Characteristic of molecular subtypes in lung adenocarcinoma based on m6A RNA methylation modification and immune microenvironment

Tumor microenvironment characteristics and prognostic role of m6A modification in lung squamous cell carcinoma

Background

It has recently been determined that N6-methyladenosine (m6A) RNA methylation regulators have prominent effects on several cancers. However, the potential role of m6A modification in lung squamous cell carcinoma (LUSC) remains unclear.

Methods

We evaluated the modification pattern of m6A and studied the biological function of m6A regulators in LUSC. Then, we constructed the m6Ascore to predict the prognosis of LUSC and analyzed the relationship between the m6Ascore and tumor mutation burden, immune cell infiltration, and immunotherapy.

Result

In the unsupervised consensus cluster analysis, three different m6Aclusters were identified, which correspond to an immune activation state, a moderate immune activation state, and an immune tolerance state. Forty-two genes related to the m6A phenotype were used to construct the m6Ascore; subsequently, multiple validations of the m6Ascore were carried out to determine the relationship between the score and immune cell infiltration and response to CTLA-4/PD-1 inhibitor treatment. Further analysis revealed that the m6Ascore could effectively predict the prognosis of LUSC and that the m6A phenotype-related genes, FAM162A and LOM4, might be potential biomarkers.

Conclusion

These findings highlight the potential role of m6A modification in the prognosis, TME, and immunotherapy of LUSC and have profound implications for developing more effective personalized treatment strategies for LUSC.

ALKBH5 promotes the development of lung adenocarcinoma by regulating the polarization of M2 macrophages through CDCA4

OBJECTIVE

Lung adenocarcinoma (LUAD) is the most common subtype of lung cancer, with high morbidity and mortality. N6-methyladenosine (m6A) is an important regulator of LUAD progression. Here, we investigated the potential biological functions of ALKBH5 (a m6A demethylated enzyme) and cell division cycle associated protein 4 (CDCA4) in the progression of LUAD.

METHODS

The expressions of CDCA4, METTL3, ALKBH5, FTO, YTHDC2 and YTHDC1 mRNA and proteins in LUAD and adjacent tissues, as well as NCI-H1299 and NCI-H157 cells were detected by RT-qPCR and western blot. Meanwhile, the role of ALKBH5 and CDCA4 in macrophage polarization was explored through tumor formation in Lewis lung carcinoma (LLC) mice and the co-culture system of NCI-H1299 and NCI-H157/THP-1 cells. Cell characterization was further analyzed. The expression of Ki-67 in tumor tissue was tested by immunohistochemistry. The scale of M1 and M2 macrophages was determined by flow cytometry.

RESULTS

CDCA4 was significantly overexpressed in NCI-H1299 and NCI-H157 cell lines compared with BEAS-2B cells. The fold enrichment of CDCA4 m6A level in the overexpression (oe)-METTL3 or short hairpin (sh)-ALKBH5 cells was enhanced. Overexpression of CDCA4 promoted the cell viability, proliferation and migration, and inhibited apoptosis, which was reversed by sh-ALKBH5 intervention. Overexpression of YTHDC2 (not YTHDC1) inhibited the effect of CDCA4 on sh-ALKBH5 cells. sh-CDCA4 inhibited tumor growth and weight of LLC cells in mice, and promoted M1/M2 ratio in LLC mice and NCI-H1299/THP-1 and NCI-H157/THP-1 co-culture systems. Oe-CDCA4 promoted the volume and weight of tumor and inhibited the M1/M2 ratio of tumor tissue in LLC mice, but was reversed by sh-ALKBH5 intervention.

CONCLUSION

m6A demethylase ALKBH5 promotes the development of LUAD through CDCA4 regulation of malignant characterization and M1/M2 macrophage polarization.

A novel disulfidptosis-related prognostic gene signature and experimental validation identify ACTN4 as a novel therapeutic target in lung adenocarcinoma

BACKGROUND

Lung adenocarcinoma (LUAD) is a prevalent form of malignancy globally. Disulfidptosis is novel programmed cell death pathway based on disulfide proteins, may have a positive impact on the development of LUAD treatment strategies.

OBJECTIVE

