scholarly journals N6-Methyladenosine-Related lncRNA Signature Predicts the Overall Survival of Colorectal Cancer Patients

Genes ◽  
2021 ◽  
Vol 12 (9) ◽  
pp. 1375
Author(s):  
Wei Song ◽  
Jun Ren ◽  
Wenzheng Yuan ◽  
Rensheng Xiang ◽  
Yuhang Ge ◽  
...  

Background: The N6-methyladenosine (m6A) RNA modification can modify long non-coding RNAs (lncRNAs), thereby affecting the tumorigenesis and progression of tumors. However, the underlying role of m6A-modified lncRNAs in colorectal cancer (CRC) remains largely unknown. Therefore, our aim was to assess the prognostic value of m6A-modified lncRNAs in CRC patients. Methods: The gene expression and clinicopathological data of CRC were extracted from The Cancer Genome Atlas (TCGA) database. Pearson correlation analysis was used to investigate the m6A-modified lncRNAs. Consensus clustering was conducted to identify molecular subtypes of CRC, and the clinical significance of molecular subtypes was identified. The least absolute shrinkage and selection operator analysis (LASSO) was applied to establish a risk signature. Finally, a prognostic nomogram with risk score and clinicopathological variables was established. Results: In total, 29 m6A-modified lncRNAs were identified as prognostic lncRNAs. Two molecular clusters were identified and significant differences were found with respect to clinicopathological features and prognosis. Cluster1 is associated with poor overall survival (OS), down-regulation of Programmed cell death ligand-1 (PD-L1) expression, lower immune score, and less immune cell infiltration. Then, an m6A-modified lncRNA signature for predicting OS was constructed in the TCGA training cohort. The signature demonstrated favorable prediction performance in both training and validation sets. Compared with low-risk patients, patients with high risk showed worse clinical outcomes, lower immune scores, and downregulated PD-L1 expression. Further analysis indicated that the signature was an independent prognostic indicator, and then a prognostic nomogram based on risk score, tumor location, and tumor stage was established. Conclusions: Our study identified a seven m6A-modified lncRNA signature and established a prognostic nomogram that reliably predicts OS in CRC. These findings may improve the understanding of m6A modifications in CRC and provide insights into the prognosis and treatment strategy of CRC.

2020 ◽  
Author(s):  
Ran Wei ◽  
Jichuan Quan ◽  
Shuofeng Li ◽  
Zhao Lu ◽  
Xu Guan ◽  
...  

Abstract Background: Cancer stem cells (CSCs), which are characterized by self-renewal and plasticity, are highly correlated with tumor metastasis and drug resistance. To fully understand the role of CSCs in colorectal cancer (CRC), we evaluated the stemness traits and prognostic value of stemness-related genes in CRC.Methods: In this study, the data from 616 CRC patients from The Cancer Genome Atlas (TCGA) were assessed and subtyped based on the mRNA expression-based stemness index (mRNAsi). The correlations of cancer stemness with the immune microenvironment, tumor mutational burden (TMB) and N6-methyladenosine (m6A) RNA methylation regulators were analyzed. Weighted gene co-expression network analysis (WGCNA) was performed to identify the crucial stemness-related genes and modules. Furthermore, a prognostic expression signature was constructed using Lasso-penalized Cox regression analysis. The signature was validated via multiplex immunofluorescence staining of tissue samples in an independent cohort of 48 CRC patients.Results: This study suggests that high mRNAsi scores are associated with poor overall survival in stage Ⅳ CRC patients. Moreover, the levels of TMB and m6A RNA methylation regulators were positively correlated with mRNAsi scores, and low mRNAsi scores were characterized by increased immune activity in CRC. The analysis identified 2 key modules and 34 key genes as prognosis-related candidate biomarkers. Finally, a 3-gene prognostic signature (PARPBP, KNSTRN and KIF2C) was explored together with specific clinical features to construct a nomogram, which was successfully validated in an external cohort. Conclusions: There is a unique correlation between CSCs and the prognosis of CRC patients, and the novel biomarkers related to cell stemness could accurately predict the clinical outcomes of these patients.


