Programming Dynamic Division of Labor Using Horizontal Gene Transfer.

The metabolic engineering of microbes has broad applications, including biomanufacturing, bioprocessing, and environmental remediation. The introduction of a complex, multistep pathway often imposes a substantial metabolic burden on the host cell, restraining the accumulation of productive biomass and limiting pathway efficiency. One strategy to alleviate metabolic burden is the division of labor (DOL) in which different subpopulations carry out different parts of the pathway and work together to convert a substrate into a final product. However, the maintenance of different engineered subpopulations is challenging due to competition and convoluted interstrain population dynamics. Through modeling, we show that dynamic division of labor (DDOL), which we define as the DOL between indiscrete populations capable of dynamic and reversible interchange, can overcome these limitations and enable the robust maintenance of burdensome, multistep pathways. We propose that DDOL can be mediated by horizontal gene transfer (HGT) and use plasmid genomics to uncover evidence that DDOL is a strategy utilized by natural microbial communities. Our work suggests that bioengineers can harness HGT to stabilize synthetic metabolic pathways in microbial communities, enabling the development of robust engineered systems for deployment in a variety of contexts.

Duplicated antibiotic resistance genes reveal ongoing selection and horizontal gene transfer in bacteria.

Horizontal gene transfer (HGT) and gene duplication are often considered as separate mechanisms driving the evolution of new functions. However, the mobile genetic elements (MGEs) implicated in HGT can copy themselves, so positive selection on MGEs could drive gene duplications. Here, we use a combination of modeling and experimental evolution to examine this hypothesis and use long-read genome sequences of tens of thousands of bacterial isolates to examine its generality in nature. Modeling and experiments show that antibiotic selection can drive the evolution of duplicated antibiotic resistance genes (ARGs) through MGE transposition. A key implication is that duplicated ARGs should be enriched in environments associated with antibiotic use. To test this, we examined the distribution of duplicated ARGs in 18,938 complete bacterial genomes with ecological metadata. Duplicated ARGs are highly enriched in bacteria isolated from humans and livestock. Duplicated ARGs are further enriched in an independent set of 321 antibiotic-resistant clinical isolates. Our findings indicate that duplicated genes often encode functions undergoing positive selection and horizontal gene transfer in microbial communities.

Programming dynamic division of labor using horizontal gene transfer.

The metabolic engineering of microbes has broad applications, including in biomanufacturing, bioprocessing, and environmental remediation. The introduction of a complex, multi-step pathway often imposes a substantial metabolic burden on the host cell, restraining the accumulation of productive biomass and limiting pathway efficiency. One strategy to alleviate metabolic burden is division of labor (DOL), in which different subpopulations carry out different parts of the pathway and work together to convert a substrate into a final product. However, the maintenance of different engineered subpopulations is challenging due to competition and convoluted inter-strain population dynamics. Through modeling, we show that dynamic division of labor (DDOL) mediated by horizontal gene transfer (HGT) can overcome these limitations and enable the robust maintenance of burdensome, multi-step pathways. We also use plasmid genomics to uncover evidence that DDOL is a strategy utilized by natural microbial communities. Our work suggests that bioengineers can harness HGT to stabilize synthetic metabolic pathways in microbial communities, enabling the development of robust engineered systems for deployment in a variety of contexts.

Idiosyncratic Purifying Selection on Metabolic Enzymes in the Long-Term Evolution Experiment with Escherichia coli.

Bacteria, Archaea, and Eukarya all share a common set of metabolic reactions. This implies that the function and topology of central metabolism has been evolving under purifying selection over deep time. Central metabolism may similarly evolve under purifying selection during long-term evolution experiments, although it is unclear how long such experiments would have to run (decades, centuries, millennia) before signs of purifying selection on metabolism appear. I hypothesized that central and superessential metabolic enzymes would show evidence of purifying selection in the long-term evolution experiment with Escherichia coli (LTEE). I also hypothesized that enzymes that specialize on single substrates would show stronger evidence of purifying selection in the LTEE than generalist enzymes that catalyze multiple reactions. I tested these hypotheses by analyzing metagenomic time series covering 62,750 generations of the LTEE. I find mixed support for these hypotheses, because the observed patterns of purifying selection are idiosyncratic and population-specific. To explain this finding, I propose the Jenga hypothesis, named after a children's game in which blocks are removed from a tower until it falls. The Jenga hypothesis postulates that loss-of-function mutations degrade costly, redundant, and non-essential metabolic functions. Replicate populations can therefore follow idiosyncratic trajectories of lost redundancies, despite purifying selection on overall function. I tested the Jenga hypothesis by simulating the evolution of 1,000 minimal genomes under strong purifying selection. As predicted, the minimal genomes converge to different metabolic networks. Strikingly, the core genes common to all 1,000 minimal genomes show consistent signatures of purifying selection in the LTEE.

