Immunomodulation of IL-1, IL-6 and IL-8 cytokines by Prosopis juliflora alkaloids during bovine sub-clinical mastitis

The current work is focused on establishing therapeutic protocol using unconventional drugs of herbal origin and studying their mechanism of action at molecular level in the treatment of bovine sub-clinical mastitis. It explores the potential of different cytokines which can be used for diagnosis, prognosis and monitoring of bovine sub-clinical mastitis. Prosopis juliflora alkaloids was administered intramammarily in 24 sub-clinically affected quarters once a day for 5 consecutive days at the rate of 10 ml of 1% formulation. In 18 disease control quarters, sterile normal saline was infused. The bacterial cultural examination, somatic cell count (SCC) and cytokines (IL-1β, IL-6, IL-8, IL-12, GM-CSF, IFN-γ, TNF-α) expression by real-time PCR were evaluated on day 7, 14, 21 and 28 post-last treatment from milk samples. Around 75.0% of treatment group quarters showed significant (p < 0.05) reduction in SCC on day 28 post-last treatment, whereas 94.4% control group quarters did not show any significant decline in SCC. 58.3% of treated quarters showed both bacteriological cure as well as significant (p < 0.05) reduction in SCC on day 28 post-last treatment. While, among control group quarters, 83.3% quarters not only remained bacteriological positive, they also did not show any significant decline in SCC. The in vitro antimicrobial activity of alkaloids of P. juliflora was evaluated. Lower concentrations of alkaloids (0.25% and 0.50%) dissolved in normal saline showed zone of inhibition against 12 out of 15 isolates, however higher concentration (1, 1.5, 2, 2.5 and 5%) showed zone of inhibition against all 15 bacterial isolates. The gene expression level of IL-1β, IL-8 and IFN-γ cytokines exhibited significant difference between healthy and sub-clinically affected quarters highlighting the potential of these cytokines in the diagnosis of bovine sub-clinical mastitis. Down-regulation of IL-1, IL-6, IL-8 and IFN-γ cytokines in treated quarters can be explored for making the prognosis and monitoring post-treatment disease progression of bovine sub-clinical mastitis. The P. juliflora alkaloid demonstrated strong in vitro and in vivo antibacterial activity, along with causing immunomodulation by enhancing post-treatment gene expression of IL-1, IL-6 and IL-8 cytokines. Therefore, P. juliflora alkaloids hold a strong claim as an effective alternative herbal therapy in bovine sub-clinical mastitis.

Identification of differentially expressed microRNAs in Sahiwal (Bos indicus) breed of cattle during thermal stress

microRNAs (miRNAs) are a class of small non-coding RNAs that play key roles in post transcriptional gene regulation that influence various fundamental cellular processes, including the cellular responses during environmental stresses. However, perusal of literatures revealed few reports on the differential expression of miRNA during thermal stress in Indian native (Bos indicus) cattle breeds. The present investigation aimed to identify differentially expressed miRNAs during thermal stress in Sahiwal (Bos indicus) dairy cattle breed of India, adapted with tropical climate over a long period of time. Stress responses of the animals were characterized by determining various physiological as well as biochemical parameters and differential expression profile of major heat shock protein genes. Ion Torrent deep sequencing and CLC-genomic analysis identified a set of differentially expressed miRNAs during summer and winter seasons. Most of the identified differentially expressed miRNAs were found to target heat shock responsive genes especially members of heat shock protein (HSP) family. Real-time quantification-based analysis of selected miRNAs revealed that bta-mir-1248, bta-mir-2332, bta-mir-2478, and bta-mir-1839 were significantly (p < 0.01) over expressed while bta-mir-16a, bta-let-7b, bta-mir-142, and bta-mir-425 were significantly (p < 0.01) under expressed during summer in comparison to winter. The present study enlists differentially expressed miRNAs at different environmental temperatures in Sahiwal (Bos indicus) that may be importance for further understanding the role of miRNAs on thermo-regulatory mechanisms.

Illumina Next Generation Sequencing for the Analysis of Eimeria Populations in Commercial Broilers and Indigenous Chickens

Eimeria species parasites can cause the enteric disease coccidiosis, most notably in chickens where the economic and welfare implications are significant. Seven Eimeria species are recognized to infect chickens, although understanding of their regional occurrence, abundance, and population structure remains limited. Reports of Eimeria circulating in chickens across much of the southern hemisphere with cryptic genotypes and the capacity to escape current anticoccidial vaccines have revealed unexpected levels of complexity. Consequently, it is important to supplement validated species-specific molecular diagnostics with new genus-level tools. Here, we report the application of Illumina MiSeq deep sequencing to partial 18S rDNA amplicons generated using Eimeria genus-specific primers from chicken caecal contents collected in India. Commercial Cobb400 broiler and indigenous Kadaknath type chickens were sampled under field conditions after co-rearing (mixed type farms, n = 150 chickens for each) or separate rearing (single type farms, n = 150 each). Comparison of MiSeq results with established Internal Transcribed Spacer (ITS) and Sequence Characterised Amplified Region (SCAR) quantitative PCR assays suggest greater sensitivity for the MiSeq approach. The caecal-dwelling Eimeria tenella and E. necatrix dominated each sample set, although all seven species which infect chickens were detected. Two of the three cryptic Eimeria genotypes were detected including OTU-X and OTU-Y, the most northern report for the latter to date. Low levels of DNA representing other Eimeria species were detected, possibly representing farm-level contamination with non-replicating oocysts or Eimeria DNA, or false positives, indicating a requirement for additional validation. Next generation deep amplicon sequencing offers a valuable resource for future Eimeria studies.

Emergence of a genotype I variant of avian infectious bronchitis virus from Northern part of India

