Somatostatin is essential for the sexual dimorphism of GH secretion, corticosteroid-binding globulin production, and corticosterone levels in mice.

GROWTH

HORMONE-SOMATOSTATIN-GRH

Somatostatin Is Essential for the Sexual Dimorphism of GH Secretion, Corticosteroid-Binding Globulin Production, and Corticosterone Levels in Mice Jessica M. Adams, Veronica Otero-Corchon,† Geoffrey L. Hammond, Johannes D. Veldhuis, Nathan Qi, and Malcolm J. Low Neuroscience Graduate Program (J.M.A.), University of Michigan, Ann Arbor, Michigan, 48109; Department of Molecular and Integrative Physiology (J.M.A., V.O.-C., M.J.L.), University of Michigan Medical School, Ann Arbor, Michigan, 48109; Department of Cellular and Physiological Sciences (G.L.H.), University of British Columbia, Vancouver, British Columbia, Canada, V6T 1Z3; Department of Internal Medicine (J.D.V.), Endocrine Research Unit, Mayo Clinic College of Medicine, Rochester, Minnesota, 55905; and Department of Internal Medicine (N.Q., M.J.L.), Division of Metabolism, Endocrinology and Diabetes, University of Michigan Medical School, Ann Arbor, Michigan, 48109 Distinct male and female patterns of pituitary GH secretion produce sexually differentiated hepatic gene expression profiles, thereby influencing steroid and xenobiotic metabolism. We used a fully automated system to obtain serial nocturnal blood samples every 15 minutes from cannulated wild-type (WT) and somatostatin knockout (Sst-KO) mice to determine the role of SST, the principal inhibitor of GH release, in the generation of sexually dimorphic GH pulsatility. WT males had lower mean and median GH values, less random GH secretory bursts, and longer trough periods between GH pulses than WT females. Each of these parameters was feminized in male Sst-KO mice, whereas female Sst-KO mice had higher GH levels than all other groups, but GH pulsatility was unaffected. We next performed hepatic mRNA profiling with high-density microarrays. Male Sst-KO mice exhibited a globally feminized pattern of GH-dependent mRNA levels, but female Sst-KO mice were largely unaffected. Among the differentially expressed female-predominant genes was Serpina6, which encodes corticosteroid-binding globulin (CBG). Increased CBG was associated with elevated diurnal peak plasma corticosterone in unstressed WT females and both sexes of Sst-KO mice compared with WT males. Sst-KO mice also had exaggerated ACTH and corticosterone responses to acute restraint stress. However, consistent with their lack of phenotypic signs of excess glucocorticoids, cerebrospinal fluid concentrations of free corticosterone in Sst-KO mice were not elevated. In summary, SST is necessary for the prolonged interpulse troughs that define masculinized pituitary GH secretion. SST also contributes to sexual dimorphism of the hypothalamic-pituitary-adrenal axis via GH-dependent regulation of hepatic CBG production. (Endocrinology 156: 1052–1065, 2015)

H is released in a pulsatile manner from somatotrophs in all species studied (1), is positively regulated by GHRH, and negatively regulated by somatostatin (SST) inputs from the hypothalamus. Ghrelin, produced primarily by the stomach, has an additional role in stimulating GH secretion (2, 3) and acts synergistically with GHRH. The specific properties of GH secretion differ between the sexes in rodents and humans, measured by frequent repeated blood sampling from indwelling venous

G

cannulae. In male rats, low basal GH levels are interrupted by large episodic bursts every 3– 4 hours (4); on the other hand, female rats exhibit higher basal GH levels and more frequent but smaller bursts of GH secretion (5). Humans exhibit a similar but less pronounced sex difference, with females secreting more GH due to more pulses (6). Additionally, in humans and rats, GH secretory patterns are more regular in males than females (7). In the only report of sex differences in GH pulsatility in mice, males and

ISSN Print 0013-7227 ISSN Online 1945-7170 Printed in U.S.A. Copyright © 2015 by the Endocrine Society Received May 21, 2014. Accepted December 23, 2014. First Published Online December 31, 2014

Abbreviations: ApEn, approximate entropy; CBC, corticosteroid-binding capacity; CBG, corticosteroid-binding globulin; CSF, cerebrospinal fluid; CV, coefficient of variability; FDR, false discovery rate; HPA, hypothalamic-pituitary-adrenal; KO, knockout; RMANOVA, repeated measures ANOVA; SST, somatostatin; WT, wild type.

1052

endo.endojournals.org

Endocrinology, March 2015, 156(3):1052–1065

doi: 10.1210/en.2014-1429

The Endocrine Society. Downloaded from press.endocrine.org by [${individualUser.displayName}] on 03 June 2015. at 11:47 For personal use only. No other uses without permission. . All rights reserved.

doi: 10.1210/en.2014-1429

females exhibited similar baseline and peak GH levels, but males had longer interpulse troughs than females (8). An important consequence of the sex difference in GH secretion is the resulting sexual dimorphism of the liver. First shown in rats (9) and later confirmed in mice (10), specific features of ultradian GH secretion dictate sex differences in the hepatic mRNA levels of many genes, including cytochrome p450 and sulfotransferase enzymes. Sex-specific hepatic gene expression is also evident in humans, although the magnitude of expression difference is much less than in rodents, and the sex predominance of many orthologous genes is reversed (11). Sexually dimorphic gene expression patterns are the result of a complex interplay of many transcription factors, the most well studied being signal transducer and activator of transcription 5b (Stat5b), which is phosphorylated in response to GH receptor activation (for review, see Ref. 12). Extensive studies in rats have determined that it is the relatively long interpulse interval of very low or undetectable plasma GH levels, rather than any change in pulse amplitude, duration, or frequency, that is critical to produce the typical masculine hepatic gene expression profile (13). The mechanism driving GH secretory dynamics and pulsatility in male vs female mammals is not fully understood. The principal driver of plasma GH pulsatility is likely GHRH, because neutralization with anti-GHRH antibodies or treatment with GHRH receptor-specific antagonists completely ablates GH pulses in rats (14, 15). On the other hand, immune-neutralization of SST does not affect GH peak amplitude, pulse frequency, or average plasma GH levels but elevates the mean trough levels (16, 17), which is critical to sexual differentiation of hepatic gene expression. Indeed, Sst mRNA levels are greater in the hypothalamus of male rats compared with female throughout development and adulthood (18) and Sst expression is regulated by testosterone (19), suggesting that SST could play a role in sexually dimorphic GH secretion. However, confirmation of this role has been challenging, because there are no available SST receptor antagonists and serial blood sampling has not been done previously in strains of Sst-knockout (Sst-KO) mice (20, 21). Therefore, the current study was designed to define the role of SST in sex-specific GH pulsatility and the subsequent downstream effects on hepatic gene expression in C57BL/6 mice. In a previous study of Sst-KO mice, we inferred the alterations in GH secretion from analysis of multiple, randomly timed GH measurements together with hepatic mRNA levels of 5 sexually dimorphic genes (20). Here, we used an automated blood sampling system in freely behaving Sst-KO and control wild-type (WT) mice of both sexes to directly measure plasma GH pulsatility patterns and correlated those results with global


1053

changes in hepatic mRNA levels measured by microarrays and specific changes in corticosteroid-binding globulin (CBG) and circulating levels of corticosterone.

Materials and Methods Animals Sst-KO mice were generated as described previously (20) and backcrossed onto the C57BL/6J (The Jackson Laboratory) genetic background for more than or equal to 10 generations. Sst-KO and WT littermates were generated from the cross of Sst-heterozygote breeders. Experimental female mice were not controlled for estrous cycle. All mice were group housed in ventilated cages on a 12-hour light, 12-hour dark cycle (lights on from 6 AM to 6 PM) in a temperature-controlled environment (22 ⫾ 1°C) with free access to standard rodent chow and water. All animal studies were approved by either the Institutional Animal Care and Use Committee at Oregon Health and Science University or the University Committee on Use and Care of Animals at University of Michigan and followed Public Health Service guidelines.

Carotid artery cannulations Five- to eight-month-old mice were anesthetized with sodium pentobarbital (50 – 60 mg/kg, ip). Under aseptic conditions, the right carotid artery was isolated. A 2-part cannula constructed from microrenathane tubing (0.635 mm outer diameter, stretched) and SILASTIC tubing (0.635 mm outer diameter) was inserted to the level of the aortic arch. The free end was tunneled sc and exteriorized at the back of the neck via a stainless steel tubing connector that was fixed sc during closure of the skin incision. Cannulae were flushed with saline and sealed with heparin-saline (200 U/mL) daily after surgery. Cannulated mice were housed individually and weighed daily. A total of 30 of 33 mice survived the procedure, and of these, cannulae remained patent in 27 (90%; n ⫽ 8 male WT, 5 female WT, 5 male Sst-KO, 6 female Sst-KO).

