Phasome analysis of pathogenic and commensal Neisseria species expands the known repertoire of phase variable genes, and highlights common adaptive strategies

100+ years of phase variation: the premier bacterial bet-hedging phenomenon

Stochastic, reversible switches in the expression of Salmonella flagella variants were first described by Andrewes in 1922. Termed phase variation (PV), subsequent research found that this phenomenon was widespread among bacterial species and controlled expression of major determinants of bacterial-host interactions. Underlying mechanisms were not discovered until the 1970s/1980s but were found to encompass intrinsic aspects of DNA processes (i.e. DNA slippage and recombination) and DNA modifications (i.e. DNA methylation). Despite this long history, discoveries are ongoing with expansions of the phase-variable repertoire into new organisms and novel insights into the functions of known loci and switching mechanisms. Some of these discoveries are somewhat controversial as the term 'PV' is being applied without addressing key aspects of the phenomenon such as whether mutations or epigenetic changes are reversible and generated prior to selection. Another 'missing' aspect of PV research is the impact of these adaptive switches in real-world situations. This review provides a perspective on the historical timeline of the discovery of PV, the current state-of-the-art, controversial aspects of classifying phase-variable loci and possible 'missing' real-world effects of this phenomenon.

Short tandem repeats delineate gene bodies across eukaryotes

Short tandem repeats (STRs) have emerged as important and hypermutable sites where genetic variation correlates with gene expression in plant and animal systems. Recently, it has been shown that a broad range of transcription factors (TFs) are affected by STRs near or in the DNA target binding site. Despite this, the distribution of STR motif repetitiveness in eukaryote genomes is still largely unknown. Here, we identify monomer and dimer STR motif repetitiveness in 5.1 billion 10-bp windows upstream of translation starts and downstream of translation stops in 25 million genes spanning 1270 species across the eukaryotic Tree of Life. We report that all surveyed genomes have gene-proximal shifts in motif repetitiveness. Within genomes, variation in gene-proximal repetitiveness landscapes correlated to the function of genes; genes with housekeeping functions were depleted in upstream and downstream repetitiveness. Furthermore, the repetitiveness landscapes correlated with TF binding sites, indicating that gene function has evolved in conjunction with cis-regulatory STRs and TFs that recognize repetitive sites. These results suggest that the hypermutability inherent to STRs is canalized along the genome sequence and contributes to regulatory and eco-evolutionary dynamics in all eukaryotes.

High-throughput phenotype-to-genotype testing of meningococcal carriage and disease isolates detects genetic determinants of disease-relevant phenotypic traits

Genome-wide association studies (GWAS) with binary or single phenotype data have successfully identified disease-associated genotypes and determinants of antimicrobial resistance. We describe a novel phenotype-to-genotype approach for a major bacterial pathogen that involves simultaneously testing for associations among multiple disease-related phenotypes and linkages between phenotypic variation and genetic determinants. High-throughput assays quantified variation among 163 Neisseria meningitidis serogroup W ST-11 clonal complex isolates for 11 phenotypic traits. A comparison of carriage and two disease subgroups detected significant differences between groups for eight phenotypic traits. Candidate genotypic testing indicated that indels in csw, a capsular biosynthesis gene, were associated with reduced survival in antibody-depleted heat-inactivated serum. GWAS testing detected 341 significant genetic variants (3 single-nucleotide polymorphisms and 338 unitigs) across all traits except serum bactericidal antibody-depleted assays. Growth traits were associated with variants of capsular biosynthesis genes, carbonic anhydrase, and an iron-uptake system while adhesion-linked variation was in pilC2, marR, and mutS. Multiple phase variation states or combinatorial phasotypes were associated with significant differences in multiple phenotypes. Controlling for group effects through regression and recursive random forest approaches detected group-independent effects for nalP with biofilm formation and fetA with a growth trait. Through random forest testing, nine phenotypes were weakly predictive of MenW:cc11 sub-lineage, original or 2013, for disease isolates while three characteristics separated carriage and disease isolates with >80% accuracy. This study demonstrates the power of combining high-throughput phenotypic testing of pathogenically relevant isolate collections with genomics for identifying genetic determinants of specific disease-relevant phenotypes and the pathobiology of microbial pathogens.IMPORTANCENext-generation sequencing technologies have led to the creation of extensive microbial genome sequence databases for several bacterial pathogens. Mining of these databases is now imperative for unlocking the maximum benefits of these resources. We describe a high-throughput methodology for detecting associations between phenotypic variation in multiple disease-relevant traits and a range of genetic determinants for Neisseria meningitidis, a major causative agent of meningitis and septicemia. Phenotypic variation in 11 disease-related traits was determined for 163 isolates of the hypervirulent ST-11 lineage and linked to specific single-nucleotide polymorphisms, short sequence variants, and phase variation states. Application of machine learning algorithms to our data outputs identified combinatorial phenotypic traits and genetic variants predictive of a disease association. This approach overcomes the limitations of generic meta-data, such as disease versus carriage, and provides an avenue to explore the multi-faceted nature of bacterial disease, carriage, and transmissibility traits.

