A Resting-state Functional Magnetic Resonance Imaging Study of Whole-brain Functional Connectivity of Voxel Levels in Patients With Irritable Bowel Syndrome With Depressive Symptoms.

BACKGROUND/AIMS: Depressive symptom is one of the most common symptoms in patients with irritable bowel syndrome (IBS), but its pathogenetic mechanisms remain unclear. As a voxel-level graph theory analysis method, degree centrality (DC) can provide a new perspective for exploring the abnormalities of whole-brain functional network of IBS with depressive symptoms (DEP-IBS).
METHODS: DC, voxel-wise image and clinical symptoms correlation and seed-based functional connectivity (FC) analyses were performed in 28 DEP-IBS patients, 21 IBS without depressive symptoms (nDEP-IBS) patients and 36 matched healthy controls (HC) to reveal the abnormalities of whole brain FC in DEP-IBS.
RESULTS: Compared to nDEP-IBS patients and HC, DEP-IBS patients showed significant decrease of DC in the left insula and increase of DC in the left precentral gyrus. The DC's z-scores of the left insula negatively correlated with depression severity in DEP-IBS patients. Compared to nDEP-IBS patients, DEP-IBS patients showed increased left insula-related FC in the left inferior parietal lobule and right inferior occipital gyrus, and decreased left insula-related FC in the left precentral gyrus, right supplementary motor area (SMA), and postcentral gyrus. In DEP-IBS patients, abstracted clusters' mean FC in the right SMA negatively correlated with depressive symptoms.
CONCLUSIONS: DEP-IBS patients have abnormal FC in brain regions associated with the fronto-limbic and sensorimotor networks, especially insula and SMA, which explains the vicious circle between negative emotion and gastrointestinal symptoms in IBS. Identification of such alterations may facilitate earlier and more accurate diagnosis of depression in IBS, and development of effective treatment strategies.

Introduction

Irritable bowel syndrome (IBS) is a highly recurrent chronic functional gastrointestinal disorder (FGID) with the typical clinical symptoms of recurrent abdominal pain or discomfort, changes in stool characteristics and bowel habits.[1] The occurrence of IBS is associated with gastrointestinal motility abnormalities, visceral hypersensitivity, intestinal inflammation, and psychological disorders, which can influence IBS symptoms through brain-gut interaction.[2] Almost half of IBS patients suffer from psychological-psychiatric disorders, such as depression, approximately 20-40% of patients accompany with depressive symptoms, which seriously affect the quality of life and work efficiency of patients.[3,4] It has been reported that the combined application of antidepressants can not only effectively improve the mood of IBS patients, but also play a central analgesic effect and weaken visceral sensitivity.[5] The influence of psychological factors on the intestinal tract may be related to the “brain-gut axis,” which is defined as the bidirectional regulatory system, including central nervous system CNS, neuro-endocrine-immune system, hypothalamic-pituitary-adrenal axis, autonomy nervous system, and enteric nervous system. It plays an important role in regulating gastrointestinal motility, visceral sensitivity, brain-gut peptide secretion, stress response, and cognitive function.[6] Abnormal emotional activities and mental states can change autonomic nervous function through the emotional stress circuit, which mainly include cerebral cortex, limbic system, locus ceruleus, vagus/hypothalamic dorsal nucleus, autonomic nervous system, and enteric nervous system, leading to gastrointestinal dysfunction and visceral perception disorders.[7] Up to now, the neural mechanism of depressive symptoms in IBS is still poorly understood and remains inconclusive. Hence, to further explore the brain function changes, such as neuronal activity in the local brain region and the strength of functional connectivity (FC) between brain regions in IBS patients with specific psychological symptoms (such as depression) will help us understand the pathophysiological mechanism of the brain-gut interaction of FGID. In clinical terms, it can provide a scientific basis for formulating personalized treatment strategies.

