
==== Front
Eur J Nucl Med Mol Imaging
Eur J Nucl Med Mol Imaging
European Journal of Nuclear Medicine and Molecular Imaging
1619-7070
1619-7089
Springer Berlin Heidelberg Berlin/Heidelberg

38052927
6542
10.1007/s00259-023-06542-4
Original Article
High-temporal resolution functional PET/MRI reveals coupling between human metabolic and hemodynamic brain response
http://orcid.org/0000-0001-9727-7580
Hahn Andreas andreas.hahn@meduniwien.ac.at

12
http://orcid.org/0000-0002-4873-608X
Reed Murray B. 12
http://orcid.org/0000-0003-2065-6093
Vraka Chrysoula 3
http://orcid.org/0000-0002-9739-0724
Godbersen Godber M. 12
http://orcid.org/0000-0001-8714-6608
Klug Sebastian 12
http://orcid.org/0000-0002-1674-4333
Komorowski Arkadiusz 12
http://orcid.org/0000-0003-4817-3499
Falb Pia 12
http://orcid.org/0000-0002-1034-876X
Nics Lukas 3
http://orcid.org/0000-0003-1118-926X
Traub-Weidinger Tatjana 3
http://orcid.org/0000-0002-4222-4083
Hacker Marcus 3
http://orcid.org/0000-0003-4641-9539
Lanzenberger Rupert rupert.lanzenberger@meduniwien.ac.at

12
1 https://ror.org/05n3x4p02 grid.22937.3d 0000 0000 9259 8492 Department of Psychiatry and Psychotherapy, Medical University of Vienna, Waehringer Guertel 18-20, 1090 Vienna, Austria
2 https://ror.org/05n3x4p02 grid.22937.3d 0000 0000 9259 8492 Comprehensive Center for Clinical Neurosciences and Mental Health (C3NMH), Medical University of Vienna, Vienna, Austria
3 https://ror.org/05n3x4p02 grid.22937.3d 0000 0000 9259 8492 Department of Biomedical Imaging and Image-Guided Therapy, Division of Nuclear Medicine, Medical University of Vienna, Vienna, Austria
6 12 2023
6 12 2023
2024
51 5 13101322
25 8 2023
22 11 2023
© The Author(s) 2023
2023
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Purpose

Positron emission tomography (PET) provides precise molecular information on physiological processes, but its low temporal resolution is a major obstacle. Consequently, we characterized the metabolic response of the human brain to working memory performance using an optimized functional PET (fPET) framework at a temporal resolution of 3 s.

Methods

Thirty-five healthy volunteers underwent fPET with [18F]FDG bolus plus constant infusion, 19 of those at a hybrid PET/MRI scanner. During the scan, an n-back working memory paradigm was completed. fPET data were reconstructed to 3 s temporal resolution and processed with a novel sliding window filter to increase signal to noise ratio. BOLD fMRI signals were acquired at 2 s.

Results

Consistent with simulated kinetic modeling, we observed a constant increase in the [18F]FDG signal during task execution, followed by a rapid return to baseline after stimulation ceased. These task-specific changes were robustly observed in brain regions involved in working memory processing. The simultaneous acquisition of BOLD fMRI revealed that the temporal coupling between hemodynamic and metabolic signals in the primary motor cortex was related to individual behavioral performance during working memory. Furthermore, task-induced BOLD deactivations in the posteromedial default mode network were accompanied by distinct temporal patterns in glucose metabolism, which were dependent on the metabolic demands of the corresponding task-positive networks.

Conclusions

In sum, the proposed approach enables the advancement from parallel to truly synchronized investigation of metabolic and hemodynamic responses during cognitive processing. This allows to capture unique information in the temporal domain, which is not accessible to conventional PET imaging.

Keywords

Simultaneous PET/MRI
Brain metabolism
Functional PET
BOLD signal
High temporal resolution
http://dx.doi.org/10.13039/501100002428 Austrian Science Fund KLI 610 KLI 1006 Hahn Andreas Lanzenberger Rupert http://dx.doi.org/10.13039/501100001821 Vienna Science and Technology Fund 10.47379/CS18039 Lanzenberger Rupert http://dx.doi.org/10.13039/501100001822 Österreichischen Akademie der Wissenschaften Austrian Science Fund (FWF)Open access funding provided by Austrian Science Fund (FWF).

issue-copyright-statement© Springer-Verlag GmbH Germany, part of Springer Nature 2024
==== Body
pmcIntroduction

The human brain has the ability to process a vast amount of information at incredibly high temporal scales, largely defined by the underlying neuronal architecture [1–3]. Therefore, characterizing these biological processes at a high temporal resolution is key to strengthen our understanding of brain function. However, current non-invasive approaches are often limited in either their temporal or spatial domain [4, 5]. In the worst case scenario, measurement parameters can only be obtained after the actual stimulation and the corresponding neuronal activity have ceased [6].

One of the most commonly used techniques is electroencephalography (EEG), which provides the highest temporal resolution in the millisecond range. However, EEG has limited spatial resolution, which is why blood-oxygen level dependent (BOLD) fMRI often represents the method of choice to assess human brain function in vivo with a good compromise between temporal (seconds) and spatial resolution (mm). As a drawback, the BOLD signal can only be indirectly related to neuronal activation [7] since it reflects a composite of blood flow, volume and oxygenation [6]. In contrast, positron emission tomography (PET) has a spatial resolution comparable to that of fMRI and provides unmatched molecular sensitivity and specificity. However, standard PET imaging is often limited to snapshots of the underlying molecular processes in the range of minutes to hours.

The recent introduction of functional PET (fPET) imaging offers a solution to this problem through the constant administration of the radioligand and repeated stimulation, similar to fMRI designs [8, 9]. By using the radiolabelled glucose analogue [18F]FDG, it enables the acquisition of baseline metabolism and task-specific changes in a single PET scan. Previous work with simultaneous assessment of fPET and fMRI has characterized the spatial agreement as well as unexpected differences between stimulation-induced changes in BOLD contrast and glucose metabolism by either acknowledging the different timescales of the signals [10] or by using a hierarchical design [11, 12]. Although the temporal resolution of single fPET time frames has been improved from minutes to 16 s [13] or even 6–12 s [14], this still impedes a direct comparison between BOLD and metabolic signals.

To overcome these issues we combined multiple recent advances in PET imaging. These include i) fPET with a bolus plus constant infusion protocol to increase the signal-to-noise ratio (SNR) of the metabolic signal [14]. ii) A novel dynamic filtering technique [15] with a sliding window was implemented, which enhances task-specific effects and enables reconstruction of high-temporal resolution fPET frames. iii) Finally, a truly synchronized acquisition of [18F]FDG fPET and BOLD fMRI data was employed during performance of the established working memory n-back paradigm. Notably, these developments are readily applicable to conventional PET scanner systems, thus enabling a widespread implementation without specific or costly hardware requirements. As a secondary aim, we evaluated the performance of a new-generation PET/CT scanner with time-of-flight imaging compared to an older, more widely available scanner in the context of fPET imaging.

Methods

Data acquisition, blood sampling, quantification of CMRGlu [10, 14, 16] and estimation of BOLD signal changes [17] were conducted analogous to our previous work, unless otherwise specified.

Experimental design

In this cross-sectional study 35 participants underwent either a simultaneous fPET/fMRI (n = 19) or separate fPET and fMRI scans (n = 16, 5.2 ± 5.8 days between scans) with the radiotracer [18F]FDG. During each of the scans an established working memory task was completed in a conventional block design. fPET started with an initial resting period of 8 min, which was followed by 12 min of task performance. The initial baseline was omitted for the MRI-only scan, but the task was otherwise identical. MRI acquisitions included a T1-weighted structural image and a BOLD sequence.

Cognitive task

The established n-back working memory task [18] was implemented in Psychtoolbox for Matlab. In this task letters were presented on a screen for 0.5 s, followed by an interstimulus interval of 1.9 s. During the 0-back control condition, a button press was required when the letter “X” appeared on the screen. For the 2-back condition, a button press was required when the current letter was the same as the one shown two letters before. Before each task block the corresponding instructions were shown for 2 s. The two conditions were shown in blocks of 36 s, each comprising 15 stimuli and 5 button presses (pseudo-randomized order). Resting periods were included at the beginning of the task (14 s), after each instruction (2 s) and after each task block (32 s). For each condition 5 blocks were presented in pseudo-randomized order, resulting in a total task duration of 12 min. The letters used in the task were phonologically similar in the German alphabet to avoid learning strategies [19, 20]. Small as well as capital letters were used equally (b, c, d, g, p, t, w). During all periods of rest, participants were instructed to look at a crosshair, relax and not to focus on anything in particular. Before data acquisition, participants completed a training session of each task condition inside the scanner to familiarize themselves with the task and the controls.

Participants

For this study 42 healthy participants were initially recruited, and data from 35 were used in the final analysis (mean age ± sd = 24.5 ± 4.4 years, 19 women). Dropout reasons included: blood sampling failures (n = 3), problems with the tracer administration (n = 3) and gross anatomical abnormalities (n = 1). Additionally, BOLD data could not be used for two participants due to technical reasons, blood glucose levels were not available for one subject and working memory performance data for another one due to failure of the response box. As no studies with such a high temporal resolution of fPET data are available (previously 6–12 s [14] and 16 s [13]), we aimed for a sample size that exceeds previous work with robust task activation by at least 50%. Each participant underwent an initial screening visit, where general health was ensured. Routine medical examinations included blood tests, electrocardiography, neurological testing and the structural clinical interview for DSM-IV. Female participants underwent pregnancy tests at the screening visit and before the PET and/or MRI scans. Further exclusion criteria were current and previous severe somatic, neurological and psychiatric conditions, substance abuse or medication, pregnancy or breastfeeding, MRI contraindications and previous study-related radiation exposure. At the screening visit, all participants provided written informed consent after detailed explanation of the study protocol. Participants were reimbursed and insured during the study. The study was approved by the ethics committee of the Medical University of Vienna (ethics numbers 1479/2015 and 2054/2020) and all procedures were carried out according to the Declaration of Helsinki.