Infectious bronchitis virus (IBV) is one of the foremost causes of a persistent economic burden to poultry industries worldwide. IBV belongs to the genus Gammacoronavirus within the family Coronaviridae. The IBV infection leads to respiratory and nephrogenic symptoms in broiler chickens. In addition, its infection leads to reduced fertility and hatchability in layer birds. We determined the first complete genome sequence of a variant IBV from an outbreak in Haryana state of the Northern part of India using next generation sequencing. On phylogenetic analysis of the IBV isolate, it clustered with genotype I lineage 1 (GI-1). The deduced amino acid sequence of S gene of IBV isolates showed a high level of identity with strains isolated from Tamil Nadu and the reference vaccine strains. Our result suggests that the IBV virus isolated from unvaccinated chicken flocks in North India might be a revertant strain originally evolved from the live attenuated vaccine strains used in the region. Determination of the complete genome sequence of additional IBV isolates from India is necessary to understand the epidemiology of IBV in India.

Microbial diversity and community composition of caecal microbiota in commercial and indigenous Indian chickens determined using 16s rDNA amplicon sequencing

BACKGROUND

The caecal microbiota plays a key role in chicken health and performance, influencing digestion and absorption of nutrients, and contributing to defence against colonisation by invading pathogens. Measures of productivity and resistance to pathogen colonisation are directly influenced by chicken genotype, but host driven variation in microbiome structure is also likely to exert a considerable indirect influence.

METHODS

Here, we define the caecal microbiome of indigenous Indian Aseel and Kadaknath chicken breeds and compare them with the global commercial broiler Cobb400 and Ross 308 lines using 16S rDNA V3-V4 hypervariable amplicon sequencing.

RESULTS

Each caecal microbiome was dominated by the genera Bacteroides, unclassified bacteria, unclassified Clostridiales, Clostridium, Alistipes, Faecalibacterium, Eubacterium and Blautia. Geographic location (a measure recognised to include variation in environmental and climatic factors, but also likely to feature varied management practices) and chicken line/breed were both found to exert significant impacts (p < 0.05) on caecal microbiome composition. Linear discriminant analysis effect size (LEfSe) revealed 42 breed-specific biomarkers in the chicken lines reared under controlled conditions at two different locations.

CONCLUSION

Chicken breed-specific variation in bacterial occurrence, correlation between genera and clustering of operational taxonomic units indicate scope for quantitative genetic analysis and the possibility of selective breeding of chickens for defined enteric microbiota.

Metagenomic data of DNA viruses of poultry affected with respiratory tract infection

The incidence and severity of respiratory diseases in commercial broiler chicken flocks have increased recently in India because of intensification of the broiler industry. Viral population are predominant in respiratory tract infections and they pose continuous economic burden to poultry industry by causing severe economic losses through decreased productivity [1], [2]. To understand viral metagenome of poultry associated with respiratory infections, we performed DNA virome sequencing and data analysis of broilers from 8 districts of Gujarat State in India. We report high quality sequencing reads and highly abundant DNA viral population present in the infected broiler birds. The raw sequencing data used to perform metagenomic analysis is available in the Sequence Read Archive (SRA) under the BioProject No. PRJNA322592 and Accession No. MAUZ00000000, MAVA00000000, MAVB00000000, MAVC00000000, MAVD00000000, MAVE00000000, MAVF00000000, MAVG00000000 (https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA322592).

Microbiota composition, gene pool and its expression in Gir cattle (Bos indicus) rumen under different forage diets using metagenomic and metatranscriptomic approaches

Zebu (Bos indicus) is a domestic cattle species originating from the Indian subcontinent and now widely domesticated on several continents. In this study, we were particularly interested in understanding the functionally active rumen microbiota of an important Zebu breed, the Gir, under different dietary regimes. Metagenomic and metatranscriptomic data were compared at various taxonomic levels to elucidate the differential microbial population and its functional dynamics in Gir cattle rumen under different roughage dietary regimes. Different proportions of roughage rather than the type of roughage (dry or green) modulated microbiome composition and the expression of its gene pool. Fibre degrading bacteria (i.e. Clostridium, Ruminococcus, Eubacterium, Butyrivibrio, Bacillus and Roseburia) were higher in the solid fraction of rumen (P<0.01) compared to the liquid fraction, whereas bacteria considered to be utilizers of the degraded product (i.e. Prevotella, Bacteroides, Parabacteroides, Paludibacter and Victivallis) were dominant in the liquid fraction (P<0.05). Likewise, expression of fibre degrading enzymes and related carbohydrate binding modules (CBMs) occurred in the solid fraction. When metagenomic and metatranscriptomic data were compared, it was found that some genera and species were transcriptionally more active, although they were in low abundance, making an important contribution to fibre degradation and its further metabolism in the rumen. This study also identified some of the transcriptionally active genera, such as Caldicellulosiruptor and Paludibacter, whose potential has been less-explored in rumen. Overall, the comparison of metagenomic shotgun and metatranscriptomic sequencing appeared to be a much richer source of information compared to conventional metagenomic analysis.

Exploring the deleterious SNPs in XRCC4 gene using computational approach and studying their association with breast cancer in the population of West India

Non-homologous end joining (NHEJ) pathway has pivotal role in repair of double-strand DNA breaks that may lead to carcinogenesis. XRCC4 is one of the essential proteins of this pathway and single-nucleotide polymorphisms (SNPs) of this gene are reported to be associated with cancer risks. In our study, we first used computational approaches to predict the damaging variants of XRCC4 gene. Tools predicted rs79561451 (S110P) nsSNP as the most deleterious SNP. Along with this SNP, we analysed other two SNPs (rs3734091 and rs6869366) to study their association with breast cancer in population of West India. Variant rs3734091 was found to be significantly associated with breast cancer while rs6869366 variant did not show any association. These SNPs may influence the susceptibility of individuals to breast cancer in this population.

Differential expression of microRNAs associated with thermal stress in Frieswal (Bos taurus x Bos indicus) crossbred dairy cattle