Automated blood sampling After surgery, mice appeared healthy with normal activity and body weight regain to more than or equal to 90% of their presurgery level. They were acclimatized for 3– 4 days to Ratturn cages as part of the Culex automated blood sampling system (Basi, Inc), and then, 44 serial blood samples were obtained from the carotid cannulae at 15-minute intervals during the 12-hour nocturnal period in complete darkness starting at 7 PM. Because of the dead volume in cannulae and tubing, 55 ␮L of whole blood was withdrawn at each time point, 10 ␮L was diluted with 50 ␮L of heparinized saline into refrigerated glass tubes, and the remaining 45 ␮L blood with 10 ␮L of heparinized saline was infused back into the mouse to maintain blood volume. The next morning, all mice appeared healthy. Collection tubes were centrifuged, and 50 ␮L of clear, diluted plasma was frozen at ⫺20°C before assay. Samples were considered to be adequate if they contained the expected pellet of red blood cells (⬃5 ␮L). Only occasional samples were excluded from mice with patent cannulae.


1054

Adams et al

SST and Sexually Dimorphic GH Secretion

Manual collection of blood and cerebrospinal fluid (CSF) For basal daytime corticosterone analysis, adult mice (n ⫽ 8 per sex, genotype, and time point) were euthanized by decapitation immediately upon removal from their cages at 9 AM or 5 PM and trunk blood was collected into tubes containing 10 ␮L of 75 mg/mL EDTA. Plasma was separated via centrifugation and stored at ⫺20°C before assay. For the initial restraint stress study, adult mice (n ⫽ 6 per sex, genotype, and time point) were either decapitated between 8 and 10 AM (time ⫽ 0) or subjected to restraint stress in 50-mL Falcon tubes for 5, 10, or 20 minutes and then decapitated. An additional group recovered for 40 minutes after the stress. Trunk blood was processed as above. For the CSF study, adult mice were either anesthetized with tribromoethanol immediately upon removal from their cages (nonstressed controls) between 4:30 and 6 PM or subjected to restraint stress for 20 minutes as described above, then anesthetized and placed in a stereotaxic frame. Under a dissecting microscope, the neck muscles were bluntly dissected through a sagital incision, and the dura mater overlying the spinal medulla and cisterna magna was visualized. The dura was punctured with a tapered glass micropipette, and 4- to 10-␮L CSF was collected within 10 minutes as previously described (22). Samples that were red were considered to be contaminated with blood and discarded. Trunk blood was processed as above, and matching plasma and CSF samples were stored at ⫺20°C. For plasma corticosteroid-binding capacity (CBC), adult mice (n ⫽ 7– 8 per sex, genotype, and time point) were euthanized by decapitation immediately upon removal from their cage at 9 AM or 5 PM, and trunk blood was collected and processed as above. A subset of the samples collected at 9 AM (n ⫽ 4 –5 per sex and genotype) were also analyzed for CBG immunoreactivity via Western blotting as described below.


using 3H-corticosterone and 10␮M corticosterone (competitor) for nonspecific binding determination (24). For Western blotting, plasma samples were diluted 1:500 in 100mM Tris (pH 7.5), boiled for 3 minutes, and 10 ␮L of each sample were assayed using antiserum specific for mouse CBG, as described previously (25).

Deconvolution analysis Plasma GH concentration time series (total 10.75 h) were analyzed by a variable-waveform deconvolution method (26). The MatLab (MathWorks, Inc) program first detrended and normalized concentrations to the unit interval [0, 1] (27), then successive potential pulse-time sets were created by an incremental smoothing process, which deleted the least significant nadir one at a time. Maximum-likelihood expectation parameter estimation was used to calculate secretion and elimination rates simultaneously for each candidate pulse-time set. The model specified basal secretion, a slow-phase half-life, secretory-burst mass, random effects on burst mass, procedural/measurement error, and a 3-parameter flexible ␥ probability distribution to embody secretory-burst waveform. The rapid phase half-life was assumed to be 1.0 minutes, and the slow half-life was calculated to be 8.0 minutes. Finally, the Akaike information criterion was applied to distinguish objectively among candidate pulse-time sets (28). Observed interpulse intervals are described by a 2-parameter Weibull renewal process. Units of parameters are burst frequency (number per 10.75 h), regularity of interpulse intervals (unitless ␥ of Weibull), basal and pulsatile secretion rates (ng/mL per 10.75h), mass secreted per burst (ng), and waveform mode (time delay to maximal secretion after burst onset, min).

Microarrays Three mice per genotype and sex were euthanized between 9 AM and 12 PM by cervical dislocation at age 4 months. Liver RNA

Plasma hormone and CBG assays Prediluted plasma samples from the Culex study were assayed in singlet using ELISA kits for rat/mouse GH (Millipore Corp) and corticosterone (Arbor Assays). Sample concentrations were calculated from the standard curves using a 4PLC sigmoidal fit in Prism 6.0 (GraphPad). Inter- and intraassay coefficients of variability (CVs) were 13% and 8%, respectively, for the GH assay. We also assayed plasma samples in 10 replicates from 2 Snell dwarf (dw/dw) mice that lack circulating GH, and the optical density of each sample was below the absorption of the blank controls. For our experimental samples, results less than the lowest point of the standard curve were assigned that value (0.4 ng/mL after correcting for dilution). The intraassay CV for the corticosterone assay was 8%, and the interassay CV was 8% as reported by the manufacturer. Occasional measurements less than the lowest point of the standard curve were assigned that value (7.8 ng/mL). Basal daytime corticosterone measurements were measured by an in-house assay at the Oregon National Primate Research Center hormone assay core as previously described (23). Plasma ACTH was measured by immunoradiometric assay (Nichols Institute). For the restraint stress and CSF collection experiment, corticosterone was measured by a double antibody RIA kit (MP Biomedicals) with both inter- and intraassay CVs of 7% as reported by the manufacturer. Measurements lower than the lowest point on the standard curve were assigned that value (6.3 ng/mL). A CBC assay was performed as described

was extracted using TRIzol reagent (Invitrogen) and its quality determined by Agilent Bioanalyzer. Total RNA samples with RNA integrity numbers between 7.4 and 7.9 were submitted for cDNA/amplified RNA synthesis using 3⬘ IVT Express kits and hybridization to GeneChip Mouse Genome 430 2.0 Arrays (Affymetrix). After quality control checks, expression values were calculated using a robust multiarray average (29), weighted (30), 3 specific contrasts computed (31), and P values adjusted using false discovery rates (FDRs) (32). Analyses were performed using the affy (33), affyPLM, and limma (34) packages of Bioconductor implemented in the R statistical environment. Expression values were entered into Cluster 3.0, and both genes and samples were hierarchically clustered using average linkage and visualized on a heatmap using Java TreeView. The complete microarray dataset has been deposited to Gene Expression Omnibus under the accession number GSE56520.

Semiquantitative real-time PCR Between 9 AM and 3 PM, 4 –10 hours after serial blood sampling, mice were euthanized and samples of liver collected and stored at ⫺80°C. Total RNA was extracted with RNeasy spin columns (QIAGEN), quantified, and quality checked by NanoDrop (Thermo Scientific), and 400 ng were used to synthesize cDNA with random hexamers (GoScript Reverse Transcription System; Promega). Real-time PCR was performed on a StepOne Real-Time PCR System (Applied Biosystems) using SYBR Green


doi: 10.1210/en.2014-1429

A

Male WT

100

GH (ng/mL)


*

80

80

60

*

*

*

40

20

20

0

0

B GH (ng/mL)

:0 21

0

:0 23

0

:0 01

0

:0 03

0

:0 05

0

:00

*

*

*

*

21

:00

23

:00

01

:00

03

:00

05

:00

Female WT 100

80

* *

*

*

* *

*

80

60

* **

* *

* *

60

* *

40

40

20

20

0

0

:00

C

21

:00

23

:00

01

:00

03

:00

05

:00

Male Sst-KO

100

GH (ng/mL)

19

Female WT

100

19

:00

21

:00

23

:00

01

:00

03

:00

05

:00

Male Sst-KO

100

80

80

*

60

* *

60

40

40

20

20

0

19

*

*

*

* * *

0

:00

D

21

:00

23

:00

01

:00

03

:00

05

:00

Female Sst-KO

100

GH (ng/mL)

* *

0

:0 19

19

:00

80

60

60

40

21

:00

23

:00

01

:00

03

:00

05

:00

*

* * *

*

*

* * *

40

*

* ** *

**

0

:0 03

20

0

:0 19

0

0

:0 21

0

:0 23

0

:0 01

0

:0 05

0

19

:00

21

:00

23

:00

Hour

01

:00

03

:00

05

:00

Hour

Figure 1. Serial plasma GH levels in C57BL/6J mice. Raw traces of plasma GH collected every 15 minutes during the nighttime in 2 representative animals of each sex (A and C, male; B and D, female) and genotype (A and B, WT; C and D, Sst-KO). Asterisks indicate the corresponding GH secretory bursts identified by deconvolution analysis.

Master mix (Life Technologies) and the relative quantification method. Expression of all genes was normalized to peptidylprolyl isomerase A (Ppia) as the internal control. Primers were designed using Primer3 to span at least one intron when possible (35) and are listed in Supplemental Table 1. Primers were used at a final concentration of 300nM and had calculated efficiencies of 85%–100% over a 1000-fold change standard curve. Relative values below the lowest point on the standard curve were as-