2020 ◽  
Vol 38 (5_suppl) ◽  
pp. 25-25
Author(s):  
Yuanyuan Shen ◽  
Justin Hummel ◽  
Isabel Cristina Trindade ◽  
Christos Papageorgiou ◽  
Chi-Ren Shyu ◽  
...  

25 Background: Low cytotoxic T lymphocyte (CTLs) infiltration in colorectal cancer (CRC) tumors is a challenge to treatment with immune checkpoint inhibitors. Consensus molecular subtypes (CMS) classify patients based on tumor attributes, and CMS1 patients include the majority of patients with high CTL infiltration and “inflamed” tumors. Epigenetic modification plays a critical role in gene expression and therapy resistance. Therefore, in this study we compared DNA methylation, gene expression, and CTL infiltration of CMS1 patients to other CMS groups to determine targets for improving immunotherapy in CRC. Methods: RNA-seq (n = 511) and DNA methylation (n = 316) from The Cancer Genome Atlas databases were used to determine gene expression and methylation profiles based on CMSs. CMS1 was used as a reference and compared to other subtypes (CMS2-4). Microenvironment Cell Populations- counter (MCPcounter) was used to determine tumor CTL infiltration. Genes with significantly different expression (p < 0.01, LogFC≥|1.5|) and difference of mean methylation β value ≥|0.25| were integrated for Pearson correlation coefficient analysis with MCPcounter score (r > |0.7|). Results: Comparing CMS1 and CMS2, ARHGAP9, TBX21, and LAG3 were differentially methylated and correlated with CTL scores. ARHGAP9 and TBX21 were decreased and hypomethylated in CMS2. Comparing CMS1 and CMS3, ARHGAP9, TBX21, FMNL1, HLA-DPB1, and STX11 were downregulated in CMS3 and highly correlated with CTL scores. ARHGAP9, FMNL1, HLA-DPB1, and STX11 were hypomethylated in CMS3 and TBX21 was methylated in both, but had a higher methylation ratio in CMS1. Comparing CMS1 and CMS4, TBX21 was the only gene downregulated, hypomethylated, and highly correlated with CTL scores in CMS4 patients. Conclusions: We found six genes differentially expressed, differentially methylated, and highly correlated with CTL infiltration when comparing CMS1 to other CMS groups. Specifically, TBX21 was the only gene highly correlated with CTL scores with differential gene expression and methylation in CMS2-4 when compared to CMS1. Thus, T-bet may be a critical regulator of T cell responses in CRC.


2020 ◽  
Author(s):  
Chen Yang ◽  
Changhao Huang ◽  
Pengwei Zeng ◽  
Heyuan Huang ◽  
Zhikang Chen ◽  
...  

Abstract Background: B3GNT6 encodes the core 3 synthase in O-glycan biosynthesis. It is commonly expressed in the GI tract, while its clinical significance in colorectal cancer remains largely unknown.Methods: We gathered mRNA transcriptomic sequencing data from 3 Gene Expression Omnibus (GEO) datasets (GSE37182, GSE39582, GSE103512) and The Cancer Genome Atlas (TCGA) to compare the B3GNT6 mRNA level between colorectal cancer tissues and normal tissues and to evaluate its value as a prognostic marker. We further validated this in protein level using online database Human Protein Atlas and with immunohistochemical staining of B3GNT6 with our own cohort. Results: B3GNT6 expression was downregulated in colorectal cancer tissue compared with that in normal tissue in both mRNA and in protein level. Downregulation of B3GNT6 was associated with poor overall survival of colorectal cancer in GSE39582 and in TCGA database. Low B3GNT6 mRNA level was significantly associated with chromosome stable (CIN negative) and KRAS mutated group colorectal cancer patient. GSEA revealed that low B3GNT6 level in colorectal cancer is associated with upregulated proteasome activity.Conclusions: Downregulated B3GNT6 was correlated with poor overall survival of colorectal cancer patients. B3GNT6 could be used as a good prognostic marker in colorectal cancer.


Author(s):  
Xibo Zhao ◽  
Shanshan Cong ◽  
Qiuyan Guo ◽  
Yan Cheng ◽  
Tian Liang ◽  
...  

