Fast growth increases the selective advantage of a mutation arising recurrently during evolution under metal limitation

Environmental modulation of global epistasis in a drug resistance fitness landscape

Interactions between mutations (epistasis) can add substantial complexity to genotype-phenotype maps, hampering our ability to predict evolution. Yet, recent studies have shown that the fitness effect of a mutation can often be predicted from the fitness of its genetic background using simple, linear relationships. This phenomenon, termed global epistasis, has been leveraged to reconstruct fitness landscapes and infer adaptive trajectories in a wide variety of contexts. However, little attention has been paid to how patterns of global epistasis may be affected by environmental variation, despite this variation frequently being a major driver of evolution. This is particularly relevant for the evolution of drug resistance, where antimicrobial drugs may change the environment faced by pathogens and shape their adaptive trajectories in ways that can be difficult to predict. By analyzing a fitness landscape of four mutations in a gene encoding an essential enzyme of P. falciparum (a parasite cause of malaria), here we show that patterns of global epistasis can be strongly modulated by the concentration of a drug in the environment. Expanding on previous theoretical results, we demonstrate that this modulation can be quantitatively explained by how specific gene-by-gene interactions are modified by drug dose. Importantly, our results highlight the need to incorporate potential environmental variation into the global epistasis framework in order to predict adaptation in dynamic environments.

Epistasis and evolution: recent advances and an outlook for prediction

As organisms evolve, the effects of mutations change as a result of epistatic interactions with other mutations accumulated along the line of descent. This can lead to shifts in adaptability or robustness that ultimately shape subsequent evolution. Here, we review recent advances in measuring, modeling, and predicting epistasis along evolutionary trajectories, both in microbial cells and single proteins. We focus on simple patterns of global epistasis that emerge in this data, in which the effects of mutations can be predicted by a small number of variables. The emergence of these patterns offers promise for efforts to model epistasis and predict evolution.

Global epistasis on fitness landscapes

Epistatic interactions between mutations add substantial complexity to adaptive landscapes and are often thought of as detrimental to our ability to predict evolution. Yet, patterns of global epistasis, in which the fitness effect of a mutation is well-predicted by the fitness of its genetic background, may actually be of help in our efforts to reconstruct fitness landscapes and infer adaptive trajectories. Microscopic interactions between mutations, or inherent nonlinearities in the fitness landscape, may cause global epistasis patterns to emerge. In this brief review, we provide a succinct overview of recent work about global epistasis, with an emphasis on building intuition about why it is often observed. To this end, we reconcile simple geometric reasoning with recent mathematical analyses, using these to explain why different mutations in an empirical landscape may exhibit different global epistasis patterns-ranging from diminishing to increasing returns. Finally, we highlight open questions and research directions. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Mutational robustness changes during long-term adaptation in laboratory budding yeast populations

As an adapting population traverses the fitness landscape, its local neighborhood (i.e., the collection of fitness effects of single-step mutations) can change shape because of interactions with mutations acquired during evolution. These changes to the distribution of fitness effects can affect both the rate of adaptation and the accumulation of deleterious mutations. However, while numerous models of fitness landscapes have been proposed in the literature, empirical data on how this distribution changes during evolution remains limited. In this study, we directly measure how the fitness landscape neighborhood changes during laboratory adaptation. Using a barcode-based mutagenesis system, we measure the fitness effects of 91 specific gene disruption mutations in genetic backgrounds spanning 8000–10,000 generations of evolution in two constant environments. We find that the mean of the distribution of fitness effects decreases in one environment, indicating a reduction in mutational robustness, but does not change in the other. We show that these distribution-level patterns result from differences in the relative frequency of certain patterns of epistasis at the level of individual mutations, including fitness-correlated and idiosyncratic epistasis.

A Periplasmic Lanthanide Mediator, Lanmodulin, in Methylobacterium aquaticum Strain 22A

Methylobacterium and Methylorubrum species oxidize methanol via pyrroloquinoline quinone-methanol dehydrogenases (MDHs). MDHs can be classified into two major groups, Ca²⁺-dependent MDH (MxaF) and lanthanide (Ln³⁺)-dependent MDH (XoxF), whose expression is regulated by the availability of Ln³⁺. A set of a siderophore, TonB-dependent receptor, and an ABC transporter that resembles the machinery for iron uptake is involved in the solubilization and transport of Ln³⁺. The transport of Ln³⁺ into the cytosol enhances XoxF expression. A unique protein named lanmodulin from Methylorubrum extorquens strain AM1 was identified as a specific Ln³⁺-binding protein, and its biological function was implicated to be an Ln³⁺ shuttle in the periplasm. In contrast, it remains unclear how Ln³⁺ levels in the cells are maintained, because Ln³⁺ is potentially deleterious to cellular systems due to its strong affinity to phosphate ions. In this study, we investigated the function of a lanmodulin homolog in Methylobacterium aquaticum strain 22A. The expression of a gene encoding lanmodulin (lanM) was induced in response to the presence of La³⁺. A recombinant LanM underwent conformational change upon La³⁺ binding. Phenotypic analyses on lanM deletion mutant and overexpressing strains showed that LanM is not necessary for the wild-type and XoxF-dependent mutant’s methylotrophic growth. We found that lanM expression was regulated by MxcQE (a two-component regulator for MxaF) and TonB_Ln (a TonB-dependent receptor for Ln³⁺). The expression level of mxcQE was altered to be negatively dependent on Ln³⁺ concentration in ∆lanM, whereas it was constant in the wild type. Furthermore, when exposed to La³⁺, ∆lanM showed an aggregating phenotype, cell membrane impairment, La deposition in the periplasm evidenced by electron microscopy, differential expression of proteins involved in membrane integrity and phosphate starvation, and possibly lower La content in the membrane vesicle (MV) fractions. Taken together, we concluded that lanmodulin is involved in the complex regulation mechanism of MDHs and homeostasis of cellular Ln levels by facilitating transport and MV-mediated excretion.

Genetic Context Significantly Influences the Maintenance and Evolution of Degenerate Pathways

Understanding the evolution of novel physiological traits is highly relevant for expanding the characterization and manipulation of biological systems. Acquisition of new traits can be achieved through horizontal gene transfer (HGT). Here, we investigate drivers that promote or deter the maintenance of HGT-driven degeneracy, occurring when processes accomplish identical functions through nonidentical components. Subsequent evolution can optimize newly acquired functions; for example, beneficial alleles identified in an engineered Methylorubrum extorquens strain allowed it to utilize a “Foreign” formaldehyde oxidation pathway substituted for its Native pathway for methylotrophic growth. We examined the fitness consequences of interactions between these alleles when they were combined with the Native pathway or both (Dual) pathways. Unlike the Foreign pathway context where they evolved, these alleles were often neutral or deleterious when moved into these alternative genetic backgrounds. However, there were instances where combinations of multiple alleles resulted in higher fitness outcomes than individual allelic substitutions could provide. Importantly, the genetic context accompanying these allelic substitutions significantly altered the fitness landscape, shifting local fitness peaks and restricting the set of accessible evolutionary trajectories. These findings highlight how genetic context can negatively impact the probability of maintaining native and HGT-introduced functions together, making it difficult for degeneracy to evolve. However, in cases where the cost of maintaining degeneracy was mitigated by adding evolved alleles impacting the function of these pathways, we observed rare opportunities for pathway coevolution to occur. Together, our results highlight the importance of genetic context and resulting epistasis in retaining or losing HGT-acquired degenerate functions.

