Impact of Hunger, Satiety, and Oral Glucose on the Association Between Insulin and Resting-State Human Brain Activity

Highlights

- Brain areas relevant for ingestive behavior are modulated as a function of hunger/satiety.

- Metabolic state, insulin levels and glucose administration interact in their effects on brain activation patterns.

- FALFF is a reliable index of spontaneous brain activity.

Introduction

Food ingestion and energy homeostasis are regulated by central nervous pathways (Schwartz and Porte, 2005; Lam, 2010; Suzuki et al., 2010) and modulated by endocrine signals (Gao and Horvath, 2008). Insulin is one of the hormones that is part of a negative feedback loop that ensures balanced energy homeostasis (Schwartz et al., 2000; Schwartz, 2001; Obici et al., 2002). Of note, both intraventricular injection and intranasal administration of insulin decreases food intake and body weight in rodents (Chavez et al., 1995; Brown et al., 2006) and humans (Benedict et al., 2008). In contrast, inactivation of insulin leads to opposite effects (McGowan et al., 1990; Brüning et al., 2000).

Task-based functional magnetic resonance imaging (fMRI) has been used to study the impact of insulin and glucose on brain responses to food-related stimuli under different homeostatic conditions on the network level (Woods et al., 1998; Liu et al., 2000) linking plasma glucose (Kroemer et al., 2013b) and fasting insulin levels (Wallner-Liebmann et al., 2010) to neural activity changes in the hypothalamus, thalamus, amygdala, insula, and superior frontal cortex, brain sites assumed to regulate appetite-related eating behavior (Kullmann et al., 2013). According to Smitha et al. (2017), task-based designs focus on a small fraction of the brain’s overall activity only. To avoid these disadvantages in the present investigation, we used resting-state fMRI in order to reveal potential links between whole brain activity and insulin and glucose levels in different metabolic states.

Resting-state fMRI (rs-fMRI) is a task-free paradigm in which participants do not perform any specific task (Biswal et al., 1995; van den Heuvel and Hulshoff Pol, 2010; Lee et al., 2013). Several studies using rs-fMRI revealed associations of insulin or glucose levels with the functional connectivity (FC) of particular brain networks or brain sites that are related to homeostatic regulation but not without contradictions. Some studies reported an increase of resting-state FC within the limbic system in response to the administration of 75 g oral glucose after overnight fasting. This increase in FC was positively correlated with plasma insulin levels (Kullmann et al., 2012; Wölnerhanssen et al., 2015). Conflicting results were found by Page et al. (2013) using MRI-cerebral blood flow (CBF) during rest: here, changes of insulin levels were negatively associated with changes of CBF signals in the caudate and the putamen in response to glucose administration. Similarly, changes in plasma insulin levels in response to a meal after overnight fasting were negatively correlated with changes of CBF signals in the insula and prefrontal cortex (Tataranni et al., 1999). Additionally, the FC between the posterior insula (PINS) and superior frontal gyrus (SFG) under hunger conditions was partially moderated by the plasma glucose levels, indicating that the PINS connectivity depends on the homeostatic energy deficit caused by fasting (Wright et al., 2016). Furthermore, a study with experimentally induced hypoglycemia reported increased FC of the default mode network (DMN) with posterior cingulate cortex (PCC) and decreased FC of insula, superior and inferior frontal gyri with temporal networks, basal ganglia, and cerebellum in healthy subjects (Bolo et al., 2015). The inconsistencies reported here might be related to the different experimental paradigms and neuroimaging modalities. In most of these studies, FC parameters were computed by using seed-based analysis or independent component analysis (ICA) to define brain networks of interest or to decompose the brain into multiple networks, respectively. This means that FC studies rely on the correlations and ignore the changes in regional brain activity under different metabolic conditions. However, it has been shown that there is a linear relationship between the amplitude of the blood-oxygen-level-dependent (BOLD) signal and brain metabolism (Tomasi et al., 2013). Therefore, in our study, we used the amplitude of the rs-fMRI signal to investigate the sensitivity of resting-state brain activity (Zou et al., 2008) to changes in metabolic states. A critical question of the present search is, whether or not the changes in the amplitude of the BOLD signal of specific brain regions are associated with changes in hormonal signals, such as insulin.

To clarify some of the issues discussed in the previous paragraph, we investigated: (i) whether changes in brain regions linked to hypoglycemia, such as the DMN, are modulated by insulin and glucose or different metabolic states; and (ii) whether changes in activity of specific brain regions, such as hypothalamus and insula, are modulated by changes of blood insulin or glucose levels after glucose administration. In the present investigation, we used a voxel-wise frequency-domain approach to measure whole brain activity by using amplitude (intensity) values of BOLD signals called fractional amplitude of low-frequency fluctuation (fALFF) for each individual per condition (Zou et al., 2008). The fALFF allow us to study the local spontaneous brain activity across the whole brain based on the magnitude of the BOLD signal in the low-frequency range (Zou et al., 2008; Egorova et al., 2017) which is higher in gray than in white matter (Zuo et al., 2010). The fALFF approach efficiently suppresses non-specific signal components, such as physiological noise (Zou et al., 2008; Cole et al., 2010). Other studies have observed that fALFF is associated with body mass index (BMI) after intranasal insulin application (Kullmann et al., 2013).

In a recent rs-fMRI study (Al-Zubaidi et al., 2018), we examined the activity and connectivity brain responses to the interaction of metabolic state (hungry vs. satiated) and glucose administration (before vs. after administration of 75 g of oral glucose). In that study, we used multimethod rs-fMRI approaches to identify brain activation patterns that are associated with changes in metabolic states and caloric intake. We showed that in contrast to other voxel-wise analyses like regional homogeneity or degree of centrality, fALFF is a more sensitive metric for identifying differences in the resting brain activity, for example, the amplitude of the SFG and PCC were increased after oral glucose treatment and in hunger conditions, respectively. However, our previous study focused on the effect of glucose treatment solely on brain activity and connectivity without taking the effect of plasma glucose and insulin levels into account when analyzing neuroimaging data (Al-Zubaidi et al., 2018). To get a better insight into hormone-dependent brain activity at rest, in the present study we used multiple linear regression analysis (MLRA) with covariates (plasma insulin and glucose levels) to investigate the dependencies between brain activity, measured as the amplitude of the BOLD signal, and hormone levels. By manipulating metabolic state (hunger/satiety) and glucose administration (before and after 75 g of oral glucose) in a factorial design, we tested the hypothesis that brain areas involved in homeostatic regulation are modulated by peripheral circulating insulin and glucose levels (Williams et al., 2000; Figlewicz, 2003; Zanchi et al., 2017).

Materials and Methods

Participants

We investigated 24 healthy male volunteers aged 20–30 years (mean ± SEM: 24.5 ± 0.6 years) with a BMI of 20–25 kg/m² (mean ± SEM in kg/m²: 23.4 ± 0.3) recruited from the local university community. All participants underwent a medical interview before participating in this study and examination assessing general health, medicines, mood, blood glucose concentration and cognitive disorders. More details on the exclusion criteria of participants can be found in Al-Zubaidi et al. (2018). Within 4 weeks before and during the trials, subjects were instructed not to participate in other studies or to donate blood. The study was approved by the ethics committee of the University of Lübeck and carried out under the Declaration of Helsinki (2000). Before participation, each participant gave written informed consent.

Experimental Design

Each subject was investigated twice, once while being in a hungry (36 h fasting) and once while being in a satiated state (standardized eating, three meals/day for 36 h). The order of the two sessions was counterbalanced across participants, with a break of at least 1 week between sessions. To control physical activity and food intake, subjects stayed at the sleep laboratory of the Department of Psychiatry of the University Hospital of Lübeck (UKSH Lübeck) for 36 h in both conditions. At the end of each 36 h stay, the functional MRI session took place. Each MRI session started with the first resting state recording. Subjects were removed from the scanner to drink a solution containing the equal of 75 g of glucose (determined by a 300 ml mixture of mono and oligosaccharides; ACCU-CHEK^® Dextro^® O.G-T., Roche, Grenzach-Wyhlen, Germany). To avoid an impact of any circadian variations the timing of the glucose drink was kept constant and started each time at 01:25 pm. Twenty minutes after glucose intake, the second resting state fMRI was recorded.

