Construction of an Immune-Related Six-lncRNA Signature to Predict the Outcomes, Immune Cell Infiltration, and Immunotherapy Response in Patients With Hepatocellular Carcinoma

Prognostic Value of an Immune Long Non-Coding RNA Signature in Liver Hepatocellular Carcinoma

In recent years, there has been a growing recognition of the important role that long non-coding RNAs (lncRNAs) play in the immunological process of hepatocellular carcinoma (LIHC). An increasing number of studies have shown that certain lncRNAs hold great potential as viable options for diagnosis and treatment in clinical practice. The primary objective of our investigation was to devise an immune lncRNA profile to explore the significance of immune-associated lncRNAs in the accurate diagnosis and prognosis of LIHC. Gene expression profiles of LIHC samples obtained from TCGA database were screened for immune-related genes. The optimal immune-related lncRNA signature was built via correlational analysis, univariate and multivariate Cox analysis. Then, the Kaplan-Meier plot, ROC curve, clinical analysis, gene set enrichment analysis, and principal component analysis were performed to evaluate the capability of the immune lncRNA signature as a prognostic indicator. Six long non-coding RNAs were identified via correlation analysis and Cox regression analysis considering their interactions with immune genes. Subsequently, tumor samples were categorized into two distinct risk groups based on different clinical outcomes. Stratification analysis indicated that the prognostic ability of this signature acted as an independent factor. The Kaplan-Meier method was employed to conduct survival analysis, results showed a significant difference between the two risk groups. The predictive performance of this signature was validated by principal component analysis (PCA). Additionally, data obtained from gene set enrichment analysis (GSEA) revealed several potential biological processes in which these biomarkers may be involved. To summarize, this study demonstrated that this six-lncRNA signature could be identified as a potential factor that can independently predict the prognosis of LIHC patients.

Long non-coding RNA NRAV in the 12q24.31 risk locus drives gastric cancer development through glucose metabolism reprogramming

Long non-coding RNAs (lncRNAs) serve as vital candidates to mediate cancer risk. Here, we aimed to identify the risk single-nucleotide polymorphisms (SNPs)-induced lncRNAs and to investigate their roles in gastric cancer (GC) development. Through integrating the differential expression analysis of lncRNAs in GC tissues and expression quantitative trait loci analysis in normal stomach tissues and GC tissues, as well as genetic association analysis based on GC genome-wide association studies and an independent validation study, we identified four lncRNA-related SNPs consistently associated with GC risk, including SNHG7 [odds ratio (OR) = 1.16, 95% confidence interval (CI): 1.09-1.23], NRAV (OR = 1.11, 95% CI: 1.05-1.17), LINC01082 (OR = 1.16, 95% CI: 1.08-1.22) and FENDRR (OR = 1.16, 95% CI: 1.07-1.25). We further found that a functional SNP rs6489786 at 12q24.31 increases binding of MEOX1 or MEOX2 at a distal enhancer and results in up-regulation of NRAV. The functional assays revealed that NRAV accelerates GC cell proliferation while inhibits GC cell apoptosis. Mechanistically, NRAV decreases the expression of key subunit genes through the electron transport chain, thereby driving the glucose metabolism reprogramming from aerobic respiration to glycolysis. These findings suggest that regulating lncRNA expression is a crucial mechanism for risk-associated variants in promoting GC development.

Screening immune-related blood biomarkers for DKD-related HCC using machine learning

Background

Diabetes mellitus is a significant health problem worldwide, often leading to diabetic kidney disease (DKD), which may also influence the occurrence of hepatocellular carcinoma (HCC). However, the relationship and diagnostic biomarkers between DKD and HCC are unclear.

Methods

Using public database data, we screened DKD secretory RNAs and HCC essential genes by limma and WGCNA. Potential mechanisms, drugs, and biomarkers for DKD-associated HCC were identified using PPI, functional enrichment, cMAP, and machine learning algorithms, and a diagnostic nomogram was constructed. Then, ROC, calibration, and decision curves were used to evaluate the diagnostic performance of the nomograms. In addition, immune cell infiltration in HCC was explored using CIBERSORT. Finally, the detectability of critical genes in blood was verified by qPCR.

Results

104 DEGs associated with HCC using WGCNA were identified. 101 DEGs from DKD were predicated on secreting into the bloodstream with Exorbase datasets. PPI analysis identified three critical modules considered causative genes for DKD-associated HCC, primarily involved in inflammation and immune regulation. Using lasso and RM, four hub genes associated with DKD-associated HCC were identified, and a diagnostic nomogram confirmed by DCA curves was established. The results of immune cell infiltration showed immune dysregulation in HCC, which was associated with the expression of four essential genes. PLVAP was validated by qPCR as a possible blood-based diagnostic marker for DKD-related HCC.

Conclusion

We revealed the inflammatory immune pathways of DKD-related HCC and developed a diagnostic nomogram for HCC based on PLVAP, C7, COL15A1, and MS4A6A. We confirmed with qPCR that PLVAP can be used as a blood marker to assess the risk of HCC in DKD patients.

Establishment of an endoplasmic reticulum stress-associated lncRNAs model to predict prognosis and immunological characteristics in hepatocellular carcinoma

BACKGROUND

The endoplasmic reticulum stress (ERS) and unfolded protein response (UPR) pathways play an essential role in the pathophysiology of hepatocellular carcinoma (HCC), and activation of the UPR pathway is strongly associated with tumor growth. However, the function of ERS-associated long non-coding RNAs (lncRNAs) in HCC is less recognized.

METHODS

We have used TCGA (The Cancer Genome Atlas) to obtain clinical and transcriptome data for HCC patients and the GSEA (Gene Set Enrichment Analysis) molecular signature database to get the ERS gene. ERS-associated prognostic lncRNA was determined using univariate Cox regression study. Then, least absolute shrinkage and selection operator and multivariate Cox regression study were used to construct ERS-associated lncRNAs risk model. Next, we use Kaplan-Meier (KM) survival study, time-dependent receiver operating characteristic (ROC) curve, univariate and multivariate Cox regression study to validate and evaluate the risk model. GSEA reveals the underlying molecular mechanism of the risk model. In addition, differences in Immune cell Infiltration Study, half-maximal inhibitory concentration (IC50) and immune checkpoints blockade (ICB) treatment between high and low risk groups were analyzed.

RESULTS

We constructed a risk model consisting of 6 ERS-associated lncRNAS (containingMKLN1-AS, LINC01224, AL590705.3, AC008622.2, AC145207.5, and AC026412.3). The KM survival study showed that the prognosis of HCC patients in low-risk group was better than that in high-risk group. ROC study, univariate and multivariate Cox regression study showed that the risk model had good predictive power for HCC patients. Our verification sample verified the aforesaid findings. GSEA suggests that several tumor- and metabolism-related signaling pathways are associated with risk groups. Simultaneously, we discovered that the risk models may help in the treatment of ICB and the selection of chemotherapeutic drugs.

CONCLUSIONS

In this article, we created an ERS-associated lncRNAs risk model to help prognostic diagnosis and personalized therapy in HCC.

Assessment of prognostic role of a novel 7-lncRNA signature in HCC patients

Background

Hepatocellular carcinoma (HCC) is characterized by extensive risk factors, high morbidity and mortality. Clinical prognostic evaluation assay assumes a nonspecific quality. Better HCC prognostics are urgently needed. Long noncoding RNAs (lncRNAs) exerts a crucial role in tumorigenesis and development. Excavating specific lncRNAs signature to ameliorate the high-risk survival prediction in HCC patients is worthwhile.

Methods

Differentially expressed lncRNAs (DElncRNAs) profile was acquired from The Cancer Genome Atlas database (TCGA). Then, the lncRNAs high-risk survival prognostic model was established using the least absolute shrinkage and selection operator (LASSO)-Cox regression algorithm. The lncRNAs were evaluated in clinical specimen by PCR. The receiver operating characteristic curve (ROC) analysis was further conducted to assess the potential prognostic value of the model. Moreover, a visible nomogram containing clinicopathological features and prognostic model was developed for prediction of survival property. Potential molecular mechanism was assessed by GO, KEGG, GSEA enrichment analysis and CIBERSORT immune infiltration analysis.

