Abstract
The adverse effects of hypercortisolism on the human brain have been highlighted in previous studies of Cushing’s disease (CD). However, the relative alterations in regional hypercortisolism in the brain remain unclear. Thus, we investigated regional volumetric alterations in CD patients. We also analyzed the associations between these volumetric changes and clinical characteristics. The study participants comprised of active CD (n = 60), short-term-remitted CD (n = 28), and long-term-remitted CD (n = 32) patients as well as healthy control subjects (n = 66). Gray matter volumes (GMVs) were measured via voxel-based morphometry. The GMVs of substructures were defined using the automated anatomical labeling (AAL) atlas. Trends toward normalization in GMV were found in most brain substructures of CD patients. Different trends, including enlarged, irreversible, and unaffected, were observed in the other subregions, such as the amygdala, thalamus, and caudate. Morphological changes in GMVs after the resolution of hypercortisolism are a complex phenomenon; the characteristics of these changes significantly differ within the brain substructures.
Introduction
Cushing’s disease (CD) is a rare disorder of chronic hypercortisolism exposure caused by an adrenocorticotrophic hormone (ACTH)-secreting pituitary adenoma (1). Hypercortisolism can induce a variety of physical manifestations and psychiatric symptoms in CD. Tumor resection can not only correct hypercortisolism but can also achieve partial recovery of these symptoms (2). Thus, CD could represent a unique human model for examining hypercortisolism-related alterations, reversibility of disease-related changes after the resolution of the condition, and the relationships between these alterations and clinical characteristics (3).
The detrimental effects of hypercortisolism on the brain were first highlighted in early autopsy and pneumoencephalography reports (4, 5). The introduction of MRI enabled a more accurate assessment of brain structures. Previous MRI studies have already demonstrated the relationship between hypercortisolism and whole-brain damage (2, 6). However, the associations between regional structural alterations and hypercortisolism remain controversial. For instance, Starkman et al. demonstrated that the hippocampal volumes of CD patients were smaller than those of normal subjects reported in the literature (7). Their longitudinal study also demonstrated significant recovery of hippocampal volumes 1 year after curative surgery (8). However, another longitudinal study found only an increase in caudate volume 12 months after the resolution of the condition (9). In addition, although volumes of the left amygdala of active CS patients negatively correlated with depression scores, differences could only be found in right (but not left) amygdala volumes when compared with controls (10).
Our total brain volume analysis demonstrated trends of a partially reversible reduction in total gray matter volume (GMV) in CD patients (11, 12). However, previous studies found only partially reversible GMVs reduction in particular subregions, such as the hippocampus (8). Furthermore, segmented GMVs of previous voxel-based morphometry (VBM) studies demonstrated inconsistent results in other subregions, such as the amygdala, frontal lobe, and caudate (2). In addition, several fMRI studies have demonstrated alterations in the frontal lobe and thalamus (13, 14, 15). Nevertheless, partial reversibility of the changes in each brain subregion after the resolution of hypercortisolism remains controversial.
Moreover, long after a biochemical cure, CD patients still have elevated cardiovascular morbidity and mortality, a higher prevalence of psychopathology, and impaired cognitive function (16). Moreover, irreversible brain alterations have long been demonstrated in previous studies of long-term-remitted CD patients (12, 17). Previous studies have tended to explore long-term reversibility after the resolution of hypercortisolism. However, a recent longitudinal study demonstrated significant reversibility of GMVs in short-term-remitted CD patients (3 months after the resolution of hypercortisolism) (17). Our research group also demonstrated trends of potential reversibility in GMV changes in short-term-remitted CD patients (mean remission time of almost 2 months) (11). To our knowledge, there is little information on the similarities and differences among active, short-term-remitted, and long-term-remitted CD patients that have been evaluated by modern, sophisticated analytical tools.
In the present study, based on the VBM technique, we retrospectively evaluated GMVs in active CD, short-term-remitted CD, and long-term-remitted CD patients and healthy controls (HCs) in greater detail using the automated anatomical labeling (AAL) atlas. We conducted an explorative whole-brain analysis to detect possible structural changes in areas outside these a priori-defined regions of interest (ROI). In addition, we aimed to explore the associations between structural changes and clinical characteristics. Detailed knowledge of the nature of morphological changes of gray matter could greatly help in determining the underlying mechanism of hypercortisolism in the human brain.
Materials and methods
Participants
CD patients were recruited during their routine visits to Rui-Jin Hospital. HCs were recruited via social media and online advertisements and matched to the CD patients based on sex, age, and education. All individuals provided written informed consent before the study participation. The participants were excluded if they had a history of drug and/or alcohol abuse, neurological problems, a history of brain lesion, and contraindications to MRI and are left-handed. A total of 186 study subjects were included, all right-handed with at least 9 years of education. According to the remission duration of CD, we included three groups of CD patients and one group of HCs. These groups comprised of 60 patients with active CD (aCD), 28 patients with short-term-remitted CD (S-rCD) (remission time, 17–135 days; mean ± s.d., 62.14 ± 35.09 days), 32 patients with long-term-remitted CD (L-rCD) (remission time, 198–1472 days; mean ± s.d., 709.75 ± 385.47 days), and 66 HCs. This study was approved by the Institutional Review Board of Rui-Jin Hospital. CD and its etiology were confirmed by an elevated 24-h urinary free cortisol (UFC) level, bilateral petrosal sinus sampling, absence of a blunted circadian rhythm of cortisol secretion, and other clinical features (18). CD remission was confirmed in all remitted CD patients by normal UFC after surgery.
Clinical characteristics
Clinical characteristics, including 24-h UFC levels, serum cortisol and ACTH levels, BMI, weight, blood pressure, and plasma glucose level, were evaluated in all CD patients. Disease duration was calculated from the earliest point at which clinical signs were noted in the patient’s history to the time of cortisol normalization. Remission duration was calculated from the date of surgery. The medical histories of all study subjects were recorded using a standardized questionnaire.
MRI acquisition and processing
The subjects were scanned using a 3.0T Philips Ingenia MRI scanner and a 3.0T GE Signa Excite HD (Chicago, IL, USA) MRI scanner. High-resolution T1-weighted images were obtained using three-dimensional brain volume imaging sequences with the following settings: repetition time, 5.552 ms; echo time, 1.752 ms; field of view, 256 × 256 mm; flip angle, 12°; 1 sagittal acquisition; matrix, 256 × 256; NEX, 1; slices, 96 thickness, 1.0 mm with no interslice gap; and bandwidth, 244.141 Hz.
The Computational Anatomy Toolbox (CAT12; http://dbm.neuro.uni-jena.de/cat12/) in MATLAB (R2012b; MathWorks, Natick, MA, USA) was used to analyze T1-weighted images. We conducted brain extraction and spatial normalization and obtained GMVs after brain segmentation according to the standard template. In addition, we corrected the influence of spatial normalization. We defined the total GMV as the summation of the GMVs of all voxels and analyzed GMVs in the ROIs (i.e. the hippocampus, amygdala, cerebellum, frontal lobe, thalamus, and caudate). To provide a complete view of alterations in the whole brain, we also used the AAL atlas to parcellate the GMV maps into all other brain subregions (ROIs) (19).
Statistical analysis
We analyzed the demographic and clinical characteristics of the four groups. Differences in age and education were determined via one-way ANOVA. Sex differences were determined using the chi-squared test. Among the three groups of CD patients, we employed one-way ANOVA to compare the patients’ BMI, weight, plasma glucose level, systolic pressure, diastolic pressure, serum cortisol levels, 24-h UFC level, ACTH level, and disease duration.
To examine gray matter changes during hypercortisolism exposure and after the resolution of the condition, the four groups’ total GMVs, as well as 116 AAL-based GMVs, were compared via one-way ANOVA. Age, sex, education level, differences in MRI, and total intracranial volume (TIV) were used as covariates. Furthermore, we calculated the partial correlation between each clinical characteristic and each GMV measurement while controlling the abovementioned covariates. All data analyses were conducted using the SPSS software (version 22.0; SPSS). For all statistical tests, a two-sided P-value of 0.05 was used, unless Bonferroni correction for multiple comparisons was indicated. Moreover, we calculated the partial correlation between each GMV measure and each clinical characteristic, including 24-h UFC levels, serum cortisol level, ACTH level, BMI, weight, blood pressure, plasma glucose level, and disease duration.
To investigate the potential reversibility of GMV changes in CD patients, one-way ANOVA of regional GMVs from longitudinal MRI data (including active and short-term remission stages for each CD patient) and data of age-, sex-, and education-matched HC subjects were calculated. A total of 11 patients and 11 HC subjects were included in the preliminary longitudinal study. Partial correlations between these GMVs and clinical characteristics were also calculated.
Results
Demographic and clinical characteristics
The demographic data, clinical characteristics, hormone levels, GMVs, and TIVs of all four groups are presented in Table 1.
Demographics and clinical characteristics of study patients and healthy control subjects.
Characteristics | aCD (n = 60) | S-rCD (n = 28) | L-rCD (n = 32) | HC (n = 66) | P value |
---|---|---|---|---|---|
Current age (years) | 42.28 ± 11.49 | 38.89 ± 14.59 | 44.50 ± 11.43 | 45.32 ± 14.59 | 0.153a |
Gender (male/female) | 10/50 | 5/23 | 3/29 | 18/48 | 0.181b |
Education (years) | 12.67 ± 3.07 | 12.29 ± 2.75 | 11.34 ± 3.62 | 13.02 ± 3.02 | 0.092a |
BMI | 26.09 ± 4.85 | 26.23 ± 4.13 | 23.73 ± 3.54 | – | 0.032a, d |
Weight (kg) | 68.84 ± 14.33 | 69.49 ± 11.18 | 60.73 ± 9.88 | – | 0.007a, e |
Plasma glucose (mmol/L) | 6.56 ± 2.29 | 5.46 ± 1.18 | 5.12 ± 1.23 | – | 0.001a, e |
Systolic pressure (mmHg) | 143.88 ± 19.57 | 129.39 ± 21.14 | 122.25 ± 20.36 | – | <0.001a, f |
Diastolic pressure (mmHg) | 90.78 ± 13.40 | 86.86 ± 13.17 | 81.59 ± 13.90 | – | 0.009a,e |
Estimated illness duration (months) | 78.78 ± 108.86 | 55.86 ± 68.05 | 59.44 ± 56.23 | – | 0.429a |
Remission time (days) | – | 62.14 ± 35.09 | 709.75 ± 385.47 | – | <0.001c, f |
Hydrocortisone replacement (mg) | – | 39.21 ± 18.49 | 16.41 ± 16.63 | – | <0.001c, f |
Urinary free cortisol (μg/24 h) | 781.49 ± 869.91 | 120.82 ± 125.28 | 103.93 ± 142.13 | – | <0.001a, f |
Plasma cortisol (08:00 h) (μg/dL) | 26.54 ± 15.24 | 2.46 ± 4.10 | 4.29 ± 7.22 | – | <0.001a, f |
Plasma cortisol (16:00 h) (μg/dL) | 21.85 ± 10.54 | 4.87 ± 4.15 | 5.31 ± 5.32 | – | <0.001a, f |
Plasma cortisol (00:00 h) (μg/dL) | 21.56 ± 14.54 | 3.82 ± 3.80 | 4.20 ± 5.66 | – | <0.001a, f |
Plasma ACTH (pg/mL) | 98.45 ± 69.03 | 14.31 ± 14.59 | 15.78 ± 16.89 | – | <0.001a, f |
Gray matter volumes (GMVs) (cm3) | 564.79 ± 30.01 | 585.69 ± 30.09 | 602.46 ± 31.23 | 614.14 ± 31.12 | <0.001a, f |
Total intracranial volumes (cm3) | 1387.68 ± 95.32 | 1421.99 ± 97.09 | 1409.68 ± 100.78 | 1446.38 ± 98.91 | 0.011a, d |
BMI | 26.09 ± 4.85 | 26.23 ± 4.13 | 23.73 ± 3.54 | – | 0.032a, d |
Weight (kg) | 68.84 ± 14.33 | 69.49 ± 11.18 | 60.73 ± 9.88 | – | 0.007a, e |
Serum glucose (mmol/L) | 6.56 ± 2.29 | 5.46 ± 1.18 | 5.12 ± 1.23 | – | 0.001a, e |
Systolic pressure (mmHg) | 143.88 ± 19.57 | 129.39 ± 21.14 | 122.25 ± 20.36 | – | <0.001a, f |
Diastolic pressure (mmHg) | 90.78 ± 13.40 | 86.86 ± 13.17 | 81.59 ± 13.90 | – | 0.009a, e |
Estimated illness duration (months) | 78.78 ± 108.86 | 55.86 ± 68.05 | 59.44 ± 56.23 | – | 0.429a |
Remission time (days) | – | 62.14 ± 35.09 | 709.75 ± 385.47 | – | <0.001c, f |
Hydrocortisone replacement (mg) | – | 39.21 ± 18.49 | 16.41 ± 16.63 | – | <0.001c, f |
Urinary free cortisol (μg/24 h) | 781.49 ± 869.91 | 120.82 ± 125.28 | 103.93 ± 142.13 | – | <0.001a, f |
Plasma cortisol (08:00 h) (μg/dL) | 26.54 ± 15.24 | 2.46 ± 4.10 | 4.29 ± 7.22 | – | <0.001a, f |
Plasma cortisol (16:00 h) (μg/dL) | 21.85 ± 10.54 | 4.87 ± 4.15 | 5.31 ± 5.32 | – | <0.001a, f |
Plasma cortisol (00:00 h) (μg/dL) | 21.56 ± 14.54 | 3.82 ± 3.80 | 4.20 ± 5.66 | – | <0.001a, f |
Plasma ACTH (pg/mL) | 98.45 ± 69.03 | 14.31 ± 14.59 | 15.78 ± 16.89 | – | <0.001a, f |
GMVs (cm3) | 564.79 ± 30.01 | 585.69 ± 30.09 | 602.46 ± 31.23 | 614.14 ± 31.12 | <0.001a, f |
Total intracranial volumes (cm3) | 1387.68 ± 95.32 | 1421.99 ± 97.09 | 1409.68 ± 100.78 | 1446.38 ± 98.91 | 0.011a, d |
BMI | 26.09 ± 4.85 | 26.23 ± 4.13 | 23.73 ± 3.54 | – | 0.032a, d |
Weight (kg) | 68.84 ± 14.33 | 69.49 ± 11.18 | 60.73 ± 9.88 | – | 0.007a, e |
Serum glucose (mmol/L) | 6.56 ± 2.29 | 5.46 ± 1.18 | 5.12 ± 1.23 | – | 0.001a, e |
Systolic pressure (mmHg) | 143.88 ± 19.57 | 129.39 ± 21.14 | 122.25 ± 20.36 | – | <0.001a, f |
Diastolic pressure (mmHg) | 90.78 ± 13.40 | 86.86 ± 13.17 | 81.59 ± 13.90 | – | 0.009a, e |
Values are expressed as mean ± s.d.
aThe P value for current age, education, BMI, weight, serum glucose level, systolic pressure, diastolic pressure, disease duration, and hormone level among groups was obtained by one-way ANOVA; bThe P value for gender distribution in the four groups was obtained by chi-square test; cThe P value for remission time and dose of hydrocortisone replacement between the two rCD groups was obtained by two sample t-test; d P < 0.05; e P < 0.01; f P < 0.001.
aCD, patient group with active CD; HC, healthy control group; L-rCD, patient group with long-term remitted CD; S-rCD, patient group with short-term remitted CD.
No significant differences were observed among the four groups in terms of age, sex, or education level. Moreover, there were no significant differences among the three CD groups in terms of disease duration.
All adenomas of CD patients were smaller than 10 mm. All CD patients were treated with transsphenoidal surgery without radiotherapy or bilateral adrenalectomy. About 8 of 60 active CD patients and 1 short-term-remitted CD patient had previously undergone surgery. One long-term-remitted CD patient used steroidogenesis inhibitors to achieve eucortisolism after surgery. CD remission was confirmed by normal UFC and adrenal insufficiency in all remitted CD patients following surgery or medical treatment.
Compared with the two remitted CD groups, significantly higher BMI, body weight, plasma glucose level, blood pressure, 24-h UFC level, serum cortisol level, and adrenocorticotropin level were observed in active CD patients. All short-term-remitted CD patients and 21 of all 32 long-term-remitted CD patients were substituted with hydrocortisone.
ANOVA and correlation analysis of total GMVs
The total GMV analysis of this study was similar to that of our previous total brain volume study, which included some of the participants of the current study (12). A significant difference (P < 0.001) was observed in the current study’s ANOVA of total GMV. HCs had the largest total GMV (614.14 ± 31.12 mm3) among the four groups, whereas active CD patients had the smallest total GM volume (564.79 ± 30.01 mm3) among the four groups. The total GM volumes of the two remitted CD groups (585.69 ± 30.09 mm3 and 602.46 ± 31.23 mm3) were between those of the HC group and the active CD group. Furthermore, significant and positive correlations were found between the total GMVs and remission time in the short-term-remitted CD group (R = 0.591, P = 0.003). These results indicated that trends of reversible changes in the total GMV occurred shortly after the resolution of hypercortisolism (Fig. 1 and Table 1).
The trends of GMV changes in ROIs
ANOVA and post hoc analysis of GM volumes among four groups were performed. The correlation analysis between remission time and GM volumes was also performed in two remitted CD groups. Four patterns could be found in brain subregions under chronic exposure to hypercortisolism. The details of the differences in the GMVs of ROIs among the four groups are presented in Table 2 and in Supplementary Tables 1, 2 and 3 (see section on supplementary materials given at the end of this article).
GM volumes of ROIs among study patients and healthy control subjects (samples).
AAL template | Regions of interest | aCD (n = 60) | S-rCD (n = 28) | L-rCD (n = 32) | HC (n = 66) | P values |
---|---|---|---|---|---|---|
Type A: trend toward normalization (number of ROIs = 87) | ||||||
003_Frontal_Sup_L | Left superior frontal gyrus | 7.368 ± 0.674 | 7.701 ± 0.677 | 7.971 ± 0.701 | 8.189 ± 0.699 | <0.001a |
Post hoc analysis | aCD << S-rCD | aCD <<< L-rCD | aCD <<< HC | S-rCD = L-rCD | S-rCD << HC | L-rCD < HC |
P = 0.008b | P < 0.001a | P < 0.001a | P = 0.334 | P = 0.001b | P = 0.039c | |
004_Frontal_Sup_R | Right superior frontal gyrus | 8.244 ± 0.736 | 8.628 ± 0.741 | 9.037 ± 0.764 | 9.155 ± 0.745 | <0.001a |
Post hoc analysis | aCD << S-rCD | aCD <<< L-rCD | aCD <<< HC | S-rCD = L-rCD | S-rCD << HC | L-rCD = HC |
P = 0.006b | P < 0.001a | P < 0.001a | P = 0.182 | P = 0.002b | P = 0.120 | |
037_Hippocampus_L | Left hippocampus | 3.204 ± 0.325 | 3.227 ± 0.328 | 3.319 ± 0.339 | 3.460 ± 0.341 | <0.001a |
Post hoc analysis | aCD = S-rCD | aCD = L-rCD | aCD <<< HC | S-rCD = L-rCD | S-rCD << HC | L-rCD < HC |
P = 0.398 | P = 0.064 | P < 0.001a | P = 0.412 | P = 0.001b | P = 0.021c | |
038_Hippocampus_R | Right hippocampus | 3.300 ± 0.313 | 3.336 ± 0.315 | 3.419 ± 0.239 | 3.508 ± 0.251 | 0.008b |
Post hoc analysis | aCD = S-rCD | aCD = L-rCD | aCD <<< HC | S-rCD = L-rCD | S-rCD < HC | L-rCD = HC |
P = 0.188 | P = 0.063 | P < 0.001a | P = 0.671 | P = 0.034c | P = 0.105 | |
Type B: trends for enlargement (number of ROIs = 8) | ||||||
071_Caudate_L | Left caudate | 3.470 ± 0.558 | 3.578 ± 0.556 | 3.929 ± 0.577 | 3.499 ± 0.545 | 0.002b |
Post hoc analysis | aCD = S-rCD | aCD <<< L-rCD | aCD = HC | S-rCD < L-rCD | S-rCD = HC | L-rCD > HC |
P = 0.173 | P < 0.001a | P = 0.134 | P = 0.043c | P = 0.858 | P = 0.021c | |
072_Caudate_R | Right caudate | 3.615 ± 0.465 | 3.805 ± 0.460 | 3.988 ± 0.481 | 3.648 ± 0.479 | 0.002b |
Post hoc analysis | aCD < S-rCD | aCD <<< L-rCD | aCD = HC | S-rCD = L-rCD | S-rCD = HC | L-rCD > HC |
P = 0.024c | P < 0.001a | P = 0.089 | P = 0.242 | P = 0.351 | P = 0.024c | |
Type C: trends of irreversible changes (number of ROIs = 10) | ||||||
077_Thalamus_L | Left thalamus | 3.908 ± 0.457 | 3.846 ± 0.460 | 4.028 ± 0.481 | 4.206 ± 0.479 | 0.001b |
Post hoc analysis | aCD = S-rCD | aCD = L-rCD | aCD <<< HC | S-rCD = L-rCD | S-rCD <<< HC | L-rCD < HC |
P = 0.931 | P = 0.153 | P < 0.001a | P = 0.206 | P < 0.001a | P = 0.036c | |
078_Thalamus_R | Right thalamus | 4.013 ± 0.534 | 4.043 ± 0.512 | 4.295 ± 0.560 | 4.437 ± 0.552 | <0.001a |
Post hoc analysis | aCD = S-rCD | aCD < L-rCD | aCD <<< HC | S-rCD = L-rCD | S-rCD << HC | L-rCD = HC |
P = 0.572 | P = 0.012c | P < 0.001a | P = 0.105 | P = 0.001b | P = 0.145 | |
Type D: trends of an uninfluenced (number of ROIs = 11) | ||||||
041_Amygdala_L | Left amygdala | 1.050 ± 0.085 | 1.093 ± 0.085 | 1.075 ± 0.085 | 1.070 ± 0.089 | 0.139 |
Post hoc analysis | aCD << S-rCD | aCD = L-rCD | aCD <<< HC | S-rCD = L-rCD | S-rCD = HC | L-rCD = HC |
P = 0.007b | P = 0.094 | P = 0.006b | P = 0.354 | P = 0.593 | P = 0.587 | |
042_Amygdala_R | Right amygdala | 1.137 ± 0.077 | 1.156 ± 0.074 | 1.158 ± 0.079 | 1.152 ± 0.081 | 0.511 |
Post hoc analysis | aCD = S-rCD | aCD = L-rCD | aCD <<< HC | S-rCD = L-rCD | S-rCD = HC | L-rCD = HC |
P = 0.072 | P = 0.109 | P = 0.004b | P = 0.837 | P = 0.638 | P = 0.482 |
Values are expressed as mean ± s.d. (cm3). The P value among groups were obtained by one-way ANOVA and post hoc analysis.
a P < 0.001; b P < 0.01; c P < 0.05.
aCD, patient group with active CD; HC, healthy control group; L-rCD, patient group with long-term remitted CD; S-rCD, patient group with short-term remitted CD.
Type A trends toward normalization
In most brain subregions (87/116), the trends for GMV change after the resolution of hypercortisolism were similar to those of the total GMVs. The aCD group had the smallest GMV among the four groups, whereas the HC group had the largest GMV. The GMVs of the two remitted CD groups were between those of the aCD patients and HCs. These results indicated that trends for reversibility of regional GMVs also occurred after the resolution of hypercortisolism.
Moreover, in 83 subregions of type A, the GMVs were significantly or nearly significantly correlated with the remission time in the short-term-remitted CD group, indicating that these trends occurred in the short term following curative surgery. In some other subregions, such as the hippocampus, significant correlations between remission time and GMVs could only be found in the long-term-remitted CD group, indicating that the trends for recovery occurred in the late stage following curative surgery in these patients (Fig. 1, Types A1 and A2).
Type B trends for enlargement
In eight subregions (bilateral caudate, etc.), the GMVs of the two remitted CD groups were not only larger than those of the active CD patients but also larger than those of the HC subjects. No correlations with remission time were observed between the remitted CD groups. These results may indicate that, after the resolution of hypercortisolism, trends for increased GMVs might occur in these subregions (Fig. 1, Type B; Supplementary Table 3).
Type C trends of irreversible changes
In the bilateral thalamus and several subregions of the cerebellum, the GMVs of the HC group were significantly larger than those of any CD group. These results may indicate that after the resolution of hypercortisolism, no significant recovery of GMVs could be found in these regions. Thus, the reductions in GMVs in these regions seem to be irreversible (Fig. 1, Type C; Supplementary Table 3).
Type D trends of unaffected
In 11 subregions (bilateral amygdala, etc.), no significant difference in GMVs could be found among the four groups. In addition, no correlations were observed between remission time and GMVs in either of the remitted CD groups. These findings may indicate that neither hypercortisolism nor the resolution of hypercortisolism could influence these subregions (Fig. 1, Type D; Supplementary Table 3).
Preliminary longitudinal analysis
Patterns similar to those found via cross-sectional analysis could be found via one-way ANOVA. However, a significant and positive correlation with remission time was observed only in AAL 079 (Heschl left). The details of one-way ANOVA and correlation analysis of longitudinal data are presented in Supplementary Fig. 3.
Correlation analysis between GMVs and other clinical features in brain subregions
Negative correlations between disease duration and the GMVs of subregions were observed in the right operculum of the inferior frontal gyrus, bilateral triangle of the inferior frontal gyrus, right rectus, bilateral insula, left hippocampus, left middle occipital gyrus, left inferior occipital gyrus, bilateral postcentral gyrus, right inferior parietal gyrus, right supramarginal gyrus, bilateral angular gyrus, and left middle temple gyrus in long-term-remitted CD patients. Interestingly, in the bilateral caudate, positive correlations were observed between disease duration and GMV in CD patients with long-term remission, indicating that the duration of hypercortisolism exposure is related to the increase in GMV in the caudate. The results of correlated analysis with disease durations are presented in Supplementary Table 4. The GMVs of 49 subregions were found to be significantly correlated with the cortisol level in CD patients. The results of the correlation analysis with the level of cortisol are presented in Supplementary Table 5.
Discussion
Major finding
After the resolution of hypercortisolism, trends of possibly reversible alterations in GMV were observed in most brain subregions of CD patients. However, in other subregions, several different trends of GMV change were found. These results may indicate that, under chronic hypercortisolism exposure and after the resolution of the condition, the regional change in GMV in CD patients was a complex phenomenon, not mere partial reversibility.
Previous studies on hypercortisolism-related brain alterations
The limbic system plays an essential role in stress-related disorders, for example major depression and posttraumatic stress disorder (PTSD). It is rich in receptors for cortisol, especially mineralocorticoid receptors (MRs), to which both mineralocorticoids and glucocorticoids can bind (20). Experimental animal studies have already demonstrated the neurotoxic effects of hypercortisolism on the limbic system (21). Previous studies of CD patients have focused on alterations in the limbic system, especially in the hippocampus, associated with depression and anxiety (2, 6, 22, 23, 24). However, alteration of the limbic system cannot explain other factors, including the psychological problem of chronic exposure to hypercortisolism, cognitive impairment, and even dementia (25). Contrary to the limited distribution of MRs, glucocorticoid receptors (GRs), another type of cortisol receptor, are widely distributed throughout the brain (21). Glucocorticoids affect all cell types in the brain, including non-neuronal cells, such as microglia, oligodendrocytes, astrocytes, and endothelial cells (26). Several recent studies of CD patients demonstrated widely distributed alterations in structure, function, and blood supply as well as diffusion in the human brain (3, 17, 27). Thus, the widespread alterations observed in the current study are not only in line with those previous studies and our previous total GM study (12) but are also in line with the different patterns of regional morphological changes observed in the current study. This indicates that the influence of hypercortisolism on the human brain is a complex phenomenon.
Trends in GMV changes after the resolution of hypercortisolism
Compared with the GMVs of the HC group, significantly smaller GMVs were observed in most brain subregions of active CD patients. ANOVA analysis revealed trends of reversibility in these regions. In addition, significant and positive correlations were observed between GMVs in these subregions and remission time in short-term-remitted CD patients. Previous studies demonstrated not only partially reversible GMV reductions in particular brain subregions, such as the hippocampus (8, 28), but also similar partial reversibility in other brain subregions and in the whole brain (17, 29). Thus, it is reasonable to assume that partial reversibility of GMVs could be found, not only in particular brain subregions but also in most brain subregions of remitted CD patients.
Moreover, we found a significant and positive correlation between remission time and GMVs of the hippocampus and adjacent structures in CD patients with long-term remission. The hippocampus and adjacent regions are rich in GRs compared with the number of GRs in any other subregions of the human brain (21). Thus, these structures exhibit increased sensitivity to cortisol, affecting both volume loss and recovery (8, 17). Previous studies of stress-related disorders and neurodegenerative disorders have often demonstrated alterations in hippocampal volumes (30). However, our VBM study of short-term-remitted CD patients did not find any difference among active CD patients, short-term-remitted CD patients, and HCs in the hippocampus (11). Contrarily, significant alterations could be observed in the hippocampus among the three patient groups when compared with the hippocampus of HC subjects. Thus, as the hippocampus has an abundance of cortisol receptors, the trends of the hippocampal recovery process after hypercortisolism resolution might be quite different from the trends in other brain regions.
The GMVs in eight subregions (the bilateral supply motor area, bilateral caudate, bilateral paracentral lobule, and bilateral putamen) of the two remitted CD groups were not only larger than those of the active CD patients but also larger than those of the HC subjects. About 12 months after the curative operation, Starkman et al. found that an increase in the volume of the right caudate was significantly associated with decreases in depression and anxiety scores (9). However, without including any HC subjects, it was impossible to draw any further conclusions about whether caudate volume decreased in the active stage of CD. Contrarily, our previous VBM study found greater GMVs in the bilateral caudate of CD patients in short-term remission than in HC subjects (11). Moreover, we found a positive correlation between disease duration and GMVs in these regions of CD patients in long-term remission. According to a recently proposed theory, 'The Scaffolding Theory of Aging and Cognition', the process of neurodegeneration is characterized by both decline and preservation (31). This means that when reductions of GMV occurred in some brain regions, structural compensation could be observed in other regions (32, 33, 34). Nevertheless, it is likely that the cause of larger GMVs in these subregions is multifactorial, including intrinsic imperfections of surgery and endocrine replacement therapy. However, the possibility of a compensatory strategy after the resolution of hypercortisolism cannot be entirely neglected.
In the bilateral thalamus, the GMVs of the HC group were greater than those of any CD group, whereas no significant differences could be found between the active group and the two remitted CD groups. Irreversible effects on the human brain caused by chronic exposure to hypercortisolism have long been observed in CD patients with long-term remission (2). Our findings could indicate that the detriment induced by hypercortisolism might be irreversible in these regions.
The GMVs of several subregions, including the bilateral amygdala, demonstrated no difference among the four groups. Although previous studies revealed functional and structural alterations in the amygdala of adolescents with hypercortisolemia (14), no differences were found in amygdala activation of adult remitted CD patients (2). In structural MRI studies, Andela et al. did not find any volume difference in the bilateral amygdala (35). Santos et al. only found that the volumes of the right (but not the left) amygdala of active CS patients were barely smaller (P = 0.043) than those of HCs (10). Our findings could indicate that hypercortisolism has at least some structural effect on some brain subregions.
Limitations
This analysis of altered GMVs in the brain subregions of CD patients has some limitations. Our study did not include neuropsychological tests or clinical data of HC subjects. These limitations and the cross-sectional design of our study prevent the determination of further conclusions with regard to alterations in regional GMVs caused by hypercortisolism, resolution of the condition, and influence of neuropsychological disorders on these effects. Thus, prospective functional MRI and structural MRI studies that include psychological measures and neuropsychological tests of CD patients and HC subjects are required.
Conclusion
Cross-sectional analysis revealed significantly altered GMVs and trends of reversibility in most brain subregions of CD patients. However, different trends of GMVs were observed in other brain subregions of these patients. These results indicated that regional morphological changes in the gray matter after the resolution of hypercortisolism in CD were due to complex phenomena. Our study enhances our understanding of the structural alterations in the human brain due to hypercortisolism (12). Furthermore, our study may facilitate the development of future studies on the influence of hypercortisolism, with comprehensive neuropsychological tests.
Supplementary materials
This is linked to the online version of the paper at https://doi.org/10.1530/EC-21-0385.
Declaration of interest
The authors declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research reported.
Funding
This study was supported by grants from the National Natural Science Foundation of China General Projects (81270856, 81770779, 81501467, and 61502059), the China Postdoctoral Science Foundation funded project (2016M592656), and the Sichuan Provence Science and Technology Program (2018JY0272).
Data availability statement
The data that support the findings of this study are available on request from the corresponding author, but they are not publicly available due to privacy restrictions.
Author contribution statement
Hong Jiang and Wen-Jie Yang contributed equally to this work. Liu-Guan Bian and Chang Liu contributed equally to this work.
Acknowledgements
The authors gratefully acknowledge the financial supports by the National Science Foundation of China. Useful suggestions given by department of Endocrinology are also acknowledged.
References
- 1↑
Newell-Price J, Bertagna X, Grossman AB, Nieman LK. Cushing’s syndrome. Lancet 2006 367 1605–1617. (https://doi.org/10.1016/S0140-6736(0668699-6)
- 2↑
Andela CD, Van Haalen FM, Ragnarsson O, Papakokkinou E, Johannsson G, Santos A, Webb SM, Biermasz NR, Van Der Wee NJA, Pereira AM. Cushing’s syndrome causes irreversible effects on the human brain: a systematic review of structural and functional magnetic resonance imaging studies.European Journal of Endocrinology 2015 173 R1–R14.(https://doi.org/10.1530/EJE-14-1101)
- 3↑
van der Werff SJA, Pannekoek JN, Andela CD, Meijer OC, van Buchem MA, Rombouts SARB, van der Mast RC, Biermasz NR, Pereira AM, van der Wee NJA. Resting-state functional connectivity in patients with long-term remission of Cushing’s disease. Neuropsychopharmacology 2015 40 1888–1898. (https://doi.org/10.1038/npp.2015.38)
- 4↑
Momose KJ, Kjellberg RN, Kliman B. High incidence of cortical atrophy of the cerebral and cerebellar hemispheres in Cushing’s disease. Radiology 1971 99 341–348. (https://doi.org/10.1148/99.2.341)
- 5↑
Trethowan WH, Cobb S. Neuropsychiatric aspects of Cushing’s syndrome. AMA Archives of Neurology and Psychiatry 1952 67 283–309. (https://doi.org/10.1001/archneurpsyc.1952.02320150016002)
- 6↑
Chen Y, Zhang J, Tan H, Li J, Yu Y. Detrimental effects of hypercortisolism on brain structure and related risk factors. Scientific Reports 2020 10 1–12. (https://doi.org/10.1038/s41598-020-68166-0)
- 7↑
Starkman MN, Gebarski SS, Berent S, Schteingart DE. Hippocampal formation volume, memory dysfunction, and Cortisoi levels in patients with Cushing’s syndrome. Biological Psychiatry 1992 32 756–765. (https://doi.org/10.1016/0006-3223(9290079-f)
- 8↑
Starkman MN, Giordani B, Gebarski SS, Berent S, Schork MA, Schteingart DE. Decrease in cortisol reverses human hippocampal atrophy following treatment of Cushing’s disease. Biological Psychiatry 1999 46 1595–1602. (https://doi.org/10.1016/s0006-3223(9900203-6)
- 9↑
Starkman MN, Giordani B, Gebarski SS, Schteingart DE. Improvement in mood and ideation associated with increase in right caudate volume. Journal of Affective Disorders 2007 101 139–147. (https://doi.org/10.1016/j.jad.2006.11.007)
- 10↑
Santos A, Granell E, Gómez-Ansón B, Crespo I, Pires P, Vives-Gilabert Y, Valassi E, Webb SM, Resmini E. Depression and anxiety scores are associated with amygdala volume in Cushing’s syndrome: preliminary study. BioMed Research International 2017 2017 2061935. (https://doi.org/10.1155/2017/2061935)
- 11↑
Jiang H, Ren J, He NY, Liu C, Sun YH, Jian FF, Bian LG, Shen JK, Yan FH & Pan SJ et al.Volumetric magnetic resonance imaging analysis in patients with short-term remission of Cushing’s disease. Clinical Endocrinology 2017 87 367–374. (https://doi.org/10.1111/cen.13381)
- 12↑
Jiang H, Liu C, Pan SJ, Ren J, He NY, Sun YH, Bian LG, Yan FH, Yang WJ, Sun QF. Reversible and the irreversible structural alterations on brain after resolution of hypercortisolism in Cushing’s disease. Steroids 2019 151 108457. (https://doi.org/10.1016/j.steroids.2019.108457)
- 13↑
Crespo I, Esther GM, Santos A, Valassi E, Yolanda VG, De Juan-Delago M, Webb SM, Gõmez-Ansõn B, Resmini E. Impaired decision-making and selective cortical frontal thinning in Cushing’s syndrome. Clinical Endocrinology 2014 81 826–833. (https://doi.org/10.1111/cen.12564)
- 14↑
Maheu FS, Mazzone L, Merke DP, Keil MF, Stratakis CA, Pine DS, Ernst M. Altered amygdala and hippocampus function in adolescents with hypercortisolemia: a functional magnetic resonance imaging study of Cushing syndrome. Development and Psychopathology 2008 20 1177–1189. (https://doi.org/10.1017/S0954579408000564)
- 15↑
Bas-Hoogendam JM, Andela CD, van der Werff SJA, Pannekoek JN, van Steenbergen H, Meijer OC, van Buchem MA, Rombouts SARB, van der Mast RC & Biermasz NR et al.Altered neural processing of emotional faces in remitted Cushing’s disease. Psychoneuroendocrinology 2015 59 134–146. (https://doi.org/10.1016/j.psyneuen.2015.05.001)
- 16↑
Pivonello R, De Martino MC, De Leo M, Simeoli C, Colao A. Cushing’s disease: the burden of illness. Endocrine 2017 56 10–18. (https://doi.org/10.1007/s12020-016-0984-8)
- 17↑
Hou B, Gao L, Shi L, Luo Y, Guo X, Young GS, Qin L, Zhu H, Lu L & Wang Z et al.Reversibility of impaired brain structures after transsphenoidal surgery in Cushing’s disease: a longitudinal study based on an artificial intelligence–assisted tool. Journal of Neurological Surgery 2020 [epub]. (https://doi.org/10.3171/2019.10.JNS191400)
- 18↑
Pivonello R, De Leo M, Cozzolino A, Colao A. The treatment of Cushing’s disease. Endocrine Reviews 2015 36 385–486. (https://doi.org/10.1210/er.2013-1048)
- 19↑
Tzourio-Mazoyer N, Landeau B, Papathanassiou D, Crivello F, Etard O, Delcroix N, Mazoyer B, Joliot M. Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage 2002 15 273–289. (https://doi.org/10.1006/nimg.2001.0978)
- 20↑
De Kloet ER, Joëls M, Holsboer F. Stress and the brain: from adaptation to disease. Nature Reviews: Neuroscience 2005 6 463–475. (https://doi.org/10.1038/nrn1683)
- 21↑
De Kloet ER, Vreugdenhil E, Oitzl MS, Joëls M. Brain corticosteroid receptor balance in health and disease. Endocrine Reviews 1998 19 269–301. (https://doi.org/10.1210/edrv.19.3.0331)
- 22↑
Bauduin SEEC, Van Der Wee NJA, Van Der Werff SJA. Structural brain abnormalities in Cushing’s syndrome.Current Opinion in Endocrinology, Diabetes and Obesity 2018 25 285–289.(https://doi.org/10.1097/MED.0000000000000414)
- 23↑
Gao L, Liu L, Shi L, Luo Y, Wang Z, Guo X, Xing B. Dynamic changes of views on the brain changes of Cushing’s syndrome using different computer-assisted tool. Reviews in Endocrine and Metabolic Disorders 2020 21 185–200. (https://doi.org/10.1007/s11154-020-09540-1)
- 24↑
Resmini E, Santos A, Webb SM. Cortisol excess and the brain. Frontiers of Hormone Research 2016 46 74–86. (https://doi.org/10.1159/000443868)
- 25↑
Forget H, Lacroix A, Bourdeau I, Cohen H. Long-term cognitive effects of glucocorticoid excess in Cushing’s syndrome. Psychoneuroendocrinology 2016 65 26–33. (https://doi.org/10.1016/j.psyneuen.2015.11.020)
- 26↑
Mifsud KR, Reul JMHM. Mineralocorticoid and glucocorticoid receptor-mediated control of genomic responses to stress in the brain. Stress 2018 21 389–402. (https://doi.org/10.1080/10253890.2018.1456526)
- 27↑
van der Werff SJA, Andela CD, Pannekoek JN, Meijer OC, van Buchem MA, Rombouts SARB, van der Mast RC, Biermasz NR, Pereira AM, van der Wee NJA. Widespread reductions of white matter integrity in patients with long-term remission of Cushing’s disease. NeuroImage: Clinical 2014 4 659–667. (https://doi.org/10.1016/j.nicl.2014.01.017)
- 28↑
Starkman MN, Giordani B, Gebarski SS, Schteingart DE. Improvement in learning associated with increase in hippocampal formation volume. Biological Psychiatry 2003 53 233–238. (https://doi.org/10.1016/s0006-3223(0201750-x)
- 29↑
Tirosh A, RaviPrakash H, Papadakis GZ, Tatsi C, Belyavskaya E, Charalampos L, Lodish MB, Bagci U, Stratakis CA. Computerized analysis of brain MRI parameters dynamics in young patients with Cushing syndrome – a case-control study. Journal of Clinical Endocrinology and Metabolism 2020 105 e2069–e2077. (https://doi.org/10.1210/clinem/dgz303)
- 30↑
Purnell JQ, Klopfenstein BA, Stevens AA, Havel PJ, Adams SH, Dunn TN, Krisky C, Rooney WD. Brain functional magnetic resonance imaging response to glucose and fructose infusions in humans. Diabetes, Obesity and Metabolism 2011 13 229–234. (https://doi.org/10.1111/j.1463-1326.2010.01340.x)
- 31↑
Goh JO, Park DC. Neuroplasticity and cognitive aging: the scaffolding theory of aging and cognition. Restorative Neurology and Neuroscience 2009 27 391–403. (https://doi.org/10.3233/RNN-2009-0493)
- 32↑
Zhang H, Lee A, Qiu A. A posterior-to-anterior shift of brain functional dynamics in aging. Brain Structure and Function 2017 222 3665–3676. (https://doi.org/10.1007/s00429-017-1425-z)
- 33↑
Park DC, Reuter-Lorenz P. The adaptive brain: aging and neurocognitive scaffolding. Annual Review of Psychology 2009 60 173–196. (https://doi.org/10.1146/annurev.psych.59.103006.093656)
- 34↑
Fuchs E, Flügge G. Adult neuroplasticity: more than 40 years of research. Neural Plastcity 2014 2014 541870.(https://doi.org/10.1155/2014/541870)
- 35↑
Andela CD, Van Der Werff SJA, Pannekoek JN, Van Den Berg SM, Meijer OC, Buchem MA, Rombouts SARB, Van Der Mast RC, Romijn JA & Tiemensma J et al.Smaller grey matter volumes in the anterior cingulate cortex and greater cerebellar volumes in patients with long-term remission of Cushing’s disease: a case-control study. European Journal of Endocrinology 2013 169 811–819. (https://doi.org/10.1530/EJE-13-0471)