Relationship between serum insulin-like growth factor-1, IGF binding protein-3 levels and body height before and after gonadotropin-releasing hormone agonist therapy

Han Sol Song; Won Bok Choi; Joon Sup Song; Il Tae Hwang; Seung Yang

doi:10.6065/apem.2014.19.4.208

Ann Pediatr Endocrinol Metab > Volume 19(4); 2014 > Article

Song, Choi, Song, Hwang, and Yang: Relationship between serum insulin-like growth factor-1, IGF binding protein-3 levels and body height before and after gonadotropin-releasing hormone agonist therapy

Original Article

Annals of Pediatric Endocrinology & Metabolism 2014;19(4): 208-213.

Published online: December 31, 2014

DOI: https://doi.org/10.6065/apem.2014.19.4.208

Relationship between serum insulin-like growth factor-1, IGF binding protein-3 levels and body height before and after gonadotropin-releasing hormone agonist therapy

Han Sol Song, MD, Won Bok Choi, MD, Joon Sup Song, MD, Il Tae Hwang, MD, Seung Yang, MD

Department of Pediatrics, Hallym University College of Medicine, Seoul, Korea.

Address for correspondence: Seung Yang, MD. Department of Pediatrics, Hallym University Kangdong Sacred Heart Hospital, 150 Seongan-ro, Gangdong-gu, Seoul 134-701, Korea.
Tel: +82-2-2224-2251, Fax: +82-2-482-8334, drsyang@hallym.or.kr

Received September 12, 2014 Revised October 06, 2014 Accepted December 16, 2014

(open-access, http://creativecommons.org/licenses/by-nc/3.0/):

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Purpose

The gonadotropin-releasing hormone agonist (GnRHa) is widely used to treat patients with precocious puberty. However, its effect on growth is often difficult to predict because of the diverse nature of its causes and presentation. This study aims to show the impact of GnRHa treatment on insulin-like growth factor-1 (IGF-I) and IGF binding protein-3 (IGFBP-3) secretion, growth, and on other parameters that may help estimate the height velocity.

Methods

Data from 60 girls (mean age, 8.8±0.7 years) treated with GnRHa were analyzed. Their height, bone age (BA), serum IGF-I, and IGFBP-3 concentrations were measured at the start and after a year of GnRHa treatment. To eliminate the confounding effect of chronological age (CA), the standard deviation scores (SDSs) of their height, IGF-I, and IGFBP-3 concentrations according to their CA at the start and after a year of GnRHa treatment were calculated. We looked for possible correlations between these variables and compared the subgroups based on their height velocities and midparental heights.

Results

During their one-year GnRHa therapy, height SDS for CA significantly decreased to 0.81±0.83 (P<0.001), but height SDS for BA increased to -0.28±0.68 (P<0.001). There was no significant change in serum IGF-I SDS, IGFBP-3 SDS, and IGF-I/IGFBP-3 ratio. The advanced BA was the factor most strongly correlated to the height velocity (R=0.265, P=0.041).

Conclusion

These findings suggest that GnRHa treatment may affect the height velocity due to mechanisms other than suppression of the IGF-I and IGFBP-3 secretory axis.

Keywords: Precocious puberty, Insulin-like growth factor I, Body height

Introduction

Precocious puberty is generally defined as the appearance of secondary sex characteristics before the age of 8 years in girls and 9 years in boys¹⁾. Central precocious puberty (CPP) is a subtype of precocious puberty, which is caused by early maturation of the hypothalamic-pituitary-gonadal axis²⁾. Studies have reported a significant increase in the incidence of CPP³⁾. For example, in 2006-2011, there was a four- to five-fold increase in the number of CPP patients among Korean girls⁴⁾. As such, the physical and psychological effects of CPP on the affected children are becoming increasingly important.

Children with true precocious puberty or early puberty tend to exhibit temporary accelerated growth due to increased sex hormones, but also a shorter growing period, which ultimately lead to a lower final adult height (FAH). Increased sex hormones shortens growing period by promoting growth plate senescence, which refers to the structural and functional changes of epiphyseal plate including decline in chondrocyte proliferation and rate of longitudinal bone growth^5,6).

