Premenstrual dysphoric disorder (PMDD) is a serious condition that causes significant suffering in 2%–8% of women of reproductive age worldwide and their families (
1,
2).
PMDD is characterized by affective, cognitive, behavioral, and somatic symptoms that occur consistently during the luteal phase of the menstrual cycle. For a diagnosis, the symptoms must be present during most menstrual cycles for at least 1 year and must be severe enough to impact daily functioning (
1). The impairment related to the affective and somatic symptoms in PMDD is particularly striking when expressed in terms of disease burden, calculated according to the World Health Organization model that determines disability-adjusted life years (i.e., the years of life lost to premature death or lived with a disability). Based on the 2000 U.S. census and a PMDD prevalence of 5%, researchers (
3,
4) estimated the number of symptomatic menstrual cycles a woman with PMDD would experience during her reproductive life, and they determined that the disability-adjusted life-years burden for the United States is 14.5 million years. Thus, there is considerable public health impetus for understanding the pathophysiology of this disorder.
While symptoms of PMDD correspond with menstrual cycle phase, no differences have been detected in plasma hormone levels—including estradiol, progesterone, luteinizing hormone (LH), and follicle-stimulating hormone (FSH)—between women with PMDD and healthy comparison subjects (
5,
6). We previously reported (
7) that PMDD symptoms remitted during ovarian suppression induced by gonadotropin-releasing hormone (GnRH) agonist in women with PMDD; when physiological levels of ovarian steroids were replaced, typical symptoms recurred in women with PMDD, but no symptoms emerged in women without PMDD under the same conditions. These findings suggest that women with PMDD have an underlying behavioral sensitivity to the normal physiological events of the menstrual cycle. Although the symptoms undoubtedly involve brain mechanisms, the neural substrate of this differential response has yet to be fully characterized.
There is substantial evidence that hormones affect brain function in healthy women. Neuroimaging studies have demonstrated that ovarian steroids modulate activity in brain regions and circuits involved in processes relevant to the symptoms of PMDD, including the functions of the prefrontal cortex, reward systems, and stress neurocircuitry (
8–
11). In contrast, few imaging studies have investigated brain function in PMDD, but some differences from comparison subjects have been observed in both the luteal (symptomatic) and follicular (nonsymptomatic) phases of the menstrual cycle (
12–
14). Similarly, abnormalities of verbal recall have been observed in women with PMDD during both luteal and follicular phases of the menstrual cycle (
15). These findings are of particular interest because some of the neurofunctional abnormalities were observed after menstruation, when women with PMDD are not typically symptomatic (
16). This observation may be consistent with an underlying trait-like neural vulnerability. However, it remains unclear whether trait or state characteristics of brain function are related to PMDD.
In order to clarify the relationship between brain function and disability in PMDD, we took a multimodal neuroimaging approach, using both fMRI and PET together with a cognitive task that accesses dorsolateral prefrontal cortex circuitry, important in both the cognitive and affective components of PMDD, to validate these associations. If the behavioral impairments of PMDD were linked to altered brain function at the trait level, then not only would abnormal task-related neural activation be seen in women with this disorder, but also the magnitude of those changes would be related to the degree of disability and disease burden, as measured independently of hormone manipulation and brain imaging. To test this hypothesis and the proposal that PMDD involves abnormal reactions to normal hormone levels, we used PET and fMRI to compare neural responses of patients with PMDD and healthy comparison subjects undergoing a 6-month pharmacological protocol that allowed us to carefully control the hormonal milieu to which each woman’s brain was exposed. Additionally, we investigated the relationship between activation in affected brain regions during this protocol and illness severity measured before initiation of the protocol.
Method
Participants
The study participants (
Table 1) were women who had regular cycles, were medication free, were not medically ill (as assessed by history, physical examination, neurological examination, MRI, gynecological examination, pap smear within the last year, laboratory tests, and electrocardiogram), and were not pregnant. They were paid for participation according to NIH volunteer guidelines. The study protocol was approved by the NIH Combined Neurosciences Institutional Review Board and Radiation Safety Committee, and all women provided written consent.
Prior to entering the pharmacological protocol, clinical characteristics were obtained by semistructured clinical interview, and all patients with PMDD confirmed the timing and severity of their mood-related symptoms prospectively with daily self-ratings for 3 months using a four-item visual analogue scale (17–19). Average self-ratings of negative moods (i.e., irritability, depression, anxiety, and mood swings) increased at least 30% (relative to the range of the scale used) for PMDD patients in the 7 days before menses relative to the 7 days after menses in at least two of the three baseline cycles. After this screening phase, but before the pharmacological protocol, each woman who met the 30% threshold completed a modification of the daily rating form (20, 21) to confirm that she also met DSM-IV criteria for PMDD. Patients with current axis I psychiatric diagnoses or any diagnoses within the past 2 years as measured by the Structured Clinical Interview for DSM-IV (SCID [
22]) were excluded.
Comparison subjects were recruited through advertisements. They had no history of menstrual-related mood or behavior disturbances, as confirmed during the 2 months before the study with the same daily self-ratings used by the patients, and they had no current or past axis I diagnoses, including alcohol and substance abuse, as confirmed by the SCID.
Hormone Manipulation Protocol
Throughout the 6-month protocol (
Figure 1), participants received monthly injections of the GnRH agonist leuprolide acetate (3.75 mg i.m.). Leuprolide suppresses ovarian function and the secretion of endogenous estradiol and progesterone. For the first 3 months, women received leuprolide only. Following this initial phase, the women entered a 3-month hormone add-back phase while continuing to receive monthly leuprolide injections. Women were randomly assigned to receive 5 weeks of transdermal 17β-estradiol, 0.1 mg/day, or progesterone vaginal suppositories, 200 mg b.i.d., in a double-blind crossover design with a 2-week washout between hormone administration periods. In addition, during the fifth week of estradiol add-back, all women received both estradiol and progesterone to induce menses. Plasma estradiol and progesterone levels were measured at each study visit and before each imaging session (
Table 2).
Rating Scales
Symptom rating forms were completed daily by all women before the study began and during the 6-month leuprolide protocol. The forms included a 12-item visual analogue scale (an extended version of the scale used during the 3-month baseline screening phase) and a modification of the daily rating form (
20,
21), both completed each evening. The daily rating form and the visual analogue scale were used to assess the severity of common symptoms of PMDD, to confirm that each woman with PMDD met DSM-IV criteria for PMDD, and to measure symptom severity in all participants during the hormone manipulation study. Finally, all women completed the Rating Scale for Premenstrual Tension Syndrome (
24) during each scanning session.
As part of the SCID, a Global Assessment of Functioning Scale (GAF) score (
25), measuring overall psychosocial/somatic disturbance and impairment (100=high function, 0=low function/suicidality), was obtained during random points in each woman’s menstrual cycle before the initiation of pharmacological treatment. These scores, obtained independently of the imaging procedures and the pharmacological protocol, were used to examine the relationship between neural activation abnormalities and PMDD-related disability in the patients.
Imaging
The [
15O]H
2O PET method was used to measure rCBF, and fMRI was used to measure the BOLD signal during each of three different hormone conditions: ovarian suppression with leuprolide only, leuprolide plus estradiol replacement, and leuprolide plus progesterone replacement (
Figure 1). Scanning took place during weeks 8–12 of the leuprolide-alone condition and during the third or fourth week of each hormone replacement condition. During both imaging modalities, the n-back working memory test (
26), a cognitive paradigm widely used in neuroimaging, was used to probe prefrontally related activation.
n-Back Task
For the n-back task (see supplemental figure 1 in the data supplement that accompanies the online edition of this article), the numbers 1–4 were displayed on a computer monitor on the points of a diamond shape and were shown in random order (one every 2 seconds). The participants responded to each trial by pressing one of four buttons arrayed on a response box in the same configuration as the stimuli on the monitor. During the 0-back sensorimotor control task, participants pressed the button corresponding to the number shown at the time of the current trial, whereas during the 2-back working memory task, they were instructed to press the button corresponding to the number shown two trials previously.
PET Acquisition and Preprocessing
During each hormone condition, rCBF measurements (10 mCi [15O]H2O/scan) were obtained with a GE Advance three-dimensional scanner (General Electric, Waukesha, Wisc.) during 14 scans of 60-seconds each. Scans from the 2-back and 0-back tasks were alternated and were separated by 6 minutes. The SPM5 software package (London, Wellcome Department of Imaging Neuroscience) was used for image attenuation and correction, reconstruction (32 planes, 6.5 mm full width at half maximum), anatomical normalization to an average template, smoothing (10 mm3 Gaussian kernel), and scaling to remove global blood flow variations. First-level single-subject activation maps (2-back > 0-back) were calculated for each scan session (one statistical map per hormone condition) at a significance threshold of p<0.05.
fMRI Acquisition and Preprocessing
During each hormone condition, participants underwent two runs of the n-back task on a GE 3-T scanner using T2*-weighted gradient echo planar imaging (36 axial slices, 4 mm thickness, 1 mm gap; repetition time=3,000 msec, echo time=35 msec, field of view=24 cm, matrix=64×64). Each run consisted of 14 blocks of 24 seconds each, switching between 2-back and 0-back tasks. Images were preprocessed in SPM5 (slice timing and motion correction, coregistration to a standard template, alignment to the first volume for each subject, and spatial normalization to the Montreal Neurological Institute T1-weighted template). Data were then smoothed with a 10 mm3 Gaussian kernel. First-level single-subject activation maps (2-back > 0-back) were calculated identically to the PET analysis.
Between-Group PET and fMRI Activation Analyses
We used the same approach for the group-level analyses of both the PET and fMRI data. First, for each woman, one first-level activation map per hormone condition was entered as a repeated measure, and diagnosis was entered as a between-group measure in a flexible factorial model in SPM5. Next, for each imaging modality, a binary mask of the main effect of task (2-back > 0-back) across all participants at a significance threshold of p<0.05 (uncorrected) was used to restrict the between-group analyses to regions relevant to task performance. Finally, within these imaging modality-specific masks, between-group differences in working memory activation (2-back > 0-back) were evaluated with a significance threshold of p<0.05 (corrected for false discovery rate). These analyses tested for main-effect between-group differences across all three hormone conditions (i.e., any identified group differences occurring regardless of hormonal state). Additionally, to assess the potential effects of task performance, these analyses were repeated with 2-back accuracy as a covariate of no interest.
Correlations Between GAF Scores and Activation
To test the relationship between activation abnormalities and GAF scores, voxel-wise correlational analyses were performed with data from the PMDD patients. Because comparison subjects had near-ceiling GAF scores with little variation, they were not included in the GAF correlation. For patients, parallel analyses were carried out for the PET and fMRI data. As in the between-group activation analyses, a binary mask of the main effect of task (2-back > 0-back) at p<0.05 (uncorrected) was created for each imaging modality to restrict the correlational analysis to regions relevant to task performance. Next, for each patient, the average activation map (2-back > 0-back) for each hormone condition (three activation maps per patient) was entered into a flexible factorial model with the hormone condition as a repeated measure. Finally, mean-centered GAF scores were entered as contrast weights to produce correlational maps, which were thresholded at p<0.05 (false discovery rate corrected).
Additional Post Hoc Analyses
To determine whether the GAF-activation correlations were consistently observed during each hormone condition, post hoc correlational analyses between the GAF scores and the activation data were carried out separately for each of the three hormone conditions. Using an independently derived dorsolateral prefrontal cortex mask (as cytoarchitechtonically defined in standard stereotaxic space in the postmortem human brain [
27]), we extracted average activation values from a 4-mm sphere surrounding the most robust voxel in the across-hormone GAF-PET and GAF-fMRI correlations (one per woman per hormone condition) and correlated these values with each woman’s GAF score within each hormone condition using Predictive Analytics Software (PASW Statistics 18.0, 2009).
We also investigated the correlation of dorsolateral prefrontal cortex activity with additional dimensions of pre-existing PMDD disease severity, including reported age at onset of PMDD, duration of PMDD, and the average change from pre- to postmenses in visual analogue scale scores for negative affective symptoms obtained during baseline screening (before leuprolide treatment). Finally, to further examine the relationship of the dorsolateral prefrontal cortex to the hormone manipulation-induced symptoms of PMDD, dorsolateral prefrontal cortex activity was correlated with symptom ratings obtained during the hormone manipulation protocol (i.e., Premenstrual Tension Syndrome Scale scores at the time of each scan and the 7-day average of irritability ratings from the daily rating form during the week of each scan).
Results
The patients and comparison subjects did not differ significantly in age, racial distribution, body mass index, handedness, or education for either imaging modality (
Table 1). Five patients with PMDD had a history of major depression, and one patient met criteria for a history of substance abuse disorder (alcohol). The patients’ GAF scores ranged from 51 to 75 (mean=63, SD=6, including two patients with scores >70 who were less severely affected). We observed no differences between GAF scores acquired during luteal and follicular phases (two-sample t test, α=0.05). GAF scores for the comparison subjects ranged from 90 to 100 (mean=92, SD=3). Reported ages at PMDD onset ranged from 11 to 45 years (mean=26.2 years, SD=11.2). The average duration of PMDD was 11 years (range=1–29 years).
Hormone Levels and Behavioral Findings
For both the PET and fMRI cohorts, blood plasma measurements confirmed ovarian suppression by leuprolide and the replacement of the appropriate ovarian steroid during each add-back condition. The hormone levels of patients and comparison subjects did not differ significantly (
Table 2).
Patients experienced a recurrence of typical PMDD symptoms while on both estradiol and progesterone replacement, but they were asymptomatic during the leuprolide-only condition. In contrast, comparison subjects experienced no mood or behavioral symptoms during any of the three hormone conditions.
Task Performance During Imaging
All participants performed well above chance (25%) on all runs across all hormone conditions. There was a main effect of hormone condition on 2-back percent accuracy (number of correct responses divided by the total number of possible responses) during both PET and fMRI scanning. Post hoc comparisons between hormone conditions using the Bonferroni correction revealed that this effect was due to lower scores during the leuprolide-only condition, likely because this testing phase occurred first (although other factors cannot be ruled out). A main effect of diagnosis was observed in the fMRI study (patients performed more poorly than comparison subjects in all three hormone conditions), but not in the PET study; patients’ scores did not differ between PET and fMRI, but comparison subjects scored higher in fMRI compared with PET. We observed no significant diagnosis-by-hormone interaction for either imaging study (
Table 2). Covarying for performance did not change the activation results.
Between-Group Differences in PET and fMRI Activation
The results were remarkably consistent across the two imaging modalities. In both data sets, both groups showed a consistent and robust pattern of activation in regions associated with working memory and cognitive control, including the dorsolateral prefrontal cortex, orbitofrontal cortex, cingulate cortex, and inferior and superior parietal lobules. However, a number of highly significant between-group differences were observed.
In both the PET and fMRI activation data, a main effect of group across all hormone conditions taken together was observed at p<0.05, false discovery rate corrected (
Figure 2A; see also supplemental table 1 online). In both imaging modalities, the patients showed abnormal prefrontal recruitment, specifically greater activation than comparison subjects throughout the dorsolateral prefrontal cortex (Brodmann’s area [BA] 9/46) bilaterally, as well as in the medial frontal gyrus (BA 8) and the cerebellum (
Figure 2A).
In contrast, regions in which patients showed less activation than comparison subjects were sparse. There were no such findings in the PET data, and only small clusters were found in the cuneus (BA 19), precuneus (BA 7), and lateral temporal cortex (BA 39) in the fMRI data. Covarying for performance demonstrated that the between-group differences were not due to task performance.
Correlations Between GAF Scores and Activation
As in the between-group activation analyses, the correlational findings from the patient PET and fMRI data sets were remarkably congruent (
Figure 2B; see also supplemental table 2 online). In both imaging modalities, patients’ 2-back > 0-back activations correlated negatively with GAF scores (the greater the overactivation, the greater the disability as indicated by lower GAF scores) throughout broad swaths of the working memory/executive function pathway bilaterally in regions where this group demonstrated abnormal working memory recruitment: most prominently in the dorsolateral prefrontal cortex (BA 9/46, p<0.005, false discovery rate corrected), as well as in the inferior and superior parietal lobules (BA 40/7). Negative correlations were also observed in the cerebellum in both imaging modalities and were observed in the posterolateral temporal cortex (BA 39) in the fMRI data.
Positive correlations with GAF scores were sparse. There were no such findings in the PET data, and activation correlated positively only in small clusters in the left superior temporal gyrus (BA 22) and the right middle frontal gyrus (BA 10) in the fMRI data.
Anatomical Convergence Between Activation Abnormalities and GAF Activation Correlations
Figure 3 depicts the overlap between maps of the patients’ overactivation and negative correlations between activations and GAF scores. In both the PET and fMRI data, when the correlational and activation findings were viewed together (p<0.05, false discovery rate corrected for both), the patients’ areas of overactivation converged with regions where GAF correlated negatively, indicating that the abnormal neural recruitment in areas involved in working memory is relevant to the degree of clinical disturbance accompanying PMDD.
Figure 3 also depicts post hoc correlational analyses between the GAF scores and the activation data separately for each of the three hormone conditions. For each imaging modality and each hormone condition, data were plotted for a 4-mm sphere around the voxel with the most robust GAF activation relationship within the independently derived dorsolateral prefrontal cortex region of interest. For both the PET and fMRI data sets, these voxels fell within the patients’ overactivated regions. Moreover, the correlations were similarly present in all of the hormone conditions.
Correlations Between Clinical Characteristics and Dorsolateral Prefrontal Cortex Activation
Both age at onset and duration of PMDD correlated with activation in the dorsolateral prefrontal cortex bilaterally (p<0.005, false discovery rate corrected), with earlier age at onset and longer duration associated with greater activations in this region. Neither age at onset nor duration of illness correlated with GAF scores; however, these two measures were highly correlated with each other (p<0.003 in all cases). To disambiguate the relationships between age at onset and duration with the neuroimaging results, we performed a stepwise regression analysis to determine the relative contributions of the two predictor variables in accounting for the dorsolateral prefrontal cortex findings. For both the PET and fMRI data, age at onset entered into the equation first and accounted for 48% and 37% of the variance, respectively. Including symptom duration in the model did not account for any additional variance in dorsolateral prefrontal cortex activation, indicating that age at onset was the primary predictor variable.
Figure 3 depicts the relationship between age at onset and the neuroimaging data, both across- and within-hormone conditions, assessed the same way as the GAF scores.
The correlation between dorsolateral prefrontal cortex activation and the average pre- to postmenses change in negative symptom scores on the visual analogue scale (before study entry) was also significant (p<0.005, false discovery rate corrected); more severe pre- and postmenses changes in symptoms were associated with greater dorsolateral prefrontal cortex activation. Irritability was the predominant symptom that contributed to these correlations. In contrast, no significant correlations between symptoms and dorsolateral prefrontal cortex activation were observed during the hormone manipulation protocol (daily irritability ratings and Premenstrual Tension Syndrome Scale scores).
Discussion
Premenstrual dysphoric disorder is characterized by somatic symptoms and impairment in affective and cognitive processing. We examined neural recruitment in brain circuits relevant to these behaviors, and we tested the relationship between these findings and disability in this disorder. Our results indicate that under controlled hormonal conditions, patients with PMDD display abnormal patterns of activation during working memory. These abnormalities occurred in the absence of significant differences between comparison subjects and patients in plasma hormone levels, thus providing strong support for the presence of a differential neurophysiological response to cognitive challenge in patients with PMDD relative to comparison subjects under equivalent hormone conditions.
Additionally, correlations between the patients’ activations (2-back > 0-back) and several clinical characteristics, measures of disease burden, and prestudy symptoms were observed with both imaging modalities. These correlations occurred within the same circuitry where activation abnormalities were observed, particularly in the dorsolateral prefrontal cortex. The greater the abnormal dorsolateral prefrontal cortex overactivation, the greater the disability and disease burden (as evidenced by lower GAF scores, earlier age at onset, and greater prestudy menses-related change in negative affect). These findings demonstrate the clinical relevance of this neural network in the manifestations of and possibly the biological substrate of risk for PMDD.
Our findings of pathophysiology in PMDD and correlations of that pathophysiology with clinically relevant measures were confirmed with two different imaging platforms that assess related but distinct neurophysiological parameters of brain function: the PET imaging data provide a direct measurement of rCBF, whereas the fMRI data reflect changes in BOLD signal. Although correlation does not permit inference about causality, the fact that our two imaging data sets were highly convergent, both in anatomical distribution and with regard to consistency across hormone conditions, strengthens the assertion that dorsolateral prefrontal cortex function plays an important role in PMDD.
This assertion is consistent with the fact that dorsolateral prefrontal cortex circuitry is crucial to the cognitive, emotional, and social functions affected in PMDD—not only with regard to cognitive control and specific cognitive operations such as working memory, but also as related to insight, judgment, and maintaining appropriate social and emotional behaviors, which are all affected in PMDD. We used GAF scores to quantify the disease burden in our PMDD patients and to measure the severity and impact of impairment across all of these prefrontally dependent domains. Since activation in the regions found to be abnormal in PMDD, notably the dorsolateral prefrontal cortex, correlated with functional disability (GAF) and with measures of disease severity (age at onset and baseline pre- to postmenses symptom changes), both determined independently of the experimental protocol, our results are consistent with the interpretation that abnormalities in the dorsolateral prefrontal cortex represent a neural substrate of vulnerability, which, when most prominent, may confer earlier age at onset or more severe clinical presentation. Since these correlations were present during all three experimental hormone conditions, our findings suggest that these persistent effects, even in the clinically asymptomatic phase of the hormone manipulation paradigm (i.e., leuprolide alone), may represent an enduring “trait-like” predisposition to this hormonally triggered disorder.
In contrast, we observed no correlations between the symptoms present during each of the three hormone conditions—a negative result that may be consistent with the possibility of a persistent neurobiological diathesis that underlies PMDD. However, this lack of correlation must be considered in light of the fact that these symptoms did not arise in response to naturally fluctuating hormones, as occurs in PMDD, but instead arose in the context of controlled hormonal manipulation. Moreover, this observation requires studies with larger sample sizes, and we cannot rule out the possibility that the pathophysiology could be a sequela of a chronic recurrent mood disorder. Since PMDD symptoms correspond to fluctuating hormone levels, suggesting that state characteristics operate in conjunction with more enduring neural vulnerabilities, future research should explore the interactions between individual hormone conditions and PMDD. Since the dorsolateral prefrontal cortex is a known target for gonadal hormones, as demonstrated by a considerable body of research documenting hormonal modulation of dorsolateral prefrontal cortex function in animals and humans (
8,
28–
31), this region and other nodes modulated by the prefrontal cortex may show hormone-specific changes and interactions not observable with our current sample size. Indeed, seminal work by Goldman et al. (
32–
34) demonstrated that the dorsolateral prefrontal cortex in nonhuman primates comes online with a time course similar to that of the pubertal surge in gonadal hormones, a developmental epoch that also may have relevance to the genesis of PMDD and should be carefully studied in this light. In sum, by identifying the prefrontal cortex, specifically the dorsolateral prefrontal cortex, as key in the pathophysiology of PMDD, our research provides insight into potential neural mechanisms contributing to this disorder.