Structural and Functional Brain Connectivity of People with Obesity and Prediction of Body Mass Index Using Connectivity (original) (raw)

Open Access

Peer-reviewed

Research Article

Structural and Functional Brain Connectivity of People with Obesity and Prediction of Body Mass Index Using Connectivity

PLOS

x

Figures

Abstract

Obesity is a medical condition affecting billions of people. Various neuroimaging methods including magnetic resonance imaging (MRI) have been used to obtain information about obesity. We adopted a multi-modal approach combining diffusion tensor imaging (DTI) and resting state functional MRI (rs-fMRI) to incorporate complementary information and thus better investigate the brains of non-healthy weight subjects. The objective of this study was to explore multi-modal neuroimaging and use it to predict a practical clinical score, body mass index (BMI). Connectivity analysis was applied to DTI and rs-fMRI. Significant regions and associated imaging features were identified based on group-wise differences between healthy weight and non-healthy weight subjects. Six DTI-driven connections and 10 rs-fMRI-driven connectivities were identified. DTI-driven connections better reflected group-wise differences than did rs-fMRI-driven connectivity. We predicted BMI values using multi-modal imaging features in a partial least-square regression framework (percent error 15.0%). Our study identified brain regions and imaging features that can adequately explain BMI. We identified potentially good imaging biomarker candidates for obesity-related diseases.

Citation: Park B-y, Seo J, Yi J, Park H (2015) Structural and Functional Brain Connectivity of People with Obesity and Prediction of Body Mass Index Using Connectivity. PLoS ONE 10(11): e0141376. https://doi.org/10.1371/journal.pone.0141376

Editor: Martin Walter, Leibniz Institute for Neurobiology, GERMANY

Received: June 16, 2015; Accepted: October 6, 2015; Published: November 4, 2015

