BACKGROUND
neisseria meningitidis  is a gram-negative, encapsulated bacterium that is present, as a commensal organism, in the nasopharyngeal cavity of five to ten percent of the adult population  <cit> .
despite its prevalence as a harmless organism some strains, for reasons not yet completely understood, can cross the epithelial barrier and enter the bloodstream, causing septicemia and life-threatening disease .
in order to maximize its fitness in the diverse and changing environments offered by the host-pathogen interplay, nm has evolved multiple and complementary adaptive strategies  <cit> .
one such mechanism is represented by simple sequence repeats , contiguous iterations of short dna motifs, generally assumed to range from  <dig> to  <dig> nucleotides in length  <cit> .
ssrs, often called contingency loci, are hyper-mutable dna sequences mediating high frequency, stochastic, reversible and heritable genotypic switching  <cit> , whose number of tandemly repeated motifs can vary with relatively high frequencies  <cit> .
a variable number ssr  located within a coding sequence can change translation by introducing frameshifts in the reading frame  or, when located in the proximity of a promoter, can modulate transcription either by switching between alternative translational start sites  <cit>  or by altering the promoter sequence .
the first nm genome-wide ssrs analyses were based on the single genome sequences determined for the mc <dig>  <cit>  and z <dig>  <cit>  strains.
putative phase-variable genes were predicted on the basis of cis factors such as ssr sequence, number of repetitions and sequence context.
for a subset of genes phase variation was confirmed experimentally.
comparative investigations, in which the genomic sequences of strains n. gonorrhoeae fa <dig> and nm fam <dig> were added to the analysis, suggested that the presence of length polymorphisms among orthologous ssrs present in different genomes is a reliable predictor of phase variation  <cit> .
overall,  <dig> putative phase variable genes were reported in the neisseria genus, of which  <dig> were specific to the nm species and included genes involved in cell adhesion , capsule formation , biosynthesis of the lipopolysaccharide layer , cell-surface receptors  and restriction/modification systems  <cit> .
among these, phase variation was confirmed experimentally for only  <dig> genes, primarily combining pcr with immunoblotting techniques  or with northern blotting and quantitative pcr .
here, we report that the limited number of sequences employed in previous vnssrs studies, together with the lack of high-throughput technologies  suitable for a large-scale validation of putative phase variable genes, have potentially led to an underestimation of the overall impact of ssr-mediated phase variation on nm phenotype and fitness.
also, the evolutionary dynamics underlying the generation of vnssrs in nm and other species has not been fully elucidated.
recently, new insights into the nm population structure and dynamics have been generated through a comparative genomic analysis based on the complete genome sequences of 20 nm strains, including multiple isolates from each of the most virulent clonal complexes  <cit> .
also, the advent of ngs technologies have allowed for an unprecedented sequencing depth across the whole genome, such that the presence of mixed populations at specific loci of a single genome can be detected with high throughput, including those variants occurring at low frequencies during bacterial replication .
in the present study we took advantage of nm pan-genome data to map, characterize and infer functional impact and evolutionary properties of ssr contingency loci, and ngs to experimentally validate in silico predictions on stochastic genotypic switching.
RESULTS
two hundred seventy-seven meningococcal genes are associated with simple sequence repeats that show inter-strain length polymorphisms
an average of  <dig> ssrs were identified in each genome, the majority of which  were represented by homopolymeric tracts .
the number of ssrs was similar across genomes  and their abundance was found to be significantly higher than random expectation .
significant deviation from neutrality was also observed for each individual ssr type.
we identified  <dig> clusters of orthologous ssrs, 35 %  of which were represented in every genome analyzed.
among these,  <dig> clusters had polymorphisms in the repeat tract length across different genomes, located either in intragenic regions  or in the associated gene’s upstream region .
these were selected for further analyses based on their potential to change gene expression.
consistent with previous findings  <cit> , homopolymeric tracts were the most abundant simple sequence repeat type, with  <dig> vnssrs clusters divided into  <dig> a/t and 59 g/c tracts.
the remaining  <dig> clusters were represented by tandem repeats of two- to nine-nucleotides motifs .
positional analysis identified  <dig>  intragenic and  <dig>  intergenic vnssrs, respectively.
interestingly, all vnssrs with a unit motif of three-, six- and nine-nucleotides, for which a length polymorphism would not result in the disruption of the reading frame, were found to be located within coding sequences .
finally, thirty-two genes were associated with more than a single vnssr, resulting in a total of  <dig> phase variable gene candidates.
ssr contingency loci show different genotype switching frequencies
based on repeat variability and intra-strain phylogenetic relationships we stratified the  <dig> vnssr-associated genes into  <dig> strong,  <dig> moderate and  <dig> weak candidates for phase variation.
most weak candidate phase variable genes  were associated with a/t homopolymeric tracts while g/c repeats were the most abundant vnssr type among the strong candidate genes .
vnssrs other than homopolymeric repeats were found in  <dig> moderate  and  <dig> strong  candidate genes, respectively .fig.
<dig> properties of identified vnssrs.
a relative contribution of different vnssrs types to the strong, moderate and weak candidate phase variable genes.
other: all non-homopolymeric repeats.
g/c: guanine or cytosine variable homopolymeric tracts.
a/t: adenine or thymine variable homopolymeric tracts.
b nei’s diversity index distribution of vnssrs associated with the  <dig> strong,  <dig> moderate and  <dig> weak candidate phase variable genes.
genes associated with intragenic and intergenic vnssrs are represented full and empty circles, respectively.
the thick horizontal lines represent distributions median.
p: wilcoxon-mann-whitney p value
intrinsic variability of vnssrs clusters, defined as the nei’s diversity index  <cit> , was then computed for the three gene sets.
vnssrs associated to strong candidate genes were significantly more variable than those associated to moderate genes, which in turn showed higher variability than those linked to weak phase variable gene candidates , providing supporting evidence that genes associated with a different likelihood of undergoing phase variation exist and that these can be identified through the analysis of a suitable genome dataset.
a positive correlation between the ssr tract length and repeat instability has been proposed  <cit> , possibly as a result of a decreased proof-reading efficiency of the dna polymerase over the longer tracts  <cit> .
we could confirm such association for the strong  but not for the moderate or weak candidate gene pools .
vnssrs variation is detectable in single colonies grown overnight under non-selective conditions
the genomes of five of the analyzed strains were re-sequenced at high depth using a next generation sequencing platform with the goal of detecting repeats length polymorphisms occurring in single colonies during an overnight growth.
analysis of the sequencing reads confirmed repeat length polymorphisms in  <dig> of the  <dig> predicted phase variable genes, distributed in  <dig> strong ,  <dig> moderate  and  <dig> weak  putative phase variable genes .
the same procedure applied to  <dig> randomly chosen, non-variable ssrs detected length polymorphisms in  <dig> loci , a significantly lower proportion compared to vnssrs .
additionally, the average frequencies of polymorphic reads  observed for vnssrs was  <dig>  %, a higher proportion compared to the  <dig> false positive ssrs  or to the allowed read mapping error rate .
frequency distributions of polymorphic reads identified in the  <dig> genomes are reported in additional file 2: figure s <dig> fig.
<dig> phase variable genes confirmed by illumina sequencing.
numbers inside bars indicate the frequency of phase variable genes present in each category.
confirmed: predicted phase variable genes confirmed by deep sequencing data.
not confirmed: predicted phase variable genes not showing evidence of variation.
control: analysis of  <dig> randomly selected ssrs that were predicted as non-variable by the genomic comparative analysis
seventy of the  <dig> validated vnssrs were present in all five re-sequenced genomes.
among these,  <dig>  were shown to be variable in all the  <dig> genomes analyzed, while  <dig>  were confirmed in four or fewer genomes .
interestingly, for any one genome there was a strong association between ssr size  and the likelihood of finding it variable.
specifically, homopolymeric tracts of  <dig> or more nucleotides were predominantly found to be variable in four or five genomes .
because there were only three non-homopolymeric tracts, it was not possible to generalize this conclusion to longer repeats.
additionally, for the  <dig> experimentally confirmed core vnssrs, the proportion of polymorphic reads was found to correlate  with the nei’s diversity index derived from the comparison of the  <dig> available genome sequences.
specifically, vnssrs associated to a nei’s index > <dig>  were mainly characterized by above-average polymorphic reads proportions , meaning that the extent of vnssrs variability can be predicted with good confidence by the comparison of multiple genomic sequences.
finally, the number of phase variable genes that would be confirmed by deep sequencing of a wider genome collection was estimated by applying an approach originally proposed for estimating the bacterial pan-genome size  <cit> .
the regression analysis based on the number of validated phase variable genes predicted that if all the  <dig> genomes were re-sequenced, the number of validated phase-variable genes would have been 146 ±  <dig> , while the analysis of  <dig> genomes could raise the number to 190 ±  <dig> .
the pool of genes whose expression is influenced by vnssrs is wider than previously hypothesized
the pool of putative phase variable genes predicted in this study was compared with the  <dig> genes previously described or predicted to be phase-variable in nm  <cit> .
our approach confirmed  <dig> of the  <dig> genes whose phase variation had been previously demonstrated in nm, as well as  <dig> of the  <dig> putative phase variable genes previously proposed as strong candidates and  <dig> of the  <dig> genes previously reported as either moderate or weak candidates .
the  <dig> genes not confirmed by our method, either didn’t satisfy the criteria for ssr search or didn’t show inter-strain length polymorphisms.
in the particular case of a poly-c tract located within the nmb <dig> coding region, the repeat showed between-strain variation, alternating between the c <dig> and c <dig> states.
however, since the c <dig> was below the minimum length cut off  our approach failed to identify this variable ssr.
remarkably, the wider genome collection allowed for the identification of  <dig> new putative phase variable genes, predicted as  <dig> strong,  <dig> moderate and  <dig> weak candidates respectively .
candidate phase variable genes related to cell surface structure or involved in dna metabolism are overrepresented
the  <dig> strong and moderate candidate genes, described in results section  <dig>  distributed across  <dig> different functional roles .
consistent with previous studies  <cit>  vnssr-associated genes were found to be primarily involved in host-microbial interplay, with “cell envelope”, “transport and binding” and “dna metabolism” functional roles being significantly overrepresented.
intragenic and intergenic vnssrs were both present in the enriched categories.
other biological processes included functions related to the peptide synthesis machinery .
additionally,  <dig> putative phase variable genes were found in functional categories that had not been previously associated to ssr-mediated regulation, such as genes encoding for transposons and prophage related functions, transcription factors and proteins participating in nucleic acids biosynthesis .fig.
<dig> tigr functional roles represented by the  <dig> strong and moderate candidate phase variable genes.
x-axis represents the proportion of putative phase variable genes present in each annotation.
numbers in brackets represent the number of putative phase variable genes over the total number of genes that are associated to a specific function.
genes associated with intragenic and intergenic vnssrs are represented in dark and light grey, respectively.
‘*’: enriched functional roles .
‘+’: functional roles not previously described to be associated with ssr-mediated regulation in nm
the  <dig> weak candidate genes were associated with  <dig> different functional roles, fourteen of which were common to the strong and moderate groups.
this subset showed a different profile, whereby genes were significantly overrepresented in the “hypothetical proteins” and “mobile and extrachromosomal element functions” functional annotations .
most vnssrs are predicted to induce phenotypic changes by modulating the level of gene expression rather than mediating on/off translational switching
we sought to gather new insights into the vnssr-mediated regulatory mechanisms by analyzing the repeat sequence context in the  <dig> strong candidate phase variable genes.
nineteen of the  <dig> vnssrs associated with strong candidate genes were surrounded by multiple in-frame 5′-atg-3′ translational initiation codons .
these vnssrs can shift the origin of translation between alternative start sites, resulting in the modulation of gene expression  <cit> .
such modulation can arise either from the switching between start sites associated with different expression levels or from the production of different protein isoforms, as may be the case for the nmb <dig>  and nmb <dig> open reading frames  , respectively.
different peptide isoforms may also result from vnssrs located in the 3′ portion of the gene  <cit> .
in this location, the vnssr can bring the c-terminal peptide portion out of frame while still allowing the translation of a functional protein.
two orfs, nmb <dig>  and nmb <dig>  were consistent with this mechanism.
in contrast, vnssrs present within the nmb <dig>  nmb <dig> and nmb <dig> reading frames were either a 3- or a 9-nucleotides repeat.
these vnssrs do not interfere with the reading frame but can still influence the cell phenotype by producing proteins with different structures and/or post-translational profiles .
also,  <dig> of the  <dig> analyzed vnssrs were predicted to interfere with the associated gene promoters, as represented by the t-homopolymeric tract located upstream of the nmb <dig> start .
the remaining  <dig> repeats were predicted to introduce frameshifts likely resulting in on-off type of transcriptional regulation, as proposed for nmb <dig> .fig.
<dig> schematic representation of five vnssrs and their sequence context.
a vnssr causing translational start site switching.
b vnssr causing the loss of a membrane-spanning domain.
c vnssr leading to the loss of the peptide c-terminal region.
d vnssr introducing changes in the peptide sequence.
e) vnssr influencing the gene promoter.
f vnssr introducing an inactivating frameshift.
dark grey arrows represent open reading frames.
black arrows marked with atg represent in-frame atg translational start sites.
light grey boxes represent the annotated functional domains.
stripped boxes represent vnssrs and the related tags indicate the repeat unit motif along with the minimum and maximum number of repetitions observed in the  <dig> analyzed genomes.
numbers below each gene indicate the position relative to the annotated translational tart site
most strong candidate phase variable genes are not specific to n. meningitidis
among the  <dig> strong candidate phase variable genes,  <dig> were identified only within nm through sequence database comparison.
the gc content of these genes largely deviates from the typical meningococcal value , suggesting acquisition through lateral gene transfer from outside the genus.
conversely, orthologous copies of the remaining  <dig> genes were also identified in other neisseria species .
among these,  <dig> were associated with the same ssrs observed in nm, indicating that these ssr-contingency loci arose before meningococcus speciation.
in the remaining four cases  the ssr associated with the gene was detected in nm only, suggesting that these contingency loci have evolved following nm speciation .
orthologous copies of  <dig> strong candidate phase variable genes were also identified in species outside of the neisseria genus , including cardiobacterium hominis, haemophilus aegyptius, haemophilus haemolyticus, haemophilus influenzae, kigella dentrificans and streptococcus pneumoniae.
among these, the nmb <dig> encoded gene identified in the three haemophilus species and the  <dig> opacity proteins identified in s. pneumoniae, were associated with the same ssr as in nm , raising the hypothesis that ssr contingency loci are transferred across species through lateral dna exchange.
vnssrs reorganization following n. meningitidis speciation primarily involved cell surface determinants
among the  <dig> vnssrs associated to the  <dig> strong candidate genes identified in multiple taxa,  <dig>  were found to be significantly longer  in nm than in other species, while  <dig> were found to be significantly shorter .
twelve and  <dig> of the repeats with increased length in nm were g/c and a/t homopolymeric tracts, respectively, along with  <dig> other non-homopolymeric ssrs.
these repeats were primarily associated with cell surface determinants, including surface receptors, surface transporters, pili and lipopolysaccharide biosynthetic enzymes.
the  <dig> vnssrs with a reduced size in nm included a poly-t and a poly-c homopolymeric tracts, respectively associated with a hypothetical protein  and to a ferric enterobactin receptor , and a 5′-cttct-3′ repeat associated to an opacity protein .
also in this case, an association  with surface exposed proteins was observed.
overall, these results suggest that evasion of host immune responses, along with changes in surface adhesion properties and tissue tropism, have likely been the main drivers of ssr contingency loci evolution within nm.
vnssrs chromosomic loci are fixed in the meningococcal population
chromosomic loci containing the  <dig> identified vnssrs were tested for signatures of selective pressures that could confirm a beneficial role  of those loci and justify their over-abundance in nm species.
only  <dig>  vnssr loci, all associated with a/t homopolymeric tracts, showed a significant deviation from neutrality with all the applied tests .
signatures of positive selection  were identified in the upstream region of the two strong candidate phase variable genes mtrf , a surface-exposed efflux-pump component involved in hydrophobic antimicrobial resistance  <cit> , and the transcription factor dksa .
a signature of negative selection  was detected in the upstream region of the moderate candidate phase variable gene sstt , which also encodes for a surface transporter.
all other signals were identified in vnssrs associated to weak candidate phase variable genes.
the remaining  <dig>  tested loci did not show signatures of selection, supporting the hypothesis that ssr contingency loci have reached an equilibrium state and are now fixed in the population.
discussion
after the first nm complete genome sequence became available  <cit> , the highly abundant ssrs were largely associated with genes involved in host adaptation, commensalism and virulence.
these initial studies on nm ssrs contingency loci were based on analysis of individual genomes, strains mc <dig>  <cit>  and z <dig>  <cit> .
with the technology available at that time, the generation of new genomic sequences and the validation of the newly proposed phase-variable genes were laborious and time-consuming tasks.
despite the challenges, a few years after the first meningococcal genome sequence was revealed, a study reported a list of 69 nm phase variable genes,  <dig> of which were experimentally confirmed  <cit> .
benefiting from the most recent sequencing technologies, the present study employed a 20 nm genomes dataset, with a resolved phylogenetic structure  <cit> , in a comparative genomic analysis aimed at the study of ssr-driven phase variation as a population-based, as opposed to single cell, adaptive strategy.
all analyzed strains were found to have more ssrs than expected in their chromosomes.
a comparable number of ssrs were also found in two non-pathogenic strains, α <dig> and ox99- <dig>  questioning the hypothesis that ssrs-mediated phase variation is an adaptation associated with the evolution of virulence.
most of the identified ssrs were homopolymeric tracts, with an excess of a/t tandem repeats, suggesting slipped-strand mispairing as the primary mechanism driving the occurrence of ssrs.
it has been proposed that slipped-strand mispairing is preferentially associated with homopolymeric tracts due to their high sequence redundancy  <cit>  and favouring a/t repeats due to the lower stability of the dna-dna hybrid during bacterial genome replication  <cit> .
the likelihood that repeat-associated genes undergo phase variation is variable.
efforts to identify molecular signatures indicative of such likelihood have defined a number of cis factors, including repeat sequence, size  and sequence context, along with the associated gene function  <cit> .
particularly, genes encoding for surface antigens were regarded as strong candidates for phase variation.
we applied a more unsupervised approach, in which only the repeat sequence context and the distribution of its variation over the sampled phylogeny were considered.
the over-representation of a/t repeats and their alleged higher propensity to cause transcriptional slippage compared to g/c repeats  <cit>  would suggest these to be the main contributor of ssr mediated phase variation.
conversely, in agreement with previous findings  <cit> , homopolymeric g/c repeats were the most frequent  among the strong candidate genes pool, while weak candidate genes were almost exclusively associated with a/t repeats .
a primary involvement of c/g homopolymeric in the mediation of phase variation has also been described in c. jejuni  <cit> .
this particular mutational pattern, together with the fact that we didn’t identify traces of selective pressure for most vnssrs, suggest that there may be molecular drivers governing the mutability of these contingency loci.
in support of this hypothesis, increased rates of phase variation have been associated with the loss of the dam dna methylase  <cit>  and the two mismatch repair proteins muts and mutl  <cit> .
regarding non-homopolymeric repeats, 3-, 4-, 5-, 6- and 9-nucelotides vnssrs were also found to be abundant.
three-, 6- and 9-nucleotides vnssrs do not disrupt the reading frame and their abundance may be justified by a reduced selective pressure for their stability.
abundance of 5-nucleotide vnssrs is due to the fact that  <dig> out of  <dig> cases are represented by the same 5′-cttct-3′ repeat associated to different paralogs of an opacity protein .
interestingly, a homolog of this opacity protein, coupled with the same ssr, was also identified in streptococcus pneumoniae genome , providing evidence that ssr-contingency loci can be transferred across different species, or strains of the same species and even different chromosomic loci of the same strain through lateral gene transfer.
conversely, abundance of 4-nucleotides vnssrs, which has been reported both in nm and in haemophilus influenzae  <cit> , has no obvious interpretation.
given the conservation of those repeat sequences and their proximity to the translational origin, the interaction with yet unidentified regulatory elements may be involved.
one such regulatory mechanism has been described for the taaa repeat associated with the nada  promoter  <cit> .
the availability of high-throughput sequencing data enabled a high throughput experimental validation of in silico predicted ssr contingency loci.
the rationale behind this approach was that dna libraries used for illumina sequencing were generated from a bacterial culture derived from a single colony.
after an 18 h growth and approximately  <dig> bacterial replication cycles , a library was expected to contain, and represent in the form of variable alleles, most of the vnssrs variation.
this test provided circumstantial evidence for phase variability of  <dig>   <dig> and 27 % of the strong, moderate and weak candidate genes respectively, leading to a total of  <dig> experimentally confirmed phase variable genes.
this approach also allowed estimating vnssrs variation frequencies.
the  <dig> validated vnssr showed an average proportion of alternative alleles of  <dig>  %, indicating an average mutation rate of once every seven cell replication cycles.
however, given the contribution of environmental, population and molecular factors to the mutability of ssrs in phase variable loci, it is likely that real frequencies may differ from those we observed in vitro.
ssr contingency loci were over-represented among genes related to cell surface functions and structure or involved in dna metabolism, primarily dna restriction and modification.
these results, which are in agreement with previous findings  <cit> , reflect a primary involvement of ssr-mediated phase variation in microbe-host interaction.
dis-regulation of those genes, in fact, alters nm ability to adhere to and invade host cells, its ability to scavenge required nutrients and to evade host immune defenses by generating antigenic variants  <cit> .
additionally, putative phase variable genes were found in several other functional classes, indicating that ssr-mediated gene regulation interferes with a set of metabolic functions wider than previously hypothesized.
our attempt to understand how vnssrs have evolved in nm revealed that most strong candidate phase variable genes  are widespread to the neisseria genus and that about half of them have undergone size change after meningococcal speciation.
this is suggestive of the fact that most ssr-associated contingency loci have developed prior to nm speciation and have subsequently evolved, possibly as a consequence of the process of adaptation to the human host environment.
in agreement with this is the fact that identified phase variable genes were enriched for cell surface determinants, suggesting that evasion of host immune responses and changes in surface adhesion properties and tissue tropism are the main drivers of ssr contingency loci evolution in nm.
regarding the remaining  <dig> genes that are specific to nm, they were characterized by a g/c content that deviates from that of meningococcus , suggesting a possible acquisition through horizontal gene transfer.
phase variable genes candidates were also identified in species outside the neisseria genus .
interestingly, all of them are human pathogens that coexist in the respiratory flora of healthy individuals.
moreover, we found that genes shared by nm, s. pneumoniae and three haemophilus members were associated with the same repeat element as in nm, supporting the possibility that ssr contingency loci are exchanged across species through lateral dna exchange.
overall, our findings suggest that ssr contingency loci have been acquired by nm either by spontaneous evolution, inheritance from an ancestral species or acquisition by means of lateral dna exchange.
in the context of a commensal organism, constantly interacting with its host, the onset of new mutations or allele variants may lead to an increased fitness , a reduced fitness  or have no effect on fitness   <cit> .
one of the most accredited hypotheses is that vnssrs have evolved as an adaptive strategy, in organisms like nm, to increase their fitness in the hostile and mutable host environment  <cit> .
we tested all vnssr-containing loci for signatures of selection that could confirm a beneficial role of those loci and justify their over-abundance.
surprisingly, the analysis revealed that most ssr-contingency loci are not under selective pressure, suggesting that nm has had enough time to adapt to its niche and has apparently reached an equilibrium state.
it must be noted, however, that substantial homologous recombination and repeated population expansions and bottlenecks, along with possible sampling biases, are likely to confound or even obscure neutrality tests results.
nonetheless, we envision that, applying this approach to a wider, unbiased strain collection, coming from a single outbreak, will have the potential of providing valuable insights into the understanding of vnssrs evolutionary processes.
despite a proportion of coding dna in nm higher than 80 %  <cit> , we observed that 51 % of the identified vnssrs are found within intergenic regions.
this bias toward intergenic vnssrs, which may be explained by a reduced selective pressure along non-coding regions, can also be indicative of a regulation of gene expression occurring at the transcriptional level.
this type of vnssrs can indeed modulate promoter efficiency by changing the relative spacing between key promoter elements  <cit> .
differently, intragenic vnssrs, which can induce frameshift mutations, are responsible for an on/of type of regulation occurring at the translational level  <cit> .
additionally, given the existence of alternative translational start codons, intragenic vnssrs located in the 5′ portion of a gene can switch between multiple active reading frames, with the potential of affecting the level of gene expression  <cit>  or producing alternative peptide isoforms  <cit> .
in the present analysis we observed that most intragenic vnssrs are located in the proximity of the translational start site and that among the  <dig> strong candidate phase variable genes, there are  <dig> cases in which the vnssr was found between two or more in-frame atg translational start sites.
twenty other vnssrs were predicted to either interact with the gene promoter or to produce alternative peptides isoforms.
overall, collected evidences suggest that ssr-contingency loci have evolved in such a way to allow the maximum degree of flexibility, reflected by their ability to finely modulate gene expression level, or the functionality of the resulting protein, rather than merely switching between the two extreme cases of a gene being expressed or silenced through the introduction of frameshifts.
CONCLUSIONS
in conclusion, this analysis, which capitalized on the most recent dna sequencing technologies, confirmed the fundamental contribution of ssr contingency loci in promoting phenotypic variation and their deep implications in nm survival and adaptation to the surrounding environment.
our unsupervised approach allowed the generation of a panel of  <dig> genes whose expression may be controlled or influenced by ssr elements and to provide corroborative evidence for the phase variability for  <dig> of them.
functional characterization of these genes highlighted an enrichment for cell surface determinants, which included members of the adhesins, evasins, lipopolisaccharide biosynthesis and nutrient scavenging protein families.
such proteins have been attributed multiple roles in attachment of bacterial cells to host membranes  <cit> , in regulating bacterial resistance to both innate and adaptive immune system  and as major determinants of meningococcal invasive disease  <cit> .
ssr-mediated phase variation has also been reported to control the expression of multiple immunogens currently included in commercial meningococcal vaccines  <cit>  and to modulate resistance to antimicrobial agents .
we therefore envision that future studies on meningococcal microbe-host interaction and studies aimed at the identification of new vaccine antigens and antimicrobial molecules will benefit from this comprehensive characterization of the meningococcal putative phase variable genes repertoire.

 neisseria meningitidis is a gram-negative, encapsulated bacterium that is present, as a commensal organism, in the nasopharyngeal cavity of five to ten percent of the adult population. despite its prevalence as a harmless organism some strains, for reasons not yet completely understood, can cross the epithelial barrier and enter the bloodstream, causing septicemia and life-threatening disease. in order to maximize its fitness in the diverse and changing environments offered by the host-pathogen interplay, nm has evolved multiple and complementary adaptive strategies. one such mechanism is represented by simple sequence repeats, contiguous iterations of short dna motifs, generally assumed to range from to nucleotides in length. ssrs, often called contingency loci, are hyper-mutable dna sequences mediating high frequency, stochastic, reversible and heritable genotypic switching, whose number of tandemly repeated motifs can vary with relatively high frequencies. a variable number ssr located within a coding sequence can change translation by introducing frameshifts in the reading frame or, when located in the proximity of a promoter, can modulate transcription either by switching between alternative translational start sites or by altering the promoter sequence. the first nm genome-wide ssrs analyses were based on the single genome sequences determined for the mc and z strains. putative phase-variable genes were predicted on the basis of cis factors such as ssr sequence, number of repetitions and sequence context. for a subset of genes phase variation was confirmed experimentally. comparative investigations, in which the genomic sequences of strains n. gonorrhoeae fa and nm fam were added to the analysis, suggested that the presence of length polymorphisms among orthologous ssrs present in different genomes is a reliable predictor of phase variation. overall, putative phase variable genes were reported in the neisseria genus, of which were specific to the nm species and included genes involved in cell adhesion, capsule formation, biosynthesis of the lipopolysaccharide layer, cell-surface receptors and restriction/modification systems. among these, phase variation was confirmed experimentally for only genes, primarily combining pcr with immunoblotting techniques or with northern blotting and quantitative pcr. here, we report that the limited number of sequences employed in previous vnssrs studies, together with the lack of high-throughput technologies suitable for a large-scale validation of putative phase variable genes, have potentially led to an underestimation of the overall impact of ssr-mediated phase variation on nm phenotype and fitness. also, the evolutionary dynamics underlying the generation of vnssrs in nm and other species has not been fully elucidated. recently, new insights into the nm population structure and dynamics have been generated through a comparative genomic analysis based on the complete genome sequences of 20 nm strains, including multiple isolates from each of the most virulent clonal complexes. also, the advent of ngs technologies have allowed for an unprecedented sequencing depth across the whole genome, such that the presence of mixed populations at specific loci of a single genome can be detected with high throughput, including those variants occurring at low frequencies during bacterial replication. in the present study we took advantage of nm pan-genome data to map, characterize and infer functional impact and evolutionary properties of ssr contingency loci, and ngs to experimentally validate in silico predictions on stochastic genotypic switching.  two hundred seventy-seven meningococcal genes are associated with simple sequence repeats that show inter-strain length polymorphisms an average of ssrs were identified in each genome, the majority of which were represented by homopolymeric tracts. the number of ssrs was similar across genomes and their abundance was found to be significantly higher than random expectation. significant deviation from neutrality was also observed for each individual ssr type. we identified clusters of orthologous ssrs, 35 % of which were represented in every genome analyzed. among these, clusters had polymorphisms in the repeat tract length across different genomes, located either in intragenic regions or in the associated gene’s upstream region. these were selected for further analyses based on their potential to change gene expression. consistent with previous findings, homopolymeric tracts were the most abundant simple sequence repeat type, with vnssrs clusters divided into a/t and 59 g/c tracts. the remaining clusters were represented by tandem repeats of two- to nine-nucleotides motifs. positional analysis identified intragenic and intergenic vnssrs, respectively. interestingly, all vnssrs with a unit motif of three-, six- and nine-nucleotides, for which a length polymorphism would not result in the disruption of the reading frame, were found to be located within coding sequences. finally, thirty-two genes were associated with more than a single vnssr, resulting in a total of phase variable gene candidates. ssr contingency loci show different genotype switching frequencies based on repeat variability and intra-strain phylogenetic relationships we stratified the vnssr-associated genes into strong, moderate and weak candidates for phase variation. most weak candidate phase variable genes were associated with a/t homopolymeric tracts while g/c repeats were the most abundant vnssr type among the strong candidate genes. vnssrs other than homopolymeric repeats were found in moderate and strong candidate genes, respectively.fig. properties of identified vnssrs. a relative contribution of different vnssrs types to the strong, moderate and weak candidate phase variable genes. other: all non-homopolymeric repeats. g/c: guanine or cytosine variable homopolymeric tracts. a/t: adenine or thymine variable homopolymeric tracts. b nei’s diversity index distribution of vnssrs associated with the strong, moderate and weak candidate phase variable genes. genes associated with intragenic and intergenic vnssrs are represented full and empty circles, respectively. the thick horizontal lines represent distributions median. p: wilcoxon-mann-whitney p value intrinsic variability of vnssrs clusters, defined as the nei’s diversity index, was then computed for the three gene sets. vnssrs associated to strong candidate genes were significantly more variable than those associated to moderate genes, which in turn showed higher variability than those linked to weak phase variable gene candidates, providing supporting evidence that genes associated with a different likelihood of undergoing phase variation exist and that these can be identified through the analysis of a suitable genome dataset. a positive correlation between the ssr tract length and repeat instability has been proposed, possibly as a result of a decreased proof-reading efficiency of the dna polymerase over the longer tracts. we could confirm such association for the strong but not for the moderate or weak candidate gene pools. vnssrs variation is detectable in single colonies grown overnight under non-selective conditions the genomes of five of the analyzed strains were re-sequenced at high depth using a next generation sequencing platform with the goal of detecting repeats length polymorphisms occurring in single colonies during an overnight growth. analysis of the sequencing reads confirmed repeat length polymorphisms in of the predicted phase variable genes, distributed in strong, moderate and weak putative phase variable genes. the same procedure applied to randomly chosen, non-variable ssrs detected length polymorphisms in loci, a significantly lower proportion compared to vnssrs. additionally, the average frequencies of polymorphic reads observed for vnssrs was  %, a higher proportion compared to the false positive ssrs or to the allowed read mapping error rate. frequency distributions of polymorphic reads identified in the genomes are reported in additional file 2: figure s fig. phase variable genes confirmed by illumina sequencing. numbers inside bars indicate the frequency of phase variable genes present in each category. confirmed: predicted phase variable genes confirmed by deep sequencing data. not confirmed: predicted phase variable genes not showing evidence of variation. control: analysis of randomly selected ssrs that were predicted as non-variable by the genomic comparative analysis seventy of the validated vnssrs were present in all five re-sequenced genomes. among these, were shown to be variable in all the genomes analyzed, while were confirmed in four or fewer genomes. interestingly, for any one genome there was a strong association between ssr size and the likelihood of finding it variable. specifically, homopolymeric tracts of or more nucleotides were predominantly found to be variable in four or five genomes. because there were only three non-homopolymeric tracts, it was not possible to generalize this conclusion to longer repeats. additionally, for the experimentally confirmed core vnssrs, the proportion of polymorphic reads was found to correlate with the nei’s diversity index derived from the comparison of the available genome sequences. specifically, vnssrs associated to a nei’s index > were mainly characterized by above-average polymorphic reads proportions, meaning that the extent of vnssrs variability can be predicted with good confidence by the comparison of multiple genomic sequences. finally, the number of phase variable genes that would be confirmed by deep sequencing of a wider genome collection was estimated by applying an approach originally proposed for estimating the bacterial pan-genome size. the regression analysis based on the number of validated phase variable genes predicted that if all the genomes were re-sequenced, the number of validated phase-variable genes would have been 146 ± , while the analysis of genomes could raise the number to 190 ± . the pool of genes whose expression is influenced by vnssrs is wider than previously hypothesized the pool of putative phase variable genes predicted in this study was compared with the genes previously described or predicted to be phase-variable in nm. our approach confirmed of the genes whose phase variation had been previously demonstrated in nm, as well as of the putative phase variable genes previously proposed as strong candidates and of the genes previously reported as either moderate or weak candidates. the genes not confirmed by our method, either didn’t satisfy the criteria for ssr search or didn’t show inter-strain length polymorphisms. in the particular case of a poly-c tract located within the nmb coding region, the repeat showed between-strain variation, alternating between the c and c states. however, since the c was below the minimum length cut off our approach failed to identify this variable ssr. remarkably, the wider genome collection allowed for the identification of new putative phase variable genes, predicted as strong, moderate and weak candidates respectively. candidate phase variable genes related to cell surface structure or involved in dna metabolism are overrepresented the strong and moderate candidate genes, described in results section distributed across different functional roles. consistent with previous studies vnssr-associated genes were found to be primarily involved in host-microbial interplay, with “cell envelope”, “transport and binding” and “dna metabolism” functional roles being significantly overrepresented. intragenic and intergenic vnssrs were both present in the enriched categories. other biological processes included functions related to the peptide synthesis machinery. additionally, putative phase variable genes were found in functional categories that had not been previously associated to ssr-mediated regulation, such as genes encoding for transposons and prophage related functions, transcription factors and proteins participating in nucleic acids biosynthesis.fig. tigr functional roles represented by the strong and moderate candidate phase variable genes. x-axis represents the proportion of putative phase variable genes present in each annotation. numbers in brackets represent the number of putative phase variable genes over the total number of genes that are associated to a specific function. genes associated with intragenic and intergenic vnssrs are represented in dark and light grey, respectively. ‘*’: enriched functional roles. ‘+’: functional roles not previously described to be associated with ssr-mediated regulation in nm the weak candidate genes were associated with different functional roles, fourteen of which were common to the strong and moderate groups. this subset showed a different profile, whereby genes were significantly overrepresented in the “hypothetical proteins” and “mobile and extrachromosomal element functions” functional annotations. most vnssrs are predicted to induce phenotypic changes by modulating the level of gene expression rather than mediating on/off translational switching we sought to gather new insights into the vnssr-mediated regulatory mechanisms by analyzing the repeat sequence context in the strong candidate phase variable genes. nineteen of the vnssrs associated with strong candidate genes were surrounded by multiple in-frame 5′-atg-3′ translational initiation codons. these vnssrs can shift the origin of translation between alternative start sites, resulting in the modulation of gene expression. such modulation can arise either from the switching between start sites associated with different expression levels or from the production of different protein isoforms, as may be the case for the nmb and nmb open reading frames, respectively. different peptide isoforms may also result from vnssrs located in the 3′ portion of the gene. in this location, the vnssr can bring the c-terminal peptide portion out of frame while still allowing the translation of a functional protein. two orfs, nmb and nmb were consistent with this mechanism. in contrast, vnssrs present within the nmb nmb and nmb reading frames were either a 3- or a 9-nucleotides repeat. these vnssrs do not interfere with the reading frame but can still influence the cell phenotype by producing proteins with different structures and/or post-translational profiles. also, of the analyzed vnssrs were predicted to interfere with the associated gene promoters, as represented by the t-homopolymeric tract located upstream of the nmb start. the remaining repeats were predicted to introduce frameshifts likely resulting in on-off type of transcriptional regulation, as proposed for nmb.fig. schematic representation of five vnssrs and their sequence context. a vnssr causing translational start site switching. b vnssr causing the loss of a membrane-spanning domain. c vnssr leading to the loss of the peptide c-terminal region. d vnssr introducing changes in the peptide sequence. e) vnssr influencing the gene promoter. f vnssr introducing an inactivating frameshift. dark grey arrows represent open reading frames. black arrows marked with atg represent in-frame atg translational start sites. light grey boxes represent the annotated functional domains. stripped boxes represent vnssrs and the related tags indicate the repeat unit motif along with the minimum and maximum number of repetitions observed in the analyzed genomes. numbers below each gene indicate the position relative to the annotated translational tart site most strong candidate phase variable genes are not specific to n. meningitidis among the strong candidate phase variable genes, were identified only within nm through sequence database comparison. the gc content of these genes largely deviates from the typical meningococcal value, suggesting acquisition through lateral gene transfer from outside the genus. conversely, orthologous copies of the remaining genes were also identified in other neisseria species. among these, were associated with the same ssrs observed in nm, indicating that these ssr-contingency loci arose before meningococcus speciation. in the remaining four cases the ssr associated with the gene was detected in nm only, suggesting that these contingency loci have evolved following nm speciation. orthologous copies of strong candidate phase variable genes were also identified in species outside of the neisseria genus, including cardiobacterium hominis, haemophilus aegyptius, haemophilus haemolyticus, haemophilus influenzae, kigella dentrificans and streptococcus pneumoniae. among these, the nmb encoded gene identified in the three haemophilus species and the opacity proteins identified in s. pneumoniae, were associated with the same ssr as in nm, raising the hypothesis that ssr contingency loci are transferred across species through lateral dna exchange. vnssrs reorganization following n. meningitidis speciation primarily involved cell surface determinants among the vnssrs associated to the strong candidate genes identified in multiple taxa, were found to be significantly longer in nm than in other species, while were found to be significantly shorter. twelve and of the repeats with increased length in nm were g/c and a/t homopolymeric tracts, respectively, along with other non-homopolymeric ssrs. these repeats were primarily associated with cell surface determinants, including surface receptors, surface transporters, pili and lipopolysaccharide biosynthetic enzymes. the vnssrs with a reduced size in nm included a poly-t and a poly-c homopolymeric tracts, respectively associated with a hypothetical protein and to a ferric enterobactin receptor, and a 5′-cttct-3′ repeat associated to an opacity protein. also in this case, an association with surface exposed proteins was observed. overall, these results suggest that evasion of host immune responses, along with changes in surface adhesion properties and tissue tropism, have likely been the main drivers of ssr contingency loci evolution within nm. vnssrs chromosomic loci are fixed in the meningococcal population chromosomic loci containing the identified vnssrs were tested for signatures of selective pressures that could confirm a beneficial role of those loci and justify their over-abundance in nm species. only vnssr loci, all associated with a/t homopolymeric tracts, showed a significant deviation from neutrality with all the applied tests. signatures of positive selection were identified in the upstream region of the two strong candidate phase variable genes mtrf, a surface-exposed efflux-pump component involved in hydrophobic antimicrobial resistance, and the transcription factor dksa. a signature of negative selection was detected in the upstream region of the moderate candidate phase variable gene sstt, which also encodes for a surface transporter. all other signals were identified in vnssrs associated to weak candidate phase variable genes. the remaining tested loci did not show signatures of selection, supporting the hypothesis that ssr contingency loci have reached an equilibrium state and are now fixed in the population. discussion after the first nm complete genome sequence became available, the highly abundant ssrs were largely associated with genes involved in host adaptation, commensalism and virulence. these initial studies on nm ssrs contingency loci were based on analysis of individual genomes, strains mc and z. with the technology available at that time, the generation of new genomic sequences and the validation of the newly proposed phase-variable genes were laborious and time-consuming tasks. despite the challenges, a few years after the first meningococcal genome sequence was revealed, a study reported a list of 69 nm phase variable genes, of which were experimentally confirmed. benefiting from the most recent sequencing technologies, the present study employed a 20 nm genomes dataset, with a resolved phylogenetic structure, in a comparative genomic analysis aimed at the study of ssr-driven phase variation as a population-based, as opposed to single cell, adaptive strategy. all analyzed strains were found to have more ssrs than expected in their chromosomes. a comparable number of ssrs were also found in two non-pathogenic strains, α and ox99- questioning the hypothesis that ssrs-mediated phase variation is an adaptation associated with the evolution of virulence. most of the identified ssrs were homopolymeric tracts, with an excess of a/t tandem repeats, suggesting slipped-strand mispairing as the primary mechanism driving the occurrence of ssrs. it has been proposed that slipped-strand mispairing is preferentially associated with homopolymeric tracts due to their high sequence redundancy and favouring a/t repeats due to the lower stability of the dna-dna hybrid during bacterial genome replication. the likelihood that repeat-associated genes undergo phase variation is variable. efforts to identify molecular signatures indicative of such likelihood have defined a number of cis factors, including repeat sequence, size and sequence context, along with the associated gene function. particularly, genes encoding for surface antigens were regarded as strong candidates for phase variation. we applied a more unsupervised approach, in which only the repeat sequence context and the distribution of its variation over the sampled phylogeny were considered. the over-representation of a/t repeats and their alleged higher propensity to cause transcriptional slippage compared to g/c repeats would suggest these to be the main contributor of ssr mediated phase variation. conversely, in agreement with previous findings, homopolymeric g/c repeats were the most frequent among the strong candidate genes pool, while weak candidate genes were almost exclusively associated with a/t repeats. a primary involvement of c/g homopolymeric in the mediation of phase variation has also been described in c. jejuni. this particular mutational pattern, together with the fact that we didn’t identify traces of selective pressure for most vnssrs, suggest that there may be molecular drivers governing the mutability of these contingency loci. in support of this hypothesis, increased rates of phase variation have been associated with the loss of the dam dna methylase and the two mismatch repair proteins muts and mutl. regarding non-homopolymeric repeats, 3-, 4-, 5-, 6- and 9-nucelotides vnssrs were also found to be abundant. three-, 6- and 9-nucleotides vnssrs do not disrupt the reading frame and their abundance may be justified by a reduced selective pressure for their stability. abundance of 5-nucleotide vnssrs is due to the fact that out of cases are represented by the same 5′-cttct-3′ repeat associated to different paralogs of an opacity protein. interestingly, a homolog of this opacity protein, coupled with the same ssr, was also identified in streptococcus pneumoniae genome, providing evidence that ssr-contingency loci can be transferred across different species, or strains of the same species and even different chromosomic loci of the same strain through lateral gene transfer. conversely, abundance of 4-nucleotides vnssrs, which has been reported both in nm and in haemophilus influenzae, has no obvious interpretation. given the conservation of those repeat sequences and their proximity to the translational origin, the interaction with yet unidentified regulatory elements may be involved. one such regulatory mechanism has been described for the taaa repeat associated with the nada promoter. the availability of high-throughput sequencing data enabled a high throughput experimental validation of in silico predicted ssr contingency loci. the rationale behind this approach was that dna libraries used for illumina sequencing were generated from a bacterial culture derived from a single colony. after an 18 h growth and approximately bacterial replication cycles, a library was expected to contain, and represent in the form of variable alleles, most of the vnssrs variation. this test provided circumstantial evidence for phase variability of and 27 % of the strong, moderate and weak candidate genes respectively, leading to a total of experimentally confirmed phase variable genes. this approach also allowed estimating vnssrs variation frequencies. the validated vnssr showed an average proportion of alternative alleles of  %, indicating an average mutation rate of once every seven cell replication cycles. however, given the contribution of environmental, population and molecular factors to the mutability of ssrs in phase variable loci, it is likely that real frequencies may differ from those we observed in vitro. ssr contingency loci were over-represented among genes related to cell surface functions and structure or involved in dna metabolism, primarily dna restriction and modification. these results, which are in agreement with previous findings, reflect a primary involvement of ssr-mediated phase variation in microbe-host interaction. dis-regulation of those genes, in fact, alters nm ability to adhere to and invade host cells, its ability to scavenge required nutrients and to evade host immune defenses by generating antigenic variants. additionally, putative phase variable genes were found in several other functional classes, indicating that ssr-mediated gene regulation interferes with a set of metabolic functions wider than previously hypothesized. our attempt to understand how vnssrs have evolved in nm revealed that most strong candidate phase variable genes are widespread to the neisseria genus and that about half of them have undergone size change after meningococcal speciation. this is suggestive of the fact that most ssr-associated contingency loci have developed prior to nm speciation and have subsequently evolved, possibly as a consequence of the process of adaptation to the human host environment. in agreement with this is the fact that identified phase variable genes were enriched for cell surface determinants, suggesting that evasion of host immune responses and changes in surface adhesion properties and tissue tropism are the main drivers of ssr contingency loci evolution in nm. regarding the remaining genes that are specific to nm, they were characterized by a g/c content that deviates from that of meningococcus, suggesting a possible acquisition through horizontal gene transfer. phase variable genes candidates were also identified in species outside the neisseria genus. interestingly, all of them are human pathogens that coexist in the respiratory flora of healthy individuals. moreover, we found that genes shared by nm, s. pneumoniae and three haemophilus members were associated with the same repeat element as in nm, supporting the possibility that ssr contingency loci are exchanged across species through lateral dna exchange. overall, our findings suggest that ssr contingency loci have been acquired by nm either by spontaneous evolution, inheritance from an ancestral species or acquisition by means of lateral dna exchange. in the context of a commensal organism, constantly interacting with its host, the onset of new mutations or allele variants may lead to an increased fitness, a reduced fitness or have no effect on fitness. one of the most accredited hypotheses is that vnssrs have evolved as an adaptive strategy, in organisms like nm, to increase their fitness in the hostile and mutable host environment. we tested all vnssr-containing loci for signatures of selection that could confirm a beneficial role of those loci and justify their over-abundance. surprisingly, the analysis revealed that most ssr-contingency loci are not under selective pressure, suggesting that nm has had enough time to adapt to its niche and has apparently reached an equilibrium state. it must be noted, however, that substantial homologous recombination and repeated population expansions and bottlenecks, along with possible sampling biases, are likely to confound or even obscure neutrality tests results. nonetheless, we envision that, applying this approach to a wider, unbiased strain collection, coming from a single outbreak, will have the potential of providing valuable insights into the understanding of vnssrs evolutionary processes. despite a proportion of coding dna in nm higher than 80 %, we observed that 51 % of the identified vnssrs are found within intergenic regions. this bias toward intergenic vnssrs, which may be explained by a reduced selective pressure along non-coding regions, can also be indicative of a regulation of gene expression occurring at the transcriptional level. this type of vnssrs can indeed modulate promoter efficiency by changing the relative spacing between key promoter elements. differently, intragenic vnssrs, which can induce frameshift mutations, are responsible for an on/of type of regulation occurring at the translational level. additionally, given the existence of alternative translational start codons, intragenic vnssrs located in the 5′ portion of a gene can switch between multiple active reading frames, with the potential of affecting the level of gene expression or producing alternative peptide isoforms. in the present analysis we observed that most intragenic vnssrs are located in the proximity of the translational start site and that among the strong candidate phase variable genes, there are cases in which the vnssr was found between two or more in-frame atg translational start sites. twenty other vnssrs were predicted to either interact with the gene promoter or to produce alternative peptides isoforms. overall, collected evidences suggest that ssr-contingency loci have evolved in such a way to allow the maximum degree of flexibility, reflected by their ability to finely modulate gene expression level, or the functionality of the resulting protein, rather than merely switching between the two extreme cases of a gene being expressed or silenced through the introduction of frameshifts.  in conclusion, this analysis, which capitalized on the most recent dna sequencing technologies, confirmed the fundamental contribution of ssr contingency loci in promoting phenotypic variation and their deep implications in nm survival and adaptation to the surrounding environment. our unsupervised approach allowed the generation of a panel of genes whose expression may be controlled or influenced by ssr elements and to provide corroborative evidence for the phase variability for of them. functional characterization of these genes highlighted an enrichment for cell surface determinants, which included members of the adhesins, evasins, lipopolisaccharide biosynthesis and nutrient scavenging protein families. such proteins have been attributed multiple roles in attachment of bacterial cells to host membranes, in regulating bacterial resistance to both innate and adaptive immune system and as major determinants of meningococcal invasive disease. ssr-mediated phase variation has also been reported to control the expression of multiple immunogens currently included in commercial meningococcal vaccines and to modulate resistance to antimicrobial agents. we therefore envision that future studies on meningococcal microbe-host interaction and studies aimed at the identification of new vaccine antigens and antimicrobial molecules will benefit from this comprehensive characterization of the meningococcal putative phase variable genes repertoire. 
