Identification and Validation of a Novel Immune-Related Four-lncRNA Signature for Lung Adenocarcinoma

Establishment of a Lactylation-Related Gene Signature for Hepatocellular Carcinoma Applying Bulk and Single-Cell RNA Sequencing Analysis

Background: Lactylation is closely involved in cancer progression, but its role in hepatocellular carcinoma (HCC) is unclear. The present work set out to develop a lactylation-related gene (LRG) signature for HCC. Methods: The lactylation score of tumor and normal groups was calculated using the gene set variation analysis (GSVA) package. The single-cell RNA sequencing (scRNA-seq) analysis of HCC was performed in the "Seurat" package. Prognostic LRGs were selected by performing univariate and least absolute shrinkage and selection operator (LASSO) Cox regression analyses to develop and validate a Riskscore model. Functional enrichment analysis was conducted by gene set enrichment analysis (GSEA) using the "clusterProfiler" package. Genomic characteristics between different risk groups were compared, and tumor mutational burden (TMB) was calculated by the "Maftools" package. Immune cell infiltration was assessed by algorithms of cell-type identification by estimating relative subsets of RNA transcript (CIBERSORT), microenvironment cell populations-counter (MCP-counter), estimating the proportions of immune and cancer cells (EPIC), tumor immune estimation resource (TIMER), and single-sample gene set enrichment analysis (ssGSEA). Immunotherapy response was predicted by the tumor immune dysfunction and exclusion (TIDE) algorithm. Drug sensitivity was analyzed using the "pRRophetic" package. A nomogram was established using the "rms" package. The expressions of the prognostic LRGs in HCC cells were verified by in vitro test, and cell counting kit-8 (CCK-8), wound healing, and transwell assays were carried out to measure the viability, migration, and invasion of HCC cells. Results: The lactylation score, which was higher in the tumor group than in the normal group, has been confirmed as an independent factor for the prognostic evaluation in HCC. Six prognostic LRGs, including two protective genes (FTCD and APCS) and four risk genes (LGALS3, C1orf43, TALDO1, and CCT5), were identified to develop a Riskscore model with a strong prognostic prediction performance in HCC. The scRNA-seq analysis revealed that LGALS3 was largely expressed in myeloid cells, while APCS, FTCD, TALDO1, CCT5, and C1orf43 were mainly expressed in hepatocytes. The high-risk group was primarily enriched in the pathways involved in tumor occurrence and development, with higher T cell infiltration. Moreover, the high-risk group was found to be less responsive to immunotherapy but was more sensitive to chemotherapeutic drugs. By integrating Riskscore and clinical features, a nomogram with a high predictive accuracy was developed. Additionally, C1orf43, CCT5, TALDO1, and LGALS3 were highly expressed in HCC cells. Silencing CCT5 inhibited the viability, migration, and invasion of HCC cells. Conclusion: The present work developed a novel LRG gene signature that could be considered a promising therapeutic target and biomarker for HCC.

Construction of an immunogenic cell death-related LncRNA signature to predict the prognosis of patients with lung adenocarcinoma

BACKGROUND

Lung adenocarcinoma (LUAD) is one of the most common malignant diseases worldwide. This study aimed to construct an immunogenic cell death (ICD)-related long non-coding RNA (lncRNA) signature to effectively predict the prognosis of LUAD.

METHODS

The RNA-sequencing and clinical data of LUAD were downloaded from The Cancer Genome Atlas (TCGA). Least absolute shrinkage and selection operator (LASSO) and stepwise multivariate Cox proportional hazard regression analysis were utilized to construct lncRNA signature. Then, the reliability of the signature was evaluated in the training, validation and whole cohorts. The differences in the immune landscape and drug sensitivity between the low- and high-risk groups were analyzed. Finally, the expression level of the selected ICD-related lncRNAs in LUAD cell lines via reverse transcription quantitative PCR (RT-qPCR). CCK-8 and transwell assays were performed to study biological function of AC245014.3.

