Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Prolactin Secretion in Healthy Adults Is Determined by Gender, Age and Body Mass Index

  • Ferdinand Roelfsema ,

    f.roelfsema@lumc.nl

    Affiliation Department of Endocrinology and Metabolic Diseases, Leiden University Medical Center, Leiden, The Netherlands

  • Hanno Pijl,

    Affiliation Department of Endocrinology and Metabolic Diseases, Leiden University Medical Center, Leiden, The Netherlands

  • Daniel M. Keenan,

    Affiliation Department of Statistics, University of Virginia, Charlottesville, Virginia, United States of America

  • Johannes D. Veldhuis

    Affiliation Endocrine Research Unit, Mayo Medical and Graduate Schools, Clinical Translational Research Center, Mayo Clinic, Rochester, Minnesota, United States of America

Correction

30 Sep 2014: The PLOS ONE Staff (2014) Correction: Prolactin Secretion in Healthy Adults Is Determined by Gender, Age and Body Mass Index. PLOS ONE 9(9): e110203. https://doi.org/10.1371/journal.pone.0110203 View correction

Abstract

Background

Prolactin (PRL) secretion is quantifiable as mean, peak and nadir PRL concentrations, degree of irregularity (ApEn, approximate entropy) and spikiness (brief staccato-like fluctuations).

Hypothesis

Distinct PRL dynamics reflect relatively distinct (combinations of) subject variables, such as gender, age, and BMI.

Location

Clinical Research Unit.

Subjects

Seventy-four healthy adults aged 22–77 yr (41 women and 33 men), with BMI 18.3–39.4 kg/m2.

Measures

Immunofluorometric PRL assay of 10-min samples collected for 24 hours.

Results

Mean 24-h PRL concentration correlated jointly with gender (P<0.0001) and BMI (P = 0.01), but not with age (overall R2 = 0.308, P<0.0001). Nadir PRL concentration correlated with gender only (P = 0.017) and peak PRL with gender (P<0.001) and negatively with age (P<0.003), overall R2 = 0.325, P<0.0001. Forward-selection multivariate regression of PRL deconvolution results demonstrated that basal (nonpulsatile) PRL secretion tended to be associated with BMI (R2 = 0.058, P = 0.03), pulsatile secretion with gender (R2 = 0.152, P = 0.003), and total secretion with gender and BMI (R2 = 0.204, P<0.0001). Pulse mass was associated with gender (P = 0.001) and with a negative tendency to age (P = 0.038). In male subjects older than 50 yr (but not in women) approximate entropy was increased (0.942±0.301 vs. 1.258±0.267, P = 0.007) compared with younger men, as well as spikiness (0.363±0.122 vs. 0463±2.12, P = 0.031). Cosinor analysis disclosed higher mesor and amplitude in females than in men, but the acrophase was gender-independent. The acrophase was determined by age and BMI (R2 = 0.186, P = 0.001).

Conclusion

In healthy adults, selective combinations of gender, age, and BMI specify distinct PRL dynamics, thus requiring balanced representation of these variables in comparative PRL studies.

Introduction

Prolactin (PRL) is a 23 kDa protein secreted by the pituitary gland and has many actions of which the lactotrophic function is best known. This hormone has various other effects on reproduction, metabolism and tumorigenicity mainly described in animals, but the precise role in the human is less well established [1], [2].

Prolactin secretion proceeds via combined pulsatile (burst-like) and basal (time-invariant) modes of release. A complicating issue in defining normative ranges even for measures of PRL secretion is that they may depend upon one or more biological or clinical factors, such as gender, age, BMI, sex-steroid concentrations, core temperature, nutrition, stress, exercise, medications and renal disease [1], [3][9]. Prolactin secretion is primarily regulated by the inhibitory action of hypothalamic dopamine, and by ultrashort autofeedback, but the physiological role of various releasing hormones is not established in man [2]. These factors would putatively determine more complex PRL dynamics, which arise physiologically from feedforward (stimulatory) and feedback (inhibitory) signals interacting in an integrative fashion, as demonstrated in detail for GH and LH in men [10]. Novel integrative measures are approximate entropy (ApEn) and spikiness, which reflect the complexity and stability of signaling interactions in homeostatic systems [10].

Other pituitary hormone systems, including GH, TSH, ACTH and LH exhibit age-, gender-, and BMI-related changes [11][13]. Although PRL levels are generally lower in men than in women, the influence of aging and adiposity is less well investigated, especially in men. The small size of most cohorts evaluated to date, the narrow age and BMI ranges encompassed, and the lack of inclusion of both genders collectively make the possibility of statistical type I or type II errors high in earlier studies. As importantly, because of correlations among age, BMI and gender, multivariate regression is needed for definitive inferences. Nonetheless, multivariate analysis also is unreliable in small cohorts. To overcome these obstacles would require investigation of a large number of healthy adults, both men and women, over wide ranges of age and BMI. In this light, the present study examines the dependencies of PRL release (mean, peak, nadir, 24-h secretion, ApEn, spikiness) on individual and/or combined clinical characteristics in 74 healthy individuals sampled frequently (every 10 min) for a sufficiently representative duration (24 h) and analyzed with a high-sensitivity PRL assay (immunofluorometric platform).

Methods

Clinical protocol

The cohort of healthy individuals studied in this project, originated from different studies, in which they served as controls, including studies on PRL secretion in patients with prolactinoma, obese subjects and patients with neurological disorders [3], [14][18]. In these studies the women and men volunteered for and completed the sampling study. Subjects originated from the same community, and were evaluated in an identical sampling paradigm and PRL assay (below). Informed written consent was obtained from the subjects in all these published studies and the studies were approved, including this retrospective analysis, by the ethics committee of the Leiden University Medical Centre. All analyses reported here used techniques not previously applied in any of the published studies. Clinical characteristics of the volunteers (41 women and 33 men) are listed in Table 1. Postmenopausal individuals studied here did not use estrogen therapy. Premenopausal women were studied in the follicular phase of the menstrual cycle.

Participants maintained conventional work and sleeping patterns and reported no recent (within 10 days) transmeridian travel, weight change (>2 kg in 6 weeks), shift work, psychosocial stress, prescription medication use, substance abuse, neuropsychiatric illness, or acute or chronic systemic disease. A complete medical history, physical examination, and screening biochemistry tests were normal. Volunteers were admitted to the Study Unit the evening before sampling for adaptation. Ambulation was permitted to the lavatory only. Vigorous exercise, daytime sleep, snacks, caffeinated beverages, and cigarette smoking were disallowed. Meals were provided at 0900, 1230 and 1730 h, and room lights were turned off between 2200 and 2400 h, depending upon individual sleeping habits. Blood samples (2.0 mL) were withdrawn at 10-min intervals for 24 h. Total blood loss was less than 360 mL. Volunteers were compensated for the time spent in the study.

