An autophagy-related long non-coding RNA prognostic signature accurately predicts survival outcomes in bladder urothelial carcinoma patients

Cuproptosis-Related lncRNAs Modulate the Prognosis of MIBC by Regulating the Expression Pattern of Immunosuppressive Molecules Within the Tumor Microenvironment

Background

Cuproptosis-related gene and long non-coding RNA (lncRNA) modulation of cancer regulation is well-established. This investigation aimed to elucidate the prognostic implications of cuproptosis-associated lncRNAs in muscle-invasive bladder cancer (MIBC).

Methods

Employing the Cancer Genome Atlas (TCGA) and IMvigor210 cohorts, bioinformatics and statistical analyses probed the prognostic relevance of cuproptosis-related lncRNAs.

Results

Co-expression analysis revealed tight associations between lncRNA expression and cuproptosis-linked genes, with 13 cuproptosis-related lncRNAs found to correlate with MIBC prognosis. Lasso regression identified a six-lncRNA prognostic signature, enabling patient stratification into high- and low-risk categories. Tissue validation substantiated differential expression of FAM13A-AS1, GHRLOS, LINC00456, OPA1-AS1, RAP2C-AS1, and UBE2Q1-AS1 between MIBC tumor and normal tissues. Comparative analyses of tumor microenvironments and immune profiles between risk groups disclosed elevated immunosuppressive molecule expression, including programmed cell death-1 (PD-L1) and T-cell immunoglobulin-3 (TIM-3), in high-risk individuals.

Conclusion

These findings suggest that cuproptosis-related lncRNAs may modulate the expression of immunosuppressive molecules, thereby influencing MIBC tumorigenesis and progression. Further exploration is warranted to unveil novel therapeutic targets for MIBC based on the expression patterns of cuproptosis-related lncRNAs and their impact on immune responses in the tumor microenvironment.

The development and experimental validation of hypoxia-related long noncoding RNAs prognostic signature in predicting prognosis and immunotherapy of cutaneous melanoma

Cutaneous melanoma (CM) is widely acknowledged as a highly aggressive form of malignancy that is associated with a considerable degree of morbidity and poor prognosis. Despite this recognition, the precise role of hypoxia-related long noncoding RNAs (HRLs) in the pathogenesis of CM remains an area of active research. This study sought to elucidate the contribution of HRLs in CM by conducting a thorough screening and extraction of hypoxia-related genes (HRGs). In particular, we conducted univariate and multivariate Cox regression analyses to assess the independence of the prognostic signature of HRLs. Our results demonstrated that a novel risk model could be established based on five prognostic HRLs. Remarkably, patients with low-risk scores exhibited significantly higher overall survival rates compared to their high-risk counterparts, as confirmed by Kaplan-Meier survival analysis. Furthermore, we utilized consensus clustering analysis to categorize CM patients into two distinct subtypes, which revealed marked differences in their prognosis and immune infiltration landscapes. Our nomogram results confirmed that the HRLs prognostic signature served as an independent prognostic indicator, offering an accurate evaluation of the survival probability of CM patients. Notably, our findings from ESTIMATE and ssGSEA analyses highlighted significant disparities in the immune infiltration landscape between low- and high-risk groups of CM patients. Additionally, IPS and TIDE results suggested that CM patients in different risk subtypes may exhibit favorable responses to immunotherapy. Enrichment analysis and GSVA results indicated that immune-related signaling pathways may mediate the role of HRLs in CM. Finally, our tumor mutation burden (TMB) results indicated that patients with low-risk scores had a higher TMB status. In summary, the establishment of a risk model based on HRLs in this study provided an accurate prognostic prediction and correlated with the immune infiltration landscape of CM, thereby providing novel insights for the future clinical management of this disease.

Unravelling the complexity of lncRNAs in autophagy to improve potential cancer therapy

Autophagy is well-known as an internal catabolic process that is evolutionarily conserved and performs the key biological function in maintaining cellular homeostasis. It is tightly controlled by several autophagy-related (ATG) proteins, which are closely associated with many types of human cancers. However, what has remained controversial is the janus roles of autophagy in cancer progression. Interestingly, the biological function of long non-coding RNAs (lncRNAs) in autophagy has been gradually understood in different types of human cancers. More recently, numerous studies have demonstrated that several lncRNAs may regulate some ATG proteins and autophagy-related signaling pathways to either activate or inhibit the autophagic process in cancer. Thus, in this review, we summarize the latest advance in the knowledge of the complicated relationships between lncRNAs and autophagy in cancer. Also, the in-depth dissection of the lncRNAs-autophagy-cancers axis involved in this review would shed new light on discovery of more potential cancer biomarkers and therapeutic targets in the future.

New algorithms based on autophagy-related lncRNAs pairs to predict the prognosis of skin cutaneous melanoma patients

Skin cutaneous melanoma (SKCM) is the most malignant skin tumor for it is enormously easy to develop invasion and metastasis. Autophagy is a process by which cellular material is degraded by lysosomes or vacuoles and recycled. Autophagy-related long non-coding RNAs (lncRNAs) have been thought to correlate with SKCM. This study aims to explore the prognostic significance of autophagy-related lncRNAs and establish a prognostic model of autophagy-related lncRNA pairs in SKCM. Firstly, the RNA-seq data and related clinical information were downloaded from the TCGA database. 446 qualified samples were enrolled. 222 autophagy-related genes were obtained from the HADb database. Pearson correlation analysis was conducted to identify autophagy-related lncRNAs (ARLs). After that, we obtained prognosis-related ARLs and autophagy-related lncRNA pairs (ARLPs). Using Lasso-Cox regression analysis, an autophagy-related lncRNA-pair prognostic signature was established. The accuracy of the signature were confirmed through a series of validations in terms of mutation profiles, immunity infiltration, and cellular pathways. And we used the random forest method to find USP30-AS1 as a key mediating factor in SKCM.

Comprehensive analysis of cuproptosis-related lncRNAs signature to predict prognosis in bladder urothelial carcinoma

BACKGROUND

Cuproptosis-related genes (CRGs) have been recently discovered to regulate the occurrence and development of various tumors by controlling cuproptosis, a novel type of copper ion-dependent cell death. Although cuproptosis is mediated by lipoylated tricarboxylic acid cycle proteins, the relationship between cuproptosis-related long noncoding RNAs (crlncRNAs) in bladder urothelial carcinoma (BLCA) and clinical outcomes, tumor microenvironment (TME) modification, and immunotherapy remains unknown. In this paper, we tried to discover the importance of lncRNAs for BLCA.

METHODS

The BLCA-related lncRNAs and clinical data were first obtained from The Cancer Genome Atlas (TCGA). CRGs were obtained through Coexpression, Cox regression and Lasso regression. Besides, a prognosis model was established for verification. Meanwhile, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, gene ontology (GO) analysis, principal component analysis (PCA), half-maximal inhibitory concentration prediction (IC50), immune status and drug susceptibility analysis were carried out.

RESULTS

We identified 277 crlncRNAs and 16 survival-related lncRNAs. According to the 8-crlncRNA risk model, patients could be divided into high-risk group and low-risk group. Progression-Free-Survival (PFS), independent prognostic analysis, concordance index (C-index), receiver operating characteristic (ROC) curve and nomogram all confirmed the excellent predictive capability of the 8-lncRNA risk model for BLCA. During gene mutation burden survival analysis, noticeable differences were observed in high- and low-risk patients. We also found that the two groups of patients might respond differently to immune targets and anti-tumor drugs.

CONCLUSION

The nomogram with 8-lncRNA may help guide treatment of BLCA. More clinical studies are necessary to verify the nomogram.

Crosstalk between autophagy and bladder transitional cell carcinoma by autophagy-related lncRNAs

The aim of this study was to investigate the crosstalk between autophagy and bladder transitional cell carcinoma (TCC) by autophagy-related long noncoding RNAs (lncRNAs). A total of 400 TCC patients from The Cancer Genome Atlas were enrolled in this study. We identified the autophagy-related lncRNA expression profile of the TCC patients and then constructed a prognostic signature using the least absolute shrinkage and selection operation and Cox regression. Risk, survival, and independent prognostic analyses were carried out. Receiver operating characteristic curve, nomogram, and calibration curves were explored. Gene Set Enrichment Analysis was employed to verify the enhanced autophagy-related functions. Finally, we compared the signature with several other lncRNA-based signatures. A 9-autophagy-related lncRNA signature was established by least absolute shrinkage and selection operation-Cox regression that was significantly associated with overall survival in TCC. Among them, 8 of the 9 lncRNAs were protective factors while the remaining was a risk factor. The risk scores calculated by the signature showed significant prognostic value in survival analysis between the high- or low-risk groups. The 5-year survival rate for the high-risk group was 26.0% while the rate for the low-risk group was 56.0% (P < .05). Risk score was the only significant risk factor in the multivariate Cox regression survival analysis (P < .001). A nomogram connecting this signature with clinicopathologic characteristics was assembled. To assess the performance of the nomogram, a C-index (0.71) was calculated, which showed great convergence with an ideal model. The Gene Set Enrichment Analysis results demonstrated 2 major autophagy-related pathways were significantly enhanced in TCC. And this signature performed a similar predictive effect as other publications. The crosstalk between autophagy and TCC is significant, and this 9 autophagy-related lncRNA signature is a great predictor of TCC.

Construction and experimental validation of a B cell-related gene signature to predict the prognosis and immunotherapeutic sensitivity in bladder cancer

BACKGROUND

B cells are essential components of tumor microenvironment and exert important functions in anti-tumor immune response. However, the prognosis value of B cell-related genes in bladder cancer (BLCA) remains obscure.

MATERIALS AND METHODS

The infiltrating levels of B cells were measured via the CD20 staining in the local samples and the computational biology analyses in the TCGA-BLCA cohort. The single-cell RNA sequencing analysis, gene-pair strategy, LASSO regression, random forest, and Cox regression were used for B cell-related signature construction. TCGA-BLCA cohort was chosen as the training cohort, and three independent cohorts from GEO and the local cohort were used for external validation. 326 B cells were adopted to explore the association between the model and B cells' biological processes. TIDE algorithm and two BLCA cohorts receiving anti-PD1/PDL1 treatment were utilized to detect its predictive ability to the immunotherapeutic response.

RESULTS

High infiltration levels of B cells heralded favorable prognosis, both in the TCGA-BLCA cohort and the local cohort (all P < 0.05). A 5-gene-pair model was established and served as a significant prognosis predictor across multiple cohorts (pooled hazard ratio = 2.79, 95% confidence interval = 2.22-3.49). The model could evaluate the prognosis effectively in 21 of 33 cancer types (P < 0.05). The signature was negatively associated with B cells' activation, proliferation, and infiltrating levels, and could serve as a potential predictor of immunotherapeutic outcomes.

CONCLUSIONS

A B cell-related gene signature was constructed to predict the prognosis and immunotherapeutic sensitivity in BLCA, helping to guide the personalized treatment.

Establishing a prognostic model based on five starvation-related long non-coding RNAs for clear cell renal cell carcinoma

BACKGROUND

Starvation-induced tumor microenvironment significantly alters genetic profiles including long non-coding RNAs (lncRNAs), further regulating the malignant biological characteristics (invasion and migration) of clear cell renal cell carcinoma (ccRCC).