With the highest case-fatality rate among women, the molecular pathological alterations of ovarian cancer (OV) are complex, depending on the diversity of genomic alterations. Increasing evidence supports that immune infiltration in tumors is associated with prognosis. Therefore, we aim to assess infiltration in OV using multiple methods to capture genomic signatures regulating immune events to identify reliable predictions of different outcomes. A dataset of 309 ovarian serous cystadenocarcinoma patients with overall survival &gt;90 days from The Cancer Genome Atlas (TCGA) was analyzed. Multiple estimations and clustering methods identified and verified two immune clusters with component differences. Functional analyses pointed out immune-related alterations underlying internal genomic variables potentially. After extracting immune genes from a public database, the LASSO Cox regression model with 10-fold cross-validation was used for selecting genes associated with overall survival rate significantly, and a risk score model was then constructed. Kaplan–Meier survival and Cox regression analyses among cohorts were performed systematically to evaluate prognostic efficiency among the risk score model and other clinical pathological parameters, establishing a predictive ability independently. Furthermore, this risk score model was compared among identified signatures in previous studies and applied to two external cohorts, showing better prediction performance and generalization ability, and also validated as robust in association with immune cell infiltration in bulk tissues. Besides, a transcription factor regulation network suggested upper regulatory mechanisms in OV. Our immune risk score model may provide gyneco-oncologists with predictive values for the prognosis and treatment management of patients with OV.


2021 ◽  
Vol 10 ◽  
Author(s):  
Liang Zhao ◽  
Jiayue Zhang ◽  
Zhiyuan Liu ◽  
Yu Wang ◽  
Shurui Xuan ◽  
...  

Alternative splicing (AS) of pre-mRNA has been widely reported to be associated with the progression of malignant tumors. However, a systematic investigation into the prognostic value of AS events in glioblastoma (GBM) is urgently required. The gene expression profile and matched AS events data of GBM patients were obtained from The Cancer Genome Atlas Project (TCGA) and TCGA SpliceSeq database, respectively. 775 AS events were identified as prognostic factors using univariate Cox regression analysis. The least absolute shrinkage and selection operator (LASSO) cox model was performed to narrow down candidate AS events, and a risk score model based on several AS events were developed subsequently. The risk score-based signature was proved as an efficient predictor of overall survival and was closely related to the tumor purity and immunosuppression in GBM. Combined similarity network fusion and consensus clustering (SNF-CC) analysis revealed two distinct GBM subtypes based on the prognostic AS events, and the associations between this novel molecular classification and clinicopathological factors, immune cell infiltration, as well as immunogenic features were further explored. We also constructed a regulatory network to depict the potential mechanisms that how prognostic splicing factors (SFs) regulate splicing patterns in GBM. Finally, a nomogram incorporating AS events signature and other clinical-relevant covariates was built for clinical application. This comprehensive analysis highlights the potential implications for predicting prognosis and clinical management in GBM.


2020 ◽  
Author(s):  
Luping Zhang ◽  
Shaokun Wang ◽  
Yachen Wang ◽  
Weidan Zhao ◽  
Yingli Zhang ◽  
...  

Abstract Background: Imbalanced nutritional supply and demand in the tumor microenvironment often leads to hypoxia. The subtle interaction between hypoxia and immune cell behavior plays an important role in tumor occurrence and development. However, the functional relationship between hypoxia and the tumor microenvironment remains unclear. Therefore, we aimed to investigate the effect of hypoxia on the intestinal tumor microenvironment.Method: We extracted the names of hypoxia-related genes from the Gene Set Enrichment Analysis (GSEA) database and screened them for those associated with the prognosis of colorectal cancer, with the final list including ALDOB, GPC1, ALDOC, and SLC2A3. Using the sum of the expression levels of these four genes, provided by The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, and the expression coefficients, we developed a hypoxia risk score model. Using the median risk score value, we divided the patients in the two databases into high- and low-risk groups.GSEA was used to compare the enrichment differences between the two groups.We used the CIBERSORT computational method to analyze immune cell infiltration.Finally,the correlation between these five genes and hypoxia was analyzed. Result: The prognosis of the two groups differed significantly, with a higher survival rate in the low-risk group than in the high-risk group.We found that the different risk groups were enriched by immune-related and inflammatory pathways. We identified activated CD4 memory T cells and M0 macrophages in TCGA and GEO databases and found that CCL2/4/5, CSF1, and CX3CL1 contributed toward the increased infiltration rate of these immune cell types. Finally, we observed a positive correlation between the five candidate genes’ expression and the risk of hypoxia, with significant differences in the level of expression of each of these genes between patient risk groups.Conclusion: Overall, our data suggest that hypoxia is associated with the prognosis and rate of immune system infiltration in patients with colorectal cancer. This finding may improve immunotherapy for colorectal cancer.


