陳柳 華清泉

[摘要] 轉(zhuǎn)錄組（RNA-Seq）是新一代測(cè)序技術(shù)，利用高通量測(cè)序技術(shù)對(duì)組織或細(xì)胞中所有RNA反轉(zhuǎn)錄而成的cDNA文庫(kù)進(jìn)行測(cè)序，從整體水平研究基因的功能及其結(jié)構(gòu)，揭示特定生物學(xué)過(guò)程以及疾病發(fā)生過(guò)程中的分子機(jī)制，而癌癥療效評(píng)估的關(guān)鍵在于腫瘤的復(fù)發(fā)率和對(duì)放化療的敏感性。新一代測(cè)序系統(tǒng)有望識(shí)別導(dǎo)致腫瘤的形成和對(duì)治療產(chǎn)生抵抗的相關(guān)靶基因，對(duì)進(jìn)一步指導(dǎo)臨床治療和診斷具有重要的理論及實(shí)踐意義。

[關(guān)鍵詞] 轉(zhuǎn)錄組學(xué)；RNA-Seq；腫瘤；測(cè)序技術(shù)

[中圖分類號(hào)] R733 [文獻(xiàn)標(biāo)識(shí)碼] A [文章編號(hào)] 1673-7210（2016）02（c）-0058-04

Application and progress of RNA-Seq in tumor

CHEN Liu HUA Qingquan

Department of Otolaryngology & Head and Neck Surgery， Renmin Hospital of Wuhan University， Hubei Province， Wuhan 430060， China

[Abstract] The transcriptome is a new generation sequencing technology， which uses high-throughput sequencing technology to sequence cDNA library formed by the reverse transcription of all RNA from tissue or cells， studies gene function and structure from the overall level， reveals the specific molecular mechanism in the process of biological processes and disease. The key point evaluating the efficacy of cancer is the tumor recurrence rate and sensitivity to radiation and chemotherapy. The new generation sequencing system is expected to identify the related target genes which lead to the formation of tumor and resistance to treatment. Further more， it has important theoretical and practical significance for guiding clinical treatment and diagnosis.

[Key words] Transcriptomics； RNA-Seq； Tumor； Sequencing technology

轉(zhuǎn)錄組是指特定組織或細(xì)胞在某一發(fā)育階段或功能狀態(tài)下轉(zhuǎn)錄出來(lái)的所有RNA的集合，包括蛋白質(zhì)編碼的信使RNA和非編碼RNA（rRNA、tRNA和其他ncRNAs）。1995年Velculescu等首次提出了關(guān)于轉(zhuǎn)錄組的概念，并在酵母細(xì)胞中獲得第1個(gè)轉(zhuǎn)錄組[1]。腫瘤的發(fā)生發(fā)展大多是多種基因改變積累造成的，包括致癌基因的激活和抑癌基因的失活。RNA-Seq技術(shù)可以從轉(zhuǎn)錄水平檢測(cè)腫瘤中的改變基因。通過(guò)相關(guān)軟件分析，鑒定相關(guān)改變基因在腫瘤發(fā)生發(fā)展中所起的作用。然后經(jīng)過(guò)刷選，進(jìn)一步挑選出相關(guān)靶基因進(jìn)行實(shí)驗(yàn)研究，有望為腫瘤的基因治療提供新的靶點(diǎn)。

1 轉(zhuǎn)錄組學(xué)的概念及進(jìn)展

1.1 轉(zhuǎn)錄組學(xué)的概念

高通量測(cè)序技術(shù)也被稱為下一代測(cè)序（NGS）技術(shù)，推進(jìn)了基因組學(xué)的研究進(jìn)展。該項(xiàng)新技術(shù)除了研究靜態(tài)基因組，近年來(lái)也開(kāi)始應(yīng)用于動(dòng)態(tài)轉(zhuǎn)錄分析，由此衍生的技術(shù)稱為RNA序列分析（RNA-seq）[2]。在20世紀(jì)90年代末和21世紀(jì)初，隨著DNA基因芯片技術(shù)的快速發(fā)展，開(kāi)創(chuàng)了癌癥近十年的高通量轉(zhuǎn)錄組研究。早期研究的顯著發(fā)現(xiàn)主要是關(guān)于黑色素瘤[3]和乳腺癌[4]。van't Veer等[5]曾使用基因芯片對(duì)乳腺癌淋巴結(jié)陽(yáng)性患者進(jìn)行相關(guān)研究，確定了一組70個(gè)基因的堿基排列順序，可以很好地預(yù)測(cè)乳腺癌患者的預(yù)后[6]。近年來(lái)，基于利用轉(zhuǎn)錄組所得到的成果，轉(zhuǎn)錄組逐漸成為研究癌癥生物學(xué)信息的一個(gè)主導(dǎo)實(shí)驗(yàn)技術(shù)。

