内质网应激相关基因能预测骨肉瘤预后并与肿瘤免疫微环境相关

以内质网应激相关基因构建骨肉瘤患者的风险模型,探索其与肿瘤免疫微环境的关系。采用生物信息学分析法,训练集的转录组数据及临床数据下载于UCSC Xena数据库,验证集的相应数据下载于GEO数据库(GSE21257,GSE39058)。采用单因素COX回归分析、LASSO回归分析及多因素COX回归分析提取风险特征基因构建风险模型,使用决策曲线分析、受试者工作特征曲线分析验证模型的准确性,随后构建列线图进一步预测骨肉瘤患者预后;根据风险评分将患者分为高、低风险组,使用Kaplan-Meier生存曲线评估高、低风险组间的生存差异,对差异表达基因(Differentially expressed genes, DEGs)进行GO/KEGG联合富集分析、基因集富集分析(Gene set enrichment analysis, GSEA)及基因集变异分析(Gene set variation analysis, GSVA);采用ESTIMATE算法、微环境种群计数器(Microenvironment cell population counter, MCP counter)方法、单样本基因集富集分析(Single sample gene set enrichment analysis, ssGSEA)进行免疫分析;最终在验证集中验证上述结果。6个风险特征基因中VEGFA、PTGIS及SERPINH1与骨肉瘤患者的不良预后相关,而TMED10、MAPK10及TOR1B与与骨肉瘤患者的良好预后相关,高、低风险组患者之间具有显著生存差异;GO/KEGG联合富集分析、GSVA、GSEA结果表明DEGs与免疫状态相关;免疫分析显示高风险组具有更低的免疫评分及免疫景观;列线图进一步准确地预测了骨肉瘤患者的预后。内质网应激相关基因构建的风险模型能准确预测骨肉瘤患者预后,并与肿瘤免疫微环境相关。     

基于单细胞测序构建胶质母细胞瘤预后模型

[目的]基于单细胞测序筛选胶质母细胞瘤特征基因并构建预后模型。[方法]分析GEO数据库单细胞RNA测序数据集GSE84465,筛选出GBM细胞分化相关的差异基因。下载TCGA数据库GBM的基因表达谱和临床数据,采用Lasso回归、Cox回归分析筛选出特征基因构建预后模型,根据独立预后因素构建列线图,GSE83300作为外部验证集。基于风险评分中位数将患者分组,比较两组生存差异。[结果]通过scRNA-seq得到492个分化差异基因,经过回归分析得到基于6个基因(PLAUR、RARRES2、G0S2、MDK、SERPINE2、CD81)的预后模型。其1、3、5年ROC曲线下面积均大于0.7;KM分析显示高低风险组预后存在差异(P<0.001),GSE83300验证结果与TCGA一致。多因素Cox回归分析表明年龄和风险评分可以作为独立影响因素(P<0.01);C-Index(0.679)、校准图显示列线图预测模型有良好的拟合度。GSEA分析示高低风险组差异基因集参与细胞因子受体相互作用、抗原处理与提呈等通路。[结论]由PLAUR、RARRES2、G0S2、MDK、SERPINE...     

铁死亡相关基因与结肠癌分型及预后的关系研究

结肠癌是一种常见的发生于结肠部位的消化道恶性肿瘤,它好发于直肠与乙状结肠交界处,恶性程度高,侵袭性强,病情发展快。本研究利用铁死亡相关基因对结肠癌进行分型且对不同分型在生存时间及临床表型方面的差异进行评估,为探索该疾病的发病机制和个性化治疗提供思路。首先,从TCGA-GDC官网下载结肠癌患者的表达数据,通过查找文献检索到60个铁死亡相关基因,筛选有显著差异表达的铁死亡相关基因对结肠癌患者进行无监督聚类分型,同时比较各个分子亚型之间在生存时间和临床特征方面的差异;运用单因素Cox分析法筛选出与预后相关性较高的基因并构建Lasso回归模型,根据回归模型对患者的风险评分将患者分为高风险组和低风险组,比较两组间生存时间的差异并确定风险评分与其他临床特征之间的关联。通过单因素独立预后分析和多因素独立预后分析,筛选出影响结肠癌预后的独立因素。通过无监督聚类将样本分为两种分子亚型,两组间的生存时间差异不显著,不同分子亚型在肿瘤分期这个临床特征间存在一定的差异。使用5个与预后显著相关的基因(FDFT1、HMGCR、CARS1、AKR1C1、ALOX12)构建了Lasso回归模型,根据Lasso回归模型...     

