Transcriptome-wide N6-methyladenosine methylation landscape of coronary artery disease

The potential of RNA methylation in the treatment of cardiovascular diseases

RNA methylation has emerged as a dynamic regulatory mechanism that impacts gene expression and protein synthesis. Among the known RNA methylation modifications, N6-methyladenosine (m6A), 5-methylcytosine (m5C), 3-methylcytosine (m3C), and N7-methylguanosine (m7G) have been studied extensively. In particular, m6A is the most abundant RNA modification and has attracted significant attention due to its potential effect on multiple biological processes. Recent studies have demonstrated that RNA methylation plays an important role in the development and progression of cardiovascular disease (CVD). To identify key pathogenic genes of CVD and potential therapeutic targets, we reviewed several common RNA methylation and summarized the research progress of RNA methylation in diverse CVDs, intending to inspire effective treatment strategies.

N6-Methyladenosine in Vascular Aging and Related Diseases: Clinical Perspectives

Aging leads to progressive deterioration of the structure and function of arteries, which eventually contributes to the development of vascular aging-related diseases. N6-methyladenosine (m6A) is the most prevalent modification in eukaryotic RNAs. This reversible m6A RNA modification is dynamically regulated by writers, erasers, and readers, playing a critical role in various physiological and pathological conditions by affecting almost all stages of the RNA life cycle. Recent studies have highlighted the involvement of m6A in vascular aging and related diseases, shedding light on its potential clinical significance. In this paper, we comprehensively discuss the current understanding of m6A in vascular aging and its clinical implications. We discuss the molecular insights into m6A and its association with clinical realities, emphasizing its significance in unraveling the mechanisms underlying vascular aging. Furthermore, we explore the possibility of m6A and its regulators as clinical indicators for early diagnosis and prognosis prediction and investigate the therapeutic potential of m6A-associated anti-aging approaches. We also examine the challenges and future directions in this field and highlight the necessity of integrating m6A knowledge into patient-centered care. Finally, we emphasize the need for multidisciplinary collaboration to advance the field of m6A research and its clinical application.

Unveiling blood pressure-associated genes in aortic cells through integrative analysis of GWAS and RNA modification-associated variants

Background

Genome-wide association studies (GWAS) have identified more than a thousand loci for blood pressure (BP). Functional genes in these loci are cell-type specific. The aim of this study was to elucidate potentially functional genes associated with BP in the aorta through the utilization of RNA modification-associated single-nucleotide polymorphisms (RNAm-SNPs).

Methods

Utilizing large-scale genetic data of 757,601 individuals from the UK Biobank and International Consortium of Blood Pressure consortium, we identified associations between RNAm-SNPs and BP. The association between RNAm-SNPs, gene expression, and BP were examined.

Results

A total of 355 RNAm-SNPs related to m6A, m1A, m5C, m7G, and A-to-I modification were associated with BP. The related genes were enriched in the pancreatic secretion pathway and renin secretion pathway. The BP GWAS signals were significantly enriched with m6A-SNPs, highlighting the potential functional relevance of m6A in physiological processes influencing BP. Notably, m6A-SNPs in CYP11B1, PDE3B, HDAC7, ACE, SLC4A7, PDE1A, FRK, MTHFR, NPPA, CACNA1D, and HDAC9 were identified. Differential methylation and differential expression of the BP genes in FTO-overexpression and METTL14-knockdown vascular smooth muscle cells were detected. RNAm-SNPs were associated with ascending and descending aorta diameter and the genes showed differential methylation between aortic dissection (AD) cases and controls. In scRNA-seq study, we identified ARID5A, HLA-DPB1, HLA-DRA, IRF1, LINC01091, MCL1, MLF1, MLXIPL, NAA16, NADK, RERG, SRM, and USP53 as differential expression genes for AD in aortic cells.

Conclusion

The present study identified RNAm-SNPs in BP loci and elucidated the associations between the RNAm-SNPs, gene expression, and BP. The identified BP-associated genes in aortic cells were associated with AD.

Transcriptome-wide N6-methyladenosine methylation profile of atherosclerosis in mice

BACKGROUND

