Genetic diversity analysis of an indigenous Chinese buffalo breed and hybrids based on microsatellite data

Chinese native buffaloes have faced the threat of extinction, along with an increase in crossbreeding with domesticated river buffaloes; consequently, conservation of local buffalo genetic resources has become a priority. A Chinese native breed, Jianghan, is often crossed intentionally and unintentionally with imported breeds from India and Pakistan, Murrah, and Nili-Ravi. A total of 128 buffaloes of the breeds Jianghan, Murrah, and Nili-Ravi and their presumed hybrid offspring were genotyped for 10 microsatellite markers. Heterozygosity and Wright's F-statistics were calculated to determine the genetic variation in those populations. The observed average heterozygosities ranged from 0.836 (Murrah) to 0.986 (Jianghan), higher than the expected heterozygosities and all the inbreeding values within the populations were negative. The genetic distances between the presumed hybrid buffaloes and the two imported river type dairy buffalo breeds (Murrah and Nili-Ravi) were lower than with the native Jianghan, indicating strong contributions of the imported breeds to this presumed hybrid buffalo population. This information will be useful for the development of rational breeding for the dairy buffalo industry and for conservation strategies for the Jianghan buffalo.     

Identification of X- and Y-chromosome bearing buffalo (Bubalus bubalis) sperm

The aim of this study was to identify the difference in DNA content characterizing the X- and Y-chromosome bearing sperm of buffalo. Sperm from six Murrah buffaloes and six Nili-Ravi buffaloes were collected and stained with Hoechst 33342 followed by flow cytometry analysis of the DNA content. Two symmetrical, separate but overlapping peaks presumed to be X- and Y-chromosome bearing sperm were detected. The difference in fluorescence intensity, which related to the DNA content, between the X- and Y-sperm was 3.59+/-0.11% for Murrah buffalo and 3.55+/-0.14% for Nili-Ravi buffalo, respectively. Significant differences were observed among males within each breed, but there were no differences between the averages of the two breeds. The results indicate that flow cytometric sorting of X- and Y-sperm of buffalo is feasible.     

Biometry of frozen-thawed sperm from eight breeds of Indian buffaloes (Bubalus bubalis)

Sperm morphometry, in combination with other objective traits, can be useful for developing a fertility index. The objective of the present study was to measure various biometric end points of frozen-thawed sperm from eight breeds of Indian buffaloes (Murrah, Surti, Tarai, Mehsana, Jaffrabadi, Bhadawari, Pandharpuri and Nili-Ravi). The sperm head of Pandharpuri buffaloes had the greatest length (10.21 microm), width (6.05 microm), area (52.31 microm(2)) and perimeter (31.86 microm). The ratio of sperm width to length was also greatest (0.61) in Pandharpuri as well as in two other breeds, viz. Nili-Ravi and Jaffrabadi. Murrah had the smallest sperm head width (4.75 microm), area (41.65 microm(2)) and perimeter (29.17 microm), but its sperm tail was longest (57.02 microm), along with that of Jaffrabadi buffaloes (56.96 microm). Based on mean values of sperm tail length, mid piece length and its width the eight buffalo breeds were categorized into three, four and five groups, respectively. Multivariate analysis and clustering put six breeds (Surti, Tarai, Mehsana, Jaffrabadi, Bhadawari and Nili-Ravi) in one cluster, whereas Murrah and Pandharpuri appeared as separate entities.     

Genetic variation and relationships among eight Indian riverine buffalo breeds

Twenty-seven microsatellite loci were used to define genetic variation and relationships among eight Indian riverine buffalo breeds. The total number of alleles ranged from 166 in the Toda breed to 194 each in the Mehsana and the Murrah. Significant departures from the Hardy-Weinberg equilibrium were observed for 26 locus-breed combinations due to heterozygote deficiency. Breed differentiation was analysed by estimation of F(ST) index (values ranging from 0.75% to 6.00%) for various breed combinations. The neighbour-joining tree constructed from chord distances, multidimensional scaling (MDS) display of F(ST) values and Bayesian clustering approach consistently identified the Toda, Jaffarabadi, and Pandharpuri breeds as one lineage each, and the Bhadawari, Nagpuri, Surati, Mehsana and Murrah breeds as admixture. Analysis of molecular variance refuted the earlier classification of these breeds proposed on the basis of morphological and geographical parameters. The Toda buffaloes, reared by a tribe of the same name, represent an endangered breed from the Nilgiri hills in South India. Divergence time of the Toda buffaloes from the other main breeds, calculated from Nei's standard genetic distances based on genotyping data on seven breeds and 20 microsatellite loci, suggested separation of this breed approximately 1800-2700 years ago. The results of the present study will be useful for development of rational breeding and conservation strategies for Indian buffaloes.     

