Promoter recruitment drives the emergence of proto-genes in a long-term evolution experiment with Escherichia coli.

The phenomenon of de novo gene birth-the emergence of genes from non-genic sequences-has received considerable attention due to the widespread occurrence of genes that are unique to particular species or genomes. Most instances of de novo gene birth have been recognized through comparative analyses of genome sequences in eukaryotes, despite the abundance of novel, lineage-specific genes in bacteria and the relative ease with which bacteria can be studied in an experimental context. Here, we explore the genetic record of the Escherichia coli long-term evolution experiment (LTEE) for changes indicative of "proto-genic" phases of new gene birth in which non-genic sequences evolve stable transcription and/or translation. Over the time span of the LTEE, non-genic regions are frequently transcribed, translated and differentially expressed, with levels of transcription across low-expressed regions increasing in later generations of the experiment. Proto-genes formed downstream of new mutations result either from insertion element activity or chromosomal translocations that fused preexisting regulatory sequences to regions that were not expressed in the LTEE ancestor. Additionally, we identified instances of proto-gene emergence in which a previously unexpressed sequence was transcribed after formation of an upstream promoter, although such cases were rare compared to those caused by recruitment of preexisting promoters. Tracing the origin of the causative mutations, we discovered that most occurred early in the history of the LTEE, often within the first 20,000 generations, and became fixed soon after emergence. Our findings show that proto-genes emerge frequently within evolving populations, can persist stably, and can serve as potential substrates for new gene formation.

No assembly required: Time for stronger, simpler publishing standards for DNA sequences.

Uniformly accessible DNA sequences are needed to improve experimental reproducibility and automation. Rather than descriptions of how engineered DNA is assembled, publishers should require complete and empirically validated sequences.

Coevolutionary phage training leads to greater bacterial suppression and delays the evolution of phage resistance.

The evolution of antibiotic-resistant bacteria threatens to become the leading cause of worldwide mortality. This crisis has renewed interest in the practice of phage therapy. Yet, bacteria's capacity to evolve resistance may debilitate this therapy as well. To combat the evolution of phage resistance and improve treatment outcomes, many suggest leveraging phages' ability to counter resistance by evolving phages on target hosts before using them in therapy (phage training). We found that in vitro, λtrn, a phage trained for 28 d, suppressed bacteria ∼1,000-fold for three to eight times longer than its untrained ancestor. Prolonged suppression was due to a delay in the evolution of resistance caused by several factors. Mutations that confer resistance to λtrn are ∼100× less common, and while the target bacterium can evolve complete resistance to the untrained phage in a single step, multiple mutations are required to evolve complete resistance to λtrn. Mutations that confer resistance to λtrn are more costly than mutations for untrained phage resistance. Furthermore, when resistance does evolve, λtrn is better able to suppress these forms of resistance. One way that λtrn improved was through recombination with a gene in a defunct prophage in the host genome, which doubled phage fitness. This transfer of information from the host genome is an unexpected but highly efficient mode of training phage. Lastly, we found that many other independently trained λ phages were able to suppress bacterial populations, supporting the important role training could play during phage therapeutic development.

Genomic evolution of antibiotic resistance is contingent on genetic background following a long-term experiment with Escherichia coli.

Antibiotic resistance is a growing health concern. Efforts to control resistance would benefit from an improved ability to forecast when and how it will evolve. Epistatic interactions between mutations can promote divergent evolutionary trajectories, which complicates our ability to predict evolution. We recently showed that differences between genetic backgrounds can lead to idiosyncratic responses in the evolvability of phenotypic resistance, even among closely related Escherichia coli strains. In this study, we examined whether a strain's genetic background also influences the genotypic evolution of resistance. Do lineages founded by different genotypes take parallel or divergent mutational paths to achieve their evolved resistance states? We addressed this question by sequencing the complete genomes of antibiotic-resistant clones that evolved from several different genetic starting points during our earlier experiments. We first validated our statistical approach by quantifying the specificity of genomic evolution with respect to antibiotic treatment. As expected, mutations in particular genes were strongly associated with each drug. Then, we determined that replicate lines evolved from the same founding genotypes had more parallel mutations at the gene level than lines evolved from different founding genotypes, although these effects were more subtle than those showing antibiotic specificity. Taken together with our previous work, we conclude that historical contingency can alter both genotypic and phenotypic pathways to antibiotic resistance.

