The Role of m6A RNA Methylation-Related lncRNAs in the Prognosis and Tumor Immune Microenvironment of Papillary Thyroid Carcinoma.
Ontology highlight
ABSTRACT:
SUBMITTER: Wang W
PROVIDER: S-EPMC8762243 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundPapillary thyroid cancer has been associated with chronic inflammation. A systematic understanding of immune cell infiltration in PTC is essential for subsequent immune research and new diagnostic and therapeutic strategies.MethodsThree different algorithms, single-sample gene set enrichment analysis (ssGSEA), immune cell marker and CIBERSORT, were used to evaluate immune cell infiltration levels (abundance and proportion) in 10 data sets (The Cancer Genome Atlas [TCGA], GSE3467, GSE3678, GSE5364, GSE27155, GSE33630, GSE50901, GSE53157, GSE58545, and GSE60542; a total of 799 PTC and 194 normal thyroid samples). Consensus unsupervised clustering divided PTC patients into low-immunity and high-immunity groups. Weighted gene coexpression network analysis (WGCNA) and gene set enrichment analysis (GSEA) were used to analyze the potential mechanisms causing differences in the immune response.ResultsCompared with normal tissues, PTC tissues had a higher overall immune level and higher abundance levels and proportions of M2 macrophages, Tregs, monocytes, neutrophils, dendritic cells (DCs), mast cells (MCs), and M0 macrophages. Compared with early PTC, advanced PTC showed higher immune infiltration and higher abundance levels and proportions of M2 macrophages, Tregs, monocytes, neutrophils, DCs, MCs, and M0 macrophages. Compared to the low-immunity group, the high-immunity group exhibited more advanced stages, larger tumor sizes, greater lymph node metastases, higher tall-cell PTCs, lower follicular PTC proportions, more BRAF mutations, and fewer RAS mutations. Epstein-Barr virus (EBV) infection was the most significantly enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway for key module genes.ConclusionsIn human PTC, M2 macrophages, Tregs, monocytes, neutrophils, DCs, MCs, and M0 macrophages appear to play a tumor-promoting role, while M1 macrophages, CD8+ T cells, B cells, NK cells, and T follicular helper (TFH) cells (including eosinophils, γδ T cells, and Th17 cells with weak supporting evidence) appear to play an antitumor role. During the occurrence and development of PTC, the overall immune level was increased, and the abundance and proportion of tumor-promoting immune cells were significantly increased, indicating that immune escape had been aggravated. Finally, we speculate that EBV may play an important role in changing the immune microenvironment of PTC tumors.
| S-EPMC7642595 | biostudies-literature
Project description:Acute myeloid leukemia (AML) is a complex hematologic malignancy. Survival rate of AML patients is low. N6-methyladenosine (m6A) and long non-coding RNAs (lncRNAs) play important roles in AML tumorigenesis and progression. However, the relationship between lncRNAs and biological characteristics of AML, as well as how lncRNAs influence the prognosis of AML patients, remain unclear. In this study. In this study, Pearson correlation analysis was used to identify lncRNAs related to m6A regulatory genes, namely m6A-related lncRNAs. And we analyzed their roles and prognostic values in AML. m6A-related lncRNAs associated with patient prognosis were screened using univariate Cox regression analysis, followed by systematic analysis of the relationship between these genes and AML clinicopathologic and biologic characteristics. Furthermore, we examined the characteristics of tumor immune microenvironment (TIME) using different IncRNA clustering models. Using LASSO regression, we identified the risk signals related to prognosis of AML patients. We then constructed and verified a risk model based on m6A-related lncRNAs for independent prediction of overall survival in AML patients. Our results indicate that risk scores, calculated based on risk-related signaling, were related to the clinicopathologic characteristics of AML and level of immune infiltration. Finally, we examined the expression level of TRAF3IP2-AS1 in patient samples through real-time polymerase chain reaction analysis and in GEO datasets, and we identified a interaction relationship between SRSF10 and TRAF3IP2-AS1 through in vitro assays. Our study shows that m6A-related lncRNAs, evaluated using the risk prediction model, can potentially be used to predict prognosis and design immunotherapy in AML patients.
| S-EPMC8810799 | biostudies-literature
Project description:BackgroundN6-methyladenosine is the most abundant RNA modification, which plays a prominent role in various biology processes, including tumorigenesis and immune regulation. Multiple myeloma (MM) is the second most frequent hematological malignancy.Materials and methodsTwenty-two m6A RNA methylation regulators were analyzed between MM patients and normal samples. Kaplan-Meier survival analysis and least absolute shrinkage and selection operator (LASSO) Cox regression analysis were employed to construct the risk signature model. Receiver operation characteristic (ROC) curves were used to verify the prognostic and diagnostic efficiency. Immune infiltration level was evaluated by ESTIMATE algorithm and immune-related single-sample gene set enrichment analysis (ssGSEA).ResultsHigh expression of HNRNPC, HNRNPA2B1, and YTHDF2 and low expression of ZC3H13 were associated with poor survival. Based on these four genes, a prognostic risk signature model was established. Multivariate Cox regression analysis demonstrated that the risk score was an independent prognostic factor of MM. Enrichment analysis showed that cell cycle, immune response, MYC, proteasome, and unfold protein reaction were enriched in high-risk MM patients. Furthermore, patients with higher risk score exhibited lower immune scores and lower immune infiltration level.ConclusionThe m6A-based prognostic risk score accurately and robustly predicts the survival of MM patients and is associated with the immune infiltration level, which complements current prediction models and enhances our cognition of immune infiltration.
| S-EPMC8599583 | biostudies-literature
Project description:Recent studies have highlighted that N6-methyladenosine (m6A) plays a significant role in tumorigenicity and progression. However, the mechanism of m6A modifications in the tumor microenvironment (TME) immune cell infiltration in cervical cancer (CC) remains unclear. Clinical and RNA sequencing data of 25 m6A RNA methylation regulators were acquired from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. LASSO Cox regression analysis was used to generate a prognostic risk signature. m6A modification patterns were identified based on the expression of 25 m6A regulators, and their correlation with TME immune cell-infiltrating characterization was analyzed. Principal component analysis was used to construct an m6A-scoring signature (m6A score) to evaluate the m6A modification patterns of individual CC samples and guide the selection of more effective immunotherapeutic strategies. Genetic and expression alterations of 25 m6A regulators were highly heterogeneous between CC and normal tissues. METTL14 and IGF2BP1 were selected to conduct the prognostic risk signature. Three m6A modification patterns were identified in 659 CC samples, which were associated with distinct clinical outcomes and biological pathways. The TME immune cell-infiltrating characterization of the three m6A modification patterns was highly consistent with 3 tumor immune phenotypes, including immune-excluded, immune-inflamed, and immune-desert phenotypes. Due to the heterogeneity of m6A modification patterns, an m6A scoring signature was established to evaluate the m6A modification patterns of individual CC samples. Univariate and multivariate Cox regression analyses revealed that the m6A score is a robust and independent prognostic biomarker for assessing the prognosis of CC patients. A low m6A score, characterized by higher somatic mutation and higher expression of proliferation-related and DNA repair-related genes, indicated poor overall survival. Activation of immune infiltration was exhibited by the high m6A score, which was likely to have a good response and clinical benefits to antiPD-1/L1 immunotherapy. This study highlights the prognostic value of 25 m6A regulators in CC. The m6A modification is related to immune regulation and the formation of TME heterogeneity and complexity. An m6A scoring signature to clarify the individual m6A modification pattern could enhance our understanding of TME immune cell-infiltrating characterization and guide immunotherapy strategies.
| S-EPMC9239609 | biostudies-literature
Project description:BackgroundCervical lymph node metastasis is an important prognostic indicator for papillary thyroid carcinoma (PTC) and affects treatment strategies for PTC. lncRNAs essentially contribute to the biological functions of tumors. This study aimed to identify the lncRNAs associated with cervical lymph node metastasis and prognosis of PTC and their potential pathophysiological mechanisms.Materials and methodsPTC-associated lncRNAs were selected from The Cancer Genome Atlas database, and correlations among lncRNAs, lymph node metastasis, tumor staging, and prognosis of PTC were analyzed in silico. These correlations were then validated through quantitative reverse transcription PCR (qRT-PCR) and immunohistochemistry (IHC).ResultsIn silico analysis showed that FAM95B1 and UCA1 were significantly correlated with cervical lymph node metastasis, tumor staging, and PTC prognosis (P<0.05). qRT-PCR analysis revealed high UCA1 expression in PTC tissues and correlations between UCA1 expression levels and cervical lymph node metastasis and tumor staging in PTC, that is, higher UCA1 expression resulted in poorer PTC prognosis. IHC analysis revealed that a high expression of UCA1 was accompanied by a high expression of metastasis-related proteins (MMP-2 and MMP-9), thereby validating the correlation of UCA1 expression with metastasis.ConclusionFAM95B1 and UCA1 expression was significantly correlated with the occurrence and progression of PTC. The expression levels of UCA1 significantly affected the prognosis of PTC patients and were significantly correlated with tumor staging and cervical lymph node metastasis.
| S-EPMC6388952 | biostudies-literature
Project description:BackgroundN6-methyladenosine (m6A) modification and long non-coding RNAs (lncRNAs) play pivotal roles in gastric cancer (GC) progression. The emergence of immunotherapy in GC has created a paradigm shift in the approaches of treatment, whereas there is significant heterogeneity with regard to degree of treatment responses, which results from the variability of tumor immune microenvironment (TIME). How the interplay between m6A and lncRNAs enrolling in the shaping of TIME remains unclear.MethodsThe RNA sequencing and clinical data of GC patients were collected from TCGA database. Pearson correlation test and univariate Cox analysis were used to screen out m6A-related lncRNAs. Consensus clustering method was implemented to classify GC patients into two clusters. Survival analysis, the infiltration level of immune cells, Gene set enrichment analysis (GSEA) and the mutation profiles were analyzed and compared between two clusters. A competing endogenous RNA (ceRNA) network and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were applied for the identification of pathways in which m6A-related lncRNAs enriched. Then least absolute shrinkage and selection operator (LASSO) COX regression was implemented to select pivotal lncRNAs, and risk model was constructed accordingly. The prognosis value of the risk model was explored. In addition, the response to immune checkpoint inhibitors (ICIs) therapy were compared between different risk groups. Finally, we performed qRT-PCR to detect expression patterns of the selected lncRNAs in the 35 tumor tissues and their paired adjacent normal tissues, and validated the prognostic value of risk model in our cohort (N = 35).ResultsThe expression profiles of 15 lncRNAs were included to cluster patients into 2 subtypes. Cluster1 with worse prognosis harbored higher immune score, stromal score, ESTIMATE score and lower mutation rates of the genes. Different immune cell infiltration patterns were also displayed between the two clusters. GSEA showed that cluster1 preferentially enriched in tumor hallmarks and tumor-related biological pathways. KEGG pathway analysis found that the target mRNAs which m6A-related lncRNAs regulated by sponging miRNAs mainly enriched in vascular smooth muscle contraction, cAMP signaling pathway and cGMP-PKG signaling pathway. Next, eight lncRNAs were selected by LASSO regression algorithm to construct risk model. Patients in the high-risk group had poor prognoses, which were consistent in our cohort. As for predicting responses to ICIs therapy, patients from high-risk group were found to have lower tumor mutation burden (TMB) scores and account for large proportion in the Microsatellite Instability-Low (MSI-L) subtype. Moreover, patients had distinct immunophenoscores in different risk groups.ConclusionOur study revealed that the interplay between m6A modification and lncRNAs might have critical role in predicting GC prognosis, sculpting TIME landscape and predicting the responses to ICIs therapy.
| S-EPMC8943990 | biostudies-literature
Project description:N6-methyladenosine (m6A) is the most prevalent modification of RNA in mammals. m6A RNA methylation levels are dynamically regulated by m6A RNA methylation regulators. While increasing evidence has suggested that m6A RNA methylation is vital in the initiation and progression of human carcinoma, little is known about the expression and effect of m6A RNA methylation regulators in differentiated thyroid carcinoma (DTC). Herein, we demonstrate that most of the thirteen main m6A RNA methylation regulators are differentially expressed in DTC tissues and normal thyroid tissues. Based on consensus clustering of m6A RNA methylation regulators, DTC cases were divided into two subgroups (TC1 and TC2). Compared with the TC1 subgroup, the TC2 subgroup was associated with a poorer prognosis, older age, higher T grade, higher N grade and higher TNM stage. The results indicated that alteration of m6A RNA methylation regulators was closely related to DTC. We further established a risk signature of four m6A RNA methylation regulators that could evaluate prognosis and clinicopathological features in DTC. Finally, the results of the TCGA analysis were verified by other cohorts from Gene Expression Omnibus (GEO) database. In conclusion, m6A RNA methylation regulators play a crucial part in the progression of DTC and are potentially useful for evaluating the prognosis and providing potential novel insights into treatment strategies.
| S-EPMC7378910 | biostudies-literature
Project description:BackgroundPapillary thyroid carcinoma (PTC) is the most common malignant endocrine tumour, and its incidence and prevalence are increasing considerably. Cellular heterogeneity in the tumour microenvironment is important for PTC prognosis. Spatial transcriptomics is a powerful technique for cellular heterogeneity study.MethodsIn conjunction with a clinical pathologist identification method, spatial transcriptomics was employed to characterise the spatial location and RNA profiles of PTC-associated cells within the tissue sections. The spatial RNA-clinical signature genes for each cell type were extracted and applied to outlining the distribution regions of specific cells on the entire section. The cellular heterogeneity of each cell type was further revealed by ContourPlot analysis, monocle analysis, trajectory analysis, ligand-receptor analysis and Gene Ontology enrichment analysis.ResultsThe spatial distribution region of tumour cells, typical and atypical follicular cells (FCs and AFCs) and immune cells were accurately and comprehensively identified in all five PTC tissue sections. AFCs were identified as a transitional state between FCs and tumour cells, exhibiting a higher resemblance to the latter. Three tumour foci were shared among all patients out of the 13 observed. Notably, tumour foci No. 2 displayed elevated expression levels of genes associated with lower relapse-free survival in PTC patients. We discovered key ligand-receptor interactions, including LAMB3-ITGA2, FN1-ITGA3 and FN1-SDC4, involved in the transition of PTC cells from FCs to AFCs and eventually to tumour cells. High expression of these patterns correlated with reduced relapse-free survival. In the tumour immune microenvironment, reduced interaction between myeloid-derived TGFB1 and TGFBR1 in tumour focus No. 2 contributed to tumourigenesis and increased heterogeneity. The spatial RNA-clinical analysis method developed here revealed prognosis-associated cellular heterogeneity in the PTC microenvironment.ConclusionsThe occurrence of tumour foci No. 2 and three enhanced ligand-receptor interactions in the AFC area/tumour foci reduced the relapse-free survival of PTC patients, potentially leading to improved prognostic strategies and targeted therapies for PTC patients.
| S-EPMC10905537 | biostudies-literature
Project description:ObjectiveThe objective of this research was to screen prognostic related genes of thyroid papillary carcinoma (PTC) by single-cell RNA sequencing (scRNA-seq), to construct the diagnostic and prognostic models based on The Cancer Genome Atlas Thyroid Cancer (TCGA-THCA) data, and to evaluate the association between tumor immune microenvironment and the prognostic model.MethodThe differentially expressed genes (DEGs) and tumor evolution were analyzed by scRNA-seq based on public databases. The potential regulatory networks of DEGs related to prognosis were analyzed by multi-omics data in the THCA. Logistic regression and Cox proportional hazards regression were utilized to construct the diagnosis and prognostic model of PTC. The performance of the diagnostic model was verified by bulk RNA sequencing (RNA-seq) of our cohort. The tumor immune microenvironment associated with the prognostic model was evaluated using multi-omics data. In addition, qRT-PCR was performed on tumor tissues and adjacent normal tissues of 20 patients to verify the expression levels of DEGs.ResultsThe DEGs screened by scRNA-seq can distinguish between tumor and healthy samples. DEGs play different roles in the evolution from normal epithelial cells to malignant cells. Three DEGs ((FN1, CLU, and ANXA1)) related to prognosis were filtered, which may be regulated by DNA methylation, RNA methylation (m6A) and upstream transcription factors. The area under curve (AUC) of the diagnostic model based on 3-gene in the validation of our RNA-seq was 1. In the prognostic model based on 3-gene, the overall survival (OS) of high-risk patients was shorter. Combined with the clinical information of patients, a nomogram was constructed by using tumor size (pT) and risk score to quantify the prognostic risk. The age and tumor size of high-risk patients in the prognostic model were greater. In addition, the increase of tumor mutation burden (TMB) and diversity of T cell receptor (TCR), and the decrease of CD8+ T cells in high-risk group suggest the existence of immunosuppressive microenvironment.ConclusionWe applied the scRNA-seq pipeline to focus on epithelial cells in PTC, simulated the process of tumor evolution, and revealed a prognostic prediction model based on 3 genes, which is related to tumor immune microenvironment.
| S-EPMC8962641 | biostudies-literature
Project description:LncRNA N6-methylandenosine (m6A) modification has been shown to be associated with the constitution of the tumor microenvironment (TME) and tumorigenesis. It's essential to understand the mechanisms of lncRNA m6A modification in hepatocellular carcinoma (HCC) and identify relative prognostic predictors to guide therapy and explore potential therapeutic targets. Pearson correlation analysis was performed to identify m6A-related lncRNAs in 374 patients with HCC. Unsupervised cluster analysis of the potential m6A-related lncRNA-based HCC subtypes was conducted, followed by the concurrent analysis of their relationship with TME characteristics, immune checkpoints, immune features, and prognosis through single sample gene set enrichment analysis and ESTIMATE algorithm. Cox regression analyses were performed to screen prognostic m6A-related lncRNA, construct an m6A-related lncRNA signature (m6A-RLRS), and establish an integrated nomogram for the prognosis of patients with HCC. We identified 61 m6A-related lncRNAs and two HCC subtypes defined by consensus cluster of m6A-related lncRNAs with distinct clinical features. Progression-free survival (PFS), three TME-related scores, 15 immune-associated gene sets, and two immune checkpoints expression were found to be significantly different among the two subtypes. Twenty-five prognostic m6A-related lncRNAs were determined, four of which were included to establish an m6A-RLRS with favorable discrimination, and the signature was validated in the validation set and an independent FAHWMU cohort (n = 60). Furthermore, a novel nomogram combining signature and clinical predictors was generated with a C-index of 0.703, and an original ceRNA regulatory network consisting of 9 lncRNAs, 28 miRNAs, and 75 target mRNAs was constructed. Finally, the differential expression of four m6A-related lncRNA was verified by qRT-PCR. In conclusion, m6A-related lncRNA prognostic signature and molecular subtype contributes to accurately predict the prognosis of HCC and provide potential novel therapeutic targets.
| S-EPMC9361634 | biostudies-literature