Development and validation of stable ferroptosis- and pyroptosis-related signatures in predicting prognosis and immune status in breast cancer

To develop and validate the predictive effects of stable ferroptosis- and pyroptosis-related features on the prognosis and immune status of breast cancer (BC). We retrieved as well as downloaded ferroptosis- and pyroptosis-related genes from the FerrDb and GeneCards databases. The minimum absolute contraction and selection operator (LASSO) algorithm in The Cancer Genome Atlas (TCGA) was used to construct a prognostic classifier combining the above two types of prognostic genes with differential expression, and the Integrated Gene Expression (GEO) dataset was used for validation. Seventeen genes presented a close association with BC prognosis. Thirteen key prognostic genes with prognostic value were considered to construct a new expression signature for classifying patients with BC into high- and low-risk groups. Kaplan–Meier analysis revealed a worse prognosis in the high-risk group. The receiver operating characteristic (ROC) curve and multivariate Cox regression analysis identified its predictive and independent features. Immune profile analysis showed that immunosuppressive cells were upregulated in the high-risk group, and this risk model was related to immunosuppressive molecules. We successfully constructed combined features of ferroptosis and pyroptosis in BC that are closely related to prognosis, clinicopathological and immune features, chemotherapy efficacy and immunosuppressive molecules. However, further experimental studies are required to verify these findings.     

肝细胞癌患者自噬相关基因的预后作用

构建由自噬相关基因组成的预后模型,预测肝细胞癌(HCC)患者的生存预后情况,为其个性化诊疗和临床研究提供依据.利用TCGA数据库中HCC的测序信息与人类自噬数据库联合,筛选差异表达的自噬相关基因,对其进行GO富集与KEGG通路分析;通过单因素与多因素Cox分析筛选与患者生存预后明显相关的风险基因,构建预后风险评分模型;...     

Cbl induced ubiquitination of HER2 mediate immune escape from HER2-targeted CAR-T

Breast cancer (BC) with high HER2 expression has higher recurrence rate and worse prognosis, and its immunotherapy is promising. Based on the high expression of HER2, develop Chimeric Antigen Receptor T-cell (CAR-T) and PDL-1 immunotherapy, and study the molecular pathways of related immune cells and recurrence. HER2-CAR-T cells were constructed using retroviruses, and their specific recognition and immune effects on HER2+ BC cells were verified by in vivo and in vitro experiments. PDL-1 was used as adjuvant immunotherapy, knocking down PDL-1 in tumor cells or dendritic cells, or depleted macrophages to study immune pathways. The negative regulation of HER2 by cbl was determined by IP, ubiquitination experiments, and segmented plasmids, elucidating the molecular mechanism of HER2+ BC recurrence after immunotherapy. HER2-CAR-T specifically recognizes HER2-positive tumor cells and inhibits tumor growth in vivo and in vitro, and anti-PDL1 treatment enhances the therapeutic effect of HER2-CAR-T on tumors. HER2-CART therapy eradicated solid tumors after PDL1 knockdown in dendritic cells. Immunotherapy of relapsed tumors lost HER2 expression by upregulating cbl. HER2-CAR-T shows specific recognition of HER2+ cells and can mediate immune response therapy with the cooperation of PDL-1.     

T-cell receptors provide potential prognostic signatures for breast cancer

Although T-cell receptors (TCRs) are related to the progression of breast cancer (BC), their prognostic values remain unclear. We downloaded the messenger RNA (mRNA) profiles and corresponding clinical information of 1413 BC patients from the Cancer Genome Atlas and Gene Expression Omnibus database, respectively. The different expression analysis of 104 TCRs in BC samples was performed, and the consensus clustering based on 104 TCRs was performed by using the K-mean method of R language. Univariate cox regression analysis was used to screen TCRs significantly associated with the prognosis of BC, and LASSO Cox analysis was applied to optimize key TCRs. The risk score was calculated using the prognostic model constructed based on six optimal TCRs, and multivariate Cox regression analysis was used to determine whether it was an independent prognostic signature. Finally, the nomogram was constructed to predict the overall survival of BC patients. Six optimal TCRs (ZAP70, GRAP2, NFKBIE, IFNG, NFKBIA, and PAK5), which were favorable for the prognosis of BC patients, were screened. Risk score could reliably predict the prognosis of BC patients as an independent prognostic signature. In addition, when bringing into two independent prognostic signatures, age and risk score, the nomogram model could better predict the overall survival of BC patients. Our results suggested that the poor prognosis of BC patients with high risk might be due to an immunosuppressive microenvironment. In summary, a prognostic risk model based on six TCRs was established and could efficiently predict the prognosis of BC patients.     

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.     

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

探讨铁死亡相关基因在肾透明细胞癌患者中的表达及其预后价值。通过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的一种新选择。     

