單細胞測序技術(shù)在乳腺癌中的研究進展

楊帥 王新恒 吳佳樂 付子彤 魏博 楊靈冰 綜述 許守平 審校

作者單位：哈爾濱醫(yī)科大學附屬腫瘤醫(yī)院乳腺外科一病區(qū)(哈爾濱 150081)

乳腺癌是女性中最常見的惡性腫瘤,2020年大約有230萬新發(fā)病例,超過68.5萬人死于乳腺癌[1]。乳腺癌是一種起源于乳腺導(dǎo)管上皮細胞高度異質(zhì)性的疾病,由于這些異質(zhì)性的存在,乳腺癌治療面臨很大的困難。單細胞測序技術(shù)是一種在單細胞水平上揭示細胞的基因組、轉(zhuǎn)錄組或表觀遺傳變化的技術(shù),從不同角度揭示細胞的功能和特征。通過構(gòu)建數(shù)據(jù)庫,可以一次檢測龐大的每個單元信息。人類對單細胞的研究主要是基于測序技術(shù)。第一代測序技術(shù)是1997年開發(fā)的Sanger測序方法,主要測序標記的單個核酸分子的堿基,科學家使用該技術(shù)在2003年繪制了人類基因的初步序列,為生命的解碼奠定了基礎(chǔ)。2005年,Roche公司開發(fā)了第二代測序(高通量平行測序),進一步解碼單細胞轉(zhuǎn)錄組、基因組、表觀遺傳學和結(jié)合多組學單細胞測序技術(shù)的細胞和分子變化,為醫(yī)學研究進展提供了良好的平臺。傳統(tǒng)高通量測序技術(shù)獲得的測序數(shù)據(jù)只是對大量組織樣本中混合細胞的平均基因測定。與傳統(tǒng)的測序技術(shù)相比,單細胞測序技術(shù)不僅能夠準確測量基因表達水平和檢測非編碼RNA的表達,還可以充分利用測序特殊樣本的優(yōu)勢,彌補特殊樣本樣本量小、不易獲得的缺點[2-4]。

近年來不斷發(fā)展的單細胞測序技術(shù)已經(jīng)能夠分離單細胞,在單細胞水平上建庫測序,捕獲單個細胞的基因表達信息。2009年,Tang等[5]首次報道了單細胞轉(zhuǎn)錄組測序技術(shù)。最初的單細胞測序只能對10～100個細胞進行研究,隨著技術(shù)的不斷完善,現(xiàn)階段人們已經(jīng)能夠在單個項目中對成千上萬個單細胞進行測序[6-8]。運用單細胞測序技術(shù),不僅能夠區(qū)分腫瘤細胞間的異質(zhì)性,也能識別不同細胞類型。本文綜述了單細胞測序技術(shù)應(yīng)用于乳腺癌圖譜、乳腺癌微環(huán)境、乳腺癌轉(zhuǎn)移與耐藥以及乳腺癌異質(zhì)性與克隆演化方面的文獻,為乳腺癌的研究提供參考。

1 單細胞測序技術(shù)應(yīng)用于乳腺癌圖譜

乳房由復(fù)雜的上皮導(dǎo)管和小葉網(wǎng)絡(luò)組成,這些上皮導(dǎo)管和小葉網(wǎng)絡(luò)嵌在結(jié)締組織和脂肪組織中[9-11]。目前大多數(shù)的研究都集中在乳腺上皮細胞類型上,而對非上皮細胞類型的研究相對較少[12-14]。Kumar等[15]構(gòu)建了具有單細胞和空間分辨率的綜合人類乳腺細胞圖譜?？蒲腥藛T通過單細胞轉(zhuǎn)錄組測序研究分析了146例女性的714 331個細胞和117 346個細胞核,得到了包括內(nèi)皮細胞、免疫細胞、基底肌上皮細胞、間充質(zhì)細胞和肥大細胞在內(nèi)的12種不同的細胞類型圖譜。這些正常乳腺細胞類型圖譜揭示了乳腺組織豐富的血管周圍樣細胞(Perivascular like, PVL)、內(nèi)皮細胞和免疫細胞群,以及高度多樣化的管腔上皮細胞群。同時這些數(shù)據(jù)為研究乳腺癌和乳腺生物學等疾病提供了正常乳腺組織的參考信息。