2021 ◽  
Vol 8 ◽  
Author(s):  
Chen Jin ◽  
Rui Li ◽  
Tuo Deng ◽  
Jialiang Li ◽  
Yan Yang ◽  
...  

Hepatocellular carcinoma (HCC) is a highly invasive malignancy prone to recurrence, and patients with HCC have a low 5-year survival rate. Long non-coding RNAs (lncRNAs) play a vital role in the occurrence and development of HCC. N6-methyladenosine methylation (m6A) is the most common modification influencing cancer development. Here, we used the transcriptome of m6A regulators and lncRNAs, along with the complete corresponding clinical HCC patient information obtained from The Cancer Genome Atlas (TCGA), to explore the role of m6A regulator-related lncRNA (m6ARlnc) as a prognostic biomarker in patients with HCC. The prognostic m6ARlnc was selected using Pearson correlation and univariate Cox regression analyses. Moreover, three clusters were obtained via consensus clustering analysis and further investigated for differences in immune infiltration, immune microenvironment, and prognosis. Subsequently, nine m6ARlncs were identified with Lasso-Cox regression analysis to construct the prognostic signature m6A-9LPS for patients with HCC in the training cohort (n = 226). Based on m6A-9LPS, the risk score for each case was calculated. Patients were then divided into high- and low-risk subgroups based on the cutoff value set by the X-tile software. m6A-9LPS showed a strong prognosis prediction ability in the validation cohort (n = 116), the whole cohort (n = 342), and even clinicopathological stratified survival analysis. Combining the risk score and clinical characteristics, we established a nomogram for predicting the overall survival (OS) of patients. To further understand the mechanism underlying the m6A-9LPS-based classification of prognosis differences, KEGG and GO enrichment analyses, competitive endogenous RNA (ceRNA) network, chemotherapeutic agent sensibility, and immune checkpoint expression level were assessed. Taken together, m6A-9LPS could be used as a precise prediction model for the prognosis of patients with HCC, which will help in individualized treatment of HCC.


BMC Cancer ◽  
2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Jingyi Chen ◽  
Yuxuan Song ◽  
Mei Li ◽  
Yu Zhang ◽  
Tingru Lin ◽  
...  

Abstract Background Competing endogenous RNA (ceRNA) represents a class of RNAs (e.g., long noncoding RNAs [lncRNAs]) with microRNA (miRNA) binding sites, which can competitively bind miRNA and inhibit its regulation of target genes. Increasing evidence has underscored the involvement of dysregulated ceRNA networks in the occurrence and progression of colorectal cancer (CRC). The purpose of this study was to construct a ceRNA network related to the prognosis of CRC and further explore the potential mechanisms that affect this prognosis. Methods RNA-Seq and miRNA-Seq data from The Cancer Genome Atlas (TCGA) were used to identify differentially expressed lncRNAs (DElncRNAs), microRNAs (DEmiRNAs), and mRNAs (DEmRNAs), and a prognosis-related ceRNA network was constructed based on DElncRNA survival analysis. Subsequently, pathway enrichment, Pearson correlation, and Gene Set Enrichment Analysis (GSEA) were performed to determine the function of the genes in the ceRNA network. Gene Expression Profiling Interactive Analysis (GEPIA) and immunohistochemistry (IHC) were also used to validate differential gene expression. Finally, the correlation between lncRNA and immune cell infiltration in the tumor microenvironment was evaluated based on the CIBERSORT algorithm. Results A prognostic ceRNA network was constructed with eleven key survival-related DElncRNAs (MIR4435-2HG, NKILA, AFAP1-AS1, ELFN1-AS1, AC005520.2, AC245884.8, AL354836.1, AL355987.4, AL591845.1, LINC02038, and AC104823.1), 54 DEmiRNAs, and 308 DEmRNAs. The MIR4435-2HG- and ELFN1-AS1-associated ceRNA subnetworks affected and regulated the expression of the COL5A2, LOX, OSBPL3, PLAU, VCAN, SRM, and E2F1 target genes and were found to be related to prognosis and tumor-infiltrating immune cell types. Conclusions MIR4435-2HG and ELFN1-AS1 are associated with prognosis and tumor-infiltrating immune cell types and could represent potential prognostic biomarkers or therapeutic targets in colorectal carcinoma.


