猕猴桃属（Actinidia）不同倍性与遗传背景下的全基因组DNA甲基化模式及基因表达差异：多倍体化后表观调控独立于剂量效应

《Current Plant Biology》：Whole-genome methylation patterns and gene-expression differences across ploidy and genetic backgrounds in Actinidia species

　　

在六倍体 A. deliciosa 中，大多数基因相对于二倍体呈下调表达；上调基因富集于激素生物合成（hormone biosynthesis）和胁迫响应（stress-response）通路；DNA甲基化（DNA methylation）随倍性水平发生动态重塑；多倍体化后的表观遗传调控（epigenetic regulation）独立于剂量效应（dosage effects）。

研究背景与问题：全基因组复制（whole-genome duplication, WGD）在植物演化中广泛存在，可提高基因组可塑性并促进生态与功能分化，但同时改变基因剂量，扰动剂量敏感调控因子和分子复合体的表达化学计量。DNA甲基化（DNA methylation）被认为参与多倍体化后重复基因组稳定和表达调控；RNA导向DNA甲基化（RNA-directed DNA methylation, RdDM）在部分多倍体系统中与恢复沉默有关，但该研究所分析的 Actinidia 材料未测试RdDM。已有盐响应甲基化等实例为检查 Actinidia 中胁迫相关通路注释提供理由，但不能证明所采样材料发生胁迫适应。转录-甲基化耦合具有物种和位点特异性，需按位点和通路解释而非统一全基因组解释。剂量敏感基因倾向于保留以维持平衡，剂量不敏感的单功能基因更易丢失；多倍体化早期可通过RdDM等暂时限制特定重复以缓解剂量失衡，长期可发生亚功能化或新功能化。尽管 Actinidia 基因组资源已阐明基因组演化与结构，并已编目 A. chinensis 中甲基化/去甲基化机器，但整合 CpG、CHG、CHH 和候选6 mA 的跨倍性与遗传背景全基因组甲基化比较仍有限。因此，研究人员开展跨材料甲基组比较，目的为描述甲基化差异、识别差异甲基化区域（differentially methylated regions, DMRs）相关候选基因和通路，并检查六倍体 A. deliciosa 中甲基化-表达关联；研究将TPM（transcripts per million，转录本丰度指标）作为转录本丰度，而非基因剂量或剂量补偿的直接测量。该研究为多倍体表观调控与剂量效应研究提供候选位点和通路，但因设计限制不能建立因果或单调倍性效应。

关键技术方法：研究人员对4份材料叶片进行Oxford Nanopore Technologies（ONT）长读长测序，包括湖北野生六倍体 A. deliciosa JZP036、野生四倍体 A. deliciosa TshT4、野生二倍体 A. chinensis TshD2及武汉植物园栽培二倍体 A. chinensis 'Hongyang' HyaD2，各1株。Guppy碱基识别、minimap2比对；Tombo检测CHH/CHG/候选6 mA，nanopolish检测CpG；SMART2 DeNovoDMR检测DMR；RNA-seq覆盖JZP036茎、叶、花、根、果，HISAT2/featureCounts计算TPM；果实跨研究比较使用已发表' Hongyang' 80 d after pollination数据（3个生物重复；Illumina HiSeq 2000；Red5_PS1_1.69.0；CRA003106）。

研究结果：2.1 果实基因表达的探索性跨研究差异（Exploratory cross-study differences in fruit gene expression）：通过比较六倍体 A. deliciosa 与二倍体 A. chinensis 'Hongyang' 果实TPM，发现60.3%的基因对在 A. deliciosa 中较低（log₂ratio<-1），5.8%较高（log₂ratio>1）；单等位和单拷贝基因高表达比例高于六等位和六拷贝基因；最高TPM比（log₂ratio>4）基因涉及35个GO（Gene Ontology，基因本体）生物学过程术语和11个KEGG（Kyoto Encyclopedia of Genes and Genomes，京都基因与基因组百科全书）通路，包括激素相关代谢与生物合成、维生素E等脂溶性维生素生物合成及对真菌、细菌、盐、低光响应。结论为这些注释模式提供候选功能差异，但受物种、遗传背景和实验处理差异限制。

