《Frontiers in Plant Science》：Big data approaches to understanding gene regulatory networks in the shoot apical meristem and de novo shoot regeneration
编辑推荐：
基因调控网络（Gene Regulatory Networks, GRNs）通过调控分生组织细胞的分裂、分化及发育命运，控制高等植物茎顶端分生组织（Shoot Apical Meristem, SAM）的发育与功能。这些网络由转录因子及其靶基因构成，同时受植物
基因调控网络（Gene Regulatory Networks, GRNs）通过调控分生组织细胞的分裂、分化及发育命运，控制高等植物茎顶端分生组织（Shoot Apical Meristem, SAM）的发育与功能。这些网络由转录因子及其靶基因构成，同时受植物激素、microRNA及表观遗传调控因子的协同调节。随着多组学与单细胞技术的深度融合，大量与SAM调控相关的数据被系统整合，为解析GRN各组分间的互作提供了重要依据。尽管基因敲除、转基因过表达及报告基因成像等传统实验手段仍对网络解析具有重要贡献，但大规模数据集驱动的GRN构建已显著推动了新假设的提出，并成功用于预测网络拓扑结构与组分互作关系。许多预测的组分与互作随后经计算机模拟（in silico）验证并在活体（in vivo）实验中得到证实，凸显了数据驱动方法在研究调控复杂发育过程的GRNs中的实用价值。本综述重点阐述了部分关键的SAM相关转录因子之间的互作，这些转录因子在正常植物发育及再生过程中的从头（de novo）茎器官形成中均作为网络枢纽发挥作用，强调了整合多种全基因组方法以构建稳健的SAM活性及de novo茎再生的GRNs与预测模型的重要性。
SAM的结构与功能区
茎顶端分生组织是位于茎尖的三维穹顶状结构，包含明确分区以界定不同细胞的功能与命运。中央区（Central Zone, CZ）位于SAM顶端，包含真正的干细胞群；其下方的组织中心（Organising Centre, OC）负责维持CZ干细胞的命运。周边区（Peripheral Zone, PZ）环绕CZ分布，由扩增细胞组成，生长素在此累积并激活侧生器官特异性基因表达程序，启动器官原基发生。SAM同时在垂直方向上分为三个克隆独立的细胞层：L1和L2层仅进行垂周分裂，共同构成表皮层；L3层兼具垂周与平周分裂，构成髓部组织。不同区域细胞群的严格区分确保了CZ干细胞不会过早并入新生器官原基，维持正常发育进程。早期研究依赖荧光报告系与空间微阵列技术定义SAM功能区，近年来逐步转向全基因组与单细胞技术。靶向纯化多核糖体mRNA（TRAP-Seq）技术通过特定启动子驱动标记核糖体，实现了对特定功能区转录组的精准解析，但该策略仍依赖已知标记基因且存在群体平均偏差。单细胞RNA测序（scRNA-seq）突破了这一限制，基于转录组相似性而非空间位置定义细胞群，成功重构了SAM组织的发育轨迹，鉴定了新的发育调控因子，并结合细胞谱系追踪绘制了从合子到成体植物的细胞演化路径。相关可视化工具如ePlant与CoNekT进一步促进了社区对SAM表达数据与共表达网络的探索。
SAM中的调控模块
SAM需协调多类发育过程，研究中常将其拆解为相对独立的调控模块，包括多能性维持、干细胞稳态与侧生器官发生，而这些模块间存在广泛的互作与冗余。多个调控因子同时参与正常发育与体外（in vitro）组织培养中的de novo茎再生过程。植物激素在其中发挥核心作用：生长素（auxin）促进PZ细胞启动侧生器官原基，并介导愈伤组织形成与根再生中的脱分化；细胞分裂素（cytokinin）则在SAM与器官原基中促进细胞增殖与茎器官形成，二者均通过经典信号通路调控下游靶基因的表达。
模块1：SAM中的多能性促进
多能性是SAM细胞的核心属性，指细胞阻滞终末分化、保留多向分化潜能的状态，主要由Ⅰ类KNOTTED1样同源盒（KNOX I）基因编码的TALE类同源域转录因子调控。拟南芥中KNOX I成员包括SHOOT MERISTEMLESS（STM）、KNAT1/BREVIPEDICELLUS（BP）、KNAT2与KNAT3，而Ⅱ类KNOX（KNOX II）基因不参与SAM发育，主要调控叶分化与次生壁合成。KNOX I基因在胚胎发育球型期即开始表达，STM受NAC结构域转录因子CUPSHAPED COTELYDONS（CUC1/2）直接激活，二者共同建立子叶间的分生组织-器官边界。STM通过直接调控KNOX I家族内其他成员、AP2家族基因PLT7/AIL7及HB25等维持多能性，同时抑制促进分化的Ⅱ类CIN-TCP转录因子（TCP3、TCP4、TCP10）表达，防止细胞过早分化。KNOX I基因还广泛调控植物激素途径：通过抑制赤霉素（Gibberellic Acid, GA）生物合成酶基因表达降低GA水平，同时激活异戊烯基转移酶（ISOPENTYL TRANSFERASE, IPT）家族基因促进细胞分裂素合成，进而诱导细胞周期增强因子CYCD3;1表达。染色质免疫沉淀结合转录组分析显示，STM直接调控大量分生组织发育相关转录因子，贝叶斯网络分析进一步揭示了STM与PLT7、CIN-TCP等因子间的调控关系，并经实验验证。
模块2：SAM中央区的干细胞维持
WUSCHEL（WUS）/CLAVATA（CLV）通路是干细胞稳态的核心调控模块，介导OC与CZ间的通讯以维持干细胞数量恒定。WUS编码WOX家族同源域转录因子，其表达受细胞分裂素通过B型ARABIDOPSIS RESPONSE REGULATORS（ARRs）诱导，同时通过抑制A型ARRs增强细胞分裂素响应。WUS在OC表达后通过胞间连丝迁移至CZ的L1与L2层，直接激活CLV3转录。CLV3编码的小分泌肽通过质外体扩散并结合CLR1与CLV2-CORYNE（CRN）富亮氨酸重复受体复合物，启动抑制WUS表达的负反馈信号通路。该反馈受WUS蛋白水平动态调控：低浓度WUS单体通过结合CLV3顺式调控元件（cis-regulatory elements, CREs）中的TAAT核心基序激活其表达，高浓度下WUS形成同源二聚体则抑制CLV3表达。此外，WUS可与HAIRY MERISTEMS1-4（HAM1-4）编码的GRAS结构域转录因子形成异源二聚体，该复合物在OC中抑制CLV3表达，但不直接结合其CREs，主要通过调控WUS的稳定性与扩散能力发挥作用。基于常微分方程（Ordinary Differential Equation, ODE）的“口袋抑制子”模型成功复现了这一调控逻辑，并解释了干细胞稳态与分生组织模式形成的机制。近期多组学研究进一步揭示，组蛋白乙酰转移酶GCN5通过在WUS与CLV3位点添加H3K9Ac活性标记促进其表达，染色质景观的动态变化是WUS-CLV系统功能实现的重要基础。STM与WUS在功能与表达域上存在广泛重叠，STM通过结合CLV3启动子区域的TGACA基序直接调控其表达，且该调控依赖于STM与WUS的二聚化，二者在不同发育阶段动态协同调控CLV3，确保干细胞稳态。
模块3：侧生器官形成
侧生器官发生依赖生长素在SAM侧翼的局部累积形成生长素最大值，该过程由PIN-FORMED（PIN）蛋白介导的极性运输与生长素的正反馈调控实现。器官原基起始伴随多能性基因的稳定抑制与器官特异性分化基因的激活。精确的原基几何定位（拟南芥中为137.5°螺旋叶序）受AP2/ERF转录因子亚家族成员AINTEGUMENTA（ANT）、AIL6/PLT3、AIL5/PLT5与AIL7/PLT7调控，这些PLT因子通过调控生长素生物合成限速酶编码基因YUCCA（YUC1/YUC4）的表达控制生长素最大值形成，且STM直接调控PLT7表达，将干细胞维持与叶序调控偶联。生长素响应因子（Auxin Response Factors, ARFs）如ETTIN、ARF4通过组蛋白去乙酰化直接抑制STM与KNAT1表达，MONOPTEROS（MP）则通过FILAMENTOUS FLOWER（FIL）间接抑制二者。WUS也通过H3K9/K14去乙酰化抑制SAM干细胞中生长素积累，维持低生长素水平以保干细胞特性。器官原基分化需要多能性拮抗因子的作用，如BLADE-ON-PETIOLE1/2（BOP1/2）与ASYMMETRIC LEAVES1/2（AS1/AS2）。STM过表达可上调BOP1/2表达，形成负反馈环。AS1-AS2异源二聚体通过结合KNOX I基因启动子并招募染色质重塑因子HIRA建立抑制性染色质状态，排除KNOX I基因在器官原基中的表达。TCP转录因子是AS2的重要互作伙伴，其中CIN亚族的CIN-TCPs受microRNA miR319转录后调控，通过抑制KNOX I与CUC基因表达促进叶分化。CIN-TCPs与KNOX II家族成员协同抑制STM-CUC模块，防止叶片过度分化为复叶形态。CUC1通过调控PIN1极性定位协调SAM与叶片发育中的生长素最大值形成，是器官边界建立的关键因子。STM与CUC1构成正反馈环路，但STM同时激活靶向CUC1的miR164c，形成不一致前馈环路。STM的细胞间移动能力是该环路空间特异性的基础：STM在CZ高浓度区同时激活CUC1与miR164c，而在器官边界低浓度区仅激活CUC1，导致CUC1表达局限于边界区域。STM的移动受FT INTERACTING PROTEIN（FTIP）3/4调控的内体运输途径介导，其缺失会导致CUC1表达扩散至整个SAM，引发器官融合缺陷。
De Novo茎再生
De novo茎再生是植物损伤响应的核心过程，允许已分化细胞重编程并获得多能性，最终形成新的茎分生组织。该过程分为两个关键阶段：第一阶段为损伤触发伤口响应基因表达，细胞脱分化并重新获得多能性，伴随细胞增殖；第二阶段为茎相关基因表达上调，驱动脱分化细胞获得茎身份，最终建立de novo茎分生组织。体外培养中，外植体先在富含生长素的愈伤组织诱导培养基（Callus Induction Medium, CIM）上形成愈伤组织，再转移至富含细胞分裂素的茎诱导培养基（Shoot Induction Medium, SIM）启动茎身份基因表达，完成再生。
茎再生的遗传调控
多个在SAM发育中发挥功能的基因同样调控de novo再生。损伤早期诱导WOUND INDUCED DEDIFFERENTIATION（WIND）基因（如WIND1）表达，进而上调ENHANCER OF SHOOT REGENERATION 1（ESR1），驱动细胞重编程。PLETHORA（PLT3、PLT5、PLT7）在第一阶段上调PLT1与PLT2，建立再生能力，该过程与侧根原基发育的转录轨迹高度相似，WOX5也参与多能性获得。细胞分裂素信号在特定阶段触发转录开关，抑制侧根发育程序并启动茎细胞命运决定，该时机的精确性决定再生效率。获得再生能力的细胞中，细胞分裂素信号联合ESR1、PLT3/5/7上调WUS、STM、CLV3、CUC1/2等茎特征基因，抑制根相关基因表达。WUS被认为是de novo茎器官发生的核心调控因子，异位表达即可诱导根中发生de novo茎再生。
拟南芥中的茎再生GRNs
早期研究通过增强酵母单杂交（enhanced Yeast One-Hybrid, eY1H）构建了首个植物细胞重编程GRN，涵盖252个转录因子与48个启动子的互作，鉴定出PLT3、ESR1与HEAT SHOCK FACTOR B1（HSFB1）为核心节点。全基因组关联分析（Genome-Wide Association Study, GWAS）利用190份拟南芥种质鉴定出5个再生主效调控因子，其中WUS启动子区的单核苷酸多态性（Single Nucleotide Polymorphism, SNP）通过改变ARR结合基序影响其对细胞分裂素的响应能力，进而调控再生效率。转录调控动力学分析结合染色质可及性与表达谱鉴定出IDD/JACKDAW（JKD）、MYC、PIF、MYB93与NGA2等新正调控因子，其功能缺失导致再生缺陷。
作物中的茎再生GRNs
番茄与小麦的研究显示，de novo再生GRN在物种间具有保守性，WUS、WIND1、ESR1/2及PLT3/7、WOX5等核心组分高度保守。番茄中鉴定出CDF3为新调控因子，可能通过代谢重编程调控植物激素稳态参与再生。小麦研究中发现DOF家族基因TaDOF5.6与TaDOF3.4可显著提升再生效率，为作物遗传转化体系优化提供了新靶点。不同物种在早期愈伤形成阶段的优势转录因子家族存在差异，拟南芥以LBD与NAC为主，小麦则以DOF与G2-like为主。
解析再生与发育中的功能重叠
由于多数再生调控基因同时参与正常发育，解析其在再生中的特异性功能仍是挑战。REGENOMICS数据库整合了大量植物再生相关的转录组数据，支持多物种、多条件的共表达与GRN分析。单细胞转录组学可克服 bulk 数据的群体平均偏差，精准追踪细胞重编程轨迹，未来结合多组学与空间信息将进一步阐明再生GRNs的动态调控逻辑。
未来展望
随着高通量测序技术的发展，多组学大数据的整合极大推动了SAM发育与再生GRNs的解析。机器学习与贝叶斯算法结合染色质免疫沉淀、开放染色质与进化保守性等多元数据，显著提升了GRN推断的准确性。当前研究的局限主要来自bulk样本的时间点“快照”设计，难以捕捉GRNs的动态连续性与跨模块互作。单细胞技术的整合正逐步解决这一问题，未来结合单细胞表观基因组与空间转录组，将实现GRNs在细胞分辨率下的精细解析。公共数据资源如植物公共RNA-seq数据库与拟南芥茎细胞图谱的共享，将进一步加速领域发展，为作物遗传改良与合成生物学设计提供理论基础。
