Development, validation and application of single molecule molecular inversion probe based novel integrated genetic screening method for 29 common lysosomal storage disorders in India.

BACKGROUND: Current clinical diagnosis pathway for lysosomal storage disorders (LSDs) involves sequential biochemical enzymatic tests followed by DNA sequencing, which is iterative, has low diagnostic yield and is costly due to overlapping clinical presentations. Here, we describe a novel low-cost and high-throughput sequencing assay using single-molecule molecular inversion probes (smMIPs) to screen for causative single nucleotide variants (SNVs) and copy number variants (CNVs) in genes associated with 29 common LSDs in India. RESULTS: 903 smMIPs were designed to target exon and exon-intron boundaries of targeted genes (n = 23; 53.7 kb of the human genome) and were equimolarly pooled to create a sequencing library. After extensive validation in a cohort of 50 patients, we screened 300 patients with either biochemical diagnosis (n = 187) or clinical suspicion (n = 113) of LSDs. A diagnostic yield of 83.4% was observed in patients with prior biochemical diagnosis of LSD. Furthermore, diagnostic yield of 73.9% (n = 54/73) was observed in patients with high clinical suspicion of LSD in contrast with 2.4% (n = 1/40) in patients with low clinical suspicion of LSD. In addition to detecting SNVs, the assay could detect single and multi-exon copy number variants with high confidence. Critically, Niemann-Pick disease type C and neuronal ceroid lipofuscinosis-6 diseases for which biochemical testing is unavailable, could be diagnosed using our assay. Lastly, we observed a non-inferior performance of the assay in DNA extracted from dried blood spots in comparison with whole blood. CONCLUSION: We developed a flexible and scalable assay to reliably detect genetic causes of 29 common LSDs in India. The assay consolidates the detection of multiple variant types in multiple sample types while having improved diagnostic yield at same or lower cost compared to current clinical paradigm.

Interplay of gene expression and regulators under salinity stress in gill of Labeo rohita.

BACKGROUND: Labeo rohita is the most preferred freshwater carp species in India. The concern of increasing salinity concentration in freshwater bodies due to climate change may greatly impact the aquatic environment. Gills are one of the important osmoregulatory organs and have direct contact with external environment. Hence, the current study is conducted to understand the gill transcriptomic response of L. rohita under hypersalinity environment. RESULTS: Comprehensive analysis of differentially expressed long non-coding RNAs (lncRNAs), microRNAs (miRNAs) and mRNAs was performed in gills of L. rohita treated with 2, 4, 6 and 8ppt salinity concentrations. Networks of lncRNA-miRNA-mRNA revealed involvement of 20, 33, 52 and 61 differentially expressed lncRNAs, 11, 13, 26 and 21 differentially expressed miRNAs in 2, 4, 6 and 8ppt groups between control and treatment respectively. These lncRNA-miRNA pairs were regulating 87, 214, 499 and 435 differentially expressed mRNAs (DE mRNAs) in 2, 4, 6 and 8ppt treatments respectively. Functional analysis of these genes showed enrichment in pathways related to ion transportation and osmolyte production to cope with induced osmotic pressure due to high salt concentration. Pathways related to signal transduction (MAPK, FOXO and phosphatidylinositol signaling), and environmental information processing were also upregulated under hypersalinity. Energy metabolism and innate immune response pathways also appear to be regulated. Protein turnover was high at 8ppt as evidenced by enrichment of the proteasome and aminoacyl tRNA synthesis pathways, along with other enriched KEGG terms such as apoptosis, cellular senescence and cell cycle. CONCLUSION: Altogether, the RNA-seq analysis provided valuable insights into competitive endogenous (lncRNA-miRNA-mRNA) regulatory network of L. rohita under salinity stress. L. rohita is adapting to the salinity stress by means of upregulating protein turnover, osmolyte production and removing the damaged cells using apoptotic pathway and regulating the cell growth and hence diverting the essential energy for coping with salinity stress.

Potential of camel rumen derived Bacillus subtilis and Bacillus velezensis strains for application in plant biomass hydrolysis.