Treatment of precocious puberty with the gonadotropin-releasing hormone agonist (GnRHa) was introduced in 1981, and it is now widely used and is known to increase FAH^7,8). GnRHa binds to the gonadotropin-releasing hormone (GnRH) receptors in the gonadotropic cells of the anterior pituitary gland, initially resulting in a brief surge of luteinizing hormone (LH) and follicle-stimulating hormone (FSH). Then it is followed by desensitization of the GnRH receptors and eventual gonadal suppression due to down-regulation of the intrinsic pulsatile secretion of the LH and FSH⁹⁾. GnRHa treatment thus suppresses bone maturation, which normalizes the height velocity according to the chronological age (CA).

However some studies have reported a drop in the height velocity to below the normal range, possibly due to suppression of the secretory axis of the insulin-like growth factor I (IGF-I) and the IGF binding protein 3 (IGFBP-3) secondary to the reduced sex hormones¹⁰⁾.

The response of precocious puberty patients to GnRHa treatment is often difficult to predict, mainly because of the diverse nature of the clinical presentation of the disease and its underlying causes. Partsch et al.¹¹⁾ suggested that CPP was a continuum of clinical presentations and rates of progression that range from normal variants or incomplete forms of pubertal development to slowly progressive or transient forms and rapidly progressive forms. The factors affected the better response to GnRHa treatment and final height outcome are: a younger CA at puberty or at the start of treatment, an advanced bone age (BA), a shorter interval between the onset of the CPP and its treatment, a greater height at the initial evaluation, the target height, the predicted adult height before the treatment, and the duration of the treatment^12,13,14,15). Other genetic factors may also affect the efficacy of GnRHa treatment¹⁶⁾.

This study aims to show the impact of GnRHa treatment on IGF-I and IGFBP-3 secretion, growth, and on other anthropometric or laboratory parameters that may help estimate the height velocity.

Materials and methods

1. Patients

The subjects were 60 girls who started GnRHa treatment at our center from May 2010 to June 2012, including 10 girls who were diagnosed as having had true precocious puberty. All the subjects experienced breast development above Tanner Stage B2 before 9 years of CA, had recently experienced rapid growth in height or body weight, and more than one year advanced BA. The subjects who were born small for their gestational age, whose magnetic resonance imaging showed an intracranial lesion, who had other endocrine diseases, and who received concurrent growth hormone therapy were excluded. Long-acting GnRHa (triptorelin or leuprorelin 3.75 mg) was administered as subcutaneous injection every 28 days.

2. Measurements

The height, BA, and serum IGF-I, IGFBP-3, LH, FSH and estradiol concentrations of all the patients were measured at the start and after a year of GnRHa treatment. The height was measured to the nearest 0.1 cm using a Harpenden Stadiometer (Holtain Ltd., Crymych, UK). The heights of both parents of each subject were also recorded on a questionnaire. BA was assessed by comparing the plain left-hand radiographs of patients with reference data set published by Greulich and Pyle¹⁷⁾.

3. Statistics

To eliminate confounding effect of age, age-adjusted height velocity was defined as the difference between standard deviation score (SDS) of height after a year of treatment and that at the diagnosis. SDSs for IGF-I and IGFBP-3 concentrations based on the CA were also calculated. The SDSs were calculated using Cole's modified LMS methods¹⁸⁾. For the height, we used the L (box-cox power), M (median), and S (coefficient of variation) values provided by the Growth Charts of Korean Children and Adolescents, which was published by the Korea Centers for Disease Control and Prevention in 2007. For IGF-I and IGFBP-3, we estimated the parameters from the percentile values for the age of Korean girls using the Levenberg-Marquardt algorithm^19,20). The IGF-I:IGFBP-3 molar ratios were calculated at the start and after a year of treatment. The respective values at the start and after a year of treatment were compared using a paired t-test, and Pearson correlation coefficient was used to measure the correlation between the age-adjusted height velocity and other measured or derived parameters (i.e., the IGF-I SDS, IGFBP-3 SDS, and IGF-I/IGFBP-3 molar ratio at the start and after a year of treatment, the difference between the two values, and the midparental height). The subjects were divided into subgroups based on the differences in their height SDSs for CA or their mid-parental heights. The same statistical analyses were performed for each subgroup. The statistical analyses were performed using the IBM SPSS Statistics ver. 21.0 (IBM Co., Armonk, NY, USA). Value of P<0.05 was considered statistically significant.

Results