The Genomic Basis of Evolutionary Novelties in a Leafhopper.

Evolutionary innovations generate phenotypic and species diversity. Elucidating the genomic processes underlying such innovations is central to understanding biodiversity. In this study, we addressed the genomic basis of evolutionary novelties in the glassy-winged sharpshooter (Homalodisca vitripennis, GWSS), an agricultural pest. Prominent evolutionary innovations in leafhoppers include brochosomes, proteinaceous structures that are excreted and used to coat the body, and obligate symbiotic associations with two bacterial types that reside within cytoplasm of distinctive cell types. Using PacBio long-read sequencing and Dovetail Omni-C technology, we generated a chromosome-level genome assembly for the GWSS and then validated the assembly using flow cytometry and karyotyping. Additional transcriptomic and proteomic data were used to identify novel genes that underlie brochosome production. We found that brochosome-associated genes include novel gene families that have diversified through tandem duplications. We also identified the locations of genes involved in interactions with bacterial symbionts. Ancestors of the GWSS acquired bacterial genes through horizontal gene transfer (HGT), and these genes appear to contribute to symbiont support. Using a phylogenomics approach, we inferred HGT sources and timing. We found that some HGT events date to the common ancestor of the hemipteran suborder Auchenorrhyncha, representing some of the oldest known examples of HGT in animals. Overall, we show that evolutionary novelties in leafhoppers are generated by the combination of acquiring novel genes, produced both de novo and through tandem duplication, acquiring new symbiotic associations that enable use of novel diets and niches, and recruiting foreign genes to support symbionts and enhance herbivory.

Electron Tomography and Machine Learning for Understanding the Highly Ordered Structure of Leafhopper Brochosomes.

Insects known as leafhoppers (Hemiptera: Cicadellidae) produce hierarchically structured nanoparticles known as brochosomes that are exuded and applied to the insect cuticle, thereby providing camouflage and anti-wetting properties to aid insect survival. Although the physical properties of brochosomes are thought to depend on the leafhopper species, the structure-function relationships governing brochosome behavior are not fully understood. Brochosomes have complex hierarchical structures and morphological heterogeneity across species, due to which a multimodal characterization approach is required to effectively elucidate their nanoscale structure and properties. In this work, we study the structural and mechanical properties of brochosomes using a combination of atomic force microscopy (AFM), electron microscopy (EM), electron tomography, and machine learning (ML)-based quantification of large and complex scanning electron microscopy (SEM) image data sets. This suite of techniques allows for the characterization of internal and external brochosome structures, and ML-based image analysis methods of large data sets reveal correlations in the structure across several leafhopper species. Our results show that brochosomes are relatively rigid hollow spheres with characteristic dimensions and morphologies that depend on leafhopper species. Nanomechanical mapping AFM is used to determine a characteristic compression modulus for brochosomes on the order of 1-3 GPa, which is consistent with crystalline proteins. Overall, this work provides an improved understanding of the structural and mechanical properties of leafhopper brochosomes using a new set of ML-based image classification tools that can be broadly applied to nanostructured biological materials.

The dynamics of molecular evolution over 60,000 generations.

