Transcriptomic analysis reveals molecular insights into lactation dynamics in Jakhrana goat mammary gland.

BACKGROUND: Goat milk is gaining popularity as a superior alternative to bovine milk due to its closer resemblance to human milk. Understanding the molecular processes underlying lactation is crucial for improving milk quality and production in goats. However, the genetic mechanisms governing lactation in goats, particularly in indigenous breeds like the Jakhrana, remain largely unexplored. RESULTS: In this study, we performed a comprehensive transcriptomic analysis of Jakhrana goat mammary glands during early and late lactation stages. We isolated milk somatic cells and conducted RNA sequencing, followed by transcript quantification and mapping against the ARS1.2 Capra hircus reference assembly. Our analysis identified differentially expressed genes (DEGs) and commonly expressed genes (CEGs) across the lactation phases. Early lactation showed enrichment of genes encoding antimicrobial peptides and lubrication proteins, while late lactation exhibited heightened expression of genes encoding major milk proteins. Additionally, DEG analysis revealed upregulation of pivotal genes, such as the ABC transporter gene MRP4, implicated in modulating milk composition and quality. CONCLUSION: Our findings provide insights into the genetic mechanisms underlying lactation dynamics in the Jakhrana goat. Understanding these mechanisms could help in improving milk production and quality in goats, benefiting both the dairy industry and consumers.

Selection of reference genes for normalizing gene expression data across seasons in spermatozoa of water buffalo (Bubalus bubalis).

Selection of the most stably expressed reference genes is key to monitoring accurate target gene expression across any tissue or cell type. The mRNA in spermatozoa stores valuable information related to changes in spermatogenesis due to variations in environmental conditions, especially during heat stress, which affects various sperm functions. Semen quality in buffalo bulls is significantly influenced by the seasons. In the study, a panel of nine genes was evaluated to identify the most stably expressed internal control gene (ICG) for the normalization of real-time gene expression data generated across various seasons for Murrah buffalo bulls' spermatozoa. Sperm cells were purified from the semen samples collected during different seasons, with temperature-humidity index (THI) ranging from 80.80 ± 1.47 (hot summer) to 55.88 ± 1.98 (winter), using the BoviPure™ gradient purification method. The RNA isolated from the purified spermatozoa fraction was quality checked prior to reverse transcription and subjected to qPCR (quantitative real-time PCR) based expression analysis. An automated 'endoGene' pipeline was employed to apply the geNorm, NormFinder, and BestKeeper algorithms for data analysis. The result indicated that GAPDH and PP1A were the most stably expressed among the gene panel, whereas ATPSF1 and ACTB were the two least stable expressed reference genes. Further, the most suitable ICGs identified were validated by normalization of real time expression data of heat stress and sperm quality genes, HSFY2 and AKAP4, respectively. The genes identified would help in generating the most reliable results for the expression profiling of the genes dictating sperm quality and heat stress cope-up mechanism in buffalo spermatozoa, collected during different seasons.

Deciphering the molecular mechanisms of heat stress tolerance in goats: Insights from transcriptome and Gene Co-expression analysis.

Climate change poses a significant threat to the sustainability of livestock production systems in developing countries, particularly impacting small ruminants like goats, which are highly susceptible to heat stress. This stressor not only reduces productivity but also undermines economic viability. This study aimed to delve into the molecular mechanisms underlying heat stress tolerance in goats by conducting a comprehensive transcriptome analysis of heat-tolerant (HT, n = 4) and heat-susceptible (HS, n = 6) Jamunapari goats. Physiological metrics, such as rectal temperature, respiratory rate, and heart rate, were meticulously monitored under extreme environmental conditions (Temperature Humidity Index >92) to effectively classify goats based on their distinct heat stress responses. Samples of blood were obtained, and peripheral blood mononuclear cells (PBMCs) were extracted for subsequent RNA extraction. RNA-Seq analysis revealed a sum of 734 differentially expressed genes (DEGs), comprising 251 upregulated and 483 downregulated genes in HT goats compared to their HS counterparts. The WGCNA revealed three key modules, darkorange (tolerance), paleturquoise (respiration rate), and darkmagenta (heart rate). Moreover, functional enrichment analysis revealed that DEGs within these modules played intricate roles in crucial biological processes and pathways, including mitochondrial function, cardiac function, immune response, genomic stability, and metabolic regulation. This research notably enhances our comprehension of the genetic underpinnings of thermo-tolerance in goats and provides invaluable guidance for formulating breeding strategies aimed at bolstering livestock resilience against the challenges of climate change.

