Expression profiles reveal parallel evolution of epistatic interactions involving the CRP regulon in Escherichia coli

The extent and nature of epistatic interactions between mutations are issues of fundamental importance in evolutionary biology. However, they are difficult to study and their influence on adaptation remains poorly understood. Here, we use a systems-level approach to examine epistatic interactions that arose during the evolution of Escherichia coli in a defined environment. We used expression arrays to compare the effect on global patterns of gene expression of deleting a central regulatory gene, crp. Effects were measured in two lineages that had independently evolved for 20,000 generations and in their common ancestor. We found that deleting crp had a much more dramatic effect on the expression profile of the two evolved lines than on the ancestor. Because the sequence of the crp gene was unchanged during evolution, these differences indicate epistatic interactions between crp and mutations at other loci that accumulated during evolution. Moreover, a striking degree of parallelism was observed between the two independently evolved lines; 115 genes that were not crp-dependent in the ancestor became dependent on crp in both evolved lines. An analysis of changes in crp dependence of well-characterized regulons identified a number of regulatory genes as candidates for harboring beneficial mutations that could account for these parallel expression changes. Mutations within three of these genes have previously been found and shown to contribute to fitness. Overall, these findings indicate that epistasis has been important in the adaptive evolution of these lines, and they provide new insight into the types of genetic changes through which epistasis can evolve. More generally, we demonstrate that expression profiles can be profitably used to investigate epistatic interactions.     

Adaptive evolution of the lactose utilization network in experimentally evolved populations of Escherichia coli

Adaptation to novel environments is often associated with changes in gene regulation. Nevertheless, few studies have been able both to identify the genetic basis of changes in regulation and to demonstrate why these changes are beneficial. To this end, we have focused on understanding both how and why the lactose utilization network has evolved in replicate populations of Escherichia coli. We found that lac operon regulation became strikingly variable, including changes in the mode of environmental response (bimodal, graded, and constitutive), sensitivity to inducer concentration, and maximum expression level. In addition, some classes of regulatory change were enriched in specific selective environments. Sequencing of evolved clones, combined with reconstruction of individual mutations in the ancestral background, identified mutations within the lac operon that recapitulate many of the evolved regulatory changes. These mutations conferred fitness benefits in environments containing lactose, indicating that the regulatory changes are adaptive. The same mutations conferred different fitness effects when present in an evolved clone, indicating that interactions between the lac operon and other evolved mutations also contribute to fitness. Similarly, changes in lac regulation not explained by lac operon mutations also point to important interactions with other evolved mutations. Together these results underline how dynamic regulatory interactions can be, in this case evolving through mutations both within and external to the canonical lactose utilization network.     

Evolution of global regulatory networks during a long-term experiment with Escherichia coli

Evolution has shaped all living organisms on Earth, although many details of this process are shrouded in time. However, it is possible to see, with one's own eyes, evolution as it happens by performing experiments in defined laboratory conditions with microbes that have suitably fast generations. The longest-running microbial evolution experiment was started in 1988, at which time twelve populations were founded by the same strain of Escherichia coli. Since then, the populations have been serially propagated and have evolved for tens of thousands of generations in the same environment. The populations show numerous parallel phenotypic changes, and such parallelism is a hallmark of adaptive evolution. Many genetic targets of natural selection have been identified, revealing a high level of genetic parallelism as well. Beneficial mutations affect all levels of gene regulation in the cells including individual genes and operons all the way to global regulatory networks. Of particular interest, two highly interconnected networks -- governing DNA superhelicity and the stringent response -- have been demonstrated to be deeply involved in the phenotypic and genetic adaptation of these experimental populations.     

Evolution after Introduction of a Novel Metabolic Pathway Consistently Leads to Restoration of Wild-Type Physiology

