整合转录组学鉴定花生幼苗根盐胁迫响应的转录调控网络及生长素相关的AT-hook转录因子AhAHL15

《Current Plant Biology》：Integrative transcriptome analysis identifies a transcriptional regulatory network and an auxin associated transcription factor AhAHL15 underlying salt tolerance in peanut seedling roots

　　

土壤盐渍化是全球农业可持续发展的关键非生物胁迫之一。为探究花生（Arachis hypogaea L.）幼苗根系耐盐性的分子机制，研究人员对太空诱变花生突变系ZHM6的根系进行了盐胁迫处理下的转录组测序（RNA-Seq）。共鉴定出54,430个表达基因和3,817个新转录本。差异表达分析筛选出4,985个差异表达基因（DEGs），其中包含209个差异表达转录因子（DETFs）。通过蛋白-蛋白互作（PPI）网络构建，167个DETFs构成了核心盐胁迫响应调控网络。研究人员进一步功能验证了AhAHL15基因，该基因编码一个含有AT-hook结构域和DUF296结构域的核定位转录因子。异源表达AhAHL15的拟南芥株系在盐胁迫下表现出显著增强的根长和更高的吲哚-3-乙酸（IAA）含量。实时定量PCR（qRT-PCR）分析表明，AhAHL15过表达上调了生长素合成基因YUCCA8和信号转导基因IAA1的表达。因此，AhAHL15通过维持生长素稳态介导盐耐受性，为花生耐盐分子育种提供了候选基因资源。

一、研究背景与科学问题

土壤盐渍化是全球农业面临的主要非生物胁迫之一，全球盐渍化土壤面积已超过10亿公顷并以每年约100万公顷的速度扩张。高浓度的Na⁺和Cl^-会破坏植物细胞渗透稳态、诱导活性氧（ROS）过量积累并抑制根系水分和养分吸收，最终导致作物生长迟缓和产量损失。花生（Arachis hypogaea L.）作为重要的油料和蛋白质作物，对盐胁迫敏感，但其作为异源四倍体作物（2n=4x=40, AABB基因组，~2.7 Gb），富含重复序列和同源基因拷贝，使得其盐胁迫响应调控网络极为复杂。此前的研究多集中于描述差异基因表达，缺乏对转录因子（TF）互作网络的系统性刻画，特别是AT-hook基序核定位（AHL）家族成员在花生盐胁迫中的功能尚属空白。本研究旨在利用太空诱变花生突变系ZHM6构建盐胁迫响应转录组，系统鉴定核心转录因子，并重点解析AHL家族成员AhAHL15的功能。

二、研究方法

研究人员采用以下关键技术：（1）对太空诱变花生突变系ZHM6（第八代突变体，M8，源自地方品种‘Huodou’经神舟十三号飞船搭载六个月空间环境诱变）进行0.5%（w/v）NaCl盐胁迫处理；（2）利用Illumina HiSeq^TM 2500平台进行根系转录组测序（RNA-Seq），使用HISAT2比对参考基因组，StringTie组装转录本，DESeq2鉴定差异表达基因（DEGs）；（3）基于PlantTFDB数据库注释转录因子，利用STRING数据库构建蛋白-蛋白互作（PPI）网络；（4）克隆AhAHL15全长编码序列（CDS），构建35S-AhAHL15-GFP融合载体进行烟草表皮细胞亚细胞定位；（5）通过农杆菌介导的花序浸染法获得AhAHL15过表达（OE）拟南芥纯合株系（T₃代），并进行梯度NaCl胁迫表型分析；（6）采用超高效液相色谱-串联质谱（UPLC-MS/MS）定量植物激素含量。

三、研究结果

3.1 盐胁迫下花生幼苗的表型变异

盐胁迫显著抑制根伸长（对照组平均根长比处理组长1.18 cm），同时增加根直径。纵向切片显示根尖分生区和伸长区结构改变。激素测定表明，盐胁迫下吲哚-3-乙酸（IAA）、脱落酸（ABA）、水杨酸（SA）、茉莉酸（JA）、二氢玉米素（DL-DZ）含量显著下降，而6-苄氨基嘌呤（6-BA）显著升高，其中IAA降幅最大。

3.2 转录组测序数据质量与比对

所有样品产生超过3.8×10¹⁰核苷酸的清洁数据，Q20≥97.47%，Q30≥93.80%。比对到参考基因组的总体映射率大于95.86%，唯一映射率大于81.82%。

3.3 基因表达定量分析与整体特征