2.2 跨 Actinidia 材料的全基因组甲基化模式（Whole-genome methylation patterns across Actinidia materials）：通过ONT甲基组测序和材料间比较，获得各材料候选甲基化位点，超过60%位点处于6 mA上下文；CHH和CHG仅在174条染色体中的30条（29个染色体组中的9个）检测到，CpG和6 mA覆盖全部染色体；四材料平均甲基化高度相关（r>0.85），TshT4与HyaD2聚类且6 mA较高；JZP036中63.20%基因无四类甲基化，2.79%具四类，六等位基因未甲基化比例高于单等位或单拷贝基因，表达-甲基化相关弱；JZP036和TshD2基因体及上下游2 kb的候选6 mA和CHH较低，HyaD2和TshT4较高；未观察到跨采样倍性水平的单调甲基化模式。DMR比较显示JZP036与TshD2比较DMR最少，多倍体比较共享更多低甲基化DMR尤其6 mA；多数DMR位于启动子和远端基因间区，重叠487个基因的674个外显子，关联基因多涉及转运和代谢，38个为WRKY、bHLH、MADS-box等转录因子（transcription factors, TF）；GO/KEGG注释包括光合电子传递、除草剂响应、有氧电子传递链、糖蛋白生物合成与糖基化、花生四烯酸代谢、角质/软木脂/蜡生物合成、二萜生物合成、托烷/哌啶/吡啶生物碱生物合成、酮体合成与降解、油菜素类固醇生物合成和光合天线蛋白。结论为DMR相关候选基因和通路提供描述性证据，但未进行区域富集正式检验，甲基化水平与TPM无清晰关联。

2.3 甲基化水平与基因表达的两两相关（Pairwise correlations between methylation levels and gene expression）：通过JZP036五组织TPM与启动子甲基化指标相关分析，发现全基因分析中启动子甲基化与表达相关较弱，487个DMR相关基因中相关高于全基因分析；茎、叶、花、根表达正相关，茎-叶和茎-花r分别为0.72和0.59；CHG-CHH甲基化正相关在全基因和DMR相关基因中r分别为0.73和0.84，其他甲基化两两相关弱；外显子重叠DMR基因的甲基化与表达无清晰两两关联；果实表达与其他组织相关在全基因分析不显著，在DMR相关基因中显著正相关，但无正式检验系数差异。结论为启动子甲基化与表达的全基因组关联弱，DMR相关基因中某些关联需谨慎解释。

2.4 JZP036叶片区域甲基化与基因表达关联（Regional methylation associations with gene expression in JZP036 leaves）：通过叶组织匹配的区域分层分析，发现关联具上下文和区域依赖：启动子Spearman系数CpG为-0.075（q<0.001）、CHG为0.015（q=0.0445）、CHH为0.106（q<0.001）；基因体CpG为0.051（q<0.001）、CHG为-0.103（q<0.001）、CHH为-0.092（q<0.001）；重复相关基因体CpG为0.010（q=0.0578）、CHG为-0.151（q<0.001）、CHH为-0.139（q<0.001）；候选6 mA在启动子、基因体和重复相关基因体分别为-0.006（q=0.1233）、0.014（q<0.001）、0.025（q<0.001）。在相同3,667个基因子集中，基因体CHG和CHH系数分别为-0.234和-0.224，重复相关分别为-0.185和-0.212（均q<0.001），CpG基因体系数由0.051变为-0.046（q=0.00683）。结论为CHG和CHH在基因体及重复相关部分呈负相关，候选6 mA信号为描述性且身份和检测准确性未验证。

讨论总结：多数六倍体 A. deliciosa 基因具可检测甲基化，仅少数具全部四类；六等位基因总体甲基化比例接近全基因组平均，但未甲基化位点比例高于单拷贝基因，兼容重复基因间异质表观状态，但不能证明低甲基化导致转录活性、功能分化或剂量缓冲。全基因组甲基化与表达无强相关，提示甲基化单独不能解释表达模式，耦合更应按位点、组织或通路解释并可能与其他调控层包括组蛋白修饰协调。四材料甲基化分布不同，但倍性与物种身份、种质、来源和栽培历史共变，不能归因于倍性。果实表达探索比较中，多数基因对在六倍体较低、少数较高，高TPM比基因涉及植物激素代谢、维生素相关生物合成和防御通路，提示候选功能差异；因数据集物种、遗传背景和处理不同，需匹配比较和表型测量。DMR相关基因涉及光合电子传递、萜类/生物碱生物合成和次生代谢防御，并与MADS-box、WRKY等TF重叠，提供候选基因。相对DMR计数和方向在三次比较中不同，天然二倍体DMR较少，多倍体比较低甲基化DMR更频繁尤其候选6 mA，提供候选区域。优先低甲基化在多倍体系统有报道，但对其材料仍是假设，需匹配二倍体-多倍体谱系、生物重复和独立功能实验检验。局限包括每材料1株、倍性与物种/背景/来源共变、共享六倍体参比可引入材料依赖检测差异、CHG/CHH覆盖染色体较少可由有效基因数和共享基因敏感性解释、叶TPM组织匹配但需重复和匹配外显子注释、跨研究果实TPM混合生物和处理差异、候选6 mA身份和准确性未验证且区域摘要为描述性候选信号。

研究结论翻译：总之，该观察性比较识别了 Actinidia 材料间甲基化模式、DMR和转录本丰度差异。JZP036的组织匹配分析显示甲基化-表达关联具有区域和序列上下文依赖性，基因体及其重复相关部分中CHG和CHH系数为负。这些观察为利用具重复、遗传匹配材料及整合调控测量研究多倍体物种甲基化提供了候选位点和通路优先级。