|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Journal of Lipid Research, Vol. 46, 1450-1456, July 2005 Determinants of variation in serum paraoxonase enzyme activity in baboons
* Department of Genetics, Southwest Foundation for Biomedical Research, San Antonio, TX 78245 Published, JLR Papers in Press, April 16, 2005. DOI 10.1194/jlr.M400473-JLR200
1 To whom correspondence should be addressed. e-mail: david{at}darwin.sfbr.org
Paraoxonase (PON), an HDL-associated enzyme, is one of many circulating antioxidants thought to play a vital protective role. To better understand the determinants of quantitative variation in serum PON activity, we assayed PON in samples from 611 pedigreed baboons fed three diets. PON was measured enzymatically; the main determinant of variation was genetic and consisted of at least three components: two loci detected by linkage analyses and a residual polygenic component. Multipoint linkage analyses gave peak log of the odds (LOD) scores on the baboon homolog of human chromosome 7q21-22 (near PON1, the structural gene) of 9.1 on the low-cholesterol, high-fat diet and 4.1 on the high-cholesterol, high-fat diet (genome-wide P values were 1 x 108 and 0.0018, respectively). Surprisingly, a second locus on the baboon homolog of human chromosome 12q13 gave a LOD score of 2.9 on the high-cholesterol, high-fat diet (genome-wide P value was 0.032). We identified several significant covariates, including age, sex, diet, and apolipoprotein A-I concentrations. We estimate that 53% of total trait variation in baboons is explained by genes and 17% by covariates, thus accounting for 70% of total variation in baboon PON. Although the generation of free radicals is influenced primarily by environmental factors, our findings suggest strong genetic regulation of one component in the antioxidant defense system that plays a major role in susceptibility to atherosclerosis.
Supplementary key words PON1 genetic linkage analysis
Aerobic organisms, including humans, derive most of their metabolic energy from the reduction of oxygen and, consequently, are susceptible to the damaging effects of small amounts of O2, OH, H2O2, and their unstable intermediates, which can react with nearly all biological macromolecules. Protection against oxidative damage has developed as a complex and poorly understood assembly of biomolecules and enzyme activities. However, protection is incomplete and variable among individuals. Oxidative damage has been invoked as a major causal factor in a number of degenerative disorders, including atherosclerosis (1, 2), diabetes (3, 4), neurological damage (5), and carcinogenesis (6) as well as in the process of aging (7). Recently, our group has begun to characterize the components of the antioxidant system that play a role in the process of atherosclerosis. A large number of proteins that have antioxidant activity are known to be associated with lipoproteins and are hypothesized to play key roles in protection from lipid oxidation and its consequences (8). Among these, paraoxonase (PON; aryldialkylphosphatase, EC 3.1.8.1) is a well-known enzyme associated specifically with HDL particles. Although, the actual physiological role of PON is controversial, several studies have demonstrated that it is capable of preventing and reversing LDL oxidation in vitro. It is speculated that the presence of PON on HDLs accounts in large part for the cardioprotective benefit of having high levels of HDL (9). Variation in the gene encoding PON (PON1) has been well studied in several human populations. Most consistently, the Q192R polymorphism is associated with major differences in enzyme activity, with the 192R form having higher and the 192Q form having lower PON activity. Other mutations in human PON1, such as M55L and C108T, also are associated with variation in enzyme activity levels (10). Baboons are a well studied model of human atherosclerosis and, in particular, of the genetic and dietary effects on risk factors for cardiovascular disease (11). In the present study, we have investigated the determinants of variation in serum PON activity in pedigreed baboons fed diets differing in levels of fat and cholesterol. We determined that PON activity levels are strongly heritable, and we obtained evidence from a genome screen for at least two distinct genes that influence baboon PON activity.
Animals and diet protocol Data were analyzed for 611 pedigreed baboons (225 males and 386 females) that averaged 9.6 years of age (SD = 6.7) and 17.9 kg in weight (SD = 8.4). Baboons were subjected to an experimental dietary challenge protocol described previously (12). Briefly, baboons were bled while consuming a basal diet low in fat and cholesterol (basal diet sample). These same baboons were then fed each of two high-fat diets for 7 weeks before obtaining blood samples. The diets were low-cholesterol, high-fat (LCHF) and high-cholesterol, high-fat (HCHF); the two diets were separated by a 7 week washout period on the basal diet. Blood samples were taken from the femoral arteries of baboons that had been fasted overnight and immobilized with ketamine (10 mg/kg) for an average of 10 min before bleeding (range, 515 min). Serum was prepared by low-speed centrifugation and frozen at 80°C in small aliquots protected from oxidation and desiccation (13). Sample aliquots were stored frozen until assayed and thus were subjected to a single freeze-thaw cycle. Animals were maintained at the Southwest Foundation for Biomedical Research, a facility certified by the Association for Assessment and Accreditation of Laboratory Animal Care International, and experimental procedures were approved by the Institutional Animal Care and Use Committee.
PON activity assay A concern with measurements of activity is the stability of the enzyme. To assess stability, we compared enzyme activity levels in two aliquots from each of eight serum samples; one aliquot was freshly thawed, and the other was left at room temperature for 30 h before assay. We found an average 3% decrease in enzyme activity (from 18.5 to 18.0 µmol/min/l) with incubation at room temperature; this decrease was not significant by ANOVA.
Baboon genotypes Data from nuclear families were used within the laboratories to identify for retyping of genotypes inconsistent with Mendelian inheritance. Genotype data were then further checked, using Markov-Chain Monte Carlo (MCMC) and descent path analysis algorithms implemented in the program SimWalk2 (17), to identify and blank two classes of genotypes: those inconsistent with Mendelian inheritance when data from the full pedigrees were considered and those that appeared to be highly unlikely (e.g., apparent double recombinants occurring within a narrow interval).
Baboon linkage map To allow multipoint linkage analyses to detect and localize quantitative trait loci (QTLs) for PON activity, we used the genotype data and information in the baboon whole genome linkage map to produce multipoint IBD matrices using the stochastic MCMC approach implemented in the computer package Loki (20, 21). This MCMC method yields multipoint IBD estimates that more closely approximate those obtained by exact methods that cannot be applied conveniently to large extended pedigrees, such as those in this study. Using the multipoint IBD matrices, we could perform a log of the odds (LOD) score evaluation at 1 cM intervals along each chromosome.
Baboon pedigree
Statistical analyses We conducted all statistical genetic analyses using a maximum likelihood-based variance decomposition approach (22) implemented in the computer program SOLAR (23). Basic univariate quantitative genetic analyses were used to simultaneously assess the additive effects of genes and the mean effects of selected covariates on PON activity data obtained from animals on each of the three diets. We estimated heritability as the proportion of the residual phenotypic variance (i.e., after accounting for the mean effects of covariates) attributable to the additive effects of genes. Using multivariate extensions to this approach (2426), we also estimated six additional parameters: the additive genetic and random environmental correlations between PON activity levels assayed in samples from animals on the three diets. Significance in all models was assessed by likelihood ratio tests (27). We used the variance decomposition approach to test for evidence of QTLs that affect variation in PON activity (23). This method entails specification of the genetic covariance between arbitrary relatives as a function of the IBD relationships at a given locus. The covariance matrix for a pedigree is modeled as the sum of the additive genetic covariance attributable to the QTL, the additive genetic covariance attributable to the effects of other loci, and the variance attributable to unmeasured environmental factors. We tested the hypothesis of linkage by comparing the likelihood of a restricted model in which variance attributable to the QTL was constrained to zero (i.e., no linkage) with that of a model in which it was estimated. The LOD score of classical linkage analysis was obtained as the quotient of the difference between the two loge-likelihoods divided by loge 10 (28). We hypothesize that the genetic contribution to variation in many complex traits is oligogenic (i.e., resulting from the actions of two or more genes with individually measurable effects, expressed over a polygenic background). Consequently, we performed whole genome linkage screens in a sequential manner to facilitate the detection of multiple QTLs for PON activity. That is, after a first screen that detected a putative QTL, we sequentially performed a series of screens in which the previously detected QTLs were fixed in each subsequent model. This was done until a screen failed to detect a putative QTL (defined as LOD > 1.9). A more detailed description of the theory and method underlying oligogenic linkage analysis was presented previously (23). To control for the overall false-positive rate given the finite marker locus density in the baboon genome linkage map, we estimated the genome-wide P value by the method of Feingold, Brown, and Siegmund (29). In this population, a LOD of 2.696 corresponded to a genome-wide P value of 0.05, which is considered significant. Using the same approach, the threshold for "suggestive" evidence of linkage, as proposed by Lander and Kruglyak (30), was a LOD of 1.461. Kurtosis can inflate LOD scores, but in no case did kurtosis exceed ±0.31 in the starting data or ±0.20 for the residual data in the null models for linkage analysis. In addition, we estimated a robust LOD score that was based on 10,000 simulations of the same trait and pedigree data with an unlinked marker (31). The simulations, implemented in SOLAR, yielded a correction factor that could be used to adjust the LOD score. However, for the three PON traits, we estimated adjustment factors >1, indicating that no adjustment was necessary.
Measurement of PON activity We assayed PON activity in serum samples from 611 pedigreed baboons that had been subjected to a dietary challenge protocol; the average value was 26.1 µmol/min/l serum. Table 1 gives average PON activity levels for male and female baboons fed each of three diets (basal, LCHF, and HCHF) and illustrates the >20% higher levels observed in females compared with males. We investigated the question of whether diet influenced enzyme activity by use of paired t-tests. These tests indicated no significant effect of increasing levels of dietary fat (i.e., basal vs. LCHF diet values; P = 0.61) but a small (3%) increase when levels of dietary cholesterol were increased (i.e., LCHF vs. HCHF diet values; P = 0.007).
Quantitative genetic analyses of PON activity on three diets Genetic models, implemented in SOLAR, included covariates for age and age squared, sex, and age-by-sex, which together accounted for 915% of total trait variance (Table 2). PON activity was significantly heritable on each of the diets [heritability (h2) for basal = 0.515; h2 for LCHF = 0.705; h2 for HCHF = 0.569]. Thus, between 51% and 71% of the residual variance was attributable to the additive effects of genes.
Multivariate genetic analyses of dietary effects Multivariate models were implemented in SOLAR, and we first investigated whether the same gene(s) influenced PON activity on each of the diets. The genetic correlations were 0.981 ± 0.041 for basal/LCHF, 0.995 ± 0.019 for LCHF/HCHF, and 0.995 ± 0.021 for basal/HCHF. In no case was the genetic correlation significantly different from unity, suggesting that all, or nearly all, of the same genes influenced enzyme activity on each of the diets. We estimated genetic variance for serum PON activity in a trivariate model that considered data from the three diets simultaneously. The genetic variation in PON activity on the LCHF diet ( G2 = 9.61) was 27% greater than that estimated for the HCHF diet ( G2 = 7.56; P = 0.055) but 44% greater than that for the basal diet ( G2 = 6.66; P = 0.039); the differences in genetic variance between the basal and HCHF diets were not significant. These latter results are consistent with a genotype-by-diet interaction in which the magnitude of the genetic effect on serum PON activity is influenced by the level of dietary fat and cholesterol.
Linkage analyses
Effects of HDL metabolism We tested three indicators of HDL metabolism for effects on PON activity by including them as covariates in the models presented in Table 2. Table 4 gives the estimates of HDL effects on PON activity and the associated P values. By far, the strongest covariate was concentration of the dominant protein component of HDL, apolipoprotein A-I (apoA-I), which explained 3.26.4% of the variance. It should be noted that apoA-I concentration data were not available for >100 animals in the study. HDL-cholesterol concentrations explained less variance than did apoA-I (2.73.7%), but all effects were significant. However, HDL median diameter, a measure of HDL particle size distributions, explained very little of the variance and was significant for only one diet. These data suggest that HDL metabolism exerts a significant effect on PON variation and that apoA-I concentration was the stronger predictor of the effect.
Components of PON variation To approximate the relative contributions of the various determinants of PON that were detected in this study, we ran models in which we estimated the effects of the covariates (age, sex, age-by-sex, age squared), QTLs on the baboon homologs of human chromosomes 7 and 12, and a residual additive genetic component. ApoA-I levels were not included in these models because >100 animals did not have this measure; thus, the effect of apoA-I was taken from the data given in Table 4. Table 5 gives a breakdown of the proportion of total trait variation explained by the various components we studied on each of the diets and also calculates an average value across diets. Overall, we were able to explain 70% of PON activity variance.
PON activity levels were reliably measured in baboon serum samples using published approaches. The coefficient of variation for a control product run with each plate was 7%. Baboon PON activity averaged 26 µmol/min/l, lower than we have observed for human samples ( 130 µmol/min/l; n = 1,284; data not shown) using the same assay approach.
By far, the predominant effect on PON activity was genetic, accounting for The remaining genetic variance (19%) was not accounted for by any single locus and is designated a "residual additive" genetic effect in Table 5. Given the fact that PON is one of the many antioxidant enzymes that help maintain a delicate balance between free radicals and antioxidants, it is not surprising that other genes affecting redox pathways might also exert significant influence on PON expression.
Of the residual, nongenetic variance, we detected effects by several measured covariates that explained
In our study, females had Finally, increases in dietary cholesterol caused significant increases in PON activity in this set of animals, but we were unable to estimate the size of the effect of diet. Because increasing dietary levels of fat and cholesterol increase HDL levels in baboons (40), it is possible that the dietary category served as a surrogate indicator of HDL metabolism. Studies of diet effects in humans have identified a number of dietary components that can significantly affect PON activity, including ethanol, certain types of fat, and fruits and vegetables (4144). Perhaps accounting for some of the reported effects of diet is a xenobiotic responsive element-like sequence within the human PON1 promoter region that appears to mediate the effects of polyphenols on PON expression in vitro (45). In contrast to the present observations, rabbits fed an atherogenic diet showed reduced PON activity. However, the atherogenic diet also caused a 70% reduction in HDL-cholesterol levels, again suggesting that some diet effects are actually exerted through effects on HDL metabolism. A possible concern with this report of two QTLs for PON activity is that they were not consistently significant on all diets (i.e., even though the diet means and variances were quite similar, the magnitude of the statistical support for each QTL differed noticeably). However, we can assert that the most significant signals on chromosomes 7 and 12 were replicated in analyses of data from each of the other diets by signals that were significant or suggestive at the genome-wide level (30). Although potential sources of error that could affect our level of confidence in a QTL include seemingly random biological variation and assay error, we believe that the apparent differences in linkage results were attributable, at least in part, to the effects of diet on gene expression. Evidence for genotype-by-diet interaction is accepted if we can determine that different genes influence the trait or if the genetic variances differ according to diet (26, 40). All three possible pair-wise genetic correlations were indistinguishable from unity, suggesting that the same set of genes influences PON variation on each diet. However, our comparisons of the diet-specific variances provided evidence that the genetic contribution to the variance in PON activity on the LCHF diet was significantly greater than on the other two diets. Thus, the results of our analyses are consistent with a genotype-by-diet interaction and, therefore, a biological basis for some of the differences in genetic models observed in this study. These results highlight the value of an animal model for characterizing gene-environment interactions that are virtually impossible to detect in free-living populations (such as humans) and yet are likely to be quite important individually. There are, however, limitations in our ability to extrapolate results from the baboon model to the human situation. First is the fact the our observations of a QTL on the homolog of human chromosome 12 have not yet been replicated in another population or species (although we report replication in the same animals fed different diets). Such replication would increase confidence that the locus might also contribute to genetic variance in humans. Another limitation of this study is that we have assessed PON activity using only the single substrate, paraoxon. It is well documented that various isoforms of human PON can show dramatically different substrate-dependent specific activities (10, 16), so it is conceivable that our results might not be true for other substrates or that we have missed critical genetic effects by choosing to assay with a single substrate. We do not anticipate that baboons will have the same polymorphisms as exert important effects on enzyme activity in humans. Therefore, until we can identify the relevant functional polymorphisms in the baboon structural gene, we cannot determine whether our results regarding baboon PON will apply directly to the human enzyme. Given the many metabolic similarities between the two species, however, it is likely that this report of a number of factors that determine PON variation will be relevant to humans.
In summary, our study has demonstrated significant genetic regulation of an important HDL-associated antioxidant, PON. Effects of at least two genes plus several covariates (age, sex, and HDL metabolism) together account for
The authors are grateful to Perry H. Moore, Jr., for technical assistance. This work was supported by grants from the National Institutes of Health (HL-28972 and RR-013986).
Submitted on
December 1, 2004
This article has been cited by other articles:
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||