The outcomes of evolution are determined by a stochastic dynamical process that governs how mutations arise and spread through a population. However, it is difficult to observe these dynamics directly over long periods and across entire genomes. Here we analyse the dynamics of molecular evolution in twelve experimental populations of Escherichia coli, using whole-genome metagenomic sequencing at five hundred-generation intervals through sixty thousand generations. Although the rate of fitness gain declines over time, molecular evolution is characterized by signatures of rapid adaptation throughout the duration of the experiment, with multiple beneficial variants simultaneously competing for dominance in each population. Interactions between ecological and evolutionary processes play an important role, as long-term quasi-stable coexistence arises spontaneously in most populations, and evolution continues within each clade. We also present evidence that the targets of natural selection change over time, as epistasis and historical contingency alter the strength of selection on different genes. Together, these results show that long-term adaptation to a constant environment can be a more complex and dynamic process than is often assumed.

pLannotate: engineered plasmid annotation.

Engineered plasmids are widely used in the biological sciences. Since many plasmids contain DNA sequences that have been reused and remixed by researchers for decades, annotation of their functional elements is often incomplete. Missing information about the presence, location, or precise identity of a plasmid feature can lead to unintended consequences or failed experiments. Many engineered plasmids contain sequences-such as recombinant DNA from all domains of life, wholly synthetic DNA sequences, and engineered gene expression elements-that are not predicted by microbial genome annotation pipelines. Existing plasmid annotation tools have limited feature libraries and do not detect incomplete fragments of features that are present in many plasmids for historical reasons and may impact their newly designed functions. We created the open source pLannotate web server so users can quickly and comprehensively annotate plasmid features. pLannotate is powered by large databases of genetic parts and proteins. It employs a filtering algorithm to display only the most relevant feature matches and also reports feature fragments. Finally, pLannotate displays a graphical map of the annotated plasmid, explains the provenance of each feature prediction, and allows results to be downloaded in a variety of formats. The webserver for pLannotate is accessible at: http://plannotate.barricklab.org/.

Gut microbe Lactiplantibacillus plantarum undergoes different evolutionary trajectories between insects and mammals.

BACKGROUND: Animals form complex symbiotic associations with their gut microbes, whose evolution is determined by an intricate network of host and environmental factors. In many insects, such as Drosophila melanogaster, the microbiome is flexible, environmentally determined, and less diverse than in mammals. In contrast, mammals maintain complex multispecies consortia that are able to colonize and persist in the gastrointestinal tract. Understanding the evolutionary and ecological dynamics of gut microbes in different hosts is challenging. This requires disentangling the ecological factors of selection, determining the timescales over which evolution occurs, and elucidating the architecture of such evolutionary patterns. RESULTS: We employ experimental evolution to track the pace of the evolution of a common gut commensal, Lactiplantibacillus plantarum, within invertebrate (Drosophila melanogaster) and vertebrate (Mus musculus) hosts and their respective diets. We show that in Drosophila, the nutritional environment dictates microbial evolution, while the host benefits L. plantarum growth only over short ecological timescales. By contrast, in a mammalian animal model, L. plantarum evolution results to be divergent between the host intestine and its diet, both phenotypically (i.e., host-evolved populations show higher adaptation to the host intestinal environment) and genomically. Here, both the emergence of hypermutators and the high persistence of mutated genes within the host's environment strongly differed from the low variation observed in the host's nutritional environment alone. CONCLUSIONS: Our results demonstrate that L. plantarum evolution diverges between insects and mammals. While the symbiosis between Drosophila and L. plantarum is mainly determined by the host diet, in mammals, the host and its intrinsic factors play a critical role in selection and influence both the phenotypic and genomic evolution of its gut microbes, as well as the outcome of their symbiosis.

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. Here 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 in which 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.

Rapid and assured genetic engineering methods applied to Acinetobacter baylyi ADP1 genome streamlining.