Polymorphism detection and characterization of sperm cells chromatin remodeling associated genes in Murrah buffalo.

Seasonal variations significantly impact buffalo bull semen production and quality, particularly during the summer months. Understanding the genetic basis of these changes is important for managing bull fertility and improving sperm quality. The present study focused on characterizing and identifying polymorphisms in chromatin remodeling genes, protamines (PRMs) and Transition Nuclear Proteins (TNPs) in Murrah buffalo bulls with varying semen quality due to seasonal effects. Our findings revealed none of the coding region variation in PRM1, PRM2, TNP1, and TNP2, these genes are highly conserved in buffalo. Two intronic variants were identified, including G16C in PRM1 intron 1 and intronic SNP in PRM2 intron 1 (G96A). The complete CDS of consensus sequence of bubaline PRM1 was 86.3% identical and 94.1% similar to the bovine PRM1. Whereas the complete CDS of consensus sequence of bubaline TNP2 was 78.2% identical and 91.0% similar to bovine TNP2. Further, no statistically significant differences in the fold change of TNP1, TNP2, PRM1, and PRM2 levels between the hot summer SNA and SA groups and the winter SNA and SA groups This study represents the first comprehensive report on the characterization of bubaline PRM1 (complete CDS), PRM2 (partial CDS), TNP1 (partial CDS), and TNP2 (complete CDS) genes in buffalo sperm cells. Results of the study, clearly indicate that the genes associated with protamine (PRM1 and TNP2) are highly conserved in Bubalus bubalis. Understanding these genetic underpinnings can have implications for improving buffalo bull fertility and semen quality.

Development and validation of most efficient RNA isolation method from buffalo bull spermatozoa.

BACKGROUND: Being highly fragmented and low in concentration, isolation of good quality RNA from sperm cells is a big challenge. Attempts have been made to evaluate various sperm RNA isolation methods from purified buffalo bull sperm cells. METHODS: Both, non-membrane and membrane-based methods have been evaluated for isolating RNA from Murrah buffalo sperms and compared for their respective efficacies. The traditional TRIzol, TRIzol-heat lysed (H-TRIzol) and cocktail of TCEP-RLT lysis buffer (Qiagen RNeasy mini kit)-TRIzol (C-TRIzol) based isopropanol isolation methods have been evaluated. RESULTS: H-TRIzol yielded best results among conventional methods. The combined T-RLT RNA isolation protocol yielded best quality and quantity compared to other membrane-based methods, due to high lytic property of cocktail of lysis reagents, necessary for complete breakdown of sperm membrane and RNA binding membrane for RNA isolation. Combined lysis performed by treatment with RLT-T and T-RLT differing in order of reagents used were also evaluated. T-RLT combination giving better results compared to RLT-T due to high gDNA contamination and membrane clogging in later protocol steps. CONCLUSION: Overall, in terms of total RNA quantity and quality per million spermatozoa, the heat-lysed TRIzol method (H-TRIzol) performs best among RNA separation techniques employed and is also quite easy to perform. This comparative evaluation of sperm RNA isolation protocols can be useful in deciding the best protocol for isolation of good quality and high concentration sperm RNA from buffalo semen, for transcriptome and other downstream studies.

Allelic diversity at BoLA DRB3 locus and association with predisposition to clinical mastitis in indicus and crossbred cattle.

