Immune-related long noncoding RNA signature for predicting survival and immune checkpoint blockade in hepatocellular carcinoma

Long noncoding RNAs (lncRNAs) show multiple functions, including immune response. Recently, the immune-related lncRNAs have been reported in some cancers. We first investigated the immune-related lncRNA signature as a potential target in hepatocellular carcinoma (HCC) survival. The training set (n = 368) and the independent external validation cohort (n = 115) were used. Immune genes and lncRNAs coexpression were constructed to identify immune-related lncRNAs. Cox regression analyses were perfumed to establish the immune-related lncRNA signature. Regulatory roles of this signature on cancer pathways and the immunologic features were investigated. The correlation between immune checkpoint inhibitors and this signature was examined. In this study, the immune-related lncRNA signature was identified in HCC, which could stratify patients into high- and low-risk groups. This immune-related lncRNA signature was correlated with disease progression and worse survival and was an independent prognostic biomarker. Our immune-related lncRNA signature was still a powerful tool in predicting survival in each stratum of age, gender, and tumor stage. This signature mediated cell cycle, glycolysis, DNA repair, mammalian target of rapamycin signaling, and immunologic characteristics (i.e., natural killer cells vs. Th1 cells down, etc). This signature was associated with immune cell infiltration (i.e., macrophages M0, Tregs, CD4 memory T cells, and macrophages M1, etc.,) and immune checkpoint blockade (ICB) immunotherapy-related molecules (i.e., PD-L1, PD-L2, and IDO1). Our findings suggested that the immune-related lncRNA signature had an important value for survival prediction and may have the potential to measure the response to ICB immunotherapy. This signature may guide the selection of the immunotherapy for HCC.     

Ferroptosis-related gene signature correlates with the tumor immune features and predicts the prognosis of glioma patients

Background: Glioma is a malignant intracranial tumor and the most fatal cancer. The role of ferroptosis in the clinical progression of gliomas is unclear.Method: Univariate and least absolute shrinkage and selection operator (Lasso) Cox regression methods were used to develop a ferroptosis-related signature (FRSig) using a cohort of glioma patients from the Chinese Glioma Genome Atlas (CGGA), and was validated using an independent cohort of glioma patients from The Cancer Genome Atlas (TCGA). A single-sample gene set enrichment analysis (ssGSEA) was used to calculate levels of the immune infiltration. Multivariate Cox regression was used to determine the independent prognostic role of clinicopathological factors and to establish a nomogram model for clinical application.Results: We analyzed the correlations between the clinicopathological features and ferroptosis-related gene (FRG) expression and established an FRSig to calculate the risk score for individual glioma patients. Patients were stratified into two subgroups with distinct clinical outcomes. Immune cell infiltration in the glioma microenvironment and immune-related indexes were identified that significantly correlated with the FRSig, the tumor mutation burden (TMB), copy number alteration (CNA), and immune checkpoint expression was also significantly positively correlated with the FRSig score. Ultimately, an FRSig-based nomogram model was constructed using the independent prognostic factors age, World Health Organization (WHO) grade, and FRSig score.Conclusion: We established the FRSig to assess the prognosis of glioma patients. The FRSig also represented the glioma microenvironment status. Our FRSig will contribute to improve patient management and individualized therapy by offering a molecular biomarker signature for precise treatment.     

Identification of m6A methyltransferase-related lncRNA signature for predicting immunotherapy and prognosis in patients with hepatocellular carcinoma