Data acquisition

Simultaneous acquisition of fPET and fMRI data was done with a mMR scanner system. Separate acquisition was carried out with a Biograph Vision PET/CT and a Prisma MRI scanner (all Siemens Healthineers, Erlangen, Germany).

Participants had to fast at least 5.5 h before the start of the fPET scan (except for unsweetened water) [21]. The radiotracer [18F]FDG was administered in a bolus (816 ml/h for 1 min) plus constant infusion protocol (114.9 ml/h for 19 min, total of 50 ml) for the entire fPET scan using a perfusion pump (Syramed µSP6000, Arcomed, Regensdorf, Switzerland), which was kept in an MR-shield (UniQUE, Arcomed).

With the PET/MR scanner, BOLD fMRI was acquired simultaneously with fPET (EPI sequence, TE/TR = 30/2000 ms, flip angle = 90°, matrix size = 80 × 80, 34 slices, voxel size = 2.5 × 2.5 × 3.325 mm, 12 min). Furthermore, a structural MRI was recorded before fPET acquisition, to rule out gross anatomical abnormalities and for spatial normalization (T1-weighted MPRAGE sequence, TE/TR = 4.21/2200 ms, TI = 900 ms, flip angle = 9°, matrix size = 240 × 256, 160 slices, voxel size 1 × 1 × 1.1 mm, 7.72 min).

For the separate MRI acquisition, corresponding imaging data were acquired (EPI sequence: TE/TR = 30/2050 ms, flip angle = 78°, matrix size = 100 × 100, 35 slices, voxel size = 2.1 × 2.1 × 3.5 mm, 12 min; T1-weighted MPRAGE sequence: TE/TR = 2.91/2000 ms, TI = 900 ms, flip angle = 9°, matrix size = 240 × 256, 192 slices, voxel size 1 × 1 × 1 mm, 8 min).

Blood sampling

Before the fPET scans, blood glucose levels were determined using the average of 3 sequential measurements (Gluplasma). During the fPET scans, arterial blood samples were taken manually from the left radial artery. This included sampling every 20 s for the first three minutes, and thereafter at 4, 5, 8, 12, 16 and 20 min. We aimed to minimize the effects of blood sampling on task performance in several ways. First, samples after 5 min were taken remotely (outside 20 mT line) with a VAMP Plus system (Edwards Lifesciences), which enables continuous and gentle blood drawing and subsequent flushing with NaCl. This sampling technique allowed for minimal interaction with the participant, avoiding manipulation and movement of the participant’s arm for blood drawing. Second, it was carried out cautiously by an experienced and carefully trained physician to minimize the ability of subjects to recognize the sampling. Third, to minimize biasing task effects, the sampling was done on the left arm during baseline, whereas the task was executed with the right hand. Notably, none of the subjects reported that blood sampling was distracting them from task performance. Whole-blood activity was measured in a gamma counter (Wizard2, Perkin Elmer). Samples obtained after 3 min were centrifuged and the plasma activity was measured. Whole-blood data were linearly interpolated to match fPET frames and multiplied with the average plasma-to-whole-blood ratio to obtain the arterial input function.

fPET preprocessing

fPET data recorded at the PET/MR were corrected for attenuation using a database approach based on the T1-weighted structural image [22] and reconstructed to frames of 3, 6 and 12 s (matrix size = 344 × 344, 127 slices). PET/CT data were attenuation-corrected with an ultra low-dose CT scan obtained in the beginning of the fPET session (voltage: 120 kV, tube current 20 mA, Care-Dose4D, CarekV). Data were reconstructed to the same time frames as for the PET/MR scans (matrix size = 220 × 220, 159 slices). Preprocessing was carried out in SPM12 (https://www.fil.ion.ucl.ac.uk/spm/), which included motion correction (quality = 1, register to mean) and spatial normalization to MNI space with the T1-weighted structural image.

We implemented a novel filtering technique in Matlab to increase the low SNR inherent to high-temporal resolution PET data since short frames contain less radioactive counts. The filter represents a dynamic non-local means (NLM) filter [15] with a modified local patch selection. Generally, NLM filters create a weighted average of neighbouring voxels within a given search window in the spatial and temporal domain. These types of filters have been shown to provide an improved contrast-to-noise ratio in dynamic PET images compared with other more conventional denoising approaches [23]. The novel aspect implemented here is a sliding-window approach, instead of using the entire PET time course [15]. This confines the denoising of each voxel both temporally and spatially, preventing the baseline [18F]FDG uptake to drive the signal and thereby specifically strengthening task-specific effects. The local neighbourhood was selected with a sliding-window approach for each voxel using a spatiotemporal patch of 3 × 3 × 3 voxels in a search window of 11 × 11 × 11 voxels. The temporal patch comprised 15, 18 and 12 s for the 3, 6 and 12 s reconstructions, respectively. The weight is given by the similarity between voxels (i.e., Gaussian distance from centre voxel) as in the original description [15]. Data were subsequently smoothed with a 5 mm Gaussian kernel in SPM12. Together with the NLM filter, this corresponds approximately to a total kernel of 8 mm, which facilitates comparison with BOLD data and previous work.

Cerebral metabolic rate of glucose (CMRGlu)

Quantification of CMRGlu was carried out voxel-wise. Data were masked to include only gray matter by using a binarized gray matter segment of the SPM12 tissue probability map, thresholded at 0.1. To further improve the SNR in the temporal domain, a low pass filter (12-order finite impulse response filter) was used with a cutoff frequency of 1/18 s. This was selected in analogy to our previous work [9, 10, 14] as it represents half the task block duration. To separate task-specific effects from baseline metabolism, the general linear model (GLM) was used, which included one regressor for the baseline, one for each task condition and motion regressors. The baseline regressor was defined individually for each participant as the average fPET signal across all gray matter voxels, excluding a union of voxels activated during each participant’s individual BOLD fMRI response (2-back vs. rest, p < 0.001 uncorrected) and those identified in a meta-analysis [18]. This combination of a meta-analytical ground and individual activation ensures that task effects do not contaminate the baseline and yields an optimal model fit [14]. The task regressors were defined as a ramp function with slope = 1 kBq/frame. For the motion regressors, the realignment parameters obtained from SPM12 were subject to a principal component analysis and a variable amount of regressors were included using the elbow method, to ensure that the majority of explained variance is covered. Using the GLM output of the different conditions and the arterial input function, the Patlak plot was employed for absolute quantification of glucose metabolism. This yields the influx constant Ki, which reflects the combination of the individual rate constants K1-k31 Ki=K1*k3/(k2+k3)

Ki was then converted to the cerebral metabolic rate of glucose (CMRGlu)2 CMRGlu=Ki*Gluplasma/LC*100

with LC being the lumped constant, which was set to 0.89 [24]. This results in task-specific maps of Ki and CMRGlu, which were used for subsequent analysis.

Blood-oxygen level dependent (BOLD) signal changes

BOLD data were processed in SPM12 [17]. Images were corrected for slice timing effects (reference = middle slice) and head motion (quality = 1, register to mean), followed by spatial normalization to MNI space and spatial smoothing with an 8 mm Gaussian kernel. Task-related changes in the BOLD signal were computed with the GLM, including one regressor for each condition (0-back, 2-back) as well as potentially confounding signals (instructions, 24 motion parameters [25], white matter and cerebrospinal fluid). The autocorrelation option was set to FAST [26]. The contrast of interest was 2-back vs. rest to facilitate comparison with fPET data.

Statistical analyses

For statistical analyses data from n = 35 participants were used for [18F]FDG glucose metabolism and n = 33 for BOLD signal analyses (n = 2 missing due to technical reasons). For correlation between the signals in the temporal domain and association with task performance, n = 17 with simultaneous PET/MRI acquisition were used.

Spatial conjunction of task-specific neuronal activation

To assess the spatial overlap between changes in the BOLD signal and metabolic demands as induced by the 2-back working memory task, a statistical conjunction analysis was carried out in SPM12 [27]. BOLD and Ki maps were individually z-scored and included in a one-way ANOVA with each modality representing a “group” and an additional factor accounting for the different scanners (i.e., PET/MR vs. separate PET/CT and MRI). Ki was used instead of CMRGlu as blood glucose data was missing for one subject. The conjunction was corrected for multiple comparisons using family-wise error rate (FWE) at the voxel level (p < 0.05) and clusters with less than 5 voxels were removed. This was carried out for the three different fPET reconstructions of 3, 6 and 12 s. The Dice coefficient was used to assess the spatial agreement between the three resulting conjunction maps. For completeness and comparison with previous work [12, 27], inference was also computed at p < 0.05 FWE-corrected cluster level after an initial voxel threshold of p < 0.001 uncorrected.

Time-domain analyses

Task-specific signals were obtained as GLM residuals after removing all effects except those of the task. The resulting BOLD signal and [18F]FDG time series were extracted from significant voxels of the above spatial conjunction (p < 0.05 FWE-corrected voxel-level). For visualization, these were averaged across 2-back blocks, participants and voxels for a time window starting 12 s before the task began until 36 s after the task (total of 84 s). Next, the BOLD and metabolic signals were correlated within each subject to evaluate the individual coupling between the two imaging modalities. These correlation values were further associated with the working memory performance, specifically the average reaction time of correct responses during the 2-back condition. This was done separately for each brain region (three clusters of the right DLPFC, left primary motor cortex, right insula, left and right intraparietal sulcus). Spearman’s rho was used to account for an outlier with substantially slower reaction time (> 3*sd of the group) and corrected for the seven brain regions using the Bonferroni method. To assess the robustness of the brain-behaviour association, the time window was also constrained to -12… + 24 s, -6… + 12 s and -6… + 6 s, followed by re-calculation of the correlation. To rule out that individual associations between metabolic and BOLD signals are driven by motion, these values were correlated with framewise displacement as obtained from the BOLD realignment parameters [28].

Analyses of the posterior cingulate cortex (PCC)

We further investigated a previously reported spatial dissociation of glucose metabolism and BOLD responses in the PCC [12]. First, separate maps with a negative BOLD response and a positive metabolic response during the 2-back working memory task were computed (p < 0.05 FWE-corrected cluster level following p < 0.001 uncorrected voxel level). Specific attention was given to the PCC, as we recently showed an opposite influence of the FPN and DAN metabolism onto the PCC [27]. Thus, the difference in task-specific Ki between FPN and DAN during the 2-back working memory condition was computed. Percent signal change was used to enable better comparison between the different networks. Based on this difference, participants were split into three groups comprising those where metabolism of FPN < DAN (lower quartile, low group), FPN = DAN (interquartile range, balanced group) and FPN > DAN (upper quartile, high group). [18F]FDG metabolic time series were then extracted for the PCC region [27] and the three groups, which were compared using a repeated measures ANOVA, where the interaction between groups by time served as the relevant outcome parameter. Similarly, BOLD signal time courses were extracted from the PCC and compared in a repeated measures ANOVA.

Scanner comparison

The secondary aim was to investigate differences between PET scanners, here we compared a widely used mMR PET/MR to a new generation Vision PET/CT with time-of-flight (both Siemens) in the feasibility to compute high-temporal resolution fPET. For the subsequent analyses, data were extracted from regions with significant activation during the 2-back condition (p < 0.05 FWE corrected voxel-level, Fig. 1a). In the temporal domain the identifiability of task-specific changes in glucose metabolism is given by t-values of the GLM analysis, i.e., the estimated effect size divided by the standard error [14]. A two-sample t-test between the scanners was performed. Next, Ki percent signal change was compared between the scanners, again with a two-sample t-test.Fig. 1 Spatial correspondence of BOLD fMRI changes and [18F]FDG fPET glucose metabolism induced by the 2-back task condition. Statistical conjunction analysis between the two signals yields robust task effects in several brain regions involved in working memory [18], such as the right dorsolateral prefrontal cortex, left primary motor cortex and bilateral anterior insula (blue: p < 0.05 FWE-corrected voxel level, yellow: p < 0.05 FWE-corrected cluster level following p < 0.001 uncorrected voxel level). The effects were stable for different temporal resolutions of the fPET data, namely 3s (a), 6s (b) and 12s (c). Dice coefficients between the activation patterns obtained with the three temporal resolutions were between 0.71 and 0.78. Voxel-level corrected regions of the 3s data (i.e., blue regions in a) were used for subsequent temporal assessment. MNI coordinates for slices are z = -1, 17 and 35 mm. Left is left

Simulations

Simulations were carried out to assess whether the observed time course of the [18F]FDG metabolic signal matches theoretical calculations. The two-tissue compartment model was implemented in Matlab with an average arterial input function and rate constants as observed in our previous work at resting state (K1 = 0.0827 ml/cm3/min, k2 = 0.0771 min−1, k3 = 0.0629 min−1, k4 = 0 min−1) [14, 16]. Task-specific changes in rate constants were introduced to simulate physiologically plausible increases in Ki (see Eq. 1). Previous work showed task-induced changes in Ki around 22% [9, 10], but for completeness also increases of 12 and 33% were introduced. These values were obtained by the following changes in rate constants i) K1 = 10% and k3 = 10%, ii) K1 = 12% and k3 = 5%, iii) K1 = 8% and k3 = 15%, iv) K1 = 15% and k3 = 12%, v) K1 = 5% and k3 = 8%. In all cases k2 = -K1 [29–32]. Please see discussion for a neurophysiological rationale of these parameter variations. In addition, the simulation was repeated with baseline rate constants of previous work (K1 = 0.1 ml/cm3/min, k2 = 0.15 min−1, k3 = 0.08 min−1, k4 = 0 min−1) [8].

