Bacterial contingency loci: the role of simple sequence DNA repeats in bacterial adaptation

Investigation of a Quadruplex-Forming Repeat Sequence Highly Enriched in Xanthomonas and Nostoc sp

In prokaryotes simple sequence repeats (SSRs) with unit sizes of 1–5 nucleotides (nt) are causative for phase and antigenic variation. Although an increased abundance of heptameric repeats was noticed in bacteria, reports about SSRs of 6–9 nt are rare. In particular G-rich repeat sequences with the propensity to fold into G-quadruplex (G4) structures have received little attention. In silico analysis of prokaryotic genomes show putative G4 forming sequences to be abundant. This report focuses on a surprisingly enriched G-rich repeat of the type GGGNATC in Xanthomonas and cyanobacteria such as Nostoc. We studied in detail the genomes of Xanthomonas campestris pv. campestris ATCC 33913 (Xcc), Xanthomonas axonopodis pv. citri str. 306 (Xac), and Nostoc sp. strain PCC7120 (Ana). In all three organisms repeats are spread all over the genome with an over-representation in non-coding regions. Extensive variation of the number of repetitive units was observed with repeat numbers ranging from two up to 26 units. However a clear preference for four units was detected. The strong bias for four units coincides with the requirement of four consecutive G-tracts for G4 formation. Evidence for G4 formation of the consensus repeat sequences was found in biophysical studies utilizing CD spectroscopy. The G-rich repeats are preferably located between aligned open reading frames (ORFs) and are under-represented in coding regions or between divergent ORFs. The G-rich repeats are preferentially located within a distance of 50 bp upstream of an ORF on the anti-sense strand or within 50 bp from the stop codon on the sense strand. Analysis of whole transcriptome sequence data showed that the majority of repeat sequences are transcribed. The genetic loci in the vicinity of repeat regions show increased genomic stability. In conclusion, we introduce and characterize a special class of highly abundant and wide-spread quadruplex-forming repeat sequences in bacteria.

OxyR-dependent formation of DNA methylation patterns in OpvABOFF and OpvABON cell lineages of Salmonella enterica

Phase variation of the Salmonella enterica opvAB operon generates a bacterial lineage with standard lipopolysaccharide structure (OpvAB(OFF)) and a lineage with shorter O-antigen chains (OpvAB(ON)). Regulation of OpvAB lineage formation is transcriptional, and is controlled by the LysR-type factor OxyR and by DNA adenine methylation. The opvAB regulatory region contains four sites for OxyR binding (OBSA-D), and four methylatable GATC motifs (GATC1-4). OpvAB(OFF) and OpvAB(ON) cell lineages display opposite DNA methylation patterns in the opvAB regulatory region: (i) in the OpvAB(OFF) state, GATC1 and GATC3 are non-methylated, whereas GATC2 and GATC4 are methylated; (ii) in the OpvAB(ON) state, GATC2 and GATC4 are non-methylated, whereas GATC1 and GATC3 are methylated. We provide evidence that such DNA methylation patterns are generated by OxyR binding. The higher stability of the OpvAB(OFF) lineage may be caused by binding of OxyR to sites that are identical to the consensus (OBSA and OBSc), while the sites bound by OxyR in OpvAB(ON) cells (OBSB and OBSD) are not. In support of this view, amelioration of either OBSB or OBSD locks the system in the ON state. We also show that the GATC-binding protein SeqA and the nucleoid protein HU are ancillary factors in opvAB control.

Epigenetic Control of Salmonella enterica O-Antigen Chain Length: A Tradeoff between Virulence and Bacteriophage Resistance

The Salmonella enterica opvAB operon is a horizontally-acquired locus that undergoes phase variation under Dam methylation control. The OpvA and OpvB proteins form intertwining ribbons in the inner membrane. Synthesis of OpvA and OpvB alters lipopolysaccharide O-antigen chain length and confers resistance to bacteriophages 9NA (Siphoviridae), Det7 (Myoviridae), and P22 (Podoviridae). These phages use the O-antigen as receptor. Because opvAB undergoes phase variation, S. enterica cultures contain subpopulations of opvABOFF and opvABON cells. In the presence of a bacteriophage that uses the O-antigen as receptor, the opvABOFF subpopulation is killed and the opvABON subpopulation is selected. Acquisition of phage resistance by phase variation of O-antigen chain length requires a payoff: opvAB expression reduces Salmonella virulence. However, phase variation permits resuscitation of the opvABOFF subpopulation as soon as phage challenge ceases. Phenotypic heterogeneity generated by opvAB phase variation thus preadapts Salmonella to survive phage challenge with a fitness cost that is transient only.

Accelerating Discovery and Functional Analysis of Small RNAs with New Technologies

Over the past decade, bacterial small RNAs (sRNAs) have gone from a biological curiosity to being recognized as a major class of regulatory molecules. High-throughput methods for sampling the transcriptional output of bacterial cells demonstrate that sRNAs are universal features of bacterial transcriptomes, are plentiful, and appear to vary extensively over evolutionary time. With ever more bacteria coming under study, the question becomes how can we accelerate the discovery and functional characterization of sRNAs in diverse organisms. New technologies built on high-throughput sequencing are emerging that can rapidly provide global insight into the numbers and functions of sRNAs in bacteria of interest, providing information that can shape hypotheses and guide research. In this review, we describe recent developments in transcriptomics (RNA-seq) and functional genomics that we expect to help us develop an integrated, systems-level view of sRNA biology in bacteria.

The external domains of the HIV-1 envelope are a mutational cold spot

In RNA viruses, mutations occur fast and have large fitness effects. While this affords remarkable adaptability, it can also endanger viral survival due to the accumulation of deleterious mutations. How RNA viruses reconcile these two opposed facets of mutation is still unknown. Here we show that, in human immunodeficiency virus (HIV-1), spontaneous mutations are not randomly located along the viral genome. We find that the viral mutation rate experiences a threefold reduction in the region encoding the most external domains of the viral envelope, which are strongly targeted by neutralizing antibodies. This contrasts with the hypermutation mechanisms deployed by other, more slowly mutating pathogens such as DNA viruses and bacteria, in response to immune pressure. We show that downregulation of the mutation rate in HIV-1 is exerted by the template RNA through changes in sequence context and secondary structure, which control the activity of apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3 (A3)-mediated cytidine deamination and the fidelity of the viral reverse transcriptase.

Mutations allow RNA virus to adapt fast but also entail fitness costs. Geller et al. show that, in HIV-1, mutations occur three times less often in the most external domains of the envelope, and that this is due to changes in RNA sequence context and structure, which control viral and host-encoded mutational mechanisms.

When natural selection gives gene function the cold shoulder

It is tempting to invoke organismal selection as perpetually optimizing the function of any given gene. However, natural selection can drive genic functional change without improvement of biochemical activity, even to the extinction of gene activity. Detrimental mutations can creep in owing to linkage with other selectively favored loci. Selection can promote functional degradation, irrespective of genetic drift, when adaptation occurs by loss of gene function. Even stabilizing selection on a trait can lead to divergence of the underlying molecular constituents. Selfish genetic elements can also proliferate independent of any functional benefits to the host genome. Here we review the logic and evidence for these diverse processes acting in genome evolution. This collection of distinct evolutionary phenomena - while operating through easily understandable mechanisms - all contribute to the seemingly counterintuitive notion that maintenance or improvement of a gene's biochemical function sometimes do not determine its evolutionary fate.

Stochastic Switching of Cell Fate in Microbes