To investigate the impact of disulfidptosis-related genes (DRGs) on the prognosis of LUAD, developed a risk model to facilitate the diagnosis and prognostication of patients. We also explored ACTN4 (DRGs) as a new therapeutic biomarker for LUAD.

METHODS

We investigated the expression patterns of DRGs in both LUAD and noncancerous tissues. To assess the prognostic value of the DRGs, we developed risk models through univariate Cox analysis and lasso regression. The expression and function of ACTN4 was evaluated by qRT-PCR, immunohistochemistry and in vitro experiments. The TIMER examined the association between ACTN4 expression and immune infiltration in LUAD.

RESULTS

Ten differentially expressed DRGs were identified. And ACTN4 was identified as potential risk factors through univariate Cox regression analysis (P< 0.05). ACTN4 expression and riskscore were used to construct a risk model to predict overall survival in LUAD, and high-risk demonstrated a significantly higher mortality rate compared to the low-risk cohort. qRT-PCR and immunohistochemistry assays indicated ACTN4 was upregulated in LUAD, and the upregulation was associated with clinicopathologic features. In vitro experiments showed the knockdown of ACTN4 expression inhibited the proliferation in LUAD cells. The TIMER analysis demonstrated a correlation between the expression of ACTN4 and the infiltration of diverse immune cells. Elevated ACTN4 expression was associated with a reduction in memory B cell count. Additionally, the ACTN4 expression was associated with m6A modification genes.

CONCLUSIONS

Our study introduced a prognostic model based on DRGs, which could forecast the prognosis of patients with LUAD. The biomarker ACTN4 exhibits promise for the diagnosis and management of LUAD, given its correlation with tumor immune infiltration and m6A modification.

Identification of molecular subgroups and establishment of risk model based on the response to oxidative stress to predict overall survival of patients with lung adenocarcinoma

OBJECTIVE

Oxidative stress is associated with the occurrence and development of lung cancer. However, the specific association between lung cancer and oxidative stress is unclear. This study aimed to investigate the role of oxidative stress in the progression and prognosis of lung adenocarcinoma (LUAD).

METHODS

The gene expression profiles and corresponding clinical information were collected from GEO and TCGA databases. Differentially expressed oxidative stress-related genes (OSRGs) were identified between normal and tumor samples. Consensus clustering was applied to identify oxidative stress-related molecular subgroups. Functional enrichment analysis, GSEA, and GSVA were performed to investigate the potential mechanisms. xCell was used to assess the immune status of the subgroups. A risk model was developed by the LASSO algorithm and validated using TCGA-LUAD, GSE13213, and GSE30219 datasets.

RESULTS

A total of 40 differentially expressed OSRGs and two oxidative stress-associated subgroups were identified. Enrichment analysis revealed that cell cycle-, inflammation- and oxidative stress-related pathways varied significantly in the two subgroups. Furthermore, a risk model was developed and validated based on the OSRGs, and findings indicated that the risk model exhibits good prediction and diagnosis values for LUAD patients.

CONCLUSION

The risk model based on the oxidative stress could act as an effective prognostic tool for LUAD patients. Our findings provided novel genetic biomarkers for prognosis prediction and personalized clinical treatment for LUAD patients.

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.

Thyroid cancer risk prediction model using m6A RNA methylation regulators: integrated bioinformatics analysis and histological validation

BACKGROUND

Epigenetic reprogramming has been reported to play a critical role in the progression of thyroid cancer. RNA methylation accounts for more than 60% of all RNA modifications, and N6-methyladenosine (m6A) is the most common modification of RNAs in higher organisms. The purpose of this study was to explore the related modification mode of m6A regulators construction and its evaluation on the clinical prognosis and therapeutic effect of thyroid cancer.

METHODS

The levels of 23 m6A regulators in The Cancer Genome Atlas (TCGA) were analyzed. Differentially expressed genes (DEGs) and survival analysis were performed based on TCGA-THCA clinicopathological and follow-up information, and the mRNA levels of representative genes were verified using clinical thyroid cancer data. In order to detect the effects of m6A regulators and their DEGs, consensus cluster analysis was carried out, and the expression of different m6A scores in Tumor Mutation Burden (TMB) and immune double antibodies (PD-1 antibody and CTLA4 antibody) were evaluated to predict the correlation between m6A score and thyroid cancer tumor immunotherapy response.

RESULTS

