Development and validation of a survival model for thyroid carcinoma based on autophagy-associated genes

Development of a thyroid cancer prognostic model based on the mitophagy-associated differentially expressed genes

BACKGROUND

The prevalence of thyroid cancer (ThyC), a frequent malignant tumor of the endocrine system, has been rapidly increasing over time. The mitophagy pathway is reported to play a critical role in ThyC onset and progression in many studies. This research aims to create a mitophagy-related survival prediction model for ThyC patients.

METHODS

Genes connected to mitophagy were found in the GeneCards database. The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases provided information on the expression patterns of ThyC-related genes. To identify differentially expressed genes (DEGs), R software was employed. The prognostic significance of each DEG was assessed using the prognostic K-M curve. The prognostic model was built using LASSO, ROC, univariate, and multivariate Cox regression analyses. Finally, a nomogram model was developed to predict the survival outcome of ThyC patients in the clinical setting.

RESULTS

Through differential analysis, functional enrichment analysis, and protein-protein interaction (PPI) network analysis, we screened 10 key genes related to mitophagy in ThyC. The risk model was eventually developed using LASSO and Cox regression analyses based on the six DEGs related to mitophagy. An altered expression level of a mitophagy-related prognostic gene, GGCT, was found to be the most significant one, according to the KM survival curve analysis. An immunohistochemical (IHC) investigation revealed that ThyC tissues expressed higher levels of GGCT than normal thyroid tissues. The ROC curve verified the satisfactory performance of the model in survival prediction. Multivariate Cox regression analysis showed that the pathological grade, residual tumor volume, and initial tumor lesion type were significantly linked to the prognosis. Finally, we created a nomogram to predict the overall survival rate of ThyC patients at 3-, 5-, and 7- year time points.

CONCLUSION

The nomogram risk prediction model was developed to precisely predict the survival rate of ThyC patients. The model was validated based on the most significant DEG GGCT gene expression in ThyC. This model may serve as a guide for the creation of precise treatment plans for ThyC patients.

A metabolism-related gene signature for predicting the prognosis in thyroid carcinoma

Metabolic reprogramming is one of the cancer hallmarks, important for the survival of malignant cells. We investigated the prognostic value of genes associated with metabolism in thyroid carcinoma (THCA). A prognostic risk model of metabolism-related genes (MRGs) was built and tested based on datasets in The Cancer Genome Atlas (TCGA), with univariate Cox regression analysis, LASSO, and multivariate Cox regression analysis. We used Kaplan-Meier (KM) curves, time-dependent receiver operating characteristic curves (ROC), a nomogram, concordance index (C-index) and restricted mean survival (RMS) to assess the performance of the risk model, indicating the splendid predictive performance. We established a three-gene risk model related to metabolism, consisting of PAPSS2, ITPKA, and CYP1A1. The correlation analysis in patients with different risk statuses involved immune infiltration, mutation and therapeutic reaction. We also performed pan-cancer analyses of model genes to predict the mutational value in various cancers. Our metabolism-related risk model had a powerful predictive capability in the prognosis of THCA. This research will provide the fundamental data for further development of prognostic markers and individualized therapy in THCA.

The role of lactate metabolism-related LncRNAs in the prognosis, mutation, and tumor microenvironment of papillary thyroid cancer

BACKGROUND

Lactate, a byproduct of glucose metabolism, is primarily utilized for gluconeogenesis and numerous cellular and organismal life processes. Interestingly, many studies have demonstrated a correlation between lactate metabolism and tumor development. However, the relationship between long non-coding RNAs (lncRNAs) and lactate metabolism in papillary thyroid cancer (PTC) remains to be explored.

METHODS

Lactate metabolism-related lncRNAs (LRLs) were obtained by differential expression and correlation analyses, and the risk model was further constructed by least absolute shrinkage and selection operator analysis (Lasso) and Cox analysis. Clinical, immune, tumor mutation, and enrichment analyses were performed based on the risk model. The expression level of six LRLs was tested using RT-PCR.

RESULTS

