Diagnostic, progressive and prognostic performance of m6A methylation RNA regulators in lung adenocarcinoma

Roles of naïve CD4+ T cells and their differentiated subtypes in lung adenocarcinoma and underlying potential regulatory pathways

BACKGROUND

Naïve CD4+ T cells and their differentiated counterparts play a significant regulatory role in the tumor immune microenvironment, yet their effects on lung adenocarcinoma (LUAD) are not fully understood.

METHODS

We utilized Mendelian randomization to assess the causal association between naïve CD4+ T cells and LUAD. Employing a modified single-sample Gene Set Enrichment Analysis (ssGSEA) algorithm with The Cancer Genome Atlas (TCGA) database, we determined the infiltration levels of naïve CD4+ T cells and their differentiation subtypes and investigated their correlation with clinical characteristics. Potential regulatory pathways of T helper cells were identified through Mantel tests and Kyoto Encyclopedia of Genes and Genomes (KEGG) database enrichment analysis.

RESULTS

Mendelian randomization analysis revealed an inhibitory effect of naïve CD4+ T cells on LUAD (false discovery rate < 0.05), which was corroborated by observational experiments using TCGA database. Specifically, T helper cell type 2 demonstrated a promotive effect on LUAD in terms of overall, disease-free, and progression-free survival (p < 0.05). Moreover, regulatory T cells exhibited a protective effect on LUAD in terms of disease-specific survival (p < 0.01). Concurrently, we explored the overall impact of naïve CD4+ T cell differentiation subtypes on LUAD, revealing upregulation in pathways such as neutrophil degranulation, MAPK family signaling pathways, and platelet activation, signaling, and aggregation.

CONCLUSION

Naïve CD4+ T cells and their differentiated counterparts play essential regulatory roles in the tumor immune microenvironment, demonstrating bidirectionality in their effects.Thus, elucidating the mechanisms and developing novel cell differentiation-inducing agents will benefit anti-cancer therapy.

Progression of m6A in the tumor microenvironment: hypoxia, immune and metabolic reprogramming

Recently, N6-methyladenosine (m6A) has aroused widespread discussion in the scientific community as a mode of RNA modification. m6A comprises writers, erasers, and readers, which regulates RNA production, nuclear export, and translation and is very important for human health. A large number of studies have found that the regulation of m6A is closely related to the occurrence and invasion of tumors, while the homeostasis and function of the tumor microenvironment (TME) determine the occurrence and development of tumors to some extent. TME is composed of a variety of immune cells (T cells, B cells, etc.) and nonimmune cells (tumor-associated mesenchymal stem cells (TA-MSCs), cancer-associated fibroblasts (CAFs), etc.). Current studies suggest that m6A is involved in regulating the function of various cells in the TME, thereby affecting tumor progression. In this manuscript, we present the composition of m6A and TME, the relationship between m6A methylation and characteristic changes in TME, the role of m6A methylation in TME, and potential therapeutic strategies to provide new perspectives for better treatment of tumors in clinical work.

N6-methyladenosine methylation regulatory pattern of pulmonary lymphoepithelioma-like carcinoma based on exosomal transcriptome analysis

Pulmonary lymphoepithelioma-like carcinoma (pLELC) is a rare malignancy that lacks specific biomarkers. N6-methyladenosine (m6 A) is the most widespread internal modification of messenger RNA (mRNA), and its dysregulation is involved in the development of many cancers. However, the expression of m6 A genes in pLELC and their roles are unknown. We obtained an exosomal transcriptome data set of patients diagnosed with pLELC and healthy controls using RNA sequencing and identified differentially expressed genes (DEGs) in the two groups using R software. The differential expression of the 37 m6 A genes in the two sets of samples was further analyzed, and receiver operating characteristic (ROC) curves were plotted for each gene to identify their grouping ability. The STRING database was used to construct a protein-protein interaction network for m6 A genes. An mRNA-miRNA regulatory network of m6 A-related DEGs was constructed using the miRNet database, and a prediction score formula was established. A nomogram was constructed based on the candidate m6 A genes and prediction scores. The expression of key genes was determined through the immunohistochemical (IHC) staining of clinical tissue sections. Using ROC curves, nine m6 A genes were revealed to have classification efficacy in both groups of samples. We screened seven m6 A-related DEGs (MAN2C1, HNRNPCL1, FUS, EIF6, DIP2A, COA3, and BUD13) that were beneficial for grouping and constructed nomogram models. Through IHC, we identified FUS and EIF6 as being possibly involved in the occurrence and development of pLELC. The m6 A gene expression patterns in pLELC-derived exosomes were significantly different from those in healthy controls. We screened several key genes to facilitate the development of diagnostic markers for pulmonary lymphoepithelioma.

Link between m6A modification and infiltration characterization of tumor microenvironment in lung adenocarcinoma

N6-methyladenosine (m6A) RNA methylation plays a pivotal role in immune responses and the onset and advancement of cancer. Nonetheless, the precise impact of m6A modification in lung adenocarcinoma (LUAD) and its associated tumor microenvironment (TME) remains to be fully elucidated. Here, we distinguished distinct m6A modification patterns within two separate LUAD cohorts using a set of 21 m6A regulators. The TME characteristics associated with these two patterns align with the immune-inflamed and immune-excluded phenotypes, respectively. We identified 2064 m6A-related genes, which were used as a basis to divide all LUAD samples into three distinct m6A gene clusters. We applied a scoring system to evaluate the m6A gene signature of the m6A modification pattern in individual patients. To authenticate the categorization significance of m6A modification patterns, we established a correlation between m6A score and TME infiltration profiling, tumor somatic mutations, and responses to immunotherapy. A high level of m6A modification may be associated with the aggressiveness and poor prognosis of LUAD. Further studies should investigate the mechanism of action of m6A regulators and m6A-related genes to improve the diagnosis and treatment of patients with LUAD.

Bayesian Approaches in Exploring Gene-environment and Gene-gene Interactions: A Comprehensive Review

Rapid advancements in high-throughput biological techniques have facilitated the generation of high-dimensional omics datasets, which have provided a solid foundation for precision medicine and prognosis prediction. Nonetheless, the problem of missing heritability persists. To solve this problem, it is essential to explain the genetic structure of disease incidence risk and prognosis by incorporating interactions. The development of the Bayesian theory has provided new approaches for developing models for interaction identification and estimation. Several Bayesian models have been developed to improve the accuracy of model and identify the main effect, gene-environment (G×E) and gene-gene (G×G) interactions. Studies based on single-nucleotide polymorphisms (SNPs) are significant for the exploration of rare and common variants. Models based on the effect heredity principle and group-based models are relatively flexible and do not require strict constraints when dealing with the hierarchical structure between the main effect and interactions (M-I). These models have a good interpretability of biological mechanisms. Machine learning-based Bayesian approaches are highly competitive in improving prediction accuracy. These models provide insights into the mechanisms underlying the occurrence and progression of complex diseases, identify more reliable biomarkers, and develop higher predictive accuracy. In this paper, we provide a comprehensive review of these Bayesian approaches.

The role of m6A demethylases in lung cancer: diagnostic and therapeutic implications

m6A is the most prevalent internal modification of eukaryotic mRNA, and plays a crucial role in tumorigenesis and various other biological processes. Lung cancer is a common primary malignant tumor of the lungs, which involves multiple factors in its occurrence and progression. Currently, only the demethylases FTO and ALKBH5 have been identified as associated with m6A modification. These demethylases play a crucial role in regulating the growth and invasion of lung cancer cells by removing methyl groups, thereby influencing stability and translation efficiency of mRNA. Furthermore, they participate in essential biological signaling pathways, making them potential targets for intervention in lung cancer treatment. Here we provides an overview of the involvement of m6A demethylase in lung cancer, as well as their potential application in the diagnosis, prognosis and treatment of the disease.

