Introduction
============

Central precocious puberty (CPP) is identified by the attainment of breast and testicle development prior to the age of 8 years in girls and 9 years in boys.  Early development of secondary sexual characteristics are driven by the premature activation of the hypothalamic-pituitary-gonadal axis \[[@b1-apem-2017-22-4-259],[@b2-apem-2017-22-4-259]\].  Gonadotropin-releasing hormone analogs (GnRHa) treatment for CPP have been established as a gold-standard since the early 1980s \[[@b3-apem-2017-22-4-259]\].  GnRHa continuously stimulating the pituitary gonadotrophs leads it to desensitization and decreases the release of follicle stimulating hormone (FSH) and luteinizing hormone (LH) \[[@b4-apem-2017-22-4-259]\].  As such, GnRHa treatment postpones development of puberty and elongates the duration of growth and one maturation, thus, height potential is preserved \[[@b3-apem-2017-22-4-259]\]. 

Although the efficacy of GnRHa is well known, the long-term side effects remain controversial, especially in association with obesity and metabolic syndrome \[[@b5-apem-2017-22-4-259],[@b6-apem-2017-22-4-259]\].  Many previous studies have reported that GnRHa treatment is associated with changes in weight and body composition in precocious pubertal girls \[[@b7-apem-2017-22-4-259]\].  In obese or altered energy balance, metabolic modulators such as leptin and ghrelin regulate kisspeptin neuron \[[@b8-apem-2017-22-4-259]\].  Then kisspeptin neurons stimulate GnRH secretion and gonadal steroid production.  In the treatment of GnRHa, obesity may affect gonadotropin suppression, and this feedback may affect the clinical finding as BA advancement and height outcome.  However, few studies have reported the effect of obesity on the GnRHa treatment outcome. 

Therefore, we aimed to investigate the effect and the clinical course of GnRHa treatment in obese girls compared to girls of normal weight with CPP. 

Materials and methods
=====================

1.  Subjects
-----------

We reviewed the medical records of 182 girls with CPP who had been treated with GnRHa in the pediatric endocrinology clinic of the Korea University Medical Center from July 2008 to March 2016. 

The inclusion criteria for the diagnosis of CPP were (1) breast development before the chronological age (CA) of 8 year, (2) accelerated growth velocity, (3) advancement of bone age (BA) by 1 or more years above CA, and (4) peak LH concentration of ≥5 IU/L on the gonadotropin-releasing hormone stimulation test.  Girls with known endocrine disease or any other chronic illness that might affect pubertal development and body composition were excluded.  The patients were divided into 2 groups according to their BMI: overweight/obesity group, BMI ≥85th percentile for age and sex; normal weight group, BMI ≥3rd percentile and BMI \<85th percentile for age and sex using 2007 Korean national data \[[@b9-apem-2017-22-4-259]\]. 

2.  Methods
----------

All patients were administered with leuprolide acetate 75--100 μg/kg depot injection every 4 weeks and regularly checked for height, weight, and degree of pubertal development.  We evaluated the changes of height, weight, BMI, and sexual maturity in both groups at baseline, after 1 year, and at the end of GnRHa treatment.  Height and BMI were expressed as the standard deviation score (SDS) for CA and BA.  The SDS for height, weight, and BMI were calculated using the 2007 Korean National Growth Charts by LMS method \[[@b9-apem-2017-22-4-259]\].  BA was evaluated every 6 months by taking a radiograph of the left hand and determined using the Greulich-Pyle method \[[@b10-apem-2017-22-4-259]\].  The difference between BA and CA (BA--CA) was evaluated to assess degree of BA advancement.  Serum LH, FSH concentrations were measured by immunoradiometric assays.  Using a GnRH stimulation test, the concentrations of LH and FSH at baseline, 30, 45, 60, and 90 minutes after the synthetic GnRH injection were measured.  The predicted adult height (PAH) was calculated using the Bayley-Pinneau method \[[@b11-apem-2017-22-4-259]\].  We also compared the changes in BMI to both groups during GnRHa treatment. 

The sexual maturity ratings were determined using the Marshall and Tanner method \[[@b12-apem-2017-22-4-259]\].  The midparental height (MPH) was calculated by subtracting 6.5 cm from the parents' mean height \[[@b13-apem-2017-22-4-259]\]. 

Of the total subjects, 46 girls were followed up after discontinuation of treatment and evaluated on reaching their near adult height.  The near adult height was considered to be reached when BA was around 13.5 years. 

