Comparative Genomic Analysis of Members of the Genera Methanosphaera and Methanobrevibacter Reveals Distinct Clades with Specific Potential Metabolic Functions.
The microbial ecology of the intestinal tract of humans and animals has been subject to extensive research in recent years, and it is becoming increasingly evident that commensal intestinal microbes have a strong impact on host physiology and well-being. One microbial group, the methanogenic archaea (methanogens), has been of particular interest in this regard. Methanogens have not only been implicated in greenhouse gas emission from livestock animals [1, 2] but several studies have also linked these microorganisms to specific bodyweight phenotypes in humans  and to differences in feed conversion efficiency in ruminants [4, 5].
The majority of the rumen and intestinal methanogens appear to belong to two of the seven known orders, the Methanobacteriales  and the recently described Methanomassiliicoccales . Two genera of the order Methanobacteriales, Methanobrevibacter and Methanosphaera, comprise a large proportion of intestinal and rumen methanogenic archaea [8, 9]. Species of the genus Methanobrevibacter have been isolated from not only a wide range of different environments, including invertebrate and vertebrate guts, but also non-host-associated environments. The majority of the members of the genus Methanobrevibacter grow primarily hydrogenotrophically , using C[O.sub.2] and [H.sub.2] as substrates. A few strains and/or species have been shown by cultivation-independent methods to be the predominant in specific intestinal environments, such as Mbb. smithii from the human gut  or Mbb. ruminantium, and Mbb. gottschalkii from the rumen of sheep and cows . The former three species have gained considerable attention in the analysis of the human and rumen microbiome, but there are also Mbb. isolates that are less well investigated by cultivation-independent studies, for example, Mbb. thaueri strain DSM [11995.sup.T] and Mbb. woesei strain DSM [11976.sup.T]. Mbb. thaueri strain DSM [11995.sup.T] was isolated from cow feces while Mbb. woesei strain DSM [11976.sup.T] had been isolated from goose feces . The 16S rRNA gene of either of the two species shares the highest sequence identity to Mbb. smithii strain DSM [861.sup.T], Mbb. gottschalkii strain DSM [11977.sup.T], and Mbb. millerae strain DSM [16643.sup.T] [12-15]. Both species have been detected in a few cultivation-independent studies on fecal and rumen samples (for Mbb. thaueri, see [16, 17] and for Mbb. woesei, see [18-20]), but besides their original description, there is very little information on the physiology and ecology of these two species.
The genus Methanosphaera is less well characterized than the genus Methanobrevibacter, and only few species have been isolated or detected by cultivation-independent methods. Two species, Msp. stadtmanae DSM [3091.sup.T] and Msp. cuniculi DSM [4103.sup.T], are currently the only formally described Methanosphaera-type species [21, 22]. The first genome-sequenced Methanosphaera species, Msp. stadtmanae [DSM3091.sup.T], was isolated from human feces and has been shown to be restricted to growth on methanol and hydrogen under in vitro conditions . This substrate restriction can to some extent be considered to be a typical trait of the genus Methanosphaera, but its genetic basis remained poorly understood until the genome of Msp. stadtmanae strain DSM [3091.sup.T] was analyzed . Comparative genome analysis revealed that most genes for biosynthesis of the molybdopterin cofactor (Moco), the cofactor of formylmethanofuran dehydrogenase (fmd), are absent from the genome of this methanogen. This explained why Msp. stadtmanae strain [DSM3091.sup.T] is not able to grow by disproportionation of methanol nor of C[O.sub.2] and hydrogen as substrates. Moreover did the analysis reveal that the genome encodes four homologues of each subunit of the methanol: coenzyme-M methyltransferase (mtaABC) . The MtaABC proteins form the key enzyme complex for the utilization of methanol and had before primarily been detected in species of the order Methanosarcinales.
The second genome of a genus Methanosphaera representative, Methanosphaera sp. WGK6, has only recently become available. This methanogen's genome sequence is to large extent nearly identical to that of Msp. stadtmanae DSM [3091.sup.T], but some of the few genomic differences result in phenotypic differences between the two species. This regards mainly two genes, putative alcohol (walC) and aldehyde (walD) dehydrogenases . These two genes were apparently acquired via horizontal gene transfer and seem to enable Msp. sp. WGK6 to utilize ethanol as substrate for methanogenesis (in addition to its ability to grow on methanol and hydrogen) . Homologues of walC and walD genes have also been detected in genomes of other methanogens but have not been investigated systematically.
Here, draft genomes of Mbb. thaueri strain DSM [11995.sup.T], Mbb. woesei strain DSM [11979.sup.T], and Msp. cuniculi strain DSM [4103.sup.T] (a draft genome of M. cuniculi DSM [4103.sup.T] of similar quality was also recently published by Gilmore et al. ) are presented and compared to those of all other currently available Methanosphaera- and Methanobrevibactertype strain genomes. The comparative analyses aim at investigating and identifying some of the distinctive features of the Methanobrevibacter and Methanosphaera genomes that allow differentiating the two closely related genera. The major emphasis was set on key genes that may have large impact on the overall physiology of the analyzed species, such as those involved in substrate utilization and nitrogen and acetate assimilation.
2. Materials and Methods
2.1. Cultivation of Microorganism and DNA Extraction. Mbb. thaueri DSM [11995.sup.T], Mbb. woesei DSM [11979.sup.T], and Msp. cuniculi DSM [4103.sup.T] were obtained from the Deutsche Sammlung von Mikroorganismen und Zellkulturen (DSMZ), Braunschweig, Germany. Genomic DNA was ordered by the DSMZ (Braunschweig) or was isolated using the MasterPure complete DNA purification kit (Epicentre, Madison, WI, USA).
2.2. Genome Sequencing. Extracted DNA was used to generate Illumina-shotgun libraries (Nextera_XT) according to the manufacturer's protocol (Illumina, San Diego, CA, USA). Sequencing was conducted using a MiSeq and Miseq Reagent kit v3 (2 x 300 bp paired end) as recommended by the manufacturer (Illumina). Trimmomatic 0.32  was used to filter low-quality reads and for clipping of adapter contaminations. The assembly was performed with the SPAdes genome assembler software 3.11 . Coverages were determined using QualiMap version 2.1 [27-29], and automatic annotation was performed using the software tool PROKKA , and data analysis was partly performed using the IMG/ER system (Integrated Microbial Genomes & Microbiomes) . The quality and the completeness of the draft genomes have been validated with CheckM .
2.3. Multilocus Sequence Analysis. For multilocus sequence analysis (MLSA), total protein sequences from 19 genomes were extracted from the corresponding GenBank files using cds_extractor.pl v0.6  and used for downstream analysis with an in-house pipeline at the Goettingen Genomics Laboratory . In detail, proteinortho version 5 (default specification: blast = blastp v2.2.24, E value = 1e - 10, alg.conn. = 0.1, coverage = 0.5, percent_identity = 50, adaptive_similarity=0.95, inc_pairs=1, inc_singles= 1, selfblast=1, and unambiguous = 0)  was used to generate clusters of orthologue groups, inparalogues were removed, and MUSCLE  used to align the remaining sequences and poorly aligned positions were automatically filtered from the alignments using Gblocks . A maximum-likelihood tree from 574 orthologues was inferred with 500 bootstraps with RAxML . A phylogenetic tree was inferred with neighbour joining and 500 bootstraps. PO_2_MLSA.py is available at GitHub (https://github.com/jvollme). Visualization was done using Proteinortho results and DNAPlotter . The Venn diagram of the different clades was plotted using p[O.sub.2]group_stats (v0.1.1) .
2.3.1. Analysis of 16S rRNA Genes. Aligned 16S rRNA gene sequences were selected from the ARB-compatible rumen and intestinal methanogen database (RIM-DB) [14, 40]. Aligned sequences were exported in PHYLIP format to construct phylogenetic trees using all available base positions. Maximum-likelihood phylogenetic trees based on aligned archaeal 16S rRNA gene sequences were generated using RAxML version 7.3.0 . Unless stated otherwise, the parameters "-m GTRGAMMA -# 500 -f a-x 2 -p 2" were used.
2.4. Nucleotide Sequence Accession Number. The annotated genomes of Mbb. thaueri DSM [11995.sup.T] and Mbb. woesei DSM [11979.sup.T] have been deposited at DDBJ/EMBL/GenBank under the accession MZGS00000000 and MZGU00000000, respectively. The versions described in this paper are versions MZGS01000000 and MZGU01000000, respectively. The annotated genome of Msp. cuniculi DSM [4103.sup.T] has been deposited at DDBJ/EMBL/GenBank under the accession LWMS00000000. The version described in this paper is version LWMS01000000.
3.1. General Genome Features. The Methanosphaera cuniculi strain DSM [4103.sup.T] draft genome has been assembled into 48 contigs, with a N50 of 111,976 bp. The GC content of the draft genome is 28%, which is similar to the G + C content of the other two published Methanosphaera genome sequences. The genomes of Mbb. thaueri strain DSM [11995.sup.T] and Mbb. woesei strain DSM [11979.sup.T] were assembled into 39 contigs and 10 contigs with N50 of 178.217 bp and 240.239 bp, respectively. Plasmids were not detected from the assembled contigs in either of the two draft genomes. General features of the of Mbb. thaueri DSM [11995.sup.T], Mbb. woesei strain DSM [11979.sup.T], and Methanosphaera cuniculi strain DSM [4103.sup.T] genomes and comparison with other Methanobrevibacter and Methanosphaera strains are shown in Table 1, and results of CheckM analysis (including all genomes used in the analysis) are shown in Table S1. Circular representation of the three genomes is shown in supporting Figures S1 and S2.
3.2. Comparative Analysis of Methanobrevibacter Genomes. Sixteen different Methanobrevibacter and three Methanosphaera genomes were included for the comparative analyses. The Methanobrevibacter strains had been isolated from a range of different intestinal environments, for example, bovine and ovine rumen [12, 41], human gut , and termite hindguts [43, 44], but include also the genome of the non-host-associated Mbb. arboriphilus strain DSM [1125.sup.T] [45, 46]. Results of the single- and multilocus analyses are shown in Figure 1. Both analyses revealed clear separation between the two genera.
The MLSA also supports distinct clades within the genus Methanobrevibacter (Figure 1). There is limited association of the observed clades with a specific host origin, for example, some species from the rumen form distinct clades while others cluster with Methanobrevibacter species isolated from other hosts. A BLAST-based analysis of each clade identifies genes specific for clades and genes that are shared between them (Figure 2). Among the clade-specific genes of clade 1 are catalase, nitrogenase, and also a tRNA-specific endonuclease (VapC), but most other clade-specific genes of clade 1 and the other clades are without functional annotation and await further characterization.
It was also investigated whether certain metabolic traits could be associated with specific clades as pointed out in Figure 1. This concerns mainly the metabolism of alcohols and acetate assimilation. Carbon monoxide/acetyl-CoA-synthetase complex (CODH-ACS) is--in addition to other enzymes--important for acetate assimilation, and its presence is one prerequisite for autotrophic growth. CODH-ACS is found in non-host-associated and autotrophic methanogens, for example, Methanothermobacter species [47-49], and is detected in four out of five clade 1 species, but gene homologues of the enzyme appear to be absent from the genomes of methanogens in clades 2-4.
The potential utilization of alcohols, specifically methanol and ethanol, is also less clearly distributed among the 19 genomes (and the four different clades) and there are currently no known Methanobrevibacter species that are restricted to growth on [H.sub.2] and methanol/ethanol . Genes for methanol utilization of the Methanosphaera/Methanosarcina type (mtaABC) are distributed broadly within the genus Methanobrevibacter as outlined in Figure 1 and have also already been reported for some strains [46,50]. Physiological characterizations of the strains that harbor mtaABC genes (Mbb. arboriphilus strains DSM [1125.sup.T] and ANOR1, Mbb. smithii strain DSM [861.sup.T], Mbb. wolinii strain DSM [11976.sup.T]) is still outstanding to confirm that these species are capable of growth on methanol/methanol-hydrogen.
The utilization of ethanol by Methanobrevibacter species has only recently gained additional attention and it has been shown that Mbb. sp. AbM4 is capable of growth in the absence of hydrogen but in presence of methanol/ethanol [51, 52]. Details regarding the exact metabolism remain speculative, but the walC and walD genes recently identified in Msp. WGK6  are also present in Mbb. sp. AbM4 and in the closely related Mbb. boviskoreani strain DSM [25824.sup.T], Mbb. olleyae strain DSM [16632.sup.T], and Mbb. wolinii strain [DSM11976.sup.T] conferring this metabolic trait potentially to the entire clade. Mbb. olleyae strain DSM [16632.sup.T] is currently the only Mbb. species outside this clade that is harboring the walCD genes.
Only few additional genome sequences of Methanobrevibacter isolates other than for Mbb. smithii are currently available, but the clade-distinguishing features appear to be conserved in these isolates as well (see Table S2 for general genome features and Table S3 for shared genes).
3.3. Shared and Distinctive Genome Features in the Genomes of Methanosphaera Species. The currently known Methanosphaera species are characterized by their limited ability to use substrates for methanogenesis, for example, Msp. stadtmanae and Msp. cuniculi are able to utilize only methanol and hydrogen for growth, Msp. sp. WKG6 can also utilize ethanol, but none of the species is capable of utilizing hydrogen and C[O.sub.2] or formate like organisms from closely related genera. The analysis of the Msp. stadtmanae genome helped explaining some of the distinctive physiological features of this methanogen, but comparative genome analysis with other closely related Methanosphaera and Methanobrevibacter species was not possible at that time due to the lack of sequenced genomes. After more than a decade, several Methanosphaera and Methanobrevibacter genomes have become available [24, 41, 46, 50, 53-58] and comparative analyses allow determining whether certain traits are species specific or shared by members of either or both genera. Comparative genome analysis reveals that some of the observations made for the Msp. stadtmanae genome appear to be consistent for all three available Methanosphaera genomes. Overall, more than 1000 genes are shared among all three Methanosphaera species (Figure 3), but the comparison with Methanobrevibacter genomes shows that several key genes are missing. Notably, this concerns the lack of genes for molybdopterin cofactor biosynthesis, energy-converting hydrogenase A (ehaA-Q), formate dehydrogenase (fdh), and methyl-coenzyme M reductase isoenzyme 1 mcrABCG (only genes encoding isoenzyme II are present (mrtABDG)). Despite the presence of formylmethanofuran dehydrogenase genes, it is not possible for Methanosphaera species to produce a functional enzyme due to the absence of the molybdopterin cofactor. The absence of functional formylmethanofuran dehydrogenase leads to the inability of these methanogens to grow hydrogenotrophically or to disproportionate methanol. The absence of Moco biosynthesis also seems to be specific to Methanosphaera species as all of the analyzed Methanobrevibacter genomes are capable of synthesizing this cofactor.
One of the noteworthy differences between the known Methanosphaera species is the capability of Msp. WGK6 to utilize ethanol as an electron donor for methanogenesis. The recently reported alcohol dehydrogenase and aldehyde dehydrogenase (NL43_02835, NL43_02830) in the genome of Methanosphaera sp. WGK6 were not detected in the Msp. stadtmanae genome (as previously reported by Hoedt et al., 2016) or the Msp. cuniculi genome , indicating that this may be a trait of Msp. WGK6 that is shared with other lessrelated methanogens.
The presented comparative genome analysis reveals insights into the genome content of different Methanobrevibacter species and how it compares to that of closely related Methanosphaera species. The analysis suggests that phenotypic variation among Methanobrevibacter strains may be larger than previously assumed. The comparative analysis indicates the presence of genes for methanol (in four Mbb. genomes) and ethanol utilization (in four Mbb. genomes) in the sixteen analyzed Methanobrevibacter genomes, with the Mbb. wolinii genome harboring both walCD and mtaABC gene homologues. More physiological analysis will be required to determine if the detected genes are conferring the ability to utilize ethanol/methanol to all these strains. The ability of Methanobrevibacter strains to utilize alcohols may not have been investigated in an in-depth manner previously as this trait is more typically associated with other orders of methanogens, but some studies suggest that the genes may in some cases also not be functional or may have unknown different functions, for example, growth of Mbb. smithii on methanol was not detected despite the presence of mtaABC genes (growth on methanol but not methanol and hydrogen was tested) . It is noteworthy that all three sequenced strains of the Mbb. wolinii clade harbor the genes for ethanol utilization, which could point to a specific ecological role of the species in this clade; however, there are only few studies that have detected significant numbers of either of these three species in a natural environment [5, 60] making it currently difficult to determine potential cooccurrences or specific syntrophic interactions with other microorganisms.
The presence of genes that could contribute to autotrophic growth, for example, nitrogenase and CODH-ACS in clade 1 containing Mbb. species from termites and in Mbb. arboriphilus strain, has been reported previously
, and these genes appear to be absent from species clades 2 to 4. It can only be speculated why Mbb. species of clade 1 did not undergo the same loss of these key genes like their counterparts in the clades that have primarily been isolated from the vertebrate intestinal tract. However, auxotrophy of some Mbb. species may be the result of a close symbiotic interaction with other microorganisms (and potentially the host) that provide favorable growth conditions as well as ammonium and acetate for the methanogen. It is also noteworthy that Mbb. curvatus in the "autotrophic clade" has nitrogenase genes, but apparently, no CODH complex genes and may therefore require externally supplied acetate. This could represent an intermediate stage between the potentially autotrophic strains and the other clades and/or may be an adaptation to a specific niche.
In addition to differences between species, there is currently only little information on strain diversity within a species but studies indicate that there may be considerable differences as shown for Mbb. smithii and Mbb. arboriphilus [46, 61]. However, undertaking pan-genome approaches, such as the one for Mbb. smithii by Hansen et al. , requires substantial cultivation efforts as there are only few strains available for each of the described methanogen species. As development of sequencing technologies continues to advance, it may also become feasible to obtain highquality closed genomes from low amounts of starting DNA or through metagenomic approaches. Having such genomes will allow greater certainty in determining the presence and absence of specific genome features and will also enable detection of small genomic differences between strains that may go unnoticed in draft genomes.
Lastly, the results of our analyses corroborate the hypothesis that the absence of molydopterin cofactor biosynthesis is a characteristic trait shared by members of the genus Methanosphaera, while all sequenced Methanobrevibacter genomes seem to encode genes for Moco biosynthesis. However, it also needs to be considered that Methanosphaera species encode the genes for formylmethanofuran dehydrogenase and other genes required for hydrogenotrophic or methylotrophic methanogenesis [62, 63]. The presence of these genes may indicate that Methanosphaera species may still be able to utilize the aforementioned methanogenic pathways, if molybdopterin cofactor is present and taken up by the methanogen (as suggested by Fricke et al. 2006). Metatranscriptional profiling could be used to determine, if fmd and other genes for hydrogenotrophic growth are expressed at all in Methanosphaera or specifically in the absence and presence of Moco or Moco-synthesizing microorganisms. It needs to be emphasized that the analyses of this study are based on the currently available three genomes of the three Methanosphaera species that have been isolated from three different host species [23-25]. It provides strong evidence that the physiology will be similar for other Msp. species in other environments, but it can also not be ruled out that other Msp. species exist that harbor genes for Moco biosynthesis and that grow hydrogenotrophically or methylotrophically. The ability of Msp. sp. WGK6 to utilize ethanol does indicate that more phenotypic variation exists among Methanosphaera species and that further isolation and characterization of new Msp. species is necessary.
The presented study includes genomes of all currently available types of strain (and other selected isolate) of Methanobrevibacter and Methanosphaera genomes and allows comprehensive insights into the genera Methanobrevibacter and Methanosphaera. The analyses reveal that distinct clades within the genus Methanobrevibacter exist and that the lack of molybdopterin cofactor biosynthesis may be a specific trait for the genus Methanosphaera. Additional isolates and further physiological and genomic analyses will be required to determine if division of the genus Methanobrevibacter in more than one genus could be justified. The primary use of the type of strains (and other isolates) for the analysis in this study does warrant access of the scientific community to most of the analyzed Methanobrevibacter and Methanosphaera isolates and will facilitate testing of the predicted physiological phenotypes.
The data used to the support the findings of this study are availabe in the supplemental information or -as in the case of the genomes- are downloadable using the accession numbers provided in Section 2.4.
Please Note: Illustration(s) are not available due to copyright restrictions.
Conflicts of Interest
The authors declare that they have no conflicts of interest. Acknowledgments
This work was funded by Temasek Life Sciences Laboratory. The authors thank Melanie Heinemann for the technical support.
Supplementary 1. Figure S1: circular display of the Mbb. thaueri strain DSM [11995.sup.T]and Mbb. woesei strain DSM [11979.sup.T] genomes and comparison with other Methanobrevibacter genomes. The Mbb. thaueri strain DSM [11995.sup.T] and the Mbb. woesei strain DSM [11979.sup.T] genomes are shown in panels A and B, respectively. The genes encoded by the leading and the lagging strand of the respective species are shown in circles 1 and 2, and rRNA (pink) and tRNA (green) encoding genes are shown in circle 3. The presence of orthologous genes (red: high similarity; orange: medium similarity; yellow: low similarity; see color code below) is indicated for the other Methanobrevibacter genomes (circles 4 to 19) in comparison to that of the respective Methanobrevibacter genome. The two innermost plots represent the GC content and the GC skew (circle 20 and 21). Color code according to E values of the blastp analysis performed using Proteinortho4.26. Gray: 1e - 20 to 1; light yellow: 1e - 21 to 1e - 50; gold: 1e - 51 to 1e - 90; light orange: 1e - 91 to 1e - 100; orange: 1e - 101 to 1e - 120; red: > 1e - 120.
Supplementary 2. Figure S2: circular display of the Msp. cuniculi strain DSM [4103.sup.T] genome and comparison with the genomes of Msp. stadtmanae strain DSM [3091.sup.T] and Msp. sp. WGK6. The genes encoded by the leading and the lagging strand of Msp. cuniculi strain DSM [4103.sup.T] are shown in circles 1 and 2, and rRNA (pink) and tRNA (green) encoding genes are shown in circle 3. The presence of orthologous genes (red: high similarity; orange: medium similarity; yellow: low similarity; see color code below) is indicated for the genomes of Msp. stadtmanae strain DSM [3091.sup.T] and Msp. sp. WGK6 (circles 4 and 5 in comparison to the Msp. cuniculi strain DSM [4103.sup.T] genome). The two innermost plots represent the GC content and the GC skew (circles 6 and 7). Color code according to E values of the blastp analysis performed using Proteinortho4.26. Gray: 1e - 20 to 1; light yellow: 1e - 21 to 1e - 50; gold: 1e - 51 to 1e - 90; light orange: 1e - 91 to 1e -100; orange: 1e - 101 to 1e - 120; red: > 1e -120.
Supplementary 3. Table S1: CheckM analysis to determine the completeness of genomes. Table S1A shows results from analysis of genomes from the main manuscript while Table S1B lists the results for genomes mainly mentioned in the supporting information. Table S2: general features of additional Methanobrevibacter genomes. General features of genomes mainly mentioned in the supporting information. Table S3: table of shared and unshared Methanobrevibacter and Methanosphaera genes. Highlighted in colors corresponding to those of Figure 1 are genes for CODH-ACS complex, nitrogenase, methanol, and ethanol utilization. Annotations were taken from the Mbb. arboriphilus DSM1125t genome.
 S. Kittelmann, C. S. Pinares-Patino, H. Seedorf et al., "Two different bacterial community types are linked with the low-methane emission trait in sheep," PLoS One, vol. 9, no. 7, article e103171, 2014.
 W. Shi, C. D. Moon, S. C. Leahy et al., "Methane yield phenotypes linked to differential gene expression in the sheep rumen microbiome," Genome Research, vol. 24, no. 9, pp. 1517-1525, 2014.
 J. K. Goodrich, J. L. Waters, A. C. Poole et al., "Human genetics shape the gut microbiome," Cell, vol. 159, no. 4, pp. 789-799, 2014.
 C. A. Carberry, S. M. Waters, D. A. Kenny, and C. J. Creevey, "Rumen methanogenic genotypes differ in abundance according to host residual feed intake phenotype and diet type," Applied and Environmental Microbiology, vol. 80, no. 2, pp. 586-594, 2014.
 M. Zhou, E. Hernandez-Sanabria, and L. L. Guan, "Assessment of the microbial ecology of ruminal methanogens in cattle with different feed efficiencies," Applied and Environmental Microbiology, vol. 75, no. 20, pp. 6524-6533, 2009.
 A. S. Bonin, D. R. Boone et al.M. Dworkin, S. Falkow, E. Rosenberg et al., "The Order Methanobacteriales," in The Prokaryotes, vol. 3pp. 231-243, Springer, 2006.
 T. Iino, H. Tamaki, S. Tamazawa et al., "Candidatus Methanogranum caenicola: a novel methanogen from the anaerobic digested sludge, and proposal of Methanomassiliicoccaceae fam. nov. and Methanomassiliicoccales ord. nov., for a methanogenic lineage of the class Thermoplasmata," Microbes and Environments, vol. 28, no. 2, pp. 244-250, 2013.
 P. B. Eckburg, E. M. Bik, C. N. Bernstein et al., "Diversity of the human intestinal microbial flora," Science, vol. 308, no. 5728, pp. 1635-1638, 2005.
 G. Henderson, F. Cox, S. Ganesh et al., "Rumen microbial community composition varies with diet and host, but a core microbiome is found across a wide geographical range," Scientific Reports, vol. 5, no. 1, 2015.
 T. L. Miller, "Methanobrevibacter, Bergey's manual of systematics of archaea and bacteria," The Archaea and Deeply Branching and Phototrophic Bacteria, vol. 1, pp. 218-226, 2001.
 H. Seedorf, S. Kittelmann, and P. H. Janssen, "Few highly abundant operational taxonomic units dominate within rumen methanogenic archaeal species in New Zealand sheep and cattle," Applied and Environmental Microbiology, vol. 81, no. 3, pp. 986-995, 2015.
 T. L. Miller and C. Lin, "Description of Methanobrevibacter gottschalkii sp. nov., Methanobrevibacter thaueri sp. nov., Methanobrevibacter woesei sp. nov. and Methanobrevibacter wolinii sp. nov," International Journal of Systematic and Evolutionary Microbiology, vol. 52, pp. 819-822, 2002.
 P. H. Janssen and M. Kirs, "Structure of the archaeal community of the rumen," Applied and Environmental Microbiology, vol. 74, no. 12, pp. 3619-3625, 2008.
 H. Seedorf, S. Kittelmann, G. Henderson, and P. H. Janssen, "RIM-DB: a taxonomic framework for community structure analysis of methanogenic archaea from the rumen and other intestinal environments," PeerJ, vol. 2, article e494, 2014.
 B. St-Pierre and A.-D. G. Wright, "Molecular analysis of methanogenic archaea in the forestomach of the alpaca (Vicugna pacos)," BMC Microbiology, vol. 12, no. 1, p. 1, 2012.
 S. L. Ishaq, M. A. Sundset, J. Crouse, and A.-D. G. Wright, "High-throughput DNA sequencing of the moose rumen from different geographical locations reveals a core ruminal methanogenic archaeal diversity and a differential ciliate protozoal diversity," Microbial Genomics, vol. 1, no. 4, article e000034, 2015.
 A. Salgado-Flores, L. H. Hagen, S. L. Ishaq et al., "Rumen and cecum microbiomes in reindeer (Rangifer tarandus tarandus) are changed in response to a lichen diet and may affect enteric methane emissions," PLoS One, vol. 11, no. 5, article e0155213, 2016.
 N. Just, P. B. Lecours, M. Marcoux-Voiselle et al., "Archaeal characterization of bioaerosols from cage-housed and floorhoused poultry operations," Canadian Journal of Microbiology, vol. 59, no. 1, pp. 46-50, 2013.
 F. Li, G. Henderson, X. Sun, F. Cox, P. H. Janssen, and L. L. Guan, "Taxonomic assessment of rumen microbiota using total RNA and targeted amplicon sequencing approaches," Frontiers in Microbiology, vol. 7, 2016.
 S. Saengkerdsub, R. C. Anderson, H. H. Wilkinson, W. K. Kim, D. J. Nisbet, and S. C. Ricke, "Identification and quantification of methanogenic archaea in adult chicken ceca," Applied and Environmental Microbiology, vol. 73, no. 1, pp. 353-356, 2006.
 B. Biavati, M. Vasta, and J. G. Ferry, "Isolation and characterization of Methanosphaera cuniculi sp. nov," Applied and Environmental Microbiology, vol. 54, no. 3, pp. 768-771, 1988.
 T. L. Miller and M. J. Wolin, "Methanosphaera stadtmaniae gen. nov., sp. nov.: a species that forms methane by reducing methanol with hydrogen," Archives of Microbiology, vol. 141, no. 2, pp. 116-122, 1985.
 W. F. Fricke, H. Seedorf, A. Henne et al., "The genome sequence of Methanosphaera stadtmanae reveals why this human intestinal archaeon is restricted to methanol and [H.sub.2] for methane formation and ATP synthesis," Journal of Bacteriology, vol. 188, no. 2, pp. 642-658, 2005.
 E. C. Hoedt, P. O. Cuiv, P. N. Evans et al., "Differences downunder: alcohol-fueled methanogenesis by archaea present in Australian macropodids," The ISME Journal, vol. 10, no. 10, pp. 2376-2388, 2016.
 S. P. Gilmore, J. K. Henske, J. A. Sexton et al., "Genomic analysis of methanogenic archaea reveals a shift towards energy conservation," BMC Genomics, vol. 18, no. 1, p. 639, 2017.
 A. M. Bolger, M. Lohse, and B. Usadel, "Trimmomatic: a flexible trimmer for Illumina sequence data," Bioinformatics, vol. 30, no. 15, pp. 2114-2120, 2014.
 A. Bankevich, S. Nurk, D. Antipov et al., "SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing," Journal of Computational Biology, vol. 19, no. 5, pp. 455-477, 2012.
 K. Okonechnikov, A. Conesa, and F. Garcia-Alcalde, "Qualimap 2: advanced multi-sample quality control for high-throughput sequencing data," Bioinformatics, vol. 32, no. 2, pp. 292-294, 2016.
 F. Garcia-Alcalde, K. Okonechnikov, J. Carbonell et al., "Qualimap: evaluating next-generation sequencing alignment data," Bioinformatics, vol. 28, no. 20, pp. 2678-2679, 2012.
 T. Seemann, "Prokka: rapid prokaryotic genome annotation," Bioinformatics, vol. 30, no. 14, pp. 2068-2069, 2014.
 V. M. Markowitz, K. Mavromatis, N. N. Ivanova, I.-M. A. Chen, K. Chu, and N. C. Kyrpides, "IMG ER: a system for microbial genome annotation expert review and curation," Bioinformatics, vol. 25, no. 17, pp. 2271-2278, 2009.
 D. H. Parks, M. Imelfort, C. T. Skennerton, P. Hugenholtz, and G. W. Tyson, "CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes," Genome Research, vol. 25, no. 7, pp. 1043-1055, 2015.
 A. Leimbach, "Bovine E. coli mastitis comparative genomics edition," Bac-genomics-scripts, 2016.
 F. S. Haack, A. Poehlein, C. Kroger et al., "Molecular keys to the Janthinobacterium and Duganella spp. interaction with the plant pathogen Fusarium graminearum," Frontiers in Microbiology, vol. 7, no. 1668, 2016.
 M. Lechner, S. Findeifi, L. Steiner, M. Marz, P. F. Stadler, and S. J. Prohaska, "Proteinortho: detection of (co-) orthologs in large-scale analysis," BMC Bioinformatics, vol. 12, no. 1, p. 124, 2011.
 R. C. Edgar, "MUSCLE: multiple sequence alignment with high accuracy and high throughput," Nucleic Acids Research, vol. 32, no. 5, pp. 1792-1797, 2004.
 J. Castresana, "Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis," Molecular Biology and Evolution, vol. 17, no. 4, pp. 540-552, 2000.
 A. Stamatakis, "RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models," Bioinformatics, vol. 22, no. 21, pp. 2688-2690, 2006.
 T. Carver, N. Thomson, A. Bleasby, M. Berriman, and J. Parkhill, "DNAPlotter: circular and linear interactive genome visualization," Bioinformatics, vol. 25, no. 1, pp. 119-120, 2009.
 W. Ludwig, O. Strunk, R. Westram et al., "ARB: a software environment for sequence data," Nucleic Acids Research, vol. 32, no. 4, pp. 1363-1371, 2004.
 J.-H. Lee, M.-S. Rhee, S. Kumar et al., "Genome sequence of Methanobrevibacter sp. strain JH1, isolated from rumen of Korean native cattle," Genome Announcements, vol. 1,no. 1, pp. e00002-e00013, 2013.
 T. L. Miller and M. J. Wolin, "Enumeration of Methanobrevibacter smithii in human feces," Archives of Microbiology, vol. 131, no. 1, pp. 14-18, 1982.
 J. R. Leadbetter and J. A. Breznak, "Physiological ecology of Methanobrevibacter cuticularis sp. nov. and Methanobrevibacter curvatus sp. nov., isolated from the hindgut of the termite Reticulitermes flavipes," Applied and Environmental Microbiology, vol. 62, no. 10, pp. 3620-3631, 1996.
 J. R. Leadbetter, L. D. Crosby, and J. A. Breznak, "Methanobrevibacter filiformis sp. nov., a filamentous methanogen from termite hindguts," Archives of Microbiology, vol. 169, no. 4, pp. 287-292, 1998.
 J. G. Zeikus and D. L. Henning, "Methanobacterium arbophilicum sp. nov. an obligate anaerobe isolated from wetwood of living trees," Antonie Van Leeuwenhoek, vol. 41, no. 1, pp. 543-552, 1975.
 A. Poehlein, R. Daniel, and H. Seedorf, "The draft genome of the non-host-associated methanobrevibacter arboriphilus strain dh1 encodes a large repertoire of adhesin-like proteins," Archaea, vol. 2017, 9 pages, 2017.
 A.-K. Kaster, M. Goenrich, H. Seedorf et al., "More than 200 genes required for methane formation from [H.sub.2] and C[O.sub.2] and energy conservation are present in Methanothermobacter marburgensis and Methanothermobacter thermautotrophicus," Archaea, vol. 2011, 23 pages, 2011.
 H. Liesegang, A. K. Kaster, A. Wiezer et al., "Complete genome sequence of Methanothermobacter marburgensis, a methanoarchaeon model organism," Journal of Bacteriology, vol. 192, no. 21, pp. 5850-5851, 2010.
 D. R. Smith, L. A. Doucette-Stamm, C. Deloughery et al., "Complete genome sequence of Methanobacterium thermoautotrophicum deltaH: functional analysis and comparative genomics," Journal of Bacteriology, vol. 179, no. 22, pp. 7135-7155, 1997.
 B. S. Samuel, E. E. Hansen, J. K. Manchester et al., "Genomic and metabolic adaptations of Methanobrevibacter smithii to the human gut," Proceedings of the National Academy of Sciences, vol. 104, no. 25, pp. 10643-10648, 2007.
 S. C. Leahy, W. J. Kelly, R. S. Ronimus, N. Wedlock, E. Altermann, and G. T. Attwood, "Genome sequencing of rumen bacteria and archaea and its application to methane mitigation strategies," Animal, vol. 7, no. s2, pp. 235-243, 2013.
 M. R. Weimar, J. Cheung, D. Dey et al., "Development of multi-well plate methods using pure cultures of methanogens to identify new inhibitors for suppressing ruminant methane emissions," Applied and Environmental Microbiology, vol. 83, no. 15, pp. e00396-e00317, 2017.
 S. C. Leahy, W. J. Kelly, D. Li et al., "The complete genome sequence of Methanobrevibacter sp. AbM4," Standards in Genomic Sciences, vol. 8, no. 2, pp. 215-227, 2013.
 S. Khelaifia, M. Garibal, C. Robert, D. Raoult, and M. Drancourt, "Draft genome sequence of a human-associated isolate of Methanobrevibacter arboriphilicus, the lowest-G+C-content archaeon," Genome Announcements, vol. 2, no. 1, pp. e01181-e01113, 2014.
 W. J. Kelly, D. Li, S. C. Lambie et al., "Draft genome sequence of the rumen methanogen Methanobrevibacter olleyae YLM1," Genome Announcements, vol. 4, no. 2, pp. e00232-e00216, 2016.
 A. Poehlein and H. Seedorf, "Draft genome sequences of Methanobrevibacter curvatus DSM11111, Methanobrevibacter cuticularis DSM11139, Methanobrevibacter filiformis DSM11501, and Methanobrevibacter oralis DSM7256," Genome Announcements, vol. 4, no. 3, pp. e00617-e00616, 2016.
 S. Khelaifia, M. Garibal, C. Robert, D. Raoult, and M. Drancourt, "Draft genome sequencing of Methanobrevibacter oralis strain JMR01, isolated from the human intestinal microbiota," Genome Announcements, vol. 2, no. 1, pp. e00073-e00014, 2014.
 S. C. Leahy, W. J. Kelly, E. Altermann et al., "The genome sequence of the rumen methanogen Methanobrevibacter ruminantium reveals new possibilities for controlling ruminant methane emissions," PLoS One, vol. 5, no. 1, article e8926, 2010.
 T. L. Miller, M. J. Wolin, E. Conway de Macario, and A. J. Macario, "Isolation of Methanobrevibacter smithii from human feces," Applied and Environmental Microbiology, vol. 43, no. 1, pp. 227-232, 1982.
 S. Federici, F. Miragoli, V. Pisacane, A. Rebecchi, L. Morelli, and M. L. Callegari, "Archaeal microbiota population in piglet feces shifts in response to weaning: Methanobrevibacter smithii is replaced with Methanobrevibacter boviskoreani," FEMS Microbiology Letters, vol. 362, no. 10, 2015.
 E. E. Hansen, C. A. Lozupone, F. E. Rey et al., "Pan-genome of the dominant human gut-associated archaeon, Methanobrevibacter smithii, studied in twins," Proceedings of the National Academy of Sciences, vol. 108, Supplement_1, pp. 4599-4606, 2011.
 R. K. Thauer, "Biochemistry of methanogenesis: a tribute to Marjory Stephenson: 1998 Marjory Stephenson Prize lecture," Microbiology, vol. 144, no. 9, pp. 2377-2406, 1998.
 R. K. Thauer, A.-K. Kaster, H. Seedorf, W. Buckel, and R. Hedderich, "Methanogenic archaea: ecologically relevant differences in energy conservation," Nature Reviews Microbiology, vol. 6, no. 8, pp. 579-591, 2008.
Anja Poehlein, (1) Dominik Schneider, (1) Melissa Soh (ID), (2) Rolf Daniel, (1) and Henning Seedorf (ID) (2, 3)
(1) Genomic and Applied Microbiology and Gottingen Genomics Laboratory, Georg-August University, Gottingen, Germany
(2) Research Link, Temasek Life Sciences Laboratory, Singapore 117604
(3) National University of Singapore, Department of Biological Sciences, Singapore 117604
Correspondence should be addressed to Henning Seedorf; firstname.lastname@example.org
Received 30 April 2018; Accepted 2 July 2018; Published 19 August 2018
Academic Editor: Stefan Spring
Caption: Figure 1: Single- and multilocus sequence analysis of Methanobrevibacter and Methanosphaera species. A maximum likelihood tree (left) of 16 Methanobrevibacter and three Methanosphaera genomes was inferred with 500 bootstraps with RAxML and visualized with Dendroscope. Phylogeny of Methanobrevibacter and Methanosphaera based on the 16S rRNA gene is shown on the right. The tree was resampled 500 times, and only bootstrap values [greater than or equal to] 70% are shown. The 16S rRNA tree was rooted with five Methanobacterium sequences. The scale bar indicates 0.10 inferred nucleotide substitutions per position. Red-colored dots indicate the presence of mtaABC genes in the species/clade, greencolored dots indicate the presence of walB and walC gene homologues (potential utilization of ethanol), blue-colored dots indicate the presence of nitrogenase genes, and black-colored dots indicate the presence of carbon monoxide/acetyl CO-DH genes (see also Table S3 for details).
Caption: Figure 2: Pan/core genome analysis of four different Methanobrevibacter clades. Venn diagram showing the numbers of orthologous genes (OGs) in the core, dispensable, and specific genome of compared strains. Ortholog detection was done with the Proteinortho software (blastp) with a similarity cut-off of 50% and an E value of 1e - 10. The total numbers of genes and paralogs are depicted under the corresponding species name. Open-reading frames that were classified as pseudogenes were not included in this analysis. See also Table S3 for details on shared genes.
Caption: Figure 3: Pan/core genome analysis of three different Methanosphaera species. Venn diagram showing the numbers of orthologous genes (OGs) in the core, dispensable, and specific genome of compared strains. Ortholog detection was done with the Proteinortho software (blastp) with a similarity cut-off of 50% and an E value of 1e - 10. The total numbers of genes and paralogs are depicted under the corresponding species name. Open-reading frames that were classified as pseudo genes were not included in this analysis. See also Table S3 for details on shared genes.
Table 1: General features of Methanobrevibacter and Methanosphaera genomes. Size (bp) GC Coding content (%) percentage (%) Methanobrevibacter 2,606,143 26.99 69.94 filiformis DSM [11501.sup.T] Methanobrevibacter 2,414,608 25.72 70.65 curvatus DSM [11111.sup.T] Methanobrevibacter 2,608,702 26.79 68.12 cuticularis DSM [11139.sup.T] Methanobrevibacter 2,445,031 25.44 74.96 arboriphilus DSM [1125.sup.T] Methanobrevibacter 2,221,072 25.53 74.02 arboriphilus ANOR1 Methanobrevibacter 2,937,203 32.64 78.12 ruminantium DSM [1093.sup.T] Methanobrevibacter 2,122,444 26.87 76.85 olleyae DSM [16632.sup.T] Methanobrevibacter 2,725,667 36.54 89.32 millerae DSM [16643.sup.T] Methanobrevibacter 2,243,115 36.87 87.52 thaueri DSM [11995.sup.T] Methanobrevibacter 1,879,371 30.02 88.01 gottschalkii DSM [11977.sup.T] Methanobrevibacter 2,110,861 27.73 84.49 oralis DSM [7256.sup.T] Methanobrevibacter 1,853,160 31.03 90.31 smithii DSM [861.sup.T] Methanobrevibacter 1,543,150 29.90 89.61 woesei DSM [11979.sup.T] Methanobrevibacter 2,041,814 24.21 75.94 wolinii DSM [11976.sup.T] Methanobrevibacter 2,045,801 28.98 78.00 boviskoreani DSM [25824.sup.T] Methanobrevibacter 1,998,189 29.04 76.7 sp. AbM4 Methanosphaera sp. 1,729,155 27.70 78.05 WGK6 Methanosphaera 1,767,403 27.63 84.10 stadtmanae DSM [3091.sup.T] Methanosphaera 1,881,497 28.08 81.36 cuniculi DSM [4103.sup.T] CDS Genes rRNA tRNA (pseudo) (pseudo) Methanobrevibacter 1933 1965 3 29 filiformis DSM [11501.sup.T] Methanobrevibacter 1969 2004 4 31 curvatus DSM [11111.sup.T] Methanobrevibacter 2061 2094 3 30 cuticularis DSM [11139.sup.T] Methanobrevibacter 1963 2005 5 37 arboriphilus DSM [1125.sup.T] Methanobrevibacter 1993 2038 7 35 arboriphilus ANOR1 Methanobrevibacter 2217 2283 (5) 8 53 ruminantium DSM [1093.sup.T] Methanobrevibacter 1813 1854 4 33 olleyae DSM [16632.sup.T] Methanobrevibacter 2383 2467 4 77 millerae DSM [16643.sup.T] Methanobrevibacter 2138 2171 2 31 thaueri DSM [11995.sup.T] Methanobrevibacter 1845 1889 7 34 gottschalkii DSM [11977.sup.T] Methanobrevibacter 2036 1994 9 32 oralis DSM [7256.sup.T] Methanobrevibacter 1795 1841 8 36 smithii DSM [861.sup.T] Methanobrevibacter 1581 1614 2 31 woesei DSM [11979.sup.T] Methanobrevibacter 1700 1747 8 36 wolinii DSM [11976.sup.T] Methanobrevibacter 1756 1799 4 36 boviskoreani DSM [25824.sup.T] Methanobrevibacter 1671 1716 49 36 sp. AbM4 Methanosphaera sp. 1456 1616 (114) 3 42 WGK6 Methanosphaera 1534 1592 (2) 12 42 stadtmanae DSM [3091.sup.T] Methanosphaera 1585 1629 5 39 cuniculi DSM [4103.sup.T] Accession Contigs/ CRISPR Original scaffolds region strain ID Methanobrevibacter LWMT00000000 237 6 RFM-3 filiformis DSM [11501.sup.T] Methanobrevibacter LWMV00000000 187 5 RFM-2 curvatus DSM [11111.sup.T] Methanobrevibacter LWMW00000000 120 6 RFM-1 cuticularis DSM [11139.sup.T] Methanobrevibacter JXMW00000000 40 9 DH1 arboriphilus DSM [1125.sup.T] Methanobrevibacter CBVX000000000 5 3 ANOR1 arboriphilus ANOR1 Methanobrevibacter CP001719 1 4 Ml ruminantium DSM [1093.sup.T] Methanobrevibacter FOTLOOOOOOOO 49 4 KM1H5-1P olleyae DSM [16632.sup.T] Methanobrevibacter CP011266 1 1 ZA-10 millerae DSM [16643.sup.T] Methanobrevibacter MZGS00000000 39 4 CW thaueri DSM [11995.sup.T] Methanobrevibacter FOAKOOOOOOOO 19 NA HO gottschalkii DSM [11977.sup.T] Methanobrevibacter LWMU00000000 99 2 ZR oralis DSM [7256.sup.T] Methanobrevibacter CP000678 1 1 PS smithii DSM [861.sup.T] Methanobrevibacter MZGU00000000 10 2 GS woesei DSM [11979.sup.T] Methanobrevibacter JHWX00000000 32 2 SH wolinii DSM [11976.sup.T] Methanobrevibacter BAGX00000000 54 2 JH1 boviskoreani DSM [25824.sup.T] Methanobrevibacter CP004050 1 6 N.A. sp. AbM4 Methanosphaera sp. JRWK00000000 37 1 N.A. WGK6 Methanosphaera CP000102 1 4 MCB-3 stadtmanae DSM [3091.sup.T] Methanosphaera LWMS00000000 48 0 1R-7 cuniculi DSM [4103.sup.T]
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Poehlein, Anja; Schneider, Dominik; Soh, Melissa; Daniel, Rolf; Seedorf, Henning|
|Date:||Jan 1, 2018|
|Previous Article:||Crystal Structures of the Putative Isocitrate Dehydrogenase from Sulfolobus tokodaii Strain 7 in the Apo and [NADP.sup.+]-Bound Forms.|
|Next Article:||Characterization of Microbial Communities in Pilot-Scale Constructed Wetlands with Salicornia for Treatment of Marine Aquaculture Effluents.|