使用scHiCAR对复杂组织中的转录组、表观基因组和三维基因组进行三模态单细胞分析

Trimodal single-cell profiling of transcriptome, epigenome and 3D genome in complex tissues with scHiCAR.

作者信息Xiaolin Wei, Yueyuan Xu, Dongchan Yang, Kyukwang Kim, Lizhi Yi, Wenzhe Luo, Xin Lin, Yu Xiang, Ashley B Williams, Xiaotao Wang, Sweta Srivas, Christabel Tan, Kan Zhang, Wei Li, Yang Eric Li, Feng Yue, Z Josh Huang, Inkyung Jung, Yarui Diao
PMID41714417
发布时间2026-02-19
DOI10.1038/s41587-026-03013-7

实验完整度

高

包含体外细胞系、小鼠脑组织、骨骼肌再生模型等多层级验证,并有功能与机制验证

主要模型

H1人胚胎干细胞 GM12878人淋巴母细胞系 小鼠胚胎干细胞 小鼠额叶皮层组织 小鼠骨骼肌再生模型

重点核对

细胞类型注释准确性:基于BC1和标记基因,三种细胞系分类准确率100% scHiCAR与bulk HiCAR开放染色质峰重叠比例:H1中59%,GM12878中81% 长程顺式相互作用比例:scHiCAR为48%,Dip-C为~6%,LiMAC为~10% scHiCAR环与bulk HiCAR环重叠:H1中49%(5kb),GM12878中69%(5kb) scDeepLUCIA在低深度数据中的AUROC:0.95,优于CovNorm(0.74)和Peakachu(0.65)

摘要

顺式调控元件(CRE)的三维(3D)组织在转录控制中至关重要。然而,从同一个单细胞中捕获转录组、表观基因组和3D基因组仍然具有挑战性。在此,我们提出了scHiCAR(单细胞Hi-C结合ATAC和RNA-seq),一种基于平板组合条形码的方法,可同时从同一细胞中分析mRNA、开放染色质和染色体构象捕获。与现有的单细胞3D基因组方法相比,scHiCAR更有效地富集锚定在候选CRE(cCRE)上的长程顺式相互作用。应用于162万个小鼠脑细胞并结合基于深度学习的环调用器,scHiCAR准确地在22种脑细胞类型中定义了细胞类型特异性转录组、可及的cCRE和5kb分辨率的增强子-启动子对。scHiCAR在骨骼肌等具有挑战性的组织中也表现稳健,能够对肌肉干细胞再生过程中的基因调控动态进行三模态单细胞水平分析。通过提供一种可扩展且成本效益高的系统,用于复杂组织中基因调控景观的单细胞三模态分析,scHiCAR揭示了3D基因组重组在转录控制中的基因座特异性调控作用。

实验结论

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

研究问题
能否从同一个单细胞中同时捕获转录组、表观基因组和3D基因组,并用于解析复杂组织中的基因调控?
核心机制
scHiCAR通过组合条形码和Tn5转座酶介导的开放染色质标记,结合原位连接和环化,特异性富集锚定在开放染色质区域的长程顺式染色质相互作用。
主要证据
在混合细胞系中验证细胞类型分类准确率100%;在小鼠额叶皮层中分析162万个细胞,定义22种细胞类型和294,395个5kb染色质环;在骨骼肌再生中揭示基因座特异性的3D基因组重组动态。
研究意义
scHiCAR提供了一种可扩展、成本效益高的三模态单细胞分析系统,能够揭示复杂组织中3D基因组重组在转录控制中的基因座特异性调控作用。

研究路径

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

1

scHiCAR技术开发与验证

验证scHiCAR同时捕获转录组、开放染色质和3D基因组的能力及准确性

在H1 hESC、GM12878和mESC混合细胞中应用scHiCAR,通过BC1条形码作为 ground truth 评估细胞类型分类准确性,并与Dip-C、LiMAC等已发表数据比较开放染色质富集和染色质接触捕获效率。

2

小鼠额叶皮层复杂组织中的应用

评估scHiCAR在复杂组织中的性能,并定义细胞类型特异性转录组、开放染色质和3D基因组特征

对小鼠额叶皮层细胞进行scHiCAR分析,基于RNA、ATAC和染色质接触数据聚类,鉴定22种主要脑细胞类型,并与BICCN snATAC-seq和bulk HiCAR数据比较。

3

大规模小鼠脑细胞图谱构建

展示scHiCAR的高通量能力,并在大规模细胞中定义细胞类型特异性调控元件和3D基因组特征

对162万个小鼠脑额叶皮层细胞核进行scHiCAR分析,使用组合条形码,鉴定22种主要细胞类型,定义554,946个开放染色质峰和260,781个独特cCRE,并与BICCN数据比较。

4

scDeepLUCIA环调用器开发与验证

开发一种对细胞丰度稳健的深度学习环调用器,用于从scHiCAR数据中识别高分辨率染色质相互作用

基于DeepLUCIA框架开发scDeepLUCIA,调整表观基因组状态模块以学习染色质可及性,并整合低分辨率染色质接触图。在 down sampled bulk HiCAR数据上评估AUROC,并与HiCCUPS、Peakachu、CovNorm等比较。

5

增强子-基因对注释与功能验证

利用scHiCAR数据注释细胞类型特异性增强子-基因对,并评估其功能相关性

