The complete genome sequence of the rumen methanogen Methanosarcina barkeri CM1

Abstract

Methanosarcina species are the most metabolically versatile of the methanogenic Archaea and can obtain energy for growth by producing methane via the hydrogenotrophic, acetoclastic or methylotrophic pathways. Methanosarcina barkeri CM1 was isolated from the rumen of a New Zealand Friesian cow grazing a ryegrass/clover pasture, and its genome has been sequenced to provide information on the phylogenetic diversity of rumen methanogens with a view to developing technologies for methane mitigation. The 4.5 Mb chromosome has an average G + C content of 39 %, and encodes 3523 protein-coding genes, but has no plasmid or prophage sequences. The gene content is very similar to that of M. barkeri Fusaro which was isolated from freshwater sediment. CM1 has a full complement of genes for all three methanogenesis pathways, but its genome shows many differences from those of other sequenced rumen methanogens. Consequently strategies to mitigate ruminant methane need to include information on the different methanogens that occur in the rumen.

Introduction

Ruminants are foregut fermenters and have evolved an efficient digestive system in which microbes ferment plant fibre and provide fermentation end-products and other nutrients for growth of the animal [1]. A variety of methanogens can be found in the rumen [2] and ruminant derived methane (CH₄) accounts for about one quarter of all anthropogenic CH₄ emissions [3], and is implicated as a driver of global climate change. In terms of their metabolism the rumen methanogens fall into three groups, hydrogenotrophs ( Methanobrevibacter , Methanomicrobium and Methanobacterium spp) which convert hydrogen and/or formate to CH₄, methylotrophs ( Methanosphaera spp and members of the order Methanomassiliicoccales ) which produce CH₄ from methyl compounds such as methanol and methylamines, and acetoclastic methanogens ( Methanosarcina ) which can utilise acetate to produce CH₄ in addition to the hydrogenotrophic and methylotrophic pathways. Obtaining representative genome sequences from each of the above organisms will be important to understanding the metabolic capacity of these archaea and how they contribute to rumen fermentation processes. Currently, genome sequences are available for five rumen methanogens including strains of Methanobrevibacter ruminantium [4], M. boviskoreani [5, 6], Methanobacterium formicicum [7] and Thermoplasmatales archaeon BRNA1 [NCBI Reference Sequence: NC_020892.1]. Development of strategies to reduce CH₄ emissions from farmed ruminant animals are currently being investigated with methanogen genome sequence information used to inform mitigation strategies based on vaccines and small-molecule inhibitors [8, 9]. Here we present the genome sequence from a rumen acetoclastic methanogen, Methanosarcina barkeri CM1.

Organism information

Classification and features

Methanosarcina sp. CM1 was isolated from the rumen of a New Zealand Friesian cow grazing a ryegrass/clover pasture [10]. CM1 grew as large cell aggregates in broth culture and showed the characteristic morphology associated with Methanosarcina barkeri [11] (Fig. 1). It was described as non-motile, and able to grow and produce methane from H₂/CO₂, acetate, methanol and methylamines. Growth occurred between 30° and 45 °C, and at pH 5.0 to 7.4. Rumen fluid was required for growth. The 16S rRNA gene from CM1 is 99 % similar to that of the Methanosarcina barkeri type strain MS (DSM 800) (Fig. 2) which was isolated from a sewage sludge digester [12, 13], and as such CM1 can be considered as a strain of M. barkeri . M. barkeri is found at high densities in anaerobic digesters and anoxic marine and freshwater sediments, but there have been several reports describing Methanosarcina from the rumen although these organisms were not characterized [14, 15]. In addition, non-rumen strains of M. barkeri have been used in co-culture studies with rumen anaerobic fungi [16] and ciliate protozoa [17, 18]. Characteristics of Methanosarcina barkeri CM1 are shown in Table 1 and Additional file 1.

Fig. 1

Morphology of Methanosarcina barkeri CM1. Micrograph showing aggregates of Methanosarcina barkeri CM1 cells captured with sectional depth scanning using an Olympus Fluoview FV1000D Spectral laser confocal scanning inverted microscope, with an UPLSAPO 60X oil objective (1.35 NA). Olympus Fluoview 10-ASW software was used to view fluorescent signals and to generate images. Emission at 635 nm wavelength shows methyl green stain incorporated into nucleic acids within cells and cell clusters. Bar is 10 μm

Fig. 2

Phylogenetic tree showing the position of CM1 relative to type strains of other Methanosarcina species. The strains and their corresponding accession numbers are shown. The evolutionary history was inferred using the Neighbor-Joining method [45] with Methanobrevibacter ruminantium used as an outgroup. The optimal tree with the sum of branch length = 0.43777587 is shown. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) are shown next to the branches [46]. The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the Kimura 2-parameter method [47] and are in the units of the number of base substitutions per site. The rate variation among sites was modeled with a gamma distribution (shape parameter = 1). The analysis involved 13 nucleotide sequences. All positions containing gaps and missing data were eliminated. There were a total of 1081 positions in the final dataset. Evolutionary analyses were conducted in MEGA5 [48]