乳腺癌的異質(zhì)性可能源于其不同的細胞來源,體現(xiàn)在不同的亞型與獨特的治療不敏感性[16]。然而科研人員對早期腫瘤形成和進展的研究常常受到不發(fā)達、偶爾相互矛盾的乳腺譜系注釋的限制。現(xiàn)有的注釋通常由少量標記定義,不能充分代表乳腺譜系內(nèi)異質(zhì)性。乳房是一個動態(tài)器官,其對生理和病理生理條件的反應(yīng)會改變其疾病易感性,但這些臨床變量對細胞狀態(tài)的具體影響仍不清楚。Gray等[17]通過整合單細胞RNA測序(single-cell RNA sequencing,scRNA-seq)、質(zhì)量流式細胞術(shù)和循環(huán)免疫熒光,呈現(xiàn)一個統(tǒng)一的高分辨率乳腺癌圖譜。他們定義了肺泡、激素感應(yīng)和基底上皮譜系中的細胞亞型,描述了不同亞型與癌癥風險因素的關(guān)聯(lián),包括年齡、胎次和BRCA2種系突變。被稱為基底腔細胞的肺泡細胞子集,帶有與基底樣乳腺癌相關(guān)的基因特征。有趣的是肺泡細胞子集表現(xiàn)出較差的轉(zhuǎn)錄譜系保真度,并且隨著年齡的增長而積累。這些數(shù)據(jù)更全面地高分辨率定義乳腺癌細胞類型。

乳腺癌是復(fù)雜的細胞生態(tài)系統(tǒng),其中異型細胞相互作用在疾病進展和對治療的反應(yīng)中發(fā)揮核心作用。然而,目前的研究對其細胞組成和組織的了解是有限的。Wu等[18]提出了人類乳腺癌的單細胞和空間轉(zhuǎn)錄組學分析,并開發(fā)了一種內(nèi)在亞型分類的單細胞方法來揭示復(fù)發(fā)性腫瘤細胞的異質(zhì)性。通過scRNA-seq使用轉(zhuǎn)錄組和表位的細胞索引進行免疫表型分析可提供高分辨率的免疫圖譜,包括與臨床結(jié)果相關(guān)的細胞程序性死亡-配體1陽性或者細胞程序性死亡-配體2陽性的巨噬細胞群。他們發(fā)現(xiàn)間充質(zhì)細胞通過在內(nèi)皮細胞、癌癥相關(guān)成纖維細胞(Cancer-associated fibroblast, CAF)和PVL三個主要譜系內(nèi)分化顯示出不同的功能和細胞表面蛋白表達的差異。該研究提供了乳腺癌細胞結(jié)構(gòu)的全面轉(zhuǎn)錄圖譜。

2 單細胞測序技術(shù)應(yīng)用于乳腺癌微環(huán)境研究

腫瘤微環(huán)境是由腫瘤細胞、免疫細胞、成纖維細胞和血管細胞組成的復(fù)雜微環(huán)境。腫瘤內(nèi)的細胞組成決定了腫瘤微環(huán)境的高度異質(zhì)和動態(tài)特性[19]。乳腺癌的微環(huán)境包括腫瘤浸潤淋巴細胞、腫瘤相關(guān)骨髓細胞[20-21]、脂質(zhì)相關(guān)巨噬細胞[22]、癌癥相關(guān)成纖維細胞(Cancer-associated fibroblast,CAF)和其他促進腫瘤細胞生長和遷移的因子。微環(huán)境中的細胞成分被證明參與乳腺癌的發(fā)生、生長和轉(zhuǎn)移。目前,乳腺癌微環(huán)境細胞的模式研究相對較少,沒有很好的評估這些微環(huán)境細胞在發(fā)病機制中的具體作用 。單細胞多組學生物技術(shù),尤其是單細胞RNA測序技術(shù)以更高的分辨率揭示單細胞表達水平,精細剖析腫瘤微環(huán)境的分子特征[23]。

