ORIGINAL RESEARCH article

Front. Neurosci., 20 September 2022

Sec. Gut-Brain Axis

Volume 16 - 2022 | https://doi.org/10.3389/fnins.2022.985190

Altered mean apparent propagator-based microstructure and the corresponding functional connectivity of the parahippocampus and thalamus in Crohn’s disease

  • 1. Department of Radiology, Renji Hospital Affiliated to Shanghai Jiao Tong University School of Medicine, Shanghai, China

  • 2. Department of Gastroenterology, Renji Hospital Affiliated to Shanghai Jiao Tong University School of Medicine, Shanghai, China

  • 3. Department of Gastroenterology, Huadong Hospital, Shanghai Medical College, Fudan University, Shanghai, China

  • 4. Shanghai Key Laboratory of Magnetic Resonance, East China Normal University School of Physics and Electronics Science, Shanghai, China

Abstract

Crohn’s disease (CD) is a chronic and relapsing inflammatory bowel disorder that has been shown to generate neurological impairments, which has the potential to signify disease activity in an underlying neurological manner. The objective of this study was to investigate the abnormalities of brain microstructure and the corresponding functional connectivity (FC) in patients with CD, as well as their associations with disease condition. Twenty-two patients with CD and 22 age-, gender-, and education-matched healthy controls (HCs) were enrolled in this study. All subjects underwent mean apparent propagator (MAP)-MRI and resting-state functional magnetic resonance imaging (MRI) (rs-fMRI) data collection. Each patient was evaluated clinically for the condition and duration of the disease. The MAP metrics were extracted and compared between two groups. Pearson’s correlation analysis was conducted to determine the relationship between disease characteristics and significantly abnormal MAP metrics in the CD group. Regions of interest (ROIs) for ROI-wise FC analysis were selected based on their correlation with MAP metrics. Results showed that multiple brain regions, including the parahippocampus and thalamus, exhibited statistically significant differences in MAP metrics between CD patients and HCs. Additionally, CD patients exhibited decreased FC between the left parahippocampus and bilateral thalamus, as well as the right parahippocampus and bilateral thalamus. The findings of this work provide preliminary evidence that structural abnormalities in the parahippocampal gyrus (PHG) and thalamus, as well as decreased FC between them, may reflect the degree of inflammatory of the disease and serve as brain biomarkers for evaluating CD activity.

Introduction

Crohn’s disease (CD) is an inflammatory bowel disease (IBD) that is chronic and relapsing and may affect any portion of the digestive tract. In recent years, the disease incidence of CD in China has surged in tandem with the urbanization of Asian countries, population mobility, and dietary structure changes (; ). The disease condition typically requires lifetime medication and supportive care on account of relapsing digestive tract inflammation and recurring intestinal or extraintestinal symptoms, greatly contributing to the poor life quality of CD patients (). It has been proved that there are white matter (WM) lesions and neurological impairments existing in CD patients’ brains, which might represent an uncommon extraintestinal manifestation of the disease (). Given that no apparent brain structure abnormalities were found by conventional magnetic resonance imaging (MRI), more sophisticated imaging techniques and methodologies are required to study brain abnormalities underlying extraintestinal symptoms and to aid clinical diagnosis based on identified biomarkers (; ; ).

Although the classical tensor model in diffusion tensor imaging (DTI) appears to be a sensitive tool for studying WM changes, it comes with a significant modeling constraint. That is, this method is predicated on the presumption that water diffusion in the brain is unrestricted and follows Gaussian distribution (). However, as a matter of fact, water molecule diffusion often follows a non-Gaussian distribution in complex biological tissues owing to the presence of cell membranes, organelles, and liquid compartments (). These microstructures cannot be adequately described by the tensor model, and biases may be introduced in the corresponding diffusion metrics in the meantime. Mean apparent propagator (MAP)-MRI, a newly developed diffusion model, has the potential to circumvent this limitation. The MAP-MRI model makes no prior presuppositions regarding the behavior of water diffusion in tissues (; ). Instead, it is based on q-space sampling () and measures the probability density function of spin displacements in complex microstructures of brain tissue to examine the dispersion distribution of water molecules (). This method offers a more comprehensive, delicate, and accurate tissue characterization than the tensor model by quantifying the non-Gaussian nature of the diffusion process, reflecting physiologically significant microstructural characteristics with a higher degree of neuroanatomical specificity (; ).

