Synonymous Genetic Variation in Natural Isolates of Escherichia coli Does Not Predict Where Synonymous Substitutions Occur in a Long-Term Experiment

Core genes can have higher recombination rates than accessory genes within global microbial populations

Recombination is essential to microbial evolution, and is involved in the spread of antibiotic resistance, antigenic variation, and adaptation to the host niche. However, assessing the impact of homologous recombination on accessory genes which are only present in a subset of strains of a given species remains challenging due to their complex phylogenetic relationships. Quantifying homologous recombination for accessory genes (which are important for niche-specific adaptations) in comparison to core genes (which are present in all strains and have essential functions) is critical to understanding how selection acts on variation to shape species diversity and genome structures of bacteria. Here, we apply a computationally efficient, non-phylogenetic approach to measure homologous recombination rates in the core and accessory genome using >100,000 whole genome sequences from Streptococcus pneumoniae and several additional species. By analyzing diverse sets of sequence clusters, we show that core genes often have higher recombination rates than accessory genes, and for some bacterial species the associated effect sizes for these differences are pronounced. In a subset of species, we find that gene frequency and homologous recombination rate are positively correlated. For S. pneumoniae and several additional species, we find that while the recombination rate is higher for the core genome, the mutational divergence is lower, indicating that divergence-based homologous recombination barriers could contribute to differences in recombination rates between the core and accessory genome. Homologous recombination may therefore play a key role in increasing the efficiency of selection in the most conserved parts of the genome.

Divergent Evolution of Mutation Rates and Biases in the Long-Term Evolution Experiment with Escherichia coli

All organisms encode enzymes that replicate, maintain, pack, recombine, and repair their genetic material. For this reason, mutation rates and biases also evolve by mutation, variation, and natural selection. By examining metagenomic time series of the Lenski long-term evolution experiment (LTEE) with Escherichia coli (Good BH, McDonald MJ, Barrick JE, Lenski RE, Desai MM. 2017. The dynamics of molecular evolution over 60,000 generations. Nature 551(7678):45–50.), we find that local mutation rate variation has evolved during the LTEE. Each LTEE population has evolved idiosyncratic differences in their rates of point mutations, indels, and mobile element insertions, due to the fixation of various hypermutator and antimutator alleles. One LTEE population, called Ara+3, shows a strong, symmetric wave pattern in its density of point mutations, radiating from the origin of replication. This pattern is largely missing from the other LTEE populations, most of which evolved missense, indel, or structural mutations in topA, fis, and dusB—loci that all affect DNA topology. The distribution of mutations in those genes over time suggests epistasis and historical contingency in the evolution of DNA topology, which may have in turn affected local mutation rates. Overall, the replicate populations of the LTEE have largely diverged in their mutation rates and biases, even though they have adapted to identical abiotic conditions.

Universal Constraints on Protein Evolution in the Long-Term Evolution Experiment with Escherichia coli

Although it is well known that abundant proteins evolve slowly across the tree of life, there is little consensus for why this is true. Here, I report that abundant proteins evolve slowly in the hypermutator populations of Lenski's long-term evolution experiment with Escherichia coli (LTEE). Specifically, the density of all observed mutations per gene, as measured in metagenomic time series covering 60,000 generations of the LTEE, significantly anticorrelates with mRNA abundance, protein abundance, and degree of protein-protein interaction. The same pattern holds for nonsynonymous mutation density. However, synonymous mutation density, measured across the LTEE hypermutator populations, positively correlates with protein abundance. These results show that universal constraints on protein evolution are visible in data spanning three decades of experimental evolution. Therefore, it should be possible to design experiments to answer why abundant proteins evolve slowly.

Comparative genomics of Staphylococcus aureus associated with subclinical and clinical bovine mastitis

