Next Article in Journal
Characterizing the Countrywide Epidemic Spread of Influenza A(H1N1)pdm09 Virus in Kenya between 2009 and 2018
Next Article in Special Issue
Effects of the Newly Isolated T4-like Phage on Transmission of Plasmid-Borne Antibiotic Resistance Genes via Generalized Transduction
Previous Article in Journal
Morphological Aspects and Viremia Analysis of BALB/c Murine Model Experimentally Infected with Dengue Virus Serotype 4
Previous Article in Special Issue
Bacteriophages fEV-1 and fD1 Infect Yersinia pestis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reduced Infection Efficiency of Phage NCTC 12673 on Non-Motile Campylobacter jejuni Strains Is Related to Oxidative Stress

1
Department of Biological Sciences, University of Alberta, Edmonton, AB T6G 2E9, Canada
2
Complex Carbohydrate Research Center, Department of Microbiology, University of Georgia, Athens, GA 30602, USA
3
Ottawa Institute of Systems Biology, University of Ottawa, Ottawa, ON K1H 8M5, Canada
4
Department of Biochemistry, Microbiology and Immunology, Faculty of Medicine, University of Ottawa, Ottawa, ON K1H 8M5, Canada
*
Author to whom correspondence should be addressed.
Viruses 2021, 13(10), 1955; https://doi.org/10.3390/v13101955
Submission received: 6 July 2021 / Revised: 8 September 2021 / Accepted: 20 September 2021 / Published: 29 September 2021
(This article belongs to the Special Issue Phage-Host Interactions 2021)

Abstract

:
Campylobacter jejuni is a Gram-negative foodborne pathogen that causes diarrheal disease and is associated with severe post-infectious sequelae. Bacteriophages (phages) are a possible means of reducing Campylobacter colonization in poultry to prevent downstream human infections. However, the factors influencing phage-host interactions must be better understood before this strategy can be predictably employed. Most studies have focused on Campylobacter phage binding to the host surface, with all phages classified as either capsule- or flagella-specific. Here we describe the characterization of a C. jejuni phage that requires functional flagellar glycosylation and motor genes for infection, without needing the flagella for adsorption to the cell surface. Through phage infectivity studies of targeted C. jejuni mutants, transcriptomic analysis of phage-resistant mutants, and genotypic and phenotypic analysis of a spontaneous phage variant capable of simultaneously overcoming flagellar gene dependence and sensitivity to oxidative stress, we have uncovered a link between oxidative stress, flagellar motility, and phage infectivity. Taken together, our results underscore the importance of understanding phage-host interactions beyond the cell surface and point to host oxidative stress state as an important and underappreciated consideration for future phage-host interaction studies.

Graphical Abstract

1. Introduction

Campylobacter jejuni is an important Gram-negative foodborne pathogen worldwide that causes diarrheal disease and is associated with severe post-infectious sequelae, such as Guillain-Barré syndrome, irritable bowel syndrome, growth stunting and arthritis [1,2,3]. Furthermore, reports of antibiotic resistance continue to increase, and thus alternative solutions are required to combat this pathogen [4,5]. Since the most common Campylobacter species are commensals in chickens, many Campylobacter infections result from mishandling or consumption of contaminated poultry products [1]. As such, strategies aimed at preventing and/or at reducing its colonization of poultry flocks are predicted to have a large impact on reducing the number of human infections [5,6].
Bacteriophages (phages), the viruses that infect bacteria, are highly abundant in all environments harbouring bacteria. Because of their natural ability to specifically recognize and kill their hosts, along with their ability to evolve and adapt alongside their hosts, the exploitation of phages has been considered as a possible antimicrobial strategy [7]. As well, phage-host interaction studies have provided valuable insights into bacterial evolution [8]. Phages have been examined as a means to reduce Campylobacter colonization on farms, particularly since phages can be added directly to chicken feed and water sources [5,9]. However, several challenges associated with the broadness of strain coverage and efficiency of phage adsorption and infection must be surmounted before this strategy can be considered feasible in large-scale agricultural settings [9,10,11]. Improved understanding of Campylobacter phage-host interactions can directly contribute to the success of this strategy.
C. jejuni coats its surface with many glycoconjugates, both protein- and lipid-linked, and these structures largely contribute to strain-strain variability in this species [12,13]. Phages in turn have been shown to encode glycan binding proteins that permit them to interact with a wide variety of bacterial glycans (reviewed in [14]). However, the study of how phages depend on bacterial glycosylation pathways, and the resultant glycans, is still in its infancy. We therefore sought to better understand how Campylobacter phages interact with C. jejuni glycans. Previous work has shown that Campylobacter phages tend to require either capsular polysaccharides (CPS) or flagellar motility for adsorption to host cells [15,16]. Taxonomic classification of Campylobacter phages has clustered most CPS-dependent phages into the Fletchervirus genus, and most flagellar-dependent phages into the Firehammervirus genus [17].
Phage NCTC 12673 is a Campylobacter jejuni-specific lytic phage of the Myoviridae family of tailed phages [18]. This phage was originally used as a Campylobacter typing phage and has been shown to require CPS for C. jejuni infection [16]. It is a member of the Fletchervirus genus of the Eucampyvirinae subfamily, of which most members characterized to date have been shown to be CPS-dependent and flagella-independent [16]. However, the NCTC 12673 phage was also shown not to plaque on the paralyzed flagella cj0390 (pflB) mutant, which expresses immobilized flagella due to a lack of a functional flagellar motor complex [16]. We thus sought to better understand the involvement of flagellar motility in C. jejuni infection by NCTC 12673.
Here we describe the observation that phage NCTC 12673 cannot plaque on isogenic mutants of at least four genes (pseCFGH) in the C. jejuni NCTC 11168 pseudaminic acid biosynthesis pathway for flagellar O-glycosylation, which are also required for filament assembly and thus motility [19,20]. We show that these gene products are not required for phage adsorption, and show using whole transcriptome sequencing that their mutation does not induce expression of anti-phage or stress response pathways. Interestingly, we identified a spontaneous NCTC 12673 phage mutant, MutC, that was able to plaque efficiently on these pse mutants and that could also plaque efficiently on motA and motB mutant strains, which express paralyzed flagella [21]. Additionally, we found that MutC also plaques more efficiently on mutants in the oxidative stress defense response (catalase, ΔkatA and alkyl-hydroxyperoxidase, ΔahpC) compared to the parent phage, and is less sensitive to oxidative stress caused by the bile salt deoxycholate. This suggests that the reduced plaquing efficiency of phage NCTC 12673 on non-motile strains may be related to the fact that non-motile mutants in C. jejuni have increased sensitivity to oxidative stress [22]. These new data support and help further explain our previously reported observations that NCTC 12673 phage infection induces upregulation of oxidative stress defense genes in C. jejuni, and that mutation of these genes leads to reduced phage infection efficiency [23]. Lastly, we found that relative to the parent phage, MutC appears to more frequently encode a truncated version of a phase-variable flagellar glycan binding protein, Gp047 (recently renamed FlaGrab, [24]). These results support a model whereby phage expression of FlaGrab during infection may exacerbate oxidative stress in the cell through interaction with flagellar glycan intermediates, suggesting that the phage may have evolved a mechanism of toggling an oxidatively-stressed state in its host, and as a result, its own infection efficiency. Taken together, our results point to a complex interplay between phage infectivity, oxidative stress and the flagellar motility pathway in C. jejuni, underscoring the importance of understanding phage-host interactions beyond the cell surface.

2. Materials and Methods

2.1. Bacterial Growth Conditions

C. jejuni NCTC 11168 (MP21) [25] and its isogenic mutants and complement strains were grown on 1.5% NZCYM (Research Products International, Mt Prospect, IL, USA) agar plates, supplemented with 50 µg/mL kanamycin or 12–25 µg/mL chloramphenicol where needed, at 37 °C under microaerobic conditions (85% N2, 10% CO2, 5% O2). Escherichia coli strains were grown on LB agar, supplemented with 100 µg/mL ampicillin, 50 µg/mL kanamycin or 25 µg/mL chloramphenicol where needed. The list of bacterial strains and phages used in this study is given in Table 1.

2.2. Mutagenesis and Complementation

The mutation of pseF in C. jejuni 11168 was generated by transferring the mutation from C.jejuni 81-176ΔpseF (kindly provided by Susan Logan) according to the method described in [29] with some modifications. The pseF:cat DNA fragment was PCR-amplified from C.jejuni 81-176ΔpseF using Phusion® high-fidelity DNA polymerase. C. jejuni 11168 was transformed with the purified PCR product, transformed colonies were selected on MH agar supplemented with chloramphenicol and successful ΔpseF mutants in 11168 were confirmed by PCR.
To generate the pseH complementation construct, pseH was PCR-amplified using the primer set CS-1158/CS-1034 (5′-GGTAGATCTTTGATAAAACTTAAAAATTTCGCAGA-3′ and 5′-ATACTCGAGTTAGCTAGGCAAGGCTTTGC-3′) and cloned into BglII/XhoI-digested pAJ23 [29]. Plasmids containing the pseH construct were PCR-verified and then used to transform C. jejuni 11168 ΔpseH cells using natural transformation [29]. The complete protocol for natural transformation is available on Protocols.io at https://doi.org/10.17504/protocols.io.magc2bw (accessed on 22 September 2021).