The critical roles of m6A RNA methylation in lung cancer: from mechanism to prognosis and therapy

Lung cancer, a highly malignant disease, greatly affects patients' quality of life. N6-methyladenosine (m6A) is one of the most common posttranscriptional modifications of various RNAs, including mRNAs and ncRNAs. Emerging studies have demonstrated that m6A participates in normal physiological processes and that its dysregulation is involved in many diseases, especially pulmonary tumorigenesis and progression. Among these, regulators including m6A writers, readers and erasers mediate m6A modification of lung cancer-related molecular RNAs to regulate their expression. Furthermore, the imbalance of this regulatory effect adversely affects signalling pathways related to lung cancer cell proliferation, invasion, metastasis and other biological behaviours. Based on the close association between m6A and lung cancer, various prognostic risk models have been established and novel drugs have been developed. Overall, this review comprehensively elaborates the mechanism of m6A regulation in the development of lung cancer, suggesting its potential for clinical application in the therapy and prognostic assessment of lung cancer.

A novel signature predicts prognosis and immunotherapy in lung adenocarcinoma based on cancer-associated fibroblasts

Background

Extensive research has established the significant correlations between cancer-associated fibroblasts (CAFs) and various stages of cancer development, including initiation, angiogenesis, progression, and resistance to therapy. In this study, we aimed to investigate the characteristics of CAFs in lung adenocarcinoma (LUAD) and develop a risk signature to predict the prognosis of patients with LUAD.

Methods

We obtained single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq data from the public database. The Seurat R package was used to process the scRNA-seq data and identify CAF clusters based on several biomarkers. CAF-related prognostic genes were further identified using univariate Cox regression analysis. To reduce the number of genes, Lasso regression was performed, and a risk signature was established. A novel nomogram that incorporated the risk signature and clinicopathological features was developed to predict the clinical applicability of the model. Additionally, we conducted immune landscape and immunotherapy responsiveness analyses. Finally, we performed in vitro experiments to verify the functions of EXO1 in LUAD.

Results

We identified 5 CAF clusters in LUAD using scRNA-seq data, of which 3 clusters were significantly associated with prognosis in LUAD. A total of 492 genes were found to be significantly linked to CAF clusters from 1731 DEGs and were used to construct a risk signature. Moreover, our immune landscape exploration revealed that the risk signature was significantly related to immune scores, and its ability to predict responsiveness to immunotherapy was confirmed. Furthermore, a novel nomogram incorporating the risk signature and clinicopathological features showed excellent clinical applicability. Finally, we verified the functions of EXP1 in LUAD through in vitro experiments.

Conclusions

The risk signature has proven to be an excellent predictor of LUAD prognosis, stratifying patients more appropriately and precisely predicting immunotherapy responsiveness. The comprehensive characterization of LUAD based on the CAF signature can predict the response of LUAD to immunotherapy, thus offering fresh perspectives into the management of LUAD patients. Our study ultimately confirms the role of EXP1 in facilitating the invasion and growth of tumor cells in LUAD. Nevertheless, further validation can be achieved by conducting in vivo experiments.

Identification of a m6A-related ferroptosis signature as a potential predictive biomarker for lung adenocarcinoma

BACKGROUND

Both N6-methyladenosine (m⁶A) and ferroptosis-related genes are associated with the prognosis of lung adenocarcinoma. However, the predictive value of m⁶A-related ferroptosis genes remains unclear. Here, we aimed to identify the prognostic value of m⁶A-related ferroptosis genes in lung adenocarcinoma.

METHODS

Lung adenocarcinoma sample data were downloaded from the University of California Santa Cruz Xena and Gene Expression Omnibus databases. Spearman's correlation analysis was used to screen for m⁶A-related ferroptosis genes. Univariate Cox regression, Kaplan-Meier, and Lasso analyses were conducted to identify prognostic m⁶A-related ferroptosis genes, and stepwise regression was used to construct a prognostic gene signature. The predictive value of the gene signature was assessed using a multivariate Cox analysis. In the validation cohort, survival analysis was performed to verify gene signature stability. The training cohort was divided into high- and low-risk groups according to the median risk score to assess differences between the two groups in terms of gene set variation analysis, somatic mutations, and tumor immune infiltration cells.

RESULTS

Six m⁶A-related ferroptosis genes were used to construct a gene signature in the training cohort and a multivariate Cox analysis was conducted to determine the independent prognostic value of these genes in lung adenocarcinoma. In the validation cohort, Kaplan-Meier and receiver operating characteristic analyses confirmed the strong predictive power of this signature for the prognosis of lung adenocarcinoma. Gene set variation analysis showed that the low-risk group was mainly related to immunity, and the high-risk group was mainly related to DNA replication. Somatic mutation analysis revealed that the TP53 gene had the highest mutation rate in the high-risk group. Tumor immune infiltration cell analysis showed that the low-risk group had higher levels of resting CD4 memory T cells and lower levels of M0 macrophages.

CONCLUSION

Our study identified a novel m⁶A-related ferroptosis-associated six-gene signature (comprising SLC2A1, HERPUD1, EIF2S1, ACSL3, NCOA4, and CISD1) for predicting lung adenocarcinoma prognosis, yielding a useful prognostic biomarker and potential therapeutic target.

Emerging role of N6-methyladenosine RNA methylation in lung diseases

In recent years, with the increase of air pollution, smoking, aging, and respiratory infection, the incidence rate and mortality of lung diseases are increasing annually, which has become a major hazard to human health. N6-methyladenosine (m⁶A) RNA methylation is the most abundant modifications in eukaryotes, and such modified RNA can be specifically recognized and combined by m⁶A recognition proteins and then mediate RNA splicing, maturation, enucleation, degradation, and translation. More and more studies have revealed that the m⁶A modification is involved in the pathogenesis and development of some diseases; however, the mechanisms of m⁶A in lung diseases are poorly understood. In this review, we summarize the latest progress in the biological function of m⁶A modifications in lung diseases and discuss the potential therapeutic and prognostic strategies. The dysregulation of global m⁶A levels and m⁶A regulators may affect the occurrence and development of asthma, chronic obstructive pulmonary disease, lung cancer, and other lung diseases through inflammation and immune function. In lung cancer, this modification has an important impact on malignant cell proliferation, migration, invasion, and drug resistance. In addition, abnormally changed m⁶A-modified proteins in lung cancer tissue samples and circulating tumor cells (CTCs) may be used as diagnostic and prognostic markers of lung cancer. Models composed of multiple m⁶A regulators can be used to evaluate the risk prediction or prognosis of asthma and pulmonary fibrosis. In general, the in-depth study of m⁶A modifications is a frontier direction in disease research. It provides novel insights for understanding of the molecular mechanisms underlying disease occurrence, development, and drug resistance, as well as for the development of effective novel therapeutics.

METTL16 predicts a favorable outcome and primes antitumor immunity in pancreatic ductal adenocarcinoma