In the hunger condition, subjects began fasting (no food or beverages, except water) from 11:00 pm the night before the examination started until the end of the second resting state recording. In the satiety condition, three standardized meals per day were provided. Standardized meals were served according to recommendations of the clinical diabetes counseling department at the University Medical Campus Schleswig-Holstein: Breakfast (25% protein, 50% carbohydrate and 25% fat), lunch (20% protein, 63% carbohydrate and 17% fat) and dinner (22% protein, 60% carbohydrate and 18% fat) were provided at 09:00 am, 12:00 pm and 07:00 pm, respectively. In both sessions, subjects arrived at the sleep lab at 08:00 am. Nurses inserted a cannula into a peripheral vein on the back of the hand to collect blood samples. The first blood sample for defining the basal blood glucose and insulin levels was taken at 08:45 am. Then the experimental protocol of the first day was completed and nine blood samples were collected across different time points. More details on the exact timing for obtaining blood samples can be found in our previous work (please see the “Experimental design” section and Figure 1 in Al-Zubaidi et al., 2018). All individuals stayed and slept overnight in the sleep laboratory. The next morning, blood samples were drawn at 08:45 am and at fixed time points throughout noon until the MRI sessions at 01:05 pm. For the resting state fMRI recording (duration 6 min), subjects were instructed to lie still inside the scanner with their eyes closed and to not engage in any particular cognitive activity. After the MRI measurements, subjects were brought back to the sleep lab and provided blood samples every 30 min until 04:45 pm on the same day. To investigate the changes in plasma glucose and insulin levels and to relate them to the fluctuation of resting brain activity, the four blood samples that were collected before i.e., 20, 100, 205 and 280 min before 75 (g oral dose of glucose) and the six blood samples that were taken after (20, 50, 80, 110, 140 and 170 min after glucose intake) oral glucose administration were used for the analyses (Figure 1). For each condition subjects rated how hungry they were 20 min before and 20 min after oral glucose administration and this was carried out using a visual analog scale from 0 (not hungry at all) to 9 (very hungry).

FIGURE 1

Figure 1. Time course of physiological parameters under hunger and satiety conditions. (A) Mean plasma concentration of insulin and (B) glucose under hunger and satiety conditions. Glucose was ingested at time point 0. The time-points of the first and second functional magnetic resonance imaging (fMRI) investigations are marked in gray. Boxes on the bottom of the graph indicate the time points of meals on the second day (B, breakfast; L, lunch at 09:00 am and 12:00 pm, respectively). The error bars represent the standard deviation.

Image Acquisition

We used a 3-T Philips Achieva scanner (Philips Medical Systems, Best, Netherlands) and a standard eight-channel phased array head coil to record functional and structural images. Functional images were acquired with a gradient-echo-EPI sequence in ascending interleaved fashion (repetition time TR = 2,000 ms; echo time TE = 28 ms; voxel dimensions 3 × 3 × 3 mm; field of view 192 × 192 mm²; 80° flip angle; 40 slices). 178 whole-brain volumes were recorded for each of the two runs per session. Further, the high-resolution structural T1 image consisted of 180 sagittal slices were acquired by applying a T1-weighted 3D turbo gradient-echo sequence with SENSE (image matrix 240 × 240; voxel dimensions 1 × 1 × 1 mm; field of view 240 × 240 mm²; 1 mm slice thickness; 9° flip angle).

Handling and Analyses of Blood Samples, Hunger Ratings

Plasma glucose concentration was measured with B-Glucose-Data-Management system (HemoCue GmbH, Grossostheim, Germany). The blood samples were, if necessary, cooled and immediately centrifuged. The supernatants were stored at −80°C until analyses. Serum insulin concentration was measured by commercial enzyme-linked immunoassay (Immulite DPC, Los Angeles, USA; insulin: intra-assay coefficient of variation (CV) <1.5% and inter-assay CV <4.9%).

The average plasma levels of glucose and insulin were calculated for each subject under hunger and satiety conditions before and after glucose treatment (Figure 2). To test the interaction between metabolic state (two levels: hunger, satiety) and the effect of glucose administration (two levels: before, and after treatment) on plasma glucose, plasma insulin, and subjective hunger rating, we applied a two-way repeated measures ANOVA for each of the dependent variables separately (SPSS, Version 22.0). Values are reported as mean (M) and standard deviation (± SD).

FIGURE 2

Figure 2. Statistical differences in physiological and behavioral parameters under hunger and satiety conditions. (A) Mean of plasma glucose, (B) plasma insulin, and (C) hunger ratings per factor: metabolic states (two levels: hunger and satiety) and oral glucose treatment (two levels: before, and after treatment) across participants. The error bars represent the standard deviation. ***Represent the significant differences (p < 0.0001).

We also ran a Pearson correlation analysis to investigate the associations between changes of the subjective feeling of hunger with changes plasma glucose and insulin as well as with changes in brain activity (i.e., clusters that survived the cluster-significance test).

Functional Image Preprocessing and Analyses