Microbial Primer: Phase variation - survival and adaptability by generation of a diverse population

Phase variation is defined as the rapid and reversible switching of gene expression, and typically occurs in genes encoding surface features in small genome bacterial pathogens. Phase variation has evolved to provide an extra survival mechanism in bacteria that lack multiple 'sense-and-respond' gene regulation systems. Many bacterial pathogens also encode DNA methyltransferases that are phase-variable, controlling systems called 'phasevarions' (phase-variable regulons). This primer will summarize the current understanding of phase variation, describing the role of major phase-variable factors, and phasevarions, in bacterial pathobiology.

Sweet complexity: O-linked protein glycosylation in pathogenic Neisseria

The genus Neisseria, which colonizes mucosal surfaces, includes both commensal and pathogenic species that are exclusive to humans. The two pathogenic Neisseria species are closely related but cause quite different diseases, meningococcal sepsis and meningitis (Neisseria meningitidis) and sexually transmitted gonorrhea (Neisseria gonorrhoeae). Although obvious differences in bacterial niches and mechanisms for transmission exists, pathogenic Neisseria have high levels of conservation at the levels of nucleotide sequences, gene content and synteny. Species of Neisseria express broad-spectrum O-linked protein glycosylation where the glycoproteins are largely transmembrane proteins or lipoproteins localized on the cell surface or in the periplasm. There are diverse functions among the identified glycoproteins, for example type IV biogenesis proteins, proteins involved in antimicrobial resistance, as well as surface proteins that have been suggested as vaccine candidates. The most abundant glycoprotein, PilE, is the major subunit of pili which are an important colonization factor. The glycans attached can vary extensively due to phase variation of protein glycosylation (pgl) genes and polymorphic pgl gene content. The exact roles of glycosylation in Neisseria remains to be determined, but increasing evidence suggests that glycan variability can be a strategy to evade the human immune system. In addition, pathogenic and commensal Neisseria appear to have significant glycosylation differences. Here, the current knowledge and implications of protein glycosylation genes, glycan diversity, glycoproteins and immunogenicity in pathogenic Neisseria are summarized and discussed.

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.

Adaptive protein evolution through length variation of short tandem repeats in Arabidopsis

Intrinsically disordered protein regions are of high importance for biotic and abiotic stress responses in plants. Tracts of identical amino acids accumulate in these regions and can vary in length over generations because of expansions and retractions of short tandem repeats at the genomic level. However, little attention has been paid to what extent length variation is shaped by natural selection. By environmental association analysis on 2514 length variable tracts in 770 whole-genome sequenced Arabidopsis thaliana, we show that length variation in glutamine and asparagine amino acid homopolymers, as well as in interaction hotspots, correlate with local bioclimatic habitat. We determined experimentally that the promoter activity of a light-stress gene depended on polyglutamine length variants in a disordered transcription factor. Our results show that length variations affect protein function and are likely adaptive. Length variants modulating protein function at a global genomic scale has implications for understanding protein evolution and eco-evolutionary biology.

Variable disruption of epithelial monolayers by Neisseria meningitidis carriage isolates of the hypervirulent MenW cc11 and MenY cc23 lineages