共鉴定54,430个非零表达基因，其中对照组特有基因1,584个，处理组特有基因1,749个。主成分分析（PCA）显示两组样品明显分离。染色体表达分布显示末端区域基因表达水平高于中部区域。KEGG富集分析表明，对照组特有基因富集于光合作用通路，处理组特有基因富集于苯丙烷生物合成、二萜生物合成等次级代谢通路。

3.4 新基因挖掘与功能注释

共鉴定3,817个新基因，其中369个为差异表达基因（228个下调，141个上调）。GO富集分析显示，分子功能（MF）类别中最显著富集的是肽酶活性，生物学过程（BP）类别中细胞内转运条目显著富集。

3.5 盐胁迫下根系基因结构变异检测

盐胁迫改变了单核苷酸多态性（SNP）编辑率的频率分布，处理组非同义SNV多于对照组。B亚基因组（BB）表达基因数多于A亚基因组（AA），盐胁迫增加了A亚基因组的相对表达贡献，降低了B亚基因组的贡献。

3.6 差异表达基因筛选与表达模式分析

共鉴定4,985个DEGs（2,268个上调，2,717个下调）。KEGG富集显示，次级代谢物生物合成和代谢途径显著富集。上调DEGs富集于谷胱甘肽代谢、二萜生物合成、苯丙烷生物合成、植物激素信号转导和氮代谢；下调DEGs富集于光合作用、昼夜节律、黄酮和异黄酮生物合成等。共筛选出209个差异表达转录因子（DETFs），其中167个纳入PPI网络。从网络内部区域筛选出12个核心转录因子，包括MYC2（茉莉酸信号）、LHY（生物钟）、WRKY40/WRKY70/WRKY6（防御调控）、AHL15（生长素相关）、BZIP8（ABA和胁迫信号）、MYB86（转录抑制）、BHLH130（胁迫响应bHLH）、SOC1和TCP4（生长发育调控）。

3.7 AhAHL15异源表达增强拟南芥耐盐性

AhAHL15（LOC112733647）编码325个氨基酸，含保守AT-hook结构域和DUF296结构域。亚细胞定位显示AhAHL15-GFP融合蛋白特异性富集于烟草叶表皮细胞核。系统发育分析表明AhAHL15与拟南芥AtAHL15和AtAHL16聚类紧密。组织表达谱显示其在花生下胚轴中表达最高，其次是茎和根。在0.125%、0.25%和0.5% NaCl处理下，过表达株系（OE1和OE2）根长均显著长于野生型（WT）（p<0.01），且随NaCl浓度增加根长减少幅度更平缓。IAA含量测定显示，各盐浓度下OE株系IAA含量显著高于WT，且下降趋势更缓慢。qRT-PCR分析表明，盐胁迫下OE株系中生长素合成基因YUCCA8和信号基因IAA1表达显著上调，而TAA1和ARF7无显著变化。启动子顺式元件分析发现AhAHL15启动子含有bZIP、AP2/RAV和B3/ARF结合位点。

四、讨论与结论

讨论部分指出：盐胁迫诱导花生根系从初级生长相关代谢向胁迫防御和适应方向转变，苯丙烷、二萜和次级代谢物生物合成通路的富集反映了这种适应性策略。新转录本的发现可能归因于异源四倍体花生的基因组复杂性，A和B亚基因组表达偏倚的观察与棉花和小麦中报道的亚基因组表达可塑性一致。167个DETFs构成的PPI网络代表了花生盐胁迫研究中规模最大的TF调控网络，支持“TF枢纽”模型，即多个TF汇聚于共享下游靶标同时调控不同通路特异性基因集。AhAHL15被确认为通过维持生长素稳态发挥作用的耐盐正调控因子，这是豆科植物中首个被功能鉴定的参与非生物胁迫耐受的AHL家族成员。

研究结论：盐胁迫诱导花生幼苗根系发生向次级代谢物生物合成的代谢转变，伴随由WRKY、NAC、bHLH和MYB家族成员构成的多枢纽TF网络。AhAHL15被鉴定并功能表征为通过YUCCA8介导的生长素生物合成途径维持生长素稳态的正向耐盐调控因子。本研究提供了目前最全面的花生根系盐胁迫响应转录组数据集和TF互作网络，并首次在豆科植物中鉴定出参与非生物胁迫耐受的AHL家族成员。这些发现表明，异源四倍体的调控特征——亚基因组表达互补、扩增的TF家族和新转录本储备库——可能共同增强了多倍体作物的胁迫适应能力。AhAHL15及鉴定的枢纽TF可作为标记辅助选择和遗传改良耐盐花生品种的分子靶标。