Integrative Analysis of Transcriptome and GWAS Data to Identify the Hub Genes Associated With Milk Yield Trait in Buffalo

Abstract

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.

Integrative genomics of the mammalian alveolar macrophage response to intracellular mycobacteria

“…Integrative multi-omics approaches to data integration are now widely used to explore and dissect the genomic architecture and physiological basis of complex traits in domestic livestock, including network-based methods to integrate functional genomics and GWAS data (Canovas et al 2014;Fang et al 2017;Cai et al 2018;Fang et al 2018;Deng et al 2019;Yan et al 2020). However, to the best of our knowledge, this study is the first that uses network biology to systematically combine transcriptomics data from M. bovis-infected macrophages with GWAS data for M. bovis infection resistance in cattle.…”

Asian water buffalo: domestication, history and genetics

“…Comparative transcriptome analyses facilitate the discovery of differentially expressed genes and microRNAs. For example, Deng et al () performed an integrated analysis of GWAS and transcriptome data of lactating mammary gland tissues of Mediterranean river buffaloes at different lactation stages. Through a weighted gene co‐expression network analysis, they identified a set of 12 hub genes ( BNIPL , TUBA1C , C2CD4B , DCP1B , MAP3K5 , PDCD11 , SRGAP1 , GDPD5 , BARX2 , SCARA3 , CTU2 and RPL27A ) involved in several pathways related to milk production.…”

Using SNP Weights Derived From Gene Expression Modules to Improve GWAS Power for Feed Efficiency in Pigs

“…Genome-wide gene expression data from thousands of studies have been accumulating and made available through public repositories such as the Gene Expression Omnibus (GEO; Edgar et al, 2002). Recently, GWAS results have been interpreted by interrogating significant SNP for associations with gene expression data in livestock (Ballester et al, 2017;Fang et al, 2017;Kommadath et al, 2017;Cai et al, 2018;Deng et al, 2019). These studies have integrated GWAS and gene expression data post-GWAS.…”