One goal of synthetic biology is to improve the efficiency and predictability of living cells by removing extraneous genes from their genomes. We demonstrate improved methods for engineering the genome of the metabolically versatile and naturally transformable bacterium Acinetobacter baylyi ADP1 and apply them to a genome streamlining project. In Golden Transformation, linear DNA fragments constructed by Golden Gate Assembly are directly added to cells to create targeted deletions, edits, or additions to the chromosome. We tested the dispensability of 55 regions of the ADP1 chromosome using Golden Transformation. The 18 successful multiple-gene deletions ranged in size from 21 to 183 kb and collectively accounted for 23.4% of its genome. The success of each multiple-gene deletion attempt could only be partially predicted on the basis of an existing collection of viable ADP1 single-gene deletion strains and a new transposon insertion sequencing (Tn-Seq) dataset that we generated. We further show that ADP1's native CRISPR/Cas locus is active and can be retargeted using Golden Transformation. We reprogrammed it to create a CRISPR-Lock, which validates that a gene has been successfully removed from the chromosome and prevents it from being reacquired. These methods can be used together to implement combinatorial routes to further genome streamlining and for more rapid and assured metabolic engineering of this versatile chassis organism.

Engineering a Culturable Serratia symbiotica Strain for Aphid Paratransgenesis.

Aphids are global agricultural pests and important models for bacterial symbiosis. To date, none of the native symbionts of aphids have been genetically manipulated, which limits our understanding of how they interact with their hosts. Serratia symbiotica CWBI-2.3T is a culturable, gut-associated bacterium isolated from the black bean aphid. Closely related Serratia symbiotica strains are facultative aphid endosymbionts that are vertically transmitted from mother to offspring during embryogenesis. We demonstrate that CWBI-2.3T can be genetically engineered using a variety of techniques, plasmids, and gene expression parts. Then, we use fluorescent protein expression to track the dynamics with which CWBI-2.3T colonizes the guts of multiple aphid species, and we measure how this bacterium affects aphid fitness. Finally, we show that we can induce heterologous gene expression from engineered CWBI-2.3T in living aphids. These results inform the development of CWBI-2.3T for aphid paratransgenesis, which could be used to study aphid biology and enable future agricultural technologies.IMPORTANCE Insects have remarkably diverse and integral roles in global ecosystems. Many harbor symbiotic bacteria, but very few of these bacteria have been genetically engineered. Aphids are major agricultural pests and an important model system for the study of symbiosis. This work describes methods for engineering a culturable aphid symbiont, Serratia symbiotica CWBI-2.3T These approaches and genetic tools could be used in the future to implement new paradigms for the biological study and control of aphids.

Innovation in an E. coli evolution experiment is contingent on maintaining adaptive potential until competition subsides.

Key innovations are disruptive evolutionary events that enable a species to escape constraints and rapidly diversify. After 15 years of the Lenski long-term evolution experiment with Escherichia coli, cells in one of the twelve populations evolved the ability to utilize citrate, an abundant but previously untapped carbon source in the environment. Descendants of these cells became dominant in the population and subsequently diversified as a consequence of invading this vacant niche. Mutations responsible for the appearance of rudimentary citrate utilization and for refining this ability have been characterized. However, the complete nature of the genetic and/or ecological events that set the stage for this key innovation is unknown. In particular, it is unclear why it took so long for citrate utilization to evolve and why it still has evolved in only one of the twelve E. coli populations after 30 years of the Lenski experiment. In this study, we recapitulated the initial mutation needed to evolve citrate utilization in strains isolated from throughout the first 31,500 generations of the history of this population. We found that there was already a slight fitness benefit for this mutation in the original ancestor of the evolution experiment and in other early isolates. However, evolution of citrate utilization was blocked at this point due to competition with other mutations that improved fitness in the original niche. Subsequently, an anti-potentiated genetic background evolved in which it was deleterious to evolve rudimentary citrate utilization. Only later, after further mutations accumulated that restored the benefit of this first-step mutation and the overall rate of adaptation in the population slowed, was citrate utilization likely to evolve. Thus, intense competition and the types of mutations that it favors can lead to short-sighted evolutionary trajectories that hide a stepping stone needed to access a key innovation from many future generations.

Arresting Evolution.

