High Throughput Method for Analysis of Repeat Number for 28 Phase Variable Loci of Campylobacter jejuni Strain NCTC11168

Novel method for prediction of combinatorial phase-variable gene expression states

Short sequence repeat mediated phase variation results in diverse phenotype presentation in many bacteria including Campylobacter and Neisseria species. Current methods for identifying the expression states of phase-variable genes involve taking a high number of single colonies. This approach is subject to bias, sampling effects and high workloads that reduce the ability to perform intermediary sampling. The use of high concentration colony sweeps provides a work around but reduces the resolution of combinatorial expression profiles (termed phasotypes). A parsimonious approach combining both single colony and sweep data was developed to overcome these limitations. The critical methodological advance is the use of an algorithm that utilises the experimental data from the two sample types and a parsimonious, iterative mathematical analysis that outputs the phasotype distribution with the highest likelihood of underpinning the experimental data sets. The advantages of this unified method are increased resolution and accuracy of gene expression state combinations as compared to conventional single colony sampling, reduced requirement for sampling large numbers of colonies leading to reduced costs, and a higher capacity for collecting samples and replicates.•Inputting of sweep and single colony data into an algorithm for a rapid determination of the combinatorial phase variation states (phasotypes) for repeat-mediated phase-variable bacterial genes•This method reduces the number of single colony samples required to produce accurate estimates of phasotypes•This method will reduce the costs of phasotype analyses and increase potential to analyse more time points or sample sites leading to an improved understanding of how phase variation contributes to bacterial host persistence and the ability to cause disease.

Phase variation in the glycosyltransferase genes of Pasteurella multocida associated with outbreaks of fowl cholera on free-range layer farms

Fowl cholera caused by Pasteurella multocida has re-emerged in Australian poultry production since the increasing adoption of free-range production systems. Currently, autogenous killed whole-cell vaccines prepared from the isolates previously obtained from each farm are the main preventative measures used. In this study, we use whole-genome sequencing and phylogenomic analysis to investigate outbreak dynamics, as well as monitoring and comparing the variations in the lipopolysaccharide (LPS) outer core biosynthesis loci of the outbreak and vaccine strains. In total, 73 isolates from two different free-range layer farms were included. Our genomic analysis revealed that all investigated isolates within the two farms (layer A and layer B) carried LPS type L3, albeit with a high degree of genetic diversity between them. Additionally, the isolates belonged to five different sequence types (STs), with isolates belonging to ST9 and ST20 being the most prevalent. The isolates carried ST-specific mutations within their LPS type L3 outer core biosynthesis loci, including frameshift mutations in the outer core heptosyltransferase gene (htpE) (ST7 and ST274) or galactosyltransferase gene (gatG) (ST20). The ST9 isolates could be separated into three groups based on their LPS outer core biosynthesis loci sequences, with evidence for potential phase variation mechanisms identified. The potential phase variation mechanisms included a tandem repeat insertion in natC and a single base deletion in a homopolymer region of gatG. Importantly, our results demonstrated that two of the three ST9 groups shared identical rep-PCR (repetitive extragenic palindromic PCR) patterns, while carrying differences in their LPS outer core biosynthesis loci region. In addition, we found that ST9 isolates either with or without the natC tandem repeat insertion were both associated with a single outbreak, which would indicate the importance of screening more than one isolate within an outbreak. Our results strongly suggest the need for a metagenomics culture-independent approach, as well as a genetic typing scheme for LPS, to ensure an appropriate vaccine strain with a matching predicted LPS structure is used.

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

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.

Stabilizing Genetically Unstable Simple Sequence Repeats in the Campylobacter jejuni Genome by Multiplex Genome Editing: a Reliable Approach for Delineating Multiple Phase-Variable Genes