Organisms cope with physiological stressors through acclimatizing mechanisms in the short-term and adaptive mechanisms over evolutionary timescales. During adaptation to an environmental or genetic perturbation, beneficial mutations can generate numerous physiological changes: some will be novel with respect to prior physiological states, while others might either restore acclimatizing responses to a wild-type state, reinforce them further, or leave them unchanged. We examined the interplay of acclimatizing and adaptive responses at the level of global gene expression in Methylobacterium extorquens AM1 engineered with a novel central metabolism. Replacing central metabolism with a distinct, foreign pathway resulted in much slower growth than wild-type. After 600 generations of adaptation, however, eight replicate populations founded from this engineered ancestor had improved up to 2.5-fold. A comparison of global gene expression in wild-type, engineered, and all eight evolved strains revealed that the vast majority of changes during physiological adaptation effectively restored acclimatizing processes to wild-type expression states. On average, 93% of expression perturbations from the engineered strain were restored, with 70% of these occurring in perfect parallel across all eight replicate populations. Novel changes were common but typically restricted to one or a few lineages, and reinforcing changes were quite rare. Despite this, cases in which expression was novel or reinforced in parallel were enriched for loci harboring beneficial mutations. One case of parallel, reinforced changes was the pntAB transhydrogenase that uses NADH to reduce NADP⁺ to NADPH. We show that PntAB activity was highly correlated with the restoration of NAD(H) and NADP(H) pools perturbed in the engineered strain to wild-type levels, and with improved growth. These results suggest that much of the evolved response to genetic perturbation was a consequence rather than a cause of adaptation and that physiology avoided “reinventing the wheel” by restoring acclimatizing processes to the pre-stressed state.     

Insertion Sequence-Driven Evolution of <Emphasis Type="Italic">Escherichia coli</Emphasis> in Chemostats

Insertion sequence (IS) elements are present in almost all bacterial genomes and are mobile enough to provide genomic tools to differentiate closely related isolates. They can be used to estimate genetic diversity and identify fitness-enhancing mutations during evolution experiments. Here, we determined the genomic distribution of eight IS elements in 120 genomes sampled from Escherichia coli populations that evolved in glucose- and phosphate-limited chemostats by comparison to the ancestral pattern. No significant differential transposition of the various IS types was detected across the environments. The phylogenies revealed substantial diversity amongst clones sampled from each chemostat, consistent with the phenotypic diversity within populations. Two IS-related changes were common to independent chemostats, suggesting parallel evolution. One of them corresponded to insertions of IS1 elements within rpoS encoding the master regulator of stress conditions. The other parallel event was an IS5-dependent deletion including mutY involved in DNA repair, thereby providing the molecular mechanism of generation of mutator clones in these evolving populations. These deletions occurred in different co-existing genotypes within single populations and were of various sizes. Moreover, differential locations of IS elements combined with their transpositional activity provided evolved clones with different phenotypic landscapes. Therefore, IS elements strongly influenced the evolutionary processes in continuous E. coli cultures by providing a way to modify both the global regulatory network and the mutation rates of evolving cells.     

cis-Regulatory changes in Kit ligand expression and parallel evolution of pigmentation in sticklebacks and humans

Dramatic pigmentation changes have evolved within most vertebrate groups, including fish and humans. Here we use genetic crosses in sticklebacks to investigate the parallel origin of pigmentation changes in natural populations. High-resolution mapping and expression experiments show that light gills and light ventrums map to a divergent regulatory allele of the Kit ligand (Kitlg) gene. The divergent allele reduces expression in gill and skin tissue and is shared by multiple derived freshwater populations with reduced pigmentation. In humans, Europeans and East Asians also share derived alleles at the KITLG locus. Strong signatures of selection map to regulatory regions surrounding the gene, and admixture mapping shows that the KITLG genomic region has a significant effect on human skin color. These experiments suggest that regulatory changes in Kitlg contribute to natural variation in vertebrate pigmentation, and that similar genetic mechanisms may underlie rapid evolutionary change in fish and humans.     

Identification of adaptive changes in an evolving population of Escherichia coli: the role of changes with regulatory and highly pleiotropic effects