Many efforts have been made to understand the pathogenesis of bovine mastitis to reduce losses and promote animal welfare. Staphylococcus aureus may cause bovine clinical mastitis, but it is mainly associated with subclinical infection, which is usually persistent and can easily reoccur. Here, we conducted a comparative genomic analysis between strains of S. aureus causing subclinical infection (Sau170, 302, 1269, 1364), previously sequenced by our group, and two well-characterized strains causing clinical mastitis (N305 and RF122) to find differences that could be linked to mastitis outcome. A total of 146 virulence-associated genes were compared and no appreciable differences were found between the bacteria. However, several nonsynonymous single nucleotide polymorphisms (SNPs) were identified in genes present in the subclinical strains when compared to RF122 and N305, especially in genes encoding host immune evasion and surface proteins. The secreted and surface proteins predicted by in silico tools were compared through multidimensional scaling analysis (MDS), revealing a high degree of similarity among the strains. The comparison of orthologous genes by OrthoMCL identified a membrane transporter and a lipoprotein as exclusive of bacteria belonging to the subclinical and clinical groups, respectively. No hit was found in RF122 and N305 for the membrane transporter using BLAST algorithm. For the lipoprotein, sequences of Sau170, 302, 1269, and 1364 with identities between 68–73% were found in the MDS dataset. A conserved region found only in the lipoprotein genes of RF122 and N305 was used for primer design. Although the polymerase chain reaction (PCR) on field isolates of S. aureus did not validate the findings for the transporter, the lipoprotein was able to separate the clinical from the subclinical isolates. These results show that sequence variation among bovine S. aureus, and not only the presence/absence of virulence factors, is an important aspect to consider when comparing isolates causing different mastitis outcomes.

The comparative population genetics of Neisseria meningitidis and Neisseria gonorrhoeae

Neisseria meningitidis and N. gonorrhoeae are closely related pathogenic bacteria. To compare their population genetics, we compiled a dataset of 1,145 genes found across 20 N. meningitidis and 15 N. gonorrhoeae genomes. We find that N. meningitidis is seven-times more diverse than N. gonorrhoeae in their combined core genome. Both species have acquired the majority of their diversity by recombination with divergent strains, however, we find that N. meningitidis has acquired more of its diversity by recombination than N. gonorrhoeae. We find that linkage disequilibrium (LD) declines rapidly across the genomes of both species. Several observations suggest that N. meningitidis has a higher effective population size than N. gonorrhoeae; it is more diverse, the ratio of non-synonymous to synonymous polymorphism is lower, and LD declines more rapidly to a lower asymptote in N. meningitidis. The two species share a modest amount of variation, half of which seems to have been acquired by lateral gene transfer and half from their common ancestor. We investigate whether diversity varies across the genome of each species and find that it does. Much of this variation is due to different levels of lateral gene transfer. However, we also find some evidence that the effective population size varies across the genome. We test for adaptive evolution in the core genome using a McDonald–Kreitman test and by considering the diversity around non-synonymous sites that are fixed for different alleles in the two species. We find some evidence for adaptive evolution using both approaches.

Genetic Diversity in FUB Genes of Fusarium oxysporum f. sp. cubense Suggests Horizontal Gene Transfer

Fusaric acid (FA) is an important secondary metabolite of many Fusarium species and involved in the wilt symptoms caused in banana by Fusarium oxysporum f. sp. cubense (Foc). To investigate the evolution characteristics of the 12 Foc FA biosynthetic genes (FUB), coding sequences of the 12 FUB genes and three housekeeping genes, EF-1α/RPB1/RPB2 (translation elongation factor-1α/RNA polymerase II subunit I/RNA polymerase II subunit II), were subjected to genetic diversity analysis, phylogenetic analysis, recombination detection, and selective pressure analysis. The results of selective pressure analysis showed that the 15 genes were mainly subjected to negative selection. However, a significantly higher number of silent mutations, which could not be simply explained by selective pressure difference, were observed in the 12 FUB genes in Foc than in the three housekeeping genes. Infraspecies phylogeny and recombination detection analysis showed that significantly more horizontal gene transfer (HGT) events (normalized) had occurred in the FUB genes than in the three housekeeping genes. In addition, many of these events involved outgroup isolates and significantly increased the genetic diversity of FUB genes in Foc. The infraspecies phylogenetic analysis suggested that the polyphyletic phylogeny proposed for Foc requires further discussion, and the divergence of race 1, race 4, and the common ancestor of several F. oxysporum (Fo) isolates pathogenic to nonbanana plants should have diverged over a short period. Finally, our results suggest that the FUB genes in Fo should have benefited from HGT to gain a relatively high genetic diversity to respond to different host plants and environments despite mainly being subject to negative selection.