Rumen inhabiting Bacillus species possesses a high genetic potential for plant biomass hydrolysis and conversion to value-added products. In view of the same, five camel rumen-derived Bacillus strains, namely B. subtilis CRN 1, B. velezensis CRN 2, B. subtilis CRN 7, B. subtilis CRN 11, and B. velezensis CRN 23 were initially assayed for diverse hydrolytic activities, followed by genome mining to unravel the potential applications. CRN 1 and CRN 7 showed the highest endoglucanase activity with 0.4 U/ml, while CRN 23 showed high ß-xylosidase activity of 0.36 U/ml. The comprehensive genomic insights of strains resolve taxonomic identity, clusters of an orthologous gene, pan-genome dynamics, and metabolic features. Annotation of Carbohydrate active enzymes (CAZymes) reveals the presence of diverse glycoside hydrolases (GH) GH1, GH5, GH43, and GH30, which are solely responsible for the effective breakdown of complex bonds in plant polysaccharides. Further, protein modeling and ligand docking of annotated endoglucanases showed an affinity for cellotrioside, cellobioside, and ß-glucoside. The finding indicates the flexibility of Bacillus-derived endoglucanase activity on diverse cellulosic substrates. The presence of the butyrate synthesis gene in the CRN 1 strain depicts its key role in the production of important short-chain fatty acids essential for healthy rumen development. Similarly, antimicrobial peptides such as bacilysin and non-ribosomal peptides (NRPS) synthesized by the Bacillus strains were also annotated in the genome. The findings clearly define the role of Bacillus sp. inside the camel rumen and its potential application in various plant biomass utilizing industry and animal health research sectors.

Genomic profiling and characteristics of a C1 degrading heterotrophic fresh-water bacterium Paracoccus sp. strain DMF.

Paracoccus species are metabolically versatile gram-negative, aerobic facultative methylotrophic bacteria showing enormous promise for environmental and bioremediation studies. Here we report, the complete genome analysis of Paracoccus sp. strain DMF (P. DMF) that was isolated from a domestic wastewater treatment plant in Kanpur, India (26.4287 °N, 80.3891 °E) based on its ability to degrade a recalcitrant organic solvent N, N-dimethylformamide (DMF). The results reveal a genome size of 4,202,269 base pairs (bp) with a G + C content of 67.9%. The assembled genome comprises 4141 coding sequences (CDS), 46 RNA sequences, and 2 CRISPRs. Interestingly, catabolic operons related to the conventional marine-based methylated amines (MAs) degradation pathway were functionally annotated within the genome of an obligated aerobic heterotroph that is P. DMF. The genomic data-based characterization presented here for the novel heterotroph P. DMF aims to improve the understanding of the phenotypic gene products, enzymes, and pathways involved with greater emphasis on facultative methylotrophic motility-based latent pathogenicity.

Nasopharyngeal microbiome of COVID-19 patients revealed a distinct bacterial profile in deceased and recovered individuals.

The bacterial co-infections in SARS-CoV-2 patients remained the least explored subject of clinical manifestations that may also determine the disease severity. Nasopharyngeal microbial community structure within SARS-CoV-2 infected patients could reveal interesting microbiome dynamics that may influence the disease outcomes. Here, in this research study, we analyzed distinct nasopharyngeal microbiome profile in the deceased (n = 48) and recovered (n = 29) COVID-19 patients and compared it with control SARS-CoV-2 negative individuals (control) (n = 33). The nasal microbiome composition of the three groups varies significantly (PERMANOVA, p-value <0.001), where deceased patients showed higher species richness compared to the recovered and control groups. Pathogenic genera, including Corynebacterium (LDA score 5.51), Staphylococcus, Serratia, Klebsiella and their corresponding species were determined as biomarkers (p-value <0.05, LDA cutoff 4.0) in the deceased COVID-19 patients. Ochrobactrum (LDA score 5.79), and Burkholderia (LDA 5.29), were found in the recovered group which harbors ordinal bacteria (p-value <0.05, LDA-4.0) as biomarkers. Similarly, Pseudomonas (LDA score 6.19), and several healthy nasal cavity commensals including Veillonella, and Porphyromonas, were biomarkers for the control individuals. Healthy commensal bacteria may trigger the immune response and alter the viral infection susceptibility and thus, may play important role and possible recovery that needs to be further explored. This research finding provide vital information and have significant implications for understanding the microbial diversity of COVID-19 patients. However, additional studies are needed to address the microbiome-based therapeutics and diagnostics interventions.