Atherosclerosis (AS) is a critical pathological event during the progression of cardiovascular diseases. It exhibits fibrofatty lesions on the arterial wall and lacks effective treatment. N6-methyladenosine (m6A) is the most common modification of eukaryotic RNA and plays an important role in regulating the development and progression of cardiovascular diseases. However, the role of m6A modification in AS remains largely unknown. Therefore, in this study, we explored the transcriptome distribution of m6A modification in AS and its potential mechanism.

METHODS

Methylation Quantification Kit was used to detect the global m6A levels in the aorta of AS mice. Western blot was used to analyze the protein level of methyltransferases. Methylated RNA immunoprecipitation with next-generation sequencing (MeRIP-seq) and RNA sequencing (RNA-seq) were used to obtain the first transcriptome range analysis of the m6A methylene map in the aorta of AS mice, followed by bioinformatics analysis. qRT-PCR and MeRIP-qRT-PCR were used to measure the mRNA and m6A levels in target genes.

RESULTS

The global m6A and protein levels of methyltransferase METTL3 were significantly increased in the aorta of AS mice. However, the protein level of demethylase ALKBH5 was significantly decreased. Through MeRIP-seq, we obtained m6A methylation maps in AS and control mice. In total, 26,918 m6A peaks associated with 13,744 genes were detected in AS group, whereas 26,157 m6A peaks associated with 13,283 genes were detected in the control group. Peaks mainly appeared in the coding sequence (CDS) regions close to the stop codon with the RRACH motif. Moreover, functional enrichment analysis demonstrated that m6A-containing genes were significantly enriched in AS-relevant pathways. Interestingly, a negative correlation between m6A methylation abundance and gene expression level was found through the integrated analysis of MeRIP-seq and RNA-seq data. Among the m6A-modified genes, a hypo-methylated but up-regulated (hypo-up) gene Fabp5 may be a potential biomarker of AS.

CONCLUSIONS

Our study provides transcriptome-wide m6A methylation for the first time to determine the association between m6A modification and AS progression. Our study lays a foundation for further exploring the pathogenesis of AS and provides a new direction for the treatment of AS.

Dynamic alterations of the transcriptome-wide m6A methylome in cytogenetically normal acute myeloid leukaemia during initial diagnosis and relapse

Accumulating studies have indicated that N6-methyladenosine (m6A) plays an important role in acute myeloid leukaemia (AML). However, little is known about the m6A methylome at a transcriptome-wide scale in AML patients. We obtained three pairs of bone marrow (BM) samples from cytogenetically normal AML patients at the timepoints of diagnosis (AML) and relapse (R_AML) and three BM samples from healthy donors used as normal controls (NCs). Methylated RNA immunoprecipitation next-generation sequencing (MeRIP-Seq) was conducted to identify differences in the m6A methylomes between AML and NC and between R_AML and AML. We identified a total of 11,076 and 11,962 differential m6A peaks in AML and R_AML group, respectively. These dysregulated m6A peaks were detected on all chromosomes, especially chr1, chr19 and chr17, and were mainly enriched in 3' untranslated regions, stop codon and coding sequence regions. Moreover, GO and KEGG analyses indicated that m6A -modified genes were significantly enriched in cancer-related biological functions and pathways. Additionally, we identified a link between the m6A methylome and RNA transcriptome via combined analyses of MeRIP-seq and RNA-seq data. In addition, 5 genes, HSPG2, HOMER3, TSPO2, CXCL12 and FUT1 regulated by m6A modification potentially, were shown to be related to the prognosis of AML patients. Additionally, we detected the mRNA expression of major m6A regulators and potential target mRNA on the leukemogenesis and found that the expression of IGF2BP2, HSPG2 and HOMER3 were upregulated in AML at the time of diagnosis. Moreover, their expression became downregulated after remission and then elevated again at relapse. Our study provides the first data on the differential m6A methylome in AML patients during initial diagnosis and relapse. This study demonstrates a novel relationship between m6A modification and AML relapse and paves the way for further studies aimed at elucidating the epigenic mechanisms involved in the relapse of AML.

