Identification of a Novel Prognostic Signature of Genome Instability-Related LncRNAs in Early Stage Lung Adenocarcinoma

Identification of Genomic Instability-Associated LncRNAs as Potential Therapeutic Targets in Lung Adenocarcinoma

BACKGROUND/OBJECTIVES

Non-small cell lung cancer (NSCLC) is the most common type of cancer, with lung adenocarcinoma (LUAD) as the predominant subtype. Despite advancements in targeted therapies, many NSCLC patients still experience poor outcomes due to treatment resistance and disease progression. Genomic instability (GI), a hallmark of cancer, defined as the increased tendency of DNA mutations and alterations, is closely linked to cancer initiation, progression, and resistance to therapy. Emerging evidence suggests that long non-coding RNAs (lncRNAs)-molecules longer than 200 nucleotides that do not encode proteins but regulate gene expression-play critical roles in cancer biology and are associated with GI. However, the relationship between GI and lncRNA expression in LUAD remains poorly understood.

METHODS

In this study, we analyzed the transcript profiles of lncRNAs and mRNAs from LUAD samples in The Cancer Genome Atlas (TCGA) database and classified them based on their Homologous Recombination Deficiency (HRD) score. The HRD score is an unweighted sum of three independent DNA-based measures of genomic instability: loss of heterozygosity, telomeric allelic imbalance, and large-scale transitions. We then performed a differential gene expression analysis to identify lncRNAs and mRNAs that were either upregulated or downregulated in samples with high HRD scores compared to those with low HRD scores. Following this, we conducted a correlation analysis to assess the significance of the association between HRD scores and the expression of both lncRNAs and mRNAs.

RESULTS

We identified 30 differentially expressed lncRNAs and 200 mRNAs associated with genomic instability. Using an RNA interactome database from sequencing experiments, we found evidence of interactions between GI-associated lncRNAs (GI-lncRNAs) and GI-associated mRNAs (GI-mRNAs). Further investigation showed that some GI-lncRNAs play regulatory and functional roles in LUAD and other diseases. We also found that GI-lncRNAs have potential as prognostic biomarkers, particularly when integrated with HRD stratification. The expression of specific GI-lncRNAs was associated with primary therapy response and immune infiltration in LUAD. Additionally, we identified existing drugs that could modulate GI-lncRNAs, offering potential therapeutic strategies to address GI in LUAD.

CONCLUSIONS

Our findings suggest that GI-associated lncRNAs could serve as valuable biomarkers for LUAD prognosis and therapeutic response. Furthermore, modulating these lncRNAs presents potential treatment avenues to address genomic instability in LUAD.

Deciphering disulfidptosis: Uncovering a lncRNA-based signature for prognostic assessment, personalized immunotherapy, and therapeutic agent selection in lung adenocarcinoma patients

BACKGROUND

Disulfidptosis, a recently identified type of regulated cell death, plays critical roles in various biological processes of cancer; however, whether they can impact the prognosis of lung adenocarcinoma (LUAD) remains to be fully elucidated. We aimed to adopt this concept to develop and validate a lncRNA signature for LUAD prognostic prediction.

METHODS

For this study, the TCGA-LUAD dataset was used as the training cohort, and multiple datasets from the GEO database were pooled as the validation cohort. Disulfidptosis regulated genes were obtained from published studies, and various statistical methods, including Kaplan-Meier (KM), Cox, and LASSO, were used to train our gene signature DISULncSig. We utilized KM analysis, COX analysis, receiver operating characteristic analysis, time-dependent AUC analysis, principal component analysis, nomogram predictor analysis, and functional assays in our validation process. We also compared DISULncSig with previous studies. We performed analyses to evaluate DISULncSig's immunotherapeutic ability, focusing on eight immune algorithms, TMB, and TIDE. Additionally, we investigated potential drugs that could be effective in treating patients with high-risk scores. Additionally qRT-PCR examined the expression patterns of DISULncSig lncRNAs, and the ability of DISULncSig in pan-cancer was also assessed.

RESULTS

DISULncSig containing twelve lncRNAs was trained and showed strong predictive ability in the validation cohort. Compared with previous similar studies, DISULncSig had more prognostic ability advantages. DISULncSig was closely related to the immune status of LUAD, and its tight relationship with checkpoints KIR2DL3, IL10, IL2, CD40LG, SELP, BTLA, and CD28 may be the key to its potential immunotherapeutic ability. For the high DISULncSig score population, we found ten drug candidates, among which epothilone-b may have the most potential. The pan-cancer analysis found that DISULncSig was a risk factor in multiple cancers. Additionally, we discovered that some of the DISULncSig lncRNAs could play crucial roles in specific cancer types.

