基于转录组及WGCNA的甘蔗干旱响应调控网络分析

doi:10.3724/SP.J.1006.2022.14121

作物学报 ›› 2022, Vol. 48 ›› Issue (7): 1583-1600.doi: 10.3724/SP.J.1006.2022.14121

• 作物遗传育种·种质资源·分子遗传学 • 上一篇下一篇

基于转录组及WGCNA的甘蔗干旱响应调控网络分析

李佩婷¹(), 赵振丽¹, 黄潮华¹, 黄国强¹, 徐良年¹, 邓祖湖¹, 张玉²^,^*(), 赵新旺¹^,³^,⁴^,^*()

¹福建农林大学国家甘蔗工程技术研究中心, 福建福州 350002
²福建农林大学生命科学学院 / 福建省植物功能生物学与绿色农业重点实验室, 福建福州 350002
³广西大学广西甘蔗生物学重点实验室, 广西南宁 530004
⁴福建农林大学 / 农业农村部福建甘蔗生物学与遗传育种重点实验室, 福建福州 350002

收稿日期:2021-07-13 接受日期:2021-10-19 出版日期:2022-07-12 网络出版日期:2021-11-25
通讯作者: 张玉,赵新旺
作者简介:E-mail: lipeiting147258@163.com
基金资助:
国家自然科学基金青年基金项目(31801423);国家自然科学基金青年基金项目(32001437)

Analysis of drought responsive regulatory network in sugarcane based on transcriptome and WGCNA

LI Pei-Ting¹(), ZHAO Zhen-Li¹, HUANG Chao-Hua¹, HUANG Guo-Qiang¹, XU Liang-Nian¹, DENG Zu-Hu¹, ZHANG Yu²^,^*(), ZHAO Xin-Wang¹^,³^,⁴^,^*()

¹National Engineering Research Center for Sugarcane, Fujian Agriculture and Forestry University, Fuzhou 350002, Fujian, China
²Fujian Key Laboratory of Plant Functional Biology and Green Agriculture / College of Life Sciences, Fujian Agriculture and Forestry University, Fuzhou 350002, Fujian, China
³Guangxi Key Laboratory of Sugarcane Biology, Guangxi University, Nanning 530004, Guangxi, China
⁴Fujian Key Laboratory of Sugarcane Biology, Genetics and Breeding, Ministry of Agriculture and Rural Affairs / Fujian Agriculture and Forestry University, Fuzhou 350002, Fujian, China

Received:2021-07-13 Accepted:2021-10-19 Published:2022-07-12 Published online:2021-11-25
Contact: ZHANG Yu,ZHAO Xin-Wang
Supported by:
Youth Fund Project of the National Natural Science Foundation of China(31801423);Youth Fund Project of the National Natural Science Foundation of China(32001437)

摘要/Abstract

摘要：

干旱是制约甘蔗产业发展的重要因素之一。前人研究表明, 斑茅具有良好的抗性基因, 可以通过远缘杂交遗传给后代。本研究以斑茅与甘蔗杂交F₁代无性系YCE96-40为材料, 对苗期干旱处理0 h和24 h后的叶和根进行转录组测序分析, 比较了根和叶在转录水平上响应干旱的差异, 鉴定出21,885个(DR vs CR:10176, DL vs CL:7907)差异表达基因(DEGs), 根中差异表达基因多于叶中, 说明根对干旱胁迫响应更为剧烈。GO功能富集分析发现, 根和叶中DEGs均富集到与脱水反应相关及激素信号转导过程相关的条目, 比如“对渗透胁迫生物过程的响应”和“对缺水生物过程的反应”等。与叶不同的是, 根中大量DEGs显著富集到与细胞膜相关的条目。在根中鉴定出多个木质素相关DEGs, 表明木质素参与了根的干旱响应。通过对所有DEGs进行WGCNA分析, 共鉴定出11个基因共表达模块, 其中与干旱处理后的根显著相关有5个模块, 与干旱处理后的叶显著相关的有2个模块。进一步筛选出了26个转录因子为甘蔗干旱响应的候选转录因子, 构建了调控网络。研究结果为进一步理解甘蔗抗旱性的分子机制及甘蔗抗旱性育种提供了理论指导。

关键词: 甘蔗, 抗旱性, 转录组, 加权基因共表达网络, 转录因子

Abstract:

Drought is one of the important factors restricting the development of sugarcane industry. Previous studies have shown that E. arundinaceus has good resistance genes, which can be inherited through distant hybridization. In this study, we used YCE96-40, a F₁ hybrid clone of E. arundinaceus and sugarcane, as material to analyze the transcriptome of leaves and roots after drought treatment for 0 and 24 hours at seedling stage. We compared the differences in response of roots and leaves to drought at transcriptional level, and identified 21,885 differentially expressed genes (DEGs) (DR vs CR: 10176, DL vs CL: 7907). The root response to drought stress was more intense. GO functional enrichment analysis showed that DEGs in roots and leaves were enriched in items related to dehydration response and hormone signal transduction process, such as “response to osmotic stress biological process” and “response to water shortage biological process”. Different from leaves, a large number of DEGs in roots were significantly enriched in the items related to cell membrane. Several lignin related DEGs were identified in roots, which indicated that lignin was involved in root drought response. Through WGCNA analysis of all DEGs, a total of 11 gene co-expression modules were identified, including five modules significantly related to root after drought treatment and two modules significantly related to leaf after drought treatment. Furthermore, 26 transcription factors were selected as candidate transcription factors for drought response in sugarcane, and the regulatory network was constructed. The results provide theoretical guidance for further understanding the molecular mechanism of sugarcane drought resistance and sugarcane drought resistance breeding.

Key words: sugarcane, drought resistance, transcriptome, weighted gene co-expression network, transcription factor

李佩婷, 赵振丽, 黄潮华, 黄国强, 徐良年, 邓祖湖, 张玉, 赵新旺. 基于转录组及WGCNA的甘蔗干旱响应调控网络分析[J]. 作物学报, 2022, 48(7): 1583-1600.

LI Pei-Ting, ZHAO Zhen-Li, HUANG Chao-Hua, HUANG Guo-Qiang, XU Liang-Nian, DENG Zu-Hu, ZHANG Yu, ZHAO Xin-Wang. Analysis of drought responsive regulatory network in sugarcane based on transcriptome and WGCNA[J]. Acta Agronomica Sinica, 2022, 48(7): 1583-1600.

图/表 22

表1

附表1

图1

图2

附表2

根中显著富集的GO条目"

GO条目GO_trem	描述 Description	条目类型Term_type	p值矫正值 FDR	DEGs数量 DEGs number
GO:0006804	过氧化物酶反应Obsolete peroxidase reaction	BP	4.02E-34	185
GO:0044710	单生物代谢过程Single-organism metabolic process	BP	2.95E-32	4319
GO:0006979	对氧化应激的反应Response to oxidative stress	BP	4.4E-32	271
GO:0009628	对非生物刺激的反应Response to abiotic stimulus	BP	4.4E-32	558
GO:1901700	对含氧化合物的反应Response to oxygen-containing compound	BP	4.4E-32	396
GO:0055114	氧化还原过程Oxidation-reduction process	BP	1.75E-31	1687
GO:0001101	对酸性化学物质的反应Response to acid chemical	BP	1.96E-31	319
GO:0044699	单生物过程Single-organism process	BP	3.51E-27	6670
GO:0042221	对化学物质的反应Response to chemical	BP	1.57E-22	812
GO:0009415	对水的反应Response to water	BP	3.42E-21	138
GO:0009737	脱落酸反应Response to abscisic acid	BP	5.87E-19	177
GO:0006970	渗透胁迫反应Response to osmotic stress	BP	1.1E-18	217
GO:0009719	对内源性刺激的反应Response to endogenous stimulus	BP	1.77E-18	359
GO:0009414	对缺水的反应Response to water deprivation	BP	1.98E-18	116
GO:0005975	碳水化合物代谢过程Carbohydrate metabolic process	BP	3.04E-17	972
GO:0009651	盐胁迫反应Response to salt stress	BP	4.57E-17	201
GO:0009725	对激素的反应Response to hormone	BP	1.54E-16	340
GO:0033993	对脂质的反应Response to lipid	BP	7.82E-16	202
GO:0010033	对有机物的反应Response to organic substance	BP	9.23E-16	406
GO:0008152	代谢过程Metabolic process	BP	4.92E-15	8200
GO:0005667	转录因子复合物Transcription factor complex	CC	9.68E-16	763
GO:0071944	细胞外围Cell periphery	CC	4.83E-15	1239
GO:0016020	膜Membrane	CC	1.29E-13	3886
GO:0031224	膜的固有成分Intrinsic component of membrane	CC	3.45E-12	1911
GO:0048046	质外体Apoplast	CC	3.98E-12	151
GO:0016021	膜的组成部分Integral component of membrane	CC	5.1E-12	1869
GO:0005773	液泡Vacuole	CC	5.26E-11	287
GO:0005886	质膜Plasma membrane	CC	1.71E-10	945
GO:0005680	后期促进复合体Anaphase-promoting complex	CC	8.84E-10	180
GO:0000152	核泛素连接酶复合物Nuclear ubiquitin ligase complex	CC	1.24E-09	180
GO:0005618	细胞壁Cell wall	CC	8.72E-09	236
GO:0044425	膜部件Membrane part	CC	1.36E-07	2276
GO:0031461	Cullin-RING泛素连接酶复合物Cullin-RING ubiquitin ligase complex	CC	1.55E-07	185
GO:0031012	细胞外基质Extracellular matrix	CC	4.46E-07	83
GO:0030312	细胞外矩阵External encapsulating structure	CC	6.09E-07	300
GO:0000151	泛素连接酶复合物Ubiquitin ligase complex	CC	9.98E-07	204
GO:0005578	蛋白质细胞外基质Proteinaceous extracellular matrix	CC	3.23E-05	69
GO:0005811	脂质颗粒Lipid particle	CC	0.000125	21
GO:0005576	胞外区Extracellular region	CC	0.000181	751
GO:0000325	植物型液泡Plant-type vacuole	CC	0.000289	51
GO:0020037	血红素结合Heme binding	MF	1.55E-34	443
GO:0004601	过氧化物酶活性Peroxidase activity	MF	4.02E-34	185
GO:0046906	四吡咯结合Tetrapyrrole binding	MF	2.48E-33	449
GO:0016491	氧化还原酶活性Oxidoreductase activity	MF	2.32E-32	1658
GO:0016684	氧化还原酶活性,作用于过氧化物作为受体 Oxidoreductase activity, acting on peroxide as acceptor	MF	7.59E-29	186
GO:0001071	核酸结合转录因子活性 Nucleic acid binding transcription factor activity	MF	1.08E-27	692
GO:0003700	转录因子活性,序列特异性DNA结合 Transcription factor activity, sequence-specific DNA binding	MF	1.08E-27	692
GO:0016705	氧化还原酶活性,作用于成对供体,分子氧掺入或减少 Oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen	MF	6.14E-27	452
GO:0003824	催化活性Catalytic activity	MF	2.78E-22	7173
GO:0005506	铁离子结合Iron ion binding	MF	7.16E-22	395
GO:0016209	抗氧化活性Antioxidant activity	MF	3.69E-19	212
GO:0051213	双加氧酶活性Dioxygenase activity	MF	6.87E-16	163
GO:0004553	水解酶活性,水解邻糖基化合物 Hydrolase activity, hydrolyzing O-glycosyl compounds	MF	1.19E-15	421
GO:0043169	阳离子结合Cation binding	MF	3.65E-15	2408
GO:0016798	水解酶活性,作用于糖基键 Hydrolase activity, acting on glycosyl bonds	MF	6.31E-15	454
GO:0004497	单加氧酶活性Monooxygenase activity	MF	1.22E-13	229
GO:0016161	β-淀粉酶活性Beta-amylase activity	MF	1.29E-13	32
GO:0005509	钙离子结合Calcium ion binding	MF	1.78E-13	313
GO:0046872	金属离子结合Metal ion binding	MF	3.49E-13	2295
GO:0043565	序列特异性DNA结合Sequence-specific DNA binding	MF	2.5E-12	375

