Longitudinal changes in growth hormone secretory process irregularity assessed transpubertally in healthy boys

Steven M. Pincus1, Johannes D. Veldhuis2,3, and Alan D. Rogol3,4,5

1 Guilford, Connecticut 06437; 2 Department of Internal Medicine, General Clinical Research Center, Center for Biomathematical Technology, 3 National Science Foundation Center for Biological Timing, and 4 Department of Pediatrics, University of Virginia, Charlottesville 22908; and 5 INSMED Pharmaceuticals, Richmond, Virginia 23219


    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Clarifying the time evolution, and underlying neuroendocrine regulation, of hormone secretion during puberty is of pivotal importance both physiologically and clinically. We sought to determine whether clinical growth and elevated growth hormone (GH) release in pubertal boys are associated with shifts in the irregularity of GH secretory patterns. We studied GH release in 23 healthy boys longitudinally at ~4-mo intervals over a 6-yr period throughout puberty, by repetitive blood sampling at 20-min intervals for 24 h on each study date. To quantify serial irregularity in the GH profiles, we utilized approximate entropy (ApEn), a scale- and model-independent quantification of the extent of sequential "randomness." Complementary statistical analyses establish that on a per-subject basis, serum GH concentration-time series show greatest secretory disorderliness (maximal ApEn) in prepuberty and mid-to-late puberty, followed by a steep decline in ApEn to maximal orderliness in postpubertal adolescence (P < 0.0001, ANOVA). Pooling all subject comparisons, we observed a persistent positive correlation between ApEn and growth velocity (GrVel), Pearson r = 0.467 (P < < 10-10). Moreover, in general, ApEn(GH) preceded GrVel evolution, with a time frame lead of 4-8 mo providing the most pronounced correlations between ApEn and GrVel. In the setting of low postpubertal growth, per-subject ApEn values fell to approximately one-half of maximal ApEn values and, on average, were in the 13th percentile from minimal (P < < 10-10) for fully mature boys. Thus, in a longitudinal transpubertal analysis, greater irregularity in GH secretion corresponds to greater linear growth in boys, which culminates in highly regular GH secretory dynamics after sexual maturity. In addition to clinical utility, there may be added experimental merit in knowing that GH dynamics typically predict accelerated linear growth in 4-8 mo.

approximate entropy; pulsatility; growth rate; somatotropic axis; humans


    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

GREATER UNDERSTANDING of the regulation of hormonal dynamics during puberty is of vital importance both physiologically and clinically. Whereas many studies have quantified the amounts of hormone secreted by the somatotropic and gonadotropic axes during this transition (1, 5, 7, 25, 30, 31), basic neuroregulatory questions remain to be clarified about the underlying control mechanisms that drive increased hormone output (3, 24). One means of assessing pulsatility changes on the basis of hormonal time-series data is approximate entropy (ApEn), a model-independent quantification of sequential irregularity (10, 21). ApEn applications have allowed the endocrinologist to monitor and quantify changes in secretory patterns from a perspective complementary to that provided by output estimates or pulse detection algorithms (26). The present study explores the hypothesis that the marked pubertal enhancement of growth hormone (GH) release in boys is accompanied by shifts in the irregularity (extent of "randomness") of GH secretory patterns, thus identifying a distinct aspect of changes in GH system control. Given that ApEn changes bear a mechanistic relationship to alterations in feedback control (discussed further below), the ensuing analyses thus allow inferential assessment of changing feedback regulation in the somatotropin release-inhibiting hormone (SRIH)/GH-releasing hormone (GHRH)-somatotropic (GH)-insulin-like growth factor I (IGF-I) axis across human puberty.

In a recent cross-sectional study of 53 boys (29), we observed that ApEn was largest (GH dynamics most irregular) at mid-to-late puberty, with a return to maximal orderliness (smallest ApEn) thereafter toward reproductive maturity. However, this study was significantly limited by its cross-sectional study design and absence of concurrent growth analysis. Thus causality of per-subject ApEn evolution could at best be inferred only indirectly and could not be related to individual growth activity. Critically from a statistical perspective, cross-subject heterogeneity often blurs longitudinal evolutionary determinations that are more apparent on a per-subject basis. Herein, we wished to determine longitudinal per-subject profiles of ApEn evolution of GH across puberty, and to assess the time association (if any) between ApEn changes and concurrently measured linear growth rate. For this purpose, we utilize the clinical model to be described in METHODS on the basis of GH profiles collected from 23 boys studied longitudinally at 4-mo intervals over a 6-yr period as they progressed through puberty. These data provide an uncommon, possibly unprecedented, 72-point GH time series per study date for each of ~18 study dates for each of the 23 boys.

