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.https://doi.org/10.7554/eLife.04279.001
Many microbes can rapidly evolve to adapt to new or extreme habitats. Most often the characteristics that develop via evolution result from individuals inheriting new combinations of genes from their parents. However, many species can also acquire new genes through a process called horizontal gene transfer, where organisms that share an environment exchange potentially useful genes. The ability of bacteria to acquire genes by horizontal gene transfer is thought to be the main reason for the rapid evolution of antibiotic-resistant bacteria such as MRSA.
A horizontally transferred gene rarely works efficiently when first transferred, and may even cause problems for its new host. The host organism must therefore evolve further after receiving a new gene, but it is difficult to trace how this occurs.
A toxic chemical called dichloromethane (or DCM) has been used in many industries since World War II, and has caused widespread contamination to the environment. Strains from several microbial genera have been identified that have adapted to break down DCM and use it as their sole energy source. Moreover, a gene responsible for breaking down DCM appears to have been shared between different species by horizontal gene transfer.
In work that was presented earlier in 2014, researchers introduced this gene into several bacterial strains from the genus Methylobacterium that had not been previously exposed to DCM. However, the resulting new strains of bacteria still had difficulties growing on DCM. There were diverse ways that the bacteria could have been prevented from growing—for example, multiple by-products generated when DCM is broken down are highly toxic.
Now Michener et al.—including the researchers that performed the 2014 work—have experimentally evolved five of the modified strains of Methylobacterium to discover how they adapt over time to living on DCM. Examining the DNA sequences of these bacteria showed several mutations that both improved the ability of the bacteria to survive on DCM and helped the bacteria to more efficiently remove chloride ions from their systems. Chloride ions accumulate as DCM is broken down; therefore, Michener et al. suggest that being exposed to too much chloride prevented the bacteria in previous experiments from growing well. This idea is further supported by mutations found in previously isolated bacteria that live on DCM. These mutations increase the bacteria's ability to excrete chloride ions from their cells, and when Michener et al. removed these mutations, the bacteria no longer grew well on DCM.
Understanding how bacteria adapt once they have acquired the DCM-degrading gene could make it easier to help bacteria to mop up DCM and other contaminants in the environment.https://doi.org/10.7554/eLife.04279.002
Microbes frequently acquire genetic material from organisms that are distant evolutionary relatives but close ecological neighbors (Ochman et al., 2000; Gogarten et al., 2002; Smillie et al., 2011). These horizontally-transferred genes and pathways play important roles in processes ranging from the spread of antibiotic resistance to ecological differentiation (Forsberg et al., 2012; Shapiro et al., 2012). When, for example, a new environmental niche opens up through the introduction of a xenobiotic compound, horizontal gene transfer (HGT) can speed the assembly and dissemination of a corresponding catabolic pathway (Springael and Top, 2004). However, newly-acquired abilities may be costly for the host until they are carefully integrated into existing metabolic and regulatory networks (Kim and Copley, 2012; Yadid et al., 2013). The challenges associated with such adaptations are highlighted by many examples from metabolic engineering, where productively transferring genes and pathways into new hosts often requires significant post–transfer refinement (Chang et al., 2007; Michener et al., 2012).
Observing this process of transfer and refinement in nature is challenging, but we can recreate the same process in the laboratory using a combination of genetic engineering and experimental evolution. Laboratory evolution has been used extensively to study natural evolutionary processes. Recent experiments have highlighted the mechanisms by which evolution can optimize existing traits (Elena and Lenski, 2003), select for the emergence of a novel ability (Blount et al., 2012), and refine a rudimentary pathway (Quandt et al., 2014). Laboratory studies focusing on genetic exchanges between microbes have described the recombination of mutations within an evolving population (Winkler and Kao, 2012) or the replacement of an endogenous metabolic pathway with an alternative route (Chou et al., 2011). Meanwhile, metabolic engineers have used adaptive evolution as an engineering tool, including the evolution of highly modified strains, but under conditions that do not reproduce HGT (Fong et al., 2005; Trinh and Srienc, 2009; Lee and Palsson, 2010; Lee et al., 2012). In this work, we combined these approaches to study how microbes evolve following the acquisition through horizontal gene transfer of a novel metabolic ability by deliberately transferring a mobile metabolic pathway into a new host and then using experimental evolution to select for its efficient use.
We investigated this process of HGT followed by evolutionary refinement using dichloromethane (DCM) catabolism in Methylobacterium strains as an example. Dichloromethane became an important industrial solvent during World War II and is still the most prevalent chlorinated solvent today (www.eurochlor.org). Following widespread environmental contamination by DCM, several microbial strains have been isolated based on their ability to grow on DCM as the sole source of carbon and energy. Most notably, Methylobacterium extorquens DM4, a methylotrophic Alphaproteobacterium, has served as a reference model to elucidate the details of bacterial growth with DCM (Gälli and Leisinger, 1985; Muller et al., 2011a). M. extorquens DM4 expresses a dehalogenase, DcmA, whose gene shows clear signs of horizontal transfer (Schmid-Appert et al., 1997; Vuilleumier et al., 2009) and whose product, formaldehyde, feeds directly into central methylotrophic metabolism (Chistoserdova, 2011).
Past efforts to understand the genetics of growth on DCM used gene knockouts to uncover DCM-specific genes (Figure 1A) (Muller et al., 2011b). We complemented these efforts using a synthetic approach involving conjugal transfer of the dcmA gene into several naïve Methylobacterium strains (Figure 1B) (Kayser et al., 2000; Michener et al., 2014). These dcmA transconjugants showed DCM dehalogenase activity but grew poorly on DCM, and we previously investigated factors that might explain their limited growth (Michener et al., 2014). Growing on DCM using the dcmA catabolic pathway is very challenging for the host, as its products and intermediates include formaldehyde, hydrochloric acid, and the alkylating agent S-chloromethylglutathione (Figure 1F) (Kayser et al., 2000; Kayser and Vuilleumier, 2001). As DCM dehalogenase activity was generally high in transconjugant strains (Michener et al., 2014), we hypothesized that one or more of these stresses was limiting growth in the transconjugants and thus required post–transfer evolutionary refinement of the host to produce highly active DCM-degrading strains.
In this work, we used experimental evolution to identify the factors that were limiting growth in dcmA-containing transconjugants. Beginning from these ancestral transconjugants, we evolved replicate populations with DCM as the sole source of carbon and energy and obtained evolved isolates with substantially higher fitness on DCM (Figure 1C). We sequenced the genomes of these evolved strains to identify the mutations that increased fitness on DCM and then reconstructed these mutations in isogenic backgrounds to measure their effect on fitness (Figure 1D). Based on these mutations, we could infer both the primary limiting stress that the pathway placed on its host as well as the biochemical mechanisms that the cells used to overcome this stress. Finally, we used this knowledge to uncover naturally-occurring refining mutations in two DCM-degrading environmental isolates and to design an improved gene cassette that enhances DCM bioremediation by preempting the need for post–transfer refinement (Figure 1E).
We previously transferred the DCM dehalogenase dcmA into three strains of M. extorquens as well as three other species of Methylobacterium (Michener et al., 2014). Each of the transconjugants was significantly less fit than the natural isolate, M. extorquens DM4 (hereafter referred to as ‘DM4’). In this work, replicate populations of five of these transconjugants were serially propagated with DCM as the sole carbon and energy source to select for increased fitness on DCM. One of these five transconjugants, M. extorquens AM1, was initially unable to grow on DCM alone and instead was propagated on a mixture of DCM and methanol. After a total of 150 generations of growth, individual clones were isolated from each replicate population. We measured the fitness of each clone in direct competition with the reference strain, DM4, and selected the most-fit clone from each population for further analysis.
For four out of the five ancestral transconjugant strains, the fitness of isolates from each evolved population was significantly higher than the fitness of the ancestor during growth on DCM (Figure 2), increasing by 1.6- to 13-fold relative to the respective ancestral transconjugant. Methylobacterium radiotolerans was the only transconjugant that did not yield improved isolates after 150 generations, and these non-improved populations were not further characterized.
The genomes of all ancestral strains had previously been sequenced (Vuilleumier et al., 2009; Marx et al., 2012), so we used a combination of whole-genome resequencing and targeted Sanger sequencing to identify putative causal mutations in the evolved isolates. In the seven isolates that we resequenced, we repeatedly found mutations to the genes for the protein translocase secY, the chloride/proton antiporter clcA, and a hypothetical protein with a domain of unknown function (DUF599, which we renamed edgA for Evolved Dichloromethane Growth), as well as a single mutation to besA, a homolog of the eukaryal bestrophin family of chloride channels (Figure 3A). Sanger sequencing confirmed these mutations in the resequenced isolates and targeted sequencing of these loci demonstrated that each remaining isolate had mutations at one of these four loci. We note that additional mutations were identified in each resequenced isolate (Figure 3—source data 1).
Repeated mutations to the same loci in independently evolved replicates strongly suggested that the mutations were causal. However, to conclusively demonstrate causality we introduced several of these mutations into the ancestral strains M. extorquens AM1 and PA1 (hereafter referred to as ‘AM1’ and ‘PA1’). We focused here upon AM1 and PA1, as these two strains incurred mutations to all four loci described above and are genetically tractable. Since the secY mutation was not identified in PA1 nor the clcA mutation in AM1, we also constructed hybrid allele exchange vectors to move these mutations into the hosts in which they were not observed. We then measured the fitness of the reconstructed mutants against the natural isolate, DM4. Each of the mutations was highly beneficial, regardless of whether the mutation arose in that host during evolution (Figure 3B). In contrast, the mutations had only minor fitness effects during growth on an alternate carbon source, succinate, indicating that their beneficial effects were specific to the challenge of growing on DCM and were not broadly beneficial during laboratory growth (Figure 3—figure supplement 1).
While each of the single mutations was beneficial, the fitness effects were not additive. For example, isolate E2 contained mutations in both clcA and besA. Each of these mutations was beneficial alone, but the double mutant was no more fit than the besAE2 single mutant (Figure 3C). Similarly, when we constructed clcAE1/secYA1 and clcAE1/edgAE3 double mutants, we found that the double mutants were less fit than the clcAE1 single mutant (Figure 3—figure supplement 2).
The protein translocase gene secY was mutated in seven of thirteen evolved isolates, totaling six amino acid mutations and one in-frame deletion. In the crystal structure of the SecY homolog from Methanocaldococcus jannaschii (Van den Berg et al., 2004), these seven mutations cluster to the same region of the protein surface, at the interface of the channel pore and its blocking plug (Figure 4A). Several of these mutations have been identified and characterized in the Escherichia coli SecY homolog (with 58% amino acid identity to M. extorquens DM4 SecY) based on their ability to secrete proteins with mistranslated sequences (Emr et al., 1981). These mutations are predicted to disrupt interactions between the plug and pore, leading to a leaky channel (Smith et al., 2005), and detailed biochemical investigation of one such mutant in E. coli demonstrated that this mutation allows the facilitated diffusion of small anions, most notably chloride (Dalal and Duong, 2009).
Similarly, the promoter of the chloride/proton antiporter clcA was mutated in four of the thirteen evolved isolates. The two evolved PA1 isolates with clcA promoter mutations showed increased levels of clcA mRNA compared to either the ancestor or the PA1 isolate with a wild type clcA promoter (Figure 4B and Figure 4—figure supplement 1). However, the E. coli ClcA antiporter (with 33% amino acid identity to DM4) is capable of importing or exporting chloride, while necessarily transporting protons in the opposite direction (Accardi and Miller, 2004). To determine the directionality of ClcA-mediated transport in our strains, we used a pH biosensor to determine the effect of clcA overexpression on proton transport during growth with DCM. We observed a transient decrease in intracellular pH upon addition of DCM that was significantly larger in the strain overexpressing clcA (Figure 4C), consistent with ClcA importing protons and exporting chloride.
Next, we tested the effect of secY and clcA mutations on the chloride sensitivity of the mutant strains, since mutations that benefit the cell by facilitating chloride diffusion during growth on DCM may also result in detrimental levels of chloride import in high salt media. We measured the growth rate of strains on succinate with varying external chloride concentrations, using sucrose as a control for the effects of osmolarity. We were unable to detect a change in the sensitivity of AM1 mutants because the strain was already highly chloride sensitive (Michener et al., 2014). In PA1, the clcAE1 mutant was no more sensitive to NaCl or sucrose than the wild type strain (Figure 4D). In contrast, the PA1 secYA1 mutant was highly sensitive to NaCl, but only slightly more sensitive to sucrose (Figure 4D).
We next asked whether our laboratory evolution experiments had recapitulated the natural process of adaptation to growth on DCM, thereby allowing the identification of refining mutations in environmental isolates. The sequences of secY, besA, and edgA in the reference isolate, DM4, did not contain any of the mutations identified in the laboratory-evolved strains. However, we found that DM4 expressed significantly higher levels of clcA mRNA than the other three ancestral strains of M. extorquens tested (Figure 5A). The sequences upstream of clcA differ between all four strains, making it difficult to identify the causal mutations in DM4. Consequently, we exchanged a 140-bp region upstream of clcA (containing a total of five nucleotide changes) between DM4 and PA1. The resulting DM4 strain with the PA1 clcA promoter had decreased clcA expression and minimal growth with DCM as the sole carbon and energy source, while the PA1 strain with the DM4 clcA promoter had increased clcA expression and was significantly more fit than the wild type PA1 strain (Figure 5B). However, complementing the DM4 clcAE0 mutant strain with a clcA/dcmA expression plasmid restored growth on DCM (Figure 6—figure supplement 1).
Having identified necessary refining mutations in DM4, we examined the other known DCM-degrading Methylobacterium isolate, M. extorquens DM17 (Firsova et al., 2010), to determine whether a parallel event had occurred in its evolutionary history. When we sequenced the secY, besA, and edgA loci in this strain, the nucleotide sequences were exactly identical to M. extorquens AM1. Compared to M. extorquens AM1, however, the clcA locus of DM17 contained an insertion sequence 52 bp upstream of the clcA promoter, and we found that clcA expression was even higher in DM17 than in DM4 (Figure 5A).
Having demonstrated that clcA overexpression increases the fitness in several strains of M. extorquens during growth on DCM, we tested whether this effect would also hold in other environmental isolates across the Methylobacterium genus. We constructed a plasmid that expressed dcmA from the native DM4 promoter as well as clcA from a constitutive promoter. We then mated either our original dcmA plasmid (pJM10) or this new dcmA/clcA plasmid (pJM83) into a collection of natural Methylobacterium isolates. We found that every strain showed higher fitness on DCM with the dual-expression dcmA/clcA plasmid than with the single-expression dcmA plasmid (Figure 6).
Previous studies in experimental evolution generally focused on a single wild type microbial strain. In contrast, we evolved five different Methylobacterium strains, each of which carried the same horizontally-transferred metabolic pathway for DCM catabolism. By sequencing the evolved isolates, we identified mutations across diverse backgrounds that improve the ability to grow on DCM. Working backward from these mutations allowed us to determine that chloride accumulation appears to be the most critical factor limiting growth, and to identify multiple mechanisms that alleviate this stress. We found parallel responses at both the genetic and physiological level, indicating that the different strains faced similar challenges when forced to use their new metabolic capability and overcame these challenges in similar ways. Horizontal transfer of dcmA is not sufficient for efficient growth on DCM and, in both the laboratory and the environment, microbes must also acquire refining mutations to optimize the newly-acquired pathway.
Several lines of evidence suggest that chloride export was the key limiting factor during growth on DCM. First, the observed mutations in secY are predicted to turn the protein translocase into a leaky chloride channel, as demonstrated previously in E. coli (Dalal and Duong, 2009). Our growth rate measurements confirm that secY mutants are more sensitive to extracellular chloride, as we would expect if the mutant SecY facilitates chloride diffusion across the cell membrane. Second, mutations to the promoter of clcA increase clcA expression. The in vivo pH measurements upon DCM addition establish that ClcA imports protons, with a necessary concomitant export of chloride, and that increased ClcA expression leads to increased proton import and chloride export. Similarly, the E. coli response to extreme acid stress uses ClcA to export chloride, with other mechanisms used for proton export (Accardi and Miller, 2004). Third, BesA is predicted to be a chloride channel (Sun et al., 2002) and while we cannot yet explain the mechanism by which the identified mutation would affect its function, we hypothesize that the mutation may also allow increased chloride export. Finally, we observed that double mutants are no more fit than the best single mutant, suggesting that all the mutations we detected are addressing the same physiological limitation. Previous experiments found that mutations affecting the same pathway were less beneficial when combined (Chou et al., 2011), in some cases to the point of being collectively deleterious (Kvitek and Sherlock, 2011; Rokyta et al., 2011; Chou et al., 2014). Due to the antagonism of edgAE3 with clcAE1, it is likely that EdgA is also involved in chloride export.
While chloride export is not an unexpected stress during growth on DCM, there was no a priori reason to expect that chloride export was the major limiting stress across all of the strains that we tested. Multiple other stresses could have been dominant, including sensitivity to internal acidification, the solvent effects of DCM, changes in formaldehyde fluxes perturbing the native host metabolism, or the mutagenic effects of glutathione-conjugated DCM (Kayser and Vuilleumier, 2001). Indeed, our previous efforts to find a single physiological difference that explained the varied fitness of the ancestral dcmA transconjugants on DCM did not identify chloride sensitivity as a predictive phenotype for fitness during growth on DCM (Michener et al., 2014). There was no reason to expect that chloride was more toxic than protons, but using ClcA to export the chloride, even at the cost of importing yet more protons, was highly beneficial. It appears that excess intracellular chloride is either more toxic or less effectively managed than excess intracellular protons.
The apparent linkage between besA and chloride export may offer insights into a human disease. Mutations to human bestrophin homologs can lead to vitelliform macular dystrophy, but little is known about these proteins aside from their roles as chloride channels (Sun et al., 2002; Tsunenari et al., 2003). Bacterial homologs of the bestrophin family are similarly uncharacterized. Select bacterial homologs have been shown to affect swarming (Inoue et al., 2007) and sporulation (Rahn-Lee et al., 2009), but these systems have not explicitly replicated the role of the human homolog as a chloride channel. However, our Methylobacterium system may offer an opportunity to characterize the chloride transporter function of this protein family in a tractable host with a strong linkage between BesA function and fitness.
By design, our analysis has focused on mutations with very large fitness effects. Mutations that alleviate stresses imposed during growth on DCM had the potential to dramatically increase fitness, and restricting our evolution experiment to 150 generations predominantly selected these mutations. Our populations undoubtedly sampled mutations that would be generally beneficial during growth in our media and laboratory conditions, as has been seen in earlier experiments with M. extorquens AM1 (Lee et al., 2009), but the fitness effects of these mutations appear to have been minor compared to the DCM-specific mutations. For example, isolate N1 lost an entire 458 kbp plasmid, which was likely due to selection but with a small fitness benefit as compared to acquiring a secY mutation (Lee and Marx, 2012). The rarity of mutations generically beneficial during laboratory growth greatly simplified our analysis of the evolved isolates.
We expect that multiple lineages within a population independently arrived at solutions to increase chloride efflux. However, the size of the fitness effects produced evolutionary dynamics that resembled periodic selection (Atwood et al., 1951), in which the first large effect mutation to arise in the population rapidly fixed. Combinations of mutations were generally no more beneficial than the best single mutation, so we predominantly observed the result of a single selective sweep (Figure 3—figure supplement 3). Population E2 is the exception, where the clcAE2 mutation fixed first, followed by the besAE2 mutation (Figure 3—figure supplement 4). The clcAE2 mutant has a fitness of 0.6 as compared to 0.8 for besAE2 or the besAE2/clcAE2 double mutant (Figure 3C). Thus, while the double mutant is no more fit than the besAE2 single mutant, the besAE2 mutation is still highly beneficial in the clcAE2 background.
Experimental populations of M. radiotolerans did not increase in fitness over 150 generations of growth on DCM. When we tested this strain with the dual-expression clcA/dcmA plasmid, the fitness increased but the increase was smaller than for other strains (Figure 6—source data 1). We hypothesize that our inability to find improved isolates of M. radiotolerans stems from a combination of fewer accessible beneficial mutations, smaller fitness effects of the beneficial mutations, and a decreased sensitivity to the mutagenic effects of growing on DCM (Perez-Pantoja et al., 2013).
Experimental evolution is increasingly used to simulate natural evolutionary processes (Kawecki et al., 2012); however, there are few examples where mutations that arise during laboratory evolution can also be identified in natural populations (Wong et al., 2012; Traverse et al., 2013). We have shown that mutations to clcA can be highly beneficial during laboratory growth, and that parallel mutations to clcA have arisen in two independent DCM-degrading environmental isolates of M. extorquens. In the reference environmental isolate M. extorquens DM4, we can directly link nucleotide mutations in the clcA promoter to clcA overexpression and increased organismal fitness during growth on DCM. Demonstrating linkages between evolutionary processes observed in the laboratory and in nature will be increasingly important for the further development and application of experimental evolution.
Having identified clcA overexpression as a beneficial evolutionary refinement following acquisition of dcmA, we used this knowledge to develop a novel bioremediation strategy and designed new expression cassettes for bioremediation that are less dependent on subsequent adaptation for effective function. Rather than introducing an exogenous microbe into a contaminated site, we propose to introduce a genetic cassette to the indigenous microflora, an approach known as genetic bioaugmentation (Top et al., 2002; Ikuma and Gunsch, 2013). Our experiments clearly demonstrated that simply introducing the degradation pathway for a contaminant can be inefficient, as the recipients may be unprepared for the stresses produced by the pathway (Kayser et al., 2002; Michener et al., 2014). Instead, we propose to identify these stresses in the laboratory and then provide both the catabolic pathway and solutions to the most common limitations. In our example of DCM degradation, transconjugants that received both the catabolic gene dcmA and the chloride exporter gene clcA were more efficient at degrading DCM than those that received dcmA alone (Figure 6).
In summary, our results highlight the important role of evolutionary refinement in the horizontal transfer of a challenging catabolic pathway. Dichloromethane degradation led to chloride accumulation and required modifications to the recipient to export the chloride. Based on recent examples from metabolic engineering, we expect that these types of host–pathway interactions are common (Kizer et al., 2008; Ro et al., 2008; Verwaal et al., 2010; Michener et al., 2012; Bernhardt and Urlacher, 2014) and, consequently, that post–transfer refinement of microbial hosts is also widespread in nature. Our approach of deliberate horizontal transfer in the laboratory followed by experimental evolution offers an opportunity to identify both the limiting stresses and the mutations that can overcome these stresses. As screens and selections for complex microbial phenotypes are further developed, we expect that this approach will find many applications, both in the analysis of natural microbial isolates and in the optimization of genetically engineered microbes and mobile genetic elements.
E. coli cultures were grown in LB containing 10 g/l NaCl except for cultures containing the dcmA/clcA dual expression plasmid (pJM83), which were salt-sensitive and instead were grown in LB containing 0.5 g/l NaCl. Unless otherwise noted, Methylobacterium strains were grown in 10 ml cultures in a 50 ml flask at 30°C and 220 rpm in M-PIPES medium supplemented with 3.5 mM succinate or 5 mM DCM (Delaney et al., 2013a). Antibiotics were added to final concentrations of 10 µg/ml streptomycin, 12.5 µg/ml tetracycline, or 50 µg/ml kanamycin. All DCM cultures were grown in gas-tight 50 ml screw-top flasks sealed with Mininert valves (Supelco, Bellefonte, PA) and teflon tape. Before use, valves were sterilized with ethanol and dried in a laminar flow hood. Unless otherwise noted, chemicals were purchased from Sigma–Aldrich (St. Louis, MO) and enzymes from New England Biolabs (Ipswich, MA). The additional Methylobacterium strains presented here were either isolated from the University of Washington campus (L1, D21, D23, D24; CJM, unpublished) or from the surroundings of Woods Hole, MA (J-4-1, C-7-2, C-7-1, C-2-3, G-1-1, M-1-1; NF Delaney, unpublished). Strains and plasmids are listed in Supplementary file 1.
Transconjugant strains containing pJM10 were streaked to single colonies on M-PIPES + succinate + kanamycin plates. Liquid cultures were inoculated from single colonies into M-PIPES + succinate + kanamycin and grown to saturation. Cultures were then diluted 100× into M-PIPES + DCM to initiate the evolution experiment, starting 3–4 replicate flasks for each transconjugant. Every 3.5 days, cultures were diluted into fresh M-PIPES + DCM. Initial dilution factors were 16×, rising to 64× by the end of the experiment. Aliquots were taken at 12, 24, 36, 48, 60, 90, 120, and 150 generations and frozen at −80°C in 8% DMSO. After 150 generations, the cultures were plated on M-PIPES + succinate plates to isolate single colonies. Individual colonies were restreaked onto M-PIPES + succinate plates, then inoculated into liquid cultures in M-PIPES + succinate. Liquid cultures were diluted 100× into M-PIPES + DCM to test for DCM growth. For each replicate population, the isolate with the highest yield on DCM as assayed by OD600 was selected for further analysis.
The ancestral AM1 transconjugant was unable to grow on DCM alone. Consequently, the initial stages of evolution were performed in DCM supplemented with 0.75 mM methanol and kanamycin. For two of the populations, A2 and A3, the populations were able to grow on DCM alone after 60 generations of evolution, and the remaining 90 generations were conducted on DCM alone. For the final population, A1, the entire evolution experiment was performed with DCM + methanol + kanamycin.
Competition experiments were conducted to compare the fitness of the tested strain against M. extorquens DM4 ΔdcmA hptA::Venus + pJM10, as described previously (Michener et al., 2014). Briefly, strains were grown to saturation in M-PIPES + succinate + kanamycin, and then diluted 100× into M-PIPES + DCM. After 3 days, cultures were mixed at defined ratios in fresh M-PIPES + DCM and allowed to grow for 3 more days. A flow cytometer was used to measure the population ratios before and after the final round of growth on DCM. The fitness of the isolate was then calculated based on the change in population ratios and the fold-growth of the mixed population.
To prepare genomic DNA for sequencing, evolved isolates were streaked to single colonies on M-PIPES + succinate, and then grown to saturation in M-PIPES + succinate. Saturated cultures were centrifuged at 4500×g for 10 min, washed in 1 ml of water, and centrifuged at 8000×g for 3 min. After discarding the remaining supernatant, cell pellets were stored at −20°C overnight. Pellets were resuspended in 570 µl of TET (10 mM Tris pH 7.5, 1 mM EDTA, and 1% Triton X-100), heated to 90°C for 1 hr, then cooled to RT. Lysozyme was added to a final concentration of 2 mg/ml, and the suspension was incubated at 37°C for 30 min. Next, 30 µl of 10% SDS and 0.2 mg/ml of proteinase K were added, and the suspension was incubated for a further 1 hr at 37°C. After addition of 0.5 mg/ml RNAse A, the suspension was incubated at 37°C for 1 hr, heated to 90°C for 10 min, and 100 µl 5 M NaCl and 80 µl CTAB/NaCl (10% CTAB in 0.7 M NaCl) were added. After a further 10 min incubation at 90°C, the mixture was extracted twice with phenol/chloroform, precipitated with isopropanol, air dried, and resuspended in 10 mM Tris pH 7.5.
DNA sequencing was performed at the Microarray and Genomic Analysis Core Facility of the University of Utah and the IBEST Genomics Resources Core of the University of Idaho using a HiSeq 2000 and yielding 90×–200× coverage. Genomic DNA libraries were constructed using a TruSeq DNA Sample Prep LT kit (Illumina, San Diego, CA) following the manufacturer's instructions. Genome sequences were analyzed using breseq-0.24 (Deatherage and Barrick, 2014). Putative mutations were confirmed by PCR amplification from the chromosome and Sanger sequencing. To determine the first appearance of a given mutation, the chromosomal loci were amplified by PCR from the mixed population aliquots that had been frozen at intermediate time points during evolution. The fraction of the population with a given allele was determined by Sanger sequencing of these mixed PCR products, with an estimated detection limit of ∼5% of the population.
To construct plasmid pJM83, plasmid pJM40 was amplified in an around-the-horn PCR to remove the pHluorin/mCherry coding frame. The clcA gene was amplified from PA1 gDNA with primers that added homology to the pJM40 expression construct (promoter and terminator). The clcA insert was then cloned into the pJM40 PCR product using Gibson assembly (Gibson et al., 2009).
To construct the allele exchange vectors, plasmid pPS04 was digested with XbaI and SacI. The appropriate chromosomal locus, including ∼500 bp on either side of the desired mutation, was amplified by PCR using primers that added homology to the pPS04 vector backbone. The chromosomal amplicon was then cloned into the pPS04 backbone using Gibson assembly (Gibson et al., 2009). In the case of plasmids pJM66, pJM67, pJM75, pJM76, pJM88, and pJM89, the mutant allele was moved into a strain in which it did not occur. Rather than cloning a single ∼1000 bp chromosomal locus, these plasmids were constructed from three separate amplicons. These amplicons included a ∼500 bp upstream homology region amplified from the recipient chromosome, a small region containing the desired mutation amplified from the donor chromosome, and a second ∼500 bp downstream homology region amplified from the recipient chromosome. The PCR primers used added the appropriate homology to combine all three amplicons into the pPS04 backbone in a single Gibson assembly reaction.
Plasmids were mated into recipient Methylobacterium strains using triparental matings as described previously (Fulton et al., 1984). Allele exchanges were performed as described (Marx, 2008). Putative allele exchange mutants were confirmed by Sanger sequencing a PCR amplification of the modified chromosomal locus.
In vivo pH measurements were conducted as described previously (Michener et al., 2014). Briefly, an expression plasmid containing both dcmA and a pHluorin-mCherry translational fusion (pJM40) was mated into the appropriate strain. The transconjugant was grown to saturation in M-PIPES + succinate + tetracycline, then diluted 100× into M-PIPES + DCM, and grown for three days. Cultures were diluted to a final optical density of 0.01 before analysis on an LSRII flow cytometer (BD, Franklin Lakes, NJ). pHluorin was excited at 488 nm and measured at 530/30 nm. mCherry was excited at 561 nm and measured at 620/40 nm. Samples were gated for forward scatter, side scatter, and mCherry fluorescence. After determining the population fluorescence, DCM was added to a final concentration of 10 mM. The culture was quickly vortexed and returned to the flow cytometer, leading to an approximately 20 s delay between DCM addition and consistent fluorescence measurements. After DCM addition, the population fluorescence was monitored for a further 10 min.
For each strain, a standard curve was constructed by diluting the DCM-grown culture to a final OD600 of 0.01 in a solution composed of 30 mM buffer, 50 mM NaCl, 3 mM KCl, 10 µM valinomycin, and 10 µM nigericin. Buffers used were MES pH 5.1, MES pH 5.3, MES pH 5.5, MES pH 5.7, MES pH 6.1, PIPES pH 6.5, PIPES pH 6.9, and PIPES pH 7.3. The population mean fluorescence ratio (pHluorin fluorescence divided by mCherry fluorescence) was measured for each combination of strain and pH, then fit to a modified Henderson–Hasselbalch equation. The internal pH of the experimental samples was calculated by finding the fluorescence ratio of each cell, dividing the timecourse into 5 s intervals, calculating the population mean of the fluorescence ratio for each interval, and comparing that ratio to the standard curve to calculate the internal pH.
For qRT-PCR measurements, cultures were grown to saturation in M-PIPES plus succinate, and then diluted into the indicated media conditions. Upon reaching mid-log phase (typically at half the optical density of a saturated culture), cultures were centrifuged for 10 min at 4500×g and 4°C, washed once with 1 ml of water, and pelleted again at 8000×g for 3 min. Cell pellets were flash frozen in liquid nitrogen and stored at −80°C overnight.
RNA was extracted using an RNeasy Kit (Qiagen, Germantown, MD) and RNAse-free DNAse (Qiagen) according to the manufacturer's instructions. Total RNA was reverse transcribed using SuperScript III (Life Technologies, Carlsbad, CA), approximately 1 µg of RNA, and gene specific primers (rpsB FWD 5′-ACCAACTGGAAGACCATCTC-3′; rpsB REV 5′-CTTCTCGAGCTTGTCCTTCTCAC-3′; clcA FWD 5′-ATCGTCACCGAGATGACCCAG-3′; clcA REV 5′-CCAAGGTGTGATAGAGGCCG-3′) according to the manufacturer's directions. cDNA was quantified by qPCR, using a CFX-96 qPCR machine (Bio-Rad, Hercules, CA) and EvaGreen qPCR mix (Biotium, Hayward, CA). Three technical replicates were performed for each biological replicate. For each replicate, the clcA concentration was normalized by the rpsB concentration. For each strain, the average normalized clcA concentration was compared to the reference strain, typically PA1.
Growth rates were measured using an automated system as described previously (Delaney et al., 2013b). Briefly, strains were grown to saturation in M-PIPES + succinate and then diluted 64× into 640 µl of M-PIPES + succinate with the appropriate concentration of osmolyte (NaCl or sucrose) in a 48-well plate. After the cultures reached saturation, they were again diluted 64× into 640 µl of the appropriate media, with three replicate wells per condition. The optical density of these assay plates was monitored every 45 min until the cultures again reached saturation. The growth rate of the culture was calculated using CurveFitter (Delaney et al., 2013a).
Cytochromes P450 as promising catalysts for biotechnological application: chances and limitationsApplied Microbiology and Biotechnology 98:6185–6203.https://doi.org/10.1007/s00253-014-5767-7
Engineering Escherichia coli for production of functionalized terpenoids using plant P450sNature Chemical Biology 3:274–277.https://doi.org/10.1038/nchembio875
Identification of mutations in laboratory-evolved microbes from next-generation sequencing data using breseqMethods in Molecular Biology 1151:165–188.https://doi.org/10.1007/978-1-4939-0554-6_12
Clarity: an open-source manager for laboratory automationJournal of Laboratory Automation 18:171–177.https://doi.org/10.1177/2211068212460237
Evolution experiments with microorganisms: the dynamics and genetic bases of adaptationNature Reviews. Genetics 4:457–469.https://doi.org/10.1038/nrg1088
Analysis of the key functional genes in new aerobic degraders of dichloromethaneNature Reviews. Microbiology 79:66–72.https://doi.org/10.1134/S0026261710010091
In silico design and adaptive evolution of Escherichia coli for production of lactic acidBiotechnology and Bioengineering 91:643–648.https://doi.org/10.1002/bit.20542
Molecular cloning of a malyl coenzyme A lyase gene from Pseudomonas sp. strain AM1, a facultative methylotrophJournal of Bacteriology 160:718–723.
Specialized bacterial strains for the removal of dichloromethane from industrial wasteConservation & Recycling 8:91–100.
Prokaryotic evolution in light of gene transferMolecular Biology and Evolution 19:2226–2238.https://doi.org/10.1093/oxfordjournals.molbev.a004046
Successful genetic bioaugmentation with Pseudomonas putida for toluene degradation in soil columnsEnvironmental Chemistry Letters 11:365–370.https://doi.org/10.1007/s10311-013-0416-4
Genome-wide screening of genes required for swarming motility in Escherichia coli K-12Journal of Bacteriology 189:950–957.https://doi.org/10.1128/JB.01294-06
Dichloromethane metabolism and C1 utilization genes in Methylobacterium strainsNature Reviews. Microbiology 148:1915–1922.
Inhibitory cross-talk upon introduction of a new metabolic pathway into an existing metabolic networkProceedings of the National Academy of Sciences of USA 109:E2856–E2864.https://doi.org/10.1073/pnas.1208509109
Application of functional genomics to pathway optimization for increased isoprenoid productionApplied and Environmental Microbiology 74:3229–3241.https://doi.org/10.1128/AEM.02750-07
Adaptive evolution of Escherichia coli K-12 MG1655 during growth on a Nonnative carbon source, L-1,2-propanediolApplied and Environmental Microbiology 76:4158–4168.https://doi.org/10.1128/AEM.00373-10
Directed evolution of xylose isomerase for improved xylose catabolism and fermentation in the yeast Saccharomyces cerevisiaeApplied and Environmental Microbiology 78:5708–5716.https://doi.org/10.1128/AEM.01419-12
Complete genome sequences of six strains of the genus MethylobacteriumJournal of Bacteriology 194:4746–4748.https://doi.org/10.1128/JB.01009-12
Identification and treatment of heme depletion attributed to overexpression of a lineage of evolved P450 monooxygenasesProceedings of the National Academy of Sciences of USA 109:19504–19509.https://doi.org/10.1073/pnas.1212287109
Phylogeny poorly predicts the utility of a challenging horizontally transferred gene in Methylobacterium strainsJournal of Bacteriology 196:2101–2107.https://doi.org/10.1128/JB.00034-14
Dichloromethane-degrading bacteria in the genomic ageResearch in Microbiology 162:869–876.https://doi.org/10.1016/j.resmic.2011.01.008
Functional genomics of dichloromethane utilization in Methylobacterium extorquens DM4Environmental Microbiology 13:2518–2535.https://doi.org/10.1111/j.1462-2920.2011.02524.x
Recursive genomewide recombination and sequencing reveals a key refinement step in the evolution of a metabolic innovation in Escherichia coliProceedings of the National Academy of Sciences of USA 111:2217–2222.https://doi.org/10.1073/pnas.1314561111
The conserved sporulation protein YneE inhibits DNA replication in Bacillus subtilisJournal of Bacteriology 191:3736–3739.https://doi.org/10.1128/JB.00216-09
Association of newly discovered IS elements with the dichloromethane utilization genes of methylotrophic bacteriaNature Reviews. Microbiology 143:2557–2567.https://doi.org/10.1099/00221287-143-8-2557
Modeling the effects of prl mutations on the Escherichia coli SecY complexJournal of Bacteriology 187:6454–6465.https://doi.org/10.1128/JB.187.18.6454-6465.2005
The vitelliform macular dystrophy protein defines a new family of chloride channelsProceedings of the National Academy of Sciences of USA 99:4008–4013.https://doi.org/10.1073/pnas.052692999
Tangled bank of experimentally evolved Burkholderia biofilms reflects selection during chronic infectionsProceedings of the National Academy of Sciences of USA 110:E250–E259.https://doi.org/10.1073/pnas.1207025110
Metabolic engineering of Escherichia coli for efficient conversion of glycerol to ethanolApplied and Environmental Microbiology 75:6696–6705.https://doi.org/10.1128/AEM.00670-09
Structure-function analysis of the bestrophin family of anion channelsThe Journal of Biological Chemistry 278:41114–41125.https://doi.org/10.1074/jbc.M306150200
X-ray structure of a protein-conducting channelNature 427:36–44.
Harnessing recombination to speed adaptive evolution in Escherichia coliMetabolic Engineering 14:487–495.https://doi.org/10.1016/j.ymben.2012.07.004
Sequestration of a highly reactive intermediate in an evolving pathway for degradation of pentachlorophenolProceedings of the National Academy of Sciences of USA 110:E2182–E2190.https://doi.org/10.1073/pnas.1214052110
Michael LaubReviewing Editor; Massachusetts Institute of Technology, United States
eLife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see review process). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.
Thank you for sending your work entitled “Effective use of a horizontally-transferred pathway for dichloromethane catabolism requires post-transfer refinement” for consideration at eLife. Your article has been favorably evaluated by Detlef Weigel (Senior editor), Michael Laub (Reviewing editor), and 2 reviewers.
The Reviewing editor and the reviewers discussed their comments before we reached this decision, and the Reviewing editor has assembled the following comments to help you prepare a revised submission.
We all generally liked this manuscript and felt that it reported an interesting story with conclusions well-justified by the data shown. There were a few issues with regard to the writing that need attention in a revised manuscript. There were also a couple of issues raised in the reviews and during our subsequent discussion that need to be addressed experimentally. These primary issues are summarized here and, with one possible exception, are straightforward experiments to add:
1) We each had questions about the nature of the naturally occurring mutations in clcA in DM4 and DM17. You should provide stronger evidence that these mutations have the anticipated effects/consequences on clcA expression levels.
2) What are the fitness costs of dcmA and refinement alleles under non-DCM selective conditions?
3) Can you include FreqSeq data or at least some sampling of intermediates from the evolved lines to probe temporal dynamics of accumulated mutations?
4) One final item that came up is whether the clcA mutations identified actually make the transporter leaky, as predicted, leading to better chloride ion export. It was not clear to us how involved this type of experiment would be, but if it's reasonably straightforward it would be good to add this as well. You reference a chloride ion transport study in your paper and perhaps the approach taken previously could be implemented to solidify your conclusions about the clcA mutants.https://doi.org/10.7554/eLife.04279.018
1) We each had questions about the nature of the naturally occurring mutations in clcA in DM4 and DM17. You should provide stronger evidence that these mutations have the anticipated effects/consequences on clcA expression levels.
We have revised Figure 5 to include new qRT-PCR data showing the clcA expression in DM17 and in the promoter exchange strains. We can now show that both of the natural isolates, DM4 and DM17, have high clcA expression. Additionally, the promoter exchange mutants have the expected trend in clcA expression.
2) What are the fitness costs of dcmA and refinement alleles under non-DCM selective conditions?
We have added a new figure (Figure 3–figure supplement 2) and explanatory sentence describing the fitness effects of these mutations during growth on succinate. In short, none of the mutations have major fitness effects during growth on succinate compared to large effects on DCM.
3) Can you include FreqSeq data or at least some sampling of intermediates from the evolved lines to probe temporal dynamics of accumulated mutations?
We have added two new supplementary figures, Figure 3–figure supplement 3 and Figure 3–figure supplement 4, to show the allele frequency dynamics of selected populations. Of note, we can connect changes in allele frequency of the M. nodulans isolates to changes in bulk properties of the evolving population, such as the yield on DCM. Additionally, we include allele frequency measurements for the single isolate, E2, in which we find two highly beneficial mutations (as shown in Figure 3C). The remaining populations have less informative allele frequency dynamics, as we only observe a single mutation rapidly sweeping to fixation.
4) One final item that came up is whether the clcA mutations identified actually make the transporter leaky, as predicted, leading to better chloride ion export. It was not clear to us how involved this type of experiment would be, but if it's reasonably straightforward it would be good to add this as well. You reference a chloride ion transport study in your paper and perhaps the approach taken previously could be implemented to solidify your conclusions about the clcA mutants.
We assume that the reviewers were referencing the mutations to secY that we propose yield a leaky anion channel. While we agree that directly measuring the chloride permeability of inverted membrane vesicles containing the various SecY mutants would strengthen our hypothesis, running the necessary experiments would be quite challenging. We estimate that it would take an additional collaborator and a year’s work from a PhD student to run these experiments as there is no precedent for doing any of this work with inverted vesicles in Methylobacterium, or indeed for many organisms outside of E. coli. Instead, we direct the reviewers’ attention to our in vivo experiments where we demonstrate increased sensitivity to external NaCl as an alternative demonstration of facilitated chloride diffusion.https://doi.org/10.7554/eLife.04279.019
- Joshua K Michener
- Aline A Camargo Neves
- Aline A Camargo Neves
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
The authors acknowledge financial support from the National Institutes of Health (F32 GM106629 to JKM), the Fundação de Amparo à Pesquisa do Estado de São Paulo, and the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (fellowships to AACN). The authors thank Y Trotsenko for providing M. extorquens DM17, C Miller for advice on testing ClcA directionality, DD Nayak for outlining several of the early concepts of this project, NF Delaney for providing six Methylobacterium isolates, and members of the Marx Laboratory for offering helpful comments on the manuscript. Genome sequencing at the University of Idaho IBEST Genomics Resources Core was supported by grants from the National Center for Research Resources (P20RR016448) and the National Institute of General Medical Sciences (P20GM103397) from the National Institutes of Health.
- Michael Laub, Massachusetts Institute of Technology, United States
© 2014, Michener et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.