Li Minyu, Fu Xiaodan, Zhou Tianhan, Han Hui
Zhejiang Chinese Medical University, Hangzhou, Zhejiang, China.
Department of Endocrinology, Affiliated Hangzhou First People's Hospital, School of medicine, Westlake University, Hangzhou, Zhejiang, China.
BMC Med Genomics. 2024 Aug 7;17(1):199. doi: 10.1186/s12920-024-01975-8.
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.
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.
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.
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.
研究表明,m6A修饰与甲状腺乳头状癌(PTC)的发生发展有关。琥珀酸代谢紊乱与多种肿瘤的发生发展相关。然而,基于m6A和琥珀酸代谢相关基因(SMRGs)在PTC中的研究较少。
从在线数据库收集TCGA-甲状腺癌(THCA)、GSE33630、1159个SMRGs和23个m6A调控因子。随后,在PTC(肿瘤)和正常样本之间筛选差异表达基因(DEGs)。将DEGs、m6A和SMRGs中的重叠基因用于筛选生物标志物。使用三种机器学习算法,基于重叠基因确定生物标志物。接下来,通过ROC曲线以及在TCGA-THCA和GSE33630中的表达分析对生物标志物进行评估。然后,比较高表达和低表达生物标志物之间的总生存期(OS)差异。最后,基于生物标志物进行免疫浸润分析、分子调控网络分析和药物预测。
在TCGA-THCA中,肿瘤样本和正常样本之间有2800个DEGs,随后获得了7个重叠基因。重要的是,已确定ADK、TNFRSF10B、CYP7B1、FGFR2和CPQ为诊断效率优异的生物标志物(AUC>0.7)。在PTC样本中,ADK和TNFRSF10B高表达,而CYP7B1、FGFR2和CPQ低表达。特别是,ADK高表达组预后较好,而CYP7B1、FGFR2和CPQ高表达组预后较差。随后,免疫浸润分析发现肿瘤样本和正常样本之间有16种免疫细胞存在浸润差异。最后,转录因子SP1可调控CYP7B1和TNFRSF10B。此外,Navitoclax是PTC患者的一种潜在药物。
总体而言,我们描述了5种与PTC不良预后相关的生物标志物,包括ADK、TNFRSF10B、CYP7B1、FGFR2和CPQ。所有这些生物标志物都参与琥珀酸代谢和RNA的m6A修饰。这组生物标志物在PTC中的诊断价值应进一步探索。有必要研究琥珀酸代谢改变和RNA途径的m6A修饰在PTC病理生理学中的作用机制。