Person:

Marx, Christopher J

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.

How do adapting populations navigate the tensions between the costs of gene expression and the benefits of gene products to optimize the levels of many genes at once? Here we combined independently-arising beneficial mutations that altered enzyme levels in the central metabolism of Methylobacterium extorquens to uncover the fitness landscape defined by gene expression levels. We found strong antagonism and sign epistasis between these beneficial mutations. Mutations with the largest individual benefit interacted the most antagonistically with other mutations, a trend we also uncovered through analyses of datasets from other model systems. However, these beneficial mutations interacted multiplicatively (i.e., no epistasis) at the level of enzyme expression. By generating a model that predicts fitness from enzyme levels we could explain the observed sign epistasis as a result of overshooting the optimum defined by a balance between enzyme catalysis benefits and fitness costs. Knowledge of the phenotypic landscape also illuminated that, although the fitness peak was phenotypically far from the ancestral state, it was not genetically distant. Single beneficial mutations jumped straight toward the global optimum rather than being constrained to change the expression phenotypes in the correlated fashion expected by the genetic architecture. Given that adaptation in nature often results from optimizing gene expression, these conclusions can be widely applicable to other organisms and selective conditions. Poor interactions between individually beneficial alleles affecting gene expression may thus compromise the benefit of sex during adaptation and promote genetic differentiation.

Methylobacterium extorquens AM1, a strain serendipitously isolated half a century ago, has become the best-characterized model system for the study of aerobic methylotrophy (the ability to grow on reduced single-carbon compounds). However, with 5 replicons and 174 insertion sequence (IS) elements in the genome as well as a long history of domestication in the laboratory, genetic and genomic analysis of M. extorquens AM1 face several challenges. On the contrary, a recently isolated strain - M. extorquens PA1- is closely related to M. extorquens AM1 (100% 16S rRNA identity) and contains a streamlined genome with a single replicon and only 20 IS elements. With the exception of the methylamine dehydrogenase encoding gene cluster (mau), genes known to be involved in methylotrophy are well conserved between M. extorquens AM1 and M. extorquens PA1. In this paper we report four primary findings regarding methylotrophy in PA1. First, with a few notable exceptions, the repertoire of methylotrophy genes between PA1 and AM1 is extremely similar. Second, PA1 grows faster with higher yields compared to AM1 on C1 and multi-C substrates in minimal media, but AM1 grows faster in rich medium. Third, deletion mutants in PA1 throughout methylotrophy modules have the same C1 growth phenotypes observed in AM1. Finally, the precision of our growth assays revealed several unexpected growth phenotypes for various knockout mutants that serve as leads for future work in understanding their basis and generality across Methylobacterium strains.

Evolutionary adaptation to a constant environment is often accompanied by specialization and a reduction of fitness in other environments. We assayed the ability of the Lenski Escherichia coli populations to grow on a range of carbon sources after 50,000 generations of adaptation on glucose. Using direct measurements of growth rates, we demonstrated that declines in performance were much less widespread than suggested by previous results from Biolog assays of cellular respiration. Surprisingly, there were many performance increases on a variety of substrates. In addition to the now famous example of citrate, we observed several other novel gains of function for organic acids that the ancestral strain only marginally utilized. Quantitative growth data also showed that strains with a higher mutation rate exhibited significantly more declines, suggesting that most metabolic erosion was driven by mutation accumulation and not by physiological tradeoffs. These reductions in growth by mutator strains were ameliorated by growth at lower temperature, consistent with the hypothesis that this metabolic erosion is largely caused by destabilizing mutations to the associated enzymes. We further hypothesized that reductions in growth rate would be greatest for substrates used most differently from glucose, and we used flux balance analysis to formulate this question quantitatively. To our surprise, we found no significant relationship between decreases in growth and dissimilarity to glucose metabolism. Taken as a whole, these data suggest that in a single resource environment, specialization does not mainly result as an inevitable consequence of adaptive tradeoffs, but rather due to the gradual accumulation of disabling mutations in unused portions of the genome.

