单细胞测序揭露小细胞肺癌化疗后瘤内异质性

当你的才华还撑不起你的野心时,请潜下心来,脚踏实地,跟着我们慢慢进步。不知不觉在单细胞转录组领域做知识分析也快两年了,通过文献速递这个栏目很幸运聚集了一些小伙伴携手共进,一起成长。

文献速递栏目通过简短介绍,扩充知识面,每天关注,希望你也能有所收获!

单细胞转录组分析综述

单细胞入门-读一篇scRNA-seq综述

文章信息

● 题目:Single-cell analyses reveal increased intratumoral heterogeneity after the onset of therapy resistance in small-cell lung cancer
● 发表日期:2020年2月17日
● 杂志:Nature Cancer
● DOI:https:/ / doi.org/10.1038/ s43018-019-0020-z

摘 要 

小细胞肺癌(small-cell lung cancer,SCLC)常常表现为化疗敏感到化疗耐药的快速演变,由于复发后组织样本的稀缺,这种转变的机制仍然不清楚。作者从 SCLC患者体内分离出循环肿瘤细胞(Circulating tumor cell ,CTC),以此构建CDX模型 (CTC-derived xenograft, CDX),然后通过对患者CTC,以及化疗耐药和化疗敏感患者CDX模型进行单细胞测序,以阐释SCLC瘤内异质性( intratumoral heterogeneity, ITH)。从中,其发现化疗耐药的细胞群体ITH普遍性增加,包括治疗靶点和潜在耐药途径的异质性表达。同样,直接从血液中提取患者CTC进行分析证实复发后ITH增加。这些发现表明,SCLC耐药的特点是细胞亚群之间的基因异质性的表达,导致肿瘤出现多重并发的耐药。这些发现强调了临床上需要对初发的SCLC进行合理的联合治疗,以最大限度地提高初始治疗效果,抑制ITH及多重耐药的出现。

样本处理

1

样本选择

化疗敏感性

MDA-SC39, MDA-SC68 , HCI-008, MDA-SC4化疗敏感;

MDA-SC16,MDA-SC49,MDA-SC55,MDA-SC75化疗不敏感。

一共构建了8个CDX模型

其中MDA-SC4、MDA-SC68 又构建了CDX耐药模型

单细胞测序取样点CDX传代次数:

MDA-SC4s: P3;  MDA-SC16r: P3;  MDA-SC39s: P2;  MDA-SC49r: P2;

MDA-SC55r: P1;  MDA-SC68s: P2;  MDA-SC75r: P1;  HCI-008: P3.

患者CTC取样:

MDA-SC55三个时间点: diagnosis, responding, relapsed;

1

流程图

所有单细胞及bulk样本数据储存在GSE138474中

1

数据产生及初步处理

  • 10x 测序平台v2试剂盒,HiSeq 4000 (Illumina)测序

  • 过滤、总细胞数量:CDX模型 Cells with fewer than 3,000 expressed genes and genes expressed in <10% cells were filtered out. The total numbers of cells in each model after filtering were as follows: MDA-SC4s: 8,998; MDA-SC16r: 7,173; MDA-SC39s: 2,053; MDA-SC49r: 2,593; MDA-SC55r: 2,853; MDA-SC68s: 2,779; MDA-SC75r: 3,959; HCI-008: 5,233. Each CDX model was downsampled to 2,000 for inter-model comparison.  单独CTCs测序数据:CTCs with fewer than 50 expressed genes and genes that were expressed in <1% cells were also filtered out.After filtering, the treatment-naïve, treatment responding and treatment relapsed samples had 637, 1,108 and 974 cells remaining, with the numbers of genes expressed ranging from 50–2,153.

  • 初步处理:Seurat V2.3流程:

    Normalize data

    Scale Data(regress out G2/M and S phase cell cycle gene scores)

    PCA(based on highly variable genes)

    Elbowplot选dims

    降维-RUNTSNE(default perplexity value of 30)

数据分析思路

1

CDX数据样本聚类及常见肿瘤相关marker展示

同一样本内肿瘤marker表达不一致如MDA-SC49r里MYC、EMTscores表达不一致,体现了瘤内异质性的存在。

1

CDX模型之间baseline ITH计算

  • 选取前20个PC,计算细胞之间的平均Euclidean distance来计算ITH score。(对应图a)

  • 通过dispersion of gene expression来衡量ITH:使用所有基因的normalized后的表达量,通过EdgeR包的estimateDisp()计算。(对应图b)

  • 基于tSNE,计算每个样本的 the optimum number of clusters(cluster 包 Silhouette score 最大值时),可发现铂耐药组细胞cluster较多,有明显异质性。

  • 差异表达基因表明铂敏感和耐铂CDX之间差异表达的基因只在单个CDX中检测到,而不是存在于在两个或更多耐药模型中。

  • GSEA分析揭示CDX模型各个cluster显著富集的通路,提示即使在一个肿瘤内部也有多个耐药相关信号通路共表达。

1

CDX模型内部各cluster之间存在基因差异表达

运用heatmap和violinplot展示与SCLC治疗和耐药相关的基因在一个样本内部的各cluster间存在差异表达(如PARP1, CHEK1, AURKA, AURKB, DLL3, VIM, 及EMT scores等)。

1

患者CTCs单细胞测序

分别从diagnosis, responding, relapsed三个时间点手机患者MDA-SC55血液中的CTC进行测序,质控降维聚类后发现两个cluster是CTCs。

提出CTCs,重新降维聚类,可得到5个不同的cluster,其中cluster1主要来自diagnosis时间点的CTC。

分别计算diagnosis时的CTC及relapsed时间点的CTC及CDX模型的ITH评分,发现relapsed时间点ITH评分更高。relapsed时间点CTC及CDX模型SCLC转移及耐药相关基因表达有所不同。

1

CDX模型模拟患者化疗耐药过程

对MDA-SC68s CDX模型分别进行对照及cisplatin处理,当肿瘤体积增至1000mm3后取肿瘤进行单细胞测序。图c可看出对照组及cisplatin耐药复发组肿瘤细胞几乎无混合。

与对照组相比,cisplatin-treated CDX ITH score增高

与对照组比,cisplatin-treated CDX里 ASCL1明显下调,与药物耐药显著相关。cluster3 EMT score升高, DLL3特异性降低。提示在MDA-SC68s中,经过铂处理后,间充质细胞取代了对铂相对敏感的ASCL1阳性细胞。

作者接着用SCLC的靶向治疗药物PARPi及CHKi进行了如上述类似的实验及分析,提出PARPi与CHKi耐药机制可能不同。

临床意义

  • 构建了可模拟SCLC疾病进展及转移的CDX模型;

  • 化疗耐药后肿瘤异质性增加。同一肿瘤内部存在转录异质性,体现多种共存的耐药机制。提示SCLC在治疗初期应采用联合化疗方法。

  • CTCs转录谱在同一病人不同时间段存在异质性表达,提示可通过病人病程中连续取样,从而捕捉患者治疗靶点表达的动态变化。

(0)

相关推荐