Rumen and fecal microbial profiles in cattle fed high lignin diets using metagenome analysis.

Coconut coir (a lignin-rich, organic material) is widely used for its commercial and economic benefits. In this study, crossbred (exotic) and Kankrej (indigenous) breeds of cattle were fed diets containing 7 or 14% coconut coir. Metagenomic analyses (16S rRNA gene amplicon and shotgun sequencing) were used to characterize the microbial community in the rumen and fecal samples along with their functional capabilities. Both amplicon and shotgun analyses revealed the predominance of bacterial phyla, Bacteroidetes, Firmicutes, Actinobacteria and Fibrobacter in ruminal liquid, ruminal solid and fecal samples. 16S rRNA gene amplicon sequencing revealed a total of 18 different bacterial taxa were found to be enriched exclusively in the animals fed with 14% coir. The shotgun analysis revealed abundance of bacterial genera, Fibrobacter, Clostridium, Prevotella, Butyrivibrio, and Ruminococcus in both liquid and solid fractions of ruminal contents, while in the fecal sample, Bacteroides, Alistipes, Plaudibacter, Parabacteroides, Porphyromonas, and Victivallis and archaeal genus, Methanocorpusculum were abundant. The functional analysis based on dbCAN database suggested that among the Glycoside hydrolase family, genes that encode oligosaccharide degrading enzymes, GH3, GH13 (p-value < 0.05), and GH43 were abundant in the feces. In ruminal solid, cellulase encoding the GH5 family was abundant. Also, lignocellulosic binding modules encoded by the CBM family, including cellulose (CBM3) and hemicellulose binding modules (CBM32 and CBM67) were abundant. Thus, the study indicated the enrichment of lignocellulosic enzymes in ruminal contents in response to feeding the coconut coir, which could be mined for potential biofuel production and other biotechnological applications.

Metagenomic analysis of bacterial communities during estrous cycle in Bos indicus.

The present study, using 16 s rRNA sequencing of the V3-V4 hypervariable region, was aimed to check diversity of vaginal microbiota throughout different stages of the estrous cycle in Bos indicus, with attention to changes in progesterone hormone and microorganism diversity. Metagenomic research was conducted on vaginal swabs obtained from nine healthy Indian Gir cows at estrus (day 0), metestrus (day 04), diestrus (day 12), and proestrus (day 16) phases of the estrous cycle. The findings revealed that the diestrus phase has a different bacterial diversity than the other three estrous cycle phases, implying that progesterone affects bacterial diversity. Proteobacteria and Firmicutes were the most abundant phyla at the phylum level, accounting for 94% of bacterial diversity. Actinobacteriota, Patescibacteria, Cyanobacteria, and Bacteroidota were among the less prevalent phyla observed in all samples. After statistical analysis, Bacillaceae, Alcaligenes, Enterobacteriaceae, and Morganellaceae families were more significant. The Enterobacteriaceae family was found to be lower in the diestrus phase than in the other three phases; in contrast, all statistically significant genera were high at the diestrus phase. The luteal stage had higher levels of Micrococcus, Stenotrophomonas, UGC-010, Massilia, and Methylobacillus than the follicular stage, and statistical analysis revealed substantial difference between the luteal and follicular stages. Lactobacillus genus was present in both the estrus and diestrus phases. This study represents an important step toward the understanding of microbial diversity within different stages of the estrous cycle of Indian cows.

Unraveling the camel rumen microbiome through metaculturomics approach for agriculture waste hydrolytic potential.

