Transcriptome-wide map of N6-methyladenosine (m6A) profiling in coronary artery disease (CAD) with clopidogrel resistance
Ruoyan Yu,Qinglin Yu,Zhenwei Li,Jiyi Li,Jin Yang,Yingchu Hu,Nan Zheng,Xiaojin Li,Yudie Song,Jiahui Li,Xiaomin Chen,Weiping Du,Jia Su
DOI: https://doi.org/10.1186/s13148-023-01602-w
2023-12-15
Abstract:Background: Clopidogrel resistance profoundly increases the risk of major cardiovascular events in coronary artery disease (CAD) patients. Here, we comprehensively analyse global m6A modification alterations in clopidogrel-resistant (CR) and non-CR patients. Methods: After RNA isolation, the RNA transcriptome expression (lncRNA, circRNA, and mRNA) was analysed via RNA-seq, and m6A peaks were identified by MeRIP-seq. The altered m6A methylation sites on mRNAs, lncRNAs, and circRNAs were identified, and then, GO and KEGG pathway analyses were performed. Through joint analysis with RNA-seq and MeRIP-seq data, differentially expressed mRNAs harbouring differentially methylated sites were identified. The changes in m6A regulator levels and the abundance of differentially methylated sites were measured by RT-PCR. The identification of m6A-modified RNAs was verified by m6A-IP-qPCR. Results: The expression of 2919 hypermethylated and 2519 hypomethylated mRNAs, 192 hypermethylated and 391 hypomethylated lncRNAs, and 375 hypermethylated and 546 hypomethylated circRNAs was shown to be altered in CR patients. The m6A peaks related to CR indicated lower mark density at the CDS region. Functional enrichment analysis revealed that inflammatory pathways and insulin signalling pathways might be involved in the pathological processes underlying CR. The expression of mRNAs (ST5, KDM6B, GLB1L2, and LSM14B), lncRNAs (MSTRG.13776.1 and ENST00000627981.1), and circRNAs (hsa_circ_0070675_CBC1, hsa-circRNA13011-5_CBC1, and hsa-circRNA6406-3_CBC1) was upregulated in CR patients, while the expression of mRNAs (RPS16 and CREG1), lncRNAs (MSTRG.9215.1), and circRNAs (hsa_circ_0082972_CBC1) was downregulated in CR patients. Moreover, m6A regulators (FTO, YTHDF3, and WTAP) were also differentially expressed. An additional combined analysis of gene expression and m6A peaks revealed that the expression of mRNAs (such as ST5, LYPD2, and RPS16 mRNAs) was significantly altered in the CR patients. Conclusion: The expression of m6A regulators, the RNA transcriptome, and the m6A landscape was altered in CR patients. These findings reveal epitranscriptomic regulation in CR patients, which might be novel therapeutic targets in future.