Continent-wide genomic analysis of the African buffalo (Syncerus caffer).

The African buffalo (Syncerus caffer) is a wild bovid with a historical distribution across much of sub-Saharan Africa. Genomic analysis can provide insights into the evolutionary history of the species, and the key selective pressures shaping populations, including assessment of population level differentiation, population fragmentation, and population genetic structure. In this study we generated the highest quality de novo genome assembly (2.65 Gb, scaffold N50 69.17 Mb) of African buffalo to date, and sequenced a further 195 genomes from across the species distribution. Principal component and admixture analyses provided little support for the currently described four subspecies. Estimating Effective Migration Surfaces analysis suggested that geographical barriers have played a significant role in shaping gene flow and the population structure. Estimated effective population sizes indicated a substantial drop occurring in all populations 5-10,000 years ago, coinciding with the increase in human populations. Finally, signatures of selection were enriched for key genes associated with the immune response, suggesting infectious disease exert a substantial selective pressure upon the African buffalo. These findings have important implications for understanding bovid evolution, buffalo conservation and population management.

Mammary fat globules as a source of mRNA to model alterations in the expression of some milk component genes during lactation in bovines.

BACKGROUND: The milk's nutritional value is determined by its constituents, including fat, protein, carbohydrates, and minerals. The mammary gland's ability to produce milk is controlled by a complex network of genes. Thereby, the fat, protein, and lactose synthesis must be boost in milk to increase milk production efficiency. This can be accomplished by fusing genetic advancements with proper management practices. Therefore, this study aimed to investigate the association between the Lipoprotein lipase (LPL), kappa casein CSN3, and Glucose transporter 1 (GLUT1) genes expression levels and such milk components as fat, protein, and lactose in different dairy breeds during different stages of lactation. METHODS: To achieve such a purpose, 94 milk samples were collected (72 samples from 36 multiparous black-white and red-white Holstein-Friesian (HF) cows and 22 milk samples from 11 Egyptian buffaloes) during the early and peak lactation stages. The milk samples were utilized for milk analysis and genes expressions analyses using non- invasive approach in obtaining milk fat globules (MFGs) as a source of Ribonucleic acid (RNA). RESULTS: LPL and CSN3 genes expressions levels were found to be significantly higher in Egyptian buffalo than Holstein-Friesian (HF) cows as well as fat and protein percentages. On the other hand, GLUT1 gene expression level was shown to be significantly higher during peak lactation than early lactation. Moreover, lactose % showed a significant difference in peak lactation phase compared to early lactation phase. Also, fat and protein percentages were significantly higher in early lactation period than peak lactation period but lactose% showed the opposite pattern of Egyptian buffalo. CONCLUSION: Total RNA can be successfully obtained from MFGs. The results suggest that these genes play a role in glucose absorption and lactose synthesis in bovine mammary epithelial cells during lactation. Also, these results provide light on the differential expression of these genes among distinct Holstein-Friesian cow breeds and Egyptian buffalo subspecies throughout various lactation phases.

CRISPR-mediated editing of ß-lactoglobulin (BLG) gene in buffalo.

Milk is a good source of nutrition but is also a source of allergenic proteins such as α-lactalbumin, ß-lactoglobulin (BLG), casein, and immunoglobulins. The Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/Cas technology has the potential to edit any gene, including milk allergens. Previously, CRISPR/Cas has been successfully employed in dairy cows and goats, but buffaloes remain unexplored for any milk trait. In this study, we utilized the CRISPR/Cas9 system to edit the major milk allergen BLG gene in buffaloes. First, the editing efficiency of designed sgRNAs was tested in fibroblast cells using the T7E assay and Sanger sequencing. The most effective sgRNA was selected to generate clonal lines of BLG-edited cells. Analysis of 15 single-cell clones, through TA cloning and Sanger sequencing, revealed that 7 clones exhibited bi-allelic (-/-) heterozygous, bi-allelic (-/-) homozygous, and mono-allelic (-/+) disruptions in BLG. Bioinformatics prediction analysis confirmed that non-multiple-of-3 edited nucleotide cell clones have frame shifts and early truncation of BLG protein, while multiple-of-3 edited nucleotides resulted in slightly disoriented protein structures. Somatic cell nuclear transfer (SCNT) method was used to produce blastocyst-stage embryos that have similar developmental rates and quality with wild-type embryos. This study demonstrated the successful bi-allelic editing (-/-) of BLG in buffalo cells through CRISPR/Cas, followed by the production of BLG-edited blastocyst stage embryos using SCNT. With CRISPR and SCNT methods described herein, our long-term goal is to generate gene-edited buffaloes with BLG-free milk.

