基于RNA-Seq技术的转录组分析揭示向日葵（Helianthus annuus L.）响应盐胁迫的分子机制

《Genes》：Transcriptomic Analysis Based on RNA-Seq Technology Reveals the Molecular Mechanisms of Sunflower (Helianthus annuus L.) Response to Salt Stress

　　

背景/目的：向日葵（Helianthus annuus L.）是全球四大油料作物之一，具有较强抗逆性，但盐胁迫仍是限制向日葵产量与品质提升的重要因素。方法：研究人员以耐盐品种P50和盐敏感品种P29为实验材料，对NaCl处理的根和叶样品进行转录组测序（RNA-Seq），随后通过组装与拼接、功能注释、差异表达分析、富集分析及转录因子（TFs）预测，揭示向日葵耐盐分子机制。结果：研究表明，两个品种分别获得54,860,184和60,601,572条高质量clean reads；组装聚类后共产生110,751条all-unigenes，其中77,536条获得功能注释；鉴定出21,332个差异表达基因（DEGs），包括10,306个上调基因和11,026个下调基因；对15个DEGs进行实时荧光定量PCR（RT-qPCR）验证，与测序数据一致性达93.33%；基因本体（GO）富集分析显示DEGs显著富集于抗氧化酶活性相关通路；京都基因与基因组百科全书（KEGG）富集分析表明DEGs主要参与15条碳水化合物代谢通路，尤其是淀粉与蔗糖代谢；此外鉴定出67个差异表达的转录因子（TFs）家族（含528个DEGs），包括bHLH、AP2/ERF-ERF、MYB、C3H、WRKY、EREBP、B3-ARF、NAC等。结论：该研究构建了向日葵响应盐胁迫的全面转录图谱，系统阐明了耐盐分子机制；耐盐品种P50通过三条核心策略建立高效盐胁迫防御系统：（i）激活抗氧化系统以快速清除过量活性氧（ROS）并减轻氧化损伤；（ii）通过淀粉与蔗糖再分配调控碳水化合物代谢，提供能量与渗透保护以抵御生理干旱；（iii）动员多个转录因子（TFs）家族建立复杂调控网络，实现对下游功能基因的精确控制。

研究背景：向日葵（Helianthus annuus L.）是世界四大油料作物之一，具有较高的经济价值，同时因耐盐、耐旱、适应贫瘠土壤等特点，被广泛用于盐碱地等边际土地的生态修复与开发利用，尤其在中国内蒙古等干旱半干旱地区成为改良盐碱土的首选作物。尽管向日葵在长期进化中形成了应对盐胁迫的内在调控机制，但其耐盐分子机制仍大多未阐明，这限制了向日葵耐盐基因资源的拓展以及通过分子育种提升耐盐性的效率。因此，从分子层面系统解析向日葵响应盐胁迫的机制具有重要理论与应用价值。研究人员在《Genes》发表了上述基于RNA-Seq的转录组研究，以耐盐品种P50与盐敏感品种P29为材料，通过比较转录组分析鉴定差异表达基因（DEGs）、进行GO与KEGG富集分析并预测差异表达的转录因子（TFs），从而系统揭示向日葵耐盐分子机制，为油料作物稳产出和盐碱地可持续利用提供理论基础与候选基因资源。

主要关键技术方法：研究人员以内蒙古农牧业科学院作物所提供的耐盐品种P50与盐敏感品种P29为样本队列；材料经水浸种、蛭石育苗、1/2 MS营养液适应后，用120 mmol/L NaCl处理，设对照；分别于处理后1、2、3天采集根、下胚轴、幼叶混合样品（对照为0天与3天），设3次生物学重复；样品提取总RNA后等量混合，构建cDNA文库，在Illumina HiSeq^TM2000平台测序；生物信息学流程包括过滤获得clean reads，用Trinity V2.0.3从头组装unigenes并经Tgicl去冗余聚类得all-unigenes；通过BLASTx比对NR、Swiss-Prot、KEGG、COG数据库预测编码序列（CDS），无hit者用ESTScan预测；采用FPKM法以|log₂Fold-change|≥1且FDR<0.001筛选DEGs；用GO与KEGG数据库进行超几何检验富集分析（q-value≤0.05显著）；用iTAK V2.0.1预测DEGs中转录因子（TFs）家族；选15个DEGs以向日葵18S rRNA为内参，通过RT-qPCR（2^?ΔΔCT法）验证表达趋势一致性。

3.1. Transcriptome Sequencing and De Novo Assembly（转录组测序与从头组装）：研究人员对P50和P29的cDNA文库测序，分别获得57,232,364和63,276,298条raw reads，过滤后得54,860,184（P50）和60,601,572（P29）条clean reads，总碱基数分别为4,937,416,560 nt和5,454,141,480 nt；单独组装P50得106,275条unigenes（平均长716 nt，N50为1301 nt），P29得119,350条（平均685 nt，N50为1263 nt）；合并组装得110,751条all-unigenes（总长96,970,017 nt，平均876 nt，N50为1409 nt），其中>2000 nt的有1494条（占1.35%），表明组装质量较高。

