The joint role of methylation and immune-related lncRNAs in ovarian cancer: Defining molecular subtypes and developing prognostic signature

IntroductionComplex outcome of ovarian cancer (OC) stems from the tumor immune microenvironment (TIME) influenced by genetic and epigenetic factors. This study aimed to comprehensively explored the subclasses of OC through lncRNAs related to both N6-methyladenosine (m6A)/N1-methyladenosine (m1A)/N7-methylguanosine (m7G)/5-methylcytosine (m5C) in terms of epigenetic variability and immune molecules and develop a new set of risk predictive systems.Material and methodsThe lncRNA data of OC were collected from TCGA. Spearman correlation analysis on lncRNA data of OC with immune-related gene expression and with m6A/m5C/m1A/m7G were respectively conducted. The m6A/m5C/m1A/m7G-related m6A/m5C/m1A/m7G related immune lncRNA subtypes were identified on the basis of the prognostic lncRNAs. Heterogeneity among subtypes was evaluated by tumor mutation analysis, tumor microenvironment (TME) component analysis, response to immune checkpoint blocked (ICB) and chemotherapeutic drugs. A risk predictive system was developed based on the results of Cox regression analysis and random survival forest analysis of the differences between each specific cluster and other clusters.ResultsThree m6A/m5C/m1A/m7G-related immune lncRNA subtypes of OC showing distinct differences in prognosis, mutation pattern, TIME components, immunotherapy and chemotherapy response were identified. A set of risk predictive system consisting of 10 lncRNA for OC was developed, according to which the risk score of samples in each OC dataset was calculated and risk type was defined.ConclusionsThis study classified three m6A/m5C/m1A/m7G-related immune lncRNA subtypes with distinct heterogeneous mutation patterns, TME components, ICB therapy and immune response, and provided a set of risk predictive system consisted of 10 lncRNA for OC.     

Identification of an m6A-related lncRNA prognostic signature in colorectal cancer

Data sets of colorectal cancer (CRC) were obtained from The Cancer Genome Atlas (TCGA), three N6-methyladenosine (m6A) subtypes were identified using 21 m6A-related long noncoding RNAs (lncRNAs) and differential m6A subtypes of different CRC tumors were determined in this study to evaluate the m6A expression and the prognosis of patients with CRC. Subsequently, eight key lncRNAs were identified based on co-expression with 21 m6A-related genes in CRC tumors using the single-factor Cox and least absolute shrinkage and selection operator. Finally, an m6A-related lncRNA risk score model of CRC tumor was established using multifactor Cox regression based on the eight important lncRNAs and found to have a better performance in evaluating the prognosis of patients in the TCGA-CRC data set. TCGA-CRC tumor samples were divided based on the risk scores: high and low. Then, the clinical characteristics, tumor mutation load, and tumor immune cell infiltration difference between the high- and low-risk-score groups were explored, and the predictive ability of the risk score was assessed for immunotherapeutic benefits. We found that the risk score model can determine the overall survival, be a relatively independent prognostic indicator, and better evaluate the immunotherapeutic benefits for patients with CRC. This study provides data support for accurate immunotherapy in CRC.     

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.     

Pyroptosis-related lncRNA pairs to estimate the molecular features and prognostic outcomes of pancreatic ductal adenocarcinoma

Pyroptosis is a form of programmed cell death associated with inflammatory alterations. However, the intrinsic mechanisms and underlying correlation of pyroptosis-related lncRNAs (PRLs) in pancreatic ductal adenocarcinoma (PDAC) remain unclear. The objective of the current research was to identify pyroptosis-related lncRNAs and a prognostic model to predict the prognosis of patients. We extracted pyroptosis-related lncRNAs to construct a risk model and validated them at Fudan University Shanghai Cancer Center. Crosstalk between lncRNA SNHG10 and GSDMD was found to regulate pyroptosis levels. A new algorithm was used to establish a 0 or 1 PRL pair matrix and prognostic model. Six pyroptosis-related lncRNA pairs were identified and utilized to construct a risk model. The low-risk groups exhibited better prognoses than the high-risk groups. The area under the curve (AUC) indicated extremely high accuracy, reaching 0.810 at 1 year, 0.850 at 2 years, and 0.850 at 3 years in the training set. Patients with different risk scores exhibited distinct metabolic, inflammatory, and immune microenvironments as well as tumor mutation landscapes. Additionally, 9 commonly used chemotherapeutic drugs exhibited different sensitivities between the high- and low-risk groups. To conclude, we propose that pyroptosis exhibits a close correlation with PDAC. Our risk model based on PRL pairs may be beneficial for the accurate estimation of prognostic outcomes, the immune microenvironment, and drug sensitivity, bringing therapeutic hope for patients with PDAC.     

