空間轉錄組+基因組數據聯(lián)合分析揭示乳腺原位癌的異質(zhì)性

文章題目：Genomic profiling reveals heterogeneous?populations of ductal carcinoma in situ of?the breast

中文題目：基因組分析揭示乳腺原位癌的異質(zhì)性群體

期刊：Communications Biology

影響因子：6.6

DOI：10.1038/s42003-021-01959-9

導讀

對大部分患者來(lái)說(shuō)，乳腺導管原位癌（DCIS）不會(huì )發(fā)展為浸潤性導管癌，在當前的臨床標準下，這些患者經(jīng)常被過(guò)度治療。雖然有各種候選標記物可用，但尚未建立用于描述風(fēng)險類(lèi)別的相關(guān)標記物。在本研究中，作者分析了431例DCIS患者的臨床特征，并對21例患者進(jìn)行了全外顯子組測序分析，對72例患者進(jìn)行了靶向深度測序分析。確定了年齡<45歲、HER2擴增和GATA3突變可能影響復發(fā)。PIK3CA突變陰性和PgR陰性也被認為是危險因素?？臻g轉錄組分析進(jìn)一步顯示GATA3功能障礙上調上皮-間充質(zhì)轉化和血管生成，隨后PgR下調。這些結果揭示了DCIS中異質(zhì)細胞群的存在，并為DCIS分類(lèi)和優(yōu)化治療提供了預測性標記。

百邁客不僅提供單細胞系列產(chǎn)品，6大核心產(chǎn)品，單細胞轉錄組、空間轉錄組、單細胞免疫組庫、單細胞ATAC-seq、單細胞核轉錄組、單細胞全長(cháng)轉錄組，而且實(shí)現了單細胞轉錄組、空間轉錄組百邁客云端APP分析，深入大數據挖掘，一鍵快速分析，單細胞和空間轉錄組的聯(lián)合分析近似實(shí)現單細胞分辨率的空間異質(zhì)性，助您進(jìn)入時(shí)空組學(xué)研究新時(shí)代。接下來(lái)我們來(lái)看看作者如何實(shí)現基因組數據和空間轉錄組數據聯(lián)合研究乳腺原位癌的異質(zhì)性。

方法

外顯子組文庫構建和全外顯子組測序

從21例DCIS（原發(fā)性腫瘤）患者和4例IDC（復發(fā)性腫瘤）患者中選擇福爾馬林固定石蠟包埋（FFPE）組織樣本。同時(shí)對每位患者制作相應的冰凍切片，染色確定DCIS和正常上皮組織中的區域。

全外顯子組測序中基因組風(fēng)險因子的選擇

為了提高基于臨床評估區分低風(fēng)險和高風(fēng)險DCI標準的準確性，作者嘗試識別基因組因素。評估了21例DCIS患者，他們是隨機選擇的，具有匹配的固有亞型和相對較高的復發(fā)率。對匹配的正常乳腺組織進(jìn)行測序（平均深度×109.1），以區分體細胞突變的種系變異。獲得的變異圖如圖1a所示。

GATA3和PIK3CA突變是隊列中最常見(jiàn)的突變，分別出現在4例（19%）和5例（24%）患者中。在4例GATA3突變患者中，3例出現IDC復發(fā)。GATA3突變陽(yáng)性與GATA3突變陰性復發(fā)的優(yōu)勢比（OR）為5.5（95%CI=0.63–76.7）。為了進(jìn)一步研究GATA3突變在復發(fā)中的意義，對復發(fā)腫瘤患者的樣本進(jìn)行額外分析，以確定與配對原發(fā)性DCI相比的突變狀態(tài)。

結果

全外顯子靶向測序

獲得的突變圖如圖2所示?？偣灿?0名（56%）和36名（50%）患者分別攜帶GATA3和PIK3CA突變。與上述結果一致，GATA3突變與復發(fā)呈正相關(guān)（OR=7.8；95%CI=1.17~88.4），而PIK3CA突變往往與復發(fā)呈負相關(guān)（OR=0.45；95%可信區間=0.12~1.7）。在測試的180個(gè)基因中，沒(méi)有一個(gè)比這兩個(gè)因素更能作為預測標記。對于PgR蛋白表達，低PgR表達再次與復發(fā)相關(guān)（OR=25.6；95%可信區間=3.64–142.2）。結果表明在11例中的8例中發(fā)現了GATA3突變，盡管ER陽(yáng)性，但PgR較低。