Total RNA was prepared from 100 mg of liver tissue per animal (n ⫽ 8 per sex/genotype) using TRIzol reagent (Invitrogen). A total of 30-␮g RNA samples were resuspended in 50% formamide, 17% vol/vol formaldehyde and 1⫻ 3(N-morpholino)propansulfonicacid,loaded in a 1.5% denaturing agarose gel, run 5 hours at 80 V, and transferred to Nylon membranes (Duralon-UV; Stratagene) overnight using 20⫻ SSC. A random priming protocol (Invitrogen) was used to label a rat CBG cDNA with 32P-dCTP (36). Hepatic CBG mRNA was detected after overnight hybridization of Northern blottings with the radiolabeled cDNA in a 50% formamide solution at 42°C, followed by washing protocol, as previously described (36, 37). The Northern blottings were then exposed in a Phosphor Imager, and the relative intensity of the 1.2-kb bands was quantified using IPGel Lab software. Sample RNA loading was controlled by hybridization of a 32P-labeled oligonucleotide (5⬘-ACGGTATCTGATCGTCTTCGAA CC-3⬘) to ribosomal 18s RNA on the same membranes as previously described (37).

Statistical analyses

Female Sst-KO

100

80

20

Northern blottings

60

*

40

19

signed that value, and occasional outliers were identified using Grubbs’ test and removed.

Male WT

100

1055

Repeated measures ANOVA (RMANOVA) followed by Tukey’s multiple comparisons test was used for analyses of serial corticosterone and ACTH measurements. Raw plasma GH values were normalized by log transformation before calculation of the means and then analyzed by two-factor ANOVA. Median GH levels and lengths of interpulse troughs were compared by the nonparametric Kruskal-Wallis test followed by Dunn’s post hoc test. Student’s unpaired t test was used for all other comparisons. For all tests, P ⬍ .05 was considered significant. All statistical analyses were performed with Prism 6.0.

Results WT C57BL/6 mice exhibit sexually dimorphic ultradian GH secretion Representative plasma GH profiles for 2 WT mice of each sex are shown in Figure 1, A and B. There were sig-


1056

Adams et al

Table 1.



Descriptive Statistics and Deconvolution Analysis of Serial GH Measurements

Raw GH (ng/mL) measurements Mean Median Range Mean maximum per group (n ⫽ 5– 8) Deconvolution analysis Total GH secretion rate (ng/mL per 10.75 h) Pulsatile GH secretion rate (ng/mL per 10.75 h) Mass of GH secreted/burst (ng/mL) Basal secretion rate (ng/mL per 10.75 h) Number of bursts/10.75 h Mode of secretory bursts (min) ApEn

Male WT

Female WT

Male Sst-KO

Female Sst-KO

n ⫽ 349 2.9 3.6 ⬍0.4 –262 71 ⫾ 80 n⫽8 696 ⫾ 520 562 ⫾ 516 99 ⫾ 108 134 ⫾ 74 6.3 ⫾ 1.7 14.3 ⫾ 2.1 0.84 ⫾ 0.18

n ⫽ 215 5.3 6.0a ⬍0.4 –201 74 ⫾ 76 n⫽5 782 ⫾ 331 500 ⫾ 337 68 ⫾ 42 282 ⫾ 203 6.8 ⫾ 1.9 16.1 ⫾ 3.5 1.11 ⫾ 0.25c

n ⫽ 217 5.0 5.2a ⬍0.4 –54 26 ⫾ 17 n⫽5 591 ⫾ 317 180 ⫾ 75 36 ⫾ 22 412 ⫾ 250c 5.8 ⫾ 2.6 14.7 ⫾ 2.7 1.08 ⫾ 0.15c

n ⫽ 258 7.8 7.3a,b ⬍0.4 –198 82 ⫾ 71 n⫽6 1248 ⫾ 890 879 ⫾ 814 108 ⫾ 89 369 ⫾ 171c 7.5 ⫾ 1.4 13.3 ⫾ 1.5 1.06 ⫾ 0.15c

For mean GH values, data were normalized by log transformation and a two-way ANOVA showed a significant effect of sex (F1,1037 ⫽ 44.14, P ⬍ .0001) and of genotype (F1,1037 ⫽ 36.99, P ⬍ .0001) but no significant interaction; values shown are back transformed. The nonparametric Kruskal-Wallis ANOVA was used to compare medians of the untransformed data. For mean maximum per group, values are mean ⫾ SD. a P ⬍ .001 compared with WT male, Dunn’s post hoc comparison test. b P ⬍ .05 compared with WT female, Dunn’s post hoc comparison test. c P ⬍ .05 compared with WT male, Student’s unpaired t test.

nificant individual variations in peak GH levels but no systematic differences between sexes for the average maximum plasma GH value per mouse or the range of values. However, both the mean and median GH levels were significantly higher in WT female than male mice (Table 1). To further visualize these results, we organized the data for every GH measurement in rank plots (Figure 2A) (38). The dissociation constant for human GH binding to its receptor is 2 ng/mL (39), and this concentration is likely close to the lower physiological limit for cellular effects of GH. A much higher percentage of GH values were less than 2 ng/mL in male WT (41%) than in female WT (13%) mice, indicating that females were more likely to have a physiologically relevant level of plasma GH than males at any given time point. Additionally, male WT mice had significantly more prolonged trough periods with consecutive GH values of less than 2 ng/mL compared with females (Figure 2B). Deconvolution analysis revealed no significant differences between male and female WT mice in total GH secretion, pulsatile GH secretion rate, GH mass per secretory burst, mode of secretory bursts, or number of secretory bursts in the 10.75-hour time frame (Table 1). There was a trend toward increased basal secretion rate in females (Table 1). Indeed, the only sex difference found for WT mice was in the approximate entropy (ApEn) parameter, which is a measure of pulse regularity, and the lower value in male mice is consistent with more reproducible (nonrandom) pulsing in males than females (Table 1). Notably, although the deconvolution algorithm identified similar numbers of secretory bursts with a range of 6 –7 for both sexes, it is obvious from inspection of the raw GH profiles (Figure 1, A and B) that in male mice, the bursts tend to occur frequently in volleys lasting up to 2 hours

separated by relatively long trough periods or valleys, whereas in female mice, they are more sharply defined and isolated with short intervening troughs. Taken together, the physiologically relevant period of complex GH secretory events separated by prolonged troughs in WT male mice appears to be in the range of 3–3.5 hours, consistent with data from other rodent studies. SST deficiency alters GH secretory patterns in a sex-specific manner Representative plasma GH profiles from Sst-KO mice are shown in Figure 1, C and D. Both male and female Sst-KO mice had significantly higher mean and median GH levels than WT mice of the respective sexes (Table 1), although the range of values was similar in all groups. The percentage of GH values of less than 2 ng/mL was significantly reduced in the male Sst-KO mice (11%) compared with male WT mice (41%) but similar to both WT (13%) and Sst-KO (12%) females (Figure 2A). In addition, the distribution of trough lengths was skewed to the right in the male Sst-KO mice and was indistinguishable from that of female WT and Sst-KO mice (Figure 2B). Further parsing of the GH rank plots revealed that Sst-KO mice had higher GH values than WT mice at the observed concentration for 5%, or OC5, which is one approximation of baseline values (Figure 2A, inset) (40). In the lower 10%– 40% ranks of the data distribution, Sst-KO mice of both sexes were similar to female WT but significantly higher than male WT mice, whereas in the upper quartile, Sst-KO females had higher rank values than all the other groups (Figure 2A). Deconvolution analysis of the GH data from the Sst-KO mice revealed no genotype differences in most of


doi: 10.1210/en.2014-1429


1057

