Amrutlal K. Patel

Comparative evaluation of rumen metagenome community using qPCR and MG-RAST

Microbial profiling of metagenome communities have been studied extensively using MG-RAST and other related metagenome annotation databases. Although, database based taxonomic profiling provides snapshots of the metagenome architecture, their reliability needs to be validated through more accurate methods. Here, we performed qPCR based absolute quantitation of selected rumen microbes in the liquid and solid fraction of the rumen fluid of river buffalo adapted to varying proportion of concentrate to green or dry roughages and compared with the MG-RAST based annotation of the metagenomes sequences of 16S r-DNA amplicons and high throughput shotgun sequencing. Animals were adapted to roughage-to-concentrate ratio in the proportion of 50:50, 75:25 and 100:00, respectively for six weeks. At the end of each treatment, rumen fluid was collected at 3xa0h post feeding. qPCR revealed that the relative abundance of Prevotella bryantii was higher, followed by the two cellulolytic bacteria Fibrobacter succinogens and Ruminococcus flavefaciens that accounted up to 1.33% and 0.78% of the total rumen bacteria, respectively. While, Selenomonas ruminantium and archaea Methanomicrobiales were lower in microbial population in the rumen of buffalo. There was no statistically significant difference between the enumerations shown by qPCR and analysis of the shotgun sequencing data by MG-RAST except for Prevotella. These results indicate the variations in abundance of different microbial species in buffalo rumen under varied feeding regimes as well as in different fractions of rumen liquor, i.e. solid and the liquid. The results also present the reliability of shotgun sequencing to describe metagenome and analysis/annotation by MG-RAST.

Characterization of the rumen microbiome of Indian Kankrej cattle (Bos indicus) adapted to different forage diet.

Present study described rumen microbiome of Indian cattle (Kankrej breed) to better understand the microbial diversity and largely unknown functional capacity of the rumen microbiome under different dietary treatments. Kankrej cattle were gradually adapted to a high-forage diet (four animals with dry forage and four with green forage) containing 50xa0% (K1), 75xa0% (K2) to 100xa0% (K3) forage, and remaining concentrate diet, each for 6xa0weeks followed by analysis of rumen fiber adherent and fiber-free metagenomic community by shotgun sequencing using ion torrent PGM platform and EBI-metagenomics annotation pipeline. Taxonomic analysis indicated that rumen microbiome was dominated by Bacteroidetes followed by Firmicutes, Fibrobacter, Proteobacteria, and Tenericutes. Functional analysis based on gene ontology classified all reads in total 157 categories based on their functional role in biological, molecular, and cellular component with abundance of genes associated with hydrolase activity, membrane, transport, transferase, and different metabolism (such as carbohydrate and protein). Statistical analysis using STAMP revealed significant differences (Pu2009<u20090.05) between solid and liquid fraction of rumen (in 65 categories), between all three treatments (in 56 categories), and between green and dry roughage (17 categories). Diet treatment also exerted significant difference in environmental gene tags (EGTs) involved in metabolic pathways for production of volatile fatty acids. EGTs for butyrate production were abundant in K2, whereas EGTs for propionate production was abundant during K1. Principal component analysis also demonstrated that diet proportion, fraction of rumen, and type of forage affected rumen microbiome at taxonomic as well as functional level.

Microbial profiles of liquid and solid fraction associated biomaterial in buffalo rumen fed green and dry roughage diets by tagged 16S rRNA gene pyrosequencing.

AbstractThe microbiome of buffalo rumen plays an important role in animal health and productivity. The rumen bacterial composition of both liquid and solid fraction was surveyed using pyrosequencing of the 16S rRNA gene. Sequences were analyzed using taxonomy-dependent clustering methods and revealed that the dominant ruminal bacteria shared by samples belonged to phyla Bacteroidetes, Firmicutes, Fibrobacteres and Proteobacteria. The core rumen microbiome of the rumen consisted of 10 phyla, 19 classes, 22 orders and 25 families. However, the relative abundance of these bacterial groups was markedly affected by diet composition as well as in type of biomaterial. In animals fed with a green and dry roughage diet, the cellulolytic bacteria, Ruminococcaceae, and Fibrobacteraceae was found in highest abundance in all biomaterials which reflected the need for enhanced fiber-digesting capacity in buffalo. The polysaccharide-degrading Prevotellaceae bacteria were most abundant in buffalo rumen. In taxonomic comparison of rumen bacteria, about 26 genera were differentially abundant among liquid and solid fraction of ruminal fluid. These results highlight the buffalo ruminal microbiome’s ability to adapt to feed with different composition.n

High Potential Source for Biomass Degradation Enzyme Discovery and Environmental Aspects Revealed through Metagenomics of Indian Buffalo Rumen