Microbes transiently differentiate into distinct, specialized cell types to generate functional diversity and cope with changing environmental conditions. Though alternate programs often entail radically different physiological and morphological states, recent single-cell studies have revealed that these crucial decisions are often left to chance. In these cases, the underlying genetic circuits leverage the intrinsic stochasticity of intracellular chemistry to drive transition between states. Understanding how these circuits transform transient gene expression fluctuations into lasting phenotypic programs will require a combination of quantitative modeling and extensive, time-resolved observation of switching events in single cells. In this article, we survey microbial cell fate decisions demonstrated to involve a random element, describe theoretical frameworks for understanding stochastic switching between states, and highlight recent advances in microfluidics that will enable characterization of key dynamic features of these circuits.

Occurrence, distribution and possible functional roles of simple sequence repeats in phytoplasma genomes

Phytoplasmas are unculturable, cell-wall-less bacteria that parasitize plants and insects. This transkingdom life cycle requires rapid responses to vastly different environments, including transitions from plant phloem sieve elements to various insect tissues and alternations among diverse plant hosts. Features that enable such flexibility in other microbes include simple sequence repeats (SSRs) — mutation-prone, phase-variable short DNA tracts that function as ‘evolutionary rheostats’ and enhance rapid adaptations. To gain insights into the occurrence, distribution and potentially functional roles of SSRs in phytoplasmas, we performed computational analysis on the genomes of five completely sequenced phytoplasma strains, ‘Candidatus Phytoplasma asteris’-related strains OYM and AYWB, ‘Candidatus Phytoplasma australiense’-related strains CBWB and SLY and ‘Candidatus Phytoplasma mali’-related strain AP-AT. The overall density of SSRs in phytoplasma genomes was higher than in representative strains of other prokaryotes. While mono- and trinucleotide SSRs were significantly overrepresented in the phytoplasma genomes, dinucleotide SSRs and other higher-order SSRs were underrepresented. The occurrence and distribution of long SSRs in the prophage islands and phytoplasma-unique genetic loci indicated that SSRs played a role in compounding the complexity of sequence mosaics in individual genomes and in increasing allelic diversity among genomes. Findings from computational analyses were further complemented by an examination of SSRs in varied additional phytoplasma strains, with a focus on potential contingency genes. Some SSRs were located in regions that could profoundly alter the regulation of transcription and translation of affected genes and/or the composition of protein products.

A biphasic epigenetic switch controls immunoevasion, virulence and niche adaptation in non-typeable Haemophilus influenzae

Non-typeable Haemophilus influenzae contains an N⁶-adenine DNA-methyltransferase (ModA) that is subject to phase-variable expression (random ON/OFF switching). Five modA alleles, modA2, modA4, modA5, modA9 and modA10, account for over two-thirds of clinical otitis media isolates surveyed. Here, we use single molecule, real-time (SMRT) methylome analysis to identify the DNA-recognition motifs for all five of these modA alleles. Phase variation of these alleles regulates multiple proteins including vaccine candidates, and key virulence phenotypes such as antibiotic resistance (modA2, modA5, modA10), biofilm formation (modA2) and immunoevasion (modA4). Analyses of a modA2 strain in the chinchilla model of otitis media show a clear selection for ON switching of modA2 in the middle ear. Our results indicate that a biphasic epigenetic switch can control bacterial virulence, immunoevasion and niche adaptation in an animal model system.

Non-typeable Haemophilus influenzae, which causes ear and lung infections, has a DNA methyltransferase encoded by alternative alleles that are subject to random ON/OFF switching. Here, Atack et al. show that this epigenetic switch controls the expression of key proteins involved in virulence.

Chlamydia trachomatis In Vivo to In Vitro Transition Reveals Mechanisms of Phase Variation and Down-Regulation of Virulence Factors

Research on the obligate intracellular bacterium Chlamydia trachomatis demands culture in cell-lines, but the adaptive process behind the in vivo to in vitro transition is not understood. We assessed the genomic and transcriptomic dynamics underlying C. trachomatis in vitro adaptation of strains representing the three disease groups (ocular, epithelial-genital and lymphogranuloma venereum) propagated in epithelial cells over multiple passages. We found genetic features potentially underlying phase variation mechanisms mediating the regulation of a lipid A biosynthesis enzyme (CT533/LpxC), and the functionality of the cytotoxin (CT166) through an ON/OFF mechanism. We detected inactivating mutations in CT713/porB, a scenario suggesting metabolic adaptation to the available carbon source. CT135 was inactivated in a tropism-specific manner, with CT135-negative clones emerging for all epithelial-genital populations (but not for LGV and ocular populations) and rapidly increasing in frequency (~23% mutants per 10 passages). RNA-sequencing analyses revealed that a deletion event involving CT135 impacted the expression of multiple virulence factors, namely effectors known to play a role in the C. trachomatis host-cell invasion or subversion (e.g., CT456/Tarp, CT694, CT875/TepP and CT868/ChlaDub1). This reflects a scenario of attenuation of C. trachomatis virulence in vitro, which may take place independently or in a cumulative fashion with the also observed down-regulation of plasmid-related virulence factors. This issue may be relevant on behalf of the recent advances in Chlamydia mutagenesis and transformation where culture propagation for selecting mutants/transformants is mandatory. Finally, there was an increase in the growth rate for all strains, reflecting gradual fitness enhancement over time. In general, these data shed light on the adaptive process underlying the C. trachomatis in vivo to in vitro transition, and indicates that it would be prudent to restrict culture propagation to minimal passages and check the status of the CT135 genotype in order to avoid the selection of CT135-negative mutants, likely originating less virulent strains.

Bioinformatic Analysis of Chlamydia trachomatis Polymorphic Membrane Proteins PmpE, PmpF, PmpG and PmpH as Potential Vaccine Antigens

Chlamydia trachomatis is the most important infectious cause of infertility in women with important implications in public health and for which a vaccine is urgently needed. Recent immunoproteomic vaccine studies found that four polymorphic membrane proteins (PmpE, PmpF, PmpG and PmpH) are immunodominant, recognized by various MHC class II haplotypes and protective in mouse models. In the present study, we aimed to evaluate genetic and protein features of Pmps (focusing on the N-terminal 600 amino acids where MHC class II epitopes were mapped) in order to understand antigen variation that may emerge following vaccine induced immune selection. We used several bioinformatics platforms to study: i) Pmps' phylogeny and genetic polymorphism; ii) the location and distribution of protein features (GGA(I, L)/FxxN motifs and cysteine residues) that may impact pathogen-host interactions and protein conformation; and iii) the existence of phase variation mechanisms that may impact Pmps' expression. We used a well-characterized collection of 53 fully-sequenced strains that represent the C. trachomatis serovars associated with the three disease groups: ocular (N=8), epithelial-genital (N=25) and lymphogranuloma venereum (LGV) (N=20). We observed that PmpF and PmpE are highly polymorphic between LGV and epithelial-genital strains, and also within populations of the latter. We also found heterogeneous representation among strains for GGA(I, L)/FxxN motifs and cysteine residues, suggesting possible alterations in adhesion properties, tissue specificity and immunogenicity. PmpG and, to a lesser extent, PmpH revealed low polymorphism and high conservation of protein features among the genital strains (including the LGV group). Uniquely among the four Pmps, pmpG has regulatory sequences suggestive of phase variation. In aggregate, the results suggest that PmpG may be the lead vaccine candidate because of sequence conservation but may need to be paired with another protective antigen (like PmpH) in order to prevent immune selection of phase variants.

Comparative analysis of the Mycoplasma capricolum subsp. capricolum GM508D genome reveals subrogation of phase-variable contingency genes and a novel integrated genetic element