female mice using the parameters of a FDR less than 5% and a greater 8 1.5 100 Male WT than 2-fold change (Supplemental Female WT 1.0 Table 2). This is similar to the numMale Sst-KO 80 6 0.5 Female Sst-KO ber of sexual dimorphic hepatic 0.0 0 2 4 6 8 10 60 genes found by other groups using 4 40 alternative microarray chips and sta2 tistical analyses (41– 44). Table 2 20 lists the genes with the highest fold 0 0 change (⬎8) between WT male and 50 60 70 80 90 100 0 10 20 30 40 50 Ranking (%) Ranking (%) female livers, and these trends are mostly in agreement with the previB ously published reports. Two highly 15 Male WT female-predominant sequences were Female WT 12 identified that have not been previMale Sst-KO Female Sst-KO ously reported, likely because they 9 are not represented by a probe on 6 Agilent microarrays. They are transcription factor 24 (Tcf24), about which little 3 has been published, and the unanno0 tated sequence AI132709. 15 30 45 60 75 90 105 120 In general, the loss of SST was asTrough Length (min) sociated with feminized hepatic gene Figure 2. Distributions of individual GH levels and lengths of GH troughs. A, Rank plots of all GH values from all mice grouped by sex and genotype, divided into lower (left) and upper (right) expression profiles in male mice halves to accommodate the difference in scale of the y-axes. The horizontal dotted line in the left while having little effect on female panel indicates the probable physiologically relevant lower limit of GH at 2 ng/mL. The inset mice. The microarray results of the shows the lowest range of GH values and the vertical dashed line indicates the OC5. A Kruskal206 sexually dimorphic genes were Wallis nonparametric ANOVA of the complete datasets revealed a highly significant (P ⬍ .0001) difference between the groups, and Dunn’s multiple comparison post hoc tests showed visualized by a heatmap (Figure 3). significant differences for every pairwise comparison except between WT females and Sst-KO After clustering by gene, a separation males. To dissect these differences further, we performed the same tests on the 3 subsets of data was readily apparent between the shown on the graphs and found a significant difference between groups (P ⬍ .0001) for each subset. In the lowest 10% of the rankings (inset), WT male and females were not different from WT female predominant (group I) each other but Sst-KO mice of both sexes were significantly different from their WT controls (P ⬍ and the WT male predominant .01). In the lower 50% rankings of the data (left panel), WT males were significantly different (group II) groups. The Sst-KO male (P ⬍ .0001) from all other groups, which were not different from each other. In the upper 50% of the rankings (right panel), Sst-KO females had significantly different (P ⬍ .0001) GH values mRNA levels more closely matched from all other groups, which were not different from each other. Seven GH values of more than both groups of female rather than 100 ng/mL are not shown on the graph but were included in the statistical analyses. B, A WT male mice. Upon closer examidistribution histogram of trough lengths, where trough length is defined by consecutive GH nation, 78% of the 73 male-predomvalues of less than 2 ng/mL at 15-minute intervals. Kruskal-Wallis nonparametric ANOVA followed by Dunn’s pairwise post hoc comparisons showed Male WT different from each other inant and 71% of the 133 femalegroup, P ⬍ .05. predominant genes showed a feminization in Sst-KO males, whereas 85% and 87% of these gene the calculated parameters. However, the basal GH secregroups, respectively, were unchanged in Sst-KO females tion rate for both sexes of Sst-KO mice was significantly higher than that of male WT but not different from female (Supplemental Table 2). Of the 28 sexually dimorphic WT mice or each other (Table 1). The ApEn for male genes whose levels were changed in Sst-KO females, only Sst-KO mice was increased to female WT and Sst-KO lev- 10 were shifted towards a more masculinized profile, els, indicating a decrease in pulse reproducibility (Table 1). whereas 18 were “hyperfeminized.” In contrast, all of the 151 genes with altered levels in Sst-KO male mice were feminized and none were “hypermasculinized.” Sexually dimorphic hepatic gene expression To validate the microarray results, we performed quanprofiles are altered by the absence of SST Our microarray analysis identified 195 annotated titative RT-PCR for a subset of the differentially expressed genes/cDNAs and 11 unannotated sequences that were genes using RNA from livers of an independent cohort of differentially expressed between the livers of WT male and animals. We chose genes from each of 5 categories (male

A

OC5

Number of Troughs

GH (ng/mL)

2.0


1058

Adams et al

Table 2.


Highly Sexually Dimorphic Hepatic Genes in Wild-Type C57BL/6J Mice Revealed by Microarray

Gene symbol Female predominant Xist Cyp2b9 Sult3a1 Slc22a26 Fmo3 Cyp2b13 Sult2a2 Hao2 AI132709 Acot3 Cyp3a41 Rad51b Slc22a27 Slc22a29 Tcf24 Abcd2 Cux2 Nipal1 Cyp17a1 Fam19a2 Cyp3a16 Nt5e Male predominant Cyp4a12 Ddx3y Hsd3b5 Eif2s3y Ugt2b38 Scara5 Nat8 Elovl3 Cyp2d9 Kdm5d C6 Cyp7b1


P value (FDR)

Male:female Previous Log2 FC reference

Inactive X-specific transcripts Cytochrome P450, family 2, subfamily b, polypeptide 9 Sulfotransferase family 3A, member 1 Solute carrier family 22 (organic cation transporter), member 26 Flavin containing monooxygenase 3 Cytochrome P450, family 2, subfamily b, polypeptide 13 Sulfotransferase family 2A, dehydroepiandrosterone (DHEA)preferring, member 2 Hydroxyacid oxidase 2 Expressed sequence AI132709 Acyl-CoA thioesterase 3 Cytochrome P450, family 3, subfamily a, polypeptide 41 RAD51 homolog B Solute carrier family 22, member 27 Solute carrier family 22. member 29 Transcription factor 24 ATP-binding cassette, subfamily D (ALD), member 2 Cut-like homeobox 2 NIPA-like domain containing 1 Cytochrome P450, family 17, subfamily a, polypeptide 1 Family with sequence similarity 19, member A2 Cytochrome P450, family 3, subfamily a, polypeptide 16 5⬘ nucleotidase, ecto

6.82E-07 2.29E-09 8.56E-08 2.51E-12 3.26E-08 3.69E-06 2.46E-04

⫺7.85 ⫺7.18 ⫺7.13 ⫺7.10 ⫺6.94 ⫺6.80 ⫺5.93

43, 44 41, 43, 44 41– 44 41– 44 41– 44 41, 43, 44 41– 44

5.86E-09 7.15E-10 1.94E-04 4.56E-04 9.28E-05 5.84E-08 2.09E-06 7.46E-10 3.82E-07 7.35E-06 9.94E-04 0.004 7.17E-07 0.016 1.14E-05

⫺5.92 ⫺5.73 ⫺4.83 ⫺4.78 ⫺4.70 ⫺4.56 ⫺4.19 ⫺4.15 ⫺4.05 ⫺3.69 ⫺3.10 ⫺3.07 ⫺3.05 ⫺3.00 ⫺3.00

41– 44 None 41, 43, 43 41, 43, 44 41, 44 41 44 None 41– 44 41– 44 41– 44 41– 44 43, 44 41– 44 41– 44

Cytochrome P450, family 4, subfamily a, polypeptide 12 DEAD (Asp-Glu-Ala-Asp) box polypeptide 3, Y-linked Hydroxy-␦-5-steroid dehydrogenase, 3 ␤- and steroid ␦-isomerase 5 Eukaryotic translation initiation factor 2, subunit 3, structural gene Y linked UDP glucuronosyltransferase 2 family, polypeptide B38 Scavenger receptor class A, member 5 (putative) N-acetyltransferase8 (GCN5 related, putative) Elongation of very long chain fatty acids (FEN1/Elo2, SUR4/Elo3, yeast)-like 3 Cytochrome P450, family 2, subfamily d, polypeptide 9 Lysine (K)-specific demethylase 5D Complement component 6 Cytochrome P450, family 7, subfamily b, polypeptide 1

2.68E-06 7.37 5.71E-12 6.75 3.24E-04 6.58

41– 44 41– 44 41– 44

9.44E-11 6.18

41, 43, 44

2.68E-06 4.71E-04 3.84E-08 0.003

5.26 4.73 4.67 4.38

41– 44 41– 44 41– 44 41– 44

0.003 5.84E-08 3.48E-06 1.07E-08

3.52 3.33 3.20 3.10

41, 43, 44 41– 44 41– 44 41– 44

Gene name

All genes with FDR ⬍ 0.05 and Log2(fold change) ⬎ I3I are shown and listed in order of fold change from highest to lowest. Genes with multiple probes on the microarrays are listed only once. The indicated references are for previous studies that reported the same pattern of expression.

or female predominant in WT and direction of change in each sex of Sst-KO mice) that exhibited robust differences in the microarrays. Relative mRNA levels of Abcd2, Cyp2b9, Cyp7b1, Cyp2b13, Cyp17a1, and Slco1a1 showed comparable sex and genotype differences and identical directions of change measured by qRT-PCR and by microarray, although the absolute fold changes among groups was sometimes larger by the microarray method (Figure 4). SST deficiency affects the hypothalamic-pituitaryadrenal axis Stress is known to affect the GH axis, and previous publications have reported significantly elevated daytime

corticosterone levels in another strain of Sst-KO compared with WT mice (21, 45), prompting us to also examine the HPA axis of the Sst-KO mice. First, we measured corticosterone once per hour to evaluate stress in the cannulated mice during blood collection. Corticosterone levels of WT mice were within the expected basal range and exhibited the typical circadian pattern of higher levels at the onset of darkness that declined until the end of the dark cycle (Figure 5A) (46, 47). RMANOVA of all 4 groups revealed a significant main effect of time. Although corticosterone levels of the male Sst-KO mice were almost double that of WT mice at several time points, these differences across the 11-hour time frame were not statisti-


doi: 10.1210/en.2014-1429

-3.00 -2.00 -1.00 0.00 1.00 2.00 3.00

Male WT Sst-KO


Female WT Sst-KO

I.

II.

Figure 3. Heatmap of the microarray data visualizing mRNA levels of the 206 sexually dimorphic genes in WT mice compared with Sst-KO mice of both sexes. Relatively high levels are indicated in red and low levels in green. Values were clustered according to gene and sample using the hierarchical clustering method in Cluster3. Femalepredominant genes clustered to branch I and male-predominant genes clustered to branch II, and the 3 biological replicates from each group clustered together.

cally significant (P ⫽ .07, Sidak post hoc test, effect of genotype). However, in a separate cohort of mice, we found elevated plasma corticosterone levels in unstressed male Sst-KO mice at both the daytime nadir and peak

1059

