Introduction
Altered volumes of the pituitary gland have been associated with increased vulnerability for psychotic disorder (Garner et al. Reference Garner, Pariante, Wood, Velakoulis, Phillips, Soulsby, Brewer, Smith, Dazzan, Berger, Yung, van den Buuse, Murray, McGorry and Pantelis2005). Patients with a first episode of psychosis (FEP) display increased pituitary volume, compared to controls (Pariante et al. Reference Pariante, Dazzan, Danese, Morgan, Brudaglio, Morgan, Fearon, Orr, Hutchinson, Pantelis, Velakoulis, Jones, Leff and Murray2005; Takahashi et al. Reference Takahashi, Suzuki, Velakoulis, Lorenzetti, Soulsby, Zhou, Nakamura, Seto, Kurachi and Pantelis2009), as do individuals with early prodromal signs and subsequent transition to psychotic disorder, compared to those who do not make the transition (Garner et al. Reference Garner, Pariante, Wood, Velakoulis, Phillips, Soulsby, Brewer, Smith, Dazzan, Berger, Yung, van den Buuse, Murray, McGorry and Pantelis2005). In samples diagnosed with psychotic disorder, an inverse association has been reported between illness duration and pituitary volume (Pariante et al. Reference Pariante, Vassilopoulou, Velakoulis, Phillips, Soulsby, Wood, Brewer, Smith, Dazzan, Yung, Zervas, Christodoulou, Murray, McGorry and Pantelis2004), although not consistently so (Upadhyaya et al. Reference Upadhyaya, El-Sheikh, MacMaster, Diwadkar and Keshavan2007). In addition, there is evidence for larger pituitary volume in first-degree relatives of patients with psychotic disorder (Mondelli et al. Reference Mondelli, Dazzan, Gabilondo, Tournikioti, Walshe, Marshall, Schulze, Murray, McDonald and Pariante2008).
The occurrence of pituitary enlargement in prodromal and antipsychotic (AP)-naïve patients, and in their relatives, suggests that the enlargement is not due to AP medication. However, cross-sectional studies show positive associations between first-generation APs (FGAs) or prolactin-elevating second-generation APs (SGAs) and pituitary volume (Pariante et al. Reference Pariante, Dazzan, Danese, Morgan, Brudaglio, Morgan, Fearon, Orr, Hutchinson, Pantelis, Velakoulis, Jones, Leff and Murray2005; MacMaster et al. Reference MacMaster, El-Sheikh, Upadhyaya, Nutche, Rosenberg and Keshavan2007) in FEP, although not all studies agree (Takahashi et al. Reference Takahashi, Suzuki, Velakoulis, Lorenzetti, Soulsby, Zhou, Nakamura, Seto, Kurachi and Pantelis2009). In addition, a longitudinal study reported a dose–response effect of prolactin-sparing SGAs on pituitary volume, higher doses predicting reduced volume (Nicolo et al. Reference Nicolo, Berger, Garner, Velakoulis, Markulev, Kerr, McGorry, Proffitt, McConchie, Pantelis and Wood2010). The mixed findings possibly reflect the fact that the pituitary gland is a highly dynamic organ, with many factors impacting on its volume over time. For example, studies have shown that women have larger pituitary volumes than men (Doraiswamy et al. Reference Doraiswamy, Potts, Axelson, Husain, Lurie, Na, Escalona, McDonald, Figiel and Ellinwood1992), and that both endogenous and exogenous estrogens are associated with pituitary volume (Grams et al. Reference Grams, Gempt, Stahl and Forschler2010).
The pituitary gland is a crucial organ involved in the hypothalamic–pituitary–adrenal (HPA) axis stress response. In depression, increased pituitary volume is thought to reflect HPA axis hyperfunctioning. Exposure to stress activates the HPA axis through the release of corticotropin-releasing hormone (CRH) in the hypothalamus, resulting in secretion of adrenocorticotropic hormone (ACTH) in the pituitary, stimulating the release of cortisol in the adrenal cortex. Cortisol, through a negative feedback mechanism, regulates HPA axis activity. Persistently elevated cortisol levels can impair the negative feedback mechanism (Walker & Diforio, Reference Walker and Diforio1997), resulting in HPA axis hyperactivity. As a result of HPA axis hyperactivity, patients become more sensitive to stressors, referred to as stress sensitization (Collip et al. Reference Collip, Myin-Germeys and van Os2008).
Increased cortisol levels have been found in FEP patients (Ryan et al. Reference Ryan, Sharifi, Condren and Thakore2004) and also in siblings of patients with psychotic disorder (Collip et al. Reference Collip, Nicolson, Lardinois, Lataster, van Os and Myin-Germeys2011), and may predict transition to a clinical psychotic outcome in young individuals at risk (Walker et al. Reference Walker, Brennan, Esterberg, Brasfield, Pearce and Compton2010). FEP patients in the acute phase also show non-suppression of cortisol secretion in the dexamethasone suppression test and in the dexamethasone/CRH test (Herz et al. Reference Herz, Fava, Molnar and Edwards1985; Lammers et al. Reference Lammers, Garcia-Borreguero, Schmider, Gotthardt, Dettling, Holsboer and Heuser1995). There is also evidence that patients with psychotic disorder are more sensitive to daily life stressors, resulting in increased levels of stress-related negative emotions (Myin-Germeys et al. Reference Myin-Germeys, van Os, Schwartz, Stone and Delespaul2001) and psychotic experiences (Myin-Germeys et al. Reference Myin-Germeys, Delespaul and van Os2005). Genetic risk for psychotic disorder may also be associated with increased cortisol reactivity (Collip et al. Reference Collip, Nicolson, Lardinois, Lataster, van Os and Myin-Germeys2011), although other evidence suggests a blunted cortisol response to metabolic and psychosocial stressors in patients (Marcelis et al. Reference Marcelis, Cavalier, Gielen, Delespaul and van Os2004; Brenner et al. Reference Brenner, Liu, Laplante, Lupien, Pruessner, Ciampi, Joober and King2009; van Venrooij et al. Reference van Venrooij, Fluitman, Lijmer, Kavelaars, Heijnen, Westenberg, Kahn and Gispen-de Wied2010).
To date, the association between pituitary volume and stress reactivity in psychotic disorder remains unclear. The present study is the first to examine pituitary volume in relation to cortisol levels and negative affect (NA) in the context of daily life stressors, using a genetically sensitive sample (patients at highest genetic risk, siblings of patients at higher than average genetic risk, and controls at average genetic risk). The hypothesis was that patients and, to a lesser extent, siblings would show increased emotional stress reactivity compared to controls, accompanied by larger pituitary volumes. Similarly, with regard to biological stress reactivity, we hypothesized that siblings would show increased cortisol reactivity, accompanied by larger pituitary volumes.
Method
Subjects
Data pertain to baseline measures of an ongoing longitudinal study in Maastricht, The Netherlands. Details on the selection procedures are described in a prior publication (Habets et al. Reference Habets, Marcelis, Gronenschild, Drukker and van Os2010). The sample consisted of 20 patients with a psychotic disorder, 37 siblings of patients with psychotic disorder and 32 controls. Inclusion criteria were: (i) age 16 to 50 years, (ii) diagnosis of non-affective psychotic disorder and (iii) good command of the Dutch language.
Diagnosis was based on DSM-IV criteria (APA, 2000), assessed with the Comprehensive Assessment of Symptoms and History (CASH) interview (Andreasen et al. Reference Andreasen, Flaum and Arndt1992). Patients' diagnoses were: schizophrenia (n=11), schizo-affective disorder (n=2), schizophreniform disorder (n=1), brief psychotic disorder (n=2), and psychotic disorder not otherwise specified (n=4). The CASH was also used to confirm the absence of a diagnosis of non-affective psychotic disorder in the siblings, and the absence of a lifetime diagnosis of any psychotic disorder or any current affective disorder in the healthy controls. Six siblings and seven controls had a history of major depressive disorder, but none presented in a current depressive episode.
Participants were screened for the following magnetic resonance imaging (MRI) exclusion criteria: (i) brain injury with unconsciousness (>1 h), (ii) neurological diseases that might have affected brain structure/function, (iii) cardiac arrhythmia requiring medical treatment (iv), severe claustrophobia, (v) metal corpora aliena, and (vi) women with intrauterine device status or (suspected) pregnancy.
The study was approved by the standing ethics committee, and all the subjects gave written informed consent in accordance with the committee's guidelines.
Measures
The Positive and Negative Syndrome Scale (PANSS; Kay et al. Reference Kay, Fiszbein and Opler1987) was used to measure psychotic symptoms over the past 2 weeks. AP medication use was determined using reports from the participant's psychiatrist. Best-estimate lifetime (cumulative) AP use was determined by multiplying the number of days of AP use by the corresponding haloperidol equivalents and summing these scores for all periods of AP use. Estrogen exposure was determined by multiplying the number of months of contraceptive medication use by micrograms of estrogen (of the corresponding contraceptive) per month. Substance use was assessed with the Composite International Diagnostic Interview (CIDI), sections B-J-L (WHO, 1990). Use of cannabis and other drugs was assessed as reported frequency of use during the past 12 months and lifetime. Alcohol intake was calculated as the average number of weekly consumptions during the past 12 months and tobacco use as the average number of daily cigarettes during the past 12 months.
Experience Sampling Method (ESM)
The ESM is a random time-sampling self-assessment technique; studies have demonstrated the feasibility, validity and reliability of the ESM in general and patient populations (Csikszentmihalyi & Larson, Reference Csikszentmihalyi and Larson1987; Myin-Germeys et al. Reference Myin-Germeys, Oorschot, Collip, Lataster, Delespaul and van Os2009). Subjects received a digital wristwatch that emitted a signal 10 times a day on 6 consecutive days, at unpredictable moments between 07:30 and 22:30 hours. After each ‘beep’, subjects completed ESM self-assessment forms concerning current context, thoughts, emotions, and psychotic experiences. Subjects were instructed to complete their reports immediately after the beep, thus minimizing memory distortions. Reports were considered valid when subjects responded within 15 min after the beep, as determined by comparing the actual beep time with the reported time of completion. For inclusion in the analyses, a valid response to at least one-third of the emitted beeps was required (Delespaul et al. Reference Delespaul, deVries and van Os2002).
ESM measures
In accordance with previous work, two different stress measures were computed: event stress and social stress (Myin-Germeys et al. Reference Myin-Germeys, van Os, Schwartz, Stone and Delespaul2001). In the event stress assessment, after each beep, the most important event that had occurred between the current and the previous beep was reported. This event was subsequently rated on a Likert scale (−3=very unpleasant, 0=neutral, 3=very pleasant). The responses were recoded to allow high scores to reflect stress. In the social stress assessment, participants were asked to evaluate the social context when other persons were present, by rating on a Likert scale (1=not true, 7=very true) the item ‘I don't like the company’ and the subsequent item ‘I would rather be alone’. The mean of these two items constituted the social stress scale. NA was assessed with five mood-related adjectives (down, guilty, insecure, lonely, anxious) scored on Likert scales (1=not at all, 7=very) after each beep. These were reduced to a single measure of mean NA with factor analyses (Cronbach's α=0.86).
Salivary cortisol sampling
After each ESM beep, siblings and controls collected a saliva sample with a cotton swab from a salivette tube (Salivette; Sarstedt, The Netherlands), placed the swab back in the salivette tube, and recorded the exact collection time. Samples were stored in subjects' home freezers until transport to the laboratory, where uncentrifuged samples were kept at −20°C until analysis. Samples collected more than 15 min after the beep were excluded from the analysis. Salivary cortisol is a reliable and non-invasive measure of the free, unbound cortisol in blood. Radioimmunoassays were run in duplicate, using a tracer solution of cortisol-3CMO coupled with 2-[125I]histamine and specific antibodies raised against cortisol-3CMO-BSA (Sulon et al. Reference Sulon, Demey-Ponsart, Beauduin and Sodoyez1978) (Dr J. Sulon, University of Liège). The lower detection limit of the assay was 0.2 nmol/l. The intra- and inter-assay coefficients of variation were <5% and <12% respectively. All samples from an individual were analyzed in the same assay to reduce sources of variability. Samples with cortisol >44 nmol/l (n=0) were considered physiologically abnormal and were excluded from the statistical analysis. As AP medication impacts cortisol values, cortisol was only measured in siblings and controls.
Stress reactivity
Stress reactivity, as described in previous work, was conceptualized as emotional (Myin-Germeys et al. Reference Myin-Germeys, van Os, Schwartz, Stone and Delespaul2001) and biological fluctuations (Collip et al. Reference Collip, Nicolson, Lardinois, Lataster, van Os and Myin-Germeys2011) associated with minor disturbances in daily life, based on ESM data. Emotional stress reactivity was defined as the change in NA in reaction to small daily life stressors; biological stress reactivity was defined as the change in cortisol level in reaction to these stressors. Within the ESM sampling framework, the level of stress refers to experience since the last beep whereas negative emotions refer to experience at the moment of the current beep; therefore, change in NA is measured in reaction to stress (stress reactivity) because stress precedes NA.
MRI acquisition and processing
Acquisition parameters of the MRI scans, acquired with a 3-T Siemens scanner, were as follows. For the modified driven equilibrium Fourier transform (MDEFT) sequence: 176 slices, 1 mm isotropic voxel size, echo time (TE) 2.4 ms, repetition time (TR) 7.92 ms, flip angle 15°, total acquisition time 12 min 51s. For the magnetization prepared rapid acquisition gradient echo (MP-RAGE) sequence developed for the Alzheimer's Disease Neuroimaging Initiative (ADNI): 192 slices, 1 mm isotropic voxel size, TE 2.6 ms, TR 2250 ms, flip angle 9°, total acquisition time 7 min 23 s. For both sequences the matrix size was 256×256 and the field of view 256×256 mm. Two sequences were used because of a scanner update during data collection. The MP-RAGE and MDEFT sequences are very similar, but to prevent systematic bias, the total proportion of MP-RAGE scans (27%) was balanced between groups.
MRI preprocessing
Scans were preprocessed with Freesurfer stable release v5.0. Technical details of these procedures are described in prior publications (Dale et al. Reference Dale, Fischl and Sereno1999; Fischl et al. Reference Fischl, Sereno and Dale1999).
Volume measures
Pituitary volumes were obtained after manual tracing using GIANT (General Image Analysis Tools developed by EHBMG; Gronenschild et al. Reference Gronenschild, Burgmans, Smeets, Vuurman, Uylings and Jolles2010), a customized software program that allows tracing of regions of interest in triplanar view and calculation of volumes of interest. All tracings were carried out on the native images. Each pituitary was traced in all coronal slices where it could be visualized (Sassi et al. Reference Sassi, Nicoletti, Brambilla, Harenski, Mallinger, Frank, Kupfer, Keshavan and Soares2001). Pituitary volume (in mm3) was calculated by summing the areas of the relevant slices and multiplying by the voxel size. All images were traced by the same rater (S.B.), who was blind to group assignment. The intra-class correlation coefficient of repeated measurements, based on a training set of 15 scans, was r=0.92.
Statistical analyses
ESM and cortisol data were analyzed using multilevel regression techniques, taking into account the hierarchical structure of the data. Repeated momentary measurements (level 1) were nested in subjects (level 2) who were part of the same family (level 3). Given these three levels, data were analyzed using the xtmixed multilevel random regression routine in Stata v. 11.0 (StataCorp, 2009). The B values are the fixed regression coefficients of the predictors in the multilevel model. Interaction terms were evaluated by Wald tests (Clayton & Hills, Reference Clayton and Hills1993). In case of significant interaction effects, stratified analyses were conducted to quantify group differences, using the Stata margin command to calculate the appropriate linear combinations from the model containing the interaction. All analyses were adjusted for the a priori hypothesized confounders age, sex, intracranial volume, estrogen exposure and scan type.
Group differences in pituitary volume were examined using multiple regression procedures, with pituitary volume as the dependent variable and group (dummy variable: 0=controls, 1=siblings, 2=patients) as the independent variable. Group differences in stress reactivity were assessed by multilevel analyses, with NA as the dependent variable and group, stress (event stress and social stress) and their interactions as the independent variables.
To test the hypothesis that group status modified the association between pituitary volume and stress reactivity, multilevel regression analyses were conducted with stress, pituitary volume and group, and their interaction terms, as independent variables, and NA as the dependent variable: NA=B0+B1(group)+B2(stress)+B3(volume)+B4(group×stress)+B5(group×volume)+B6(stress×volume)+B7(group×stress×volume). The group×stress×volume interaction was fitted with controls as the reference category. Pituitary volume was entered both as a linear variable and as dummy variables representing the distribution of volume calculated at the 33rd and 66th percentiles of pituitary volume of the controls, allowing visualization of dose–response: 1 corresponds to <33rd percentile (low pituitary volume), 2 to ⩾33rd and <66th percentile (medium pituitary volume), and 3 to ⩾66th percentile (high pituitary volume).
The same pattern of analyses as described above was applied to the log transformed cortisol values, with ln cort as the dependent variable. Additional confounders were: time of cortisol sample, time squared, and recent consumption of food or tobacco (since previous beep). To test whether the mean cortisol level differed between the two groups, a regression model was estimated with ln cort as the dependent variable and group (0=controls, 1=siblings) as the independent variable.
The association between pituitary volume and both exogenous estrogen and lifetime AP was investigated using multiple regression analysis.
Power analysis
Power calculations for the three-way interaction analyses were carried out by empirical statistical simulation in Stata (www.stata.com/support/faqs/stat/power.html), as described previously (Habets et al. Reference Habets, Marcelis, Gronenschild, Drukker and van Os2011). Effect sizes were based on previously published work in this area. The effect sizes for the association between group and NA and between stress and NA, and the interaction between group and stress, were adapted from Myin-Germeys et al. (Reference Myin-Germeys, van Os, Schwartz, Stone and Delespaul2001) . Similarly, effect sizes in the model of cortisol were adapted from Collip et al. (Reference Collip, Nicolson, Lardinois, Lataster, van Os and Myin-Germeys2011) . The effect sizes for the other interactions were not known. Therefore, effect sizes of 0.2 standard deviation (s.d.) (small) were used for the patient versus control simulations, and effect sizes of 0.1 s.d. were used for the sibling versus control simulations. At an α of 0.05, the three-way interaction analysis had a power of 56% to detect a significant difference between patients and controls and a power of 18% to detect a significant difference between siblings and controls. Therefore, the results of the three-way interactions measuring differences between siblings and controls are only described exploratively so as to generate hypotheses for future research and help to generate realistic power calculations based on published effect sizes in these groups.
Sensitivity analysis
Studies report an association between pituitary volume and illness duration in schizophrenia (Pariante et al. Reference Pariante, Vassilopoulou, Velakoulis, Phillips, Soulsby, Wood, Brewer, Smith, Dazzan, Yung, Zervas, Christodoulou, Murray, McGorry and Pantelis2004; Tournikioti et al. Reference Tournikioti, Tansella, Perlini, Rambaldelli, Cerini, Versace, Andreone, Dusi, Balestrieri, Malago, Gasparini and Brambilla2007), with increased pituitary volume for FEP patients and decreased pituitary volume for patients with established illness (EIP). Therefore, explorative analyses were carried out separately in patients with an illness duration of ⩽5 years (FEP; n=10) and >5 years (EIP; n=10).
Results
Descriptive analyses
Groups were well matched on most demographic variables (Table 1). Patients smoked more cigarettes and cannabis and used more hard drugs (lifetime) than siblings and controls, with no difference between the latter two groups. Eighteen patients received AP medication (SGA: n=2; FGA: n=16). Furthermore, three patients and two controls used antidepressants and one patient and one control used benzodiazepines. Patients reported higher PANSS scores than controls and siblings, with no difference between the latter two groups (Table 1).
Table 1. Subject demographics