Assays

Plasma PRL concentrations were measured with a sensitive time-resolved fluoroimmunoassay (Wallac Oy, Turku, Finland). The limit of detection (defined as the value 2 SD above the mean value of the zero standard) was 0.04 µg/L. The assay was calibrated against the 3rd WHO standard 84/500. The intra-assay coefficient of variation (CV) varied from 3.0–5.2% in the assay range 0.1–250 µg/L, with corresponding interassay CV's of 3.4–6.2%. Serum estradiol concentrations were assayed by a sensitive RIA (Spectria Estradiol Sensitive RIA, Orion Diagnostica, Espoo, Finland).The detection limit of the assay is 5 pmol/L. The intraassay CV was 21% at concentrations below 30 pmol, 4.5% at 85 nmol/L and 1.7% at 200 pmol/L. Testosterone was measured by RIA (Siemens Healthcare Diagnostics, Deerfield, Il, USA). The detection limit is 0.2 nmol/L. The intraassay CV at 1 nmol/L is 20% and at 14 nmol/L 12%. Free thyroxine was measured by electrochemoluminescence immunoassay (Elecsys 2010, Roche Diagnostics, Almere, The Netherlands). The detection limit is 0.6 pmol/L and the intraassay CV range amounts 5–8%. Serum IGF-I concentration was measured with the Immulite 2500 system (Diagnostic Products Corporation, Los Angeles, CA, USA). The detection limit is 1.5 nmol/L. The intra-assay variation was 5.0 and 7.5% at levels of 8 and 75 nmol/L, respectively. Simple mean, peak and nadir PRL concentrations Mean (24-h average), peak (single daily maximum) and nadir (single daily minimum) PRL concentrations were determined in each subject.

Deconvolution analysis

Prolactin concentration time series were analyzed via a recently developed automated deconvolution method, empirically validated using hypothalamo-pituitary sampling and simulated pulsatile time series [19][22]. The Matlab-based algorithm first detrends the data and normalizes concentrations to the unit interval [0, 1]. Second, the program creates multiple successive potential pulse-time sets, each containing one burst less via a smoothing process (a nonlinear adaptation of the heat-diffusion equation). Third, a maximum-likelihood expectation estimation method computes all secretion and elimination parameters simultaneously conditional on each of the multiple candidate pulse-time sets. Deconvolution parameters comprise basal secretion (β0), two half-lives (α1,α2), secretory-burst mass (η0, η1), random effects on burst mass (σA), measurement error (σε), and a three-parameter flexible Gamma-secretory-burst waveform (β1, β2, β3). The unit-area normalized shape of secretory bursts (plot of rate of secretion over time) was permitted to differ in the day and night, thus constituting a dual-waveform model of secretion. Two change point times were estimated to demarcate onset of the day and onset of the nighttime waveforms within each 24-h pulse train. For PRL, the fast half-life was represented as 3.5 min constituting 37% of the decay amplitude. The slow half-life was estimated as an unknown variable between 20 and 50 min [23], [24]. All candidate pulse-time sets were deconvolved. Statistical model selection was then performed to distinguish among the independently framed fits of the multiple candidate pulse-time sets using the Akaike information criterion. Observed interpulse intervals were described by a two-parameter Weibull process (more general form of a Poisson process, which uncouples the mean from the variance).The parameters (and units) are frequency (number of bursts per total sampling period, lambda of Weibull distribution), regularity of interpulse intervals (unitless gamma of Weibull), slow half-life (minutes), basal and pulsatile secretion rates (concentration units/session), mass secreted per burst (concentration units), and waveform shape (mode, or time delay to maximal secretion after objectively estimated burst onset, in min). A typical example of part of the graphic output of the deconvolution calculations is shown in Fig. 1.

thumbnail
Figure 1. Part of the graphic output of the deconvolution analysis in a male control subject, the upper panel shows the original serum PRL concentrations (µg/L) and the fitted concentration curve (interrupted line).

Asterisks denote pulse onsets, and the rhomboids the time of waveform switch. Time 0 min is 0900 hr. The second panel shows the secretion rate in µg/L.min. The third panel represents the secretion rate within bursts (normalized secretion over time) for the daytime (interrupted line) and nighttime (continuous line), stressing the difference in time at which the maximal secretion rate is reached. The lowest plot shows the statistical distribution of the interpulse delays.

https://doi.org/10.1371/journal.pone.0031305.g001

ApEn

Approximate Entropy (ApEn) is a sensitive and specific statistic for discriminating insidious differences in serial dynamics. ApEn is calculated for any time series as a single nonnegative number, with zero denoting perfect orderliness, as for a sine wave, and larger ApEn values corresponding to more apparently irregular dynamics [25]. The ApEn metric evaluates the consistency of recurrent subordinate (non-pulsatile) patterns in the data, and thus yields information distinct from and complementary to deconvolution (pulse) analyses [10], [25]. In typical biological applications, ApEn calculations are normalized against the SD of the data series by defining a pattern-reproducibility threshold value of r = 0.2 SD validated for data lengths n≥60 samples [26]. This choice of r limits random effects of low measurement variability (typically ≤0.065 SD), thus allowing discrimination between fine gradations in the orderliness of the underlying process. Validation studies have established the suitability of the m = 1 as the pattern-recurrence length for time series comprising 60<n<300 points, as would be true for many endocrine profiles [26]. For this (m, r) pair, there is quantifiably greater regularity (lower ApEn) of nocturnal GH secretion in adult male than female rats castrated prepubertally, as well as in normal men and women [27]. ApEn is translation- and scale-independent mathematically, which means that adding or multiplying each data value by a fixed number does not alter ApEn [28]. This feature ensures valid comparisons between different mean concentrations, overall variation, or secretion rates due to age, gender, physiological state, and pathology. For example, more irregular (higher ApEn, less orderly) GH secretion occurs in patients with either hypersomatotropism due to GH-secreting pituitary tumors or hyposomatotropism due to hypopituitarism despite 1000-fold differences in GH production [29], [30].