Results

Spatial domain

Statistical conjunction analysis was used to assess whether working memory performance induces similar spatial activations between metabolic and BOLD signal changes. Increases across both signals were observed in the right dorsolateral prefrontal cortex, the left primary motor cortex as well as the intraparietal sulcus and anterior insula bilaterally during the 2-back condition (p < 0.05 FWE corrected). This activation pattern was robust for all temporal resolutions of 3, 6 and 12 s with Dice coefficients between 0.71 and 0.78 (Fig. 1). Across these regions the average increase was Ki = 0.011 ± 0.005 ml/cm3/min and CMRGlu = 6.74 ± 2.78 µmol/100 g/min, which corresponds to a change of 27.15 ± 9.91%.

In a similar conjunction analysis we observed a dissociation between signals with a negative BOLD response and increased glucose metabolism in the left frontal cortex (x/y/z = -18/54/26 mm, p < 0.05 FWE corrected). This closely matches previous work, showing a dissociation at the borders between default mode and frontoparietal networks [12]. The inverse effect (i.e., positive BOLD response and decreased metabolism) was observed in the cerebellar vermis.

Temporal domain

The time courses of the [18F]FDG metabolic and BOLD hemodynamic signals were then extracted from significant brain regions (Fig. 1, p < 0.05 FWE corrected voxel level) and averaged across all 2-back task blocks and participants. As expected, the BOLD signal changes exhibited a slightly delayed plateau with respect to task performance (Fig. 2a). The high-temporal resolution [18F]FDG signal showed a continuous increase during task execution due to the irreversible binding of the radioligand. Interestingly, also the metabolic signal returned to baseline within 10 s, which was most visible for the 3 s reconstruction (Fig. 2b). Therefore, the highest temporal resolution was used for all subsequent computations, further enabling direct comparison with the BOLD signal. Additional simulations were carried out to assess the temporal dynamics of brain metabolism in response to task-induced changes in the corresponding rate constants of the two-tissue compartment model. The simulated time courses showed a remarkable match with the observed task response (Fig. 2c). Specifically, simulations yielded a linear increase of the signal, which returned to baseline after stimulation ceased. The simulation results were robust across a physiologically plausible range of different kinetics, which included changes in Ki of 12, 22 and 33%, elicited by variation of K1, k2 and k3 between 5 and 15% (Fig. 2d). The different simulations only changed the amplitude of the time course, but not the overall shape. This was also the case when using rate constants from previous work [8].Fig. 2 Concurrent temporal changes in BOLD fMRI and [18F]FDG fPET signals during working memory performance across all active regions. a) Task-induced changes in the BOLD signal show a plateau of activation as expected from the block design. b) Task-specific time course of [18F]FDG glucose metabolism resulted in a constant increase during task performance (due to the irreversible binding of the radiotracer) and rapid return to baseline afterwards (i.e., a steady activity level). The time course was generally similar for the different temporal resolutions of 3 s (black), 6 s (blue) and 12 s (green), but the 3 s data showed highest ability to reconstruct the return to baseline. c) Simulation of task-related changes in the [18F]FDG signal was in close agreement with the shape of the experimental observation in b across different variations in rate constants. Colors of lines match those of symbols in d. d) Matrix of changes in Ki induced by variation of rate constants K1, k2 and k3 from 0 to 20%. Changes in rate constants were chosen to yield physiological increases in the net influx constant Ki around 22% (black line as well as black, yellow and magenta circles) as observed in this and previous work [9]. In addition, rate constants yielding Ki = 12 and 33% (i.e., orthogonal to the expected 22%-line) were investigated (orange and red x, respectively). For all simulations k2 = -K1. Gray areas in a-c indicate the time of the task performance. Data in a and b represent the average across participants, 2-back task blocks and regions obtained from the statistical conjunction analysis (blue regions in brain slices, p < 0.05 FWE-corrected voxel level, identical to Fig. 1a). Data in b and c show task effects after subtraction of baseline metabolism