Table 1 Classification and general features of Methanosarcina barkeri CM1

Genome sequencing information

Genome project history

Methanosarcina barkeri CM1 was selected for genome sequencing on the basis of its phylogenetic position relative to other methanogens isolated from the rumen. A summary of the genome project information is shown in Table 2.

Table 2 Project information

Growth conditions and genomic DNA preparation

Methanosarcina barkeri CM1 was grown in BY medium [19] with added SL10 Trace Elements solution (1 ml added l⁻¹ [20], 20 mM sodium acetate, 60 mM sodium formate and Vitamin 10 solution (0.1 ml added to 10 ml culture before inoculation) [4]. H₂ was supplied as the energy source by pumping the culture vessels to 180 kPa over pressure with an 80:20 mixture of H₂:CO₂. Genomic DNA was extracted from freshly grown cells using a modified version of a liquid N₂ and grinding method as described previously [6].

Genome sequencing and assembly

The complete genome sequence of CM1 was determined using pyrosequencing of 3Kb mate paired-end sequence libraries using a 454 GS FLX platform with Titanium chemistry (Macrogen, Korea). Pyrosequencing reads provided 97× coverage of the genome and were assembled using the Newbler assembler version 2.7 (Roche 454 Life Sciences, USA). The Newbler assembly resulted in 85 contigs across 9 scaffolds. Gap closure was managed using the Staden package [21] and gaps were closed using additional Sanger sequencing by standard and inverse PCR based techniques. In addition, CM1 genomic DNA was sequenced using the Illumina HiSeq 2000 platform (Beijing Genomics Institute, China) which provided 223× genome coverage. Illumina reads were assembled using the Spades assembler version 3.0 [22] and combined with the Newbler assembly using the Staden package. Assembly validation was confirmed by pulsed-field gel electrophoresis as described previously [6].

Genome annotation

The procedure for genome annotation was as described previously for Methanobrevibacter sp. [4, 6], and the CM1 genome sequence was prepared for NCBI submission using Sequin. The adenine residue of the start codon of the Cdc6-1 replication initiation protein A (MCM1_0001) gene was chosen as the first base for the CM1 genome. The nucleotide sequence of the Methanosarcina barkeri CM1 chromosome has been deposited in Genbank under accession number CP008746.

Genome properties

The genome of Methanosarcina barkeri CM1 consists of a single 4,501,171 basepair (bp) circular chromosome with an average G + C content of 39 %. A total of 3656 genes were predicted, 3523 of which were protein-coding genes, representing 70 % of the total genome sequence. A COG category was assigned to 2267 of the protein-coding genes. The properties and statistics of the genome are summarized in Tables 3 and 4. As with the other sequenced Methanosarcina strains CM1 has dual origins of replication (MCM1_001 and MCM1_3593, 95 kb apart) surrounded by conserved genes [23]. The CM1 genome has neither plasmid nor prophage sequences, but does contain three clusters of CRISPR genes associated with CRISPR repeat regions, and three type I restriction/modification systems.

Table 3 Genome statistics

Table 4 Number of genes associated with the 25 general COG functional categories

Insights from the genome sequence

The genome of Methanosarcina barkeri CM1 is compared with genomes of other sequenced methanogens from the genus Methanosarcina in Table 5. Overall, the gene content of the CM1 genome is very similar to that of Methanosarcina barkeri Fusaro, but gene organization shows very little synteny.

Table 5 Genomes of Methanosarcina species from various anaerobic environments

Methanogenesis

Methanosarcina species are the most metabolically versatile of the methanogenic archaea [24] and can obtain energy for growth by producing methane via three different pathways (Fig. 3). Methane can be derived from the reduction of CO₂ with hydrogen (hydrogenotrophic pathway), from the methyl group of acetate (acetoclastic pathway), or from the methyl group of methanol, methylamines or methylthiols (methylotrophic pathway). Each pathway culminates in the transfer of a methyl group to coenzyme M and the subsequent reduction to methane. The bioenergetics of aceticlastic methanogens have been recently reviewed [25, 26], and a metabolic reconstruction presented for M. barkeri Fusaro [27]. There is evidence that the genes essential to both the acetoclastic and methyoltrophic pathways were horizontally acquired during evolution of the Methanosarcinaceae [28–30]. Although acetoclastic methanogenesis contributes approximately two-thirds of the methane in the biosphere [31], acetate is not metabolized to methane to any significant extent in the rumen [2]. CM1 has a full complement of genes for all three methanogenesis pathways and as with other Methanosarcina species several genes are present as multiple copies [32]. Unlike many hydrogenotrophic methanogens, CM1 does not have the gene for [Fe]-hydrogenase dehydrogenase (hmd), or the genes that encode methyl coenzyme M reductase II (mrt), however it does have genes for formate dehydrogenase (MCM1_3047-3048) although CM1 and other M. barkeri strains are unable to use formate [10, 23]. The pathway for coenzyme M biosynthesis differs from that found in other sequenced rumen methanogens which belong to the order Methanobacteriales [33]. Consequently methanogen inhibitors targeting coenzyme M biosynthesis would not be expected to work against all rumen methanogen species. M. barkeri is the organism in which the 22nd amino acid (pyrrolysine) was discovered during examination of the methyltransferases required for methane formation from methylamines [34]. Biosynthesis of this amino acid requires specialized enzymes together with a specific aminoacyl-tRNA synthetase [35], and the genes encoding these (pylSBCD) are found in CM1 (MCM1_2535-2538).

Fig. 3

The three methanogenesis pathways inferred from the genome of Methanosarcina barkeri CM1

Cell envelope

The majority of rumen methanogens belong to the family Methanobacteriaceae and have a characteristic pseudomurein-containing cell wall [4]. The cell surface of Methanosarcina sp. is different and electron microscopy shows a protein S-layer surrounding the cytoplasmic membrane. A major S-layer protein has been identified in three Methanosarcina species and used to define a family of proteins [36, 37]. All paralogs possess signal peptides and one or two DUF1608 (Pfam accession PF07752) domains. Both M. barkeri Fusaro and CM1 possess nine proteins containing this domain. Gene expression studies show that a single DUF1608 domain-containing protein is abundantly expressed in both M. acetivorans (MA0829) and M. barkeri (Mbar_A1758), and is among the most highly expressed of all proteins in the cell [37]. A similar protein is predicted from the CM1 genome (MCM1_2018, 84 % amino acid identity to the product of Mbar_A1758), and likely to be the major S-layer protein for this strain. The DUF1608 domain from MA0829 has been crystallised and was used to propose an elegant model of the Methanosarcina S-layer structure [38]. Methanosarcina cells can form large aggregates (Fig. 1) mediated by the production of methanochondroitin, a polymer composed of uronic acid and N-acetylgalactosamine residues [39]. While the steps in methanochondroitin biosynthesis have been determined, the genes involved have yet to be identified. These aggregates are observed to separate into single cells and CM1 encodes eleven proteins with disaggregatase-related domains (Pfam accession PF08480). The CM1 genome has four oligosaccharyl transferase genes, three of which (MCM1_1841-1843) are at the end of a large gene cluster that contains 14 glycosyl transferases and is likely to be involved in polysaccharide biosynthesis (MCM1_1841-1889). The fourth oligosaccharyl transferase is associated with a smaller gene cluster containing glycosyl transferases, methyltransferases and transporters (MCM1_2113-2123). A third cluster of polysaccharide biosynthesis genes is found at MCM1_2831-2857. CM1 also encodes a secreted protein (MCM_2974) containing a glycoside hydrolase family 18 (chitinase) domain that is not found in M. barkeri Fusaro which may be involved in mediating interaction with rumen anaerobic fungi. Like many other archaea, CM1 has an identifiable archaella (archaeal flagella) operon (FlaB-FlaJ, MCM1_1947-1953), together with a cluster of chemotaxis genes (MCM1_3655-3662) [40, 41]. However, motility has never been observed in any Methanosarcina species and thus the function of these genes remains unknown.

Methanosarcina barkeri has been reported to fix nitrogen [42] and sets of nitrogenase genes are found in Methanosarcina genomes. CM1 contains two different nif operons comprising nitrogenase and nitrogenase cofactor biosynthesis genes that match to those reported from M. barkeri strain 227 [43]. These are a molybdenum- and iron-containing nitrogenase (MCM1_2924-2930) and a vanadium- and iron-containing nitrogenase (MCM1_1063-1072). However, it does not have the genes for the third type, the iron-only nitrogenase that is found in M. acetivorans and M. barkeri Fusaro. Electron micrographs of M. barkeri log phase cells [11] show the presence of numerous electron-dense granules in the cytoplasm. In M. thermophila similar granules were found to contain glycogen [44], and CM1 has several genes predicted to encode the enzymes necessary for the biosynthesis and degradation of this reserve polysaccharide. CM1 does not have the genes for gas vesicle biosynthesis that are found in the M. barkeri Fusaro genome [23], but it does have genes for the two-subunit acetyl-CoA synthetase (MCM1_1658 and 2708) that have been lost from the Fusaro strain.

Conclusion

The genome of Methanosarcina barkeri CM1 is very similar to that of the freshwater sediment isolate M. barkeri Fusaro, but markedly different from the dominant rumen methanogens, most of which are members of the family Methanobacteriaceae . CM1 has a much larger genome and its sequence provides new insights into the metabolic versatility of rumen methanogens. With its ability to use three different methanogenesis pathways, M. barkeri appears to be a generalist able to occupy a range of different environments but is not particularly at home in the rumen. This is in contrast to the more specialised rumen methanogens, such as the Methanobrevibacter species which dominate the rumen environment. Analysis of the methanogenesis pathway and the cell envelope have been important for the design of methane mitigation strategies targeting rumen methanogens, but differences highlighted from the CM1 genome stress the need to include information from all rumen methanogens in the design of mitigation approaches.