(Figure 5B). In addition, male Sst-KO mice had enlarged adrenal glands (Figure 5C) due to hypertrophy of the zona fasiculata (data not shown). Next, we measured plasma ACTH and corticosterone at several time points during and after acute restraint stress. ACTH was doubled in Sst-KO males compared with WT males after the 20-minute stress (Figure 5D). Sst-KO mice of both sexes had exaggerated elevations in corticosterone during and after stress compared with their WT controls (Figure 5E). Despite elevated basal corticosterone levels, enlarged adrenals, and exaggerated hormonal stress responses, Sst-KO mice do not exhibit behavioral signs of anxiety or depression (V. Otero-Corchon and M.J. Low, unpublished results) (20, 21), and their thymus glands do not involute (29.6 ⫾ 4.4 vs 27.1 ⫾ 4.0 mg [mean ⫾ SD], female WT and Sst-KO, respectively, P ⫽ .15; and 21.1 ⫾ 3.7 vs 27.1 ⫾ 4.0 mg, male WT and Sst-KO, respectively, P ⫽ .20). Most circulating corticosterone is biologically unavailable to tissues, because it is bound to CBG, encoded by the Serpina6 gene. CBG does not cross the blood brain barrier (48), and corticosterone in the brain has been shown to reflect free corticosterone levels in the plasma (49). We therefore estimated circulating free corticosterone levels by measuring corticosterone in the CSF of unstressed and stressed mice. In contrast to the genotype differences found in total plasma corticosterone levels, free corticosterone levels in CSF did not differ between WT and Sst-KO mice (Figure 6A), consistent with the lack of evidence for functional glucocorticoid excess in the mutant animals. We did not initially identify Serpina6 as a sexually dimorphic gene due to the stringent restrictions on statistical analysis of the microarray data, but using quantitative RT-PCR, we found it to be female predominant (Figure 6B), in agreement with previous literature (25, 42– 44). Furthermore, hepatic mRNA levels of Serpina6 were increased in Sst-KO compared with WT males, suggesting a feminized expression pattern in Sst-KO males (Figure 6B). Northern blottings further confirmed femalepredominant Serpina6 mRNA levels in WT mice, feminization in Sst-KO males, and “ultrafeminization” in Sst-KO females (Figure 6C). The same differential expression patterns held true for plasma CBG protein levels (Figure 6D) and CBC (Figure 6E), indicating that the group differences in Serpina6 mRNA and CBG protein levels parallel the observed differences in total corticosterone levels.

Discussion Although serial blood sampling from indwelling cannulae is considered the gold standard method to measure GH pulsatility, it has only been performed in mice by one other


1060

Adams et al



convolution analysis in the present nighttime study are very similar to *** ** 150 150 150 *** ** those reported in the published day*** 100 100 100 time studies using WT male mice of 50 50 50 n.s. comparable age (3, 50, 56, 57), sug0 0 0 gesting that GH pulsatility is not diurnally regulated, at least in male mice. D 200 E F The only other studies that have * ** 200 600 ** *** 500 examined sex differences in GH pul150 150 *** *** ** 400 satility in mice collected serial blood 100 100 300 ** 200 samples during the daytime from ** 50 50 100 outbred CD-1 mice (8, 58). Using a 0 0 0 similar but older deconvolution algorithm to ours, MacLeod et al (8) found more frequent GH peaks and Figure 4. qRT-PCR validation of the microarray results. mRNA levels relative to a control group significantly shorter nadir trough pe(either WT male or female, depending on the sex predominance of the gene) as measured by qRT-PCR are shown by the open bars. The open circles indicate the 3 biological replicates from riods between peaks, resulting in a the corresponding microarray results. A–C, Three “typical” genes were unaffected in Sst-KO shortened total cycle length, in the females but feminized in Sst-KO males. D–F, Three “atypical” genes showed changes in both female mice. There were no sex difsexes of Sst-KO mice. RT-qPCR values shown are the mean ⫹ SEM (n ⫽ 5–11) and were analyzed by unpaired Student’s t test; *, P ⬍ .05; **, P ⬍ .01; ***, P ⬍ .001. Statistical analyses of the ferences in the range of GH values or microarray results are presented in Supplemental Table 2. the width and area of GH release in each pulse. Our data analysis is in group to our knowledge (8), due primarily to the technical agreement that the defining sex difference in GH secretion challenges of the surgery. An alternative method was de- in mice is the prolonged interpulse trough period in males veloped recently involving the collection of 4-␮L blood compared with females and suggests that any possible cirsamples from the tail tip every 10 minutes for 6 hours from cadian GH variation does not affect the sex difference in awake, but gently restrained, mice (50). Intrinsic to this this parameter. However, our analysis also revealed a sigmethod is frequent handling of the mice during bleeding, nificantly lower ApEn for GH secretion in male mice, inwhich is likely to activate the stress response. Although the dicating that GH secretory bursts are less random in male handling stress was claimed to be negligible, reported cor- compared with female mice. This difference is reflected in the ticosterone levels were relatively high at 370 ng/mL (51). pattern of GH secretory events. We found that the GH seThus, our automated method has the advantage of ob- cretory bursts in female C57BL/6 mice usually occurred as taining serial blood samples from individual mice for ex- isolated events lasting 15–30 minutes, whereas the bursts in tended periods of time during the day or night for the male C57BL/6 mice frequently occurred in overlapping or analysis of hormone patterns in the absence of activating tightly related volleys of 2 or more bursts. In contrast, the GH the stress axis. secretory peaks described previously in CD-1 mice were isoHuman males secrete more GH at night than during the lated and equally narrow in both sexes (8, 58). These differday (52), and nocturnal GH pulses are tightly correlated ences might reflect the different mouse strains. with slow-wave sleep (53), leading to the hypothesis that A possible confounding factor of our results is that the rodents might secrete more GH in the daytime when they female mice could have been in different stages of the esare sleeping. Early studies in rats found no difference in trous cycle at the time of sampling. Early studies indicated male GH pulsatility patterns in the night vs the day (4, 5) no difference in GH secretory patterns across estrous stage but increased pulsatility in females at night (54), suggest- in rats (54, 59). However, reports in humans have been ing that diurnal GH patterns may be more related to light contradictory. Some groups report no change in GH pulthan to sleep-wake patterns. In support of this, a more satility across the menstrual cycle (60, 61), whereas others recent study found that light stimulation in the middle of report elevated mean GH in the late follicular phase (62) the night rapidly down-regulated spontaneous GH secre- and increased pulse frequency in the ovulatory phase (63), tion (55). Together, these studies also suggest the possi- both of which positively correlate with plasma estrogen. bility of circadian variation in GH pulsatility in mice, but The effect of estrous stage on GH pulsatility in mice has no previous study has measured nocturnal GH pulsatility. never been studied, but it is possible that some parameters However, the GH secretory parameters determined by de- that we measured could have been influenced by estrous

B

O

t-K

O

al

e

e

Ss

Ss

t-K

W T

al e

m

al

Fe

O

Ss

al e

Fe m

M

al e

Ss

t-K

t-K

W T

al e

W T

Fe m

t-K

Fe m

al e

Ss

Ss

e

al

O

Slco1a1 mRNA Level (% of Male WT)

O

O

t-K

W T

M

Fe m

al e

W T

M al e

t-K

Ss

e

al

m

Fe

M

Fe

Cyp17a1 mRNA Level (% of Female WT)

O

O

Ss

e

al

M

Fe m

al e

t-K

W T

W T

Fe m

t-K

Ss

al

m

al

M

Fe

Cyp2b13 mRNA Level (% of Female WT)

M al e

W T

O

O

t-K

e

Ss

e

Fe m

al e

al e

M

Ss

e

W T

W T

Cyp7b1 mRNA Level (% of Male WT)

200

Cyp2b9 mRNA Level (% of Female WT)

O

t-K

t-K

Ss

al

m

M

al

e

Fe m

al e

W T

W T

O

Abcd2 mRNA Level (% of Female WT)

M al e

C

n.s.

200

***

M al e

n.s.

200

M al e

A


doi: 10.1210/en.2014-1429


1061