Environmental temperature is one of the important abiotic factors that influence the normal physiological function and productive performance of dairy cattle. Temperature stress evokes complex responses that are essential for safeguarding of cellular integrity and animal health. Post-transcriptional regulation of gene expression by miRNA plays a key role cellular stress responses. The present study investigated the differential expression of miRNA in Frieswal (Holstein Friesian × Sahiwal) crossbred dairy cattle that are distinctly adapted to environmental temperature stress as they were evolved by using the temperate dairy breed Holstein Friesian. The results indicated that there was a significant variation in the physiological and biochemical indicators estimated under summer stress. The differential expression of miRNA was observed under heat stress when compared to the normal winter season. Out of the total 420 miRNAs, 65 were differentially expressed during peak summer temperatures. Most of these miRNAs were found to target heat shock responsive genes especially members of heat shock protein (HSP) family, and network analysis revealed most of them having stress-mediated effects on signaling mechanisms. Being greater in their expression profile during peak summer, bta-miR-2898 was chosen for reporter assay to identify its effect on the target HSPB8 (heat shock protein 22) gene in stressed bovine PBMC cell cultured model. Comprehensive understanding of the biological regulation of stress responsive mechanism is critical for developing approaches to reduce the production losses due to environmental heat stress in dairy cattle.

Culture independent assessment of human milk microbial community in lactational mastitis

Breastfeeding undoubtedly provides important benefits to the mother-infant dyad and should be encouraged. Mastitis, one of the common but major cause of premature weaning among lactating women, is an inflammation of connective tissue within the mammary gland. This study reports the influence of mastitis on human milk microbiota by utilizing 16 S rRNA gene sequencing approach. We sampled and sequenced microbiome from 50 human milk samples, including 16 subacute mastitis (SAM), 16 acute mastitis (AM) and 18 healthy-controls. Compared to controls, SAM and AM microbiota were quite distinct and drastically reduced. Genera including, Aeromonas, Staphylococcus, Ralstonia, Klebsiella, Serratia, Enterococcus and Pseudomonas were significantly enriched in SAM and AM samples, while Acinetobacter, Ruminococcus, Clostridium, Faecalibacterium and Eubacterium were consistently depleted. Further analysis of our samples revealed positive aerotolerant odds ratio, indicating dramatic depletion of obligate anaerobes and enrichment of aerotolerant bacteria during the course of mastitis. In addition, predicted functional metagenomics identified several gene pathways related to bacterial proliferation and colonization (e.g. two-component system, bacterial secretion system and motility proteins) in SAM and AM samples. In conclusion, our study confirmed previous hypothesis that mastitis women have lower microbial diversity, increased abundance of opportunistic pathogens and depletion of commensal obligate anaerobes.

Molecular cloning, heterologous expression, and functional characterization of a cellulolytic enzyme (Cel PRII) from buffalo rumen metagenome

A cellulase encoding gene, Cel PRII, was identified from Mehsani buffalo rumen metagenome, and cloned and expressed in Escherichia coli BL21(DE3)pLysS. The 1170 bp full length gene encodes a 389 residue polypeptide (Cel PRII) containing a catalytic domain belonging to glycosyl hydrolase (GH) 5 family. The fusion protein consisting of the Cel PRII, thioredoxin tag and 6x Histidine tag with predicted molecular weight of 63 kDa when recovered from inclusion bodies under denaturing conditions, exhibited cellulolytic activity against carboxymethyl cellulose (CMC). Recombinant Cel PRII was stable in the pH range 4.0-10.0 with pH optima 6.0. The optimal reaction temperature of Cel PRII was 30 °C with more than 50% of its activity retained at the temperatures ranging from 0 to 50 °C. Cel PRII exhibited enhanced enzymatic activity in the presence of Mn²⁺ ions and was inhibited in the presence of chelating agent EDTA. The K _m and V _max values for CMC were found to be 166 mg/mL and 1292 IU/mg, respectively. Cel PRII identified in the present study may act as an excellent candidate for industrial applications, and may aid in lignocellulosic biomass conversion because of its potential cellulolytic activity, thermostability, and excellent pH stability.

Characterization of NPHS2 gene polymorphisms associated to steroid resistance nephrotic syndrome in Indian children

Nephrotic syndrome (NS) is the common glomerular disease in children. These children are treated with steroids, depending upon their behavior. They are either steroid sensitive (SSNS) or steroid resistant (SRNS). NPHS2 gene mutants are linked to the risk of autosomal recessive SRNS and in some cases to SSNS. The present study has been performed to screen single nucleotide polymorphisms (SNPs) of the NPHS2 gene in a group of 90 Indian children suffering with NS (30 SSNS, 30 SRNS and 30 Controls) by PCR method followed by direct exon sequencing. Effect of SNPs on fold expression changes at transcript level of podocin was checked using quantitative real time PCR (qRT-PCR). SNPs identified through sequencing helps to carry out in-silico analysis. Overall 17 SNPs were identified in NPHS2 gene where 6 were found novel. Three missense SNPs p.R299Q, p.P20L and p.G35D were also identified in this population where SNP, p.G35D was found novel. In addition to sequencing analysis, results of in silico analysis shows that a mutant with these three missense SNPs has least ligand binding efficiency compared to native model. Moreover the significant observation of this study included two intronic SNPs c.451+23C>T and c.451+58A>T present in SRNS group of patients. These SNPs has shown high level of clinical significance within genomic and allelic frequency along with haplotypes and linkage disequilibrium count. The qRT-PCR analysis shows, down expression of podocin protein at transcript level in SRNS patients compared to SSNS patients. All these results support the fact that SNPs present in this population could affect the protein structural stability. Thus it is concluded that the polymorphisms predicted in this study might be disease causing in the NPHS2 gene and may have influence on the therapeutic response of NS patients.

Complete mitochondrial genome of threatened mahseer Tor tor (Hamilton 1822) and its phylogenetic relationship within Cyprinidae family