Spikiness

Spikiness was defined as the ratio of the SD of the first-differenced (incremental) time series to the SD of the original series [31]. Spikiness quantifies the extent of sharp, brief, staccato-like unpatterned fluctuations.

Cosinor analysis

The diurnal variation of PRL was analyzed by a non-linear cosine approximation. Measures are the mesor (average level around which the 24-hour oscillation occurs), amplitude (half of the difference between the highest and lowest values) and acrophase (time of the maximum).

Statistical analysis

Stepwise forward-selection multivariate linear regression analysis of untransformed PRL measures was used to examine correlations between preselected PRL parameters (individual dependent variables) and one or more of age, BMI, gender and serum hormone concentrations (independent variables). Statistical comparisons by gender were carried out via the nonparametric Kolmogorov-Smirnov test. Significant contrasts were confirmed by unpaired two-tailed Student's t-test of log-transformed PRL-deconvolution measures. In addition, ANOVA was used for comparisons of more than 2 groups. Data are given as median and absolute range, and as mean and standard deviation or standard error. Analyses used Systat, version 11 (SPSS Inc., Chicago, IL, USA). P<0.05 was considered significant.

Results

Table 1 shows the median (absolute range) subject characteristics for the 33 men and 41 women. Their ages were similar, but the BMI was higher in women than men. Estradiol and testosterone concentrations were as expected for the genders, but free thyroxine was higher in men than women and within normal limits, while IGF-I concentrations did not differ between genders. Mean 24-h, peak and nadir PRL concentrations were all larger in women than in men. ApEn (regularity) and spikiness (brief sharp elevations) were similar in men and women (Table 1).

The serum PRL profiles across the 24 hour cycle are displayed in Fig. 2, showing the marked gender differences, especially during the phase with lights off.

thumbnail
Figure 2. Twenty-four hour serum PRL concentration profiles in 41 healthy women and 33 healthy men.

Blood samples were drawn every 10 min. Data are shown as mean, and the bars represent the SEM. The period with lights off (shown as the black horizontal bar) was between 2300 and 0700 h.

https://doi.org/10.1371/journal.pone.0031305.g002

Stepwise forward-selection multivariate regression analysis was employed to assess the association of individual PRL measures with gender, age and BMI. Mean 24-h PRL concentration was associated with gender (female>male, P<0.0001) and BMI (P = 0.01), but not with age (R2 = 0.308, ANOVA P<0.0001) (see Fig. 3). Nadir PRL concentration correlated with gender only (R2 = 0.077, P = 0.017). However, peak PRL concentration correlated with gender (P<0.0001) and negatively with age (P<0.0001), overall R2 = 0.325, P<0.0001 (see Fig. 4). PRL ApEn, a measure of secretory regularity, and spikiness, a metric of brief, staccato-like increases in secretion, were gender-, age- and BMI-invariant.

thumbnail
Figure 3. Multiple linear regression between age, gender and mean serum PRL concentration.

Data were obtained in 74 healthy subjects, who underwent 24-h blood sampling at 10-min intervals. Male subjects are shown as triangles, female subjects as circles. Data points above the regression plane are open, below they are closed.

https://doi.org/10.1371/journal.pone.0031305.g003

thumbnail
Figure 4. Multiple linear regression between age, gender and maximal PRL concentration in the 24-h serum profile.

Data were obtained in 74 healthy subjects, who underwent 24-h blood sampling at 10-min intervals. Male subjects are shown as triangles, female subjects as circles. Data points above the regression plane are open, below they are closed.

https://doi.org/10.1371/journal.pone.0031305.g004

Based on deconvolution analysis and on unpaired statistical comparisons, gender determined pulsatile PRL secretion (P<0.0001), total secretion (P<0.0001), but not basal secretion. Pulsatile secretion was amplified by 1.5 fold due to increased burst mass in women, with unchanged pulse frequency (Table 2). Forward-selection multivariate regression of PRL deconvolution results demonstrated that basal (nonpulsatile) secretion tended to be associated with BMI (R2 = 0.058, P = 0.03), pulsatile secretion with gender (R2 = 0.152, P = 0.003), and total secretion with gender and BMI (R2 = 0.204, P<0.0001, Fig. 5). Pulse mass was associated with gender (P = 0.001) and with a negative tendency to age (P = 0.038) (Fig. 6).

thumbnail
Figure 5. Multiple linear regression between age, gender and total logarithmically transformed PRL secretion, calculated by deconvolution analysis.

Data were obtained in 74 healthy subjects, who underwent 24-h blood sampling at 10-min intervals. Male subjects are shown as triangles, female subjects as circles. Data points above the regression plane are open, below they are closed.

https://doi.org/10.1371/journal.pone.0031305.g005

thumbnail
Figure 6. Multiple linear regression between age, gender and mean PRL pulse mass, calculated by deconvolution analysis.

Data were obtained in 74 healthy subjects, who underwent 24-h blood sampling at 10-min intervals. Male subjects are shown as triangles, female subjects as circles. Data points above the regression plane are open, below they are closed.

https://doi.org/10.1371/journal.pone.0031305.g006

thumbnail
Table 2. Prolactin deconvolution parameters in men and women.

https://doi.org/10.1371/journal.pone.0031305.t002

We also compared results in premenopausal women with postmenopausal women. In postmenopausal women the following parameters were smaller than in women with a regular menstrual cycle: mean 24 h PRL, peak PRL, fasting PRL, nadir PRL, pulsatile secretion, mean pulse mass and total secretion (see Table 3). However, these parameters were not different in postmenopausal women and men (P-values lying between 0.27 and 0.75) (Table 3). In male subjects older than 50 yr, ApEn (0.942±0.301 vs. 1.258±0.267, P = 0.007) and spikiness (0.363±0.122 vs. 0.463±0.212, P = 0.031) were increased compared with values in younger men. Other parameters were statistically not different.

thumbnail
Table 3. Prolactin secretion characteristics in premenopausal and postmenopausal women, and men.

https://doi.org/10.1371/journal.pone.0031305.t003

In the linear regression analyses of the PRL secretion parameters with age, BMI, gender and serum estradiol concentration, the last was an independent predictor of total PRL secretion (P = 0.002), basal secretion (P = 0.01), mean 24-h concentration (P = 0.002), minimal concentration (P = 0.002), but not of pulsatile secretion (P = 0.51), and maximal concentration (P = 0.32). Since we had testosterone values only in men, and with a few exceptions not in women, the regression analysis of the influence of testosterone was restricted to men. None of the parameters was related to the serum testosterone concentration.