METHODS

Transcriptome RNA-sequencing data of 539 ccRCC tumors and 72 normal tissues were acquired from the TCGA and paired clinical samples of 50 ccRCC patients. In vitro experiments, such as qPCR, migration and invasion assays were applied to reveal the clinical relevance of LINC-PINT, AC108449.2 and AC007637.1.

RESULTS

170 lncRNAs were verified as starvation-related lncRNAs (SR-LncRs), of which 25 lncRNAs were associated with overall survival in ccRCC patients. Furthermore, a starvation-related risk score model (SRSM) was built based on the expression levels of LINC-PINT, AC108449.2, AC009120.2, AC008702.2 and AC007637.1. ccRCC patients with high level of LINC-PINT expression were divided into high-risk group and led to higher mortality, but AC108449.2 and AC007637.1 were contrary. Analogously, LINC-PINT was highly expressed in ccRCC cell lines and tumor tissues, especially in patients with advanced stage, T-stage and M-stage, while AC108449.2 and AC007637.1 showed the opposite results. In addition, the increased levels of AC108449.2 and AC007637.1 were significantly correlated with grade. Silencing LINC-PINT reduced the invasion and migration characteristics of ccRCC cells. SiR-AC108449.2 and siR-AC007637.1 enhanced the ability of invasion and migration in ccRCC cells.

CONCLUSIONS

In this study, we find the clinical significance of LINC-PINT, AC108449.2 and AC007637.1 in predicting the prognosis of ccRCC patients and verify their correlation with various clinical parameters. These findings provide an advisable risk score model for ccRCC clinical decision-making.

T2DB: A Web Database for Long Non-Coding RNA Genes in Type II Diabetes

Type II diabetes (T2D) is a growing health problem worldwide due to increased levels of obesity and can lead to other life-threatening diseases, such as cardiovascular and kidney diseases. As the number of individuals diagnosed with T2D rises, there is an urgent need to understand the pathogenesis of the disease in order to prevent further harm to the body caused by elevated blood glucose levels. Recent advances in long non-coding RNA (lncRNA) research may provide insights into the pathogenesis of T2D. Although lncRNAs can be readily detected in RNA sequencing (RNA-seq) data, most published datasets of T2D patients compared to healthy donors focus only on protein-coding genes, leaving lncRNAs to be undiscovered and understudied. To address this knowledge gap, we performed a secondary analysis of published RNA-seq data of T2D patients and of patients with related health complications to systematically analyze the expression changes of lncRNA genes in relation to the protein-coding genes. Since immune cells play important roles in T2D, we conducted loss-of-function experiments to provide functional data on the T2D-related lncRNA USP30-AS1, using an in vitro model of pro-inflammatory macrophage activation. To facilitate lncRNA research in T2D, we developed a web application, T2DB, to provide a one-stop-shop for expression profiling of protein-coding and lncRNA genes in T2D patients compared to healthy donors or subjects without T2D.

A Novel Aging-Related Prognostic lncRNA Signature Correlated with Immune Cell Infiltration and Response to Immunotherapy in Breast Cancer

Breast cancer (BC) is among the most universal malignant tumors in women worldwide. Aging is a complex phenomenon, caused by a variety of factors, that plays a significant role in tumor development. Consequently, it is crucial to screen for prognostic aging-related long non-coding RNAs (lncRNAs) in BC. The BC samples from the breast-invasive carcinoma cohort were downloaded from The Cancer Genome Atlas (TCGA) database. The differential expression of aging-related lncRNAs (DEarlncRNAs) was screened by Pearson correlation analysis. Univariate Cox regression, LASSO-Cox analysis, and multivariate Cox analysis were performed to construct an aging-related lncRNA signature. The signature was validated in the GSE20685 dataset from the Gene Expression Omnibus (GEO) database. Subsequently, a nomogram was constructed to predict survival in BC patients. The accuracy of prediction performance was assessed through the time-dependent receiver operating characteristic (ROC) curves, Kaplan-Meier analysis, principal component analyses, decision curve analysis, calibration curve, and concordance index. Finally, differences in tumor mutational burden, tumor-infiltrating immune cells, and patients' response to chemotherapy and immunotherapy between the high- and low-risk score groups were explored. Analysis of the TCGA cohort revealed a six aging-related lncRNA signature consisting of MCF2L-AS1, USP30-AS1, OTUD6B-AS1, MAPT-AS1, PRR34-AS1, and DLGAP1-AS1. The time-dependent ROC curve proved the optimal predictability for prognosis in BC patients with areas under curves (AUCs) of 0.753, 0.772, and 0.722 in 1, 3, and 5 years, respectively. Patients in the low-risk group had better overall survival and significantly lower total tumor mutational burden. Meanwhile, the high-risk group had a lower proportion of tumor-killing immune cells. The low-risk group could benefit more from immunotherapy and some chemotherapeutics than the high-risk group. The aging-related lncRNA signature can provide new perspectives and methods for early BC diagnosis and therapeutic targets, especially tumor immunotherapy.

Integrative analysis to screen novel pyroptosis-related LncRNAs for predicting clinical outcome of glioma and validation in tumor tissue

BACKGROUND

Pyroptosis, also known as inflammatory necrosis, is a programmed cell death that manifests itself as a continuous swelling of cells until the cell membrane breaks, leading to the liberation of cellular contents, which triggers an intense inflammatory response. Pyroptosis might be a panacea for a variety of cancers, which include immunotherapy and chemotherapy-insensitive tumors such as glioma. Several findings have observed that long non-coding RNAs (lncRNAs) modulate the bio-behavior of tumor cells by binding to RNA, DNA and protein. Nevertheless, there are few studies reporting the effect of lncRNAs in pyroptosis processes in glioma.

METHODS

The principal goal of this study was to identify pyroptosis-related lncRNAs (PRLs) utilizing bioinformatic algorithm and to apply PCR techniques for validation in human glioma tissues. The second goal was to establish a prognostic model for predicting the overall survival patients with glioma. Predict algorithm was used to construct prognosis model with good diagnostic precision for potential clinical translation.

RESULTS

Noticeably, molecular subtypes categorized by the PRLs were not distinct from any previously published subtypes of glioma. The immune and mutation landscapes were obviously different from previous subtypes of glioma. Analysis of the sensitivity (IC50) of patients to 30 chemotherapeutic agents identified 22 agents as potential therapeutic agents for patients with low riskscores.

CONCLUSIONS

We established an exact prognostic model according to the expression profile of PRLs, which may facilitate the assessment of patient prognosis and treatment patterns and could be further applied to clinical.

Cytoskeleton regulator RNA expression on cancer-associated fibroblasts is associated with prognosis and immunotherapy response in bladder cancer

Background

Dysregulation of long noncoding RNAs (lncRNAs) has been reported to be associated with multiple tumors where they act as tumor suppressors or accelerators. The lncRNA CYTOR was identified as an oncogene involved in many cancers, such as gastric cancer, colorectal cancer, hepatocellular carcinoma, and renal cell carcinoma. However, the role of CYTOR in bladder cancer (BCa) has rarely been reported.

Methods

Using cancer datasets from The Cancer Genome Atlas (TCGA) program, we analyzed the association between CYTOR expression and prognostic value, oncogenic pathways, antitumor immunity and immunotherapy response in BCa. The influence of CYTOR on the immune infiltration pattern in the urothelial carcinoma microenvironment was further verified in our dataset. Single-cell analysis revealed the role of CYTOR in the tumor microenvironment (TME) of BCa. Finally, we evaluated the expression of CYTOR in BCa in the Peking University First Hospital (PKU-BCa) dataset and its correlation with the malignant phenotype of BCa in vitro and in vivo.

Results

The results indicated that CYTOR was highly expressed in multiple cancer samples, including BCa, and increased CYTOR expression contributed to poor overall survival (OS). Additionally, elevated CYTOR expression was significantly correlated with clinicopathological features of BCa, such as female sex, advanced TNM stage, high histological grade and non-papillary subtype. Functional characterization revealed that CYTOR may be involved in immune-related pathways and the epithelial mesenchymal transformation (EMT) process. Moreover, CYTOR had a significant association with infiltrating immune cells, including M2 macrophages and regulatory T cells (Tregs). CYTOR facilitates the crosstalk between cancer-associated fibroblasts (CAFs) and macrophages, and mediates M2 polarization of macrophages. Correlation analysis revealed a positive correlation between CYTOR expression and programmed cell death-1 (PD-1)/programmed death ligand 1 (PD-L1)/expression and other targets for specific immunotherapy in BCa, which are recognized to predict the efficacy of immunotherapy.

Conclusions

These results suggest that CYTOR serves as a potential biomarker for predicting survival outcome, TME cell infiltration characteristics and immunotherapy response in BCa.

Necroptosis-Related LncRNA Signatures for Prognostic Prediction in Uterine Corpora Endometrial Cancer

Necroptosis is one of the common modes of apoptosis, and it has an intrinsic association with cancer prognosis. However, the role of the necroptosis-related long non-coding RNA LncRNA (NRLncRNAs) in uterine corpora endometrial cancer (UCEC) has not yet been fully elucidated at present. Therefore, the present study is designed to investigate the potential prognostic value of necroptosis-related LncRNAs in UCEC. In the present study, the expression profiles and clinical data of UCEC patients were downloaded from TCGA database to identify the differentially expressed NRLncRNAs associated with overall survival. A LncRNA risk model was constructed via Cox regression analysis, and its prognostic value was evaluated. We have also further evaluated the relationships between the LncRNA features and the related cellular function, related pathways, immune status, and immune checkpoints m6A-related genes. Seven signatures, including PCAT19, CDKN2B-AS1, LINC01936, LINC02178, BMPR1B-DT, LINC00237, and TRPM2-AS, were established to assess the overall survival (OS) of the UCEC in the present study. Survival analysis and ROC curves indicated that the correlated signature has good predictable performance. The normogram could accurately predict the overall survival of the patients with an excellent clinical practical value. Enrichment analysis of gene sets indicated that risk signals were enriched in several immune-related pathways. In addition, the risk characteristics were significantly correlated with immune cells, immune function, immune cell infiltration, immune checkpoints, and some m6A-related genes. This study has identified seven necroptosis-related LncRNA signatures for the first time, providing a valuable basis for a more accurate prognostic prediction of UCEC.

Cuproptosis-related LncRNAs are correlated with immunity and predict prognosis in HNSC independent of TMB

Aims: Cuproptosis is a novel cell death pathway, and the regulatory mechanism in head and neck squamous cell carcinoma (HNSC) remains to be explored. We determined whether cuproptosis-related lncRNAs (CRLs) could predict prognosis in HNSC. Methods and Results: First, we identified 10 prognostic CRLs by Pearson correlation and univariate Cox regression analyses. Next, we constructed the CRLs prognostic model based on 5 CRLs screened by the least absolute shrinkage and selection operator (LASSO) Cox analysis. Following this, we calculated the risk score for HNSC patients and divided patients into high- and low-risk groups. In our prognostic model, HNSC patients with higher risk scores had poorer outcomes. Based on several prognostic features, a predictive nomogram was established. Furthermore, we investigated principal component analysis to distinguish two groups, and functional enrichment analysis of 176 differentially expressed genes (DEGs) between risk groups was performed. Finally, we analyzed relationships between tumor mutation burden (TMB) and risk scores. Conclusion: Cuproptosis-related lncRNAs can be applied to predict HNSC prognosis independent of TMB, which is closely correlated with tumor immunity.

