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

Contrasting genes conferring short- and long-term biofilm adaptation in Listeria

Listeria monocytogenes is an opportunistic food-borne bacterium that is capable of infecting humans with high rates of hospitalization and mortality. Natural populations are genotypically and phenotypically variable, with some lineages being responsible for most human infections. The success of L. monocytogenes is linked to its capacity to persist on food and in the environment. Biofilms are an important feature that allow these bacteria to persist and infect humans, so understanding the genetic basis of biofilm formation is key to understanding transmission. We sought to investigate the biofilm-forming ability of L. monocytogenes by identifying genetic variation that underlies biofilm formation in natural populations using genome-wide association studies (GWAS). Changes in gene expression of specific strains during biofilm formation were then investigated using RNA sequencing (RNA-seq). Genetic variation associated with enhanced biofilm formation was identified in 273 genes by GWAS and differential expression in 220 genes by RNA-seq. Statistical analyses show that the number of overlapping genes flagged by either type of experiment is less than expected by random sampling. This novel finding is consistent with an evolutionary scenario where rapid adaptation is driven by variation in gene expression of pioneer genes, and this is followed by slower adaptation driven by nucleotide changes within the core genome.

Genomic Characterization of Campylobacter jejuni Associated with Perimyocarditis: A Family Case Report

Campylobacter spp. is the leading cause of foodborne gastrointestinal infections in humans worldwide. This study reports the first case of four family members who had contact with the same source of Campylobacter jejuni contamination with different results. Only the little siblings were infected by the same C. jejuni strain, but with different symptoms. Whereas the daughter was slightly affected with mild enteritis, the son suffered a longer campylobacteriosis followed with a perimyocarditis. This is the first case of the youngest patient affected by C. jejuni-related perimyocarditis published to date. The genomes of both strains were characterized by whole-genome sequencing and compared with the C. jejuni NCTC 11168 genome to gain insights into the molecular features that may be associated with perimyocarditis. Various comparison tools were used for the comparative genomics analysis, including the identification of virulence and antimicrobial resistance genes, phase variable (PV) genes, and single nucleotide polymorphisms (SNPs) identification. Comparisons of the strains identified 16 SNPs between them, which constituted small but significant changes mainly affecting the ON/OFF state of PV genes after passing through both hosts. These results suggest that PV occurs during human colonization, which modulates bacteria virulence through human host adaptation, which ultimately is related to complications after a campylobacteriosis episode depending on the host status. The findings highlight the importance of the relation between host and pathogen in severe complications of Campylobacter infections.

Identification of Novel Phage Resistance Mechanisms in Campylobacter jejuni by Comparative Genomics

Phages infecting Campylobacter jejuni are considered a promising intervention strategy at broiler farms, yet phage sensitivity of naturally occurring poultry isolates is not well studied. Here, we investigated phage sensitivity and identified resistance mechanisms of C. jejuni strains originating from Danish broilers belonging to the most prevalent MLST (ST) types. Determining plaque formation of 51 phages belonging to Fletchervirus or Firehammervirus showed that 21 out of 31 C. jejuni strains were susceptible to at least one phage. While C. jejuni ST-21 strains encoded the common phase variable O-methyl phosphoramidate (MeOPN) receptor of the Fletchervirus and were only infected by these phages, ST-45 strains did not encode this receptor and were exclusively infected by Firehammervirus phages. To identify internal phage resistance mechanism in ST-21 strains, we performed comparative genomics of two strains, CAMSA2002 sensitive to almost all Fletchervirus phages and CAMSA2038, resistant to all 51 phages. The strains encoded diverse clustered regularly interspaced short palindromic repeats (CRISPR) spacers but none matched the tested phages. Sequence divergence was also observed in a predicted SspE homolog and putative restriction modification systems including a methyl-specific McrBC endonuclease. Furthermore, when mcrB was deleted, CAMSA2038 became sensitive to 17 out of 43 phages, three being Firehammervirus phages that otherwise did not infect any ST-21 strains. Yet, 16 phages demonstrated significantly lower efficiencies of plating on the mcrB mutant suggesting additional resistance mechanism still restricting phage propagation in CAMSA2038. Thus, our work demonstrates that C. jejuni isolates originating from broilers may have acquired several resistance mechanisms to successfully prevent phage infection in their natural habitat.

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.

