整合RNA-Seq（RNA测序）、转录因子注释与WGCNA（加权基因共表达网络分析）鉴定陆地棉幼苗短期干旱胁迫下关键候选基因

《Genes》：Integrating RNA-Seq, Transcription Factor Annotation, and WGCNA Identifies Key Candidate Genes in Upland Cotton Seedlings Under Short-Term Drought Stress

　　

背景：干旱是影响并限制棉花生长与产量的主要非生物胁迫之一，然而耐旱与敏旱陆地棉品系间的转录组差异在很大程度上仍属未知。方法与结果：研究人员以耐旱品种XLZ80与敏感品种XLZ61为材料，在干旱胁迫下进行比较表型与转录组分析。表型评估显示，胁迫后XLZ80仅表现轻度叶片萎蔫，而XLZ61出现严重萎蔫。RNA-Seq分析揭示，XLZ80的差异表达基因（DEG）特异性富集于磷脂酰肌醇信号、苯丙氨酸代谢、MAPK（丝裂原活化蛋白激酶）信号及甜菜碱生物合成等通路，而XLZ61的DEG主要参与基础代谢过程。两品系间及内部共鉴定9302个核心DEG，归为8个动态表达簇，含841个转录因子（TF）。WGCNA进一步识别3个与耐旱性相关的关键模块。12个枢纽基因（含GH_D02G2153（MADS-box）与GH_A05G1087（bZIP））被确定为核心调控因子。qRT-PCR（定量实时聚合酶链反应）验证确认这些基因在耐旱品系中诱导更快更强。综上，该研究深化了棉花干旱胁迫响应的转录水平理解，并为抗旱育种提供了宝贵基因资源。

研究背景与立项依据：干旱作为首要非生物胁迫，随气候变化频发重发，严重制约陆地棉（Gossypium hirsutum）苗期生长、产量与纤维品质。传统育种侧重产量与品质，耐旱机制在转录层面尤其在耐旱与敏旱品系间差异尚未明晰。为挖掘苗期根系快速应答的核心调控节点，研究人员以新疆农科院提供的耐旱种XLZ80与敏旱种XLZ61为对照，设置PEG6000模拟干旱0、2、4、6 h多时间点取样，旨在从动态转录组中区分“有效防御”与“紊乱代谢”的遗传基础，研究发表于《Genes》。

主要关键技术方法：研究人员构建XLZ80与XLZ61三叶期苗期根系PEG6000急性干旱时间序列样本队列（每时间点三生物学重复）；采用Illumina HiSeq 2500平台RNA-Seq，以ZJU-AD1_v2.1为参考基因组，用HISAT2+StringTie定量FPKM，DESeq2筛选FDR<0.01且|log₂FC|>1的DEG；利用PlantTFDB注释转录因子；通过WGCNA（软阈值β=8，minModuleSize=30）构建共表达网络并以模块特征基因相关性锁定耐旱模块，取kME最高前4基因为枢纽基因；最后用qRT-PCR（内参GhUBQ7，2^?ΔΔCt法）验证候选基因表达趋势。

研究结果：

3.1 表型分析显示XLZ80较XLZ61更耐旱：PEG6000处理下XLZ80各时点仅轻度萎蔫，XLZ61自2 h即现症状并持续恶化，证实基因型耐旱差异。

3.2 RNA-Seq数据可靠性确认：总clean data 351.31 Gb，Q30>91.24%，同处理重复Pearson r>0.91，PCA按基因型分离，qRT-PCR随机10基因与RNA-Seq趋势一致（r=0.87，p<0.01）。

3.3 时间序列DEG特征：XLZ61共18410 DEG（多于XLZ80的7881），XLZ61特有13901 DEG富集甘油磷脂代谢、光传导等基础/光周期过程；XLZ80特有3372 DEG富集缬亮异亮氨酸合成、磷脂酰肌醇信号、苯丙氨酸代谢、MAPK信号、甜菜碱合成等防御相关通路。

3.4 品系间DEG随干旱持续扩增：共13895品系间DEG，GO富集乙烯合成、角质合成、昼夜节律、蔗糖分解、叶衰老等，KEGG涉光传导、卟啉代谢、泛素介导蛋白水解、黄酮合成、MAPK等。

3.5 9302核心重叠DEG分8簇：簇8（1306基因）在XLZ80持续上调6 h达峰、XLZ61略降，提示该簇为耐旱正向调控群；整体KEGG覆盖钙信号、乙醛酸代谢、碳固定、糖酵解、MAPK。

3.6 841个DEG转录因子分属ERF/MYB/bHLH/NAC/WRKY等11族，分3簇：簇1在XLZ61剧升XLZ80稳；簇3在XLZ80明确上调含bHLH/ERF/MYB/C2H2/NAC，反映耐旱品系稳态型TF程序。

3.7 WGCNA获7模块，brown（XLZ80 6 h正相关）、green（XLZ80 2 h正相关）、yellow（XLZ61 0 h相关）为关键；brown富组氨酸/β-丙氨酸/肌醇磷酸/苯丙氨酸代谢，green富甘油磷脂/类黄酮/胡萝卜素，yellow富钙信号/光传导/昼夜节律。每模块取kME前4得12枢纽基因：GH_A03G1290（NAD(P)氧化还原酶）、GH_D11G0858（YABBY）、GH_D10G0069（肉桂醇脱氢酶）、GH_A06G2102（C2H2，ABA激活信号）、GH_D02G2153（MADS-box）、GH_A05G1087（bZIP）、GH_A07G1345（钙依赖蛋白激酶CDPK）、GH_A04G0096（甘油磷酰二酯磷酸二酯酶）、GH_A12G0213（UDP-糖基转移酶）、GH_D10G0049（几丁质酶）、GH_A12G0644（泛酰巯基乙胺磷酸腺苷转移酶）、GH_D05G3699（抗坏血酸氧化酶）。

3.8 qRT-PCR验证：12基因中多数在XLZ80于2 h即显著诱导且幅值高于XLZ61，TF类（YABBY/C2H2/MADS/bZIP）早诱早峰，证实耐旱品系具备更快更强转录启动能力。

讨论部分总结：研究人员指出XLZ61大量DEG反映剧烈但失衡的转录扰动，偏向基础代谢与细胞周期维持，后期致ROS累积与光合崩溃；XLZ80则经MAPK级联、钙信号、磷脂代谢实现稳态防御。841 TF中NAC/WRKY/bZIP等家族在棉花已有功能佐证，本研C2H2与CDPK枢纽位处网络中心，契合ABA与非ABA途径交叉。局限在于仅苗期根、急旱、双品系，未做过表达/突变体因果验证，后续需扩品系与多时期多器官并联合表观/代谢层。

结论翻译：通过整合苗期干旱表型与RNA-Seq，研究人员发现耐旱XLZ80形态损伤轻，其转录响应有别于敏旱XLZ61，由基因型特异调控主导；XLZ80的DEG集中于磷脂酰肌醇信号系统、MAPK信号通路及次生代谢，XLZ61偏基础代谢。结合WGCNA锁定3个耐旱核心模块与12枢纽基因。该研究为解析棉花耐旱分子机制提供基因资源与理论基础，为基因工程与标记辅助育种增强耐旱性奠基。