Results

A novel 7-lncRNA risk model (AL161937.2, LINC01063, AC145207.5, POLH-AS1, LNCSRLR, MKLN1-AS, AC105345.1) was constructed and validated for HCC prognosis prediction. Kaplan-Meier analysis revealed that patients in the high-risk group suffered a poor prognosis (p = 1.813 × 10^-8). These genes were detected by PCR, and the expression trend was in accordance with TCGA database. Interestingly, the risk score served as an independent risk factor for HCC patients (HR: 1.166, 95% CI:1.119-1.214, p < 0.001). The nomogram was established, and the predictive accuracy in the nomogram was prior to the TNM stage according to the ROC curve analysis. Cell proliferation related pathway, decreased CD4⁺ T cell, CD8⁺ T cell, NK cell and elevated Neutrophil, Macrophage M0 were observed in high-risk group. Besides, suppression of MKLN1-AS expression inhibited cell proliferation of HCC cells by CCK8 assay in vitro.

Conclusion

The 7-lncRNA signature may exert a particular prognostic prediction role in HCC and provide new insight in HCC carcinogenesis.

Integrated Analysis of N1-Methyladenosine Methylation Regulators-Related lncRNAs in Hepatocellular Carcinoma

N1-methyladenosine (m1A) and long non-coding RNAs (lncRNAs) play significant roles in tumor progression in hepatocellular carcinoma (HCC). However, their association with HCC is still unclear. In this study, lncRNAs related to m1A were extracted from the mRNA expression matrix in The Cancer Genome Atlas (TCGA) database. Five m1A-related lncRNAs (AL031985.3, NRAV, WAC-AS1, AC026412.3, and AC099850.4) were identified based on lasso Cox regression and they generated a prognostic signature of HCC. The prognostic signature was identified as an independent prognosis factor in HCC patients. Moreover, the prognostic signature achieved better performance than TP53 mutation status or tumor mutational burden (TMB) scores in the stratification of patient survival. The immune landscape indicated that most immune checkpoint genes and immune cells were distributed differently between both risk groups. A higher IC50 of chemotherapeutics (sorafenib, nilotinib, sunitinib, and gefitinib) was observed in the high-risk group, and a lower IC50 of gemcitabine in the low-risk group, suggesting the potential of the prognostic signature in chemosensitivity. In addition, fifty-five potential small molecular drugs were found based on drug sensitivity and NRAV expression. Together, five m1A-related lncRNAs generated a prognostic signature that could be a promising prognostic prediction approach and therapeutic response assessment tool for HCC patients.

SnoRNA and lncSNHG: Advances of nucleolar small RNA host gene transcripts in anti-tumor immunity

The small nucleolar RNA host genes (SNHGs) are a group of genes that can be transcript into long non-coding RNA SNHG (lncSNHG) and further processed into small nucleolar RNAs (snoRNAs). Although lncSNHGs and snoRNAs are well established to play pivotal roles in tumorigenesis, how lncSNHGs and snoRNAs regulate the immune cell behavior and function to mediate anti-tumor immunity remains further illustrated. Certain immune cell types carry out distinct roles to participate in each step of tumorigenesis. It is particularly important to understand how lncSNHGs and snoRNAs regulate the immune cell function to manipulate anti-tumor immunity. Here, we discuss the expression, mechanism of action, and potential clinical relevance of lncSNHGs and snoRNAs in regulating different types of immune cells that are closely related to anti-tumor immunity. By uncovering the changes and roles of lncSNHGs and snoRNAs in different immune cells, we aim to provide a better understanding of how the transcripts of SNHGs participate in tumorigenesis from an immune perspective.

A passage-dependent network for estimating the in vitro senescence of mesenchymal stromal/stem cells using microarray, bulk and single cell RNA sequencing

Long-term in vitro culture of human mesenchymal stem cells (MSCs) leads to cell lifespan shortening and growth stagnation due to cell senescence. Here, using sequencing data generated in the public domain, we have established a specific regulatory network of "transcription factor (TF)-microRNA (miRNA)-Target" to provide key molecules for evaluating the passage-dependent replicative senescence of mesenchymal stem cells for the quality control and status evaluation of mesenchymal stem cells prepared by different procedures. Short time-series expression miner (STEM) analysis was performed on the RNA-seq and miRNA-seq databases of mesenchymal stem cells from various passages to reveal the dynamic passage-related changes of miRNAs and mRNAs. Potential miRNA targets were predicted using seven miRNA target prediction databases, including TargetScan, miRTarBase, miRDB, miRWalk, RNA22, RNAinter, and TargetMiner. Then use the TransmiR v2.0 database to obtain experimental-supported transcription factor for regulating the selected miRNA. More than ten sequencing data related to mesenchymal stem cells or mesenchymal stem cells reprogramming were used to validate key miRNAs and mRNAs. And gene set variation analysis (GSVA) was performed to calculate the passage-dependent signature. The results showed that during the passage of mesenchymal stem cells, a total of 29 miRNAs were gradually downregulated and 210 mRNA were gradually upregulated. Enrichment analysis showed that the 29 miRNAs acted as multipotent regulatory factors of stem cells and participated in a variety of signaling pathways, including TGF-beta, HIPPO and oxygen related pathways. 210 mRNAs were involved in cell senescence. According to the target prediction results, the targets of these key miRNAs and mRNAs intersect to form a regulatory network of "TF-miRNA-Target" related to replicative senescence of cultured mesenchymal stem cells, across 35 transcription factor, 7 miRNAs (has-mir-454-3p, has-mir-196b-5p, has-mir-130b-5p, has-mir-1271-5p, has-let-7i-5p, has-let-7a-5p, and has-let-7b-5p) and 7 predicted targets (PRUNE2, DIO2, CPA4, PRKAA2, DMD, DDAH1, and GATA6). This network was further validated by analyzing datasets from a variety of mesenchymal stem cells subculture and lineage reprogramming studies, as well as qPCR analysis of early passages mesenchymal stem cells versus mesenchymal stem cells with senescence morphologies (SA-β-Gal⁺). The "TF-miRNA-Target" regulatory network constructed in this study reveals the functional mechanism of miRNAs in promoting the senescence of MSCs during in vitro expansion and provides indicators for monitoring the quality of functional mesenchymal stem cells during the preparation and clinical application.

Revealing Prognostic and Immunotherapy-Sensitive Characteristics of a Novel Cuproptosis-Related LncRNA Model in Hepatocellular Carcinoma Patients by Genomic Analysis

Immunotherapy has shown strong anti-tumor activity in a subset of patients. However, many patients do not benefit from the treatment, and there is no effective method to identify sensitive immunotherapy patients. Cuproptosis as a non-apoptotic programmed cell death caused by excess copper, whether it is related to tumor immunity has attracted our attention. In the study, we constructed the prognostic model of 9 cuproptosis-related LncRNAs (crLncRNAs) and assessed its predictive capability, preliminarily explored the potential mechanism causing treatment sensitivity difference between the high-/low-risk group. Our results revealed that the risk score was more effective than traditional clinical features in predicting the survival of HCC patients (AUC = 0.828). The low-risk group had more infiltration of immune cells (B cells, CD8+ T cells, CD4+ T cells), mainly with anti-tumor immune function (p < 0.05). It showed higher sensitivity to immune checkpoint inhibitors (ICIs) treatment (p < 0.001) which may exert the effect through the AL365361.1/hsa-miR-17-5p/NLRP3 axis. In addition, NLRP3 mutation-sensitive drugs (VNLG/124, sunitinib, linifanib) may have better clinical benefits in the high-risk group. All in all, the crLncRNAs model has excellent specificity and sensitivity, which can be used for classifying the therapy-sensitive population and predicting the prognosis of HCC patients.