Pulsatile PRL secretion was dependent on the diurnal cycle. During day time PRL secretion in men amounted 51.0 (14–239) µg/liter.24 h and in women 79.0 (26–358) µg/liter.24 h, P = 0.22. During the period with lights off these values were 61.0 (16–111) µg/liter.24 h in men and 111.0 (26–421) µg/liter.24 h in women, P = 0.001 (Fig. 7).

thumbnail
Figure 7. Box plots of the pulsatile PRL secretion during day time and during the period with lights off (2300 -0700 h).

Differences between men and women were calculated with the Kolmogorov-Smirnov test.

https://doi.org/10.1371/journal.pone.0031305.g007

All subjects had a significant diurnal PRL rhythm. The detailed results are displayed in Table 4. The mesor and amplitude were larger in women than men, but the time at which the maximal value was obtained did not differ between genders. The acrophase correlated negatively with BMI and age (R2 = 0.186, P = 0.003), see Fig. 8.

thumbnail
Figure 8. Multiple linear regression between age, BMI and the acrophase of the PRL rhythm.

Data points above the regression plane are open, below the plane they are closed. There were no gender differences in acrophase.

https://doi.org/10.1371/journal.pone.0031305.g008

thumbnail
Table 4. Cosinor analysis of the 24-h serum prolactin profiles in healthy subjects.

https://doi.org/10.1371/journal.pone.0031305.t004

Discussion

This study evaluated simple and complex measures of prolactin secretion in relation to gender (women vs men), age and BMI. Simple measures of prolactin secretion, e.g. nadir, peak and mean levels, as used in clinical settings, showed a strong relation to gender, age (decreased peak value only) and BMI (increased mean and nadir). Deconvolution analysis showed that pulsatile, but not basal (non-pulsatile), secretion was larger in women than men and caused by amplified mean burst mass with unchanged pulse frequency. Basal secretion was positively dependent on BMI and serum estradiol concentration.

Recent publications on 24-h secretion profiles of pituitary hormones in the human, with a 10-min sampling scheme, measured with robust and sensitive assays, and analyzed with operator-independent tools have demonstrated different impacts of gender, age and BMI in multivariate regression analysis. For instance, GH secretion in 100 adult subjects was determined independently by age, BMI and gender, thus decreasing by advancing age and increasing BMI, and with a higher secretion rate in women than men [11]. In contrast, ACTH secretion is larger in men than women, and positively correlated with BMI [13]. On the other hand, TSH secretion is gender invariant and only age-dependent in women [12].

The present study shows that PRL secretion is especially determined by gender and BMI. When analyzed separately postmenopausal women exhibited a 40% decrease in PRL secretion compared with premenopausal women in the follicular phase of their menstrual cycle. Healthy men did not show a decrease after age fifty, but rather exhibited a slight 18% increase, suggesting that age per se may not be the cause of the decrease in women. Therefore, studies on PRL physiology and pathophysiology should include a carefully balanced control group for meaningful comparisons.

One obvious cause for the decrease in PRL in women is the lack of estrogens in the menopausal phase, as suggested by the regression analysis in women. However the role of estrogens in regulating PRL secretion is certainly not straightforward. During the menstrual cycle no influence of different estradiol concentrations on PRL is demonstrable, but in this study only part of the 24 h cycle was analyzed [32]. Furthermore, the PRL increase after TRH injection is independent of the stage of the menstrual cycle [33]. However the value of such studies is limited because of the non-physiological TRH dose, which achieved PRL levels higher than those in spontaneous PRL secretory patterns. On the other hand, high dose orally administered estradiol in postmenopausal women, leading to greatly increased estradiol levels increases PRL. This does not occur after transdermal administration [34]. Other studies have even shown a decrease in PRL after transdermal estrogen administration in postmenopausal women [35], or no change [36]. However, it is conceivable that long-term effects of estrogens regulate the prolactin cell mass, thus explaining hypothetically the remarkable decrease in serum PRL concentration after menopause and possibly also the moderate increase in elderly men by enhanced bioavailable estradiol [37].

Interestingly, estradiol correlated with basal PRL secretion and the minimal 24-h concentration, and secondary to this the mean 24-h concentration and total secretion, but not with pulsatile secretion and the maximal value. This observation suggests a diminished restraint on secretion, but the mechanisms involved are currently not known. In addition, this observation also underlines a marked difference with GH secretion, where estrogens increase specifically pulsatile secretion [38], [39].The present result, however, should be confirmed by estradiol measurements with the more precise high performance liquid chromatography-tandem mass spectrometry [40] and free estradiol. Collectively, the impact of estrogens on PRL secretion is still not completely clear, partly because of nonphysiological experiments and partly because of correlation studies that do not prove a causal role. Potentially, blocking of the estrogen receptor by specific drugs may provide more insight into the physiological role of estrogens on PRL secretion.

Leptin is one of the various factors modulating PRL secretion and the administration restores lactation in the leptin-deficient ob/ob mice [41]. Infusion of leptin raises serum PRL concentration in fasted rats to levels present in normal fed animals [42]. Furthermore, a direct effect on PRL secretion by leptin on the pituitary in vitro has also been demonstrated [43]. Therefore, it is possible that age- and gender-dependent differences in serum leptin concentrations, i.e higher levels in women than men, and a 20% decrease after menopause, modulate PRL secretion as we describe in this study [44]. Unfortunately, leptin levels were measured only in a limited number of subjects, so that we could not investigate the relationship between this hormone and PRL secretion parameters.

Differences in PRL secretion between genders have rarely been investigated in relation to the diurnal cycle and not yet with deconvolution techniques. In this study the gender difference was only present during the sleep period but not during day time. This finding suggests that the dopamine restraint and putative stimulatory factors exhibit circadian properties, but also reveals the impact of gender.

