- Letter to the Editor
- Open Access
The landscape of m6A regulators in small cell lung cancer: molecular characteristics, immuno-oncology features, and clinical relevance
Molecular Cancer volume 20, Article number: 122 (2021)
Lung cancer continues to be the world’s most common and deadly cancer . Small cell lung cancer (SCLC) is the most lethal lung cancer and accounts for ~ 15% of total lung cancer cases. SCLC is an aggressive high-grade neuroendocrine tumor and characterized by a short doubling time, rapid growth, and early metastatic spread . Most patients with SCLC rapidly develop drug resistance and their 5-year survival is low (5–6%), even in cases with a good initial response to standard chemotherapy . The addition of immune checkpoint inhibitors to conventional chemotherapy for SCLCs is promising; however, their absolute prolonged benefit is moderate . The complex mechanisms underlying widespread SCLC metastases and recurrence require elucidation to extend the durable benefit of chemotherapy and immunotherapy to more patients.
N6-methyladenosine (m6A), the most abundant and prevalent RNA modification in eukaryotic RNA, is an important component in cancer biology . M6A modification—similar to other epigenetic modifications—is a dynamic reversible process regulated by methyltransferases, RNA binding proteins, and demethylases. Accumulating evidence strongly indicates that the m6A modification is a novel decisive factor for tumor metastases, recurrence, and therapy resistance, particularly anti-PD-1/PD-L1 monotherapies resistance . Multiple epigenetic abnormalities have been proved tightly associated with cancerous phenotype and aggressiveness in SCLC . However, as the most common RNA epigenetic modification, m6A has never been linked to SCLC progression.
This was the first characterization of the landscape of m6A regulators—including their molecular characteristics, immuno-oncology features, and clinical relevance. Our data highlight the importance of m6A regulators on cancer pathogenesis and shaping of the tumor immune microenvironment (TIME), and lay a rational foundation and buttress for developing therapeutic strategies for targeting the m6A modification in patients with SCLC.
Results and discussion
The genetic variation, expression patterns, and therapeutic potential of m6A regulators in SCLC
We curated a catalog of 30 m6A regulators, including 11 writers, two erasers, and 17 readers (Additional file 1: Table S1). The multiple functions of m6A modification in cancers are summarized in Fig. 1a. Somatic mutations of 30 m6A regulators were explored in 110 SCLC samples. Twenty-eight samples exhibited mutations of m6A regulators, with a mutation frequency of 25.5% (Fig. 1b). Of all the regulators, readers showed relatively high mutation frequencies and FMR1 showed the highest mutation frequency. In contrast, no erasers exhibited mutations. We also noted co-occurrent mutations between METTL3 and YTHDC2 and between IGF2BP2 and YTHDC2 (Additional file 2: Figure S1). Then, we investigated the copy number variations (CNVs) of these regulators within 53 SCLC cell lines from the Cancer Cell Line Encyclopedia (CCLE). As shown in Fig. 1c, CNV alterations of m6A regulators in SCLC were universal. Most readers (11/17) showed widespread CNV amplification, while all erasers had more loss frequency. The chromosomal CNV alteration mutations are displayed in Fig. 1d. As the genetic alterations of m6A regulators are common in SCLC, we should determine if these changes affect expression patterns. The panorama of m6A regulator expression in normal lung and SCLC samples were evaluated through principal component analysis, which revealed a dramatically different distribution pattern (Fig. 1e). The expression details of the regulators between normal and SCLC samples are displayed in Fig. 1f. Nearly all the writers and readers were significantly upregulated in SCLCs; nevertheless, eraser expression tended to decline, indicating that SCLC is associated with abundant m6A modification. IGF2BP3 has been found at significantly higher levels in multiple solid tumors , this research is the first to extend this phenomenon to SCLC. Combining the CNV results with the expression pattern of m6A regulators, we infer that CNV alterations may contribute to perturbations in m6A regulator expression.
To further investigate the relationship between m6A expression profile and SCLC metastasis, we compared the distributions of 30 regulators in 50 SCLC cell lines from CCLE with different pathological origins. Most of the writers and readers showed an upward distribution tendency in metastasis-derived cell lines compared to primary-derived cell lines (Additional file 2: Figure S2), indicating that m6A modification may help promote SCLC metastasis. Of note, four regulators—METTL3, HNRNPA2B1, HNRNPC, and RBMX—were significantly upregulated in metastasis-derived cell lines. Could these molecules be therapy targets for metastatic SCLC? We first sought to determine which of the four regulators acted as pan-cancer fitness genes (Additional file 2: Figure S3) . We found that HNRNPC and RBMX were essential genes almost in all different cancer types. Moreover, DepMap data also reminded that HNRNPC and RBMX were necessary for pan-cancer cells (Additional file 2: Figure S4). Therefore, we assumed that HNRNPC and RBMX may also affect the malignant biological behavior of SCLC. We carried out several in vitro and in vivo experiments using two metastasis-derived cell lines and found that HNRNPC and RBMX knockdown drastically inhibited SCLC proliferation and migration and promoted apoptosis (Fig. 1g, h, Additional file 1: Table S2, Additional file 2: Figures S5–S7). We also found the protein expressions of HNRNPC and RBMX were related to SCLC staging in clinical samples (Additional file 2: Figure S8).
Relationship between m6A regulators and cancer pathways and immune characteristics in SCLC
To determine how m6A regulators participated in SCLC, we explored correlations between m6A regulators and 50 cancer hallmark-related pathways. (Additional file 3: methods). Regulator expression was closely related to different activation statuses of multiple oncogenic pathways (Fig. 2a). With the highest number of connections with activated pathways—including G2M checkpoint, DNA repair, PI3K/AKT/mTOR pathway, et cetera—HNRNPC and METTL5 showed the closest relationship with cancer promotion in SCLC (Fig. 2a, Additional file 2: Figure S9). In contrast, HNRNPA2B1 was correlated with the most inhibited pathways, suggesting that HNRNPA2B1 may be protective for SCLC progression (Additional file 2: Figure S9). In vitro experiments affirmed that METTL5 was a pro-tumoral regulator, and HNRNPA2B1 had no significant effect on SCLC cells (Additional file 2: Figure S10). Of note, m6A regulators are well known for their collaboration in cancer progression. We therefore tested the correlation between these 30 regulators in SCLC (Fig. 2b). Several highly correlated relationships were observed between writers and readers, or readers and readers (Additional file 2: Figure S11, R > 0.5). These relationships indicated that the cooperation between the writers and readers may promote SCLC progression. Among these relationships, VIRMA and YTHDF3 showed the highest correlation coefficient (Fig. 2b, R = 0.81). Additionally, protein–protein interaction networks revealed that these regulators interacted with every single other frequently (Fig. 2c). Notably, writers formed more-frequent connections, which is in line with that the writers are working together as protein complex. Therefore, m6A regulator cross-talk and cooperation may participate in SCLC progression.
We also investigated the relationships between regulators and the immune phenotype. Results from seven clusters of metagenes indicated that most regulators were negatively associated with multiple inflammatory actives (Fig. 2d). For instance, ELAVL1 and HNRNPA2B1 showed consistent negative associations with HCK, LCK, MHC-I, STAT1, MHC_ II, Interferon, and Ig G. We also used the ESTIMATE and MCP-counter methods to explore the distributions of stromal and immune cells with different regulators expression patterns (Fig. 2e). Frequent negative correlations were found between the writers and readers with immune-stromal cells, indicating that high m6A modification is closely related to immune-desert phenotype in SCLC. Meanwhile, we analyzed the immune checkpoint profiles and found many negative correlations, particularly in writers and readers [especially METTL3, HNRNPA2B1, and ELAVL1 (Fig. 2f)]. These results suggest that m6A regulators may be closely related to the heterogeneity and complexity of TIME in SCLCs.
Clinical significance of the m6A regulators in SCLC
Observed connections between m6A regulators and immuno-oncological features prompted us to explore the clinical significance of these molecules in SCLC. Most regulators (21/30) exhibited clinical relevance (Additional file 2: Figures S12 and S13a, Additional file 1: Table S3), inspiring us to building a m6A regulator-based signature (m6A score) for prediction. We assembled three independent cohorts to construct the m6A score, including the International Cohort (n = 77), the Shanghai Cohort (n = 48), and the National Cancer Center (NCC) Cohort (n = 152) (Additional file 1: Table S4). A least absolute shrinkage and selection operator model was introduced and a 11-regulator based m6A score was established (Additional file 2: Figure S13b, c). The m6A score-high patients had significantly worse overall survival (OS) in the International Cohort (Fig. 2g, P < 0.001). We confirmed the repeatability and robustness of m6A score with the Shanghai cohort (Fig. 2h, P = 0.001). To determine if the m6A score could be used in clinical practice, we collected 152 SCLC samples with qPCR data from our institute (NCC Cohort) (Additional file 1: Table S5). As expected, the classifier significantly stratified patients into high- vs low-m6A score groups for OS (Fig. 2i, P < 0.001) and progression-free survival (PFS) (Additional file 2: Figure S14a, P = 0.001). The m6A score was well verified in different clinical subgroups (Additional file 2: Figure S15) and performed well in patients who received chemotherapy, regardless of cohort (Additional file 2: Figure S14b, c & d, P < 0.001). Importantly, after adjusting for clinical and pathologic characteristics, the m6A score remained an independent prognostic factor (Additional file 2: Figure S14e). To our best knowledge, the m6A score is the first multicenter molecular model for SCLC with a large total sample size. Large-scale investigations are rare because of challenges related to obtaining tumor specimens during standard-of-care treatment regimens. In summary, m6A regulators are promising prognostic factors and their future application may help determine the optimal treatment strategy for SCLC.
This study describes, for the first time, the abnormal expression patterns, specific immuno-oncology features, and clinical relevance of m6A regulators in SCLCs. Comprehensive evaluation of m6A regulators in SCLCs will help enhance our understanding of tumorigenesis and the remodeled TIME and allow for more effective chemo- and immune-therapeutic options.
Availability of data and materials
The genetic and expression data of m6A regulators in the International Cohort were retrieved from cBioPortal (https://www.cbioportal.org). The expression details of m6A regulators in small cell lung cancer cell lines were downloaded from CCLE (https://portals.broadinstitute.org/ccle). The data extracted from GSE40275 and GSE4824 can be downloaded from GEO dataset (https://www.ncbi.nlm.nih.gov/geo/). The other data used and analyzed during this study are available from the corresponding author on reasonable request.
Small cell lung cancer
Tumor immune microenvironment
Copy number variation
Cancer Cell Line Encyclopedia
National Cancer Center
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49. https://doi.org/10.3322/caac.21660.
Rudin CM, Brambilla E, Faivre-Finn C, Sage J. Small-cell lung cancer. Nat Rev Dis Primers. 2021;7(1):3.
Khan P, Siddiqui JA, Maurya SK, Lakshmanan I, Jain M, Ganti AK, Salgia R, Batra SK, Nasser MW. Epigenetic landscape of small cell lung cancer: small image of a giant recalcitrant disease. Semin Cancer Biol. 2020. https://doi.org/10.3322/caac.21660.
Iams WT, Porter J, Horn L. Immunotherapeutic approaches for small-cell lung cancer. Nat Rev Clin Oncol. 2020;17(5):300–12.
Huang H, Weng H, Chen J. m(6)A modification in coding and non-coding RNAs: roles and therapeutic implications in cancer. Cancer Cell. 2020;37(3):270–88.
Nombela P, Miguel-López B, Blanco S. The role of m(6)A, m(5)C and Ψ RNA modifications in cancer: novel therapeutic opportunities. Mol Cancer. 2021;20(1):18.
Li Y, Xiao J, Bai J, Tian Y, Qu Y, Chen X, Wang Q, Li X, Zhang Y, Xu J. Molecular characterization and clinical relevance of m(6)A regulators across 33 cancer types. Mol Cancer. 2019;18(1):137.
Behan FM, Iorio F, Picco G, Gonçalves E, Beaver CM, Migliardi G, Santos R, Rao Y, Sassi F, Pinnelli M, et al. Prioritization of cancer therapeutic targets using CRISPR-Cas9 screens. Nature. 2019;568(7753):511–6.
All authors would like to thank the specimen donors used in this study and research groups for the International Cohort, the Shanghai Cohort, CCLE, GSE40275, and GSE4824. All authors would also like to thank CapitalBio Technology for their kind help.
This work was supported by the CAMS Innovation Fund for Medical Sciences (2017-I2M-1–005, 2016-I2M-1–001), the National Key R&D Program of China (2016YFC1303201), the National Natural Science Foundation of China (81802299, 81502514), the Fundamental Research Funds for the Central Universities (3332018070), and the National Key Basic Research Development Plan (2018YFC1312105).
Ethics approval and consent to participate
Patient data we used were acquired by publicly available datasets that were collected with patients’ informed consent. The clinical samples used in this study was approved by the Ethics Committee of our institute (the National Cancer Center/Cancer Hospital of the Chinese Academy of Medical Sciences).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The descriptions of 30 m6A regulators enrolled in this study. Supplementary Table S2. Primer sequences of cell line samples for qPCR. Supplementary Table S3. The optimum cutoff survival analysis of 30 m6A regulators in the International Cohort. Supplementary Table S4. Clinical characteristics of the patients from multiple institutions. Supplementary Table S5. Primer sequences of FFPE samples for qPCR.
Co-occurrence of genetic alterations of the m6A regulators in small cell lung cancer. Supplementary Figure S2. The expression pattern of 30 m6A regulators between 18 primary- and 32 metastasis-derived small cell lung cancer cell lines from CCLE. Supplementary Figure S3. The essential gene proportion of METTL3 (a), HNRNPA2B1 (b), HNRNPC (c), and RBMX (d) in pan-cancer cell lines. Supplementary Figure S4. The landscape of therapeutic potential of HNRNPC and RBMX in pan-cancer cell lines from the DepMap Portal. Supplementary Figure S5. HNRNPC is upregulated in small cell lung cancer and promotes cell proliferation and inhibits cell apoptosis. Supplementary Figure S6. RBMX is upregulated in small cell lung cancer and promotes cell proliferation and inhibits cell apoptosis. Supplementary Figure S7. HNRNPC and RBMX promote SCLC cells metastasis in vivo. Supplementary Figure S8. The protein level of HNRNPC and RBMX are positively related to SCLC staging. Supplementary Figure S9. The number of m6A regulator-related pathways. Supplementary Figure S10. Experimental exploration of HNRNPA2B1 and METTL5 in SCLC. Supplementary Figure S11. Positive correlations between them6A regulators in small cell lung cancer. Supplementary Figure S12. Kaplan–Meier survival analysis of SCLCs grouped by the expression of m6A regulators in the International Cohort. Supplementary Figure S13. The clinical analyze of m6A regulators in small cell lung cancer. Supplementary Figure S14. The clinical significance of m6A regulators in small cell lung cancer. Supplementary Figure S15. The performance of the m6A score in different clinical subgroups from different cohorts.
Supplementary materials and methods.
About this article
Cite this article
Zhang, C., Zhang, Z., Zhang, Z. et al. The landscape of m6A regulators in small cell lung cancer: molecular characteristics, immuno-oncology features, and clinical relevance. Mol Cancer 20, 122 (2021). https://doi.org/10.1186/s12943-021-01408-5