[文獻(xiàn)解讀]-全轉(zhuǎn)錄組RNA測(cè)序?qū)毙运柘蛋籽〉木C合診斷

百邁客一直秉承“生物科技長(zhǎng)信，服務(wù)社會(huì)，造福人民”的企業(yè)使命，致力于讓生物科技更快，更好的提高人類生活質(zhì)量。通過(guò)整合高通量測(cè)序技術(shù)、生物信息技術(shù)與云計(jì)算、大數(shù)據(jù)等新興IT技術(shù)，為用戶開(kāi)啟生物科技服務(wù)2.0新時(shí)代。各個(gè)測(cè)序的經(jīng)驗(yàn)也都是十分豐富，今天給大家?guī)?lái)一篇全轉(zhuǎn)錄組對(duì)白血病治療輔助的研究報(bào)道，了解應(yīng)用的新方向。詳解見(jiàn)下文：

英文題目：Comprehensive diagnostics of acute myeloid leukemia by whole transcriptome RNA sequencing

中文題目：全轉(zhuǎn)錄組RNA測(cè)序?qū)毙运柘蛋籽〉木C合診斷

期刊：Leukemia

影響因子：11.528

發(fā)表日期：2020.3.3

摘要

急性髓細(xì)胞白血病(AML)是由基因變異引起的，基因變異也控制著患者的預(yù)后，并指導(dǎo)風(fēng)險(xiǎn)適應(yīng)和靶向治療。急性髓系白血病的遺傳變異在結(jié)構(gòu)上是多樣的，目前通過(guò)不同的診斷分析來(lái)檢測(cè)。本研究旨在建立全轉(zhuǎn)錄組核糖核酸測(cè)序作為AML診斷的單一、全面和靈活的平臺(tái)。我們開(kāi)發(fā)了HAMLET(人類急性髓細(xì)胞白血病快速轉(zhuǎn)錄組學(xué))作為生物信息學(xué)管道，用于同時(shí)檢測(cè)fusion genes、small variants、tandem duplications和gene expression，所有信息都集中在一個(gè)被注釋的、簡(jiǎn)易的輸出文件中。對(duì)100例急性髓系白血病患者進(jìn)行了全轉(zhuǎn)錄組RNA測(cè)序，并通過(guò)reference assays和targeted resequencing驗(yàn)證了HAMLET結(jié)果。數(shù)據(jù)顯示，HAMLET z確地檢測(cè)到所有融合基因和過(guò)表達(dá)EVI1與3q26突變無(wú)關(guān)。此外，AML中經(jīng)常突變的13個(gè)基因的small variants有99.2%的敏感性和100%的特異性，以及FLT3和KMT2A的tandem duplications以100%的靈敏度和97.1%的特異性通過(guò)了新算法檢測(cè)。總之，HAMLET有潛力成為在單一技術(shù)平臺(tái)上提供與AML分類、風(fēng)險(xiǎn)評(píng)估和靶向治療相關(guān)的z確全面的診斷信息的平臺(tái)。

背景介紹

急性髓系白血病(AML)是由功能互補(bǔ)的基因突變引起的，這些突變導(dǎo)致髓系前體細(xì)胞不斷增殖和成熟被抑制。顯性 AML 平均攜帶13個(gè)非同義突變。造血干細(xì)胞中AML起始突變的獲得可能比AML診斷早幾十年。2016版的WHO classification of hematologic malignancies區(qū)分臨床和預(yù)后重要性的九種AML亞型。在50-55%的AML病例中進(jìn)行對(duì)遺傳變異的風(fēng)險(xiǎn)評(píng)估的z確分類，并可能指導(dǎo)靶向治療。在80-85%的AML病例對(duì)重復(fù)突變基因的擴(kuò)展組的基因組進(jìn)行分類。