In addition to clarifying the evolution of ApEn transpubertally, we examine the association between ApEn values and a classical primary clinical marker of pubertal change, growth velocity (GrVel), considering both lagged and unlagged comparisons. GrVel is a well-established and relevant biological response (22), which can be measured quite accurately via standard auxological methods.


    METHODS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Clinical Protocols

This study was approved by the University of Virginia Human Investigation Committee and performed in the General Clinical Research Center. Written informed consent was obtained from a parent or guardian, and written assent was obtained from each subject. A longitudinal study was conducted of 23 prepubertal and early pubertal normally growing boys. Each volunteer was observed for >= 6 yr as he entered and progressed through pubertal development, i.e., from prepuberty, as defined by Tanner Stage I for pubic hair and genitalia, to early postpuberty (clinical Tanner Stage V for the same characteristics). At entry, each boy was within 2 SD of height for age and was growing normally. Each boy also had a normal weight for height. None of the boys had chronic illnesses or was receiving medications chronically. Data were gathered from each boy approximately every 4 mo. Height, from which GrVel was calculated, was taken mainly by a single experienced technician by use of the mean of three early morning measurements with a Harpenden stadiometer. All measurements were made to the nearest millimeter. All volunteers underwent blood sampling at 20-min intervals for 24 h on each study date (n = 72 point time series). Serum GH concentrations were assayed in duplicate by immunoradiometric assay, or IRMA (Nichols Laboratories, San Juan Capistrano, CA) (6). The minimum detectable concentration was 0.5 µg/l. A detailed description of the mean intra-assay coefficients of variation (CV) across a full range of serum GH concentrations has been reported previously (4).

Quantification of Irregularity

To quantify irregularity, we utilize ApEn, a model-independent statistic defined in Ref. 10, with further mathematical properties and representative biological (including endocrinological) applications given in Refs. 8, 9, 12, 13, 15, 17, 19-21, 23, and 28. ApEn assigns a nonnegative number to a time series, with larger values corresponding to greater apparent process randomness or serial irregularity and smaller values corresponding to more instances of recognizable features or patterns in the data. Two input parameters, a run length m and a tolerance window r, must be specified to compute ApEn. Briefly, ApEn measures the logarithmic likelihood that runs of patterns that are close (within r) for m contiguous observations remain close (within the same tolerance width r) on next incremental comparisons; the precise mathematical definition is given in Ref. 10.

For this study, we calculated ApEn values for all data sets, m = 1 and r = 20% of the standard deviation (SD) of the individual subject time series. Normalizing r to each time- series SD gives ApEn a translation and scale invariance to absolute serum concentration levels (12). Additionally, ApEn provides a direct barometer of feedback system change in many coupled systems (11, 18). Further technical discussion of mathematical and statistical properties of ApEn can be found elsewhere (14, 16).

Analytic Strategy

Challenges in forming statistical tests. The goal of this study is to deduce typical and general characteristics of longitudinal per-subject ApEn(GH) curves transpubertally. This requires a careful analytic strategy for several reasons. First, the primary per-subject inputs are ApEn and GrVel curves (longitudinal point sets), rather than single numerical values; thus a means to aggregate these correlated data for suitable statistical testing is not straightforward. Reducing each per-subject curve to a several-parameter set via a low-order best-fit polynomial (or other elementary function) seems inappropriate, because many of the per-subject longitudinal data sets are highly nonmonotonic; i.e., they show numerous changes in direction (as in Fig. 1). Virtually no low-parameter models are characterized by such frequent shifts in direction (manifesting many first derivative sign changes) across time. Second, pooling the data is not so straightforward, even upon z-transformation standardization (see below), because subjects were not sampled on common dates and had differing start and end study dates relative to their date of maximal GrVel and unequal numbers of studies. Third, GrVel itself does not smoothly increase to a maximal value, thereafter subsequently showing monotonic decrease (see Fig. 1), so that comparisons of ApEn values to the clinical reference GrVel must at least accommodate this irregular variation.


View larger version (31K):
[in this window]
[in a new window]
 
Fig. 1.   Representative longitudinal approximate entropy [ApEn; run length (m) = 1; tolerance window (r) = 20% SD] of 24-h growth hormone (GH) time series and concomitant growth velocity (GrVel) values for 4 boys as each progressed through puberty. Data were obtained approximately every 4 mo.

To address these issues, we chose to analyze the data via several different and complementary, although thematically related, assessments to establish robustness and generality to the findings.