Analysis of culling probability in dairy buffalo using survival models

In order to contribute to the genetic breeding programs of buffaloes, this study aimed to determine the influence of environmental effects on the stayability (ST) of dairy female Murrah buffalo in the herd. Data from 1016 buffaloes were used. ST was defined as the ability of the female to remain in the herd for 1, 2, 3, 4, 5 or 6 years after the first calving. Environmental effects were studied by survival analysis, adjusted to the fixed effects of farm, year and season of birth, class of first-lactation milk yield and age at first calving. The data were analyzed using the LIFEREG procedure of the SAS program that fits parametric models to failure time data (culling or ST = 0), and estimates parameters by maximum likelihood estimation. Breeding farm, year of birth and first-lactation milk yield significantly influenced (P < 0.0001) the ST to the specific ages (1 to 6 years after the first calving). Buffaloes that were older at first calving presented higher probabilities of being culled 1 year after the first calving, without any effect on culling at older ages. Buffaloes with a higher milk yield at first calving presented a lower culling probability and remained for a longer period of time in the herd. The effects of breeding farm, year of birth and first-lactation milk yield should be included in models used for the analysis of ST in buffaloes.     

The role of MC1R gene in buffalo coat color

Melanocortin-1 receptor (MC1R) plays a major role in pigmentation in many species. To investigate if the MC1R gene is associated with coat color in water buffalo, the coding region of MC1R gene of 216 buffalo samples was sequenced, which included 49 black river buffalo (Murrah and Nili-Ravi), 136 swamp buffalo (Dehong, Diandongnan, Dechang, Guizhou, and Xilin) with white and gray body, and 31 hybrid offspring of river buffalo Nili-Ravi (or Murrah) and swamp buffalo. Among the three variation sites found, SNP684 was synonymous, while SNP310 and SNP384 were nonsynonymous, leading to p.S104G and p.I128M changes, respectively. Only Individuals carrying homozygote EBR/EBR were black. The genotype and phenotype analysis of the hybrid offspring of black river buffalo and gray swamp buffalo further revealed that the river buffalo type allele EBR or the allele carrying the amino acid p.104S was important for the full function of MC1R. The in silico functional analysis showed that the amino acid substitutions p.G104S and p.M128I had significant impact on the function of MC1R. Above results indicate that the allele EBR or the allele carrying the amino acid p.104S was associated with the black coat color in buffalo.     

Evolution and diversity studies of innate immune genes in Indian buffalo (<Emphasis Type="Italic">Bubalus bubalis</Emphasis>) breeds using next generation sequencing

Indian buffalo (Bubalus bubalis) breeds are acknowledged for their disease resistance power, hardiness and considerable milk production. River buffaloes are more disease resistance despite of tropical environment of India than cattle breeds of the world. The genetic polymorphisms of innate immune genes play a prominent role in disease susceptibility of livestock. Here in our study, we mainly focused on SNPs in ten selected candidate innate immune genes viz CHGA, CHGB, CHGC, Slc11a1, Slc11a2, DEFB1, BNBD4, BNBD5, LAP and TAP among three most economically important Indian buffalo breeds viz., Jaffrabadi, Mehsani and Murrah using amplicon sequencing on Ion torrent PGM sequencing platform. Polymorphism identification using by GATK tool of Unified Genotyper revealed 274, 245 and 224 high quality SNPs in Jaffrabadi, Mehsani and Murrah buffalo breeds, respectively as well as 168 common SNPs among these three breeds and uniquely detected SNPs were 31, 15 and 37 SNPs in Jaffrabadi, Mehsani and Murrah breeds respectively. Consequences of non-synonymous SNPs were analysed using several computational tools viz., I-Mutant, SIFT, PolyPhen and PROVEAN which identified three deleterious nsSNPs. This wealth of sequence information will be productive for conservation studies, selective breeding and designing future strategies for identifying disease associations.     

