阿尔茨海默病细胞类型特征在不同人群中的共有性

Cell-type signatures of Alzheimer's disease shared across population groups

作者信息Tain Luquez, Jonathan Algoo, Rebecca Chiu, Jason A Mares, Archana Yadav, Matti Lam, Pallavi Gaur, Xiaoying Lai, Dylan I Lee, Fahad Paryani, Rafe Batchelor, Irla Belli, Jordan Henry, Bat Hoter-Ishay, Courteney Mattison, Lindsey Starr, Tsering Lama, Berke Karaahmet, Wenqing Cao, Philip L De Jager, Mariko Taga, Lisa L Barnes, David X Marquez, David A Bennett, Ya Zhang, Vilas Menon
PMID42457956
期刊Nature
发布时间2026-09
DOI10.1038/s41586-026-10793-0

实验完整度

研究包含多人群(非裔美国人、拉丁裔、白人)死后脑组织的单核多组学(RNA+ATAC)分析,并用Xenium空间转录组验证了关键细胞比例关联,属于多证据层级综合。

主要模型

死后脑组织样本(背外侧前额叶皮层、颞上回、尾状核) 单核RNA测序数据(约190万核) 单核ATAC测序数据 Xenium原位杂交组织切片(33位供体)

重点核对

样本来源:ROSMAP、MARS和Latino Core研究队列的自报种族/族裔分组(非裔美国人、拉丁裔、白人) 脑区:背外侧前额叶皮层、颞上回、前尾状核 数据类型:snRNA-seq和snATAC-seq,共399个样本 关联分析:使用ANCOM-BC、quasibinomial logistic regression、LINDA和sccomp四种方法,meta分析FDR调整P<0.1,I2<50%,至少两种方法一致 验证:Xenium平台33个供体的组织验证,标记基因面板包含259个基因

摘要

单细胞分辨率的基因组研究已识别出几种与阿尔茨海默病临床和病理特征相关的细胞类型,但尚未检验在不同人群中共享的关联。为弥合这一差距,我们使用单核RNA测序和转座酶可及染色质测序,对拉丁裔、白人(不含拉丁裔)和非裔美国人(不含拉丁裔)个体的死后脑组织样本中的皮层和皮层下区域进行了分析。通过离散和连续的分子程序解析,我们识别出在所有三个种群组中,与阿尔茨海默病相关的细胞类型特异性簇,这些簇具有区域特异性,包括小胶质细胞(GPNMB+和CD74+亚群)、星形胶质细胞(SERPINH1+、CD44+和WIF1+亚群)和神经元(SST+ GABA能神经元和浅层谷氨酸能神经元)特征。我们还报告了星形胶质细胞和少突胶质细胞中的连续基因表达因子,这些因子未被离散簇分配捕获,但与疾病表型显示出强烈的关联;这些因子富集了与注释功能(如脂质处理和神经递质再摄取)相关的基因。最后,我们发现分子程序揭示了六个不同的认知障碍个体亚组,这些亚组跨越所有三个种群,未被神经病理学捕获,而是由并非普遍存在但与死前障碍相关的分子特征所区分。总的来说,我们的研究确定了与阿尔茨海默病相关的关键细胞类型和基因程序,这些在不同的种群组中是共享的,并强调了代表性采样如何捕获共享特征和疾病异质性,从而能够更好地优先考虑关键细胞类型以供进一步研究。

实验结论

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

研究问题
在非裔美国人、拉丁裔和白人群体中,哪些细胞类型特异性的分子特征与阿尔茨海默病的病理和认知表型共同关联?
核心机制
研究未明确建立分子机制,但识别出与阿尔茨海默病相关的细胞类型特异性的共享特征,包括小胶质细胞(GPNMB+,CD74+)、星形胶质细胞(SERPINH1+,CD44+,WIF1+)和神经元(SST+ GABAergic,浅层谷氨酸能)亚群,以及少突胶质细胞中的连续基因表达因子,这些因子富集于脂质处理和神经递质再摄取等功能。
主要证据
通过单核RNA测序和ATAC测序,对167名个体的死后脑组织样本进行分析,并结合Xenium原位杂交验证,发现GPNMB+小胶质细胞在痴呆和认知衰退个体中比例增加,与tau病理相关;并发现六个分子亚组与认知障碍相关,主要通过细胞比例和连续的基因因子分析得到。
研究意义
该研究强调代表性采样能够识别跨种群的共享分子特征,有助于优先考虑关键细胞类型进行进一步研究,并为针对多种族背景的阿尔茨海默病治疗靶点的识别奠定基础。

