Bioinformatic analyses of the role of m6A RNA methylation regulators in abdominal aortic aneurysm.

Background: N6-methyladenosine (m6A) is found in almost all nuclear RNAs of eukaryotes, playing an important and diverse role in many biological processes. Nonetheless, the roles of m6A regulators in abdominal aortic aneurysm (AAA) unknown. Therefore, there is a pressing need to identify m6A RNA methylation regulators in the diagnosis of AAA, determination of individualized risk, discovery of therapeutic targets, and improve understanding of pathogenesis.
Methods: The GSE98278 dataset were obtained from the Gene Expression Omnibus (GEO) database to perform differential analysis of m6A-related regulators between elective stable abdominal aortic aneurysms (eAAA) and abdominal aortic aneurysm ruptured (rAAA). The random forest model was used to screen candidate m6A regulators to predict the risk of rAAA. The single sample gene set enrichment analysis (ssGSEA) method was then used to evaluate the abundance of 23 immune cells in AAA. The m6A RNA Methylation Quantification Kit was used to measure the total m6A levels of AAA and normal abdominal aorta. The terminal deoxynucleotidyl transferase mediated dUTP Nick-End Labeling (TUNEL) apoptosis assay kit was used to detect the human aortic smooth muscle cells (HASMCs) apoptotic after RBM15 knockdown. Mechanically, RBM15 knockdown was found to reduce the expression of CASP3 in an m6A-dependent manner by Western blotting, RNA immunoprecipitation (RIP), and methylated RNA immunoprecipitation-quantitative polymerase chain reaction (MeRIP-RT-PCR).
Results: RBM15, WTAP, ALKBH5, and IGFBP3 were highly expressed in rAAA. In contrast, RBM15B showed opposite results (P<0.05). The high m6A level in the rAAA compared with eAAA and normal abdominal aorta (P<0.05). The random forest model was used to screen 5 candidate m6A regulators to predict the risk of rAAA. Expression of the 5 m6A methylation regulators was validated in AAA samples (P<0.05). RBM15 knockdown inhibited the apoptosis of HASMCs. RBM15 knockdown reduced the expression of CASP3 in an m6A-dependent manner. A strong correlation between the five m6A methylation regulators and immune cell infiltration was identified. Conclusions: In summary, m6A regulators play nonnegligible roles in the occurrence of rAAA. Our investigation of m6A patterns may be able to guide future immunotherapy strategies for AAA. 2022 Annals of Translational Medicine. All rights reserved.

Introduction

Abdominal aortic aneurysm (AAA) is a permanent focal dilation of the abdominal aorta, and aorta aneurysm rupture has an incidence of approximately 80% (1-3). Ruptured abdominal aortic aneurysm (rAAA) is responsible for approximately 200,000 deaths annually (1,4). As an inflammatory and chronic vascular degenerative disease, AAA is characterized by remodeling and progressive dilation of the vessel wall (2,5-7). Lymphocytes and macrophages infiltrate the vascular wall, proteases destroy elastin and collagen in the middle and outer membrane, and there is smooth muscle cell apoptosis (8,9). Inflammatory cells activate and secrete various cytokines and chemokines to induce activation of matrix metalloproteinase, resulting in extracellular matrix degradation and smooth muscle cell apoptosis, leading to the formation and rupture of AAA (7,10-12). AAA usually develops asymptomatically and unpredictably in normal conditions until its rupture. Moreover, elective stable abdominal aortic aneurysms (eAAAs) can be fatal. In spite of extensive literature on the pathogenesis, evolution, and rupture of AAA (13-15), prevention or suppression of AAA in a clinical setting still remains an unanswered question.

It is well known that epigenetic regulation of gene expression directs and maintains cellular phenotypes (16,17). N6-methyladenosine (m6A) is the most abundant and reversible internal epigenetic modification in eukaryotic organisms (18,19). In the past few years, m6A modification of messenger RNA (mRNA) or non-coding RNA has been reported to play a vital role in almost all major normal biological processes, including the self-renewal and differentiation of embryonic stem cells, tissue development, response to heat shock or DNA damage, biological clock regulation, spermatogenesis, and transition between maternal and zygote (20,21). Using the “limma” R package (22), Li et al. noted low expression of METTL14 and HNRNPC, and reported that RBM15B was the upregulated methylation transferase used in human AAA. They also found that METTL14, HNRNPC, and RBM15B levels were correlated with immune infiltration (23). We hypothesized that m6A modification may play an important role in AAA. Nonetheless, the roles of m6A regulators in AAA unknown. Therefore, there is a pressing need to identify m6A RNA methylation regulators in the diagnosis of AAA, determination of individualized risk, discovery of therapeutic targets, and improve understanding of pathogenesis.