附表2

附表3

叶片中显著富集的GO条目"

GO条目GO_trem	描述 Description	条目类型 Term_type	P值矫正值FDR	DEGs数量 DEGs number
GO:0015979	光合作用Photosynthesis	BP	4.87E-29	201
GO:0055114	氧化还原法Oxidation-reduction process	BP	5.7E-24	946
GO:0015995	叶绿素生物合成过程Chlorophyll biosynthetic process	BP	7.37E-18	44
GO:0009628	对非生物刺激的反应Response to abiotic stimulus	BP	7.85E-16	295
GO:0044699	单生物过程Single-organism process	BP	4.78E-14	3540
GO:0009725	对激素的反应Response to hormone	BP	6.56E-14	203
GO:0009409	对寒冷的反应Response to cold	BP	2.02E-13	87
GO:0009719	对内源性刺激的反应Response to endogenous stimulus	BP	2.33E-13	207
GO:0009416	对光刺激的反应Response to light stimulus	BP	8.04E-13	119
GO:0009737	脱落酸反应Response to abscisic acid	BP	9.82E-13	103
GO:0009637	对蓝光的反应Response to blue light	BP	7.79E-12	26
GO:0010033	对有机物的反应Response to organic substance	BP	9.31E-12	233
GO:0046148	色素生物合成过程Pigment biosynthetic process	BP	9.39E-12	98
GO:0009415	对水的反应Response to water	BP	1.08E-11	76
GO:0001101	对酸性化学物质的反应Response to acid chemical	BP	1.12E-11	156
GO:0009314	对辐射的反应Response to radiation	BP	1.78E-11	119
GO:1901700	对含氧化合物的反应 Response to oxygen-containing compound	BP	1.78E-11	193
GO:0097305	对酒精的反应Response to alcohol	BP	1.82E-11	109
GO:0010114	对红灯的反应Response to red light	BP	3.16E-11	23
GO:0008152	代谢过程Metabolic process	BP	4.1E-11	4395
GO:0042651	类囊体膜Thylakoid membrane	CC	1.86E-54	199
GO:0044436	类囊体部分Thylakoid part	CC	3.33E-41	267
GO:0009535	叶绿体类囊体膜Chloroplast thylakoid membrane	CC	8.11E-41	138
GO:0055035	质体类囊体膜Plastid thylakoid membrane	CC	8.11E-41	138
GO:0044434	叶绿体部分Chloroplast part	CC	2.83E-40	329
GO:0009579	类囊体Thylakoid	CC	3.32E-40	280
GO:0034357	光合膜Photosynthetic membrane	CC	3.64E-40	246
GO:0009534	叶绿体类囊体Chloroplast thylakoid	CC	6.45E-40	144
GO:0031976	质体类囊体Plastid thylakoid	CC	6.45E-40	144
GO:0044435	质体部分Plastid part	CC	1.09E-39	332
GO:0031984	细胞器小室Organelle subcompartment	CC	1.61E-36	154
GO:0009507	叶绿体Chloroplast	CC	7.17E-32	520
GO:0009536	质体Plastid	CC	1.27E-30	547
GO:0009941	叶绿体包膜Chloroplast envelope	CC	2.35E-19	156
GO:0009526	质体包膜Plastid envelope	CC	4.84E-19	161
GO:0009654	光系统II析氧复合体 Photosystem II oxygen evolving complex	CC	1.37E-16	51
GO:0009532	质体基质Plastid stroma	CC	7.47E-16	158
GO:0009521	光系统Photosystem	CC	1.45E-14	122
GO:0009570	叶绿体基质Chloroplast stroma	CC	1.76E-14	151
GO:0031977	类囊体腔Thylakoid lumen	CC	7.77E-13	39
GO:0016491	氧化还原酶活性Oxidoreductase activity	MF	2.17E-25	934
GO:0048037	辅因子结合Cofactor binding	MF	1.16E-14	364
GO:0003824	催化活性Catalytic activity	MF	1.86E-12	3817
GO:0016705	氧化还原酶活性,作用于成对供体,分子氧掺入或减少 Oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen	MF	5.66E-12	233
GO:0001071	核酸结合转录因子活性Nucleic acid binding transcription factor activity	MF	8.34E-12	356
GO:0003700	转录因子活性,序列特异性DNA结合 Transcription factor activity, sequence-specific DNA binding	MF	8.34E-12	356
GO:0003959	NADPH脱氢酶活性NADPH dehydrogenase activity	MF	3.16E-11	13
GO:0050662	辅酶结合Coenzyme binding	MF	2.01E-10	287
GO:0005506	铁离子结合Iron ion binding	MF	3.4E-10	205
GO:0016851	镁螯合酶活性Magnesium chelatase activity	MF	8.95E-09	12
GO:0051002	连接酶活性,形成氮金属键 Ligase activity, forming nitrogen-metal bonds	MF	8.95E-09	12
GO:0051003	连接酶活性,形成氮金属键,形成配合物 Ligase activity, forming nitrogen-metal bonds, forming coordination complexes	MF	8.95E-09	12
GO:0004497	单加氧酶活性Monooxygenase activity	MF	1.82E-08	127
GO:0043565	序列特异性DNA结合Sequence-specific DNA binding	MF	1.85E-08	210
GO:0003883	CTP合成酶活性CTP synthase activity	MF	2.32E-08	24
GO:0004753	糖精嘌呤脱氢酶活性Saccharopine dehydrogenase activity	MF	4.17E-07	7
GO:0047130	糖精嘌呤脱氢酶(NADP+,L-赖氨酸形成)活性 Sccharopine dehydrogenase (NADP+, L-lysine-forming) activity	MF	4.17E-07	7
GO:0050661	NADP结合NADP binding	MF	4.63E-07	52
GO:0004553	水解酶活性,水解邻糖基化合物 Hydrolase activity, hydrolyzing O-glycosyl compounds	MF	5.58E-07	217
GO:0004471	苹果酸脱氢酶(脱羧)(NAD+)活性 Malate dehydrogenase (decarboxylating) (NAD+) activity	MF	5.58E-07	22

附表3

图3

附表4

根中显著富集的KEGG通路"

通路编号 Pathway ID	通路 Pathway	DEGs数量 DEGs number	p值矫正值 FDR
DR vs CR UP
ko04075	植物激素信号转导Plant hormone signal transduction	206	3.16E-20
ko00052	半乳糖代谢Galactose metabolism	114	2.82E-16
ko00906	类胡萝卜素生物合成Carotenoid biosynthesis	63	1.25E-15
ko00592	α-亚麻酸代谢α-Linolenic acid metabolism	86	1.67E-15
ko00591	亚油酸Linoleic acid metabolism	41	6.81E-11
ko00500	淀粉和蔗糖代谢Starch and sucrose metabolism	182	5.68E-09
ko00564	甘油磷脂代谢Glycerophospholipid metabolism	120	4.38E-07
ko00250	丙氨酸、天冬氨酸和谷氨酸代谢Alanine, aspartate and glutamate metabolism	77	1.84E-06
ko00561	甘油酯代谢Glycerolipid metabolism	84	2.57E-06
ko00600	鞘脂代谢Sphingolipid metabolism	50	1.35E-05
ko00910	氮代谢Nitrogen metabolism	46	1.46E-05
ko00480	谷胱甘肽代谢Glutathione metabolism	102	0.000107079
ko00760	烟酸和烟酰胺代谢Nicotinate and nicotinamide metabolism	32	0.00017074
ko00945	二苯乙烯类、二芳基庚烷类和姜辣素生物合成 Stilbenoid, diarylheptanoid and gingerol biosynthesis	34	0.000260611
ko00710	光合生物的固碳作用Carbon fixation in photosynthetic organisms	85	0.000380085
ko00010	糖酵解/糖异生Glycolysis / Gluconeogenesis	132	0.00063726
ko00590	花生四烯酸代谢Arachidonic acid metabolism	23	0.001180007
ko00562	磷酸肌醇代谢Inositol phosphate metabolism	62	0.001180007
ko00330	精氨酸和脯氨酸代谢Arginine and proline metabolism	58	0.001345202
ko00430	牛磺酸和低牛磺酸代谢Taurine and hypotaurine metabolism	21	0.002349213
ko00565	乙醚脂质代谢Ether lipid metabolism	39	0.002572574
ko00410	β-丙氨酸代谢beta-Alanine metabolism	48	0.002572574
ko00904	β-丙氨酸代谢Diterpenoid biosynthesis	31	0.002789186
ko04626	植物-病原相互作用Plant-pathogen interaction	141	0.002931207
ko00603	鞘糖脂生物合成-globo系列Glycosphingolipid biosynthesis - globo series	17	0.007304209
ko00280	缬氨酸、亮氨酸和异亮氨酸降解Valine, leucine and isoleucine degradation	49	0.015423467
ko00402	苯并恶嗪类生物合成Benzoxazinoid biosynthesis	5	0.020686992
ko01040	不饱和脂肪酸的生物合成Biosynthesis of unsaturated fatty acids	33	0.025968338
ko00400	苯丙氨酸、酪氨酸和色氨酸生物合成 Phenylalanine, tyrosine and tryptophan biosynthesis	44	0.025968338
ko00908	玉米素生物合成Zeatin biosynthesis	20	0.026518753
DR vs CR DOWN
ko00940	苯丙烷生物合成Phenylpropanoid biosynthesis	124	1.19E-31
ko00195	光合作用Photosynthesis	39	4.25E-13
ko00196	光合作用天线蛋白Photosynthesis - antenna proteins	18	4.61E-09
ko00630	乙醛酸和二羧酸代谢Glyoxylate and dicarboxylate metabolism	43	0.00000338
ko00460	氰基氨基酸代谢Cyanoamino acid metabolism	36	0.00000901
ko00730	硫胺素代谢Thiamine metabolism	12	0.000256332
ko04075	植物激素信号转导Plant hormone signal transduction	54	0.000477217
ko00941	类黄酮生物合成Flavonoid biosynthesis	19	0.006976581
ko00260	甘氨酸、丝氨酸和苏氨酸代谢Glycine, serine and threonine metabolism	29	0.007551386
ko00860	卟啉与叶绿素代谢Porphyrin and chlorophyll metabolism	18	0.012474933
ko00905	油菜素类固醇生物合成Brassinosteroid biosynthesis	8	0.030624012

