Dynamic Changes of Functional Pain Connectome in Women with Primary Dysmenorrhea.

Primary dysmenorrhea (PDM) is the most prevalent gynecological problem. Many key brain systems are engaged in pain processing. In light of dynamic communication within and between systems (or networks) in shaping pain experience and behavior, the intra-regional functional connectivity (FC) in the hub regions of the systems may be altered and the functional interactions in terms of inter-regional FCs among the networks may be reorganized to cope with the repeated stress of menstrual pain in PDM. Forty-six otherwise healthy PDM subjects and 49 age-matched, healthy female control subjects were enrolled. Intra- and inter-regional FC were assessed using regional homogeneity (ReHo) and ReHo-seeded FC analyses, respectively. PDM women exhibited a trait-related ReHo reduction in the ventromedial prefrontal cortex, part of the default mode network (DMN), during the periovulatory phase. The trait-related hypoconnectivity of DMN-salience network and hyperconnectivity of DMN-executive control network across the menstrual cycle featured a dynamic transition from affective processing of pain salience to cognitive modulation. The altered DMN-sensorimotor network may be an ongoing representation of cumulative menstrual pain. The findings indicate that women with long-term PDM may develop adaptive neuroplasticity and functional reorganization with a network shift from affective processing of salience to the cognitive modulation of pain.

Primary dysmenorrhea (PDM), or menstrual pain without organic causes, is the most prevalent menstrual complaint of young females. Between 40 to 90% of female adolescents have experienced PDM and 10–20% describe their suffering as so severe and distressing that it requires absence from school or work1. Notably, otherwise healthy women with PDM often display a higher prevalence of incidental brain findings2, particularly of normal variants, when compared with control subjects. Additionally, genetic factors, such as the brain-derived neurotrophic factor (BDNF) Val66Met polymorphism, may account for the susceptibility of an individual to PDM3.

We have previously reported that long-term PDM is associated with alterations in brain metabolism and functional connectivity (FC) in pain modulatory systems, as well as in state- and trait-related structures4567. State-related changes are menstrual pain-primed, whereas trait-related changes exist even without symptoms. Notably, pain-related regions, such as the medial prefrontal cortex (mPFC), posterior cingulate cortex and insula, exhibit abnormal functional and structural changes in otherwise healthy women with PDM468. These regions are parts of the pain matrix9 and are associated with cognitive control of emotions, self-referential processing10, contextualizing1112 and pain generation and modulation1314.

Many key brain systems are engaged for menstrual pain processing; many key regions of these systems are also neural substrates of pain matrix that participate substantially in central processing of pain. Recent studies have demonstrated that the FC of large-scale brain systems/networks may be dynamically altered in patients with PDM and other chronic pain disorders and that these alterations can be either adaptive or maladaptive71516. Such altered cross-network connectivity and the imbalance between the systems have been suggested to be a principal pain-related brain abnormality in chronic pain patients17. The default mode network (DMN) plays an important role in pain disorders1518. It not only participates in episodic memory and in intrinsic attention but also modulates the experience of pain1819. The DMN and sensorimotor network (SMN) reorganize in patients with arm amputation that can be with and without phantom pain20. The mPFC, a central node within the DMN, acts as a cardinal coordinator by interacting with the dorsal attention network (DAN), executive control network (ECN) and salience network (SN) during stress and pain adaptation and modulation152122. The ECN balances self-referential awareness with attention to salient information23 and modulates pain anticipation24. The SN, which is related to anxiety25, is activated when attention is maintained on pain15. The DAN provides a direct corticocortical pathway for the attentional modulation of pain22. The DMN, ECN and SN networks constitute the core neural systems in the human brain26. To address the resilient crosstalk between these systems under clinical pain condition, PDM as a genuine condition of chronic pain7 with spontaneous alternation of painful and pain-free periods may serve as a good clinical model of choice. As pain experiences in females are usually more sophisticated than males27 and PDM subjects are vulnerable to anxiety28, it is plausible the FC alterations in PDMs may be manifested as trait-related changes in response to different affective states. Furthermore, PDM may feature functional interactions between large-scale brain networks.

In this study, we implemented a hierarchical approach to study the FC constructs of PDM. We reasoned that FC changes could be manifested as intra-regional (local) and inter-regional (long- distance) FC changes. To study the neuroplasticity of intra-regional FC, we used regional homogeneity (ReHo) analyses to estimate the local synchronization of resting brain activity within a segregated brain region29. An increased local synchrony indicates increased coherence of spontaneous neuronal activity30, which implicates functional integration within a region. Conversely, desynchronized local activity indicates divergent local functions or increased intra-regional functional segregation31. To study inter-regional FC alterations among the aforementioned networks, we performed seed-based FC analyses of large-scale networks in the significant regions that were identified in the ReHo analysis. Incorporating ReHo in the seed-based FC analyses may improve the sensitivity of FC analyses and reduce the uncertainty of seed selection32. We considered the regions of ReHo alterations the key functional hubs for PDM, and the changes of inter-regional FCs as seeded on these ReHo-hubs may serve as critical functional correlates for brain resilience under long-term PDM. Collectively, we hypothesized that the local synchronization within hub region(s) of the key systems, particularly the DMN, would be altered and that the functional interactions among the networks (DMN, SN, ECN, and DAN) would be dynamically and functionally reorganized in patients with PDM to cope with the repeated stress of menstrual pain.

Results

Demographic and psychological data