In this study, we firstly detected the m6A levels of total RNAs in AAA tissues and normal abdominal aortas using a m6A quantitative kit, and found high m6A levels in AAA. Furthermore, compared with eAAA, the m6A level in rAAA was upregulated. We compared the predictive performance of the support vector machine (SVM) and the random forest (RF) methods, which illustrated that RF has a substantially high accuracy in predicting the probability of AAA rupture. Moreover, we assessed the significance of differentially-expressed m6A methylation regulators in immune cell infiltration. We present the following article in accordance with the TRIPOD reporting checklist (available at https://atm.amegroups.com/article/view/10.21037/atm-22-1891/rc).

Materials

Data access

The GSE98278 dataset containing 31 eAAAs and 17 rAAAs patients were obtained from the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/). Twenty-one m6A methylation regulators were extracted from multiple literatures to identify different m6A modification patterns mediated by 21 m6A regulators. These regulators consisted of two erasers (ALKBH5 and FTO), 11 readers (YTHDC1, FMR1, ELAVL1, YTHDF1, IGF2BP1, YTHDC2, YTHDF2, HNRNPA2B1, YTHDF3, LRPPRC, and HNRNPC), and eight writers (RBM15, RBM15B, CBLL1, WTAP, METTL3, ZC3H13, METTL14, and KIAA1429).

The feature genes were screened based on the RF analysis

A flowchart for the data processing as in Figure S1. Differentially-expressed m6A methylation regulators between eAAA and rAAA were screened using the “limma” R package, and the criterion for screening differential genes was P<0.05. Pearson correlation coefficients were used to analyze differential m6A methylation regulators (r>0.3, P<0.05). The SVM and RF methods were adopted to construct a training model to predict the occurrence of rAAA. We used “Boxplots of residual”, “Reverse cumulative distribution of residual”, and a receiver operating characteristic (ROC) curve to evaluate the accuracy of the model. The RF method was used to differentially express m6A methylation regulators to predict the occurrence of rAAA, and the calculations were performed using the R library “randomForest”. In this study, mtry and ntree were set to 3 and 500, respectively. The optimal ntree was selected according to minimum cross-validation error in 10-fold cross-validation. We assessed the significance of differentially-expressed m6A methylation regulators with the optimal ntree. Nomogram construction was performed using the “rms” package. The consistency between the actual observed values and the predicted values was evaluated using calibration curves. Finally, we conducted a clinical impact curve and decision curve analysis to evaluate the clinical benefit of our model.

Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) pathway enrichment analyses

Differentially-expressed genes between eAAA and rAAA were screened using the “limma” package in R, and the criterion for screening differential genes was P<0.05. Meanwhile, KEGG pathway and GO enrichment analyses were performed using the R package, including “clusterProfiler”, “org.Hs.eg.db”, “enrichplot”, “ggplot2”, “RColorBrewer”, “dplyr”, and “ComplexHeatmap”.

m6A cluster analysis and estimation of immune cell infiltration

Using the “ConsensusClusterPlus” R package, consensus clustering was used to identify m6A clusters based on m6A methylation regulators. The single sample gene set enrichment analysis (ssGSEA) method was then used to evaluate the abundance of 23 immune cells in AAA to further investigate the correlation, and a bar plot was constructed to visualize the features.

Cell culture and transfection

Human aortic smooth muscle cells (HASMCs) from The First Hospital of China Medical University were isolated by an enzyme digestion method from the human aorta and cultured. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). The present study was approved by the Ethics Committee of The First Hospital of China Medical University (No. AF-SOP-07-1.1-01), and informed consent was taken from all the patients. HASMCs were cultured in a smooth muscle cell growth medium (SMCM), which included smooth muscle basal medium, penicillin and streptomycin (1%), Fetal bovine serum (FBS) (5%), and smooth muscle cell growth supplement (1%), at 37 ℃ in 5% carbon dioxide (CO2). HASMCs were seeded into six-well plates (1×105 cells) and then transfected with small interfering Ribonucleic Acid (siRNA). The siRNA sequence was as follows: 5'-GGUGAUAGUUGGGCAUAUA-3'.