Few areas of science have benefited more from the expansion in sequencing capability than the study of microbial communities. Can sequence data, besides providing hypotheses of the functions the members possess, detect the evolutionary and ecological processes that are occurring? For example, can we determine if a species is adapting to one niche, or if it is diversifying into multiple specialists that inhabit distinct niches? Fortunately, adaptation of populations in the laboratory can serve as a model to test our ability to make such inferences about evolution and ecology from sequencing. Even adaptation to a single niche can give rise to complex temporal dynamics due to the transient presence of multiple competing lineages. If there are multiple niches, this complexity is augmented by segmentation of the population into multiple specialists that can each continue to evolve within their own niche. For a known example of parallel diversification that occurred in the laboratory, sequencing data gave surprisingly few obvious, unambiguous signs of the ecological complexity present. Whereas experimental systems are open to direct experimentation to test hypotheses of selection or ecological interaction, the difficulty in “seeing ecology” from sequencing for even such a simple system suggests translation to communities like the human microbiome will be quite challenging. This will require both improved empirical methods to enhance the depth and time resolution for the relevant polymorphisms and novel statistical approaches to rigorously examine time-series data for signs of various evolutionary and ecological phenomena within and between species.

The most powerful genome-scale framework to model metabolism, flux balance analysis (FBA), is an evolutionary optimality model. It hypothesizes selection upon a proposed optimality criterion in order to predict the set of internal fluxes that would maximize fitness. Here we present a direct test of the optimality assumption underlying FBA by comparing the central metabolic fluxes predicted by multiple criteria to changes measurable by a 13C-labeling method for experimentally-evolved strains. We considered datasets for three Escherichia coli evolution experiments that varied in their length, consistency of environment, and initial optimality. For ten populations that were evolved for 50,000 generations in glucose minimal medium, we observed modest changes in relative fluxes that led to small, but significant decreases in optimality and increased the distance to the predicted optimal flux distribution. In contrast, seven populations evolved on the poor substrate lactate for 900 generations collectively became more optimal and had flux distributions that moved toward predictions. For three pairs of central metabolic knockouts evolved on glucose for 600–800 generations, there was a balance between cases where optimality and flux patterns moved toward or away from FBA predictions. Despite this variation in predictability of changes in central metabolism, two generalities emerged. First, improved growth largely derived from evolved increases in the rate of substrate use. Second, FBA predictions bore out well for the two experiments initiated with ancestors with relatively sub-optimal yield, whereas those begun already quite optimal tended to move somewhat away from predictions. These findings suggest that the tradeoff between rate and yield is surprisingly modest. The observed positive correlation between rate and yield when adaptation initiated further from the optimum resulted in the ability of FBA to use stoichiometric constraints to predict the evolution of metabolism despite selection for rate.

Methylobacterium extorquens strains are the best-studied methylotrophic model system, and their metabolism of single carbon compounds has been studied for over 50 years. Here we develop a new system for high-throughput batch culture of M. extorquens in microtiter plates by jointly optimizing the properties of the organism, the growth media and the culturing system. After removing cellulose synthase genes in M. extorquens strains AM1 and PA1 to prevent biofilm formation, we found that currently available lab automation equipment, integrated and managed by open source software, makes possible reliable estimates of the exponential growth rate. Using this system, we developed an optimized growth medium for M. extorquens using response surface methodologies. We found that media that used EDTA as a metal chelator inhibited growth and led to inconsistent culture conditions. In contrast, the new medium we developed with a PIPES buffer and metals chelated by citrate allowed for fast and more consistent growth rates. This new Methylobacterium PIPES (‘MP’) medium was also robust to large deviations in its component ingredients which avoided batch effects from experiments that used media prepared at different times. MP medium allows for faster and more consistent growth than other media used for M. extorquens.

Understanding evolutionary dynamics within microbial populations requires the ability to accurately follow allele frequencies through time. Here we present a rapid, cost-effective method (FREQ-Seq) that leverages Illumina next-generation sequencing for localized, quantitative allele frequency detection. Analogous to RNA-Seq, FREQ-Seq relies upon counts from the >105 reads generated per locus per time-point to determine allele frequencies. Loci of interest are directly amplified from a mixed population via two rounds of PCR using inexpensive, user-designed oligonucleotides and a bar-coded bridging primer system that can be regenerated in-house. The resulting bar-coded PCR products contain the adapters needed for Illumina sequencing, eliminating further library preparation. We demonstrate the utility of FREQ-Seq by determining the order and dynamics of beneficial alleles that arose as a microbial population, founded with an engineered strain of Methylobacterium, evolved to grow on methanol. Quantifying allele frequencies with minimal bias down to 1% abundance allowed effective analysis of SNPs, small in-dels and insertions of transposable elements. Our data reveal large-scale clonal interference during the early stages of adaptation and illustrate the utility of FREQ-Seq as a cost-effective tool for tracking allele frequencies in populations.