Abstract

STUDY QUESTION

Do mental distress and mood fluctuations in women undergoing GnRH agonist and GnRH antagonist protocols for assisted reproductive technology (ART) differ depending on protocol and the personality trait, neuroticism?

SUMMARY ANSWER

ART treatment did not induce elevated levels of mental distress in either GnRH antagonist or agonist protocols but neuroticism was positively associated with increased mental distress, independent of protocols.

WHAT IS KNOWN ALREADY

ART treatment may increase mental distress by mechanisms linked to sex hormone fluctuations. General psychological characteristics, such as personality traits indexing negative emotionality, e.g. neuroticism, are likely to affect mental distress during ART treatment.

STUDY DESIGN, SIZE, DURATION

A total of 83 women undergoing their first ART cycle were consecutively randomized 1:1 to GnRH antagonist (n = 42) or GnRH agonist (n = 41) protocol. The study population was a subgroup of a larger ongoing Danish clinical randomized trial and was established as an add-on in the period 2010–2012.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Women in the GnRH antagonist protocol received daily injections with recombinant follicle-stimulating hormone, Puregon® and subcutaneous injections with GnRH antagonist, Orgalutran®. Women in the GnRH agonist protocol received nasal administration of the GnRH agonist, Synarela® and subcutaneous injections with FSH, Puregon®. The study design did not allow for a blinding procedure. All women self-reported the Profile of Mood States, the Perceived Stress Scale, the Symptom Checklist-92-Revised, and the Major Depression Inventory questionnaires, at baseline, at ART cycle day 35, on the day of oocyte pick-up, and on the day of hCG testing. Also, a series of Profile of Mood States were reported daily during pharmacological treatment to monitor mood fluctuations. The personality trait Neuroticism was assessed at baseline by the self-reported NEO-PI-R questionnaire.

MAIN RESULTS AND THE ROLE OF CHANCE

ART did not induce within- or between-protocol changes in any of the applied measures of mental distress. However, the GnRH antagonist protocol was associated with more pronounced median mood fluctuations during the stimulation phase (antagonist, 11.0 SD, [IQR = 21.1–6.1]; agonist, 8.9 SD, [IQR = 11.3–5.7], P = 0.025). This association became non-significant after applying a Bonferroni–Holm correction. Neuroticism was highly positively associated with increased levels of mental distress throughout treatment independent of protocols (all P-values <0.006), and cross-sectional analysis revealed that women with high or low Neuroticism scores at baseline showed a significant trend towards lower chances of a positive pregnancy test (P-value =0.028).

LIMITATIONS, REASONS FOR CAUTION

Information on prognostic factors such as preceding length of infertility, number of retrieved oocytes and number of prior insemination treatments was not accounted for in the analyses. The stratification of protocols by age in the subgroups of women included in this study was suboptimal. Women with prior or current use of antidepressant medication were excluded from our study.

WIDER IMPLICATIONS

Our results imply that mental distress emerging during ART treatment is not causally linked to hypogonadism per se or to the choice of protocol. Rather, our data highlight the potential importance of (i) rapid increases in ovarian steroids and (ii) addressing personality traits indexing negative emotionality, i.e. Neuroticism, in women undergoing ART treatment, to optimize both emotional adjustment and, possibly, the chances of obtaining pregnancy.

STUDY FUNDING/COMPETING INTEREST(S)

The Danish Research Council for Independent Research and MSD, Denmark kindly supported the study. The authors declare no competing financial interests.

TRIAL REGISTRATION NUMBER

EudraCT – 2008-005452-24.

Introduction

