《Epigenomes》：Cold Stress-Induced DNA Methylation Profiling and Its Association with Terpenoid Biosynthesis in Germinating Mikania micrantha Seeds
编辑推荐：
摘要专业翻译
薇甘菊（Mikania micrantha）是一种高度入侵性杂草，在中国华南地区，由于冬季低温的限制，其向北扩展的范围受到制约。尽管DNA甲基化日益被认为是植物逆境适应的介导因子，但在入侵植物中，表观遗传调控是否将耐冷性与化感次生代谢物的产生相偶
摘要专业翻译
薇甘菊（Mikania micrantha）是一种高度入侵性杂草，在中国华南地区，由于冬季低温的限制，其向北扩展的范围受到制约。尽管DNA甲基化日益被认为是植物逆境适应的介导因子，但在入侵植物中，表观遗传调控是否将耐冷性与化感次生代谢物的产生相偶联仍属未知。本研究通过整合全基因组重亚硫酸盐测序（whole-genome bisulfite sequencing, WGBS）与RNA测序（RNA-Seq），在三种温度条件（14 °C、25 °C和30 °C）下全面表征了萌发胚中的全基因组DNA甲基化动态及其转录后果。暴露于14 °C显著延迟了萌发进程但并未使其停滞。全局甲基化水平在CG位点最高（80.21%–84.63%），其次为CHG（63.67%–66.77%）和CHH（6.63%–7.51%）序列背景。整合分析鉴定出1103个差异甲基化区域（differentially methylated region, DMR）相关基因，其中434个同时差异表达。值得注意的是，对这434个重叠基因的KEGG富集分析揭示了倍半萜和三萜生物合成途径的显著过度表达，且甲羟戊酸（mevalonate, MVA）途径的关键基因在14 °C下表现出与DMR相关的表达改变。由于这是一项整合性甲基化组–转录组关联研究，未进行代谢物定量或功能验证，因此这些结果与以下假说一致：冷诱导的表观遗传重编程与萜类生物合成的转录变化相关联，这可能有助于低温胁迫下的抗氧化和膜保护功能；化感性倍半萜内酯的可能积累及其在冷季维持竞争优势中的任何作用仍属待实验验证的工作假说。本研究结果为薇甘菊入侵的管理与防控提供了新的表观遗传学视角和支持数据。
论文解读
一、研究背景与科学问题
薇甘菊（Mikania micrantha）以其极快的繁殖速度著称，常被形容为“一分钟一英里”杂草，这一特性主要归因于其广泛的生理耐受性、旺盛的无性繁殖能力、化感作用以及大量优质种子的产生。在中国南方，薇甘菊的花期和瘦果产生期从8月持续至次年5月，而幼苗建立的高峰期通常为9月至11月，与约25–30 °C的最适萌发温度相吻合。尽管有报道称薇甘菊在15 °C或35 °C下也能萌发，但明显偏离最适温度的条件会显著抑制其萌发和幼苗建立。在中国南方薇甘菊的入侵区域，冬季温度经常降至15 °C以下。结合中国境内分布记录，该物种主要分布于广东、广西、湖南和福建等省份，最北端可达福州（福建）且有进一步向北扩展的趋势。这些证据表明，低温生长环境是限制其在中国向北扩散的主要制约因素之一。然而，关于低温胁迫对薇甘菊种子发育影响的研究仍然匮乏，尤其是在萌发过程中分子遗传变异与表观遗传调控的整合分析方面。
低温是限制植物生长、生产力及地理分布的主要环境因素。植物已进化出精密的生理和分子适应机制来应对冷胁迫，包括CBF/DREB1转录级联、相容性渗透调节物质的积累以及膜脂组成的调节。在这些已充分表征的遗传途径之外，包括DNA甲基化、组蛋白修饰、染色质重塑和非编码RNA在内的表观遗传机制在调控胁迫响应基因表达中发挥关键作用，并能提供快速、可逆且可能可遗传的适应性响应。冷诱导的DNA甲基化变化已在多种植物中被记录。表观遗传在生物入侵中的作用日益受到关注。入侵植物物种往往在奠基者效应导致的低遗传多样性条件下仍能建立和扩散，提示非遗传来源的表型变异可能尤为重要。在薇甘菊中，群体水平研究揭示了与环境异质性相关的显著表观遗传分化，且前期研究已证明冷依赖性DNA去甲基化与抗氧化调控共同促进低温下种子萌发的恢复。然而，冷诱导表观遗传重编程在该物种中的功能靶标，特别是甲基化变化是否调控与胁迫耐受性和入侵性均相关的次生代谢途径，仍基本未被探索。
冷胁迫诱导活性氧（reactive oxygen species, ROS）的过量产生，导致脂质、蛋白质和核酸的氧化损伤。除酶促抗氧化系统外，植物还积累多种多样的专化（次生）代谢物作为非酶抗氧化剂和膜稳定剂。其中，萜类化合物作为植物专化代谢物中最大且结构最多样化的类别，在非生物胁迫耐受中发挥突出作用。挥发性异戊二烯类化合物，包括异戊二烯、单萜和倍半萜，能够直接清除ROS、淬灭脂质过氧化链式反应，并在温度胁迫下稳定类囊体膜。冷胁迫还可上调檀香中萜类生物合成基因，并在番茄中按胁迫严重程度比例定量增强单萜和倍半萜的释放。类黄酮（包括花青素）构成另一大类冷响应代谢物，薇甘菊本身也会积累花青素作为可见的冬季适应特征。重要的是，次生代谢的表观遗传调控已有充分文献记载：DNA甲基化控制黄花蒿中青蒿素（倍半萜）的生物合成、长春花中单萜吲哚生物碱的产生以及丹参中丹参酮和酚酸的积累。这些发现确立了DNA甲基化可作为专化代谢途径的调控开关，但在入侵植物冷胁迫过程中是否存在此类表观遗传调控，以及其是否靶向萜类、类黄酮或生物碱途径，此前尚未被研究。
尽管取得了上述进展，尚无整合性甲基化组–转录组研究考察冷诱导的表观遗传变化是否重编程入侵植物的次生代谢，也未解决此类代谢重编程在胁迫耐受与化感作用中的潜在双重功能。基于前期关于薇甘菊冷胁迫萌发伴随DNA甲基化动态和抗氧化调控的发现，本研究采用全基因组重亚硫酸盐测序和RNA-Seq技术，在14 °C、25 °C和30 °C下对萌发胚进行分析，旨在：（i）表征CG、CHG和CHH三种序列背景下全基因组甲基化变化；（ii）鉴定DMR相关的差异表达基因及其功能富集；（iii）确定次生代谢途径（特别是萜类生物合成）是否在冷胁迫下受表观遗传调控。
二、研究内容与主要结论
研究人员在三种恒定温度（14 °C、25 °C和30 °C）下进行种子萌发试验。吸水曲线显示，30 °C下种子在0–18 h内快速吸胀，48–72 h间40%种子破种皮成功萌发；25 °C下呈现相似模式，0–24 h快速吸水，48–72 h间40%萌发。因此选择72 h（第3天）作为30 °C和25 °C的萌发里程碑时间点。14 °C下萌发显著延迟，吸胀至第5天才完成，平台期从第5天延伸至第8天，第9天观察到萌发。根据萌发胚长度测量，14 °C第7天的胚长度与30 °C/25 °C第3天相当，14 °C第10天的胚长度与30 °C/25 °C第7天相当。最终采样方案包含六个时间点：30 °C第3天和第7天；25 °C第3天和第7天；14 °C第7天和第10天。
RNA-Seq共生成约14.1亿条原始读段，质控后保留96.82%，获得约13.6亿条干净读段，平均每个文库7580万条。Q20和Q30值分别落在97.86%–98.39%和93.51%–94.90%范围内。配对末端干净读段比对至薇甘菊参考基因组，获得40,088个注释基因，平均比对率为90.31%。外显子碱基占比对碱基的84.0%–89.1%。生物学重复间的Pearson相关系数超过0.8，证实了表达数据的可靠性。差异表达分析共进行九对比较。30 °C下第3天与第7天之间鉴定出7956个差异表达基因（differentially expressed genes, DEGs）；25 °C对应比较获得2274个DEGs；14 °C第7天与第10天之间检测到8679个DEGs，为最大数量。主成分分析显示生物学重复良好聚类，PC1和PC2分别占总方差的27.26%和16.53%。
全基因组重亚硫酸盐测序数据质量控制显示，总检出位点数为1,343,423,094至1,369,305,764个，全基因组位点平均覆盖深度为13.14×至16.76×。≥1×测序深度的碱基占全基因组的87.64%–89.33%，≥5×和≥10×分别占73.42%–77.79%和56.65%–64.56%。所有基因组C位点的整体平均甲基化水平为23.55%–25.54%。三种序列背景中，CG背景平均甲基化水平最高（80.21%–84.63%），其次为CHG（63.67%–66.77%），CHH背景最低（6.63%–7.51%）。整合分析显示，14 °C萌发胚中基因侧翼区域CHH背景甲基化水平升高，而25 °C萌发胚中基因体CG和CHG背景甲基化水平略高。25 °C与14 °C比较中，共鉴定出1757个DEGs和1103个DMR相关基因，其中434个基因重叠。KEGG功能富集分析显示，这434个重叠基因中倍半萜和三萜生物合成途径显著富集。对MVA途径关键基因（MVK、HMGCS、MVD和FDPS）的整合分析显示，HMGCS在14 °C下FPKM明显更高，MVK和FDPS较低，MVD基本不变。
三、关键技术方法
本研究采用的主要技术方法包括：全基因组重亚硫酸盐测序（WGBS）用于全基因组DNA甲基化谱分析，使用Bismark进行比对和甲基化胞嘧啶检出，DSS用于差异甲基化区域鉴定；RNA-Seq转录组测序用于基因表达定量，使用HISAT2比对参考基因组，featureCounts获得读段计数并转换为FPKM值，DESeq2用于差异表达基因鉴定；整合甲基化组–转录组分析，通过维恩图取DMR相关基因与DEGs的交集，clusterProfiler进行KEGG途径富集分析。植物材料采集自中国广东省深圳市内伶仃岛。种子经表面消毒后，在MS培养基上于14 °C、25 °C和30 °C三种温度条件下萌发。
四、研究结果
2.1 种子萌发预实验、RNA提取与质量评估。 通过吸水曲线确定了三种温度下的萌发阶段划分和采样时间点。RNA质量检测显示所有样品RIN值大于6.0，满足cDNA文库构建要求。三种温度条件下MS培养基上培育的幼苗表现出明显的形态差异。
2.2 冷胁迫下萌发胚的转录重编程。 通过RNA-Seq技术生成约14.1亿条原始读段，质控后获得约13.6亿条干净读段，平均比对率为90.31%，确认了转录组数据的可靠性。外显子碱基占比对碱基的84.0%–89.1%。
2.3 基因表达定量与样本相关性分析。 FPKM密度分布显示三种温度条件下表达模式总体可比，14 °C处理的基因整体表达水平略高于其他两种温度。生物学重复间Pearson相关系数超过0.8，证实了表达数据的可靠性。
2.4 差异基因表达分析与功能富集。 九对比较共鉴定出大量DEGs，其中14 °C第7天与第10天之间数量最多（8679个）。PCA分析显示生物学重复良好聚类。
2.5 重亚硫酸盐测序数据质量控制。 全基因组平均覆盖深度为13.14×至16.76×，至少56.65%的胞嘧啶位点达到≥10×测序深度。CG背景甲基化水平最高（80.21%–84.63%），CHH背景最低（6.63%–7.51%）。
2.6 表观遗传与转录组整合分析。 14 °C萌发胚中基因侧翼区域CHH背景甲基化水平升高，25 °C萌发胚中基因体CG和CHG背景甲基化水平略高。基因体甲基化与基因表达总体呈正相关。
2.7 DMR相关基因功能注释分析。 25 °C与14 °C比较中鉴定出434个DEGs与DMR相关基因的重叠基因，KEGG富集显示倍半萜和三萜生物合成途径显著富集。MVA途径中HMGCS在14 °C下表达明显更高。
五、讨论与结论
讨论部分指出，环境因素已知可调节植物的表观遗传性状和适应性进化过程。胞嘧啶甲基化（5-mC）是与程序化发育和胁迫响应相关的表观遗传标记，DNA甲基化已被证明可调节植物对环境挑战的响应。在植物中，DNA甲基化常见于CG、CHG和CHH三种序列背景。物种内产生的甲基化模式相对稳定且可跨代遗传，但在胁迫暴露下，植物基因组DNA甲基化可发生改变，这被认为是应对不利环境的适应机制之一。结合薇甘菊低温萌发胚中观察到的甲基化和表达谱，该物种可能通过增加胁迫响应基因的表达来增强耐冷性。萜类代谢代表了植物抵御冷诱导氧化胁迫的额外防御层。冷胁迫可主动触发植物中萜类积累，这些专化代谢物具有内在抗氧化活性，可直接清除低温条件下产生的过量ROS，从而减轻膜脂过氧化并保持细胞膜完整性。整合甲基化组和转录组分析揭示了倍半萜和三萜生物合成途径的显著富集，HMGCS在14 °C萌发胚中表达明显更高，提示改变的萜类生物合成可能构成表观遗传调控的适应策略。DNA甲基化对基因表达的调控并非孤立发生，而是与组蛋白修饰、染色质重塑和非编码RNA等多种表观遗传机制协调运作。薇甘菊作为全球性入侵杂草，其显著的环境适应性可能部分归因于表观遗传调控的灵活性和可塑性。研究也承认了两项局限性：采样虽按胚长度匹配以比较发育等效阶段，但无法完全分离温度效应与时间年龄依赖的表观遗传变化；本研究为相关性研究，缺乏功能验证，关于萜类生物合成表观遗传调控的结论仍属初步，需通过qRT-PCR验证DEGs、位点特异性重亚硫酸盐验证DMRs、GC-MS/LC-MS萜类定量及化感生物测定加以确认。
研究结论部分指出，对14 °C、25 °C和30 °C下萌发薇甘菊胚的整合甲基化组和转录组分析表明，冷胁迫延迟但不阻止萌发，并伴随序列背景和位置依赖的DNA甲基化重编程。冷萌发条件下基因上游2 kb（启动子）和下游侧翼区域CHH甲基化显著升高，而基因体甲基化（尤其是CG和CHG背景）在25 °C下略高，且与基因表达总体呈正相关。在434个DMR相关的差异表达基因中，倍半萜和三萜生物合成途径显著富集，MVA途径关键基因（AACT、HMGCS、MVD、MVK、FDPS）在14 °C下表现出与DMR相关的差异表达。这些发现与以下假说一致：冷相关的表观遗传变化与萜类生物合成的转录调控相关联，可能通过抗氧化和膜保护性萜类有助于耐冷性；化感性倍半萜内酯的可能产生及其生态后果仍属未验证的工作假说。这种通过表观遗传调控将胁迫耐受性与入侵化学联系起来的潜在机制，可能代表了薇甘菊在亚热带和暖温带边缘生境定殖的重要代谢基础，并为气候变化下入侵植物管理提供了新视角，但因果性表观遗传调控的建立仍需功能验证。