A Characterization and Functional Analysis of Peroxisome Proliferator-Activated Receptor Gamma Splicing Variants in the Buffalo Mammary Gland.

Peroxisome proliferator-activated receptor γ (PPARG) has various splicing variants and plays essential roles in the regulation of adipocyte differentiation and lipogenesis. However, little is known about the expression pattern and effect of the PPARG on milk fat synthesis in the buffalo mammary gland. In this study, we found that only PPARG-X17 and PPARG-X21 of the splicing variant were expressed in the buffalo mammary gland. Amino acid sequence characterization showed that the proteins encoded by PPARG-X17 and PPARG-X21 are endonuclear non-secreted hydrophilic proteins. Protein domain prediction found that only the PPARG-X21-encoded protein had PPAR ligand-binding domains (NR_LBD_PPAR), which may lead to functional differences between the two splices. RNA interference (RNAi) and the overexpression of PPARG-X17 and PPARG-X21 in buffalo mammary epithelial cells (BMECs) were performed. Results showed that the expression of fatty acid synthesis-related genes (ACACA, CD36, ACSL1, GPAT, AGPAT6, DGAT1) was significantly modified (p < 0.05) by the RNAi and overexpression of PPARG-X17 and PPARG-X21. All kinds of FAs detected in this study were significantly decreased (p < 0.05) after RNAi of PPARG-X17 or PPARG-X21. Overexpression of PPARG-X17 or PPARG-X21 significantly decreased (p < 0.05) the SFA content, while significantly increased (p < 0.05) the UFA, especially the MUFA in the BMECs. In conclusion, there are two PPARG splicing variants expressed in the BMECs that can regulate FA synthesis by altering the expression of diverse fatty acid synthesis-related genes. This study revealed the expression characteristics and functions of the PPARG gene in buffalo mammary glands and provided a reference for further understanding of fat synthesis in buffalo milk.

[Authenticity and adulteration identification of Bubali Cornu based on DNA fingerprints].

Whether adulteration exists is a difficult problem in the identification of traditional Chinese medicine(TCM). Bubali Cornu is mainly available in the medicinal material market in the form of buffalo horn silk or buffalo horn powder but lacks obvious identification characteristics, so there is a risk of adulteration. However, the method of identification of adulteration in Bubali Cornu is lacking at present. In order to ensure authenticity and identify adulteration of TCM Bubali Cornu, control the quality of TCM Bubali Cornu, and ensure the authenticity of clinical use, the DNA fingerprints of 43 batches of samples from pharmaceutical companies and medicinal material markets were identified, and the amplification primers of fluorescent DNA fingerprints of Bubali Cornu and Bovis Grunniens Cornu were screened. The DNA fingerprints of Bubali Cornu were obtained by fluorescent capillary typing. The identification effect of fluorescent capillary typing on different adulteration ratios was also tested. Two pairs of fluorescent STR typing primers, namely 16Sa and CRc, which can distinguish Bubali Cornu and Bovis Grunniens Cornu, were screened out, and a DNA fingerprint identification method was established. The 16Sa migration peaks of Bovis Grunniens Cornu and Bubali Cornu were 223.4-223.9 bp and 225.5-226.1 bp. The CRc migration peaks of Bovis Grunniens Cornu and Bubali Cornu were 518.8-524.8 bp and 535.9-542.5 bp. The peak height of the migration peak could be used for preliminary quantification of the adulterants with an adulteration ratio below 50%, and the quantitative results were similar to the adulteration ratio. In this study, a simple and quick universal DNA fingerprint method was established for the identification of Bubali Cornu and its adulterants, which could realize the identification of TCM Bubali Cornu and the semi-quantitative identification of the adulterants.

