《Epigenomes》：Low Depth Epigenetic Mapping of Maturation Versus Retrodifferentiation in HepaRG Cells
编辑推荐：
摘要：背景：长读长、单CpG分辨率测序正在重新定义表观基因组学中信息量与测序深度之间的关系。尽管传统甲基组分析通常需要较高覆盖度，研究人员提出了一种可扩展的分析流程，旨在从浅层测序数据中提取高密度的调控逻辑。方法：研究人员利用具有祖细胞样特征的HepaRG细胞
摘要：背景：长读长、单CpG分辨率测序正在重新定义表观基因组学中信息量与测序深度之间的关系。尽管传统甲基组分析通常需要较高覆盖度，研究人员提出了一种可扩展的分析流程，旨在从浅层测序数据中提取高密度的调控逻辑。方法：研究人员利用具有祖细胞样特征的HepaRG细胞系作为肝脏可塑性模型，在两条分化方向不同的发育轨迹中验证该框架：肝细胞成熟以及球体诱导的逆分化。技术策略结合了以CpG为中心的富集分析与区域甲基化聚合分析，从稀疏数据中重建调控景观。研究人员采用长读长Nanopore测序，对5-甲基胞嘧啶（5mC）和5-羟甲基胞嘧啶（5hmC）动态进行了绘制。结果：该流程揭示，这两条轨迹并非彼此相反的过程，而是采用了不同的表观遗传策略。肝细胞成熟的特征是5hmC积累，其部分靶向抑制性异染色质标记H3K9me3、H4K20me3以及先导因子FOXA2。相较之下，逆分化增加了5mC，可能通过与Polycomb相关的网络沉默HNF1A等成人期调控因子。此外，基于聚合的分析能够区分广泛存在的局灶性扰动与一小部分将表观遗传变化转化为区域可及性的转录因子。结论：本研究提供了一个用于研究细胞命运转变的可扩展计算框架，证明即使在降低测序深度的条件下，仍可获得高价值的表观遗传学见解。
该文发表于《Epigenomes》，围绕HepaRG细胞在肝细胞成熟与球体诱导逆分化两种命运转变中的表观遗传调控逻辑展开。研究背景在于，肝脏具有显著的细胞可塑性，成熟肝细胞既可参与组织再生，也可在慢性应激下逆分化为祖细胞样状态，这种双向转换与肝细胞癌（HCC，肝细胞癌）中的干性维持、化疗耐受和复发密切相关。尽管既往研究已提示DNA修饰和代谢重编程参与这些过程，但支撑两种转变的具体表观遗传程序尚未完全阐明。尤其是在长读长Nanopore测序能够直接检测天然DNA上5-甲基胞嘧啶（5mC）与5-羟甲基胞嘧啶（5hmC）的背景下，浅覆盖度数据造成的CpG位点稀疏性，常使依赖连续高密度覆盖的差异甲基化区域分析难以稳定识别关键调控事件。因此，开展本研究的意义在于建立一套适用于低深度测序的可扩展分析框架，从有限数据中提取细胞命运转换的核心调控信息。
研究人员以具祖细胞样特征的人HepaRG细胞系为模型，比较二维分化培养所得成熟样肝细胞状态（Diff）与三维球体培养诱导的逆分化状态（Spheres）之间的5mC和5hmC变化。研究结论表明，这两条轨迹并不是简单互逆过程，而是受本质上不对称的表观遗传程序控制：肝细胞成熟更突出表现为5hmC积累，提示主动去甲基化和精细调控增强；逆分化则以5mC增加为主，倾向于通过更稳定的抑制性网络巩固祖细胞样状态。该研究的重要意义在于证明，即使在长读长浅层测序条件下，通过适当的平滑、功能富集与区域聚合分析，仍可有效解析细胞命运转变中的高价值表观遗传规律，并为研究肝脏可塑性及肿瘤耐药相关机制提供可推广的方法学模板。
在技术方法方面，研究人员使用独立三重复的HepaRG-Diff与HepaRG-Spheres样本，采用Oxford Nanopore长读长平台进行天然DNA修饰测序，以Guppy完成碱基识别，并用Remora同步判定5mC与5hmC。随后，利用DSS基于β-二项分布模型进行差异甲基化位点（DML）分析；以KnowYourCG（KYCG）对差异CpG在转录因子结合位点（TFBS）、染色质状态、组蛋白修饰和CpG岛注释中的富集进行评估；再通过Methylation-Based Inference of Regulatory Activity（MIRA，基于甲基化的调控活性推断）对区域甲基化足迹进行聚合，推断转录因子调控活性。
在结果部分，论文首先以“2.1. Global Imbalance in 5mC and 5hmC in Diff vs. Sphere HepaRG Culture”为题，分析了Diff与Spheres之间5mC和5hmC的整体失衡。研究人员在满足最低覆盖要求的约171,000个CpG位点上开展单碱基差异分析，发现两种状态之间存在显著的5mC重塑，共鉴定约1936个5mC差异甲基化位点，而5hmC差异位点仅159个，数量明显更少。进一步比较显示，Diff相较Spheres更常出现5mC升高，同时5hmC变化更加不对称，绝大多数5hmC差异位点也集中表现为Diff中升高。少数同时发生5mC与5hmC显著变化的位点呈反向相关，说明这两类胞嘧啶修饰在局部位点可能存在动态转换关系。由此得出的结论是：成熟与逆分化之间既有广泛5mC重塑，也存在偏向成熟状态的5hmC积累。
在“2.2. 5mC in Polycomb Subunits Is Differentially Enriched in Diff vs. Sphere HepaRG Culture”部分，研究人员通过KYCG对5mC差异位点的调控注释富集进行分析，发现Diff和Spheres均有大量富集特征，但靶向对象并不相同。Diff中的高甲基化5mC显著富集于CpG岛（CGI）及其邻近区域，并呈现明显的Polycomb Repressive Complex 2（PRC2，Polycomb抑制复合体2）特征，包括H3K27me3、SUZ12和EZH2富集，同时涉及双价转录起始位点。这提示成熟过程中的5mC重塑与PRC2介导的抑制性染色质环境密切相关。相比之下，Spheres中的高甲基化5mC虽也富集于CpG岛相关区域，但更突出地关联H3K9me3、SETDB1、异染色质状态以及PRC1相关成分CBX7、BMI1、PCGF2和RING1，并涉及双价增强子及侧翼区域。该部分结果支持这样一个结论：在两条命运轨迹中，Polycomb复合体的组装与作用模式不同，成熟偏向PRC2/H3K27me3相关程序，而逆分化偏向PRC1/H3K9me3相关的更强抑制性甲基化程序。
在“2.3. Hepatocyte TFs Are Dynamically Modulated by 5hmC in Differentiation Conditions”部分，研究人员重点考察5hmC富集所指向的调控网络。结果显示，Diff中的5hmC富集特征极其丰富，而Spheres中仅有极少数富集项，说明5hmC的动态调控主要发生在成熟条件下。Diff中5hmC不仅指向H3K9me3和H4K20me3等深度抑制性异染色质标记，还富集于ATRX、PARP1及多种锌指蛋白相关区域。更重要的是，在与肝细胞分化相关的转录因子层面，Diff中的5hmC显著富集于POU5F1（OCT4）、NANOG、SOX2等干性主调控因子靶位点，也富集于SOX10、TP63等分化因子靶位点。进一步筛选肝谱系决定相关主转录因子后发现，HNF1A、HNF4A、HNF4G、FOXA1、FOXA2及PROX1等核心肝脏调控因子靶位点在Diff中均偏向获得5hmC。相反，Spheres则对其中部分相同位点，如HNF1A和PROX1，转而出现5mC升高，提示逆分化可能借助更持久的DNA甲基化沉默成人肝调控网络。该部分说明，5hmC在肝细胞成熟中并非被动标记，而是与肝细胞特异性转录网络重塑紧密相关。
在“2.4. TF Activity Inferred from 5mC Aggregation”部分，研究人员利用MIRA基于5mC聚合足迹推断转录因子活性，从138个高可信度人类转录因子中仅识别出6个在两种状态间具有显著活性差异的因子：PAX6、FOXO1、FOXK1、LEF1、PRDM1和ATF4。其中，FOXO1和ATF4在Spheres中的推断活性更高，而PAX6、FOXK1、LEF1和PRDM1在Diff中更高。尽管这些因子并非KYCG分析中的最强富集项，但它们同时也出现在Diff的5hmC富集区域中，提示MIRA检测到的5mC足迹变化可能与Diff环境下活跃的DNA去甲基化过程有关。研究人员据此认为，广泛的局部CpG扰动最终只转化为少数关键转录因子的区域开放与功能输出，这些因子分别关联细胞命运调控以及代谢和应激稳态。
在讨论部分，研究人员综合指出，浅层长读长测序虽然在全基因组尺度上存在明显稀疏性，但通过KYCG和MIRA这类功能聚合策略，仍可从单分子修饰数据中提取具有解释力的全局调控模式。研究表明，HepaRG的成熟与逆分化不是简单镜像，而是依赖不对称的表观遗传机制。成熟过程可能与TET介导的5hmC沉积有关，且这一过程定向发生于既有异染色质区域，提示局部位点存在5mC与5hmC并存的动态周转状态；与此同时，整体5mC可能用于压制替代谱系程序，而局部5hmC则有助于解除特定肝脏增强子的限制。逆分化方面，球体环境通过定向新生5mC沉积沉默HNF1A等成熟肝主调控因子，并在PRC相关位点形成更稳定的抑制状态，从而维持祖细胞样特征。FOXO1和ATF4在Spheres中的活性升高，则提示逆分化细胞还伴随代谢重编程和应激适应。研究人员同时指出，本研究缺乏正交实验验证，所识别的具体DML、转录因子和标志物仍需借助重亚硫酸盐特异性PCR、RT-qPCR或染色质免疫沉淀等方法进一步确认。
研究结论部分可译为：通过结合近期发展的生物信息学算法，可以从浅层长读长测序中提取具有功能意义的5mC/5hmC信息。平滑处理（DSS）、聚合分析（KYCG）和分箱策略（MIRA）为体外两条对比鲜明的分化路径提供了全局视角。研究人员得以检测到5mC/5hmC周转在分化精细调控中的更强影响，Polycomb机器在Spheres中可能被导向巩固性与终末性沉默，以及主导性转录因子活性如何指向逆分化相关化疗耐受的来源。更广泛地看，肝脏可塑性似乎并非由静态的“开/关”状态支配，而是由甲基化—去甲基化偶联循环构成的高度动态区域所调控。
