Abstract
Background
Thyroid cancer (TC) is the most common endocrine malignancy worldwide, and immunotherapy is a new cancer treatment that stimulates and enhances the natural ability of the immune system to fight cancer cells. The role of RNA N6-methyladenosine (m6A) related genes in these challenges has recently become a research hotspot, but he potential role of m6A modifications in tumor microenvironment (TME) cell infiltration remains unknown.
Purpose
There is growing evidence that m6A plays a critical role in the regulation of gene expression by participating in important biological processes. A comprehensive analysis of the m6A regulator-mediated infiltration characteristics of the TME will help advance the understanding of immune regulation in thyroid tumors.
Methods
This study assessed m6A modification modes in 510 thyroid cancer samples from the Cancer Genome Atlas (TCGA) databases according to a comprehensive set of 24 m6A regulators. In this study, we analyzed the biological characteristics and m6A methylation modification patterns. Based on this, we constructed m6A signatures and analyzed m6A modification features in tumor somatic mutations and TCGA molecular subtypes.
Results
These modification modes were systematically linked to TME cell infiltration signatures. m6A modification patterns were comprehensively assessed and correlated with immune cell infiltration features in the TME. An unsupervised clustering approach was applied and three distinct m6A modification subtypes and three m6A-associated gene subtypes were identified. Additionally, three distinct m6A methylation modification modes were identified in the thyroid cancer samples. The TME profiles of the identified genetic subtypes were strongly congruent with the immuno-heat and immuno-cold phenotypes.
Conclusions
The results revealed that m6A modifications play an integral role in the diversity and complexity of thyroid carcinomas. Evaluating the m6A modification patterns of individual tumors will create more efficient immunotherapeutic strategies. A comprehensive analysis of the role of TME in thyroid cancer provides a research idea for studying the effect of m6A epigenetics on thyroid tumors and their immune microenvironment.
Similar content being viewed by others
Availability of data and materials
The datasets analyzed during the current study are available in the TCGA repository (https:// portal. gdc. cancer. gov) and UCSC Xena (https://xena.ucsc.edu/) database.
Abbreviations
- TCGA:
-
The Cancer Genome Atlas
- TC:
-
Thyroid cancer
- PTC:
-
Papillary thyroid carcinoma
- m6A:
-
RNA N6-methyladenosine
- TME:
-
Tumor microenvironment
- ZC3H13:
-
Zinc finger CCCH-type containing 13
- WTAP:
-
Wilms' tumor 1 associated protein
- VIRMA:
-
Vir like m6A methyltransferase associated
- ALKBH5:
-
AlkB homolog 5, RNA demethylase
- RBM15:
-
RNA binding motif protein 15
- METTL3/14/16:
-
Methyltransferase 3/14/16
- RBM15B:
-
RNA binding motif protein 15B
- LRPPRC:
-
Leucine rich pentatricopeptide repeat containing
- FTO:
-
Fat mass and obesity-associated gene
- FMR1:
-
Fragile X messenger ribonucleoprotein 1
- RBMX:
-
RNA binding motif protein X-linked
- HNRNPC:
-
Heterogeneous nuclear ribonucleoprotein C
- EIF3A:
-
Eukaryotic translation initiation factor 3 subunit A
- YTHDF1/2/3:
-
YTH N6-methyladenosine RNA binding protein 1/2/3
- YTHDC1/2:
-
YTH domain containing 1/2
- HNRNPA2B1:
-
Heterogeneous nuclear ribonucleoprotein A2/B1
- IGFBP1:
-
Insulin like growth factor binding protein 1/2/3, /2/3
- ssGSEA:
-
Single-Sample Gene Set Enrichment Analysis
- CNV:
-
Copy number variation
- GDC:
-
Genomic Data Commons
- DEGs:
-
Differential expression genes
- TCIA:
-
The Cancer Tumor Immune Atlas
- GSVA:
-
Gene Set Variation Analysis
- UCSC:
-
University of California, Santa Cruz
- PD-1:
-
Programmed death-1
- KEGG:
-
Kyoto Encyclopedia of Gene and Genomes
- CTLA4:
-
Cytotoxic T-lymphocyte associated protein 4
- GO:
-
Gene Ontology
- CXCR4:
-
C-X-C motif chemokine receptor 4
- IPS:
-
Immunophenoscore
- SOX10:
-
SRY-box transcription factor 10
References
Cabanillas ME, McFadden DG, Durante C. Thyroid cancer. Lancet. 2016;388:2783–95. https://doi.org/10.1016/S0140-6736(16)30172-6.
Siegel RL, et al. Colorectal cancer statistics, 2020. CA Cancer J Clin. 2020;70:145–64. https://doi.org/10.3322/caac.21601.
Wang X, et al. Identification and validation of m(6)A RNA methylation regulators with clinical prognostic value in Papillary thyroid cancer. Cancer Cell Int. 2020;20:203. https://doi.org/10.1186/s12935-020-01283-y.
Desrosiers R, Friderici K, Rottman F. Identification of methylated nucleosides in messenger RNA from Novikoff hepatoma cells. Proc Natl Acad Sci USA. 1974;71:3971–5. https://doi.org/10.1073/pnas.71.10.3971.
Zaccara S, Ries RJ, Jaffrey SR. Reading, writing and erasing mRNA methylation. Nat Rev Mol Cell Biol. 2019;20:608–24. https://doi.org/10.1038/s41580-019-0168-5.
Chen XY, Zhang J, Zhu JS. The role of m(6)A RNA methylation in human cancer. Mol Cancer. 2019;18:103. https://doi.org/10.1186/s12943-019-1033-z.
Zhen D, et al. m(6)A reader: epitranscriptome target prediction and functional characterization of N (6)-methyladenosine (m(6)A) readers. Front Cell Dev Biol. 2020;8:741. https://doi.org/10.3389/fcell.2020.00741.
Nombela P, Miguel-Lopez B, Blanco S. The role of m(6)A, m(5)C and Psi RNA modifications in cancer: novel therapeutic opportunities. Mol Cancer. 2021;20:18. https://doi.org/10.1186/s12943-020-01263-w.
Jiang X, et al. The role of m6A modification in the biological functions and diseases. Signal Transduct Target Ther. 2021;6:74. https://doi.org/10.1038/s41392-020-00450-x.
Chong W, et al. m(6)A regulator-based methylation modification patterns characterized by distinct tumor microenvironment immune profiles in colon cancer. Theranostics. 2021;11:2201–17. https://doi.org/10.7150/thno.52717.
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:270–88. https://doi.org/10.1016/j.ccell.2020.02.004.
Quan C, et al. N(6)-methyladenosine in cancer immunotherapy: an undervalued therapeutic target. Front Immunol. 2021;12: 697026. https://doi.org/10.3389/fimmu.2021.697026.
Wang Y, et al. The emerging role of m6a modification in regulating the immune system and autoimmune diseases. Front Cell Dev Biol. 2021;9:755691. https://doi.org/10.3389/fcell.2021.755691.
Liu C, et al. Potential roles of N6-methyladenosine (m6A) in immune cells. J Transl Med. 2021;19:251. https://doi.org/10.1186/s12967-021-02918-y.
Quail DF, Joyce JA. The microenvironmental landscape of brain tumors. Cancer Cell. 2017;31:326–41. https://doi.org/10.1016/j.ccell.2017.02.009.
Schulz M, Salamero-Boix A, Niesel K, Alekseeva T, Sevenich L. Microenvironmental regulation of tumor progression and therapeutic response in brain metastasis. Front Immunol. 2019;10:1713. https://doi.org/10.3389/fimmu.2019.01713.
Qin S, et al. Novel immune checkpoint targets: moving beyond PD-1 and CTLA-4. Mol Cancer. 2019;18:155. https://doi.org/10.1186/s12943-019-1091-2.
Ruffo E, Wu RC, Bruno TC, Workman CJ, Vignali DAA. Lymphocyte-activation gene 3 (LAG3): the next immune checkpoint receptor. Semin Immunol. 2019;42: 101305. https://doi.org/10.1016/j.smim.2019.101305.
Yu C, et al. Combination of immunotherapy with targeted therapy: theory and practice in metastatic melanoma. Front Immunol. 2019;10:990. https://doi.org/10.3389/fimmu.2019.00990.
Pardoll DM. The blockade of immune checkpoints in cancer immunotherapy. Nat Rev Cancer. 2012;12:252–64. https://doi.org/10.1038/nrc3239.
Aghajani M, et al. Clinicopathologic and prognostic significance of programmed cell death ligand 1 expression in patients with non-medullary thyroid cancer: a systematic review and meta-analysis. Thyroid. 2018;28:349–61. https://doi.org/10.1089/thy.2017.0441.
Ulisse S, et al. PD-1 ligand expression in epithelial thyroid cancers: potential clinical implications. Int J Mol Sci. 2019. https://doi.org/10.3390/ijms20061405.
Gunda V, et al. Combinations of BRAF inhibitor and anti-PD-1/PD-L1 antibody improve survival and tumour immunity in an immunocompetent model of orthotopic murine anaplastic thyroid cancer. Br J Cancer. 2018;119:1223–32. https://doi.org/10.1038/s41416-018-0296-2.
Hanzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinform. 2013;14:7. https://doi.org/10.1186/1471-2105-14-7.
Li HB, et al. m(6)A mRNA methylation controls T cell homeostasis by targeting the IL-7/STAT5/SOCS pathways. Nature. 2017;548:338–42. https://doi.org/10.1038/nature23450.
Yao Y, et al. METTL3-dependent m(6)A modification programs T follicular helper cell differentiation. Nat Commun. 2021;12:1333. https://doi.org/10.1038/s41467-021-21594-6.
Chen DS, Mellman I. Elements of cancer immunity and the cancer-immune set point. Nature. 2017;541:321–30. https://doi.org/10.1038/nature21349.
Turley SJ, Cremasco V, Astarita JL. Immunological hallmarks of stromal cells in the tumour microenvironment. Nat Rev Immunol. 2015;15:669–82. https://doi.org/10.1038/nri3902.
Gajewski TF, et al. Cancer immunotherapy strategies based on overcoming barriers within the tumor microenvironment. Curr Opin Immunol. 2013;25:268–76. https://doi.org/10.1016/j.coi.2013.02.009.
Zhang B, et al. m(6)A regulator-mediated methylation modification patterns and tumor microenvironment infiltration characterization in gastric cancer. Mol Cancer. 2020;19:53. https://doi.org/10.1186/s12943-020-01170-0.
Acknowledgements
We are grateful to the Center for Translational Medicine of the First Affiliated Hospital of Zhengzhou University for its support in providing analytical tools. We would also like to thank the English editors of Editage (www.editage.cn).
Funding
This work was supported by a grant from Thermal Ablation of Thyroid Nodules International Joint Laboratory (Henan Province) (YUKEWAI【2016】NO.11).
Author information
Authors and Affiliations
Contributions
FHJ and ZY contributed to the conception, study design, and collected data, CS and SL analyzed the data and drafted the manuscript. And XGQ were responsible for revising and correcting the manuscript. All authors reviewed and approved the manuscript.
Corresponding author
Ethics declarations
Conflict of interest
The authors declare no competing interests.
Ethics approval
The study was conducted in accordance with the Declaration of Helsinki Principles, approval was obtained from the Ethics Committee of the First Affiliated Hospital of Zhengzhou University (ethical review number: 2021-KY-0202).
Consent to participate
Informed consent was obtained from all patients.
Consent for publication
All authors have given their consent for publication.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Springer Nature or its licensor holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Ji, FH., yang, Z., Sun, C. et al. Characterization of m6A methylation modifications and tumor microenvironment infiltration in thyroid cancer. Clin Transl Oncol 25, 269–282 (2023). https://doi.org/10.1007/s12094-022-02940-6
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s12094-022-02940-6