然而，急性髓細(xì)胞白血病的z確分類需要不同的技術(shù)來(lái)捕獲結(jié)構(gòu)多樣的基因突變，如較大的插入和缺失(indels)、融合基因和除小變異外的結(jié)構(gòu)性染色體變異。此外，結(jié)構(gòu)正常基因的表達(dá)水平可對(duì)預(yù)后產(chǎn)生決定性影響。因此，在當(dāng)前的技術(shù)平臺(tái)下，急性髓細(xì)胞白血病診斷及風(fēng)險(xiǎn)評(píng)估仍然繁瑣、貴，且不完整。為了解決這一需求，我們實(shí)施了全轉(zhuǎn)錄組RNA測(cè)序(mRNAseq)，作為獲取與當(dāng)前和未來(lái)分類和風(fēng)險(xiǎn)分類相關(guān)的所有遺傳信息的單一平臺(tái)。在這種診斷范式中，一個(gè)經(jīng)認(rèn)可的mRNAseq協(xié)議獲取診斷樣本的全面數(shù)據(jù)。對(duì)于單次運(yùn)行數(shù)據(jù)分析，我們開(kāi)發(fā)了一個(gè)指定的集成生物信息管道(人類急性髓細(xì)胞白血病快速轉(zhuǎn)錄組學(xué)(HALLET)；Fig. 1)。HAMLET具有靈活的適應(yīng)性和個(gè)性化，可以查詢一組選定基因的序列變異和表達(dá)水平，以及易位情況。

實(shí)驗(yàn)材料

從萊頓大學(xué)醫(yī)學(xué)中心(LUMC)的血液學(xué)生物庫(kù)中選擇了100份經(jīng)機(jī)構(gòu)審查委員會(huì)批準(zhǔn)獲得的冷凍保存的急性髓系白血病樣本。所有樣本都進(jìn)行了基因分型NPM1突變和FLT3內(nèi)部串聯(lián)復(fù)制(FLT3-ITD)由LUMC分子診斷實(shí)驗(yàn)室認(rèn)證。

實(shí)驗(yàn)方法

全轉(zhuǎn)錄組測(cè)序，HAMLET輸出的驗(yàn)證

結(jié)果

經(jīng)認(rèn)證的mRNAseq在100個(gè)AML樣本上產(chǎn)生了4885萬(wàn)個(gè)reads(9700-17000萬(wàn)個(gè)reads)，平均插入大小為149-177個(gè)堿基。質(zhì)量控制后的平均閱讀長(zhǎng)度為123–124個(gè)堿基。
融合基因檢測(cè)

HAMLET在24例病例中z確檢測(cè)到中期細(xì)胞遺傳學(xué)預(yù)測(cè)的所有復(fù)發(fā)融合基因(Fig. 2a)。另外3例除t(9;11)(p21;q23)外的11q23易位病例，HAMLET正確識(shí)別了KMT2A融合相關(guān)基因MLLT1,MLLT4，和MLLT6(圖1，2)。此外，HAMLET在8例病例中正確識(shí)別了融合轉(zhuǎn)錄本，這些轉(zhuǎn)錄本不能通過(guò)中期核型分析預(yù)測(cè)，但與急性髓系白血病的預(yù)后相關(guān)(圖2b)。FUS-ERG 檢測(cè)到融合基因，其對(duì)應(yīng)的t(16；21)(p11；q22)僅在其復(fù)雜核型中被鑒定。攜帶ETV6-LYN融合基因的AML(2-020)的核型為add(8)(q24)和del(12)(p13)，但缺乏特征ins(12;8)(p13;q11q21)。三個(gè)AML (2-003，2-031，2-032)發(fā)現(xiàn)NUP98-NSD1融合轉(zhuǎn)錄本。由于兩個(gè)參與基因的端粒位置NUP98-NSD1源于t(5；11)(q35；p15.5)，因此目前未被AML分類，盡管其與不良預(yù)后相關(guān)。在三個(gè)病例(2-053、3-021、3-024)中發(fā)現(xiàn)PIM3-SCO2發(fā)生了基因融合。此前已在兒童急性髓系白血病中發(fā)現(xiàn)，可起源于慢性中性粒細(xì)胞白血病的0.6 Mb號(hào)染色體倒位。

HAMLET提供了融合轉(zhuǎn)錄物的詳細(xì)信息，包括差異剪接(KMT2A-MLLT3、RUNX1-RUNX1T1、PML-RARA、CBFB-MYH11和NUP98-NSD1)和來(lái)自平衡易位的相互轉(zhuǎn)錄物，例如，RARA-PML, ERG-FUS, NUP214-DEK, 和 NSD1-NUP98?？傊琀AMLET檢測(cè)到54種不同的融合轉(zhuǎn)錄本，這些轉(zhuǎn)錄物不能通過(guò)中期核型分析預(yù)測(cè)?？偟膩?lái)說(shuō)，28例急性髓系白血病攜帶兩個(gè)或多個(gè)融合，在一個(gè)復(fù)雜的AML中*多有9個(gè)不同的融合基因。