A large number of studies have reported that FGID patients have functional abnormalities in the limbic system and sensorimotor networks, which are partly associated with depressive symptoms.[8-13] Resting-state functional magnetic resonance imaging (fMRI) can detect the activities of brain neurons in the resting state, which was one of the effective methods to study the neural network mechanism related to brain-gut interaction. Using the amplitude of low-frequency fluctuations[14,15] and regional homogeneity[16] analysis methods, the changes in the activity of local brain regions in the prefrontal lobe-limbic system, basal ganglia and default mode network were showed in patients with IBS, which partly associated with depressive symptoms. Studies based on seed-point and whole-brain functional connection analysis methods have shown that IBS patients had disconnection of emotions, limbic system, sensorimotor, and default mode networks, which are partly associated with depressive symptoms.[17,18] Our previous study found that patients with IBS with depressive symptoms (DEP-IBS) have abnormalities in the microstructure of brain regions related to pain and emotion in the prefrontal-limbic circuit, which further causes and aggravates depressive symptoms (or gastrointestinal symptoms) through the intermediation of gastrointestinal symptoms (or depressive symptoms).[19] It is speculated from the results of the above studies that some regions in the frontal lobe-limbic system is involved not only in the perception of visceral pain, but also in the emotional experience triggered by various stresses in the emotional activity system. It can mediate unpleasantness or depression caused by abdominal pain, and it can also affect the perception of abdominal pain in patients with IBS when depression and other emotional disorders occur. Depression can directly cause the activation of the limbic system, thereby affecting the transmission of visceral pain, leading to visceral hyperesthesia, abdominal pain, and related emotional abnormalities. The results help to understand more clearly why the gastrointestinal symptoms of IBS patients with psychological abnormalities are more serious than those of patients without psychological abnormalities.[3] The inclusion of psychosocial factors may be crucial to elucidate the neural network mechanism of IBS.

The above studies focus on local spontaneous brain activity or analyze neural networks between selected brain regions based on a prior assumption. In recent years, graph-based analysis methods applied to the analysis of complex neural networks, such as the attribute analysis of “small-world” brain networks.[20] According to graph theory, the network considered as a pair relationship between system elements, which consists of a set of edges connecting a set of nodes.[21] Network analysis provides a new conceptual framework for the study of network biology in neurodegenerative and psychosocial diseases.[21,22] Among measures based on the voxel level, degree centrality (DC) refers to the number of connections that are strongly correlated with nodes of a given voxel or binary graph. It can analyze the brain at the voxel level and quantify the importance of a node to other parts of the brain, thus avoiding the bias caused by the selection of brain regions based on prior assumptions.[23] The DC method can reflect the topological properties of nodes or brain regions in the whole brain network and is a measure to describe the importance of network nodes.[20,23] In the present study, the DC method for whole-brain FC analysis of voxel-level based on a graph theoretical approach and a secondary seed-based FC analysis were used to comprehensively explore the intrinsic abnormality of the whole-brain functional network in DEP-IBS patient from the node and network attributes. We hypothesized that long-term negative emotions in IBS patients lead to abnormal FC in certain brain areas, especially those involved in emotional, limbic, basal ganglia and sensorimotor networks, resulting in abnormal regulation of visceral pain perception (such as visceral hypersensitivity), further aggravate gastrointestinal symptoms and negative emotions, thus forming a vicious circle.

Materials and Methods

Subjects

Forty-nine right-handed patients with IBS were recruited from the digestive clinic, as well as 36 HC with matched age, sex and educational level, were enrolled from the physical examination center. This study protocol conforms to the ethical guidelines of the 1975 Declaration of Helsinki as reflected in a priori approval by the institution’s human research committee (Approval No. HZNU20160427), and written informed consent was obtained from all participants before data collection.

All patients fulfilled the following criteria: (1) diagnoses were confirmed on the basis of the Rome III symptom criteria,[24] (2) were not taken antidepressant and prokinetic drugs before enrollment, (3) were without any previous or current neurological disorders and head injury, and (4) were right-handed. A complete medical history taking and physical examination, blood profiles, and colonoscopy were performed in all patients. HCs were no history of chronic pain, bowel disturbance, or visceral disease. Other criteria were the same as those applied to the patients. The enrolled IBS patients were further divided based on the presence or absence of depressive symptoms according to the 17-item Hamilton Depression Rating Scale (HAMD-17),[25] forming 2 groups: DEP-IBS (mean ± SD, 12.75 ± 1.71) and IBS with non-depressive symptoms (nDEP-IBS) (mean ± SD, 3.71 ± 1.38). HAMD-17 is the most commonly used scale for clinical assessment of depressive status. The total score is 24 points. Patients with a score higher than 7 were considered depressed. Higher scores indicate more severe depression. Their demographic and clinical information are shown in Table 1.

Table 1

Demographic and Clinical Characteristics of Participants

