Demographic data and general PDM experience
Forty-six subjects were recruited, but 11 were excluded because of abnormal ultrasonography findings, taking hormone medicine, having a chronic pain history, and withdrawal of consent. Among the 35 subjects included, 19 were allocated to the verum acupuncture group and 16 to the sham acupuncture group. During the intervention session, one subject in the verum group withdrew her consent, leaving 18 patients in the verum group and 16 patients in the sham group who completed the intervention (Fig. 1).
Demographic data were not significantly different between the verum and sham groups in terms of age, handedness, age of menarche, gynecologic age, and length of menstrual cycle. No significant intergroup difference was noted in the age of PDM onset, duration of PDM history, duration of menstrual pain per cycle, menstrual pain experience, and menstrual pain intensity (Table 1). All subjects reported that they believed they were receiving the verum acupuncture intervention.
Table 1
Clinical and demographic data of women with primary dysmenorrhea
| | | Verum (N = 18) | Sham (N = 16) | p-Value |
Age (year) | | 24.89 ± 4.59 | 26.13 ± 4.54 | 0.44 |
Headedness (-1 ~ 1) | 0.84 ± 0.16 | 0.91 ± 0.17 | 0.23 |
The age of menarche (year) | 11.72 ± 1.36 | 11.88 ± 1.36 | 0.75 |
Gynecologic age (year) | 13.17 ± 4.78 | 14.25 ± 4.81 | 0.52 |
Length of menstrual cycle (day) | 31.06 ± 2.34 | 30.06 ± 1.95 | 0.19 |
The age of dysmenorrhea onset (year) | 14.50 ± 2.73 | 14.69 ± 2.77 | 0.84 |
Dysmenorrhea history (year) | 10.39 ± 5.80 | 11.44 ± 5.46 | 0.59 |
The duration of menstrual pain per cycle (day) | 1.89 ± 0.68 | 1.63 ± 0.62 | 0.25 |
MPQ | | | |
PRI (0 ~ 78) | 35.83 ± 13.28 | 36.31 ± 12.14 | 0.91 |
PPI (0 ~ 5) | 3.17 ± 1.20 | 3.25 ± 1.00 | 0.83 |
MPQ: McGill Pain Questionnaire; PRI: Pain Rating Index; PPI: Present Pain Intensity. All values are presented as mean ± SD. |
Changes in FC maps across acupuncture sessions
Due to scanner scheduling issue, three and four subjects in the verum and sham groups, respectively, were unable to undergo the rfMRI on week 4. Thus, only 15 and 12 subjects, respectively, were included in the imaging analysis.
After the 8-week acupuncture intervention, in the verum group, increased FC was found in the right middle frontal gyrus, right inferior occipital gyrus, right cuneus, right parahippocampal gyrus, midbrain, and cerebellum, and decreased FC was evident in the right inferior frontal gyrus, right middle temporal gyrus, right anterior cingulate gyrus, left precentral gyrus, left inferior parietal lobule, and left angular gyrus (Fig. 2, Table 3). In the sham group, increased FC was found in the right superior parietal lobule, left temporal gyrus, midbrain, and cerebellum, whereas decreased FC was evident in the right rectal gyrus, right superior frontal gyrus, right inferior frontal gyrus, left inferior temporal gyrus, and bilateral inferior parietal lobule (Fig. 3, Table 3). After the 8-week intervention, the interaction analysis revealed that compared with the sham group, the verum group had a higher increase in FC in the right middle frontal gyrus, right caudate body, and cerebellum and a higher decrease in FC in the left precentral gyrus, bilateral precuneus, and brainstem (Table 3).
Table 3
Changes in functional connectivity map after acupuncture intervention
Week 0 > Week 8 | | Week 0 < Week 8 |
| | | | Coordinate | | | | | | Coordinate |
Anatomical Area | BA | Size | Zmax | x | y | z | | Anatomical Area | BA | Size | Zmax | x | y | z |
Verum | | | | | | | | | | | | | | |
R Ant Cingulate G | 32 | 112 | 3.86 | 18 | 24 | 42 | | Bil Midbrain | | 285 | 4.20 | -4 | -28 | -16 |
L Inf Parietal Lob | 40 | 124 | 3.68 | -44 | -44 | 56 | | | | | 3.21 | 6 | -42 | -6 |
L Precentral G | 4/6 | 347 | 3.56 | -24 | -30 | 70 | | R Cuneus | 18 | 912 | 3.71 | 4 | -64 | 2 |
L Angular G | 39 | 419 | 3.44 | -40 | -78 | 30 | | R Mid Frontal G | 10 | 182 | 3.59 | 40 | 50 | 14 |
L Paracentral Lob | 31 | 47 | 3.37 | -6 | -24 | 46 | | R Inf Occipital G | 18 | 45 | 3.57 | 26 | -92 | -22 |
R Inf Frontal G | 9 | 47 | 3.36 | 64 | 6 | 30 | | Bil Cerebellum Uvula | | 248 | 3.48 | 8 | -76 | -44 |
R Mid Temporal G | 39 | 87 | 3.24 | 42 | -60 | 24 | | | | | 2.86 | -2 | -86 | -34 |
| | | | | | | | L Cerebellum Culmen | | 45 | 3.09 | -34 | -50 | -36 |
| | | | | | | | L Cerebellum Vermis | | 79 | 3.08 | 0 | -78 | -20 |
| | | | | | | | R Parahippocampal G | 30 | 46 | 2.67 | 18 | -44 | 8 |
Sham | | | | | | | | | | | | | | |
R Rectal G | 25 | 63 | 4.12 | 14 | 12 | -22 | | L Midbrain | | 52 | 4.06 | -8 | -26 | -16 |
| 11 | 68 | 3.35 | 14 | 30 | -22 | | L Parahippocampal G | | 53 | 3.50 | -26 | -22 | -14 |
R Sup Frontal G | 8 | 70 | 3.45 | 10 | 26 | 54 | | R Cerebellum Vermis | | 142 | 3.34 | 0 | -50 | -36 |
Bil Inf Parietal Lob | 40 | 147 | 3.30 | -52 | -48 | 36 | | R Sup Parietal Lob | 7 | 60 | 3.29 | 34 | -82 | 44 |
| | 43 | 3.03 | 54 | -42 | 28 | | L Mid Temporal G | 21 | 60 | 3.09 | -44 | -6 | -20 |
L Inf Temporal G | 21 | 148 | 3.26 | -66 | -28 | -20 | | L Cerebellar Tonsil | | 47 | 3.06 | -12 | -46 | -46 |
R Inf Frontal G | 47 | 121 | 3.19 | 50 | 38 | 0 | | | | | | | | |
| 44 | 77 | 3.04 | 60 | 18 | 16 | | | | | | | | |
Verum > Sham | | | | | | | | | | | | | | |
R Brainstem | | 43 | 3.67 | 20 | -22 | -50 | | Bil Cerebellum Culmen | | 81 | 3.77 | -4 | -36 | -6 |
Bil Precuneus | 19 | 112 | 3.64 | 34 | -82 | 44 | | | | | 3.65 | 4 | -40 | -6 |
| 19 | 102 | 3.21 | -34 | -76 | 34 | | R Mid Frontal G | 46/10 | 322 | 3.52 | 44 | 48 | 10 |
L Precentral G | 4/6 | 168 | 3.33 | -24 | -30 | 68 | | R Caudate Body | | 60 | 3.50 | 12 | -6 | 24 |
BA: Brodmann area; Size: number of voxels in the cluster; Zmax: peak Z value; L: left; R: right; Bil: bilateral; Sup: superior; Inf: inferior; Mid: middle; Med: medial; Ant: anterior; G: gyrus; Lob: lobule; N: nucleus.
After the early stage of intervention, in the verum group, increased FC was evident in the left medial frontal gyrus, left superior temporal gyrus, left parahippocampal gyrus, and cerebellum, and decreased FC was evident in the right superior frontal gyrus, left precentral gyrus, left paracentral lobule, and bilateral precuneus. In the sham group, increased FC was identified in the right middle occipital gyrus, left precentral gyrus, and cerebellum, whereas decreased FC was apparent in the right medial frontal gyrus, right superior temporal gyrus, right middle temporal gyrus, and cerebellum. The interaction analysis revealed that compared with the sham group, the verum group had a higher increase in FC in the right middle temporal gyrus, left superior temporal gyrus, and cerebellum and a higher decrease in FC in the right middle occipital gyrus and left superior frontal gyrus (Table S1).
After the late stage of intervention, in the verum group, increased FC was apparent in the left middle frontal gyrus, left middle occipital gyrus, left precuneus, bilateral cuneus, and bilateral parahippocampal gyrus, whereas decreased FC was evident in the right inferior parietal lobule, left medial frontal gyrus, left supramarginal gyrus, and bilateral middle temporal gyrus. In the sham group, increased FC was identified in the right inferior parietal lobule, left precentral gyrus, left superior temporal gyrus, bilateral posterior cingulate gyrus, and cerebellum, and decreased FC was apparent in the right middle frontal gyrus, right inferior frontal gyrus, left inferior parietal lobule, left fusiform gyrus, bilateral superior frontal gyrus, and bilateral anterior cingulate gyrus. The interaction analysis revealed that compared with the sham group, the verum group had a higher increase in FC in the left paracentral lobule, left fusiform gyrus, left caudate body, and bilateral superior frontal gyrus and a higher decrease in FC in the right interior parietal lobule, left precentral gyrus, left middle temporal gyrus, and cerebellum (Table S2).
Correlation between menstrual pain experience and the FC of PAG
In the verum group, menstrual pain experience was positively correlated with the FC of PAG in the right middle frontal gyrus, right subcallosal gyrus, left precentral gyrus, and bilateral caudate head and negatively correlated with the right inferior temporal gyrus, right lingual gyrus, left postcentral gyrus, and left superior parietal lobule (Fig. 4, Table 4). In the sham group, menstrual pain experience was positively correlated with the FC of PAG in the right inferior frontal gyrus, left middle temporal gyrus, left fusiform gyrus, left cingulate gyrus, bilateral cuneus, and cerebellum and negatively correlated with the right inferior temporal gyrus, left caudate body, bilateral thalamus, and cerebellum (Fig. 4, Table 4). The interaction analysis revealed that compared with the sham group, the menstrual pain experience in the verum group was more positively correlated with the FC of PAG in the right middle frontal gyrus, right postcentral gyrus, right inferior temporal gyrus, right caudate head, left precentral gyrus, and left middle temporal gyrus and more negatively correlated in the left cuneus, left anterior cingulate gyrus, and cerebellum (Table 4).
Table 4
Correlation between menstrual pain experience and functional connectivity of PAG
Positive | | Negative |
| | | | Coordinate | | | | | | Coordinate |
Anatomical Area | BA | Size | Zmax | x | y | z | | Anatomical Area | BA | Size | Zmax | x | y | z |
Verum | | | | | | | | | | | | | | |
Bil Caudate Head | | 463 | 3.92 | 12 | 18 | 6 | | L Postcentral G | 2 | 59 | 4.01 | -30 | -42 | 72 |
| | 44 | 2.83 | -14 | 18 | 4 | | R Inf Temporal G | 20 | 90 | 3.84 | 58 | -28 | -28 |
R Subcallosal G | 13 | 463 | 3.32 | 18 | 16 | -14 | | L Sup Parietal Lob | 7 | 79 | 3.51 | -40 | -60 | 60 |
R Mid Frontal G | 6 | 137 | 3.28 | 34 | 6 | 54 | | R Lingual G | 18 | 87 | 3.28 | 12 | -78 | -14 |
L Precentral G | 6 | 88 | 3.18 | -26 | -20 | 54 | | | | | | | | |
Sham | | | | | | | | | | | | | | |
Bil Cuneus | 19 | 927 | 4.65 | 32 | -86 | 20 | | R Inf Temporal G | 20 | 89 | 4.47 | 36 | -8 | -44 |
| 19 | 635 | 3.68 | -16 | -98 | 20 | | R Cerebellum Culmen | | 75 | 3.37 | 8 | -60 | -26 |
L Mid Temporal G | 21 | 62 | 3.26 | -54 | 8 | -16 | | L Caudate Body | | 69 | 3.21 | -20 | -6 | 24 |
L Ant Cingulate G | 32 | 61 | 3.26 | -16 | 14 | 32 | | Bil Thalamus | | 73 | 2.99 | -4 | -2 | 4 |
L Fusiform G | 37 | 74 | 3.20 | -38 | -62 | -16 | | | | | 2.82 | 2 | -8 | 4 |
L Cerebellum Lob | | 48 | 3.08 | -4 | -68 | -52 | | | | | | | | |
R Inf Frontal G | 47 | 43 | 3.06 | 28 | 12 | -22 | | | | | | | | |
Verum > Sham | | | | | | | | | | | | | | |
R Inf Temporal G | 20 | 64 | 3.85 | 36 | -6 | -40 | | Bil Cuneus | 19 | 432 | 3.85 | -26 | -88 | 34 |
R Mid Frontal G | 6 | 212 | 3.51 | 34 | 6 | 56 | | | 19 | 356 | 3.72 | 34 | -90 | 18 |
R Caudate Head | | 187 | 3.48 | 12 | 18 | 6 | | R Cerebellum Tuber | | 76 | 3.68 | 30 | -90 | -38 |
L Precentral G | 6 | 87 | 3.37 | -26 | -20 | 56 | | L Cerebellum Lob | | 91 | 3.21 | -26 | -74 | -58 |
L Mid Temporal G | 21 | 62 | 3.21 | -36 | 6 | -38 | | L Cingulate G | 24 | 52 | 3.13 | -14 | 14 | 34 |
R Postcentral G | 3 | 46 | 3.00 | 22 | -30 | 62 | | | | | | | | |
BA: Brodmann area; Size: number of voxels in the cluster; Zmax: peak Z value; L: left; R: right; Bil: bilateral; Sup: superior; Mid: middle; Inf: inferior; Ant: anterior; G: gyrus; Lob: lobule.