Cellulose is the most abundant natural polymer present on Earth in the form of agriculture waste. Hydrolysis of agriculture waste for simple fermentable reducing sugars is the bottleneck in the area of biofuel generation and other value-added products. The present study aims to utilize the camel rumen as a bioreactor for potent cellulolytic and hemicellulolytic bacteria by altering the feed types with varying cellulosic concentrations. A total of 6716 bacterial cultures were subjected to three layers of screening, where plate zymography and chromophoric substrate screening served as primary screening method for cellulolytic and hemicellulolytic potential. The potential isolates were genetically grouped using RAPD, and 51 representative isolates from each group were subjected to molecular identification through 16S rDNA sequencing, followed by quantification of various cellulolytic and hemicellulolytic enzymes. Out of 51 potent isolates, 5 isolates had high endoglucanase activity ranging from 0.3 to 0.48 U/ml. The selected five key isolates identified as Pseudomonas, Paenibacillus, Citrobacter, Bacillus subtilis, and Enterobacter were employed for hydrolyzing the various agriculture residues and resulted in approximately 0.4 mg/ml of reducing sugar. Furthermore, the metaculturomics approach was implemented to deduce the total cultured diversity through 16S rRNA amplicon library sequencing. The metaculturomics data revealed the dominance of proteobacteria and unidentified bacterial population in all four feed types, which indicates the possibility of culturing novel cellulose-deconstructing bacteria. Moreover, the presence of diverse hydrolytic enzymes in cultured isolates supports the usage of these bacteria in bio-processing of agriculture waste residues and obtaining the biofuels and other value-added products.

Isolation, purification and characterization of a novel esterase from camel rumen metagenome.

Bacterial esterases are gaining the importance in pharmaceuticals and agrochemical industries due to their excellent biocatalytic properties and a wide range of applications. In the present study, a novel gene encoding an esterase (designated as Est-CR) was identified from shotgun metagenomic sequencing data of camel rumen (Camelus dromedarius) liquor. The open reading frame consisted of 1,224bp, which showed 84.03% sequence identity to Bacteroidales bacterium, corresponding to a protein of 407 amino acids and has a catalytic domain belonging to an esterase. Est-CR belonged to family V with GLSMG domain. The purified enzyme with a molecular mass of 62.64 kDa was checked on SDS-PAGE, and its expression was confirmed by western blotting. The enzyme was active and stable over a broad range of temperature (35-65 °C), displayed the maximum activity at 50 °C and pH 7.0. Individually all metal ions inhibited the enzyme activity, while in combination, K2+, Ca2+, Mg2+ and Mn2+ metal ions enhanced the enzyme activity. The detergents strongly inhibited the activity, while EDTA (10 mM) increased the activity of the Est-CR enzyme. The enzyme showed specificity to short-chain substrates and displayed an optimum activity against butyrate ester. This novel enzyme might serve as a promising candidate to meet some harsh industrial processes enzymatic needs.

Mutational landscape of TRPC6, WT1, LMX1B, APOL1, PTPRO, PMM2, LAMB2 and WT1 genes associated with Steroid resistant nephrotic syndrome.

BACKGROUND: Nephrotic syndrome appears as a group of symptoms like proteinuria, edema and hyperlipidemia. Identification of monogenic forms revealed the physiology and pathogenesis of the SRNS. METHODS AND RESULTS: We performed Illumina panel sequencing of seven genes in 90 Indian patients to determine the role of these genetic mutations in nephrotic syndrome prognosis. Samtool was used for variants calling, and SnpEff and Snpsift did variants annotation. Clinical significance and variant classification were performed by the ClinVar database. In SSNS and SRNS patients, we found 0.78% pathogenic and 3.41% likely pathogenic mutations. Pathogenic mutations were found in LAMB2, LMX1B and WT1 genes, while likely pathogenic mutations were found in (6/13) LAMB2, (2/13) LMX1B, (2/13) TRPC6, (2/13) PTPRO and (1/13) PMM2 genes. Approximately 46% likely pathogenic mutations were contributed to the LAMB2 gene in SSNS and SRNS patients. We also detect 30 VUS (variants of uncertain significance), which were found (17/30) pathogenic and (13/30) likely pathogenic by different prediction tools. CONCLUSIONS: Multigene panels were used for genetic screening of heterogeneous disorders like nephrotic syndrome in the Indian population. We found pathogenic, likely pathogenic and certain VUS, which were responsible for the pathogenesis of the disease. Therefore, mutational analysis of SSNS and SRNS is necessary to avoid adverse effects of corticosteroids, modify the intensity of immunosuppressing agents, and prevent the disease's progression.

Unravelling the early warning capability of wastewater surveillance for COVID-19: A temporal study on SARS-CoV-2 RNA detection and need for the escalation.

