SUPPLEMENT ARTICLE |

* Department of Psychology, University of California at Davis, Davis, California
Childrens Hospital Los Angeles, Los Angeles, California
| ABSTRACT |
|---|
|
|
|---|
Methods. The participants were the 413 children and their mothers from the International Maternal PKU Collaborative Study.
Results. Results supported a nonlinear relation between prenatal Phe exposure and offspring cognitive outcomes, with damage to the developing fetus if average Phe levels are above approximately 360 µmol/L. Moreover, prenatal Phe exposure had a strong effect on offspring outcomes at 1 year of age and was the only one of the background, pregnancy-related, or perinatal variables to influence directly offspring outcomes at 2, 4, and 7 years of age.
Conclusion. The present study was able to document the importance of prenatal exposure to Phe for predicting offspring cognitive outcomes in the presence of other predictors of these outcomes.
Key Words: prenatal effects phenylalanine PKU maternal PKU intelligence teratogenic effects socioeconomic status maternal intelligence maternal age pregnancy birth outcomes infancy childhood regression spline models structural equation models longitudinal study
Abbreviations: Phe, phenylalanine SES, socioeconomic status PKU, phenylketonuria SD, standard deviation PAH, phenylalanine hydroxylase MDI, Mental Development Index PDI, Psychomotor Development Index GCI, General Cognitive Index WISC-R, Wechsler Intelligence Scale for ChildrenRevised
Despite considerable research demonstrating detrimental effects of prenatal phenylalanine (Phe) exposure,1 much remains to be determined with regard to the effects of prenatal exposure to Phe on offspring developmental outcomes. One key question is the form of the relation between prenatal Phe exposure and offspring outcomes. If this relation is linear, then any prenatal exposure to Phe might cause damage to the developing fetus and thus be teratogenic, although Phe exposure might have to reach some specifiable level (eg, 600 µmol/L) before clinically significant problems can be documented. Conversely, if the relation is nonlinear, then different functional forms could be considered, and critical thresholds might be identified beyond which teratogenic effects occur.
One candidate model for the relation between prenatal Phe exposure and offspring outcomes is a 2-piece linear spline, which is characterized by 2 linear functions that intersect at a "knot" point. The underlying assumption of this model is that the relation between the predictor and outcome variables is a linear function with a given slope below the knot point and a linear function with a different slope above the knot point.
A 2-piece linear spline model with 0 slope below the knot point and negative slope above the knot point is consistent with a threshold of prenatal Phe exposure; although hypothesized by many researchers, there is considerable variability in estimates of threshold values required for teratogenic effects to occur. For example, Perry et al2 argued that maternal Phe levels >600 µmol/L were likely to be necessary before damage to the fetus occurs, but both Johnson3 and Hsia4 claimed that even higher levels of maternal Phe might be required to confirm a teratogenic effect of exposure. Consistent with this latter claim, Levy and Waisbren5 found in a sample of 53 offspring that mental retardation occurred, with 1 exception, only when the mothers Phe level exceeded 1080 µmol/L. However, the survey by Lenke and Levy1 documented mental retardation in 21% of offspring of mothers with Phe levels between 180 and 600 µmol/L.
A second important question involves the relative impact of prenatal exposure to Phe in the context of other variables. One of the most common correlates of offspring cognitive outcomes in infancy and childhood is parental intelligence,6 because there are both genetic and environmental contributions by natural parents to their offspring that are thought to influence intelligence. Recent reviews of the behavioral genetic literature7 conclude that genetic sources explain more than half of the variance in intelligence and that shared environmental factors explain small amounts of variance. In addition to parental intelligence, parental socioeconomic status (SES)8 and parental education9 are strong correlates of offspring intelligence throughout infancy, childhood, and adolescence. Advocates of genetic explanations of intelligence often argue that these presumed effects of SES and education of the parents on childhood outcomes are largely indirect effects of genetic sources of variance: parents with more favorable genetic endowments not only have higher intelligence but also have higher socioeconomic and education levels and provide more stimulating home environments.
Prenatal exposure to Phe is yet another factor that may influence offspring cognitive outcomes. Most research on maternal phenylketonuria (PKU) has focused on the relationship between maternal blood Phe levels during pregnancy and child outcomes, sometimes adjusting for other treatment and/or demographic factors, but has failed to estimate these effects within larger or more inclusive models for child developmental outcomes. Because mothers with higher prenatal Phe levels may have lower levels of intelligence, SES, and education, it is possible that maternal demographic factors may be more important than Phe levels during pregnancy when predicting offspring cognitive outcomes.
Thus, 2 general goals motivated the current study. The first goal was to determine the relative utility of linear and nonlinear regression models for representing the relation between prenatal Phe exposure and offspring cognitive outcomes. The second goal was to examine the relative effects of prenatal Phe exposure on offspring cognitive outcomes in the context of a range of maternal background variables, pregnancy-related variables, and indicators of perinatal status.
| METHODS |
|---|
|
|
|---|
Measures
Descriptive statistics for variables used in these analyses, presented in Table 1, include number of patients with valid scores, mean, standard deviation (SD), minimum, and maximum.
|
Maternal intelligence, assessed by the Verbal IQ and Performance IQ scores from the Wechsler Adult Intelligence ScaleRevised10 fell, on average, approximately 1 SD below the population mean for intelligence.
The final background variable, reflecting the severity of the genetic mutation on the phenylalanine hydroxylase (PAH) gene, was assessed by 2 variables: 1) assigned PHE level, the maximum blood Phe level on an unrestricted diet, and 2) a score developed by Güttler and Guldberg11 reflecting the severity of both mutations on the PAH gene, using DNA identification of the mutation. The Güttler score can range from 2 to 16, with higher numbers reflecting milder mutations.
Pregnancy-Related Variables
Five pregnancy-related variables were used in analyses reported in this article. Weight gain during pregnancy was assessed using 2 measures: 1) total weight gain in pounds during pregnancy and 2) percentage of recommended weight gain achieved at term. The second pregnancy-related variable was variability in Phe levels during pregnancy, assessed by the SD of all Phe levels during the pregnancy for an individual. The third was Phe level, assessed using 3 indicators: 1) the average of all Phe levels recorded during the pregnancy, 2) the week during pregnancy when all remaining Phe levels were <600 µmol/L, and 3) the week during pregnancy when all subsequent Phe levels were <360 µmol/L.
The fourth pregnancy-related variable was duration of the pregnancy, which ranged from approximately 28 to 43 weeks gestation. The final pregnancy-related variable was average daily protein intake during pregnancy.
Birth Variables
Three birth variables were included in the models: offspring length, weight, and head circumference.
Measures of Cognitive Status
The measures of offspring cognitive status used in the present study were obtained at ages 1, 2, 4, and 7 years. At ages 1 and 2 years, the Bayley Scales of Infant Development12 were administered, yielding Mental Development Index (MDI) and Psychomotor Development Index (PDI) scores at each age. When children were 4 years of age, the McCarthy Scales of Childrens Abilities13 were administered, which yielded General Cognitive Index (GCI) scores, as was the Test of Language Development,14 from which the Composite Spoken Language Quotient was used. Finally, at 7 years of age, the Wechsler Intelligence Scale for ChildrenRevised (WISC-R)15 was used to assess intelligence, from which the Verbal IQ and Performance IQ scores were used.
Analytic Procedures
To model the relation between prenatal Phe exposure and offspring intellectual outcomes, we used the REG and NLIN procedures within the SAS/STAT Software, Version 8 of the SAS system for PCs (SAS Institute, Inc, Cary, NC) to fit linear and 2-piece linear spline models to the data. Ordinary least-squares estimation was used under REG, and the multivariant secant method of estimation (option DUD) was used under NLIN assuming a slope of 0 for the segment representing lower maternal Phe levels. Because it is not possible to formulate statistical tests of the difference in fit of linear and 2-piece spline models, the models were evaluated with regard to their squared multiple correlation and the meaningfulness of their parameter estimates and associated standard errors.
When fitting structural equation models to the prenatal, perinatal, and offspring data, the LISREL 8.52 program16 was used. Because of the presence of missing data for some participants, we used full information maximum likelihood estimation of model parameters. Full information maximum likelihood estimation has been found to be an efficient and unbiased method of estimation in the presence of missing data that are missing at random and is one of the best alternatives to use when data do not meet the "missing at random" assumption. Model fit was evaluated using both the likelihood ratio
2 goodness of fit statistical test and practical fit indices, as described by Jöreskog and Sörbom.16
| RESULTS |
|---|
|
|
|---|
60 µmol/L). Relevant values are reported in Table 2 in mg/dL, in both systems of units.
|
The parameters of the 2-piece linear spline model also seem to have a more reasonable interpretation than those from the standard linear model. The standard linear model had an intercept of 119.5 and a regression slope of 2.47. These estimates suggest that predicted offspring MDI scores at the lowest observed maternal Phe levels (1.3 mg/dL or 78 µmol/L) would approach 120 and that MDI scores decrease linearly approximately 2.5 points for every 1 mg/dL (60 µmol/L) increase in the mothers average Phe level. Thus, offspring of mothers with average Phe levels of 10 mg/dL (600 µmol/L) had a mean predicted MDI score of 95.3, or (120 2.47 x 10). The predicted MDI of 116 at the lowest maternal Phe is extremely unlikely in this sample, given the mean IQ of the mothers, which fell approximately 1 SD below the population mean.
In the 2-piece linear spline model, the intercept of 108.1 provides constant estimate of average offspring MDI scores for maternal Phe levels up to the estimated knot point of 6.10 mg/dL (366 µmol/L). Above the knot point, predicted MDI scores decrease more rapidly than in the linear model, 3.08 points for every additional 1 mg/dL (60 µmol/L) in maternal Phe. Thus, offspring of mothers with average Phe levels of 10 mg/dL (600 µmol/L) had a mean predicted MDI score of 96.1, or (108 3.08 x 3.9). The 2-piece linear spline model supports the conclusion that average Phe exposure below 6.1 mg/dL (366 µmol/L) is not teratogenic and has no measurable effect but that substantial fetal damage occurs when average maternal Phe exposure levels exceed 6.1 mg/dL (366 µmol/L), with greater damage as Phe levels rise higher above the knot point.
The results for the PDI at 1 year of age, shown in Table 2, are similar to those for the MDI. Although the linear model explained considerable variance in PDI scores, the spline model explained approximately 2% more variance, all parameter estimates in both equations were statistically significant, and the parameter estimates from the spline model were more interpretable. The estimated knot point for PDI scores was 8.1 mg/dL (486 µmol/L), suggesting that higher levels of Phe exposure are required to affect PDI scores negatively than was the case for MDI scores.
Outcomes at 2 Years of Age
Table 2 report results for Bayley scores at 2 years of age. The linear model for the MDI explained slightly more variance than did the 2-piece linear spline model, a surprising outcome because the latter model has 1 more parameter than does the former model. Despite having slightly superior fit, the standard linear model had less interpretable parameter estimates, with predicted MDI scores approaching 120 at the lowest maternal Phe levels (ie, the intercept of 121.3). The 2-piece spline model had predicted MDI scores of approximately 110 for offspring whose mothers maintained average Phe levels below the knot point of 3.9 mg/dL (234 µmol/L) and a decrease of 3.3 MDI score units for every 1 mg/dL (60 µmol/L) increase in average Phe exposure beyond the knot point.
For Bayley PDI scores at 2 years of age, the 2-piece linear spline model explained more variance than did the standard linear model and also had more interpretable parameter estimates. The knot point of 5.6 mg/dL (336 µmol/L) for PDI scores was larger than that for MDI scores, suggesting that higher levels of Phe exposure are required to affect PDI scores negatively than was the case for MDI scores, but note that the knot points at 2 years of age (3.9 for MDI and 5.6 for PDI) were >2 mg/dL (120 µmol/L) lower than the estimates at 1 year of age (6.1 for MDI and 8.1 for PDI).
Outcomes at 4 Years of Age
Table 2 provides results from modeling McCarthy GCI scores. The 2-piece linear spline fit the data only slightly better than the standard linear model, but the parameter estimates for the spline model were more interpretable. In particular, predicted GCI scores were close to the population mean for offspring of mothers with Phe levels below the knot point of 3.9 mg/dL (234 µmol/L), after which GCI scores decreased approximately 3.1 points for every 1 mg/dL (60 µmol/L) increase in average Phe exposure. The knot point for GCI at age 4 was virtually identical to the knot point for MDI scores at 2 years of age.
Table 2 also shows results for Composite Spoken Language Quotient scores. The results closely approximated those for the GCI, with a better fit for the spline model and a knot point estimate of 4.5 mg/dL (270 µmol/L), after which teratogenic effects of Phe exposure occur.
Outcomes at 7 Years of Age
Table 2 show results for modeling WISC-R IQ scores at 7 years of age. Mirroring previous results, the 2-piece linear spline fit the data better than did the standard linear regression, explaining almost 3% more variance. Although all parameter estimates in both equations were statistically significant, the parameter estimates from the 2-piece linear spline were more easily interpreted. Specifically, predicted Verbal IQ scores were close to the population mean for offspring whose mothers had average Phe levels below the knot point of 6.8 mg/dL (408 µmol/L), above which Verbal IQ scores fell >4 points for every 1 mg/dL (60 µmol/L) increase in Phe exposure.
The results for Performance IQ scores closely parallel those for the Verbal IQ, with the 2-piece linear spline explaining >4% more variance than the standard linear regression model. The estimates from the 2-piece linear spline model predicted offspring Performance IQ near the population mean for maternal Phe below the knot point at approximately 6.8 mg/dL (408 µmol/L) and an estimated decrease in IQ scores of 4.25 points for every 1 mg/dL (60 µmol/L) increase after the knot point.
Modeling Influences of Prenatal and Perinatal Variables on Offspring Intellectual Outcomes
To model the relationship of maternal background, pregnancy-related, and birth variables to offspring cognitive outcomes, we formulated a large structural equation model that embodied 4 sets of variables. The first set of variables consisted of 5 background latent variables: 1) mothers SES, with Hollingshead scores as measured variable; 2) mothers education, with education ratings as its indicator; 3) mothers IQ, with the mothers Verbal IQ and Performance IQ scores from the Wechsler Adult Intelligence ScaleRevised as indicators; 4) mothers age, with chronological age at conception of child as indicator; and 5) severity of mutation, with assigned PHE level and Güttler score as the 2 indicators. In this analysis, the "direction" of variables for which higher values were less desirable (SES score and education code) was mathematically transformed so that higher values indicated better maternal status.
Five pregnancy-related latent variables were included: 1) weight gain, with total weight gain and percentage of recommended weight gain as indicators; 2) PHE variability, with SD of Phe levels during pregnancy as indicator; 3) PHE level, with average Phe level during pregnancy, weeks gestation after which all subsequent Phe levels were <600 µmol/L, and weeks gestation after which all subsequent Phe levels were <360 µmol/L as the 3 indicators; 4) weeks gestation, with number of weeks gestation at birth as indicator; and 5) protein intake, with average daily protein intake (in grams) as indicator.
Two birth-related latent variables were specified: 1) birth size, with birth weight and birth length as its 2 indicators, and 2) birth head circumference, with measured head circumference at birth as indicator.
Finally, 4 offspring cognitive outcome latent variables were specified: 1) Bayley MDI at 1 year; 2) Bayley MDI at 2 years; 3) McCarthy GCI at 4 years; and 4) WISC IQ, with Verbal IQ and Performance IQ as its 2 indicators. The first 3 of the offspring cognitive outcome latent variables were single-indicator latent variables.
Given limitations of space, the modeling of the relations among the preceding variables cannot be described in detail. Three latent variablesweight gain, protein intake, and birth sizewere excluded from the final model because the background latent variables had negligible effects on these variables and because the variables had negligible or 0 effects on offspring cognitive outcomes. In the final model, all single-indicator latent variables were assumed to be error-free, and indicators on multiple-indicator latent variables had moderate to high loadings (0.600.95, in standardized metric) on their respective latent variables.
The final structural equation model exhibited adequate levels of fit, with a root mean square error of approximation16 of 0.055. The final structural model, showing all directed paths in the model, is shown in Fig 1. All path coefficients shown in Fig 1 are statistically significant (P < .05), and all path coefficients with absolute magnitude >0.20 were significant past the P < .01 level.
|
The single birth latent variable remaining in the model was birth head circumference. Three variables had significant effects on head circumference, although the effect of mothers SES was of minor importance (ß = 0.11). The 2 strong effects were those of PHE level during pregnancy (ß = 0.46), with higher PHE levels leading to smaller head circumference, and weeks gestation (ß = 0.40), with longer pregnancies leading to newborns with larger head circumferences.
Turning to the 4 offspring cognitive outcome latent variables, we first specified autoregressive paths between these variables. That is, Bayley MDI at 1 year was allowed to influence Bayley MDI at 2 years, Bayley MDI at 2 years influenced McCarthy GCI at 4 years, and McCarthy GCI at 4 years influenced WISC-R IQ at 7 years. We then tested the effects of the maternal background, pregnancy-related, and birth-related latent variables on the cognitive outcome latent variables. Only PHE level during pregnancy had significant direct effects on the offspring cognitive outcome latent variables. That is, once PHE level during pregnancy was allowed to predict the cognitive outcome latent variables, none of the remaining potential predictorsnot even mothers IQ or birth head circumferencehad significant effects on the cognitive outcome latent variables. Moreover, the effects of PHE level on the cognitive outcomes were pervasive. As shown in Fig 1, PHE level had a strong effect on Bayley MDI at 1 year (ß = 0.57) and also had strong direct effects on Bayley MDI at 2 years (ß = 0.39), McCarthy GCI at 4 years (ß = 0.41), and WISC-R IQ at 7 years (ß = 0.39), despite the presence of autoregressive paths (ranging between 0.26 and 0.47) between the successive cognitive outcome measures.
| DISCUSSION |
|---|
|
|
|---|
The knot point in the 2-piece linear spline model is a key parameter, as exposure to levels falling below the knot point presumably have no teratogenic effects whereas exposure to levels above the knot point result in damage to the offspring. The estimates of the knot point varied somewhat across offspring age levels and cognitive outcomes, ranging between 234 and 486 µmol/L. Averaging across the several estimates of the knot point leads to an estimate of 330 to 360 µmol/L as the best estimate of the critical threshold of Phe exposure beyond which damage to the offspring occurs. A critical threshold in the range from 330 to 360 µmol/L is consistent with the data presented by Lenke and Levy,1 which showed elevated rates of mental retardation even in offspring with maternal Phe levels between 180 and 600 µmol/L, and is inconsistent with the much higher threshold values discussed by many researchers.2,5
The second goal was to estimate the effects of prenatal exposure to Phe within the context of a model that included an array of maternal background, pregnancy-related, and birth-related variables. The results of this modeling were concise. The severity of the mutation on the mothers PAH gene had stronger influences than the other background variables on pregnancy-related variables, including on both level and variability of maternal Phe levels during pregnancy. Thus, mothers with more severe mutations had higher Phe levels during pregnancy and exhibited greater variability in their Phe levels.
Of the pregnancy-related latent variables, only PHE level had notable effects on birth-related and offspring cognitive outcome variables, and all of these effects were large. PHE level had a large effect on birth head circumference and then had moderate to large direct effects on all 4 of the cognitive outcome latent variables. Importantly, PHE level mediated the effects of maternal background variables on offspring cognitive outcomes. That is, all 5 of the maternal background variables had significant influences on PHE level, and PHE level then significantly affected all remaining birth-related and cognitive outcomes. Moreover, with 1 minor exception (the effect of mothers SES on birth head circumference), the maternal background variables had no significant direct effects on the birth-related and cognitive outcomes once the effects of PHE level on these latent variables were estimated.
Taken together, these results point to the crucial nature of Phe level during pregnancy with regard to offspring cognitive outcomes. The model suggests that the influence of the PAH mutations on offspring cognitive outcomes could be prevented if Phe levels during pregnancy were controlled at low levels. If all mothers with PKU could keep their average Phe levels during pregnancy below the crucial threshold value of 330 to 360 µmol/L, then teratogenic effects of Phe exposure during pregnancy would likely be negligible or nonexistent, and the offspring of these mothers would have every expectation of having cognitive skills within the normal range.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
F. Maillot, M. Lilburn, J. Baudin, D. W Morley, and P. J Lee Factors influencing outcomes in the offspring of mothers with phenylketonuria during pregnancy: the importance of variation in maternal blood phenylalanine Am. J. Clinical Nutrition, September 1, 2008; 88(3): 700 - 705. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||