In brief, 46 PDMs and 49 controls (CON) were enrolled in the current neuroimaging study. There were no significant between-group differences regarding age (PDM: 23.33 ± 2.42 years of age, CON: 23.80 ± 2.47 years of age, p = 0.351), age at menarche (PDM: 11.93 ± 1.32 years of age, CON: 12.31 ± 1.08 years of age, p = 0.14), years of menstruating (PDM: 11.39 ± 2.84 years, CON: 11.49 ± 2.86, p = 0.867), or average duration of one menstrual cycle (PDM: 29.50 ± 1.28 days, CON: 29.52 ± 1.19 days, p = 0.936). For the menstrual pain experience, PDM subjects had in average 9.22 ± 2.85 years of pain history; 65.2% of them had absenteeism from school or works and 41.3% of them had taken analgesics occasionally. The PDM group had been experiencing moderate cyclic menstrual pain (present total scores of McGill Pain Questionnaire-pain rating index: 30.11 ± 11.99, with sensory score: 16.54 ± 6.01, affective score: 3.91 ± 2.85, evaluation score: 2.54 ± 2.02, miscellaneous score: 7.11 ± 3.62 in each category; present pain intensity: 2.72 ± 0.96).

For the psychological inventories, significant group (PDM vs. CON)-phase (menstruation [MENS] vs. periovulatory [POV]) interactions were found in the total (F = 5.921, R = 0.060, p = 0.017) and helplessness scores (F = 8.498, R = 0.084, p = 0.004) of the Pain Catastrophizing Scale (PCS), Beck Anxiety Inventory (BAI; F = 43.802, R = 0.320, p < 0.001), the state anxiety (F = 20.055, R = 0.177, p < 0.001) of Spielberger State-Trait Anxiety Inventory (STAI), and Beck Depression Inventory (BDI; F = 31.910, R = 0.255, p < 0.001). No significant group-phase interactions were found in the rumination (F = 1.569, R = 0.017, p = 0.214) and magnification (F = 2.021, R = 0.021, p = 0.158) scores of the PCS, and trait anxiety of STAI (F = 2.892, R = 0.030, p = 0.092). Significant main effects of group were found in the total (F = 71.504, R = 0.435, p < 0.001), rumination (F = 75.315, R = 0.447, p < 0.001), helplessness (F = 68.698, R = 0.425, p < 0.001) and magnification (F = 31.026, R = 0.250, p < 0.001) scores of the PCS, BAI (F = 57.994, R = 0.384, p < 0.001), the state anxiety (F = 16.437, R = 0.150, p < 0.001) and trait anxiety (F = 8.312, R = 0.083, p = 0.005) of STAI, and BDI (F = 14.514, R = 0.135, p < 0.001). Significant main effects of phase were found in the BAI (F = 27.334, R = 0.227, p < 0.001), state anxiety (F = 23.071, R = 0.199, p < 0.001) of STAI and BDI (F = 20.724, R = 0.182, p < 0.001). No main effects of phase were found in the total (F = 3.276, R = 0.034, p = 0.074), rumination (F = 2.544, R = 0.027, p = 0.114), helplessness (F = 3.713, R = 0.038, p = 0.057) and magnification (F = 0.264, R = 0.003, p = 0.609) scores of the PCS, and the trait anxiety (F = 1.518, R = 0.016, p = 0.221) of STAI (see Supplementary Table S1).

Serum hormone levels

There were no significant between-group differences regarding the levels of estradiol, progesterone or testosterone in MENS phase and POV phase, respectively (see Supplementary Table S2).

Intra-regional connectivity (ReHo)

For the interaction effect, the between-group comparison of phase differences revealed significant greater ReHo alterations in the left declive and right orbitofrontal cortex (OFC), and lesser ReHo alterations in the left aIPS/SPL in the PDM group (Table 1; Fig. 1A). No main effects of group and phase were noted. For the between-group planned contrast, there was no difference during the MENS phase. The bilateral ventromedial prefrontal cortex (vmPFC, BA10) of DMN network, exhibited a trait-related reduction of ReHo during the POV phase in the PDM group vs. CON group (Table 1; Figs 1A and S1). For the between-phase planned contrast, PDM group demonstrated state-related reduction of the ReHo in the anterior part of the left intraparietal sulcus (aIPS, BA40) and left superior parietal lobe (SPL) of DAN network during the MENS phase vs. the POV phase (Table 1; Figs 2A and S1). The CON group demonstrated state-related reduction of ReHo in the left declive during the MENS phase vs. the POV phase.

Table 1

The main effect of group and phase, interaction effect, between-group differences, and within-group differences in regional homogeneity.

Contrast	Region	BA	Size	tscore	Coordinate
					x	y	z
Main effects of group (PDM vs. CON)
PDM > CON	NS
CON > PDM	NS
Main effects of phase (MENS vs. POV)
MENS > POV	NS
POV > MENS	NS
Interaction
PDMΔ(MENS-POV) > CONΔ(MENS-POV)	declive	–	415	4.62	−22	−68	−26
	OFC	10	512	3.94	38	52	−12
CONΔ(MENS-POV) > PDMΔ(MENS-POV)	aIPS/SPL	40	429	4.58	−42	−38	52
Between-group planned contrast
MENS: PDM > CON	NS
MENS: CON > PDM	NS
POV: PDM > CON	NS
POV: CON > PDM	vmPFC	10	377	3.89	6	62	22
				3.6	0	66	12
				3.51	−6	50	4