Whole transcriptome analysis reveals that immune infiltration- lncRNAs are related to cellular apoptosis in liver transplantation

Introduction

In most instances, liver transplantation (LT) is the only available treatment for end-stage liver diseases. However, LT could also induce serious liver diseases or injury, and the underlying mechanisms of LT-induced complications remain largely unknown, especially the mechanisms of the dysfunction of the immune system mediated by long noncoding RNAs (lncRNAs).

Methods

In this study, we globally analyzed the proportion of immune cells by using the transcriptome sequencing data (RNA-seq) of needle-core liver biopsies from pre- and post-transplantation recipients. Dysregulated lncRNAs were found to be correlated with the altered fractions of immune cells. We finally explored the potential targets of dysregulated lncRNAs and analyzed their functions in LT.

Results

We found that in the samples, some immune cells changed significantly after LT, including CD4 T cells, NK cells and mast cells. The proportion of macrophages in different polarization states also changed significantly, with M0 macrophages increasing and M2 macrophages decreasing. Through weighted gene co-expression network analysis (WGCNA), 7 gene expression modules related to LT were identified. These modules were related to changes in the proportion of different immune cells. The functions of these modules represent the response modes of different functional genes after LT. Among these modules, MEtan and MEyellow modules were primarily enriched in apoptosis and inflammatory pathways. Twelve immunity-related lncRNAs were identified for the first time, and the regulatory network co-changing with immune cells was also identified. The co-expressed genes of these lncRNAs were highly enriched in apoptosis-related pathways. Many apoptosis-related genes were found to be up-regulated after LT.

Discussion

In summary, we speculated that the expression and regulation of these apoptotic genes may be related to the changes in the proportion of immune cells. Some of these lncRNAs and apoptosis-related genes have been reported to be related to cell proliferation and apoptosis. They are also potential biomarkers or therapeutic targets.

Hsa_circ_0015278 Regulates FLT3-ITD AML Progression via Ferroptosis-Related Genes

AML with the FLT3-ITD mutation seriously threatens human health. The mechanism by which circRNAs regulate the pathogenesis of FLT3-ITD mutant-type AML through ferroptosis-related genes (FerRGs) remains unclear. Differentially expressed circRNAs and mRNAs were identified from multiple integrated data sources. The target miRNAs and mRNAs of the circRNAs were predicted using various databases. The PPI network, ceRNA regulatory network, GO, and KEGG enrichment analyses were performed. The "survival" and the "pROC" R packages were used for K-M and ROC analysis, respectively. GSEA, immune infiltration analysis, and clinical subgroup analysis were performed. Finally, circRNAs were validated by Sanger sequencing and qRT-PCR. In our study, 77 DECircs-1 and 690 DECircs-2 were identified. Subsequently, 11 co-up-regulated DECircs were obtained by intersecting DECircs-1 and DECircs-2. The target miRNAs of the circRNAs were screened by CircInteractome, circbank, and circAtlas. Utilizing TargetScan, ENCORI, and miRWalk, the target mRNAs of the miRNAs were uncovered. Ultimately, 73 FerRGs were obtained, and the ceRNA regulatory network was constructed. Furthermore, MAPK3 and CD44 were significantly associated with prognosis. qRT-PCR results confirmed that has_circ_0015278 was significantly overexpressed in FLT3-ITD mutant-type AML. In summary, we constructed the hsa_circ_0015278/miRNAs/FerRGs signaling axis, which provides new insight into the pathogenesis and therapeutic targets of AML with FLT3-ITD mutation.

Prognostic analysis of RAS-related lncRNAs in liver hepatocellular carcinoma

Background

Liver hepatocellular carcinoma (LIHC), whose incidence is increasing globally, is one of the most prevalent malignant cancers. RAS-related pathways are involved in the cell proliferation, migration, apoptosis, and metabolism in LIHC. Long noncoding RNAs (lncRNAs) also play important roles in the progression and prognosis of LIHC. However, the clinical role, prognostic significance, and immune regulation of RAS-related lncRNAs in LIHC remains unclear. Our study aims to construct and validate a RAS-related lncRNA prognostic risk signature that can estimate the prognosis and response to immunotherapy in LIHC.

Methods

The clinical information and corresponding messenger RNA (mRNA)/lncRNA expression profiles were obtained from The Cancer Genome Atlas (TCGA) database, and 502 RAS-related lncRNAs were identified by Pearson correlation analysis. A prognostic risk signature with 5 RAS-related lncRNAs was then developed based on the Cox regression and least absolute shrinkage and selection operator (LASSO) algorithm analyses. Subsequently, Kaplan-Meier survival curve, receiver operating characteristic (ROC) curve, and the nomogram were established to evaluate the predictive accuracy of the signature. In addition, the immune microenvironment, tumor mutation burden, and drug sensitivity associated with the signature were also analyzed in LIHC.

Results

Compared with the low-risk groups, the high-risk groups had an unfavorable outcome. Multivariate regression analysis revealed that the risk score signature was the independent prognostic factor superior to the other clinical variables. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses demonstrated that the risk score was highly associated with the nuclear division, DNA replication, and immune response. The group with high risk tended to hold a lower immune escape rate and better immunotherapy efficacy, while the group with low risk was more sensitive to some small molecular targeted drugs.

Conclusions

We developed a RAS-related lncRNA risk signature that was highly associated with the prognosis and response to immunotherapy and targeted drugs and which provided novel mechanistic insights into the personalized treatment and potential drug selection for patients with LIHC.

Development of a novel immune-related lncRNA prognostic signature for patients with hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is the most common neoplasm and the major cause of cancer-associated death worldwide. The high mortality rate of HCC is mainly attributed to its widespread prevalence and the lack of effective treatment. Immunotherapy as a promising, innovative approach has revolutionised the treatment of solid tumours. However, owing to the heterogeneity and complex tumour microenvironment of HCC, an efficient biomarker for immunotherapy has yet to be identified. We investigated the role of immune-related long non-coding RNAs (lncRNAs) as prognostic biomarkers in patients with HCC from The Cancer Genome Atlas (TCGA) database. Spearman correlation, univariate and multivariate Cox, and lasso regression analyses were utilised to screen lncRNAs associated with prognosis. Four lncRNAs were filtered out to develop an immune-associated lncRNA prognostic signature in TCGA training as well as validation cohorts. Patients with HCC were then categorised into low- and high-risk groups according to the median value of the risk scores to evaluate the ability of the prognostic model between training and validation cohorts. A nomogram (based on risk score and stage) was constructed to appraise the general overall survival (OS) of patients with HCC. Differences in immune cell infiltration, immune checkpoint inhibitor (ICI) treatment response, gene mutation, and drug sensitivity were observed between the two groups. Thus, the lncRNA prognostic signature can serve as a sensitive prognostic biomarker with potential in individualised immunotherapy for HCC patients.

Identification of senescence-associated long non-coding RNAs to predict prognosis and immune microenvironment in patients with hepatocellular carcinoma