Reasons for 2-furaldehyde and 5-hydroxymethyl-2-furaldehyde resistance in Saccharomyces cerevisiae: current state of knowledge and perspectives for further improvements

Common toxic compounds 2-furaldehyde (furfural) and 5-hydroxymethyl-2-furaldehyde (HMF) are formed from dehydration of pentose and hexose, respectively, during decomposition of lignocellulosic biomass polymers. Furfural and HMF represent a major class of aldehyde toxic chemicals that inhibit microbial growth and interfere with subsequent fermentation for production of renewable fuels and chemicals. Understanding mechanisms of yeast tolerance aids development of tolerant strains as the most economic means to overcome the toxicity. This review updates current knowledge on yeast resistance to these toxic chemicals obtained from rapid advances in the past few years. Findings are largely exemplified by an adapted strain NRRL Y-50049 compared with its progenitor, the industrial yeast Saccharomyces cerevisiae type strain NRRL Y-12632. Newly characterized molecular phenotypes distinguished acquired resistant components of Y-50049 from innate stress response of its progenitor Y-12632. These findings also raised important questions on how to address more deeply ingrained changes in addition to local renovations for yeast adaptation. An early review on understandings of yeast tolerance to these inhibitory compounds is available and its contents omitted here to avoid redundancy. Controversial and confusing issues on identification of yeast resistance to furfural and HMF are further clarified aiming improved future research. Propositions and perspectives on research understanding molecular mechanisms of yeast resistance and future improvements are also presented. KEY POINTS: • Distinguished adapted resistance from innate stress response in yeast. • Defined pathway-based molecular phenotypes of yeast resistance. • Proposed genomic insight and perspectives on yeast resistance and adaptation.

A glimpse of potential transposable element impact on adaptation of the industrial yeast Saccharomyces cerevisiae

The adapted industrial yeast strain Saccharomyces cerevisiae NRRL Y-50049 is able to in situ detoxify major toxic aldehyde compounds derived from sugar conversion of lignocellulosic biomass while producing ethanol. Pathway-based studies on its mechanisms of tolerance have been reported previously, however, little is known about transposable element (TE) involvement in its adaptation to inhibitory compounds. This work presents a comparative dynamic transcription expression analysis in response to a toxic treatment between Y-50049 and its progenitor, an industrial type strain NRRL Y-12632, using a time-course study. At least 77 TEs from Y-50049 showed significantly increased expression compared with its progenitor, especially during the late lag phase. Sequence analysis revealed significant differences in TE sequences between the two strains. Y-50049 was also found to have a transposons of yeast 2 (Ty2) long terminal repeat-linked YAT1 gene showing significantly higher copy number changes than its progenitor. These results raise awareness of potential TE involvement in the adaptation of industrial yeast to the tolerance of toxic chemicals.

Lanthanide-dependent alcohol dehydrogenases require an essential aspartate residue for metal coordination and enzymatic function

The lanthanide elements (Ln3+), those with atomic numbers 57-63 (excluding promethium, Pm3+), form a cofactor complex with pyrroloquinoline quinone (PQQ) in bacterial XoxF methanol dehydrogenases (MDHs) and ExaF ethanol dehydrogenases (EDHs), expanding the range of biological elements and opening novel areas of metabolism and ecology. Other MDHs, known as MxaFIs, are related in sequence and structure to these proteins, yet they instead possess a Ca2+-PQQ cofactor. An important missing piece of the Ln3+ puzzle is defining what features distinguish enzymes that use Ln3+-PQQ cofactors from those that do not. Here, using XoxF1 MDH from the model methylotrophic bacterium Methylorubrum extorquens AM1, we investigated the functional importance of a proposed lanthanide-coordinating aspartate residue. We report two crystal structures of XoxF1, one with and another without PQQ, both with La3+ bound in the active-site region and coordinated by Asp320 Using constructs to produce either recombinant XoxF1 or its D320A variant, we show that Asp320 is needed for in vivo catalytic function, in vitro activity, and La3+ coordination. XoxF1 and XoxF1 D320A, when produced in the absence of La3+, coordinated Ca2+ but exhibited little or no catalytic activity. We also generated the parallel substitution in ExaF to produce ExaF D319S and found that this variant loses the capacity for efficient ethanol oxidation with La3+ These results provide evidence that a Ln3+-coordinating aspartate is essential for the enzymatic functions of XoxF MDHs and ExaF EDHs, supporting the notion that sequences of these enzymes, and the genes that encode them, are markers for Ln3+ metabolism.

Establishment of CRISPR interference in Methylorubrum extorquens and application of rapidly mining a new phytoene desaturase involved in carotenoid biosynthesis

The methylotrophic bacterium Methylorubrum extorquens AM1 holds a great potential of a microbial cell factory in producing high value chemicals with methanol as the sole carbon and energy source. However, many gene functions remain unknown, hampering further rewiring of metabolic networks. Clustered regularly interspaced short palindromic repeat interference (CRISPRi) has been demonstrated to be a robust tool for gene knockdown in diverse organisms. In this study, we developed an efficient CRISPRi system through optimizing the promoter strength of Streptococcus pyogenes-derived deactivated cas9 (dcas9). When the dcas9 and sgRNA were respectively controlled by medium PR/tetO and strong PmxaF-g promoters, dynamic repression efficacy of cell growth through disturbing a central metabolism gene glyA was achieved from 41.9 to 96.6% dependent on the sgRNA targeting sites. Furthermore, the optimized CRISPRi system was shown to effectively decrease the abundance of exogenous fluorescent protein gene mCherry over 50% and to reduce the expression of phytoene desaturase gene crtI by 97.7%. We then used CRISPRi technology combined with 26 sgRNAs pool to rapidly discover a new phytoene desaturase gene META1_3670 from 2470 recombinant mutants. The gene function was further verified through gene deletion and complementation as well as phylogenetic tree analysis. In addition, we applied CRISPRi to repress the transcriptional level of squalene-hopene cyclase gene shc involved in hopanoid biosynthesis by 64.9%, which resulted in enhancing 1.9-fold higher of carotenoid production without defective cell growth. Thus, the CRISPRi system developed here provides a useful tool in mining functional gene of M. extorquens as well as in biotechnology for producing high-valued chemicals from methanol. KEY POINTS: Developing an efficient CRISPRi to knockdown gene expression in C1-utilizing bacteria CRISPRi combined with sgRNAs pool to rapidly discover a new phytoene desaturase gene Improvement of carotenoid production by repressing a competitive pathway.

Patterns and Mechanisms of Diminishing Returns from Beneficial Mutations

Diminishing returns epistasis causes the benefit of the same advantageous mutation smaller in fitter genotypes and is frequently observed in experimental evolution. However, its occurrence in other contexts, environment dependence, and mechanistic basis are unclear. Here, we address these questions using 1,005 sequenced segregants generated from a yeast cross. Under each of 47 examined environments, 66-92% of tested polymorphisms exhibit diminishing returns epistasis. Surprisingly, improving environment quality also reduces the benefits of advantageous mutations even when fitness is controlled for, indicating the necessity to revise the global epistasis hypothesis. We propose that diminishing returns originates from the modular organization of life where the contribution of each functional module to fitness is determined jointly by the genotype and environment and has an upper limit, and demonstrate that our model predictions match empirical observations. These findings broaden the concept of diminishing returns epistasis, reveal its generality and potential cause, and have important evolutionary implications.