2.3. Phage Propagation and Titration

Phage NCTC 12673 and its propagating strain C. jejuni 12661 [34] were obtained from the National Collection of Type Cultures (NCTC; Salisbury, UK). Phage MutC was also propagated on C. jejuni 12661 following initial isolation on pseC cells. Phage propagation and titration were performed following the methods described in [23].

2.4. Efficiency of Plating (EOP) Assays

EOP assays were done by spotting serial dilutions of NCTC 12673 phage onto different strains and determining the proportion of plaque forming units formed on mutant strains compared to the corresponding wild type strain, as described in [23]. Briefly, overnight bacterial cultures were harvested in NZCYM broth and set to an OD600 of 0.35. A 5-mL aliquot of this suspension was transferred to a standard sized empty Petri dish and incubated at 37 °C without shaking for 4 h under microaerobic conditions. The suspension was then set to an OD600 of 0.5, and 200 µL of this was mixed with 5 mL sterile 0.6% molten NZCYM agar at 45 °C. This suspension was poured onto the surface of a pre-warmed NZCYM plate containing 1.5% agar. Plates were allowed to solidify for 15 min and then 10 µL of serial dilutions of a phage suspension (starting at 107 PFU/mL) was spotted onto the agar surface and allowed to completely soak into the agar (15 min) before inverting the plate and incubating at 37 °C under microaerobic conditions. Plaques were counted after 18–24 h and converted to PFU/mL by multiplying countable numbers by the total dilution factor. Percent infectivity was determined by dividing PFU/mL for each test strain by the same value for the wild type reference strain. This method is available on Protocols.io (https://doi.org/10.17504/protocols.io.mahc2b6) (accessed on 22 September 2021).

2.5. Deoxycholate Assay

The effect of deoxycholate on the ability of phage to infect C. jejuni was tested under conditions previously shown to induce oxidative stress [35] with some modifications. Overnight cultures were harvested and set to OD600 of 0.5 in 10 mL of BHI broth in a polystyrene petri dish. Deoxycholate was added to a final concentration of 0.05% and the cultures were incubated for 18 h in microaerobic conditions. After growth in deoxycholate, cultures were collected, washed with fresh broth and the EOP assay was performed as described above.

2.6. Adsorption Assays

Adsorption assays were done according to [36] with some modifications. Briefly, cells were harvested as detailed above for growth curves after washing pellets three times in 1 mL NZCYM broth and then set to an OD600 of 0.4 in 2 mL NZCYM broth in a small empty Petri plate. Then, 20 μL of 0.22-μm-filtered phage lysate at 5.6 × 105 PFU/mL was added and mixed by pipetting. The number of unadsorbed phages at t = 0 was determined by removing 500 μL immediately after mixing phage and host and centrifuging at 15,000 rpm for 2 min at 4 °C. Then, 100 μL from the supernatant was removed and kept on ice until all samples were ready for titration. After the t = 0 sample was taken, phage-cell suspensions were incubated aerobically with 100 rpm shaking at 37 °C for 90 min, at which point the sampling process was repeated. Samples from each time point were serially diluted in SM buffer and titered for PFU/mL to determine the number of unadsorbed phages over time. Percent adsorption was determined by dividing PFU/mL of unadsorbed phage particles in each test strain by the same value for the wild type reference strain.

2.7. Total RNA Extraction

C. jejuni 11168 wild type, ΔpseC and ΔpseF mutant cells were harvested from overnight NZCYM plate cultures, pelleted and washed once in NZCYM broth and set to an OD600 of 0.05 (2 × 108 colony forming units per mL (CFU/mL)) in 20 mL NZCYM broth in 125-mL Erlenmeyer flasks, each containing a 1-inch sterile magnetic stir-bar. Cells were grown under microaerobic conditions and magnetically stirred at 200 rpm. After 4.5 h incubation (mid-log phase, cell counts were approximately 5 × 108 CFU/mL), the entire contents of each flask was transferred to a pre-prepared tube containing 2.6 mL (0.1 volume) ice cold 10% buffered phenol in 100% ethanol to stabilize RNA followed by immediate mixing and storage on ice until all samples were collected [37]. RNA was extracted from each sample using a hot phenol method [37]. RNA samples were sequentially DNAse-treated (37 °C for 30 min) using RNAse-free DNAse I (Epicentre, Madison, WI, USA) and cleaned using the Zymo RNA Clean & Concentrator (Zymo Research, Irvine, CA, USA). PCR was used to confirm the absence of residual DNA. Total RNA quality was assessed using an Agilent Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) and RNA was stored at −80 °C until further use. Samples were extracted in biological triplicate.

2.8. RNA-Sequencing

Total mRNA libraries from all replicates were generated. Samples were depleted of rRNA using the RiboZero bacterial kit (Illumina, San Diego, CA, USA) according to the manufacturer’s instructions. Successful rRNA depletion was confirmed using the Agilent Bioanalyzer RNA 6000 Pico Kit (Agilent Technologies, Santa Clara, CA, USA). Strand-specific barcoded sequencing libraries were constructed using the Ion Total RNA-seq kit (Thermo Fisher Scientific, Waltham, MA, USA). Libraries were quality-checked and quantified using the Bioanalyzer High Sensitivity DNA kit (Agilent Technologies, Santa Clara, CA, USA) and pooled together in equimolar amounts. The pooled libraries were templated using the Ion PI Hi-Q kit (Thermo Fisher Scientific, Waltham, MA, USA) and sequenced on an Ion Torrent Proton using the Ion PI Hi-Q sequencing 200 kit (Thermo Fisher Scientific, Waltham, MA, USA) on a single Proton V2 chip.
The raw sequencing reads were demultiplexed by the Ion Torrent suite software (version 5.2.2., Thermo Fisher Scientific, Waltham, MA, USA) and sequentially mapped to the host (NCTC 11168) and phage (NCTC 12673) genomes using Spliced Transcripts Alignment to a Reference (STAR) [38] (Table S1). Reads aligning to coding regions were counted using HT-seq using the default settings. The raw demultiplexed sequencing reads have been deposited at the NCBI Sequence Read Archive (SRA) archive under accession number PRJNA743365. DESeq2 was used to identify differentially expressed transcripts between each strain [39]. Genes with a fold change +/− 1.5 and false discovery rate (FDR)-corrected p-value < 0.05 were considered differentially expressed.

2.9. Isolation of MutC

MutC was isolated as a clear plaque following an EOP assay (as described above) whereby serial dilutions of phage NCTC 12673 were spotted onto a plate of C. jejuni NCTC 11168 pseC mutant cells. The plaque was isolated and purified according to the methods of [40] and propagated on C. jejuni 12661, the propagating strain for the parent phage, NCTC 12673.

2.10. Whole Genome Sequencing and Analysis

Phage DNA (MutC and NCTC 12673) was prepared using hot phenol-chloroform extraction following pre-amplification with phi29 polymerase, as described in [18]. Genomic DNA was submitted for Illumina library preparation and whole-genome sequencing to The Applied Genomics Core (Edmonton, AB, Canada). Genomes of NCTC 12673 (resequenced during this work) and MutC were aligned against the NCTC 12673 reference strain in GenBank (accession: GU296433) [18].
Genetic differences between NCTC 12673 and MutC phages, as compared to the reference strain, were identified by variant calling using the Illumina sequencing reads using Geneious (version 8).

3. Results

3.1. Phage NCTC 12673 Requires a Functional Pseudaminic Acid Biosynthetic Pathway for Infection

To determine whether flagellar glycosylation is important for NCTC 12673 infection, we examined phage infection of ΔpseC, ΔpseF, ΔpseG, ΔpseH mutants and found that compared to wild type C. jejuni NCTC 11168 infection (EOP set to 100%), infectivity was decreased on all mutants (Figure 1). No plaquing was observed on ΔpseC or ΔpseH, and infectivity relative to wild type was drastically reduced upon infection of pseF (EOP = 1.7%; p = 0.0001) and pseG (EOP = 3.3%; p = 0.0001). We complemented ΔpseH with a functional copy of the gene and found that infectivity was significantly improved in the ΔpseH/+pseH complemented strain (EOP = 46.6%; p = 0.0044) compared to the pseH mutant. As a negative control, we tested infectivity on kpsM, which lacks the capsule required for NCTC 12673 binding and is known not to be infected by the NCTC 12673 phage [16].

3.2. Phage NCTC 12673 Adsorbs to Pseudaminic Acid Pathway Mutants at Wild Type Levels

It has been well established that inactivation of pseB, pseC, pseE, pseF, pseG, pseH or pseI, which all encode enzymes required for biosynthesis and transfer of pseudaminic acid to flagella in C. jejuni, prevents flagellar biogenesis and motility [20,42] (Figure S1). As other Campylobacter phages have been shown to depend on flagellar motility for infection, and this has been attributed to the use of flagella as a surface receptor for the phage [36], we first sought to determine whether the reason NCTC 12673 infectivity was impaired on pse mutants was due to a lack of adsorption to these aflagellate mutants. We titered the number of free phages remaining in the supernatant following incubation with C. jejuni NCTC 11168 wild type, ΔpseC, ΔpseF, ΔpseG and ΔpseH mutants, a ΔpseH/ + pseH complement strain, and a ΔkpsM mutant as a negative control. With the exception of ΔpseH, which showed a slight, yet statistically significant decrease in adsorption compared to wild type cells (p = 0.002), we observed no differences in adsorption for most mutant strains relative to wild type cells (Figure 2). In line with previous results, NCTC 12673 phage did not adsorb to ΔkpsM mutant cells (p = 0.0002).