ANLN Serves as an Oncogene in Bladder Urothelial Carcinoma via Activating JNK Signaling Pathway

INTRODUCTION

To understand the significance of ANLN (anillin, actin-binding protein)-mediated c-Jun N-terminal kinase (JNK) signal pathway on the progression of bladder urothelial carcinoma (BLCA).

METHODS

The Cancer Genome Atlas (TCGA) database was utilized to perform the clinical significance of ANLN in BLCA. Then, ANLN expression was determined in human normal primary bladder epithelial cells (BdEC) and BLCA cells. Later, ANLN knockdown was performed in BLCA cells, where the expression of MAPK8, MAPK9, and p-JNK/JNK was detected. BLCA cells were divided into the Mock, siNC, siANLN, SP600125 (a selective JNK inhibitor), and ANLN + SP600125 group, followed by measurements of real-time quantitative polymerase chain reaction, 3-4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide, Annexin V-FITC/PI, Wound-healing, Transwell, and immunofluorescence assays.

RESULTS

ANLN was upregulated in the BLCA tissues, which showed a relation with the stage of patients. Besides, BLCA patients with high expression of ANLN had a worse prognosis than those with low expression of ANLN. Besides, the expression of ANLN in the BLCA tissues was positively correlated with MAPK8 and MAPK9. SP600125 suppressed the JNK signal pathway, reduced the proliferation, and increased BLCA cell apoptosis, with the reductions in the invasion and migration and the upregulation of phospho-histone H3 Ser-10 (pHH3), which was abolished by the overexpression of ANLN.

CONCLUSION

ANLN, as an oncogene of BLCA, may associate with the activation of JNK signal pathway. Inhibiting ANLN could deactivate the JNK signal pathway, thereby suppressing the proliferation, invasion, and migration while promoting the apoptosis of BLCA cells.

Development and validation of an endoplasmic reticulum stress long non-coding RNA signature for the prognosis and immune landscape prediction of patients with lung adenocarcinoma

Background: Lung adenocarcinoma (LUAD), the most common histotype of lung cancer, may have variable prognosis due to molecular variations. This work investigated long non-coding RNA (lncRNA) related to endoplasmic reticulum stress (ERS) to predict the prognosis and immune landscape for LUAD patients. Methods: RNA data and clinical data from 497 LUAD patients were collected in the Cancer Genome Atlas database. Pearson correlation analysis, univariate Cox regression, least absolute shrinkage and selection operator regression analyses, as well as the Kaplan-Meier method, were used to screen for ERS-related lncRNAs associated with prognosis. The risk score model was developed using multivariate Cox analysis to separate patients into high- and low-risk groups and a nomogram was constructed and evaluated. Finally, we explore the potential functions and compared the immune landscapes of two groups. Quantitative real-time PCR was used to verify the expression of these lncRNAs. Results: Five ERS-related lncRNAs were shown to be strongly linked to patients' prognosis. A risk score model was built by using these lncRNAs to categorize patients based on their median risk scores. For LUAD patients, the model was found to be an independent prognostic predictor (p < 0.001). The signature and clinical variables were then used to construct a nomogram. With 3-year and 5-year OS' AUC of 0.725 and 0.740, respectively, the nomogram's prediction performance is excellent. The 5-lncRNA signature was associated with DNA replication, epithelial-mesenchymal transition, and the pathway of cell cycle, P53 signaling. Between the two risk groups, immune responses, immune cells, and immunological checkpoints were found to be considerably different. Conclusion: Overall, our findings indicate that the 5 ERS-related lncRNA signature was an excellent prognostic indicator and helped to predict the immunotherapy response for patients with LUAD.

A prognostic risk prediction model based on ferroptosis-related long non-coding RNAs in bladder cancer: A bulk RNA-seq research and scRNA-seq validation

BACKGROUND

To construct a prognostic risk model of bladder cancer (BC) from the perspective of long non-coding RNAs (lncRNAs) and ferroptosis, in order to guide clinical prognosis and identify potential therapeutic targets.

METHODS

In-hours BC samples were collected from 4 patients diagnosed with BC, who underwent radical cystectomy. Single cell transcriptome sequencing was performed and Seurat package were used for quality control and secondary analysis. LncRNAs expression profiles of BC samples were extracted from The Cancer Genome Atlas database. And sex, age, tumor, node, metastasis stage and other clinical data was downloaded at the same time. Ferroptosis-related lncRNAs were identified by co-expression analysis. We constructed a risk model by Cox regression and least absolute shrinkage and selection operator regression analyses. The predictive strength of the risk model for overall survival (OS) of patients with BC was evaluated by the log-rank test and Kaplan-Meier method. Finally, the enrichment analysis was performed and visualized.

RESULTS

We identified and included 15 prognostic ferroptosis-related lncRNAs (AL356740.1, FOXC2AS1, ZNF528AS1, LINC02535, PSMB8AS1, AL590428.1, AP000347.2, OCIAD1-AS1, AP001347.1, AC104986.2, AC018926.2, LINC00867, AC099518.4, USP30-AS1, and ARHGAP5-AS1), to build our ferroptosis-related lncRNAs risk model. Using this risk model, BC patients were divided into high and low-risk groups, and their respective survival lengths were calculated. The results showed that the OS of the low-risk group was significantly longer than that of the high-risk group. A nomogram was utilized to predict the survival rate of BC patients. As indicated in the nomogram, risk score was the most important indicator of OS in patients with BC. The ferroptosis-related lncRNAs risk model is an independent tool for prognostic risk assessment in patients with BC. Single cell transcriptome sequencing suggests that ferroptosis-related lncRNAs express specifically in BC tumor microenvironment. AL356740.1, LINC02535 and LINC00867 were mainly expressed in tumor cells.

CONCLUSION

The risk model based on the ferroptosis-related lncRNAs and the genomic clinico-pathological nomogram could be used to accurately predict the prognosis of patients with BC. The lncRNAs used to build this model might become potential therapeutic targets in the future.

A novel autophagy-related long non-coding RNAs prognostic risk score for clear cell renal cell carcinoma

BACKGROUND

As the main histological subtype of renal cell carcinoma, clear cell renal cell carcinoma (ccRCC) places a heavy burden on health worldwide. Autophagy-related long non-coding RNAs (ARlncRs) have shown tremendous potential as prognostic signatures in several studies, but the relationship between them and ccRCC still has to be demonstrated.

METHODS

The RNA-sequencing and clinical characteristics of 483 ccRCC patients were downloaded download from the Cancer Genome Atlas and International Cancer Genome Consortium. ARlncRs were determined by Pearson correlation analysis. Univariate and multivariate Cox regression analyses were applied to establish a risk score model. A nomogram was constructed considering independent prognostic factors. The Harrell concordance index calibration curve and the receiver operating characteristic analysis were utilized to evaluate the nomogram. Furthermore, functional enrichment analysis was used for differentially expressed genes between the two groups of high- and low-risk scores.

RESULTS

A total of 9 SARlncRs were established as a risk score model. The Kaplan-Meier survival curve, principal component analysis, and subgroup analysis showed that low overall survival of patients was associated with high-risk scores. Age, M stage, and risk score were identified as independent prognostic factors to establish a nomogram, whose concordance index in the training cohort, internal validation, and external ICGC cohort was 0.793, 0.671, and 0.668 respectively. The area under the curve for 5-year OS prediction in the training cohort, internal validation, and external ICGC cohort was 0.840, 0.706, and 0.708, respectively. GO analysis and KEGG analysis of DEGs demonstrated that immune- and inflammatory-related pathways are likely to be critically involved in the progress of ccRCC.

CONCLUSIONS

We established and validated a novel ARlncRs prognostic risk model which is valuable as a potential therapeutic target and prognosis indicator for ccRCC. A nomogram including the risk model is a promising clinical tool for outcomes prediction of ccRCC patients and further formulation of individualized strategy.

Prognostic implications of necroptosis-related long noncoding RNA signatures in muscle-invasive bladder cancer

Background: Bladder cancer (BLCA) is the sixth most common cancer in men, with an increasing incidence of morbidity and mortality. Necroptosis is a type of programmed cell death and plays a critical role in the biological processes of bladder cancer (BLCA). However, current studies focusing on long noncoding RNA (lncRNA) and necroptosis in cancer are limited, and there is no research about necroptosis-related lncRNAs (NRLs) in BLCA. Methods: We obtained the RNA-seq data and corresponding clinical information of BLCA from The Cancer Genome Atlas (TCGA) database. The seven determined prognostic NLRs were analyzed by several methods and verified by RT-qPCR. Then, a risk signature was established based on the aforementioned prognostic NLRs. To identify it, we evaluated its prognostic value by Kaplan-Meier (K-M) survival curve and receiver operating characteristics (ROC) curve analysis. Moreover, the relationships between risk signature and clinical features, functional enrichment, immune landscape, and drug resistance were explored as well. Results: We constructed a signature based on seven defined NLRs (HMGA2-AS1, LINC02489, ETV7-AS1, EMSLR, AC005954.1, STAG3L5P-PVRIG2P-PILRB, and LINC02178). Patients in the low-risk cohort had longer survival times than those in the high-risk cohort, and the area under the ROC curve (AUC) value of risk signature was higher than other clinical variables. Functional analyses, the infiltrating level of immune cells and functions, ESTIMATE score, and immune checkpoint analysis all indicated that the high-risk group was in a relatively immune-activated state. In terms of treatments, patients in the high-risk group were more sensitive to immunotherapy, especially anti-PD1/PD-L1 immunotherapy and conventional chemotherapy. Conclusion: The novel NLR signature acts as an invaluable tool for predicting prognosis, immune microenvironment, and drug resistance in muscle-invasive bladder cancer (MIBC) patients.

A New Prognostic Signature Constructed with Necroptosis-Related lncRNA in Bladder Cancer

BACKGROUND

Bladder cancer (BC) accounts for the most common urologic malignancy, leading to a heavy social burden over the world. We aim to search for a novel prognostic biomarker with necroptosis-related lncRNAs of bladder cancer in this study.

METHODS

We download the RNA-sequencing data and corresponding clinical information of BC patients from TCGA. We performed Pearson correlation analysis to identify necroptosis-related lncRNAs (NRlncRNAs). Then, we used univariate Cox regression, Lasso Cox analysis, and multivariate Cox regression to construct the optimal prognostic model. Next, we used Kaplan-Meier curves, Cox regression, receiver operating characteristic (ROC) curves, nomogram, and stratified survival analysis to evaluate the capacity of the prognostic signature. Furthermore, gene set enrichments in the signature and the correlation between prognostic signature and necroptosis genes, tumor microenvironment, immune infiltration, and immune checkpoints of BC were also explored.

RESULTS

