油棕（Elaeis guineensis Jacq.）幼苗冷胁迫生理损伤时序特征与冷响应miRNA表达谱解析

《Frontiers in Plant Science》：Physiological characterization and cold-responsive miRNA profiling of oil palm (Elaeis guineensis Jacq.) under cold stress

　　

引言：冷胁迫严重制约中国油棕产业发展。微小RNA（miRNA）是一类小非编码RNA（small non-coding RNAs），作为基因表达的关键转录后调控因子，在植物生长发育和逆境适应中发挥重要作用。尽管miRNA已在多种植物非生物逆境中被广泛表征，但其在油棕（Elaeis guineensis Jacq.）冷胁迫响应中的调控功能仍缺乏系统研究。方法：研究人员以薄壳Tenera油棕基因型幼苗为材料，将苗龄一致的幼苗置于8 °C下进行冷处理，并于0（26 °C）、0.5、1、2、4和8 h采样，每个时间点设三个生物学重复。研究人员评估了丙二醛（MDA）、电解质渗漏和叶绿素含量等生理损伤指标；在Illumina HiSeq Xten平台上对18个独立小RNA文库进行测序，鉴定冷响应miRNA及其预测靶基因；并通过定量PCR（qPCR）进一步验证17个候选miRNA的表达。结果与讨论：生理测定显示，随着冷胁迫加剧，叶片中MDA积累和电解质渗漏率显著升高，叶绿素含量持续下降，表明细胞膜完整性逐步受损。为此，研究人员构建了18个独立小RNA测序文库，用于小RNA表达谱分析，包括一个常温对照（26 °C）和五个冷处理时间点（8 °C下0.5、1、2、4和8 h），每个实验组三个生物学重复。测序初步鉴定出191个未过滤的已知miRNA候选和98个预测的新miRNA候选。功能富集分析显示，差异表达miRNA的靶基因主要富集于转录调控通路以及倍半萜和三萜生物合成。qPCR确认了17个冷响应miRNA的表达动态，其中基于表达相关性观察到miR166家族与其预测靶基因ATHB9-like和ATHB12-like之间的负调控关系。研究数据支持以下推断：miRNA依赖的转录后重编程参与油棕冷胁迫适应，而非解析完整分子机制。进一步分析表明，冷响应miRNA的时序表达变化与植物形态、生理和生化指标的动态变化一致。对五个冷处理比较的交集分析定义了保守和新miRNA的核心冷响应集合与阶段特异性集合；无偏Mfuzz时间序列聚类将所有冷响应miRNA分为六个时序表达模块。

研究背景与问题：油棕（Elaeis guineensis Jacq.）是多年生热带油料作物，主要分布于南纬10°至北纬15°，在中国主要种植于海南、云南等热带地区。油棕适于温暖湿润气候，冬季低温可导致花败育和产量下降，是限制其生长发育和棕榈油生产的关键因素。已有研究报道油棕可通过形态变化、渗透调节、抗氧化酶系统和激素信号通路适应冷胁迫，SSR标记和转录组分析也鉴定了部分冷响应基因，但其冷响应调控机制仍待深入。miRNA作为转录后调控因子，可通过切割冷响应基因或抑制翻译参与植物耐冷性调控，且许多miRNA靶基因编码转录因子。尽管部分保守miRNA及其靶基因已在油棕中被鉴定，并涉及盐胁迫和花器官发育，但油棕冷胁迫下miRNA表达谱的系统研究仍有限。为此，研究人员以薄壳Tenera油棕幼苗为对象，结合生理指标、小RNA测序、靶基因预测和qPCR验证，解析冷胁迫下miRNA表达变化，旨在提供优先候选miRNA（包括Eg-miR166、Eg-miR393、Eg-miR156、Eg-miR444）供后续功能验证和耐冷分子育种，并为理解该热带作物miRNA相关冷胁迫响应提供转录水平基础资料。