Effects of Beneficial Mutations in pykF Gene Vary over Time and across Replicate Populations in a Long-Term Experiment with Bacteria

The fitness effects of mutations can depend on the genetic backgrounds in which they occur and thereby influence future opportunities for evolving populations. In particular, mutations that fix in a population might change the selective benefit of subsequent mutations, giving rise to historical contingency. We examine these effects by focusing on mutations in a key metabolic gene, pykF, that arose independently early in the history of 12 Escherichia coli populations during a long-term evolution experiment. Eight different evolved nonsynonymous mutations conferred similar fitness benefits of ∼10% when transferred into the ancestor, and these benefits were greater than the one conferred by a deletion mutation. In contrast, the same mutations had highly variable fitness effects, ranging from ∼0% to 25%, in evolved clones isolated from the populations at 20,000 generations. Two mutations that were moved into these evolved clones conferred similar fitness effects in a given clone, but different effects between the clones, indicating epistatic interactions between the evolved pykF alleles and the other mutations that had accumulated in each evolved clone. We also measured the fitness effects of six evolved pykF alleles in the same populations in which they had fixed, but at seven time points between 0 and 50,000 generations. Variation in fitness effects was high at intermediate time points, and declined to a low level at 50,000 generations, when the mean fitness effect was lowest. Our results demonstrate the importance of genetic context in determining the fitness effects of different beneficial mutations even within the same gene.

Power law fitness landscapes and their ability to predict fitness

Whether or not evolution by natural selection is predictable depends on the existence of general patterns shaping the way mutations interact with the genetic background. This interaction, also known as epistasis, has been observed during adaptation (macroscopic epistasis) and in individual mutations (microscopic epistasis). Interestingly, a consistent negative correlation between the fitness effect of beneficial mutations and background fitness (known as diminishing returns epistasis) has been observed across different species and conditions. We tested whether the adaptation pattern of an additional species, Schizosaccharomyces pombe, followed the same trend. We used strains that differed by the presence of large karyotype differences and observed the same pattern of fitness convergence. Using these data along with published datasets, we measured the ability of different models to describe adaptation rates. We found that a phenotype-fitness landscape shaped like a power law is able to correctly predict adaptation dynamics in a variety of species and conditions. Furthermore we show that this model can provide a link between the observed macroscopic and microscopic epistasis. It may be very useful in the development of algorithms able to predict the adaptation of microorganisms from measures of the current phenotypes. Overall, our results suggest that even though adaptation quickly slows down, populations adapting to lab conditions may be quite far from a fitness peak.

Instability of the 16S rRNA methyltransferase-encoding npmA gene: why have bacterial cells possessing npmA not spread despite their high and broad resistance to aminoglycosides?

The NpmA bacterial 16S rRNA methyltransferase, which is identified from Escherichia coli strains, confers high resistance to many types of aminoglycoside upon its host cells. But despite its resistance-conferring ability, only two cases of its isolation from E. coli (14 years apart) have been reported to date. Here, we investigated the effect of the npmA gene on aminoglycoside resistance in Pseudomonas aeruginosa and Klebsiella pneumoniae and its stability in E. coli cells by comparing it with armA, another 16S rRNA methyltransferase gene currently spreading globally. As a result, we found that npmA conferred resistance to all types of aminoglycoside antibiotics we tested (except streptomycin) in both P. aeruginosa and K. pneumoniae, as well in E. coli. In addition, co-expression of armA and npmA resulted in an additive effect for the resistance. However, in return for the resistance, we also observed that the growth rates and the cell survivability of the strains transformed with the npmA-harboring plasmids were inferior than those of the control strains and that these plasmids were easily disrupted by IS10, IS1, and IS5 insertion sequences. We discuss these data in the context of the threat posed by pathogenic strains possessing npmA.

Different adaptive strategies in E. coli populations evolving under macronutrient limitation and metal ion limitation

BACKGROUND

Adaptive responses to nutrient limitation involve mutations that increase the efficiency of usage or uptake of the limiting nutrient. However, starvation of different nutrients has contrasting effects on physiology, resulting in different evolutionary responses. Most studies performed to understand these evolutionary responses have focused only on macronutrient limitation. Hence our understanding of adaptation under limitation of other forms of nutrients is limited. In this study, we compared the evolutionary response in populations evolving under growth-limiting conditions for a macronutrient and a major cation.

RESULTS

We evolved eight populations of E. coli in nutrient-limited chemostats for 400 generations to identify the genetic basis of the mechanisms involved in efficient usage of two nutrients: nitrogen and magnesium. Our population genomic sequencing work, based on this study and previous work, allowed us to identify targets of selection under these nutrient limiting conditions. Global transcriptional regulators glnGL were targets of selection under nitrogen starvation, while proteins involved in outer-membrane biogenesis (genes from the lpt operon) were targets of selection under magnesium starvation. The protein involved in cell-cycle arrest (yhaV) was a target of selection in both environments. We re-constructed specific mutants to analyze the effect of individual mutations on fitness in nutrient limiting conditions in chemostats and in batch cultures. We further demonstrated that adaptation to nitrogen starvation proceeds via a nutrient specific mechanism, while that to magnesium starvation involves a more general mechanism.

CONCLUSIONS

Our results show two different forms of adaptive strategies under limitation of nutrients that effect cellular physiology in different ways. Adaptation to nitrogen starvation proceeds by upregulation of transcriptional regulator glnG and subsequently of transporter protein amtB, both of which results in increased nitrogen scavenging ability of the cell. On the other hand, adaptation to magnesium starvation proceeds via the restructuring of the cell outer-membrane, allowing magnesium to be redistributed to other biological processes. Also, adaptation to the chemostat environment involves selection for loss of function mutations in genes that under nutrient-limiting conditions interfere with continuous growth.

Hitchhiking and epistasis give rise to cohort dynamics in adapting populations

Beneficial mutations are the driving force of adaptive evolution. In asexual populations, the identification of beneficial alleles is confounded by the presence of genetically linked hitchhiker mutations. Parallel evolution experiments enable the recognition of common targets of selection; yet these targets are inherently enriched for genes of large target size and mutations of large effect. A comprehensive study of individual mutations is necessary to create a realistic picture of the evolutionarily significant spectrum of beneficial mutations. Here we use a bulk-segregant approach to identify the beneficial mutations across 11 lineages of experimentally evolved yeast populations. We report that nearly 80% of detected mutations have no discernible effects on fitness and less than 1% are deleterious. We determine the distribution of driver and hitchhiker mutations in 31 mutational cohorts, groups of mutations that arise synchronously from low frequency and track tightly with one another. Surprisingly, we find that one-third of cohorts lack identifiable driver mutations. In addition, we identify intracohort synergistic epistasis between alleles of hsl7 and kel1, which arose together in a low-frequency lineage.

Transposable Elements Mediate Adaptive Debilitation of Flagella in Experimental Escherichia coli Populations