Statistical tests: ApEn(GH) evolution. On the basis of previous cross-sectional study (29), we hypothesized that per-subject ApEn values would be at a midlevel prepubertally, would tend to increase toward maximal ApEn at or near the time of maximal GrVel, and then generally and steadily decrease to a very low value at maturity, markedly smaller than the initial study value, indicating very regular GH dynamics at the final test date. To evaluate this notion in a robust manner that did not require model-based assumptions for the ApEn longitudinal curves, we performed the following calculations. First, for a given subject, denote the ApEn(GH) value on the ith study date as ApEni and denote the GrVel value on the ith study date as GrVeli. Denote the study date for which GrVel is maximal for the subject as tmax, with corresponding ApEn value ApEnvmax; also, denote by t1 and tfin the subject study start and final dates, respectively, with corresponding ApEn values denoted ApEn1 and ApEnfin. Last, define tearly as the study date midway between t1 and tmax, and define tlate as the study date midway between tmax and tfin. Thus the five study dates t1, tearly, tmax, tlate, and tfin can be seen as a per-subject increasing time sequence of markers segmenting the start date-maximal growth date-final date continuum.

One-way ANOVA was then applied to evaluate the null hypothesis that the ApEn values were statistically unrelated to pubertal development across the five epochs t1, tearly, tmax, tlate, and tfin. ANOVA was applied to both the raw ApEn values and to the standardized (z-transformed) ApEn values.

To further appraise the longitudinal ApEn trajectory, we also calculated 1) the percentage of dates i < tmax for which ApEni > ApEn1, and 2) the percentage of dates i > tmax for which ApEni < ApEnvmax. Both of these calculations were performed per subject and were subsequently pooled. Calculation 1 assesses the extent and consistency of increases in GH irregularity from study onset to the time of maximal GrVel, and calculation 2 assesses the extent and consistency of decreases in GH irregularity from the time of maximal GrVel to study completion.

Statistical tests: joint ApEn(GH), GrVel associations. We performed several sets of calculations to determine linear correlation between ApEn and GrVel. First, both per subject and subsequently pooled, we calculated the correlation between ApEni and GrVeli..Next, we performed led and lagged correlations and calculated the correlation between ApEni and GrVeli-N for each of N = 1, 2, 3, -1, -2, and -3.

The above overall correlation calculations have the limitation that they are affected by cross-subject variation and heterogeneity. One major hypothesis of this study is that, on a per-subject basis, ApEn correlates strongly and positively with GrVel. To assess this hypothesis by using the power of the pooled data base, yet in a manner free of cross-subject heterogeneity, we also calculated pooled overall correlations between standardized (z-transformed) ApEni and GrVeli values. The z-transformations were performed on each subject's ApEn and GrVel values. We also performed led and lagged correlations of the pooled z-transformed data, identically to the nontransformed comparisons, for the lead and/or lag times described.

Low final date ApEn values. We tested the postulate that, and determined the degree to which, final date ApEn values are low, relative to the maximal growth date on a per-subject analysis, in a robust, nonparametric manner as follows. First, rank for each subject the final test date ApEn value among all (N) ApEni values for that subject, denoting this as Aprankfin, which has a value between 1 (lowest) and N (largest). Then normalize this ranking by nAprankfin = (Aprankfin - 1)/(N - 1), producing a range between 0 and 1. Under the null hypothesis of identical ApEn distributions on all dates, the mean nAprankfin should have a mean of 0.5 (i.e., the final date ApEn value would on average equal the median subject ApEn value). A value of nAprankfin substantially below 0.5 would indicate a final date ApEn value that is among the lowest of that subject's ApEn values.

We also calculated both the percentage of subjects for which ApEn1 > ApEnfin and the percentage of subjects for which ApEnvmax > ApEnfin, as another pair of assessments to evaluate the hypothesized very regular GH dynamics at the final test date.

Statistical Analysis

All statistical comparisons below for discrimination between two groups employed the two-sided t-test with unknown variance. Linear correlation was assessed via the Pearson r statistic. One-way ANOVA was performed as we have discussed.


    RESULTS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Representative longitudinal ApEn(GH) and GrVel plots for four boys are shown in Fig. 1. In Fig. 2, representative GH serum concentration-time series are shown for the five indicated study dates for the boy whose longitudinal ApEn values are depicted in Fig. 1A.


View larger version (26K):
[in this window]
[in a new window]
 
Fig. 2.   Representative serum GH concentration (µg/l) time series over 24 h measured in blood collected at 20-min intervals for the 5 study dates t1, tearly, tmax, tlate, and tfin for the individual boy whose longitudinal ApEn values are depicted in Fig. 1A. ApEn values indicated are ApEn(1, 20% SD): A: t1, age 10.1 yr, Tanner stage 1, ApEn = 0.582, GrVel = 2.88 cm/yr; B: tearly, age 11.6 yr, Tanner stage 1.5, ApEn = 1.030, GrVel = 6.70 cm/yr; C: tmax, age 12.9 yr, Tanner stage 3.5, ApEn = 0.776, GrVel = 10.34 cm/yr; D: tlate, age 15.2 yr, Tanner stage 5, ApEn = 0.365, GrVel = 3.70 cm/yr; E: tfin, age 16.6 yr, Tanner stage 5, ApEn = 0.088, GrVel = 1.30 cm/yr.

