Qi Bao1,2, Xiaolan Zhang1,2, Pengjia Bao1,2, Chunnian Liang1,2, Xian Guo1,2, Min Chu1,2, Ping Yan3,4. 1. Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, 730050, China. 2. Key Laboratory of Yak Breeding Engineering, Lanzhou, 730050, China. 3. Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, 730050, China. pingyanlz@163.com. 4. Key Laboratory of Yak Breeding Engineering, Lanzhou, 730050, China. pingyanlz@163.com.
Abstract
BACKGROUND: As a mammal living at the highest altitude in the world, the yak has strong adaptability to the harsh natural environment (such as low temperature, scarce food, especially low oxygen) of Qinghai-Tibet Plateau (QTP) after a long process of natural selection. OBJECTIVE: Here, we used Weighted Correlation Network Analysis (WGCNA), a systematic biology method, to identify hypoxic adaptation-related modules and hub genes. The research of the adaptability of yak against hypoxia is of great significance to identify the genetic characteristics and yak breeding. METHODS: Based on the transcriptome sequencing data (PRJNA362606), the R package DESeq2 and WGCNA were conducted to analyze differentially expressed genes (DEGs) and construct the gene co-expression network. The module hub genes were identified and characterized by the correlation of gene and trait, module membership (kME). In addition, GO and KEGG enrichment analyses were used to explore the functions of hub genes. RESULTS: Our results revealed that 1098, 1429, and 1645 DEGs were identified in muscle, spleen, and lung, respectively. Besides, a total of 13 gene co-expression modules were detected, of which two hypoxic adaptation-related modules (saddlebrown and turquoise) were found. We identified 39 and 150 hub genes in these two modules. Functional enrichment analyses showed that 12 GO terms and 18 KEGG pathways were enriched in the saddlebrown module while 85 GO terms and 22 KEGG pathways were enriched in the turquoise module. The significant pathways related to hypoxia adaptation include FoxO signaling pathway, Thermogenesis pathway, and Retrograde endocannabinoid signaling pathway, etc. CONCLUSIONS: In this study, we obtained two hypoxia-related specific modules and identified hub genes based on the connectivity by constructing a weighted gene co-expression network. Function enrichment analysis of two modules revealed mitochondrion is the most important organelle for hypoxia adaptation. Moreover, the insulin-related pathways and thermogenic-related pathways played a major role. The results of this study provide theoretical guidance for further understanding the molecular mechanism of yak adaptation to hypoxia.
BACKGROUND: As a mammal living at the highest altitude in the world, the yak has strong adaptability to the harsh natural environment (such as low temperature, scarce food, especially low oxygen) of Qinghai-Tibet Plateau (QTP) after a long process of natural selection. OBJECTIVE: Here, we used Weighted Correlation Network Analysis (WGCNA), a systematic biology method, to identify hypoxic adaptation-related modules and hub genes. The research of the adaptability of yak against hypoxia is of great significance to identify the genetic characteristics and yak breeding. METHODS: Based on the transcriptome sequencing data (PRJNA362606), the R package DESeq2 and WGCNA were conducted to analyze differentially expressed genes (DEGs) and construct the gene co-expression network. The module hub genes were identified and characterized by the correlation of gene and trait, module membership (kME). In addition, GO and KEGG enrichment analyses were used to explore the functions of hub genes. RESULTS: Our results revealed that 1098, 1429, and 1645 DEGs were identified in muscle, spleen, and lung, respectively. Besides, a total of 13 gene co-expression modules were detected, of which two hypoxic adaptation-related modules (saddlebrown and turquoise) were found. We identified 39 and 150 hub genes in these two modules. Functional enrichment analyses showed that 12 GO terms and 18 KEGG pathways were enriched in the saddlebrown module while 85 GO terms and 22 KEGG pathways were enriched in the turquoise module. The significant pathways related to hypoxia adaptation include FoxO signaling pathway, Thermogenesis pathway, and Retrograde endocannabinoid signaling pathway, etc. CONCLUSIONS: In this study, we obtained two hypoxia-related specific modules and identified hub genes based on the connectivity by constructing a weighted gene co-expression network. Function enrichment analysis of two modules revealed mitochondrion is the most important organelle for hypoxia adaptation. Moreover, the insulin-related pathways and thermogenic-related pathways played a major role. The results of this study provide theoretical guidance for further understanding the molecular mechanism of yak adaptation to hypoxia.
Authors: Jamie R Crait; Henry D Prange; Noah A Marshall; Henry J Harlow; Clark J Cotton; Merav Ben-David Journal: J Exp Biol Date: 2012-01-15 Impact factor: 3.312
Authors: Jonathan P DiNitto; Anna Delprato; Meng-Tse Gabe Lee; Thomas C Cronin; Shaohui Huang; Adilson Guilherme; Michael P Czech; David G Lambright Journal: Mol Cell Date: 2007-11-30 Impact factor: 17.970