Skip to main content

Analysis on the polymorphisms of site RS4977574, and RS1333045 in region 9p21 and the susceptibility of coronary heart disease in Chinese population

Abstract

Background

Rs4977574 (A > G) and Rs1333045 (C > T) are both single nucleotide polymorphisms (SNPs) related with coronary artery disease, locating on chromosome 9p21.3. The study aimed to identify the correlation between rs4977574 and rs1333045 polymorphism genotypes and coronary heart disease (CHD) in a Chinese population.

Methods

Blood samples were collected from 855 subjects. A case-control study was used in this experiment, and 598 cases in the CHD group and 257 subjects in the control group were enrolled. Genotyping was identified by the Agena MassARRAY system. Statistical analysis was conducted by SPSS (Ver 16.0) and plink (Ver. 1.07, Shaun Purcell). Haplotype analysis was performed using Haploview software.

Results

Association analysis by plink indicated a significant difference in the allele distribution for single nucleotide polymorphisms between cases and controls (rs4977574 P = 0.003, rs1333045 P = 0.035). Fisher’s exact test by plink proved that allele G may be associated with a higher risk of CHD (P = 0.003, odds ratio (OR= 1.371) and the T allele was likely to reduce the risk of coronary events (P = 0.035, OR = 0.798). The serum levels of apolipoprotein A (ApoA) were higher in subjects with the AG + AA genotype of rs4977574 compared to those with the GG genotype (P = 0.028). In the dominant model of rs1333045, the levels of ApoA were higher and LDL levels were lower in the TC + TT genotype than in the CC genotype.

Conclusions

The present study examined the association between the 9p21 chromosome rs4977574 and rs1333045 polymorphism genotypes and CHD in a population of Chinese patients. The G allele of rs4977574 and the C allele of rs1333045 are the susceptibility sites of CHD.

Peer Review reports

Background

Coronary heart disease (CHD) is coronary atherosclerotic heart disease. CHD often occurs due to coronary atherosclerosis, and the rupture of the atherosclerotic plaque contributes to thrombosis, resulting in myocardial infarction, myocardial remodeling, heart failure, and even death. CHD is the top cause of human death in the United States [1] and Asian countries such as China [2] and Japan [3]. CHD has been a hot topic in medical research because of its high incidence and high mortality. New research shows that CHD is a complex disease caused by a combination of genetic and environmental factors [4]. In addition to traditional risk factors such as age, smoking, high blood pressure and abnormal blood lipids, the combination of disease susceptibility genes and environmental factors affects the development of CHD. Identifying the susceptibility gene for coronary heart disease can help in early diagnosis and treatment. However, the exact mechanism of the disease is still unclear. In recent years, genetic predisposition has been thought to be closely associated with CHD and has been widely studied [5, 6].

Based on genome-wide association studies (GWAS), one of the important loci associated with CHD is the 9p21 locus, which has been introduced by the deCODE database. CDKN2B-AS1 of antisense noncoding RNA in the INK4 locus (ANIRL) is an antisense long noncoding RNA (lncRNA) mapped on the 9p21 locus [7] consisting of 19 exons, spanning 126.3 kb in the genome. The gene contributes to epigenetic inhibition of transcription by regulating Polycomb proteins [8]. According to the published literature, this gene is involved in developing cardiovascular diseases, diabetes, Alzheimer’s disease, and several types of cancer [9,10,11]. Regulation of cardiac CDKN2B-AS1 expression has been found to play a pivotal role in the development of CHD by altering the dynamics of vascular cell proliferation [12]. Rs4977574 is a nonprotein coding SNP (A > G) that is located in proximity to the cyclin-dependent kinase inhibitor 2A and B genes on chromosome 9p21.3. Rs4977574 has been recently found to be associated with the early onset of CHD, which is differentially expressed in a variety of tissues, such as vascular endothelial cells and smooth coronary muscle cells [13,14,15,16], and rs4977574 is characterized by a guanine nucleotide (G) instead of an adenine nucleotide (A) [17]. rs1333045 (C > T) is also an artery disease susceptibility SNP located in a conserved region in CDKN2B-AS1 that has been shown to have enhancer activity in a reporter gene experiment [13]. However, the underlying mechanisms in the regulation of rs4977574 and rs1333045 to CHD are unknown, and there is a lack of genetic evidence describing the pathogenesis in greater detail. Giving its limitations, coronary arteriography remains the gold standard for documenting the extent and severity of CHD.

Therefore, to provide a more detailed description of the correlation between Rs4977574 (A > G) and Rs1333045 (C > T) and CHD in a Chinese population, we conduct a case and control study by analysing the datasets collected from 855 subjects. It demonstrates that the G allele of rs4977574 and the C allele of rs1333045 are the susceptibility sites of CHD, laying a foundation for further study about the mechanisms in the regulation of rs4977574 and rs1333045 to CHD.

Methods

Study subjects

All selected subjects signed informed consent, and the study was approved by the ethics committee of the First Affiliated Hospital of Nanjing Medical University and the Friendship Hospital of Ili Kazakh Autonomous Prefecture in China. There were 855 cases, including 629 males and 226 females, with an average age of 60 years. Patients with spastic angina pectoris, infections within the previous 2 weeks, heart failure, adrenal dysfunction, and thyroid dysfunction were excluded from the study.

In the CHD group, we selected a total of 598 patients (465 males and 133 females) aged 61 (53–70) years who were clinically diagnosed with CHD in the Friendship Hospital of Ili Kazakh Autonomous Prefecture in China from March 1, 2010 to April 31, 2015. The patients were confirmed by coronary angiography, which was performed by at least two experienced doctors simultaneously. In addition, 257 subjects (164 males and 93 females) aged 58 (49.5–66) years with negative angiography were selected as the control group. Coronary arteries were cannulated using either the Judkins technique [18] or a radial artery approach with 6F catheters. CHD subjects were defined as at least one major epicardial vessel with > 50% stenosis; control subjects were defined as all of the major epicardial vessels with < 50% stenosis [19]. All subjects were unrelated Chinese to ensure the consistency of genetic background.

Laboratory investigation

In this study, four milliliters of peripheral blood was drawn from the subjects’ veins after 12 h of fasting to perform biochemical assays on the second day of hospitalization. Each assay index was the reported clinical parameter associated with coronary heart disease [20,21,22,23,24,25,26]. Total cholesterol (TC, mmol/L), triglyceride (TG, mmol/L), creatinephosphokinase isoenzyme (CKMB, U/L), fasting high-density lipoprotein cholesterol (HDL-C, mmol/L), fasting low-density lipoprotein cholesterol (LDL-C, mmol/L), apolipoprotein A (ApoA, g/L), apolipoprotein B (ApoB, g/L), glucose (Glu, mmol/L), blood urea nitrogen (Bun, mmol/L) and serum creatinine (Cre, μmol/L) were detected by the enzymatic method with an automated autoanalyzer (AU 2700 Olympus, 1st Chemical Ltd., Japan). Blood pressure was measured in the bare right arm with the participant seated, and the results were averaged over three measurements.

Smoking and drinking

For smoking, patients were classified as either nonsmokers (never smoking) or smokers (smoking now or in the past). Patients who smoked one cigarette a day were classified as current smokers. For drinking, patients were classified as nondrinkers (never having a drink) or drinkers (drinking now or in the past). Patients who drank at least 50 g of alcohol a week were considered drinkers.

Hypertension

Patients were defined as having hypertension if they had been previously diagnosed with hypertension and took antihypertensive drugs or if their systolic blood pressure was 140 mmHg or their diastolic blood pressure was 90 mmHg.

Scoring of coronary angiogram

Coronary angiography results were scored according to Gensini’s scoring system. In Gensini’s scoring system [27], the degree of coronary stenosis in the range of 0–25% was scored as 1 point, stenosis in the range of 25–50% was scored as 2 points, 50–75% was scored as 4 points, 75–90% was scored as 8 points, 90–99% was scored as 16 points, and completely occluded was scored as 32 points. Every stenosed segment was divided into 0.5 to 5 according to the importance of the segment’s function. According to Gensini’s standard, each major vessel and segment were graded, and the total score was calculated to show the severity of coronary lesion.

DNA extraction and determination of SNP genotypes

The experimental process of MassARRAY SNP typing is shown in Fig. 1. Gene polymorphisms were identified on the Agena MassARRAY system (Agena/Sequenom Inc., San Diego, CA, USA) according to the manufacturer’s instructions. The literature and AssayDesigner3.1 software (Sequenom Inc., San Diego, CA, USA) were both used to design polymerase chain reaction (PCR) and single base extension primers. The primers were diluted for backup analysis by a professional biotech company.

Fig. 1
figure1

The experimental flow graph of MassARRAY SNP typing

The DNA in blood samples, tissues, cells and saliva was extracted using the finished product kit [Axygene Biotechnology (Hangzhou) Limited, Hangzhou City, China]. All DNA quality checks were qualified by using the NanoDrop2000 instrument (Thermo, Wilmington, DE, USA) for OD value detection and 1.25% agarose gel electrophoresis detection. All samples were transferred to a 96-well plate and stored for reserve at − 20 °C. DNA samples were amplified using standard PCR. The PCR master mix was configured in a 1.5 ml EPtube and oscillated at low speed centrifugation. Eight or 12 channel shifters were chosen, and 4 μL of PCR master mix were added to each well in 384 well plate, along with 1 μL of template DNA (20 ng/μL). The 384-well plate was carefully sealed and pressed over each well to prevent evaporation during the PCR procedure. The plate was centrifuged at 1000 RPM for 1 min. The reaction environment and procedures were set as 94 °C, 5 min; 94 °C, 20 s; 56 °C, 30 s; 72 °C, 1 min; 72 °C, 3 min. Forty-five amplification cycles were performed, and the reactants were stored at 4 °C. The PCR products were treated with shrimp alkaline phosphatase (SAP) to remove the free dNTPs in the system. Five microliters of PCR products and 2 μL of SAP mix were mixed in 384-well plates. The SAP reaction program was set to 37 °C for 20 min and heated to 85 °C for 5 min after centrifugation. The products were stored at 4 °C. Two microliters of extend mix and 7 μL of PCR + SAP reaction reagent were mixed in each well of 384-well plates. The products were transferred into 384-well spectroCHIP bioarray by MassARRAY Nanodispenser RS1000 sampler (Agena, Inc.). The gene chip was analyzed by mass spectrometry (MALDI-TOF-MS). The original data and genotypes were obtained using MassARRAY TYPER 4.0 software (Agena, Inc).

Statistical analysis

All statistical analyses were performed with the Statistics Package for Social Sciences statistical software (Ver. 16.0, Chicago, USA) and plink (Ver. 1.07, Shaun Purcell). Quantitative variables are presented as the mean ± standard deviation, while qualitative variables are expressed as percentages. Skewed data are presented as the median (interquartile range). Smoking status, drinking status, hypertension, gender and the SNP site were examined by the chi-square test. The serum level of ApoB was examined by independent samples t test, and the remaining baseline characteristics were examined by the Mann-Whitney test for the demographics of the study population. To detect clinical parameters distributed in different genotypes, the serum level of ApoB was examined by one-way ANOVA, and the remaining baseline characteristics were examined by the Kruskal-Wallis test. Fisher’s exact test was used for the distribution of alleles between the two groups. Logistic regression analysis was used to assess the significance of the association between SNP rs4977574, rs1333045 and CHD. Comparison of the genotype and environmental factors between cases and controls was determined by multivariable logistic regression analysis. The odds ratio (OR) and 95% confidence intervals (CI) were also calculated. The Hardy-Weinberg law tests the equilibrium of allele frequency and genotype frequency in two groups. Pairwise linkage disequilibrium (LD) between sites was assessed by calculating D’ and squared correlation (r2) using Haploview software. P < 0.05 was considered statistically significant between the two groups.

The effects of the combination of gene and environment include not only the effects of both genes but also the superposition of gene and environmental effects and the multiplicative effect of gene and environmental effects. Through different models, we can distinguish the existence and size of interactions between two factors based on different models. Based on the additive model, there are several indicators for calculating interaction: (1) the interaction index (the synergy index S, S): when S = 1, no interaction; when S > 1, there is a positive interaction between the two factors; when S < 1, there is a negative interaction between the two factors. (2) The attribution ratio of interaction (attributable proportion of interaction AP): it shows how much of the total effect is attributable to two-factor interactions. (3) The attribution ratio of pure interaction: AP* indicates the proportion of interaction between two factors, which is caused by two factors of gene and environment. (4) The interaction excess relative risk degree (relative excess risk of interaction, RERI): this indicator indicates the difference between the combined effect of two factors and the sum of their individual effects and indicates the relationship between interaction and the effect of factors other than the two factors. If the unknown factor plays a great role, the interaction of the study becomes very minor and meaningless. Here, RERI is the two-factor interaction value based on the additive model.

Results

Hardy-Weinberg equilibrium

Hardy-Weinberg equilibrium (HWE) analysis revealed that the two study groups showed a remarkable genetic disequilibrium for both genotypes of rs4977574 and rs1333045 polymorphisms in this study, and the two groups were representative of the group.

Characteristics of the subjects

All 855 Chinese participants were genotyped for the rs4977574 and rs1333045 SNPs and included in the study. In this study, the CHD group included 465 males and 133 females, and the control group included 164 males and 93 females. The clinical characteristics of the participants are shown in Table 1. There were significant differences in age (P<0.001) and gender (P<0.001) between the CHD and control groups; older or male subjects are presented in the CHD group. Smokers (P = 0.031) are more susceptible to CHD, and the incidence of hypertension (P = 0.008) in the CHD group was significantly increased. The serum levels of CKMB (P = 0.005) in the CHD group were higher than in the control group, while the levels of HDL (P = 0.015) were much lower. In addition, the content of Glu (P<0.001) and Cre (P = 0.003) in CHDs was significantly higher than that in the controls. These variables are likely to be associated with CHD. The levels of TC and TG revealed no strong differences between CHDs and controls, which may be related to the use of lipid-lowering drugs in hospitalized patients.

Table 1 Demographics of the study population

Clinical parameter distribution in different genotypes

According to Table 1, the subjects were divided into three groups according to different genotypes: GG represents the homozygote of minor alleles, AG represents the heterozygote, and AA represents the homozygote of major alleles in rs4977574. TT represents the homozygote of minor alleles, TC represents the heterozygote, and CC represents the homozygote of major alleles in rs1333045. The three groups of genotypes were compared regarding clinical parameters, and the results are summarized in Table 2. Accordingly, no significant difference was observed between the frequencies of rs4977574 and rs1333045 in the distribution of age, gender, smoking status, drinking status, hypertension, CKMB, TC, TG, HDL, LDL, ApoA, ApoB, Glu, Bun and Cre (P > 0.05). However, the Gensini score did differ significantly among the rs4977574 genotypes (P = 0.017). This indicates that the G allele was associated with the severity of CHD. The allele frequency in subjects is all consistent with HWE (P = 0.596 for rs4977574, P = 0.566 for rs1333045).

Table 2 Clinical parameters distributed in different genotypes

Serological biomarkers in the dominant model and recessive model

Table 3 shows the serological biomarkers in the dominant model (rs4977574, AA vs. AG + GG; rs1333045, CC vs. TC + TT) and recessive model (rs4977574, GG vs. AA+AG; rs1333045, TT vs. TC + CC). Table 3 shows that the serum levels of ApoA in the GG genotype were significantly lower than those in the AA+AG genotype (P = 0.028). In the dominant model of rs1333045, the serum levels of LDL in the CC genotype were higher than those in the TC + TT genotype (P = 0.032). In contrast, the serum levels of ApoA in the CC genotype were significantly lower than those in the TC + TT genotype (P = 0.049). In addition, there were no differences between the two groups in other indicators. Moreover, significant difference of the serological biomarkers was observed neither in overdominant models of RS4977574 or RS1333045 (details seen in Additional file 2: Table S2).

Table 3 Comparisons of serological biomarkers in the dominant model and recessive model

The analysis of the rs4977574 and rs1333045 gene polymorphism with CHD

Fisher’s exact test by plink proved that allele G may be associated with a higher risk of CHD (P = 0.003, OR = 1.371, 95% CI = 1.113–1.689) and allele T was likely to reduce the incidence of coronary events (P = 0.035, OR = 0.798, 95% CI = 0.649–0.982) (Table 4). The chi-square test showed consistent results, allele G of rs4977574 (P = 0.003, OR = 1.371, 95% CI = 1.113–1.689) and allele T of rs1333045 (P = 0.033, OR = 0.798, 95% CI = 0.649–0.982) (Table 5).

Table 4 Fisher’s exact test of alleles
Table 5 The chi-square test of alleles

The genotypic distribution and allele frequencies of rs4977574 and rs1333045 in CHDs and controls are shown in Table 6. The rs4977574-A allele frequencies were 0.498 in CHDs and 0.576 in controls, and rs4977574-G allele frequencies were 0.502 in CHDs and 0.424 in controls. The rs1333049-C allele frequency was 0.525 in the CHD group and 0.469 in the control group, and rs1333045-T allele frequencies were 0.475 in CHDs and 0.531 in controls.

Table 6 Logistic analysis of the association between SNPs rs4977574 and rs1333045 and CHD risk

The association between the genotypes of SNP rs4977574 and CHD risk was investigated by logistic regression analysis. For rs4977574, the CHD group included 149 individuals with the AA allele, 297 individuals with the AG allele, and 152 individuals with the GG allele. In the control group, 87, 122, and 48 subjects had AA, AG, and GG alleles, respectively. The study results revealed a significant difference between the CHDs and controls in the rs4977574 polymorphism (P = 0.012). The frequency of genotype GG (P = 0.042, OR = 1.421, 95% CI = 1.014–1.993) and AG (P = 0.004, OR = 1.849, 95% CI = 1.217–2.810) were increased in the CHD group compared with the control group. In comparison to genotype CC, genotype TC significantly reduced the incidence of CHD (P = 0.035, OR = 0.641 and 95% CI = 0.424–0.970). As reported in Table 6, the subjects with AG + GG (P = 0.008, OR = 1.542, 95% CI = 1.122–2.119) genotypes had a higher risk of CHD than those with AA genotype. Corresponding to this, in a recessive model (GG vs. (AA+AG)), rs4977574 may decrease the incidence of CHD (P = 0.033, OR = 0.674 and 95% CI = 0.468–0.970). Moreover, rs1333045 was more likely to be a protective variant of CHD under the dominant model (CC vs. TC + TT, P = 0.048, AOR = 0.700, 95% CI = 0.492–0.997) when the P value was adjusted for age and gender. The gene polymorphisms of two sites in the CHD group and control group were in Hardy-Weinberg equilibrium.

The relationship between the clinical characteristics and CHD

For further study, we used the receiver operating characteristic (ROC) to describe the clinical indicators associated with CHD. The five indicators of significance are age, CKMB, HDL, Glu and Cre. The area under curve (AUC) is 0.583 for age (P<0.001, 95% CI = 0.541–0.624); 0.560 for CKMB (P = 0.005, 95% CI = 0.520–0.601); 0.552 for HDL (P = 0.015, 95% CI = 0.511–0.593), 0.617 for Glu (P<0.001, 95% CI = 0.578–0.657); and 0.563 for Cre (P = 0.003, 95% CI = 0.521–0.605) (as shown in Table 7).

Table 7 Receiver operating characteristic curve analyses for predicting CHD prevalence

The CHD prevalence and gene-environment interactions on the risk of CHD

In the logistic regression of interactions between CHD, genotypes and variables, crossover analysis was used to reveal the relationship between CHD prevalence and gene-environment interactions on the risk of CHD (details seen in Additional file 1: Table S1).

In the dominant model of rs4977574, the AG + GG genotype showed an increased prevalence of CHD in participants who were ≥ 59.5 years of age (P < 0.001, OR = 2.304, 95% CI = 1.448–3.666) compared to the reference group who had AA genotype who were < 59.5. Females had a lower risk of CHD than males (AA, P = 0.004, OR = 0.415, 95% CI = 0.229–0.753). In addition, the AA genotype without smoking or drinking of the reference group, the AA genotype with smoking (P = 0.029, OR = 1.831, 95% CI = 1.063–3.154), the AG + GG genotype with smoking (P < 0.001, OR = 2.249, 95% CI = 1.448–3.493), drinking (P = 0.034, OR = 1.803, 95% CI = 1.044–3.113) or without smoking (P = 0.005, OR = 1.824, 95% CI = 1.201–2.77), and without drinking (P = 0.004, OR = 1.661, 95% CI = 1.178–2.341) were more likely to have CHD. Patients who do not have hypertension (AG + GG, P = 0.004, OR = 1.914, 95% CI = 1.230–2.977), and who have hypertension (AA, P = 0.010, OR = 2.028, 95% CI = 1.184–3.473; AG + GG, P < 0.001, OR = 2.487, 95% CI = 1.591–3.887) have an increase in CHD prevalence. For the subjects had CKMB< 25.9 (AG + GG, P = 0.016, OR = 1.515, 95% CI = 1.080–2.123), CKMB≥25.9 (AA, P = 0.012, OR = 3.262, 95% CI = 1.296–8.212; AG + GG, P < 0.001, OR = 6.075, 95% CI = 2.982–12.375) existed as a risk factor compared with the AA genotype where CKMB< 25.9. Similarly, the AG + GG genotype with Glu < 5.445 (P = 0.004, OR = 1.750, 95% CI = 1.194–2.564) and the AA or AG + GG genotype with Glu ≥ 5.445 (AA, P < 0.001, OR = 3.551, 95% CI = 1.910–6.603; AG + GG, P < 0.001, OR = 4.399, 95% CI = 2.721–7.113) may have a greater chance of having CHD compared to the reference group (Glu < 5.445, AA). Moreover, the AG + GG genotype increased the risk of CHD in subjects with Cre ≥ 73.2 (P < 0.001, OR = 2.453, 95% CI = 1.565–3.844).

In the recessive model of rs4977574, females had a lower risk of CHD (AG + AA, P < 0.001, OR = 0.333, 95% CI = 0.201–0.551). No drinking also contributed to health (AG + AA, P = 0.019, OR = 0.625, 95% CI = 0.422–0.925). The AG + AA genotype exists as a positive factor for health in subjects with CKMB< 25.9 (AG + AA, P = 0.031, OR = 0.654, 95% CI = 0.444–0.962), Glu < 5.445 (AG + AA, P = 0.003, OR = 0.521, 95% CI = 0.338–0.802), and HDL > 1.325 (AG + AA, P = 0.001, OR = 0.405, 95% CI = 0.236–0.697; GG, P = 0.05, OR = 0.516, 95% CI = 0.267–1.000). However, CKMB≥25.9 (AG + AA, P = 0.005, OR = 2.675, 95% CI = 1.336–5.356), Glu ≥ 5.445 (AG + AA, P = 0.034, OR = 1.741, 95% CI = 1.041–2.911) or Cre ≥ 73.2 (GG, P = 0.038, OR = 2.122, 95% CI = 1.041–4.328) increased the prevalence of CHD.

For rs1333045, in the dominant model, female (CC, P = 0.037, OR = 0.501, 95% CI = 0.262–0.958; TC + TT, P < 0.001, OR = 0.365, 95% CI = 0.226–0.590) or no drinking (TC + TT, P = 0.042, OR = 0.680, 95% CI = 0.469–0.985) possibly reduced the risk of CHD. There are other protective factors, TC + TT genotype with HDL < 1.325 (CC, P = 0.023, OR = 0.530, 95% CI = 0.307–0.916), HDL ≥ 1.325 (CC, P = 0.005, OR = 0.405, 95% CI = 0.217–0.757; TC + TT, P = 0.001, OR = 0.394, 95% CI = 0.232–0.670) or Glu < 5.445 (TC + TT, P = 0.006, OR = 0.559, 95% CI = 0.371–0.843). In contrast, all other interactions increased the risk of CHD. Subjects with age ≥ 59.5 (CC, P = 0.006, OR = 2.400, 95% CI = 1.292–4.459), CKMB≥25.9 (CC, P = 0.040, OR = 3.118, 95% CI = 1.052–9.243; TC + TT, P = 0.003, OR = 2.825, 95% CI = 1.429–5.588), Glu ≥ 5.445 (TC + TT, P = 0.009, OR = 1.949, 95% CI = 1.179–3.225) and Cre ≥ 73.2 (CC, P = 0.017, OR = 2.221, 95% CI = 1.157–4.265) may be associated with CHD.

In the recessive model of rs1333045, compared to the reference group, the likelihood of onset is reduced in females (TT, P = 0.003, OR = 0.377, 95% CI = 0.198–0.717). However, when other factors, such as age ≥ 59.5 (TC + CC, P = 0.003, OR = 2.042, 95% CI = 1.268–3.291), smoking (TT, P = 0.029, OR = 1.916, 95% CI = 1.068–3.435; TC + CC, P = 0.003, OR = 1.983, 95% CI = 1.256–3.130), hypertension (TC + CC, P = 0.002, OR = 2.048, 95% CI = 1.291–3.247), CKMB≥25.9 (TT, P = 0.017, OR = 3.102, 95% CI = 1.226–7.846; TC + CC, P < 0.001, OR = 5.521, 95% CI = 2.684–11.357), Glu ≥ 5.445 (TT, P < 0.001, OR = 3.550, 95% CI = 1.785–7.061; TC + CC, P < 0.001, OR = 3.654, 95% CI = 2.263–5.900), Cre ≥ 73.2 (TT, P = 0.047, OR = 1.812, 95% CI = 1.007–3.260; TC + CC, P < 0.001, OR = 2.218, 95% CI = 1.399–3.517) exist, the risk of CHD increased.

The synergistic effect between dominant model, recessive model and classical risk factors

The analysis of the associations between gene and classical risk factors is shown in Table 8. We found a positive correlation between the dominant model of rs4977574 and age (S = 3.34, AP = 0.40, and AP* = 0.70), and the proportion of CHD attributable to the interaction between the dominant model, and age was as high as 70%. The same positive correlation was found between the dominant model and CKMB (S = 1.83, AP = 0.38, and AP* = 0.45). There is a positive interaction between Glu (S = 1.03, AP = 0.02, and AP* = 0.03), Cre (S = 1.91, AP = 0.28, and AP* = 0.48) and the dominant model of rs4977574. In the recessive model of rs4977574 and rs1333045, CKMB shows a positive interaction between them (rs4977574, S = 1.23, AP = 0.12, and AP* = 0.19; rs1333045, S = 1.86, AP = 0.38, and AP* = 0.46). In contrast, a negative interaction was found between HDL and the recessive model of rs4977574 (S = 0.67, AP = 0.73, and AP* = − 0.50) and between HDL and the dominant model of rs1333045 (S = 0.57, AP = 1.16, and AP* = − 0.76).

Table 8 The indexes of the synergistic effect between the dominant model, recessive model and risk factors

Haplotype analysis between rs4977574 and rs1333045

The LD plot on the left indicates the linkage of the gene SNP. R2 refers to the statistical correlation between the two sites. An LD plot was constructed based on the pairwise correlation between the two sites on the 9p21 locus, but unfortunately, it shows that there is no linkage disequilibrium between two sites (D’ = 0.92 and r2 = 0.75) (Fig. 2).

Fig. 2
figure2

Haplotype analysis between the two groups of rs4977574 and rs1333045

Discussion

CHD is a major cause of morbidity and mortality worldwide and represents a tremendous social and economic burden on society [28]. Recent genome-wide association studies have indicated that there is an association between increased susceptibility to coronary artery disease and specific SNPs within the genome that play a role in conjunction with other known traditional CHD risk factors, but the exact mechanism is unclear. Clinical observation has found that atherosclerosis is one of the major pathophysiological mechanisms of CHD [29]. Endothelial cells are monolayer continuous cells covering the inner surface of vessels with significant biological functions, including the regulation of thrombosis and coagulation, dilation of vascular smooth muscle, suppression of platelet adhesion and aggregation [30, 31]. The 9p21.3 risk allele in CHD is associated with altered expression of the CDKN2B-AS1 gene in blood. Evidence has demonstrated that CDKN2B-AS1 is a new susceptibility gene for CHD [32,33,34]. Both rs4977574 and rs1333045 are ANRILs of the CDKN2B-AS1 gene. Long noncoding RNAs have been shown to have a regulatory role in increasing cell proliferation and decreasing apoptosis in addition to participating in the inflammatory response.

In the present study, we recruited 855 inpatients, 598 cases with CHD confirmed by angiography as the CHD group and 257 subjects with normal coronary artery as the control group. In associating the two groups with SNPs, we analyzed their basic demographic information, including age, gender, smoking status, drinking status, and hypertension, and the following variable values, such as CKMB, TC, TG, HDL, LDL, ApoA, ApoB, Glu, Bun and Cre, which may be related to CHD.

The results of allelic and genotype logic regression showed that the G allele may be associated with a higher risk of CHD (P = 0.003, OR = 1.371, 95% CI = 1.113–1.689), and allele T was likely to reduce the occurrence of coronary events (P = 0.035, OR = 0.798, 95% CI = 0.649–0.982). The frequencies of genotypes GG (P = 0.042, OR = 1.421, 95% CI = 1.014–1.993) and AG (P = 0.004, OR = 1.849, 95% CI = 1.217–2.810) were increased in the CHD group compared with the control group. In comparison to genotype CC, genotype TC significantly reduced the incidence of CHD (P = 0.035, OR = 0.641 and 95% CI = 0.424–0.970). We also evaluated the association between the severity of coronary atherosclerosis and SNP polymorphisms. We found that the degree of vascular stenosis estimated by the Gensini score differed in relation to the rs4977574 genotype. We found that a significantly higher Gensini score appeared in the heterozygous AG genotype and the homozygous GG genotype compared with the AA genotype. This proves that allele G is positively correlated with the severity of CHD. In addition, we analyzed the distribution of clinical parameters in different models. In the recessive model of rs4977574, the genotype AA+AG carriers had significantly higher levels of ApoA (P = 0.028) than the GG genotype. In the dominant model of rs1333045, the serum levels of LDL in the CC genotype were higher than in the TC + TT genotype (P = 0.032), and the serum levels of ApoA in the CC genotype were significantly lower than in the TC + TT genotype (P = 0.049). The conclusion that allele A and allele T are protective factors of CHD should be confirmed more powerfully. Therefore, we suspected that this phenomenon was related to genes affecting the development of CHD by regulating lipid metabolism. Furthermore, we conducted an analysis of CHD prevalence and gene-environment interactions with the occurrence of CHD. It can be seen that genotype GG and AG are more susceptible to CHD than genotype AA when other variables are present. In comparison to genotype CC, genotype TC significantly reduced the incidence of CHD. When traditional risk factors and susceptibility genes coexist, the incidence of CHD is greatly increased.

Nevertheless, the other downside is that in haplotype analysis, it shows that there is no linkage disequilibrium between two sites, but the data showed a degree of association between the two sites that on the same gene. The deeper connection between the two sites needs to be explored further.

In the present study, due to the limited collection of clinical variables, further mechanisms exploring variables were not evaluated in the model building process. At the same time, the sample size of the case-control study was comparatively small, and the age and gender were not well matched. Moreover, the ratio of case and control was not 1:1 due to the consecutiveness of the recruited subjects and the invasiveness of coronary angiography. Therefore, we could not exclude a chance of random positive finding of the genotype association. Further studies utilizing a well-matched group with a larger sample size are warranted to increase the confidence of our findings. However, we can draw some conclusions from this study, and it also lays the foundation for further research in the future.

Conclusions

We identified that the SNPs tested in the case-control study were significantly associated with CHD. We proved that allele G of rs4977574 may be associated with a higher risk of CHD, and allele T of rs1333045 was likely to reduce the incidence of coronary events. When susceptibility factors (age ≥ 59.5, male, smoking, hypertension, CKMB≥25.9, Glu ≥ 5.445, Cre ≥ 73.2) and susceptibility genes (allele G or allele C) coexist, the incidence of CHD is greatly increased. However, further studies are needed to explain SNPs rs4977574 and rs1333045 in the specific pathogenesis of CHD.

Availability of data and materials

The datasets used or analysed during the current study are available from the corresponding author upon reasonable request.

Abbreviations

ANIRL:

Antisense noncoding RNA in the INK4 locus

ApoA:

Apolipoprotein A

ApoB:

Apolipoprotein B

AUC:

Area under curve

Bun:

Blood urea nitrogen

CHD:

Coronary heart disease

CI:

Confidence intervals

CKMB:

Creatinephosphokinase isoenzyme

Cre:

Creatinine

Glu:

Glucose

GWAS:

Genome-wide association studies

HDL-C:

High-density lipoprotein cholesterol

HWE:

Hardy-Weinberg equilibrium

LD:

Linkage disequilibrium

LDL-C:

Low-density lipoprotein cholesterol

LncRNA:

Long noncoding RNA

OR:

Odds ratio

PCR:

Polymerase chain reaction

RERI:

Relative excess risk of interaction

ROC:

Receiver operating characteristic

SAP:

Shrimp alkaline phosphatase

SE:

Standard error

SNPs:

Single nucleotide polymorphisms

TC:

Total cholesterol

TG:

Triglyceride

References

  1. 1.

    Owan TE, Roe MT, Messenger JC, Dai D, Michaels AD. Contemporary use of adjunctive thrombectomy during primary percutaneous coronary intervention for ST-elevation myocardial infarction in the United States. Catheter Cardiovasc Interv. 2012;80:1173–80.

    Article  Google Scholar 

  2. 2.

    Chien KL, Hsu HC, Su TC, Chang WT, Chen PC, Sung FC, et al. Constructing a point-based prediction model for the risk of coronary artery disease in a Chinese community: a report from a cohort study in Taiwan. Int J Cardiol. 2012;157:263–8.

    Article  Google Scholar 

  3. 3.

    Tada N, Maruyama C, Koba S, Tanaka H, Birou S, Teramoto T, et al. Japanese dietary lifestyle and cardiovascular disease. J Atheroscler Thromb. 2011;18:723–34.

    Article  Google Scholar 

  4. 4.

    Kangas-Kontio T, Huotari A, Ruotsalainen H, Herzig KH, Tamminen M, Ala-Korpela M, et al. Genetic and environmental determinants of total and high-molecular weight adiponectin in families with low HDL-cholesterol and early onset coronary heart disease. Atherosclerosis. 2010;210:479–85.

    CAS  Article  Google Scholar 

  5. 5.

    Feng L, Nian S-Y, Hao Y-L, Xu W-B, Ye D, Zhang X-F, et al. A single nucleotide polymorphism in the stromal cell-derived factor 1 gene is associated with coronary heart disease in Chinese patients. Int J Mol Sci. 2014;15:11054–63.

    CAS  Article  Google Scholar 

  6. 6.

    Kanu JS, Gu Y, Zhi S, Yu M, Lu Y, Cong Y, et al. Single nucleotide polymorphism rs3774261 in the AdipoQ gene is associated with the risk of coronary heart disease (CHD) in Northeast Han Chinese population: a case-control study. Lipids Health Dis. 2016;15:6.

    Article  Google Scholar 

  7. 7.

    Cunnington MS, Koref MS, Mayosi BM, Burn J, Keavney B. Chromosome 9p21 SNPs associated with multiple disease phenotypes correlate with ANRIL expression. PLoS Genet. 2010;6:e1000899.

    Article  Google Scholar 

  8. 8.

    Yap KL, Li S, Munoz-Cabello AM, Raguz S, Zeng L, Mujtaba S, et al. Molecular interplay of the noncoding RNA ANRIL and methylated histone H3 lysine 27 by polycomb CBX7 in transcriptional silencing of INK4a. Mol Cell. 2010;38:662–74.

    CAS  Article  Google Scholar 

  9. 9.

    Kotake Y, Nakagawa T, Kitagawa K, Suzuki S, Liu N, Kitagawa M, et al. Long non-coding RNA ANRIL is required for the PRC2 recruitment to and silencing of p15(INK4B) tumor suppressor gene. Oncogene. 2011;30:1956–62.

    CAS  Article  Google Scholar 

  10. 10.

    McPherson R, Pertsemlidis A, Kavaslar N, Stewart A, Roberts R, Cox DR, et al. A common allele on chromosome 9 associated with coronary heart disease. Science. 2007;316:1488–91.

    CAS  Article  Google Scholar 

  11. 11.

    Scott LJ, Mohlke KL, Bonnycastle LL, Willer CJ, Li Y, Duren WL, et al. A genome-wide association study of type 2 diabetes in Finns detects multiple susceptibility variants. Science. 2007;316:1341–5.

    CAS  Article  Google Scholar 

  12. 12.

    Psaty BM, O'Donnell CJ, Gudnason V, Lunetta KL, Folsom AR, Rotter JI, et al. Cohorts for heart and aging research in genomic epidemiology (CHARGE) consortium: design of prospective meta-analyses of genome-wide association studies from 5 cohorts. Circ Cardiovasc Genet. 2009;2:73–80.

    Article  Google Scholar 

  13. 13.

    Jarinova O, Stewart AF, Roberts R, Wells G, Lau P, Naing T, et al. Functional analysis of the chromosome 9p21.3 coronary artery disease risk locus. Arterioscler Thromb Vasc Biol. 2009;29:1671–7.

    CAS  Article  Google Scholar 

  14. 14.

    Liu Y, Sanoff HK, Cho H, Burd CE, Torrice C, Mohlke KL, et al. INK4/ARF transcript expression is associated with chromosome 9p21 variants linked to atherosclerosis. PLoS One. 2009;4:e5027.

    Article  Google Scholar 

  15. 15.

    Broadbent HM, Peden JF, Lorkowski S, Goel A, Ongen H, Green F, et al. Susceptibility to coronary artery disease and diabetes is encoded by distinct, tightly linked SNPs in the ANRIL locus on chromosome 9p. Hum Mol Genet. 2008;17:806–14.

    CAS  Article  Google Scholar 

  16. 16.

    Liu Y, Sanoff HK, Cho H, Burd CE, Torrice C, Ibrahim JG, et al. Expression of p16(INK4a) in peripheral blood T-cells is a biomarker of human aging. Aging Cell. 2009;8:439–48.

    CAS  Article  Google Scholar 

  17. 17.

    Sakalar C, Gurbuz E, Kalay N, Kaya MG. Higher frequency of rs4977574 (the G Allele) on chromosome 9p21.3 in patients with myocardial infarction as revealed by PCR-RFLP analysis. Tohoku J Exp Med. 2013;230:171–6.

    CAS  Article  Google Scholar 

  18. 18.

    Judkins MP. A percutaneous transfemoral technique. Radiology. 1967;89:819–21.

    Article  Google Scholar 

  19. 19.

    Sun X, Zhang M, Sanagawa A, Mori C, Ito S, Iwaki S, et al. Circulating microRNA-126 in patients with coronary artery disease: correlation with LDL cholesterol. Thromb J. 2012;10:16.

    Article  Google Scholar 

  20. 20.

    Lawler PR, Filion KB, Dourian T, Atallah R, Garfinkle M, Eisenberg MJ. Anemia and mortality in acute coronary syndromes: a systematic review and meta-analysis. Am Heart J. 2013;165:143–53.e5.

    Article  Google Scholar 

  21. 21.

    O'Connor NT, Cederholm-Williams S, Copper S, Cotter L. Hypercoagulability and coronary artery disease. Br Heart J. 1984;52:614–6.

    CAS  Article  Google Scholar 

  22. 22.

    Kannel WB, McGee DL. Diabetes and glucose tolerance as risk factors for cardiovascular disease: the Framingham study. Diabetes Care. 1979;2:120–6.

    CAS  Article  Google Scholar 

  23. 23.

    Hulley SB, Rosenman RH, Bawol RD, Brand RJ. Epidemiology as a guide to clinical decisions. The association between triglyceride and coronary heart disease. N Engl J Med. 1980;302:1383–9.

    CAS  Article  Google Scholar 

  24. 24.

    May AL, Kuklina EV, Yoon PW. Prevalence of cardiovascular disease risk factors among US adolescents, 1999-2008. Pediatrics. 2012;129:1035–41.

    Article  Google Scholar 

  25. 25.

    Boekholdt SM, Arsenault BJ, Mora S, Pedersen TR, LaRosa JC, Nestel PJ, et al. Association of LDL cholesterol, non-HDL cholesterol, and apolipoprotein B levels with risk of cardiovascular events among patients treated with statins: a meta-analysis. Jama. 2012;307:1302–9.

    CAS  Article  Google Scholar 

  26. 26.

    Tsimikas S, Hall JL. Lipoprotein(a) as a potential causal genetic risk factor of cardiovascular disease: a rationale for increased efforts to understand its pathophysiology and develop targeted therapies. J Am Coll Cardiol. 2012;60:716–21.

    CAS  Article  Google Scholar 

  27. 27.

    Gensini GG. A more meaningful scoring system for determining the severity of coronary heart disease. Am J Cardiol. 1983;51:606.

    CAS  Article  Google Scholar 

  28. 28.

    Zou JG, Ma YT, Xie X, Yang YN, Pan S, Adi D, et al. The association between CYP1A1 genetic polymorphisms and coronary artery disease in the Uygur and Han of China. Lipids Health Dis. 2014;13:145.

    Article  Google Scholar 

  29. 29.

    Madjid M, Willerson JT. Inflammatory markers in coronary heart disease. Br Med Bull. 2011;100:23–38.

    CAS  Article  Google Scholar 

  30. 30.

    Tasaki H. Low-density lipoprotein apheresis in the prevention of recurrent coronary heart disease: a review. Ther Apher Dial. 2003;7:408–12.

    CAS  Article  Google Scholar 

  31. 31.

    Campo G, Pavasini R, Pollina A, Tebaldi M, Ferrari R. Coagulation factors and recurrence of ischemic and bleeding adverse events in patients with acute coronary syndromes. Thromb Res. 2013;132:151–7.

    CAS  Article  Google Scholar 

  32. 32.

    Holdt LM, Beutner F, Scholz M, Gielen S, Gabel G, Bergert H, et al. ANRIL expression is associated with atherosclerosis risk at chromosome 9p21. Arterioscler Thromb Vasc Biol. 2010;30:620–7.

    CAS  Article  Google Scholar 

  33. 33.

    Peden JF, Hopewell JC. A genome-wide association study in Europeans and[ ]outh Asians identifies five new loci for coronary artery disease. Nat Genet. 2011;43:339–44.

    CAS  Article  Google Scholar 

  34. 34.

    Harismendy O, Notani D, Song X, Rahim NG, Tanasa B, Heintzman N, et al. 9p21 DNA variants associated with coronary artery disease impair interferon-gamma signalling response. Nature. 2011;470:264–8.

    CAS  Article  Google Scholar 

Download references

Acknowledgements

This study received support from the National Natural Science Foundations of China (No. 81170180, 30400173 and 30971257) and the Priority Academic Program Development of Jiangsu Higher Education Institutions. Dr. En-Zhi Jia is an Assistant Fellow at the Collaborative Innovation Center for Cardiovascular Disease Translational Medicine.

Funding

This study received support from the National Natural Science Foundations of China (No. 81170180, 30400173 and 30971257) and the Priority Academic Program Development of Jiangsu Higher Education Institutions.

Author information

Affiliations

Authors

Contributions

As the guarantor, LH, JXY, and EZJ conceived the study. SH, CHZ and QWJ initially drafted the paper. JZ, FHA, ZHC and LHL collected and analysed the data. LSW, WZM and GXX revised the manuscript. We ensure that all authors have read and approved the content of the manuscript.

Corresponding author

Correspondence to En-Zhi Jia.

Ethics declarations

Ethics approval and consent to participate

The study was performed in accordance with the protocol approved by the ethics committee of the First Affiliated Hospital of Nanjing Medical University and the Friendship Hospital of Ili Kazakh Autonomous Prefecture in China and all selected subjects signed informed consent.

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:

Table S1. CHD incidence by interactions with environmental factors.

Additional file 2:

Table S2. Serological biomarkers in the overdominant model.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Hua, L., Yuan, J., He, S. et al. Analysis on the polymorphisms of site RS4977574, and RS1333045 in region 9p21 and the susceptibility of coronary heart disease in Chinese population. BMC Med Genet 21, 36 (2020). https://doi.org/10.1186/s12881-020-0965-x

Download citation

Keywords

  • Coronary heart disease
  • SNP
  • Chromosome 9p21