Search

ABSTRACT Methylothon is an inquiry-based high school learning module in microbial ecology, molecular biology, and bioinformatics that centers around pink-pigmented plant-associated methylotrophic bacteria. Here, we present an overview of the module’s learning goals, describe course resources (available for public use at http://methylothon.com), and relate lessons learned from adapting Methylothon for remote learning during the pandemic in spring of 2021. This curriculum description is intended not only for instructors but also for microbial ecology researchers with an interest in conducting K-12 outreach. The original in-person version of the module allows students to isolate their own strains of methylotrophic bacteria from plants they sample from the environment, to identify these using PCR, sequencing, and phylogenetic analysis, and to contribute their strains to original research in a university lab. The adapted version strengthens the focus on bioinformatics and increases its flexibility and accessibility by making the lab portion optional and adopting free web-based tools. Student feedback and graded assignments from spring 2021 revealed that the lesson was especially effective at introducing the concepts of BLAST and phylogenetic trees and that students valued and felt inspired by the opportunity to conduct hands-on work and to participate in community science.

Combination antimicrobial therapy has been considered a promising strategy to combat the evolution of antimicrobial resistance. Francisella tularensis is the causative agent of tularemia and in addition to being found in the nature, is recognized as a threat agent that requires vigilance. We investigated the evolutionary outcome of adapting the Live Vaccine Strain (LVS) of F. tularensis subsp. holarctica to two non-interacting drugs, ciprofloxacin and doxycycline, individually, sequentially, and in combination. Despite their individual efficacies and independence of mechanisms, evolution to the combination arose on a shorter time scale than evolution to the two drugs sequentially. We conducted a longitudinal mutational analysis of the populations evolving to the drug combination, genetically reconstructed the identified evolutionary pathway, and carried out biochemical validation. We discovered that, after the appearance of an initial weak generalist mutation (FupA/B), each successive mutation alternated between adaptation to one drug or the other. In combination, these mutations allowed the population to more efficiently ascend the fitness peak through a series of evolutionary switch-backs. Clonal interference, weak pleiotropy, and positive epistasis also contributed to combinatorial evolution. This finding suggests that the use of this non-interacting drug pair against F. tularensis may render both drugs ineffective because of mutational switch-backs that accelerate evolution of dual resistance.

Pink-pigmented facultative methylotrophs have long been studied for their ability to grow on reduced single-carbon (C1) compounds. The C1 groups that support methylotrophic growth may come from a variety of sources. Here, we describe a group of Methylobacterium strains that can engage in methoxydotrophy: they can metabolize the methoxy groups from several aromatic compounds that are commonly the product of lignin depolymerization. Furthermore, these organisms can utilize the full aromatic ring as a growth substrate, a phenotype that has rarely been described in Methylobacterium. We demonstrated growth on p-hydroxybenzoate, protocatechuate, vanillate, and ferulate in laboratory culture conditions. We also used comparative genomics to explore the evolutionary history of this trait, finding that the capacity for aromatic catabolism is likely ancestral to two clades of Methylobacterium, but has also been acquired horizontally by closely related organisms. In addition, we surveyed the published metagenome data to find that the most abundant group of aromatic-degrading Methylobacterium in the environment is likely the group related to Methylobacterium nodulans, and they are especially common in soil and root environments. The demethoxylation of lignin-derived aromatic monomers in aerobic environments releases formaldehyde, a metabolite that is a potent cellular toxin but that is also a growth substrate for methylotrophs. We found that, whereas some known lignin-degrading organisms excrete formaldehyde as a byproduct during growth on vanillate, Methylobacterium do not. This observation is especially relevant to our understanding of the ecology and the bioengineering of lignin degradation.

ABSTRACT Methylobacterium is a prevalent bacterial genus of the phyllosphere. Despite its ubiquity, little is known about the extent to which its diversity reflects neutral processes like migration and drift, versus environmental filtering of life history strategies and adaptations. In two temperate forests, we investigated how phylogenetic diversity within Methylobacterium is structured by biogeography, seasonality, and growth strategies. Using deep, culture-independent barcoded marker gene sequencing coupled with culture-based approaches, we uncovered a considerable diversity of Methylobacterium in the phyllosphere. We cultured different subsets of Methylobacterium lineages depending upon the temperature of isolation and growth (20°C or 30°C), suggesting long-term adaptation to temperature. To a lesser extent than temperature adaptation, Methylobacterium diversity was also structured across large (>100 km; between forests) and small (

The recalcitrance of complex organic polymers such as lignocellulose is one of the major obstacles to sustainable energy production from plant biomass, and the generation of toxic intermediates can negatively impact the efficiency of microbial lignocellulose degradation. Here, we describe the development of a model microbial consortium for studying lignocellulose degradation, with the specific goal of mitigating the production of the toxin formaldehyde during the breakdown of methoxylated aromatic compounds. Included are Pseudomonas putida, a lignin degrader; Cellulomonas fimi, a cellulose degrader; and sometimes Yarrowia lipolytica, an oleaginous yeast. Unique to our system is the inclusion of Methylorubrum extorquens, a methylotroph capable of using formaldehyde for growth. We developed a defined minimal “Model Lignocellulose” growth medium for reproducible coculture experiments. We demonstrated that the formaldehyde produced by P. putida growing on vanillic acid can exceed the minimum inhibitory concentration for C. fimi, and, furthermore, that the presence of M. extorquens lowers those concentrations. We also uncovered unexpected ecological dynamics, including resource competition, and interspecies differences in growth requirements and toxin sensitivities. Finally, we introduced the possibility for a mutualistic interaction between C. fimi and M. extorquens through metabolite exchange. This study lays the foundation to enable future work incorporating metabolomic analysis and modeling, genetic engineering, and laboratory evolution, on a model system that is appropriate both for fundamental eco-evolutionary studies and for the optimization of efficiency and yield in microbially-mediated biomass transformation.

The potency and indiscriminate nature of formaldehyde reactivity upon biological molecules make it a universal stressor. However, some organisms such as Methylorubrum extorquens possess means to rapidly and effectively mitigate formaldehyde-induced damage. EfgA is a recently identified formaldehyde sensor predicted to halt translation in response to elevated formaldehyde as a means to protect cells. Herein, we investigate growth and changes in gene expression to understand how M. extorquens responds to formaldehyde with and without the EfgA-formaldehyde-mediated translational response, and how this mechanism compares to antibiotic-mediated translation inhibition. These distinct mechanisms of translation inhibition have notable differences: they each involve different specific players and in addition, formaldehyde also acts as a general, multi-target stressor and a potential carbon source. We present findings demonstrating that in addition to its characterized impact on translation, functional EfgA allows for a rapid and robust transcriptional response to formaldehyde and that removal of EfgA leads to heightened proteotoxic and genotoxic stress in the presence of increased formaldehyde levels. We also found that many downstream consequences of translation inhibition were shared by EfgA-formaldehyde- and kanamycin-mediated translation inhibition. Our work uncovered additional layers of regulatory control enacted by functional EfgA upon experiencing formaldehyde stress, and further demonstrated the importance this protein plays at both transcriptional and translational levels in this model methylotroph.

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.

Methylamine plays an important role in the global carbon and nitrogen budget; microorganisms that grow on reduced single carbon compounds, methylotrophs, serve as a major biological sink for methylamine in aerobic environments. Two non-orthologous, functionally degenerate routes for methylamine oxidation have been studied in methylotrophic Proteobacteria: Methylamine dehydrogenase and the N-methylglutamate pathway. Recent work suggests the N-methylglutamate (NMG) pathway may be more common in nature than the well-studied methylamine dehydrogenase (MaDH, encoded by the mau gene cluster). However, the distribution of these pathways across methylotrophs has never been analyzed. Furthermore, even though horizontal gene transfer (HGT) is commonly invoked as a means to transfer these pathways between strains, the physiological barriers to doing so have not been investigated. We found that the NMG pathway is both more abundant and more universally distributed across methylotrophic Proteobacteria compared to MaDH, which displays a patchy distribution and has clearly been transmitted by HGT even amongst very closely related strains. This trend was especially prominent in well-characterized strains of the Methylobacterium extroquens species, which also display significant phenotypic variability during methylamine growth. Strains like Methylobacterium extorquens PA1 that only encode the NMG pathway grew on methylamine at least five-fold slower than strains like Methylobacterium extorquens AM1 that also possess the mau gene cluster. By mimicking a HGT event through the introduction of the M. extorquens AM1 mau gene cluster into the PA1 genome, the resulting strain instantaneously achieved a 4.5-fold increase in growth rate on methylamine and a 11-fold increase in fitness on methylamine, which even surpassed the fitness of M. extorquens AM1. In contrast, when three replicate populations of wild type M. extorquens PA1 were evolved on methylamine as the sole carbon and energy source for 150 generations neither fitness nor growth rate improved. These results suggest that the NMG pathway permits slow growth on methylamine and is widely distributed in methylotrophs; however, rapid growth on methylamine can be achieved quite readily through acquisition of the mau cluster by HGT.

The human population is growing and, globally, we must meet the challenge of increased protein needs required to feed this population. Single cell proteins (SCP), when coupled to aquaculture production, offer a means to ensure future protein needs can be met without direct competition with food for people. To demonstrate a given type of SCP has potential as a protein source for use in aquaculture feed, a number of steps need to be validated including demonstrating that the SCP is accepted by the species in question, leads to equivalent survival and growth, does not result in illness or other maladies, is palatable to the consumer, is cost effective to produce and can easily be incorporated into diets using existing technology. Here we examine white shrimp (Litopenaeus vannamei) growth and consumer taste preference, smallmouth grunt (Haemulon chrysargyreum) growth, survival, health and gut microbiota, and Atlantic salmon (Salmo salar) digestibility when fed diets that substitute the bacterium Methylobacterium extorquens at a level of 30% (grunts), 100% (shrimp), or 55% (salmon) of the fishmeal in a compound feed. In each of these tests, animals performed equivalently when fed diets containing M. extorquens as when fed a standard aquaculture diet. This transdisciplinary approach is a first validation of this bacterium as a potential SCP protein substitute in aquafeeds. Given the ease to produce this SCP through an aerobic fermentation process, the broad applicability for use in aquaculture indicates the promise of M. extorquens in leading toward greater food security in the future.

The interspecies exchange of metabolites plays a key role in the spatiotemporal dynamics of microbial communities. This raises the question of whether ecosystem-level behavior of structured communities can be predicted using genome-scale metabolic models for multiple organisms. We developed a modeling framework that integrates dynamic flux balance analysis with diffusion on a lattice and applied it to engineered communities. First, we predicted and experimentally confirmed the species ratio to which a two-species mutualistic consortium converges and the equilibrium composition of a newly engineered three-member community. We next identified a specific spatial arrangement of colonies, which gives rise to what we term the “eclipse dilemma”: does a competitor placed between a colony and its cross-feeding partner benefit or hurt growth of the original colony? Our experimentally validated finding that the net outcome is beneficial highlights the complex nature of metabolic interactions in microbial communities while at the same time demonstrating their predictability. Video Abstract:

Chloromethane is an ozone-depleting gas, produced predominantly from natural sources, that provides an important carbon source for microbes capable of consuming it. Chloromethane catabolism has been difficult to study owing to the challenging genetics of its native microbial hosts. Since the pathways for chloromethane catabolism show evidence of horizontal gene transfer, we reproduced this transfer process in the laboratory to generate new chloromethane-catabolizing strains in tractable hosts. We demonstrate that six putative accessory genes improve chloromethane catabolism, though heterologous expression of only one of the six is strictly necessary for growth on chloromethane. In contrast to growth of Methylobacterium strains with the closely-related compound dichloromethane, we find that chloride export does not limit growth on chloromethane and, in general, that the ability of a strain to grow on dichloromethane is uncorrelated with its ability to grow on chloromethane. This heterologous expression system allows us to investigate the components required for effective chloromethane catabolism and the factors that limit effective catabolism after horizontal transfer.

Adaptation under similar selective pressure often leads to comparable phenotypes. A longstanding question is whether such phenotypic repeatability entails similar (parallelism) or different genotypic changes (convergence). To better understand this, we characterized mutations that optimized expression of a plasmid-borne metabolic pathway during laboratory evolution of a bacterium. Expressing these pathway genes was essential for growth but came with substantial costs. Starting from overexpression, replicate populations founded by this bacterium all evolved to reduce expression. Despite this phenotypic repetitiveness, the underlying mutational spectrum was highly diverse. Analysis of these plasmid mutations identified three distinct means to modulate gene expression: (1) reducing the gene copy number, (2) lowering transcript stability, and (3) integration of the pathway-bearing plasmid into the host genome. Our study revealed diverse molecular changes beneath convergence to a simple phenotype. This complex genotype-phenotype mapping presents a challenge to inferring genetic evolution based solely on phenotypic changes.

Complete genome sequences are now available for many bacterial species that lack sophisticated genetic tools. We describe the development of a broad-host-range crelox system that allows antibiotic marker recycling in a variety of Gram-negative bacteria. This system consists of an allelic exchange vector bearing a kanamycin cassette flanked by loxP sites and a tetracycline-resistant IncP plasmid that provides expression of the Cre recombinase. We demonstrate this system by generating unmarked deletions of genes in two different bacteria, Methylobacterium extorquens AM1 and Burkholderia fungorum LB400. This new antibiotic marker recycling system offers the possibility of creating unmarked mutants in a wide variety of Gram-negative bacteria. Furthermore, marker recycling allows the generation of strains bearing multiple genetic manipulations in organisms for which few antibiotic markers are currently available.

Epistasis between mutations in two genes is thought to reflect an interdependence of their functions. While sometimes epistasis is predictable using mechanistic models, its roots seem, in general, hidden in the complex architecture of biological networks. Here, we ask how epistasis can be quantified based on the mathematical dependence of a system-level trait (e.g. fitness) on lower-level traits (e.g. molecular or cellular properties). We first focus on a model in which fitness is the difference between a benefit and a cost trait, both pleiotropically affected by mutations. We show that despite its simplicity, this model can be used to analytically predict certain properties of the ensuing distribution of epistasis, such as a global negative bias, resulting in antagonism between beneficial mutations, and synergism between deleterious ones. We next extend these ideas to derive a general expression for epistasis given an arbitrary functional dependence of fitness on other traits. This expression demonstrates how epistasis relative to fitness can emerge despite the absence of epistasis relative to lower level traits, leading to a formalization of the concept of independence between biological processes. Our results suggest that epistasis may be largely shaped by the pervasiveness of pleiotropic effects and modular organization in biological networks. [ABSTRACT FROM AUTHOR]

[This corrects the article DOI: 10.1371/journal.pgen.1008458.].

Normal cellular processes give rise to toxic metabolites that cells must mitigate. Formaldehyde is a universal stressor and potent metabolic toxin that is generated in organisms from bacteria to humans. Methylotrophic bacteria such as Methylorubrum extorquens face an acute challenge due to their production of formaldehyde as an obligate central intermediate of single-carbon metabolism. Mechanisms to sense and respond to formaldehyde were speculated to exist in methylotrophs for decades but had never been discovered. Here, we identify a member of the DUF336 domain family, named efgA for enhanced formaldehyde growth, that plays an important role in endogenous formaldehyde stress response in M. extorquens PA1 and is found almost exclusively in methylotrophic taxa. Our experimental analyses reveal that EfgA is a formaldehyde sensor that rapidly arrests growth in response to elevated levels of formaldehyde. Heterologous expression of EfgA in Escherichia coli increases formaldehyde resistance, indicating that its interaction partners are widespread and conserved. EfgA represents the first example of a formaldehyde stress response system that does not involve enzymatic detoxification. Thus, EfgA comprises a unique stress response mechanism in bacteria, whereby a single protein directly senses elevated levels of a toxic intracellular metabolite and safeguards cells from potential damage.

While microbiologists often make the simplifying assumption that genotype determines phenotype in a given environment, it is becoming increasingly apparent that phenotypic heterogeneity (in which one genotype generates multiple phenotypes simultaneously even in a uniform environment) is common in many microbial populations. The importance of phenotypic heterogeneity has been demonstrated in a number of model systems involving binary phenotypic states (e.g., growth/non-growth); however, less is known about systems involving phenotype distributions that are continuous across an environmental gradient, and how those distributions change when the environment changes. Here, we describe a novel instance of phenotypic diversity in tolerance to a metabolic toxin within wild-type populations of Methylobacterium extorquens, a ubiquitous phyllosphere methylotroph capable of growing on the methanol periodically released from plant leaves. The first intermediate in methanol metabolism is formaldehyde, a potent cellular toxin that is lethal in high concentrations. We have found that at moderate concentrations, formaldehyde tolerance in M. extorquens is heterogeneous, with a cell's minimum tolerance level ranging between 0 mM and 8 mM. Tolerant cells have a distinct gene expression profile from non-tolerant cells. This form of heterogeneity is continuous in terms of threshold (the formaldehyde concentration where growth ceases), yet binary in outcome (at a given formaldehyde concentration, cells either grow normally or die, with no intermediate phenotype), and it is not associated with any detectable genetic mutations. Moreover, tolerance distributions within the population are dynamic, changing over time in response to growth conditions. We characterized this phenomenon using bulk liquid culture experiments, colony growth tracking, flow cytometry, single-cell time-lapse microscopy, transcriptomics, and genome resequencing. Finally, we used mathematical modeling to better understand the processes by which cells change phenotype, and found evidence for both stochastic, bidirectional phenotypic diversification and responsive, directed phenotypic shifts, depending on the growth substrate and the presence of toxin.

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.

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.

Multi-species microbial communities play a critical role in human health, industry, and waste remediation. Recently, the evolution of synthetic consortia in the laboratory has enabled adaptation to be addressed in the context of interacting species. Using an engineered bacterial consortium, we repeatedly evolved cooperative genotypes and examined both the predictability of evolution and the phenotypes that determine community dynamics. Eight Salmonella enterica serovar Typhimurium strains evolved methionine excretion sufficient to support growth of an Escherichia coli methionine auxotroph, from whom they required excreted growth substrates. Non-synonymous mutations in metA, encoding homoserine trans-succinylase (HTS), were detected in each evolved S. enterica methionine cooperator and were shown to be necessary for cooperative consortia growth. Molecular modeling was used to predict that most of the non-synonymous mutations slightly increase the binding affinity for HTS homodimer formation. Despite this genetic parallelism and trend of increasing protein binding stability, these metA alleles gave rise to a wide range of phenotypic diversity in terms of individual versus group benefit. The cooperators with the highest methionine excretion permitted nearly two-fold faster consortia growth and supported the highest fraction of E. coli, yet also had the slowest individual growth rates compared to less cooperative strains. Thus, although the genetic basis of adaptation was quite similar across independent origins of cooperative phenotypes, quantitative measurements of metabolite production were required to predict either the individual-level growth consequences or how these propagate to community-level behavior.

How do bacteria regulate their cellular physiology in response to starvation? Here, we present a detailed characterization of Escherichia coli growth and starvation over a time-course lasting two weeks. We have measured multiple cellular components, including RNA and proteins at deep genomic coverage, as well as lipid modifications and flux through central metabolism. Our study focuses on the physiological response of E. coli in stationary phase as a result of being starved for glucose, not on the genetic adaptation of E. coli to utilize alternative nutrients. In our analysis, we have taken advantage of the temporal correlations within and among RNA and protein abundances to identify systematic trends in gene regulation. Specifically, we have developed a general computational strategy for classifying expression-profile time courses into distinct categories in an unbiased manner. We have also developed, from dynamic models of gene expression, a framework to characterize protein degradation patterns based on the observed temporal relationships between mRNA and protein abundances. By comparing and contrasting our transcriptomic and proteomic data, we have identified several broad physiological trends in the E. coli starvation response. Strikingly, mRNAs are widely down-regulated in response to glucose starvation, presumably as a strategy for reducing new protein synthesis. By contrast, protein abundances display more varied responses. The abundances of many proteins involved in energy-intensive processes mirror the corresponding mRNA profiles while proteins involved in nutrient metabolism remain abundant even though their corresponding mRNAs are down-regulated.

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.

When microbes acquire new abilities through horizontal gene transfer, the genes and pathways must function under conditions with which they did not coevolve. If newly-acquired genes burden the host, their utility will depend on further evolutionary refinement of the recombinant strain. We used laboratory evolution to recapitulate this process of transfer and refinement, demonstrating that effective use of an introduced dichloromethane degradation pathway required one of several mutations to the bacterial host that are predicted to increase chloride efflux. We then used this knowledge to identify parallel, beneficial mutations that independently evolved in two natural dichloromethane-degrading strains. Finally, we constructed a synthetic mobile genetic element carrying both the degradation pathway and a chloride exporter, which preempted the adaptive process and directly enabled effective dichloromethane degradation across diverse Methylobacterium environmental isolates. Our results demonstrate the importance of post–transfer refinement in horizontal gene transfer, with potential applications in bioremediation and synthetic biology.

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.

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.

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.

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 (13)C-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.

Recent studies have suggested that plant genomes have undergone potentially rampant horizontal gene transfer (HGT), especially in the mitochondrial genome. Parasitic plants have provided the strongest evidence of HGT, which appears to be facilitated by the intimate physical association between the parasites and their hosts. A recent phylogenomic study demonstrated that in the holoparasite Rafflesia cantleyi (Rafflesiaceae), whose close relatives possess the world's largest flowers, about 2.1% of nuclear gene transcripts were likely acquired from its obligate host. Here, we used next-generation sequencing to obtain the 38 protein-coding and ribosomal RNA genes common to the mitochondrial genomes of angiosperms from R. cantleyi and five additional species, including two of its closest relatives and two host species. Strikingly, our phylogenetic analyses conservatively indicate that 24%-41% of these gene sequences show evidence of HGT in Rafflesiaceae, depending on the species. Most of these transgenic sequences possess intact reading frames and are actively transcribed, indicating that they are potentially functional. Additionally, some of these transgenes maintain synteny with their donor and recipient lineages, suggesting that native genes have likely been displaced via homologous recombination. Our study is the first to comprehensively assess the magnitude of HGT in plants involving a genome (i.e., mitochondria) and a species interaction (i.e., parasitism) where it has been hypothesized to be potentially rampant. Our results establish for the first time that, although the magnitude of HGT involving nuclear genes is appreciable in these parasitic plants, HGT involving mitochondrial genes is substantially higher. This may represent a more general pattern for other parasitic plant clades and perhaps more broadly for angiosperms.

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.

Measureable rates of genome evolution are well documented in human pathogens but are less well understood in bacterial pathogens in the wild, particularly during and after host switches. Mycoplasma gallisepticum (MG) is a pathogenic bacterium that has evolved predominantly in poultry and recently jumped to wild house finches (Carpodacus mexicanus), a common North American songbird. For the first time we characterize the genome and measure rates of genome evolution in House Finch isolates of MG, as well as in poultry outgroups. Using whole-genome sequences of 12 House Finch isolates across a 13-year serial sample and an additional four newly sequenced poultry strains, we estimate a nucleotide diversity in House Finch isolates of only ∼2% of ancestral poultry strains and a nucleotide substitution rate of 0.8-1.2×10(-5) per site per year both in poultry and in House Finches, an exceptionally fast rate rivaling some of the highest estimates reported thus far for bacteria. We also found high diversity and complete turnover of CRISPR arrays in poultry MG strains prior to the switch to the House Finch host, but after the invasion of House Finches there is progressive loss of CRISPR repeat diversity, and recruitment of novel CRISPR repeats ceases. Recent (2007) House Finch MG strains retain only ∼50% of the CRISPR repertoire founding (1994-95) strains and have lost the CRISPR-associated genes required for CRISPR function. Our results suggest that genome evolution in bacterial pathogens of wild birds can be extremely rapid and in this case is accompanied by apparent functional loss of CRISPRs.

Genome reduction has been observed in many bacterial lineages that have adapted to specialized environments. The extreme genome degradation seen for obligate pathogens and symbionts appears to be dominated by genetic drift. In contrast, for free-living organisms with reduced genomes, the dominant force is proposed to be direct selection for smaller, streamlined genomes. Most variation in gene content for these free-living species is of "accessory" genes, which are commonly gained as large chromosomal islands that are adaptive for specialized traits such as pathogenicity. It is generally unclear, however, whether the process of accessory gene loss is largely driven by drift or selection. Here we demonstrate that selection for gene loss, and not a shortened genome, per se, drove massive, rapid reduction of accessory genes. In just 1,500 generations of experimental evolution, 80% of populations of Methylobacterium extorquens AM1 experienced nearly parallel deletions removing up to 10% of the genome from a megaplasmid present in this strain. The absence of these deletion events in a mutation accumulation experiment suggested that selection, rather than drift, has dominated the process. Reconstructing these deletions confirmed that they were beneficial in their selective regimes, but led to decreased performance in alternative environments. These results indicate that selection can be crucial in eliminating unnecessary genes during the early stages of adaptation to a specialized environment.

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 >10(5) 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.

Understanding the evolution of biological systems requires untangling the molecular mechanisms that connect genetic and environmental variations to their physiological consequences. Metal limitation across many environments, ranging from pathogens in the human body to phytoplankton in the oceans, imposes strong selection for improved metal acquisition systems. In this study, we uncovered the genetic and physiological basis of adaptation to metal limitation using experimental populations of Methylobacterium extorquens AM1 evolved in metal-deficient growth media. We identified a transposition mutation arising recurrently in 30 of 32 independent populations that utilized methanol as a carbon source, but not in any of the 8 that utilized only succinate. These parallel insertion events increased expression of a novel transporter system that enhanced cobalt uptake. Such ability ensured the production of vitamin B(12), a cobalt-containing cofactor, to sustain two vitamin B(12)-dependent enzymatic reactions essential to methanol, but not succinate, metabolism. Interestingly, this mutation provided higher selective advantages under genetic backgrounds or incubation temperatures that permit faster growth, indicating growth-rate-dependent epistatic and genotype-by-environment interactions. Our results link beneficial mutations emerging in a metal-limiting environment to their physiological basis in carbon metabolism, suggest that certain molecular features may promote the emergence of parallel mutations, and indicate that the selective advantages of some mutations depend generically upon changes in growth rate that can stem from either genetic or environmental influences.

BACKGROUND:Methylotrophy describes the ability of organisms to grow on reduced organic compounds without carbon-carbon bonds. The genomes of two pink-pigmented facultative methylotrophic bacteria of the Alpha-proteobacterial genus Methylobacterium, the reference species Methylobacterium extorquens strain AM1 and the dichloromethane-degrading strain DM4, were compared. METHODOLOGY/PRINCIPAL FINDINGS:The 6.88 Mb genome of strain AM1 comprises a 5.51 Mb chromosome, a 1.26 Mb megaplasmid and three plasmids, while the 6.12 Mb genome of strain DM4 features a 5.94 Mb chromosome and two plasmids. The chromosomes are highly syntenic and share a large majority of genes, while plasmids are mostly strain-specific, with the exception of a 130 kb region of the strain AM1 megaplasmid which is syntenic to a chromosomal region of strain DM4. Both genomes contain large sets of insertion elements, many of them strain-specific, suggesting an important potential for genomic plasticity. Most of the genomic determinants associated with methylotrophy are nearly identical, with two exceptions that illustrate the metabolic and genomic versatility of Methylobacterium. A 126 kb dichloromethane utilization (dcm) gene cluster is essential for the ability of strain DM4 to use DCM as the sole carbon and energy source for growth and is unique to strain DM4. The methylamine utilization (mau) gene cluster is only found in strain AM1, indicating that strain DM4 employs an alternative system for growth with methylamine. The dcm and mau clusters represent two of the chromosomal genomic islands (AM1: 28; DM4: 17) that were defined. The mau cluster is flanked by mobile elements, but the dcm cluster disrupts a gene annotated as chelatase and for which we propose the name "island integration determinant" (iid). CONCLUSION/SIGNIFICANCE:These two genome sequences provide a platform for intra- and interspecies genomic comparisons in the genus Methylobacterium, and for investigations of the adaptive mechanisms which allow bacterial lineages to acquire methylotrophic lifestyles.

Genome-scale analysis of predicted metabolic pathways has revealed the common occurrence of apparent redundancy for specific functional units, or metabolic modules. In many cases, mutation analysis does not resolve function, and instead, direct experimental analysis of metabolic flux under changing conditions is necessary. In order to use genome sequences to build models of cellular function, it is important to define function for such apparently redundant systems. Here we describe direct flux measurements to determine the role of redundancy in three modules involved in formaldehyde assimilation and dissimilation in a bacterium growing on methanol. A combination of deuterium and (14)C labeling was used to measure the flux through each of the branches of metabolism for growth on methanol during transitions into and out of methylotrophy. The cells were found to differentially partition formaldehyde among the three modules depending on the flux of methanol into the cell. A dynamic mathematical model demonstrated that the kinetic constants of the enzymes involved are sufficient to account for this phenomenon. We demonstrate the role of redundancy in formaldehyde metabolism and have uncovered a new paradigm for coping with toxic, high-flux metabolic intermediates: a dynamic, interconnected metabolic loop.