CAF大量存在于幾乎所有腫瘤微環(huán)境中,并影響腫瘤進展。雖然人們逐漸地了解到它們的功能和異質(zhì)性,但對CAF的起源知之甚少。Houthuijzen等[24]通過對小鼠進行的各種移植研究,發(fā)現(xiàn)浸潤性小葉乳腺癌和三陰性乳腺癌中的CAF均源自乳腺組織常駐的正常成纖維細胞(Normal fibroblasts,NF)。單細胞轉(zhuǎn)錄組學、體內(nèi)和體外研究揭示了CD26+和CD26-NF細胞群分別轉(zhuǎn)變?yōu)檠仔訡AF和肌成纖維細胞CAF。功能性共培養(yǎng)實驗表明,CD26+NF轉(zhuǎn)變?yōu)榇倌[瘤炎性CAF,后者以趨化因子12依賴性方式募集骨髓細胞,并通過基質(zhì)金屬蛋白酶活性增強腫瘤細胞侵襲。他們的結(jié)果表明CD26+和CD26-NF在乳腺癌小鼠模型中轉(zhuǎn)化為不同的CAF亞群。Bartoschek等[25]利用遺傳工程改造的乳腺癌小鼠模型的768例間充質(zhì)細胞轉(zhuǎn)錄組的單細胞RNA測序,定義了三個不同的CAF亞群。并且在轉(zhuǎn)錄和蛋白水平上揭示了不同來源的CAF亞類的空間分布,包括PVL、乳腺脂肪墊和轉(zhuǎn)化的上皮細胞。每種CAF亞型的基因譜與其獨特的功能程序相關(guān),CAF三種亞群的分辨率可作為生物標志物為CAF精確靶向藥物的開發(fā)提供了可能性。

3 單細胞測序技術(shù)應(yīng)用于乳腺癌轉(zhuǎn)移研究

作為臨床預(yù)后的指標,乳腺癌的淋巴結(jié)轉(zhuǎn)移引起了人們的廣泛關(guān)注。許多報道揭示了轉(zhuǎn)移性乳腺癌細胞的特征,但仍然需要更多乳腺癌相關(guān)淋巴結(jié)微環(huán)境成分和空間轉(zhuǎn)錄組圖譜的研究。Xu等[26]他們的研究中,通過整合scRNA-seq和空間轉(zhuǎn)錄組學,研究了6例手術(shù)切除淋巴結(jié)樣本的轉(zhuǎn)錄譜和1例陽性淋巴結(jié)的組織。結(jié)果發(fā)現(xiàn)高表達CD68和CD163的破骨細胞樣巨細胞。CD68和CD163是腫瘤相關(guān)巨噬細胞的生物標志物。通過空間轉(zhuǎn)錄組學方法,他們發(fā)現(xiàn)破骨細胞樣巨細胞分散在轉(zhuǎn)移性乳腺癌細胞中。在乳腺癌細胞浸潤的淋巴結(jié)微環(huán)境中,腫瘤相關(guān)巨噬細胞富含包括NF-κB信號通路和NOD樣受體信號通路在內(nèi)的前腫瘤通路。進一步的子聚類展示了巨噬細胞從活躍的趨化因子產(chǎn)生狀態(tài)發(fā)展為活躍的淋巴細胞激活狀態(tài)的潛在分化軌跡。該研究首次將scRNA-seq和空間轉(zhuǎn)錄組學整合到腋窩淋巴結(jié)腫瘤微環(huán)境中,為深入研究乳腺癌淋巴結(jié)轉(zhuǎn)移提供了系統(tǒng)方法。

盡管轉(zhuǎn)移仍然是大多數(shù)癌癥相關(guān)死亡的原因,但對遠端組織轉(zhuǎn)移播種的機制知之甚少。Davis等[27]建立了一種可靠的方法,使用單細胞RNA測序和患者衍生的乳腺癌異種移植模型,用于在過程中識別罕見轉(zhuǎn)移細胞的整體轉(zhuǎn)錄組變化。他們發(fā)現(xiàn)原發(fā)性腫瘤和微轉(zhuǎn)移都顯示出轉(zhuǎn)錄異質(zhì)性,但微轉(zhuǎn)移具有獨特的轉(zhuǎn)錄組程序,該轉(zhuǎn)錄組程序在患者衍生的異種移植模型中保守,可以高度預(yù)測患者的不良生存率。通路分析顯示線粒體氧化磷酸化是微轉(zhuǎn)移中上調(diào)的主要通路,而原發(fā)性腫瘤細胞中糖酵解酶水平較高,并且通過流式細胞術(shù)和代謝組學分析證實了這一點。通過抑制氧化磷酸化可以顯著減弱肺部的轉(zhuǎn)移播種,證明了氧化磷酸化在轉(zhuǎn)移中的重要性,并突出了其作為預(yù)防乳腺癌患者轉(zhuǎn)移擴散的治療靶點的潛力。

4 單細胞測序技術(shù)應(yīng)用于乳腺癌耐藥研究

耐藥性是乳腺癌治療的主要障礙。為了解決這個問題,Prieto-Vila等[28]對腔面型乳腺癌細胞系MCF7及其耐藥細胞進行了廣泛的單細胞基因表達譜分析。在耐藥細胞中觀察到上皮-間質(zhì)轉(zhuǎn)化和干性相關(guān)基因的上調(diào)和細胞周期相關(guān)基因的下調(diào)。這些基因主要受淋巴樣增強因子1調(diào)節(jié)。親代細胞中的少數(shù)細胞表現(xiàn)出與耐藥細胞相似的基因表達譜,表明未經(jīng)處理的親代細胞已經(jīng)含有耐藥的細胞亞群。

雌激素受體陽性乳腺癌的主要治療方法是內(nèi)分泌治療,臨床上大約90%的轉(zhuǎn)移性乳腺癌患者最終對內(nèi)分泌治療產(chǎn)生耐藥,至少33%的早期乳腺癌患者將產(chǎn)生內(nèi)分泌耐藥[29-31]。細胞周期蛋白依賴性激酶(Cyclin-dependent kinase,CDK)抑制劑聯(lián)合內(nèi)分泌治療能夠改善晚期轉(zhuǎn)移性雌激素受體陽性乳腺癌患者的預(yù)后,但其在早期患者中的價值尚不清楚。Griffiths等[32]運用scRNA-seq對單獨使用內(nèi)分泌治療或者聯(lián)用CDK抑制劑的內(nèi)分泌治療活體組織進行測序。研究者發(fā)現(xiàn)單獨用藥和聯(lián)合用藥均出現(xiàn)了共同的耐藥表型。聯(lián)合用藥治療的耐藥腫瘤顯示出雌激素信號的加速喪失,同時通過生長因子受體對JNK信號進行集中上調(diào)。相比之下,在單一或聯(lián)合治療期間維持雌激素信號傳導(dǎo)的癌細胞通過ERBB4信號傳導(dǎo),顯示CDK4/6激活和ERK上調(diào)的增強。這些結(jié)果表明,早期ER陽性乳腺癌的聯(lián)合治療通過從雌激素轉(zhuǎn)變?yōu)樘娲L信號介導(dǎo)的增殖而導(dǎo)致耐藥性的出現(xiàn)。

