腫瘤組織是由正常細(xì)胞和腫瘤細(xì)胞組成的復(fù)雜組織，在以往基因組或轉(zhuǎn)錄組的研究中，很難對(duì)其中各種細(xì)胞間關(guān)系和組織微環(huán)境進(jìn)行系統(tǒng)全面的解釋。那么，如何將其與單細(xì)胞技術(shù)結(jié)合，根據(jù)單細(xì)胞水平的基因表達(dá)情況對(duì)每一個(gè)細(xì)胞進(jìn)行精確分類呢？

在單細(xì)胞測(cè)序中，對(duì)腫瘤組織細(xì)胞進(jìn)行分類的方法通?；?strong>無(wú)監(jiān)督聚類分析，根據(jù)已知Marker基因的表達(dá)對(duì)細(xì)胞類型進(jìn)行判定，以及基于拷貝數(shù)變異(copy number variation, CNV)從單細(xì)胞的表達(dá)數(shù)據(jù)推斷細(xì)胞類型。前者基于聚類分析的方法可能受到多種因素影響，包括聚類算法、用于注釋的Marker基因選擇等，準(zhǔn)確判定需要結(jié)合大量過(guò)往生物學(xué)知識(shí)；而后者基于cnv的方法則需要進(jìn)行大量密集型計(jì)算，相對(duì)耗時(shí)。

因此，Mariano Alvarez 團(tuán)隊(duì)提出了一種名為基因組不穩(wěn)定性分析(genomic instability Analysis, GIA)的新方法，運(yùn)行速度相比inferCNV快至少5~6倍，能夠快速評(píng)估單個(gè)細(xì)胞內(nèi)的基因組變化情況，輔助腫瘤細(xì)胞的判定。

Genomic Instability原理

該軟件使用aREA算法，對(duì)單細(xì)胞基因表達(dá)譜中的連續(xù)編碼基因(位點(diǎn)塊)進(jìn)行富集分析，從而定量估計(jì)出基因表達(dá)與染色體位置之間的關(guān)聯(lián)性（inferCNV）。

genomicInstability軟件主要整合了三大功能函數(shù)：

①inferCNV：量化表達(dá)譜與位點(diǎn)塊之間的關(guān)聯(lián)函數(shù)；

②genomicInstabilityScore：估計(jì)每個(gè)細(xì)胞的基因組不穩(wěn)定性評(píng)分；

③giLikelihood：估計(jì)每個(gè)細(xì)胞基因組不穩(wěn)定的可能性。

基因組不穩(wěn)定性分析

2.1 不穩(wěn)定性打分和高斯建模

首先根據(jù)物種信息（人/小鼠），對(duì)表達(dá)數(shù)據(jù)的基因和蛋白ID進(jìn)行匹配，得到與數(shù)據(jù)庫(kù)基因位點(diǎn)塊對(duì)應(yīng)的表達(dá)數(shù)據(jù)，從而進(jìn)行關(guān)聯(lián)性分析。軟件用歸一化富集分?jǐn)?shù)（NES）表示每個(gè)細(xì)胞在每個(gè)位點(diǎn)塊的富集程度，該位點(diǎn)塊富集到的基因表達(dá)范圍越廣，NES離散度越大，該區(qū)域的基因組不穩(wěn)定性也越高。接著根據(jù)富集分?jǐn)?shù)計(jì)算基因組不穩(wěn)定性評(píng)分（GIS），并對(duì)整體細(xì)胞建立多個(gè)高斯模型。最后，根據(jù)GIS由高斯模型確定細(xì)胞基因組不穩(wěn)定概率（gi_likelihood_score）。

表1 不同細(xì)胞基因組不穩(wěn)定性評(píng)分示

從表1可以看出，control組較case組的細(xì)胞基因組不穩(wěn)定概率更低，可以較明顯的區(qū)分開(kāi)腫瘤和非腫瘤細(xì)胞。

2.2 模型可視化和CNV熱圖

分析過(guò)程中會(huì)默認(rèn)建立三個(gè)高斯模型，在不同高斯模型中，具有最低GIS均值的細(xì)胞被定義為基因組穩(wěn)定細(xì)胞，其他模型則代表基因組不穩(wěn)定細(xì)胞（即腫瘤細(xì)胞）。由圖1可以看出，灰色/藍(lán)色模型中更多的是GIS在-1 ~1.5 的非腫瘤細(xì)胞，而腫瘤細(xì)胞主要集中在右側(cè)GIS大于1.5的橙色模型中。

圖1 基因組不穩(wěn)定性評(píng)分分布以及高斯模型示例

根據(jù)基因組不穩(wěn)定性評(píng)分（GIS）和不穩(wěn)定概率，可以繪制染色體水平熱圖和featureplot，從而更直觀地看出基因組、各clusters細(xì)胞的不穩(wěn)定情況。圖2熱圖中，紅色表示位點(diǎn)塊區(qū)域存在擴(kuò)增，藍(lán)色則表示存在缺失或刪除，其中Unstable表示細(xì)胞gi_likelihood_score 大于0.8， Stable 小于0.2。

圖2 染色體水平的位點(diǎn)塊富集熱圖示例

圖3 featureplot圖中，紅色區(qū)域表示細(xì)胞有較高的不穩(wěn)定性，可以看出右側(cè)一小簇細(xì)胞有更高的不穩(wěn)定性。

圖3 基因組不穩(wěn)定性featureplot可視化示例

genomicInstability軟件能夠根據(jù)單細(xì)胞基因表達(dá)數(shù)據(jù)，快速預(yù)估每個(gè)細(xì)胞是腫瘤細(xì)胞或正常細(xì)胞的概率，為腫瘤細(xì)胞的判定提供參考。

參考文獻(xiàn)：

[1] Patel, A. P. et al. Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma. Science 344, 1396–1401 (2014).

[2] Puram, S. V. et al. Single-Cell Transcriptomic Analysis of Primary and Metastatic Tumor Ecosystems in Head and Neck Cancer. Cell 171, 1611–1624 (2017).

[3] Muller, S. & Diaz, A. Single-Cell mRNA Sequencing in Cancer Research: Integrating the Genomic Fingerprint. Front. Genet. 8, 73 (2017).

[4] Fan, J., Slowikowski, K. & Zhang, F. Single-cell transcriptomics in cancer: computational challenges and opportunities. Exp. Mol. Med. 52, 1452–1465 (2020).

[5] Fan, J. et al. Linking transcriptional and genetic tumor heterogeneity through allele Analysis of single-cell RNA-seq data. Genome Res. 28, 1217–1227 (2018).

[6] Serin Harmanci, A., Harmanci, A. O. & Zhou, X. CaSpER identifies and visualizes CNV events by integrative Analysis of single-cell or bulk RNA-sequencing data. Nat. Commun. 11, 89 (2020)

[7] Alvarez, M. J. et al. Functional characterization of somatic mutations in cancer using network-based inference of protein activity. Nat. Genet. 48, 838–47 (2016).

本文系歐易生物原創(chuàng)