轻度干旱下多年生黑麦草耐旱基因型与敏感基因型的转录组差异及生长维持机制研究

《Journal of Plant Interactions》：Growth under mild drought: transcriptomics comparison of perennial ryegrass genotypes with differential tolerance to water deficit

摘要 轻度夏季干旱很少威胁多年生黑麦草的存活，但敏感品种可能过早停止生长，导致产量下降。本研究考察了两种叶片生长行为截然不同的基因型在轻度干旱（停止浇水5天）及复水条件下的转录组响应。敏感基因型表现出显著更多的差异表达基因（Differentially Expressed Genes, DEGs），其中包括上调的胁迫相关转录因子（bZIP、DREB、MYB和WRKY）以及下调的细胞生长基因，这与其早期生长停滞和自噬激活相一致。耐旱基因型则表现出延迟的胁迫响应，并上调了与根系发育和水分保持相关的基因（TSO1和WSD1），从而通过生理适应维持持续生长。两种基因型在复水后均激活了恢复通路，尤其是脂肪酸和核苷酸合成，但敏感基因型表现出干旱响应的突然关闭，限制了恢复能力。这些发现揭示了基因型特异的分子策略，并为培育耐旱多年生黑麦草提供了候选基因。 关键政策要点： •轻度干旱（除极端水分胁迫外）对饲草生产力有显著影响，理解植物在该胁迫水平下的响应对于发展气候适应型草地系统至关重要。 •多年生黑麦草基因型在有限水分条件下维持叶片生长的能力存在差异，这凸显了选择能在中度干旱期间维持生物量生产的品种的重要性。 •将生理生长测量与转录组分析相结合，可识别与耐旱性相关的候选分子通路和基因网络，并可能支持改进育种策略的开发。 •本研究中鉴定的生长相关候选基因为开发分子标记提供了潜力，可用于旨在增强耐旱性并在日益多变的气候条件下维持饲草产量的育种项目。 关键词：多年生黑麦草（Lolium perenne）；叶片生长；RNA-seq；轻度干旱；差异表达基因（DEGs）

论文解读：轻度干旱下多年生黑麦草耐旱与敏感基因型的转录组比较研究

一、研究背景与意义

气候变化正重塑温带地区的农业格局，对饲草作物生产力产生深远影响。尽管延长生长季和升温可能增加生物量，但干旱胁迫日益成为北欧-波罗的海地区农业面临的严重问题。多年生黑麦草（Lolium perenne L.）作为全球重要的饲草作物，具有较高的水分需求，其生长和干旱响应存在显著遗传变异。在轻度干旱下，耐旱基因型往往能通过适度的气孔调节、渗透调节和维持光合作用来维持生长，而敏感基因型可能过早停止生长，导致产量下降。然而，以往针对多年生黑麦草的转录组研究多集中于重度干旱，且未将叶片生长停滞与基因表达模式联系起来。本研究旨在填补这一空白，通过比较轻度干旱下生长行为截然不同的两种基因型的转录组响应，揭示与叶片生长维持相关的分子机制，为培育气候适应型多年生黑麦草提供候选基因和分子标记。

二、主要技术方法

研究人员选用立陶宛来源的两种基因型——耐旱基因型'3177'和敏感基因型'3575'，每种基因型通过无性繁殖获得3个克隆植株作为生物学重复。实验在生长室中进行，设置16/8小时光暗周期，停止浇水103小时以模拟轻度干旱。研究人员利用Tri-Phase模型测定叶片伸长速率（LER）以量化基因型对水分亏缺的响应差异，并测定相对含水量（RWC）和光系统II最大光化学效率（Fv/Fm）作为生理指标。在干旱第1、3、4、5天的上午和下午以及复水后2小时采集冠部组织（含叶片分生组织）进行RNA-seq测序。数据分析流程包括FastQC质量控制、fastp修剪、hisat2比对至多年生黑麦草参考基因组（MPB_Lper_Kyuss_1697, ensembl 60）、featureCounts提取计数矩阵、DIAMOND进行功能注释、edgeR的quasi-likelihood方法鉴定差异表达基因（DEGs，p.adj<0.05且fold change≥log₂(1.5)），以及clusterProfiler进行GO和KEGG富集分析。此外，使用WGCNA和multiWGCNA进行共表达网络分析，鉴定生长相关模块中的hub基因。

三、研究结果

3.1 干旱条件下的叶片生长测定（LER）

通过Tri-Phase模型分析，耐旱基因型的叶片伸长速率在土壤水势（Ψ）为3.7 log₁₀(hPa)时开始减缓（Σ），在4.3 log₁₀(hPa)时停止（σ）；而敏感基因型的减缓点和停止点分别为1.6和1.8 log₁₀(hPa)。敏感基因型在更早、更低的土壤水势水平下即减缓并完全停止生长，耐旱基因型在生长减缓点表现出显著更高的LER值。