Between-phase planned contrast
PDM: MENS > POV	NS
PDM: POV > MENS	aIPS/SPL/VPL	40	520	5.11	−42	−36	48
				3.5	−30	−26	44
				3.39	−36	−46	56
CON: MENS > POV	NS
CON: POV > MENS	declive	−	409	5.22	−22	−68	−26

Peak coordinates refer to Montreal Neurological Institute space. Significance was set at the uncorrected voxel level p = 0.005 followed by the FWE-corrected cluster level p = 0.05. BA, Brodmann area; OFC, orbitofrontal cortex; vmPFC, ventromedial prefrontal cortex; VPL, ventral parietal lobe; aIPS, anterior part of the intraparietal sulcus; SPL, superior parietal lobe; MENS, menstruation phase; POV, periovulatory phase; PDM, primary dysmenorrhea group; CON, control group; NS, non-significant.

Figure 1

Statistical maps of ReHo, ReHo-seeded (Lt vmPFC) FC and the interaction effects.

(A) ReHo maps. PDM group shows decreased ReHo of the Lt vmPFC during the POV phase when comparing to the CON group. For the interaction effect, the between-group (PDM vs. CON) comparison of phase differences (MENS vs. POV) reveals significant greater ReHo alterations in the Lt declive and right orbitofrontal cortex (OFC), and lesser ReHo alterations in the Lt anterior intraparietal sulcus and superior parietal lobule. The bar charts show adjusted ReHo values at the peak voxel of the vmPFC. The error bar corresponds to a 90% confidence interval. (B) FC maps. For ReHo-seeded (Lt vmPFC) FC between-group comparisons, hypoconnectivity in the Lt vmPFC-dorsal anterior cingulate cortex FC and hyperconnectivity in the Lt vmPFC-dorsomedial prefrontal cortex FC are noted in the PDM group when compared to the CON group during both phases. Hypoconnectivity in the Lt vmPFC-pregenual anterior cingulate cortex FC and hyperconnectivity in the Lt vmPFC-dorsolateral prefrontal cortex and -ventral parietal lobe FCs are observed in the PDM group during the POV phase. The PDM group (vs. CON group) shows greater between-phase (MENS vs. POV) Lt vmPFC-supplementary motor area (SMA) FC. Red regions are associated with increased FC and greater phase-difference of ReHo or FC in the PDM group. Blue regions are associated with decreased ReHo, decreased FC and lesser phase-difference of ReHo or FC in the PDM group. Significance thresholded at the uncorrected voxel level p = 0.005 followed by the FWE-corrected cluster level p = 0.05. ReHo, regional homogeneity; PDM, primary dysmenorrhea; CON, control; MENS, menstrual phase; POV, periovulatory phase; FC, functional connectivity; vmPFC, ventromedial prefrontal cortex; Lt, left. *Denotes the contrast(s) that are significant in between-group comparisons.

Figure 2

The state-related ReHo and FC alterations in PDM.

(A) ReHo maps. PDM group shows decreased ReHo of the anterior part of the intraparietal sulcus (aIPS) and superior parietal lobe (SPL) during the MENS phase when compared to the POV phase. The bar charts indicate the adjusted ReHo values at the peak voxel of the aIPS/SPL. The error bar corresponds to a 90% confidence interval. (B) FC maps. For ReHo-seeded (aIPS/SPL) FC between-phase comparisons, hypoconnectivity in the aIPS-bilateral posterior insula and -right postcentral gyrus FCs are observed during the MENS phase in the PDM group. Blue regions are associated with either decreased ReHo or decreased FC. Significance thresholded at the uncorrected voxel level p = 0.005 followed by the FWE-corrected cluster level p = 0.05. *Denotes the contrast(s) that are significant in between-phase comparisons. For abbreviations, see the legend of Fig. 1.

Inter-regional connectivity (FC)

vmPFC (DMN network)

For the interaction effect, the PDM group as compared to the CON group exhibited greater left vmPFC (x = −6, y = 50, z = 4)-supplementary motor area (SMA) FC during MENS phase relative to POV phase (Table 2; Fig. 1B). No interaction was found when the right vmPFC (x = 6, y = 62, z = 22) was used as the seed. For the main effects, please refer to the Supplementary Table S3. For the between-group planned contrast, the PDM group as compared to the CON group exhibited hypoconnectivity of the left vmPFC-bilateral dorsal anterior cingulate cortex (dACC; substrate of SN network) circuit and hyperconnectivity of the left vmPFC-left dorsomedial prefrontal cortex (dmPFC; substrate of ECN network) circuit during both phases. In addition, the PDM group as compared to the CON group exhibited hypoconnectivity of the left vmPFC-bilateral pregenual anterior cingulate cortex (pACC; substrate of SN network) circuit and hyperconnectivity of the left vmPFC-left dorsolateral prefrontal cortex (DLPFC; substrate of ECN network) and left ventral parietal lobe (VPL; substrate of ECN network) circuits during the POV phase (Fig. 1B; Table 2; Figs S2, S3). Furthermore, both the right vmPFC-right inferior frontal gyrus (IFG; substrate of SN network) circuit and the right vmPFC-pACC circuit (DMN-DMN network) exhibited hypoconnectivity in the PDM group as compared to CON group during the POV phase (Table 2; Fig. S4). No between-group difference was detected during the MENS phase when the right vmPFC was used as the seed. For the between-phase planned contrast, no difference was found in the PDM group when the left vmPFC was seeded. In the CON group, hypoconnectivity of the left vmPFC-SMA circuit was found during the MENS phase as compared to the POV phase. In the PDM group, hypoconnectivity of the right vmPFC-pons circuit was found during the MENS phase as compared to the POV phase (Table 2). There was no between-phase difference in the CON group when the right vmPFC was seeded.