Biomarkers for detection and prognosis of breast cancer identified by a functional hypermethylome screen

Breast cancer (BC) is a disease with diverse tumor heterogeneity, which challenges conventional approaches to develop biomarkers for early detection and prognosis. To identify effective biomarkers, we performed a genome-wide screen for functional methylation changes in BC, i.e., genes silenced by promoter hypermethylation, using a functionally proven gene expression approach. A subset of candidate hypermethylated genes were validated in primary BCs and tested as markers for detection and prognosis prediction of BC. We identified 33 cancer specific methylated genes and, among these, two categories of genes: (1) highly frequent methylated genes that detect early stages of BC. Within that category, we have identified the combination of NDRG2 and HOXD1 as the most sensitive (94%) and specific (90%) gene combination for detection of BC; (2) genes that show stage dependent methylation frequency pattern, which are candidates to help delineate BC prognostic signatures. For this category, we found that methylation of CDO1, CKM, CRIP1, KL and TAC1 correlated with clinical prognostic variables and was a significant prognosticator for poor overall survival in BC patients. CKM [Hazard ratio (HR) = 2.68] and TAC1 (HR = 7.73) were the strongest single markers and the combination of both (TAC1 and CKM) was associated with poor overall survival independent of age and stage in our training (HR = 1.92) and validation cohort (HR = 2.87). Our study demonstrates an efficient method to utilize functional methylation changes in BC for the development of effective biomarkers for detection and prognosis prediction of BC.     

Systematic profiling of invasion-related gene signature predicts prognostic features of lung adenocarcinoma

Due to the high heterogeneity of lung adenocarcinoma (LUAD), molecular subtype based on gene expression profiles is of great significance for diagnosis and prognosis prediction in patients with LUAD. Invasion-related genes were obtained from the CancerSEA database, and LUAD expression profiles were downloaded from The Cancer Genome Atlas. The ConsensusClusterPlus was used to obtain molecular subtypes based on invasion-related genes. The limma software package was used to identify differentially expressed genes (DEGs). A multi-gene risk model was constructed by Lasso-Cox analysis. A nomogram was also constructed based on risk scores and meaningful clinical features. 3 subtypes (C1, C2 and C3) based on the expression of 97 invasion-related genes were obtained. C3 had the worst prognosis. A total of 669 DEGs were identified among the subtypes. Pathway enrichment analysis results showed that the DEGs were mainly enriched in the cell cycle, DNA replication, the p53 signalling pathway and other tumour-related pathways. A 5-gene signature (KRT6A, MELTF, IRX5, MS4A1 and CRTAC1) was identified by using Lasso-Cox analysis. The training, validation and external independent cohorts proved that the model was robust and had better prediction ability than other lung cancer models. The gene expression results showed that the expression levels of MS4A1 and KRT6A in tumour tissues were higher than in normal tissues, while CRTAC1 expression in tumour tissues was lower than in normal tissues. The 5-gene signature prognostic stratification system based on invasion-related genes could be used to assess prognostic risk in patients with LUAD.     

Systematic construction and validation of an immune prognostic model for lung adenocarcinoma

Lung adenocarcinoma (LUAD), the most common non‐small‐cell lung cancer, is characterized by a dense lymphocytic infiltrate, which indicates that the immune system plays an active role in the development and growth of this cancer. However, no investigations to date have proposed robust models for predicting survival outcome for patients with LUAD in terms of tumour immunology. A total of 761 LUAD patients were included in this study, in which the database of The Cancer Genome Atlas (TCGA) was utilized for discovery, and the Gene Expression Omnibus (GEO) database was utilized for validation. Bioinformatics analysis and R language tools were utilized to construct an immune prognostic model and annotate biological functions. Lung adenocarcinoma showed a weakened immune phenotype compared with adjacent normal tissues. Immune‐related gene sets were profiled, an immune prognostic model based on 2 immune genes (ANLN and F2) was developed with the TCGA database to distinguish cases as having a low or high risk of unfavourable prognosis, and the model was verified with the GEO database. The model was prognostically significant in stratified cohorts, including stage I‐II, stage III‐IV and epidermal growth factor receptor (EGFR) mutant subsets, and was considered to be an independent prognostic factor for LUAD. Furthermore, the low‐ and high‐risk groups showed marked differences in tumour‐infiltrating leucocytes, tumour mutation burden, aneuploidy and PD‐L1 expression. In conclusion, an immune prognostic model was proposed for LUAD that is capable of independently identifying patients at high risk for poor survival, suggesting a relationship between local immune status and prognosis.     

Modulation of Mutational Landscape in HER2-Positive Breast Cancer after Neoadjuvant Chemotherapy