The mahseers (Tor, Neolissochilus and Naziritor) are an important group of fishes endemic to Asia with the conservation status of most species evaluated as threatened. Conservation plans to revive these declining wild populations are hindered by unstable taxonomy. Molecular phylogeny studies with mitochondrial genome have been successfully used to reconstruct the phylogenetic tree and to resolve taxonomic ambiguity. In the present study, complete mitochondrial genome of Tor tor has been sequenced using ion torrent next-generation sequencing platform with coverage of more than 1000 x. Comparative mitogenome analysis shows higher divergence value at ND1 gene than COI gene. Further, occurrence of a distinct genetic lineage of T. tor is revealed. The phylogenetic relationship among mahseer group has been defined as Neolissochilus hexagonolepis ((T. sinensis (T. putitora, T. tor), (T. khudree, T. tambroides)).

Transcriptomic comparison of primary bovine horn core carcinoma culture and parental tissue at early stage

Aim:

Squamous cell carcinoma or SCC of horn in bovines (bovine horn core carcinoma) frequently observed in Bos indicus affecting almost 1% of cattle population. Freshly isolated primary epithelial cells may be closely related to the malignant epithelial cells of the tumor. Comparison of gene expression in between horn’s SCC tissue and its early passage primary culture using next generation sequencing was the aim of this study.

Materials and Methods:

Whole transcriptome sequencing of horn’s SCC tissue and its early passage cells using Ion Torrent PGM were done. Comparative expression and analysis of different genes and pathways related to cancer and biological processes associated with malignancy, proliferating capacity, differentiation, apoptosis, senescence, adhesion, cohesion, migration, invasion, angiogenesis, and metabolic pathways were identified.

Results:

Up-regulated genes in SCC of horn’s early passage cells were involved in transporter activity, catalytic activity, nucleic acid binding transcription factor activity, biogenesis, cellular processes, biological regulation and localization and the down-regulated genes mainly were involved in focal adhesion, extracellular matrix receptor interaction and spliceosome activity.

Conclusion:

The experiment revealed similar transcriptomic nature of horn’s SCC tissue and its early passage cells.

Functional gene profiling through metaRNAseq approach reveals diet-dependent variation in rumen microbiota of buffalo (Bubalus bubalis)

Recent advances in next generation sequencing technology have enabled analysis of complex microbial community from genome to transcriptome level. In the present study, metatranscriptomic approach was applied to elucidate functionally active bacteria and their biological processes in rumen of buffalo (Bubalus bubalis) adapted to different dietary treatments. Buffaloes were adapted to a diet containing 50:50, 75:25 and 100:0 forage to concentrate ratio, each for 6 weeks, before ruminal content sample collection. Metatranscriptomes from rumen fiber adherent and fiber-free active bacteria were sequenced using Ion Torrent PGM platform followed by annotation using MG-RAST server and CAZYmes (Carbohydrate active enzymes) analysis toolkit. In all the samples Bacteroidetes was the most abundant phylum followed by Firmicutes. Functional analysis using KEGG Orthology database revealed Metabolism as the most abundant category at level 1 within which Carbohydrate metabolism was dominating. Diet treatments also exerted significant differences in proportion of enzymes involved in metabolic pathways for VFA production. Carbohydrate Active Enzyme(CAZy) analysis revealed the abundance of genes encoding glycoside hydrolases with the highest representation of GH13 CAZy family in all the samples. The findings provide an overview of the activities occurring in the rumen as well as active bacterial population and the changes occurring through different dietary treatments.

De novo Transcriptome Sequencing to Dissect Candidate Genes Associated with Pearl Millet-Downy Mildew (Sclerospora graminicola Sacc.) Interaction

Understanding the plant-pathogen interactions is of utmost importance to design strategies for minimizing the economic deficits caused by pathogens in crops. With an aim to identify genes underlying resistance to downy mildew, a major disease responsible for productivity loss in pearl millet, transcriptome analysis was performed in downy mildew resistant and susceptible genotypes upon infection and control on 454 Roche NGS platform. A total of ~685 Mb data was obtained with 1 575 290 raw reads. The raw reads were pre-processed into high-quality (HQ) reads making to ~82% with an average of 427 bases. The assembly was optimized using four assemblers viz. Newbler, MIRA, CLC and Trinity, out of which MIRA with a total of 14.10 Mb and 90118 transcripts proved to be the best for assembling reads. Differential expression analysis depicted 1396 and 936 and 1000 and 1591 transcripts up and down regulated in resistant inoculated/resistant control and susceptible inoculated/susceptible control respectively with a common of 3644 transcripts. The pathways for secondary metabolism, specifically the phenylpropanoid pathway was up-regulated in resistant genotype. Transcripts up-regulated as a part of defense response included classes of R genes, PR proteins, HR induced proteins and plant hormonal signaling transduction proteins. The transcripts for skp1 protein, purothionin, V type proton ATPase were found to have the highest expression in resistant genotype. Ten transcripts, selected on the basis of their involvement in defense mechanism were validated with qRT-PCR and showed positive co-relation with transcriptome data. Transcriptome analysis evoked potentials of hypersensitive response and systemic acquired resistance as possible mechanism operating in defense mechanism in pearl millet against downy mildew infection.

VDAP-GUI: a user-friendly pipeline for variant discovery and annotation of raw next-generation sequencing data

Even though next-generation sequencing (NGS) has become an invaluable tool in molecular biology, several laboratories with NGS facilities lack trained Bioinformaticians for data analysis. Here, focusing on the variant detection application of NGS analysis, we have developed a fully automated pipeline, namely Variant Discovery and Annotation Tool-Graphical User Interface (VDAP-GUI), which detects and annotates single nucleotide polymorphisms and insertions/deletions from raw sequence reads. VDAP-GUI consolidates several proven methods in each step such as quality control, trimming, mapping, variant detection and annotation. It supports multiple NGS platforms and has four methodological choices for variant detection. Further, it can re-analyze existing data with alternate thresholds and generates easily interpretable reports in html and tab-delimited formats. Using VDAP-GUI, we have analyzed a publically available human whole-exome sequence dataset. VDAP-GUI is developed using Perl/Tk programming, and is available for free download and use at http://sourceforge.net/projects/vdapgui/.

Mining of Ruminant Microbial Phytase (RPHY1) from Metagenomic Data of Mehsani Buffalo Breed: Identification, Gene Cloning, and Characterization

Phytases have been widely used as animal feed supplements to increase the availability of digestible phosphorus, especially in monogastric animals fed cereal grains. The present study describes the identification of a full-length phytase gene of <i>Prevotella</i> species present in Mehsani buffalo rumen. The gene, designated as RPHY1, consists of 1,251 bp and is expressed into protein with 417 amino acids. A homology search of the deduced amino acid sequence of the RPHY1 phytase gene in a nonredundant protein database showed that it shares 92% similarity with the histidine acid phosphatase domain. Subsequently, the RPHY1 gene was expressed using a pET32a expression vector in <i>Escherichia coli </i>BL21 and purified using a His60 Ni-NTA gravity column. The mass of the purified RPHY1 was estimated to be approximately 63 kDa by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE). The optimal RPHY1 enzyme activity was observed at 55°C (pH 5) and exhibited good stability at 5°C and within the acidic pH range. Significant inhibition of RPHY1 activity was observed for Mg²⁺ and K⁺ metal ions, while Ca²⁺, Mn²⁺, and Na⁺ slightly inhibited enzyme activity. The RPHY1 phytase was susceptible to SDS, and it was highly stimulated in the presence of EDTA. Overall, the observed comparatively high enzyme activity levels and characteristics of the RPHY1 gene mined from rumen prove its promising candidature as a feed supplement enzyme in animal farming.