Discovery of positive and purifying selection in metagenomic time series of hypermutator microbial populations.

A general method to infer both positive and purifying selection during the real-time evolution of hypermutator pathogens would be broadly useful. To this end, we introduce a Simple Test to Infer Mode of Selection (STIMS) from metagenomic time series of evolving microbial populations. We test STIMS on metagenomic data generated by simulations of bacterial evolution, and on metagenomic data spanning 62,750 generations of Lenski's long-term evolution experiment with Escherichia coli (LTEE). This benchmarking shows that STIMS detects positive selection in both nonmutator and hypermutator populations, and purifying selection in hypermutator populations. Using STIMS, we find strong evidence of ongoing positive selection on key regulators of the E. coli gene regulatory network, even in some hypermutator populations. STIMS also detects positive selection on regulatory genes in hypermutator populations of Pseudomonas aeruginosa that adapted to subinhibitory concentrations of colistin-an antibiotic of last resort-for just twenty-six days of laboratory evolution. Our results show that the fine-tuning of gene regulatory networks is a general mechanism for rapid and ongoing adaptation. The simplicity of STIMS, together with its intuitive visual interpretation, make it a useful test for positive and purifying selection in metagenomic data sets that track microbial evolution in real-time.

Intra- and interpopulation transposition of mobile genetic elements driven by antibiotic selection.

The spread of genes encoding antibiotic resistance is often mediated by horizontal gene transfer (HGT). Many of these genes are associated with transposons, a type of mobile genetic element that can translocate between the chromosome and plasmids. It is widely accepted that the translocation of antibiotic resistance genes onto plasmids potentiates their spread by HGT. However, it is unclear how this process is modulated by environmental factors, especially antibiotic treatment. To address this issue, we asked whether antibiotic exposure would select for the transposition of resistance genes from chromosomes onto plasmids and, if so, whether antibiotic concentration could tune the distribution of resistance genes between chromosomes and plasmids. We addressed these questions by analysing the transposition dynamics of synthetic and natural transposons that encode resistance to different antibiotics. We found that stronger antibiotic selection leads to a higher fraction of cells carrying the resistance on plasmids because the increased copy number of resistance genes on multicopy plasmids leads to higher expression of those genes and thus higher cell survival when facing antibiotic selection. Once they have transposed to plasmids, antibiotic resistance genes are primed for rapid spread by HGT. Our results provide quantitative evidence for a mechanism by which antibiotic selection accelerates the spread of antibiotic resistance in microbial communities.

Maintenance of Metabolic Plasticity despite Relaxed Selection in a Long-Term Evolution Experiment with Escherichia coli.

AbstractTraits that are unused in a given environment are subject to processes that tend to erode them, leading to reduced fitness in other environments. Although this general tendency is clear, we know much less about why some traits are lost while others are retained and about the roles of mutation and selection in generating different responses. We addressed these issues by examining populations of a facultative anaerobe, Escherichia coli, that have evolved for >30 years in the presence of oxygen, with relaxed selection for anaerobic growth and the associated metabolic plasticity. We asked whether evolution led to the loss, improvement, or maintenance of anaerobic growth, and we analyzed gene expression and mutational data sets to understand the outcomes. We identified genomic signatures of both positive and purifying selection on aerobic-specific genes, while anaerobic-specific genes showed clear evidence of relaxed selection. We also found parallel evolution at two interacting loci that regulate anaerobic growth. We competed the ancestor and evolved clones from each population in an anoxic environment, and we found that anaerobic fitness had not decayed, despite relaxed selection. In summary, relaxed selection does not necessarily reduce an organism's fitness in other environments. Instead, the genetic architecture of the traits under relaxed selection and their correlations with traits under positive and purifying selection may sometimes determine evolutionary outcomes.