RESULTS

A signature consisting of 5 ICD-related lncRNAs was constructed. Kaplan Meier (K-M) survival analysis showed shorter overall survival (OS) in high-risk group. The receiver operating characteristic (ROC) curves and Multivariate Cox regression analysis showed the signature was good predictive and independent prognostic factor in LUAD. Moreover, the high-risk group had a lower level of antitumor immunity and was less sensitive to some chemotherapeutics and targeted drugs. Finally, the expression level of selected ICD-related lncRNAs was validated in LUAD cell lines by RT-qPCR. Knockdown of AC245014.3 significantly suppressed LUAD proliferation, migration and invasion.

CONCLUSIONS

In this study, an ICD-related lncRNA signature was constructed, which could accurately predict the prognosis of LUAD patients and guide clinical treatment.

Screening of Immune-related lncRNAs in Lung Adenocarcinoma and Establishing a Survival Prognostic Risk Prediction Model

OBJECTIVE

This study aimed to improve lung adenocarcinoma (LUAD) prognosis prediction based on a signature of immune-related long non-coding RNAs (lncRNAs).

METHODS

LUAD samples from the TCGA database were divided into the immunity_H group and the immunity_L group. Differentially expressed RNAs (DERs) between the two groups were identified. Optimized immune-related lncRNAs combination was obtained using LASSO Cox regression. A prognostic risk prediction (RS) model was built and further validated in the training and validation datasets. A network among lncRNAs in the RS model, their co-expressed DERs, and the related KEGG pathways were established. Critical lncRNAs were validated in LUAD tissue samples.

RESULTS

In total, 255 DERs were obtained, and 11 immune-related lncRNAs were significantly related to prognosis. Six lncRNAs were demonstrated as an optimal combination for building the RS model, including LINC00944, LINC00930, LINC00607, LINC00582, LINC00543, and LINC00319. The KM curve and ROC curve revealed the RS model to be a reliable indicator for LUAD prognosis. LINC00944 and LINC00582 showed a co-expression relationship with the MS4A1. LINC00944, LINC00582, and MS4A1 were successfully validated in LUAD samples.

CONCLUSION

We have established a promising LUAD patient survival prediction model based on six immune-related lncRNAs. For LUAD patients, this prognostic model could guide personalized treatment.

Single-cell transcriptome analysis reveals the clinical implications of myeloid-derived suppressor cells in head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSC) is the most common malignant tumor that arises in the epithelium of the head and neck regions. Myeloid-derived suppressor cells (MDSCs) are one of the tumor-infiltrating immune cell populations, which play a powerful role in inhibiting anti-tumor immune response. Herein, we employed a single-cell RNA sequencing (scRNA-seq) dataset to dissect the heterogeneity of myeloid cells. We found that SPP1 ⁺ tumor-associated macrophages (TAMs) and MDSCs were the most abundant myeloid cells in the microenvironment. By cell cluster deconvolution from bulk RNA-seq datasets of larger patient groups, we observed that highly-infiltrated MDSC was a poor prognostic marker for patients' overall survival (OS) probabilities. To better apply the MDSC OS prediction values, we identified a set of six MDSC-related genes (ALDOA, CD52, FTH1, RTN4, SLC2A3, and TNFAIP6) as the prognostic signature. In both training and test cohorts, MDSC-related prognostic signature showed a promising value for predicting patients' prognosis outcomes. Further parsing the ligand-receptor pairs of intercellular communications by CellChat, we found that MDSCs could frequently interact with cytotoxic CD8 ⁺ T cells, SPP1 ⁺ TAMs, and endothelial cells. These interactions likely contributed to the establishment of an immunosuppressive microenvironment and the promotion of tumor angiogenesis. Our findings suggest that targeting MDSCs may serve as an alternative and promising target for the immunotherapy of HNSC.

Prognostic significance of lncRNA AP004608.1 in prostate cancer