N6-methyladenosine in macrophage function: a novel target for metabolic diseases

N6-methyladenosine (m⁶A) is one of the most prevalent internal transcriptional modifications. Evidence has highlighted changes in m⁶A in metabolic disorders and various metabolic diseases. However, the precise mechanisms of these m⁶A changes in such conditions are not understood. Macrophages are crucial for the innate immune system and exert either beneficial or harmful roles in metabolic disease. Notably, m⁶A was found to be closely related to macrophage phenotype and dysfunction. In this review, we summarize m⁶A in macrophage function from the perspective of macrophage development, activation, and polarization, pyroptosis, and metabolic disorders. Furthermore, we discuss how m⁶A-mediated macrophage function affects metabolic diseases, including atherosclerosis and nonalcoholic fatty liver disease (NAFLD). Finally, we discuss challenges and prospects for m⁶A in macrophage and metabolic diseases with the aim of providing guidance for the treatment of metabolic diseases.

The Mechanism and Role of N6-Methyladenosine (m6A) Modification in Atherosclerosis and Atherosclerotic Diseases

N6-methyladenosine (m6A) modification is a newly discovered regulatory mechanism in eukaryotes. As one of the most common epigenetic mechanisms, m6A's role in the development of atherosclerosis (AS) and atherosclerotic diseases (AD) has also received increasing attention. Herein, we elucidate the effect of m6A on major risk factors for AS, including lipid metabolism disorders, hypertension, and hyperglycemia. We also describe how m6A methylation contributes to endothelial cell injury, macrophage response, inflammation, and smooth muscle cell response in AS and AD. Subsequently, we illustrate the m6A-mediated aberrant biological role in the pathogenesis of AS and AD, and analyze the levels of m6A methylation in peripheral blood or local tissues of AS and AD, which helps to further discuss the diagnostic and therapeutic potential of m6A regulation for AS and AD. In summary, studies on m6A methylation provide new insights into the pathophysiologic mechanisms of AS and AD, and m6A methylation could be a novel diagnostic biomarker and therapeutic target for AS and AD.

Exploration of N6-Methyladenosine Profiles of mRNAs and the Function of METTL3 in Atherosclerosis

Objectives: N6-methylladenosine (m6A) modification has not been fully studied in atherosclerosis. The objectives of this study were to investigate differentially expressed m6A methylated peaks and mRNAs, along with the regulatory role of methyltransferase 3 (METTL3) in pathological processes of atherosclerosis. Methods: The pathological models of human coronary artery smooth muscle cells (HCASMCs) were induced in vitro. The differentially expressed mRNAs and m6A peaks were identified by RNA-Seq and meRIP-Seq. The potential mechanisms were analyzed via bioinformatic assays. Methylases expression was tested by quantitative real-time polymerase chain reaction (qRT-PCR) and Western blotting (WB) in HCASMCs, and by immunohistochemical assays in 40 human coronary arteries. The knockdown of METTL3 expression in cells was performed by siRNA transfection, and cell proliferation and migration were detected after transfection. Results: We identified 5121 m6A peaks and 883 mRNAs that were expressed differentially in the pathological processes of HCASMCs. Bioinformatic analyses showed that the different m6A peaks were associated with cell growth and cell adhesion, and the 883 genes showed that the extracellular matrix and PI3K/AKT pathway regulate the processes of HCASMCs. Additionally, 10 hub genes and 351 mRNAs with differential methylation and expression levels were found. METTL3 was upregulated in the arteries with atherosclerotic lesions and in the proliferation and migration model of HCASMCs, and pathological processes of HCASMCs could be inhibited by the knockdown of METTL3. The mechanisms behind regulation of migration and proliferation reduced by siMETTL3 are concerned with protein synthesis and energy metabolism. Conclusions: These results revealed a new m6A epigenetic method to regulate the progress of atherosclerosis, which suggest approaches for potential therapeutic interventions that target METTL3 for the prevention and treatment of coronary artery diseases.

Coronary artery disease risk factors affected by RNA modification-related genetic variants

Background