Wastewater-based Epidemiological (WBE) surveillance offers a promising approach to assess the pandemic situation covering pre-symptomatic and asymptomatic cases in highly populated area under limited clinical tests. In the present study, we analyzed SARS-CoV-2 RNA in the influent wastewater samples (n = 43) from four wastewater treatment plants (WWTPs) in Gandhinagar, India, during August 7th to September 30th, 2020. A total of 40 samples out of 43 were found positive i.e. having at least two genes of SARS-CoV-2. The average Ct values for S, N, and ORF 1 ab genes were 32.66, 33.03, and 33.95, respectively. Monthly variation depicted a substantial rise in the average copies of N (~120%) and ORF 1 ab (~38%) genes in the month of September as compared to August, while S-gene copies declined by 58% in September 2020. The SARS-CoV-2 genome concentration was higher in the month of September (~924.5 copies/L) than August (~897.5 copies/L), corresponding to a ~2.2-fold rise in the number of confirmed cases during the study period. Further, the percentage change in genome concentration level on a particular date was found in the lead of 1-2 weeks of time with respect to the official confirmed cases registered based on clinical tests on a temporal scale. The results profoundly unravel the potential of WBE surveillance to predict the fluctuation of COVID-19 cases to provide an early warning. Our study explicitly suggests that it is the need of hour that the wastewater surveillance must be included as an integral part of COVID-19 pandemic monitoring which can not only help the water authorities to identify the hotspots within a city but can provide up to 2 weeks of time lead for better tuning the management interventions.

Transcriptomic analysis to affirm the regulatory role of long non-coding RNA in horn cancer of Indian zebu cattle breed Kankrej (Bos indicus).

Long non-coding RNA (lncRNA) was previously considered as a non-functional transcript, which now established as part of regulatory elements of biological events such as chromosome structure, remodeling, and regulation of gene expression. The study presented here showed the role of lncRNA through differential expression analysis on cancer-related coding genes in horn squamous cell carcinoma of Indian zebu cattle. A total of 10,360 candidate lncRNAs were identified and further analyzed for its coding potential ability using three tools (CPC, CPAT, and PLEK) that provide 8862 common lncRNAs. Pfam analysis of these common lncRNAs gave 8612 potential candidates for lncRNA differential expression analysis. Differential expression analysis showed a total of 59 significantly differentially expressed genes and 19 lncRNAs. Pearson's correlation analysis was used to identify co-expressed mRNA-lncRNAs to established relation of the regulatory role of lncRNAs in horn cancer. We established a positive relation of seven upregulated (XLOC_000016, XLOC_002198, XLOC_002851, XLOC_ 007383, XLOC_010701, XLOC_010272, and XLOC_011517) and one downregulated (XLOC_011302) lncRNAs with eleven genes that are related to keratin family protein, keratin-associated protein family, cornifelin, corneodesmosin, serpin family protein, and metallothionein that have well-established role in squamous cell carcinoma through cellular communication, cell growth, cell invasion, and cell migration. These biological events were found to be related to the MAPK pathway of cell cycle regulation indicating the role of lncRNAs in manipulating cell cycle regulation during horn squamous cell carcinomas that will be useful in identifying molecular portraits related to the development of horn cancer.

Exploring the eukaryotic diversity in rumen of Indian camel (Camelus dromedarius) using 18S rRNA amplicon sequencing.

In addition to a wide variety of anaerobic and facultative anaerobic bacteria, camel rumen also harbors a diverse of eukaryotic organisms. In the present study, the eukaryotic communities of camel rumen were characterized using 18S rRNA amplicon sequencing. Metagenomic DNA was isolated from rumen samples of fourteen adult Bikaneri and Kachchhi breeds of camel fed different diets containing Jowar, Bajra, Maize, and Guar. Illumina sequencing generated 27,161,904 number of reads corresponding to 1543 total operational taxonomic units (OTUs). Taxonomic classification of community metagenome sequences from all the samples revealed the presence of 92 genera belonging to 16 different divisions, out of which Ciliophora (73%), Fungi (13%) and Streptophyta (9%) were found to be the most dominant. Notably, the abundance of Ciliophora was significantly higher in the case of Guar feed, while Fungi was significantly higher in the case of Maize feed, indicating the influence of cellulose and hemicellulose content of feedstuff on the composition of eukaryotes. The results suggest that the camel rumen eukaryotes are highly dynamic and depend on the type of diet given to the animal. Pearson's correlation analysis suggested the ciliate protozoa and fungi were negatively correlated with each other. To the best of our knowledge, this is first systematic study to characterize camel rumen eukaryotes, which has provided newer information regarding eukaryotic diversity patterns amongst camel fed on different diets.