m6A quantification

For m6A quantitation, the change in global m6A mRNA levels was measured using a m6A RNA Methylation Quantification Kit (Colorimetric, Epigentek, USA). Briefly, total RNA was extracted from AAA samples with an RNA extraction kit (TIANGEN, Beijing, China). Next, the RNA was incubated with m6A antibody in wells, and m6A levels were measured at 450 nm using a microplate reader (BioTek microplate reader, Colmar, France).

Protein extraction and Western blot analysis

A total protein extraction kit (KeyGen Biotech, Nanjing, China) was used for total protein extraction, and a bicinchoninic acid (BCA) Protein Assay Kit (KeyGen Biotech, Nanjing, China) was used for total protein quantification. Total protein was subjected to SDS-polyacrylamide gel electrophoresis (SDS-PAGE) and Western blot analysis. Primary antibodies were incubated overnight at 4 ℃. The following antibodies were used: RBM15, Caspase-3, and Actin (Cell Signaling Technology, MA, USA). Membranes were washed three times in Tris-buffered saline containing 0.05% Tween-20 (TBST) buffer and incubated with the appropriate secondary antibodies (Cell Signaling Technology, MA, USA) for 2h. The chemiluminescence method was used to detect protein bands.

Quantitative real-time polymerase chain reaction (RT-qPCR)

Total RNA was extracted using the Total RNA Extraction Kit (Solarbio, Beijing, China) according to the manufacturer’s instructions. The RNA was reverse-transcribed using the TaKaRa RNA PCR kit (TaKaRa, Dalian, China). SYBR Premix Ex Taq (TaKaRa, Dalian, China) was used for the RT-qPCR experiments. RT-qPCR was performed using an ABS 7500 Realtime PCR System (Thermo Fisher, Waltham, MA, USA) according to the manufacturer’s protocol. The qPCR primer sequences are listed in Table S1.

Terminal deoxynucleotidyl transferase mediated dUTP nick-end labeling (TUNEL) staining assay

The apoptotic cells were detected according to the TUNEL apoptosis assay kit instructions (in situ cell death detection kit-POD system, Roche, USA). HASMCs were seeded into six-well plates (1×105 cells) with 100 µm hydrogen peroxide (H2O2) treatment, and then fixed, permeabilized, and incubated with the TUNEL reagent.

RNA immunoprecipitation (RIP) assay

A RIP assay was carried out using the Magna RIP RNA-binding protein immunoprecipitation kit (Millipore, USA) according to the manufacturer’s protocol. Briefly, total RNA was extracted using the total RNA extraction kit (Solarbio, Beijing, China) according to the manufacturer’s instructions. Fragmented RNA was incubated with an antibody. The beads were washed with RIP buffer, followed by RNA purification and qRT-PCR analysis.

Methylated RNA immunoprecipitation-quantitative polymerase chain reaction (MeRIP-qPCR)

The m6A modifications of individual genes were detected using a MeRIP-qPCR assay according to the manufacturer’s instructions (Millipore, Germany). We used purified RNA incubated with m6A antibody in immunoprecipitation buffer containing RNasin Plus RNase inhibitor (Promega, Madison WI, USA). Precipitated RNA was reverse-transcribed and quantified by qPCR for m6A enrichment.

Statistical analysis

Linear regression analysis and the Pearson correlation coefficient (r) were used to determine the correlation between the gene expression levels. Non-parametric one-way analysis of variance (ANOVA) was utilized to compare variables between different groups. Two-tailed t-tests were used for all analyses comparing two groups. All statistical analysis was performed in R (Version 4.1.0). Significant significance was found when the p-values were less than 0.05.

Results

Relationship between m6A RNA methylation and rAAA

To elucidate the m6A methylation changes associated with rAAA, a differential m6A methylation regulator expression analysis was performed. The GSE98278 dataset was downloaded from the GEO database. We analyzed differentially-expressed m6A methylation regulators in the eAAA (n=31) and rAAA (n=17), and a heat map was constructed to display the expression of differential genes (). Specifically, RBM15, WTAP, ALKBH5, and IGFBP3 were highly expressed in eAAA. In contrast, RBM15B showed opposite results (). Correlation analysis also revealed a positive correlation between MTTL3 and FTO expression levels based on the GSE98278 dataset (). The m6A levels of total RNAs in AAA tissues and normal abdominal aortas were detected using a m6A quantitative kit showing the high m6A levels in AAA. Moreover, the m6A level in rAAA was upregulated compared with that in eAAA (). The changes in m6A methylation regulator expression in rAAA were validated by qRT-PCR (). Taken together, a high m6A modification methylation level was related to rAAA.