At the start of the GnRHa treatment, the mean chronological and BAs were 8.77±0.70 years and 10.38±0.76 years, respectively. Serum LH and estradiol concentrations decreased significantly, indicating adequate gonadal suppression (1.62±2.07 to 1.31±1.82, P<0.001; and 5.17±5.06 to 3.32±3.11, P=0.024, respectively). Serum FSH concentrations had no significant difference. The height SDS according to the CA was 1.01±0.84, and the height SDS according to the BA was -0.64±0.88. The SDSs for the serum IGF-I and IGFBP-3 concentrations were 0.48±0.89 and 2.96±0.94, respectively, and the IGF-I/IGFBP-3 ratio was 59.3±13.7. The mean midparental height was 159.6±3.7 cm. After a year of GnRHa treatment, the BA became 10.89±0.49 years. The difference between the BA and the CA significantly decreased from 1.61±0.82 to 1.12±0.72 (P<0.001), which indicates adequate suppression of the growth plate senescence. The height SDS according to the CA significantly decreased to 0.81±0.83 (P<0.001), but the height SDS according to the BA significantly increased to -0.28±0.68 (P<0.001). There was no significant change in the serum IGF-I, IGFBP-3 SDSs, and IGF-I/IGFBP-3 ratio (Table 1).

Girls whose BA was more advanced than their CA tended to have a greater age-adjusted height velocity for both their CA (R=0.265, P=0.041) and BA (R=0.642, P<0.001). Girls with a higher BA at the start of their treatment also had a greater age-adjusted height velocity for their BA (R=0.753, P<0.001). Age-adjusted height velocity for both CA and BA did not correlate with IGF-I, IGFBP-3, IGF-I:IGFBP-3 ratio, or sex hormones (Figs. 1 and 2). No other significant correlation between the variables was evident.

The mean age-adjusted height velocity for CA during the treatment was -0.20, and 33 subjects had a greater difference than the mean (Table 2). Between the two subgroups (that with a below-the-mean difference in their height SDS, and that with an above-the-mean difference in their height SDS), the above-the-mean group had a higher serum IGF-I concentration SDS for their CA after their treatment (0.52±0.96 vs. 0.25±0.60, P=0.048). The below-the-mean group had a higher IGF-I/IGFBP-3 ratio after their treatment (60.7±15.7 vs. 58.9±9.6, P=0.018). However, the ratio increased during the treatment in the above-the-mean group, whereas the below-the-mean group showed a decrease (2.0±9.0 vs. -1.6±15.9, P=0.037). No significant correlation between the variables was found in each subgroup.

Subgroup analysis was also performed between the below-the-mean and above-the-mean groups for their midparental heights (below and above: 159.6 cm; n=32 and 28, respectively) to elucidate the potential effect of genetic factors on the height velocity. The two subgroups showed no significant difference, and no significant correlation between the variables was evident in each subgroup.

Discussion

The role of IGF-I as a modulator of maturational timing is being increasingly recognized, and several investigators have examined the effect of gonadal suppression with GnRHa on the growth hormone axis and the height velocity with conflicting results. Significant decreases in the serum somatomedin-C concentration and in the mean growth hormone level after treatment with GnRHa have been reported^21,22). A previous analysis of the height velocities and the serum IGF-I and IGFBP-3 levels of 26 Korean girls also showed significant decreases in the serum IGF-I and IGFBP-3 concentrations after a year of GnRHa treatment, but unchanged normal height velocities²³⁾. Other studies have reported no significant change in the IGF-I and IGFBP-3 concentrations despite a decrease in the height velocity^24,25).

In this study, there was a significant decrease in the height velocity but not in the mean IGF-I and IGFBP-3 concentrations. These results imply that gonadal suppression may influence growth via direct pathways independent of the circulating IGF-I. The exact mechanism of the effect of sex steroids on the growth plate still needs investigation. Kanety et al.²⁴⁾ suggested that sex steroids have a role in inducing an increase in the IGF-I and IGFBP-3 axis, but may not be important in the subsequent maintenance of the levels of these hormones. This conclusion may be applicable to other groups of patients, since male patients with virilizing congenital adrenal hyperplasia have shown that hydrocortisone therapy suppresses testosterone without suppressing IGF-I^24,26).

Lack of change in IGF-I and IGFBP-3 with decrease in sex hormones level and height velocity signifies direct effect of sex hormones on growth. One possible explanation is that the sex hormones directly affect bone growth. In vitro and animal studies have shown that sex steroids may act via the locally produced IGF-I in the target tissues without significantly raising the circulating IGF-I concentrations, with the bone as a target organ⁹⁾. The exact quantitative role of growth and the sex hormones on the height velocity needs further investigation.