Part of preprocessing on functional images was initially carried out using FSLv5.0 (available at http://fsl.fmrib.ox.ac.uk/fsl/) to implement ICA-based strategy for automatic removal of motion artifacts (ICA-AROMA) for head motion correction (Pruim et al., 2015a). It has been shown that ICA-AROMA enhances the sensitivity and specificity of rs-fMRI activation and connectivity analyses (Pruim et al., 2015a). To improve inter-subject alignment (Klein et al., 2009; McLaren et al., 2010), the spatial preprocessing of the data were performed with the statistical parametric mapping 12b (SPM12b; available at http://www.fil.ion.ucl.ac.uk/spm/) in MATLAB (MathWorks, Natick, MA, USA) and data processing assistant for resting-state fMRI toolbox (DPARSF advanced edition, version 3.2, available at http://rfmri.org/DPARSF).

The rs-fMRI images were preprocessed as the following: (i) we discarded the first seven functional volumes from each participant’s four runs to allow steady-state tissue magnetization; (ii) we manually reoriented all functional volumes to the anterior commissure; (iii) we implemented head movement correction during data acquisition by volume-realignment to the middle volume using MCFLIRT (Jenkinson et al., 2002). The MCFLIRT results in estimated realignment parameters for each experimental condition were summarized in Figure 3 for motion translation and Supplementary Figure S1 for motion rotation parameters; (iv) we applied ICA-AROMA to the data in order to identify and remove motion-related components using four spatial and temporal features as summarized in Al-Zubaidi et al. (2018) and detailed in Pruim et al. (2015b). Then the denoised functional and structural images were further preprocessed with DPARSF toolbox as following: (v) we co-registered the T1 structural image to the mean functional image; (vi) we ran a segmentation protocol to distinguish gray matter, white matter and cerebrospinal fluid (CSF); (vii) we applied bias correction and spatial normalization of the T1 structural image and adjusted them to the Montreal Neurological Institute (MNI) template using diffeomorphic anatomical registration through exponentiated Liealgebra (DARTEL) algorithm (Ashburner and Friston, 2005); (viii) we performed nuisance regression (including white matter and CSF signals) to reduce the influence of unspecific physiological effects on BOLD signal (Liu, 2016); (ix) we spatially normalized functional images to the MNI-template using the normalization parameters estimated by DARTEL algorithm with a 3 mm isotropic voxels size; (x) we performed spatial smoothing with a 6 mm full width at half maximum (FWHM) Gaussian kernel; and (xi) finally, we masked all functional images with a gray matter mask, which was calculated by averaging the gray matter masks of all subjects (Al-Zubaidi et al., 2018).

FIGURE 3

Figure 3. Summary of MCFLIRT estimated translation parameters for each experimental condition.

We carried out a first-level analysis on each subject and each run by calculating the fALFF on the preprocessed data using DPARSF toolbox (Yan and Zang, 2010). DPARSF has built-in Fast Fourier Transform (FFT) to convert time series for each voxel to the frequency domain and compute the power spectrum. This procedure estimates the amplitude of each frequency as the square root of the power spectrum. Last, the total amplitude of the low-frequency range (0.01–0.08 Hz) was divided by that of the entire frequency range 0–0.25 Hz (TR = 2 s). Put simply, fALFF is the ratio of rs-fMRI signal fluctuation in the low-frequency range relative to the entire frequency range (Zou et al., 2008). Individual fALFF maps were transformed to z-scores to reduce the effect of individual variability (Al-Zubaidi et al., 2018).

Statistical Analysis of Functional Images

To capture the modulatory impact of glucose administration on the association between the activity of brain regions and plasma insulin or glucose levels under different metabolic states, we performed two multiple linear regression analyses (MLRA) using SPM12, the first “before” and the second “after” oral glucose treatment. Each MLRA was designed with two covariates (plasma insulin and glucose levels). Briefly, MLRA is used to describe how a “linear combination” of multiple variables, called independent or explanatory variables, to predict a single response variable, referred to as the dependent or target variable. The contribution of each independent variable to the model is assessed by looking at the regression coefficients (Nathans et al., 2012). In this study, we used MLRA to figure out the contribution of plasma insulin (an independent variable) to the prediction of brain activity (the dependent variable) when taking the effect of plasma glucose (another independent variable) into account (Kiebel and Holmes, 2007; Field, 2014), and vice versa. Accordingly, our multiple linear regression model is:

$Y_j=\mu +x_{j1}{\beta }_1+x_{j2}{\beta }_2+{\in }_j,(1)$

where Y_j is a dependent variable (such as fALFF at a particular voxel) and j = 1,…,J indexes the observation. The regression coefficient μ represents a constant term (the mean of the dependent variable when all predictors are zero), while β₁ represents the regression slope, which quantifies the association of Y with x₁ (such as plasma insulin), adjusting for the effect of x₂ (plasma glucose) on Y and vice versa for β₂ and ∈ is the error associated with the regression (the variance of the dependent variable from its mean when all predictors are zero). The parameters were estimated by using the least squares method. To find significant voxels whose activity was affected by hunger vs. satiety or by glucose or insulin, we used one-sample t-tests for each regression coefficient on the voxel level per MLRA (see section below). The resulting statistics indicate whether the regression coefficient of a particular voxel is significantly different from the error in that estimate (Field, 2014). To correct for multiple comparisons, the topological false discovery rate (FDR) q = 0.05 was employed with a cluster defining threshold of p < 0.001 for the t-tests (Chumbley et al., 2010).

In the 1st MLRA (before glucose administration), we contrasted hunger vs. satiety states while controlling the moderating fluctuation of glucose and insulin. We calculated differences [delta (Δ) = hunger-satiety] of fALFF maps before oral glucose treatment. For glucose and insulin, the area under the concentration-time curve (AUC) of the first four samples (Figures 1A,B) was calculated and the Δ of the AUC was obtained. The AUC has been used as an index to assess the regulation of glucose and insulin (Tzagournis and Skillman, 1970; Owen et al., 2012). It was calculated using the standard trapezoid method, which is computing the AUC with respect to ground (see formula 2 in Pruessner et al., 2003). The group level analysis was performed using Δ AUC of glucose and insulin as covariates (independent variables) and the Δ fALFF maps as the dependent variable in the MLRA. Notably, the Δ AUC of glucose and insulin were not significantly correlated (|r| = 0.3, p = 0.1). To check whether the AUC of plasma glucose and insulin influence our findings, we also recomputed the 1st MLRA by including the Δ plasma glucose and insulin based on the value 20 min before glucose ingestion (Figure 1) as independent variables. The results of this model (see Supplementary Table S1 and Supplementary Figure S2 in the supplementary material) were similar to the results of the model with AUC (see Table 1 and Figure 4A). The AUC provides an overview of plasma glucose and insulin level profiles under diet or standard meal vs. time (Johnson et al., 2006). Also, we believe that the changes in brain activity before glucose ingestion may be related to profile change more than single glucose and insulin values. Therefore, we will report the AUC model results only.

TABLE 1

Table 1. Changes and associations of fractional amplitude of low-frequency fluctuation (fALFF) with food conditions and hormone levels.

FIGURE 4

Figure 4. Hunger vs. satiety effects on brain activity. (A) Mean fractional amplitude of low-frequency fluctuation (fALFF) value of all voxels of a significant cluster per condition across participants. (B) Results of the first model (before glucose administration). (C) Mean fALFF value of all voxels of a significant cluster per condition across participants. (D) Results of the second model (after glucose administration). Abbreviations: PCC, posterior cingulate cortex; PCUN, precuneus; IPG, inferior parietal gyrus.

The 2nd MLRA (after glucose administration) had a similar design as first MLRA. Differences (Δ) of amplitude rs-fMRI signals were calculated by subtracting the fALFF maps of hunger condition from satiety condition. The Δ calculated for glucose and insulin was based on the value 20 min after glucose ingestion (Figure 1) only to keep the temporal association of endocrine and functional imaging data as clear as possible. Δ glucose and Δ insulin were used as covariates in the MLRA. Again, glucose and insulin were not significantly correlated (|r| = 0.14, p = 0.5). Additional glucose and insulin samples collected after oral glucose intake were not included in this analysis, because they were taken after the second rs-fMRI recording.

To ensure that the correlation between specific brain regions and glucose or insulin values is not biased (Esterman et al., 2010), we extracted fALFF values by averaging across voxels in each cluster that survived the cluster-significance test. Then, we performed the leave-one-out cross-validation (LOOCV) of Pearson correlation to evaluate the relationship between averaged ΔfALFF values of each brain region with Δ of the plasma glucose and insulin. Finally, we performed full and partial correlation analyses between clusters linked to plasma insulin levels to investigate the association between them and the effects of glucose and insulin values on that association.

Finally, to investigate the acute effect of glucose administration on the interaction between brain activity and physiological changes, we performed two additional MLRA. The 3rd and 4th MLRAs estimate the correlations between changes in fALFF (Δ fALFF = pre-post glucose ingestion) and changes in plasma glucose/insulin (Δ of plasma glucose/insulin were calculated based on the value of 20 min before and after glucose ingestion; Figure 1) under satiety and hunger conditions separately, respectively.

We used the automated anatomical labeling (AAL) atlas (Tzourio-Mazoyer et al., 2002) included in the xjView toolbox¹ to label the anatomical location of significant clusters. Thalamic nuclei were identified by applying the MNI-based Morel Atlas (Jakab et al., 2012).

Results

Physiological and Behavioral Effects

The ANOVA of plasma glucose concentration (Figure 2A) revealed a significant metabolic state (hunger/satiety) and glucose administration interaction (F_(1,23) = 53.6, p < 0.001). The main effects metabolic state (F_(1,23) = 23.1, p < 0.001) and glucose treatment (F_(1,23) = 256, p < 0.001) were significant as well. The interaction was driven by an increased area under curve for glucose following the exogenous glucose challenge in the hunger (M = 7.6 mmol/l, SD = 1.1) compared to the satiety condition (M = 5.9 mmol/l, SD = 0.6) whereas this was reversed before glucose treatment with higher plasma glucose level in the satiety (M = 4.6 mmol/L, SD = 0.5) compared to the hunger condition (M = 3.9 mmol/L, SD = 0.5).

For insulin, a significant metabolic state and glucose treatment interaction was found (F_(1,23) = 101.8, p < 0.0001). In addition, a significant main effect glucose treatment (F_(1,23) = 106, p < 0.001) was revealed, but the main effect metabolic state was not significant (F_(1,23) = 0.65, p = 0.4). As shown in Figure 2B, before glucose treatment the plasma insulin level for the satiated state (M = 198.8 pmol/L, SD = 65.2) was higher compared to the hunger state (M = 22.5 pmol/L, SD = 10.1). In contrast, after glucose treatment, the plasma insulin level was greater in the hunger (M = 457.9 pmol/L, SD = 185.6) compared to the satiated state (M = 314.6 pmol/L, SD = 123.5), indicating reduced responsiveness to circulating insulin.

The analysis of the subjective hunger ratings revealed higher hunger ratings in the fasting condition compared to satiety (F_(1,23) = 28.9, p < 0.001), confirming the success of our fasting treatment (Figure 2C). There was neither a main effect of glucose treatment and nor an interaction between both factors. In addition, we found no significant correlations between subjective feeling of hunger with changes plasma glucose and insulin as well as with changes brain activity (i.e., clusters that survived the cluster-significance test). These results are shown in Supplementary Figure S3.

Resting State fMRI: Hunger vs. Satiety Effects

The 1st MLRA (before glucose administration) revealed differences of resting brain activity as a function of hunger in the PCC and the left anterior precuneus (PCUN), as well as in the left inferior parietal gyrus (IPG; Figures 4A,B, Table 1). The 2nd MLRA (after glucose administration) showed an increase of the fALFF signal in the hunger compared to the satiety condition in the left thalamus and the left PCC (Figures 4C,D, Table 1).

Resting State fMRI: Before vs. After Glucose Treatment

The 3rd (under satiety condition) and 4th (under hunger condition) MLRA results showed that the resting brain activity of supplementary motor area (SMA) was significantly decreased after glucose ingestion (Figures 5A,B, Table 2). In the hunger condition (4th MLRA), a reduced resting brain activity was evident in precentral gyrus (PreCG) and the postcentral gyrus (PoCG) after glucose administration (Figure 5B, Table 2).

FIGURE 5

Figure 5. Before vs. after glucose administration on brain activity. (A) Results of the third model (under satiety condition). (B) Results of the fourth model (under hunger condition). Abbreviations: SMA, supplementary motor area; PreCG, precentral gyrus; PoCG, postcentral gyrus.

TABLE 2

Table 2. Influences of oral glucose administration on brain activity under hunger and satiety conditions.

Correlations Between Physiological and Neural Effects

For the 1st, 3rd and 4th MLRAs, we found no correlation between changes in resting brain activity and changes in glucose and insulin passing the correction for multiple comparisons.

For 2nd MLRA (i.e., hunger-satiety), after glucose administration and while controlling for effects of Δ insulin we found a correlation between Δ fALFF and Δ plasma glucose in the left fusiform gyrus (Figures 6A,B, Table 1). A LOOCV between the average Δ fALFF values of the fusiform cluster and Δ plasma glucose after data being adjusted for Δ plasma insulin level (Figure 6C) revealed a significant positive correlation (r = 0.75).

FIGURE 6

Figure 6. Glucose-associated low-frequency BOLD fluctuations. (A) The Δ (hunger-satiety) of fALFF in Fusiform connected with Δ glucose plasma after glucose administration. The statistical image was assessed for cluster-wise significance using a cluster defining threshold P < 0.001, 0.05 false discovery rate (FDR) corrected at the cluster level. (B) Amplitude of BOLD signal in the fusiform gyrus and plasma glucose levels in the hunger and satiety states. (C) The LOOCV showed that a positive correlation (r = 0.75) between the average Δ fALFF value of fusiform cluster and Δ plasma glucose, adjusted for Δ plasma insulin level. LOOCV, leave-one-out cross-validation.

Next, when controlling for Δ plasma glucose we found a connection between Δ fALFF and Δ plasma insulin after glucose administration in the left posterior insula and SFG (Figures 7A,B). A significant negative correlation (r = −0.7) was found between average Δ fALFF of the posterior Insula cluster and Δ plasma insulin (Figure 7C; red dots and line), as well as a significant negative correlation (r = −0.8) between average Δ fALFF SFG and Δ plasma insulin (Figure 7C; blue dots and line).

FIGURE 7

Figure 7. Insulin-associated low-frequency BOLD fluctuations. (A) The Δ (hunger-satiety) of fALFF in the left posterior (pos.) insula and superior frontal gyrus (SFG) connected with Δ plasma insulin after glucose administration. The results were assessed for cluster-wise significance using a cluster defining threshold P < 0.001, 0.05 FDR corrected at the cluster level. (B) SFG and pos. insula amplitude values of the BOLD signal and plasma insulin levels in the hunger and satiety states. (C) Scatter plot of the correlation between Δ plasma insulin and average Δ fALFF value of left pos. Insula (LOOCV of r = 0.7; red dots and line) and left SFG (LOOCV of r = 0.8; blue dots and line), adjusted for Δ plasma glucose level. Panels (D,E) represent a significant positive correlation between the average Δ fALFF values of left pos. insula and left SFG when the data was either not adjusted or adjusted for Δ plasma glucose, respectively. (F) The correlation between the average Δ fALFF value of left pos. insula and left SFG was no longer significant when the data was adjusted for Δ plasma insulin level, which suggests the effects were driven by plasma insulin. LOOCV, leave-one-out cross-validation.

A potential influence of plasma insulin and plasma glucose on the association between Δ fALFF posterior insula and Δ fALFF SFG was tested by using the average Δ fALFF in these clusters and revealed a significant positive correlation (r = 0.5, p = 0.004, Figure 7D) which survived when adjusting for Δ plasma glucose(r = 0.5, p = 0.004, Figure 7E) but not when adjusting for Δ plasma insulin level (r = −0.06, p = 0.7, Figure 7F).

Discussion

In the present study, we investigated the interaction of hunger and satiety with plasma glucose and plasma insulin levels before and after glucose administration and explored this interaction’s impact on changes in resting state fMRI. We used fALFF, a measure reflecting the temporal fluctuation of the BOLD signal, to study neural activity and connectivity. Multiple linear regressions analyses (MLRA) with two covariates (glucose, insulin) was used to assess the moderating influence of these covariates on the differences between hunger and satiety conditions.

Expectedly, plasma glucose and insulin levels differed between hunger and satiety conditions and were also differentially influenced by glucose administration (Figures 1, 2). In the satiety condition, a major peak in insulin and glucose levels around 3 h before glucose ingestion is related to the breakfast consumed in this condition. Generally, before glucose administration, insulin and glucose levels are higher in the satiety compared to fasting condition, which is the expected result (Johnson et al., 2006). After ingestion of glucose, there is a massive rise of both insulin and glucose levels as expected. This curve is higher and wider for both parameters in the fasting condition, which can be explained by post-fasting glucose intolerance (Johnson et al., 2006). Furthermore, the delay in the time till the maximum is reached in the fasting condition compared to the satiety condition can be explained with a delayed gastric emptying time for fasted individuals (Corvilain et al., 1995). The analyses of the fALFF data revealed an increased BOLD signal amplitude in the PCC in hunger relative to satiety that was independent of glucose and insulin plasma levels. Brain activity reflecting the difference between hunger vs. satiety was modulated by changes in glucose and insulin plasma levels in fusiform and PINS after oral glucose intake, respectively.

Physiological and Behavioral Effects

The results for glucose (Figure 2) suggest that hunger is stimulated, at least in part, by changing the blood glucose level (Pannacciulli et al., 2007). Also, a decrease in insulin concentrations in fasting compared to the fed state might function as a trigger to increase food intake (i.e., carbohydrate intake) in lean men (Mars et al., 2005). In addition, a meta-analysis showed that high insulin levels in response to a meal suppress appetite (Flint et al., 2007). Therefore, glucose and insulin may act as biomarkers of the satiety (de Graaf et al., 2004). Hunger ratings have been shown to increase with blood glucose concentration of 4.4 mmol/L (85 mg/dL) or lower (Ciampolini and Bianchi, 2006). In our study, the mean blood glucose level was 3.6 mmol/L under hunger condition (Figure 1B), and the feeling of hunger was higher under hunger relative to satiety condition before the fMRI scans (Figure 2C). These findings suggest that the experimental manipulation worked.

Resting-State fMRI: Hunger vs. Satiety Effects

The amplitude of the rs-fMRI signal was increased for hunger vs. satiety in the left IPG (IPG) before oral glucose treatment and in the PCC before and after oral glucose treatment (Figure 4). Several studies using visual stimuli showed enhanced activity of IPG, PCC, precuneus and insula during food choice vs. non-food choice under hunger conditions (Führer et al., 2008; Charbonnier et al., 2015). Moreover, PCC and precuneus are core hubs of the DMN (Cavanna and Trimble, 2006; Benedict et al., 2008). The DMN is activated during the rs-fMRI paradigm and deactivated during task-based fMRI, when attention is oriented towards internal rather than external stimuli (Gerozissis, 2003; Buckner et al., 2008). In line with that, the DMN, especially in the PCC, has been shown to be involved in self-referential processes (Buckner and Carroll, 2007). A study using ICA approach to investigate the FC of obese/overweight subjects found a reduction of DMN activity associated with a decrease in hunger ratings and fat mass after 6 months of exercise (McFadden et al., 2013). In the current experiment, the increased fALFF in PCC region seems to be driven by metabolic state (hunger vs. satiety; Figure 4) independently of changes in the plasma glucose and insulin levels. Consistent with our results, DMN activity has been shown to be positively correlated with lower blood glucose levels (Bolo et al., 2015) and with hunger ratings (McFadden et al., 2013). Our findings in the current and previous study (Al-Zubaidi et al., 2018) about increased activation of the DMN (i.e., PCC) under hunger condition might be indicative of increased self-focus during food deprivation.

We found a significant increase of fALFF in the left mediodorsal thalamus under hunger conditions and a decrease in satiety after (but not before) glucose treatment indicating an interaction between metabolic state and glucose administration in this region. Thalamic brain activity has been reported to vary as a function of hunger (Tataranni et al., 1999; Zhao et al., 2016), ghrelin application (Higgins et al., 2007) and glucose infusion (Jones et al., 2012; Little et al., 2014). A previous study (Kroemer et al., 2013a) reported a positive correlation between fasting levels of ghrelin and thalamic reactivity to food images, though this effect was not further modulated by glucose intake (75 g glucose). The thalamus is central for the integration of proprioceptive information stemming from the gastrointestinal (GI) tract (Kelley et al., 2005; Little et al., 2014). Specifically, afferent information from this tract reaches the thalamus via the vagus nerve (Coss-Adame and Rao, 2014). Consistent with this observation, Rolls (2005, 2006) postulated that the mediodorsal thalamus impacts short-term eating behavior. Our results agree with these proposals in that mediodorsal thalamus activity is partly dependent on metabolic state and associated with glucose metabolism.

Resting State fMRI: Before vs. After Glucose Treatment

FALFF was lower in the left SMA after glucose ingestion (Figure 5). The SMA is crucial for planning movements (Cunnington et al., 2002; Nachev et al., 2008) and observation of actions (Grèzes and Decety, 2002). The response to glucose ingestion in the form of lower activity in the SMA could potentially lead to suppressed action planning or initiation because of an alteration in the available energy resources of the body. Therefore, plasma glucose might have an effect on brain regions controlling attention, food evaluation and motor planning. Further research is needed to back up this interpretation in the form of task-related fMRI or behavioral studies that show such functional differences in dependency of blood glucose levels.

Correlations Between Physiological and Neural Effects

We observed that changes in plasma insulin levels (hunger vs. satiety) after oral glucose administration were negatively associated with changes in the amplitude of the BOLD signal in the left PINS and the left SFG). Furthermore, we found a significant positive correlation between left PINS and left SFG when the ΔfALFF values were adjusted for the Δplasma glucose, while this relationship disappeared when the ΔfALFF values were adjusted for the Δplasma insulin (Figure 7). In contrast, changes in plasma glucose levels were positively correlated with changes in the fALFF in the left fusiform after oral glucose administration (Figure 6). The PINS is involved in sensory, motor and interoception taste intensity (Kurth et al., 2010; Kelly et al., 2012; Nieuwenhuys, 2012; Tang et al., 2012) but is also activated during hunger (Tataranni et al., 1999), during food craving (Siep et al., 2012), and on receiving an appetitive drink (Bohon and Stice, 2011). The SFG is activated in response to appetitive stimuli in fasted subjects (Malik et al., 2011; Martens et al., 2013), and it is frequently involved in inhibiting approach behavior towards food (Gautier et al., 2000; McCaffery et al., 2009; Batterink et al., 2010). The fusiform gyrus harbors high-level visual processes (e.g., face recognition; Hadjikhani et al., 2004), but has been shown to be activated in response of food cues after glucose administration (Heni et al., 2014), to high vs. low caloric food pictures in lean subjects (Kroemer et al., 2013b), and as a function of increasing hunger (Zhao et al., 2016).

