The mammary gland is the production organ in mammals that is of great importance for milk production and quality. However, characterization of the buffalo mammary gland transcriptome and identification of the valuable candidate genes that affect milk production is limited. Here, we performed the differential expressed genes (DEGs) analysis of mammary gland tissue on day 7, 50, 140, and 280 after calving and conducted gene-based genome-wide association studies (GWAS) of milk yield in 935 Mediterranean buffaloes. We then employed weighted gene co-expression network analysis (WGCNA) to identify specific modules and hub genes related to milk yield based on gene expression profiles and GWAS data. The results of the DEGs analysis showed that a total of 1,420 DEGs were detected across different lactation points. In the gene-based analysis, 976 genes were found to have genome-wide association (P ≤ 0.05) that could be defined as the nominally significant GWAS geneset (NSGG), 9 of which were suggestively associated with milk yield (P < 10−4). Using the WGCNA analysis, 544 and 225 genes associated with milk yield in the turquoise module were identified from DEGs and NSGG datasets, respectively. Several genes (including BNIPL, TUBA1C, C2CD4B, DCP1B, MAP3K5, PDCD11, SRGAP1, GDPD5, BARX2, SCARA3, CTU2, and RPL27A) were identified and considered as the hub genes because they were involved in multiple pathways related to milk production. Our findings provide an insight into the dynamic characterization of the buffalo mammary gland transcriptome, and these potential candidate genes may be valuable for future functional characterization of the buffalo mammary gland.

Collagens, as extracellular matrix proteins, support cells for structural integrity and contribute to support mammary basic structure and development. This study aims to perform the genomic identification, evolution, and expression analyses of the collagen gene family in water buffalo (Bubalus bubalis) during lactation. A total of 128 buffalo collagen protein sequences were deduced from the 45 collagen genes identified in silico from buffalo genome, which classified into six groups based on their phylogenetic relationships, conserved motifs, and gene structure analyses. The identified collagen sequences were unequally distributed on 16 chromosomes. The tandem duplicated genes were found within three chromosomes, while only one segmental event occurred between Chr3 and Chr8. Collinearity analysis revealed that a total of 36 collagen gene pairs were orthologous between buffalo and cattle genomes despite having different chromosome numbers. Comparative transcription analyses revealed that a total of 23 orthologous collagen genes were detected in the milk samples at different lactation periods between the two species. Notably, the duplicated gene pair of COL4A1-COL4A2 during lactation had a higher mRNA expression level than that of cattle, while a higher expression level of COL6A1-COL6A2 pair was found in cattle compared with that of buffalo. The present study provides useful information for investigating the potential functions of the collagen family in buffalo during lactation and helps in the functional characterization of collagen genes in additional research.

The water buffalo is an important dual-purpose livestock that is widespread throughout central and southern China. However, there has been no characterization of the population genetics of Chinese buffalo. Using an Axiom buffalo genotyping array (Thermo Fisher Scientific, Wilmington, DE), we analyzed the genetic diversity, linkage disequilibrium pattern, and signature of selection in 176 Chinese buffaloes from 13 breeds. A total of 35,547 SNP passed quality control and were used for further analyses. Population genetic analysis revealed a clear separation between swamp and river types. Ten Chinese indigenous breeds were clustered into the swamp group, the Murrah and Nili-Ravi breeds were clustered into the river group, and the crossbred breed was closer to the river group. Genetic diversity analysis showed that the swamp group had a lower average expected heterozygosity. Linkage disequilibrium decay distance was much shorter in the swamp group compared with the river group, with an average square of correlation coefficient value of 0.2 of approximately 50 kb. Analysis of runs of homozygosity indicated extensive remote and recent inbreeding within swamp and river groups, respectively. Moreover, one genomic region under selection was detected between the river and swamp groups. Our findings contribute to our understanding of the characterization of population genetics in Chinese buffaloes, which in turn may be used in buffalo breeding programs.

No abstract

Cytochrome P450 aromatase 19A1 (CYP19A1) is a critical enzyme in estrogen synthesis. However, the effect of CYP19A1 on cell growth and hormone secretion of buffalo follicular granulosa cells (BFGCs) is poorly understood. The objective of this study was to assess the role of CYP19A1 in cell proliferation and hormone secretion of BFGCs by knocking down CYP19A1 mRNA expression. The mRNA expression level of CYP19A1 gene was knocked down in BFGCs using the siCYP19A1-296 fragment with the best interference efficiency of 72.63%, as affirmed by real-time quantitative PCR (qPCR) and cell morphology analysis. The CYP19A1 knockdown promoted the proliferation of BFGCs through upregulating the mRNA expression levels of six proliferation-related genes (CCND1, CCNE1, CCNB1, CDK2, CDKN1A, and CDKN1B). Moreover, CYP19A1 knockdown increased (P < 0.05) the concentrations of progesterone secretion (P4) in BFGCs through increasing the mRNA expression levels of three steroidogenic genes (HSD17B1, HSD17B7, and CYP17A1). Our data further found that the FSH could inhibit the mRNA expression level of CYP19A1 in BFGCs, while LH obtains the opposite effect. These findings showed that the CYP19A1 knockdown had a regulatory role in the hormone secretion and cell proliferation in BFGCs.