附表4

附表5

叶中显著富集的KEGG通路"

通路编号 Pathway ID	通路 Pathway	DEGs数量 DEGs_Num	p值矫正值 FDR
DL vs CL UP
ko04075	植物激素信号转导Plant hormone signal transduction	141	5.93E-31
ko00052	半乳糖代谢Galactose metabolism	63	2.57E-13
ko00603	鞘糖脂生物合成-globo系列Glycosphingolipid biosynthesis - globo series	24	1.61E-11
ko00500	淀粉和蔗糖代谢Starch and sucrose metabolism	96	3.42E-08
ko00561	甘油酯代谢Glycerolipid metabolism	52	4.28E-08
ko00600	鞘脂代谢Sphingolipid metabolism	33	2.72E-07
ko00250	丙氨酸、天冬氨酸和谷氨酸代谢Alanine, aspartate and glutamate metabolism	43	3.16E-06
ko00906	类胡萝卜素生物合成Carotenoid biosynthesis	25	3.20E-06
ko00430	牛磺酸和低牛磺酸代谢Taurine and hypotaurine metabolism	17	1.23E-05
ko00280	缬氨酸、亮氨酸和异亮氨酸降解Valine, leucine and isoleucine degradatio	36	1.23E-05
ko00330	精氨酸和脯氨酸代谢Arginine and proline metabolism	35	0.000131771
ko00100	甾体生物合成Steroid biosynthesis	31	0.000139675
ko00592	α-亚麻酸代谢alpha-Linolenic acid metabolism	30	0.000209326
ko00909	倍半萜和三萜生物合成Sesquiterpenoid and triterpenoid biosynthesis	16	0.002175606
ko00650	丁酸代谢Butanoate metabolism	17	0.002261337
ko00904	二萜生物合成Diterpenoid biosynthesis	18	0.002565069
ko00760	烟酸和烟酰胺代谢Nicotinate and nicotinamide metabolism	16	0.003174526
ko00402	苯并恶嗪类生物合成Benzoxazinoid biosynthesis	4	0.00822359
ko00232	咖啡因代谢Caffeine metabolism	6	0.014517862
ko00590	花生四烯酸代谢Arachidonic acid metabolism	11	0.020800302
ko00564	甘油磷脂代谢Glycerophospholipid metabolism	46	0.026196513
ko00071	脂肪酸降解Fatty acid degradation	27	0.039857743
ko00620	丙酮酸代谢Pyruvate metabolism	36	0.043040988
DL vs CL DOWN
ko00195	光合作用Photosynthesis	77	4.61E-52
ko00860	卟啉与叶绿素代谢Porphyrin and chlorophyll metabolism	59	3.90E-35
ko00196	光合作用天线蛋白Photosynthesis - antenna proteins	34	2.00E-28
ko00710	光合生物的固碳作用Carbon fixation in photosynthetic organisms	48	9.80E-13
ko00910	氮代谢Nitrogen metabolism	19	2.05E-05
ko00630	乙醛酸和二羧酸代谢Glyoxylate and dicarboxylate metabolism	33	4.86E-05
ko04712	植物昼夜节律Circadian rhythm - plant	19	9.61E-05
ko00740	核黄素代谢Riboflavin metabolism	9	0.000113929
ko00030	戊糖磷酸途径Pentose phosphate pathway	25	0.000347306
ko00902	单萜生物合成Monoterpenoid biosynthesis	10	0.00191122
ko00960	托烷、哌啶和吡啶生物碱的生物合成Tropane, piperidine and pyridine alkaloid biosynthesis	11	0.01412646
ko00190	氧化磷酸化Oxidative phosphorylation	30	0.02482389
ko00480	谷胱甘肽代谢Glutathione metabolism	27	0.0322657

附表5

附表6

木质素生物合成过程涉及的DEGs"

基因编号 Gene ID	CL	DL	CR	DR	注释 Annotation
Cluster-14906.71489	3.070	9.323	14.043	29.253	苯丙氨酸解氨酶Phenylalanine ammonia-lyase
Cluster-14906.139395	0.553	0.923	3.330	6.990	苯丙氨酸解氨酶Phenylalanine ammonia-lyase
Cluster-14906.56879	1.883	1.333	16.627	55.193	苯丙氨酸解氨酶Phenylalanine ammonia-lyase
Cluster-14906.98337	26.197	79.867	63.683	102.497	苯丙氨酸解氨酶Phenylalanine ammonia-lyase
Cluster-14906.76036	0.907	1.157	7.933	24.347	苯丙氨酸解氨酶Phenylalanine ammonia-lyase
Cluster-14906.166499	3.867	30.673	31.877	78.123	反式肉桂酸4-单加氧酶Trans-cinnamate 4-monooxygenase
Cluster-14906.161634	1.840	17.440	10.913	33.137	反式肉桂酸4-单加氧酶Trans-cinnamate 4-monooxygenase
Cluster-14906.183751	2.063	17.230	9.277	20.610	反式肉桂酸4-单加氧酶Trans-cinnamate 4-monooxygenase
Cluster-14906.205461	1.340	32.647	5.007	16.92	反式肉桂酸4-单加氧酶Trans-cinnamate 4-monooxygenase
Cluster-14906.183750	1.547	9.180	3.297	12.877	反式肉桂酸4-单加氧酶Trans-cinnamate 4-monooxygenase
Cluster-14906.161532	12.513	111.977	74.997	213.130	反式肉桂酸4-单加氧酶Trans-cinnamate 4-monooxygenase
Cluster-14906.182827	19.500	159.120	117.697	293.707	反式肉桂酸4-单加氧酶Trans-cinnamate 4-monooxygenase
Cluster-14906.104971	1.913	6.293	4.903	9.997	反式肉桂酸4-单加氧酶Trans-cinnamate 4-monooxygenase
Cluster-14906.135776	2.060	10.270	5.010	14.483	反式肉桂酸4-单加氧酶Trans-cinnamate 4-monooxygenase
Cluster-14906.41002	0.090	0.660	3.777	8.323	4 -香豆酸-辅酶A连接酶4-coumarate--CoA ligase
Cluster-14906.172556	2.697	6.117	7.847	11.320	4 -香豆酸-辅酶A连接酶4-coumarate--CoA ligase
Cluster-14906.244214	0.250	0.320	3.243	8.450	4 -香豆酸-辅酶A连接酶4-coumarate--CoA ligase
Cluster-14906.98045	1.410	1.550	7.030	11.090	4 -香豆酸-辅酶A连接酶4-coumarate--CoA ligase
Cluster-14906.159381	27.833	62.880	74.293	109.827	4 -香豆酸-辅酶A连接酶4-coumarate--CoA ligase
Cluster-14906.35439	0.940	27.317	30.183	60.867	4 -香豆酸-辅酶A连接酶4-coumarate--CoA ligase
Cluster-14906.123210	1.273	0.233	8.260	2.420	4 -香豆酸-辅酶A连接酶4-coumarate--CoA ligase
Cluster-14906.159383	0.007	0.000	1.207	0.233	4 -香豆酸-辅酶A连接酶4-coumarate--CoA ligase
Cluster-14906.60211	0.150	0.507	1.057	6.523	肉桂酰辅酶A还原酶Cinnamoyl-CoA reductase
Cluster-14906.172984	0.070	0.300	0.523	1.110	肉桂酰辅酶A还原酶Cinnamoyl-CoA reductase
Cluster-14906.60213	1.100	4.167	6.277	34.983	肉桂酰辅酶A还原酶Cinnamoyl-CoA reductase
Cluster-14906.172314	1.140	11.927	9.780	35.253	肉桂酰辅酶A还原酶Cinnamoyl-CoA reductase
Cluster-35239.0	0.020	0.023	1.727	3.883	肉桂酰辅酶A还原酶Cinnamoyl-CoA reductase
Cluster-14906.60210	0.040	0.253	0.147	0.873	肉桂酰辅酶A还原酶Cinnamoyl-CoA reductase
Cluster-14906.60209	0.053	0.543	0.433	5.577	肉桂酰辅酶A还原酶Cinnamoyl-CoA reductase
Cluster-14906.42614	0.000	1.003	20.833	3.030	肉桂酰辅酶A还原酶Cinnamoyl-CoA reductase
Cluster-14561.1	0.000	0.000	0.037	2.430	肉桂醇脱氢酶Cinnamyl-alcohol dehydrogenase
Cluster-14561.0	0.000	0.000	0.333	2.470	肉桂醇脱氢酶Cinnamyl-alcohol dehydrogenase
Cluster-14906.55023	0.220	0.290	505.470	115.510	肉桂醇脱氢酶Cinnamyl-alcohol dehydrogenase
Cluster-14906.52428	0.000	0.000	5.457	0.263	肉桂醇脱氢酶Cinnamyl-alcohol dehydrogenase
Cluster-14906.52638	0.000	0.090	55.993	11.630	肉桂醇脱氢酶Cinnamyl-alcohol dehydrogenase
Cluster-14906.48450	0.000	0.000	26.770	5.013	肉桂醇脱氢酶Cinnamyl-alcohol dehydrogenase
Cluster-14906.122136	0.500	0.253	133.550	27.537	肉桂醇脱氢酶Cinnamyl-alcohol dehydrogenase
Cluster-14906.151146	337.100	109.617	58.610	21.810	肉桂醇脱氢酶Cinnamyl-alcohol dehydrogenase
Cluster-14906.148566	0.030	0.730	58.430	11.743	肉桂醇脱氢酶Cinnamyl-alcohol dehydrogenase
Cluster-14906.162254	17.813	166.080	79.383	705.197	过氧化物酶Peroxidase
Cluster-14906.99962	1.957	1.903	251.467	363.743	过氧化物酶Peroxidase
Cluster-14906.51580	0.107	1.793	103.743	314.013	过氧化物酶Peroxidase
Cluster-14906.163436	24.420	148.720	140.720	196.070	过氧化物酶Peroxidase
Cluster-14906.112155	26.703	127.633	111.990	171.277	过氧化物酶Peroxidase
Cluster-14906.162250	7.010	52.423	19.460	158.550	过氧化物酶Peroxidase
Cluster-14906.167188	14.620	104.527	41.543	140.513	过氧化物酶Peroxidase
Cluster-14906.51123	0.407	2.883	34.840	126.203	过氧化物酶Peroxidase
Cluster-14906.75719	0.000	3.357	7.083	121.560	过氧化物酶Peroxidase
Cluster-14906.149942	241.280	43.240	75.630	12.810	过氧化物酶Peroxidase
Cluster-14906.96512	1.513	46.580	54.063	5.037	过氧化物酶Peroxidase
Cluster-14906.122006	1.790	0.703	6.870	1.363	过氧化物酶Peroxidase
Cluster-14906.180977	0.000	1.697	25.100	8.317	过氧化物酶Peroxidase
Cluster-14906.51983	0.000	0.000	16.220	0.277	过氧化物酶Peroxidase
Cluster-14906.242928	0.000	0.000	6.163	1.570	过氧化物酶Peroxidase
Cluster-14906.38206	0.000	0.000	12.043	0.470	过氧化物酶Peroxidase
Cluster-14906.237061	0.303	2.687	284.187	118.177	过氧化物酶Peroxidase
Cluster-14906.60176	2.130	21.483	36.497	2.127	过氧化物酶Peroxidase
Cluster-14906.48834	0.303	9.257	9.500	35.240	莽草酸O-邻羟基肉桂酰转移酶 Shikimate O-hydroxycinnamoyltransferase
Cluster-14906.234705	0.463	1.297	0.083	2.147	莽草酸O-邻羟基肉桂酰转移酶 Shikimate O-hydroxycinnamoyltransferase
Cluster-14906.234707	1.097	2.443	1.410	5.907	莽草酸O-邻羟基肉桂酰转移酶 Shikimate O-hydroxycinnamoyltransferase
Cluster-14906.84135	0.327	2.637	3.597	9.643	莽草酸O-邻羟基肉桂酰转移酶 Shikimate O-hydroxycinnamoyltransferase
Cluster-14906.234709	0.053	0.510	0.317	3.823	莽草酸O-邻羟基肉桂酰转移酶 Shikimate O-hydroxycinnamoyltransferase
Cluster-14906.33828	0.310	2.977	0.703	6.697	莽草酸O-邻羟基肉桂酰转移酶 Shikimate O-hydroxycinnamoyltransferase
Cluster-14906.235478	3.173	1.997	12.850	4.297	莽草酸O-邻羟基肉桂酰转移酶 Shikimate O-hydroxycinnamoyltransferase
Cluster-14906.52653	0.640	0.037	1.370	0.033	莽草酸O-邻羟基肉桂酰转移酶 Shikimate O-hydroxycinnamoyltransferase
Cluster-14906.246300	1.347	0.000	4.387	0.343	莽草酸O-邻羟基肉桂酰转移酶 Shikimate O-hydroxycinnamoyltransferase
Cluster-14906.236102	0.237	0.157	2.400	0.433	莽草酸O-邻羟基肉桂酰转移酶 Shikimate O-hydroxycinnamoyltransferase
Cluster-19452.0	0.013	1.097	2.707	0.757	莽草酸O-邻羟基肉桂酰转移酶 Shikimate O-hydroxycinnamoyltransferase
Cluster-14906.186135	0.000	0.073	6.267	0.997	莽草酸O-邻羟基肉桂酰转移酶 Shikimate O-hydroxycinnamoyltransferase
Cluster-14906.183828	1.253	1.263	1.163	0.070	莽草酸O-邻羟基肉桂酰转移酶 Shikimate O-hydroxycinnamoyltransferase
Cluster-14906.183825	1.657	0.703	8.633	0.897	莽草酸O-邻羟基肉桂酰转移酶 Shikimate O-hydroxycinnamoyltransferase
Cluster-14906.77787	1.153	0.597	2.630	0.257	莽草酸O-邻羟基肉桂酰转移酶 Shikimate O-hydroxycinnamoyltransferase
Cluster-14906.77785	4.030	2.130	19.020	7.520	莽草酸O-邻羟基肉桂酰转移酶 Shikimate O-hydroxycinnamoyltransferase