A 7-NRlncRNAs signature comprising FKBP14-AS1, AL731567.1, LINC02178, AC011503.2, LINC02195, AC068196.1, and AL136084.2 was constructed to predict the prognosis of BC in this research. Cox regression analysis showed that the signature could be an independent prognostic factor for BC patients (P < 0.001). Compared to other clinicopathological characteristics, this signature displayed a better capacity of prediction with the area under the curve (AUC) of 0.745. Stratified analysis using various clinical variables demonstrated that the prognostic signature has good clinical fitness. GSEA showed that focal adhesion and the WNT signaling pathway were enriched in the high-risk group. Immune infiltration analysis indicated that the signature was significantly inversely correlated with infiltration of CD8⁺ T cells and CD4⁺ T cells while positively correlated with macrophages and cancer associated fibroblasts. Immune checkpoint analysis revealed that the expressions of protective factors were significantly lower in the high-risk group, while expressions of cancer promotors were significantly higher in this group. The gene expression analysis displayed that necroptosis genes such as FADD, FAS, MYC, STAT3, PLK1, LEF1, EGFR, RIPK3, CASP8, BRAF, ID1, GATA3, MYCN, CD40, and TNFRSF21 were significantly different between the two groups.

CONCLUSIONS

The 7-NRlncRNAs signature can predict the overall survival of BC and may provide help for the individualized treatment of BC patients.

Development and validation of a novel hypoxia-related signature for prognostic and immunogenic evaluation in head and neck squamous cell carcinoma

Hypoxia plays a critical role in head and neck squamous cell carcinoma (HNSCC) prognosis. However, till now, robust and reliable hypoxia-related prognostic signatures have not been established for an accurate prognostic evaluation in HNSCC patients. This article focused on establishing a risk score model to evaluate the prognosis and guide treatment for HNSCC patients. RNA-seq data and clinical information of 502 HNSCC patients and 44 normal samples were downloaded from The Cancer Genome Atlas (TCGA) database. 433 samples from three Gene Expression Omnibus (GEO) datasets were incorporated as an external validation cohort. In the training cohort, prognostic-related genes were screened and LASSO regression analyses were performed for signature establishment. A scoring system based on SRPX, PGK1, STG1, HS3ST1, CDKN1B, and HK1 showed an excellent prediction capacity for an overall prognosis for HNSCC patients. Patients were divided into high- and low-risk groups, and the survival status of the two groups exhibited a statistically significant difference. Subsequently, gene set enrichment analysis (GSEA) was carried out to explore the underlying mechanisms for the prognosis differences between the high- and low-risk groups. The tumor immune microenvironment was evaluated by CIBERSORT, ESTIMATE, TIDE, and xCell algorithm, etc. Then, we explored the relationships between this prognostic model and the levels of immune checkpoint-related genes. Cox regression analysis and nomogram plot indicated the scoring system was an independent predictor for HNSCC. Moreover, a comparison of predictive capability has been made between the present signature and existing prognostic signatures for HNSCC patients. Finally, we detected the expression levels of proteins encoded by six-HRGs via immunohistochemical analysis in tissue microarray. Collectively, a novel integrated signature considering both HRGs and clinicopathological parameters will serve as a prospective candidate for the prognostic evaluation of HNSCC patients.

M7G-Related lncRNAs predict prognosis and regulate the immune microenvironment in lung squamous cell carcinoma

Background

N7-Methylguanosine (m7G) and long non-coding RNAs (lncRNAs) have been widely studied in cancer and have been found to be useful for assessing tumor progression. However, the role of m7G-related lncRNAs in lung squamous cell carcinoma (LUSC) remains unclear. Thus, it is crucial to identify m7G-associated lncRNAs with definitive prognostic value. This study aimed to investigate the prognostic value, correlation with tumor mutation burden, and impact on the tumor immune microenvironment of m7G-related lncRNAs in LUSC. 

Methods

LUSC transcriptome data and clinical data were downloaded from The Cancer Genome Atlas, and an m7G-related lncRNA-mRNA co-expression network was constructed using Pearson’s correlation analysis. Cox regression analyses were used to determine a risk model for m7G-associated lncRNAs with prognostic value. The risk signature was verified using the Kaplan–Meier method, receiver operating characteristic curve analysis, and principal component analysis. A nomogram based on risk scores and clinical characteristics was then developed. Gene set enrichment analysis was used for functional annotation to analyze the risk signature. The association among the risk signature, tumor mutational burden, and tumor-infiltrating immune cells was then analyzed. RT-qPCR was used to investigate the expression of 6 m7G-related lncRNAs in LUSC cells. The cytological function of SRP14-AS1 was verified by wound-healing assay and transwell assay.

Results

A total of 293 m7G-related lncRNAs were identifed, 27 candidate m7G-related lncRNAs were signifcantly associated with overall survival (OS). Six of these lncRNAs (CYP4F26P, LINC02178, MIR22HG, SRP14-AS1, TMEM99, PTCSC2) were selected for establishment of the risk model. The OS of patients in the low-risk group was higher than that of patients in the high-risk group (p < 0.001). Multivariate cox regression analysis indicated that the model could be an independent prognostic factor for LUSC (HR = 1.859; 95% CI 1.452–2.380, p < 0.001). The ROC curve analysis revealed that the AUCs for OS in the 3-, and 5-year were 0.682, 0.657, respectively. GSEA analysis revealed that the risk model was closely related to immune-related pathways. Compared with normal lung epithelial cells, four m7G-related lncRNAs were higher expressed in cancer cells and two were lower expressed, among which knockdown of SRP14-AS1 promoted the proliferation and migration of LUSC cells.

Conclusion

A risk model based on six m7G-related lncRNAs with prognostic value may be a promising prognostic tool in LUSC and guide individualized patient treatment.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-10232-z.

Identification of novel cuproptosis-related lncRNA signatures to predict the prognosis and immune microenvironment of breast cancer patients

Background

Cuproptosis is a new modality of cell death regulation that is currently considered as a new cancer treatment strategy. Nevertheless, the prognostic predictive value of cuproptosis-related lncRNAs in breast cancer (BC) remains unknown. Using cuproptosis-related lncRNAs, this study aims to predict the immune microenvironment and prognosis of BC patients. and develop new therapeutic strategies that target the disease.

Methods

The Cancer Genome Atlas (TCGA) database provided the RNA-seq data along with the corresponding clinical and prognostic information. Univariate and multivariate Cox regression analyses were performed to acquire lncRNAs associated with cuproptosis to establish predictive features. The Kaplan-Meier method was used to calculate the overall survival rate (OS) in the high-risk and low-risk groups. High risk and low risk gene sets were enriched to explore functional discrepancies among risk teams. The mutation data were analyzed using the “MAFTools” r-package. The ties of predictive characteristics and immune status had been explored by single sample gene set enrichment analysis (ssGSEA). Last, the correlation between predictive features and treatment condition in patients with BC was analyzed. Based on prognostic risk models, we assessed associations between risk subgroups and immune scores and immune checkpoints. In addition, drug responses in at-risk populations were predicted.

Results

We identified a set of 11 Cuproptosis-Related lncRNAs (GORAB-AS1, AC 079922.2, AL 589765.4, AC 005696.4, Cytor, ZNF 197-AS1, AC 002398.1, AL 451085.3, YTH DF 3-AS1, AC 008771.1, LINC 02446), based on which to construct the risk model. In comparison to the high-risk group, the low-risk patients lived longer (p < 0.001). Moreover, cuproptosis-related lncRNA profiles can independently predict prognosis in BC patients. The AUC values for receiver operating characteristics (ROC) of 1-, 3-, and 5-year risk were 0.849, 0.779, and 0.794, respectively. Patients in the high-risk group had lower OS than those in the low-risk group when they were divided into groups based on various clinicopathological variables. The tumor burden mutations (TMB) correlation analysis showed that high TMB had a worse prognosis than low-TMB, and gene mutations were found to be different in high and low TMB groups, such as PIK3CA (36% versus 32%), SYNE1 (4% versus 6%). Gene enrichment analysis indicated that the differential genes were significantly concentrated in immune-related pathways. The predictive traits were significantly correlated with the immune status of BC patients, according to ssGSEA results. Finally, high-risk patients showed high sensitivity in anti-CD276 immunotherapy and conventional chemotherapeutic drugs such as imatinib, lapatinib, and pazopanib.

Conclusion

We successfully constructed of a cuproptosis-related lncRNA signature, which can independently predict the prognosis of BC patients and can be used to estimate OS and clinical treatment outcomes in BRCA patients. It will serve as a foundation for further research into the mechanism of cuproptosis-related lncRNAs in breast cancer, as well as for the development of new markers and therapeutic targets for the disease.

Emerging role of lncRNAs in drug resistance mechanisms in head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) originates in the squamous cell lining the mucosal surfaces of the head and neck region, including the oral cavity, nasopharynx, tonsils, oropharynx, larynx, and hypopharynx. The heterogeneity, anatomical, and functional characteristics of the patient make the HNSCC a complex and difficult-to-treat disease, leading to a poor survival rate and a decreased quality of life due to the loss of important physiologic functions and aggressive surgical injury. Alteration of driver-oncogenic and tumor-suppressing lncRNAs has recently been recently in HNSCC to obtain possible biomarkers for diagnostic, prognostic, and therapeutic approaches. This review provides current knowledge about the implication of lncRNAs in drug resistance mechanisms in HNSCC. Chemotherapy resistance is a major therapeutic challenge in HNSCC in which lncRNAs are implicated. Lately, it has been shown that lncRNAs involved in autophagy induced by chemotherapy and epithelial–mesenchymal transition (EMT) can act as mechanisms of resistance to anticancer drugs. Conversely, lncRNAs involved in mesenchymal–epithelial transition (MET) are related to chemosensitivity and inhibition of invasiveness of drug-resistant cells. In this regard, long non-coding RNAs (lncRNAs) play a pivotal role in both processes and are important for cancer detection, progression, diagnosis, therapy response, and prognostic values. As the involvement of more lncRNAs is elucidated in chemoresistance mechanisms, an improvement in diagnostic and prognostic tools could promote an advance in targeted and specific therapies in precision oncology.

Ferroptosis-Related lncRNA Signature Correlates with the Prognosis, Tumor Microenvironment, and Therapeutic Sensitivity of Esophageal Squamous Cell Carcinoma

Esophageal squamous cell carcinoma (ESCC) is the most prevalent form of esophageal cancer in China and is closely associated with malignant biological characteristics and poor survival. Ferroptosis is a newly discovered iron-dependent mode of cell death that plays an important role in the biological behavior of ESCC cells. The clinical significance of ferroptosis-related long noncoding RNAs (FRLs) in ESCC remains unknown and warrants further research. The current study obtained RNA sequencing profiles and corresponding clinical data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, and FRLs were obtained through coexpression analysis. Consensus clustering was employed to divide the subjects into clusters, and immune-associated pathways were identified by functional analysis. The current study observed significant differences in the enrichment scores of immune cells among different clusters. Patients from TCGA-ESCC database were designated as the training cohort. A ten-FRL prediction signature was established using the least absolute shrinkage and selection operator Cox regression model and validated using the GEO cohort and our own independent validation database. Real-time quantitative polymerase chain reaction was used to verify the expression of the ten FRLs, and the ssGSEA analysis was employed to evaluate their function. In addition, the IMvigor database was used to assess the predictive value of the signature in terms of immunotherapeutic responses. Multivariate Cox and stratification analyses revealed that the ten-FRL signature was an independent predictor of the overall survival (OS). Patients with ESCC in the high-risk group displayed worse survival, a characteristic tumor immune microenvironment, and low immunotherapeutic benefits compared to those in the low-risk group. Collectively, the risk model established in this study could serve as a promising predictor of prognosis and immunotherapeutic response in patients with ESCC.