As summarized in Table 1 and displayed in Fig. 3A, the comparison of raw ApEn values at the five study dates t1,tearly,tmax, tlate, and tfin showed the largest irregularity at tearly and a minimum ApEn (greatest orderliness) at tfin (ANOVA, P < 0.0001). The corresponding analyses based on the standardized (z-transformed) ApEn values show the same general pattern of evolution (Table 1, Fig. 3B), with the same significance (P < 0.0001). In addition, the standardized data analysis indicates a slight, although not statistically significant, rise in ApEn values at tmax compared with t1, following the maximum at tearly.

                              
View this table:
[in this window]
[in a new window]
 
Table 1.   GH ApEn values associated with 5 study times



View larger version (12K):
[in this window]
[in a new window]
 
Fig. 3.   Histograms of mean GH ApEn(1, 20% SD) values associated with the 5 study dates t1, tearly, tmax, tlate, and tfin (see text for definition of these per-subject dates), based on raw ApEn values (A) or standardized ApEn values (B). Precise values are indicated in Table 1. By ANOVA, these 5 dates do not have identical ApEn distributions (P < 0.0001). Also, for each plot, columns identified with different super- or subscripts differ significantly (P < 0.0001).

To further evaluate longitudinal curve characteristics, assessing increases in GH irregularity from study onset to the time of maximal GrVel showed that, pooled, ApEni > ApEn1 for 64% of dates i < tmax. Assessment of decreases in GH irregularity from the time of maximal GrVel to study completion indicated that, pooled, ApEni < ApEnvmax for 90% of dates i > tmax.

Table 2 indicates the linear correlation between ApEni and GrVeli-N for each of N = 1, 2, 3, -1, -2, and -3 for the pooled ensemble data, both raw and standardized. All the indicated correlations are positive, with high significance (P < 10-10 for all but the lag +3 correlations). For both standardized and raw data comparisons, the largest positive correlation between ApEni and GrVeli-N was seen for the lag -1 data (standardized Pearson r = 0.510, raw data Pearson r = 0.494), with next largest positive correlation seen for the lag -2 data (standardized Pearson r = 0.487, raw data Pearson r = 0.467). Interestingly, both lag -1 and lag -2 data show more pronounced positive correlation than lag 0, and, importantly, all negative lag correlations (ApEn leading GrVel comparisons) were larger than any positive lag correlation.

                              
View this table:
[in this window]
[in a new window]
 
Table 2.   Lagged correlation strengths (Pearson r) for the relationship between GH ApEn and growth velocity

In general, on a per-subject basis, correlation is generally moderately large and positive for each of the studied lags. Table 3 summarizes the per-subject correlations of ApEn to GrVel at the indicated lags. For time lags -3, -2, -1, and 0, 91% of the correlations were positive; for time lag 0, 87% were positive; for time lag +2, 83% were positive; and for time lag +3, 65% of the per-subject correlations were positive. From an alternative viewpoint, the percentages of subjects with linear correlations that were >= 0.4 (and positive) were 74% for lags -1 and -2; 65% for lags 0 and +1; 56% for lag -3; 48% for lag +2; and 13% for lag +3.

                              
View this table:
[in this window]
[in a new window]
 
Table 3.   Summary of per-subject lagged correlations between GH ApEn and growth velocity in boys

As another calculation pertinent to determining whether ApEn leads or follows GrVel, define tApmax as the study date with maximal ApEn value. We found that tApmax < tmax for 18 of 23 subjects, indicating that maximal ApEn values occurred earlier than maximal GrVel in 78% of subjects. As well, for three of the remaining five subjects, tApmax is either only one or two test periods (4-8 mos) delayed from tmax.

Final date ApEn values were quite low. As seen in Table 1, cross-sectional ApEnfin = 0.380 ± 0.152 is approximately one-half of ApEn at tearly = 0.786 ± 0.186. Thus the difference in the relative orderliness of GH release between early puberty to midpuberty and maturity is great. This is suggested by and consistent with the longitudinal plots shown as Fig. 1. However, the degree to which final-date ApEn values are low (GH secretion quite regular) relative to maximal growth date is most evident by per-subject analysis via calculation of the nAprankfin values described above. The mean nAprankfin (over the 23 subjects) of 0.1269, i.e., the 13th percentile, is smaller (P < 10-10) than the expected mean nAprankfin of 0.5 under the null hypothesis. Moreover, for 9 of the 23 subjects, the final date had the lowest ApEn, whereas for five additional subjects, the final date was either second or third lowest.