Selection Maintains Protein Interactome Resilience in the Long-Term Evolution Experiment with Escherichia coli.

Most cellular functions are carried out by a dynamic network of interacting proteins. An open question is whether the network properties of protein interactomes represent phenotypes under natural selection. One proposal is that protein interactomes have evolved to be resilient, such that they tend to maintain connectivity when proteins are removed from the network. This hypothesis predicts that interactome resilience should be maintained by natural selection during long-term experimental evolution. I tested this prediction by modeling the evolution of protein-protein interaction (PPI) networks in Lenski's long-term evolution experiment with Escherichia coli (LTEE). In this test, I removed proteins affected by nonsense, insertion, deletion, and transposon mutations in evolved LTEE strains, and measured the resilience of the resulting networks. I compared the rate of change of network resilience in each LTEE population to the rate of change of network resilience for corresponding randomized networks. The evolved PPI networks are significantly more resilient than networks in which random proteins have been deleted. Moreover, the evolved networks are generally more resilient than networks in which the random deletion of proteins was restricted to those disrupted in LTEE. These results suggest that evolution in the LTEE has favored PPI networks that are, on average, more resilient than expected from the genetic variation across the evolved strains. My findings therefore support the hypothesis that selection maintains protein interactome resilience over evolutionary time.

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.

Large-scale discovery of protein interactions at residue resolution using co-evolution calculated from genomic sequences.

Increasing numbers of protein interactions have been identified in high-throughput experiments, but only a small proportion have solved structures. Recently, sequence coevolution-based approaches have led to a breakthrough in predicting monomer protein structures and protein interaction interfaces. Here, we address the challenges of large-scale interaction prediction at residue resolution with a fast alignment concatenation method and a probabilistic score for the interaction of residues. Importantly, this method (EVcomplex2) is able to assess the likelihood of a protein interaction, as we show here applied to large-scale experimental datasets where the pairwise interactions are unknown. We predict 504 interactions de novo in the E. coli membrane proteome, including 243 that are newly discovered. While EVcomplex2 does not require available structures, coevolving residue pairs can be used to produce structural models of protein interactions, as done here for membrane complexes including the Flagellar Hook-Filament Junction and the Tol/Pal complex.

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.

Genomic and phenotypic evolution of Escherichia coli in a novel citrate-only resource environment.

Evolutionary innovations allow populations to colonize new ecological niches. We previously reported that aerobic growth on citrate (Cit+) evolved in an Escherichia coli population during adaptation to a minimal glucose medium containing citrate (DM25). Cit+ variants can also grow in citrate-only medium (DM0), a novel environment for E. coli. To study adaptation to this niche, we founded two sets of Cit+ populations and evolved them for 2500 generations in DM0 or DM25. The evolved lineages acquired numerous parallel mutations, many mediated by transposable elements. Several also evolved amplifications of regions containing the maeA gene. Unexpectedly, some evolved populations and clones show apparent declines in fitness. We also found evidence of substantial cell death in Cit+ clones. Our results thus demonstrate rapid trait refinement and adaptation to the new citrate niche, while also suggesting a recalcitrant mismatch between E. coli physiology and growth on citrate.

Genome-wide discovery of epistatic loci affecting antibiotic resistance in Neisseria gonorrhoeae using evolutionary couplings.

Genome analysis should allow the discovery of interdependent loci that together cause antibiotic resistance. In practice, however, the vast number of possible epistatic interactions erodes statistical power. Here, we extend an approach that has been successfully used to identify epistatic residues in proteins to infer genomic loci that are strongly coupled. This approach reduces the number of tests required for an epistatic genome-wide association study of antibiotic resistance and increases the likelihood of identifying causal epistasis. We discovered 38 loci and 240 epistatic pairs that influence the minimum inhibitory concentrations of 5 different antibiotics in 1,102 isolates of Neisseria gonorrhoeae that were confirmed in a second dataset of 495 isolates. Many known resistance-affecting loci were recovered; however, the majority of associations occurred in unreported genes, such as murE. About half of the discovered epistasis involved at least one locus previously associated with antibiotic resistance, including interactions between gyrA and parC. Still, many combinations involved unreported loci and genes. While most variation in minimum inhibitory concentrations could be explained by identified loci, epistasis substantially increased explained phenotypic variance. Our work provides a systematic identification of epistasis affecting antibiotic resistance in N. gonorrhoeae and a generalizable approach for epistatic genome-wide association studies.