The Interferon Gamma-Related Long Noncoding RNA Signature Predicts Prognosis and Indicates Immune Microenvironment Infiltration in Colon Adenocarcinoma

Colon adenocarcinoma (COAD) is one of the most common clinically malignant tumours of the digestive system, with high incidence and mortality and poor prognosis. Interferon-gamma (IFN-γ) and long noncoding RNAs (lncRNAs) have prognostic values and were closely associated with immune microenvironment in COAD. Thus, identifying IFN-γ-related lncRNAs may be valuable in predicting the survival of patients with COAD. In this study, we identified IFN-γ-related lncRNAs and divided COAD patients from the Cancer Genome Atlas (TCGA) database into training and validation sets. Pearson’s correlation analysis and least absolute shrinkage and selection operator (LASSO) Cox regression were performed to select IFN-γ-related lncRNA-associated prognoses. Thirteen lncRNAs (AC025165.8, AC091633.3, FENDRR, LINC00882, LINC01828, LINC01829, MYOSLID, RP11-154H23.4, RP11-20J15.3, RP11-324L17.1, RP11-342A23.2, RP11-805I24.3, SERTAD4-AS1) were identified to construct an IFN-γ-related lncRNA prognostic signature in TCGA training (n =213) and validation (n =213) cohorts. COAD patient risk scores were calculated and classified into high- and low-risk groups based on the median value of the risk scores in each dataset. We compared the overall survival (OS) of patients stratified by age, gender, and stage. The OS in the high-risk group was significantly shorter than that in the low-risk group. In addition, the clinical nomogram incorporating the prognostic signature and clinical features showed a high concordance index of 0.78 and accurately predicted 1-, 3-, and 5-year survival times among COAD patients in the high- and low-risk groups. Based on the risk model, the high- and low-risk groups exhibited distinct differences in the immune system by gene set enrichment analysis (GSEA) functional annotation, and differentially expressed genes (DEGs) between the high- and low-risk groups were subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. We investigated the expression of multiple immune checkpoint genes in the high- and low-risk groups and plotted Kaplan-Meier survival curves, indicating that immune checkpoint genes, such as LAG3 and PD. L1, STING and TIM 3, were also expressed differently between the two risk groups. Subsequently, there were dramatic differences in mutated genes, SNV (single nucleotide variants) classes, variant types and variant allele frequencies between low- and high-risk patients with COAD. Patients stratified by risk scores had different sensitivities to common chemotherapeutic agents. Finally, we used quantitative real-time polymerase chain reaction (qRT-PCR) assays to demonstrate that three lncRNAs were significantly differentially expressed in COAD tissues and adjacent normal tissues. Considered together, a thirteen-lncRNA prognostic signature has great potential to be a prognostic biomarker and could play an essential role in the immune microenvironment of COAD.

Clinicopathological and Prognostic Value of Necroptosis-Associated lncRNA Model in Patients with Kidney Renal Clear Cell Carcinoma

Background. Necroptosis, a recently identified type of programmed necrotic cell death, is closely related to the tumorigenesis and development of cancer. However, it remains unclear whether necroptosis-associated long noncoding RNAs (lncRNAs) can be used to predict the prognosis of kidney renal clear cell carcinoma (KIRC). This work was designed to probe the possible prognostic worth of necroptosis-associated lncRNAs along with their impact on the tumor microenvironment (TME) in KIRC. Methods. The Cancer Genome Atlas (TCGA) database was used to extract KIRC gene expression and clinicopathological data. Pearson correlation analysis was used to evaluate necroptosis-associated lncRNAs against 159 known necroptosis-associated genes. To define molecular subtypes, researchers used univariate Cox regression analysis and consensus clustering, as well as clinical significance, TME, and tumor immune cells in each molecular subtype. We develop the necroptosis-associated lncRNA prognostic model using univariate Cox regression analysis and least absolute shrinkage and selection operator (LASSO) regression analysis. Patients were divided into high- and low-risk groups according to prognostic model. Moreover, comprehensive analyses, including prognostic value, gene set enrichment analysis (GSEA), immune infiltration, and immune checkpoint gene expression, were performed between the two risk groups. Finally, anticancer drug sensitivity analyses were employed for assessing associations for necroptosis-associated lncRNA expression profile and anticancer drug chemosensitivity. Results. Through univariate analysis, sixty-nine necroptosis-associated lncRNAs were found to have a significant relationship with KIRC prognosis. Two molecular clusters were identified, and significant differences were found with respect to clinicopathological features and prognosis. The segregation of patients into two risk groups was done by the constructed necroptosis-associated lncRNA model. The survival prognosis, clinical features, degree of immune cell infiltration, and expression of immune checkpoint genes of high-risk and low-risk groups were all shown to vary. Conclusions. Our study identified a model of necroptosis-associated lncRNA signature and revealed its prognostic role in KIRC. It is expected to provide a reference for the screening of KIRC prognostic markers and the evaluation of immune response.

Establishment of a Cell Necroptosis Index to Predict Prognosis and Drug Sensitivity for Patients With Triple-Negative Breast Cancer

Background: Necroptosis has been an alternatively identified mechanism of programmed cancer cell death, which plays a significant role in cancer. However, research about necroptosis-related long noncoding RNAs (lncRNAs) in cancer are still few. Moreover, the potentially prognostic value of necroptosis-related lncRNAs and their correlation with the immune microenvironment remains unclear. The present study aimed to explore the potential prognostic value of necroptosis-related lncRNAs and their relationship to immune microenvironment in triple-negative breast cancer (TNBC).

Methods: The RNA expression matrix of patients with TNBC was obtained from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases. Finally, 107 patients of GSE58812, 159 patients of TCGA, and 143 patients of GSE96058 were included. Necroptosis-related lncRNAs were screened by Cox regression and Pearson correlation analysis with necroptosis-related genes. By LASSO regression analysis, nine necroptosis-related lncRNAs were employed, and a cell necroptosis index (CNI) was established; then, we evaluated its prognostic value, clinical significance, pathways, immune infiltration, and chemotherapeutics efficacy.

Results: Based on the CNI value, the TNBC patients were divided into high- and low-CNI groups, and the patients with high CNI had worse prognosis, more lymph node metastasis, and larger tumor (p < 0.05). The receiver operating characteristic (ROC) analysis showed that the signature performed well. The result of the infiltration proportion of different immune cell infiltration further explained that TNBC patients with high CNI had low immunogenicity, leading to poor therapeutic outcomes. Moreover, we found significant differences of the IC50 values of various chemotherapeutic drugs in the two CNI groups, which might provide a reference to make a personalized chemotherapy for them.

Conclusion: The novel prognostic marker CNI could not only precisely predict the survival probability of patients with TNBC but also demonstrate a potential role in antitumor immunity and drug sensitivity.

Multi-omics reveals that Bifidobacterium breve M-16V may alleviate the immune dysregulation caused by nanopolystyrene

There is a growing attention regarding the toxic effect of microplastics pollutants. However, comprehensive phenotyping- and omics-based strategies for the toxicity evaluation of microplastics on the host remain to be established. To this end, we designed an encompassing phenotyping and multi-omics analysis method to detect the molecular interference of nanopolystyrene (PS)-exposed mice. The exposure time was 28 days with 1000 μg/L PS. We found that PS induced microbial alteration and metabolic disorders, which was closely related to immune disturbances. In addition, the altered expression of some genes related to immune dysregulation was observed. Interestingly, Bifidobacterium breve M-16V (B. breve M-16V) significantly inhibited Th2 and Th17 lymphocyte subset. Simultaneously, B.breve M-16V may activate MyD88 expression and promote Th1-related cytokine IL-12 production. In addition, B. breve M-16V may partially restore the gut microbiota dysbiosis. In summary, we demonstrated that the combined phenotyping and omics-based profiling established a practical framework that allowed us to gain a deeper understanding of the maladaptive consequences of PS exposure. It can be utilized to evaluate the toxicity of other environmental microplastics pollutants. Meanwhile, we found that B. breve M-16V has certain anti-inflammatory and immunomodulatory functions through host-microbiome interactions.

A novel prognostic model associated with the overall survival in patients with breast cancer based on lipid metabolism-related long noncoding RNAs

BACKGROUND

Lipid metabolism is closely related to the occurrence and development of breast cancer. Our purpose was to establish a novel model based on lipid metabolism-related long noncoding RNAs (lncRNAs) and evaluate the potential clinical value in predicting prognosis for patients suffering from breast cancer.

METHODS

RNA data and clinical information for breast cancer were obtained from the cancer genome atlas (TCGA) database. Lipid metabolism-related lncRNAs were identified via the criteria of correlation coefficient |R² | > 0.4 and p < 0.001, and prognostic lncRNAs were identified to establish model through Cox regression analysis. The training set and validation set were established to certify the feasibility, and all samples were separated into high-risk group or low-risk group. Gene Ontology (GO) and Gene Set Enrichment Analysis (GSEA) were conducted to evaluate the potential biological functions, and the immune infiltration levels were explored through Cibersortx database.

RESULTS

A total of 14 lncRNAs were identified as protective genes (AC022150.4, AC061992.1, AC090948.3, AC092794.1, AC107464.3, AL021707.8, AL451085.2, AL606834.2, FLJ42351, LINC00926, LINC01871, TNFRSF14-AS1, U73166.1 and USP30-AS1) with HRs < 1 while 10 lncRNAs (AC022150.2, AC090948.1, AC243960.1, AL021707.6, ITGB2-AS1, OTUD6B-AS1, SP2-AS1, TOLLIP-AS1, Z68871.1 and ZNF337-AS1) were associated with increased risk with HRs >1. A total of 24 prognostic lncRNAs were selected to construct the model. The patients in low-risk group were associated with better prognosis in both training set (p < 0.001) and validation set (p < 0.001). The univariate and multivariate Cox regression analyses revealed that risk score was an independent prognostic factors in both training set (p < 0.001) and validation set (p < 0.001). GO and GSEA analyses revealed that these lncRNAs were related to metabolism-related signal pathway and immune cells signal pathway. Risk score was negatively correlated with B cells (r = -0.097, p = 0.002), NK cells (r = -0.097, p = 0.002), Plasma cells (r = -0.111, p = 3.329e-04), T-cells CD4 (r = -0.064, p = 0.039) and T-cells CD8 (r = -0.322, p = 2.357e-26) and positively correlated with Dendritic cells (r = 0.077, p = 0.013) and Monocytes (r = 0.228, p = 1.107e-13).

CONCLUSION

The prognostic model based on lipid metabolism lncRNAs possessed an important value in survival prediction of breast cancer patients.