Effect of oestrus synchronization methods on oestrus behaviour, timing of ovulation and pregnancy rate during the breeding and low breeding seasons in Nili-Ravi buffaloes

The objective of this study was to determine the effect of oestrous synchronization methods on oestrous behaviour, timing of ovulation and pregnancy rate during the breeding and low breeding seasons in Nili-Ravi buffaloes. In Experiment 1, oestrous behaviour and timing of ovulation were determined from (n=34) oestruses. The mean (+/- S.E.M.) time of ovulation after the onset of standing oestrus was greater (P<0.05) in PGF(2alpha)-induced luteolysis (30.6+/-1.5h) compared to Ovsynch buffaloes (15.0+/-0.8h). In Experiment 2, pregnancy rates were compared between two methods of synchronization (detected oestrus and Ovsynch protocol) during the breeding and low breeding seasons. Pregnancy rates of buffaloes bred at detected oestrus (62.5%) or by the Ovsynch protocol (36.3%) during the breeding season did not differ significantly (P>0.05) from those which were inseminated during the low breeding season (55.5%) and (30.4%), respectively. This study demonstrates clearly that (1) timing of ovulation in Nili-Ravi buffalo is about 30h after the onset of standing oestrus and (2) buffaloes can be successfully synchronized with optimum fertility using either PGF(2alpha) alone (detected oestrus) or using (Ovsynch protocol) during low breeding season, to calve during the period when milk availability is short.     

Mitochondrial sequence-based evolutionary analysis of riverine–swamp hybrid buffaloes of India indicates novel maternal differentiation and domestication patterns

In this study, mitochondrial D-loop sequence data on riverine, swamp and hybrid buffaloes from India have been generated and compared with other reported Indian riverine, Chinese and Bangladeshi swamp buffalo populations. Sequence analysis revealed the presence of 132 haplotypes, with a haplotype diversity of 0.9611 ± 0.0045 and a nucleotide diversity of 0.04801 ± 0.00126. For the first time, the existence of riverine–swamp hybrids among the Indian Chilika buffalo population has been recorded, having 49 chromosomes, which was also confirmed by mitochondrial haplotype sharing between Chilika and Indian swamp as well as Chinese swamp buffalo populations in the network analysis. Phylogenetic analysis documents the sharing of reported pre-domestication haplogroups ‘SA1’, ‘SA2’, ‘SA3’ and ‘SB1’ between the Chilika and swamp buffalo populations of India, China and Bangladesh, an indication of the migration of swamp buffaloes towards Bangladesh and adjoining lower parts of India and north towards Chinese domestication sites. The results have also been supplemented by multidimension scaling, grouping Indian and Chinese swamp buffaloes more closely together with Bangladeshi buffaloes, but into a separate quadrant, whereas Chilika grouped away from other riverine as well as swamp buffaloes. These findings thus confirm the previous reports that the northeast region of India, close to the Indo-China border, is the point of evolution of swamp buffaloes with multiple sites of domestication.     

Investigation of transferability of BovineSNP50 BeadChip from cattle to water buffalo for genome wide association study

Cattle and water buffalo belong to the same subfamily Bovinae and share chromosome banding and gene order homology. In this study, we used genome-wide Illumina BovineSNP50 BeadChip to analyze 91 DNA samples from three breeds of water buffalo (Nili-Ravi, Murrah and their crossbred with local GuangXi buffalos in China), to demonstrate the genetic divergence between cattle and water buffalo through a large single nucleotide polymorphism (SNP) transferability study at the whole genome level, and performed association analysis of functional traits in water buffalo as well. A total of 40,766 (75.5 %) bovine SNPs were found in the water buffalo genome, but 49,936 (92.5 %) were with only one allele, and finally 935 were identified to be polymorphic and useful for association analysis in water buffalo. Therefore, the genome sequences of water buffalo and cattle shared a high level of homology but the polymorphic status of the bovine SNPs varied between these two species. The different patterns of mutations between species may associate with their phenotypic divergence due to genome evolution. Among 935 bovine SNPs, we identified a total of 9 and 7 SNPs significantly associated to fertility and milk production traits in water buffalo, respectively. However, more works in larger sample size are needed in future to verify these candidate SNPs for water buffalo.     

8个亚洲水牛群体的遗传结构分析