Our data showed that oral glucose administration modulates the insulin-dependent association between PINS and SFG, as well as the association between left fusiform gyrus activity with plasma glucose. Our results suggest that the fusiform, PINS, and SFG activity is modulated by an intricate interplay of glucose and insulin levels, most likely to regulate ingestive behavior by differential engagement of attentional, executive and effective processes.

These findings dovetail nicely with results from rs-fMRI studies demonstrating that body weight and insulin levels influence the FC within and between cognitive and homeostatic brain regions (Kullmann et al., 2012, 2013).

Limitations

Although our study was tightly controlled, there are some limiting factors that should be taken into account. First, although the hypothalamus is a central region for the regulation of energy homeostasis (Reno et al., 2016; Rosario et al., 2016), it was not modulated in the current study, which may be due to insufficient resolution of this small structure in fMRI imaging (De Silva et al., 2012) and the use of cluster-level correction (number of voxels per cluster; Woo et al., 2014). Specific scanning protocols targeted at subcortical structures in combination with region-of-interest analyses might be needed to detect hypothalamic effects. Second, we only investigated healthy young male participants. Therefore, our results may not be readily generalizable to other populations. Second, we only investigated healthy young male participants. Therefore, our results may not be readily generalizable to other populations. For instance, Haase et al. (2011) reported gender differences in response to sweet taste under hunger and satiety conditions in several brain regions. Also, the effects of insulin signaling on the human brain show sex differences (Benedict et al., 2008; Ghasemi et al., 2013). There have been multiple previous studies addressing differences in brain activity between lean and obese subjects (Hogenkamp et al., 2016). It is thus conceivable that the rs-fMRI measures used in the current study might also differ between lean and obese subjects. Furthermore, the correlation analysis after glucose administration was only performed with endocrine data narrowly associated with the fMRI data collection to keep this connection as clear as possible. Finally, our study focused on revealing correlations between plasma insulin and glucose concentrations and brain activity under different metabolic states. To extend these results, future studies might employ a longitudinal design with multiple fMRI sessions during the course of the fasting and control conditions and collection of the endocrine data. If the shown correlations persist or build up in a sensible way over time these results would be a strong addition to the correlational findings of this study. As this study was designed to confirm previous findings and to show the validity of these potential correlations in a physiological setting with an oral glucose intake intervention we only performed two fMRI sessions. To establish a causal relationship in contrast to our correlational findings, an experimental manipulation using insulin clamp and/or glucose clamp techniques is needed.