附表6

图4

图5

附表7

激素信号转导涉及的DEGs"

基因编号 Gene ID	CL	DL	CR	DR	注释 Annotation
Cluster-14906.221505	0.557	1.390	1.167	3.040	DELLA蛋白DELLA protein
Cluster-14906.251505	0.250	1.577	1.577	4.297	DELLA蛋白DELLA protein
Cluster-14906.212188	0.790	2.627	2.263	6.180	DELLA蛋白DELLA protein
Cluster-14906.36214	0.213	0.893	0.387	1.240	DELLA蛋白DELLA protein
Cluster-14906.136410	7.867	103.147	43.327	94.863	赤霉素受体GID1 Gibberellin receptor GID1
Cluster-14906.159942	4.393	72.967	27.303	66.610	赤霉素受体GID1 Gibberellin receptor GID1
Cluster-14906.136409	14.547	322.227	75.783	180.640	赤霉素受体GID1 Gibberellin receptor GID1
Cluster-14906.175210	19.627	493.157	104.097	319.060	赤霉素受体GID1 Gibberellin receptor GID1
Cluster-14906.175209	2.593	55.220	20.887	63.517	赤霉素受体GID1 Gibberellin receptor GID1
Cluster-14906.245283	0.767	0.023	0.370	0.233	F-box蛋白GID2 F-box protein GID2
Cluster-14906.85071	0.590	2.543	1.283	4.213	光敏色素相互作用因子3 Phytochrome-interacting factor 3
Cluster-14906.245085	1.693	8.603	4.467	11.933	光敏色素相互作用因子3 Phytochrome-interacting factor 3
Cluster-14906.171581	2.620	6.073	1.313	0.120	光敏色素相互作用因子3 Phytochrome-interacting factor 3
Cluster-14906.84753	0.657	1.453	2.650	0.227	光敏色素相互作用因子4 phytochrome-interacting factor 4
Cluster-14906.196922	2.823	5.563	2.803	0.527	光敏色素相互作用因子3 Phytochrome-interacting factor 3
Cluster-14906.100642	0.217	0.670	0.317	1.140	光敏色素相互作用因子3 Phytochrome-interacting factor 3
Cluster-14906.94012	0.133	0.430	0.160	1.100	光敏色素相互作用因子3 Phytochrome-interacting factor 3
Cluster-14906.87877	0.547	2.797	0.820	0.300	光敏色素相互作用因子3 Phytochrome-interacting factor 3
Cluster-14906.104785	1.133	1.300	1.580	0.253	光敏色素相互作用因子4 Phytochrome-interacting factor 4
Cluster-14906.156785	51.100	3.237	1.983	0.543	光敏色素相互作用因子4 Phytochrome-interacting factor 4
Cluster-14906.150558	46.077	4.493	2.690	0.983	光敏色素相互作用因子4 Phytochrome-interacting factor 4
Cluster-14906.190884	0.430	0.273	0.447	0.047	生长素反应因子Auxin response factor
Cluster-14906.108900	4.437	12.997	7.723	14.110	生长素反应因子Auxin response factor
Cluster-14906.97388	3.767	3.137	14.040	23.287	生长素内流载体(AUX1-LAX家族) Auxin influx carrier (AUX1 LAX family)
Cluster-14906.167834	4.077	5.457	17.403	29.683	生长素内流载体(AUX1-LAX家族) Auxin influx carrier (AUX1 LAX family)
Cluster-14906.118371	1.780	2.737	10.503	18.667	生长素内流载体(AUX1-LAX家族) Auxin influx carrier (AUX1 LAX family)
Cluster-14906.167832	1.733	2.737	7.550	16.573	生长素内流载体(AUX1-LAX家族) Auxin influx carrier (AUX1 LAX family)
Cluster-14906.124564	1.593	2.467	6.840	16.883	生长素内流载体(AUX1-LAX家族) Auxin influx carrier (AUX1 LAX family)
Cluster-14906.192689	0.120	0.690	1.370	3.903	生长素内流载体(AUX1-LAX家族) Auxin influx carrier (AUX1 LAX family)
Cluster-14906.193369	0.460	2.440	2.243	6.810	生长素内流载体(AUX1-LAX家族) Auxin influx carrier (AUX1 LAX family)
Cluster-14906.11436	0.017	0.217	0.947	0.130	生长素反应性GH3基因家族 Auxin responsive GH3 gene family
Cluster-14906.42600	0.000	0.537	8.167	14.530	生长素反应性GH3基因家族 Auxin responsive GH3 gene family
Cluster-14906.50539	0.260	1.067	7.483	29.600	生长素反应性GH3基因家族 Auxin responsive GH3 gene family
Cluster-14906.50538	0.037	0.253	1.813	7.353	生长素反应性GH3基因家族 Auxin responsive GH3 gene family
Cluster-14906.50540	0.033	0.317	1.067	4.763	生长素反应性GH3基因家族 Auxin responsive GH3 gene family
Cluster-14906.92468	0.087	0.237	0.793	3.657	生长素反应性GH3基因家族 Auxin responsive GH3 gene family
Cluster-14906.213311	1.277	6.033	2.257	19.243	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.205367	2.700	13.417	26.930	84.550	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.140785	11.743	78.487	34.137	237.477	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.140784	5.783	47.150	18.533	138.790	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.175488	0.627	5.257	2.377	9.377	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.209309	0.073	0.620	1.483	0.140	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.69886	1.063	10.667	7.183	41.047	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.214614	0.760	10.533	10.563	89.317	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.52693	0.087	1.150	0.087	0.733	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.234130	0.423	7.470	0.350	3.190	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.181232	0.447	8.980	1.613	5.237	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.175039	0.237	6.097	0.733	3.297	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.95904	0.173	4.337	0.200	1.607	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.181231	1.123	31.097	5.303	19.280	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.28541	0.133	4.007	1.117	6.690	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.211052	0.787	25.223	1.893	9.567	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.171815	0.417	14.793	2.193	8.610	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.211051	1.017	36.383	2.087	14.030	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.83296	0.100	4.510	0.237	2.447	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.37981	0.050	2.787	2.807	55.770	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.101613	11.710	13.753	24.937	44.163	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.30073	0.190	5.890	2.510	4.873	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.123161	1.320	4.520	4.563	9.457	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.261113	0.083	0.643	1.003	2.327	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.104522	1.497	5.473	10.070	26.157	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.205368	3.603	8.410	17.550	49.360	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.55479	4.720	2.567	12.250	36.047	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.104520	1.233	4.450	6.840	22.983	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.66377	0.670	4.227	2.167	8.437	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.51751	0.047	0.167	0.220	0.830	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.236180	0.330	1.427	1.627	6.770	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.224347	0.000	0.480	0.867	4.703	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.84443	0.310	4.580	1.033	9.260	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.252340	0.880	0.427	0.193	1.953	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.3150	0.000	0.110	2.310	28.077	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.3151	0.000	0.057	0.580	9.260	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.175040	0.000	1.167	0.037	0.870	生长素反应蛋白IAA Auxin-responsive protein IAA
Cluster-14906.89778	1.313	8.323	2.383	14.960	SAUR家族蛋白SAUR family protein
Cluster-14906.74371	1.107	25.770	3.400	17.370	SAUR家族蛋白SAUR family protein
Cluster-14906.74563	0.887	3.477	4.590	9.077	SAUR家族蛋白SAUR family protein
Cluster-14906.74564	2.703	2.243	4.177	10.027	SAUR家族蛋白SAUR family protein
Cluster-14906.74562	3.727	6.370	8.090	21.447	SAUR家族蛋白SAUR family protein
Cluster-14906.61729	2.470	4.837	6.120	19.717	SAUR家族蛋白SAUR family protein
Cluster-14906.61742	2.210	6.630	4.213	16.503	SAUR家族蛋白SAUR family protein
Cluster-14906.60146	4.073	30.210	6.277	43.960	SAUR家族蛋白SAUR family protein
Cluster-14906.243068	1.013	31.953	4.980	51.570	SAUR家族蛋白SAUR family protein
Cluster-14906.216156	1.837	0.753	0.080	0.970	SAUR家族蛋白SAUR family protein
Cluster-14906.2322	0.040	0.047	0.077	1.673	SAUR家族蛋白SAUR family protein
Cluster-14906.121117	7.640	6.140	6.117	9.467	调节蛋白NPR1 Regulatory protein NPR1
Cluster-14906.121119	3.090	3.567	3.033	5.747	调节蛋白NPR1 Regulatory protein NPR1
Cluster-14906.96141	0.993	1.450	0.677	1.667	调节蛋白NPR1 Regulatory protein NPR1
Cluster-14906.91484	4.457	19.817	7.900	21.490	调节蛋白NPR1 Regulatory protein NPR1
Cluster-14906.174587	0.063	0.113	30.947	2.330	致病相关蛋白1 Pathogenesis-related protein 1
Cluster-14906.189409	0.000	0.160	28.003	3.120	致病相关蛋白1 Pathogenesis-related protein 1
Cluster-14906.129257	3.483	3.117	61.843	11.880	致病相关蛋白1 Pathogenesis-related protein 1
Cluster-14906.35029	0.000	0.007	0.833	0.017	转录因子TGA Transcription factor TGA
Cluster-14906.35028	0.000	0.000	1.623	0.210	转录因子TGA Transcription factor TGA
Cluster-14906.58847	0.113	0.153	1.743	0.227	转录因子TGA Transcription factor TGA
Cluster-14906.235475	0.437	0.970	1.680	0.287	转录因子TGA Transcription factor TGA
Cluster-14906.204994	0.640	0.147	3.060	0.647	转录因子TGA Transcription factor TGA
Cluster-14906.51255	0.000	0.007	2.253	0.520	转录因子TGA Transcription factor TGA
Cluster-14906.77294	0.407	0.613	1.957	0.463	转录因子TGA Transcription factor TGA
Cluster-14906.213586	1.370	1.773	4.797	1.300	转录因子TGA Transcription factor TGA
Cluster-14906.69390	2.027	3.407	8.877	2.373	转录因子TGA Transcription factor TGA
Cluster-14906.126707	3.753	1.147	25.537	6.910	转录因子TGA Transcription factor TGA
Cluster-14906.126704	4.027	1.530	17.170	4.670	转录因子TGA Transcription factor TGA
Cluster-14906.213952	0.497	0.170	7.107	1.987	转录因子TGA Transcription factor TGA
Cluster-14906.214468	2.427	3.947	9.313	3.327	转录因子TGA Transcription factor TGA
Cluster-14906.184928	0.740	3.387	7.423	5.863	转录因子TGA Transcription factor TGA
Cluster-14906.156973	0.783	4.343	3.070	2.990	转录因子TGA Transcription factor TGA
Cluster-14906.80433	0.657	6.800	18.270	10.683	转录因子TGA Transcription factor TGA
Cluster-14906.17984	0.033	1.503	0.603	1.437	丝氨酸/苏氨酸蛋白激酶 CTR1Serine/threonine-protein kinase CTR1
Cluster-14906.88790	0.113	1.777	1.197	1.940	丝氨酸/苏氨酸蛋白激酶 CTR1Serine/threonine-protein kinase CTR1
Cluster-14906.88795	0.313	3.543	1.733	3.230	丝氨酸/苏氨酸蛋白激酶 CTR1Serine/threonine-protein kinase CTR1
Cluster-14906.243327	0.650	9.103	4.270	8.493	丝氨酸/苏氨酸蛋白激酶 CTR1Serine/threonine-protein kinase CTR1
Cluster-14906.88794	0.643	12.410	4.250	8.900	丝氨酸/苏氨酸蛋白激酶 CTR1Serine/threonine-protein kinase CTR1
Cluster-14906.89642	0.540	2.887	1.983	4.750	丝氨酸/苏氨酸蛋白激酶 CTR1Serine/threonine-protein kinase CTR1
Cluster-14906.163398	4.007	16.237	39.750	9.367	EIN3结合F-box蛋白EIN3-binding F-box protein
Cluster-14906.255933	0.070	0.000	0.030	0.253	乙烯不敏感蛋白2 Ethylene-insensitive protein 2
Cluster-14906.76809	2.777	27.703	7.583	26.123	乙烯不敏感蛋白3 Ethylene-insensitive protein 3
Cluster-14906.39924	0.023	1.180	0.160	0.050	乙烯不敏感蛋白3 Ethylene-insensitive protein 3
Cluster-14906.45915	0.420	0.633	1.433	3.353	乙烯反应转录因子1 Ethylene-responsive transcription factor 1
Cluster-14906.76670	0.067	0.230	1.003	0.073	乙烯受体Ethylene receptor
Cluster-14906.123629	3.440	6.847	17.543	6.620	乙烯受体Ethylene receptor

