Abstract
Coronary artery disease (CAD), a common cardiovascular disease, has become a vital cause of mortality worldwide. Genetic microRNA (miR) polymorphisms might contribute to CAD susceptibility. In this study, we selected miR-146a, miR-196a2, and miR-499 single nucleotide polymorphisms and conducted a case–control study. In total, 505 CAD cases and 1109 controls were recruited. We used SNPscan™ genotyping assay to obtain genotyping of miR rs2910164, rs11614913, and rs3746444 variants. We found that miR-196a2 rs11614913 T > C decreased the susceptibility of myocardial infarction (MI) (TC vs. TT: adjusted p = 0.007 and CC/TC vs. TT: adjusted p = 0.012). In female subgroup, our results indicated that miR-196a2 rs11614913 T > C variants might also decrease the susceptibility of CAD (TC vs. TT: adjusted p = 0.017 and TC/CC vs. TT: adjusted p = 0.015). In summary, these results suggest that miR-196a2 rs11614913 T > C locus decreases the susceptibility of CAD in female and MI subgroups. However, further studies are needed to validate the potential associations of miR-196a2 rs11614913 locus with CAD.
Introduction
Coronary artery disease (CAD) manifests in middle and old age populations worldwide (Newton et al., 2015). The etiology of CAD remains unclear. Traditional factors, such as smoking, drinking, obesity, overweight, hypertension, type 2 diabetes mellitus (T2DM), and hyperlipidemia, have contributed to the susceptibility and progression of CAD (Moon et al., 2015; Li et al., 2016; Albrektsen et al., 2017; Bokma et al., 2017; Lerman et al., 2017; Lim, 2017; Lui et al., 2017; Zhang et al., 2017). Nevertheless, some potential susceptibility factors promoting CAD development are needed to be studied. Recently, accumulating molecular epidemiological investigations have highlighted that hereditary factors are important for CAD, and the individual's single nucleotide polymorphisms (SNPs) are attracting more attention (Gui et al., 2014).
microRNAs (miRs), a group of noncoding RNAs, are ∼22 nucleotides in length (Noonan et al., 2010; Panarelli and Yantiss, 2011) and regulate gene expression at the post-translational level (Vasudevan et al., 2007). miRs are correlated with human diseases, including CAD, diabetes, cancer, and so on. Accumulating evidences have indicated the vital roles of miRs in regulating fundamental biological processes of CAD, such as proliferation, differentiation and stress resistance, and abnormal apoptosis (Klingelhoffer et al., 2016; Brunquell et al., 2017; Lu et al., 2017; Wang et al., 2017; Yu et al., 2017). In addition, several studies identified that miRs were implicated in the regulation of energy metabolism, oxidative stress, inflammation, and proliferation of smooth muscle cells (Maves et al., 1989; McDonald et al., 2015; Afzal et al., 2016; Jin et al., 2017), which are associated with the development of CAD.
SNPs located on pre-miR genes and certain miR-binding sites in target genes may influence the interaction of specific miR and mRNA target site and then lead to the abnormal expression of these genes (Sun et al., 2010). In view of the potential roles of miRs in regulating translation, genetic miR SNPs were thought to contribute to CAD susceptibility. Previous studies have focused on the correlation of some common miR SNPs (miR-146a rs2910164 C>G, miR-196a2 rs11614913 T>C, and miR-499 rs3746444 A>G) with CAD development. Recently, these important polymorphisms in miR were under investigation to explore the possible association to CAD, but the results remained conflicting. Thus, in this study, we tried to identify the correlation of miR SNPs with the susceptibility to CAD. We selected rs2910164, rs11614913, and rs3746444 SNPs and conducted a case–control study in Chinese population.
Materials and Methods
Study population and ethical approval
A total of 505 CAD patients who had presented to the Fujian Union Hospital from October 2014 to May 2016 were enrolled in the study. And 1109 controls were recruited at the same time. The criterion of included CAD cases was coronary stenosis ≥50%, which was confirmed based on coronary angiography (Kucukhuseyin et al., 2013). The major criteria of healthy controls were as follows: (1) without history of CAD, (2) without any symptom of CAD, and (3) without myocardial ischemia in electrocardiogram test (Kucukhuseyin et al., 2013).
The related risk factors and demography data, such as smoking, alcohol consumption, weight, height, body mass index (BMI), sex, age, fasting blood glucose (FBG), history of T2DM, hypertension, low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), triglyceride (TG), and total cholesterol (TC), were collected. The definition of smoker and drinker was presented in our previous studies (Tang et al., 2014c, 2016). The diagnosis of T2DM was based on WHO1999 criteria. The systolic pressure ≥140 mmHg and/or diastolic pressure ≥90 mmHg or history of antihypertensive therapy was defined as hypertension (Soderstrom et al., 2013). BMI ≥24 was defined as overweight or obese (Zhai et al., 2010). Hyperlipidemia was diagnosed according to Chinese Medical Association 2016 criteria (Luo et al., 2019).
The written informed consent was signed by each participant. The present study was approved by Ethics Committee of Fujian Medical University (Approval No. 2016KY006).
miR SNP selection
To analyze the correlation of miR SNPs with the risk of CAD, we chose rs2910164, rs11614913, and rs3746444 SNPs according to the previous publications; these polymorphisms were found to be associated with cancer (Chen et al., 2016; Yin et al., 2017a; Tang et al., 2019; Thakur et al., 2019), ischemic stroke (Huang et al., 2015b; Zhu et al., 2017a), T2DM (Alipoor et al., 2016; Kaidonis et al., 2016; Ciccacci et al., 2018), autoimmune diseases (Ciccacci et al., 2016; Park et al., 2016; Toraih et al., 2016; Assmann et al., 2017; Shaker et al., 2018), and CAD (Zhi et al., 2012; Chen et al., 2014, 2017b; Xiong et al., 2014; Li et al., 2015; Bastami et al., 2016; Fragoso et al., 2019).
Isolation of DNA and genotyping
Two-milliliter venous blood specimens were drawn into a test tube containing ethylene diamine tetraacetic acid. Genomic DNA was extracted from venous blood samples using the DNA Kit (Promega, Madison). We used SNPscan™ genotyping assay (Genesky Biotechnologies, Inc., Shanghai, China) to obtain genotyping of rs2910164, rs11614913, and rs3746444 polymorphisms. To ensure the accuracy of our results, 65 randomly selected DNA samples were tested repeatedly, and the results were unchanged.
Statistical analyses
The continuous variables (e.g., FBG, age, height, weight, BMI, TG, TC, LDL-C, and HDL-C) are presented as the mean ± standard deviation. Differences in age, sex, hypertension, T2DM, hyperlipidemia, smoking, drinking, BMI, and genotypes of selected polymorphisms were analyzed using χ
2 test. Student's t-test was harnessed to determine the potential difference between two groups. The relationship of rs2910164, rs11614913, and rs3746444 genotypes with CAD susceptibility was determined by odds ratios and their 95% confidence intervals. Hardy–Weinberg equilibrium (HWE) was calculated by an Internet HWE-analysis software (online at
Results
Population characteristics
The frequency distributions of drinking, BMI, age, sex, smoking, hypertension, T2DM, and hyperlipidemia information for the 505 CAD patients (64.57 ± 9.91 years) and 1109 healthy controls (64.75 ± 10.20 years) are presented in Table 1. The gender and age had no significant difference between two groups (p = 0.704 and 0.942, respectively). The differences of BMI, hypertension, T2DM, smoking, drinking, and hyperlipidemia between CAD patients and controls were significant (p < 0.05). Among the CAD patients, 147 were myocardial infarction (MI) cases and 358 were non-MI cases.
Distribution of Selected Demographic Variables and Risk Factors in Coronary Artery Disease Cases and Controls
Bold values are statistically significant (p < 0.05).
Two-sided χ 2 test and Student's t-test.
BMI, body mass index; CAD, coronary artery disease; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol; MI, myocardial infarction; TC, total cholesterol; TG, triglyceride.
Primary information of miR SNPs
The primary information of rs2910164, rs11614913, and rs3746444 SNPs is listed in Table 2. These SNPs were noncoding variants. Genotype distributions of three miR SNPs in controls were in accord with HWE (p > 0.05). All of the genotyping successful rates were 99.50% (Table 2). All data were presented in Supplementary Data.
Primary Information for miR- 146 a rs2910164 C>G, miR- 196 a 2 rs11614913 T>C, and miR- 499 rs3746444 A>G Polymorphisms
HWE, Hardy–Weinberg equilibrium; MAF, minor allele frequency; miR, microRNA; SNPs, single nucleotide polymorphisms.
Association between miR SNPs and CAD risk
The genotype frequencies of miR polymorphisms in different groups are summarized in Table 3. The frequencies of rs2910164 CC, CG, and GG genotypes were 38.80%, 49.20%, and 12.00% in CAD cases and 39.42%, 46.65%, and 13.92% in controls, respectively. Adjusting for BMI, hypertension, smoking, drinking, age, sex, T2DM, and hyperlipidemia, we identified that miR-146a rs2910164 was not associated with the risk of CAD (CG vs. CC: adjusted p = 0.709; GG vs. CC: adjusted p = 0.732; CG/GG vs. CC: adjusted p = 0.843; GG vs. CC/CG: adjusted p = 0.607; Table 4).
The Frequencies of miR- 146 a rs2910164 C>G, miR- 196 a 2 rs11614913 T>C, and miR- 499 rs3746444 A>G Polymorphisms in Coronary Artery Disease Patients and Controls
Overall and Stratified Analyses of miR- 146 a rs2910164 C>G, miR- 196 a2 rs11614913 T>C, and miR- 499 rs3746444 A>G Polymorphisms with Coronary Artery Disease
Bold values are statistically significant (p < 0.05).
Adjusted for age, sex, smoking, drinking, BMI, hypertension, T2DM, and hyperlipidemia.
95% CI, 95% confidence interval; OR, odds ratio; T2DM, type 2 diabetes mellitus.
The frequencies of rs3746444 AA, AG, and GG variants were 71.40%, 27.60%, and 1.00% in CAD group and 72.33%, 25.23%, and 2.44% in controls, respectively. When we adjusted BMI, hypertension, smoking, drinking, age, sex, T2DM, and hyperlipidemia, the null association of miR-499 rs3746444 A>G locus with CAD risk was also found (AG vs. AA: adjusted p = 0.685; GG vs. AA: adjusted p = 0.157; AG/GG vs. AA: adjusted p = 0.944; GG vs. AA/AG: adjusted p = 0.149; Table 4).
The frequencies of rs11614913 TT, TC, and CC genotypes were 32.60%, 49.40%, and 18.00% in CAD group and 30.29%, 52.08%, and 17.63% in controls, respectively. When we adjusted BMI, hypertension, smoking, drinking, age, sex, T2DM, and hyperlipidemia, the results suggested that rs11614913 genotype was not different between CAD cases and controls (TC vs. TT: adjusted p = 0.239; CC vs. TT: adjusted p = 0.903; TC/CC vs. TT: adjusted p = 0.335; CC vs. TT/TC: adjusted p = 0.620; Table 4).
Association of miR SNPs with CAD risk in different CAD types
The genotype frequencies of rs2910164, rs11614913, and rs3746444 polymorphisms in different stratification groups are summarized in Table 3. To determine whether the role of rs2910164, rs11614913, and rs3746444 variants was influenced by CAD type, we performed a stratified analysis. We identified that rs11614913 SNP was associated with a decreased susceptibility of MI (TC vs. TT: adjusted p = 0.007 and CC/TC vs. TT: adjusted p = 0.012). However, we found that rs2910164 and rs3746444 loci were not associated with the susceptibility of CAD (Table 4).
Association of miR SNPs with CAD risk in a stratification analysis
Table 5 summarized rs11614913 variant frequencies in different drinking, BMI, hypertension, age, smoking, T2DM, sex, and hyperlipidemia subgroups. In female subgroup, we identified that rs11614913 SNP decreased the risk of CAD (TC vs. TT: adjusted p = 0.017 and TC/CC vs. TT: adjusted p = 0.015). However, no difference was shown in distribution of rs2910164 and rs3746444 variants among CAD patients and controls in other subgroups (Tables 6 and 7).
Stratified Analyses Between miR- 196 a2 rs11614913 T > C Polymorphism and Coronary Artery Disease Risk by Sex, Age, Body Mass Index, Smoking Status, and Alcohol Consumption
Bold values are statistically significant (p < 0.05).
For miR-196a2 rs11614913 T>C, the genotyping was successful in 500 (99.01%) CAD cases and 1106 (99.73%) controls.
Adjusted for age, sex, smoking, drinking, BMI, hypertension, T2DM, and hyperlipidemia (besides stratified factors accordingly) in a multiple logistic regression model.
Stratified Analyses Between miR- 146 a rs2910164 C > G Polymorphism and Coronary Artery Disease Risk by Sex, Age, Body Mass Index, Smoking Status, and Alcohol Consumption
For miR-146a rs2910164 C>G, the genotyping was successful in 500 (99.01%) CAD cases and 1106 (99.73%) controls.
Adjusted for age, sex, smoking, drinking, BMI, hypertension, T2DM, and hyperlipidemia (besides stratified factors accordingly) in a multiple logistic regression model.
Stratified Analyses Between miR- 499 rs3746444 A>G Polymorphism and Coronary Artery Disease Risk by Sex, Age, Body Mass Index, Smoking Status, and Alcohol Consumption
For miR-499 rs3746444 A>G, the genotyping was successful in 500 (99.01%) CAD cases and 1106 (99.73%) controls.
Adjusted for age, sex, smoking, drinking, BMI, hypertension, T2DM, and hyperlipidemia (besides stratified factors accordingly) in a multiple logistic regression model.
Combination analysis of three miR-polymorphisms (rs11614913, rs2910164, and rs3746444)
Table 8 showed the combined analysis of rs11614913, rs2910164, and rs3746444 polymorphisms. Three different combinations, rs11614913/rs2910164, rs11614913/rs3746444, and rs2910164/rs3746444, were made to explore the interaction of different genotypes of miR-SNPs and their potential roles on the development of CAD. However, we did not find any association between combined genotypes and the risk of CAD.
Combination Analysis of microRNA Polymorphisms (rs11614913, rs2910164, and rs3746444) in Coronary Artery Disease Cases and Controls
Discussion
Several case–control studies have regarded miR rs2910164, rs11614913, and rs3746444 polymorphisms as promising candidates for CAD. However, the association between these SNPs and CAD was conflicting. The present study in Eastern Chinese Han populations failed to confirm this relationship in overall comparison. However, we found that rs11614913 locus might decrease the susceptibility of CAD in MI and female subgroups. To our knowledge, this is the first study suggesting the protective role of miR-196a2 rs11614913 T > C polymorphism to MI.
Luthra et al. (2008) indicated that miR-196a2 could regulate annexin A1 (ANXA1). ANXA1 is important for controlling inflammation. An increased expression of ANXA1 might be a protective factor for CAD (Bergstrom et al., 2017). Thus, miR-196a2 might be associated with the development of CAD. Mir-196a2 rs11614913 T → C variant locates in the 3p strand of mature miR regions, which could influence the pre-miR maturation and the interaction between 3p mature miRs and target mRNAs (Landgraf et al., 2007). Sung et al. (2016) suggested that rs11614913 T > C polymorphism was a protective factor for CAD. However, the others reported null association (Zhi et al., 2012; Chen et al., 2014; Xiong et al., 2014; Huang et al., 2015a).
Recently, a quantitative assessment of the relationship between miR-196a2 rs11614913 T > C polymorphism and susceptibility of CAD in Asians indicated that this SNP may contribute to a decreased susceptibility of CAD (Wang et al., 2018). In this study, our findings were similar to these reports. Some functional studies reported that rs11614913 T → C variant might be associated with a higher expression miR-196a2 in tumor tissues (Zhao et al., 2016; Yin et al., 2017b). In glioma tissue, a previous study found that the rs11614913 T > C SNP could influence the expression of target gene HOXC8 (Sibin et al., 2017). In liver cirrhosis patients with hepatopulmonary syndrome, compared to the patients carrying rs11614913 TT genotype, the expression of miR-196a increased in individuals with rs11614913CT/CC genotype (Chen et al., 2017a).
In addition, Xu et al. (2009) suggested that miR-196a2 rs11614913 CC could promote the expression of mature miR-196a in cardiac tissue, and the higher expression of miR-196a might lead to a decreased mRNA target of homeobox B8. Thus, miR-196a2 rs11614913 T > C polymorphism may influence the risk of CAD by increasing the expression of miR-196a2 and decreasing mRNA target of homeobox B8. It is noticeable, however, that the relatively few studies were performed on the relationship of rs11614913 with CAD risk, and the sample size of these studies was moderate. In addition, we only found the association between rs11614913 with decreased CAD risk in MI and female subgroups. It could be explained that the association of rs11614913 T > C SNP with the risk of CAD might be influenced by environment factor and subtype of disease.
Considering a fact that most of the SNPs make a small-to-moderate contribution to the development of CAD, the potential role of the rs11614913 T > C locus in miR-196a2 may be diluted or masked by interaction of gene–environment factors. In the future, more studies with detailed environment factors are needed to support our findings.
Nevertheless, some potential limitations should be considered. First, the CAD cases and controls were both recruited from local hospital, and these participants could not be well representative of the Chinese, and bias might have occurred. Second, the SNPs studied in the present study were selected based on literatures and could not obtain an extensive view of the heritable variations of miRs. Further case–control studies are needed to explore the mechanism of CAD development by a functional study. Third, for the moderate sample size in MI and female subgroup, the power might be insufficient in these genetic models. Finally, we did not collect detailed information on cardiovascular event, which restricted further assess of the roles of the miR rs2910164, rs11614913, and rs3746444 polymorphisms in CAD progression and prognosis.
To sum up, this study highlights that rs11614913 locus might decrease the susceptibility of CAD in MI and female subgroups. However, further studies are needed to confirm these potential associations between rs11614913 polymorphism and CAD development.
Footnotes
Acknowledgments
The authors appreciate all subjects who participated in this study. The authors thank Dr. Yan Liu (Genesky Biotechnologies, Inc., Shanghai, China) for technical support.
Disclosure Statement
No competing financial interests exist.
Funding Information
The project was supported by Research Foundation for Senior Talents of Jiangsu University (Grant no. 16JDG066) and Zhenjiang Social Development Science and Technology Project (SH2014087).
Supplementary Material
Supplementary Data
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