3.3. ΔpseC and ΔpseF Mutants Display No Evidence of Stress Response or Phage Defense, but Downregulate Many Flagellar Genes

We next sought to determine whether downstream effects of pse mutation led to stress responses or expression of any anti-phage defense genes, which might help explain why the NCTC 12673 phage is unable to efficiently infect these strains. To test this, we used RNA-seq to compare gene expression between C. jejuni NCTC 11168 wild type, ΔpseC and ΔpseF mutant cells during the mid-log phase of growth (Figure S2). We identified widespread downregulation of many genes in both ΔpseC compared to wild type cells (Figure 3A, Table S1) and ΔpseF compared to wild type cells (Figure 3B, Table S1). For both mutant strains relative to wild type cells, most of the significantly downregulated genes were related to flagellar biogenesis (e.g., flgE2, flaB, flgE, flgD, flgI, flgK), with the next most altered category being amino acid biosynthesis (e.g., gltB, glnA). cj0501, a pseudogene [43], was downregulated in ΔpseC and ΔpseF. Interestingly, relative to wild type and ΔpseF, the ΔpseC mutant showed substantial upregulation of metA and metB genes (encoding enzymes for methionine biosynthesis). We also observed the slight upregulation of the CPS biosynthesis genes (cj1422c, hddA) in ΔpseC mutant cells compared to wild type cells. Direct comparison of ΔpseF to ΔpseC showed limited differences but did highlight that metB and metA were downregulated in ΔpseF cells compared to ΔpseC cells, as were pseG, pseH, and cj1422c, while asd, cj1295, and cj1022c were expressed more highly in ΔpseF than ΔpseC (Figure 3C, Table S1).

3.4. Phage NCTC 12673 Is Unable to Infect Cells in the Absence of the Flagellar Motor Proteins MotA or MotB

To better understand the possible NCTC 12673 dependence on flagellar motility, and to verify previous results that showed that this phage did not plaque on a “paralyzed” (full-length but non-motile) flagella ΔpflB mutant, we tested NCTC 12673 phage plaquing efficiency on flagellar motor mutants ΔmotA and ΔmotB [21,24]. Interestingly, we found that NCTC 12673 phage infection was essentially abolished on both mutants (Figure 1A).

3.5. Oxidative Stress Sensitivity of Non-Motile Mutant Strains May Explain Reduced NCTC 12673 Plaquing Efficiency

We recently showed that the NCTC 12673 phage displays reduced plaquing on mutants in the oxidative stress defense genes ΔkatA, ΔahpC and ΔsodB [23]. As well, Flint et al. previously showed that many non-motile C. jejuni NCTC 11168 mutants, including ΔflgK, ΔflgH, ΔflgD, ΔflgI, ΔflgR, ΔpseB, and ΔmotAB, are hypersensitive to oxidative stress, presumably through disruption of the proton gradient/electron leakage through the electron transport chain increasing reactive oxygen species [22]. With the exception of ΔmotA and ΔmotB, we observed that all of these genes highlighted to lead to oxidative stress sensitivity by Flint et al. were among the most significantly downregulated genes across the transcriptome in both ΔpseC and ΔpseF mutant strains (Figure 3A,B). We therefore hypothesized that the observed dependence of NCTC 12673 on flagellar motility might be linked to increased oxidative stress sensitivity in these strains.

3.6. A spontaneous NCTC 12673 Mutant Phage, “MutC”, Efficiently Plaques on Both Non-Motile and Oxidative Stress Defense Mutants

To better understand the mechanism for motility dependence by NCTC 12673, we sought to characterize an NCTC 12673 phage escape mutant, MutC, which we isolated as a clear plaque on ΔpseC mutant cells. Interestingly, we found that MutC displayed a significantly higher EOP on all of the Δpse and Δmot mutants tested compared to the parent phage, NCTC 12673 (Figure 4). While the parent phage displayed undetectable plaquing on ΔpseC, ΔpseH, ΔmotA, and ΔmotB, and low efficiency on ΔpseF and ΔpseG (Figure 4), phage MutC plaqued at a statistically significantly higher efficiency than the parent phage on ΔpseC (15.7%; p = 0.001), ΔpseF (26.4%; p = 0.002), ΔpseG (48.4%; p = 0.003), ΔpseH (9.8%; p = 0.01), ΔmotA (75.6%; p = 0.005), and ΔmotB (57.3%; p = 0.02), as determined by t-test. Similar to the parent phage, MutC did not infect ΔkpsM cells, suggesting that it retained dependence on C. jejuni CPS.
To determine whether the ability to overcome dependence on motility by MutC also allowed this phage to overcome dependence on oxidative stress defense genes, we next tested MutC infection of mutants in the oxidative stress defense genes ΔkatA (catalase), ΔahpC (alkyl-hydroxyperoxide reductase), and ΔsodB (superoxide dismutase) [31]. We found that MutC more efficiently infected ΔkatA (83.0% vs. 28.6%; p = 0.01) and ΔahpC (79.2% vs. 3.7%; p = 0.00005) mutant cells compared to phage NCTC 12673 (Figure 4). Conversely, we did not observe a difference between MutC and NCTC 12673 infection of the ΔsodB mutant (75.0% vs. 61.2%; p = 0.62).

3.7. Phage MutC Is Less Impacted by Exposure to the Oxidative Stress-Inducing Agent Deoxycholate Than Phage NCTC 12673

To further demonstrate that phage MutC had gained an improved ability to infect oxidatively-stressed cells, we compared MutC and the parent phage NCTC 12673 infection efficiency on C. jejuni 11168 following exposure to 0.05% deoxycholate, a bile salt previously shown to induce the production of intracellular reactive oxygen species in C. jejuni when supplemented in the growth medium [35]. We found that MutC indeed had a significantly improved (5.8-fold) infection efficiency following growth in the presence of deoxycholate compared to NCTC 12673 (73.5% vs. 12.6%, respectively; p = 0.001) (Figure 4).

3.8. Genomic Comparison between Phages NCTC 12673 and MutC Predicts Differences in Several Proteins, including FlaGrab, a Flagellar Glycan-Binding Protein

In order to understand the mechanism for the observed differences in NCTC 12673 and MutC plaquing on non-motile and ΔkatAahpC mutants, we performed whole-genome sequencing of both phages and assessed differences in their genomes. Although the two phages displayed nucleotide changes at approximately a dozen locations across their 135-kbp genomes (Table S2), amino acid-level consequences were predicted in only six genes: gp047/flagrab, which encodes the flagellar glycan-binding protein FlaGrab [24,29,30]; gp041, annotated as a baseplate wedge protein homologue; gp058 and gp167, both annotated as homing endonucleases; and gp114 and gp116, both hypothetical protein-encoding genes (Table 2). We chose to focus on the flagrab variation for follow-up analysis since our previous studies documented associations between FlaGrab and the C. jejuni flagellar motility pathway, leading us to surmise that it might be implicated in the NCTC 12673/MutC flagellar mutant infectivity differences observed.

3.9. Differences in Variant Frequency within Two Poly-Adenosine Nucleotide Tracts Suggests MutC Phage Expresses a Truncated FlaGrab Protein More Frequently Than NCTC 12673 Phage

FlaGrab (encoded by gp047/flagrab) is a flagellar glycan binding protein encoded by all Campylobacter phages characterized to date [17], which we have previously shown binds to acetamidino-modified pseudaminic acid glycans on C. jejuni flagella [29], inhibits C. jejuni growth [24,30] and reduces expression of energy metabolism genes following binding to motile flagella [24]. Within the NCTC 12673 phage genome, we identified the presence of two poly-A tracts within flagrab. Variation in the length of either or both poly-A tracts is predicted to lead to frameshift mutations that would severely truncate the protein; when both tracts contain 7 As, the predicted FlaGrab protein is 1365 amino acids in length, whereas encoding 6 As at one or both tracts results in a protein of less than or equal to 83 amino acids. By aligning the reads generated by short-read whole genome sequencing (Illumina) data for MutC and NCTC 12673 phages to the NCTC 12673 reference genome in GenBank (accession: GU296433) [18], we analyzed the length of each of the two poly-A tracts within the flagrab gene for both phages and compared the frequency of 6-A vs. 7-A tracts to determine whether we would expect differences in FlaGrab length expressed by the two phages (Table 2). At the first poly-A tract (chromosomal position: 40,821), we observed that MutC encoded the 6-A allele at a frequency of 68.80%, whereas NCTC 12673 did not show variability, encoding only the 7-A allele. At the second tract (chromosomal position: 40,935), we observed the 6-A allele in similar proportions for both phages (NCTC 12673: 33.20%, MutC: 25.30%). These data suggest that while both phages would be predicted to express a truncated FlaGrab as a result of variation at the second tract approximately one third of the time, phage MutC is predicted to express a truncated FlaGrab most of the time, and thus far more frequently than NCTC 12673, due to variation at both tracts.
To verify the flagrab variability observed upon whole-genome sequencing, we designed primers to PCR-amplify the poly-A tract regions of flagrab. We then performed Sanger sequencing on the resultant amplicons using NCTC 12673 and MutC phage lysates of C. jejuni NCTC 12661 cells (the strain routinely used to propagate these phages) as templates. Indeed, we found variability in sequence data (represented by multiple chromatogram peaks) starting directly downstream of the flagrab poly-A tracts for both phages, supporting the variability predicted from the variant frequency analysis described above.

