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

Machine learning-based prognostic model of lactylation-related genes for predicting prognosis and immune infiltration in patients with lung adenocarcinoma

BACKGROUND

Histone lactylation is a novel epigenetic modification that is involved in a variety of critical biological regulations. However, the role of lactylation-related genes in lung adenocarcinoma has yet to be investigated.

METHODS

RNA-seq data and clinical information of LUAD were downloaded from TCGA and GEO datasets. Unsupervised consistent cluster analysis was performed to identify differentially expressed genes (DEGs) between the two clusters, and risk prediction models were constructed by Cox regression analysis and LASSO analysis. Kaplan-Meier (KM) survival analysis, ROC curves and nomograms were used to validate the accuracy of the models. We also explored the differences in risk scores in terms of immune cell infiltration, immune cell function, TMB, TIDE, and anticancer drug sensitivity. In addition, single-cell clustering and trajectory analysis were performed to further understand the significance of lactylation-related genes. We further analyzed lactate content and glucose uptake in lung adenocarcinoma cells and tissues. Changes in LUAD cell function after knockdown of lactate dehydrogenase (LDHA) by CCK-8, colony formation and transwell assays. Finally, we analyzed the expression of KRT81 in LUAD tissues and cell lines using qRT-PCR, WB, and IHC. Changes in KRT81 function in LUAD cells were detected by CCK-8, colony formation, wound healing, transwell, and flow cytometry. A nude mouse xenograft model and a KrasLSL-G12D in situ lung adenocarcinoma mouse model were used to elucidate the role of KRT81 in LUAD.

RESULTS

After identifying 26 lactylation-associated DEGs, we constructed 10 lactylation-associated lung adenocarcinoma prognostic models with prognostic value for LUAD patients. A high score indicates a poor prognosis. There were significant differences between the high-risk and low-risk groups in the phenotypes of immune cell infiltration rate, immune cell function, gene mutation frequency, and anticancer drug sensitivity. TMB and TIDE scores were higher in high-risk score patients than in low-risk score patients. MS4A1 was predominantly expressed in B-cell clusters and was identified to play a key role in B-cell differentiation. We further found that lactate content was abnormally elevated in lung adenocarcinoma cells and cancer tissues, and glucose uptake by lung adenocarcinoma cells was significantly increased. Down-regulation of LDHA inhibits tumor cell proliferation, migration and invasion. Finally, we verified that the model gene KRT81 is highly expressed in LUAD tissues and cell lines. Knockdown of KRT81 inhibited cell proliferation, migration, and invasion, leading to cell cycle arrest in the G0/G1 phase and increased apoptosis. KRT81 may play a tumorigenic role in LUAD through the EMT and PI3K/AKT pathways. In vivo, KRT81 knockdown inhibited tumor growth.

CONCLUSION

We successfully constructed a new prognostic model for lactylation-related genes. Lactate content and glucose uptake are significantly higher in lung adenocarcinoma cells and cancer tissues. In addition, KRT81 was validated at cellular and animal levels as a possible new target for the treatment of LUAD, and this study provides a new perspective for the individualized treatment of LUAD.

Construction and validation of a prognostic model based on oxidative stress-related genes in non-small cell lung cancer (NSCLC): predicting patient outcomes and therapy responses

Background

Non-small cell lung cancer (NSCLC) is a significant health concern. The prognostic value of oxidative stress (OS)-related genes in NSCLC remains unclear. The study aimed to explore the prognostic significance of OS-genes in NSCLC using extensive datasets from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO).

Methods

The research used the expression data and clinical information of NSCLC patients to develop a risk-score model. A total of 74 OS-related differentially expressed genes (DEGs) were identified by comparing NSCLC and control samples. Univariate Cox and least absolute shrinkage and selection operator (LASSO) regression analyses were employed to identify the prognostic biomarkers. A risk-score model was constructed and validated with receiver operating characteristic (ROC) curves in TCGA and GSE72094 datasets. The model's accuracy was further verified by univariate and multivariate Cox regression.

Results

The identified biomarkers, including lactate dehydrogenase A (LDHA), protein tyrosine phosphatase receptor type N (PTPRN), and transient receptor potential cation channel subfamily A (TRPA1) demonstrated prognostic significance in NSCLC. The risk-score model showed good predictive accuracy, with 1-year area under the curves (AUC) of 0.661, 3-year AUC of 0.648, and 5-year AUC of 0.634 in the TCGA dataset, and 1-year AUC of 0.643, 3-year AUC of 0.648, and 5-year AUC of 0.662 in the GSE72094 dataset. A nomogram integrating risk score and tumor node metastasis (TNM) stage was developed. The signature effectively distinguished between patient responses to immunotherapy. High-risk groups were characterized by an immunosuppressive microenvironment and an increased tumor mutational burden (TMB), marked by a higher incidence of mutations in genes such as TP53, DCP1B, ELN, and MAGI2. Organoid drug sensitivity testing revealed that NSCLC patients with a low-risk score responded better to chemotherapy.

Conclusions

