关于《整合miRNA‑mRNA‑lncRNA分析揭示香椿根蘖克隆繁殖调控网络》研究的学术报告

本研究于2026年发表在《植物分子生物学报道》（Plant Molecular Biology Reporter）期刊第44卷第53期。研究团队来自中国多所科研单位，主要作者包括 Jiaming Tie（第一作者）、Qiangqiang Li、Qiangqiang Cheng、Qiwei Zhong、Zhonghua Tu、Rongxi Sun、Yao Xiao、Chunce Guo、Lu Zhang 以及通讯作者 Jikai Ma。研究机构包括江西农业大学江西省亚热带森林资源培育重点实验室、江西省本土树种改良与高效利用重点实验室、江西省蔬菜栽培与利用重点实验室，以及江西省林业科学院江西省植物生物技术重点实验室。

研究学术背景 本研究的科学领域属于植物分子生物学、发育生物学与森林遗传学的交叉领域，聚焦于植物无性繁殖（克隆生长）的遗传调控机制。研究对象为香椿（Toona fargesii），一种中国特有且具有高经济价值的濒危树种。香椿主要通过“根蘖”（root sprouting, RS）进行再生，即从根部产生新的芽苗。这种无性繁殖策略虽有利于个体生存和快速占据生态位，但长期依赖可能导致种群遗传多样性下降，削弱物种对环境变化的适应能力，这对于濒危物种的保护尤为不利。尽管根蘖现象在约10%的陆生植物中存在，但其背后的遗传调控网络，特别是非编码RNA（如长链非编码RNA, lncRNA；微小RNA, miRNA）在其中扮演的角色，尚缺乏深入研究。此前研究表明，激素（尤其是生长素与细胞分裂素的比例失衡）以及经典调控模块（如miR156-SPL）在芽的再生与发育中起关键作用，但对lncRNA在根蘖过程中的功能知之甚少。因此，本研究旨在通过全转录组测序技术，整合分析mRNA、lncRNA和miRNA的表达谱，系统揭示香椿根蘖发育背后的潜在遗传调控网络，以阐明其克隆繁殖的分子基础，并为理解这一濒危树种的繁殖策略提供新见解。

详细研究流程 本研究流程严谨，从表型观察到分子机制解析，共包含五个主要阶段。

第一阶段：根蘖现象的动态观察与样品准备。 研究团队首先对江西省南昌市梅岭山脉一片八年生香椿人工林进行了为期一年（2023年1月至12月）的定点观测。他们系统记录了根蘖发生的时间动态、数量以及所附着根系的直径。结果显示，根蘖主要发生在5月至8月，其中6月为高峰期，10月后完全停止。约46.68%的根蘖发生在直径为4.0-6.0毫米的根系上。这些细致的表型数据为后续确定采样关键期和选择合适的实验材料提供了关键依据。为进一步明确根蘖的形态发生特征，研究者使用扫描电子显微镜（SEM）对早期根蘖进行了观察，发现芽原基从根表皮突出，并被密集的毛状体包裹。基于形态学观察，研究团队严格区分了根蘖（RS）组织样品和非根蘖（NRS，即正常根部组织）样品。为满足测序的生物学重复要求，于6月至8月期间，分别采集了RS和NRS样品各约100份，随后等分为两部分，一部分用于构建mRNA和lncRNA测序文库，另一部分用于构建小RNA（sRNA，主要用于miRNA分析）测序文库。每个处理组设置三个生物学重复。

第二阶段：转录组测序与文库构建。 使用RNAprep Pure植物试剂盒提取总RNA，并检测其质量与完整性。对于mRNA和lncRNA测序，首先使用Ribo-Zero rRNA去除试剂盒去除核糖体RNA，然后将RNA片段化并反转录构建cDNA文库。对于sRNA测序，则通过连接接头、反转录和PCR扩增来构建文库。所有文库均在Illumina HiSeq 4000平台上进行高通量测序，最终获得了约81.15 GB的mRNA/lncRNA数据和3.32 GB的sRNA数据。数据质量评估显示Q30均高于90%，表明测序数据可靠。

