Skip to main content

MiR-146a rs2910164 (G/C) polymorphism is associated with the development and prognosis of acute coronary syndromes: an observational study including case control and validation cohort

Abstract

Background

Polymorphisms in microRNAs (miRNAs) play an important role in acute coronary syndromes (ACS). The purpose of this study was to assess the association of miR-146a rs2910164 and miR-34b rs4938723 polymorphisms with the development and prognosis of ACS and to explore the underlying mechanisms.

Methods

A case–control study of 1171 subjects was included to determine the association of miR-146a rs2910164 and miR-34b rs4938723 polymorphisms with ACS risk. An additional 612 patients with different miR-146a rs2910164 genotypes, who underwent percutaneous coronary intervention (PCI) were included in the validation cohort and followed for 14 to 60 months. The endpoint was major adverse cardiovascular events (MACE). A luciferase reporter gene assay was used to validate the interaction of oxi-miR-146a(G) with the IKBA 3'UTR. Potential mechanisms were validated using immunoblotting and immunostaining.

Results

The miR-146a rs2910164 polymorphism was significantly associated with the risk of ACS (Dominant model: CG + GG vs. CC, OR = 1.270, 95% CI (1.000–1.613), P = 0.049; Recessive model: GG vs. CC + CG, OR = 1.402, 95% CI (1.017–1.934), P = 0.039). Serum inflammatory factor levels were higher in patients with the miR-146a rs2910164 G allele than in those with the C allele. MiR-146a rs2910164 polymorphism in dominant model was associated with the incidence of MACE in post-PCI patients (CG + GG vs. CC, HR = 1.405, 95% CI (1.018–1.939), P = 0.038). However, the miR-34b rs4938723 polymorphism was not associated with the prevalence and prognosis of ACS. The G allele of miR-146a rs2910164 tends to be oxidized in ACS patients. The miRNA fractions purified from monocytes isolated from ACS patients were recognized by the 8OHG antibody. Mispairing of Oxi-miR-146a(G) with the 3'UTR of IKBA results in decreased IκBα protein expression and activation of the NF-κB inflammatory pathway. P65 expression was higher in atherosclerotic plaques from patients carrying the miR-146a rs2910164 G allele.

Conclusion

The variant of miR-146a rs2910164 is closely associated with the risk of ACS in Chinese Han population. Patients carrying miR-146a rs2910164 G allele may have worse pathological change and poorer post-PCI prognosis, partly due to the oxidatively modified miR-146a mispairing with 3′UTR of IKBA and activating NF-κB inflammatory pathways.

Introduction

Despite significant advances in the diagnosis and treatment of coronary artery disease (CAD), cardiovascular disease remains the leading cause of death worldwide, with nearly half of these deaths due to ischemic heart disease [1, 2]. As the most important type of ischemic heart diseases, acute coronary syndromes (ACS), characterized pathologically by unstable atherosclerotic lesion, commonly begin with atherosclerotic plaque rupture and intracoronary thrombus formation and result in cardiovascular disease events [3].

MicroRNAs (miRNAs) are approximately 22 nucleotide RNAs that negatively regulate gene expression by inhibiting mRNA translation or promoting mRNA degradation [4] thereby participate in physiological processes such as cell differentiation and proliferation, apoptosis, and organ development [5, 6]. It has been demonstrated that miRNAs are closely related to the pathogenesis of cardiovascular disease [7]. MiR-146a and miR-34b are significantly up-regulated in human atherosclerotic plaques [8]. In addition, circulating miR-146a levels were significantly increased in ACS patients [9]. Furthermore, the inhibition of the miR-34 family reduces ventricular remodeling and cardiac dysfunction due to myocardial infarction (MI) or pressure overload, and increases angiogenesis [10].

The polymorphism of miRNAs may alter their target selection and change their original biological function, leading to gene susceptibility as one of the risk factors in many diseases [11]. Polymorphisms of miR-146a rs2910164 and miR-34b rs4938723 have been reported to be associated with the development and prognosis of Crohn’s disease and cancer [12, 13]. Recent studies have further demonstrated that miR-146a rs2910164 is associated with the pathogenesis of ischemic stroke [14]. However, no studies have elaborated the roles of miR-146a rs2910164 and miR-34b rs4938723 in the development of ACS and their underlying mechanisms.

Our present work demonstrated that the miR-146a rs2910164, not miR-34b rs4938723, polymorphism is a risk factor for ACS and that patients carrying the G allele have higher mortality after percutaneous coronary intervention (PCI). Our data suggest a mechanism by which oxidatively modified miR-146a rs2910164 mismatches with the 3ʹUTR of IKBA, thereby activating NF-κB inflammatory pathways that leads to cardiovascular events.

Methods

Study population

This research consists of two independent studies: a case–control and a validation cohort. In the case–control study, we screened a total of 1605 CAD patients and finally included 1171 Han patients aged 24–80 years. Patients with stable angina (SA) were included in the control group (n = 552). The ACS group (n = 619) included patients with acute ST-segment elevation myocardial infarction (STEMI), acute non-ST-segment elevation myocardial infarction (NSTEMI), and unstable angina (UA).

The prospective cohort study screened a total of 697 patients with ACS undergoing primary percutaneous coronary intervention (PCI), 53 patients were excluded according to the exclusion criteria, 32 patients (4.5%) were lost to follow-up, and finally 612 patients were included in the study. Based on the genetic analysis, they were divided into miR-146a rs2910164 CC, CG and GG groups. Subjects with different genotypes were followed for 14 to 60 months (median follow-up time is 42 months) to assess the incidence of mortality and MACE. The endpoint was major adverse cardiovascular events (MACE), including all-cause death, nonfatal acute myocardial infarction, revascularization, and stroke.