Table 2

Significant interaction effect, between-group and within-group differences in ReHo-seed (vmPFC) functional connectivity analyses.

Contrast	Region	BA	Size	tscore	Coordinate
					x	y	z
left side ReHo-seed (−6, 50, 4)
Interaction
PDMΔ(MENS-POV) > CONΔ(MENS-POV)	SMA/MCC	6	539	3.95	2	−6	50
Between-group planned contrast
MENS: PDM > CON	SFG	6	388	3.97	−14	24	58
	dmPFC	9	384	3.94	−10	48	28
MENS: CON > PDM	dACC	32	374	4.2	2	20	32
POV: PDM > CON	dmPFC/DLPFC	8	952	4.24	−12	30	40
	VPL/angular	39	372	3.85	−42	−68	48
POV: CON > PDM	pACC	24	823	4.28	0	38	10
					2	22	28#
Between-phase planned contrast
CON: POV > MENS	SMA	6	580	4.46	−12	−12	54
right side ReHo-seed (6, 62, 22)
Between-group planned contrast
POV: CON > PDM	IFG	–	335	5.06	46	20	6
	pACC	32	691	3.99	4	40	8
Between-phase planned contrast
PDM: POV > MENS	pons	–	325	4.58	12	−38	−34

Peak coordinates refer to Montreal Neurological Institute space. Significance thresholded at the uncorrected voxel level p = 0.005 followed by the FWE-corrected cluster level p = 0.05.

SMA, supplementary motor area; MCC, midcingulate cortex; SFG, superior frontal gyrus; dmPFC, dorsomedial prefrontal cortex; dACC, dorsal anterior cingulate cortex; DLPFC, dorsolateral prefrontal cortex; VPL, ventral parietal lobe; pACC, pregenual anterior cingulate cortex; IFG, inferior frontal gyrus (see Table 1 for other abbreviation details) #the second highest peak, almost the same location with the highest peak in the MENS phase.

aIPS/SPL (DAN network)

For the interaction effect, no group-phase interaction was found when the left aIPS/SPL (x = −42, y = −36, z = 48) was seeded. For the main effects, please refer to the Supplementary Table S4. For the between-group planned contrast, the PDM group as compared to the CON group exhibited hyperconnectivity of the aIPS/SPL-MTG circuit and hypoconnectivity of the aIPS/SPL-anterior lobe of the cerebellum circuit during the MENS phase. The PDM group as compared to the CON group exhibited hyperconnectivity of the aIPS/SPL-MTG and -middle occipital gyrus circuits during the POV phase. For the between-phase planned contrast, PDM group exhibited hypoconnectivity of the aIPS-bilateral posterior insula (DAN-pain matrix) and -right postcentral gyrus circuits (DAN-pain matrix) during the MENS phase as compared to the POV phase (Fig. 2B; Table 3). No significant between-phase difference was found in the CON group.

Table 3

Significant between-group and within-group differences in ReHo-seed (aIPS/SPL) functional connectivity analyses.

Contrast	Region	BA	Size	tscore	Coordinate
					x	y	z
Between-group planned contrast
MENS: PDM > CON	MTG	39	323	3.71	40	−60	20
MENS: CON > PDM	cerebellum anterior lobe	31	378	4.34	6	−48	−34
POV: PDM > CON	STG/MTG	22	996	4.54	56	−56	6
	MOG	19	469	4.27	−34	−80	−2
Between-phase planned contrast
PDM: POV > MENS	claustrum	–	351	4.01	−36	−8	0
	supramarginal	40	498	4.35	54	−26	18
	insula	13	434	4.16	42	2	6

Peak coordinates refer to Montreal Neurological Institute space. Significance was set at the uncorrected voxel level p = 0.005 followed by the FWE-corrected cluster level p = 0.05. MTG, middle temporal gyrus; STG, superior temporal gyrus; MOG, middle occipital gyrus.

Correlation between FC and psychological measurements

In the PDM group, the FC of the left vmPFC-bilateral posterior-mid insula (DMN-pain matrix connectivity) and -dACC (DMN-pain matrix connectivity) were more correlated with BAI during the POV phase when compared to the MENS phase (Fig. 3; Table 4). The FCs of the left vmPFC-bilateral midcingulate cortex and -right middle frontal gyrus positively correlated with BAI, and the FC of the left vmPFC-bilateral declive negatively correlated with BAI during the POV phase. The FC of the right vmPFC-right inferior temporal gyrus and -left declive circuit negatively correlated with BAI during the POV phase. The FCs of the aIPS-bilateral rectal gyrus and -right cerebellum connectivity and BAI were significantly more correlated during the POV phase when compared with the MENS phase. The FC of the aIPS-bilateral rectal gyrus positively correlated with BAI and the FC of the aIPS-left precuneus negatively correlated with BAI during the POV phase. For brain-behavioral correlations, there was no significant correlation during the MENS phase (Table 4).

Figure 3

Correlation between FC and anxiety score inverses across phases in PDM.

The Lt vmPFC-Rt insula, -Lt insula and -bilateral dorsal anterior cingulate cortex (dACC) FCs are positively correlated with BAI during the POV phase but negatively correlated with BAI during the MENS phase. Original FC values were used for the correlation to avoid artificial and inadvertent zeroing. Blue regions are associated with lower correlation with BAI during MENS phase. Significance thresholded at the uncorrected voxel level p = 0.005 followed by the FWE-corrected cluster level p = 0.05. BAI, Beck Anxiety Inventory; Lt, left; Rt, right.