4. Discussion

Phage NCTC 12673 is classified as a member of the Fletchervirus genus of Campylobacter phages, members of which tend to target CPS and not flagella [15,17,18]. Coward et al. (2006) screened this phage against a C. jejuni mutant library and showed that acapsular mutants were consistently resistant to the phage [16], and we have confirmed that acapsular mutants are not infected or bound by the phage [29] (this work). However, Coward et al. also found that there was one non-motile mutant that NCTC 12673 could not infect: the paralyzed flagellar mutant ΔpflB [16]. This observation suggested that in addition to dependence on CPS for infection, motile flagella may also be important for successful infection of C. jejuni cells by phage NCTC 12673.
We sought to elucidate whether host cell motility was indeed important for NCTC 12673 infection, and if so, why. We analyzed whether phage NCTC 12673 could plaque on several mutants of the pseudaminic acid (pse) flagellar glycosylation pathway, which are known to be aflagellate and non-motile [44] (Figure S1). We found that mutagenesis of genes in the pse pathway (pseC, pseF, pseH, or pseG) drastically reduced NCTC 12673 phage ability to infect C. jejuni NCTC 11168 cells, and that complementation of pseH in trans restored flagellar motility and improved phage infection efficiency. These results supported our hypothesis that NCTC 12673 phage may indeed have a dependency on C. jejuni flagellar motility.
One of the simplest explanations for why a phage might require flagellar motility is that it uses motile flagella as a host cell surface receptor and binds to the flagella as part of its initial interaction with the host; indeed, others have shown that this is the case for the flagella-requiring Campylobacter phage F341 [36]. To determine if this was the case for NCTC 12673, we assessed adsorption of the NCTC 12673 phage to our set of pse mutants. Surprisingly, we found that the phage could still adsorb to the pse mutants at wild type levels, ruling this out as an explanation. We thus hypothesized that pse mutation may impact phage infection at a later point in the phage lifecycle such as during the intracellular phase of phage infection. We thus sought to determine whether we could detect any other impacts of pse mutations on C. jejuni cells that might help explain why the phage could not infect these mutants, such as host stress responses or anti-phage defense that might be detrimental to phage replication, or other unexpected downstream effects of pse mutagenesis.
Since downstream effects on gene expression in response to gene mutagenesis are common, we were interested in understanding the total effects of pse gene inactivation to identify possible explanations for the lack of phage infectivity of pse mutants. To test this, we performed RNA-seq on ΔpseC and ΔpseF mutant cells. We did not observe signs of general stress response activation in the ΔpseC and ΔpseF mutants relative to wild type cells, or increased expression of any known anti-phage defense genes [23,45], but we did observe widespread downregulation of many flagellar genes. This downregulation was not surprising, since C. jejuni is known to tightly regulate flagellar expression at multiple levels including hierarchical expression of genes through distinct sigma factors (σ70, σ54, σ28), feedback inhibition by σ28 expression, global transcriptional regulation by CsrA-FliW [46,47], PseB inhibition by Pse glycans and pse gene requirement for filament biogenesis [41,48,49,50,51]. It is therefore expected that the cell would downregulate σ54- and σ28-regulated flagellar biogenesis genes in response to the lack of σ28-regulated pse gene expression, which is what we observed.
Interestingly, Flint et al. (2014) showed that several non-motile C. jejuni mutants in genes downregulated in our ΔpseC and ΔpseF mutants are hypersensitive to oxidative stress [22]. Additionally, our recent transcriptomic analysis of NCTC 12673 phage during lytic infection of C. jejuni 11168 showed that upon infection, C. jejuni upregulates the oxidative stress defense gene katA as well as the heme acquisition genes chuABCDZ, which import the heme co-factor required for KatA activity [23]. This led to the finding that ΔkatA, as well as ΔahpC and ΔsodB, which constitute other oxidative stress response genes, are important for efficient phage infection by NCTC 12673, suggesting that this phage cannot efficiently infect oxidatively-stressed cells [23]. In support of this, Chatterjee et al. (2020) have recently shown that oxidative stress response genes are also important for Enterococcus faecalis phage infection, suggesting that oxidative stress may play an important role in dictating phage-host dynamics more broadly [52]. We therefore hypothesized that increased sensitivity to oxidative stress of pse mutants, brought on by their lack of motility, might explain the inability of NCTC 12673 phage to plaque on the pse mutants.
During our experiments plaquing phage NCTC 12673 on the pseC mutant, we isolated a spontaneous phage variant that could form clear plaques on ΔpseC mutant cells. We decided to characterize this variant, which we named MutC, to see if it might provide clues into why pse mutation impacted NCTC 12673 infection efficiency. We tested MutC’s ability to plaque on the other pse mutants and found that in addition to plaquing on ΔpseC cells, MutC could efficiently infect mutants in ΔpseC, ΔpseF, ΔpseG and ΔpseH, as well as the flagellar filament mutant ΔflaA and the paralyzed flagellar mutants ΔmotA and ΔmotB.
To validate the hypothesis that increased oxidative stress was the reason the pse and mot mutants were not efficiently infected by phage NCTC 12673, we next sought to determine whether MutC could plaque on C. jejuni mutants lacking oxidative stress defense genes. Interestingly, we found that MutC was able to plaque more efficiently than the parent phage on the oxidative stress defense mutants ΔkatA and ΔahpC. In other words, the mutation(s) in MutC that allowed it to overcome dependence on pse genes also allowed it to overcome dependence on oxidative stress defense genes, providing a clear link between pse mutation and oxidative stress defense. This result supported our hypothesis that NCTC 12673′s inability to infect pse and mot mutants may be caused by oxidative stress sensitivity of these mutants. Of note, NCTC 12673 and MutC similarly infected a ΔsodB mutant, suggesting that the phages may be sensitive to specific oxidative stresses, as hydrogen peroxides/alkyl peroxides are detoxified by KatA/AhpC while superoxides are detoxified by SodB [31]. Alternatively, this could indicate the involvement of other regulators known to impact oxidative stress regulation, such as the Campylobacter oxidative stress regulator, CosR [53].
To further demonstrate that phage MutC had gained an improved ability to infect oxidatively-stressed cells, we compared phage infection efficiency on C. jejuni 11168 of MutC and the parent phage NCTC 12673 following exposure to 0.05% deoxycholate, a bile salt known to induce the production of reactive oxygen species in C. jejuni [35]. We found that MutC indeed had 5.8-fold higher infection efficiency following deoxycholate exposure compared to NCTC 12673. This provided further support to the hypothesis that the sensitivity of phage NCTC 12673 to oxidative stress explains its apparent dependence on a functional flagellar pathway.
Reanalyzing our transcriptomic data given these new findings about oxidative stress sensitivity of the pse mutants, we found it noteworthy that the ΔpseC mutant displayed higher expression of the methionine biosynthesis genes metA and metB compared with wild type or the ΔpseF mutant. Others have proposed that increased methionine biosynthesis may be necessary to regenerate oxidized methionines that exceed the capacity of the methionine sulfoxide reductase repair system, and/or for the increased production of oxidative repair enzymes such as KatA and AhpC, which are known for their high methionine content/quenching abilities in the related pathogen Helicobacter pylori [54,55]. It should also be noted that C. jejuni does not encode genes for the production of glutathione [56], a key redox buffer involved in the repair of oxidized methionine/cysteine [57]. Of note, our results suggest that mutation of early pse pathway genes (pseC, pseH) is more detrimental to phage infection compared to mutation of late pse pathway genes (pseG, pseF), but the reason for this remains unknown. It is possible that differences in metA, metB, CPS biosynthesis and/or amino acid biosynthesis may explain some of the difference between NCTC 12673 infectivity of these strains.
To understand the molecular mechanism driving the observed differences in host requirements between the two phages, we sought to identify the genetic differences between NCTC 12673 and MutC. While we identified several genes with variations between the two phages, only five showed amino acid-level changes. Interestingly, one of the amino acid-level differences between the two phages was in FlaGrab, a conserved Campylobacter phage protein that we previously showed is a flagellar glycan-binding protein specific for acetamidino-modified pseudaminic acid [29] that causes growth inhibition in C. jejuni, presumably also through disruption of the PMF leading to oxidative stress [24]. We chose to focus on the FlaGrab variation for follow-up analysis, as its documented association with the C. jejuni flagellar motility pathway led us to surmise that it might be implicated in the NCTC 12673/MutC flagellar mutant infectivity differences observed.
Interestingly, the difference we observed in the flagrab gene between the two phages was in the frequency of variation at a 6–7-nucleotide poly-A tract, predicted to lead to a frameshift mutation in flagrab more frequently in MutC than in NCTC 12673. In fact, phase variability through hypermutable poly-G tracts (a rapid adaptation mechanism that provides a stochastic on/off switch for gene expression via homopolymeric nucleotide stretches encoded within open reading frames [58]) has recently been reported in Campylobacter phages [59]. As C. jejuni expresses at least 30 phase-variable genes, many of them governing surface structure expression [26,58,60], it is not surprising to find that Campylobacter phages have evolved similar mechanisms to ensure evolutionary dominance over such a variable host. Tracts of only 6–7 nucleotides in length would not be predicted to switch on and off as frequently as longer tracts, given the slipped-strand mispairing mechanism that leads to variation at these loci, but would still be reasonably expected to represent a mechanism of variation that could provide an adaptive advantage for the encoding biological entity. The differences in variant frequency at the poly-A tracts in flagrab that we observed suggest that a much larger proportion of the MutC population expresses a severely truncated (or “off-switched”) version of FlaGrab compared to the parent phage NCTC 12673. Although we have extensively described the effects of adding recombinantly purified FlaGrab to C. jejuni cells, namely that binding of this protein to C. jejuni flagella leads to growth inhibition [24,30], the function of FlaGrab in the phage life cycle has remained elusive. And yet, every Campylobacter phage characterized to date has been shown to encode a FlaGrab homologue, suggesting it is of high importance to any phage infecting this species [17]. We are still exploring how FlaGrab may or may not be involved in the phenotypes observed, but based on our results to date, we hypothesize that full-length FlaGrab may prevent NCTC 12673 phage from lysing its host under conditions of oxidative stress, while MutC, encoding a truncated FlaGrab, is able to constitutively lyse its host. How this occurs is currently purely speculative, but our hypothesis is that FlaGrab may have another function intracellularly, as many other phage proteins have been shown to have [61]. Our future directions will seek to establish whether a definitive link can be made between the length of FlaGrab and phage infection efficiency.
Although most Campylobacter phages studied to date are lytic and not lysogenic, a ‘carrier state’ association between two NCTC 12673-related phages and C. jejuni has been described [62]. This is defined as a longer-term association between the phage and host, where the phage associates with the host both intracellularly (episomally) and extracellularly (by associating with the surface) and waits to lyse the cell for a period of time, without integrating into the genome. The carrier state in C. jejuni has been shown to lead to downregulation of flagellar motility, suggesting that motility and the phage life cycle may be linked for phages targeting this organism [62]. However, the mechanism for how C. jejuni phages might accomplish this has not been elucidated. Given our results here, which point toward FlaGrab having a role in blocking the lytic phage life cycle under non-motile conditions, it is tempting to speculate that the lack of phage infection efficiency we have described here is in fact representative of the carrier state, with FlaGrab representing the switch that dictates the lysis/no-lysis decision.
Together our data suggest a possible role for flagrab in dictating NCTC 12673 phage interactions with oxidatively stressed C. jejuni cells. However, the data do not rule out possible contributions of the other observed genetic changes between the two phage variants, and the effects of these changes should also be examined. Our preliminary work suggests the L to S amino acid change in MutC gp041, predicted to encode a baseplate wedge subunit protein, is reproducible (could be confirmed by Sanger sequencing of phage genomic DNA), suggesting that this change could also play a role in the phenotypic differences between the phages. In addition, analysis of gp114 and gp116 has shown that these genes exhibit phase-variability via poly-nucleotide tracts. However, we found these tracts to be highly variable following both ΔpseC and wild type infection for both NCTC 12673 and MutC phages, thus reducing the likelihood that on/off switching in these genes is responsible for the stable difference in phenotypes observed between the two phages. Finally, our efforts to confirm the observed differences in the homing endonucleases gp058 and gp167 between NCTC 12673 and MutC phages showed that the nucleotide changes observed upon whole genome sequencing were not reproducible, suggesting against a role for these changes in explaining the phenotypic differences between NCTC 12673 and MutC.

