急性髓系白血病單細胞層次的研究

題目:A single-cell survey of cellular hierarchy in acute myeloid leukemia
發(fā)表期刊:Journal of Hematology & Oncology (IF:23.2)
發(fā)表年份:2020-09

摘要

Background:急性髓系白血病 (AML) 是一種致命的造血惡性腫瘤螟凭,其預后與遺傳復雜性相關(guān)。目前尚未有從單細胞水平對不同 AML 亞型進行綜合分析掂榔。
Methods:使用Microwell-seq平臺弯予,分析了來自40名患者和3名健康供體的骨髓樣本的單細胞轉(zhuǎn)錄組測序。還使用single-cell single-molecule real-time (SMRT)測序來研究AML細胞的克隆異質(zhì)性。
Results:通過對191727個AML細胞的數(shù)據(jù)分析籍滴,建立了一個單細胞層面的AML“景觀”,并鑒定了一群具有新型AML marker的AML祖細胞群體,推斷出兩種具有不同臨床結(jié)局的AML亞型榴啸。聯(lián)合 Microwell-seq 和 SMRT sequencing探究AML中的線粒體突變孽惰。且提出存在一種表型“cancer attractor”, 這可能有助于定義AML祖細胞具有的共同表型。最后插掂,通過比較AML landscape和Human Cell Landscape來探索潛在的藥物靶點灰瞻。
Conclusion:他們鑒定出了一群的AML祖細胞群腥例。核糖體蛋白基因水平高表明提示病人預后不良辅甥。還推斷出具有不同臨床結(jié)局的AML亞型,并探索了潛在的藥物靶點燎竖。這些結(jié)果表明AML存在cancer attractor璃弄。

實驗設計及質(zhì)控

樣本: 40 病人和3個健康對照人的骨髓單核細胞(BMMC,Bone marrow mononuclear cell)
細胞:

結(jié)果

1.正常人BMMC單細胞分析

使用了Microwell-seq對三個健康供體BMMC進行分析(Figure 1A and B).
Figure 1A and B

鑒定了3種主要的細胞群體:淋巴細胞、紅細胞和髓系細胞(Figure1 A 和C).
Figure 1C
中性粒細胞再細分為6類細胞亞群, neutrophil A, neutrophil B, neutrophil C, neutrophil Dneutrophil Eneutrophil F 和neutrophil G(Figure S2).
Figure S2
為了進行l(wèi)ineage trajectory分析构回,從既往的研究中(數(shù)據(jù)ID號:GSE92274)整合了額外的2000個造血干/祖細胞(HSPC, hematopoietic stem/progenitor cells)和2719個外周血單核細胞(PBMC, peripheral blood mononuclear cells)夏块,總共獲得13280個健康對照人的細胞,使用partition-based graph abstraction (PAGA),構(gòu)建健康人的造血分化軌跡景觀(Figure 1D,E and F)
Figure 1D,E and F

2.鑒定初發(fā)AML的造血細胞群體

對40個AML樣本和3個健康對照樣本進行批次矯正后纤掸,得到20個亞群脐供。中性粒細胞,單核細胞借跪,紅系細胞和淋巴細胞都可在AML樣本和正常樣本中找到(Figure 2A and B).但是中央有一團細胞無特異性的基因表達(Figure 2C).
Figure 2A, B and C

網(wǎng)絡圖提示這群團細胞與髓系細胞很靠近政己,相關(guān)性圖提示這群基因與HSPC相關(guān),因此這群細胞定義為AML progenitor cells.
Figure 2D and E

3.AML progenitor cells的特征