Bovine lymphocyte antigen (BoLA) DRB3 locus in healthy and mastitis affected cattle has been genotyped by a polymerase chain reaction and restriction fragment length polymorphisms (PCR-RLFP) using RsaI restriction enzyme, followed by sequencing. In 130 farm animals, 25 BoLA DRB3 alleles have been detected by PCR-RFLP. Three distinct allelic patterns significantly associated with mastitis in Karan Fries crossbred and Sahiwal indicus cattle have been identified, whereas, four other allelic patterns were significantly high in frequency among healthy animals. Sequencing of RFLP genotypes revealed 25 and 47 alleles among healthy Sahiwal and Karan Fries, respectively, while 17 and 38 patterns observed in mastitis affected Sahiwal and Karan Fries animals, respectively. From Tajima's D-test of neutrality, it was concluded that alleles associated with mastitis were expanding in the population, whereas those of healthy were under contraction. Phylogenetic analysis carried out to delineate the evolutionary relationship of the farm and field animals at DRB3 locus, differentiating allelic patterns into six different clusters. Among the phylogenetic lineages, five patterns DRB3*028:01, DRB3*011:03, DRB3*031:01, DRB3*001:01 and DRB3*043:01, were previously reported, whereas one novel allelic variant was observed in indicus and crossbred cattle. This information will help in further exploring the association between BoLA-DRB3 genetic diversity and disease resistance in distinct cattle breeds, important in designing breeding strategies for increasing the distribution of favorable alleles.

Sequence characterization and comparative expression profile of buffalo WNT10B gene in adult and fetal tissues.

In this study, Wingless-type MMTV (mouse mammary tumor virus) integration site family member (WNT10B) gene was sequence characterized in the Indian water buffalo. Sequence analysis revealed an open reading frame of 1176 nucleotides in buffalo, encoding 391 amino acids long protein. Nineteen nucleotide variations were observed between cattle and buffalo resulting in six amino acid changes. Phylogenetic analysis showed the clustering of ruminant species together. Real-time expression analysis of WNT10B in tissues collected from different organs of fetal and adult buffalo, revealed, the gene being abundantly expressed in the rumen and liver of the fetus. The fetal ovary, heart, kidney, lung, testis and mammary gland showed moderate expression, while in adult tissues, expression was high in the ovary, testis, brain, kidney, small intestine and liver, whereas lower expression was observed in the adult rumen. Significant differences in WNT10B expression levels were found for the brain, small intestine, testes, kidney, heart, rumen, and ovary when adult and fetal tissues were compared. A moderate level of genetic variation was found between cattle and buffalo WNT10B and expression patterns in a variety of tissues in adult buffalo implies that in addition to possible roles in adipogenesis and hematopoiesis, the WNT10B gene might be playing a significant role in other regulatory pathways as well.

Genetic admixture and population structure analysis of Indian water buffaloes (Bubalus bubalis) using STR markers.

BACKGROUND: India has a vast riverine and swamp buffalo diversity adapted to various agro-ecological conditions. In the present study, genetic diversity data for 10 different buffalo populations of India, using 20 highly polymorphic microsatellite markers has been generated for the genetic diversity analysis. The buffalo populations of Eastern Odisha state, were the primary focus. METHODS AND RESULTS: The minimal spanning network based on Bruvo's distance, PCA (Principal Component Analysis) based on the Fst (Fixation Index) values, and genetic admixture analysis using both the STRUCTURE and 'snapclust' were performed. The analysis could identify the Manda population as distinct from other Odisha buffalo breeds as well as adjoining Chhattisgarhi buffalo breeds. The total observed number of alleles ranged between 143 (Manda) and 301 (Paralakhemundi) with an average of 204 alleles per breed. The Sambhalpuri buffalo population also clustered into two separate subpopulations, half of the unique sub-population located geographically south-wards, displayed no admixture with any of the adjacent buffalo populations. The Manda buffalo population has shown sufficient allelic richness and heterozygosity under random mating being practiced in the field conditions. CONCLUSIONS: The study has led to the identification of the Manda as a distinct buffalo population, and the germplasm has been registered as a new Indian buffalo breed. Whereas, the Sambhalpuri population requires elaborate analysis to confirm the existence of two distinct sub-populations.

Alternate PCR assays for screening of JH1 mutation associated with embryonic death in Jersey cattle.

Jersey haplotype (JH) 1, a stop-gain lethal mutation in the CWC15 gene, causes embryonic losses in Jersey cattle. Two PCR based assays using Amplification Refractory Mutation System (T-ARMS-PCR) and restriction fragment length polymorphism (PCR-RFLP) were developed for screening of the JH1 in cattle. During the screening, seven among 30 Indian Jersey bulls were identified as carriers of the mutant JH1 allele, the first time in the country. These PCR assays are economical, rapid and accurate; and can be used separately or in combination for screening and cross-validation of the JH1 carriers in Jersey cattle.