Mycoplasma capricolum subspecies capricolum is both a pathogen of small ruminants and a model recipient organism for gene transplantation and synthetic biology. With the availability of the complete genome of the type strain California kid (released in 2005), a draft genome of strain GM508D was determined to investigate genomic variation in this subspecies. Differences in mobile genetic element location and complement, catabolic pathway genes, contingency loci, surface antigen genes and type II restriction-modification systems highlight the plasticity and diversity within this taxon.

Populations adapt to fluctuating selection using derived and ancestral allelic diversity

Populations can adapt to changing environments by using allelic diversity, yet whether diversity is recently derived or ancestral is often debated. Although evolution could productively use both types of diversity in a changing environment, their relative frequency has not been quantified. We address this question experimentally using budding yeast strains that harbor a tandem repeat containing URA3 gene, which we expose to cyclical selection and counterselection. We characterize and quantify the dynamics of frameshift events in the URA3 gene in eight populations over 12 cycles of selection and find that ancestral alleles account for 10-20% of all adaptive events. Using a general model of fluctuating selection, we determine how these results depend on mutation rates, population sizes, and fluctuation timescales. We quantify the contribution of derived alleles to the adaptation process using the de novo mutation rate along the population's ancestral lineage, a novel measure that is applicable in a wide range of settings. We find that the adaptive dynamics undergoes a sharp transition from selection on ancestral alleles to selection on derived alleles as fluctuation timescales increase. Our results demonstrate that fluctuations can select between different modes of adaptation over evolutionary timescales.

Specificity of the ModA11, ModA12 and ModD1 epigenetic regulator N(6)-adenine DNA methyltransferases of Neisseria meningitidis

Phase variation (random ON/OFF switching) of gene expression is a common feature of host-adapted pathogenic bacteria. Phase variably expressed N(6)-adenine DNA methyltransferases (Mod) alter global methylation patterns resulting in changes in gene expression. These systems constitute phase variable regulons called phasevarions. Neisseria meningitidis phasevarions regulate genes including virulence factors and vaccine candidates, and alter phenotypes including antibiotic resistance. The target site recognized by these Type III N(6)-adenine DNA methyltransferases is not known. Single molecule, real-time (SMRT) methylome analysis was used to identify the recognition site for three key N. meningitidis methyltransferases: ModA11 (exemplified by M.NmeMC58I) (5'-CGY M6A: G-3'), ModA12 (exemplified by M.Nme77I, M.Nme18I and M.Nme579II) (5'-AC M6A: CC-3') and ModD1 (exemplified by M.Nme579I) (5'-CC M6A: GC-3'). Restriction inhibition assays and mutagenesis confirmed the SMRT methylome analysis. The ModA11 site is complex and atypical and is dependent on the type of pyrimidine at the central position, in combination with the bases flanking the core recognition sequence 5'-CGY M6A: G-3'. The observed efficiency of methylation in the modA11 strain (MC58) genome ranged from 4.6% at 5'-GCGC M6A: GG-3' sites, to 100% at 5'-ACGT M6A: GG-3' sites. Analysis of the distribution of modified sites in the respective genomes shows many cases of association with intergenic regions of genes with altered expression due to phasevarion switching.

Mutations in pneumococcal cpsE generated via in vitro serial passaging reveal a potential mechanism of reduced encapsulation utilized by a conjunctival isolate

The polysaccharide capsule of Streptococcus pneumoniae is required for nasopharyngeal colonization and for invasive disease in the lungs, blood, and meninges. In contrast, the vast majority of conjunctival isolates are acapsular. The first serotype-specific gene in the capsule operon, cpsE, encodes the initiating glycosyltransferase and is one of the few serotype-specific genes that can tolerate null mutations. This report characterizes a spontaneously arising TIGR4 mutant exhibiting a reduced capsule, caused by a 6-nucleotide duplication in cpsE which results in duplication of Ala and Ile at positions 45 and 46. This strain (AI45dup) possessed more exposed phosphorylcholine and was hypersusceptible to C3 complement deposition compared to the wild type. Accordingly, the mutant was significantly better at forming abiotic biofilms and binding epithelial cells in vitro but was avirulent in a sepsis model. In vitro serial passaging of the wild-type strain failed to reproduce the AI45dup mutation but instead led to a variety of mutants with reduced capsule harboring single nucleotide polymorphisms (SNPs) in cpsE. A single passage in the sepsis model after high-dose inoculation readily yielded revertants of AI45dup with restored wild-type capsule level, but the majority of SNP alleles of cpsE could not revert, suppress, or bypass. Analysis of cpsE in conjunctival isolates revealed a strain with a single missense mutation at amino acid position 377, which was responsible for reduced encapsulation. This study supports the hypothesis that spontaneous, nonreverting mutations in cpsE serve as a form of adaptive mutation by providing a selective advantage to S. pneumoniae in niches where expression of capsule is detrimental.

IMPORTANCE

While the capsule of Streptococcus pneumoniae is required for colonization and invasive disease, most conjunctival isolates are acapsular by virtue of deletion of the entire capsular operon. We show that spontaneous acapsular mutants isolated in vitro harbor mostly nonrevertible single nucleotide polymorphism (SNP) null mutations in cpsE, encoding the initiating glycosyltransferase. From a small collection of acapsular conjunctival isolates, we identified one strain with a complete capsular operon but containing a SNP in cpsE that we show is responsible for the acapsular phenotype. We propose that acapsular conjunctival isolates may arise initially from such nonreverting SNP null mutations in cpsE, which can be followed later by deletion of portions or all of the cps operon.

Neisseria gonorrhoeae MutS affects pilin antigenic variation through mismatch correction and not by pilE guanine quartet binding

Many pathogens use homologous recombination to vary surface antigens to avoid immune surveillance. Neisseria gonorrhoeae achieves this in part by changing the properties of its surface pili in a process called pilin antigenic variation (AV). Pilin AV occurs by high-frequency gene conversion reactions that transfer silent pilS sequences into the expressed pilE locus and requires the formation of an upstream guanine quartet (G4) DNA structure to initiate this process. The MutS and MutL proteins of the mismatch correction (MMC) system act to correct mismatches after replication and prevent homeologous (i.e., partially homologous) recombination, but MutS orthologs can also bind to G4 structures. A previous study showed that mutation of MutS resulted in a 3-fold increase in pilin AV, which could be due to the loss of MutS antirecombination properties or loss of G4 binding. We tested two site-directed separation-of-function MutS mutants that are both predicted to bind to G4s but are not able to perform MMC. Pilus phase variation assays and DNA sequence analysis of pilE variants produced in these mutants showed that all three mutS mutants and a mutL mutant had similar increased frequencies of pilin AV. Moreover, the mutS mutants all showed similar increased levels of pilin AV-dependent synthetic lethality. These results show that antirecombination by MMC is the reason for the effect that MutS has on pilin AV and is not due to pilE G4 binding by MutS.

IMPORTANCE

Neisseria gonorrhoeae continually changes its outer surface proteins to avoid recognition by the immune system. N. gonorrhoeae alters the antigenicity of the pilus by directed recombination between partially homologous pilin copies in a process that requires a guanine quartet (G4) structure. The MutS protein of the mismatch correction (MMC) system prevents recombination between partially homologous sequences and can also bind to G4s. We confirmed that loss of MMC increases the frequency of pilin antigenic variation and that two MutS mutants that are predicted to separate the two different functions of MutS inhibit pilin variation similarly to a complete-loss-of-function mutant, suggesting that interaction of MutS with the G4 structure is not a major factor in this process.

Courting disaster: How diversification rate affects fitness under risk