Conclusion

Our results suggest that plasma glucose and insulin respond similarly to oral glucose intake depending on metabolic states (hunger vs. satiety) and that these responses are related to different neural processing in the brain. Changes in plasma glucose were associated with changes of activation patterns in the fusiform gyrus, while changes in plasma insulin enhanced connectivity between the posterior insula and SFG when added as a covariate in the MLRA, indicating that changes in plasma insulin levels were at least partially responsible for the augmented connectivity. This connectivity appears to be related to alleviating an acute hemostatic energy deficit. These results contribute to identifying the neural mechanisms through which insulin regulates food intake (Kullmann et al., 2017). All in all, our findings expand existing neural models of homeostatic energy and highlight the complex nature of food intake and hormone-relationships in humans.

Ethics Statement

The study was approved by the ethics committee of the University of Lübeck and carried out under the Declaration of Helsinki (2000). Before participation, each participant gave written informed consent.

Author Contributions

KJ-C and TM designed the study and wrote the protocol. AA-Z, MH, JN and KJ-C participated in the data collection. AA-Z managed the literature searches and wrote the first draft of the manuscript. AA-Z and AM performed data processing and statistical analyses. MH, KJ-C, GB and TM helped with interpretation of data. All authors contributed to and had approved the final manuscript.

Funding

This research was supported by a grant of the German Research Foundation to the Research Training Group 1957 “Adipocyte-Brain Crosstalk” and the SFB TR 134, project C1. We acknowledge financial support by Land Schleswig-Holstein within the funding programme Open Access Publication Funds.

Conflict of Interest Statement

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.

Acknowledgments

AA-Z thanks Macià Buades Rotger for helpful comments on this work.

Abbreviations

APCUNS, anterior precuneus; AUC, area under the concentration-time curve; BMI, body mass index; BOLD, Blood-oxygen-level-dependent; CSF, cerebrospinal fluid; DARTEL, diffeomorphic anatomical registration through exponentiated Liealgebra; CBF, cerebral blood flow; CNS, central nervous system; Δ, delta = hunger-satiety; DMN, default mode network; DPARSFA, data processing assistant for resting-state fMRI advanced edition; EPI, echo-planar imaging; fALFF, fractional amplitude of low-frequency fluctuations; FC, functional connectivity; fMRI, functional magnetic resonance imaging; FSL, FMRIB Software Library; FDR, false discovery rate; ICA-AROMA, ICA-based strategy for automatic removal of motion artifacts; IPG, inferior parietal gyrus; K, number of voxel per cluster; LOOCV, leave-one-out cross validation; M, mean; MLRA, multiple linear regression analysis; MNI, Montreal Neurological Institute; PCC, posterior cingulate cortex; PET, positron emission tomography; PINS, posterior insula; rs-fMRI, resting-state functional magnetic resonance imaging; rm-ANOVA, repeated measures analysis of variance; SD, standard deviation; SPM, statistical parametric mapping; SFG, superior frontal gyrus; TE, echo time; TR, repetition time.

Footnotes

1. ^ http://www.alivelearn.net/xjview8

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnhum.2019.00162/full#supplementary-material

References

Al-Zubaidi, A., Heldmann, M., Mertins, A., Jauch-Chara, K., and Münte, T. F. (2018). Influences of hunger, satiety and oral glucose on functional brain connectivity: a multimethod resting-state fMRI study. Neuroscience 382, 80–92. doi: 10.1016/j.neuroscience.2018.04.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashburner, J., and Friston, K. J. (2005). Unified segmentation. Neuroimage 26, 839–851. doi: 10.1016/j.neuroimage.2005.02.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Batterink, L., Yokum, S., and Stice, E. (2010). Body mass correlates inversely with inhibitory control in response to food among adolescent girls: an fMRI study. Neuroimage 52, 1696–1703. doi: 10.1016/j.neuroimage.2010.05.059

PubMed Abstract | CrossRef Full Text | Google Scholar

Benedict, C., Kern, W., Schultes, B., Born, J., and Hallschmid, M. (2008). Differential sensitivity of men and women to anorexigenic and memory-improving effects of intranasal insulin. J. Clin. Endocrinol. Metab. 93, 1339–1344. doi: 10.1210/jc.2007-2606

PubMed Abstract | CrossRef Full Text | Google Scholar

Biswal, B., Yetkin, F. Z., Haughton, V. M., and Hyde, J. S. (1995). Functional connectivity in the motor cortex of resting human brain using. Magn. Reson. Med. 34, 537–541. doi: 10.1002/mrm.1910340409

PubMed Abstract | CrossRef Full Text | Google Scholar

Bohon, C., and Stice, E. (2011). Reward abnormalities among women with full and subthreshold bulimia nervosa: a functional magnetic resonance imaging study. Int. J. Eat. Disord. 44, 585–595. doi: 10.1002/eat.20869

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolo, N. R., Musen, G., Simonson, D. C., Nickerson, L. D., Flores, V. L., Siracusa, T., et al. (2015). Functional connectivity of insula, basal ganglia and prefrontal executive control networks during hypoglycemia in type 1 diabetes. J. Neurosci. 35, 11012–11023. doi: 10.1523/JNEUROSCI.0319-15.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Brown, L. M., Clegg, D. J., Benoit, S. C., and Woods, S. C. (2006). Intraventricular insulin and leptin reduce food intake and body weight in C57BL/6J mice. Physiol. Behav. 89, 687–691. doi: 10.1016/j.physbeh.2006.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Brüning, J. C., Gautam, D., Burks, D. J., Gillette, J., Schubert, M., Orban, P. C., et al. (2000). Role of brain insulin receptor in control of body weight and reproduction. Science 289, 2122–2125. doi: 10.1126/science.289.5487.2122

PubMed Abstract | CrossRef Full Text | Google Scholar

Buckner, R. L., Andrews-Hanna, J. R., and Schacter, D. L. (2008). The brain’s default network: anatomy, function, and relevance to disease. Ann. N Y Acad. Sci. 1124, 1–38. doi: 10.1196/annals.1440.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Buckner, R. L., and Carroll, D. C. (2007). Self-projection and the brain. Trends Cogn. Sci. 11, 49–57. doi: 10.1016/j.tics.2006.11.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Cavanna, A. E., and Trimble, M. R. (2006). The precuneus: a review of its functional anatomy and behavioural correlates. Brain 129, 564–583. doi: 10.1093/brain/awl004