MetaRNAseq analysis of surti buffalo rumen content reveals that transcriptionally active microorganisms need not be abundant.

The present study describes rumen microbiota composition and their functional profiles in Indian Surti buffaloes by metagenomic (MG) and metatranscriptomic (MT) approaches. The study compares samples from buffaloes fed three different proportion of roughages; green and dry type of roughage; and different rumen liquor fractions. Irrespective of sample, Bacteroidetes and Firmicutes were the most predominant bacterial phyla, followed by Proteobacteria, Fibrobacteres and Actinobacteria while, Prevotella, Bacteroides, Ruminococcus and Clostridium were the most abundant genera. Different proportions of taxa were observed in both MG and MT approaches indicating the differences in organisms present and organisms active in the rumen. Higher proportions of fungal taxa were observed in MT while important organisms like Fibrobacter and Butyrivibrio and abundant organisms like Bacteroides and Prevotella were underrepresented in MT data. Functionally, higher proportions of genes involved in Carbohydrate metabolism, Amino acid metabolism and Translation were observed in both data. Genes involved in Metabolism were observed to be underrepresented in MT data while, those involved in Genetic information processing were overrepresented in MT data. Further, genes involved in Carbohydrate metabolism were overexpressed compared to genes involved in Amino acid metabolism in MT data compared to MG data which had higher proportion of genes involved in Amino acid metabolism than Carbohydrate metabolism. In all significant differences were observed between both approaches, different fractions of rumen liquor (liquid and solid) and different proportions of roughage in diet.

The road less travelled: The efficacy of canine pluripotent stem cells.

The path-breaking discovery of induced pluripotent stem cells has fuelled the scientific advancements of stem cells. Nevertheless, the need to ensure the safety of stem cell therapy at translational level is still at large, prompting scientists to use animal models which are genetically and anatomically homologous to that of humans. Dogs, being genomically and physiologically more similar to humans serve as better models in mimicking human diseases as compared to rodents. The heterogeneity in canine breeds offers an excellent opportunity to comprehend the complexities of many genetic diseases, making them exceptional tools for stem cell therapies. Various canine gene therapy models have paved the foundation for strategizing therapies for humans. But a similar progress is lacking in utilizing canine stem cells for stem cell-based therapies in both dogs and humans. This review attempts to bridge the gap, by articulating the key differences in canine pluripotency pathways, based on the recent derivation of canine embryonic stem cells (cESCs) and canine induced pluripotent stem cells (ciPSCs), thereby attempting to position dog in the reprogramming landscape. The potential clinical application of canine iPSCs also offers great hope to canine patients and might lead to significant contributions in veterinary medicine.

Molecular characterization and expression profiling of 17-beta-hydroxysteroid dehydrogenase 2 and spermatogenesis associated protein 2 genes in endangered catfish, Clarias magur (Hamilton, 1822).

The 17-beta-hydroxysteroid dehydrogenase 2 (17ß-HSD2) enzyme regulates steroid levels by the inactivation of estrogen and androgens. Spermatogenesis associated protein 2 (SPATA2) plays a vital role in spermatogenesis in vertebrates including fish. We report cloning and characterization of full cds of 17ß-HSD2 and SPATA2 genes in Clarias magur. The full-length cDNA sequences of 17ß-HSD2 and SPATA2 were 1187 bp (ORF 1125 bp) and 1806 bp (ORF 1524 bp) encoding 375 and 508 amino acids, respectively. Signal peptide analysis revealed SPATA2 is nonsecretory, while 17ß-HSD2 is a secretory protein. Hydropathy profiles showed both proteins are hydrophilic in nature. Tissue distribution of both the genes revealed high mRNA level of SPATA2 in all tissues examined indicating its wide range of expression. 17ß-HSD2 indicated higher expression in preparatory phase compared to spawning phase in ovary while it was opposite in case of testis. SPATA2 showed significantly higher expression in preparatory phase compared to spawning phase in both ovary and testis. Administration of OvatideTM (GnRH analog) resulted in upregulation of SPATA2 expression at 6 and 16 h post-injection while 17ß-HSD2 showed upregulation only at 6 h post-injection. To the best of our knowledge, this is a first report on characterization of 17ß-HSD2 and SPATA2 full-length cDNA in catfish.