Life is full of risk. To deal with this uncertainty, many organisms have evolved bet-hedging strategies that spread risk through phenotypic diversification. These rates of diversification can vary by orders of magnitude in different species. Here we examine how key characteristics of risk and organismal ecology affect the fitness consequences of variation in diversification rate. We find that rapid diversification is strongly favored when the risk faced has a wide spatial extent, with a single disaster affecting a large fraction of the population. This advantage is especially great in small populations subject to frequent disaster. In contrast, when risk is correlated through time, slow diversification is favored because it allows adaptive tracking of disasters that tend to occur in series. Naturally evolved diversification mechanisms in diverse organisms facing a broad array of environmental risks largely support these results. The theory presented in this article provides a testable ecological hypothesis to explain the prevalence of slow stochastic switching among microbes and rapid, within-clutch diversification strategies among plants and animals.

Evolution in changing environments: modifiers of mutation, recombination, and migration

The production and maintenance of genetic and phenotypic diversity under temporally fluctuating selection and the signatures of environmental changes in the patterns of this variation have been important areas of focus in population genetics. On one hand, periods of constant selection pull the genetic makeup of populations toward local fitness optima. On the other, to cope with changes in the selection regime, populations may evolve mechanisms that create a diversity of genotypes. By tuning the rates at which variability is produced--such as the rates of recombination, mutation, or migration--populations may increase their long-term adaptability. Here we use theoretical models to gain insight into how the rates of these three evolutionary forces are shaped by fluctuating selection. We compare and contrast the evolution of recombination, mutation, and migration under similar patterns of environmental change and show that these three sources of phenotypic variation are surprisingly similar in their response to changing selection. We show that the shape, size, variance, and asymmetry of environmental fluctuation have different but predictable effects on evolutionary dynamics.

A random six-phase switch regulates pneumococcal virulence via global epigenetic changes

Streptococcus pneumoniae (the pneumococcus) is the world’s foremost bacterial pathogen in both morbidity and mortality. Switching between phenotypic forms (or ‘phases’) that favour asymptomatic carriage or invasive disease was first reported in 1933. Here, we show that the underlying mechanism for such phase variation consists of genetic rearrangements in a Type I restriction-modification system (SpnD39III). The rearrangements generate six alternative specificities with distinct methylation patterns, as defined by single-molecule, real-time (SMRT) methylomics. The SpnD39III variants have distinct gene expression profiles. We demonstrate distinct virulence in experimental infection and in vivo selection for switching between SpnD39III variants. SpnD39III is ubiquitous in pneumococci, indicating an essential role in its biology. Future studies must recognize the potential for switching between these heretofore undetectable, differentiated pneumococcal subpopulations in vitro and in vivo. Similar systems exist in other bacterial genera, indicating the potential for broad exploitation of epigenetic gene regulation.

Pneumococci can alternate between harmless and highly virulent forms. Here the authors show that such variation may be due to random rearrangements in a genetic locus encoding a restriction-modification system, resulting in epigenetic changes that affect expression of many genes.

The genetics of Neisseria species

Neisseria gonorrhoeae and Neisseria meningitidis are closely related organisms that cause the sexually transmitted infection gonorrhea and serious bacterial meningitis and septicemia, respectively. Both species possess multiple mechanisms to alter the expression of surface-exposed proteins through the processes of phase and antigenic variation. This potential for wide variability in surface-exposed structures allows the organisms to always have subpopulations of divergent antigenic types to avoid immune surveillance and to contribute to functional variation. Additionally, the Neisseria are naturally competent for DNA transformation, which is their main means of genetic exchange. Although bacteriophages and plasmids are present in this genus, they are not as effective as DNA transformation for horizontal genetic exchange. There are barriers to genetic transfer, such as restriction-modification systems and CRISPR loci, that limit particular types of exchange. These host-restricted pathogens illustrate the rich complexity of genetics that can help define the similarities and differences of closely related organisms.

ModM DNA methyltransferase methylome analysis reveals a potential role for Moraxella catarrhalis phasevarions in otitis media

Moraxella catarrhalis is a significant cause of otitis media and exacerbations of chronic obstructive pulmonary disease. Here, we characterize a phase-variable DNA methyltransferase (ModM), which contains 5'-CAAC-3' repeats in its open reading frame that mediate high-frequency mutation resulting in reversible on/off switching of ModM expression. Three modM alleles have been identified (modM1-3), with modM2 being the most commonly found allele. Using single-molecule, real-time (SMRT) genome sequencing and methylome analysis, we have determined that the ModM2 methylation target is 5'-GAR(m6)AC-3', and 100% of these sites are methylated in the genome of the M. catarrhalis 25239 ModM2 on strain. Proteomic analysis of ModM2 on and off variants revealed that ModM2 regulates expression of multiple genes that have potential roles in colonization, infection, and protection against host defenses. Investigation of the distribution of modM alleles in a panel of M. catarrhalis strains, isolated from the nasopharynx of healthy children or middle ear effusions from patients with otitis media, revealed a statistically significant association of modM3 with otitis media isolates. The modulation of gene expression via the ModM phase-variable regulon (phasevarion), and the significant association of the modM3 allele with otitis media, suggests a key role for ModM phasevarions in the pathogenesis of this organism.

Differential distributions of mononucleotide repeat sequences in 256 viral genomes and its potential implications

Mononucleotide repeats (MNRs) have been systematically investigated in the genomes of eukaryotic and prokaryotic organisms. However, detailed information on the distribution of MNRs in viral genomes is limited. In this study, we examined the distributions of MNRs in 256 fully sequenced virus genomes which showed extensive variations across viral genomes, and is significantly influenced by both genome size and CG content. Furthermore, the ratio of the observed to the expected number of MNRs (O/E ratio) appears to be influenced by both the host range and genome type of a particular virus. Additionally, the densities and frequencies of MNRs in genic regions are lower than in non-coding regions, suggesting that selective pressure acts on viral genomes. We also discuss the potential functional roles that these MNR loci could play in virus genomes. To our knowledge, this is the first analysis focusing on MNRs in viruses, and our study could have potential implications for a deeper understanding of virus genome stability and the co-evolution that occurs between a virus and its host.

Mathematical and live meningococcal models for simple sequence repeat dynamics - coherent predictions and observations

Evolvability by means of simple sequence repeat (SSR) instability is a feature under the constant influence of opposing selective pressures to expand and compress the repeat tract and is mechanistically influenced by factors that affect genetic instability. In addition to direct selection for protein expression and structural integrity, other factors that influence tract length evolution were studied. The genetic instability of SSRs that switch the expression of antibiotic resistance ON and OFF was modelled mathematically and monitored in a panel of live meningococcal strains. The mathematical model showed that the SSR length of a theoretical locus in an evolving population may be shaped by direct selection of expression status (ON or OFF), tract length dependent (α) and tract length independent factors (β). According to the model an increase in α drives the evolution towards shorter tracts. An increase in β drives the evolution towards a normal distribution of tract lengths given that an upper and a lower limit are set. Insertion and deletion biases were shown to skew allelic distributions in both directions. The meningococcal SSR model was tested in vivo by monitoring the frequency of spectinomycin resistance OFF→ON switching in a designed locus. The instability of a comprehensive panel of the homopolymeric SSRs, constituted of a range of 5-13 guanine nucleotides, was monitored in wildtype and mismatch repair deficient backgrounds. Both the repeat length itself and mismatch repair deficiency were shown to influence the genetic instability of the homopolymeric tracts. A possible insertion bias was observed in tracts ≤G10. Finally, an inverse correlation between the number of tract-encoded amino acids and growth in the presence of ON-selection illustrated a limitation to SSR expansion in an essential gene associated with the designed model locus and the protein function mediating antibiotic resistance.

Mutations in the control of virulence sensor gene from Streptococcus pyogenes after infection in mice lead to clonal bacterial variants with altered gene regulatory activity and virulence

