- Research
- Open access
- Published:
Differential brainstem connectivity according to sex and menopausal status in healthy male and female individuals
Biology of Sex Differences volume 16, Article number: 25 (2025)
Abstract
Background
Brainstem nuclei play a critical role in both ascending monoaminergic modulation of cortical function and arousal, and in descending bulbospinal pain modulation. Even though sex-related differences in the function of both systems have been reported in animal models, a complete understanding of sex differences, as well as menopausal effects, in brainstem connectivity in humans is lacking. This study evaluated resting-state connectivity of the dorsal raphe nucleus, right and left locus coeruleus complex (LCC), and periaqueductal gray (PAG) according to sex and menopausal status in healthy individuals. In addition, relationships between systemic estrogen levels and brainstem-network connectivity were examined in a subset of participants.
Methods
Resting-state fMRI was performed in 47 healthy male (age, 31.2 ± 8.0 years), 53 healthy premenopausal female (age, 24.7 ± 7.3 years; 22 in the follicular phase, 31 in the luteal phase), and 20 postmenopausal female participants (age, 54.6 ± 7.2 years). Permutation Analysis of Linear Models (5000 permutations) was used to evaluate differences in brainstem-network connectivity according to sex and menopausal status, controlling for age. In 10 males and 17 females (9 premenopausal; 8 postmenopausal), estrogen and estrogen metabolite levels in plasma and stool were determined by liquid chromatography-mass spectrometry/mass spectrometry. Relationships between estrogen levels and brainstem-network connectivity were evaluated by partial least squares analysis.
Results
Left LCC-executive control network connectivity showed an overall sex difference (p = 0.02), with higher connectivity in females than in males; however, this was mainly due to differences between males and premenopausal females (p = 0.008). Additional sex differences were dependent on menopausal status: PAG-default mode network (DMN) connectivity was higher in postmenopausal females than in males (p = 0.04), and PAG-sensorimotor network (SMN) connectivity was higher in premenopausal females than in males (p = 0.03) and postmenopausal females (p = 0.007). Notably, higher free 2-hydroxyestrone levels in stool were reliably associated with higher PAG-SMN and PAG-DMN connectivity in premenopausal females (p < 0.01).
Conclusions
Healthy females show higher brainstem-network connectivity involved in cognitive control, sensorimotor function, and self-relevant processes than males, dependent on their menopausal status. Further, 2-hydroxyestrone, implicated in pain, may modulate PAG connectivity in premenopausal females. These findings may relate to differential vulnerabilities to chronic stress-sensitive disorders at different life stages.
Plain English summary
Males and females show differences in the likelihood to develop chronic pain disorders, as well as anxiety and depression, which often accompany chronic pain, across the lifespan. Brainstem nuclei are small structures in the brain yet show powerful influences across the entire brain. Here, we investigated differences in the connectivity of several brainstem nuclei with major brain networks involved in stress responsiveness and pain modulation between healthy males and healthy premenopausal and postmenopausal females, as well as their association with estrogen levels in the body. We found higher connectivity between the locus coeruleus and the executive control network in females, especially premenopausal females, than in males, which may relate to differences in preferred information processing strategies and anxiety symptoms. We also found higher connectivity between the periaqueductal gray and the default mode network in postmenopausal females and the sensorimotor network in premenopausal females than in males. Further, higher levels of 2-hydroxyestrone, an estrogen metabolite that plays a role in visceral pain sensitivity, was related to periaqueductal gray connectivity in premenopausal females. These latter findings may relate to differences in vulnerability to visceral pain disorders at different life stages. These results help in identifying potential contributors to sex differences in vulnerability to chronic disorders that reduce the quality of life, before a disorder develops.
Highlights
-
We evaluated differences in the connectivity of several brainstem nuclei with major brain networks between healthy premenopausal and postmenopausal females and males, as well as relationships between connectivity and estrogen levels in plasma and stool.
-
Premenopausal females showed (1) higher connectivity between the left locus coeruleus complex and executive control network compared to that in males, and (2) higher connectivity between the periaqueductal gray and sensorimotor network compared to that in males and postmenopausal females.
-
Postmenopausal females showed higher connectivity between the periaqueductal gray and default mode network compared to that in males.
-
Levels of free 2-hydroxyestrone, an estrogen metabolite implicated in visceral pain, in plasma and stool were associated with periaqueductal gray connectivity in premenopausal females.
-
These results may relate to differential vulnerabilities to chronic pain and stress-sensitive disorders across the lifespan.
Background
Brainstem nuclei, including the locus coeruleus complex (LCC), dorsal raphe nucleus (DRN), and periaqueductal gray (PAG), play critical roles in ascending monoaminergic modulation of brain and vital functions, as well as in endogenous descending pain modulation. Alterations in these modulations have been demonstrated in chronic pain disorders, as well as in anxiety and depression, which are often comorbid with each other [1,2,3,4,5]. The LCC is the primary source of noradrenergic innervation of the forebrain and exerts a powerful modulatory role over cognitive and affective functions via widespread cortical and subcortical projections [6]. The DRN is the primary serotonergic nucleus in the central nervous system and modulates several vital functions, including mood, appetite, and sleep, through ascending projections to many cortical and subcortical brain regions [7]. The PAG is involved in integrated descending modulation of pain, as well as in autonomic and behavioral responses to threat. It has reciprocal connections with prefrontal and emotion-regulation regions and receives top-down input from the orbitofrontal cortex and insula [4].
Sex differences in LCC, DRN, and PAG structure and function exist [8, 9]. For instance, animal studies suggest that the LCC arousal system is more sensitive to corticotropin-releasing factor, which is involved in stress, in females than in males [10, 11]. In addition, neuroimaging studies suggest that the LCC has stronger connectivity with the hippocampus in males than in females [12], consistent with animal research demonstrating greater LCC noradrenergic input to the hippocampus in males than in females [13]. Animal research indicates PAG/DRN involvement in sex differences in pain-related behaviors, promoting anti-nociception in males and pain-related locomotor behaviors in females [9]. Females have greater risk for numerous chronic pain disorders, as well as greater risk for anxiety and depression [14]. Sex differences in the connectivity of brainstem nuclei may be related to differential vulnerability to chronic pain conditions and comorbid mood disorders.
Estrogens are known to impact LCC, DRN, and PAG function. Estrogens modulate LCC output, generally increasing noradrenergic levels in target regions [15]. Estrogens also increase the expression of tryptophan hydroxylase, a rate-limiting enzyme for serotonergic synthesis, in the DRN, reducing anxiety and increasing active coping behavior in animal studies [16,17,18]. The PAG contains a large population of estrogen receptors, contributing to known sex differences in response to morphine administration, with a greater antinociceptive effect in males than in females [19].
Menopause is associated with increased risk for anxiety and depression [20], as well as for some chronic pain conditions, such as fibromyalgia, migraine, and back pain [21]. However, the role of estrogen in the risk for and severity of symptoms in disorders of gut brain interaction (DGBI) in premenopausal and postmenopausal females is incompletely understood. We previously found that postmenopausal females with irritable bowel syndrome (IBS) have greater overall IBS symptom severity and worse health-related quality of life [22]. However, in a large global epidemiology study in individuals with IBS and other DGBI, premenopausal females reported a greater frequency of gastrointestinal symptoms compared to males and postmenopausal females [23]. In addition, another study found that menopause may be associated with less risk for DGBI such as IBS [24]. Changes in brainstem connectivity during menopause may contribute to these changes in symptom severity in postmenopausal females. However, a complete understanding of sex and menopausal effects on brainstem connectivity, especially in humans, is lacking.
In the current study, we aimed to evaluate sex and menopausal status effects on resting-state connectivity of the DRN, left and right LCC, and PAG with major brain networks involved in stress responsiveness, pain modulation, and emotion regulation, including the central autonomic network (CAN), default mode network (DMN), emotional arousal network (EAN), executive control network (ECN), salience network (SAL), and sensorimotor network (SMN), in healthy individuals, as well as relationships between affected connectivity and anxiety and somatic symptoms. In addition, relationships between estrogen levels and brainstem-network connectivity were examined in a subset of participants.
Methods
Participants
Healthy male and female participants were recruited from the Los Angeles area through advertisements and local clinics. Exclusion criteria were as follows: chronic pain disorder, major neurological condition or vascular disease, current or past psychiatric disorder, substance use disorder, use of centrally acting medications, pregnant or breastfeeding, weight > 400 lbs, and MRI contraindications (e.g. metal implants). In addition, individuals with < 8 min of low-motion resting-state fMRI data (with low motion defined as framewise displacement < 0.2 mm) were excluded.
All participants underwent a medical history and physical examination. Participants also underwent the Mini-International Neuropsychiatric Interview to assess past/current psychiatric disorders [25]. Sex was self-reported as the sex assigned at birth; prior to October 2023, male and female were the only options presented to the participants, after which additional options including intersex and ‘none of these’ were presented. For consistency, only those indicating male or female at birth were included. Menopausal status was determined by the following criteria: females who had regular menses within the previous 12 months were premenopausal; females who did not have menses within the previous 12 months, with decreased serum estradiol (E2) and increased serum follicle-stimulating hormone levels based on normal laboratory values were considered postmenopausal, in accordance with STRAW + 10 criteria [26]. Among premenopausal female participants, menstrual cycle phase was determined using urine ovulation kits. Pregnancy or childbirth within the past 12 months were exclusion criteria for all female participants.
This study was approved by the Institutional Review Board at the University of California, Los Angeles’s Office of Protection for Research Subjects (Nos. 20–000540 and 20–000515). All participants provided written informed consent.
Questionnaires
Participants completed questionnaires on anxiety and somatic symptoms to enable the examination of the behavioral correlates of implicated brainstem connectivity. These included the Hospital Anxiety and Depression (HAD) scale, as a measure of current anxiety symptoms (total score: 0–21) [27]; State-Trait Anxiety Inventory (STAI), as a measure of trait anxiety (total score: 20–80) [28]; Perceived Stress Scale (PSS), as a measure of ongoing stress burden (total score: 0–40) [29]; and the Patient Health Questionnaire-15 (PHQ-15; total score: 0–15) [30] and Pennebaker Inventory of Limbic Languidness (PILL; total score: 0–54) as measures of the bothersomeness and frequency of common somatic symptoms, respectively [31].
Imaging acquisition and preprocessing
Participants underwent neuroimaging on a 3.0 T Siemens Prisma MRI scanner (Siemens, Erlangen, Germany). Structural MRI (T1-weighted and T2-weighted) and resting-state fMRI scans were performed in accordance with Human Connectome Project (HCP) protocols (version 4.3), with a field of view optimized for the brainstem. Specifically, acquisition parameters for high-resolution T1-weighted scans were as follows: echo time, 1.81 ms; repetition time, 2500 ms; slice thickness, 0.8 mm; number of slices, 208; voxel matrix, 320 × 300; and voxel size, 1.0 × 1.0 × 0.8 mm. Parameters for the T2-weighted scans were as follows: echo time, 564 ms; repetition time, 3200 ms; slice thickness, 0.8 mm; number of slices, 208; voxel matrix, 320 × 300; and voxel size, 1.0 × 1.0 × 0.8 mm. Resting-state scans were obtained in anterior–posterior and posterior-anterior directions (8 min each; total duration, 16 min), with eyes open. Parameters for the resting-state scan were as follows: echo time, 37 ms; repetition time, 1000 ms; flip angle, 52 deg; slice thickness, 2 mm; number of slices, 88; voxel matrix, 104 × 104; and voxel size, 2.0 × 2.0 × 2.0 mm. Spin echo fieldmaps were also acquired in anterior–posterior and posterior-anterior directions for distortion correction.
Images were preprocessed using the ABCD-HCP pipeline (https://github.com/DCAN-Labs/abcd-hcp-pipeline) [32], which is based on the HCP pipeline [33]. Briefly, structural images underwent bias field correction, volume segmentation and cortical surface reconstruction using FreeSurfer 6.0 [34], and MNI registration; resting-state images underwent distortion and bias field correction, realignment, MNI registration via T1-weighted registration, and intensity normalization. Additionally, resting-state images were denoised by regressing out head motion parameters and white matter and cerebrospinal fluid signals, which were refined using respiratory motion filtering (13.2–18.6 breaths per minute), and band-pass filtering (0.008, 0.09 Hz) [35].
Denoised resting-state images were parcellated using the Destrieux atlas for cortical regions [36], Harvard–Oxford atlas for subcortical regions, and ascending arousal network atlas for brainstem regions [37]. Fisher-transformed connectivity matrices were created using 8 min of low-movement data (defined as framewise displacement < 0.20 mm). The correlation between each brainstem region of interest (DRN, left LCC, right LCC, and PAG) and each network of interest (CAN, DMN, EAN, ECN, SAL, and SMN) was calculated as the average of the pairwise correlations between the brainstem region and all of regions belonging to the network. Regions included in each network of interest are indicated in Table 1.
Determination of free and total estrogens in plasma and stool
In a subset of participants (10 males; 9 premenopausal females, with 4 and 5 in follicular and luteal phases, respectively, during scanning and sample collection; and 8 postmenopausal females), the free and total (total = free + conjugated) levels of 15 estrogens and estrogen metabolites in plasma and stool were determined by liquid chromatography-mass spectrometry/mass spectrometry (LC–MS/MS). Estrogen species analyzed included the three main types of estrogen, estrone (E1), estradiol (E2), estriol (E3), as well as 12 additional estrogen metabolites: 2-hydroxyestrone (2OHE1), 2-methoxyestrone (2MeOE1), 2-hydroxyestradiol (2OHE2), 2-methoxyestradiol (2MeOE2), 2-hydroxyestrone-3-methyl ether (3MeOE1), 4-hydroxyestrone (4OHE1), 4-methoxyestrone (4MeOE1), 4-methoxyestradiol (4MeOE2), 16α-hydroxyestrone (16aOHE1), 17-epiestriol (17epiE3), 16-ketoestradiol (16ketoE2), and 16-epiestriol (16epiE3). Both sample preparation and LC–MS/MS methods followed previously published methods [38,39,40], with minor modifications.
Plasma sample preparation
Plasma was prepared with and without enzymatic hydrolysis. For the determination of free estrogen, 100 uL plasma was combined with stable heavy isotope internal standards and 300 uL of basic reaction buffer consisting of 0.15 M acetate buffer pH 4.5 containing 1.0 mg/mL L-ascorbic acid. For the determination of total estrogen, 100 uL plasma was prepared as above, but 15 uL of β-glucuronidase/aryl-sulfatase (from Helix pomatia) was also added and samples were kept gently rocking at 37 °C overnight. For both free and total samples, isopropanol was added (5% v/v) and lipids were extracted with 400 cc supported liquid extraction cartridges (Biotage LLC) using 3 × 1.5 mL dichloromethane extractions and dried under argon.
Stool sample preparation
Stool for both free and total estrogen (100 mg each) was combined with internal standards and 1 mL basic reaction buffer and subjected to bead mill homogenization (Biotage Lysera). Stool for total estrogen was combined with 15 uL deconjugation enzyme and reacted overnight as above. Both sample sets were combined with 500 ul acetonitrile; re-homogenized; and centrifuged. Supernatant was loaded into 2 cc supported liquid extraction cartridges (Biotage) and extracted with 3 × 2 mL DCM followed by drying under argon.
Sample derivatization
Dried plasma and stool samples were combined with 100 uL of 1 mg/mL-acetone dansyl chloride and 100 uL 0.1 M sodium bicarbonate buffer pH 9.2; reacted for 20 min at 60 °C; centrifuged; and transferred to LCMS vials.
LC–MS/MS
LC–MS/MS analysis was performed on an Agilent 1290/SCIEX QTrap 5500 system using principles of tuning, method validation, quality control, calibration, and quantitation that we previously described [41]. MS/MS settings for derivatized estrogens and stable heavy isotope internal standards were determined using reference calibration and internal standards (Steraloids); values were comparable to those previously reported [38, 39]. Liquid chromatography was performed using a Kinetex C18 1.7 um 2.1 × 150 mm column (Phenomenex), while the gradient transitioned from 90% water/0.1% formic acid to 95% acetonitrile/0.1% formic acid across 20 min. Concentration of all estrogens was determined as ng/mL-plasma or ng/5 mg-stool.
Statistical analysis
Group differences in characteristics were evaluated by analyses of variance, with the exception of PHQ-15 and PILL scores. Because of strongly skewed distributions, PHQ-15 and PILL scores were dichotomized by a median split, such that no/minimal somatic burden was defined as a score of 0 or 1 on the PHQ-15/PILL and more than minimal somatic burden was defined as a score of 2 or more PHQ-15/PILL. Group differences in the frequency of more than minimal somatic burden were evaluated by chi-squared analysis.
Permutation analysis of linear models with non-parametric combination (NPC) was used to evaluate brainstem-network connectivity according to sex and menopausal status, with 5000 permutations [42, 43]. NPC combines the test statistics of separate analyses into a single joint statistic, the significance of which is assessed through synchronized permutations for each of the separate tests [43, 44]. In this case, the test statistics of separate analyses of males vs females in the follicular phase, males vs. females in the luteal phase and males vs postmenopausal females, controlling for age, were combined into a joint statistic representing male vs female participants across menstrual status/phase categories (i.e., an overall sex difference) and joint statistics representing differences according to menopausal status (i.e., premenopausal females vs males, premenopausal females vs postmenopausal females, and postmenopausal females vs males). Familywise error (FWE)-corrected p-values < 0.05 were considered significant.
Partial least squares correlation (PLSC) analysis was applied to examine relationships between estrogen/behavioral data and brainstem-network connectivity with significant sex/menopausal status effects, using plscmd in Matlab (http://www.rotman-baycrest.on.ca/pls) [45]. PLSC analysis is a multivariate analytical technique that identifies weighted patterns of variables in two blocks of variables that maximally covary with each other (i.e. latent variables) and is appropriate for data with multicollinearity. As metabolite pathways are interrelated, multicollinearity in the estrogen data was expected. In the present study, Block 1 comprised brainstem-network data and Block 2 comprised estrogen data (free and total levels in plasma and stool) or behavioral data (HAD anxiety, STAI trait anxiety, PSS, PHQ-15, and PILL). Permutation testing (5000 permutations) was used to assess latent variable significance and bootstrap estimation was applied (5000 samples) to assess the reliability of individual saliences within the latent variable. In the present study, we report latent variables with any statistically reliable saliences according to the bootstrap ratio (BSR). For behavioral analyses, BSRs of magnitude 1.96 or greater (corresponding to p < 0.05) were considered statistically reliable. For estrogen analyses, given the limited sample size, a more stringent cutoff was adopted; BSRs of magnitude 2.58 or greater (corresponding to p < 0.01) were considered statistically reliable.
Results
Participant characteristics
From the 50 male and 75 female participants (21 postmenopausal), 3 males and 2 females (1 postmenopausal) were excluded from analysis due to insufficient low-motion resting-state data. Accordingly, 47 healthy male participants (mean age: 31.2 ± 8.0 years), 53 healthy premenopausal female participants (mean age: 24.7 ± 7.3 years), and 20 healthy postmenopausal female participants (mean age: 54.6 ± 7.2 years) were included in the analysis. Among the premenopausal female participants, 22 were scanned during the follicular phase and 31 were scanned during the luteal phase of their menstrual cycle. None of the postmenopausal female participants were taking hormone replacement therapy at the time of scanning. In addition, there were no cases of induced menopause (e.g. no cases of bilateral oophorectomy). Although HAD anxiety scores were low in this healthy population, analysis of variance revealed significant differences in HAD anxiety (F(2,117) = 4.09, p = 0.02), with significantly higher scores in premenopausal female participants than in male participants (p = 0.005). In addition, there was a significant difference in the frequency of more than minimal PHQ-15 scores, with higher frequency in female participants than in male participants (χ2 = 6.7, p = 0.03) (Table 2).
Sex and menopausal status effects on brainstem connectivity
Significant differences in brainstem-network connectivity according to sex and menopausal status are summarized in Fig. 1.
Sex and menopausal status effects on brainstem connectivity. a Brainstem-network connectivity showing a significant overall sex difference (all females vs males). Left LCC connectivity with the ECN was higher in female participants than in male participants (pfwe = 0.02); however, a subgroup analysis indicated that this difference was mainly driven by premenopausal female participants (pfwe = 0.008). b Brainstem-network connectivity showing a significant difference between premenopausal and postmenopausal female participants and/or a sex difference dependent on menopausal status. PAG connectivity with the SMN was higher in premenopausal female participants than in male (pfwe = 0.03) and postmenopausal female participants (pfwe = 0.007). In contrast, PAG connectivity with the DMN was significantly higher in postmenopausal female participants than in male participants (pfwe = 0.03). This simplified figure does not show the distributed and somewhat overlapping nature of the regions in the networks (which are listed in Table 1). DMN default mode network, ECN executive control network, pfwe familywise error-corrected p-value, LCC locus coeruleus complex, PAG periaqueductal gray, SMN sensorimotor network
NPC analysis revealed a significant overall sex difference in left LCC-ECN connectivity (pfwe = 0.02), with higher connectivity in all female participants than in male participants. However, subgroup analysis according to menopausal status revealed that the overall sex difference was mainly due to significantly higher connectivity in premenopausal female participants than in male participants (pfwe = 0.008). No significant difference was observed between postmenopausal female and male participants (pfwe = 0.69). There was a trend toward higher connectivity in premenopausal than in postmenopausal female participants (pfwe = 0.051).
NPC analysis revealed significant sex differences in PAG connectivity dependent on menopausal status. Specifically, PAG-DMN connectivity was significantly higher in postmenopausal female participants compared to male participants, a difference not seen between premenopausal female and male participants (pfwe = 0.03; pfwe = 0.35, respectively). Additionally, PAG-SMN connectivity was significantly higher in premenopausal female participants than in male (pfwe = 0.03) and postmenopausal female participants (pfwe = 0.007).
Relationships between estrogens and brainstem connectivity
As indicated above, left LCC-ECN, PAG-DMN, and PAG-SMN connectivity showed significant group differences and were, thus, submitted to PLSC analysis in a subset of participants with estrogen data to evaluate relationships between connectivity and estrogen levels within each group. Specifically, estrogen and connectivity data for the three connections were simultaneously submitted to PLSC analysis in each participant group (i.e. males, premenopausal females, postmenopausal females). A summary of estrogen levels is provided in Additional file 1; as expected, levels were lower, with a more restricted range, in male and postmenopausal female participants compared to that in premenopausal female participants.
The first latent variable in the PLSC analysis of connectivity-estrogen relationships in male participants accounted for 47.3% of the cross-block variance (p = 0.05) and revealed that higher levels of free 2MeOE2 in plasma, and free and total 3MeOE1 in stool, were reliably associated with higher left LCC-ECN and PAG-DMN connectivity on bootstrap testing (Fig. 2). No other latent variables had statistically reliable saliences on bootstrap testing.
Relationships between estrogen levels and connectivity in males. In male participants, the first latent variable reflected estrogen levels reliably associated with left LCC-ECN and PAG-DMN connectivity; no other associations were statistically reliable. Heatmaps indicate the BSR of each evaluated estrogen and estrogen metabolite for the connectivity pattern shown in the corresponding bar graph, with statistically reliable positive ratios (at p < 0.01, i.e. BSR > 2.58) highlighted in red. Additional estrogens and estrogen metabolites with a BSR meeting a lower threshold of p < 0.05 are shown in light blue (BSR < − 1.96) or pink (BSR > 1.96). The plasma total levels of 16epiE3 and 17epiE3 were removed from analysis due to false positives (black bars). The bars in the bar graphs indicate the correlation between each connection and the estrogen pattern shown in the corresponding heatmap and the whiskers indicate the 95% confidence interval (confidence intervals entirely above or below zero were considered to indicate significance). This simplified figure does not show the distributed and somewhat overlapping nature of the regions in the networks (which are listed in Table 1). BSR bootstrap ratio, DMN default mode network, ECN executive control network, LCC locus coeruleus complex, PAG periaqueductal gray, SMN sensorimotor network, E1 estrone, E2 estradiol, E3 estriol: 2OHE1 2-hydroxyestrone, 2MeOE1 2-methoxyestrone, 2OHE2 2-hydroxyestradiol, 2MeOE2 2-methoxyestradiol, 3MeOE1 2-hydroxyestrone-3-methyl ether, 4OHE1 4-hydroxyestrone, 4MeOE1 4-methoxyestrone, 4MeOE2 4-methoxyestradiol, 16aOHE1 16α-hydroxyestrone, 17epiE3 17-epiestriol, 16ketoE2 16-ketoestradiol, 16epiE3 16-epiestriol, Total, summation of all free or total (free + conjugated) estrogens and estrogen metabolites, dependent on category
The first latent variable in the PLSC analysis of connectivity-estrogen relationships in premenopausal female participants accounted for 43.9% of the cross-block variance (p = 0.49) and revealed that higher levels of free E1 and 2OHE1 in plasma and free 2OHE1 and 4OHE1 in stool, were reliably associated with higher PAG-DMN and PAG-SMN connectivity on bootstrap testing. The second latent variable accounted for 36.8% of the cross-block variance, (p = 0.47) and revealed that lower levels of total 2OHE2 and 4MeOE1 in plasma, and lower levels of free 4MeOE2 in stool, were reliably associated with higher left LCC-ECN connectivity on bootstrap testing (Fig. 3).
Relationships between estrogen levels and connectivity in premenopausal females. In premenopausal female participants, the first latent variable reflected estrogen levels reliably associated with PAG-SMN and PAG-DMN connectivity, while the second variable reflected estrogen levels reliably associated with left LCC-ECN connectivity. Heatmaps indicate the BSR of each evaluated estrogen and estrogen metabolite for the connectivity pattern shown in the corresponding bar graph, with statistically reliable ratios (at p < 0.01) highlighted in dark blue (BSR < − 2.58) or red (BSR > 2.58). Additional estrogens and estrogen metabolites with a BSR meeting a lower threshold of p < 0.05 are shown in light blue (BSR < − 1.96) or pink (BSR > 1.96). The plasma total levels of 16epiE3 and 17epiE3 were removed from analysis due to false positives (black bars). The bars in the bar graphs indicate the correlation between each connection and the estrogen pattern shown in the corresponding heatmap and the whiskers indicate the 95% confidence interval (confidence intervals entirely above or below zero were considered to indicate significance). Note: This simplified figure does not show the distributed and somewhat overlapping nature of the regions in the networks (which are listed in Table 1). BSR bootstrap ratio, DMN default mode network, ECN executive control network, LCC locus coeruleus complex, PAG periaqueductal gray, SMN sensorimotor network, E1 estrone, E2 estradiol, E3 estriol: 2OHE1 2-hydroxyestrone, 2MeOE1 2-methoxyestrone, 2OHE2 2-hydroxyestradiol, 2MeOE2 2-methoxyestradiol, 3MeOE1 2-hydroxyestrone-3-methyl ether, 4OHE1 4-hydroxyestrone, 4MeOE1 4-methoxyestrone, 4MeOE2 4-methoxyestradiol, 16aOHE1 16α-hydroxyestrone, 17epiE3 17-epiestriol, 16ketoE2 16-ketoestradiol, 16epiE3 16-epiestriol; Total, summation of all free or total (free + conjugated) estrogens and estrogen metabolites, dependent on category
The first latent variable in the PLSC analysis of connectivity-estrogen relationships in postmenopausal female participants accounted for 43.8% of the cross-block variance (p = 0.05) and revealed that higher levels of total E2 and lower levels of free 16aOHE1 in plasma, and higher levels of free 16epiE3, total 3MeOE1, and total-free and total-total estrogens (summation all free estrogens and estrogen metabolites and summation of all free + conjugated estrogens and estrogen metabolites, respectively) in stool, were reliably associated with higher left LCC-ECN connectivity on bootstrap testing. The second latent variable accounted for 30.2% of the cross-block variance (p = 0.21) and revealed that higher levels of free E2, total 16aOHE1, and total-free and total-total estrogens in plasma were reliably associated with higher PAG-SMN connectivity on bootstrap testing (Fig. 4).
Relationships between estrogen levels and connectivity in postmenopausal females. In postmenopausal female participants, the first latent variable reflected estrogen levels reliably associated with left LCC-ECN connectivity, while the second variable reflected estrogen levels reliably associated with PAG-SMN connectivity. Heatmaps indicate the BSR of each evaluated estrogen and estrogen metabolite for the connectivity pattern shown in the corresponding bar graph, with statistically reliable ratios (at p < 0.01) highlighted in dark blue (BSR < − 2.58) or red (BSR > 2.58). Additional estrogens and estrogen metabolites with a BSR meeting a lower threshold of p < 0.05 are shown in light blue (BSR < − 1.96) or pink (BSR > 1.96). The plasma total levels of 16epiE3 and 17epiE3 were removed from analysis due to false positives (black bars). The bars in the bar graphs indicate the correlation between each connection and the estrogen pattern shown in the corresponding heatmap and the whiskers indicate the 95% confidence interval (confidence intervals entirely above or below zero were considered to indicate significance). Note: This simplified figure does not show the distributed and somewhat overlapping nature of the regions in the networks (which are listed in Table 1). BSR bootstrap ratio, DMN default mode network, ECN executive control network, LCC locus coeruleus complex, PAG periaqueductal gray, SMN sensorimotor network, E1 estrone, E2 estradiol, E3 estriol: 2OHE1 2-hydroxyestrone, 2MeOE1 2-methoxyestrone, 2OHE2 2-hydroxyestradiol, 2MeOE2 2-methoxyestradiol, 3MeOE1 2-hydroxyestrone-3-methyl ether, 4OHE1 4-hydroxyestrone, 4MeOE1 4-methoxyestrone, 4MeOE2 4-methoxyestradiol, 16aOHE1 16α-hydroxyestrone, 17epiE3 17-epiestriol, 16ketoE2 16-ketoestradiol, 16epiE3 16-epiestriol; Total, summation of all free or total (free + conjugated) estrogens and estrogen metabolites, dependent on category
Relationships between anxiety and somatic symptoms and brainstem connectivity
Relationships between left LCC-ECN, PAG-DMN, and PAG-SMN connectivity and anxiety and somatic symptoms were similarly examined in the total sample and in each participant group (i.e. males, premenopausal females, postmenopausal females).
The first latent variable in the PLSC analysis of relationships between brainstem connectivity and anxiety/somatic symptoms in the total sample accounted for 44.4% of the cross-block variance (p = 0.10) and revealed that higher HAD anxiety (BSR = 1.97, p = 0.049) was reliably associated higher left LCC-ECN connectivity on bootstrap testing (r = 0.10, 95% confidence interval [95% CI] 0.02–0.26). No other latent variables had statistically reliable saliences on bootstrap testing.
The first latent variable in the PLSC analysis of relationships between brainstem connectivity and anxiety and somatic symptoms in male participants accounted for 65.8%of the cross-block variance (p = 0.002) and revealed that higher perceived stress (BSR = 2.78, p = 0.005) and trait anxiety (BSR = 4.05, p < 0.001) were reliably associated with higher PAG-DMN connectivity on bootstrap testing (r = 0.40, 95% CI 0.21–0.59). No other latent variables had statistically reliable saliences on bootstrap testing.
The first latent variable in the PLSC analysis of relationships between brainstem connectivity and anxiety and somatic symptoms in premenopausal female participants accounted for 54.5% of the cross-block variance (p = 0.06) and revealed that higher PHQ-15 scores (BSR = 2.29, p = 0.02) were reliably associated with higher left LCC-ECN connectivity on bootstrap testing (r = 0.21, 95% CI 0.08–0.44). In addition, the second latent variable accounted for 36.1% of the cross-block variance (p = 0.36) and revealed that higher PILL scores (BSR = 2.40, p = 0.02) were reliably associated with higher PAG-SMN connectivity on bootstrap testing (r = 0.28, 95% CI 0.02–0.46).
PLSC analysis of relationships between brainstem connectivity and anxiety and somatic symptoms in postmenopausal female participants did not reveal any statistically reliable associations.
Discussion
The present study evaluated differences in resting-state connectivity between specific brainstem nuclei (DRN, left and right LCC, and PAG) and key brain networks (CAN, DMN, EAN, ECN, SAL, and SMN) implicated in stress responsiveness and pain modulation according to sex and menopausal status in healthy individuals. We found a significant overall sex difference in left LCC-ECN connectivity, with generally higher connectivity in female participants than in male participants, mainly driven by higher connectivity in premenopausal female participants. There was also significantly higher PAG-SMN connectivity in premenopausal female participants than in male participants and in premenopausal female participants than in postmenopausal female participants. In contrast, significantly higher PAG-DMN connectivity was present in postmenopausal female participants than in male participants. Further, relationships between left LCC-ECN, PAG-SMN, and PAG-DMN and estrogen levels in plasma and stool, as well as anxiety and somatic symptoms, were observed.
Sex differences in LCC resting-state connectivity may relate to known differences in information processing priorities
The left LCC showed higher connectivity with the ECN in female participants, especially premenopausal females, than in male participants. The ECN comprises lateral prefrontal and parietal regions and supports executive functions such as working memory, selective attention, and cognitive control [46]. Animal and human studies suggest that enhanced functional coupling of the LCC with the ECN is associated with increased goal-directed attention and decreased impulsivity [47, 48].
LCC noradrenergic signaling biases perception, attention, and memory toward more salient stimuli by selectively amplifying the activity of priority mechanisms operating at the moment [49]. Substantial evidence suggests that while males show a greater preference for allocentric knowledge (e.g., describing objects/others independent of one’s own perspective), with a greater reliance on hippocampal-based strategies [8], females show a greater preference for egocentric knowledge (e.g., describing objects in terms of one’s own spatial perspective and using more privileged information in inference), with a greater reliance on working memory processes mediated by frontal regions, such as those in the ECN [8, 50]. We speculate that these previously reported sex-related biases in information processing priorities are operating under resting-state condition and underlie the observed difference in left LCC connectivity between males and females. Further, anxiety and egocentricism are related, with greater reliance on egocentric perspective-taking/mentalizing in those experiencing anxiety [51]. Thus, under our speculation, the observed greater LCC-ECN connectivity in premenopausal females than in males may be related to the higher current anxiety symptom scores in premenopausal female participants than in male participants. Consistent with this, higher HAD anxiety scores were associated with left LCC-ECN connectivity in the total sample, with a trend (p = 0.08) towards association among premenopausal female participants. In addition, left LCC-ECN connectivity was associated with higher PHQ-15 in premenopausal female participants. Although the PHQ-15 is not specifically focused on anxiety, it includes some items on somatic symptoms of anxiety (e.g. heart racing, shortness of breath, nausea, dizziness, stomachaches) and had shared variance with HAD anxiety in the present sample of premenopausal females (r = 0.48).
Interestingly, unlike the left LCC, the right LCC did not show any sex differences in connectivity. A recent mixed-sex study of LCC connectivity gradients found greater relationships between age, anxiety/depression symptoms, and cognitive performance for the left LCC than for the right LCC [52]. Additional research suggests that neurodegenerative disorders affect the left LCC more than the right LCC [53, 54]. Thus, the left LCC may be more pliable or sensitive than the right LCC. However, further research is needed.
In premenopausal female participants, decreased levels of mainly metabolites in the 2-hydroxylation pathway of estrogen metabolites, including plasma 2OHE2, were associated with higher left LCC-ECN connectivity. Metabolites of the 2-hydroxylation pathway are generally considered to have weak estrogenic activity; however, 2OHE2 has structural similarities to catecholamines and can compete with noradrenaline in the brain [55]. Thus, in premenopausal female participants, circulating 2OHE2 may interact with LC noradrenergic output, modulating left LCC-ECN connectivity. In contrast, in male participants, higher levels of methylated estrogen metabolites, including 2MeOE2 in plasma, and in postmenopausal female participants, higher levels of total E2 in plasma, were associated with higher left LCC-ECN connectivity (i.e., connectivity more similar to that of premenopausal female participants). 2MeOE2 was one of the more abundant estrogen metabolites in plasma in the male participants in the present study and E2 is a major endogenous estrogen. These results suggest that estrogens may affect LCC-ECN connectivity in males and postmenopausal females, but the effects may be more apparent for metabolites with relatively higher levels, as levels were generally lower in these individuals.
Sex differences in PAG-SMN resting-state connectivity may relate to known sex differences in pain processing
The PAG showed higher connectivity with the SMN in premenopausal female participants than in male and postmenopausal female participants. The present results are consistent with a previous neuroimaging study that reported greater PAG connectivity with sensorimotor-related brain regions in healthy females than in healthy males [56], and add to the literature by showing a dependence on menopausal status.
The SMN comprises sensorimotor, mid-cingulate and superior frontal cortices, as well as the posterior insula, thalamus, and basal ganglia [57, 58]. The SMN is involved in central processing and modulation of visceral and somatic sensory information and both the PAG and SMN are involved in pain processing. Previous studies in patient populations indicate that increased connectivity of the PAG with the SMN, or specific regions within the SMN, may be associated with an increased risk of the development of chronic pain following mild traumatic brain injury [59] and increased central sensitization symptoms in patients with fibromyalgia [60]. In the present study, we found that higher PAG-SMN connectivity was associated with higher PILL scores, suggesting a relationship between PAG-SMN connectivity and increased somatic symptom burden, in premenopausal female participants. Numerous studies suggest that somatic symptom burden contributes to an increased risk for the development of chronic pain conditions [61, 62]. Thus, our findings align with an increased vulnerability to chronic pain conditions in premenopausal females than in males [23, 63]. However, a mixed-sex study in healthy individuals reported that resting-state connectivity between the PAG and SMN positively correlated with conditioned pain modulation, suggesting more efficient endogenous pain modulation with increased connectivity [64]. Primary somatosensory cortex modulates sensory gain and nociception, with outputs originating from layer 5 of the cortex and projecting to subcortical targets, including the PAG, comprising an anti-nociceptive pathway, and outputs originating from layer 6 of the cortex and projecting to the thalamus, which is also a component of the SMN and interacts with the PAG, comprising a pro-nociceptive pathway [65]. Thus, the interpretation of increased PAG-SMN connectivity is complicated and may require a finer-grained analysis.
However, one notable finding in the analysis of connectivity-estrogen relationships was that, in premenopausal female participants, increased plasma and stool free 2OHE1 was associated with increased PAG-SMN connectivity. 2OHE1 has been shown to increase nociceptor activation via transient receptor potential ankyrin 1Â (TRPA1) and transient receptor potential vanilloid type-1 (TRPV1) channels in a mouse model of uterine pain [66]. Peripheral TRPA1 and TRPV1 sensitization have been implicated in visceral pain disorders, including DGBI [67,68,69]. However, central factors may be necessary for the persistence of visceral hypersensitivity [70]. TRPV1 receptors are expressed in the PAG and antagonists applied to the dorsolateral PAG can reduce anxiety-like and nociceptive behavior in animal models [71, 72]. Thus, estrogen metabolites may modulate the increased PAG-SMN connectivity in premenopausal females, potentially heightening somatic symptom burden and risk for chronic pain disorders such as DGBI.
Sex differences in PAG-DMN resting-state connectivity may relate to known sex differences in risk for posttraumatic stress disorder and stress-induced analgesia
In the present study, PAG-DMN connectivity was significantly lower in male participants than in postmenopausal female participants; however, there was a positive association between PAG-DMN connectivity and trait anxiety (assessed by the STAI) and perceived stress (assessed by the PSS) in male participants, suggesting that, in males, acute and chronic feelings of worry and distress is associated with connectivity similar to that of premenopausal female participants. The DMN comprises the medial prefrontal cortex, posterior cingulate cortex, precuneus, inferior parietal cortices, and lateral temporal cortices, and is involved in self-referential processes [73, 74]. Additionally, building upon animal research, human brain imaging studies implicate DMN regions in threat processing, with the posterior cingulate cortex and medial prefrontal cortex involved in evaluating threat cues and modulating responses to threat, respectively [75]. A mixed-sex neuroimaging study reported increased PAG-DMN connectivity two weeks after a car accident as predictive of the development of posttraumatic stress disorder within 6Â months [76]. In female patients with temporomandibular disorder, higher connectivity between PAG and posterior cingulate and medial prefrontal cortices is associated with pain rumination [77]. Thus, we speculate that increased PAG-DMN connectivity in postmenopausal females, and males with higher trait anxiety, may be associated with maladaptive coping strategies that contribute to risk for posttraumatic stress disorder [78]. Generally, females have higher risk for posttraumatic stress disorder than males [79]. Some studies suggest that risk for posttraumatic stress disorder in females peaks around menopause [80], which may explain why the sex difference in PAG-DMN connectivity was more robust for postmenopausal female participants than for premenopausal female participants.
Limitations
The present study has several limitations. The number of postmenopausal female participants was relatively small, limiting the power to detect differences between postmenopausal and premenopausal females and sex differences that emerge or reverse after menopause. Additionally, although we controlled for age, bias may exist in comparisons between males and postmenopausal females because of limited overlap in age; thus, future studies with a greater representation of older males is needed to confirm the robustness of our results. Further, in the PLSC analyses, follicular and luteal females were combined into a single group due to the limited sample size, resulting in heterogeneity that may have contributed to non-significance on permutation testing. Therefore, we focused on salience reliability in our reporting. A larger sample, allowing follicular and luteal females to represented in separate groups, may result in improved power to detect significant latent variables that more precisely explain the variance. Additionally, although sex was treated as binary in the present analysis, sex is not strictly binary, as variation exists (e.g., in chromosomes, reproductive organs, hormones) that defies binary classification. Other limitations arise from the small size of brainstem nuclei, which can be difficult to delineate. We used a brainstem atlas based mostly on postmortem data to delineate the LCC, DRN, and PAG; however, imprecision in nuclei boundaries may have affected the connectivity estimates. Additionally, we investigated the overall connectivity of the LCC, DRN, and PAG, without consideration of differential connectivity within each of these brainstem regions. Although small, these regions show variations in connectivity, with a rostral-caudal connectivity gradient in the LC and subregions in the DRN and PAG with differential connectivity supporting various functions [81]. This may have contributed to the lack of significant sex differences in DRN connectivity. However, a finer-grained analysis is beyond the scope of the present study. Additionally, as a major limitation of correlational studies, the causality or directionality of interactions could not be addressed. Finally, the present study was exploratory in nature and estrogen data were available in a limited subset of participants; thus, the results should be interpreted with caution and further research is required to rigorously confirm the present findings.
Conclusions
The present study expands the limited research on sex and menopausal effects on brainstem connectivity, and their relationships with various estrogens, in humans. We found that healthy females show higher left LCC and PAG connectivity with networks involved in cognitive control, and sensorimotor function and self-relevant processes, respectively, than males, dependent on their menopausal status. Although such differences may show benefits under optimal conditions, they may also relate to differential vulnerabilities to chronic pain and stress-sensitive disorders at different life stages. In particular, PAG connectivity with the SMN may be modulated by circulating 2OHE1 and associated with somatic symptoms in premenopausal females. Given the known role of 2OHE1 in peripheral and central sensitization processes, we speculate that this contributes to an increased risk for chronic pain disorders such as DGBI in premenopausal females. However, future studies are needed in patients with chronic pain.
Availability of data and materials
Data are available from the corresponding authors on reasonable request.
References
Tache Y, Martinez V, Million M, Rivier J. Corticotropin-releasing factor and the brain-gut motor response to stress. Can J Gastroenterol. 1999;13:18A-25A. https://doiorg.publicaciones.saludcastillayleon.es/10.1155/1999/375916.
Hale MW, Shekhar A, Lowry CA. Stress-related serotonergic systems: implications for symptomatology of anxiety and affective disorders. Cell Mol Neurobiol. 2012;32(5):695–708. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s10571-012-9827-1.
Mayer EA, Tillisch K. The brain-gut axis in abdominal pain syndromes. Annu Rev Med. 2011;62:381–96. https://doiorg.publicaciones.saludcastillayleon.es/10.1146/annurev-med-012309-103958.
George DT, Ameli R, Koob GF. Periaqueductal gray sheds light on dark areas of psychopathology. Trends Neurosci. 2019;42(5):349–60. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.tins.2019.03.004.
Peng WH, Kan HW, Ho YC. Periaqueductal gray is required for controlling chronic stress-induced depression-like behavior. Biochem Biophys Res Commun. 2022;593:28–34. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.bbrc.2022.01.025.
Mouton PR, Pakkenberg B, Gundersen HJ, Price DL. Absolute number and size of pigmented locus coeruleus neurons in young and aged individuals. J Chem Neuroanat. 1994;7(3):185–90. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/0891-0618(94)90028-0.
Michelsen KA, Schmitz C, Steinbusch HW. The dorsal raphe nucleus–from silver stainings to a role in depression. Brain Res Rev. 2007;55(2):329–42. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.brainresrev.2007.01.002.
Ycaza Herrera A, Wang J, Mather M. The gist and details of sex differences in cognition and the brain: how parallels in sex differences across domains are shaped by the locus coeruleus and catecholamine systems. Prog Neurobiol. 2019;176:120–33. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.pneurobio.2018.05.005.
Yu W, Pati D, Pina MM, Schmidt KT, Boyt KM, Hunker AC, et al. Periaqueductal gray/dorsal raphe dopamine neurons contribute to sex differences in pain-related behaviors. Neuron. 2021;109(8):1365–80. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.neuron.2021.03.001.
Bangasser DA, Curtis A, Reyes BA, Bethea TT, Parastatidis I, Ischiropoulos H, et al. Sex differences in corticotropin-releasing factor receptor signaling and trafficking: potential role in female vulnerability to stress-related psychopathology. Mol Psychiatry. 2010;15(9):96–904. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/mp.2010.66.
Bangasser DA, Eck SR, Ordones SE. Sex differences in stress reactivity in arousal and attention systems. Neuropsychopharmacology. 2019;44(1):129–39. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41386-018-0137-2.
Zhang S, Hu S, Chao HH, Li CS. Resting-state functional connectivity of the locus coeruleus in humans: in comparison with the ventral tegmental area/substantia Nigra pars compacta and the effects of age. Cereb Cortex. 2016;26(8):3413–27. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/cercor/bhv172.
Sun P, Wang J, Zhang M, Duan X, Wei Y, Xu F, et al. Sex-related differential whole-brain input atlas of locus coeruleus noradrenaline neurons. Front Neural Circuits. 2020;14:53. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fncir.2020.00053.
Osborne NR, Davis KD. Sex and gender differences in pain. Int Rev Neurobiol. 2022;164:277–307. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/bs.irn.2022.06.013.
Bangasser DA, Wiersielis KR, Khantsis S. Sex differences in the locus coeruleus-norepinephrine system and its regulation by stress. Brain Res. 2016;1641:177–88. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.brainres.2015.11.021.
McEwen B. Estrogen actions throughout the brain. Recent Prog Horm Res. 2002;57:357–84. https://doiorg.publicaciones.saludcastillayleon.es/10.1210/rp.57.1.357.
Donner N, Handa RJ. Estrogen receptor beta regulates the expression of tryptophan-hydroxylase 2 mRNA within serotonergic neurons of the rat dorsal raphe nuclei. Neuroscience. 2009;163(2):705–18. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.neuroscience.2009.06.046.
Campos GV, de Noronha SR, de Souza AA, Lima PM, Abreu AR, Chianca-Jr D, de Menezes RC. Estrogen receptor beta activation within dorsal raphe nucleus reverses anxiety-like behavior induced by food restriction in female rats. Behav Brain Res. 2019;357–358:57–64. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.bbr.2018.03.029.
Averitt DL, Eidson LN, Doyle HH, Murphy AZ. Neuronal and glial factors contributing to sex differences in opioid modulation of pain. Neuropsychopharmacology. 2019;44(1):155–65. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41386-018-0127-4.
Alblooshi S, Taylor M, Gill N. Does menopause elevate the risk for developing depression and anxiety? Results from a systematic review. Australas Psychiatry. 2023;31(2):165–73. https://doiorg.publicaciones.saludcastillayleon.es/10.1177/10398562231165439.
Gibson CJ, Li Y, Bertenthal D, Huang AJ, Seal KH. Menopause symptoms and chronic pain in a national sample of midlife women veterans. Menopause. 2019;26(7):708–13. https://doiorg.publicaciones.saludcastillayleon.es/10.1097/GME.0000000000001312.
Labrenz F, Merz CJ, Icenhour A. Connecting dots in disorders of gut-brain interaction: the interplay of stress and sex hormones in shaping visceral pain. Front Psychiatry. 2023;14:1204136. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fpsyt.2023.1204136.
Sarnoff RP, Hreinsson JP, Sperber AD, Palsson OS, Bangdiwala SI, Chang L. Tu1029 effect of menopause and menses on gastrointestinal symptoms among individuals with DGBI in the global population. Gastroenterology. 2024;166(5):S1228-9. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/S0016-5085(24)03273-6.
Cain KC, Jarrett ME, Burr RL, Rosen S, Hertig VL, Heitkemper MM. Gender differences in gastrointestinal, psychological, and somatic symptoms in irritable bowel syndrome. Dig Dis Sci. 2009;54(7):1542–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s10620-008-0516-3.
Sheehan DV, Lecrubier Y, Sheehan KH, Amorim P, Janavs J, Weiller E, et al. The mini-international neuropsychiatric interview (M.I.N.I.): the development and validation of a structured diagnostic psychiatric interview for DSM-IV and ICD-10. J Clin Psychiatry. 1998;59(Suppl 20):22–33.
Harlow SD, Gass M, Hall JE, Lobo R, Maki P, Rebar RW, et al. Executive summary of the stages of reproductive aging workshop +10: addressing the unfinished agenda of staging reproductive aging. Climacteric. 2012;15(2):105–14. https://doiorg.publicaciones.saludcastillayleon.es/10.3109/13697137.2011.650656.
Zigmond AS, Snaith RP. The hospital anxiety and depression scale. Acta Psychiatr Scand. 1983;67(6):361–70. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1600-0447.1983.tb09716.x.
Spielberg CD, Gorsuch RL, Lushene R, Vagg PR, Jacobs GA. Manual for the state-trait anxiety inventory. Washington, DC: Consulting Psychologists Press; 1983.
Cohen S, Kamarck T, Mermelstein R. A global measure of perceived stress. J Health Soc Behav. 1983;24(4):385–96.
Kroenke K, Spitzer RL, Williams JB. The PHQ-15: validity of a new measure for evaluating the severity of somatic symptoms. Psychosom Med. 2002;64(2):258–66.
Pennebaker JW. The psychology of physical symptoms. NewYork: Springer; 1982.
Barch DM, Albaugh MD, Baskin-Sommers A, Bryant BE, Clark DB, Dick AS, et al. Demographic and mental health assessments in the adolescent brain and cognitive development study: updates and age-related trajectories. Dev Cogn Neurosci. 2021;52:101031. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.dcn.2021.101031.
Glasser MF, Smith SM, Marcus DS, Andersson JL, Auerbach EJ, Behrens TE, et al. The human connectome project’s neuroimaging approach. Nat Neurosci. 2016;19(9):1175–87. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/nn.4361.
Dale AM, Fischl B, Sereno MI. Cortical surface-based analysis. I. Segmentation and surface reconstruction. Neuroimage. 1999;9(2):179–94. https://doiorg.publicaciones.saludcastillayleon.es/10.1006/nimg.1998.0395.
Fair DA, Miranda-Dominguez O, Snyder AZ, Perrone A, Earl EA, Van AN, et al. Correction of respiratory artifacts in MRI head motion estimates. Neuroimage. 2020;208:116400. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.neuroimage.2019.116400.
Destrieux C, Fischl B, Dale A, Halgren E. Automatic parcellation of human cortical gyri and sulci using standard anatomical nomenclature. Neuroimage. 2010;53(1):1–15. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.neuroimage.2010.06.010.
Edlow BL, Takahashi E, Wu O, Benner T, Dai G, Bu L, et al. Neuroanatomic connectivity of the human ascending arousal system critical to consciousness and its disorders. J Neuropathol Exp Neurol. 2012;71(6):531–46. https://doiorg.publicaciones.saludcastillayleon.es/10.1097/NEN.0b013e3182588293.
Xu X, Roman JM, Issaq HJ, Keefer LK, Veenstra TD, Ziegler RG. Quantitative measurement of endogenous estrogens and estrogen metabolites in human serum by liquid chromatography-tandem mass spectrometry. Anal Chem. 2007;79(20):7813–21. https://doiorg.publicaciones.saludcastillayleon.es/10.1021/ac070494j.
Xu X, Veenstra TD, Fox SD, Roman JM, Issaq HJ, Falk R, et al. Measuring fifteen endogenous estrogens simultaneously in human urine by high-performance liquid chromatography-mass spectrometry. Anal Chem. 2005;77(20):6646–54. https://doiorg.publicaciones.saludcastillayleon.es/10.1021/ac050697c.
Flores R, Shi J, Fuhrman B, Xu X, Veenstra TD, Gail MH, et al. Fecal microbial determinants of fecal and systemic estrogens and estrogen metabolites: a cross-sectional study. J Transl Med. 2012;10:253. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/1479-5876-10-253.
Meriwether D, Sulaiman D, Volpe C, Dorfman A, Grijalva V, Dorreh N, et al. Apolipoprotein A-I mimetics mitigate intestinal inflammation in COX2-dependent inflammatory bowel disease model. J Clin Invest. 2019;129(9):3670–85. https://doiorg.publicaciones.saludcastillayleon.es/10.1172/JCI123700.
Smith SM, Nichols TE. Threshold-free cluster enhancement: addressing problems of smoothing, threshold dependence and localisation in cluster inference. Neuroimage. 2009;44(1):83–98. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.neuroimage.2008.03.061.
Winkler AM, Ridgway GR, Webster MA, Smith SM, Nichols TE. Permutation inference for the general linear model. Neuroimage. 2014;92(100):381–97. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.neuroimage.2014.01.060.
Winkler AM, Webster MA, Brooks JC, Tracey I, Smith SM, Nichols TE. Non-parametric combination and related permutation tests for neuroimaging. Hum Brain Mapp. 2016;37(4):1486–511. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/hbm.23115.
McIntosh AR, Lobaugh NJ. Partial least squares analysis of neuroimaging data: applications and advances. Neuroimage. 2004;23(Suppl 1):S250–63. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.neuroimage.2004.07.020.
Menon V. Large-scale brain networks and psychopathology: a unifying triple network model. Trends Cogn Sci. 2011;15(10):483–506. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.tics.2011.08.003.
Bari A, Xu S, Pignatelli M, Takeuchi D, Feng J, Li Y, Tonegawa S. Differential attentional control mechanisms by two distinct noradrenergic coeruleo-frontal cortical pathways. Proc Natl Acad Sci USA. 2020;117(46):29080–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1073/pnas.2015635117.
Grueschow M, Kleim B, Ruff CC. Functional coupling of the locus coeruleus is linked to successful cognitive control. Brain Sci. 2022;12:3. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/brainsci12030305.
Mather M, Clewett D, Sakaki M, Harley CW. Norepinephrine ignites local hotspots of neuronal excitation: how arousal amplifies selectivity in perception and memory. Behav Brain Sci. 2016;39:e200. https://doiorg.publicaciones.saludcastillayleon.es/10.1017/S0140525X15000667.
Torromino G, Maggi A, De Leonibus E. Estrogen-dependent hippocampal wiring as a risk factor for age-related dementia in women. Prog Neurobiol. 2021;197:101895. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.pneurobio.2020.101895.
Todd AR, Forstmann M, Burgmer P, Brooks AW, Galinsky AD. Anxious and egocentric: how specific emotions influence perspective taking. J Exp Psychol Gen. 2015;144(2):374–91. https://doiorg.publicaciones.saludcastillayleon.es/10.1037/xge0000048.
Vereb D, Mijalkov M, Canal-Garcia A, Chang YW, Gomez-Ruis E, Gerboles BZ, et al. Age-related differences in the functional topography of the locus coeruleus: implications for cognitive and affective functions. medRxiv. 2023. https://doiorg.publicaciones.saludcastillayleon.es/10.1101/2023.02.25.23286442.
Sun J, Ma J, Gao L, Wang J, Zhang D, Chen L, et al. Disruption of locus coeruleus-related functional networks in Parkinson’s disease. NPJ Parkinsons Dis. 2023;9(1):81. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41531-023-00532-x.
Jacobs HI, Wiese S, van de Ven V, Gronenschild EH, Verhey FR, Matthews PM. Relevance of parahippocampal-locus coeruleus connectivity to memory in early dementia. Neurobiol Aging. 2015;36(2):618–26. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.neurobiolaging.2014.10.041.
Paden CM, McEwen BS, Fishman J, Snyder L, DeGroff V. Competition by estrogens for catecholamine receptor binding in vitro. J Neurochem. 1982;39(2):512–20. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1471-4159.1982.tb03974.x.
Coulombe MA, Erpelding N, Kucyi A, Davis KD. Intrinsic functional connectivity of periaqueductal gray subregions in humans. Hum Brain Mapp. 2016;37(4):1514–30. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/hbm.23117.
Biswal B, Yetkin FZ, Haughton VM, Hyde JS. Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magn Reson Med. 1995;34(4):537–41. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/mrm.1910340409.
Jones LM, Fontanini A, Sadacca BF, Miller P, Katz DB. Natural stimuli evoke dynamic sequences of states in sensory cortical ensembles. Proc Natl Acad Sci U S A. 2007;104(47):18772–7. https://doiorg.publicaciones.saludcastillayleon.es/10.1073/pnas.0705546104.
Bosak N, Branco P, Kuperman P, Buxbaum C, Cohen RM, Fadel S, et al. Brain connectivity predicts chronic pain in acute mild traumatic brain injury. Ann Neurol. 2022;92(5):819–33. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/ana.26463.
Soldatelli M, Franco AO, Picon F, Duarte JA, Scherer R, Bandeira J, et al. Primary somatosensory cortex and periaqueductal gray functional connectivity as a marker of the dysfunction of the descending pain modulatory system in fibromyalgia. Korean J Pain. 2023;36(1):113–27. https://doiorg.publicaciones.saludcastillayleon.es/10.3344/kjp.22225.
Fillingim RB, Ohrbach R, Greenspan JD, Sanders AE, Rathnayaka N, Maixner W, Slade GD. Associations of psychologic factors with multiple chronic overlapping pain conditions. J Oral Facial Pain Headache. 2020;34(Suppl):s85–100. https://doiorg.publicaciones.saludcastillayleon.es/10.11607/ofph.2584.
Fillingim RB, Ohrbach R, Greenspan JD, Knott C, Diatchenko L, Dubner R, et al. Psychological factors associated with development of TMD: the OPPERA prospective cohort study. J Pain. 2013;14(12 Suppl):T75-90. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.jpain.2013.06.009.
Sperber AD, Bangdiwala SI, Drossman DA, Ghoshal UC, Simren M, Tack J, et al. Worldwide prevalence and burden of functional gastrointestinal disorders, results of rome foundation global study. Gastroenterology. 2021;160(1):99–114. https://doiorg.publicaciones.saludcastillayleon.es/10.1053/j.gastro.2020.04.014.
Harrison R, Gandhi W, van Reekum CM, Salomons TV. Conditioned pain modulation is associated with heightened connectivity between the periaqueductal grey and cortical regions. Pain Rep. 2022;7(3):e999. https://doiorg.publicaciones.saludcastillayleon.es/10.1097/PR9.0000000000000999.
Ziegler K, Folkard R, Gonzalez AJ, Burghardt J, Antharvedi-Goda S, Martin-Cortecero J, et al. Primary somatosensory cortex bidirectionally modulates sensory gain and nociceptive behavior in a layer-specific manner. Nat Commun. 2023;14(1):2999. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41467-023-38798-7.
Xie Z, Feng J, Cai T, McCarthy R, Eschbach MD 2nd, Wang Y, et al. Estrogen metabolites increase nociceptor hyperactivity in a mouse model of uterine pain. JCI Insight. 2022;7:10. https://doiorg.publicaciones.saludcastillayleon.es/10.1172/jci.insight.149107.
Peng HY, Chang HM, Lee SD, Huang PC, Chen GD, Lai CH, et al. TRPV1 mediates the uterine capsaicin-induced NMDA NR2B-dependent cross-organ reflex sensitization in anesthetized rats. Am J Physiol Renal Physiol. 2008;295(5):F1324–35. https://doiorg.publicaciones.saludcastillayleon.es/10.1152/ajprenal.00126.2008.
Ford AC, Vanner S, Kashyap PC, Nasser Y. Chronic visceral pain: new peripheral mechanistic insights and resulting treatments. Gastroenterology. 2024;166(6):976–94. https://doiorg.publicaciones.saludcastillayleon.es/10.1053/j.gastro.2024.01.045.
Balemans D, Boeckxstaens GE, Talavera K, Wouters MM. Transient receptor potential ion channel function in sensory transduction and cellular signaling cascades underlying visceral hypersensitivity. Am J Physiol Gastrointest Liver Physiol. 2017;312(6):G635–48. https://doiorg.publicaciones.saludcastillayleon.es/10.1152/ajpgi.00401.2016.
van Wanrooij SJ, Wouters MM, Van Oudenhove L, Vanbrabant W, Mondelaers S, Kollmann P, et al. Sensitivity testing in irritable bowel syndrome with rectal capsaicin stimulations: role of TRPV1 upregulation and sensitization in visceral hypersensitivity? Am J Gastroenterol. 2014;109(1):99–109. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/ajg.2013.371.
Terzian AL, Aguiar DC, Guimaraes FS, Moreira FA. Modulation of anxiety-like behaviour by transient receptor potential vanilloid type 1 (TRPV1) channels located in the dorsolateral periaqueductal gray. Eur Neuropsychopharmacol. 2009;19(3):188–95. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.euroneuro.2008.11.004.
Madasu MK, Okine BN, Olango WM, Rea K, Lenihan R, Roche M, Finn DP. Genotype-dependent responsivity to inflammatory pain: a role for TRPV1 in the periaqueductal grey. Pharmacol Res. 2016;113:44–54. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.phrs.2016.08.011.
Greicius MD, Krasnow B, Reiss AL, Menon V. Functional connectivity in the resting brain: a network analysis of the default mode hypothesis. Proc Natl Acad Sci U S A. 2003;100(1):253–8. https://doiorg.publicaciones.saludcastillayleon.es/10.1073/pnas.0135058100.
Buckner RL, Andrews-Hanna JR, Schacter DL. The brain’s default network: anatomy, function, and relevance to disease. Ann N Y Acad Sci. 2008;1124:1–38. https://doiorg.publicaciones.saludcastillayleon.es/10.1196/annals.1440.011.
Wen Z, Pace-Schott EF, Lazar SW, Rosen J, Ahs F, Phelps EA, et al. Distributed neural representations of conditioned threat in the human brain. Nat Commun. 2024;15(1):2231. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41467-024-46508-0.
Webb EK, Huggins AA, Belleau EL, Taubitz LE, Hanson JL, deRoon-Cassini TA, Larson CL. Acute posttrauma resting-state functional connectivity of periaqueductal gray prospectively predicts posttraumatic stress disorder symptoms. Biol Psychiatry Cogn Neurosci Neuroimaging. 2020;5(9):891–900. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.bpsc.2020.03.004.
Kucyi A, Moayedi M, Weissman-Fogel I, Goldberg MB, Freeman BV, Tenenbaum HC, Davis KD. Enhanced medial prefrontal-default mode network functional connectivity in chronic pain and its association with pain rumination. J Neurosci. 2014;34(11):3969–75. https://doiorg.publicaciones.saludcastillayleon.es/10.1523/JNEUROSCI.5055-13.2014.
Jenness JL, Jager-Hyman S, Heleniak C, Beck AT, Sheridan MA, McLaughlin KA. Catastrophizing, rumination, and reappraisal prospectively predict adolescent PTSD symptom onset following a terrorist attack. Depress Anxiety. 2016;33(11):1039–47. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/da.22548.
Olff M. Sex and gender differences in post-traumatic stress disorder: an update. Eur J Psychotraumatol. 2017. https://doiorg.publicaciones.saludcastillayleon.es/10.1080/20008198.2017.1351204.
Ditlevsen DN, Elklit A. The combined effect of gender and age on post traumatic stress disorder: do men and women show differences in the lifespan distribution of the disorder? Ann Gen Psychiatry. 2010;9:32. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/1744-859X-9-32.
Appels R, Eversole K, Feuillet C, Keller B, et al. Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science. 2018. https://doiorg.publicaciones.saludcastillayleon.es/10.1126/science.aar7191.
Acknowledgements
The authors would like to acknowledge support from the Neuroimaging Core, the Clinical Studies and Database Core, the Integrated Bioinformatics and Biostatistics Core, and the Biorepository Core of the Goodman-Luskin Microbiome Center (GLMC) at the University of California Los Angeles, CA, USA for help with data collection, data processing, data storage, dataset creation, and analyses.
Funding
NIH U54 DK123755 (EAM/LC).
Author information
Authors and Affiliations
Contributions
L.A.K. performed the statistical analyses, prepared the figures, and drafted the manuscript; A.C. was a major contributor to the study design and writing of the manuscript; D.M. provided the estrogen data and contributed to the writing of the manuscript; S.M-J. contributed to the study design and writing of the manuscript, V.W.L. contributed to estrogen data analysis and the writing of the manuscript, J.S. contributed to data collection and writing of the manuscript, J.R. contributed to data collection and writing of the manuscript, J.S.L. contributed to the study design and writing of the manuscript, T.D. contributed to the study design and writing of the manuscript, J.P.J. contributed to the study design and writing of the manuscript, B.D.N. contributed to the study design and writing of the manuscript, L.C. procured funding and was a major contributor in the study design and writing of the manuscript, E.A.M. procured funding and was a major contributor in the study design and writing of the manuscript. All authors reviewed the manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
This study was approved by the Institutional Review Board at the University of California, Los Angeles’s Office of Protection for Research Subjects (Nos. 20–000540 and 20–000515). All participants provided written informed consent.
Consent for publication
Not applicable.
Competing interests
A.C. is a research consultant for YAMAHA. E.A.M. is a member of the scientific advisory boards of Danone, Axial Therapeutics, Amare, Mahana Therapeutics, Pendulum, Bloom Biosciences, and APC Microbiome Ireland. L.C. serves as an advisory board member or consultant for Ardelyx, Arena Pharmaceuticals, Bausch Health, Immunic, Ironwood Pharmaceuticals, Inc., Mauna Kea Technologies, and Trellus; and receives grant support from AnX Robotica, Arena Pharmaceuticals, and Ironwood Pharmaceuticals. All other authors declare no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
13293_2025_709_MOESM1_ESM.pdf
Supplementary Material 1. Mean estrogen levels, with standard errors, in males, premenopausal females, and postmenopausal females. Plasma data is provided in units of ng/mL-plasma, while stool data is provided in units of ng/5 mL-stool. E1, estrone; E2, estradiol; E3, estriol: 2OHE1, 2-hydroxyestrone; 2MeOE1, 2-methoxyestrone; 2OHE2, 2-hydroxyestradiol; 2MeOE2, 2-methoxyestradiol; 3MeOE1, 2-hydroxyestrone-3-methyl ether; 4OHE1, 4-hydroxyestrone; 4MeOE1, 4-methoxyestrone; 4MeOE2, 4-methoxyestradiol; 16aOHE1, 16α-hydroxyestrone; 17epiE3, 17-epiestriol; 16ketoE2, 16-ketoestradiol; 16epiE3, 16-epiestriol; Total, summation of all free or totalestrogens and estrogen metabolites, dependent on category.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Kilpatrick, L.A., Church, A., Meriwether, D. et al. Differential brainstem connectivity according to sex and menopausal status in healthy male and female individuals. Biol Sex Differ 16, 25 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13293-025-00709-4
Received:
Accepted:
Published:
DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13293-025-00709-4