Evolution in the form of selective breeding has long been harnessed as a useful tool by humans. However, rapid evolution can also be a danger to our health and a stumbling block for biotechnology. Unwanted evolution can underlie the emergence of drug and pesticide resistance, cancer, and weeds. It makes live vaccines and engineered cells inherently unreliable and unpredictable, and therefore potentially unsafe. Yet, there are strategies that have been and can possibly be used to stop or slow many types of evolution. We review and classify existing population genetics-inspired methods for arresting evolution. Then, we discuss how genome editing techniques enable a radically new set of approaches to limit evolution.

Directed evolution of Escherichia coli with lower-than-natural plasmid mutation rates.

Unwanted evolution of designed DNA sequences limits metabolic and genome engineering efforts. Engineered functions that are burdensome to host cells and slow their replication are rapidly inactivated by mutations, and unplanned mutations with unpredictable effects often accumulate alongside designed changes in large-scale genome editing projects. We developed a directed evolution strategy, Periodic Reselection for Evolutionarily Reliable Variants (PResERV), to discover mutations that prolong the function of a burdensome DNA sequence in an engineered organism. Here, we used PResERV to isolate Escherichia coli cells that replicate ColE1-type plasmids with higher fidelity. We found mutations in DNA polymerase I and in RNase E that reduce plasmid mutation rates by 6- to 30-fold. The PResERV method implicitly selects to maintain the growth rate of host cells, and high plasmid copy numbers and gene expression levels are maintained in some of the evolved E. coli strains, indicating that it is possible to improve the genetic stability of cellular chassis without encountering trade-offs in other desirable performance characteristics. Utilizing these new antimutator E. coli and applying PResERV to other organisms in the future promises to prevent evolutionary failures and unpredictability to provide a more stable genetic foundation for synthetic biology.

Specificity of genome evolution in experimental populations of Escherichia coli evolved at different temperatures.

Isolated populations derived from a common ancestor are expected to diverge genetically and phenotypically as they adapt to different local environments. To examine this process, 30 populations of Escherichia coli were evolved for 2,000 generations, with six in each of five different thermal regimes: constant 20 °C, 32 °C, 37 °C, 42 °C, and daily alternations between 32 °C and 42 °C. Here, we sequenced the genomes of one endpoint clone from each population to test whether the history of adaptation in different thermal regimes was evident at the genomic level. The evolved strains had accumulated ∼5.3 mutations, on average, and exhibited distinct signatures of adaptation to the different environments. On average, two strains that evolved under the same regime exhibited ∼17% overlap in which genes were mutated, whereas pairs that evolved under different conditions shared only ∼4%. For example, all six strains evolved at 32 °C had mutations in nadR, whereas none of the other 24 strains did. However, a population evolved at 37 °C for an additional 18,000 generations eventually accumulated mutations in the signature genes strongly associated with adaptation to the other temperature regimes. Two mutations that arose in one temperature treatment tended to be beneficial when tested in the others, although less so than in the regime in which they evolved. These findings demonstrate that genomic signatures of adaptation can be highly specific, even with respect to subtle environmental differences, but that this imprint may become obscured over longer timescales as populations continue to change and adapt to the shared features of their environments.

Bioassay for Determining the Concentrations of Caffeine and Individual Methylxanthines in Complex Samples.