Pancreatic carcinogenesis is a complicated and multi-step process. It is substantially assisted by N6-methyladenosine (m6A) RNA modification, especially when mutations of driver genes (KRAS, TP53, CDKN2A, and SMAD4) occur. However, the underlying mechanism remains obscure. In this research, we identified m6A regulators as potential biomarkers when mutations of driver genes occur, and investigated the role of these m6A candidates in pancreatic ductal adenocarcinoma (PDA). We first estimated the abnormal expression patterns of potential m6A regulators when all the driver genes are mutated, using The Cancer Genome Atlas and Gene Expression Omnibus databases. METTL16, an m6A“writer,” was chosen as a unique candidate of PDA, owing to its markedly differential expression under mutations of all driver genes (KRAS, TP53, CDKN2A, and SMAD4) and its favorable prognostic value. Moreover, METTL16 was under-expressed in PDA tissues and cell lines. Consistently, gain- and loss-of-function experiments indicated that it had a tumor suppressor role in vitro and in vivo. Further, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses revealed that METTL16 may have an effect on the tumor microenvironment. Notably, a markedly positive association between METTL16 expression and infiltration of B cells and CD8+ T cells was observed according to the CIBERSORT and TIMER databases. Enhanced expression of immune checkpoints and cytokines was elicited in patients with over-expression of METTL16. Notably, decreased expression of PD-L1 was observed when upregulation of METTL16 expression occurred in MIA PaCa-2 cells, while increased expression of PD-L1 existed when downregulation of METTL16 happened in HPAF-II cells. Collectively, these findings highlight the prognostic value of METTL16, and indicate that it is a potential immunotherapy target that could be used to regulate the tumor microenvironment and promote antitumor immunity in PDA.

Gene amplification-driven RNA methyltransferase KIAA1429 promotes tumorigenesis by regulating BTG2 via m6A-YTHDF2-dependent in lung adenocarcinoma

Background

Epigenetic alterations have been shown to contribute immensely to human carcinogenesis. Dynamic and reversible N6‐methyladenosine (m6A) RNA modification regulates gene expression and cell fate. However, the reasons for activation of KIAA1429 (also known as VIRMA, an RNA methyltransferase) and its underlying mechanism in lung adenocarcinoma (LUAD) remain largely unexplored. In this study, we aimed to clarify the oncogenic role of KIAA1429 in the tumorigenesis of LUAD.

Methods

Whole‐genome sequencing and transcriptome sequencing of LUAD data were used to analyze the gene amplification of RNA methyltransferase. The in vitro and in vivo functions of KIAA1429 were investigated. Transcriptome sequencing, methylated RNA immunoprecipitation sequencing (MeRIP‐seq), m6A dot blot assays and RNA immunoprecipitation (RIP) were performed to confirm the modified gene mediated by KIAA1429. RNA stability assays were used to detect the half‐life of the target gene.

Results

Copy number amplification drove higher expression of KIAA1429 in LUAD, which was correlated with poor overall survival. Manipulating the expression of KIAA1429 could regulate the proliferation and metastasis of LUAD. Mechanistically, the target genes of KIAA1429‐mediated m6A modification were confirmed by transcriptome sequencing and MeRIP‐seq assays. We also revealed that KIAA1429 could regulate BTG2 expression in an m6A‐dependent manner. Knockdown of KIAA1429 significantly decreased the m6A levels of BTG2 mRNA, leading to enhanced YTH m6A RNA binding protein 2 (YTHDF2, the m6A “reader”)‐dependent BTG2 mRNA stability and promoted the expression of BTG2; thus, participating in the tumorigenesis of LUAD.

Conclusions

Our data revealed the activation mechanism and important role of KIAA1429 in LUAD tumorigenesis, which may provide a novel view on the targeted molecular therapy of LUAD.

An analysis of the role of HnRNP C dysregulation in cancers

Heterogeneous nuclear ribonucleoproteins C (HnRNP C) is part of the hnRNP family of RNA-binding proteins. The relationship between hnRNP C and cancers has been extensively studied, and dysregulation of hnRNP C has been found in many cancers. According to existing public data, hnRNP C could promote the maturation of new heterogeneous nuclear RNAs (hnRNA s, also referred to as pre-mRNAs) into mRNAs and could stabilize mRNAs, controlling their translation. This paper reviews the regulation and dysregulation of hnRNP C in cancers. It interacts with some cancer genes and other biological molecules, such as microRNAs (miRNAs), long noncoding RNAs (lncRNAs), and double-stranded RNAs (dsRNAs). Even directly binds to them. The effects of hnRNP C on biological processes such as alternative cleavage and polyadenylation (APA) and N6-methyladenosine (m6A) modification differ among cancers. Its main function is regulating stability and level of translation of cancer genes, and the hnRNP C is regarded as a candidate biomarker and might be valuable for prognosis evaluation.

Insight into the structure, physiological function, and role in cancer of m6A readers—YTH domain-containing proteins

YT521-B homology (YTH) domain-containing proteins (YTHDF1-3, YTHDC1-2) are the most crucial part of N6-methyladenosine (m6A) readers and play a regulatory role in almost all stages of methylated RNA metabolism and the progression of various cancers. Since m6A is identified as an essential post-transcriptional type, YTH domain-containing proteins have played a key role in the m6A sites of RNA. Hence, it is of great significance to study the interaction between YTH family proteins and m6A-modified RNA metabolism and tumor. In this review, their basic structure and physical functions in RNA transcription, splicing, exporting, stability, and degradation as well as protein translation are introduced. Then we discussed the expression regulation of YTH domain-containing proteins in cancers. Furthermore, we introduced the role of the YTH family in cancer biology and systematically demonstrated their functions in various aspects of tumorigenesis and development. To provide a more institute understanding of the role of YTH family proteins in cancers, we summarized their functions and specific mechanisms in various cancer types and presented their involvement in cancer-related signaling pathways.

Comprehensive Analysis Revealed the Potential Implications of m6A Regulators in Lung Adenocarcinoma

Background: The biological significance of RNA N6-methyladenosine (m6A) decoration in tumorigenicity and progression has been highlighted in recent studies, but whether m6A modification plays a potential role in tumor microenvironment (TME) formation and immune regulation in lung adenocarcinoma (LUAD) remains unclear.

Methods: m6A modification features were evaluated by analyzing the multi-omics features of 17 m6A regulators in over 1900 LUAD samples, and at the same time, the correlation between these modification patterns and TME characteristics was analyzed. An m6A score signature–based principal component analysis (PCA) algorithm was constructed to assess the prognosis and responses of individual patients to immunotherapeutic and targeted therapies.

Results: Three different m6A modification patterns were determined in 1901 LUAD samples, which were found to be related to diverse clinical outcomes via different biological pathways. Based on the m6A score extracted from the m6A-associated signature genes, LUAD patients were separated into high- and low-m6A score groups. It was discovered that patients with high m6A scores had longer survival, lower tumor mutation loads, and low PD-L1/PDCD1/CTLA4/TAG3 expression level. In addition, LUAD patients with high m6A scores displayed lower IC₅₀ to some targeted drugs, including nilotinib, erlotinib, imatinib, and lapatinib.

Conclusion: m6A modification was significantly associated with the TME and clinical outcomes. These findings may help gain more insights into the role of m6A decoration in the molecular mechanism of LUAD, thus facilitating the development of more effective personalized treatment strategies.

The molecular mechanism of METTL3 promoting the malignant progression of lung cancer

