首先針對部分粉絲的需求統(tǒng)一回答下，弗雷賽斯目前只有2種課題服務(wù)模式：

1.提供完整技術(shù)路線的中文分析報告（包括原始數(shù)據(jù)，圖片及解讀，用戶可自行翻譯后投稿）

2.深度合作，完成從選題到發(fā)表的全過程，未達到協(xié)議目標，弗雷賽斯將全額退款！

言歸正傳，代謝一直是生命科學(xué)中久火不衰的研究方向，原因就在于任何生命活動都離不開代謝的變化，就不用說腫瘤了。有關(guān)腫瘤代謝的研究在pubmed上多如牛毛，人們把腫瘤的代謝改變叫做“重編程”。 567 那么可以想象，特定腫瘤狀態(tài)會有特定的重編程狀態(tài)。

今天ZY就為大家解讀一篇文獻有關(guān)代謝組學(xué)重編程得的文章，發(fā)表在4月份的

JOURNAL OF CLINICAL ENDOCRINOLOGY & METABOLISM（IF=5.605）上， 期刊詳情。作者以 去分化甲狀腺癌（DDTC）為研究對象，研究其代謝相關(guān)基因的表達模式。

研究背景：

甲狀腺乳頭狀癌（PTC）是去分化的甲狀腺癌（DDTC）中最常見的組織學(xué)類型[1]，大多數(shù)的PTC是不易復(fù)發(fā)且能夠被治愈的，但有5%-20%的患者會出現(xiàn)復(fù)發(fā)的情況[2]，有一些還會進一步惡化甚至死于該疾病。在BenMa等人之前發(fā)表的工作中，他們發(fā)現(xiàn)甲狀腺乳頭狀癌的去分化是導(dǎo)致其預(yù)后差的原因[3]。同時，雖然有幾項關(guān)于從基因組層面對差預(yù)后甲狀腺癌(PDTC)和甲狀腺未分化癌（ATC）的研究，這些研究表明分化良好的甲狀腺癌是通過積累一些重要基因的突變，如TERT和TP53突變，進而從逐步腫瘤進展（WDTC）到PDTC最終發(fā)展成ATC的[4-6]，但PTC去分化的分子機制仍未澄清。

之前的研究在甲狀腺癌中發(fā)現(xiàn)了代謝組學(xué)特征的改變、一些氧化應(yīng)激反應(yīng)的介質(zhì)和異常的脂質(zhì)代謝。代謝組學(xué)的分析可以有助于對惡性甲狀腺結(jié)節(jié)的準確診斷和實現(xiàn)對PTC遠端轉(zhuǎn)移和癌癥干細胞的預(yù)測。因此，在上海復(fù)旦大學(xué)癌癥研究中心頭頸外科和腫瘤科的BenMa等研究員的研究中，他們試圖去找到與PTC分化相關(guān)的代謝基因，且進一步識別和PTC預(yù)后相關(guān)的代謝基因，該工作的研究思路如Fig1所示。

Fig1. 流程圖

材料與方法

甲狀腺癌數(shù)據(jù)：FUSCC數(shù)據(jù)集作為訓(xùn)練數(shù)據(jù)集，包括5個PDTC，5個PTC和6個NT（normalthyroid）；GEO數(shù)據(jù)庫中下載甲狀腺癌芯片數(shù)據(jù)，即：GSE29265，GSE33630，GSE53157，GSE65144和GSE76039；在UCSCXena數(shù)據(jù)庫下載TCGA甲狀腺癌表達數(shù)據(jù)（n=502）。

GEO平臺數(shù)據(jù)處理：對芯片表達譜進行背景校正和RMA標準化處理，Combat去批次，探針注釋到基因，當(dāng)多個探針注釋到同一個基因，取均值。

識別與DDTC相關(guān)的代謝基因：從KEGG數(shù)據(jù)庫中收集70條代謝相關(guān)通路，整合通路中的基因，共得到1466個候選代謝基因。在FUSCC數(shù)據(jù)集中，通過分別計算PDTCv.s. PTC以及PDTC v.s. NT的Foldchange值和p值，得到在這兩個比較集中都顯著上調(diào)和顯著下調(diào)的基因（|FC|≥1.5& p<0.05），并進一步和在GSE29265和GSE33630兩套數(shù)據(jù)集中的失調(diào)基因取交集，得到在DDTC中顯著失調(diào)的基因。

識別預(yù)后相關(guān)的signature：基于單因素cox回歸分析，首先識別和預(yù)后相關(guān)的代謝基因和臨床以及其他的突變特征（如BRAFV600E等特征）；再通過多因素回歸分析進一步識別可以作為獨立預(yù)后因素的代謝基因。將在多因素cox檢驗中可以作為獨立的預(yù)后標志物的代謝基因挑選出來，作為影響甲狀腺癌預(yù)后的特征基因。

風(fēng)險模型的構(gòu)建：結(jié)合TCGA甲狀腺癌，利用上述預(yù)后相關(guān)基因構(gòu)建風(fēng)險模型，采取中位數(shù)值作為閾值將樣本劃分高低風(fēng)險組。在TCGA數(shù)據(jù)集以及整合的GEO數(shù)據(jù)集（包括GSE29265，GSE33630，GSE53157，GSE65144和GSE76039）中分別去刻畫高低風(fēng)險組，并計算該預(yù)后模型分別在TCGA數(shù)據(jù)集以及整合的GEO數(shù)據(jù)集中的曲線下面積AUC值，進而評價該預(yù)后模型的性能。