Identifying Drivers of Parallel Evolution: A Regression Model Approach

Parallel evolution, defined as identical changes arising in independent populations, is often attributed to similar selective pressures favoring the fixation of identical genetic changes. However, some level of parallel evolution is also expected if mutation rates are heterogeneous across regions of the genome. Theory suggests that mutation and selection can have equal impacts on patterns of parallel evolution; however, empirical studies have yet to jointly quantify the importance of these two processes. Here, we introduce several statistical models to examine the contributions of mutation and selection heterogeneity to shaping parallel evolutionary changes at the gene-level. Using this framework, we analyze published data from forty experimentally evolved Saccharomyces cerevisiae populations. We can partition the effects of a number of genomic variables into those affecting patterns of parallel evolution via effects on the rate of arising mutations, and those affecting the retention versus loss of the arising mutations (i.e., selection). Our results suggest that gene-to-gene heterogeneity in both mutation and selection, associated with gene length, recombination rate, and number of protein domains drive parallel evolution at both synonymous and nonsynonymous sites. While there are still a number of parallel changes that are not well described, we show that allowing for heterogeneous rates of mutation and selection can provide improved predictions of the prevalence and degree of parallel evolution.

Parallels between experimental and natural evolution of legume symbionts

The emergence of symbiotic interactions has been studied using population genomics in nature and experimental evolution in the laboratory, but the parallels between these processes remain unknown. Here we compare the emergence of rhizobia after the horizontal transfer of a symbiotic plasmid in natural populations of Cupriavidus taiwanensis, over 10 MY ago, with the experimental evolution of symbiotic Ralstonia solanacearum for a few hundred generations. In spite of major differences in terms of time span, environment, genetic background, and phenotypic achievement, both processes resulted in rapid genetic diversification dominated by purifying selection. We observe no adaptation in the plasmid carrying the genes responsible for the ecological transition. Instead, adaptation was associated with positive selection in a set of genes that led to the co-option of the same quorum-sensing system in both processes. Our results provide evidence for similarities in experimental and natural evolutionary transitions and highlight the potential of comparisons between both processes to understand symbiogenesis.

It is unclear if experimental evolution is a good model for natural processes. Here, Clerissi et al. find parallels between the evolution of symbiosis in rhizobia after horizontal transfer of a plasmid over 10 million years ago and experimentally evolved symbionts.

Analysis of bacterial genomes from an evolution experiment with horizontal gene transfer shows that recombination can sometimes overwhelm selection

Few experimental studies have examined the role that sexual recombination plays in bacterial evolution, including the effects of horizontal gene transfer on genome structure. To address this limitation, we analyzed genomes from an experiment in which Escherichia coli K-12 Hfr (high frequency recombination) donors were periodically introduced into 12 evolving populations of E. coli B and allowed to conjugate repeatedly over the course of 1000 generations. Previous analyses of the evolved strains from this experiment showed that recombination did not accelerate adaptation, despite increasing genetic variation relative to asexual controls. However, the resolution in that previous work was limited to only a few genetic markers. We sought to clarify and understand these puzzling results by sequencing complete genomes from each population. The effects of recombination were highly variable: one lineage was mostly derived from the donors, while another acquired almost no donor DNA. In most lineages, some regions showed repeated introgression and others almost none. Regions with high introgression tended to be near the donors' origin of transfer sites. To determine whether introgressed alleles imposed a genetic load, we extended the experiment for 200 generations without recombination and sequenced whole-population samples. Beneficial alleles in the recipient populations were occasionally driven extinct by maladaptive donor-derived alleles. On balance, our analyses indicate that the plasmid-mediated recombination was sufficiently frequent to drive donor alleles to fixation without providing much, if any, selective advantage.