A population of Escherichia coli initiated with a single clone developed extensive morphological and physiological polymorphism after being maintained for 773 generations in glucose-limited continuous culture. To understand the mechanisms of adaptation to this environment, total protein patterns of four adaptive clones and of the parent strains were examined by two-dimensional gel electrophoresis. Approximately 20% of the proteins (approximately 160 in absolute numbers) showed significantly different levels of expression in pairwise comparisons of parent and adapted clones. The extent of these changes points to the importance of mutations with regulatory and/or highly pleiotropic effects in the adaptive process. The four evolved clones all expressed fewer proteins than did the parent strain, supporting the hypothesis of energy conservation during evolutionary change. Forty-two proteins that could be assigned to known cellular functions were identified. The changes in some of them indicated that the evolved clones developed different adaptive mechanisms to glucose-limited environment. Changes were observed in the expression levels of proteins associated with translation, membrane composition, shock response, and active transport. A fraction of the changes could not be either explained or predicted from a consideration of the nature of the environment in which the clones evolved.     

Altered regulation of the OmpF porin by Fis in Escherichia coli during an evolution experiment and between B and K-12 strains

The phenotypic plasticity of global regulatory networks provides bacteria with rapid acclimation to a wide range of environmental conditions, while genetic changes in those networks provide additional flexibility as bacteria evolve across long time scales. We previously identified mutations in the global regulator-encoding gene fis that enhanced organismal fitness during a long-term evolution experiment with Escherichia coli. To gain insight into the effects of these mutations, we produced two-dimensional protein gels with strains carrying different fis alleles, including a beneficial evolved allele and one with an in-frame deletion. We found that Fis controls the expression of the major porin-encoding gene ompF in the E. coli B-derived ancestral strain used in the evolution experiment, a relationship that has not been described before. We further showed that this regulatory connection evolved over two different time scales, perhaps explaining why it was not observed before. On the longer time scale, we showed that this regulation of ompF by Fis is absent from the more widely studied K-12 strain and thus is specific to the B strain. On a shorter time scale, this regulatory linkage was lost during 20,000 generations of experimental evolution of the B strain. Finally, we mapped the Fis binding sites in the ompF regulatory region, and we present a hypothetical model of ompF expression that includes its other known regulators.     

Spontaneously arising mutL mutators in evolving Escherichia coli populations are the result of changes in repeat length

Over the course of thousands of generations of growth in a glucose-limited environment, 3 of 12 experimental populations of Escherichia coli spontaneously and independently evolved greatly increased mutation rates. In two of the populations, the mutations responsible for this increased mutation rate lie in the same region of the mismatch repair gene mutL. In this region, a 6-bp repeat is present in three copies in the gene of the wild-type ancestor of the experimental populations but is present in four copies in one of the experimental populations and two copies in the other. These in-frame mutations either add or delete the amino acid sequence LA in the MutL protein. We determined that the replacement of the wild-type sequence with either of these mutations was sufficient to increase the mutation rate of the wild-type strain to a level comparable to that of the mutator strains. Complementation of strains bearing the mutator mutations with wild-type copies of either mutL or the mismatch repair gene uvrD rescued the wild-type mutation rate. The position of the mutator mutations-in the region of MutL known as the ATP lid-suggests a possible deficiency in MutL's ATPase activity as the cause of the mutator phenotype. The similarity of the two mutator mutations (despite the independent evolutionary histories of the populations that gave rise to them) leads to a discussion of the potential adaptive role of DNA repeats.     

Insertion-sequence-mediated mutations isolated during adaptation to growth and starvation in Lactococcus lactis

We studied the activity of three multicopy insertion sequence (IS) elements in 12 populations of Lactococcus lactis IL1403 that evolved in the laboratory for 1000 generations under various environmental conditions (growth or starvation and shaken or stationary). Using RFLP analysis of single-clone representatives of each population, nine IS-mediated mutations were detected across all environmental conditions and all involving IS981. When it was assumed that these mutations were neutral, their frequency was higher under shaken than under stationary conditions, possibly due to oxygen stress. We characterized seven of the nine mutations at the molecular level and studied their population dynamics where possible. Two were simple insertions into new positions and the other five were recombinational deletions (of <1->10 kb) among existing and new copies of IS981; in all but one case these mutations disrupted gene functions. The best candidate beneficial mutations were two deletions of which similar versions were detected in two populations each. One of these two parallel deletions, affecting a gene involved in bacteriophage resistance, showed intermediate rearrangements and may also have resulted from increased local transposition rates.     

