文章标题:《Entinostat decreases immune suppression to promote anti-tumor responses in a HER2+ breast tumor microenvironment》
发表日期和杂志:2022年发表在Cancer Immunol Research上
在线阅读链接:https://doi.org/10.1158%2F2326-6066.CIR-21-0170
Eninostat简介
通过BioinfoArk提供的中国区chatGPT查询到药物的作用:
Entinostat是一种组蛋白去乙酰化酶抑制剂,也被称为HDAC抑制剂。它被用作抗癌药物,主要用于治疗某些类型的癌症,如淋巴瘤和乳腺癌。
Entinostat的作用机制是通过抑制组蛋白去乙酰化酶(HDAC),从而增加组蛋白的乙酰化水平。这可以导致基因的表达模式发生改变,包括抑制肿瘤细胞的增殖、诱导细胞凋亡和促进免疫应答。
在HER2阳性乳腺肿瘤微环境中,Entinostat可以少免疫抑制,从而增强免疫系统对肿瘤的反应。它可以改善肿瘤微环境中的免疫细胞功能,增加抗肿瘤免疫细胞的浸润和活性,并促进肿瘤细胞的免疫识别和杀伤。
单细胞实验设计
对使用Eninostat和ICIS治疗的小鼠的HER2过表达的乳腺肿瘤进行了单细胞RNA测序,以全面描述TME内多种细胞类型的变化。
从20个完整肿瘤中提取RNA,分别为载体对照组(V)、恩替诺他处理组(E)、恩替诺他与抗pd -1 (EP)、恩替诺他与抗ctla -4 (EC)、恩替诺他与抗pd -1和抗ctla -4 (EPC), 5个实验组每组4个生物重复。
数据链接是:https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE166321
文章数据是包含了scRNA-seq以及bulk RNA-seq数据,我们只选择下载scRNA-seq进行分析
#samples
GSM5068516_wtile1_barcodes.tsv.gz 12.5 Kb
GSM5068516_wtile1_features.tsv.gz 221.0 Kb
GSM5068516_wtile1_matrix.mtx.gz 17.9 Mb
GSM5068517_wtile2_barcodes.tsv.gz 9.8 Kb
GSM5068517_wtile2_features.tsv.gz 221.0 Kb
GSM5068517_wtile2_matrix.mtx.gz 10.5 Mb
GSM5068518_wtilv1_barcodes.tsv.gz 6.9 Kb
GSM5068518_wtilv1_features.tsv.gz 212.7 Kb
GSM5068518_wtilv1_matrix.mtx.gz 15.3 Mb
GSM5068519_wtilv2_barcodes.tsv.gz 18.3 Kb
GSM5068519_wtilv2_features.tsv.gz 221.0 Kb
GSM5068519_wtilv2_matrix.mtx.gz 26.8 Mb
提供的是10X格式的标准三个文件,选择下载我们需要的scRNA数据之后需要对数据进行整理,将三个文件分别整理到对应的文件夹中。
#整理文件
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("./mice/outputs/")
samples
dir <- file.path('./mice/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整合以及单细胞细分亚群定义等。
通过对RNA表达谱进行聚类并评估规范和细胞类型标记物的表达,共将56,731个细胞注释为细胞类型。
将细胞分为5种细胞类型,其中包括24,798个癌细胞(Erbb2+, Cdh1+), 6,645个CAFs (Col12a1+, Mmp2+), 1,419个淋巴样细胞(Cd3e+), 13,301个单核/巨噬细胞(Adgre1+, Itgam+, Cd14+)和10,568个MDSCs (S100a8/9+, Itgam+, Cd14+)。
通过BioinfoArk提供的中国区chatGPT查询到各个细分亚群的解释:
集群之间的差异表达基因
使用CoGAPS非负矩阵分解进行了无监督学习,比较了是否有任何细胞水平的特征(称为模式)与CoGAPS中直接获得的单细胞集群及其各自排名最高的基因相对应。
观察到总共36种CoGAPS模式,可以区分免疫细胞、髓细胞、上皮/成纤维细胞。
Eninostat对肿瘤细胞的影响
对肿瘤细胞进行了差异表达和通路分析,比较了Eninostat治疗和赋形剂。
Eninostat+ICIS对乳腺TME内未成熟髓系细胞表型的调节
MDSCs和巨噬细胞中Cd274和Pdcd1lg2 (PD-L1/2检查点配体)的主要髓系TME特征和表达,所以将分析重点放在Eninostat治疗诱导的分子变化上
首先收集了用于鉴定特定髓细胞亚群的典型标记和差异表达基因的代表性列表
用仅在髓系细胞上执行的第二次CoGAPS实施获得的额外模式交叉验证了这些亚群,分析确定了18种髓系细胞模式,其中7种对应于髓系亚群,其余的对应于具有共同功能或技术噪音的更广泛的细胞组
最终分析显示5,069个M-MDSCs(模式14),5,499个G-MDSCs(模式5),368个树突状细胞(DC,模式11),2,471个单核细胞(模式16),430个增殖的髓系细胞(模式6),8,079个M1样(模式9)和1,953个M2样(模式3)TAM
为了评估联合治疗对这些髓系细胞内功能可塑性的影响,在TAM和MDSC簇中进行了伪时间分析。
由于G-MDSCs被认为比M-MDSCs抑制更少,更不可能进一步分化为更特化的髓系细胞,而且M1样巨噬细胞被认为有助于更强大的髓系抗肿瘤反应,这两种细胞表型的变化都表明TME抑制较少,这可能有助于ICIS治疗的增敏。