The adaptation of Escherichia coli cells grown in simulated microgravity for an extended period is both phenotypic and genomic

Microorganisms impact spaceflight in a variety of ways. They play a positive role in biological systems, such as waste water treatment but can be problematic through buildups of biofilms that can affect advanced life support. Of special concern is the possibility that during extended missions, the microgravity environment will provide positive selection for undesirable genomic changes. Such changes could affect microbial antibiotic sensitivity and possibly pathogenicity. To evaluate this possibility, Escherichia coli (lac plus) cells were grown for over 1000 generations on Luria Broth medium under low-shear modeled microgravity conditions in a high aspect rotating vessel. This is the first study of its kind to grow bacteria for multiple generations over an extended period under low-shear modeled microgravity. Comparisons were made to a non-adaptive control strain using growth competitions. After 1000 generations, the final low-shear modeled microgravity-adapted strain readily outcompeted the unadapted lac minus strain. A portion of this advantage was maintained when the low-shear modeled microgravity strain was first grown in a shake flask environment for 10, 20, or 30 generations of growth. Genomic sequencing of the 1000 generation strain revealed 16 mutations. Of the five changes affecting codons, none were neutral. It is not clear how significant these mutations are as individual changes or as a group. It is concluded that part of the long-term adaptation to low-shear modeled microgravity is likely genomic. The strain was monitored for acquisition of antibiotic resistance by VITEK analysis throughout the adaptation period. Despite the evidence of genomic adaptation, resistance to a variety of antibiotics was never observed.

Core Genes Evolve Rapidly in the Long-Term Evolution Experiment with Escherichia coli

Bacteria can evolve rapidly under positive selection owing to their vast numbers, allowing their genes to diversify by adapting to different environments. We asked whether the same genes that evolve rapidly in the long-term evolution experiment (LTEE) with Escherichia coli have also diversified extensively in nature. To make this comparison, we identified ∼2000 core genes shared among 60 E. coli strains. During the LTEE, core genes accumulated significantly more nonsynonymous mutations than flexible (i.e., noncore) genes. Furthermore, core genes under positive selection in the LTEE are more conserved in nature than the average core gene. In some cases, adaptive mutations appear to modify protein functions, rather than merely knocking them out. The LTEE conditions are novel for E. coli, at least in relation to its evolutionary history in nature. The constancy and simplicity of the environment likely favor the complete loss of some unused functions and the fine-tuning of others.

Enigmas in tumor resistance to kinase inhibitors and calculation of the drug resistance index for cancer (DRIC)

Darwinian selection is also applicable when antibiotics, the immune system or other host factors shape the repertoire of microorganisms, and similarly, clonal selection is the hallmark of tumor evolution. The ongoing revolution in new anti-cancer treatment modalities, combined with an unprecedented precision in characterizing malignant clones at the level below one percent, profoundly improves the understanding of repertoire-tuning mechanisms. There is no fundamental difference between selection of the tumor cells in the presence, or absence, of therapy. However, under treatment the influence of a single agent can be measured, simplifying the analysis. Because of their beneficial and selective therapeutic effect, the focus in this review is set on protein kinase inhibitors (PKIs), predominantly tyrosine kinase inhibitors (TKIs). This is one of the most rapidly growing families of novel cancer medicines. In order to limit the number of drugs, the following representative target kinases are included: ALK, BCR-ABL, BRAF, BTK, and EGFR. A key therapeutic challenge is how to reduce tumor growth after treatment, since this is rate-limiting for the generation and expansion of more malignant escape mutants. Thus, upon efficient treatment, tumor cell loss often enables a profoundly increased growth rate among resistant cells. Strategies to reduce this risk, such as concomitant, competitive outgrowth of non-transformed cells, are described. Seven parameters: 1. Drug type, 2. tumor type, 3. presence of metastases or phenotypic change, 4. tumor cell number, 5. net growth rate (proliferation minus cell death), 6. inherited genetic- and 7. epigenetic- variations are crucial for drug responses. It is envisaged that it might become possible to calculate a clinically relevant Drug Resistance Index for Cancer (DRIC) for each patient.

