亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于TCGA數(shù)據(jù)庫前列腺癌預(yù)后相關(guān)mRNA風(fēng)險模型的構(gòu)建

        2021-06-25 12:19:06張漢榮林聯(lián)拯于忠英李金雨
        關(guān)鍵詞:差異基因前列腺癌精度

        張漢榮,林聯(lián)拯,趙 力,于忠英,李金雨

        (中國人民解放軍聯(lián)勤保障部隊第909醫(yī)院(廈門大學(xué)附屬東南醫(yī)院)泌尿外科,福建 漳州363000)

        0 引言

        前列腺癌是一種常見而復(fù)雜的惡性腫瘤,根據(jù)最新的全球癌癥統(tǒng)計數(shù)據(jù),全球每年約有130萬例前列腺癌新發(fā)病例,死亡病例達36萬例,占所有腫瘤新發(fā)病例的15%,是男性第二常見的腫瘤類型和第五大癌癥死亡原因[1-2].目前,前列腺癌的治療方法主要有手術(shù)治療、放療、化療和激素治療[3-4].選擇何種治療手段與腫瘤的分期、分級及危險度等因素密切相關(guān)[5].前列腺癌的預(yù)后模型的構(gòu)建、或利用構(gòu)建的預(yù)測模型指導(dǎo)前列腺癌患者分層治療是前列腺癌的研究熱點之一.目前的預(yù)測模型大多數(shù)是以患者的年齡、腫瘤大小、TNM分期等臨床資料為影響因素,預(yù)測模型的效率和精度上仍存在不足[6-8].隨著高通量基因測序技術(shù)的發(fā)展和成熟,利用癌癥基因組或轉(zhuǎn)錄組數(shù)據(jù)構(gòu)建的癌癥預(yù)后模型顯示出良好的預(yù)測性能[9-10].本研究的目的是,利用公共數(shù)據(jù)庫TCGA中提供的前列腺癌RNA-Seq數(shù)據(jù),篩選潛在的前列腺癌相關(guān)基因,并利用這些候選基因構(gòu)建預(yù)測模型,為臨床前列腺癌的風(fēng)險評估、分層治療提供參考依據(jù).

        1. 資料與方法

        1.1 研究資料

        從the cancer genome atlas(TCGA)數(shù)據(jù)庫(http://www.tcga.org/)下載前列腺癌及對應(yīng)癌旁組織的mRNA表達數(shù)據(jù)與臨床資料.TCGA數(shù)據(jù)庫共收錄前列腺癌癌組織492例,癌旁組織52例.

        1.2 前列腺癌組與對應(yīng)癌旁組織差異基因的篩選與富集分析

        根據(jù)TCGA下載的數(shù)據(jù),利用R軟件的“edge”包篩選出癌組織與癌旁組織間差異表達的mRNA.篩選標準為:|Log2FC|≥1,Padj≤0.05.利用Volcano Plot包對差異基因數(shù)據(jù)繪制火山圖.利用clusterProfiler包對前列腺癌組織的上調(diào)基因分別進行GO富集和KEGG富集分析,以了解前列腺癌與癌旁組織差異基因的基因功能和信號通路的富集情況.

        1.3 差異基因?qū)η傲邢侔┗颊哳A(yù)后的影響

        以上調(diào)最明顯的10個基因為候選基因,分別以上述基因mRNA表達量中值為分界線,將492例前列腺癌患者分為高表達組和低表達組,并做生存分析.比較差異基因?qū)η傲邢侔╊A(yù)后的影響.生存分析采用Log-rank檢驗方法,P<0.05為差異有統(tǒng)計學(xué)意義.

        1.4 COX風(fēng)險回歸模型的構(gòu)建

        將上述生存分析有統(tǒng)計學(xué)意義的基因納入COX回歸分析,以各基因的表達量為自變量,患者的生存時間為因變量.COX回歸模型計算各患者的風(fēng)險得分(Risk score),最后根據(jù)患者的風(fēng)險得分將492例患者分為高危組和低危組,進一步進行生存分析,比較高危組和低危組的生存差異.ROC曲線用于驗證該風(fēng)險模型的預(yù)測精度.

        2. 結(jié)果

        2.1 前列腺癌對應(yīng)癌旁組織的差異基因及富集分析

        基因表達差異分析結(jié)果如圖1A所示,共篩選表達上調(diào)基因1978個,下調(diào)基因1644個.其中,上調(diào)最明顯的基因為:PCA3、AMACR、MTND4P12、RNY3P8、DLX1、OR51E2、PCAT14、GOLM1、HPN、GLYATL1.下調(diào)最明顯的基因為:ADAM33、MEG3、GSTM2、GABRE、CYP3A5、CORO6、GOLGA8A、DIO3OS、SLC26A10、HIF3A.

        對上調(diào)基因做GO功能富集分析,結(jié)果如圖1B所示,上調(diào)的基因主要與carboxylic acid biosynthetic process,organic acid biosynthetic process及sulfur compound metabolic process等代謝過程有關(guān).KEGG富集分析(圖1C)顯示,上調(diào)基因主要集中在Phagosome、Cell adhesion molecules和Hematopopoietic aureus infection等信號通路.

        圖1 差異基因火山圖及GO功能、KEGG通路富集分析結(jié)果

        2.2 差異基因?qū)η傲邢侔┗颊哳A(yù)后的影響

        上調(diào)最明顯的10個基因?qū)︻A(yù)后的影響分別如圖2所示,生存分析有顯著差異的基因有6個,分別為PCA3、MTND4P12、RNY3P8、OR51E2、PCAT14、GOLM1.且上述基因的生存分析中,均表現(xiàn)為低表達組的總生存率較高表達組高.因此,考慮上述6個基因可能是前列腺癌發(fā)生、發(fā)展的關(guān)鍵基因.

        圖2 10個差異基因?qū)?92例前列腺癌患者預(yù)后的影響

        2.3 COX風(fēng)險回歸模型的構(gòu)建

        以各基因的表達量為自變量,生存時間為因變量,得到風(fēng)險得分公式為:Risk score=0.5503×(PCA3)+0.6777×(MTND4P12)+0.5605×(RNY3P8)-0.365×(OR51E2)+0.5358×(PCAT14)-0.1972×(GOLM1).計算492例前列腺癌患者的風(fēng)險得分,并以風(fēng)險得分的中值將患者分為高風(fēng)險組與低分險組(圖3A),該模型的ROC曲線(圖3B)曲線下面積為0.801,故認為該模型的預(yù)測精度較好.高、低風(fēng)險組的生存分析結(jié)果如圖3C所示,結(jié)果提示,高風(fēng)險組的總生存時間明顯較低風(fēng)險組低,P<0.0001.從高、低風(fēng)險組的生存曲線可看出,多個基因表達水平構(gòu)建的預(yù)測模型,較單一基因表達水平的生存分析差異更大,提示多基因表達構(gòu)建的預(yù)測模型比單一基因構(gòu)建的模型預(yù)測精度更高.

        圖3 COX風(fēng)險回歸模型的構(gòu)建

        3. 討論

        與其他腫瘤相比,前列腺癌具有更明顯的腫瘤異質(zhì)性,不同個體間的差異大[11].良好的預(yù)后模型可以促進臨床咨詢和指導(dǎo)醫(yī)生制定治療和隨訪計劃.目前臨床上用于評估前列腺癌患者預(yù)后的模型主要有前列腺特異抗原(PAS)水平,Gleason評分腫瘤臨床分期等[12-14].然而,上述用于評估前列腺癌患者危險度的模型的精確度仍不能滿足準確的患者分類和分層的需求.因此,臨床上亟需一種可精確預(yù)測患者危險度的預(yù)測模型,以指導(dǎo)前列腺癌患者的指導(dǎo)治療,實現(xiàn)個體化治療.目前已有大量基因單個基因表達量構(gòu)建的前列腺癌預(yù)后模型,如Mee YK等[15]的研究提示,高表達FOXM1基因與前列腺癌預(yù)后差相關(guān).除FOXM1基因外,近年來,ARHGAP10[16],SNW1[17],NCAPH[18]等基因也有類似報道.但是基于單基因構(gòu)建的預(yù)測模型預(yù)測效能及模型的穩(wěn)定性相對較差.L Schmidt等[19]通過對接受根治性前列腺切除術(shù)的123名男性的前列腺癌組織樣本進行全基因組miRNA表達譜分析,篩選出6個與前列腺癌預(yù)后最相關(guān)的miRNA,并利用其中的4個miRNA的表達量比值構(gòu)建預(yù)測模型.結(jié)果顯示該預(yù)測模型對前列腺癌具有良好的預(yù)測作用.

        本研究通過對TCGA數(shù)據(jù)庫收錄的492例前列腺癌組織和52例癌旁組織的全基因組測序結(jié)果分析,篩選出在癌組織中高表達的基因.以上調(diào)最明顯的10個基因為候選基因,分別以各基因的表達水平中值為界,比較高表達組與低表達組間患者預(yù)后的差異.結(jié)果提示,10個候選基因,有6個(PCA3、MTND4P12、RNY3P8、OR51E2、PCAT14、GOLM1)基因的表達量與預(yù)后相關(guān),均表現(xiàn)為高表達提示預(yù)后不良.進一步對聯(lián)合上述6個基因做COX風(fēng)險回歸模型.以上述6個基因的表達量為自變量,患者生存資料為因變量.結(jié)果提示,聯(lián)合6個基因構(gòu)建的風(fēng)險預(yù)測模型具有良好的預(yù)測精度(AUC=0.801).根據(jù)模型可利用各患者上述6個基因的表達量計算該患者的風(fēng)險得分.本研究進一步將492例患者風(fēng)險得分的中值為分界線,比較高、低分險組生存差異.結(jié)果提示,高風(fēng)險組的生存時間明顯較低風(fēng)險組短(P<0.00001).聯(lián)合6個基因構(gòu)建的風(fēng)險模型做的生存分析的P值明顯小于單一基因生存分析的P值,提示,該風(fēng)險模型具有更好的預(yù)測精度.

        綜上所述,本研究利用生物信息學(xué)方法,成功構(gòu)建了一個以6基因(PCA3、MTND4P12、RNY3P8、OR51E2、PCAT14、GOLM1)mRNA表達水平為自變量的風(fēng)險模型.該模型顯示出良好的預(yù)測精度.

        猜你喜歡
        差異基因前列腺癌精度
        ICR鼠肝和腎毒性損傷生物標志物的篩選
        前列腺癌復(fù)發(fā)和轉(zhuǎn)移的治療
        關(guān)注前列腺癌
        認識前列腺癌
        基于RNA 測序研究人參二醇對大鼠心血管內(nèi)皮細胞基因表達的影響 (正文見第26 頁)
        前列腺癌,這些蛛絲馬跡要重視
        基于DSPIC33F微處理器的采集精度的提高
        電子制作(2018年11期)2018-08-04 03:25:38
        GPS/GLONASS/BDS組合PPP精度分析
        改進的Goldschmidt雙精度浮點除法器
        SSH技術(shù)在絲狀真菌功能基因篩選中的應(yīng)用
        蜜桃视频网站在线观看一区| 国产久视频| 用力草我小逼视频在线播放| 日本午夜剧场日本东京热| 成人亚洲一区二区三区在线| 久久99精品国产99久久6男男| 久久亚洲日本免费高清一区| 美女性色av一区二区三区| 强奸乱伦影音先锋| 风间由美性色一区二区三区| 中文字幕一区二区三区在线不卡| 黄片亚洲精品在线观看| 亚洲一区二区三区内裤视| 午夜无码片在线观看影视| 国产精品无码无片在线观看3D| 国产av熟女一区二区三区蜜臀| 最近免费中文字幕中文高清6| 亚洲av无码潮喷在线观看| 久久男人av资源网站无码| 最好的99精品色视频大全在线| 久久婷婷色香五月综合缴缴情| 亚洲乱码国产一区三区| 国产真实伦视频在线视频| 色婷婷久久综合中文蜜桃| 免费观看羞羞视频网站| 欧美三级不卡视频| 亚洲精品一区二区视频| 国产精品女老熟女一区二区久久夜 | 日韩国产成人无码av毛片蜜柚| 国产精品一区二区韩国AV| 92自拍视频爽啪在线观看| 精品国产天堂综合一区在线| 成人亚洲性情网站www在线观看 | 丰满人妻被黑人中出849 | 亚洲国产精品无码久久一线| 双乳被一左一右吃着动态图| 韩国无码精品人妻一区二| 久久av粉嫩一区二区| 高潮又爽又无遮挡又免费| 四虎成人精品国产一区a| 国产护士一区二区三区|