Phase Variation During Host Colonization and Invasion by Campylobacter jejuni and Other Campylobacter Species

Phase variation (PV) is a phenomenon common to a variety of bacterial species for niche adaption and survival in challenging environments. Among Campylobacter species, PV depends on the presence of intergenic and intragenic hypermutable G/C homopolymeric tracts. The presence of phase-variable genes is of especial interest for species that cause foodborne or zoonotic infections in humans. PV influences the formation and the structure of the lipooligosaccharide, flagella, and capsule in Campylobacter species. PV of components of these molecules is potentially important during invasion of host tissues, spread within hosts and transmission between hosts. Motility is a critical phenotype that is potentially modulated by PV. Variation in the status of the phase-variable genes has been observed to occur during colonization in chickens and mouse infection models. Interestingly, PV is also involved in bacterial survival of attack by bacteriophages even during chicken colonization. This review aims to explore and discuss observations of PV during model and natural infections by Campylobacter species and how PV may affect strategies for fighting infections by this foodborne pathogen.

Revisiting Campylobacter jejuni Virulence and Fitness Factors: Role in Sensing, Adapting, and Competing

Campylobacter jejuni is the leading cause of bacterial foodborne gastroenteritis world wide and represents a major public health concern. Over the past two decades, significant progress in functional genomics, proteomics, enzymatic-based virulence profiling (EBVP), and the cellular biology of C. jejuni have improved our basic understanding of this important pathogen. We review key advances in our understanding of the multitude of emerging virulence factors that influence the outcome of C. jejuni–mediated infections. We highlight, the spatial and temporal dynamics of factors that promote C. jejuni to sense, adapt and survive in multiple hosts. Finally, we propose cohesive research directions to obtain a comprehensive understanding of C. jejuni virulence mechanisms.

Selective and non-selective bottlenecks as drivers of the evolution of hypermutable bacterial loci

Bottlenecks reduce the size of the gene pool within populations of all life forms with implications for their subsequent survival. Here, we examine the effects of bottlenecks on bacterial commensal-pathogens during transmission between, and dissemination within, hosts. By reducing genetic diversity, bottlenecks may alter individual or population-wide adaptive potential. A diverse range of hypermutable mechanisms have evolved in infectious agents that allow for rapid generation of genetic diversity in specific genomic loci as opposed to the variability arising from increased genome-wide mutation rates. These localised hypermutable mechanisms include multi-gene phase variation (PV) of outer membrane components, multi-allele PV of restriction systems and recombination-driven antigenic variation. We review selected experimental and theoretical (mathematical) models pertaining to the hypothesis that localised hypermutation (LH) compensates for fitness losses caused by bottlenecks and discuss whether bottlenecks have driven the evolution of hypermutable loci.

Meningococcal core and accessory phasomes vary by clonal complex

Neisseria meningitidis is a Gram-negative human commensal pathogen, with extensive phenotypic plasticity afforded by phase-variable (PV) gene expression. Phase variation is a stochastic switch in gene expression from an ON to an OFF state, mediated by localized hypermutation of simple sequence repeats (SSRs). Circulating N. meningitidis clones vary in propensity to cause disease, with some clonal complexes (ccs) classified as hypervirulent and others as carriage-associated. We examined the PV gene repertoires, or phasome, of these lineages in order to determine whether phase variation contributes to disease propensity. We analysed 3328 genomes representative of nine circulating meningococcal ccs with PhasomeIt, a tool that identifies PV genes by the presence of SSRs and homologous gene clusters. The presence, absence and functions of all identified PV gene clusters were confirmed by annotation or blast searches within the Neisseria PubMLST database. While no significant differences were detected in the number of PV genes or the core, conserved phasome content between hypervirulent and carriage lineages, individual ccs exhibited major variations in PV gene numbers. Phylogenetic clusters produced by phasome or core genome analyses were similar, indicating co-evolution of PV genes with the core genome. While conservation of PV clusters is high, with 76 % present in all meningococcal isolates, maintenance of an SSR is variable, ranging from conserved in all isolates to present only in a single cc, indicating differing evolutionary trajectories for each lineage. Diverse functional groups of PV genes were present across the meningococcal lineages; however, the majority directly or indirectly influence bacterial surface antigens and could impact on future vaccine development. Finally, we observe that meningococci have open pan phasomes, indicating ongoing evolution of PV gene content and a significant potential for adaptive changes in this clinically relevant genus.

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.