CONCLUSION

The current study established a powerful prognostic DISULncSig signature for LUAD that was also valid for most pan-cancers. This signature could serve as a potential target for immunotherapy and might help the more efficient application of drugs to specific populations.

Comprehensive landscape and future perspective of long noncoding RNAs in non-small cell lung cancer: it takes a village

Long noncoding RNAs (lncRNAs) are a distinct subtype of RNA that lack protein-coding capacity but exert significant influence on various cellular processes. In non-small cell lung cancer (NSCLC), dysregulated lncRNAs act as either oncogenes or tumor suppressors, contributing to tumorigenesis and tumor progression. LncRNAs directly modulate gene expression, act as competitive endogenous RNAs by interacting with microRNAs or proteins, and associate with RNA binding proteins. Moreover, lncRNAs can reshape the tumor immune microenvironment and influence cellular metabolism, cancer cell stemness, and angiogenesis by engaging various signaling pathways. Notably, lncRNAs have shown great potential as diagnostic or prognostic biomarkers in liquid biopsies and therapeutic strategies for NSCLC. This comprehensive review elucidates the significant roles and diverse mechanisms of lncRNAs in NSCLC. Furthermore, we provide insights into the clinical relevance, current research progress, limitations, innovative research approaches, and future perspectives for targeting lncRNAs in NSCLC. By summarizing the existing knowledge and advancements, we aim to enhance the understanding of the pivotal roles played by lncRNAs in NSCLC and stimulate further research in this field. Ultimately, unraveling the complex network of lncRNA-mediated regulatory mechanisms in NSCLC could potentially lead to the development of novel diagnostic tools and therapeutic strategies.

Genome instability-derived genes as a novel prognostic signature for lung adenocarcinoma

Background: An increasing number of patients are being diagnosed with lung adenocarcinoma, but there remains limited progress in enhancing prognostic outcomes and improving survival rates for these patients. Genome instability is considered a contributing factor, as it enables other hallmarks of cancer to acquire functional capabilities, thus allowing cancer cells to survive, proliferate, and disseminate. Despite the importance of genome instability in cancer development, few studies have explored the prognostic signature associated with genome instability for lung adenocarcinoma. Methods: In the study, we randomly divided 397 lung adenocarcinoma patients from The Cancer Genome Atlas database into a training group (n = 199) and a testing group (n = 198). By calculating the cumulative counts of genomic alterations for each patient in the training group, we distinguished the top 25% and bottom 25% of patients. We then compared their gene expressions to identify genome instability-related genes. Next, we used univariate and multivariate Cox regression analyses to identify the prognostic signature. We also performed the Kaplan-Meier survival analysis and the log-rank test to evaluate the performance of the identified prognostic signature. The performance of the signature was further validated in the testing group, in The Cancer Genome Atlas dataset, and in external datasets. We also conducted a time-dependent receiver operating characteristic analysis to compare our signature with established prognostic signatures to demonstrate its potential clinical value. Results: We identified GULPsig, which includes IGF2BP1, IGF2BP3, SMC1B, CLDN6, and LY6K, as a prognostic signature for lung adenocarcinoma patients from 42 genome instability-related genes. Based on the risk score of the risk model with GULPsig, we successfully stratified the patients into high- and low-risk groups according to the results of the Kaplan-Meier survival analysis and the log-rank test. We further validated the performance of GULPsig as an independent prognostic signature and observed that it outperformed established prognostic signatures. Conclusion: We provided new insights to explore the clinical application of genome instability and identified GULPsig as a potential prognostic signature for lung adenocarcinoma patients.

Application of an angiogenesis-related genes risk model in lung adenocarcinoma prognosis and immunotherapy

