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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Exome-wide comparative analyses revealed differentiating genomic regions for performance traits in Indian native buffaloes.

In India, 20 breeds of buffalo have been identified and registered, yet limited studies have been conducted to explore the performance potential of these breeds, especially in the Indian native breeds. This study is a maiden attempt to delineate the important variants and unique genes through exome sequencing for milk yield, milk composition, fertility, and adaptation traits in Indian local breeds of buffalo. In the present study, whole exome sequencing was performed on Chhattisgarhi (n = 3), Chilika (n = 4), Gojri (n = 3), and Murrah (n = 4) buffalo breeds and after stringent quality control, 4333, 6829, 4130, and 4854 InDels were revealed, respectively. Exome-wide FST along 100-kb sliding windows detected 27, 98, 38, and 35 outlier windows in Chhattisgarhi, Chilika, Gojri, and Murrah, respectively. The comparative exome analysis of InDels and subsequent gene ontology revealed unique breed specific genes for milk yield (CAMSAP3), milk composition (CLCN1, NUDT3), fertility (PTGER3) and adaptation (KCNA3, TH) traits. Study provides insight into mechanism of how these breeds have evolved under natural selection, the impact of these events on their respective genomes, and their importance in maintaining purity of these breeds for the traits under study. Additionally, this result will underwrite to the genetic acquaintance of these breeds for breeding application, and in understanding of evolution of these Indian local breeds.

Genetic parameters and genome-wide association studies for mozzarella and milk production traits, lactation length, and lactation persistency in Murrah buffaloes.

Genetic and genomic analyses of longitudinal traits related to milk production efficiency are paramount for optimizing water buffaloes breeding schemes. Therefore, this study aimed to (1) compare single-trait random regression models under a single-step genomic BLUP setting based on alternative covariance functions (i.e., Wood, Wilmink, and Ali and Schaeffer) to describe milk (MY), fat (FY), protein (PY), and mozzarella (MZY) yields, fat-to-protein ratio (FPR), somatic cell score (SCS), lactation length (LL), and lactation persistency (LP) in Murrah dairy buffaloes (Bubalus bubalis); (2) combine the best functions for each trait under a multiple-trait framework; (3) estimate time-dependent SNP effects for all the studied longitudinal traits; and (4) identify the most likely candidate genes associated with the traits. A total of 323,140 test-day records from the first lactation of 4,588 Murrah buffaloes were made available for the study. The model included the average curve of the population nested within herd-year-season of calving, systematic effects of number of milkings per day, and age at first calving as linear and quadratic covariates, and additive genetic, permanent environment, and residual as random effects. The Wood model had the best goodness of fit based on the deviance information criterion and posterior model probabilities for all traits. Moderate heritabilities were estimated over time for most traits (0.30 ± 0.02 for MY; 0.26 ± 0.03 for FY; 0.45 ± 0.04 for PY; 0.28 ± 0.05 for MZY; 0.13 ± 0.02 for FPR; and 0.15 ± 0.03 for SCS). The heritability estimates for LP ranged from 0.38 ± 0.02 to 0.65 ± 0.03 depending on the trait definition used. Similarly, heritabilities estimated for LL ranged from 0.10 ± 0.01 to 0.14 ± 0.03. The genetic correlation estimates across days in milk (DIM) for all traits ranged from -0.06 (186-215 DIM for MY-SCS) to 0.78 (66-95 DIM for PY-MZY). The SNP effects calculated for the random regression model coefficients were used to estimate the SNP effects throughout the lactation curve (from 5 to 305 d). Numerous relevant genomic regions and candidate genes were identified for all traits, confirming their polygenic nature. The candidate genes identified contribute to a better understanding of the genetic background of milk-related traits in Murrah buffaloes and reinforce the value of incorporating genomic information in their breeding programs.

Analysis of OCT4 and PGP9.5 gene expression in prenatal and postnatal buffalo (Bubalus bubalis) testes.

This study aimed to investigate and characterize the spermatogonial stem cells (SSCs) in buffaloes at different stages of development, including prenatal, neonatal, prepubertal, and adult testes. We sought a comprehensive understanding of these cells through a combination of histological, immunohistochemical, and ultrastructural analyses. Specifically, we examined changes in the expression of two potential SSC markers, OCT4 and PGP9.5, using immunohistochemistry. Additionally, we conducted a real-time quantitative polymerase chain reaction (RT-qPCR) to assess the relative gene expression of OCT4 and PGP9.5. The relative expression of the OCT4 gene was down-regulated in the adult testes compared to its expression during prepubertal and neonatal life. The relative expression of the PGP9.5 gene was up-regulated in the neonatal testes and down-regulated in the prepubertal and adult testes. The spermatogonia were round, oval-to-ellipsoidal cells lying over the basement membrane (BM) with a round-to-oval nucleus. Based on the immunoexpression of the putative SSC markers, OCT4 and PGP9.5, we concluded that the proportion of stem cells was highest during the neonatal stage, followed by the prepubertal and prenatal stages. This finding sheds light on the dynamics of spermatogonial stem cells in buffalo testes at different developmental stages, providing valuable insights into these cells' regulation and potential applications.