N6-methyladenosine (m6A) methyltransferase has been shown to be an oncogene in a variety of cancers. Nevertheless, the relationship between the long non-coding RNAs (lncRNAs) and hepatocellular carcinoma (HCC) remains elusive. We integrated the gene expression data of 371 HCC and 50 normal tissues from The Cancer Genome Atlas (TCGA) database. Differentially expressed protein-coding genes (DE-PCGs)/lncRNAs (DE-lncRs) analysis and univariate regression and Kaplan–Meier (K–M) analysis were performed to identify m6A methyltransferase-related lncRNAs. Three prognostic lncRNAs were selected by univariate and LASSO Cox regression analyses to construct the m6A methyltransferase-related lncRNA signature. Multivariate Cox regression analyses illustrated that this signature was an independent prognostic factor for overall survival (OS) prediction. The Gene Set Enrichment Analysis (GSEA) suggested that the m6A methyltransferase-related lncRNAs were involved in the immune-related biological processes (BPs) and pathways. Besides, we discovered that the lncRNAs signature was correlated with the tumor microenvironment (TME) and the expression of critical immune checkpoints. Tumor Immune Dysfunction and Exclusion (TIDE) analysis revealed that the lncRNAs could predict the clinical response to immunotherapy. Our study had originated a prognostic signature for HCC based on the potential prognostic m6A methyltransferase-related lncRNAs. The present study had deepened the understanding of the TME status of HCC patients and laid a theoretical foundation for the choice of immunotherapy.     

A robust six-miRNA prognostic signature for head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) remains a major health problem worldwide. We aimed to identify a robust microRNA (miRNA)-based signature for predicting HNSCC prognosis. The miRNA expression profiles of HNSCC were obtained from The Cancer Genome Atlas (TCGA) database. The TCGA HNSCC cohort was randomly divided into the discovery and validation cohort. A miRNA-based prognostic signature was built up based on TGCA discovery cohort, and then further validated. The downstream targets of prognostic miRNAs were subjected to functional enrichment analyses. The role of miR-1229-3p, a prognosis-related miRNA, in tumorigenesis of HNSCC was further evaluated. A total of 305 significantly differentially expressed miRNAs were found between HNSCC samples and normal tissues. A six-miRNA prognostic signature was constructed, which exhibited a strong association with overall survival (OS) in the TCGA discovery cohort. In addition, these findings were successfully confirmed in TCGA validation cohort and our own independent cohort. The miRNA-based signature was demonstrated as an independent prognostic indicator for HNSCC. A risk signature-based nomogram model was constructed and showed good performance for predicting the OS for HNSCC. The functional analyses revealed that the downstream targets of these prognostic miRNAs were closely linked to cancer progression. Mechanistically, in vitro analysis revealed that miR-1229-3p played a tumor promoting role in HNSCC. In conclusion, our study has developed a robust miRNA-based signature for predicting the prognosis of HNSCC with high accuracy, which will contribute to improve the therapeutic outcome.     

Survival prediction and response to immune checkpoint inhibitors: A prognostic immune signature for hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the most common cancers all over the world. Several studies have explored if immune-related genes and tumor immune microenvironment could play roles in HCC prognoses. This study is aimed at developing a prognostic signature of HCC based on immune-related genes or tumor immune microenvironment to predict survival and response to immune checkpoint inhibitors (ICIs). We constructed a prognostic signature using bioinformatics method and validated its predictive capability. The mechanisms of the signature prediction were explored with The Cancer Immunome Atlas (TCIA) and mutation analysis. We also explored the association between the signature and immunophenoscore (IPS), which is the marker of ICIs response. A 6 immune-related-gene (6-IRG) signature was developed. It was revealed in a multivariate analysis that the 6-IRG signature was an independent prognostic factor of overall survival and progression-free interval among HCC patients. In the high-risk group of 6-IRG signature score, macrophage M0 cells and regulatory T cells, which are observed associated with poor overall survival in our study, were higher. The low-risk group had a higher IPS, which meant a better response to ICIs. Taken together, we constructed a reliable 6-IRG signature for prediction of survival and response to ICIs. The signature needs further testing for clinical application.     

An Immune-Related Six-lncRNA Signature to Improve Prognosis Prediction of Glioblastoma Multiforme