1.2 轉(zhuǎn)錄組學(xué)的相關(guān)進(jìn)展

復(fù)合性基因融合現(xiàn)象常見(jiàn)于與血液學(xué)關(guān)聯(lián)的惡性腫瘤和罕見(jiàn)的骨骼及軟組織腫瘤中。近年來(lái)研究表明，常見(jiàn)的實(shí)體腫瘤中也有復(fù)發(fā)性基因融合現(xiàn)象。融合基因是最普遍的癌基因，是由染色體重組造成的。Maher等[7]綜合運(yùn)用轉(zhuǎn)錄組技術(shù)，在慢性粒細(xì)胞性白血病細(xì)胞系中發(fā)現(xiàn)融合基因BCR-ABL1以及在前列腺的癌細(xì)胞系和組織中發(fā)現(xiàn)融合基因TMPRSS2-ERG。同時(shí)Huang等[8]研究提示RNA-Seq技術(shù)可以有效識(shí)別移位的染色體，并發(fā)現(xiàn)融合基因PARP14-TFE3。另一方面，可變剪接是導(dǎo)致人類蛋白質(zhì)多樣性的主要來(lái)源。單個(gè)基因經(jīng)過(guò)可變剪接，形成多種不同的mRNAs。據(jù)目前情況估計(jì)，基于全基因組技術(shù)的應(yīng)用，研究表明超過(guò)90%的人類基因都進(jìn)行過(guò)可變剪接[9-10]。在許多癌癥中發(fā)現(xiàn)可變剪接。許多與腫瘤生物學(xué)功能相關(guān)的關(guān)鍵蛋白質(zhì)經(jīng)過(guò)了與癌癥相關(guān)的可變剪接，這些腫瘤生物學(xué)功能包括細(xì)胞凋亡、細(xì)胞周期調(diào)控、入侵和轉(zhuǎn)移[11-13]。近年來(lái)，全基因組方法的顯著發(fā)展有望發(fā)現(xiàn)腫瘤中改變的可變剪接基因，有助于發(fā)現(xiàn)腫瘤細(xì)胞中發(fā)生變化的具體信號(hào)通路。

2 RNA測(cè)序（RNA-Seq）

2.1 RNA-Seq的技術(shù)平臺(tái)

NGS已經(jīng)徹底改變了既往對(duì)研究細(xì)胞通路的系統(tǒng)分析方法。癌癥是一種復(fù)雜的疾病，涉及多種異常的基因表達(dá)調(diào)控。識(shí)別這種驅(qū)動(dòng)腫瘤生物學(xué)發(fā)展的表達(dá)模式為我們理解癌癥提供了很好的分子基礎(chǔ)[14]。NGS的兩個(gè)主要平臺(tái)是全外顯子組測(cè)序（WES）和RNA-Seq。WES主要是用于發(fā)現(xiàn)DNA的變異，RNA-Seq主要用于衡量基因的表達(dá)。兩種技術(shù)都可以用于檢測(cè)基因變異，尤其是單核苷酸變異（SNVs）[15]。鑒于持續(xù)降低成本和改進(jìn)數(shù)據(jù)分析方法的優(yōu)勢(shì)，NGS已經(jīng)開(kāi)始直接影響生物醫(yī)學(xué)的發(fā)展[16]，其測(cè)序步驟依賴于NGS技術(shù)。目前使用最廣泛并且廣受歡迎的應(yīng)用于RNA-seq的NGS平臺(tái)有3個(gè)，包括454年焦磷酸測(cè)序系統(tǒng)、AB固體系統(tǒng)和Illummina基因組分析儀。具體來(lái)說(shuō)，Illumina/Solexa Illumina公司（美國(guó)圣地亞哥Illumina公司Inc. CA）平臺(tái)是基于可逆染料終結(jié)合測(cè)序原理，而454年和固體系統(tǒng)采用創(chuàng)新型的乳液聚合酶鏈反應(yīng)機(jī)制。