Isolation and characterization of novel multifunctional recombinant family 26 glycoside hydrolase from Mehsani buffalo rumen metagenome

Rumen microbiota harbor a diverse set of carbohydrate-active enzymes (CAZymes), which play a crucial role in the degradation of a complex plant polysaccharide thereby providing metabolic energy to the host animals. Earlier, we reported CAZYme analysis from the buffalo rumen metagenome by high throughput shotgun sequencing. Among the various CAZymes, glycoside hydrolase family 26 (GH26) enzymes have a number of industrial applications including in paper, oil, biofuel, food, feed, pharmaceutical, coffee, and detergent industries. Here, we report isolation and characterization of GH26 enzyme from the buffalo rumen metagenome. A novel GH26 gene composed of 1,119 base pairs was successfully amplified using the gene-specific primers inferred based on the contig generated from metagenome sequence assembly and cloned in a pET32a (+) expression vector as an N-terminal histidine tag fusion protein. A novel GH26 protein from an unknown rumen microorganism shared a maximum of 68% identity with the Prevotella ruminicola 23 encoded carbohydrate esterase family 7 and 46% with Bacteroides sp. 2_1_33B encoded mannan endo-1, 4-β-mannosidase. The recombinant GH26-histidine tag fusion protein was expressed in Escherichia coli and purified using Ni-NTA affinity chromatography. The purified enzyme displayed multifunctional activities against various carbohydrate substrates including locust bean gum, beechwood xylan, pectin, and carboxymethyl cellulose suggesting mannanase, xylanase, pectin esterase, and endoglucanase activities, respectively.

Effect of roughage on rumen microbiota composition in the efficient feed converter and sturdy Indian Jaffrabadi buffalo (Bubalus bubalis)

BACKGROUND

The rumen microbiota functions as an effective system for conversion of dietary feed to microbial proteins and volatile fatty acids. In the present study, metagenomic approach was applied to elucidate the buffalo rumen microbiome of Jaffrabadi buffalo adapted to varied dietary treatments with the hypothesis that the microbial diversity and subsequent in the functional capacity will alter with diet change and enhance our knowledge of effect of microbe on host physiology. Eight adult animals were gradually adapted to an increasing roughage diet (4 animals each with green and dry roughage) containing 50:50 (J1), 75:25 (J2) and 100:0 (J3) roughage to concentrate proportion for 6 weeks. Metagenomic sequences of solid (fiber adherent microbiota) and liquid (fiber free microbiota) fractions obtained using Ion Torrent PGM platform were analyzed using MG-RAST server and CAZymes approach.

RESULTS

Taxonomic analysis revealed that Bacteroidetes was the most abundant phylum followed by Firmicutes, Fibrobacter and Proteobacteria. Functional analysis revealed protein (25-30 %) and carbohydrate (15-20 %) metabolism as the dominant categories. Principal component analysis demonstrated that roughage proportion, fraction of rumen and type of forage affected rumen microbiome at taxonomic as well as functional level. Rumen metabolite study revealed that rumen fluid nitrogen content reduced in high roughage diet fed animals and pathway analysis showed reduction in the genes coding enzymes involved in methanogenesis pathway. CAZyme annotation revealed the abundance of genes encoding glycoside hydrolases (GH), with the GH3 family most abundant followed by GH2 and GH13 in all samples.

CONCLUSIONS

Results reveals that high roughage diet feed improved microbial protein synthesis and reduces methane emission. CAZyme analysis indicated the importance of microbiome in feed component digestion for fulfilling energy requirements of the host. The findings help determine the role of rumen microbes in plant polysaccharide breakdown and in developing strategies to maximize productivity in ruminants.

Deep insights into carbohydrate metabolism in the rumen of Mehsani buffalo at different diet treatments

Ruminants are known to harbour a vast and diverse microbial community that functions in utilizing the fibrous and starchy feedstuffs. The microbial fermentation of fibrous and starchy feed is carried out by different groups of microbiota, which function in synergistic mechanism. The exploration of the shift in carbohydrate utilizing microbial community with the change in diet will reveal the efficient role of that group of microbial community in particular carbohydrate utilization. The present study explains the shifts in microbial enzymes for carbohydrate utilization with the change in the feed proportions and its correlation with the microbial community abundance at that particular treatment. The sequencing data of the present study is submitted to NCBI SRA with experiment accession IDs (ERX162128, ERX162129, ERX162130, ERX162131, ERX162139, ERX162134, ERX162140, ERX162141, ERX197218, ERX197219, ERX197220, ERX197221, ERX162158, ERX162159, ERX162160, ERX162161, ERX162176, ERX162164, ERX162165, ERX162166, ERX162167, ERX162168, ERX162169, ERX162177).

