Abstract

Purpose. Single nucleotide polymorphisms of the CYBA gene may modify the risk of coronary artery disease (CAD). The aim of the present study was to investigate whether the 49A>G (rs7195830) polymorphism is associated with CAD. Materials and Methods. CYBA gene 49A>G polymorphism was determined in 481 subjects: 242 patients with premature CAD and 239 age and sex matched controls using the fluorescently labeled allele-specific oligonucleotides method. Results. The frequency of the 49G allele carrier state was significantly higher in patients than in controls (84.8% versus 76.6%, resp., ), as well as the frequency of the 49G allele (62.2% versus 54.0%, ). Both factors were associated with CAD in the analyzed population (OR = 1.70, 95% CI: 1.04–2.76 for GG+AG versus AA and OR = 1.40, 95% CI: 1.08–1.83 for 49G versus  49A). Carrier state of the 49G allele was a stronger and independent risk factor for CAD among women (OR = 4.35, 95% CI: 1.50–13.20, ), as well as the 49G allele (OR = 2.25, 95% CI: 1.34–3.77, ). The 49G allele carrier state was also associated with left ventricular hypertrophy in patients with coronary artery disease (). Conclusion. The CYBA gene 49A>G polymorphism modifies the risk of coronary artery disease.

1. Introduction

NADPH oxidases (nicotinamide adenine dinucleotide phosphate oxidases) are one of the main sources of reactive oxygen species (ROS) in human body and especially in vasculature. NADPH oxidases are expressed in endothelium, vascular smooth muscle cells, fibroblasts, cardiomyocytes, and mainly phagocytes. Free radicals produced by NADPH oxidases are involved in the regulation of ROS-dependent intracellular signaling pathways, affecting expression of many classes of genes. In physiological conditions, ROS can stimulate growth, differentiation, migration, division, and apoptosis of cells of the vessel wall, to regulate vascular remodeling processes [15]. When synthesized in excessive levels, ROS may intensify the progress and development of atherosclerosis. Many data suggest the role of NADPH oxidases in the pathogenesis of cardiovascular disorders with atherosclerotic background, including coronary artery disease (CAD) [3, 69].

NADPH oxidases are multisubunit enzymatic complexes, which produce superoxide anion from molecular oxygen [6, 10]. Their active complexes are composed of NOX (NADPH oxidase) catalytic protein and many auxiliary subunits. There are few types of NOX proteins, namely, NOX2 (gp91phox) and its homologues (NOX1, NOX3–5). NOX2-based oxidases have the most important role in vascular superoxide production. Membrane-bound p22phox and cytoplasmic p47phox, Noxo1 (p47phox homologue), p67, Noxa1 (p67phox homologue), p40phox, and Rac GTP-ases are additional subunits of NADPH oxidases active complexes. p22phox and NOX form a heterodimer named cytochrome b558. During assembly of an active complex, cytosolic subunits migrate to the membrane-bound cytochrome b558 which becomes activated and generates superoxide anions [1117].

The p22phox is common element of NADPH oxidases complexes based on NOX1–NOX4 proteins [9]. It stabilizes the NOX-p22phox heterodimer formation and acts as a link between NOX and cytoplasmic subunits. p22phox is encoded by CYBA (cytochrome b245 alpha) gene (16q24). Opposite to genes encoding other subunits of NADPH oxidases, the CYBA gene is characterized by a relatively high polymorphic variability. It is probably due to the location of the CYBA gene, near a telomeric end of chromosome 16, and is a result of higher frequency of recombination than in loci of genes for other components of NADPH oxidases [8, 9].

In our previous studies, we demonstrated that single nucleotide polymorphisms (SNPs) of the CYBA gene may increase the risk of CAD in the Polish population [18] or intensify the effects caused by influence of traditional risk factors like cigarette smoking [1820], hypercholesterolemia [19, 20], or overweight/obesity [18]. The current study presents results from the next step of research on CYBA gene polymorphisms in the context of CAD and focuses on A>G (rs7195830) SNP. The A>G (NC_000016.9:g.88709712A>G, NM_000101.2:c. T>C) polymorphism lies in the 3′ UTR (3′ untranslated region) of the CYBA gene. To date, very little is known about this polymorphism and its potential functionality. Lack of previous studies with regard to its association with coronary artery disease prompted us to undertake such research. The aim of the present study was to investigate whether the A>G polymorphism shows an association with CAD and its clinical phenotype and whether it modifies the risk of disease associated with exposure to the traditional risk factors.