The cluster of virulence sensor (CovS)/responder (CovR) two-component operon (CovRS) regulates ∼15% of the genes of the Group A Streptococcal pyogenes (GAS) genome. Bacterial clones containing inactivating mutations in the covS gene have been isolated from patients with virulent invasive diseases. We report herein an assessment of the nature and types of covS mutations that can occur in both virulent and nonvirulent GAS strains, and assess whether a nonvirulent GAS can attain enhanced virulence through this mechanism. A group of mice were infected with a globally-disseminated clonal M1T1 GAS (isolate 5448), containing wild-type (WT) CovRS (5448/CovR⁺S⁺), or less virulent engineered GAS strains, AP53/CovR⁺S⁺ and Manfredo M5/CovR⁺S⁺. SpeB negative GAS clones from wound sites and/or from bacteria disseminated to the spleen were isolated and the covS gene was subjected to DNA sequence analysis. Numerous examples of inactivating mutations were found in CovS in all regions of the gene. The mutations found included frame-shift insertions and deletions, and in-frame small and large deletions in the gene. Many of the mutations found resulted in early translation termination of CovS. Thus, the covS gene is a genomic mutagenic target that gives GAS enhanced virulence. In cases wherein CovS⁻ was discovered, these clonal variants exhibited high lethality, further suggesting that randomly mutated covS genes occur during the course of infection, and lead to the development of a more invasive infection.

Fitness is strongly influenced by rare mutations of large effect in a microbial mutation accumulation experiment

Our understanding of the evolutionary consequences of mutation relies heavily on estimates of the rate and fitness effect of spontaneous mutations generated by mutation accumulation (MA) experiments. We performed a classic MA experiment in which frequent sampling of MA lines was combined with whole genome resequencing to develop a high-resolution picture of the effect of spontaneous mutations in a hypermutator (ΔmutS) strain of the bacterium Pseudomonas aeruginosa. After ∼644 generations of mutation accumulation, MA lines had accumulated an average of 118 mutations, and we found that average fitness across all lines decayed linearly over time. Detailed analyses of the dynamics of fitness change in individual lines revealed that a large fraction of the total decay in fitness (42.3%) was attributable to the fixation of rare, highly deleterious mutations (comprising only 0.5% of fixed mutations). Furthermore, we found that at least 0.64% of mutations were beneficial and probably fixed due to positive selection. The majority of mutations that fixed (82.4%) were base substitutions and we failed to find any signatures of selection on nonsynonymous or intergenic mutations. Short indels made up a much smaller fraction of the mutations that were fixed (17.4%), but we found evidence of strong selection against indels that caused frameshift mutations in coding regions. These results help to quantify the amount of natural selection present in microbial MA experiments and demonstrate that changes in fitness are strongly influenced by rare mutations of large effect.

Predicted configurations of oligosaccharide extensions in the lipooligosaccharide of nontypeable Haemophilus influenzae isolates

Lipooligosaccharide configurations were predicted in nontypeable Haemophilus influenzae isolates based on the presence of seven oligosaccharide extension-initiating genes (or alleles). Predicted configurations with 2 to 3 oligosaccharide extensions were more prevalent among middle ear than throat strains. In addition, strains with these configurations averaged higher levels of serum resistance than strains with other configurations.

Predicting the virulence of MRSA from its genome sequence

Microbial virulence is a complex and often multifactorial phenotype, intricately linked to a pathogen’s evolutionary trajectory. Toxicity, the ability to destroy host cell membranes, and adhesion, the ability to adhere to human tissues, are the major virulence factors of many bacterial pathogens, including Staphylococcus aureus. Here, we assayed the toxicity and adhesiveness of 90 MRSA (methicillin resistant S. aureus) isolates and found that while there was remarkably little variation in adhesion, toxicity varied by over an order of magnitude between isolates, suggesting different evolutionary selection pressures acting on these two traits. We performed a genome-wide association study (GWAS) and identified a large number of loci, as well as a putative network of epistatically interacting loci, that significantly associated with toxicity. Despite this apparent complexity in toxicity regulation, a predictive model based on a set of significant single nucleotide polymorphisms (SNPs) and insertion and deletions events (indels) showed a high degree of accuracy in predicting an isolate’s toxicity solely from the genetic signature at these sites. Our results thus highlight the potential of using sequence data to determine clinically relevant parameters and have further implications for understanding the microbial virulence of this opportunistic pathogen.

Phase variation mediates reductions in expression of surface proteins during persistent meningococcal carriage

Asymptomatic and persistent colonization of the upper respiratory tract by Neisseria meningitidis occurs despite elicitation of adaptive immune responses against surface antigens. A putative mechanism for facilitating host persistence of this bacterial commensal and pathogen is alterations in expression of surface antigens by simple sequence repeat (SSR)-mediated phase variation. We investigated how often phase variation occurs during persistent carriage by analyzing the SSRs of eight loci in multiple isolates from 21 carriers representative of 1 to 6 months carriage. Alterations in repeat number were detected by a GeneScan analysis and occurred at 0.06 mutations/gene/month of carriage. The expression states were determined by Western blotting and two genes, fetA and nadA, exhibited trends toward low expression states. A critical finding from our unique examination of combinatorial expression states, "phasotypes," was for significant reductions in expression of multiple phase-variable surface proteins during persistent carriage of some strains. The immune responses in these carriers were examined by measuring variant-specific PorA IgG antibodies, capsular group Y IgG antibodies and serum bactericidal activity in concomitant serum samples. Persistent carriage was associated with high levels of specific IgG antibodies and serum bactericidal activity while recent strain acquisition correlated with a significant induction of antibodies. We conclude that phase-variable genes are driven into lower expression states during long-term persistent meningococcal carriage, in part due to continuous exposure to antibody-mediated selection, suggesting localized hypermutation has evolved to facilitate host persistence.

Bacterial genome instability

Bacterial genomes are remarkably stable from one generation to the next but are plastic on an evolutionary time scale, substantially shaped by horizontal gene transfer, genome rearrangement, and the activities of mobile DNA elements. This implies the existence of a delicate balance between the maintenance of genome stability and the tolerance of genome instability. In this review, we describe the specialized genetic elements and the endogenous processes that contribute to genome instability. We then discuss the consequences of genome instability at the physiological level, where cells have harnessed instability to mediate phase and antigenic variation, and at the evolutionary level, where horizontal gene transfer has played an important role. Indeed, this ability to share DNA sequences has played a major part in the evolution of life on Earth. The evolutionary plasticity of bacterial genomes, coupled with the vast numbers of bacteria on the planet, substantially limits our ability to control disease.

Selection for phase variation of LOS biosynthetic genes frequently occurs in progression of non-typeable Haemophilus influenzae infection from the nasopharynx to the middle ear of human patients

Surface structures in Haemophilus influenzae are subject to rapid ON/OFF switching of expression, a process termed phase variation. We analyse tetranucleotide repeats controlling phase variation in lipo-oligosaccharide (LOS) genes of H. influenzae in paired isolates from both the nasopharynx and middle ears of paediatric patients with chronic or recurrent otitis media. A change in expression of at least one of the seven phase variable LOS biosynthesis genes was seen in 12 of the 21 strain pairs. Several strains showed switching of expression in multiple LOS genes, consistent with a key role for phase variable LOS biosynthetic genes in human infection.

Genomic and global approaches to unravelling how hypermutable sequences influence bacterial pathogenesis