Metagenomic characterisation of ruminal bacterial diversity in buffaloes from birth to adulthood using 16S rRNA gene amplicon sequencing.

Microbial colonisation in the forestomach of a ruminant is one of the most crucial factors in determining many of its physiological developments and digestive capabilities. The present study attempts to identify establishment pattern of microbes in relation to food, age and rumen development in the buffalo calves at every fortnight interval from birth to 6 months of age, followed by every month till animals became 1 year of age. Diversity study based on 16S rRNA gene sequencing identified rapidly changing bacterial population during initial 60 days of life, which got assemblage as rumen became physiologically mature with increasing age of animals. A lactate fermenting aerobic to facultative anaerobic genera found during initial 30 days of life were expeditiously replaced by strict anaerobic cellulolytic bacterial population with increasing age. The study confirms that initial colonisation mainly depends on the oral cavity and skin of the mother, followed by the surrounding environment and feed offered, which is reversed in order once animal gets older. Some of the well-described genera based on culture-dependent studies like Ruminococcus spp. were found to be in lesser proportion suggesting an additional role of other microbes or niche in cellulose degradation. We report the presence of Porphyromonas spp. and Mannheimia glucosidal for the first time in bovine infants.

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.

Influence of Diet Composition on Cattle Rumen Methanogenesis: A Comparative Metagenomic Analysis in Indian and Exotic Cattle.

Comparative metagenomics approach has been used in this study to discriminate colonization of methanogenic population in different breeds of cattle. We compared two Indian cattle breeds (Gir and Kankrej) and two exotic cattle (Holstein and Jersey) breeds. Using a defined dietary plan for selected Indian varieties, the diet dependent shifts in microbial community and abundance of the enzymes associated with methanogenesis were studied. This data has been compared with the available rumen metagenome data from Holstein and Jersey dairy cattle. The abundance of genes for methanogenesis in Holstein and Jersey cattle came from Methanobacteriales order whereas, majority of the enzymes for methanogenesis in Gir and Kankrej cattle came from Methanomicrobiales order. The study suggested that by using slow/less digestible feed, the propionate levels could be controlled in rumen; and in turn, this would also help in further reducing the hydrogenotrophic production of methane. The study proposes that with the designed diet plan the overall methanogenic microbial pool or the individual methanogens could be targeted for development of methane mitigation strategies.

A genome-wide approach to screen for genetic variants in broilers (Gallus gallus) with divergent feed conversion ratio.

Feed conversion ratio (FCR) is an economically important trait in broilers and feed accounts for a significant proportion of the costs involved in broiler production. To explore the contribution of functional variants to FCR trait, we analyzed coding and non-coding single-nucleotide variants (SNVs) across the genome by exome sequencing in seven pairs of full-sibs broilers with divergent FCR and with a sequence coverage at an average depth of fourfold. We identified 192,119 high-quality SNVs, including 30,380 coding SNVs (cSNVs) in the experimental population. We discovered missense SNVs in PGM2, NOX4, TGFBR3, and TMX4, and synonymous SNVs in TSNAX, ITA, HSP90B1, and COL18A1 associated with FCR. Haplotype analyses of genome-wide significant SNVs in PGM2, PHKG1, DGKZ, and SOD2 were also observed with suggestive evidence of haplotype association with FCR. Single-variant and FCR QTL-related genes-based association analyses of SNVs identified newly associated genes for FCR in the regions subjected to targeted exome sequencing. The top seven SNVs were next evaluated in independent replication data sets where SNV chr. 3: 13,990,160 (c. 961G>C) at TMX4 was replicated (p < 0.05). Collectively, we have detected SNVs associated with FCR in broiler as well as identification of SNVs in known FCR QTL region. These findings should facilitate the discovery of causative variants for FCR and contribute to marker-assisted selection.