Table 4

ReHo-seeded functional connectivity significantly covarying with the Beck Anxiety Index (BAI).

ReHo-seed	Positive correlation							Negative correlation
	Region	BA	Size	tscore	Coordinate			Region	BA	Size	tscore	Coordinate
phase contrast					X	Y	Z					X	Y	Z
left vmPFC (−6, 50, 4)
POV
	MFG	10	388	4.69	28	48	12	declive	–	265	3.99	−12	−68	−22
	MCC	24	354	4.04	2	2	38
MENS > POV
	NS							insula/dACC	–	1292	4.84	36	−10	−2
											4.04	−6	22	26
											3.82	2	36	8
								insula	–	618	4.26	−28	4	14
right vmPFC (6, 62, 22)
POV
	NS							ITG	21	273	4.34	58	−10	−18
								declive	–	431	3.9	−10	−74	−28
aIPS/SPL (−42, −36, 48)
POV
	rectal gyrus	11	536	5.09	−4	28	−26	precuneus/postcentral gyrus	5/7	874	4.76	−6	−50	66
MENS > POV
	NS							rectal gyrus	11	803	5	4	18	−26
								cerebellum	–	351	4.59	50	−70	−36

Peak coordinates refer to Montreal Neurological Institute space. Significance was set at the uncorrected voxel level p = 0.005 followed by the FWE-corrected cluster level p = 0.05.

MFG, middle frontal gyrus; MCC, midcingulate cortex; ITG, inferior temporal gyrus; dACC, dorsal anterior cingulate cortex; NS, non-significant. (see Table 1 for other abbreviation details).

Discussion

Using ReHo and ReHo-seeded FC analyses, the current study demonstrated state- and trait-related intra- and inter-regional FC alterations that collectively may result in subclinical central modulations of pain in otherwise healthy women with PDM. In our study, we observed a decreased ReHo, or a reduced local synchrony, in the vmPFC, which indicates increased intra-regional functional segregations in reaction to long-term menstrual pain. We also measured changes of inter-regional FC (both intrinsic and extrinsic network connectivity) in the PDM group, which may represent network reorganization in the context of strengthened DMN-ECN FC and attenuated DMN-SN FC when compared to the CON group. In addition, functional interaction between aIPS/SPL and pain-related regions (DAN-pain matrix connectivity) might be associated with a top-down regulation of pain transmission. Finally, we observed an inverse correlation between the vmPFC-posterior/mid insula FC (DMN-pain matrix connectivity) and anxiety scores across different phases. This suggests an adaptive neuroplasticity that is related to pain-primed mental state changes.

The vmPFC of PDM group manifested a relative lack of resilience (Fig. 1A) across the menstrual cycle as reflected in the significantly lower ReHo in the PDMs vs. CONs. The functional expression of DMN throughout the menstrual cycle differed between PDM and CON groups. The stereotyped desynchronized local activity (similar ReHo value across the menstrual cycle) in the PDMs indicates sustained and constrained neuroplasticity for divergent local functionalities or enhanced intra-regional functional segregation31. Being a part of the DMN and a functional hub of the brain, the vmPFC integrates and regulates conceptual information, including affective sensory, self-perception, social cognition, and emotion1033. The vmPFC is also involved in the anticipatory coping mechanism34. It has also been suggested that the vmPFC and the linked nucleus accumbens constitute a system that mediates the effects of self-regulation on pain evaluation1435. Of note, the DMN is a key network altered in chronic pain diseases17 and abnormal resting-state FC of it can be associated with pain intensity and rumination363738. The ReHo reduction herein could possibly be attributed to various factors, such as the self-regulation of emotion, behavioral coping, appraisal/re-appraisal of self-relevant social information, and adaptation for pain modulation in relevant networks in response to PDM.

The functional expression and engagement of DAN throughout the menstrual cycle differed between PDM and CON groups, as reflected by the interaction effect of group and phase. The aIPS/SPL services attention shifting and saliently driven sensory discrimination and is a cardinal substrate of the DAN that participates in the top-down attention to pain22. Therefore, the relative lower ReHo of aIPS/SPL during the MENS phase vs. the POV phase in the PDM group (Figs 1A and 2A) might implicate an active disengagement of pain attention as part of salience modulation. The relative higher ReHo of the aIPS/SPL during the POV phase might connote an overt integration of the functionality for salience detection for the forthcoming stressful events.

Greater ReHo (MENS vs. POV) of the declive and the orbitofrontal cortex (OFC) in the PDM as compared to the CON group (Fig. 1A) could be adaptive to long-term menstrual pain. The OFC plays an important role in visceral sensory/nociceptive processing, and is activated during MENS phase of PDMs6. The OFC is crucial to adaptively cope with physiological challenges, regulate affective states and modulate pain39. Of note, we also found greater left vmPFC-SMA FC during MEMS (vs. POV phase) in PDM group as compared to CON group (Fig. 1B). The medial brain motor area has been suggested to be a visceromotor area and has underrecognized modulatory influence on visceral processes40. In addition, many chronic pain conditions are associated with increased DMN-SMN (including SMA) FC1720. The group-phase interaction of vmPFC-SMA FC (DMN-SMN network) indicates that the cross-network communication is altered in the PDM group. Similar to our recent finding of periaqueductal gray matter (PAG) FC changes to the DMN7, the altered DMN-SMN network may be an ongoing representation of cumulative menstrual pain, which might be generated in part by FC dysfunction between medial brain motor regions and both DMN and PAG40.

