Ming-Chun Lee

Methylobacterium Genome Sequences: A Reference Blueprint to Investigate Microbial Metabolism of C1 Compounds from Natural and Industrial Sources

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.

Repeated, Selection-Driven Genome Reduction of Accessory Genes in Experimental Populations

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.

Asymmetric, Bimodal Trade-Offs During Adaptation of Methylobacterium to Distinct Growth Substrates

Trade-offs between selected and nonselected environments are often assumed to exist during adaptation. This phenomenon is prevalent in microbial metabolism, where many organisms have come to specialize on a narrow breadth of substrates. One wellstudied example is methylotrophic bacteria that can use single-carbon (C1) compounds as their sole source of carbon and energy, but generally use few, if any, multi-C compounds. Here, we use adaptation of experimental populations of the model methylotroph, Methylobacterium extorquens AM1, to C1 (methanol) or multi-C (succinate) compounds to investigate specialization and trade-offs between these two metabolic lifestyles. We found a general trend toward trade-offs during adaptation to succinate, but this was neither universal nor showed a quantitative relationship with the extent of adaptation. After 1500 generations, succinate-evolved strains had a remarkably bimodal distribution of fitness values on methanol: either an improvement comparable to the strains adapted on methanol or the complete loss of the ability to grow on C1 compounds. In contrast, adaptation to methanol resulted in no such trade-offs. Based on the substantial, asymmetric loss of C1 growth during growth on succinate, we suggest that the long-term maintenance of C1 metabolism across the genus Methylobacterium requires relatively frequent use of C1 compounds to prevent rapid loss.

Complete Genome Sequences of Six Strains of the Genus Methylobacterium

The complete and assembled genome sequences were determined for six strains of the alphaproteobacterial genus Methylobacterium, chosen for their key adaptations to different plant-associated niches and environmental constraints.

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

Methylobacterium extorquens strains are the best-studied methylotrophic model system, and their metabolism of single carbon compounds has been studied for over 50 years. Here we develop a new system for high-throughput batch culture of M. extorquens in microtiter plates by jointly optimizing the properties of the organism, the growth media and the culturing system. After removing cellulose synthase genes in M. extorquens strains AM1 and PA1 to prevent biofilm formation, we found that currently available lab automation equipment, integrated and managed by open source software, makes possible reliable estimates of the exponential growth rate. Using this system, we developed an optimized growth medium for M. extorquens using response surface methodologies. We found that media that used EDTA as a metal chelator inhibited growth and led to inconsistent culture conditions. In contrast, the new medium we developed with a PIPES buffer and metals chelated by citrate allowed for fast and more consistent growth rates. This new 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.

Synchronous Waves of Failed Soft Sweeps in the Laboratory: Remarkably Rampant Clonal Interference of Alleles at a Single Locus

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

FREQ-Seq: A Rapid, Cost-Effective, Sequencing-Based Method to Determine Allele Frequencies Directly from Mixed Populations

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

OASIS: an automated program for global investigation of bacterial and archaeal insertion sequences

Insertion sequences (ISs) are simple transposable elements present in most bacterial and archaeal genomes and play an important role in genomic evolution. The recent expansion of sequenced genomes offers the opportunity to study ISs comprehensively, but this requires efficient and accurate tools for IS annotation. We have developed an open-source program called OASIS, or Optimized Annotation System for Insertion Sequences, which automatically annotates ISs within sequenced genomes. OASIS annotations of 1737 bacterial and archaeal genomes offered an unprecedented opportunity to examine IS evolution. At a broad scale, we found that most IS families are quite widespread; however, they are not present randomly across taxa. This may indicate differential loss, barriers to exchange and/or insufficient time to equilibrate across clades. The number of ISs increases with genome length, but there is both tremendous variation and no increase in IS density for genomes >2 Mb. At the finer scale of recently diverged genomes, the proportion of shared IS content falls sharply, suggesting loss and/or emergence of barriers to successful cross-infection occurs rapidly. Surprisingly, even after controlling for 16S rRNA sequence divergence, the same ISs were more likely to be shared between genomes labeled as the same species rather than as different species.

Sign epistasis limits evolutionary trade-offs at the confluence of single- and multi-carbon metabolism in Methylobacterium extorquens AM1.

Adaptation of one set of traits is often accompanied by attenuation of traits important in other selective environments, leading to fitness trade‐offs. The mechanisms that either promote or prevent the emergence of trade‐offs remain largely unknown, and are difficult to discern in most systems. Here, we investigate the basis of trade‐offs that emerged during experimental evolution of Methylobacterium extorquens AM1 to distinct growth substrates. After 1500 generations of adaptation to a multi‐carbon substrate, succinate (S), many lineages had lost the ability to use one‐carbon compounds such as methanol (M), generating a mixture of M+ and M− evolved phenotypes. We show that trade‐offs in M− strains consistently arise via antagonistic pleiotropy through recurrent selection for loss‐of‐function mutations to ftfL (formate‐tetrahydrofolate ligase), which improved growth on S while simultaneously eliminating growth on M. But if loss of FtfL was beneficial, why were M trade‐offs not found in all populations? We discovered that eliminating FtfL was not universally beneficial on S, as it was neutral or even deleterious in certain evolved lineages that remained M+. This suggests that sign epistasis with earlier arising mutations prevented the emergence of mutations that drove trade‐offs through antagonistic pleiotropy, limiting the evolution of metabolic specialists in some populations.