Although insertion sequence (IS) elements are generally considered genomic parasites, they can mediate adaptive genetic changes in bacterial genomes. We discovered that among 12 laboratory-evolved Escherichia coli populations, three had experienced at least six different IS1-mediated deletions of flagellar genes. These deletions all involved the master flagellar regulator flhDC, and as such completely incapacitate motility. Two lines of evidence strongly suggest that these deletions were adaptive in our evolution experiment: (1) parallel evolution in three independent populations is highly unlikely just by chance, and (2) one of these deletion mutations swept to fixation within ~1000 generations, which is over two million times faster than expected if this deletion was instead selectively neutral and thus evolving by genetic drift. Because flagella are energetically expensive to synthesize and operate, we suspect that debilitating their construction conferred a fitness advantage in our well-stirred evolution experiment. These findings underscore the important role that IS elements can play in mediating adaptive loss-of-function mutations in bacteria.

Identification of the potentiating mutations and synergistic epistasis that enabled the evolution of inter-species cooperation

Microbes often engage in cooperation through releasing biosynthetic compounds required by other species to grow. Given that production of costly biosynthetic metabolites is generally subjected to multiple layers of negative feedback, single mutations may frequently be insufficient to generate cooperative phenotypes. Synergistic epistatic interactions between multiple coordinated changes may thus often underlie the evolution of cooperation through overproduction of metabolites. To test the importance of synergistic mutations in cooperation we used an engineered bacterial consortium of an Escherichia coli methionine auxotroph and Salmonella enterica. S. enterica relies on carbon by-products from E. coli if lactose is the only carbon source. Directly selecting wild-type S. enterica in an environment that favored cooperation through secretion of methionine only once led to a methionine producer, and this producer both took a long time to emerge and was not very effective at cooperating. On the other hand, when an initial selection for resistance of S. enterica to a toxic methionine analog, ethionine, was used, subsequent selection for cooperation with E. coli was rapid, and the resulting double mutants were much more effective at cooperation. We found that potentiating mutations in metJ increase expression of metA, which encodes the first step of methionine biosynthesis. This increase in expression is required for the previously identified actualizing mutations in metA to generate cooperation. This work highlights that where biosynthesis of metabolites involves multiple layers of regulation, significant secretion of those metabolites may require multiple mutations, thereby constraining the evolution of cooperation.

Regulation of solvent tolerance in Pseudomonas putida S12 mediated by mobile elements

Organic solvent‐tolerant bacteria are outstanding and versatile hosts for the bio‐based production of a broad range of generally toxic aromatic compounds. The energetically costly solvent tolerance mechanisms are subject to multiple levels of regulation, involving among other mobile genetic elements. The genome of the solvent‐tolerant Pseudomonas putida S12 contains many such mobile elements that play a major role in the regulation and adaptation to various stress conditions, including the regulation of expression of the solvent efflux pump SrpABC. We recently sequenced the genome of P. putida S12. Detailed annotation identified a threefold higher copy number of the mobile element ISS12 in contrast to earlier observations. In this study, we describe the mobile genetic elements and elaborate on the role of ISS12 in the establishment and maintenance of solvent tolerance in P. putida. We identified three different variants of ISS12 of which a single variant exhibits a high translocation rate. One copy of this variant caused a loss of solvent tolerance in the sequenced strain by disruption of srpA. Solvent tolerance could be restored by applying selective pressure, leading to a clean excision of the mobile element.

Hopanoid-free Methylobacterium extorquens DM4 overproduces carotenoids and has widespread growth impairment

Hopanoids are sterol-like membrane lipids widely used as geochemical proxies for bacteria. Currently, the physiological role of hopanoids is not well understood, and this represents one of the major limitations in interpreting the significance of their presence in ancient or contemporary sediments. Previous analyses of mutants lacking hopanoids in a range of bacteria have revealed a range of phenotypes under normal growth conditions, but with most having at least an increased sensitivity to toxins and osmotic stress. We employed hopanoid-free strains of Methylobacterium extorquens DM4, uncovering severe growth defects relative to the wild-type under many tested conditions, including normal growth conditions without additional stressors. Mutants overproduce carotenoids-the other major isoprenoid product of this strain-and show an altered fatty acid profile, pronounced flocculation in liquid media, and lower growth yields than for the wild-type strain. The flocculation phenotype can be mitigated by addition of cellulase to the medium, suggesting a link between the function of hopanoids and the secretion of cellulose in M. extorquens DM4. On solid media, colonies of the hopanoid-free mutant strain were smaller than wild-type, and were more sensitive to osmotic or pH stress, as well as to a variety of toxins. The results for M. extorquens DM4 are consistent with the hypothesis that hopanoids are important for membrane fluidity and lipid packing, but also indicate that the specific physiological processes that require hopanoids vary across bacterial lineages. Our work provides further support to emerging observations that the role of hopanoids in membrane robustness and barrier function may be important across lineages, possibly mediated through an interaction with lipid A in the outer membrane.

DNA Transposition at Work

DNA transposons are defined segments of DNA that are able to move from one genomic location to another. Movement is facilitated by one or more proteins, called the transposase, typically encoded by the mobile element itself. Here, we first provide an overview of the classification of such mobile elements in a variety of organisms. From a mechanistic perspective, we have focused on one particular group of DNA transposons that encode a transposase with a DD(E/D) catalytic domain that is topologically similar to RNase H. For these, a number of three-dimensional structures of transpososomes (transposase-nucleic acid complexes) are available, and we use these to describe the basics of their mechanisms. The DD(E/D) group, in addition to being the largest and most common among all DNA transposases, is the one whose members have been used for a wide variety of genomic applications. Therefore, a second focus of the article is to provide a nonexhaustive overview of transposon applications. Although several non-transposon-based approaches to site-directed genome modifications have emerged in the past decade, transposon-based applications are highly relevant when integration specificity is not sought. In fact, for many applications, the almost-perfect randomness and high frequency of integration make transposon-based approaches indispensable.

High-Throughput Identification of Adaptive Mutations in Experimentally Evolved Yeast Populations

High-throughput sequencing has enabled genetic screens that can rapidly identify mutations that occur during experimental evolution. The presence of a mutation in an evolved lineage does not, however, constitute proof that the mutation is adaptive, given the well-known and widespread phenomenon of genetic hitchhiking, in which a non-adaptive or even detrimental mutation can co-occur in a genome with a beneficial mutation and the combined genotype is carried to high frequency by selection. We approximated the spectrum of possible beneficial mutations in Saccharomyces cerevisiae using sets of single-gene deletions and amplifications of almost all the genes in the S. cerevisiae genome. We determined the fitness effects of each mutation in three different nutrient-limited conditions using pooled competitions followed by barcode sequencing. Although most of the mutations were neutral or deleterious, ~500 of them increased fitness. We then compared those results to the mutations that actually occurred during experimental evolution in the same three nutrient-limited conditions. On average, ~35% of the mutations that occurred during experimental evolution were predicted by the systematic screen to be beneficial. We found that the distribution of fitness effects depended on the selective conditions. In the phosphate-limited and glucose-limited conditions, a large number of beneficial mutations of nearly equivalent, small effects drove the fitness increases. In the sulfate-limited condition, one type of mutation, the amplification of the high-affinity sulfate transporter, dominated. In the absence of that mutation, evolution in the sulfate-limited condition involved mutations in other genes that were not observed previously—but were predicted by the systematic screen. Thus, gross functional screens have the potential to predict and identify adaptive mutations that occur during experimental evolution.