Single nucleotide polymorphisms that affect RNA modification (RNAm-SNPs) may have functional roles in coronary artery disease (CAD). The aim of this study was to identify RNAm-SNPs in CAD susceptibility loci and highlight potential risk factors.

Methods

CAD-associated RNAm-SNPs were identified in the CARDIoGRAMplusC4D and UK Biobank genome-wide association studies. Gene expression and circulating protein levels affected by the RNAm-SNPs were identified by QTL analyses. Cell experiments and Mendelian randomization (MR) methods were applied to test whether the gene expression levels were associated with CAD.

Results

We identified 81 RNAm-SNPs that were associated with CAD or acute myocardial infarction (AMI), including m6A-, m1A-, m5C-, A-to-I- and m7G-related SNPs. The m6A-SNPs rs3739998 in JCAD, rs148172130 in RPL14 and rs12190287 in TCF21 and the m7G-SNP rs186643756 in PVT1 were genome-wide significant. The RNAm-SNPs were associated with gene expression (e.g., MRAS, DHX36, TCF21, JCAD and SH2B3), and the expression levels were associated with CAD. Differential m6A methylation and differential expression in FTO-overexpressing human aorta smooth muscle cells and peripheral blood mononuclear cells of CAD patients and controls were detected. The RNAm-SNPs were associated with circulating levels of proteins with specific biological functions, such as blood coagulation, and the proteins (e.g., cardiotrophin-1) were confirmed to be associated with CAD and AMI in MR analyses.

Conclusion

The present study identified RNAm-SNPs in CAD susceptibility genes, gene expression and circulating proteins as risk factors for CAD and suggested that RNA modification may play a role in the pathogenesis of CAD.

m6A regulator-mediated RNA methylation modification patterns are involved in immune microenvironment regulation of coronary heart disease

Background

Although the roles of m6A modification in the immune responses to human diseases have been increasingly revealed, their roles in immune microenvironment regulation in coronary heart disease (CHD) are poorly understood.

Methods

The GSE20680 and GSE20681 datasets related to CHD were acquired from the Gene Expression Omnibus (GEO) database. A total of 30 m6A regulators were used to perform LASSO regression to identify the significant genes involved in CHD. Unsupervised clustering analysis was conducted using the m6A regulators to distinguish the m6A RNA methylation patterns in patients with CHD. The differentially expressed genes (DEGs) and biological characteristics, including GO and KEGG enrichment results, were assessed for the different m6A patterns to analyse the impacts of m6A regulators on CHD. Hub genes were identified, and subsequent microRNAs-mRNAs (miRNAs–mRNAs) and mRNAs-transcriptional factors (mRNA-TFs) interaction networks were constructed by the protein and protein interaction (PPI) network method using Cytoscape software. The infiltrating proportion of immune cells was assessed by ssGSEA and the CIBERSORT algorithm. Quantitative real-time PCR (qRT-PCR) was performed to detect the expression of the significant m6A regulators and hub genes.

Results

Four of 30 m6A regulators (HNRNPC, YTHDC2, YTHDF3, and ZC3H13) were identified to be significant in the development of CHD. Two m6A RNA methylation clusters were distinguished by unsupervised clustering analysis based on the expression of the 30 m6A regulators. A total of 491 genes were identified as DEGs between the two clusters. A PPI network including 308 mRNAs corresponding to proteins was constructed, and 30 genes were identified as hub genes that were enriched in the bioprocesses of peptide cross-linking, keratinocyte differentiation. Twenty-seven hub genes were found to be related to miRNAs, and seven hub genes were found to be related to TFs. Moreover, among the 30 hub genes, eight genes were found to be upregulated in CHD, and three were found to be downregulated in CHD compared to the normal people. The high m6A modification pattern was associated with a higher infiltrated abundance of immune cells.

Conclusion

Our findings demonstrated that m6A modification plays crucial roles in the diversity and complexity of the immune microenvironment in CHD.

Validated Impacts of N6-Methyladenosine Methylated mRNAs on Apoptosis and Angiogenesis in Myocardial Infarction Based on MeRIP-Seq Analysis