The role of lncFABP4 in modulating adipogenic differentiation in buffalo intramuscular preadipocytes.

Intramuscular fat (IMF) is a crucial determinant of meat quality and is influenced by various regulatory factors. Despite the growing recognition of the important role of long noncoding RNAs (lncRNAs) in IMF deposition, the mechanisms underlying buffalo IMF deposition remain poorly understood. In this study, we identified and characterized a lncRNA, lncFABP4, which is transcribed from the antisense strand of fatty acid-binding protein 4 (FABP4). lncFABP4 inhibited cell proliferation in buffalo intramuscular preadipocytes. Moreover, lncFABP4 significantly increased intramuscular preadipocyte differentiation, as indicated by an increase in the expression of the adipogenic markers peroxisome proliferator-activated receptor gamma (PPARG), CCAAT enhancer binding protein alpha (C/EBPα), and FABP4. Mechanistically, lncFABP4 was found to have the potential to regulate downstream gene expression by participating in protein-protein interaction pathways. These findings contribute to further understanding of the intricate mechanisms through which lncRNAs modulate intramuscular adipogenesis in buffaloes.

Regulatory impacts of PPARGC1A gene expression on milk production and cellular metabolism in buffalo mammary epithelial cells.

The PPARGC1A gene plays a fundamental role in regulating cellular energy metabolism, including adaptive thermogenesis, mitochondrial biogenesis, adipogenesis, gluconeogenesis, and glucose/fatty acid metabolism. In a previous study, our group investigated seven SNPs in Mediterranean buffalo associated with milk production traits, and the current study builds on this research by exploring the regulatory influences of the PPARGC1A gene in buffalo mammary epithelial cells (BuMECs). Our findings revealed that knockdown of PPARGC1A gene expression significantly affected the growth of BuMECs, including proliferation, cell cycle, and apoptosis. Additionally, we observed downregulated triglyceride secretion after PPARGC1A knockdown. Furthermore, the critical genes related to milk production, including the STATS, BAD, P53, SREBF1, and XDH genes were upregulated after RNAi, while the FABP3 gene, was downregulated. Moreover, Silencing the PPARGC1A gene led to a significant downregulation of ß-casein synthesis in BuMECs. Our study provides evidence of the importance of the PPARGC1A gene in regulating cell growth, lipid, and protein metabolism in the buffalo mammary gland. In light of our previous research, the current study underscores the potential of this gene for improving milk production efficiency and overall dairy productivity in buffalo populations.

Characterization and expression profiling of buffalo IFN-lambda family.

Interferon lambda (IFN-λ) is an important type III interferon triggered mainly by viral infection. IFN-λ binds to their heterodimeric receptors and signals through JAK-STAT pathways similar to type I IFN. In this study, we deduced the buffalo IFN-λ sequences through the polymerase chain reaction, and then studied IFN-λ's expression patterns in different tissues, and post induction with poly I:C and live MRSA using RT-qPCR. The full-length sequences of buffalo IFN-λ3, IFN-λ receptors, and a transcript variant of IFN-λ4 were determined. IFN-λ1 is identified as a pseudogene. Virus response elements and a recombination hotspot factor was observed in the regulatory region of IFN-λ. The IFN-λ3 expressed highest in lungs and monocytes but IFN-λ4 did not. The expression of Interferon Lambda Receptor 1 was tissue specific, while Interleukin 10 Receptor subunit beta was ubiquitous. Following poly I:C induction, IFN-λ3 expression was primarily observed in epithelial cells as opposed to fibroblasts, displaying cell type-dependent expression. The cytosolic RNA sensors were expressed highest in endometrial epithelial cells, whereas the endosomal receptor was higher in fibroblasts. 2',5'-oligoadenylate synthetase expressed higher in fibroblasts, myxoma resistance protein 1 and IFN-stimulated gene 56 in epithelial cells, displaying cell-specific antiviral response of the interferon stimulated genes (ISGs). The endometrial epithelial cells expressed IFN-λ3 after live S. aureus infection indicating its importance in bacterial infection. The induction of IFN-λ3 was S. aureus isolate specific at the same multiplicity of infection (MOI). This study elucidates the IFN-λ sequences, diverse expression patterns revealing tissue specificity, and specificity in response to poly I:C and bacterial stimuli, emphasising its crucial role in innate immune response modulation.