Hypermutable simple sequence repeats (SSRs) are major drivers of phase variation in Campylobacter jejuni. The presence of multiple SSR-mediated phase-variable genes encoding enzymes that modify surface structures, including capsular polysaccharide (CPS) and lipooligosaccharide (LOS), generates extreme cell surface diversity within bacterial populations, thereby promoting adaptation to selective pressures in host environments. Therefore, genetically controlling SSR-mediated phase variation can be important for achieving stable and reproducible research on C. jejuni. Here, we show that natural "cotransformation" is an effective method for C. jejuni genome editing. Cotransformation is a trait of naturally competent bacteria that causes uptake/integration of multiple different DNA molecules, which has been recently adapted to multiplex genome editing by natural transformation (MuGENT), a method for introducing multiple mutations into the genomes of these bacteria. We found that cotransformation efficiently occurred in C. jejuni. To examine the feasibility of MuGENT in C. jejuni, we "locked" different polyG SSR tracts in strain NCTC11168 (which are located in the biosynthetic CPS/LOS gene clusters) into either the ON or OFF configurations. This approach, termed "MuGENT-SSR," enabled the generation of all eight edits within 2 weeks and the identification of a phase-locked strain with a highly stable type of Penner serotyping, a CPS-based serotyping scheme. Furthermore, extensive genome editing of this strain by MuGENT-SSR identified a phase-variable gene that determines the Penner serotype of NCTC11168. Thus, MuGENT-SSR provides a platform for genetic and phenotypic engineering of genetically unstable C. jejuni, making it a reliable approach for elucidating the mechanisms underlying phase-variable expression of specific phenotypes. IMPORTANCE Campylobacter jejuni is the leading bacterial cause of foodborne gastroenteritis in developed countries and occasionally progresses to the autoimmune disease Guillain-Barré syndrome. A relatively large number of hypermutable simple sequence repeat (SSR) tracts in the C. jejuni genome markedly decreases its phenotypic stability through reversible changes in the ON or OFF expression states of the genes in which they reside, a phenomenon called phase variation. Thus, controlling SSR-mediated phase variation can be important for achieving stable and reproducible research on C. jejuni. In this study, we developed a feasible and effective approach for genetically manipulate multiple SSR tracts in the C. jejuni genome using natural cotransformation, a trait of naturally transformable bacterial species that causes the uptake and integration of multiple different DNA molecules. This approach will greatly help to improve the genetic and phenotypic stability of C. jejuni to enable diverse applications in research and development.

Binding of Phage-Encoded FlaGrab to Motile Campylobacter jejuni Flagella Inhibits Growth, Downregulates Energy Metabolism, and Requires Specific Flagellar Glycans

Many bacterial pathogens display glycosylated surface structures that contribute to virulence, and targeting these structures is a viable strategy for pathogen control. The foodborne pathogen Campylobacter jejuni expresses a vast diversity of flagellar glycans, and flagellar glycosylation is essential for its virulence. Little is known about why C. jejuni encodes such a diverse set of flagellar glycans, but it has been hypothesized that evolutionary pressure from bacteriophages (phages) may have contributed to this diversity. However, interactions between Campylobacter phages and host flagellar glycans have not been characterized in detail. Previously, we observed that Gp047 (now renamed FlaGrab), a conserved Campylobacter phage protein, binds to C. jejuni flagella displaying the nine-carbon monosaccharide 7-acetamidino-pseudaminic acid, and that this binding partially inhibits cell growth. However, the mechanism of this growth inhibition, as well as how C. jejuni might resist this activity, are not well-understood. Here we use RNA-Seq to show that FlaGrab exposure leads C. jejuni 11168 cells to downregulate expression of energy metabolism genes, and that FlaGrab-induced growth inhibition is dependent on motile flagella. Our results are consistent with a model whereby FlaGrab binding transmits a signal through flagella that leads to retarded cell growth. To evaluate mechanisms of FlaGrab resistance in C. jejuni, we characterized the flagellar glycans and flagellar glycosylation loci of two C. jejuni strains naturally resistant to FlaGrab binding. Our results point toward flagellar glycan diversity as the mechanism of resistance to FlaGrab. Overall, we have further characterized the interaction between this phage-encoded flagellar glycan-binding protein and C. jejuni, both in terms of mechanism of action and mechanism of resistance. Our results suggest that C. jejuni encodes as-yet unidentified mechanisms for generating flagellar glycan diversity, and point to phage proteins as exciting lenses through which to study bacterial surface glycans.

Mode and Tempo of Microsatellite Length Change in a Malaria Parasite Mutation Accumulation Experiment