Recent studies have demonstrated the utility and superiority of long non-coding RNAs (lncRNAs) as novel biomarkers for cancer diagnosis, prognosis, and therapy. In the present study, the prognostic value of lncRNAs in glioblastoma multiforme was systematically investigated by performing a genome-wide analysis of lncRNA expression profiles in 419 glioblastoma patients from The Cancer Genome Atlas (TCGA) project. Using survival analysis and Cox regression model, we identified a set of six lncRNAs (AC005013.5, UBE2R2-AS1, ENTPD1-AS1, RP11-89C21.2, AC073115.6, and XLOC_004803) demonstrating an ability to stratify patients into high- and low-risk groups with significantly different survival (median 0.899 vs. 1.611 years, p = 3.87e?09, log-rank test) in the training cohort. The six-lncRNA signature was successfully validated on independent test cohort of 219 patients with glioblastoma, and it revealed superior performance for risk stratification with respect to existing lncRNA-related signatures. Multivariate Cox and stratification analysis indicated that the six-lncRNA signature was an independent prognostic factor after adjusting for other clinical covariates. Further in silico functional analysis suggested that the six-lncRNA signature may be involved in the immune-related biological processes and pathways which are very well known in the context of glioblastoma tumorigenesis. The identified lncRNA signature had important clinical implication for improving outcome prediction and guiding the tailored therapy for glioblastoma patients with further prospective validation.     

Development and validation of a 14-gene signature for prognosis prediction in hepatocellular carcinoma

Worldwide, hepatocellular carcinoma (HCC) remains a crucial medical problem. Precise and concise prognostic models are urgently needed because of the intricate gene variations among liver cancer cells. We conducted this study to identify a prognostic gene signature with biological significance. We applied two algorithms to generate differentially expressed genes (DEGs) between HCC and normal specimens in The Cancer Genome Atlas cohort (training set included) and performed enrichment analyses to expound on their biological significance. A protein-protein interactions network was established based on the STRING online tool. We then used Cytoscape to screen hub genes in crucial modules. A multigene signature was constructed by Cox regression analysis of hub genes to stratify the prognoses of HCC patients in the training set. The prognostic value of the multigene signature was externally validated in two other sets from Gene Expression Omnibus (GSE14520 and GSE76427), and its role in recurrence prediction was also investigated. A total of 2000 DEGs were obtained, including 1542 upregulated genes and 458 downregulated genes. Subsequently, we constructed a 14-gene signature on the basis of 56 hub genes, which was a good predictor of overall survival. The prognostic signature could be replicated in GSE14520 and GSE76427. Moreover, the 14-gene signature could be applied for recurrence prediction in the training set and GSE14520. In summary, the 14-gene signature extracted from hub genes was involved in some of the HCC-related signalling pathways; it not only served as a predictive signature for HCC outcome but could also be used to predict HCC recurrence.     

Comprehensive analysis of prognostic immune-related genes in the tumor microenvironment of cutaneous melanoma

Cutaneous malignant melanoma (hereafter called melanoma) is one of the most aggressive cancers with increasing incidence and mortality rates worldwide. In this study, we performed a systematic investigation of the tumor microenvironmental and genetic factors associated with melanoma to identify prognostic biomarkers for melanoma. We calculated the immune and stromal scores of melanoma patients from the Cancer Genome Atlas (TCGA) using the ESTIMATE algorithm and found that they were closely associated with patients’ prognosis. Then the differentially expressed genes were obtained based on the immune and stromal scores, and prognostic immune-related genes further identified. Functional analysis and the protein–protein interaction network further revealed that these genes enriched in many immune-related biological processes. In addition, the abundance of six infiltrating immune cells was analyzed using prognostic immune-related genes by TIMER algorithm. The unsupervised clustering analysis using immune-cell proportions revealed eight clusters with distinct survival patterns, suggesting that dendritic cells were most abundant in the microenvironment and CD8⁺ T cells and neutrophils were significantly related to patients’ prognosis. Finally, we validated these genes in three independent cohorts from the Gene Expression Omnibus database. In conclusion, this study comprehensively analyzed the tumor microenvironment and identified prognostic immune-related biomarkers for melanoma.     

Nucleolar accumulation of RNA binding proteins induced by Actinomycin D is functional in Trypanosoma cruzi and Leishmania mexicana but not in T. brucei