PubMed Abstract | CrossRef Full Text | Google Scholar

Charbonnier, L., van der Laan, L. N., Viergever, M. A., Smeets, P. A. M., and Tregellas, J. R. (2015). Functional MRI of challenging food choices: forced choice between equally liked high-and low-calorie foods in the absence of hunger. PLoS One 10:e0131727. doi: 10.1371/journal.pone.0131727

PubMed Abstract | CrossRef Full Text | Google Scholar

Chavez, M., Kaiyala, K., Madden, L. J., Schwartz, M. W., and Woods, S. C. (1995). Intraventricular insulin and the level of maintained body weight in rats. Behav. Neurosci. 109, 528–531. doi: 10.1037/0735-7044.109.3.528

PubMed Abstract | CrossRef Full Text | Google Scholar

Chumbley, J., Worsley, K., Flandin, G., and Friston, K. (2010). Topological FDR for neuroimaging. Neuroimage 49, 3057–3064. doi: 10.1016/j.neuroimage.2009.10.090

PubMed Abstract | CrossRef Full Text | Google Scholar

Ciampolini, M., and Bianchi, R. (2006). Training to estimate blood glucose and to form associations with initial hunger. Nutr. Metab. 3:42. doi: 10.1186/1743-7075-3-42

PubMed Abstract | CrossRef Full Text | Google Scholar

Cole, M. W., Pathak, S., and Schneider, W. (2010). Identifying the brain’s most globally connected regions. Neuroimage 49, 3132–3148. doi: 10.1016/j.neuroimage.2009.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Corvilain, B., Abramowicz, M., Fery, F., Schoutens, A., Verlinden, M., Balasse, E., et al. (1995). Effect of short-term starvation on gastric emptying in humans: relationship to oral glucose tolerance. Am. J. Physiol. Liver Physiol. 269, G512–G517. doi: 10.1152/ajpgi.1995.269.4.g512

PubMed Abstract | CrossRef Full Text | Google Scholar

Coss-Adame, E., and Rao, S. S. C. (2014). Brain and gut interactions in irritable bowel syndrome: new paradigms and new understandings. Curr. Gastroenterol. Rep. 16:379. doi: 10.1007/s11894-014-0379-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Cunnington, R., Windischberger, C., Deecke, L., and Moser, E. (2002). The preparation and execution of self-initiated and externally-triggered movement: a study of event-related fMRI. Neuroimage 15, 373–385. doi: 10.1006/nimg.2001.0976

PubMed Abstract | CrossRef Full Text | Google Scholar

de Graaf, C., Blom, W. A. M., Smeets, P. A. M., Stafleu, A., and Hendriks, H. F. J. (2004). Biomarkers of satiation and satiety. Am. J. Clin. Nutr. 79, 946–961. doi: 10.1093/ajcn/79.6.946

PubMed Abstract | CrossRef Full Text | Google Scholar

De Silva, A., Salem, V., Matthews, P. M., and Dhillo, W. S. (2012). The use of functional MRI to study appetite control in the CNS. Exp. Diabetes Res. 2012:764017. doi: 10.1155/2012/764017

PubMed Abstract | CrossRef Full Text | Google Scholar

Egorova, N., Veldsman, M., Cumming, T., and Brodtmann, A. (2017). Fractional amplitude of low-frequency fluctuations (fALFF) in post-stroke depression. Neuroimage Clin. 16, 116–124. doi: 10.1016/j.nicl.2017.07.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Esterman, M., Tamber-Rosenau, B. J., Chiu, Y.-C., and Yantis, S. (2010). Avoiding non-independence in fMRI data analysis: leave one subject out. Neuroimage 50, 572–576. doi: 10.1016/j.neuroimage.2009.10.092

PubMed Abstract | CrossRef Full Text | Google Scholar

Field, A. (2014). Discovering Statistics Using IBM SPSS Statistics. 4th Edn. London: Sage Publications Ltd.

Google Scholar

Figlewicz, D. P. (2003). Adiposity signals and food reward: expanding the CNS roles of insulin and leptin. Am. J. Physiol. Regul. Integr. Comp. Physiol. 284, R882–R892. doi: 10.1152/ajpregu.00602.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

Flint, A., Gregersen, N. T., Gluud, L. L., Møller, B. K., Raben, A., Tetens, I., et al. (2007). Associations between postprandial insulin and blood glucose responses, appetite sensations and energy intake in normal weight and overweight individuals: a meta-analysis of test meal studies. Br. J. Nutr. 98, 17–25. doi: 10.1017/s000711450768297x

PubMed Abstract | CrossRef Full Text | Google Scholar

Führer, D., Zysset, S., and Stumvoll, M. (2008). Brain activity in hunger and satiety: an exploratory visually stimulated FMRI study. Obesity 16, 945–950. doi: 10.1038/oby.2008.33

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, Q., and Horvath, T. L. (2008). Neuronal control of energy homeostasis. FEBS Lett. 582, 132–141. doi: 10.1016/j.febslet.2007.11.063

PubMed Abstract | CrossRef Full Text | Google Scholar

Gautier, J. F., Chen, K., Salbe, A. D., Bandy, D., Pratley, R. E., Heiman, M., et al. (2000). Differential brain responses to satiation in obese and lean men. Diabetes 49, 838–846. doi: 10.2337/diabetes.49.5.838

PubMed Abstract | CrossRef Full Text | Google Scholar

Gerozissis, K. (2003). Brain insulin: regulation, mechanisms of action and functions. Cell. Mol. Neurobiol. 23, 1–25. doi: 10.1023/A:1022598900246

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghasemi, R., Haeri, A., Dargahi, L., Mohamed, Z., and Ahmadiani, A. (2013). Insulin in the brain: sources, localization and functions. Mol. Neurobiol. 47, 145–171. doi: 10.1007/s12035-012-8339-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Grèzes, J., and Decety, J. (2002). Does visual perception of object afford action? Evidence from a neuroimaging study. Neuropsychologia 40, 212–222. doi: 10.1016/s0028-3932(01)00089-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Haase, L., Green, E., and Murphy, C. (2011). Males and females show differential brain activation to taste when hungry and sated in gustatory and reward areas. Appetite 57, 421–434. doi: 10.1016/j.appet.2011.06.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Hadjikhani, N., Joseph, R. M., Snyder, J., Chabris, C. F., Clark, J., Steele, S., et al. (2004). Activation of the fusiform gyrus when individuals with autism spectrum disorder view faces. Neuroimage 22, 1141–1150. doi: 10.1016/j.neuroimage.2004.03.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Heni, M., Kullmann, S., Ketterer, C., Guthoff, M., Bayer, M., Staiger, H., et al. (2014). Differential effect of glucose ingestion on the neural processing of food stimuli in lean and overweight adults. Hum. Brain Mapp. 35, 918–928. doi: 10.1002/hbm.22223

PubMed Abstract | CrossRef Full Text | Google Scholar

Higgins, S. C., Gueorguiev, M., and Korbonits, M. (2007). Ghrelin, the peripheral hunger hormone. Ann. Med. 39, 116–136. doi: 10.1080/07853890601149179

PubMed Abstract | CrossRef Full Text | Google Scholar

Hogenkamp, P. S., Zhou, W., Dahlberg, L. S., Stark, J., Larsen, A. L., Olivo, G., et al. (2016). Higher resting-state activity in reward-related brain circuits in obese versus normal-weight females independent of food intake. Int. J. Obes. 40, 1687–1692. doi: 10.1038/ijo.2016.105

PubMed Abstract | CrossRef Full Text | Google Scholar

Jakab, A., Blanc, R., Berényi, E. L., and Székely, G. (2012). Generation of individualized thalamus target maps by using statistical shape models and thalamocortical tractography. AJNR Am. J. Neuroradiol. 33, 2110–2116. doi: 10.3174/ajnr.a3140

PubMed Abstract | CrossRef Full Text | Google Scholar

Jenkinson, M., Bannister, P., Brady, M., and Smith, S. (2002). Improved optimization for the robust and accurate linear registration and motion correction of brain images. Neuroimage 17, 825–841. doi: 10.1016/s1053-8119(02)91132-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, N. A., Stannard, S. R., Rowlands, D. S., Chapman, P. G., Thompson, C. H., O’Connor, H., et al. (2006). Effect of short-term starvation versus high-fat diet on intramyocellular triglyceride accumulation and insulin resistance in physically fit men. Exp. Physiol. 91, 693–703. doi: 10.1113/expphysiol.2006.033399

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, R. B., McKie, S., Astbury, N., Little, T. J., Tivey, S., Lassman, D. J., et al. (2012). Functional neuroimaging demonstrates that ghrelin inhibits the central nervous system response to ingested lipid. Gut 61, 1543–1551. doi: 10.1136/gutjnl-2011-301323