附表7

图6

附表8

附表9

图7

附表10

表2

抗旱相关特异性模块中核心转录因子的功能注释"

样本 Sample	模块 Module	甘蔗基因ID Gene ID in sugarcane	基因名称 Gene name	基因功能 Gene function
根Root	red	Cluster-14906.26480	ABI4	编码一个包含AP2结构域的假定转录因子。对ABA、渗透胁迫、糖胁迫和干旱响应。 Encodes a putative transcription factor containing an AP2 domain. Expressed in response to ABA, osmotic stress, sugar stress and drought.
		Cluster-14906.49886	ABR1	编码ERF/AP2转录因子家族的ERF (乙烯反应因子)亚家族B-6的一个成员。 Encodes a member of the ERF (ethylene response factor) subfamily B-6 of ERF/AP2 transcription factor family.
		Cluter-14906.198691	BIM1	编码基本螺旋-环-螺旋(bHLH)家族蛋白BIM1 (BES1-相互作用的MYC-LIKE 1), 参与油菜素类固醇信号传导。 Encodes a basic helix-loop-helix (bHLH) family protein BIM1 (BES1-INTERACTING MYC-LIKE 1), involved in brassinosteroid signaling.
		Cluster-14906.14127	bZIP60	bZIP60由一个bZIP DNA结合域和一个假定的跨膜域组成, 在盐胁迫中发挥作用, 但这种反应通过一个信号通路进行。 bZIP60 consists of a bZIP DNA binding domain followed by a putative transmembrane domain, play a role in salt stress but this response goes through a signaling pathway.
		Cluster-14906.17595	ERF27	编码ERF/AP2转录因子家族的ERF (乙烯反应因子)亚家族的一个成员。 Encodes a member of the ERF (ethylene response factor) subfamily of ERF/AP2 transcription factor family.
		Cluster-14906.82449	ERF8	编码ERF/AP2转录因子家族的ERF (乙烯反应因子)亚家族成员, 参与ABA和免疫信号传导。 Encodes a member of the ERF (ethylene response factor) subfamily of ERF/AP2 transcription factor family, involved in both ABA and immune signaling.
		Cluster-14906.17611	LOB4	植物特异性LATERAL ORGAN BOUNDARIES (LOB)结构域基因家族的成员, 参与土壤中植物根系形态的形成。 Members of the plant-specific LATERAL ORGAN BOUNDARIES (LOB) domain gene family, participate in the formation of root morphology of plants in soil.
		Cluter-14906.257000	NAC71	NAC转录因子家族的成员, 通过上调脱落酸(ABA)信号基因的表达来增强对干旱和盐胁迫的抗性。 A member of NAC transcription factor family, which enhances the resistance to drought and salt stress by up regulating the expression of abscisic acid (ABA) signal gene.
样本 Sample	模块 Module	甘蔗基因ID Gene ID in sugarcane	基因名称 Gene name	基因功能 Gene function
根Root	red	Cluster-14906.44675	TIFY 11c-like	它包含一种由非生物胁迫诱导的锌指结构域tify。 It contains tify domain, a kind of zinc-finger domain, which is induced by abiotic stress.
		Cluster-14906.17972	TIFY9	它包含一种由非生物胁迫诱导的锌指结构域tify。 It contains tify domain, a kind of zinc-finger domain, which is induced by abiotic stress.
		Cluster-14906.81425	DIVARICATA	MYB转录因子的一个成员。 A member of MYB transcription factor.
		Cluster-14906.86529	C2H2-1	C2H2型锌指转录因子。 C2H2 type of zinc finger transcription factors.
根Root	pink	Cluter-14906.231981	ERF1	编码ERF/AP2转录因子家族的ERF (乙烯反应因子)亚家族的一个成员。 Encodes a member of the ERF (ethylene response factor) subfamily of ERF/AP2 transcription factor family.
		Cluter-14906.214614	IAA24	与生长素反应元件结合的转录因子。 A transcription factor that binds to an auxin response element.
		Cluter-14906.117403	MYBS3	一种单一的DNA结合重复序列MYB转录因子, 参与植物对冷胁迫的反应。 A single DNA-binding repeat MYB transcription factor, it is involved in plant response to cold stress.
		Cluster-14906.45120	RAP2-13	作为AP2/ERF转录因子的一员, 与缺氧、氧化和渗透胁迫有关。 A member of AP2 / ERF transcription factors, it is related to hypoxia, oxidation and osmotic stress.
根Root	brown	Cluster-14906.37143	CSA	油菜素类固醇诱导的MYB结构域蛋白CSA参与花粉和种子发育。 The MYB domain protein CSA, induced by brassinosteroid, participates in pollen and seed development.
		Cluster-14906.47480	DREB2B	编码脱水反应元件结合(DREB)转录因子的基因。 Genes encoding the dehydration-responsive element-binding (DREB) transcription factors.
		Cluster-14906.47622	HOX29	同源结构域亮氨酸拉链家族成员。 Member of homeodomain-leucine zipper family.
		Cluter-14906.146281	NAC92	NAC转录因子家族的成员。 Member of NAC transcription factor family.
叶Leaf	greenyellow	Cluter-14906.136592	WRKY50	WRKY蛋白。 WRKY proteins.
		Cluter-14906.182530	MADS51	MADS-box基因, 与植物根系发育有关。 MADS-box genes, it is related to the development of plant root system.
		Cluter-14906.111752	Hsf C-1b	热应激转录因子C-1b。 Heat stress transcription factor C-1b.
		Cluter-14906.123506	ERF5	编码ERF/AP2转录因子家族的ERF (乙烯反应因子)亚家族成员, 参与ABA和免疫信号传导。 Encodes a member of the ERF (ethylene response factor) subfamily of ERF/AP2 transcription factor family, involved in both ABA and immune signaling.
		Cluster-14906.17935	HSB1	热休克因子结合蛋白1。 Heat shock factor-binding protein 1.
		Cluster-14906.41976	ODORANT1	R2R3型MYB家族的成员。 A member of the R2R3-type MYB family.

表2

图8

附表11

图9

参考文献 62

