Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Both Paraoxonase-1 Genotype and Activity Do Not Predict the Risk of Future Coronary Artery Disease; the EPIC-Norfolk Prospective Population Study

  • Rakesh S. Birjmohun ,

    r.s.birjmohun@amc.uva.nl

    Affiliation Department of Vascular Medicine, Academic Medical Center, Amsterdam, The Netherlands

  • Menno Vergeer,

    Affiliation Department of Vascular Medicine, Academic Medical Center, Amsterdam, The Netherlands

  • Erik S. G. Stroes,

    Affiliation Department of Vascular Medicine, Academic Medical Center, Amsterdam, The Netherlands

  • Manjinder S. Sandhu,

    Affiliation Department of Public Health and Primary Care, Institute of Public Health, University of Cambridge, Cambridge, United Kingdom

  • Sally L. Ricketts,

    Affiliation Department of Public Health and Primary Care, Institute of Public Health, University of Cambridge, Cambridge, United Kingdom

  • Michael W. Tanck,

    Affiliation Department of Clinical Epidemiology and Biostatistics, Academic Medical Center, Amsterdam, The Netherlands

  • Nicholas J. Wareham,

    Affiliation Medical Research Council Epidemiology Unit, Cambridge, United Kingdom

  • J. Wouter Jukema,

    Affiliation Department of Cardiology, Leiden University Medical Center, Leiden, The Netherlands

  • John J. P. Kastelein,

    Affiliation Department of Vascular Medicine, Academic Medical Center, Amsterdam, The Netherlands

  • Kay-Tee Khaw,

    Affiliation Department of Public Health and Primary Care, Institute of Public Health, University of Cambridge, Cambridge, United Kingdom

  • S. Matthijs Boekholdt

    Affiliation Department of Cardiology, Academic Medical Center, Amsterdam, The Netherlands

Abstract

Background

Paraoxonase-1 (PON1) is an antioxidant enzyme, that resides on high-density lipoprotein (HDL). PON1-activity, is heavily influenced by the PON1-Q192R polymorphism. PON1 is considered to protect against atherosclerosis, but it is unclear whether this relation is independent of its carrier, HDL. In order to evaluate the atheroprotective potential of PON1, we assessed the relationships among PON1-genotype, PON1-activity and risk of future coronary artery disease (CAD), in a large prospective case-control study.

Methodology/Principal Findings

Cases (n = 1138) were apparently healthy men and women aged 45–79 years who developed fatal or nonfatal CAD during a mean follow-up of 6 years. Controls (n = 2237) were matched by age, sex and enrollment time. PON1-activity was similar in cases and controls (60.7±45.3 versus 62.6±45.8 U/L, p = 0.3) and correlated with HDL-cholesterol levels (r = 0.16, p<0.0001). The PON1-Q192R polymorphism had a profound impact on PON1-activity, but did not predict CAD risk (Odds Ratio [OR] per R allele 0.98[0.84–1.15], p = 0.8). Using conditional logistic regression, quartiles of PON1-activity showed a modest inverse relation with CAD risk (OR for the highest versus the lowest quartile 0.77[0.63–0.95], p = 0.01; p-trend = 0.06). PON1-activity adjusted for Q192R polymorphism correlated better with HDL-cholesterol (r = 0.26, p<0.0001) and more linearly predicted CAD risk (0.79[0.64–0.98], p = 0.03; p-trend = 0.008). However, these relationships were abolished after adjustment for HDL (particles-cholesterol-size) and apolipoproteinA-I (0.94[0.74–1.18], p-trend = 0.3).

Conclusions/Significance

This study, shows that PON1-activity inversely relates to CAD risk, but not independent of HDL, due to its close association with the HDL-particle. These data strongly suggest that a low PON1-activity is not a causal factor in atherogenesis.

Introduction

Oxidation of low-density lipoprotein (LDL) and the subsequent generation of lipid hydroperoxides has been shown to accelerate the development of atherosclerosis and to adversely affect cardiovascular outcome [1][3]. Consequently, a robust number of antioxidants that potently inhibit LDL-oxidation in vitro [4], [5] have been tested in clinical trials of impressive size, but most of these studies have yielded negative results [6][9]. Consequently, attention has shifted towards other physiological targets that may play a role in LDL-oxidation in vivo. In this respect, paraoxonase1 (PON1), a plasma enzyme that resides on the high-density lipoprotein (HDL) particle, has been shown to inhibit LDL-oxidation and is therefore also a strong candidate as a protective enzyme in the development of atherosclerosis [10][12].