研究人员开展的研究、结论与意义：研究人员对8 °C冷处理0、0.5、1、2、4和8 h的油棕幼苗进行破坏性单时点采样，测定MDA、电解质渗漏、叶绿素和CIELAB颜色指标，并构建18个独立小RNA文库进行测序，鉴定保守和新miRNA，预测靶基因，进行GO和KEGG富集、Upset交集分析和Mfuzz时序聚类，并以qPCR验证17个冷响应miRNA。研究结论是：冷胁迫引起渐进性膜脂过氧化和叶绿素降解；miRNA表达谱在冷暴露后快速重塑；冷响应miRNA包括核心集合和阶段特异性集合，并可分为六个时序模块；靶基因显著富集于转录调控及倍半萜和三萜生物合成等通路；miR166家族与ATHB9-like、ATHB12-like转录本呈负相关。该研究的意义在于提供油棕冷响应miRNA的时序框架和候选基因资源，但原文明确指出其不能完全解析油棕耐冷分子机制，靶基因调控关系仍需降解组测序或5′-RLM-RACE等验证。

关键技术方法概括（不超过250字）：研究人员以薄壳Tenera（Dura × Pisifera F₁）油棕幼苗为材料，幼苗来自中国热带农业科学院椰子研究所种质苗圃（海南文昌，19.63°N，110.75°E）；选取10月龄、大小一致且具8片完全展开功能叶的幼苗，18株用于破坏性单时点采样，常温对照3株，8 °C冷处理0.5、1、2、4、8 h各3株，取枪叶。关键技术包括：测定MDA、电解质渗漏、叶绿素及CIELAB颜色指标；构建18个独立小RNA文库并用Illumina HiSeq Xten测序；采用miRDeep-P2等预测新miRNA，psRNATarget、TargetFinder和miRanda预测靶基因；qPCR验证候选miRNA及靶转录本；进行GO/KEGG富集、Upset交集分析和Mfuzz时序聚类。

研究结果：（1）冷胁迫下油棕叶片叶绿素、MDA和电导率变化：通过8 °C时间梯度处理结合生理和颜色指标测定，发现叶片由轻微黄化发展至严重褐变，MDA和相对电导率随时间显著升高，叶绿素含量下降，总色差ΔE_ab与MDA和电导率正相关、与总叶绿素负相关，表明冷胁迫导致渐进性细胞膜损伤，颜色指标可定量反映冷害程度。（2）油棕小RNA文库构建与下一代测序：18个独立文库测序获得高质量clean reads，sRNA长度主要分布于21—25 nt；相关性和主成分分析显示常温对照与冷处理样本明显分离，冷处理0.5 h后全局miRNA谱已快速重塑。（3）油棕冷响应sRNA鉴定与特征：通过与miRBase等数据库比对注释sRNA，发现冷胁迫下24-nt sRNA比例总体下降，21-nt sRNA比例先降后升，表明冷胁迫扰动异染色质siRNA稳态并重塑21-nt转录后调控小RNA库。（4）冷胁迫下油棕保守miRNA鉴定：鉴定到191个未过滤已知成熟miRNA候选，并报告269个保守miRNA条目（含isoform/variant）及6,179个靶基因；严格过滤去除24-nt siRNA后，六组共鉴定5,032个已知植物来源成熟miRNA，286个为六组共有；长度以21/22 nt为主，具5′尿苷偏好。（5）冷胁迫下油棕新miRNA鉴定：基于参考基因组和miRDeep-P2等预测，六组共发现95个新miRNA序列，21-nt为主导长度类；预测获得7,978个对应新miRNA的靶基因，提示物种特异miRNA参与冷响应调控网络。（6）冷条件下油棕miRNA表达变化谱：共鉴定48个已知成熟miRNA和11个新miRNA差异表达；差异表达miRNA数量随处理时间先增后减，早期至中期以上调为主；Upset分析定义核心保守冷响应集合（包括上调的miR167-x、miR162-x、miR11108-x、miR3627-z、miR393-z、miR444-y、miR7767-x和下调的miR157-y、miR1863-z）及核心新miRNA集合，并鉴定阶段特异性成员。（7）冷响应miRNA的GO和KEGG通路分析：靶基因GO注释覆盖生物过程、细胞组分和分子功能；KEGG富集显示倍半萜和三萜生物合成通路在五个冷处理比较中均显著富集，8 h富集最强；早期涉及抗氧化和防御相关代谢，中期涉及膜系统与能量代谢，晚期涉及DNA修复；HMGCR被miR159和miR396靶向，SQS被miR172靶向，β-AS被miR396靶向。（8）冷条件下油棕miRNA表达谱：qPCR验证17个冷响应miRNA，分为高于对照和低于对照两组；Eg-miR393呈线性上升，Eg-miR444早期升高后下降，Eg-miR166a、Eg-miR166-x、Eg-miR391-x、Eg-miR395-z和Eg-miR398-y低于对照且8 h最低；miR166成员与ATHB9-like、ATHB12-like转录本呈负相关；qPCR与small RNA-seq的log₂(fold-change)高度相关，R²为0.9441—0.9960。（9）油棕响应冷胁迫的阶段特异性miRNA：结合qPCR表达轨迹将冷响应miRNA归为早期瞬时、中期、晚期持续、持续上调和持续下调等模式；Mfuzz对317个冷响应miRNA（266个新miRNA和51个保守miRNA）进行无偏聚类，解析出M1—M6六个时序模块，M1为0.5 h早期瞬时上调模块，包含miR6255-z和miR168-z等。