Colonization of mucosal tissues by Neisseria meningitidis requires adhesion mediated by the type IV pilus and multiple outer-membrane proteins. Penetration of the mucosa and invasion of epithelial cells are thought to contribute to host persistence and invasive disease. Using Calu-3 cell monolayers grown at an air-liquid interface, we examined adhesion, invasion and monolayer disruption by carriage isolates of two clonal complexes of N. meningitidis. Carriage isolates of both the serogroup Y cc23 and the hypervirulent serogroup W cc11 lineages exhibited high levels of cellular adhesion, and a variable disruption phenotype across independent isolates. Inactivation of the gene encoding the main pilus sub-unit in multiple cc11 isolates abrogated both adhesive capacity and ability to disrupt epithelial monolayers. Contrastingly, inactivation of the phase-variable opa or nadA genes reduced adhesion and invasion, but not disruption of monolayer integrity. Adherence of tissue-disruptive meningococci correlated with loss of staining for the tight junction protein, occludin. Intriguingly, in a pilus-negative strain background, we observed compensatory ON switching of opa genes, which facilitated continued adhesion. We conclude that disruption of epithelial monolayers occurs in multiple meningococcal lineages but can vary during carriage and is intimately linked to pilus-mediated adhesion.

Insights from targeting transferrin receptors to develop vaccines for pathogens of humans and food production animals

While developing vaccines targeting surface transferrin receptor proteins in Gram-negative pathogens of humans and food production animals, the common features derived from their evolutionary origins has provided us with insights on how improvements could be implemented in the various stages of research and vaccine development. These pathogens are adapted to live exclusively on the mucosal surfaces of the upper respiratory or genitourinary tract of their host and rely on their receptors to acquire iron from transferrin for survival, indicating that there likely are common mechanisms for delivering transferrin to the mucosal surfaces that should be explored. The modern-day receptors are derived from those present in bacteria that lived over 320 million years ago. The pathogens represent the most host adapted members of their bacterial lineages and may possess factors that enable them to have strong association with the mucosal epithelial cells, thus likely reside in a different niche than the commensal members of the bacterial lineage. The bacterial pathogens normally lead a commensal lifestyle which presents challenges for development of relevant infection models as most infection models either exclude the early stages of colonization or subsequent disease development, and the immune mechanisms at the mucosal surface that would prevent disease are not evident. Development of infection models emulating natural horizontal disease transmission are also lacking. Our aim is to share our insights from the study of pathogens of humans and food production animals with individuals involved in vaccine development, maintaining health or regulation of products in the human and animal health sectors.

Genetic, physiological, and cellular heterogeneities of bacterial pathogens in food matrices: Consequences for food safety

In complex food systems, bacteria live in heterogeneous microstructures, and the population displays phenotypic heterogeneities at the single-cell level. This review provides an overview of spatiotemporal drivers of phenotypic heterogeneity of bacterial pathogens in food matrices at three levels. The first level is the genotypic heterogeneity due to the possibility for various strains of a given species to contaminate food, each of them having specific genetic features. Then, physiological heterogeneities are induced within the same strain, due to specific microenvironments and heterogeneous adaptative responses to the food microstructure. The third level of phenotypic heterogeneity is related to cellular heterogeneity of the same strain in a specific microenvironment. Finally, we consider how these phenotypic heterogeneities at the single-cell level could be implemented in mathematical models to predict bacterial behavior and help ensure microbiological food safety.

The Host-Pathogen Interactions and Epicellular Lifestyle of Neisseria meningitidis

Neisseria meningitidis is a gram-negative diplococcus and a transient commensal of the human nasopharynx. It shares and competes for this niche with a number of other Neisseria species including N. lactamica, N. cinerea and N. mucosa. Unlike these other members of the genus, N. meningitidis may become invasive, crossing the epithelium of the nasopharynx and entering the bloodstream, where it rapidly proliferates causing a syndrome known as Invasive Meningococcal Disease (IMD). IMD progresses rapidly to cause septic shock and meningitis and is often fatal despite aggressive antibiotic therapy. While many of the ways in which meningococci survive in the host environment have been well studied, recent insights into the interactions between N. meningitidis and the epithelial, serum, and endothelial environments have expanded our understanding of how IMD develops. This review seeks to incorporate recent work into the established model of pathogenesis. In particular, we focus on the competition that N. meningitidis faces in the nasopharynx from other Neisseria species, and how the genetic diversity of the meningococcus contributes to the wide range of inflammatory and pathogenic potentials observed among different lineages.

Modelling evolutionary pathways for commensalism and hypervirulence in Neisseria meningitidis