RNA结合蛋白对肝癌预后和免疫浸润的评估价值

[目的]探索RNA结合蛋白(RBP)在肝癌预后以及免疫浸润中的评估价值。[方法]下载TCGA数据库中肝癌(LIHC)数据，筛选出差异表达的RBP基因。利用LASSO和多因素COX构建风险模型，Kaplan-Meier方法比较高低风险组患者的预后。对高低风险组患者的差异基因进行功能富集分析。利用GSVA分析高低风险组免疫浸润，并比较两组患者药物敏感性和免疫治疗反应性。[结果]总共得到330个差异基因，LASSO-COX回归模型最终纳入6个基因(UPF3B、MRPL54、ZC3H13、DHX58、PPARGC1A、EIF2AK4)。高风险组患者总体生存率劣于低风险组患者(P<0.05)。风险模型预测患者总体生存率的AUROC在TCGA和ICGC中分别为0.756和0.781。免疫浸润分析提示高风险组调节T细胞富集量高于低风险组，而NK细胞则较低。两组患者对索拉非尼、多西他赛、顺铂敏感性有明显差异(P<0.05),且低风险组对免疫治疗更敏感。[结论] RBP风险模型可有效评估肝癌患者的预后、免疫浸润状态、化疗药物及免疫治疗的敏感性，可能为肝癌患者的风险分级和免疫治疗带来新的思路。     

结合自噬相关基因表达谱和临床因素的肺腺癌预后预测模型

对肺腺癌自噬相关基因进行生物信息学分析,结合多基因预后标志和临床参数构建能够预测肺腺癌患者预后的模型。首先,对TCGA肺腺癌数据中的938个自噬相关基因进行差异分析,获得了82个差异自噬相关基因,使用单因素Cox比例风险回归模型从差异自噬相关基因中筛选出候选基因,通过 lasso回归进一步筛选出预后相关基因,分别是ARNTL2、NAPSA、ATG9B、CAPN12、MAP1LC3C和KRT81。通过多因素Cox回归分析以构建风险评分模型,根据最优cutoff值将患者分为高低风险组,生存曲线显示高低风险组之间生存差异显著,ROC曲线显示风险评分的预测能力良好,并在内、外验证集中得到验证。同时对传统的临床因素进行单因素和多因素Cox回归分析,结果显示Stage、复发和风险评分能够独立预测预后,结合这三个独立的预后参数以构建列线图模型,使用一致性指数、校准曲线评估列线图的预测能力,结果显示预测结果与实际结果之间具有良好的一致性。通过与Stage和风险评分的比较发现,列线图的预测能力表现最佳。基于肺腺癌相关的自噬基因和临床参数构建了一个列线图模型来预测肺腺癌患者的预后生存,这可能为临床医生提供了一种可靠的预后评估工具。     

急性髓系白血病铁死亡相关基因筛选及预后模型的建立

