子宫肌瘤与配对肌层中非编码RNA启动子甲基化状态的比较分析

《International Journal of Molecular Sciences》：A Comparative Analysis of the Methylation Status of Non-Coding RNA Promoters in Fibroid and Matched Myometrium

子宫肌瘤表现出非编码RNA（ncRNA）的表达失调，但其潜在机制仍不完全清楚。研究人员探究了启动子DNA甲基化及其与子宫肌瘤中ncRNA表达的关系。采用MeDIP–chip分析8对配对子宫肌瘤和匹配肌层组织来源的基因组DNA，以鉴定差异甲基化的ncRNA启动子。选取候选位点在16对样本中用甲基化特异性PCR（MSP）验证，并在68–94对标本中用qRT–PCR评估转录本表达。MeDIP–chip鉴定出538个lncRNA和61个miRNA具有差异启动子甲基化，包括300个高甲基化和238个低甲基化lncRNA，以及47个高甲基化和14个低甲基化miRNA。启动子甲基化与转录本表达无显著相关（r = ?0.1224）。MSP证实LINC–PINT和MIR9–3高甲基化，WT1–AS和TTLL10–AS1低甲基化。相应地，子宫肌瘤中LINC–PINT和MIR9–3表达降低，而WT1–AS和TTLL10–AS1表达升高。然而，LINC–PINT和TTLL10–AS1的甲基化与MeDIP–chip结果未完全对应。这些结果揭示了子宫肌瘤中广泛的ncRNA启动子甲基化改变，但表明全基因组甲基化并不能一致预测转录本表达，突出了ncRNA表观遗传调控的复杂性及位点特异性验证的重要性。

研究背景与立项依据

子宫肌瘤影响近70%育龄女性，其分子驱动机制尚未完全阐明。除蛋白编码基因外，非编码RNA（ncRNA，包括长链非编码RNA lncRNA和微RNA miRNA）被证实受表观遗传调控并参与肌瘤发病。启动子DNA甲基化与转录活性呈反比虽是经典原则，但全基因组甲基化能否一致预测ncRNA表达在肌瘤中未知。既往研究多聚焦蛋白编码基因甲基化或单个ncRNA，缺乏ncRNA启动子甲基化的系统图谱及其与表达整合分析。该研究由此开展，旨在绘制子宫肌瘤ncRNA启动子甲基化谱并评估其与表达的关系，识别候选表观遗传调控ncRNA，为生物标志物和非激素治疗提供依据。论文发表于《International Journal of Molecular Sciences》。

主要技术方法与样本队列

研究收集美国Harbor UCLA Medical Center妇科手术获得的68–94对肌壁间子宫肌瘤（2–5 cm）及配对肌层组织，供者30–49岁、术前至少3个月未用激素。发现阶段取8对样本做MeDIP–chip（Arraystar人ncRNA启动子芯片，覆盖27248个lncRNA启动子及622个miRNA启动子）；同批6对样本有既往RNA–seq数据用于甲基化–表达整合。验证阶段以16对样本做亚硫酸氢盐转化后甲基化特异性PCR（MSP），以68–94对样本做qRT–PCR验证LINC–PINT、WT1–AS、TTLL10–AS1、MIR9–3表达；大队列中行MED12 exon 2 Sanger测序分型。统计采用非参数检验（Wilcoxon配对、Mann–Whitney）。

研究结果

2.1 MeDIP–Chip分析揭示肌瘤与配对肌层间lncRNA和miRNA启动子差异甲基化

8对样本MeDIP–chip经标准化与分层聚类（?log₁₀ p ≥ 2）鉴定538个lncRNA启动子条目与61个miRNA启动子差异甲基化。其中lncRNA高甲基300（55.8%）、低甲基238（44.2%）；miRNA高甲基47（77.0%）、低甲基14（23.0%）。按CpG密度分：高CpG密度启动子（HCP）、中CpG密度启动子（ICP）、低CpG密度启动子（LCP）在高低甲基组中各占约35%、11%、19%。结论为肌瘤存在广泛ncRNA启动子甲基化改变。

2.2 肌瘤中差异启动子甲基化lncRNA和miRNA通路分析揭示受影响通路

对300高甲基/238低甲基lncRNA及47高甲基/14低甲基miRNA做GO、KEGG、Reactome富集。低甲基lncRNA富集于“癌症中miRNA”“癌症蛋白聚糖”“细胞衰老”；高甲基lncRNA富集于“HPV感染”“癌症中miRNA”“蛋白聚糖”。低甲基miRNA连向鞘脂信号、PI3K/AKT、组蛋白乙酰转移酶结合；高甲基miRNA连向癌症中miRNA、FoxO、核外雌激素信号。表明差异甲基化ncRNA关联癌症相关、激素应答及PI3K/AKT网络。

2.3 qPCR和甲基化特异性PCR验证芯片数据

整合MeDIP–chip与同6对RNA–seq（基因水平），Spearman r = ?0.1224（p = 0.0883），无显著相关。在69个满足RNA–seq差异表达（倍数≥1.5或≤0.67）的位点中，40个（58.0%）呈典型反比（低甲基高表达29、高甲基低表达11），29个（42.0%）非反比。选4个高GC/CpG候选做MSP（n=16对）与qRT–PCR（n=68–94对）：肌瘤中WT1–AS、TTLL10–AS1表达升，LINC–PINT、MIR9–3表达降；按MED12突变分层，LINC–PINT在突变瘤更低、WT1–AS更高。MSP确证LINC–PINT与MIR9–3高甲基、WT1–AS与TTLL10–AS1低甲基；但LINC–PINT与TTLL10–AS1的MeDIP–chip与MSP方向不完全一致。说明启动子甲基化可参与部分ncRNA调控但不统一预测表达。

讨论与结论翻译

讨论指出，该研究通过MeDIP–chip结合靶向验证给出肌瘤ncRNA启动子甲基化全景，证实广泛改变且MSP验证WT1–AS/TTLL10–AS1低甲基、LINC–PINT/MIR9–3高甲基，但LINC–PINT与TTLL10–AS1的芯片与MSP不一致源于MeDIP–chip检区域富集、MSP检特异CpG位点及异质性。全队列甲基化–表达无显著相关，69个位点中仅58.0%呈经典反比，提示ncRNA表达受染色质重塑、组蛋白修饰、转录因子、增强子及转录后机制共同塑造。所验四位ncRNA（LINC–PINT抑癌、WT1–AS涉WT1/p53网络、TTLL10–AS1乏表征、MIR9–3沉默调NF?κB）可能经Wnt/β?catenin、NF?κB等促增殖促ECM积累。局限含发现队列小、MeDIP分辨率低于亚硫酸氢盐测序、MED12分层仅验证队列可做、芯片用预设阈值非FDR校正、无功能因果实验。

结论（翻译）：研究发现全基因组启动子甲基化谱不能一致预测子宫肌瘤中ncRNA表达，突出超越DNA甲基化本身的表观遗传调控复杂性；尽管启动子甲基化可参与选定ncRNA调控，但染色质重塑、组蛋白修饰、转录因子活性及转录后调控共同塑造肌瘤ncRNA格局。结果强调解读全基因组甲基化数据时需位点特异性验证，并为未来整合表观遗传网络及评估ncRNA作为肌瘤生物标志物或治疗靶点奠基。