This study aimed to screen and determine the value of AP004608.1 expression as a biomarker for Prostate cancer (PCa) survival. We investigated the expression and prognosis of AP004608.1 through bioinformatics analysis. Low AP004608.1 expression predicted favorable Overall survival (OS) and Progression-free survival (PFS) in PCa patients, according to the Cancer Genome Atlas (TCGA) database. Cox regression demonstrated that low AP004608.1 expression were in-dependent biomarkers for OS. Moreover, Gene Expression Omnibus (GEO) database was utilized to verify the prognostic role of AP004608.1 in PCa, and the similar results were reached. A meta-analysis revealed that low AP004608.1 expression was closely relevant to better OS. AP004608.1 could constitute a promising prognostic biomarker, and probably plays an important role in PCa.

Necroptosis-related lncRNA in lung adenocarcinoma: A comprehensive analysis based on a prognosis model and a competing endogenous RNA network

Background: Necroptosis, an innovative type of programmed cell death, involves the formation of necrosomes and eventually mediates necrosis. Multiple lines of evidence suggest that necroptosis plays a major role in the development of human cancer. However, the role of necroptosis in lung adenocarcinoma (LUAD) remains unclear. In this study, we aimed to construct an NRL-related prognostic model and comprehensively analyze the role of NRL in LUAD.

Methods: A necroptosis-related lncRNA (NRL) signature was constructed in the training cohort and verified in the validation and all cohorts based on The Cancer Genome Atlas database. In addition, a nomogram was developed. The tumor microenvironment (TME), checkpoint, human leukocyte antigen, and m6A methylation levels were compared between low-risk and high-risk groups. Then, we identified five truly prognostic lncRNAs (AC107021.2, AC027117.1, FAM30A, FAM83A-AS1, and MED4-AS1) and constructed a ceRNA network, and four hub genes of downstream genes were identified and analyzed using immune, pan-cancer, and survival analyses.

Results: The NRL signature could accurately predict the prognosis of patients with LUAD, and patients with low risk scores were identified with an obvious “hot” immune infiltration level, which was strongly associated with better prognosis. Based on the ceRNA network, we postulated that NRLs regulated the TME of patients with LUAD via cyclin-dependent kinase (CDK) family proteins.

Conclusion: We constructed an NRL signature and a ceRNA network in LUAD and found that NRLs may modulate the immune microenvironment of LUAD via CDK family proteins.

Construction of a prognostic immune-related lncRNA model and identification of the immune microenvironment in middle- or advanced-stage lung squamous carcinoma patients

Background

Methods

Results

Conclusion

Identification of a chromatin regulator signature and potential candidate drugs for bladder cancer

Background

Bladder cancer (BLCA) is a malignant tumor with a dismay outcome. Increasing evidence has confirmed that chromatin regulators (CRs) are involved in cancer progression. Therefore, we aimed to explore the function and prognostic value of CRs in BLCA patients.

Methods

Chromatin regulators (CRs) were acquired from the previous top research. The mRNA expression and clinical information were downloaded from TCGA and GEO datasets. Cox regression analysis and least absolute shrinkage and selection operator (LASSO) regression analysis were performed to select the prognostic gene and construct the risk model for predicting outcome in BLCA. The Kaplan-Meier analysis was used to assess the prognosis between high- and low-risk groups. We also investigated the drug sensitivity difference between high- and low-risk groups. CMAP dataset was performed to screen the small molecule drugs for treatment.

Results

We successfully constructed and validated an 11 CRs-based model for predicting the prognosis of patients with BLCA. Moreover, we also found 11 CRs-based model was an independent prognostic factor. Functional analysis suggested that CRs were mainly enriched in cancer-related signaling pathways. The CR-based model was also correlated with immune cells infiltration and immune checkpoint. Patients in the high-risk group were more sensitive to several drugs, such as mitomycin C, gemcitabine, cisplatin. Eight small molecule drugs could be beneficial to treatment for BLCA patients.