RNA-seq也被稱為整個(gè)轉(zhuǎn)錄組鳥(niǎo)槍測(cè)序（wts）。RNA-seq包含三個(gè)步驟：測(cè)序文庫(kù)的建立，在特定NGS平臺(tái)中測(cè)序和生物信息學(xué)分析。RNA-Seq是一種基于NGS技術(shù)的新的轉(zhuǎn)錄分析方法。采用NGS技術(shù)對(duì)從感興趣的RNA樣本序列反轉(zhuǎn)錄得到的cDNA進(jìn)行測(cè)序[17-19]。簡(jiǎn)而言之，就是對(duì)剪切變體、轉(zhuǎn)錄后RNA編輯、內(nèi)含子表達(dá)水平和特定等位基因的表達(dá)情況進(jìn)行定性和定量分析。

2.2 RNA-Seq原理

在對(duì)基因組測(cè)序和分析研究的同時(shí)，復(fù)雜的轉(zhuǎn)錄組研究也廣泛發(fā)展起來(lái)。利用高通量測(cè)序技術(shù)平臺(tái)分析轉(zhuǎn)錄組的結(jié)構(gòu)和表達(dá)水平情況，更能挖掘未知轉(zhuǎn)錄本和稀有轉(zhuǎn)錄本信息，精確地識(shí)別RNA 的選擇性剪切以及編碼序列的單核苷酸多態(tài)性（SPN），進(jìn)一步解析復(fù)雜的轉(zhuǎn)錄組信息[20]。最初的轉(zhuǎn)錄組研究主要以基因芯片微陣列技術(shù)為基礎(chǔ)，而基因芯片技術(shù)的檢測(cè)范圍取決于芯片上的探針信息，所以只能檢測(cè)已知序列的特征，缺少發(fā)現(xiàn)新基因的能力。而高通量測(cè)序技術(shù)可以很好地彌補(bǔ)基因芯片技術(shù)在這方面的不足。因此，現(xiàn)階段轉(zhuǎn)錄組的研究是借助于NGS[21]來(lái)完成的。通過(guò)構(gòu)建cDNA文庫(kù)并對(duì)cDNA 進(jìn)行高通量測(cè)序，分析測(cè)序結(jié)果進(jìn)而解析轉(zhuǎn)錄組學(xué)中的復(fù)雜變化。目前，以基因測(cè)序技術(shù)為核心的新技術(shù)平臺(tái)支撐體系已經(jīng)相對(duì)成熟和完善，例如：Illumina公司的Solexa測(cè)序技術(shù)、羅氏公司的454測(cè)序技術(shù)、ABI公司的SOLID測(cè)序技術(shù)以及美國(guó)螺旋生物科學(xué)公司的新型納米孔測(cè)序技術(shù)等[13]。

隨著NGS技術(shù)的不斷更新和提高而日益成熟完善，RNA測(cè)序技術(shù)平臺(tái)在測(cè)序通量、測(cè)序長(zhǎng)度、錯(cuò)配率、堿基配對(duì)讀取能力等測(cè)序性能方面技術(shù)的提高均有利于轉(zhuǎn)錄組的研究。RNA-Seq技術(shù)的不斷更新和創(chuàng)新，為人類逐步全面了解真核生物和原核生物的轉(zhuǎn)錄組信息提供了新的定性和定量的生物信息學(xué)方法。

2.3 RNA-Seq技術(shù)優(yōu)勢(shì)

