## 中文标题 北美短叶松（*Pinus banksiana*）铜相关甲基化组变异：与铜抗性关联的DNA甲基化特征

《Plant Stress》：Copper-Associated Methylome Variation in Pinus banksiana: DNA Methylation Signatures Linked to Copper Resistance

　　

## 摘要专业翻译

铜毒性是一种重要的环境胁迫因子，可能与植物DNA甲基化和生理变化相关。本研究利用简并代表性亚硫酸氢盐测序（reduced-representation bisulfite sequencing，RRBS）表征了经表型筛选的铜抗性北美短叶松（*Pinus banksiana*）幼苗的甲基化模式。经铜处理的抗性幼苗表现出独特的甲基化模式，主要特征为向更高甲基化水平区间偏移。在注释基因区域及其侧翼序列之间观察到相反的甲基化模式：基因相关区域甲基化升高，而相邻侧翼区域甲基化降低。这些相反趋势表明，铜暴露后与所筛选铜抗性表型相关的DNA甲基化发生了空间重分布。铜处理的抗性样本在CG、CHG和CHH甲基化背景下还表现出更高的共定位程度，表明这些背景下的甲基化变化发生在重叠的基因组区域内。与对照相比，铜处理的抗性幼苗呈现948个高甲基化差异甲基化区域（differentially methylated regions，DMRs）和279个低甲基化DMRs，大多数位点分别集中于约93%和30%的甲基化水平。本分析中鉴定的所有基因相关DMR均位于注释基因边界内。DMRs与候选胁迫响应基因相关，包括一个MATE1外排转运蛋白和一个patatin样磷脂酶A酰基水解酶（pPLA）。总体而言，本研究鉴定了铜暴露后经表型筛选的铜抗性*P. banksiana*幼苗相关的DNA甲基化特征，并为未来功能研究提供了候选基因组区域。

## 论文解读文章

### 一、研究背景与目的

铜是植物必需微量元素，参与光合作用、营养转运、细胞呼吸和抗氧化防御等多种生理过程，但过量铜会通过Haber-Weiss和Fenton反应产生ROS，造成细胞损伤并诱导DNA甲基化模式改变。已有少数研究报道铜对植物全基因组甲基化的影响，但针对松科针叶树的研究极为有限。北美短叶松（*Pinus banksiana*）具有铜抗性，且在所有已报道物种中具有最高水平的基线甲基化之一，其基因组庞大并富含重复序列和转座元件。此前对该物种的转录组分析已发现铜抗性幼苗中金属转运和螯合相关基因广泛上调，但铜诱导的全基因组甲基化特征仍不清楚。本研究首次利用RRBS技术对铜处理的抗性*P. banksiana*进行全基因组甲基化分析，旨在评估过量铜对全基因组甲基化的影响，并鉴定与铜暴露相关的差异甲基化区域（DMRs）及其关联基因。

### 二、研究内容与技术方法

研究人员从加拿大安大略省萨德伯里VALE植物生产设施获取六个月苗龄的*P. banksiana*幼苗，在受控环境生长室中以1,300 mg CuSO₄·5H₂O kg^-1土壤（10 mmol Cu）处理14天，并设置去离子水阴性对照和等摩尔浓度K₂SO₄硫酸根对照。依据1–9级损伤评分进行表型筛选，选取损伤评分为1、无可见症状的三株最强抗性幼苗作为独立生物学重复。提取针叶总DNA，经MspI酶切构建RRBS文库，在Illumina HiSeq X Ten平台测序。测序数据经TrimGalore和FastQC质控后，使用Bismark软件比对到*Pinus taeda*参考基因组（ptaeda.v1.0），计算CG、CHG和CHH三种序列背景的甲基化水平。采用Metilene软件基于二元分割算法识别DMR，以q值≤0.05且平均甲基化差异≥10%为阈值；通过Pearson相关、层级聚类和主成分分析（PCA）评估样本间甲基化谱关系；利用Gene3D、PANTHER、InterPro等数据库对DMR关联基因进行功能注释。

### 三、主要研究结果

**3.1 RRBS测序质量与全基因组甲基化胞嘧啶比例**：测序深度和覆盖率均较高，MspI消化效率超过95%，亚硫酸氢盐转化效率超过99%。铜处理样本总甲基化胞嘧啶比例最高，主要由mCHG和mCG增加驱动，mCHH比例降低。