The next goal of our study was to identify the role of SST in mediating 200 200 WT WT sex-specific GH secretory patterns. Sst-KO Sst-KO 150 150 Overall mean and median GH values were increased above WT in both 100 100 sexes of Sst-KO mice, and further analysis revealed sex-specific effects 50 50 of Sst deletion in the distribution of 0 0 GH values over different ranges. The major effect of SST deficiency, howHour Hour ever, was to greatly reduce the fraction of individual GH levels less than *** *** B C 2 ng/mL in male Sst-KO mice, conMale WT 6 400 Female WT sequently eliminating the masculin5 Male Sst-KO 300 ized pattern of prolonged GH Female Sst-KO 4 troughs. Theoretically, SST could in3 200 hibit GH secretion through any of 5 2 * known SST receptors (Sst1–Sst5). 100 1 GHRH neurons in the arcuate nu0 0 cleus express Sst2A in a sexually di09:00 17:00 morphic distribution (64), and GH D E administration in Sst2-KO mice Male WT *** 800 500 Female WT failed to decrease the activity of 400 Male Sst-KO * 600 GHRH neurons (65), suggesting that Female Sst-KO 300 *** SST might contribute to sexually di*** 400 ** morphic GH-mediated negative 200 feedback via Sst2 at the level of the 200 100 hypothalamus. The anterior pitu0 0 itary expresses all 5 SST receptors, 0 10 20 30 40 50 60 0 10 20 30 40 50 60 but Sst2 and Sst5 are thought to be Time (min) Time (min) most important in regulating GH reFigure 5. Characterization of the HPA axis in Sst-KO mice. A, Hourly overnight plasma lease from somatotrophs and are corticosterone levels of mice from the Culex study. RMANOVA on all 4 groups showed no difference between groups but a significant effect of time (F10,200 ⫽ 6.12, P ⬍ .0001). more highly expressed in male vs feRMANOVA restricted to males (left panel) revealed a significant effect of time (F10,110 ⫽ 5.01, male mice (66). However, knockout P ⬍ .001) and a near significant effect of genotype (F1,11 ⫽ 5.06, P ⫽ .07), whereas RMANOVA mice for each of the SST receptors restricted to females showed only a significant effect of time (F10,90 ⫽ 2.69, P ⬍ .01). RMANOVA between WT male and WT female mice revealed a significant effect of time (F10,110 ⫽ 4.17, P ⬍ have been generated with no reports .001) and a near significant effect of sex (F1,11 ⫽ 4.57, P ⫽ .06). B, Corticosterone of elevated circulating GH, suggestmeasurements from unstressed mice at the diurnal nadir (9 AM) and peak (5 PM). Data were ing either functional redundancy analyzed by unpaired Student’s t test, *, P ⬍ .05; ***, P ⬍ .001. C, Adrenal weights analyzed by unpaired Student’s t test, ***, P ⬍ .001. D and E, Hormone responses to 20 minutes of acute among the receptors or developmenrestraint stress starting at time t ⫽ 0. D, RMANOVA of the ACTH response to stress revealed a tal compensations in the knockout significant effect of time (F3,60 ⫽ 25.7, P ⬍ .001) and group (F3,20 ⫽ 5.5, P ⬍ .01) but no models. Future experiments are interaction. Tukey’s multiple comparisons test found a significant difference between male WT needed to characterize GH pulsatiland male Sst-KO mice at 20 minutes (*, P ⬍ .05). E, RMANOVA of the corticosterone response to stress found a significant effect of time (F3,60 ⫽ 62.0, P ⬍ .001), group (F3,20 ⫽ 24.1, P ⬍ ity in mice with individual or com.0001) and an interaction (F9,60 ⫽ 7.8, P ⬍ .001). Tukey’s multiple comparisons test identified pound deletions of the SST receptors pairwise differences between Sst-KO and WT males at 5 and 20 minutes after stress and and to determine the relative imporbetween Sst-KO and WT females at 20 and 60 minutes after stress (**, P ⬍ .01; ***, P ⬍ .001). All data shown are the mean ⫾ SEM. tance of SST effects in the hypothalamus and pituitary gland. Although GH pulsatility in Sst-KO mice has not been stage, leading to increased variability within females and thus a lack of power to detect a difference from males. directly studied previously, Luque et al (67) examined sevNevertheless, collapsing all females regardless of estrous eral other aspects of the GH axis using an alternative strain stage, we still found several significant sex differences that of Sst-KO mice, also shown to have elevated random GH levels in both sexes (21, 67). They found elevated plasma are in agreement with previous literature.

A

Female

0

04

:0

0

:0

02

00

:0

0

0

:0

0

22

Adrenal Weight (mg)

20

:0

18

:0

04

:0

0

Corticosterone (ng/mL)

0

0

02

:0

0

:0

00

:0

0

0

22

:0

20


ACTH (pg/mL)


18

:0

0


Male



25

10 5 0

Plasma CBG IR Intensity (AU)

* 0.4 0.2

0.0 Serpina6 18s

50

O t-K

O

Ss e

Ss

al m Fe

M al e

Ss al

Fe

m

t- K

W T

W T

O t-K

O e

Ss e al M

E

250

* ***

150

**

100

***

1000

200

Plasma CBC (nM)

*** 0.6

t- K

W T al e

W T M

Ss

Fe m

M

m

al

al

e

e

Ss

D

*

0.8

100

0

al e

t-K

O t- K

W T al e

Fe m

al e

O

0

**

150

al e

200

15

M

400

20

#

200

al e

600

Restraint Stress Basal

Fe m

CSF corticosterone (ng/mL)

*

M

50 0

***

800 600

***

400 200

O Ss

t-K

al m Fe

M

al

e

e

Ss

e al m Fe

t-K

O

T W

W T M al e

O Ss al

m Fe

M

al

e

e

al

Ss

e

t-K

t-K

O

T W

W T m

e m

al

e Fe

al M

Fe

Ss

t-K Ss

e al m Fe

t-K

O

T W

W T M al e

O

0

M al e

Serpina6 mRNA/ 18s RNA Intensity (AU)

C

B

* ***

800

W T

Plasma corticosterone (ng/mL)

A


Serpina6 mRNA Level (% of Female WT)

Adams et al

Fe

1062