2021 ◽  
Author(s):  
Oscar Brück ◽  
Moon Hee Lee ◽  
Riku Turkki ◽  
Ilona Uski ◽  
Patrick Penttilä ◽  
...  

AbstractWhile the abundance and phenotype of tumor-infiltrating lymphocytes are linked with clinical survival, their spatial coordination and its clinical significance remain unclear. Here, we investigated the immune profile of intratumoral and peritumoral tissue of clear cell renal cell carcinoma patients (n = 64). We trained a cell classifier to detect lymphocytes from hematoxylin and eosin stained tissue slides. Using unsupervised classification, patients were further classified into immune cold, hot and excluded topographies reflecting lymphocyte abundance and localization. The immune topography distribution was further validated with The Cancer Genome Atlas digital image dataset. We showed association between PBRM1 mutation and immune cold topography, STAG1 mutation and immune hot topography and BAP1 mutation and immune excluded topography. With quantitative multiplex immunohistochemistry we analyzed the expression of 23 lymphocyte markers in intratumoral and peritumoral tissue regions. To study spatial interactions, we developed an algorithm quantifying the proportion of adjacent immune cell pairs and their immunophenotypes. Immune excluded tumors were associated with superior overall survival (HR 0.19, p = 0.02) and less extensive metastasis. Intratumoral T cells were characterized with pronounced expression of immunological activation and exhaustion markers such as granzyme B, PD1, and LAG3. Immune cell interaction occurred most frequently in the intratumoral region and correlated with CD45RO expression. Moreover, high proportion of peritumoral CD45RO+ T cells predicted poor overall survival. In summary, intratumoral and peritumoral tissue regions represent distinct immunospatial profiles and are associated with clinicopathologic characteristics.


2021 ◽  
Vol 12 ◽  
Author(s):  
Ruiyu Li ◽  
Yangzhige He ◽  
Hui Zhang ◽  
Jing Wang ◽  
Xiaoding Liu ◽  
...  

BackgroundPancreatic ductal adenocarcinoma (PDAC) remains treatment refractory. Immunotherapy has achieved success in the treatment of multiple malignancies. However, the efficacy of immunotherapy in PDAC is limited by a lack of promising biomarkers. In this research, we aimed to identify robust immune molecular subtypes of PDAC to facilitate prognosis prediction and patient selection for immunotherapy.MethodsA training cohort of 149 PDAC samples from The Cancer Genome Atlas (TCGA) with mRNA expression data was analyzed. By means of non-negative matrix factorization (NMF), we virtually dissected the immune-related signals from bulk gene expression data. Detailed immunogenomic and survival analyses of the immune molecular subtypes were conducted to determine their biological and clinical relevance. Validation was performed in five independent datasets on a total of 615 samples.ResultsApproximately 31% of PDAC samples (46/149) had higher immune cell infiltration, more active immune cytolytic activity, higher activation of the interferon pathway, a higher tumor mutational burden (TMB), and fewer copy number alterations (CNAs) than the other samples (all P &lt; 0.001). This new molecular subtype was named Immune Class, which served as an independent favorable prognostic factor for overall survival (hazard ratio, 0.56; 95% confidence interval, 0.33-0.97). Immune Class in cooperation with previously reported tumor and stroma classifications had a cumulative effect on PDAC prognostic stratification. Moreover, programmed cell death-1 (PD-1) inhibitors showed potential efficacy for Immune Class (P = 0.04). The robustness of our immune molecular subtypes was further verified in the validation cohort.ConclusionsBy capturing immune-related signals in the PDAC tumor microenvironment, we reveal a novel molecular subtype, Immune Class. Immune Class serves as an independent favorable prognostic factor for overall survival in PDAC patients.


Sign in / Sign up

Export Citation Format

Share Document