It is salient to acknowledge that the vast majority of previous research has focused on either brain structure or brain function. To our knowledge, very few studies provide information on both brain structure and function in the CD sample (; ), mostly using separate analysis approaches. While beyond doubt useful, such approaches may overlook vital sources of intermodal information. Yet, such information is essential for directly detecting potential links and similarities between patterns of brain structure and function in a specific population. Using concurrent multimodal MRI techniques, such as advanced diffusion model establishment and resting-state functional MRI (rs-fMRI) analysis, it is possible to detect co-altered patterns of brain changes that may be partially missed by conventional separate analyses.

In the meantime, few researches have examined the relationship between brain structure/function regulation and the inflammatory activity level of the disease, mostly focusing instead on association with cognitive function, mood, or pain scale. The dysfunction of the brain-gut axis is believed to participate in the dysregulation of CD (; ; ). The brain-gut axis refers to a complex bi-directional interactive system comprised of multiple interconnections between the neuroendocrine pathways, the autonomous nervous system, and the gastrointestinal tract. Behavior and brain function may be affected by intestinal inflammation, and vice versa. Therefore, it is of great value to explore the connection between neural changes and inflammatory activity.

The present research used a novel approach of combining MAP-MRI data with regions of interest (ROI)-wise functional connectivity (FC) analysis to uncover brain alterations in CD patients and their relationships with disease activity conditions. We hypothesized that on account of the inseparability of brain structure and function, there were both MAP parameter and FC changes in specific brain regions between these two groups. We also hypothesized that the altered microstructural and functional features would be associated with the levels of inflammatory activity of CD and, hence, may serve as sensitive clinical biomarkers for evaluating disease activity in CD patients.

Materials and methods

Participants

This study was approved by the Research Ethics Committee of Renji Hospital Affiliated to Shanghai Jiao Tong University School of Medicine. Written informed consent was obtained from each subject. A total of fifty patients with CD were recruited from April 2020 to January 2022. Inclusion criteria were as follows: (1) right-handed; (2) 18–55 years old; and (3) education level greater than 6 years. Exclusion criteria: (1) intestinal-related abdominal surgery; (2) recent use of glucocorticoids, biologics, or psychotropic medications; (3) organic brain lesions or obvious WM degeneration; and (4) MRI scan-related contraindications such as claustrophobia or metallic implants. Twenty-two patients were finally enrolled in this study after being evaluated and screened by a gastroenterologist with extensive expertise. The medical records of the patients were reviewed throughout the evaluation to acquire endoscopic and hematological data, as well as supplementary information such as disease duration in months.

Twenty-two healthy control subjects (HCs) were recruited via advertisements and were matched for age, gender, education level, and handedness. They had no digestive or pain-related disorders, and their colonoscopy examinations were negative. The inclusion and exclusion criteria for the HC group were identical to those for the patient group.

Clinical examinations

The serum C-reactive protein (CRP), erythrocyte sedimentation rate (ESR), fecal calprotectin (Calpro), and the Crohn’s Disease Activity Index (CDAI) () were collected to assess the disease condition of each patient, including the degree of inflammatory activity and severity. The CDAI is a globally accepted index for accurately evaluating the disease severity and tentatively predicting the therapeutic outcome of CD patients. CDAI ≤ 150 was considered remission, and CDAI > 150 was considered an active disease. In addition, the disease duration was recorded in months.

Magnetic resonance imaging data acquisition

All participants were scanned using a 3.0 Tesla MR system (SIEMENS MAGNETOM Prisma) equipped with a 64-channel phase-array head coil. Participants were advised to stay awake throughout the scanning, with their eyes closed and ears muffled, and to avoid thinking about anything in particular. The three-dimensional T1-weighted anatomical images were obtained in the sagittal orientation using the following parameters: TR = 1, 800 ms, TE = 2.28 ms, slice thickness = 1 mm, flip angle = 8°, field of view = 256 × 256 mm2, matrix = 256 × 256, and number of slices = 160. The rs-fMRI data were collected using an echo planar imaging sequence with the following parameters (multi-band, acceleration factor = 2): TR = 2,000 ms, TE = 30 ms, slice thickness = 2 mm, flip angle = 90°, field of view = 230 × 230 mm2, matrix = 64 × 64, number of slices = 70, and total volume = 220 which was acquired in 7 min and 33 s. The diffusion-weighted images were obtained with the following parameters (multi-band, acceleration factor = 2): TR = 6,800 ms, TE = 75 ms, field of view = 192 × 192 mm2, matrix = 128 × 128, number of slices = 92, slice thickness = 1.5 mm, gap = 0. Diffusion-sensitizing gradients were applied along 100 non-collinear directions with a b-value of 3, 000 s/mm2, and a reference image was acquired at b = 0 s/mm2.

Magnetic resonance imaging data preprocessing