Copyright: © 2015 Park et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: This study (H.P.) was supported by the Basic Science Research Program through the National Research Foundation of Korea (grant number 20100023233, http://www.nrf.re.kr/nrf_eng_cms/). This study (H.P.) was also supported in part by the Institute for Basic Science (grant number IBS-R015-D1, http://www.ibs.re.kr/). Data were provided by the Human Connectome Project, WU-Minn Consortium (NIH 1U54MH091657, http://www.humanconnectome.org/). There were no other sources of funding than those listed above. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Obesity is a medical condition affecting billions of people worldwide and is considered a serious health issue in the 21st century [1]. Obesity is a key factor for many cardiovascular diseases such as diabetes, hypertension, and stroke [2]. People with obesity have a high body mass index (BMI), a measurement based on weight and height [3]. While BMI is not the most accurate measure of accumulated body fat, it is a convenient and widely adopted measure.

Many studies have adopted neuroimaging tools, including various flavors of magnetic resonance imaging (MRI), single photon emission computed tomography (SPECT), and positron emission tomography (PET), to extract information related to obesity [2,4–6]. MRI is especially useful for assessing both structural and functional brain information. Neuroimaging studies have revealed that specific brain systems are associated with obesity [7–9]. Altered brain activity in brain systems related to reward, motor, cognition, control, and attention has been observed [7]. The reward system has been repeatedly identified as crucial for obesity, and the hypothalamus was shown to play an important role [10–15]. Hormones regulating eating behavior such as leptin and ghrelin convey appetite-related information to the hypothalamus and suppress or promote a person’s food intake [16,17]. Alterations in levels of these hormones confuse the brain’s reward system and might promote weight gain [16,17]. Brain regions controlling the reward system are functionally or structurally related to the hypothalamus and affect food regulation [18].

Multi-modal images of a given brain can be obtained due to recent advances in neuroimaging. Different imaging modalities might provide complementary information and thus lead to better overall information to assess disease-specific changes. We chose two different MRI methods, resting state functional MRI (rs-fMRI) and diffusion tensor imaging (DTI), to explore the brains of non-healthy weight (HW) subjects. Raw data from neuroimaging need to be processed using computer algorithms to extract relevant information [19–23].

One type of analysis algorithm, known as connectivity analysis, focuses on how activities in one region correlate with activities in another region. Connectivity analysis allows observation of the whole brain as a complex interconnected network [24–26]. Resting-state fMRI measures local brain activity using the blood oxygen level-dependent effect (BOLD), while DTI provides in vivo neuronal fiber information using anisotropic water diffusion. DTI has several limitations and cannot distinguish between efferent and afferent connections; however, it is the only practical option for assessing in vivo fiber information. Connectivity derived from fMRI, referred to as functional connectivity, reflects functional correlations between time series data of different regions. Connectivity derived from DTI is known as structural connectivity, because DTI reflects actual structural connections via neuronal fibers between regions.

Previous studies have explored functional brain connectivity of non-HW subjects and have reported decreased connectivity in the right anterior cingulate cortex and left insula, and increased connectivity in the bilateral precuneus, putamen, and posterior cingulate cortex [27,28]. Increased BOLD signals were observed in insula and orbitofrontal cortex [10]. Another study explored structural brain connectivity of non-HW subjects and reported low structural integrity of the connection between frontal and temporal lobes [29]. Cortical volume and thickness reductions were observed in the insula, thalamus, and orbitofrontal cortex [11,13,15]. These studies generally used a single imaging modality to extract connectivity information. In the present study, we used multi-modal imaging to explore connectivity differences between HW and non-HW subjects.

Group-wise differences from connectivity analysis might be used as biomarker of important clinical variables. We attempted to predict BMI values from multi-modal imaging features derived from group-wise connectivity analysis. The aim of this study was to assess differences in connectivity between HW and non-HW subjects using multi-modal neuroimaging.

We obtained rs-fMRI and DTI data from the Human Connectome Project (HCP) database [30]. Connectivity analysis was applied to the following brain regions, which were reported by other studies to be involved in obesity and to be related to the reward system; thalamus, insula, putamen, and orbitofrontal cortex [10–15]. We identified regions that showed group-wise structural connection differences. Structural connections and functional connectivity that was correlated with structural connections were retained for further regression analyses. Regions and associated features were used to predict BMI scores in a partial least-square regression (PLSR) framework. In summary, we identified regions and features from rs-fMRI and DTI responsible for group-wise differences, and used the results to predict BMI values.

Materials and Methods

Subjects and imaging data

Institutional Review Board (IRB) of Sungkyunkwan University approved our retrospective study. Our study did not require participants’ consent as we analyzed anonymized data. Participant data were anonymized prior to analysis. We obtained DTI and rs-fMRI data from the HCP database [30]. This study included 120 subjects divided into a healthy weight (HW) group (n = 60) and a non-HW group (n = 60). Groups were classified using BMI scores. All individuals in the HW group had a BMI less than 25; the mean BMI was 22.53 kg/m2 with a standard deviation (SD) of 1.37. All individuals in the non-HW group had a BMI greater or equal to 25; the mean BMI was 32.55 kg/m2 with a SD of 3.60. BMI scores were significantly different between the two groups (p < 0.05). The mean age of the HW group was 28.93 years (SD 3.48 years), while the mean age of the non-HW group was 29.60 years (SD 3.60 years). The ratio of males and females was 24:36 in the HW group and 25:35 in the non-HW group. Age and gender were not significantly different between the two groups (p > 0.05). No subjects had mental disorders according to the diagnostic and statistical manual of mental disorders fourth edition (DSM-IV) [31].

The following imaging parameters were used for DTI on a Siemens Skyra 3T scanner: image matrix = 168 x 144; pixel resolution = 1.25 mm isotropic; b-value = 3000 s/mm2; gradient = 270; b0 images = 18; repetition time (TR) = 5520 ms; echo time (TE) = 89.5 ms; field of view (FOV) = 201x180 mm; flip angle = 78 degrees. For rs-fMRI, the following imaging parameters were used: image matrix = 104 x 90; number of slices = 1200; pixel resolution = 2.0 mm isotropic; slice thickness = 2.0 mm; TR = 720 ms; TE = 33.1 ms; FOV = 208 x 180 mm; flip angle = 52 degrees. HCP rs-fMRI data provided two different phase encodings that differed from those of the conventional phase encoding: “left-to-right” and “right-to-left.” We concatenated the two phase encodings and considered the concatenated data as one-time series data. Each rs-fMRI acquisition took approximately 15 minutes. The HCP team quality checked imaging data, including excessive head movement, as described in [30]. Data used in this study had no outliers in terms of head movement.

Image pre-processing: DTI

The HCP database provides pre-processed DTI and rs-fMRI data for many commonly performed pre-processing steps using FSL and FreeSurfer software [20,21,32]. Pre-processing steps and procedures applied to analyze the DTI data are described briefly below. The intensity of mean b0 images was normalized across six diffusion series. Echo planar (EPI) distortions were removed using the TOPUP algorithm in FSL. Distortion-corrected images were not different between the two groups. Eddy current-induced distortions and subject motion were corrected with the EDDY algorithm in FSL. All b0 images were registered to T1-weighted images using FSL’s FLIRT algorithm, and the registration results were refined using FreeSurfer’s Bbregister algorithm. Registered images were imported into standard Montreal Neurological Institute (MNI) structural space.

Image pre-processing: rs-fMRI

The following procedures were applied to the rs-fMRI data: distortions due to gradient non-linearity were corrected with the gradient_nonlin_unwarp package in FreeSurfer. Head motions were corrected by realigning all volumes to the single-band reference image acquired at the beginning of data acquisition using FSL’s FLIRT. Motion-corrected images were co-registered onto T1-weighted structural images using FLIRT in FSL and then FreeSurfer’s Bbregister for fine tuning. Co-registered images were registered onto the standard MNI space with non-linear transformation using FNIRT in FSL. Intensity normalization of total volumes was performed with a whole brain global mean of 10,000. Non-brain extraction was performed by applying the standard MNI brain mask to individual subject spaces. Temporal pre-processing of rs-fMRI data was performed using AFNI software [19]. Nuisance variables, such as six rigid motion parameters and white matter and cerebrospinal contributions were regressed out using the 3dDeconvolve package in AFNI. Band-pass filtering ranging from 0.009–0.08 Hz was applied to remove slow drift using the 3dFourier package in AFNI.

Tractography and structural network construction

Neuronal fiber densities were extracted to estimate group-wise structural connection differences between the HW and non-HW groups. Many studies have adopted the number of fibers or fiber density to effectively differentiate between healthy and diseased groups, thus we considered fiber density to be a reliable metric [33–37]. Neuronal fiber information was extracted using the Diffusion Toolkit [23]. DTI images were reconstructed with the corresponding gradient tables. Linear least-squares fitting method was used for tensor estimation. Fiber tracking was performed with the Fiber Assignment using the Continuous Tracking (FACT) algorithm, which propagates a line from the center of a seed voxel following the vector orientation, the largest eigenvector, to the next voxel. Tracking was complete when the fiber direction changed rapidly (i.e., angle threshold greater than 35°). A spline filter to smooth tracks and a cutoff filter 20–500 mm in length were applied. Connectivity analysis requires regions of interest (ROIs) to be specified so that correlations among them can be investigated. We transferred information from a pre-defined atlas onto individual subject image spaces to specify the ROIs via registration. The Automated Anatomical Labeling (AAL) atlas containing 116 ROIs was adopted [38]. Nodes were assigned to ROIs transferred from the atlas. We computed fiber density values to construct a structural network. Fiber density between two regions is the number of fiber tracks per unit surface and is computed as , where S u and S v denote the cortical surface areas of the two regions involved, respectively, and l(f) denotes the length of fiber f [36]. Edge values were fiber density values that were entered into the matrix as elements. This matrix was referred to as the structural connectivity matrix. The fiber tracking algorithm in this study and many other studies cannot distinguish between efferent and afferent connections [39–41]. We adopted a simple undirected network model. Our tractography method estimated the number of fibers starting from one brain region and ending at the other brain region. Not all brain regions are interconnected, thus we expected that some pairs of regions would not be connected.

Functional network construction

A functional network was constructed in a similar manner to that used for structural network construction. The same atlas was used for ROIs and network node specification. Each edge value was the Pearson correlation between pre-processed time-series data of two ROIs. Edge values were entered into the functional connectivity matrix, which was further normalized using Fisher’s r-to-z transformation.

Combining structural and functional connectivity

All subjects had respective structural connectivity matrices, and we searched for connections (i.e., elements in the matrix) in the whole brain that distinguished between HW and non-HW subjects (p < 0.05, corrected). The multiple comparison issue was corrected using permutation tests randomly assigning HW and non-HW cases 5,000 times to identify statistically significant structural connections [42]. Some subjects had no fiber connecting two ROIs, and thus were excluded from the analysis to identify significant connections. Such removal of subjects might cause bias, therefore we confirmed that the proportion of removed subjects was relatively low (i.e., less than 20%). Among the structural connections that differed between groups, we focused on regions in the reward system previously identified as important to obesity research [10–15]: thalamus, insula, putamen, and orbitofrontal cortex regions. The identified structural connections were further analyzed with functional connectivity. A structural connection involves two ROIs, and mean functional connectivity values of both ROI were computed. The mean functional connectivity of a region was computed as the mean of correlation values between the given region and the rest of the brain. For example, a structural connection between region A and region B were first tested with the functional connectivity of region A and then tested with the functional connectivity of region B. We adopted the mean functional connectivity value, as this is a scaled version of degree centrality. Degree centrality quantifies the importance of a given node (i.e., region) in terms of the information flow in a network [24,43,44]. If the structural connection was significantly correlated (p < 0.05) with the functional connectivity of either of the two regions, then the structural connection was retained. We focused on regions whose functional connectivity was correlated with structural connection. For example, if we identified 12 structural connections from DTI, then up to 24 ROIs were considered for functional connectivity. The identified ROIs and the associated features, structural or functional connectivity, were further used to predict BMI values. Strong functional connectivity (via rs-fMRI) does not imply actual structural connection (via DTI). We considered regions with significant functional and structural connectivity in terms of group-wise difference so that a limited number of independent variables would be entered into the regression framework.

BMI prediction

The identified ROIs and the associated features were used as independent variables in a partial least-square regression (PLSR) framework to predict BMI values [45–47]. First, each independent variable was tested using a simple linear regression model (i.e., one independent variable) in terms of its R-squared value to quantify the contribution of each feature to BMI prediction. Then, PLSR was applied using many independent variables. PLSR is especially useful in a regression framework when there are many independent variables. PLSR is a combination of principal component analysis (PCA) and multiple linear regression [45]. The PCA portion of PLSR ensures the avoidance of the multi-collinearity problem [45]. A major parameter in PLSR is the number of latent variables (LVs). We made sure that we did not overfit the data by determining the number of LVs from the predicted residual estimated sum of squares (PRESS) [45]. The prediction quality of PLSR is measured using PRESS values [45]. The number of LVs was increased until the value of PRESS no longer improved [45,47]. The error value was computed to quantify differences between the predicted BMI and actual BMI. Two types of error, root mean squared (RMS) error and percent error, were computed.

Statistical tests

Leave-one-subject-out cross validation was applied by assigning 119 (= 120–1) subjects to the training set and the remaining subject to the test set. Differences between HW subjects and non-HW subjects were assessed within the training set using permutation tests randomly assigning HW and non-HW cases 5,000 times to identify statistically significant structural connections (p < 0.05, corrected) [42]. Each identified structural connection was correlated with functional connectivity from rs-fMRI using standard Pearson’s correlation. Structural connectivity of a given ROI was retained only if the correlation was significant (p < 0.05). We quantified the abilities of the chosen ROI and associated features to explain BMI values using R-squared (for simple linear regression) and explained variance (for PLSR) statistics. The explained variance increases as the number of LVs increased. However, we avoided overfitting the data by determining the number of LVs from PRESS. All statistical analyses were performed using MATLAB (Mathworks Inc., USA).

Results

Structural connectivity and group-wise difference

Group-wise structural connection differences between HW and non-HW subjects were quantified using leave-one-subject-out cross validation. The distribution of the fiber lengths was given in the supplement (Table A and Figure A in S1 File). The mean length of the fiber was 26.2mm, and it corresponded to the typical length of the neuronal fibers [48,49]. We repeated the process 120 times to identify features with significant group-wise differences. The identified features changed very little among training sets. We report the most representative set of features in this section. We found 10 structural connections showing significant (p < 0.05, corrected) group-wise differences between HW and non-HW subjects (Table 1 and Fig 1A). The regions related to the reward system are the thalamus, insula, putamen, and orbitofrontal cortex. The structural connection between the right medial orbitofrontal cortex and right olfactory had the lowest p-value.

thumbnail

Fig 1.

(A) Structural (fiber density) connections with significant group-wise differences. (B) Structural (fiber density) connections that showed significant correlations with fMRI. Green dots and lines represent nodes and edges showing the most significant correlation between structural and functional connectivity, respectively.

https://doi.org/10.1371/journal.pone.0141376.g001

Combined connectivity and group-wise difference.

Each of the identified structural connection from the previous step was correlated with functional connectivity of two ROIs, and structural connectivity was retained only if the correlation was significant (p < 0.05). Four structural connections failed to achieve meaningful correlations with rs-fMRI connectivity, and thus were removed. Only six (of 10) structural connections were retained after correlation with functional connectivity information (Table 2 and Fig 1B). All structural connections showed a significant correlation with rs-fMRI except for the connection between the right thalamus and right caudate (Table 2 and Fig 2). The connection between the right superior orbitofrontal cortex and the right parahippocampal showed the highest correlation value. Six structural connections and 10 functional connectivity values were retained. Among 12 possible functional connectivity candidates, one was removed due to poor correlation with structural connectivity, and another were removed because it was a duplicate of a region already identified. In more detail, structural connections involving the right putamen occurred twice (row 3 in Table 2, right putamen and right hippocampus; row 4 in Table 2, right putamen and right Heschl gyrus); thus, we had only one region of functional connectivity from two structural connections. Similar to what we reported in the preceding section, the identified features changed very little from different training sets. We reported the most representative set of features in this section.

thumbnail

Table 2. Structural (fiber density) connections that showed significant correlation with fMRI.

r1 and p1 represent the correlation value and p-value of the first region of the connection listed in the leftmost column, respectively. r2 and p2 represent the correlation value and p-value of the second region of the connection, respectively.

https://doi.org/10.1371/journal.pone.0141376.t002

BMI prediction

The identified regions and the associated features from multi-modal connectivity were used as independent variables in a PLSR framework to determine BMI values using leave-one-subject-out cross validation. The number of LVs in the PLSR framework changed very little for different training sets within the context of leave-one-subject-out cross validation. The most frequently used number of LVs was nine when using all identified regions and features (from both DTI and rs-fMRI), four when using only identified regions and features from DTI, and six when using only identified regions and features from rs-fMRI. Using all identified regions and features (from both DTI and rs-fMRI), the independent variables adequately accounted for BMI values (explained variance = 0.5732). Using only identified regions and features from DTI resulted in an explained variance of 0.3591. Using only identified regions and features from rs-fMRI resulted in an explained variance of 0.3479. Each independent variable was tested using a simple linear regression model (i.e., one independent variable) in terms of R-squared value to quantify the contribution of each feature to BMI (Table 3). Structural connections between the right thalamus and right caudate, right putamen and right Heschl gyrus, and left superior orbitofrontal cortex and left amygdala contributed significantly (p < 0.05) to explaining the BMI values (Table 3 and Fig 3). Among these, the connection between the right thalamus and right caudate was the most significant factor. Functional connectivity of the right superior orbitofrontal cortex and right thalamus contributed significantly (p < 0.05) to explaining the BMI values (Table 3 and Fig 3). Among these, functional connectivity of the right thalamus was the most significant factor.

thumbnail

Table 3. Identified regions and features and their contributions to BMI values.

R-squared and p-values are reported. Features with two regions are associated with structural (fiber density) connections, and features with one region are associated with functional correlations. Each structural connection has two regions, and the first region is used as the seed region. Structural and functional features that best explained BMI are shown in italicized bold font.

https://doi.org/10.1371/journal.pone.0141376.t003

The coefficients of PLSR were computed from the training data, and the model was tested on test data. Predicted BMI values using identified regions and features and actual BMI values showed a meaningful correlation (r = 0.4414, p < 0.001; Fig 4). The mean RMS error between predicted and actual BMI was 5.26 with a SD of 5.26. The percent error mean was 15.0% with a SD of 11.1%.

Discussion

We considered multi-modal connectivity information from rs-fMRI and DTI to distinguish between HW and non-HW subjects. Significant group-wise differences in structural connections were observed in the thalamus, insula, putamen, and orbitofrontal cortex. Identified connections were correlated with functional connectivity. Statistically significant regions and associated features were used as regressors in a PLSR framework to predict BMI, a widely accepted measure for assessing obesity.

Use of multi-modal features increased the amount of BMI in terms of explained variance (combined 0.5732_vs_. rs-fMRI only 0.3479 and DTI only 0.3591). Rs-fMRI and DTI features complemented each other and better explain BMI. Our BMI prediction model was able to predict BMI values in HW and non-HW subjects. We explicitly avoided overfitting the BMI data by determining LVs of PLSR from PRESS.

Structural connections between the right thalamus and right caudate contributed most to explaining BMI among the structural factors evaluated. This connection was previously identified in research demonstrating that the thalamus and caudate comprise a brain circuit and co-activate on several tasks such as memory and inhibition [50]. Functional connectivity of the right thalamus contributed most to explaining BMI among the functional connectivity-derived factors. Previous studies demonstrated decreased functional connectivity and altered activations in the thalamus in obesity [51–53].

Our study identified brain regions whose structural and functional connectivity measures could predict BMI values well. Obesity is conventionally assessed with body fat percentage, BMI, and other similar measures. Stratifying patients using these conventional measures does not consider the rich neuroimaging data available for many obesity-related diseases, including hypertension and diabetes [2]. Stratifying patients using neuroimaging might allow us to consider brain alterations such as cognitive and executive dysfunctions so that many obesity-related diseases could be better assessed [2,54,55]. Stratifying patients using neuroimaging features (e.g., structural connection between the right thalamus and right caudate within the reward system) could also allow better characterization of obesity-related diseases. The same neuroimaging features could be used to assess the differential impacts of obesity-related hormones such as leptin and ghrelin [16,17]. Further research is needed explore the use of neuroimaging to assess obesity.

We adopted a PLSR framework instead of multiple linear regression to predict BMI values from regressors. We chose PLSR among many possible machine learning-based approaches as we felt that avoiding the multi-collinearity problem among independent variables was important [45–47]. We left exploration of support-vector regression and other machine learning approaches for future studies [56]. Other advanced algorithms exist to analyze structural and functional connectivity. A measure such as betweenness centrality could be used to identify discriminating regions and associated features [26,43]. Two flavors of MRI, DTI and rs-fMRI, were adopted in this study. Several other imaging modalities, including PET, are also sensitive to obesity-related changes [6]. The addition of another modality could provide complementary information and thus lead to better assessment of obesity-related changes.

Supporting Information

Acknowledgments

This study (H.P.) was supported in part by the Institute for Basic Science (grant number IBS-R015-D1, http://www.ibs.re.kr/). Data were provided in part by the Human Connectome Project, WU-Minn Consortium (NIH 1U54MH091657, http://www.humanconnectome.org/).

Author Contributions

Conceived and designed the experiments: BP HP. Performed the experiments: BP. Analyzed the data: HP. Contributed reagents/materials/analysis tools: BP HP JS JY. Wrote the paper: JS JY.

References

  1. 1.World Health Organization: Obesity and overweight. 2015.
  2. 2.Raji CA, Ho AJ, Parikshak N, Becker JT, Lopez OL, Kuller LH, et al. Brain Structure and Obesity. Hum Brain Mapp. 2010;31: 353–364. pmid:19662657
  3. 3.Rothman KJ. BMI-related errors in the measurement of obesity. Int J Obes. 2008;32: S56–59.
  4. 4.Burger KS, Berner LA. A functional neuroimaging review of obesity, appetitive hormones and ingestive behavior. Physiol Behav. Elsevier B.V.; 2014;136: 121–7.
  5. 5.Willeumier KC, Taylor DV, Amen DG. Elevated BMI is associated with decreased blood flow in the prefrontal cortex using SPECT imaging in healthy adults. Obesity. Nature Publishing Group; 2011;19: 1095–7.
  6. 6.Albe ARDS, Ratley RIEP, Awson MIL, Eiman ERICMR, Avussin ERICR. Neuroanatomical correlates of hunger and satiation in humans using positron emission tomography. Proc Natl Acad Sci U S A. 1999;96: 4569–74. pmid:10200303
  7. 7.Carnell S, Gibson C, Benson L, Ochner CN, Geliebter A. Neuroimaging and obesity: current knowledge and future directions. Obes Rev. 2012;13: 43–56. pmid:21902800
  8. 8.Stice E, Spoor S, Zald DH. Relation of Obesity to Consummatory and Anticipatory Food Reward. Physiol Behav. 2009;97: 551–60. pmid:19328819
  9. 9.Stice E, Spoor S, Bohon C. Relation of Reward from Food Intake and Anticipated Food Intake to Obesity: A Functional Magnetic Resonance Imaging Study. J Abnorm Psychol. 2008;117: 924–35. pmid:19025237
  10. 10.Opel N, Redlich R, Grotegerd D, Dohm K, Haupenthal C, Heindel W, et al. Enhanced neural responsiveness to reward associated with obesity in the absence of food-related stimuli. Hum Brain Mapp. 2015;36: 2330–2337. pmid:25704752
  1. 11.Marqués-Iturria I, Pueyo R, Garolera M, Segura B, Junqué C, García-García I, et al. Frontal cortical thinning and subcortical volume reductions in early adulthood obesity. Psychiatry Res Neuroimaging. Elsevier; 2013;214: 109–115.
  1. 12.Carnell S, Gibson C, Benson L, Ochner CN, Geliebter A. Neuroimaging and obesity: current knwoledge and future directions. Obes Rev. 2012;13: 43–56. pmid:21902800
  1. 13.Shott ME, Cornier M, Mittal VA, Pryor TL, Orr JM, Brown MS, et al. Orbitofrontal cortex volume and brain reward response in obesity. Int J Obes. Nature Publishing Group; 2015;39: 214–221.
  1. 14.Stoeckel LE, Weller RE, Cook EW, Twieg DB, Knowlton RC, Cox JE. Widespread reward-system activation in obese women in response to pictures of high-calorie foods. Neuroimage. 2008;41: 636–47. pmid:18413289
  1. 15.Opel N, Redlich R, Grotegerd D, Dohm K, Heindel W, Kugel H, et al. Obesity and major depression: Body-mass index (BMI) is associated with a severe course of disease and specific neurostructural alterations. Psychoneuroendocrinology. Elsevier Ltd; 2015;51: 219–226.
  1. 16.Monteleone P, Maj M. Dysfunctions of leptin, ghrelin, BDNF and endocannabinoids in eating disorders: Beyond the homeostatic control of food intake. Psychoneuroendocrinology. Elsevier Ltd; 2013;38: 312–330.
  1. 17.Meier U, Gressner AM. Endocrine regulation of energy metabolism: review of pathobiochemical and clinical chemical aspects of leptin, ghrelin, adiponectin, and resistin. Clin Chem. 2004;50: 1511–1525. pmid:15265818
  1. 18.Barbas H. Flow of information for emotions through temporal and orbitofrontal pathways. J Anat. 2007;211: 237–49. pmid:17635630
  1. 19.Cox RW. AFNI: Software for Analysis and Visualization of Functional Magnetic Resonance Neuroimages. Comput Biomed Res. 1996;29: 162–73. pmid:8812068
  1. 20.Fischl B. FreeSurfer. Neuroimage. 2012;62: 774–781. pmid:22248573
  1. 21.Jenkinson M, Beckmann CF, Behrens TEJ, Woolrich MW, Smith SM. Fsl. Neuroimage. 2012;62: 782–790. pmid:21979382
  1. 22.Daducci A, Gerhard S, Griffa A, Lemkaddem A, Cammoun L, Gigandet X, et al. The connectome mapper: an open-source processing pipeline to map connectomes with MRI. PLoS One. 2012;7: 1–9.
  1. 23.Wang R, Benner T, Sorensen AG, Wedeen VJ. Diffusion Toolkit: A Software Package for Diffusion Imaging Data Processing and Tractography. Proc Intl Soc Mag Reson Med. 2007;15: 3720.
  1. 24.Bullmore E, Sporns O. Complex brain networks: graph theoretical analysis of structural and functional systems. Neuroscience. 2009;10: 186–198. pmid:19190637
  1. 25.Watts D, Strogatz S. Collective dynamics of “small-world”networks. Nature. 1998;393: 440–442. Available: http://www.nature.com/nature/journal/v393/n6684/abs/393440a0.html pmid:9623998
  1. 26.Stam CJ, Reijneveld JC. Graph theoretical analysis of complex networks in the brain. Nonlinear Biomed Phys. 2007;1: 1–19.
  1. 27.García-García I, Jurado MÁ, Garolera M, Segura B, Sala-Llonch R, Marqués-Iturria I, et al. Alterations of the salience network in obesity: a resting-state fMRI study. Hum Brain Mapp. 2013;34: 2786–97. pmid:22522963
  1. 28.Kullmann S, Heni M, Veit R, Ketterer C, Schick F, Häring H-U, et al. The obese brain: association of body mass index and insulin sensitivity with resting state network functional connectivity. Hum Brain Mapp. 2012;33: 1052–61. pmid:21520345
  1. 29.Bolzenius JD, Laidlaw DH, Cabeen RP, Conturo TE, McMichael AR, Lane EM, et al. Brain structure and cognitive correlates of body mass index in healthy older adults. Behav Brain Res. Elsevier B.V.; 2015;278: 342–347.
  1. 30.Van Essen DC, Smith SM, Barch DM, Behrens TEJ, Yacoub E, Ugurbil K. The WU-Minn Human Connectome Project: an overview. Neuroimage. Elsevier Inc.; 2013;80: 62–79.
  1. 31.Association AP. Diagnostic and Statistical Manual of Mental Disorders. 4th ed. American Psychiatric Press, Washington, DC; 1994.
  2. 32.Glasser MF, Sotiropoulos SN, Wilson JA, Coalson TS, Fischl B, Andersson JL, et al. The minimal preprocessing pipelines for the Human Connectome Project. Neuroimage. Elsevier Inc.; 2013;80: 105–124.
  1. 33.Lee S-H, Seo J, Lee J, Park H. Differences in early and late mild cognitive impairment tractography using a diffusion tensor MRI. Neuroreport. 2014;25: 1393–1398. pmid:25325351
  1. 34.Lo C-Y, Wang P-N, Chou K-H, Wang J, He Y, Lin C-P. Diffusion tensor tractography reveals abnormal topological organization in structural cortical networks in Alzheimer’s disease. J Neurosci. 2010;30: 16876–16885. pmid:21159959
  1. 35.Zhou Y, Dougherty JH, Hubner KF, Bai B, Cannon RL, Hutson RK. Abnormal connectivity in the posterior cingulate and hippocampus in early Alzheimer’s disease and mild cognitive impairment. Alzheimer’s Dement. 2008;4: 265–270.
  1. 36.Hagmann P, Cammoun L, Gigandet X, Meuli R, Honey CJ, Van Wedeen J, et al. Mapping the structural core of human cerebral cortex. PLoS Biol. 2008;6: 1479–1493.
  1. 37.Torgerson CM, Irimia A, Leow AD, Bartzokis G, Moody TD, Jennings RG, et al. DTI tractography and white matter fiber tract characteristics in euthymic bipolar I patients and healthy control subjects. Brain Imaging Behav. 2013;7: 129–139. pmid:23070746
  1. 38.Tzourio-Mazoyer N, Landeau B, Papathanassiou D, Crivello F, Etard O, Delcroix N, et al. Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage. 2002;15: 273–289. pmid:11771995
  1. 39.Gong G, He Y, Concha L, Lebel C, Gross DW, Evans AC, et al. Mapping anatomical connectivity patterns of human cerebral cortex using in vivo diffusion tensor imaging tractography. Cereb Cortex. 2009;19: 524–536. pmid:18567609
  1. 40.Yan C, Gong G, Wang J, Wang D, Liu D, Zhu C, et al. Sex- and brain size-related small-world structural cortical networks in young adults: A DTI tractography study. Cereb Cortex. 2011;21: 449–458. pmid:20562318
  1. 41.Iturria-Medina Y, Sotero RC, Canales-Rodríguez EJ, Alemán-Gómez Y, Melie-García L. Studying the human brain anatomical network via diffusion-weighted MRI and Graph Theory. Neuroimage. 2008;40: 1064–1076. pmid:18272400
  1. 42.Smith SM, Vidaurre D, Beckmann CF, Glasser MF, Jenkinson M, Miller KL, et al. Functional connectomics from resting-state fMRI. Trends Cogn Sci. 2013;17: 666–682. pmid:24238796
  1. 43.Xia S, Foxe JJ, Sroubek AE, Branch C, Li X. Topological organization of the “small-world” visual attention network in children with attention deficit/hyperactivity disorder (ADHD). Front Hum Neurosci. 2014;8: 1–14.
  1. 44.Rubinov M, Sporns O. Complex network measures of brain connectivity: uses and interpretations. Neuroimage. Elsevier Inc.; 2010;52: 1059–1069.
  1. 45.Krishnan A, Williams LJ, McIntosh AR, Abdi H. Partial Least Squares (PLS) methods for neuroimaging: a tutorial and review. Neuroimage. Elsevier Inc.; 2011;56: 455–475.
  1. 46.Abdi H. Partial least squares regression and projection on latent structure regression (PLS Regression). Wiley Interdiscip Rev Comput Stat. 2010;2: 97–106.
  1. 47.Park H, Yang J, Seo J, Choi Y, Lee K, Lee J. Improved explanation of human intelligence using cortical features with second order moments and regression. Comput Biol Med. Elsevier; 2014;47: 139–146.
  1. 48.Gao J, Cheung RTF, Chan YS, Chu LW, Mak HKF, Lee TMC. The relevance of short-range fibers to cognitive efficiency and brain activation in aging and dementia. PLoS One. 2014;9: 1–9.
  1. 49.Englander ZA, Pizoli CE, Batrachenko A, Sun J, Worley G, Mikati MA, et al. Diffuse reduction of white matter connectivity in cerebral palsy with specific vulnerability of long range fiber tracts. NeuroImage Clin. The Authors; 2013;2: 440–447.
  1. 50.Hart J, Maguire MJ, Motes M, Mudar RA, Chiang HS, Womack KB, et al. Semantic memory retrieval circuit: role of pre-SMA, caudate, and thalamus. Brain Lang. Elsevier Inc.; 2013;126: 89–98.
  1. 51.Musen G, Jacobson AM, Bolo NR, Simonson DC, Shenton ME, McCartney RL, et al. Resting-state brain functional connectivity is altered in type 2 diabetes. Diabetes. 2012;61: 2375–2379. pmid:22664957
  1. 52.Matsuda M, Liu Y, Mahankali S, Pu Y, Mahankali A, Wang J, et al. Altered hypothalamic function in response to glucose ingestion in obese humans. Diabetes. 1999;48: 1801–06. pmid:10480611
  1. 53.Gautier JF, Chen K, Salbe AD, Bandy D, Pratley RE, Heiman M, et al. Differential brain responses to satiation in obese and lean men. Diabetes. 2000;49: 838–46. pmid:10905495
  1. 54.Dahl AK, Hassing LB. Obesity and cognitive aging. Epidemiol Rev. 2013;35: 22–32. pmid:23258415
  1. 55.Liang J, Matheson BE, Kaye WH, Boutelle KN. Neurocognitive correlates of obesity and obesity-related behaviors in children and adolescents. Int J Obes. Nature Publishing Group; 2014;38: 494–506.
  1. 56.Mourão-Miranda J, Bokde A LW, Born C, Hampel H, Stetter M. Classifying brain states and determining the discriminating activation patterns: Support Vector Machine on functional MRI data. Neuroimage. 2005;28: 980–995. pmid:16275139