Experimental evolution allows us to observe evolution in real time. New advances in genome sequencing make it trivial to discover the mutations that have arisen in evolved cultures; however, linking those mutations to particular adaptive traits remains difficult. We evaluated the fitness impacts of thousands of single-gene losses and amplifications in yeast. We discovered that only a fraction of the hundreds of possible beneficial mutations were actually detected in evolution experiments performed previously. Our results provide evidence that 35% of the mutations identified in experimentally evolved populations are advantageous and that the distribution of beneficial fitness effects depends on the genetic background and the selective conditions. Furthermore, we show that it is possible to select for alternative mutations that improve fitness by blocking particularly high-fitness routes to adaptation.

Difference in C3-C4 metabolism underlies tradeoff between growth rate and biomass yield in Methylobacterium extorquens AM1

BACKGROUND

Two variants of Methylobacterium extorquens AM1 demonstrated a trade-off between growth rate and biomass yield. In addition, growth rate and biomass yield were also affected by supplementation of growth medium with different amounts of cobalt. The metabolism changes relating to these growth phenomena as well as the trade-off were investigated in this study. (13)C metabolic flux analysis was used to generate a detailed central carbon metabolic flux map with both absolute and normalized flux values.

RESULTS

The major differences between the two variants occurred at the formate node as well as within C3-C4 inter-conversion pathways. Higher relative fluxes through formyltetrahydrofolate ligase, phosphoenolpyruvate carboxylase, and malic enzyme led to higher biomass yield, while higher relative fluxes through pyruvate kinase and pyruvate dehydrogenase led to higher growth rate. These results were then tested by phenotypic studies on three mutants (null pyk, null pck mutant and null dme mutant) in both variants, which agreed with the model prediction.

CONCLUSIONS

In this study, (13)C metabolic flux analysis for two strain variants of M. extorquens AM1 successfully identified metabolic pathways contributing to the trade-off between cell growth and biomass yield. Phenotypic analysis of mutants deficient in corresponding genes supported the conclusion that C3-C4 inter-conversion strategies were the major response to the trade-off.

Selection Maintains Apparently Degenerate Metabolic Pathways due to Tradeoffs in Using Methylamine for Carbon versus Nitrogen

Microorganisms often encode multiple non-orthologous metabolic modules that catalyze the same reaction. However, little experimental evidence actually demonstrates a selective basis for metabolic degeneracy. Many methylotrophs-microorganisms that grow on reduced single-carbon compounds-like Methylobacterium extorquens AM1 encode two routes for methylamine oxidation: the periplasmic methylamine dehydrogenase (MaDH) and the cytoplasmic N-methylglutamate (NMG) pathway. In Methylobacterium extorquens AM1, MaDH is essential for methylamine growth, but the NMG pathway has no known physiological role. Here, we use experimental evolution of two isolates lacking (or incapable of using) MaDH to uncover the physiological challenges that need to be overcome in order to use the NMG pathway for growth on methylamine as a carbon and energy source. Physiological characterization of the evolved isolates revealed regulatory rewiring to increase expression of the NMG pathway and novel mechanisms to mitigate cytoplasmic ammonia buildup. These adaptations led us to infer and validate environmental conditions under which the NMG pathway is advantageous compared to MaDH. The highly expressed MaDH enables rapid growth on high concentrations of methylamine as the primary carbon and energy substrate, whereas the energetically expensive NMG pathway plays a pivotal role during growth with methylamine as the sole nitrogen source, which we demonstrate is especially true under limiting concentrations (<1 mM). Tradeoffs between cellular localization and ammonium toxicity lead to selection for this apparent degeneracy as it is beneficial to facultative methylotrophs that have to switch between using methylamine as a carbon and energy source or just a nitrogen source.

Adaptation to Parasites and Costs of Parasite Resistance in Mutator and Nonmutator Bacteria

Parasitism creates selection for resistance mechanisms in host populations and is hypothesized to promote increased host evolvability. However, the influence of these traits on host evolution when parasites are no longer present is unclear. We used experimental evolution and whole-genome sequencing of Escherichia coli to determine the effects of past and present exposure to parasitic viruses (phages) on the spread of mutator alleles, resistance, and bacterial competitive fitness. We found that mutator alleles spread rapidly during adaptation to any of four different phage species, and this pattern was even more pronounced with multiple phages present simultaneously. However, hypermutability did not detectably accelerate adaptation in the absence of phages and recovery of fitness costs associated with resistance. Several lineages evolved phage resistance through elevated mucoidy, and during subsequent evolution in phage-free conditions they rapidly reverted to nonmucoid, phage-susceptible phenotypes. Genome sequencing revealed that this phenotypic reversion was achieved by additional genetic changes rather than by genotypic reversion of the initial resistance mutations. Insertion sequence (IS) elements played a key role in both the acquisition of resistance and adaptation in the absence of parasites; unlike single nucleotide polymorphisms, IS insertions were not more frequent in mutator lineages. Our results provide a genetic explanation for rapid reversion of mucoidy, a phenotype observed in other bacterial species including human pathogens. Moreover, this demonstrates that the types of genetic change underlying adaptation to fitness costs, and consequently the impact of evolvability mechanisms such as increased point-mutation rates, depend critically on the mechanism of resistance.

Buried Treasure: Evolutionary Perspectives on Microbial Iron Piracy

Host–pathogen interactions provide valuable systems for the study of evolutionary genetics and natural selection. The sequestration of essential iron has emerged as a crucial innate defense system termed nutritional immunity, leading pathogens to evolve mechanisms of ‘iron piracy’ to scavenge this metal from host proteins. This battle for iron carries numerous consequences not only for host–pathogen evolution but also microbial community interactions. Here we highlight recent and potential future areas of investigation on the evolutionary implications of microbial iron piracy in relation to molecular arms races, host range, competition, and virulence. Applying evolutionary genetic approaches to the study of microbial iron acquisition could also provide new inroads for understanding and combating infectious disease.

The battle between microbes and their hosts for nutrient iron is emerging as a new front of evolutionary genetic conflict.

Molecular arms races can emerge between host iron-binding proteins and microbial ‘iron piracy’ factors that steal this nutrient for growth. Such rapid evolution may also contribute to the host range of pathogenic microbes.

Iron acquisition plays an important role in evolutionary interactions between microbes, both in the environment and within the host. Competition for iron can prevent infection by pathogens, while genetic changes in iron acquisition systems can enhance microbial virulence.

Evolutionary conflicts for nutrient iron are revealing potential new genetic mechanisms of disease resistance as well as avenues for therapeutic development.

ISQuest: finding insertion sequences in prokaryotic sequence fragment data

MOTIVATION

Insertion sequences (ISs) are transposable elements present in most bacterial and archaeal genomes that play an important role in genomic evolution. The increasing availability of sequenced prokaryotic genomes offers the opportunity to study ISs comprehensively, but development of efficient and accurate tools is required for discovery and annotation. Additionally, prokaryotic genomes are frequently deposited as incomplete, or draft stage because of the substantial cost and effort required to finish genome assembly projects. Development of methods to identify IS directly from raw sequence reads or draft genomes are therefore desirable. Software tools such as Optimized Annotation System for Insertion Sequences and IScan currently identify IS elements in completely assembled and annotated genomes; however, to our knowledge no methods have been developed to identify ISs from raw fragment data or partially assembled genomes. We have developed novel methods to solve this computationally challenging problem, and implemented these methods in the software package ISQuest. This software identifies bacterial ISs and their sequence elements-inverted and direct repeats-in raw read data or contigs using flexible search parameters. ISQuest is capable of finding ISs in hundreds of partially assembled genomes within hours, making it a valuable high-throughput tool for a global search of IS elements. We tested ISQuest on simulated read libraries of 3810 complete bacterial genomes and plasmids in GenBank and were capable of detecting 82% of the ISs and transposases annotated in GenBank with 80% sequence identity.