According to Weise et al.²⁷⁾, the height velocity is positively affected by the current estrogen levels, and negatively by growth plate senescence, which reflects prior estrogen exposure. The authors hypothesized that impaired linear growth due to excessive senescence is unmasked by treatment with GnRHa, which removes the stimulatory effect of estrogen on growth. In this study, the height SDS for the CA decreased after the GnRHa treatment, but the height SDS for the BA increased. The height velocity showed a stronger correlation with the BA than with the CA, which supports the hypothesis. A study of 38 Korean girls with CPP who had received GnRHa treatment also reported a significant decrease in the girls' height SDS for their CA and a significant increase in their height SDS for their BA²⁸⁾. Our findings were consistent with those of these studies. Furthermore, the increases in the height SDSs were greater in the girls whose BAs were more advanced than their CAs.

One of limitations of this study is that majority of the subjects were diagnosed as early puberty, and older than 8 years at the start of treatment. Although there was no evident correlation between CA at initial treatment and difference in height SDS in this study, previous reports have shown that GnRHa treatment is less effective in girls beyond 8 years of age^29,30). Younger patients with true precocious puberty may yield different results.

The dose of GnRHa administered during this study is relatively lower than that of previous studies with similar settings in general. Although sex steroids level does not correlate with dose of GnRHa as long as adequate gonadal suppression is achieved, the dose-effect of GnRHa on IGF-I/IGFBP-3 and growth needs further investigation³¹⁾. Another limitation is that all samples are collected before and during treatment, thus may not fully reflect long-term effects. A follow-up study is required to show longterm effect of GnRHa on IGF-I/IGFBP-3 axis, especially after cessation of treatment.

In summary, gonadal suppression with GnRHa results in a lower height velocity but no change in the serum IGF-I and IGFBP-3 concentrations. The decrease in the height velocity independent of the circulating IGF-I during GnRHa treatment suggests the presence of more direct mechanisms through which sex hormones affect bone growth. Additional studies are required to elucidate other factors that affect the height velocity during and after GnRHa treatment.

Notes

No potential conflict of interest relevant to this article was reported.

References

1. Grumbach MM, Styne DM. Puberty : ontogeny, neuroendocrinology, physiology, and disorders. Wilson JD, Foster DW, Kronenberg HM, Larsen PRet al., editors. Williams textbook of endocrinology. 9th ed. Philadelphia: W.B. Saunders. 1998;pp 1509–1625.

2. Chalumeau M, Chemaitilly W, Trivin C, Adan L, Breart G, Brauner R. Central precocious puberty in girls: an evidence-based diagnosis tree to predict central nervous system abnormalities. Pediatrics 2002;109:61–67. PMID: 11773542.

3. Mogensen SS, Aksglaede L, Mouritsen A, Sorensen K, Main KM, Gideon P, et al. Diagnostic work-up of 449 consecutive girls who were referred to be evaluated for precocious puberty. J Clin Endocrinol Metab 2011;96:1393–1401. PMID: 21346077.

4. Kim HK, Kee SJ, Seo JY, Yang EM, Chae HJ, Kim CJ. Gonadotropin-releasing hormone stimulation test for precocious puberty. Korean J Lab Med 2011;31:244–249. PMID: 22016677.

5. Weise M, De-Levi S, Barnes KM, Gafni RI, Abad V, Baron J. Effects of estrogen on growth plate senescence and epiphyseal fusion. Proc Natl Acad Sci U S A 2001;98:6871–6876. PMID: 11381135.

6. Byers S, Moore AJ, Byard RW, Fazzalari NL. Quantitative histomorphometric analysis of the human growth plate from birth to adolescence. Bone 2000;27:495–501. PMID: 11033444.

7. Carel JC, Lahlou N, Roger M, Chaussain JL. Precocious puberty and statural growth. Hum Reprod Update 2004;10:135–147. PMID: 15073143.

8. Crowley WF Jr, Comite F, Vale W, Rivier J, Loriaux DL, Cutler GB Jr. Therapeutic use of pituitary desensitization with a long-acting lhrh agonist: a potential new treatment for idiopathic precocious puberty. J Clin Endocrinol Metab 1981;52:370–372. PMID: 6780592.

9. Attie KM, Ramirez NR, Conte FA, Kaplan SL, Grumbach MM. The pubertal growth spurt in eight patients with true precocious puberty and growth hormone deficiency: evidence for a direct role of sex steroids. J Clin Endocrinol Metab 1990;71:975–983. PMID: 2401720.