人表皮生長因子受體2(Human epidermal growth factor receptor2,HER2)陽性乳腺癌患者受益于HER2靶向治療,其中最常用的是曲妥珠單抗。然而,獲得性耐藥通常在接受靶向治療的一年內(nèi)發(fā)生,這種耐藥機制尚不清楚。Guoyu等[33]在兩種HER2陽性乳腺癌細胞系SKBR-3和BT474中通過耐藥篩選,構(gòu)建了曲妥珠單抗耐藥細胞。通過外顯子組測序技術(shù),檢測抗性誘導(dǎo)過程中不同時間點的細胞,用來研究基因組隨時間的變化。其中單細胞靶向測序被用于識別與耐藥性相關(guān)的基因突變。結(jié)果發(fā)現(xiàn)了拷貝數(shù)變異區(qū)域的快速增加,以及單核苷酸變異的逐漸積累。單細胞靶向測序發(fā)現(xiàn)在SKBR-3和BT474細胞中多個基因位點的突變導(dǎo)致腫瘤細胞對曲妥珠單抗敏感性降低。

5 單細胞測序技術(shù)應(yīng)用于乳腺癌異質(zhì)性與克隆演化研究

乳腺癌的異質(zhì)性使得不同患者之間(腫瘤間異質(zhì)性)、甚至在每個個體腫瘤內(nèi)(腫瘤內(nèi)異質(zhì)性)差異很大。腫瘤異質(zhì)性發(fā)生在形態(tài)學、基因組、轉(zhuǎn)錄組和蛋白質(zhì)組水平,給診斷和治療帶來挑戰(zhàn)[34]。