A synonymous single nucleotide polymorphism (g.36417726C > A) in the Lama2 gene influencing fat deposition is associated with post-partum anestrus interval in Murrah buffalo.

Postpartum absence of estrus exhibition known as postpartum anestrus interval (PPAI) for more than 90 days after calving is a concerning issue for dairy buffalo farmers' economy. The PPAI duration is influenced by both management practices and animal genetics. Investigating genetic markers associated with PPAI is crucial for incorporating them into marker-assisted selection programs. Towards this goal, our study focused on exploring potential genetic markers from early postpartum adipose tissue gene networks. We successfully identified 24 Single Nucleotide Polymorphisms (SNPs) within 9 candidate genes. In our initial analysis involving 100 buffaloes, we detected a significant association (P = 0.02267) between a specific synonymous SNP within the Lama2 gene (g.36417726C > A) and PPAI. This finding was subsequently validated (P = 0.02937) in a larger cohort of 415 buffaloes, where the SNP explained 1.36 % of the genetic variance. Intriguingly, buffaloes with the CC genotype of this SNP exhibited a PPAI that was 12.71 ± 3.21 days longer compared to buffaloes with AA and CA genotypes. To gain insight into the functional relevance of this SNP, a computational analysis was performed which indicated that the C allele of the SNP (g.36417726C > A) increased the stability of LAMA2 mRNA compared to the A allele. This computational prediction was corroborated by observing a significant increase (P = 0.01798) in Lama2 gene expression (greater than 8-fold) and higher fat percentage (P < 0.05) in adipose tissue of CC genotypes (48.78 ± 1.87 %) compared to AA genotypes (33.59 ± 4.5 %). Furthermore, we noted a significant (P < 0.05) upregulation of C/ebpß, Pparγ, Fasn, C/ebpα, and Pnpla2 genes, along with the downregulation of Bmp2 and Ptch1 in CC genotypes as opposed to AA genotypes. This observation suggests the involvement of the Pparγ-mediated pathway in both adipogenesis and lipolysis within CC genotypes. In summary, our comprehensive analysis involving association and functional validation underscores the potential of the SNP (g.36417726C > A) within the Lama2 gene as a promising genetic marker against extended PPAI in Murrah buffalo.

Identification of hub genes and their expression profiling for predicting buffalo (Bubalus bubalis) semen quality and fertility.

Sperm transcriptomics provide insights into subtle differences in sperm fertilization competence. For predicting the success of complex traits like male fertility, identification of hub genes involved in various sperm functions are essential. The bulls from the transcriptome profiled samples (n = 21), were grouped into good and poor progressive motility (PM), acrosome integrity (AI), functional membrane integrity (FMI) and fertility rate (FR) groups. The up-regulated genes identified in each group were 87, 470, 1715 and 36, respectively. Gene networks were constructed using up- and down-regulated genes from each group. The top clusters from the upregulated gene networks of the PM, AI, FMI and FR groups were involved in tyrosine kinase (FDR = 1.61E-11), apoptosis (FDR = 1.65E-8), translation (FDR = 2.2E-16) and ribosomal pathway (FDR = 1.98E-21), respectively. From the clusters, the hub genes were identified and validated in a fresh set of semen samples (n = 12) using RT-qPCR. Importantly, the genes (fold change) RPL36AL (14.99) in AI, EIF5A (54.32) in FMI, and RPLP0 (8.55) and RPS28 (13.42) in FR were significantly (p < 0.05) up-regulated. The study suggests that the expression levels of MAPK3 (PM), RPL36AL + RPS27A or RPL36AL + EXT2 (AI), RPL36AL or RPS27A (FMI) and RPS18 + RPS28 (FR) are potential markers for diagnosing the semen quality and fertility status of bulls which can be used for the breeding program.

TWAS revealed significant causal loci for milk production and its composition in Murrah buffaloes.

