Identification and verification of m7G-Related genes as biomarkers for prognosis of sarcoma

Front Genet. 2023 Feb 3:14:1101683. doi: 10.3389/fgene.2023.1101683. eCollection 2023.

Abstract

Background: Increasing evidence indicates a crucial role for N7-methylguanosine (m7G) methylation modification in human disease development, particularly cancer, and aberrant m7G levels are closely associated with tumorigenesis and progression via regulation of the expression of multiple oncogenes and tumor suppressor genes. However, the role of m7G in sarcomas (SARC) has not been adequately evaluated. Materials and methods: Transcriptome and clinical data were gathered from the TCGA database for this study. Normal and SARC groups were compared for the expression of m7G-related genes (m7GRGs). The expression of m7GRGs was verified using real-time quantitative PCR (RT-qPCR) in SARC cell lines. Then, differentially expressed genes (DEGs) were identified between high and low m7GRGs expression groups in SARC samples, and GO enrichment and KEGG pathways were evaluated. Next, prognostic values of m7GRGs were evaluated by Cox regression analysis. Subsequently, a prognostic model was constructed using m7GRGs with good prognostic values by Lasso regression analysis. Besides, the relationships between prognostic m7GRGs and immune infiltration, clinical features, cuproptosis-related genes, and antitumor drugs were investigated in patients with SARC. Finally, a ceRNA regulatory network based on m7GRGs was constructed. Results: The expression of ten m7GRGs was higher in the SARC group than in the control group. DEGs across groups with high and low m7GRGs expression were enriched for adhesion sites and cGMP-PKG. Besides, we constructed a prognostic model that consists of EIF4A1, EIF4G3, NCBP1, and WDR4 m7GRGs for predicting the survival likelihood of sarcoma patients. And the elevated expression of these four prognostic m7GRGs was substantially associated with poor prognosis and elevated expression in SARC cell lines. Moreover, we discovered that these four m7GRGs expressions were negatively correlated with CD4+ T cell levels, dendritic cell level and tumor purity, and positively correlated with tumor mutational burden, microsatellite instability, drug sensitivity and cuproptosis-related genes in patients with sarcomas. Then, a triple regulatory network of mRNA, miRNA, and lncRNA was established. Conclusion: The current study identified EIF4A1, EIF4G3, NCBP1, and WDR4 as prognostic genes for SARC that are associated with m7G.These findings extend our knowledge of m7G methylation in SARC and may guide the development of innovative treatment options.

Keywords: ceRNA regulatory network; m7G; prognostic model; sarcoma; tumor immune microenvironment.

Grants and funding

This research was continuously funded by National Natural Science Foundation of China (No. 82102568; 82172432 and 82001319), National & Local Joint Engineering Research Center of Orthopaedic Biomaterials (XMHT20190204007), Shenzhen High-level Hospital Construction Fund, Shenzhen Key Medical Discipline Construction Fund (No. SZXK023), Shenzhen “San-Ming” Project of Medicine (No. SZSM201612092), Research and Development Projects of Shenzhen (No. Z2021N054), Shenzhen Science and Technology Program (No. JCYJ20210324110214040 and JCYJ20190809152409606), Guangdong Basic and Applied Basic Research Foundation (No. 2021A1515012586), Bethune Charitable Foundation and CSPC Osteoporosis Research Foundation Project (No. G-X-2020–1107-21), and The Scientific Research Foundation of Peking University Shenzhen Hospital (No. KYQD2021099).