Figure 6. Plasma and CSF corticosterone, liver Serpina6 and plasma CBG levels. A, Plasma corticosterone (left) and CSF corticosterone (right) both basally (gray squares) and after restraint stress (black circles). The dotted line represents the assay sensitivity of 6.3 ng/mL. B, Hepatic Serpina6 mRNA levels as measured by qRT-PCR (open bars) and microarrays (open circles; #, P ⬍ .001, see Supplemental Table 2). C, Northern blot analysis of Serpina6 mRNA normalized to 18s. D, Western blottings of plasma CBG immunoreactivity (IR). E, CBC was measured at 9 AM and 5 PM, but no diurnal rhythm was found and results displayed are collapsed across time points. All data shown, except microarray data in B, are mean ⫾ SEM (n ⫽ 4 –11) and were analyzed by unpaired Student’s t test; *, P ⬍ .05; **, P ⬍ .01; ***, P ⬍ .001.

IGF-I, increased pituitary Gh, Ghrhr, and Ghsr expression, and decreased hypothalamic Ghrh expression only in female Sst-KO mice, leading them to conclude that endogenous SST plays a predominant role in females in the regulation of the GH axis (67). Those findings are consistent with our observations that female Sst-KO mice exhibited the highest levels of mean GH, a skewed distribution of ranked GH levels, and the largest GH secretion rate compared with the other 3 groups. However, our results suggest that a lack of SST more profoundly affects the GH interpulse trough interval in males than it does in females. It is possible that in normal females, estrogen signaling interacts with SST signaling in the hypothalamus and the pituitary to organize the GHRH-GH signaling axis. In males, testosterone directly activates SST neurons, which act as a brake on GH release from somatotrophs. Thus, the sex steroid milieu dictates a differential role for SST and, while in females, a lack of SST affects expression of several key genes in the GH axis, SST is more physiologically relevant in males because of the profound effects of Sst deletion on male GH secretory patterns and GH-regulated hepatic mRNA levels. An intriguing observation from our analysis of ultradian GH secretion and sexually dimorphic hepatic gene expression was the elevation of CBG production in both

sexes of Sst-KO mice and in WT female compared with WT male mice. Although higher Serpina6 mRNA levels have been noted previously in the livers of female mice compared with males (25, 42– 44), and serum CBG levels have been shown to be dependent on GH pulse patterns (68), to our knowledge, a potential mechanistic relationship between sexually dimorphic GH secretion and sexual dimorphism of the HPA axis has never been proposed. Previous reports of elevated corticosterone levels in Sst-KO mice attributed the increase to a putative loss of SST inhibition of ACTH release from pituitary corticotrophs (45). That study demonstrated in vitro effects of SST to inhibit basal ACTH release and to attenuate CRHinduced ACTH release from primary cultured pituitary cells of both rats and nonhuman primates. Our in vivo data showing exaggerated ACTH and corticosterone responses to acute restraint stress in Sst-KO mice and increased adrenal gland weight in male Sst-KO mice are consistent with this possibility. Paradoxically, the mutant mice do not exhibit phenotypic signs of glucocorticoid excess, and free corticosterone levels in CSF are the same in both sexes and genotypes of mice after acute restraint stress, despite their differences in total corticosterone levels in plasma.


doi: 10.1210/en.2014-1429

Therefore, we propose that SST has 2 distinct and parallel actions in regulation of the HPA axis. SST indirectly contributes to the sexual dimorphism of the HPA axis through its shaping of sexually dimorphic GH secretion and hepatic gene expression. Sex differences in the levels of circulating CBG produced by the liver contribute to the higher total corticosterone levels (bound and free) measured in WT female mice. Due to feminized GH pulsatility, Sst-KO male mice have feminized levels of CBG and hence elevated total corticosterone levels but do not exhibit a Cushingoid phenotype, because their biologically available free corticosterone remains normal. In addition, SST can inhibit ACTH release from the pituitary (45). A normally higher pituitary SST tone in male mice, proposed as a mechanism for the sexual dimorphism of GH secretion, may therefore lead to chronically slightly increased basal ACTH levels in male Sst-KO mice, which we have been unable to confirm due to assay limitations (data not shown). This could explain the hypertrophied adrenal zona fasciculata in male Sst-KO mice. However, SST inhibition of ACTH release may be more physiologically relevant in the context of simultaneous CRH stimulation as occurs during stress, as we have shown. In summary, we used an automated blood sampling system to characterize the sex differences in GH pulsatility in WT C57BL/6J mice and to determine the necessity of SST for this sexual dimorphism. Carotid artery cannulation and automated sampling were well tolerated by the mice, as reflected by normal corticosterone levels, and the results accurately portrayed physiological changes of plasma hormones. We found that male mice have lower mean and median GH, more nonrandom GH secretory bursts, and more consecutive nadir GH values of less than 2 ng/mL than female mice during the nocturnal period. In the absence of SST, male mice exhibited a feminized pattern of GH secretion, associated with a global feminization of sexually dimorphic hepatic mRNA levels. SST deficiency in female mice also increased mean and median GH levels but produced no measurable change in GH pulsatility and only affected hepatic mRNA levels in a handful of sexually dimorphic genes, including Serpina6. Together, these results suggest that the default GH secretory pattern is feminine and that SST is necessary to masculinize the hypothalamic-pituitary-liver axis. As a consequence, decreased Serpina6 expression and plasma CBG contribute to the lower circulating total corticosterone levels in WT male compared with female mice.

Acknowledgments We thank Richard Miller for providing us with plasma from dw/dw mice, Craig Jaffe for insightful discussion, Amanda


1063

Sharpe for CSF collections from the cisterna magna, and Courtney Attard and Jeffrey Jian for technical assistance and the University of Michigan Center for Statistical Consultation and Research for advice on statistical analyses. Address all correspondence and requests for reprints to: Malcolm J. Low, MD, PhD, Department of Molecular and Integrative Physiology, University of Michigan Medical School, 6116 Brehm Tower, 1000 Wall Street, Ann Arbor, MI 48105. E-mail: [email protected]. † Deceased, 18 February, 2013. This work was supported by the National Institutes of Health (NIH) Early Stage Training in the Neurosciences Training Grant T32-NS076401 (to J.M.A.). G.L.H. was supported by a Canada Research Chair in Reproductive Health and an operating grant from Canadian Institutes of Health Research (MOP-111102). M.J.L. was supported by a NIH Grant R01-DK066604. This work used core services provided by the University of Michigan Animal Phenotyping Core supported by the Michigan Nutrition and Obesity Research Center and the Michigan Diabetes Research Center (NIH Grants P30-DK089503 and P30DK020572) and the Molecular Biology Core supported by the Michigan Gastrointestinal Peptide Research Center (NIH Grant P30-DK034933). Disclosure Summary: The authors have nothing to disclose.

References 1. Gahete MD, Durán-Prado M, Luque RM, et al. Understanding the multifactorial control of growth hormone release by somatotropes: lessons from comparative endocrinology. Ann NY Acad Sci. 2009; 1163:137–153. 2. Kojima M, Hosoda H, Date Y, Nakazato M, Matsuo H, Kangawa K. Ghrelin is a growth-hormone-releasing acylated peptide from stomach. Nature. 1999;402(6762):656 – 660. 3. Hassouna R, Zizzari P, Tomasetto C, et al. An early reduction in GH peak amplitude in preproghrelin-deficient male mice has a minor impact on linear growth. Endocrinology. 2014;155(9):3561–3571. 4. Tannenbaum GS, Martin JB. Evidence for an endogenous ultradian rhythm governing growth hormone secretion in the rat. Endocrinology. 1976;98(3):562–570. 5. Edén S. Age- and sex-related differences in episodic growth hormone secretion in the rat. Endocrinology. 1979;105(2):555–560. 6. Hartman ML, Veldhuis JD, Vance ML, Faria AC, Furlanetto RW, Thorner MO. Somatotropin pulse frequency and basal concentrations are increased in acromegaly and are reduced by successful therapy. J Clin Endocrinol Metab. 1990;70(5):1375–1384. 7. Pincus SM, Gevers EF, Robinson I, et al. Females secrete growth hormone with more process irregularity than males in both humans and rats. Am J Physiol. 1996;270(1):E107–E115. 8. MacLeod JN, Pampori NA, Shapiro BH. Sex differences in the ultradian pattern of plasma growth hormone concentrations in mice. J Endocrinol. 1991;131(3):395–399. 9. Mode A, Norstedt G, Simic B, Eneroth P, Gustafsson JA. Continuous infusion of growth hormone feminizes hepatic steroid metabolism in the rat. Endocrinology. 1981;108(6):2103–2108. 10. Norstedt G, Palmiter R. Secretory rhythm of growth hormone regulates sexual differentiation of mouse liver. Cell. 1984;36(4):805– 812. 11. Zhang Y, Klein K, Sugathan A, et al. Transcriptional profiling of


1064

12.

13.

14.

15.

16.

17.

18.

19.

20.

21.

22.

23.

24.

25.

26.

27.

28.

Adams et al


human liver identifies sex-biased genes associated with polygenic dyslipidemia and coronary artery disease. Plos One. 2011;6(8). Waxman DJ, Holloway MG. Sex differences in the expression of hepatic drug metabolizing enzymes. Mol Pharmacol. 2009;76(2): 215–228. Waxman DJ, Pampori NA, Ram PA, Agrawal AK, Shapiro BH. Interpulse interval in circulating growth hormone patterns regulates sexually dimorphic expression of hepatic cytochrome p450. Proc Natl Acad Sci USA. 1991;88(15):6868 – 6872. Wehrenberg WB, Brazeau P, Luben R, Böhlen P, Guillemin R. Inhibition of the pulsatile secretion of growth hormone by monoclonal antibodies to the hypothalamic growth hormone releasing factor (GRF). Endocrinology. 1982;111(6):2147–2148. Lumpkin MD, McDonald JK. Blockade of growth hormone-releasing factor (GRF) activity in the pituitary and hypothalamus of the conscious rat with a peptidic GRF antagonist. Endocrinology. 1989; 124(3):1522–1531. Tannenbaum GS, Epelbaum J, Colle E, Brazeau P, Martin JB. Antiserum to somatostatin reverses starvation-induced inhibition of growth hormone but not insulin secretion. Endocrinology. 1978; 102(6):1909 –1914. Terry LC, Martin JB. The effects of lateral hypothalamic-medial forebrain stimulation and somatostatin antiserum on pulsatile growth hormone secretion in freely behaving rats: evidence for a dual regulatory mechanism. Endocrinology. 1981;109(2):622– 627. Argente J, Chowen JA, Zeitler P, Clifton DK, Steiner RA. Sexual dimorphism of growth hormone-releasing hormone and somatostatin gene expression in the hypothalamus of the rat during development. Endocrinology. 1991;128(5):2369 –2375. Chowen-Breed JA, Steiner RA, Clifton DK. Sexual dimorphism and testosterone-dependent regulation of somatostatin gene expression in the periventricular nucleus of the rat brain. Endocrinology. 1989; 125(1):357–362. Low MJ, Otero-Corchon V, Parlow AF, et al. Somatostatin is required for masculinization of growth hormone-regulated hepatic gene expression but not of somatic growth. J Clin Invest. 2001; 107(12):1571–1580. Zeyda T, Diehl N, Paylor R, Brennan MB, Hochgeschwender U. Impairment in motor learning of somatostatin null mutant mice. Brain Res. 2001;906(1–2):107–114. Liu L, Herukka SK, Minkeviciene R, van Groen T, Tanila H. Longitudinal observation on CSF A␤42 levels in young to middle-aged amyloid precursor protein/presenilin-1 doubly transgenic mice. Neurobiol Dis. 2004;17(3):516 –523. Gruenewald DA, Hess DL, Wilkinson CW, Matsumoto AM. Excessive testicular progesterone secretion in aged male Fischer 344 rats: a potential cause of age-related gonadotropin suppression and confounding variable in aging studies. J Gerontol. 1992;47(5): B164 –B170. Hammond GL, Lähteenmäki PL. A versatile method for the determination of serum cortisol binding globulin and sex hormone binding globulin binding capacities. Clin Chim Acta. 1983;132(1):101– 110. Scrocchi LA, Orava M, Smith CL, Han VK, Hammond GL. Spatial and temporal distribution of corticosteroid-binding globulin and its messenger ribonucleic acid in embryonic and fetal mice. Endocrinology. 1993;132(2):903–909. Veldhuis JD, Keenan DM, Pincus SM. Motivations and methods for analyzing pulsatile hormone secretion. Endocr Rev. 2008;29(7): 823– 864. Keenan DM, Iranmanesh A, Veldhuis JD. Analytical construct of reversible desensitization of pituitary-testicular signaling: illustrative application in aging. Am J Physiol Regul Integr Comp Physiol. 2011;300(2):R349 –R360. Akaike H. A new look at statistical model identification. IEEE Trans. 1974;AC19(6):716 –723.


29. Irizarry RA, Hobbs B, Collin F, et al. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003;4(2):249 –264. 30. Ritchie ME, Diyagama D, Neilson J, et al. Empirical array quality weights in the analysis of microarray data. BMC Bioinformatics. 2006;7:261. 31. Smyth GK. Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:Article 3. 32. Benjamini Y, Hochberg Y. Controlling the false discovery rate - A practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol. 1995;57(1):289 –300. 33. Gautier L, Cope L, Bolstad BM, Irizarry RA. affy–analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004; 20(3):307–315. 34. Smyth GK. Limma: linear models for microarray data. In: Gentleman R, Carey V, Huber W, Irizarry R, Dudoit S, eds. Bioinformatics and Computational Biology Solution Using R and Bioconductor. Springer, New York; 2005:397– 420. 35. Rozen S, Skaletsky HJ. Primer3. Code available at http://wwwgenome.wi.mit.edu/genome_software/other/primer3.html. 1998. 36. Smith CL, Hammond GL. Rat corticosteroid binding globulin: primary structure and messenger ribonucleic acid levels in the liver under different physiological conditions. Mol Endocrinol. 1989; 3(2):420 – 426. 37. Smith CL, Hammond GL. Ontogeny of corticosteroid-binding globulin biosynthesis in the rat. Endocrinology. 1991;128(2):983–988. 38. Xu J, Bekaert AJ, Dupont J, et al. Exploring endocrine GH pattern in mice using rank plot analysis and random blood samples. J Endocrinol. 2011;208(2):119 –129. 39. Fuh G, Cunningham BC, Fukunaga R, Nagata S, Goeddel DV, Wells JA. Rational design of potent antagonists to the human growth hormone receptor. Science. 1992;256(5064):1677–1680. 40. Matthews DR, Hindmarsh PC, Pringle PJ, Brook CGD. A distribution method for analyzing the baseline of pulsatile endocrine signals as exemplified by 24-hour growth hormone profiles Clin Endocrinol (Oxf). 1991;35(3):245–252. 41. Clodfelter KH, Holloway MG, Hodor P, Park SH, Ray WJ, Waxman DJ. Sex-dependent liver gene expression is extensive and largely dependent upon signal transducer and activator of transcription 5b (STAT5b): STAT5b-dependent activation of male genes and repression of female genes revealed by microarray analysis. Mol Endocrinol. 2006;20(6):1333–1351. 42. Yang X, Schadt EE, Wang S, et al. Tissue-specific expression and regulation of sexually dimorphic genes in mice. Genome Res. 2006; 16(8):995–1004. 43. Gatti DM, Zhao N, Chesler EJ, et al. Sex-specific gene expression in the BXD mouse liver. Physiol Genomics. 2010;42(3):456 – 468. 44. Conforto TL, Zhang Y, Sherman J, Waxman DJ. Impact of CUX2 on the female mouse liver transcriptome: activation of female-biased genes and repression of male-biased genes. Mol Cell Biol. 2012; 32(22):4611– 4627. 45. Luque RM, Gahete MD, Hochgeschwender U, Kineman RD. Evidence that endogenous SST inhibits ACTH and ghrelin expression by independent pathways. Am J Physiol. 2006;291(2):E395–E403. 46. Saito M, Bray GA. Diurnal rhythm for corticosterone in obese (ob/ ob) diabetes (db/db) and gold-thioglucose-induced obesity in mice. Endocrinology. 1983;113(6):2181–2185. 47. Zaki MS, Katayama T, Murata T, Konishi H, Shiota K, Takahashi M. The regulation of food intake and correlated energy balance in mice. J Vet Med Sci. 1991;53(2):249 –254. 48. Schwarz S, Pohl P. Steroid hormones and steroid hormone binding globulins in cerebrospinal fluid studied in individuals with intact and with disturbed blood-cerebrospinal fluid barrier. Neuroendocrinology. 1992;55(2):174 –182. 49. Qian X, Droste SK, Gutièrrez-Mecinas M, et al. A rapid release of corticosteroid-binding globulin from the liver restrains the gluco-


doi: 10.1210/en.2014-1429

50.

51.

52.

53.

54.

55.

56.

57.

58.

corticoid hormone response to acute stress. Endocrinology. 2011; 152(10):3738 –3748. Steyn FJ, Huang L, Ngo ST, et al. Development of a method for the determination of pulsatile growth hormone secretion in mice. Endocrinology. 2011;152(8):3165–3171. Steyn FJ, Leong JW, Huang L, et al. GH does not modulate the early fasting-induced release of free fatty acids in mice. Endocrinology. 2012;153(1):273–282. Jaffe CA, Ocampo-Lim B, Guo W, et al. Regulatory mechanisms of growth hormone secretion are sexually dimorphic. J Clin Invest. 1998;102(1):153–164. Vancauter E, Kerkhofs M, Caufriez A, Vanonderbergen A, Thorner MO, Copinschi G. A quantitative estimation of growth hormone secretion in normal man - Reproducibility and relation to sleep and time of day J Clin Endocrinol Metab. 1992;74(6):1441–1450. Clark RG, Carlsson LM, Robinson IC. Growth hormone secretory profiles in conscious female rats. J Endocrinol. 1987;114(3):399 – 407. Davies JS, Carter DA, Wells T. Photic stimulation inhibits growth hormone secretion in rats: a hypothalamic mechanism for transient entrainment. Endocrinology. 2004;145(6):2950 –2958. Steyn FJ, Lee K, Fogarty MJ, et al. Growth hormone secretion is correlated with neuromuscular innervation rather than motor neuron number in early-symptomatic male amyotrophic lateral sclerosis mice. Endocrinology. 2013;154(12):4695– 4706. Tan HY, Huang L, Simmons D, Veldhuis JD, Steyn FJ, Chen C. Hypothalamic distribution of somatostatin mRNA expressing neurones relative to pubertal and adult changes in pulsatile growth hormone secretion in mice. J Neuroendocrinol. 2013;25(10):910 – 919. Pampori NA, Shapiro BH. Effects of neonatally administered monosodium glutamate on the sexually dimorphic profiles of circulating growth hormone regulating murine hepatic monooxygenases. Biochem Pharmacol. 1994;47(7):1221–1229.


1065

59. Takahashi S, Gottschall PE, Quigley KL, Goya RG, Meites J. Growth hormone secretory patterns in young, middle-aged and old female rats. Neuroendocrinology. 1987;46(2):137–142. 60. Klein NA, Battaglia DE, Miller PB, Soules MR. Circulating levels of growth hormone, insulin-like growth factor-I and growth hormone binding protein in normal women of advanced reproductive age. Clin Endocrinol (Oxf). 1996;44(3):285–292. 61. Jaffe CA, Ocampo-Lim B, Guo W, et al. Growth hormone secretory dynamics over the menstrual cycle. Endocr J. 2000;47(5):549 –556. 62. Faria ACS, Bekenstein LW, Booth RA, et al. Pulsatile growth hormone release in normal women during the menstrual cycle. Clin Endocrinol (Oxf). 1992;36(6):591–596. 63. Ovesen P, Vahl N, Fisker S, Veldhuis JD, Christiansen JS, Jørgensen JO. Increased pulsatile, but not basal, growth hormone secretion rates and plasma insulin-like growth factor I levels during the periovulatory interval in normal women. J Clin Endocrinol Metab. 1998;83(5):1662–1667. 64. Bouyer K, Loudes C, Robinson IC, Epelbaum J, Faivre-Bauman A. Sexually dimorphic distribution of sst2A somatostatin receptors on growth hormone-releasing hormone neurons in mice. Endocrinology. 2006;147(6):2670 –2674. 65. Zheng H, Bailey A, Jiang MH, et al. Somatostatin receptor subtype 2 knockout mice are refractory to growth hormone-negative feedback on arcuate neurons. Mol Endocrinol. 1997;11(11):1709 – 1717. 66. Cordoba-Chacon J, Gahete MD, Castano JP, Kineman RD, Luque RM. Somatostatin and its receptors contribute in a tissue-specific manner to the sex-dependent metabolic (fed/fasting) control of growth hormone axis in mice. Am J Physiol. 2011;300(1):E46 –E54. 67. Luque RM, Kineman RD. Gender-dependent role of endogenous somatostatin in regulating growth hormone-axis function in mice. Endocrinology. 2007;148(12):5998 – 6006. 68. Jansson JO, Oscarsson J, Mode A, Ritzén EM. Plasma growth hormone pattern and androgens influence the levels of corticosteroidbinding globulin in rat serum. J Endocrinol. 1989;122(3):725–732.


Sexual dimorphism of pyridostigmine potentiation of growth hormone (GH)-releasing hormone-induced GH release in humans.

IGF-1 axis contribute to skeletal sexual dimorphism? Evidence from mouse studies.

Reliability of stimulated and spontaneous growth hormone (GH) levels for identifying the child with low GH secretion.

Sexual dimorphism of somatostatin and growth hormone-releasing factor signaling in the control of pulsatile growth hormone secretion in the rat.

Is Sexual Size Dimorphism Inherent in the Scallop Patinopecten yessoensis?

The role of sexual steroids in the modulation of growth hormone (GH) secretion in humans.

Sexual dimorphism in autoimmunity.

The dilemma of choosing a reference character for measuring sexual size dimorphism, sexual body component dimorphism, and character scaling: cryptic dimorphism and allometry in the scorpion Hadrurus arizonensis.

Impaired inhibitory effects of somatostatin on growth hormone (GH)-releasing hormone stimulation of GH secretion after short term infusion.

Sexual dimorphism of growth hormone-releasing hormone and somatostatin gene expression in the hypothalamus of the rat during development.

Parasitism and the expression of sexual dimorphism.

Masculinization of growth hormone (GH) secretory pattern by dihydrotestosterone is associated with augmentation of hypothalamic somatostatin and GH-releasing hormone mRNA levels in ovariectomized adult rats.

Growth hormone (GH) autofeedback action in the neonatal rat: involvement of GH-releasing hormone and somatostatin.

corticosterone secretion.

Sexual dimorphism in the weight of chicks.

NMDA receptor-deficient mice display sexual dimorphism in the onset and severity of behavioural abnormalities.

Sexual dimorphism in the preferential secretion of unsaturated lysophosphatidylcholine by rat hepatocytes but no secretion by sheep hepatocytes.

Survivin is essential for fertile egg production and female fertility in mice.

Sexual dimorphism in medulloblastoma features.

Human chromosomes. Evidence for autosomal sexual dimorphism.

Sexual dimorphism in the vertebrate nervous system.

Autophagy is essential for hearing in mice.

Embryonic differentiation of sexual dimorphism in vasotocin and mesotocin levels in chickens.

The earliest fossil evidence for sexual dimorphism in primates.