Diversity analysis at MHC class II DQA locus in buffalo (Bubalus bubalis) indicates extensive duplication and trans-species evolution.

Variation at MHC Class II-DQA locus in riverine and swamp buffaloes (Bubu) has been explored in this study. Through sequencing of buffalo DQA, 48 nucleotide variants identified from 17 individuals, reporting 42 novel alleles, including one pseudogene. Individual animal displayed two to seven variants, suggesting the presence of more than two Bubu-DQA loci, as an evidence of extensive duplication. dN values were found to be higher than dS values at peptide binding sites, separately for riverine and swamp buffaloes, indicating locus being under positive selection. Evolutionary analysis revealed numerous trans-species polymorphism with alleles from water buffalo assigned to at least three different loci (Bubu-DQA1, DQA2, DQA3). Alleles of both the sub-species intermixed within the cluster, showing convergent evolution of MHC alleles in bovines. The results thus suggest that both riverine and swamp buffaloes share con-current arrangement of DQA region, comparable to cattle in terms of copy number and population polymorphism.

Current status and unique attributes of Indian Chilika buffalo for adaptation to brackish water ecology.

Chilika buffalo is native to the Eastern coast of India and well adapted to the largest coastal brackish water lagoon of Asia, Chilika Lake. We present here a report on the Chilika buffalo breed emphasizing the conservational urgency based on unique biochemical and molecular evidence related to liver and kidney functions while comparing it with tropically adapted other water buffalo breeds (Bubalus bubalis) of India. It is found that the Chilika buffalo breed has a better ability to withstand a long dehydration period as evident from its better glomerular filtration and higher expression of the ion transport channel. Mitochondrial D-loop sequencing results have shown these buffaloes being closer to swamp-type buffaloes of Bangladesh and northeast India and represent a unique "hybrid zone" on the eastern coast of India. Conservation of such uniquely adapted germplasm is crucial owing to the current global trend, where the introduction of exotic breeds has negatively impact "sui-generis" germplasm and they require higher managerial resource consumption for maintaining higher productivity. Further, the introduction of unconventional fisheries activities has proved detrimental to the lagoon ecosystem, potentially causing more threat to the buffalo's population.

Evaluation of therapeutic potential of recombinant buffalo lactoferrin N-lobe expressed in E. coli.

Lactoferrin (Lf) is a multifunctional bi-lobate iron-binding glycoprotein belonging to transferrin family with a mass of approximately 80 kD. Being ubiquitously present in almost all biological secretions, it performs important biological functions. One of the earliest and very well-documented functions of Lf is the antibacterial effect against broad spectrum Gram-negative and Gram-positive bacteria. In this study, buffalo Lf N-lobe cDNA was amplified, cloned and expressed as a fusion protein in Escherichia coli cells using pQE30 expression vector. After post-induction confirmation of expressed protein by SDS-PAGE, purification of recombinant protein using Ni-NTA was attempted and the yield of recombinant buffalo N-lobe Lf was estimated to be 1 mg/ml. Antibacterial activity of recombinant buffalo Lf N-lobe was assessed on pathogenic E. coli and Staphylococcus aureus strains. Peptic digest of recombinant N-lobe buffalo Lf showed antibacterial activity comparable to commercially available bovine Lf. The successful expression and characterization of functional recombinant N-lobe of buffalo Lf expressed in E. coli opens new vistas for developing alternate therapeutics, particularly against the diseases caused by Gram-negative microbes such as septicemia and diarrhea in newborn calves and mastitis in dairy animals.

Matrix-based three-dimensional culture of buffalo mammary epithelial cells showed higher induction of genes related to milk protein and fatty acid metabolism.