Lung adenocarcinoma (LUAD) is an essential pathological subtype of non-small cell lung cancer and offers a severe problem for worldwide public health. There is mounting proof that angiogenesis is a crucial player in LUAD progression. Consequently, the purpose of this research was to construct a novel LUAD risk assessment model based on genetic markers related to angiogenesis. We accessed The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases for LUAD mRNA sequencing data and clinical information. Based on machine algorithms and bioinformatics, angiogenic gene-related risk scores (RS) were calculated. Patients in the high-risk category had a worse prognosis (p < 0.001) in the discovery TCGA cohort, and the results were confirmed by these three cohorts (validation TCGA cohort, total TCGA cohort, and GSE68465 cohort). Moreover, risk scores for genes involved in angiogenesis were independent risk factors for lung cancer in all four cohorts. The low-risk group was associated with better immune status and lower tumor mutational load. In addition, the somatic mutation study revealed that the low-risk group had a lower mutation frequency than the high-risk group. According to an analysis of tumor stem cell infiltration, HLA expression, and TIDE scores, the low-risk group had higher TIDE scores and HLA expression levels than the high-risk group, and the amount of tumor stem cell infiltration correlated with the risk score. In addition, high-risk groups may benefit from immune checkpoint inhibitors and targeted therapies. In conclusion, we developed an angiogenesis-related gene risk model to predict the prognosis of LUAD patients, which may aid in the classification of patients with LUAD and select medications for LUAD patients.

Novel GIRlncRNA Signature for Predicting the Clinical Outcome and Therapeutic Response in NSCLC

Background: Non–small cell lung cancer (NSCLC) is highly malignant with driver somatic mutations and genomic instability. Long non-coding RNAs (lncRNAs) play a vital role in regulating these two aspects. However, the identification of somatic mutation-derived, genomic instability-related lncRNAs (GIRlncRNAs) and their clinical significance in NSCLC remains largely unexplored.

Methods: Clinical information, gene mutation, and lncRNA expression data were extracted from TCGA database. GIRlncRNAs were screened by a mutator hypothesis-derived computational frame. Co-expression, GO, and KEGG enrichment analyses were performed to investigate the biological functions. Cox and LASSO regression analyses were performed to create a prognostic risk model based on the GIRlncRNA signature (GIRlncSig). The prediction efficiency of the model was evaluated by using correlation analyses with mutation, driver gene, immune microenvironment contexture, and therapeutic response. The prognostic performance of the model was evaluated by external datasets. A nomogram was established and validated in the testing set and TCGA dataset.

Results: A total of 1446 GIRlncRNAs were selected from the screen, and the established GIRlncSig was used to classify patients into high- and low-risk groups. Enrichment analyses showed that GIRlncRNAs were mainly associated with nucleic acid metabolism and DNA damage repair pathways. Cox analyses further identified 19 GIRlncRNAs to construct a GIRlncSig-based risk score model. According to Cox regression and stratification analyses, 14 risk lncRNAs (AC023824.3, AC013287.1, AP000829.1, LINC01611, AC097451.1, AC025419.1, AC079949.2, LINC01600, AC004862.1, AC021594.1, MYRF-AS1, LINC02434, LINC02412, and LINC00337) and five protective lncRNAs (LINC01067, AC012645.1, AL512604.3, AC008278.2, and AC089998.1) were considered powerful predictors. Analyses of the model showed that these GIRlncRNAs were correlated with somatic mutation pattern, immune microenvironment infiltration, immunotherapeutic response, drug sensitivity, and survival of NSCLC patients. The GIRlncSig risk score model demonstrated good predictive performance (AUCs of ROC for 10-year survival was 0.69) and prognostic value in different NSCLC datasets. The nomogram comprising GIRlncSig and tumor stage exhibited improved robustness and feasibility for predicting NSCLC prognosis.

Conclusion: The newly identified GIRlncRNAs are powerful biomarkers for clinical outcome and prognosis of NSCLC. Our study highlights that the GIRlncSig-based score model may be a useful tool for risk stratification and management of NSCLC patients, which deserves further evaluation in future prospective studies.

Mitochondrial Homeostasis–Related lncRNAs are Potential Biomarkers for Predicting Prognosis and Immune Response in Lung Adenocarcinoma

The prognosis of the most common histological subtype of lung cancer, lung adenocarcinoma (LUAD), is relatively poor. Mitochondrial homeostasis depends to a great extent on the coordination between mitophagy and mitochondrial biogenesis, the deregulation of which causes various human diseases, including cancer. There is accumulating evidence that long noncoding RNAs (lncRNAs) are critical in predicting the prognosis and immune response in carcinoma. Therefore, it is critical to discern lncRNAs related to mitochondrial homeostasis in LUAD patients. In this study, we identified mitochondrial homeostasis–related lncRNAs (MHRlncRNAs) by coexpression analysis. In order to construct a prognostic signature composed of three MHRlncRNAs, univariate and multivariate Cox regression analyses were performed. Kaplan–Meier analysis, stratification analysis, principal component analysis (PCA), receiver operating characteristic (ROC) curve, gene set enrichment analysis (GSEA), and nomogram were applied to evaluate and optimize the risk model. Subsequently, we identified the mitochondrial homeostasis–related lncRNA signature (MHLncSig) as an independent predictive factor of prognosis. Based on the LUAD subtypes regrouped by this risk model, we further investigated the underlying tumor microenvironment, tumor mutation burden, and immune landscape behind different risk groups. Likewise, individualized immunotherapeutic strategies and candidate compounds were screened to aim at different risk subtypes of LUAD patients. Finally, we validated the expression trends of lncRNAs included in the risk model using quantitative real-time polymerase chain reaction (qRT-PCR) assays. The established MHLncSig may be a promising tool for predicting the prognosis and guiding individualized treatment in LUAD.

