Significance of RNA N6-methyladenosine regulators in the diagnosis and subtype classification of coronary heart disease using the Gene Expression Omnibus database

Front Cardiovasc Med. 2023 Oct 19:10:1185873. doi: 10.3389/fcvm.2023.1185873. eCollection 2023.

Abstract

Background: Many investigations have revealed that alterations in m6A modification levels may be linked to coronary heart disease (CHD). However, the specific link between m6A alteration and CHD warrants further investigation.

Methods: Gene expression profiles from the Gene Expression Omnibus (GEO) databases. We began by constructing a Random Forest model followed by a Nomogram model, both aimed at enhancing our predictive capabilities on specific m6A markers. We then shifted our focus to identify distinct molecular subtypes based on the key m6A regulators and to discern differentially expressed genes between the unique m6A clusters. Following this molecular exploration, we embarked on an in-depth analysis of the biological characteristics associated with each m6A cluster, revealing profound differences between them. Finally, we delved into the identification and correlation analysis of immune cell infiltration across these clusters, emphasizing the potential interplay between m6A modification and the immune system.

Results: In this research, 37 important m6Aregulators were identified by comparing non-CHD and CHD patients from the GSE20680, GSE20681, and GSE71226 datasets. To predict the risk of CHD, seven candidate m6A regulators (CBLL1, HNRNPC, YTHDC2, YTHDF1, YTHDF2, YTHDF3, ZC3H13) were screened using the logistic regression model. Based on the seven possible m6A regulators, a nomogram model was constructed. An examination of decision curves revealed that CHD patients could benefit from the nomogram model. On the basis of the selected relevant m6A regulators, patients with CHD were separated into two m6A clusters (cluster1 and cluster2) using the consensus clustering approach. The Single Sample Gene Set Enrichment Analysis (ssGSEA) and CIBERSORT methods were used to estimate the immunological characteristics of two separate m6A Gene Clusters; the results indicated a close association between seven candidate genes and immune cell composition. The drug sensitivity of seven candidate regulators was predicted, and these seven regulators appeared in numerous diseases as pharmacological targets while displaying strong drug sensitivity.

Conclusion: m6A regulators play crucial roles in the development of CHD. Our research of m6A clusters may facilitate the development of novel molecular therapies and inform future immunotherapeutic methods for CHD.

Keywords: Gene Expression Omnibus database; consensus clustering; coronary heart disease; immune cell infiltration; m6A RNA methylation regulator; risk prediction model.

Grants and funding

This work was supported by the Key S&T Special Projects funded by the Yunnan Provincial Department of Science and Technology (No. 202102AA310003-19).