Diffusion-weighted images underwent eddy current and motion correction using the DiffusionKit eddy tool (). MAP-MRI parameters were calculated using NeuDiLab, a software based on an open-source tool named DIPY (Diffusion Imaging in Python) (). In a nutshell, the MAP was fit with a radial order of six. A positivity constraint and Laplacian regularization with a weighting of 0.05 were utilized for denoizing during fitting. The ultimate generating parameters were recommended for robust fitting by the developers of the DIPY MAP-MRI toolbox. Finally, maps of MAP-MRI parameters were constructed, encompassing mean square displacement (MSD), non-Gaussianity (NG), non-Gaussianity axial (NGAx), non-Gaussianity vertical (NGRad), Q-space inverse variance (QIV), return to the origin probability (RTOP), return to the axis probability (RTAP), and return to the plane probability (RTPP). They were all normalized to the standard Montreal Neurological Institute (MNI) space in two steps using MRIcron1 and then smoothed with a 6 mm full width at half-maximum (FWHM) Gaussian kernel with a reslicing resolution of 2 × 2 × 2 mm3.

The rs-fMRI data were preprocessed in MATLAB (MathWorks, Natick, Massachusetts, United States) using the Data Processing and Analysis for Brain Imaging (DPABI) toolkit package (version v4.3).2 The main procedures were brain extraction, slice timing correction, rigid-body motion correction, spatial smoothing using a 6 mm FWHM Gaussian kernel, and 150 s high-pass temporal filtering. Moreover, 24 head-motion parameters (6 head motion parameters at baseline, 6 head motion parameters at 1 time point prior, and the 12 associated squared items) were utilized as nuisance variables for regression, along with cerebrospinal fluid and WM signal.

Group comparison and correlation analysis of mean apparent propagator parameters

A two-sample independent t-test was used to compare all of the standardized MAP parameter maps between CD patients and HCs, using the Statistical Parametric Mapping 12 (SPM12)3 toolbox in MATLAB. Given the relatively small sample size, a permutation test with a threshold free cluster enhancement (number of permutations = 1,000, p < 0.01) was used to extract the brain areas with significant inter-group differences. The automated anatomical labeling (AAL) atlas () was used to identify the above-mentioned brain regions. The mean values of MAP parameters, including MSD, NG, NGAx, NGRad, QIV, RTAP, RTOP, and RTPP in the CD group, were calculated from the regions exhibiting significant group differences. After adjusting for age, gender, education level, and body mass index (BMI), Pearson’s partial correlation analyses were performed to assess the relationship between clinical characteristics and MAP parameters using SPSS v25 (IBM, Armonk, NY, United States). All p-values < 0.05 were considered statistically significant, after false discovery rate (FDR) correction for multiple comparisons.

Regions of interest-wise functional connectivity analysis

ROIs were defined as brain regions where MAP values were significantly correlated with clinical characteristics. We conducted an ROI-to-ROI FC analysis using the DPABI toolkit. The mean time series for each ROI were extracted using the preprocessed rs-fMRI data. The pairwise FC strength was then computed using Pearson’s correlation coefficients on the time series and Fisher’s r-to-z transformation to convert the correlation coefficients into z-scores. CD-related changes in FC were evaluated by the group comparison between patients and HCs. All p-values < 0.05 were considered statistically significant after FDR correction for multiple comparisons.

Results

Demographic and clinical characteristics

Table 1 shows the demographic and clinical characteristics of the two groups. No significant differences were found in age, gender, or level of education between patients and HCs (two-sample t-tests and a chi-squared test, all p-values > 0.1).

TABLE 1

Patients (n = 22)HCs (n = 22)P-value
Age, years33.60 ± 14.6737.55 ± 7.430.292
Gender (male/female)12/1012/101.00
Education level, years12.06 ± 3.2312.80 ± 1.320.487
BMI19.15 ± 2.93/
Disease duration, months59.78 ± 72.67/
CRP, mg/dL24.52 ± 24.41/
ESR, mm/h38.00 ± 27.70/
Calpro, mg/L323.90 ± 151.18/
CDAI171.12 ± 68.58/

Demographic and clinical characteristics of all participants.

Values are represented as the mean ± standard deviation, except for the gender distribution. Two-sample t-tests were performed to assess group differences for age and education, and a chi-squared test for gender. P-value < 0.05 was considered to be statistically significant.

HC, healthy control; BMI, body mass index; CRP, C-reactive protein; ESR, erythrocyte sedimentation rate; Calpro, calprotectin; CDAI, Crohn’s Disease active index.

Group comparison of the mean apparent propagator parameters