Malaria parasites have small extremely AT-rich genomes: microsatellite repeats (1–9 bp) comprise 11% of the genome and genetic variation in natural populations is dominated by repeat changes in microsatellites rather than point mutations. This experiment was designed to quantify microsatellite mutation patterns in Plasmodium falciparum. We established 31 parasite cultures derived from a single parasite cell and maintained these for 114–267 days with frequent reductions to a single cell, so parasites accumulated mutations during ∼13,207 cell divisions. We Illumina sequenced the genomes of both progenitor and end-point mutation accumulation (MA) parasite lines in duplicate to validate stringent calling parameters. Microsatellite calls were 99.89% (GATK), 99.99% (freeBayes), and 99.96% (HipSTR) concordant in duplicate sequence runs from independent sequence libraries, whereas introduction of microsatellite mutations into the reference genome revealed a low false negative calling rate (0.68%). We observed 98 microsatellite mutations. We highlight several conclusions: microsatellite mutation rates (3.12 × 10⁻⁷ to 2.16 × 10⁻⁸/cell division) are associated with both repeat number and repeat motif like other organisms studied. However, 41% of changes resulted from loss or gain of more than one repeat: this was particularly true for long repeat arrays. Unlike other eukaryotes, we found no insertions or deletions that were not associated with repeats or homology regions. Overall, microsatellite mutation rates are among the lowest recorded and comparable to those in another AT-rich protozoan (Dictyostelium). However, a single infection (>10¹¹ parasites) will still contain over 2.16 × 10³ to 3.12 × 10⁴ independent mutations at any single microsatellite locus.

Random sorting of Campylobacter jejuni phase variants due to a narrow bottleneck during colonization of broiler chickens

Phase variation (PV), involving stochastic switches in gene expression, is exploited by the human pathogen Campylobacter jejuni to adapt to different environmental and host niches. Phase-variable genes of C. jejuni modulate expression of multiple surface determinants, and hence may influence host colonization. Population bottlenecks can rapidly remove the diversity generated by PV, and strict single-cell bottlenecks can lead to propagation of PV states with highly divergent phenotypes. Using a combination of high-throughput fragment size analysis and comparison with in vivo and in silico bottleneck models, we have characterized a narrow population bottleneck during the experimental colonization of broiler chickens with C. jejuni strain 81-176. We identified high levels of variation in five PV genes in the inoculum, and subsequently, massively decreased population diversity following colonization. Each bird contained a dominant five-gene phasotype that was present in the inoculum indicative of random sorting through a narrow, non-selective bottleneck during colonization. These results are evidence of the potential for confounding effects of PV on in vivo studies of Campylobacter colonization factors and poultry vaccine studies. Our results are also an argument for population bottlenecks as mediators of stochastic variability in the propensity to survive through the food chain and cause clinical human disease.

PhasomeIt: an 'omics' approach to cataloguing the potential breadth of phase variation in the genus Campylobacter

Hypermutable simple sequence repeats (SSRs) are drivers of phase variation (PV) whose stochastic, high-frequency, reversible switches in gene expression are a common feature of several pathogenic bacterial species, including the human pathogen Campylobacter jejuni. Here we examine the distribution and conservation of known and putative SSR-driven phase variable genes – the phasome – in the genus Campylobacter. PhasomeIt, a new program, was specifically designed for rapid identification of SSR-mediated PV. This program detects the location, type and repeat number of every SSR. Each SSR is linked to a specific gene and its putative expression state. Other outputs include conservation of SSR-driven phase-variable genes and the ‘core phasome’ – the minimal set of PV genes in a phylogenetic grouping. Analysis of 77 complete Campylobacter genome sequences detected a ‘core phasome’ of conserved PV genes in each species and a large number of rare PV genes with few, or no, homologues in other genome sequences. Analysis of a set of partial genome sequences, with food-chain-associated metadata, detected evidence of a weak link between phasome and source host for disease-causing isolates of sequence type (ST)-828 but not the ST-21 or ST-45 complexes. Investigation of the phasomes in the genus Campylobacter provided evidence of overlapping but distinctive mechanisms of PV-mediated adaptation to specific niches. This suggests that the phasome could be involved in host adaptation and spread of campylobacters. Finally, this tool is malleable and will have utility for studying the distribution and genic effects of other repetitive elements in diverse bacterial species.

C-Type Lectin Receptor (CLR)-Fc Fusion Proteins As Tools to Screen for Novel CLR/Bacteria Interactions: An Exemplary Study on Preselected Campylobacter jejuni Isolates