[1]	Fang Y, Xiong L. General mechanisms of drought response and their application in drought resistance improvement in plants. Cell Mol Life Sci, 2015, 72: 673-689. doi: 10.1007/s00018-014-1767-0
[2]	Ma Y, Dias M C, Freitas H. Drought and salinity stress responses and microbe-induced tolerance in plants. Front Plant Sci, 2020, 11: 591911. doi: 10.3389/fpls.2020.591911
[3]	曾廷儒.玉米转高粱SbSNAC1基因抗旱性鉴定及其分子机制解析. 中国农业科学院博士学位论文, 北京, 2020.
	Zeng T R.Identification of Drought-resistant Lines and Analysis of Its Molecular Mechanism of Transgenic Sorghum SbSNAC1 Gene in Maize. PhD Dissertation of Chinese Academy of Agricultural Sciences, Beijing, China, 2020. (in Chinese with English abstract)
[4]	Gentile A, Dias L I, Mattos R S, Ferreira T H, Menossi M. MicroRNAs and drought responses in sugarcane. Front Plant Sci, 2015, 6: 58.
[5]	Zhang Y B, Yang S L, Dao J M, Deng J, Shahzad A N, Fan X, Li R D, Quan Y J, Bukhari S H, Zeng Z H. Drought-induced alterations in photosynthetic, ultrastructural and biochemical traits of contrasting sugarcane genotypes. PLoS One, 2020, 15: e0235845. doi: 10.1371/journal.pone.0235845
[6]	罗晟昇, 廖韦卫, 韦海球, 何洪良, 蒋亚琴, 唐利球. 干旱胁迫对甘蔗生长研究进展. 中国热带农业, 2020, (5): 70-73.
	Luo S S, Liao W W, Wei H Q, He H L, Jiang Y Q, Tang L Q. Research progress of drought stress on sugarcane growth. China Trop Agric, 2020, (5): 70-73. (in Chinese with English abstract)
[7]	李海碧, 桂意云, 张荣华, 韦金菊, 杨荣仲, 张小秋, 李杨瑞, 周会, 刘昔辉. 甘蔗抗旱性及抗旱育种研究进展. 分子植物育种, 2019, 17: 3406-3415.
	Li H B, Gui Y Y, Zhang R H, Wei J J, Yang R Z, Zhang X Q, Li Y R, Zhou H, Liu X H. Research progress on drought resistance and drought-resistant breeding of sugarcane. Mol Plant Breed, 2019, 17: 3406-3415. (in Chinese with English abstract)
[8]	Farooq M, Wahid A, Kobayashi N, Fujita D, Basra S J F S D, Plant drought stress: effects, mechanisms and management. Agron Sustain Dev, 2009, 29: 185-212. doi: 10.1051/agro:2008021
[9]	Sallam A, Alqudah A M, Dawood M F A, Baenziger P S, Börner A. Drought stress tolerance in wheat and barley: advances in physiology, breeding and genetics research. Int J Mol Sci, 2019, 20: 3137. doi: 10.3390/ijms20133137
[10]	刘励蔚. 玉米根系发育与抗旱栽培的研究. 吉林农业, 2019, (22): 35.
	Liu L W. Study on root development and drought resistant cultivation of maize. Jilin Agric, 2019, (22): 35. (in Chinese with English abstract)
[11]	Marcos F C C, Silveira N M, Mokochinski J B, Sawaya A, Marchiori P E R, Machado E C, Souza G M, Landell M G A, Ribeiro R V. Drought tolerance of sugarcane is improved by previous exposure to water deficit. J Plant Physiol, 2018, 223: 9-18. doi: 10.1016/j.jplph.2018.02.001
[12]	温琦, 赵文博, 张幽静, 梁塔娜, 张艳欣, 李丽丽, 黄凤兰. 植物干旱胁迫响应的研究进展. 江苏农业科学, 2020, 48(12): 11-15.
	Wen Q, Zhao W B, Zhang Y J, Liang T N, Zhang Y X, Li L L, Huang F L. Research progress of plant response to drought stress. Jiangsu Agric Sci, 2020, 48(12): 11-15. (in Chinese with English abstract)
[13]	Takahashi F, Kuromori T, Sato H, Shinozaki K. Regulatory gene networks in drought stress responses and resistance in plants. Adv Exp Med Biol, 2018, 1081: 189-214. doi: 10.1007/978-981-13-1244-1_11 pmid: 30288711
[14]	Gahlaut V, Jaiswal V, Kumar A, Gupta P K.Transcription factors involved in drought tolerance and their possible role in developing drought tolerant cultivars with emphasis on wheat (Triticum aestivum L.). Theor Appl Genet, 2016, 129: 2019-2042. pmid: 27738714
[15]	Thirumalaikumar V P, Devkar V, Mehterov N, Ali S, Ozgur R, Turkan I, Mueller-Roeber B, Balazadeh S. NAC transcription factor JUNGBRUNNEN1 enhances drought tolerance in tomato. Plant Biotechnol J, 2018, 16: 354-366. doi: 10.1111/pbi.12776 pmid: 28640975
[16]	El-Esawi M A, Alayafi A A. Overexpression of StDREB2 transcription factor enhances drought stress tolerance in cotton(Gossypium barbadense L.) . Genes (Basel), 2019, 10: 142. doi: 10.3390/genes10020142
[17]	Zhao J, Ren W, Zhi D, Wang L, Xia G. Arabidopsis DREB1A/CBF3 bestowed transgenic tall fescue increased tolerance to drought stress. Plant Cell Rep, 2007, 26: 1521-1528. doi: 10.1007/s00299-007-0362-3
[18]	张余. 水稻抗逆相关基因OsEBP89和OsRMT1功能研究. 华中农业大学博士学位论文 , 湖北武汉, 2020.
	Zhang Y. Functional Characterization of Two Stress-related Genes OsEBP89 and OsRMT1 in Rice . PhD Dissertation of Huazhong Agricultural University, Hubei, Wuhan, China, 2020. (in Chinese with English abstract)
[19]	翟楠鑫, 迟会, 夏玥琳, 刘彩月, 裴新梧, 袁潜华. 海南山栏稻抗旱基因转录组分析. 生物技术通报, 2020, 36(12): 12-20. doi: 10.13560/j.cnki.biotech.bull.1985.2020-0422
	Zhai N X, Chi H, Xia Y L, Liu C Y, Pei X W, Yuan Q H. Transcriptome analysis of drought-resistant genes in Hainan Shanlan Upland Rice. Biotechnol Bull, 2020, 36(12): 12-20. (in Chinese with English abstract)
[20]	张文云, 张建诚, 姚景珍. 氮胁迫下小麦叶片转录组分析. 中国农业科技导报, 2020, 22(11): 26-34.
	Zhang W Y, Zhang J C, Yao J Z. Transcriptome analysis of wheat leaves under nitrogen stress. J Agric Sci Technol, 2020, 22(11): 26-34. (in Chinese with English abstract)
[21]	Henry R J, Furtado A, Rangan P. Wheat seed transcriptome reveals genes controlling key traits for human preference and crop adaptation. Curr Opin Plant Biol, 2018, 45: 231-236. doi: 10.1016/j.pbi.2018.05.002
[22]	Nakabayashi R, Yonekura-Sakakibara K, Urano K, Suzuki M, Yamada Y, Nishizawa T, Matsuda F, Kojima M, Sakakibara H, Shinozaki K, Michael A J, Tohge T, Yamazaki M, Saito K. Enhancement of oxidative and drought tolerance in Arabidopsis by overaccumulation of antioxidant flavonoids. Plant J, 2014, 77: 367-379. doi: 10.1111/tpj.12388
[23]	Rodrigues F A, Laia M L D, Zingaretti S M. Analysis of gene expression profiles under water stress in tolerant and sensitive sugarcane plants. Plant Sci, 2009, 176: 286. doi: 10.1016/j.plantsci.2008.11.007
[24]	Belesini A A, Carvalho F M S, Telles B R, De Castro G M, Giachetto P F, Vantini J S, Carlin S D, Cazetta J O, Pinheiro D G, Ferro M I T. De novo transcriptome assembly of sugarcane leaves submitted to prolonged water-deficit stress. Genet Mol Res, 2017, 16: 1-20
[25]	吴嘉云.甘蔗与斑茅后代染色体遗传分析及抗性初步评价. , 福建农林大学博士学位论文, 福建福州 2013.
	Wu J Y.The Chromosome Inheritance for the Hybrid Progeny of S. officinarum L. and E. arundinaceum and Germplasm Resistance Preliminary Evaluation. PhD Dissertation of Fujian Agriculture and Forestry University, Fuzhou, Fujian, China, 2013. (in Chinese with English abstract)
[26]	刀志学, 鄢家俊, 张建波, 常丹, 白史且, 陈智华, 李达旭, 游明鸿, 张玉, 张昌兵, 张劲, 严旭, 胡超, 张瑜. 斑茅野生种质资源的考察与收集. 植物遗传资源学报, 2013, 14: 816-820.
	Dao Z X, Yan J J, Zhang J B, Chang D, Bai S Q, Chen Z H, Li D X, You M H, Zhang Y, Zhang C B, Zhang J, Yan X, Hu C, Zhang Y. Investigation and ppub of wild Erianthus arundinaceum germplasm resources. J Plant Genet Resour, 2013, 14: 816-820. (in Chinese with English abstract)
[27]	黄永吉, 吴嘉云, 刘少谋, 邓祖湖, 符成, 林彦铨, 李奇伟, 黄忠兴, 陈如凯. 基于GISH的甘蔗与斑茅F₁染色体遗传与核型分析. 植物遗传资源学报, 2014, 15: 394-398.
	Huang Y J, Wu J Y, Liu S M, Deng Z H, Fu C, Lin Y Q, Li Q W, Huang Z X, Chen R K. Chromosome genetic and karyotype analysis of F₁ hybrids between Saccharum officinarum and Erianthus. J Plant Genet Resour, 2014, 15: 394-398. (in Chinese with English abstract)