This study was conducted at the First Affiliated Hospital of Xi'an Jiaotong University from January 2016 to October 2021. Patients were excluded from the case–control study and the cohort study with the following cases: acute infections, post-revascularization, serious liver or kidney disease, cancer or acute stroke. Written informed consent was obtained from all study participants, and this study was approved by the ethnic committee approval of the First Affiliated Hospital of Xi’an Jiaotong University.

Clinical measurements

Well-trained interviewers and investigators gathered all samples and clinical data. A history of smoking was defined as either past or current history of cigarettes. Hypertension was defined as either the systolic blood pressure higher than 140 mmHg or diastolic blood pressure higher than 90 mmHg. Diabetes mellitus (DM) was defined as fasting plasma glucose ≥ 7.0 mmol/L or use of anti-diabetic therapies, or HbA1c ≥ 6.5%. The body weight (kg) and height (m) of each subject were acquired during the initial visit. Body mass index (BMI) was calculated by the equation of body weight (kg)/height square (m2).

Genetic analysis

Genomic DNA was extracted from white blood cells using the commercially available DNA isolation kit (Tiangen Biotech, Beijing, China) according to the manufacturer’s instructions. Genotypes for individual DNA samples were genotyped using the ABI PRISM- Snapshot method (Applied Bio system, CA, USA). More detail information was shown in Additional files. Genotypes of 2 SNPs were identified by capillary electrophoresis (ABI PRISM3730 DNA Sequencer; Applied Biosystems). The results were analyzed with GeneMapper 3.0 software (Applied Biosystems). All SNapShot and PCR primers were listed in Additional file 1: Table S2. Genotyping quality control procedures for study have been described in detail [15]. In order to evaluate the multiplex SNaPShot results, 40 samples were randomly selected and regenotyped by direct sequencing using a BigDye terminator (Applied Biosystem). All assays were 100% concordant.

Measurement of serum cytokines

Overnight after admission, all blood samples (5 mL) were collected into tubes containing anticoagulant and then were centrifuged for 10 min at 3000 rpm. The concentrations of serum high sensitivity C-reactive protein (hs-CRP), fasting glucose, high density lipoprotein cholesterol(HDL-C) and triglyceride(TG) were determined at the laboratory in the hospital. Low-density lipoprotein cholesterol (LDL-C) was calculated using the Friedewald formula. Interleukin-1β (IL-1β), interleukin-6 (IL-6), and tumor necrosis factor-α (TNF-α) in serum were determined by ELISA kits (R&D Systems, Minneapolis, MN, USA).

Calculation of Gensini’s score

The Gensini score was calculated for each patient from the coronary angiogram. Gensini score grades a severity score to each coronary stenosis as 1 for 1 to 25% narrowing, 2 for 26 to 50%, 4 for 51 to 75%, 8 for 76% to 90%, 16 for 91% to 99%, and 32 for a completely occluded artery. The score is then multiplied by a factor according to the importance of the coronary artery. The multiplication factor is 5 for a left main coronary artery, 2.5 for proximal left anterior descending artery and proximal circumflex artery, 1.5 for a mid left anterior descending artery, and 1 for distal left anterior descending artery, mid or distal circumflex artery, and right coronary artery. The multiplication factor for any other branch is 0.5. Gensini score was showed as the sum of the scores for all coronary arteries.

Statistical analysis

Categorical variables are presented (%) and continuous variables are presented as mean ± SD. Data from all groups were compared using theχ2 test for categorical variables and a Student t test for continuous variables. The odds ratio and 95% confidence intervals (CI) were calculated using logistic regression to achieve the relative risks of various genotypes for ACS. The association between different genotypes and the risks of MACE was estimated by Kaplan–Meier survival and Cox regression analysis. Hazard ratios (HR) were presented with 95% CI to reveal the risk of an event when this factor is appeared. All calculations and graphs were performed using R version 4.2.2 (Vienna, Austria) and Graphpad prism version 8.0 softwares. Statistically significant differences were defined as P < 0.05 using a two-tailed test.

Results

MiR-146a and miR-34b, two microRNAs with rs2910164 and rs4938723 polymorphisms, are highly expressed in CAD patients