讨论部分总结：讨论将冷害进程概括为早期感知（0.5 h）、膜损伤起始（1—2 h）和严重冷害/细胞损伤（4—8 h），并认为2 h为生理转折点；CIELAB总色差可作为叶绿素降解和膜氧化损伤的高通量代理。小RNA测序显示冷胁迫快速重塑sRNA和miRNA谱，24-nt异染色质siRNA比例下降而21-nt调控性小RNA动态变化。保守和新miRNA共同构成冷响应网络，新miRNA数量随冷胁迫持续增加，提示物种特异调控成分；21/22 nt长度和5′尿苷偏倚符合植物miRNA特征。miR393、miR166-HD-ZIP模块和miR156等与其他单子叶植物冷响应具有部分一致性，且miRNA家族成员存在表达分化。靶基因富集提示倍半萜和三萜生物合成参与冷响应，但原文将其表述为可检验工作假设，因缺乏降解组和代谢组直接证据；HMGCR、SQS和β-AS为预测靶基因。整合模型认为miRNA依赖的转录后重编程在可见冷害前已快速启动，随后与细胞损伤信号叠加。局限性包括：靶基因均为计算预测，缺乏降解组测序或5′-RLM-RACE切割证据；仅使用单一冷敏感Tenera基因型，不能区分耐冷调控机制与一般冷害响应；种子繁殖F₁存在残余遗传异质性，结合破坏性单时点采样使个体遗传差异与处理效应混杂；缺少与冷处理采样时间匹配的常温对照，不能完全排除昼夜节律影响；4—8 h连续光照可能加剧冷害褐变。

研究结论翻译：本研究表征了油棕冷胁迫下差异表达miRNA及其靶基因。冷响应miRNA的靶基因显著富集于转录调控和次级代谢通路，包括倍半萜和三萜生物合成。17个冷响应miRNA经qPCR验证并按表达模式分组。相关性表达分析表明，这些miRNA的时序变化与冷暴露过程中发生的生理生化改变（膜脂过氧化、叶绿素降解和形态损伤）相一致。尽管该工作不能完全解析油棕耐冷的完整分子调控机制，但识别出优先候选miRNA（如Eg-miR166、Eg-miR393、Eg-miR156、Eg-miR444）用于功能验证和耐冷油棕材料的分子育种，并为未来热带油棕幼苗miRNA介导冷胁迫响应的机制研究提供转录水平基础资源。全基因组Upset交集分析和无偏Mfuzz时间序列聚类进一步定义了核心冷响应miRNA集合和六个时序表达模块，为油棕冷胁迫响应提供了全基因组时序框架。