Rapid adaptation to fluctuations in the host milieu contributes to the host persistence and virulence of bacterial pathogens. Adaptation is frequently mediated by hypermutable sequences in bacterial pathogens. Early bacterial genomic studies identified the multiplicity and virulence-associated functions of these hypermutable sequences. Thus, simple sequence repeat tracts (SSRs) and site-specific recombination were found to control capsular type, lipopolysaccharide structure, pilin diversity and the expression of outer membrane proteins. We review how the population diversity inherent in the SSR-mediated mechanism of localised hypermutation is being unlocked by the investigation of whole genome sequences of disease isolates, analysis of clinical samples and use of model systems. A contrast is presented between the problematical nature of analysing simple sequence repeats in next generation sequencing data and in simpler, pragmatic PCR-based approaches. Specific examples are presented of the potential relevance of this localized hypermutation to meningococcal pathogenesis. This leads us to speculate on the future prospects for unravelling how hypermutable mechanisms may contribute to the transmission, spread and persistence of bacterial pathogens.

Engineering reduced evolutionary potential for synthetic biology

The field of synthetic biology seeks to engineer reliable and predictable behaviors in organisms from collections of standardized genetic parts. However, unlike other types of machines, genetically encoded biological systems are prone to changes in their designed sequences due to mutations in their DNA sequences after these devices are constructed and deployed. Thus, biological engineering efforts can be confounded by undesired evolution that rapidly breaks the functions of parts and systems, particularly when they are costly to the host cell to maintain. Here, we explain the fundamental properties that determine the evolvability of biological systems. Then, we use this framework to review current efforts to engineer the DNA sequences that encode synthetic biology devices and the genomes of their microbial hosts to reduce their ability to evolve and therefore increase their genetic reliability so that they maintain their intended functions over longer timescales.

Single cell stochastic regulation of pilus phase variation by an attenuation-like mechanism

The molecular triggers leading to virulence of a number of human-adapted commensal bacteria such as Streptococcus gallolyticus are largely unknown. This opportunistic pathogen is responsible for endocarditis in the elderly and associated with colorectal cancer. Colonization of damaged host tissues with exposed collagen, such as cardiac valves and pre-cancerous polyps, is mediated by appendages referred to as Pil1 pili. Populations of S. gallolyticus are heterogeneous with the majority of cells weakly piliated while a smaller fraction is hyper piliated. We provide genetic evidences that heterogeneous pil1 expression depends on a phase variation mechanism involving addition/deletion of GCAGA repeats that modifies the length of an upstream leader peptide. Synthesis of longer leader peptides potentiates the transcription of the pil1 genes through ribosome-induced destabilization of a premature stem-loop transcription terminator. This study describes, at the molecular level, a new regulatory mechanism combining phase variation in a leader peptide-encoding gene and transcription attenuation. This simple and robust mechanism controls a stochastic heterogeneous pilus expression, which is important for evading the host immune system while ensuring optimal tissue colonization.

Streptococcus gallolyticus (formely known as S. bovis biotype I) is an emerging cause of septicemia and endocarditis in the elderly. Intriguingly, epidemiological studies revealed a strong association, up to 65%, between endocarditis due to S. gallolyticus and colorectal malignancies. Whether S. gallolyticus infection is a cause or a consequence of colon cancer remains to be investigated. We previously showed that colonization of damaged cardiac valves with exposed collagen is mediated by the Pil1 pilus in S. gallolyticus. In the present work, we report that Pil1 is heterogeneously expressed at the single cell level, giving rise to two distinct bacterial subpopulations, a majority of weakly piliated cells and a minority of hyper-piliated cells. We have characterized, at the molecular level, a novel regulatory mechanism responsible for Pil1 heterogeneous expression combining phase variation in the leader peptide and transcriptional attenuation. Pili are highly immunogenic proteins proposed as vaccine candidate in pathogenic streptococci whose expression involves a fitness cost due to the selective pressure of host immune responses. Hence, this robust and simple system mitigates susceptibility to immune defenses while ensuring optimal colonization of host tissues.

A variable homopolymeric G-repeat defines small RNA-mediated posttranscriptional regulation of a chemotaxis receptor in Helicobacter pylori

Phase variation of hypermutable simple sequence repeats (SSRs) is a widespread and stochastic mechanism to generate phenotypic variation within a population and thereby contributes to host adaptation of bacterial pathogens. Although several examples of SSRs that affect transcription or coding potential have been reported, we now show that a SSR also impacts small RNA-mediated posttranscriptional regulation. Based on in vitro and in vivo analyses, we demonstrate that a variable homopolymeric G-repeat in the leader of the TlpB chemotaxis receptor mRNA of the human pathogen Helicobacter pylori is directly targeted by a small RNA (sRNA), RepG (Regulator of polymeric G-repeats). Whereas RepG sRNA is highly conserved, the tlpB G-repeat length varies among diverse H. pylori strains, resulting in strain-specific RepG-mediated tlpB regulation. Based on modification of the G-repeat length within one strain, we demonstrate that the G-repeat length determines posttranscriptional regulation and can mediate both repression and activation of tlpB through RepG. In vitro translation assays show that this regulation occurs at the translational level and that RepG influences tlpB translation dependent on the G-repeat length. In contrast to the digital ON-OFF switches through frame-shift mutations within coding sequences, such modulation of posttranscriptional regulation allows for a gradual control of gene expression. This connection to sRNA-mediated posttranscriptional regulation might also apply to other genes with SSRs, which could be targeting sites of cis- or trans-encoded sRNAs, and thereby could facilitate host adaptation through sRNA-mediated fine-tuning of virulence gene expression.

Simple sequence repeats together with mismatch repair deficiency can bias mutagenic pathways in Pseudomonas aeruginosa during chronic lung infection

Pseudomonas aeruginosa is an opportunistic pathogen that chronically infects the airways of cystic fibrosis (CF) patients and undergoes a process of genetic adaptation based on mutagenesis. We evaluated the role of mononucleotide G:C and A:T simple sequence repeats (SSRs) in this adaptive process. An in silico survey of the genome sequences of 7 P. aeruginosa strains showed that mononucleotide G:C SSRs but not A:T SSRs were greatly under-represented in coding regions, suggesting a strong counterselection process for G:C SSRs with lengths >5 bp but not for A:T SSRs. A meta-analysis of published whole genome sequence data for a P. aeruginosa strain from a CF patient with chronic airway infection showed that G:C SSRs but not A:T SSRs were frequently mutated during the infection process through the insertion or deletion of one or more SSR subunits. The mutation tendency of G:C SSRs was length-dependent and increased exponentially as a function of SSR length. When this strain naturally became a stable Mismatch Repair System (MRS)-deficient mutator, the degree of increase of G:C SSRs mutations (5-fold) was much higher than that of other types of mutation (2.2-fold or less). Sequence analysis of several mutated genes reported for two different collections, both containing mutator and non-mutator strains of P. aeruginosa from CF chronic infections, showed that the proportion of G:C SSR mutations was significantly higher in mutators than in non-mutators, whereas no such difference was observed for A:T SSR mutations. Our findings, taken together, provide genome-scale evidences that under a MRS-deficient background, long G:C SSRs are able to stochastically bias mutagenic pathways by making the genes in which they are harbored more prone to mutation. The combination of MRS deficiency and virulence-related genes that contain long G:C SSRs is therefore a matter of concern in P. aeruginosa CF chronic infection.

Natural selection promotes antigenic evolvability