PubMed Abstract | CrossRef Full Text | Google Scholar

Kelley, A. E., Baldo, B. A., Pratt, W. E., and Will, M. J. (2005). Corticostriatal-hypothalamic circuitry and food motivation: integration of energy, action and reward. Physiol. Behav. 86, 773–795. doi: 10.1016/j.physbeh.2005.08.066

PubMed Abstract | CrossRef Full Text | Google Scholar

Kelly, C., Toro, R., Di Martino, A., Cox, C. L., Bellec, P., Castellanos, F. X., et al. (2012). A convergent functional architecture of the insula emerges across imaging modalities. Neuroimage 61, 1129–1142. doi: 10.1016/j.neuroimage.2012.03.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Kiebel, S. J., and Holmes, A. P. (2007). “The general linear modal,” in Statistical Parametric Mapping: The Analysis of Functional Brain Images, eds K. J. Friston, J. Ashburner, S. J. Kiebel, T. E. Nichols and W. D. Penny (London: Academic Press), 101–125.

Klein, A., Andersson, J., Ardekani, B. A., Ashburner, J., Avants, B., Chiang, M.-C., et al. (2009). Evaluation of 14 nonlinear deformation algorithms applied to human brain MRI registration. Neuroimage 46, 786–802. doi: 10.1016/j.neuroimage.2008.12.037

PubMed Abstract | CrossRef Full Text | Google Scholar

Kroemer, N. B., Krebs, L., Kobiella, A., Grimm, O., Pilhatsch, M., Bidlingmaier, M., et al. (2013a). Fasting levels of ghrelin covary with the brain response to food pictures. Addict. Biol. 18, 855–862. doi: 10.1111/j.1369-1600.2012.00489.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kroemer, N. B., Krebs, L., Kobiella, A., Grimm, O., Vollstädt-Klein, S., Wolfensteller, U., et al. (2013b). (Still) longing for food: insulin reactivity modulates response to food pictures. Hum. Brain Mapp. 34, 2367–2380. doi: 10.1002/hbm.22071

PubMed Abstract | CrossRef Full Text | Google Scholar

Kullmann, S., Frank, S., Heni, M., Ketterer, C., Veit, R., Häring, H. U., et al. (2013). Intranasal insulin modulates intrinsic reward and prefrontal circuitry of the human brain in lean women. Neuroendocrinology 97, 176–182. doi: 10.1159/000341406

PubMed Abstract | CrossRef Full Text | Google Scholar

Kullmann, S., Heni, M., Veit, R., Ketterer, C., Schick, F., Häring, H. U., et al. (2012). The obese brain: association of body mass index and insulin sensitivity with resting state network functional connectivity. Hum. Brain Mapp. 33, 1052–1061. doi: 10.1002/hbm.21268

PubMed Abstract | CrossRef Full Text | Google Scholar

Kullmann, S., Heni, M., Veit, R., Scheffler, K., Machann, J., Häring, H.-U., et al. (2017). Intranasal insulin enhances brain functional connectivity mediating the relationship between adiposity and subjective feeling of hunger. Sci. Rep. 7:1627. doi: 10.1038/s41598-017-01907-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Kurth, F., Zilles, K., Fox, P. T., Laird, A. R., and Eickhoff, S. B. (2010). A link between the systems: functional differentiation and integration within the human insula revealed by meta-analysis. Brain Struct. Funct. 214, 519–534. doi: 10.1007/s00429-010-0255-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Lam, T. K. T. (2010). Neuronal regulation of homeostasis by nutrient sensing. Nat. Med. 16, 392–395. doi: 10.1038/nm0410-392

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, M. H., Smyser, C. D., and Shimony, J. S. (2013). Resting-state fMRI: a review of methods and clinical applications. AJNR Am. J. Neuroradiol. 34, 1866–1872. doi: 10.3174/ajnr.A3263

PubMed Abstract | CrossRef Full Text | Google Scholar

Little, T. J., McKie, S., Jones, R. B., D’Amato, M., Smith, C., Kiss, O., et al. (2014). Mapping glucose-mediated gut-to-brain signalling pathways in humans. Neuroimage 96, 1–11. doi: 10.1016/j.neuroimage.2014.03.059

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, T. T. (2016). Noise contributions to the fMRI signal: an overview. Neuroimage 143, 141–151. doi: 10.1016/j.neuroimage.2016.09.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y., Gao, J. H., Liu, H. L., and Fox, P. T. (2000). The temporal response of the brain after eating revealed by functional MRI. Nature 405, 1058–1062. doi: 10.1038/35016590

PubMed Abstract | CrossRef Full Text | Google Scholar

Malik, S., McGlone, F., and Dagher, A. (2011). State of expectancy modulates the neural response to visual food stimuli in humans. Appetite 56, 302–309. doi: 10.1016/j.appet.2011.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Mars, M., de Graaf, C., de Groot, L. C., and Kok, F. J. (2005). Decreases in fasting leptin and insulin concentrations after acute energy restriction and subsequent compensation in food intake. Am. J. Clin. Nutr. 81, 570–577. doi: 10.1093/ajcn/81.3.570

PubMed Abstract | CrossRef Full Text | Google Scholar

Martens, M. J. I., Born, J. M., Lemmens, S. G. T., Karhunen, L., Heinecke, A., Goebel, R., et al. (2013). Increased sensitivity to food cues in the fasted state and decreased inhibitory control in the satiated state in the overweight. Am. J. Clin. Nutr. 97, 471–479. doi: 10.3945/ajcn.112.044024

PubMed Abstract | CrossRef Full Text | Google Scholar

McCaffery, J. M., Haley, A. P., Sweet, L. H., Phelan, S., Raynor, H. A., Del Parigi, A., et al. (2009). Differential functional magnetic resonance imaging response to food pictures in successful weight-loss maintainers relative to normal-weight and obese controls. Am. J. Clin. Nutr. 90, 928–934. doi: 10.3945/ajcn.2009.27924

PubMed Abstract | CrossRef Full Text | Google Scholar

McFadden, K. L., Cornier, M.-A., Melanson, E. L., Bechtell, J. L., and Tregellas, J. R. (2013). Effects of exercise on resting-state default mode and salience network activity in overweight/obese adults. Neuroreport 24, 866–871. doi: 10.1097/WNR.0000000000000013

PubMed Abstract | CrossRef Full Text | Google Scholar

McGowan, M. K., Andrews, K. M., Kelly, J., and Grossman, S. P. (1990). Effects of chronic intrahypothalamic infusion of insulin on food intake and diurnal meal patterning in the rat. Behav. Neurosci. 104, 373–385. doi: 10.1037//0735-7044.104.2.373

PubMed Abstract | CrossRef Full Text | Google Scholar

McLaren, D. G., Kosmatka, K. J., Kastman, E. K., Bendlin, B. B., and Johnson, S. C. (2010). Rhesus macaque brain morphometry: a methodological comparison of voxel-wise approaches. Methods 50, 157–165. doi: 10.1016/j.ymeth.2009.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Nachev, P., Kennard, C., and Husain, M. (2008). Functional role of the supplementary and pre-supplementary motor areas. Nat. Rev. Neurosci. 9, 856–869. doi: 10.1038/nrn2478

PubMed Abstract | CrossRef Full Text | Google Scholar

Nathans, L. L., Oswald, F. L., and Nimon, K. (2012). Interpreting multiple linear regression: a guidebook of variable importance. Pract. Assess. Res. Eval. 17, 1–19.

Google Scholar

Nieuwenhuys, R. (2012). The Insular Cortex: A Review. 1st Edn. Elsevier B.V.

Google Scholar

Obici, S., Zhang, B. B., Karkanias, G., and Rossetti, L. (2002). Hypothalamic insulin signaling is required for inhibition of glucose production. Nat. Med. 8, 1376–1382. doi: 10.1038/nm798

PubMed Abstract | CrossRef Full Text | Google Scholar

Owen, L., Scholey, A. B., Finnegan, Y., Hu, H., and Sünram-Lea, S. I. (2012). The effect of glucose dose and fasting interval on cognitive function: a double-blind, placebo-controlled, six-way crossover study. Psychopharmacology 220, 577–589. doi: 10.1007/s00213-011-2510-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Page, K. A., Chan, O., Arora, J., Belfort-Deaguiar, R., Dzuira, J., Roehmholdt, B., et al. (2013). Effects of fructose vs. glucose on regional cerebral blood flow in brain regions involved with appetite and reward pathways. JAMA 309, 63–70. doi: 10.1001/jama.2012.116975

PubMed Abstract | CrossRef Full Text | Google Scholar

Pannacciulli, N., Ortega, E., Koska, J., Salbe, A. D., Bunt, J. C., and Krakoff, J. (2007). Glucose response to an oral glucose tolerance test predicts weight change in non-diabetic subjects. Obesity 15, 632–639. doi: 10.1038/oby.2007.572

PubMed Abstract | CrossRef Full Text | Google Scholar