DCIS患者的分子解剖分析

為了確定所識別的高危標記物的分子病因，調查了另外三名具有代表性的新鮮冰凍標本患者?；颊逜患有帶有GATA3突變（圖2，*1）和微侵襲的DCIS，病變被認為是真正的DCIS?；颊連患有PgR低表達的DCIS（圖2，*2）和微侵襲，病變被視為真正的DCIS?；颊逤患有PIK3CA突變的DCIS（圖2，*3），沒(méi)有任何微侵襲，病變被認為可能是假DCIS（低風(fēng)險DCIS）。

含有GATA3突變的DCIS病變的STseq

首先，利用空間轉錄組學(xué)技術(shù)——10X基因組學(xué)Visium平臺，通過(guò)STseq分析了患者A標本的空間基因表達。該患者在術(shù)前病理診斷中被診斷為DCIS，但在術(shù)后病理診斷中發(fā)現有微侵襲部位。因此，DCIS病變被確定為IDC的真正前兆。平均深度為×2391.4的面板測序分析顯示該樣本含有GATA3突變（外顯子4:c.865dupG:p。C288fs，變異等位基因頻率[VAF]=7%）。在較高的VAF（外顯子5:c.T1035A:p。N345K，VAF=26.1%），表明新出現的GATA3突變用PIK3CA突變覆蓋了DCIS的基本特征。

PgR低表達但無(wú)GATA3突變的DCIS病變的STseq

接下來(lái)，作者將患者B作為PgR低表達的例外患者進(jìn)行研究，盡管沒(méi)有GATA3突變（圖2，*2）。從平均深度×3030.6的靶向測序中，檢測到PIK3CA突變（外顯子21:c.A3140T:p。H1047L，VAF=40%），而即使在這個(gè)足夠的測序深度，也沒(méi)有檢測到GATA3突變。因此，從分子水平來(lái)看，該患者的預后良好。然而，術(shù)后病理檢查顯示微侵襲，表明DCIS病變是IDC的真正前兆。免疫組化分析顯示該患者ER陽(yáng)性，PgR陰性。

DCIS和PIK3CA突變

最后，將患者C的病例分析為可能的假DCIS病例，不太可能發(fā)展為IDC（低風(fēng)險DCIS）。該患者特征為PIK3CA突變，GATA3突變或下調缺失，無(wú)微侵襲，術(shù)后無(wú)復發(fā)。該患者的DCIS病變包含兩個(gè)PIK3CA突變（外顯子10:c.G1633A:p。E545K，外顯子10:c.A1634G:p。E545G，兩種突變的VAF=25%）。