Background: Cellular senescence plays a complicated and vital role in cancer development because of its divergent effects on tumorigenicity. However, the long non-coding RNAs (lncRNAs) associated with tumor senescence and their prognostic value in hepatocellular carcinoma (HCC) remain unexplored. Methods: The trans-cancer oncogene-induced senescence (OIS) signature was determined by gene set variation analysis (GSVA) in the cancer genome atlas (TCGA) dataset. The OIS-related lncRNAs were identified by correlation analyses. Cox regression analyses were used to screen lncRNAs associated with prognosis, and an optimal predictive model was created by regression analysis of the least absolute shrinkage and selection operator (LASSO). The performance of the model was evaluated by Kaplan-Meier survival analyses, nomograms, stratified survival analyses, and receiver operating characteristic curve (ROC) analyses. Gene set enrichment analysis (GSEA) and cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT) were carried out to explore the functional relevance and immune cell infiltration, respectively. Results: Firstly, we examined the pan-cancer OIS signature, and found several types of cancer with OIS strongly associated with the survival of patients, including HCC. Subsequently, based on the OIS signature, we identified 76 OIS-related lncRNAs with prognostic values in HCC. We then established an optimal prognostic model based on 11 (including NRAV, AC015908.3, MIR100HG, AL365203.2, AC009005.1, SNHG3, LINC01138, AC090192.2, AC008622.2, AL139423.1, and AC026356.1) of these lncRNAs by LASSO-Cox regression analysis. It was then confirmed that the risk score was an independent and potential risk indicator for overall survival (OS) (HR [95% CI] = 4.90 [2.74-8.70], p < 0.001), which outperforms those traditional clinicopathological factors. Furthermore, patients with higher risk scores also showed more advanced levels of a proinflammatory senescence-associated secretory phenotype (SASP), higher infiltration of regulatory T (Treg) cells and lower infiltration of naïve B cells, suggesting the regulatory effects of OIS on immune microenvironment. Additionally, we identified NRAV as a representative OIS-related lncRNA, which is over-expressed in HCC tumors mainly driven by DNA hypomethylation. Conclusion: Based on 11 OIS-related lncRNAs, we established a promising prognostic predictor for HCC patients, and highlighted the potential immune microenvironment-modulatory roles of OIS in HCC, providing a broad molecular perspective of tumor senescence.

Identification of Prognostic Fatty Acid Metabolism lncRNAs and Potential Molecular Targeting Drugs in Uveal Melanoma

Background

The aim of this study was to identify prognostic fatty acid metabolism lncRNAs and potential molecular targeting drugs in uveal melanoma through integrated bioinformatics analysis.

Methods

In the present study, we obtained the expression matrix of 309 FAM-mRNAs and identified 225 FAM-lncRNAs by coexpression network analysis. We then performed univariate Cox analysis, LASSO regression analysis, and cross-validation and finally obtained an optimized UVM prognosis prediction model composed of four PFAM-lncRNAs (AC104129.1, SOS1-IT1, IDI2-AS1, and DLGAP1-AS2).

Results

The survival curves showed that the survival time of UVM patients in the high-risk group was significantly lower than that in the low-risk group in the train cohort, test cohort, and all patients in the prognostic prediction model (P < 0.05). We further performed risk prognostic assessment, and the results showed that the risk scores of the high-risk group in the train cohort, test cohort, and all patients were significantly higher than those of the low-risk group (P < 0.05), patient survival decreased and the number of deaths increased with increasing risk scores, and AC104129.1, SOS1-IT1, and DLGAP1-AS2 were high-risk PFAM-lncRNAs, while IDI2-AS1 were low-risk PFAM-lncRNAs. Afterwards, we further verified the accuracy and the prognostic value of our model in predicting prognosis by PCA analysis and ROC curves.

Conclusion

We identified 24 potential molecularly targeted drugs with significant sensitivity differences between high- and low-risk UVM patients, of which 13 may be potential targeted drugs for high-risk patients. Our findings have important implications for early prediction and early clinical intervention in high-risk UVM patients.

Identification and validation of a prognostic model of necroptosis-related lncRNAs in hepatocellular carcinoma

Background: The study focused on establishing a prognostic survival model with six necroptosis-related lncRNAs to predict overall survival (OS) in patients with hepatocellular carcinoma (HCC).

Methods: The data of gene expression and clinical information of HCC patients were obtained from The Cancer Genome Atlas (TCGA). Cox regression with LASSO was used for constructing a necroptosis-related lncRNA survival model, which we further validated with qRT-PCR in vitro. The relative bioinformatics analysis and consensus cluster analysis were performed based on six differentially expressed lncRNAs.

Results: The survival prognostic model was constructed by using data from TCGA. Receiver operating characteristic (ROC) curves showed a good survival prediction by this model. GSEA showed that several signaling pathways were related to HCC progression. Immune-related functional analysis showed that aDCs, macrophages, Th2 cells, and Tregs have stronger correlation with the high-risk group. The consensus cluster analysis further validated the 6-lncRNA prognostic model.

Conclusion: A novel 6-lncRNA (AL606489.1, NRAV, LINC02870, DUXAP8, “ZFPM2-AS1,” and AL031985.3) prognostic model had an accurately predictive power in HCC prognosis, which might be worthy of clinical application.

Non-Invasive Biomarkers for Immunotherapy in Patients with Hepatocellular Carcinoma: Current Knowledge and Future Perspectives

Simple Summary

The search for non-invasive biomarkers is a hot topic in modern oncology, since a tissue biopsy has significant limitations in terms of cost and invasiveness. The treatment perspectives have been significantly improved after the approval of immunotherapy for patients with hepatocellular carcinoma; therefore, the quick identification of responders is crucial to define the best therapeutic strategy. In this review, the current knowledge on the available non-invasive biomarkers of the response to immunotherapy is described.

Abstract

The treatment perspectives of advanced hepatocellular carcinoma (HCC) have deeply changed after the introduction of immunotherapy. The results in responders show improved survival compared with Sorafenib, but only one-third of patients achieve a significant benefit from treatment. As the tumor microenvironment exerts a central role in shaping the response to immunotherapy, the future goal of HCC treatment should be to identify a proxy of the hepatic tissue condition that is easy to use in clinical practice. Therefore, the search for biomarkers that are accurate in predicting prognosis will be the hot topic in the therapeutic management of HCC in the near future. Understanding the mechanisms of resistance to immunotherapy may expand the patient population that will benefit from it, and help researchers to find new combination regimens to improve patients’ outcomes. In this review, we describe the current knowledge on the prognostic non-invasive biomarkers related to treatment with immune checkpoint inhibitors, focusing on serological markers and gut microbiota.

A novel LUAD prognosis prediction model based on immune checkpoint-related lncRNAs

Lung adenocarcinoma (LUAD) is a malignant disease with an extremely poor prognosis, and there is currently a lack of clinical methods for early diagnosis and precise treatment and management. With the deepening of tumor research, more and more attention has been paid to the role of immune checkpoints (ICP) and long non-coding RNAs (lncRNAs) regulation in tumor development. Therefore, this study downloaded LUAD patient data from the TCGA database, and finally screened 14 key ICP-related lncRNAs based on ICP-related genes using univariate/multivariate COX regression analysis and LASSO regression analysis to construct a risk prediction model and corresponding nomogram. After multi-dimensional testing of the model, the model showed good prognostic prediction ability. In addition, to further elucidate how ICP plays a role in LUAD, we jointly analyzed the immune microenvironmental changes in LAUD patients and performed a functional enrichment analysis. Furthermore, to enhance the clinical significance of this study, we performed a sensitivity analysis of common antitumor drugs. All the above works aim to point to new directions for the treatment of LUAD.

Development and validation of a novel cellular senescence-related prognostic signature for predicting the survival and immune landscape in hepatocellular carcinoma

Background: Cellular senescence is a typical irreversible form of life stagnation, and recent studies have suggested that long non-coding ribonucleic acids (lncRNA) regulate the occurrence and development of various tumors. In the present study, we attempted to construct a novel signature for predicting the survival of patients with hepatocellular carcinoma (HCC) and the associated immune landscape based on senescence-related (sr) lncRNAs.

Method: Expression profiles of srlncRNAs in 424 patients with HCC were retrieved from The Cancer Genome Atlas database. Lasso and Cox regression analyses were performed to identify differentially expressed lncRNAs related to senescence. The prediction efficiency of the signature was checked using a receiver operating characteristic (ROC) curve, Kaplan–Meier analysis, Cox regression analyses, nomogram, and calibration. The risk groups of the gene set enrichment analysis, immune analysis, and prediction of the half-maximal inhibitory concentration (IC50) were also analyzed. Quantitative real-time polymerase chain reaction (qPCR) was used to confirm the levels of AC026412.3, AL451069.3, and AL031985.3 in normal hepatic and HCC cell lines.