基于scHiCAR环支持的CRE-基因对,计算22种脑细胞类型中 distal cCRE可及性与mRNA水平的相关性,鉴定20,270个高置信度增强子-基因对,并与VISTA验证的增强子、H3K27ac数据等比较。

6

骨骼肌再生中的三模态单细胞动态分析

在单细胞水平解析肌肉干细胞再生过程中的基因调控动态,包括3D基因组重组

对未损伤和BaCl2诱导损伤后第5、7天的小鼠骨骼肌进行scHiCAR分析,获得41,578个细胞核,基于RNA、ATAC和染色质接触数据鉴定8种主要细胞类型,并对19,255个Pax7+ MuSCs及其后代进行拟时间分析。

研究方法

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

产品清单

实验环节名称品牌货号
基质胶Corning354230
mTeSR Plus培养基STEMCELL Technologies#05825
RPMI 1640培养基----
胎牛血清----
组织研磨器VWR62400-595
二硫苏糖醇SigmaD9779-250MG
无EDTA蛋白酶抑制剂Sigma5056489001
SUPERase RNA酶抑制剂InvitrogenAM2694
RNaseOUT重组RNA酶抑制剂Invitrogen10777019
CellTrics过滤器(30微米)Sysmex04-004-2326
水平转子离心机----
37%甲醛----
2.5 M甘氨酸----
条形码寡核苷酸(cHiC-R12-#1~#48)----
Tn5Merev寡核苷酸----
裸Tn5转座酶----
条形码oligo dT引物(cHiC-RNA-#1~#48)----
条形码随机六聚体(cHiC-RNA-Hex#1~#48)----
条形码连接接头(cHiC-R22-#1~#96)----
连接寡核苷酸(cHiC-linker0102-30bp-2)----
NPB缓冲液(含5% BSA、0.2% IGEPAL-CA630、1 mM DTT、1X无EDTA蛋白酶抑制剂、0.4 U/μl RNaseOUT的PBS)----
牛血清白蛋白SigmaA4503-100G
IGEPAL-CA630Alfa AesarJ61055
二硫苏糖醇Apex16-204
1X转座缓冲液(含33 mM Tris-乙酸、66 mM乙酸钾、10 mM乙酸镁、16% DMF)----
乙酸钾ThermoFisherP171-500
乙酸镁ThermoFisherBP215-500
台盼蓝----
涡旋振荡器(ika MS 3)Agilent--
Maxima H Minus逆转录酶thermoEP0753
T4 DNA连接酶NEBM0202L
DRAQ7染料cell signaling7406
Sony SH800分选仪Sony--
CviQI限制性内切酶----
夹板寡核苷酸(bind-rev-3)----
糖原Invitrogen10814010
NlaIII限制性内切酶NEBR0125L
DNA清洁与浓缩试剂盒ZymoD4013
Q5聚合酶NEBm0491L
KAPA高保真DNA聚合酶roche9420398001
DNA结合缓冲液zymoD4004-1-L
DNA回收试剂盒Zymo#D4002
DAPI染料----
木瓜蛋白酶----
Miltenyi髓鞘去除磁珠Miltenyi--
Miltenyi抗CD11b磁珠Miltenyi--
Miltenyi MACS抗ACSA-2微珠试剂盒Miltenyi--
MseI限制性内切酶----
T4 DNA连接酶----
NlaIII限制性内切酶----

关键环节

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

环节核对要点
scHiCAR细胞系验证
细胞类型:H1 hESC、GM12878、mESC;每个细胞类型12个孔;BC1条形码分配;测序深度:深度测序
阅读提示:Methods 'scHiCAR Protocol'及Figure 1b图注
小鼠额叶皮层scHiCAR
样本:小鼠额叶皮层;细胞数:5,313个细胞(两个生物学重复);测序:深度测序;参考数据集:BICCN snATAC-seq、bulk HiCAR
阅读提示:Results 'scHiCAR captures high quality mRNA, open chromatin and 3D genome contacts from mouse brain cells'及Figure 2图注
大规模小鼠脑细胞scHiCAR
细胞数:1.62百万个细胞核;批次:四个批次实验;测序:浅测序(~1,758接触/细胞)
阅读提示:Results 'scHiCAR accurately identifies the cell-type-specific transcriptome, accessible cCREs, and key 3D genome signatures in 22 cell types from 1.62 million mouse brain frontal cortex cells'及Figure 3图注
scDeepLUCIA训练与测试
训练数据:深度测序的bulk HiCAR数据集(小鼠和人类细胞系或原代细胞);测试数据:down sampled小鼠星形胶质细胞bulk HiCAR(30M reads);分辨率:5kb、10kb、50kb;评估指标:AUROC
阅读提示:Methods 'scDeepLUCIA'及Supplementary Methods;Figure 4图注和Extended Data Fig. 5
骨骼肌再生模型scHiCAR
模型:成年小鼠骨骼肌BaCl2诱导急性损伤;时间点:0、5、7天 post-injury;每个时间点两只小鼠;细胞数:41,578个细胞核;测序:深度测序
阅读提示:Methods 'Nuclei isolation from mouse frontal cortex and muscle for scHiCAR'及Figure 6a图注
增强子-基因对注释
增强子-基因对:20,270对(PCC FDR < 0.05);验证数据集:VISTA Enhancer Browser(182个验证增强子)、Paired-tag H3K27ac数据
阅读提示:Results 'scHiCAR accurately identifies functional enhancers and annotates their distal target genes'及Figure 5图注