第三阶段：转录本识别与差异表达分析。 这是数据处理的核心环节，运用了多个生物信息学工具和流程。首先，使用fastp软件对原始测序数据进行质控，获得干净读数（clean reads）。随后，使用HISAT2将干净读数比对到研究团队之前已发表的香椿参考基因组上。接着，使用StringTie进行转录本拼接和定量，并利用其“merge”功能合并所有样本的转录本，过滤掉长度小于200个核苷酸及链方向不确定的转录本。为了从海量转录本中准确识别出lncRNA，研究采用了“多工具交叉验证”的策略：使用CPC2（编码潜力计算器）、CNCI（编码-非编码指数）和Pfam（蛋白质家族数据库）三个独立的工具来预测转录本的编码潜力。只有被三个工具同时预测为非编码的转录本，才被最终确定为新的lncRNA（novel_lncRNAs）。这种严格的筛选流程有效降低了假阳性。最终，共鉴定出14,544条新的lncRNA，并根据其与已知蛋白质编码基因的位置关系分为反义、基因间、正义重叠和正义内含子等四类。对于miRNA，则将测序得到的小RNA标签通过Bowtie比对到参考序列，并参照miRBase 20.0、PmiREN和sRNAanno数据库，利用改进的miRDeep2软件进行已知miRNA的鉴定和新miRNA的预测。在基因表达定量方面，使用StringTie计算每个基因的FPKM值。差异表达分析则采用DESeq2 R包进行，筛选标准为校正后的p值（padj）小于0.05且表达倍数变化（log2FoldChange）绝对值大于1。通过此流程，共鉴定出2,603个差异表达基因（DEGs）、135个差异表达lncRNA（DELs，此处“135”可能指筛选后重点关注的数量或注释到靶基因的数量，文中另一处提及1,692个DELs）和39个差异表达miRNA（DEMs）。

第四阶段：功能富集分析与调控网络构建。 为理解差异表达分子可能参与的生物学过程，研究进行了基因本体（GO）和京都基因与基因组百科全书（KEGG）通路富集分析。DEGs显著富集于“光合作用”和“ABC转运蛋白”等通路，后者提示激素介导的调控（如细胞分裂素转运蛋白ABCG11/14）可能参与根蘖过程。差异表达lncRNA的预测靶基因则在“玉米素生物合成”（与细胞分裂素合成相关）等通路中富集。基于已有文献报道和差异表达分析结果，研究筛选出两个关键枢纽基因：SQUAMOSA启动子结合蛋白类似物10（SPL10）和BRANCHED1（BRC1）。随后，通过多步骤预测构建了miRNA-mRNA-lncRNA共表达调控网络：1）计算DELs与DEGs之间的皮尔逊相关系数（>0.5， FDR<0.01），确定共表达对；2）预测lncRNA在基因组上、下游10 kb范围内的顺式作用靶基因；3）综合以上两种结果，确定lncRNA与mRNA的潜在相互作用；4）使用psRobot软件预测枢纽基因（如BRC1， SPL10）可能靶向它们的miRNA（如miR319， miR156）。最终，利用Cytoscape软件将这些相互作用整合并可视化。

第五阶段：实验验证。 为验证转录组测序结果的可靠性，研究者选取了包括SPL10、BRC1、miR319b、miR156g-3p以及部分lncRNA在内的候选分子，进行了逆转录定量聚合酶链式反应（RT-qPCR）验证。由于香椿基因组近期经历过全基因组复制事件，存在大量同源基因，研究特别设计了特异性引物以避免交叉反应。RT-qPCR的结果与转录组测序的表达趋势高度一致，证实了测序数据的准确性。例如，miR319b在RS中低表达，而在NRS中高表达，与其预测靶基因BRC1的表达模式（RS中高表达）相反，符合miRNA负调控靶基因的预期。