基于急性髓系白血病(Acute Myeloid Leukemia,AML)临床大数据及多组学数据库探讨铁死亡相关基因在AML中的作用,并建立铁死亡基因表达相关预后模型。整合TCGA数据库中151例AML患者和GTEx数据库中337例正常人外周血的临床和转录组数据。将Wilcoxon检验和单因素Cox分析结果取交集,筛选出预后相关差异表达基因(Differential Expression Genes, DEGs),使用Lasso回归建立基因标志物预后模型,利用受试者工作特征曲线(Receiver Operating Characteristic Curve,ROC曲线)评价预测价值,Kaplan-Meier法进行生存分析,对AML患者临床数据进行单因素和多因素Cox回归分析,使用差异基因表达分析等方法比较高、低风险患者间的组学差异,最后,利用BeatAML数据库对基因标志物进行验证。将差异基因表达分析和单因素分析结果取交集,得到13个预后相关DEGs。构建了8个基因标志物的预后评分模型,并将患者分为高、低风险两组;ROC曲线分析证实了模型良好的预测性能;生存分析提示高、低风险组患者的生存率具有显著差异;单因素分析显示年龄和风险评分与患者整体生存显著相关,多因素分析显示,年龄和风险评分是独立预后指标。在2个风险组之间筛选出384个DEGs,GO富集分析结果显示,富集的基因大多与中性粒细胞和白细胞的趋化与迁移等免疫相关分子和通路显著相关,KEGG富集通路主要与TNF信号通路、细胞因子与细胞因子受体相互作用相关。BeatAML数据库验证结果显示,5个基因与预后显著相关。铁死亡相关基因在AML中显著表达,且高风险患者预后较差,该研究对AML铁死亡相关潜在生物标志物的发现和应用奠定了一定的基础。     

肾透明细胞癌患者铁死亡相关基因的预后作用

探讨铁死亡相关基因在肾透明细胞癌患者中的表达及其预后价值。通过TCGA数据库下载KIRC的相关测序数据与检索到的铁死亡相关基因取交集,进行铁死亡相关基因的差异分析。之后利用单变量和多变量Cox回归分析,筛选具有预后价值的基因,构建预测患者生存情况的风险评分模型,并对模型进行验证。对高低风险组进行GO与KEGG通路富集,探讨风险差异的可能原因;通过ssGSEA分析,评估高低风险组间的免疫浸润情况。在KIRC患者的肿瘤组织和正常组织中,共得到21个差异的铁死亡相关基因;通过单因素Cox回归分析,获得 28 个与KIRC预后相关的基因;之后进行Lasso回归与多因素Cox回归分析,结果显示有10个基因被纳入模型,计算公式为:风险值(Risk score)=(0.024 5)×ALOX5表达值+(0.126 0)×CBS表达值+(0.199 5)×CD44表达值+(0.218 3)×CHAC1表达值+(-0.295 9)×HMGCR表达值+(0.036 7)×MT1G表达值+(0.061 4)×SLC7A11表达值+(-0.080 7)×FDFT1表达值+(0.160 3)×PEBP1表达值+(-0.220 5)×GOT1表达值。生存状态图表明,高风险组死亡病例数多于低风险组;ROC曲线表明风险评分模型具备一定预测能力;K-M生存分析显示,高风险组总体生存率低于低风险组(P=5.73×10^-13)。GO与KEGG富集分析提示,高低风险组间免疫情况及IL-17信号通路存在显著差异;进一步的ssGSEA富集显示,高低风险组间大部分免疫细胞的评分存在显著差异。基于铁死亡相关基因的预后风险评分模型可用于KIRC的预后预测,针对铁死亡相关基因设计靶点可能是治疗KIRC的一种新选择。     

免疫相关基因FGF13、GZMB、FLT3、CRLF3可能与急性髓系白血病患者预后密切相关

为了构建急性髓系白血病（AML）患者的预后风险模型,本研究通过IMMPORT数据库提取免疫基因,通过基因表达综合（GEO）数据库和AML相关转录组RNA测序数据集筛选差异表达的免疫基因,基于癌症基因组图谱（TCGA）数据库下载AML相关的表达谱数据集（TCGA-LAML）及临床数据,对差异表达的免疫基因进一步进行单因素COX回归分析和多因素COX回归分析,最终确定了关键免疫基因,随后根据免疫基因进行单基因预后分析以及免疫浸润分析。结果显示：与正常对照组相比,在AML患者的骨髓样本中,有55个免疫基因显著上调或下调,其中由4个免疫基因（FGF13、GZMB、FLT3、CRLF3）构建的风险模型能预测AML患者的预后[曲线下面积（AUC）=0.741];高风险组和低风险组之间免疫细胞的含量存在显著差异,其中FGF13、CRLF3为低风险基因,GZMB、FLT3为高风险基因。这为4个免疫基因FGF13、GZMB、FLT3、CRLF3及其构建的风险模型用于监测AML患者的预后和免疫浸润情况提供了理论基础。     