Interestingly, two genetic variants of PON1, the Q192R and L55M polymorphisms (rs662 and rs854560, respectively) profoundly affect the activity of serum paraoxonase. The PON1-192 QQ gene product exhibits low hydrolytic activity towards paraoxon, whereas that of the RR-genotype shows high activity towards paraoxon [13], [14]. Conversely, the isoenzyme corresponding to this RR-genotype exhibits low hydrolytic activity towards lipid hydroperoxides. Likewise, gene products of the PON1-55 LL-genotype, LM-genotype and MM-genotype exhibit high, intermediate and low hydrolytic activity towards paraoxon, respectively. Although many association studies have reported on PON1-gene polymorphisms and CAD risk, these studies have yielded conflicting results; these reports were statistically underpowered, studied different ethnic and patient populations, utilized different sampling strategies as well as different genotyping procedures. These impediments have been discussed in a meta-analysis of association studies on PON1-genotype and CAD risk [15]. In contrast to PON1-genotype, only few studies have assessed the relationship between PON1-activity and CAD risk. In a cross-sectional study, hydrolysis rates of PON1 towards paraoxon and diazoxon were significantly lower in people with carotid disease than in controls [13]. Consistently, the prospective Caerphilly study showed that low PON1-activity towards paraoxon was associated with an increased risk of future CAD [14]. In line, PON1 activity and concentration were significantly lower in patients with severe coronary artery stenosis [16], whereas another study found no association between PON1 activity and coronary artery calcification [17]. In addition to this, elevated levels of PON1 activity healthy middle-aged women were shown to be associated with increased CAD risk [18]. These conflicting data have made it difficult to fully comprehend the role of PON1 in atherosclerosis. In fact, it remains uncertain whether the relation between PON1-activity and CAD risk is causal and whether it is independent of other CAD risk factors such as HDL cholesterol [15]. We used the concept of Mendelian randomization to evaluate whether PON1-activity may play a causal role in determining CAD risk [19]. We quantified the activity of PON1 towards paraoxon, determined two prominent genetic variants affecting PON1-activity, and assessed the consistency and proportionality of relationships between genotype, phenotype and risk of future CAD. We studied these relationships in a case-control study nested in the large, prospective EPIC (European Prospective Investigation into Cancer and Nutrition)–Norfolk study cohort.

Results

Serum samples from 1138 cases and 2237 controls were included in the analysis from which 1099 cases were matched to two controls and 39 cases were matched to one control. The Q192R and L55M polymorphisms were available for 1077 cases (61 missing cases). The numbers in tables 1, 2, 3 may not add up to the total number because of missing data. Matching ensured that the distribution of age and sex was comparable between cases and controls. As expected, individuals who developed CAD during follow-up were more likely than controls to smoke and to have diabetes mellitus (Table 1). The use of alcohol was comparable between cases and controls. Likewise, in both men and women, levels of total cholesterol, LDL-cholesterol, triglycerides (non-fasting), apoB, systolic and diastolic blood pressure, BMI, waist circumference and CRP were significantly higher in cases than in controls, whereas HDL-cholesterol and apoA-I were significantly lower in cases than controls. Serum PON1-activity was similar in cases and controls (60.7±45.3 versus 62.6±45.8 U/L, P = 0.3). Sex-specific analyses revealed similar results. Interestingly, women had higher PON1-activity than men, consistent with higher HDL cholesterol and apoA-I levels in women (see data supplement for sex-specific data: Table S1). The distribution of PON1-Q192R and PON1-L55M polymorphisms did not differ between cases and controls (Table 1).

thumbnail
Table 2. Associations between PON1 activity, PON1 genotype and risk factors.

https://doi.org/10.1371/journal.pone.0006809.t002

thumbnail
Table 3. Distribution of cardiovascular risk factors by quartiles of PON1 activity.

https://doi.org/10.1371/journal.pone.0006809.t003

Table 2 shows the correlations between PON1-activity and risk factors and associations between both PON1-genotypes and risk factors. In brief, PON1-activity was positively correlated with HDL-particle number, HDL-cholesterol, HDL-size, apoA-I, LDL-cholesterol and vitamin-C concentration. There was a negative but weak correlation with levels of CRP and MPO. Similar trends were found among men and women in the sex-specific analyses (see data supplement for sex-specific data: Table S2).

Although both genotypes were significantly associated with PON1-activity, we found a stronger association for the PON1-Q192R polymorphism. Carriership of more PON1-192R alleles or more PON1-55M alleles was associated with a higher serum PON1-activity, concordant with a gene-dosing effect. Distribution of PON1-activity over genotypes was triphasic according to the PON1-Q192R polymorphism, illustrating the strong influence of this variant on the PON1-activity phenotype (Figure 1). Statistical adjustment for the PON1-Q192R polymorphism yielded a normally distributed PON1-activity (Figure 2), which serves as a proxy for PON1-concentration. PON1-activity adjusted for the PON1-Q192R polymorphism was similar in cases and controls (0.49±1.09 versus 0.76±0.79 AU/L, P = 0.8).

thumbnail
Figure 1. Histogram showing the triphasic distribution of PON1-activity according to the PON1-Q192R polymorphism.