Lung cancer remains one of the major causes of cancer-related death globally. Recent studies have shown that aberrant m⁶A levels caused by METTL3 are involved in the malignant progression of various tumors, including lung cancer. The m⁶A modification, the most abundant RNA chemical modification, regulates RNA stabilization, splicing, translation, decay, and nuclear export. The methyltransferase complex plays a key role in the occurrence and development of many tumors by installing m⁶A modification. In this complex, METTL3 is the first identified methyltransferase, which is also the major catalytic enzyme. Recent findings have revealed that METTL3 is remarkably associated with different aspects of lung cancer progression, influencing the prognosis of patients. In this review, we will focus on the underlying mechanism of METT3 in lung cancer and predict the future work and potential clinical application of targeting METTL3 for lung cancer therapy.

RNA-Binding Protein IGF2BP1 Associated With Prognosis and Immunotherapy Response in Lung Adenocarcinoma

N6-methyladenosine (m⁶A) is the most common modification in eukaryotic RNAs and plays a vital role in the tumorigenesis and metastasis of various cancers. However, a comprehensive study of m⁶A methylation regulators in lung adenocarcinoma (LUAD) is still lacking. The present study aimed to systematically explore the role of m⁶A methylation regulators in LUAD. RNA sequencing data of 20 m⁶A methylation regulators and clinical data of LUAD patients were downloaded from The Cancer Genome Atlas (TCGA) database. The prognosis value of m⁶A methylation regulators in LUAD was evaluated using the Gene Expression Profiling Interactive Analysis (GEPIA) and PrognoScan database. The correlation between IGF2BP1 and immune infiltrates in LUAD was investigated via CIBERSORT and Tumor Immune Estimation Resource (TIMER). A total of 15 m⁶A modification regulators were significantly abnormally expressed in LUAD tissues. Survival analysis revealed that four genes (HNRNPC, HNRNPA2B1, IGF2BP1, and IGF2BP3) were significantly associated with poor prognosis in LUAD. Multivariate Cox regression analysis showed that only IGF2BP1 was an independent predictor of LUAD after adjusting common clinical parameters. The mutation rates of m⁶A modification regulators in LUAD were less than 10%. Further analysis revealed that IGF2BP1 expression was significantly correlated with immune infiltration, the expression of immune checkpoints, and tumor mutational burden (TMB) in LUAD. Our findings suggest that IGF2BP1 is an independent predictor and related to immunotherapy response in LUAD, which maybe a potential novel biomarker for LUAD prognosis and the status of tumor immunity.

Prognostic implications of N6-methyladenosine RNA regulators in breast cancer

The significance of N⁶-methyladenosine (m6A) RNA modifications in the progression of breast cancer (BC) has been recognised. However, their potential role and mechanism of action in the tumour microenvironment (TME) and immune response has not been demonstrated. Thus, the role of m6A regulators and their downstream target gene components in BC remain to be explored. In this study, we used a series of bioinformatics methods and experiments to conduct exploratory research on the possible role of m6A regulators in BC. First, two regulatory modes of immune activation and inactivation were determined by tumour classification. The TME, immune cell infiltration, and gene set variation analysis results confirmed the reliability of this pattern. The prognostic model of the m6A regulator was established by the least absolute shrinkage and selection operator and univariate and multivariate Cox analyses, with the two regulators most closely related to survival verified by real-time quantitative reverse transcription polymerase chain reaction. Next, the prognostic m6A regulator identified in the model was crossed with the differential copy number of variant genes in invasive BC (IBC), and it was determined that YTHDF1 was a hub regulator. Subsequently, single-cell analysis revealed the expression patterns of m6A regulators in different IBC cell populations and found that YTHDF1 had significantly higher expression in immune-related IBC cells. Therefore, we selected the intersection of the BC differential expression gene set and the differential expression gene set of a cell line with knocked-down YTHDF1 in literature to identify downstream target genes of YTHDF1, in which we found IFI6, EIR, and SPTBN1. A polymerase chain reaction was conducted to verify the results. Finally, we confirmed the role of YTHDF1 as a potential prognostic biomarker through pan-cancer analysis. Furthermore, our findings revealed that YTHDF1 can serve as a new molecular marker for BC immunotherapy.

A novel 14-gene signature for overall survival in lung adenocarcinoma based on the Bayesian hierarchical Cox proportional hazards model

There have been few investigations of cancer prognosis models based on Bayesian hierarchical models. In this study, we used a novel Bayesian method to screen mRNAs and estimate the effects of mRNAs on the prognosis of patients with lung adenocarcinoma. Based on the identified mRNAs, we can build a prognostic model combining mRNAs and clinical features, allowing us to explore new molecules with the potential to predict the prognosis of lung adenocarcinoma. The mRNA data (n = 594) and clinical data (n = 470) for lung adenocarcinoma were obtained from the TCGA database. Gene set enrichment analysis (GSEA), univariate Cox proportional hazards regression, and the Bayesian hierarchical Cox proportional hazards model were used to explore the mRNAs related to the prognosis of lung adenocarcinoma. Multivariate Cox proportional hazard regression was used to identify independent markers. The prediction performance of the prognostic model was evaluated not only by the internal cross-validation but also by the external validation based on the GEO dataset (n = 437). With the Bayesian hierarchical Cox proportional hazards model, a 14-gene signature that included CPS1, CTPS2, DARS2, IGFBP3, MCM5, MCM7, NME4, NT5E, PLK1, POLR3G, PTTG1, SERPINB5, TXNRD1, and TYMS was established to predict overall survival in lung adenocarcinoma. Multivariate analysis demonstrated that the 14-gene signature (HR 3.960, 95% CI 2.710–5.786), T classification (T₁, reference; T₃, HR 1.925, 95% CI 1.104–3.355) and N classification (N₀, reference; N₁, HR 2.212, 95% CI 1.520–3.220; N₂, HR 2.260, 95% CI 1.499–3.409) were independent predictors. The C-index of the model was 0.733 and 0.735, respectively, after performing cross-validation and external validation, a nomogram was provided for better prediction in clinical application. Bayesian hierarchical Cox proportional hazards models can be used to integrate high-dimensional omics information into a prediction model for lung adenocarcinoma to improve the prognostic prediction and discover potential targets. This approach may be a powerful predictive tool for clinicians treating malignant tumours.

Genomics and Prognosis Analysis of N6-Methyladenosine Regulators in Lung Adenocarcinoma

Objective: N⁶-methyladenosine (m⁶A) modification is involved in modulating various biological processes in human cancers. But the implication of m⁶A modification in lung adenocarcinoma (LUAD) is still unclear. Hence, this study conducted a comprehensive analysis of the expression and clinical implication of m⁶A regulators in LUAD.

Methods: Consensus clustering analysis of 502 LUAD samples in the TCGA dataset was presented based on the expression profiles of 20 m⁶A regulators using ConsensusClusterPlus package. Overall survival (OS), activation of signaling pathways and tumor immunity (immune/stromal score, tumor purity, expression of HLA and immune checkpoints, and immune cell infiltration) were compared between m⁶A modification patterns. The m⁶A-related genes between patterns were identified and prognostic m⁶A-related genes were imported into LASSO-cox regression analysis. The m⁶A risk score was developed and its prognostic implication was evaluated and externally verified in the GSE30219 and GSE72094 dataset. Furthermore, a nomogram that contained independent prognostic indicators was established, followed by external verification.

Results: Two m⁶A modification patterns were clustered across LUAD based on the expression similarity of the m⁶A regulators via consensus clustering analysis, with distinct OS, activation of signaling pathways and tumor immunity. Totally, 213 m⁶A-related genes that were identified by comparing two patterns were significantly related to LUAD prognosis. By LASSO method, we constructed the m⁶A risk score that was a reliable and independent prognostic factor for LUAD. Patients with low m⁶A risk score displayed a prominent survival advantage. After incorporating independent clinical features, we developed the prognostic nomogram that exhibited high predictive accuracy and the best clinical net benefit for OS.