Identification of putative SNPs in progressive retinal atrophy affected Canis lupus familiaris using exome sequencing

Progressive retinal atrophy (PRA) is one of the major causes of retinal photoreceptor cell degeneration in canines. The inheritance pattern of PRA is autosomal recessive and genetically heterogeneous. Here, using targeted sequencing technology, we have performed exome sequencing of 10 PRA-affected (Spitz=7, Cocker Spaniel=1, Lhasa Aphso=1 and Spitz-Labrador cross breed=1) and 6 normal (Spitz=5, Cocker Spaniel=1) dogs. The high-throughput sequencing using 454-Roche Titanium sequencer generated about 2.16 Giga bases of raw data. Initially, we have successfully identified 25,619 single nucleotide polymorphisms (SNPs) that passed the stringent SNP calling parameters. Further, we performed association study on the cohort, and the highly significant (0.001) associations were short-listed and investigated in-depth. Out of the 171 significant SNPs, 113 were previously unreported. Interestingly, six among them were non-synonymous coding (NSC) SNPs, which includes CPPED1 A>G (p.M307V), PITRM1 T>G (p.S715A), APP G>A (p.T266M), RNF213 A>G (p.V1482A), C>A (p.V1456L), and SLC46A3 G>A (p.R168Q). On the other hand, 35 out of 113 unreported SNPs were falling in regulatory regions such as 3'-UTR, 5'-UTR, etc. In-depth bioinformatics analysis revealed that majority of NSC SNPs have damaging effect and alter protein stability. This study highlighted the genetic markers associated with PRA, which will help to develop genetic assay-based screening in effective breeding.

Isolation of chitinolytic Clostridium sp. NCR from Mehsani buffalo rumen, its genomic analysis and potential role in rumen

Genomic analysis of Clostridium sp. NCR, an anaerobic Gram positive bacterium which was isolated from rumen fluid of Mehsani breed of buffalo revealed presence of various environmental gene tags (EGTs) involved in pathways for utilizing a wide range of substrates. Here we report the sequence of this rumen isolate, its whole genome sequence has been deposited in DDBJ/EMBL/GenBank under the accession number JQHY00000000. The genome comprises of a 3.62-Mb draft genome with a G + C content of 28.10%, which encodes a total of 3126 proteins. Functional analysis provides information about the microbe's role in maintaining host homeostasis and its fiber degradation potential.

Geriatric Respondents and Non-Respondents to Probiotic Intervention Can be Differentiated by Inherent Gut Microbiome Composition

Scope

Probiotic interventions are known to have been shown to influence the composition of the intestinal microbiota in geriatrics. The growing concern is the apparent variation in response to identical strain dosage among human volunteers. One factor that governs this variation is the host gut microbiome. In this study, we attempted to define a core gut metagenome, which could act as a predisposition signature marker of inherent bacterial community that can help predict the success of a probiotic intervention.

Methods and results

To characterize the geriatric gut microbiome, we designed primers targeting the 16S rRNA hypervariable region V2–V3 followed by semiconductor sequencing using Ion Torrent PGM. Among respondents and non-respondents, the chief genera of phylum Firmicutes that showed significant differences are Lactobacillus, Clostridium, Eubacterium, and Blautia (q < 0.002), while in the genera of phylum Proteobacteria included Shigella, Escherichia, Burkholderia and Camphylobacter (q < 0.002).

Conclusion

We have identified potential microbial biomarkers and taxonomic patterns that correlate with a positive response to probiotic intervention in geriatric volunteers. Future work with larger cohorts of geriatrics with diverse dietary influences could reveal the potential of the signature patterns of microbiota for personalized nutrition.

Functional Characterization Reveals Novel Putative Coding Sequences in Prevotella ruminicola Genome Extracted from Rumen Metagenomic Studies

AIM

To reassemble Prevotella ruminicola genome from rumen metagenomic data of cattle and buffalo and compare with the published reference genome.

METHOD

Rumen microbial communities from Mehsani buffaloes (n = 8) and Kankrej cattle (n = 8), each adapted to different proportions of a dry or green roughage diet, were subjected to metagenomic sequencing by Ion Torrent PGM, and subsequent reads were analyzed by MG-RAST. Using reference-guided assembly of the sequences against the published P. ruminicola strain 23, draft genomes of 2.56 and 2.46 Mb were reconstructed from Mehsani buffalo and Kankrej cows, respectively. The genomes were annotated using the RAST Server and carbohydrate active enzyme (CAZyme) analysis.

RESULTS

Taxonomic analysis by MG-RAST revealed P. ruminicola to be the most abundant species present among the rumen microflora. Functional annotation of reconstructed genomes using the RAST Server depicted the maximum assignment of coding sequences involved in the subsystems amino acid and derivatives and carbohydrate metabolism. CAZyme profiling revealed the glycoside hydrolases (GH) family to be the most abundant. GH family subclassification revealed that the extracted genomes had more sequence hits for GH2, GH3, GH92 and GH97 as compared to the reference.

CONCLUSION

The results reflect the metabolic significance of rumen-adapted P. ruminicola in utilizing a coarse diet for animals based on acquisition of novel genetic elements.

Low-depth shotgun sequencing resolves complete mitochondrial genome sequence of Labeo rohita

Labeo rohita, popularly known as rohu, is a widely cultured species in whole Indian subcontinent. In the present study, we used in-silico approach to resolve complete mitochondrial genome of rohu. Low-depth shotgun sequencing using Roche 454 GS FLX (Branford, Connecticut, USA) followed by de novo assembly in CLC Genomics Workbench version 7.0.4 (Aarhus, Denmark) revealed the complete mitogenome of L. rohita to be 16 606 bp long (accession No. KR185963). It comprised of 13 protein-coding genes, 22 tRNAs, 2 rRNAs and 1 putative control region. The gene order and organization are similar to most vertebrates. The mitogenome in the present investigation has 99% similarity with that of previously reported mitogenomes of rohu and this is also evident from the phylogenetic study using maximum-likelihood (ML) tree method. This study was done to determine the feasibility, accuracy and reliability of low-depth sequence data obtained from NGS platform as compared to the Sanger sequencing. Thus, NGS technology has proven to be competent and a rapid in-silico alternative to resolve the complete mitochondrial genome sequence, thereby reducing labors and time.