Results: We identified 3 srlncRNAs (AC026412.3, AL451069.3, and AL031985.3) and constructed a new risk model. The results of the ROC curve and Kaplan–Meier analysis suggested that it was concordant with the prediction. Furthermore, a nomogram model was constructed to accurately predict patient prognosis. The risk score also correlated with immune cell infiltration status, immune checkpoint expression, and chemosensitivity. The results of qPCR revealed that AC026412.3 and AL451069.3 were significantly upregulated in hepatoma cell lines.

Conclusion: The novel srlncRNA (AC026412.3, AL451069.3, and AL031985.3) signatures may provide insights into new therapies and prognosis predictions for patients with HCC.

Analysis of the role of m6A and lncRNAs in prognosis and immunotherapy of hepatocellular carcinoma

Objective

To study the role of m6A and lncRNAs in the prognosis and immunotherapy of hepatocellular carcinoma, construct the risk score of overall survival of hepatocellular carcinoma, and search for new therapeutic targets and drugs.

Methods

The data used in this study are obtained from The Cancer Genome Atlas (TCGA) database. A total of 424 HCC samples were included. The co-expression of lncRNAs and M6A-related genes in HCC was analyzed, COX regression analysis was conducted to construct the risk score for HCC prognosis, and the model's validity was further verified in different clinical trials subtypes and principal component analysis. GO enrichment analysis and immune function analysis were performed for the differential genes in the high-risk group and the low-risk group divided by risk score and analyzed the prognostic effect of TMB on the two groups. Based on the results, potential therapeutic agents for HCC were screened.

Results

The risk score can better predict the prognosis of HCC, the area under the ROC curve is 0.727. Differential genes were mainly located in the extracellular matrix and chromosomal regions and may play regulatory roles in binding sites and catalytic enzymes, thereby affecting the chromosome division and cell proliferation of cells. Type Ⅱ IFN response, type Ⅰ IFN response and MHC class Ⅰ were the three most different functions in terms of immune function between the high-risk group and the low-risk group. Type II IFN response, type I IFN response was significantly down-regulated in the high-risk group, while MHC class I was up-regulated. 14 potential therapeutic drugs were screened out.

Conclusions

The risk score constructed with NRAV and AL031985.3 had a good predictive effect on the prognosis of HCC. Differences in genes and immune function between high-risk and low-risk groups promoted the occurrence and progression of HCC.

Construction and validation of a necroptosis-related lncRNAs prognosis signature of hepatocellular carcinoma

Background: Immunotherapy has achieved remarkable success in treating advanced liver cancer. Current evidence shows that most of the available immune checkpoint inhibitor (ICB) treatments are suboptimal, and specific markers are needed for patients regarded as good candidates for immunotherapy. Necroptosis, a type of programmed cell death, plays an important role in hepatocellular carcinoma (HCC) progression and outcome. However, studies on the necroptosis-related lncRNA in HCC are scarce. In this view, the present study investigates the link among necroptosis-related lncRNA, prognosis, immune microenvironment, and immunotherapy response.

Methods: Gene transcriptome and clinical data were retrieved from The Cancer Genome Atlas database. Pearson correlation analysis of necroptosis-related genes was performed to identify necroptosis-related lncRNAs. The Wilcoxon method was used to detect differentially expressed genes, and prognostic relevant lncRNAs were obtained by univariate Cox regression analysis. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis were utilized to perform functional enrichment analysis. Lasso–Cox stepwise regression analysis was employed to calculate risk score, which was involved in analyzing immune cells infiltration, immune checkpoints expression, and predicting immunotherapeutic efficacy. Quantitative RT-PCR (qRT-PCR) was performed to detect the expression pattern of lncRNA in cell lines.

Results: The 10 lncRNAs generated in this study were used to create a prognostic risk model for HCC and group patients into groups based on risk. High-risk patients with HCC have a significantly lower OS rate than low-risk patients. Multivariate Cox regression analysis showed that risk score is an independent risk factor for HCC with high accuracy. Patients in the high-risk group exhibited a weaker immune surveillance and higher expression level of immune checkpoint molecules. In terms of drug resistance, patients in the low-risk group were more sensitive to sorafenib. The OS-related nomogram was constructed to verify the accuracy of our model. Finally, quantitative RT-PCR experiments were used to verify the expression patterns of candidate genes.

Conclusion: The lncRNA signature established herein, encompassing 10 necroptosis-related lncRNAs, is valuable for survival prediction and holds promise as prognostic markers for HCC.

Identification of HM13 as a prognostic indicator and a predictive biomarker for immunotherapy in hepatocellular carcinoma

Background

Histocompatibility minor 13 (HM13) is a signal sequence stubbed intramembrane cleavage catalytic protein that is essential for cell signaling, intracellular communication, and cancer. However, the expression of HM13 and its prognostic value, association with tumor-infiltrating immune cells (TIICs) in the microenvironment, and potential to predict immunotherapeutic response in HCC are unknown.

Methods

The HM13 expression, clinicopathology analysis, and its influence on survival were analyzed in multiple public databases and further verified in collected HCC and normal tissues by qRT-PCR and immunohistochemistry staining assay (IHC). Furthermore, the lentivirus vector encoding HM13-shRNA to manipulate HM13 expression was selected to investigate whether HM13 could influence the malignant growth and metastasis potential of HCC cells. Finally, significant impacts of HM13 on the HCC tumor microenvironment (TME) and reaction to immune checkpoint inhibitors were analyzed.

Results

Upregulated HM13 was substantially correlated with poor prognosis in patients with HCC, and could facilitate the proliferation and migratory potential of HCC cells. Additionally, patients with high HM13 expression might be more sensitive to immunotherapy.

Conclusions

HM13 might be a prognostic biomarker and potential molecular therapeutic target for HCC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09987-2.

The Pyroptosis-Related Long Noncoding RNA Signature Predicts Prognosis and Indicates Immunotherapeutic Efficiency in Hepatocellular Carcinoma

Pyroptosis was recently demonstrated to be an inflammatory form of gasdermin-regulated programmed cell death characterized by cellular lysis and the release of several proinflammatory factors and participates in tumorigenesis. However, the effects of pyroptosis-related long noncoding RNAs (lncRNAs) on hepatocellular carcinoma (HCC) have not yet been completely elucidated. Based on the regression coefficients of ZFPM2-AS1, KDM4A-AS1, LUCAT1, NRAV, CRYZL2P-SEC16B, AL031985.3, SNHG4, AL049840.5, AC008549.1, MKLN1-AS, AC099850.3, and LINC01224, HCC patients were classified into a low- or high-risk group. The high-risk score according to pyroptosis-related lncRNA signature was significantly associated with poor overall survival even after adjusting for age and clinical stage. Receiver operating characteristic curves and principal component analysis further supported the accuracy of the model. Our study revealed that a higher pyroptosis-related lncRNA risk score was significantly associated with tumor staging, pathological grade, and tumor-node-metastasis stages. The nomogram incorporating the pyroptosis-related lncRNA risk score and clinicopathological factors demonstrated good accuracy. Furthermore, we observed distinct tumor microenvironment cell infiltration characteristics between high- and low-risk tumors. Notably, based on the risk model, we found that the risk score is closely related to the expression of immune checkpoint genes, immune subtypes of tumors, and the sensitivity of HCC to chemotherapy drugs and immunotherapy. In conclusion, our novel risk score of pyroptosis-related lncRNA can serve as a promising prognostic biomarker for HCC patients and provide help for HCC patients to guide precision drug treatment and immunotherapy.

Identification of a novel prognostic signature for HCC and analysis of costimulatory molecule-related lncRNA AC099850.3