The hypothesis that evolvability - the capacity to evolve by natural selection - is itself the object of natural selection is highly intriguing but remains controversial due in large part to a paucity of direct experimental evidence. The antigenic variation mechanisms of microbial pathogens provide an experimentally tractable system to test whether natural selection has favored mechanisms that increase evolvability. Many antigenic variation systems consist of paralogous unexpressed 'cassettes' that recombine into an expression site to rapidly alter the expressed protein. Importantly, the magnitude of antigenic change is a function of the genetic diversity among the unexpressed cassettes. Thus, evidence that selection favors among-cassette diversity is direct evidence that natural selection promotes antigenic evolvability. We used the Lyme disease bacterium, Borrelia burgdorferi, as a model to test the prediction that natural selection favors amino acid diversity among unexpressed vls cassettes and thereby promotes evolvability in a primary surface antigen, VlsE. The hypothesis that diversity among vls cassettes is favored by natural selection was supported in each B. burgdorferi strain analyzed using both classical (dN/dS ratios) and Bayesian population genetic analyses of genetic sequence data. This hypothesis was also supported by the conservation of highly mutable tandem-repeat structures across B. burgdorferi strains despite a near complete absence of sequence conservation. Diversification among vls cassettes due to natural selection and mutable repeat structures promotes long-term antigenic evolvability of VlsE. These findings provide a direct demonstration that molecular mechanisms that enhance evolvability of surface antigens are an evolutionary adaptation. The molecular evolutionary processes identified here can serve as a model for the evolution of antigenic evolvability in many pathogens which utilize similar strategies to establish chronic infections.

In Darwinian evolution, feedback from natural selection leads to biased mutations

Natural selection provides feedback through which information about the environment and its recurring challenges is captured, inherited, and accumulated within genomes in the form of variations that contribute to survival. The variation upon which natural selection acts is generally described as "random." Yet evidence has been mounting for decades, from such phenomena as mutation hotspots, horizontal gene transfer, and highly mutable repetitive sequences, that variation is far from the simplifying idealization of random processes as white (uniform in space and time and independent of the environment or context).  This paper focuses on what is known about the generation and control of mutational variation, emphasizing that it is not uniform across the genome or in time, not unstructured with respect to survival, and is neither memoryless nor independent of the (also far from white) environment. We suggest that, as opposed to frequentist methods, Bayesian analysis could capture the evolution of nonuniform probabilities of distinct classes of mutation, and argue not only that the locations, styles, and timing of real mutations are not correctly modeled as generated by a white noise random process, but that such a process would be inconsistent with evolutionary theory.

The role of variable DNA tandem repeats in bacterial adaptation

DNA tandem repeats (TRs), also designated as satellite DNA, are inter- or intragenic nucleotide sequences that are repeated two or more times in a head-to-tail manner. Because TR tracts are prone to strand-slippage replication and recombination events that cause the TR copy number to increase or decrease, loci containing TRs are hypermutable. An increasing number of examples illustrate that bacteria can exploit this instability of TRs to reversibly shut down or modulate the function of specific genes, allowing them to adapt to changing environments on short evolutionary time scales without an increased overall mutation rate. In this review, we discuss the prevalence and distribution of inter- and intragenic TRs in bacteria and the mechanisms of their instability. In addition, we review evidence demonstrating a role of TR variations in bacterial adaptation strategies, ranging from immune evasion and tissue tropism to the modulation of environmental stress tolerance. Nevertheless, while bioinformatic analysis reveals that most bacterial genomes contain a few up to several dozens of intra- and intergenic TRs, only a small fraction of these have been functionally studied to date.

Natural variation in the promoter of the gene encoding the Mga regulator alters host-pathogen interactions in group a Streptococcus carrier strains

Humans commonly carry pathogenic bacteria asymptomatically, but the molecular factors underlying microbial asymptomatic carriage are poorly understood. We previously reported that two epidemiologically unassociated serotype M3 group A Streptococcus (GAS) carrier strains had an identical 12-bp deletion in the promoter of the gene encoding Mga, a global positive gene regulator. Herein, we report on studies designed to test the hypothesis that the identified 12-bp deletion in the mga promoter alters GAS virulence, thereby potentially contributing to the asymptomatic carrier phenotype. Using allelic exchange, we introduced the variant promoter into a serotype M3 invasive strain and the wild-type promoter into an asymptomatic carrier strain. Compared to strains with the wild-type mga promoter, we discovered that strains containing the promoter with the 12-bp deletion produced significantly fewer mga and Mga-regulated gene transcripts. Consistent with decreased mga transcripts, strains containing the variant mga promoter were also significantly less virulent in in vivo and ex vivo models of GAS disease. Further, we provide evidence that the pleiotropic regulator protein CodY binds to the mga promoter and that the 12-bp deletion in the mga promoter reduces CodY-mediated mga transcription. We conclude that the naturally occurring 12-bp deletion in the mga promoter significantly alters the pathogen-host interaction of these asymptomatic carrier strains. Our findings provide new insight into the molecular basis of the carrier state of an important human pathogen.

Adaptation of Enterococcus faecalis to daptomycin reveals an ordered progression to resistance

With increasing numbers of hospital-acquired antibiotic resistant infections each year and staggering health care costs, there is a clear need for new antimicrobial agents, as well as novel strategies to extend their clinical efficacy. While genomic studies have provided a wealth of information about the alleles associated with adaptation to antibiotics, they do not provide essential information about the relative importance of genomic changes, their order of appearance, or potential epistatic relationships between adaptive changes. Here we used quantitative experimental evolution of a single polymorphic population in continuous culture with whole-genome sequencing and allelic frequency measurements to study daptomycin (DAP) resistance in the vancomycin-resistant clinical pathogen Enterococcus faecalis S613. Importantly, we sustained both planktonic and nonplanktonic (i.e., biofilm) populations in coculture as the concentration of antibiotic was raised, facilitating the development of more ecological complexity than is typically observed in laboratory evolution. Quantitative experimental evolution revealed a clear order and hierarchy of genetic changes leading to resistance, the signaling and metabolic pathways responsible, and the relative importance of these mutations to the evolution of DAP resistance. Despite the relative simplicity of this ex vivo approach compared to the ecological complexity of the human body, we showed that experimental evolution allows for rapid identification of clinically relevant adaptive molecular pathways and new targets for drug design in pathogens.

Prevalence and phase variable expression status of two autotransporters, NalP and MspA, in carriage and disease isolates of Neisseria meningitidis

Neisseria meningitidis is a human nasopharyngeal commensal capable of causing life-threatening septicemia and meningitis. Many meningococcal surface structures, including the autotransporter proteins NalP and MspA, are subject to phase variation (PV) due to the presence of homopolymeric tracts within their coding sequences. The functions of MspA are unknown. NalP proteolytically cleaves several surface-located virulence factors including the 4CMenB antigen NhbA. Therefore, NalP is a phase-variable regulator of the meningococcal outer membrane and secretome whose expression may reduce isolate susceptibility to 4CMenB-induced immune responses. To improve our understanding of the contributions of MspA and NalP to meningococcal-host interactions, their distribution and phase-variable expression status was studied in epidemiologically relevant samples, including 127 carriage and 514 invasive isolates representative of multiple clonal complexes and serogroups. Prevalence estimates of >98% and >88% were obtained for mspA and nalP, respectively, with no significant differences in their frequencies in disease versus carriage isolates. 16% of serogroup B (MenB) invasive isolates, predominately from clonal complexes ST-269 and ST-461, lacked nalP. Deletion of nalP often resulted from recombination events between flanking repetitive elements. PolyC tract lengths ranged from 6–15 bp in nalP and 6–14 bp in mspA. In an examination of PV status, 58.8% of carriage, and 40.1% of invasive nalP-positive MenB isolates were nalP phase ON. The frequency of this phenotype was not significantly different in serogroup Y (MenY) carriage strains, but was significantly higher in invasive MenY strains (86.3%; p<0.0001). Approximately 90% of MenB carriage and invasive isolates were mspA phase ON; significantly more than MenY carriage (32.7%) or invasive (13.7%) isolates. This differential expression resulted from different mode mspA tract lengths between the serogroups. Our data indicates a differential requirement for NalP and MspA expression in MenB and MenY strains and is a step towards understanding the contributions of phase-variable loci to meningococcal biology.