This study found several lncRNAs linked to lactate metabolism in both The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) datasets. Using Cox regression analysis, 303 lactate LRLs were found to be substantially associated with prognosis. Lasso was done on the TCGA cohort. Six LRLs were identified as independent predictive indicators for the development of a PTC prognostic risk model. The cohort was separated into two groups based on the median risk score (0.39717 -0.39771). Subsequently, Kaplan-Meier survival analysis and multivariate Cox regression analysis revealed that the high-risk group had a lower survival probability and that the risk score was an independent predictive factor of prognosis. In addition, a nomogram that can easily predict the 1-, 3-, and 5-year survival rates of PTC patients was established. Furthermore, the association between PTC prognostic factors and tumor microenvironment (TME), immune escape, as well as tumor somatic mutation status was investigated in high- and low-risk groups. Lastly, gene expression analysis was used to confirm the differential expression levels of the six LRLs.

CONCLUSION

In conclusion, we have constructed a prognostic model that can predict the prognosis, mutation status, and TME of PTC patients. The model may have great clinical significance in the comprehensive evaluation of PTC patients.

A systematic pan-cancer analysis reveals the clinical prognosis and immunotherapy value of C-X3-C motif ligand 1 (CX3CL1)

It is now widely known that C-X3-C motif ligand 1 (CX3CL1) plays an essential part in the process of regulating pro-inflammatory cells migration across a wide range of inflammatory disorders, including a number of malignancies. However, there has been no comprehensive study on the correlation between CX3CL1 and cancers on the basis of clinical features. In order to investigate the potential function of CX3CL1 in the clinical prognosis and immunotherapy, I evaluated the expression of CX3CL1 in numerous cancer types, methylation levels and genetic alterations. I found CX3CL1 was differentially expressed in numerous cancer types, which indicated CX3CL1 may plays a potential role in tumor progression. Furthermore, CX3CL1 was variably expressed in methylation levels and gene alterations in most cancers according to The Cancer Genome Atlas (TCGA). CX3CL1 was robustly associated with clinical characteristics and pathological stages, suggesting that it was related to the degree of tumor malignancy and the physical function of patients. As determined by the Kaplan-Meier method of estimating survival, high CX3CL1 expression was associated with either favorable or unfavorable outcomes depending on the different types of cancer. It suggests the correlation between CX3CL1 and tumor prognosis. Significant positive correlations of CX3CL1 expression with CD4⁺ T cells, M1 macrophage cells and activated mast cells have been established in the majority of TCGA malignancies. Which indicates CX3CL1 plays an important role in tumor immune microenvironment. Gene Ontology (GO) terms and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis suggested that the chemokine signaling pathway may shed light on the pathway for CX3CL1 to exert function. In a conclusion, our study comprehensively summarizes the potential role of CX3CL1 in clinical prognosis and immunotherapy, suggesting that CX3CL1 may represent a promising pharmacological treatment target of tumors.

Establishment and validation of individualized clinical prognostic markers for LUAD patients based on autophagy-related genes

There is considerable heterogeneity in the genomic drivers of lung adenocarcinoma, which has a dismal prognosis. Bioinformatics analysis was performed on lung adenocarcinoma (LUAD) datasets to establish a multi-autophagy gene model to predict patient prognosis. LUAD data were downloaded from The Cancer Genome Atlas (TCGA) database as a training set to construct a LUAD prognostic model. According to the risk score, a Kaplan-Meier cumulative curve was plotted to evaluate the prognostic value. Furthermore, a nomogram was established to predict the three-year and five-year survival of patients with LUAD based on their prognostic characteristics. Two genes (ITGB1 and EIF2AK3) were identified in the autophagy-related prognostic model, and the multivariate Cox proportional risk model showed that risk score was an independent predictor of prognosis in LUAD patients (HR=3.3, 95%CI= 2.3 to 4.6, P< 0.0001). The Kaplan-Meier cumulative curve showed that low-risk patients had significantly better overall (P<0.0001). The validation dataset GSE68465 further confirmed the nomogram’s robust ability to assess the prognosis of LUAD patients. A prognosis model of autophagy-related genes based on a LUAD dataset was constructed and exhibited diagnostic value in the prognosis of LUAD patients. Moreover, real-time qPCR confirmed the expression patterns of EIF2AK3 and ITGB1 in LUAD cell lines. Two key autophagy-related genes have been suggested as prognostic markers for lung adenocarcinoma.