3.2 轻度干旱下RWC和Fv/Fm的评估

两种基因型的RWC均随干旱进程下降，但敏感基因型在第4天和第5天即显著低于对照，而耐旱基因型仅在第5天才显著下降。敏感基因型的Fv/Fm在第4天和第5天显著降低，而耐旱基因型在整个实验期间保持相对稳定。所有Fv/Fm值均未低于0.75的PSII损伤阈值，表明胁迫程度为轻度。

3.3 功能注释、比对与探索性数据分析

比对率为70%-81%（平均78%），92%的基因获得功能注释，过滤后保留27,342个基因用于差异表达分析。PCA分析显示PC1（46%变异）和PC2（~16%变异）沿对角线分离了基因型和采样天数，PC3（~10%变异）分离了敏感基因型的上午和下午样本。

3.4 轻度干旱和复水期间的差异基因表达

敏感基因型在轻度干旱和复水期间共鉴定出13,703个DEGs，耐旱基因型为3,161个，其中2,673个为两基因型共有。敏感基因型在所有对比中（除D1M vs D3A外）均表现出更多的上调和下调基因。多种已知干旱响应基因和转录因子（如drought-induced 19、ERD15、ERD7、LEA蛋白、HSF A-2c、HSPs、MYBS3 TF等）在两种基因型中均差异表达。KEGG富集分析显示，敏感基因型在干旱期间上调了MAPK信号通路、半乳糖/蔗糖/淀粉代谢和植物激素信号转导等通路，下调了DNA复制、核糖体生物发生和马达蛋白通路。自噬相关通路仅在敏感基因型干旱最后一天富集。

3.5 轻度干旱下基因型间差异

比较转录组分析鉴定出16,215个基因在两种基因型间差异表达。bZIP、DREB、MYB和WRKY转录因子家族相关基因在敏感基因型中高表达。细胞生长和分裂相关基因（如DNA连接酶、细胞周期蛋白和组蛋白编码基因）在敏感基因型中下调，而在耐旱基因型中保持稳定。耐旱基因型中TSO1（抑制MYB3R转录因子以限制根分生组织活性）和WSD1（蜡酯合成酶，参与角质层蜡生物合成）在干旱后期上调，而敏感基因型中TSO1下调。敏感基因型中DHN3（脱水蛋白）、DREB2B和PAT1（ABA依赖性胁迫响应基因产物的降解因子）显著上调。

3.6 共表达分析

共表达分析将DEGs聚类为67个模块，其中6个生长相关模块在基因型间表现出显著差异表达。从每个模块中鉴定出hub基因，包括TUBB1（微管蛋白β-1链）、RHD3（根毛缺陷3）、CDKB2-1（周期蛋白依赖性激酶B2-1）、CSLA1（葡甘聚糖4-β-甘露糖基转移酶1）、XOAT6（木聚糖O-乙酰转移酶6）和TRS130（运输蛋白颗粒复合物II特异性亚基130）。这些基因在敏感基因型中普遍下调，而在耐旱基因型中保持稳定或恢复表达。

四、讨论与结论

讨论部分指出，敏感基因型对轻度干旱表现出更早、更广泛的转录重编程，包括胁迫相关通路的早期激活和生长相关通路的早期抑制，这与叶片生长减缓（σ）和停止（Σ）的提前以及RWC和Fv/Fm的降低一致。耐旱基因型则表现出延迟的胁迫响应，通过上调WSD1介导的蜡生物合成和TSO1介导的根系发育维持水分保持和持续生长。共表达分析揭示的hub基因下调模式进一步表明敏感基因型中细胞分裂、微管动力学和细胞壁生物合成受到更广泛的抑制。从育种角度看，WSD1和PAT1等基因可作为改良耐旱性的潜在靶点。

研究结论部分：对两种研究基因型的整合转录组和生理分析表明，敏感基因型对轻度干旱表现出比耐旱基因型更早、更具破坏性的响应。这一结论得到叶片生长测量的支持，敏感基因型在较高水分可利用性下即表现出更早的生长减缓与停止。这些生长响应与生理指标一致，包括RWC降低和Fv/Fm下降。共表达网络分析进一步揭示，敏感基因型中这种明显的胁迫响应与细胞分裂和细胞壁组织等生长相关基因的协调下调有关。相比之下，耐旱基因型表现出生长相关基因表达的稳定或恢复，使其能在轻度干旱下维持持续生长。这些发现揭示了所研究基因型间不同的干旱响应策略，并鉴定了可用于改善多年生黑麦草耐旱性的分子特征，但需在更多基因型中验证。