Complementary evidence of the apparently very regular GH dynamics at the final test date is that ApEn1 > ApEnfin for 96% of subjects, whereas ApEnvmax > ApEnfin for 100% of subjects.


    DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Summarizing the results, we have established that on a per-subject basis, 24-h serum GH concentration time series show pronounced variation in their sequential irregularity transpubertally, with maximal secretory disorderliness (ApEn) in pre- and mid-to-late puberty, followed by a steep decline in ApEn to maximal orderliness in postpubertal adolescence. Moreover, we determined that, in general, ApEn(GH) elevations precede GrVel increases, with a time-frame lead of one to two study periods (4-8 mo) providing the most pronounced correlations, as seen in Tables 2 and 3. This time delay is also consistent with visual inference given by Fig. 1, and by the finding that for 78% of boys, the date of maximal ApEn preceded that of maximal GrVel. As discussed below, there may be added clinical and experimental utility in knowing that high ApEn(GH) values typically predict high GrVel values in 4-8 mo, beyond the strictly diagnostic or confirmatory perspective that positive unlagged correlation conveys.

Statistically, per-subject standardization sharpens the (lagged and led) correlations between ApEn and GrVel by a modest, albeit significant, amount compared with raw data comparisons. More notably, one critical utility of the longitudinal study design applied here is to clarify that the hormonal dynamics precede the (accelerated) physical growth. In contrast, whereas the cross-sectional study in Ref. 29 gives qualitatively relatively similar results regarding ApEn(GH) evolution with increasing age (compare Fig. 1 of Ref. 29 with Fig. 3 here), the data available in Ref. 29 preclude any direct comparison of hormonal dynamics, e.g., ApEn, with either GrVel or an alternative measure of physical growth.

Of course, it is imperative to note that growth is a multifactor-dependent phenomenon, and that the present analysis reflects only one correlate (GH), albeit an important one. We anticipate that, ultimately, forecasts based on multiple physiological (including hormonal) characteristics will provide enhanced accuracy in predicting periods of very rapid linear growth. To reinforce this point, in modeling the growth velocity of these boys, we have found marked effects of testosterone concentration in addition to the mean GH level (32). In addition to age (either chronological or bone), the interaction between the mean GH level and testosterone concentration has a significant effect on the GrVel. The maximal rise in circulating testosterone concentration (and perhaps its conversion to biologically effective estrogens) precedes the maximal mean GH concentration and peak GrVel by ~3 mo (32), commensurate with the findings of the present study.

The longitudinal analyses allow a first-order, clinically applicable, predictive relationship between ApEn and GrVel. Across all subjects, an ApEn value of >0.8 for the 24-h GH time series (which occurred in ~20% of measurements) produced a mean GrVel in the subsequent measurement period of 7.33 ± 2.95 cm/yr, which is significantly larger than the overall mean GrVel, which is 5.43 ± 3.36 cm/yr (P < 0.0001). The GrVel of 7.33 cm/yr is the 68th percentile of pooled GrVel values. The decision rule that an ApEn value >0.8 generally predicts a high GrVel in the next 4-mo period is relatively easy to implement in potential applications, in that it is an overall population rule (rather than a per-subject determinant), thus incorporating cross-subject variation.

ApEn values of GH time series on each boy's final test date (0.380 ± 0.152) were insignificantly different (P > 0.4) from comparable 20-min sampled GH ApEn values (0.425 ± 0.183) in healthy young men (mean age 25 yr; age range 22-28 yr) reported in Ref. 13. Thus we infer that, from a secretory irregularity perspective, in the present analysis, healthy postpubertal boys exhibit a fully mature GH subnetwork by the mean per-subject final study age of 17.3 (±0.9) yr.

One issue should be clarified concerning statistical study design. For statistical reasons, we chose to avoid fitting unimodal functions (or splines) to either the ApEn or GrVel longitudinal data. In principle, if such best-fit functions were established, then characteristics of each curve, such as the date of maximal growth rate from such a GrVel curve, could be inferred algebraically. However, given the nonmonotonicity of each of ApEn and GrVel as a function of time (see comments in Analytical Strategy), we expect that such fitting procedures could produce somewhat equivocal or arbitrary results. Indeed, the ApEn-GrVel correlations as assessed herein provided a robust, and we believe biologically realistic, quantification of the association of ApEn and transpubertal growth.

Neuroendocrine Considerations