Demanding transcriptomic studies in livestock animal species could be replaced by good in vitro models mimicking the function of mammary gland. Mammary epithelial cells (MEC) are the functional unit of the mammary gland. Extracellular matrix is known to be a key factor providing normal homeostasis in three-dimensional (3D) environment as important signals are lost when cells are cultured in two-dimensional (2D) environment. The aims of this study were to establish a buffalo mammary epithelial cells (BMECs) in 3D culture using extracellular matrix and to determine whether such a 3D culture model has different expression pattern than 2D counterpart. The purified MEC generated after several passages were used to establish 3D culture using Geltrex matrix. The expression of milk casein genes viz., alpha S1-casein (CSN1S1), alpha S2-casein (CSN1S2), beta-casein (CSN2), kappa-casein (CSN3); and fatty acid metabolism genes viz., butyrophilin (BTN1A1), glycerol-3-phosphate acyltransferase (GPAM), fatty acid-binding protein 3 (FABP3), and stearoyl-CoA desaturase (SCD) was assessed in 3D culture in comparison to traditional monolayer culture using qRT-PCR. Notable morphological differences were observed for BMECs grown in 3D culture in comparison to 2D culture. Morphologically, epithelial structures grown in Geltrex matrix (3D) environment showed enhanced functional differentiation in comparison to 2D culture. In 3D culture, lumen and dome-like structures were formed by day 5, whereas polarized acinus-like structure were formed within 15 days of culturing. The expression data showed higher mRNA induction of milk casein and fatty acid metabolism genes in 10-day-old 3D BMECs culture in comparison to 2D monolayer culture. The result suggests that 3D organization of epithelial cells has favorable effect on induction of milk and fatty acid metabolism-related genes. Therefore, matrix-based 3D culture of MEC that recapitulate the structural and functional context of normal tissues could provide a better in vitro model to understand the mammary gland functioning of buffaloes.

Identification of polymorphism in fatty acid binding protein 3 (FABP3) gene and its association with milk fat traits in riverine buffalo (Bubalus bubalis).

The fatty acid binding protein 3 (FABP3) gene, known to be associated with fat percentage of milk and meat in bovines, was screened among swamp and riverine buffaloes for polymorphism detection and further association with milk fat contents. An SNP g.307C > T was identified in the intron 2 (+53 exon 2) region of FABP3 gene of Indian buffaloes. The SNP identified was genotyped in 692 animals belonging to 15 riverine, swamp and hybrid (riverine × swamp) buffalo populations of diverse phenotypes and utilities, by PCR-RFLP. A marked contrast was observed between the C and T allele frequencies in three types of buffaloes. The frequency of C allele ranged from 0.67 to 0.96 in pure swamp buffalo populations, with the highest in Mizoram (0.96). Whereas the frequency of T allele was high across all the Indian riverine buffalo breeds, ranging from 0.57 to 0.96. None of the genotypes at FABP3 g.307C > T locus was found to have significant association with milk fat and other production traits in Mehsana dairy buffalo breed. Our study revealed marked differences in the allele frequencies between riverine and swamp buffaloes at FABP3 g.307C > T locus, without any significant association with different milk traits in riverine buffaloes.

Genome-wide 5'-C-phosphate-G-3' methylation patterns reveal the effect of heat stress on the altered semen quality in Bubalus bubalis.

Semen production and quality are closely correlated with different environmental factors in bovines, particularly for the buffalo (Bubalus bubalis) bulls reared under tropical and sub-tropical conditions. Factors including DNA methylation patterns, an intricate process in sperm cells, have an impact on the production of quality semen in buffalo bulls under abiotic stress conditions. The present study was conducted to identify DNA methylome signatures for semen quality in Murrah buffalo bulls, acclaimed as a major dairy breed globally, under summer heat stress. Based on semen quality parameters that significantly varied between the two groups over the seasons, the breeding bulls were classified into seasonally affected (SA = 6) and seasonally non-affected (SNA = 6) categories. DNA was isolated from purified sperm cells and sequenced using the RRBS (Reduced Representation Bisulfite Sequencing) technique for genome-wide methylome data generation. During the hot summer months, the physiological parameters such as scrotal surface temperature, rectal temperature, and respiration rate for both the SA and SNA bulls were significantly higher in the afternoon than in the morning. Whereas, the global CpG% of SA bulls was positively correlated with the afternoon's scrotal surface and rectal temperature. The RRBS results conveyed differentially methylated cytosines in the promoter region of the genes encoding the channels responsible for Ca2+ exchange, NPTN, Ca2+ activated chloride channels, ANO1, and a few structure-related units such as septins (SEPT4 and SEPT6), SPATA, etc. Additionally, the hypermethylated set of genes in SA was significantly enriched for pathways such as the FOXO signaling pathway and oocyte meiosis. The methylation patterns suggest promoter methylation in the genes regulating the sperm structure as well as surface transporters, which could contribute to the reduced semen quality in the Murrah buffalo bulls during the season-related heat stress.