[28]	Grabherr M G, Haas B J, Yassour M, Levin J Z, Thompson D A, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, Chen Z, Mauceli E, Hacohen N, Gnirke A, Rhind N, Di Palma F, Birren B W, Nusbaum C, Lindblad-Toh K, Friedman N, Regev A. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol, 2011, 29: 644-652. doi: 10.1038/nbt.1883 pmid: 21572440
[29]	Li B, Dewey C N. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinf, 2011, 12: 323. doi: 10.1186/1471-2105-12-323
[30]	Love M I, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol, 2014, 15: 550. doi: 10.1186/s13059-014-0550-8
[31]	Young M D, Wakefield M J, Smyth G K, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol, 2010, 11: R14. doi: 10.1186/gb-2010-11-2-r14
[32]	Mao X, Cai T, Olyarchuk J G, Wei L. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics, 2005, 21: 3787-3793. doi: 10.1093/bioinformatics/bti430
[33]	Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinf, 2008, 9: 559. doi: 10.1186/1471-2105-9-559
[34]	Shannon P, Markiel A, Ozier O, Baliga N S, Wang J T, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res, 2003, 13: 2498-2504. doi: 10.1101/gr.1239303 pmid: 14597658
[35]	Livak K J, Schmittgen T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods, 2001, 25: 402-408. doi: 10.1006/meth.2001.1262 pmid: 11846609
[36]	Begcy K, Mariano E D, Lembke C G, Zingaretti S M, Souza G M, Araújo P, Menossi M. Overexpression of an evolutionarily conserved drought-responsive sugarcane gene enhances salinity and drought resilience. Ann Bot, 2019, 124: 691-700. doi: 10.1093/aob/mcz044
[37]	Kim Y, Chung Y S, Lee E, Tripathi P, Heo S, Kim K H. Root response to drought stress in rice (Oryza sativa L.) . Int J Mol Sci 2020, 21: 1513. doi: 10.3390/ijms21041513
[38]	Uga Y, Sugimoto K, Ogawa S, Rane J, Ishitani M, Hara N, Kitomi Y, Inukai Y, Ono K, Kanno N, Inoue H, Takehisa H, Motoyama R, Nagamura Y, Wu J, Matsumoto T, Takai T, Okuno K, Yano M. Control of root system architecture by DEEPER ROOTING 1 increases rice yield under drought conditions. Nat Genet, 2013, 45: 1097-1102. doi: 10.1038/ng.2725
[39]	邱丽丽, 赵琪, 张玉红, 戴绍军. 植物质膜蛋白质组的逆境应答研究进展. 植物学报, 2017, 52: 128-147. doi: 10.11983/CBB16001
	Qiu L L, Zhao Q, Zhang Y H, Dai S J. Research progress on stress response of plant plasma membrane proteome. Chin Bull Bot, 2017, 52: 128-147. (in Chinese with English abstract)
[40]	Liu X, Zhang Y, Yang H, Liang Y, Li X, Oliver M J, Zhang D.Functional aspects of early light-induced protein (ELIP) genes from the desiccation-tolerant moss Syntrichia caninervis. Int J Mol Sci, 2020, 21: 1411. doi: 10.3390/ijms21041411
[41]	Chaves M M, Flexas J, Pinheiro C. Photosynthesis under drought and salt stress: regulation mechanisms from whole plant to cell. Ann Bot, 2009, 103: 551-560. doi: 10.1093/aob/mcn125
[42]	Reddy A R, Chaitanya K V, Vivekanandan M. Drought-induced responses of photosynthesis and antioxidant metabolism in higher plants. italic>J Plant Physiol, 2004, 161: 1189-1202.
[43]	Challabathula D, Zhang Q, Bartels D. Protection of photosynthesis in desiccation-tolerant resurrection plants. J Plant Physiol, 2018, 227: 84-92. doi: 10.1016/j.jplph.2018.05.002
[44]	马文静. 外源NO对干旱胁迫下紫花苜蓿幼苗苯丙烷类代谢及miRNA调控的影响. 甘肃农业大学硕士学位论文, 甘肃兰州, 2019.
	Ma W J. Effects of Exogenous Nitric Oxide on Phenylpropanoid Metabolism and Regulation of miRNAs in Alfalfa Seedlings under Drought Stress. MS Thesis of Gansu Agricultural University, Lanzhou, Gansu, China, 2019. (in Chinese with English abstract)
[45]	Ortega-García F, Peragón J. The response of phenylalanine ammonia-lyase, polyphenol oxidase and phenols to cold stress in the olive tree (Olea europaea L. cv. Picual). J Sci Food Agric, 2010, 89: 1565-1573. doi: 10.1002/jsfa.3625
[46]	刘威.薄皮甜瓜肉桂醇脱氢酶(CAD)在非生物胁迫下参与木质素合成的功能探究. 沈阳农业大学博士学位论文,辽宁沈阳, 2019.
	Liu W.Functional Research of CmCADs in Lignin Biosynthesis in Oriental Melon (Cucumis melo L.) under Abiotic Stresses. PhD Dissertation of Shenyang Agricultural University, Shenyang, Liaoning, China, 2019. (in Chinese with English abstract)
[47]	黄汉唐.木质素代谢相关基因在柚果实发育过程中的作用机制初探. 福建农林大学硕士学位论文, 福建福州, 2020.
	Huang H T. Preliminary Exploration on Mechanism of Lignin Metabolism Related Genes of Pomelo during Fruit Development Stages. MS Thesis of Fujian Agriculture and Forestry University, Fuzhou, Fujian, China, 2020. (in Chinese with English abstract)
[48]	刘少锋.甘蓝型油菜甘油三磷酸脱氢酶(BnaGPDH)基因cDNA的克隆与功能分析. 湖南农业大学硕士学位论文, 湖南长沙, 2018.
	Liu S F. Cloning and Function Analysis of Glycerol-3-Phosphate Dehydrogenase (BnaGPDH) cDNA from Brassica napus. MS Thesis of Hunan Agricultural University,Changsha, Hunan, China 2018. (in Chinese with English abstract)
[49]	Zhao Y, Liu M, He L, Li X, Wang F, Yan B, Wei J, Zhao C, Li Z, Xu J. A cytosolic NAD(+)-dependent GPDH from maize (ZmGPDH1) is involved in conferring salt and osmotic stress tolerance. BMC Plant Biol, 2019, 19: 16. doi: 10.1186/s12870-018-1597-6
[50]	齐岩, 李杰, 李莹, 刘西燕, 冀好布套, 朱延明, 柏锡, 才华. 转OsMAPK4基因烟草的抗旱性研究与遗传分析. 遗传, 2007, 29: 1144-1148.
	Qi Y, Li J, Li Y, Liu X Y, Ji Hao B T, Zhu Y M, Bo X, Cai H. Drought resistance and heredity of transgenic tobacco with OsMAPK4 gene. Hereditas, 2007, 29: 1144-1148. (in Chinese with English abstract)
[51]	Spitzer-Rimon B, Farhi M, Albo B, Cna'ani A, Ben Zvi M M, Masci T, Edelbaum O, Yu Y, Shklarman E, Ovadis M, Vainstein A. The R2R3-MYB-like regulatory factor EOBI, acting downstream of EOBII, regulates scent production by activating ODO1 and structural scent-related genes in petunia. Plant Cell, 2012, 24: 5089-5105. doi: 10.1105/tpc.112.105247
[52]	Wei Q, Luo Q, Wang R, Zhang F, He Y, Zhang Y, Qiu D, Li K, Chang J, Yang G, He G. A wheat R2R3-type MYB transcription factor TaODORANT1 positively regulates drought and salt stress responses in transgenic tobacco plants. Front Plant Sci, 2017, 8: 1374. doi: 10.3389/fpls.2017.01374
[53]	Xiang L, Liu C, Luo J, He L, Deng Y, Yuan J, Wu C, Cai Y. A tuber mustard AP2/ERF transcription factor gene, BjABR1, functioning in abscisic acid and abiotic stress responses, and evolutionary trajectory of the ABR1 homologous genes in Brassica species. Peer J, 2018, 6: e6071. doi: 10.7717/peerj.6071
[54]	Chen Y S, Ho T D, Liu L, Lee D H, Lee C H, Chen Y R, Lin S Y, Lu C A, Yu S M. Sugar starvation-regulated MYBS2 and 14-3-3 protein interactions enhance plant growth, stress tolerance, and grain weight in rice. Proc Natl Acad Sci USA, 2019, 116: 21925-21935. doi: 10.1073/pnas.1904818116
[55]	Kuang J F, Wu C J, Guo Y F, Walther D, Shan W, Chen J Y, Chen L, Lu W J. Deciphering transcriptional regulators of banana fruit ripening by regulatory network analysis. Plant Biotechnol J, 2020, 19: 477-489. doi: 10.1111/pbi.13477
[56]	王爽, 李海英. 植物E3泛素连接酶与非生物胁迫相关研究进展. 中国农学通报, 2020, 36(29): 47-53.
	Wang S, Li H Y. Plant E3 ubiquitin ligase and abiotic stress: research progress. Chin Agric Sci Bull, 2020, 36(29): 47-53. (in Chinese with English abstract)
[57]	Qin Q, Wang Y, Huang L, Du F, Zhao X, Li Z, Wang W, Fu B. A U-box E3 ubiquitin ligase OsPUB67 is positively involved in drought tolerance in rice. Plant Mol Biol, 2020, 102: 89-107. doi: 10.1007/s11103-019-00933-8
[58]	Lim C W, Park C, Kim J H, Joo H, Hong E, Lee S C. Pepper CaREL1, a ubiquitin E3 ligase, regulates drought tolerance via the ABA-signalling pathway. Sci Rep, 2017, 7: 477. doi: 10.1038/s41598-017-00490-4
[59]	Qanmber G, Liu J, Yu D, Liu Z, Lu L, Mo H, Ma S, Wang Z, Yang Z. Genome-wide identification and characterization of the PERK gene family in Gossypium hirsutum reveals gene duplication and functional divergence. Int J Mol Sci, 2019, 20: 1750. doi: 10.3390/ijms20071750
[60]	Silva N F, Goring D R. The proline-rich, extensin-like receptor kinase-1 (PERK1) gene is rapidly induced by wounding. Plant Mol Biol, 2002, 50: 667-685. doi: 10.1023/A:1019951120788
[61]	Wang X, Li B B, Ma T T, Sun L Y, Tai L, Hu C H, Liu W T, Li W Q, Chen K M. The NAD kinase OsNADK1 affects the intracellular redox balance and enhances the tolerance of rice to drought. BMC Plant Biol, 2020, 20: 11. doi: 10.1186/s12870-019-2234-8
[62]	Muhammad T, Zhang J, Ma Y, Li Y, Zhang F, Zhang Y, Liang Y. Overexpression of a mitogen-activated protein kinase SlMAPK3 positively regulates tomato tolerance to cadmium and drought stress. Molecules, 2019, 24: 556. doi: 10.3390/molecules24030556

Metrics

Viewed

Full text

Abstract

Cited

Shared

Discussed

基因ID Gene ID	正向引物 Forward sequence (5°-3°)	反向引物 Reverse sequence (5°-3°)
Cluster-14906.81942	GTGGCGGCCCTTACGAATGA	GCCCGTACCATGCATTGCGA
Cluster-14906.153942	TTCTGCACGCCTGAGAGCTG	CACGCCATCAGCCCATCTGT
Cluster-14906.115135	GCCCATATTGGACCAGAGCCA	GCTCTAGCAGTGGATGTCGCA
Cluster-14906.152933	GATCTGAGGCTGTGCGTGCT	AGGCGCTGCTGCTGATTCTT
Cluster-14906.154818	AGCGAGTGGCTCCTCTGACA	TGGACGTGCCAATTCCAGTTTGT
Cluster-14906.156090	CTTCAACCCGCTGCACAGGA	TCGGCGATGTCGTGTAACCG
Cluster-14906.142139	CTGACGCAGCAGAAGAGGCA	GTCGATCGCTCCACTGCCAT
Cluster-14906.147361	GCCGACGAGCAGAAGGTGAA	ACAGGCCCGAGAACATGCTG
Cluster-14906.157951	GATGTCGGTGCTGGTAAAAT	CAACGACATTATCACACCAT
Cluster-14906.136592	CGAGGAACTACTACCGGTGC	CCGAGTATGGAGGCGAGAAG
Cluster-14906.146281	GGCACACGTGACCTGCTTCT	CGGCGTCTGCTGGTTGAAGA
Cluster-14906.123506	TGGAGCCCACCACGAGGTAT	CTTCCGAACACGCCCTCCAG
GAPDH	CACGGCCACTGGAAGCA	TCCTCAGGGTTCCTGATGCC