Evolution in microbes

This review presents a broad survey of experimental microbial evolution, covering diverse topics including trade-offs, epistasis, fluctuating conditions, spatial dynamics, cooperation, aging, and stochastic switching. Emphasis is placed on examples that highlight key conceptual points or address theoretical predictions. Experimental evolution is discussed from two points of view. First, population trajectories are described as adaptive walks on a fitness landscape, whose genetic structure can be probed by experiments. Second, populations are viewed from a physiological perspective, and their nongenetic heterogeneity is examined. Bringing together these two viewpoints remains a major challenge for the future.

Genome-wide association study identifies vitamin B5 biosynthesis as a host specificity factor in Campylobacter

Genome-wide association studies have the potential to identify causal genetic factors underlying important phenotypes but have rarely been performed in bacteria. We present an association mapping method that takes into account the clonal population structure of bacteria and is applicable to both core and accessory genome variation. Campylobacter is a common cause of human gastroenteritis as a consequence of its proliferation in multiple farm animal species and its transmission via contaminated meat and poultry. We applied our association mapping method to identify the factors responsible for adaptation to cattle and chickens among 192 Campylobacter isolates from these and other host sources. Phylogenetic analysis implied frequent host switching but also showed that some lineages were strongly associated with particular hosts. A seven-gene region with a host association signal was found. Genes in this region were almost universally present in cattle but were frequently absent in isolates from chickens and wild birds. Three of the seven genes encoded vitamin B5 biosynthesis. We found that isolates from cattle were better able to grow in vitamin B5-depleted media and propose that this difference may be an adaptation to host diet.

Complete genome sequence and transcription profiles of the rock bream iridovirus RBIV-C1

The family Iridoviridae consists of 5 genera of double-stranded DNA viruses, including the genus Megalocytivirus, which contains species that are important fish pathogens. In a previous study, we isolated the first rock bream iridovirus from China (RBIV-C1) and identified it as a member of the genus Megalocytivirus. In this report, we determined the complete genomic sequence of RBIV-C1 and examined its in vivo expression profiles. The genome of RBIV-C1 is 112333 bp in length, with a GC content of 55% and a coding density of 92%. RBIV-C1 contains 4584 simple sequence repeats, 89.8% of which are distributed among coding regions. A total of 119 potential open reading frames (ORFs) were identified in RBIV-C1, including the 26 core iridovirus genes; 41 ORFs encode proteins that are predicted to be associated with essential biological functions. RBIV-C1 exhibits the highest degree of sequence conservation and colinear arrangement of genes with orange-spotted grouper iridovirus (OSGIV) and rock bream iridovirus (RBIV). The pairwise nucleotide identities are 99.49% between RBIV-C1 and OSGIV and 98.69% between RBIV-C1 and RBIV. Compared to OSGIV, RBIV-C1 contains 11 insertions, 13 deletions, and 103 single nucleotide mutations. Whole-genome transcription analysis showed that following experimental infection of rock bream with RBIV-C1, all but 1 of the 119 ORFs were expressed at different time points and clustered into 3 hierarchical groups based on their expression patterns. These results provide new insights into the genetic nature and gene expression features of megalocytiviruses.

Selection-driven accumulation of suppressor mutants in bacillus subtilis: the apparent high mutation frequency of the cryptic gudB gene and the rapid clonal expansion of gudB(+) suppressors are due to growth under selection

Soil bacteria like Bacillus subtilis can cope with many growth conditions by adjusting gene expression and metabolic pathways. Alternatively, bacteria can spontaneously accumulate beneficial mutations or shape their genomes in response to stress. Recently, it has been observed that a B. subtilis mutant lacking the catabolically active glutamate dehydrogenase (GDH), RocG, mutates the cryptic gudB^CR gene at a high frequency. The suppressor mutants express the active GDH GudB, which can fully replace the function of RocG. Interestingly, the cryptic gudB^CR allele is stably inherited as long as the bacteria synthesize the functional GDH RocG. Competition experiments revealed that the presence of the cryptic gudB^CR allele provides the bacteria with a selective growth advantage when glutamate is scarce. Moreover, the lack of exogenous glutamate is the driving force for the selection of mutants that have inactivated the active gudB gene. In contrast, two functional GDHs are beneficial for the cells when glutamate was available. Thus, the amount of GDH activity strongly affects fitness of the bacteria depending on the availability of exogenous glutamate. At a first glance the high mutation frequency of the cryptic gudB^CR allele might be attributed to stress-induced adaptive mutagenesis. However, other loci on the chromosome that could be potentially mutated during growth under the selective pressure that is exerted on a GDH-deficient mutant remained unaffected. Moreover, we show that a GDH-proficient B. subtilis strain has a strong selective growth advantage in a glutamate-dependent manner. Thus, the emergence and rapid clonal expansion of the active gudB allele can be in fact explained by spontaneous mutation and growth under selection without an increase of the mutation rate. Moreover, this study shows that the selective pressure that is exerted on a maladapted bacterium strongly affects the apparent mutation frequency of mutational hot spots.

Allelic variation in a simple sequence repeat element of neisserial pglB2 and its consequences for protein expression and protein glycosylation

Neisseria species express an O-linked glycosylation system in which functionally distinct proteins are elaborated with variable glycans. A major source of glycan diversity in N. meningitidis results from two distinct pglB alleles responsible for the synthesis of either N,N'-diacetylbacillosamine or glyceramido-acetamido trideoxyhexose that occupy the reducing end of the oligosaccharides. Alternative modifications at C-4 of the precursor UDP-4-amino are attributable to distinct C-terminal domains that dictate either acetyltransferase or glyceramidotransferase activity, encoded by pglB and pglB2, respectively. Naturally occurring alleles of pglB2 have homopolymeric tracts of either 7 or 8 adenosines (As) bridging the C-terminal open reading frame (ORF) and the ORF encompassing the conserved N-terminal domain associated with phosphoglycosyltransferase activity. In the work presented here, we explored the consequences of such pglB2 allele variation and found that, although both alleles are functional vis-à-vis glycosylation, the 7A form results in the expression of a single, multidomain protein, while the 8A variant elicits two single-domain proteins. We also found that the glyceramidotransferase activity-encoding domain is essential to protein glycosylation, showing the critical role of the C-4 modification of the precursor UDP-4-amino in the pathway. These findings were further extended and confirmed by examining the phenotypic consequences of extended poly(A) tract length variation. Although ORFs related to those of pglB2 are broadly distributed in eubacteria, they are primarily found as two distinct, juxtaposed ORFs. Thus, the neisserial pglB2 system provides novel insights into the potential influence of hypermutability on modular evolution of proteins by providing a unique snapshot of the progression of ongoing gene fusion.

Reverse vaccinology in the 21st century: improvements over the original design

Reverse vaccinology (RV), the first application of genomic technologies in vaccine research, represented a major revolution in the process of discovering novel vaccines. By determining their entire antigenic repertoire, researchers could identify protective targets and design efficacious vaccines for pathogens where conventional approaches had failed. Bexsero, the first vaccine developed using RV, has recently received positive opinion from the European Medicines Agency. The use of RV initiated a cascade of changes that affected the entire vaccine development process, shifting the focus from the identification of a list of vaccine candidates to the definition of a set of high throughput screens to reduce the need for costly and labor intensive tests in animal models. It is now clear that a deep understanding of the epidemiology of vaccine candidates, and their regulation and role in host-pathogen interactions, must become an integral component of the screening workflow. Far from being outdated by technological advancements, RV still represents a paradigm of how high-throughput technologies and scientific insight can be integrated into biotechnology research.