small variants檢測(cè)

HAMLET的VARSCAN組件在13個(gè)AML基因中鑒定出246個(gè)小變異體。HAMLET對(duì)NPM1和DNMT3A的FLT3-ITD和熱點(diǎn)區(qū)域進(jìn)行了敏感性評(píng)估，對(duì)其余的418個(gè)被HAMLET稱為野生型的樣本進(jìn)行了靶向重測(cè)序反應(yīng)。VARSCAN只識(shí)別了七個(gè)相當(dāng)短的FLT3-34人中的ITDFLT3-ITD攜帶AML。在這種情況下，HAMLET都正確地識(shí)別了一個(gè)ASXL1(3-030)的30bp C末端插入和CEBPA (3-003)中一個(gè)N端易位。HAMLET檢測(cè)到的小變異體的總體類型和分布與在COSMIC中報(bào)告的分布非常一致。

HAMLET檢測(cè)的變異在COSMIC中列出并位于定義的遺傳熱點(diǎn)內(nèi)，很容易被解釋為A類變異(n = 114;46.3%)。COSMIC中報(bào)告了60種變異(24.4%)，但沒(méi)有位于熱點(diǎn)。這些B類變異大部分在已知突變發(fā)生在整個(gè)編碼序列中的基因中被檢測(cè)到(CEBPA,TET2，RUNX1,WT1和TP53).其余變異并未報(bào)告(C類變異體:n= 52;21.1%)或未在COSMIC中列出但被描述為SNP(D類變異體，n= 20;8.1%).D類變異可能代表真正的低等位基因頻率的遺傳多態(tài)性。C類變異體可以根據(jù)其預(yù)測(cè)的致病作用來(lái)解釋。按這種規(guī)定進(jìn)行解釋，90%以上的A類變異，50%的B類和C類變異體和小于1%的D類變異體是髓系惡性腫瘤的復(fù)發(fā)突變。

在TCGA AML隊(duì)列中，基因組DNA中存在的小變體已被證明在相應(yīng)的RNA序列分析中是不存在的。因此，我們驗(yàn)證了小變異的HAMLET結(jié)果RUNX1,TET2,TP53,DNMT3A,KIT1，和IDH1通過(guò)靶向NGS對(duì)從56例急性髓系白血病分離的基因組DNA。選擇這些病例是因?yàn)門P53，RUNX1和TET2的遺傳特征與基因突變?cè)黾拥目赡苄杂嘘P(guān)，在這56例AML病例中，HAMLET檢測(cè)到67個(gè)變異，包括18個(gè)易位，它們是無(wú)義介導(dǎo)的RNA衰變的潛在靶點(diǎn)。在這67個(gè)變體中，有一個(gè)TET2三例急性髓系白血病病例中的變異體(c.4317dup)被目標(biāo)NGS遺漏。這個(gè)A插入發(fā)生在6個(gè)其他A核苷酸的延伸中，顯然很難在Ion GeneStudio S5系統(tǒng)上進(jìn)行排序。原始數(shù)據(jù)中存在這種變異的Reads，等位基因頻率為2-3%。還有兩個(gè)變異被目標(biāo)NGS發(fā)現(xiàn)，但被HAMLET錯(cuò)過(guò)了。這些變異包括RUNX1是已被靶向NGS檢測(cè)到的錯(cuò)義變異體，具有10%的低等位基因頻率和IDH1?？傊?，我們的數(shù)據(jù)顯示，HAMLETz確地調(diào)用了TET2,TP53,RUNX1,KIT,IDH1和DNMT3A包括由轉(zhuǎn)錄物編碼的所有變異，這些轉(zhuǎn)錄物是無(wú)義介導(dǎo)的衰變的潛在目標(biāo)。我們還確定了HAMLET-RNAseq檢測(cè)所需的測(cè)序深度，并對(duì)100個(gè)急性髓系白血病進(jìn)行了20萬(wàn)、30萬(wàn)、40萬(wàn)和4750萬(wàn)個(gè)讀對(duì)的流水線操作。在HAMLET所稱的246個(gè)變異中，有53個(gè)變異是無(wú)義突變或移碼突變。HAMLET檢測(cè)這些變異的z確度從100%、96%、85%下降到70%，而其余變異體的z確度從100%、98%、96%下降到81%。這些數(shù)據(jù)表明，約測(cè)5000萬(wàn)個(gè)reads的高測(cè)序深度對(duì)于確保檢測(cè)所有變異類型是必不可少的。