Gain-of-function experiments with bacteriophage lambda uncover residues under diversifying selection in nature.

Viral gain-of-function mutations frequently evolve during laboratory experiments. Whether the specific mutations that evolve in the lab also evolve in nature and whether they have the same impact on evolution in the real world is unknown. We studied a model virus, bacteriophage λ, that repeatedly evolves to exploit a new host receptor under typical laboratory conditions. Here, we demonstrate that two residues of λ's J protein are required for the new function. In natural λ variants, these amino acid sites are highly diverse and evolve at high rates. Insertions and deletions at these locations are associated with phylogenetic patterns indicative of ecological diversification. Our results show that viral evolution in the laboratory mirrors that in nature and that laboratory experiments can be coupled with protein sequence analyses to identify the causes of viral evolution in the real world. Furthermore, our results provide evidence for widespread host-shift evolution in lambdoid viruses.

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.

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.

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.

Haplotype-Phased Synthetic Long Reads from Short-Read Sequencing.

Next-generation DNA sequencing has revolutionized the study of biology. However, the short read lengths of the dominant instruments complicate assembly of complex genomes and haplotype phasing of mixtures of similar sequences. Here we demonstrate a method to reconstruct the sequences of individual nucleic acid molecules up to 11.6 kilobases in length from short (150-bp) reads. We show that our method can construct 99.97%-accurate synthetic reads from bacterial, plant, and animal genomic samples, full-length mRNA sequences from human cancer cell lines, and individual HIV env gene variants from a mixture. The preparation of multiple samples can be multiplexed into a single tube, further reducing effort and cost relative to competing approaches. Our approach generates sequencing libraries in three days from less than one microgram of DNA in a single-tube format without custom equipment or specialized expertise.

Sustained fitness gains and variability in fitness trajectories in the long-term evolution experiment with Escherichia coli.

Many populations live in environments subject to frequent biotic and abiotic changes. Nonetheless, it is interesting to ask whether an evolving population's mean fitness can increase indefinitely, and potentially without any limit, even in a constant environment. A recent study showed that fitness trajectories of Escherichia coli populations over 50 000 generations were better described by a power-law model than by a hyperbolic model. According to the power-law model, the rate of fitness gain declines over time but fitness has no upper limit, whereas the hyperbolic model implies a hard limit. Here, we examine whether the previously estimated power-law model predicts the fitness trajectory for an additional 10 000 generations. To that end, we conducted more than 1100 new competitive fitness assays. Consistent with the previous study, the power-law model fits the new data better than the hyperbolic model. We also analysed the variability in fitness among populations, finding subtle, but significant, heterogeneity in mean fitness. Some, but not all, of this variation reflects differences in mutation rate that evolved over time. Taken together, our results imply that both adaptation and divergence can continue indefinitely--or at least for a long time--even in a constant environment.

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

Synonymous genetic differences vary by more than 20-fold among genes in natural isolates of Escherichia coli. One hypothesis to explain this heterogeneity is that genes with high levels of synonymous variation mutate at higher rates than genes with low synonymous variation. If so, then one would expect to observe similar mutational patterns in evolution experiments. In fact, however, the pattern of synonymous substitutions in a long-term evolution experiment with E. coli does not support this hypothesis. In particular, the extent of synonymous variation across genes in that experiment does not reflect the variation observed in natural isolates of E. coli. Instead, gene length alone predicts with high accuracy the prevalence of synonymous changes in the experimental populations. We hypothesize that patterns of synonymous variation in natural E. coli populations are instead caused by differences across genomic regions in their effective population size that, in turn, reflect different histories of recombination, horizontal gene transfer, selection, and population structure.