An exosome-related lncRNA signature correlates with prognosis,immune microenvironment,and therapeutic responses in hepatocellular carcinoma

BackgroundExosomes act as essential modulators of cancer development and progression in hepatocellular carcinoma. However, little is known about the potential prognostic value and underlying molecular features of exosome-related long non-coding RNAs.MethodsGenes associated with exosome biogenesis, exosome secretion, and exosome biomarkers were collected. Exosome-related lncRNA modules were identified using PCA and WGCNA analysis. A prognostic model based on data from the TCGA, GEO, NODE, and ArrayExpress was developed and validated. A comprehensive analysis of the genomic landscape, functional annotation, immune profile, and therapeutic responses underlying the prognostic signature was performed on multi-omics data, and bioinformatics methods were also applied to predict potential drugs for patients with high risk scores. qRT-PCR was used to validate the differentially expressed lncRNAs in normal and cancer cell lines.ResultsTwenty-six hub lncRNAs were identified as highly correlated with exosomes and overall survival and were used for prognosis modeling. Three cohorts consistently showed higher scores in the high-risk group, with an AUC greater than 0.7 over time. These higher scores implied poorer overall survival, higher genomic instability, higher tumor purity, higher tumor stemness, pro-tumor pathway activation, lower anti-tumor immune cell and tertiary lymphoid structure infiltration, and poor responses to immune checkpoint blockade therapy and transarterial chemoembolization therapy.ConclusionThrough developing an exosome-related lncRNA predictor for HCC patients, we revealed the clinical relevance of exosome-related lncRNAs and their potential as prognostic biomarkers and therapeutic response predictors.     

Development of a tumor microenvironment-related prognostic signature in glioma to predict immune landscape and potential therapeutic drugs

The involvement of the tumor microenvironment (TME) in the biology of gliomas has expanded, while it is yet uncertain its potential of supporting diagnosis and therapy choices. According to immunological characteristics and overall survival, cohorts of glioma patients from public databases were separated into two TME-relevant clusters in this analysis. Based on differentially expressed genes between TME clusters and correlative regression analysis, a 21-gene molecular classifier of TME-related prognostic signature (TPS) was constructed. Afterward, the prognostic efficacy and effectiveness of TPS were assessed in the training and validation groups. The outcome demonstrated that TPS might be utilized alone or in conjunction with other clinical criteria to act as a superior prognostic predictor for glioma. Also, high-risk glioma patients classified by TPS were considered to associate with enhanced immune infiltration, greater tumor mutation, and worse general prognosis. Finally, possible treatment medicines specialized for different risk subgroups of TPS were evaluated in drug databases.     

An eight-long noncoding RNA expression signature for colorectal cancer patients’ prognosis

Long noncoding RNAs (lncRNAs) have recently emerged as important biomarkers of cancer progression. Here, we proposed to develop a lncRNA-based signature with a prognostic value for colorectal cancer (CRC) overall survival (OS). Through mining microarray datasets, we analyzed the lncRNA expression profiles of 122 patients with CRC from Gene Expression Omnibus. Associations between lncRNA and CRC OS were firstly evaluated through univariate Cox regression analysis. A random survival forest method was applied for further screening of the lncRNA signature, which resulted in eight lncRNAs, including PEG3-AS1, LOC100505715, MINCR, DBH-AS1, LINC00664, FAM224A, LOC642852, and LINC00662. Combination of the eight lncRNAs weighted by their multivariate Cox regression coefficients formed a prognostic signature, through which, we could divide the 122 patients with CRC into two subgroups with significantly different OS. Good robustness of the lncRNA signature's prognostic value was verified through an independent data set consisting of 55 patients with CRC. In addition, gene set enrichment analysis indicated the potential association between high prognostic value and oxygen metabolism-related processes. This result should indicate that lncRNAs could be a useful signature for CRC prognosis.     

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.     