Conclusion: Collectively, our study may provide a clinically useful tool for precise prognostic management and optimization of immunotherapeutic strategies for LUAD patients.

The roles and mechanisms of the m6A reader protein YTHDF1 in tumor biology and human diseases

YTHDF1 is the most versatile and powerful reader protein of N⁶-methyladenosine (m6A)-modified RNA, and it can recognize both G(m6A)C and A(m6A)C RNAs as ligands without sequence selectivity. YTHDF1 regulates target gene expression by different mechanisms, such as promoting translation or regulating the stability of mRNA. Numerous studies have shown that YTHDF1 plays an important role in tumor biology and nontumor lesions by mediating the protein translation of important genes or by affecting the expression of key factors involved in many important cell signaling pathways. Therefore, in this review we focus on some of the roles of YTHDF1 in tumor biology and diseases.

Survival-associated N6-adenosine methyltransferase signatures in lung squamous cell carcinoma and clinical verification

Background

N⁶-methyladenine (m⁶A) is the most common modification of mRNA and IncRNA in higher organisms. m⁶A has been confirmed to be related to the formation and progression of tumors and m⁶A-related genes can be used as prognostic biomarkers in a variety of tumors. However, there have been no similar studies on lung squamous cell carcinoma. The main purpose of this study was aimed to explore the differential expression of m⁶A-related genes in lung squamous cell carcinoma tissues and its relationship with patient clinical prognosis.

Methods

We integrated three m⁶A writers that catalyze the methylation of adenine on mRNA molecules. The training set including 501 patients with LUSC was collected from The Cancer Genome Atlas (TCGA) database and the test set including 181 patients with LUSC was collected from the Gene Expression Omnibus (GEO) database. Based on the expression level of the m⁶A methylase gene, we established a tumor subgroup and risk-prognosis model to quantify the risk index and long-term patient prognosis, which were confirmed by principal component analysis (PCA) and receiver operating characteristic (ROC) curve analysis. After lung squamous cell carcinoma tissue specimens were obtained during surgery, immunohistochemistry (IHC) was used to verify the results in vitro.

Results

The results of the study showed that the expression of the three m⁶A methylases in tumor tissues and normal tissues was significantly different (P < 0.05). The survival-prognostic model based on METTL3 gene expression showed better predictive performance (AUC: 0.706). Patients in the high-risk and low-risk groups exhibited significant differences in terms of survival time and 5-year and 10-year survival rates. Immunohistochemistry revealed that patients with high METTL3 expression exhibited a longer survival time than those with low METTL3 expression.

Conclusions

Our study showed that the molecular phenotype based on the expression of METTL3 may be an independent risk factor affecting the prognosis of lung squamous cell carcinoma. These findings not only prove the important role of m⁶A methylase in lung squamous cell carcinoma, but are also expected to provide more accurate prognostic assessment and individualized treatment for patients with lung squamous cell carcinoma.

m6A regulator expression profile predicts the prognosis, benefit of adjuvant chemotherapy, and response to anti-PD-1 immunotherapy in patients with small-cell lung cancer

BACKGROUND

Small cell lung cancer (SCLC) is lethal and possesses limited therapeutic options. Platinum-based chemotherapy-with or without immune checkpoint inhibitors (anti-PDs)-is the current first-line therapy for SCLCs; however, its associated outcomes are heterogeneous. N⁶-methyladenosine (m⁶A) is a novel and decisive factor in tumour progression, chemotherapy resistance, and immunotherapy response. However, m⁶A modification in SCLC remains poorly understood.

METHODS

We systematically explored the molecular features and clinical significance of m⁶A regulators in SCLC. We then constructed an m⁶A regulator-based prognostic signature (m⁶A score) based on our examination of 256 cases with limited-stage SCLC (LS-SCLC) from three different cohorts-including an independent cohort that contained 150 cases with qPCR data. We additionally evaluated the relationships between the m⁶A score and adjuvant chemotherapy (ACT) benefits and the patients' responses to anti-PD-1 treatment. Immunohistochemical (IHC) staining and the HALO digital pathological platform were used to calculate CD8+ T cell density.

RESULTS

We observed abnormal somatic mutations and expressions of m⁶A regulators. Using the LASSO Cox model, a five-regulator-based (G3BP1, METTL5, ALKBH5, IGF2BP3, and RBM15B) m⁶A score was generated from the significant regulators to classify patients into high- and low-score groups. In the training cohort, patients with high scores had shorter overall survival (HR, 5.19; 2.75-9.77; P < 0.001). The prognostic accuracy of the m⁶A score was well validated in two independent cohorts (HR 4.6, P = 0.006 and HR 3.07, P < 0.001). Time-dependent ROC and C-index analyses found the m⁶A score to possess superior predictive power than other clinicopathological parameters. A multicentre multivariate analysis revealed the m⁶A score to be an independent prognostic indicator. Additionally, patients with low scores received a greater survival benefit from ACT, exhibited more CD8+ T cell infiltration, and were more responsive to cancer immunotherapy.

CONCLUSIONS

Our results, for the first time, affirm the significance of m⁶A regulators in LS-SCLC. Our multicentre analysis found that the m⁶A score was a reliable prognostic tool for guiding chemotherapy and immunotherapy selections for patients with SCLC.

An integrated model of FTO and METTL3 expression that predicts prognosis in lung squamous cell carcinoma patients

Background

Lung squamous cell carcinoma (LUSC) approximately accounts for a third of lung cancers. However, the role of N6-methyladenosine (m6A) in LUSC remains largely unknown according to previous studies.

Methods

In this study, we investigated the mutations, copy number variants (CNVs), expression of 20 m6A RNA methylation regulators, and clinical data from The Cancer Genome Atlas-LUSC (TCGA-LUSC). These data were used for the training cohort of screening potential biomarkers. The prognostic model of m6A RNA methylation regulators was constructed. A receiver operating characteristic (ROC) analysis was undertaken to determine the area under the curves (AUCs) (for 3- and 5-year survival) for the model. Additionally, the accuracy of the two-gene model was confirmed with external data verifications. Combined two-gene model and clinincal information were performed to construct a nomogram to predict patient’s prognostic risk assessment.

Results

Fat mass- and obesity-associated protein (FTO) and methyltransferase-like 3 (METTL3) were identified as potential prognostic biomarkers to evaluate benign and malignant tumors and prognosticate. The following prognostic model of m6A RNA methylation regulators was constructed: risk score = 0.162 × FTO − 0.069 × METTL3. Patients in low-risk group [median overall survival (mOS), 43.4 months] had longer survival than those with high-risk (mOS, 67.3 months) with P=0.0023. The smoking grade and risk score could be independent prognostic factors (P=0.00098 and P=0.0014, respectively). Ultimately, a nomogram was developed to assist clinicians to predict clinical outcomes.

Conclusions

FTO and METTL3 are potential prognostic biomarkers of LUSC. The two-gene model’s use of prognostic risk scores may provide guidance in the selection of therapeutic strategies.

LncRNA LINC00337 sponges mir-1285-3p to promote proliferation and metastasis of lung adenocarcinoma cells by upregulating YTHDF1

BACKGROUND

Emerging studies have shown that long noncoding RNAs (lncRNAs) predominantly function in the carcinogenesis of multiple developing human tumors. The current study aimed to investigate the underlying mechanisms of LINC00337 in lung adenocarcinoma.

METHODS