To determine which microRNAs are highly expressed in patients with coronary artery disease (CAD), we downloaded the non-coding RNA dataset GSE64566 from the NCBI_GEO database (https://www.ncbi.nlm.nih.gov/gds) [16]. Epicardial adipose tissues (EAT) from control (n = 8) and CAD (n = 11) patients was used as samples for differential gene expression analysis. Principal component analysis (PCA) revealed significant differences in gene expression between the control and CAD groups (Fig. 1A). Volcano plots and heatmaps showed increased expression of 13 microRNAs in CAD patients (Fig. 1B, C). See Additional file 2: Excel1_Different miRNA for details. We then performed a secondary analysis of the data on microRNA expression in atherosclerotic plaques provided by Emma Raitoharju’s article[8] and found that seven microRNAs were highly expressed in atherosclerotic plaques (Additional file 1: Table S1). By comparison, we found that miR-146a and miR-34b were increased in both epicardial fat and atherosclerotic plaques in CAD patients (Fig. 1C–E, Additional file 1: Table S1). Next, All SNPs of miR-146a and miR-34b were then screened using the online tool miRNASNP V3 (http://bioinfo.life.hust.edu.cn/miRNASNP#!/). We excluded SNPs with minor allele frequency (MAF) < 0.05 and finally targeted two SNPs: miR-146a rs2910164 and miR-34b rs4938723. Several literatures [14, 17, 18] have reported that miR-146a rs2910164 and miR-34b rs4938723 polymorphisms are associated with CAD and ischemic stroke. However, there are few reports on the association of the two SNPs with acute coronary syndrome (ACS). Therefore, we hypothesize that miR-146a rs2910164 and miR-34b rs4938723 may also be involved in the development of ACS.

Fig. 1
figure 1

MiR-146a and miR-34b are highly expressed in CAD patients. A Principal component analysis (PCA) of GSE64566 dataset from the NCBI_GEO database (https://www.ncbi.nlm.nih.gov/gds/); B Volcano plot. Red represents up-regulated genes; blue represents down-regulated genes; black represents genes with no significant difference in expression, each point in the figure represents a gene, and the screened differential genes are genes within the range of p-value < 0.05, |FC|> 1.2; C Heatmap. Each small square in the figure represents 1 gene in 1 sample, and its color represents the expression value of the gene. The darker the color, the higher the gene expression level (red for up-regulated genes, blue for down-regulated genes). Differentially expressed genes with p-value < 0.05 and |FC|≥ 1.2 are shown in the figure; D, E. Boxplot of miR-146a and miR-34b expression

MiR-146a rs2910164, but not miR-34b rs4938723, polymorphism is a risk factor for ACS

We screened 1605 CAD patients and eventually identified 1171 patients (Fig. 2A). Among them, patients with stable angina pectoris (SA) were included in the control group (n = 552). Patients with acute ST-segment elevation myocardial infarction (STEMI), acute non-ST-segment elevation myocardial infarction (NSTEMI), and unstable angina pectoris (UA) were included in the ASC group (n = 619). Comparing baseline data between the two groups, we found that the ACS group was older, had higher rates of smoking, hypertension, and family history, and had worse cardiac and renal function (Table 1). Furthermore, serum levels of inflammatory factors IL-1β and IL-6 were significantly higher in the ACS group than in the control group (p < 0.001). Although TNF-α levels were also higher in the ACS group than in the control group, there was no statistical difference. (Fig. 2B–D).

Fig. 2
figure 2

MiR-146a rs2910164 polymorphism is a risk factor for ACS. A Flow chart of case control study; B–D Serum IL-1β, IL-6 and TNF-α levels in control group and ACS group. The unit is pg/ml. ***p < 0.001 versus Control; E–G Among all 1171 individuals, serum IL-1β, IL-6 and TNF-α levels in different genotypes of miR-146a rs2910164. The unit is pg/ml. ***p < 0.001 versus CC genotype, ### p < 0.001 versus CG genotype; H Forest plot

Table 1 Basic information of control subjects and ACS patients

The Hardy–Weinberg law of equilibrium (HWE) is the most important principle of population inheritance. In Table 2, after the HWE test was performed on both the control and the ACS group, it was confirmed that the gene distribution of miR-146a rs2910164 and miR-34b rs4938723 loci in the two groups were in line with the Hardy–Weinberg equilibrium (rs2910164: p of control group = 0.614, p of ACS group = 0.975; rs4938723: p of control group = 0.478, p of ACS group = 0.111). Subsequently, we performed a chi-square test on the miR-146a rs2910164 and miR-34b rs4938723 SNPs between the two groups and found that the miR-146a rs2910164 SNP was a risk factor for ACS in both dominant and recessive modes (Dominant model: OR = 1.270, 95%CI (1.000, 1.613), p = 0.049; Recessive model: OR = 1.402, 95%CI (1.017, 1.934), p = 0.039). Whereas, the difference of miR-34b rs4938723 SNP between the two groups was not statistically significant in both dominant and recessive modes (Table 2). In addition, as shown in Fig. 2E–G, among all 1171 patients, serum IL-1β, IL-6 and TNF-α levels were significantly higher in patients carrying the miR-146a rs2910164 G allele (CG + GG genotype) than in those carrying the miR-146a rs2910164 C allele (CC genotype). However, there was no difference in serum inflammatory factor levels among patients with miR-34b rs4938723 TT, TC and CC genotypes (Additional file 1: Figure S1 A-C). Finally, Age, gender, smoking, hypertension, diabetes, LVEF, serum CRP, TC, LDLC and HDLC levels, rs2910164 CG and rs2910164 GG were included to construct a multifactorial logistic regression equation. The results found that: men were more likely to develop ACS than women (OR = 1.438, 95%CI (1.032, 2.004), p = 0.032); diabetes increased the risk of ACS (OR = 1.602, 95%CI (1.109, 2.314), p = 0.012); higher serum CRP and TC levels were associated with a higher risk of ACS (CRP: OR = 2.348, 95%CI (2.010, 2.626), p < 0.001; TC: OR = 1.486, 95%CI (1.125, 1.963), p = 0.005); higher serum HDLC levels were associated with a lower risk of ACS (OR = 0.138, 95%CI (0.058, 0.331), p < 0.001); and most importantly rs2910164 GG remained a risk factor for ACS (OR = 1.820, 95%CI (1.119, 2.961), p = 0.016) (Fig. 2H). These results suggest that miR-146a rs2910164 polymorphism is a risk factor for ACS, and patients carrying miR-146a rs2910164 G allele have higher levels of inflammation.

Table 2 Association of miRNA polymorphisms with ACS Risk

MiR-146a rs2910164 polymorphism is associated with inflammatory factor levels and the severity of coronary lesions in the ACS group

We examined the relationship between the above two SNPs and inflammatory factor levels within the control and ACS groups, respectively, and found that the miR-146a rs2910164 polymorphism was weakly associated with inflammatory factor levels in the control group (Fig. 3A–C). Control individuals with GG genotype had slightly higher TNF-α serum level than CC and CG genotypes. However, IL-1β and IL-6 serum levels did not differ among individuals with the three genotypes (Fig. 3A–C). Surprisingly, within the ACS group, the three serum inflammatory factor levels were significantly higher in patients with the CG and GG genotypes than in patients with the CC genotype. And the TNF-α level of patients with GG genotype was higher than that of CG genotype patients (Fig. 3D–F). The miR-34b rs4938723 polymorphism was not associated with inflammatory factor levels in either the control or ACS groups (Additional file 1: Figure S1D–I).

Fig. 3
figure 3

MiR-146a rs2910164 polymorphism is associated with inflammatory factor levels and the severity of coronary lesions in the ACS group. A–C Within control group, serum IL-1β, IL-6 and TNF-α levels in different genotypes of miR-146a rs2910164. The unit is pg/ml. ***p < 0.001 versus CC genotype, ### p < 0.001 versus CG genotype; D, E Within ACS group, serum IL-1β, IL-6 and TNF-α levels in different genotypes of miR-146a rs2910164. The unit is pg/ml. ***p < 0.001 versus CC genotype, ### p < 0.001 versus CG genotype; F–G Within ACS group, Gensini scores of patients with different genotypes of miR-146a rs2910164. **p < 0.01 and ***p < 0.001 versus CC genotype, ### p < 0.001 versus CG genotype

Gensini score takes full account of the number, location and extent of coronary lesions and is a commonly used and relatively scientific evaluation standard based on angiography to assess the severity of coronary atherosclerosis. Within the ACS group, the mean Gensini scores of patients with miR-146a rs2910164 CC, CG, and GG genotypes were 52.46 ± 2.11, 63.31 ± 2.35, and 84.89 ± 4.00, respectively. As shown in Fig. 3G, H, patients carrying the G allele had significantly higher Gensini scores: the score of patients with CG genotype was significantly higher than that of patients with CC genotype, and it was further increased in patients with GG genotype. Taken together, it appears that patients carrying the G allele had higher levels of inflammation and more severe coronary lesions.

Polymorphism of miR-146a rs2910164 was associated with the outcome of post-PCI patients

The above study found that miR-146a rs2910164 polymorphism was associated with levels of inflammatory factors only in the ACS group, and that patients carrying the G allele had more severe coronary lesions, so we hypothesized that it might be related to the prognosis of post-PCI patients. We screened 697 post-PCI patients, and 612 patients were finally included in the validation cohort due to inclusion and exclusion criteria and loss to follow-up (Fig. 4A). They were categorized into CC, CG, and GG groups according to miR-146a rs2910164 genotypes, and were followed up for 14 to 60 months (median follow-up time is 42 months). Baseline data showed that patients with CG genotype were older than patients with CC genotype. GG genotype patients had higher lipoprotein a (Lpa) levels than CC genotype patients, and their fibrinogen (FIB) levels were higher than CG genotype patients (Table 3). The rest of the medical history, biochemical indicators and medication history did not differ significantly among patients with different genotypes (Table 3). Proportional hazards model was used to detect the correlation between miR-146a rs2910164 polymorphism and MACE. Table 4 showed that in the dominant model post-PCI patients with G allele were more prone to have MACE (HR = 1.405 95% CI (1.018–1.939), P = 0.038). Kaplan–Meier survival curves also indicated that patients with the GG genotype were most likely to have MACE compared to those with the CG and CC genotypes (Fig. 4B, C, Additional file 1: Figure S2A, B). In addition, Kaplan–Meier curves for mortality indicated that patients carrying G allele had higher mortality compared with patients carrying C allele (Fig. 4D, E Additional file 1: Figure S2C, D). In brief, these data indicate that miR-146a rs2910164 polymorphism was associated with the outcome of post-PCI patients.

Fig. 4
figure 4

Kaplan–Meier curve of clinical outcomes of post-PCI patients. A Flow chart of validation cohort; B, C Kaplan–Meier survival curve for MACE judging by miR-146a rs2910164 polymorphism. MACE: major adverse cardiovascular events; D, E Kaplan–Meier survival curve for death judging by miR-146a rs2910164 polymorphism

Table 3 Basic information of patients with different miR-146a rs2910164 genotypes
Table 4 Association of miR-146a rs2910164 polymorphisms with MACE

The G allele of miR-146a rs2910164 is susceptible to oxidation

Guanine can be hydroxylated by reactive oxygen species (ROS) to yield 8-oxo-7,8-dihydro-2'-deoxyguanosine (8OHdG) in DNA and 8-oxo-7,8-dihydroguanosine (8OHG) in RNA [19,20,21]. Wang, J.X et al. reported that miRNAs can also be oxidized in response to oxidative stress [22]. Oxidative stress plays a pivotal role in the process of atherosclerosis. Under such pathological conditions, ROS, including hydrogen peroxide (H2O2) and hydroxyl radicals (·OH), cause cell-damaging effects through oxidative modification of biomolecules such as proteins and nucleotides (DNA and RNA). Therefore, we infer that the G allele of miR-146a rs2910164 is susceptible to oxidation in ACS patients. To test whether the G allele can be oxidatively modified, we used an anti-8-OHG antibody that recognizes oxidatively modified 8-OHG. To avoid interference from other RNAs, miRNA fractions were purified from monocytes isolated from ACS patients and control individuals. As shown in Fig. 5A, we observed an increase in oxidatively modified miRNAs in monocytes from ACS patients. To verify whether the G allele of miR-146a rs2910164 could be oxidized to 8-OHG, we synthesized wild-type miR-146a(miR-146a(C)) and miR-146a mutated to G at rs2910164 (miR-146a(G)). They were then oxidized by the Fenton reaction and identified by HPLC–UV spectroscopy. Figure 5B shows that miR-146a(G) can be modified by oxidation. Together, these data suggest that miRNAs in ACS patients are more susceptible to oxidative modification, and miR-146a(G) can be oxidized to oxi-miR-146a(G).

Fig. 5
figure 5

Oxidized miR-146a rs2910164 mismatches with the 3’UTRs of IKBA and activates the NF-κB inflammatory pathway. A MiRNA fractions were purified from monocytes isolated from controls and ACS patients, and 8-OHG levels were determined by Northwestern blotting using an anti-8-OHG antibody; B Analysis of 8-OHG in mir-146(C) and mir-146(G) by HPLC–UV spectroscopy; C Outline of luciferase reporter assay for validating the interaction of oxi-miR-146a(G) with IKBA 3′ UTR. The miRNA response elements (MREs) of oxi-miR-146a(G) in human IKBA 3′ UTR were predicted by the bioinformatics program RNAhybrid. The red color indicates the potential mismatching sites in which 8-OHG matches with A. M1, MRE1; M2, MRE2; D Oxidized miR-146a(G) reduces the luciferase activities of IKBA 3′ UTR. HEK293 cells were transfected with the luciferase constructs of the wild-type IKBA 3′ UTR, along with miR-NC, miR-146a(C), miR-146a(G), or oxi-146a(G). The empty vector pGL3 served as a control. ∗ p < 0.05 versus IKBA 3′UTR. The results are expressed as means ± SEM of at least three independent experiments; E–F THP-1 cells were transfected with miR-146a(C) and miR-146a(G), respectively, and then treated with ox-LDL or its control. IκBα(E) and P65 (F) protein levels were analyzed by immunoblot. Quantitative analysis of IκBα and P65 levels is shown in the lower panels. *p < 0.05 and ** p < 0.01 versus Control; G, H P65 immunostaining of atherosclerotic plaque sections of CC, CG, and GG genotypes. Left, overview of plaque; middle and right, higher magnification

Oxidized miR-146a rs2910164 mismatches with the 3’UTR of IKBA and activates the NF-κB inflammatory pathway

8-OHG has been well characterized as pairing with A instead of C, resulting in changes in the targeting mRNA [23]. Through Targetscan, miRanda, miRMap, miTarBase and miRDB database predictions [24,25,26,27,28], it was found that miR-146a may interact with mRNAs related to the NF-κB signaling pathway (Additional file 1: Figure S3A, B, Additional file 3:  Excel2_miR 146a target). Besides, the above results suggest that patients carrying rs2910164 G allele have higher levels of inflammation. Therefore, we speculated that oxi-miR-146a(G) interact with mRNAs of the IKB protein family. After alignment, it was found that oxi-miR-146a(G) may be complementary to two sites in the 3'UTR of IKBA (Fig. 5C). To investigate whether these potential sites can be targeted by oxi-miR-146a(G), we made a luciferase construct of wild-type IKBA 3'UTR. The luciferase reporter assay indicated that oxi-miR-146a(G) caused a significant decrease in the luciferase activity of the IKBA 3' UTR, whereas miR-146a(C) and miR-146a(G) had no effect on the IKBA 3' UTR (Fig. 5D). Next, we transfected THP-1 with miR-146(C) and miR-146(G). In monocytes transfected with miR-146(G), a decrease in the protein expression level of IκBα and an increase in the protein expression of P65 were observed after ox-LDL treatment (Fig. 5E, F). This implies that ROS stimulated by ox-LDL oxidatively modified miR-146(G), which in turn inhibited the protein expression of IκBα, leading to an increased level of inflammation. Immunohistochemical staining of carotid plaques from atherosclerotic patients obtained after atherectomy showed that the expression of P65 protein was increased in the plaques of patients with CG and GG genotypes (Fig. 5G). This result further demonstrated that in patients with the rs2910164 G allele, miR-146a rs2910164 mismatches with the IKBA 3ʹUTRs and the intra-plaque NF-κB inflammatory pathways were activated.

Discussion

This research includes a case–control study and a validation cohort. We first identified miR-146 rs2910164, but not miR-34b rs4938723, associated with the development of ACS using a case–control approach, and then followed with a validation cohort to determine a causal relationship between miR-146 rs2910164 and outcomes of post- PCI patient. The G allele of miR-146a rs2910164 was an independent predictor of MACE for ACS patients after PCI. The results provided new evidence for the relationship between miRNAs polymorphisms and both the onset of ACS and long-term prognosis of ACS patients.

The C allele of miR-146a rs2910164 has been reported to be a risk allele for nasopharyngeal carcinoma, while the G allele is a risk allele for asthma and vascular disease [29,30,31]. In our current case–control study, patients carrying the G allele of miR-146a rs2910164 had 1.270-fold (dominant model) or 1.402-fold (recessive model) higher risk of ACS than those with the C allele, indicating that the G allele of miR-146a rs2910164 is a susceptible allele for ACS. In addition, compared to C allele carriers, G allele carriers had a 1.405-fold(dominant model) higher incidence of MACE after PCI during the 5-year follow-up period. Moreover, the G allele of miR-146a rs2910164 was associated with serum inflammation factor levels.

Although several studies have shown that the miR-34 family is associated with the development of cardiovascular diseases such as coronary atherosclerosis and acute myocardial infarction [32,33,34]. However, our study revealed that polymorphism of miR-34b rs4938723 was not associated with the incidence of ACS and serum inflammatory factor levels.

Several articles have demonstrated that miR-146a inhibits inflammation levels by targeting the upstream NF-κB signaling proteins: interleukin-1 receptor-associated kinase 1 (IRAK1) and TNF receptor-associated factor 6 (TRAF6) through a negative feedback mechanism [35,36,37]. However, we observed that patients with miR-146a rs2910164 variant had higher levels of inflammation, so we speculate that perhaps the rs2910164 polymorphism alters the original target of miR146a and its physiological function.

Guanines in miRNAs can be hydroxylated by reactive oxygen species (ROS) to generate 8-oxo-7,8-dihydroguanosine (8OHG) [19,20,21,22]. 8-OHG pairs with A instead of C, resulting in changes in the miRNA-mRNA pairing [23].We first detected higher levels of miRNA oxidation, including miR146a, in monocytes from ACS patients. Secondly, it was demonstrated by HPLC–UV spectroscopy that mutated mir-146(miR-146a(G)), but not wild-type 146(miR-146a(C)), can be oxidatively modified. Oxi-miR-146a(G) activates the NF-κB pathway by mismatching with 3ʹUTR of IKBA. These results suggest that the pathogenicity of the miR-146a rs2910164 G allele is mediated in part by its mismatch with IKBA and its regulatory role in the inflammatory response. This conclusion is further supported by immunostaining for P65 in human atherosclerotic plaques.

Interestingly, some miRNAs were downregulated in CAD patients during our bioinformatics analysis. The downregulated genes are sorted according to the fold difference from large to small: miR-589, miR-520h, miR-144, miR-518f, miR-379, miR-512-5p, miR-668 and miR-554. Most of these miRNAs play roles in cancer, while there are few studies on them related to cardiovascular diseases.

MiR-520h can promote tumor growth and progression in hepatocellular carcinoma cell and breast cancer. Meanwhile, miR-520 h can target smad7/EMT, promoting the metastasis of bladder cancer cells [38]. Similarly, miR-518f-5p can reduce the expression of the tetraspanin CD9 in prostate cancer cells [39]. Danielle R Bond et al. also verified that the transfection of miR-518f-5p significantly decreased CD9 protein expression and increased breast cell migration in vitro [40]. These studies show that miR-518f-5p can enhance the migration and invasion ability of prostate cancer cells and breast cancer cells. As of miR-554, Yasuo Takashima et al. found that miR-554 is designated as a miRNA predictor in PCNSL (primary central nervous system lymphoma), which is associated with all categories such as Th-1/Th-2 and T-reg differentiation status, and stimulatory and inhibitory immune checkpoints [41]. MiR-554 has been identified as a reliable prognostic predictor of PCNSL.

However, miR-512-5p plays a significant inhibitory role in tumor progression. MiR-512-5p can target a variety of different genes, inhibit the expression of tumor related genes, and reduce the tumor cell’s ability of proliferation, migration, and invasion. Yixiang Huang et al. found that in breast cancer, circular RNA circRPPH1 acted as a miRNA sponge of miR-512-5p, which inhibits the combination of miR-512-5p and 3 '- UTR of STAT1, reduces the expression of STAT1, and promotes the he proliferation, migration and colony formation ability of Breast cancer cells [42].miR-668 plays a protective role against renal injury after ischemia–reperfusion (IRI). In ischemic acute renal injury, HIF1 induces miR-668 during the ischemic reperfusion (IRI) interval. This microRNA can inhibit the expression of mitochondrial protein 18 kDa (MTP18), maintaining mitochondrial dynamics and promoting the survival of renal tubular cells after IRI [43].

In acute myocardial infarction (AMI), miR-379 binds tumor necrosis factor-α-induced protein 8 (TNFAIP8) and inhibits its activity significantly, leading to the reduce of the level of apoptosis and the inhibition NK-κB signal pathway in H9c2 cells. This study shows that miR-379 can play a protective role in the heart [44].

MiR-144 can be used as a sequence specific inhibitor of ABCA1/G1 (ATP binding cassette transporter) mRNA to control gene expression. MiR-144 binds to the 3'UTR region site of ABCA1/G1.This will reduce the outflow of cholesterol to apolipoprotein A1 (ApoA1) mediated by the ABCA1/G1 of hepatocytes and macrophages, promoting inflammatory activation and atherosclerosis [45].

Studies have shown that COX2 promoter RNA contains two miR-589-5p binding sites. MiR-589 is an endogenous activator in the nucleus that can regulate the expression of COX-2 and PLA2G4A genes at the transcriptional level. Activation of the COX-2 gene by miR-589 requires AGO2 and GW182. COX-2 catalyzes the conversion of arachidonic acid and further metabolizes it to prostaglandins and other arachidonic acids. These arachidonic acids mediate many biological processes, including inflammation, cellular immunity, and cancer [46].

There are several limitations in our present study. Firstly, our data did not provide the expression levels of 3 different isoforms of miR-146a under the current circumstances. Secondly, the content of serum cytokines was measured at the time of admission and were not monitored dynamically. Thirdly, the validation cohort should be collected at different hospitals, or even multiple hospitals in different regions.

Conclusion

The rs2910164 variant of miR-146a is closely associated with the risk of ACS in Chinese Han population. The patients carrying miR-146a with G allele inclines to be oxidized and may have a more serious pathological change and poorer post-PCI prognosis, partly due to the oxidatively modified miR-146a mispairing with 3’UTR of IKBA and activating NF-κB inflammatory pathways.

Availability of data and materials

The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

ACS:

Acute coronary syndromes

PCI:

Percutaneous coronary intervention

miRNAs:

MicroRNAs

MACE:

Major adverse cardiovascular events

MI:

Myocardial infarction

SA:

Stable angina

STEMI:

ST-segment elevation myocardial infarction

NSTEMI:

Non-ST-segment elevation myocardial infarction

UA:

Unstable angina

CI:

Confidence intervals

HR:

Hazard ratio

OR:

Odds ratio

CAD:

Coronary artery disease

EAT:

Epicardial adipose tissues

PCA:

Principal component analysis

HWE:

Hardy–Weinberg law of equilibrium

SNPs:

Single nucleotide polymorphisms

8OHdG:

8-Oxo-7,8-dihydro-2'-deoxyguanosine

8OHG:

8-Oxo-7,8-dihydroguanosine

ROS:

Reactive oxygen species

H2O2:

hydrogen peroxide

·OH:

Hydroxyl radical

IRAK1:

Interleukin-1 receptor-associated kinase 1

TRAF6:

TNF receptor-associated factor 6

BMI:

Body mass index

SBP:

Systolic blood pressure

DBP:

Diastolic blood pressure

LVEF:

Left ventricular ejection fraction

FBG:

Fasting blood glucose

CK:

Creatine kinase

CKMB:

Creatine kinase isoenzymes

FIB:

Fibrinogen

FDP:

Fructose-1, 6-diphosphate

CRP:

C-reactive protein

BNP:

Brain natriuretic peptide

TC:

Total cholesterol

TG:

Triglyceride

HDL:

High-density lipoprotein cholesterol

LDL:

Low-density lipoprotein cholesterol

Apo:

Apolipoprotein

Lpa:

Lipoprotein a

HGB:

Hemoglobin

WBC:

White blood cell count

ACEI:

Angiotensin converting enzyme inhibitor

ARB:

Angiotensin receptor blocker

CCB:

Calcium channel blocker

References

  1. Collaborators, G.B.D.C.o.D. Global regional, and national age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980–2017: a systematic analysis for the global burden of disease study 2017. Lancet. 2018;392(10159):1736–88.

    Article  Google Scholar 

  2. Diseases GBD, Injuries C. Global burden of 369 diseases and injuries in 204 countries and territories, 1990–2019: a systematic analysis for the global burden of disease study 2019. Lancet. 2020;396(10258):1204–22.

    Article  Google Scholar 

  3. Fuster V, et al. The pathogenesis of coronary artery disease and the acute coronary syndromes (2). N Engl J Med. 1992;326(5):310–8.

    Article  CAS  PubMed  Google Scholar 

  4. Wandzioch E, Zaret KS. Dynamic signaling network for the specification of embryonic pancreas and liver progenitors. Science. 2009;324(5935):1707–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Plasterk RH. Micro RNAs in animal development. Cell. 2006;124(5):877–81.

    Article  CAS  PubMed  Google Scholar 

  6. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116(2):281–97.

    Article  CAS  PubMed  Google Scholar 

  7. Jones Buie JN, et al. The role of miRNAs in cardiovascular disease risk factors. Atherosclerosis. 2016;254:271–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Raitoharju E, et al. miR-21, miR-210, miR-34a, and miR-146a/b are up-regulated in human atherosclerotic plaques in the tampere vascular study. Atherosclerosis. 2011;219(1):211–7.

    Article  CAS  PubMed  Google Scholar 

  9. Oerlemans MI, et al. Early assessment of acute coronary syndromes in the emergency department: the potential diagnostic value of circulating microRNAs. EMBO Mol Med. 2012;4(11):1176–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Zhang X, et al. Expression profiles of six circulating microRNAs critical to atherosclerosis in patients with subclinical hypothyroidism: a clinical study. J Clin Endocrinol Metab. 2014;99(5):E766–74.

    Article  CAS  PubMed  Google Scholar 

  11. Macauda A, et al. Identification of miRSNPs associated with the risk of multiple myeloma. Int J Cancer. 2017;140(3):526–34.

    Article  CAS  PubMed  Google Scholar 

  12. Gazouli M, et al. Association study of genetic variants in miRNAs in patients with inflammatory bowel disease: preliminary results. Dig Dis Sci. 2013;58(8):2324–8.

    Article  CAS  PubMed  Google Scholar 

  13. Xu Y, et al. A potentially functional polymorphism in the promoter region of miR-34b/c is associated with an increased risk for primary hepatocellular carcinoma. Int J Cancer. 2011;128(2):412–7.

    Article  PubMed  Google Scholar 

  14. Bao MH, et al. Meta-Analysis of miR-146a polymorphisms association with coronary artery diseases and ischemic stroke. Int J Mol Sci. 2015;16(7):14305–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Song H, et al. A genome-wide association study identifies a new ovarian cancer susceptibility locus on 9p22.2. Nat Genet. 2009;41(9):996–1000.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Vacca M, et al. Integrative miRNA and whole-genome analyses of epicardial adipose tissue in patients with coronary atherosclerosis. Cardiovasc Res. 2016;109(2):228–39.

    Article  CAS  PubMed  Google Scholar 

  17. Wang Y, et al. Two single nucleotide polymorphisms (rs2431697 and rs2910164) of miR-146a are associated with risk of coronary artery disease. Int J Environ Res Public Health. 2017;14(5):514.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Cen J, et al. Pri-miR-34b/c rs4938723 polymorphism contributes to coronary artery disease susceptibility. Arch Med Res. 2019;50(4):170–4.

    Article  CAS  PubMed  Google Scholar 

  19. Kasai H, et al. 8-Hydroxyguanine, a DNA adduct formed by oxygen radicals: its implication on oxygen radical-involved mutagenesis/carcinogenesis. J Toxicol Sci. 1991;16(Suppl 1):95–105.

    Article  CAS  PubMed  Google Scholar 

  20. Maki H, Sekiguchi M. MutT protein specifically hydrolyses a potent mutagenic substrate for DNA synthesis. Nature. 1992;355(6357):273–5.

    Article  CAS  PubMed  Google Scholar 

  21. Taddei F, et al. Counteraction by MutT protein of transcriptional errors caused by oxidative damage. Science. 1997;278(5335):128–30.

    Article  CAS  PubMed  Google Scholar 

  22. Wang JX, et al. Oxidative modification of miR-184 enables It to target Bcl-xL and Bcl-w. Mol Cell. 2015;59(1):50–61.

    Article  CAS  PubMed  Google Scholar 

  23. Cheng KC, et al. 8-hydroxyguanine, an abundant form of oxidative DNA damage, causes G––T and A––C substitutions. J Biol Chem. 1992;267(1):166–72.

    Article  CAS  PubMed  Google Scholar 

  24. Agarwal V, et al. Predicting effective microRNA target sites in mammalian mRNAs. Elife. 2015. https://doi.org/10.7554/eLife.05005.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Vejnar CE, Zdobnov EM. MiRmap: comprehensive prediction of microRNA target repression strength. Nucleic Acids Res. 2012;40(22):11673–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Enright AJ, et al. MicroRNA targets in drosophila. Genome Biol. 2003;5(1):R1.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Huang HY, et al. miRTarBase update 2022: an informative resource for experimentally validated miRNA-target interactions. Nucleic Acids Res. 2022;50(D1):D222–30.

    Article  CAS  PubMed  Google Scholar 

  28. Chen Y, Wang X. miRDB: an online database for prediction of functional microRNA targets. Nucleic Acids Res. 2020;48(D1):D127–31.

    Article  CAS  PubMed  Google Scholar 

  29. Jeon YJ, et al. Association of the miR-146a, miR-149, miR-196a2, and miR-499 polymorphisms with ischemic stroke and silent brain infarction risk. Arterioscler Thromb Vasc Biol. 2013;33(2):420–30.

    Article  CAS  PubMed  Google Scholar 

  30. Lung RW, et al. A single nucleotide polymorphism in microRNA-146a is associated with the risk for nasopharyngeal carcinoma. Mol Carcinog. 2013;52(Suppl 1):E28-38.

    Article  CAS  PubMed  Google Scholar 

  31. Su XW, et al. Association between single-nucleotide polymorphisms in pre-miRNAs and the risk of asthma in a Chinese population. DNA Cell Biol. 2011;30(11):919–23.

    Article  CAS  PubMed  Google Scholar 

  32. Bernardo BC, et al. Therapeutic inhibition of the miR-34 family attenuates pathological cardiac remodeling and improves heart function. Proc Natl Acad Sci USA. 2012;109(43):17615–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Iekushi K, et al. Regulation of cardiac microRNAs by bone marrow mononuclear cell therapy in myocardial infarction. Circulation. 2012;125(14):1765–73.

    Article  CAS  PubMed  Google Scholar 

  34. Rokavec M, et al. The p53/miR-34 axis in development and disease. J Mol Cell Biol. 2014;6(3):214–30.

    Article  CAS  PubMed  Google Scholar 

  35. Bhaumik D, et al. Expression of microRNA-146 suppresses NF-kappaB activity with reduction of metastatic potential in breast cancer cells. Oncogene. 2008;27(42):5643–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Garo LP, et al. MicroRNA-146a limits tumorigenic inflammation in colorectal cancer. Nat Commun. 2021;12(1):2419.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Su YL, et al. Myeloid cell-targeted miR-146a mimic inhibits NF-kappaB-driven inflammation and leukemia progression in vivo. Blood. 2020;135(3):167–80.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Ren L, et al. Nitroxoline suppresses metastasis in bladder cancer via EGR1/circNDRG1/miR-520h/smad7/EMT signaling pathway. Int J Biol Sci. 2022;18(13):5207–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Bond DR, et al. miR-518f-5p decreases tetraspanin CD9 protein levels and differentially affects non-tumourigenic prostate and prostate cancer cell migration and adhesion. Oncotarget. 2018;9(2):1980–91.

    Article  PubMed  Google Scholar 

  40. Bond DR, et al. Tetraspanin CD9 is regulated by miR-518f-5p and functions in breast cell migration and in vivo tumor growth. Cancers. 2020;12(4):e0229577.

    Article  Google Scholar 

  41. Takashima Y, et al. miR-101, miR-548b, miR-554, and miR-1202 are reliable prognosis predictors of the miRNAs associated with cancer immunity in primary central nervous system lymphoma. PLoS ONE. 2020;15(2):e0229577.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Huang Y, et al. Circular RNA circRPPH1 promotes breast cancer progression via circRPPH1-miR-512-5p-STAT1 axis. Cell Death Discov. 2021;7(1):376.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Wei Q, et al. MicroRNA-668 represses MTP18 to preserve mitochondrial dynamics in ischemic acute kidney injury. J Clin Invest. 2018;128(12):5448–64.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Dai H, et al. MiR-379 relieves myocardial injury after acute myocardial infarction by regulating tumor necrosis factor-alpha-induced protein 8. Panminerva Med. 2022;64(1):31–7.

    Article  PubMed  Google Scholar 

  45. Ramirez CM, et al. Control of cholesterol metabolism and plasma high-density lipoprotein levels by microRNA-144. Circ Res. 2013;112(12):1592–601.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Matsui M, et al. Promoter RNA links transcriptional regulation of inflammatory pathway genes. Nucleic Acids Res. 2013;41(22):10086–109.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This study was financially supported by the National Natural Science Foundation of China (81870330). Funds were used to cover investigator costs, determine the genotypes of miR-146a rs2910164 and miR-34b rs4938723, and to purchase reagents and consumables.

Author information

Authors and Affiliations

Authors

Contributions

TC and DZ designed this project. X-RQ was responsible for statistical analysis and cellular experiments. TZ collected patients’ data. YX and XY were in charge of detecting serum inflammatory factors. TC, YW and ZY were responsible for reviewing and validating all tables, figures and reagents. X-RQ and DZ were major contributors in writing the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Dong Zhou or Tao Chen.

Ethics declarations

Ethics approval and consent to participate

Written informed consent was obtained from all study participants, and this study was approved by the ethnic committee approval of the First Affiliated Hospital of Xi’an Jiaotong University.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Supplementary material.

Additional file 2.

 Excel1_Different miRNA.

Additional file 3.

Excel2_miR 146a target.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Qiao, XR., Zheng, T., Xie, Y. et al. MiR-146a rs2910164 (G/C) polymorphism is associated with the development and prognosis of acute coronary syndromes: an observational study including case control and validation cohort. J Transl Med 21, 325 (2023). https://doi.org/10.1186/s12967-023-04140-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-023-04140-4

Keywords