10. Cara JF, Burstein S, Cuttler L, Moll GW Jr, Rosenfield RL. Growth hormone deficiency impedes the rise in plasma insulin-like growth factor I levels associated with precocious puberty. J Pediatr 1989;115:64–68. PMID: 2738797.

11. Partsch CJ, Heger S, Sippell WG. Management and outcome of central precocious puberty. Clin Endocrinol (Oxf) 2002;56:129–148. PMID: 11874402.

12. Adan L, Chemaitilly W, Trivin C, Brauner R. Factors predicting adult height in girls with idiopathic central precocious puberty: implications for treatment. Clin Endocrinol (Oxf) 2002;56:297–302. PMID: 11940040.

13. Paul D, Conte FA, Grumbach MM, Kaplan SL. Long-term effect of gonadotropin-releasing hormone agonist therapy on final and near-final height in 26 children with true precocious puberty treated at a median age of less than 5 years. J Clin Endocrinol Metab 1995;80:546–551. PMID: 7852518.

14. Allali S, Lemaire P, Couto-Silva AC, Prete G, Trivin C, Brauner R. Predicting the adult height of girls with central precocious puberty. Med Sci Monit 2011;17:PH41–PH48. PMID: 21629198.

15. Kauli R, Galatzer A, Kornreich L, Lazar L, Pertzelan A, Laron Z. Final height of girls with central precocious puberty, untreated versus treated with cyproterone acetate or GnRH analogue. A comparative study with re-evaluation of predictions by the Bayley-Pinneau method. Horm Res 1997;47:54–61. PMID: 9030968.

16. Oerter KE, Manasco P, Barnes KM, Jones J, Hill S, Cutler GB Jr. Adult height in precocious puberty after long-term treatment with deslorelin. J Clin Endocrinol Metab 1991;73:1235–1240. PMID: 1955504.

17. Greulich WW, Pyle SI. Radiographic atlas of skeletal development of hand wrist. Stanford, CA: Stanford University Press. 1971.

18. Cole TJ. The LMS method for constructing normalized growth standards. Eur J Clin Nutr 1990;44:45–60. PMID: 2354692.

19. Korea Centers for Disease Control and Prevention, Division of Chronic Disease Surveillance, Committee for the Development of Growth Standard for Korean Children and Adolescents. Korean Pediatric Society, Committee for School Health and Public Health Statistics. 2007 Korean children and adolescents growth standard (commentary for the development of 2007 growth chart). Cheongwon: Korea Centers for Disease Control and Prevention, Division of Chronic Disease Surveillance. 2007.

20. Hyun SE, Lee BC, Suh BK, Chung SC, Ko CW, Kim HS, et al. Reference values for serum levels of insulin-like growth factor-I and insulin-like growth factor binding protein-3 in Korean children and adolescents. Clin Biochem 2012;45:16–21. PMID: 22032863.

21. Harris DA, Van Vlie G, Egli CA, Grumbach MM, Kaplan SL, Styne DM, et al. Somatomedin-C in normal puberty and in true precocious puberty before and after treatment with a potent luteinizing hormone-releasing hormone agonist. J Clin Endocrinol Metab 1985;61:152–159. PMID: 3889035.

22. DiMartino-Nardi J, Wu R, Fishman K, Saenger P. The effect of long-acting analog of luteinizing hormone-releasing hormone on growth hormone secretory dynamics in children with precocious puberty. J Clin Endocrinol Metab 1991;73:902–906. PMID: 1909708.

23. Yi KH. Serum IGF-1 and IGFBP-3 levels in central precocious puberty girls treated with gonadotropin releasing hormone agonist (GnRHa). J Korean Soc Pediatr Endocrinol 2011;16:20–23.

24. Kanety H, Karasik A, Pariente C, Kauschansky A. Insulin-like growth factor-I and IGF binding protein-3 remain high after GnRH analogue therapy in girls with central precocious puberty. Clin Endocrinol (Oxf) 1996;45:7–12. PMID: 8796132.

25. Sklar CA, Rothenberg S, Blumberg D, Oberfield SE, Levine LS, David R. Suppression of the pituitary-gonadal axis in children with central precocious puberty: effects on growth, growth hormone, insulin-like growth factor-I, and prolactin secretion. J Clin Endocrinol Metab 1991;73:734–738. PMID: 1909701.

26. Belgorosky A, Rivarola M. Serum IGF-I during treatment of three different types of abnormal puberty, possible role of IGF-I in pubertal maturation [abstract 551] In: Proceedings of the 76th Annual Meeting of the Endocrine Society; 1994 June 15-18; Anaheim, CA, USA. Bethesda, MD, The Endocrine Society. 1994.