We have recently shown in T. cruzi that a group of RNA Binding Proteins (RBPs), involved in mRNA metabolism, are accumulated into the nucleolus in response to Actinomycin D (ActD) treatment. In this work, we have extended our analysis to other members of the trypanosomatid lineage. In agreement with our previous study, the mechanism seems to be conserved in L. mexicana, since both endogenous RBPs and a transgenic RBP were relocalized to the nucleolus in parasites exposed to ActD. In contrast, in T. brucei, neither endogenous RBPs (TbRRM1 and TbPABP2) nor a transgenic RBP from T. cruzi were accumulated into the nucleolus under such treatment. Interestingly, when a transgenic TbRRM1 was expressed in T. cruzi and the parasites exposed to ActD, TbRRM1 relocated to the nucleolus, suggesting that it contains the necessary sequence elements to be targeted to the nucleolus. Together, both experiments demonstrate that the mechanism behind nucleolar localization of RBPs, which is present in T. cruzi and L. mexicana, is not functional in T. brucei, suggesting that it has been lost or retained differentially during the evolution of the trypanosomatid lineage.     

Integrated analysis of gene expression and DNA methylation datasets identified key genes and a 6-gene prognostic signature for primary lung adenocarcinoma

Lung adenocarcinoma (LUAD) is the main subtype of non-small cell lung cancer with a poor survival prognosis. In our study, gene expression, DNA methylation, and clinicopathological data of primary LUAD were utilized to identify potential prognostic markers for LUAD, which were recruited from The Cancer Genome Atlas (TCGA) database. Univariate regression analysis showed that there were 21 methylation-associated DEGs related to overall survival (OS), including 9 down- and 12 up-regulated genes. The 12 up-regulated genes with hypomethylation may be risky genes, whereas the other 9 down-regulated genes with hypermethylation might be protective genes. By using the Step-wise multivariate Cox analysis, a methylation-associated 6-gene (consisting of CCL20, F2, GNPNAT1, NT5E, B3GALT2, and VSIG2) prognostic signature was constructed and the risk score based on this gene signature classified patients into high- or low-risk groups. Patients of the high-risk group had shorter OS than those of the low-risk group in both the training and validation cohort. Multivariate Cox analysis and the stratified analysis revealed that the risk score was an independent prognostic factor for LUAD patients. The methylation-associated gene signature may serve as a prognostic factor for LUAD patients and the represent hypermethylated or hypomethylated genes might be potential targets for LUAD therapy.     

An atlas of posttranslational modifications on RNA binding proteins

RNA structure and function are intimately tied to RNA binding protein recognition and regulation. Posttranslational modifications are chemical modifications which can control protein biology. The role of PTMs in the regulation RBPs is not well understood, in part due to a lacking analysis of PTM deposition on RBPs. Herein, we present an analysis of posttranslational modifications (PTMs) on RNA binding proteins (RBPs; a PTM RBP Atlas). We curate published datasets and primary literature to understand the landscape of PTMs and use protein–protein interaction data to understand and potentially provide a framework for understanding which enzymes are controlling PTM deposition and removal on the RBP landscape. Intersection of our data with The Cancer Genome Atlas also provides researchers understanding of mutations that would alter PTM deposition. Additional characterization of the RNA–protein interface provided from in-cell UV crosslinking experiments provides a framework for hypotheses about which PTMs could be regulating RNA binding and thus RBP function. Finally, we provide an online database for our data that is easy to use for the community. It is our hope our efforts will provide researchers will an invaluable tool to test the function of PTMs controlling RBP function and thus RNA biology.     

