- Open Access
Interaction between polyphenols intake and PON1 gene variants on markers of cardiovascular disease: a nutrigenetic observational study
Journal of Translational Medicinevolume 14, Article number: 186 (2016)
Paraoxonase 1 (PON1) gene polymorphisms and polyphenols intake have been reported independently associated to lipid profile and susceptibility to atherosclerosis and cardiovascular disease. However, the interaction between these factors remains to be investigated. We performed an observational nutrigenetic study to examine whether the interaction between polyphenols and anthocyanins intake and PON1 genetic variants can modulate biomarkers of cardiovascular health in an Italian healthy population.
We recruited 443 healthy volunteers who participated in the EC funded ATHENA project (AnThocyanin and polyphenols bioactive for Health Enhancement through Nutritional Advancement). Data collection included detailed demographic, clinical, dietary, lifestyle, biochemical and genetic data. Polyphenols and anthocyanins intake was measured by 24 h dietary recall repeated three times a year in order to get seasonal variations. We tested the interaction between 18 independent tagging SNPs in PON1 gene and polyphenols intake on HDL, LDL, cholesterol, triglycerides and atherogenic index of plasma.
Without considering the genetic background, we could not observe significant differences in the lipid profile between high and low polyphenols and anthocyanins intake. Using a nutrigenetic approach, we identified protective genotypes in four independent polymorphisms that, at Bonferroni level (p ≤ 0.0028), present a significant association with increased HDL level under high polyphenols and anthocyanins intake, compared to risk genotypes (rs854549, Beta = 4.7 per C allele; rs854552, Beta = 5.6 per C allele; rs854571, Beta = 3.92 per T allele; rs854572, Beta = 3.94 per C allele).
We highlight the protective role of genetic variants in PON1 towards cardiovascular risk under high polyphenols and anthocyanins consumption. PON1 variants could represent novel biomarkers to stratify individuals who might benefit from targeted dietary recommendation for health promotion and strategies of preventive medicine.
Cardiovascular disease (CVD) is responsible of approximately one-third of deaths in the world and its risk factors continue to increase [1, 2]. Since 2007, the ‘Grand challenges in chronic non-communicable diseases’ outlined that promoting salutary lifestyle and increasing the availability and consumption of healthy food have a key role in the protection against chronic diseases, a global epidemic and economic burden for our society [3, 4].
The physiological response to diet can be modulated by the genetic background, with common DNA variants having the potential to affect specific functional pathways and making subjects more or less responder to a specific diet . For instance, changes in plasma lipid concentration in response to dietary fat intake are largely genetically controlled [6, 7].
Epidemiological studies support the preventive effect of anthocyanins and polyphenols towards the onset of CVD [8–11] through their great antioxidant and antiatherosclerotic activity [12, 13]. Animal and in vitro studies are in favor of their potential to influence lipid profile, a commonly used biomarker of cardiovascular risk . They can slow or inhibit the absorption of lipids and glucose in the intestine and inhibit cholesterol synthesis , which results in a decrease in serum triglyceride, total cholesterol and non-HDL cholesterol and in increase in serum HDL. Moreover, they promote reverse cholesterol transport, which also contributes to their effect on lipid profile [12, 15, 16].
Polyphenols are the most abundant antioxidants in the diet. Their dietary intake can be as high as 1 g/dl, much higher than that of all other known dietary antioxidants. Their main dietary sources are fruits and plant-derived beverages (fruit juices, tea, coffee and red wine), vegetables, cereals, chocolate and dry legumes.
Anthocyanins are plant pigments members of the flavonoid family of polyphenols, providing the bright red–orange to blue–violet colors present in many fruits and vegetables.
Interventional studies have consistently reported that polyphenols and athocyanins-rich diets can modulate paraoxonase 1 (PON1) activity and/or the level of PON1 expression [17–19]. PON1 is a member of a gene family, which also comprises PON2 and PON3, all clustered in tandem on the long arm of human chromosome 7 (q21.22). PON1 enzyme is a glycoprotein with hydrolytic activity , which associates with HDL in the circulation. It metabolizes a broad variety of substrates, and is a primary determinant of the antioxidant and anti-inflammatory activity of HDL, promoting the HDL-mediated macrophage cholesterol efflux [15, 16]. As such, PON1 plays a relevant role in determining susceptibility to atherosclerosis and cardiovascular disease [17, 21, 22]. It is recognized that PON1 activity is genetically regulated, with Single Nucleotide Polymorphisms (SNPs) exhibiting strong association with its arylesterase and paraoxanase activities [17, 20, 23]. By the way, the wide range of serum PON1 activity among individuals is only partially explained by genetic polymorphisms. The presence of additional factors known to modulate PON1 activity and HDL (e.g. dietary factors, lifestyle, statins, etc.) also needs to be taken into account.
The EC funded the Athena project (FP7-KBBE-2009-3) “Anthocyanin and polyphenols bioactive for health enhancement through nutritional advancement”. Within this project, we conducted an observational study with the aim to explore the relationship between SNPs in PON1 and lipid profile as biomarker of cardiovascular health, taking into consideration polyphenols and anthocyanins consumption. We selected HDL, LDL, total cholesterol, triglycerides and the atherogenic index of plasma (AIP) as common biomarkers, broadly used in clinical practice and with limited inter- and intra assay variability.
Study design and participants
We performed an observational nutrigenetic study on 500 volunteers (age range 20–85 years), as defined in the ATHENA project. Participants were recruited in two healthcare centers in Milano (Italy), ASST Santi Paolo e Carlo, Division of Nephrology, and IRCCS Galeazzi Orthopedic Institute, from June 2012 to December 2013. The subjects participating to the study were required to be in good health, defined according to individual medical history and to routine biochemical, urine and instrumental examination (electrocardiogram, ambulatory blood pressure monitoring, echocardiography, carotid artery and renal ultrasound). Exclusion criteria included evidence of cardiovascular or coronary artery disease, neoplasm, psychosis and diabetes mellitus. The ethic committees of both centres approved the study, (Register 441 signed off on 9 June 2011 and Register 45-1250-ATHENA signed off on 22 February 2013, respectively), in accordance with principles of the Declaration of Helsinki and all participants gave written informed consent.
Volunteers underwent a rigorous clinical examination. At first visit, we collected baseline anthropometric and metabolic measurements and information on diet, personal medical history and lifestyle. A blood sample was also taken for the genetic analysis. We repeated the nutritional interview four times throughout the year in order to get seasonal diet variations. Subjects with less than 3 dietary recall interviews were excluded from the analysis. We then calculated the average annual dietary anthocyanins and polyphenols intake for each participant, assuming that this reflects the individual dietary habit. We divided participants in two groups according to their low or high antioxidant intake, considering anthocyanin and polyphenols separately.
Nutritional data collection
Experienced dieticians carried out the dietary assessment through 24 h recall interviews. Food portions were defined using a Food Atlas (edited by Scotti-Bassani), including 99 tables with photographs. 24 h recall data were recorded at least three times in a year, once during the first clinical visit and then every 3–4 month.
For data collection and dietary nutrient composition estimate we used the Diet Monitoring Solution (DMS) software, designed and developed by KOS genetic s.r.l., in the context of the Athena project . DMS calculates macro- and micro-nutrients composition using composition tables from the Food Composition Database for Epidemiological Studies in Italy (BDA-IEO), integrated with tables from the USDA Database for the Flavonoid Content and from the Phenol Explorer 2.0 Database.
Serum lipids analysis
Determination of HDL, LDL, total cholesterol and triglycerides was carried out on automatic clinical chemistry instrumentation (Abbott Architect 8000). Samples were analyzed in batches to limit the analytical variations, after storage at −40 °C for a time not exceeding one month. AIP was calculated as the logarithmically transformed ratio of triglycerides and HDL concentration in plasma .
Genotyping and imputing
Genomic DNA was extracted from peripheral blood using the commercial NucleoSpin Blood Kit by Macherey–Nagel (Macherey–Nagel, Düren, Germany) following the manufacturer’s protocol, and stored at −20°. Forty hundred and seventy-one samples were genotyped using the Illumina HumanCore array (Illumina Inc, San Diego, CA, USA) that contains about 300 K highly informative genome-wide tag-SNPs including indels and updated exome-focused markers. Genotypes were assigned all simultaneously using Genome Studio software (version V2011.1 genotyping module 1.9.4). Genome-wide imputation was performed using MINIMAC software [26, 27] and HapMap CEU haplotypes (release 22) as reference. Measured SNPs with call rate ≥99 % and minor allele frequency (MAF) ≥1 % were included in the dataset. Imputed SNPs with low imputation quality (r2 <0.7) were not used in the association analysis.
PON1 gene maps to chromosome 7 (chr7:94764924-94791780). To select SNPs for the analyses, we first reviewed all SNPs in this gene, including SNPs mapping 1 Kb up-stream and down-stream the 5′ and the 3′ of the gene. We selected 18 independent tagging SNPs in high linkage disequilibrium (r2 ≥0.80) with 18 neighboring SNPs (Table 1). All SNPs complied with Hardy–Weinberg equilibrium. The 18 selected SNPs covered the entire gene and the 3′ and 5′ flanking regions.
The genotyped sample underwent quality control in accordance with the protocol written by Anderson et al. . The degree of recent shared ancestry for a pair of individuals was estimated using genome-wide identity by descent (IBD) as implemented in PLINK software v1.07 . We assessed population stratification with principal component analysis (PCA) as implemented in SNP and Variation Suite v8.x (Golden Helix, Inc., Bozeman http://www.goldenhelix.com) . We selected the first 10 PCs to include as covariates in the linear regression model.
We performed a quantitative trait interaction analysis (gene per environment analysis, GxE) testing PON1 SNPs association with HDL, LDL, Cholesterol, Triglycerides and AIP phenotypes in anthocyanins and polyphenols intake environments using PLINK v1.07 . The two environmental groups (high and low intake) were defined using the two extreme tertiles of anthocyanins and total polyphenol intake distributions in order to maximize the differences between the two groups. Quantitative traits were analysed as residuals, adjusted for sex, age, BMI and the first 10 principal components calculated using R software [31–35]. For each trait, individuals with missing data were excluded from the analysis. To correct for multiple testing we defined a Bonferroni threshold of p = 0.0028. Furthermore, in order to identify the genotypes that are associated to CVD protection (high HDL, low non-HDL cholesterol, triglycerides and AIP) in high anthocyanins/polyphenols intake, we performed a multivariate regression analysis in high and low environment separately using STATA (Stata corp.2015. Stata Statistical Software: Release 14. College Station, TX: StataCorp LP). Analyses were adjusted for BMI, age, gender and principal components.
Four hundred and forty-three subjects, 175 males and 268 females, were fully eligible for the analysis after quality control of genetic and phenotypic data (Fig. 1).
Among the 500 enrolled participants, 29 subjects with less than three dietary recall interviews were excluded, leaving 471 subjects for genotyping and downstream analysis.
After genotyping, one subject was excluded for low call rate (≤0.95) and eight subjects were filtered out for a reduced proportion of heterozygosity. We then removed 15 related individuals and four genetic outliers after PCA, defined as individuals that exceed six standard deviations from the whole sample along any of the principal components (Additional file 1: Figure S1).
Anthropometric characteristics of the overall sample and of each sub-group identified according to anthocyanins or polyphenols intake are described in Table 2. Anthocyanins intake subgroup included 283 individuals while polyphenols subgroup included 295 individuals. Subjects with missing biochemical data were excluded from each specific analysis (Table 2).
Daily anthocyanins intake ranged between 0 and 5.98 mg/day in the lower and 25.7 and 614.4 mg/day in the higher tertile. Polyphenols intake ranged between 99.4 and 804.5 mg/day in the lower and between 1288.0 and 4342.2 mg/day in the higher tertile.
High and low sub-groups were significantly different for age in anthocyanins and for sex (Chi squared test p < 0.05) and BMI in polyphenols environment (t test p < 0.05) (Table 2). BMI did not significantly correlate with anthocyanins intake (r = −0.05, p = 0.3) whereas it weakly correlated with polyphenols intake (r = −0.2, p < 0.0001) (Additional file 2: Figure S2). Interestingly, without considering the genetic background, no significant difference was observed for HDL, LDL, cholesterol, triglycerides and AIP between high and low anthocyanins and polyphenols intake. Evaluating the frequency of subjects out of the physiological ranges for each lipid parameter, we observed that only the percentage of subjects with HDL <45 in females and <35 in males is significantly higher in the low polyphenols intake sub-group compared to the high sub-group (odds ratio 4.4, CI 95 % 1.6–12.1, p = 0.004 for low intake) (Table 2).
Quantitative trait interaction results for the 18 selected SNPs in PON1 and CVD biomarkers, considering anthocyanins and polyphenols intake as environment, are presented in Table 3. According to the number of analyzed SNPs, the Bonferroni threshold is p = 0.0028. A significant p value indicates a difference in association with the lipid profiles between the two environments. The interaction analysis for HDL showed that 2 SNPs had a significant p value in interaction with anthocyanins (rs854549 flanking the gene at the 3′, p = 0.0008 and rs854552, at the 3′ untranslated region (3′ UTR) p = 0.001) and 2 SNPs in interaction with polyphenols (rs854571, p = 0.0021 and rs854572 p = 0.0020, both in the promoter region of the gene). A significant interaction model was found for rs854551, located at 3′ UTR (p = 0.0022) in relation to AIP and anthocyanins and for the intronic rs3917477 (p = 0.0026) in relation to total cholesterol and polyphenols.
Among the significant GxE interactions with HDL, cholesterol and AIP, we then selected the genotypes associated to CVD protection in the high anthocyanins and polyphenols intake. Figure 2 shows the trend of lipid profiles across the three genotypes for each significant GxE SNP in low or high environment. In particular, in high anthocyanins intake, HDL concentration was higher in carriers of the major CC genotype for rs854549 (p = 0.001, Beta = 4.7) and in the minor CC genotype for rs854552 (p = 0.001, Beta = 5.6). AIP was lower in minor AA carriers for rs854551 (p = 0.034, Beta = −0.07). Considering high polyphenols intake, HDL concentration was higher in minor TT carriers for rs854571 (p = 0.026, Beta = 3.92) and higher in minor CC compared to GG for rs854572 (p = 0.025, Beta = 3.94). Rs3917477 was not significantly associated to cholesterol in high polyphenols intake.
Within the ATHENA project, we performed a nutrigenetic observational study to determine whether SNPs that describe the genetic variability in PON1 gene can influence the response of cardiovascular health biomarkers to polyphenols and anthocyanins. We collected genetic, dietary, environmental, lifestyle data and laboratory measurements in 443 healthy Italians. As protective biomarkers of cardiovascular health we considered high HDL, low total cholesterol, LDL, triglycerides and AIP [14, 25].
In our analysis, high and low antioxidant intakes did not exert any beneficial effect on the target phenotypes if the genetic background related to PON1 gene was not considered. These results are in line with previous interventional studies that showed discordant findings on the effect of anthocyanins on common biomarkers of CVD . On the contrary, using a nutrigenetic approach, we could identify 5 SNPs significant at Bonferroni level (rs854549, rs854551, rs854552, rs854571, rs854572) and for each SNP we pointed out the genotype with a significant cardiovascular protective effect under high antioxidants intake.
In high anthocyanins intake, carriers of the C protective allele at rs854549 experienced an increase in HDL levels of 4.7 mg/dl (p = 0.001) while carriers of the C protective allele at rs854552 showed an increase of 5.6 mg/dl (p value 0.001). Considering high polyphenols intake, HDL levels were 3.92 mg/dl higher in T carriers for rs854571 (p = 0.026) and 3.94 mg/dl higher in C carriers for rs854572 (p = 0.025). These effects are of note if we consider that Boes  estimated that an increase of 1 mg/dl of HDL levels is associated with a 2 and 3 % reduction of the risk for coronary artery disease in men and women, respectively. AIP was lower in A carriers for rs854551, with a decrease of 0.07 (p = 0.034) in high anthocyanins intake.
PON1 gene is associated with several human diseases, related to oxidative stress including cardiovascular disease, Parkinson’s disease and cancer  and is inversely associated to the risk of CVD, particularly to atherosclerosis .
PON1 enzyme is tightly associated with HDL particles and protects both LDL and HDL from oxidation, a major step in the progression of atherosclerosis, the underlying pathophysiologic factor for the majority of cardiovascular diseases [37–39]. HDL contributes to PON1 enzyme stabilization, furnishes a hydrophobic environment that could be important for PON1 function and is a key player in the reverse cholesterol transport, which shuttles cholesterol from peripheral cells (e.g. macrophages) to the liver or other tissues.
As lifestyle determinants such as smoking, alcohol intake and exogenous or endogenous oxidants can modify PON1 levels and activity, several strategies were used to test if antioxidant supplementation, including polyphenols and anthocyanins, could improve PON1 function. It has been demonstrated that anthocyanins and polyphenols promote antioxidant activity and cholesterol efflux capacity of HDL. They also enhance PON1 stabilization, its association with HDL and catalytic activity [15, 16, 40].
Huen et al. reported rs854551 and rs854552, both mapping at the 3′ UTR of the gene, as significantly associated to paraoxonase activity of PON1 in a study on Mexican Newborn and Mothers .
Rs854571 and rs854572 map in the promoter region of PON1 gene and have been independently reported to produce an approximatively two-fold change in PON1 expression levels in human hepatoma cell line HepG2 [38, 43]. Leviev  reported a significant increase in activity of PON1 promoter related to the T allele of rs854571, the same allele that we found associated to higher level of HDL in high anthocyanins intake. Brophy  showed a significant increase in activity of PON1 promoter related to the G allele of rs854572, that in our sample was associated to lower HDL levels in high anthocyanins intake. A recent genome-wide association analysis demonstrated that PON1 polymorphisms are strongly associated to PON1 function, and especially rs854572 is the best predictor of arylesterase activity .
On the contrary, associations between PON1 SNPs and lipid profiles have been mostly controversial [21–23, 36, 44–46]. Nus et al.  evaluated the effect of walnut enriched meat consumption on lipid profile, in subjects at increased risk of CVD, according to rs662 and rs854560 genotypes. They observed that rs662 TT carriers had lower HDL, LDL and triglycerides levels compared to C carriers and that the difference increased if the SNPs were considered in interaction (rs662 TT + rs854560 TT + AT versus the other genotypes). However, the lipid profile did not differ between rs662 and rs854560 after consumption of either walnut enriched or control meat diet. De Souza et al. in a pharmacogenomics study on statin response in Brazilians, found that T allele carriers of rs854560 reached their HDL lipid target (HDL >1.55 mmol/L) more often than patients who were homozygous AA. Rios et al. . showed rs662 as associated to triglycerides and HDL levels in male Caucasian-Brazilians. Van Himbergen et al.  found that genetic variants in PON1 associated with high levels and activity of the enzyme, were also associated with higher HDL levels. That study, however, was performed in patients with familial hypercholesterolemia with a high prevalence of atherosclerosis in the population. Van Aalst-Cohen et al.  in a familial hypercholesterolemia study and Blatter Garin et al.  in a case control study on coronary artery disease could not find any association between the two coding polymorphisms rs854560 and rs662 and HDL level. In our study, rs854560 was significantly associated to HDL in interaction with anthocyanins (p = 0.0083) but the significance was not confirmed after correction for multiple testing. There was a tendency for TT homozygous to higher HDL levels in high anthocyanins intake. Rs662 was not associated to any of the phenotype tested.
Though a limitation of our study could be the small sample size, its strength is an accurately phenotyped cohort that includes subjects with well-established long-term dietary habits regarding antioxidant intake. To address bias due to self-reported data from the nutritional interview, trained dieticians conducted the 24 h dietary recall using pictures (i.e. validated Food Atlas) to help people identifying food portions. In addition, as seasonality affects anthocyanins levels that are higher in summer and lower in winter, we repeated the 24 h recall interview at least three times a year and used the average yearly intake for all the analyses [47, 48]. This has increased the precision in evaluating the chronic effect of the diet on lipid profile. Moreover, the DMS software, integrating different databases, returns a detailed estimate of the anthocyanins and polyphenols content for all the recorded foods. A validation on an independent sample would allow confirming causality of the SNPs identified in the present study.
This nutrigenetic study highlights the role of PON1 as susceptibility gene for cardiovascular health under high antioxidant intake. We identified protective genotypes in 4 independent polymorphisms, presenting a significant association with increased HDL level under high polyphenols and anthocyanins intake. The identified alleles could be used to stratify subjects who can benefit from high antioxidant consumption to improve their HDL levels.
atherogenic index of plasma
anthocyanin and polyphenols bioactive for health enhancement through nutritional advancement
Banca Dati Alimentare-Istituto Europeo Oncologico
diet monitoring solution
gene per environment analysis
human hepatoma cell line
identity by descent
minor allele frequency
principal component analysis
single nucleotide polymorphism
Fuster V. Global burden of cardiovascular disease: time to implement feasible strategies and to monitor results. J Am Coll Cardiol. 2014;64:520–2.
Deaton C, Froelicher ES, Wu LH, Ho C, Shishani K, Jaarsma T. The global burden of cardiovascular disease. Eur J Cardiovasc Nurs. 2011;10(Suppl 2):S5–13.
Daar AS, Singer PA, Persad DL, Pramming SK, Matthews DR, Beaglehole R, et al. Grand challenges in chronic non-communicable diseases. Nature. 2007;450:494–6.
Niebylski ML, Redburn KA, Duhaney T, Campbell NR. Healthy food subsidies and unhealthy food taxation: a systematic review of the evidence. Nutrition. 2015;31:787–95.
de Roos B. Personalised nutrition: ready for practice? Proc Nutr Soc. 2013;72:48–52.
Kim DS, Marsillach J, Furlong CE, Jarvik GP. Pharmacogenetics of paraoxonase activity: elucidating the role of high-density lipoprotein in disease. Pharmacogenomics. 2013;14:1495–515.
Robitaille J, Houde A, Lemieux S, Gaudet D, Pérusse L, Vohl MC. The lipoprotein/lipid profile is modulated by a gene-diet interaction effect between polymorphisms in the liver X receptor-alpha and dietary cholesterol intake in French-Canadians. Br J Nutr. 2007;97:11–8.
Jennings A, Welch AA, Fairweather-Tait SJ, Kay C, Minihane A-M, Chowienczyk P, et al. Higher anthocyanin intake is associated with lower arterial stiffness and central blood pressure in women. Am J Clin Nutr. 2012;96:781–8.
McCullough ML, Peterson JJ, Patel R, Jacques PF, Shah R, Dwyer JT. Flavonoid intake and cardiovascular disease mortality in a prospective cohort of US adults. Am J Clin Nutr. 2012;95:454–64.
Cassidy A, O’Reilly ÉJ, Kay C, Sampson L, Franz M, Forman JP, et al. Habitual intake of flavonoid subclasses and incident hypertension in adults. Am J Clin Nutr. 2011;93:338–47.
Du G, Sun L, Zhao R, Du L, Song J, Zhang L, et al. Polyphenols: potential source of drugs for the treatment of ischaemic heart disease. Pharmacol Ther. 2016;162:23–34.
Qin Y, Xia M, Ma J, Hao Y, Liu J. Anthocyanin supplementation improves serum LDL-and HDL-cholesterol concentrations associated with the inhibition of cholesteryl ester transfer protein in dyslipidemic subjects. Am J Clin Nutr. 2009;90:485–92.
Wallace TC. Anthocyanins in cardiovascular disease 1. Adv Nutr. 2011;2:1–7.
Wallace T, Slavin M, Frankenfeld C. Systematic review of anthocyanins and markers of cardiovascular disease. Nutrients. 2016;8:32.
Zhu Y, Huang X, Zhang Y, Wang Y, Liu Y, Sun R, et al. Anthocyanin supplementation improves HDL-Associated paraoxonase 1 activity and enhances cholesterol efflux capacity in subjects with hypercholesterolemia. J Clin Endocrinol Metab. 2014;99:561–9.
Farrell N, Norris G, Lee SG, Chun OK, Blesso CN. Anthocyanin-rich black elderberry extract improves markers of HDL function and reduces aortic cholesterol in hyperlipidemic mice. Food Funct. 2015;6:1278–87.
Costa LG, Giordano G, Furlong CE. NIH Public Access. 2012;81:337–44.
Aviram M, Rosenblat M, Gaitini D, Nitecki S, Hoffman A, Dornfeld L, et al. Pomegranate juice consumption for 3 years by patients with carotid artery stenosis reduces common carotid intima-media thickness, blood pressure and LDL oxidation. Clin Nutr. 2004;23:423–33.
Aviram M, Rosenblat M. Pomegranate for your cardiovascular health. Rambam Maimonides Med J. 2013;4:e0013.
Tang WHW, Hartiala J, Fan Y, Wu Y, Stewart AFR, Erdmann J, et al. Clinical and genetic association of serum paraoxonase and arylesterase activities with cardiovascular risk. Arterioscler Thromb Vasc Biol. 2012;32:2803–12.
Blatter Garin MC, Moren X, James RW. Paraoxonase-1 and serum concentrations of HDL-cholesterol and apoA-I. J Lipid Res. 2006;47:515–20.
de Souza JA, Menin A, Lima LO, Smiderle L, Hutz MH, Van Der Sand CR, et al. PON1 polymorphisms are predictors of ability to attain HDL-C goals in statin-treated patients. Clin Biochem. 2015;48:1039–44.
Nus M, Frances F, Librelotto J, Canales A, Corella D, Sanchez-Montero JM, et al. Arylesterase Activity and antioxidant status depend on PON1-Q192R and PON1-L55M Polymorphisms in subjects with increased risk of cardiovascular disease consuming walnut-enriched meat. J Nutr. 2007;137:1783–8.
Conti C, Rossi E, Marceglia S, Tauro V, Rizzi F, Lazzaroni M, et al. An integrated diet monitoring solution for nutrigenomic research. Stud Health Technol Inform. 2015;210:632–6.
Dobiásová M, Frohlich J, Sedová M, Cheung MC, Brown BG. Cholesterol esterification and atherogenic index of plasma correlate with lipoprotein size and findings on coronary angiography. J Lipid Res. 2011;52:566–71.
Delaneau O, Marchini J, Zagury J-F. A linear complexity phasing method for thousands of genomes. Nat Methods. 2012;9:179–81.
Li Y, Willer CJ, Ding J, Scheet P, Abecasis GR. MaCH: using sequence and genotype data to estimate haplotypes and unobserved genotypes. Genet Epidemiol. 2010;34:816–34.
Anderson CA, Pettersson FH, Clarke GM, Cardon LR, Morris AP, Zondervan KT. Data quality control in genetic case-control association studies. Nat Protoc. 2010;5:1564–73.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–75.
Bozeman MT. Golden Helix, Inc. SNP and Variation Suite™ (Version 8.x) [Software]. http://www.goldenhelix.com. Accessed Sept 2015.
R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2014.
Sarkar D. Lattice. New York: Springer; 2008.
Lenth RV. Least-squares means: the R Package lsmeans. J Stat Softw. 2016;69(1):1–33.
Auguie B. gridExtra: miscellaneous functions for “Grid” graphics. R package version 2.2.1. 2016. https://CRAN.R-project.org/package=gridExtra. Accessed Jan 2016.
Aulchenko YS, Ripke S, Isaacs A, van Duijn CM. GenABEL: an R library for genome-wide association analysis. Bioinformatics. 2007;23:1294–6.
Boes E, Coassin S, Kollerits B, Heid IM, Kronenberg F. Cholesterol levels: a systematic in-depth review. Genet Epidemiol. 2010;44:136–60.
Rosenblat M, Vaya J, Shih D, Aviram M. Paraoxonase 1 (PON1) enhances HDL-mediated macrophage cholesterol efflux via the ABCA1 transporter in association with increased HDL binding to the cells: a possible role for lysophosphatidylcholine. Atherosclerosis. 2005;179:69–77.
Brophy VH, Hastings MD, Clendenning JB, Richter RJ, Jarvik GP, Furlong CE. Polymorphisms in the human paraoxonase (PON1) promoter. Pharmacogenetics. 2001;11:77–84.
Sahebkar A, Simental-Mendía LE, Giorgini P, Ferri C, Grassi D. Lipid profile changes after pomegranate consumption: a systematic review and meta-analysis of randomized controlled trials. Phytomedicine. 2016;16:2–7.
Rock W, Rosenblat M, Miller-Lotan R, Levy AP, Elias M, Aviram M. Consumption of wonderful variety pomegranate juice and extract by diabetic patients increases paraoxonase 1 association with high-density lipoprotein and stimulates its catalytic activities. J Agric Food Chem. 2008;56:8704–13.
Carlson CS, Heagerty PJ, Hatsukami TS, Richter RJ, Ranchalis J, Lewis J, et al. TagSNP analyses of the PON gene cluster: effects on PON1 activity, LDL oxidative susceptibility, and vascular disease. J Lipid Res. 2006;47:1014–24.
Huen K, Barcellos L, Beckman K, Rose S, Eskenazi B, Holland N. Effects of PON polymorphisms and haplotypes on molecular phenotype in Mexican-American mothers and children. Environ Mol Mutagen. 2011;52:105–16.
Leviev I, James RW. Promoter polymorphisms of human paraoxonase PON1 gene and serum paraoxonase activities and concentrations. Arterioscler Thromb Vasc Biol. 2000;20:516–21.
Rios DLS, D’Onofrio LO, Cerqueira CCS, Bonfim-Silva R, Carvalho HG, Santos-Filho A, et al. Paraoxonase 1 gene polymorphisms in angiographically assessed coronary artery disease: evidence for gender interaction among Brazilians. Clin Chem Lab Med. 2007;45:874–8.
van Himbergen TM, Roest M, de Graaf J, Jansen EHJM, Hattori H, Kastelein JJP, et al. Indications that paraoxonase-1 contributes to plasma high density lipoprotein levels in familial hypercholesterolemia. J Lipid Res. 2005;46:445–51.
van Aalst-Cohen ES, Jansen ACM, Boekholdt SM, Tanck MWT, Fontecha MR, Cheng S, et al. Genetic determinants of plasma HDL-cholesterol levels in familial hypercholesterolemia. Eur J Hum Genet. 2005;13:1137–42.
Kassim A, Poette J, Paterson A, Zait D, McCallum S, Woodhead M, et al. Environmental and seasonal influences on red raspberry anthocyanin antioxidant contents and identification of quantitative traits loci (QTL). Mol Nutr Food Res. 2009;53:625–34.
Malcher ET, Carvalho JCT. The influence of seasonality on the anthocyanin concentrations in the açaí fruit (Euterpe oleracea Mart.) from the brazilian amazon. Int J Plant Anim Env Sci. 2011;1:224–32.
The authors’ responsibilities were as follows: CB, CC, DC, FR, ES: developed the overall research plan; MT and FLR: were involved in managing the observational study and collecting clinical and nutritional data at IRCCS Galeazzi Orthopedic Institute, Milano, Italy; ED, AT, AM, CAB, FP, LS were involved in managing the observational study and collecting clinical and nutritional data at ASST Santi Paolo e Carlo, Division of Nephrology, Milano, Italy; ML conducted the biochemical sample analysis; DB, FD performed genotyping; CC and MB managed the clinical and nutritional database and performed the phenotype quality control; SL, ES, FR, MC performed the statistical analysis; FR, CC, CB, and ES wrote the manuscript. All authors read and approved the final manuscript.
We thank Klizia Cavalletti, Amina Remli, Sara Palladino, Sara Paleologo, Giuseppe Beretta, Laura Scaroni, Emmanuele Giglio for their help with the clinical data collection.
The authors declare that they have no competing interests.
Availability of data and materials
Additional analyses are ongoing on other possible candidate genes, so the dataset is not available.
The owner of DMS software is a private company, KOS Genetic S.R.L, Milano and it is not available for sharing. Detailed information on software characteristics can be found in the paper by Conti et al. .
Ethics approval and consent to participate
The ethic committees of ASST Santi Paolo e Carlo and IRCCS Galeazzi Orthopedic Institute approved the study (Register 441 signed off on 9 June 2011 and Register 45-1250-ATHENA signed off on 22 February 2013, respectively), in accordance with principles of the Declaration of Helsinki and all participants gave written informed consent.
European Union, Grant number FP7-KBBE-2009-3 (ATHENA) and MIUR-CNR PB05 InterOmics Italian Flagship Project, supported this work. Elena Dogliotti was supported by a Fondazione Umberto Veronesi Fellowship.
Federica Rizzi and Costanza Conti are joint first Authors