Caffeine and other methylxanthines are stimulant molecules found in formulated beverages, including sodas and energy drinks, and in brewed beverages, such as coffee and teas. Previously, we developed a bioassay for caffeine that involves monitoring the growth of a ΔguaB mutant of Escherichia coli defective in de novo guanine biosynthesis. When supplemented with a plasmid expressing the genes for an N-demethylation pathway from Pseudomonas putida CBB5, these bacteria demethylate caffeine (1,3,7-trimethylxanthine) and other methylxanthines into xanthine, which is then converted into guanine to support cell growth. A major limitation of this bioassay was that it could only measure the total concentration of all methylxanthines in a mixture. Therefore, it could not be used to measure the caffeine content of beverages like teas, which contain substantial quantities of multiple methylxanthines. To overcome this limitation, we created seven new plasmids containing all subsets of the three demethylase genes (ndmA, ndmB, and ndmC). We show that strains of ΔguaBE. coli containing each plasmid are able to demethylate specific subsets of methylxanthines and that they can be used to determine the concentrations of individual methylxanthines in complex mixtures containing multiple methylxanthines, including coffee doped with an additional methylxanthine. While validating this assay, we also discovered an unexpected demethylation event at the 1-methyl position when NdmB and NdmC were expressed in the absence of NdmA. The improved cell-based bioassay is inexpensive, is easy to use, and gives results comparable to standard high-performance liquid chromatography methods for measuring methylxanthine concentrations.IMPORTANCE Caffeine (1,3,7-trimethylxanthine) is the dominant neurostimulant found in coffee, teas, sodas, and energy drinks. Measuring the amount of caffeine and other methylxanthines in these beverages is important for quality assurance and safety in food science. Methylxanthines are also used in medicine and as performance-enhancing drugs, two contexts in which accurately determining their concentrations in bodily fluids is important. Liquid chromatography is the standard method for measuring methylxanthine concentrations in a sample, but it requires specialized equipment and expertise. We improved a previous bioassay that links E. coli growth to methylxanthine demethylation so that it can now be used to determine the amounts of individual methylxanthines in complex mixtures or beverages, such as coffee.

Genome dynamics during experimental evolution.

Evolutionary changes in organismal traits may occur either gradually or suddenly. However, until recently, there has been little direct information about how phenotypic changes are related to the rate and the nature of the underlying genotypic changes. Technological advances that facilitate whole-genome and whole-population sequencing, coupled with experiments that 'watch' evolution in action, have brought new precision to and insights into studies of mutation rates and genome evolution. In this Review, we discuss the evolutionary forces and ecological processes that govern genome dynamics in various laboratory systems in the context of relevant population genetic theory, and we relate these findings to evolution in natural populations.

Genomic analysis of a key innovation in an experimental Escherichia coli population.

Evolutionary novelties have been important in the history of life, but their origins are usually difficult to examine in detail. We previously described the evolution of a novel trait, aerobic citrate utilization (Cit(+)), in an experimental population of Escherichia coli. Here we analyse genome sequences to investigate the history and genetic basis of this trait. At least three distinct clades coexisted for more than 10,000 generations before its emergence. The Cit(+) trait originated in one clade by a tandem duplication that captured an aerobically expressed promoter for the expression of a previously silent citrate transporter. The clades varied in their propensity to evolve this novel trait, although genotypes able to do so existed in all three clades, implying that multiple potentiating mutations arose during the population's history. Our findings illustrate the importance of promoter capture and altered gene regulation in mediating the exaptation events that often underlie evolutionary innovations.

Large-scale analysis of post-translational modifications in E. coli under glucose-limiting conditions.

BACKGROUND: Post-translational modification (PTM) of proteins is central to many cellular processes across all domains of life, but despite decades of study and a wealth of genomic and proteomic data the biological function of many PTMs remains unknown. This is especially true for prokaryotic PTM systems, many of which have only recently been recognized and studied in depth. It is increasingly apparent that a deep sampling of abundance across a wide range of environmental stresses, growth conditions, and PTM types, rather than simply cataloging targets for a handful of modifications, is critical to understanding the complex pathways that govern PTM deposition and downstream effects. RESULTS: We utilized a deeply-sampled dataset of MS/MS proteomic analysis covering 9 timepoints spanning the Escherichia coli growth cycle and an unbiased PTM search strategy to construct a temporal map of abundance for all PTMs within a 400 Da window of mass shifts. Using this map, we are able to identify novel targets and temporal patterns for N-terminal N α acetylation, C-terminal glutamylation, and asparagine deamidation. Furthermore, we identify a possible relationship between N-terminal N α acetylation and regulation of protein degradation in stationary phase, pointing to a previously unrecognized biological function for this poorly-understood PTM. CONCLUSIONS: Unbiased detection of PTM in MS/MS proteomics data facilitates the discovery of novel modification types and previously unobserved dynamic changes in modification across growth timepoints.