Milk yield is the most complex trait in dairy animals, and mapping all causal variants even with smallest effect sizes has been difficult with the genome-wide association study (GWAS) sample sizes available in geographical regions with small livestock holdings such as Indian sub-continent. However, Transcriptome-wide association studies (TWAS) could serve as an alternate for fine mapping of expression quantitative trait loci (eQTLs). This is a maiden attempt to identify milk production and its composition related genes using TWAS in Murrah buffaloes (Bubalus bubalis). TWAS was conducted on a test (N = 136) set of Murrah buffaloes genotyped through ddRAD sequencing. Their gene expression level was predicted using reference (N = 8) animals having both genotype and mammary epithelial cell (MEC) transcriptome information. Gene expression prediction was performed using Elastic-Net and Dirichlet Process Regression (DPR) model with fivefold cross-validation and without any cross-validation. DPR model without cross-validation predicted 80.92% of the total genes in the test group of Murrah buffaloes which was highest compared to other methods. TWAS in test individuals based on predicted gene expression, identified a significant association of one unique gene for Fat%, and two for SNF% at Bonferroni corrected threshold. The false discovery rates (FDR) corrected P-values of the top ten SNPs identified through GWAS were comparatively higher than TWAS. Gene ontology of TWAS-identified genes was performed to understand the function of these genes, it was revealed that milk production and composition genes were mainly involved in Relaxin, AMPK, and JAK-STAT signaling pathway, along with CCRI, and several key metabolic processes. The present study indicates that TWAS offers a lower false discovery rate and higher significant hits than GWAS for milk production and its composition traits. Hence, it is concluded that TWAS can be effectively used to identify genes and cis-SNPs in a population, which can be used for fabricating a low-density genomic chip for predicting milk production in Murrah buffaloes.

Follistatin (FST) is expressed in buffalo (Bubalus bubalis) ovarian follicles and promotes oocyte maturation and early embryonic development.

Follistatin (FST), a member of the transforming growth factor-ß (TGF-ß) superfamily, has been identified as an inhibitor of follicle-stimulating hormone. Previous studies showed that it plays an important role in animal reproduction. Therefore, this study aims to investigate its effect on the maturation of buffalo oocytes in vitro, and the underlying mechanism of FST affecting oocyte maturation was also explored in buffalo cumulus cells. Results showed that FST was enriched in the ovary and expressed at different stages of buffalo ovarian follicles as well as during oocyte maturation and early embryo development. The FST expression level was up-regulated in MII buffalo oocytes compared with the GV stage (p < .05). To study the effects of FST on buffalo oocytes' maturation and early embryonic development, we added the pcD3.1 skeleton vector and PCD3.1-EGFP-FST vector into the maturation fluid of buffalo oocytes, respectively. It was demonstrated that FST promoted the in vitro maturation rate of buffalo oocytes and the blastocyst rate of embryos cultured in vitro (p < .05). By interfering with FST expression, we discovered that FST in cumulus cells plays a crucial role in oocyte maturation. Interference with the FST expression during the buffalo oocyte maturation did not affect the first polar body rate of buffalo oocyte (p > .05). In contrast, the location of mitochondria in oocytes was abnormal, and the cumulus expansion area was reduced (p < .05). After parthenogenetic activation, the cleavage and blastocyst rates of the FST-interfered group were reduced (p < .05). Furthermore, RT-qPCR was performed to investigate further the underlying mechanism by which FST enhances oocyte maturation. We found that overexpression of FST could up-regulate the expression level of apoptosis suppressor gene Bcl-2 and TGF-ß/SMAD pathway-related genes TGF-ß, SMAD2, and SMAD3 (p < .05). In contrast, the expression levels of SMAD4 and pro-apoptotic gene BAX were significantly decreased (p < .05). The FST gene could affect buffalo oocyte maturation by regulating the oocyte mitochondria integrity, the cumulus expansion, cumulus cell apoptosis, and the expression levels of TGF-ß/SMAD pathway-related genes.

Role of Different Members of the AGPAT Gene Family in Milk Fat Synthesis in Bubalus bubalis.

During triacylglycerol synthesis, the acylglycerol-3-phosphate acyltransferase (AGPAT) family catalyzes the conversion of lysophosphatidic acid to phosphatidic acid and the acylation of sn-2 fatty acids. However, the catalytic activity of different AGPAT members is different. Therefore, this study aimed to investigate the mechanism through which different AGPATs affect the efficiency of TAG synthesis and fatty acid composition. The conservation of amino acid sequences and protein domains of the AGPAT family was analyzed, and the functions of AGPAT1, AGPAT3, and AGPAT4 genes in buffalo mammary epithelial cells (BMECs) were studied using RNA interference and gene overexpression. Prediction of the protein tertiary structure of the AGPAT family demonstrated that four conservative motifs (motif1, motif2, motif3, and motif6) formed a hydrophobic pocket in AGPAT proteins, except AGPAT6. According to cytological studies, AGPAT1, AGPAT3, and AGPAT4 were found to promote the synthesis and fatty acid compositions of triacylglycerol, especially UFA compositions of triacylglycerol, by regulating ACSL1, FASN, GPAM, DGAT2, and PPARG gene expression. This study provides new insights into the role of different AGPAT gene family members involved in TAG synthesis, and a reference for improving the fatty acid composition of milk.