《International Journal of Molecular Sciences》：Deciphering Cancer-Specific H3K27ac Landscapes Through Systematic Multi-Omics Integration and Development of Open-Access Genomic Annotation Platform
编辑推荐：
H3K27ac（组蛋白H3第27位赖氨酸乙酰化）通过表观遗传方式调控基因表达，并可作为癌症状态的指示标志。尽管异常H3K27ac已被广泛研究，但癌症特异性H3K27ac模式的全面景观仍然缺乏。研究人员通过整合覆盖48种组织和细胞类型的209套H3K27ac染色
H3K27ac（组蛋白H3第27位赖氨酸乙酰化）通过表观遗传方式调控基因表达，并可作为癌症状态的指示标志。尽管异常H3K27ac已被广泛研究，但癌症特异性H3K27ac模式的全面景观仍然缺乏。研究人员通过整合覆盖48种组织和细胞类型的209套H3K27ac染色质免疫沉淀测序（ChIP-seq）数据集，定义了3612个普遍保守区域和大量癌症特异性H3K27ac改变。引人注目的是，尽管癌症特异性H3K27ac获得（gains）在白血病、肝细胞癌（HCC）、乳腺癌（BRCA）、肺腺癌（LUAD）和淋巴瘤中普遍存在，但丢失（losses）却极为罕见。这些癌症特异性H3K27ac获得主要富集于非编码区域，并与染色质可及性增加和染色质相互作用增强相关，从而使致癌转录因子（TF）被募集以激活靶基因。相反，癌症特异性H3K27ac丢失优先发生在启动子，并与染色质致密化和转录抑制相关。白血病中的机制性多组学分析进一步证明了这种表观遗传-结构偶联。白血病特异性H3K27ac获得建立了允许性染色质环境，促进致癌性GATA1结合及其靶基因TAL1的激活。在此背景下，TAL1与STIL的相互作用增加，STIL是已知参与白血病发生的融合伴侣。为促进这些发现的更广泛应用，研究人员开发了一个开放获取的基因组注释网络服务器（web server）。总之，该研究提供了癌症特异性H3K27ac修饰的全面景观，将其定义为致癌性染色质重塑的功能性表观遗传特征，并为理解癌症基因组调控提供了注释资源。
研究背景与问题：全球癌症负担沉重，2022年全球新发癌症超过2000万例、癌症相关死亡970万例，中国新发482万例、死亡257万例。癌症由遗传变异、表观遗传修饰和环境因素复杂交互所致，表观遗传失调是连接遗传变异与环境暴露的重要机制之一。启动子、增强子和基因体上不同组蛋白修饰（HMs）组合定义疾病相关染色质状态并贡献癌症亚型异质性。H3K27ac（组蛋白H3第27位赖氨酸乙酰化）是活跃染色质的标志，可通过中和组蛋白尾部电荷降低核小体DNA亲和力，增加转录因子（TF）结合所需的染色质可及性。异常H3K27ac与多种癌症癌基因激活相关，例如BRCA中NBR2调控增强子H3K27ac激活TBX1，结直肠癌中H3K27ac介导STEAP3激活并与ATF3形成复合物；Zhou et al.和Cai et al.的全基因组H3K27ac分析识别了超级增强子及相关致癌靶点。然而，既往研究多关注癌变过程中的差异H3K27ac区域，对仅见于特定癌种、其他癌种和正常样本缺失的癌症特异性H3K27ac事件关注有限，且这些修饰如何影响染色质结构和转录仍 largely uncharacterized。为此，研究人员构建多癌H3K27ac图谱，系统识别高度保守区域和各癌种癌症特异性改变，整合多组学数据并在白血病模型验证，开发公开网络服务器CCHMR（癌症特异性或跨组织/细胞类型保守的组蛋白H3第27位赖氨酸乙酰化修饰区域），为癌症基因组调控提供注释资源。
关键技术方法概括：研究人员从ENCODE和GEO收集209套H3K27ac ChIP-seq（染色质免疫沉淀测序）数据，覆盖48种细胞或组织类型，并纳入24套H3K4me1（组蛋白H3第4位赖氨酸单甲基化）数据；整合112套RNA-seq、22套DNase-seq（DNase I超敏位点测序）、27套ATAC-seq（转座酶可及性染色质测序）、8套细胞Hi-C（染色体构象捕获）和10套白血病患者Hi-C、10套白血病患者RNA-seq，以及GATA1、MEF2B和CTCF（CCCTC结合因子）ChIP-seq数据。扰动数据包括CRISPR介导p300（组蛋白乙酰转移酶p300）耗竭后2套H3K27ac ChIP-seq和K562细胞经SGC-CBP30或C646处理的8套ATAC-seq。生存分析收集GEO和TARGET数据库9个基因表达队列（包括GSE71014和TARGET等），并使用BART、GENIE3、LASSO-Cox（最小绝对收缩和选择算子-Cox回归）等分析。
研究结果：
2.1 全基因组识别与表征保守及癌症特异性H3K27ac模式。2.1.1 全球H3K27ac分析定义核心保守H3K27ac修饰区域：对209套数据peak calling获得9,619,866个H3K27ac峰，中位数41,192个/套；依据相邻峰summit距离约150 bp拐点合并为678,446个非重叠union区域；以占用频率≥80%（≥168套数据）定义3612个保守区域，并通过幂律模型显示观测/期望比值达3.26。2.1.2 广泛H3K27ac分析揭示癌症特异性H3K27ac景观：比较各癌种与其他样本及匹配正常对照的占用和信号，识别获得与丢失区域；白血病中1322个获得、32个丢失，获得区域信号较正常升高5.55–30.64倍，丢失区域平均仅0.24倍，效应量分析显示获得区域效应更大。2.1.3 癌症特异性H3K27ac区域具有高特异性和空间偏好：跨癌种Jaccard指数平均为获得0.002、丢失0.01，远低于全体H3K27ac区域的0.12–0.24；获得区域主要富集于远端非编码和内含子区域，丢失区域在全部癌种中>90%位于启动子。
2.2 癌症特异性H3K27ac修饰事件与癌变过程中染色质结构动态变化相关。2.2.1 癌症特异性H3K27ac修饰伴随局部染色质状态改变：DNase-seq和独立ATAC-seq显示获得区域可及性升高、丢失区域降低、保守区域稳定；K562细胞经SGC-CBP30和C646处理后，获得区域染色质可及性显著降低，保守和丢失区域相对稳定。2.2.2 癌症特异性H3K27ac获得和丢失平行于高级染色质重塑：CTCF结合在获得区域附近升高，H3K4me1信号与H3K27ac变化一致；Hi-C显示获得和丢失区域侧翼相互作用频率分别增加和降低，白血病中57%获得区域在1000 kb内维持升高相互作用，62.5%丢失区域在该尺度持续降低。
2.3 癌症特异性H3K27ac改变与转录重编程相关，获得区域对p300耗竭更敏感：H3K27ac获得与基因上调相关、丢失与基因下调相关，保守区域变化少；癌症特异性区域尤其丢失区域与基因表达变化相关性最高；CRISPR介导p300耗竭后，获得区域H3K27ac信号显著降低，丢失区域基本不受影响。
2.4 癌症特异性H3K27ac获得可能促进致癌TF募集和癌症相关转录重编程。2.4.1 癌症特异性H3K27ac获得区域可能促进致癌TF募集：BART推断TF富集，淋巴瘤获得区域最显著富集MEF2B，丢失区域富集KDM5C；BRD4在白血病、HCC和LUAD丢失区域中排名靠前。2.4.2 GATA1–TAL1轴可能驱动白血病发生并预测患者生存：GATA1在白血病特异性获得区域富集最强，其结合信号在获得区域高于保守和丢失区域；GATA1靶基因主要参与红系和髓系发育调控；交集获得区域相关基因与GATA1靶基因得到12个基因，LASSO-Cox保留7个，可将1153例白血病患者分为高低风险组，高风险组表达更高，TARGET队列结果类似；TAL1为GENIE3权重最高的GATA1调控基因，白血病TAL1位点附近两个获得区域与STIL等邻近基因空间相互作用增强，并在GSE152136队列验证。
2.5 CCHMR：功能性H3K27ac注释交互平台：研究人员开发CCHMR网络服务器，用户可快速判断目标基因组区域是否与3612个高度保守H3K27ac区域或癌症特异性H3K27ac区域重叠，支持手动输入和BED格式文件批量上传。
讨论部分总结：H3K27ac不仅参与癌症，也参与阿尔茨海默病和心肌梗死等疾病，但既往研究多聚焦癌变过程中差异H3K27ac信号，未充分考虑疾病特异性H3K27ac模式及其调控效应。该研究构建多癌H3K27ac图谱，识别5种癌症的癌症特异性区域及48种细胞和组织类型的保守区域，整合组蛋白修饰和TF ChIP-seq、RNA-seq、ATAC-seq、DNase-seq及原位Hi-C，发现癌症特异性H3K27ac获得与染色质可及性增加、长程染色质相互作用和致癌转录程序相关，获得区域可作为转录允许性枢纽，例如白血病中GATA1富集于白血病特异性获得区域，其靶基因TAL1在高风险病例中上调，并获得与STIL等邻近基因增强的空间邻近；癌症特异性H3K27ac丢失主要位于启动子，与染色质致密化和基因抑制相关。局限性包括：主要基于多来源公共数据再分析，扰动验证主要在白血病，数据质量、测序深度、样本量和实验设计差异可能影响一致性；仅聚焦高度癌症特异性改变，未系统研究广泛共享或较不动态区域；多数分析为bulk水平，可能掩盖细胞类型特异性异质性；H3K27ac改变、染色质重塑和转录调控之间的关联不必然建立直接因果关系，未来需单细胞多组学、靶向功能扰动和体内验证。
研究结论翻译：总之，该研究建立了高分辨率、多癌种H3K27ac修饰区域图谱，成功将定义恶性身份的调控区域与高度保守区域解耦。整合多组学分析提示，癌症特异性H3K27ac获得与动态调控枢纽相关，这些枢纽以致癌TF募集和三维染色质重塑为特征。通过以GATA1–TAL1轴作为这种表观遗传-结构偶联的示例，研究结果为理解局部组蛋白修饰与染色质组织更广泛变化之间的潜在关系提供了框架。此外，CCHMR网络服务器为更广泛研究群体提供了可转化资源，以在不同生物学背景中优先评估功能性调控元件。最终，研究结果描绘了癌症表观基因型的全面景观，为发现下一代表观遗传生物标志物和治疗靶点提供了稳健框架。