We analyzed TCGA and GTEx datasets and chose LINC00337 as the research object. Cell proliferation, cell apoptosis, cell cycle, migration, and invasion were detected in the gain and loss experiments of LINC00337 both in vitro and in vivo. Moreover, RNA pull-down, luciferase reporter assays, western blotting analysis, and rescue experiments were performed to investigate the underlying molecular mechanisms of LINC00337 function.

RESULTS

LINC00337 expression was remarkably upregulated in lung adenocarcinoma. In addition, LINC00337 knockdown was shown to repress cell migration, invasion, and proliferation, as well as the cell cycle, and gear up apoptosis in lung adenocarcinoma in vitro and in vivo. With respect to the mechanism, LINC00337 knockdown boosted miR-1285-3p expression and then restrained YTHDF1 expression post-transcriptionally. Crucially, both miR-1285-3p decrement and YTHDF1 overexpression successfully reversed the influence on cell proliferation, migration, invasion, and apoptosis caused by LINC00337 shRNA.

CONCLUSIONS

These results suggest that LINC00337 acts as an oncogenic lncRNA, targeting miR-1285-3p and regulating YTHDF1 expression, to promote the progression of lung adenocarcinoma.

Development and validation of a novel N6-methyladenosine (m6A)-related multi- long non-coding RNA (lncRNA) prognostic signature in pancreatic adenocarcinoma

Accumulating evidence has unveiled the pivotal roles of N6-methyladenosine (m6A) in pancreatic adenocarcinoma (PAAD). However, there are not many researches to predict the prognosis of PAAD using m6A-related long non-coding RNAs (lncRNAs). Raw data from The Cancer Genome Atlas (TCGA), International Cancer Genome Consortium (ICGC), and the Genotype-Tissue Expression project (GTEx) were utilized to comprehensively analyze the expression and prognostic performances of 145 m6A-related lncRNAs in PAAD and to develop and validate a novel m6A-related multi-lncRNA prognostic signature (m6A-LPS) for PAAD patients. In total, 57 differentially expressed m6A-related lncRNAs with prognostic values were identified. Based on LASSO-Cox regression analysis, m6A-LPS was constructed and verified by using five-lncRNA expression profiles for TCGA and ICGC cohorts. PAAD patients were then divided into high- and low-risKBIE_A_1933868k subgroups with different clinical outcomes according to the median risk score; this was further verified by time-dependent receiver operating characteristic curves. Risk scores were significantly associated with clinical parameters such as histological grade and cancer status among PAAD patients. A nomogram consisting of risk score, grade, and cancer status was generated to predict the survival probability of PAAD patients, as also demonstrated by calibration curves. Discrepancies in cellular processes, signaling pathways, and immune status between the high- and low-risk subgroups were investigated by functional and single-sample gene set enrichment analyses. In conclusion, the novel m6A-LPS for PAAD patients was developed and validated, which might provide new insight into clinical decision-making and precision medicine.

Identification of the Signature Associated With m6A RNA Methylation Regulators and m6A-Related Genes and Construction of the Risk Score for Prognostication in Early-Stage Lung Adenocarcinoma

Background

N6-methyladenosine (m⁶A) RNA modification is vital for cancers because methylation can alter gene expression and even affect some functional modification. Our study aimed to analyze m⁶A RNA methylation regulators and m⁶A-related genes to understand the prognosis of early lung adenocarcinoma.

Methods

The relevant datasets were utilized to analyze 21 m⁶A RNA methylation regulators and 5,486 m⁶A-related genes in m⁶Avar. Univariate Cox regression analysis, random survival forest analysis, Kaplan-Meier analysis, Chi-square analysis, and multivariate cox analysis were carried out on the datasets, and a risk prognostic model based on three feature genes was constructed.

Results

Respectively, we treated GSE31210 (n = 226) as the training set, GSE50081 (n = 128) and TCGA data (n = 400) as the test set. By performing univariable cox regression analysis and random survival forest algorithm in the training group, 218 genes were significant and three prognosis-related genes (ZCRB1, ADH1C, and YTHDC2) were screened out, which could divide LUAD patients into low and high-risk group (P < 0.0001). The predictive efficacy of the model was confirmed in the test group GSE50081 (P = 0.0018) and the TCGA datasets (P = 0.014). Multivariable cox manifested that the three-gene signature was an independent risk factor in LUAD. Furthermore, genes in the signature were also externally validated using the online database. Moreover, YTHDC2 was the important gene in the risk score model and played a vital role in readers of m⁶A methylation.

Conclusion

The findings of this study suggested that associated with m⁶A RNA methylation regulators and m⁶A-related genes, the three-gene signature was a reliable prognostic indicator for LUAD patients, indicating a clinical application prospect to serve as a potential therapeutic target.

Pan-Cancer Prognostic, Immunity, Stemness, and Anticancer Drug Sensitivity Characterization of N6-Methyladenosine RNA Modification Regulators in Human Cancers

N6-methyladenosine RNA modification plays a significant role in the progression of multiple tumorigenesis. Our study identified the imperative role of m6A regulators in the tumor immune microenvironment, survival, stemness score, and anticancer drug sensitivity of pan-cancer. The Wilcox test was to identify the differential expression between 17 m6A regulators across 33 TCGA cancer types and their normal tissues from UCSC Xena GDC pan-cancer. Survival analysis of m6A-related regulators in 33 TCGA cancer types was identified using the "survival" and "survminer" package. The Spearman correlation test and Pearson correlation test were used to identify the correlation relationship between m6A regulators expression and tumor microenvironment, tumor stem cell score, and drug sensitivity of anticancer drugs. ConsensusPathDB was used for exploring m6A regulators functional enrichment. The 17 (METTL3, WTAP, METTL14, RBM15, RBM15B, VIRMA, HNRNPC, HNRNPA2B1, YTHDC1, ZC3H13, YTHDF1, YTHDC2, YTHDF2, IGF2BP3, IGF2BP1, FTO, and ALKBH5) m6A regulators were differentially expressed in 18 TCGA cancer types and adjacent normal tissues. Correlation analysis indicated that the relationship between the expression of 17 m6A regulators and tumor microenvironment indicated that the higher expression of m6A regulators, the higher the degree of tumor stem cells. The anticancer drug sensitivity analysis indicated that ZC3H13 expression had a positive relationship with anticancer drugs such as selumetinib, dabrafenib, cobimetinib, trametinib, and hypothemycin (p < 0.001). YTHDF2 expression was significantly negatively correlated with the anticancer drug dasatinib (p < 0.001). The pan-cancer immune subtype analysis showed that the 17 m6A regulators were significantly different in immune subtype C1 (wound healing), C3 (inflammatory), C2 (IFN-gamma dominant), C5 (immunological quiet), C4 (lymphocyte depleted), and C6 (TGF-beta dominant) (p < 0.001). Our study provides a comprehensive insight for revealing the significant role of m6A regulators in the tumor immune microenvironment, stemness score, and anticancer drug sensitivity of human cancers.

Importance of N6-methyladenosine RNA modification in lung cancer (Review)

The N6-methyladenosine (m6A) modification is the most common mRNA modification in eukaryotes and exerts biological functions by affecting RNA metabolism. The m6A modification is installed by m6A methyltransferases, removed by demethylases and recognized by m6A-binding proteins. The interaction between these three elements maintains the dynamic equilibrium of m6A in cells. Accumulating evidence indicates that m6A RNA methylation has a significant impact on RNA metabolism and is involved in the pathogenesis of cancer. Lung cancer is the leading cause of cancer-related deaths worldwide. The treatment options for lung cancer have developed considerably over the past few years; however, the survival rate of patients with lung cancer still remains very low. Although diagnostic methods and targeted therapies have been rapidly developed in recent years, the underlying mechanism and importance of m6A RNA methylation in the pathogenesis of lung cancer remains ambiguous. The current review summarized the biological functions of m6A modification and considers the potential roles of m6A regulators in the occurrence and development of lung cancer.