相比基因芯片和標(biāo)記的堿基序列的測(cè)序方法，RNA-seq具有一定的優(yōu)越性。與基因芯片相比，RNA-Seq不僅可以在更高的分辨率下用來(lái)測(cè)量基因的表達(dá)水平情況，還可以揭示未知的轉(zhuǎn)錄組和拼接亞型，并為選擇性拼接亞型提供定量分析[22]。相對(duì)于傳統(tǒng)的芯片雜交平臺(tái)，轉(zhuǎn)錄組測(cè)序無(wú)需預(yù)先針對(duì)已知序列設(shè)計(jì)探針，即可對(duì)任意物種的整體轉(zhuǎn)錄活動(dòng)進(jìn)行檢測(cè)，提供更精確的數(shù)字化信號(hào)、更高的檢測(cè)通量以及更廣泛的檢測(cè)范圍，是目前深入研究復(fù)雜性轉(zhuǎn)錄組的強(qiáng)大工具。鑒于這些優(yōu)勢(shì)，RNA-Seq很快就被應(yīng)用到癌癥相關(guān)研究的多個(gè)方面。Brooks等[23]研究發(fā)現(xiàn)，相比于qRT-PCR，NGS能夠更加全面以及更加準(zhǔn)確地定量和定性評(píng)價(jià)細(xì)胞或組織中的mRNA。NGS可以對(duì)細(xì)胞或組織中的RNA分子的所有亞型進(jìn)行量化和綜合評(píng)價(jià)。同時(shí)，RNA-seq技術(shù)可以檢測(cè)到低水平表達(dá)、未識(shí)別的轉(zhuǎn)錄子和新拼接亞型基因[24]。另一方面，與基因芯片相比，RNA-Seq有一個(gè)更廣的檢測(cè)動(dòng)態(tài)范圍，可以研究編碼和非編碼RNA，更易于基因網(wǎng)絡(luò)的構(gòu)建、發(fā)現(xiàn)選擇性剪切轉(zhuǎn)錄物，檢測(cè)到等位基因的具體表達(dá)方式，也可以用來(lái)提取基因型信息。在RNA-Seq單細(xì)胞研究[25]中，一個(gè)主要優(yōu)勢(shì)是可以分析整個(gè)轉(zhuǎn)錄序列，而不是有限制數(shù)量的基因，并且，研究DNA和RNA的方法并不是相互排斥的，他們可以相互結(jié)合使用，產(chǎn)生更多的生物學(xué)上重要的信息[26]。

3 RNA-Seq在腫瘤中的應(yīng)用

癌癥的生物起源、發(fā)展、轉(zhuǎn)移與轉(zhuǎn)錄組中的許多變化息息相關(guān)?；虻倪x擇性剪接是一個(gè)主要的轉(zhuǎn)錄后調(diào)控機(jī)制，參與許多惡性腫瘤的發(fā)生發(fā)展。NGS應(yīng)用于RNA層面上的研究，為研究轉(zhuǎn)錄組提供了一種新技術(shù)。同時(shí)，實(shí)體腫瘤與生物學(xué)特征及臨床行為息息相關(guān)，在治療方面越來(lái)越凸顯出個(gè)體化治療的重要性，而RNA-Seq是一種比較合適的技術(shù)來(lái)評(píng)估個(gè)性化的治療效果[27]。

3.1 RNA測(cè)序與食道鱗狀細(xì)胞癌

食管癌是一種常見(jiàn)的惡性腫瘤。陳友芳等[28]對(duì)6對(duì)食管鱗癌（ESCC）組織及同一患者的正常食管黏膜組織進(jìn)行RNA-Seq分析，發(fā)現(xiàn)142個(gè)差異表達(dá)基因。Ma等[29]對(duì)3例確診為ESCC患者的癌組織和癌旁組織進(jìn)行RNA-Seq分析。在癌組織中，PTK6均為低表達(dá)，并運(yùn)用免疫組化、RT-PCR、流式細(xì)胞儀等技術(shù)進(jìn)行驗(yàn)證，發(fā)現(xiàn)PTK6可能是通過(guò)影響Akt/GSK3/β-catenin信號(hào)通路，從而促進(jìn)腫瘤的形成和轉(zhuǎn)移。Jiang等[30]通過(guò)對(duì)一個(gè)確診為ESCC患者的癌組織和癌旁組織進(jìn)行RNA-Seq分析，確定新型融合基因HLA-E和HLA-B。該融合基因有望為ESCC的診斷和治療提供新的靶點(diǎn)，使我們更加了解ESCC的形成機(jī)制。

3.2 RNA測(cè)序與肺癌