应用13个微卫星标记结合荧光–多重PCR技术, 对德昌水牛、兴隆水牛、富钟水牛、温州水牛、东流水牛、福安水牛和两个引进品种摩拉水牛、尼里－拉菲水牛进行遗传结构分析。结果表明: 8个水牛群体在13个微卫星座位中共检测到157个等位基因, 其中7个群体具有各自的特有等位基因, 其和为23; 8个群体的有效等位基因数在2.2908～4.2308之间, 杂合度在0.4951～0.7194之间, 多态信息含量在0.4495～0.6776之间; 有11个座位为高度多态座位, 是适合分析水牛遗传多样性的多态标记; 聚类分析表明富钟水牛和东流水牛先聚在一起, 再与兴隆水牛聚在一起, 然后与温州水牛和福安水牛聚在一起, 德昌水牛独自聚为一类; 两个引进品种聚为一类。     

Analysis of Rumen Methanogen Diversity in Water Buffaloes (Bubalus bubalis) Under Three Different Diets

The water buffalo (Bubalus bubalis) is a prominent livestock species for the production of milk and meat in many countries. We investigated the diversity of rumen methanogens in Mediterranean water buffaloes maintained in Brazil under different diets: corn silage, grazing pasture, or sugar cane. A total of 467 clones were isolated from three methanogen 16S rRNA gene clone libraries that each represented a distinct feed type. The 467 clones were assigned to 19 species-level operational taxonomic units (OTUs). Four OTUs were represented in all three libraries, eight OTUs were library-specific, six OTUs were found in only the corn silage and pasture grazing libraries, and one OTU was shared only between pasture grazing and sugar cane libraries. We found that Methanobrevibacter-related sequences were the most abundant in the water buffaloes sampled for our analysis, in contrast to previously reported studies showing that Methanomicrobium mobile-like methanogens were the most abundant methanogens in water buffaloes of Murrah and Surti breeds sampled in India. Considering the worldwide distribution of water buffaloes and the likely wide variety of diets provided, our results combined with studies from other groups support that larger scope analyses of microbiomes for this livestock species would provide great insight into the contribution of geographical location, breed, and diet in determining the population structure of rumen microorganisms.     

Evaluating the potential for reintroducing the endangered wild water buffalo (Bubalus arnee) in Kanha National Park,central India

Reintroducing megafauna to their historic range is an effective strategy to halt their extinctions and restore ecosystems. Wild water buffalo (Bubalus arnee) is an endangered megaherbivore that is lost from 95% of its range. About 90% of its global population (less than 4000) resides within India, in two isolated populations: northeast and central India. The central Indian population is on the verge of extinction and warrants urgent conservation interventions. We assess the potential and provide a strategy for reintroducing buffaloes in Kanha National Park, India. Habitat suitability using the global occurrence of buffalo revealed low-lying grasslands with least human pressure found in Kanha (390 km²) to be suitable. Within this suitable range, we evaluated vegetation composition, forage biomass, and potential carrying capacity. Multidimensional ordination classified these suitable sites into moist and dry grassland clusters. Moist grasslands were found to better suit reintroduction due to lower grazing pressure, higher productivity, and availability of perennial water sources. Distinct matriarchal clades of swamp and river buffalo were observed. Within the river buffalo clade, Indian wild buffaloes formed a distinct cluster with close proximity between northeast and central Indian samples, suggesting northeast buffaloes could be sourced for a founding population in Kanha. Following IUCN guidelines, we discuss the reintroduction strategy that could sustain approximately 200 buffaloes in the area and subsequently expand to other suitable habitats in central India. If implemented as proposed, recovery of this lost ecosystem engineer will help to restore grasslands, and swamp habitats, as well as contribute to its global conservation efforts.     

Analysis of mitochondrial D-loop region casts new light on domestic water buffalo (Bubalus bubalis) phylogeny

The phylogeny of water buffaloes (Bubalus bubalis) is still a matter of discussion, especially if the two types of domestic water buffalo (swamp and river) derived from different domestication events or if they are products of human selection. To obtain more insight, we analyzed the entire mitochondrial D-loop region of 80 water buffaloes of four different breeds, i.e., 19 swamp buffaloes (Carabao) and 61 river buffaloes (Murrah, Jafarabadi, and Mediterranean), sampled in Brazil and Italy. We detected 36 mitochondrial haplotypes with 128 polymorphic sites. Pooled with published data of South-East Asian and Australian water buffaloes and based on comprehensive median-joining network and population demography analyses we show evidence that both river and swamp buffaloes decent from one domestication event, probably in the Indian subcontinent. However, the today swamp buffaloes have an unravelled mitochondrial history, which can be explained by introgression of wild water buffalo mtDNA into domestic stocks. We are also discussing indications for an independent domestication of buffaloes in China.     