The Role of RNA Modifications and RNA-modifying Proteins in Cancer Therapy and Drug Resistance

The advent of new genome-wide sequencing technologies has uncovered abnormal RNA modifications and RNA editing in a variety of human cancers. The discovery of reversible RNA N6-methyladenosine (RNA: m⁶A) by fat mass and obesity-associated protein (FTO) demethylase has led to exponential publications on the pathophysiological functions of m⁶A and its corresponding RNA modifying proteins (RMPs) in the past decade. Some excellent reviews have summarized the recent progress in this field. Compared to the extent of research into RNA: m⁶A and DNA 5-methylcytosine (DNA: m⁵C), much less is known about other RNA modifications and their associated RMPs, such as the role of RNA: m⁵C and its RNA cytosine methyltransferases (RCMTs) in cancer therapy and drug resistance. In this review, we will summarize the recent progress surrounding the function, intramolecular distribution and subcellular localization of several major RNA modifications, including 5' cap N7-methylguanosine (m7G) and 2'-O-methylation (Nm), m⁶A, m⁵C, A-to-I editing, and the associated RMPs. We will then discuss dysregulation of those RNA modifications and RMPs in cancer and their role in cancer therapy and drug resistance.

Suppression of Heterogeneous Nuclear Ribonucleoprotein C Inhibit Hepatocellular Carcinoma Proliferation, Migration, and Invasion via Ras/MAPK Signaling Pathway

Hepatocellular carcinoma (HCC), the most common malignant tumor, has high fatality and recurrence rates. Accumulating evidence shows that heterogeneous nuclear ribonucleoprotein C (HNRNPC), which is mainly involved in RNA splicing, export, and translation, promotes progression and metastasis of multiple tumor types; however, the effects of HNRNPC in HCC are unknown. In the present study, high levels of HNRNPC were detected in tumor tissues compared with para-tumor tissues by immunohistochemical and western blot assays. Furthermore, Cox proportional hazards regression models, the Kaplan–Meier method, and clinicopathologic features analysis showed that HNRNPC was not only an independent prognostic factor for both overall and disease-free survival in HCC but also a predictor of large tumor size and advanced tumor stage. Functional experiments revealed that silencing of HNRNPC not only led to arrest of more HCC cells at G0/G1 phase to inhibit their proliferation, but also suppressed EMT process to block their invasion, and migration in vitro; this was related to the Ras/MAPK signaling pathway. In addition, blocking of HCC cell proliferation regulated by HNRNPC silencing was observed in vivo. Finally, rescue tests showed that after recovery of Ras/MAPK signaling pathway activity by treatment with Ras agonists, the proliferation, migration, and invasion suppression of Huh-7 and Hep 3B cell lines caused by HNRNPC knockdown was partially reversed. Taken together, these results indicate that HNRNPC knockdown inhibits HCC cell proliferation, migration and invasion, in part via the Ras/MAPK signaling pathway. Thus, HNRNPC may have an important role in the progression of HCC and represents a promising biomarker for evaluation of prognosis and a potential therapeutic target in HCC patients.

A five-m6A regulatory gene signature is a prognostic biomarker in lung adenocarcinoma patients

We analyzed the prognostic value of N6-methyladenosine (m6A) regulatory genes in lung adenocarcinoma (LADC) and their association with tumor immunity and immunotherapy response. Seventeen of 20 m6A regulatory genes were differentially expressed in LDAC tissue samples from the TCGA and GEO databases. We developed a five-m6A regulatory gene prognostic signature based on univariate and Lasso Cox regression analysis. Western blot analysis confirmed that the five prognostic m6A regulatory proteins were highly expressed in LADC tissues. We constructed a nomogram with five-m6A regulatory gene prognostic risk signature and AJCC stages. ROC curves and calibration curves showed that the nomogram was well calibrated and accurately distinguished high-risk and low-risk LADC patients. Weighted gene co-expression analysis showed significant correlation between prognostic risk signature genes and the turquoise module enriched with cell cycle genes. The high-risk LADC patients showed significantly higher PD-L1 levels, increased tumor mutational burden, and a lower proportion of CD8+ T cells in the tumor tissues and improved response to immune checkpoint blockade therapy. These findings show that this five-m6A regulatory gene signature is a prognostic biomarker in LADC and that immune checkpoint blockade is a potential therapeutic option for high-risk LADC patients.

The Important Role of N6-methyladenosine RNA Modification in Non-Small Cell Lung Cancer

N6-methyladenosine (m⁶A) is one of the most prevalent epigenetic modifications of eukaryotic RNA. The m⁶A modification is a dynamic and reversible process, regulated by three kinds of regulator, including m⁶A methyltransferases, demethylases and m⁶A-binding proteins, and this modification plays a vital role in many diseases, especially in cancers. Accumulated evidence has proven that this modification has a significant effect on cellular biological functions and cancer progression; however, little is known about the effects of the m⁶A modification in non-small cell lung cancer (NSCLC). In this review, we summarized how various m⁶A regulators modulate m⁶A RNA metabolism and demonstrated the effect of m⁶A modification on the progression and cellular biological functions of NSCLC. We also discussed how m⁶A modification affects the treatment, drug resistance, diagnosis and prognosis of NSCLC patients.

m6A RNA Methylation Regulators Act as Potential Prognostic Biomarkers in Lung Adenocarcinoma

N⁶-methyladenosine [m(6)A/m⁶A] methylation is one of the most common RNA modifications in eukaryotic cell mRNA and plays an important regulatory role in mRNA metabolism, splicing, translocation, stability, and translation. Previous studies have demonstrated that the m⁶A modification is highly associated with tumor cell proliferation, migration, and invasion. In the present study, five m⁶A regulatory factors have been revealed, namely heterogeneous nuclear ribonucleoprotein A2/B1(HNRNPA2B1), heterogeneous nuclear ribonucleoprotein C (HNRNPC), Vir like m⁶A methyltransferase associated protein (KIAA1429/VIRMA), RNA binding motif protein 15 (RBM15) and methyltransferase like 3 (METTL3), which are closely related to the overall survival (OS) of patients with lung adenocarcinoma (LUAD). These five m⁶A regulatory factors exhibited potential prognostic value for the 1, 3, and 5-years survival outcomes of LUAD patients. Our findings revealed that several signaling pathways, such as cell cycle, DNA replication, RNA degradation, RNA polymerase, nucleotide excision repair and basal transcription factors, are activated in the high-risk group of LUAD patients.

Gene signatures of 6-methyladenine regulators in women with lung adenocarcinoma and development of a risk scoring system: a retrospective study using the cancer genome atlas database

Although the emergence of new treatments has improved the prognosis of women with lung adenocarcinoma (LUAD), the emergence of drug resistance limits their clinical efficacy. Therefore, there is an urgent need to identify new targets and develop a risk scoring system to evaluate the prognosis of patients. 6-methyladenine (M6A), as the most common methyl modification in RNA modification, its clinicopathological features, diagnosis and prognostic value in lung cancer, especially in LUAD remain to be discussed. We analyzed the clinical and sequencing data of the female LUAD cohort from The Cancer Genome Atlas (TCGA), evaluated the expression profiles of 16 M6A regulation-related genes in the cohort and the relationships between genetic changes and clinical characteristics, developed an M6A-related risk scoring system using Cox analysis. Finally, the copy number variations (CNVs) of the related genes in the samples were analyzed and verified using the cBioPortal platform. Compared with other clinical factors, this risk scoring system showed a higher predictive sensitivity and specificity. The M6A-related risk scoring system developed in this study may help to improve the screening of female patients at high risk of LUAD and provides important theoretical bioinformatics support for evaluating the prognosis of such patients.