Next, we investigated the individual association between metabolic and hemodynamic signals and their relationship with task performance for the sample that underwent truly simultaneous PET/MR imaging (n = 17). Across all task regions, we observed a moderate agreement between both modalities (|r|= 0.31 ± 0.19). Moreover, the association of BOLD and metabolic signals in the left primary motor cortex correlated with the reaction time of correct 2-back button presses, except for one outlier whose reaction time was more than 3 standard deviations higher than the group average (rho = -0.69, p = 0.02 corrected, time window = -12 to + 36 s with respect to task performance, Fig. 3d). The correlation with task performance was independent of the time windows (-12 to + 24 s, -6 to + 12 s and -6 to + 6 s) used for the multimodal signal correlations with rho = -0.68…-0.65 (p = 0.024…0.043 corrected). The corresponding time series are exemplarily shown in Fig. 3a–c, indicating a similar temporal profile of the metabolic and BOLD signals with individual delays. Importantly, the association between the two signals was not driven by motion artifacts (rho = -0.15, p = 0.57).Fig. 3 Temporal associations between signals of fMRI BOLD (blue) and fPET glucose metabolism (CMRGlu, red) in the left primary motor cortex. A–c depict three representative examples of individual time courses of the BOLD signal and [18F]FDG glucose metabolism, averaged across the 2-back working memory task blocks. For CMRGlu, the first derivative is shown to account for the irreversible accumulation of the radiotracer. The three time courses show that a) the metabolic signal may either precede the BOLD signal, b) both signal changes occur simultaneously or c) that the BOLD signal precedes that of CMRGlu. Data in a–c were z-transformed to account for different amplitudes of the signal. Gray areas indicate the time of the task performance. d) The individual temporal correlation between the CMRGlu and BOLD signals was associated with the corresponding response time of correct button presses during the 2-back working memory task (rho = -0.69, p = 0.02 corrected)

Finally, we investigated the temporal domain in relation to previous findings of a spatial dissociation between BOLD changes and glucose metabolism in the posteromedial default mode network (DMN) during working memory [12]. Consistent with those findings we observed a negative BOLD response in the posterior cingulate cortex (PCC) during the 2-back working memory execution. However, no change in metabolic demands was found when averaged across the entire sample (Fig. 4b, p < 0.05 FWE-corrected cluster level following p < 0.001 uncorrected voxel level). This dissociation was also evident when using the PCC as region of interest (BOLD beta = -0.93 ± 0.29, p < 10–17; Ki = 0.0006 ± 0.0072 ml/cm3/min, p = 0.6).Fig. 4 fPET metabolic and fMRI BOLD signals in the posterior cingulate cortex (PCC) during the 2-back working memory task. a) Participants were split into three groups according to significant changes in CMRGlu in the frontoparietal (FPN) and dorsal attention network (DAN, p < 0.05 FWE-corrected voxel level, Fig. 1a) as these networks drive CMRGlu changes in the PCC [27]. By definition, the quartile split separates 50% of participants with similar changes in these networks (green, balanced group), while 25% had lower task-specific CMRGlu in FPN than DAN (blue, low) and 25% had higher CMRGlu in FPN than DAN (red, high). b) Task-specific decreases in the BOLD signal (blue) and increases in CMRGlu (red) during working memory performance (p < 0.05 FWE-corrected cluster level following p < 0.001 uncorrected voxel level). Replicating previous work [12], BOLD decreases where pronounced in the PCC, while CMRGlu did not change significantly. The white area indicated by the cross-hair represents the PCC cluster of our recent study [27], which was used for c and d. c) The time course of [18F]FDG glucose metabolism in the PCC was significantly different between the three groups (group*time interaction p = 0.0016). Participants with lower CMRGlu in FPN than DAN exhibited a decrease in PCC CMRGlu (blue) during working memory, while participants with higher CMRGlu in FPN than DAN showed a delayed positive deflection in the metabolic signal. d) BOLD signals in the PCC showed similar decreases during the task for all three groups (interaction p = 0.9). Colors in a, c and d represent the same groups of participants. Boxplots in a indicate median values (center line), upper and lower quartiles (box limits) and 1.5*interquartile range (whiskers). Data in c show task effects after subtraction of baseline metabolism. Gray areas indicate the time of the task performance

Given the corresponding metabolic influence of frontoparietal (FPN) and dorsal attention networks (DAN) on the PCC [27], participants were split into quartiles according to the difference in task-specific Ki between the FPN and DAN (Fig. 4a). This resulted in three distinct groups with the metabolism of FPN < DAN (low, 25% of participants), FPN = DAN (balanced, 50%) and FPN > DAN (high, 25%). Evaluating the temporal profile of the PCC between these groups revealed an interaction effect (p = 0.0016, Fig. 4c), with significant post-hoc differences between high vs. low participants (p < 10–5) and low vs. balanced (p = 0.025). In detail, the low group (i.e., with lower metabolism in FPN than DAN) was characterized by a constant decrease in PCC metabolism during task execution, followed by an increase afterwards. In contrast, the high group only exhibited a delayed increase in the metabolic signal, whereas the balanced group showed almost no variation in their metabolic response. Again, the temporal difference in the metabolic profile between the three groups was consistent across all time windows (p = 0.0005…0.040). However, the BOLD time series of the PCC showed robust decreases for all three groups, without any significant difference (Fig. 4d, group-by-time interaction p = 0.9).

Comparison of PET systems

In the temporal domain, there was no significant difference in the ability to identify task-specific changes in glucose metabolism between the scanners (PET/MR t = 3.13 ± 0.98, PET/CT t = 3.71 ± 1.49, p = 0.17). Similarly, Ki percent signal change of regions involved in the task (Fig. 1a, p < 0.05 FWE-corrected voxel level) was not significantly different between the scanners (PET/MR = 29.79 ± 10.22%, PET/CT = 25.54 ± 10.87%, p = 0.24).

Discussion

In this work we employed high temporal resolution fPET with 3 s to investigate the dynamics of glucose metabolism during working memory performance. The metabolic response of the [18F]FDG signal was characterized by a constant increase during task execution and followed by a rapid return to baseline, which was consistent with simulations. Using fully synchronized fPET/fMRI we further demonstrate the feasibility to assess the individual coupling of the hemodynamic and metabolic response, whose effect in the primary motor cortex was related to working memory reaction time. Finally, analysis of PCC signals revealed that only specific subgroups exhibited decreased metabolism during the task while others showed an increased response, which was dependent on the corresponding task-positive metabolic demands. Together, these findings highlight the value of assessing brain energy metabolism in the temporal domain, which provides unique information not accessible to conventional PET imaging.

Regardless of the reconstructed frame length (3, 6 and 12 s), we observed spatially robust activations induced by working memory performance which were in line with a meta-analysis of BOLD signal changes [18]. Moreover, the average increase in glucose metabolism was consistent with previous work using lower, more typical, temporal resolutions of 30 to 60 s, both in terms of absolute units of Ki and CMRGlu as well as percent signal changes [9, 10]. This indicates that the lower SNR inherent to shorter time frames is compensated by the increased number of data points for the general linear model estimation, with the additional benefit to investigate fast metabolic changes. Altogether, these observations provide confidence for the validity to quantify task-induced changes in glucose metabolism at a high temporal resolution of 3 s. Compared to our previous work [14], this is now possible through further improvement in the protocol, specifically the combination of tracer administration as bolus + constant infusion and the novel implementation of the dynamic filter. In contrast, the use of constant infusion only (i.e., without an initial bolus) [9, 12, 13] does not seem to provide sufficient SNR to reconstruct data with such a high temporal resolution (see e.g., Fig. 2b in [33]). Regarding the different PET systems, the new-generation PET/CT scanner with time-of-flight imaging showed similar performance compared to the PET/MR. This may be related to the similar sensitivity of the two systems (16.4 kcps/MBq [34] and 15.0 kcps/MBq [35], respectively) and the fact that the benefit of time-of-flight is less pronounced for smaller volumes such as the head [36], which improved the spatial SNR but had little to no effect on the temporal resolution.

Our improved protocol allowed us to depict the corresponding time course of the activated brain regions in unprecedented detail. This was characterized by fast changes during cognitive performance, namely an expected constant increase in the [18F]FDG signal during the task but also a rapid return to baseline within 10 s once stimulation ceased. These fast metabolic changes are supported by animal studies. Specifically, almost immediate differences in extracellular availability of ATP [37] and glucose [38, 39] upon neuronal activation have been reported across different animals and approaches. Although these methods differ with respect to the return to baseline, two of these studies also matched our work in this respect with 12 ± 3 s [39] and a range between 10 and 100 s [37]. As these differences are likely related to the specific approach and species, the exact metabolic response obtained with [18F]FDG as well as its variation across participants and paradigms need to be determined in future studies. Nevertheless, the time course of the task-induced metabolic signal was supported by previous [40] and our own simulations. In contrast to other work simulating slow metabolic changes due to variation of k3 only [8, 12], we employed a complete simulation of the irreversible two-tissue compartment model with changes in all three relevant rate constants. In particular, the variation of K1 is feasible as this parameter is closely related to cerebral blood flow (CBF). CBF in turn is well known to change with neuronal activation and also represents a major driver of the BOLD signal [41, 42]. Furthermore, glucose as well as [18F]FDG are subject to transport across the blood brain barrier (BBB) that is facilitated by the glucose transporter 1 (GLUT1) carrier protein. In this transport system, however, influx from blood plasma to the brain (K1) will occupy the carrier protein and thus block its availability for the reverse transport (k2), i.e., influx and efflux rate constants are inversely related and affect each other [29–32]. Although k3 most closely resembles hexokinase activity and the first step in glycolysis, glucose transport across the BBB reflected by K1 and k2 represents an essential (but often neglected) aspect to meet increased metabolic demands of neuronal activation [43–45]. Furthermore, all three rate constants contribute to the [18F]FDG signal and the final quantified CMRGlu by definition, thus, supporting their combined variation in the model [40] and the resulting temporal profile of glucose metabolism.