A novel necroptosis-related gene signature associated with immune landscape for predicting the prognosis of papillary thyroid cancer

Background: Necroptosis, a type of programmed cell death, has been implicated in a variety of cancer-related biological processes. However, the roles of necroptosis-related genes in thyroid cancer yet remain unknown.

Methods: A necroptosis-related gene signature was constructed using the least absolute shrinkage and selection operator (LASSO) regression analysis and Cox regression analysis. The predictive value of the prognostic signature was validated in an internal cohort. Additionally, the single-sample gene set enrichment analysis (ssGSEA) was used to examine the relationships between necroptosis and immune cells, immunological functions, and immune checkpoints. Next, the modeled genes expressions were validated in 96 pairs of clinical tumor and normal tissue samples. Finally, the effects of modeled genes on PTC cells were studied by RNA interference approaches in vitro.

Results: In this study, the risk signature of seven necroptosis-related genes was created to predict the prognosis of papillary thyroid cancer (PTC) patients, and all patients were divided into high- and low-risk groups. Patients in the high-risk group fared worse in terms of overall survival than those in the low-risk group. The area under the curve (AUC) of the receiving operating characteristic (ROC) curves proved the predictive capability of created signature. The risk score was found to be an independent risk factor for prognosis in multivariate Cox analysis. The low-risk group showed increased immune cell infiltration and immunological activity, implying that they might respond better to immune checkpoint inhibitor medication. Next, GEO database and qRT-PCR in 96 pairs of matched tumorous and non-tumorous tissues were used to validate the expression of the seven modeled genes in PTCs, and the results were compatible with TCGA database. Finally, overexpression of IPMK, KLF9, SPATA2 could significantly inhibit the proliferation, invasion and migration of PTC cells.

Conclusion: The created necroptosis associated risk signature has the potential to have prognostic capability in PTC for patient outcome. The findings of this study could pave the way for further research into the link between necroptosis and tumor immunotherapy.

Identification of immune-related ferroptosis prognostic marker and in-depth bioinformatics exploration of multi-omics mechanisms in thyroid cancer

Background: Factors such as variations in thyroid carcinoma (THCA) gene characteristics could influence the clinical outcome. Ferroptosis and immunity have been verified to play an essential role in various cancers, and could affect the cancer patients’ prognosis. However, their relationship to the progression and prognosis of many types of THCA remains unclear.

Methods: First, we extracted prognosis-related immune-related genes and ferroptosis-related genes from 2 databases for co-expression analysis to obtain prognosis-related differentially expressed immune-related ferroptosis genes (PR-DE-IRFeGs), and screened BID and CDKN2A for building a prognostic model. Subsequently, multiple validation methods were used to test the model’s performance and compare its performance with other 4 external models. Then, we explored the mechanism of immunity and ferroptosis in the occurrence, development and prognosis of THCA from the perspectives of anti-tumor immunity, CDKN2A-related competitive endogenous RNA regulatory, copy number variations and high frequency gene mutation. Finally, we evaluated this model’s clinical practice value.

Results: BID and CDKN2A were identified as prognostic risk and protective factors, respectively. External data and qRT-PCR experiment also validated their differential expression. The model’s excellent performance has been repeatedly verified and outperformed other models. Risk scores were significantly associated with most immune cells/functions. Risk score/2 PR-DE-IRFeGs expression was strongly associated with BRAF/NRAS/HRAS mutation. Single copy number deletion of CDKN2A is associated with upregulation of CDKN2A expression and worse prognosis. The predicted regulatory network consisting of CYTOR, hsa-miRNA-873-5p and CDKN2A was shown to significantly affect prognosis. The model and corresponding nomogram have been shown to have excellent clinical practice value.

