筛选了100篇文献，终于找到这篇，文章所有复现代码都提供了，单细胞、蛋白质组，学这一篇就够了

科研不掉发快来这个地表最强的生信神仙网站中国银河生信云平台 立即访问https://usegalaxy.cn成年后的神经胶质细胞是否真的已经彻底摆脱了发育期留下的分子印记如果这些印记仍然留在染色质层面它们是帮助细胞在损伤后快速响应还是也可能成为肿瘤发生时被重新利用的“旧程序”2026 年 3 月 19 日Nature Neuroscience发表了 Karolinska Institutet 等团队的研究。作者通过单核染色质开放性、组蛋白修饰和三维染色质构象联合分析发现成人人类中枢神经系统中的少突胶质细胞谱系仍保留部分发育阶段的表观遗传记忆。今天我们就来拆解一下这篇生信文章Single-nucleus epigenomic profiling of the adult human central nervous system unveils epigenetic memory of developmental programs。研究概述这项研究关注的是一个基础但长期缺少直接证据的问题成人人类中枢神经系统不同细胞类型的转录异质性背后究竟由怎样的染色质调控层次支撑。作者以成人人类初级运动皮层、小脑和颈段脊髓为对象系统绘制了单核水平的染色质开放性和 H3K27ac、H3K27me3 图谱。研究的核心结论有两点一是少突胶质细胞谱系存在新的 SOX10 远端增强子二是脊髓来源的少突胶质细胞和星形胶质细胞在 HOX 位点保留了类似发育期的 primed 染色质状态而这种状态与成熟细胞中的稳定高表达并不一致提示成人细胞中存在“表观遗传记忆”。实验设计样本部分作者收集了 20 例成人供体的尸后冷冻白质样本年龄 34 至 74 岁男女各半来自初级运动皮层、 小脑和颈段脊髓三个区域。关键实验主要分为四部分。首先作者对 60 份组织样本进行了snATAC-seq用于分析单核染色质开放性。其次在来自 4 名供体的脊髓和皮层样本中作者建立了可同时检测H3K27ac与H3K27me3的nanoCUTTag多模态数据。第三作者在颈段脊髓和皮层样本中开展10x Multiome联合分析 ATAC 与 RNA以比较开放性与转录输出是否一致。第四作者在 iPSC 来源的人 OPC 中进行高分辨率Micro-C并结合CRISPRa/CRISPRi对候选增强子进行功能验证。分析流程上snATAC 数据用于细胞聚类、细胞类型注释、差异开放区分析、co-accessibility 推断和 TF motif 分析nanoCUTTag 数据用于比较活性和抑制性组蛋白修饰分布Multiome 用于检验开放性与 RNA 是否耦合Micro-C 则用于解析 HOX 基因簇周围的三维结构和远端增强子联系。研究结果图 1作者首先构建了成人人类 CNS 的 snATAC-seq 与 H3K27ac/H3K27me3 nanoCUTTag 图谱成功识别主要神经元和胶质细胞类型并显示 H3K27ac 比单纯开放性更能区分细胞类型特异的调控状态。图 2通过 co-accessibility、footprinting 和 CRISPRa 验证作者在少突胶质细胞谱系中鉴定出一个位于SOX10上游约 91 kb 的新远端增强子并证明其可特异性提升SOX10表达。图 3基于 ATAC 与 H3K27ac 数据构建的核心转录因子调控网络显示不同神经细胞谱系具有不同的主导 TF而 HOX 家族 TF 在 OPC 中的调控强度高于成熟少突胶质细胞提示其在祖细胞状态中更活跃。图 4区域比较发现颈段脊髓来源的少突胶质细胞在大量 HOX 基因位点具有更高开放性但 Multiome 结果显示这些基因只有很弱的 RNA 表达说明开放性和转录输出发生了解耦。图 5在脊髓来源少突胶质细胞中HOXA 和 HOXD 基因簇呈现 3 端开放并富集 H3K27ac、5 端富集 H3K27me3 的梯度结构形成 inactive、primed、silenced 三类染色质域显示典型的发育样表观构型。图 6iPSC 来源 hOPC 的 Micro-C 表明 HOX 基因簇的三维染色质组织与发育期 c-Dom/t-Dom 结构相似且针对MIR10B/LINC01116相关增强子的 CRISPRi/a 可调控同一结构域内的 3 HOXD 基因表达。图 7成人脊髓少突胶质细胞中 primed 的 HOX 位点在桥脑 H3K27M 高级别胶质瘤中对应为被重新激活的位点提示正常细胞中的发育记忆可能为特定脑肿瘤提供易感的染色质背景。数据分析生信分析文中涉及的组学技术包括snATAC-seq、单细胞nanoCUTTagH3K27ac 与 H3K27me3、10x MultiomeATAC RNA、Micro-C以及用于外部参照的scRNA-seq、公开 Hi-C/ChIP-seq 数据。1. 单核染色质开放性组snATAC 原始数据使用cellranger-atac count处理并以cellranger-atac aggr聚合。作者用自定义脚本构建 2-kb genome bin count matrix使用ArchR计算 TSS enrichment只保留TSSe 7且唯一片段数大于 3,000 的细胞。下游分析中作者保留前 100,000 个特征使用 TF-IDF 归一化、SVD 降维、kNN 图构建和 Leiden 聚类并通过Harmony做批次校正。差异开放区分析同时使用episcanpy rank_features与diffxpy并将性别、年龄和组织区域纳入协变量多重检验使用 Benjamini-Hochberg 校正。细胞类型注释方面作者构建了基于启动子附近 5 kb 区域的 gene activity matrix并用MAGIC平滑再结合既往 snRNA-seq 数据中每类细胞前 50 个差异基因构建 metagene score从而完成细胞类型标注。2. 组蛋白修饰组作者在归档人脑组织上建立了可同时检测H3K27ac与H3K27me3的nanoCUTTag。下游将 H3K27ac 数据与 snATAC 数据整合后做 label transfer从而将组蛋白修饰信号映射到具体细胞类型。作者还基于 10-kb 基因组分箱将 ATAC、H3K27ac、H3K27me3 三种模态的信号输入deepTools computeMatrix对全基因组 bin 做层次聚类和 Pearson 相关分析用来描述不同细胞类型和不同模态之间的全局对应关系。对 H3K27ac/H3K27me3 的基因层面信号作者用deepTools中k-means将 1-kb padded TSS 的信号分为 10 类再结合gget enrichr进行细胞类型富集解释。3. 转录组与多组学联合分析Multiome 数据用于联合比较染色质开放性和 RNA 表达。作者将 Multiome ATAC 与 snATAC 图谱整合并用配对 RNA 信息完成注释从而直接比较 OLG marker 与 HOX 基因在不同区域中的“启动子开放性-转录表达”关系。此外研究还整合了既往公开的人脑scRNA-seq数据集和作者自建的前额叶scRNA-seq数据进行标签转移所用方法包括Seurat、Signac、FindTransferAnchors和TransferData。4. 调控元件与 TF 网络分析候选增强子分析中作者使用Cicero进行 co-accessibility 推断分别以 0.25 和 0.5 作为显著交互和高可信交互阈值。TF motif 层面作者使用chromVAR分析 motif accessibility 差异并使用Regulatory Genomics Toolbox中的HINT做 footprinting识别活跃 TF 结合位点。核心 TF 调控网络则通过Core Regulatory Circuit框架构建输入包括 H3K27ac 或 ATAC 的 BAM/bigwig 和 super-enhancer 信息。TF 强度定义为出边数减入边数用于评估某一 TF 在网络中的调控影响力。5. 三维基因组与多组学联合分析Micro-C 数据使用 Dovetail pipeline 预处理之后用cooler做矩阵平衡与归一化并生成 observed/expected contact matrix。hicPCA用于 compartment 分析coolpuppy用于 aggregate pileupcooltools用于 insulation score 与边界强度分析mustache用于 loop callinghicPlotViewpoint用于 virtual 4C。作者将 Micro-C 结果与成人少突胶质细胞的 ATAC、H3K27ac、H3K27me3 图谱叠加比较 HOX 基因簇周围的 TAD、子结构域边界和远端增强子联系并进一步与桥脑、丘脑 H3K27M 高级别胶质瘤以及 PFA ependymoma 的公开三维结构数据比较。统计分析统计信息主要由图注和方法部分提供。分组比较中差异开放区分析使用双侧t test并进行Benjamini-Hochberg多重校正图 2 的 qPCR 比较使用双侧t test图 6 的 CRISPRi/a 实验同时使用单因素ANOVA加Tukey事后检验以及双侧t test。文中还使用线性回归展示 OPC 与 MOL 中共享 TF 强度之间的关系。对于部分 qPCR 实验作者说明默认假设数据服从正态分布但未进行正式检验。总体来说统计方法与多重校正策略在关键分析处有交代但并非所有图形的统计细节都展开到非常充分的程度。总结研究意义这项研究补上了成人人类 CNS 单核表观组学中的一个关键空白即在单细胞分辨率下同时描绘染色质开放性和活性/抑制性组蛋白修饰并将这些信息与三维染色质结构联系起来。从生物学上看文章说明成人少突胶质细胞和部分星形胶质细胞并没有完全抹去发育时期的 HOX 调控框架而是以一种低表达、但可被再次调用的 primed 状态将其保留下来。这为理解胶质细胞在损伤修复、再髓鞘化以及肿瘤易感性中的区域差异提供了新的线索。从转化角度看研究提示正常组织中的表观遗传记忆既可能支持再生也可能在 PRC2 失衡等病理背景下被肿瘤重新利用因此 HOX 调控域及其远端增强子可能成为后续研究脑胶质瘤细胞起源和区域特异性的切入点。文章复现原始数据已提供。•EGAD50000000410scATAC-seq、nanoCT-seq、sc multiome 和 hOPCs Micro-C 数据•EGAD50000001542人前额叶活检 scRNA-seq 数据•EGAD50000001535B cells Micro-C 数据代码已提供。•https://github.com/mkabbe/snATACnanoCT_AdultHumanCNS中国银河生信云平台精品课程中国银河生信云平台UseGalaxy.cn致力于生信平权。海量云端算力、8000生信工具结合AI推动生信进入3.0时代数据分析从本地到云端从手工到 AI。加入交流群免费领取学习资料。左手代码右手云平台。特色生信培训助你丝滑发顶刊单细胞数据分析培训班Python/Galaxy可选不怕学不会转录组数据分析实战Galaxy| 直播回放咨询小助手usegalaxy