人类肠道染色质可及性的单细胞分析识别调控程序并阐明克罗恩病的遗传关联

Single-cell analysis of chromatin accessibility in the human intestine identifies regulatory programs and clarifies genetic associations in Crohn's disease

作者信息Yu Zhao, Ran Zhou, Zepeng Mu, Peter Carbonetto, Xiaoyuan Zhong, Bingqing Xie, Kaixuan Luo, Zhiwei Jiang, Jianqiao Liu, Candace M Cham, Jason Koval, Xin He, Andrew W Dahl, Xuanyao Liu, Eugene B Chang, Anindita Basu, Sebastian Pott
PMID42791370
期刊Nat Genet
发布时间2026-10
DOI10.1038/s41588-026-02755-z

实验完整度

高

包含scATAC–seq、scRNA-seq、caQTL、sLDSC和精细定位等多个独立证据层级,且有功能验证和机制验证

主要模型

23例克罗恩病患者及16例健康对照的末端回肠和升结肠活检组织 178,030个细胞进行scATAC–seq 29种细胞类型聚类

重点核对

scATAC–seq数据保留178,030个细胞,鉴定29个细胞类型聚类 topic模型识别45个topic,其中topic 7和topic 9与炎症相关 caQTL分析在19种细胞类型中检测到21-8,533个caQTL sLDSC分析显示10种免疫细胞类型与IBD风险富集相关 精细定位识别30个PIP≥0.2的变异位于cCREs内

摘要

克罗恩病是一种复杂的炎症性肠病,由遗传、微生物和环境因素相互作用引起。细胞类型特异性的克罗恩病病因和遗传风险贡献尚未完全理解。在这里,我们建立了一个综合性的细胞类型解析染色质可及性图谱,包括来自23名活动性和非活动性克罗恩病患者及16名健康对照的末端回肠和升结肠中的557,310个候选顺式调控元件(cCREs)。我们识别了细胞类型特异性、解剖位置特异性和情境特异性的cCREs,并表征了克罗恩病患者肠道黏膜中炎症反应背后的调控程序。这些细胞类型解析数据证实,克罗恩病遗传力主要在适应性免疫细胞,特别是T细胞,以及先天免疫细胞,包括中性粒细胞中富集。使用精细定位的非编码克罗恩病变异,我们识别了30个位于cCREs内的变异。我们的图谱为研究克罗恩病和健康中的基因调控效应提供了全面资源,并突出了克罗恩病风险背后的细胞复杂性。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
克罗恩病中细胞类型特异性的染色质可及性调控程序是什么,以及哪些细胞类型和变异介导遗传风险?
核心机制
炎症相关调控程序涉及AP-1因子(如FOSL1/FOSL2)和干扰素调控因子(如IRF2、IRF3、IRF9)等转录因子,以及JAK-STAT信号通路;caQTL分析显示TF结合基序破坏影响染色质可及性,如TCF3在浆细胞和CDX2在结肠细胞中。
主要证据
通过scATAC–seq构建557,310个cCREs图谱,topic模型识别炎症相关topic 7和9,caQTL分析显示基序破坏与可及性相关,sLDSC显示免疫细胞富集IBD风险,精细定位将30个变异映射到cCREs。
研究意义
提供了研究克罗恩病和健康中基因调控效应的综合资源,并强调了克罗恩病风险背后的细胞复杂性,有助于生成关于非编码变异在基因调控和疾病中作用的机制假设。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

样本收集与scATAC–seq图谱构建

获取肠道黏膜细胞类型解析的染色质可及性图谱

从23例克罗恩病患者和16例健康对照的末端回肠和升结肠收集活检,进行scATAC–seq,质控后保留178,030个细胞,聚类为29种细胞类型,鉴定557,310个cCREs。

2

主题模型分析识别调控程序

识别细胞类型特异性和炎症相关调控程序

使用FastTopics对scATAC–seq数据进行主题建模,选择K=45,识别细胞类型特异性和炎症相关topic(如topic 7和9),进行基序富集和通路富集分析。

3

炎症相关差异可及性区域(iDARs)鉴定

发现炎症条件下细胞类型特异性的染色质可及性变化

使用dreamlet进行伪批量差异可及性测试,比较炎症样本与对照样本,识别iDARs,进行基序富集分析。

4

caQTL定位和基序破坏分析

探究遗传变异对染色质可及性的影响及TF基序破坏效应

使用RASQUAL在19种细胞类型中调用caQTLs,使用motifbreakR评估基序破坏,计算caQTL效应大小与基序破坏分数的相关性。

5

遗传力富集分析

评估细胞类型和topic特异性染色质可及性对IBD风险的富集

使用sLDSC评估IBD及其亚型在细胞类型和topic特异性开放染色质中的遗传力富集,包括比较其他自身免疫和对照性状。

6

精细定位和变异优先级排序

优先考虑非编码克罗恩病风险变异并关联到细胞类型

使用SuSiE对克罗恩病GWAS信号进行精细定位,获得PIPs,优先考虑PIP≥0.2且位于cCREs内的变异,使用topic加载注释细胞类型和上下文。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
Cell Ranger ATAC10x Genomics--
ArchR----
MACS2----
GLIMPSE2----
RASQUAL----
SuSiE----

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
样本收集
队列规模:23例克罗恩病患者和16例健康对照;组织类型:末端回肠和升结肠活检
阅读提示:Methods中的Ethics statement和Patient information
scATAC–seq
细胞保留数量:178,030个细胞;细胞类型数:29种
阅读提示:Results中的A regulatory atlas of the intestinal mucosa和图2
主题模型
主题数量:K=45
阅读提示:Methods中的Topic modeling and pathway enrichment
caQTL定位
细胞类型数:19种;FDR<10%和MAF>5%
阅读提示:Methods中的CaQTL calling
精细定位
PIP≥0.2;使用UKBB参考面板
阅读提示:Methods中的Statistical fine-mapping