5. Conclusions

Overall, this work provides novel insights into NCTC 12673 phage interactions with C. jejuni 11168. We have shown that the apparent dependence of phage NCTC 12673 on C. jejuni motility is not related to adsorption to the flagella, but instead can be explained by the phage’s reduced infection efficiency under conditions of oxidative stress, previously shown to be the case for non-motile C. jejuni mutants. We have also shown evidence that phage-encoded FlaGrab, a flagellar glycan binding protein encoded by all Campylobacter phages with an unknown role in the phage life cycle, is phase-variable, and that variable expression of this protein may be linked to a mechanism by which the phage can tune its infectivity according to the oxidative stress state of the cell. Taken together, our results point to a complex interplay between phage infectivity, oxidative stress and flagellar motility in C. jejuni, underscoring the importance of understanding phage-host interactions beyond the cell surface.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/v13101955/s1, Figure S1: Transmission electron micrographs showing presence of flagella in wild type C. jejuni 11168 cells (A), absence of flagella in pseC-H mutants (BE), and restoration of flagella in the pseH/+pseH complement strain. Figure S2: Principal Component Analysis plot for C. jejuni NCTC 11168 wild type and mutant cells: wild type, ΔpseC, and ΔpseF. Table S1: Differentially expressed genes for C. jejuni NCTC 11168 wild type and mutant cells: wild type, ΔpseC, and ΔpseF. Table S2: Complete list of genetic polymorphisms and variant frequencies between NCTC 12673 phage (parent) and MutC, a spontaneous mutant phage isolated on C. jejuni 11168 ΔpseC mutant cells.

Author Contributions

Conceptualization, J.C.S., M.A.J. and C.M.S.; Data curation, J.B.; Formal analysis, J.C.S. and J.B.; Funding acquisition, C.M.S. and A.S.; Investigation, J.C.S., M.A.J., A.F. and C.S.C.; Methodology, J.C.S., M.A.J., J.B. and A.F.; Resources, A.S. and C.M.S.; Supervision, A.S. and C.M.S.; Writing—original draft, J.C.S.; Writing—review & editing, J.C.S., A.F., J.B., A.S. and C.M.S. All authors have read and agreed to the published version of the manuscript.

Funding

J.C.S. is a recipient of an NSERC Alexander Graham Bell Canada Graduate Student Scholarship. C.M.S. was an Alberta Innovates Strategic Chair in Bacterial Glycomics. Research in Alain Stintzi′s laboratory is supported by CIHR grant MOP#84224, the Government of Canada through Genome Canada and the Ontario Genomics Institute (OGI-149), the Canadian Institutes of Health Research (ECD-144627), and the Ontario Ministry of Economic Development and Innovation (Project 13440).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw demultiplexed RNA-sequencing reads have been deposited at the NCBI Sequence Read Archive (SRA) archive under accession number PRJNA743365.

Acknowledgments