Profiles of immune-related genes and immune cell infiltration in the tumor microenvironment of diffuse lower-grade gliomas

The tumor microenvironment is highly correlated with tumor occurrence, progress, and prognosis. We aimed to investigate the immune-related gene (IRG) expression and immune infiltration pattern in the tumor microenvironment of lower-grade glioma (LGG). We employed the Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) algorithm to calculate immune and stromal scores and identify prognostic IRG based on The Cancer Genome Atlas data set. The potential molecular functions of these genes were explored with the help of functional enrichment analysis and the protein–protein interaction network. Remarkably, three cohorts that were downloaded from the Chinese Glioma Genome Atlas database were analyzed to further verify the prognostic values of these genes. Moreover, the Tumor IMmune Estimation Resource (TIMER) algorithm was used to estimate the abundance of infiltrating immune cells and explore the immune infiltration pattern in LGG. And unsupervised cluster analysis determined three clusters of the immune infiltration pattern and indicated that CD8⁺ T cells and macrophages were significantly associated with LGG outcomes. Altogether, our study identified a list of prognostic IRGs and provided a perspective to explore the immune infiltration pattern in LGG.     

Typical tumor immune microenvironment status determine prognosis in lung adenocarcinoma

BackgroundImmune cells, vital components of tumor microenvironment, regulate tumor survival and progression. Lung adenocarcinoma (LUAD), the tumor with the highest mortality rate worldwide, reconstitutes tumor immune microenvironment (TIME) to avoid immune destruction. Data have shown that TIME influences LUAD prognosis and predicts immunotherapeutic efficacy. The related information about the role of TIME's characteristics in LUAD is limited.MethodsWe performed unsupervised consensus clustering via machine-learning techniques to identify TIME clusters among 1906 patients and gathered survival data. The characteristics of TIME clusters of LUAD were visualized by multi-omics analysis, pseudo-time dynamic analysis, and enrichment analysis. TIME score model was constructed by principal component analysis. Comprehensive analysis and validation were conducted to test the prognostic efficacy and immunotherapeutic response of TIME score.ResultsTIME clusters (A, B and C) were constructed and exhibited different immune infiltration states. Multi-omics analyses included significant mutated genes (SMG), copy number variation (CNV) and cancer stemness that were significantly different among the three clusters. TIME cluster A had a lower SMG, lower CNV, and lower stemness but a higher immune infiltration level compared to TIME clusters B and C. TIME score showed that patients in low TIME score group had higher overall survival rates, higher immune infiltration level and high expression of immune checkpoints. In validation cohorts, low TIME score subgroup had better drug sensitivity and favorable immunotherapeutic response.ConclusionWe constructed a stable model of LUAD immune microenvironment characteristics that may improve the prognostic accuracy of patients, provide improved explanations of LUAD responses to immunotherapy, and provide new strategies for LUAD treatment.     

Identification of autophagy-related long non-coding RNA prognostic signature for breast cancer

Autophagy-related long non-coding RNAs (lncRNAs) disorders are related to the occurrence and development of breast cancer. The purpose of this study is to explore whether autophagy-related lncRNA can predict the prognosis of breast cancer patients. The autophagy-related lncRNAs prognostic signature was constructed by Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression. We identified five autophagy-related lncRNAs (MAPT-AS1, LINC01871, AL122010.1, AC090912.1, AC061992.1) associated with prognostic value, and they were used to construct an autophagy-related lncRNA prognostic signature (ALPS) model. ALPS model offered an independent prognostic value (HR = 1.664, 1.381-2.006), where this risk score of the model was significantly related to the TNM stage, ER, PR and HER2 status in breast cancer patients. Nomogram could be utilized to predict survival for patients with breast cancer. Principal component analysis and Sankey Diagram results indicated that the distribution of five lncRNAs from the ALPS model tends to be low-risk. Gene set enrichment analysis showed that the high-risk group was enriched in autophagy and cancer-related pathways, and the low-risk group was enriched in regulatory immune-related pathways. These results indicated that the ALPS model composed of five autophagy-related lncRNAs could predict the prognosis of breast cancer patients.     