對于Visium分析，收集了799133428個(gè)測序讀數，測序飽和度為86.2%。分析點(diǎn)的數量（直徑55μm）為2043個(gè)，其中包含7469個(gè)唯一分子標識符（UMI）讀數的中位數和每個(gè)點(diǎn)2928個(gè)檢測基因的中位數。如圖3A的中間板所示，癌細胞斑點(diǎn)被分為三組，通過(guò)非分層k-均值聚類(lèi)（K＝9），表明DCIS病變中存在異質(zhì)性。此外，可能代表癌細胞周?chē)h(huán)境的非癌點(diǎn)被分為四類(lèi)。在代表病理學(xué)鑒定的癌細胞的總共422個(gè)點(diǎn)中，通過(guò)分析Visium讀數在46個(gè)點(diǎn)中檢測到GATA3突變讀數（圖3b）。幸運的是，在該患者中，GATA3突變位于轉錄本的3′端，因此可以通過(guò)Visium讀數來(lái)表示。為了研究攜帶GATA3突變（spots）的細胞是否是克隆性的，檢測了46個(gè)攜帶GATA3突變的點(diǎn)的VISUM讀數上的GATA3突變位點(diǎn)。突變完全存在于同一位點(diǎn)（外顯子4:c.865dupG:p。C288fs），表明它們是單克隆來(lái)源。使用TCC R檢測了GATA3突變的斑點(diǎn)之間的差異表達基因（DEG）?？偣矙z測到1468個(gè)DEG（錯誤發(fā)現率[FDR]<0.05，圖3c，左面板）。如圖1d所示，在具有GATA3突變的斑點(diǎn)中觀(guān)察到PgR表達下調（P=0.02901；t檢驗；圖3c，右面板）。為了預測檢測到的DEG的功能結果，使用Metascape進(jìn)行了路徑分析(http://metascape.org/)28.發(fā)現了9條與11種癌癥特征相關(guān)的關(guān)鍵途徑。

與突變陰性點(diǎn)的結果相比，GATA3突變的癌點(diǎn)受影響（圖3d）。特別是，關(guān)鍵基因和關(guān)鍵途徑包括EMT和血管生成。重要的是，具有表達變化的關(guān)鍵基因是那些被鑒定為對異常GATA3功能作出反應的基因，如VIM和FN1（圖3e）。這些結果表明GATA3突變在DCIS進(jìn)展過(guò)程中出現，并伴有惡性特征。據報道，增加的VIM表達發(fā)生在管腔癌細胞的GATA3突變下游24；因此，在這些斑點(diǎn)中觀(guān)察到的GATA3突變可能代表惡性發(fā)展過(guò)程中先前的基因改變。相反，在沒(méi)有GATA3突變的斑點(diǎn)中，DEG主要富集于雌激素反應、緊密連接和mTORC1的信號通路，表明這些斑點(diǎn)中的細胞獲得了允許細胞轉化的最小變化。

Visium分析發(fā)現了異質(zhì)基因表達，這取決于癌區和非癌區細胞的位置。PgR表達下調，盡管觀(guān)察到高ER表達（圖4a，上圖），與免疫染色結果一致（圖4a，下圖）。針對DCIS細胞中基因表達的變化，手動(dòng)選擇了188個(gè)位于導管內區域的形態(tài)學(xué)斑點(diǎn)。這些斑點(diǎn)的無(wú)監督分層聚類(lèi)確定了三個(gè)明顯的聚類(lèi)（圖4b，左面板）。在空間上，不同的簇對應于不同的區域（圖4b，中間面板）。重要的是，簇1的48個(gè)斑點(diǎn)（紅色）與DCIS細胞的位置重疊，DCIS細胞根據其形態(tài)即將入侵（由圖4b右面板中的紅色箭頭指示）。同時(shí)，簇2（有色綠色）的101個(gè)斑點(diǎn)位于同一導管的中心，并被視為不侵入基質(zhì)的細胞（由圖4B的中間板中的綠色斑點(diǎn)指示）。簇1和簇2之間的差異表達分析確定了2747個(gè)DEG（FDR<0.05，圖4c，左）。有趣的是，GATA3表達在簇1中顯著(zhù)下調（圖4c，右圖）。盡管該基因未發(fā)生突變，但該患者的GATA3功能異?？赡馨l(fā)生在基因水平。根據基因富集和通路分析的結果，該患者始終觀(guān)察到以GATA3為中心的基因表達變化，如患者A所觀(guān)察到的（圖4d，e）。共有21條途徑和18個(gè)癌癥特征，包括EMT和血管生成途徑受到影響。如對患者A所觀(guān)察到的，在未侵入間質(zhì)的第2組中，DEG主要代表雌激素反應相關(guān)基因?？傊?，在這名患者身上，我們得出結論，GATA3通過(guò)其表達的變化發(fā)揮關(guān)鍵作用。即使沒(méi)有基因組突變，GATA3的異常表達也可能在某些患者中導致相同的結果。

通過(guò)全外顯子組測序，平均深度為×134.8。Visium分析揭示了一種單調的基因表達模式，與該癌的形態(tài)學(xué)單克隆結構一致（圖5a，左面板）。事實(shí)上，這些斑點(diǎn)通過(guò)非層次聚類(lèi)大致分為兩個(gè)簇（圖5a，中間和右側面板）。它們幾乎與形態(tài)學(xué)確定的癌細胞和非癌細胞（基質(zhì)細胞）完全重疊。因此，使用DEG分析比較了癌細胞中的347個(gè)斑點(diǎn)（藍色斑點(diǎn)）和基質(zhì)細胞中的151個(gè)斑點(diǎn)（橙色斑點(diǎn)），并確定了癌斑點(diǎn)中508個(gè)上調的DEG（FDR<0.05）。加權基因共表達網(wǎng)絡(luò )分析（WGCNA）顯示，與非癌點(diǎn)相比，癌點(diǎn)中只有兩個(gè)主要節點(diǎn)（模塊）發(fā)生改變（圖5b，上面板）。進(jìn)行了子網(wǎng)絡(luò )分析，以確定樞紐基因（圖5b，中間和下部面板）。hub基因的途徑富集分析顯示，參與雌激素反應和p53途徑的基因在模塊1中富集（圖5c，補充數據5），這與在患者A和B中觀(guān)察到的良性出現斑點(diǎn)的特征相同，表明該病變中沒(méi)有發(fā)生惡性轉化。

討論

在本研究中，將GATA3突變確定為DCI分類(lèi)的潛在標記。ER陽(yáng)性DCIS患者的PIK3CA突變陰性和PgR蛋白陰性也被認為是危險因素。傳統的臨床病理因素對于預測IDC發(fā)生的風(fēng)險也至關(guān)重要，包括年齡(≥45歲與<45歲）和HER2狀態(tài)（陰性與陽(yáng)性）。在缺乏HER2靶向治療的情況下，HER2擴增是眾所周知的IDC預后不良因素。然而，HER2擴增的頻率在DCIS和IDC33–37之間有所不同，HER2擴增是否是DCIS的風(fēng)險因素仍有爭議。HER2陽(yáng)性率/年齡的HR≥45歲與2歲消極/年齡≥45年為1.54（95%可信區間=0.31-7.65，P=0.595）。HER2陰性/年齡<45歲與HER2陰性/年齡的HR比較≥45為2.43（95%可信區間=0.78-7.52，P=0.125）。相反，HER2陽(yáng)性/年齡<45歲與HER2陰性/年齡的HR≥45為11.03（95%可信區間=3.54-34.4，P<0.0001）。盡管交互試驗的P值為0.121，但這些結果可能表明年齡與HER2擴增之間存在交互作用。還發(fā)現PIK3CA突變的DCIS細胞沒(méi)有任何導致惡性轉化的基因組或轉錄組改變。未來(lái)，DCIS的臨床研究(包括正在進(jìn)行的試驗12 14)可能會(huì )提供一種準確的算法，使用包括本研究提出的標記物在內的標志物來(lái)區分高風(fēng)險和低風(fēng)險DCIS，從而避免不必要的治療。

參考文獻

Nagasawa S, Kuze Y, Maeda I,?et al. Genomic profiling reveals heterogeneous populations of ductal carcinoma in situ of the breast. Commun Biol. 2021 Apr 1;4(1):438.

 

文獻下載

https://international.biocloud.net/zh/article/detail/33795819

?

推薦閱讀

文獻分享|【綜述】單細胞技術(shù)開(kāi)啟對肝臟病理學(xué)及發(fā)病機制的新認知

單細胞轉錄組助力構建不同胃腺癌亞型的轉錄組圖譜

百邁客：10x單細胞轉錄組與空間轉錄組聯(lián)合分析一睹為快

空間轉錄組-在生長(cháng)發(fā)育中的應用2

空間轉錄組-在癌癥中的應用1

單細胞轉錄組應用4——結直腸癌研究

單細胞&空間轉錄組應用3-動(dòng)植物研究“遍地開(kāi)花”

單細胞轉錄組應用2——食管鱗狀細胞癌研究

單細胞轉錄組應用1——衰老研究