3.  Statistical analysis
-----------------------

Data analysis was performed using IBM SPSS Statistics ver.  20.0 (IBM Co., Armonk, NY, USA).  Comparison of the mean values of the 2 groups were performed by the Student t-test.  The paired *t*-test was applied to compare the values at baseline and after GnRHa treatment in each group.  Multiple linear regression analysis was performed to evaluate the several clinical factors influencing on changes in PAH.  Data are represented as mean±standard deviation.  Differences were considered to be statistically significant at *P*\<0.05. 

4.  Ethics Statement
-------------------

This study was approved by the Institutional Review Board of Korea University Anam Hospital (approval number: ED16208) and informed consent was waived by the IRB. 

Results
=======

1.  Comparison of clinical parameters of the normal weight and overweight/obesity groups at baseline
---------------------------------------------------------------------------------------------------

One hundred eighty-two female patients were included in the study: 108 patients in the normal weight group and 74 patients in the overweight/obesity group.  The auxological characteristics of both groups before and after GnRHa treatment are described in [Table 1](#t1-apem-2017-22-4-259){ref-type="table"}.  The mean CA and BA ate the start of GnRHa treatment were similar in both groups.  The MPH, PAH and height SDS were not significantly different between the 2 groups.  The mean BMIs of the normal weight group and the overweight/obese group were 16.87±1.28 kg/m^2^ and 20.81±1.36 kg/m^2^, respectively (*P*\<0.001).  The BMI SDS were 0.08±0.60 and 1.55±0.36 in each group (*P*\<0.001). 

2.  Comparison of clinical parameters of the normal weight and overweight/obesity groups at the end of the treatment
-------------------------------------------------------------------------------------------------------------------

All patients received GnRHa treatment for more than 1 year.  The 2 groups had a similar mean treatment duration (2.43±0.60 years vs. 2.45±0.59 years, *P*=0.787) ([Table 2](#t2-apem-2017-22-4-259){ref-type="table"}).  At the end of treatment, The CA was similar between the groups after the last GnRHa treatment, whereas the BA was more advanced and the height SDS for CA was higher in the overweight/obesity group.  A significant decrease in BA--CA during GnRHa treatment was similar between groups.  The height SDS for BA and PAH revealed similar results for both groups. 

3.  Changes in clinical parameters on growth outcome during GnRHa treatment
--------------------------------------------------------------------------

We followed up parameters indicating a treatment effect at 1 year and at the finish on GnRHa treatment ([Fig.  1](#f1-apem-2017-22-4-259){ref-type="fig"}).  In both groups, the BA--CA significantly and gradually decreased over the treatment period compared with baseline.  The height SDS for CA significantly decreased and the height SDS for BA increased during treatment in both groups, respectively.  At the end of treatment, PAH in both groups was significantly increased compared with baseline PAH.  The mean PAH increased by 1.61±2.87 cm in the normal weight group and by 1.86±3.01 cm in the overweight/obesity group, before and after treatment (all *P*\<0.001).  During the treatment of GnRHa, changes in parameters were compared between the 2 groups ([Table 3](#t3-apem-2017-22-4-259){ref-type="table"}).  The change in height SDS or BA and BA--CA after 1 year and end of the treatment between both groups are not different. 

4.  Correlation between the auxological factors and the △PAH at the end of treatment
-----------------------------------------------------------------------------------

Multiple regression analysis was performed to assess the factors influencing ΔPAH at the end of treatment ([Table 4](#t4-apem-2017-22-4-259){ref-type="table"}).  In both groups, ΔPAH was statistically correlated with baseline BA, baseline height SDS for BA and growth velocity during the 1st year.  Baseline BMI for CA was not significantly associated with ΔPAH at the end of treatment in both groups. 

5.  Changes in BMI during GnRHa treatment
----------------------------------------

Before and after GnRHa treatment, the mean BMI SDS in the normal weight group significantly increased ([Fig.  2](#f2-apem-2017-22-4-259){ref-type="fig"}).  The mean BMI SDS also increased significantly 0.27 at 1 year and 0.40 at the end of treatment, compared to baseline (all *P*\<0.001).  However, mean BMI-SDS in overweight/obesity group did not show significant change during this period.  At 1 year and the end of treatment, mean BMI SDS in the overweight/obesity group was 1.51 and 1.55 (*P*=0.201 and *P*=0.994, respectively). 

6.  Comparison of treatment efficacy according to change in BMI during GnRHa treatment
-------------------------------------------------------------------------------------

Of the total 182 patients, the BMI SDS was increased in 120 patients (65.9%) at the end of treatment ([Table 5](#t5-apem-2017-22-4-259){ref-type="table"}).  The baseline BMI SDS of the increased BMI SDS group was lower than that of the decreased BMI SDS group (*P*\<0.001).  There were no significant difference in ΔPAH, ΔBA--CA, Δheight SDS for CA, Δheight SDS for BA between the 2 groups. 

7.  Follow-up after discontinuation of GnRHa treatment
-----------------------------------------------------

After discontinuation of treatment, the near adult height of 46 girls, reaching approximately 13.5 years of BA, was subsequently evaluated ([Table 6](#t6-apem-2017-22-4-259){ref-type="table"}).  The near adult height of 29 girls in the normal weight group and of 17 girls in the overweight/obesity group were not significantly different (160.17±3.72 cm and 161.04±3.91 cm, respectively, *P*=0.569).  Their initial PAH was 157.91±3.08 cm and 158.31±2.43 cm, respectively.  The near adult height of 2 groups was significantly greater than their initial PAH (2.26±3.70 cm vs. 2.72±3.85 cm, respectively, *P*\<0.001).  Normal weight group revealed more increase in BMISDS compared to those of overweight/obesity group. 

Discussion
==========

In the current study we compared the clinical course and treatment effect between normal weight and overweight/obese girls during GnRHa treatment for CPP and GnRHa treatment improved the height outcome and delayed progression of bone maturation in both normal weight girls and obese girls with CPP.  GnRHa treatment efficacy was not affected by the obesity in female CPP patients.  The BMI SDS significantly increased only in normal weight female CPP patients during long-term GnRHa treatment. 

The number of CPP patients is increasing worldwide, and in Korea \[[@b14-apem-2017-22-4-259],[@b15-apem-2017-22-4-259]\].  As the report of the Health Insurance Review and Assessment Service, CPP patients increased annually about 44.4% from 2006 to 2010 in Korea \[[@b16-apem-2017-22-4-259]\].  In addition, a marked increase of obesity prevalence in childhood has been observed around the world \[[@b17-apem-2017-22-4-259],[@b18-apem-2017-22-4-259]\].  Many studies have been performed to investigate the association between precocious puberty and obesity \[[@b19-apem-2017-22-4-259]\].  A Swedish longitudinal study reported that an increase of BMI in childhood was associated with an earlier timing of pubertal onset, by as much as 0.7 years in girls and 0.6 years in boys \[[@b20-apem-2017-22-4-259]\].  Since recent reports have primarily focused on whether GnRHa treatment has increased BMI or not, we investigated the clinical parameters to evaluate whether the obesity affects the treatment outcome in CPP girls. 

In this study, the overall course of treatment in both groups was similar.  We found that BA-CA in both groups significantly decreased during follow-up.  It means that GnRHa also worked well in postponing bone maturity in the overweight/obesity group.  Both groups similarly revealed distinct improvement of PAH and height SDS for BA.  At the end of treatment, PAH increased 1.61±2.87 cm vs. 1.86±3.01 cm in the normal weight and overweight/obese groups, respectively, compared to the beginning of treatment.  Near adult height of 46 patients showed height gain of 2.26±3.70 cm in the normal weight group and 2.72±3.85 cm in the overweight/obese group, respectively, compared to their baseline PAH.  Near adult height was comparable to their MPH, although gain of height is smaller than that of previous other studies.  Thus, it is considered that GnRHa treatment can be effective in height outcome in our study.  According to previous study about GnRHa treatment effects in obese CPP girls, Yoon et al.  \[[@b21-apem-2017-22-4-259]\] reported that PAH SDS significantly increased in overweight/obesity group compared to in the normal weight group, and BA--CA similarly decreased in both groups during 1-year of GnRHa treatment in CPP girls. 

We performed a multiple regression analysis and compared the factors predicting ΔPAH in both groups, respectively.  Our results suggest that greater baseline BA and height SDS for BA and growth velocity for the first year of treatment indicate good growth outcomes, regardless of the treatment group.  However, body weight and BMI did not exhibit any statistical correlation with PAH at the end of treatment.  Obesity did not contribute significantly to the height outcome.  In previous other studies, Partsch et al.  \[[@b22-apem-2017-22-4-259]\] reported that initial BA advancement and treatment duration were factors that explained 68% of height gain.  Other studies found that MPH, CA, pretreatment height, height at the end of treatment, and growth rate during treatment were factors influencing near final height \[[@b23-apem-2017-22-4-259],[@b24-apem-2017-22-4-259]\]. 

Several previous studies have evaluated alterations in BMI before and after GnRHa treatment in CPP.  Whether GnRHa treatment affects BMI and aggravates obesity remain controversial \[[@b25-apem-2017-22-4-259]\].  Shiasi Arani and Heidari \[[@b26-apem-2017-22-4-259]\] reported that girls treated for CPP revealed no change in BMI compared with girls not treated for CPP, during 1 year of follow up.  In contrast, Paterson et al.  \[[@b27-apem-2017-22-4-259]\] reported a marked increase in mean BMI SDS after treatment from 0.93 to 1.2 in girls with early puberty and CPP.  On the other hand in our study, BMI SDS significantly increased in the normal weight group only during GnRHa treatment.  The reasons explaining an increase in BMI SDS in normal weight group are unclear.  Paterson et al.  \[[@b27-apem-2017-22-4-259]\] found that the mean BMI SDS returned to the pretreatment state 3 years after finishing GnRHa treatment course, suggesting that the weight increase was a transient phenomenon.  Wolters et al.  \[[@b28-apem-2017-22-4-259]\] explained the beneficial effect of GnRHa on weight status by reducing sex hormones that lead to weight gain through increase in body fat and muscle mass during GnRHa treatment.  However, this is insufficient to account for the difference in the impact of GnRHa on weight status between overweight and normal weight children.  In our study, 19 of the normal weight groups showed increased BMI and they were turned into overweight or obese at the end of treatment, whereas, 89 patients maintained normal weight persistently.  Patients who turned into overweight/obesity revealed higher BMI SDS at baseline and more increase in BMI SDS compared to patients with normal BMI SDS at the end of treatment (0.51 vs. -0.01, *P*\<0.001 and 0.74 vs. 0.23, *P*\<0.001, respectively).  On the other hand, overweight/obese patients at baseline showed a tendency to maintain their BMI SDS, and this finding is consistent with that in previous studies \[[@b28-apem-2017-22-4-259],[@b29-apem-2017-22-4-259]\].  Pich et al.  \[[@b30-apem-2017-22-4-259]\] reported that obese patients may have tendency to recognize their obese status and modify risk factors such as diet, sleep and physical activity.  However, our study cannot investigate these factors.  A prospective further study is warranted to confirm the tendency. 

The strength of our study is that all subjects were had CPP, except for those in early puberty and that they were separated into 2 groups; namely overweight and normal weight girls.  However, in our retrospective study, there were some limitations.  There were no clinical measures of body adiposity, such as waist circumference, skinfold thickness or dual-energy x-ray absorptiometry.  Although BMI is a valid tool to determine obesity, it has limitations as an indirect measure of fat mass.  A long-term follow-up following discontinuation of GnRHa treatment is lacking in our study.  To confirm the effects of obesity on the treatment outcome, additional analysis of the final height, BMI SDS and menarche age is required. 

In conclusion, obese girls with CPP had an improved height outcome similar to normal weight girls during GnRHa treatment.  Obesity may not affect the efficacy of GnRHa in girls with CPP. 

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

![ Changes in treatment-effect parameters during gonadotropin-releasing hormone agonist treatment.  BA--CA (A), height SDS for CA (B), and height SDS for BA (C).  BA, bone age; CA, chronological age; BA--CA, difference between BA and CA; SDS, standard deviation score.  \**P*\<0.001 compared to baseline of treatment.](apem-2017-22-4-259f1){#f1-apem-2017-22-4-259}

![ Changes in the BMI SDS during gonadotropin-releasing hormone agonist treatment.  BMI, body mass index; SDS, standard deviation score. \* * P*\<0.001 compared to baseline of treatment.](apem-2017-22-4-259f2){#f2-apem-2017-22-4-259}

###### 

Auxological characteristics of normal weight and overweight/obesity groups at baseline of gonadotropin-releasing hormone analogs treatment

  Characteristic                Normal (n=108)   Overweight/obesity (n=74)   *P*-value
  ----------------------------- ---------------- --------------------------- --------------------------------------------------------------
  CA (yr)                       8.55±0.47        8.47±0.49                   0.226
  BA (yr)                       10.21±0.66       10.28±0.71                  0.492
  Mid parental height (cm)      160.01±3.48      159.16±3.61                 0.114
  Predicted adult height (cm)   157.58±2.82      158.02±2.66                 0.297
  Tanner stage                                                               0.006^[\*](#tfn1-apem-2017-22-4-259){ref-type="table-fn"}^
   B2                           76 (70.4)        36 (48.6)                   
   B3                           29 (26.9)        35 (47.3)                   
   B4                           3 (2.8)          3 (4.1)                     
  Height (cm)                   135.27±5.54      136.03±5.08                 0.350
  Height-SDS for CA             1.14±0.84        1.35±0.70                   0.080
  Height-SDS for BA             -0.47±0.50       -0.42±0.48                  0.488
  Weight (kg)                   30.95±3.69       38.60±4.19                  \<0.001^[\*](#tfn1-apem-2017-22-4-259){ref-type="table-fn"}^
  Weight-SDS                    1.82±0.61        1.69±0.42                   \<0.001^[\*](#tfn1-apem-2017-22-4-259){ref-type="table-fn"}^
  BMI (kg/m^2^)                 16.87±1.28       20.81±1.36                  \<0.001^[\*](#tfn1-apem-2017-22-4-259){ref-type="table-fn"}^
  BMI-SDS                       0.08±0.60        1.55±0.36                   \<0.001^[\*](#tfn1-apem-2017-22-4-259){ref-type="table-fn"}^
  Basal LH                      0.72±0.73        0.74±0.63                   0.870
  Peak LH                       16.37±12.97      12.53±11.79                 0.043^[\*](#tfn1-apem-2017-22-4-259){ref-type="table-fn"}^
  Basal FSH                     3.14±2.44        2.99±2.34                   0.697
  Peak FSH                      16.20±8.38       15.65±7.89                  0.658

Values are presented as mean±standard deviation or number (%). 

CA, chronological age; BA, bone age; SDS, standard deviation score; BMI, body mass index; LH, luteinizing hormone; FSH, follicle stimulating hormone. 

*P*\<0.05, statistically significant differences. 

###### 

Auxological characteristics of normal weight and overweight/obesity groups at the end of gonadotropin-releasing hormone analogs treatment

  Characteristic            Normal (n=108)   Overweight/obesity (n=74)   *P*-value
  ------------------------- ---------------- --------------------------- --------------------------------------------------------------
  CA (yr)                   10.98±0.36       10.89±0.42                  0.134
  BA (yr)                   11.77±0.44       11.94±0.38                  0.005^[\*](#tfn2-apem-2017-22-4-259){ref-type="table-fn"}^
  BA--CA (yr)               -0.78±0.50       -1.05±0.55                  0.001^[\*](#tfn2-apem-2017-22-4-259){ref-type="table-fn"}^
  Treatment duration (yr)   2.43±0.60        2.45±0.59                   0.787
  PAH (cm)                  159.19±3.25      159.88±3.41                 0.169
  ΔPAH (cm)                 1.61±2.87        1.86±3.01                   0.577
  Height (cm)               148.31±4.90      150.34±4.40                 0.011^[\*](#tfn2-apem-2017-22-4-259){ref-type="table-fn"}^
  Height-SDS for CA         0.68±0.78        1.00±0.75                   0.007^[\*](#tfn2-apem-2017-22-4-259){ref-type="table-fn"}^
  Height-SDS for BA         -0.04±0.61       0.08±0.62                   0.204
  BMI (kg/m^2^)             19.25±1.88       23.22±1.71                  \<0.001^[\*](#tfn2-apem-2017-22-4-259){ref-type="table-fn"}^
  BMI-SDS                   0.40±0.64        1.55±0.42                   \<0.001^[\*](#tfn2-apem-2017-22-4-259){ref-type="table-fn"}^

Values are presented as mean±standard deviation. 

CA, chronological age; BA, bone age; PAH, predicted adult height; SDS, standard deviation score; BMI, body mass index. 

*P*\<0.05, statistically significant differences. 

###### 

Comparison of clinical parameters between normal weight and overweight/obesity groups during gonadotropinreleasing hormone analogs treatment

  Variable                  Normal (n=108)   Overweight/obesity (n=74)   *P*-value
  ------------------------- ---------------- --------------------------- -----------
  ΔHeight-SDS for CA, 1yr   -0.12±0.29       -0.06±0.15                  0.088
  ΔHeight-SDS for CA, end   -0.46±0.30       -0.35±0.26                  0.012\*
  ΔHeight-SDS for BA, 1yr   0.36±0.42        0.35±0.40                   0.815
  ΔHeight-SDS for BA, end   0.43±0.49        0.50±0.46                   0.356
  ΔBA--CA, 1yr              -0.46±0.41       -0.42±0.39                  0.482
  ΔBA--CA, end              -0.87±0.52       -0.77±0.45                  0.163
  ΔBMI-SDS, 1yr             0.19±0.38        -0.04±0.28                  \<0.001\*
  ΔBMI-SDS, end             0.32±0.46        0.00±0.37                   \<0.001\*

Values are presented as mean±standard deviation. 

SDS, standard deviation score; CA, chronological age; BA, bone age; BA--CA, difference between BA and CA; BMI, body mass index. 

###### 

A multiple regression analysis on ΔPAH according to obesity groups

  Variable                          Total    Normal   Overweight/obesity                                               
  --------------------------------- -------- -------- -------------------- -------- -------- ------- -------- -------- -------
  Baseline BA                       1.265    0.275    \<0.001              1.315    0.374    0.001   1.187    0.423    0.007
  Baseline height-SDS for BA        4.411    0.932    \<0.001              3.553    1.163    0.003   5.614    1.635    0.001
  Growth velocity during 1st year   68.982   15.449   \<0.001              62.201   18.233   0.001   82.757   31.659   0.011
  Baseline BMI-SDS                  0.141    0.212    0.507                0.342    0.402    0.397   -0.621   0.856    0.471

Total group: *R*^2^=0.317, *P*\<0.001; normal group: *R*^2^=0.316, *P*\<0.001; overweight/obesity group; *R*^2^=0.343, *P*\<0.001. 

PAH, predicted adult height; SE, standard error; BA, bone age; SDS, standard deviation score; BMI, body mass index. 

###### 

Comparison of 2 subgroups according to change in BMI during gonadotropin-releasing hormone analogs treatment

  Variable                     ΔBMI-SDS\>0 (n=120)   ΔBMI-SDS\<0 (n=62)   *P*-value
  ---------------------------- --------------------- -------------------- -----------
  Baseline BMI-SDS             0.50±0.86             1.03±0.85            \<0.001\*
  End BMI-SDS                  0.94±0.76             0.73±0.84            0.087
  ΔBMI-SDS                     0.45±0.32             -0.30±0.21           \<0.001\*
  Baseline PAH                 157.73±2.72           157.82±2.86          0.825
  ΔPAH                         1.80±2.78             1.54±3.20            0.572
  Baseline BA--CA              1.74±0.52             1.69±0.43            0.617
  ΔBA--CA                      -0.81±0.50            -0.87±0.49           0.507
  Baseline height-SDS for CA   1.27±0.77             1.16±0.83            0.377
  ΔHeight-SDS for CA           -0.41±0.25            -0.44±0.34           0.451
  Baseline height-SDS for BA   -0.43±0.49            -0.50±0.49           0.425
  ΔHeight-SDS for BA           0.47±0.45             0.45±0.52            0.811

Values are presented as mean±standard deviation. 

BMI, body mass index; SDS, standard deviation score; PAH, predicted adult height; BA--CA, difference between BA and CA; CA, chronological age; BA, bone age. 

###### 

Final data of 46 girls after discontinuation of gonadotropin-releasing hormone analogs treatment

  Variable                  Normal (n=29)   Overweight/obesity (n=17)   *P*-value
  ------------------------- --------------- --------------------------- -----------
  CA (yr)                   12.85±0.88      12.62±1.23                  0.480
  BA (yr)                   13.91±0.70      14.06±0.90                  0.544
  Midparental height (cm)   160.01±3.48     159.16±3.61                 0.114
  Baseline PAH (cm)         157.91±3.08     158.31±2.43                 0.645
  Near adult height (cm)    160.17±3.72     161.04±3.91                 0.456
  Baseline BMI-SDS          0.09±0.68       1.42±0.50                   \<0.001\*
  ΔBMI-SDS, 1yr             0.09±0.43       -0.12±0.30                  0.084
  ΔBMI-SDS, end             0.40±0.56       -0.06±0.36                  0.003\*

Values are presented as mean±standard deviation. 

CA, chronological age; BA, bone age; PAH, predicted adult height; BMI, body mass index; SDS, standard deviation score.