Dominance of Methanomicrobium phylotype in methanogen population present in Murrah buffaloes (Bubalus bubalis)

Aims:  To study the diversity of rumen methanogens in Murrah buffaloes ( Bubalus bubalis ) from North India by using 16S rRNA gene libraries obtained from the pooled rumen content from four animals and using suitable software analysis.
Methods and Results:  Genomic DNA was isolated and PCR was set up by using specific primers. Amplified product was cloned into a suitable vector and the positive clones were selected on the basis of blue–white screening and sequenced. The resulting nucleotide sequences were arranged in the phylogenetic tree. A total of 108 clones were examined, revealing 17 different 16S rRNA gene sequences or phylotypes. Of the 17 phylotypes, 15 (102 of 108 clones) belonged to the genus Methanomicrobium , indicating that the genus Methanomicrobium is the most dominant component of methanogen populations in Murrah buffaloes ( Bubalus bubalis ) from North India. The largest group of clones (102 clones) was more than 98% similar to Methanomicrobium mobile . BLAST analysis of the rumen contents from individual animals also revealed 17 different phylotypes with a range of 3–10 phylotypes per animal.
Conclusion:  Methanomicrobium phylotype is the most dominant phylotype of methanogens present in Murrah buffaloes ( Bubalus bubalis ).
Significance and Impact of the Study:  Effective strategies can be made to inhibit the growth of Methanomicrobium phylotype to reduce the methane emission from rumen contents and thus help in preventing global warming.     

Serum Electrolyte and Mineral Variations During Pregnancy and Lactation in Nili-Ravi Buffalo

The aim of present study was to determine the changes of sodium, potassium, chloride, copper, and zinc in serum of Nili-Ravi buffalo during pregnancy and lactation. The study was carried out on 25 Nili-Ravi buffaloes during March 2008 to February 2009 at Buffalo Research Institute, Pattoki, District Kasur. Blood samples were taken from each buffalo during early pregnancy (Phase EP), i.e., between 1 and 3 months gestation, mid pregnancy (Phase MP), i.e., between 4 and 6 months gestation, late pregnancy (Phase LP), i.e., between 7 and 10 months gestation, and during lactation (Phase LT). Serum sodium, potassium, and chloride were determined on a clinical chemistry analyzer whereas copper and zinc were determined on an atomic absorption spectrophotometer. The mean serum sodium, potassium, and chloride concentrations varied non-significantly (P > 0.05) during pregnancy and lactation. With advancing pregnancy, copper increased and was significantly high (P < 0.05) during late pregnancy whereas during lactation copper concentrations dropped significantly (P < 0.05). The zinc concentrations increased non-significantly during early and mid pregnancy. There were significant differences (P < 0.05) in zinc concentrations during lactation compared with late pregnancy. This study demonstrates that concentrations of serum sodium, potassium, and chloride remain unchanged during pregnancy and lactation whereas serum copper concentrations increased and zinc concentrations decreased with advancing pregnancy because of the demand of the fetus at different periods of pregnancy.     

Molecular phylogeny of Pakistani riverine buffalo based on genetic variability of mitochondrial cytochrome b gene

Mitochondrial cytochrome b gene is considered to be one of the best markers for breed characterization as well as studying the ancestry in the vertebrates due to its exclusive maternal inheritance. DNA fingerprinting by single nucleotide polymorphism is most reliable and widely used molecular technique in modern forensics and is being considered in this study. Partial sequencing of 1,061?bp of aforementioned gene from 14580 to 15643 was conducted in two famous Pakistani buffalo breeds named Nili-Ravi and Kundi. In which we explore seven haplotypes within earlier and none in the latter breed. Nili-Ravi is polymorphic at four codons of this gene, and the protein translation is also different from the reference sample while monomorphic at three codons with no amino acid replacement. Haplotypes frequency distribution of these four haplotypes named NR3, NR4, NR5, NR7 revealed that the prevalence of each haplotype is 0.04?% in the Pakistani buffalo population of this Nili-Ravi breed while complete homoplasmy was observed in the Kundi breed population. Nili-Ravi breed of buffalo is genetically more variable than the Kundi breed as far as the gene in subject is concerned. It means later breed has spent more time to propagate its wild type haplotype which make this breed more ancestral as compare to Nili-Ravi. Secondly both breeds share their common ancestors with regional water buffalo rather than the swamp one.     