PANSS, Positive and Negative Syndrome Scale; MDEFT, modified driven equilibrium Fourier transform; ADNI, Alzheimer's Disease Neuroimaging Initiative; s.d., standard deviation.
F/χ2 and p values refer to between-group differences.
Group differences in pituitary volume and emotional stress reactivity
Mean pituitary volume did not differ significantly between groups (patients: 602.6 mm3, s.e.=23.4; siblings: 624.5 mm3, s.e.=17.5; controls: 624.4 mm3, s.e.=18.0; B=5.34, p=0.76). There was a trend significant group×event stress interaction in the model of NA (Table 2), indicating that stress reactivity differed between groups. Stratified analyses revealed that patients had significantly higher levels of stress reactivity compared to both controls and siblings. For social stress, there was no significant interaction effect with group in the model of NA (Table 2).
Table 2. Multilevel regression estimates of emotional and biological stress reactivity

NA, Negative affect; C, Controls; S, siblings; P, patients.
* p<0.05.
χ2 and p values represent the results of the Wald test, testing the significance of stratified effects calculated from the model with the group×stress interaction on cortisol and NA respectively.
The B coefficients represent the effect sizes of level of stress sensitivity compared to the reference level.
Interaction between genetic risk, emotional stress reactivity and pituitary volume
In the model of NA, there was a significant social stress×volume (linear variable) interaction (χ2=23.03, p=0.00) and also a significant group×social stress×volume, the latter indicating that the association between pituitary volume and stress reactivity differed between groups (χ2=20.13, p=0.00). Stratified analyses showed that, with increasing pituitary volume, stress reactivity was progressively greater in patients compared to both controls and siblings, with no difference between the latter two groups. Visualizing the effect over the three tertile groups of pituitary volume revealed progressively higher stress reactivity with progressively higher pituitary volume in all three groups, but most pronounced in patients (Table 3; Fig. 1). Stress reactivity was significantly higher in the high-volume compared to the low-volume group in controls (χ2=7.17, p=0.01), siblings (χ2=5.34, p=0.02) and patients (χ2=12.72, p=0.00). In patients, stress reactivity was also significantly increased in the medium-volume compared to the low-volume group (χ2=30.62, p=0.00).