CONTACT

abiswas@cs.odu.edu.

Efficient reduction of the formation of by-products and improvement of production yield of 2,3-butanediol by a combined deletion of alcohol dehydrogenase, acetate kinase-phosphotransacetylase, and lactate dehydrogenase genes in metabolically engineered Klebsiella oxytoca in mineral salts medium

Klebsiella oxytoca KMS005 (∆adhE∆ackA-pta∆ldhA) was metabolically engineered to improve 2,3-butanediol (BDO) yield. Elimination of alcohol dehydrogenase E (adhE), acetate kinase A-phosphotransacetylase (ackA-pta), and lactate dehydrogenase A (ldhA) enzymes allowed BDO production as a primary pathway for NADH re-oxidation, and significantly reduced by-products. KMS005 was screened for the efficient glucose utilization by metabolic evolution. KMS005-73T improved BDO production at a concentration of 23.5±0.5 g/L with yield of 0.46±0.02 g/g in mineral salts medium containing 50 g/L glucose in a shake flask. KMS005-73T also exhibited BDO yields of about 0.40-0.42 g/g from sugarcane molasses, cassava starch, and maltodextrin. During fed-batch fermentation, KMS005-73T produced BDO at a concentration, yield, and overall and specific productivities of 117.4±4.5 g/L, 0.49±0.02 g/g, 1.20±0.05 g/Lh, and 27.2±1.1 g/gCDW, respectively. No acetoin, lactate, and formate were detected, and only trace amounts of acetate and ethanol were formed. The strain also produced the least by-products and the highest BDO yield among other Klebsiella strains previously developed.

Parallel and Divergent Evolutionary Solutions for the Optimization of an Engineered Central Metabolism in Methylobacterium extorquens AM1

Bioengineering holds great promise to provide fast and efficient biocatalysts for methanol-based biotechnology, but necessitates proven methods to optimize physiology in engineered strains. Here, we highlight experimental evolution as an effective means for optimizing an engineered Methylobacterium extorquens AM1. Replacement of the native formaldehyde oxidation pathway with a functional analog substantially decreased growth in an engineered Methylobacterium, but growth rapidly recovered after six hundred generations of evolution on methanol. We used whole-genome sequencing to identify the basis of adaptation in eight replicate evolved strains, and examined genomic changes in light of other growth and physiological data. We observed great variety in the numbers and types of mutations that occurred, including instances of parallel mutations at targets that may have been "rationalized" by the bioengineer, plus other "illogical" mutations that demonstrate the ability of evolution to expose unforeseen optimization solutions. Notably, we investigated mutations to RNA polymerase, which provided a massive growth benefit but are linked to highly aberrant transcriptional profiles. Overall, we highlight the power of experimental evolution to present genetic and physiological solutions for strain optimization, particularly in systems where the challenges of engineering are too many or too difficult to overcome via traditional engineering methods.

Transhydrogenase promotes the robustness and evolvability of E. coli deficient in NADPH production

Metabolic networks revolve around few metabolites recognized by diverse enzymes and involved in myriad reactions. Though hub metabolites are considered as stepping stones to facilitate the evolutionary expansion of biochemical pathways, changes in their production or consumption often impair cellular physiology through their system-wide connections. How does metabolism endure perturbations brought immediately by pathway modification and restore hub homeostasis in the long run? To address this question we studied laboratory evolution of pathway-engineered Escherichia coli that underproduces the redox cofactor NADPH on glucose. Literature suggests multiple possibilities to restore NADPH homeostasis. Surprisingly, genetic dissection of isolates from our twelve evolved populations revealed merely two solutions: (1) modulating the expression of membrane-bound transhydrogenase (mTH) in every population; (2) simultaneously consuming glucose with acetate, an unfavored byproduct normally excreted during glucose catabolism, in two subpopulations. Notably, mTH displays broad phylogenetic distribution and has also played a predominant role in laboratory evolution of Methylobacterium extorquens deficient in NADPH production. Convergent evolution of two phylogenetically and metabolically distinct species suggests mTH as a conserved buffering mechanism that promotes the robustness and evolvability of metabolism. Moreover, adaptive diversification via evolving dual substrate consumption highlights the flexibility of physiological systems to exploit ecological opportunities.

Methylobacterium extorquens: methylotrophy and biotechnological applications

Methylotrophy is the ability to use reduced one-carbon compounds, such as methanol, as a single source of carbon and energy. Methanol is, due to its availability and potential for production from renewable resources, a valuable feedstock for biotechnology. Nature offers a variety of methylotrophic microorganisms that differ in their metabolism and represent resources for engineering of value-added products from methanol. The most extensively studied methylotroph is the Alphaproteobacterium Methylobacterium extorquens. Over the past five decades, the metabolism of M. extorquens has been investigated physiologically, biochemically, and more recently, using complementary omics technologies such as transcriptomics, proteomics, metabolomics, and fluxomics. These approaches, together with a genome-scale metabolic model, facilitate system-wide studies and the development of rational strategies for the successful generation of desired products from methanol. This review summarizes the knowledge of methylotrophy in M. extorquens, as well as the available tools and biotechnological applications.

Genetic and phenotypic comparison of facultative methylotrophy between Methylobacterium extorquens strains PA1 and AM1

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 C₁ 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 C₁ 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.

Genome information of Methylobacterium oryzae, a plant-probiotic methylotroph in the phyllosphere

Pink-pigmented facultative methylotrophs in the Rhizobiales are widespread in the environment, and many Methylobacterium species associated with plants produce plant growth-promoting substances. To gain insights into the life style at the phyllosphere and the genetic bases of plant growth promotion, we determined and analyzed the complete genome sequence of Methylobacterium oryzae CBMB20T, a strain isolated from rice stem. The genome consists of a 6.29-Mb chromosome and four plasmids, designated as pMOC1 to pMOC4. Among the 6,274 coding sequences in the chromosome, the bacterium has, besides most of the genes for the central metabolism, all of the essential genes for the assimilation and dissimilation of methanol that are either located in methylotrophy islands or dispersed. M. oryzae is equipped with several kinds of genes for adaptation to plant surfaces such as defense against UV radiation, oxidative stress, desiccation, or nutrient deficiency, as well as high proportion of genes related to motility and signaling. Moreover, it has an array of genes involved in metabolic pathways that may contribute to promotion of plant growth; they include auxin biosynthesis, cytokine biosynthesis, vitamin B12 biosynthesis, urea metabolism, biosorption of heavy metals or decrease of metal toxicity, pyrroloquinoline quinone biosynthesis, 1-aminocyclopropane-1-carboxylate deamination, phosphate solubilization, and thiosulfate oxidation. Through the genome analysis of M. oryzae, we provide information on the full gene complement of M. oryzae that resides in the aerial parts of plants and enhances plant growth. The plant-associated lifestyle of M. oryzae pertaining to methylotrophy and plant growth promotion, and its potential as a candidate for a bioinoculant targeted to the phyllosphere and focused on phytostimulation are illuminated.

Microbial evolution. Global epistasis makes adaptation predictable despite sequence-level stochasticity