Increased age is characterized by advance (earlier in day) shifting of acrophases of rhythms and diminished amplitude, possibly related to sleep fragmentation and earlier sleep onset and awakening [45]. In this study, age was weakly negatively correlated with PRL pulse mass and the maximal value attained during the 24-h rhythm. Studies which compared nocturnal PRL secretion in elderly and young subjects established decreased pulse amplitude in older age [46], [47]. On the other hand, PRL secretion during the daytime was age-invariant [46].The latter observation suggests that the decrease in nocturnal PRL secretion is not the result of decreased lactotrope cell mass. In another study in men aged 30–96 yr, non-fasting early afternoon serum PRL concentration increased slightly with advancing age, which was attributable to subjects older than 75 yr, whose data are not included in other investigations with a lesser age span [48]. Indeed, another analysis comparing basal and TRH-stimulated PRL concentrations in two age groups (23–45 vs. 56–75 yr) found comparable basal levels, but a negative correlation between age and stimulated PRL concentration [49]. Collectively, these observations are consistent with altered (increased) dopaminergic tonus in the elderly during the sleep period, and which can be restored by metoclopramide administration and diminished responsiveness of the prolactin-secreting cell to TRH [5]. Body composition is an important modulator of hormone secretion, as found for instance for insulin, leptin, TSH, ACTH and cortisol secretion. In this study, total PRL secretion was correlated with BMI in both genders. This result corroborates the notion of enhanced PRL secretion in obesity in most studies [3], [50], but not all [51]. Weight reduction by very low calorie diet or bariatric surgery diminished PRL secretion [52][54]. However, in a study by Ernst and colleagues no change in basal PRL levels was observed after massive weight loss (average 50 kg) in their patients [50]. Whether this finding is the result of differences in experimental design, i.e, frequent blood sampling with deconvolution versus a single fasting morning specimen, is not known.

Approximate Entropy of PRL was increased in elderly men compared with subjects younger than 50 yr, but not in women. A common denominator of irregularity is attenuation of negative feedback compared with feedforward. Thus, pathophysiologies that impair feedback elevate ApEn (process randomness). Well-established conditions of high ApEn include primary failure of a target gland like the thyroid, testis, and ovary, autonomous endocrine tumors, and PCOS [55][61]. Additionally, excessive feedforward enforces irregular patterns. Parathyroid hormone secretion is more disorderly in hyperparathyroidism [62], as is aldosterone secretion in primary and secondary hyperaldosteronism [63]. In keeping with these observations, irregularity can be induced experimentally by muting negative feedback by testosterone, cortisol, and IGF-I, which normally maintain regularity of their upstream hormones (LH, ACTH, GH), as well as by augmenting feedforward by GnRH or GHRH on downstream hormones (LH and GH) [64]. In old age hormone secretory regularity is generally diminished, as described for GH, cortisol, and LH, but not ACTH and TSH [11][13], and here in elderly men for PRL, but not in women. Nonetheless, one previous study in 10 young and 10 old men found a non-significant 20% rise in PRL ApEn in elderly men [65]. Clearly, more elderly healthy subjects, with ages extending into the nineties are required to confirm ApEn changes of PRL secretion with age.

Spikiness is a measure of sharp staccato-like excursions in serial measurements, putatively reflecting acute stimuli from outside normative feedback loops [31]. In elderly men, spikiness was increased compared with younger males. In another large study by us, GH spikiness was determined jointly by gender (higher in women) and BMI (positively), accounting for 29% of inter individual variability. Spikiness also marks other endocrine physiology and pathophysiology, such as greater survival when applied to serial glucose data in patients with protracted critical illness, older age in women when applied to ACTH time series, sharper TSH excursions in hypothyroidism and thyrotropinoma compared with normal, and unstable premenstrual mood-rating dynamics in dysphoria and response to therapy in the latter group [13], [66][68]. Recognition of increased spikiness in older adults is important in planning investigations in this age group because greater spikiness would decrease statistical power by accentuating single-sample variability in aging individuals.

As expected, and concordant with the literature, the mesor and amplitude were larger in women than in men [69]. No gender difference was found for the acrophase of the PRL rhythm. A novel finding in this study was the dependence of the acrophase on BMI and age, both causing an advance shift, thus leading to an earlier maximum. This result is comparable with the advanced shift of cortisol in elderly subjects reported in other studies [70][72]. The mechanism(s) behind these shifts are not known, although changes in the sleep quality have been advocated [70]. However, no detailed studies on the effect of adiposity on hormone acrophases are available.

This study has several limitations. First, the retrospective design could have introduced a selection bias, because the volunteers participated in several independent studies as mentioned in the Introduction. However all subjects were recruited by advertisements from the local area. In addition all subjects underwent a structured medical history assessment, physical examination and routine clinical chemistry. The sampling protocol and chemical methods were identical in these studies and all studies were uneventful. Although we cannot completely rule out stress or anxiety in the volunteers, we did not detect this in the individual serum hormone patterns. Second, as mentioned above we did not use the more sensitive and specific high performance liquid chromatography-tandem mass spectrometry, but a reasonably sensitive, but necessarily less accurate estradiol RIA. Third, we have not included healthy volunteers aged over eighty years. Future inclusion of this cohort into the analysis may strengthen conclusions.

In summary, in healthy adults, selective combinations of gender, age, and BMI specify distinct PRL dynamics, thus requiring balanced representation of these variables in comparative PRL studies.

Author Contributions

Conceived and designed the experiments: FR HP DMK JDV. Performed the experiments: FR HP. Analyzed the data: FR HP DMK JDV. Contributed reagents/materials/analysis tools: DMK. Wrote the paper: FR JDV.