Clinical Value Analysis of Xiaozheng Decoction Combined with Bladder Perfusion for Postoperative Treatment of Bladder Cancer and Its Effect on Serum miR-143 and miR-92a

Objective

To study the clinical value of Xiaozheng decoction combined with bladder perfusion treating bladder cancer after the operation and its effect on serum miR-143 and miR-92a.

Methods

The patients in the control group were treated with gemcitabine bladder instillation, and patients in the study group were treated with the combination of gemcitabine bladder instillation + Xiaozheng decoction. The clinical efficacy, postoperative adverse effects, and recurrence between the two groups were compared. miR-143 and miR-92a levels, immune function levels, and tumor factor levels were compared before and after treatment. The relationship between patient prognosis and miR-143/miR-92a expression was analyzed.

Results

The overall effective rate of treatment in the study group (86.67%) was significantly increased, and the occurrence of adverse reactions and recurrence were significantly decreased (P < 0.05). After treatment, serum miR-143 and miR-92a levels, CD4+, CD4+/CD8+, and NK levels increased in both groups (P<0.05). CD8+ levels, BTA, NMP, and UBC levels decreased in both groups (P < 0.05). Analysis of survival results indicated that the two-year survival rates of patients with miR-143 and miR-92a high expressions were significantly higher than patients with low expressions (P < 0.05).

Conclusion

The efficacy of Xiaozheng decoction combined with bladder perfusion in treating postoperative patients with bladder cancer was significant, which could reduce the incidence of adverse reactions and postoperative recurrence rate, improve serum tumor marker levels, and enhance patients' immunity with a good prognosis.

LncRNA FAM13A-AS1 Promotes Renal Carcinoma Tumorigenesis Through Sponging miR-141-3p to Upregulate NEK6 Expression

Long non-coding RNAs are a diverse catalog of RNAs that have been implicated in various aspects of tumorigenesis. Emerging evidence indicates that they play crucial roles in tumor growth, disease progression, and drug resistance. However, the clinical significance of lncRNAs in tumor behavior prediction and disease prognosis as well as the underlying mechanism in renal cell carcinoma (RCC) remains elusive. By analyzing the gene expression profiles of 539 RCC patients from the TCGA cohort and 40 RCC patients from an independent cohort, we identified FAM13A-AS1, a poorly studied lncRNA, upregulated in RCC patients. Knockdown experiments revealed that FAM13A-AS1 promotes cell proliferation, migration, and invasion by interacting with miR-141-3p. FAM13A-AS1 regulates the expression of NEK6 by decoying miR-141-3p. In addition, there was a strong positive correlation between the expression of FAM13A-AS1 and NEK6 in RCC patients. In summary, our results demonstrate the oncogenic role of FAM13A-AS1 in RCC and suggest that it promotes tumorigenesis by upregulating the expression of NEK6 by competitively binding to miR-141-3p.

lncRNA USP30-AS1 sponges miR-765 and modulates the progression of colon cancer

Background

The incidence and mortality of colon cancer is increasing recently. It is necessary to identify effective biomarkers for the progression and prognosis of colon cancer. To assess the potential of lncRNA USP30-AS1 (USP30-AS1) in serving as the biomarker of colon cancer and unearth the underlying mechanism.

Methods

There were 123 colon cancer patients enrolled. The expression of USP30-AS1 was evaluated with PCR in tissue and cell samples. The clinical significance of USP30-AS1 was assessed with a series of statistical methods, while the CCK8 and Transwell assay were conducted to estimate its biological effect on the colon cancer cellular processes. In mechanism, the interaction of USP30-AS1 with miR-765 was evaluated with the dual-luciferase reporter assay.

Results

In colon cancer tissues, the USP30-AS1 downregulation and the miR-765 upregulation were observed, and there was a negative correlation between the USP30-AS1 expression level and the miR-765 expression level. The downregulation of USP30-AS1 related to the malignant progression and served as an adverse prognostic indicator of colon cancer. The overexpression of USP30-AS1 dramatically suppressed colon cancer cellular processes, which was alleviated by miR-765.

Conclusions

USP30-AS1 predicts the malignancy and prognosis of colon cancer patients. USP30-AS1 suppressed the progression of colon cancer through modulating miR-765.

Identification of Novel Key Genes and Pathways in Multiple Sclerosis Based on Weighted Gene Coexpression Network Analysis and Long Noncoding RNA-Associated Competing Endogenous RNA Network

Objective

Multiple sclerosis (MS) is an autoimmune disease of the central nervous system characterized by chronic inflammation and demyelination. This study is aimed at identifying crucial genes and molecular pathways involved in MS pathogenesis.

Methods

Raw data in GSE52139 were collected from the Gene Expression Omnibus. The top 50% expression variants were subjected to weighted gene coexpression network analysis (WGCNA), and the key module associated with MS occurrence was identified. A long noncoding RNA- (lncRNA-) associated competing endogenous RNA (ceRNA) network was constructed in the key module. The hub gene candidates were subsequently verified in an individual database.

Results

Of the 18 modules obtained, the cyan module was designated as the key module. The established ceRNA network was composed of seven lncRNAs, 45 mRNAs, and 21 microRNAs (miRNAs), and the FAM13A-AS1 was the lncRNA with the highest centrality. Functional assessments indicated that the genes in the cyan module primarily gathered in ribosome-related functional terms. Interestingly, the targeted mRNAs of the ceRNA network enriched in diverse categories. Moreover, highly expressed CYBRD1, GNG12, and SMAD1, which were identified as hub genes, may be associated with “valine leucine and isoleucine degradation,” “base excision repair,” and “fatty acid metabolism,” respectively, according to the results of single gene-based genomes and gene set enrichment analysis (GSEA).

Conclusions

Combined with the WGCNA and ceRNA network, our findings provide novel insights into the pathogenesis of MS. The hub genes discovered herein might also serve as novel biomarkers that correlate with the development and management of MS.

Identification of a novel immune-related long noncoding RNA signature to predict the prognosis of bladder cancer

Tumour immune regulation has attracted widespread attention, and long noncoding RNAs (lncRNAs) play an important role in this process. Therefore, we evaluated patient prognosis by exploring the relationship between bladder cancer (BLCA) and immune-related lncRNAs (IRlncRNAs). Transcriptome data and immune-related genes were analysed for coexpression, and then, the IRlncRNAs were analysed to determine the differentially expressed IRlncRNAs (DEIRlncRNAs) between normal and tumour samples in The Cancer Genome Atlas. The screened lncRNAs were pairwise paired and combined with clinical data, and finally, a signature was constructed by Lasso regression and Cox regression in 13 pairs of DEIRlncRNAs. According to the Akaike information criterion (AIC) values of the 1-year receiver operating characteristic curve, BLCA patients were stratified into high- or low-risk groups. The high-risk group had a worse prognosis. A comprehensive analysis showed that differences in risk scores were associated with the immune status of BLCA-infiltrated patients. The identified signature was correlated with the expression of immune checkpoint inhibitor-related molecules and sensitivity to chemotherapeutic drugs. We also identified three BLCA clusters with different immune statuses and prognoses that are also associated with immunotherapy response and drug sensitivity. In conclusion, we constructed a powerful predictive signature with high accuracy and validated its prognostic value.

A novel prognostic model for hepatocellular carcinoma based on 5 microRNAs related to vascular invasion

BACKGROUND

Hepatocellular carcinoma (HCC) is prevalent worldwide with a high mortality rate. Prognosis prediction is crucial for improving HCC patient outcomes, but effective tools are still lacking. Characteristics related to vascular invasion (VI), an important process involved in HCC recurrence and metastasis, may provide ideas on prognosis prediction.

METHODS

Tools, including R 4.0.3, Funrich version 3, Cytoscape 3.8.2, STRING 11.5, Venny 2.1.0, and GEPIA 2, were used to perform bioinformatic analyses. The VI-related microRNAs (miRNAs) were identified using Gene Expression Omnibus HCC miRNA dataset GSE67140, containing 81 samples of HCC with VI and 91 samples of HCC without VI. After further evaluated the identified miRNAs based on The Cancer Genome Atlas database, a prognostic model was constructed via Cox regression analysis. The miRNAs in this model were also verified in HCC patients. Moreover, a nomogram was developed by integrating risk score from the prognostic model with clinicopathological parameters. Finally, a potential miRNA-mRNA network related to VI was established through weighted gene co-expression network analysis of HCC mRNA dataset GSE20017, containing 40 samples of HCC with VI and 95 samples of HCC without VI.

RESULTS

A prognostic model of 5 VI-related miRNAs (hsa-miR-126-3p, hsa-miR-148a-3p, hsa-miR-15a-5p, hsa-miR-30a-5p, hsa-miR-199a-5p) was constructed. The area under receiver operating characteristic curve was 0.709 in predicting 5-year survival rate, with a sensitivity of 0.74 and a specificity of 0.63. The nomogram containing risk score could also predict prognosis. Moreover, a VI-related miRNA-mRNA network covering 4 miRNAs and 15 mRNAs was established.

CONCLUSION

The prognostic model and nomogram might be potential tools in HCC management, and the VI-related miRNA-mRNA network gave insights into how VI was developed.

Identification and Validation of an EMT-Related LncRNA Signature for HNSCC to Predict Survival and Immune Landscapes

Long noncoding RNAs (lncRNAs) are increasingly recognized as decisive factors in the progression of head and neck squamous cell carcinoma (HNSCC), and they participate in the epithelial–mesenchymal transformation (EMT) of HNSCC. LncRNAs are closely related to the prognosis of patients with HNSCC; thus, it is essential to identify EMT-related lncRNAs with prognostic value for HNSCC. The coexpression network of EMT-related lncRNAs was constructed using The Cancer Genome Atlas (TCGA). An EMT-related eight-lncRNA-based prognostic signature was constructed using LASSO Cox regression and Cox proportional hazards analyses. Univariate and multivariate analyses and stratified prognosis confirmed that the prognostic signature was an independent predictive factor. Subsequently, we performed immune cell infiltration analysis, gene set enrichment analysis (GSEA), and single-sample GSEA (ssGSEA) pathway enrichment analysis to uncover the potential molecular mechanisms of prognostic differences in the high- and low-risk groups. Next, we discussed the relationship between the prognostic signature and immune checkpoint-related genes, their TIDE scores, and the sensitivity of common chemotherapeutics. Finally, we further verified the expression differences in lncRNAs that were included in our signature via RT–qPCR in eighteen paired tissues. In summary, this prognostic signature provides powerful prognostic biomarkers for HNSCC and could serve as a predictor for the sensitivity of common chemotherapeutics and immunotherapy responses as well as providing a reference for further personalized treatment.

Bioinformatic Analyses of the Ferroptosis-Related lncRNAs Signature for Ovarian Cancer