基于m5C相关基因构建三阴性乳腺癌预后预测模型及药物敏感性分析

为了探讨5-甲基胞嘧啶(5-methylcytosine,m5C)相关基因在三阴性乳腺癌(triple negative breast cancer,TNBC)患者治疗及预后中的潜在价值,构建了基于m5C相关基因的预后预测模型,用于评估TNBC患者的预后和生存状况。从基因表达总库(gene expression omnibus,GEO)数据库和癌症基因组图谱(the cancer genome atlas,TCGA)数据库中下载TNBC基因表达谱和相应的临床数据。通过Pearson分析确定了99个m5C相关基因,进一步采用单因素Cox分析鉴定出5个与预后有关的m5C相关基因(SLC6A14、BCL11A、UGT8、LMO4、PSAT1)并构建了风险评分(risk score)预测模型,根据风险评分中位值将患者划分为高风险组和低风险组。使用Kaplan-Meier(K-M)生存分析、受试者工作特征(receiver operating characteristic,ROC)曲线、多变量Cox回归分析、构建列线图和校准曲线评估了模型的预测效能。训练集和验证集的K-M生存曲线、受试者工作特征...     

肺腺癌转移相关的免疫特征开发及潜在中药预测

本研究旨在确定具有免疫相关基因的可靠预后特征,该特征可以预测预后并对肺腺癌(lung adenocarcinoma, LUAD)患者的个体化管理提供帮助。从癌症基因组图谱(The Cancer Genome Atlas, TCGA)数据库下载LUAD患者的mRNA表达谱和相应的临床数据;使用单因素COX和LASSO模型来构建预后模型;使用基于风险评分的方法开发预后特征;通过Kaplan-Meier分析比较高风险患者和低风险患者之间的总生存期(overall survival, OS), OS的独立预测因子通过单变量和多变量COX分析确定;单样本基因集富集分析(single sample gene set enrichment analysis, ssGSEA)用于评估免疫细胞浸润程度;通过LASSO和COX回归分析构建生存预后特征。根据预后特征,在OS方面将患者显著分层为高风险组和低风险组,与低风险组相比,高风险组的LUAD患者OS显著降低。通过ROC曲线分析证实了预后基因标记的预测能力。多因素COX分析显示,风险评分是OS的独立预测因子。通过免疫分析,发现了肺腺癌转移组与未转移组的不...     

The prediction of tumor and normal tissues based on the DNA methylation values of ten key sites

Abnormal DNA methylation can alter the gene expression to promote or inhibit tumorigenesis in colon adenocarcinoma (COAD). However, the finding important genes and key sites of abnormal DNA methylation which result in the occurrence of COAD is still an eventful task. Here, we studied the effects of DNA methylation in the 12 types of genomic features on the changes of gene expression in COAD, the 10 important COAD-related genes and the key abnormal DNA methylation sites were identified. The effects of important genes on the prognosis were verified by survival analysis. Moreover, it was shown that the important genes were participated in cancer pathways and were hub genes in a co-expression network. Based on the DNA methylation levels in the ten sites, the least diversity increment algorithm for predicting tumor tissues and normal tissues in seventeen cancer types are proposed. The better results are obtained in jackknife test. For example, the predictive accuracies are 94.17 %, 91.28 %, 89.04 % and 88.89 %, respectively, for COAD, rectum adenocarcinoma, pancreatic adenocarcinoma and cholangiocarcinoma. Finally, by computing enrichment score of infiltrating immunocytes and the activity of immune pathways, we found that the genes are highly correlated with immune microenvironment.     

肺腺癌预后相关miRNA生物信息学筛选及其临床意义