Women and men with fertility problems often experience their situation as highly stressful (Cousineau and Domar, 2007). The strain of undergoing assisted reproductive technologies (ART) is also well-documented (Hammarberg et al., 2001; Verhaak et al., 2007), reflected by the fact that psychological burden is a major reason for discontinuing treatment (Olivius et al., 2004; Gameiro et al., 2012). Such adverse psychological effects of ART treatment could have a negative impact on both mother and infant future health (O'hara and Swain, 1996; Eugster and Vingerhoets, 1999; Brouwers et al., 2001; Hjelmstedt et al., 2003; Alder et al., 2007), and currently it is not clear if mental distress experienced by women undergoing ART treatment is associated with chances of obtaining pregnancy (Boivin et al., 2011; Matthiesen et al., 2011). ART treatment may increase mental distress in women through a complex set of biological and psychological factors; controlled ovarian stimulation in ART is achieved by the pharmacological induction of sex steroid hormone changes, which could influence levels of mental distress as sex steroid hormones may play a role in the regulation of mood and the pathophysiology of mood disorders (Payne, 2003; Douma et al., 2005; Freeman et al., 2006, 2014; Munk-Olsen et al., 2006; Deecher et al., 2008). Especially, the administration of gonadotrophin-releasing hormone (GnRH) agonists has been coupled to negative mood symptoms, typically attributed to medically induced hypogonadism (Warnock et al., 1998; Patten and Barbui, 2004). Interestingly, adverse responses to changes in sex hormone levels appear only to affect a subgroup of women undergoing ART treatment (Van den Broeck et al., 2010; Bloch et al., 2011). Therefore identifying individual markers of susceptibility is pivotal to enable targeted prevention. The psychological construct of Neuroticism indexes individual differences in the tendency to experience negative emotions, impulsiveness, anxiety and angry hostility, and vulnerability to stress (McCrae and Costa, 2003). High Neuroticism scores have been associated with risk for major depression (Kendler and Myers, 2010), also in the context of stress (Jacobs et al., 2006), psychopathology in general (Ormel et al., 2013), and increased release of the stress hormone cortisol (Portella et al., 2005; Nater et al., 2010). However, the role of Neuroticism for ART-induced mental distress is at present unknown. In this study, we tested the following hypotheses: (i) women in the GnRH agonist protocol exhibit increased levels of mental distress compared with women in the GnRH antagonist protocol; (ii) mood fluctuations are more pronounced in the GnRH agonist protocol compared with the GnRH antagonist protocol; and (iii) neuroticism interacts with the protocols and independently affects levels of mental distress and mood fluctuations.

Method

Participants

A total of 83 eligible women (mean age: 33.1 ± 4.8, range: 22–39 years) undergoing their first ART cycle were consecutively randomized 1:1 to a GnRH antagonist (n = 42) or GnRH agonist (n = 41) protocol by a project nurse, stratified by age (≤36 years and >36 years) and the need for ICSI or general IVF. The study population was a subgroup of a larger ongoing Danish clinical randomized trial evaluating treatment outcome of the two protocols and was established as an add-on to the overall study in the period 2010–2012. Exclusion criteria were: prior IVF treatment, uterine anomalies, testicular sperm aspiration (TESA) needed, allergy to the ingredients used in the pharmacological treatment, reduced kidney or liver function, women >40 years of age, or prior or current use of antidepressant medication. The study was registered as a clinical trial (EudraCT – 2008-005452-24) and approved by the Ethics Committee for the Capital Region of Denmark (H-B-2008-109). All participants signed an informed consent form.

Intervention

Women in the GnRH antagonist protocol received daily injections with the recombinant follicle-stimulating hormone (rFSH) analogue, Puregon®, to initiate ovarian stimulation (150, i.e. ≤36 years, and 225, i.e. >36 years, respectively) starting at cycle day 2–3. After 5 days of stimulation treatment, the women received additional daily subcutaneous injections with the GnRH antagonist, Orgalutran® (1 × 0.25 mg). Women in the GnRH agonist protocol received daily nasal administration of the GnRH agonist, Syranela®, to suppress ovarian hormone production (200 mg × 3 daily) starting at cycle day 21. After 14 days of GnRH agonist administration (cycle day 35), the women received additional daily injections with Puregon® (150, i.e. ≤36 years and 225, i.e. >36 years, respectively). Nasal administration of Syranela® was continued (200 mg × 2 daily) until the day of oocyte pick-up. In both protocols, ovulation induction was induced by subcutaneous injection with Ovitrelle® (6500, i.e.), when the three largest ovarian follicles had a diameter ≥17 mm. Oocyte retrieval was performed 36–38 h later.

Measures

Time points for data collection

Figure 1 is an overview of the study design. Participating women completed their questionnaires using a secure online survey system (https://survey.nru.dk/). Baseline questionnaires (T0) were completed at home on the day of inclusion and randomization. A computer was set up at the fertility clinic, which the women accessed to complete their questionnaires when routinely attending the clinic for treatment at three time points: when ovarian stimulation was initiated (T1) [at cycle day 35 for women in the agonist protocol]; on the day of oocyte pick-up (T2); and on the day of hCG-testing (T3). In addition, mood was reported daily during the pharmacological treatment from home. At T2, questionnaires were completed before the aspiration of follicles, and likewise, questionnaires were completed prior to hCG-testing at T3. Baseline plasma measures of FSH, luteinizing hormone (LH) and plasma for determining estradiol concentrations were collected on the day of inclusion (T0). Subsequent plasma samples were collected on the days described above as attendance points for questionnaires (T1, T2 and T3).

Overview of the treatment course, cycle days and time points of assessment in the two protocols. POMS: Profile of Mood States; Serial POMS: Series of daily reported mood disturbances scores (POMS) during pharmacological treatment; CD: cycle day.
Figure 1

Overview of the treatment course, cycle days and time points of assessment in the two protocols. POMS: Profile of Mood States; Serial POMS: Series of daily reported mood disturbances scores (POMS) during pharmacological treatment; CD: cycle day.

The NEO Personality Inventory

At baseline (T0), all women completed the Danish version (Skovdahl Hansen et al., 2003) of the NEO Personality Inventory (NEO PI-R) (Costa and McCrae, 1992). The NEO PI-R is a self-report inventory, which measures five major domains of personality with six sub-facets for each domain. To serve the purpose of this study, the Neuroticism domain was used. The Neuroticism domain comprises 48 items (e.g. ‘I am not a worrier’ or ‘I often get angry at the way people treat me’) rated on a 5-point Likert scale from 0 (strongly disagree) to 4 (strongly agree). Individuals who score high on Neuroticism tend to experience difficulties in coping with stress and interpret the world around them as threatening and frustrating, which make them more prone to experience anxiety, anger, guilt, stress or sadness. In the present sample, internal consistency (Cronbach's alpha, α) for the Neuroticism domain was high, α = 0.91.

The Profile of Mood States

The Profile of Mood States (POMS) is a psychological rating scale used to assess transient, distinct mood states (McNair et al., 1992). It consists of six factors and a total score of mood disturbance rated by 65 adjectives (e.g. ‘Furious’, ‘Hopeless’ and ‘Carefree’) on a 5-point Likert scale from 1 (not at all) to 5 (extremely) based on the recollection of ‘the last 24 h’ for serial daily reports at home and ‘past week including present day’ for attendance at the clinic, respectively. For the purpose of this study, only the total mood disturbance scores were used. The POMS was completed at T0, T1, T2 and T3 and daily during the intervention. Internal consistency for the total mood disturbance score was high, α = 0.85.

The Perceived Stress Scale

The Perceived Stress Scale (PSS) is a rating scale, which provides an overall estimation of the degree to which respondents experience their lives as unpredictable, uncontrollable and overloaded (Cohen et al., 1983). It consists of 10 stress-related items (e.g. ‘How often have you felt that you were unable to control the important things in your life?’) rated on a 5-point Likert scale from 0 (never) to 4 (very often) based on the recollection of the last 2 weeks. The PSS was completed at T0, T1, T2 and T3. Internal consistency for the PSS was high, α = 0.91.

The Symptom Checklist Revised

The Danish version of the Symptom Checklist-90-Revised (SCL-92-R) (Olsen et al., 2006) was used to assess severity of mental distress. The SCL-92-R comprises 92 items (e.g. ‘difficulty making decision’ or ‘feeling afraid to go out of your house alone’) rated on a 5-point Likert scale of distress from 0 (none) to 4 (extreme). Nine primary symptom scales and three global indices of distress are derived. For the purpose of this study, only the Global Severity Index (GSI) was used. The SCL-92-R was completed at T0, T1, T2 and T3. Internal consistency for the GSI was high, α = 0.95.

The Major Depression Inventory

The Major Depression Inventory (MDI) is a rating scale of depressive symptoms according to DSM-IV and ICD-10 diagnostic criteria (Bech et al., 2001). It comprises 10 items (e.g. ‘Have you felt very restless?’ or ‘Have you suffered from reduced appetite?’) rated on a 6-point Likert scale from 0 (never) to 5 (all the time) based on the recollection of the last 2 weeks. The MDI was completed at T0, T1, T2 and T3. Internal consistency for the MDI was high, α = 0.88.

Reported measures

We use ‘mental distress’ as an umbrella term covering mood disturbances, perceived stress, global symptom distress and depressive symptoms, except when otherwise explicitly stated. Mood fluctuations are not included in this term and refer to the standard deviation (SD) of the serial reported mood disturbances scores for each woman during pharmacological treatment (see also under statistical analyses).

Plasma estradiol

Venous blood was drawn from all women at T0, T1, T2 and T3. Plasma was kept at a temperature between 2 and 8°C until analysed on a routine basis at the hospital laboratory after a maximum of 7 days. Estradiol concentrations were determined by electrochemiluminescence immunoassays (ECLIA) on a Cobas E601 Immunology Analyzer (Roche, Mannheim, Germany) with a lower detection limit of 0.02 nmol/l.

Statistical analyses

Main analyses

Demographic data were analysed using Student's t-tests or Mann–Whitney's U-tests for continuous data and Chi-square tests (χ2) for categorical data. Friedman's test was used to examine ART-induced changes in mental distress and plasma estradiol within the protocols and Kruskal–Wallis one-way analysis of variance of the pair wise change (from T0 – T2) was used to examine differences in ART-induced changes in mental distress between the protocols. For interaction analyses, we used a general linear model adjusted for baseline measures of mental distress, age and BMI. Likewise, we used a general linear model for our cross-sectional analyses of Neuroticism effects on mental distress adjusted for age and BMI. Mann–Whitney's U-test was used to compare protocol differences in mood fluctuations. Levels of mood disturbances over the course of treatment fluctuated around a quite stationary mean in both protocols. This was evaluated using generalized estimating equations with independence working correlation structure showing no significant linear trends (all P-values >0.22). We therefore assigned the SD around each woman's own mean as an index of fluctuations in serial daily reported mood disturbances scores (POMS) for the phases shown in Fig. 1. We controlled for the family-wise error-rate using the Bonferroni–Holm procedure (Holm, 1979). All presented P-values in the manuscript and the tables are uncorrected P-values.

Follow-up analysis

We conducted a follow-up cross-sectional analysis to examine whether baseline Neuroticism predicted probability of positive pregnancy test (as defined by hCG > 50 ui/l). This was analysed using logistic regression allowing a non-linear function of Neuroticism and measures of mental distress modelled flexibly using a natural cubic spline with a single knot at the median of the predictor adjusted for age, protocol and BMI (Harrell, 2001). All analyses were performed in SPSS (version 20.0) and R (version 3.0.1) (Team, 2013). We used a significance level of ≤0.05 and all hypothesis-tests were considered two-sided.

Results

Baseline data

Descriptive data for women in the GnRH antagonist and GnRH agonist protocols are shown in Table I. Median age differed significantly between the two protocols (antagonist, 31.2 years [IQR = 35.5–28.4]; agonist, 36.4 years [IQR = 37.6–32.7], P = 0.001). No women presented with clinical levels of psychopathology at baseline, according to established Danish criteria for clinical cut off scores on the SCL-92-R (Olsen, Mortensen and Bech, 2006), and no significant protocol differences were observed on measures of mental distress (Table I). Plasma concentrations of estradiol were not associated with measures of mental distress or Neuroticism at baseline (all P-values > 0.251); however, Neuroticism scores were, as expected, positively associated with mental distress at baseline (all P-values <0.001, two-tailed).

Table I

Descriptive baseline data.

Descriptive variablesAntagonistAgonistP-value
Age in years31.2
[35.5–28.4]
36.4
[37.6–32.7]
0.001
Body mass index (kg/m2)24.9 ± 4.725.0 ± 5.60.817
Smoking (cigarettes per day)0.50.60.705
Alcohol use (units per week)11.20.433
Neuroticism scores86.9 ± 22.084.7 ± 21.10.506
Mental distress
 POMS8.0
[44.0–0.5]
8.0
[21.00–5.00]
0.197
 PSS13.5
[20.3–7.0]
12.0
[17.00–5.00]
0.298
 MDI6.0
[16.3–3.8]
6.0
[9.00–4.00]
0.407
 SCL-92-R0.3
[0.6–0.2]
0.3
[0.5–0.2]
0.529
Descriptive variablesAntagonistAgonistP-value
Age in years31.2
[35.5–28.4]
36.4
[37.6–32.7]
0.001
Body mass index (kg/m2)24.9 ± 4.725.0 ± 5.60.817
Smoking (cigarettes per day)0.50.60.705
Alcohol use (units per week)11.20.433
Neuroticism scores86.9 ± 22.084.7 ± 21.10.506
Mental distress
 POMS8.0
[44.0–0.5]
8.0
[21.00–5.00]
0.197
 PSS13.5
[20.3–7.0]
12.0
[17.00–5.00]
0.298
 MDI6.0
[16.3–3.8]
6.0
[9.00–4.00]
0.407
 SCL-92-R0.3
[0.6–0.2]
0.3
[0.5–0.2]
0.529

Information obtained at baseline with medians and interquartile ranges in square brackets or means ± standard deviations and P-values of protocol differences.

POMS, Profile of Mood States; PSS, Perceived Stress Scale; MDI, Major Depression Inventory; SCL-92-R, Symptom Checklist 92-Revised.

Table I

Descriptive baseline data.

Descriptive variablesAntagonistAgonistP-value
Age in years31.2
[35.5–28.4]
36.4
[37.6–32.7]
0.001
Body mass index (kg/m2)24.9 ± 4.725.0 ± 5.60.817
Smoking (cigarettes per day)0.50.60.705
Alcohol use (units per week)11.20.433
Neuroticism scores86.9 ± 22.084.7 ± 21.10.506
Mental distress
 POMS8.0
[44.0–0.5]
8.0
[21.00–5.00]
0.197
 PSS13.5
[20.3–7.0]
12.0
[17.00–5.00]
0.298
 MDI6.0
[16.3–3.8]
6.0
[9.00–4.00]
0.407
 SCL-92-R0.3
[0.6–0.2]
0.3
[0.5–0.2]
0.529
Descriptive variablesAntagonistAgonistP-value
Age in years31.2
[35.5–28.4]
36.4
[37.6–32.7]
0.001
Body mass index (kg/m2)24.9 ± 4.725.0 ± 5.60.817
Smoking (cigarettes per day)0.50.60.705
Alcohol use (units per week)11.20.433
Neuroticism scores86.9 ± 22.084.7 ± 21.10.506
Mental distress
 POMS8.0
[44.0–0.5]
8.0
[21.00–5.00]
0.197
 PSS13.5
[20.3–7.0]
12.0
[17.00–5.00]
0.298
 MDI6.0
[16.3–3.8]
6.0
[9.00–4.00]
0.407
 SCL-92-R0.3
[0.6–0.2]
0.3
[0.5–0.2]
0.529

Information obtained at baseline with medians and interquartile ranges in square brackets or means ± standard deviations and P-values of protocol differences.

POMS, Profile of Mood States; PSS, Perceived Stress Scale; MDI, Major Depression Inventory; SCL-92-R, Symptom Checklist 92-Revised.

Changes in mental distress

Contrary to our hypotheses, we observed no significant ART-induced, within protocol changes in mental distress from baseline (T0) to the day of oocyte pick-up (T2), including the suppression phase for women in GnRH agonist protocol (T1) (Table II). Furthermore, we observed no differences in changes of mental distress from baseline (T0) to the day of oocyte pick-up (T2) between the agonist and antagonist protocols (all P-values >0.125). We conducted a post hoc Friedman's test, which further included the day of hCG-testing (T3) for the subgroup of women who had embryos transferred (Fig. 1). Here, a significant effect was observed for women in the GnRH antagonist protocol (Table II). Analyses of pair wise change (Kruskal–Wallis) showed a decrease in median perceived stress (PSS) from the day of oocyte pick-up (T2) to the day of hCG-testing (T3) for these women (T2, 13.0 PSS [IQR = 18.0–7.0]; T3, 10.0 PSS [IQR = 16.0–5.8], P = 0.036). Changes in estradiol during pharmacological treatment were assessed by calculating a change score (Δ) for each woman (T3 − T0 for all women and T1 − T0 for women in the agonist protocol only) (Δantagonist, 3.30 ± 2.57 nmol/l; Δagonist, 3.74 ± 2.57 and −0.11 ± 0.09 nmol/l, respectively). Neither the absolute concentrations of estradiol nor the magnitude of changes during pharmacological treatment were associated with levels of mental distress (all P-values >0.203).

Table II

Mental distress and estradiol during ART treatment.

Mental distressFrom baseline to the end of hormone stimulation
Including day of hCG-testing
ProtocolBaseline T0Hypogonadal T1Oocyte pick-up T2X2P-valuehCG-test T3X2P-value
POMSAntagonist8.0
[44.0–0.5]

11.0
[33.5–1.5]
0.6760.5069.0
[54.0–1.0]
0.8380.658
Agonist8.0
[21.0–5.0]
1.0
[14.8–8.0]
0.0
[15.0–8.0]
1.8600.40910.0
[23.0–5.0]
3.9130.271
PSSAntagonist13.5
[20.3–7.0]

13.0
[18.0–8.0]
0.2730.7319.0
[16.0–6.0]
7.7070.029
Agonist12.0
[17.0–5.0]
9.0
[15.0–6.3]
10.0
[12.0–6.0]
0.5950.7499.0
[15.0–6.0]
1.5410.673
MDIAntagonist6.0
[16.3–3.8]

7.0
[13.5–4.0]
0.1110.8647.0
[14.0–3.0]
0.7030.704
Agonist6.0
[9.0–4.0]
6.5
[11.8–3.0]
6.0
[11.0–4.0]
1.7500.4178.0
[11.0–5.0]
2.2860.515
SCL-92-RAntagonist0.3
[0.6–0.2]

0.3
[0.6–0.2]
0.0271.0000.3
[0.7–0.2]
1.4610.482
Agonist0.3
[0.5–0.2]
0.3
[0.5–0.1]
0.2
[0.4–0.1]
0.2680.8870.3
[0.4–0.2]
3.2710.352
Estradiol (nmol/l)Antagonist0.14
[0.24–0.11]

3.40
[4.22–2.19]
36.100<0.0010.35
[2.50–0.12]
30.260<0.001
Agonist0.12
[0.18–0.09]
0.04
[0.06–0.00]
3.12
[5.11–2.07]
63.510<0.0010.18
[4.21–0.12]
67.070<0.001
Mental distressFrom baseline to the end of hormone stimulation
Including day of hCG-testing
ProtocolBaseline T0Hypogonadal T1Oocyte pick-up T2X2P-valuehCG-test T3X2P-value
POMSAntagonist8.0
[44.0–0.5]

11.0
[33.5–1.5]
0.6760.5069.0
[54.0–1.0]
0.8380.658
Agonist8.0
[21.0–5.0]
1.0
[14.8–8.0]
0.0
[15.0–8.0]
1.8600.40910.0
[23.0–5.0]
3.9130.271
PSSAntagonist13.5
[20.3–7.0]

13.0
[18.0–8.0]
0.2730.7319.0
[16.0–6.0]
7.7070.029
Agonist12.0
[17.0–5.0]
9.0
[15.0–6.3]
10.0
[12.0–6.0]
0.5950.7499.0
[15.0–6.0]
1.5410.673
MDIAntagonist6.0
[16.3–3.8]

7.0
[13.5–4.0]
0.1110.8647.0
[14.0–3.0]
0.7030.704
Agonist6.0
[9.0–4.0]
6.5
[11.8–3.0]
6.0
[11.0–4.0]
1.7500.4178.0
[11.0–5.0]
2.2860.515
SCL-92-RAntagonist0.3
[0.6–0.2]

0.3
[0.6–0.2]
0.0271.0000.3
[0.7–0.2]
1.4610.482
Agonist0.3
[0.5–0.2]
0.3
[0.5–0.1]
0.2
[0.4–0.1]
0.2680.8870.3
[0.4–0.2]
3.2710.352
Estradiol (nmol/l)Antagonist0.14
[0.24–0.11]

3.40
[4.22–2.19]
36.100<0.0010.35
[2.50–0.12]
30.260<0.001
Agonist0.12
[0.18–0.09]
0.04
[0.06–0.00]
3.12
[5.11–2.07]
63.510<0.0010.18
[4.21–0.12]
67.070<0.001

Friedman's test of within protocol changes at the various time points with medians and interquartile ranges in square brackets (see also Fig. 1).

POMS, Profile of Mood States; PSS, Perceived Stress Scale; MDI, Major Depression Inventory; SCL-92-R, Symptom Checklist 92-Revised.

Table II

Mental distress and estradiol during ART treatment.

Mental distressFrom baseline to the end of hormone stimulation
Including day of hCG-testing
ProtocolBaseline T0Hypogonadal T1Oocyte pick-up T2X2P-valuehCG-test T3X2P-value
POMSAntagonist8.0
[44.0–0.5]

11.0
[33.5–1.5]
0.6760.5069.0
[54.0–1.0]
0.8380.658
Agonist8.0
[21.0–5.0]
1.0
[14.8–8.0]
0.0
[15.0–8.0]
1.8600.40910.0
[23.0–5.0]
3.9130.271
PSSAntagonist13.5
[20.3–7.0]

13.0
[18.0–8.0]
0.2730.7319.0
[16.0–6.0]
7.7070.029
Agonist12.0
[17.0–5.0]
9.0
[15.0–6.3]
10.0
[12.0–6.0]
0.5950.7499.0
[15.0–6.0]
1.5410.673
MDIAntagonist6.0
[16.3–3.8]

7.0
[13.5–4.0]
0.1110.8647.0
[14.0–3.0]
0.7030.704
Agonist6.0
[9.0–4.0]
6.5
[11.8–3.0]
6.0
[11.0–4.0]
1.7500.4178.0
[11.0–5.0]
2.2860.515
SCL-92-RAntagonist0.3
[0.6–0.2]

0.3
[0.6–0.2]
0.0271.0000.3
[0.7–0.2]
1.4610.482
Agonist0.3
[0.5–0.2]
0.3
[0.5–0.1]
0.2
[0.4–0.1]
0.2680.8870.3
[0.4–0.2]
3.2710.352
Estradiol (nmol/l)Antagonist0.14
[0.24–0.11]

3.40
[4.22–2.19]
36.100<0.0010.35
[2.50–0.12]
30.260<0.001
Agonist0.12
[0.18–0.09]
0.04
[0.06–0.00]
3.12
[5.11–2.07]
63.510<0.0010.18
[4.21–0.12]
67.070<0.001
Mental distressFrom baseline to the end of hormone stimulation
Including day of hCG-testing
ProtocolBaseline T0Hypogonadal T1Oocyte pick-up T2X2P-valuehCG-test T3X2P-value
POMSAntagonist8.0
[44.0–0.5]

11.0
[33.5–1.5]
0.6760.5069.0
[54.0–1.0]
0.8380.658
Agonist8.0
[21.0–5.0]
1.0
[14.8–8.0]
0.0
[15.0–8.0]
1.8600.40910.0
[23.0–5.0]
3.9130.271
PSSAntagonist13.5
[20.3–7.0]

13.0
[18.0–8.0]
0.2730.7319.0
[16.0–6.0]
7.7070.029
Agonist12.0
[17.0–5.0]
9.0
[15.0–6.3]
10.0
[12.0–6.0]
0.5950.7499.0
[15.0–6.0]
1.5410.673
MDIAntagonist6.0
[16.3–3.8]

7.0
[13.5–4.0]
0.1110.8647.0
[14.0–3.0]
0.7030.704
Agonist6.0
[9.0–4.0]
6.5
[11.8–3.0]
6.0
[11.0–4.0]
1.7500.4178.0
[11.0–5.0]
2.2860.515
SCL-92-RAntagonist0.3
[0.6–0.2]

0.3
[0.6–0.2]
0.0271.0000.3
[0.7–0.2]
1.4610.482
Agonist0.3
[0.5–0.2]
0.3
[0.5–0.1]
0.2
[0.4–0.1]
0.2680.8870.3
[0.4–0.2]
3.2710.352
Estradiol (nmol/l)Antagonist0.14
[0.24–0.11]

3.40
[4.22–2.19]
36.100<0.0010.35
[2.50–0.12]
30.260<0.001
Agonist0.12
[0.18–0.09]
0.04
[0.06–0.00]
3.12
[5.11–2.07]
63.510<0.0010.18
[4.21–0.12]
67.070<0.001

Friedman's test of within protocol changes at the various time points with medians and interquartile ranges in square brackets (see also Fig. 1).

POMS, Profile of Mood States; PSS, Perceived Stress Scale; MDI, Major Depression Inventory; SCL-92-R, Symptom Checklist 92-Revised.

Mood fluctuations

To allow comparison between the phases in the treatment period, the GnRH agonist protocol was divided into a suppression and stimulation phase. During the stimulation phase, women in the GnRH antagonist protocol exhibited significantly more pronounced median mood fluctuation scores (SD's) relative to women in the GnRH agonist protocol (antagonist, 11.0 SD, [IQR = 21.1–6.1]; agonist, 8.9 SD, [IQR = 11.3–5.7], P = 0.025), which reflects a difference between protocols of ∼20%. No difference was observed between the stimulation phase in GnRH antagonist protocol and the suppression phase in GnRH agonist protocol (Table III). Applying a Bonferroni–Holm correction rendered the observed difference in mood fluctuations non-significant.

Table III

Mood fluctuations during pharmacological treatment.

Protocol and phaseSuppression phase
Stimulation phase
P-value
DaysSDDaysSD
Antagonist (N = 39)1011.0
[21.1–6.1]
Agonist (N = 37)1410.3
[14.9–7.0]
128.9
[11.3–5.7]
Antagonist Stim. versus Agonist Supp.0.531
Antagonist Stim. versus Agonist Stim.0.025
Protocol and phaseSuppression phase
Stimulation phase
P-value
DaysSDDaysSD
Antagonist (N = 39)1011.0
[21.1–6.1]
Agonist (N = 37)1410.3
[14.9–7.0]
128.9
[11.3–5.7]
Antagonist Stim. versus Agonist Supp.0.531
Antagonist Stim. versus Agonist Stim.0.025

Mean duration (days) and median standard deviations (SD) and interquartile ranges in square brackets of daily mood disturbances scores (POMS) during pharmacological treatment with P-values of protocol differences.

POMS, Profile of Mood States; Stim, stimulation phase; Supp., suppression phase.

Table III

Mood fluctuations during pharmacological treatment.

Protocol and phaseSuppression phase
Stimulation phase
P-value
DaysSDDaysSD
Antagonist (N = 39)1011.0
[21.1–6.1]
Agonist (N = 37)1410.3
[14.9–7.0]
128.9
[11.3–5.7]
Antagonist Stim. versus Agonist Supp.0.531
Antagonist Stim. versus Agonist Stim.0.025
Protocol and phaseSuppression phase
Stimulation phase
P-value
DaysSDDaysSD
Antagonist (N = 39)1011.0
[21.1–6.1]
Agonist (N = 37)1410.3
[14.9–7.0]
128.9
[11.3–5.7]
Antagonist Stim. versus Agonist Supp.0.531
Antagonist Stim. versus Agonist Stim.0.025

Mean duration (days) and median standard deviations (SD) and interquartile ranges in square brackets of daily mood disturbances scores (POMS) during pharmacological treatment with P-values of protocol differences.

POMS, Profile of Mood States; Stim, stimulation phase; Supp., suppression phase.

Neuroticism

Neuroticism significantly predicted mental distress throughout treatment (all P-values < 0.006), but did not interact with protocols (all P-values > 0.115). Higher Neuroticism scores at baseline were also significantly associated with more pronounced mood fluctuations during the stimulation phases across protocols (P = 0.035). In the cross-sectional follow-up analysis, a non-linear inversed U-shaped association was found, such that women with high or low Neuroticism scores at baseline showed a significant trend (P-value = 0.028) towards a lower probability of pregnancy (defined as hCG > 50 ui/l) (Fig. 2).

A non-linear association was observed, such that women with high or low Neuroticism scores at baseline showed a significant trend (P = 0.028) towards lower chances of pregnancy (defined as HCG > 50). This was examined using a logistic regression model with the effects of Neuroticism modelled using a natural cubic spline basis with a single knot located at the median. Although similar trends were observed using different statistical approaches, significance levels varied in both directions with the number and placement of knots, polynomial degree, and type of test (Wald/Likelihood Ratio Test). The grey area indicates the 95% confidence interval.
Figure 2

A non-linear association was observed, such that women with high or low Neuroticism scores at baseline showed a significant trend (P = 0.028) towards lower chances of pregnancy (defined as HCG > 50). This was examined using a logistic regression model with the effects of Neuroticism modelled using a natural cubic spline basis with a single knot located at the median. Although similar trends were observed using different statistical approaches, significance levels varied in both directions with the number and placement of knots, polynomial degree, and type of test (Wald/Likelihood Ratio Test). The grey area indicates the 95% confidence interval.

Discussion

We observed no significant changes in mental distress from baseline to the end of hormone stimulation, or a protocol × Neuroticism interaction, within or between the GnRH agonist and antagonist protocols. However, mood fluctuations were more pronounced during the stimulation phase for women in the GnRH antagonist protocol versus the GnRH agonist protocol. High Neuroticism scores were highly significantly associated with increased mental distress throughout treatment independent of protocols and with mood fluctuations during the stimulation phase. Also, a follow-up analysis suggested that women with high or low Neuroticism scores showed a significant trend towards a lower chance of pregnancy.

ART protocols and changes in mental distress

Contrary to our a priori hypothesis, the initial suppression phase in the GnRH agonist protocol was not associated with elevated symptoms of mental distress compared with baseline, suggesting that the initial suppression phase per se does not induce mental distress. Consistent with our findings, no exacerbations of mood symptoms during the hypogonadal phase were observed in women undergoing IVF treatment compared with baseline (Bloch et al., 2011). However, these findings contrast with an earlier study, where elevated symptoms of depression were observed during the hypogonadal phase in women undergoing their first cycle of conventional IVF treatment when compared with a mild ovarian stimulation with single embryo transfer protocol (de Klerk et al., 2006), and a smaller study, where symptoms of depression and anxiety progressively increased from baseline to the day of oocyte pick-up (Toren et al., 1996). The discrepant findings may to some extent reflect the use of different psychometric tools across studies. De Klerk et al. (2006) found elevated depressive symptoms during pituitary down-regulation using the self-reported Hospital Anxiety and Depression Scale (HADS) in women undergoing ART, unlike the present study and the study by Bloch et al. (2011), where no elevated levels of distress were found using SCL-92-R based measures. Toren et al. (1996) used the interview-based Hamilton Rating Scale and found elevated depressive symptoms during pituitary down-regulation corroborating de Klerk et al.'s (2006)' findings. The HADS assesses only non-physical symptoms of anxiety and depression, while the SCL-92-R assesses a broad range of psychopathological symptoms and associated levels of global distress, making the latter more suitable for a global characterization of mental distress and the former more suitable for focusing on anxiety and depression in a medical population. Therefore it is possible that the SCL-92-R is less specific in detecting changes in mood and anxiety symptoms in women undergoing ART when compared with HADS. However, the women studied here presented with reproductive difficulties but were otherwise healthy and as de Klerk et al. (2006) did not exclude women with prior or current use of antidepressant medication, direct comparison between studies is difficult.

Mood fluctuations during ART treatment

Our data are consistent with the hypothesis that controlled ovarian hormone stimulation, in the absence of prior suppression of ovarian hormone production, is associated with more pronounced mood fluctuations. Rodent and non-human primate models suggest that both withdrawal from and exposure to ovarian hormones can potently affect brain functions of relevance to emotional regulation and mood, such as serotonergic signalling (Bethea et al., 2002; Lu et al., 2003; Suda et al., 2008). However, it is far from clear how rapid hormone fluctuations may affect mental distress and through which mechanisms or which phases of fluctuations are the more critical (Ben Dor et al., 2013). Based on the current observations, we therefore speculate that prior treatment with GnRH agonists, at the time suppression is established, may dampen the potential adverse effects of subsequent controlled ovarian stimulation. However, correcting for multiple comparisons rendered the observation non-significant and replication is needed in order to validate this possible inter-protocol difference in mood fluctuations using a design with serial daily reports. The stimulation phase coincides with initiation of treatment in GnRH antagonist as opposed to GnRH agonist protocol. Therefore factors not induced pharmacologically (e.g. anxiety and expectation pressures related to the procedure) may also contribute to the observed differences in mood fluctuations.

The role of neuroticism

Neuroticism (i.e. negative emotions and stress vulnerability, impulsive and labile mood dispositions) was highly positively associated with increased levels of mental distress and with more pronounced mood fluctuations during the stimulation phase in both protocols. The latter suggests that higher levels of neuroticism combined with rapid increases in ovarian hormones may elicit adverse mood symptoms, such as unstable and labile mood. Corroborating this, similar personality constructs to Neuroticism have been coupled to subtle menstrual cycle disturbances (Demyttenaere et al., 1994), levels of ovarian reproductive hormones during the menstrual cycle (Ziomkiewicz et al., 2012), and infertility-related distress (Van den Broeck et al., 2010). High or low Neuroticism scores at baseline showed a significant trend towards subsequent negative pregnancy test. As Neuroticism is thought to be a risk marker for mental distress and depression, this is consistent with a meta-analysis on emotional distress and outcome of ART treatment in which perceived stress and trait/state anxiety showed a negative association with clinical pregnancy, and clinical depression showed a non-significant negative trend towards lower pregnancy rates (Matthiesen et al., 2011). Likewise, an adverse effect of Neuroticism and a positive effect of reducing stress on probability of pregnancy and live birth have been observed (Hämmerli et al., 2009; Volgsten et al., 2010). However, a recent meta-analysis did not report any association between pretreatment emotional distress and outcome of ART treatment (Boivin et al., 2011). The relatively few observations in the lower end of the Neuroticism spectrum in this study make interpretation of the association between low Neuroticism scores and chances of pregnancy somewhat inconclusive, as also reflected by the broad confidence interval (Fig. 2).

Of particular relevance for clinicians handling the treatment of subfertile women, personality appears to be a stronger predictor than choice of protocol for the mental distress experienced by women undergoing ART. If replicated, the ease of administration and the clinical relevance make ratings of Neuroticism an applicable tool for health care professionals in evaluating ART-related vulnerability to mental distress. Our tentative results also warrant further longitudinal investigation of the effect of personality and mental distress for successful ART treatment outcome, ideally with clinical end-points, such as number of live children born, and infant health measures.

Methodological considerations

The strength of this study is that it is part of a controlled randomized and prospective study with an extensive set-up of well-validated and reliable psychometric instruments covering the entire treatment cycle, including daily reports during pharmacological treatment. The study also took into account general psychological characteristics associated with increased risk of mood disorders, namely the personality trait, Neuroticism. However, the present findings should be interpreted under the following potentially important methodological limitations. First, even though only women undergoing their first ART-cycle were included, information on prognostic factors such as preceding length of infertility and number of prior inseminations was not accounted for in the analyses. Second, the stratification of protocols by age (≤36 years and >36 years) in this subgroup of women was not successful. Hence, women in the GnRH antagonist protocol were younger on average than those in the GnRH agonist protocol. Moreover, the different doses of Puregon® administered across protocols (150, i.e. ≤36 years and 225, i.e. >36 years, respectively) may have biased our results. However, age, group and BMI were included as covariates in our cross-sectional analyses. Third, women with prior or current use of antidepressant medication were excluded from our study. This is likely to have biased our study population towards low-risk women with lower Neuroticism scores. We may therefore have underestimated the effects of the stimulation phase on mood fluctuations and the predictive value of Neuroticism for mental distress and probability of pregnancy.

In conclusion, ART treatment did not induce elevated levels of mental distress and no differences between agonist and antagonist protocols or interactions with Neuroticism were observed. However, in the absence of prior ovarian hormone suppression (i.e. GnRH antagonist protocol), more pronounced mood fluctuations during the stimulation phase were observed, which was rendered non-significant after correcting for multiple comparisons. Neuroticism predicted levels of mental distress and mood fluctuations during treatment independent of protocols and was associated with probability of achieving pregnancy.

Authors’ roles

D.S.S.: contributed to conception and design of the study and the acquisition of data; managed literature searches; organized and conducted the statistical analyses; wrote the first draft of the manuscript; contributed substantially to analyses and interpretation of the results; and revised the manuscript critically for important intellectual content. M.T.: contributed to conception and design of the study and the acquisition of data; carried out medical examinations of the participants; was involved in the statistical analyses; contributed substantially to analyses and interpretation of the results; and revised the manuscript critically for important intellectual content. L.V.H.: contributed to conception and design of the study; was involved in the statistical analyses; contributed substantially to analyses and interpretation of the results; and revised the manuscript critically for important intellectual content. P.S.J.: contributed to conception and design of the study; contributed considerably to analyses and interpretation of the results; and revised the manuscript critically for important intellectual content. K.K.H.: contributed considerably to analyses and interpretation of the results; conducted all modelling of the applied non-linear models; and revised the manuscript critically for important intellectual content. T.B.: contributed to conception and design of the study; carried out medical examinations of the participants; contributed considerably to analyses and interpretation of the results; and revised the manuscript critically for important intellectual content. T.H.: recruited all participants; handled data collection; and revised the manuscript critically for important intellectual content. J.B.: carried out medical examinations of the participants; substantially contributed to analyses and interpretation of data; and revised the manuscript critically for important intellectual content. A.P.: carried out medical examinations of the participants; substantially contributed to analyses and interpretation of data; and revised the manuscript critically for important intellectual content. P.H.: carried out medical examinations of the participants; substantially contributed to analyses and interpretation of data; and revised the manuscript critically for important intellectual content. V.G.F.: conceptualized and designed the study in collaboration with co-authors; participated in acquisition of data; managed literature searches; organized and contributed to the statistical analyses; contributed substantially to analyses and interpretation of the results; participated in manuscript drafting; and revised the manuscript critically for important intellectual content.

Funding

The Danish Research Council for Independent Research, The Lundbeck Foundation, The Capital Region of Denmark, Foundation for Health Research, and MSD, Denmark kindly supported the study.

Conflict of interest

None declared.

Acknowledgements

The authors thank the participating women.

References

Alder
J
Fink
N
Bitzer
J
Hösli
I
Holzgreve
W
Depression and anxiety during pregnancy: a risk factor for obstetric, fetal and neonatal outcome? A critical review of the literature
J Matern Fetal Neonatal Med
2007
, vol. 
20
 (pg. 
189
-
209
)
Bech
P
Rasmussen
NA
Olsen
LR
Noerholm
V
Abildgaard
W
The sensitivity and specificity of the Major Depression Inventory, using the Present State Examination as the index of diagnostic validity
J Affect Disord
2001
, vol. 
66
 (pg. 
159
-
164
)
Ben Dor
R
Harsh
VL
Fortinsky
P
Koziol
DE
Rubinow
DR
Schmidt
PJ
Effects of pharmacologically induced hypogonadism on mood and behavior in healthy young women
Am J Psychiatry
2013
, vol. 
170
 (pg. 
426
-
433
)
Bethea
CL
Lu
NZ
Gundlah
C
Streicher
JM
Diverse actions of ovarian steroids in the serotonin neural system
Front Neuroendocrinol
2002
, vol. 
23
 (pg. 
41
-
100
)
Bloch
M
Azem
F
Aharonov
I
Ben Avi
I
Yagil
Y
Schreiber
S
Amit
A
Weizman
A
GnRH-agonist induced depressive and anxiety symptoms during in vitro fertilization-embryo transfer cycles
Fertil Steril
2011
, vol. 
95
 (pg. 
307
-
309
)
Boivin
J
Griffiths
E
Venetis
CA
Emotional distress in infertile women and failure of assisted reproductive technologies: meta-analysis of prospective psychosocial studies
BMJ
2011
pg. 
342
 
Brouwers
EPM
van Baar
AL
Pop
VJM
Maternal anxiety during pregnancy and subsequent infant development
Infant Behav Dev
2001
, vol. 
24
 (pg. 
95
-
106
)
Cohen
S
Kamarck
T
Mermelstein
R
A global measure of perceived stress
J Health Soc Behav
1983
, vol. 
24
 (pg. 
385
-
396
)
Costa
PT
McCrae
RR
Professional Manual for Revised NEO Personality Inventory
1992
Odessa, Florida
Psychological Assessment Resources
Cousineau
TM
Domar
AD
Psychological impact of infertility
Best Pract Res Clin Obstet Gynaecol
2007
, vol. 
21
 (pg. 
293
-
308
)
Deecher
D
Andree
TH
Sloan
D
Schechter
LE
From menarche to menopause: exploring the underlying biology of depression in women experiencing hormonal changes
Psychoneuroendocrinology
2008
, vol. 
33
 (pg. 
3
-
17
)
de Klerk
C
Heijnen
EM
Macklon
NS
Duivenvoorden
HJ
Fauser
BC
Passchier
J
Hunfeld
JA
The psychological impact of mild ovarian stimulation combined with single embryo transfer compared with conventional IVF
Hum Reprod
2006
, vol. 
21
 (pg. 
721
-
727
)
Demyttenaere
K
Nijs
P
Evers-Kiebooms
G
Koninckx
PR
Personality characteristics, psychoendocrinological stress and outcome of IVF depend upon the etiology of infertility
Gynecol Endocrinol
1994
, vol. 
8
 (pg. 
233
-
240
)
Douma
SL
Husband
C
O'Donnell
ME
Barwin
BN
Woodend
AK
Estrogen-related mood disorders: reproductive life cycle factors
ANS Adv Nurs Sci
2005
, vol. 
28
 (pg. 
364
-
375
)
Eugster
A
Vingerhoets
AJJM
Psychological aspects of in vitro fertilization: a review
Soc Sci Med
1999
, vol. 
48
 (pg. 
575
-
589
)
Freeman
EW
Sammel
MD
Lin
H
Nelson
DB
Associations of hormones and menopausal status with depressed mood in women with no history of depression
Arch Gen Psychiatry
2006
, vol. 
63
 (pg. 
375
-
382
)
Freeman
EW
Sammel
MD
Boorman
DW
Zhang
R
Longitudinal pattern of depressive symptoms around natural menopause
JAMA Psychiatry
2014
, vol. 
71
 (pg. 
36
-
43
)
Gameiro
S
Boivin
J
Peronace
L
Verhaak
CM
Why do patients discontinue fertility treatment? A systematic review of reasons and predictors of discontinuation in fertility treatment
Hum Reprod Update
2012
, vol. 
18
 (pg. 
652
-
669
)
Hammarberg
K
Astbury
J
Baker
H
Women's experience of IVF: a follow-up study
Hum Reprod
2001
, vol. 
16
 (pg. 
374
-
383
)
Hämmerli
K
Znoj
H
Barth
J
The efficacy of psychological interventions for infertile patients: a meta-analysis examining mental health and pregnancy rate
Hum Reprod Update
2009
, vol. 
15
 (pg. 
279
-
295
)
Harrell
FE
Jr
Regression Modeling Strategies: With Applications to Linear Models, Logistic Regression and Survival Analysis
2001
New York
Springer
Hjelmstedt
A
Widström
A-M
Wramsby
H
Matthiesen
A-S
Collins
A
Personality factors and emotional responses to pregnancy among IVF couples in early pregnancy: a comparative study
Acta Obstet Gynecol Scand
2003
, vol. 
82
 (pg. 
152
-
161
)
Holm
S
A simple sequentially rejective multiple test procedure
Scand J Stat
1979
, vol. 
6
 (pg. 
65
-
70
)
Jacobs
N
Kenis
G
Peeters
F
Derom
C
Vlietinck
R
van Os
J
Stress-related negative affectivity and genetically altered serotonin transporter function: evidence of synergism in shaping risk of depression
Arch Gen Psychiatry
2006
, vol. 
63
 (pg. 
989
-
996
)
Kendler
KS
Myers
J
The genetic and environmental relationship between major depression and the five-factor model of personality
Psychol Med
2010
, vol. 
40
 (pg. 
801
-
806
)
Lu
NZ
Eshleman
AJ
Janowsky
A
Bethea
CL
Ovarian steroid regulation of serotonin reuptake transporter (SERT) binding, distribution, and function in female macaques
Mol Psychiatry
2003
, vol. 
8
 (pg. 
353
-
360
)
Matthiesen
SMS
Frederiksen
Y
Ingerslev
HJ
Zachariae
R
Stress, distress and outcome of assisted reproductive technology (ART): a meta-analysis
Hum Reprod
2011
, vol. 
26
 (pg. 
2763
-
2776
)
McCrae
RR
Costa
PT
Personality in Adulthood : A Five-Factor Theory Perspective
2003
2nd edn
New York
Guilford Press
McNair
DM
Lorr
M
Droppleman
L
Manual for the Profile of Mood States
1992
San Diego, CA
EdITS/Educational and Industrial Testing Service
Munk-Olsen
T
Laursen
T
Pedersen
C
Mors
O
Mortensen
P
New parents and mental disorders: a population-based register study
JAMA
2006
, vol. 
296
 (pg. 
2582
-
2589
)
Nater
UM
Hoppmann
C
Klumb
PL
Neuroticism and conscientiousness are associated with cortisol diurnal profiles in adults—role of positive and negative affect
Psychoneuroendocrinology
2010
, vol. 
35
 (pg. 
1573
-
1577
)
O'hara
MW
Swain
AM
Rates and risk of postpartum depression—a meta-analysis
Int Rev Psychiatry
1996
, vol. 
8
 (pg. 
37
-
54
)
Olivius
C
Friden
B
Borg
G
Bergh
C
Why do couples discontinue in vitro fertilization treatment? A cohort study
Fertil Steril
2004
, vol. 
81
 (pg. 
258
-
261
)
Olsen
LR
Mortensen
EL
Bech
P
Mental distress in the Danish general population
Acta Psychiatr Scand
2006
, vol. 
113
 (pg. 
477
-
484
)
Ormel
J
Jeronimus
BF
Kotov
R
Riese
H
Bos
EH
Hankin
B
Rosmalen
JG
Oldehinkel
AJ
Neuroticism and common mental disorders: meaning and utility of a complex relationship
Clin Psychol Rev
2013
, vol. 
33
 (pg. 
686
-
697
)
Patten
SB
Barbui
C
Drug-induced depression: a systematic review to inform clinical practice
Psychother Psychosom
2004
, vol. 
73
 (pg. 
207
-
215
)
Payne
JL
The role of estrogen in mood disorders in women
Int Rev Psychiatry
2003
, vol. 
15
 (pg. 
280
-
290
)
Portella
MJ
Harmer
CJ
Flint
J
Cowen
P
Goodwin
GM
Enhanced early morning salivary cortisol in neuroticism
Am J Psychiatry
2005
, vol. 
162
 (pg. 
807
-
809
)
Skovdahl Hansen
H
Mortensen
EL
Schiøtz
HK
NEO PI-R. Manual – Clinical. Danish Version
2003
Copenhagen
Dansk Psykologisk Forlag
Suda
S
Segi-Nishida
E
Newton
SS
Duman
RS
A postpartum model in rat: behavioral and gene expression changes induced by ovarian steroid deprivation
Biol Psychiatry
2008
, vol. 
64
 (pg. 
311
-
319
)
Team
RC
R: A Language and Environment for Statistical Computing
2013
Vienna, Austria
R Foundation for Statistical Computing
Toren
P
Dor
J
Mester
R
Mozes
T
Blumensohn
R
Rehavi
M
Weizman
A
Depression in women treated with a gonadotropin-releasing hormone agonist
Biol Psychiatry
1996
, vol. 
39
 (pg. 
378
-
382
)
Van den Broeck
U
D'Hooghe
T
Enzlin
P
Demyttenaere
K
Predictors of psychological distress in patients starting IVF treatment: infertility-specific versus general psychological characteristics
Hum Reprod
2010
, vol. 
25
 (pg. 
1471
-
1480
)
Verhaak
CM
Smeenk
JM
Nahuis
MJ
Kremer
JA
Braat
DD
Long-term psychological adjustment to IVF/ICSI treatment in women
Hum Reprod
2007
, vol. 
22
 (pg. 
305
-
308
)
Volgsten
H
Ekselius
L
Poromaa
IS
Svanberg
AS
Personality traits associated with depressive and anxiety disorders in infertile women and men undergoing in vitro fertilization treatment
Acta Obstet Gynecol Scand
2010
, vol. 
89
 (pg. 
27
-
34
)
Warnock
JK
Bundren
JC
Morris
DW
Depressive symptoms associated with gonadotropin-releasing hormone agonists
Depress Anxiety
1998
, vol. 
7
 (pg. 
171
-
177
)
Ziomkiewicz
A
Wichary
S
Bochenek
D
Pawlowski
B
Jasienska
G
Temperament and ovarian reproductive hormones in women: evidence from a study during the entire menstrual cycle
Horm Behav
2012
, vol. 
61
 (pg. 
535
-
540
)