References

  1. 1. Ben-Jonathan N, LaPensee CR, LaPensee EW (2008) What can we learn from rodents about prolactin in humans? Endocr Rev 29: 1–41. er.2007-0017 [pii];10.1210/er.2007-0017 [doi].
  2. 2. Molitch ME (1995) Prolactin. In: Melmed S, editor. The Pituitary. Cambridge MA: Blackwell Science.
  3. 3. Kok P, Roelfsema F, Frolich M, Meinders AE, Pijl H (2004) Prolactin release is enhanced in proportion to excess visceral fat in obese women. J Clin Endocrinol Metab 89: 4445–4449. 10.1210/jc.2003-032184 [doi];89/9/4445 [pii].
  4. 4. Veldhuis JD, Iranmanesh A, Wilkowski MJ, Samojlik E (1994) Neuroendocrine alterations in the somatotropic and lactotropic axes in uremic men. Eur J Endocrinol 131: 489–498.
  5. 5. Greenspan SL, Klibanski A, Rowe JW, Elahi D (1990) Age alters pulsatile prolactin release: influence of dopaminergic inhibition. Am J Physiol 258: E799–E804.
  6. 6. Christensen SE, Jorgensen O, Moller J, Moller N, Orskov H (1985) Body temperature elevation, exercise and serum prolactin concentrations. Acta Endocrinol (Copenh) 109: 458–462.
  7. 7. Corenblum B, Taylor PJ (1981) Mechanisms of control of prolactin release in response to apprehension stress and anesthesia-surgery stress. Fertil Steril 36: 712–715.
  8. 8. Quigley ME, Ropert JF, Yen SS (1981) Acute prolactin release triggered by feeding. J Clin Endocrinol Metab 52: 1043–1045.
  9. 9. Franchimont P, Dourcy C, Legros JJ, Reuter A, Vrindts-Gevaert Y, et al. (1976) Prolactin levels during the menstrual cycle. Clin Endocrinol (Oxf) 5: 643–650.
  10. 10. Veldhuis JD, Keenan DM, Pincus SM (2008) Motivations and methods for analyzing pulsatile hormone secretion. Endocr Rev 29: 823–864. er.2008-0005 [pii];10.1210/er.2008-0005 [doi].
  11. 11. Veldhuis JD, Roelfsema F, Keenan DM, Pincus S (2011) Gender, age, body mass index, and IGF-I individually and jointly determine distinct GH dynamics: analyses in one hundred healthy adults. J Clin Endocrinol Metab 96: 115–121. jc.2010-1669 [pii];10.1210/jc.2010-1669 [doi].
  12. 12. Roelfsema F, Pereira AM, Veldhuis JD, Adriaanse R, Endert E, et al. (2009) Thyrotropin secretion profiles are not different in men and women. J Clin Endocrinol Metab 94: 3964–3967. jc.2009-1155 [pii];10.1210/jc.2009-1155 [doi].
  13. 13. Veldhuis JD, Roelfsema F, Iranmanesh A, Carroll BJ, Keenan DM, et al. (2009) Basal, pulsatile, entropic (patterned), and spiky (staccato-like) properties of ACTH secretion: impact of age, gender, and body mass index. J Clin Endocrinol Metab 94: 4045–4052. jc.2009-1143 [pii];10.1210/jc.2009-1143 [doi].
  14. 14. Donjacour CE, Aziz NA, Roelfsema F, Frolich M, Overeem S, et al. (2011) Effect of sodium oxybate on growth hormone secretion in narcolepsy patients and healthy controls. Am J Physiol Endocrinol Metab 300: E1069–E1075. ajpendo.00623.2010 [pii];10.1152/ajpendo.00623.2010 [doi].
  15. 15. Aziz NA, Pijl H, Frolich M, Roelfsema F, Roos RA (2010) Altered thyrotropic and lactotropic axes regulation in Huntington's disease. Clin Endocrinol (Oxf) 73: 540–545. CEN3836 [pii];10.1111/j.1365-2265.2010.03836.x [doi].
  16. 16. van der Klaauw AA, Pereira AM, van Thiel SW, Frolich M, Iranmanesh A, et al. (2007) Attenuated pulse size, disorderly growth hormone and prolactin secretion with preserved nyctohemeral rhythm distinguish irradiated from surgically treated acromegaly patients. Clin Endocrinol (Oxf) 66: 489–498. CEN2757 [pii];10.1111/j.1365-2265.2006.02757.x [doi].
  17. 17. Kok SW, Meinders AE, Overeem S, Lammers GJ, Roelfsema F, et al. (2002) Reduction of plasma leptin levels and loss of its circadian rhythmicity in hypocretin (orexin)-deficient narcoleptic humans. J Clin Endocrinol Metab 87: 805–809.
  18. 18. Veldman RG, Frolich M, Pincus SM, Veldhuis JD, Roelfsema F (2001) Basal, pulsatile, entropic, and 24-hour rhythmic features of secondary hyperprolactinemia due to functional pituitary stalk disconnection mimic tumoral (primary) hyperprolactinemia. J Clin Endocrinol Metab 86: 1562–1567.
  19. 19. Liu PY, Keenan DM, Kok P, Padmanabhan V, O'Byrne KT, et al. (2009) Sensitivity and specificity of pulse detection using a new deconvolution method. Am J Physiol Endocrinol Metab 297: E538–E544. 00071.2009 [pii];10.1152/ajpendo.00071.2009 [doi].
  20. 20. Keenan DM, Roelfsema F, Biermasz N, Veldhuis JD (2003) Physiological control of pituitary hormone secretory-burst mass, frequency, and waveform: a statistical formulation and analysis. Am J Physiol Regul Integr Comp Physiol 285: R664–R673. 10.1152/ajpregu.00195.2003 [doi];00195.2003 [pii].
  21. 21. Keenan DM, Chattopadhyay S, Veldhuis JD (2005) Composite model of time-varying appearance and disappearance of neurohormone pulse signals in blood. J Theor Biol 236: 242–255. S0022-5193(05)00117-7 [pii];10.1016/j.jtbi.2005.03.008 [doi].
  22. 22. Chattopadhyay SKDM, Veldhuis JD (2008) Probabilistic recovery of neuroendocrine pulsatile secretory and kinetic structure: an alternating discrete and continuous scheme. Quarterly Appl Math 66: 401–421.
  23. 23. Sievertsen GD, Lim VS, Nakawatase C, Frohman LA (1980) Metabolic clearance and secretion rates of human prolactin in normal subjects and in patients with chronic renal failure. J Clin Endocrinol Metab 50: 846–852.
  24. 24. Cooper DS, Ridgway EC, Kliman B, Kjellberg RN, Maloof F (1979) Metabolic clearance and production rates of prolactin in man. J Clin Invest 64: 1669–1680. 10.1172/JCI109629 [doi].
  25. 25. Pincus SM (1991) Approximate entropy as a measure of system complexity. Proc Natl Acad Sci U S A 88: 2297–2301.
  26. 26. Pincus SM, Hartman ML, Roelfsema F, Thorner MO, Veldhuis JD (1999) Hormone pulsatility discrimination via coarse and short time sampling. Am J Physiol 277: E948–E957.
  27. 27. Pincus SM, Gevers EF, Robinson IC, van den Berg G, Roelfsema F, et al. (1996) Females secrete growth hormone with more process irregularity than males in both humans and rats. Am J Physiol 270: E107–E115.
  28. 28. Pincus SM (2000) Irregularity and asynchrony in biologic network signals. Methods Enzymol 321: 149–182. S0076-6879(00)21192-0 [pii].
  29. 29. Roelfsema F, Biermasz NR, Veldhuis JD (2002) Pulsatile, nyctohemeral and entropic characteristics of GH secretion in adult GH-deficient patients: selectively decreased pulsatile release and increased secretory disorderliness with preservation of diurnal timing and gender distinctions. Clin Endocrinol (Oxf) 56: 79–87. 1433 [pii].
  30. 30. Hartman ML, Pincus SM, Johnson ML, Matthews DH, Faunt LM, et al. (1994) Enhanced basal and disorderly growth hormone secretion distinguish acromegalic from normal pulsatile growth hormone release. J Clin Invest 94: 1277–1288. 10.1172/JCI117446 [doi].
  31. 31. Pincus SM, Schmidt PJ, Palladino-Negro P, Rubinow DR (2008) Differentiation of women with premenstrual dysphoric disorder, recurrent brief depression, and healthy controls by daily mood rating dynamics. J Psychiatr Res 42: 337–347. S0022-3956(07)00003-9 [pii];10.1016/j.jpsychires.2007.01.001 [doi].
  32. 32. Genazzani AD, Petraglia F, Volpogni C, Forti G, Surico N, et al. (1994) The duration of prolactin secretory bursts from the pituitary is independent from both prolactin and gonadal steroid plasma levels in women and in men. J Endocrinol Invest 17: 83–89.
  33. 33. Sawin CT, Hershman JM, Boyd AE III, Longcope C, Bacharach P (1978) The relationship of changes in serum estradiol and progesterone during the menstrual cycle to the thyrotropin and prolactin responses to thyrotropin-releasing hormone. J Clin Endocrinol Metab 47: 1296–1302.
  34. 34. Stanosz S, Zochowska E, Safranow K, Sieja K, Stanosz M (2009) Influence of modified transdermal hormone replacement therapy on the concentrations of hormones, growth factors, and bone mineral density in women with osteopenia. Metabolism 58: 1–7. S0026-0495(08)00267-9 [pii];10.1016/j.metabol.2008.07.016 [doi].
  35. 35. Castelo-Branco C, Martinez de Osaba MJ, Fortuny A, Iglesias X, Gonzalez-Merlo J (1995) Circulating hormone levels in menopausal women receiving different hormone replacement therapy regimens. A comparison. J Reprod Med 40: 556–560.
  36. 36. Perrone G, Falaschi P, Capri O, Pastore R, Galoppi P, et al. (1994) Hormonal and metabolic effects of transdermal estradiol/progestagen administration in postmenopausal women. Int J Fertil Menopausal Stud 39: 202–207.
  37. 37. Vermeulen A, Kaufman JM, Goemaere S, van Pelt I (2002) Estradiol in elderly men. Aging Male 5: 98–102.
  38. 38. Ho KY, Evans WS, Blizzard RM, Veldhuis JD, Merriam GR, et al. (1987) Effects of sex and age on the 24-hour profile of growth hormone secretion in man: importance of endogenous estradiol concentrations. J Clin Endocrinol Metab 64: 51–58.
  39. 39. Hudson SB, Schroeder DR, Bailey JN, Mielke KL, Erickson D, et al. (2010) Pre- versus postmenopausal age, estradiol, and peptide-secretagogue type determine pulsatile growth hormone secretion in healthy women: studies using submaximal agonist drive and an estrogen clamp. J Clin Endocrinol Metab 95: 353–360. jc.2009-1769 [pii];10.1210/jc.2009-1769 [doi].
  40. 40. Soldin SJ, Soldin OP (2009) Steroid hormone analysis by tandem mass spectrometry. Clin Chem 55: 1061–1066. clinchem.2007.100008 [pii];10.1373/clinchem.2007.100008 [doi].
  41. 41. Chehab FF (1997) The reproductive side of leptin. Nat Med 3: 952–953.
  42. 42. Watanobe H, Suda T, Wikberg JE, Schioth HB (1999) Evidence that physiological levels of circulating leptin exert a stimulatory effect on luteinizing hormone and prolactin surges in rats. Biochem Biophys Res Commun 263: 162–165. 10.1006/bbrc.1999.1331 [doi];S0006-291X(99)91331-1 [pii].
  43. 43. Yu WH, Kimura M, Walczewska A, Karanth S, McCann SM (1997) Role of leptin in hypothalamic-pituitary function. Proc Natl Acad Sci U S A 94: 1023–1028.
  44. 44. Isidori AM, Strollo F, More M, Caprio M, Aversa A, et al. (2000) Leptin and aging: correlation with endocrine changes in male and female healthy adult populations of different body weights. J Clin Endocrinol Metab 85: 1954–1962.
  45. 45. Copinschi G, Van Cauter E (1995) Effects of ageing on modulation of hormonal secretions by sleep and circadian rhythmicity. Horm Res 43: 20–24.
  46. 46. Iranmanesh A, Mulligan T, Veldhuis JD (1999) Mechanisms subserving the physiological nocturnal relative hypoprolactinemia of healthy older men: dual decline in prolactin secretory burst mass and basal release with preservation of pulse duration, frequency, and interpulse interval–a General Clinical Research Center study. J Clin Endocrinol Metab 84: 1083–1090.
  47. 47. Blackman MR, Kowatch MA, Wehmann RE, Harman SM (1986) Basal serum prolactin levels and prolactin responses to constant infusions of thyrotropin releasing hormone in healthy aging men. J Gerontol 41: 699–705.
  48. 48. Arnetz BB, Lahnborg G, Eneroth P (1986) Age-related differences in the pituitary prolactin response to thyrotropin-releasing hormone. Life Sci 39: 135–139.
  49. 49. McTiernan A, Wu L, Chen C, Chlebowski R, Mossavar-Rahmani Y, et al. (2006) Relation of BMI and physical activity to sex hormones in postmenopausal women. Obesity (Silver Spring) 14: 1662–1677. 14/9/1662 [pii];10.1038/oby.2006.191 [doi].
  50. 50. Ernst B, Thurnheer M, Schultes B (2009) Basal serum prolactin levels in obesity–unrelated to parameters of the metabolic syndrome and unchanged after massive weight loss. Obes Surg 19: 1159–1162. 10.1007/s11695-009-9856-0 [doi].
  51. 51. Mingrone G, Manco M, Iaconelli A, Gniuli D, Bracaglia R, et al. (2008) Prolactin and insulin ultradian secretion and adipose tissue lipoprotein lipase expression in severely obese women after bariatric surgery. Obesity (Silver Spring) 16: 1831–1837. oby2008297 [pii];10.1038/oby.2008.297 [doi].
  52. 52. Friedrich N, Rosskopf D, Brabant G, Volzke H, Nauck M, et al. (2010) Associations of anthropometric parameters with serum TSH, prolactin, IGF-I, and testosterone levels: results of the study of health in Pomerania (SHIP). Exp Clin Endocrinol Diabetes 118: 266–273. 10.1055/s-0029-1225616 [doi].
  53. 53. Camastra S, Manco M, Frascerra S, Iaconelli A, Mingrone G, et al. (2009) Daylong pituitary hormones in morbid obesity: effects of bariatric surgery. Int J Obes (Lond) 33: 166–172. ijo2008226 [pii];10.1038/ijo.2008.226 [doi].
  54. 54. Veldhuis JD, Roemmich JN, Richmond EJ, Bowers CY (2006) Somatotropic and gonadotropic axes linkages in infancy, childhood, and the puberty-adult transition. Endocr Rev 27: 101–140. er.2005-0006 [pii];10.1210/er.2005-0006 [doi].
  55. 55. Pincus SM, Mulligan T, Iranmanesh A, Gheorghiu S, Godschalk M, et al. (1996) Older males secrete luteinizing hormone and testosterone more irregularly, and jointly more asynchronously, than younger males. Proc Natl Acad Sci U S A 93: 14100–14105.
  56. 56. van den Berg G, Pincus SM, Veldhuis JD, Frolich M, Roelfsema F (1997) Greater disorderliness of ACTH and cortisol release accompanies pituitary-dependent Cushing's disease. Eur J Endocrinol 136: 394–400.
  57. 57. van den Berg G, Pincus SM, Frolich M, Veldhuis JD, Roelfsema F (1998) Reduced disorderliness of growth hormone release in biochemically inactive acromegaly after pituitary surgery. Eur J Endocrinol 138: 164–169.
  58. 58. van Aken MO, Pereira AM, van Thiel SW, van den Berg G, Frolich M, et al. (2005) Irregular and frequent cortisol secretory episodes with preserved diurnal rhythmicity in primary adrenal Cushing's syndrome. J Clin Endocrinol Metab 90: 1570–1577. jc.2004-1281 [pii];10.1210/jc.2004-1281 [doi].
  59. 59. Pincus SM, Veldhuis JD, Mulligan T, Iranmanesh A, Evans WS (1997) Effects of age on the irregularity of LH and FSH serum concentrations in women and men. Am J Physiol 273: E989–E995.
  60. 60. Veldhuis JD, Pincus SM, Garcia-Rudaz MC, Ropelato MG, Escobar ME, et al. (2001) Disruption of the synchronous secretion of leptin, LH, and ovarian androgens in nonobese adolescents with the polycystic ovarian syndrome. J Clin Endocrinol Metab 86: 3772–3778.
  61. 61. Veldhuis JD, Zwart AD, Iranmanesh A (1997) Neuroendocrine mechanisms by which selective Leydig cell castration unleashes increased pulsatile LH release. Am J Physiol 272: R464–R474.
  62. 62. Schmitt CP, Locken S, Mehls O, Veldhuis JD, Lehnert T, et al. (2003) PTH pulsatility but not calcium sensitivity is restored after total parathyroidectomy with heterotopic autotransplantation. J Am Soc Nephrol 14: 407–414.
  63. 63. Siragy HM, Vieweg WV, Pincus S, Veldhuis JD (1995) Increased disorderliness and amplified basal and pulsatile aldosterone secretion in patients with primary aldosteronism. J Clin Endocrinol Metab 80: 28–33.
  64. 64. Liu PY, Takahashi PY, Roebuck PD, Veldhuis JD (2006) Age or factors associated with aging attenuate testosterone's concentration-dependent enhancement of the regularity of luteinizing hormone secretion in healthy men. J Clin Endocrinol Metab 91: 4077–4084. jc.2005-2811 [pii];10.1210/jc.2005-2811 [doi].
  65. 65. Veldhuis JD, Iranmanesh A, Mulligan T, Pincus SM (1999) Disruption of the young-adult synchrony between luteinizing hormone release and oscillations in follicle-stimulating hormone, prolactin, and nocturnal penile tumescence (NPT) in healthy older men. J Clin Endocrinol Metab 84: 3498–3505.
  66. 66. Roelfsema F, Pereira AM, Adriaanse R, Endert E, Fliers E, et al. (2010) Thyrotropin secretion in mild and severe primary hypothyroidism is distinguished by amplified burst mass and Basal secretion with increased spikiness and approximate entropy. J Clin Endocrinol Metab 95: 928–934. jc.2009-1959 [pii];10.1210/jc.2009-1959 [doi].
  67. 67. Roelfsema F, Pereira AM, Keenan DM, Veldhuis JD, Romijn JA (2008) Thyrotropin secretion by thyrotropinomas is characterized by increased pulse frequency, delayed diurnal rhythm, enhanced basal secretion, spikiness, and disorderliness. J Clin Endocrinol Metab 93: 4052–4057. jc.2008-1145 [pii];10.1210/jc.2008-1145 [doi].
  68. 68. Pincus SM, Alam S, Rubinow DR, Bhuvaneswar CG, Schmidt PJ (2011) Predicting response to leuprolide of women with premenstrual dysphoric disorder by daily mood rating dynamics. J Psychiatr Res 45: 386–394. S0022-3956(10)00218-9 [pii];10.1016/j.jpsychires.2010.07.006 [doi].
  69. 69. Latta F, Leproult R, Tasali E, Hofmann E, L'Hermite-Baleriaux M, et al. (2005) Sex differences in nocturnal growth hormone and prolactin secretion in healthy older adults: relationships with sleep EEG variables. Sleep 28: 1519–1524.
  70. 70. Sherman B, Wysham C, Pfohl B (1985) Age-related changes in the circadian rhythm of plasma cortisol in man. J Clin Endocrinol Metab 61: 439–443.
  71. 71. Van Cauter E, Leproult R, Kupfer DJ (1996) Effects of gender and age on the levels and circadian rhythmicity of plasma cortisol. J Clin Endocrinol Metab 81: 2468–2473.
  72. 72. Haus E, Nicolau G, Lakatua DJ, Sackett-Lundeen L, Petrescu E (1989) Circadian rhythm parameters of endocrine functions in elderly subjects during the seventh to the ninth decade of life. Chronobiologia 16: 331–352.