The coupling between the BOLD signal and CMRGlu has previously been shown to increase during task execution at the spatial level [10]. The current work extends this finding to the time domain, where stronger temporal coupling between metabolic and BOLD time series was associated with faster working memory performance. Both of these signals are related to neuronal activity via different yet intertwined mechanisms [44]. The BOLD contrast results from disproportional increases in oxygen consumption (CMRO2) and CBF [41, 42]. The rise in CBF is related to glutamate release during neuronal activation [46, 47] to maintain the supply of nutrients, termed neurovascular coupling. On the other hand, neurometabolic coupling links glutamate release to an increase in glucose consumption to meet the energy demands for the reversal of ion gradients [48–50]. These changes in CBF, CMRO2, CMRGlu and thus also the neurovascular and the neurometabolic responses have previously been considered as parallel processes driven by neuronal activation, instead of being a serial connection of events [39, 51]. This may explain the observed delay between metabolic and BOLD signals in the primary motor cortex, which varied across participants. We therefore put forward the hypothesis that the individual differences in the coupling between BOLD and [18F]FDG signals potentially informs us about the integration of the metabolic and hemodynamic interplay in neuronal processing. Possible underlying reasons for these differences could be inter-individual differences in the hemodynamic response [52] as well as regional variability in the neurovascular unit [53, 54] and BBB glucose transport [55], which probably further translates to individual differences across participants and may also be present in a similar fashion for neurometabolic coupling. On the other hand, dissociation between the signals may provide promising information about alterations in different brain disorders. For instance, in healthy aging increased BOLD-derived activations were not matched by glucose metabolism [56]. Moreover, patients with Alzheimer’s disease suffer from both impaired neurovascular coupling [57] and widespread decreased metabolism [58].

Our analysis in the temporal domain also provides further insight into the divergence between glucose metabolism and BOLD deactivations in the PCC [12]. In line with our previous work in the spatial domain [27], we observed that the metabolic time course is dependent on energy demands of the corresponding task-positive networks involved in the cognitive process. More specifically, individuals with lower CMRGlu in FPN than DAN exhibited consistent negative deflections in both metabolic and BOLD time series. In contrast, higher CMRGlu in FPN than DAN did not lead to relevant changes in the metabolic signal during working memory. Since BOLD signal decreases without affecting metabolism seem unlikely [50], we speculate that metabolically expensive suppression [12, 27] and downregulation of neuronal signalling (presumably mediated by GABA and glutamate, respectively) [27] occur simultaneously, thus resulting in BOLD deactivations but a net metabolism change around zero. Interestingly, right after task execution both groups showed an increase in glucose metabolism. This may indicate that the actual return of the BOLD signal to baseline also requires energy. The effect may potentially be related to increased glutamate signalling, which leads to increases in both glucose metabolism and the BOLD signal (see above).

Since the BOLD signal in the PCC showed a negative response throughout, we assume that the ratio between CBF and CMRO2 is also similar for all participants. As a consequence, the variation of glucose metabolism across individuals and time further implies that also the oxygen-to-glucose index (OGI) differs in the same manner. In particular, those individuals without a decreased metabolic response (and thus higher CMRGlu in FPN than DAN) would exhibit a decreased OGI in the PCC during task performance, indicating that energy is supplied by aerobic glycolysis [59]. This is supported by recent work suggesting that immediate energy demands are preferably met by aerobic glycolysis [60]. However, to verify our hypothesis acquisition of the actual CMRO2 time series is required to fully understand the mechanisms of neurometabolic coupling [51].

In conclusion, being the first study that employs such a high temporal resolution of fPET data, our work does not claim to provide definitive answers to the observed phenomena. Rather, the approach aims to spark future discussions on the relationship between hemodynamic and metabolic signals of the human brain. Due to the relatively small sample of simultaneously acquired PET/MR data, more cases are required to thoroughly assess the associations between hemodynamics and metabolic demands. Future studies should aim to identify the underlying coupling mechanisms, particularly when the signals diverge, the delay between BOLD and metabolic time series as well as potential dissociations in brain disorders. Furthermore, the performance of challenging cognitive tasks may not be feasible in certain brain disorders with impaired cognitive and/or physical abilities, such as genetic disorders with mental retardation, motor or neurological disorders, schizophrenia, Alzheimer’s disease or traumatic brain injury, as this may bias the outcome of brain imaging. However, our approach may allow identifying minor impairments at the early stages of the disease, e.g., prodromal stage of schizophrenia or mild cognitive impairment, which is of paramount clinical and diagnostic relevance. In addition, the presented technique allows assessment of [18F]FDG fPET and BOLD fMRI signals with a variety of task settings or difficulties and also during resting state, i.e., where no explicit cognitive abilities are required. Thus, the direct comparison of these signals in the temporal domain offers numerous possibilities to investigate these aspects in health and disease, thereby introducing a previously unseen perspective to our understanding of brain function.

Acknowledgements

We thank the graduated team members and the diploma students of the Neuroimaging Lab (NIL, head: R. Lanzenberger) as well as the clinical colleagues from the Department of Psychiatry and Psychotherapy for clinical and/or administrative support. In detail, we would like to thank S. Kasper, K. Papageorgiou, L. Silberbauer, C. Schmidt, B. Eggersdorfer, J. Unterholzner and V. Popper for medical support, L. Rischka for acquisition and analysis support, V. Ritter and C. Wotawa for subject recruitment. We are further grateful to W. Wadsak, V. Pichler, G. Karanikas, W. Langsteger and the radioligand synthesis team from the Department of Biomedical Imaging and Image-guided Therapy, Division of Nuclear Medicine for acquisition support and supervision. The scientific project was performed with the support of the Medical Imaging Cluster of the Medical University of Vienna.

Author contribution

Study design: A.H., R.L., T. T-W., M.H.

Data acquisition: C.V., L.N., G.M.G., S.K., A.K., M.B.R., A.H.

Methods: M.B.R., A.H., P.F.

Data analysis: M.B.R., A.H.

Manuscript preparation: A.H., M.B.R., G.M.G., S.K., R.L.

All authors discussed the implications of the findings and approved the final version of the manuscript.

Funding

Open access funding provided by Austrian Science Fund (FWF). This research was funded in whole, or in part, by the Austrian Science Fund (FWF, KLI 610, PI: A. Hahn; and KLI 1006, PI: R. Lanzenberger), and the Vienna Science and Technology Fund (WWTF) [10.47379/CS18039], Co-PI: R. Lanzenberger. For the purpose of open access, the author has applied a CC BY public copyright license to any Author Accepted Manuscript version arising from this submission. M.B. Reed is a recipient of a DOC fellowship of the Austrian Academy of Sciences at the Department of Psychiatry and Psychotherapy, Medical University of Vienna.

Data and code availability

Raw data will not be publicly available due to reasons of data protection. Processed data and custom code can be obtained from the corresponding author with a data-sharing agreement, approved by the departments of legal affairs and data clearing of the Medical University of Vienna.

Declarations

Ethics approval

Participants were reimbursed and insured during the study. The study was approved by the ethics committee of the Medical University of Vienna (ethics numbers 1479/2015 and 2054/2020) and all procedures were carried out according to the Declaration of Helsinki.

Consent to participate

At the screening visit, all participants provided written informed consent after detailed explanation of the study protocol.

Conflict of interest