Costimulatory molecules are involved in initiation of anti-tumor immune responses while long non‐coding RNAs (lncRNAs) regulate the development of various cancers. However, the roles of lncRNA in hepatocellular carcinoma (HCC) have not been fully established. In this study, we aimed at identifying lncRNAs-related costimulatory molecules in HCC and to construct a prognostic signature for predicting the clinical outcomes for HCC patients. Data were downloaded from The Cancer Genome Atlas database for bioinformatics analyses. Costimulatory molecules were obtained from published literature. The R software, SPSS, and GraphPad Prism were used for statistical analyses. A risk model that is based on five costimulatory molecule-related lncRNAs was constructed using lasso and Cox regression analyses. Multivariate regression analysis revealed that the risk score could predict the prognostic outcomes for HCC. Samples in high- and low-risk groups exhibited significant differences in gene set enrichment and immune infiltration levels. Through colony formation and CCK8 assays, we found that AC099850.3 was strongly associated with HCC cell proliferation. We identified and validated a novel costimulatory molecule-related survival model. In addition, AC099850.3 was found to be closely associated with clinical stages and proliferation of HCC cells, making it a potential target for HCC treatment.

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

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

Non-Coding RNAs Implicated in the Tumor Microenvironment of Colorectal Cancer: Roles, Mechanisms and Clinical Study

Colorectal cancer (CRC) is one of the most common malignant tumors. The morbidity and mortality rates have been increasing all over the world. It is critical to elucidate the mechanism of CRC occurrence and development. However, tumor microenvironment (TME) includes immune cells, fibroblasts, endothelial cells, cytokines, chemokines and other components that affect the progression of CRC and patients’ prognosis. Non-coding RNAs (ncRNAs) including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), circular RNAs (circRNAs) without protein-coding ability have been shown to engage in tumor microenvironment-mediated angiogenesis and metastasis. Therefore, clarifying the mechanism of ncRNAs regulating the microenvironment is very important to develop the therapeutic target of CRC and improve the survival time of patients. This review focuses on the role and mechanism of ncRNAs in the CRC microenvironment and puts forward possible clinical treatment strategies.

Comprehensive Analysis of N6-Methyladenosine-Related Long Noncoding RNA Prognosis of Acute Myeloid Leukemia and Immune Cell Infiltration

N6-Methyladenosine-related long noncoding RNAs play an essential role in many cancers’ development. However, the relationship between m6A-related lncRNAs and acute myelogenous leukemia (AML) prognosis remains unclear. We systematically analyzed the association of m6A-related lncRNAs with the prognosis and tumor immune microenvironment (TME) features using the therapeutically applicable research to generate effective treatment (TARGET) database. We screened 315 lncRNAs associated with AML prognosis and identified nine key lncRNAs associated with m6A by the LASSO Cox analysis. A model was established based on these nine lncRNAs and the predictive power was explored in The Cancer Genome Atlas (TCGA) database. The areas under the ROC curve of TARGET and TCGA databases for ROC at 1, 3, and 5 years are 0.701, 0.704, and 0.696, and 0.587, 0.639, and 0.685, respectively. The nomogram and decision curve analysis (DCA) showed that the risk score was more accurate than other clinical indicators in evaluating patients’ prognoses. The clusters with a better prognosis enrich the AML pathways and immune-related pathways. We also found a close correlation between prognostic m6A-related lncRNAs and tumor immune cell infiltration. LAG3 expression at the immune checkpoint was lower in the worse prognostic cluster. In conclusion, m6A-related lncRNAs partly affected AML prognosis by remodeling the TME and affecting the anticarcinogenic ability of immune checkpoints, especially LAG3 inhibitors. The prognostic model constructed with nine key m6A-related lncRNAs can provide a method to assess the prognosis of AML patients in both adults and children.

Construction of an algorithm based on oncosis-related LncRNAs comprising the molecular subtypes and a risk assessment model in lung adenocarcinoma

Background

As an important non‐apoptotic cell death method, oncosis has been reported to be closely associated with tumors in recent years. However, few research reported the relationship between oncosis and lung cancer.

Methods

In this study, we established an oncosis‐based algorithm comprised of cluster grouping and a risk assessment model to predict the survival outcomes and related tumor immunity of patients with lung adenocarcinomas (LUAD). We selected 11 oncosis‐related lncRNAs associated with the prognosis (CARD8‐AS1, LINC00941, LINC01137, LINC01116, AC010980.2, LINC00324, AL365203.2, AL606489.1, AC004687.1, HLA‐DQB1‐AS1, and AL590226.1) to divide the LUAD patients into different clusters and different risk groups. Compared with patients in clsuter1, patients in cluster2 had a survival advantage and had a relatively more active tumor immunity. Subsequently, we constructed a risk assessment model to distinguish between patients into different risk groups, in which low‐risk patients tend to have a better prognosis. GO enrichment analysis revealed that the risk assessment model was closely related to immune activities. In addition, low‐risk patients tended to have a higher content of immune cells and stromal cells in tumor microenvironment, higher expression of PD‐1, CTLA‐4, HAVCR2, and were more sensitive to immune checkpoint inhibitors (ICIs), including PD‐1/CTLA‐4 inhibitors. The risk score had a significantly positive correlation with tumor mutation burden (TMB). The survival curve of the novel oncosis‐based algorithm suggested that low‐risk patients in cluster2 have the most obvious survival advantage.

Conclusion

The novel oncosis‐based algorithm investigated the prognosis and the related tumor immunity of patients with LUAD, which could provide theoretical support for customized individual treatment for LUAD patients.

Inflammatory Response-Related Long Non-Coding RNA Signature Predicts the Prognosis of Hepatocellular Carcinoma

Background

Hepatocellular carcinoma (HCC) is a high mortality malignant tumor with genetic and phenotypic heterogeneity, making predicting prognosis challenging. Meanwhile, the inflammatory response is an indispensable player in the tumorigenesis process and regulates the tumor microenvironment, which can affect the prognosis of tumor patients.

Methods

Using HCC samples in the TCGA-LIHC dataset, we explored lncRNA expression profiles associated with the inflammatory response. The inflammatory response-related lncRNA signature was constructed by univariate Cox regression, LASSO regression, and multivariate Cox regression methods based on inflammatory response-related differentially expressed lncRNAs in HCC.

Results

Seven inflammatory response-related lncRNA signatures were identified in predicting HCC prognosis. Kaplan–Meier (K-M) survival analysis indicated that high-risk group HCC patients were associated with poor prognosis. The utility of the inflammatory response-related lncRNA signatures was proved by the AUC and DCA analysis. The nomogram further confirmed the accuracy of the novel signature in predicting HCC patients' prognoses. In validation, our novel signature is more accurate than traditional clinicopathological performance for prognosis prediction of HCC patients. GSEA analysis further elucidated the underlying mechanisms and pathways of HCC progression in the low- and high-risk groups. Moreover, immune cells infiltration responses and immune function analyses revealed a significant difference between high- and low-risk groups in cytolytic activity, MHC class I, type I INF response, type II INF response, inflammation-promoting, and T cell coinhibition. Finally, HHLA2, NRP1, CD276, TNFRSF9, TNFSF4, CD80, and VTCN1 were expressed higher in high-risk groups in the immune checkpoint analysis.

Conclusions

A novel inflammatory response-related lncRNA signature (AC145207.5, POLHAS1, AL928654.1, MKLN1AS, AL031985.3, PRRT3AS1, and AC023157.2) is capable of predicting the prognosis of HCC patients and providing new immune targeted therapies insight.

Pyroptosis-Related LncRNA Signature Predicts Prognosis and Is Associated With Immune Infiltration in Hepatocellular Carcinoma

Pyroptosis is an inflammatory form of programmed cell death that is involved in various cancers, including hepatocellular carcinoma (HCC). Long non-coding RNAs (lncRNAs) were recently verified as crucial mediators in the regulation of pyroptosis. However, the role of pyroptosis-related lncRNAs in HCC and their associations with prognosis have not been reported. In this study, we constructed a prognostic signature based on pyroptosis-related differentially expressed lncRNAs in HCC. A co-expression network of pyroptosis-related mRNAs-lncRNAs was constructed based on HCC data from The Cancer Genome Atlas. Cox regression analyses were performed to construct a pyroptosis-related lncRNA signature (PRlncSig) in a training cohort, which was subsequently validated in a testing cohort and a combination of the two cohorts. Kaplan-Meier analyses revealed that patients in the high-risk group had poorer survival times. Receiver operating characteristic curve and principal component analyses further verified the accuracy of the PRlncSig model. Besides, the external cohort validation confirmed the robustness of PRlncSig. Furthermore, a nomogram based on the PRlncSig score and clinical characteristics was established and shown to have robust prediction ability. In addition, gene set enrichment analysis revealed that the RNA degradation, the cell cycle, the WNT signaling pathway, and numerous immune processes were significantly enriched in the high-risk group compared to the low-risk group. Moreover, the immune cell subpopulations, the expression of immune checkpoint genes, and response to chemotherapy and immunotherapy differed significantly between the high- and low-risk groups. Finally, the expression levels of the five lncRNAs in the signature were validated by quantitative real-time PCR. In summary, our PRlncSig model shows significant predictive value with respect to prognosis of HCC patients and could provide clinical guidance for individualized immunotherapy.

An Effective Hypoxia-Related Long Non-Coding RNA Assessment Model for Prognosis of Lung Adenocarcinoma

Background: Lung adenocarcinoma (LUAD) represents one of the highest incidence rates worldwide. Hypoxia is a significant biomarker associated with poor prognosis of LUAD. However, there are no definitive markers of hypoxia-related long non-coding RNAs (lncRNAs) in LUAD.

Methods: From The Cancer Genome Atlas (TCGA) and the Molecular Signatures Database (MSigDB), we acquired the expression of hypoxia-related lncRNAs and corresponding clinical information of LUAD patients. The hypoxia-related prognostic model was constructed by univariable COX regression analysis, least absolute shrinkage and selection operator (LASSO), and multivariable Cox regression analysis. To assess the performance of the model, the Kaplan–Meier (KM) survival and receiver operating characteristic (ROC) curve analyses were performed.

Results: We found seven lncRNAs, AC022613.1, AC026355.1, GSEC, LINC00941, NKILA, HSPC324, and MYO16-AS1, as biomarkers of the potential hypoxia-related prognostic signature. In the low-risk group, patients had a better overall survival (OS). In addition, the results of ROC analysis indicated that the risk score predicted LUAD prognosis exactly. Furthermore, combining the expression of lncRNAs with clinical features, two predictive nomograms were constructed, which could accurately predict OS and had high clinical application value.

Conclusion: In summary, the seven-lncRNA prognostic signature related to hypoxia might be useful in predicting clinical outcomes and provided new molecular targets for the research of LUAD patients.

Identification of a Novel Survival-Related circRNA–miRNA–mRNA Regulatory Network Related to Immune Infiltration in Liver Hepatocellular Carcinoma

Increasing studies have reported that circular RNAs (circRNAs) play critical roles in tumorigenesis and cancer progression. However, the underlying regulatory mechanisms of circRNA-related competing endogenous RNA (ceRNA) in liver hepatocellular carcinoma (LIHC) are still unclear. In the present study, we discovered dysregulated circRNAs through Gene Expression Omnibus (GEO) analysis and validated the expression of the top seven circRNAs with upregulated expression by qRT–PCR and Sanger sequencing. Then, the Cancer-Specific CircRNA Database (CSCD) was used to predict the downstream miRNAs of seven circRNAs, and expression and survival analyses through The Cancer Genome Atlas (TCGA) were performed to identify the key miRNA in LIHC. Thereafter, the hsa_circ_0017264-hsa-miR-195–5p subnetwork was successfully constructed. Subsequently, we predicted downstream target genes of hsa-miR-195–5p with TargetScan, miRDB, and mirtarbase and overlapped them with differentially expressed mRNAs to obtain 21 target genes. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed to predict the biological and functional roles of these target genes. Finally, with Pearson correlation and prognostic value analysis, a survival-related hsa_circ_0017264-hsa-miR-195-5p-CHEK1/CDC25A/FOXK1 axis was established. Gene set enrichment analysis (GSEA) was performed to determine the function of CHEK1/CDC25A/FOXK1 in the ceRNA network. Moreover, immune infiltration analysis revealed that the ceRNA network was markedly associated with the levels of multiple immune cell infiltrates, immune cell biomarkers and immune checkpoints. Overall, the hsa_circ_0017264-hsa-miR-195-5p-CHEK1/CDC25A/FOXK1 network might provide novel insights into the potential mechanisms underlying LIHC onset and progression.

Construction and Verification of a Combined Hypoxia and Immune Index for Clear Cell Renal Cell Carcinoma

Clear cell renal cell carcinoma (ccRCC) is one of the most aggressive malignancies in humans. Hypoxia-related genes are now recognized as a reflection of poor prognosis in cancer patients with cancer. Meanwhile, immune-related genes play an important role in the occurrence and progression of ccRCC. Nevertheless, reliable prognostic indicators based on hypoxia and immune status have not been well established in ccRCC. The aims of this study were to develop a new gene signature model using bioinformatics and open databases and to validate its prognostic value in ccRCC. The data used for the model structure can be accessed from The Cancer Genome Atlas database. Univariate, least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression analyses were used to identify the hypoxia- and immune-related genes associated with prognostic risk, which were used to develop a characteristic model of prognostic risk. Kaplan-Meier and receiver-operating characteristic curve analyses were performed as well as independent prognostic factor analyses and correlation analyses of clinical characteristics in both the training and validation cohorts. In addition, differences in tumor immune cell infiltrates were compared between the high and low risk groups. Overall, 30 hypoxia- and immune-related genes were identified, and five hypoxia- and immune-related genes (EPO, PLAUR, TEK, TGFA, TGFB1) were ultimately selected. Survival analysis showed that the high-risk score on the hypoxia- and immune-related gene signature was significantly associated with adverse survival outcomes. Furthermore, clinical ccRCC samples from our medical center were used to validate the differential expression of the five genes in tumor tissue compared to normal tissue through quantitative real-time polymerase chain reaction (qRT-PCR). However, more clinical trials are needed to confirm these results, and future experimental studies must verify the potential mechanism behind the predictive value of the hypoxia- and immune-related gene signature.

Immune-Related lncRNA Pairs Clinical Prognosis Model Construction for Hepatocellular Carcinoma

Background

Long non-coding RNA (lncRNA) plays an essential regulatory role in the occurrence and development of hepatocellular carcinoma (HCC). This paper aims to establish an immune-related lncRNA (irlncRNA) pairs model independent of expression level for risk assessment and prognosis prediction of HCC.

Methods

Transcriptome data and corresponding clinical data were downloaded from TCGA. HCC patients were randomly divided into training group and test group. Univariate Cox regression analysis, LASSO regression analysis, and stepwise multiple Cox regression analysis were used to establish a prognostic model. The prediction ability of the model was verified by ROC curves. Next, the patients were divided into low-risk and high-risk groups. We compared the differences between the two groups in survival rate, clinicopathological characteristics, tumor immune cell infiltration status, chemotherapeutic drug sensitivity and immunosuppressive molecules.

Results

A prognosis prediction model was established based on 7 irlncRNA pairs, namely irlncRNA pairs (IRLP). ROC curves of the training group and test group showed that the IRLP model had high sensitivity and specificity for survival prediction. Kaplan–Meier analysis showed that the survival rate of the high-risk group was significantly lower than that of the low-risk group. Immune cell infiltration analysis showed that the high-risk group was significantly correlated with various immune cell infiltration. Finally, there were statistically significant differences in chemosensitivity and molecular marker expression between the two groups.

Conclusion

The prognosis prediction model established by irlncRNA pairs has a certain guiding significance for the prognosis prediction of HCC. It may provide valuable clinical applications in antitumor immunotherapy.

Identification of Novel Tumor Microenvironment-Related Long Noncoding RNAs to Determine the Prognosis and Response to Immunotherapy of Hepatocellular Carcinoma Patients

Introduction: Hepatocellular carcinoma (HCC) is one of the most common malignant tumors with poor prognosis. The tumor microenvironment (TME) plays a vital role in HCC progression. Thus, this research was designed to analyze the correlation between the TME and the prognosis of HCC patients and to construct a TME-related long noncoding RNA (lncRNA) signature to determine HCC patients’ prognosis and response to immunotherapy.

Methods: We assessed the stromal–immune–estimate scores within the HCC microenvironment using the ESTIMATE (Estimation of Stromal and Immune Cells in Malignant Tumor Tissues Using Expression Data) algorithm based on The Cancer Genome Atlas database, and their associations with survival and clinicopathological parameters were also analyzed. Thereafter, differentially expressed lncRNAs were filtered out according to the immune and stromal scores. Cox regression analysis was performed to build a TME-related lncRNA risk signature. Kaplan–Meier analysis was used to explore the prognostic value of the risk signature. Furthermore, we explored the biological functions and immune microenvironment features in the high- and low-risk groups. Lastly, we probed the association of the risk model with treatment responses to immune checkpoint inhibitors (ICIs) in HCC.

Results: The stromal, immune, and estimate scores were obtained utilizing the ESTIMATE algorithm for patients with HCC. Kaplan–Meier analysis showed that high scores were significantly correlated with better prognosis in HCC patients. Six TME-related lncRNAs were screened to construct the prognostic model. The Kaplan–Meier curves suggested that HCC patients with low risk had better prognosis than those with high risk. Receiver operating characteristic (ROC) curve and Cox regression analyses indicated that the risk model could predict HCC survival exactly and independently. Functional enrichment analysis revealed that some tumor- and immune-related pathways were activated in the high-risk group. We also revealed that some immune cells, which were important in enhancing immune responses toward cancer, were significantly increased in the low-risk group. In addition, there was a close correlation between ICIs and the risk signature, which can be used to predict the treatment responses of HCC patients.

Conclusion: We analyzed the influence of the stromal, immune, and estimate scores on the prognosis of HCC patients. A novel TME-related lncRNA risk model was established, which could be effectively applied as an independent prognostic biomarker and predictor of ICIs for HCC patients.

Construction of a novel immune-related lncRNA signature and its potential to predict the immune status of patients with hepatocellular carcinoma

Background

The accuracy of existing biomarkers for predicting the prognosis of hepatocellular carcinoma (HCC) is not satisfactory. It is necessary to explore biomarkers that can accurately predict the prognosis of HCC.

Methods

In this study, original transcriptome data were downloaded from The Cancer Genome Atlas (TCGA) database. Immune-related long noncoding ribonucleic acids (irlncRNAs) were identified by coexpression analysis, and differentially expressed irlncRNA (DEirlncRNA) pairs were distinguished by univariate analysis. In addition, the least absolute shrinkage and selection operator (LASSO) penalized regression was modified. Next, the cutoff point was determined based on the area under the curve (AUC) and Akaike information criterion (AIC) values of the 5-year receiver operating characteristic (ROC) curve to establish an optimal model for identifying high-risk and low-risk groups of HCC patients. The model was then reassessed in terms of clinicopathological features, survival rate, tumor-infiltrating immune cells, immunosuppressive markers, and chemotherapy efficacy.

Results

A total of 1009 pairs of DEirlncRNAs were recognized in this study, 30 of these pairs were included in the Cox regression model for subsequent analysis. After regrouping according to the cutoff point, we could more effectively identify factors such as aggressive clinicopathological features, poor survival outcomes, specific immune cell infiltration status of tumors, high expression level of immunosuppressive biomarkers, and low sensitivity to chemotherapy drugs in HCC patients.

Conclusions

The nonspecific expression level signature involved with irlncRNAs shows promising clinical value in predicting the prognosis of HCC patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-09059-x.

Using Immune-Related Long Non-coding Ribonucleic Acids to Develop a Novel Prognosis Signature and Predict the Immune Landscape of Colon Cancer

Purpose: This study aimed to construct a novel signature to predict the survival of patients with colon cancer and the associated immune landscape, based on immune-related long noncoding ribonucleic acids (irlncRNAs).

Methods: Expression profiles of irlncRNAs in 457 patients with colon cancer were retrieved from the TCGA database (https://portal.gdc.cancer.gov). Differentially expressed (DE) irlncRNAs were identified and irlncRNA pairs were recognized using Lasso regression and Cox regression analyses. Akaike information criterion (AIC) values of receiver operating characteristic (ROC) curve were calculated to identify the ideal cut-off point for dividing patients into two groups and constructing the prognosis signature. Quantitative real-time polymerase chain reaction (qRT-PCR) was performed to validate the expression of LINC02195 and SCARNA9 in colon cancer.

Results: We identified 22 irlncRNA pairs and patients were divided into high-risk and low-risk groups based on the calculated risk score using these 22 irlncRNA pairs. The irlncRNA pairs were significantly related to patient survival. Low-risk patients had a significantly longer survival time than high-risk patients (p < 0.001). The area under the curve of the signature to predict 5-year survival was 0.951. The risk score correlated with tumor stage, infiltration depth, lymph node metastasis, and distant metastasis. The risk score remained significant after univariate and multivariate Cox regression analyses. A nomogram model to predict patient survival was developed based on the results of Cox regression analysis. Immune cell infiltration status, expression of some immune checkpoint genes, and sensitivity to chemotherapeutics were also related to the risk score. The results of qRT-PCR revealed that LINC02195 and SCARNA9 were significantly upregulated in colon cancer tissues.

Conclusion: The constructed prognosis signature showed remarkable efficiency in predicting patient survival, immune cell infiltration status, expression of immune checkpoint genes, and sensitivity to chemotherapeutics.

A novel lncRNA-miRNA-mRNA competing endogenous RNA regulatory network in lung adenocarcinoma and kidney renal papillary cell carcinoma

Background

GPRIN1 may be a novel tumor regulator, but its role and mechanism in tumors are still unclear.

Methods

First, a pan‐cancer correlation analysis was conducted on the expression and prognosis of GPRIN1 based on the data downloaded from The Cancer Genome Atlas (TCGA) database. Second, the Starbase database was used to predict the upstream miRNAs and lncRNAs of GPRIN1, and the expression analysis, survival analysis, and correlation analysis were performed to screen the microRNA (miRNAs)/long non‐coding RNAs (lncRNAs) that had a correlation with kidney renal papillary cell carcinoma (KIRP) or lung adenocarcinoma (LUAD). Third, the CIBERSORT algorithm was employed to calculate the proportion of various types of immune cells, and then the R packages were used for evaluating the relation between GPRIN1 expression and tumor immune cell infiltration as well as between GPRIN1 and the immune cell biomarker. Finally, the correlation analysis was made on GPRIN1 and immune checkpoints (CD274, CTLA4, and PDCD1).

Results

The pan‐cancer analysis suggested that GPRIN1 was up‐expressed in KIRP and LUAD, and it correlated with poor prognosis. LINC00894/MMP25‐AS1/SNHG1/LINC02298/MIR193BHG‐miR‐140‐3p was likely to be the most promising upstream regulation pathway of GPRIN1. Upexpression of LINC00894/MMP25‐AS1/SNHG1/LINC02298/MIR193BHG and downexpression of miR‐140‐3p were found relevant with poor outcomes of KIRP and LUAD. GPRIN1 expression was significantly correlated with tumor immune cell infiltration, immune cell biomarkers, and immune checkpoints.

Conclusions

The competitive endogenous (ceRNA) of miR‐140‐3p‐GPRIN1 axis and its upstream lncRNAs are closely related to KIRP and LUAD, and might affect the prognosis and therapeutic effect of KIRP and LUAD.