Figure 1

Relationship between m6A RNA methylation and rAAA. (A) Heat map showing the significantly different expression of five differential m6A methylation regulators (P<0.05). (B) Differential gene expression shown by box plots between eAAA and rAAA. *P<0.05, **P<0.01, ***P<0.001. (C) Analysis of the associations between the expression of METTL3 and FTO. (D) m6A quantitative analysis demonstrated the percentage of m6A content in eAAA and rAAA, *P<0.05. (E) RT-PCR analysis showed the expression of five differential m6A methylation regulators in eAAA and rAAA. *P<0.05. rAAA, ruptured abdominal aortic aneurysms; eAAA, elective stable abdominal aortic aneurysms; RT-PCR, real-time polymerase chain reaction.

The SVM and RF methods were used to construct the rAAA predictive models

We compared the prediction performance of the SVM and RF methods. “Boxplots of residual” (), “reverse cumulative distribution of residual” (), and a ROC curve () showed that RF has substantially high prediction accuracy. All of the samples in the model had minimum residuals, indicating that the model is optimal. In this study, mtry and ntree were set to 3 and 500, respectively. The optimal ntree was selected according to the minimum cross-validation error in 10-fold cross-validation (). After ranking the genes according to their importance, we identified the five m6A methylation regulators ().

Figure 2

The SVM and RF methods were used to construct rAAA predictive models. (A,B) Boxplot of the residual distribution (A) and reverse cumulative distribution of residual (B) as a function of the values of observed sensitivity between RF and SVM. (C) ROC curves showing the predictions of the two models: SVM and RF. (D) RF: prediction error curves based on 10-fold cross-validation. (E) The importance of the five m6A regulators based on the RF model. (F) Nomogram graph of the predictive model based on three m6A regulators. (G) The calibration curves showed that the nomogram model may be an ideal predictive model for rAAA. (H,I) DCA (H) and clinical impact plot (I) were used to determine the clinical utility of risk prediction nomograms. SVM, support vector machine; RF, random forest; rAAA, ruptured abdominal aortic aneurysms; ROC, receiver operating characteristic; DCA, decision curve analysis.

The nomogram evaluation model was constructed based on three m6A methylation regulators to predict the probability of AAA rupture (). The calibration curves showed that the nomogram model may be an ideal predictive model for rAAA (). Decision curve analysis (DCA) and clinical impact plots were used to determine the clinical utility of the risk prediction nomograms (). The nomogram showed remarkable predictive power for rAAA.

GO and KEGG pathway enrichment analyses

To further understand the differences between eAAA and rAAA, we screened the differential genes from the GSE98278 dataset. A total of 2,899 differential genes were screened using the “limma” package in R and displayed via a heat map ().

Figure 3

GO and KEGG pathway enrichment analyses. (A) Heat map diagram of the differential gene expression in eAAA and rAAA. (B,C) GO analysis (B) and KEGG pathway enrichment analysis (C) of differential genes. (D-H) Correlation analysis between RBM15 and CASP3, HSPA4, NFKB1, PPP3CC, and TNFRSF10B. GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; rAAA, ruptured abdominal aortic aneurysms; eAAA, elective stable abdominal aortic aneurysms.

Meanwhile, the GO and KEGG pathway enrichment analyses of differential genes were also performed. The top 12 GO results based on the P value are shown in . KEGG pathway enrichment analysis showed that the most significantly enriched signaling pathway was the lipid and atherosclerosis signaling pathway (). The previous result showed that RBM15 is the most crucial m6A methylation regulator in rAAA (). Correlation analysis was carried out to explore further the relationship between RBM15 and the lipid and atherosclerosis signaling pathway-related genes. A total of 22 genes were found to be positively correlated with RBM15 expression, and the most significant were CASP3, HSPA4, NFKB1, PPP3CC, and TNFRSF10B ().

RBM15 regulates the expression of CASP3 in an m6A-dependent manner