等位基因的CEBPA突變?yōu)楹诵驼５腁ML帶來(lái)了良好的預(yù)后并且將一個(gè)等位基因上編碼序列的第一個(gè)357 bp內(nèi)的N-末端移碼突變與c.834-1074之間的C-末端框內(nèi)插入或缺失相結(jié)合，從而破壞另一個(gè)等位基因上的DNA結(jié)合。通過(guò)Illumina平臺(tái)上通常可獲得的reads長(zhǎng)度來(lái)證明它們?cè)谙喾慈旧w上的位置。然而，雙等位基因CEBPA突變與特定的基因表達(dá)譜有關(guān)?；趍RNAseq的基因表達(dá)譜顯示了特征性的雙等位基因CEBPA攜帶兩個(gè)樣本的所有四個(gè)急性髓系白血病樣本的突變相關(guān)特征CEBPA突變(2-009、2-039、2-045和3-003)，包括其C端突變被HAMLET遺漏的病例3-003(圖3)。11例急性髓系白血病中，有一例CEBPA變異，一個(gè)病例(2-025)的特征性N-末端移碼突變與雙聚在一起CEBPA變異。因此，將變異檢測(cè)與基因表達(dá)譜結(jié)合分析，可以單獨(dú)解決變異檢測(cè)的不確定性。由于CEBPA 雙突變體的基因標(biāo)記不是定量的，所以尚未在HAMLET中實(shí)施臨床應(yīng)用。

檢測(cè)中的串聯(lián)重復(fù)FLT3和KMT2A

外顯子中常見(jiàn)串聯(lián)重復(fù)的可靠檢測(cè)FLT3在急性髓系白血病診斷中尤為重要，因?yàn)镕LT3-ITD是酪氨酸激酶抑制劑治療的適應(yīng)癥。除了對(duì)FLT3-ITD，其位置、大小和突變與等位基因比率對(duì)急性髓系白血病的預(yù)后有重要影響。為了改善對(duì)不充分的檢測(cè)FLT3-ITD通過(guò)VARSCAN，我們開(kāi)發(fā)了ReSCU作為一種新的算法，該算法基于在相互基因位置發(fā)生的部分對(duì)齊、SC reads的調(diào)用。

ReSCU確定了36例有FLT3-ITD，包括由VARSCAN檢測(cè)到的所有7例small ITD病例和經(jīng)認(rèn)證診斷確定的所有34例病例(圖4a)。NGS的有針對(duì)性的重新排序進(jìn)一步證實(shí)了FLT3-除了兩種情況之外，所有情況下，ITD都處于軟限幅方法確定的位置(表SXIV).在兩個(gè)由ReSCU調(diào)用但不是由標(biāo)準(zhǔn)診斷調(diào)用的不一致案例中，SC reads占其各自位置總覆蓋率的小于1%。仔細(xì)審查原始數(shù)據(jù)顯示FLT3-在這兩種情況下，ITD信號(hào)都低于檢測(cè)閾值，表明HAMLET的靈敏度更高(圖S6).突變型與正常型的等位基因比FLT3片段與ReSCU調(diào)用的SC reads百分比有很好的相關(guān)性(圖4b).

KMT2A-PTD外顯子2–13，這對(duì)于將AML 分類到染色質(zhì)剪接組，導(dǎo)致蛋白質(zhì)延長(zhǎng)，對(duì)AML的臨床結(jié)果產(chǎn)生不利影響。ReSCU識(shí)別的KMT2A-PTD外顯子2-8和2-10分別出現(xiàn)在4和2病例中，外顯子2、2-3和3-6的重復(fù)出現(xiàn)在單個(gè)病例中(圖4c)。全部九個(gè)KMT2A-PTD通過(guò) RT-PCR和測(cè)序進(jìn)行驗(yàn)證。

基因過(guò)表達(dá)的檢測(cè)