The functional interplay between the three core networks might reorganize after stress events41. In the present study, we observed reduced connectivity in the vmPFC-dACC circuit (DMN-SN network) in the PDM group. The dACC is a critical neural substrate of the SN network42. Hong et al.43 suggested that females allocate more neural resources to interoceptive awareness. The authors demonstrated hypoconnectivity between the SN- and DMN-related regions in patients with irritable bowel syndrome (a chronic pain condition that is often comorbid with PDM). Such reduced DMN-SN FC can be adaptive for stress and anxiety44 and perhaps a reactive modulation to PDM. Conversely, we found hyperconnectivity in the vmPFC-VPL (DMN-ECN network) and vmPFC-DLPFC circuits (DMN-ECN network) in the PDM group. These two regions of the ECN are associated with pain attention and modulation2245. The dynamic transition from the DMN-SN to DMN-ECN FC in individuals with PDM might implicate a shift from affective processing of pain salience to the cognitive modulation of long-term menstrual pain. DMN engagement is essential to an improved performance in self-related tasks46, including the reappraisal (involving both ECN and DMN47) of potentially harmful inputs that can intervene relatively early in the process of emotion-generation and the recruitment of executive cognitive control processes. This cross-network coordination plausibly implicates a functional plasticity and a coping mechanism underlying a normal protective adaptation in females with PDM as a result of repeated cyclic stress.

The right anterior insula and right inferior frontal gyrus exhibited trait-related decrease of FC with right vmPFC (Fig. S4, Table 2), and are involved in the SN and right-lateralized ventral attention network (VAN), respectively; possibly encoding the prolonged salience of pain48. The hypoconnectivity might connote a possible coping in active inhibition of interoceptive awareness for pain saliency and attention to pain. The hypoconnectivity of the aIPS-right postcentral gyrus and -bilateral posterior-mid insula FCs (DAN-pain matrix connectivity) can also be adaptive in patients with chronic PDM. The postcentral gyrus and posterior-mid insula service the discriminative dimension of pain experience4950. The posterior insula integrates viscerosensory and motor processing, and participate in the interaction between motor cortical areas controlling painful body regions and distant non-motor cortical brain regions51. As a neural substrate of the DAN, the aIPS/SPL has an important role in top-down attention regulation52, which modulates pain perception22. Diverting attention from pain can decrease the activity of the ascending pain system53 and can be mechanistic for cognitive–behavioral therapy and meditation54. Grant et al.55 reported that meditators who are less sensitive to pain feature a decoupled FC between the dACC and the dorsolateral prefrontal cortex. This disengagement from higher-order evaluative processes has an important role in the regulation of pain and emotion as well as in cognitive control53.

As a part of the pain matrix50, the insula integrates interoception into the subjective social-emotion56 and has an important role in allocating attention and evaluating context5758. Beside, dACC and mPFC are associated with interoceptive awareness, and play an important role in the appraisal and expression of fear and anxiety10. Increased DMN/mPFC-insula FC has been reported in several chronic pain disorders1837. The alteration of mPFC-insula FC is also present in the unpleasant itch condition59. Chronic pain and itch may invoke guarding or scratching behaviors and elicit anxiety concomitantly. Hence, the mPFC-insula FC could be a state-specific biomarker for chronic pain perception38. This proposition has been corroborated by a study that found that opioids exert their analgesic efficacy by disrupting the vmPFC/rACC-insula FC60. State anxiety is positively correlated with DMN-insula FC in healthy subjects61. Although pain and anxiety are based on distinguishable neural mechanisms62, anxiety is often accompanied by menstrual pain356. In the present study, a trend of positive correlation between mPFC-posterior-mid insula FC (DMN-pain matrix connectivity) and anxiety was observed during the POV phase. This finding is consistent with the view that neuroticism can be a personality trait in patients with PDM63. The positive correlation shifted to a negative correlation when menstrual pain was present, thereby implicating a possible adaptation to pain-primed anxiety.

Pain is a subjective experience with wide individual differences. One possible explanation is that females with PDM have a propensity for anxiety that can be attributed to genetic factors, such as the BDNF Val66Met polymorphism3. BDNF genotypes may exhibit cognitive and neurobiological differences64. In addition, menstrual pain is accompanied by complex and subtle changes in emotion and cognition. The associated coping strategies and emotional and cognitive regulation vary across individuals. It can be heuristic to explore the interaction between the BDNF Val66Met polymorphism and its neuronal regulation to stress in individuals with PDM.

The neuroimaging is a sensitive measurement that can provide subclinical or preclinical evidence of significant alterations in brain function even without conspicuous behavioral manifestations65. Our results suggest that PDM subjects, after long-term cyclic menstrual pain, may develop neuroplasticity and brain reorganization as manifested by intra- and inter-regional FC changes (intrinsic and extrinsic network connectivity) of neural networks. The interoceptive awareness and switching roles of the SN network for pain saliency can be substituted by cognitive control function of the ECN network. Such resilient transition can be adaptive in long-term PDM subjects and may contribute to the understanding of system dynamics of functional pain connectome in other acute and chronic pain conditions15.

Materials and Methods

Subjects

