- 1Institute of Neuroscience and Medicine – 4, Forschungszentrum Jülich, Jülich, Germany
- 2Department of Nuclear Medicine, Medical Faculty RWTH Aachen, Aachen, Germany
- 3Department of Nuclear medicine, University Hospital Halle, Halle (Saale), Germany
- 4Department of Radiology and Nuclear Medicine, Maastricht University Medical Centre, Maastricht, Netherlands
- 5Institute of Neuroscience and Medicine – 5, Forschungszentrum Jülich, Jülich, Germany
- 6Mind, Brain Imaging and Neuroethics Research Unit, Institute of Mental Health, Royal Ottawa Mental Health Centre and University of Ottawa, Ottawa, ON, Canada
- 7Institute of Neuroscience and Medicine – 11, Forschungszentrum Jülich, Jülich, Germany
- 8JARA – BRAIN – Translational Medicine, Aachen, Germany
- 9Department of Neurology, RWTH Aachen University, Aachen, Germany
Introduction: The fMRI BOLD response to emotional stimuli highlighting the role of the medial prefrontal cortex (MPFC) has been thoroughly investigated. Recently, the relationship between emotion processing and GABA levels has been studied using MPFC proton magnetic resonance spectroscopy (1H-MRS). However, the role of GABAA receptors in the MPFC during emotion processing remains unexplored.
Methods: Using [11C]flumazenil PET, we investigated the relationship between the binding potential of GABAA receptors and emotion processing as measured using simultaneous fMRI BOLD. We hypothesized a correlation between the percent signal change in the BOLD signal and the binding potential of GABAA receptors in the MPFC. In a combined simultaneous fMRI and [11C]flumazenil-PET study, we analyzed the data from 15 healthy subjects using visual emotional stimuli. Our task comprised two types of emotional processing: passive viewing and appraisal. Following the administration of a bolus plus infusion protocol, PET and fMRI data were simultaneously acquired in a hybrid 3 T MR-BrainPET.
Results: We found a differential correlation of BOLD percent signal change with [11C]flumazenil binding potential in the MPFC. Specifically, [11C]flumazenil binding potential in the ventromedial prefrontal cortex (vMPFC) correlated with passive viewing of emotionally valenced pictures. In contrast, the [11C]flumazenil binding potential and the BOLD signal induced by picture appraisal did show a correlation in the paracingulate gyrus.
Conclusion: Our data deliver first evidence for a relationship between MPFC GABAA receptors and emotion processing in the same region. Moreover, we observed that GABAA receptors appear to play different roles in emotion processing in the vMPFC (passive viewing) and paracingulate gyrus (appraisal).
1. Introduction
The basis of any brain function results from neural signaling and the associated interplay of neurotransmitters that are organized in large-scale networks. These networks cause fluctuations in the cerebral metabolic rate of oxygen consumption, inducing changes in the blood content of deoxyhaemoglobin, as well as changes in blood flow and volume, representing the basis of the fMRI BOLD signal (Lauritzen et al., 2012). One of the biggest challenges in neuroimaging is to understand how these processes (i.e., neurochemistry and neurophysiology) are related to each other on various functional levels.
Gamma-aminobutyric acid (GABA) is considered the principal inhibitory neurotransmitter in the central nervous system (CNS) (de Leon and Tadi, 2021), and many studies have addressed its influence on neural activity and brain hemodynamic oxygen response on a microscopic level (Lauritzen et al., 2012; Kiemes et al., 2021). However, in order to address this relationship at a macroscopic level, different imaging modalities, such as functional magnetic resonance imaging (fMRI), proton magnetic resonance spectroscopy (1H-MRS) and positron emission tomography (PET), have been employed in multimodal imaging studies (Duncan et al., 2014a). To date, studies combining 1H-MRS and fMRI have delivered evidence of associations between local GABA levels and fMRI activation in visual processing in the occipital lobe (Kiemes et al., 2021). Moreover, they have found a similar relationship between local MRS-based GABA levels and fMRI activation and emotional processing in the MPFC, and, more predominately, in the anterior cingulate (ACC). Many fMRI studies have demonstrated the role of the MPFC in emotional processing, in particular, in the processing of emotional valences (the emotional dimension referring to the extent to which an emotion is positive or negative) (Phan et al., 2002; Heinzel et al., 2005; Grimm et al., 2006; Willinger et al., 2019). Recently, the processing of different emotional valences has been linked to GABA levels in the ACC (Levar et al., 2017). However, it remains unclear how the local GABA levels are related to the anatomical distribution of GABA receptors.
Receptor imaging with [11C]flumazenil using positron emission tomography (PET) is a powerful technique used to probe the binding potential of GABAA receptors in humans non-invasively (Murrell et al., 2020). However, the use of combined fMRI and [11C]flumazenil PET in human studies remains rather limited. Furthermore, the existing literature is predominantly based on the less sophisticated technique of serial measurements, in contrast to the more advanced approach employed in this study. Applying aversive electrical stimuli, Hayes et al. (2013) showed that individual differences in sensorimotor cortex signal changes could be predicted by GABAA receptors in the MPFC and sensorimotor cortex. Other combined studies have found that GABAA receptors are involved in differential processing of internally and externally guided awareness (Wiebking et al., 2012) as well as in neural variability modulation (Qin et al., 2018).
In this study, we investigated the relationship of the GABAA receptor binding potential with the percent signal change of the fMRI BOLD signal during emotional valence processing using a unique approach to obtain simultaneous fMRI and [11C]flumazenil-PET measurements in healthy participants. We hypothesized that changes in the BOLD percent signal would correlate with the binding potential of GABAA receptors in the MPFC. In order to further investigate the specific type of emotional valence processing, we used different tasks requiring either passive viewing or emotional appraisal. In addition to our hypothesis-driven approach regarding the MPFC, we additionally examined other regions in an exploratory whole-brain approach.
2. Materials and methods
2.1. Participants
Twenty healthy volunteers participated in the study. None of the participants had a history of neurological or psychiatric disease [assessed by the Mini-International Neuropsychiatry Interview (MINI)] (Sheehan et al., n.d.), and all participants completed the German version of the Beck Depression Inventory (BDI) (Kühner et al., 2007), the German version of the Toronto Alexithymia Scale (Bach et al., 1996), and the Positive and Negative Affect Schedule (PANAS) (Watson et al., 1988). None of the participants were using medication at the time of scanning. Standard MRI exclusion criteria were applied.
All measurements were approved by the local ethics committee/federal authorities and were conducted in accordance with the declaration of Helsinki. All participants gave written informed consent prior to the measurement.
2.2. MR and PET data acquisition and reconstruction
A detailed description of the synthesis and administration of [11C]flumazenil using a bolus plus infusion scheme, along with the details of data acquisition and reconstruction, can be found in (Mauler et al., 2020). Briefly, the measurements were performed using a Siemens (Erlangen, Germany) hybrid 3 T MR-BrainPET scanner at the Forschungszentrum Jülich, Germany (Herzog et al., 2011). Following the administration of the [11C]flumazenil bolus and the initiation of the infusion pump, list-mode data were acquired for 100 min. MR data were acquired concurrently with PET data.
Shortly after a bolus injection with a duration of 3 s, structural, T1-weighted MPRAGE (Magnetization Prepared RApid Gradient Echo) MR data were measured (repetition time TR = 2.25 s, echo time TE = 3.03 ms, in-plane field of view FOV = 256 × 192, resolution 1 mm isotropic, 256 sagittal slices acquired with zero gaps). After 60 min, when the beginning of the steady-state interval was reached, functional scans were performed with an echo-planar imaging (EPI) sequence for 50 min (TR = 3 s, TE = 30 ms, matrix size 64 × 64, in-plane resolution 3 mm, 36 slices in ascending order, slice thickness 3 mm, and inter-slice gap 0.75 mm).
The PET images were iteratively reconstructed using PRESTO (120 iterations) (Scheins et al., 2011, 2015) and a framing of 20 × 300 s. Corrections were included to obtain quantitative images, i.e., for normalization, dead-time, attenuation, scattered and random events (Mauler et al., 2020). In addition, motion correction and registration to the T1-weighted MPRAGE data using the Pmod software package (Version 3.5) was also performed.
2.3. Paradigm for fMRI
Subjects were presented with the fMRI paradigm during the time interval in which the [11C]flumazenil steady state was reached (Figure 1A). Since the combined analysis of fMRI and PET data did not focus on binding potential changes at distinct fMRI event times, the fMRI measurement was continued for 10 min after the end of the PET measurement to gain a higher SNR.
Figure 1. (A) Schematic of the experimental design. (B) Schematic depiction of a single trial out of the sequence of 114 trials of the fMRI paradigm. The stimulus picture was presented for 3 s, and followed by a fixation cross with variable duration (3–5 s). The stimulus was either shown again with a question referring to the type of judgment and terminated by the fixation cross for 10–14 s or a new picture without question was shown.
The viewing task comprised two types of emotional processing: passive viewing of emotional stimuli and cognitive appraisal of these stimuli with regard to emotional valence (i.e., positive or negative) or self-relevance (i.e., high or low personal relevance of the picture content) (Northoff et al., 2004; Heinzel et al., 2005; Northoff, 2016). Accordingly, participants were asked to first passively view and then to judge visual stimuli from the International Affective Picture System (IAPS) (Lang et al., 2008). Since the focus of this paper is valence processing, the selected IAPS images were matched for arousal levels according to the values stated in Libkuman et al. (2007). Pictures with normative valence ratings ranging between four and six were categorized as neutral, while those with ratings lower than four were considered negative, and those with ratings above six were treated as positive (Grühn and Scheibe, 2008; Gong and Wang, 2016). To this end, 64 emotional (38 negative) stimuli and 50 neutral stimuli were selected for inclusion, with mean arousal values of 3.70 for neutral stimuli and 3.77 for emotional stimuli (negative 3.68; positive 3.90). The mean valence values of the negative stimuli were 2.53, 5.04 for neutral stimuli, and 6.64 for positive stimuli. A schematic representation of the fMRI task is shown in Figure 1B. The IAPS pictures were presented for 3 s (passive viewing), followed by a fixation cross with a variable duration (3.0, 3.5, 4.0, 4.5, and 5.0). Subsequently, the same IAPS picture was presented for a second time, but this time with written instructions indicating the type of judgment (i.e., emotional valence or self-relevance each). The trial was terminated by the presentation of the fixation cross with a variable duration (10.0, 10.5, 11.0, 11.5, 12.0, 12.5, 13.0, 13.5, and 14.0). Altogether, 94 of these trials were presented (47 for each type of judgment) to each participant. In addition, in 20 trials, the passive viewing task was not followed by a judgment task but by another passive viewing task with a different IAPS picture. The judgment was performed by pressing a button on a response device (LUMItouch™, Photon Control Inc., BC, Canada) with the index finger of the left or right hand. Reaction times were recorded. The fixation cross served as the baseline condition (Newman et al., 2001; Northoff et al., 2004). The different types of IAPS passive viewing and judgment tasks were presented in a randomized order. Before the experimental session, participants familiarized themselves with the paradigm using a test run with different stimuli.
The pictures were generated using the stimulus presentation and response recording tools from the Presentation Version 14.5 (Neurobehavioral Systems, Inc., www.neurobs.com/) software package. During the fMRI measurements, the pictures were presented on a monitor inside the scanner room, and participants viewed the monitor through a mirror positioned on the head coil in a supine position.
The fMRI task-based functional data were preprocessed and analyzed using statistical parametric mapping 12 (SPM12, Wellcome Department of Imaging Neuroscience, London, https://www.fil.ion.ucl.ac.uk/spm). In order to permit magnetic field saturation, the first three volumes of each condition were discarded. The remaining images were realigned to the mean image, slice-time corrected, and co-registered using the individual 3D-MPRAGE images. All images were normalized to the Montreal Neurological Institute (MNI) space, resampled to the resolution of 2 × 2 × 2 mm, and spatially smoothed using a 5 mm full-width-at-half-maximum Gaussian kernel. In order to obtain the signal changes of the hemodynamic response for the subsequent correlation analysis, intra-individual first-level analyses using the general linear model (Friston et al., 1995) were conducted throughout. The design matrix included regressors encoding neutral, negative, and positive IAPS stimuli, as well as for valence judgment and self-relevance judgment. After defining the first-level models, subject-specific activations were calculated. The contrasts were defined such that they corresponded to percent signal change courses in relation to the following events: passive viewing against baseline for negative, positive, neutral and all pictures; picture appraisal against baseline for negative, positive, neutral and all pictures. Voxels showing possible motor activation upon pressing the buttons on the response devices were excluded from the subsequent combined fMRI and PET data analysis. For this purpose, the two-sample t-tests “picture appraisal minus passive viewing” were conducted for every valence class and converted into masks.
The results were then subjected to the combined fMRI and PET data analysis.
2.3.1. Combined fMRI and PET data analysis
The flumazenil data were registered to the MNI space (iso 2 mm) using SPM12, and the uptake C in every voxel was expressed as the non-displaceable binding potential BPND = (C-Cref)/Cref with the activity Cref from the pure white matter compartment of the pons as the reference region. In order to exclude any (potential) effects from the analysis that were unrelated to the experiment, the binding potential was calculated for exactly those timepoints (i.e., for the corresponding PET frames) at which both conditions were fulfilled: the emotional stimulation paradigm was presented and the steady state was reached (SD from mean calculated over the stimulation interval was less than 10%). In addition, the binding potential was calculated for the pre-stimulation frames – however, it was only calculated for those frames in which the signal was within the same mean ± SD range determined for the stimulation interval. This allowed the comparison of the binding potentials before and during emotional stimulation. The PET signal from the reference region, the pons, was exported and numerically decomposed (Mauler et al., 2020) into the contributions from white and grey matter. The uptake in the pure white matter compartment was used as the reference region uptake, Cref. The unwanted contribution from uptake in neuronal nuclei was captured in the value of the pure grey matter compartment and discarded. In order to level out any remaining small signal fluctuations, Cref and C were calculated as the mean values over the included frames. This, and the following calculations, were performed by an in-house developed Matlab program.
Since the definition of the vMPFC varies to some extent between different authors, anatomical VOI masks were defined according to the Harvard-Oxford Cortical Atlas (Desikan et al., 2006). In addition, in order to ensure comparability with other studies, the analysis of the vMPFC was calculated in terms of the regions defined by this atlas (Figure 2). To this end, the vMPFC was defined to include the ventral medial frontal pole, frontal medial cortex (FMC), ventral paracingulate gyrus, ventral anterior cingulate gyrus, and the subcallosal cortex, limited dorsally by the plane z = +10, and laterally by the planes x = ±20 (MNI152 coordinates) (Santos et al., 2011).
Figure 2. Delineation of the ventromedial prefrontal cortex (vMPFC) composed by regions of the Harvard Oxford Cortical Structures Atlas. It is projected on the MNI iso 2 mm brain. The green cross indicates the Frontal Medial Cortex as part of vMPFC.
To minimize the influence of partial volume effects on the analysis, the partial volume of the grey matter compartment was assessed voxel-wisely by means of the tissue probability maps, which were computed from the MPRAGE data by using SPM12. To match the resolution of the individual tissue probability maps to the effective PET data resolution, the maps were smoothed with an iso 1 mm Gauss kernel following registration to the iso 2 mm MNI space (Mauler et al., 2020). Due to the partial volume effect, voxels with a grey matter content of less than 60% were excluded from further analysis. In addition, voxels indicating potential motor activations were excluded using the masks previously described. The percentage signal change P = event signal(s)/mean signal * 100 in the BOLD signal and the flumazenil binding potential values of all remaining voxels from the previously described VOIs were exported. The values from both modalities were subject-individually z-transformed to compensate for the large inter-individual signal differences, VOI-wise averaged and analyzed for the Pearson correlation coefficient between percent signal change and flumazenil binding potential.
3. Results
3.1. Overview
Twenty participants were measured, of which five subjects were excluded because the acquisition of one modality failed, the measurement was interrupted, or the quality of the steady state was insufficient, leaving data from 15 volunteers (five male) with an average age of 33 ± 13 years [mean ± standard deviation (SD)] for inclusion in the analysis. On average, 383 ± 31 MBq [11C]flumazenil was administered per participant using the bolus plus infusion protocol with a mean kbol value of 67 min [53–100 min (minimum – maximum)].
3.2. Behavioral results
For the BDI, the mean score was 2.33, with a SD of 2.38. For the TAS-20, the mean score was 48.73, with a SD of 7.01. The mean PANAS score of negative affect was 1.41 (SD 0.4), and the mean score of positive affect was 3.45 (SD 0.51). The reaction times exceeded the maximum time interval of 3 s in 2.9 ± 1.9 out of 94 trials (1–7 trials).
3.3. fMRI-PET relationship
The [11C]flumazenil binding potential used for the correlation analysis was calculated as the mean of the uptake values over a time interval of 38 ± 8 min (corresponding to 7.6 ± 1.6 frames), beginning at 64 ± 10 min after the bolus injection. The [11C]flumazenil binding potential before stimulation was compared to the binding potential upon emotional stimulation and showed no difference. The binding potential before stimulation was calculated from a time interval of 20 ± 12 min (4.1 ± 2.5 frames). The following describes the results of the correlation analysis between the static property of GABA receptor availability (i.e., baseline flumazenil BPND) and the dynamic percent signal change in the BOLD signal as performed on z-transformed data.
3.3.1. Passive viewing
The results of the correlation analysis are shown in Table 1A and Figures 3A,B and, in line with our a priori hypothesis, indicate that percent signal changes in the BOLD signal in the FMC correlate moderately positively with [11C]flumazenil binding potential. The underlying absolute BPND values (i.e., without z-transformation) in the FMC were 4.34 ± 1.8 (1.51–6.32) (the individual absolute BPND and BOLD percent signal change values for all subjects are listed in Supplementary Table S1). However, while this was found to be true for passive viewing of emotionally valenced pictures (negative and/or positive), it was not the case for neutral pictures. Considering the vMPFC as a whole, similar results were found for passive viewing of emotionally valenced pictures, and significant correlations (p < 0.05) were noted for positive pictures. Thus, the vMPFC, and the FMC in particular, showed a specific moderate correlation between the percent signal change and the [11C]flumazenil binding potential while passively viewing emotional stimuli.
Table 1. Results from the correlation analyses for the contrasts related to passive viewing of visual stimuli in part A and picture appraisal (i.e., emotional judgment) in part B.
Figure 3. Results after stimulation with the event class “passive viewing.” (A) Scatterplots illustrating the VOI-based correlation analyses of the flumazenil binding potential after z-transformation and the z-transformed percent signal change (BOLD) for different types of passive viewing of emotional stimuli (Pearson correlation coefficient and the corresponding p-value) in the frontal medial cortex. (B) VOI-specific correlation coefficients of the flumazenil binding potential and the percentage signal change in the BOLD signal, both after z-transformation, which were mapped to the planum temporale (r = −0.57, p = 0.03) and the frontal medial cortex (r = 0.62, p = 0.01) of the Harvard Oxford Cortical Structures Atlas. Positive stimuli.
The whole-brain analysis revealed other correlations that were not included in our a priori hypothesis. These regions comprise the insular cortex, planum temporale, Heschl’s gyrus, and the planum polare. In contrast to the MPFC, these regions all showed negative (rather than positive) correlations in the BOLD signal changes and the [11C]flumazenil binding potential (Table 1). The highest absolute correlation coefficients were found for the insular cortex (passive viewing of positive stimuli) and the planum temporale (combined analysis of passive viewing of negative, positive, and neutral stimuli) (Table 1A). The most correlations were found for the passive viewing of positive stimuli.
3.3.2. Picture appraisal
The results from the correlation analysis shown in Table 1B and Figures 4A,B demonstrate that changes in the BOLD signal correlate significantly (p < 0.05) with [11C]flumazenil binding potential in the paracingulate gyrus. The correlations were found to be significant for the combined analysis of picture appraisal of negative, positive, and neutral stimuli, as well as separate analyses of neutral stimuli. The largest moderately sized correlation was found for neutral stimuli. No correlation was found in the vMPFC. The underlying absolute BPND values (i.e., without z-transformation) in the paracingulate gyrus were 3.99 ± 1.58 (1.4–5.68) (the individual absolute BPND and BOLD percent signal change values for all subjects are listed in Supplementary Table S2).
Figure 4. Results after stimulation with the event class “picture appraisal.” (A) Scatterplots of the paracingulate gyrus illustrating the Pearson correlation coefficients and the corresponding value of ps of the flumazenil binding potential and the percent signal change (BOLD), both after z-transform, for different appraisals of the stimuli. (B) Representation of the VOIs of the Harvard Oxford Cortical Structures Atlas. The color code represents the correlation between the z-transformed flumazenil binding potential and the z-transformed percentage signal change in the BOLD signal in the paracingulate gyrus (r = 0.60, p = 0.019) and the planum temporale (r = 0.71, p = 0.003). The subjects assessed neutral and negative visual stimuli in terms of self-relatedness and emotional valence.
The whole-brain analysis revealed additional significant correlations (p < 0.05) in the parietal operculum cortex, Heschl’s gyrus, and the planum temporale. In contrast to the results from passive viewing, all regions except the parietal operculum showed positive correlations between the percent change in the BOLD signal and the [11C]flumazenil binding potential (Table 1B).
In order to further characterize the type of judgments, correlations of valence-related and self-related judgment were also investigated. A significant positive correlation with the paracingulate gyrus (r = 0.52, p < 0.05) was found for self-related judgments. The planum temporale showed a positive correlation (r = 0.61, p < 0.05) for the valence-related judgment task. This suggests, albeit tentatively, that the above-mentioned correlations were driven by emotional valence and the self-specificity of our emotional stimuli.
4. Discussion
A differential correlation of the BOLD signal with [11C]flumazenil binding potential was found within the MPFC. Furthermore, the binding potential of [11C]flumazenil in the vMPFC, and particularly in the FMC, correlated moderately with the passive viewing of negatively and positively valenced pictures. No significant correlation was found with passive viewing of neutrally valenced pictures.
In contrast, although the binding potential of [11C]flumazenil in the paracingulate gyrus correlated with the BOLD signal induced by picture appraisal, no such correlation was found in the vMPFC. Moreover, the results indicate that changes in BOLD signal induced by picture appraisal are related to both emotional valence judgment and to self-related judgment.
Various functional imaging studies have demonstrated that the MPFC is commonly activated in emotional processing, including general processes, emotion evaluation, experience, and regulation (Phan et al., 2002; Northoff et al., 2004; Phan et al., 2004; Heinzel et al., 2005; Grimm et al., 2006; Willinger et al., 2019). In addition, studies have shown that the MPFC is specifically involved in the processing of emotional valence (Heinzel et al., 2005; Grimm et al., 2006; Hayes et al., 2014; Yang et al., 2020). The fMRI-PET correlation analysis conducted in this study extends these findings by indicating a specific correlation between visual emotional valence processing and [11C]flumazenil binding potential. To date, only 1H-MRS has been used to investigate the relationship between GABA levels and emotional processing in the MPFC, and meta-analyses investigating GABA levels in the MPFC have focused predominately on the ACC. The results of these studies have been mixed, with some suggesting a negative relationship with neural activity (Duncan et al., 2014b; Kiemes et al., 2021) and others contradicting this view (Kiemes et al., 2021). Although we did not find a negative correlation between the [11C]flumazenil binding potential and the ACC, a moderate positive correlation with the paracingulate gyrus and the FMC was established. However, the role of the experimental setup must be considered when comparing GABA levels resulting from MRS with [11C]flumazenil binding potential, as MRS measures global GABA levels (intra- and extra-cellular) generated from the average tissue concentration of GABA, whereas PET permits the detection of alterations in GABA levels in the extracellular space (Frankle et al., 2015). The experimental design used did not measure rapid changes in BPND under emotional stimulation. Instead, it measured the availability of the GABA receptor and the associated percent signal change in the BOLD signal induced by different types of emotional processing.
In addition, the results from the literature show initial evidence for a negative relationship between GABA levels (measured with 1H-MRS) and [11C]flumazenil binding potential (Persson et al., 2021) as well as with plasma GABA and [11C]flumazenil binding potential (Klumpers et al., 2010). This was interpreted as a form of homeo-static regulation reflected by a regional balance between the transmitter and receptor of the GABA system. Thus, more availability of GABA might be associated with less GABAA receptor availability (i.e., due to downregulation) in the brain (Persson et al., 2021). Therefore, it can be speculated that the positive correlation of GABAA receptor availability and BOLD signal in the vMPFC is associated with less available GABA, as measured with 1H-MRS.
Animal experiment data have shown a positive BOLD response in excitatory neurons, whereas the stimulation of inhibitory neurons resulted in a biphasic response at the stimulation site and a negative BOLD response at downstream sites (Moon et al., 2021). The presented VOI-based analysis captures the dominating correlation but does not differentiate between the spatially intermingled excitatory, inhibitory and downstream neurons in the VOI. This effect may have contributed to the finding of negative correlation coefficients in the planum temporale, Heschl’s gyrus, planum polare and parietal operculum cortex, and the positive correlation coefficients in the regions of the MPFC, which were partially in contrast with literature findings.
Moreover, we identified correlations in two different regions within the MPFC. In contrast to emotional appraisal, which was found to be related to the paracingulate gyrus, significant correlations for passive viewing were located in the vMPFC and, in particular, proved to be moderately large in the FMC. Based on functional imaging studies mainly using fMRI, the vMPFC has been specifically related to emotional experience, which fits in with our results for the task involving passive viewing of emotional stimuli (Heinzel and Northoff, 2009, 2014). As noted above, no correlation was found with regard to the passive viewing of neutral stimuli, indicating that the correlation with [11C]flumazenil binding potential in the vMPFC might be specific to the emotional component of picture viewing. Our findings are further supported by previous studies indicating an association between the paracingulate gyrus and reality monitoring, cognitive functions, and the salience network, which is responsible for determining and assigning saliency to events and stimuli (Fornito et al., 2004; Buda et al., 2011; Garrison et al., 2015; Levar et al., 2019). Picture appraisal, as applied in our experimental paradigm, contains cognitive elements in the context of the evaluation of emotional stimuli with regard to valence or self-relatedness.
In addition to our a priori hypotheses, the whole-brain analysis revealed that the planum temporale and the Heschl’s gyrus also showed correlations with passive viewing and judgment of the stimuli. More specifically, a negative correlation was found for passive viewing, and a positive correlation was found for emotional appraisal in both regions. The planum temporale has been related to various functions, including visual motion perception and visual memory (Antal et al., 2008; Isato et al., 2021), and, based on our results, it is possible to speculate that it has an even more general role in visual processing, including the visual processing of emotional and neutral stimuli. The association of Heschl’s gyrus with inner speech (Hurlburt et al., 2016) may also indicate that it plays a role in decision-making during emotional judgment. However, since these results are not based on our a priori hypotheses, they should be interpreted with caution.
The analysis used the calculation of a modified distribution volume that is determined at tracer equilibrium in order to calculate the [11C]flumazenil binding potential (Lassen, 1992). This method was preferred over modeling with a reference tissue model because it was easy to apply – especially with relatively low SNR and when possible changes in the binding potential upon emotional stimulation could not be completely ruled out. However, for comparison with kinetic modeling, data were fitted with the simplified reference tissue model (Lammertsma and Hume, 1996); based on the assumption that modeling the changes in the BP during stimulation is unnecessary) and the linear simplified reference tissue model (Alpert et al., 2003) (results not shown), both of which were published based on a bolus without infusion experiment. The test of the first model revealed that the underlying system of equations was rank deficient and ill-conditioned. The test of the second model showed that it was not able to fit the steady state. Therefore, using these reference tissue models did not result in more robust outcomes than the steady state distribution volume method.
When evaluating the results of this study, some further limitations should also be considered. Perhaps most notably, GABA levels were not assessed using 1H-MRS. In order to further investigate the relationship between GABA levels, BOLD response, and [11C]flumazenil binding potential in the MPFC, studies should be performed that measure all three parameters. The range of raw BPND values includes the range of mean values found for [11C]flumazenil by other groups (Koeppe et al., 1991; Abadie et al., 1992; Klumpers et al., 2008; Miederer et al., 2009). In contrast to the results from the literature, the binding potential values show a larger standard deviation. These comparisons are only for orientation since the VOIs differ to some extent in position and size. Moreover, other groups frequently carried through kinetic modeling on data acquired immediately after the bolus injection when the SNR is high. Our analysis is based on data that were acquired when the steady state was reached after three half-lives of 11C. This may have contributed to the expanded BP range. The literature (Feng et al., 2016) shows elevated standard deviations of the steady state distribution volumes, which results in an increased propagated scatter of BP values. Relatively small deviations in the denominator of the distribution volume ratio cause a large deviation in the result. A possible approach to estimate BPND with higher accuracy in future studies could be to perform arterial blood sampling so that kinetic modeling with an arterial input function can be applied.
The applied event-related design contained a presentation of the stimuli for 3 s, followed by a relatively long period devoid of stimulus presentation, which – if at all – leads to relatively rapid fluctuations in GABA binding. Given the sensitivity and spatial resolution of the PET scanner used (Herzog et al., 2011), rapid changes in the binding potential after single and short stimuli are not measurable. The analysis focuses on the correlation between the static [11C]flumazenil binding potential and the percent signal change in the BOLD signal.
An inherent problem of correlating variables derived from methods with low signal strength, such as fMRI, is that larger sample sizes may be required in order to ensure robust contrast to noise ratio. Large groups of >100 participants have been suggested (Turner et al., 2018). A possible solution to overcome the enormous financial and logistical barriers might be to set up multi-center studies (Duncan et al., 2019).
Finally, our results are based on predefined regions from the Harvard-Oxford atlas. Since this method considers each region as a whole, the exact sizes of the significant clusters inside the predefined regions remain unclear.
5. Conclusion
These results provide initial evidence for a direct relationship between MPFC GABAA receptor availability and emotional processing in the same region. Moreover, differential correlations of the BOLD signal and [11C]flumazenil binding potential were found for active appraisal and passive viewing of emotional stimuli, indicating the role of the vMPFC in passive viewing and the paracingulate gyrus in appraisal.
Data availability statement
The datasets presented in this article are not readily available because the ethical approval of the study does not allow the data to be made available to outsiders. Requests to access the datasets should be directed to j.mauler@fz-juelich.de.
Ethics statement
The studies involving humans were approved by Ethikkommission an der Med. Fakultät der HHU Düsseldorf Gebäude 14.82.01 Moorenstr. 5 D-40225 Düsseldorf Germany. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
AH, GN, HH, and NS designed the study. AH, JM, HH, JS, FB, CL, FM, BN, and K-JL contributed to data collection and data analysis. AH, JM, NS, and GN drafted the manuscript. All authors contributed to the article and approved the submitted version.
Funding
The study was funded from the institute’s budget (Institute of Neuroscience and Medicine - 4, Forschungszentrum Jülich (member of the Helmholtz Association), Jülich, Germany), which is funded by federal and state funds. Open access publication funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – 491111487.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2023.1027697/full#supplementary-material
References
Abadie, P., Baron, J. C., Bisserbe, J. C., Boulenger, J. P., Rioux, P., Travère, J. M., et al. (1992). Central benzodiazepine receptors in human brain: estimation of regional Bmax and KD values with positron emission tomography. Eur. J. Pharmacol. 213, 107–115. doi: 10.1016/0014-2999(92)90239-Z
Alpert, N. M., Badgaiyan, R. D., Livni, E., and Fischman, A. J. (2003). A novel method for noninvasive detection of neuromodulatory changes in specific neurotransmitter systems. NeuroImage 19, 1049–1060. doi: 10.1016/S1053-8119(03)00186-1
Antal, A., Baudewig, J., Paulus, W., and Dechent, P. (2008). The posterior cingulate cortex and planum temporale/parietal operculum are activated by coherent visual motion. Vis. Neurosci. 25, 17–26. doi: 10.1017/s0952523808080024
Bach, M., Bach, D., de Zwaan, M., Serim, M., and Böhmer, F. (1996). Validation of the German version of the 20-item Toronto Alexithymia Scale in normal persons and psychiatric patients. Psychother. Psychosom. Med. Psychol. 46, 23–28.
Buda, M., Fornito, A., Bergström, Z. M., and Simons, J. S. (2011). A specific brain structural basis for individual differences in reality monitoring. J. Neurosci. 31, 14308–14313. doi: 10.1523/jneurosci.3595-11.2011
de Leon, A. S., and Tadi, P. (2021). “Biochemistry, Gamma Aminobutyric Acid” in StatPearls (Treasure Island (FL): StatPearls Publishing)
Desikan, R. S., Ségonne, F., Fischl, B., Quinn, B. T., Dickerson, B. C., Blacker, D., et al. (2006). An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest, Neuroimage. 31, 968–80 doi: 10.1016/j.neuroimage.2006.01.021
Duncan, N. W., Wiebking, C., Muñoz-Torres, Z., and Northoff, G. (2014a). How to investigate neuro-biochemical relationships on a regional level in humans? Methodological considerations for combining functional with biochemical imaging. J. Neurosci. Methods 221, 183–188. doi: 10.1016/j.jneumeth.2013.10.011
Duncan, N. W., Wiebking, C., and Northoff, G. (2014b). Associations of regional GABA and glutamate with intrinsic and extrinsic neural activity in humans—a review of multimodal imaging studies. Neurosci. Biobehav. Rev. 47, 36–52. doi: 10.1016/j.neubiorev.2014.07.016
Duncan, N. W., Zhang, J., Northoff, G., and Weng, X. (2019). Investigating GABA concentrations measured with macromolecule suppressed and unsuppressed MEGA-PRESS MR spectroscopy and their relationship with BOLD responses in the occipital cortex. J. Magn. Reson. Imaging 50, 1285–1294. doi: 10.1002/jmri.26706
Feng, L., Svarer, C., Madsen, K., Ziebell, M., Dyssegaard, A., Ettrup, A., et al. (2016). Design of Infusion Schemes for Neuroreceptor Imaging: Application to [11C]Flumazenil-PET Steady-State Study. Biomed. Res. Int. 2016, 1–8. doi: 10.1155/2016/9132840
Fornito, A., Yücel, M., Wood, S., Stuart, G. W., Buchanan, J. A., Proffitt, T., et al. (2004). Individual differences in anterior cingulate/paracingulate morphology are related to executive functions in healthy males. Cereb. Cortex 14, 424–431. doi: 10.1093/cercor/bhh004
Frankle, W. G., Cho, R. Y., Prasad, K. M., Mason, N. S., Paris, J., Himes, M. L., et al. (2015). In vivo measurement of GABA transmission in healthy subjects and schizophrenia patients. Am. J. Psychiatry 172, 1148–1159. doi: 10.1176/appi.ajp.2015.14081031
Friston, K. J., Frith, C. D., Turner, R., and Frackowiak, R. S. (1995). Characterizing evoked hemodynamics with fMRI. NeuroImage 2, 157–165. doi: 10.1006/nimg.1995.1018
Garrison, J. R., Fernyhough, C., McCarthy-Jones, S., Haggard, M., and Simons, J. S. (2015). Paracingulate sulcus morphology is associated with hallucinations in the human brain. Nat. Commun. 6:8956. doi: 10.1038/ncomms9956
Gong, X., and Wang, D. (2016). Applicability of the international affective picture system in chinese older adults: a validation study. Psych. J. 5, 117–124. doi: 10.1002/pchj.131
Grimm, S., Schmidt, C. F., Bermpohl, F., Heinzel, A., Dahlem, Y., Wyss, M., et al. (2006). Segregated neural representation of distinct emotion dimensions in the prefrontal cortex-an fMRI study. NeuroImage 30, 325–340. doi: 10.1016/j.neuroimage.2005.09.006
Grühn, D., and Scheibe, S. (2008). Age-related differences in valence and arousal ratings of pictures from the International Affective Picture System (IAPS): do ratings become more extreme with age? Behav. Res. Methods 40, 512–521. doi: 10.3758/brm.40.2.512
Hayes, D. J., Duncan, N. W., Wiebking, C., Pietruska, K., Qin, P., Lang, S., et al. (2013). GABAA receptors predict aversion-related brain responses: an fMRI-PET investigation in healthy humans. Neuropsychopharmacology 38, 1438–1450. doi: 10.1038/npp.2013.40
Hayes, D. J., Duncan, N. W., Xu, J., and Northoff, G. (2014). A comparison of neural responses to appetitive and aversive stimuli in humans and other mammals. Neurosci. Biobehav. Rev. 45, 350–368. doi: 10.1016/j.neubiorev.2014.06.018
Heinzel, A., Bermpohl, F., Niese, R., Pfennig, A., Pascual-Leone, A., Schlaug, G., et al. (2005). How do we modulate our emotions? Parametric fMRI reveals cortical midline structures as regions specifically involved in the processing of emotional valences. Brain Res. Cogn. Brain Res. 25, 348–358. doi: 10.1016/j.cogbrainres.2005.06.009
Heinzel, A., and Northoff, G. (2009). Emotional feeling and the orbitomedial prefrontal cortex: theoretical and empirical considerations. Philos. Psychol. 22, 443–464. doi: 10.1080/09515080903153592
Heinzel, A., and Northoff, G. (2014). The Relationship of Self-Relatedness and Emotional Processing. J. Conscious. Stud. 21, 30–48.
Herzog, H., Langen, K. J., Weirich, C., Rota Kops, E., Kaffanke, J., Tellmann, L., et al. (2011). High resolution BrainPET combined with simultaneous MRI. Nuklearmedizin 50, 74–82. doi: 10.3413/Nukmed-0347-10-09
Hurlburt, R. T., Alderson-Day, B., Kühn, S., and Fernyhough, C. (2016). Exploring the Ecological Validity of Thinking on Demand: Neural Correlates of Elicited vs Spontaneously Occurring Inner Speech. PloS one 11, –e0147932. doi: 10.1371/journal.pone.0147932
Isato, A., Suhara, T., and Yamada, M. (2021). Resting-state functional connectivity relates to interindividual variations in positive memory. BioRxiv :2021.2006.2020.448902. doi: 10.1101/2021.06.20.448902
Kiemes, A., Davies, C., Kempton, M. J., Lukow, P. B., Bennallick, C., Stone, J. M., et al. (2021). GABA, Glutamate and Neural Activity: A Systematic Review With Meta-Analysis of Multimodal (1)H-MRS-fMRI Studies. Front. Psych. 12:644315. doi: 10.3389/fpsyt.2021.644315
Klumpers, U. M. H., Veltman, D. J., Boellaard, R., Comans, E. F., Zuketto, C., Yaqub, M., et al. (2008). Comparison of plasma input and reference tissue models for analysing [(11)C]flumazenil studies. J. Cereb. Blood Flow Metab. 28, 579–587. doi: 10.1038/sj.jcbfm.9600554
Klumpers, U. M., Veltman, D. J., Drent, M. L., Boellaard, R., Comans, E. F., Meynen, G., et al. (2010). Reduced parahippocampal and lateral temporal GABAA-[11C]flumazenil binding in major depression: preliminary results. Eur. J. Nucl. Med. Mol. Imaging 37, 565–574. doi: 10.1007/s00259-009-1292-9
Koeppe, R. A., Holthoff, V. A., Frey, K. A., Kilbourn, M. R., and Kuhl, D. E. (1991). Compartmental analysis of [11C]flumazenil kinetics for the estimation of ligand transport rate and receptor distribution using positron emission tomography. J. Cereb. Blood Flow Metab. 11, 735–744. doi: 10.1038/jcbfm.1991.130
Kühner, C., Bürger, C., Keller, F., and Hautzinger, M. (2007). Reliability and validity of the Revised Beck Depression Inventory (BDI-II). Results from German samples. Nervenarzt 78, 651–656. doi: 10.1007/s00115-006-2098-7
Lammertsma, A. A., and Hume, S. P. (1996). Simplified reference tissue model for PET receptor studies. NeuroImage 4, 153–158. doi: 10.1006/nimg.1996.0066
Lang, P. J., Bradley, M. M., and Cuthbert, B. N. (2008). International affective picture system (IAPS): affective ratings of pictures and instruction manual. Technical Report A-8. Gainesville: University of Florida.
Lassen, N. A. (1992). Neuroreceptor quantitation in vivo by the steady-state principle using constant infusion or bolus injection of radioactive tracers. J. Cereb. Blood Flow Metab. 12, 709–716. doi: 10.1038/jcbfm.1992.101
Lauritzen, M., Mathiesen, C., Schaefer, K., and Thomsen, K. J. (2012). Neuronal inhibition and excitation, and the dichotomic control of brain hemodynamic and oxygen responses. NeuroImage 62, 1040–1050. doi: 10.1016/j.neuroimage.2012.01.040
Levar, N., Van Doesum, T. J., Denys, D., and Van Wingen, G. A. (2019). Anterior cingulate GABA and glutamate concentrations are associated with resting-state network connectivity. Sci. Rep. 9:2116. doi: 10.1038/s41598-018-38078-1
Levar, N., van Leeuwen, J. M. C., Denys, D., and van Wingen, G. A. (2017). Divergent influences of anterior cingulate cortex GABA concentrations on the emotion circuitry. NeuroImage 158, 136–144. doi: 10.1016/j.neuroimage.2017.06.055
Libkuman, T. M., Otani, H., Kern, R., Viger, S. G., and Novak, N. (2007). Multidimensional normative ratings for the International Affective Picture System. Behav. Res. Methods 39, 326–334. doi: 10.3758/bf03193164
Mauler, J., Heinzel, A., Matusch, A., Herzog, H., Neuner, I., Scheins, J., et al. (2020). Bolus Infusion Scheme for the Adjustment of Steady State [11C]Flumazenil Levels in the Grey Matter and in the Blood Plasma for Neuroreceptor Imaging. NeuroImage 221:117160. doi: 10.1016/j.neuroimage.2020.117160
Miederer, I., Ziegler, S. I., Liedtke, C., Spilker, M. E., Miederer, M., Sprenger, T., et al. (2009). Kinetic modeling of [11C]flumazenil using data-driven methods. Eur. J. Nucl. Med. Mol. Imaging 36, 659–670. doi: 10.1007/s00259-008-0990-z
Moon, H. S., Jiang, H., Vo, T. T., Jung, W. B., Vazquez, A. L., and Kim, S. G. (2021). Contribution of Excitatory and Inhibitory Neuronal Activity to BOLD fMRI. Cereb. Cortex 31, 4053–4067. doi: 10.1093/cercor/bhab068
Murrell, E., Pham, J. M., Sowa, A. R., Brooks, A. F., Kilbourn, M. R., Scott, P. J. H., et al. (2020). Classics in Neuroimaging: Development of Positron Emission Tomography Tracers for Imaging the GABAergic Pathway. ACS Chem. Neurosci. 11, 2039–2044. doi: 10.1021/acschemneuro.0c00343
Newman, S. D., Twieg, D. B., and Carpenter, P. A. (2001). Baseline conditions and subtractive logic in neuroimaging. Hum. Brain Mapp. 14, 228–235. doi: 10.1002/hbm.1055
Northoff, G. (2016). Is the self a higher-order or fundamental function of the brain? The "basis model of self-specificity" and its encoding by the brain's spontaneous activity. Cogn. Neurosci. 7, 203–222. doi: 10.1080/17588928.2015.1111868
Northoff, G., Heinzel, A., Bermpohl, F., Niese, R., Pfennig, A., Pascual-Leone, A., et al. (2004). Reciprocal modulation and attenuation in the prefrontal cortex: an fMRI study on emotional-cognitive interaction. Hum. Brain Mapp. 21, 202–212. doi: 10.1002/hbm.20002
Persson, J., Wall, A., Weis, J., Gingnell, M., Antoni, G., Lubberink, M., et al. (2021). Inhibitory and excitatory neurotransmitter systems in depressed and healthy: A positron emission tomography and magnetic resonance spectroscopy study. Psychiatry Res. Neuroimaging 315:111327. doi: 10.1016/j.pscychresns.2021.111327
Phan, K. L., Wager, T., Taylor, S. F., and Liberzon, I. (2002). Functional neuroanatomy of emotion: a meta-analysis of emotion activation studies in PET and fMRI. NeuroImage 16, 331–348. doi: 10.1006/nimg.2002.1087
Phan, K. L., Wager, T. D., Taylor, S. F., and Liberzon, I. (2004). Functional neuroimaging studies of human emotions. CNS Spectr. 9, 258–266. doi: 10.1017/s1092852900009196
Qin, P., Duncan, N. W., Chen, D. Y., Chen, C. J., Huang, L. K., Huang, Z., et al. (2018). Vascular-metabolic and GABAergic Inhibitory Correlates of Neural Variability Modulation. A Combined fMRI and PET Study. Neuroscience 379, 142–151. doi: 10.1016/j.neuroscience.2018.02.041
Santos, J. P., Seixas, D., Brandão, S., and Moutinho, L. (2011). Investigating the role of the ventromedial prefrontal cortex in the assessment of brands. Front. Neurosci. 5:77. doi: 10.3389/fnins.2011.00077
Scheins, J. J., Herzog, H., and Shah, N. J. (2011). Fully-3D PET image reconstruction using scanner-independent, adaptive projection data and highly rotation-symmetric voxel assemblies. IEEE Trans. Med. Imaging 30, 879–892. doi: 10.1109/tmi.2011.2109732
Scheins, J. J., Vahedipour, K., Pietrzyk, U., and Shah, N. J. (2015). High performance volume-of-intersection projectors for 3D-PET image reconstruction based on polar symmetries and SIMD vectorisation. Phys. Med. Biol. 60, 9349–9375. doi: 10.1088/0031-9155/60/24/9349
Sheehan, D. V., Lecrubier, Y., Sheehan, K. H., Amorim, P., Janavs, J., Weiller, E., et al. The Mini-International Neuropsychiatric Interview (M.I.N.I.): the development and validation of a structured diagnostic psychiatric interview for DSM-IV and ICD-10, The Journal of Clinical Psychiatry. 20, 22–33.
Turner, B. O., Paul, E. J., Miller, M. B., and Barbey, A. K. (2018). Small sample sizes reduce the replicability of task-based fMRI studies. Commun. Biol. 1:62. doi: 10.1038/s42003-018-0073-z
Watson, D., Clark, L. A., and Tellegen, A. (1988). Development and validation of brief measures of positive and negative affect: the PANAS scales. J. Pers. Soc. Psychol. 54, 1063–1070. doi: 10.1037//0022-3514.54.6.1063
Wiebking, C., Duncan, N. W., Qin, P., Hayes, D. J., Lyttelton, O., Gravel, P., et al. (2012). External awareness and GABA-A multimodal imaging study combining fMRI and [(18) F]flumazenil-PET. Hum. Brain Mapp. 35, 173–184. doi: 10.1002/hbm.22166
Willinger, D., Karipidis, I. I., Beltrani, S., Di Pietro, S. V., Sladky, R., Walitza, S., et al. (2019). Valence-Dependent Coupling of Prefrontal-Amygdala Effective Connectivity during Facial Affect Processing. eNeuro 6:ENEURO.0079-0019.2019. doi: 10.1523/ENEURO.0079-19.2019
Yang, M., Tsai, S. J., and Li, C. R. (2020). Concurrent amygdalar and ventromedial prefrontal cortical responses during emotion processing: a meta-analysis of the effects of valence of emotion and passive exposure versus active regulation. Brain Struct. Funct. 225, 345–363. doi: 10.1007/s00429-019-02007-3
Keywords: simultaneous fMRI/PET, GABAA receptor, emotion, [11C]flumazenil, BOLD signal, binding potential
Citation: Heinzel A, Mauler J, Herzog H, Boers F, Mottaghy FM, Langen K-J, Scheins J, Lerche C, Neumaier B, Northoff G and Shah NJ (2023) GABAA receptor availability relates to emotion-induced BOLD responses in the medial prefrontal cortex: simultaneous fMRI/PET with [11C]flumazenil. Front. Neurosci. 17:1027697. doi: 10.3389/fnins.2023.1027697
Edited by:
Rajendra D. Badgaiyan, The University of Texas Health Science Center at San Antonio, United StatesReviewed by:
Vassiliy Tsytsarev, University of Maryland, United StatesChristin Y. Sander, Massachusetts General Hospital and Harvard Medical School, United States
Copyright © 2023 Heinzel, Mauler, Herzog, Boers, Mottaghy, Langen, Scheins, Lerche, Neumaier, Northoff and Shah. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alexander Heinzel, a.heinzel@fz-juelich.de
‡ORCID: Alexander Heinzel, https://orcid.org/0000-0002-2430-4557
Jörg Mauler, https://orcid.org/0000-0003-4937-0355
Hans Herzog, https://orcid.org/0000-0002-6751-7165
Frank Boers, https://orcid.org/0000-0002-5379-9368
Felix M. Mottaghy, https://orcid.org/0000-0002-7212-6521
Karl-Josef Langen, https://orcid.org/0000-0003-1101-5075
Jürgen Scheins, https://orcid.org/0000-0003-3474-430X
Christoph Lerche, https://orcid.org/0000-0003-2749-2108
Bernd Neumaier, https://orcid.org/0000-0001-5425-3116
Georg Northoff, https://orcid.org/0000-0002-5236-0951
N. Jon Shah, https://orcid.org/0000-0002-8151-6169
†These authors have contributed equally to this work