Conclusion:

In conclusion, our study provided novel insights into the function of CRs in BLCA. We identified a reliable prognostic biomarker for the survival of patients with BLCA.

Construction of a novel risk model based on the random forest algorithm to distinguish pancreatic cancers with different prognoses and immune microenvironment features

Immune-related long noncoding RNAs (irlncRNAs) are actively involved in regulating the immune status. This study aimed to establish a risk model of irlncRNAs and further investigate the roles of irlncRNAs in predicting prognosis and the immune landscape in pancreatic cancer. The transcriptome profiles and clinical information of 176 pancreatic cancer patients were retrieved from The Cancer Genome Atlas (TCGA). Immune-related genes (irgenes) downloaded from ImmPort were used to screen 1903 immune-related lncRNAs (irlncRNAs) using Pearson's correlation analysis (R > 0.5; p < 0.001). Random survival forest (RSF) and survival tree analysis showed that 9 irlncRNAs were highly correlated with overall survival (OS) according to the variable importance (VIMP) and minimal depth. Next, Cox regression analysis was used to establish a risk model with 3 irlncRNAs (LINC00462, LINC01887, RP11-706C16.8) that was evaluated by Kaplan-Meier analysis, the areas under the curve (AUCs) of the receiver operating characteristics and the C-index. Additionally, we performed Cox regression analysis to establish the clinical prognostic model, which showed that the risk score was an independent prognostic factor (p < 0.001). A nomogram and calibration plots were drawn to visualize the clinical features. The Wilcoxon signed-rank test and Pearson's correlation analysis further explored the irlncRNA signatures and immune cell infiltration, as well as the immunotherapy response.

Construction of a Ferroptosis-Related Long Non-coding RNA Prognostic Signature and Competing Endogenous RNA Network in Lung Adenocarcinoma

Ferroptosis-related genes play an important role in the progression of lung adenocarcinoma (LUAD). However, the potential function of ferroptosis-related lncRNAs in LUAD has not been fully elucidated. Thus, to explore the potential role of ferroptosis-related lncRNAs in LUAD, the transcriptome RNA-seq data and corresponding clinical data of LUAD were downloaded from the TCGA dataset. Pearson correlation was used to mine ferroptosis-related lncRNAs. Differential expression and univariate Cox analysis were performed to screen prognosis related lncRNAs. A ferroptosis-related lncRNA prognostic signature (FLPS), which included six ferroptosis-related lncRNAs, was constructed by the least absolute shrinkage and selection operator (LASSO) Cox regression. Patients were divided into a high risk-score group and low risk-score group by the median risk score. Receiver operating characteristic (ROC) curves, principal component analysis (PCA), and univariate and multivariate Cox regression were performed to confirm the validity of FLPS. Enrichment analysis showed that the biological processes, pathways and markers associated with malignant tumors were more common in high-risk subgroups. There were significant differences in immune microenvironment and immune cells between high- and low-risk groups. Then, a nomogram was constructed. We further investigated the relationship between six ferroptosis-related lncRNAs and tumor microenvironment and tumor stemness. A competing endogenous RNA (ceRNA) network was established based on the six ferroptosis-related lncRNAs. Finally, we detected the expression levels of ferroptosis-related lncRNAs in clinical samples through quantitative real-time polymerase chain reaction assay (qRT-PCR). In conclusion, we identified the prognostic ferroptosis-related lncRNAs in LUAD and constructed a prognostic signature which provided a new strategy for the evaluation and prediction of prognosis in LUAD.

Relationships of N6-Methyladenosine-Related Long Non-Coding RNAs With Tumor Immune Microenvironment and Clinical Prognosis in Lung Adenocarcinoma

