Biomarkers related to m6A and succinic acid metabolism in papillary thyroid carcinoma

Abstract

Abstract Background Studies have shown that m6A modification is related to the occurrence and development of papillary thyroid carcinoma (PTC). The disorder of succinic acid metabolism is associated with the occurrence and development of various tumors. However, there are few studies based on m6A and succinate metabolism-related genes (SMRGs) in PTC. Methods The TCGA-Thyroid carcinoma (THCA), GSE33630, 1159 SMRGs, and 23 m6A regulatory factors were collected from the online databases. Subsequently, the differentially expressed genes (DEGs) were selected between PTC (Tumor) and Normal samples. The overlapping genes among the DEGs, m6A, and SMRGs were applied to screen the biomarkers. Using the 3 machine-learning algorithms, the biomarkers were determined based on the overlapping genes. Next, the biomarkers were evaluated by the ROC curve and expression analysis in TCGA-THCA and GSE33630. Then, the overall survival (OS) differences were compared between the high-and low-expression biomarkers. Finally, immune infiltration analysis, molecular regulatory network, and drug prediction were performed based on the biomarkers. Results In TCGA-THCA, there were 2800 DEGs between and Normal samples, and then 7 overlapping genes were obtained. Importantly, ADK, TNFRSF10B, CYP7B1, FGFR2, and CPQ were determined as biomarkers with excellent diagnostic efficiency (AUC > 0.7). In PTC samples, ADK and TNFRSF10B were high-expressed while CYP7B1, FGFR2, and CPQ were low-expressed. Especially, the high-expression groups of ADK had a better prognosis, while the high-expression groups of CYP7B1, FGFR2, and CPQ had a worse prognosis. Afterward, immune infiltration analysis found that 16 immune cells had infiltration differences between the Tumor and Normal samples. Finally, transcription factor SP1 could regulate CYP7B1 and TNFRSF10B. Moreover, Navitoclax was a potential drug for PTC patients. Conclusion Overall, we described 5 biomarkers associated with adverse prognosis of PTC, including ADK, TNFRSF10B, CYP7B1, FGFR2, and CPQ. All these biomarkers were involved in succinate metabolism and m6A modification of RNA. This set of biomarkers should be explored further for their diagnostic value in PTC. Investigations into the mechanistic role of alteration of succinate metabolism and m6A modification of RNA pathways in the pathophysiology of PTC are warranted.

Keywords

• m6A
• Papillary thyroid carcinoma
• Succinic acid metabolism
• Bioinformatics