3.2. Prediction of CDS（编码序列预测）：研究人员将all-unigenes按NR>Swiss-Prot>KEGG>COG优先级做BLASTx（E-value<0.00001），从匹配unigenes中提取71,260条CDS并翻译为肽序列；无数据库hit的unigenes用ESTScan预测得9970条CDS并翻译；共获81,230条CDS，其长度分布经作图展示。

3.3. Identification of DEGs（差异表达基因鉴定）：研究人员以|log₂Ratio|≥1且FDR≤0.001为标准，在P50与P29间鉴定出21,332个DEGs，其中10,306个上调（P50相对P29高表达）、11,026个下调；表达变化倍数（log₂尺度）范围为?17至16，分布图展示了变化频率特征。

3.4. RT-qPCR Validation（RT-qPCR验证）：研究人员选取15个盐胁迫响应DEGs（如CL10185.Contig2_All等）进行RT-qPCR验证，除Unigene12292_All外，其余14个表达趋势与RNA-Seq一致，一致性达93.33%，证实转录组数据的可靠性。

3.5. GO and KEGG Enrichment Analysis of DEGs（DEGs的GO与KEGG富集分析）：研究人员发现DEGs显著富集于55个GO条目（细胞组分17个、分子功能16个、生物过程22个）；分子功能中1733个DEGs（23.19%）显著富集于氧化还原酶、GTP酶、抗氧化、过氧化物酶、L-抗坏血酸氧化酶活性等抗氧化相关条目。KEGG富集显示DEGs显著富集于128条代谢通路（涵盖20个功能类），其中最显著的是碳水化合物代谢（15条通路），以淀粉和蔗糖代谢居首，其次为戊糖与葡萄糖醛酸转换、糖酵解/糖异生等，表明碳水化合物代谢重编程在耐盐中起关键作用。

3.6. Prediction of TFs（转录因子预测）：研究人员从21,332个DEGs中鉴定出67个转录因子（TFs）家族共528个DEGs（302个上调、226个下调）；数量最多的为bHLH家族（44个DEGs：32上调、12下调），其后依次为AP2/ERF-ERF（37个：19上调、18下调）、MYB（36个：21上调、15下调）、C3H（30个：16上调、14下调）、WRKY（25个：14上调、11下调）、EREBP（25个：14上调、11下调）、B3-ARF（25个：13上调、12下调）、NAC（15个：9上调、6下调）等；同一家族内既有上调又有下调成员，显示精细的特异性调控。

讨论部分总结：研究人员指出高质量转录组（N50>1200 nt）与RT-qPCR验证（93.33%一致）保障了分析可靠性。盐胁迫破坏活性氧（ROS）产生与清除平衡，导致氧化损伤；GO富集显示大量DEGs涉及氧化还原酶、抗氧化、过氧化物酶、L-抗坏血酸氧化酶等，P50中多个抗氧化相关基因上调幅度显著高于P29，与前期生理观测的SOD、POD活性增幅更高一致，说明P50更高效激活过氧化物酶、抗坏血酸氧化酶等基因以清除ROS，增强抗氧化防御是向日葵耐盐核心策略之一；GTP酶相关基因富集暗示其在抗氧化防御与调控通路中可作重要信号分子。KEGG显示DEGs显著富集于碳水化合物代谢特别是淀粉与蔗糖代谢，推测盐胁迫下通过加速淀粉降解为可溶性糖，一方面供能，另一方面作为渗透保护剂维持渗透平衡、保护膜系统与蛋白结构，缓解生理干旱，该策略与Prunellae Spica、木豆等研究一致；糖也可能作为信号分子协调整合胁迫响应基因表达。转录因子（TFs）作为基因表达分子开关，研究人员鉴定出67个家族528个DEGs，bHLH数量最多且多数上调于P50，文献表明bHLH可通过维持渗透稳态、促进ROS清除、调控抗性次生代谢物提升耐盐性，故bHLH是核心TF家族之一；AP2/ERF-ERF、MYB、C3H、WRKY、EREBP、B3-ARF、NAC等亦显著差异表达，这些家族普遍参与植物非生物胁迫响应，可能协同构建向日葵盐胁迫调控网络；同一家族内上调与下调并存，说明并非整个家族单纯激活或抑制，而是通过特定成员精确调控下游靶标，形成高效有序防御网络，与棉花等研究结果类似。研究人员也指出局限：仅用根叶混合样品，未来需组织特异性转录组；仅转录水平，需整合蛋白组、代谢组多维组学；关键候选基因需用转基因或基因编辑做功能验证。

结论翻译：该研究构建了向日葵响应盐胁迫的全面转录图谱，系统阐明了耐盐分子机制；耐盐向日葵品种P50通过三条核心策略实现高效盐胁迫防御系统：（i）激活抗氧化系统以快速清除过量活性氧（ROS）并减轻氧化损伤；（ii）通过淀粉与蔗糖再分配调控碳水化合物代谢，提供能量与渗透保护以抵御生理干旱；（iii）动员多个转录因子（TFs）家族建立复杂调控网络，实现对下游功能基因的精确控制。