27. Weise M, Flor A, Barnes KM, Cutler GB Jr, Baron J. Determinants of growth during gonadotropin-releasing hormone analog therapy for precocious puberty. J Clin Endocrinol Metab 2004;89:103–107. PMID: 14715835.

28. Lee SJ, Yang EM, Seo JY, Kim CJ. Effects of gonadotropinreleasing hormone agonist therapy on body mass index and height in girls with central precocious puberty. Chonnam Med J 2012;48:27–31. PMID: 22570812.

29. Bouvattier C, Coste J, Rodrigue D, Teinturier C, Carel JC, Chaussain JL, et al. Lack of effect of GnRH agonists on final height in girls with advanced puberty: a randomized long-term pilot study. J Clin Endocrinol Metab 1999;84:3575–3578. PMID: 10522998.

30. Cassio A, Cacciari E, Balsamo A, Bal M, Tassinari D. Randomised trial of LHRH analogue treatment on final height in girls with onset of puberty aged 7.5-8.5 years. Arch Dis Child 1999;81:329–332. PMID: 10490438.

31. Badaru A, Wilson DM, Bachrach LK, Fechner P, Gandrud LM, Durham E, et al. Sequential comparisons of one-month and three-month depot leuprolide regimens in central precocious puberty. J Clin Endocrinol Metab 2006;91:1862–1867. PMID: 16449344.

Fig. 1

Scatter plot of the relationship between the changes in the IGF-I SDS and the height SDS during the GnRHa treatment. IGF-I, insulin-like growth factor I; SDS, standard deviation score; GnRHa, gonadotropin-releasing hormone agonist.

Fig. 2

Scatter plot of the relationship between the changes in the IGF-I:IGFBP-3 ratio and the height SDS during the gonadotropin-releasing hormone agonist treatment. IGF-I, insulin-like growth factor I; IGFBP-3, insulin-like growth factor binding protein 3; SDS, standard deviation score.

Table 1.

Characteristics of the 60 girls treated with the GnRH agonist

Characteristic	At the diagnosis	After a year
CA (yr)	8.77±0.70	-
BA (yr)	10.38±0.76	10.89±0.49
BA–CA (yr)	1.61±0.82	1.12±0.72^*
Height SDS for the CA	1.01±0.84	0.81±0.83^*
Height SDS for the BA	–0.64±0.88	0.28±0.68^*
IGF-I SDS	0.48±0.89	0.40±0.82
IGFBP-3 SDS	2.96±0.94	3.17±0.89
IGF-I/IGFBP-3	59.3±13.7	59.6±12.7
Midparental height (cm)	159.6±3.7	-

Values are presented as mean±standard deviation.

GnRH, gonadotropin-releasing hormone; CA, chronological age; BA, bone age; SDS, standard deviation score; IGF-I, insulinlike growth factor I; IGFBP-3, insulin-like growth factor binding protein 3.

^* P<0.05 compared to before the treatment.

Table 2.

Comparison of the two subgroups grouped according to the amount of change in their height SDS during their one-year treatment with GnRHa

Variable	Group A (n=33)^a)	Group B (n=27)^b)	P-value
Midparental height (cm)	160.2	158.8	0.281
CA before the treatment (yr)	8.64	8.91	0.553
Bone age (yr)
Before	10.31	10.46	0.926
After	10.84	10.94	0.563
Difference	0.53	0.48	0.937
Height SDS
Before	0.89	1.15	0.437
After	0.81	0.80	0.347
Difference	–0.08	–0.35	0.189
IGF-I SDS
Before	0.37	0.61	0.633
After^*	0.52	0.25	0.048
Difference	0.15	–0.35	0.275
IGFBP-3 SDS
Before	2.97	2.94	0.854
After	3.09	3.28	0.165
Difference	0.13	0.34	0.112
IGF-I/IGFBP-3
Before	57.3	62.3	0.189
After^*	58.9	60.7	0.018
Difference^*	2	–1.6	0.037

SDS, standard deviation score; GnRH, gonadotropin-releasing hormone; CA, chronological age; SDS, standard deviation score; IGF-I, insulin-like growth factor I; IGFBP-3, insulin-like growth factor binding protein 3; Δ, difference between the value after a year of treatment and that at the diagnosis.

^a) ΔHeight SDS≥–0.20.

^b) ΔHeight SDS<–0.20.

^* P<0.05 between the two subgroups.