Double-Tube Multiplex TaqMan Real-Time PCR for the Detection of Eight Animal-Derived Dairy Ingredients.

Milk and dairy products represent important sources of nutrition in our daily lives. The identification of species within dairy products holds importance for monitoring food adulteration and ensuring traceability. This study presented a method that integrated double-tube and duplex real-time polymerase chain reaction (PCR) with multiplex TaqMan probes to enable the high-throughput detection of animal-derived ingredients in milk and dairy products. The detection system utilized one pair of universal primers, two pairs of specific primers, and eight animal-derived specific probes for cow, buffalo, goat, sheep, camel, yak, horse, and donkey. These components were optimized within a double-tube and four-probe PCR multiplex system. The developed double-tube detection system could simultaneously identify the above eight targets with a detection limit of 10-0.1 pg/µL. Validation using simulated adulterated milk samples demonstrated a detection limit of 0.1%. The primary advantage of this method lies in the simplification of the multiplex quantitative real-time PCR (qPCR) system through the use of universal primers. This method provides an efficient approach for detecting ingredients in dairy products, providing powerful technical support for market supervision.

Research Progress and Applications of Bovine Genome in the Tribe Bovini.

Various bovine species have been domesticated and bred for thousands of years, and they provide adequate animal-derived products, including meat, milk, and leather, to meet human requirements. Despite the review studies on economic traits in cattle, the genetic basis of traits has only been partially explained by phenotype and pedigree breeding methods, due to the complexity of genomic regulation during animal development and growth. With the advent of next-generation sequencing technology, genomics projects, such as the 1000 Bull Genomes Project, Functional Annotation of Animal Genomes project, and Bovine Pangenome Consortium, have advanced bovine genomic research. These large-scale genomics projects gave us a comprehensive concept, technology, and public resources. In this review, we summarize the genomics research progress of the main bovine species during the past decade, including cattle (Bos taurus), yak (Bos grunniens), water buffalo (Bubalus bubalis), zebu (Bos indicus), and gayal (Bos frontalis). We mainly discuss the development of genome sequencing and functional annotation, focusing on how genomic analysis reveals genetic variation and its impact on phenotypes in several bovine species.

Principal component analysis of phenotypic and breeding value data for semen traits in Egyptian buffalo bulls.

Buffalo bull semen traits are economically important traits that influence farm fertility and profitability. Genetic improvement of semen characteristics is an important detail of the genetic improvement. This study was conducted to assess the relationship between the breeding values as well as the phenotypic values for semen traits (VOL, MM, LS, AS and CONC) of the Egyptian buffalo bulls. A total of 7761 normal semen ejaculates were collected and characterized at ILMTC laboratory from 26 bulls from 2009 to 2019. For VOL, MM, LS, AS, and CONC, the actual means were 3.89 mL, 62.37%, 60.64%, 3.94%, and 0.67 × 109 sperm/mL, respectively. The prediction of breeding values for semen traits was estimated using a Bayesian procedure. The estimated standardized EBVs and phenotypic values were used in the principal component analysis (PCA). Of five PCs, one PC (PC1) had > 1 eigenvalues that was responsible for 87.19% of the total variation of SEBV, and two PCs had > 1 eigenvalues that were responsible for 59.61% and 21.35% of the total variation of the phenotypic values. Together, PC1 and PC2 accounted for 97.97% of the total variance of SEBV and 80.96% of the total variance of phenotypic values. A graphs of the first two components showed the traits separated into two different directions by group. This indicates each group was under similar genetic influence. Therefore, selection can be done separately for each group without influencing the other. Principal component analysis reduced variables to describe the key information in buffalo semen data.

Effect of heat exposure on prostaglandin production and expression of COX-2, PGES, PGFS, ITGAV and LGALS15 mRNAs in endometrial epithelial cells of buffalo (Bubalus bubalis).