We thank David Simpson, Stephen Pollo, and Zachary Lewis for assistance with and advice on phage genome assembly and analysis. We thank Hayley Reynolds for help carrying out phage experiments, Stephen Andersen for help with mutagenesis and complementation of C. jejuni mutants, and Susan Logan for the C. jejuni 81-176ΔpseF strain. This work was adapted from work published in J.C.S.’s doctoral thesis [63].

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kaakoush, N.O.; Castaño-Rodríguez, N.; Mitchell, H.M.; Man, S.M. Global epidemiology of Campylobacter infection. Clin. Microbiol. Rev. 2015, 28, 687–720. [Google Scholar] [CrossRef] [Green Version]
  2. Yuki, N.; Susuki, K.; Koga, M.; Nishimoto, Y.; Odaka, M.; Hirata, K.; Taguchi, K.; Miyatake, T.; Furukawa, K.; Kobata, T.; et al. Carbohydrate mimicry between human ganglioside GM1 and Campylobacter jejuni lipooligosaccharide causes Guillain-Barre syndrome. Proc. Natl. Acad. Sci. USA 2004, 101, 11404–11409. [Google Scholar] [CrossRef] [Green Version]
  3. Amour, C.; Gratz, J.; Mduma, E.; Svensen, E.; Rogawski, E.T.; McGrath, M.; Seidman, J.C.; McCormick, B.J.J.J.; Shrestha, S.; Samie, A.; et al. Etiology, Risk Factors, and Interactions of Enteric Infections and Malnutrition and the Consequences for Child Health and Development Project (MAL-ED) Network Investigators Epidemiology and Impact of Campylobacter Infection in Children in 8 Low-Resource Settings: Results from the MAL-ED Study. Clin. Infect. Dis. 2016, 63, 1171–1179. [Google Scholar] [CrossRef] [Green Version]
  4. Hampton, T. Report Reveals Scope of US Antibiotic Resistance Threat. JAMA 2013, 310, 1661. [Google Scholar] [CrossRef]
  5. Johnson, T.J.; Shank, J.M.; Johnson, J.G. Current and Potential Treatments for Reducing Campylobacter Colonization in Animal Hosts and Disease in Humans. Front. Microbiol. 2017, 8, 487. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Wassenaar, T.M. Following an imaginary Campylobacter population from farm to fork and beyond: A bacterial perspective. Lett. Appl. Microbiol. 2011, 53, 253–263. [Google Scholar] [CrossRef] [PubMed]
  7. Abedon, S.; Kuhl, S.J.; Blasdel, B.G.; Kutter, E.M. Phage treatment of human infections. Bacteriophage 2011, 1, 66–85. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. De Smet, J.; Hendrix, H.; Blasdel, B.G.; Danis-Wlodarczyk, K.; Lavigne, R. Pseudomonas predators: Understanding and exploiting phage-host interactions. Nat. Rev. Microbiol. 2017, 15, 517–530. [Google Scholar] [CrossRef]
  9. Connerton, P.L.; Timms, A.R.; Connerton, I.F. Campylobacter bacteriophages and bacteriophage therapy. J. Appl. Microbiol. 2011, 111, 255–265. [Google Scholar] [CrossRef]
  10. El-Shibiny, A.; Connerton, P.L.; Connerton, I.F. Campylobacter succession in broiler chickens. Vet. Microbiol. 2007, 125, 323–332. [Google Scholar] [CrossRef]
  11. El-Shibiny, A.; Scott, A.; Timms, A.; Metawea, Y.; Connerton, P.; Connerton, I. Application of a group II Campylobacter bacteriophage to reduce strains of Campylobacter jejuni and Campylobacter coli colonizing broiler chickens. J. Food Prot. 2009, 72, 733–740. [Google Scholar] [CrossRef]
  12. Karlyshev, A.V.; Ketley, J.M.; Wren, B.W. The Campylobacter jejuni glycome. FEMS Microbiol. Rev. 2005, 29, 377–390. [Google Scholar]
  13. Nothaft, H.; Szymanski, C.M. Protein glycosylation in bacteria: Sweeter than ever. Nat. Rev. Microbiol. 2010, 8, 765–778. [Google Scholar] [CrossRef] [PubMed]
  14. Simpson, D.J.; Sacher, J.C.; Szymanski, C.M. Exploring the interactions between bacteriophage-encoded glycan binding proteins and carbohydrates. Curr. Opin. Struct. Biol. 2015, 34, 69–77. [Google Scholar] [CrossRef] [PubMed]
  15. Sørensen, M.C.H.; Gencay, Y.E.; Birk, T.; Baldvinsson, S.B.; Jäckel, C.; Hammerl, J.A.; Vegge, C.S.; Neve, H.; Brøndsted, L. Primary isolation strain determines both phage type and receptors recognised by Campylobacter jejuni bacteriophages. PLoS ONE 2015, 10, e0116287. [Google Scholar] [CrossRef] [Green Version]
  16. Coward, C.; Grant, A.J.; Swift, C.; Philp, J.; Towler, R.; Heydarian, M.; Frost, J.A.; Maskell, D.J. Phase-variable surface structures are required for infection of Campylobacter jejuni by bacteriophages. Appl. Environ. Microbiol. 2006, 72, 4638–4647. [Google Scholar] [CrossRef] [Green Version]
  17. Javed, M.A.; Ackermann, H.W.; Azeredo, J.; Carvalho, C.M.; Connerton, I.; Evoy, S.; Hammerl, J.A.; Hertwig, S.; Lavigne, R.; Singh, A.; et al. A suggested classification for two groups of Campylobacter myoviruses. Arch. Virol. 2014, 159, 181–190. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Kropinski, A.M.; Arutyunov, D.; Foss, M.; Cunningham, A.; Ding, W.; Singh, A.; Pavlov, A.R.; Henry, M.; Evoy, S.; Kelly, J.; et al. Genome and proteome of Campylobacter jejuni bacteriophage NCTC 12673. Appl. Environ. Microbiol. 2011, 77, 8265–8271. [Google Scholar] [CrossRef] [Green Version]
  19. Goon, S.; Kelly, J.F.; Logan, S.M.; Ewing, C.P.; Guerry, P. Pseudaminic acid, the major modification on Campylobacter flagellin, is synthesized via the Cj1293 gene. Mol. Microbiol. 2003, 50, 659–671. [Google Scholar] [CrossRef]
  20. Guerry, P.; Szymanski, C.M. Campylobacter sugars sticking out. Trends Microbiol. 2008, 16, 428–435. [Google Scholar] [CrossRef]
  21. Beeby, M.; Ribardo, D.A.; Brennan, C.A.; Ruby, E.G.; Jensen, G.J.; Hendrixson, D.R. Diverse high-torque bacterial flagellar motors assemble wider stator rings using a conserved protein scaffold. Proc. Natl. Acad. Sci. USA 2016, 113, E1917–E1926. [Google Scholar] [CrossRef] [Green Version]
  22. Flint, A.; Sun, Y.Q.; Butcher, J.; Stahl, M.; Huang, H.; Stintzi, A. Phenotypic screening of a targeted mutant library reveals Campylobacter jejuni defenses against oxidative stress. Infect. Immun. 2014, 82, 2266–2275. [Google Scholar] [CrossRef] [Green Version]
  23. Sacher, J.C.; Flint, A.; Butcher, J.; Blasdel, B.; Reynolds, H.M.; Lavigne, R.; Stintzi, A.; Szymanski, C.M. Transcriptomic analysis of the Campylobacter jejuni response to T4-like phage NCTC 12673 infection. Viruses 2018, 10, 332. [Google Scholar] [CrossRef] [Green Version]
  24. Sacher, J.C.; Shajahan, A.; Butcher, J.; Patry, R.T.; Flint, A.; Hendrixson, D.R.; Szymanski, C.M. Binding of phage-encoded FlaGrab to motile Campylobacter jejuni flagella inhibits growth, downregulates energy metabolism, and requires specific flagellar glycans. Front. Microbiol. 2020, 11, 397. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Sørensen, M.C.H.; van Alphen, L.B.; Harboe, A.; Li, J.; Christensen, B.B.; Szymanski, C.M.; Brøndsted, L. Bacteriophage F336 recognizes the capsular phosphoramidate modification of Campylobacter jejuni NCTC11168. J. Bacteriol. 2011, 193, 6742–6749. [Google Scholar] [CrossRef] [Green Version]
  26. Parkhill, J.; Wren, B.W.; Mungall, K.; Ketley, J.M.; Churcher, C.; Basham, D.; Chillingworth, T.; Davies, R.M.; Feltwell, T.; Holroyd, S.; et al. The genome sequence of the food-borne pathogen Campylobacter jejuni reveals hypervariable sequences. Nature 2000, 403, 665–668. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Michael, F.S.; Szymanski, C.M.; Li, J.; Chan, K.H.; Khieu, N.H.; Larocque, S.; Wakarchuk, W.W.; Brisson, J.R.; Monteiro, M.A. The structures of the lipooligosaccharide and capsule polysaccharide of Campylobacter jejuni genome sequenced strain NCTC 11168. Eur. J. Biochem. 2002, 269, 5119–5136. [Google Scholar] [CrossRef] [PubMed]
  28. Irons, J.; Sacher, J.C.; Szymanski, C.M.; Downs, D.M. Cj1388 is a RidA homolog and is required for flagella biosynthesis and/or function in Campylobacter jejuni. Front. Microbiol. 2019, 10, 2058. [Google Scholar] [CrossRef] [PubMed]
  29. Javed, M.A.; van Alphen, L.B.; Sacher, J.; Ding, W.; Kelly, J.; Nargang, C.; Smith, D.F.; Cummings, R.D.; Szymanski, C.M. A receptor-binding protein of Campylobacter jejuni bacteriophage NCTC 12673 recognizes flagellin glycosylated with acetamidino-modified pseudaminic acid. Mol. Microbiol. 2015, 95, 101–115. [Google Scholar] [CrossRef] [PubMed]
  30. Javed, M.A.; Sacher, J.C.; van Alphen, L.B.; Patry, R.T.; Szymanski, C.M. A flagellar glycan-specific protein encoded by Campylobacter phages inhibits host cell growth. Viruses 2015, 7, 6661–6674. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Palyada, K.; Sun, Y.Q.; Flint, A.; Butcher, J.; Naikare, H.; Stintzi, A. Characterization of the oxidative stress stimulon and PerR regulon of Campylobacter jejuni. BMC Genom. 2009, 10, 481. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Dwivedi, R.; Nothaft, H.; Garber, J.; Szymanski, C.M. L-Fucose influences chemotaxis and biofilm formation in Campylobacter jejuni. Mol. Microbiol. 2016, 101, 575–589. [Google Scholar] [CrossRef] [Green Version]
  33. Grajewski, B.A.; Kusek, J.W.; Gelfand, H.M. Development of bacteriophage typing system for Campylobacter jejuni and Campylobacter coli. J. Clin. Microbiol. 1985, 22, 13–18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Frost, J.A.; Kramer, J.M.; Gillanders, S.A. Phage typing of Campylobacter jejuni and Campylobacter coli and its use as an adjunct to serotyping. Epidemiol. Infect. 1999, 123, 47–55. [Google Scholar] [CrossRef]
  35. Negretti, N.M.; Gourley, C.R.; Clair, G.; Adkins, J.N.; Konkel, M.E. The food-borne pathogen Campylobacter jejuni responds to the bile salt deoxycholate with countermeasures to reactive oxygen species. Sci. Rep. 2017, 7, 1–11. [Google Scholar] [CrossRef] [PubMed]
  36. Baldvinsson, S.B.; Holst Sørensen, M.C.; Vegge, C.S.; Clokie, M.R.J.; Brøndsted, L. Campylobacter jejuni motility is required for infection of the flagellotropic bacteriophage F341. Appl. Environ. Microbiol. 2014, 80, 7096–7106. [Google Scholar] [CrossRef] [Green Version]
  37. Palyada, K.; Threadgill, D.; Stintzi, A. Iron acquisition and regulation in Campylobacter jejuni. J. Bacteriol. 2004, 186, 4714–4729. [Google Scholar] [CrossRef] [Green Version]
  38. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
  39. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 1–21. [Google Scholar] [CrossRef] [Green Version]
  40. Sørensen, M.C.; Gencay, Y.E.; Brøndsted, L. Methods for Initial Characterization of Campylobacter jejuni Bacteriophages. Methods Mol. Biol. 2017, 1512, 91–105. [Google Scholar] [CrossRef]
  41. McNally, D.J.; Schoenhofen, I.C.; Houliston, R.S.; Khieu, N.H.; Whitfield, D.M.; Logan, S.M.; Brisson, J.R. CMP-Pseudaminic Acid is a Natural Potent Inhibitor of PseB, the First Enzyme of the Pseudaminic Acid Pathway in Campylobacter jejuni and Helicobacter pylori. ChemMedChem Chem. Enabling Drug Discov. 2008, 3, 55–59. [Google Scholar]
  42. Szymanski, C.M.; Logan, S.M.; Linton, D.; Wren, B.W. Campylobacter—A tale of two protein glycosylation systems. Trends Microbiol. 2003, 11, 233–238. [Google Scholar] [CrossRef]
  43. Hofreuter, D.; Tsai, J.; Watson, R.O.; Novik, V.; Altman, B.; Benitez, M.; Galán, J.E. Unique features of a highly pathogenic Campylobacter jejuni strain. Infect. Immun. 2006, 74, 4694–4707. [Google Scholar] [CrossRef] [Green Version]
  44. Logan, S.M. Flagellar glycosylation—A new component of the motility repertoire? Microbiology 2006, 152, 1249–1262. [Google Scholar] [CrossRef] [Green Version]
  45. Hooton, S.; Connerton, I.F. Campylobacter jejuni acquire new host-derived CRISPR spacers when in association with bacteriophages harboring a CRISPR-like Cas4 protein. Front. Microbiol. 2015, 5, 744. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Radomska, K.A.; Ordoñez, S.R.; Wösten, M.M.; Wagenaar, J.A.; van Putten, J.P. Feedback control of Campylobacter jejuni flagellin levels through reciprocal binding of FliW to flagellin and the global regulator CsrA. Mol. Microbiol. 2016, 102, 207–220. [Google Scholar] [CrossRef]
  47. Dugar, G.; Svensson, S.L.; Bischler, T.; Wäldchen, S.; Reinhardt, R.; Sauer, M.; Sharma, C.M. The CsrA-FliW network controls polar localization of the dual-function flagellin mRNA in Campylobacter jejuni. Nat. Commun. 2016, 7, 1–18. [Google Scholar] [CrossRef] [Green Version]
  48. Boll, J.M.; Hendrixson, D.R. A regulatory checkpoint during flagellar biogenesis in Campylobacter jejuni initiates signal transduction to activate transcription of flagellar genes. MBio 2013, 4, e00432-3. [Google Scholar] [CrossRef] [Green Version]
  49. Gilbreath, J.J.; Cody, W.L.; Merrell, D.S.; Hendrixson, D.R. Change is good: Variations in common biological mechanisms in the epsilonproteobacterial genera Campylobacter and Helicobacter. Microbiol. Mol. Biol. Rev. 2011, 75, 84–132. [Google Scholar] [CrossRef] [Green Version]
  50. Guerry, P. Campylobacter flagella: Not just for motility. Trends Microbiol. 2007, 15, 456–461. [Google Scholar] [CrossRef]
  51. Carrillo, C.D.; Taboada, E.; Nash, J.H.; Lanthier, P.; Kelly, J.; Lau, P.C.; Szymanski, C.M. Genome-wide expression analyses of Campylobacter jejuni NCTC11168 reveals coordinate regulation of motility and virulence by flhA. J. Biol. Chem. 2004, 279, 20327–20338. [Google Scholar] [CrossRef] [Green Version]
  52. Chatterjee, A.; Willett, J.L.; Nguyen, U.T.; Monogue, B.; Palmer, K.L.; Dunny, G.M.; Duerkop, B.A. Parallel genomics uncover novel enterococcal-bacteriophage interactions. MBio 2020, 11, e03120-19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Kim, J.C.; Oh, E.; Kim, J.; Jeon, B. Regulation of oxidative stress resistance in Campylobacter jejuni, a microaerophilic foodborne pathogen. Front. Microbiol. 2016, 6, 751. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Benoit, S.L.; Bayyareddy, K.; Mahawar, M.; Sharp, J.S.; Maier, R.J. Alkyl hydroperoxide reductase repair by Helicobacter pylori methionine sulfoxide reductase. J. Bacteriol. 2015, 195, 5396–5401. [Google Scholar] [CrossRef] [Green Version]
  55. Benoit, S.L.; Maier, R.J. Helicobacter Catalase Devoid of Catalytic Activity Protects the Bacterium against Oxidative Stress. J. Biol. Chem. 2016, 291, 23366–23373. [Google Scholar] [CrossRef] [Green Version]
  56. Hofreuter, D.; Novik, V.; Galán, J.E. Metabolic diversity in Campylobacter jejuni enhances specific tissue colonization. Cell Host Microbe 2008, 4, 425–433. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Ezraty, B.; Gennaris, A.; Barras, F.; Collet, J.F. Oxidative stress, protein damage and repair in bacteria. Nat. Rev. Microbiol. 2017, 15, 385–396. [Google Scholar] [CrossRef] [PubMed]
  58. Lango-Scholey, L.; Aidley, J.; Woodacre, A.; Jones, M.A.; Bayliss, C.D. High throughput method for analysis of repeat number for 28 phase variable loci of Campylobacter jejuni strain NCTC11168. PLoS ONE 2016, 11, e0159634. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Sørensen, M.C.H.; Vitt, A.; Neve, H.; Soverini, M.; Ahern, S.J.; Klumpp, J.; Brøndsted, L. Campylobacter phages use hypermutable polyG tracts to create phenotypic diversity and evade bacterial resistance. Cell Rep. 2021, 35, 109214. [Google Scholar] [CrossRef]
  60. Aidley, J.; Holst Sørensen, M.C.; Bayliss, C.D.; Brøndsted, L. Phage exposure causes dynamic shifts in the expression states of specific phase-variable genes of Campylobacter jejuni. Microbiology 2017, 163, 911–919. [Google Scholar] [CrossRef]
  61. Tormo-Más, M.Á.; Mir, I.; Shrestha, A.; Tallent, S.M.; Campoy, S.; Lasa, Í.; Penadés, J.R. Moonlighting bacteriophage proteins derepress staphylococcal pathogenicity islands. Nature 2010, 465, 779–782. [Google Scholar] [CrossRef] [Green Version]
  62. Siringan, P.; Connerton, P.L.; Cummings, N.J.; Connerton, I.F. Alternative bacteriophage life cycles: The carrier state of Campylobacter jejuni. Open Biol. 2014, 4, 130200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Sacher, J.C. Insights into the Role of the Flagellar Glycosylation System in Campylobacter jejuni Phage-Host Interactions. Ph.D. Thesis, University of Alberta, Edmonton, AB, Canada, 2018. [Google Scholar]