結(jié)構(gòu)正常基因的過(guò)度表達(dá)會(huì)影響急性髓系白血病患者的預(yù)后BAALC, ERG, MN1, DNMT3B, SPARC, and EVI1。mRNAseq測(cè)量過(guò)表達(dá)的可靠性被評(píng)估為預(yù)示AML預(yù)后非常差的EVI1轉(zhuǎn)錄本(圖5a)。過(guò)表達(dá)EVI1可由inv(3)(q21q26)或t(3;3)(q21;q26)，但是存在替代的分子機(jī)制。mRNAseq檢測(cè)到的EVIl外顯子1的表達(dá)與定量RT-PCR檢測(cè)結(jié)果有很好的相關(guān)性(圖5b).兩例(1- 005,3 -029)攜帶inv(3)(q21q26)和1例(1-008)攜帶t(3;8)(q26;q24)均可見(jiàn)高EVI1表達(dá)。在另外4例EVI1高表達(dá)的病例中，1例有der(3)t(1;3)(q3?1;q2?5)，沒(méi)有可識(shí)別的MECOM1位點(diǎn)參與(病例3-022)。這3例EVIl過(guò)表達(dá)而無(wú)3q26異常的病例均攜帶KMT2A融合轉(zhuǎn)錄本(case 1-001和3-008:KMT2A-MLLT3；case2-047:KMT2A-MLLT4)。這種遺傳變異的組合表明AML的預(yù)后極差。

Fig. 5 Detection of EVI1 overexpression. a Schematic representation of two mRNA transcripts from the MECOM locus on chromosome 3. One transcript contains exon 1–2 of MDS1 fused to exon 2–15 of EVI1 (MDS1-EVI1 transcript), whereas the other transcript contains exon 1–15 of EVI1 (EVI1 transcript). b Comparison of EVI1 expression by RNAseq and quantitative RT-PCR. Xaxis: Expression of the first exon of EVI1 normalized for expression of the PBGD housekeeping gene by quantitative RT-PCR (log2 EVI1/PBGD). Y-axis: Sum of base coverage of the first exon of EVI1 per kb transcript and one million mapped reads (log2 BPKM) by RNAseq

克隆進(jìn)化

HAMLET揭示了克隆進(jìn)化，在診斷和復(fù)發(fā)時(shí)采集的所有四對(duì)樣本中獲得和丟失了遺傳變異。盡管基因型明顯不同，但I(xiàn)DH2突變作為起始事件的持續(xù)存在確立了新生AML和后續(xù)治療相關(guān)AML的共同起源(case 3-021、3-038)。因此，一個(gè)獨(dú)立的、與治療相關(guān)的AML的診斷被糾正為原始AML的復(fù)發(fā)，盡管兩個(gè)樣本將被分類為不同的AML亞型。

急性髓系白血病的分類和預(yù)測(cè)

根據(jù)HAMLET和中期核型分析結(jié)果，99例病例可根據(jù)世衛(wèi)組織2016年和基因組分類進(jìn)行重新分類(圖6)。HAMLET信息促進(jìn)了6個(gè)案例的風(fēng)險(xiǎn)重新評(píng)估，值得注意的是通過(guò)NUP98-NSD1融合而無(wú)分類明確病變的AML，或通過(guò)EVIl過(guò)表達(dá)而無(wú)inv(3)(q21q26)來(lái)分配風(fēng)險(xiǎn)。

 

原文鏈接：

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

百邁客轉(zhuǎn)錄組優(yōu)勢(shì)：

超高性價(jià)比：一次測(cè)序，4組數(shù)據(jù)，多篇文章

樣本全：組織、細(xì)胞、粘膜、毛囊均可建庫(kù)測(cè)序

經(jīng)驗(yàn)豐富：常規(guī)建庫(kù)起始量低；低起始量的珍惜樣本也可建庫(kù)成功

強(qiáng)大生信分析：可提供基本分析和高階分析，滿足個(gè)性化調(diào)整

全套服務(wù)：組織樣本取樣、建庫(kù)、測(cè)序、數(shù)據(jù)分析全套服務(wù)

云分析平臺(tái)：百邁客云平臺(tái)已經(jīng)研發(fā)快速高效的云平臺(tái)，可高效呈現(xiàn)高質(zhì)量文章圖表和對(duì)大數(shù)據(jù)進(jìn)行挖掘；

百邁客的測(cè)序服務(wù)結(jié)合強(qiáng)大的云分析平臺(tái)，可以根據(jù)您的研究目的提供個(gè)性化分析，針對(duì)性地挖掘出豐富多彩的結(jié)果圖譜。