This study successfully developed a robust model for predicting patient prognosis in NSCLC, highlighting the critical prognostic value of OS-genes. These findings hold significant potential to refine treatment strategies, and enhance survival outcomes for NSCLC patients. By enabling a personalized therapeutic approach tailored to individual risk scores, this model may facilitate more precise decisions concerning immunotherapy and chemotherapy, thereby optimizing patient management and treatment efficacy.

Identification of oxidative stress signatures of lung adenocarcinoma and prediction of patient prognosis or treatment response with single-cell RNA sequencing and bulk RNA sequencing data

Lung adenocarcinoma (LUAD), the most common subtype of lung cancer globally, has seen improved prognosis with advancements in diagnostic, surgical, radiotherapy, and molecular therapy techniques, while its 5-year survival rate remains low. Molecular biomarkers provide prognostic value. Oxidative stress factors, such as reactive nitrogen species and ROS, are crucial in various stages of tumor progression, influencing cell transformation, proliferation, angiogenesis, and metastasis. ROS demonstrate dual roles, affecting tumor cells, hypoxia sensitivity, and the microenvironment. Comprehensive analysis of oxidative stress in LUAD has not been conducted to date. Therefore, we systematically investigated the regulatory patterns of oxidative stress in LUAD based on oxidative stress-related genes and correlated these patterns with cellular infiltration characteristics of the tumor immune microenvironment. The model utilizes single-factor Cox analysis to screen key differential genes with prognostic value and employs least absolute shrinkage and selection operator (LASSO) penalized Cox regression analysis to construct a prognostic-related prediction model. Ten candidate genes were selected based on this model. The risk score was constructed using the coefficients and expression levels of these ten genes. Furthermore, the impact of this risk score on overall survival (OS) was determined. Two genes with the most significant differential expression, SFTPB and S100P, were selected through qRT-PCR. Cell experiments including CCK-8, Edu, transwell assays confirmed their effects on lung cancer cells growth, consistent with the results of bioinformatics analysis. These findings suggested that this model held potential clinical value for evaluating the prognosis of lung adenocarcinoma.

CRTAC1 identified as a promising diagnosis and prognostic biomarker in lung adenocarcinoma

CRTAC1, one of the pyroptosis-related genes, has been identified as a protective factor in certain kinds of cancer, such as gastric adenocarcinoma and bladder cancer. The study aimed to investigate the role of CRTAC1 in lung adenocarcinoma (LUAD). LUAD datasets were obtained from Gene Expression Omnibus (GEO) database and The Cancer Genome Atlas (TCGA), pyroptosis-related genes from GeneCard. Limma package used to find differentially expressed genes (DEGs), least absolute shrinkage and selection operator (LASSO) regression and weighted genes co-expression network analysis (WGCNA) to identify CRTAC1 as hub gene. CRTAC1 expression was confirmed in a real-world cohort using quantitative polymerase chain reaction (qPCR) and Western Blot (WB) analyses. Cellular experiments were conducted to investigate CRTAC1's potential oncogenic mechanisms. CRTAC1 mRNA expression was significantly lower in LUAD tissues (p < 0.05) and showed high accuracy in diagnosing LUAD. Reduced CRTAC1 expression was associated with a poor prognosis. Higher CRTAC1 expression correlated with increased immune cell infiltration. Individuals with high CRTAC1 expression showed increased drug sensitivity. Additionally, qPCR and WB analyses showed that CRTAC1 expression was lower in tumor tissue compared to adjacent normal tissue at both the RNA and protein levels. Upregulation of CRTAC1 significantly inhibited LUAD cell proliferation, invasion, and migration in cellular experiments. CRTAC1 has the potential to serve as a diagnostic and prognostic biomarker in LUAD.

The Tumor Stemness Indice mRNAsi can Act as Molecular Typing Tool for Lung Adenocarcinoma

Due to the high heterogeneity, lung adenocarcinoma (LUAD) cannot be distinguished into precise molecular subtypes, thereby resulting in poor therapeutic effect and low 5-year survival rate clinically. Although the tumor stemness score (mRNAsi) has been shown to accurately characterize the similarity index of cancer stem cells (CSCs), whether mRNAsi can serve as an effective molecular typing tool for LUAD isn't reported to date. In this study, we first demonstrate that mRNAsi is significantly correlated with the prognosis and disease degree of LUAD patients, i.e., the higher the mRNAsi, the worse the prognosis and the higher the disease degree. Second, we identify 449 mRNAsi-related genes based on both weighted gene co-expression network analysis (WGCNA) and univariate regression analysis. Third, our results display that 449 mRNAsi-related genes can accurately distinguish the LUAD patients into two molecular subtypes: ms-H subtype (with high mRNAsi) and ms-L subtype (with low mRNAsi), particularly the ms-H subtype has a worse prognosis. Remarkably, significant differences in clinical characteristics, immune microenvironment, and somatic mutation exist between the two molecular subtypes, which might lead to the poorer prognosis of the ms-H subtype patients than that of the ms-L subtype ones. Finally, we establish a prognostic model containing 8 mRNAsi-related genes, which can effectively predict the survival rate of LUAD patients. Taken together, our work provides the first molecular subtype related to mRNAsi in LUAD, and reveals that these two molecular subtypes, the prognostic model and marker genes may have important clinical value for effectively monitoring and treating LUAD patients.