2. Materials and Methods

2.1. Subjects

We studied 481 Polish Caucasians, inhabitants of Upper Silesia. There were 242 patients with angiographically documented premature coronary artery disease in the CAD group (69 females and 173 males), aged 27–55 years (mean ). The second group included 239 blood donors (BD) matched by age and gender, with no signs of CAD and with negative familial history of CAD as inclusion criteria. CAD subjects were classified by the same cardiologist from (1) patients admitted to the 1st Department and Clinic of Cardiology at the Upper Silesian Center of Cardiology in Katowice and (2) patients admitted to the 1st Department of Cardiac Surgery at the Upper Silesian Center of Cardiology in Katowice. Blood donors were used as a control group. They were recruited from the Regional Centers of Blood Donation and Blood Treatment in Katowice and Raciborz. Blood was collected only from subjects with systolic blood pressure (SBP) < 140 mmHg and diastolic blood pressure (DBP) < 90 mmHg on the day of blood collection, according to recommendations of the Polish Centers of Blood Donation and Blood Treatment.

The criteria of inclusion and exclusion to the study as well as criteria for coronary artery disease, myocardial infarction (MI), and common risk factors of CAD were described previously [20].

The study was approved by the Ethics Committee of the Medical University of Silesia in Katowice. All submitted cases signed an informed consent prior to the study.

2.2. Methods
2.2.1. Serum Lipids Measurement

Serum lipid parameters like triglycerides (TG), total cholesterol (TC), and HDL cholesterol (HDL-C) were measured using colorimetric methods (Analco, Warsaw, Poland). The levels of LDL cholesterol (LDL-C) were calculated using Friedewald formula [21] only in subjects with TG levels < 4.4 mmol/L.

2.2.2. DNA Extraction and Genotyping

Genomic DNA was extracted using the MasterPure genomic DNA purification kit (Epicentre Technologies, Madison, USA). The CYBA gene A>G polymorphism (rs7195830) was genotyped using the TaqMan® Predesigned SNP Genotyping Assay (Applied Biosystems, Foster City, USA). Primer-probe sets were designed and prepared by the manufacturer. 7300 Real-Time PCR System (Applied Biosystems, Foster City, USA) was used to genotyping the rs7195830 polymorphism. Correctness of the genotyping was verified by regenotyping 15% of the samples and the repeatability was 100%.

2.2.3. Statistical Analysis

Statistica 10.0 software (STATSOFT, Tulsa, USA) was used in statistical analyses. Normality of distribution of quantitative data was assessed by the Shapiro-Wilk test and then a comparison of data was performed by Student’s t-test (variables with normal distribution) or the Mann-Whitney U test (variables with nonnormal distribution).

The χ2 test was used to test the Hardy-Weinberg equilibrium in all groups. Genotypes and alleles frequencies were also compared by χ2 test as well as the presence of possible associations between genotype variants of the A>G polymorphism and characteristics of clinical phenotype of CAD. Odds ratios (OR) and their 95% confidence intervals (CI) were estimated using a univariate analysis. In univariate analyses statistical significance was accepted at . Multiple logistic regression analysis was performed after adjustment for age, gender, and traditional risk factors of CAD. Only variables with values of less than 0.1 in univariate analysis were classified to the multiple logistic regression model. Risk ratio (RR) values with 95% CI were used when the number of subjects in any of subgroups was 0. Epi Info™ 7.1.1.0 (Centers for Disease Control and Prevention, Atlanta, USA) was used to compute the effective sample size and statistical power.

To determine possible interactions between CYBA genotypes and traditional risk factors of CAD (synergism/antagonism), the 4 × 2 table approach was used. Multiplicative synergy indices (SIMs) were calculated on the basis of OR values taken from 4 × 2 table comparisons [22]. 95% confidence intervals for SIMs were calculated using the script described by Cortina-Borja et al. [23]. Additive interactions were analyzed independently by a χ2 test using ORs as measures of interaction.

3. Results

3.1. Clinical Characteristics and Traditional Risk Factors of CAD

There were 72.5% cases who had suffered from MI in the CAD group. 65.8% of patients had multivessel disease (stenoses in at least two coronary vessels). Critical stenoses (>90%) in coronary vessels were diagnosed in 66.9% patients. 12% of patients with CAD suffered from a left ventricular hypertrophy (LVH). Other clinical and biochemical parameters of patients and controls are shown in Table 1. CAD patients had increased levels of TC, LDL-C, and TG and a higher BMI. HDL-C levels were significantly lower in CAD patients. Hypertension, cigarette smoking, and diabetes mellitus were significantly more frequent in CAD subjects (Table 1).

