PDF Version 1.5 MB
Background: Paraoxonase 1 (PON1) is an enzyme that detoxifies activated organophosphorus pesticides (OPs) and is also involved in oxidative stress pathways.
Objectives: PON1 activity in newborns is lower than in adults, but the ontogeny of PON1 activity is poorly characterized in young children. We examined the effects of age and PON1 genotype on enzyme activity in a birth cohort of Mexican-American children.
Methods: We determined three substrate-specific measures of PON1 activity in 1,143 plasma samples collected longitudinally from 458 children at five time points from birth through 7 years of age, and genotyped PON1 polymorphisms at positions 192 and –108 in these children.
Results: Contrary to previous reports that PON1 activities plateau by 2 years of age, we observed an age-dependent increase in all three PON1 measures from birth through 7 years of age (p < 0.0001). The PON1192 genotype significantly modified the effect of age on paraoxonase (POase) activity (p < 0.0001) such that increases in enzyme activity with age were influenced by the number of R alleles in a dose-dependent manner. Children with the PON1-108CC192RR diplotype had significantly higher mean PON1 activities and also experienced steeper increases of POase activity over time compared with children with the PON1-108TT192QQ diplotype.
Conclusions: Lower levels of the PON1 enzyme, which is involved in protection against OPs and oxidative stress, persist in young children past 2 years of age through at least 7 years of age. Future policies addressing pesticide exposure in children should take into account that the window of vulnerability to OPs in young children may last beyond infancy.
Citation: Huen K, Harley K, Brooks J, Hubbard A, Bradman A, Eskenazi B, Holland N. 2009. Developmental Changes in PON1 Enzyme Activity in Young Children and Effects of PON1 Polymorphisms. Environ Health Perspect 117:1632–1638; http://dx.doi.org/10.1289/ehp.0900870
Address correspondence to N. Holland, 50 University Hall, School of Public Health, UC Berkeley, Berkeley, CA 94720-7360 USA. Telephone: (510) 455-0561. Fax: (510) 643-5426. E-mail: firstname.lastname@example.org
We gratefully acknowledge Center for Health Assessment in Mothers and Children of Salinas (CHAMACOS) staff, community partners, and especially the CHAMACOS participants. We appreciate the insightful comments of C. Furlong and the assistance with genotyping from K. Beckman.
This publication was made possible by grants from the U.S. Environmental Protection Agency (EPA) (R826886, R82670901) and the National Institute of Environmental Health Science (NIEHS) (R01ESO12503-03, P01 ES009605). The contents of this article are solely the responsibility of the authors and do not necessarily represent the official views of the NIEHS or the U.S. EPA.
The authors declare they have no competing financial interests.
Received: 3 April 2009
Accepted: 9 June 2009
Advance Publication: 9 June 2009
Final Publication: 1 October 2009
Supplemental Material 188 KB
The enzyme paraoxonase 1 (PON1; EC 220.127.116.11) detoxifies activated (oxon) organophosphorus pesticides (OPs), which are neurotoxic (Costa et al. 2005a; Furlong et al. 2002). PON1 also is an antioxidant, and several studies indicate that it inhibits oxidation of low-density lipoprotein, a marker of oxidative stress (James 2006; Li et al. 2003). PON1 genotypes and enzyme activities have been associated with numerous health outcomes related to OP exposure and oxidative stress, including Alzheimer’s disease (Erlich et al. 2006), Parkinson’s disease (Zintzaras and Hadjigeorgiou 2004), cardiovascular disease (Bhattacharyya et al. 2008), and preterm birth (Chen et al. 2004; Lawlor et al. 2006). Studies of in utero OP exposures have also demonstrated potential gene–environment interactions with PON1. Berkowitz et al. (2004) found that among pregnant women with measureable levels of chlorpyrifos in their blood, those with lower PON1 levels had children with smaller infant head circumferences at birth. The same research group also reported that higher maternal OP exposure (measured by urinary metabolites) was significantly associated with decreased child birth weight among mothers with the PON1192QQ genotype (Wolff et al. 2007). PON1 enzyme activities and levels vary broadly in humans (Costa et al. 2005b), and the determinants of PON1 variation, including genetics and age, may further explain the role this enzyme plays in relation to exposures and diseases.
Several common polymorphisms in the coding and promoter regions of the PON1 gene influence substrate-specific PON1 enzyme activities (Brophy et al. 2001; Ferre et al. 2003), particularly single-nucleotide polymorphisms (SNPs) at positions 192 and –108 (James et al. 2005). The SNP at codon 192 results in an amino acid substitution from the active-site residue glutamine (Q) to arginine (R), and the catalytic efficiency of the PON1192R alloform toward the oxon derivatives of OP pesticides parathion and chlorpyrifos is greater than that of the PON1192Q alloform. Recent studies have demonstrated that the PON1192 genotype accounts for a large portion of the variation of PON1 activity toward the OP substrate paraoxon in Caucasian and African-American (59%) adults (Bhattacharyya et al. 2008) and also in a Mexican-American (48%) population (Rainwater et al. 2008). Of the promoter polymorphisms, the PON1-108 SNP exerts the most noticeable effect on arylesterase (AREase) activity, a measure of PON1 quantity (22.4%), such that the C allele is associated with increased PON1 levels compared with the T allele (Brophy et al. 2001; Deakin et al. 2003). Although genetic polymorphisms account for a large portion of PON1 variation, it is not sufficient in epidemiologic studies to consider PON1 genotypes alone (Jarvik et al. 2000; Richter and Furlong 1999). Interindividual variability of PON1 phenotypes can range widely even among individuals with the same PON1 genotypes because enzyme levels can also be different. Therefore, studies that measure PON1 activities may be more informative than studies that rely solely on PON1 genotype data.
Several animal studies have demonstrated that low PON1 activity during development results in a reduced metabolic capacity in the young and therefore an increased susceptibility to OP toxicity (Moser et al. 1998; Timchalk et al. 2007). Similarly in humans, PON1 levels in infants are much lower than levels in adults (Chen et al. 2003; Holland et al. 2006). Furthermore, preterm babies have 24% lower activities compared with term babies (Ecobichon and Stephens 1973), suggesting that PON1 levels in the fetus may be even lower. A few studies have measured PON1 activities in healthy young children at single points in time (Garces et al. 2008; Karikas et al. 2006; Sumegova et al. 2007; Vincent-Viry et al. 2003), but only one study has described the ontogenetic profile of PON1 activity in the same individuals over time. This study (Cole et al. 2003) of nine children demonstrated high inter-individual variation of PON1 from birth to 2 years of age; PON1 activities increased over time and reached a plateau between 6 months and 2 years of age in most, but not all, individuals. PON1 activity is relatively stable in adults (Costa et al. 2005b); however, in the elderly, it decreases steadily over time (Seres et al. 2004). To date, longitudinal changes in PON1 activity in children after 2 years of age and potential differences between genotypes have not yet been characterized.
Previously, we reported a relatively high level of OP exposure in the Salinas Valley, California, the location of our longitudinal birth cohort (Bradman et al. 2005, 2007). Additionally, we observed associations between in utero OP exposure and decreases in gestational duration (Eskenazi et al. 2004), abnormal reflexes in neonates (Young et al. 2005), poorer mental development at 2 years of age (Eskenazi et al. 2008), and increases in maternal report of pervasive development disorder (Eskenazi et al. 2008) in these children. Initially, we reported PON1 activity in 130 mothers and their newborns and demonstrated much lower activities in the newborns (Furlong et al. 2006; Holland et al. 2006). In the present study, we measured PON1 activity in 458 children at several time points spanning birth through 7 years of age to determine the ontogeny of PON1 phenotypes. In the same children, we also genotyped two PON1 SNPs (–108 and 192) known to influence phenotypic variation and examined how PON1 genotypes affect children’s PON1 activities at different ages.
Materials and Methods
Study subjects. The Center for Health Assessment in Mothers and Children of Salinas (CHAMACOS) study is a longitudinal birth cohort study of the effects of pesticide and other environmental exposures on neurodevelopment, growth, and respiratory disease in children from primarily Mexican-American families (Eskenazi et al. 2003). Located in Monterey County, California, the Salinas Valley is an area of heavy agricultural production where approximately 200,000 kg of OPs are applied annually (Department of Pesticide Regulation 2005). Women and children in this community are at risk for higher exposures to pesticides through direct exposure if they live near agricultural fields or if the mother works in the field, and through the take-home pathway. A total of 601 pregnant women were enrolled in the CHAMACOS study, and 528 delivered newborns. Mothers in the CHAMACOS cohort were primarily young (mean ± SD, 25.6 ± 5.3 years), married, low-income, Mexican-born, Spanish-speaking women who were farmworkers themselves (44%) or lived with farmworkers at the time of enrollment (84%). Ethnicity of children and mothers was based on mothers’ self-report. For this particular study, we limited analyses to include only subjects who were of Mexican descent (97%). More than half the women in CHAMACOS came from three states in Mexico: Michoacán (23%), Guanajuato (21%), and Jalisco (11%). Study protocols were approved by the University of California–Berkeley human subjects review committee. Written informed consent was obtained from all mothers, which included consent for their children to participate.
Blood collection and processing. Blood specimens were collected from 458 CHAMACOS children (50.1% boys) for the measurement of PON1 enzyme activity. Umbilical cord blood samples (n = 336) were collected by delivery room staff once the baby was safely delivered. Blood specimens were also collected from children at approximately 1 year of age (n = 52; mean ± SD, 1.07 ± 0.19 years), 2 years of age (n = 249; 2.01 ± 0.09 years), 5 years of age (n = 225; 5.10 ± 0.22 years), and 7 years of age (n = 281; 7.10 ± 0.23 years). We were able to measure PON1 enzyme activity at four or more time points for 108 of these children.
Heparinized whole blood was collected in BD Vacutainers (Becton, Dickinson and Company, Franklin Lakes, NJ), centrifuged, and then divided into plasma, buffy coats, and red blood cells and stored at −80°C. Vacutainers without anticoagulant were used to collect serum and clot. DNA was isolated from blood clots as described previously (Holland et al. 2006).
Determination of PON1 levels and activities. We measured PON1 enzyme activity toward three substrates [paraoxon (PO), phenyl acetate (aryl ester; ARE), and chlorpyrifos-oxon (CPO)] in plasma samples using spectrophotometric methods as described previously (Richter and Furlong 1999). The AREase assay is considered a measure of PON1 enzyme quantity because rates of phenyl acetate hydrolysis do not differ between PON1192 Q and R alloforms as they do for paraoxon hydrolysis (POase). Methods based on ELISA and Western blot using PON1 antibodies confirm a high correlation (r > 0.85) between measured PON1 quantity and AREase activity (Connelly et al. 2008; Kujiraoka et al. 2000). In contrast, POase and chlorpyrifos-oxonase (CPOase) substrate-specific assays reflect both quantity and catalytic efficiency of the enzyme. In this article, we refer to enzyme quantity measured by AREase activity as PON1 levels. CPOase and POase activities (quantity and substrate-specific efficiency) are referred to as enzyme activity. All assays were performed in triplicate. Quality assurance, described in more detail by Huen et al. (2009), included assessment of repeat samples (separate aliquots of the same sample run on different days), internal controls (aliquots of the same sample run on all assay plates), and concurrent analyses of specimens from different time points (samples from different time points run on the same plates). Repeated analysis of 3% of samples showed a high degree of concordance. For the three substrate-specific assays, the average coefficient of variation (CV) for repeated samples ranged from 6% to 9%, and the correlation coefficients between repeated runs ranged from 0.91 to 0.98. Interassay variability, as measured by the average CV for internal controls samples, was between 7 and 9%.
Determination of PON1 genotypes. The coding polymorphism, PON1192, was genotyped using the TaqMan real-time polymerase chain reaction (PCR) method. Briefly, primers for the nucleotide sequence flanking the SNP and probes specific for the SNP were custom designed by Applied Biosystems, Inc. (Foster City, CA). The promoter SNP, PON1-108, was genotyped using a fluorogenic allele-specific genotyping assay (Amplifluor; Flowgen Biosciences Ltd., Nottingham, UK). This assay required a two-part nested PCR strategy, where the region surrounding the SNP was preamplified using nonallelic flanking primers. The amplicon was then diluted and used as the template for the Amplifluor assay. Quality assurance procedures for genotyping both PON1 SNPs included assessment of randomly distributed blank samples in each plate and duplicates of randomly selected samples with independently isolated DNA from the same subjects. Repeated analysis (4% of samples) in several runs showed a high degree (> 99%) of concordance. All discrepancies were resolved with additional genotyping.
Statistical analysis. To model longitudinal changes of PON1 enzyme activity, we used generalized estimating equations (GEEs) with an exchangeable correlation structure because we had repeated measures on individual subjects. This method is often advantageous in longitudinal models because inferences based on robust estimators are less dependent on the assumption that the data are normally distributed (Hox 2002). Correlations of measurements on individuals at various ages were similar. We created separate models for each of the three outcomes, AREase, CPOase, and POase activity. Independent variables in the model included age (years), PON1 genotype (either PON1192 or PON1-108), and an inter-action term between age and PON1 genotype (age × PON1192 or age × PON1-108). Because the number of specimens collected at 1 year of age was much smaller than at other time points, we also ran GEE models without this time point. However, because results did not change significantly, we show only the GEE models using all time points.
We previously established that storage duration can affect PON1 activity, particularly in specimens that had been stored for > 2 years (Huen et al. 2009). PON1 activity decreased on average by 17.1%, 39.4%, and 37.6% for AREase, CPOase, and POase, respectively after 5 years of storage (p < 0.001 for all three assays). This decline of PON1 enzymatic activity over time has also been reported in other studies (Brackley et al. 1983; Stenzel et al. 2007). In this study, specimens collected at earlier ages were sometimes stored for > 2 years, so age and storage duration were inversely correlated with each other (Pearson’s correlation coefficient r = 0.986). To estimate the effects of age while taking into account the potential effects of storage duration, we generated a storage duration correction factor using pilot data from a previous study (Huen et al. 2009). Briefly, parallel aliquots of specimens from 95 subjects were assayed after 2 years and 7 years of storage at −80°C. Measured AREase, CPOase, and POase activities were significantly lower after 7 years of storage compared with 2 years of storage. Using these data, we constructed linear models to predict the percent change in PON1 activity as a function of years of storage duration and storage duration × assay temperature. The models including an interaction between storage duration and assay temperature had the lowest Akaike Information Criterion and were thus selected to be used to determine the correction factors. The intercept was constrained to equal zero. We then used the resultant coefficients from these models as correction factors to predict the PON1 activity (AREase, CPOase, POase) at a storage time of zero years, and used the storage duration–corrected PON1 activity measures in subsequent GEE models to determine the effects of age on PON1 activity as described above.
Because we could not completely eliminate the residual confounding by storage duration in our GEE models, we conducted additional analyses to provide evidence that the association between age and PON1 substrate-specific measures is not attributable solely to the strong correlation of age with storage duration. We performed a likelihood ratio test comparing a full model, including terms for age, storage, PON1 genotype, and the interaction between age and genotype, with a restricted model without the age terms. It indicated that age should be retained in the model and that the fit of the full model was significantly better (p < 0.0001). Additionally, because the correlation between age and storage duration was considerably weaker within shorter time intervals (r = −0.3 between 1- and 2-year-olds and r = −0.93 between 5- and 7-year-olds), we also ran GEE models including a term for duration of storage in data sets restricted a) only to 1- and 2-year-olds and b) only to 5- and 7-year-olds. In the models for 1- and 2-year-olds, an interaction term for PON1 genotype and age was not included because there was not sufficient power with the limited sample size.
All analyses were performed in STATA (version 10.0; StataCorp, College Station, TX) and R (version 2.7.0; R Development Core Team 2008).
PON1 polymorphisms. Genotype distributions did not deviate significantly from Hardy–Weinberg equilibrium. For the PON1192 polymorphism, the allele frequencies of the Q and R alleles were 50.2% and 49.8%, respectively, and genotype frequencies were 23.2%, 53.1%, and 23.7% for QQ, QR, and RR genotypes (Table 1), respectively. Genotype frequencies for PON1-108 were 28.3%, 53.3%, and 18.4% for CC, CT, and TT, respectively. The allele frequency of the major allele C was 54.7%.
Effect of age on enzyme activity. Table 2 lists descriptive statistics for AREase, CPOase, and POase activity at birth and 1, 2, 5, and 7 years of age. Mean PON1 enzyme activity in CHAMACOS children was lowest in newborns for all three substrate-specific assays and highest at 7 years of age, and there were no differences in PON1 activity by sex (analysis of variance, F < 1.14; p > 0.30 for all three PON1 assays at all time points). Mean ± SD AREase activity, which is considered a measure of PON1 levels (quantity), approximately doubled from 33.1 ± 14.5 U/mL at birth to 66.5 ± 30.4 U/mL at 1 year of age and quadrupled (compared with birth) to 121.5 ± 30.3 U/mL at 7 years of age. Although PON1 activity did not increase appreciably between 2 and 5 years of age for AREase and CPOase, all three substrate specific measures grew noticeably between 5 and 7 years of age. At birth, differences in overall enzymatic activity between newborns with the lowest and highest values were 26-, 34-, and 136-fold for AREase, CPOase, and POase, respectively. In 7-year-old children, the fold differences were somewhat smaller (17-, 13-, and 21-fold for AREase, CPOase, and POase, respectively).
Effects of PON1 genotype on enzyme activity. Within each age group, the PON1 enzyme levels were influenced by both PON1192 and PON1-108. Mean AREase values differed by the PON1192 polymorphism (p for trend < 0.05 for increasing number of R alleles). At birth, average AREase levels were 15.8% higher in RR children than QQ children (Table 1). In contrast, at 2, 5, and 7 years of age, mean AREase activity was 9.8–14.1% higher in QQ children than in RR children. Mean AREase activity also varied by PON1-108 genotypes, and these differences were larger at all ages than for PON1192. At birth, mean AREase activity was 18.9 U/mL or about 44.9% lower in TT newborns than in CC newborns (Figure 1D). PON1 levels continued to be lower in CC children (p for trend < 0.001 for all time points except at year 1), with differences ranging from 18.8% to 25.4% in 2-, 5-, and 7-year-old children.
The PON1192 SNP explained a relatively small amount of the variance of AREase activity at different ages (1.4–3.6%), whereas the promoter polymorphism PON1-108 accounted for 8.2–24.4% of AREase activity [see Supplemental Material, Table 1 (available online at doi:10.1289/ehp.0900870.S1 via http://dx.doi.org)]. Although 4.4% of the variance of CPOase activity was explained by the PON1192 SNP at birth, the coding region SNP made a seemingly smaller contribution at other time points (0.1–1.2%), and the role of PON1-108 seemed much larger (4.1–27.3%). Both SNPs accounted for a significant portion of POase activity variation. The contribution of PON1192, which was 39.8% at birth, was generally higher in older children and highest at 7 years of age (70.7%).
Effects of age and PON1 genotype on enzyme activity. We used GEE augmentation of linear models to determine the effects of age and genotype on PON1 activity in children (Table 3). Point estimates in the models using storage duration–corrected data were not markedly different from those generated using raw data, although 95% confidence intervals (CIs) were slightly wider in storage duration–corrected models. Because the results for both models were similar, only results from the storage duration–corrected models are shown. Regression lines were generated for AREase, CPOase, and POase activity by age and are plotted with respect to PON1192 (Figure 1A–C) and PON1-108 (Figure 1D–F) genotypes. All three substrate-specific PON1 enzyme activities were positively associated with age (p < 0.0001 for all models). It appeared that PON1 activities continued to increase past the age of five because paired t-tests showed that AREase, CPOase, and POase were all significantly higher in 7-year-olds than in 5-year-olds (p < 0.0005).
For AREase activity, there was a small but statistically significant interaction (p = 0.016) between age and PON1192 genotype in which the slopes for age were lowest (8.3 U/mL per year) for RR children and highest (10.3 U/mL per year) for QQ children (Figure 1A). The promoter polymorphism, PON1-108, had a more pronounced effect on AREase activity than did PON1192. Furthermore, we observed a small but significant interaction (p = 0.027) between age and PON1-108 genotype in which mean AREase activity increased 10.0, 9.1, and 8.2 U/mL per year increase in age for CC, CT, and TT children, respectively.
Similar to AREase, CPOase activity was modestly associated (p = 0.017) with PON1192 genotype (Table 4, Figure 1B). Mean CPOase activity was 591.0 U/mL (95% CI, 104.6–1,077.4 U/mL) higher in RR children and 295.5 U/mL (95% CI, 52.3–538.7 U/mL) higher in QR children than in QQ children. There was no interaction between age and PON1192 genotype. In contrast, we did observe a small but significant interaction between age and PON1-108 genotype (β-estimate for age × PON1-108 genotype = −54.6, p = 0.022).
We also observed significant effect modification by the PON1192 genotype (p < 0.0001) on the relationship between age and POase activity (Figure 1C). POase activity increased with age for all PON1192 genotypes, but the slopes for age were much steeper for RR children. For example, mean activity increased by only 21.8 U/L (95% CI, 15.3–28.3 U/L) per year in QQ children, whereas in RR children mean activity increased 119.2 U/L or almost 5.5-fold more per year (95% CI, 97.7–140.9 U/L). The effect is so large that even at 7 years of age the average POase activity in QQ children was lower than the average POase activity in QR and RR children at 2 years of age. Although not as strikingly as the coding polymorphism PON1192, the promoter SNP, PON1-108, also influenced POase activity significantly (Figure 1F). The effect of age on POase activity was modified by PON1-108 genotype (p < 0.0001) such that slopes for mean POase activity per year were highest in CC children and lowest in TT children.
To determine the effects of combinations of PON1192 and PON1-108 genotypes on enzyme activity, we performed GEE models comparing children with the PON1-108CC192RR homozygous diplotype (reference group) with those with the PON1-108TT192QQ homozygous diplotype (Table 4). Children with the PON1-108TT192QQ diplotype had significantly lower mean AREase (14.1 U/mL; 95% CI, 3.7–24.4 U/mL), CPOase (1525.9 U/L; 95% CI, 835.0–2216.3 U/L), and POase (627.4 U/L; 95% CI, 534.5–720.4 U/L) activities compared with PON1-108CC192RR children. Additionally, we observed a statistically significant interaction between age and diplotype on POase activity (p < 0.0001). The slope for mean POase activity over age was much steeper for PON1-108CC192RR children (135.3 U/L per year; 95% CI, 115.0–155.7 U/L) than for PON1-108TT192QQ children, which was 112.6 U/L per year (95% CI, 91.3–133.8 U/L) lower.
Effects of age and PON1 genotype on enzyme activity in restricted data sets. Because the correlation between age and storage duration was appreciably smaller within shorter time intervals, we were able to include a storage duration term (rather than use the storage duration–corrected data) in GEE models where the time interval was restricted to a) 1- and 2-year-olds only or b) 5- and 7-year-olds only. The results are presented in Table 4. For 1- and 2-year-olds, age was significantly associated with all three substrate-specific measures of PON1 even after adjusting for potential cofounding by storage duration (p < 0.001 for AREase and CPOase and p = 0.016 for POase). Further, storage duration was only significantly associated with POase activity (p = 0.04), and the magnitude of the effect (−55.8 U/L) was much smaller than the effect of age (128.0 U/L). Similarly, for 5- and 7-year-olds, age was strongly associated with AREase, CPOase, and POase (p < 0.001) activities, and storage duration was not. For POase, the prominent interaction between age and PON1192 genotype observed in the storage duration–corrected models was also significant in these models (p = 0.02).
Correlation analysis. Correlations of each of the PON1 substrate-specific enzyme activities in children at different time points were highly significant except at 1 year of age, when there was a weak correlation with 7-year-olds and no correlation with cord blood samples at birth. This may be attributable to the relatively small sample size at 1 year of age. Further, because the cord blood measurement accounts for both fetal and placental expression, the correlation between these two time points (birth and 1 year of age) may be weaker than with other age groups. For pairwise comparisons among all other age groups (birth and 2, 5, and 7 years of age), Pearson’s correlation coefficients ranged from 0.37 to 0.62 for AREase and from 0.41 to 0.73 for CPOase assays. POase measurements were even more highly correlated, with correlation coefficients ranging from 0.70 to 0.93 [see Supplemental Material, Table 2 (available online at doi:10.1289/ehp.0900870.S1)]. Within each time point, all three PON1 activities were also significantly correlated with each other (p < 0.001). The correlation coefficients tended to be highest between CPOase and AREase, ranging from 0.84 to 0.97. The pairwise correlation coefficients between POase and CPOase ranged from 0.44 to 0.78. Between POase and AREase, the correlation coefficients were lower, ranging from 0.23 to 0.65 [see Supplemental Material, Table 3 (available online at doi:10.1289/ehp.0900870.S1)].
In this study, we examined the ontogeny of PON1 phenotypes in a longitudinal birth cohort of 458 children who were followed from birth through 7 years of age. Our data demonstrated a dramatic range of PON1 variability in young children, particularly in newborns. We observed an age-dependent increase in measures of PON1 levels and enzymatic activity that continued through 7 years of age. Furthermore, for POase activity, the rates of increase were allele specific, such that children with PON1192R alleles experienced a much more pronounced increase in POase activity over time. Similarly, children with the PON1-108CC192RR diplotype had significantly higher mean PON1 activities and also experienced steeper increases of POase activity over time compared with children with the PON1-108TT192QQ diplotype.
Earlier studies (Chen et al. 2003; Holland et al. 2006 ) demonstrated that newborns have lower PON1 levels than adults. Cole et al. (2003) reported that PON1 levels (AREase) and activities (POase) increased noticeably in nine children during the first year of life and appeared to have plateaued by 2 years of age in some but not all of them. In contrast, our study, which followed a much larger number of children for a longer period of time, demonstrated that children’s PON1 levels and activities continue to increase past 2 years of age until at least 7 years of age. These data suggest that children’s PON1 activities remain lower than their expected adult levels longer than previously thought. Although this phenomenon has not yet been reported in other epidemiologic studies, similar observations in mice have also shown that PON1 levels are low at birth and increase significantly with age (Zemke et al. 2009). Because PON1 can detoxify activated OPs and also has antioxidant properties, our findings provide evidence that the time in which infants and young children remain more vulnerable than adults to OP exposures and oxidative stress may extend from birth through 7 years of age and possibly even later.
In addition to the age-dependent increases in PON1 measures, our data in this large Mexican-American cohort corroborate allele and genotype frequencies of the PON1 SNPs at positions –108 and 192 in other studies of Mexican and Mexican-American populations (Rainwater et al. 2008; Rojas-Garcia et al. 2005). Consistent with other studies, we found that the Q allele of the PON1192 coding polymorphism was associated with decreased POase activities and that the C allele of the PON1-108 promoter polymorphism was associated with increased activity in all three PON1 measures (Deakin and James 2004; Furlong et al. 2005). In transgenic mouse models, the PON1192Q alloform does not protect against CPO exposure as well as the PON1192R alloform (Cole et al. 2005). Furthermore, mice with lower PON1 levels were more sensitive to OP toxicity as measured by cholinesterase inhibition in the brain (Li et al. 2000). These studies demonstrated that because of their influence on PON1 levels and catalytic efficiency, the PON1-108 and PON1192 genotypes are significant determinants of sensitivity to OP exposure. Furlong et al. (2006) used these data to estimate that the range of sensitivity to CPO in CHAMACOS mothers and their newborns is very wide (131- to 164-fold). Although age and PON1 genotype appear to affect this variability significantly, it is likely that other factors such as environmental exposures and additional genetic polymorphisms may account for some of this variability.
To our knowledge, no other studies have examined the effects of common polymorphisms on the ontogeny of PON1 phenotypes in young children. In our study, children with the PON1192R and PON1-108C alleles experienced much steeper increases in POase activity with age in a dose-dependent manner (dose of allele), indicating significant effect modification by PON1192 and PON1-108 genotypes. A similar pattern was reported by Parker-Katiraee et al. (2008), who found that PON1 expression in murine liver tissue increased during embryonic development in an allele-specific manner. These data suggest that PON1 ontogeny is influenced by genetic polymorphisms and children with certain genotypes (PON1192Q and PON1-108T alleles) may be more susceptible to the adverse effects of OP exposure and oxidative stress throughout childhood. In addition to PON1 genotypes, PON1 developmental expression may also be affected by prenatal exposure to OPs via epigenetic mechanisms. It will be interesting to determine whether OP exposure may be associated with PON1 ontogeny in future studies.
In our analyses, we performed GEE models using raw data and also using data corrected for storage duration. Previous studies, including our own, have reported a decline in POase, AREase, and diazoxonase activities over time, particularly in samples stored for several years (Brackley et al. 1983; Huen et al. 2009; Stenzel et al. 2007). To adjust for storage duration, we created a correction factor by constructing linear models from a data set of 95 subjects whose PON1 activities were measured at two time points. Although it is possible that the relationship between storage duration and PON1 activity may not be linear, it was the most reasonable assumption to use with only two time points available. This assumption would tend to overestimate the effect of storage duration, thus biasing the effect of age on PON1 activity toward the null. We also assumed that variance in measurements at the two time points was due primarily to storage duration and not other unmeasured factors. Although the storage duration–corrected data introduces a degree of uncertainty into our models, we were still able to identify highly significant age-dependent increases in all three PON1 measures using this data. In addition to the models using storage duration data, we were able to apply GEE models to subsets of the data (smaller time intervals) in which the correlation between age and storage duration was weaker. In these models, where a term for storage duration was included, we confirmed our findings that age is positively associated with PON1 activity after adjustment for storage duration. Despite the potential limitations of our analysis, this is the largest longitudinal study examining PON1 activity during childhood development. Use of a cross-sectional study design in future studies may help to confirm our results while avoiding the potential bias introduced by storage duration effects on PON1 activity.
From our data, it is not yet apparent when children’s PON1 levels and activities reach mature levels. Furthermore, levels did not increase appreciably between 2 and 5 years of age, but then continued to increase between 5 and 7 years of age. It is possible that the approach to puberty may affect levels of this physiologically important enzyme. Compared with younger children, average PON1 levels and activities in 7-year-olds were much closer to those we previously reported in pregnant CHAMACOS mothers (Holland et al. 2006). Whether children with different genotypes reach mature levels at different ages remains unclear. Because this study was performed in a cohort of Mexican-American children with relatively high levels of OP exposure, generalization of the age-dependent increase of PON1 activity to other populations should be done with caution. Future studies of PON1 ontogeny should include children as they grow older (past 7 years of age) and also additional populations.
Our present study provides additional evidence that children with the PON1192QQ and PON1-108TT genotypes may represent vulnerable populations. This may be of particular concern in the Mexican-American population, including our CHAMACOS subjects, because both genotypes are relatively common (23.2% QQ, 18.4% TT, and 7.5% QQ and TT). The PON1192Q and PON1-108T alleles are also quite frequent in Caucasian populations (73% Q, 62% T) but much less so in African-American populations (27% Q, 15% T) (Chen et al. 2003). Thus, sensitivity to OP toxicity and oxidative stress may vary between populations.
In summary, we followed a cohort of Mexican-American children from a farm-worker community from birth through 7 years of age and found that PON1 levels and activities continue to increase past 2 years of age and possibly even past 7 years of age. Additionally, PON1 polymorphisms modify the effect of age on enzyme activity such that some children experience much steeper increases over time than others. These findings suggest that children’s PON1 activities continue to increase until much later than previously believed, and therefore, the window of increased susceptibility to OP exposures and oxidative stress may also last much longer in young children. Future policies regarding protection from pesticide exposures should take into account that children may be particularly vulnerable to OP exposures and those with certain PON1 genotypes may be at even greater risk because of lower PON1 levels and activities.
Berkowitz GS, Wetmur JG, Birman-Deych E, Obel J, Lapinski RH, Godbold JH, et al. 2004. In utero pesticide exposure, maternal paraoxonase activity, and head circumference. Environ Health Perspect 112:388–391.
Bhattacharyya T, Nicholls SJ, Topol EJ, Zhang R, Yang X, Schmitt D, et al. 2008. Relationship of paraoxonase 1 (PON1) gene polymorphisms and functional activity with systemic oxidative stress and cardiovascular risk. JAMA 299(11):1265–1276.
Brackley M, Carro-Ciampi G, Stewart DJ, Lowden JA, Ray AK, Kalow W. 1983. Stability of the paraoxonase phenotyping ratio in collections of human sera with differing storage times. Res Commun Chem Pathol Pharmacol 41(1):65–78.
Bradman A, Eskenazi B, Barr DB, Bravo R, Castorina R, Chevrier J, et al. 2005. Organophosphate urinary metabolite levels during pregnancy and after delivery in women living in an agricultural community. Environ Health Perspect 113:1802–1807.
Bradman A, Whitaker D, Quiros L, Castorina R, Henn BC, Nishioka M, et al. 2007. Pesticides and their metabolites in the homes and urine of farmworker children living in the Salinas Valley, CA. J Expo Sci Environ Epidemiol 17(4):331–349.
Cole TB, Walter BJ, Shih DM, Tward AD, Lusis AJ, Timchalk C, et al. 2005. Toxicity of chlorpyrifos and chlorpyrifos oxon in a transgenic mouse model of the human paraoxonase (PON1) Q192R polymorphism. Pharmacogenet Genomics 15(8):589–598.
Connelly PW, Maguire GF, Picardo CM, Teiber JF, Draganov D. 2008. Development of an immunoblot assay with infrared fluorescence to quantify paraoxonase 1 in serum and plasma. J Lipid Res 49(1):245–250.
Deakin SP, Leviev I, Brulhart-Meynet MC, James RW. 2003. Paraoxonase-1 promoter haplotypes and serum paraoxonase: a predominant role for polymorphic position –107, implicating the Sp1 transcription factor. Biochem J 372(pt 2):643–649.
Eskenazi B, Harley K, Bradman A, Weltzien E, Jewell NP, Barr DB, et al. 2004. Association of in utero organophosphate pesticide exposure and fetal growth and length of gestation in an agricultural population. Environ Health Perspect 112:1116–1124.
Ferre N, Camps J, Fernandez-Ballart J, Arija V, Murphy MM, Ceruelo S, et al. 2003. Regulation of serum paraoxonase activity by genetic, nutritional, and lifestyle factors in the general population. Clin Chem 49(9):1491–1497.
Furlong CE, Cole TB, Jarvik GP, Pettan-Brewer C, Geiss GK, Richter RJ, et al. 2005. Role of paraoxonase (PON1) status in pesticide sensitivity: genetic and temporal determinants. Neurotoxicology 26(4):651–659.
Furlong CE, Holland N, Richter RJ, Bradman A, Ho A, Eskenazi B. 2006. PON1 status of farmworker mothers and children as a predictor of organophosphate sensitivity. Pharmacogenet Genomics 16(3):183–190.
Garces C, Lopez-Simon L, Rubio R, Benavente M, Cano B, Ortega H, et al. 2008. High-density lipoprotein cholesterol and paraoxonase 1 (PON1) genetics and serum PON1 activity in prepubertal children in Spain. Clin Chem Lab Med 46(6):809–813.
Holland N, Furlong C, Bastaki M, Richter R, Bradman A, Huen K, et al. 2006. Paraoxonase polymorphisms, haplotypes, and enzyme activity in Latino mothers and newborns. Environ Health Perspect 114:985–991.
James RW, Kalix B, Bioletto S, Brulhart-Meynet MC. 2005. Paraoxonase-1 promoter polymorphism C–107T and serum apolipoprotein AI interact to modulate serum paraoxonase-1 status. Pharmacogenet Genomics 15(6):441–446.
Jarvik GP, Rozek LS, Brophy VH, Hatsukami TS, Richter RJ, Schellenberg GD, et al. 2000. Paraoxonase (PON1) phenotype is a better predictor of vascular disease than is PON1(192) or PON1(55) genotype. Arterioscler Thromb Vasc Biol 20(11):2441–2447.
Karikas GA, Kriebardis A, Samara I, Schulpis K, Papachristodoulou M, Fytou-Pallikari A. 2006. Serum homocysteine levels and paraoxonase 1 activity in pre-school aged children in Greece. Clin Chem Lab Med 44(5):623–627.
Lawlor DA, Gaunt TR, Hinks LJ, Davey Smith G, Timpson N, Day IN, et al. 2006. The association of the PON1 Q192R polymorphism with complications and outcomes of pregnancy: findings from the British Women’s Heart and Health cohort study. Paediatr Perinat Epidemiol 20(3):244–250.
Li WF, Costa LG, Richter RJ, Hagen T, Shih DM, Tward A, et al. 2000. Catalytic efficiency determines the in-vivo efficacy of PON1 for detoxifying organophosphorus compounds. Pharmacogenetics 10(9):767–779.
Moser VC, Chanda SM, Mortensen SR, Padilla S. 1998. Age- and gender-related differences in sensitivity to chlorpyrifos in the rat reflect developmental profiles of esterase activities. Toxicol Sci 46(2):211–222.
Parker-Katiraee L, Bousiaki E, Monk D, Moore GE, Nakabayashi K, Scherer SW. 2008. Dynamic variation in allele-specific gene expression of Paraoxonase-1 in murine and human tissues. Hum Mol Genet 17(21):3263–3270.
Rojas-Garcia AE, Solis-Heredia MJ, Pina-Guzman B, Vega L, Lopez-Carrillo L, Quintanilla-Vega B. 2005. Genetic polymorphisms and activity of PON1 in a Mexican population. Toxicol Appl Pharmacol 205(3):282–289.
Stenzel J, Worek F, Eyer P. 2007. Preparation and characterization of dialkylphosphoryl-obidoxime conjugates, potent anticholinesterase derivatives that are quickly hydrolyzed by human paraoxonase (PON1192Q). Biochem Pharmacol 74(9):1390–1400.
Timchalk C, Kousba AA, Poet TS. 2007. An age-dependent physiologically based pharmacokinetic/pharmacodynamic model for the organophosphorus insecticide chlorpyrifos in the preweanling rat. Toxicol Sci 98(2):348–365.
Young JG, Eskenazi B, Gladstone EA, Bradman A, Pedersen L, Johnson C, et al. 2005. Association between in utero organophosphate pesticide exposure and abnormal reflexes in neonates. Neurotoxicology 26(2):199–209.
Zemke AC, Snyder JC, Brockway BL, Drake JA, Reynolds SD, Kaminski N, et al. 2009. Molecular staging of epithelial maturation using secretory cell-specific genes as markers. Am J Respir Cell Mol Biol 40(3):340–348.