Fertility is a complex trait that has a major impact on the development of the buffalo industry. Genome-wide association study (GWAS) has increased the ability to detect genes influencing complex traits, and many important genes related to reproductive traits have been identified in ruminants. However, reproductive traits are influenced by many factors. The development of the follicle is one of the most important internal processes affecting fertility. Genes found by GWAS to be associated with follicular development may directly affect fertility. The present study combined GWAS and RNA-seq of follicular granulosa cells to identify important genes which may affect fertility in the buffalo. The 90 K Affymetrix Axiom Buffalo SNP Array was used to identify the SNPs, genomic regions, and genes that were associated with reproductive traits. A total of 40 suggestive loci (related to 28 genes) were identified to be associated with six reproductive traits (first, second and third calving age, calving interval, the number of services per conception and open days). Interestingly, the mRNA expressions of 25 of these genes were also observed in buffalo follicular granulosa cells. The IGFBP7 gene showed high level of expression during whole antral follicle growth. The knockdown of IGFBP7 in buffalo granulosa cells promoted cell apoptosis and hindered cell proliferation, and increased the production of progesterone and estradiol. Furthermore, a notable signal was detected at 2.3–2.7 Mb on the equivalent of bovine chromosome 5 associated with age at second calving, calving interval, and open days. The genes associated with buffalo reproductive traits in this study may have effect on fertility by regulating of follicular growth. These results may have important implications for improving buffalo breeding programs through application of genomic information.

Background: The 90K Axiom Buffalo SNP Array is expected to improve and speed up various genomic analyses for the buffalo (Bubalus bubalis). Genomic prediction is an effective approach in animal breeding to improve selection and reduce costs. As buffalo genome research is lagging behind that of the cow and production records are also limited, genomic prediction performance will be relatively poor. To improve the genomic prediction in buffalo, we introduced a new approach (pGBLUP) for genomic prediction of six buffalo milk traits by incorporating QTL information from the cattle milk traits in order to help improve the prediction performance for buffalo. Results: In simulations, the pGBLUP could outperform BayesR and the GBLUP if the prior biological information (i.e., the known causal loci) was appropriate; otherwise, it performed slightly worse than BayesR and equal to or better than the GBLUP. In real data, the heritability of the buffalo genomic region corresponding to the cattle milk trait QTLs was enriched (fold of enrichment > 1) in four buffalo milk traits (FY270, MY270, PY270, and PM) when the EBV was used as the response variable. The DEBV as the response variable yielded more reliable genomic predictions than the traditional EBV, as has been shown by previous research. The performance of the three approaches (GBLUP, BayesR, and pGBLUP) did not vary greatly in this study, probably due to the limited sample size, incomplete prior biological information, and less artificial selection in buffalo. Conclusions: To our knowledge, this study is the first to apply genomic prediction to buffalo by incorporating prior biological information. The genomic prediction of buffalo traits can be further improved with a larger sample size, higher-density SNP chips, and more precise prior biological information.

Background: The diacylglycerol acyltransferases (DGAT) are a vital group of enzymes in catalyzing triacylglycerol biosynthesis. DGAT genes like DGAT1 and DGAT2 , have been identified as two functional candidate genes affecting milk production traits, especially for fat content in milk. Buffalo milk is famous for its excellent quality, which is rich in fat and protein content. Therefore, this study aimed to characterize DGAT family genes in buffalo and to find candidate markers or DGAT genes influencing lactation performance. Results: We performed a genome-wide detection and identified eight DGAT genes in buffalo. All the DGAT genes classified into two distinct clades (DGAT1 and DGAT2 subfamily) based on their phylogenetic relationships and structural features. Chromosome localization displayed eight buffalo DGAT genes distributed on five chromosomes. Collinearity analysis revealed that the DGAT family genes were extensive homologous between buffalo and cattle. Afterward, we discovered genetic variants loci within the genomic regions that DGAT genes located in buffalo. Eight haplotype blocks were constructed and were associated with buffalo milk production traits. Single marker association analyses revealed four most significant single nucleotide polymorphisms (SNPs) mainly affecting milk protein percentage or milk fat yield in buffalo. Genes functional analysis indicated that these DGAT family genes could influence lactation performance in the mammal through regulating lipid metabolism. Conclusion: In the present study, we performed a comprehensive analysis for the DGAT family genes in buffalo, which including identification, structural characterization, phylogenetic classification, chromosomal distribution, collinearity analysis, association analysis, and functional analysis. These findings provide useful information for an in-depth study to determine the role of DGAT family gens play in the regulation of milk production and milk quality improvement in buffalo. Keywords: DGAT family, DGAT1, DGAT2 , Milk production traits, Buffalo genome