利用WGCNA鉴定花生主茎生长基因共表达模块  

Identification of gene co-expression modules of peanut main stem growth by WGCNA

作　　者：汪颖 高芳[1] 刘兆新[1] 赵继浩 赖华江 潘小怡 毕晨 李向东[1] 杨东清[1] WANG Ying;GAO Fang;LIU Zhao-Xin;ZHAO Ji-Hao;LAI Hua-Jiang;PAN Xiao-Yi;BI Chen;LI Xiang-Dong;YANG Dong-Qing(College of Agronomy,Shandong Agricultural University/State Key Laboratory of Crop Biology,Tai’an 271018,Shandong,China)

机构地区：[1]山东农业大学农学院/作物生物学国家重点实验室,山东泰安271018

出　　处：《作物学报》2021年第9期1639-1653,共15页Acta Agronomica Sinica

基　　金：国家重点研发计划项目“大田经济作物优质丰产的生理基础与调控”(2018YFD1000900);山东省重大科技创新工程项目(2018YFJH0601-3);山东省农业重大应用技术创新项目(SD2019ZZ11);山东省现代农业产业技术体系花生创新团队首席专家专项基金(SDAIT-04-01)资助。

摘　　要：以不同主茎高花生品种为材料,利用转录组测序技术分析茎秆转录组基因表达的异同,并结合加权基因共表达网络分析(WGCNA),深入挖掘与主茎生长相关基因,深入认识花生茎秆形态建成的分子机制。结果表明,矮秆型Df216与高秆型花育33号相比较共有5872个差异基因;Df216与中间型山花108相比较共有6662个差异基因,这些差异基因涉及细胞壁和次生细胞壁的生物起源及调控过程、苯丙烷生物合成及代谢过程、木质素生物合成过程、纤维素合酶活性等分子功能。WGCNA鉴定到5个与主茎高呈极显著相关的共表达模块。编码咖啡酰辅酶A-O-甲基转移酶、转录因子ATAF2、WAT1、GDSL脂肪酶等基因是模块内的核心基因。通过筛选权重值构建核心基因的局部网络发现,Grey60模块的核心基因ADRL3L与编码莽草酸香豆酯/奎酸酯3’-羟化酶、4-香豆酸辅酶A连接酶、羟基肉桂酰辅酶A莽草酸/奎尼酸羟基肉桂酰转移酶、以及快速碱化因子、锌指蛋白、类COBRA蛋白等基因有较高互作网络关系;Brown模块核心基因TZB0A2则与编码b-1,4-半乳糖基转移酶、果胶乙酰酯酶、类受体丝氨酸/苏氨酸蛋白激酶、富含亮氨酸重复序列的伸展蛋白等基因有较高互作网络关系。相关模块与核心基因的挖掘以及基因生物学功能和互作网络的解析有助于揭示花生主茎生长的遗传基础。This study was investigated the difference of transcriptome using three different peanut varieties with high main stem by RNA-seq.Transcriptomics combined with weighted gene co-expression network analysis(WGCNA)was used to explore the hub genes related to main stem growth and the molecular mechanisms of morphological formation of peanut stems.Results showed that 5872 differential expressed genes(DEGs)were detected in the Df216 and Huayu 33 comparation group,while 6662 DEGs were detected in the Df216 and Shanhua 108 comparation group.GO analysis suggested that these DEGs were mainly in-volved in molecular function and biological process,including the primary and secondary cell wall organization and biogenesis,phenylpropanoid biosynthetic and metabolic process,lignin biosynthetic process,and cellulose synthase activity,respectively.There were 33 modules were identified by WGCNA,among which five modules(Grey60,Cyan,Darkolivegreen,Brown,and Blue)were highly significant association with main stem height.According to the connectivity of genes in modules,caffeoyl-CoA O-methyltransferase,transcription factor ATAF2,WAT1(walls are thin1),and GDSL esterase/lipase were the hub genes,respec-tively.The results of hub gene networks by weighted values indicated that coumaroylquinate 3’-monooxygenase,4-coumarate-CoA ligase,shikimate O-hydroxycinnamoyltransferase,rapid alkalinization factor,COBRA-like protein,and zinc finger protein had high connections with ADRL3L in the Grey60 module,whileβ-1,4-galactosyltransferase,LRR receptor-like serine/threonine-protein kinase,pectin acetylesterase,leucine-rich repeat extensin-like protein had high connections with TZB0A2 in the Brown module.The identification of co-expression modules and their hub genes,and the analysis of gene function and gene networks of key genes will be helpful for revealing the genetic basis of the main height in peanut.

关 键 词：花生 主茎高 转录组 加权基因共表达网络 核心基因 

分 类 号：S565.2[农业科学—作物学]