目的：探讨肺腺癌预后相关miRNA组学特征及其临床意义。方法：应用癌症基因组图谱（TCGA）数据库,检索人肺腺癌miRNA表达谱数据,进行差异分析,再利用Cox风险回归模型筛选预后相关miRNA;利用mirwalk分析平台,对筛选出的miRNA进行靶向调控基因预测、KEGG功能富集分析,最后,预测出预后相关miRNA的功能。结果：共筛选肺腺癌差异miRNA46个,其中,上调19个、下调27个;通过Cox生存分析筛选到预后相关miRNA有6个,即hsa-mir-21、hsa-mir-142、hsa-mir-200a高表达,hsa-mir-101、hsa-let-7c、hsa-mir-378e低表达,其中hsa-mir-21、hsa-mir-378e与肺腺癌患者不良预后有关,生存期显著缩短（P<0.05,AUC=0.618）。KEGG分析上述预后相关miRNA靶向调控基因与免疫反应通路、miRNA与癌症通路、代谢通路等有关。结论：hsa-mir-21、hsa-mir-378e与肺腺癌预后不良有关,未来经进一步临床验证有可能作为肺腺癌预后相关的分子标记物。     

Profiles of immune infiltration in lung adenocarcinoma and their clinical significant: A gene-expression–based retrospective study

Lung cancer is one of the fatal tumors. The tumor microenvironment plays a key role in regulating tumor progression. To figure out the role of tumor microenvironment in lung adenocarcinoma (LUAD), ESTIMATE algorithm was used to evaluate the immune scores in LUAD. Patients with low immune scores had a worse overall survival (OS) compared with high immune scores. Using RNA-Seq data of 489 patients in The Cancer Genome Atlas (TCGA), differentially expressed genes (DEGs) were identified between high- and low-immune score groups. Based on the DEGs, nine-gene signature was constructed by the least absolute shrinkage and selection operator Cox regression model in TCGA set. The signature demonstrated significant prognostic value in both TCGA and Gene Expression Omnibus database. Multivariate Cox regression analyses indicated that nine-genes signature was an independent prognostic factor. Subgroup analysis also revealed a robust prognostic ability of nine-gene signature. A nomogram with a C-index of 0.722 had a favorable power for predicting 3-, 5-, and 10-year survival for clinical use by integrating nine-gene signature and other clinical features. Co-expression and functional enrichment analysis showed that nine-gene signature was significantly associated with immune response and provided potential profound molecules for revealing the mechanism of tumor initiation and progression. In conclusion, we revealed the significance of immune infiltration and built a novel nine-gene signature as a reliable prognostic factor for patients with LUAD.     

Survival prediction in patients with colon adenocarcinoma via multiomics data integration using a deep learning algorithm

The present study proposed a deep learning (DL) algorithm to predict survival in patients with colon adenocarcinoma (COAD) based on multiomics integration. The survival-sensitive model was constructed using an autoencoder for DL implementation based on The Cancer Genome Atlas (TCGA) data of patients with COAD. The autoencoder framework was compared with PCA, NMF, t-SNE, and univariable Cox-PH model for identifying survival-related features. The prognostic robustness of the inferred survival risk groups was validated using three independent confirmation cohorts. Differential expression analysis, Pearson’s correlation analysis, construction of miRNA–target gene network, and function enrichment analysis were performed. Two risk groups with significant survival differences were identified in TCGA set using the autoencoder-based model (log-rank P-value = 5.51e⁻⁰⁷). The autoencoder framework showed superior performance compared with PCA, NMF, t-SNE, and the univariable Cox-PH model based on the C-index, log-rank P-value, and Brier score. The robustness of the classification model was successfully verified in three independent validation sets. There were 1271 differentially expressed genes, 10 differentially expressed miRNAs, and 12 hypermethylated genes between the survival risk groups. Among these, miR-133b and its target genes (GNB4, PTPRZ1, RUNX1T1, EPHA7, GPM6A, BICC1, and ADAMTS5) were used to construct a network. These genes were significantly enriched in ECM–receptor interaction, focal adhesion, PI3K–Akt signaling pathway, and glucose metabolism-related pathways. The risk subgroups obtained through a multiomics data integration pipeline using the DL algorithm had good robustness. miR-133b and its target genes could be potential diagnostic markers. The results would assist in elucidating the possible pathogenesis of COAD.     

