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.

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.

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.

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.

Identification of stably expressed Internal Control Genes (ICGs) for normalization of expression data in liver of C57BL/6 mice injected with beta casomorphins.

In recent years, beta-casomorphin peptides (BCM7/BCM9) derived from the digestion of cow milk have drawn a lot of attention world over because of their proposed impact on human health. In order to evaluate the transcriptional modulation of target genes through RT-qPCR in response to these peptides, availability of appropriate reference or internal control genes (ICGs) will be the key. The present study was planned to identify a panel of stable ICGs in the liver tissue of C57BL/6 mice injected with BCM7/BCM9 cow milk peptides for 3 weeks. A total of ten candidate genes were evaluated as potential ICGs by assessing their expression stability using software suites; geNorm, NormFinder and BestKeeper. The suitability of the identified ICGs was validated by assessing the relative expression levels of target genes, HP and Cu/Zn SOD. Based on geNorm, PPIA and SDHA gene pair was identified to be most stably expressed in liver tissue during the animal trials. Similarly, NormFinder analysis also identified PPIA as the most stable gene. BestKeeper analysis showed crossing point SD value for all the genes in the acceptable range that is closer to 1. Overall, the study identified a panel of stable ICGs for reliable normalization of target genes expression data in mice liver tissues during BCM7/9 peptides trial.

ASIP gene polymorphism associated with black coat and skin color in Murrah buffalo.

The melanogenesis pathway regulates pigmentation through the synergic action of various genes. We are interested in analyzing the genetic variations in the ASIP which determine eumelanin production in the dermis layer. In the present study, the ASIP gene was characterized in buffalo and 268 genetically unrelated buffaloes belonging to 10 different populations were genotyped for the non-synonymous SNP (c.292C>T) identified in the exon 3 region of the gene using Tetra-ARMS-PCR. The TT genotype occurred at a higher rate in Murrah, followed by Nili Ravi, Tripura, and Paralakhemundi (42.63%, 19.30%, 3.45%, and 3.33%). These results convey the association of the black coat color of Murrah with the ASIP gene TT genotype and the lighter shades of black coat (brown and grayish-black) color phenotype in other breeds with the CC genotype.

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.

Selection of species specific panel of reference genes in peripheral blood mononuclear cells of native livestock species adapted to trans-Himalayan region of Leh-Ladakh.

The identification of appropriate references genes is an integral component of any gene expression-based study for getting accuracy and reliability in data interpretation. In this study, we evaluated the expression stability of 10 candidate reference genes (GAPDH, RPL4, EEF1A1, RPS9, HPRT1, UXT, RPS23, B2M, RPS15, ACTB) in peripheral blood mononuclear cells of livestock species that are adapted to high altitude hypoxia conditions of Leh-Ladakh. A total of 37 PBMCs samples from six native livestock species of Leh-Ladakh region such as Ladakhi cattle, Ladakhi yak, Ladakhi donkey, Chanthangi goat, Double hump cattle and Zanskar ponies were included in this study. The commonly used statistical algorithms such as geNorm, Normfinder, BestKeeper and RefFinder were employed to assess the stability of these RGs in all the livestock species. Our study has identified different panel of reference genes in each species; for example, EEF1A1, RPL4 in Ladakhi cattle; GAPDH, RPS9, ACTB in Ladakhi yak; HPRT1, B2M, ACTB in Ladakhi donkey; HPRT1, B2M, ACTB in Double hump camel, RPS9, HPRT1 in Changthangi goat, HPRT1 and ACTB in Zanskar ponies. To the best of our knowledge, this is the first systematic attempt to identify panel of RGs across different livestock species types adapted to high altitude hypoxia conditions. In future, the findings of the present study would be quite helpful in conducting any transcriptional studies to understand the molecular basis of high altitude adaptation of native livestock population of Leh-Ladakh.

Identification of differential methylome signatures of white pigmented skin patches in Nili Ravi buffalo of India.

The DNA methylation events mark a major epigenetic change in the genome, reflecting non-genetic disease developments and varied phenotypes. The water buffalo is a dairy production animal with wide agro-climatic distribution in India. Breed-wise the coat color of water buffalo varies from ash-gray to jet black. A typical pigmentation pattern is found in one of the breeds of North India, Nili Ravi, with variedly distributed white patches. The DNA methylation pattern could potentially reveal the epigenetic factors responsible for the pigmentation patterns. To address this question, the DNA isolated from the skin tissues of Nili Ravi with varied white pigmentation and black Murrah buffaloes was subjected to reduced representation bisulfite sequencing. DNA methylation analysis revealed, 68.44%, 63.39%, and 47.94% of the promoter regions were hypermethylated in Nili Ravi over-white versus Murrah, Nili Ravi under-white versus Murrah, and Nili Ravi under-white versus Nili Ravi over-white, respectively. Major genes identified to be differentially methylated among over-white and under-white skin tissues in Nili Ravi included TBX2, SNAI2, HERC2, and CITED1. Overall the results have indicated differential methylation patterns to be potentially involved in hyper or hypopigmentation in Nili Ravi and Murrah buffaloes.

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.

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.

