Abstract Background Mounting evidence suggests that circulating microRNAs (miRNAs) are critical indicators of cardiovascular disease. However, prospective studies linking circulating miRNAs to incident acute coronary syndrome (ACS) are limited, and the underlying effect of associated miRNA on incident ACS remains unknown. Methods and Results Based on a 2‐stage prospective nested case–control design within the Dongfeng‐Tongji cohort, we profiled plasma miRNAs from 23 pairs of incident ACS cases and controls by microarray and validated the candidate miRNAs in 572 incident ACS case–control pairs using quantitative real‐time polymerase chain reaction. We observed that plasma miR‐4286 was associated with higher risk of ACS (adjusted odds ratio according to an interquartile range increase, 1.26 [95% CI, 1.07–1.48]). Further association analysis revealed that triglyceride was positively associated with plasma miR‐4286, and an interquartile range increase in triglyceride was associated with an 11.04% (95% CI, 3.77%–18.83%) increase in plasma miR‐4286. In addition, the Mendelian randomization analysis suggested a potential causal effect of triglyceride on plasma miR‐4286 (β coefficients: 0.27 [95% CI, 0.01–0.53] and 0.27 [95% CI, 0.07–0.47] separately by inverse variance‐weighted and Mendelian randomization‐pleiotropy residual sum and outlier tests). Moreover, the causal mediation analysis indicated that plasma miR‐4286 explained 5.5% (95% CI, 0.7%–17.0%) of the association of triglyceride with incident ACS. Conclusions Higher level of plasma miR‐4286 was associated with an increased risk of ACS. The upregulated miR‐4286 in plasma can be attributed to higher triglyceride level and may mediate the effect of triglyceride on incident ACS. Keywords: acute coronary syndrome, microRNA, miR‐4286, prospective study, triglyceride Subject Categories: Epidemiology, Acute Coronary Syndromes, Epigenetics __________________________________________________________________ Nonstandard Abbreviations and Acronyms FRS Framingham Risk Score IR insulin resistant miRNA microRNA MR Mendelian randomization UA unstable angina Clinical Perspective What Is New? * In our study, we found that plasma miR‐4286 was associated with higher risk of ACS based on a prospective nested case–control design. * Further Mendelian randomization analysis revealed that triglyceride was causally associated with a higher level of plasma miR‐4286. * Plasma miR‐4286 explained 5.5% of the association of triglyceride with incident ACS, suggesting the mediation effect of miR‐4286 in ACS development. What Are the Clinical Implications? * These findings identified a novel risk contributor to incident ACS at the circulating miRNA level and revealed its mediation potential in ACS development, suggesting a potential target for intervention or prevention of ACS. Coronary heart disease (CHD) is a major disease burden worldwide.[60] ^1 , [61]^2 Acute coronary syndrome (ACS), including acute myocardial infarction (AMI) and unstable angina (UA), is the main fatal subtype of CHD.[62] ^3 To date, the precise prevention and management of ACS remains a challenge because of its complicated pathogenic and physiological processes. Previous studies have demonstrated that the development of endothelial dysfunction and subsequent vulnerable plaque rupture or erosion is a major pathophysiological process of ACS.[63] ^4 However, the original regulator in ACS progression remains to be elucidated, especially from an epigenetic aspect. MicroRNAs (miRNAs) are a class of small noncoding RNAs that are involved in posttranscriptional regulation by binding target mRNAs and play important roles in cardiac remodeling, angiogenesis, and vessel development.[64] ^5 Given the feature of high stability and sensitivity, circulating miRNAs have been widely characterized as diagnostic and prognostic biomarkers for cardiovascular disease (CVD).[65] ^6 , [66]^7 However, prospective clues on circulating miRNAs and risk of ACS are limited and based on Europeans.[67] ^8 , [68]^9 Over the past decades, with the change of lifestyles and westernization of the diet, CVD has been the leading cause of death in China, and among cardiovascular deaths, the majority are attributed to ACS.[69] ^10 Thereby, prospective associations of circulating miRNAs with incident ACS and its risk factors among the Chinese need to be investigated. It is now well‐established that hyperlipidemia is a major risk factor for fatal and nonfatal CHD.[70] ^11 , [71]^12 With the development of the genome‐wide association study, Mendelian randomization (MR) has become a useful approach to estimating the causal relationship between exposure and outcome.[72] ^13 Recent MR studies have demonstrated that abnormal plasma lipid levels were causally associated with CHD, indicating that the genetic variants of lipid traits played a key role in the development of atherosclerosis and vulnerable plaque.[73] ^14 , [74]^15 Previous reviews have summarized that liver‐ or vessel‐enriched miRNAs are involved in regulating cholesterol homeostasis.[75] ^16 , [76]^17 Additionally, current studies have identified several circulating miRNAs that are related to lipids or are impacted by lipid‐related dietary interventions.[77] ^18 , [78]^19 However, there is still uncertainty regarding whether the relationships are causal between aberrant circulating miRNAs and lipid levels, as well as the mediation potential of associated miRNAs in ACS progression. Herein, we conducted a 2‐stage nested case–control study in the Dongfeng‐Tongji cohort to investigate the associations of plasma miRNAs with incident ACS, and of the associated target miRNA with cardiovascular traits and risk factors. For lipid traits associated with target miRNA, we further evaluated their causal relationships and estimated the mediation effects of target miRNA on the associations between lipids and incident ACS. An overview of the study design is presented in Figure [79]S1. Methods Details of methods are provided in Data [80]S1. Data and materials supporting the conclusions of the article are available from the corresponding author on reasonable request. Study Design and Participants This 2‐stage prospective nested case–control study was performed in the Dongfeng‐Tongji cohort,[81] ^20 which recruited 27 009 retirees of Dongfeng Motor Corporation during September 2008 and June 2010 in Shiyan, China. The first follow‐up was completed in October 2013 and 14 120 individuals were newly recruited. All of the study participants were free of ACS at the baseline of the study (2013). During follow‐up, 595 cases were newly diagnosed ACS during follow‐up through December 31, 2016, including 147 AMI and 448 UA cases. We randomly selected a control who was free of CVD and cancer at the time of the case event to match each incident ACS case using the incidence density sampling method.[82] ^21 Each pair was matched on age (±3 years), sex, and blood drawing time (±1 month). All individuals had completed questionnaires and physical examinations, and provided blood samples in 2013. In the discovery stage, we randomly selected 23 pairs of plasma samples to conduct miRNA microarray detection. In the validation stage, the remaining 572 pairs of plasma samples were used for validation by quantitative real‐time polymerase chain reaction. A detailed sampling flowchart is provided in Figure [83]S2A. The study was approved by the Ethics and Human Subject Committee of Tongji Medical College (2012‐10) and Dongfeng General Hospital. Each participant provided written informed consent. Discovery and Validation of Candidate miRNAs in Plasma We used 400 µL plasma to extract RNA by the miRNeasy Serum/Plasma kit (Qiagen, Germany) according to the instruction handbook. Using the Agilent Human miRNA Microarray, Release 21.0, 8x60K (Agilent Technologies Inc, USA), we selected differentially expressed miRNAs by paired t tests in the discovery stage. Candidate miRNAs reaching the significant level (P<0.05 and absolute fold change >1.3, no multiple testing correction) were included for quantitative real‐time polymerase chain reaction validation. Given that the present study was prospective and the cases were relatively healthy status at baseline, and the expression pattern of circulating miRNAs was complex because of its various sources, the miRNA profiles might not be dramatically altered in comparison with controls at baseline. Therefore, we used the uncorrected P value as screening criterion in line with previously published prospective or circulating miRNA‐based studies.[84] ^9 , [85]^22 We quantitated candidate miRNAs by TaqMan Advanced miRNA assay according to the protocols of the manufacturer (Applied Biosystems Inc, USA), and the reactions were carried out by QuantStudio 7 Flex Real‐Time PCR System (Applied Biosystems Inc, USA). The miRNA levels were measured by the 2^−ΔCt method. Bidirectional 2‐Sample MR Analyses We screened 39 and 58 independent loci separately associated with triglyceride and high‐density lipoprotein cholesterol (HDL‐C) with weak linkage disequilibrium (r ^2<0.2) as genetic instrumental variables based on the previous reports on Asians (P<5E‐08). For instrumental variables of miR‐4286, we selected 24 single nucleotide polymorphisms of miR‐4286 based on the genome‐wide association study findings from 340 participants in the Dongfeng‐Tongji study. A detailed sampling flowchart is provided in Figure [86]S2B. The proportion of exposure explained by the genetic instruments (adjusted R ^2) and F statistics were computed using linear models that fitted the weighted genetic risk score of exposure and adjusted for age and sex. To estimate causal relationships between lipids and miR‐4286, we used the bidirectional 2‐sample MR analyses by the “TwoSampleMR” and “MRPRESSO” packages, which provide inverse variance–weighted, MR‐Egger, weighted median, and MR‐PRESSO methods.[87] ^23 , [88]^24 Statistical Analysis The power estimation for the discovery stage was calculated by the PASS software (version 11.0.10). In the situation of setting 23 pairs of samples, 81% power was obtained when anticipating to achieve minimum significant condition (fold change=1.3, P=0.05) based on the test of eligible 408 miRNAs on the array after data filtering. Baseline characteristics of the case and control groups in 2 stages were compared by t tests or Mann–Whitney U tests for continuous variables, and χ ^2 tests or Fisher exact tests for categorical variables. Data of plasma miRNAs were natural log‐transformed to improve normality. To investigate the association of candidate miRNAs with incident ACS, we categorized the individual miRNAs into tertiles based on their distribution among the controls and fitted conditional logistic regressions to estimate the odds ratio (OR) and 95% CI for incident ACS. Similarly, we also calculated the OR (95% CI) for incident ACS per interquartile range (IQR) increase in the miRNA level. The multivariable model adjusted for established potential confounders including age, body mass index (BMI), smoking status, drinking status, education level, metabolic equivalent, family history of CHD, diabetes mellitus, hypertension, HDL‐C, low‐density lipoprotein cholesterol, triglyceride, and use of lipid‐lowering medication. Missing data of categorical covariates were replaced with an extra category to indicate missingness. Linear trend P value was estimated by assigning the median value of miRNA to each tertile and using this as a continuous variable in the logistic regression model. We defined the statistical significance as a false discovery rate (FDR) <0.1 using the Benjamini‐Hochberg method.[89] ^25 We further estimated the risk difference that represented the absolute risk change for ACS attributed to plasma miR‐4286 using generalized linear regression models assuming an identity link function.[90] ^26 We also used elastic net regression models with 10‐fold cross‐validation to test the robustness between miRNAs and incident ACS, which combines the Lasso and Ridge penalties to reduce the occurrence of false‐positives.[91] ^27 Sensitivity analysis was conducted by excluding ACS events that had been diagnosed within 1 year after collection of blood samples in 2013. In addition, we conducted restricted cubic splines to estimate the associations between target miRNA and incident ACS. We also performed stratified analyses according to age, sex, BMI, current smokers, current drinkers, hypertension, diabetes mellitus, fasting glucose level, and hyperlipidemia. To evaluate the predictive value of ACS‐associated miRNA, we compared the area under the receiver‐operating characteristic curves between the Framingham Risk Score (FRS) for hard CHD and addition of selected miRNA to the FRS, and also calculated the net reclassification improvement and the integrated discrimination improvement to assess the improvement of predictive ability. The net reclassification improvement was calculated according to low (<10%), moderate (10%–30%), and high (>30%) levels of predicted risk of ACS. To estimate the heterogeneity of association between ACS subtypes, we fitted a multinomial logistic regression model with the outcomes (AMI, UA, or non‐ACS) as the dependent variable and miRNA level as independent variable, and adjusted for age, sex, BMI, smoking status, drinking status, education levels, metabolic equivalent, diabetes mellitus, hypertension, family history of CHD, HDL‐C, low‐density lipoprotein cholesterol, triglyceride, and use of lipid‐lowering medication. P value for heterogeneity was obtained from multinomial logistic regression model in the comparison between AMI and UA. The relationships between lipid traits and target miRNA were estimated by linear regression. Plasma levels of HDL‐C, total cholesterol, triglyceride, and triglyceride/HDL‐C ratio were natural log‐transformed and included in linear models and further analyses. Multivariable models adjusted for age, sex, BMI, smoking status, drinking status, education level, metabolic equivalent, and use of lipid‐lowering drugs. Calculation of linear trend P value, and definition of statistical significance were the same as above. We presented the percent difference in miRNA level when categorizing the participants according to tertiles of each lipid level, which represented the percent change in miRNA in relation to the upper lipid group versus lower lipid group. The percent difference (95% CI) in miRNA level was also estimated in relation to per IQR increase in each lipid level.[92] ^28 Furthermore, we used the “mediation” R package,[93] ^29 which contained mediator models and outcome models, to estimate the mediation effects of target miRNA between lipids and incident ACS. The mediated proportion reflected the average mediation effect of target miRNA between lipids and incident ACS, estimated by the proportion of indirect effect in total effect. Analyses were performed with R software (version 3.4.3, R Core Team) or SAS program (version 9.4, SAS Institute). Results Characteristics of the Study Participants The cases were newly diagnosed with ACS within 0.5 to 3.6 years since the 2013 baseline, and the median follow‐up was 1.8 years. Among the participants, mean±SD age was 67.2±7.6 years, and 57.4% were men. In the validation stage, ACS cases were more likely to have a history of hypertension, with higher BMI, fasting glucose, low‐density lipoprotein cholesterol, triglyceride, and lower HDL‐C levels. By contrast, a consistent but weak group difference was detected between ACS cases and controls in the discovery stage ([94]Table). In comparison with respective controls, both AMI and UA cases had higher blood pressure, fasting glucose, and triglyceride levels. In addition, AMI cases were more likely to be less educated, and physically inactive, while UA cases had higher BMI and lower HDL‐C levels (Table [95]S1). Table 1. Baseline Characteristics of the 2‐Stage Nested Case–Control Participants Variables Discovery Stage Validation Stage ACS Cases (n=23) Controls (n=23) P Value ACS Cases (n=572) Controls (n=572) P Value Age, y 68.1±6.7 68.1±6.7 0.96 67.2±7.7 67.2±7.6 0.91 Male, n (%) 14 (60.9) 14 (60.9) 1.00 327 (57.2) 327 (57.2) 1.00 BMI, kg/m^2 24.60±3.45 23.90±3.58 0.50 24.71±3.39 24.21±3.04 0.009 SBP, mm Hg 146.83±26.00 140.22±22.18 0.36 148.21±24.44 142.78±22.37 <0.001 DBP, mm Hg 79.04±11.33 79.65±12.54 0.86 83.28±13.26 80.67±12.58 <0.001 Smoking status, n (%) 1.00 0.08 Current smoker 5 (21.7) 5 (21.7) 136 (23.8) 105 (18.4) Former smoker 4 (17.4) 4 (17.4) 119 (20.8) 110 (19.2) Never smoker 14 (60.9) 14 (60.9) 315 (55.1) 355 (62.1) Drinking status, n (%) 1.00 0.11 Current drinker 9 (39.1) 10 (43.5) 158 (27.6) 185 (32.3) Former drinker 2 (8.7) 1 (4.3) 75 (13.1) 56 (9.8) Never drinker 12 (52.2) 12 (52.2) 336 (58.7) 330 (57.8) Education level, n (%) 0.87 0.22 Primary school or below 7 (30.4) 8 (34.8) 175 (30.6) 154 (26.9) Middle school 10 (43.5) 8 (34.8) 218 (38.1) 212 (37.1) High school or beyond 6 (26.1) 7 (30.4) 174 (30.4) 196 (34.3) Physical activity, MET‐h/wk 21.00 (12.00, 44.00) 21.00 (14.25, 42.00) 0.61 21.00 (9.25, 42.00) 23.31 (13.88, 42.00) 0.05 Family history of CHD, n (%) 2 (8.7) 1 (4.3) 0.61 39 (6.8) 41 (7.2) 0.91 Lipid‐lowering medication, n (%) 1 (4.3) 1 (4.3) 1.00 85 (14.9) 82 (14.3) 0.87 Antihypertensive medication, n (%) 12 (52.2) 7 (30.4) 0.23 231 (40.4) 193 (33.7) 0.02 Antidiabetic medication, n (%) 4 (17.4) 3 (13.0) 1.00 76 (13.3) 57 (10.0) 0.10 Hyperlipidemia, n (%) 11 (47.8) 6 (26.1) 0.22 380 (66.4) 374 (65.4) 0.76 Hypertension, n (%) 16 (69.6) 15 (65.2) 1.00 449 (78.5) 409 (71.5) 0.01 Diabetes mellitus, n (%) 9 (39.1) 4 (17.4) 0.19 160 (28.0) 134 (23.4) 0.09 FG, mmol/L 5.90 (5.11, 6.54) 5.68 (5.50, 6.15) 0.92 5.80 (5.30, 6.50) 5.69 (5.20, 6.30) 0.003 HDL‐C, mmol/L 1.29 (1.17, 1.66) 1.44 (1.08, 1.60) 0.84 1.31 (1.16, 1.53) 1.40 (1.19, 1.65) <0.001 LDL‐C, mmol/L 2.67±1.02 2.60±0.84 0.79 2.89±0.82 2.78±0.86 0.03 TC, mmol/L 4.75 (4.12, 5.39) 4.34 (3.91, 5.29) 0.46 4.89 (4.29, 5.46) 4.84 (4.17, 5.45) 0.15 triglyceride, mmol/L 1.37 (0.92, 2.19) 1.36 (1.00, 1.70) 0.79 1.40 (1.04, 2.02) 1.23 (0.89, 1.68) <0.001 [96]Open in a new tab ACS indicates acute coronary syndrome; BMI, body mass index; CHD, coronary heart disease; DBP, diastolic blood pressure; FG, fasting glucose; HDL‐C, high‐density lipoprotein cholesterol; LDL‐C, low‐ density lipoprotein cholesterol; MET, metabolic equivalent; SBP, systolic blood pressure; and TC, total cholesterol. Selection of Candidate miRNAs in the Discovery Stage We screened 33 miRNAs to be significantly associated with incident ACS by paired t tests (P<0.05, fold change>1.3). Via comparing the 33 candidate miRNAs to the miRmine database, we removed 11 extremely low abundance miRNAs in plasma. Finally, we selected 18 candidate miRNAs that could be detected in pooled plasma samples for quantitative real‐time polymerase chain reaction validation (Table [97]S2). Hierarchical clustering showed that the cases and controls were grouped distinctly across the heat map of 18 candidate miRNAs (Figure [98]S3). Identification of Plasma miR‐4286 Associated With Incident ACS in the Validation Stage Among 14 miRNAs that were selected for further analyses after quality control, most of them were significantly correlated with each other (Spearman’s rank correlation coefficients range from −0.44 to 0.83; Figure [99]S4). After multivariate adjustment, we observed a strong positive association between miR‐4286 and incident ACS; participants with plasma miR‐4286 in the highest tertile of the distribution had a 1.80‐fold increased risk of ACS compared with participants with plasma miR‐4286 in the lowest tertile (OR, 1.80; 95% CI, 1.28–2.53; P‐trend<0.001, FDR=0.01). Corresponding ORs for an IQR increase in plasma miR‐4286 were 1.26 (95% CI, 1.07–1.48; Figure [100]1). In addition, compared with the lowest tertile, the greatest risk difference was 10.56% (95% CI, 3.46%–17.66%) for the highest tertile of plasma miR‐4286. Meanwhile, an IQR increase in plasma miR‐4286 was associated with an absolute risk increase in incident ACS of 4.43% (95% CI, 0.94%–7.91%; Figure [101]S5). Elastic net model including all 14 candidate miRNAs showed that only miR‐4286 had a nonzero coefficient (β=0.02), indicating a strong association between miR‐4286 and incident ACS. Figure 1. Associations of candidate miRNAs with incident ACS in the validation stage. Figure 1 [102]Open in a new tab Plasma miRNA levels were normalized to miR‐26b‐5p and were expressed as 2^−∆Ct. Adjusted OR (95% CI) for incident ACS was obtained using multivariable conditional logistic regression model with adjustment for age, body mass index, smoking status, drinking status, education levels, metabolic equivalent, diabetes mellitus, hypertension, family history of coronary heart disease, high‐density lipoprotein cholesterol, low‐density lipoprotein cholesterol, triglyceride, and use of lipid‐lowering medication. P‐trend was estimated by assigning the median value of miRNA to each tertile and using this as a continuous variable in the logistic regression model. ACS indicates acute coronary syndrome; FDR, false discovery rate; IQR, interquartile range; miRNAs, microRNAs; and OR, odds ratio. Of note, the significant result for miR‐4286 was also confirmed by the global mean normalization strategy, and the association remained significant after excluding cases with ACS diagnosed within 1 year (Tables [103]S3 and S4). Cubic spline regression analysis demonstrated a significant nonlinear association between plasma miR‐4286 and incident ACS (nonlinear P=0.002; Figure [104]S6). According to the stratified analyses, no significant interaction was observed between miR‐4286 and ACS‐related risk factors (Table [105]S5). The area under the receiver‐operating characteristic curves for prediction of incident ACS was 0.55 (95% CI, 0.52–0.59, P=0.002), with a sensitivity and specificity being 65.0% and 43.9% at 0.09 of plasma miR‐4286 as a cut‐off value. In comparison with FRS alone, the addition of miR‐4286 increased area under the receiver‐operating characteristic curves from 0.64 (95% CI, 0.61–0.68) to 0.66 (95% CI, 0.63–0.69; P difference=0.07), the net reclassification improvement (95% CI) was 3.2% (0.9%–5.4%, P=0.006), and the integrated discrimination improvement (95% CI) was 1.1% (0.5%–1.8%, P<0.001; Figure [106]S7). In addition, we observed that miR‐4286 was associated with higher risk of AMI and UA, with no significant heterogeneity observed for AMI and UA (P [heterogeneity]=0.49; Table [107]S6). It suggested that the associations of plasma miR‐4286 on ACS were consistent across different subtypes. Consequently, we selected miR‐4286 as target miRNA associated with incident ACS to further explore its influencing factors and its role in ACS progression. Associations of Cardiovascular Traits and Risk Factors With Plasma miR‐4286 We observed a higher miR‐4286 level in females and the significant correlations between lipid traits and plasma miR‐4286 (Tables [108]S7 and S8). Decreased HDL‐C level, increased triglyceride level, and elevated triglyceride/HDL‐C ratio were significantly associated with higher levels of plasma miR‐4286 after multivariable adjustment among all participants and incident ACS cases. According to the results among all participants, comparing the lowest tertiles, the highest tertiles of HDL‐C, triglyceride, and triglyceride/HDL‐C ratio were significantly associated with 18.08% decrease (95% CI, −28.21% to −6.51%; P‐trend=0.003, FDR=0.005), 23.89% increase (95% CI, 8.67%–41.24%; P‐trend=0.001, FDR=0.004), and 23.63% increase (95% CI, 8.37%–41.04%; P‐trend=0.001, FDR=0.004) in plasma miR‐4286, respectively. In addition, an IQR increase in HDL‐C, triglyceride, and triglyceride/HDL‐C ratio was significantly associated with an 11.05% decrease (95% CI, −16.65% to −5.08%), 11.04% increase (95% CI, 3.77%–18.83%), and 15.01% increase (95% CI, 7.16%–to 23.42%) in plasma miR‐4286, respectively (Figure [109]2). Moreover, we also found the significant associations of triglyceride, HDL‐C, and triglyceride/HDL‐C ratio with incident ACS (Table [110]S9). Figure 2. Associations of lipid traits with plasma miR‐4286 in the validation stage. Figure 2 [111]Open in a new tab Forest plots of A, B, and C separately show the associations of lipid traits with plasma miR‐4286 among all participants, incident ACS cases, and controls. Adjusted PD (95% CI) for plasma miR‐4286 was obtained using multivariable linear regression model with adjustment for age, sex, body mass index, smoking status, drinking status, education levels, metabolic equivalent, and use of lipid‐lowering medication. P‐trend was estimated by assigning the median value of lipid to each tertile and using this as a continuous variable in the linear regression model. ACS indicates acute coronary syndrome; FDR, false discovery rate; HDL‐C, high‐density lipoprotein cholesterol; IQR, interquartile range; LDL‐C, low‐density lipoprotein cholesterol; PD, percent difference; TC, total cholesterol; and TG, triglyceride. Causal Relationships Between Lipids and Plasma miR‐4286 The genome‐wide association study results of miR‐4286 are presented in Table [112]S10 and Figure [113]S8, and details of triglyceride and HDL‐C associated instrumental variables information are shown in Tables [114]S11 and S12. The F statistics of triglyceride, HDL‐C, and miR‐4286 were 76.50, 83.65, and 105.41, respectively (Table [115]S13). In the MR analyses, triglyceride was causally associated with an increased level of miR‐4286, and the adjusted β (95% CI) according to per 1‐unit increase was 0.27 (0.01–0.53; P=0.04) by conventional inverse variance–weighted test. Although the significance in MR‐Egger and weighted median tests was marginal, the significant effect was validated by the MR‐PRESSO test (0.27 [0.07–0.47]; P=0.01). The HDL‐C was negatively associated with miR‐4286 in MR‐PRESSO test, the adjusted β (95% CI) according to per 1‐unit increase was −0.20 (−0.39 to −0.01; P=0.04; Figure [116]3 and Figure [117]S9). Notably, there was no suggestion of a reverse causal effect of miR‐4286 on triglyceride or HDL‐C (Figures [118]S10 and S11). Figure 3. Causal effects of TG and HDL‐C on plasma miR‐4286. Figure 3 [119]Open in a new tab The causal effect (95% CI) of TG and HDL‐C on miR‐4286 was estimated by Mendelian randomization analysis, using 39 TG SNPs and 58 HDL‐C SNPs derived from Asian reports, respectively. Pleiotropy P value derived from the intercept of MR‐Egger test or MR‐PRESSO Global test, a small P value indicates an unbalanced pleiotropy. HDL‐C, high‐density lipoprotein cholesterol; IVW, inverse variance‐weighted; SNPs, single nucleotide polymorphisms; and TG, triglyceride. Causal Mediation Effects of Plasma miR‐4286 on the Associations of Lipids and Incident ACS The positive associations of triglyceride and triglyceride/HDL‐C ratio with incident ACS were partly mediated through an increasing level of plasma miR‐4286, and the mediation proportions (95% CI) were 5.5% (0.7%–17.0%; P [mediation]=0.02) and 6.6% (0.9%–19.5%; P [mediation]=0.02), respectively. Nevertheless, no significant mediation effect was detected for HDL‐C as exposure, presumably because of the weak association between HDL‐C and incident ACS (Figure [120]4). Figure 4. Mediation effects of plasma miR‐4286 on the associations of lipids with incident ACS. Figure 4 [121]Open in a new tab A, Path diagram: Total effect (c)=direct effect (c′)+indirect effect (ab). B, Total effects and the proportion of mediation effects were obtained by quasi‐Bayesian Monte Carlo simulation for 10 000 times in the R package “mediation.” Indirect effects were derived from multiple linear regression for the association of lipid with miRNA level adjusting for age, sex, body mass index, smoking status, drinking status, education levels, metabolic equivalent, and use of lipid‐lowering medication. Direct effects were calculated by logistic regression for the association of lipid with ACS additional adjusting for miRNA level, diabetes mellitus, hypertension, family history of coronary heart disease, and low‐density lipoprotein cholesterol. ACS indicates acute coronary syndrome; HDL‐C, high‐density lipoprotein cholesterol; miRNA, microRNA; and TG, triglyceride. Discussion Based on a 2‐stage nested case–control design, we found that plasma miR‐4286 was positively associated with incident ACS and causally related to triglyceride and HDL‐C. Importantly, an increase in plasma miR‐4286 also partially mediated the positive association of triglyceride and triglyceride/HDL‐C ratio with incident ACS. To summarize, our study identified a novel risk contributor of incident ACS and indicated its mediation potential in ACS development. According to the prospective findings from the Bruneck cohort and HUNT (Nord‐Trøndelag Health Study) cohort, circulating miRNAs associated with risk of AMI were nonoverlapping.[122] ^8 , [123]^9 Interestingly, in the present study, we also found none of the predicted biomarkers as abovementioned to be associated with incident ACS. In comparison with the Bruneck cohort including 47 incident AMI cases and 773 controls, we performed the nested case–control design that used age‐ and sex‐matched cases and controls, which enhanced statistical power and reduced heterogeneity.[124] ^8 Although the HUNT study also performed the nested case–control design, they have additionally matched on risk factors including BMI, total cholesterol, and HDL‐C. These could be potential explanations for the discrepancy, given that the miRNA identified in the present study was correlated with lipid levels.[125] ^9 In addition, the different platform and normalized strategy for miRNA analysis might lead to technical differences. In line with the HUNT study, we also confirmed our result by global mean normalization strategy. Of note, the association between miR‐4286 and incident ACS was moderate in this method. This weaker association could be explained by the fact that the closely associated miRNAs cluster might bias the average of cycle thresholds.[126] ^8 As for the predictive potential of plasma miR‐4286 for incident ACS, we observed a slight but significant improvement for net reclassification with the FRS plus miR‐4286, suggesting that miR‐4286 might be a novel ACS risk predictor independent of traditional risk factors. Nevertheless, given that adding plasma miR‐4286 to the FRS model could not significantly increase the area under the receiver‐operating characteristic curves for incident ACS prediction, the optimal combination of biomarkers including miRNAs or clinical traits in ACS risk prediction should be further explored. Participants with abnormal levels of plasma triglyceride and HDL‐C remain at a high risk of CVD even with normal low‐density lipoprotein cholesterol level.[127] ^30 Higher triglyceride and lower HDL‐C levels, as previously reported, were symbolic of dyslipidemia of insulin resistance (IR), which played an important role in the development of CVD.[128] ^31 Other studies have demonstrated that the triglyceride/HDL‐C ratio was positively associated with insulin sensitivity, and thus was considered a useful surrogate index to estimate IR.[129] ^32 Herein, we identified potentially causal associations of triglyceride and HDL‐C with plasma miR‐4286, as well as a positive association between triglyceride/HDL‐C ratio and miR‐4286. Accordingly, we suggested that an increase in plasma miR‐4286 might be caused by abnormal triglyceride and HDL‐C levels accompanied by impaired insulin sensitivity in the development of ACS. The impact of systematic chronic inflammation caused by impaired insulin sensitivity and abnormal lipid levels might account for the observation, which might stimulate the expression and release of miRNAs involved in the metabolic disturbance and ACS progression.[130] ^33 According to pathway enrichment analysis, our result suggested that the highly expressed target genes of miR‐4286 in the cardiovascular system were most enriched in the IR pathway (Figure [131]S12). Previous studies have demonstrated the crucial role of IR in endothelial dysfunction.[132] ^34 Indeed, proliferation, migration, and dysfunction of endothelial cells is a key point in the development of atherosclerosis and vulnerable plaque.[133] ^35 In this study, miR‐4286 could partially mediate the impact of triglyceride and triglyceride/HDL‐C on incident ACS. It seems possible that endothelium might be a potential target for the mediation effect of miR‐4286 on triglyceride‐associated ACS. Although the mechanism role of miR‐4286 in the IR‐related pathway in endothelial cells was lacking, functional study has found that miR‐4286 played an important part in cell proliferation and migration via regulating the PI3K/AKT pathway in the development of lung cancer.[134] ^36 From a clinical perspective, our findings inferred a possible novel target for miRNA‐based therapeutics.[135] ^37 Further experimental studies are required to investigate the functional role of miR‐4286 in endothelial cells, especially in insulin‐related pathways. As mentioned in a review,[136] ^38 the pool of circulating miRNAs is the sum of secreted and leaked miRNAs from different cells. MiRNAs derived from active secretion via microvesicles or binding with HDL or Argonaute2 remain stable in the extracellular environment. On the contrary, passively leaked miRNAs from injured cells were susceptible to degradation. According to the EVmiRNA database[137] ^39 —an up‐to‐date extracellular vesicle miRNA database including 462 sequencing data sets from 17 tissues—we found that miR‐4286 was undetectable in the secreted exosomes from endothelial cells or epithelial cells. However, it remains to be investigated whether miR‐4286 was passively released from endothelial cells or derived from other cardiovascular cells. In addition, based on the abovementioned database, miR‐4286 was also enriched in exosomes from B lymphoblastoid cell lines and lymph. It suggested that lymphocytes might be one of the sources of plasma miR‐4286. There are several strengths in the present study. First, our study was a prospective nested case–control design based on a well‐established large prospective cohort, which could provide a relatively strong basis for the associations between plasma miRNAs and incident ACS. Second, we adjusted for traditional risk factors and lipid‐lowering medication used in ACS in all association analyses to reduce potential confounding bias. Third, we further investigated the internal causal relationships of plasma miR‐4286 and lipids through MR analyses and explored the causal mediation effect of miR‐4286 on triglyceride‐incident ACS association. Several limitations in this study should also be acknowledged. First, compared with previous prospective studies, the follow‐up period of 3 years was relatively short.[138] ^8 , [139]^9 However, we have excluded incident ACS cases that were diagnosed within 6 months after blood collections, and our findings were largely unchanged after further excluding cases with a follow‐up period <1 year. Second, the sample size of the microarray exploration data was limited, which may contribute to a relatively high false‐positive rate. Although we used a 1:1 matched strategy to improve the statistical power and verified the association by correcting multiple testing in the validation stage, further long‐term prospective studies with a larger sample size are needed to replicate our results. Third, there is a lack of established endogenous control to normalize the circulating miRNA levels. In response, we chose 3 eligible miRNAs as candidate endogenous controls and used acknowledged NormFinder to estimate the optimal endogenous control and repeated a global mean normalized strategy to check the stability of our results. Fourth, expression patterns and function of miR‐4286 in tissues remain largely unclear. Although there is a lack of evidence for the underlying mechanism, we inferred a potential role of miR‐4286 in regulating insulin sensitivity via functional pathway analysis and provided a basis for further mechanism research. Finally, we could not determine the source of miR‐4286, given that dysregulation of miRNAs in different tissues might be affected by different factors. Further systematic exploration of the source of circulating miR‐4286, as well as its regulatory mechanism, needs to be conducted. In conclusion, the present study identified a positive association of plasma miR‐4286 with incident ACS in a Chinese population. We indicated the causal effect of triglyceride on plasma miR‐4286, and highlighted a possible mediation role of plasma miR‐4286 in triglyceride‐incident ACS association. Further prospective and functional studies are warranted to systemically elucidate the biological mechanisms behind these observations and translate these results into clinical practice for early and precise prevention and treatment. Sources of Funding This work was supported by the National Key Research and Development Program of China (2016YFC0900800), the National Natural Science Foundation of China (81930092), the Fundamental Research Funds for the Central Universities (2019kfyXMBZ015), the 111 Project, and the Program for Changjiang Scholars and Innovative Research Team in University. Disclosures None. Supporting information Supplemental Methods Tables S1–S13 Figures S1–S12 References [140]40 , [141]41 , [142]42 , [143]43 , [144]44 , [145]45 ,