https://doi.org/10.1371/journal.pone.0006809.g001

thumbnail
Figure 2. Histogram showing the normal distribution of PON1-activity after adjusting for the PON1-Q192R polymorphism.

https://doi.org/10.1371/journal.pone.0006809.g002

Table 3 shows the distribution of cardiovascular risk factors as stratified by PON1-activity quartiles (Table 3) and the PON1-Q192R polymorphism adjusted PON1-activity quartiles (Table 4). In brief, in both men and women, the PON1-Q192R polymorphism adjusted PON1-activity showed a strong association with HDL-cholesterol, HDL-particle number and apoA-I. There was a modest inverse, but significant linear relationship between the PON1-Q192R polymorphism adjusted PON1-activity and waist circumference, CRP, MPO and age, whereas a modest positive linear relationship was observed for gender, diabetes mellitus, total cholesterol, LDL-cholesterol, plasma vitamin-C concentration and HDL-size. Similar, although weaker, associations were found for the unadjusted PON1-activity quartiles (Table 4).

thumbnail
Table 4. Distribution of cardiovascular risk factors by quartiles of PON1-activity adjusted for PON1-192 genotype.

https://doi.org/10.1371/journal.pone.0006809.t004

Table 5A shows the unadjusted and adjusted ORs for future CAD by quartiles of PON1-activity. Serum PON1-activity was modestly, although not linearly related to risk of future CAD. Adjustment for HDL-parameters (HDL-cholesterol, HDL-particle number, HDL-size and apoA-I), did not considerably affect the risk estimate. Further adjustment for risk factors (alcohol use, BMI, CRP, diabetes mellitus, fasting time, LDL-cholesterol, MPO, smoking, systolic blood pressure, plasma vitamin-C, vitamin supplement use, waist circumference and triglycerides (non-fasting) did not substantially alter the observed relationship.

thumbnail
Table 5. Risk of future coronary artery disease by quartiles of PON1 activity and by quartiles of PON-1 activity adjusted for PON1-192 genotype.

https://doi.org/10.1371/journal.pone.0006809.t005

Table 5B shows, however, that the PON1-Q192R polymorphism adjusted PON1-activity was strongly inversely related with the risk of future CAD (OR 0.79 for the highest versus the lowest quartile, 95% CI 0.64 to 0.98, P for linearity = 0.008). This relation was abolished by adjustment for HDL related parameters (HDL-cholesterol, HDL-particle number, HDL-size and apoA-I). Additional adjustments for cardiovascular risk factors did not affect the results. Similar trends were found among men and women in the gender-specific analyses (see data supplement for gender-specific data: Table S3).

Backward stepwise regression analysis for PON1-activity, PON1-genotype as well as PON1-activity adjusted for PON1-192 and PON1-55 genotype revealed no significant effect of these variables and CAD risk (see data supplement: Table S4, S5, S6).

Haplotypes of PON1-192 and PON1-55 genotypes where equally distributed between cases and controls. Further the effect of these different haplotypes on PON1 activity, PON1-genotypes and PON1-activity adjusted for PON1-192 genotype are reported in the supplementary table (see data supplement: Table S7). No significant relation between combined haplotypes and CAD risk was found.

Discussion

In this large prospective case-control study among apparently healthy men and women, we observed that PON1-activity was modestly and inversely related with the risk of future CAD. However, given the fact that the Q192R polymorphism strongly affects lifelong PON1-activity, the absence of a relation between the Q192R polymorphism and CAD risk suggests that PON1-activity is not a causal factor in atherogenesis. The apparent protective function of PON1 can only be appreciated against the background of its close association with HDL.

Whereas PON1-L55M genotype was modestly associated with PON1-activity in our study, the PON1-Q192R polymorphism was the strongest determinant of PON1-activity against paraoxon. Although higher PON1-activity against paraoxon tended to be associated with a lower CAD risk, this relation was not linear. More importantly, PON1-Q192R genotype had no effect on CAD risk. In the context of Mendelian randomization [19], these findings suggest that the relation between PON1-activity against paraoxon and CAD risk is not causal. Had this been the case, one would expect consistent and proportional relationships among genotype, phenotype and CAD risk. We observed that carriership of a 192R allele was associated with a 62.5 U/l increase in PON1-activity per allele. In turn, an increase of 62.5 U/l in PON1-activity corresponds approximately to the difference between the first and third PON1 activity quartile, which was associated with an OR for CAD of 0.81(95% CI [0.67–0.99]). However, carriership of a 192R allele was not associated with a proportional decrease in CAD risk. In fact, we observed a null relationship (OR 0.98[0.84–1.15], p = 0.8). The present findings imply that previous findings of an inverse relation between PON1 activity and CAD risk were troubled by confounding. In this respect, several issues deserve closer attention. PON1 is located predominantly on HDL-particles. Since HDL has a strong inverse relation to CAD risk [24], other protective effects mediated by or associated with HDL [25] are likely to contribute to the PON1–CAD relation. People with more circulating HDL-particles can be expected to have a higher plasma PON1-concentration. Because PON1-activity is a function of PON1-concentration and PON1-genotype [26][27], we could obtain a proxy for PON1-concentration by statistically adjusting PON1-activity for the PON1-Q192R polymorphism. Indeed, the PON1-Q192R polymorphism adjusted PON1-activity correlated better with HDL-related parameters than unadjusted PON1-activity did. Similarly, it was a better predictor of CAD risk than unadjusted PON1-activity. However, statistical adjustment for HDL-related parameters abolished the inverse relation between PON1-activity adjusted for the PON1-Q192R polymorphism and CAD risk in the present cohort (Table 4B). These findings lend further support to the concept that the relation between PON1-activity and CAD risk to a relevant extent confounded by HDL-cholesterol levels and is not based on causality.

PON1 and CAD risk

Numerous studies have reported associations between the PON1-Q192R and L55M polymorphisms and risk of CAD [28][34], whereas other studies found no such relation [35][39]. However, a meta-analysis reported no significant association for the L55M variant [15]. In addition, it observed that the significant association for the Q192R variant was largely explained by small studies while larger studies which are less prone to statistical bias, showed no significant association. Interestingly, recent genome-wide association studies support our findings and have not revealed any associations between the PON1-genotype and risk of CAD [40], [41]. In addition, Jarvik et al. have shown that PON1-activity towards paraoxon is a better predictor of cardiovascular risk than PON1 genotype [13]. In the present study, we substantiate this finding, since no relationship was observed between the PON1-Q192R or PON1-L55M polymorphisms and risk of future CAD events. In this respect, it should be noted that the present study was a large prospective case-control cohort, comprising 1138 cardiovascular events, whereas other observational studies have used a small sample size and/or a retrospective design.

Finally, the role of inflammation in atherosclerosis has been substantiated in numerous studies. However, in the present study, inflammatory markers such as CRP and MPO did not affect the relation between PON1 and CAD risk. An explanation for this finding might be the modest relation between both CRP and MPO versus PON1 activity.

PON1 and Vitamins

Vitamin C intake has previously been shown to be associated with increased PON1 activity [42]. In line, the present study demonstrated that PON1 activity is primarily genotype dependent, and is just modestly influenced by plasma vitamin C concentration. It has been suggested that a reduction in oxidative stress related to vitamin C intake may preserve PON1 activity. In addition, it has been shown that those subjects using vitamin supplements are more likely to engage in healthy lifestyles, of which the factors are multivariate [43]. Thus, it is possible that some factor correlated with vitamin intake is actually influencing PON1 activity.

Study Limitations

The results of the present study should be considered in light of its potential limitations. First, case ascertainment is an issue in the design of every prospective study, including this one. However, a validation study indicated that case ascertainment in our study was at least equivalent to that of other large prospective cohort studies [44]. Second, the findings of the present study apply to Caucasians and should be carefully extrapolated to non-Caucasian populations, especially since PON1-genotype frequencies have been reported to vary substantially between Caucasian and East Asian populations (PON1-192RR homozygosity occurs in only about 10% of Caucasians, but in about 40% of East Asians). Also, activities of PON1 towards other substrates than paraoxon, e.g. lipid hydroperoxides, might have a distinctly different impact on atheroprotection. Therefore we cannot exclude that other tests for PON1 activity would yield different results. However, since the PON1-192 RR genotype, apart from high reactivity towards paraoxon, is also characterized by low reactivity towards lipid hydroperoxides, the absence of a relationship between the PON1-Q192R polymorphism and CAD risk suggests that this particular phenotype can not strongly modulate CAD risk either. Third, the in vitro hydrolysis of paraoxon reflecting PON1 activity may be a suboptimal method to measure its actual in vivo antioxidant capacity. Methods to more accurately PON1 in vivo antioxidant capacity are currently not at hand.

Fourth, a common limitation of case-control studies, even when they are of a prospective nature is that subjects who contract CAD at a future date are bound to have more classical risk factors than subjects who do not. This may limit the contribution of PON1-activity and genotype on CAD risk. Despite these limitations we tried to minimize the influence of these stronger risk factors on CAD risk to an acceptable degree (see stepwise backward regression analysis in supplementary tables S4-S6). Finally, the Mendellian randomization analyses in the present study should be interpreted with caution, since the PON1-activity and CAD risk association is nonlinear. The concept of Mendelian randomization is based on several assumptions, including a linear relation between the trait of interest and disease risk.

Conclusion

We demonstrate a modest and inverse relationship between PON1-activity and risk of future CAD among apparently healthy men and women. The PON1-Q192R polymorphism strongly affects PON1-activity but was not associated with risk of future CAD. The inconsistency of the relationships among PON1-genotype, PON1-activity and CAD risk suggests that PON1 is not a causal factor in the development of CAD. Concurrently, we observed that the relationship between a proxy for PON1-concentration and CAD risk was abolished upon adjustment for HDL-related parameters, indicating that this relationship is confounded by HDL-cholesterol.

Materials and Methods

The EPIC-Norfolk cohort study is a prospective population study of 25663 male and female inhabitants of Norfolk, United Kingdom, aged between 45 and 79 years old. EPIC-Norfolk is part of the 10-country collaborative EPIC study designed to investigate dietary and other determinants of cancer. Additional data were obtained to enable assessment of determinants of other diseases. Participants were enrolled between 1993 and 1997, completed a baseline questionnaire survey and attended a clinic visit. The study cohort was similar to UK population samples with regard to many characteristics, including anthropometry, blood pressure, and lipids, but with a lower proportion of smokers. Participants were recruited by mail from age-gender registers of general practices. At the baseline survey between 1993 and 1997, participants completed a detailed health and lifestyle questionnaire, and additional data collection was performed by trained nurses at a clinic visit as described previously. All individuals have been flagged for mortality at the UK Office of National Statistics, with vital status ascertained for the entire cohort. Death certificates of individuals enrolled in the study were used to classify the deaths. The death certificates were coded by trained nosologists according to the International Classification of Diseases (ICD) 9th revision. Death was considered due to CAD if the underlying cause was coded as ICD 410 to 414. In addition, participants admitted to a hospital were identified by their unique National Health Service number by data linkage with ENCORE (East Norfolk Health Authority database), which identifies all hospital contacts throughout England and Wales for Norfolk residents. Participants were identified as having CAD during follow-up if they had a hospital admission and/or died with CAD listed as an underlying cause. We report results with follow-up up to January 2003, an average of 6 years. The Norwich District Health Authority Ethics Committee approved the study, and all participants gave signed informed consent.

Participants

For the present nested case-control study, we identified study participants who did not report a history of heart attack or stroke at the baseline clinic visit. A total of 1138 apparently healthy individuals developed fatal or nonfatal CAD during follow-up. A total of 2237 controls were apparently healthy study participants who remained free of CAD during follow-up. 1099 cases were matched to two controls and 39 cases were matched to one control by sex, age (within 5 years), and date of visit (within 3 months). Data on PON1-activity was available in all study subjects, but data on PON1-genotype was missing in 61 cases.

Biochemical Analyses

Non-fasting blood samples were taken by vein puncture into serum tubes. Serum levels of total cholesterol, HDL cholesterol, and triglycerides were measured on fresh samples with the RA 1000 auto-analyzer (Bayer Diagnostics, Basingstoke, United Kingdom). LDL-cholesterol levels were calculated with the Friedewald formula [20]. Blood samples were also stored at minus 80° Celsius. Samples were later thawed and processed for analysis. Serum levels of apolipoproteinA-I (apoA-I) and B (apoB) were measured by rate immunonephelometry (Behring Nephelometer BNII, Marburg, Germany) with calibration traceable to the International Federation of Clinical Chemistry primary standards [21]. The interassay coefficients of variation of the apoA-I and apoB measurements were 5% and 3%, respectively. HDL-particle number and HDL-size were measured with an automated nuclear magnetic resonance spectroscopic assay as described previously [22]. C-reactive protein (CRP) levels were measured as described previously [23]. Serum PON1-activity toward paraoxon was analyzed as previously described [14]. Serum concentration of myeloperoxidase (MPO) was measured by use of a commercially available ELISA (CardioMPO Test, Prognostix, Cleveland, Ohio). Samples were analyzed in random order to avoid systematic bias. Researchers and laboratory personnel had no access to identifiable information and could identify samples by number only.

PON1 Genotyping

We investigated the two most common single nucleotide polymorphisms (SNP) in the PON1-gene; two substitutions in PON1 [Q192R, rs662, within exon 6 of the PON1 locus, and L55M, rs854560, within exon 3 of the PON1 locus) were analyzed. Genotyping for rs662 and rs854560 was conducted by KBioscience (http://www.kbioscience.co.uk) using KASPar technology. Cases and controls were randomly allocated across DNA study plates with two duplicate samples and two water controls in each 96-well plate. Concordancy between duplicate samples was 100% for both SNPs. Both SNPs were in Hardy-Weinberg equilibrium in control participants (P>0.4).

Statistical Analysis

Baseline characteristics were compared between cases and controls with a mixed-effects model for continuous variables or conditional logistic regression for categorical variables. Triglycerides (non-fasting) and CRP levels had a skewed distribution, and therefore both variables were log-transformed before statistical analysis. Since the PON1-Q192R polymorphism affects serum PON1-activity, but not its concentration, we adjusted PON1-activity for the PON1-Q192R genotype to obtain a proxy for PON1-concentration. Analyses were performed for both (unadjusted) PON1-activity and PON1-activity adjusted for the PON1-Q192R polymorphism.

Associations between PON1-genotype, PON1-activity and cardiovascular risk factors were assessed with Pearson correlations for parametric variables and Spearman correlations for non-parametric variables.

Subsequently, PON1-activity and PON1-activity adjusted for the PON1-Q192R polymorphism were categorized into (sex-specific) quartiles based on the distribution in the controls. Mean levels of cardiovascular risk factors were calculated per quartile. Further, conditional logistic regression was used to calculate odds ratios (ORs) and corresponding 95% confidence intervals (CIs) per quartiles of PON1-activity (±adjustment for PON1-Q192R), as an estimate of the relative risk of CAD with the lowest quartile as the reference category. Conditional logistic regression took into account the matching for sex, age and enrollment time, and ORs were additionally adjusted for alcohol use, apolipoprotein A-I, body mass index (BMI), CRP, diabetes mellitus, fasting time, HDL-cholesterol, HDL-particle number, HDL-size, LDL-cholesterol, MPO, smoking, systolic blood pressure, plasma vitamin-C, vitamin supplement use, waist circumference and non-fasting triglycerides. ORs were also calculated after additional adjustment for PON1-L55M genotype. Finally, we performed stepwise (backward) cox regression analysis for PON1-activity, PON1-genotype as well as PON1-activity adjusted for PON1-genotype to predict the best model for risk of CAD. Statistical analyses were performed with SPSS software (version 12.0.1).

Supporting Information

Table S1.

Sex-specific characteristics of study participants.

https://doi.org/10.1371/journal.pone.0006809.s001

(0.08 MB DOC)

Table S2.

Sex-specific associations of PON1 activity, PON1 genotype and risk factors.

https://doi.org/10.1371/journal.pone.0006809.s002

(0.07 MB DOC)

Table S3.

Sex-specific odds ratios of future coronary artery disease by quartiles of PON1 activity and by quartiles of PON-1 activity adjusted for PON1-192 genotype.

https://doi.org/10.1371/journal.pone.0006809.s003

(0.06 MB DOC)

Table S4.

Backward Stepwise Cox Regression Analysis: HDL-cholesterol with all and excluded variables.

https://doi.org/10.1371/journal.pone.0006809.s004

(0.08 MB DOC)

Table S5.

Backward Stepwise Cox Regression Analysis: HDL-particles with all and excluded variables.

https://doi.org/10.1371/journal.pone.0006809.s005

(0.08 MB DOC)

Table S6.

Backward Stepwise Cox Regression Analysis: apolipoproteinA-I with all and excluded variables.

https://doi.org/10.1371/journal.pone.0006809.s006

(0.08 MB DOC)

Table S7.

Distribution of combined PON1-haplotypes and its effect on HDL-cholesterol, HDL-particles, PON1-acitivity and PON1-activity adjusted for PON1-192 genotype.

https://doi.org/10.1371/journal.pone.0006809.s007

(0.04 MB DOC)

Acknowledgments

We thank the participants, general practitioners, and staff in EPIC-Norfolk. Lipoprotein NMR spectroscopy measurements were kindly performed by Liposcience Inc.

Author Contributions

Conceived and designed the experiments: RSB ESGS MSS NJW JJPK KTK SMB. Performed the experiments: RSB SLR SMB. Analyzed the data: RSB MV MWT SMB. Contributed reagents/materials/analysis tools: RSB MSS NJW KTK SMB. Wrote the paper: RSB MV ESGS JJPK KTK SMB.

References

  1. 1. Tsimikas S, Brilakis ES, Miller ER, McConnell JP, Lennon RJ, et al. (2005) Oxidized phospholipids, Lp(a) lipoprotein, and coronary artery disease. N Engl J Med 353: 46–57.
  2. 2. Diaz MN, Frei B, Vita JA, Keaney JF Jr (1997) Antioxidants and atherosclerotic heart disease. N Engl J Med 337: 408–416.
  3. 3. Witztum JL (1994) The oxidation hypothesis of atherosclerosis. Lancet 344: 793–795.
  4. 4. Frei B, Keaney JF Jr, Retsky KL, Chen K (1996) Vitamins C and E and LDL oxidation. Vitam Horm 52: 1–34.
  5. 5. Li D, Devaraj S, Fuller C, Bucala R, Jialal I (1996) Effect of alpha-tocopherol on LDL oxidation and glycation: in vitro and in vivo studies. J Lipid Res 37: 1978–1986.
  6. 6. Heart Protection Study Collaborative Group (2002) MRC/BHF Heart Protection Study of antioxidant vitamin supplementation in 20,536 high-risk individuals: a randomised placebo-controlled trial. Lancet 360: 23–33.
  7. 7. Lonn E, Bosch J, Yusuf S, Sheridan P, Pogue J, et al. (2005) HOPE and HOPE-TOO Trial Investigators. Effects of long-term vitamin E supplementation on cardiovascular events and cancer: a randomized controlled trial. JAMA 293: 1338–1347.
  8. 8. Lonn E, Yusuf S, Arnold MJ, Sheridan P, Pogue J, et al. (2006) Heart Outcomes Prevention Evaluation (HOPE) 2 Investigators. Homocysteine lowering with folic acid and B vitamins in vascular disease. N Engl J Med 354: 1567–1577. Erratum in:N Engl J Med 355 : 746.
  9. 9. Thomson MJ, Puntmann V, Kaski JC (2007) Atherosclerosis and oxidant stress: the end of the road for antioxidant vitamin treatment? Cardiovasc Drugs Ther 21: 195–210.
  10. 10. Watson AD, Berliner JA, Hama SY, La Du BN, Faull KF, et al. (1995) Protective effect of high density lipoprotein associated paraoxonase. Inhibition of the biological activity of minimally oxidized low density lipoprotein. J Clin Invest 96: 2882–2891.
  11. 11. Mackness MI, Arrol S, Durrington PN (1991) Paraoxonase prevents accumulation of lipoperoxides in low-density lipoprotein. FEBS Lett 286: 152–154.
  12. 12. Aviram M, Rosenblat M, Bisgaier CL, Newton RS, Primo-Parmo SL, et al. (1998) Paraoxonase inhibits high-density lipoprotein oxidation and preserves its functions. A possible peroxidative role for paraoxonase. J Clin Invest 101: 1581–1590.
  13. 13. Jarvik GP, Rozek LS, Brophy VH, Hatsukami TS, Richter RJ, et al. (2000) Paraoxonase (PON1) phenotype is a better predictor of vascular disease than is PON1(192) or PON1(55) genotype. Arterioscler Thromb Vasc Biol 20: 2441–2447.
  14. 14. Mackness B, Durrington P, McElduff P, Yarnell J, Azam N, et al. (2003) Low paraoxonase activity predicts coronary events in the Caerphilly Prospective Study. Circulation 107: 2775–2779.
  15. 15. Wheeler JG, Keavney BD, Watkins H, Collins R, Danesh J (2004) Four paraoxonase gene polymorphisms in 11212 cases of coronary heart disease and 12786 controls: meta-analysis of 43 studies. Lancet 363: 689–695.
  16. 16. Granér M, James RW, Kahri J, Nieminen MS, Syvänne M, et al. (2006) Association of paraoxonase-1 activity and concentration with angiographic severity and extent of coronary artery disease. J Am Coll Cardiol 20: 2429–2435.
  17. 17. Thyagarajan B, Jacobs DR Jr, Carr JJ, Alozie O, Steffes MW, et al. (2008) Factors associated with paraoxonase genotypes and activity in a diverse, young, healthy population: the Coronary Artery Risk Development in Young Adults (CARDIA) study. Clin Chem 54: 738–746.
  18. 18. van Himbergen TM, van der Schouw YT, Voorbij HA, van Tits LJ, Stalenhoef AF, et al. (2008) Paraoxonase (PON1) and the risk for coronary heart disease and myocardial infarction in a general population of Dutch women. Atherosclerosis 199: 408–414.
  19. 19. Davey Smith G, Ebrahim S (2008) ‘Mendelian randomization’: can genetic epidemiology help redress the failures of observational epidemiology? Int J Epidemiol 37: 414–416.
  20. 20. Friedewald WT, Levy RI, Fredrickson DS (1972) Estimation of the concentration of LDL cholesterol in plasma, without use of the preparative ultracentrifuge. Clin Chem 18: 499–502.
  21. 21. Albers JJ, Marcovina SM, Kennedy H (1992) International Federation of Clinical Chemistry standardization project for measurements of apolipoproteins A-I and B. II. Evaluation and selection of candidate reference materials. Clin Chem 38: 658–662.
  22. 22. Jeyarajah EJ, Cromwell WC, Otvos JD (2006) Lipoprotein particle analysis by nuclear magnetic resonance spectroscopy. Clin Lab Med 26: 847–870.
  23. 23. Bruins P, Velthuis Ht, Yazdanbakhsh AP, Jansen PGM, van Hardevelt FWJ, et al. (1997) Activation of the Complement System During and After Cardiopulmonary Bypass Surgery : Postsurgery Activation Involves C-Reactive Protein and Is Associated With Postoperative Arrhythmia. Circulation 96: 3542–3548.
  24. 24. Gordon T, Castelli WP, Hjortland MC, Kannel WB, Dawber TR (1977) High density lipoprotein as a protective factor against coronary heart disease. The Framingham Study. Am J Med 62: 707–714.
  25. 25. Barter P (2004) HDL: a recipe for longevity. Atheroscler Suppl 5: 25–31.
  26. 26. Mackness B, Mackness MI, Arrol S, Turkie W, Julier K, et al. (1998) Serum paraoxonase (PON1) 55 and 192 polymorphism and paraoxonase activity and concentration in non-insulin dependent diabetes mellitus. Atherosclerosis 139: 341–349.
  27. 27. Blatter Garin MC, Abbott C, Messmer S, Mackness M, Durrington P, et al. (1994) Quantification of human serum paraoxonase by enzyme-linked immunoassay: population differences in protein concentrations. Biochem J 304: 549–554.
  28. 28. Bhattacharyya T, Nicholls SJ, Topol EJ, Zhang R, Yang X, et al. (2008) Relationship of paraoxonase 1 (PON1) gene polymorphisms and functional activity with systemic oxidative stress and cardiovascular risk. JAMA 299: 1265–1276.
  29. 29. Odawara M, Tachi Y, Yamashita K (1997) Paraoxonase polymorphism (Gln192-Arg) is associated with coronary heart disease in Japanese noninsulin-dependent diabetes mellitus. J Clin Endocrinol Metab 82: 2257–2260.
  30. 30. Leus FR, Wittekoek ME, Prins J, Kastelein JJ, Voorbij HA (2000) Paraoxonase gene polymorphisms are associated with carotid arterial wall thickness in subjects with familial hypercholesterolemia. Atherosclerosis 149: 371–377.
  31. 31. Watzinger N, Schmidt H, Schumacher M, Schmidt R, Eber B, et al. (2002) Human paraoxonase 1 gene polymorphisms and the risk of coronary heart disease: a community-based study. Cardiology 98: 116–122.
  32. 32. Wang X, Fan Z, Huang J, Su S, Yu Q, et al. (2003) Extensive association analysis between polymorphisms of PON gene cluster with coronary heart disease in Chinese Han population. Arterioscler Thromb Vasc Biol 23: 328–334.
  33. 33. Fortunato G, Rubba P, Panico S, Trono D, Tinto N, et al. (2003) A paraoxonase gene polymorphism, PON 1 (55), as an independent risk factor for increased carotid intima-media thickness in middle-aged women. Atherosclerosis 167: 141–148.
  34. 34. Oliveira SA, Mansur AP, Ribeiro CC, Ramires JA, Annichino-Bizzacchi JM (2004) PON1 M/L55 mutation protects high-risk patients against coronary artery disease. Int J Cardiol 94: 73–77.
  35. 35. Cascorbi I, Laule M, Mrozikiewicz PM, Mrozikiewicz A, Andel C, et al. (1999) Mutations in the human paraoxonase 1 gene: frequencies, allelic linkages, and association with coronary artery disease. Pharmacogenetics 9: 755–761.
  36. 36. Herrmann SM, Blanc H, Poirier O, Arveiler D, Luc G, et al. (1996) The Gln/Arg polymorphism of human paraoxonase (PON 192) is not related to myocardial infarction in the ECTIM Study. Atherosclerosis 126: 299–303.
  37. 37. Gardemann A, Philipp M, Hess K, Katz N, Tillmanns H, et al. (2000) The paraoxonase Leu-Met54 and Gln-Arg191 gene polymorphisms are not associated with the risk of coronary heart disease. Atherosclerosis 152: 421–431.
  38. 38. Ombres D, Pannitteri G, Montali A, Candeloro A, Seccareccia F, et al. (1998) The gln-Arg192 polymorphism of human paraoxonase gene is not associated with coronary artery disease in Italian patients. Arterioscler Thromb Vasc Biol 18: 1611–1616.
  39. 39. Troughton JA, Woodside JV, Yarnell JW, Arveiler D, Amouyel P, Ferrières J, et al. (2008) Paraoxonase activity and coronary heart disease risk in healthy middle-aged males: The PRIME study. Atherosclerosis 197: 556–563.
  40. 40. Samani NJ, Erdmann J, Hall AS, et al. (2007) Genomewide association analysis of coronary artery disease. N Engl J Med 357: 443–453.
  41. 41. Rosenzweig A (2007) Scanning the genome for coronary risk. N Engl J Med 357: 497–499.
  42. 42. Jarvik GP, Tsai NT, McKinstry LA, et al. (2002) Vitamin C and E intake is associated with increased paraoxonase activity. Arterioscler Thromb Vasc Biol 22: 1329–1333.
  43. 43. Ascherio A, Rimm EB, Hernán MA, et al. (1999) Relation of consumption of vitamin E, vitamin C, and carotenoids to risk for stroke among men in the United States. Ann Intern Med 130: 963–970.
  44. 44. Boekholdt SM, Peters RJG, Day NE, Luben R, Bingham SA, et al. (2004) Macrophage migration inhibitory factor and the risk of myocardial infarction or death due to coronary artery disease in adults without prior myocardial infarction or stroke: The EPIC-Norfolk Prospective Population study. The American Journal of Medicine 117: 390–397.