The brain regions with significantly different MAP parameters are shown in Table 2. They are as follows: bilateral parahippocampal gyrus (PHG), bilateral thalamus (THA), bilateral insula (INS), left hippocampus (HIP.L), left putamen (PUT. L), left amygdala (AMYG.L), left temporal pole: superior temporal gyrus (TPOsup.L), left Rolandic operculum (ROL.L), left fusiform gyrus (FFG.L), right middle frontal gyrus (MFG.R), right medial superior frontal gyrus (SFGmed.R), and right anterior cingulate and paracingulate gyri (ACG.R).

TABLE 2

MAP parametersBrain regionsBAMNI
t-valueVoxels
XYZ
MAP-MSDPHG.L28−14−26−524.43283
PHG.R3526−14−305.29119
HIP.L−30−10−164.4683
PUT.L34−282−144.7397
MFG.R104458124.37188
TPOsup.L22−581004.5171
ROL.L13−44−425.651,140
SFGmed.R9456384.01169
MAP-NGTHA.R8−1812−4.13109
ROL.L13−38−220−4.58178
MAP-NGAxMFG.R46485484.1287
THA.R8−1812−4.20126
ROL.L13−38−220−4.77218
MAP-NGRadMFG.R10485484.0577
THA.R16−1814−4.0975
INS.L13−38−218−4.0974
MAP-QIVPHG.L35−14−18−484.37395
PHG.R3526−12−324.5075
HIP.L35−30−10−164.3270
PUT.L34−282−125.2389
MFG.R104458124.3478
ROL.L13−36−2165.54884
MAP-RTAPAMYG.L34−264−16−5.3591
FFG.L36−22−40−12−5.5182
ROL.L13−42−62−5.50994
MAP-RTOPAMYG.L34−264−16−5.08122
FFG.L36−22−40−12−5.1286
ROL.L13−42−42−5.441,009
THA.L−14−242−3.8065
THA.R10−186−4.4291
MAP-RTPPMFG.R10485484.1790
PHG.R3526−16−28−4.8298
PHG.L28−22−16−26−4.68130
THA.R2710−26−8−4.51449
HIP.R1340−22−6−4.32973
INS.L13−44−42−5.78835
INS.R1332−2026−3.9682
ACG.R2462420−4.9877

Significant differences in the MAP parameters between CD patients and HCs.

The statistical threshold was set at p < 0.01.

BA, Brodmann area; CD, Crohn’s Disease; HC, healthy controls; PHG, parahippocampal gyrus; THA, thalamus; INS, insula; HIP.L, left hippocampus; PUT. L, left putamen; AMYG.L, left amygdala; TPOsup.L, left temporal pole: superior temporal gyrus; ROL.L, left Rolandic operculum; FFG.L, left fusiform gyrus; MFG.R, right middle frontal gyrus; SFGmed.R, right medial superior frontal gyrus; ACG.R, right anterior cingulate and paracingulate gyri.

Associations between clinical characteristics and mean apparent propagator parameters

Table 3 shows the brain regions where MAP parameters were significantly correlated with clinical characteristics in patients. The MSD value in the left parahippocampal gyrus (PHG.L) was negatively correlated with CDAI, while the NG value in the right thalamus (THA.R) and the QIV value in the PHG.L were both positively correlated with CDAI (Figure 1). The values of NG, NGAx, RTOP, and RTAP in the ROL.L, and the NGRad value in the left insula (INS.L) were both positively correlated with disease duration. The RTOP value in the left thalamus (THA.L) was positively correlated with CRP. The RTPP value in the right parahippocampal gyrus (PHG.R) was negatively correlated with both ESR and Calpro. (All p-values < 0.05, after FDR correction).

TABLE 3

Disease durationCRPESRCalproCDAI
MAP-MSD in PHG.Lr−0.03−0.33−0.210.01−0.61
P0.930.260.470.990.02*
MAP-NG in ROL.Lr0.770.150.240.090.36
p0.001**0.610.400.770.21
MAP-NG in THA.Rr0.080.360.300.050.53
p0.760.210.290.870.019*
MAP-NGAx in ROL.Lr0.760.090.190.080.36
p0.001**0.770.520.780.21
MAP-NGRad in INS.Lr0.740.340.320.170.32
p0.001**0.230.270.560.27
MAP-QIV in PHG.Lr−0.030.300.250.030.70
p0.910.290.390.920.005**
MAP-RTAP in ROL.Lr0.62−0.14−0.060.210.06
p0.01*0.630.850.460.84
MAP-RTOP in ROL.Lr0.63−0.040.050.180.04
p0.009**0.880.860.540.90
MAP-RTOP in THA.Lr−0.100.600.330.100.02
p0.710.023*0.250.750.94
MAP-RTPP in PHG.Rr0.190.250.56−0.570.23
p0.470.390.037*0.032*0.44