To further explore the role of RBM15 in AAA, the expression of RBM15 was measured in eAAA and rAAA by Western blotting. RBM15 was found to be more highly expressed in rAAA compared with eAAA (). Next, we examined the effects of RBM15 knockdown on CASP3, HSPA4, NFKB1, PPP3CC, and TNFRSF10B expression in HASMCs. The knockdown efficiency of RBM15 was assessed by Western blotting (). RT-qPCR and Western blotting showed that the mRNA level of CASP3 was decreased upon RBM15 knockdown of HASMCs ().

Figure 4

RBM15 regulates the expression of CASP3 in a m6A-dependent manner. (A,B) RBM15 protein expression was assessed by Western blotting. (C) RT-PCR analysis of CASP3, HSPA4, NFKB1, PPP3CC, and TNFRSF10B expression after RBM15 knockdown. *P<0.05. (D) Protein expressions of CASP3 were detected by Western blotting. (E) RIP assay demonstrated that RBM15 is directly bound to CASP3. *P<0.05. (F) Quantification of m6A level in total RNA in the HASMCs after RBM15 Knockdown. *P<0.05. (G) The m6A levels of CASP3 by MeRIP-qPCR in HASMCs. *P<0.05. (H) TUNEL assay showing apoptosis in HASMCs after RBM15 knockdown. *P<0.05. RT-PCR, real-time polymerase chain reaction; RIP, RNA immunoprecipitation; rAAA, ruptured abdominal aortic aneurysms; eAAA, elective stable abdominal aortic aneurysms.

Furthermore, we investigated the regulatory mechanism of CASP3 expression by RBM15. The RIP assay revealed an interaction between RBM15 and mRNA-CASP3 (). In addition, we wondered whether RBM15 could affect m6A levels, and detected m6A levels after RBM15 knockdown. RBM15 knockdown significantly decreased the m6A level (). We also measured the m6A levels of CASP3 by MeRIP-qPCR in HASMCs, which showed that m6A modification of CASP3 mRNA decreased relative to the normal control (NC) group ().

We then examined the effects of RBM15 knockdown on apoptosis by TUNEL staining. The result showed that RBM15 knockdown inhibited apoptosis (). As mentioned above, RBM15 knockdown decreased the expression of CASP3 in an m6A-dependent manner and inhibited HASMC apoptosis.

Three distinct m6A clusters were identified and immune cell infiltration

Consensus clustering was performed using the “ConsensusClusterPlus” R package, which identified three m6A clusters (m6Acluster.A, m6Acluster.B, and m6Acluster.C) based on the five m6A methylation regulators (). The heat map showing the differential gene expression in three m6A clusters is displayed in . Principal component analysis (PCA) was used to verify the three distinct m6A clusters divided by consensus clustering of the five m6A RNA methylation regulators (). Moreover, a significant difference in ALKBH5 expression was identified between the three m6A clusters ().

Figure 5

Identification of three distinct m6A clusters and immune cell infiltration. (A) Consensus clustering matrix of AAA samples for k=3. (B) The AAA were divided into three clusters by the consensus clustering method (clusters A, B and C). Heat map showing the differential gene expression in the three m6A clusters. (C) PCA was used to verify the three distinct m6A clusters. (D) ALKBH5 expression was identified in the three m6A clusters. **P<0.01. (E) Box plot showing the infiltrating immune cells in the three m6A clusters. *P<0.05. AAA, abdominal aortic aneurysms; PCA, principal component analysis.

Given that AAA has been shown to be infiltrated with a large number of immune cells, we further analyzed the infiltrating immune cells between the three m6A clusters. A significant difference was found in cluster difference 56 (CD56) bright natural killer cells and immature dendritic cells between the different m6A clusters (). Moreover, we analyzed the relationship between the expression of the five m6A methylation regulators and immune cell infiltration using the ssGSEA method (). The results were as follows: (I) ALKBH5 was positively correlated with CD56 bright natural killer cells and T follicular helper cells (); (II) RBM15 was negatively correlated with macrophages (); (III) RBM15B was positively correlated with monocytes and type 1 T helper cells (); (IV) WTAP was positively correlated with type 17 T helper cells (); and (V) IGFBP3 levels were not statistically significantly correlated with immune cell infiltration ().

Figure 6