however, in addition we identified an unexpectedly broad spectrum of other metabolic functions, and most of the variable ssrs were predicted to induce phenotypic changes by modulating gene expression at a transcriptional level or by producing different protein isoforms rather than mediating on/off translational switching through frameshifts. investigation of the evolutionary history of ssr contingency loci revealed that these loci were inherited from a nm ancestor, evolved independently within nm, or were acquired by nm through lateral dna exchange. deep-coverage dna sequencing of single colonies grown overnight under non-selective conditions confirmed the presence of high-frequency, stochastic variation in of them, providing circumstantial evidence for their phase variability. overall, our results have identified a broader and qualitatively different phenotypic diversification of ssrs-mediated stochastic variation than previously documented, including its impact on central nm metabolism. the neisseria meningitidis chromosome shows a high abundance of simple sequence dna repeats that undergo stochastic, reversible mutations at high frequency. we confirmed previous observations of a predominance of variable ssrs within genes for components located on the cell surface or dna metabolism. here we present a population-scale comparative genomic analysis that identified genes and classified them into strong, moderate and weak candidates for phase variation. neisseria meningitidisphase variationcontingency locicomparative genomicshost-pathogen interactionissue-copyright-statement© the author 2016
 the neisseria meningitidis chromosome shows a high abundance of simple sequence dna repeats that undergo stochastic, reversible mutations at high frequency. this mechanism is reflected in an extensive phenotypic diversity that facilitates nm adaptation to dynamic environmental changes. to date, phase-variable phenotypes mediated by ssrs variation have been experimentally confirmed for 26 nm genes.  here we present a population-scale comparative genomic analysis that identified genes and classified them into strong, moderate and weak candidates for phase variation. deep-coverage dna sequencing of single colonies grown overnight under non-selective conditions confirmed the presence of high-frequency, stochastic variation in of them, providing circumstantial evidence for their phase variability. we confirmed previous observations of a predominance of variable ssrs within genes for components located on the cell surface or dna metabolism. however, in addition we identified an unexpectedly broad spectrum of other metabolic functions, and most of the variable ssrs were predicted to induce phenotypic changes by modulating gene expression at a transcriptional level or by producing different protein isoforms rather than mediating on/off translational switching through frameshifts. investigation of the evolutionary history of ssr contingency loci revealed that these loci were inherited from a nm ancestor, evolved independently within nm, or were acquired by nm through lateral dna exchange.  overall, our results have identified a broader and qualitatively different phenotypic diversification of ssrs-mediated stochastic variation than previously documented, including its impact on central nm metabolism. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords neisseria meningitidisphase variationcontingency locicomparative genomicshost-pathogen interactionissue-copyright-statement© the author 2016 