Novel insight into m6A regulator-mediated methylation modification patterns and immune characteristics in intracranial aneurysm

Growing evidence demonstrated that m6A modification in cardiovascular diseases. However, how it is involved in the intracranial aneurysm (IA) is still unclear. This study aimed to identify the role of m6A modification in IA.

The m6A modification may shape the IAs microenvironment and participates in the formation and rupture of IAs by regulating immune infiltration.

Three datasets downloaded from the Gene Expression Omnibus (GEO) database were used, including GSE122897, GSE15629, and GSE3679. The landscapes of 24 m6A regulators were depicted using the STRING database, Pearson's correlation analysis, and Wilcoxon test. The targets of differentially expressed m6A (DEm6A) were predicted in the m6A2Target database and the modification m6A sites of hub targets were identified in SRAMP online tool. A diagnostic model based on DEm6A was constructed and verified in training and test databases. A consensus clustering algorithm was performed to classify IA patients into distinct m6A-related clusters. Functional analyses including gene ontology, Kyoto Encyclopedia of Genes and Genomes (KEGG), gene set variation analysis, and gene set enrichment analysis analyses were conducted to elucidate the underlying mechanisms. ssGSEA algorithm was performed to uncover the immune characteristics. A PCA method was adopted to quantify the m6A score.

Nine DEm6A (IGF2BP1, IGF2BP3, YTHDF2, ZNF217, RBM15, YTHDF3, YTHDC1, FTO, and LRPPRC) significantly differed between IA and controls. Biological annotations showed that immune-related pathways (such as complement activation, inflammatory response, and interleukin signaling) and apoptosis were more enriched in IAs than in controls. Immune analyses indicate that the abundance of immune cells, immune responses, and HLA gene expression were elevated in IA samples than in controls. PCA results showed that IA has a lower m6A score than controls. An immune/apoptosis-related network modified by DEm6A was constructed. The m6A sites of six hub targets (CDK1, ASPM, AURKB, BUB1B, MKI67, and TPX2) were predicted with very high confidence. A diagnostic model with four genes (LRPPRC, YTHDF3, IGF2BP1, and ZNF217) was constructed and verified. Two m6A modification subtypes were identified with unsupervised cluster analysis. Immune infiltration analysis revealed that cluster 1 had higher immune activation than cluster 2. Further study showed that cluster 1 had a larger proportion of ruptured IAs.