Structural and functional annotation of the porcine immunome.

BACKGROUND: The domestic pig is known as an excellent model for human immunology and the two species share many pathogens. Susceptibility to infectious disease is one of the major constraints on swine performance, yet the structure and function of genes comprising the pig immunome are not well-characterized. The completion of the pig genome provides the opportunity to annotate the pig immunome, and compare and contrast pig and human immune systems. RESULTS: The Immune Response Annotation Group (IRAG) used computational curation and manual annotation of the swine genome assembly 10.2 (Sscrofa10.2) to refine the currently available automated annotation of 1,369 immunity-related genes through sequence-based comparison to genes in other species. Within these genes, we annotated 3,472 transcripts. Annotation provided evidence for gene expansions in several immune response families, and identified artiodactyl-specific expansions in the cathelicidin and type 1 Interferon families. We found gene duplications for 18 genes, including 13 immune response genes and five non-immune response genes discovered in the annotation process. Manual annotation provided evidence for many new alternative splice variants and 8 gene duplications. Over 1,100 transcripts without porcine sequence evidence were detected using cross-species annotation. We used a functional approach to discover and accurately annotate porcine immune response genes. A co-expression clustering analysis of transcriptomic data from selected experimental infections or immune stimulations of blood, macrophages or lymph nodes identified a large cluster of genes that exhibited a correlated positive response upon infection across multiple pathogens or immune stimuli. Interestingly, this gene cluster (cluster 4) is enriched for known general human immune response genes, yet contains many un-annotated porcine genes. A phylogenetic analysis of the encoded proteins of cluster 4 genes showed that 15% exhibited an accelerated evolution as compared to 4.1% across the entire genome. CONCLUSIONS: This extensive annotation dramatically extends the genome-based knowledge of the molecular genetics and structure of a major portion of the porcine immunome. Our complementary functional approach using co-expression during immune response has provided new putative immune response annotation for over 500 porcine genes. Our phylogenetic analysis of this core immunome cluster confirms rapid evolutionary change in this set of genes, and that, as in other species, such genes are important components of the pig's adaptation to pathogen challenge over evolutionary time. These comprehensive and integrated analyses increase the value of the porcine genome sequence and provide important tools for global analyses and data-mining of the porcine immune response.

Sequence and topological characterization of Toll-like receptor 8 gene of Indian riverine buffalo (Bubalus bubalis).

In this study, buffalo (Bubalus bubalis) Toll-like receptor 8 (TLR8) gene has been characterized by sequence analysis and detecting polymorphism. Complete ORF of buffalo TLR8 gene was amplified using the RNA isolated from spleen tissue, which was found to be 3,102 nucleotides long encoding a 1,033 amino acid protein. Buffalo TLR8 had 10 nucleotide changes as compared to other livestock species resulting in six unique amino acid changes, four of them lying within leucine-rich repeat (LRR) domains. As compared to cattle (Bos indicus and Bos taurus), out of fifteen cysteine residues, fourteen were conserved and Cys at position 521 was replaced by Arg. Nine of the LRR domains had no amino acid change as compared to cattle, whereas LRR-C-terminus had maximum, five amino acid changes. Sequence characterization of 12 riverine and swamp buffaloes revealed presence of four polymorphic nucleotides, two of them were non-synonymous, one synonymous and one site in 3'UTR. PCR-RFLP genotyping of non-synonymous SNP 2758A>G (ILeu920Val) in Toll-interleukin-1 receptor domain of 463 swamp and riverine buffaloes showed a higher frequency of allele A in swamp (95 %) as compared to riverine (9.84 %) buffaloes.