- Research article
- Open Access
- Open Peer Review
Prediction models for short children born small for gestational age (SGA) covering the total growth phase. Analyses based on data from KIGS (Pfizer International Growth Database)
BMC Medical Informatics and Decision Making volume 11, Article number: 38 (2011)
Mathematical models can be developed to predict growth in short children treated with growth hormone (GH). These models can serve to optimize and individualize treatment in terms of height outcomes and costs. The aims of this study were to compile existing prediction models for short children born SGA (SGA), to develop new models and to validate the algorithms.
Existing models to predict height velocity (HV) for the first two and the fourth prepubertal years and during total pubertal growth (TPG) on GH were applied to SGA children from the KIGS (Pfizer International Growth Database) - 1st year: N = 2340; 2nd year: N = 1358; 4th year: N = 182; TPG: N = 59. A new prediction model was developed for the 3rd prepubertal year based upon 317 children by means of the all-possible regression approach, using Mallow's C(p) criterion.
The comparison between the observed and predicted height velocity showed no significant difference when the existing prediction models were applied to new cohorts. A model for predicting HV during the 3rd year explained 33% of the variability with an error SD of 1.0 cm/year. The predictors were (in order of importance): HV previous year; chronological age; weight SDS; mid-parent height SDS and GH dose.
Models to predict growth to GH from prepubertal years to adult height are available for short children born SGA. The models utilize easily accessible predictors and are accurate. The overall explained variability in SGA is relatively low, due to the heterogeneity of the disorder. The models can be used to provide patients with a realistic expectation of treatment, and may help to identify compliance problems or other underlying causes of treatment failure.
Small for gestational age (SGA) is a working diagnostic term used to describe foetuses or newborn infants who have a lower weight and/or length than what is normal for their gestational age, in the absence of any other specific diagnosis or reason for their small stature . Although in the majority of children born SGA catch-up growth occurs by 2 years of age, in about 10% catch-up growth does not occur. Without treatment, these children remain short and constitute some 20-25% of adults whose final height is below -2 SD scores . Although, these individuals are not growth hormone (GH) deficient, recent long-term studies have shown that treatment with recombinant human GH is able to improve height. Today GH has become an approved treatment in promoting catch-up growth  in short children born SGA. Experience with SGA children treated with GH to adult height is still very limited and studies have shown that a considerable fraction of the children treated does not achieve an adult height within the normal range. These observations demand an improved and individualized therapy with GH in short children with SGA. Recently, growth prediction models for the optimization of GH in SGA have been published based on patients documented within KIGS [4, 5]. The aim of this study is to test the existing prediction models and to develop additional growth prediction model from the KIGS database, to allow an optimisation and individualization of GH treatment from pre-pubertal onset to the end of the growth phase in SGA.
Patients included in this analysis were receiving recombinant human GH (Genotropin®; Pfizer Inc.) during follow-up and were documented in the pharmacoepidemiological survey - KIGS (Pfizer International Growth Database) by July 18, 2010. The KIGS survey is conducted in accordance with the Declaration of Helsinki . Diagnosis was made according to the KIGS aetiology classification list . Additional inclusion criteria for the patients were a birth weight and/or birth length for gestational age below -2.0 SD scores. Pre-pubertal stage was defined as a mean testes volume ≤ 3 ml in boys or Tanner breast stage B1 in girls. All patients received six or seven injections of GH per week. In order to avoid the effects of initial catch-up growth on total pubertal growth (TPG) only patients who were treated with GH for a total of at least 5 years and who were treated at least two years before puberty onset were included in the analysis of TPG.
The height standards used for normal children were those of Tanner et al  or - as indicated - those of Prader et al.  and the weight standards were those of Freeman et al. . Birth weight for gestational age was transformed to an SD score based on the standards of Niklasson et al . The MPH SD score was calculated as: (father's height SD score + mother's height SD score) ÷ 1.61 . Bone ages, calculated according to the method of Greulich and Pyle , were taken as reported by the treating physician.
The existing prediction models for the first two prepubertal years for short children born SGA  were tested on children documented within KIGS by July 18, 2010. The characteristics of the patients used for the development of prediction models for the first two prepubertal years are listed in Table 1. The existing prediction models for total pubertal growth (TPG) developed for boys and girls with idiopathic GHD  were tested for applicability to SGA children.
Development of the prediction model
Growth responses (annualized height velocities) during GH therapy were correlated, by multiple regression analysis, with potentially relevant variables as published before . It was attempted to develop prediction models for the third and fourth pre-pubertal year on GH. The variables tested were (i) status at birth: sex, weight SD score, length SD score, ponderal index, mode of delivery and APGAR score; (ii) genetic background: height SD score of the mother, height SD score of the father, and midparental height (MPH) SD score; (iii) treatment modality: GH dose (per kg body weight and per kg ideal body weight [weight for height]), frequency of GH injections, and accumulated years of GH treatment; (iv) patient variables at the beginning of the treatment period: age, bone age, height SD score, weight SD score, height SD score minus MPH SD score, the peak serum GH concentration during stimulation testing and the height velocity during the previous year of GH treatment.
Variables are reported as mean and standard deviation. SD scores were calculated as follows: SD score = (patient value - mean value for age- and sex-matched normal subjects) ÷ SD of the value for age- and sex-matched normal subjects. For comparison of patient groups, Student's t-test was used if data had a Gaussian distribution, and Wilcoxon's rank-sum test was used otherwise. Significance was considered at the 1% level (P < 0.01), unless otherwise specified. SAS® version 8 for Sun Solarix (SAS Institute, Cary, North Carolina) was used for all statistical analyses. The prediction models were developed by means of a multiple linear regression analysis fitted by least squares and the REG procedure in the SAS computer program (version 8.02, SAS Institute, Inc., Cary, NC). A hierarchy of predictive factors was derived by the all-possible regression approach, using Mallow's C (p) criterion for ordering predictive factors, as described previously [15, 16]. Differences between observed and predicted height velocities were expressed in terms of Studentized Residuals. The residual is calculated as the observed height velocity minus the predicted height velocity for each observation, and the Studentized Residual is the residual divided by its standard error, which is equivalent to a SD score.
Response to GH: 1st prepubertal year
A total of 2340 short children (62% boys) born SGA were treated a full first year with GH. The characteristics at start of GH treatment and the responses to GH are listed in Table 2. The characteristics of the patients and the parameters relevant for the developed first-year prediction model  are listed in Table 1. The patients were started on GH at an age of 7.2 (+/- 3.2) years and a height of -3.4 (+/- 1.0) SDS. The height velocity was 8.5 (+/- 2.1) cm/year and the delta height was 0.75 (+/- 0.44) SDS. The predicted height velocity was 8.8 (+/- 1.1) cm/year, which resulted in an SR of -0.27 (+/- 1.3 (see Table 2).
Response to GH: 2nd prepubertal year
A total of 1358 short stature children (65% boys) born SGA were treated a full second year with GH. The characteristics at start of GH treatment and the responses to GH are listed in Table 2. The characteristics of the patients and the parameters relevant for the developed second-year prediction model  are listed in Table 1. The patients were started on GH at an age of 6.7 (+/- 2.8) years and a height of -3.5 (+/- 1.0) SDS. The height velocity during the second year was 7.2 (+/- 1.5) cm/year and the delta height was 0.45 (+/- 0.28) SDS. The predicted height velocity was 7.1 (+/- 0.7) cm/year, which resulted in an SR of 0.04 (+/- 1.2).
Response to GH: 3rd prepubertal year
A total of 317 short stature children (72% boys) born SGA were treated a full third year with GH. The characteristics at start of GH treatment and the responses to GH are listed in Table 2. The patients were started on GH at an age of 5.5 (+/- 2.1) years and a height of -3.7 (+/- 1.1) SDS. The height velocity during the third year was 6.6 (+/- 1.2) cm/year and the delta height was 0.32 (+/- 0.22) SDS. The group was used to develop a third-year prediction model for SGA (see below).
Response to GH: 4th prepubertal year
A total of 182 short stature children (62% boys) born SGA were treated a full fourth year with GH. The characteristics at start of GH treatment and the responses to GH are listed in Table 2. The patients were started on GH at an age of 3.7 (+/- 0.9) years and a height of -4.2 (+/- 1.3) SDS. The height velocity during the fourth year was 6.1 (+/- 1.2) cm/year and the delta height was 0.20 (+/- 0.25) SDS. The predicted height velocity based on the fourth year prediction model of GHD patients  was 6.2 (+/- 0.8) cm/year, which resulted in an SR of -0.10 (+/- 1.2).
Response to GH: Total pubertal growth (TPG)
In 35 male patients whose puberty started at an age of 12.9 (+/- 1.1) years and a height of 140.1 (+/- 6.8) cm, and who received a mean GH dose of 0.26 (+/-0.1) mg/kg/week TPG was 23.0 (+/-6.2) cm. In 24 female patients whose puberty started at an age of 11.8 (+/- 1.8) years and a height of 135.4 (+/- 7.2) cm, and who received a mean GH dose of 0.25 (+/-0.1) mg/kg/week TPG was 17.3 (+/-6.3) cm. Based on the prediction models developed for TPG of males and females with GHD  the predicted TPG for boys and girls were 26.3 (+/-5.1) and 18.0 (+/-5.7) cm. These predicted numbers for height gain during puberty are not different from those actually observed.
3rd year prepubertal prediction model for short children born SGA
The results of the development of the prediction model for the 3rd prepubertal year based on 317 children (see Table 2) are listed in Table 3. Table 3 also gives the rank order of importance of the variables as predictors, the overall correlation coefficients of the prediction models (R2) and the error SD of the prediction in centimetres. The equation describing the predicted height velocity (PHV) for the third year of GH therapy (from Table 3) is as follows: PHV (cm/year) = 6.2 + [-0.18 × age at start (years)] + [0.19 × weight SD score at start] + [1.21 × GH dose (mg/kg/week)] + [0.12 × MPH SD score] + [0.26 × HV 2nd year (cm/year)] ± 1.0. The equation explains 33% of the variability.
The difference between observed and predicted height velocity during the first to fourth prepubertal year and of observed and predicted total growth during puberty expressed in terms Studentized are illustrated in Figures 1 a-e.
Validation of 3rd year prepubertal prediction model for SGA children
A subgroup of 34 children was randomly assigned from the total cohort identified for the validation of the derived model. The characteristics of these patients were not different compared to model group. The correlation between the predicted and the observed height velocity in the validation group was: R = 0.53 (p = 0.001) and the mean of Studentized Residual was 0.0 (SD = 1.0) [Figure 2a, b].
Most short stature children born SGA catch up within the first two years of life; however, about 10-15% remain permanently short . Obviously, since the causes of intrauterine as well as permanent growth impairment are often unknown, the terminology used to define these children is descriptive rather than causal . GH treatment in short children born SGA has been investigated over many years. Both the US Food and Drug Administration (FDA) and the European Agency for the Evaluation of Medical Products (EMEA) have approved GH for the treatment of short stature in children born SGA and Silver-Russell syndrome, provided that other causes of short stature are excluded, such as growth-inhibiting medication, chronic diseases, endocrine disorders, emotional deprivation, or syndromes associated with short stature. GH treatment for short stature in SGA children is standard practice today [18, 19]. So far only a few studies on SGA have documented the growth response to GH until adult height [20–28]. The data suggest that a considerable fraction of the short stature children born SGA treated with GH does not reach a normal height. Thus, an optimising of the treatment approach with GH is mandatory.
The response to GH can be expressed in different modes. Traditionally the response has been expressed in terms of annual height velocity (cm/year) which can be compared to normal age-related references either directly or transformed into SD scores. During puberty total height gain (cm) from the onset of puberty to the end of growth can be described. In addition height velocity can be compared with empirical somatograms derived from patients during GH therapy versus age and sex. Such height velocity targets  have been published for the first two prepubertal years for SGA children . Moreover, the response can be predicted based on mathematical algorithms derived from the analyses of large cohorts of children with the same diagnosis treated with GH . The observed and the predicted growth can be compared and expressed as a Studentized Residual; - an" index of responsiveness".
We and others [32–34, 17] ) have developed algorithms (models) for the prediction of the response to GH and have shown that the most important determinants of first pre-pubertal year growth on GH in short children born SGA including Silver Russell Syndrome (SRS) are the dose of GH (the higher, the greater) and the age of the children (the younger, the greater). We were also able to show that the height velocity observed during the first year is the major determinant of the second pre-pubertal year growth response to GH in SGA . Unlike in GHD there are no SGA growth prediction models existing for the subsequent pre-pubertal treatment years. The present analysis confirms that the SGA prediction models developed for the first two prepubertal years are suited to accurately predict height velocity in other patients documented within KIGS, who by the nature of the database are very heterogeneous. During the third prepubertal year of GH treatment again the observed previous height velocity and age are the most important predictors of height velocity. During the fourth year of GH treatment in SGA growth is only marginally above the spontaneous normal growth rate. Thus, most of the catch-up growth occurs during the first three years on GH. Since the number of children during the fourth prepubertal treatment year within KIGS did not exceed two hundred it was not appropriate to calculate an independent prediction model for SGA. However, the fact that the fourth-year prediction model for children with GHD  could be applied to the SGA cohort indicates, that after most of catch-up growth has occurred, similar mechanisms appear to be involved in the response to GH of both disorders at that stage of treatment.
Likewise, the number of individuals with all documented information required for the development of a disease-specific mathematical algorithm to predict TPG was presently still too small within KIGS. However, the models developed for male and female adolescents with GHD  could be applied to adolescents with SGA and gave accurate results. The models for TPG show that the anthropometrical characteristics at puberty onset are more important for TPG than the dose of GH, which is of course the only factor which can be modified by the treating physician.
The presently available growth prediction models derived from KIGS data are suited to predict height development during the most relevant phases of growth in short children born SGA. This will serve to adapt treatment modalities in such a way that an optimal height outcome can be achieved. The prediction models can be used to provide patients with a realistic expectation of treatment, and may help to identify compliance problems or other underlying causes of treatment failure.
Bryan SM, Hindmarsh PC: Normal and abnormal fetal growth. Horm Res. 2006, 65 (Suppl 3): 19-27.
Albertsson-Wikland K, Karlberg J: Postnatal growth of children born small for gestational age. Acta Paediatr. 1997, 423 (Suppl): 193-195.
Hokken-Koelega AC, De Ridder MA, Lemmen RJ, Den Hartog H, De Muinck Keizer-Schrama SM, Drop SL: Children born small for gestational age: do they catch up?. Pediatr Res. 1995, 38: 267-271. 10.1203/00006450-199508000-00022.
Ranke MB, Lindberg A, Cowell CT, Wikland KA, Reiter EO, Wilton P, Price DA, KIGS International Board: Prediction of response to growth hormone treatment in short children born small for gestational age: analysis of data from KIGS (Pharmacia International Growth Database). J Clin Endocrinol Metab. 2003, 88: 125-131. 10.1210/jc.2002-020867.
Ranke MB, Lindberg A: Height at Start, First-Year Growth Response and Cause at Birth Are Major Determinants of Adult Height Outcomes of Short Children Born Small for Gestational Age and Silver-Russell Syndrome Treated with Growth Hormone: Analysis of Data from KIGS. Horm Res Paediatr. 2010, 74: 259-266. 10.1159/000289570.
Declaration of Helsinki: Recommendations guiding doctors in clinical research. Adopted by the World Medical Association in 1964. Wiss Med J. 1967, 66: 25-6.
Ranke MB: The Kabi Pharmacia International Growth Study: aetiology classification list with comments. Acta Paediatr Scand. 1991, 379 (Suppl): 87-92.
Tanner JM, Whitehouse RH, Takaishi M: Standards from birth to maturity for height, height velocity and weight velocity: British children, 1965. Part II. Arch Dis Child. 1966, 41: 613-635. 10.1136/adc.41.220.613.
Prader A, Largo RH, Molinari L, Issler C: Physical growth of Swiss children from birth to 20 years of age. First Zurich longitudinal study of growth and development. Helv Paediatr Acta. 1989, 52 (Suppl): 1-125.
Freeman JV, Cole TJ, Chinn S, Jones PR, White EM, Preece MA: Cross sectional stature and weight reference curves for UK (1990). Arch Dis Child. 1995, 73: 17-24. 10.1136/adc.73.1.17.
Niklasson A, Ericson A, Fryer JG, Karlberg J, Lawrence C, Karlberg P: An update of the Swedish reference standards for weight, length and head circumference at birth for given gestational age (1977-1981). Acta Paediatr Scand. 1991, 80: 756-762. 10.1111/j.1651-2227.1991.tb11945.x.
Lindberg A, Ranke MB: Data analyses within KIGS. Growth Hormone Therapy in pediatrics - 20 Years of KIGS. Edited by: Ranke MB, Price DA, Reiter EO. 2007, Basel, Karger, 23-28.
Greulich WW, Pyle: Radiographic atlas of skeletal development of the hand and wrist. 1959, Stanford: Stanford University Press, 2
Ranke MB, Lindberg A, Chatelain P, Wilton P, Cutfield W, Albertsson-Wikland K, Price DA: Derivation and validation of A mathematical model for predicting the response to exogenous recombinant human growth hormone (GH) in prepubertal children with idiopathic GH deficiency. J Clin Endocrinol Metab. 1999, 84: 1174-1183. 10.1210/jc.84.4.1174.
Cook RD, Weisberg S: Residuals and influence in regression. 1982, Chapman and Hall, New York
Weisberg S: Applied linear regression. 1985, Wiley and Sons, Chichester
Ranke MB, Lindberg A, KIGS International Board: Observed and predicted total pubertal growth (TPG) during treatment with growth hormone in adolescents with idiopathic growth hormone deficiency (GHD), Turner syndrome (TS), short stature born small for gestational age (SGA) and idiopathic short stature (ISS). KIGS analysis and review. Horm Res Ped. 2011, 75: 423-432. 10.1159/000324117.
Saenger P, Czernichow P, Hughes I, Reiter EO: Small for gestational age: short stature and beyond. Endocr Rev. 2007, 28: 219-251.
Clayton PE, Cianfarani S, Czernichow P, Johannsson G, Rapaport R, Rogol A: Management of the child born small for gestational age through to adulthood: a consensus statement of the International Societies of Pediatric Endocrinology and the Growth Hormone Research Society. J Clin Endocrinol Metab. 2007, 92: 804-810.
Albanese A, Stanhope R: GH treatment induces sustained catch-up growth in children with intrauterine growth retardation: 7-year results. Horm Res. 1997, 48: 173-177. 10.1159/000185509.
Ranke MB, Lindberg A: Growth hormone treatment of short children born small for gestational age or with Silver-Russell syndrome: results from KIGS (Kabi International Growth Study), including the first report on final height. Acta Paediatr. 1996, 417 (Suppl): 18-26.
Coutant R, Carel JC, Letrait M, Bouvattier C, Chatelain P, Coste J, Chaussain JL: Short stature associated with intrauterine growth retardation: final height of untreated and growth hormone-treated children. J Clin Endocrinol Metab. 1998, 83: 1070-1074. 10.1210/jc.83.4.1070.
Rosilio M, Carel JC, Ecosse E, Chaussain JL: Adult height of prepubertal short children born small for gestational age treated with GH. Eur J Endocrinol. 2005, 152: 835-843. 10.1530/eje.1.01916.
Zucchini S, Cacciari E, Balsamo A, Cicognani A, Tassinari D, Barbieri E, Gualandi S: Final height of short subjects of low birth weight with and without growth hormone treatment. Arch Dis Child. 2001, 84: 340-343. 10.1136/adc.84.4.340.
van Pareren Y, Mulder P, Houdijk M, Jansen M, Reeser M, Hokken-Koelega A: Adult height after long-term, continuous growth hormone (GH) treatment in short children born small for gestational age: results of a randomized, double-blind, dose-response GH trial. J Clin Endocrinol Metab. 2003, 88: 3584-3590. 10.1210/jc.2002-021172.
Carel JC, Chatelain P, Rochiccioli P, Chaussain JL: Improvement in adult height after growth hormone treatment in adolescents with short stature born small for gestational age: results of a randomized controlled study. J Clin Endocrinol Metab. 2003, 88: 1587-1593. 10.1210/jc.2002-021123.
de Zegher F, Hokken-Koelega A: Growth hormone therapy for children born small for gestational age: height gain is less dose dependent over the long term than over the short term. Pediatr. 2005, 115: e458-e462. 10.1542/peds.2004-1934.
Dahlgren J, Wikland KA: Final height in short children born small for gestational age treated with growth hormone. Pediatr Res. 2005, 57: 216-222. 10.1203/01.PDR.0000148716.71231.81.
Bakker B, Frane J, Anhalt H, Lippe B, Rosenfeld R: Height velocity targets from the national cooperative growth survey for first-year growth hormone responses in short children. J Clin Endocrinol Metab. 2008, 93: 352-357. 10.1210/jc.2007-1581.
Ranke MB, Lindberg A, KIGS International Board: Observed and predicted growth responses in prepubertal children in children with growth disorders. Guidance of growth hormone treatment by empirical variables. J Clin Endocrinol Metab. 2010, 95: 1229-1237. 10.1210/jc.2009-1471.
Ranke MB, Lindberg A: Predicting growth in response to growth hormone therapy. GH IGF Res. 2009, 19: 1-11. 10.1016/j.ghir.2008.08.001.
Boguszewski M, Albertsson-Wikland K, Aronsson S, Gustafsson J, Hagenas L, Westgren U, Westphal O, Lipsanen-Nyman M, Sipila I, Gellert P, Muller J, Madsen B: Growth hormone treatment of short children born small-for-gestational-age: the Nordic Multicentre Trial. Acta Paediatr. 1998, 87: 257-263. 10.1111/j.1651-2227.1998.tb01434.x.
Albertsson-Wikland K, Kristrom B, Rosberg S, Svensson B, Nierop AF: Validated multivariate models predicting the growth response to GH treatment in individual short children with a broad range in GH secretion capacities. Pediatr Res. 2000, 48: 475-484. 10.1203/00006450-200010000-00010.
Schönau E, Westermann F, Rauch F, Stabrey A, Wassmer G, Keller E, Brämswig J, Blum WF, German Lilly Growth Response Study Group: A new and accurate prediction model for growth response to growth hormone treatment in children with growth hormone deficiency. Eur J Endocrinol. 2001, 144: 13-20. 10.1530/eje.0.1440013.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1472-6947/11/38/prepub
The authors express their thanks to all those KIGS investigators involved in providing data. Editorial assistance was provided by Margaretha Lindell. KIGS is sponsored by Pfizer Inc.
MBR is a consultant to KIGS and receives honoraria from Pfizer, NovoNordisk, Eli Lilly and Ipsen. MBR was a paid consultant to Pfizer in connection with the development of this manuscript.
AL is a full-time employee of Pfizer Endocrine Care, KIGS/KIMS/ACROSTUDY, Stockholm, Sweden
MBR has made major contributions in designing the analyses and in interpretation of data, in drafting the manuscript and revising it critically for important intellectual content. AL participated in the design of the analyses and performed the statistical analyses. Both authors read and approved the final manuscript.
Michael B Ranke and Anders Lindberg contributed equally to this work.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Ranke, M.B., Lindberg, A. Prediction models for short children born small for gestational age (SGA) covering the total growth phase. Analyses based on data from KIGS (Pfizer International Growth Database). BMC Med Inform Decis Mak 11, 38 (2011) doi:10.1186/1472-6947-11-38
- Growth Hormone
- Growth Hormone Treatment
- Growth Hormone Therapy
- Height Velocity
- Short Child