Background: Lung adenocarcinoma (LUAD) is the major subtype of lung cancer and is associated with very high mortality. Emerging studies have shown that N6-methyladenosine (m6A)-related long non-coding (lnc) RNAs play crucial roles in tumor prognosis and the tumor immune microenvironment (TME). We aimed to explore the expression patterns of different m6A-related lncRNAs concerning patient prognosis and construct an m6A-related lncRNA prognostic model for LUAD. Methods: The prognostic value of m6A-related lncRNAs was investigated in LUAD samples from The Cancer Genome Atlas (TCGA). Potential prognostic m6A-related lncRNAs were selected by Pearson's correlation and univariate Cox regression analysis. Patients were divided into clusters using principal component analysis and the m6A-related lncRNA prognostic signature was calculated using least absolute shrinkage and selection operator (LASSO) Cox regression analysis. Results: Based on 91 prognostic m6A-related lncRNAs, we identified two m6A-related-lncRNA pattern clusters with different overall survival (OS) and different TMEs. We subsequently verified our findings multidimensionally by constructing a 13 m6A-related lncRNA prognostic signature (m6A-LPS) to calculate the risk score, which was robust in different subgroups. The receiver operating characteristic (ROC) curves and concordance index demonstrated that m6A-LPS harbored a promising ability to predict OS in TCGA data set and independent GSE11969 cohort. The risk score was also related to OS, TME, and clinical stage, and the risk score calculated by our model was also identified as independent prognostic predictive factors for LUAD patients after adjustment for age, smoking, gender, and stage. Enrichment analysis indicated that malignancy and drug resistance-associated pathways were more common in cluster2 (LUAD-unfavorable m6A-LPS). Furthermore, the results indicated that the signaling pathway enriched by the target gene of 13 m6A-related lncRNAs may be associated with metastasis and progression of cancer according to current studies. Conclusion: The current results indicated that different m6A-related-lncRNA patterns could affect OS and TME in patients with LUAD, and the prognostic signature based on 13 m6A-related lncRNAs may help to predict the prognosis in LUAD patients.

Characterization of Immune-Related Long Non-coding RNAs to Construct a Novel Signature and Predict the Prognosis and Immune Landscape of Soft Tissue Sarcoma

Background: Increasing evidence has demonstrated that immune-related long non-coding RNAs (irlncRNAs) are critically involved in tumor initiation and progression and associated with the prognosis of various cancers. However, their role in soft tissue sarcoma (STS) remains significantly uninvestigated.

Materials and Methods: Gene expression profiles were extracted from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) for the identification of irlncRNAs. Univariate analysis and modified least absolute shrinkage and selection operator (LASSO) penalized regression analysis were employed to determine differently expressed irlncRNA (DEirlncRNA) pairs of prognostic value, and subsequently, a risk signature based on DEirlncRNA pairs was established. Furthermore, Kaplan–Meier analysis and the area under the receiver operating characteristic curve (AUC) were used to assess survival differences and the predictive accuracy of the risk signature, respectively. Lastly, the correlation of irlncRNAs with immune characteristics and chemosensitivity in patients with STS were further investigated.

Results: A total of 1088 irlncRNAs were identified, and 311 irlncRNAs were distinguished as DEirlncRNAs. A total of 130 DEirlncRNA pairs were further identified as prognostic markers, and 14 pairs were selected for establishing a risk signature. The irlncRNA-based risk signature functioned as an independent prognostic marker for STS. Compared with the patients in the high-risk group, those in the low-risk group exhibited a better prognosis and were more sensitive to several chemotherapeutic agents. In addition, the irlncRNA-based risk signature was significantly associated with immune scores, infiltrating immune cells, and the expression of several immune checkpoints.

Conclusion: In conclusion, our data revealed that the irlncRNA-based risk signature resulted in reliable prognosis, effectively predicted the immune landscape of patients with STS and was significantly correlated with chemosensitivity, thus providing insights into the potential role of irlncRNAs as prognostic biomarkers and novel therapeutic targets for STS.

Construction of a Novel Immune-Related lncRNA Pair Signature with Prognostic Significance for Kidney Clear Cell Renal Cell Carcinoma

