Data information
Twenty healthy volunteers (10 women and 10 men) with ages ranging from 24 to 40 (mean 32.1 ± 4.1) years were prospectively enrolled in this study. All subjects were right-handed without a history of any psychiatric, neurologic, or medical illness. Study procedures were approved by the Institutional Review Board of Chang Gung Medical Foundation, and informed written consent was obtained from each participant before enrollment.
To track how the human brain responds to water-related psychological changes, we designed an experiment that simultaneously collected MRI scans and diverse psychological data under different hydration statuses within 13 h in a longitudinal manner. As shown in Fig. 2a, our experiment set three fixed time points in the timeline for data collection, each representing one hydration status. At the first time point, the participants were in a normohydration status without any limitations of food or fluid intake in the past several hours. The second time point was set at 12 h after the first one, during which the participants were prohibited from eating or drinking anything with a water content for 12 h to achieve a dehydration status. The third time point represents the rehydration status, which was set to one hour after the participants were asked to drink 1.0 L water for rehydration.
The MRI data were collected with a 3-T Siemens Verio MRI system (Siemens Medical System, Erlangen, Germany) using a 32-channel head coil. All subjects were instructed to stay awake and relaxed with their eyes closed during the scanning. The MRI protocol was performed in all three hydration statuses, including baseline, dehydration, and rehydration. Rs-fMRI data were acquired with a gradient EPI sequence using the following parameters: TR = 2500 ms, TE = 27 ms, FOV = 220 mm, matrix = 64×64×36, voxel size = 3.4×3.4×4 mm3, and total volumes = 240 (10 min 7 sec). 3D MP-RAGE anatomical images were acquired using a gradient echo sequence with the following parameters: TR = 1900 ms, TE = 2.98 ms, FOV = 230 mm, matrix = 220×256, slice number = 160, and voxel size = 0.9×0.9×0.9 mm3.
The categorical rating thirsty scale was used to assess the sensation of thirst [39]. Briefly, all subjects were asked to answer the question of “How thirsty do you feel now?” by placing a mark on a 10-cm horizontal line anchored by the phrases “not at all” and “very thirsty” at bilateral extremes. The physiological status of hydration was evaluated with three serum parameters (B) [i.e., glucose, sodium (Na (B)) and osmolality (osmolality (B))] and four urine parameters (U) [i.e., sodium (Na (U)), osmolality (osmolality (U)), creatinine (creatinine (U)) and specific gravity (SP gravity)]. Body water composition, including TBW, ICW, and ECW, was measured by a bioelectrical impedance body composition analyzer (X-SCAN PLUS II, Jawon Medical., Korea).
Data preprocessing
The rs-fMRI data were preprocessed with AFNI [40]. Allowing for the equilibration of the magnetic field, the first 10 volumes were discarded for each scan of each subject. The remaining rs-fMRI data were preprocessed with the following steps: (1) slice timing; (2) head motion correction; (3) normalization to the MNI atlas with a resolution of 2×2×2 mm3; (4) spatial smoothing with a 6-mm full-width half-maximum (FWHM) Gaussian kernel; (5) mean signal removal (including white matter, cerebrospinal fluid (CSF), whole-brain averaged signal, and six head motion parameters) by using the linear regression model; and (6) temporal bandpass filtering (0.01–0.1 Hz). Additionally, to further remove head motion effects, we scrubbed the data and deleted one volume before and after each bad frame to control the global measure of framewise displacement (FD > 0.5 mm) [41, 42]. After scrubbing, subjects with less than 120 volumes were excluded, and one subject was removed from further analysis. Finally, the spike was removed from the blood oxygen level-dependent (BOLD) time series by using 3dDespike in AFNI.
All collected physiological indices of hydration were normalized to z scores and then pooled into the PCA model to extract the factors that primarily contribute to the overall hydration physiological status. In our work, the factors (i.e., eigenvectors) with eigenvalue scores larger than 1 were retained as PCs, and all physiological indices were projected on the PCs to generate the new physiological representation (i.e., PC scores) with low dimensionality. Notably, all of the following experiments on physiological indices were conducted on the PC scores instead of the original physiological data to reduce the analysis complexity.
Functional connectivity map construction
All rs-fMRI data (total 57 scans: 19 subjects × 3 scans) were used to generate resting-state FC network maps by using the gICA. The number of components was set to 15 based on previous studies [43, 44]. To investigate the individual-level and group-level FC networks simultaneously, dual regression was specifically adopted in our work to identify the subject-specific spatial distribution map of each available component [45], and a one-sample t test was then performed on all individual-level maps of each component to identify the corresponding group-level spatial distribution. Based on all detected group-level FC components, three well-known spontaneous-related networks were selected by visual inspection, which included the DMN, FPCN, and SN [17].
As a hub region for thirst regulation, the MnPO was also included as a region of interest (ROI) in our work to explore its FC alterations with three identified thought-related functional subnetworks under three hydration statuses. To this end, two spherical seeds with diameter = 3 mm were first manually placed in the bilateral MnPO (see Fig. 1) in the Montreal Neurological Institute coordinates of [±3.5, 0.6, −13.2] according to [46]. MnPO-based FC maps were then constructed by estimating the Pearson correlation of BOLD time courses between the MnPO seeds and voxels in the DMN, FPCN, and SN.
Comparisons of psychological state and brain functional connectivity under different hydration statuses
To explore whether psychological states and functional connectivity are altered under different hydration statuses, repeated-measures ANOVA was performed on (1) thirstiness scores, (2) psychological PC scores, (3) FCs within three thought-related functional subnetworks (i.e., DMN, FPCN, and SN), and (4) FCs between the three networks and bilateral MnPO. For thirst scores and physiological data, the significance levels were set to p < 0.05 after false discovery rate (FDR) correction. For FCs, the significance levels were set to p
corrected < 0.05 (uncorrected p = 0.01, cluster size = 600 mm3) with 3Dclustsim correction in AFNI.
Uncovering functional pathways for thirst regulation with a multiple mediation analysis model
In addition to the above cross-sectional comparisons of physiological state and FCs between different hydration statuses, we further investigated how the brain integrates diverse regions to respond to physiological changes. To this end, we first calculated the direct relationships between two physiological indices (i.e., DCI and DAI) and VAS. In our work, the LMER model [47] was used instead of the general linear regression model to estimate the correlation coefficient between two variables for its ability to handle longitudinal data. In each LMER model, age and gender were added as the covariates, and random intercept and subject effects were included as the characterization of the temporal correlation. Second, the serial multiple mediation analysis model was adopted to estimate the indirect relationship between the physiological state and thirsty scale with the MnPO and the three identified thought-related networks as mediations [23]. In this model (see Fig. 5), two physiological PC scores (i.e., DCI and DAI) were entered as independent variables X, and VAS was entered as dependent variable Y. The significantly changing FCs between the MnPO and the three intrinsic networks across the difference groups were entered as the first-level mediator variables M
1, and the significantly changing FCs within DMN, FPCN, and SN across the difference groups were entered as the second-level mediator variables M
2. Similarly, the LMER model was used to evaluate the correlation coefficient between the input and output variables. For each pathway (Fig. 5), we utilized the bootstrapping test 1000 times to test the significance. Finally, the mediation effect of each significant path was evaluated using the Sobel test [23]. The z value (Sobel test) is estimated:
$$z\ value=\frac{a_1\times {b}_2\times d}{\sqrt{{a_1}^2\times {b_2}^2\times SE(d)+{a_1}^2\times {d}^2\times SE\left({b}_2\right)+{b_2}^2\times {d}^2\times SE\left({a}_1\right)}},$$
where a
1 is the mean coefficient of the path from X to M
1, d is the mean coefficient of the path from M
1 to M
2, b
2 is the mean coefficient of the path from M
2 to Y, and SE represents the standard error.