Seroprevalence of Neospora caninum antibodies in cattle and water buffaloes in India

Neospora caninum is now recognized as a major cause of abortion in cattle worldwide, but there is no report of N. caninum infection in cattle in India. Serum samples from 427 dairy cattle and 32 dairy water buffaloes from 7 organized dairy farms located in Punjab, India, were tested for N. caninum antibodies using a commercial monoclonal antibody-based competitive enzyme-linked immunosorbent assay (ELISA). Antibodies to N. caninum were found in 35 of 427 cattle from 6 of the 7 farms; 9.6% of cows, 5.1% of heifers, and 5.0% of calves were seropositive, suggesting postnatal transmission of N. caninum on the farm. Antibodies to N. caninum were found in 16 of 32 buffaloes tested from 2 dairy farms. In total, 64 cattle and 16 buffalo sera already tested by ELISA were also evaluated by an indirect fluorescent antibody test (IFAT) to verify ELISA results. Of the 64 cattle samples, 29 sera were negative by both tests and of the 35 ELISA-positive sera, 12 had IFAT titers of 1:100 or higher (1 had IFAT titer of 100, 2 had IFAT titer of 200, and 9 had IFAT titers of 400 or higher). Of the 16 buffalo sera positive by ELISA, 1 had an IFAT titer of 1:400. Thus, antibodies to N. caninum were demonstrated in cattle sera by 2 serologic methods. To our knowledge this is the first report of N. caninum infection in cattle and buffaloes in India.     

Identification of novel single nucleotide polymorphisms in the DGAT1 gene of buffaloes by PCR-SSCP

Diacylglycerol O-acyltransferase 1 (DGAT1) is a microsomal enzyme that catalyzes the final step of triglyceride synthesis. The DGAT1 gene is a strong functional candidate for determining milk fat content in cattle. In this work, we used PCR-SSCP (polymerase chain reaction-single-strand conformation polymorphism) and DNA sequencing to examine polymorphism in the region spanning exon 7 to exon 9 of the DGAT1 gene in Murrah and Pandharpuri buffaloes. Three alleles (A, B and C) and four novel single-nucleotide polymorphisms were identified in the buffalo DGAT1 gene. The frequencies of the alleles differed between the two buffalo breeds, with allele C being present in Murrah but not in Pandharpuri buffalo. The allele variation detected in this work may influence DGAT1 expression and function. The results described here could be useful in examining the association between the DGAT1 gene and milk traits in buffalo.     

Methyl coenzyme M reductase (<Emphasis Type="Italic">mcrA</Emphasis>) gene based phylogenetic analysis of methanogens population in Murrah buffaloes (<Emphasis Type="Italic">Bubalus bubalis</Emphasis>)

The aim of the present study was to decipher the diversity of methanogens in rumen of Murrah buffaloes so that effective strategies can be made in order to mitigate methane emission from these methanogens. In the present study diversity of rumen methanogens in Murrah buffaloes (Bubalus bubalis) from North India was evaluated by using mcr-A gene library obtained from the pooled PCR product from four animals and by using MEGA4 software. A total of 104 clones were examined, revealing 26 different mcr-A gene sequences or phylotypes. Of the 26 phylotypes, 16 (64 of 104 clones) were less than 97% similar to any of the cultured strain of methanogens. Seven clone sequences were clustered with Methanomicrobium mobile and three clone sequences were clustered with Methanobrevibacter gottschalkii during the phylogenetic analysis. Uncultured group of methanogens comes out to be the major component of the methanogens community structure in Murrah buffaloes. Methanomicrobium phylotype comes out to be major phylotype among cultured methanogens followed by Methanobrevibacter phylotype. These results help in making effective strategies to check the growth of dominant communities in the rumen of this animal which in turn help in the reduction of methane emission in the environment and ultimately helps us in fighting with the problem of global warming.