3.2. Association of the CYBA Gene 49A>G Polymorphism with CAD

Genotype frequencies of the CYBA gene A>G polymorphism were compatible with the Hardy-Weinberg equilibrium in CAD group and controls (, , resp.), as well as in subgroups divided by gender ( in CAD females, in CAD males and in BD females, in BD males).

The frequency of the G allele carriers (GG+AG genotypes) was significantly higher in patients than in controls as well as the frequency of the G allele (Table 2). The GG genotype was also more frequent in the CAD group than in the controls; however the difference obtained in univariate model lied on the border of statistical significance (). Both the G allele carrier state and GG homozygosity were associated with premature CAD in multivariate model, after adjustment for sex, age, TC, LDL-C, HDL-C, TG, BMI, diabetes mellitus, cigarette smoking status, hypertension, and familial history of coronary artery disease (Table 2).

Further analyses showed that the differences in the frequencies of genotypes and alleles of the A>G CYBA polymorphism were more significant in female subgroups of CAD patients and blood donors (Table 3). The G allele carrier state and GG homozygosity were associated with premature CAD in multivariate model, also in this case (Table 3). These associations were not observed in male subgroups (data not shown).

3.3. The CYBA Gene 49A>G Polymorphism and Traditional Risk Factors of CAD

The presence of possible interactions (synergy/antagonism, additive and/or multiplicative effects) between genotype variants of the CYBA gene A>G polymorphism and traditional risk factors of CAD (cigarette smoking, BMI, elevated levels of TC, and LDL) was checked. We did not find any statistically significant interactions between A>G polymorphism genotypes and traditional risk factors of CAD (data not shown).

We also checked whether the A>G polymorphic variants are associated with concomitant phenotypes of CAD like hypertension, dyslipidemias, and diabetes mellitus, using case-control model. We found an association between rs7195830 SNP variants and hypertension analysing subgroups of hypertensive CAD patients () and nonhypertensive blood donors (). The frequencies of GG genotype, G allele carriers (subjects with GG+AG genotypes), and the G allele were significantly higher in hypertensive CAD patients than in controls (Table 4). There were no statistically significant differences in alleles and genotypes distribution between nonhypertensive CAD patients and nonhypertensive blood donors (data not shown).

The presence of possible associations between variants of the A>G polymorphism and characteristics of clinical phenotype of CAD was also tested. The G allele was present in 100% of CAD patients with left ventricular hypertrophy and its frequency in patients without LVH was 82.6%. This difference was statistically significant (Table 5). There was no association between genotype variants of the A>G polymorphism and MI or the severity of atherosclerosis estimated on the basis of the number of coronary stenoses or critical arterial occlusions observed during a coronary angiography (data not shown).

4. Discussion

The CYBA gene encodes p22phox, a subunit of NADPH oxidases, enzymes critical in vascular production of superoxide anion. In the current work the A>G CYBA gene polymorphism was analyzed in the context of coronary artery disease in the classic case-control model. The influence of the A>G SNP on CAD risk and possible interactions with traditional risk factors of CAD were searched. We demonstrated statistically significant differences in genotype and allele frequencies of the A>G polymorphism between analyzed groups. G allele carriers were more frequent in CAD patients than in the control group and G allele carrier state was associated with CAD. We also showed that carrier state of the G allele is a stronger and independent risk factor for CAD among women. Present work also demonstrated that the G allele carrier state is associated with left ventricular hypertrophy in patients with coronary artery disease.

This last observation prompted us to check whether the rs7195830 polymorphism was associated with hypertension, as the left ventricular hypertrophy is most often the result of high blood pressure. We observed that the GG genotype, G allele, and its carrier state were independent risk factors of hypertension, when we compared extracted subgroups of CAD patients (only hypertensive subjects) and blood donors (only nonhypertensive controls). This effect was stronger than that observed between the analyzed polymorphism and CAD, which suggests that the rs7195830 is a risk factor for hypertension rather than for CAD.