Identification of 6 gene markers for survival prediction in osteosarcoma cases based on multi-omics analysis

Multiple-omics sequencing information with high-throughput has laid a solid foundation to identify genes associated with cancer prognostic process. Multiomics information study is capable of revealing the cancer occurring and developing system according to several aspects. Currently, the prognosis of osteosarcoma is still poor, so a genetic marker is needed for predicting the clinically related overall survival result. First, Office of Cancer Genomics (OCG Target) provided RNASeq, copy amount variations information, and clinically related follow-up data. Genes associated with prognostic process and genes exhibiting copy amount difference were screened in the training group, and the mentioned genes were integrated for feature selection with least absolute shrinkage and selection operator (Lasso). Eventually, effective biomarkers received the screening process. Lastly, this study built and demonstrated one gene-associated prognosis mode according to the set of the test and gene expression omnibus validation set; 512 prognosis-related genes (P < 0.01), 336 copies of amplified genes (P < 0.05), and 36 copies of deleted genes (P < 0.05) were obtained, and those genes of the mentioned genomic variants display close associations with tumor occurring and developing mechanisms. This study generated 10 genes for candidates through the integration of genomic variant genes as well as prognosis-related genes. Six typical genes (i.e. MYC, CHIC2, CCDC152, LYL1, GPR142, and MMP27) were obtained by Lasso feature selection and stepwise multivariate regression study, many of which are reported to show a relationship to tumor progressing process. The authors conducted Cox regression study for building 6-gene sign, i.e. one single prognosis-related element, in terms of cases carrying osteosarcoma. In addition, the samples were able to be risk stratified in the training group, test set, and externally validating set. The AUC of five-year survival according to the training group and validation set reached over 0.85, with superior predictive performance as opposed to the existing researches. Here, 6-gene sign was built to be new prognosis-related marking elements for assessing osteosarcoma cases’ surviving state.     

IDUA,NDST1, SAP30L,CRYBA4, and SI as novel prognostic signatures clear cell renal cell carcinoma

Carcinoma of the kidney is one of the most prevalent carcinoma worldwide. The majority types of carcinoma are clear cell renal cell carcinoma (CCRCC), which consist more than 80% of the cases. As a genetically diverse disease, identification of prognosis-related genes has utmost importance in the early diagnosis and prognosis of the CCRCC. In this study, we performed gene expression profiling to identify prognosis-related genes for CCRCC. In addition, we developed and validated a gene signature-based risk score to comprehensively assess the prognostic function of differentially expressed genes. Furthermore, we performed a ROC analysis to identify the optimal cut-off point for classification risk level of the patients. Univariate Cox regression models were used to assess the association between differentially expressed genes in relation to the prognosis of patients with different stages of CCRCC. Five genes were identified significantly differentially expressed in CCRCC and associated with their survival time, namely: IDUA, NDST1, SAP30L, CRYBA4, and SI. A 5-gene signature-based risk score was developed based on the Cox coefficient of the individual genes. The prognostic value of this risk score was validated in an internal testing data set. In summary, a gene-based risk score was identified and validated, which can predict CCRCC patient survival. The potential functions of this gene expression signature and individual differentially expressed gene as prognostic targets of CCRCC were revealed by this study. Furthermore, these findings may have important implications in the understanding of the potential therapeutic method for the CCRCC patients.     

Identification of immune microenvironment subtypes that predicted the prognosis of patients with ovarian cancer

