| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Dutch Growth Foundation (M.A.J.d.R., A.C.S.H.-K.), 3001 KB Rotterdam, The Netherlands; Department of Epidemiology and Biostatistics (M.A.J.d.R.) Erasmus Medical Center-University Medical Center, 3000 CA Rotterdam, The Netherlands; Department of Pediatrics, Division of Endocrinology (A.C.S.H.-K.), Sophia Childrens Hospital, Erasmus Medical Center-University Medical Center Rotterdam, 3000 CB Rotterdam, The Netherlands; and Department of Medical Statistics (T.S.), Leiden University Medical Center, 2300 RC Leiden, The Netherlands
Address all correspondence and requests for reprints to: M. de Ridder, Dutch Growth Foundation, P.O. Box 23068, 3001 KB Rotterdam, The Netherlands. E-mail: m.deridder{at}erasmusmc.nl.
| Abstract |
|---|
|
|
|---|
Objective: Our objective was to develop a model for prediction of height at the onset of puberty and of adult height (AH).
Design and Setting: Two GH studies were performed in short SGA children.
Patients/Intervention: A total of 150 SGA children with height SD scores (SDS) less than –2, age 3 yr or older, no signs of catch-up growth, available height at the onset of puberty, and at least 1 yr of GH treatment before the onset of puberty were studied. In one study, patients were randomly assigned to either 0.033 or 0.067 mg/kg·d; in the other study all received 0.033 mg/kg·d. In 71 children, AH was reached.
Main Outcome Measures: Height SDS at the onset of puberty and AH SDS were calculated.
Results: Determinants positively related to height SDS at the onset of puberty were: height SDS at the start; target height SDS; and GH dose, whereas age at the start and female gender were negatively related. Positively related to AH SDS were: height SDS and chronological age – bone age at the start; target height SDS; and GH dose, whereas serum IGF binding protein (IGFBP)-3 SDS at the start was negatively related. There was a significant interaction between GH dose and IGFBP-3 SDS, indicating a smaller GH dose effect for higher levels of IGFBP-3. The final model explained 57% of the variance in height SDS at the onset of puberty and 41% of AH SDS.
Conclusions: The prediction model for height SDS at the onset of puberty and AH SDS of short SGA children treated with GH provides useful information about the expected long-term growth. Because GH dosage is one of the determinants, the model aids in determining the optimal GH dose for each child.
| Introduction |
|---|
|
|
|---|
In Europe, the recommended dose for the approved indication (The European Agency for the Evaluation of Medicinal Products, 2003) is 0.035 mg/kg·d. In the United States, the Food and Drug Administration-approved indication describes a dose of 0.070 mg/kg·d. Several studies found a dose effect on growth response during the first treatment years (7, 8, 10, 11, 12, 13, 14, 15). For the long-term growth response, the dose effect was smaller but still significant (4, 5). The optimal dose for individual short SGA children is not yet established. Some investigators stated that a dose of 0.033 mg/kg·d results in significant gains in long-term growth, with IGF-I levels in the normal range and at lower costs (16). Others argued that the lower dose might be sufficient for children without extremely short stature (above –3 SDS) but that shorter or older children might better start with a higher dose (
0.050 mg/kg·d), with tapering of the dose per kg when the absolute GH dose (in milligrams) is maintained over the years (5).
In the present study, we developed a model to predict height at the onset of puberty and adult height (AH) for short children born SGA who will start GH treatment. Determinants were various baseline characteristics and GH dose. The predictions from this model can be used as information about the expected AH for an individual child to decide on the prescribed GH dose.
| Subjects and Methods |
|---|
|
|
|---|
We used data from two GH trials in short SGA children (8, 17). These studies included children with birth length SDS for gestational age below –2.00 (18), height SDS for CA at the start below –2.00 (19), and height velocity SDS for CA below zero (19, 20), to exclude children with spontaneous catch-up growth. In both studies, CA at the start had to be 3 yr or older. Study 1 (8) included prepubertal children, defined as Tanner breast stage I for girls and a testicular volume less than 4 ml for boys (21), and age younger than 12 yr in boys and younger than 10 yr in girls. Study 2 (17) comprised a prepubertal group, aged 3–8 yr, and a group aged 8–14 yr. Both studies excluded children who had had a complicated neonatal period, endocrine or metabolic disorders, chromosomal disorders, growth failure caused by other disorders (emotional deprivation, severe chronic illness, or chondrodysplasia), or syndromes (except for Silver-Russell syndrome), as well as children who used or had used drugs that could interfere with GH treatment.
In study 1, children were randomly assigned to either 0.033 mg GH/kg·d (= 1 mg/m2·d, n = 41) or 0.067 mg GH/kg·d (2 mg/m2·d, n = 38). The prepubertal children in study 2 were randomized into either a GH group or a control group that remained untreated for 3 yr and started treatment afterward. In the older age group, all children were treated from the start. The GH dose in study 2 was 0.033 mg/kg·d. The inclusion period for study 1 was between April 1991 and January 1993, and for study 2 between October 1996 and December 1998.
For the present analysis, we selected children who were prepubertal for at least 1 yr after the start of GH treatment and had started or completed puberty at the time of analysis (December 2006). Birth year for boys had to be before 1993 and for girls before 1994. AH was used of boys born before 1988 and girls born before 1989.
Assays
Serum IGF-I and IGF binding protein (IGFBP)-3 levels were measured in one central laboratory, using a specific RIA (22). The intraassay coefficient of variation was 4%, and the interassay coefficient of variation was 6%. IGF-I and IGFBP-3 values were converted into SDS (23).
Measurements
Standing height was measured every 3 months using a Harpenden stadiometer. The mean of four measurements was used for analysis. Heights were converted into height SDS (19). AH SDS was computed using the reference values of adults (age > 20 yr). BA at the start was determined according to the Tanner and Whitehouse (24) radius, ulna, short bones score. At each visit, pubertal stages were assessed according to the method of Tanner and Whitehouse (21). The onset of puberty was defined as a breast stage II for girls and a testicular volume more than or equal to 4 ml for boys.
AH was defined as the height reached when height velocity had decreased less than 0.5 cm during the previous 6 months, and a BA of 15 yr or older for girls and 16.5 yr or older for boys.
Imputation of missing values and truncation of extreme values
Missing values were imputed using multiple imputation (25, 26). We generated five imputed data sets using the procedure SAS Proc MI (27). Variables with a nonnormal distribution were transformed to normality during the imputation procedure.
To restrict the influence of outlying values, the lowest values of the outcome and of the determinants were truncated to the first percentile and the highest values to the 99th percentile.
Development of the prediction model
The potential determinants were: 1) initial characteristics: sex, birth length SDS, and birth weight SDS (18), gestational age, and TH SDS (19); and 2) characteristics at the start of GH treatment: CA, height SDS, weight SDS, body mass index SDS (19), BA and CA – BA, maximum GH response to GH stimulation tests, serum IGF-I SDS, and IGFBP-3 SDS and GH dose.
We first developed separate models for height SDS at the onset of puberty and for AH SDS, using forward selection with an inclusion criterion of P < 0.05.
Next, we constructed a model for both outcomes (height SDS at the onset of puberty and AH SDS), using repeated measurements analysis. We started with all determinants selected in the two separate models as covariables, and the interaction terms for each determinant with time (0 = onset of puberty, 1 = AH), to allow for different effects on the two outcomes. Nonsignificant (P > 0.05) terms were excluded stepwise. After this selection, possible interactions between GH dose and each determinant in the model were tested.
Analysis of residuals
Using the model obtained by the previously described procedure, we calculated predicted outcomes and residuals (observed outcome minus predicted outcome). A possible relation between these values was assessed by examining the scatter plots, and by fitting the linear regression with the absolute values of the residuals as dependent and the predictions as independent variable.
Internal validation
For internal validation of the derived model, we used bootstrap techniques (28, 29). This is an important step in the procedure of development of a prediction model, needed to make the model less dependent from the data set. The predictive performance of a model applied to other data sets than the set on which it is derived, will be lower. Bootstrapping evaluates this difference by taking many random samples from the original data set. The models derived on these samples are consequently applied to the sample and to the original data set, and the performances are compared. The mean of the differences, called the optimism, is used to correct the residual SD and the R2 of the original model. Another result from bootstrapping is the shrinkage factor for the estimated coefficients of the model, which can be used to obtain the final prediction formula, corrected for overoptimism. For the internal validation of our model, we used 200 bootstrap samples.
| Results |
|---|
|
|
|---|
|
So, for the repeated measurements analysis, fitting one model for both outcomes, we started with candidate predictors TH SDS, GH dose, gender, age, CA – BA, height SDS, IGF-I SDS, and IGFBP-3 SDS (all measured at the start of treatment), and the interactions of these variables with time. In the backward selection, IGFBP-3 SDS was significant, but IGF-I SDS was removed because it had no significant contribution to the model. IGFBP-3 SDS was correlated with IGF-I SDS (r = 0.47; P < 0.0001) but proved to be a stronger determinant. The interactions of height SDS at the start, TH SDS, and GH dose effects with time (onset of puberty or AH) were not significant, indicating that the effects of these determinants are equal for both outcomes. The effects of age at the start and gender were only significant for height SDS at the onset of puberty, whereas IGFBP-3 SDS was only significant for AH SDS. For CA – BA, the effect on height at the onset of puberty was much smaller than the effect on AH SDS. There was a significant interaction between GH dose and IGFBP-3 SDS. This indicated that there was not one constant dose effect, but the dose effect was depending on the value of IGFBP-3 SDS at the start. A higher IGFBP-3 was related to a smaller effect of GH dose. The relation between the GH dose effect and the level of IGFBP-3 is illustrated in Fig. 1
. We plotted the outcome (height SDS at the onset of puberty in Fig. 1A
and AH SDS in Fig. 1B
), adjusted for gender, TH SDS, height SDS, age, and CA – BA at the start, against the IGFBP-3 SDS value. Through these scatter plots, separate splines were drawn for cases treated with dose 1 and cases treated with dose 2. The distance between the splines, which is decreasing with increasing value of IGFBP-3 SDS, represents the dose effect on the outcome. The dose effect is plotted at the bottom of each figure. The plots also show that there is only a minority (15%) of children with IGFBP-3 below –2.5 SDS, in which the GH dose effect is substantial.
|
|
In Table 3
and Fig. 2
, examples of predictions for four children are shown, presuming that either 0.033 mg GH/kg·d (dose 1) or 0.067 mg GH/kg·d (dose 2) is given. The predicted height SDS at the onset of puberty and the predicted AH SDS were plotted, both with 95% prediction interval. A reference line was drawn at –2 SDS, and the TH range (defined as the TH SDS ± 1.3) of the child was marked. For each child, AH SDS is lower than height SDS at the onset of puberty. This decrease in height SDS during puberty is mainly related to age at the start (for late starters the decrease is lower) and, to a lesser extent, to gender (more decrease for males) and IGFBP-3 SDS (more decrease for higher levels). At the right of each plot, the total gain in height SDS during treatment is plotted. By applying the prediction formula, it is expected that child 1 will achieve an AH above –2 SDS, with any of the GH dosages. For both dosages, predicted AH lies within the TH range, and the total gain in height SDS during treatment will be more than 1.5. Treatment with a GH dose of 0.033 mg/kg·d is likely to be sufficient for this child. Child 2 has a predicted AH below –2 SDS, even when dose 2 is given. However, with any dose, the predicted AH lies within the TH range. Dose 1 will result into a gain in height SDS of 1.34, whereas for dose 2, the predicted gain is 1.49 SDS. The small difference (1 cm) in the predicted growth response between the two doses is related to the relatively high serum level of IGFBP-3 (0 SDS). Therefore, for child 2, it does not seem reasonable to prescribe the higher dose. On the contrary, for child 3, the dose effect is expected to be substantial. With dose 1, an AH below –2 SDS is predicted, whereas the prediction is –1.23 SDS with dose 2, so within the normal range and within the TH range. This may justify prescribing a dose of 0.067 mg GH/kg·d. For child 4, aged 8 yr, the prospects are rather low, and the effect of a higher dose is limited, again because of the relatively high level of IGFBP-3 (0.5 SDS). The predicted gain in height SDS is less than one, with either dose. One may doubt whether GH treatment will be beneficial for this child.
|
|
| Discussion |
|---|
|
|
|---|
Our model includes generally available factors as TH SDS and height SDS and BA at the start of treatment. Serum IGFBP-3 levels are not always available for short children born SGA. When we removed IGFBP-3 SDS as a potential determinant, the model included IGF-I SDS as a significant determinant (estimated coefficient = –0.08; P = 0.03). This model had an 11% lower explained variance compared with our final model including IGFBP-3 SDS. There was no significant interaction between IGF-I SDS and GH dose.
We found a significant interaction between GH dose and IGFBP-3 SDS, with a negative coefficient. This means that the effect of a higher dose is depending on the level of IGFBP-3 of the child: a lower IGFBP-3 SDS at the start is related to a larger dose effect. For a child with an IGFBP-3 SDS of –2, the expected difference in AH is 0.65 SDS (
4.6 cm) when treated with 0.067 mg GH/kg·d compared with 0.033 mg GH/kg·d. For a child with an IGFBP-3 SDS of zero, this difference is only 0.11 SDS (0.8 cm). So, estimating the effect of a higher dose of GH treatment for an individual child requires a measurement of serum IGFBP-3 at the start of treatment. Therefore, we recommend the measurement of IGFBP-3 and calculation of its SDS in SGA children for whom GH treatment is considered.
When a model was constructed only for the outcome height at the onset of puberty (the first step in our procedure), IGF-I SDS was selected as a significant predictor. However, in the model only for the outcome AH, IGFBP-3 SDS was superior over IGF-I SDS as a predictor. There was a significant correlation between IGFBP-3 SDS and IGF-I SDS. Also in the repeated measurement model, in which we applied backward selection and started with both IGF-I SDS and IGFBP-3 SDS as candidate predictors, IGF-I SDS was not significant and, therefore, not included in the final model.
Internal validation of the prediction model was performed by bootstrapping. With this procedure, the model is corrected for overfitting and will give valid predictions with correct confidence limits for new patients. It is not necessary to perform an external validation. Applying the model to data of an independent cohort is still interesting.
Ranke et al. (7) developed prediction models for short-term growth response to GH treatment in short children born SGA. For both first and second-year responses (in cm/yr), the predicting factors were: age at the start (negatively related), weight SDS at the start, midparental height SDS, and GH dose (all positively related). In a Swedish study (9), a regression model for AH SDS was presented with paternal height SDS, height SDS at the start (both positively related), age at the start, and maximal GH response during provocation tests (both negatively related) as predictors. The Swedish study group differed from our group with respect to the percentage of children with GH deficiency (37% with a maximal GH response < 5.3 µg/liter vs. 21% in our study group) and the age range (2.5–15.1 vs. 3.0–11.2 yr in our study group). In addition, dosage of GH treatment was only randomized during puberty, and no significant influence of GH dose on the AH was found.
The interpretation of a predicted value for AH SDS depends on the objective of GH treatment. First, the goal might be to achieve an AH within the normal range (above –2 SDS). Second, the aim could be an AH in the TH range, usually defined as the TH ± 1.3 SDS. Third, a substantial gain in height SDS during treatment might be regarded as a reasonable goal. For an individual child, when the AH predicted for a GH dose of 0.033 mg/kg·d is unsatisfactory, a higher dose should be considered. Our study included only two fixed dosages (0.033 or 0.066 mg/kg·d). Therefore, for estimating results of intermediate doses, we have to specify the type of relation between GH dose and the effect on height SDS. Two assumptions are reasonable: a linear relation, meaning that each extra milligram of GH has the same additive effect; or a log linear relation, used in some studies (30, 31), meaning that the effect of an extra mg GH decreases with increasing dose. For example, the effect of 0.040 mg/kg·d compared with 0.035 mg/kg·d is larger than the effect of 0.065 mg/kg·d compared with 0.070 mg/kg·d. According to our model, for a child with IGFBP-3 SDS at the start of –2, the effect of 0.066 mg/kg·d compared with 0.033 mg/kg·d is 0.65 SDS. If the effect of GH dose is linear, the effect of a dose in the middle between these two doses, so 0.050 mg/kg·d, is half of 0.65, so 0.325 SDS. However, if the effect of GH dose is log linear, the effect of a dose of 0.05 mg/kg·d is 0.38 SDS. However, for doses between 0.033 and 0.066 mg/kg·d, the difference between the two assumptions (linear or log linear) is small, and, for simplicity, a linear relation can be assumed. Only for doses outside this range (extrapolation of dose effect), the difference might become large. Therefore, we recommend the following. For a child born SGA, the predicted AH SDS for a GH dose of 0.033 mg/kg·d should be computed. If this prediction is satisfactory, according to the considerations mentioned previously, the child is treated with this dose. If not, the predicted AH SDS for dose 0.066 mg/kg·d is computed. If this prediction is still relatively low, the benefit of any GH treatment should be doubted. If this prediction is in the desired range, the child is treated with this double dose. If the prediction is too far above the aimed AH SDS, the dose can be computed using the following formula: GH dose (mg/kg·d) = 0.033 *[aimed AH SDS – 0.11 – 0.66 *height SDS – 0.12 *TH SDS + 0.11 *IGFBP-3 SDS – 0.21 *(CA – BA)]/(0.15 – 0.27 *IGFBP-3 SDS) + 1.
For example, a child with a height SDS of –3.86, age 7.35 yr, BA 6 yr, IGFBP-3 SDS –2.46, and TH SDS –1.57 has a predicted AH SDS for dose 0.033 mg/kg·d of –2.07, so still outside the normal range. If a dose of 0.066 mg/kg·d will be used, the predicted AH SDS is –1.26. If the aimed AH SDS is –1.5, the required GH dose is computed as 0.056 mg/kg·d.
The prediction formulas given in Table 2
, as well as the formula for calculation of the suitable GH dose, described previously, will be implemented in the next update of the Growth Analyser (www.growthanalyser.org).
In conclusion, the presented model predicts the expected AH of a short child born SGA who will be treated with GH. This may help in determining the optimal GH dose for each child.
| Footnotes |
|---|
First Published Online November 13, 2007
Abbreviations: AH, Adult height; BA, bone age; CA, chronological age; IGFBP, IGF binding protein; SDS, SD score; SGA, small for gestational age; TH, target height.
Received June 21, 2007.
Accepted November 5, 2007.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
B. Kristrom, A. S. Aronson, J. Dahlgren, J. Gustafsson, M. Halldin, S. A. Ivarsson, N.-O. Nilsson, J. Svensson, T. Tuvemo, and K. Albertsson-Wikland Growth Hormone (GH) Dosing during Catch-Up Growth Guided by Individual Responsiveness Decreases Growth Response Variability in Prepubertal Children with GH Deficiency or Idiopathic Short Stature J. Clin. Endocrinol. Metab., February 1, 2009; 94(2): 483 - 490. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Endocrinology | Endocrine Reviews | J. Clin. End. & Metab. |
| Molecular Endocrinology | Recent Prog. Horm. Res. | All Endocrine Journals |