Correlations between clinical characteristics and MAP parameters.

The statistical threshold was set at p < 0.05.

*p < 0.05, **p < 0.01.

PHG.L, left parahippocampal gyrus; THA.R, right thalamus; THA.L, left thalamus; INS.L, left insula; PHG.R, right parahippocampal gyrus; ROL.L, left Rolandic operculum.

FIGURE 1

Functional connectivity analysis

Four ROI to ROI FC were significantly lower in CD patients compared to HCs. They were THA.L-PHG.L (t = −3.117, p = 0.034), THA.L-PHG.R (t = −3.407, p = 0.021), THA.R-PHG.R (t = −2.959, p = 0.029), and THA.R-PHG.L (t = −4.485, p = 0.006) FC (Figure 2).

FIGURE 2

Discussion

In summary, our study found that the majority of MAP parameters in the parahippocampus and thalamus were significantly different in CD patients from those in the HCs, and that the corresponding FC between the parahippocampus and thalamus were significantly decreased in CD patients.

The MAP-MRI diffusion model, which was recently developed, provides several novel, measurable parameters for capturing previously concealed inherent properties of nervous tissue microstructure. The MSD can be obtained using diffusion propagator models to determine the average amount of diffusion. Zero displacement probabilities, involving the RTOP, RTAP, and RTPP, can be calculated to quantify various features of the three-dimensional diffusion process (; ; ; ). The RTOP is a metric estimated from three-dimensional q-space data, which reflects the probability of water molecules experiencing zero net displacements between the applications of two diffusion sensitizing gradients (; ). The RTAP and RTPP are two-dimensional and one-dimensional variants of RTOP, respectively. RTAP is decomposed from RTOP perpendicular to the direction of the primary eigenvector and reflects the presence of restrictive barriers in the radial orientation, while RTPP is decomposed parallel to the direction of the primary eigenvector and demonstrates the presence of restrictive obstacles in the axial direction (; ). Previous studies have shown that MAP indices are clinically applicable in a variety of neuroradiological subspecialties (; ; ), and MAP-MRI outperformed the conventional DTI in diagnosing and evaluating Parkinson’s disease (). The present study discovered significant differences in MAP parameters between CD patients and HCs in specific brain regions, suggesting the potential value of MAP-MRI in the clinical diagnosis of CD.

Previous studies employing structural and functional MRI have reported some substantial abnormalities in the brain structure and activity of CD patients. showed that CD patients exhibited lower gray matter volume in the frontal lobe and cingulate cortex. , however, found no significant structural difference between CD patients and HCs using voxel-based morphometry analysis, but rather increased FC between the frontoparietal network and the salience network. Prior task-based and rs-fMRI studies have discovered abnormal functional alterations in specific brain areas. In a task-related study, patients with CD had decreased activity in the medial temporal lobe, insula, putamen, and cerebellum (). In addition, another study using the stress-evoking task found that CD patients had increased activation in the midcingulate cortex (; ). A study using intra-network FC analysis with rs-fMRI data showed significantly increased connectivity within the executive control network and default mode network (). According to , the intrinsic FC between the amygdala and insula, parahippocampus, and anterior middle cingulate cortex significantly decreased and was correlated with the disease duration.

In our study, patients with CD displayed microstructural abnormalities in the PHG and thalamus, as well as decreased FC between them. The PHG and thalamus are two knots of the limbic system that have long been well-established as critical components in affect and emotion processing. Some evidence indicates that CD patients often suffer from mood disorders, overreaction to stress, deficiency in concentration, and gut flora imbalance (), precipitating mental symptoms such as anxiety and depression. Holland and Goadsby proposed that the limbic system structures had an effect on the pain matrix and were particularly involved in nociception inhibition (). The PHG is involved in the familiarity and retrieval process of memory. The thalamus acts as a relay station, receiving sensory input and relaying it to the cortex. The decreased FC between these two regions might be a distinctive feature of CD in pain or pain context recollection. Along with the decreased FC between the PHG and thalamus, demonstrated that CD patients with significantly higher self-reported anxiety or depression had increased FC between the amygdala, thalamus, and orbitofrontal cortex. Both previous research and our findings indicate that abnormal FC within the limbic system could contribute to emotional symptoms in CD patients ().