IntroductionIn early-stage HER2 positive breast cancer (BC) patients, tumor response to neoadjuvant chemotherapy (NACT) predict survival outcomes. Patients achieving less than pathological complete response (pCR) have a worse prognosis, however, this group is heterogeneous. Nowadays limited data on predictive/prognostic biomarkers in patients with residual cancer disease are available.MethodsUsing next-generation sequencing technology, we evaluated a panel of 21 cancer genes in a group of HER2 positive BC patients with residual disease after NACT. A control group of patients who achieved the pCR was selected too. The BC mutational profile was analyzed on both the tumor diagnostic biopsy and matched residual disease.ResultsOverall, the detection rate of mutations was 79% in the No-pCR group versus 90% in the pCR cohort and 98% in the residual BC. The most mutated genes were TP53 and PIK3CA. No correlations between single gene mutations and survival outcomes were found. In no-pCR cohort, 52% of patients had different mutational profile after NACT, 69% of them had an increased in the number of mutated genes. Mutational profile changes from diagnostic biopsy to residual BC were a negative prognostic factor in term of relapse free survival: recurrence probability in different gene profile sub-group was 42% vs 0% in the same profile one (P = .019).ConclusionsTreatment selective pressure on tumor cells due to NACT changed the gene mutational profile in more than half of BC patient with residual tumor disease. Treatment-induced gene mutations significantly increase the risk of relapse. Profiling primary and residual BC is a major step in order to further personalized adjuvant treatment strategy.     

Identification and Validation of a Prognostic Model Based on Three MVI-Related Genes in Hepatocellular Carcinoma

MVI has significant clinical value for treatment selection and prognosis evaluation in hepatocellular carcinoma (HCC). We aimed to construct a model based on MVI-Related Genes (MVIRGs) for risk assessment and prognosis prediction in patients with HCC. This study utilized various statistical analysis methods for prognostic model construction and validation in the Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) cohorts, respectively. In addition, immunohistochemistry and qRT-PCR were used to analyze and identify the value of the model in our cohort. After the analyses, 153 differentially expressed MVIRGs were identified, and three key genes were selected to construct a prognostic model. The high-risk group showed significantly lower overall survival (OS), and this trend was observed in all subgroups: different age groups, genders, stages, and grades. Risk score was a risk factor independent of age, gender, stage, and grade. Moreover, the ICGC cohort validated the prognostic value of the model corresponding to the TCGA. In our cohort, qRT-PCR and immunohistochemistry showed that all three genes had higher expression levels in HCC samples than in normal controls. High expression levels of genes and high-risk scores showed significantly lower recurrence-free survival (RFS) and OS, especially in MVI-positive HCC samples. Therefore, the prognostic model constructed by three MVIRGs can reliably predict the RFS and OS of patients with HCC and is valuable for guiding clinical treatment selection and prognostic assessment of HCC.     

Identification of prognostic biomarkers for breast cancer brain metastases based on the bioinformatics analysis

PurposeThe prognosis of breast cancer (BC) patients who develop into brain metastases (BMs) is very poor. Thus, it is of great significance to explore the etiology of BMs in BC and identify the key genes involved in this process to improve the survival of BC patients with BMs.Patients and methodsThe gene expression data and the clinical information of BC patients were downloaded from TCGA and GEO database. Differentially expressed genes (DEGs) in TCGA-BRCA and GSE12276 were overlapped to find differentially expressed metastatic genes (DEMGs). The protein-protein interaction (PPI) network of DEMGs was constructed via STRING database. ClusterProfiler R package was applied to perform the gene ontology (GO) enrichment analysis of DEMGs. The univariate Cox regression analysis and the Kaplan-Meier (K-M) curves were plotted to screen DEMGs associated with the overall survival and the metastatic recurrence survival, which were identified as the key genes associated with the BMs in BC. The immune infiltration and the expressions of immune checkpoints for BC patients with brain relapses and BC patients with other relapses were analyzed respectively. The correlations among the expressions of key genes and the differently infiltrated immune cells or the differentially expressed immune checkpoints were calculated. The gene set enrichment analysis (GSEA) of each key gene was conducted to investigate the potential mechanisms of key genes involved in BC patients with BMs. Moreover, CTD database was used to predict the drug-gene interaction network of key genes.ResultsA total of 154 DEGs were identified in BC patients at M0 and M1 in TCGA database. A total of 667 DEGs were identified in BC patients with brain relapses and with other relapses. By overlapping these DEGs, 17 DEMGs were identified, which were enriched in the cell proliferation related biological processes and the immune related molecular functions. The univariate Cox regression analysis and the Kaplan-Meier curves revealed that CXCL9 and GPR171 were closely associated with the overall survival and the metastatic recurrence survival and were identified as key genes associated with BMs in BC. The analyses of immune infiltration and immune checkpoint expressions showed that there was a significant difference of the immune microenvironment between brain relapses and other relapses in BC. GSEA indicated that CXCL9 and GPR171 may regulate BMs in BC via the immune-related pathways.ConclusionOur study identified the key genes associated with BMs in BC patients and explore the underlying mechanisms involved in the etiology of BMs in BC. These findings may provide a promising approach for the treatments of BC patients with BMs.     

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