Background

Renal cell carcinoma (RCC) is one of the most common aggressive malignant tumors in the urinary system, among which the clear cell renal cell carcinoma (ccRCC) is the most common subtype. The immune-related long noncoding ribonucleic acids (irlncRNAs) which are abundant in immune cells and immune microenvironment (IME) have potential significance in evaluating the prognosis and effects of immunotherapy. The signature based on irlncRNA pairs and independent of the exact expression level seems to have a latent predictive significance for the prognosis of patients with malignant tumors but has not been applied in ccRCC yet.

Method

In this article, we retrieved The Cancer Genome Atlas (TCGA) database for the transcriptome profiling data of the ccRCC and performed coexpression analysis between known immune-related genes (ir-genes) and lncRNAs to find differently expressed irlncRNA (DEirlncRNA). Then, we adopted a single-factor test and a modified LASSO regression analysis to screen out ideal DEirlncRNAs and constructed a Cox proportional hazard model. We have sifted 28 DEirlncRNA pairs, 12 of which were included in this model. Next, we compared the area under the curve (AUC), found the cutoff point by using the Akaike information criterion (AIC) value, and distinguished the patients with ccRCC into a high-risk group and a low-risk group using this value. Finally, we tested this model by investigating the relationship between risk score and survival, clinical pathological characteristics, cells in tumor immune microenvironment, chemotherapy, and targeted checkpoint biomarkers.

Results

A novel immune-related lncRNA pair signature consisting of 12 DEirlncRNA pairs was successfully constructed and tightly associated with overall survival, clinical pathological characteristics, cells in tumor immune microenvironment, and reactiveness to immunotherapy and chemotherapy in patients with ccRCC. Besides, the efficacy of this signature was verified in some commonly used clinicopathological subgroups and could serve as an independent prognostic factor in patients with ccRCC.

Conclusions

This signature was proven to have a potential predictive significance for the prognosis of patients with ccRCC and the efficacy of immunotherapy.

Construction of a Novel Lung Adenocarcinoma Immune-Related lncRNA Pair Signature

Background

A growing number of studies have demonstrated that immune-related long noncoding ribonucleic acids (irlncRNAs) are potential prognostic factors for lung adenocarcinoma. Two-gene combination patterns could improve the sensitivity of prognostic models, providing us a novel signature construction concept that we applied to lung adenocarcinoma.

Methods

Gene expression and clinical data were downloaded from the Lung Adenocarcinoma project of The Cancer Genome Atlas (TCGA) database. We applied a co-expression analysis with immune genes obtained from the ImmPort database to recognize irlncRNA. The matrix of irlncRNA pairs was established by a cyclic comparison of each lncRNA pair expression level. Univariate and multivariate Cox regressions and Lasso penalized regression analysis were applied to construct the risk model. Patients with lung adenocarcinoma were divided into high- and low-risk groups, according to the Akaike Information Criterion (AIC) values of the receiver operating characteristic (ROC) curve. Then, we evaluated our signature under various clinical settings: clinical-pathological characteristics, tumor-infiltrating immune cells, checkpoint-related biomarkers, targeted therapy, and chemotherapy.

Results

Based on the 239 differently expressed irlncRNAs, we constructed an 11-irlncRNA pair signature. The area under the curve (AUC) of the ROC curve for the signature to predict the 4-year survival rate was 0.819, and the cut-off point was recognized as 1.003. Subsequent analysis showed that our signature can effectively distinguish unfavorable survival outcomes, prognostic clinic-pathological characteristics, and specify tumor infiltration status. Highly expressed immune checkpoint-related genes, as well as higher chemosensitivity, were correlated to the low-risk group.

Conclusion

We constructed a novel lung adenocarcinoma irlncRNA signature with promising prognostic value using the TCGA database, based on paired irlncRNAs and not relying on lncRNAs special expression levels.