近年來,scRNA-seq在探索腫瘤細胞群體的異質(zhì)性以及腫瘤發(fā)生和轉(zhuǎn)移相關(guān)的稀有細胞類型方面,比傳統(tǒng)測序方法具有更大的優(yōu)勢[35]。scRNA-seq技術(shù)將具有不同分子亞型的乳腺癌細胞群體聚集在一起,來鑒定可能與不良預(yù)后和耐藥性相關(guān)的不同群體。scRNA-seq技術(shù)還用于鑒定潛在免疫治療靶點的不同免疫細胞亞群。此外,除了探索異質(zhì)性,scRNA-seq在乳腺癌研究中還具有多種應(yīng)用,包括分析細胞間信號通路、調(diào)節(jié)單細胞狀態(tài)和免疫細胞分布等。由于scRNA-seq能夠定義具有潛在治療靶點細胞亞群的能力,因此有助于乳腺癌的個體化治療。

異質(zhì)性和復(fù)雜的細胞結(jié)構(gòu)對乳腺癌的進展以及對治療的反應(yīng)具有重要的影響。然而,破譯腫瘤亞型及其空間組織仍然具有挑戰(zhàn)性。Liu等[36]將scRNA-seq與基于微陣列的空間轉(zhuǎn)錄組學相結(jié)合,來識別細胞群及其在乳腺癌組織中的空間分布。腫瘤細胞聚集成不同的亞群,這些細胞簇具有不同的起源、特征和功能。他們發(fā)現(xiàn)這些亞群分布在不同的組織區(qū)域,觀察到基質(zhì)細胞類型的差異富集,并且通過對大型乳腺癌RNA-seq隊列的去卷積推斷了這些腫瘤亞群的豐度,揭示了與患者生存和治療反應(yīng)的差異關(guān)聯(lián)。這項研究為乳腺癌的細胞結(jié)構(gòu)和潛在的治療策略提供了新的見解。

確定腫瘤異質(zhì)性及其對藥物反應(yīng)的影響能更好地對患者進行分層,幫助制定患者個性化治療方案。為了評估腫瘤細胞異質(zhì)性對藥物反應(yīng)的影響,Gambardella等[37]對來自32例乳腺癌細胞系的35 276個單個細胞進行轉(zhuǎn)錄分析,得到32個單細胞圖譜。研究者設(shè)計了一種單細胞轉(zhuǎn)錄譜映射到乳腺癌圖譜的算法,通過這種算法可以從腫瘤活檢的基因表達譜中確定細胞系組成,對患者來源的細胞系進行分層,從而預(yù)測藥物反應(yīng)。最后,他們將細胞系體外藥物篩選的結(jié)果與單細胞數(shù)據(jù)聯(lián)系起來進行比對,進而預(yù)測單細胞譜與藥物敏感性的關(guān)系。結(jié)果發(fā)現(xiàn)非遺傳的轉(zhuǎn)錄異質(zhì)性使具有不同藥物敏感性的細胞在相同的細胞系中也能共存。這對確定細胞系組成和藥物反應(yīng)方面的腫瘤異質(zhì)性提供了一個框架。