The subjects of this study were a subset of the participants from our previous genetic association and behavioral studies of PDM who were eligible for neuroimaging studies3. Originally, 326 PDMs and 209 CONs were registered during August 2011 to March 2014 as intended for the integrated behavioral and multimodal imaging genetics studies (magnetic resonance imaging and magnetoencephalography). Two hundred and twenty-nine PDMs and 103 CONs were then excluded after rigorous screening according to the exclusion criteria and one or a combination of the following factors: an irregular menstrual cycle, a prolonged or shortened menstrual cycle, inconsistent pain intensity, pelvic abnormalities by ultrasonography as examined after entering the program, or being unwilling to complete the entire series of genetic, hormonal, behavioral, and multimodal neuroimaging studies (see also Lee et al.3). After entering the actual multimodal neuroimaging experiments, 11 PDMs and 3 CONs were excluded due to incidental brain findings (see also Li et al.2); 36 PDMs and 38 CONs excluded from subsequent analyses due to drop-out that ended in incomplete dataset; 2 PDMs and 10 CONs excluded owing to abnormal hormone level (possibly due to sampling and technical errors); and 2 PDMs and 6 CONs excluded owing to the head motion (>2 mm) or rotation (>2°) during the scan. Eventually, 46 PDMs and 49 CONs were eligible for neuroimaging analyses in this study (whole protocols completed), as also reported in our recently published functional connectivity study of PDM (see Wei et al.7).

The inclusion criteria for subjects with PDM included: 1) 20–30-year-old Taiwanese (Asian) females; 2) a regular menstrual cycle of approximately 27–32 days; 3) a history of menstrual pain longer than 6 months; 4) averaged menstrual pain under regular treatment with a rating at least higher than 4 on a verbal numerical scale (VNS, 0 = not at all, 10 = the worst imaginable pain) in the last 6 months; and 5) right-handedness, as confirmed by the Edinburgh Handedness Inventory. All subjects in the PDM group underwent pelvic ultrasonography to exclude secondary dysmenorrhea that was caused by organic pelvic diseases, such as endometriosis or adenomyosis. All participants were clinically examined and diagnosed in the gynecology clinic by the same gynecologist. The inclusion criteria for the healthy female control subjects were similar to those for the PDM group, with the exception that the members of the control group must have no pain during menses (VNS = 0). The exclusion criteria for all the participants were: 1) using oral contraceptives, hormonal supplements, Chinese herbal medicine, or any central-acting medication (e.g., opioid, anti-epileptics) within 6 months prior to the study; 2) pathological pituitary gland disease; 3) organic pelvic disease; 4) any psychiatric or neurological disorders; 5) brain trauma or brain surgery; 6) immediate plans for pregnancy or a positive pregnancy test; 7) a history of childbirth; and 8) having a metal/pacemaker implant, claustrophobia, or any contraindications to magnetic resonance imaging (MRI). No analgesics were used within 24 hours of the scans. The study was conducted in accordance with the Declaration of Helsinki and was approved by the Institutional Review Board of Taipei Veterans General Hospital. Written informed consents were obtained from all participants.

Experimental design

MRI scanning was individually scheduled according to the commencement day of menstruation for each subject. Psychological assessments, blood samples for gonadal hormone assays and MRI images (T1 and resting-state fMRI images) were obtained at two time points during the menstrual cycle: the MENS phase (days 1–3 of the menstrual cycle) and the POV phase (days 12–16 of the menstrual cycle). Ovulation was confirmed using a urinary luteinizing hormone test (Han Chiun Proper LH Rapid Test).

Psychological measurements

The patients in the PDM group completed the PCS and recounted their overall menstrual pain by using the McGill Pain Questionnaire during the inception interview. Then, during the MENS phase, subjects were assessed for their present experience of menstrual pain using the McGill Pain Questionnaire. All participants completed the BDI, BAI and STAI during the MENS and POV phases to assess their psychological status.

Serum gonadal hormone measurements

The serum that was extracted during both phases was stored for batch analysis using commercialized assays (UniCel DxC 800 Synchron Clinical Systems, Beckman Coulter, Inc., Brea, CA). The total serum concentrations were assayed using a chemiluminescence immunoassay technique for estradiol and progesterone and a radioimmunoassay technique for the testosterone. As gonadal hormones may affect resting-state FC66, hormone fluctuations were regressed out as covariates of non-interest in the following image processing.

Image acquisition

Resting-state functional MRI images were acquired using a 3.0 Tesla MRI scanner (Magnetom Trio Tim, Siemens, Erlangen, Germany) with a 12-channel head coil. High-resolution T1-weighted 3-dimensional structural images using a magnetization-prepared rapid-acquired gradient echo sequence (MPRAGE; [TR]/[TE] = 2530 ms/3.03 ms, flip angle = 70°, field-of-view = 224 × 256 × 192 mm3, in-plane matrix size = 224 × 256 × 192, in-plane resolution = 1 mm) and T2*-weighted gradient echo sequence ([TR]/[TE] = 2500 ms/30 ms, flip angle = 90°, field-of-view =220 × 220 × 136 mm3, in-plane matrix size = 64 × 64 × 40, in-plane resolution = 3.4 mm [round-out], and 204 volumes per run) were conducted to obtain high-resolution anatomical T1 images and functional MRI images. The first four functional scans of each resting-state functional MRI series were discarded for signal saturation and magnetic field stabilization. Participants remained awake during the scan and were instructed to maintain open eyes and relaxed, still heads without thinking about any topic in particular. Head cushions and earplugs were applied to reduce head motion and noise.

Resting-state analysis

Regional homogeneity (ReHo) analysis