**3.2 铜诱导的平均全基因组甲基化水平变化**：虽然各处理组间平均甲基化水平差异较小，但铜是唯一在所有序列背景下均提高平均全基因组甲基化水平的处理；钾处理未显示此效应。

**3.3 甲基化水平跨序列背景的再分布**：铜处理使CG位点从未甲基化（0%）和完全甲基化（100%）状态向中间甲基化区间（25–99%）移动；而CHG和CHH甲基化增加主要来自完全甲基化位点的富集，说明铜促进非CG背景的二元化甲基化。

**3.4 DNA甲基化在基因相关区域的空间分布**：铜是唯一增加基因体甲基化占比的处理，同时降低基因侧翼区域的甲基化比例；下游侧翼区域甲基化水平普遍高于上游，表明铜诱导的甲基化重编程具有空间组织性。

**3.5 金属胁迫下CG、CHG和CHH甲基化的共定位**：铜处理样本在mCG与mCHG、mCHH之间表现出最高的共定位程度，且共享甲基化位点数目增加，而钾处理与水对照相似，说明铜特异性地增强不同序列背景甲基化的空间重叠。

**3.6 全基因组CpG甲基化谱与样本聚类**：Pearson相关显示所有样本间全基因组CpG甲基化整体相似性高；层级聚类中部分铜样本与水对照聚集，但PCA可将铜处理样本与对照分离。其中铜处理样本C7形成独立簇，提示抗性个体间存在甲基化异质性。

**3.7 差异甲基化区域的鉴定与特征**：基于CG背景数据，铜与水对照比较共鉴定出1,227个DMR（其中高甲基化948个，低甲基化279个），平均长度约160 bp；铜处理样本中较短DMR频率更高。

**3.8 差异甲基化区域的基因组分布**：各处理比较中，基因相关DMR主要位于基因体内，其次为内含子；铜与水比较中侧翼区域DMR数多于外显子、启动子和转录起始位点，而铜与钾比较中基因外显子DMR比例更高。

**3.9 差异甲基化区域的甲基化谱**：铜样本中高甲基化DMR的甲基化水平集中于93%附近，而对照对应位点约为30%；低甲基化DMR则相反，铜样本约30%，对照约93%，表明铜暴露同时产生局部强烈高甲基化和位点特异性低甲基化。

**3.10 共有与处理特异的差异甲基化区域**：铜和钾处理与水对照的比较中，共有DMR仅占总数的15.81%–17.37%，大部分为基因间区高甲基化DMR，提示两种金属存在部分共同的表观遗传响应靶点，但多数响应为处理特异。

**3.11 铜响应DMRs的功能注释**：基因相关DMR分布于14个染色体支架上，所有支架均含基因体DMR，多数含内含子DMR；启动子DMR较少，转录起始位点未检出DMR。候选基因包括MATE1外排转运蛋白、patatin样磷脂酶A酰基水解酶（pPLA）、WRKY1转录因子、DEAD-box ATP依赖性RNA解旋酶、丝氨酸/苏氨酸蛋白激酶等胁迫响应相关基因。

### 四、讨论总结与结论

研究讨论指出，铜处理与基因体甲基化升高、侧翼区域甲基化降低、CG/CHG/CHH甲基化共定位增强以及转座子富集区域差异甲基化相关，这些模式可能构成铜抗性表型相关的表观遗传特征。但研究也强调，由于未整合转录组数据且仅选取抗性表型，甲基化变化与基因表达或抗性机制之间的因果关系尚不能确立。钾处理作为硫酸根对照也表现出部分甲基化效应，提示铜响应中可能包含与钾相关的共同生理过程。结论部分可概括为：本研究鉴定出与铜处理的抗性*P. banksiana*相关的DNA甲基化特征，支持铜暴露与基因体甲基化增加、侧翼调控区域甲基化减少、CG–CHG–CHH甲基化协同动态及转座子富集区域差异甲基化关联的模型。该工作首次提供了针叶树铜抗性相关的全基因组甲基化图谱，为金属污染环境下植物表观遗传适应机制研究奠定了基础。