CD, as a chronic IBD, primarily affects the gastrointestinal tract. Some recent studies have proposed that the microbiota-gut-brain axis may be a critical contributor in understanding the concomitant mental health issues in IBD and related disorders (; ; ; ). Animal-model research showed that chronic gastrointestinal inflammation might induce anxiety-like behavior and alter biochemistry of the central nervous system in mice (). Furthermore, gastrointestinal inflammation is associated with hippocampal changes in inflammatory gene expression and glutamatergic transmission, resulting in increased brain excitability and cognitive impairments (). Serotonin is suspected to play a vital role in gut-to-brain transmission. An immunohistochemical study noted that serotonergic fibers connecting the thalamus to the cortex may exert a considerable influence on affective function (). All of these factors might account for our observations of altered microstructure in the hippocampus and thalamus of CD patients.

There were several limitations to our study. First, the sample size is too small, and it would be preferable to incorporate more clinical information, including manifestations and treatments, and further studies can be conducted to explore the relationship between the clinical manifestations and brain function changes, such as whether the results are consistent between active and remission stages of the disease, or whether they are reversible between courses. Second, now that the thalamus and PHG are pivotal parts of the classical hippocampal circuit within the limbic system, further studies using multimodal MRI can be conducted to investigate the structural and functional network (e.g., small world network) alterations related to emotion processing. Third, since a lot of new brain function templates have admittedly emerged in recent years, additional research is necessary to use other brain templates or atlases including but not limited to the High Angular Resolution Diffusion Imaging (HARDI) template () based on the International Consortium for Brain Mapping (ICBM)-152 space, to concentrate on more complex neuronal microarchitecture, which may yield different and more intriguing findings.

In summary, the altered microstructure and FC of the brain reveal abnormal brain activity in CD patients, which may be related to the pathophysiological mechanism of the disease through the brain-gut axis. The impairment of parahippocampus gyrus and thalamus may serve as a neurological biomarker for disease activity.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving human participants were reviewed and approved by the Research Ethics Committee of Renji Hospital Affiliated to Shanghai Jiao Tong University School of Medicine. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author contributions