Five key lncRNAs considered as prognostic targets for predicting pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) has a poor prognosis, and the 5‐year survival rate was only 7.7%. To improve prognosis, a screening biomarker for early diagnosis of pancreatic cancer is in urgent need. Long non‐coding RNA (lncRNA) expression profiles as potential cancer prognostic biomarkers play critical roles in development of tumorigenesis and metastasis of cancer. However, lncRNA signatures in predicting the survival of a patient with PDAC remain unknown. In the current study, we try to identify potential lncRNA biomarkers and their prognostic values in PDAC. LncRNAs expression profiles and corresponding clinical information for 182 cases with PDAC were acquired from The Cancer Genome Atlas (TCGA). A total of 14 470 lncRNA were identified in the cohort, and 175 PDAC patients had clinical variables. We obtained 108 differential expressed lncRNA via R packages. Univariate and multivariate Cox proportional hazards regression, lasso regression was performed to screen the potential prognostic lncRNA. Five lncRNAs have been recognized to significantly correlate with OS. We established a linear prognostic model of five lncRNA (C9orf139, MIR600HG, RP5‐965G21.4, RP11‐436K8.1, and CTC‐327F10.4) and divided patients into high‐ and low‐risk group according to the prognostic index. The five lncRNAs played independent prognostic biomarkers of OS of PDAC patients and the AUC of the ROC curve for the five lncRNAs signatures prediction 5‐year survival was 0.742. In addition, targeted genes of MIR600HG, C9orf139, and CTC‐327F10.4 were explored and functional enrichment was also conducted. These results suggested that this five‐lncRNAs signature could act as potential prognostic biomarkers in the prediction of PDAC patient's survival.     

Identification of immune-related hub genes and construction of an immune-related gene prognostic index for low-grade glioma

Low-grade glioma (LGG) poses significant management challenges and has a dismal prognosis. While immunotherapy has shown significant promise in cancer treatment, its progress in glioma has confronted with challenges. In our study, we aimed to develop an immune-related gene prognostic index (IRGPI) which could be used to evaluate the response and efficacy of LGG patients with immunotherapy. We included a total of 529 LGG samples from TCGA database and 1152 normal brain tissue samples from the GTEx database. Immune-related differentially expressed genes (DEGs) were screened. Then, we used weighted gene co-expression network analysis (WGCNA) to identify immune-related hub genes in LGG patients and performed Cox regression analysis to construct an IRGPI. The median IRGPI was used as the cut-off value to categorize LGG patients into IRGPI-high and low subgroups, and the molecular and immune mechanism in IRGPI-defined subgroups were analysed. Finally, we explored the relationship between IRGPI-defined subgroups and immunotherapy related indicators in patients after immunotherapy. Three genes (RHOA, NFKBIA and CCL3) were selected to construct the IRGPI. In a survival analysis using TCGA cohort as a training set, patients in the IRGPI-low subgroup had a better OS than those in IRGPI-high subgroup, consistent with the results in CGGA cohort. The comprehensive results showed that IRGPI-low subgroup had a more abundant activated immune cell population and lower TIDE score, higher MSI, higher TMB score, lower T cell dysfunction score, more likely benefit from ICIs therapy. IRGPI is a promising biomarker in the field of LGG ICIs therapy to distinguish the prognosis, the molecular and immunological characteristics of patients.     

Genome-wide screening and cohorts validation identifying novel lncRNAs as prognostic biomarkers for clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is the main subtype of renal cell carcinoma with varied prognosis. We aimed to identify and assess the possible prognostic long noncoding RNA (lncRNA) biomarkers. LncRNAs expression data and corresponding clinical information of 619 ccRCC patients were downloaded from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) databases. Differentially expressed genes analysis, univariate Cox regression, the least absolute shrinkage and selection operator Cox regression model were utilized to identify hub lncRNAs. Multivariate Cox regression was used to establish the risk model. Statistical analysis was performed using R 3.5.3. The expression value of five lncRNAs and the risk-score levels were significantly associated with a survival prognosis of ccRCC patients (all P < .001). In the TCGA validation cohort, the area under the curve (AUC) for the integrated nomogram was 0.905 and 0.91 for 3-, 5-year prediction separately. The AUC reached up to 0.757 in an independent ICGC cohort. Besides, the calibration plots also illustrated well curve-fitting between observation values and predictive values. Weighted gene co-expression network analysis and subsequent pathway analysis revealed that the PI3K-Akt-mTOR and hypoxia-inducible factor signaling crosstalk might function as the most essential mechanisms related to the five-lncRNAs signature. Our study suggested that lncRNA AC009654.1, AC092490.2, LINC00524, LINC01234, and LINC01885 were significantly associated with ccRCC prognosis. The prognostic model based on this five lncRNA may predict the overall survival of ccRCC.     

