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

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

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 studies 3 . 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. 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). 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. 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 FC 66 , hormone fluctuations were regressed out as covariates of non-interest in the following image processing. 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 mm 3 , 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 mm 3 , 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. 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 studies 68 . 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 mm 3 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. 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 data 69 . However, the neuroscientific interpretation of anti-correlation has been challenged 70 , and global signal regression may cause a negative shift in the distribution of correlations 71 72 ; therefore, we implemented a mask and addressed positive connectivity only 7 . 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 studies 7 73 . 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 brain 74 . 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. 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 previously 75 . Basically, all the reported network assignments were affirmed by this procedure. 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 . 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 FCs 5 7 . 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. Since studies have reported that state anxiety correlate with mPFC-insula FC 61 , 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.

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).

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 segregation 31 . 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 emotion 10 33 . The vmPFC is also involved in the anticipatory coping mechanism 34 . 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 evaluation 14 35 . Of note, the DMN is a key network altered in chronic pain diseases 17 and abnormal resting-state FC of it can be associated with pain intensity and rumination 36 37 38 . 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 pain 22 . 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 PDMs 6 . The OFC is crucial to adaptively cope with physiological challenges, regulate affective states and modulate pain 39 . 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 processes 40 . In addition, many chronic pain conditions are associated with increased DMN-SMN (including SMA) FC 17 20 . 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 DMN 7 , 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 PAG 40 . The functional interplay between the three core networks might reorganize after stress events 41 . 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 network 42 . 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 anxiety 44 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 modulation 22 45 . 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 tasks 46 , including the reappraisal (involving both ECN and DMN 47 ) 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 pain 48 . 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 experience 49 50 . 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 regions 51 . As a neural substrate of the DAN, the aIPS/SPL has an important role in top-down attention regulation 52 , which modulates pain perception 22 . Diverting attention from pain can decrease the activity of the ascending pain system 53 and can be mechanistic for cognitive–behavioral therapy and meditation 54 . 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 control 53 . As a part of the pain matrix 50 , the insula integrates interoception into the subjective social-emotion 56 and has an important role in allocating attention and evaluating context 57 58 . Beside, dACC and mPFC are associated with interoceptive awareness, and play an important role in the appraisal and expression of fear and anxiety 10 . Increased DMN/mPFC-insula FC has been reported in several chronic pain disorders 18 37 . The alteration of mPFC-insula FC is also present in the unpleasant itch condition 59 . 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 perception 38 . This proposition has been corroborated by a study that found that opioids exert their analgesic efficacy by disrupting the vmPFC/rACC-insula FC 60 . State anxiety is positively correlated with DMN-insula FC in healthy subjects 61 . Although pain and anxiety are based on distinguishable neural mechanisms 62 , anxiety is often accompanied by menstrual pain 3 5 6 . 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 PDM 63 . 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 polymorphism 3 . BDNF genotypes may exhibit cognitive and neurobiological differences 64 . 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 manifestations 65 . 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 conditions 15 .