Comprehensive analysis of novel three-long noncoding RNA signatures as a diagnostic and prognostic biomarkers of human triple-negative breast cancer

Currently, traditional predictors of prognosis (tumor size, nodal status, progesterone receptor [PR], estrogen receptor [ER], or human epidermal growth factor receptor-2 [HER2]) are insufficient for precise survival prediction for triple-negative breast cancer (TNBC). Long noncoding RNAs (lncRNAs) have been observed to exert critical functions in cancer, including in TNBC. Nevertheless, systematically tracking expression-based lncRNA biomarkers based on the sequence data for the prediction of prognosis in TNBC has not yet been investigated. To ascertain whether biomarkers exist that can distinguish TNBC from adjacent normal tissue or nTNBC, we implemented a comprehensive analysis of lncRNA expression profiles and clinical data of 1097 BC samples from The Cancer Genome Atlas database. A total of 1510 differentially expressed lncRNAs in normal and TNBC samples were extracted. Similarly, 672 differentially expressed lncRNAs between nTNBC and TNBC samples were detected. The receiver operating characteristic curve analysis indicated that three upregulated lncRNAs (AC091043.1, AP000924.1, and FOXCUT) may be of strong diagnostic value for predicting the existence of TNBC in the training and validation sets (area under the curve (AUC > 0.85). Kaplan-Meier analysis demonstrated that the other three lncRNAs (AC010343.3, AL354793.1, and FGF10-AS1) were associated with the prognosis of TNBC patients (P < 0.05). We used the three overall survival (OS)-related lncRNAs to establish a three-lncRNA signature. Multivariate Cox regression analysis suggested that the three-lncRNA signature was a prognostic factor independent of other clinical variables ( P < 0.01) for predicting OS in TNBC patients that could be utilized to classify patients into high- or low-risk subgroups. Our results might provide efficient signatures for clinical diagnosis and prognostic evaluation of TNBC.     

Hypoxia-related LncRNAs signature predicts prognosis and is associated with immune infiltration and progress of head and neck squamous cell carcinoma

BackgroundDisclosing prognostic information is necessary to enable good treatment selection and improve patient outcomes. Previous studies suggest that hypoxia is associated with an adverse prognosis in patients with HNSCC and that long non-coding RNAs (lncRNAs) show functions in hypoxia-associated cancer biology. Nevertheless, the understanding of lncRNAs in hypoxia related HNSCC progression remains confusing.MethodsData were downloaded from TCGA and GEO database. Bioinformatic tools including R packages GEOquery, limma, pheatmap, ggplot2, clusterProfiler, survivalROC and survcomp and LASSO cox analysis were utilized. Si-RNA transfection, CCK8 and real-time quantified PCR were used in functional study.ResultsGEO data (GSE182734) revealed that lncRNA regulation may be important in hypoxia related response of HNSCC cell lines. Further analysis in TCGA data identified 314 HRLs via coexpression analysis between differentially expressed lncRNAs and hypoxia-related mRNAs. 23 HRLs were selected to build the prognosis predicting model using lasso Cox regression analyses. Our model showed excellent performance in predicting survival outcomes among patients with HNSCC in both the training and validation sets. We also found that the risk scores were related to tumor stage and to tumor immune infiltration. Moreover, LINC01116 were selected as a functional study target. The knockdown of LINC01116 significantly inhibited the proliferation of HNSCC cells and effected the hypoxia induced immune and the NF-κB/AKT signaling.ConclusionsData analysis of large cohorts and functional experimental validation in our study suggest that hypoxia related lncRNAs play an important role in the progression of HNSCC, and its expression model can be used for prognostic prediction.     

Establishment and experimental validation of an immune miRNA signature for assessing prognosis and immune landscape of patients with colorectal cancer

As essential regulators of gene expression, miRNAs are engaged in the initiation and progression of colorectal cancer (CRC), including antitumour immune response. In this study, we proposed an integrated algorithm, ImmuMiRNA, for identifying miRNA modulators of immune-associated pathways. Based on these immune-associated miRNAs, we applied the LASSO algorithm to develop a reliable and individualized signature for evaluating overall survival (OS) and inflammatory landscape of CRC patients. An external public data set and qRT-PCR data from 40 samples were further utilized to validate this signature. As a result, an immune-associated miRNA prognostic signature (IAMIPS) consisting of three miRNAs (miR-194-3P, miR-216a-5p and miR-3677-3p) was established and validated. Patients in the high-risk group possessed worse OS. After stratification for clinical factors, the signature remained a powerful independent predictor for OS. IAMIPS displayed much better accuracy than the traditional clinical stage in assessing the prognosis of CRC. Further analysis revealed that patients in the high-risk group were characterized by inflammatory response, abundance immune cell infiltration, and higher immune checkpoint profiles and tumour mutation burden (TMB). In conclusion, the IAMIPS is highly predictive of OS in patients with CRC, which may serve as a powerful prognostic tool to further optimize immunotherapies for cancer.     

SPP1在头颈部鳞状细胞癌的免疫浸润及临床相关性分析

探讨分泌型磷蛋白1 (Secreted Phosphoprotein 1,SPP1)在头颈部鳞状细胞癌(Head and neck squamous cell carcinoma, HNSC)中与免疫浸润及临床的相关性,明确SPP1在HNSC预后和个体化治疗中的潜在价值。使用癌症基因组图谱(The Cancer Genome Atlas, TCGA)HNSC数据分析SPP1表达。使用来自TCGA的临床生存数据评估SPP1的临床预后价值。使用R语言的clusterProfiler包进行SPP1相关的富集分析。使用R语言的CIBERSORT函数评估22种肿瘤浸润免疫细胞在HNSC中的浸润情况,分析肿瘤浸润免疫细胞与SPP1表达之间的关联。差异表达分析发现SPP1在HNSC中高表达(P<0.001),临床相关性分析发现SPP1表达与T分期(P=0.001)、临床分期(P=0.013)相关,SPP1高表达患者的总生存期明显短于低表达患者(P=0.020 4)。基因富集分析发现SPP1在HNSC中与免疫学功能及免疫相关通路有关联。肿瘤浸润免疫细胞分析发现在高SPP1表达组中,M2巨噬细胞(...     

Hypoxia-induced PTTG3P contributes to colorectal cancer glycolysis and M2 phenotype of macrophage

Long noncoding RNAs (lncRNAs) play critical factors in tumor progression and are ectopically expressed in malignant tumors. Until now, lncRNA pituitary tumor-transforming 3, pseudogene (PTTG3P) biological function in colorectal cancer (CRC) further needs to be clarified. qRT-PCR was used to measure the PTTG3P level and CCK-8, glucose uptake, lactate assay, adenosine triphosphate (ATP) assay, extracellular acidification rate (ECAR) assay, and xenograft mice model were adopted to evaluate the glycolysis and proliferation, and macrophage polarization were determined in CRC cells. Xenograft experiments were utilized to analyze tumor growth. Ectopic expression of PTTG3P was involved in CRC and related to dismal prognosis. Through gain- and loss-of-function approaches, PTTG3P enhanced cell proliferation and glycolysis through YAP1. Further, LDHA knockdown or glycolysis inhibitor (2-deoxyglucose (2-DG), 3-BG) recovered from PTTG3P-induced proliferation. And PTTG3P overexpression could facilitate M2 polarization of macrophages. Silenced PTTG3P decreased the level of inflammatory cytokines TNF-α, IL-1β and IL-6, and low PTTG3P expression related with CD8⁺ T, NK, and TFH cell infiltration. Besides, hypoxia-inducible factor-1α (HIF1A) could increase PTTG3P expression by binding to the PTTG3P promoter region. Hypoxia-induced PTTG3P contributes to glycolysis and M2 phenotype of macrophage, which proposes a novel approach for clinical treatment.