ORIGINAL CONTRIBUTION Relationship of Paraoxonase 1 (PON1) Gene Polymorphisms and Functional Activity With Systemic Oxidative Stress and Cardiovascular Risk Tamali Bhattacharyya, MD, MS Stephen J. Nicholls, MBBS, PhD Eric J. Topol, MD Renliang Zhang, MD, PhD Xia Yang, PhD David Schmitt, BA Xiaoming Fu, MS Mingyuan Shao, MS Danielle M. Brennan, MS Stephen G. Ellis, MD Marie-Luise Brennan, PhD Hooman Allayee, PhD Aldons J. Lusis, PhD Stanley L. Hazen, MD, PhD S YSTEMIC OXIDATIVE STRESS PRO- motes a number of key events in the development of atherosclerosis, and the search to identify effective therapeutic antioxidant strategies is of potential clinical interest. For example, in addition to promoting reverse cholesterol transport, highdensity lipoprotein (HDL) is also reported to promote systemic antiinflammatory and antioxidant effects. The HDL associated esterase/lactonase paraoxonase 1 (PON1) is implicated in contributing to the anti-inflammatory and antioxidant activities of the lipoprotein,1 although direct demonstration of this in humans is lacking. Several lines of evidence suggest that PON1 promotes antioxidant and atheroprotective effects.2,3 PON1 inhibits oxidation of low-density lipoprotein (LDL) in vitro.2,4 Genetic deletion of PON1 is associated with increased sus- Context Paraoxonase 1 (PON1) is reported to have antioxidant and cardioprotective properties. The relationship between PON1 genotypes and functional activity with systemic measures of oxidative stress and cardiovascular disease (CVD) risk in humans has not been systematically investigated. Objective To investigate the relationship of genetic and biochemical determinants of PON1 activity with systemic measures of oxidative stress and CVD risk in humans. Design, Setting, and Participants The association between systemic PON1 activity measures and a functional polymorphism (Q192R) resulting in high PON1 activity with prevalent CVD and future major adverse cardiac events (myocardial infarction, stroke, or death) was evaluated in 1399 sequential consenting patients undergoing diagnostic coronary angiography between September 2002 and November 2003 at the Cleveland Clinic. Patients were followed up until December 2006. Systemic levels of multiple structurally defined fatty acid oxidation products were also measured by mass spectrometry in 150 age-, sex-, and race-matched patients and compared with regard to PON1 genotype and activity. Main Outcome Measures Relationship between a functional PON1 polymorphism and PON1 activity with global indices of systemic oxidative stress and risk of CVD. Results The PON1 genotype demonstrated significant dose-dependent associations (QQ192⬎QR192⬎RR192) with decreased levels of serum PON1 activity and with increased levels of systemic indices of oxidative stress. Compared with participants with either the PON1 RR192 or QR192 genotype, participants with the QQ192 genotype demonstrated an increased risk of all-cause mortality (43/681 deaths [6.75%] in RR192 and QR192 and 62/584 deaths [11.1%] in QQ192; adjusted hazard ratio, 2.05; 95% confidence interval [CI], 1.32-3.18) and of major adverse cardiac events (88/681 events [13.6%] in RR192 and QR192 and 102/584 events [18.0%] in QQ192; adjusted hazard ratio, 1.48; 95% CI, 1.09-2.03; P=.01). The incidence of major adverse cardiac events was significantly lower in participants in the highest PON1 activity quartile (23/315 [7.3%]) and 235/324 [7.7%] for paraoxonase and arylesterase, respectively) compared with those in the lowest activity quartile (78/311 [25.1%] and 75/319 [23.5%]; P⬍.001 for paraoxonase and arylesterase, respectively). The adjusted hazard ratios for major adverse cardiac events between the highest and lowest PON1 activity quartiles were, for paraoxonase, 3.4 (95% CI, 2.1-5.5; P⬍.001) and for arylesterase, 2.9 (95% CI, 1.8-4.7; P⬍.001) and remained independent in multivariate analysis. Conclusion This study provides direct evidence for a mechanistic link between genetic determinants and activity of PON1 with systemic oxidative stress and prospective cardiovascular risk, indicating a potential mechanism for the atheroprotective function of PON1. www.jama.com JAMA. 2008;299(11):1265-1276 ceptibility of LDL to oxidation ex vivo, increased measures of macrophage oxidative stress, and increased lesion size in animal models of atherosclerosis; ©2008 American Medical Association. All rights reserved. Author Affiliations are listed at the end of this article. Corresponding Author: Stanley L. Hazen, MD, PhD, Center for Cardiovascular Diagnostics and Prevention, Cleveland Clinic, 9500 Euclid Ave, NE-10, Cleveland, OH 44195 ([email protected]). (Reprinted) JAMA, March 19, 2008—Vol 299, No. 11 Downloaded from jama.ama-assn.org at University of California - Los Angeles on January 6, 2012 1265 RELATIONSHIP OF PARAOXONASE 1 WITH SYSTEMIC OXIDATIVE STRESS AND CVD RISK conversely, overexpression of the human PON1 transgene in mice results in reduced aortic lesion size and corresponding decreases in epitopes recognized by antibodies specific for oxidized lipid-protein adducts.5-8 Despite evidence that PON1 prevents atherosclerosis in animal models, it remains to be established whether PON1 possesses atheroprotective and antioxidant properties in humans. Several studies have suggested that PON1 may play an atheroprotective role but the simultaneous associations between PON1 polymorphisms and enzyme activity with cardiovascular disease (CVD) risk have been reported in only some population studies.1-3,9-13 Further confounding the relationship in humans has been the discovery of multiple PON1 polymorphisms in both the coding region of the protein and the promoter, some of which reportedly influence overall systemic activity levels.14 In a recent meta-analysis of 43 studies examining multiple single-nucleotide polymorphisms (SNPs) for PON1, the most promising of the genotypes was the Q192R variant, although this SNP demonstrated only a weak overall association with coronary heart disease of uncertain relevance.15 However, each of the contributing studies in the metaanalyses involved small cohorts of patients, with variable evaluation of functional activity and no quantitative examination of systemic indices of oxidative stress. The major conclusions of the meta-analyses therefore emphasized the need for both much larger and mechanistic investigations of the role of PON1 in human CVD, particularly since no significant association for the Q192R polymorphism was noted among the larger studies in the meta-analyses, and accompanying simultaneous functional and genetic studies were typically lacking.15 Thus, it remains to be established whether genetic and biochemical determinants of PON1 are linked to oxidative stress and CVD risk in humans. The Q192R polymorphism involves a mutation from glutamine (Q, wild type) to arginine (R, variant) at amino acid po- sition 192 of the protein sequence. Functional PON1 activity can be measured by its ability to hydrolyze exogenous substrates such as paraoxon and phenylacetate, reflecting so-called paraoxonase and arylesterase activity, respectively. Functional differences have been observed in hydrolysis rates of the Q192 vs R192 alloenzymes using paraoxon as substrate, although no difference has been reported regarding their ability to hydrolyze phenylacetate.16 In this present large, prospective clinical study, we report a comprehensive in vivo investigation of a mechanistic link between the functional PON1 Q192R polymorphism, serum PON1 activity using dual enzyme activity measurements to account for the differential rate of hydrolysis of the alloenzymes, multiple indices of systemic oxidative stress, and risk of both prevalent atherosclerotic CVD and nearterm incident major adverse cardiovascular events (myocardial infarction [MI], stroke, and death). METHODS Study Design and Sample Collection PON1 activity and functional polymorphisms were determined in serum and DNA samples of 1399 sequential consenting patients who participated in the GeneBank study between September 2002 and November 2003. GeneBank is a single-site (Cleveland Clinic, Cleveland, Ohio) sample repository generated from patients undergoing elective diagnostic coronary angiography with extensive clinical and laboratory characterization and longitudinal observation. For systemic measures of oxidative stress, whole blood collected in EDTA tubes was immediately spun, plasma and buffy coat were isolated, and plasma was stored under argon atmosphere with antioxidant cocktail supplement, as previously described.17 Serum for PON1 activity measures was obtained from serum separator tubes after 30 to 60 minutes of clotting time at room temperature. All specimens were stored at −80°C until time of analysis. Pa- 1266 JAMA, March 19, 2008—Vol 299, No. 11 (Reprinted) tients were followed up on an annual basis for adjudicated incident major adverse cardiac events and mortality until December 14, 2006. The GeneBank study was approved by the institutional review board of the Cleveland Clinic. All patients provided written informed consent prior to being enrolled in the study. Clinical Diagnosis and Definition of Outcomes Information regarding demographics, medical history, and medication use was obtained by patient interview and confirmed by chart review. Race information used in analyses was prespecified prior to the study and was based on self-report. All clinical outcomes data were verified by source documentation. Mortality was assessed using the Social Security Death Index.18 Cardiovascular disease was defined by the presence of coronary artery disease or peripheral arterial disease. Coronary artery disease included adjudicated diagnoses of stable or unstable angina, MI (adjudicated definition based on defined electrocardiographic changes or elevated cardiac enzymes), or angiographic evidence of at least 50% stenosis of 1 or more epicardial vessels. Peripheral artery disease was defined as the presence of any extracoronary atherosclerosis and included obstructive disease (including a history of intermittent claudication), amaurosis fugax, history of cerebrovascular accident, or evidence of either arterial stenosis or aneurysmal disease in the thoracic limbs or abdominal aorta on Doppler ultrasound. Prospective cardiovascular risk was assessed by the incidence of major adverse cardiovascular events (MACE), which included nonfatal and fatal MI, nonfatal and fatal stroke, and allcause mortality. Nonfatal events were defined as MI or stroke in patients who survived at least for 48 hours following the onset of symptoms. We also assessed the risk of true incidence or the first cardiovascular event in participants without a history of CVD at enrollment (baseline) and the risk for re- ©2008 American Medical Association. All rights reserved. Downloaded from jama.ama-assn.org at University of California - Los Angeles on January 6, 2012 RELATIONSHIP OF PARAOXONASE 1 WITH SYSTEMIC OXIDATIVE STRESS AND CVD RISK current cardiovascular events in participants with an established diagnosis of CVD at baseline according to levels of PON1 activity. Determination of PON1 Activity Serum arylesterase and paraoxonase activities were independently measured by UV spectrophotometry in a 96-well plate format (Spectramax 384 Plus, Molecular Devices, Sunnyvale, California) using phenyl acetate or paraoxon (SigmaAldrich, St Louis, Missouri) as substrates, respectively.19 Briefly, for arylesterase assays, initial hydrolysis rates were determined at 270 nm in 50-fold diluted serum (final) in reactions mixtures composed of 3.4mM phenylacetate, 9mM Tris hydrocholoride, pH 8, and 0.9mM calcium chloride at 24°C. An extinction coefficient (at 270 nm) of 1310M−1 · cm−1 was used for calculating units of arylesterase activity, which are expressed as the amount of phenyl acetate hydrolyzed in micromoles per minute per milliliter of serum. For paraoxonase activity assays, rate of generation of para-nitrophenol was determined at 405 nm in 40-fold diluted serum (final) in reaction mixtures composed of 1.5mM paraoxon, 10mM Tris hydrocholoride, pH 8, 1M sodium chloride, and 2mM calcium chloride at 24°C. An extinction coefficient (at 405 nm) of 17 000 M−1 · cm−1 was used for calculating units of paraoxonase activity, which are expressed as the amount of para-nitrophenol produced in nanomoles per minute per milliliter of serum. Paraoxonase and arylesterase assays for each sample were performed in duplicates, with average measurements of enzyme activity for each sample calculated. Each 96-well plate included blank samples to monitor spontaneous hydrolysis of substrates and aliquots of serum samples of 3 pooled calibrators (low, mid, and high levels) with known activity levels to ensure assay quality of each plate based on established acceptability criterion. The intra-assay and interassay coefficients of variance for performance of arylesterase were 1.2% and 3.9%, respectively, and the intra-assay and interas- say coefficients of variance for performance of paraoxonase activity assays were 2.0% and 5.6%, respectively, on 20 replicates performed on 10 different days. To establish normal ranges, serum arylesterase and paraoxonase activities were also determined on 100 apparently healthy volunteers (50 men and 50 women) aged 55 years or older (mean, 64 [SD, 4] years) responding to local advertisements. Systemic arylesterase activity in this healthy, middleaged population ranged from 169.3 to 814.0 µmol/min/mL, with median levels of 605.8 µmol/min/mL (interquartile range [IQR], 517.6-666.8 µmol/ min/mL) of serum. Systemic paraoxonase activity in this healthy, middle-aged population ranged from 436.4 to 4025 nmol/min/mL, with median levels of 1264 nmol/min/mL (IQR, 647.0-1865 nmol/min/mL) of serum. Genotyping Of the 1399 participants, 1386 DNA samples from the GeneBank cohort were available for genotyping for the PON1 Q192R polymorphism (SNP rs662). Primers for amplifying the sequences containing the SNP were designed using Primer 3 (http://www.genome.wi.mit .edu/cgi-bin/primer/primer3_www .cgi). Sense and antisense probes used for fluorescence polarization–single base extension detection of the SNP were designed using Primer PREMIER (PREMIER Biosoft International, Palo Alto, California). Each genotyping reaction consisted of a single-plex polymerase chain reaction followed by a single basepair extension reaction using either the sense or antisense probe and dideoxynucleotides labeled with TAMRA or R110 for the alternative alleles. The fluorescence signal was detected by fluorescence polarization using the analyst HT (Molecular Devices) and genotypes were determined based on the plot of TAMRA vs R110 signal values. Mass Spectrometry Assays Plasma levels of structurally specific species of hydroxyeicosatetraenoic acids (HETEs), hydroxyoctadecadienoic acids (HODEs), the 8-isopros- ©2008 American Medical Association. All rights reserved. tane prostaglandin F2␣ (8-isoPGF2␣), and their precursor fatty acids (arachidonic or linoleic acids) were quantified in samples of participants using stable isotope dilution high-performance liquid chromatography with online electrospray ionization tandem mass spectrometry. Prostaglandin F2␣-d4 (Cayman Chemicals, Ann Arbor, Michigan) was used as internal standard for calibration of 8-isoPGF 2␣ and 15HETE-d8 (Cayman Chemicals) for other oxidized fatty acids. Quantification of total plasma levels of each analyte (free plus esterified) were performed following addition of the appropriate isotopically labeled internal standards and base-catalyzed hydrolysis, as previously described.17 Analyses were performed on plasma specimens collected from 50 age-, sex-, and race-matched participants of each genotype (QQ192, QR192, and RR192) who were randomly selected from among the 1386 patients. To avoid bias, 50 participants with the genotype of the lowest frequency in our cohort (RR192) were initially randomly selected among the 136 participants who carried that genotype (50% male; 100% white; mean age, 63 [SD, 2] years). An additional 50 participants from each of the QQ192 and QR192 cohorts, matched with the QQ192 participants for age, sex, and race, were also selected for mass spectrometry analyses. Statistical Analyses Clinical diagnosis, outcome definition, determination of PON1 activity, genotyping, and mass spectrometry analyses were each performed by investigators who were blinded to CVD status and other measurements. Continuous variables are presented as mean (SD) or median (IQR) for non–normally distributed data and categorical variables as numbers and percentages. Regression analysis was used for calculating the variance (R2) explained by the PON1 polymorphism and activity measures. Quartile cut points were determined from the enzyme activity levels of all 1399 study participants at baseline. Analysis of variance or the Kruskal-Wallis test (for non– (Reprinted) JAMA, March 19, 2008—Vol 299, No. 11 Downloaded from jama.ama-assn.org at University of California - Los Angeles on January 6, 2012 1267 RELATIONSHIP OF PARAOXONASE 1 WITH SYSTEMIC OXIDATIVE STRESS AND CVD RISK normally distributed data) was used to test the difference in mean oxidized fatty acid levels according to PON1 genotype and tertiles of PON1 activity. Kaplan-Meier methods were used to plot time-to-event curves for PON1 activity quartiles and genotypes and the logrank test was performed to assess differences between curves. Logistic regression analysis was used to calculate adjusted odds ratios (ORs) for the association between activity levels and prevalent CVD status after adjusting for the traditional risk factors and selected classes of medications (statins and aspirin) as described herein. Cox proportional hazard models were performed to determine if PON1 is an independent predictor of future cardiac events. Unadjusted hazard ratios (HRs) for clinical events were calculated with reference to the highest quartile (corresponding to lowest risk). The models were adjusted for all traditional cardiac risk factors, including the Framingham ATP-III risk score (including diabetes status), log Creactive protein, body mass index, and medication use (statins and aspirin). Models were created separately for paraoxonase and arylesterase. All variables used in the models met the proportional hazards assumption by testing them as time-dependent covariates in the multivariate model. Receiver operating characteristic curves were plotted to estimate the C index for MACE with and without PON1 activity as predictors.20 The C index is analogous to the area under the receiver operating characteristic curve but takes into account right censoring.21 To evaluate the contribution of PON1 activity as a predictive marker, we calculated the concordance indices22 with and without each PON1 activity measurement in separate multivariate models that included the variables described herein. The improvement in predictability with the addition of each of these wellestablished risk factors was assessed by the difference in the concordance indices.23 The difference of the indices was bias-corrected and bootstrapping was used to generate 95% confidence inter- vals (CIs). A 1-sample t test was performed to determine if the difference was equal to zero. Hazard ratios for clinical events were calculated for participants with the PON1 Q192R, polymorphism and were adjusted after controlling for differences in the aforementioned traditional cardiac risk factors and medications. An additive model of inheritance was also created to assess the change in risk for having 1 and 2 copies of the Q allele. The genotypes were given values of 0, 1, and 2 and entered in a Cox proportional hazards model. All statistical analyses were performed using SPSS, version 11 (SPSS Inc, Chicago, Illinois) and verified on SAS, version 8.2 (SAS Institute Inc, Cary, North Carolina). All P values are 2-sided, with P⬍ .05 considered significant. RESULTS Clinical, Laboratory, and Demographic Characteristics The clinical characteristics of participants stratified according to a diagnosis of CVD at time of enrollment are summarized in TABLE 1. Participants with CVD were older, with a greater prevalence of history of hyperlipidemia, diabetes, and hypertension and greater use of aspirin and statins. While 66% of patients were taking statins, less than 10% of the cohort were taking fenofibrates. Participants with CVD had lower levels of HDL cholesterol (HDL-C) and LDL cholesterol (LDL-C) (the latter presumably because of increased statin use) and higher levels of triglycerides and C-reactive protein. In the entire cohort, 46.3% (642/1386) had the QQ192 genotype, 43.9% (608/ 1386)hadtheQR192genotype,and9.8% (136/1386) had the RR192 genotype. Whiteparticipantsdemonstratedagreater frequency of the QQ192 genotype (621/ 1254 [49.5%] vs 17/126 [13.5%]; P⬍.001) and lower frequency of the RR192 genotype (92/1254 [7.4%] vs 44/ 126 [34.9%]; P⬍.001) compared with nonwhite participants. Irrespective of race, paraoxonase and arylesterase activity levels were significantly lower in participants with CVD in 1268 JAMA, March 19, 2008—Vol 299, No. 11 (Reprinted) comparison with non-CVD patients: for paraoxonase, the median was 860.6 (IQR, 442.6-1599) nmol/min/mL vs 1164 (IQR, 473.5-1825) nmol/min/mL (P⬍.001), and for arylesterase, the median was 334.6 (IQR, 279.5-395.8) µmol/ min/mL vs 356.9 (IQR, 296-424.4) µmol/ min/mL (P⬍.001) (Table 1). Similarly, systemic levels of serum paraoxonase and arylesterase activity were lower in participants with CVD compared with healthy middle-aged volunteers: for paraoxonase, median, 860.6 (IQR, 442.61599) nmol/min/mL vs 1264 (IQR, 647.0-1864.6) nmol/min/mL (P⬍.001), and for arylesterase, 334.6 (IQR, 279.5395.8) µmol/min/mL vs 605.8 (IQR, 517.6-666.8) µmol/min/mL (P⬍.001). During follow-up (mean, 44 [SD, 7] months; minimum of 3 years), 13.8% of participants (193/1399) experienced at least 1 MACE (death, 7.6% [106/1399]; nonfatal MI, 5.7% [80/ 1399]; and stroke, 1.1% [16/1399]). PON1 Genotype and Paraoxonase Activity The relationships between PON1 Q192R genotypes and activity measures throughout the entire population are summarized in (TABLE 2). The PON1 Q192R polymorphism genotypes were in Hardy-Weinberg equilibrium. Participants with the QQ192 genotype demonstrated significantly less paraoxonase activity than participants with either the QR192 genotype (median, 453.5 [IQR, 360.4-560.3] nmol/min/mL vs 1436 [IQR, 11291800] nmol/min/mL; P ⬍.001) or the RR192 genotype (median, 453.5 [IQR, 360.4-560.3] nmol/min/mL vs 2374 [IQR, 1978-2851] nmol/min/mL; P ⬍.001) (Table 2). While 49.4% of participants with the wild-type QQ192 genotype (317/642) were in the lowest quartile of PON1 activity, only 4.4% of participants with the variant RR192 genotype (6/136) were in the lowest activity quartile. In contrast, 3.9% of participants with the QQ192 genotype (25/642) and 83.8% with the RR192 genotype (114/136) were in the highest activity quartile (Table 2). ©2008 American Medical Association. All rights reserved. Downloaded from jama.ama-assn.org at University of California - Los Angeles on January 6, 2012 RELATIONSHIP OF PARAOXONASE 1 WITH SYSTEMIC OXIDATIVE STRESS AND CVD RISK In the lowest activity quartile, 92.2% of study participants (317/344) had the QQ192 genotype and only 1.7% (6/ 344) had the RR192 genotype. The reverse trend was observed in the highest activity quartile, with 7.2% of participants (25/348) having the QQ192 genotype and 32.8% (114/ 348) having the RR192 genotype (TABLE 3). Regression analysis confirmed that the PON1 Q192R polymorphism accounted for 58.5% (R2 =0.585; P⬍.001) of the variation in serum paraoxonase activity levels throughout the population. PON1 Genotype, Paraoxonase Activity, and Systemic Oxidative Stress Table 1. Baseline Characteristics and Prospective Events Among Participants Undergoing Elective Diagnostic Cardiac Catheterization (N = 1399) a With CVD at Baseline (n = 1116) Characteristics Race White P Value 65.1 (10.9) 57.2 (11.8) ⬍.001 799/1116 (71.6) 135/283 (47.7) ⬍.001 1023/1111 (92.1) 243/282 (86.2) .002 71/1111 (6.4) 38/282 (13.5) ⬍.001 Age, mean (SD), y Male Without CVD at Baseline (n = 283) African American Hyperlipidemia b 940/1093 (86.0) 168/277 (60.6) ⬍.001 Hypertension 843/1100 (76.6) 150/282 (53.2) ⬍.001 Diabetes mellitus c 433/1095 (39.5) 49/271 (18.1) ⬍.001 Current smoking 163/1116 (14.6) 33/283 (11.7) .20 Statin use 717/1088 (65.9) 80/270 (29.6) ⬍.001 Aspirin use 873/1098 (79.5) 176/272 (64.7) ⬍.001 Body mass index d Blood pressure, mean (SD), mm Hg Systolic Diastolic 29.7 (5.9) 29.7 (6.3) .58 134.3 (21.2) 134.8 (21.3) .96 74.0 (13.2) 75.5 (12.2) .20 46.1 (13.0) 53.9 (15.9) ⬍.001 Plasma levels of multiple structurally specific oxidized fatty acids were quantified in participants and analyzed for their relationships with the PON1 genotype (TABLE 4) and paraoxonase activity (Table 4). Participants with the RR192 genotype had lower levels of all measured systemic indices of oxidative stress compared with age-, sex-, and race-matched participants possessing the QQ192 genotype (P ⬍.001 for all comparisons). Consistent with these findings, serum paraoxonase activity levels were inversely correlated with multiple direct systemic indices of oxidative stress in a dose-dependent fashion (Table 4). Collectively, these results provide both genetic and biochemical support for the notion that the PON1 Q192R variant strongly influences quantitative measures of systemic oxidative stress in humans. Prospective event rates, No. (Kaplan-Meier %) [95% CI] MI Association of the PON1 Q192R Polymorphism With Prevalent CVD and CVD Outcomes Table 2. Distribution of Paraoxonase Activity Quartiles in Each PON1 Q192R Genotype a An increased prevalence of coronary artery disease was observed in participants withthePON1QQ192genotype(461/962 [47.9%] with vs 169/405 [41.7%] without coronary artery disease; P=.04). In contrast, participants with the RR192 genotype showed the opposite tendency with lower prevalence of coronary artery disease (TABLE 5). Those with the PON1 QQ192genotypeweresimilarlyobserved HDL-C, mean (SD), mg/dL 112.4 (36.6) ⬍.001 Triglycerides, median (IQR), mg/dL 138.0 (99.0-202.0) 115.5 (77.0-173.0) ⬍.001 e Total cholesterol, mean (SD), mg/dL 178.4 (46.4) 193.9 (43.5) ⬍.001 LDL-C, mean (SD), mg/dL 99.0 (37.5) C-reactive protein, median (IQR), mg/L Framingham ATP-III score, mean (SD) 3.0 (1.5-7.0) .005 e 2.4 (1.2-5.9) 13.3 (3.5) ⬍.001 11.4 (5.0) Paraoxonase, median (IQR), nmol/min/mL 860.6 (442.6-1599.0) 1164.0 (473.5-1825.0) Arylesterase, median (IQR), µmol/min/mL 334.6 (279.5-395.8) 356.9 (296.0-424.4) .001 e ⬍.001 69 (7.2) [5.6-8.9] 11 (4.4) [1.8-6.9] MI/CVA 84 (8.9) [7.1-10.7] 12 (4.7) [2.1-7.4] .14 .05 Death 99 (10.3) [8.3-12.2] 7 (2.9) [0.8-5.0] ⬍.001 MACE 175 (17.9) [15.5-20.3] 18 (7.2) [4.0-10.4] ⬍.001 Abbreviations: CI, confidence interval; CVA, cerebrovascular accident; CVD, cardiovascular disease; HDL-C, highdensity lipoprotein cholesterol; IQR, interquartile range; MACE, major adverse cardiac events (ie, myocardial infarction, stroke, or death); MI, myocardial infarction; LDL-C, low-density lipoprotein cholesterol. SI conversions: To convert HDL-C, LDL-C, and total cholesterol to mmol/L, multiply by 0.0259; to convert triglycerides to mmol/L, multiply by 0.0113; to convert C-reactive protein to nmol/L, multiply by 9.524. a Data are presented as No./total (%) unless otherwise indicated. b Hyperlipidemia was defined as a total cholesterol level of higher than 200 mg/dL or treatment with an antihyperlipidemic agent. c Diabetes mellitus was defined as a fasting plasma glucose level of at least 125 mg/dL or treatment with hypoglycemic agents. d Calculated as weight in kilograms divided by height in meters squared. e By Wilcoxon rank-sum test. PON1 Q192R Genotype Paraoxonase Quartile, nmol/min/mL Median (IQR) QQ192 (n = 642) QR192 (n = 608) RR192 (n = 136) 453.5 (360.4-560.3) 1436 (1129-1800) 2374 (1978-2851) Quartile 4 (⬎1640) 25 (3.9) 209 (34.4) 114 (83.8) Quartile 3 (1640-899.1) 14 (2.2) 321 (52.8) 14 (10.3) Quartile 2 (899-450) 286 (44.5) 57 (9.4) 2 (1.5) Quartile 1 (⬍450) 317 (49.4) 21 (3.5) 6 (4.4) Abbreviations: CVD, cardiovascular disease; IQR, interquartile range; PON1, paraoxonase 1 gene. a Data are presented as No. (%) unless otherwise indicated. Genotypes: QQ192, wild type; QR192, heterozygous; RR192, mutant homozygous. ©2008 American Medical Association. All rights reserved. (Reprinted) JAMA, March 19, 2008—Vol 299, No. 11 Downloaded from jama.ama-assn.org at University of California - Los Angeles on January 6, 2012 1269 RELATIONSHIP OF PARAOXONASE 1 WITH SYSTEMIC OXIDATIVE STRESS AND CVD RISK to have an increased likelihood of having a history of coronary artery bypass graft surgery (P=.03) and a tendency toward increased history of percutaneous coronary intervention (P=.09) (data not shown). In contrast, those with the PON1 RR192 genotype demonstrated opposite tendencies (diminished history of revascularization;P=.09)(datanotshown).IndividualscarryingthePON1QR192genotype showed no observed differences in prevalent peripheral artery disease within the cohort (Table 5). Of the 1386 participants with genotype information, follow-up data were available in 1265 participants (QQ192 genotype, n = 584; QR192 genotype, n = 563; and RR192 genotype, n = 118). The PON1 functional polymorphism Q192R was significantly associated with CVD outcomes within the cohort. For example, all adverse cardiovascular outcomes monitored were observed to a significantly lesser degree in patients who carried 1 R allele. Kaplan-Meier estimates of all-cause mortality and MACE revealed a significant trend in our study participants. Homozygous and heterozygous participants (RR192 and QR192 genotypes) had significantly lower event rates for allcause mortality (43/681 [Kaplan-Meier estimate, 6.75%]) compared with participants with the QQ192 genotype (62/ 584 [Kaplan-Meier estimate, 11.10%]; P = .006). Similarly, the event rate for MACE over the ensuing 3-year period was lower with the RR192 and QR192 genotypes (88/681 [Kaplan-Meier estimate, 13.6%]) vs the QQ192 genotype (102/584 [Kaplan-Meier estimate, 18.0%]; P=.03) (TABLE 6 and FIGURE). In a multivariate model including age, Framingham ATP-III risk score, race, log C-reactive protein, body mass index, and medication (statin and aspirin) use, the addition of the risk allele Q192 showed no association with incident MI and stroke over the ensuing 3 years following enrollment but was significantly associated with in- creased likelihood of death and MACE (Table 6). Compared with participants with either the RR192 or QR192 genotype, participants with the QQ192 genotype demonstrated an adjusted HR of 2.05 (95% CI, 1.32-3.18) for allcause mortality (P = .001). Furthermore, MACE were also more likely to occur over the ensuing 3-year period in participants with the QQ192 genotype compared with participants with either the RR192 or QR192 genotype (adjusted HR, 1.48; 95% CI, 1.092.03; P=.01) (Table 6). In separate analyses, an additive model of inheritance was created to assess the change in risk for having a copy of the Q allele. An adjusted HR for 3-year all-cause mortality of 1.58 (95% CI, 1.12-2.24; P=.01) was observed for having a Q allele. Furthermore, in this model the presence of a Q allele was also associated with an increased risk of having a MACE over the ensuing 3-year period (HR, 1.32; 95% CI, 1.04-1.69; P =.03). Table 3. Distribution of PON1 Q192R Genotypes in Each Paraoxonase Activity Quartile a Paraoxonase Quartile, nmol/min/mL PON1 Q192R Genotype QQ192 QR192 RR192 Quartile 4 (⬎1640) (n = 348) Quartile 3 (1640-899.1) (n = 349) 25 (7.2) 209 (60.1) 114 (32.8) 14 (4.0) 312 (92.0) 14 (4.1) Quartile 2 (899-450) (n = 345) 286 (82.9) 57 (16.5) 2 (0.6) Quartile 1 (⬍450) (n = 344) 317 (92.2) 21 (6.1) 6 (1.7) Abbreviations: CVD, cardiovascular disease; IQR, interquartile range; PON1, paraoxonase 1 gene. a Data are presented as No./total (%). Genotypes: QQ192, wild type; QR192, heterozygous; RR192, mutant homozygous. Table 4. Total Plasma Oxidized Fatty Acid Levels According to PON1 Q192R Genotype and Tertile of Serum Paraoxonase Activity a PON1 Q192R Genotype Oxidized Fatty Acid 5-HETE 8-HETE 9-HETE 11-HETE 12-HETE 15-HETE 9-HODE 13-HODE 8-isoPGF2␣ QQ192 (n = 50) QR192 (n = 50) Paraoxonase Tertile RR192 (n = 50) 20.1 (14.3-24.9) 15.5 (12.6-19.2) 9.1 (7.6-12.8) 3.0 (2.5-4.4) 2.5 (1.8-2.7) 1.5 (1.1-2.1) 42.7 (30.3-63.1) 30.1 (22.0-35.2) 17.2 (15.1-25.3) 5.8 (4.8-7.1) 7.6 (5.9-10.4) 24.3 (20.6-30.8) 35.5 (30.2-43.7) 38.7 (30.4-49.0) 11.1 (4.4-15.6) 4.5 (3.7-6.1) 6.1 (4.9-8.1) 3.7 (3.3-5.1) 4.4 (3.3-6.4) 20.2 (16.2-26.3) 13.1 (10.5-19.5) 28.8 (22.6-37.1) 20.3 (15.5-27.9) 29.6 (23.0-38.9) 22.2 (17.6-28.3) 11.5 (3.0-28.5) 2.8 (1.1-6.7) P Value Tertile 1 (⬍630) Tertile 2 (630.1-1928) Tertile 3 (⬎1928) (n = 50) (n = 50) (n = 50) P Value ⬍.001 ⬍.001 ⬍.001 21.8 (16.0-25.3) 3.3 (2.6-4.6) 45.2 (35.2-66.3) 16.1 (13.1-20.7) 2.5 (2.2-2.7) 30.4 (25.3-37.5) 8.4 (6.0-9.7) 1.2 (0.9-1.7) 15.8 (11.5-17.9) ⬍.001 ⬍.001 ⬍.001 ⬍.001 ⬍.001 5.9 (5.0-8.2) 7.8 (6.2-11.3) 4.6 (4.1-6.3) 6.4 (5.4-8.7) 3.6 (2.5-4.1) 3.8 (2.7-4.8) ⬍.001 ⬍.001 ⬍.001 ⬍.001 ⬍.001 ⬍.001 24.7 (21.2-30.8) 36.3 (30.9-45.7) 39.8 (31.5-52.6) 12.1 (4.5-28.2) 22.0 (17.4-30.7) 33.0 (24.8-38.8) 30.5 (25.5-43.4) 14.2 (7.8-29.4) 11.3 (8.4-14.8) 18.3 (10.1-23.1) 19.3 (10.1-26.0) 1.6 (0.1-5.1) ⬍.001 ⬍.001 ⬍.001 ⬍.001 Abbreviations: HETE, hydroxyeicosatetraenoic acid; HODE, hydroxyloctadecadienoic acid; PON1, paraoxonase 1 gene; QQ192, wild type subjects; QR192, heterozygous subjects; RR192, mutant homozygous subjects; 8-iso PGF2␣, 8-isoprostane prostaglandin F2␣. a Data are presented as median (interquartile range) of oxidized fatty acid in picomoles per milliliter. 1270 JAMA, March 19, 2008—Vol 299, No. 11 (Reprinted) ©2008 American Medical Association. All rights reserved. Downloaded from jama.ama-assn.org at University of California - Los Angeles on January 6, 2012 RELATIONSHIP OF PARAOXONASE 1 WITH SYSTEMIC OXIDATIVE STRESS AND CVD RISK Association of PON1 Activity With Prevalent CVD and CVD Outcomes A higher prevalence of CVD was also observed in patients with low serum activity levels of either paraoxonase or arylesterase. Following adjustment for differences in known CVD risk factors, this association was not significant for paraoxonase (OR, 1.5; 95% CI, 0.97-2.2; P =.07 for comparison between upper and lower quartiles) but was significant for arylesterase (OR, 2.1; 95% CI, 1.4-3.1; P=.001). Low systemic levels of enzyme activity were significantly associated with the presence of coronary artery disease (paraoxonase: OR, 1.5; 95% CI, 1.03-2.3; P = .03; arylesterase: OR, 2.0; 95% CI, 1.4-3.1; P = .001) or in combination with peripheral artery disease (paraoxonase: OR, 1.4; 95% CI, 0.9-2.0; P = .11; arylesterase: OR, 1.2; 95% CI, 1.0-2.11; P =.05). Event rates for all prospective cardiovascular events were significantly lower in participants in the highest PON1 activity quartile compared with participants in the lowest quartile (P⬍.001 for all prospective events comparisons) (TABLE 7 and Figure). The frequency of incident cases of nonfatal MI or stroke among participants within the highest quartiles of PON1 activity was 2.5% (8/315) based on paraoxonase activity and 2.8% (9/ 324) using arylesterase activity. Higher rates of incident nonfatal MI and stroke were observed in participants within the lowest quartile of paraoxonase activity (37/311 [11.9%]) and arylesterase activity (40/319 [12.5%]). Similarly, lower frequency of 3-year incident all-cause mortality was observed in participants within the highest quartiles of paraoxonase activity (17/ 315 [5.4%]) and arylesterase activity (16/324 [4.9%]) compared with participants within the lowest quartiles of paraoxonase activity (37/311 [11.9%]) and arylesterase activity (41/319 [12.9%]). Lower event rates were also noted for MACE among study participants in the highest quartiles of paraoxonase activity (23/315 [7.3%]) and Table 5. Prevalent Disease by PON1 Q192R Genotype a Prevalent Disease by PON1 Q192R Genotype Coronary artery disease QQ192 QR192 RR192 Peripheral artery disease QQ192 QR192 RR192 Disease Present, No. (%) n = 962 461 (47.9) 410 (42.6) 91 (9.5) n = 363 171 (47.1) 164 (45.2) 28 (7.7) Disease Absent, No. (%) n = 405 169 (41.7) 191 (47.2) 45 (11.1) n = 1023 471 (46.0) 444 (43.3) 108 (10.6) P Value .03 .12 .35 .73 .56 .12 Abbreviation: PON1, paraoxonase 1 gene. a QQ , wild type; QR , heterozygotes; RR , mutant homozygous. Of the 1399 participants, 1386 had genotype 192 192 192 information and 1367 had information on history of coronary artery disease. arylesterase activity (25/324 [7.7%]) compared with event rates for MACE among study participants in the lowest quartiles of paraoxonase activity (78/ 311 [25.1%]) and arylesterase activity (75/319 [23.5%]) (Table 7 and Figure). Following multivariate analysis, serum PON1 activity measures remained independently associated with prospective risk of cardiac events (Table 7). The lowest quartiles of both paraoxonase and arylesterase activity were associated with a greater incident risk of nonfatal MI or stroke (paraoxonase: HR, 4.4; 95% CI, 2.0-9.6; P⬍.001; arylesterase: HR, 4.5; 95% CI, 2.2-9.4; P⬍.001) during the 3-year follow-up interval. The risk of all-cause mortality was also greatest in participants in the lowest quartiles of either paraoxonase activity or arylesterase activity (paraoxonase: HR, 2.4; 95% CI, 1.3-4.4; P =.004; arylesterase: HR, 2.2; 95% CI, 1.2-4.2; P = .01). This translated to a greater incidence of MACE in participants in the lowest quartiles of paraoxonase and arylesterase activity (paraoxonase: HR, 3.4; 95% CI, 2.15.5; P⬍.001; arylesterase: HR, 2.9; 95% CI, 1.8-4.7; P ⬍.001). In further analyses, low systemic levels of arylesterase activity were associated with an increased risk of having a first cardiovascular event (true incidence) among participants without either a history of CVD or angiographic evidence of significant coronary artery disease (defined as ⱖ50% stenosis) at baseline (adjusted HR, 5.8; 95% CI, 1.2-28.6; P=.03) (TABLE 8). In- ©2008 American Medical Association. All rights reserved. creased risk of having a recurrent nonfatal MI or stroke, all-cause mortality, and MACE was observed with low levels of paraoxonase and arylesterase activity measurements in participants with an established diagnosis of CVD at enrollment (TABLE 9). The adjusted HR for recurrent nonfatal MI or stroke was 4.5 (95% CI, 1.9-11.0; P=.001) for paraoxonase and 4.2 (95% CI, 1.9-9.6; P = .001) for arylesterase between the highest and lowest activity quartiles. Participants in the lowest activity quartile for either paraoxonase or arylesterase were more likely to have a recurrent MACE compared with participants in the highest activity quartile (paraoxonase: HR, 3.4; 95% CI, 2.0-5.9; P⬍.001; arylesterase: HR, 2.4; 95% CI, 1.4-3.9; P ⬍.001). Prognostic Value of PON1 The C index for the outcomes of nonfatal MI and stroke during a minimum 3-year follow-up period was 0.59 (with traditional risk factors as predictors), and significantly increased to either 0.66 (P=.007) or 0.69 (P=.001) in comparisons with traditional risk factors when including systemic measures of either paraoxonase or arylesterase activity, respectively. For the composite MACE outcome, addition of either paraoxonase or arylesterase activity to traditional risk factors also significantly increased the predictive value of the model during a mean of 44 (SD, 7) months. For example, the C index of 0.67 with traditional risk factors alone increased to 0.71 (P = .007) and 0.70 (Reprinted) JAMA, March 19, 2008—Vol 299, No. 11 Downloaded from jama.ama-assn.org at University of California - Los Angeles on January 6, 2012 1271 RELATIONSHIP OF PARAOXONASE 1 WITH SYSTEMIC OXIDATIVE STRESS AND CVD RISK Table 6. Relationship Between PON1 Q192R Genotype and CVD Outcomes PON1 Q192R Genotypes CVD Outcomes RR192 ⫹ QR192 (n = 681) Nonfatal MI/CVA No./total (Kaplan-Meier %) 51/681 (7.97) Adjusted HR (95% CI) a 1 [Reference] All-cause mortality No./total (Kaplan-Meier %) 43/681 (6.75) Adjusted HR (95% CI) a 1 [Reference] MACE No./total (Kaplan-Meier %) 88/681 (13.59) Adjusted HR (95% CI) a 1 [Reference] QQ192 (n = 584) 43/584 (7.88) 1.01 (0.65-1.57) 62/584 (11.10) 2.05 (1.32-3.18) 102/584 (18.04) 1.48 (1.09-2.03) P Value .95 .96 .006 .001 .03 .01 Abbreviations: CI, confidence interval; CVA, cerebrovascular accident; CVD, cardiovascular disease; HR, hazard ratio; MACE, major adverse cardiac events (ie, myocardial infarction, stroke, or death); MI, myocardial infarction; PON1, paraoxonase 1 gene. a Adjusted HRs were calculated for each clinical end point by adding the risk allele Q192 to the multivariate model that included the Framingham ATP-III risk score of individual participants (including diabetes and smoking status), race, log C-reactive protein, body mass index, and use of s2tatins and aspirin. (P=.01) for risk factors plus paraoxonase or arylesterase, respectively. COMMENT Oxidative stress is thought to have a pivotal role in the pathogenesis of a number of chronic inflammatory disease processes, including atherosclerosis. The failure of ␣-tocopherol supplementation studies with alleged antioxidant properties to prevent cardiovascular events has brought the oxidation hypothesis of atherosclerosis into question.24 Therefore, the search continues to identify effective strategies that promote systemic antioxidant effects in humans and to determine whether they have a beneficial influence on the rate of CVD. The present genetic and biochemical studies demonstrate that the HDLassociated protein PON1 promotes pronounced systemic antioxidant effects in humans, with coincident links to reduction in coronary artery disease prevalence and prospective risks of MACE. Our results demonstrate that both the PON1 Q192R polymorphism and serum PON1 activity are associated with both prevalent coronary artery disease and incident adverse cardiovascular events. It has been speculated that PON1 contributes to the atheroprotective property of HDL via promotion of a systemic antioxidant effect.1,25 However, there has been no definitive in vivo evi- dence that PON1 promotes systemic antioxidant effects in humans. A recent meta-analysis of 43 studies examining the relationship between PON1 and clinical outcomes15 concluded that considerable uncertainty remained, since published studies to date primarily involve small cohorts of participants. Moreover, they often failed to simultaneously examine PON1 genotype, systemic PON1 activity measures, and cardiovascular outcomes, and none examined systemic quantitative indices of oxidative stress. In this study, we investigated the potential relationship between PON1 activity, PON1 genotype, systemic oxidative stress, and both coronary and peripheral artery disease in a large, prospective cohort of patients to more fully interrogate the aforementioned relationships. We demonstrate that the PON1 Q192R polymorphism is functional, resulting in increased enzymatic activity. In parallel, we demonstrate that elevated systemic levels of multiple structurally distinct fatty acid oxidation products that are increased in both atherosclerotic plaque and plasma of participants with CVD26,27 are associated with low systemic levels of PON1 activity and the PON1 QQ192 genotype. Importantly, the plasma samples analyzed in the present study were collected and processed under conditions designed to prevent artificial oxidation of lipids during both storage and analysis. The find- 1272 JAMA, March 19, 2008—Vol 299, No. 11 (Reprinted) ing that levels of 9-HETE, an isomer of arachidonic acid oxidation produced exclusively by free radical–mediated processes, are higher in participants with low levels of paraoxonase activity suggests that PON1 can influence oxidative events beyond the cyclooxygenase and lipoxygenase pathways. While the mechanism(s) for PON1mediated systemic antioxidant effects remains to be determined, the present findings strongly support a role for this HDL-associated protein in modulating systemic oxidative stress in humans. Moreover, the present study suggests an important mechanistic link among PON1, systemic oxidative stress, and risk of development of atherosclerotic heart disease and its acute complications. In a recent publication, amino acid position 192 of PON1 was suggested to participate in HDL binding.28 The Q192 alloenzyme was shown to bind to the HDL particle with 3-fold lower affinity than the R192 alloenzyme and, consequently, exhibited lower stability, lipolactonase activity, and modulatory effect on macrophage cholesterol efflux. The findings of the present clinical study complement these results, demonstrating that individuals with the arginine (R) mutation at position 192 have higher serum levels of PON1 activity, lower systemic indices of systemic oxidative stress, and corresponding reductions in both prevalent coronary artery disease and prospective cardiac events. This is, to our knowledge, the first large, prospective study that comprehensively examines the genetics and biochemical activity of PON1 using dual enzyme measurements to predict prevalent disease risks, as well as prospective risk of MACE, while simultaneously also examining whether a potential mechanism for CVD associations is linked to systemic oxidative stress measures in humans. Compared with prior studies that often focused on lower-risk populations,15 the present cohort is also substantially enriched with patients with both coronary and peripheral artery disease, enabling us to better study the association of the PON1 Q192R variant on ©2008 American Medical Association. All rights reserved. Downloaded from jama.ama-assn.org at University of California - Los Angeles on January 6, 2012 RELATIONSHIP OF PARAOXONASE 1 WITH SYSTEMIC OXIDATIVE STRESS AND CVD RISK Figure. Relationship Between Functional Genetic and Biochemical Indices of PON1 Activity and Incident Cardiovascular Risk MI/CVA Major Adverse Cardiac Events Q192R Polymorphism Q192R Polymorphism 30 Genotype QQ RQ RR 10 25 Events, % Events, % 15 5 20 15 10 5 Log-rank P = .95 Log-rank P = .02 0 0 0 183 365 549 732 915 1095 0 183 365 Days No. at Risk QQ RQ RR 584 563 118 544 522 112 527 517 110 501 492 101 732 915 1095 487 477 101 446 440 91 435 429 90 Days 487 477 101 446 440 91 435 429 90 584 563 118 544 522 112 Arylesterase Activity 527 517 110 501 492 101 Arylesterase Activity 30 15 Quartile 1 2 3 4 10 25 Events, % Events, % 549 5 20 15 10 5 Log-rank P <.001 Log-rank P <.001 0 0 0 183 365 549 732 915 1095 0 183 365 Days No. at Risk Quartile 1 Quartile 2 Quartile 3 Quartile 4 319 317 318 324 289 290 300 310 278 284 296 307 263 267 277 297 549 732 915 1095 253 259 270 293 234 233 249 270 223 228 245 267 Days 253 259 270 293 234 233 249 270 223 228 245 267 319 317 318 324 289 290 300 310 Paraoxonase Activity 278 284 296 307 263 267 277 297 Paraoxonase Activity 15 30 10 Events, % Events, % 25 5 20 15 10 5 Log-rank P <.001 0 0 183 365 549 732 915 1095 Log-rank P <.001 0 0 183 365 Days No. at Risk Quartile 1 Quartile 2 Quartile 3 Quartile 4 311 326 326 315 280 303 305 301 267 299 301 298 254 288 283 279 549 732 915 1095 288 290 275 277 283 259 253 257 279 256 247 252 Days 288 290 275 277 283 259 253 257 279 256 247 252 311 326 326 315 280 303 305 301 267 299 301 298 254 288 283 279 The paraoxonase 1 (PON1) Q192R genotypes are as follows: RR192 (mutant homozygous); QR192 (heterozygous); and QQ192 (wild type). In the top panels, log-rank P values are shown for the at-risk genotype QQ192 vs RR192⫹QR192. PON1 activity (paraoxonase and arylesterase) were categorized into quartiles; Q1: lowest activity quartile; Q4: highest activity quartile. For paraoxonase, Q4⬎1640, Q3=1640-899.1, Q2=899-450, and Q1⬍450 nmol/min/mL. For arylesterase, Q4⬎403.9, Q3=403.9-338.5, Q2=338.4-283, and Q1⬍283 µmol/min/mL. In the middle and bottom panels, log-rank P values across activity quartiles are shown. Days indicates number of days from enrollment to first cardiac event. Event rates were calculated at 6-month intervals. Y-axis scales in blue indicate range from 0% to 15%. MI indicates myocardial infarction; CVA, cerebrovascular accident. ©2008 American Medical Association. All rights reserved. (Reprinted) JAMA, March 19, 2008—Vol 299, No. 11 Downloaded from jama.ama-assn.org at University of California - Los Angeles on January 6, 2012 1273 RELATIONSHIP OF PARAOXONASE 1 WITH SYSTEMIC OXIDATIVE STRESS AND CVD RISK Table 7. Risk of Prospective CVD Events in Relation to PON1 Activity Among Entire Cohort Paraoxonase Activity Quartile, nmol/min/mL CVD Events Nonfatal MI/CVA No./total (%) Unadjusted HR (95% CI) Adjusted HR (95% CI) a All-cause mortality No./total (%) Unadjusted HR (95% CI) Adjusted HR (95% CI) a MACE No./total (%) Unadjusted HR (95% CI) Adjusted HR (95% CI) a Quartile 4 (⬎ 1640) (n = 315) Quartile 3 (1640-899.1) (n = 325) Quartile 2 (899-450) (n = 326) Quartile 1 (⬍ 450) (n = 311) 8/315 (2.5) 1 [Reference] 1 [Reference] 23/325 (7.4) 27/326 (8.3) 37/311 (11.9) 3.0 (1.3-6.7) 3.3 (1.5-7.3) 5.0 (2.3-10.8) 2.9 (1.3-6.4) 3.1 (1.4-7.0) 4.4 (2.0-9.6) 17/315 (5.4) 1 [Reference] 1 [Reference] 21/325 (6.5) 27/326 (8.3) 37/311 (11.9) 1.2 (0.6-2.3) 1.3 (0.7-2.5) 2.7 (1.5-4.7) 1.1 (0.6-2.2) 1.3 (0.7-2.6) 2.4 (1.3-4.4) 23/315 (7.3) 1 [Reference] 1 [Reference] 44/325 (13.5) 48/326 (14.7) 78/311 (25.1) 1.9 (1.2-3.2) 2.1 (1.2-3.4) 3.7 (2.3-5.9) 1.9 (1.1-3.2) 2.0 (1.2-3.4) 3.4 (2.1-5.5) Arylesterase Activity Quartile, µmol/min/mL Quartile 4 (⬎ 403.9) (n = 324) Nonfatal MI/CVA No./total (%) Unadjusted HR (95% CI) Adjusted HR (95% CI) a All-cause mortality No./total (%) Unadjusted HR (95% CI) Adjusted HR (95% CI) a MACE No./total (%) Unadjusted HR (95% CI) Adjusted HR (95% CI) a 9/324 (2.8) 1 [Reference] 1 [Reference] Quartile 3 (403.9-338.5) (n = 318) Quartile 2 (338.4-283) (n = 316) Quartile 1 (⬍ 283) (n = 319) 15/318 (4.7) 32/316 (10.1) 40/319 (12.5) 1.7 (0.8-4.0) 3.8 (1.8-8.0) 4.8 (2.3-9.9) 1.5 (0.7-3.6) 3.4 (1.6-7.2) 4.5 (2.2-9.4) 16/324 (4.9) 1 [Reference] 1 [Reference] 22/318 (6.9) 27/316 (8.5) 41/319 (12.9) 1.4 (0.8-2.7) 1.8 (1.0-3.3) 2.7 (1.5-4.8) 1.4 (0.7-2.8) 1.8 (0.9-3.4) 2.2 (1.2-4.2) 25/324 (7.7) 1 [Reference] 1 [Reference] 37/318 (11.6) 56/316 (17.7) 75/319 (23.5) 1.6 (0.9-2.6) 2.4 (1.5-3.9) 3.3 (2.1-5.2) 1.5 (0.9-2.5) 2.3 (1.4-3.7) 2.9 (1.8-4.7) Abbreviations: CI, confidence interval; CVA, cerebrovascular accident; CVD, cardiovascular disease; HR, hazard ratio; MACE, major adverse cardiac events (ie, myocardial infarction, stroke, or death); MI, myocardial infarction; PON1, paraoxonase 1. a Adjusted HRs were calculated including Framingham ATP-III risk score (including diabetes status), log C-reactive protein, body mass index, and use of statins and aspirin. Table 8. Incidence of First CVD Event in Relation to PON1 Activity Among Patients With No CVD at Baseline Paraoxonase Activity Quartile, nmol/min/mL First CVD Events MACE No./total (%) Unadjusted HR (95% CI) Adjusted HR (95% CI) a Quartile 4 (⬎ 1640) (n = 80) Quartile 3 (1640-899.1) (n = 61) Quartile 2 (899-450) (n = 57) Quartile 1 (⬍ 450) (n = 56) 4/80 (5.0) 1 [Reference] 1 [Reference] 3/61 (4.9) 1.0 (0.2-4.4) 0.6 (0.1-3.1) 3/57 (5.3) 1.1 (0.2-4.7) 1.2 (0.3-5.3) 8/56 (14.3) 3.0 (0.9-9.8) 1.7 (0.5-6.4) Arylesterase Activity Quartile, µmol/min/mL MACE No./total (%) Unadjusted HR (95% CI) Adjusted HR (95% CI) a Quartile 4 (⬎ 403.9) (n = 91) Quartile 3 (403.9-338.5) (n = 61) 2/91 (2.2) 1 [Reference] 1 [Reference] 3/61 (4.9) 2.2 (0.4-13.4) 2.1 (0.3-12.8) Quartile 2 (338.4-283) (n = 51) Quartile 1 (⬍ 283) (n = 51) 3/51 (5.9) 10/51 (19.6) 2.7 (0.4-16.0) 9.7 (2.1-44.1) 2.4 (0.4-14.4) 5.8 (1.2-28.6) Abbreviations: CI, confidence interval; CVD, cardiovascular disease; HR, hazard ratio; MACE, major adverse cardiac events (ie, myocardial infarction, stroke, or death); PON1, paraoxonase 1. a Adjusted HRs were calculated including the Framingham ATP-III risk score (including diabetes status), log C-reactive protein, body mass index, and use of statins and aspirin. 1274 JAMA, March 19, 2008—Vol 299, No. 11 (Reprinted) coronary and extracoronary atherosclerosis after controlling for established risk factors. However, in the present study, diminished systemic arylesterase activity levels remained significantly associated with 3-year incident MACE (composite of MI, stroke, or death) in the subgroup of participants without history or clinical evidence of either coronary or peripheral artery disease and without significant angiographic evidence of coronary artery disease. Furthermore, the present study also suggests a potential prognostic value of PON1 activity measurement in patients. The addition of systemic PON1 activity measures to traditional risk factors and C-reactive protein provided a significant incremental improvement in the ability to predict clinical outcome during a 3-year period. Inrecentgenomewideassociationstudies,29,30 evidence for association of the PON1 gene with coronary artery disease or MI was not observed. One likely explanation is that genomewide association studies are not ideal for assessing specific candidate genes since the panel of SNPs placed on the chip may not capture all of the genetic variation for any particular gene.31 Furthermore, genomewide association studies tend to identify genes with the strongest genetic effect because of the stringency in considering what is statistically significant. In this regard, a weak associationwithPON1mayhavebeenobserved but not reported since it did not exceedthethresholdforsignificancebased on the thousands of statistical tests performed. Moreover, one of the phenotypes wehavestudied,namely,mortalitywithin thecohort,wasnottheexactsameasthose studied in the genomewide association studies,whichcouldalsoaccountforwhy an association with PON1 was not reported. ThePON1geneislocatedincloseproximity to 2 other members of the paraoxonasegenefamily,PON2andPON3,which raises the possibility that the association we have observed could be due to linkage disequilibrium with a variant in either PON2orPON3.However,anexamination of the HAPMAP database for the white population(whichmatchestheGeneBank ©2008 American Medical Association. All rights reserved. Downloaded from jama.ama-assn.org at University of California - Los Angeles on January 6, 2012 RELATIONSHIP OF PARAOXONASE 1 WITH SYSTEMIC OXIDATIVE STRESS AND CVD RISK cohort)showsthatthePON1Q192Rvariant is located within a haplotype block in which the SNPs are in moderate linkage disequilibriumwitheachother.Moreover, this haplotype block only covers PON1 and does not extend to the PON2 and PON3genes.Thus,whileitispossiblethat variants of PON2 and PON3 could also contribute to altered PON1 activity and CVDrisk,theavailableHAPMAPdatasuggest that, in whites, the association observed with the PON1 Q192R variant is not due to linkage disequilibrium with other variants in the adjacent PON2 and PON3genes,butlikelyresultsfromaltered PON1function.Althoughitispossiblethat other PON1 variants that are in linkage disequilibrium with the Q192R SNP are the causal alleles, recent studies using recombinantPON1mutantsatposition192 have shown that this residue is important for HDL binding28 and suggest that the Q192R SNP is the causal variant underlying the clinical associations that we observe. Further studies to address this issue are required. A number of caveats to the present study should be noted. All participants presented for elective diagnostic coronary angiography, which limits the generalization of the present findings and raises the possibility of selection bias. However, measurements of both paraoxonase and arylesterase activity in healthy volunteers are higher than in both the CVD and non-CVD participants evaluated in the present study. Given that the majority of participants were white, it also remains to be determined whether the same relationships are observed in large cohorts of other racial/ethnic groups. The lower percentage of nonwhite participants in our cohort (7.9%) limited our ability to analyze the mechanistic link between PON1 and global indices of systemic oxidative stress and thereby its association with CVD in participants of other racial/ethnic backgrounds. It is also uncertain whether concomitant medical problems influenced PON1 activity; however, inclusion of traditional cardiac risk factors including age, sex, race, diabetes, hypertension, smoking, lipids, C- reactive protein, and medication use (statins and aspirin) in the multivariate analyses failed to alter the results. It also is possible that the reduction in PON1 activity may result from the presence of vascular disease rather than be the direct cause of future CVD events. However, arylesterase activity levels in participants with minimal angiographic evidence of coronary artery disease still remained independently associated with incident cardiovascular events over the ensuing 3-year interval following participant enrollment, consistent with an association of PON1 in early macrovascular atherosclerotic disease processes. Regardless, the current observation of a relationship among the PON1 Q192R polymorphism and PON1 activity, oxidative stress, and CVD outcomes is consistent with findings in murine models and provides evidence that the PON1 protein protects against the development and propagation of CVD. CONCLUSION The current findings provide direct prospective evidence of an important mechanistic link between the PON1 gene and PON1 systemic activity measures with both multiple quantitative indices of oxidative stress and atherosclerotic heart disease development in humans. Paraoxonase 1 is almost exclusively found to be associated with HDL particles within the circulation and has been argued to promote some of the Table 9. Risk of Recurrent CVD Events in Relation to PON1 Activity Among Patients With CVD at Baseline Paraoxonase Activity Quartile, nmol/min/mL Recurrent Events Nonfatal MI/CVA No./total (%) Unadjusted HR (95% CI) Adjusted HR (95% CI) a All-cause mortality No./total (%) Unadjusted HR (95% CI) Adjusted HR (95% CI) a MACE No./total (%) Unadjusted HR (95% CI) Adjusted HR (95% CI) a Quartile 4 (⬎ 1640) (n = 235) Quartile 3 (1640-899.1) (n = 265) Quartile 2 (899-450) (n = 269) Quartile 1 (⬍ 450) (n = 255) 6/235 (2.6) 1 [Reference] 1 [Reference] 22/265 (8.3) 24/269 (8.9) 40/255 (12.5) 3.4 (1.4-8.3) 3.6 (1.5-8.7) 5.3 (2.2-12.8) 3.2 (1.3-7.9) 3.3 (1.3-8.1) 4.5 (1.9-11.0) 15/235 (6.4) 1 [Reference] 1 [Reference] 20/265 (7.5) 24/269 (8.9) 40/255 (15.7) 1.2 (0.6-2.3) 1.4 (0.7-2.6) 2.6 (1.4-4.6) 1.1 (0.6-2.3) 1.4 (0.7-2.7) 2.4 (1.3-4.5) 19/235 (8.1) 1 [Reference] 1 [Reference] 41/265 (15.5) 45/269 (16.7) 70/255 (27.5) 2.0 (1.2-3.4) 2.1 (1.2-3.6) 3.7 (2.2-6.2) 2.0 (1.1-3.5) 2.1 (1.2-3.7) 3.4 (2.0-5.9) Arylesterase Activity Quartile, µmol/min/mL Quartile 4 (⬎ 403.9) (n = 233) Nonfatal MI/CVA No./total (%) Unadjusted HR (95% CI) Adjusted HR (95% CI) a All-cause mortality No./total (%) Unadjusted HR (95% CI) Adjusted HR (95% CI) a MACE No./total (%) Unadjusted HR (95% CI) Adjusted HR (95% CI) a Quartile 3 (403.9-338.5) (n = 257) Quartile 2 (338.4-283) (n = 266) Quartile 1 (⬍ 283) (n = 268) 7/233 (3.0) 1 [Reference] 1 [Reference] 12/257 (4.7) 30/266 (11.3) 35/268 (13.1) 1.6 (0.6-4.1) 4.0 (1.8-9.1) 4.6 (2.0-10.4) 1.4 (0.5-3.6) 3.4 (1.5-7.9) 4.2 (1.9-9.6) 16/233 (6.9) 1 [Reference] 1 [Reference] 22/257 (8.6) 26/266 (9.8) 35/268 (13.1) 1.3 (0.7-2.4) 1.5 (0.8-2.8) 1.9 (1.1-3.5) 1.3 (0.6-2.5) 1.4 (0.7-2.8) 1.7 (0.9-3.2) 23/233 (9.9) 1 [Reference] 1 [Reference] 34/257 (13.2) 53/266 (19.9) 65/268 (24.3) 1.4 (0.8-2.4) 2.2 (1.3-3.5) 2.6 (1.6-4.2) 1.3 (0.7-2.3) 2.0 (1.2-3.3) 2.4 (1.4-3.9) Abbreviations: CI, confidence interval; CVA, cerebrovascular accident; CVD, cardiovascular disease; HR, hazard ratio; MACE, major adverse cardiac events (ie, myocardial infarction, stroke, or death); MI, myocardial infarction; PON1, paraoxonase 1. a Adjusted HRs were calculated including the Framingham ATP-III risk score (including diabetes status), log C-reactive protein, body mass index, and use of statins and aspirin. ©2008 American Medical Association. All rights reserved. (Reprinted) JAMA, March 19, 2008—Vol 299, No. 11 Downloaded from jama.ama-assn.org at University of California - Los Angeles on January 6, 2012 1275 RELATIONSHIP OF PARAOXONASE 1 WITH SYSTEMIC OXIDATIVE STRESS AND CVD RISK anti-inflammatory and antioxidant effects attributed to HDL. Thus, the present studies also provide further support for the concept that functional properties beyond the ability of HDL and its associated proteins to promote reverse cholesterol transport contribute to the overall ability of this lipoprotein to reduce or prevent development of atherosclerosis. Author Affiliations: Department of Cell Biology (Drs Bhattacharyya, Nicholls, Zhang, M.-L. Brennan, and Hazen and Messrs Schmitt and Fu), Center for Cardiovascular Diagnostics and Prevention (Drs Bhattacharyya, Nicholls, Zhang, M.-L. Brennan, and Hazen and Messrs Schmitt and Fu), and Department of Cardiovascular Medicine (Drs Nicholls, Ellis, and Hazen, Mr Shao, and Ms D. M. Brennan), Cleveland Clinic, Cleveland, Ohio; Scripps Research Institute and Scripps Clinic, La Jolla, California (Dr Topol); Departments of Medicine (Drs Yang and Lusis), Human Genetics (Dr Lusis), and Microbiology and Immunology (Dr Lusis), University of California, and Department of Preventive Medicine, University of Southern California (Dr Allayee), Los Angeles. Author Contributions: Dr Hazen had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. Study concept and design: Lusis, Hazen. Acquisition of data: Bhattacharyya, Topol, Zhang, Yang, Schmitt, Fu, M.-L. Brennan, Allayee, Lusis, Hazen. Analysis and interpretation of data: Bhattacharyya, Nicholls, Topol, Yang, Shao, D. M. Brennan, Ellis, Lusis, Hazen. Drafting of the manuscript: Bhattacharyya, Nicholls, Lusis, Hazen. Critical revision of the manuscript for important intellectual content: Bhattacharyya, Zhang, Schmitt, Fu, Shao, D. M. Brennan, M.-L. Brennan, Allayee, Lusis, Hazen. Statistical analysis: Bhattacharyya, Nicholls, Yang, Shao, D. M. Brennan, Lusis, Hazen. Obtained funding: Topol, Allayee, Lusis, Hazen. Administrative, technical, or material support: Bhattacharyya, Topol, Zhang, M.-L. Brennan, Allayee, Lusis, Hazen. Study supervision: Nicholls, Topol, Ellis, M.-L. Brennan, Hazen. Financial Disclosures: Dr Nicholls reports that he has received speaking honoraria from Pfizer, AstraZeneca, Merck, Schering-Plough, and Takeda and consulting fees from Roche, AstraZeneca, Pfizer, and NovoNordisk. Dr Ellis reports that he is named as coinventor on pending patents filed by the Cleveland Clinic referring to use of genetic markers in cardiovascular disease. Dr Ellis also reports that he has received research support from Lilly, Centocor, Boston Scientific, and Cardlodx and has received consulting fees from Celera, Cordis, Boston Scientific, and Abbott Vascular. Dr Allayee reports that he is named as coinventor on pending and approved patents filed by the Univer- sity of California that are related to identification of 5-lipoxygenase as a gene involved in cardiovascularand diabetes-related traits. Dr Allayee also reports that he has received honoraria/consulting fees from Reliant Pharmaceuticals and Corautus Genetics. Dr Lusis reports that he has received research grant support from Merck, Anthera, and Bristol-Myers Squibb. Dr Hazen reports that he is named as coinventor on pending and approved patents filed by the Cleveland Clinic that refer to the use of biomarkers to inflammatory and cardiovascular diseases. Dr Hazen also reports that he is the scientific founder of PrognostiX Inc; has received research grant support from Abbott Diagnostics, Pfizer, Merck, PrognostiX Inc, Hawaii Biotech, ArgiNOx, Sanofi, and Takeda; and has received honoraria and consulting fees from Abbott Diagnostics, BioSite, Merck, Lilly, Pfizer, PrognostiX Inc, Wyeth, BioPhysical, and AstraZeneca. No other disclosures were reported. Funding/Support: This study was supported by National Institutes of Health grants P01 HL076491, P01 HL077107, P01 HL087018, P01 HL30568, and RO1 HL079353, the Case Western Reserve University/ Cleveland Clinic Clinical and Translational Science Award (grant 1KL2RR024990), and AHA grant 0435223N. A portion of this work was conducted in a facility constructed with support from Research Facilities Improvement Program grant C06 (RR1060001, CA62528-01, RR14514-01) from the National Center for Research Resources. Role of the Sponsors: The funding organizations played no role in the design and conduct of the study; in the collection, analysis, and interpretation of the data; or in the preparation, review, or approval of the manuscript. Human paraoxonase 1 gene polymorphisms and the risk of coronary heart disease: a community-based study. Cardiology. 2002;98(3):116-122. 12. Robertson KS, Hawe E, Miller GJ, et al. Human paraoxonase gene cluster polymorphisms as predictors of coronary heart disease risk in the prospective Northwick Park Heart Study II. Biochim Biophys Acta. 2003;1639(3):203-212. 13. Lawlor DA, Day IN, Gaunt TR, et al. The association of the PON1 Q192R polymorphism with coronary heart disease: findings from the British Women’s Heart and Health cohort study and a meta-analysis. BMC Genet. 2004;5:17. 14. Brophy VH, Hastings MD, Clendenning JB, et al. Polymorphisms in the human paraoxonase (PON1) promoter. Pharmacogenetics. 2001;11(1):77-84. 15. Wheeler JG, Keavney BD, Watkins H, et al. Four paraoxonase gene polymorphisms in 11 212 cases of coronary heart disease and 12 786 controls: metaanalysis of 43 studies. Lancet. 2004;363(9410):689695. 16. Mackness B, Davies GK, Turkie W, et al. Paraoxonase status in coronary heart disease: are activity and concentration more important than genotype? Arterioscler Thromb Vasc Biol. 2001;21(9):14511457. 17. Zhang R, Brennan ML, Shen Z, et al. Myeloperoxidase functions as a major enzymatic catalyst for initiation of lipid peroxidation at sites of inflammation. J Biol Chem. 2002;277(48):46116-46122. 18. Newman TB, Brown AN. Use of commercial record linkage software and vital statistics to identify patient deaths. J Am Med Inform Assoc. 1997;4(3): 233-237. 19. Eckerson HW, Wyte CM, La Du BN. The human serum paraoxonase/arylesterase polymorphism. Am J Hum Genet. 1983;35(6):1126-1138. 20. Pencina MJ, D’Agostino RB. Overall C as a measure of discrimination in survival analysis: model spe- cific population value and confidence interval estimation. Stat Med. 2004;23(13):2109-2123. 21. Kattan MW. Validating a prognostic model. Cancer. 2006;107(11):2523-2524. 22. Harrell FE, Califf RM, Pryor DB, et al. Evaluating the yield of medical tests JAMA. 1982;247(18):25432546. 23. Kattan MW. Evaluating a new marker’s predictive contribution. Clin Cancer Res. 2004;10(3):822824. 24. Vivekananthan DP, Penn MS, Sapp SK, et al. Use of antioxidant vitamins for the prevention of cardiovascular disease: meta-analysis of randomised trials. Lancet. 2003;361(9374):2017-2023. 25. Lusis AJ, Fogelman AM, Fonarow GC. Genetic basis of atherosclerosis, I: new genes and pathways. Circulation. 2004;110(13):1868-1873. 26. Thukkani AK, McHowat J, Hsu FF, et al. Identification of alpha-chloro fatty aldehydes and unsaturated lysophosphatidylcholine molecular species in human atherosclerotic lesions. Circulation. 2003;108 (25):3128-3133. 27. Shishehbor MH, Zhang R, Medina H, et al. Systemic elevations of free radical oxidation products of arachidonic acid are associated with angiographic evidence of coronary artery disease. Free Radic Biol Med. 2006;41(11):1678-1683. 28. Gaidukov L, Rosenblat M, Aviram M, et al. The 192R/Q polymorphs of serum paraoxonase PON1 differ in HDL binding, lipolactonase stimulation, and cholesterol efflux. J Lipid Res. 2006;47(11):2492-2502. 29. Samani NJ, Erdmann J, Hall AS, et al. Genomewide association analysis of coronary artery disease. N Engl J Med. 2007;357(5):443-453. 30. Rosenzweig A. Scanning the genome for coronary risk. N Engl J Med. 2007;357(5):497-499. 31. Shriner D, Vaughn LK, Padilla MA, Tiwari HK. Problems with genome-wide association studies. Science. 2007;316(5833):1840-1842. REFERENCES 1. Mackness M, Mackness B. Paraoxonase 1 and atherosclerosis: is the gene or the protein more important? Free Radic Biol Med. 2004;37(9):1317-1323. 2. Li HL, Liu DP, Liang CC. Paraoxonase gene polymorphisms, oxidative stress, and diseases. J Mol Med. 2003;81(12):766-779. 3. Durrington PN, Mackness B, Mackness MI. Paraoxonase and atherosclerosis. Arterioscler Thromb Vasc Biol. 2001;21(4):473-480. 4. Mackness MI, Arrol S, Durrington PN. Paraoxonase prevents accumulation of lipoperoxides in lowdensity lipoprotein. FEBS Lett. 1991;286(1-2):152154. 5. Shih DM, Xia YR, Wang XP, et al. Combined serum paraoxonase knockout/apolipoprotein E knockout mice exhibit increased lipoprotein oxidation and atherosclerosis. J Biol Chem. 2000;275(23):1752717535. 6. Shih DM, Gu L, Xia YR, et al. Mice lacking serum paraoxonase are susceptible to organophosphate toxicity and atherosclerosis. Nature. 1998;394(6690): 284-287. 7. Tward A, Xia YR, Wang XP, et al. Decreased atherosclerotic lesion formation in human serum paraoxonase transgenic mice. Circulation. 2002;106(4): 484-490. 8. Rozenberg O, Rosenblat M, Coleman R, et al. Paraoxonase (PON1) deficiency is associated with increased macrophage oxidative stress: studies in PON1knockout mice. Free Radic Biol Med. 2003;34(6): 774-784. 9. Reddy ST, Waldeigh DJ, Grijalva V, et al. Human paraoxonase-3 is an HDL-associated enzyme with biological activity similar to paraoxonase-1 protein but is not regulated by oxidized lipids. Arterioscler Thromb Vasc Biol. 2001;21(4):542-547. 10. Hegele RA. Paraoxonase genes and disease. Ann Med. 1999;31(3):217-224. 11. Watzinger N, Schmidt H, Schumacher M, et al. 1276 JAMA, March 19, 2008—Vol 299, No. 11 (Reprinted) ©2008 American Medical Association. All rights reserved. Downloaded from jama.ama-assn.org at University of California - Los Angeles on January 6, 2012
© Copyright 2026 Paperzz