肺癌的5年生存率不到15%。Beane等[31]利用RNA測(cè)序技術(shù)分析了抽煙對(duì)肺癌的影響，揭示了抽煙與肺癌的關(guān)系以及肺癌的發(fā)生機(jī)制。Riccardo等[32]在K-rasLA1/p53R172HΔg小鼠非小細(xì)胞肺癌模型中，運(yùn)用RNA測(cè)序技術(shù)發(fā)現(xiàn)GM-CSF和ADORA3的高表達(dá)均與腫瘤的轉(zhuǎn)移相關(guān)。ADORA3參與肺癌細(xì)胞系中p53誘導(dǎo)凋亡過(guò)程。在p53突變的腫瘤中，ADORA3有望成為潛在免疫治療的靶點(diǎn)。lncRNAs的破壞在非小細(xì)胞肺癌的發(fā)生發(fā)展中起著關(guān)鍵作用，Ricciuti等[33]揭示lncRNAs是早期非小細(xì)胞肺癌患者潛在生物標(biāo)志物，可以預(yù)測(cè)腫瘤患者對(duì)化療和靶向治療的敏感性。

3.3 RNA測(cè)序與前列腺癌

前列腺癌是男性泌尿系統(tǒng)中最常見(jiàn)的惡性腫瘤之一，同時(shí)也是男性中因癌癥死亡的第二大原因。Palanisamy等[34]對(duì)前列腺癌（PCa）組織標(biāo)本進(jìn)行RNA-Seq分析來(lái)探索新的嵌合體融合轉(zhuǎn)錄片段，最后分析發(fā)現(xiàn)了7種腫瘤特異性融合基因——其中有兩個(gè)是ETS基因家族中的ETl和ERG基因，4個(gè)是非ETS基因如CDKNlA（p21），CD9和IKBκB（IKK-beta）。該研究首次把非ETS基因納入到PCa融合基因中。Eswaran等[35]分析了15例前列腺癌樣本的RNA-Seq結(jié)果，驗(yàn)證了RAF信號(hào)途徑在癌癥發(fā)病過(guò)程中起到了非常重要的作用，同時(shí)也表明了RAF信號(hào)途徑中的融合基因有潛力成為抗腫瘤治療與抗腫瘤藥物篩選的靶點(diǎn)。Xu等[36]利用RNA測(cè)序數(shù)據(jù)對(duì)前列腺癌的體細(xì)胞突變進(jìn)行篩查，得到一組可用于前列腺癌的診斷和治療的候選基因。

3.4 RNA-Seq面臨的挑戰(zhàn)及發(fā)展前景

目前，RNA-seq技術(shù)仍處于起步階段，然而它的優(yōu)越性明顯超過(guò)傳統(tǒng)的芯片雜交測(cè)序方法。在過(guò)去的幾年中，RNA-seq在研究人類癌癥的轉(zhuǎn)錄注釋、基因轉(zhuǎn)錄結(jié)構(gòu)的確定、不同轉(zhuǎn)錄生物過(guò)程的定量分析中做出了不可磨滅的貢獻(xiàn)。隨著RNA-seq技術(shù)不斷發(fā)展和成本持續(xù)下降，在未來(lái)幾年中，將RNA-seq技術(shù)應(yīng)用于人類腫瘤的研究中，相信很有可能會(huì)產(chǎn)生更加令人振奮的發(fā)現(xiàn)。

此外，RNA-Seq實(shí)驗(yàn)還可以揭示出那些基因芯片不能檢測(cè)出的新型的轉(zhuǎn)錄子、非編碼RNA和其他小分子RNA 。我們認(rèn)識(shí)到，拼接變異以及其他基因組變化是腫瘤重要的驅(qū)動(dòng)因素[34]。RNA-Seq可以有效捕獲內(nèi)含子水平的亞型和非編碼基因。在尋找與癌癥預(yù)后、診斷及治療靶點(diǎn)的生物標(biāo)志上，RNA-Seq是一種很有潛力的實(shí)驗(yàn)技術(shù)，為研究腫瘤的發(fā)生、發(fā)展提供了一個(gè)前所未有的機(jī)遇。