Morphometric and microsatellite-based comparative genetic diversity analysis in Bubalus bubalis from North India.

To understand the similarities and dissimilarities of a breed structure among different buffalo breeds of North India, it is essential to capture their morphometric variation, genetic diversity, and effective population size. In the present study, diversity among three important breeds, namely, Murrah, Nili-Ravi and Gojri were studied using a parallel approach of morphometric characterization and molecular diversity. Morphology was characterized using 13 biometric traits, and molecular diversity through a panel of 22 microsatellite DNA markers recommended by FAO, Advisory Group on Animal Genetic Diversity, for diversity studies in buffaloes. Canonical discriminate analysis of biometric traits revealed different clusters suggesting distinct genetic entities among the three studied populations. Analysis of molecular variance revealed 81.8% of genetic variance was found within breeds, while 18.2% of the genetic variation was found between breeds. Effective population sizes estimated based on linkage disequilibrium were 142, 75 and 556 in Gojri, Nili-Ravi and Murrah populations, respectively, indicated the presence of sufficient genetic variation and absence of intense selection among three breeds. The Bayesian approach of STRUCTURE analysis (at K = 3) assigned all populations into three clusters with a degree of genetic admixture in the Murrah and Nili-Ravi buffalo populations. Admixture analysis reveals introgression among Murrah and Nili-Ravi breeds while identified the Gojri as unique buffalo germplasm, indicating that there might be a common origin of Murrah and Nili-Ravi buffaloes. The study provides important insights on buffalo breeds of North India that could be utilized in designing an effective breeding strategy, with an appropriate choice of breeds for upgrading local non-descript buffaloes along with conservation of unique germplasm.

Identification of Internal Reference Genes in Peripheral Blood Mononuclear Cells of Cattle Populations Adapted to Hot Arid Normoxia and Cold Arid Hypoxia Environments.

To estimate gene expression in a reliable manner, quantitative real-time polymerase chain reaction data require normalisation using a panel of stably expressed reference genes (RGs). To date, information on an appropriate panel of RGs in cattle populations reared at cold arid high-altitude hypoxia and hot arid tropical normoxia environments is not available. Therefore, the present study was carried out to identify a panel of stably expressed RGs from 10 candidate genes (GAPDH, RPL4, EEF1A1, RPS9, HPRT1, UXT, HMBS, B2M, RPS15, and ACTB) in peripheral blood mononuclear cells (PBMCs) of cattle populations reared at cold arid high-altitude hypoxia and hot arid normoxia environments. Four different statistical algorithms: geNorm, NormFinder, BestKeeper, and RefFinder were used to assess the stability of these genes. A total of 30 blood samples were collected: six adult heifers each of Ladakhi (LAC) and Holstein Frisian crosses (HFX) and 4 Jersey (JYC) cows from cold arid high-altitude hypoxia environments (group I) and five adult heifers each of Sahiwal (SAC), Karan Fries (KFC), and Holstein Friesian (HFC) cows from hot arid normoxia environments (group II). Combined analysis of group I and group II resulted in identification of a panel of RGs like RPS9, RPS15, and GAPDH that could act as a useful resource to unravel the accurate transcriptional profile of PBMCs from diverse cattle populations adapted to distinct altitudes.

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.

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.

Transcriptome Analysis of Circulating PBMCs to Understand Mechanism of High Altitude Adaptation in Native Cattle of Ladakh Region.

Ladakhi cattle is native population of Leh and Ladakh region and constantly exposed to hypobaric hypoxia over many generations. In present study, transcriptome signatures of cattle from Ladakh region (~5500 m) and Sahiwal cattle from tropical regions were evaluated using Agilent 44 K microarray chip. The top up-regulated genes in Ladakhi cows were INHBC, ITPRI, HECA, ABI3, GPR171, and HIF-1α involved in hypoxia and stress response. In Sahiwal cows, the top up-regulated genes eEF1A1, GRO1, CXCL2, DEFB3 and BOLA-DQA3 were associated with immune function and inflammatory response indicating their strong immune potential to combat the pathogens prevalent in the tropical conditions. The molecular pathways highly impacted were MAPK signaling, ETC, apoptosis, TLR signaling and NF- kB signaling pathway indicating signatures of adaptive evolution of these two cattle types in response to diverse environments. Further, qPCR analysis revealed increased expression of DEGs such as HIF-1, EPAS-1, VEGFA, NOS2, and GLUT-1/SLC2A1 in cattle types from high altitude suggesting their pivotal role in association with high altitude adaptation. Based on data generated, native cattle of Ladakh region was found to be genetically distinct from native cattle adapted to the tropical region of India.

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.