基于急性髓系白血病(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铁死亡相关潜在生物标志物的发现和应用奠定了一定的基础。     

Ferroptosis-related gene signature as a prognostic marker for lower-grade gliomas

Ferroptosis is a newly discovered form of programmed cell death, which has unique biological effects on metabolism and redox biology. In this study, the prognostic value of ferroptosis-related genes was investigated in lower-grade gliomas (LGG). We downloaded the ferroptosis-related genes from the FerrDb dataset. Univariate Cox and LASSO regression analyses were applied to identify genes correlated with overall survival (OS). Subsequently, 12 ferroptosis-related genes were screened to establish the prognostic signature using stepwise multivariate Cox regression. According to the median value of risk scores, patients were divided into low- and high-risk subgroups. The Kaplan-Meier curves showed the high-risk group had a lower OS. The predictive power of the risk model was validated using the CGGA. Functional analysis revealed that the terms associated with plasma membrane receptor complex, immune response and glutamate metabolic process were primarily related to the risk model. Moreover, we established a nomogram that had a strong forecasting ability for the 1-, 3- and 5-year OS. In addition, we compared the risk scores between different clinical features. We also detected infiltration of macrophages and monocytes in different subgroups. Overall, our study identified the prognostic signature of 12 ferroptosis-related genes, which has the potential to predict the prognosis of LGG.     

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

以内质网应激相关基因构建骨肉瘤患者的风险模型,探索其与肿瘤免疫微环境的关系。采用生物信息学分析法,训练集的转录组数据及临床数据下载于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与免疫状态相关;免疫分析显示高风险组具有更低的免疫评分及免疫景观;列线图进一步准确地预测了骨肉瘤患者的预后。内质网应激相关基因构建的风险模型能准确预测骨肉瘤患者预后,并与肿瘤免疫微环境相关。     

Screening of significant biomarkers related with prognosis of liver cancer by lncRNA-associated ceRNAs analysis

This study aimed to identify significant biomarkers related to the prognosis of liver cancer using long noncoding RNA (lncRNA)-associated competing endogenous RNAs (ceRNAs) analysis. Differentially expressed mRNA and lncRNAs between liver cancer and paracancerous tissues were screened, and the functions of these mRNAs were predicted by gene ontology and pathway enrichment analyses. A ceRNA network consisting of differentially expressed mRNAs and lncRNAs was constructed. LncRNA FENDRR and lncRNA HAND2-AS1 were hub nodes in the ceRNA network. A risk score assessment model consisting of eight genes (PDE2A, ESR1, FBLN5, ALDH8A1, AKR1D1, EHHADH, ADRA1A, and GNE) associated with prognosis were developed. Multivariate Cox regression suggested that both pathologic_T and risk group could be regarded as independent prognostic factors. Furthermore, a nomogram model consisting of pathologic_T and risk group showed a good prediction ability for predicting the survival rate of liver cancer patients. The nomogram model consisting of pathologic_T and a risk score assessment model could be regarded as an independent factor for predicting prognosis of liver cancer.     

Novel gene signatures for prognosis prediction in ovarian cancer

Ovarian cancer (OV) is one of the leading causes of cancer deaths in women worldwide. Late diagnosis and heterogeneous treatment result to poor survival outcomes for patients with OV. Therefore, we aimed to develop novel biomarkers for prognosis prediction from the potential molecular mechanism of tumorigenesis. Eight eligible data sets related to OV in GEO database were integrated to identify differential expression genes (DEGs) between tumour tissues and normal. Enrichment analyses discovered DEGs were most significantly enriched in G2/M checkpoint signalling pathway. Subsequently, we constructed a multi‐gene signature based on the LASSO Cox regression model in the TCGA database and time‐dependent ROC curves showed good predictive accuracy for 1‐, 3‐ and 5‐year overall survival. Utility in various types of OV was validated through subgroup survival analysis. Risk scores formulated by the multi‐gene signature stratified patients into high‐risk and low‐risk, and the former inclined worse overall survival than the latter. By incorporating this signature with age and pathological tumour stage, a visual predictive nomogram was established, which was useful for clinicians to predict survival outcome of patients. Furthermore, SNRPD1 and EFNA5 were selected from the multi‐gene signature as simplified prognostic indicators. Higher EFNA5 expression or lower SNRPD1 indicated poorer outcome. The correlation between signature gene expression and clinical characteristics was observed through WGCNA. Drug‐gene interaction was used to identify 16 potentially targeted drugs for OV treatment. In conclusion, we established novel gene signatures as independent prognostic factors to stratify the risk of OV patients and facilitate the implementation of personalized therapies.