Different expression patterns of m6A regulatory factors were detected in thyroid cancer tumors and normal tissues, and several prognoses related m6A genes were obtained. Two different m6A modification patterns were determined by consensus cluster analysis. Two different subgroups were established by screening overlapping DEGs between two m6A clusters, with cluster A having the best prognosis. According to the m6A score extracted from DEGs, thyroid cancer patients can be divided into high and low score subgroups. Patients with lower m6A score have longer survival time and better clinical features. The relationship between m6A score and Tumor Mutation Burden (TMB) and its correlation with the expression of PD-1 antibody and CTLA4 antibody proved that m6A score could be used as a potential predictor of the efficacy of immunotherapy in thyroid cancer patients.

CONCLUSIONS

We screened DEGs from cluster m6A and constructed a highly predictive model with prognostic value by dividing TCGA-THCA into two different clusters and performing m6A score analysis. This study will help clarify the overall impact of m6A modification patterns on thyroid cancer progression and formulate more effective immunotherapy strategies.

Aberrant Expression and Prognostic Potential of IL-37 in Human Lung Adenocarcinoma

Interleukin-37 (IL-37) is a relatively new IL-1 family cytokine that, due to its immunoregulatory properties, has lately gained increasing attention in basic and translational biomedical research. Emerging evidence supports the implication of this protein in any human disorder in which immune homeostasis is compromised, including cancer. The aim of this study was to explore the prognostic and/or diagnostic potential of IL-37 and its receptor SIGIRR (single immunoglobulin IL-1-related receptor) in human tumors. We utilized a series of bioinformatics tools and -omics datasets to unravel possible associations of IL-37 and SIGIRR expression levels and genetic aberrations with tumor development, histopathological parameters, distribution of tumor-infiltrating immune cells, and survival rates of patients. Our data revealed that amongst the 17 human malignancies investigated, IL-37 exhibits higher expression levels in tumors of lung adenocarcinoma (LUAD). Moreover, the expression profiles of IL-37 and SIGIRR are associated with LUAD development and tumor stage, whereas their high mRNA levels are favorable prognostic factors for the overall survival of patients. What is more, IL-37 correlates positively with a LUAD-associated transcriptomic signature, and its nucleotide changes and expression levels are linked with distinct infiltration patterns of certain cell subsets known to control LUAD anti-tumor immune responses. Our data indicate the potential value of IL-37 and its receptor SIGIRR to serve as biomarkers and/or immune-checkpoint therapeutic targets for LUAD patients. Further, the data highlight the urgent need for further exploration of this cytokine and the underlying pathogenetic mechanisms to fully elucidate its implication in LUAD development and progression.

The regulation and potential roles of m6A modifications in early embryonic development and immune tolerance at the maternal-fetal interface

The immune microenvironment at the maternal-fetal interface was determined by the crosstalk between the trophoblast and maternal-derived cells, which dynamically changed during the whole gestation. Trophoblasts act as innate immune cells and dialogue with maternal-derived cells to ensure early embryonic development, depending on the local immune microenvironment. Therefore, dysfunctions in trophoblasts and maternal decidual cells contribute to pregnancy complications, especially recurrent pregnancy loss in early pregnancy. Since many unknown regulatory factors still affect the complex immune status, exploring new potential aspects that could influence early pregnancy is essential. RNA methylation plays an important role in contributing to the transcriptional regulation of various cells. Sufficient studies have shown the crucial roles of N6-methyladenosine (m6A)- and m6A-associated- regulators in embryogenesis during implantation. They are also essential in regulating innate and adaptive immune cells and the immune response and shaping the local and systemic immune microenvironment. However, the function of m6A modifications at the maternal-fetal interface still lacks wide research. This review highlights the critical functions of m6A in early embryonic development, summarizes the reported research on m6A in regulating immune cells and tumor immune microenvironment, and identifies the potential value of m6A modifications in shaping trophoblasts, decidual immune cells, and the microenvironment at the maternal-fetal interface. The m6A modifications are more likely to contribute to embryogenesis, placentation and shape the immune microenvironment at the maternal-fetal interface. Uncovering these crucial regulatory mechanisms could provide novel therapeutic targets for RNA methylation in early pregnancy.

Characterization of Necroptosis-Related Molecular Subtypes and Therapeutic Response in Lung Adenocarcinoma