样本编号Sample ID	总读数(百万) Total reads (in millions)	总数据 Total data (Gb)	GC (%)	Q30 (%)	比对率 Number of aligned reads to splicing transcriptome (%)
CL_1	71146938	10.67	55.66	92.55	68.12
CL_2	64776612	9.72	55.9	90.08	68.12
CL_3	62901160	9.44	56.77	90.59	67.98
DL_1	58972730	8.85	55.35	91.34	65.55
DL_2	61055990	9.16	54.95	91.87	66.39
DL_3	75358428	11.3	54.26	90.49	65.4
CR_1	61805904	9.27	53.74	92.03	65.93
CR_2	68885630	10.33	54.58	91.75	65.02
CR_3	69318686	10.4	53.87	91.84	66.56
DR_1	72940638	10.94	54.98	89.98	64.31
DR_2	62329192	9.35	55.32	92.12	66.25
DR_3	75884494	11.38	56.23	89.81	63.35

基因编号 Gene ID	CL	DL	CR	DR	注释 Annotation
Cluster-14906.196840	0.000	10.697	0.047	19.380	DREB1C, AP2-EREBP
Cluster-14906.17595	0.107	13.857	1.057	18.227	ERF027, AP2-EREBP
Cluster-14906.14127	0.000	3.050	0.047	15.543	bZIP60, bZLP
Cluster-14906.26480	0.000	15.983	0.463	106.220	ABI4, AP2-EREBP
Cluster-14906.17611	0.840	5.330	2.497	16.907	LBD4, LOB
Cluster-14906.44675	0.567	14.780	2.590	42.597	TIFY3, Tify
Cluster-14906.81425	0.420	7.670	1.930	20.953	MYBR13, MYB
Cluster-14906.47480	0.000	3.107	7.757	30.450	ERF054, AP2-EREBP
Cluster-14906.117403	0.293	2.217	2.930	11.377	MYBS3, MYB
Cluster-14906.214614	0.760	10.533	10.563	89.317	IAA24, AUX/IAA
Cluster-14906.45120	0.000	10.573	9.077	46.850	RAP2-13, AP2-EREBP
Cluster-14906.231981	0.070	2.330	1.720	8.390	ERF1, AP2-EREBP
Cluster-14906.82449	1.070	7.210	4.790	27.327	ERF8, AP2-EREBP
Cluster-14906.257000	0.377	2.343	1.367	4.593	NAC71, NAC
Cluster-14906.198691	0.323	4.170	2.593	10.763	BIM1, bHLH
Cluster-14906.47622	0.183	1.800	32.773	80.773	HOX29, HB
Cluster-14906.44261	0.717	0.750	2.337	327.897	MYB2, MYB
Cluster-14906.146281	0.410	1.570	6.573	52.243	NAC92, NAC
Cluster-14906.209309	0.073	0.620	1.483	0.140	IAA25, AUX/Laa
Cluster-14906.106436	0.170	19.527	0.233	8.913	ERF061, AP2-EREBP
Cluster-14906.111752	0.593	30.990	0.237	5.283	HSF C-1b, HSF
Cluster-14906.182530	1.333	10.267	0.000	0.393	MDA-box 51, MADS
Cluster-14906.123506	0.260	5.080	0.400	1.080	ERF5, AP2-EREBP
Cluster-14906.41976	0.053	30.537	0.127	2.207	ODORANT1, MYB

基于转录组及WGCNA的甘蔗干旱响应调控网络分析

Analysis of drought responsive regulatory network in sugarcane based on transcriptome and WGCNA

RichHTML

PDF

摘要/Abstract

引用本文

使用本文

图/表 22

参考文献 62

相关文章 15

Metrics

本文评价

推荐阅读 10

关于本刊

在线期刊

作者中心

相关单位

联系我们

转录因子家族 Transcription factor family	基因数量 Genes number	DEGs数量DEGs number	转录因子家族Transcription factor family	基因数量 Genes number	DEGs数量 DEGs number
MYB	555	144	Alfin-like	40	4
bHLH	426	87	CPP	40	4
Orphans	416	116	TAZ	40	8
AP2-EREBP	394	117	PLATZ	39	7
WRKY	376	84	IWS1	38	2
bZIP	364	116	BSD	37	4
NAC	339	77	ARID	34	0
C2H2	317	30	GeBP	33	3
FAR1	305	10	zf-HD	33	8
C3H	290	29	EIL	27	10
GRAS	253	56	BES1	25	6
HB	236	53	Sigma70-like	25	11
AUX/IAA	234	50	Rcd1-like	24	0
ABI3VP1	221	22	C2C2-YABBY	21	1
G2-like	208	23	SRS	21	2
TRAF	203	23	TIG	20	1
SET	193	12	CSD	19	4
CCAAT	175	29	DBP	17	10
PHD	169	8	GRF	17	1
SNF2	164	5	PBF-2-like	14	0
GNAT	146	27	SWI/SNF-SWI3	14	2
mTERF	130	0	ARR-B	13	4
ARF	128	13	LUG	13	2
MADS	125	22	MED6	13	11
TUB	112	25	MBF1	12	1
C2C2-GATA	104	10	SOH1	12	1
HSF	90	37	DDT	10	0
Jumonji	87	18	Pseudo ARR-B	10	4
LOB	87	13	C2C2-CO-like	8	1
C2C2-Dof	85	20	CAMTA	8	2
FHA	85	14	VOZ	8	1
Trihelix	83	13	LIM	7	0
Tify	82	39	Coactivator p15	6	0
SBP	81	12	HRT	6	0
OFP	76	12	ULT	6	6
HMG	73	10	RB	5	0
SWI/SNF-BAF60b	67	9	BBR/BPC	4	0
E2F-DP	61	10	MED7	4	0
TCP	56	17	Zn-clus	2	0
RWP-RK	41	3	S1Fa-like	1	0

模块 Module	DEGs数量 DEGs number
green	1476
cyan	96
greenyellow	260
magenta	1097
salmon	216
red	601
blue	978
brown	813
black	569
midnightblue	39
pink	264

基因编号 Gene ID	连通性 kME	模块 Module
Cluster-14906.26480	0.994853	red
Cluster-14906.49886	0.9924	red
Cluster-14906.14127	0.9841	red
Cluster-14906.231981	0.9881	pink
Cluster-14906.47622	0.9801	brown
Cluster-14906.117403	0.9730	pink
Cluster-14906.146281	0.9802	brown
Cluster-14906.45120	0.9729	pink
Cluster-14906.44675	0.9779	red
Cluster-14906.136592	0.961263	greenyellow
Cluster-14906.41976	0.973491	greenyellow
Cluster-14906.123506	0.958913	greenyellow

[1]	张天宇, 王越, 刘影, 周婷, 岳彩鹏, 黄进勇, 华营鹏. 油菜脯氨酸代谢基因家族的生物信息学分析与核心成员鉴定[J]. 作物学报, 2022, 48(8): 1977-1995.
[2]	韩尚玲, 霍轶琼, 李辉, 韩华蕊, 侯思宇, 孙朝霞, 韩渊怀, 李红英. 基于WGCNA发掘谷子穗部类黄酮合成途径调控关键基因[J]. 作物学报, 2022, 48(7): 1645-1657.
[3]	柯丹霞, 霍娅娅, 刘怡, 李锦颖, 刘晓雪. 大豆TGA转录因子基因GmTGA26在盐胁迫中的功能分析[J]. 作物学报, 2022, 48(7): 1697-1708.
[4]	李旭娟, 李纯佳, 吴转娣, 田春艳, 胡鑫, 丘立杭, 吴建明, 刘新龙. 甘蔗HTD2基因的表达特征及基因多态性分析[J]. 作物学报, 2022, 48(7): 1601-1613.
[5]	郭楠楠, 刘天策, 史硕, 胡心亭, 牛亚丹, 李亮. 长链非编码RNA (LncRNA)在印度梨形孢促进大麦根部生长发育中的调控作用[J]. 作物学报, 2022, 48(7): 1625-1634.
[6]	王兴荣, 李玥, 张彦军, 李永生, 汪军成, 徐银萍, 祁旭升. 青稞种质资源成株期抗旱性鉴定及抗旱指标筛选[J]. 作物学报, 2022, 48(5): 1279-1287.
[7]	朱峥, 王田幸子, 陈悦, 刘玉晴, 燕高伟, 徐珊, 马金姣, 窦世娟, 李莉云, 刘国振. 水稻转录因子WRKY68在Xa21介导的抗白叶枯病反应中发挥正调控作用[J]. 作物学报, 2022, 48(5): 1129-1140.
[8]	肖健, 陈思宇, 孙妍, 杨尚东, 谭宏伟. 不同施肥水平下甘蔗植株根系内生细菌群落结构特征[J]. 作物学报, 2022, 48(5): 1222-1234.
[9]	周慧文, 丘立杭, 黄杏, 李强, 陈荣发, 范业赓, 罗含敏, 闫海锋, 翁梦苓, 周忠凤, 吴建明. 甘蔗赤霉素氧化酶基因ScGA20ox1的克隆及功能分析[J]. 作物学报, 2022, 48(4): 1017-1026.
[10]	陈悦, 孙明哲, 贾博为, 冷月, 孙晓丽. 水稻AP2/ERF转录因子参与逆境胁迫应答的分子机制研究进展[J]. 作物学报, 2022, 48(4): 781-790.
[11]	孔垂豹, 庞孜钦, 张才芳, 刘强, 胡朝华, 肖以杰, 袁照年. 不同施肥水平下丛枝菌根真菌对甘蔗生长及养分相关基因共表达网络的影响[J]. 作物学报, 2022, 48(4): 860-872.
[12]	杨宗桃, 刘淑娴, 程光远, 张海, 周营栓, 商贺阳, 黄国强, 徐景升. 甘蔗类泛素蛋白UBL5应答SCMV侵染及其与SCMV-6K2的互作[J]. 作物学报, 2022, 48(2): 332-341.
[13]	张海燕, 解备涛, 姜常松, 冯向阳, 张巧, 董顺旭, 汪宝卿, 张立明, 秦桢, 段文学. 不同抗旱性甘薯品种叶片生理性状差异及抗旱指标筛选[J]. 作物学报, 2022, 48(2): 518-528.
[14]	李玲红, 张哲, 陈永明, 尤明山, 倪中福, 邢界文. 普通小麦颖壳蜡质缺失突变体glossy1的转录组分析[J]. 作物学报, 2022, 48(1): 48-62.
[15]	汪颖, 高芳, 刘兆新, 赵继浩, 赖华江, 潘小怡, 毕晨, 李向东, 杨东清. 利用WGCNA鉴定花生主茎生长基因共表达模块[J]. 作物学报, 2021, 47(9): 1639-1653.