单细胞毛细胞型星形细胞瘤细胞中的促有丝分裂与祖细胞基因程序

Mitogenic and progenitor gene programmes in single pilocytic astrocytoma cells.

作者信息Zachary J Reitman, Brenton R Paolella, Guillaume Bergthold, Kristine Pelton, Sarah Becker, Robert Jones, Claire E Sinai, Hayley Malkin, Ying Huang, Leslie Grimmet, Zachary T Herbert, Yu Sun, Jessica L Weatherbee, John A Alberta, John F Daley, Orit Rozenblatt-Rosen, Alexandra L Condurat, Kenin Qian, Prasidda Khadka, Rosalind A Segal, Daphne Haas-Kogan, Mariella G Filbin, Mario L Suva, Aviv Regev, Charles D Stiles, Mark W Kieran, Liliana Goumnerova, Keith L Ligon, Alex K Shalek, Pratiti Bandopadhayay, Rameen Beroukhim
PMID31427603
发布时间2019-08
DOI10.1038/s41467-019-11493-2

实验完整度

高

含6例PA患者样本scRNA-seq、荧光免疫组化、RNA原位杂交、bulk表达反卷积及小鼠神经干细胞过表达BRAF的scRNA-seq功能验证

主要模型

6例毛细胞型星形细胞瘤患者肿瘤组织 表达KIAA1549-BRAF、BRAF-V600E或空载体的小鼠神经干细胞(mNSC) 151例儿童低级别胶质瘤bulk表达谱

重点核对

BRAF融合检测:采用SMART-seq2联合BRAF 3'寡核苷酸spike-in及针对KIAA1549:BRAF融合接头的qPCR,检测率由0–0.7%提高至29–54% 细胞分选:按A2B5免疫标记状态分选单细胞,A2B5+细胞中79%检出BRAF融合 质量控制:931个细胞通过质控(表达基因>1000且线粒体基因<4%),基因需在≥10个细胞中表达 亚群划分:共享最近邻聚类(k=150,前7个主成分)及NMF鉴定5个细胞簇 程序评分与相关性:用MSigDB C2细胞周期和衰老基因集计算评分,Spearman相关并以Bonferroni校正

摘要

毛细胞型星形细胞瘤(PA)是儿童最常见的脑肿瘤,是一种由单一BRAF重排驱动的低级别胶质瘤。本文中,我们对六个PA进行了scRNAseq,所用方法能够检测该重排。与更高级别胶质瘤相比,显著更高比例的PA癌细胞表现出分化的、星形细胞样表型。较小比例的细胞表现出具有增殖证据的祖细胞样表型。这些细胞表达一种在更高级别胶质瘤中不存在的丝裂原活化蛋白激酶(MAPK)程序。免疫细胞,尤其是小胶质细胞,占PA所有细胞的40%,并解释了肿瘤位置和亚型之间bulk表达谱的差异。这些数据表明,MAPK信号仅限于PA中相对未分化的癌细胞,这对针对该通路的在研治疗具有意义。

实验结论

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

研究问题
PA癌细胞在单细胞水平存在何种基因表达程序与发育层级,MAPK通路是否在所有癌细胞中均匀激活,以及肿瘤相关免疫细胞如何影响bulk表达谱。
核心机制
PA癌细胞呈现从OPC样细胞到成熟星形细胞样细胞的发育分化层级,MAPK信号程序局限于相对未分化的祖细胞样亚群,并与增殖或衰老程序互斥;致癌性BRAF表达可抑制成熟胶质基因程序。
主要证据
6例PA的scRNA-seq鉴定出肿瘤簇与免疫细胞簇,肿瘤细胞富集OPC特征,AC样与MAPK程序分布于不同亚群;荧光免疫组化和RNA原位杂交验证OLIG2、GFAP、JUN、APOD的分区表达;表达KIAA1549-BRAF或BRAF-V600E的小鼠神经干细胞中AC样程序评分显著低于对照;反卷积分析显示小胶质细胞中位贡献约30%。
研究意义
为PA的生物学与治疗研究提供框架,并提出可解释MEK抑制剂反应异质性的临床假设供后续研究检验。

研究路径

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

1

建立PA癌细胞识别与单细胞测序流程

在低级别肿瘤缺乏大范围CNV和表达突变的情况下,将PA癌细胞与肿瘤相关细胞区分开。

采用快速解离、A2B5免疫标记与流式分选富集A2B5+胶质祖细胞,结合SMART-seq2、BRAF 3'寡核苷酸spike-in、KIAA1549:BRAF融合接头qPCR及表达胶质肿瘤标志物聚类鉴定癌细胞。

2

单细胞转录组聚类与癌细胞/免疫细胞解卷积

确定癌细胞与非癌细胞的转录组分组,并量化各细胞类型对bulk表达谱的贡献。

对scRNA-seq数据进行t-SNE、PCA、共享最近邻聚类和NMF,比较A2B5与BRAF融合状态;用差异基因富集评分估计151例儿童低级别胶质瘤bulk表达谱中各细胞类型贡献。

3

鉴定PA癌细胞正常脑细胞类型相似性