BACKGROUND: Early embryonic mortality is one of the major intriguing factors of reproductive failure that causes considerable challenge to the mammalian cell biologists. Heat stress is the major factor responsible for reduced fertility in farm animals. The present study aimed to investigate the influence of heat stress on prostaglandin production and the expression of key genes, including COX-2, PGES, PGFS, ITGAV and LGALS15, in buffalo endometrial epithelial cells. METHODS AND RESULTS: Buffalo genitalia containing ovaries with corpus luteum (CL) were collected immediately post-slaughter. The stages of the estrous cycle were determined based on macroscopic observations of the ovaries. Uterine lumens of the mid-luteal phase (days 6-10 of the estrous cycle) were washed and treated with trypsin to isolate epithelial cells, which were then cultured at control temperature (38.5 °C for 24 h) or exposed to elevated temperatures [38.5 °C for 6 h, 40.5 °C for 18 h; Heat Stressed (HS)]. The supernatant and endometrial epithelial cells were collected at various time points (0, 3, 6, 12, and 24 h) from both the control and treatment groups. Although heat stress (40.5 °C) significantly (P < 0.05) increased COX-2, PGES, and PGFS transcripts in epithelial cells but it did not affect the in vitro production of PGF2α and PGE2. The expression of ITGAV and LGALS15 mRNAs in endometrial epithelial cells remained unaltered under elevated temperature conditions. CONCLUSION: It can be concluded that elevated temperature did not directly modulate prostaglandin production but, it promoted the expression of COX-2, PGES and PGFS mRNA in buffalo endometrial epithelial cells.

Weighted single step GWAS reveals genomic regions associated with economic traits in Murrah buffaloes.

The objective of the present study was to identify genomic regions influencing economic traits in Murrah buffaloes using weighted single step Genome Wide Association Analysis (WssGWAS). Data on 2000 animals, out of which 120 were genotyped using a double digest Restriction site Associated DNA (ddRAD) sequencing approach. The phenotypic data were collected from NDRI, India, on growth traits, viz., body weight at 6M (month), 12M, 18M and 24M, production traits like 305D (day) milk yield, lactation length (LL) and dry period (DP) and reproduction traits like age at first calving (AFC), calving interval (CI) and first service period (FSP). The biallelic genotypic data consisted of 49353 markers post-quality check. The heritability estimates were moderate to high, low to moderate, low for growth, production, reproduction traits, respectively. Important genomic regions explaining more than 0.5% of the total additive genetic variance explained by 30 adjacent SNPs were selected for further analysis of candidate genes. In this study, 105 genomic regions were associated with growth, 35 genomic regions with production and 42 window regions with reproduction traits. Different candidate genes were identified in these genomic regions, of which important are OSBPL8, NAP1L1 for growth, CNTNAP2 for production and ILDR2, TADA1 and POGK for reproduction traits.

Genome-wide diversity, population structure and signatures of inbreeding in the African buffalo in Mozambique.

The African buffalo, Syncerus caffer, is a key species in African ecosystems. Like other large herbivores, it plays a fundamental role in its habitat acting as an ecosystem engineer. Over the last few centuries, African buffalo populations have declined because of range contraction and demographic decline caused by direct or indirect human activities. In Mozambique, historically home to large buffalo herds, the combined effect of colonialism and subsequent civil wars has created a critical situation that urgently needs to be addressed. In this study, we focused on the analysis of genetic diversity of Syncerus caffer caffer populations from six areas of Mozambique. Using genome-wide SNPs obtained from ddRAD sequencing, we examined the population structure across the country, estimated gene flow between areas under conservation management, including national reserves, and assessed the inbreeding coefficients. Our results indicate that all studied populations of Syncerus caffer caffer are genetically depauperate, with a high level of inbreeding. Moreover, buffaloes in Mozambique present a significant population differentiation between southern and central areas. We found an unexpected genotype in the Gorongosa National Park, where buffaloes experienced a dramatic population size reduction, that shares a common ancestry with southern populations of Catuane and Namaacha. This could suggest the past occurrence of a connection between southern and central Mozambique and that the observed population structuring could reflect recent events of anthropogenic origin. All the populations analysed showed high levels of homozygosity, likely due to extensive inbreeding over the last few decades, which could have increased the frequency of recessive deleterious alleles. Improving the resilience of Syncerus caffer caffer in Mozambique is essential for preserving the ecosystem integrity. The most viable approach appears to be facilitating translocations and re-establishing connectivity between isolated herds. However, our results also highlight the importance of assessing intraspecific genetic diversity when considering interventions aimed at enhancing population viability such as selecting suitable source populations.