Investigating subtypes of lung adenocarcinoma by oxidative stress and immunotherapy related genes

Lung adenocarcinoma (LUAD) is one of the most widespread and fatal types of lung cancer. Oxidative stress, resulting from an imbalance in the production and accumulation of reactive oxygen species (ROS), is considered a promising therapeutic target for cancer treatment. Currently, immune checkpoint blockade (ICB) therapy is being explored as a potentially effective treatment for early-stage LUAD. In this research, we aim to identify distinct subtypes of LUAD patients by investigating genes associated with oxidative stress and immunotherapy. Additionally, we aim to propose subtype-specific therapeutic strategies. We conducted a thorough search of the Gene Expression Omnibus (GEO) datasets. From this search, we pinpointed datasets that contained both expression data and survival information. We selected genes associated with oxidative stress and immunotherapy using keyword searches on GeneCards. We then combined expression data of LUAD samples from both The Cancer Genome Atlas (TCGA) and 11 GEO datasets, forming a unified dataset. This dataset was subsequently divided into two subsets, Dataset_Training and Dataset_Testing, using a random bifurcation method, with each subset containing 50% of the data. We applied consensus clustering (CC) analysis to identify distinct LUAD subtypes within the Dataset_Training. Molecular variances associated with oxidative stress levels, the tumor microenvironment (TME), and immune checkpoint genes (ICGs) were then investigated among these subtypes. Employing feature selection combined with machine learning techniques, we constructed models that achieved the highest accuracy levels. We validated the identified subtypes and models from Dataset_Training using Dataset_Testing. A hub gene with the highest importance values in the machine learning model was identified. We then utilized virtual screening to discover potential compounds targeting this hub gene. In the unified dataset, we integrated 2,154 LUAD samples from TCGA-LUAD and 11 GEO datasets. We specifically selected 1,311 genes associated with immune and oxidative stress processes. The expression data of these genes were then employed for subtype identification through CC analysis. Within Dataset_Training, two distinct subtypes emerged, each marked by different levels of immune and oxidative stress pathway values. Consequently, we named these as the OX+ and IM+ subtypes. Notably, the OX+ subtype showed increased oxidative stress levels, correlating with a worse prognosis than the IM+ subtype. Conversely, the IM+ subtype demonstrated enhanced levels of immune pathways, immune cells, and ICGs compared to the OX+ subtype. We reconfirmed these findings in Dataset_Testing. Through gene selection, we identified an optimal combination of 12 genes for predicting LUAD subtypes: ACP1, AURKA, BIRC5, CYC1, GSTP1, HSPD1, HSPE1, MDH2, MRPL13, NDUFS1, SNRPD1, and SORD. Out of the four machine learning models we tested, the support vector machine (SVM) stood out, achieving the highest area under the curve (AUC) of 0.86 and an accuracy of 0.78 on Dataset_Testing. We focused on HSPE1, which was designated as the hub gene due to its paramount importance in the SVM model, and computed the docking structures for four compounds: ZINC3978005 (Dihydroergotamine), ZINC52955754 (Ergotamine), ZINC150588351 (Elbasvir), and ZINC242548690 (Digoxin). Our study identified two subtypes of LUAD patients based on oxidative stress and immunotherapy-related genes. Our findings provided subtype-specific therapeutic strategies.

CRTAC1 enhances the chemosensitivity of non-small cell lung cancer to cisplatin by eliciting RyR-mediated calcium release and inhibiting Akt1 expression

Sensitivity to platinum-based combination chemotherapy is associated with a favorable prognosis in patients with non-small cell lung cancer (NSCLC). Here, our results obtained from analyses of the Gene Expression Omnibus database of NSCLC patients showed that cartilage acidic protein 1 (CRTAC1) plays a role in the response to platinum-based chemotherapy. Overexpression of CRTAC1 increased sensitivity to cisplatin in vitro, whereas knockdown of CRTAC1 decreased chemosensitivity of NSCLC cells. In vivo mouse experiments showed that CRTAC1 overexpression increased the antitumor effects of cisplatin. CRTAC1 overexpression promoted NFAT transcriptional activation by increasing intracellular Ca2+ levels, thereby inducing its regulated STUB1 mRNA transcription and protein expression, accelerating Akt1 protein degradation and, in turn, enhancing cisplatin-induced apoptosis. Taken together, the present results indicate that CRTAC1 overexpression increases the chemosensitivity of NSCLC to cisplatin treatment by inducing Ca2+-dependent Akt1 degradation and apoptosis, suggesting the potential of CRTAC1 as a biomarker for predicting cisplatin chemosensitivity. Our results further reveal that modulating the expression of CRTAC1 could be a new strategy for increasing the efficacy of cisplatin in chemotherapy of NSCLC patients.

Integrative analysis of TP53 mutations in lung adenocarcinoma for immunotherapies and prognosis

BACKGROUND

The TP53 tumor suppressor gene is one of the most mutated genes in lung adenocarcinoma (LUAD) and plays a vital role in regulating the occurrence and progression of cancer. We aimed to elucidate the association between TP53 mutations, response to immunotherapies and the prognosis of LUAD.