De Novo Evolution of Complex, Global and Hierarchical Gene Regulatory Mechanisms

Gene regulatory networks exhibit complex, hierarchical features such as global regulation and network motifs. There is much debate about whether the evolutionary origins of such features are the results of adaptation, or the by-products of non-adaptive processes of DNA replication. The lack of availability of gene regulatory networks of ancestor species on evolutionary timescales makes this a particularly difficult problem to resolve. Digital organisms, however, can be used to provide a complete evolutionary record of lineages. We use a biologically realistic evolutionary model that includes gene expression, regulation, metabolism and biosynthesis, to investigate the evolution of complex function in gene regulatory networks. We discover that: (i) network architecture and complexity evolve in response to environmental complexity, (ii) global gene regulation is selected for in complex environments, (iii) complex, inter-connected, hierarchical structures evolve in stages, with energy regulation preceding stress responses, and stress responses preceding growth rate adaptations and (iv) robustness of evolved models to mutations depends on hierarchical level: energy regulation and stress responses tend not to be robust to mutations, whereas growth rate adaptations are more robust and non-lethal when mutated. These results highlight the adaptive and incremental evolution of complex biological networks, and the value and potential of studying realistic in silico evolutionary systems as a way of understanding living systems.     

Fitness effects of fixed beneficial mutations in microbial populations

Beneficial mutations are intuitively relevant to understanding adaptation, yet not all beneficial mutations are of consequence to the long-term evolutionary outcome of adaptation. Many beneficial mutations-mostly those of small effect-are lost due either to (1) genetic drift or to (2) competition among clones carrying different beneficial mutations, a phenomenon called the "Hill-Robertson effect" for sexual populations and "clonal interference" for asexual populations. Competition among clones becomes more prevalent with increasing genetic linkage and increasing population size, and it is thus generally characteristic of microbial populations. Together, these two phenomena suggest that only those beneficial mutations of large fitness effect should achieve fixation, despite the fact that most beneficial mutations produced are predicted to have very small fitness effects. Here, we confirm this prediction-both empirically and theoretically-by showing that fitness effects of fixed beneficial mutations follow a distribution whose mode is positive.     

Genome Analysis of a Transmissible Lineage of Pseudomonas aeruginosa Reveals Pathoadaptive Mutations and Distinct Evolutionary Paths of Hypermutators

Genome sequencing of bacterial pathogens has advanced our understanding of their evolution, epidemiology, and response to antibiotic therapy. However, we still have only a limited knowledge of the molecular changes in in vivo evolving bacterial populations in relation to long-term, chronic infections. For example, it remains unclear what genes are mutated to facilitate the establishment of long-term existence in the human host environment, and in which way acquisition of a hypermutator phenotype with enhanced rates of spontaneous mutations influences the evolutionary trajectory of the pathogen. Here we perform a retrospective study of the DK2 clone type of P. aeruginosa isolated from Danish patients suffering from cystic fibrosis (CF), and analyze the genomes of 55 bacterial isolates collected from 21 infected individuals over 38 years. Our phylogenetic analysis of 8,530 mutations in the DK2 genomes shows that the ancestral DK2 clone type spread among CF patients through several independent transmission events. Subsequent to transmission, sub-lineages evolved independently for years in separate hosts, creating a unique possibility to study parallel evolution and identification of genes targeted by mutations to optimize pathogen fitness (pathoadaptive mutations). These genes were related to antibiotic resistance, the cell envelope, or regulatory functions, and we find that the prevalence of pathoadaptive mutations correlates with evolutionary success of co-evolving sub-lineages.The long-term co-existence of both normal and hypermutator populations enabled comparative investigations of the mutation dynamics in homopolymeric sequences in which hypermutators are particularly prone to mutations. We find a positive exponential correlation between the length of the homopolymer and its likelihood to acquire mutations and identify two homopolymer-containing genes preferentially mutated in hypermutators. This homopolymer facilitated differential mutagenesis provides a novel genome-wide perspective on the different evolutionary trajectories of hypermutators, which may help explain their emergence in CF infections.