Both ferroptosis and lncRNAs are significant for ovarian cancer (OC). Whereas, the study of ferroptosis-related lncRNAs (FRLs) still few in ovarian cancer. We first constructed an FRL-signature for patients with OC in the study. A total of 548 FRLs were identified for univariate Cox regression analysis, and 21 FRLs with significant prognosis were identified. The prognostic characteristics of nine FRLs was constructed and validated, showing opposite prognosis in two subgroups based on risk scores. The multivariate Cox regression analysis and nomogram further verified the prognostic value of the risk model. By calculating ferroptosis score through ssGSEA, we found that patients with higher risk scores exhibited higher ferroptosis scores, and high ferroptosis score was a risk factor. There were 40 microenvironment cells with significant differences in the two groups, and the difference of Stromal score between the two groups was statistically significant. Six immune checkpoint genes were expressed at different levels in the two groups. In addition, five m6A regulators (FMR1, HNRNPC, METTL16, METTL3, and METTL5) were higher expressed in the low-risk group. GSEA revealed that the risk model was associated with tumor-related pathways and immune-associated pathway. We compared the sensitivity of chemotherapy drugs between the two risk groups. We also explored the co-expression, ceRNA relation, cis and trans interaction of ferroptosis-related genes and lncRNAs, providing a new idea for the regulatory mechanisms of FRLs. Moreover, the nine FRLs were selected for detecting their expression levels in OC cells and tissues.

An Unfolded Protein Response Related Signature Could Robustly Predict Survival Outcomes and Closely Correlate With Response to Immunotherapy and Chemotherapy in Bladder Cancer

Background: The unfolded protein response (UPR) plays a significant role in maintaining protein hemostasis in tumor cells, which are crucial for tumor growth, invasion, and resistance to therapy. This study aimed to develop a UPR-related signature and explore its correlation with immunotherapy and chemotherapy in bladder cancer. Methods: The differentially expressed UPR-related genes were put into Lasso regression to screen out prognostic genes, which constituted the UPR signature, and were incorporated into multivariate Cox regression to generate risk scores. Subsequently, the predictive performance of this signature was estimated by receiver operating characteristic (ROC) curves. The CIBERSORTx, the maftool, and Gene set enrichment analysis (GSEA) were applied to explore infiltrated immune cells, tumor mutational burden (TMB), and enriched signaling pathways in both risk groups, respectively. Moreover, The Cancer Immunome Atlas (TCIA) and Genomics of Drug Sensitivity in Cancer (GDSC) databases were used to predict responses to chemotherapy and immunotherapy. Results: Twelve genes constituted the UPR-related signature. Patients with higher risk scores had worse overall survival (OS) in training and three validation sets. The UPR-related signature was closely correlated with clinicopathologic parameters and could serve as an independent prognostic factor. M0 macrophages showed a significantly infiltrated difference in both risk groups. TMB analysis showed that the risk score in the wild type and mutation type of FGFR3 was significantly different. GSEA indicated that the immune-, extracellular matrix-, replication and repair associated pathways belonged to the high risk group and metabolism-related signal pathways were enriched in the low risk group. Prediction of immunotherapy and chemotherapy revealed that patients in the high risk group might benefit from chemotherapy, but had a worse response to immunotherapy. Finally, we constructed a predictive model with age, stage, and UPR-related risk score, which had a robustly predictive performance and was validated in GEO datasets. Conclusion: We successfully constructed and validated a novel UPR-related signature in bladder cancer, which could robustly predict survival outcomes and closely correlate with the response to immunotherapy and chemotherapy in bladder cancer.

The new horizon of biomarker in melanoma patients: A study based on autophagy-related long non-coding RNA

Autophagy-related long non-coding RNAs (arlncRNAs) play a crucial role in the pathogenesis and development of the tumor. However, there is a lack of systematic analysis of arlncRNAs in melanoma patients.Melanoma data for analysis were obtained from The Cancer Genome Atlas (TCGA) database. By establishing a co-expression network of autophagy-related mRNAs-lncRNAs, we identified arlncRNAs in melanoma patients. We evaluated the prognostic value of arlncRNAs by univariate and multivariate Cox analysis and constructed an arlncRNAs risk model. Patients were divided into high- and low-risk groups based on the arlncRNAs risk score. This model was evaluated by Kaplan-Meier (K-M) analysis, univariate-multivariate Cox regression analysis, and receiver operating characteristic (ROC) curve analysis. Characteristics of autophagy genes and co-expressive tendency were analyzed by principal component analysis and Gene Set Enrichment Analysis (GSEA) functional annotation.Nine arlncRNAs (USP30-AS1, LINC00665, PCED1B-AS1, LINC00324, LINC01871, ZEB1-AS1, LINC01527, AC018553.1, and HLA-DQB1-AS1) were identified to be related to the prognosis of melanoma patients. Otherwise, the 9 arlncRNAs constituted an arlncRNAs prognostic risk model. K-M analysis and ROC curve analysis showed that the arlncRNAs risk model has good discrimination. Univariate and multivariate Cox regression analysis showed that arlncRNAs risk model was an independent prognostic factor in melanoma patients. Principal component analysis and GSEA functional annotation showed different autophagy and carcinogenic status in the high- and low-risk groups.This novel arlncRNAs risk model plays an essential role in predicting of the prognosis of melanoma patients. The model reveals new prognosis-related biomarkers for autophagy, promotes precision medicine, and provides a lurking target for melanoma's autophagy-related treatment.

Molecular Signatures Correlated With Poor IVF Outcomes: Insights From the mRNA and lncRNA Expression of Endometriotic Granulosa Cells

The outcomes of in vitro fertilization (IVF) for endometriotic women are significantly worse than for patients without ovarian endometriosis (OEM), as shown by fewer retrieved oocytes. However, the exact pathophysiological mechanism is still unknown. Thus, we conducted a prospective study that analyzed mRNA and lncRNA transcriptome between granulosa cells (GCs) from patients with fewer retrieved oocytes due to OEM and GCs from controls with male factor (MF) infertility using an RNA sequencing approach. We found a group of significantly differentially expressed genes (DEGs), including NR5A2, MAP3K5, PGRMC2, PRKAR2A, DEPTOR, ITGAV, KPNB1, GPC6, EIF3A, and SMC5, which were validated to be upregulated and negatively correlated with retrieved oocyte numbers in GCs of patients with OEM, while DUSP1 demonstrated the opposite. The molecular functions of these DEGs were mainly enriched in pathways involving mitogen-activated protein kinase (MAPK) signaling, Wnt signaling, steroid hormone response, apoptosis, and cell junction. Furthermore, we performed lncRNA analysis and identified a group of differentially expressed known/novel lncRNAs that were co-expressed with the validated DEGs and correlated with retrieved oocyte numbers. Co-expression networks were constructed between the DEGs and known/novel lncRNAs. These distinctive molecular signatures uncovered in this study are involved in the pathological regulation of ovarian reserve dysfunction in OEM patients.

LncRNA USP30-AS1 promotes the survival of acute myeloid leukemia cells by cis-regulating USP30 and ANKRD13A

Acute myeloid leukemia (AML) is a malignant tumor derived from leukemia stem cells, with complicated pathogenesis. LncRNAs play an important role in tumors genesis and progression. According to results from bioinformatics analysis, lncRNA USP30-AS1 is highly expressed in AML and both the high expression of USP30-AS1 and low methylation level at Cg03124318 locus of USP30-AS1 gene promoter are associated with poor prognosis of AML. This study knocked down and overexpressed USP30-AS1 to determine the roles in AML cell lines. High-throughput sequencing was performed to explore the genes regulated by USP30-AS1. Results showed that USP30-AS1 promoted AML cell viability and inhibited apoptosis. Genes regulated by USP30-AS1 are mainly related to genetic regulation and immune system. Among them, USP30 and ANKRD13A genes are close to USP30-AS1 gene in chromosome. Knockdown of USP30, but not ANKRD13A, abolished the cancer-promoting effects of USP30-AS1. ANKRD13A recognizes Lys-63-linked polyubiquitin chain in HLA-I. USP30-AS1 induced HLA-I internalization from the cell membrane by up-regulating ANKRD13A, which might induce the immune escape of AML cells. ChIP analysis revealed that the regulatory effects of USP30-AS1 on USP30 and ANKRD13A are associated with H3K4me3 and H3K27Ac. In summary, USP30-AS1 probably promotes AML cell survival by cis-regulating USP30 and ANKRD13A.

Immunotherapeutic Significance of a Prognostic Alternative Splicing Signature in Bladder Cancer

Objectives: Bladder cancer is the fourth most common malignancy in men in the United States. Aberrant alternative splicing (AS) events are involved in the carcinogenesis, but the association between AS and bladder cancer remains unclear. This study aimed to construct an AS-based prognostic signature and elucidate the role of the tumor immune microenvironment (TIME) and the response to immunotherapy and chemotherapy in bladder cancer. Methods: Univariate Cox regression analysis was performed to detect prognosis-related AS events. The least absolute shrinkage and selection operator (LASSO) and multivariate Cox analyses were employed to build prognostic signatures. Kaplan-Meier survival analysis, multivariate Cox regression analysis, and receiver operating characteristic (ROC) curves were conducted to validate the prognostic signatures. Then, the Estimation of Stromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) and tumor immune estimation resource (TIMER) databases were searched and the single-sample gene set enrichment analysis (ssGSEA) algorithm and CIBERSORT method were performed to uncover the context of TIME in bladder cancer. The Tumor Immune Dysfunction and Exclusion (TIDE) web tool and pRRophetic algorithm were used to predict the response to immunotherapy and chemotherapy. Finally, we constructed a correlation network between splicing factors (SFs) and survival-related AS events. Results: A total of 4684 AS events were significantly associated with overall survival in patients with bladder cancer. Eight prognostic signatures of bladder cancer were established, and a clinical survival prediction model was built. In addition, the consolidated prognostic signature was closely related to immune infiltration and the response to immunotherapy and chemotherapy. Furthermore, the correlation identified EIF3A, DDX21, SDE2, TNPO1, and RNF40 as hub SFs, and function analysis found ubiquitin-mediated proteolysis is correlated most significantly with survival-associated AS events. Conclusion: Our findings highlight the prognostic value of AS for patients with bladder cancer and reveal pivotal players of AS events in the context of TIME and the response to immunotherapy and chemotherapy, which may be important for patient management and treatment.

Recognition of Glycometabolism-Associated lncRNAs as Prognosis Markers for Bladder Cancer by an Innovative Prediction Model

The alteration of glycometabolism is a characteristic of cancer cells. Long non-coding RNAs (lncRNAs) have been documented to occupy a considerable position in glycometabolism regulation. This research aims to construct an effective prediction model for the prognosis of bladder cancer (BC) based on glycometabolism-associated lncRNAs (glyco-lncRNAs). Pearson correlation analysis was applied to get glyco-lncRNAs, and then, univariate cox regression analysis was employed to further filtrate survival time-associated glyco-lncRNAs. Multivariate cox regression analysis was utilized to construct the prediction model to divide bladder cancer (BC) patients into high- and low-risk groups. The overall survival (OS) rates of these two groups were analyzed using the Kaplan-Meier method. Next, gene set enrichment analysis and Cibersortx were used to explore the enrichment and the difference in immune cell infiltration, respectively. pRRophetic algorithm was applied to explore the relation between chemotherapy sensitivity and the prediction model. Furthermore, reverse transcriptase quantitative polymerase chain reaction was adopted to detect the lncRNAs constituting the prediction signature in tissues and urine exosomal samples of BC patients. A powerful model including 6 glyco-lncRNAs was proposed, capable of suggesting a risk score for each BC patient to predict prognosis. Patients with high-risk scores demonstrated a shorter survival time both in the training cohort and testing cohort, and the risk score could predict the prognosis without depending on the traditional clinical traits. The area under the receiver operating characteristic curve of the risk score was higher than that of other clinical traits (0.755 > 0.640, 0.485, 0.644, or 0.568). The high- and low-risk groups demonstrated very distinct immune cells infiltration conditions and gene set enriched terms. Besides, the high-risk group was more sensitive to cisplatin, docetaxel, and sunitinib. The expression of lncRNA AL354919.2 featured with an increase in low-grade patients and a decrease in T3-4 and Stage III-IV patients. Based on the experiment results, lncRNA AL355353.1, AC011468.1, and AL354919.2 were significantly upregulated in tumor tissues. This research furnishes a novel reference for predicting the prognosis of BC patients, assisting clinicians with help in the choice of treatment.