Doppler ultrasonographic scan, gene expression and serum profile of immune, APPs and antioxidant markers in Egyptian buffalo-cows with clinical endometritis.

The objective of this study was to elaborate Doppler ultrasonographic scan, genetic resistance and serum profile of markers associated with endometritis susceptibility in Egyptian buffalo-cows. The enrolled animals were designed as; twenty five apparently healthy buffalo-cows considered as a control group and twenty five infected buffalo with endometritis. There were significant (p < 0.05) increased of cervical diameter, endometrium thickness, uterine horn diameter, TAMEAN, TAMAX and blood flow through middle uterine artery with significant decrease of PI and RI values in endometritis buffalo-cows. Gene expression levels were considerably higher in endometritis-affected buffaloes than in resistant ones for the genes A2M, ADAMTS20, KCNT2, MAP3K4, MAPK14, FKBP5, FCAMR, TLR2, IRAK3, CCl2, EPHA4, and iNOS. The RXFP1, NDUFS5, TGF-ß, SOD3, CAT, and GPX genes were expressed at substantially lower levels in endometritis-affected buffaloes. The PCR-DNA sequence verdicts of healthy and affected buffaloes revealed differences in the SNPs in the amplified DNA bases related to endometritis for the investigated genes. However, MAP3K4 elicited a monomorphic pattern. There was a significant decrease of red blood cells (RBCs) count, Hb and packed cell volume (PCV) with neutrophilia, lymphocytosis and monocytosis in endometritis group compared with healthy ones. The serum levels of Hp, SAA, Cp, IL-6, IL-10, TNF-α, NO and MDA were significantly (P˂0.05) increased, along with reduction of CAT, GPx, SOD and TAC in buffalo-cows with endometritis compared to healthy ones. The variability of Doppler ultrasonographic scan and studied genes alongside alterations in the serum profile of investigated markers could be a reference guide for limiting buffalo endometritis through selective breeding of natural resistant animals.

Molecular and Structural Insights into Buffalo Interleukin-17A.

Interleukin-17A is a pro-inflammatory cytokine that plays a key role in the immune response to many pathogens and implicated in autoimmune diseases. This molecule is also involved in providing protection to many bacterial and fungal infections of gastro-intestinal tract and respiratory mucosa. Although molecular aspect of IL-17A has been studied in few species, no data are available for buffalo, which is one of the major sources of milk production in India. Therefore, in the present study, IL-17A gene of Indian Murrah Buffalo origin was cloned, expressed, and analyzed using bioinformatic tools. The coding sequence of buffalo IL-17A gene was cloned in prokaryotic expression vector (pET-28a) followed by its expression, purification, and characterization. A computational analysis was performed to understand the sequence, structure, and evolutionary relationship of buIL-17A. It revealed that the length of buIL-17A sequence without signal peptide is 132 amino acids as in cattle. However, sequence identity is found to be 99% due to one amino substitution difference between buffalo and cattle. After analysis, it can be concluded that buIL-17A recombinant protein can be used as a potential immunobiological reagent for diagnostic and therapeutic purpose.

Oocyte competence and gene expression in parthenogenetic produced embryos from repeat breeder and normally fertile buffaloes (Bubalus bubalis) raised in sub-humid tropical climate.

