- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Article Processing Charges ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Recently Accepted Articles ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Genetics Research International
Volume 2013 (2013), Article ID 980191, 7 pages
Molecular Genetic Diversity and Quantitation of Methanogen in Ruminal Fluid of Buffalo (Bubalus bubalis) Fed Ration (Wheat Straw and Concentrate Mixture Diet)
1P. G. Department of Genetics, ARIBAS, New V V Nagar, Anand, Gujarat 388121, India
2Department of Animal Biotechnology, College of Veterinary Science and Animal Husbandry, Anand Agricultural University, Anand, Gujarat 388 001, India
3Animal Nutrition Research Station, AAU, Anand, Gujarat 388001, India
4Department of Microbiology, Christ College, Rajkot, Gujarat 360 005, India
Received 27 February 2013; Revised 13 May 2013; Accepted 16 May 2013
Academic Editor: Jorge H. Leitão
Copyright © 2013 K. M. Singh et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
High roughage diet causes more methane emissions; however, the total methanogen abundance is not influenced by roughage proportion. Technologies to reduce methane emissions are lacking, and development of inhibitors and vaccines that mitigate rumen-derived methane by targeting methanogens relies on present knowledge of the methanogens. In this work, we have investigated molecular diversity of rumen methanogens of Surti buffalo. DNA from rumen fluid was extracted, and 16S rRNA encoding genes were amplified using methanogen specific primer to generate 16S rDNA clone libraries. Seventy-six clones were randomly selected and analysed by RFLP resulting in 21 operational taxonomic units (OTUs). BLAST analysis with available sequences in database revealed sequences of 13 OTUs (55 clones) showing similarity with Methanomicrobium sp, 3 OTUs (15 clones) with Methanobrevibacter sp. The remaining 5 OTUs (6 clones) belonged to uncultured archaea. The phylogenetic analysis indicated that methanogenic communities found in the library were clustered in the order of Methanomicrobiales (18 OTUs) and Methanobacteriales (3 OTUs). The population of Methanomicrobiales, Methanobacteriales, and Methanococcales were also observed, accounting for 1.94%, 0.72%, and 0.47% of total archaea, respectively.
Methanogens, members of the domain Archaea, fall within the kingdom euryarchaeota . They are obligate anaerobes and can be unmistakably differentiated from other organisms since they all produce methane as a major catabolic end-product . The most common species of methanogens isolated from the rumen are Methanobrevibacter ruminantium, Methanomicrobium mobile, Methanobacterium formicicum, Methanobacterium smithii, Methanobacterium olleyae, Methanobacterium bryantii, Methanosarcina barkeri, and Methanoculleus olentangyi [3–9]. The diversity of archaea found in the rumen has been recently reviewed [10–12].
Interest in methanogens from ruminants has resulted from the role of methane in global warming and from the fact that enteric methane emission is a major source of greenhouse gas in agriculture sector. Currently, India possesses the world’s largest livestock population of 485 million, which accounts for 13% of the global livestock population (Intergovernmental Panel on Climate Change, 2001). It has 57% of the world’s buffalo and 16% of the cattle population. Buffalo contributes to 42% of the total methane emission by livestock in India . Reducing enteric methane emissions has been identified as one way of lowering global methane emissions. However, the effectiveness of any strategy that reduces greenhouse gas emissions and also increases production or nutritional efficiency will likely depend upon having an understanding of the numbers and/or distribution of methanogen species among ruminant livestock. In the present study, diversity analysis of methanogen consortium in ruminal fluid of buffalo (Bubalus bubalis) fed wheat straw and compound concentrate mixture was carried out.
2. Materials and Methods
2.1. Sampling and DNA Extraction
The permission of the Committee for the Purpose of Control and Supervision of Experiments on Animals (CPCSEA) was obtained prior to initiation of the study. The experiments were carried out on three young Surti buffaloes, approximately two years of age, which were reared at the Department of Animal Nutrition, College of Veterinary Science and A.H., Anand. All the animals were maintained under uniform feeding regime for minimum 30 days. The diet comprised wheat straw and compound concentrate mixture (20% crude protein, 65% total digestible nutrients). Approximately 500 mL of rumen fluid was collected 4 hrs after feeding . About 100 mL rumen fluid was passed through four layers of cheese cloth to remove particulate matter. Remaining rumen fluid was stored at −80°C for further study. Total DNA was extracted separately by using a commercially available kit according to the manufacturer’s instructions (QIAGEN Stool kit; QIAGEN, CA) and finally pooled all the DNA extracts. The total DNA was used as a template in PCR experiments to amplify 16S rDNA.
2.2. PCR Primers and Amplification
The 16S r DNA was amplified by PCR using metagenomic DNA as template and primers meth 86f (5′-ACAGGCCTAACACATGCAAGTC-3′) and meth1360r (5′-AGGGCGG(AT)GTGTACAAGGC-3′) . A total of 25 μL of reaction mixtures consisted of 10 pmol of each primer, 75–100 ng of template DNA, and 12.5 μL of master mix (Fermentas, UK). The reaction mixtures were subjected to initial denaturation at 95°C for 5 minutes following 30 cycles of each denaturation at 94°C for 1 minute, annealing at 58°C for 45 seconds, and extension at 72°C for 1 minute with final extension at 72°C for 10 minutes using thermal cycler (ABI, USA). The anticipated product of approximately 1.2 kb was purified using Qiagen DNA Gel Extraction Kits (QIAGEN, CA) in accordance with the manufacturer’s instructions.
2.3. Cloning, Screening, and RFLP
The purified PCR products were cloned in pTZR57T plasmid using the InstaT/A cloning kit (Fermentas, UK) according to the manufacturer’s protocol and transformed into Escherichia coli DH5α. The blue and white colonies were screened on Luria-Bertani plates with ampicillin (100 μg/mL), X-gal (20 mg/mL), and IPTG (100 mM). A total of 76 clones were examined from pooled metagenomes by colony PCR using vector specific primers M13f (5′-GTAAAACGACGGCCAG-3′) and M13r (5′-CAGGAAACAGCTATGAC-3′). Plasmids were extracted and amplified by PCR with primers M13f and M13r. The total 25 μL volume of the reaction mixture contained 10 pmol of each primer, 12.5 μL of master mix (Fermentas, UK), and 2.5 μL plasmid DNA. PCR conditions were denaturation for 5 min at 94°C, 30 cycles of denaturation at 94°C for 1 min, annealing at 50°C for 1 min, extension at 72°C for 1 min, and a final extension at 72°C for 10 min. Aliquots (10 μL) of all successfully recovered clones were digested with 0.5 U of HaeIII restriction endonuclease at 37°C overnight as described by . The digested DNA was visualized after electrophoresis on 1.5% agarose gels. One of each RFLP profile was assumed to indicate a different sequence and each different RFLP profile was purified and used for sequencing.
2.4. Sequences and Phylogenetic Analysis
Sequencing performed for all the clones with an ABI Prism 310 Genetic Analyser (Applied Biosystems Inc., CA) using BigDye Terminator (version 3.1) cycle sequencing kit (ABI, USA) at the Animal Biotechnology Laboratory, AAU, Anand, Gujarat, India.
All reference sequences were obtained from the Genbank/EMBL . Sequences from the current study were mainly trimmed to remove low-quality base calls from the start and end of DNA sequences and further analysed by the CHECK_CHIMERA program  to remove any chimera rDNA clone. The similarity searches for sequences were carried out by BLAST (http://www.ncbi.nlm.nih.gov/BLAST/Blast.cgi), and alignment was perform using CLUSTAL W (http://www.ebi.ac.uk/Tools/clustalw2/index.html).
The phylogenetic relatedness was estimated using the neighbor-joining method . All positions containing gaps and missing data were eliminated from the dataset (complete deletion option). One thousand bootstrap replications were performed to place the confidence estimates on the major groups resolved in the tree. The bootstrap consensus tree inferred from 1,000 replicates represents the evolutionary history of the sequences analyzed . The phylogenetic analysis was carried out using MEGA software version 4.0 .
2.5. Real-Time PCR Analysis
Plasmid DNA containing the total archea, methanomicrobiales and methanobacteriales specific sequences, used as the standard DNA in real-time PCR, was obtained by PCR cloning with the specific primer sets already described . After the confirmation of a single band of the correct size with respective pair of primers on an agarose gel, the PCR products were excised from the gel. The PCR products were purified using the Qiagen gel Purification Kit (Qiagen, CA) and then ligated into pTZR57T/A cloning vector (Fermentas, UK). The ligated products were transformed to competent E. coli DH5α cells by heat shock. Plasmids were purified from positive clones using a QIAprep spin miniprep kit (Qiagen, USA), and the plasmids containing the correct insert were screened out by PCR amplification with respective primer sets. Tenfold dilution series ranging from 108 to 10 copies were prepared for each target. Real-time PCR was performed with ABI system (ABI 7500). The Qiagen 2X SYBR Green master mix was used for PCR reaction. The optimal amplification conditions for each primer set were obtained with 10 pmol of each primer with the combination of annealing temperature and extension time as described by . The 10-fold dilution series of the standard plasmid for the respective target was run along with the corresponding samples in duplicate. The copy numbers of 16S rRNA genes of targeted methanogens per mL rumen fluid were calculated using the following equation: , where QM is the quantitative mean of the copy number, C is the DNA concentration of each sample, DV is dilution volume of extracted DNA, S is the DNA amount (ng) subjected to analysis, and V is the rumen fluid volume subjected to DNA extraction . In the reaction, nearly perfect linear regressions to 0.9995 and slope (−3.3 to −4.5) were obtained between threshold cycle and quantities of standard for all targets, and data generated from the reaction were used for further analysis.
2.6. Nucleotide Sequence Accession Numbers
The nucleotide sequences of the 16S RNA gene from the representative clones (based on RFLP) were deposited in NCBI under the Accession nos. HM566228–HM566248.
3.1. Comparison of Cloned Sequences with Sequences Deposited in Databases
A total of 76 clones from sample were screened and 21 sequences (OTUs) were generated based on PCR-RFLP patterns. Twenty-one sequences (OTUs) were subjected to similarity analysis using BLAST search  and SIMILARIY_RANK program . In the library, a Methanomicrobium like clone accounted for almost 72.36% of the clones (55 of 76 clones), and 17.9% of the cloned sequences (15 clones) were aligned to Methanobrevibacter sp with 96%–98% identity. 7.89% of the cloned sequences (6 clones) Also, belonged to uncultured archaeon with 94%–98% similarity (Table 1).
3.2. Phylogenetic Analysis of Sequences and Quantitation of Methanogens
The results of phylogenetic analysis of sequences are shown in Figure 1. Rumen methanogens formed a tree mainly divided into clusters named Methanomicrobiales-II and Methanosarcinales, Methanomicrobiales-I, and Methanosarcinales. In cluster Methanomicrobiales-I, ten OTUs (26 clones) were phylogenetically placed within genus Methanomicrobia with only one species Methanomicrobia mobile. However Methanomicrobiales-I shows to be more related with Methanobacteriales. In the cluster Methanomicrobiales-II, eight OTUs (43 clones) were also phylogenetically placed within genus Methanomicrobia and the sequences were 94–98% identical to those of Methanomicrobia mobile. In cluster Methanobacteriales, three OTUs (7 clones) were phylogenetically affiliated to genus Methanobrevibacter sp with 96%–98% similarity (Table 1). Total archaea were detected with copies per mL of ruminal fluid. The number of 16S rRNA gene copies of Methanomicrobiales, Methanobacteriales, and Methanococcales was detected with , , and copies per mL of ruminal fluid and accounting for 1.94%, 0.72%, and 0.47% of total archaea, respectively.
Methanogens have been classified into more than 100 species distributed by more than 20 genera . Interestingly, few methanogens have been isolated from the rumen. So far, cultured methanogens obtained from rumen have been assigned to Methanobrevibacter ruminantium , Methanobrevibacter olleyae , Methanomicrobium mobile [9, 26], and Methanoculleus olentangyi . Methanobacterium formicicum  as well as Methanosarcina spp have also been cultured from the rumen . However, in the present study, there were only a few genera identified in our library, which may be due to the ruminant host, diet, DNA extraction methods, or PCR primers used . In addition, it is possible that the DNA from some methanogens could not be extracted or the extracted DNA of some species was too low for amplification. PCR amplification bias, primers, and geological distribution could be other reasons for the true diversity of the composition of rumen methanogens. Ramos et al.  and Nadais et al.  have developed a method for PCR analysis to characterize any microbial population in anaerobic sludge blanket (UASB) reactor sample based on the 16S r RNA-encoding genes and clones analysis using frequently cutting enzymes.
Methanogens in the present study, PCR-retrieved methanogenic 16S rRNA gene libraries, were established from the rumen fluid. The results showed that the rRNA gene clones generated in this study exhibited a high degree of sequence similarity to two methanogenic genus, Methanomicrobium and Methanobrevibacter. The most predominant species of methanogens in the study were related to the genera Methanomicrobium. This is an agreement with our previous report  of rumen methanogens of buffalo fed green fodder and concentrates, which indicates that diet composition does not affect the methanogen population. Predominance of sequences highly similar to those of the Methanomicrobium genus in the rumen of Indian buffalo has also been reported [34, 35]. In addition, Shin et al.  and Tajima et al.  have also reported the dominance of Methanomicrobium sp. rumen of cows and found that Methanobrevibacter sp. was either not detected or was not the dominant clade affiliated to genus Methanobrevibacter sp. However, in our study only few sequences were affiliated to genus Methanobrevibacter. In contrast, Whitford et al.  found that Methanobrevibacter sp. was the dominant archaea (58.5%) in the cow rumen, where Methanosphaera sp. (26.8%) and Methanimicrococcus spp. (14.6%) were also abundant. Wright et al.  found that Methanobrevibacter sp. dominated in feedlot cattle in two geographic locations in Canada while sheep in Australia and Venezuela shared very similar archaeal communities. Tatsuoka et al.  and Denman et al.  have reported the diversity of the methyl-coenzyme M reductase (mcrA) gene in the rumen of cattle and found that most of sequences belonged to Methanobrevibacter sp.
This study reveals the phylogenetic diversity of the rumen methanogens in Surti buffalo rumen by analyzing methanogenic specific16S rRNA gene sequences in a culture-independent manner. Sixteen out of the 21 sequences (OTUs) were associated with two methanogenic genera (Table 1): Methanomicrobium (13 sequences: HM566229, HM566230, HM566232, HM566234, HM566238, HM566239, HM566241, HM566242, HM566243, HM566244, HM566245, HM566246, and HM566247) and Methanobrevibacter (3 Sequences: HM566228, HM566240 and HM566248) (Figure 1). The remaining 5 sequences (OTUs) belonged to uncultured methanogenic archaeon. BLAST search showed that all sequences shared a degree of similarity ranging from 90% to 99% with ruminal archaeon 16S rDNA sequences deposited in NCBI GenBank database (Table 1).
In conclusion, the knowledge of the ruminal methanogen community is of critical importance for the development of strategies to mitigate rumen methane production. This study has revealed a predominance of Methanomicrobium in the rumen of Surti buffalo and indicates the need to better understand the diet dependent association study of methanogen with methane emission. Indeed on the basis of the information of these predominant methanogen species existing in their specific niche, this research provides basic information for further studies on elucidating the symbiotic relationship of methanogens and other microbes, the mechanism of methanogenesis, and the regulation of methane emissions in ruminants. Reducing enteric methane emissions is likely to be one of the key mitigation strategies for the reduction of greenhouse gas emissions in the agricultural sector.
Conflict of Interests
The authors declare there is no conflict of interests.
Financial support provided by the Department of Biotechnology, Government. of India, New Delhi, to conduct the study reported here is acknowledged with respect and gratitude.
- C. R. Woese, O. Kandler, and M. L. Wheelis, “Towards a natural system of organisms: proposal for the domains Archaea, Bacteria, and Eucarya,” Proceedings of the National Academy of Sciences of the United States of America, vol. 87, no. 12, pp. 4576–4579, 1990.
- D. H. Bergey, “Themethanogens,” in Bergey'S Manual of Determinative Bacteriology Baltimore, J. G. Holt, Ed., pp. 719–725, Williams and Wilkins, 1994.
- L. Raskin, L. K. Poulsen, D. R. Noguera, B. E. Rittmann, and D. A. Stahl, “Quantification of methanogenic groups in anaerobic biological reactors by oligonucleotide probe hybridization,” Applied and Environmental Microbiology, vol. 60, no. 4, pp. 1241–1248, 1994.
- C. S. Stewart, et al., “The rumen bacteria,” in The Rumen Microbial Ecosystem, P. N. Hobson and C. S. Stewart, Eds., pp. 10–72, Blackie Academic and Professional, London, UK, 1997.
- M. J. Wolin, T. L. Miller, and C. S. Stewart, “Microbe-microbe interactions,” in The Rumen Microbial Ecosystem, P. N. Hobson and C. S. Stewart, Eds., pp. 467–491, Blackie Academic and Professional, New York, NY, USA, 1997.
- D. R. Lovley, R. C. Greening, and J. G. Ferry, “Rapidly growing rumen methanogenic organism that synthesizes coenzyme M and has a high affinity for formate,” Applied and Environmental Microbiology, vol. 48, no. 1, pp. 81–87, 1984.
- T. L. Miller, M. J. Wolin, Z. Hongxue, and M. P. Bryant, “Characteristics of methanogens isolated from bovine rumen,” Applied and Environmental Microbiology, vol. 51, no. 1, pp. 201–202, 1986.
- M. F. Whitford, R. M. Teather, and R. J. Forster, “Phylogenetic analysis of methanogens from the bovine rumen,” BMC Microbiology, vol. 1, article 5, 2001.
- S. Kumar, S. S. Dagar, and A. K. Puniya, “Isolation and characterization of methanogens from rumen of Murrah buffalo,” Annals of Microbiology, vol. 62, no. 1, pp. 345–350, 2012.
- 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.
- D. P. Morgavi, E. Forano, C. Martin, and C. J. Newbold, “Microbial ecosystem and methanogenesis in ruminants,” Animal, vol. 4, no. 7, pp. 1024–1036, 2010.
- S. E. Hook, A.-D. G. Wright, and B. W. McBride, “Methanogens: methane producers of the rumen and mitigation strategies,” Archaea, vol. 2010, Article ID 945785, 11 pages, 2010.
- A. Chhabra, K. R. Manjunath, S. Panigrahy, and J. S. Parihar, “Spatial pattern of methane emissions from Indian livestock,” Current Science, vol. 96, no. 5, pp. 683–689, 2009.
- S. Khampa, M. Wanapat, C. Wachirapakorn, N. Nontaso, and M. Wattiaux, “Effects of urea level and sodium DL-malate in concentrate containing high cassava chip on ruminal fermentation efficiency, microbial protein synthesis in lactating dairy cows raised under tropical condition,” Asian-Australasian Journal of Animal Sciences, vol. 19, no. 6, pp. 837–844, 2006.
- A.-D. G. Wright and C. Pimm, “Improved strategy for presumptive identification of methanogens using 16S riboprinting,” Journal of Microbiological Methods, vol. 55, no. 2, pp. 337–349, 2003.
- D. A. Benson, I. Karsch-Mizrachi, D. J. Lipman, J. Ostell, and D. L. Wheeler, “GenBank,” Nucleic Acids Research, vol. 35, no. 1, pp. D21–D25, 2007.
- B. L. Maidak, J. R. Cole, T. G. Lilburn et al., “The RDP-II (Ribosomal Database Project),” Nucleic Acids Research, vol. 29, no. 1, pp. 173–174, 2001.
- N. Saitou and M. Nei, “The neighbor-joining method: a new method for reconstructing phylogenetic trees,” Molecular Biology and Evolution, vol. 4, no. 4, pp. 406–425, 1987.
- J. Felsenstein, “Confidence limits on phylogenies: an approach using the bootstrap,” Evolution, vol. 39, pp. 783–791, 1985.
- K. Tamura, J. Dudley, M. Nei, and S. Kumar, “MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0,” Molecular Biology and Evolution, vol. 24, no. 8, pp. 1596–1599, 2007.
- Y. Yu, C. Lee, J. Kim, and S. Hwang, “Group-specific primer and probe sets to detect methanogenic communities using quantitative real-time polymerase chain reaction,” Biotechnology and Bioengineering, vol. 89, no. 6, pp. 670–679, 2005.
- K. M. Singh, A. K. Tripathi, P. R. Pandya, D. N. Rank, R. K. Kothari, and C. G. Joshi, “Dasytricha dominance in Surti buffalo rumen revealed by 18S rRNA sequences and real-time PCR assay,” Current Microbiology, vol. 63, no. 3, pp. 281–288, 2011.
- S. F. Altschul, T. L. Madden, A. A. Schäffer et al., “Gapped BLAST and PSI-BLAST: a new generation of protein database search programs,” Nucleic Acids Research, vol. 25, no. 17, pp. 3389–3402, 1997.
- G. M. Garrity, T. G. Lilburn, J. R. Cole, S. H. Harrison, J. Euzeby, and B. J. Tindall, “Taxonomic outline of the Bacteria and Archaea,” Part1: The Archaea, phyla Crenarchaeota and Euryarchaeota. Release 7.7. Michigan State University, 2007.
- P. H. Smith and R. E. Hungate, “Isolation and characterization of Methanobacterium ruminantium n. sp,” Journal of Bacteriology, vol. 75, no. 6, pp. 713–718, 1958.
- M. J. Paynter and R. E. Hungate, “Characterization of Methanobacterium mobilis, sp. n., isolated from the bovine rumen,” Journal of Bacteriology, vol. 95, no. 5, pp. 1943–1951, 1968.
- K. N. Joblin, “Methanogenic archaea,” in Methods in Gut Microbial Ecology For Ruminants, H. P. S. Makkar and C. S. McSweeney, Eds., pp. 47–53, Springer, Amsterdam, The Netherlands, 2005.
- R. A. Oppermann, W. O. Nelson, and R. E. Brown, “In vitro studies on methanogenic rumen bacteria,” Journal of Dairy Science, vol. 40, pp. 779–788, 1957.
- W. H. Beijer, “Methane fermentation in the rumen of cattle,” Nature, vol. 170, no. 4327, pp. 576–577, 1952.
- A.-D. G. Wright, C. H. Auckland, and D. H. Lynn, “Molecular diversity of methanogens in feedlot cattle from Ontario and Prince Edward Island, Canada,” Applied and Environmental Microbiology, vol. 73, no. 13, pp. 4206–4210, 2007.
- C. G. Ramos, A. M. Grilo, S. A. Sousa, M. L. Barbosa, H. Nadais, and J. H. Leitão, “A new methodology combining PCR, cloning, and sequencing of clones discriminated by RFLP for the study of microbial populations: application to an UASB reactor sample,” Applied Microbiology and Biotechnology, vol. 85, no. 3, pp. 801–806, 2010.
- H. Nadais, M. Barbosa, I. Capela et al., “Enhancing wastewater degradation and biogas production by intermittent operation of UASB reactors,” Energy, vol. 36, no. 4, pp. e2164–e2168, 2011.
- K. M. Singh, A. K. Tripathi, P. R. Pandya et al., “Methanogen diversity in the rumen of Indian Surti buffalo (Bubalus bubalis), assessed by 16S rDNA analysis,” Research in Veterinary Science, vol. 92, no. 3, pp. 451–455, 2012.
- P. P. Chaudhary and S. K. Sirohi, “Dominance of Methanomicrobium phylotype in methanogen population present in Murrah buffaloes (Bubalus bubalis),” Letters in Applied Microbiology, vol. 49, no. 2, pp. 274–277, 2009.
- P. P. Chaudhary, S. K. Sirohi, and J. Saxena, “Diversity analysis of methanogens in rumen of Bubalus bubalis by 16S riboprinting and sequence analysis,” Gene, vol. 493, no. 1, pp. 13–17, 2012.
- E. C. Shin, B. R. Choi, W. J. Lim et al., “Phylogenetic analysis of archaea in three fractions of cow rumen based on the 16S rDNA sequence,” Anaerobe, vol. 10, no. 6, pp. 313–319, 2004.
- K. Tajima, T. Nagamine, H. Matsui, M. Nakamura, and R. I. Aminov, “Phylogenetic analysis of archaeal 16S rRNA libraries from the rumen suggests the existence of a novel group of archaea not associated with known methanogens,” FEMS Microbiology Letters, vol. 200, no. 1, pp. 67–72, 2001.
- N. Tatsuoka, N. Mohammed, M. Mitsumori, K. Hara, M. Kurihara, and H. Itabashi, “Phylogenetic analysis of methyl coenzyme-M reductase detected from the bovine rumen,” Letters in Applied Microbiology, vol. 39, no. 3, pp. 257–260, 2004.
- S. E. Denman, N. W. Tomkins, and C. S. McSweeney, “Quantitation and diversity analysis of ruminal methanogenic populations in response to the antimethanogenic compound bromochloromethane,” FEMS Microbiology Ecology, vol. 62, no. 3, pp. 313–322, 2007.