Variables	DEP-IBS (n = 28)	nDEP-IBS (n = 21)	HC (n = 36)	P-values
Female/Male	11/17	9/12	10/26	0.456a
Age (yr)	36.36 ± 7.31	32.29 ± 9.96	31.67 ± 8.85	0.087b
Education (years)	11.82 ± 3.09	12.95 ± 2.58	13.28 ± 3.32	0.164b
Duration of IBS (months)	20.64 ± 9.01	18.57 ± 4.00	-	0.331c
GSRS	42.11 ± 11.76	23.43 ± 5.42	-	< 0.001c
HAMD-17	12.75 ± 1.71d	3.71 ± 1.38e,f	0.42 ± 0.77	< 0.001b

IBS, irritable bowel syndrome; DEP-IBS, IBS patients with depressive symptoms; nDEP-IBS, IBS patients without depressive symptoms; HC, healthy controls; GSRS, gastrointestinal symptoms rating scale; HAMD, Hamilton depressive rating scale.

Data are shown as mean ± SD.

aχ2; bP-values of one-way ANOVA and post hoc comparisons (DEP-IBS vs HC, dP<0.001; nDEP vs HC, eP<0.001; DEP vs nDEP, fP<0.001); cUnpaired t test.

Magnetic Resonance Imaging Data Acquisition

Brain imaging was performed by a whole-body 3T GE Discovery MR-750 scanner and an 8-channel head coil (GE Healthcare, Waukesha, WI, USA). All subjects were instructed to remain still and awake, close their eyes, and breathe steadily during the image acquisition process. Special non-magnetic headphones and spongy pads were used to reduce noise and limit head movement. Functional images were acquired using blood oxygenation level dependent–gradient echo–echo planar imaging sequence with parameters: flip angle, 90°; repetition time, 2000 msec; echo time, 30 msec; matrix, 64 × 64; field of view, 192 mm × 192 mm; in-plane resolution, 3.44 mm × 3.44 mm; 37 slices covering the whole brain; and slice thickness, 4 mm. The resting-state sequence scanning time was 8 minutes, and 240 volumes were transversely acquired. After the resting-state scan, structural images were also collected for each subject, which acquired using the spoiled gradient-recalled pulse sequence with parameters: flip angle, 8°; TR, 8.16 msec; echo time, 3.18 msec; matrix, 256 × 256; thickness, 1 mm; and 176 sagittal slices covering the whole brain. In addition, routine clinical MRI sequences such as T1-weighted image and T2-weighted-fluid-attanuated inversion recovery (T2-FLAIR) were collected to exclude anatomical abnormalities and organic brain disease.

Magnetic Resonance Imaging Data Processing

Prior to data processing, T1- and T2-FLAIR weighted images of each participant were examined by 2 senior radiologists, and 2 subjects with abnormal signs such as mass lesions and cerebral infarction were excluded.