Fig. 1. Stress reactivity (expressed as B coefficients) as a function of pituitary volume, stratified by group. Stress reactivity was significantly higher in the high-volume group when compared to the low volume in controls (p=0.01), siblings (p=0.02) and patients (p=0.00). When comparing middle volume with low volume, only patients showed increased stress reactivity (p=0.00). Patients showed higher stress reactivity when compared to controls in the medium pituitary volume group (p=0.00) and in the high pituitary volume group (p=0.02).
Table 3. Negative affect (NA) and cortisol as a function of daily life stress and pituitary volume, stratified by group

CI, Confidence interval.
The B values represent the effect sizes of stress reactivity compared to the reference level.
χ2 and p values represent the results of the Wald test, testing the significance of stratified effects calculated from the group×stress×volume interaction in the models of NA and cortisol.
The event stress×volume interaction effect in the model of NA was not significant, nor was the group×event stress×volume (linear variable) interaction (χ2=1.15, p=0.56) (Table 3); however, effect sizes were directionally similar to those observed for social stress. The direction and significance of the results remained the same when estrogen exposure was removed as a confounding factor (results available upon request).
Group differences in cortisol levels and biological stress reactivity
Siblings had significantly higher cortisol levels than controls (B=0.30, p=0.00). The group×stress interaction in the model of cortisol was not significant for either social or event-related stress (Table 2).
Interaction between genetic risk, biological stress reactivity and pituitary volume
The stress×volume interaction and the group×stress×volume interaction in the model of cortisol were not significant for either social (χ2=1.04, p=0.31) or event-related stress (χ2=1.49, p=0.22) (Table 3).
AP medication use and estrogen exposure in relation to pituitary volume
Lifetime AP use did not predict pituitary volume (B=1.04, p=0.19), nor did estrogen exposure (B=0.00, p=0.48). There was no suggestive or significant interaction between group and estrogen exposure in the model of pituitary volume (χ2=0.00, p=0.99). Women had larger pituitary volumes than men (B=63.94, p=0.04), but no group×gender interaction was apparent in the model of pituitary volume.
Sensitivity analysis
Although the results were not statistically significant, FEP patients had higher pituitary volumes [mean =641.1 mm3, standard error (s.e.)=33.79] than controls, and EIP patients had lower pituitary volumes (mean=564 mm3, s.e.=91.08) than controls. Post-hoc analyses revealed that the group×social stress×volume (linear variable) interaction was attributable predominantly to the FEP (χ2=22.66, p=0.00), whereas in the EIP this interaction was no longer seen, although effect sizes were directionally similar (χ2=1.65, p=0.44).
Discussion
There was no mean difference in pituitary volume among the three groups. Patients with psychotic disorder showed higher levels of emotional stress reactivity than siblings and controls. Emotional stress reactivity was moderated not only by group but also by pituitary volume. In all three groups, increased emotional stress reactivity was associated with higher pituitary volume, but this association was much more pronounced in the patient group, suggesting a process of progressive sensitization mediating clinical outcome.
Emotional stress reactivity
The finding of increased emotional stress reactivity in patients is in line with previous studies that have shown that minor stressors in the flow of daily life were associated with an increase in NA in individuals with a psychotic disorder (Myin-Germeys et al. Reference Myin-Germeys, van Os, Schwartz, Stone and Delespaul2001, Reference Myin-Germeys, Oorschot, Collip, Lataster, Delespaul and van Os2009). Although there was no difference in pituitary volume between the groups, emotional stress reactivity was associated with a larger pituitary volume in all three groups, with a much more pronounced increase in the patients than in the other two groups (Table 3). More specifically, each unit increase in stress reactivity was associated with an increase in pituitary volume of 0.20 s.d. in the patients. This phenomenon may indicate a mechanism of stress sensitization in the patients and mediate the onset of psychotic disorder.
Stress sensitization occurs when early and/or frequent exposures to stressors alter the stress response system, sensitizing individuals to later stress, leading to increased biological and behavioral response to that stressor, even when less severe than the original stressor. Given that many environmental risk factors associated with psychotic disorder may be linked to stress, stress sensitization may represent a common mechanism linking multiple environmental exposures (Collip et al. Reference Collip, Myin-Germeys and van Os2008). For example, evidence suggests that patients with psychotic disorder who report previous exposure to childhood trauma (Lardinois et al. Reference Lardinois, Lataster, Mengelers, van Os and Myin-Germeys2010) or other stressful life events (Myin-Germeys et al. Reference Myin-Germeys, Krabbendam, Delespaul and van Os2003) show increased stress reactivity in daily life. Higher pituitary volume may be the result of differential functional states (e.g. increased stress reactivity/hypercortisolism) or, vice versa, contribute causally to increased stress reactivity. However, as this is a cross-sectional study, no definite conclusion can be drawn regarding the direction of this association. Post-hoc analyses revealed that FEP had higher and EIP had lower pituitary volumes than controls, in line with previous studies (Pariante et al. Reference Pariante, Dazzan, Danese, Morgan, Brudaglio, Morgan, Fearon, Orr, Hutchinson, Pantelis, Velakoulis, Jones, Leff and Murray2005; Upadhyaya et al. Reference Upadhyaya, El-Sheikh, MacMaster, Diwadkar and Keshavan2007). In addition, the association between increased stress reactivity and pituitary volume was confined to the FEP group. These observations are in agreement with the hypothesis of two biologically distinct pathways resulting in mean group differences between (i) type 1 psychotic disorder, characterized by high levels of positive symptoms and an episodic course (Myin-Germeys & van Os, Reference Myin-Germeys and van Os2007; Lataster et al. Reference Lataster, Collip, Lardinois, van Os and Myin-Germeys2010), in addition to high stress reactivity and associated larger pituitary volume, and (ii) type 2 psychotic disorder, characterized by high levels of negative symptoms, cognitive impairments and a chronic course (Robins & Guze, Reference Robins and Guze1970; Crow, Reference Crow1980; Murray et al. Reference Murray, O'Callaghan, Castle and Lewis1992; van Os et al. Reference van Os, Gilvarry, Bale, Van Horn, Tattan, White and Murray1999), in addition to low stress reactivity and smaller pituitary volumes. A larger study sample and a longitudinal design are needed to draw more definitive conclusions about underlying mechanisms and direction of effects.
Biological stress reactivity
Mean cortisol levels were higher in siblings than in controls, consistent with another study from our group (Collip et al. Reference Collip, Nicolson, Lardinois, Lataster, van Os and Myin-Germeys2011). It has been suggested that elevated cortisol levels in patients with a diagnosis of psychotic disorder (Ryan et al. Reference Ryan, Sharifi, Condren and Thakore2004; Mondelli et al. Reference Mondelli, Dazzan, Hepgul, Di Forti, Aas, D'Albenzio, Di Nicola, Fisher, Handley, Marques, Morgan, Navari, Taylor, Papadopoulos, Aitchison, Murray and Pariante2010) and their siblings reflect disturbances in the negative feedback regulation of the HPA axis (Phillips et al. Reference Phillips, McGorry, Garner, Thompson, Pantelis, Wood and Berger2006). A comparable observation in depression is that HPA axis hyperactivity is thought to reflect a lack of negative inhibitory feedback by circulating cortisol on the HPA axis, particularly at the level of pituitary cells producing ACTH, leading to an increase in the size and number of these cells. Persistently elevated cortisol secretion could also lead to resistance of the glucocorticoid receptors for cortisol (Pariante, Reference Pariante2006).
Cortisol fluctuations in response to daily life stress, however, did not show a significant difference between siblings and controls, which does not agree with a previous study (Collip et al. Reference Collip, Nicolson, Lardinois, Lataster, van Os and Myin-Germeys2011). The present analyses, although based on the same cohort as described by Collip et al. (Reference Collip, Nicolson, Lardinois, Lataster, van Os and Myin-Germeys2011), were carried out in individuals for whom MRI data were available, resulting in a sample size reduction. Thus, statistical power was insufficient to detect small differences in cortisol reactivity between siblings and controls and/or as a function of pituitary volume.
In conclusion, although some evidence for HPA axis hyperactivity was found in siblings (increased mean cortisol values), cortisol reactivity was not mediated by group status or pituitary volume in this study. Larger studies are needed before more definitive conclusions can be drawn. Nevertheless, the biological and emotional stress reactivity pathways seemed to be associated, as similar patterns in both biological and emotional stress reactivity were found in siblings and controls.
Methodological considerations
Some authors have cast doubt on the reliability and subject compliance in paper-and-pencil ESM studies, favoring the use of electronic devices (Stone et al. Reference Stone, Shiffman, Schwartz, Broderick and Hufford2002). However, in a comparative study, Green et al. (Reference Green, Rafaeli, Bolger, Shrout and Reis2006) concluded that both methods yielded similar results.
Two different scanning sequences were used in this study, which could in theory have influenced the volume measurements. However, it has been shown that within-scanner measurements were reliable for subcortical volumes across scan sessions (Jovicich et al. Reference Jovicich, Czanner, Han, Salat, van der Kouwe, Quinn, Pacheco, Albert, Killiany, Blacker, Maguire, Rosas, Makris, Gollub, Dale, Dickerson and Fischl2009). Furthermore, no large or significant interaction was found between scan type and group on pituitary volume (χ2=0.09, p=0.96). Similarly, adjusting the analyses for scanning sequence did not affect the direction or significance of the results.
Numerous external factors can influence pituitary volume, such as the use of APs. According to the study by Nicolo et al. (Reference Nicolo, Berger, Garner, Velakoulis, Markulev, Kerr, McGorry, Proffitt, McConchie, Pantelis and Wood2010) in FEP, SGAs reduce pituitary gland volume in a dose-dependent manner. Other studies report that prolactin-elevating medication increases pituitary volume in comparison with prolactin-sparing medication (Pariante et al. Reference Pariante, Dazzan, Danese, Morgan, Brudaglio, Morgan, Fearon, Orr, Hutchinson, Pantelis, Velakoulis, Jones, Leff and Murray2005; Mondelli et al. Reference Mondelli, Dazzan, Gabilondo, Tournikioti, Walshe, Marshall, Schulze, Murray, McDonald and Pariante2008), but not all studies have found this association (Pariante et al. Reference Pariante, Vassilopoulou, Velakoulis, Phillips, Soulsby, Wood, Brewer, Smith, Dazzan, Yung, Zervas, Christodoulou, Murray, McGorry and Pantelis2004; Takahashi et al. Reference Takahashi, Suzuki, Velakoulis, Lorenzetti, Soulsby, Zhou, Nakamura, Seto, Kurachi and Pantelis2009). The present study did not find an association between cumulative lifetime exposure to AP medication and pituitary volume, but the groups were too small to examine differential effects of the AP group (SGA versus FGA). Adding lifetime AP medication as a confounder did not change the results (data available upon request).
Exogenous estrogen exposure has also been associated with pituitary volume. Reduced pituitary volumes have been found in women taking oral contraceptives (Grams et al. Reference Grams, Gempt, Stahl and Forschler2010) but, by contrast, elevated pituitary volumes have been found in menopausal women using estrogens (Abech et al. Reference Abech, Moratelli, Leite and Oliveira2005). The present study did not find an association between lifetime estrogen exposure and pituitary volume. Larger studies investigating the association between estrogens and pituitary volume longitudinally are warranted. Nevertheless, there was a gender effect, in that pituitary volumes were higher in women than in men, possibly reflecting differences in endogenous estrogen levels (Doraiswamy et al. Reference Doraiswamy, Potts, Axelson, Husain, Lurie, Na, Escalona, McDonald, Figiel and Ellinwood1992). The role of estrogen in psychotic disorder remains unknown. Evidence suggests that estrogens may exert a neuroprotective effect and possibly have AP properties (Boerma et al. Reference Boerma, van der Stel, van Amelsvoort, Linszen and de Haan2010). However, if estrogen could also lead to pituitary enlargement, this may potentially be less beneficial in terms of vulnerability for psychotic disorder. It is advisable to take estrogen measures into account when examining pituitary volumes.
Acknowledgments
This work was sponsored by the Netherlands Organization for Scientific Research (NOW), the European Community's Seventh Framework Programme under grant agreement HEALTH-F2-2009-241909 (EU-GEI consortium), and by a 2006 National Alliance for Research on Schizophrenia and Depression (NARSAD) Young Investigator Award and by the Dutch Medical Research Council (VENI and VIDI grant) to I. Myin-Germeys.
We thank T. Driesen, I. Crolla and F. van Goethem for their assistance in this project as well as the G.R.O.U.P. investigators (R. Kahn, D. Linszen, J. van Os, D. Wiersma, R. Bruggeman, W. Cahn, L. de Haan, L. Krabbendam and I. Myin-Germeys).
Declaration of Interest
J. van Os is or has been an unrestricted research grant holder with, or has received financial compensation as an independent symposium speaker from, Eli Lilly, BMS, Lundbeck, Organon, Janssen, GlaxoSmithKline, AstraZeneca, Pfizer, and Servier. M. Marcelis has received financial compensation as an independent symposium speaker from Eli Lilly and Janssen.