Potential functional gene diversity involved in methanogenesis and methanogenic community structure in Indian buffalo (Bubalus bubalis) rumen.

Understanding the methanogen community structure and methanogenesis from Bubalus bubalis in India may be beneficial to methane mitigation. Our current understanding of the microbial processes leading to methane production is incomplete, and further advancement in the knowledge of methanogenesis pathways would provide means to manipulate its emission in the future. In the present study, we evaluated the methanogenic community structure in the rumen as well as their potential genes involved in methanogenesis. The taxonomic and metabolic profiles of methanogens were assessed by shotgun sequencing of rumen metagenome by Ion Torrent semiconductor sequencing. The buffalo rumen contained representative genera of all the families of methanogens. Members of Methanobacteriaceae were found to be dominant, followed by Methanosarcinaceae, Methanococcaceae, Methanocorpusculaceae, and Thermococcaceae. A total of 60 methanogenic genera were detected in buffalo rumen. Methanogens related to the genera Methanobrevibacter, Methanosarcina, Methanococcus, Methanocorpusculum, Methanothermobacter, and Methanosphaera were predominant, representing >70 % of total archaeal sequences. The metagenomic dataset indicated the presence of genes involved in the methanogenesis and acetogenesis pathways, and the main functional genes were those of key enzymes in the methanogenesis. Sequences related to CoB--CoM heterodisulfide reductase, methyl coenzyme M reductase, f420-dependent methylenetetrahydromethanopterin reductase, and formylmethanofuran dehydrogenase were predominant in rumen. In addition, methenyltetrahydrofolate cyclohydrolase, methylenetetrahydrofolate dehydrogenase, 5,10-methylenetetrahydrofolate reductase, and acetyl-coenzyme A synthetase were also recovered.