發(fā)表雜志:JTransl Med.
影響因子:5.530
請(qǐng)掃描下方二維碼
研究背景
流程圖
分析解讀:
數(shù)據(jù)收集
①TCGA數(shù)據(jù)庫下載轉(zhuǎn)錄組、DNA甲基化、突變數(shù)據(jù)、拷貝數(shù)變異(CNV)和臨床數(shù)據(jù)。
②CRC肝轉(zhuǎn)移數(shù)據(jù)下載自GSE28814數(shù)據(jù)集和GSE6988數(shù)據(jù)集。
③CRC治療數(shù)據(jù)來自GSE36864數(shù)據(jù)集。
④分子特征數(shù)據(jù)庫篩選出166個(gè)EMT相關(guān)基因。
EMT相關(guān)基因和lncRNA差異表達(dá)分析
①從TCGA數(shù)據(jù)集中篩選出EMT-RDGs,利用WGCNA共表達(dá)分析獲得EMT相關(guān)的lncRNAs(EMT-RlncRNAs)并進(jìn)行差異表達(dá)分析。
基因富集與功能分析
結(jié)果:
下圖A-C:CRC中EMT-RDGs的GO富集分析。
下圖D:CRC中EMT-RDGs的KEGG信號(hào)通路分析。
下圖E:預(yù)后相關(guān)基因的功能富集分析。
基于EMT-RDGs的預(yù)后模型構(gòu)建
①在TCGA數(shù)據(jù)中通過單變量cox和Lasso回歸分析篩選出9個(gè)與預(yù)后相關(guān)的EMT-RDGs,并計(jì)算每個(gè)基因的得分。通過基因風(fēng)險(xiǎn)評(píng)分構(gòu)建CRC風(fēng)險(xiǎn)預(yù)后模型。
結(jié)果:
下圖A:高風(fēng)險(xiǎn)和低風(fēng)險(xiǎn)組的生存曲線。
下圖B:風(fēng)險(xiǎn)評(píng)分的ROC曲線。
下圖C:風(fēng)險(xiǎn)評(píng)分分布和風(fēng)險(xiǎn)評(píng)分與生存時(shí)間關(guān)系的散點(diǎn)圖。
下圖D:高風(fēng)險(xiǎn)和低風(fēng)險(xiǎn)組的預(yù)后EMT-RDG和臨床參數(shù)的熱圖。
下圖E:臨床特征的單變量cox回歸森林圖。
下圖F:臨床特征的多變量cox回歸森林圖。
下圖G:多變量cox回歸分析的列線圖基線。
下圖H:3年生存率的ROC曲線。
下圖I:5年生存率的ROC曲線。
下圖J:22種免疫細(xì)胞在高危組和低危組中的浸潤水平。
結(jié)果:GSE40967數(shù)據(jù)中9個(gè)預(yù)后EMT-RGD的風(fēng)險(xiǎn)預(yù)后模型驗(yàn)證
下圖A:GSE40967數(shù)據(jù)集中風(fēng)險(xiǎn)評(píng)分和其他臨床特征的ROC曲線。
下圖B:高風(fēng)險(xiǎn)和低風(fēng)險(xiǎn)組的生存曲線。
下圖C:高風(fēng)險(xiǎn)和低風(fēng)險(xiǎn)組的預(yù)后EMT-RDG和臨床參數(shù)的熱圖。
下圖D:臨床特征的單變量cox森林圖。
下圖E:臨床特征的多變量cox森林圖。
下圖F:多變量cox回歸分析的列線圖基線。
下圖G:3年生存率的ROC曲線。
下圖H:5年生存率的ROC曲線。
下圖I:22種免疫細(xì)胞在高危組和低危組中的浸潤水平。
預(yù)后相關(guān)EMT-RDGs與免疫的關(guān)系
①ESTIMATE計(jì)算腫瘤微環(huán)境中基質(zhì)細(xì)胞和免疫細(xì)胞的純度。
②CIBERSORT計(jì)算22種免疫細(xì)胞的浸潤水平。
③利用TIMER數(shù)據(jù)庫中xCELL法計(jì)算免疫細(xì)胞富集水平,EPIC法計(jì)算免疫細(xì)胞浸潤程度,QUANTISEQ對(duì)免疫細(xì)胞進(jìn)行定量分析。
④在TIMER中獲得免疫浸潤與基因表達(dá)之間的關(guān)聯(lián)、免疫浸潤與突變狀態(tài)之間的關(guān)聯(lián)、體細(xì)胞CNV和臨床結(jié)果之間的關(guān)聯(lián)。
⑤通過ssGSEAScore計(jì)算免疫細(xì)胞和功能的得分。
結(jié)果:
下圖A:TCGA數(shù)據(jù)中ssGSEAScore比較高危組和低危組的免疫細(xì)胞評(píng)分。
下圖B:TCGA數(shù)據(jù)中ssGSEAScore比較高危組和低危組的免疫功能評(píng)分。
下圖C:GSE40967數(shù)據(jù)中ssGSEAScore比較高危組和低危組的免疫細(xì)胞得分。
下圖D:GSE40967數(shù)據(jù)中ssGSEAScore比較高危組和低危組的免疫功能評(píng)分。
預(yù)后相關(guān)EMT-RDGs與甲基化的相關(guān)性
①測(cè)試9個(gè)與預(yù)后相關(guān)的EMT-RDG的甲基化狀態(tài)和突變水平對(duì)CRC預(yù)后的影響,以及它們與藥物反應(yīng)和耐藥性的相關(guān)性。
結(jié)果:
下圖A:預(yù)后EMT-RGD的甲基化和基因表達(dá)之間的相關(guān)性與來自GCSC數(shù)據(jù)庫的spearman分析。
下圖B:GCSC數(shù)據(jù)庫中腫瘤組織與正常組織之間關(guān)于預(yù)后的EMT-RGDs甲基化水平的差異分析。
下圖C:來自GCSC數(shù)據(jù)庫的TIAM1高甲基化和低甲基化之間的總體生存差異。
下圖D:高甲基化對(duì)來自cBioPortal數(shù)據(jù)庫的預(yù)后EMT-RGD生存風(fēng)險(xiǎn)的影響。
下圖E:預(yù)后EMT-RGD的甲基化和基因表達(dá)之間的相關(guān)性與來自DNMIVD數(shù)據(jù)庫的Pearson分析。
下圖F:來自DNMIVD數(shù)據(jù)庫的腫瘤與正常組織之間關(guān)于預(yù)后的EMT-RGD的甲基化水平的差異分析。
下圖G:來自DNMIVD數(shù)據(jù)庫的九個(gè)預(yù)后EMT-RGD的重要CpG甲基化位點(diǎn)。
下圖H:Cox比例風(fēng)險(xiǎn)回歸模型基于來自DNMIVD數(shù)據(jù)庫的九個(gè)預(yù)后EMT-RGD的重要CpG甲基化位點(diǎn)。
下圖I:比較高危組和低危組的生存曲線。
預(yù)后相關(guān)EMT-RDGs表達(dá)與突變的關(guān)系
①分析CRC基因突變率和9個(gè)與預(yù)后相關(guān)的EMT-RDG之間的關(guān)聯(lián)。
結(jié)果:
下圖A:在CRC中突變Top基因的oncoplot圖。
下圖B:在CRC中突變的9個(gè)預(yù)后EMT-RGD的oncoplot圖。
下圖C:來自cBioPortal數(shù)據(jù)庫的3種類型CRC中9種預(yù)后EMT-RGD的突變類型分布。
下圖D:來自cBioPortal數(shù)據(jù)庫的CRC預(yù)后EMT-RGD的突變類型分布。
下圖E:來自cBioPortal數(shù)據(jù)庫的SNV類的比較。
下圖F:受基因突變影響的通路比例。
結(jié)果:CRC中9種預(yù)后EMT-RGD的CNV分析
下圖A:CNV與預(yù)后性EMT-RGD基因表達(dá)的相關(guān)性與來自GCSC數(shù)據(jù)庫的spearman分析。
下圖B:GCSC數(shù)據(jù)庫中結(jié)腸腺癌和直腸腺癌中9種預(yù)后EMT-RGD的雜合擴(kuò)增和缺失。
下圖C:GCSC數(shù)據(jù)庫中CRC中9種預(yù)后EMT-RGD的CNV類型分布。
預(yù)后相關(guān)的EMT-RDGs與CRC轉(zhuǎn)移的關(guān)系
①在GSE6988數(shù)據(jù)集中分析有或沒有肝轉(zhuǎn)移的CRC之間EMT-RDG的表達(dá)水平。
結(jié)果:
下圖A:GSE6988數(shù)據(jù)集中9個(gè)預(yù)后EMT-RGD在CRC肝轉(zhuǎn)移和非轉(zhuǎn)移之間的表達(dá)水平。
下圖B:GSE28814組(0:非轉(zhuǎn)移,1:轉(zhuǎn)移)中CRC轉(zhuǎn)移和非轉(zhuǎn)移之間9個(gè)預(yù)后EMT-RGD的表達(dá)水平。
預(yù)后EMT-RDGs與CRC治療的關(guān)系
①在GSE36864數(shù)據(jù)集中分析EMT-RDGs表達(dá)水平與晚期CRC治療(卡培他濱組、卡培他濱+伊立替康組、XELOX(卡培他濱+奧沙利鉑組)+貝伐珠單抗組)的關(guān)系。
結(jié)果:
下圖A:CARE數(shù)據(jù)庫中9個(gè)與預(yù)后CRC EMT-RDG耐藥性的計(jì)算分析。
下圖B:CTRP數(shù)據(jù)庫中9個(gè)CRC預(yù)后EMT-RGD的藥物篩選結(jié)果。
下圖C:GDSC數(shù)據(jù)庫中9種CRC預(yù)后的EMT-RGDs藥物篩選結(jié)果。
驗(yàn)證EMT-RDGs的表達(dá)水平
①在Oncomine數(shù)據(jù)庫、人類蛋白質(zhì)圖譜(HPA)數(shù)據(jù)庫和GEPIA數(shù)據(jù)庫驗(yàn)證EMT-RDGs的表達(dá)水平。
結(jié)果:
下圖A:GEPIA數(shù)據(jù)庫中九個(gè)預(yù)后EMT-RGD的表達(dá)驗(yàn)證。
下圖B:Oncomine數(shù)據(jù)庫中九個(gè)預(yù)后EMT-RGD的表達(dá)驗(yàn)證。
下圖C:HPA數(shù)據(jù)庫中9個(gè)預(yù)后EMT-RGD的免疫組織化學(xué)圖片。
小結(jié):
基因集泛癌分析文章列表
想發(fā)高分純生信文章,這也許就是你的機(jī)會(huì)
11.6分,最熱復(fù)合物泛癌純生信分析!想發(fā)高分純生信不容錯(cuò)過
11分+SCI 基因集泛癌純生信分析,兩個(gè)多月接受,可重復(fù)!
25天接受的基因集泛癌分析,先人一步的優(yōu)勢(shì)!
7.6分基因集泛癌純生信分析,從投稿到接收僅一個(gè)月!
單基因泛癌分析文章列表
對(duì)比兩篇單基因文章,一篇7.5分,一篇未發(fā)表
對(duì)比剛發(fā)表的2區(qū)6.5分生信文章和我們的泛癌分析內(nèi)容
被致謝!2021年7分+單基因泛癌純生信,投稿到接受30天!
適合發(fā)純生信或者生信加表達(dá)驗(yàn)證的基因
最新7.5分SCI 單基因單腫瘤生信+簡(jiǎn)單實(shí)驗(yàn)
8分+ 單基因干濕結(jié)合,非常規(guī)思路,值得學(xué)習(xí)
單基因干濕結(jié)合發(fā)7.5分SCI,2個(gè)月接收,有實(shí)驗(yàn)條件的可輕易模仿
7.5分,2個(gè)月接收的單基因生信+組化驗(yàn)證 可重復(fù)!
6.6分 不到兩個(gè)月接收 單基因生信+實(shí)驗(yàn)可以這樣做
7.7分nature子刊,單基因泛癌分析與免疫治療的關(guān)聯(lián)
6.2分 生信+實(shí)驗(yàn) 經(jīng)典套路,永不過時(shí)(一篇單基因泛癌的文章可以拆分成多個(gè)腫瘤單獨(dú)發(fā))
6.2分SCI,從投稿到接收僅38天!剛見刊的單基因泛癌純生信
非腫瘤生信文章列表
非腫瘤生信+簡(jiǎn)單pcr驗(yàn)證 快速發(fā)3分+SCI 自噬,鐵死亡,焦亡等都可以
又見非腫瘤7.5分SCI,生信分析+組化驗(yàn)證 可重復(fù)
最新7.5分非腫瘤純生信,分析一點(diǎn)不難,這種思路你想到?jīng)]有?
3+經(jīng)典非腫瘤生信套路,值得收藏學(xué)習(xí)!
聯(lián)系客服