METHODS

Genomic, transcriptomic, and clinical data of LUAD were downloaded from The Cancer Genome Atlas (TCGA) dataset. Gene ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, gene set enrichment analysis (GSEA). Gene set variation analysis (GSVA) were performed to determine the differences in biological pathways. A merged protein-protein interaction (PPI) network was constructed and analyzed. MSIpred was used to analyze the correlation between the expression of the TP53 gene, tumor mutation burden (TMB) and tumor microsatellite instability (MSI). CIBERSORT was used to calculate the abundance of immune cells. Univariate and multivariate Cox regression analyses were used to determine the prognostic value of TP53 mutations in LUAD.

RESULTS

TP53 was the most frequently mutated in LUAD, with a mutational frequency of 48%. GO and KEGG enrichment analysis, GSEA, and GSVA results showed a significant upregulation of several signaling pathways, including PI3K-AKT mTOR (P < 0.05), Notch (P < 0.05), E2F target (NES = 1.8, P < 0.05), and G2M checkpoint (NES = 1.7, P < 0.05). Moreover, we found a significant correlation between T cells, plasma cells, and TP53 mutations (R² < 0.01, P = 0.040). Univariate and multivariate Cox regression analyses revealed that the survival prognosis of LUAD patients was related to TP53 mutations (Hazard Ratio (HR) = 0.72 [95% CI, 0.53 to 0.98], P < 0.05), cancer status (P < 0.05), and treatment outcomes (P < 0.05). Lastly, the Cox regression models showed that TP53 exhibited good power in predicting three- and five-year survival rates.

CONCLUSIONS

TP53 may be an independent predictor of response to immunotherapy in LUAD, and patients with TP53 mutations have higher immunogenicity and immune cell infiltration.

The oncogenic role of TFAP2A in bladder urothelial carcinoma via a novel long noncoding RNA TPRG1-AS1/DNMT3A/CRTAC1 axis

BACKGROUND

Overexpression of TFAP2A has been linked to increased lymph node metastasis in basal-squamous bladder cancer. However, its downstream targets in bladder urothelial carcinoma (BLCA), the most malignant cancer of the urinary tract, remain unclear. In the current study, we aim to explore the function and mechanism of TFAP2A in BLCA.

METHODS

TFAP2A expression and the prognostic significance in BLCA was analyzed using TCGA and GTEX projects. TFAP2A was knocked-down in BLCA cells to study its impact on glucose uptake, lactate and ATP production, expression of HK2, and the number of vascular meshes formed by HUVEC. The target long noncoding RNAs (lncRNAs) of TFAP2A were predicted by bioinformatics tools, followed by ChIP-qPCR and luciferase assays. The downstream targets of TPRG1-AS1 were analyzed by microarray analysis. Rescue experiments were conducted for validation.

RESULTS

TFAP2A upregulation in BLCA predicted dismal survival of patients. Loss of TFAP2A inhibited glycolysis (as evidenced by reduced glucose uptake, lactate, ATP production, and the expression of HK2) and angiogenesis (decreased number of vascular meshes formed by HUVEC). TFAP2A promoted the transcription of TPRG1-AS1. TPRG1-AS1 reversed the inhibitory effect of TFAP2A knockdown on glycolysis and angiogenesis in BLCA cells. TPRG1-AS1 inhibited the transcription of CRTAC1 by recruiting a DNA methyltransferase to the promoter of CRTAC1 and increasing the DNA methylation of its promoter. CRTAC1 inhibited glycolysis and angiogenesis in BLCA cells. TFAP2A silencing curbed tumor growth in vivo via the TPRG1-AS1/CRTAC1 axis.

CONCLUSION

TFAP2A reduces CRTAC1 expression by promoting TPRG1-AS1 transcription, thereby expediting BLCA glycolysis and angiogenesis.

CD8+ T cell-associated genes MS4A1 and TNFRSF17 are prognostic markers and inhibit the progression of colon cancer

Background

Colon cancer (CC) is among the top three diseases with the highest morbidity and mortality rates worldwide. Its increasing incidence imposes a major global health burden. Immune checkpoint inhibitors, such as anti-PD-1 and anti-PD-L1, can be used for the treatment of CC; however, most patients with CC are resistant to immunotherapy. Therefore, identification of biomarkers that can predict immunotherapy sensitivity is necessary for selecting patients with CC who are eligible for immunotherapy.

Methods

Differentially expressed genes associated with the high infiltration of CD8+ T cells were identified in CC and para-cancerous samples via bioinformatic analysis. Kaplan–Meier survival analysis revealed that MS4A1 and TNFRSF17 were associated with the overall survival of patients with CC. Cellular experiments were performed for verification, and the protein expression of target genes was determined via immunohistochemical staining of CC and the adjacent healthy tissues. The proliferation, migration and invasion abilities of CC cells with high expression of target genes were determined via in vitro experiments.

Results

