Integrative analysis of transcriptome and gwas data to identify the hub genes associated with milk yield trait in Buffalo

被引:0
作者
Deng T. [1 ,2 ]
Liang A. [1 ]
Liang S. [2 ]
Ma X. [2 ]
Lu X. [2 ]
Duan A. [2 ]
Pang P. [2 ]
Hua G. [1 ]
Liu S. [1 ]
Campanile G. [3 ]
Salzano A. [3 ]
Gasparrini B. [3 ]
Neglia G. [3 ]
Liang X. [2 ]
Yang L. [1 ]
机构
[1] Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, Huazhong Agricultural University, Wuhan
[2] Guangxi Provincial Key Laboratory of Buffalo Genetics, Breeding and Reproduction Technology, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning
[3] Department of Veterinary Medicine and Animal Productions, University of Naples Federico II, Naples
来源
Frontiers in Physics | 2019年 / 7卷 / FEB期
关键词
Buffalo; Genome-wide association studies; Hub genes; Milk yield; RNA-seq; WGCNA;
D O I
10.3389/fphy.2019.00036
中图分类号
学科分类号
摘要
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. © 2019 Baumgart, Schneider and Schütz.
引用
收藏
相关论文
共 66 条
[1]  
Akersr M.R., Bauman D.E., Capuco A.V., Goodman G.T., Tucker A.H., Prolactin regulation of milk secretion and biochemical differentiation of mammary epithelial cells in periparturient cows, Endocrinology, 109, pp. 23-30, (1981)
[2]  
Bionaz M., Loor J.J., Identification of reference genes for quantitative real-time PCR in the bovine mammary gland during the lactation cycle, Physiol. Genomics, 29, pp. 312-319, (2007)
[3]  
Bionaz M., Periasamy K., Rodriguez-Zas S.L., Everts R.E., Lewin H.A., Hurley W.L., Et al., Old and new stories: revelations from functional analysis of the bovine mammary transcriptome during the lactation cycle, PLoS One, 7, (2012)
[4]  
Bush W.S., Moore J.H., Chapter 11: genome-wide association studies, PLoS Comput. Biol., 8, (2012)
[5]  
Canovas A., Rincon G., Bevilacqua C., Islastrejo A., Brenaut P., Hovey R.C., Et al., Comparison of five different RNA sources to examine the lactating bovine mammary gland transcriptome using RNA-Sequencing, Sci. Rep., 4, (2014)
[6]  
Cao M.D., Cheng M., Rizwan A., Jiang L., Krishnamachary B., Bhujwalla Z.M., Et al., Targeting choline phospholipid metabolism: GDPD5 and GDPD6 silencing decrease breast cancer cell proliferation, migration, and invasion, NMR Biomed., 29, pp. 1098-1107, (2016)
[7]  
Cellurale C., Girnius N., Jiang F., Cavanagh-Kyros J., Lu S., Garlick D.S., Et al., Role of JNK in mammary gland development and breast cancer, Cancer Res., 72, pp. 472-481, (2012)
[8]  
Cui X., Hou Y., Yang S., Xie Y., Zhang S., Zhang Y., Et al., Transcriptional profiling of mammary gland in Holstein cows with extremely different milk protein and fat percentage using RNA sequencing, BMC Genomics, 15, (2014)
[9]  
Dai W.T., Zou Y.X., White R.R., Liu J.X., Liu H.Y., Transcriptomic profiles of the bovine mammary gland during lactation and the dry period, Funct. Integ. Genomics, 18, pp. 1-16, (2017)
[10]  
De Camargo G., Aspilcueta-Borquis R.R., Fortes M., Porto-Neto R., Cardoso D.F., Santos D., Et al., Prospecting major genes in dairy buffaloes, BMC Genomics, 16, (2015)