杨树根蘖发生的分子机制：生长素-细胞分裂素-氮信号协同调控研究学术报告
一、 研究作者、机构及发表信息
本研究的主要作者为Hongying Pang、Wanwan Lyu、Yajuan Chen、Liping Ding、Lin Zheng和Hongzhi Wang。通讯作者为Lin Zheng和Hongzhi Wang。所有作者均来自中国北京市农林科学院生物技术研究所，隶属于北京市农业基因资源与生物技术重点实验室和北京市作物分子设计与智能育种重点实验室。
该研究成果以题为“Coordinated auxin–cytokinin–nitrogen signaling orchestrates root suckering in populus”的学术论文形式，于2025年12月18日在线发表于国际期刊 International Journal of Molecular Sciences (Int. J. Mol. Sci.) 2025年第26卷第12172页。论文于2025年11月15日收稿，经修订后于2025年12月11日被接受。
二、 学术背景与研究目的
本研究属于植物发育生物学与林木分子遗传学交叉领域，聚焦于木本植物无性繁殖的关键过程——根蘖发生。白杨组（Populus section Leuce）树种，如山杨，通常难以通过扦插生根繁殖，因此根蘖繁殖成为其主要的无性繁殖方式之一，具有成本低、操作简单、能复壮无性系并产生健壮后代的优势。根蘖是从根上产生的由不定芽发育而成的萌条，是北美颤杨等树种在干扰后重建林分的主要机制。尽管已知遗传背景、内源激素水平、养分有效性及环境条件影响根蘖发生，但其潜在的分子机制尚不明确。阐明这一机制对于在林木育种和栽培中，有目的地促进理想再生或抑制不想要的萌蘖（包括转基因逃逸）具有重要指导意义。
先前生理学研究表明，外源生长素（Auxin）应用抑制根蘖萌发，而外源细胞分裂素（Cytokinin）则促进其发生。在拟南芥等模式植物中，生长素与细胞分裂素的互作以及核心分生组织调节因子如WUSCHEL (WUS) 和 SHOOT MERISTEMLESS (STM) 在茎顶端分生组织（Shoot Apical Meristem， SAM）的从头形成与维持中的核心作用已被深入解析。然而，这些已知的调控网络如何在木本植物根上被重编程以启动根蘖发生，特别是其动态的分子时序图谱以及氮信号如何整合到这一过程中，仍属未知。
因此，本研究旨在：1）建立一种高效、快速的杨树根蘖诱导体系；2）利用高时间分辨率的转录组学，解析根蘖发生过程中的动态基因表达变化，划分发育阶段；3）鉴定调控根蘖发生的关键通路与候选基因；4）通过药理学实验验证激素功能；5）最终揭示生长素、细胞分裂素与氮信号协同调控杨树根蘖发生的分子途径。
三、 详细研究流程与方法
本研究包含一系列连贯的实验流程，从材料建立到分子验证，逻辑严密。
流程一：快速根蘖诱导系统建立与细胞学观察 * 研究对象与样本量：使用三年生山新杨（Populus davidiana × P. bolleana）的根段作为研究材料。在秋季落叶后采集根系，切成约10厘米长的根段。 * 处理与实验方法：将根段垂直插入湿润的基质中（埋入约5-7厘米，部分露出地面），置于25°C、70±5%相对湿度的温室黑暗条件下培养。此方法可在6天内高效诱导出大量根蘖，相比传统方法需要数月时间，是一个显著的优化。 * 细胞学分析：在培养第3天（D3），取出现白色椭圆状凸起的根段组织，进行固定、LR White树脂包埋、半薄切片（3 µm）和甲苯胺蓝O染色，利用光学显微镜观察根蘖原基的解剖结构。 * 数据分析流程：通过显微镜成像，描述原基的细胞学特征，特别是与根维管形成层的连接以及顶端分生组织细胞簇的出现。
流程二：时间序列转录组测序与动态分析 * 样本采集：在根段从田间移至温室培养后，每天（D0至D6）上午11点取样。在肉眼可见芽出现前，取根段维管形成层远端的皮层组织；出现原基后，则取整个不定芽原基。每个时间点从9株个体植物取样，混合为3个生物学重复。 * 实验方法：使用植物总RNA提取试剂盒提取总RNA，质检合格（RIN ≥ 8.5）后，构建链特异性cDNA文库，在DNBSEQ-T7平台上进行150 bp双端测序。 * 生物信息学分析： 1. 数据质控与比对：使用FastQC进行原始数据质控，使用HISAT2将高质量读段比对到毛果杨（Populus trichocarpa）v4.1参考基因组。 2. 基因表达量化：使用featureCounts进行基因水平计数，标准化为FPKM值。 3. 样本关系评估：计算样本间皮尔逊相关系数并进行主成分分析（PCA），以评估重复性和识别转录状态。 4. 差异表达基因（DEGs）鉴定：使用edgeR软件包，以错误发现率（FDR）< 0.05 且 |log2倍变化| ≥ 1 为标准，将每个时间点（D1-D6）与起始点（D0）进行比较，鉴定DEGs。 5. 加权基因共表达网络分析（WGCNA）：使用TBtools软件的WGCNA插件，对在所有样本中表达的15,000个基因构建共表达网络。通过自动软阈值选择构建有符号网络，并利用动态树切割法识别基因模块（最小模块大小设为100）。计算模块特征基因与发育阶段的相关性，筛选出与特定阶段显著相关（|r| ≥ 0.44， p < 0.05）的模块进行后续分析。 6. 功能富集分析：对WGCNA筛选出的阶段特异性模块中的基因，进行基因本体（GO）和京都基因与基因组百科全书（KEGG）通路富集分析，以揭示其涉及的生物学过程和代谢通路。
流程三：关键基因表达模式验证 * 实验方法：从用于转录组测序的相同样本中提取RNA，反转录成cDNA。针对转录组分析中揭示的、在激素信号（如ARF5, AHK4）和氮代谢（如ASP1, ASP2）中起关键作用的基因，设计特异性引物，进行逆转录定量聚合酶链式反应（RT-qPCR）验证。 * 数据分析：使用SYBR Green化学法在Bio-Rad CFX系统上进行RT-qPCR，通过熔解曲线分析确认扩增特异性。将RT-qPCR得到的表达模式与RNA-seq数据进行对比，以验证测序结果的可靠性。
流程四：激素功能的药理学验证 * 研究对象与样本量：使用相同的山新杨根段材料，每种处理至少13个生物学重复。 * 实验方法：将根段分别浸入以下无菌溶液中处理：0.5 mg L⁻¹ 或 1.0 mg L⁻¹ 的萘乙酸（NAA，一种生长素类似物），或5.0 mg L⁻¹ 的洛伐他汀（Lovastatin，一种细胞分裂素生物合成抑制剂）。以未处理的根段作为对照。所有处理在25°C黑暗条件下培养6天。 * 数据记录与分析：培养结束后，统计每个根段上萌发的根蘖数量，并拍照记录。使用单因素方差分析（ANOVA）及Duncan’s检验比较不同处理间根蘖数量的差异显著性（p < 0.05）。
流程五：核心分生组织调节因子的表达动态分析 * 分析方法：从时间序列转录组数据中，提取已知在SAM建立和维持中起核心作用的基因（如WUS, STM, CLV1, BAM1, SHR1, WOX4-1/2）的表达谱。 * 数据呈现：通过热图展示这些基因在D0至D6期间的表达变化模式，并结合RT-qPCR对部分基因（如STM, SHR1, WOX4-1）进行表达验证，阐述其在根蘖发生不同阶段的可能功能。
四、 主要研究结果
结果一：快速诱导系统成功建立并揭示根蘖发生的三个连续阶段 解剖学观察显示，培养3天（D3）的根段凸起物由木栓层包裹，内部包含一个根生茎原基，原基基部有两个原形成层束与根维管形成层相连，顶端有一簇分生组织细胞， destined to become the functional SAM。宏观上，D3出现白色椭圆状凸起，D5-6形成具幼叶的芽，D14长出旺盛的枝条。 时间序列转录组PCA分析将样本清晰地分为四个转录状态簇：D0、D1、D2-4和D5-6。据此推断出根蘖发生的三个连续发育阶段：原基起始（D0-1）、SAM建立（D1-4）和器官分化与生长（D4-6）。差异表达分析发现，与D0相比，每个时间点均有约9000个DEGs，其中4913个为所有时间点共有，表明整个过程中存在持续且动态的转录重编程。
结果二：WGCNA揭示阶段特异性核心通路与基因 WGCNA共识别出14个基因模块，其中6个与特定发育阶段显著相关。对这6个模块的功能富集分析揭示了各阶段的分子特征： * 原基起始阶段（D0-D1）：相关模块（绿黄、红、青模块）富集于核糖体生物发生、TCA循环、氧化磷酸化、氨基酸（如色氨酸）生物合成等通路，表明从休眠到高代谢活性的快速转变。关键基因包括色氨酸合成酶基因（TSB1, TSB2, PAT1）、细胞周期相关基因（APC4, APC8, APC10）以及生长素响应因子ARF5，这些基因在早期被强烈诱导，暗示局部生长素合成和细胞周期重启是形成根蘖原基的最早分子事件。 * SAM建立阶段（D1-D4）：相关模块（棕、绿黄模块）显著富集于氮化合物代谢过程。天冬氨酸和谷氨酸代谢关键酶基因（ASP1/2, CTPS, GAD, AK2）被强烈上调。同时，激素响应基因也富集，包括多个生长素响应因子（ARFs）和SAURs，以及细胞分裂素响应基因（RR9-1, RR9-2）和受体基因AHK4。这表明，在SAM建立过程中，生长素-细胞分裂素的互作是在氮代谢的背景下协同进行的。 * 器官分化与生长阶段（D4-D6）：相关模块（蓝、粉、绿黄、红模块）富集于光合作用、碳水化合物生物合成、叶片发育等通路。氮代谢依然突出，尤其是谷氨酰胺合成酶基因（GS1-1, GS1-2, GLN2）在此阶段表达达到峰值。同时，细胞分裂素信号组分（受体AHK1，B型响应调节因子RR4, RR5, RR9）也达到最大表达量，说明持续的细胞分裂素信号和氮代谢共同支撑SAM的维持和持续的器官发生。 RT-qPCR对ARF5、AHK4、ASP1、ASP2等基因的表达验证，支持了RNA-seq数据的可靠性。
结果三：生长素与细胞分裂素通路基因的动态表达及功能验证 对激素合成、代谢和运输相关基因的表达分析显示： * 细胞分裂素：生物合成基因（IPTs, LOGs, CYP735As）在诱导后立即上调，其中CYP735A2（负责催化形成反式玉米素型细胞分裂素）在SAM建立阶段（D2-D4）特异性强烈诱导。而细胞分裂素氧化酶基因CKX5、输出蛋白ABCG14-1等则在整个过程中被抑制，这有利于根中细胞分裂素的积累。LOG1和CKX7等基因在后期（D3-D6）上调，可能参与SAM成熟和维持阶段的局部细胞分裂素水平微调。 * 生长素：生物合成基因TAR2在起始阶段快速诱导，而YUC10则在D6逐渐达到峰值。PIN1和PIN3的多个旁系同源基因差异表达，表明生长素分布的时空控制非常严格。 药理学实验直接验证了这两种激素的功能：外施NAA（0.5和1.0 mg L⁻¹）使根蘖数量分别减少47.5%和59.9%，证实生长素抑制根蘖发育；而施用细胞分裂素生物合成抑制剂洛伐他汀使根蘖萌发减少60%，证明细胞分裂素对根蘖起始是必需的。
结果四：核心分生组织调节因子在根蘖SAM从头形成过程中的重编程 关键分生组织调节基因的表达动态分为三个簇： * 簇I：包含WUS和STM，在D2被快速上调，这与它们作为干细胞身份决定因子在SAM形成中的作用一致。 * 簇II：包含CLV1、BAM1和SHR1，均在D6达到峰值。CLV1/BAM1与CLV3的负反馈环路在拟南芥中已知用于维持分生组织稳态，而SHR1参与木质部模式和形成层活动，暗示在根蘖发生后期及建立后生长中，干细胞稳态与维管发育的协调调控。 * 簇III：包含WOX4-1和WOX4-2，在整个SAM建立和器官分化期（D2-D6）持续诱导。WOX4在杨树中是关键的形成层调节因子，其持续表达表明其在原形成层束分化和新生SAM与母根之间的维管整合中发挥重要作用。
五、 研究结论与意义
本研究通过结合高效的6天根蘖诱导系统与高时间分辨率的转录组分析，首次提供了杨树根蘖发生的整合性分子图谱。研究得出结论：根蘖发生经历原基起始、SAM建立和器官分化三个连续阶段，这一过程由生长素-细胞分裂素-氮信号的协同互作所 orchestrate（精心编排）。 这种互作依次激活了WUS/STM、WOX4-1/2和CLV1/BAM1–SHR1网络，同时驱动了茎顶端分生组织的从头形成以及新生枝条与母根之间的维管整合。
科学价值： 1. 机制解析：首次系统揭示了木本植物根蘖发生的动态分子程序，将已知的激素和分生组织调控网络在时间和空间上进行了精细解析，并突出了氮代谢信号在此过程中的整合作用，提出了“氮信号-细胞分裂素”模块可能通过调节细胞分裂素合成来影响分生组织形成的新假设。 2. 技术方法：建立的快速温室诱导体系，为后续研究提供了高效可靠的实验系统。 3. 提供候选靶点：鉴定出一系列在根蘖发生关键阶段特异性表达的核心基因（如CYP735A2, ARF5, WOX4-1⁄2, ASP1/2, GS1-2等），为通过分子育种手段操纵根蘖发生（增强克隆繁殖效率或抑制人工林中的 unwanted suckers）提供了潜在靶点。
六、 研究亮点
七、 其他有价值内容
研究在讨论部分进行了深入的比较生物学分析。例如，指出虽然本研究提示氮代谢基因（如GS1-2, ASN1）的诱导可能在杨树根蘖中具有重要作用，但在拟南芥和水稻中，同源基因的功能背景存在差异（如拟南芥gs1;1 gs1;2双突变体在低氮下侧根减少，但无明显SAM缺陷），这暗示了木本植物可能存在独特的调控回路。这为后续通过CRISPR/Cas9等基因编辑技术在杨树中进行功能验证指明了方向，以确认这些候选基因在根蘖发生中的具体贡献。此外，研究也指出了未来可利用单细胞和空间转录组学来追踪根蘖干细胞的确切细胞起源，这是对当前研究的有益延伸和展望。