Differential gene expression, weighted gene co-expression and survival analyses revealed that patients with CC with high expression of MS4A1 and TNFRSF17 had longer overall survival. The expression of these two genes was lower in CC tissues than in healthy colon tissues and was remarkably associated with the infiltration of various immune cells, including CD8+ T cells, in the tumour microenvironment (TME) of CC. Patients with CC with high expression of MS4A1 and TNFRSF17 were more sensitive to immunotherapy. Quantitative reverse transcription-polymerase chain reaction, western blotting and immunohistochemical staining validated the differential expression of MS4A1 and TNFRSF17. In addition, Cell Counting Kit-8, wound healing and transwell assays revealed that the proliferation, migration and invasion abilities of CC cells were weakened after overexpression of MS4A1 and TNFRSF17.

Conclusions

The core genes MS4A1 and TNFRSF17 can be used as markers to predict the sensitivity of patients with CC to immunotherapy and have potential applications in gene therapy to inhibit CC progression.

Molecular Typing Based on Oxidative Stress Genes and Establishment of Prognostic Characteristics of 7 Genes in Lung Adenocarcinoma

Oxidative stress could maintain different biological processes in human cancer. However, the effect of oxidative stress on lung adenocarcinoma (LUAD) should be studied. This study analyzed the expression and clinical importance of oxidative stress in LUAD in detail. The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) were employed for obtaining LUAD expression profiles. Based on oxidative stress-related genes, molecular subtypes substantially correlated with the LUAD prognosis were discovered with ConsensusClusterPlus. Differentially expressed genes (DEGs) among subtypes were found using the Limma software package. Least absolute shrinkage and selection operator- (Lasso-) Cox analysis was employed to create the polygenic risk model. RiskScore and clinically relevant features were used to create nomograms. By utilizing oxidative stress-related genes and reliable clustering, stable molecular subtypes were first discovered. The prognosis, clinical characteristics, route characteristics, and immunological characteristics of these three molecular subtypes were all different. Subsequently, by using differential expression genes among molecular subtypes and Lasso, 7 main genes linked with the oxidative stress phenotype were discovered. A prognostic risk model was also built on the basis of major genes associated with the oxidative stress phenotype. The model demonstrated a high level of resilience and was unaffected by clinical-pathological features. It played a stable predictive role in independent datasets. Ultimately, to improve the prognosis model and survival prediction, RiskScore (RS) was combined with clinicopathological variables, and a decision tree model was used. The model exhibited a high prediction accuracy as well as the ability to predict survival. This research found that oxidative stress-related genes have a major involvement in the onset and progression of LUAD and that they may influence LUAD susceptibility to immunotherapy and standard chemotherapy. Furthermore, the identified risk models for 7 genes linked with oxidative stress exhibited could assist clinical treatment decisions and prognosis prediction. The classifier could be used as a molecular diagnostic tool for assessing LUAD patients' prognosis risk.

Identification of N7-methylguanosine related signature for prognosis and immunotherapy efficacy prediction in lung adenocarcinoma

Background

Lung adenocarcinoma (LUAD) is one of the most frequent causes of tumor-related mortality worldwide. Recently, the role of N7-methylguanosine (m7G) in tumors has begun to receive attention, but no investigation on the impact of m7G on LUAD. This study aims to elucidate the significance of m7G on the prognosis and immunotherapy in LUAD.

Methods

Consensus clustering was employed to determine the molecular subtype according to m7G-related regulators extracted from The Cancer Genome Atlas (TCGA) database. Survival, clinicopathological features and tumor mutational burden (TMB) analysis were applied to research molecular characteristics of each subtype. Subsequently, “limma” package was used to screen differentially expressed genes (DEGs) between subtypes. In the TCGA train cohort (n = 245), a prognostic signature was established by univariate Cox regression, lasso regression and multivariate Cox regression analysis according to DEGs and survival analysis was employed to assess the prognosis. Then the prognostic value of the signature was verified by TCGA test cohort (n = 245), TCGA entire cohort (n = 490) and GSE31210 cohort (n = 226). Moreover, the association among immune infiltration, clinical features and the signature was investigated. The immune checkpoints, TMB and tumor immune dysfunction and exclusion (TIDE) were applied to predict the immunotherapy response.

Results

Two novel molecular subtypes (C1 and C2) of LUAD were identified. Compared to C2 subtype, C1 subtype had poorer prognosis and higher TMB. Subsequently, the signature (called the “m7G score”) was constructed according to four key genes (E2F7, FAM83A, PITX3, and HOXA13). The distribution of m7G score were significantly different between two molecular subtypes. The patients with lower m7G score had better prognosis in TCGA train cohort and three verification cohort. The m7G score was intensively related to immune infiltration. Compared with the lower score, the higher m7G score was related to remarkable upregulation of the PD-1 and PD-L1, the higher TMB and the lower TIDE score.

Conclusion

This study established a m7G-related signature for predicting prognosis and immunotherapy in LUAD, which may contribute to the development of new therapeutic strategies for LUAD.

The prognostic value of immune-related genes AZGP1, SLCO5A1, and CTF1 in Uveal melanoma

Objective

Uveal melanoma (UM) is an aggressive malignancy with a poor prognosis and no available effective treatment. Therefore, exploring a potential prognostic marker for UM could provide new possibilities for early detection, recurrence, and treatment.

Methods