12 Water buffalo holds the tremendous potential of milk and meat that widespread 13 throughout central and southern China. However, characterization of the population 14 genetics of Chinese buffalo is poorly understood. Using Axiom ® buffalo genotyping 15 array, we performed the genetic diversity, linkage disequilibrium (LD) pattern and 16 signature of selection in the 176 Chinese buffaloes from thirteen breeds. A total of 17 35,547 SNPs passed quality control and were used for further analyses. Population 18 genetic analysis revealed a clear separation between the swamp and river types. Ten 19 Chinese indigenous breeds clustered into the swamp group, Murrah and Nili-Ravi 20 breeds were the river group, and the crossbred breed was closer to the river group.21 Genetic diversity analysis showed that the swamp group had a lower average expected 22 heterozygosities compared to the river group. LD decay distance was much shorter in 23 the swamp group compared with the river group with value of approximately 50 2 0.2 24 Kb. Analysis of runs of homozygosity indicated that extensive remote and recent 25 inbreeding activity was respectively found within swamp and river groups. Moreover, 26 a total of 12 genomic regions under selection were detected between river and swamp 27 groups. Further, 12 QTL regions were found associated with buffalo milk production 28 traits. Some candidate genes within these QTLs were predicted to be involved in the 2 29 cell structure and function, suggesting that these genes might play vital roles in the 30 buffalo milk performance. Our data contribute to our understanding of the 31 characterization of population genetics in Chinese buffaloes, which in turn may be 32 utilized in buffalo breeding programs.33 Author Summary 34 Identifying the causal genes or markers associated with important economic traits in 35 livestock is critical to increasing the production level on the species. However, current 36 understanding of the genetic basis for milk production traits in buffalo is limited.37 Here, we confirmed the divergent evolution, distinct population structure, and LD 38 extent among Chinese buffalo breeds. We also identified 12 QTL regions associated 39 with milk production traits in buffaloes using the selective sweeps and haplotype 40 analysis. Further, a total of 7 genes involved in the cell structure and function were 41 predicted within the identified QTLs. These findings suggested that these genes can 42 serve as the candidate genes associated with buffalo milk production, which hold a 43 vital role in the milk trait improvement of dairy buffalo industry. Introduction45 Water buffalo (Bubalus bubalis), primarily raised for milk, meat and draught power, 46 is an essential part of the agricultural economy of many countries around the world. 47 Broadly, buffalo mainly consists of two types: the river (B. bubalis bubalis; n=50) and 48 swamp (B. bubalis carabensis; n=48), which are primarily distinguished based on 49 their distinct morphology and chromosomal karyotypes [1]. In China, the indi...

The organic anion transporter (OAT) family is the subfamily of the solute carrier (SLC) superfamily, which plays a vital role in regulating essential nutrients in milk. However, little is known about the members’ identification, evolutionary basis, and function characteristics of OAT genes associated with milk performance in buffalo. Comparative genomic analyses were performed to identify the potential role of buffalo OAT genes in milk performance in this study. The results showed that a total of 10 and 7 OAT genes were identified in river buffalo and swamp buffalo, respectively. These sequences clustered into three groups based on their phylogenetic relationship and had similar motif patterns and gene structures in the same groups. Moreover, the river-specific expansions and homologous loss of OAT genes occurred in the two buffalo subspecies during the evolutionary process. Notably, the duplicated SLCO3A1 gene specific to river buffalo showed higher expression level in mammary gland tissue than that of swamp buffalo. These findings highlight some promising candidate genes that could be potentially utilized to accelerate the genetic progress in buffalo breeding programs. However, the identified candidate genes require further validation in a larger cohort for use in the genomic selection of buffalo for milk production.

Acylglycerophosphate acyltransferases (AGPATs) are the rate-limiting enzymes for the de novo pathway of triacylglycerols (TAG) synthesis. Although AGPATs have been extensively explored by evolution, expression and functional studies, little is known on functional characterization of how many members of the AGPAT family are involved in TAG synthesis and their impact on the cell proliferation and apoptosis. Here, 13 AGPAT genes in buffalo were identified, of which 12 AGPAT gene pairs were orthologous between buffalo and cattle. Comparative transcriptomic analysis and real-time quantitative reverse transcription PCR (qRT-PCR) further showed that both AGPAT1 and AGPAT6 were highly expressed in milk samples of buffalo and cattle during lactation. Knockdown of AGPAT1 or AGPAT6 significantly decreased the TAG content of buffalo mammary epithelial cells (BuMECs) and bovine mammary epithelial cells (BoMECs) by regulating lipogenic gene expression (p < 0.05). Knockdown of AGPAT1 or AGPAT6 inhibited proliferation and apoptosis of BuMECs through the expression of marker genes associated with the proliferation and apoptosis (p < 0.05). Our data confirmed that both AGPAT1 and AGPAT6 could regulate TAG synthesis and growth of mammary epithelial cells in buffalo. These findings will have important implications for understanding the role of the AGPAT gene in buffalo milk performance.