Data preprocessing was implemented using the toolbox of Data Processing and Analysis for (Resting-State) Brain Imaging (DPABI) (http://www.restfmri.net/forum)[26] in Statistical Parametric Mapping software version 8 (available from URL: http://www.fil.ion.ucl.ac.uk/spm/software/spm8/). The software runs in MATLAB R2013b (The Mathworks, Natick, MA, USA). The Digital Imaging and Communications in Medicine (DICOM) format of MRI data of all subjects was first converted to a standard format (Neuroimaging Informatics Technology Initiative). The image data of the first 10 time points of the resting state scan were eliminated to enable the subjects to adapt to the scanning environment. Slice-timing was performed to eliminate the interval scanning time difference. The head movement signal was corrected and the head movement information of each subject was obtained. Subjects whose head movements were shifted by more than 3 mm in the x, y, and z directions or rotated by more than 3° were excluded (all the subjects met this standard). Gray matter, white matter, and cerebrospinal fluid of each individual structural image were segmented, and then resampled (3 mm × 3 mm × 3 mm) in the Montreal Neurological Institute (MNI) space. Detrending was applied to remove the linear drift of the baseline signal. The data were bandpass filtered (0.01-0.08 Hz) to eliminate physiological noise. The head motion parameters, cerebrospinal fluid signal, white matter signal, and mean value signal of the whole brain were regressed as covariables to improve the accuracy of the results.

Based on preprocessing, the DC was calculated using DPABI software in a voxel-wise manner with a threshold r > 0.25.[23] Based on the graph theory method, each voxel in the whole brain was defined as a node, and the linear correlation of neural activity signals between voxels was regarded as an edge. For each subject, the correlation between any voxel pairs was calculated within the DPABI brain default template. Then, the whole-brain functional network was constructed by thresholding each correlation, and the DC value at the individual level was divided by the mean value of the whole brain DC, to obtain the standardized DC value. To improve the normality, Fisher’s R to Z and Z transformation were used to convert the individual related data into Fisher’s Z value. Finally, the Z-valued DC distribution map of each subject was obtained. The weighted DC graph was obtained by summation of the positive correlation coefficients (Z values) between a given voxel and all other voxels. The resulting images were smoothed with 6 mm full-width half maximum.

Statistical Methods

We analyzed demographic and clinical data using a one-way analysis of variance (ANOVA), χ2-test, and post hoc unpaired t tests. DC comparisons among the 3 groups were performed using a one-way analysis of covariance (ANCOVA) in DPABI with nuisance covariates, including age, sex, education level, and head movement parameters. These areas were extracted as a mask. A 2-sample post hoc t test within this mask was performed between each pair of the 3 groups (AlphaSim correction, P < 0.05). Partial correlation analysis (age, sex, education level, and head movement as covariables) was used to assess the correlation between DC values of brain regions with intergroup differences and HAMD and gastrointestinal symptoms rating scale (GSRS) scores. Finally, brain regions associated with HAMD scores were used as seed points (spherical regions with a radius of 6 mm centered on the peak coordinates of MNI in the DC map) for FC analysis based on seed points. Besides, to explore the relationship between seed-based FC abnormalities and depressive symptoms in DEP-IBS patients, the partial correlation analyses between FC abnormalities clusters and HAMD scores of the DEP-IBS patients were also performed. As we know, IBS has typical sex differences,[27] so we further use a 2 (sex) × 2 (group) mixed model ANOVA statistical method to analyzed the data results according to sex.

Results

Demographic Characteristics and Clinical Results

Demographic details of all subjects and clinical characteristics of patients are summarized in Table 1. There were no statistically significant differences among the 3 groups in sex, age, and education level. Both patient groups had higher HAMD scores (P < 0.001) than did the HC group. Moreover, the GSRS and HAMD scores were significantly higher in the DEP-IBS group than in the nDEP-IBS group (P < 0.001). The 2 patient groups did not differ significantly in IBS disease duration (P = 0.330). No significant sex specificity was found in the gastrointestinal symptoms and depressive symptoms of the 2 patient groups (DEP-IBS group: GSRS P = 0.673, HAMD P = 0.956; nDEP-IBS group: GSRS P = 0.521, HAMD P = 0.895).

Degree Centrality Analysis

One-way ANCOVA results (P < 0.05, AlphaSim multiple comparisons corrected) showed that the brain regions with differences in DC values among the 3 groups included the left insula (Table 2 and Fig. 1A), the left medial superior frontal gyrus (Table 2 and Fig. 1B), and the left precentral gyrus (Table 2 and Fig. 1C). Post hoc 2-sample t test comparison results (Fig. 1D-F) showed that compared with nDEP-IBS patients and HC, DEP-IBS patients showed reduced DC in the left insula (Fig. 1D), and increased the DC in the left precentral gyrus (Fig. 1F). There was no significant difference in the DC value of the above 2 brain regions between nDEP-IBS group and HC group. The DC value of left medial superior frontal gyrus in both groups was significantly lower than that in the HC group, while there was no statistically significant difference between the 2 groups of IBS patients (Fig. 1E). The ANOVA results showed that the DC values in the left insula and in the left medial superior frontal gyrus did not have a significant group-sex interaction effect (F = 0.924, P = 0.401; F = 1.203, P = 0.306). There was a significant group-sex interaction effect on the DC values in the left precentral gyrus (F = 4.800, P = 0.011). In the 2 patient groups, the DC values in the left precentral gyrus of female patients was significantly lower than that of male patients in each group (DEP-IBS group: P = 0.011; nDEP-IBS group: P = 0.040). There was no sex difference in the DC values of the left precentral gyrus in the HC group (P = 0.243).

Table 2

Brain Regions Showing Significant Differences in Degree Centrality Among 3 Groups

Brain area	R/L	MNI coordinates			Cluster size (voxels)	PeakF-values
		x	y	z
Insula	L	–33	18	0	33	11.37
Medial superior frontal gyrus	L	–6	51	12	36	9.23
Precentral gyrus	L	–54	12	–18	53	12.71

R, right; L, left; MIN, Montreal Neurological Institute.

One-way ANCOVA results (P < 0.05, AlphaSim multiple comparisons corrected) showed that the brain regions with differences in degree centrality values among the irritable bowel syndrome with depressive symptoms, irritable bowel syndrome without depressive symptoms, and healthy controls.

Figure 1

One-way analysis of covariance (ANCOVA) and post hoc 2-sample t test results of degree centrality among irritable bowel syndrome (IBS) with depressive symptoms (DEP-IBS) group, IBS without depressive symptoms (nDEP-IBS) group, and healthy controls (HC) group (P < 0.05, Alphasim correction). Brain regions with different degree centrality (DC) values among groups and sex distribution in each group: (A) left insula, (B) left medial superior frontal gyrus, and (C) left precentral gyrus. The distribution and comparison of DC values of brain regions (A-C) in the 3 groups (D-F). L.INS, left insula; L.medial SFG, left medial superior frontal gyrus; L.PreCG, left precentral gyrus.

Correlation Between Degree Centrality Value and Hamilton Depression Rating Scale Scores

Partial correlation analysis showed that the DC value of the left insula of DEP-IBS patients was significantly negatively correlated with HAMD score (r = –0.524, P = 0.009), as shown in Figure 2. Since no sex difference was found in the left insula DC values in the DEP-IBS group, the results of this correlation analysis were not further analyzed by sex.

Figure 2

Scatter plots depicting partial correlation between degree centrality (DC) and the clinical depression severity for irritable bowel syndrome with depressive symptoms patients. HAMD, Hamilton depression rating scale.

Secondary Seed-based Functional Connectivity Analysis

The left insula cluster (peak MNI: x = –33, y = 18, z = 0) was further used as a seed for the secondary seed-based FC analysis. Figure 3A shows the one-way ANCOVA results, and Figure 3B shows the post hoc 2-sample t test results. Compared with HC, patients with DEP-IBS showed increased left insula-related FC in the right inferior occipital gyrus, and decreased left insula-related FC in the right supplementary motor area (SMA) and the right postcentral gyrus. The above differences were not shown between the nDEP-IBS and HC groups. Compared with the nDEP-IBS group, patients with DEP-IBS showed the same pattern of FC changes as described above, and there were additional increased left insula-related FC in the left inferior parietal lobule and decreased left insula-related FC in the bilateral anterior cingulate cortex (ACC) (Fig. 3B). In addition, FC intensity between left insula and right SMA negatively correlated with the score of HAMD in DEP-IBS patients (r = –0.438, P = 0.037) (Fig. 3C).

Figure 3

One-way analysis of covariance (ANCOVA) (A) and post hoc 2-sample t test (B) results of left insula seed functional connectivity analyses differences among irritable bowel syndrome (IBS) with depressive symptoms (DEP-IBS) group, IBS without depressive symptoms (nDEP-IBS) group, and healthy controls (HC) group. AlphaSim multiple comparisons corrected is P < 0.05. Scatter plots depicting partial correlation between the left insula-related functional connectivity in the right supplementary motor area (SMA) and the Hamilton depression rating scale (HAMD) score for DEP-IBS patients (C). L, left; R, right; IPL, inferior parietal lobule; PostCG, postcentral gyrus; IOG, inferior occipital gyrus; ACC, anterior cingulate cortex; INS, insula.

Discussion

To the best of our knowledge, this study is the first to explore voxel-level whole-brain FC abnormalities in DEP-IBS patients using both DC and seed-based FC approaches. The main findings were as follows: (1) compared with nDEP-IBS patients and HC, DEP-IBS patients showed decreased DC in the left insula and increased DC in the left precentral gyrus. The DC values of the left medial superior frontal gyrus in both IBS groups were significantly lower than that in the HC group; (2) the DC abnormalities in the left insula correlated with depression severity in the DEP-IBS patients; and (3) changes of the left insula-related FC in the right SMA, the right postcentral gyrus, the right inferior occipital gyrus, the inferior parietal lobule and the bilateral ACC were further observed. In particular, the left insula-related FC abnormalities in the right SMA of DEP-IBS patients were associated with HAMD scores.

The insula, as part of the limbic system, is involved in the integration of visceral sensory functions and is responsible for sending pain signals to the frontal cortex, playing a crucial role in modulation of visceral sensation and movement, and pain-related emotional processing.[28] A large number of functional MRI (fMRI) studies have reported the functional and structural abnormalities of the insula in IBS patients, and there were significant correlations with variables such as gastrointestinal symptoms, anxiety and depression symptoms.[29-31] A positron emission tomography–computed tomography study based on fluorine-18-deoxyglucose found that after controlling gastrointestinal symptoms, the metabolic abnormalities of the insular and prefrontal cortex (PFC) of patients with functional dyspepsia remained, and there was a significant correlation with anxiety/depressive symptoms.[32] These results suggested that anxiety/depression may interfere with glucose metabolism in the insula and prefrontal lobes of FGID patients. A recent study found abnormal insular activity in patients with IBS in response to visceral stimulation and endocrine stress, which is associated with alexithymia.[33] Our study showed that DEP-IBS patients had decreased DC in the left insula and negatively correlated with the HAMD score. This negative correlation suggest that decreased DC was associated with the severity of pain-related negative emotions. A meta-analysis of task-state fMRI studies related to pain stimuli found that brain regions positively related to pain involve ACC, PFC, primary and secondary somatosensory cortex, insular cortex, and thalamus. The recognition of pain intensity activates the ventral pathway from the insular cortex to the bilateral PFC, while the recognition of pain space involves the dorsal pathway from the posterior parietal cortex to the right dorsolateral PFC. Both of the above tasks lead to the activation of ACC.[34] The ACC is closely connected with the SMA, premotor and primary motor cortex, which together form a module that generates and transmits the final decision about the nature of the sensory features of the discriminated pain stimuli. A secondary seed-based FC analysis in our study further revealed that changes in the left insula-related FC exist in patients with chronic visceral pain, including ACC, somatosensory cortex (postcentral gyrus) and posterior parietal cortex. In particular, abnormalities in the FC between insula and SMA were associated with depression severity in DEP-IBS patients. Growing evidence also suggests that SMA is not only a motor brain region, but also plays an important role in the integration of affective and cognitive functions.[35] Patients with depression show disrupted functional integration in SMA, which may be due to psychomotor retardation, and is considered one of the vital neurophysiological indicators of depression.[36] Numerous studies have found that functional changes in SMA networks can be used as predictors of antidepressant treatment outcomes, and may serve as a biomarker of treatment response.[37-39] Our previous study found that DEP-IBS patients had decreased gray matter volume in the insula and sensorimotor cortex.[19] Therefore, we speculate that the insula and SMA may be the key pathogenesis of DEP-IBS patients, and may provide a certain guiding value for the formulation of clinical treatment strategies.

The precentral gyrus is a sensorimotor region, consisting of the premotor cortex and the SMA, receiving information projected from the basal ganglia. Bhatt et al[40] found that the FC between the precentral gyrus and caudate nucleus was enhanced in IBS patients. The basal ganglia is involved in the integration of information with the cortical, thalamic regions and 3 specific pain processing regions, which is an important part of the “pain matrix.”[41] The caudate nucleus is the largest nucleus in the basal ganglia, which can produce analgesic effects when stimulated.[42] Our study showed that increased DC in the left precentral gyrus in DEP-IBS patients, which may be related to the depression that further strengthens the patients’ perception of visceral pain. More interestingly, our data showed that female patients with IBS had lower DC than male patients in the left precentral gyrus but not in the HC group, suggesting precentral’s ability in integrating and receiving visceral hypersensitivity information can be influenced by sex.

This study showed decreased DC in the left medial superior frontal gyrus in both IBS groups. MRI studies in chronic pain patients more consistently reported decreased gray matter density in the medial PFC (mPFC) and decreased functional connections between mPFC and other brain regions.[43,44] The activity of PFC is related to anti-hyperalgesia.[45] Studies suggest that mPFC can release endogenous opioids during tension pain.[46] Therefore, the decreased DC in mPFC may be related to the weakened function of the pain suppression system.

This study has several shortcomings. First, the sample size of IBS is relatively small, which may affect the credibility and universality of the results. Second, the study did not analyze other concomitant symptoms or factors, such as anxiety symptoms or traumatic experiences. Thirdly, the HAMD score of DEP-IBS patients was within the range of 10 to 17 points, which only belonged to the state of mild depression and did not meet the clinical diagnostic criteria for depression. In the future, the sample size should be further expanded to refine the concomitant symptoms and sex differences of IBS patients, and comparative analysis should be conducted with patients with simple depression for mutual verification.

The pathogenesis of depressive symptoms in IBS was associated with dysfunctions of fronto-limbic circuit and sensorimotor network, which were reflected by the abnormal FC pattern of these networks in our study. Especially, the insula and SMA, could be taken as a characteristic feature for depression in IBS. Identification of such alterations may facilitate earlier and more accurate diagnosis of depression in IBS, and development of effective treatment strategies.