文章标题:《A transcriptomic map of murine and human alopecia areata》
发表日期和杂志:2020年发表在JCI insight上
在线阅读链接:https://doi.org/10.1172%2Fjci.insight.137424
疾病介绍
斑秃(AA)是一种常见的自身免疫性疾病,最初表现为脱发,没有其他明显的皮肤变化。皮肤表面的不起眼的外观与毛囊中复杂的免疫活动形成了鲜明的对比。AA的发病机制是由于毛囊失去免疫豁免权,导致自身免疫攻击。
Alopecia areata可以影响任何年龄段的人,包括儿童和成年人。脱发通常以圆形或椭圆形的斑块开始,可能会逐渐扩大或合并成更大的区域。在某些情况下,脱发可能会进一步发展,导致全身性脱发,包括头发、眉毛、睫毛和体毛。
单细胞实验设计
使用单细胞测序来揭示小鼠再生障碍性贫血中免疫细胞的不同表达谱,小鼠再生障碍性贫血数据被用来生成对人类再生障碍性贫血的高度预测模型。
数据链接是:https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE145095 数据情况:
#samples
GSM4306442 Mouse_AA_skin
GSM4306443 Mouse_AA_lymphnode
GSM4306444 Mouse_Con_skin
GSM4306445 Mouse_Con_lymphnode
GSM4306446 Human_AA_skin
GSM4745080 Human_con_skin
提供的是10X格式的标准三个文件,选择下载数据之后需要对数据进行整理,将三个文件分别整理到对应的文件夹中。
#整理文件
fs=list.files('./','features')
fs
samples1= gsub('.tsv.gz','',gsub('features.','',fs))
samples1
samples2 = samples1
lapply(1:length(samples2), function(i){
x=samples2[i]
y=fs[i]
dir.create(x,recursive = T)
file.copy(from= y ,
to=file.path(x, 'features.tsv.gz' ))
file.copy(from= gsub('features','matrix',gsub('tsv','mtx',y)),
to= file.path(x, 'matrix.mtx.gz' ) )
file.copy(from= gsub('features','barcodes',y),
to= file.path(x, 'barcodes.tsv.gz' ))
})
然后使用Read10X函数将数据读取进来即可进行后续的标准分析
#指定数据存放位置
samples=list.files("./GSE145095_RAW-Human/outputs/")
samples
dir <- file.path('./GSE145095_RAW-Human/outputs/',samples)
names(dir) <- samples
#读取数据创建Seurat对象
counts <- Read10X(data.dir = dir)
sce.all = CreateSeuratObject(counts,
min.cells = 5,
min.features = 300 )
dim(sce.all) #查看基因数和细胞总数
as.data.frame(sce.all@assays$RNA$counts[1:10, 1:2])
table(sce.all@meta.data$orig.ident) #查看每个样本的细胞数
head(sce.all@meta.data)
对读取进来的数据进行质控、harmony整合以及单细胞细分亚群定义等。
因为人和鼠的数据给的都是10X的标准文件,所以按照一样的方法读取之后再处理即可!
为了阐明小鼠再生障碍性贫血的免疫细胞组成和转录异质性,从6332例不稳定性再障和4173例再生障碍性贫血的皮肤和淋巴结中分离和测序共10,505个免疫细胞
使用UMAP可视化,并鉴定了15个免疫细胞簇
其中大多数细胞由淋巴结(集群0、3、4、5、10和12)、皮肤(集群2、8、9、13和14)或混合(集群1、6、7、11)细胞组成
利用每个簇的中位基因表达量,通过2种方法将每个簇分配到一个细胞谱系中:
(a)免疫基因组计划获得的小鼠纯细胞基因特征的相关性;
(b) T细胞(Cd3d、Cd28、Cd4、Cd8)、APCs (Itgax [CD11c]、Itgam [CD11b]、Xcr1、Cd207 [Langerin])和b细胞(Cd79a、Cd19)的典型标记表达模式分析。
通过BioinfoArk提供的中国区chatGPT查询到各个细分亚群的解释:
抗原呈递细胞(APCs)细胞亚群细分——小鼠再生障碍性贫血的APC向CD11b+DC极化,具有活跃的促炎信号。
第8和第13簇由714个APC组成,包括276个AA和438个UA细胞
使用来自MSigDB文库(21)和先前获得的髓系特征(22)的基因特征进行了单样本基因集富集分析(ssGSEA)。基于富集分数进行主成分分析(PCA),发现AA和UA淋巴结的apc有明显的重叠(灰色椭圆)
除了细胞类型分化外,ssGSEA还显示小鼠AA APCs中血管生成、CD40、IFN-γ、JAK/STAT和缺氧信号的显著增加,支持该群体在AA中的促炎特征。
由于目前缺乏现成的小鼠APC数据,因此使用了人类APC特征。重新分析了数据,以便根据不同簇的特征基因表达特征标记不同的簇。在校正簇间的细胞周期状态后,对APC进行重新聚类,并检查APC的典型标记(图2,D-F)
使用典型标记,6个新的小鼠APC集群标记如下:
APC成分在疾病状态和组织部位之间存在显著差异。簇M0、moDC2、M3和Lc4主要在皮肤中发现,而cDc1和Lc5主要在淋巴结中发现。在再生障碍性贫血中,有更大比例的moDC,而不稳定型小鼠皮肤APC以Trem1+巨噬细胞为主。
确定再生障碍性贫血小鼠皮肤和淋巴结细胞的T细胞分化
文章还研究了小鼠AA和UA样本中T细胞的基因表达谱。T细胞由PAN-T细胞标志物的表达和相应的TCR测序信息来定义。在UMAP图中最初确定的15个簇中,1、2、3、4、6和7簇被发现对应于T细胞
使用典型和差异T细胞标记来分配集群识别
集群3和4分别被定义为CD4+CCR7lo和CD4+CCR7+细胞,并且在小鼠AA细胞的另一个生物复制中发现了类似的淋巴结集群相关性(集群4)。集群7被定义为Treg群体,具有高水平的Foxp3, CD25 (Il2ra), Helios (Ikzf2), neuropilin (Nrp1), GITR (Tnfrsf18), Tigit和Ctla4
人再生障碍性贫血皮肤T淋巴细胞单细胞测序显示类似的基因表达动态和途径富集。
使用建模来翻译我们的小鼠再生障碍性贫血的发现,并对从健康对照组和再生障碍性贫血患者皮肤中分离的T淋巴细胞进行单细胞RNA测序。经处理后,共回收2416个细胞(AA n=1664,对照组n=752),形成9个不同的簇
在鉴定的聚类中,T2、T4、T5、T6和T7是AA为主的聚类,AA样本的相对贡献大于50%。
使用上述纯免疫细胞群和典型/功能性T细胞标记物的双重相关性分析方法,鉴定了CD4+ (T1、T3、T4、T5和T7)、CD8+ (T6和T8)和混合CD4+/CD8+ (T2和T9) T细胞群。
执行ssGSEA,以确定聚类中可能存在的功能差异,发现IL2/STAT5信号的增加与CD4+ (T5)和CD8+ (T6)的一个子集以及调节性T细胞(T7)相关,这与在小鼠数据集中的发现相一致。
50% AA和对照细胞的T8细胞毒性淋巴细胞和促炎基因集的富集程度最高,利用TCR序列信息,还检查了单细胞队列内的克隆型动态,发现在前100个克隆的库空间中,AA内的前10个重复克隆型占39.7%,而在对照T细胞中这一比例为24.4%
在所有皮肤T淋巴细胞中,通过克隆型组进行ssGSEA检测激活和信号通路,其中CD4+ TCR激活、IL2/STAT5信号传递和一般TCR信号传递基因集通过增加克隆型而增加。虽然在正常皮肤克隆型组中没有看到减少,但与人AA富集相比,正常皮肤基因集富集明显减少
使用单细胞测序来揭示小鼠AA中免疫细胞的不同表达谱。发现CD4+和CD8+ T细胞的克隆扩增,在不同的转录状态下具有共享的克隆型。小鼠AA数据用于生成人类AA病的高度预测模型、,对人AA中T细胞的单细胞测序重现了预测模型的克隆型发现和基因表达。