Multiphyletic origins of methylotrophy in Alphaproteobacteria,exemplified by comparative genomics of Lake Washington isolates

We sequenced the genomes of 19 methylotrophic isolates from Lake Washington, which belong to nine genera within eight families of the Alphaproteobacteria, two of the families being the newly proposed families. Comparative genomic analysis with a focus on methylotrophy metabolism classifies these strains into heterotrophic and obligately or facultatively autotrophic methylotrophs. The most persistent metabolic modules enabling methylotrophy within this group are the N‐methylglutamate pathway, the two types of methanol dehydrogenase (MxaFI and XoxF), the tetrahydromethanopterin pathway for formaldehyde oxidation, the serine cycle and the ethylmalonyl‐CoA pathway. At the same time, a great potential for metabolic flexibility within this group is uncovered, with different combinations of these modules present. Phylogenetic analysis of key methylotrophy functions reveals that the serine cycle must have evolved independently in at least four lineages of Alphaproteobacteria and that all methylotrophy modules seem to be prone to lateral transfers as well as deletions.     

Modularity of methylotrophy, revisited

Methylotrophy is a metabolic capability possessed by microorganisms that allows them to build biomass and to obtain energy from organic substrates containing no carbon-carbon bonds (C1 compounds, such as methane, methanol, etc.). This phenomenon in microbial physiology has been a subject of study for over 100 years, elucidating a set of well-defined enzymatic systems and pathways enabling this capability. The knowledge gained from the early genetic and genomic approaches to understanding methylotrophy pointed towards the existence of alternative enzymes/pathways for the specific metabolic goals. Different combinations of these systems in different organisms suggested that methylotrophy must be modular in its nature. More recent insights from genomic analyses, including the genomes representing novel types of methylotrophs, seem to reinforce this notion. This review integrates the new findings with the previously developed concept of modularity of methylotrophy.     

Investigating genomic and phenotypic parallelism between piscivorous and planktivorous lake trout (Salvelinus namaycush) ecotypes by means of RADseq and morphometrics analyses

Repeated adaptive ecological diversification has commonly been reported in fish and has often been associated with trophic niche diversity. The main goal of this study was to investigate the extent of parallelism in the genomic and phenotypic divergence between piscivorous and planktivorous lake trout ecotypes from Laurentian Shield lakes, Canada. This was achieved by documenting the extent of morphological differentiation using geometric morphometrics and linear measurements as well as the pattern of genomic divergence by means of RADseq genotyping (3925 filtered SNPs) in 12 lakes. Our results indicate that the two ecotypes evolved distinct body shape and several linear measurements in parallel. Neutral genetic differentiation was pronounced between all isolated populations (Mean F_ST = 0.433), indicating no or very limited migration and pronounced genetic drift. Significant genetic differentiation also suggested partial reproductive isolation between ecotypes in the two lakes where they are found in sympatry. Combining different outlier detection methods, we identified 48 SNPs putatively under divergent selection between ecotypes, among which 10 could be annotated and related to functions such as developmental processes and ionic regulation. Finally, our results indicate that parallel morphological divergence is accompanied by both parallel and nonparallel genomic divergence, which is associated with the use of different trophic niches between ecotypes. The results are also discussed in the context of management and conservation of this highly exploited species throughout northern North America.     

Widespread distribution of prophages signaling the potential for adaptability and pathogenicity evolution of Ralstonia solanacearum species complex

Integrated bacteriophages (prophages) can impact host cells, affecting their lifestyle, genomic diversity, and fitness. However, many basic aspects of how these organisms affect the host cell remain poorly understood. Ralstonia solanacearum is a gram-negative plant pathogenic bacterium that encompasses a great diversity of ecotypes regarded as a species complex (R. solanacearum Species Complex - RSSC). RSSC genomes have a mosaic structure containing numerous elements, signaling the potential for its evolution through horizontal gene transfer. Here, we analyzed 120 Ralstonia spp. genomes from the public database to identify prophage sequences. In total, 379 prophage-like elements were found in the chromosome and megaplasmid of Ralstonia spp. These elements encode genes related to host fitness, virulence factors, antibiotic resistance, and niche adaptation, which might contribute to RSSC adaptability. Prophage-like elements are widespread into the complex in different species and geographic origins, suggesting that the RSSC phages are ancestrally acquired. Complete prophages belonging to the families Inoviridae, Myoviridae, and Siphoviridae were found, being the members of Inoviridae the most abundant. Analysis of CRISPR-Cas spacer sequences demonstrated the presence of prophages sequences that indicate successive infection events during bacterial evolution. Besides complete prophages, we also demonstrated 14 novel putative prophages integrated into Ralstonia spp. genomes. Altogether, our results provide insights into the diversity of prophages in RSSC genomes and suggest that these elements may deeply affect the virulence and host adaptation and shaping the genomes among the strains of this important pathogen.     

Genomes of three methylotrophs from a single niche reveal the genetic and metabolic divergence of the methylophilaceae

The genomes of three representatives of the family Methylophilaceae, Methylotenera mobilis JLW8, Methylotenera versatilis 301, and Methylovorus glucosetrophus SIP3-4, all isolated from a single study site, Lake Washington in Seattle, WA, were completely sequenced. These were compared to each other and to the previously published genomes of Methylobacillus flagellatus KT and an unclassified Methylophilales strain, HTCC2181. Comparative analysis revealed that the core genome of Methylophilaceae may be as small as approximately 600 genes, while the pangenome may be as large as approximately 6,000 genes. Significant divergence between the genomes in terms of both gene content and gene and protein conservation was uncovered, including the varied presence of certain genes involved in methylotrophy. Overall, our data demonstrate that metabolic potentials can vary significantly between different species of Methylophilaceae, including organisms inhabiting the very same environment. These data suggest that genetic divergence among the members of this family may be responsible for their specialized and nonredundant functions in C? cycling, which in turn suggests means for their successful coexistence in their specific ecological niches.     

Investigating the extent of parallelism in morphological and genomic divergence among lake trout ecotypes in Lake Superior

Understanding the emergence of species through the process of ecological speciation is a central question in evolutionary biology which also has implications for conservation and management. Lake trout (Salvelinus namaycush) is renowned for the occurrence of different ecotypes linked to resource and habitat use throughout North America. We aimed to unravel the fine genetic structure of the four lake trout ecotypes in Lake Superior. A total of 486 individuals from four sites were genotyped at 6822 filtered SNPs using RADseq technology. Our results revealed different extent of morphological and genetic differentiation within the different sites. Overall, genetic differentiation was weak but significant and was on average three times higher between sites (mean F_ST = 0.016) than between ecotypes within sites (mean F_ST = 0.005) indicating higher level of gene flow or a more recent shared ancestor between ecotypes within each site than between populations of the same ecotype. Evidence of divergent selection was also found between ecotypes and/or in association with morphological variation. Outlier loci found in genes related to lipid metabolism and visual acuity were of particular interest in this context of ecotypic divergence. However, we did not find clear indication of parallelism at the genomic level, despite the presence of phenotypic parallelism among some ecotypes from different sampling sites. Overall, the occurrence of different levels of both genomic and phenotypic differentiation between ecotypes within each site with several differentiated loci linked to relevant biological functions supports the presence of a continuum of divergence in lake trout.     

Putative chromosomal rearrangements are associated primarily with ecotype divergence rather than geographic separation in an intertidal,poorly dispersing snail

Littorina saxatilis is becoming a model system for understanding the genomic basis of ecological speciation. The parallel formation of crab‐adapted ecotypes that exhibit partial reproductive isolation from wave‐adapted ecotypes has enabled genomic investigation of conspicuous shell traits. Recent genomic studies suggest that chromosomal rearrangements may enable ecotype divergence by reducing gene flow. However, the genomic architecture of traits that are divergent between ecotypes remains poorly understood. Here, we use 11,504 single nucleotide polymorphism (SNP) markers called using the recently released L. saxatilis genome to genotype 462 crab ecotype, wave ecotype and phenotypically intermediate Spanish L. saxatilis individuals with scored phenotypes. We used redundancy analysis to study the genetic architecture of loci associated with shell shape, shape corrected for size, shell size and shell ornamentation, and to compare levels of co‐association among different traits. We discovered 341 SNPs associated with shell traits. Loci associated with trait divergence between ecotypes were often located inside putative chromosomal rearrangements recently characterized in Swedish L. saxatilis. In contrast, we found that shell shape corrected for size varied primarily by geographic site rather than by ecotype and showed little association with these putative rearrangements. We conclude that genomic regions of elevated divergence inside putative rearrangements were associated with divergence of L. saxatilis ecotypes along steep environmental axes—consistent with models of adaptation with gene flow—but were not associated with divergence among the three geographical sites. Our findings support predictions from models indicating the importance of genomic regions of reduced recombination allowing co‐association of loci during ecological speciation with ongoing gene flow.     

Expressed Genome of Methylobacillus flagellatus as Defined through Comprehensive Proteomics and New Insights into Methylotrophy

In recent years, techniques have been developed and perfected for high-throughput identification of proteins and their accurate partial sequencing by shotgun nano-liquid chromatography-tandem mass spectrometry (nano-LC-MS/MS), making it feasible to assess global protein expression profiles in organisms with sequenced genomes. We implemented comprehensive proteomics to assess the expressed portion of the genome of Methylobacillus flagellatus during methylotrophic growth. We detected a total of 1,671 proteins (64% of the inferred proteome), including all the predicted essential proteins. Nonrandom patterns observed with the nondetectable proteins appeared to correspond to silent genomic islands, as inferred through functional profiling and genome localization. The protein contents in methylamine- and methanol-grown cells showed a significant overlap, confirming the commonality of methylotrophic metabolism downstream of the primary oxidation reactions. The new insights into methylotrophy include detection of proteins for the N-methylglutamate methylamine oxidation pathway that appears to be auxiliary and detection of two alternative enzymes for both the 6-phosphogluconate dehydrogenase reaction (GndA and GndB) and the formate dehydrogenase reaction (FDH1 and FDH4). Mutant analysis revealed that GndA and FDH4 are crucial for the organism''s fitness, while GndB and FDH1 are auxiliary.Methylotrophy is the ability of some microorganisms to grow on compounds containing no carbon-carbon bonds (C₁ compounds) as sole sources of carbon and energy. This type of metabolism requires special enzyme systems and pathways that accomplish three principal metabolic goals: (i) energy-producing primary oxidation of a C₁ substrate to formaldehyde or methyl transfer to produce methyl-tetrahydrofolate (H₄F), (ii) energy-producing oxidation of formaldehyde or methyl-H₄F to CO₂, and (iii) energy-consuming assimilation of formaldehyde, methylene-H₄F, and/or CO₂ into biomass (1, 22). In the course of microbial evolution, multiple pathways for each step have been assembled, and, theoretically, any combination of these should enable methylotrophy (9). Moreover, some organisms possess multiple pathways for certain metabolic tasks. For example, methane can be oxidized to methanol by either soluble or particulate methane monooxygenase, and many methanotrophs possess both enzymes (25). The presence of multiple formate dehydrogenases is typical among methylotrophs (10). Some methylotrophs, such as Methylococcus capsulatus, encode multiple pathways for C₁ assimilation (37). Methylotrophs employing the ribulose monophosphate (RuMP) cycle for formaldehyde assimilation can also oxidize formaldehyde via this pathway by employing a single additional reaction, catalyzed by 6-phosphogluconate dehydrogenase (Gnd) (1, 22). In addition, most RuMP cycle methylotrophs also possess a linear pathway for formaldehyde oxidation employing tetrahydromethanopterin (H₄MPT) as a cofactor (35).Methylobacillus flagellatus is a typical representative of RuMP cycle methylotrophs. It has a very limited substrate repertoire, growing robustly only on methanol or methylamine (16). Genomic analysis has revealed specific lesions in pathways for utilization of multicarbon compounds, confirming that, indeed, this organism must rely exclusively on methylotrophy to sustain its growth (11). At the same time, redundant methylotrophy pathways have been deduced from the genome. For example, in addition to the bona fide methanol dehydrogenase, four homologs of the large subunit are encoded (11). Besides the well-characterized methylamine dehydrogenase, an alternative system for methylamine oxidation is encoded, consisting of N-methylglutamate synthase and N-methylglutamate dehydrogenase (N-methylglutamate pathway) (21). Both cyclic and linear pathways for formaldehyde oxidation are encoded, along with two formate dehydrogenases. In addition, multiple terminal cytochrome oxidases are encoded (11).It has been assumed that metabolism of both methanol and methylamine, with the exception of the respective specific primary oxidation step, is carried out by M. flagellatus (and other methylotrophs) in exactly the same fashion (1, 22). However, this assumption has not been experimentally tested. The relative contributions of each of the redundant pathways encoded in the genome (for example, linear versus cyclic oxidation of formaldehyde) also remained unclear (7). These questions can be assessed through analysis of the expressed portion of the genome under specific growth conditions and especially precisely through a comprehensive analysis of the protein content in the cell (proteome) (3).In recent years, techniques have been developed and perfected for high-throughput identification of proteins and their accurate partial sequencing by shotgun nano-liquid chromatography-tandem mass spectrometry (nanoLC-MS/MS), making it feasible to assess relatively complete global protein expression profiles in prokaryotic organisms with sequenced genomes (2, 17). In this study, we employed comprehensive proteomics of M. flagellatus to determine what portion of the encoded proteome is detectable during methylotrophic growth under defined laboratory conditions and to obtain new insights into the physiology of the organism through this analysis.     

Genome of Methylobacillus flagellatus, molecular basis for obligate methylotrophy, and polyphyletic origin of methylotrophy

Along with methane, methanol and methylated amines represent important biogenic atmospheric constituents; thus, not only methanotrophs but also nonmethanotrophic methylotrophs play a significant role in global carbon cycling. The complete genome of a model obligate methanol and methylamine utilizer, Methylobacillus flagellatus (strain KT) was sequenced. The genome is represented by a single circular chromosome of approximately 3 Mbp, potentially encoding a total of 2,766 proteins. Based on genome analysis as well as the results from previous genetic and mutational analyses, methylotrophy is enabled by methanol and methylamine dehydrogenases and their specific electron transport chain components, the tetrahydromethanopterin-linked formaldehyde oxidation pathway and the assimilatory and dissimilatory ribulose monophosphate cycles, and by a formate dehydrogenase. Some of the methylotrophy genes are present in more than one (identical or nonidentical) copy. The obligate dependence on single-carbon compounds appears to be due to the incomplete tricarboxylic acid cycle, as no genes potentially encoding alpha-ketoglutarate, malate, or succinate dehydrogenases are identifiable. The genome of M. flagellatus was compared in terms of methylotrophy functions to the previously sequenced genomes of three methylotrophs, Methylobacterium extorquens (an alphaproteobacterium, 7 Mbp), Methylibium petroleiphilum (a betaproteobacterium, 4 Mbp), and Methylococcus capsulatus (a gammaproteobacterium, 3.3 Mbp). Strikingly, metabolically and/or phylogenetically, the methylotrophy functions in M. flagellatus were more similar to those in M. capsulatus and M. extorquens than to the ones in the more closely related M. petroleiphilum species, providing the first genomic evidence for the polyphyletic origin of methylotrophy in Betaproteobacteria.     

Genomic Architecture of Rapid Parallel Adaptation to Fresh Water in a Wild Fish

Rapid adaptation to novel environments may drive changes in genomic regions through natural selection. However, the genetic architecture underlying these adaptive changes is still poorly understood. Using population genomic approaches, we investigated the genomic architecture that underlies rapid parallel adaptation of Coilia nasus to fresh water by comparing four freshwater-resident populations with their ancestral anadromous population. Linkage disequilibrium network analysis and population genetic analyses revealed two putative large chromosome inversions on LG6 and LG22, which were enriched for outlier loci and exhibited parallel association with freshwater adaptation. Drastic frequency shifts and elevated genetic differentiation were observed for the two chromosome inversions among populations, suggesting that both inversions would undergo divergent selection between anadromous and resident ecotypes. Enrichment analysis of genes within chromosome inversions showed significant enrichment of genes involved in metabolic process, immunoregulation, growth, maturation, osmoregulation, and so forth, which probably underlay differences in morphology, physiology and behavior between the anadromous and freshwater-resident forms. The availability of beneficial standing genetic variation, large optimum shift between marine and freshwater habitats, and high efficiency of selection with large population size could lead to the observed rapid parallel adaptive genomic change. We propose that chromosomal inversions might have played an important role during the evolution of rapid parallel ecological divergence in the face of environmental heterogeneity in C. nasus. Our study provides insights into the genomic basis of rapid adaptation of complex traits in novel habitats and highlights the importance of structural genomic variants in analyses of ecological adaptation.     

A genomic island linked to ecotype divergence in Atlantic cod

The genomic architecture underlying ecological divergence and ecological speciation with gene flow is still largely unknown for most organisms. One central question is whether divergence is genome‐wide or localized in ‘genomic mosaics’ during early stages when gene flow is still pronounced. Empirical work has so far been limited, and the relative impacts of gene flow and natural selection on genomic patterns have not been fully explored. Here, we use ecotypes of Atlantic cod to investigate genomic patterns of diversity and population differentiation in a natural system characterized by high gene flow and large effective population sizes, properties which theoretically could restrict divergence in local genomic regions. We identify a genomic region of strong population differentiation, extending over approximately 20 cM, between pairs of migratory and stationary ecotypes examined at two different localities. Furthermore, the region is characterized by markedly reduced levels of genetic diversity in migratory ecotype samples. The results highlight the genomic region, or ‘genomic island’, as potentially associated with ecological divergence and suggest the involvement of a selective sweep. Finally, we also confirm earlier findings of localized genomic differentiation in three other linkage groups associated with divergence among eastern Atlantic populations. Thus, although the underlying mechanisms are still unknown, the results suggest that ‘genomic mosaics’ of differentiation may even be found under high levels of gene flow and that marine fishes may provide insightful model systems for studying and identifying initial targets of selection during ecological divergence.     

Genetic and Phenotypic Comparison of Facultative Methylotrophy between Methylobacterium extorquens Strains PA1 and AM1

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

Expanding our view of genomic diversity in Candidatus Accumulibacter clades

Enhanced biological phosphorus removal (EBPR) is an important industrial wastewater treatment process mediated by polyphosphate‐accumulating organisms (PAOs). Members of the genus Candidatus Accumulibacter are one of the most extensively studied PAO as they are commonly enriched in lab‐scale EBPR reactors. Members of different Accumulibacter clades are often enriched through changes in reactor process conditions; however, the two currently sequenced Accumulibacter genomes show extensive metabolic similarity. Here, we expand our understanding of Accumulibacter genomic diversity through recovery of eight population genomes using deep metagenomics, including seven from phylogenetic clades with no previously sequenced representative. Comparative genomic analysis revealed a core of shared genes involved primarily in carbon and phosphorus metabolism; however, each Accumulibacter genome also encoded a substantial number of unique genes (> 700 genes). A major difference between the Accumulibacter clades was the type of nitrate reductase encoded and the capacity to perform subsequent steps in denitrification. The Accumulibacter clade IIF genomes also contained acetaldehyde dehydrogenase that may allow ethanol to be used as carbon source. These differences in metabolism between Accumulibacter genomes provide a molecular basis for niche differentiation observed in lab‐scale reactors and may offer new opportunities for process optimization.     

Investigating niche and lineage diversification in widely distributed taxa: phylogeography and ecological niche modeling of the Peromyscus maniculatus species group

The relationship between lineage formation and variation in the ecological niche is a fundamental evolutionary question. Two prevailing hypotheses reflect this relationship: niche conservatism and niche divergence. Niche conservatism predicts a pattern where sister taxa will occupy similar niche spaces; whereas niche divergence predicts that sister taxa will occupy different niche spaces. Widely distributed species often show distinct phylogeographic structure, but little research has been conducted on how the environment may be related to these phylogenetic patterns. We investigated the relationship between lineage divergence and environmental space for the closely related species Peromyscus maniculatus and P. polionotus utilizing phylogenetic techniques and ecological niche modeling (ENM). We estimated the phylogenetic relationship among individuals based on complete cytochrome b sequences that represent individuals from a majority of the species ranges. Niche spaces that lineages occupy were estimated by using 12 environmental layers. Differences in niche space were tested using multivariate statistics based on location data, and ENMs were employed using maximum entropy algorithms. Two similarity indices estimated significant divergence in environmental space based on the ENM. Six geographically structured lineages were identified within P. maniculatus. Nested within P. maniculatus we found that P. polionotus recently diverged from a clade occupying central and western United States. We estimated that the majority of the genetic lineages occupy distinct environmental niches, which supports a pattern of niche divergence. Two sister taxa showed niche divergence and represent different ecomorphs, suggesting morphological, genetic and ecological divergence between the two lineages. Two other sister taxa were observed in the same environmental space based on multivariate statistics, suggesting niche conservatism. Overall our results indicate that a widely distributed species may exhibit both niche conservatism and niche divergence, and that most lineages seem to occupy distinct environmental niches.     

Methylotrophs in natural habitats: current insights through metagenomics

The focus of this review is on the recent data from the omics approaches, measuring the presence of methylotrophs in natural environments. Both Bacteria and Archaea are considered. The data are discussed in the context of the current knowledge on the biochemistry of methylotrophy and the physiology of cultivated methylotrophs. One major issue discussed is the recent metagenomic data pointing toward the activity of “aerobic” methanotrophs, such as Methylobacter, in microoxic or hypoxic conditions. A related issue of the metabolic distinction between aerobic and “anaerobic” methylotrophy is addressed in the light of the genomic and metagenomic data for respective organisms. The role of communities, as opposed to single-organism activities in environmental cycling of single-carbon compounds, such as methane, is also discussed. In addition, the emerging issue of the role of non-traditional methylotrophs in global metabolism of single-carbon compounds and the role of methylotrophy pathways in non-methylotrophs is briefly mentioned.     

甲醇的生物利用与转化*

甲醇作为一种来源广泛、价格低廉、还原度高的非粮原料有望成为下一代生物制造的关键原料。利用合成生物学技术构建能够高效利用甲醇的重组微生物以实现从甲醇到高值化学品的生物转化已成国内外研究热点,但由于甲醇代谢过程的特殊性及复杂性,目前人工设计的甲基营养菌还难以实现以甲醇为唯一碳源进行生长及产物合成。基于对天然甲基营养菌甲醇代谢过程的分析,从甲醇脱氢酶的筛选与改造、甲醛同化途径的重构与优化、甲醇到化学品的生物转化几个方面对合成型甲基营养菌的构建策略及面临的挑战进行总结与分析,以期为今后合成型甲基营养菌的人工设计和利用提供一定的借鉴。     

Ecotype Evolution in Glossina palpalis Subspecies,Major Vectors of Sleeping Sickness

BackgroundThe role of environmental factors in driving adaptive trajectories of living organisms is still being debated. This is even more important to understand when dealing with important neglected diseases and their vectors.Conclusions/SignificanceThese results provide an opportunity to examine whether new tsetse fly ecotypes might display different behaviour, dispersal patterns, host preferences and vectorial capacities. This work also urges a revision of taxonomic status of Glossina palpalis subspecies and highlights again how fast ecological divergence can be, especially in host-parasite-vector systems.     

Patterns and implications of gene gain and loss in the evolution of Prochlorococcus

Prochlorococcus is a marine cyanobacterium that numerically dominates the mid-latitude oceans and is the smallest known oxygenic phototroph. Numerous isolates from diverse areas of the world's oceans have been studied and shown to be physiologically and genetically distinct. All isolates described thus far can be assigned to either a tightly clustered high-light (HL)-adapted clade, or a more divergent low-light (LL)-adapted group. The 16S rRNA sequences of the entire Prochlorococcus group differ by at most 3%, and the four initially published genomes revealed patterns of genetic differentiation that help explain physiological differences among the isolates. Here we describe the genomes of eight newly sequenced isolates and combine them with the first four genomes for a comprehensive analysis of the core (shared by all isolates) and flexible genes of the Prochlorococcus group, and the patterns of loss and gain of the flexible genes over the course of evolution. There are 1,273 genes that represent the core shared by all 12 genomes. They are apparently sufficient, according to metabolic reconstruction, to encode a functional cell. We describe a phylogeny for all 12 isolates by subjecting their complete proteomes to three different phylogenetic analyses. For each non-core gene, we used a maximum parsimony method to estimate which ancestor likely first acquired or lost each gene. Many of the genetic differences among isolates, especially for genes involved in outer membrane synthesis and nutrient transport, are found within the same clade. Nevertheless, we identified some genes defining HL and LL ecotypes, and clades within these broad ecotypes, helping to demonstrate the basis of HL and LL adaptations in Prochlorococcus. Furthermore, our estimates of gene gain events allow us to identify highly variable genomic islands that are not apparent through simple pairwise comparisons. These results emphasize the functional roles, especially those connected to outer membrane synthesis and transport that dominate the flexible genome and set it apart from the core. Besides identifying islands and demonstrating their role throughout the history of Prochlorococcus, reconstruction of past gene gains and losses shows that much of the variability exists at the “leaves of the tree,” between the most closely related strains. Finally, the identification of core and flexible genes from this 12-genome comparison is largely consistent with the relative frequency of Prochlorococcus genes found in global ocean metagenomic databases, further closing the gap between our understanding of these organisms in the lab and the wild.     

Comparative genomics of 274 Vibrio cholerae genomes reveals mobile functions structuring three niche dimensions

Background

Vibrio cholerae is a globally dispersed pathogen that has evolved with humans for centuries, but also includes non-pathogenic environmental strains. Here, we identify the genomic variability underlying this remarkable persistence across the three major niche dimensions space, time, and habitat.

Results

Taking an innovative approach of genome-wide association applicable to microbial genomes (GWAS-M), we classify 274 complete V. cholerae genomes by niche, including 39 newly sequenced for this study with the Ion Torrent DNA-sequencing platform. Niche metadata were collected for each strain and analyzed together with comprehensive annotations of genetic and genomic attributes, including point mutations (single-nucleotide polymorphisms, SNPs), protein families, functions and prophages.

Conclusions

Our analysis revealed that genomic variations, in particular mobile functions including phages, prophages, transposable elements, and plasmids underlie the metadata structuring in each of the three niche dimensions. This underscores the role of phages and mobile elements as the most rapidly evolving elements in bacterial genomes, creating local endemicity (space), leading to temporal divergence (time), and allowing the invasion of new habitats. Together, we take a data-driven approach for comparative functional genomics that exploits high-volume genome sequencing and annotation, in conjunction with novel statistical and machine learning analyses to identify connections between genotype and phenotype on a genome-wide scale.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-654) contains supplementary material, which is available to authorized users.     

Mixing and matching methylotrophic enzymes to design a novel methanol utilization pathway in E. coli

One-carbon (C1) compounds, such as methanol, have recently gained attention as alternative low-cost and non-food feedstocks for microbial bioprocesses. Considerable research efforts are thus currently focused on the generation of synthetic methylotrophs by transferring methanol assimilation pathways into established bacterial production hosts. In this study, we used an iterative combination of dry and wet approaches to design, implement and optimize this metabolic trait in the most common chassis, E. coli. Through in silico modelling, we designed a new route that “mixed and matched” two methylotrophic enzymes: a bacterial methanol dehydrogenase (Mdh) and a dihydroxyacetone synthase (Das) from yeast. To identify the best combination of enzymes to introduce into E. coli, we built a library of 266 pathway variants containing different combinations of Mdh and Das homologues and screened it using high-throughput ¹³C-labeling experiments. The highest level of incorporation of methanol into central metabolism intermediates (e.g. 22% into the PEP), was obtained using a variant composed of a Mdh from A. gerneri and a codon-optimized version of P. angusta Das. Finally, the activity of this new synthetic pathway was further improved by engineering strategic metabolic targets identified using omics and modelling approaches. The final synthetic strain had 1.5 to 5.9 times higher methanol assimilation in intracellular metabolites and proteinogenic amino acids than the starting strain did. Broadening the repertoire of methanol assimilation pathways is one step further toward synthetic methylotrophy in E. coli.