In this study, we used “ConsensusClusterPlus” to classify patients with UM into subgroups, screened for significant differences in immune prognostic factors between subgroups, selected three genes using LASSO (Least absolute shrinkage and selection operator) regression to construct a risk model, and performed tumor immune cell infiltration analysis on the risk model. infiltration analysis, and then verified the heterogeneous role of the 3 core genes in other cancers by pan-cancer analysis and validate its expression by RT-qPCR in normal and tumor cells.

Results

We consistently categorized 80 UM patients into two subgroups after the immunogenetic set, where the UM1 subgroup had a better prognosis than the UM2 subgroup, and used 3 immune-related genes AZGP1, SLCO5A1, and CTF1 to derive risk scores as independent prognostic markers and predictors of UM clinicopathological features. We found significant differences in overall survival (OS) between low- and high-risk groups, and prognostic models were negatively correlated with B cell and myeloid dendritic cell and positively correlated with CD8+ T cell AZGP1 and CTF1 were significantly upregulated in UM cells compared with normal UM cells.

Conclusion

Immunogens are significantly associated with the prognosis of UM, and further classification based on genetic characteristics may help to develop immunotherapeutic strategies and provide new approaches to develop customized treatment strategies for patients.

Clustering analysis and prognostic signature of lung adenocarcinoma based on the tumor microenvironment

Because of immunotherapy failure in lung adenocarcinoma, we have tried to find new potential biomarkers for differentiating different tumor subtypes and predicting prognosis. We identified two subtypes based on tumor microenvironment-related genes in this study. We used seven methods to analyze the immune cell infiltration between subgroups. Further analysis of tumor mutation load and immune checkpoint expression among different subgroups was performed. The least absolute shrinkage and selection operator Cox regression was applied for further selection. The selected genes were used to construct a prognostic 14-gene signature for LUAD. Next, a survival analysis and time-dependent receiver operating characteristics were performed to verify and evaluate the model. Gene set enrichment analyses and immune analysis in risk groups was also performed. According to the expression of genes related to the tumor microenvironment, lung adenocarcinoma can be divided into cold tumors and hot tumors. The signature we built was able to predict prognosis more accurately than previously known models. The signature-based tumor microenvironment provides further insight into the prediction of lung adenocarcinoma prognosis and may guide individualized treatment.

A novel prognostic signature of metastasis-associated genes and personalized therapeutic strategy for lung adenocarcinoma patients

Lung adenocarcinoma (LUAD) is a highly invasive and metastatic malignant tumor with high morbidity and mortality. This study aimed to construct a prognostic signature for LUAD patients based on metastasis-associated genes (MAGs). RNA expression profiles were downloaded from the Cancer Genome Atlas (TCGA) database. RRA method was applied to identify differentially expressed MAGs. A total of 192 significantly robust MAGs were determined among seven GEO datasets. MAGs were initially selected through the Lasso Cox regression analysis and 6 MAGs were included to construct a prognostic signature model. Transcriptome profile, patient prognosis, correlation between the risk score and clinicopathological features, immune cell infiltration characteristics, immunotherapy sensitivity and chemotherapy sensitivity differed between low- and high-risk groups after grouping according to median risk score. The reliability and applicability of the signature were further validated in the GSE31210, GSE50081 and GSE68465 cohort. CMap predicted 62 small molecule drugs on the base of the prognostic MAGs. Targeted drug staurosporine had hydrogen bonding with Gln-172 of SLC2A1, which is one of MAGs. Staurosporine could inhibit cell migration in A549 and H1299. We further verified mRNA and protein expression of 6 MAGs in A549 and H1299. The signature can serve as a promising prognostic tool and may provide a novel personalized therapeutic strategy for LUAD patients.

Novel Pyroptosis-Related Gene Signatures Identified as the Prognostic Biomarkers for Bladder Carcinoma

Background

Bladder carcinoma (BLCA) is a common malignant tumor with high morbidity and mortality in the urinary system. Pyroptosis is a pattern of programmed cell death that is closely associated with progression of tumors. Therefore, it is significant to probe the expression of pyroptosis-related genes (PRGs) in BLCA.

Methods

The differentially expressed genes in normal and BLCA tissues were first obtained from the Cancer Genome Atlas (TCGA) database analysis, as well as PRGs from the National Center for Biotechnology Information (NCBI) database, intersecting to obtain differentially expressed pyroptosis-related genes (DEPRGs) in BLCA. With the construction of a prognostic model of pyroptosis by regression analysis, we derived and validated key genes, which were ascertained as a separate prognostic marker by individual prognostic and clinical relevance analysis. In addition, we gained six immune cells from the Tumor Immune Evaluation Resource (TIMER) website and analyzed the relationship between pyroptosis prognostic genes and immune infiltration.

Result

Our results revealed that 31 DEPRGs were available by comparing normal and BLCA tissues with |log2 (fold change, FC)| &gt; 0.5 and FDR &lt;0.05. Four key genes (CRTAC1, GSDMB, AIM2, and FOXO3) derived from the pyroptosis prognostic model were experimentally validated for consistent expression in BLCA patients. Following risk scoring, the low-risk group of BLCA patients had noticeably higher overall survival (OS) than the high-risk group (p &lt; 0.001). Risk score was still an independent prognostic factor (HR = 1.728, 95% CI =1.289–2.315, p &lt; 0.001). In addition, we found remarkable correlations among the expression of pyroptosis-related prognostic genes and the immune infiltration of CD4+ T cells, CD8+ T cells, B cells, dendritic cells, macrophages, and neutrophils.