Potentiation of Phase Variation in Multiple Outer-Membrane Proteins During Spread of the Hyperinvasive Neisseria meningitidis Serogroup W ST-11 Lineage

BACKGROUND

Since 2009, increases in the incidence of invasive meningococcal disease have occurred in the United Kingdom due to a sublineage of the Neisseria meningitidis serogroup W ST-11 clonal complex (hereafter, the "original UK strain"). In 2013, a descendent substrain (hereafter, the "2013 strain") became the dominant disease-causing variant. Multiple outer-membrane proteins of meningococci are subject to phase-variable switches in expression due to hypermutable simple-sequence repeats. We investigated whether alterations in phase-variable genes may have influenced the relative prevalence of the original UK and 2013 substrains, using multiple disease and carriage isolates.

METHODS

Repeat numbers were determined by either bioinformatics analysis of whole-genome sequencing data or polymerase chain reaction amplification and sizing of fragments from genomic DNA extracts. Immunoblotting and sequence-translation analysis was performed to identify expression states.

RESULTS

Significant increases in repeat numbers were detected between the original UK and 2013 strains in genes encoding PorA, NadA, and 2 Opa variants. Invasive and carriage isolates exhibited similar repeat numbers, but the absence of pilC gene expression was frequently associated with disease.

CONCLUSIONS

Elevated repeat numbers in outer-membrane protein genes of the 2013 strain are indicative of higher phase-variation rates, suggesting that rapid expansion of this strain was due to a heightened ability to evade host immune responses during transmission and asymptomatic carriage.

Domestication of Campylobacter jejuni NCTC 11168

Reference and type strains of well-known bacteria have been a cornerstone of microbiology research for decades. The sharing of well-characterized isolates among laboratories has run in parallel with research efforts and enhanced the reproducibility of experiments, leading to a wealth of knowledge about trait variation in different species and the underlying genetics. Campylobacter jejuni strain NCTC 11168, deposited at the National Collection of Type Cultures in 1977, has been adopted widely as a reference strain by researchers worldwide and was the first Campylobacter for which the complete genome was published (in 2000). In this study, we collected 23 C. jejuni NCTC 11168 reference isolates from laboratories across the UK and compared variation in simple laboratory phenotypes with genetic variation in sequenced genomes. Putatively identical isolates, identified previously to have aberrant phenotypes, varied by up to 281 SNPs (in 15 genes) compared to the most recent reference strain. Isolates also display considerable phenotype variation in motility, morphology, growth at 37 °C, invasion of chicken and human cell lines, and susceptibility to ampicillin. This study provides evidence of ongoing evolutionary change among C. jejuni isolates as they are cultured in different laboratories and highlights the need for careful consideration of genetic variation within laboratory reference strains. This article contains data hosted by Microreact.

Determination of Repeat Number and Expression States of Phase-Variable Loci Through Next Generation Sequencing and Bioinformatic Analysis

Phase variation (PV) enables high frequency, reversible switches in expression of genetic loci across numerous species of bacteria. A major mechanism of PV in bacteria is the use of slipped strand mispairing across simple sequence repeats (SSRs). The generation and online availability of genomic datasets enables a comprehensive analysis of the distribution and composition of SSRs across multiple bacterial genomes of a species. PhasomeIt is a program that was developed to rapidly identify SSRs, to determine whether these SSRs mediate PV and to find homologous PV loci across multiple genomes. We describe use of this program for analysis of neisserial genomes. We further describe a method to reassemble specific PV loci to allow analysis of large repeat tracts which are often poorly assembled due to inherent drawbacks of the Illumina next generation sequencing (NGS) platform. These methodologies allow for rapid analysis of a major mechanism of PV across numerous species of Neisseria and other bacterial species.

Mutation and Selection in Bacteria: Modelling and Calibration

Temporal evolution of a clonal bacterial population is modelled taking into account reversible mutation and selection mechanisms. For the mutation model, an efficient algorithm is proposed to verify whether experimental data can be explained by this model. The selection–mutation model has unobservable fitness parameters, and, to estimate them, we use an Approximate Bayesian Computation algorithm. The algorithms are illustrated using in vitro data for phase variable genes of Campylobacter jejuni.