Objectives: N6-methyladenosine (m⁶A) is hypothesized to play a role in the regulation of pathogenesis of myocardial infarction (MI). This study was designed to compare m⁶A-tagged transcript profiles to identify mRNA-specific changes on pathophysiological variations after MI.

Methods: N6-methyladenosine methylated RNA immunoprecipitation sequencing (MeRIP-seq) and RNA sequencing (RNA-seq) were interacted to select m⁶A-modified mRNAs with samples collected from sham operated and MI rat models. m⁶A methylation regulated mRNAs were interacted with apoptosis/angiogenesis related genes in GeneCards. Afterwards, MeRIP-quantitative real-time PCR (MeRIP-qRT-PCR) was performed to measure m⁶A methylation level of hub mRNAs. m⁶A methylation variation was tested under different oxygen concentration or hypoxic duration in H9c2 cells and HUVECs. In addition, Western blot and qRT-PCR were employed to detect expression of hub mRNAs and relevant protein level. Flow cytometry and Tunel assay were conducted to assess apoptotic level. CCK-8, EdU, and tube formation assay were performed to measure cell proliferation and tube formation ability.

Results: Upregulation of Mettl3 was firstly observed in vivo and in vitro, followed by upregulation of m⁶A methylation level. A total of 567 significantly changed m⁶A methylation peaks were identified, including 276 upregulated and 291 downregulated peaks. A total of 576 mRNAs were upregulated and 78 were downregulated. According to combined analysis of MeRIP-seq and RNA-seq, we identified 26 significantly hypermethylated and downregulated mRNAs. Based on qRT-PCR and interactive analysis, Hadh, Kcnn1, and Tet1 were preliminarily identified as hub mRNAs associated with apoptosis/angiogenesis. MeRIP-qRT-PCR assay confirmed the results from MeRIP-seq. With the inhibition of Mettl3 in H9c2 cells and HUVECs, downregulated m⁶A methylation level of total RNA and upregulated expression of hub mRNAs were observed. Increased m⁶A level was verified in the gradient context in terms of prolonged hypoxic duration and decreased oxygen concentration. Under simulated hypoxia, roles of Kcnn1 and Tet1 in angiogenesis and Hadh, Tet1, and Kcnn1 in apoptosis were further confirmed with our validation experiments.

Conclusion: Roles of m⁶A-modified mRNA transcripts in the context of MI were preliminarily verified. In the context of m⁶A methylation, three hub mRNAs were validated to impact the process of apoptosis/angiogenesis. Our study provided theoretical basis and innovative targets for treatment of MI and paved the way for future investigations aiming at exploring upstream epigenetic mechanisms of pathogenesis after MI.

Epitranscriptomics of Ischemic Heart Disease-The IHD-EPITRAN Study Design and Objectives

Epitranscriptomic modifications in RNA can dramatically alter the way our genetic code is deciphered. Cells utilize these modifications not only to maintain physiological processes, but also to respond to extracellular cues and various stressors. Most often, adenosine residues in RNA are targeted, and result in modifications including methylation and deamination. Such modified residues as N-6-methyl-adenosine (m6A) and inosine, respectively, have been associated with cardiovascular diseases, and contribute to disease pathologies. The Ischemic Heart Disease Epitranscriptomics and Biomarkers (IHD-EPITRAN) study aims to provide a more comprehensive understanding to their nature and role in cardiovascular pathology. The study hypothesis is that pathological features of IHD are mirrored in the blood epitranscriptome. The IHD-EPITRAN study focuses on m6A and A-to-I modifications of RNA. Patients are recruited from four cohorts: (I) patients with IHD and myocardial infarction undergoing urgent revascularization; (II) patients with stable IHD undergoing coronary artery bypass grafting; (III) controls without coronary obstructions undergoing valve replacement due to aortic stenosis and (IV) controls with healthy coronaries verified by computed tomography. The abundance and distribution of m6A and A-to-I modifications in blood RNA are charted by quantitative and qualitative methods. Selected other modified nucleosides as well as IHD candidate protein and metabolic biomarkers are measured for reference. The results of the IHD-EPITRAN study can be expected to enable identification of epitranscriptomic IHD biomarker candidates and potential drug targets.