主要研究结果 1. 根蘖表型特征明确： 系统观测明确了香椿根蘖发生的季节性规律（夏季为主）和空间偏好性（中等粗细根系），SEM揭示了其从根表皮发生的早期形态结构。 2. 转录组全景图绘制成功： 高质量完成了香椿根蘖与非根蘖组织的mRNA、lncRNA和miRNA全转录组测序，首次在香椿中系统鉴定了大量lncRNA，并分析了它们的分类和基本特征（如转录本长度、外显子数目通常小于mRNA）。 3. 差异表达谱全面解析： 鉴定出大量与根蘖过程相关的DEGs、DELs和DEMs。功能富集分析将根蘖与光合作用、激素（尤其是细胞分裂素相关）转运与合成通路联系起来，为表型观察（如根蘖多发生于裸露岩石上的根系，可能接触更多光照）提供了分子层面的线索。 4. 关键调控基因与模块识别： 差异表达分析成功定位了两个核心调控基因：BRC1（在RS中上调）和SPL10（在RS中上调）。BRC1是控制芽活性和分枝的关键转录因子，其表达受细胞分裂素正调控，这与研究前期观察到的RS中细胞分裂素水平升高相符。SPL10则是miR156的经典靶标，参与植物幼年性维持和再生过程。研究发现miR156g-3p在RS中表达下调，而SPL10上调，符合其调控关系。 5. 多层次调控网络构建： 通过生物信息学预测，成功构建了一个以BRC1和SPL10为枢纽的miRNA-mRNA-lncRNA互作网络。网络显示，miR319可能靶向负调控BRC1；而SPL10不仅受miR156调控，还与多个lncRNA（如LINC594， LINC7200）存在潜在的相互作用。这些发现串联起了不同层次的调控因子。 6. 核心调控模块提出： 综合所有结果，研究提出了一个名为“miR319-BRC1-miR156-SPL10-lncRNAs”的潜在遗传调控模块，认为该模块协同作用，共同调控香椿根蘖的起始与发育。

结论 本研究通过整合转录组学分析，系统揭示了香椿根蘖克隆繁殖背后的复杂分子调控网络。研究不仅提供了根蘖发育过程的详细表型记录和高通量分子表达谱数据，更重要的是识别出BRC1和SPL10等关键调控基因，并首次在根蘖研究中构建了包含miRNA和lncRNA的多层次调控网络。所提出的“miR319-BRC1-miR156-SPL10-lncRNAs”调控模块为理解植物无性繁殖的遗传机制提供了新的理论框架。该研究加深了我们对香椿这一濒危树种繁殖策略的认识，其发现对于评估克隆繁殖对种群遗传多样性和适应潜力的影响具有重要科学价值，也为未来通过分子手段调控植物繁殖方式（例如在林木无性快繁或控制入侵植物扩散方面）提供了潜在的靶点。

研究亮点 1. 研究对象的独特性与重要性： 聚焦于中国特有濒危树种香椿及其关键的根蘖繁殖现象，将基础研究与物种保护需求相结合，具有明确的生态学和保护生物学意义。 2. 研究技术的系统性与前沿性： 采用了涵盖mRNA、lncRNA和miRNA的全转录组测序策略，并进行了整合分析，能够从多维度、多层次揭示调控机制，是当前非编码RNA功能研究的前沿方法。 3. 分析流程的严谨性： 在lncRNA鉴定上采用多工具交叉验证，极大提高了预测准确性；在调控网络构建中，结合了共表达分析和顺式作用预测，增加了网络的可信度。 4. 创新性发现： 首次在香椿根蘖中系统鉴定并分析了lncRNA，并成功构建了包含lncRNA的根蘖调控网络。提出的“miR319-BRC1-miR156-SPL10-lncRNAs”调控模块是一个新颖的整合性假设，将已知的激素响应、miRNA调控和lncRNA潜在功能联系起来。 5. 扎实的实验验证： 通过RT-qPCR对关键候选分子进行了实验验证，确保了生物信息学预测结果的可靠性，形成了完整的“假设生成-验证”研究闭环。

其他有价值内容 研究在讨论部分也坦诚指出了局限性：样品来自天然林而非可控环境，采样基于形态而非精确的细胞标记；香椿遗传转化体系尚未建立，因此未能进行基因功能的体内验证。这些为未来研究指明了方向。此外，研究联系到全球变暖背景下，与香椿近缘的澳洲红椿（亦具根蘖能力）种群衰退的现象，警示长期依赖克隆繁殖可能削弱物种的气候变化适应力，从而将分子机制研究与宏观生态学问题关联起来，提升了研究的格局。最后，研究提供的详细方法学流程和数据（如lncRNA鉴定步骤、引物设计考虑同源基因问题等）对其他从事非模式植物非编码RNA研究的科研人员具有很高的参考价值。