Hypertension is one of the main causes of left ventricular hypertrophy. It seems that there is positive feedback between hypertension-induced myocardial hypertrophy and NADPH oxidases. On the one hand NADPH oxidases of the neuronal, renal, and vascular systems are strongly implicated in pathological signaling leading to hypertension (discussed in Datla and Griendling [24]). It is well documented that NOX are a primary source of superoxide in angiotensin II-induced neuronal and vascular activity and superoxide induces vascular dysfunction in hypertension by its well-described interaction with NO. On the other hand, hypertension is known activator of NADPH oxidases. NADPH oxidases may be involved in cardiac hypertrophy, because their expression and activity are upregulated by pressure overload, in both cardiomyocytes and endothelial cells of rodents [25, 26].

Until today, the CYBA gene A>G polymorphism was not related to hypertension; however, other CYBA gene SNPs, namely, -930A>G and -675A>T, were previously associated with this phenotype [3, 8, 27, 28]. It was also shown that the T allele of the 214C>T polymorphism is associated with higher left ventricular mass in Brazilian hypertensive patients [29].

Lack of functional studies and insufficient number of clinical trials greatly hinder possibility of comprehensive interpretation of the results of our work. Although the function of the A>G polymorphism is not yet known, Fan et al. [30] suggested that the rs7195830 lies in a miRNA binding region of the 3′ UTR of CYBA. Because binding of miRNA to mRNA is essential for mRNA level regulation and protein expression, it is possible that allelic variants of the A>G polymorphism may differentially affect miRNA-mediated efficiency and a rate of these processes. Although this is only a hypothesis, the neighboring polymorphism of the 3′ UTR of the CYBA gene, namely, rs1049255 (G>A), lies within the miRNA binding site [31]. The A allele increases the potential of hybridization between CYBA mRNA and miRNA (miR-320a) so lower levels of CYBA expression and NOX activity observed in A allele may result from enhanced regulation of miRNA [31]. It is possible that an association between G allele and CAD observed in the present work is a consequence of the interactions between these two 3′ UTR polymorphisms, especially as the distance between the G>A and A>G SNPs is only 25 base pairs. It should be noted, however, that we did not observe preferential coexistence of alleles of both polymorphisms and they did not create haplotype blocks in our studied groups (data not shown). In our previous study, however, we documented that the G allele influences predispositions to CAD through interactions with cigarette smoking and hypercholesterolemia with no direct association with the disease [20].

To date the A>G CYBA polymorphism was analyzed only in the aspects of artery elasticity [32], chronic kidney disease [30], and cervical cancer progression [33]. In the study analyzing large and small artery elasticity of the apparently healthy Chinese individuals the GG and AG genotypes had higher small artery elasticity compared with the AA genotype (), in physically active participants. There was no such association in less physically active participants [32]. Next work determined whether the A>G SNP influences the estimated glomerular filtration rate (eGFR) and its role in the pathogenesis of chronic kidney disease in Chinese population [30]. The results indicated that the AA homozygosity of the A>G polymorphism was associated with significantly lower eGFRs () compared with the G allele carrier state. The A allele carrier state was associated with an increased risk of intraepithelial neoplasia grade 3 or cervical cancer, in subjects with papillomavirus infections [33].

As seen, it is difficult to relate our findings to the results of the previous studies, mainly due to the difference of studied phenotypes/disease entities and ethnic differences between the studied populations. The latter factor appears to be particularly important. In some cases, the presence of differences in allele and genotype frequencies of SNPs results in contradictory findings of association studies from populations with different ethnicity. The results of the meta-analyses of another often tested CYBA polymorphism, namely, 242C>T (rs4673), can be an example of the foregoing. In Caucasians, the TT genotype was associated with an increased risk of ischemic cerebrovascular disease (OR = 2.13, 95% CI: 1.06–4.26, ) [34] and marginally with CAD (OR = 1.22, 95% CI: 0.99–1.49, ) [35]. In Asians, the T allele has a protective effect to coronary artery disease [36, 37] and ischemic cerebrovascular disease [34].

Summing up the results of our current work, carrier state of the G allele of the CYBA gene is associated with coronary artery disease, hypertension, and left ventricular hypertrophy in patients with coronary artery disease. The A>G polymorphism probably does not modify the risk of CAD through interactions with known, traditional risk factors, unlike other CYBA gene polymorphisms [1820]. The results of the current work tend to indicate that the A>G polymorphism may be particularly strong risk factor in females. Since neither current nor previous studies explain the role of the A>G polymorphism and its relationship with hypertension and atherosclerotic CAD, further functional researches and population surveys are needed.

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

Acknowledgments

This research was supported by a grant from the Medical University of Silesia KNW-1-014/N/4/0.