C-type lectin receptors (CLRs) are carbohydrate-binding receptors that recognize their ligands often in a Ca²⁺-dependent manner. Upon ligand binding, myeloid CLRs in innate immunity trigger or inhibit a variety of signaling pathways, thus initiating or modulating effector functions such as cytokine production, phagocytosis, and antigen presentation. CLRs bind to various pathogens, including viruses, fungi, parasites, and bacteria. The bacterium Campylobacter jejuni (C. jejuni) is a very frequent Gram-negative zoonotic pathogen of humans, causing severe intestinal symptoms. Interestingly, C. jejuni expresses several glycosylated surface structures, for example, the capsular polysaccharide (CPS), lipooligosaccharide (LOS), and envelope proteins. This “Methods” paper describes applications of CLR–Fc fusion proteins to screen for yet unknown CLR/bacteria interactions using C. jejuni as an example. ELISA-based detection of CLR/bacteria interactions allows a first prescreening that is further confirmed by flow cytometry-based binding analysis and visualized using confocal microscopy. By applying these methods, we identified Dectin-1 as a novel CLR recognizing two selected C. jejuni isolates with different LOS and CPS genotypes. In conclusion, the here-described applications of CLR–Fc fusion proteins represent useful methods to screen for and identify novel CLR/bacteria interactions.

Coadministration of the Campylobacter jejuni N-Glycan-Based Vaccine with Probiotics Improves Vaccine Performance in Broiler Chickens

Source attribution studies report that the consumption of contaminated poultry is the primary source for acquiring human campylobacteriosis. Oral administration of an engineered Escherichia coli strain expressing the Campylobacter jejuni N-glycan reduces bacterial colonization in specific-pathogen-free leghorn chickens, but only a fraction of birds respond to vaccination. Optimization of the vaccine for commercial broiler chickens has great potential to prevent the entry of the pathogen into the food chain. Here, we tested the same vaccination approach in broiler chickens and observed similar efficacies in pathogen load reduction, stimulation of the host IgY response, the lack of C. jejuni resistance development, uniformity in microbial gut composition, and the bimodal response to treatment. Gut microbiota analysis of leghorn and broiler vaccine responders identified one member of Clostridiales cluster XIVa, Anaerosporobacter mobilis, that was significantly more abundant in responder birds. In broiler chickens, coadministration of the live vaccine with A. mobilis or Lactobacillus reuteri, a commonly used probiotic, resulted in increased vaccine efficacy, antibody responses, and weight gain. To investigate whether the responder-nonresponder effect was due to the selection of a C. jejuni "supercolonizer mutant" with altered phase-variable genes, we analyzed all poly(G)-containing loci of the input strain compared to nonresponder colony isolates and found no evidence of phase state selection. However, untargeted nuclear magnetic resonance (NMR)-based metabolomics identified a potential biomarker negatively correlated with C. jejuni colonization levels that is possibly linked to increased microbial diversity in this subgroup. The comprehensive methods used to examine the bimodality of the vaccine response provide several opportunities to improve the C. jejuni vaccine and the efficacy of any vaccination strategy.IMPORTANCE Campylobacter jejuni is a common cause of human diarrheal disease worldwide and is listed by the World Health Organization as a high-priority pathogen. C. jejuni infection typically occurs through the ingestion of contaminated chicken meat, so many efforts are targeted at reducing C. jejuni levels at the source. We previously developed a vaccine that reduces C. jejuni levels in egg-laying chickens. In this study, we improved vaccine performance in meat birds by supplementing the vaccine with probiotics. In addition, we demonstrated that C. jejuni colonization levels in chickens are negatively correlated with the abundance of clostridia, another group of common gut microbes. We describe new methods for vaccine optimization that will assist in improving the C. jejuni vaccine and other vaccines under development.

The impact of bottlenecks on microbial survival, adaptation, and phenotypic switching in host-pathogen interactions

Microbial pathogens and viruses can often maintain sufficient population diversity to evade a wide range of host immune responses. However, when populations experience bottlenecks, as occurs frequently during initiation of new infections, pathogens require specialized mechanisms to regenerate diversity. We address the evolution of such mechanisms, known as stochastic phenotype switches, which are prevalent in pathogenic bacteria. We analyze a model of pathogen diversification in a changing host environment that accounts for selective bottlenecks, wherein different phenotypes have distinct transmission probabilities between hosts. We show that under stringent bottlenecks, such that only one phenotype can initiate new infections, there exists a threshold stochastic switching rate below which all pathogen lineages go extinct, and above which survival is a near certainty. We determine how quickly stochastic switching rates can evolve by computing a fitness landscape for the evolutionary dynamics of switching rates, and analyzing its dependence on both the stringency of bottlenecks and the duration of within-host growth periods. We show that increasing the stringency of bottlenecks or decreasing the period of growth results in faster adaptation of switching rates. Our model provides strong theoretical evidence that bottlenecks play a critical role in accelerating the evolutionary dynamics of pathogens.