Immune cell infiltration. (A) Heat map of the correlation between the expression of the five m6A methylation regulators and immune cells infiltration by the ssGSEA method. (B) ALKBH5 was positively correlated with CD56 bright natural killer cells and T follicular helper cells, *P<0.05. (C) RBM15 was negatively correlated with macrophages. *P<0.05. (D) RBM15B was positively correlated with monocytes and type 1 T helper cells. *P<0.05. (E) WTAP was positively correlated with type 17 T helper cells. *P<0.05. (F) IGFBP3 levels were not statistically significantly correlated with immune cell infiltration.

Discussion

AAA is a degenerative disease that may lead to aortic rupture (2,5). m6A modification is the most ubiquitous mRNA modification in eukaryotes, and the role of m6A modification has been well defined in mRNA (18,21,23). Reversible m6A modification regulates a series of biological progressions, such as RNA stability, alternative splicing, and mRNA translation (19). Abnormal expression of m6A modified protein or m6A binding protein or abnormal m6A modified level can lead to abnormal metabolism of m6A-modified related RNA and affect gene expression, which is closely related to the occurrence and development of many diseases (18,19,23,24). He et al. reported a significantly high m6A level in AAA tissues compared to normal aortic tissues, and high m6A levels are correlated with a greater risk of AAA rupture (25).

Similarly, we confirmed the high m6A level in eAAA, and the m6A level in rAAA was further upregulated compared with eAAA. In contrast to He et al., who used a multivariable model to screen risk genes, the SVM and RF Methods were applied in this study (25). RF exhibited substantially high prediction accuracy compared to “Boxplots of residual”, Reverse cumulative distribution of residual”, and the ROC curve. We compared the predictive performance of the SVM and RF. As a result, we identified the five m6A methylation regulators after ranking these genes according to their importance by RF. RBM15 is the most crucial gene.

The five m6A methylation regulators expression in rAAA were validated by RT-PCR, and RBM15 was found to be highly expressed in rAAA compared with eAAA by Western blotting. These results indicated a high m6A modification methylation level was related to rAAA. Due to the small sample size, further research is needed. Thus, we finally selected RBM15 for further analysis in this study. RBM15 is a member of the split-end (SPEN) family of proteins, which can bind to RNA by interacting with spliceosome components (26,27). RBM15/RBM15B recruit the WTAP-METTL3-METTLE14 RNA methyltransferase complex to install m6A (28). To further explore the role of RBM15 in AAA, RBM15 was depleted by siRNA-mediated knockdown. We found that RBM15 knockdown inhibited apoptosis in HASMCs.

Meanwhile, GO and KEGG pathway enrichment analyses of differential genes was performed based on the GSE98278 dataset, and the lipid and atherosclerosis signaling pathway was the most significant. Therefore, the relationship between RBM15 and the lipid and atherosclerosis signaling pathway was further investigated. A total of 22 genes were positively correlated with RBM15 expression, and the most significant were CASP3, HSPA4, NFKB1, PPP3CC, and TNFRSF10B.

Moreover, the mRNA level of CASP3 was decreased upon RBM15 knockdown of HASMCs. Mechanically, RBM15 knockdown decreased the expression of CASP3 in an m6A-dependent manner by RIP and MeRIP-qPCR. Furthermore, consensus clustering identified three m6A clusters based on the five m6A methylation regulators. Consider that AAA has been shown to be infiltrated with a large number of immune cells, we further analyzed the infiltrating immune cells between the three m6A clusters. A significant difference was observed in CD56 bright natural killer cells and immature dendritic cells between the different m6A clusters. Moreover, we determined the relationship between the expression of the five m6A methylation regulators and immune cell infiltration by the ssGSEA method.

Taken together, m6A plays a significant role in eAAA. RBM15 knockdown decreased the expression of CASP3 in an m6A-dependent manner and inhibited HASMC apoptosis. A strong correlation between the five m6A methylation regulators and immune cell infiltration was identified.

Conclusions

In this study, we observed a significantly higher m6A level in AAA tissues compared to normal aortic tissues, and high m6A levels increase the risk of AAA rupture. After ranking these genes according to their importance by RF, we identified five m6A methylation regulators, of which RBM15 was the most crucial factor. RBM15 knockdown decreased the expression of CASP3 in an m6A-dependent manner and inhibited HASMC apoptosis. Therefore, m6A may be a potential therapeutic target for inhibiting AAA rupture.

The article’s supplementary files as