Comprehensive Analyses of Mutation-Derived Long-Chain Noncoding RNA Signatures of Genome Instability in Kidney Renal Papillary Cell Carcinoma

Background: The role of long-chain noncoding RNA (lncRNA) in genomic instability has been demonstrated to be increasingly importance. Therefore, in this study, lncRNAs associated with genomic instability were identified and kidney renal papillary cell carcinoma (KIRP)-associated predictive features were analysed to classify high-risk patients and improve individualised treatment.

Methods: The training (n = 142) and test (n = 144) sets were created using raw RNA-seq and patient’s clinical data of KIRP obtained from The Cancer Genome Atlas (TCGA).There are 27 long-chain noncoding RNAs (lncRNAs) that are connected with genomic instability, these lncRNAs were identified using the ‘limma’ R package based on the numbers of somatic mutations and lncRNA expression profiles acquired from KIRP TCGA cohort. Furthermore, Cox regression analysis was carried out to develop a genome instability-derived lncRNA-based gene signature (GILncSig), whose prognostic value was confirmed in the test cohort as well as across the entire KIRP TCGA dataset.

Results: A GILncSig derived from three lncRNAs (BOLA3-AS1, AC004870, and LINC00839), which were related with poor KIRP survival, was identified, which was split up into high- and low-risk groups. Additionally, the GILncSig was found to be an independent prognostic predictive index in KIRP using univariate and multivariate Cox analysis. Furthermore, the prognostic significance and characteristics of GilncSig were confirmed in the training test and TCGA sets. GilncSig also showed better predictive performance than other prognostic lncRNA features.

Conclusion: The function of lncRNAs in genomic instability and the genetic diversity of KIRP were elucidated in this work. Moreover, three lncRNAs were screened for prediction of the outcome of KIRP survival and novel insights into identifying cancer biomarkers related to genomic instability were discussed.

Diagnostic and prognostic value of genomic instability-derived long non-coding RNA signature of endometrial cancer

OBJECTIVE

To investigate whether genomic instability (GI)-derived long non-coding RNAs (lncRNAs) have a prognostic impact on the patients with endometrial cancer.

MATERIAL AND METHODS

Patients with Uterine Corpus Endometrial Carcinoma (UCEC) were selected from The Cancer Genome Atlas (TCGA) database. Systematic bioinformatics analyses were performed, including Pearson correlations, GO and KEGG enrichment analysis, bivariate and multiple logistic regression analysis, and Kaplan-Meier (KM) method.

RESULTS

A total of 552 UCEC samples were included in the study. The differentially expressed lncRNAs (DELs) were identified, including 79 down-regulated lncRNAs and 31 up-regulated lncRNAs. Bivariate logistic regression analysis showed that 19 GI-derived lncRNAs were prognostic factors. By further multivariate logistic regression analysis, AC005256.1 (estimated coefficient = -0.474), AC026336.3 (estimated coefficient = -0.030), AL161618.1 (estimated coefficient = -1.661), and BX322234.1 (estimated coefficient = 1.511) were used to construct a prognostic risk model. In the train set and test set, the risk model was shown to have both a high prognostic and a diagnostic value.

CONCLUSION

We developed a novel GI-derived 4-lncRNA signature for the diagnosis and prognosis of patients with endometrial cancer. These findings offered a novel perspective in the clinical management of endometrial cancer.

Prognosis and Immunotherapy Response With a Novel Golgi Apparatus Signature-Based Formula in Lung Adenocarcinoma