A Novel RNA Binding Protein-Related Prognostic Signature for Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is a highly malignant and aggressive cancer with high recurrence rates and mortality. Some studies have illustrated that RNA binding proteins (RBPs) were involved in the carcinogenesis and development of multiple cancers, but the roles in HCC were still unclear. We downloaded the RNA-seq and corresponding clinical information of HCC from The Cancer Genome Atlas (TCGA) database, and 330 differentially expressed RBPs were identified between normal and HCC tissues. Through series of the univariate, the least absolute shrinkage selection operator (LASSO), and the stepwise multivariate Cox regression analyses, six prognosis-related key RBPs (CNOT6, UPF3B, MRPL54, ZC3H13, IFIT5, and PPARGC1A) were screened out from DE RBPs, and a six-RBP gene risk score signature was constructed in training set. Survival analysis indicated that HCC patients with high-risk scores had significantly worse overall survival than low-risk patients, and furthermore, the signature can be used as an independent prognostic indicator. The good accuracy of this prognostic signature was confirmed by the ROC curve analysis and was further validated in the International Cancer Genome Consortium (ICGC) HCC cohort. Besides, a nomogram based on six RBP genes was established and internally validated in the TCGA cohort. Gene set enrichment analysis demonstrated some cancer-related phenotypes were significantly gathered in the high-risk group. Overall, our study first identified an RBP-related six-gene prognostic signature, which could serve as a promising prognostic biomarker and provide some potential therapeutic targets for HCC.

Molecular characterization, biological function, tumor microenvironment association and clinical significance of m6A regulators in lung adenocarcinoma

N6-methyladenosine (m6A) modification can regulate a variety of biological processes. However, the implications of m6A modification in lung adenocarcinoma (LUAD) remain largely unknown. Here, we systematically evaluated the m6A modification features in more than 2400 LUAD samples by analyzing the multi-omics features of 23 m6A regulators. We depicted the genetic variation features of m6A regulators, and found mutations of FTO and YTHDF3 were linked to worse overall survival. Many m6A regulators were aberrantly expressed in tumors, among which FTO, IGF2BP3, YTHDF1 and RBM15 showed consistent alteration features across 11 independent cohorts. Besides, the regulator-pathway interaction network demonstrated that m6A modification was associated with various biological pathways, including immune-related pathways. The correlation between m6A regulators and tumor microenvironment was also assessed. We found that LRPPRC was negatively correlated with most tumor-infiltrating immune cells. On the other hand, we established a scoring tool named m6Sig, which was positively correlated with PD-L1 expression and could reflect both the tumor microenvironment characterization and prognosis of LUAD patients. Comparison of CNV between high and low m6Sig groups revealed differences on chromosome 7. Application of m6Sig on an anti-PD-L1 immunotherapy cohort confirmed that the high m6Sig group demonstrated therapeutic advantages and clinical benefits. Our study indicated that m6A modification is involved in many aspects of LUAD and contributes to tumor microenvironment formation. A better understanding of m6A modification will provide more insights into the molecular mechanisms of LUAD and facilitate developing more effective personalized treatment strategies. A web application was built along with this study (http://www.bioinfo-zs.com/luadexpress/).

N6-methyladenosine (m6A) RNA modification in cancer stem cells

Cancer stem cells (CSCs), a unique subset of undifferentiated cells with stem cell-like properties, have emerged as driving forces in mediating tumor growth, metastasis, and therapeutic resistance. Recent advances have highlighted that N6-methyladenosine (m6A) RNA modification plays an important role in cancer biology and CSCs. Dynamic m6A decoration has been demonstrated to be involved in CSC generation and maintenance, governing cancer progression and therapeutic resistance. In this review, we provide the first overview of the current knowledge of m6A modification implicated in CSCs and their impact on CSC properties, tumor progression, and responses to treatment. Finally, we also highlight the potential of m6A machinery as novel targets for cancer therapeutics. The involvement of m6A modification in CSCs provides a new direction for exploring cancer pathogenesis and inspires the development of effective strategies to fully eliminate both cancer cells and CSCs.

Downregulation of m6 A reader YTHDC2 promotes tumor progression and predicts poor prognosis in non-small cell lung cancer

Background

m⁶A modification affects the pathological progress of many diseases by affecting RNA stability and translocation. YTHDC2, a m⁶A reader, is associated with multiple cancers; however, little is known of its role in non‐small cell lung cancer (NSCLC).

Methods

The GEPIA, Oncomine and GEO databases were analyzed to assess expression of YTHDC2 in NSCLC patients. Quantitative polymerase chain reaction, western blot and immunohistochemistry were used to detect YTHDC2 expression in different NSCLC cell lines (H1299, H460, H292 and A549) and patients. The effects of YTHDC2 on NSCLC cell lines (A549 and H1299) proliferation and migration were employed using CCK8 and transwell assays. The potential target RNAs of YTHDC2 were obtained from the POSTAR database. Functional enrichment analysis of YTHDC2 targeted RNAs was performed using the Metascape database.

Results

GEPIA, Oncomine and GEO databases showed low expression of YTHDC2 in lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) patients. YTHDC2 expression was significantly decreased in different NSCLC cell lines and our clinical samples. Moreover, low expression of YTHDC2 was significantly associated with poor differentiation, lymph node metastasis, tumor size and stage. In addition, YTHDC2 could suppress the proliferation and migration ability of A549 and H1299 cell lines. Kaplan‐Meier Plotter database analysis revealed that patients with low level of YTHDC2 had a significantly poor prognosis. Finally, functional enrichment analysis of YTHDC2 targeted RNAs indicated several enriched pathways related to cancer.

Conclusions

These findings elucidate that YTHDC2 suppresses tumorigenesis in NSCLC, indicating that YTHDC2 may be a promising therapeutic target for NSCLC.

Key points

Significant findings of the study

This study demonstrated that the downregulation of YTHDC2 promotes tumor progression and predicts poor prognosis in non‐small cell lung cancer (NSCLC).

What this study adds

YTHDC2 might be a promising therapeutic target for non‐small cell lung cancer (NSCLC).

Crosstalk between RNA m6A Modification and Non-coding RNA Contributes to Cancer Growth and Progression

N6-methyladenosine (m6A) is the most common RNA modification and has an important role in normal development and tumorigenesis. The abnormal expression of m6A regulators can lead to an imbalance in m6A levels in cancer cells, leading to the dysregulated expression of oncogenes and tumor suppressor genes that may contribute to cancer development, patient response to chemoradiotherapy, and clinical prognosis. Recent studies demonstrate that non-coding RNAs are involved in epigenetic modification of both DNA and RNA in tumor cells, and may also affect the development and progression of cancer by targeting m6A regulators. In this review, we describe the functional crosstalk between m6A and non-coding RNAs, particularly microRNA, long non-coding RNA, and circular RNA, and illustrate their roles in tumor regulation. Finally, we discuss the significance of non-coding RNA and m6A modification in the diagnosis, treatment, and prognosis of cancer patients, as well as potential future research directions.