Figure 1. (A) Percent efficiency of plating (EOP) of phage NCTC 12673 on C. jejuni NCTC 11168 wild type (WT) and isogenic mutants in genes encoding enzymes involved in pseudaminic acid biosynthesis (pse), flagellar motor function (mot), and kpsM (capsular polysaccharide biosynthesis; included as a negative control). Bars represent the mean of three to six biological replicates. One-way ANOVA showed all mutant strains to be statistically significantly different from wild type (p < 0.0001). The results of an unpaired t-test for ΔpseH vs. pseH complement strain is indicated by asterisks (p = 0.0044). (B) Schematic of the pseudaminic acid pathway for flagellar O-glycosylation in C. jejuni, which is required for flagellar motility [20]. Dotted line shows PseB is inhibited by CMP-Pse [41]. OM: outer membrane, PG: peptidoglycan, IM: inner membrane, Pse5Ac7Ac: pseudaminic acid, Pse5Ac7Am: acetamidino-modified pseudaminic acid.
Figure 1. (A) Percent efficiency of plating (EOP) of phage NCTC 12673 on C. jejuni NCTC 11168 wild type (WT) and isogenic mutants in genes encoding enzymes involved in pseudaminic acid biosynthesis (pse), flagellar motor function (mot), and kpsM (capsular polysaccharide biosynthesis; included as a negative control). Bars represent the mean of three to six biological replicates. One-way ANOVA showed all mutant strains to be statistically significantly different from wild type (p < 0.0001). The results of an unpaired t-test for ΔpseH vs. pseH complement strain is indicated by asterisks (p = 0.0044). (B) Schematic of the pseudaminic acid pathway for flagellar O-glycosylation in C. jejuni, which is required for flagellar motility [20]. Dotted line shows PseB is inhibited by CMP-Pse [41]. OM: outer membrane, PG: peptidoglycan, IM: inner membrane, Pse5Ac7Ac: pseudaminic acid, Pse5Ac7Am: acetamidino-modified pseudaminic acid.
Viruses 13 01955 g001
Figure 2. Percent phage NCTC 12673 adsorption to C. jejuni NCTC 11168 wild type and ΔpseC, ΔpseF ΔpseG, ΔpseH, ΔpseH/ + pseH, and ΔkpsM mutants. Bars represent the mean of three biological replicates, and error bars represent standard deviations. Statistical significance as determined by one-way ANOVA is indicated by asterisks as follows: p < 0.05 (*) and p < 0.0001 (****). All other mutant strains were shown to have no significant difference from wild type.
Figure 2. Percent phage NCTC 12673 adsorption to C. jejuni NCTC 11168 wild type and ΔpseC, ΔpseF ΔpseG, ΔpseH, ΔpseH/ + pseH, and ΔkpsM mutants. Bars represent the mean of three biological replicates, and error bars represent standard deviations. Statistical significance as determined by one-way ANOVA is indicated by asterisks as follows: p < 0.05 (*) and p < 0.0001 (****). All other mutant strains were shown to have no significant difference from wild type.
Viruses 13 01955 g002
Figure 3. Volcano plots showing the most significantly differentially expressed genes between C. jejuni NCTC 11168 wild type, ΔpseC and ΔpseF cells. Significantly differentially expressed genes between ΔpseC and wild type (A), between ΔpseF and wild type (B) and between ΔpseF and ΔpseC (C) are shown. The negative log of the false discovery rate (FDR)-adjusted p-value vs. the log2 fold change between the conditions indicated is plotted for each gene. Genes with a fold change +/− 1.5 and FDR-corrected p-value < 0.05 were considered differentially expressed. Selected significantly differentially expressed genes are represented by different colours and symbols according to their predicted or known function.
Figure 3. Volcano plots showing the most significantly differentially expressed genes between C. jejuni NCTC 11168 wild type, ΔpseC and ΔpseF cells. Significantly differentially expressed genes between ΔpseC and wild type (A), between ΔpseF and wild type (B) and between ΔpseF and ΔpseC (C) are shown. The negative log of the false discovery rate (FDR)-adjusted p-value vs. the log2 fold change between the conditions indicated is plotted for each gene. Genes with a fold change +/− 1.5 and FDR-corrected p-value < 0.05 were considered differentially expressed. Selected significantly differentially expressed genes are represented by different colours and symbols according to their predicted or known function.
Viruses 13 01955 g003
Figure 4. Percent efficiency of plating (EOP) of phages NCTC 12673 and MutC on C. jejuni NCTC 11168 wild type (WT) and isogenic mutants in pseudaminic acid biosynthesis (pseCFGH), flagellar motor (motAB), and oxidative stress defense genes (katA, ahpC, sodB), as well as WT cells in the presence of 0.05% deoxycholate (DC), an oxidative stress-inducing agent for Campylobacter [35]. Bars represent the mean of three to six biological replicates (with the exception of WT+DC, which represents the mean of two biological replicates). Statistical significance between NCTC 12673 and MutC phage for each strain was determined by unpaired t-test and is indicated by asterisks as follows: p < 0.05 (*), p < 0.01 (**), p < 0.001 (***). n.s: no significant difference between NCTC 12673 and MutC.
Figure 4. Percent efficiency of plating (EOP) of phages NCTC 12673 and MutC on C. jejuni NCTC 11168 wild type (WT) and isogenic mutants in pseudaminic acid biosynthesis (pseCFGH), flagellar motor (motAB), and oxidative stress defense genes (katA, ahpC, sodB), as well as WT cells in the presence of 0.05% deoxycholate (DC), an oxidative stress-inducing agent for Campylobacter [35]. Bars represent the mean of three to six biological replicates (with the exception of WT+DC, which represents the mean of two biological replicates). Statistical significance between NCTC 12673 and MutC phage for each strain was determined by unpaired t-test and is indicated by asterisks as follows: p < 0.05 (*), p < 0.01 (**), p < 0.001 (***). n.s: no significant difference between NCTC 12673 and MutC.
Viruses 13 01955 g004
Table 1. List of strains and phages used in this study.
Table 1. List of strains and phages used in this study.
StrainDescription/PhenotypeReference
C. jejuni NCTC 11168Human enteropathy isolate, capsular, motile[26]
C. jejuni NCTC 11168∆motANon-motile (paralyzed flagella)[24]
C. jejuni NCTC 11168∆motBNon-motile (paralyzed flagella)[24]
C. jejuni NCTC 11168∆kpsMAcapsular[27]
C. jejuni NCTC 11168∆pseCNon-motile (aflagellate)[28]
C. jejuni NCTC 11168∆pseHNon-motile (aflagellate)[29]
C. jejuni NCTC 11168∆pseH/+pseHMotileThis work
C. jejuni NCTC 11168∆pseGNon-motile (aflagellate)[30]
C. jejuni NCTC 11168∆pseFNon-motile (aflagellate)This work
C. jejuni NCTC 11168∆katAHypersensitive to oxidative stress (lacks catalase)[31]
C. jejuni NCTC 11168∆ahpCHypersensitive to oxidative stress (lacks alkyl hydroxyperoxide reductase)[31]
C. jejuni NCTC 11168∆sodBHypersensitive to oxidative stress (lacks superoxide dismutase)[31]
C. jejuni NCTC 11168∆flaANon-motile (aflagellate)[32]
NCTC 12673UK phage typing scheme phage 1[18,33]
MutCSpontaneous variant of NCTC 12673This work
Table 2. Variant frequency table showing amino acid-level genetic differences between NCTC 12673 and MutC phages. For each gene displaying one or more variants between the two phages, the variants are listed, along with nucleotide change(s) relative to the NCTC 12673 reference strain in GenBank (accession: GU296433) [18], the corresponding amino acid change, the variant frequency, and the corresponding phage genome for each. (A)7 → (A)6: AAAAAAA → AAAAAA. gp114/gp115: change observed in the region where these genes overlap.
Table 2. Variant frequency table showing amino acid-level genetic differences between NCTC 12673 and MutC phages. For each gene displaying one or more variants between the two phages, the variants are listed, along with nucleotide change(s) relative to the NCTC 12673 reference strain in GenBank (accession: GU296433) [18], the corresponding amino acid change, the variant frequency, and the corresponding phage genome for each. (A)7 → (A)6: AAAAAAA → AAAAAA. gp114/gp115: change observed in the region where these genes overlap.
Gene.Predicted FunctionVariantPhagePosition within GenomeNucleotide ChangeAmino Acid ChangeVariant Frequency
gp041Gp6 baseplate wedge subunit1MutC31,844A → GL → S72.70%
2MutC31,885C → TM → I25.60%
gp047/flagrabFlaGrab, flagellar glycan binding protein1MutC40,821(A)7 → (A)6Frame Shift68.80%
2NCTC 1267340,935(A)7 → (A)6Frame Shift33.20%
MutC40,935(A)7 → (A)6Frame Shift25.30%
gp058Hef59 homing endonuclease1MutC50,606AA → GGE → G69.8%
gp114/gp115Hypothetical protein1MutC100,737(CC)4 → (CC)5Frame Shift37.00%
NCTC 12673100,737(C)9 → (C)10Frame Shift31.20%
MutC100,737(C)9 → (C)10Frame Shift43.70%
gp116Hypothetical protein1NCTC 12673101,692(C)11 → (C)12Frame Shift28.90%
2NCTC 12673101,691(C)11 → (C)10Frame Shift27.10%
MutC101,691(C)11 → (C)10Frame Shift62.80%
3MutC101,202G → AT → M68.70%
4NCTC 12673101,112G → TT → N73.80%
gp167Hef168 homing endonuclease1MutC132,316C → TD → N64.00%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sacher, J.C.; Javed, M.A.; Crippen, C.S.; Butcher, J.; Flint, A.; Stintzi, A.; Szymanski, C.M. Reduced Infection Efficiency of Phage NCTC 12673 on Non-Motile Campylobacter jejuni Strains Is Related to Oxidative Stress. Viruses 2021, 13, 1955. https://doi.org/10.3390/v13101955

AMA Style

Sacher JC, Javed MA, Crippen CS, Butcher J, Flint A, Stintzi A, Szymanski CM. Reduced Infection Efficiency of Phage NCTC 12673 on Non-Motile Campylobacter jejuni Strains Is Related to Oxidative Stress. Viruses. 2021; 13(10):1955. https://doi.org/10.3390/v13101955

Chicago/Turabian Style

Sacher, Jessica C., Muhammad Afzal Javed, Clay S. Crippen, James Butcher, Annika Flint, Alain Stintzi, and Christine M. Szymanski. 2021. "Reduced Infection Efficiency of Phage NCTC 12673 on Non-Motile Campylobacter jejuni Strains Is Related to Oxidative Stress" Viruses 13, no. 10: 1955. https://doi.org/10.3390/v13101955

APA Style

Sacher, J. C., Javed, M. A., Crippen, C. S., Butcher, J., Flint, A., Stintzi, A., & Szymanski, C. M. (2021). Reduced Infection Efficiency of Phage NCTC 12673 on Non-Motile Campylobacter jejuni Strains Is Related to Oxidative Stress. Viruses, 13(10), 1955. https://doi.org/10.3390/v13101955

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop