Genome sequence of a mesophilic hydrogenotrophic methanogen Methanocella paludicola, the first cultivated representative of the order Methanocellales

We report complete genome sequence of a mesophilic hydrogenotrophic methanogen Methanocella paludicola, the first cultured representative of the order Methanocellales once recognized as an uncultured key archaeal group for methane emission in rice fields. The genome sequence of M. paludicola consists of a single circular chromosome of 2,957,635 bp containing 3004 protein-coding sequences (CDS). Genes for most of the functions known in the methanogenic archaea were identified, e.g. a full complement of hydrogenases and methanogenesis enzymes. The mixotrophic growth of M. paludicola was clarified by the genomic characterization and re-examined by the subsequent growth experiments. Comparative genome analysis with the previously reported genome sequence of RC-I(MRE50), which was metagenomically reconstructed, demonstrated that about 70% of M. paludicola CDSs were genetically related with RC-I(MRE50) CDSs. These CDSs included the genes involved in hydrogenotrophic methane production, incomplete TCA cycle, assimilatory sulfate reduction and so on. However, the genetic components for the carbon and nitrogen fixation and antioxidant system were different between the two Methanocellales genomes. The difference is likely associated with the physiological variability between M. paludicola and RC-I(MRE50), further suggesting the genomic and physiological diversity of the Methanocellales methanogens. Comparative genome analysis among the previously determined methanogen genomes points to the genome-wide relatedness of the Methanocellales methanogens to the orders Methanosarcinales and Methanomicrobiales methanogens in terms of the genetic repertoire. Meanwhile, the unique evolutionary history of the Methanocellales methanogens is also traced in an aspect by the comparative genome analysis among the methanogens.     

Higher-level classification of the Archaea: evolution of methanogenesis and methanogens

We used a phylogenetic approach to analyze the evolution of methanogenesis and methanogens. We show that 23 vertically transmitted ribosomal proteins do not support the monophyly of methanogens, and propose instead that there are two distantly related groups of extant archaea that produce methane, which we have named Class I and Class II. Based on this finding, we subsequently investigated the uniqueness of the origin of methanogenesis by studying both the enzymes of methanogenesis and the proteins that synthesize its specific coenzymes. We conclude that hydrogenotrophic methanogenesis appeared only once during evolution. Genes involved in the seven central steps of the methanogenic reduction of carbon dioxide (CO(2)) are ubiquitous in methanogens and share a common history. This suggests that, although extant methanogens produce methane from various substrates (CO(2), formate, acetate, methylated C-1 compounds), these archaea have a core of conserved enzymes that have undergone little evolutionary change. Furthermore, this core of methanogenesis enzymes seems to originate (as a whole) from the last ancestor of all methanogens and does not appear to have been horizontally transmitted to other organisms or between members of Class I and Class II. The observation of a unique and ancestral form of methanogenesis suggests that it was preserved in two independent lineages, with some instances of specialization or added metabolic flexibility. It was likely lost in the Halobacteriales, Thermoplasmatales and Archaeoglobales. Given that fossil evidence for methanogenesis dates back 2.8 billion years, a unique origin of this process makes the methanogenic archaea a very ancient taxon.     

从宏基因组推测扎龙湿地未培养甲烷古菌Rice Cluster II的代谢途径与可能的盐碱适应性

【背景】湿地是重要的甲烷排放源,因为其中栖息着各种产甲烷古菌。已知未培养甲烷古菌Rice Cluster Ⅱ (RCⅡ)类群广泛分布于低温酸性和北方泥炭藓湿地、淡水湿地及草本沼泽等环境,但它们在低温盐碱湿地中的分布及代谢途径尚未知。【目的】分析扎龙盐碱湿地未培养甲烷古菌RCⅡ类群的多样性、推测产甲烷代谢途径及其潜在的盐碱适应机制。【方法】16S rRNA基因扩增子测序分析扎龙湿地土壤中甲烷古菌群组成;构建16S rRNA基因克隆文库分析扎龙湿地土壤RCⅡ的多样性;宏基因组分析推测RCⅡ古菌编码的产甲烷途径及与耐盐碱相关物质的合成基因。【结果】16SrRNA基因高通量测序发现未培养甲烷古菌的RCⅡ类群占扎龙盐碱湿地总甲烷古菌的13.280%±0.019%;系统发育学分析表明该湿地的RCⅡ由3个分支组成;宏基因组分析组装了2个优势的未培养RCⅡ的基因组,均含完整的氢还原二氧化碳产甲烷途径的基因,并编码海藻糖的转运与合成基因。【结论】扎龙盐碱湿地土壤富含未培养RCⅡ甲烷古菌,推测它们通过氢还原二氧化碳产甲烷,利用细胞内高的海藻糖适应盐碱环境。     

Methanogenesis by Methanosarcina acetivorans involves two structurally and functionally distinct classes of heterodisulfide reductase

Biochemical studies have revealed two distinct classes of Coenzyme B‐Coenzyme M heterodisulfide (CoB‐S‐S‐CoM) reductase (Hdr), a key enzyme required for anaerobic respiration in methane‐producing archaea. A cytoplasmic HdrABC enzyme complex is found in most methanogens, whereas a membrane‐bound HdrED complex is found exclusively in members of the order Methanosarcinales. Unexpectedly, genomic data indicate that multiple copies of both Hdr classes are found in all sequenced Methanosarcinales genomes. The Methanosarcina acetivorans hdrED1 operon is constitutively expressed and required for viability under all growth conditions examined, consistent with HdrED being the primary Hdr. HdrABC appears to be specifically involved in methylotrophic methanogenesis, based on reduced growth and methanogenesis rates of an hdrA1C1B1 mutant on methylotrophic substrates and downregulation of the genes during growth on acetate. This conclusion is further supported by phylogenetic analysis showing that the presence of hdrA1 in an organism is specifically correlated with the presence of genes for methylotrophic methanogenesis. Examination of mRNA abundance in methanol‐grown ΔhdrA1C1B1 strains relative to wild‐type revealed upregulation of genes required for synthesis of (di)methylsulfide and for transport and biosynthesis of CoB‐SH and CoM‐SH, suggesting that the mutant has a defect in electron transfer from ferredoxin to CoB‐S‐S‐CoM that causes cofactor limitation.     

产甲烷菌研究进展

产甲烷菌是重要的环境微生物,在自然界的碳素循环中起重要作用。迄今已有5种产甲烷菌基因组测序完成。基因组信息使人们对产甲烷茵的细胞结构、进化、代谢及环境适应性有了更深的理解。目前已知的甲烷生物合成途径有3种,它们以乙酸、甲基化合物、氢/二氧化碳为起始,通过不同的反应途径都形成了甲基辅酶M,在甲基辅酶M还原酶的催化下最终形成甲烷。     

The genome characteristics and predicted function of methyl-group oxidation pathway in the obligate aceticlastic methanogens, Methanosaeta spp

In this work, we report the complete genome sequence of an obligate aceticlastic methanogen, Methanosaeta harundinacea 6Ac. Genome comparison indicated that the three cultured Methanosaeta spp., M. thermophila, M. concilii and M. harundinacea 6Ac, each carry an entire suite of genes encoding the proteins involved in the methyl-group oxidation pathway, a pathway whose function is not well documented in the obligately aceticlastic methanogens. Phylogenetic analysis showed that the methyl-group oxidation-involving proteins, Fwd, Mtd, Mch, and Mer from Methanosaeta strains cluster with the methylotrophic methanogens, and were not closely related to those from the hydrogenotrophic methanogens. Quantitative PCR detected the expression of all genes for this pathway, albeit ten times lower than the genes for aceticlastic methanogenesis in strain 6Ac. Western blots also revealed the expression of fwd and mch, genes involved in methyl-group oxidation. Moreover, (13)C-labeling experiments suggested that the Methanosaeta strains might use the pathway as a methyl oxidation shunt during the aceticlastic metabolism. Because the mch mutants of Methanosarcina barkeri or M. acetivorans failed to grow on acetate, we suggest that Methanosaeta may use methyl-group oxidation pathway to generate reducing equivalents, possibly for biomass synthesis. An fpo operon, which encodes an electron transport complex for the reduction of CoM-CoB heterodisulfide, was found in the three genomes of the Methanosaeta strains. However, an incomplete protein complex lacking the FpoF subunit was predicted, as the gene for this protein was absent. Thus, F(420)H(2) was predicted not to serve as the electron donor. In addition, two gene clusters encoding the two types of heterodisulfide reductase (Hdr), hdrABC, and hdrED, respectively, were found in the three Methanosaeta genomes. Quantitative PCR determined that the expression of hdrED was about ten times higher than hdrABC, suggesting that hdrED plays a major role in aceticlastic methanogenesis.     

Natural selection on functional modules, a genome-wide analysis

Classically, the functional consequences of natural selection over genomes have been analyzed as the compound effects of individual genes. The current paradigm for large-scale analysis of adaptation is based on the observed significant deviations of rates of individual genes from neutral evolutionary expectation. This approach, which assumed independence among genes, has not been able to identify biological functions significantly enriched in positively selected genes in individual species. Alternatively, pooling related species has enhanced the search for signatures of selection. However, grouping signatures does not allow testing for adaptive differences between species. Here we introduce the Gene-Set Selection Analysis (GSSA), a new genome-wide approach to test for evidences of natural selection on functional modules. GSSA is able to detect lineage specific evolutionary rate changes in a notable number of functional modules. For example, in nine mammal and Drosophilae genomes GSSA identifies hundreds of functional modules with significant associations to high and low rates of evolution. Many of the detected functional modules with high evolutionary rates have been previously identified as biological functions under positive selection. Notably, GSSA identifies conserved functional modules with many positively selected genes, which questions whether they are exclusively selected for fitting genomes to environmental changes. Our results agree with previous studies suggesting that adaptation requires positive selection, but not every mutation under positive selection contributes to the adaptive dynamical process of the evolution of species.     

DETECTING EVOLUTIONARY TRANSFER OF GENES USING PhIGs1

Organisms have acquired plastids by convoluted paths that have provided multiple opportunities for gene transfer into a host nucleus from intracellular organisms, including the cyanobacterial ancestor of plastids, the proteobacterial ancestor of mitochondria, and both green and red algae whose engulfment has led to secondary acquisition of plastids. These gene movements are most accurately demonstrated by building phylogenetic trees that identify the evolutionary origin of each gene, and one effective tool for this is “PhIGs” (Phylogenetically Inferred Groups; http://PhIGs.org ), a set of databases and computer tools with a Web interface for whole‐genome evolutionary analysis. PhIGs takes as input gene sets of completely sequenced genomes, builds clusters of genes using a novel, graph‐based approach, and reconstructs the evolutionary relationships among all gene families. The user can view and download the sequence alignments, compare intron‐exon structures, and follow links to functional genomic databases. Currently, PhIGs contains 652,756 genes from 45 genomes grouped into 61,059 gene families. Graphical displays show the relative positions of these genes among genomes. PhIGs has been used to detect the evolutionary transfer of hundreds of genes from cyanobacteria and red algae into oömycete nuclear genomes, revealing that even though they have no plastids, their ancestors did, having secondarily acquired them from an intracellular red alga. A great number of genomes are soon to become available that are relevant to our broader understanding of the movement of genes among intracellular compartments after engulfing other organisms, and PhIGs will be an effective tool to interpret these gene movements.     

Shifts in methanogen community structure and function associated with long-term manipulation of sulfate and salinity in a hypersaline microbial mat

Methanogenesis was characterized in hypersaline microbial mats from Guerrero Negro, Baja California Sur, Mexico both in situ and after long-term manipulation in a greenhouse environment. Substrate addition experiments indicate methanogenesis to occur primarily through the catabolic demethylation of non-competitive substrates, under field conditions. However, evidence for the coexistence of other metabolic guilds of methanogens was obtained during a previous manipulation of sulfate concentrations. To fully characterize methanogenesis in these mats, in the absence of competition for reducing equivalents with sulfate-reducing microorganisms, we maintained microbial mats for longer than 1 year under conditions of lowered sulfate and salinity levels. The goal of this study was to assess whether observed differences in methane production during sulfate and salinity manipulation were accompanied by shifts in the composition of methanogen communities. Culture-independent techniques targeting methyl coenzyme M reductase genes ( mcrA ) were used to assess the dynamics of methanogen assemblages. Clone libraries from mats sampled in situ or maintained at field-like conditions in the greenhouse were exclusively composed of sequences related to methylotrophic members of the Methanosarcinales . Increases in pore water methane concentrations under conditions of low sulfate correlated with an observed increase in the abundance of putatively hydrogenotrophic mcrA , related to Methanomicrobiales . Geochemical and molecular data provide evidence of a significant shift in the metabolic pathway of methanogenesis from a methylotroph-dominated system in high-sulfate environments to a mixed community of methylotrophic and hydrogenotrophic methanogens under low sulfate conditions.     

Bacteriophage genomics

The past three years have seen an escalation in the number of sequenced bacteriophage genomes with more than 500 now in the NCBI phage database, representing a more than threefold increase since 2005. These span at least 70 different bacterial hosts, with two-thirds of the sequenced genomes of phages representing only eight bacterial hosts. Three key features emerge from the comparative analysis of these genomes. First, they span a very high degree of genetic diversity, suggesting early evolutionary origins. Second, the genome architectures are mosaic, reflecting an unusually high degree of horizontal genetic exchange in their evolution. Third, phage genomes contain a very high proportion of novel genetic sequences of unknown function, and probably represent the largest reservoir of unexplored genes. With an estimated 10(31) bacterial and archael viruses in the biosphere, our view of the virosphere will draw into sharper focus as further bacteriophage genomes are characterized.     

In silico description of cobalt and nickel assimilation systems in the genomes of methanogens

Methanogens are a diverse group of organisms that can live in a wide range of environments. Herein, cobalt and tungsten assimilation pathways have proposed to be established in the genomes of Methanococcus maripaludies C5 and Methanosarcina mazei Go1, respectively. All of the proteins involved in the proposed pathways were identified from public domain databases and then complied manually to reconstruct the pathways. The function of proteins with unknown function was assigned by a combined prediction approach. Totally, 17 proteins were identified to cobalt transport and assimilation processes whereas 7 proteins reported to tungsten assimilation system. Phylogenetic analysis of this study revealed that heavy metal transporter of methanogens could be evolved from closely related members in the different genera of methanogens. Nevertheless, genes encoding for metal resistance proteins could be originated from thermophilic and sulfur reducing bacteria. Many metalloenzymes in methanogens were very unique to the species of methanogens. It implied that these metal ions were utilized to produce the precursors for energy driven processes of methanogens. This study suggested that in combination of systems models and evolutionary inference can only correlate metabolic fluxes and physiological changes in methanogens. In silico models of this study will provide insights to design experiments for heavy metal assimilation processes of methanogens growing under heavy metal-rich environments and or in a laboratory condition.     

Reticulate representation of evolutionary and functional relationships between phage genomes

Bacteriophage genomes show pervasive mosaicism, indicating the importance of horizontal gene exchange in their evolution. Phage genomes represent unique combinations of modules, each of them with a different phylogenetic history. The traditional classification, based on a variety of criteria such as nucleic acid type (single/double-stranded DNA/RNA), morphology, and host range, appeared inconsistent with sequence analyses. With the genomic era, an ever increasing number of sequenced phages cannot be classified, in part due to a lack of morphological information and in part to the intrinsic incapability of tree-based methods to efficiently deal with mosaicism. This problem led some virologists to call for a moratorium on the creation of additional taxa in the order Caudovirales, in order to let virologists discuss classification schemes that might better suit phage evolution. In this context, we propose a framework for a reticulate classification of phages based on gene content. Starting from gene families, we built a weighted graph, where nodes represent phages and edges represent phage-phage similarities in terms of shared genes. We then apply various measures of graph topology to analyze the resulting graph. Most double-stranded DNA phages are found in a single component. The values of the clustering coefficient and closeness distinguish temperate from virulent phages, whereas chimeric phages are characterized by a high betweenness coefficient. We apply a 2-step clustering method to this graph to generate a reticulate classification of phages: Each phage is associated with a membership vector, which quantitatively characterizes its membership to the set of clusters. Furthermore, we cluster genes based on their "phylogenetic profiles" to define "evolutionary cohesive modules." In virulent phages, evolutionary modules span several functional categories, whereas in temperate phages they correspond better to functional modules. Moreover, despite the fact that modules only cover a fraction of all phage genes, phage groups can be distinguished by their different combination of modules, serving the bases for a higher level reticulate classification. These 2 classification schemes provide an automatic and dynamic way of representing the relationships within the phage population and can be extended to include newly sequenced phage genomes, as well as other types of genetic elements.     

The Origin and Evolution of Metabolic Pathways: Why and How did Primordial Cells Construct Metabolic Routes?

The emergence and evolution of metabolic pathways represented a crucial step in molecular and cellular evolution. In fact, the exhaustion of the prebiotic supply of amino acids and other compounds that were likely present on the primordial Earth imposed an important selective pressure, favoring those primordial heterotrophic cells that became able to synthesize those molecules. Thus, the emergence of metabolic pathways allowed primitive organisms to become increasingly less dependent on exogenous sources of organic compounds. Comparative analyses of genes and genomes from organisms belonging to Archaea, Bacteria, and Eukarya reveal that, during evolution, different forces and molecular mechanisms might have driven the shaping of genomes and the emergence of new metabolic abilities. Among these gene elongations, gene and operon duplications played a crucial role since they can lead to the (immediate) appearance of new genetic material that, in turn, might undergo evolutionary divergence, giving rise to new genes coding for new metabolic abilities. Concerning the mechanisms of pathway assembly, both the analysis of completely sequenced genomes and directed evolution experiments strongly support the patchwork hypothesis, according to which metabolic pathways have been assembled through the recruitment of primitive enzymes that could react with a wide range of chemically related substrates. However, the analysis of the structure and organization of genes belonging to ancient metabolic pathways, such as histidine biosynthesis, suggests that other different hypothesis, i.e., the retrograde hypothesis, may account for the evolution of some steps within metabolic pathways.     

Comparative Genomic Analysis of 60 Mycobacteriophage Genomes: Genome Clustering, Gene Acquisition, and Gene Size

Mycobacteriophages are viruses that infect mycobacterial hosts. Expansion of a collection of sequenced phage genomes to a total of 60—all infecting a common bacterial host—provides further insight into their diversity and evolution. Of the 60 phage genomes, 55 can be grouped into nine clusters according to their nucleotide sequence similarities, 5 of which can be further divided into subclusters; 5 genomes do not cluster with other phages. The sequence diversity between genomes within a cluster varies greatly; for example, the 6 genomes in Cluster D share more than 97.5% average nucleotide similarity with one another. In contrast, similarity between the 2 genomes in Cluster I is barely detectable by diagonal plot analysis. In total, 6858 predicted open-reading frames have been grouped into 1523 phamilies (phams) of related sequences, 46% of which possess only a single member. Only 18.8% of the phams have sequence similarity to non-mycobacteriophage database entries, and fewer than 10% of all phams can be assigned functions based on database searching or synteny. Genome clustering facilitates the identification of genes that are in greatest genetic flux and are more likely to have been exchanged horizontally in relatively recent evolutionary time. Although mycobacteriophage genes exhibit a smaller average size than genes of their host (205 residues compared with 315), phage genes in higher flux average only 100 amino acids, suggesting that the primary units of genetic exchange correspond to single protein domains.     

Reconstruction of metabolic pathways for the cattle genome

Background  

Metabolic reconstruction of microbial, plant and animal genomes is a necessary step toward understanding the evolutionary origins of metabolism and species-specific adaptive traits. The aims of this study were to reconstruct conserved metabolic pathways in the cattle genome and to identify metabolic pathways with missing genes and proteins. The MetaCyc database and PathwayTools software suite were chosen for this work because they are widely used and easy to implement.     

Whole genome-based phylogenetic analysis of free-living microorganisms.

A phylogenetic 'tree of life' has been constructed based on the observed presence and absence of families of protein-encoding genes observed in 11 complete genomes of free-living microorganisms. Past attempts to reconstruct the evolutionary relation-ships of microorganisms have been limited to sets of genes rather than complete genomes. Despite apparent rampant lateral gene transfer among microorganisms, these results indicate a single robust underlying evolutionary history for these organisms. Broadly, the tree produced is very similar to the small subunit rRNA tree although several additional phylogenetic relationships appear to be resolved, including the relationship of Archaeoglobus to the methanogens studied. This result is in contrast to notions that a robust phylogenetic reconstruction of microorganisms is impossible due to their genomes being composed of an incomprehensible amalgam of genes with complicated histories and suggests that this style of genome-wide phylogenetic analysis could become an important method for studying the ancient diversification of life on Earth. Analyses using informational and operational subsets of the genes showed that this 'tree of life' is not dependent on the phylogenetically more consistent informational genes.     

Evolution of nodal and nodal‐related genes and the putative composition of the heterodimers that trigger the nodal pathway in vertebrates

Nodal is a signaling molecule that belongs to the transforming growth factor‐β superfamily that plays key roles during the early stages of development of animals. In vertebrates Nodal forms an heterodimer with a GDF1/3 protein to activate the Nodal pathway. Vertebrates have a paralog of nodal in their genomes labeled Nodal‐related, but the evolutionary history of these genes is a matter of debate, mainly because of the presence of a variable numbers of genes in the vertebrate genomes sequenced so far. Thus, the goal of this study was to investigate the evolutionary history of the Nodal and Nodal‐related genes with an emphasis in tracking changes in the number of genes among vertebrates. Our results show the presence of two gene lineages (Nodal and Nodal‐related) that can be traced back to the ancestor of jawed vertebrates. These lineages have undergone processes of differential retention and lineage‐specific expansions. Our results imply that Nodal and Nodal‐related duplicated at the latest in the ancestor of gnathostomes, and they still retain a significant level of functional redundancy. By comparing the evolution of the Nodal/Nodal‐related with GDF1/3 gene family, it is possible to infer that there are several types of heterodimers that can trigger the Nodal pathway among vertebrates.     

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.