在克隆演化方面,拷貝數(shù)改變(Copy number alterations,CNA)在塑造乳腺癌基因組中發(fā)揮重要作用,并證明在臨床上可用于治療和預(yù)后。Balsan等[38]使用單細胞測序技術(shù),全面繪制出乳腺癌腫瘤隊列中拷貝數(shù)改變異質(zhì)性的特點。分析揭示:腫瘤塊內(nèi)非腫瘤細胞(即基質(zhì))的遺傳異質(zhì)性;拷貝數(shù)異質(zhì)性對乳腺癌基因組的影響程度以及亞克隆事件的基因組位置和劑量的重要性;染色體擴增的遺傳異質(zhì)性普遍存在;以及拷貝數(shù)異質(zhì)性與臨床和生物學參數(shù)(如多倍體和雌激素受體陰性狀態(tài))的關(guān)聯(lián)。他們的數(shù)據(jù)突出了單細胞基因組學在剖析多種形式的CNA腫瘤內(nèi)遺傳異質(zhì)性方面的特點,CNA異質(zhì)性對乳腺癌基因組的影響程度,以及CNA異質(zhì)性對腫瘤耐藥和復(fù)發(fā)的影響。

Kacavas等[6]對來自六個原發(fā)性三陰性乳腺癌患者超過1 500個細胞進行了scRNA-seq。結(jié)果發(fā)現(xiàn)每個腫瘤內(nèi)基因表達程序的細胞間異質(zhì)性是可變的,并且在很大程度上與推斷的基因組拷貝數(shù)變化的克隆性相關(guān)?；虮磉_譜的聚類確定了多個腫瘤共有的不同惡性細胞亞群,包括與治療抗性和轉(zhuǎn)移的多個特征相關(guān)的單個亞群,其功能特征是激活鞘糖脂代謝和相關(guān)的先天免疫途徑。定義該亞群的新特征預(yù)測了大型隊列中三陰性乳腺癌患者的長期結(jié)果。這項分析揭示了三陰性乳腺癌的功能異質(zhì)性及其與基因組進化的關(guān)聯(lián),并揭示了導(dǎo)致這種疾病預(yù)后不良的未預(yù)料到的生物學原理。

6 小結(jié)與展望

近年來,對單個細胞水平的分析正在改變?nèi)藗儗膊〉恼J識,腫瘤研究者開始逐漸認識到腫瘤組織和細胞群體的復(fù)雜性和異質(zhì)性,傳統(tǒng)的混合組織測序分析不能很好地表征生物學的復(fù)雜性。然而,在單細胞研究廣泛開展的大背景下,單細胞測序技術(shù)應(yīng)用于腫瘤研究也存在一些局限性。首先,單細胞測序技術(shù)的成本很高,需要一定的科研經(jīng)費支持;另外單細胞分析的常規(guī)方法還不完善,技術(shù)方面受限于單細胞的分離和捕獲、細胞膜的完整性及活性,從而使下游分析中的細胞失去了其在機體內(nèi)環(huán)境下的真實狀態(tài);其次,單細胞測序基因組學在單個實驗中測量成千上萬單細胞基因組信息,所產(chǎn)生的龐大數(shù)據(jù)的分析處理及對生物學信息的解析是一個巨大的挑戰(zhàn)。因此,與生物信息學和計算機生物學的跨學科合作將極大提高單細胞測序?qū)嶒烆I(lǐng)域的進展,對腫瘤學的研究具有重要的意義。