The Golgi apparatus (GA) is a cellular organelle that participates in the packaging, modification, and transport of proteins and lipids from the endoplasmic reticulum to be further fabricated before being presented to other cellular components. Recent studies have demonstrated that GA facilitates numerous cellular processes in cancer development. Therefore, this study aimed to establish a novel lung adenocarcinoma (LUAD) risk evaluation model based on GA gene signatures. In this study, we used TCGA-LUAD (n = 500) as the training cohort and GSE50081 (n = 127), GSE68465 (442), and GSE72094 (398) as the validation cohorts. Two immunotherapy datasets (GSE135222 and GSE126044) were also obtained from a previous study. Based on machine algorithms and bioinformatics methods, a GA gene-related risk score (GARS) was established. We found that the GARS independently predicted the prognosis of LUAD patients and remained effective across stages IA to IIIA. Then, we identified that the GARS was highly correlated with mutations in P53 and TTN. Further, this study identified that GARS is related to multiple immune microenvironmental characteristics. Furthermore, we investigated GSE135222 and GSE126044 and found that a lower GARS may be indicative of an improved therapeutic effect of PD-1/PD-L1 therapy. We also found that high GARS may lead to a better response to multiple anticancer drugs. Finally, we established a nomogram to better guide clinical application. To our knowledge, this is the first study to demonstrate a novel GA signature-based risk score formula to predict clinical prognosis and guide the treatment of LUAD patients.

Comprehensive Analysis of Cell Cycle-Related Genes in Patients With Prostate Cancer

This study aimed to identify critical cell cycle-related genes (CCRGs) in prostate cancer (PRAD) and to evaluate the clinical prognostic value of the gene panel selected. Gene set variation analysis (GSVA) of dysregulated genes between PRAD and normal tissues demonstrated that the cell cycle-related pathways played vital roles in PRAD. Patients were classified into four clusters, which were associated with recurrence-free survival (RFS). Moreover, 200 prognostic-related genes were selected using the Kaplan-Meier (KM) survival analysis and univariable Cox regression. The prognostic CCRG risk score was constructed using random forest survival and multivariate regression Cox methods, and their efficiency was validated in Memorial Sloan Kettering Cancer Center (MSKCC) and GSE70770. We identified nine survival-related genes: CCNL2, CDCA5, KAT2A, CHTF18, SPC24, EME2, CDK5RAP3, CDC20, and PTTG1. Based on the median risk score, the patients were divided into two groups. Then the functional enrichment analyses, mutational profiles, immune components, estimated half-maximal inhibitory concentration (IC50), and candidate drugs were screened of these two groups. In addition, the characteristics of nine hub CCRGs were explored in Oncomine, cBioPortal, and the Human Protein Atlas (HPA) datasets. Finally, the expression profiles of these hub CCRGs were validated in RWPE-1 and three PRAD cell lines (PC-3, C4-2, and DU-145). In conclusion, our study systematically explored the role of CCRGs in PRAD and constructed a risk model that can predict the clinical prognosis and immunotherapeutic benefits.

Expression Landscape and Functional Roles of HOXA4 and HOXA5 in Lung Adenocarcinoma

BACKGROUND

The role of HOXA family genes in the occurrence and progression of a variety of human cancers has been scatteredly reported. However, there is no systematic study on the differential expression, prognostic significance and potential molecular mechanism of HOXA4 and HOXA5 in LUAD.

METHODS

In-house immunohistochemistry (IHC), multi-center microarrays, RT-qPCR and RNA-seq data were incorporated for comprehensively evaluating the expression and prognostic value of HOXA4 and HOXA5 in LUAD. The mechanism of HOXA4 and HOXA5 in the formation and development of LUAD was analyzed from multiple aspects of immune correlations, upstream transcriptional regulation, functional states of single cells and co-expressed gene network. The functional roles of HOXA4 and HOXA5 in LUAD were validated by in vitro experiments.

RESULTS

As a result, in 3201 LUAD samples and 2494 non-cancer lung samples, HOXA4 and HOXA5 were significantly downexpressed (P < 0.05). The aberrant expression of HOXA5 was significantly correlated with the clinical progression of LUAD (P < 0.05). HOXA5 showed remarkable prognostic value for LUAD patients (P < 0.05). The expression of HOXA4 and HOXA5 in LUAD were negatively correlated with tumor purity and positively correlated with the infiltration of various immune cells such as B cells, T cells and macrophages. HOXA4 and HOXA5 overexpression had notable inhibitory effect on the proliferation, migration and invasion of LUAD cells.

CONCLUSIONS

In conclusion, the identified downexpressed HOXA4 and HOXA5 had significant distinguishing ability for LUAD samples and affected the cellular functions of LUAD cells. The low expression of HOXA5 indicated worse overall survival of LUAD patients. Therefore, the two HOXA family genes especially HOXA5 may serve as potential biomarkers for LUAD.