Pruessner, J. C., Kirschbaum, C., Meinlschmid, G., and Hellhammer, D. H. (2003). Two formulas for computation of the area under the curve represent measures of total hormone concentration versus time-dependent change. Psychoneuroendocrinology 28, 916–931. doi: 10.1016/s0306-4530(02)00108-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Pruim, R. H. R., Mennes, M., Buitelaar, J. K., and Beckmann, C. F. (2015a). Evaluation of ICA-AROMA and alternative strategies for motion artifact removal in resting state fMRI. Neuroimage 112, 278–287. doi: 10.1016/j.neuroimage.2015.02.063

PubMed Abstract | CrossRef Full Text | Google Scholar

Pruim, R. H. R., Mennes, M., van Rooij, D., Llera, A., Buitelaar, J. K., and Beckmann, C. F. (2015b). ICA-AROMA: a robust ICA-based strategy for removing motion artifacts from fMRI data. Neuroimage 112, 267–277. doi: 10.1016/j.neuroimage.2015.02.064

PubMed Abstract | CrossRef Full Text | Google Scholar

Reno, C. M., Puente, E. C., Sheng, Z., Daphna-Iken, D., Bree, A. J., Routh, V. H., et al. (2016). Brain GLUT4 knockout mice have impaired glucose tolerance, decreased insulin sensitivity, and impaired hypoglycemic counterregulation. Diabetes 66, 587–597. doi: 10.2337/db16-0917

PubMed Abstract | CrossRef Full Text | Google Scholar

Rolls, E. T. (2005). Taste, olfactory, and food texture processing in the brain and the control of food intake. Physiol. Behav. 85, 45–56. doi: 10.1016/j.physbeh.2005.04.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Rolls, E. T. (2006). Brain mechanisms underlying flavour and appetite. Philos. Trans. R. Soc. Lond. B Biol. Sci. 361, 1123–1136. doi: 10.1098/rstb.2006.1852

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosario, W., Singh, I., Wautlet, A., Patterson, C., Flak, J., Becker, T. C., et al. (2016). The brain to pancreatic islet neuronal map reveals differential glucose regulation from distinct hypothalamic regions. Diabetes 65, 2711–2723. doi: 10.2337/db15-0629

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwartz, M. W. (2001). Brain pathways controlling food intake and body weight. Exp. Biol. Med. 226, 978–981. doi: 10.1177/153537020122601103

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwartz, M. W., and Porte, D. Jr. (2005). Diabetes, obesity, and the brain. Science 307, 375–379. doi: 10.1126/science.1104344

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwartz, M. W., Woods, S. C., Porte, D. Jr., Seeley, R. J., and Baskin, D. G. (2000). Central nervous system control of food intake. Nature 404, 661–671. doi: 10.1038/35007534

PubMed Abstract | CrossRef Full Text | Google Scholar

Siep, N., Roefs, A., Roebroeck, A., Havermans, R., Bonte, M., and Jansen, A. (2012). Fighting food temptations: the modulating effects of short-term cognitive reappraisal, suppression and up-regulation on mesocorticolimbic activity related to appetitive motivation. Neuroimage 60, 213–220. doi: 10.1016/j.neuroimage.2011.12.067

PubMed Abstract | CrossRef Full Text | Google Scholar

Smitha, K. A., Akhil Raja, K., Arun, K. M., Rajesh, P. G., Thomas, B., Kapilamoorthy, T. R., et al. (2017). Resting state fMRI: a review on methods in resting state connectivity analysis and resting state networks. Neuroradiol. J. 30, 305–317. doi: 10.1177/1971400917697342

PubMed Abstract | CrossRef Full Text | Google Scholar

Suzuki, K., Simpson, K. A., Minnion, J. S., Shillito, J. C., and Bloom, S. R. (2010). The role of gut hormones and the hypothalamus in appetite regulation. Endocr. J. 57, 359–372. doi: 10.1507/endocrj.k10e-077

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, D. W., Fellows, L. K., Small, D. M., and Dagher, A. (2012). Food and drug cues activate similar brain regions: a meta-analysis of functional MRI studies. Physiol. Behav. 106, 317–324. doi: 10.1016/j.physbeh.2012.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Tataranni, P. A., Gautier, J.-F., Chen, K., Uecker, A., Bandy, D., Salbe, A. D., et al. (1999). Neuroanatomical correlates of hunger and satiation in humans using positron emission tomography. Proc. Natl. Acad. Sci. U S A 96, 4569–4574. doi: 10.1073/pnas.96.8.4569

PubMed Abstract | CrossRef Full Text | Google Scholar

Tomasi, D., Wang, G., and Volkow, N. (2013). Energetic cost of brain functional connectivity. Proc. Natl. Acad. Sci. U S A 110, 13642–13647. doi: 10.1073/pnas.1303346110

PubMed Abstract | CrossRef Full Text | Google Scholar

Tzagournis, M., and Skillman, T. G. (1970). Glucose intolerance mechanism after starvation. Metab. Clin. Exp. 19, 170–178. doi: 10.1016/0026-0495(70)90109-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Tzourio-Mazoyer, N., Landeau, B., Papathanassiou, D., Crivello, F., Etard, O., Delcroix, N., et al. (2002). Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage 15, 273–289. doi: 10.1006/nimg.2001.0978

PubMed Abstract | CrossRef Full Text | Google Scholar

van den Heuvel, M. P., and Hulshoff Pol, H. E. (2010). Exploring the brain network: a review on resting-state fMRI functional connectivity. Eur. Neuropsychopharmacol. 20, 519–534. doi: 10.1016/j.euroneuro.2010.03.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Wallner-Liebmann, S., Koschutnig, K., Reishofer, G., Sorantin, E., Blaschitz, B., Kruschitz, R., et al. (2010). Insulin and hippocampus activation in response to images of high-calorie food in normal weight and obese adolescents. Obesity 18, 1552–1557. doi: 10.1038/oby.2010.26

PubMed Abstract | CrossRef Full Text | Google Scholar

Williams, G., Harrold, J. A., and Cutler, D. J. (2000). The hypothalamus and the regulation of energy homeostasis: lifting the lid on a black box. Proc. Nutr. Soc. 59, 385–396. doi: 10.1017/s0029665100000434

PubMed Abstract | CrossRef Full Text | Google Scholar

Wölnerhanssen, B. K., Meyer-Gerspach, A. C., Schmidt, A., Zimak, N., Peterli, R., Beglinger, C., et al. (2015). Dissociable behavioral, physiological and neural effects of acute glucose and fructose ingestion: a pilot study. PLoS One 10:e0130280. doi: 10.1371/journal.pone.0130280

PubMed Abstract | CrossRef Full Text | Google Scholar

Woo, C. W., Krishnan, A., and Wager, T. D. (2014). Cluster-extent based thresholding in fMRI analyses: pitfalls and recommendations. Neuroimage 91, 412–419. doi: 10.1016/j.neuroimage.2013.12.058

PubMed Abstract | CrossRef Full Text | Google Scholar

Woods, S. C., Seeley, R. J., Porte, D. Jr., and Schwartz, M. W. (1998). Signals that regulate food intake and energy homeostasis. Science 280, 1378–1383. doi: 10.1126/science.280.5368.1378

PubMed Abstract | CrossRef Full Text | Google Scholar

Wright, H., Li, X., Fallon, N. B., Crookall, R., Giesbrecht, T., Thomas, A., et al. (2016). Differential effects of hunger and satiety on insular cortex and hypothalamic functional connectivity. Eur. J. Neurosci. 43, 1181–1189. doi: 10.1111/ejn.13182

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, C., and Zang, Y. (2010). DPARSF: a MATLAB toolbox for “pipeline” data analysis of resting-state fMRI. Front. Syst. Neurosci. 4:13. doi: 10.3389/fnsys.2010.00013

PubMed Abstract | CrossRef Full Text | Google Scholar

Zanchi, D., Depoorter, A., Egloff, L., Haller, S., Mählmann, L., Lang, U. E., et al. (2017). The impact of gut hormones on the neural circuit of appetite and satiety: a systematic review. Neurosci. Biobehav. Rev. 80, 457–475. doi: 10.1016/j.neubiorev.2017.06.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, J., Li, M., Zhang, Y., Song, H., von Deneen, K. M., Shi, Y., et al. (2016). Intrinsic brain subsystem associated with dietary restraint, disinhibition and hunger: an fMRI study. Brain Imaging Behav. 11, 264–277. doi: 10.1007/s11682-015-9491-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Zou, Q. H., Zhu, C. Z., Yang, Y., Zuo, X. N., Long, X. Y., Cao, Q. J., et al. (2008). An improved approach to detection of amplitude of low-frequency fluctuation (ALFF) for resting-state fMRI: fractional ALFF. J. Neurosci. Methods 172, 137–141. doi: 10.1016/j.jneumeth.2008.04.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Zuo, X. N., Di Martino, A., Kelly, C., Shehzad, Z. E., Gee, D. G., Klein, D. F., et al. (2010). The oscillating brain: complex and reliable. Neuroimage 49, 1432–1445. doi: 10.1016/j.neuroimage.2009.09.037

PubMed Abstract | CrossRef Full Text | Google Scholar