Lung adenocarcinoma (LUAD) is one of the most common malignant tumors with high morbidity and mortality and is usually associated with therapeutic resistance and poor prognosis because of individual biological heterogeneity. There is an unmet need to screen for reliable parameters, especially immunotherapy-related biomarkers to predict the patient’s outcomes. Necroptosis is a special caspase-independent form of necrotic cell death associated with the pathogenesis, progression, and prognosis of multiple tumors but the potential connection between necroptosis-related genes (NRGs) and LUAD still remains unclear. In this study, we expounded mutational and transcriptional alterations of 67 NRGs in 522 LUAD samples and proposed a consensus-clustering subtype of these patients into two cohorts with distinct immunological and clinical prognosis characteristics. Cluster B patients were associated with a better prognosis and characterized by relatively lower expression of NRGs, higher immune scores in the tumor microenvironment (TME), more mild clinical stages, and downregulated expression of immunotherapy checkpoints. Subsequently, the NRG score was further established to predict the overall survival (OS) of LUAD patients using univariate Cox, LASSO, and multivariate Cox regression analyses. The immunological characteristics and potential predictive capability of NRG scores were further validated by 583 LUAD patients in external datasets. In addition to better survival and immune-activated conditions, low-NRG-score cohorts exhibited a significant positive correlation with the mRNA stem index (mRNAsi) and tumor mutation burden (TMB) levels. Combined with classical clinical characteristics and NRG scores, we successfully defined a novel necroptosis-related nomogram to accurately predict the 1/3/5-year survival rate of individual LUAD patients, and the potential predictive capability was further estimated and validated in multiple test datasets with high AUC values. Integrated transcriptomic analysis helps us seek vital NRGs and supplements a novel clinical application of NRG scores in predicting the overall survival and therapeutic benefits for LUAD patients.

Construction of a Novel Prognostic Model in Lung Adenocarcinoma Based on 7-Methylguanosine-Related Gene Signatures

Increasing evidence has implicated the modification of 7-methylguanosine (m7G), a type of RNA modification, in tumor progression. However, no comprehensive analysis to date has summarized the predicted role of m7G-related gene signatures in lung adenocarcinoma (LUAD). Herein, we aimed to develop a novel prognostic model in LUAD based on m7G-related gene signatures. The LUAD transcriptome profiling data and corresponding clinical data were acquired from the Cancer Genome Atlas (TCGA) and two Gene Expression Omnibus datasets. After screening, we first obtained 29 m7G-related genes, most of which were upregulated in tumor tissues and negatively associated with overall survival (OS). According to the expression similarity of m7G-related genes, the combined samples from the TCGA-LUAD and GSE68465 datasets were further classified as two clusters that exhibit distinct OS rates and genetic heterogeneity. Then, we constructed a novel prognostic model involving four genes by using 130 differentially expressed genes among the two clusters. The combined samples were randomly divided into a training cohort and an internal validation cohort in a 1:1 ratio, and the GSE72094 dataset was used as an external validation cohort. The samples were divided into high- and low-risk groups. We demonstrated that a higher risk score was an independent negative prognostic factor and predicted poor OS. A nomogram was further constructed to better predict the survival of LUAD patients. Functional enrichment analyses indicated that cell cycle and DNA replication-related biological processes and pathways were enriched in the high-risk group. More importantly, the low-risk group had greater infiltration and enrichment of most immune cells, as well as higher ESTIMATE, immune, and stromal scores. In addition, the high-risk group had a lower TIDE score and higher expressions of most immune checkpoint-related genes. We finally noticed that patients in the high-risk group were more sensitive to chemotherapeutic agents commonly used in LUAD. In conclusion, we herein summarized for the first time the alterations and prognostic role of m7G-related genes in LUAD and then constructed a prognostic model based on m7G-related gene signatures that could accurately and stably predict survival and guide individualized treatment decision-making in LUAD patients.

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.

Integrative Analysis of m6A RNA Methylation Regulators and the Tumor Immune Microenvironment in Non-Small-Cell Lung Cancer

Background

Non-small-cell lung cancer (NSCLC) is a major component of lung cancer and is significantly correlated with poor prognosis. N6-methyladenosine (m6A) RNA methylation is closely related to the occurrence, progression, and prognosis of cancer. The potential biological functions and mechanisms of m6A RNA methylation in the immune microenvironment are still unclear.

