Analyses of m6A regulatory genes and subtype classification in atrial fibrillation

Abstract

Objective: To explore the role of m6A regulatory genes in atrial fibrillation (AF), we classified atrial fibrillation patients into subtypes by two genotyping methods associated with m6A regulatory genes and explored their clinical significance.

Methods: We downloaded datasets from the Gene Expression Omnibus (GEO) database. The m6A regulatory gene expression levels were extracted. We constructed and compared random forest (RF) and support vector machine (SVM) models. Feature genes were selected to develop a nomogram model with the superior model. We identified m6A subtypes based on significantly differentially expressed m6A regulatory genes and identified m6A gene subtypes based on m6A-related differentially expressed genes (DEGs). Comprehensive evaluation of the two m6A modification patterns was performed.

Results: The data of 107 samples from three datasets, GSE115574, GSE14975 and GSE41177, were acquired from the GEO database for training models, comprising 65 AF samples and 42 sinus rhythm (SR) samples. The data of 26 samples from dataset GSE79768 comprising 14 AF samples and 12 SR samples were acquired from the GEO database for external validation. The expression levels of 23 regulatory genes of m6A were extracted. There were correlations among the m6A readers, erasers, and writers. Five feature m6A regulatory genes, ZC3H13, YTHDF1, HNRNPA2B1, IGFBP2, and IGFBP3, were determined (p < 0.05) to establish a nomogram model that can predict the incidence of atrial fibrillation with the RF model. We identified two m6A subtypes based on the five significant m6A regulatory genes (p < 0.05). Cluster B had a lower immune infiltration of immature dendritic cells than cluster A (p < 0.05). On the basis of six m6A-related DEGs between m6A subtypes (p < 0.05), two m6A gene subtypes were identified. Both cluster A and gene cluster A scored higher than the other clusters in terms of m6A score computed by principal component analysis (PCA) algorithms (p < 0.05). The m6A subtypes and m6A gene subtypes were highly consistent.

Conclusion: The m6A regulatory genes play non-negligible roles in atrial fibrillation. A nomogram model developed by five feature m6A regulatory genes could be used to predict the incidence of atrial fibrillation. Two m6A modification patterns were identified and evaluated comprehensively, which may provide insights into the classification of atrial fibrillation patients and guide treatment.

Keywords: atrial fibrillation; bioinformatics; m6A regulatory genes; m6A subtypes; predictive model.