A risk model based on autophagy-related lncRNAs for predicting prognosis and efficacy of immunotherapy and chemotherapy in gastric cancer patients

Long non-coding RNAs (lncRNAs) are a class of non-protein-coding RNAs essential to the occurrence and development of gastric cancer (GC). We aimed to identify critical lncRNA pairs to construct a prognostic model and assess its performances in prognosis and efficacy prediction in GC patients receiving immunotherapy and chemotherapy. We searched transcriptome and clinical data of GC patients from The Cancer Genome Atlas (TCGA) database. Autophagy-related lncRNAs were identified using co-expression network analysis, and lncRNA pairs with prognostic value were selected using pairwise transcriptome analysis. The gene pairs were subjected to LASSO algorithm for identification of optimal gene pairs for risk model construction. Patients were classified into the low-risk and high-risk groups with the RiskScore as a cutoff. Finally, 9 optimal gene pairs were identified in the LASSO algorithm model for construction of a lncRNA prognostic risk model. For predictive performances, it successfully predicted a shorter survival of high-risk patients than that obtained in low-risk individuals (P < 0.001). It showed moderate AUC (area under the curve) values for 1-, 2-, and 3-year overall survival prediction of 0.713 and could serve as an independent predictor for GC prognosis. Compared to the low-risk group, high-risk patients had higher expressions of marker genes for immune checkpoint inhibitors (ICIs) and showed higher sensitivity to the chemotherapy agents, rapamycin, bexarotene, and bicalutamide. Our findings demonstrate a robust prognostic model based on nine autophagy-related lncRNA pairs for GC. It acts as an independent predictor for survival and efficacy prediction of immunotherapy and chemotherapy in GC patients.

A New Ferroptosis-Related lncRNA Signature Predicts the Prognosis of Bladder Cancer Patients

Background: Ferroptosis is closely related to the occurrence and development of cancer. An increasing number of studies have induced ferroptosis as a treatment strategy for cancer. However, the predictive value of ferroptosis-related lncRNAs in bladder cancer (BC) still need to be further elucidated. The purpose of this study was to construct a predictive signature based on ferroptosis-related long noncoding RNAs (lncRNAs) to predict the prognosis of BC patients.

Methods: We downloaded RNA-seq data and the corresponding clinical and prognostic data from The Cancer Genome Atlas (TCGA) database and performed univariate and multivariate Cox regression analyses to obtain ferroptosis-related lncRNAs to construct a predictive signature. The Kaplan-Meier method was used to analyze the overall survival (OS) rate of the high-risk and low-risk groups. Gene set enrichment analysis (GSEA) was performed to explore the functional differences between the high- and low-risk groups. Single-sample gene set enrichment analysis (ssGSEA) was used to explore the relationship between the predictive signature and immune status. Finally, the correlation between the predictive signature and the treatment response of BC patients was analyzed.

Results: We constructed a signature composed of nine ferroptosis-related lncRNAs (AL031775.1, AL162586.1, AC034236.2, LINC01004, OCIAD1-AS1, AL136084.3, AP003352.1, Z84484.1, AC022150.2). Compared with the low-risk group, the high-risk group had a worse prognosis. The ferroptosis-related lncRNA signature could independently predict the prognosis of patients with BC. Compared with clinicopathological variables, the ferroptosis-related lncRNA signature has a higher diagnostic efficiency, and the area under the receiver operating characteristic curve was 0.707. When patients were stratified according to different clinicopathological variables, the OS of patients in the high-risk group was shorter than that of those in the low-risk group. GSEA showed that tumor- and immune-related pathways were mainly enriched in the high-risk group. ssGSEA showed that the predictive signature was significantly related to the immune status of BC patients. High-risk patients were more sensitive to anti-PD-1/L1 immunotherapy and the conventional chemotherapy drugs sunitinib, paclitaxel, cisplatin, and docetaxel.

Conclusion: The predictive signature can independently predict the prognosis of BC patients, provides a basis for the mechanism of ferroptosis-related lncRNAs in BC and provides clinical treatment guidance for patients with BC.

USP30-AS1 contributes to mitochondrial quality control in glioblastoma cells

Glioblastoma is the most serious type of brain cancer with poor prognosis. Here, using the publicly available glioma database, we identified that USP30-AS1, an antisense lncRNA locating on the opposite strand of USP30 locus, is upregulated in human gliomas, particularly in high grade glioma. High level of USP30-AS1 is correlated with poor survival in both primary and recurrent glioma patients. USP30-AS1 regulates mitochondrial homeostasis and mitophagy in glioblastoma cells. Knockdown of USP30-AS1 decreases mitochondrial protein expression and mitochondrial mass, promotes mitochondrial uncoupler-induced mitophagy. However, USP30-AS1 does not regulate USP30 expression in a cis-regulatory manner. In summary, this study proposed that USP30-AS1 may serve as a valuable prognostic marker for gliomas. USP3-AS1 is a negative regulator of mitophagy and the regulatory effect is USP30-independent. USP30-AS1 mediated repression of mitophagy may contribute to the loss of mitochondrial homeostasis and tumor development in glioma.

N6-Methyladenosine-Related Long Noncoding RNAs as Potential Prognosis Biomarkers for Endometrial Cancer

Purpose

Endometrial cancer (EC) is a common gynaecologic malignancy with an increasing incidence rate and mortality in recent years. N6-methylandenosine (m6A)-related long noncoding RNA (lncRNA) plays a vital role in EC, emerging as one of the most abundant RNA modifications.

Materials and Methods

The Cancer Genome Atlas (TCGA) database and UCSC Xena were used to download data related to EC. Survival and univariate and multifactorial prognostic analyses were performed for m6A-related lncRNAs. The expression levels of the three lncRNAs were verified using q-PCR. A nomogram was used to create a clinical tool to assess overall survival. To investigate the relationship between m6A-related lncRNA and EC, we downloaded differential genes related to EC from the TCGA database and mined three m6A-related lncRNAs, namely SCARNA9, TRAF3IP2-AS1, and AL133243.2. The data were categorized into high- and low-risk groups based on m6A-associated lncRNA.

Results

Survival analysis revealed that the high-risk group had a lower survival rate. Survival analysis of three m6A-associated lncRNAs revealed that cases with high expression of SCARNA9 tended to have a poorer prognosis, whereas the opposite was true for TRAF3IP2-AS1, AL133243.2. Univariate and multifactorial prognostic analyses suggested statistical differences in patients’ age, FIGO stage, pathological grade, risk score, and prognosis of EC, which was confirmed by results of the separate prognostic factor analysis for the three lncRNAs. Risk status was validated as an independent prognostic indicator, and the prognostic nomogram combined patient age, pathological stage, and FIGO classification to assess 3–5-year survival. Cases from high- and low-risk groups were analysed for the tumour microenvironment and immune cell scores, and stromal cell scores were found to be lower in the high-risk group. Correlations were analysed using different databases for immune cell classification.

Conclusion

m6A-related lncRNAs may play a key role in the diagnosis and treatment of EC as targets of prognosis and the immune microenvironment.

Comprehensive analysis of the prognosis and biological significance for IFIT family in skin cutaneous melanoma

Interferon-induced protein with tetratricopeptide repeats (IFITs) genes, consisting of four members named IFIT1, IFIT2, IFIT3 and IFIT5, are involved in the progression of multiple cancer types, but their roles in skin cutaneous melanoma (SKCM) are still largely unknown. The TCGA-SKCM dataset, GSE15605 dataset and GSE100508 dataset were obtained in our study, and multiple online databases were used for data analysis and visualization, including GEPIA, GSCALite, MethSurv, DAVID, starBase and TIMER database. The mRNA expressing levels of all the four members included in IFIT family were elevated in SKCM tissues. In addition, ROC curve showed that the combined IFITs had a higher tumor prediction performance. Kaplan-Meier survival analysis revealed that the low expression of IFIT1/2/3/5 was associated with poor overall survival (OS) and disease-specific survival (DSS) in SKCM patients. Moreover, univariate and multivariate Cox regression analysis suggested that the low expression of IFIT2/3/5 was an independent risk factor for the prognosis of SKCM patients. Besides, cancer pathway activity analysis certified that the IFITs were involved in the apoptosis pathways, epithelial-mesenchymal transition (EMT) and cell cycle. Furthermore, drug sensitivity analysis indicated that the high expression of IFIT1/2/3 was sensitive to dasatinib drug. Additionally, the expressing levels of IFITs were found to be positively correlated with the level of immune cell infiltrates, immune biomarkers and m6A regulators. Finally, using bioinformatics analysis, we predicted that PAX8-AS1/Z83843.1-miR-92a-3p-IFIT2 axis might play crucial roles in the development and progression of SKCM. This study explored the prognostic values and biological significance of the IFITs in SKCM microenvironment. IFITs may serve as novel biomarkers for the diagnosis and prognosis of melanoma and potential immunotherapeutic targets.

Identification and validation of an autophagy-related signature for predicting survival in lower-grade glioma

Abnormal levels of autophagy have been implicated in the pathogenesis of multiple diseases, including cancer. However, little is known about the role of autophagy-related genes (ARGs) in low-grade gliomas (LGG). Accordingly, the aims of this study were to assess the prognostic values of ARGs and to establish a genetic signature for LGG prognosis. Expression profile data from patients with and without primary LGG were obtained from The Cancer Genome Atlas (TCGA) and Genome Tissue Expression databases, respectively, and consensus clustering was used to identify clusters of patients with distinct prognoses. Nineteen differentially expressed ARGs were selected with threshold values of FDR < 0.05 and |log2 fold change (FC)| ≥ 2, and functional analysis revealed that these genes were associated with autophagy processes as expected. An autophagy-related signature was established using a Cox regression model of six ARGs that separated patients from TCGA training cohort into high- and low-risk groups. Univariate and multivariate Cox regression analysis indicated that the signature-based risk score was an independent prognostic factor. The signature was successfully validated using the TCGA testing, TCGA entire, and Chinese Glioma Genome Atlas cohorts. Stratified analyses demonstrated that the signature was associated with clinical features and prognosis, and gene set enrichment analysis revealed that autophagy- and cancer-related pathways were more enriched in high-risk patients than in low-risk patients. The prognostic value and expression of the six signature-related genes were also investigated. Thus, the present study constructed and validated an autophagy-related prognostic signature that could optimize individualized survival prediction in LGG patients.