與髓系細胞相比掏愁,AML progenitor cells和HSPC有許多共同的上調(diào)基因歇由,尤其是核糖體蛋白基因(RP,ribosomal protein genes,Figure 3A).通路富集分析顯示AML祖細胞具有較高的RB基因轉(zhuǎn)錄活性(圖卵牍。(Figure 3B).
Figure 3A and B

盡管存在RP基因表達的相似性,但是一些基因在HSPC沦泌,糊昙,AML progenitor cells和髓系細胞中的表達差異還是比較大的。

基因表達模式顯示AML progenitor cells處于從HSPC到分化的髓系細胞的“中間狀態(tài)”谢谦。具體來說释牺,AML progenitor cells中的GATA2、SPI1(和MPO的表達水平處于中間水平(Figure 3D)他宛。然而船侧,一些參與primitive state 的基因在AML progenitor cells表達最高,如SOX4厅各、FOS和ITM2A (Figure 3E).此外镜撩,CD99、CFD队塘、RACK1袁梗、FTL、B2M和ADA在AML progenitor cells中表達量較高憔古,先前已經(jīng)有研究發(fā)現(xiàn)這些基因與AML存在關(guān)系(Figure 3F).而且有基因一些與實體瘤相關(guān)的基因遮怜,這些基因與AML的關(guān)系目前尚無報道,例如TMSB10鸿市、SH3BGRL3锯梁、MGST1、MRPL33和MARCKSL1基因(Figure 3G)
Figure 3D, E and F
接著用VIPER(Virtual Inference of Protein-activity by Enriched Regulon)分析焰情,以揭示AML progenitor cells的轉(zhuǎn)錄因子(TF)活性陌凳。與HSPC相比,髓系細胞中差異表達的TF因子包括SPI1内舟、KLF5和JDP2合敦,而這些基因都參與造血發(fā)育的過程(Figure S7B).(F
Figure S7B
然而,參與惡性腫瘤的TF验游,如SMARCC1充岛、HOXA9和HOXB3,在AML progenitor cells是活躍的(Figure 3H).
igure 3H

4.AML progenitor cells的瘤內(nèi)異質(zhì)性可預測預后

AML progenitor cells比例在病人異質(zhì)性很大(Figure S5A).
Figure S5A

將AML progenitor cells分為16個亞群耕蝉,并將其歸納為4大類崔梗。Clusters1、2垒在、3蒜魄、4和11為RP gene high clusters。其余3種包括neutrophil-like, monocyte-like, and myeloid cell-like clusters (Figure 4A and B)
Figure 4A and B
一般來說,這種聚類方法與FAB(French–American–British)分類一致.但是基于單細胞轉(zhuǎn)錄組數(shù)據(jù)揭示AML progenitor cells的功能狀態(tài)會更詳細权悟。例如砸王,M5患者可繼續(xù)分為4個clusters,C5峦阁、C9谦铃、C12和C1(Figure S8A).
Figure S8A

富集分析表明,在RP gene high clusters中榔昔,C1和C2在調(diào)節(jié)肌動蛋白細胞骨架中起作用驹闰。C4與P53通路有關(guān)。C11富含癌蛋白聚糖撒会。在neutrophil-like種嘹朗,C5和C15與IL-17信號通路和癌癥翻譯失調(diào)有關(guān)。C7和C8富含MYC通路(Figure 4C).
Figure 4C
Figure S8D

使用WGCNA探索AML progenitor cells的基因調(diào)控網(wǎng)絡诵肛,鑒定了18個模塊(Figure S8E)屹培。其中,一個模塊由RP基因主導(Figure S8F).
Figure S8E and F
接著展示了這個模塊的基因網(wǎng)絡怔檩。該模塊中的所有RP基因都可以在RP gene high clusters中找到(Figure 4D)
Figure 4D
該模塊中的基因與各種腫瘤通路密切相關(guān)褪秀,提示這些基因參與AML的進展(Figure 4E).
Figure 4E

5.AML progenitor cells對于從AML的診斷和復發(fā)的臨床意義

接著觀察治療前后的細胞比例變化,將N02作為對照薛训,在緩解樣本中(P-extra 1隊列)媒吗,治療后的大多數(shù)細胞與正常細胞的重疊,其中AML progenitor cells和增殖細胞的數(shù)量減少乙埃,而中性粒細胞增加(Figure 5A-B). 髓系亞群基因在正常細胞中高表達闸英,而AML progenitor cells 基因主要在P-extra 1隊列治療前的隊列表達(Figure 5J). 在緩解的樣本中,RP和primitive genes在治療前表達量高介袜,治療之后卻比較低甫何,中性粒細胞基因則表現(xiàn)出相反的趨勢(Figure 5J).在復發(fā)樣本中(P20 隊列),骨髓中AML progenitor cells 占主導地位米酬,RP和primitive genes在治療前和治療后表達量均高(Figure 5K).
Figure 5A-F
Figure 5J-L

6.單核細胞白血病的瘤內(nèi)異質(zhì)性

然后沛豌,將分析重點放在一種亞型AML-M5上趋箩,并整合了15種單核細胞白血病的數(shù)據(jù)赃额。在此研究中,P06叫确、P14和P16是難治型AML患者跳芳;其余12例為非難治性患者。
AML祖細胞在患者中的比例差異很大(Figure S11A)竹勉。隨機選擇了12000個細胞難治型樣本(C5)和非難治性患者(C13)分析(Figure S11B-E)飞盆。GSEA表明與C13相比,MYC、SRC吓歇、RELA孽水、原癌基因MTOR相關(guān)通路在C5中過表達,表明C13惡性程度更高(Figure S11F)城看。TCGA-AML中女气,SRC和MTOR的高表達水平提示預后較差(FigureS11G)。GO分析揭示了
C5中的TF测柠,如CEBPG炼鞠、GATA2、MAX和
JUN基因集在C5中處于激活狀態(tài)轰胁,而CEBPE谒主、GATA1、MAZ赃阀、费尽,

和JUND基因集在C13處于激活狀態(tài)(FigureS11 H and I).
Figure S11

7.單細胞水平SMRT測序與遺傳突變的關(guān)系

由于覆蓋范圍的限制,3′端的單細胞轉(zhuǎn)錄組分析未能識別基因突變磕诊。因此扼脐,將 long-read sequencing與Microwell seq實驗相結(jié)合。我們使用靶向上游引物和通用下游引物從單細胞cDNA中擴增特定基因進行SMRT測序肖抱,可以獲得mutation sites and cell barcodes(Figure 6a).既往研究報道备典,mtDNA突變可提供 innate and natural barcodes來推斷細胞的克隆進化。他們選擇線粒體基因MT-ND4意述,使用SMRT測序進行譜系追蹤提佣。從P25、P10和P17樣本的單個細胞中鑒定到了MT-ND4的異質(zhì)變體荤崇。過濾后拌屏,以MT-ND5轉(zhuǎn)錄本為參考,得到平均3406個讀數(shù)术荤。突變位點的平均測序深度為3041倚喂。其中,UMAP聚類后瓣戚,每個樣本約有711個read端圈。接著將這些細胞分為幾個亞克隆(Figure S12A-C)子库。采用monocle3包舱权,描繪了髓系、紅系和淋巴系的分化軌跡(Figure 6B and Figure S12D and E)仑嗅。這些亞克隆與細胞表型的相關(guān)性很差宴倍。在不同的細胞亞群中檢測到所有MT突變的亞克隆张症,這表明多個細胞克隆可能有助于AML分級,不同的細胞可能在AML中屬于相似類型的“attractor”(Figure 6C and Figure S12 D-E).

接著使用ARACNe和VIPER的算法來進行基因網(wǎng)絡分析鸵贬。Figure 6D and E顯示了在P25患者中俗他,與HSPC細胞相比,AML祖細胞阔逼、AML(E)和正常(f)髓系細胞的attractor網(wǎng)絡核心基因拯辙。MYB、RUNX2和YY1是P25的AML祖細胞獨有的attractor颜价,可能是誘導癌癥的決定因素涯保。P10和P17也有自己的attractor網(wǎng)絡.(Figure S12 F-I).(
Figure 6
Figure S12

8.結(jié)合HCL 數(shù)據(jù)庫探究AML靶點

分析AML landscape和Human Cell Landscape兩者的數(shù)據(jù)后,差異分析得到了AML中高表達基因周伦。在前10個基因中夕春,F(xiàn)LT3是眾所周知的,POU4F1也有報道专挪。其中有一半基因是lncRNA及志。其他是PRSS21、CCL1和DNTT(Figure 7A)寨腔。相關(guān)分析顯示速侈,這前10個基因與AML中最相關(guān)的基因一起屬于兩個網(wǎng)絡(Figure 7B)。RP基因在這兩個網(wǎng)絡中都起著重要作用迫卢。接著倚搬,研究了AML祖細胞群中的前5個相關(guān)蛋白編碼基因,這些基因與腫瘤或髓系分化相關(guān)(Figure 7C).
Figure 7

結(jié)論

研究中使用Microwell-seq分析了來自40名AML患者的191727細胞乾蛤,以建立單細胞層面的AML景觀每界。鑒定出具有新型AML標志物的惡性AML祖細胞亞群。祖細胞RP基因高表達的患者的緩解率較低家卖。推斷出兩種具有不同臨床結(jié)局的AML眨层。我們表明,單細胞分析可用于預測和評估AML治療效果上荡。最后趴樱,通過將Microwell-seq與SMRT測序相結(jié)合,追蹤了AML景觀中的線粒體突變酪捡,并證明了AML克隆與轉(zhuǎn)錄組表型之間缺乏關(guān)聯(lián)叁征。總體研究表明“cancer attractor”的存在可能有助于定義AML祖細胞的共同表型沛善。

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末航揉,一起剝皮案震驚了整個濱河市塞祈,隨后出現(xiàn)的幾起案子金刁,更是在濱河造成了極大的恐慌帅涂,老刑警劉巖,帶你破解...
    沈念sama閱讀 206,013評論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件尤蛮,死亡現(xiàn)場離奇詭異媳友,居然都是意外死亡,警方通過查閱死者的電腦和手機产捞,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,205評論 2 382
  • 文/潘曉璐 我一進店門醇锚,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人坯临,你說我怎么就攤上這事焊唬。” “怎么了看靠?”我有些...
    開封第一講書人閱讀 152,370評論 0 342
  • 文/不壞的土叔 我叫張陵赶促,是天一觀的道長。 經(jīng)常有香客問我挟炬,道長鸥滨,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 55,168評論 1 278
  • 正文 為了忘掉前任谤祖,我火速辦了婚禮婿滓,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘粥喜。我一直安慰自己凸主,他們只是感情好,可當我...
    茶點故事閱讀 64,153評論 5 371
  • 文/花漫 我一把揭開白布额湘。 她就那樣靜靜地躺著秕铛,像睡著了一般。 火紅的嫁衣襯著肌膚如雪缩挑。 梳的紋絲不亂的頭發(fā)上但两,一...
    開封第一講書人閱讀 48,954評論 1 283
  • 那天,我揣著相機與錄音供置,去河邊找鬼谨湘。 笑死,一個胖子當著我的面吹牛芥丧,可吹牛的內(nèi)容都是我干的紧阔。 我是一名探鬼主播,決...
    沈念sama閱讀 38,271評論 3 399
  • 文/蒼蘭香墨 我猛地睜開眼续担,長吁一口氣:“原來是場噩夢啊……” “哼擅耽!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起物遇,我...
    開封第一講書人閱讀 36,916評論 0 259
  • 序言:老撾萬榮一對情侶失蹤乖仇,失蹤者是張志新(化名)和其女友劉穎憾儒,沒想到半個月后,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體乃沙,經(jīng)...
    沈念sama閱讀 43,382評論 1 300
  • 正文 獨居荒郊野嶺守林人離奇死亡起趾,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 35,877評論 2 323
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現(xiàn)自己被綠了警儒。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片训裆。...
    茶點故事閱讀 37,989評論 1 333
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖蜀铲,靈堂內(nèi)的尸體忽然破棺而出边琉,到底是詐尸還是另有隱情,我是刑警寧澤记劝,帶...
    沈念sama閱讀 33,624評論 4 322
  • 正文 年R本政府宣布艺骂,位于F島的核電站,受9級特大地震影響隆夯,放射性物質(zhì)發(fā)生泄漏钳恕。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點故事閱讀 39,209評論 3 307
  • 文/蒙蒙 一蹄衷、第九天 我趴在偏房一處隱蔽的房頂上張望忧额。 院中可真熱鬧,春花似錦愧口、人聲如沸睦番。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,199評論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽托嚣。三九已至,卻和暖如春厚骗,著一層夾襖步出監(jiān)牢的瞬間示启,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 31,418評論 1 260
  • 我被黑心中介騙來泰國打工领舰, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留夫嗓,地道東北人。 一個月前我還...
    沈念sama閱讀 45,401評論 2 352
  • 正文 我出身青樓冲秽,卻偏偏與公主長得像舍咖,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子锉桑,可洞房花燭夜當晚...
    茶點故事閱讀 42,700評論 2 345

推薦閱讀更多精彩內(nèi)容