Epistatic interactions between mutations can make evolutionary trajectories contingent on the chance occurrence of initial mutations. We used experimental evolution in Saccharomyces cerevisiae to quantify this contingency, finding differences in adaptability among 64 closely related genotypes. Despite these differences, sequencing of 104 evolved clones showed that initial genotype did not constrain future mutational trajectories. Instead, reconstructed combinations of mutations revealed a pattern of diminishing-returns epistasis: Beneficial mutations have consistently smaller effects in fitter backgrounds. Taken together, these results show that beneficial mutations affecting a variety of biological processes are globally coupled; they interact strongly, but only through their combined effect on fitness. As a consequence, fitness evolution follows a predictable trajectory even though sequence-level adaptation is stochastic.

Mapping the fitness landscape of gene expression uncovers the cause of antagonism and sign epistasis between adaptive mutations

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.

The pace and outcome of a series of adaptive steps in an evolving lineage depends upon how well different beneficial mutations stack on top of each other. We found that independent beneficial mutations that affected gene expression for a metabolic pathway did not work well together, and were often jointly deleterious. The most beneficial mutations interacted the most poorly with others, which was a trend we found common in other biological systems. Through generating a model that accounted for enzymatic benefits and expression costs, we uncovered that this antagonism was caused by a phenotype to fitness mapping that had an intermediate peak. This allowed us to predict the fitness effect of double mutants and to uncover that the single winning mutations tended to move straight to the peak in a single step. These findings demonstrate the importance of considering the phenotypic changes that cause nonlinear interactions between mutations upon fitness, and thus influence how populations evolve.

Metabolic erosion primarily through mutation accumulation, and not tradeoffs, drives limited evolution of substrate specificity in Escherichia coli

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.

Laboratory divergence of Methylobacterium extorquens AM1 through unintended domestication and past selection for antibiotic resistance

Background

A common assumption of microorganisms is that laboratory stocks will remain genetically and phenotypically constant over time, and across laboratories. It is becoming increasingly clear, however, that mutations can ruin strain integrity and drive the divergence or “domestication” of stocks. Since its discovery in 1960, a stock of Methylobacterium extorquens AM1 (“AM1”) has remained in the lab, propagated across numerous growth and storage conditions, researchers, and facilities. To explore the extent to which this lineage has diverged, we compared our own “Modern” stock of AM1 to a sample archived at a culture stock center shortly after the strain’s discovery. Stored as a lyophilized sample, we hypothesized that this Archival strain would better reflect the first-ever isolate of AM1 and reveal ways in which our Modern stock has changed through laboratory domestication or other means.

Results

Using whole-genome re-sequencing, we identified some 29 mutations – including single nucleotide polymorphisms, small indels, the insertion of mobile elements, and the loss of roughly 36 kb of DNA - that arose in the laboratory-maintained Modern lineage. Contrary to our expectations, Modern was both slower and less fit than Archival across a variety of growth substrates, and showed no improvement during long-term growth and storage. Modern did, however, outperform Archival during growth on nutrient broth, and in resistance to rifamycin, which was selected for by researchers in the 1980s. Recapitulating selection for rifamycin resistance in replicate Archival populations showed that mutations to RNA polymerase B (rpoB) substantially decrease growth in the absence of antibiotic, offering an explanation for slower growth in Modern stocks. Given the large number of genomic changes arising from domestication (28), it is somewhat surprising that the single other mutation attributed to purposeful laboratory selection accounts for much of the phenotypic divergence between strains.

Conclusions

These results highlight the surprising degree to which AM1 has diverged through a combination of unintended laboratory domestication and purposeful selection for rifamycin resistance. Instances of strain divergence are important, not only to ensure consistency of experimental results, but also to explore how microbes in the lab diverge from one another and from their wild counterparts.

Metagenomic analysis of size-fractionated picoplankton in a marine oxygen minimum zone

Marine oxygen minimum zones (OMZs) support diverse microbial communities with roles in major elemental cycles. It is unclear how the taxonomic composition and metabolism of OMZ microorganisms vary between particle-associated and free-living size fractions. We used amplicon (16S rRNA gene) and shotgun metagenome sequencing to compare microbial communities from large (>1.6 μm) and small (0.2-1.6 μm) filter size fractions along a depth gradient in the OMZ off Chile. Despite steep vertical redox gradients, size fraction was a significantly stronger predictor of community composition compared to depth. Phylogenetic diversity showed contrasting patterns, decreasing towards the anoxic OMZ core in the small size fraction, but exhibiting maximal values at these depths within the larger size fraction. Fraction-specific distributions were evident for key OMZ taxa, including anammox planctomycetes, whose coding sequences were enriched up to threefold in the 0.2-1.6 μm community. Functional gene composition also differed between fractions, with the >1.6 μm community significantly enriched in genes mediating social interactions, including motility, adhesion, cell-to-cell transfer, antibiotic resistance and mobile element activity. Prokaryotic transposase genes were three to six fold more abundant in this fraction, comprising up to 2% of protein-coding sequences, suggesting that particle surfaces may act as hotbeds for transposition-based genome changes in marine microbes. Genes for nitric and nitrous oxide reduction were also more abundant (three to seven fold) in the larger size fraction, suggesting microniche partitioning of key denitrification steps. These results highlight an important role for surface attachment in shaping community metabolic potential and genome content in OMZ microorganisms.

Environment determines epistatic patterns for a ssDNA virus

Despite the accumulation of substantial quantities of information about epistatic interactions among both deleterious and beneficial mutations in a wide array of experimental systems, neither consistent patterns nor causal explanations for these interactions have yet emerged. Furthermore, the effects of mutations depend on the environment in which they are characterized, implying that the environment may also influence epistatic interactions. Recent work with beneficial mutations for the single-stranded DNA bacteriophage ID11 demonstrated that interactions between pairs of mutations could be understood by means of a simple model that assumes that mutations have additive phenotypic effects and that epistasis arises through a nonlinear phenotype-fitness map with a single intermediate optimum. To determine whether such a model could also explain changes in epistatic patterns associated with changes in environment, we measured epistatic interactions for these same mutations under conditions for which we expected to find the wild-type ID11 at different distances from its phenotypic optimum by assaying fitnesses at three different temperatures: 33°, 37°, and 41°. Epistasis was present and negative under all conditions, but became more pronounced as temperature increased. We found that the additive-phenotypes model explained these patterns as changes in the parameters of the phenotype-fitness map, but that a model that additionally allows the phenotypes to vary across temperatures performed significantly better. Our results show that ostensibly complex patterns of fitness effects and epistasis across environments can be explained by assuming a simple structure for the genotype-phenotype relationship.

A novel pair of inducible expression vectors for use in Methylobacterium extorquens

BACKGROUND

Due to the ever increasing use of diverse microbial taxa in basic research and industrial settings, there is a growing need for genetic tools to alter the physiology of these organisms. In particular, there is a dearth of inducible expression systems available for bacteria outside commonly used γ-proteobacteria, such as Escherichia coli or Pseudomonas species. To this end, we have sought to develop a pair of inducible expression vectors for use in the α-proteobacterium Methylobacterium extorquens, a model methylotroph.

FINDINGS