In silico analysis of consequences of non-synonymous SNPs of Slc11a2 gene in Indian bovines

The aim of our study was to analyze the consequences of non-synonymous SNPs in Slc11a2 gene using bioinformatic tools. There is a current need of efficient bioinformatic tools for in-depth analysis of data generated by the next generation sequencing technologies. SNPs are known to play an imperative role in understanding the genetic basis of many genetic diseases. Slc11a2 is one of the major metal transporter families in mammals and plays a critical role in host defenses. In this study, we performed a comprehensive analysis of the impact of all non-synonymous SNPs in this gene using multiple tools like SIFT, PROVEAN, I-Mutant and PANTHER. Among the total 124 SNPs obtained from amplicon sequencing of Slc11a2 gene by Ion Torrent PGM involving 10 individuals of Gir cattle and Murrah buffalo each, we found 22 non-synonymous. Comparing the prediction of these 4 methods, 5 nsSNPs (G369R, Y374C, A377V, Q385H and N492S) were identified as deleterious. In addition, while tested out for polar interactions with other amino acids in the protein, from above 5, Y374C, Q385H and N492S showed a change in interaction pattern and further confirmed by an increase in total energy after energy minimizations in case of mutant protein compared to the native.

•

22 nsSNPs were predicted to decrease the stability of protein based on I-Mutant.

•

From these SNPs, 5 was identified as deleterious by SIFT, PROVEAN, and PANTHER.

•

Y374C, Q385H and N492S were found to be damaging.

Genetic polymorphism of toll-like receptors 4 gene by polymerase chain reaction-restriction fragment length polymorphisms, polymerase chain reaction-single-strand conformational polymorphism to correlate with mastitic cows

Aim:

An attempt has been made to study the toll-like receptors 4 (TLR4) gene polymorphism from cattle DNA to correlate with mastitis cows.

Materials and Methods:

In present investigation, two fragments of TLR4 gene named T4CRBR1 and T4CRBR2 of a 316 bp and 382 bp were amplified by polymerase chain reaction (PCR), respectively from Kankrej (22) and Triple cross (24) cattle. The genetic polymorphisms in the two populations were detected by a single-strand conformational polymorphism in the first locus and by digesting the fragments with restriction endonuclease Alu I in the second one.

Results:

Results showed that both alleles (A and B) of two loci were found in all the two populations and the value of polymorphism information content indicated that these were highly polymorphic. Statistical results of χ² test indicated that two polymorphism sites in the two populations fit with Hardy–Weinberg equilibrium (p<0.05). Meanwhile, the effect of polymorphism of TLR4 gene on the somatic cell score (SCS) indicated the cattle with allele a in T4CRBR1 showed lower SCS than that of allele B (p<0.05). Thus, the allele A might play an important role in mastitis resistance in cows.

Conclusion:

The relationship between the bovine mastitis trait and the polymorphism of TLR4 gene indicated that the bovine TLR4 gene may play an important role in mastitis resistance.

Analysis of consequences of non-synonymous SNP in feed conversion ratio associated TGF-β receptor type 3 gene in chicken

The recent advances in high throughput sequencing technology accelerate possible ways for the study of genome wide variation in several organisms and associated consequences. In the present study, mutations in TGFBR3 showing significant association with FCR trait in chicken during exome sequencing were further analyzed. Out of four SNPs, one nsSNP p.Val451Leu was found in the coding region of TGFBR3. In silico tools such as SnpSift and PANTHER predicted it as deleterious (0.04) and to be tolerated, respectively, while I-Mutant revealed that protein stability decreased. The TGFBR3 I-TASSER model has a C-score of 0.85, which was validated using PROCHECK. Based on MD simulation, mutant protein structure deviated from native with RMSD 0.08 Å due to change in the H-bonding distances of mutant residue. The docking of TGFBR3 with interacting TGFBR2 inferred that mutant required more global energy. Therefore, the present study will provide useful information about functional SNPs that have an impact on FCR traits.

•

Investigated functional nsSNP p.Val451Leu (rs312979494) in feed conversion ratio (FCR) associated TGFBR3 of chicken

•

Computational tools (SIFT and I-Mutant 2.0) predicted that this nsSNP was deleterious.

•

Mutant structure of TGFBR3 showed high energies and RMS deviations compared to native using MD simulation.

•

Molecular docking of TGFBR3 with interacting protein TGFBR2 showed an increase in global energy of mutant compared to native.

•

We have predicted that functional SNP has an impact on TGFBR3 of chicken and thus can be treated as candidate SNP in screening.

Deciphering the microbiota of Tuwa hot spring, India using shotgun metagenomic sequencing approach

Here, we report metagenome from the Tuwa hot spring, India using shotgun sequencing approach. Metagenome consisted of 541,379 sequences with 98.7 Mbps size with 46% G + C content. Metagenomic sequence reads were deposited into the EMBL database under accession number ERP009321. Community analysis presented 99.1% sequences belong to bacteria, 0.3% of eukaryotic origin, 0.2% virus derived and 0.05% from archea. Unclassified and unidentified sequences were 0.4% and 0.07% respectively. A total of 22 bacterial phyla include 90 families and 201 species were observed in the hot spring metagenome. Firmicutes (97.0%), Proteobacteria (1.3%) and Actinobacteria (0.4%) were reported as dominant bacterial phyla. In functional analysis using Cluster of Orthologous Group (COG), 21.5% drops in the poorly characterized group. Using subsystem based annotation, 4.0% genes were assigned for stress responses and 3% genes were fit into the metabolism of aromatic compounds. The hot spring metagenome is very rich with novel sequences affiliated to unclassified and unidentified lineages, suggesting the potential source for novel microbial species and their products.