Neisseria meningitidis, the meningococcus, resides exclusively in humans and causes invasive meningococcal disease (IMD). The population of N. meningitidis is structured into stable clonal complexes by limited horizontal recombination in this naturally transformable species. N. meningitidis is an opportunistic pathogen, with some clonal complexes, such as cc53, effectively acting as commensal colonizers, while other genetic lineages, such as cc11, are rarely colonizers but are over-represented in IMD and are termed hypervirulent. This study examined theoretical evolutionary pathways for pathogenic and commensal lineages by examining the prevalence of horizontally acquired genomic islands (GIs) and loss-of-function (LOF) mutations. Using a collection of 4850 genomes from the BIGSdb database, we identified 82 GIs in the pan-genome of 11 lineages (10 hypervirulent and one commensal lineage). A new computational tool, Phaser, was used to identify frameshift mutations, which were examined for statistically significant association with genetic lineage. Phaser identified a total of 144 frameshift loci of which 105 were shown to have a statistically significant non-random distribution in phase status. The 82 GIs, but not the LOF loci, were associated with genetic lineage and invasiveness using the disease carriage ratio metric. These observations have been integrated into a new model that infers the early events of the evolution of the human adapted meningococcus. These pathways are enriched for GIs that are involved in modulating attachment to the host, growth rate, iron uptake and toxin expression which are proposed to increase competition within the meningococcal population for the limited environmental niche of the human nasopharynx. We surmise that competition for the host mucosal surface with the nasopharyngeal microbiome has led to the selection of isolates with traits that enable access to cell types (non-phagocytic and phagocytic) in the submucosal tissues leading to an increased risk for IMD.

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.

Allelic polymorphisms in a glycosyltransferase gene shape glycan repertoire in the O-linked protein glycosylation system of Neisseria

Glycosylation of multiple proteins via O-linkage is well documented in bacterial species of Neisseria of import to human disease. Recent studies of protein glycosylation (pgl) gene distribution established that related protein glycosylation systems occur throughout the genus including nonpathogenic species. However, there are inconsistencies between pgl gene status and observed glycan structures. One of these relates to the widespread distribution of pglG, encoding a glycosyltransferase that in Neisseria elongata subsp. glycolytica is responsible for the addition of di-N-acetyl glucuronic acid at the third position of a tetrasaccharide. Despite pglG residing in strains of N. gonorrhoeae, N. meningitidis and N. lactamica, no glycan structures have been correlated with its presence in these backgrounds. Moreover, PglG function in N. elongata subsp. glycolytica minimally requires UDP-glucuronic acid (GlcNAcA), and yet N. gonorrhoeae, N. meningitidis and N. lactamica lack pglJ, the gene whose product is essential for UDP-GlcNAcA synthesis. We examined the functionality of pglG alleles from species spanning the Neisseria genus by genetic complementation in N. elongata subsp. glycolytica. The results indicate that select pglG alleles from N. meningitidis and N. lactamica are associated with incorporation of an N-acetyl-hexosamine at the third position and reveal the potential for an expanded glycan repertoire in those species. Similar experiments using pglG from N. gonorrhoeae failed to find any evidence of function suggesting that those alleles are missense pseudogenes. Taken together, the results are emblematic of how allelic polymorphisms can shape bacterial glycosyltransferase function and demonstrate that such alterations may be constrained to distinct phylogenetic lineages.

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.

Localized Hypermutation is the Major Driver of Meningococcal Genetic Variability during Persistent Asymptomatic Carriage