Conclusion

Genes (CRTAC1, GSDMB, AIM2, and FOXO3) associated with pyroptosis are potential BLCA prognostic biomarkers that act as an essential part in the predictive prognosis of survival and immunotherapy of BLCA.

Predicting Chemo-Radiotherapy Sensitivity With Concordant Survival Benefit in Non-Small Cell Lung Cancer via Computed Tomography Derived Radiomic Features

Background

To identify a computed tomography (CT) derived radiomic signature for the options of concurrent chemo-radiotherapy (CCR) in patients with non-small cell lung cancer (NSCLC).

Methods

A total of 226 patients with NSCLC receiving CCR were enrolled from public dataset, and allocated to discovery and validation sets based on patient identification number. Using CT images of 153 patients in the discovery dataset, we pre-selected a list of radiomic features significantly associated with 5-year survival rate and adopted the least absolute shrinkage and selection operator regression to establish a predictive radiomic signature for CCR treatment. We performed transcriptomic analyzes of the signature, and evaluated its association with molecular lesions and immune landscapes in a dataset with matched CT images and transcriptome data. Furthermore, we identified CCR resistant genes positively correlated with resistant scores of radiomic signature and screened essential resistant genes for NSCLC using genome-scale CRIPSR data. Finally, we combined DrugBank and Genomics of Drug Sensitivity in Cancer databases to excavate candidate therapeutic agents for patients with CCR resistance, and validated them using the Connectivity Map dataset.

Results

The radiomic signature consisting of nine features was established, and then validated in the dataset of 73 patients receiving CCR log-rank P = 0.0005, which could distinguish patients into resistance and sensitivity groups, respectively, with significantly different 5-year survival rate. Furthermore, the novel proposed radiomic nomogram significantly improved the predictive performance (concordance indexes) of clinicopathological factors. Transcriptomic analyzes linked our signature with important tumor biological processes (e.g. glycolysis/glucoseogenesis, ribosome). Then, we identified 36 essential resistant genes, and constructed a gene-agent network including 10 essential resistant genes and 35 candidate therapeutic agents, and excavated AT-7519 as the therapeutic agent for patients with CCR resistance. The therapeutic efficacy of AT-7519 was validated that significantly more resistant genes were down-regulated induced by AT-7519, and the degree gradually increased with the enhanced doses.

Conclusions

This study illustrated that radiomic signature could non-invasively predict therapeutic efficacy of patients with NSCLC receiving CCR, and indicated that patients with CCR resistance might benefit from AT-7519 or CCR treatment combined with AT-7519.

A Multi-Omics Study on the Effect of Helicobacter Pylori-Related Genes in the Tumor Immunity on Stomach Adenocarcinoma

Background

Helicobacter pylori (HP), a gram-negative spiral-shaped microaerophilic bacterium, colonizes the stomach of approximately 50% of the world’s population, which is considered a risk factor for gastritis, peptic ulcers, gastric cancer, and other malignancies. HP is also considered carcinogenic since it involves the mutation and damage of multiple HP-related genes. Stomach adenocarcinoma (STAD) is a common stom5ach cancer with a poor prognosis and high risk of metastasis in the advanced stage. Therefore, an early diagnosis and targeted therapies are needed to ensure a better prognosis. In this study, a scoring system was constructed based on three HP infection–related candidate genes to enable a more accurate prediction of tumor progression and metastasis and response to immunotherapies.

Methods

HP infection–induced mutation patterns of STAD samples from six cohorts were comprehensively assessed based on 73 HP-related genes, which were then correlated with the immune cell–infiltrating characteristics of the tumor microenvironment (TME). The risk signature was constructed to quantify the influence of HP infection on individual tumors. Subsequently, an accurate nomogram was generated to improve the clinical applicability of the risk signature. We conducted immunohistochemical experiments and used the Affiliated Hospital of Youjiang Medical University for Nationalities (AHYMUN) cohort data set with survival information to further verify the clinical value of this risk signature.

Results

Two distinct HP-related mutation patterns with different immune cell–infiltrating characteristics (ICIC) and survival possibility were identified. We demonstrated that the evaluation of HP infection–induced mutation patterns of tumor could assist the prediction of stages, phenotypes, stromal activity, genetic diversity, and patient prognosis. A low risk score involved an increased mutation burden and activation of immune responses, with a higher 5-year survival rate and enhanced response to anti-PD-1/L1 immunotherapy, while a high risk score involved stromal activation and poorer survival. The efficiency of the risk signature was further evidenced by the nomogram.

Conclusions

STAD patients with a low risk score demonstrated significant therapeutic advantages and clinical benefits. HP infection–induced mutations play a nonnegligible role in STAD development. Quantifying the HP-related mutation patterns of individual tumors will contribute to phenotype classification, guide more effective targeted and personalized therapies, and enable more accurate predictions of metastasis and prognosis.