Tempo and mode of genome evolution in a 50,000-generation experiment

Adaptation by natural selection depends on the rates, effects, and interactions of many mutations, making it difficult to determine what proportion of mutations in an evolving lineage are beneficial. We analysed 264 complete genomes from 12 Escherichia coli populations to characterize their dynamics over 50,000 generations. The populations that retained the ancestral mutation rate support a model where most fixed mutations are beneficial, the fraction of beneficial mutations declines as fitness rises, and neutral mutations accumulate at a constant rate. We also compared these populations to mutation-accumulation lines evolved under a bottlenecking regime that minimizes selection. Nonsynonymous mutations, intergenic mutations, insertions, and deletions are overrepresented in the long-term populations, further supporting the inference that most mutations that reached high frequency were favoured by selection. These results illuminate the shifting balance of forces that govern genome evolution in populations adapting to a new environment.

Gene flow in microbial communities could explain unexpected patterns of synonymous variation in the Escherichia coli core genome

Researchers contest the importance of gene flow in bacterial core genomes, as traditionalists view microbes as predominantly clonal, asexually reproducing organisms. Contrary to the traditional perspective, Escherichia coli core genes vary greatly in their levels of synonymous genetic diversity. This observation indicates that the relative importance of evolutionary forces such as mutation, selection, and recombination varies from gene to gene. In this paper, I highlight why the synonymous diversity observation is broadly relevant to researchers interested in the evolutionary dynamics of microbial populations and communities. I explain how a model of evolution called the coalescent relates neutral diversity (i.e. mutations with negligible fitness effects) to mutation rates, evolutionary time, and a parameter called effective population size. I then describe the possible ways in which mutation, selection, and recombination can explain observed patterns of synonymous diversity in E. coli. Finally, I describe a model for E. coli genome evolution in which different loci are subject to varying levels of gene flow among co-occurring microbes and viruses in the environment. Researchers can falsify the gene flow hypothesis by sequencing genes and strains isolated from stable microbiomes or by carrying out evolution experiments that trace gene genealogies in real-time.

Can the experimental evolution programme help us elucidate the genetic basis of adaptation in nature?

There have been a variety of approaches taken to try to characterize and identify the genetic basis of adaptation in nature, spanning theoretical models, experimental evolution studies and direct tests of natural populations. Theoretical models can provide formalized and detailed hypotheses regarding evolutionary processes and patterns, from which experimental evolution studies can then provide important proofs of concepts and characterize what is biologically reasonable. Genetic and genomic data from natural populations then allow for the identification of the particular factors that have and continue to play an important role in shaping adaptive evolution in the natural world. Further to this, experimental evolution studies allow for tests of theories that may be difficult or impossible to test in natural populations for logistical and methodological reasons and can even generate new insights, suggesting further refinement of existing theories. However, as experimental evolution studies often take place in a very particular set of controlled conditions--that is simple environments, a small range of usually asexual species, relatively short timescales--the question remains as to how applicable these experimental results are to natural populations. In this review, we discuss important insights coming from experimental evolution, focusing on four key topics tied to the evolutionary genetics of adaptation, and within those topics, we discuss the extent to which the experimental work compliments and informs natural population studies. We finish by making suggestions for future work in particular a need for natural population genomic time series data, as well as the necessity for studies that combine both experimental evolution and natural population approaches.