RL received investigator-initiated research funding from Siemens Healthcare regarding clinical research using PET/MR. He is a shareholder of the start-up company BM Health GmbH since 2019. M. Hacker received consulting fees and/or honoraria from Bayer Healthcare BMS, Eli Lilly, EZAG, GE Healthcare, Ipsen, ITM, Janssen, Roche, and Siemens Healthineers. All other authors report no conflict of interest in relation to this study.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Ramsey NF Jansma JM Jager G Van Raalten T Kahn RS Neurophysiological factors in human information processing capacity Brain 2004 127 517 525 10.1093/brain/awh060 14691061
Ramsey NF, Jansma JM, Jager G, Van Raalten T, Kahn RS. Neurophysiological factors in human information processing capacity. Brain. 2004;127:517–25. 10.1093/brain/awh060.14691061 10.1093/brain/awh060
2. Marois R Ivanoff J Capacity limits of information processing in the brain Trends Cogn Sci 2005 9 296 305 10.1016/j.tics.2005.04.010 15925809
Marois R, Ivanoff J. Capacity limits of information processing in the brain. Trends Cogn Sci. 2005;9:296–305. 10.1016/j.tics.2005.04.010.15925809 10.1016/j.tics.2005.04.010
3. Li T Zheng Y Wang Z Zhu DC Ren J Liu T Brain information processing capacity modeling Sci Rep 2022 12 2174 10.1038/s41598-022-05870-z 35140253
Li T, Zheng Y, Wang Z, Zhu DC, Ren J, Liu T, et al. Brain information processing capacity modeling. Sci Rep. 2022;12:2174. 10.1038/s41598-022-05870-z.35140253 10.1038/s41598-022-05870-z
4. Shah NJ, Oros-Peusquens AM, Arrubla J, Zhang K, Warbrick T, Mauler J, et al. Advances in multimodal neuroimaging: hybrid MR-PET and MR-PET-EEG at 3 T and 9.4 T. J Magn Reson. 2013;229:101–15. 10.1016/j.jmr.2012.11.027.
5. Sarraf S Sun J Advances in Functional Brain Imaging: A Comprehensive Survey for Engineers and Physical Scientists Int J Adv Res 2016 4 640 660 10.21474/ijar01/1259
Sarraf S, Sun J. Advances in Functional Brain Imaging: A Comprehensive Survey for Engineers and Physical Scientists. Int J Adv Res. 2016;4:640–60. 10.21474/ijar01/1259.10.21474/ijar01/1259
6. Norris DG Principles of magnetic resonance assessment of brain function J Magn Reson Imaging 2006 23 794 807 10.1002/jmri.20587 16649206
Norris DG. Principles of magnetic resonance assessment of brain function. J Magn Reson Imaging. 2006;23:794–807. 10.1002/jmri.20587.16649206 10.1002/jmri.20587
7. Heeger DJ Ress D What does fMRI tell us about neuronal activity? Nat Rev Neurosci 2002 3 142 151 10.1038/nrn730 11836522
Heeger DJ, Ress D. What does fMRI tell us about neuronal activity? Nat Rev Neurosci. 2002;3:142–51.11836522 10.1038/nrn730
8. Villien M Wey HY Mandeville JB Catana C Polimeni JR Sander CY Dynamic functional imaging of brain glucose utilization using fPET-FDG Neuroimage 2014 100 192 199 10.1016/j.neuroimage.2014.06.025 24936683
Villien M, Wey HY, Mandeville JB, Catana C, Polimeni JR, Sander CY, et al. Dynamic functional imaging of brain glucose utilization using fPET-FDG. Neuroimage. 2014;100:192–9.24936683 10.1016/j.neuroimage.2014.06.025
9. Hahn A Gryglewski G Nics L Hienert M Rischka L Vraka C Quantification of Task-Specific Glucose Metabolism with Constant Infusion of 18F-FDG J Nucl Med 2016 57 1933 1940 10.2967/jnumed.116.176156 27390156
Hahn A, Gryglewski G, Nics L, Hienert M, Rischka L, Vraka C, et al. Quantification of Task-Specific Glucose Metabolism with Constant Infusion of 18F-FDG. J Nucl Med. 2016;57:1933–40. 10.2967/jnumed.116.176156.27390156 10.2967/jnumed.116.176156
10. Hahn A Breakspear M Rischka L Wadsak W Godbersen GM Pichler V Reconfiguration of functional brain networks and metabolic cost converge during task performance Elife 2020 9 e52443 10.7554/elife.52443 32314956
Hahn A, Breakspear M, Rischka L, Wadsak W, Godbersen GM, Pichler V, et al. Reconfiguration of functional brain networks and metabolic cost converge during task performance. Elife. 2020;9: e52443. 10.7554/elife.52443.32314956 10.7554/elife.52443
11. Jamadar SD Ward PG Li S Sforazzini F Baran J Chen Z Simultaneous task-based BOLD-fMRI and [18-F] FDG functional PET for measurement of neuronal metabolism in the human visual cortex Neuroimage 2019 189 258 266 10.1016/j.neuroimage.2019.01.003 30615952
Jamadar SD, Ward PG, Li S, Sforazzini F, Baran J, Chen Z, et al. Simultaneous task-based BOLD-fMRI and [18-F] FDG functional PET for measurement of neuronal metabolism in the human visual cortex. Neuroimage. 2019;189:258–66. 10.1016/j.neuroimage.2019.01.003.30615952 10.1016/j.neuroimage.2019.01.003
12. Stiernman LJ, Grill F, Hahn A, Rischka L, Lanzenberger R, Panes Lundmark V, et al. Dissociations between glucose metabolism and blood oxygenation in the human default mode network revealed by simultaneous PET-fMRI. Proc Natl Acad Sci U S A. 2021;118. 10.1073/pnas.2021913118.
13. Jamadar SD Ward PGD Liang EX Orchard ER Chen Z Egan GF Metabolic and Hemodynamic Resting-State Connectivity of the Human Brain: A High-Temporal Resolution Simultaneous BOLD-fMRI and FDG-fPET Multimodality Study Cereb Cortex 2021 31 2855 2867 10.1093/cercor/bhaa393 33529320
Jamadar SD, Ward PGD, Liang EX, Orchard ER, Chen Z, Egan GF. Metabolic and Hemodynamic Resting-State Connectivity of the Human Brain: A High-Temporal Resolution Simultaneous BOLD-fMRI and FDG-fPET Multimodality Study. Cereb Cortex. 2021;31:2855–67. 10.1093/cercor/bhaa393.33529320 10.1093/cercor/bhaa393
14. Rischka L Gryglewski G Pfaff S Vanicek T Hienert M Klobl M Reduced task durations in functional PET imaging with [(18)F]FDG approaching that of functional MRI Neuroimage 2018 181 323 330 10.1016/j.neuroimage.2018.06.079 29966719
Rischka L, Gryglewski G, Pfaff S, Vanicek T, Hienert M, Klobl M, et al. Reduced task durations in functional PET imaging with [(18)F]FDG approaching that of functional MRI. Neuroimage. 2018;181:323–30. 10.1016/j.neuroimage.2018.06.079.29966719 10.1016/j.neuroimage.2018.06.079
15. Gal Y Mehnert AJ Bradley AP McMahon K Kennedy D Crozier S Denoising of dynamic contrast-enhanced MR images using dynamic nonlocal means IEEE Trans Med Imaging 2010 29 302 310 10.1109/TMI.2009.2026575 19605318
Gal Y, Mehnert AJ, Bradley AP, McMahon K, Kennedy D, Crozier S. Denoising of dynamic contrast-enhanced MR images using dynamic nonlocal means. IEEE Trans Med Imaging. 2010;29:302–10. 10.1109/TMI.2009.2026575.19605318 10.1109/TMI.2009.2026575
16. Klug S Godbersen GM Rischka L Wadsak W Pichler V Klobl M Learning induces coordinated neuronal plasticity of metabolic demands and functional brain networks Communications biology 2022 5 428 10.1038/s42003-022-03362-4 35534605
Klug S, Godbersen GM, Rischka L, Wadsak W, Pichler V, Klobl M, et al. Learning induces coordinated neuronal plasticity of metabolic demands and functional brain networks. Communications biology. 2022;5:428. 10.1038/s42003-022-03362-4.35534605 10.1038/s42003-022-03362-4
17. Reed MB Vanicek T Seiger R Klobl M Spurny B Handschuh P Neuroplastic effects of a selective serotonin reuptake inhibitor in relearning and retrieval Neuroimage 2021 236 118039 10.1016/j.neuroimage.2021.118039 33852940
Reed MB, Vanicek T, Seiger R, Klobl M, Spurny B, Handschuh P, et al. Neuroplastic effects of a selective serotonin reuptake inhibitor in relearning and retrieval. Neuroimage. 2021;236: 118039. 10.1016/j.neuroimage.2021.118039.33852940 10.1016/j.neuroimage.2021.118039
18. Mencarelli L Neri F Momi D Menardi A Rossi S Rossi A Stimuli, presentation modality, and load-specific brain activity patterns during n-back task Hum Brain Mapp 2019 40 3810 3831 10.1002/hbm.24633 31179585
Mencarelli L, Neri F, Momi D, Menardi A, Rossi S, Rossi A, et al. Stimuli, presentation modality, and load-specific brain activity patterns during n-back task. Hum Brain Mapp. 2019;40:3810–31. 10.1002/hbm.24633.31179585 10.1002/hbm.24633
19. Harvey PO Fossati P Pochon JB Levy R Lebastard G Lehericy S Cognitive control and brain resources in major depression: an fMRI study using the n-back task Neuroimage 2005 26 860 869 10.1016/j.neuroimage.2005.02.048 15955496
Harvey PO, Fossati P, Pochon JB, Levy R, Lebastard G, Lehericy S, et al. Cognitive control and brain resources in major depression: an fMRI study using the n-back task. Neuroimage. 2005;26:860–9.15955496 10.1016/j.neuroimage.2005.02.048
20. Schoning S Zwitserlood P Engelien A Behnken A Kugel H Schiffbauer H Working-memory fMRI reveals cingulate hyperactivation in euthymic major depression Hum Brain Mapp 2009 30 2746 2756 10.1002/hbm.20702 19086021
Schoning S, Zwitserlood P, Engelien A, Behnken A, Kugel H, Schiffbauer H, et al. Working-memory fMRI reveals cingulate hyperactivation in euthymic major depression. Hum Brain Mapp. 2009;30:2746–56. 10.1002/hbm.20702.19086021 10.1002/hbm.20702
21. Guedj E Varrone A Boellaard R Albert NL Barthel H van Berckel B EANM procedure guidelines for brain PET imaging using [(18)F]FDG, version 3 Eur J Nucl Med Mol Imaging 2022 49 632 651 10.1007/s00259-021-05603-w 34882261
Guedj E, Varrone A, Boellaard R, Albert NL, Barthel H, van Berckel B, et al. EANM procedure guidelines for brain PET imaging using [(18)F]FDG, version 3. Eur J Nucl Med Mol Imaging. 2022;49:632–51. 10.1007/s00259-021-05603-w.34882261 10.1007/s00259-021-05603-w
22. Burgos N Cardoso MJ Thielemans K Modat M Pedemonte S Dickson J Attenuation correction synthesis for hybrid PET-MR scanners: application to brain studies IEEE Trans Med Imaging 2014 33 2332 2341 10.1109/TMI.2014.2340135 25055381
Burgos N, Cardoso MJ, Thielemans K, Modat M, Pedemonte S, Dickson J, et al. Attenuation correction synthesis for hybrid PET-MR scanners: application to brain studies. IEEE Trans Med Imaging. 2014;33:2332–41. 10.1109/TMI.2014.2340135.25055381 10.1109/TMI.2014.2340135
23. Dutta J Leahy RM Li Q Non-local means denoising of dynamic PET images PLoS ONE 2013 8 e81390 10.1371/journal.pone.0081390 24339921
Dutta J, Leahy RM, Li Q. Non-local means denoising of dynamic PET images. PLoS ONE. 2013;8: e81390. 10.1371/journal.pone.0081390.24339921 10.1371/journal.pone.0081390
24. Graham MM Muzi M Spence AM O'Sullivan F Lewellen TK Link JM The FDG lumped constant in normal human brain J Nucl Med 2002 43 1157 1166 12215553
Graham MM, Muzi M, Spence AM, O’Sullivan F, Lewellen TK, Link JM, et al. The FDG lumped constant in normal human brain. J Nucl Med. 2002;43:1157–66.12215553
25. Friston KJ Williams S Howard R Frackowiak RS Turner R Movement-related effects in fMRI time-series Magn Reson Med 1996 35 346 355 10.1002/mrm.1910350312 8699946
Friston KJ, Williams S, Howard R, Frackowiak RS, Turner R. Movement-related effects in fMRI time-series. Magn Reson Med. 1996;35:346–55.8699946 10.1002/mrm.1910350312
26. Olszowy W Aston J Rua C Williams GB Accurate autocorrelation modeling substantially improves fMRI reliability Nat Commun 2019 10 1220 10.1038/s41467-019-09230-w 30899012
Olszowy W, Aston J, Rua C, Williams GB. Accurate autocorrelation modeling substantially improves fMRI reliability. Nat Commun. 2019;10:1220. 10.1038/s41467-019-09230-w.30899012 10.1038/s41467-019-09230-w
27. Godbersen GM, Klug S, Wadsak W, Pichler V, Raitanen J, Rieckmann A, et al. Task-evoked metabolic demands of the posteromedial default mode network are shaped by dorsal attention and frontoparietal control networks. eLife. 2023;12:e84683.
28. Power JD Barnes KA Snyder AZ Schlaggar BL Petersen SE Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion Neuroimage 2012 59 2142 2154 10.1016/j.neuroimage.2011.10.018 22019881
Power JD, Barnes KA, Snyder AZ, Schlaggar BL, Petersen SE. Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage. 2012;59:2142–54. 10.1016/j.neuroimage.2011.10.018.22019881 10.1016/j.neuroimage.2011.10.018
29. Duarte JM Gruetter R Characterization of cerebral glucose dynamics in vivo with a four-state conformational model of transport at the blood-brain barrier J Neurochem 2012 121 396 406 10.1111/j.1471-4159.2012.07688.x 22324542
Duarte JM, Gruetter R. Characterization of cerebral glucose dynamics in vivo with a four-state conformational model of transport at the blood-brain barrier. J Neurochem. 2012;121:396–406. 10.1111/j.1471-4159.2012.07688.x.22324542 10.1111/j.1471-4159.2012.07688.x
30. Gruetter R Ugurbil K Seaquist ER Steady-state cerebral glucose concentrations and transport in the human brain J Neurochem 1998 70 397 408 10.1046/j.1471-4159.1998.70010397.x 9422387
Gruetter R, Ugurbil K, Seaquist ER. Steady-state cerebral glucose concentrations and transport in the human brain. J Neurochem. 1998;70:397–408. 10.1046/j.1471-4159.1998.70010397.x.9422387 10.1046/j.1471-4159.1998.70010397.x
31. Cunningham VJ Hargreaves RJ Pelling D Moorhouse SR Regional blood-brain glucose transfer in the rat: a novel double-membrane kinetic analysis J Cereb Blood Flow Metab 1986 6 305 314 10.1038/jcbfm.1986.53 3711158
Cunningham VJ, Hargreaves RJ, Pelling D, Moorhouse SR. Regional blood-brain glucose transfer in the rat: a novel double-membrane kinetic analysis. J Cereb Blood Flow Metab. 1986;6:305–14. 10.1038/jcbfm.1986.53.3711158 10.1038/jcbfm.1986.53
32. Barros LF Bittner CX Loaiza A Porras OH A quantitative overview of glucose dynamics in the gliovascular unit Glia 2007 55 1222 1237 10.1002/glia.20375 17659523
Barros LF, Bittner CX, Loaiza A, Porras OH. A quantitative overview of glucose dynamics in the gliovascular unit. Glia. 2007;55:1222–37. 10.1002/glia.20375.17659523 10.1002/glia.20375
33. Jamadar SD, Liang EX, Zhong S, Ward PGD, Carey A, McIntyre R, et al. Monash DaCRA fPET-fMRI: A dataset for comparison of radiotracer administration for high temporal resolution functional FDG-PET. GigaScience. 2022;11. 10.1093/gigascience/giac031.
34. van Sluis J de Jong J Schaar J Noordzij W van Snick P Dierckx R Performance Characteristics of the Digital Biograph Vision PET/CT System J Nucl Med 2019 60 1031 1036 10.2967/jnumed.118.215418 30630944
van Sluis J, de Jong J, Schaar J, Noordzij W, van Snick P, Dierckx R, et al. Performance Characteristics of the Digital Biograph Vision PET/CT System. J Nucl Med. 2019;60:1031–6. 10.2967/jnumed.118.215418.30630944 10.2967/jnumed.118.215418
35. Delso G Furst S Jakoby B Ladebeck R Ganter C Nekolla SG Performance measurements of the Siemens mMR integrated whole-body PET/MR scanner J Nucl Med 2011 52 1914 1922 10.2967/jnumed.111.092726 22080447
Delso G, Furst S, Jakoby B, Ladebeck R, Ganter C, Nekolla SG, et al. Performance measurements of the Siemens mMR integrated whole-body PET/MR scanner. J Nucl Med. 2011;52:1914–22. 10.2967/jnumed.111.092726.22080447 10.2967/jnumed.111.092726
36. Vandenberghe S Mikhaylova E D'Hoe E Mollet P Karp JS Recent developments in time-of-flight PET EJNMMI Phys 2016 3 3 10.1186/s40658-016-0138-3 26879863
Vandenberghe S, Mikhaylova E, D’Hoe E, Mollet P, Karp JS. Recent developments in time-of-flight PET. EJNMMI Phys. 2016;3:3. 10.1186/s40658-016-0138-3.26879863 10.1186/s40658-016-0138-3
37. Mann K Deny S Ganguli S Clandinin TR Coupling of activity, metabolism and behaviour across the Drosophila brain Nature 2021 593 244 248 10.1038/s41586-021-03497-0 33911283
Mann K, Deny S, Ganguli S, Clandinin TR. Coupling of activity, metabolism and behaviour across the Drosophila brain. Nature. 2021;593:244–8. 10.1038/s41586-021-03497-0.33911283 10.1038/s41586-021-03497-0
38. Li B Freeman RD Neurometabolic coupling between neural activity, glucose, and lactate in activated visual cortex J Neurochem 2015 135 742 754 10.1111/jnc.13143 25930947
Li B, Freeman RD. Neurometabolic coupling between neural activity, glucose, and lactate in activated visual cortex. J Neurochem. 2015;135:742–54. 10.1111/jnc.13143.25930947 10.1111/jnc.13143
39. Roche R Salazar P Martin M Marcano F Gonzalez-Mora JL Simultaneous measurements of glucose, oxyhemoglobin and deoxyhemoglobin in exposed rat cortex J Neurosci Methods 2011 202 192 198 10.1016/j.jneumeth.2011.07.003 21801749
Roche R, Salazar P, Martin M, Marcano F, Gonzalez-Mora JL. Simultaneous measurements of glucose, oxyhemoglobin and deoxyhemoglobin in exposed rat cortex. J Neurosci Methods. 2011;202:192–8. 10.1016/j.jneumeth.2011.07.003.21801749 10.1016/j.jneumeth.2011.07.003
40. Huang SC Phelps ME Hoffman EJ Kuhl DE Error sensitivity of fluorodeoxyglucose method for measurement of cerebral metabolic rate of glucose J Cereb Blood Flow Metab 1981 1 391 401 10.1038/jcbfm.1981.43 6976976
Huang SC, Phelps ME, Hoffman EJ, Kuhl DE. Error sensitivity of fluorodeoxyglucose method for measurement of cerebral metabolic rate of glucose. J Cereb Blood Flow Metab. 1981;1:391–401. 10.1038/jcbfm.1981.43.6976976 10.1038/jcbfm.1981.43
41. Fox PT Raichle ME Focal physiological uncoupling of cerebral blood flow and oxidative metabolism during somatosensory stimulation in human subjects Proc Natl Acad Sci U S A 1986 83 1140 1144 10.1073/pnas.83.4.1140 3485282
Fox PT, Raichle ME. Focal physiological uncoupling of cerebral blood flow and oxidative metabolism during somatosensory stimulation in human subjects. Proc Natl Acad Sci U S A. 1986;83:1140–4. 10.1073/pnas.83.4.1140.3485282 10.1073/pnas.83.4.1140
42. Raichle ME Behind the scenes of functional brain imaging: a historical and physiological perspective Proc Natl Acad Sci U S A 1998 95 765 772 10.1073/pnas.95.3.765 9448239
Raichle ME. Behind the scenes of functional brain imaging: a historical and physiological perspective. Proc Natl Acad Sci U S A. 1998;95:765–72. 10.1073/pnas.95.3.765.9448239 10.1073/pnas.95.3.765
43. Hargreaves RJ Planas AM Cremer JE Cunningham VJ Studies on the relationship between cerebral glucose transport and phosphorylation using 2-deoxyglucose J Cereb Blood Flow Metab 1986 6 708 716 10.1038/jcbfm.1986.127 3793806
Hargreaves RJ, Planas AM, Cremer JE, Cunningham VJ. Studies on the relationship between cerebral glucose transport and phosphorylation using 2-deoxyglucose. J Cereb Blood Flow Metab. 1986;6:708–16. 10.1038/jcbfm.1986.127.3793806 10.1038/jcbfm.1986.127
44. Leybaert L Neurobarrier coupling in the brain: a partner of neurovascular and neurometabolic coupling? J Cereb Blood Flow Metab 2005 25 2 16 10.1038/sj.jcbfm.9600001 15678108
Leybaert L. Neurobarrier coupling in the brain: a partner of neurovascular and neurometabolic coupling? J Cereb Blood Flow Metab. 2005;25:2–16. 10.1038/sj.jcbfm.9600001.15678108 10.1038/sj.jcbfm.9600001
45. Rothman DL Dienel GA Behar KL Hyder F DiNuzzo M Giove F Glucose sparing by glycogenolysis (GSG) determines the relationship between brain metabolism and neurotransmission J Cereb Blood Flow Metab 2022 42 844 860 10.1177/0271678X211064399 34994222
Rothman DL, Dienel GA, Behar KL, Hyder F, DiNuzzo M, Giove F, et al. Glucose sparing by glycogenolysis (GSG) determines the relationship between brain metabolism and neurotransmission. J Cereb Blood Flow Metab. 2022;42:844–60. 10.1177/0271678X211064399.34994222 10.1177/0271678X211064399
46. Attwell D Buchan AM Charpak S Lauritzen M Macvicar BA Newman EA Glial and neuronal control of brain blood flow Nature 2010 468 232 243 10.1038/nature09613nature09613[pii] 21068832
Attwell D, Buchan AM, Charpak S, Lauritzen M, Macvicar BA, Newman EA. Glial and neuronal control of brain blood flow. Nature. 2010;468:232–43. 10.1038/nature09613nature09613[pii].21068832 10.1038/nature09613nature09613[pii]
47. Mishra A Reynolds JP Chen Y Gourine AV Rusakov DA Attwell D Astrocytes mediate neurovascular signaling to capillary pericytes but not to arterioles Nat Neurosci 2016 19 1619 1627 10.1038/nn.4428 27775719
Mishra A, Reynolds JP, Chen Y, Gourine AV, Rusakov DA, Attwell D. Astrocytes mediate neurovascular signaling to capillary pericytes but not to arterioles. Nat Neurosci. 2016;19:1619–27. 10.1038/nn.4428.27775719 10.1038/nn.4428
48. Harris JJ Jolivet R Attwell D Synaptic energy use and supply Neuron 2012 75 762 777 10.1016/j.neuron.2012.08.019 22958818
Harris JJ, Jolivet R, Attwell D. Synaptic energy use and supply. Neuron. 2012;75:762–77. 10.1016/j.neuron.2012.08.019.22958818 10.1016/j.neuron.2012.08.019
49. Magistretti PJ Allaman I A cellular perspective on brain energy metabolism and functional imaging Neuron 2015 86 883 901 10.1016/j.neuron.2015.03.035 25996133
Magistretti PJ, Allaman I. A cellular perspective on brain energy metabolism and functional imaging. Neuron. 2015;86:883–901. 10.1016/j.neuron.2015.03.035.25996133 10.1016/j.neuron.2015.03.035
50. Raichle ME Mintun MA Brain work and brain imaging Annu Rev Neurosci 2006 29 449 476 10.1146/annurev.neuro.29.051605.112819 16776593
Raichle ME, Mintun MA. Brain work and brain imaging. Annu Rev Neurosci. 2006;29:449–76. 10.1146/annurev.neuro.29.051605.112819.16776593 10.1146/annurev.neuro.29.051605.112819
51. Buxton RB Interpreting oxygenation-based neuroimaging signals: the importance and the challenge of understanding brain oxygen metabolism Frontiers in neuroenergetics 2010 2 8 10.3389/fnene.2010.00008 20616882
Buxton RB. Interpreting oxygenation-based neuroimaging signals: the importance and the challenge of understanding brain oxygen metabolism. Frontiers in neuroenergetics. 2010;2:8. 10.3389/fnene.2010.00008.20616882 10.3389/fnene.2010.00008
52. Handwerker DA Ollinger JM D'Esposito M Variation of BOLD hemodynamic responses across subjects and brain regions and their effects on statistical analyses Neuroimage 2004 21 1639 1651 10.1016/j.neuroimage.2003.11.029 15050587
Handwerker DA, Ollinger JM, D’Esposito M. Variation of BOLD hemodynamic responses across subjects and brain regions and their effects on statistical analyses. Neuroimage. 2004;21:1639–51.15050587 10.1016/j.neuroimage.2003.11.029
53. Schaeffer S Iadecola C Revisiting the neurovascular unit Nat Neurosci 2021 24 1198 1209 10.1038/s41593-021-00904-7 34354283
Schaeffer S, Iadecola C. Revisiting the neurovascular unit. Nat Neurosci. 2021;24:1198–209. 10.1038/s41593-021-00904-7.34354283 10.1038/s41593-021-00904-7
54. Iadecola C The Neurovascular Unit Coming of Age: A Journey through Neurovascular Coupling in Health and Disease Neuron 2017 96 17 42 10.1016/j.neuron.2017.07.030 28957666
Iadecola C. The Neurovascular Unit Coming of Age: A Journey through Neurovascular Coupling in Health and Disease. Neuron. 2017;96:17–42. 10.1016/j.neuron.2017.07.030.28957666 10.1016/j.neuron.2017.07.030
55. Hawkins RA Mans AM Davis DW Hibbard LS Lu DM Glucose availability to individual cerebral structures is correlated to glucose metabolism J Neurochem 1983 40 1013 1018 10.1111/j.1471-4159.1983.tb08086.x 6834033
Hawkins RA, Mans AM, Davis DW, Hibbard LS, Lu DM. Glucose availability to individual cerebral structures is correlated to glucose metabolism. J Neurochem. 1983;40:1013–8. 10.1111/j.1471-4159.1983.tb08086.x.6834033 10.1111/j.1471-4159.1983.tb08086.x
56. Stiernman L Grill F McNulty C Bahrd P Panes Lundmark V Axelsson J Widespread fMRI BOLD Signal Overactivations during Cognitive Control in Older Adults Are Not Matched by Corresponding Increases in fPET Glucose Metabolism J Neurosci 2023 43 2527 2536 10.1523/JNEUROSCI.1331-22.2023 36868855
Stiernman L, Grill F, McNulty C, Bahrd P, Panes Lundmark V, Axelsson J, et al. Widespread fMRI BOLD Signal Overactivations during Cognitive Control in Older Adults Are Not Matched by Corresponding Increases in fPET Glucose Metabolism. J Neurosci. 2023;43:2527–36. 10.1523/JNEUROSCI.1331-22.2023.36868855 10.1523/JNEUROSCI.1331-22.2023
57. Tarantini S Tran CHT Gordon GR Ungvari Z Csiszar A Impaired neurovascular coupling in aging and Alzheimer's disease: Contribution of astrocyte dysfunction and endothelial impairment to cognitive decline Exp Gerontol 2017 94 52 58 10.1016/j.exger.2016.11.004 27845201
Tarantini S, Tran CHT, Gordon GR, Ungvari Z, Csiszar A. Impaired neurovascular coupling in aging and Alzheimer’s disease: Contribution of astrocyte dysfunction and endothelial impairment to cognitive decline. Exp Gerontol. 2017;94:52–8. 10.1016/j.exger.2016.11.004.27845201 10.1016/j.exger.2016.11.004
58. Teipel S Drzezga A Grothe MJ Barthel H Chetelat G Schuff N Multimodal imaging in Alzheimer's disease: validity and usefulness for early detection Lancet Neurol 2015 14 1037 1053 10.1016/S1474-4422(15)00093-9 26318837
Teipel S, Drzezga A, Grothe MJ, Barthel H, Chetelat G, Schuff N, et al. Multimodal imaging in Alzheimer’s disease: validity and usefulness for early detection. Lancet Neurol. 2015;14:1037–53. 10.1016/S1474-4422(15)00093-9.26318837 10.1016/S1474-4422(15)00093-9
59. Vaishnavi SN Vlassenko AG Rundle MM Snyder AZ Mintun MA Raichle ME Regional aerobic glycolysis in the human brain Proc Natl Acad Sci U S A 2010 107 17757 17762 10.1073/pnas.1010459107 20837536
Vaishnavi SN, Vlassenko AG, Rundle MM, Snyder AZ, Mintun MA, Raichle ME. Regional aerobic glycolysis in the human brain. Proc Natl Acad Sci U S A. 2010;107:17757–62.20837536 10.1073/pnas.1010459107
60. Theriault JE Shaffer C Dienel GA Sander CY Hooker JM Dickerson BC A functional account of stimulation-based aerobic glycolysis and its role in interpreting BOLD signal intensity increases in neuroimaging experiments PsyArXiv 2021 10.31234/osf.io/pkzr8
Theriault JE, Shaffer C, Dienel GA, Sander CY, Hooker JM, Dickerson BC, et al. A functional account of stimulation-based aerobic glycolysis and its role in interpreting BOLD signal intensity increases in neuroimaging experiments. PsyArXiv. 2021. 10.1016/j.neubiorev.2023.105373.10.31234/osf.io/pkzr8