Development of a novel immune-related lncRNA signature as a prognostic classifier for endometrial carcinoma

Endometrial carcinoma (EnCa) is one of the deadliest gynecological malignancies. The purpose of the current study was to develop an immune-related lncRNA prognostic signature for EnCa. In the current research, a series of systematic bioinformatics analyses were conducted to develop a novel immune-related lncRNA prognostic signature to predict disease-free survival (DFS) and response to immunotherapy and chemotherapy in EnCa. Based on the newly developed signature, immune status and mutational loading between high‑ and low‑risk groups were also compared. A novel 13-lncRNA signature associated with DFS of EnCa patients was ultimately developed using systematic bioinformatics analyses. The prognostic signature allowed us to distinguish samples with different risks with relatively high accuracy. In addition, univariate and multivariate Cox regression analyses confirmed that the signature was an independent factor for predicting DFS in EnCa. Moreover, a predictive nomogram combined with the risk signature and clinical stage was constructed to accurately predict 1-, 2-, 3-, and 5-year DFS of EnCa patients. Additionally, EnCa patients with different levels of risk had markedly different immune statuses and mutational loadings. Our findings indicate that the immune-related 13-lncRNA signature is a promising classifier for prognosis and response to immunotherapy and chemotherapy for EnCa.     

MicroRNA signature predicts survival in papillary thyroid carcinoma

Papillary thyroid cancer (PTC) accounts for the majority of malignant thyroid tumors. Recently, several microRNA (miRNA) expression profiling studies have used bioinformatics to suggest miRNA signatures as potential prognostic biomarkers in various malignancies. However, a prognostic miRNA biomarker has not yet been established for PTC. The aim of the present study was to identify miRNAs with prognostic value for the overall survival (OS) of patients with PTC by analyzing high-throughput miRNA data and their associated clinical characteristics downloaded from The Cancer Genome Atlas database. From our dataset, 150 differentially expressed miRNAs were identified between tumor and nontumor samples; of these miRNAs, 118 were upregulated and 32 were downregulated. Among the 150 differentially expressed miRNAs, a four miRNA signature was identified that reliably predicts OS in patients with PTC. This miRNA signature was able to classify patients into a high-risk group and a low-risk group with a significant difference in OS (P < .01). The prognostic value of the signature was validated in a testing set ( P < .01). The four miRNA signature was an independent prognostic predictor according to the multivariate analysis and demonstrated good performance in predicting 5-year disease survival with an area under the receiver operating characteristic curve area under the curve (AUC) score of 0.886. Thus, this signature may serve as a novel biomarker for predicting the survival of patients with PTC.     

Developing a risk scoring system based on immune-related lncRNAs for patients with gastric cancer

The immune system and the tumor interact closely during tumor development. Aberrantly expressed long non-coding RNAs (lncRNAs) may be potentially applied as diagnostic and prognostic markers for gastric cancer (GC). At present, the diagnosis and treatment of GC patients remain a formidable clinical challenge. The present study aimed to build a risk scoring system to improve the prognosis of GC patients. In the present study, ssGSEA was used to evaluate the infiltration of immune cells in GC tumor tissue samples, and the samples were split into a high immune cell infiltration group and a low immune cell infiltration group. About 1262 differentially expressed lncRNAs between the high immune cell infiltration group and the low immune cell infiltration group. About 3204 differentially expressed lncRNAs between GC tumor tissues and paracancerous tissues were identified. Then, 621 immune-related lncRNAs were screened using a Venn analysis based on the above results, and 85 prognostic lncRNAs were identified using a univariate Cox analysis. We constructed a prognostic signature using LASSO analysis and evaluated the predictive performance of the signature using ROC analysis. GO and KEGG enrichment analyses were performed on the lncRNAs using the R package, ‘clusterProfiler’. The TIMER online database was used to analyze correlations between the risk score and the abundances of the six types of immune cells. In conclusion, our study found that specific immune-related lncRNAs were clinically significant. These lncRNAs were used to construct a reliable prognostic signature and analyzed immune infiltrates, which may assist clinicians in developing individualized treatment strategies for GC patients.