Host persistence of bacteria is facilitated by mutational and recombinatorial processes that counteract loss of genetic variation during transmission and selection from evolving host responses. Genetic variation was investigated during persistent asymptomatic carriage of Neisseria meningitidis Interrogation of whole-genome sequences for paired isolates from 25 carriers showed that de novo mutations were infrequent, while horizontal gene transfer occurred in 16% of carriers. Examination of multiple isolates per time point enabled separation of sporadic and transient allelic variation from directional variation. A comprehensive comparative analysis of directional allelic variation with hypermutation of simple sequence repeats and hyperrecombination of class 1 type IV pilus genes detected an average of seven events per carrier and 2:1 bias for changes due to localized hypermutation. Directional genetic variation was focused on the outer membrane with 69% of events occurring in genes encoding enzymatic modifiers of surface structures or outer membrane proteins. Multiple carriers exhibited directional and opposed switching of allelic variants of the surface-located Opa proteins that enables continuous expression of these adhesins alongside antigenic variation. A trend for switching from PilC1 to PilC2 expression was detected, indicating selection for specific alterations in the activities of the type IV pilus, whereas phase variation of restriction modification (RM) systems, as well as associated phasevarions, was infrequent. We conclude that asymptomatic meningococcal carriage on mucosal surfaces is facilitated by frequent localized hypermutation and horizontal gene transfer affecting genes encoding surface modifiers such that optimization of adhesive functions occurs alongside escape of immune responses by antigenic variation.IMPORTANCE Many bacterial pathogens coexist with host organisms, rarely causing disease while adapting to host responses. Neisseria meningitidis, a major cause of meningitis and septicemia, is a frequent persistent colonizer of asymptomatic teenagers/young adults. To assess how genetic variation contributes to host persistence, whole-genome sequencing and hypermutable sequence analyses were performed on multiple isolates obtained from students naturally colonized with meningococci. High frequencies of gene transfer were observed, occurring in 16% of carriers and affecting 51% of all nonhypermutable variable genes. Comparative analyses showed that hypermutable sequences were the major mechanism of variation, causing 2-fold more changes in gene function than other mechanisms. Genetic variation was focused on genes affecting the outer membrane, with directional changes in proteins responsible for bacterial adhesion to host surfaces. This comprehensive examination of genetic plasticity in individual hosts provides a significant new platform for rationale design of approaches to prevent the spread of this pathogen.

Neisseria meningitidis: using genomics to understand diversity, evolution and pathogenesis

Meningococcal disease remains an important cause of morbidity and death worldwide despite the development and increasing implementation of effective vaccines. Elimination of the disease is hampered by the enormous diversity and antigenic variability of the causative agent, Neisseria meningitidis, one of the most variable bacteria in nature. These features are attained mainly through high rates of horizontal gene transfer and alteration of protein expression through phase variation. The recent availability of whole-genome sequencing (WGS) of large-scale collections of N. meningitidis isolates from various origins, databases to facilitate storage and sharing of WGS data and the concomitant development of effective bioinformatics tools have led to a much more thorough understanding of the diversity of the species, its evolution and population structure and how virulent traits may emerge. Implementation of WGS is already contributing to enhanced epidemiological surveillance and is essential to ascertain the impact of vaccination strategies. This Review summarizes the recent advances provided by WGS studies in our understanding of the biology of N. meningitidis and the epidemiology of meningococcal disease.

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.

Phase-variable bacterial loci: how bacteria gamble to maximise fitness in changing environments

Phase-variation of genes is defined as the rapid and reversible switching of expression - either ON-OFF switching or the expression of multiple allelic variants. Switching of expression can be achieved by a number of different mechanisms. Phase-variable genes typically encode bacterial surface structures, such as adhesins, pili, and lipooligosaccharide, and provide an extra contingency strategy in small-genome pathogens that may lack the plethora of 'sense-and-respond' gene regulation systems found in other organisms. Many bacterial pathogens also encode phase-variable DNA methyltransferases that control the expression of multiple genes in systems called phasevarions (phase-variable regulons). The presence of phase-variable genes allows a population of bacteria to generate a number of phenotypic variants, some of which may be better suited to either colonising certain host niches, surviving a particular environmental condition and/or evading an immune response. The presence of phase-variable genes complicates the determination of an organism's stably expressed antigenic repertoire; many phase-variable genes are highly immunogenic, and so would be ideal vaccine candidates, but unstable expression due to phase-variation may allow vaccine escape. This review will summarise our current understanding of phase-variable genes that switch expression by a variety of mechanisms, and describe their role in disease and pathobiology.

Med15: Glutamine-Rich Mediator Subunit with Potential for Plasticity

The Mediator complex is required for basal activity of the RNA polymerase (Pol) II transcriptional apparatus and for responsiveness to some activator proteins. Med15, situated in the Mediator tail, plays a role in transmitting regulatory information from distant DNA-bound transcription factors to the transcriptional apparatus poised at promoters. Yeast Med15 and its orthologs share an unusual, glutamine-rich amino acid composition. Here, we discuss this sequence feature and the tendency of polyglutamine tracts to vary in length among strains of Saccharomyces cerevisiae, and we propose that different polyglutamine tract lengths may be adaptive within certain domestication habitats.

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.

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.