Identification of m6A methylation-related genes in cerebral ischaemia‒reperfusion of Breviscapus therapy based on bioinformatics methods

Abstract

Abstract Background Cerebral ischaemia‒reperfusion (I/R) frequently causes late-onset neuronal damage. Breviscapine promotes autophagy in microvascular endothelial cells in I/R and can inhibit oxidative damage and apoptosis. However, the mediation mechanism of breviscapine on neuronal cell death is unclear. Methods First, transcriptome sequencing was performed on three groups of mice: the neuronal normal group (Control group), the oxygen-glucose deprivation/ reoxygenation group (OGD/R group) and the breviscapine administration group (Therapy group). Differentially expressed genes (DEGs) between the OGD/R and control groups and between the Therapy and OGD/R groups were obtained by the limma package. N6-methyladenosine (m6A) methylation-related DEGs were selected by Pearson correlation analysis. Then, prediction and confirmation of drug targets were performed by Swiss Target Prediction and UniProt Knowledgebase (UniProtKB) database, and key genes were obtained by Pearson correlation analysis between m6A-related DEGs and drug target genes. Next, gene set enrichment analysis (GSEA) and Ingenuity pathway analysis (IPA) were used to obtain the pathways of key genes. Finally, a circRNA-miRNA‒mRNA network was constructed based on the mRNAs, circRNAs and miRNAs. Results A total of 2250 DEGs between the OGD/R and control groups and 757 DEGs between the Therapy and OGD/R groups were selected by differential analysis. A total of 7 m6A-related DEGs, including Arl4d, Gm10653, Gm1113, Kcns3, Olfml2a, Stk26 and Tfcp2l1, were obtained by Pearson correlation analysis. Four key genes (Tfcp2l1, Kcns3, Olfml2a and Arl4d) were acquired, and GSEA showed that these key genes significantly participated in DNA repair, e2f targets and the g2m checkpoint. IPA revealed that Tfcp2l1 played a significant role in human embryonic stem cell pluripotency. The circRNA-miRNA‒mRNA network showed that mmu_circ_0001258 regulated Tfcp2l1 by mmu-miR-301b-3p. Conclusions In conclusion, four key genes, Tfcp2l1, Kcns3, Olfml2a and Arl4d, significantly associated with the treatment of OGD/R by breviscapine were identified, which provides a theoretical basis for clinical trials.