Preprocessing was performed by using the Data Processing Assistant for Resting-State fMRI (DPARSF) V2.3 Advanced Edition (State Key Laboratory of Cognitive Neuroscience and Learning, Beijing Normal University, China)67 with statistical parametrical mapping 8 (SPM 8; Wellcome Trust Center for Neuro-imaging, University College London, London, UK) in Matlab 2013b (The MathWorks, Inc., Natick, MA, USA). The preprocessing protocol was adapted from previous studies68. In brief, all functional images were slice-time corrected and realigned for head motion correction. Subjects with head motion of any volume more than 2 mm or 2° were excluded from further processing. Images were normalized using the SPM’s standard echo planar imaging template in Montreal Neurological Institute (MNI) space and re-sampled to an isotropic 2 mm3 voxel size. The resulting time series in each voxel was then linearly detrended and band-pass filtered (0.01–0.08 Hz) to extract the low-frequency oscillations. The following nuisance variables were regressed out: 1) the six head movement parameters computed based on rigid body translation and rotation during the realignment in SPM8; 2) the mean signal within the lateral ventricles for cerebral spinal fluid; and 3) the mean signal within a deep white matter region (centrum ovale). Brain activity was activated in clusters, or contiguous voxels, rather than in a few voxels. ReHo maps were generated by calculating the Kendall’s coefficient of concordance (KCC) of the time series between a given voxel with its nearest neighbors (26 voxels) in a voxel-wise manner. The ReHo maps were spatially smoothed using a 3D Gaussian kernel of 6 mm full-width at half-maximum (FWHM). For standardization purposes, each individual subject’s ReHo map was divided by its own global mean. The peaks of significant clusters were then selected as the ReHobased seeds.

Functional connectivity analysis of ReHo-based seeds

The preprocessing of the functional connectivity analysis protocol was as same as in the ReHo analysis, except that the images were 1) coregistered to individuals’ anatomical image and then normalized to the standard T1 MNI template; 2) spatially smoothed using a 3D Gaussian kernel of 8 mm FWHM; and 3) treated with the global mean signal for additional nuisance variables. We performed global signal regression because it has been shown to maximize the specificity of positive resting-state correlations when using real and simulated data69. However, the neuroscientific interpretation of anti-correlation has been challenged70, and global signal regression may cause a negative shift in the distribution of correlations7172; therefore, we implemented a mask and addressed positive connectivity only7. The seeded regions that were selected for FC analysis were those wherein the ReHo significantly differed in the between-group and between-phase comparisons. The mean time series activity was extracted within the 3mm radius spherical ReHoseeded regions. There is no strict universal standard or any fixed parameter for the seed definition in various modalities of neuroimaging. We had used 3 mm as the radius for the seeds in our published functional studies773. Three-mm was chosen out of the concern to avoid additional partial volume effect since our original fMRI resolution is 3.4 × 3.4 × 3.4 mm and the cortical thickness as reported is about 2 to 4 mm across the brain74. The individual FC maps were computed by the Pearson’s correlation coefficient (r) between the seeds and the related brain regions. After calculating the correlation between the reference time course and the time course of each voxel in the brain, r values were converted into z-values using Fisher’s r-to-z transformation to normalize the distribution. Finally, we only used positive FC maps in subsequent analyses to minimize the possible influences of artifacts of anticorrelations.

Network assignments for the altered ReHo and FC regions

For the network assignments for the significantly altered ReHo and FC regions, we conducted the inverse validation by means of FC analyses, using peaks of significant clusters as the seeds. The FC was considered significant if family-wise error level was p < 0.05 at voxel level. The FC map of each seed was visually compared with the FC map of the key brain systems as published previously75. Basically, all the reported network assignments were affirmed by this procedure.

Statistical analyses

Demographic and psychological data

SPSS Statistics 20.0 (SPSS Inc., Chicago, IL) was used for all analyses. Two-sample t-tests were conducted for the between-group (PDM vs. CON) differences in demographic characteristics and Edinburgh Handedness Inventory scores. To test the group-phase interactions, we used mixed design two-way ANOVA. The significance was thresholded at p = 0.05. For more comprehensive statistical analyses and results, please refer to our previous paper by Lee et al.3.

Image data: ReHo and ReHo-seeded FC

For the interaction effects of group and phase, we compared the phase difference (MENS phase vs. POV phase) of mean ReHo values and ReHo-seeded FCs between the PDM and CON groups by the factorial design module of SPM8. To elucidate the respective trait- and state-related neuroplasticity, we focused on the planned contrasts of between-group differences of respective phase and between-phase differences of respective group rather than the main effects for ReHo values and ReHo-seeded FCs57. The planned between-group and between-phase comparisons of imaging data were conducted using 2-sample and paired t test, respectively. Main effects analyses were also conducted for reference only. Three gonadal hormone levels were treated as non-interest covariates. Significance was thresholded at the uncorrected voxel level p = 0.005 followed by the FWE-corrected cluster level p = 0.05.

Image data: brain-behavior correlation

Since studies have reported that state anxiety correlate with mPFC-insula FC61, we also applied regression analyses between ReHo-seeded FC and BAI for the PDM group. Three gonadal hormone levels were treated as non-interest covariates. Significance was thresholded at the uncorrected voxel level p = 0.005 followed by the FWE-corrected cluster level p = 0.05.

Additional Information

How to cite this article: Wu, T.-H. et al. Dynamic Changes of Functional Pain Connectome in Women with Primary Dysmenorrhea. Sci. Rep. 6, 24543; doi: 10.1038/srep24543 (2016).