Ovarian cancer (OC) is associated with high mortality rate. However, the correlation between immune microenvironment and prognosis of OC remains unclear. This study aimed to explore prognostic significance of OC tumour microenvironment. The OC data set was selected from the cancer genome atlas (TCGA), and 307 samples were collected. Hierarchical clustering was performed according to the expression of 756 genes. The immune and matrix scores of all immune subtypes were determined, and Kruskal-Wallis test was used to analyse the differences in the immune and matrix scores between OC samples with different immune subtypes. The model for predicting prognosis was constructed based on the expression of immune-related genes. TIDE platform was applied to predict the effect of immunotherapy on patients with OC of different immune subtypes. The 307 OC samples were classified into three immune subtypes A-C. Patients in subtype B had poorer prognosis and lower survival rate. The infiltration of helper T cells and macrophages in microenvironment indicated significant differences between immune subtypes. Enrichment analyses of immune cell molecular pathways showed that JAK–STAT3 pathway changed significantly in subtype B. Furthermore, predictive response to immunotherapy in subtype B was significantly higher than that in subtype A and C. Immune subtyping can be used as an independent predictor of the prognosis of OC patients, which may be related to the infiltration patterns of immune cells in tumour microenvironment. In addition, patients in immune subtype B have superior response to immunotherapy, suggesting that patients in subtype B are suitable for immunotherapy.     

Development and validation of a CTNNB1-associated metabolic prognostic model for hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a heterogeneous malignancy closely related to metabolic reprogramming. We investigated how CTNNB1 mutation regulates the HCC metabolic phenotype and thus affects the prognosis of HCC. We obtained the mRNA expression profiles and clinicopathological data from The Cancer Genome Atlas (TCGA), the International Cancer Genomics Consortium (ICGC) and the Gene Expression Omnibus database ( GSE14520 and GSE116174 ). We conducted gene set enrichment analysis on HCC patients with and without mutant CTNNB1 through TCGA dataset. The Kaplan-Meier analysis and univariate Cox regression analysis assisted in screening metabolic genes related to prognosis, and the prognosis model was constructed using the Lasso and multivariate Cox regression analysis. The prognostic model showed good prediction performance in both the training cohort (TCGA) and the validation cohorts (ICGC, GSE14520 , GSE116174 ), and the high-risk group presented obviously poorer overall survival compared with low-risk group. Cox regression analysis indicated that the risk score can be used as an independent predictor for the overall survival of HCC. The immune infiltration in different risk groups was also evaluated in this study to explore underlying mechanisms. This study is also the first to describe an metabolic prognostic model associated with CTNNB1 mutations and could be implemented for determining the prognoses of individual patients in clinical practice.     

Screening of immunosuppressive cells from colorectal adenocarcinoma and identification of prognostic markers

Background: Colorectal cancer (CRC) is the most common type of gastrointestinal malignant tumour. Colorectal adenocarcinoma (COAD) – the most common type of CRC – is particularly dangerous. The role of the immune system in the development of tumour-associated inflammation and cancer has received increasing attention recently.Methods: In the present study, we compiled the expression profiles of 262 patients with complete follow-up data from The Cancer Genome Atlas (TCGA) database as an experimental group and selected 65 samples from the Gene Expression Omnibus (GEO) dataset (of which 46 samples were with M0) as a verification group. First, we screened the immune T helper 17 (Th17) cells related to the prognosis of COAD. Subsequently, we identified Th17 cells-related hub genes by utilising Weighted Gene Co-expression Network Analysis (WGCNA) and Least Absolute Shrinkage and Selector Operation (LASSO) regression analysis. Six genes associated with the prognosis in patients with COAD were identified, including: KRT23, ULBP2, ASRGL1, SERPINA1, SCIN, and SLC28A2. We constructed a clinical prediction model and analysed its predictive power.Results: The identified hub genes are involved in developing many diseases and closely linked to digestive disorders. Our results suggested that the hub genes could influence the prognosis of COAD by regulating Th17 cells’ infiltration.Conclusions: These newly discovered hub genes contribute to clarifying the mechanisms of COAD development and metastasis. Given that they promote COAD development, they may become new therapeutic targets and biomarkers of COAD.