We found that the P(R) promoter from rhizobial phage 16-3 was active in M. extorquens and engineered the promoter to be inducible by either p-isopropyl benzoate (cumate) or anhydrotetracycline. These hybrid promoters, P(R/cmtO) and P(R/tetO), were found to have high levels of expression in M. extorquens with a regulatory range of 10-fold and 30-fold, respectively. Compared to an existing cumate-inducible (10-fold range), high-level expression system for M. extorquens, P(R/cmtO) and P(R/tetO) have 33% of the maximal activity but were able to repress gene expression 3 and 8-fold greater, respectively. Both promoters were observed to exhibit homogeneous, titratable activation dynamics rather than on-off, switch-like behavior. The utility of these promoters was further demonstrated by complementing loss of function of ftfL--essential for growth on methanol--where we show P(R/tetO) is capable of not only fully complementing function but also producing a conditional null phenotype. These promoters have been incorporated into a broad-host-range backbone allowing for potential use in a variety of bacterial hosts.

CONCLUSIONS

We have developed two novel expression systems for use in M. extorquens. The expression range of these vectors should allow for increased ability to explore cellular physiology in M. extorquens. Further, the P(R/tetO) promoter is capable of producing conditional null phenotypes, previously unattainable in M. extorquens. As both expression systems rely on the use of membrane permeable inducers, we suspect these expression vectors will be useful for ectopic gene expression in numerous proteobacteria.

Development of an optimized medium, strain and high-throughput culturing methods for Methylobacterium extorquens

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 MethylobacteriumPIPES (‘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.

The 380 kb pCMU01 plasmid encodes chloromethane utilization genes and redundant genes for vitamin B12- and tetrahydrofolate-dependent chloromethane metabolism in Methylobacterium extorquens CM4: a proteomic and bioinformatics study

Chloromethane (CH3Cl) is the most abundant volatile halocarbon in the atmosphere and contributes to the destruction of stratospheric ozone. The only known pathway for bacterial chloromethane utilization (cmu) was characterized in Methylobacterium extorquens CM4, a methylotrophic bacterium able to utilize compounds without carbon-carbon bonds such as methanol and chloromethane as the sole carbon source for growth. Previous work demonstrated that tetrahydrofolate and vitamin B12 are essential cofactors of cmuA- and cmuB-encoded methyltransferases of chloromethane dehalogenase, and that the pathway for chloromethane utilization is distinct from that for methanol. This work reports genomic and proteomic data demonstrating that cognate cmu genes are located on the 380 kb pCMU01 plasmid, which drives the previously defined pathway for tetrahydrofolate-mediated chloromethane dehalogenation. Comparison of complete genome sequences of strain CM4 and that of four other M. extorquens strains unable to grow with chloromethane showed that plasmid pCMU01 harbors unique genes without homologs in the compared genomes (bluB2, btuB, cobA, cbiD), as well as 13 duplicated genes with homologs of chromosome-borne genes involved in vitamin B12-associated biosynthesis and transport, or in tetrahydrofolate-dependent metabolism (folC2). In addition, the presence of both chromosomal and plasmid-borne genes for corrinoid salvaging pathways may ensure corrinoid coenzyme supply in challenging environments. Proteomes of M. extorquens CM4 grown with one-carbon substrates chloromethane and methanol were compared. Of the 49 proteins with differential abundance identified, only five (CmuA, CmuB, PurU, CobH2 and a PaaE-like uncharacterized putative oxidoreductase) are encoded by the pCMU01 plasmid. The mainly chromosome-encoded response to chloromethane involves gene clusters associated with oxidative stress, production of reducing equivalents (PntAA, Nuo complex), conversion of tetrahydrofolate-bound one-carbon units, and central metabolism. The mosaic organization of plasmid pCMU01 and the clustering of genes coding for dehalogenase enzymes and for biosynthesis of associated cofactors suggests a history of gene acquisition related to chloromethane utilization.

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.

Acclimatizing and adaptive (evolutionary) processes allow organisms to thrive despite cellular and environmental perturbations. Our work examined whether adaptation restores stress responses towards wild-type (pre-stressed) versus novel physiological states during adaptation by studying a bacterium (Methylobacterium extorquens AM1) that was experimentally engineered and evolved with a novel central metabolism. The engineered strain was much slower and less fit than wild-type, but eight replicate populations evolved for six hundred generations showed substantial improvements. We found that changes in gene expression during adaptation consistently restored acclimatizing processes to the wild-type state, often in 8/8 evolved lines. Novel changes were common and largely restricted to one lineage; however, highly parallel novel changes revealed loci harboring beneficial mutations. Even rarer were reinforced changes, such as pntAB transhydrogenase, which increased beyond immediate acclimation during evolution to restore NAD(P)(H) metabolism and improve growth. Overall, a few novel or reinforcing changes drove the mass-restoration of physiology back to wild-type.

Can you sequence ecology? Metagenomics of adaptive diversification

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.

Synchronous waves of failed soft sweeps in the laboratory: remarkably rampant clonal interference of alleles at a single locus

It has increasingly been recognized that adapting populations of microbes contain not one, but many lineages continually arising and competing at once. This process, termed "clonal interference," alters the rate and dynamics of adaptation and biases winning mutations toward those with the largest selective effect. Here we uncovered a dramatic example of clonal interference between multiple similar mutations occurring at the same locus within replicate populations of Methylobacterium extorquens AM1. Because these mutational events involved the transposition of an insertion sequence into a narrow window of a single gene, they were both readily detectable at low frequencies and could be distinguished due to differences in insertion sites. This allowed us to detect up to 17 beneficial alleles of this type coexisting in a single population. Despite conferring a large selective benefit, the majority of these alleles rose and then fell in frequency due to other lineages emerging that were more fit. By comparing allele-frequency dynamics to the trajectories of fitness gains by these populations, we estimated the fitness values of the genotypes that contained these mutations. Collectively across all populations, these alleles arose upon backgrounds with a wide range of fitness values. Within any single population, however, multiple alleles tended to rise and fall synchronously during a single wave of multiple genotypes with nearly identical fitness values. These results suggest that alleles of large benefit arose repeatedly in failed "soft sweeps" during narrow windows of adaptation due to the combined effects of epistasis and clonal interference.

The basis of antagonistic pleiotropy in hfq mutations that have opposite effects on fitness at slow and fast growth rates

Mutations beneficial in one environment may cause costs in different environments, resulting in antagonistic pleiotropy. Here, we describe a novel form of antagonistic pleiotropy that operates even within the same environment, where benefits and deleterious effects exhibit themselves at different growth rates. The fitness of hfq mutations in Escherichia coli affecting the RNA chaperone involved in small-RNA regulation is remarkably sensitive to growth rate. E. coli populations evolving in chemostats under nutrient limitation acquired beneficial mutations in hfq during slow growth (0.1 h(-1)) but not in populations growing sixfold faster. Four identified hfq alleles from parallel populations were beneficial at 0.1 h(-1) and deleterious at 0.6 h(-1). The hfq mutations were beneficial, deleterious or neutral at an intermediate growth rate (0.5 h(-1)) and one changed from beneficial to deleterious within a 36 min difference in doubling time. The benefit of hfq mutations was due to the greater transport of limiting nutrient, which diminished at higher growth rates. The deleterious effects of hfq mutations at 0.6 h(-1) were less clear, with decreased viability a contributing factor. The results demonstrate distinct pleiotropy characteristics in the alleles of the same gene, probably because the altered residues in Hfq affected the regulation of expression of different genes in distinct ways. In addition, these results point to a source of variation in experimental measurement of the selective advantage of a mutation; estimates of fitness need to consider variation in growth rate impacting on the magnitude of the benefit of mutations and on their fitness distributions.