研究路径

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

1

样本收集与多组学数据生成

获取并处理来自不同种族/族裔个体的死后脑组织样本,生成单核RNA测序和ATAC测序数据,以覆盖多种细胞类型。

从ROSMAP、MARS和Latino核心研究队列中选取167名个体,解剖出三个脑区(DLPFC、STG、AC)的组织样本,提取细胞核后进行snRNA-seq和snATAC-seq建库和测序。

2

聚类与细胞类型注释

识别主要的细胞类型和亚群,并确定其在不同脑区和个体中的组成。

对snRNA-seq数据使用Seurat进行迭代聚类,首先区分主要细胞类别,然后对每个主要类别进行二次聚类以确定亚群。最终确定了61个亚群,包括神经元、胶质细胞和血管细胞等。

3

细胞比例与疾病表型的关联分析

评估不同细胞亚群的比例与阿尔茨海默病病理和认知表型之间的关联,并确定跨种族/族裔群体的共享特征。

对每个种簇,分别使用ANCOM-BC、quasibinomial logistic regression、LINDA和sccomp四种方法检验其比例与临床和病理变量的关联,然后进行元分析以确定跨群体的一致性。

4

连续表达因子分析

识别离散聚类未能捕获的、与疾病表型相关的连续基因表达模式,特别是少突胶质细胞和星形胶质细胞中的因子。

使用scHPF算法在每个主要细胞类别中提取连续表达因子,并检验因子得分与临床和病理变量的关联。

5

snATAC-seq差异可及性与驱动转录因子推测

探索染色质可及性差异,识别可能驱动聚类或因子身份的转录因子,并评估AD风险位点的开放性。

对snATAC-seq峰进行差异分析,比较不同转录组定义的亚群间的染色质开放性,并富集转录因子结合基序。在AD相关GWAS位点计算染色质开放性,比较不同群体间的差异。

6

完整组织切片的验证

使用Xenium原位杂交技术验证snRNA-seq中发现的细胞类型比例与疾病表型的关联。

设计包含259个基因的针对面板,对来自33个供体的死后DLPFC组织切片进行原位杂交。将snRNA-seq聚类标签映射到细胞,计算组织中的细胞比例,并使用ANCOM-BC检验其与临床和病理变量的关联。

7

外部数据集映射与共识评估

通过与以往两项主要包含白人个体的snRNA-seq研究(Mathys et al. 2023和Green et al. 2024)的结果进行比较,验证本研究的发现。

将本研究的聚类标签用MapMyCells工具映射到两个外部数据集,检验关键聚类(如GPNMB+小胶质细胞)的比例关联是否一致。同时,将外部聚类的标签映射回本研究,验证其关联。

8

个体分子亚组鉴定

在认知障碍个体中探索基于细胞组成和基因表达因子的分子亚组,并分析其与临床和病理特征的关系。

基于所有细胞亚群比例和因子得分,对MCI或AD个体的STG样本进行层次聚类,确定六个亚组,并评估其组成与临床、病理、种群和遗传特征的关联。

研究方法

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

产品清单

实验环节名称品牌货号
杜恩斯匀浆器SigmaD8938
RNA酶抑制剂Takara2313B
牛血清白蛋白NEBB9000S
预分离过滤器(20μm)Miltenyi130-101-812
Cellometer Vision细胞计数仪Nexcelom--
细胞计数板NexcelomCHT4-SD100-002
Chromium Next GEM单细胞Multiome ATAC+基因表达试剂盒10x GenomicsPN-1000283
Chromium Next GEM J芯片10x GenomicsPN-2000264
NovaSeq 6000测序仪Illumina--
CellRanger ARC v2.0.210x Genomics--
CellBender v0.2.0----
Seurat v.4----
popscle (freemuxlet)----
Signac----
Seurat v.5----
RCTD(稳健细胞类型分解)----
MapMyCellsAllen Institute for Brain Science--
Xenium分析仪----
DESeq2----
XSTREME----