Methods

We assessed m6A RNA methylation modification patterns in 1326 NSCLC patient samples based on 20 m6A regulators, linking these clusters to the tumor microenvironment and immune cell infiltration. The m6Ascore was created to quantify the m6A modification patterns of individual tumors. We then assessed the value of NSCLC patients in terms of clinical prognosis and immunotherapy response.

Results

According to different mRNA expression levels, two different m6A clusters were identified. m6A aggregation was significantly associated with clinical prognostic characteristics, the tumor microenvironment, and immune-related biological processes. Fifteen differential genes were screened based on these two m6A clusters, and to further investigate the mechanisms of action of these differential genes, they were subjected to unsupervised clustering analysis, which classified them into four different genomic isoforms. Prognostic analysis indicated that the survival advantage of the m6A gene cluster A modification mode was significantly prominent. We continued to construct the m6Ascore, which was used as a scoring tool to evaluate tumor typing, immunity, and prognosis. Patients with a low m6Ascore showed a significant survival advantage, and the group with a low m6Ascore had a better prognosis predicted by immunotherapy. The anti-PD-1/L1 immunotherapy cohort showed that a lower m6Ascore was associated with higher efficacy of immunotherapy.

Conclusions

The results suggest that m6A RNA methylation regulators make an important difference in the tumor immune microenvironment of patients with NSCLC. m6A gene characterization and the construction of the m6Ascore provide us with a richer understanding of m6A RNA methylation modification patterns in NSCLC patients and help to predict clinical prognosis and immunotherapeutic response.

Prognosis and Immunotherapy Response With a Novel Golgi Apparatus Signature-Based Formula in Lung Adenocarcinoma

The Golgi apparatus (GA) is a cellular organelle that participates in the packaging, modification, and transport of proteins and lipids from the endoplasmic reticulum to be further fabricated before being presented to other cellular components. Recent studies have demonstrated that GA facilitates numerous cellular processes in cancer development. Therefore, this study aimed to establish a novel lung adenocarcinoma (LUAD) risk evaluation model based on GA gene signatures. In this study, we used TCGA-LUAD (n = 500) as the training cohort and GSE50081 (n = 127), GSE68465 (442), and GSE72094 (398) as the validation cohorts. Two immunotherapy datasets (GSE135222 and GSE126044) were also obtained from a previous study. Based on machine algorithms and bioinformatics methods, a GA gene-related risk score (GARS) was established. We found that the GARS independently predicted the prognosis of LUAD patients and remained effective across stages IA to IIIA. Then, we identified that the GARS was highly correlated with mutations in P53 and TTN. Further, this study identified that GARS is related to multiple immune microenvironmental characteristics. Furthermore, we investigated GSE135222 and GSE126044 and found that a lower GARS may be indicative of an improved therapeutic effect of PD-1/PD-L1 therapy. We also found that high GARS may lead to a better response to multiple anticancer drugs. Finally, we established a nomogram to better guide clinical application. To our knowledge, this is the first study to demonstrate a novel GA signature-based risk score formula to predict clinical prognosis and guide the treatment of LUAD patients.

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.

N6-Methyladenosine RNA Modification in the Tumor Immune Microenvironment: Novel Implications for Immunotherapy

N6-methyladenosine (m⁶A) methylation is one of the most common modifications of RNA in eukaryotic cells, and is mainly regulated by m⁶A methyltransferases (writers), m⁶A demethylases (erasers), and m⁶A binding proteins (readers). Recently, accumulating evidence has shown that m⁶A methylation plays crucial roles in the regulation of the tumor immune microenvironment, greatly impacting the initiation, progression, and metastasis processes of various cancers. In this review we first briefly summarizes the m⁶A-related concepts and detection methods, and then describes in detail the associations of m⁶A methylation modification with various tumor immune components especially immune cells (e.g., regulatory T cells, dendritic cells, macrophages, and myeloid-derived suppressor cells) in a variety of cancers. We discuss the relationship between m⁶A methylation and cancer occurrence and development with the involvement of tumor immunity highlighted, suggesting novel markers and potential targets for molecular pathological diagnosis and immunotherapy of various cancers.