确定PA癌细胞最接近的正常脑细胞类型。

将PA癌细胞与发育中脑、成人皮层、发育皮层等单细胞图谱的细胞类型特征进行富集评分比较,并与H3K27M高级别胶质瘤、IDH突变星形细胞瘤和少突胶质细胞瘤比较。

4

解析PA癌细胞亚群基因程序

在癌细胞内部鉴定不同基因程序及其细胞状态分布。

对癌细胞进行PCA,按PC1和PC2排序基因,定义MAPK信号、星形细胞样(AC-like)和少突胶质细胞样(OC-like)程序,计算程序评分,并用荧光免疫组化和RNA原位杂交验证标志基因。

5

推断癌细胞发育轨迹并分析增殖与衰老

检验PA癌细胞是否存在发育过程,并确定增殖细胞与衰老细胞的分布。

用Monocle v2.6.4的反向图嵌入构建细胞轨迹并计算拟时间,分析MAPK、AC样程序评分随拟时间变化;计算细胞周期和衰老程序评分与MAPK程序、BRAF表达的相关性。

6

验证致癌性BRAF对细胞状态的影响

检验致癌性BRAF表达是否影响PA癌细胞基因程序。

对表达KIAA1549-BRAF、BRAF-V600E或空载体的小鼠神经干细胞进行scRNA-seq,比较差异表达基因与PA来源基因程序的重叠及AC样程序评分。

研究方法

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

产品清单

实验环节名称品牌货号
GentleMACS C管Miltenyi Biotec--
脑肿瘤解离试剂盒Miltenyi Biotec--
抗A2B5-APC抗体(克隆105HB29)Miltenyi Biotec--
钙黄绿素蓝AMInvitrogenC1429
LIVE/DEAD可固定近红外死细胞染色试剂盒--L-34974
SH800S荧光激活细胞分选仪Sony Biotechnology--
MoFlo Astrios EQ细胞分选仪Beckman Coulter--
SmartSeq2 v4试剂盒Takara--
Nextera XT反应试剂盒----
Agilent TapeStation 2200Agilent--
D5k筛选胶带Agilent--
Illumina NextSeq500测序仪Illumina--
Illumina HiSeq2500测序仪Illumina--
Agilent SureSelect杂交捕获试剂盒Agilent--
RNAscope人APOD探针ACD Biotech445171
RNAscope人JUN探针ACD Biotech470541
Olig2抗体R&D SystemAF2418
GFAP抗体(克隆GA5)Cell Signaling Technology3670
Ki-67抗体(克隆MIB-1)DakoM7240
Perkin Elmer Vectra 3成像系统PerkinElmer, Inc.--
DMEM/F12培养基----
无维生素A的B27添加剂----
表皮生长因子----
neurobasal-A培养基Life Technologies--
层粘连蛋白----
红细胞裂解液ThermoFisher--

关键环节

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

环节核对要点
肿瘤组织获取与解离
样本来源与伦理方案:波士顿儿童医院Dana Farber/Harvard Cancer Center IRB方案10–417;解离酶为胰蛋白酶或木瓜蛋白酶;6例肿瘤类型(5例经典KIAA1549-BRAF 16:9易位,1例非经典BRAF重复)
阅读提示:Methods第1段Tumor acquisition and preparation;Supplementary Data 1和2
A2B5免疫标记与单细胞分选
抗A2B5-APC用量10 µL/100 µL、Calcein Blue AM 1.5 µL/100 µL、LIVE/DEAD试剂0.5 µL;按A2B5状态大致等比例分选至96孔板;阳性对照为小鼠神经干细胞
阅读提示:Methods中Fluorescence-activated cell sorting小节
BRAF融合检测
BRAF 3'寡核苷酸终浓度13 µM;Taqman qPCR针对15–9、16–9或16–11断点;BT827和BT618 cDNA耗尽未做qPCR,BT646因复杂重排不适用该qPCR
阅读提示:Methods中Library preparation and sequencing with BRAF spike-in及Detection of BRAF fusion from single cell cDNA libraries小节
数据质控与聚类
剔除表达基因<1000或线粒体基因>4%的细胞;剔除在<10个细胞中表达的基因;剔除14个与免疫细胞聚在一起的A2B5+细胞及表达>1.5 log-normalized GZMB的细胞;聚类k=150、前7个主成分
阅读提示:Methods中Data quality control, dimension reduction, and clustering小节
小鼠神经干细胞实验
细胞来源为胚胎第14天小鼠内侧神经节隆起;培养条件为DMEM/F12(50:50)+1×B27(无维生素A)+20 ng/mL EGF;BRAF突变细胞培养不含生长因子;质控后487个细胞
阅读提示:Methods中KIAA1549-BRAF expressing neural stem cells及Mouse neural stem cell scRNA-seq小节
荧光免疫组化与RNA原位杂交
每种肿瘤分析20个高倍视野(20×);Olig2、GFAP、Ki-67抗体克隆与货号;RNAscope探针货号;IHC阳性阈值由20张图像交叉复核确定
阅读提示:Methods中Fluorescence immunohistochemistry及RNA in situ staining小节;Fig. 6