Conclusion: The model can effectively predict the THCA patients’ prognosis and guide clinical treatment. Ferroptosis and immunity may be involved in the THCA’s progression through antitumor immunity and BRAF/NRAS/HRAS mutation. CYTOR-hsa-miRNA-873-5p-CDKN2A regulatory networks and single copy number deletion of CDKN2A may also affect THCA′ progression and prognosis.

An autophagy-related lncRNA prognostic risk model for thyroid cancer

PURPOSE

Thyroid cancer (TC) is the most common malignancy of the endocrine system and its incidence is gradually rising. Research has demonstrated a close link between autophagy and thyroid cancer. We constructed a prognostic model of autophagy-related long non-coding RNA (lncRNA) in thyroid cancer and explored its prognostic value.

METHODS

The data used in this study were all obtained from The Cancer Genome Atlas (TCGA) database and the Human Autophagy Database (HADb). We construct a co-expression network by autophagy-related genes and lncRNA to obtain autophagy-related lncRNAs. After univariate Cox regression analysis and multivariate Cox regression analysis, autophagy-related lncRNAs significantly associated with prognosis were identified. Based on the risk score of lncRNA, thyroid cancer patients are divided into high-risk group and low-risk group.

RESULTS

A total of 14,142 lncRNAs and 212 autophagy-related genes (ATGs) were obtained from the TCGA database and the HADb, respectively. We performed lncRNA-ATGs correlation analysis and finally obtained 1,166 autophagy-associated lncRNAs. Subsequently, we conducted univariate Cox regression analysis and multivariate Cox regression analysis, nine autophagy-related lncRNAs (AC092279.1, AC096677.1, DOCK9-DT, LINC02454, AL136366.1, AC008063.1, AC004918.3, LINC02471 and AL162231.2) significantly associated with prognosis were identified. Based on these autophagy-related lncRNAs, a risk model was constructed. The area under the curve (AUC) of the risk score was 0.905, proving that the accuracy of risk signature was superior. In addition, multiple regression analysis showed that risk score was a significant independent prognostic risk factor for thyroid cancer.

CONCLUSION

In this study, nine autophagy-related lncRNAs in thyroid cancer were established to predict the prognosis of thyroid cancer patients.

Identification of the EMT-Related Genes Signature for Predicting Occurrence and Progression in Thyroid Cancer

BACKGROUND

The detection rate of thyroid cancer (TC) has been continuously improved due to the development of detection technology. Epithelial-mesenchymal transition (EMT) is thought to be closely related to the malignant progression of tumors. However, the relationship between EMT-related genes (ERGs) characteristics and the diagnosis and prognosis of TC patients has not been studied.

METHODS

Four datasets from Gene Expression Omnibus (GEO) were used to perform transcriptomic profile analysis. The overlapping differentially expressed ERGs (DEERGs) were analyzed using the R package "limma". Then, the hub genes, which had a higher degree, were identified by the protein-protein interaction (PPI) network. Gene expression analysis between the TC and normal data, the disease-free survival (DFS) analysis of TC patients from The Cancer Genome Atlas Thyroid Cancer (TCGA-THCA) cohort, function analysis, and immunohistochemistry (IHC) were performed to verify the importance of the hub genes. Finally, a prognostic risk scoring was constructed to predict DFS in patients with the selected genes.

RESULTS

A total of 43 DEERGs were identified and 10 DEERGs were considered hub ERGs, which had a high degree of connectivity in the PPI network. Then, the differential expressions of FN1, ITGA2, and KIT between TC and normal tissues were verified in the TCGA-THCA cohort and their protein expressions were also verified by IHC. DFS analysis indicated upregulations of FN1 expression (P<0.01) and ITGA2 expression (P<0.01) and downregulation of KIT expression (P=0.01) increased risks of decreased DFS for TCGA-THCA patients. Besides, by building a prognostic risk scoring model, we found that the DFS of TCGA-THCA patients was significantly worse in high-risk groups.

CONCLUSION

In summary, these hub ERGs were potential biomarkers for diagnosis and prognosis of TC, which can provide a basis for further exploring the efficacy of EMT in patients with TC.