惡性腫瘤嚴(yán)重威脅人類健康，是世界各國(guó)醫(yī)學(xué)科學(xué)領(lǐng)域中的重大科研課題。目前每年大約有200萬(wàn)例新發(fā)病例，130萬(wàn)人死于惡性腫瘤，是目前中國(guó)人的第三大死因。早治療、早診斷可以提高患者生活質(zhì)量。我們相信，基于網(wǎng)絡(luò)的生物學(xué)方法可以為解釋腫瘤的發(fā)病機(jī)制和發(fā)現(xiàn)腫瘤的新藥物靶點(diǎn)奠定有力基礎(chǔ)；為腫瘤患者的預(yù)診斷及治療提供理論及實(shí)驗(yàn)依據(jù)，有望進(jìn)一步提高腫瘤患者的生活質(zhì)量。RNA-Seq可以有效捕獲內(nèi)含子水平的亞型和非編碼基因[19]。在尋找與癌癥預(yù)后、診斷及治療靶點(diǎn)的生物標(biāo)志上，RNA-Seq是一種很有潛力的實(shí)驗(yàn)技術(shù)，為研究腫瘤的發(fā)生、發(fā)展提供了一個(gè)前所未有的機(jī)遇。

[參考文獻(xiàn)]

[1] Bawa P，Zackaria S，Verma M，et al. Integrative analysis of normal long intergenic non-coding RNAs in prostate cancer [J]. PLoS One，2015，10（5）：1-14.

[2] Qian X，Ba Y，Zhuang Q，et al. RNA-Seq technology and its application in fish transcriptomics [J]. OMICS，2014，18（2）：98-110.

[3] DeRisi J，Penland L，Brown PO，et al. Use of a cDNA microarray to analyse gene expression patterns in human cancer [J]. Nature Genetics，1996，14（4）：457-460.

[4] Kononen J，Bubendorf L，Kallioniemi A，et al. Tissue microarrays for high-throughtput molecular profiling of tumor specimens [J]. Nat Med，1998，4（7）：844-847.

[5] van't Veer LJ，Dai H，van de Vijver MJ，et al. Gene expression profiling predicts clinical outcome of breast cancer [J]. Nature，2002，415（6871）：530-536.

[6] Roukos DH. Next-generation sequencing and epigenometechnologies：potential medical applications [J]. Expert Rev Med Devices，2010，7（6）：723-726.

[7] Maher CA，Kumar-Sinha C，Cao X，et al. Transcriptome sequencing to detect gene fusions in cancer [J]. Nature，2009，458（7234）：97-101.

[8] Huang W，Goldfischer M，Babyeva S，et al. Identification of a novel PARP14-TFE3 gene fusion from 10-year-old FFPE tissue by RNA-seq [J]. Genes Chromosomes Cancer，2015. [Epub ahead of print]

[9] Wang ET，Sandberg R，Luo S，et al. Alternative isoform regulation in human tissue transcriptomes [J]. Nature，2008，456（7221）：470-476.

[10] Pan Q，Shai O，Lee LJ，et al. Deep surveying of alternative splicing complexity in the human transcriptome by high-throughput sequencing [J]. Nat Genet，2008，40（12）：1413-1415.

[11] Venables JP. Aberrant and alternative splicing in cancer [J]. Cancer Res，2004，64（21）：7647-7654.

[12] Liu S，Cheng C. Alternative RNA splicing and cancer [J]. Wiley Interdiscip Rev RNA，2013，4（5）：547-566.

[13] David CJ，Manley JL. Alternative pre-mRNA splicing regulation in cancer：pathways and programs unhinged [J]. Genes Dev，2010，24（21），2343-2364.

[14] Danan-Gotthold M，Golan-Gerstl R，Eisenberg E，et al. Identification of recurrent regulated alternative splicing events across human solid tumors [J]. Nucleic Acids Res，2015，43（10）：5130-5144.

[15] O'Brien TD，Jia P，Xia J，et al. Inconsistency and features of single nucleotide variants detected in whole exome sequencing versus transcriptome sequencing：a case study in lung cancer [J]. Methods，2015，83：118-127.

[16] Ding L，Wendl MC，Koboldt DC，et al. Analysis of next generation genomic data in cancer：accomplishments and challenges [J]. Hum Mol Genet，2010，19（R2）：188-196.

[17] Morozova O，Hirst M，Marra MA. Applications of new sequencing technologies for transcriptome analysis [J]. Annu Rev Genomics Hum Genet，2009，10：135-151.