The complex microbiomes of the rumen functions as an effective system for plant cell wall degradation, and biomass utilization provide genetic resource for degrading microbial enzymes that could be used in the production of biofuel. Therefore the buffalo rumen microbiota was surveyed using shot gun sequencing. This metagenomic sequencing generated 3.9u2009GB of sequences and data were assembled into 137270 contiguous sequences (contigs). We identified potential 2614 contigs encoding biomass degrading enzymes including glycoside hydrolases (GH: 1943 contigs), carbohydrate binding module (CBM: 23 contigs), glycosyl transferase (GT: 373 contigs), carbohydrate esterases (CE: 259 contigs), and polysaccharide lyases (PE: 16 contigs). The hierarchical clustering of buffalo metagenomes demonstrated the similarities and dissimilarity in microbial community structures and functional capacity. This demonstrates that buffalo rumen microbiome was considerably enriched in functional genes involved in polysaccharide degradation with great prospects to obtain new molecules that may be applied in the biofuel industry.

Short hairpin RNA-induced myostatin gene silencing in caprine myoblast cells in vitro.

Myostatin (MSTN) belongs to the transforming growth factor (TGF)-β superfamily and is a potent negative regulator of skeletal muscle development and growth. Dysfunction of MSTN gene either by natural mutation or induced through genetic manipulation (knockout or knockdown) has been reported to increase the remarkable muscle mass in mammalian species. RNA interference (RNAi) is the most promising method for inhibition of gene expression that can be utilized for MSTN gene knockdown by developing short hairpin RNA (shRNA) construct against it. We utilized three antisense RNA expressing vectors with six constructs to knockdown MSTN gene in in vitro caprine myoblast cell culture system. We observed that all six shRNA constructs were successful in MSTN silencing with efficiency ranging from 7 to 46xa0% by quantitative real-time PCR and up to 19xa0% by western blotting. The significant upregulation of interferon response gene OAS1 (5- to 11-fold) in cells transfected with shRNA constructs were indicative of induction of interferon response. This RNAi-based method of increasing muscle mass could provide an alternative strategy to gene knockout methods for improving the production traits and economic properties of livestock.

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 >70xa0% 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.

Insights into resistome and stress responses genes in Bubalus bubalis rumen through metagenomic analysis

Buffalo rumen microbiota experience variety of diets and represents a huge reservoir of mobilome, resistome and stress responses. However, knowledge of metagenomic responses to such conditions is still rudimentary. We analyzed the metagenomes of buffalo rumen in the liquid and solid phase of the rumen biomaterial from river buffalo adapted to varying proportion of concentrate to green or dry roughages, using high-throughput sequencing to know the occurrence of antibiotics resistance genes, genetic exchange between bacterial population and environmental reservoirs. A total of 3914.94xa0MB data were generated from all three treatments group. The data were analysed with Metagenome rapid annotation system tools. At phyla level, Bacteroidetes were dominant in all the treatments followed by Firmicutes. Genes coding for functional responses to stress (oxidative stress and heat shock proteins) and resistome genes (resistance to antibiotics and toxic compounds, phages, transposable elements and pathogenicity islands) were prevalent in similar proportion in liquid and solid fraction of rumen metagenomes. The fluoroquinolone resistance, MDR efflux pumps and Methicillin resistance genes were broadly distributed across 11, 9, and 14 bacterial classes, respectively. Bacteria responsible for phages replication and prophages and phage packaging and rlt-like streptococcal phage genes were mostly assigned to phyla Bacteroides, Firmicutes and proteaobacteria. Also, more reads matching the sigma B genes were identified in the buffalo rumen. This study underscores the presence of diverse mechanisms of adaptation to different diet, antibiotics and other stresses in buffalo rumen, reflecting the proportional representation of major bacterial groups.

Myostatin knockdown and its effect on myogenic gene expression program in stably transfected goat myoblasts

Myostatin, a negative regulator of skeletal muscle mass, is a proven candidate to modulate skeletal muscle mass through targeted gene knockdown approach. Here, we report myostatin (MSTN) knockdown in goat myoblasts stably expressing small hairpin RNA (shRNAs) against MSTN gene through lentivirus vector-mediated integration. We observed 72% (pu2009=u20090.003) and 54% (pu2009=u20090.022) downregulation of MSTN expression with sh2 shRNA compared to empty vector control and untransduced myoblasts, respectively. The knockdown of MSTN expression was accompanied with concomitant downregulation of myogenic regulatory factor MYOD (77%, pu2009=u20090.001), MYOG (94%, pu2009=u20090.000), and MYF5 (36%, pu2009=u20090.000), cell cycle regulator p21 (62%, pu2009=u20090.000), MSTN receptor ACVR2B (23%, pu2009=u20090.061), MSTN antagonist follistatin (81%, pu2009=u20090.000), and downstream signaling mediators SMAD2 (20%, pu2009=u20090.060) and SMAD3 (49%, pu2009=u20090.006). However, the expression of MYF6 was upregulated by 14% compared to control lentivirus-transduced myoblasts (pu2009=u20090.354) and 79% compared to untransduced myoblasts (pu2009=u20090.018) in sh2 shRNA-transduced goat myoblasts cells. Although, MSTN knockdown led to sustained cell proliferation of myoblasts, the myoblasts fusion was suppressed in both MSTN knocked down and control lentivirus-transduced myoblasts. The expression of interferon response gene OAS1 was significantly upregulated in control lentivirus (10.86-fold; pu2009=u20090.000)- and sh2 (1.71-fold; pu2009=u20090.002)-integrated myoblasts compared to untransduced myoblasts. Our study demonstrates stable knockdown of MSTN in goat myoblasts cells and its potential for use in generation of transgenic goat by somatic cell nuclear transfer.