We have previously linked ApEn changes to a mechanistic understanding via theoretical modeling studies (11, 18). In particular, for diverse mathematical network models, greater regularity (lower ApEn) typically corresponds to greater component and subsystem autonomy and, conversely, greater irregularity corresponds to increased external influences, increased coupling strength, and/or accelerated positive feedback. On the basis of the present analysis, we thus propose that the increased GH irregularity observed before and near peak GrVel reflects more critical interacting factors within the GH feedback axis, or a greater intensity of particular interactions, e.g., among GHRH-SRIH-GH-IGF-I. Specific, yet distinct, mechanistic explanations consistent with this interpretation are proposed and discussed in Ref. 29, although as noted there, further clinical and interventional studies will be required to distinguish among several theoretically plausible alternatives. Additionally, the very low ApEn values seen at maturity are consistent with a relatively closed GH subnetwork, in which some of the inter-system interactions that are enhanced in puberty have diminished substantially.

Two additional clinical investigations (29) support the hypothesis that estradiol in particular, or testosterone acting after its aromatization to estradiol, can govern the feedback regulation of GH release during puberty and in young adulthood. In the first study, clinical experiments in prepubertal girls with gonadal dysgenesis (Turner's syndrome) showed that treatment with a small dose of ethinyl estradiol orally for 5 wk increases GH ApEn significantly. In the second study, ApEn of GH increased in response to short-term parenteral testosterone administration to clinically prepubertal boys with constitutionally delayed adolescence, yet showed no increase upon treatment with a nonaromatizable androgen, alpha -dihydrotestosterone (DHT).

A minority of children born with intrauterine growth retardation have a propensity for the metabolic syndrome as adults (2). Theoretically, the disordered intrauterine environment had reprogrammed intermediary metabolism. Large doses of GH given for 2 yr leads to marked catch-up growth in the majority of those who fail to catch up on their own. Those that do catch up to their target height percentile often remain along that growth trajectory even without additional GH. One interpretation of this observation is that the metabolic syndrome is again reprogrammed toward normalcy by this treatment. It is likely that altering GH secretion is pivotal in this extrauterine treatment phase. Analyses such as those discussed in this report may then provide explicit numerical support of the hypothesis that feedback mechanisms have been altered in this metabolic reprogramming.

Potential Applications and Follow-Up Studies

Some follow-up studies are suggested by the present findings. For example, clinically one would like to predict, as well as possible, stature at maturity, and as an important intermediate point, the time of maximal linear growth. In addition to GH input, one might develop a multivariate predictive model utilizing concurrent (e.g., testosterone, IGF-I, and/or estradiol) measurements when they are available. Furthermore, if time-series data per study date became available for two or more hormones, e.g., GH and testosterone, one could calculate cross-ApEn values (19, 21) on paired hormonal time series to determine changes in bihormonal synchrony across the peripubertal years. Such information could aid in determining effective intervention modalities, e.g., whether and when treatment with testosterone would be necessary to augment linear growth.

Additionally, there may be practical applicability in determining whether administering GH on the basis of some secretory irregularity threshold, e.g., current ApEn(GH) >0.8, achieves more efficient growth per unit GH dosage compared with dosages given either at random times or uniformly over a longer regimen.

Whether the ApEn variations in GH across pubertal development are unique to the somatotropic axis, or also apply to other hormonal subsystems transpubertally, is not yet known. An earlier finding is that aging well beyond the young adult years results in a slowly increasing rise in GH ApEn (28), indicating a progressively more disorderly GH release pattern, also seen for luteinizing hormone/testosterone (19), ACTH/cortisol (27), and insulin (8). Hence, one could hypothesize that ApEn changes in hormonal secretory patterns may in part detect broader alterations in neuroregulatory control.

The present study was conducted in healthy boys only. Accordingly, the extent to which analogous secretory irregularity changes are present in boys with growth-related pathophysiology, in young girls in the peripubertal time frame, and in nonhuman species remains a very important, yet open, area for future research.


    ACKNOWLEDGEMENTS

We thank Sandra Jackson and the nursing staff at the General Clinical Research Center for their expert patient care, Melanie Bishop Harlow and Margaret Wood Ball for research coordination, and Catherine Kern and Ginger Bauler for technical assistance with the assays. Our long-term colleague Dr. Robert M. Blizzard is acknowledged for foresight and for development of the concept and the patient protocol. We thank a large number of pediatric endocrine fellows at the University of Virginia for the clinical care of the boy subjects. Without the remarkable generosity of these boys and their parents, the studies could not have been completed.


    FOOTNOTES

This work was supported in part by National Institutes of Health (NIH) NCRR Grant RR-00847 (to the Clinical Research Center of the University of Virginia); NIH National Institute of Child Health and Human Development Research Career Development Award 1-KO4-HD-00634 (to J. D. Veldhuis); the University of Virginia Pratt Foundation; the National Science Foundation (NSF) Center for Biological Timing (NSF Grant DIR-89-20162); NIH U-54 Specialized Cooperative Centers Program for Reproduction Research (NICHD) HD-28934; Postdoctoral Research Training in Diabetes and Hormone Grants 5T32-DK-07320 and DK-07642 (to the University of Virginia); NIH National Institute on Aging R01 AG-14799 and R03 AG-14873 (to J. D. Veldhuis); and NIH R01 HD-32631 (to A. D. Rogol).

Address for reprint requests and other correspondence: S. M. Pincus, 990 Moose Hill Road, Guilford, CT 06437 (E-mail:stevepincus{at}alum.mit.edu).

The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. §1734 solely to indicate this fact.

Received 15 December 1999; accepted in final form 21 March 2000.


    REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

1.   Attie, KM, Ramirez NR, Conte FA, Kaplan SL, and 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 71: 975-983, 1990[Abstract].

2.   Barker, DJ, Hales CN, Fall CH, Osmond C, Phipps K, and Clark PM. Type 2 (non-insulin-dependent) diabetes mellitus, hypertension, and hyperlipidaemia (syndrome X): relation to reduced fetal growth. Diabetologia 36: 62-67, 1993[ISI][Medline].

3.   Giustina, A, and Veldhuis JD. Pathophysiology of the neuroregulation of GH sceretion in experimental animals and the human. Endocr Rev 19: 717-797, 1998[Abstract/Free Full Text].

4.   Hartman, ML, Veldhuis JD, Vance ML, Faria ACS, Furlanetto RW, and Thorner MO. Somatotropin pulse frequency, basal concentrations are increased in acromegaly and are reduced by successful therapy. J Clin Endocrinol Metab 70: 1375-1384, 1990[Abstract].

5.   Keenan, BS, Richards GE, Ponder SW, Dallas JS, Nagamani M, and Smith ER. Androgen-stimulated pubertal growth: the effects of testosterone and dihydrotestosterone on growth hormone and insulin-like growth factor-I in the treatment of short stature and delayed puberty. J Clin Endocrinol Metab 76: 996-1001, 1993[Abstract].

6.   Martha, PM, Jr, Goorman KM, Blizzard RM, Rogol AD, and Veldhuis JD. Endogenous growth hormone secretion and clearance rates in normal boys as determined by deconvolution analysis: relationship to age, pubertal status and body mass. J Clin Endocrinol Metab 74: 336-344, 1992[Abstract].

7.   Mauras, N, Rogol AD, Haymond MW, and Veldhuis JD. Sex steroids, growth hormone, IGF-1: neuroendocrine and metabolic regulation in puberty. Horm Res 45: 74-80, 1996[ISI][Medline].

8.   Meneilly, GS, Ryan AS, Veldhuis JD, and Elahi D. Increased disorderliness of basal insulin release, attenuated insulin secretory burst mass, and reduced ultradian rhythmicity of insulin secretion in older individuals. J Clin Endocrinol Metab 82: 4088-4093, 1997[Abstract/Free Full Text].

9.   Morrison, S, and Newell KM. Inter- and intra-limb coordination in arm tremor. Exp Brain Res 110: 455-464, 1996[ISI][Medline].

10.   Pincus, SM. Approximate entropy as a measure of system complexity. Proc Natl Acad Sci USA 88: 2297-2301, 1991[Abstract].

11.   Pincus, SM. Greater signal regularity may indicate increased system isolation. Math Biosci 122: 161-181, 1994[ISI][Medline].

12.   Pincus, SM, Cummins TR, and Haddad GG. Heart rate control in normal and aborted-SIDS infants. Am J Physiol Regulatory Integrative Comp Physiol 264: R638-R646, 1993[Abstract/Free Full Text].

13.   Pincus, SM, Gevers EF, Robinson ICAF, van den Berg G, Roelfsema F, Hartman ML, and Veldhuis JD. Females secrete growth hormone with more process irregularity than males in both humans and rats. Am J Physiol Endocrinol Metab 270: E107-E115, 1996[Abstract/Free Full Text].

14.   Pincus, SM, and Goldberger AL. Physiological time-series analysis: what does regularity quantify? Am J Physiol Heart Circ Physiol 266: H1643-H1656, 1994[Abstract/Free Full Text].

15.   Pincus, SM, Hartman ML, Roelfsema F, Thorner MO, and Veldhuis JD. Hormone pulsatility discrimination via coarse and short time sampling. Am J Physiol Endocrinol Metab 277: E948-E957, 1999[Abstract/Free Full Text].

16.   Pincus, SM, and Huang WM. Approximate entropy: statistical properties and applications. Commun Statist-Theory Methodol 21: 3061-3077, 1992.

17.   Pincus, S, and Kalman RE. Not all (possibly) "random" sequences are created equal. Proc Natl Acad Sci USA 94: 3513-3518, 1997[Abstract/Free Full Text].

18.   Pincus, SM, and Keefe DL. Quantification of hormone pulsatility via an approximate entropy algorithm. Am J Physiol Endocrinol Metab 262: E741-E754, 1992[Abstract/Free Full Text].

19.   Pincus, SM, Mulligan T, Iranmanesh A, Gheorghiu S, Godschalk M, and Veldhuis JD. Older males secrete luteinizing hormone (LH) and testosterone more irregularly, and jointly more asynchronously, than younger males. Proc Natl Acad Sci USA 93: 14100-14105, 1996[Abstract/Free Full Text].

20.   Pincus, SM, Padmanabhan V, Lemon W, Randolph J, and Midgley AR. FSH is secreted more irregularly than LH in both human and sheep. J Clin Invest 101: 1318-1324, 1998[Abstract/Free Full Text].

21.   Pincus, S, and Singer BH. Randomness and degrees of irregularity. Proc Natl Acad Sci USA 93: 2083-2088, 1996[Abstract/Free Full Text].

22.   Preece, MA, and Baines MJ. A new family of mathematical models describing the human growth curve. Ann Hum Biol 5: 1-24, 1978[ISI][Medline].

23.   Radhakrishnan, N, and Gangadhar BN. Estimating regularity in epileptic seizure time-series data. IEEE Eng Med Biol 17: 89-94, 1998[ISI].

24.   Rogol, AD, Martha PM, Jr, Johnson ML, Veldhuis JD, and Blizzard RM. Growth hormone secretory dynamics during puberty. In: The Somatotrophic Axis and the Reproductive Process in Health and Disease, edited by Adashi EY, and Thorner MO. New York: Springer-Verlag, 1996, p. 69-82.

25.   Thompson, RG, Rodriguez A, Kowarski AA, Migeon CJ, and Blizzard RM. Integrated concentrations of growth hormone correlated with plasma testosterone and bone age in preadolescent and adolescent males. J Clin Endocrinol Metab 62: 1341-1344, 1978.

26.   Urban, RJ, Evans WS, Rogol AD, Kaiser DL, Johnson ML, and Veldhuis JD. Contemporary aspects of discrete peak-detection algorithms. 1. The paradigm of the luteinizing hormone pulse signal in men. Endocr Rev 9: 3-37, 1988[ISI][Medline].

27.   Van den Berg, G, Pincus SM, Veldhuis JD, Frolich M, and Roelfsema F. Greater disorderliness of ACTH and cortisol release accompanies pituitary-dependent Cushing's disease. Eur J Endocrinol 136: 394-400, 1997[ISI][Medline].

28.   Veldhuis, JD, Liem AY, South S, Weltman A, Weltman J, Clemmons DA, Abbott R, Mulligan T, Johnson ML, Pincus S, Straume M, and Iranmanesh A. Differential impact of age, sex-steroid hormones, and obesity on basal vs. pulsatile growth hormone secretion in men as assessed in an ultrasensitive chemiluminescence assay. J Clin Endocrinol Metab 80: 3209-3222, 1995[Abstract].

29.   Veldhuis, JD, Metzger DL, Martha PM, Mauras N, Kerrigan JR, Keenan B, Rogol AD, and Pincus SM. Estrogen and testosterone, but not a nonaromatizable androgen, direct network integration of the hypothalamo-somatotrope (GH)-IGF-1 axis in the human: evidence from pubertal pathophysiology and sex-steroid hormone replacement. J Clin Endocrinol Metab 82: 3414-3420, 1997[Abstract/Free Full Text].

30.   Wehrenberg, WB, and Giustina A. Basic counterpoint: mechanisms and pathways of gonadal steroid modulation of growth hormone secretion. Endocr Rev 13: 299-308, 1992[Abstract].

31.   Zachmann, M. Interrelations between growth hormones and sex hormones: physiology and therapeutic consequences. Horm Res 1: 1-8, 1992.

32.  Zhang J, Peddada SD, Malina RM, and Rogol AD. Longitudinal assessment of hormonal and physical alterations during normal puberty in boys. VI. Modeling of growth velocity, mean growth hormone (GH mean) and serum testosterone concentrations. Am J Hum Biol In press.


Am J Physiol Endocrinol Metab 279(2):E417-E424
0193-1849/00 $5.00 Copyright © 2000 the American Physiological Society