Genomic analysis of host - Peste des petits ruminants vaccine viral transcriptome uncovers transcription factors modulating immune regulatory pathways

Peste des petits ruminants (PPR), is an acute transboundary viral disease of economic importance, affecting goats and sheep. Mass vaccination programs around the world resulted in the decline of PPR outbreaks. Sungri 96 is a live attenuated vaccine, widely used in Northern India against PPR. This vaccine virus, isolated from goat works efficiently both in sheep and goat. Global gene expression changes under PPR vaccine virus infection are not yet well defined. Therefore, in this study we investigated the host-vaccine virus interactions by infecting the peripheral blood mononuclear cells isolated from goat with PPRV (Sungri 96 vaccine virus), to quantify the global changes in the transcriptomic signature by RNA-sequencing. Viral genome of Sungri 96 vaccine virus was assembled from the PPRV infected transcriptome confirming the infection and demonstrating the feasibility of building a complete non-host genome from the blood transcriptome. Comparison of infected transcriptome with control transcriptome revealed 985 differentially expressed genes. Functional analysis showed enrichment of immune regulatory pathways under PPRV infection. Key genes involved in immune system regulation, spliceosomal and apoptotic pathways were identified to be dysregulated. Network analysis revealed that the protein - protein interaction network among differentially expressed genes is significantly disrupted in infected state. Several genes encoding TFs that govern immune regulatory pathways were identified to co-regulate the differentially expressed genes. These data provide insights into the host - PPRV vaccine virus interactome for the first time. Our findings suggested dysregulation of immune regulatory pathways and genes encoding Transcription Factors (TFs) that govern these pathways in response to viral infection.

Exploring genetic polymorphism in innate immune genes in Indian cattle (Bos indicus) and buffalo (Bubalus bubalis) using next generation sequencing technology

Activation of innate immunity initiates various cascades of reactions that largely contribute to defense against physical, microbial or chemical damage, prompt for damage repair and removal of causative organisms as well as restoration of tissue homeostasis. Genetic polymorphism in innate immune genes plays prominent role in disease resistance capabilities in various breeds of cattle and buffalo. Here we studied single nucleotide variations (SNP/SNV) and haplotype structure in innate immune genes viz CHGA, CHGB, CHGC, NRAMP1, NRAMP2, DEFB1, BNBD4, BNBD5, TAP and LAP in Gir cattle and Murrah buffalo. Targeted sequencing of exonic regions of these genes was performed by Ion Torrent PGM sequencing platform. The sequence reads obtained corresponding to coding regions of these genes were mapped to reference genome of cattle BosTau7 by BWA program using genome analysis tool kit (GATK). Further variant analysis by Unified Genotyper revealed 54 and 224 SNPs in Gir and Murrah respectively and also 32 SNVs was identified. Among these SNPs 43, 36, 11,32,81,21 and 22 variations were in CHGA, CHGB, CHGC, NRAMP1, NRAMP2, DEFB1 and TAP genes respectively. Among these identified 278 SNPs, 24 were found to be reported in the dbSNP database. Variant analysis was followed by structure formation of haplotypes based on multiple SNPs using SAS software revealed a large number of haplotypes. The SNP discovery in innate immune genes in cattle and buffalo breeds of India would advance our understanding of role of these genes in determining the disease resistance/susceptibility in Indian breeds. The identified SNPs and haplotype data would also provide a wealth of sequence information for conservation studies, selective breeding and designing future strategies for identifying disease associations involving samples from distinct populations.

•

Sequencing of exonic regions of 10 selected innate immune genes, 895.60 Mb data was generated.

•

Variant analysis using GATK pipeline revealed 278 SNPs, 32 SNVs and from these 24 were found to be reported.

•

43, 36 and 11 SNPs were observed in chromagranin genes viz. CHGA, CHGB and CHGC respectively.

•

32 and 81 SNPs were identified in NRAMP genes viz. NRAMP1 and NRAMP2 respectively.

•

21 and 22 SNPs were observed in β-defensin genes viz. DEFB1 and TAP respectively.

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.

Evaluation of a topical herbal drug for its in-vivo immunomodulatory effect on cytokines production and antibacterial activity in bovine subclinical mastitis

BACKGROUND

Antibiotics have been in use in the treatment of bovine mastitis since decades; however, their use is associated with cost issues and human health concern. Use of herbal drugs does not generally carry these disadvantages. Many plants/herbs have been evaluated in the treatment of bovine mastitis with additional property of immunomodulation in affected mammary gland.

AIM

To evaluate a topical herbal drug in two breeds of cattle for its in-vivo immunomodulatory effect on cytokines production and antibacterial activity in bovine subclinical mastitis.

MATERIALS AND METHODS

The response to treatment was evaluated by enumerating somatic cell count (SCC), determining total bacterial load, and studying the expression of different cytokines (interleukin [IL]-6, IL-8, IL-12, granulocyte macrophage-colony stimulating factor, interferon (IFN)-γ and tumor necrosis factor [TNF]-α).

RESULTS

The pre- and post-treatment SCC in mastitic quarters statistically did not differ significantly, however, total bacterial load declined significantly from day 0 onwards in both the breeds. Highly significant differences (P < 0.01) were observed in all the cytokines on day 0, 5, and 21 postlast treatment in both the breeds. The expression level of all the cytokines showed a significant increase on day 5, while a decrease was noticed on day 21 in both the breeds of cattle. The comparison of cytokine expression profiles between crossbred and Gir cattle revealed a significant difference in expression of IL-6 and TNF-α. However, other cytokines exhibited a similar pattern of expression in both breeds, which was non-significant.

CONCLUSION

The topical herbal drug exhibited antibacterial and immunomodulatory activities in subclinical mastitis and thus the work supports its use as alternative herbal therapy against subclinical udder infection in bovines.