Goat activin receptor type IIB knockdown by artificial microRNAs in vitro.

Activin receptor type IIB (ACVR2B) has been known to negatively regulate the muscle growth through mediating the action of transforming growth factor beta superfamily ligands. Recently, the artificial microRNAs (amiRNAs) which are processed by endogenous miRNA processing machinery have been proposed as promising approach for efficient gene knockdown. We evaluated amiRNAs targeting goat ACVR2B in HEK293T and goat myoblasts cells. The amiRNAs were designed based on the miR-155 backbone and cloned in 5′- and 3′-UTR of GFP reporter gene under the CMV promoter. Although both 5′- and 3′-UTR-amiRNAs vectors showed efficient synthesis of GFP transcripts, amiRNAs in 5′-UTR drastically affected GFP protein synthesis in transfected goat myoblast cells. Among the four amiRNAs targeting ACVR2B derived from either 5′- or 3′-UTR, ami318 showed highest silencing efficiency against exogenously co-expressed ACVR2B in both 293T and goat myoblast cells whereas ami204 showed highest silencing efficiency against endogenous ACVR2B in goat myoblasts cells. The 3′-UTR-derived amiRNA exerted higher knockdown efficiency against endogenous ACVR2B at transcript level whereas 5′-UTR-derived amiRNAs exerted higher knockdown efficiency at protein level. The expression of ACVR2B showed positive correlation with the expression of MYOD (ru2009=u20090.744; pu2009=u20090.009) and MYOG (ru2009=u20090.959; pu2009=u20090.000) in the amiRNA-transfected myoblasts. Although both 5′- and 3′-UTR–amiRNA vectors led to substantial induction of interferon response, the magnitude of the response was found to be higher with the 3′-UTR–amiRNA vectors.

Buffalo alpha S1-casein gene 5′-flanking region and its interspecies comparison

The expression of milk protein genes is tightly regulated in a spatio-temporal manner through the combinatorial interaction of lactogenic hormones and a set of transcription factors mediating developmental and tissue-specific gene expression. The recruitment of a unique set of transcription factors is determined by the cis-regulatory motifs present in the gene promoter region. Here, we report the isolation, sequencing, structural analysis and interspecies comparison of the 5′cis-regulatory region of the buffalo alpha S1 (αS1)-casein gene. The proximal promoter region of the buffalo αS1-casein gene harbored the insertion of a 72-bp fragment of long interspersed nuclear element of the L1_BT retrotransposon family. Among the core and vertebrate-specific promoter elements, the motifs for the binding of Brn POU domain factors (BRNF), Lim homeodomain factors (LHXF), NK6 homeobox transcription factors (NKX6), nuclear factor kappa B/c-rel (NFKB), AT-rich interactive domain factor (ARID), Brn POU domain factor 5 (BRN5), pancreatic and intestinal homeodomain transcription factor (PDX1), Distal-less homeodomain transcription factors (DLXF), T-cell factor/lymphoid enhancer-binding factor-1 (LEFF) and GHF-1 pituitary-specific POU domain transcription factor (PIT1) were over-represented in the αS1-casein gene regulatory region (Z score >4.0). The Multiple EM for Motif elicitation predicted three motifs which consisted of the sequences known to bind mammary gland factor/signal transducer and activator of transcription 5 (MGF/STAT5), estrogen receptor-related alpha (ERα), steroidogenic factor 1 (SF1) and glucocorticoid receptor (GR), indicating their potential role in the mammary gland-specific gene expression. The interspecies comparison of the proximal promoter region revealed conserved sequences for TATA boxes and MGF/STAT5 in all species, whereas activator protein 1 (AP1), pregnancy-specific mammary nuclear factor (PMF), CCAAT/enhancer binding protein (C/EBP), double-stranded and single-stranded DNA-binding protein 1 (DS1 and SS), ying and yang factor 1 (YY1), and GR half-sites were among ruminants. The functional significance of the L1_BT retrotransposon insertion on the buffalo αS1-casein gene expression needs to be experimentally validated.