The reproductive management of the buffalo species still faces several unresolved problems, which directly affect the productivity of the herd, one of them being the presence of repeat breeder females. Given this scenario, this study aimed to verify the developmental competence of oocytes obtained from repeat breeder females and submitted to parthenogenetic activation. In addition, embryo gene expression was compared to normally fertile females. Murrah buffaloes were divided into two groups: repeat breeder (RB, n = 8) and normally fertile or control (CR, n = 7). Cumulus-oocyte complexes (COCs) were aspirated by transvaginal ovum pick-up from estrus synchronized females. The COCs were submitted to IVM for 24 h, and subsequently, the oocytes were activated using ionomycin, followed by 6-DMAP. Afterwards, the presumptive parthenotes were cultured for six or seven days in a microenvironment of 5 % CO2, 5 % O2, and 90 % N2 at 38.5 °C. The expression of OCT4, GLUT1, BCL2 and TFAM genes from blastocysts was evaluated. The overall COCs recovery rate was 70.9 % (190/268). The maturation (57.8 vs 71.1), cleavage (45.2 vs 62.2) and blastocyst (30.1 vs 45.9) rates did not differ (P > 0.05) between RB and CR females, respectively. Similarly, no significant difference (P > 0.05) was observed for the expression of studied genes in both RB and CR females. In conclusion, oocytes obtained from RB were as developmentally competent as those collected from CR females, with similar energy metabolism and in vitro development capacity. Thus, the low fertility rate of repeat breeder buffaloes, when compared to normal cyclic females, must be due to subsequent events to the blastocyst stage.

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.

Estimation of additive and maternal covariance of production traits in Murrah buffalo.

The study was done to determine additive, maternal and common permanent environmental effects and best-suited model for some production traits using six univariate animal models that differed in the (co)variance components fitted to assess the importance of maternal effect using likelihood ratio test in Murrah buffaloes. Data from 614 Murrah buffaloes related to production traits were collected from history pedigree sheets maintained at the buffalo farm, Department of Livestock Production and Management (LPM), LUVAS, Hisar. The production traits under this study were 305 days milk yield (305DMY), peak yield (PY), lactation length (LL), dry period (DP), lactation milk yield (LMY) and wet average (WA). The heritability estimates were in the range of 0.33-0.44 for 305DMY, 0.25-0.51 for PY, 0.05-0.13 for LL, 0.03-0.23 for DP, 0.17-0.40 for LMY and 0.37-0.66 for WA. Model 1 was considered best for most of the traits, viz., 305DMY, PY, LL, LMY and WA followed by model 2 for DP. Covariance and correlated values within the traits caused inflation of heritability in model 3 and model 6. The maximum covariance between the additive and maternal effect was found in trait LMY, which was 14,183.90 in model 3 and the minimum value was also reported in the same trait for model 6, valued at -3522.37. Multivariate analysis showed that all production traits were moderate to high and positively correlated with each other except for DP, which was low and negative genetic and phenotypic correlated. Spearman's rank correlation coefficients of breeding value among all six models were high and significant, ranged from 0.78 to 1.00 for all the traits except DP, therefore any of the models could be taken into account depending upon the availability of data.

Diversity assessment of a lesser known buffalo population from Central India and its comparative evaluation reveals presence of sufficient genetic variation and absence of selection.

Planned breeding and conservation strategies for a lesser-known population require an assessment of complete genetic diversity and population structure analysis in addition to its morphometric characteristics. In the present study, a comparative analysis of the genetic structure of a rare buffalo population, namely Chhattisgarhi, was extensively studied using a panel of FAO-recommended microsatellite markers along with well-established breeds namely Murrah, Nili-Ravi, Gojri, Kalahandi, and Nagpuri. Mode shift analysis indicated the absence of genetic bottleneck in the recent past. Assessment of genetic diversity indices across all loci indicated the presence of sufficient genetic variation within and between populations. Analysis of molecular variance between the six different buffalo populations attributed 19.05% of the variations to between-population differentiation. Cluster analyses using DAPC and Bayesian approach along with the phylogenetic tree based on UPGMA grouped six populations into three groups. The Chhattisgarhi population was revealed to be genetically closer to Nagpuri and Kalahandi populations. The study reveals the presence of sufficient genetic diversity within the Chhattisgarhi population and indicates the absence of a systematic selection program. We suggest improvement and conservation programs should be planned for this breed in the near future through short-term selection.