[18] Wilhelm BT，Marguerat S，Watt S，et al. Dynamic repertoire of a eukaryotic transcriptome surveyed at single-nucleotide resolution [J]. Nature，2008，453（7199）：1239-1243.

[19] Feng H，Qin Z，Zhang X. Opportunities and methods for studying alternative splicing in cancer with RNA-Seq [J]. Cancer Lett，2013，340（2）：179-191.

[20] Ding X，Zhu L，Ji T，et al. Long intergenic non-coding RNAs（LincRNAs）identified by RNA-Seq in breast cancer [J]. PLoS One，2014，9（8）：1-10.

[21] Iyer MK，Niknafs YS，Malik R. The landscape of long noncoding RNAs in the human transcriptome [J]. Nat Genet，2015，47（3）：199-208.

[22] Peng Z，Cheng Y，Tan BC. Comprehensive analysis of RNA-Seq data reveals extensive RNA editing in a human transcriptome [J]. Nat Biotechnol，2012，30（3）：253-260.

[23] Brooks MJ，Rajasimha HK，Roger JE，et al. Next-generation sequencing facilitates quantitative analysis of wild-type and Nrl-/- retinal transcriptomes [J]. Mol Vis，2011， 17：3034-3054.

[24] Trapnell C，Williams BA，Pertea G，et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation [J]. Nat Biotechnol，2010，28（5）：511-515.

[25] Hitzemann R，Bottomly D，Darakjian P，et al. Genes，behavior and next-generation RNA sequencing [J].Genes Brain Behav，2013，12（1）：1-12.

[26] Lee MC，Lopez-Diaz FJ，Khan SY，et al. Single-cell analyses of transcriptional heterogeneity during drug tolerance transition in cancer cells by RNA sequencing [J]. Proc Natl Acad Sci，2014，111（44）：4726-4735.

[27] Meiβner T，F(xiàn)isch KM，Gioia L，et al. OncoRep：an n-of-1 reporting tool to support genome-guided treatment for breast cancer patients using RNA-sequencing [J]. BMC Med Genomics，2015，8（24）：1-8.

[28] 陳友芳，曹訓(xùn)，黃志良，等.RNA測(cè)序技術(shù)在篩選食管鱗癌差異表達(dá)基因中的應(yīng)用[J].廣東醫(yī)學(xué)，2014，35（6）：844-847.

[29] Ma S，Bao JY，Kwan PS，et al. Identification of PTK6，via RNA sequencing analysis，as a suppressor of esophageal squamous cell carcinoma [J]. Gastroenterology，2012，143（3）：675-686.

[30] Jiang YZ，Li QH，Zhao JQ，et al. Identification of a novel fusion gene（HLA-E and HLA-B）by RNA-seq analysis in esophageal squamous cell carcinoma [J]. Asian Pac J Cancer Prev，2014，15（5）：2309-2312.

[31] Beane J，Vick J，Schembri F，et al. Characterizing the impact of smoking and lung cancer on the airway transcriptome using RNA-Seq [J]. Cancer Prev Res（Phila），2011，4（6）：803-817.

[32] Riccardo F，Arigoni M，Buson G，et al. Characterization of a genetic mouse model of lung cancer：a promise to identify non-small cell lung cancer therapeutic targets and biomarkers [J]. BMC Genomics，2014，15（Suppl 3）：S1.

[33] Ricciuti B，Mencaroni C，Paglialunga L，et al. Long noncoding RNAs：new insights into non-small cell lung cancer biology，diagnosis and therapy [J]. Med Oncol，2016， 33（2）：1-18.

[34] Palanisamy N，Ateeq B. Rearrangements of the RAF kinase pathway in prostate cancer，gastric cancer and melanoma [J]. Nat Med，2010，16（7）：793-798.

[35] Eswaran J，Horvath A，Godbole S，et al. RNA sequencing of cancer reveals novel splicing alterations [J]. Sci Rep，2013，3：1689.

[36] Xu YM，Liu HZ. Related biomarkers in the diagnosis of prostate cancer [J]. Zhonghua Nan Ke Xue，2015，21（10）：937-940.

（收稿日期：2015-10-02 本文編輯：張瑜杰）