Phage exposure causes dynamic shifts in the expression states of specific phase-variable genes of Campylobacter jejuni

Phase variation (PV) creates phenotypic heterogeneity at high frequencies and in a reversible manner. This phenomenon allows bacteria to adapt to a variety of different environments and selective pressures. In Campylobacterjejuni this reversible adaptive process is mediated by mutations in homopolymeric G/C tracts. Many C. jejuni-specific phages are dependent on phase-variable surface structures for successful infection. We previously identified the capsular polysaccharide (CPS) moiety, MeOPN-GalfNAc, as a receptor for phage F336 and showed that phase-variable expression of the transferase for this CPS modification, cj1421, and two other phase-variable CPS genes generated phage resistance in C. jejuni. Here we investigate the population dynamics of C. jejuni NCTC11168 when exposed to phage F336 in vitro using a newly described method - the 28-locus-CJ11168 PV analysis. Dynamic switching was observed in the ON/OFF states of three phase-variable CPS genes, cj1421, cj1422 and cj1426, during phage F336 exposure, with the dominant phage-resistant phasotype differing between cultures. Although loss of the phage receptor was predominately observed, several other PV events also led to phage resistance, a phenomenon that increases the chance of phage-resistant subpopulations being present in any growing culture. No other PV genes were affected and exposure to phage F336 resulted in a highly specific response, only selecting for phase variants of cj1421, cj1422 and cj1426. In summary, C. jejuni may benefit from modification of the surface in multiple ways to inhibit or reduce phage binding, thereby ensuring the survival of the population when exposed to phages.

Nonselective Bottlenecks Control the Divergence and Diversification of Phase-Variable Bacterial Populations

Phase variation occurs in many pathogenic and commensal bacteria and is a major generator of genetic variability. A putative advantage of phase variation is to counter reductions in variability imposed by nonselective bottlenecks during transmission. Genomes of Campylobacter jejuni, a widespread food-borne pathogen, contain multiple phase-variable loci whose rapid, stochastic variation is generated by hypermutable simple sequence repeat tracts. These loci can occupy a vast number of combinatorial expression states (phasotypes) enabling populations to rapidly access phenotypic diversity. The imposition of nonselective bottlenecks can perturb the relative frequencies of phasotypes, changing both within-population diversity and divergence from the initial population. Using both in vitro testing of C. jejuni populations and a simple stochastic simulation of phasotype change, we observed that single-cell bottlenecks produce output populations of low diversity but with bimodal patterns of either high or low divergence. Conversely, large bottlenecks allow divergence only by accumulation of diversity, while interpolation between these extremes is observed in intermediary bottlenecks. These patterns are sensitive to the genetic diversity of initial populations but stable over a range of mutation rates and number of loci. The qualitative similarities of experimental and in silico modeling indicate that the observed patterns are robust and applicable to other systems where localized hypermutation is a defining feature. We conclude that while phase variation will maintain bacterial population diversity in the face of intermediate bottlenecks, narrow transmission-associated bottlenecks could produce host-to-host variation in bacterial phenotypes and hence stochastic variation in colonization and disease outcomes.

Transmission and within-host spread of pathogenic organisms are associated with selective and nonselective bottlenecks that significantly reduced population diversity. In several bacterial pathogens, hypermutable mechanisms have evolved that mediate high-frequency reversible switching of specific phenotypes, such as surface structures, and hence counteract bottleneck-associated reductions in population diversity. Here, we investigated how combinations of hypermutable simple sequence repeats interact with nonselective bottlenecks by using a stochastic computer model and experimental data for Campylobacter jejuni, a food-borne pathogen. We find that bottleneck size qualitatively alters the output populations, with large bottlenecks maintaining population diversity while small bottlenecks produce dramatic shifts in the prevalence of particular variants. We conclude that narrow bottlenecks are capable of producing host-to-host variation in repeat-controlled bacterial phenotypes, leading to a potential for stochastic person-to-person variations in disease outcome for C. jejuni and other organisms with similar hypermutable mechanisms.