关键环节

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

环节核对要点
研究参与者选择
参与者自报种族/族裔(非洲裔美国人、拉丁裔、白人),来源队列(ROSMAP、MARS、Latino Core),纳入标准为拥有至少两个脑区的冷冻组织且具有完整的病理表征
阅读提示:Methods 的 Study participants, ethics and clinicopathological characterization
脑区取样与组织处理
脑区(DLPFC、STG、AC),组织重量(约50-100mg灰质),裂解缓冲液成分,匀浆次数(25次A、25次B),离心条件(500g,5分钟),过滤孔径(20μm)
阅读提示:Methods 的 Generation of snRNA-seq and snATAC-seq (multiome) data 下的 Tissue preparation
单核多组学建库和测序
每批混合三个样本(来自不同脑区和种族/族裔),目标细胞数(约20,000个核),使用的试剂盒(Chromium Next GEM Single Cell Multiome ATAC+Gene Expression Reagent Bundle)和芯片(Chromium Next GEM Chip J),RT反应条件(45分钟37°C,30分钟25°C),测序平台和覆盖度(NovaSeq 6000,约645M reads/样本)
阅读提示:Methods 的 Generation of snRNA-seq and snATAC-seq (multiome) data 下的 Library preparation and sequencing
数据预处理和质量控制
比对和定量使用CellRanger ARC v2.0.2默认参数,使用CellBender去除背景,过滤标准(UMI > 500,线粒体UMI百分比 < 5%),聚类参数(PC数、分辨率),迭代聚类流程,双细胞去除方法(freemuxlet)
阅读提示:Methods 的 Primary processing of sequencing data, Demultiplexing of sample barcodes and assignment to donors, and Clustering of nuclei
差异丰度分析模型
统计方法(ANCOM-BC, quasibinomial, LINDA, sccomp),模型公式(比例 ~ 变量 + 年龄 + 性别 + 尸检间隔,认知变量另加教育),元分析方法(固定效应模型),显著性标准(FDR调整P < 0.1,方向一致,I2 < 50%,至少两种方法)
阅读提示:Methods 的 Associations of cell-type proportions with clinical and pathological traits 和 Meta-analysis across populations
连续因子分析
scHPF的因子数量选择标准,基因选择(权重最高前50个基因),因子得分计算(每个供体/区域的中位数),关联分析用线性模型,协变量同前,显著性标准同聚类分析
阅读提示:Methods 的 Identifying cell-type factors through scHPF
染色质可及性分析
snATAC-seq峰的计算和统一,差异峰分析使用DESeq2,亚群比较的协变量(批次、尸检间隔、年龄),AD风险位点的选择(GWAS队列),染色质开放性的归一化方法,群体间差异的统计检验(Kruskal-Wallis,事后Mann-Whitney)
阅读提示:Methods 的 snATAC-seq visualization and differential peak analysis 和 Open chromatin quantification for AD risk loci
空间验证实验
Xenium平台,基因面板包含259个基因,组织切片厚度(10μm),探针杂交条件(50°C过夜,浓度10 nM),连接条件(37°C 2小时),扩增条件(1小时4°C,2小时37°C),图像采集循环数(15轮),信号处理参数(默认),细胞分割参数(XOA v3.0默认)
阅读提示:Methods 的 Xenium multiplexed ISH experiments
个体分子亚组分析
聚类特征(细胞比例和因子得分),输入样本(MCI或AD供体的STG样本),层次聚类方法(距离1-Pearson相关,Ward最小方差),最佳簇数通过轮廓系数,评估稳定性方法(添加噪声和扰乱,计算Jaccard相似度)
阅读提示:Methods 的 Donor clustering based on cell subcluster frequencies and median factor scores
遗传分析
WGS数据使用条件,质量控制标准(变异缺失<2%,MAF>1%,HWE P>1e-6),LD剪除参数(r2>0.2,窗口200 SNP,步进50),PC计算,ADMIXTURE分析的k值选择(k=3),多基因风险评分计算使用Bellenguez et al.汇总统计和PLINK
阅读提示:Methods 的 Genetic ancestry estimation and modelling 和 Polygenic risk scoring