功能富集分析（GSEA）：為驗證這風(fēng)險模型中所構(gòu)成的代謝基因signature的功能，根據(jù)上述的風(fēng)險得分公式將TCGA樣本所劃分成的高低風(fēng)險標簽作為phenotypelabel，70個代謝相關(guān)的KEGG通路作為背景基因集，利用GSEA軟件隨機1000次去完成基因功能的富集分析。

代謝基因?qū)ι娴挠绊懀?/span>在Gepia數(shù)據(jù)庫（http://gepia./index.html）驗證這些代謝相關(guān)基因的高低表達對甲狀腺癌病人的DFS和OS的影響。

結(jié)果展示：

1.識別和驗證在DDTC相關(guān)的代謝基因：

在FUSCC數(shù)據(jù)集中，通過|FC|≥1.5& p<0.05，我們獲取在PDTC v.s. PTC和PDTCv.s. NT中共顯著失調(diào)（上調(diào)和下調(diào)）的121個基因（81個下調(diào)基因和40個上調(diào)基因）（如圖Fig 2b），并進一步和GSE29265和GSE33630兩套數(shù)據(jù)集中識別到的上調(diào)和下調(diào)基因取交集，最終得到了在DDTC中27個上調(diào)和38個下調(diào)基因，共65個失調(diào)基因（如圖Fig 2c所示）。

Fig2. 失調(diào)代謝基因

2.構(gòu)建預(yù)后風(fēng)險模型：基于單因素分析，我們從上述的65個失調(diào)代謝基因中找到49個和生存相關(guān)的候選基因集以及T3/T4stage、LNM和BRAFV600Emutation。進一步整合這49個基因和T3/T4stage、LNM和BRAFV600Emutation這些特征構(gòu)建多因素cox回歸模型，發(fā)現(xiàn)LPCAT2、ACOT7、HSD17B8、PDE8B以及ST3GAL1這5個基因可以作為獨立的預(yù)后標志物，并利用這5個基因構(gòu)建風(fēng)險評估模型（如下所示）。我們發(fā)現(xiàn)該風(fēng)險模型和甲狀腺癌差的預(yù)后相關(guān)（如圖Fig 3d）。通過在TCGA數(shù)據(jù)集以及整合的GEO數(shù)據(jù)集（包括GSE29265，GSE33630，GSE53157，GSE65144和GSE76039）中計算AUC值去評價該預(yù)后模型，在TCGA數(shù)據(jù)集中的AUC=0.9，在整合數(shù)據(jù)集中的AUC=0.79（如Fig3所示）。

Fig3. LPCAT2、ACOT7、HSD17B8、PDE8B以及ST3GAL1這5個基因構(gòu)成的標志物與PTC差的預(yù)后相關(guān)

3.高低風(fēng)險組病人特征刻畫：對依據(jù)上述預(yù)后模型劃分的高低風(fēng)險組病人進行諸如年齡、性別、Tstage，TERT突變特征以及相應(yīng)LPCAT2、ACOT7、HSD17B8、PDE8B以及ST3GAL1這5個基因表達分布的刻畫，通過圖中的展示，發(fā)現(xiàn)高風(fēng)險組的病人顯著發(fā)生BRAFV600E的突變。

Fig4. TCGA數(shù)據(jù)集中高低風(fēng)險組病人相關(guān)臨床和分子特征的刻畫

4.基因功能的刻畫：在TCGA樣本表達譜中，映射得到LPCAT2、ACOT7、HSD17B8、PDE8B以及ST3GAL1這5個基因表達譜數(shù)據(jù)，70條代謝相關(guān)的KEGG通路作為背景基因集，預(yù)后風(fēng)險模型對樣本劃分的高低風(fēng)險作為phenotypelabel，循環(huán)1000次，去看這5個基因所參與的生物學(xué)功能（如圖Fig 5所示）。

Fig5. 代謝基因標志物的功能.a)在TCGA數(shù)據(jù)集中GSEA刻畫代謝基因所參與的功能；b-d)LPCAT2，ACOT7和HSD17B8所參與的代謝功能。

由此，大家應(yīng)該已經(jīng)了解腫瘤代謝組學(xué)技術(shù)路線的要點了，主要分為以下四步，包括：獲得代謝基因，預(yù)后模型構(gòu)建，臨床樣本分層，基因功能探索。在實際應(yīng)用中，可以進一步靈活使用，比如代謝組學(xué)和之前講過的免疫分型，單基因技術(shù)路線的結(jié)合等等...

部分參考文獻：

AschebrookkilfoyB, Ward M H, Sabra M M, et al. Thyroid Cancer Incidence Patterns in the UnitedStates by Histologic Type, 1992–2006[J]. Thyroid, 2011, 21(2): 125-134.

SchlumbergerM. Papillary and follicular thyroid carcinoma.[J]. The New England Journal ofMedicine, 1998, 338(5): 297-306.

MaB, Xu W, Wei W, et al. Clinicopathological and Survival Outcomes ofWell-Differentiated Thyroid Carcinoma Undergoing Dedifferentiation: ARetrospective Study from FUSCC[J]. International Journal of Endocrinology,2018.

......

生信學(xué)習(xí)需要持之以恒，大家可以到目前國內(nèi)最大的生信學(xué)習(xí)社區(qū)（https:///）逛逛，尤其是神秘的SangerBox可視化生信分析軟件，相信對大家玩轉(zhuǎn)生信一定有益。生信套路千千萬，不變的是生信的精髓:差異。生信及其他科研視頻，請直接點擊“閱讀全文”。