Risk stratification of lung adenocarcinoma using a nomogram combined with ferroptosis-related LncRNAs and subgroup analysis with immune and N6-methyladenosine modification

Background

Determining the prognosis of lung adenocarcinoma (LUAD) is challenging. The present study aimed to identify prognostic ferroptosis-related long noncoding RNAs (FRLs) and construct a prognostic model. Moreover, differential analysis of immune and N6-methyladenosine (m6A)-related genes was systematically conducted.

Methods

A total of 504 patients selected from a dataset from The Cancer Genome Atlas were included. The patients with LUAD were randomly divided into a training group and a test group at a ratio of 1:1. Pearson correlation analysis and univariate Cox regression analysis were used to identify the prognostic FRLs. Then, a prognostic model was constructed from the optimized subset of prognostic FRLs based on the least absolute shrinkage and selection operator (LASSO) algorithm. Subsequently, the receiver operating characteristic (ROC) curve and survival analysis were used to evaluate the performance of the model. The risk score based on the prognostic model was analyzed using Cox regression analysis. Moreover, gene set enrichment analysis and differential analysis of immune- and m6A-related genes were conducted.

Results

After univariate Cox regression analysis and LASSO algorithm analysis, a total of 19 prognostic FRLs were selected to construct the final model to obtain the risk score. The area under the ROC curve of the prognostic model for 1-year, 3-year, and 5-year overall survival (OS) was 0.763, 0.745, and 0.778 in the training set and 0.716, 0.724, and 0.736 in the validation set, respectively. Moreover, the OS of the high-risk group was significantly worse than that of the low-risk group in the training group (P < 0.001) and in the test group (P < 0.001). After univariate and multivariate Cox regression analysis, the risk score [hazard ratio (HR) = 1.734; P < 0.001] and stage (HR = 1.557; P < 0.001) were both considered significant prognostic factors for LUAD. A nomogram was constructed based on clinical features and risk score. The expression of 34 checkpoint genes and 13 m6A-related genes varied significantly between the two risk groups.

Conclusion

This study constructed a prognostic model to effectively predict the OS of patients with LUAD, and these OS-related FRLs might serve as potential therapeutic targets of LUAD.

Gene Expression Profiles of Multiple Synchronous Lesions in Lung Adenocarcinoma

Many studies support a stepwise continuum of morphologic changes between atypical adenomatous hyperplasia (AAH) and lung adenocarcinoma (ADC). Here we characterized gene expression patterns and the association of differentially expressed genes and immune tumor microenvironment behaviors in AAH to ADC during ADC development. Tumor tissues from nine patients with ADC and synchronous multiple ground glass nodules/lesions (GGN/Ls) were analyzed using RNA sequencing. Using clustering, we identified genes differentially and sequentially expressed in AAH and ADC compared to normal tissues. Functional enrichment analysis using gene ontology terms was performed, and the fraction of immune cell types was estimated. We identified up-regulated genes (ACSL5 and SERINC2) with a stepwise change of expression from AAH to ADC and validated those expressions by quantitative PCR and immunohistochemistry. The immune cell profiles revealed increased B cell activities and decreased natural killer cell activities in AAH and ADC. A stepwise change of differential expression during ADC development revealed potential effects on immune function in synchronous precursors and in tumor lesions in patients with lung cancer.

Gene signatures with therapeutic value: emerging perspective for personalized immunotherapy in renal cancer

Renal cancer is one of the deadliest urogenital diseases. In recent years, the advent of immunotherapy has led to significant improvement in the management of patients with renal cancer. Although cancer immunotherapy and its combinations had benefited numerous patients, several challenges need to be addressed. Apart from the high costs of treatment, the lack of predictive biomarkers and toxic side-effects have impeded its wider applicability. To address these issues, new biomarkers are required to predict responsiveness and design personalized treatment strategies. Recent advances in the field of single-cell sequencing and multi-dimensional spatial transcriptomics have identified clinically relevant subtypes of renal cancer. Furthermore, there is emerging potential for gene signatures based on immune cells, non-coding RNAs, and pathways such as metabolism and RNA modification. In this review article, we have discussed recent progress in the identification of gene signatures with predictive and prognostic potential in renal cancer.

A DNA methylation profile of long non-coding RNAs can predict OS in prostate cancer

Prostate cancer (PCa) is the most common male reproductive tract malignant tumor, accurate evaluation of PCa characterization and prognostic prediction at diagnosis are vital for the effective administration of the disease, especially at the molecular level. In this study, 48 CpG sites with differential methylation associated with overall survival (OS) were screened out between PCa and normal adjacent tissues. 16 CpG sites were selected by the least absolute shrinkage and selection operator (LASSO) and the risk score formula for methylated-based classifier was established. For 16-lncRNAs-CpG-classifier, the area under the curve (AUC) were 0.890, 0.917, and 0.932 at 3 years, 5 years and 7 years, respectively. Kaplan–Meier curves indicated that patients with high-risk scores had worse OS than those with low-risk scores. Prognostic methylation model of lncRNAs was identified from the whole genome in patients with PCa. This novel finding provides a novel insight for screening biomarkers of a prognosis for PCa.

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.