YQ and QL contributed to the conceptualization, formal analysis, investigation, methodology, and writing – original draft. DW contributed to the review and editing of the manuscript. YiZ performed the formal analysis. JC performed data curation and investigation. ZC and YaZ contributed to the supervision of the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This research was supported by the Shanghai Science and Technology Committee Project (Natural Science Funding; 20ZR1433200 and 17411964300), the Shanghai Science and Technology Commission Research Project (14441903103), the Shanghai Municipal Key Clinic Specialty, the National Natural Science Foundation of China (82171885), and the Shanghai Science and Technology Committee Project (the Explorer Project Funding; 21TS1400700).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

  • 1

    AgostiniA.BallottaD.RighiS.MorettiM.BertaniA.ScarcelliA.et al (2017). Stress and brain functional changes in patients with Crohn’s disease: A functional magnetic resonance imaging study.Neurogastroenterol. Motil.29110. 10.1111/nmo.13108

  • 2

    AgostiniA.BenuzziF.FilippiniN.BertaniA.ScarcelliA.FarinelliV.et al (2013a). New insights into the brain involvement in patients with Crohn’s disease: A voxel-based morphometry study.Neurogastroenterol. Motil.25147e82. 10.1111/nmo.12017

  • 3

    AgostiniA.FilippiniN.BenuzziF.BertaniA.ScarcelliA.LeoniC.et al (2013b). Functional magnetic resonance imaging study reveals differences in the habituation to psychological stress in patients with Crohn’s disease versus healthy controls.J. Behav. Med.36477487. 10.1007/s10865-012-9441-1

  • 4

    AvramA.SarllsJ.BarnettA.ÖzarslanE.ThomasC.IrfanogluM.et al (2016). Clinical feasibility of using mean apparent propagator (MAP) MRI to characterize brain tissue microstructure.Neuroimage127422434. 10.1016/j.neuroimage.2015.11.027

  • 5

    BaoC.LiuP.LiuH.JinX.ShiY.WuL.et al (2018). Difference in regional neural fluctuations and functional connectivity in Crohn’s disease: A resting-state functional MRI study.Brain Imaging Behav.1217951803. 10.1007/s11682-018-9850-z

  • 6

    BestW.BecktelJ.SingletonJ. (1979). Rederived values of the eight coefficients of the Crohn’s disease activity index (CDAI).Gastroenterology77843846.

  • 7

    BonazB.BernsteinC. (2013). Brain-gut interactions in inflammatory bowel disease.Gastroenterology1443649. 10.1053/j.gastro.2012.10.003

  • 8

    Boscolo GalazzoI.BrusiniL.ObertinoS.ZucchelliM.GranzieraC.MenegazG. (2018). On the viability of diffusion MRI-based microstructural biomarkers in ischemic stroke.Front. Neurosci.12:92. 10.3389/fnins.2018.00092

  • 9

    ChenQ.WangY.QiuY.WuX.ZhouY.ZhaiG. (2020). A deep learning-based model for classification of different subtypes of subcortical vascular cognitive impairment with FLAIR.Front. Neurosci.14:557. 10.3389/fnins.2020.00557

  • 10

    FanY.BaoC.WeiY.WuJ.ZhaoY.ZengX.et al (2019). Altered functional connectivity of the amygdala in Crohn’s disease.Brain Imaging Behav.14, 20972106. 10.1007/s11682-019-00159-8

  • 11

    FanY.BaoC.WeiY.WuJ.ZhaoY.ZengX.et al (2020). Altered functional connectivity of the amygdala in Crohn’s disease.Brain imaging Behav.1420972106. 10.1007/s11682-019-00159-8

  • 12

    FickR.WassermannD.CaruyerE.DericheR. (2016). MAPL: Tissue microstructure estimation using laplacian-regularized MAP-MRI and its application to HCP data.Neuroimage134365385. 10.1016/j.neuroimage.2016.03.046

  • 13

    GeisslerA.AndusT.RothM.KullmannF.CaesarI.HeldP.et al (1995). Focal white-matter lesions in brain of patients with inflammatory bowel disease.Lancet345897898. 10.1016/s0140-6736(95)90013-6

  • 14

    GracieD.GuthrieE.HamlinP.FordA. (2018). Bi-directionality of brain-gut interactions in patients with inflammatory bowel disease.Gastroenterology15416351646.e3. 10.1053/j.gastro.2018.01.027

  • 15

    GracieD.HamlinP.FordA. (2019). The influence of the brain-gut axis in inflammatory bowel disease and possible implications for treatment.Lancet Gastroenterol. Hepatol.4632642. 10.1016/s2468-1253(19)30089-5

  • 16

    HollanderD. (2003). Inflammatory bowel diseases and brain-gut axis.J. Physiol. Pharmacol.54 (Suppl 4)183190.

  • 17

    HosseinborA.ChungM.WuY.AlexanderA. (2013). Bessel fourier orientation reconstruction (BFOR): An analytical diffusion propagator reconstruction for hybrid diffusion imaging and computation of q-space indices.Neuroimage64650670. 10.1016/j.neuroimage.2012.08.072

  • 18

    HouJ.MohantyR.NairVaDoddK.Beniwal-PatelP.SahaS.et al (2019). Alterations in resting-state functional connectivity in patients with Crohn’s disease in remission.Sci. Rep.9:7412. 10.1038/s41598-019-43878-0

  • 19

    KhairnarA.Ruda-KucerovaJ.SzabóN.DrazanovaE.ArabA.Hutter-PaierB.et al (2017). Early and progressive microstructural brain changes in mice overexpressing human α-Synuclein detected by diffusion kurtosis imaging.Brain Behav. Immun.61197208. 10.1016/j.bbi.2016.11.027

  • 20

    KornelsenJ.WilsonA.LabusJ. S.WitgesK.MayerE. A.BernsteinC. N. (2020). Brain resting-state network alterations associated with Crohn’s disease.Front. Neurol.11:48. 10.3389/fneur.2020.00048

  • 21

    LeH.ZengW.ZhangH.LiJ.WuX.XieM.et al (2020). Mean apparent propagator MRI Is better than conventional diffusion tensor imaging for the evaluation of Parkinson’s disease: A prospective pilot study.Front. Aging Neurosci.12:563595. 10.3389/fnagi.2020.563595

  • 22

    LeeM.KimY.JahngG.KwonO. (2018). Angular resolution enhancement technique for diffusion-weighted imaging (DWI) using predicted diffusion gradient directions.Neuroimage183836846. 10.1016/j.neuroimage.2018.08.072

  • 23

    LvK.FanY.XuL.XuM. (2017). Brain changes detected by functional magnetic resonance imaging and spectroscopy in patients with Crohn’s disease.World J. Gastroenterol.2336073614. 10.3748/wjg.v23.i20.3607

  • 24

    LvK.SongW.TangR.PanZ.ZhangY.XuY.et al (2018). Neurotransmitter alterations in the anterior cingulate cortex in Crohn’s disease patients with abdominal pain: A preliminary MR spectroscopy study.Neuroimage. Clin.20793799. 10.1016/j.nicl.2018.09.008

  • 25

    MayerE.TillischK. (2011). The brain-gut axis in abdominal pain syndromes.Annu. Rev. Med.62381396. 10.1146/annurev-med-012309-103958

  • 26

    MogilevskiT.BurgellR.AzizQ.GibsonP. (2019). Review article: The role of the autonomic nervous system in the pathogenesis and therapy of IBD.Aliment. Pharmacol. Ther.50720737. 10.1111/apt.15433

  • 27

    NgS.ShiH.HamidiN.UnderwoodF.TangW.BenchimolE.et al (2017). Worldwide incidence and prevalence of inflammatory bowel disease in the 21st century: A systematic review of population-based studies.Lancet39027692778. 10.1016/s0140-6736(17)32448-0

  • 28

    NingL.WestinC.RathiY. (2015). Estimating diffusion propagator and its moments using directional radial basis functions.IEEE Trans. Med. Imaging3420582078. 10.1109/tmi.2015.2418674

  • 29

    OlsonD.ArpinarV.MuftulerL. (2019). Optimization of q-space sampling for mean apparent propagator MRI metrics using a genetic algorithm.Neuroimage199237244. 10.1016/j.neuroimage.2019.05.078

  • 30

    ÖzarslanE.KoayC.ShepherdT.KomloshM.İrfanoğluM.PierpaoliC.et al (2013). Mean apparent propagator (MAP) MRI: A novel diffusion imaging method for mapping tissue microstructure.Neuroimage781632. 10.1016/j.neuroimage.2013.04.016

  • 31

    ThomannA.ReindlW.WüstenbergT.KmucheD.EbertM.SzaboK.et al (2019). Aberrant brain structural large-scale connectome in Crohn’s disease.Neurogastroenterol. Motil.31:e13593. 10.1111/nmo.13593

  • 32

    ThomannA. K.SchmitgenM. M.KmucheD.EbertM. P.ThomannP. A.SzaboK.et al (2021). Exploring joint patterns of brain structure and function in inflammatory bowel diseases using multimodal data fusion. Neurogastroenterol. Motil.33:e14078. 10.1111/nmo.14078

  • 33

    TorresJ.MehandruS.ColombelJ.Peyrin-BirouletL. (2017). Crohn’s disease.Lancet38917411755. 10.1016/s0140-6736(16)31711-1

  • 34

    Tzourio-MazoyerN.LandeauB.PapathanassiouD.CrivelloF.EtardO.DelcroixN.et al (2002). Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain.Neuroimage15273289. 10.1006/nimg.2001.0978

  • 35

    VarentsovaA.ZhangS.ArfanakisK. (2014). Development of a high angular resolution diffusion imaging human brain template.Neuroimage91177186.

  • 36

    WangP.WengL.XieS.HeJ.MaX.LiB.et al (2021). Primary application of mean apparent propagator-MRI diffusion model in the grading of diffuse glioma.Eur. J. Radiol.138:109622. 10.1016/j.ejrad.2021.109622

  • 37

    XieS.ChenL.ZuoN.JiangT. (2016). DiffusionKit: A light one-stop solution for diffusion MRI data analysis.J. Neurosci. Methods273107119. 10.1016/j.jneumeth.2016.08.011

  • 38

    ZhangS.ChenF.WuJ.LiuC.YangG.PiaoR.et al (2022). Altered structural covariance and functional connectivity of the insula in patients with Crohn’s disease.Quant. Imaging Med. Surg.1210201036. 10.21037/qims-21-509

  • 39

    YeungA. W. K. (2021). Structural and functional changes in the brain of patients with Crohn’s disease: an activation likelihood estimation meta-analysis. Brain Imaging Behav. 15, 807818. 10.1007/s11682-020-00291-w

  • 40

    ZucchelliM.BrusiniL.Andrés MéndezC.DaducciA.GranzieraC.MenegazG. (2016). What lies beneath? Diffusion EAP-based study of brain tissue microstructure.Med. Image Anal.32145156. 10.1016/j.media.2016.03.008

Summary

Keywords

mean apparent propagator, diffusion-weighted imaging, resting-state magnetic resonance imaging, functional connectivity, Crohn’s disease

Citation

Qiu Y, Li Q, Wu D, Zhang Y, Cheng J, Cao Z and Zhou Y (2022) Altered mean apparent propagator-based microstructure and the corresponding functional connectivity of the parahippocampus and thalamus in Crohn’s disease. Front. Neurosci. 16:985190. doi: 10.3389/fnins.2022.985190

Received

03 July 2022

Accepted

18 August 2022

Published

20 September 2022

Volume

16 - 2022

Edited by

Shengxiang Liang, Fujian University of Traditional Chinese Medicine, China

Reviewed by

Sheng Hu, Anhui University of Chinese Medicine, China; Yangding Li, Hunan Normal University, China

Updates

Copyright

*Correspondence: Zhijun Cao, Yan Zhou,

These authors share first authorship

This article was submitted to Gut-Brain Axis, a section of the journal Frontiers in Neuroscience

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics