Validation of cardiac image-derived input functions for functional PET quantification.

PURPOSE: Functional PET (fPET) is a novel technique for studying dynamic changes in brain metabolism and neurotransmitter signaling. Accurate quantification of fPET relies on measuring the arterial input function (AIF), traditionally achieved through invasive arterial blood sampling. While non-invasive image-derived input functions (IDIF) offer an alternative, they suffer from limited spatial resolution and field of view. To overcome these issues, we developed and validated a scan protocol for brain fPET utilizing cardiac IDIF, aiming to mitigate known IDIF limitations. METHODS: Twenty healthy individuals underwent fPET/MR scans using [18F]FDG or 6-[18F]FDOPA, utilizing bed motion shuttling to capture cardiac IDIF and brain task-induced changes. Arterial and venous blood sampling was used to validate IDIFs. Participants performed a monetary incentive delay task. IDIFs from various blood pools and composites estimated from a linear fit over all IDIF blood pools (3VOI) and further supplemented with venous blood samples (3VOIVB) were compared to the AIF. Quantitative task-specific images from both tracers were compared to assess the performance of each input function to the gold standard. RESULTS: For both radiotracer cohorts, moderate to high agreement (r: 0.60-0.89) between IDIFs and AIF for both radiotracer cohorts was observed, with further improvement (r: 0.87-0.93) for composite IDIFs (3VOI and 3VOIVB). Both methods showed equivalent quantitative values and high agreement (r: 0.975-0.998) with AIF-derived measurements. CONCLUSION: Our proposed protocol enables accurate non-invasive estimation of the input function with full quantification of task-specific changes, addressing the limitations of IDIF for brain imaging by sampling larger blood pools over the thorax. These advancements increase applicability to any PET scanner and clinical research setting by reducing experimental complexity and increasing patient comfort.

Non-invasive assessment of stimulation-specific changes in cerebral glucose metabolism with functional PET.

PURPOSE: Functional positron emission tomography (fPET) with [18F]FDG allows quantification of stimulation-induced changes in glucose metabolism independent of neurovascular coupling. However, the gold standard for quantification requires invasive arterial blood sampling, limiting its widespread use. Here, we introduce a novel fPET method without the need for an input function. METHODS: We validated the approach using two datasets (DS). For DS1, 52 volunteers (23.2 ± 3.3 years, 24 females) performed Tetris® during a [18F]FDG fPET scan (bolus + constant infusion). For DS2, 18 participants (24.2 ± 4.3 years, 8 females) performed an eyes-open/finger tapping task (constant infusion). Task-specific changes in metabolism were assessed with the general linear model (GLM) and cerebral metabolic rate of glucose (CMRGlu) was quantified with the Patlak plot as reference. We then estimated simplified outcome parameters, including GLM beta values and percent signal change (%SC), and compared them, region and whole-brain-wise. RESULTS: We observed higher agreement with the reference for DS1 than DS2. Both DS resulted in strong correlations between regional task-specific beta estimates and CMRGlu (r = 0.763…0.912). %SC of beta values exhibited strong agreement with %SC of CMRGlu (r = 0.909…0.999). Average activation maps showed a high spatial similarity between CMRGlu and beta estimates (Dice = 0.870…0.979) as well as %SC (Dice = 0.932…0.997), respectively. CONCLUSION: The non-invasive method reliably estimates task-specific changes in glucose metabolism without blood sampling. This streamlines fPET, albeit with the trade-off of being unable to quantify baseline metabolism. The simplification enhances its applicability in research and clinical settings.

Case report: Interstitial pneumonitis after initiation of lamotrigine.

The second-generation anticonvulsant lamotrigine is widely used in the psychiatric field as a mood stabilizer or antidepressant augmentation therapy. Although particularly older anticonvulsants are known for their potential to cause hypersensitivity syndromes, newer antiepileptic drugs do hold a certain risk as well. Presenting a case of a 32-year-old male inpatient of African ethnicity suffering from a primary severe depressive episode in the course of a recurrent major depressive disorder, we report the occurrence of a rapid-onset drug-induced pneumonitis. Herewith, the interstitial pneumonitis occurred after the initiation of 25 mg lamotrigine as an augmentation therapy. Except for the clear temporal correlation between the administration of lamotrigine and the onset of pneumonitis, we did not reveal any further potentially causal diagnostic hints. Importantly, no relevant genetic variations of metabolizing enzymes or drug interactions resulting in lamotrigine overdosage as a potential cause of toxicity were identified. Our experience with a potentially life-threatening adverse drug reaction shortly after the initiation of the largely well-tolerated lamotrigine suggests a potential side effect under the second-generation anticonvulsant although similar adverse events are deemed to be very rare.

Comparison of cardiac image-derived input functions for quantitative whole body [18F]FDG imaging with arterial blood sampling.

Introduction: Dynamic positron emission tomography (PET) and the application of kinetic models can provide important quantitative information based on its temporal information. This however requires arterial blood sampling, which can be challenging to acquire. Nowadays, state-of-the-art PET/CT systems offer fully automated, whole-body (WB) kinetic modelling protocols using image-derived input functions (IDIF) to replace arterial blood sampling. Here, we compared the validity of an automatic WB kinetic model protocol to the reference standard arterial input function (AIF) for both clinical and research settings. Methods: Sixteen healthy participants underwent dynamic WB [18F]FDG scans using a continuous bed motion PET/CT system with simultaneous arterial blood sampling. Multiple processing pipelines that included automatic and manually generated IDIFs derived from the aorta and left ventricle, with and without motion correction were compared to the AIF. Subsequently generated quantitative images of glucose metabolism were compared to evaluate performance of the different input functions. Results: We observed moderate to high correlations between IDIFs and the AIF regarding area under the curve (r = 0.49-0.89) as well as for the cerebral metabolic rate of glucose (CMRGlu) (r = 0.68-0.95). Manual placing of IDIFs and motion correction further improved their similarity to the AIF. Discussion: In general, the automatic vendor protocol is a feasible approach for the quantification of CMRGlu for both, clinical and research settings where expertise or time is not available. However, we advise on a rigorous inspection of the placement of the volume of interest, the resulting IDIF, and the quantitative values to ensure valid interpretations. In protocols requiring longer scan times or where cohorts are prone to involuntary movement, manual IDIF definition with additional motion correction is recommended, as this has greater accuracy and reliability.

Whole-body metabolic connectivity framework with functional PET.

The nervous and circulatory system interconnects the various organs of the human body, building hierarchically organized subsystems, enabling fine-tuned, metabolically expensive brain-body and inter-organ crosstalk to appropriately adapt to internal and external demands. A deviation or failure in the function of a single organ or subsystem could trigger unforeseen biases or dysfunctions of the entire network, leading to maladaptive physiological or psychological responses. Therefore, quantifying these networks in healthy individuals and patients may help further our understanding of complex disorders involving body-brain crosstalk. Here we present a generalized framework to automatically estimate metabolic inter-organ connectivity utilizing whole-body functional positron emission tomography (fPET). The developed framework was applied to 16 healthy subjects (mean age ± SD, 25 ± 6 years; 13 female) that underwent one dynamic 18F-FDG PET/CT scan. Multiple procedures of organ segmentation (manual, automatic, circular volumes) and connectivity estimation (polynomial fitting, spatiotemporal filtering, covariance matrices) were compared to provide an optimized thorough overview of the workflow. The proposed approach was able to estimate the metabolic connectivity patterns within brain regions and organs as well as their interactions. Automated organ delineation, but not simplified circular volumes, showed high agreement with manual delineation. Polynomial fitting yielded similar connectivity as spatiotemporal filtering at the individual subject level. Furthermore, connectivity measures and group-level covariance matrices did not match. The strongest brain-body connectivity was observed for the liver and kidneys. The proposed framework offers novel opportunities towards analyzing metabolic function from a systemic, hierarchical perspective in a multitude of physiological pathological states.

Effects of bilateral sequential theta-burst stimulation on functional connectivity in treatment-resistant depression: First results.

BACKGROUND: Previous studies suggest that transcranial magnetic stimulation exerts antidepressant effects by altering functional connectivity (FC). However, knowledge about this mechanism is still limited. Here, we aimed to investigate the effect of bilateral sequential theta-burst stimulation (TBS) on FC in treatment-resistant depression (TRD) in a sham-controlled longitudinal study. METHODS: TRD patients (n = 20) underwent a three-week treatment of intermittent TBS of the left and continuous TBS of the right dorsolateral prefrontal cortex (DLPFC). Upon this trial's premature termination, 15 patients had received active TBS and five patients sham stimulation. Resting-state functional magnetic resonance imaging was performed at baseline and after treatment. FC (left and right DLPFC) was estimated for each participant, followed by group statistics (t-tests). Furthermore, depression scores were analyzed (linear mixed models analysis) and tested for correlation with FC. RESULTS: Both groups exhibited reductions of depression scores, however, there was no significant main effect of group, or group and time. Anticorrelations between DLPFC and the subgenual cingulate cortex (sgACC) were observed for baseline FC, corresponding to changes in depression severity. Treatment did not significantly change DLPFC-sgACC connectivity, but significantly reduced FC between the left stimulation target and bilateral anterior insula. CONCLUSIONS: Our data is compatible with previous reports on the relevance of anticorrelation between DLPFC and sgACC for treatment success. Furthermore, FC changes between left DLPFC and bilateral anterior insula highlight the effect of TBS on the salience network. LIMITATIONS: Due to the limited sample size, results should be interpreted with caution and are of exploratory nature.

ABCB1 variants and sex affect serotonin transporter occupancy in the brain.

Strategies to personalize psychopharmacological treatment promise to improve efficacy and tolerability. We measured serotonin transporter occupancy immediately after infusion of the widely prescribed P-glycoprotein substrate citalopram and assessed to what extent variants of the ABCB1 gene affect drug target engagement in the brain in vivo. A total of 79 participants (39 female) including 31 patients with major depression and 48 healthy volunteers underwent two PET/MRI scans with the tracer [11C]DASB and placebo-controlled infusion of citalopram (8 mg) in a cross-over design. We tested the effect of six ABCB1 single nucleotide polymorphisms and found lower SERT occupancy in ABCB1 rs2235015 minor allele carriers (n = 26, MAF = 0.18) compared to major allele homozygotes (t73 = 2.73, pFWE < 0.05) as well as in men compared to women (t73 = 3.33, pFWE < 0.05). These effects were robust to correction for citalopram plasma concentration, age and diagnosis. From occupancy we derived the ratio of occupied to unoccupied SERT, because in theory this measure is equal to the product of drug affinity and concentration at target sites. A model combining genotype with basic clinical variables, predicted that, at the same dosage, occupied to unoccupied SERT ratio was -14.48 ± 5.38% lower in rs2235015 minor allele carriers, +19.10 ± 6.95% higher in women, -4.83 ± 2.70% lower per 10 kg bodyweight, and -2.68 ± 3.07% lower per 10 years of age. Our results support the exploration of clinical algorithms with adjustment of initial citalopram dosing and highlight the potential of imaging-genetics for precision pharmacotherapy in psychiatry.

Biodistribution and dosimetry of the GluN2B-specific NMDA receptor PET radioligand (R)-[11C]Me-NB1.

BACKGROUND: The NMDA receptor (NMDAR) plays a key role in the central nervous system, e.g., for synaptic transmission. While synaptic NMDARs are thought to have protective characteristics, activation of extrasynaptic NMDARs might trigger excitotoxic processes linked to neuropsychiatric disorders. Since extrasynaptic NMDARs are typically GluN2B-enriched, the subunit is an interesting target for drug development and treatment monitoring. Recently, the novel GluN2B-specific PET radioligand (R)-[11C]Me-NB1 was investigated in rodents and for the first time successfully translated to humans. To assess whether (R)-[11C]Me-NB1 is a valuable radioligand for (repeated) clinical applications, we evaluated its safety, biodistribution and dosimetry. METHODS: Four healthy subjects (two females, two males) underwent one whole-body PET/MR measurement lasting for more than 120 min. The GluN2B-specific radioligand (R)-[11C]Me-NB1 was administered simultaneously with the PET start. Subjects were measured in nine passes and six bed positions from head to mid-thigh. Regions of interest was anatomically defined for the brain, thyroid, lungs, heart wall, spleen, stomach contents, pancreas, liver, kidneys, bone marrow and urinary bladder contents, using both PET and MR images. Time-integrated activity coefficients were estimated to calculate organ equivalent dose coefficients and the effective dose coefficient. Additionally, standardized uptake values (SUV) were computed to visualize the biodistribution. RESULTS: Administration of the radioligand was safe without adverse events. The organs with the highest uptake were the urinary bladder, spleen and pancreas. Organ equivalent dose coefficients were higher in female in almost all organs, except for the urinary bladder of male. The effective dose coefficient was 6.0 µSv/MBq. CONCLUSION: The GluN2B-specific radioligand (R)-[11C]Me-NB1 was well-tolerated without reported side effects. Effective dose was estimated to 1.8 mSv when using 300 MBq of presented radioligand. The critical organ was the urinary bladder. Due to the low effective dose coefficient of this radioligand, longitudinal studies for drug development and treatment monitoring of neuropsychiatric disorders including neurodegenerative diseases are possible. Trial registration Registered on 11th of June 2019 at https://www.basg.gv.at (EudraCT: 2018-002933-39).

The Impact of Theta-Burst Stimulation on Cortical GABA and Glutamate in Treatment-Resistant Depression: A Surface-Based MRSI Analysis Approach.

Background: Theta burst stimulation (TBS) belongs to one of the biological antidepressant treatment options. When applied bilaterally, excitatory intermittent TBS (iTBS) is commonly targeted to the left and inhibitory continuous TBS (cTBS) to the right dorsolateral prefrontal cortex. TBS was shown to influence neurotransmitter systems, while iTBS is thought to interfere with glutamatergic circuits and cTBS to mediate GABAergic neurotransmission. Objectives: We aimed to expand insights into the therapeutic effects of TBS on the GABAergic and glutamatergic system utilizing 3D-multivoxel magnetic resonance spectroscopy imaging (MRSI) in combination with a novel surface-based MRSI analysis approach to investigate changes of cortical neurotransmitter levels in patients with treatment-resistant depression (TRD). Methods: Twelve TRD patients (five females, mean age ± SD = 35 ± 11 years) completed paired MRSI measurements, using a GABA-edited 3D-multivoxel MEGA-LASER sequence, before and after 3 weeks of bilateral TBS treatment. Changes in cortical distributions of GABA+/tNAA (GABA+macromolecules relative to total N-acetylaspartate) and Glx/tNAA (Glx = mixed signal of glutamate and glutamine), were investigated in a surface-based region-of-interest (ROI) analysis approach. Results: ANCOVAs revealed a significant increase in Glx/tNAA ratios in the left caudal middle frontal area (p corr. = 0.046, F = 13.292), an area targeted by iTBS treatment. Whereas, contralateral treatment with cTBS evoked no alterations in glutamate or GABA concentrations. Conclusion: This study demonstrates surface-based adaptions in the stimulation area to the glutamate metabolism after excitatory iTBS but not after cTBS, using a novel surface-based analysis of 3D-MRSI data. The reported impact of facilitatory iTBS on glutamatergic neurotransmission provides further insight into the neurobiological effects of TBS in TRD.

Serotonergic modulation of effective connectivity in an associative relearning network during task and rest.

An essential core function of one's cognitive flexibility is the use of acquired knowledge and skills to adapt to ongoing environmental changes. Animal models have highlighted the influence serotonin has on neuroplasticity. These effects have been predominantly demonstrated during emotional relearning which is theorized as a possible model for depression. However, translation of these mechanisms is in its infancy. To this end, we assessed changes in effective connectivity at rest and during associative learning as a proxy of neuroplastic changes in healthy volunteers. 76 participants underwent 6 weeks of emotional or non-emotional (re)learning (face-matching or Chinese character-German noun matching). During relearning participants either self-administered 10 mg/day of the selective serotonin reuptake inhibitor (SSRI) escitalopram or placebo in a double-blind design. Associative learning tasks, resting-state and structural images were recorded before and after both learning phases (day 1, 21 and 42). Escitalopram intake modulated relearning changes in a network encompassing the right insula, anterior cingulate cortex and right angular gyrus. Here, the process of relearning during SSRI intake showed a greater decrease in effective connectivity from the right insula to both the anterior cingulate cortex and right angular gyrus, with increases in the opposite direction when compared to placebo. In contrast, intrinsic connections and those at resting-state were only marginally affected by escitalopram. Further investigation of gray matter volume changes in these functionally active regions revealed no significant SSRI-induced structural changes. These findings indicate that the right insula plays a central role in the process of relearning and SSRIs further potentiate this effect. In sum, we demonstrated that SSRIs amplify learning-induced effective connections rather than affecting the intrinsic task connectivity or that of resting-state.

Escitalopram modulates learning content-specific neuroplasticity of functional brain networks.

Learning-induced neuroplastic changes, further modulated by content and setting, are mirrored in brain functional connectivity (FC). In animal models, selective serotonin reuptake inhibitors (SSRIs) have been shown to facilitate neuroplasticity. This is especially prominent during emotional relearning, such as fear extinction, which may translate to clinical improvements in patients. To investigate a comparable modulation of neuroplasticity in humans, 99 healthy subjects underwent three weeks of emotional (matching faces) or non-emotional learning (matching Chinese characters to unrelated German nouns). Shuffled pairings of the original content were subsequently relearned for the same time. During relearning, subjects received either a daily dose of the SSRI escitalopram or placebo. Resting-state functional magnetic resonance imaging was performed before and after the (re-)learning phases. FC changes in a network comprising Broca's area, the medial prefrontal cortex, the right inferior temporal and left lingual gyrus were modulated by escitalopram intake. More specifically, it increased the bidirectional connectivity between medial prefrontal cortex and lingual gyrus for non-emotional and the connectivity from medial prefrontal cortex to Broca's area for emotional relearning. The context dependence of these effects together with behavioral correlations supports the assumption that SSRIs in clinical practice improve neuroplasticity rather than psychiatric symptoms per se. Beyond expanding the complexities of learning, these findings emphasize the influence of external factors on human neuroplasticity.

Predicting Antidepressant Citalopram Treatment Response via Changes in Brain Functional Connectivity After Acute Intravenous Challenge.

Introduction: The early and therapy-specific prediction of treatment success in major depressive disorder is of paramount importance due to high lifetime prevalence, and heterogeneity of response to standard medication and symptom expression. Hence, this study assessed the predictability of long-term antidepressant effects of escitalopram based on the short-term influence of citalopram on functional connectivity. Methods: Twenty nine subjects suffering from major depression were scanned twice with resting-state functional magnetic resonance imaging under the influence of intravenous citalopram and placebo in a randomized, double-blinded cross-over fashion. Symptom factors were identified for the Hamilton depression rating scale (HAM-D) and Beck's depression inventory (BDI) taken before and after a median of seven weeks of escitalopram therapy. Predictors were calculated from whole-brain functional connectivity, fed into robust regression models, and cross-validated. Results: Significant predictive power could be demonstrated for one HAM-D factor describing insomnia and the total score (r = 0.45-0.55). Remission and response could furthermore be predicted with an area under the receiver operating characteristic curve of 0.73 and 0.68, respectively. Functional regions with high influence on the predictor were located especially in the ventral attention, fronto-parietal, and default mode networks. Conclusion: It was shown that medication-specific antidepressant symptom improvements can be predicted using functional connectivity measured during acute pharmacological challenge as an easily assessable imaging marker. The regions with high influence have previously been related to major depression as well as the response to selective serotonin reuptake inhibitors, corroborating the advantages of the current approach of focusing on treatment-specific symptom improvements.

Reinforcement and Punishment Shape the Learning Dynamics in fMRI Neurofeedback.

INTRODUCTION: Neurofeedback (NF) using real-time functional magnetic resonance imaging (fMRI) has proven to be a valuable neuroscientific tool for probing cognition and promising therapeutic approach for several psychiatric disorders. Even though learning constitutes an elementary aspect of NF, the question whether certain training schemes might positively influence its dynamics has largely been neglected. METHODS: To address this issue, participants were trained to exert control on their subgenual anterior cingulate cortex (sgACC) blood-oxygenation-level-dependent signal, receiving either exclusively positive reinforcement (PR, "positive feedback") or also positive punishment (PP, "negative feedback"). The temporal dynamics of the learning process were investigated by individually modeling the feedback periods and trends, offering the possibility to assess activation changes within and across blocks, runs and sessions. RESULTS: The results show faster initial learning of the PR + PP group by significantly lower deactivations of the sgACC in the first session and stronger regulation trends during the first runs. Independent of the group, significant control over the sgACC could further be shown with but not without feedback. CONCLUSION: The beneficial effect of PP is supported by previous findings of multiple research domains suggesting that error avoidance represents an important motivational factor of learning, which complements the reward spectrum. This hypothesis warrants further investigation with respect to NF, as it could offer a way to generally facilitate the process of gaining volitional control over brain activity.

Modeling the acute pharmacological response to selective serotonin reuptake inhibitors in human brain using simultaneous PET/MR imaging.

Pharmacological imaging of the effects of selective serotonin reuptake inhibitors (SSRI) may aid the clarification of their mechanism of action and influence treatment of highly prevalent neuropsychiatric conditions if the detected effects could be related to patient outcomes. In a randomized double-blind design, 38 healthy participants received a constant infusion of 8 mg citalopram or saline during either their first or second of two PET/MR scans. Resting-state functional MRI (fMRI) was acquired simultaneously with PET data on the binding of serotonin transporters (5-HTT) using [11C]DASB. Three different approaches for modeling of pharmacological fMRI response were tested separately. These relied on the use of regressors corresponding to (1) the drug infusion paradigm, (2) time courses of citalopram plasma concentrations and (3) changes in 5-HTT binding measured in each individual, respectively. Furthermore, the replication of results of a widely used model-free analysis method was attempted which assesses the deviation of signal in discrete time bins of fMRI data acquired after start of drug infusion. Following drug challenge, average 5-HTT occupancy was 69±7% and peak citalopram plasma levels were 111.8 ±â€¯21.1 ng/ml. None of the applied methods could detect significant differences in the pharmacological response between SSRI and placebo scans. The failed replication of SSRI effects reported in the literature despite a threefold larger sample size highlights the importance of appropriate correction for family-wise error in order to avoid spurious results in pharmacological imaging. This calls for the development of analysis methods which take regional specialization and the dynamics of brain activity into account.

Structural changes in amygdala nuclei, hippocampal subfields and cortical thickness following electroconvulsive therapy in treatment-resistant depression: longitudinal analysis.

BACKGROUND: Electroconvulsive therapy (ECT) is the treatment of choice for severe mental illness including treatment-resistant depression (TRD). Increases in volume of the hippocampus and amygdala following ECT have consistently been reported.AimsTo investigate neuroplastic changes after ECT in specific hippocampal subfields and amygdala nuclei using high-resolution structural magnetic resonance imaging (MRI) (trial registration: clinicaltrials.gov - NCT02379767). METHOD: MRI scans were carried out in 14 patients (11 women, 46.9 years (s.d. = 8.1)) with unipolar TRD twice before and once after a series of right unilateral ECT in a pre-post study design. Volumes of subcortical structures, including subfields of the hippocampus and amygdala, and cortical thickness were extracted using FreeSurfer. The effect of ECT was tested using repeated-measures ANOVA. Correlations of imaging and clinical parameters were explored. RESULTS: Increases in volume of the right hippocampus by 139.4 mm3 (s.d. = 34.9), right amygdala by 82.3 mm3 (s.d. = 43.9) and right putamen by 73.9 mm3 (s.d. = 77.0) were observed. These changes were localised in the basal and lateral nuclei, and the corticoamygdaloid transition area of the amygdala, the hippocampal-amygdaloid transition area and the granule cell and molecular layer of the dentate gyrus. Cortical thickness increased in the temporal, parietal and insular cortices of the right hemisphere. CONCLUSIONS: Following ECT structural changes were observed in hippocampal subfields and amygdala nuclei that are specifically implicated in the pathophysiology of depression and stress-related disorders and retain a high potential for neuroplasticity in adulthood.Declaration of interestS.K. has received grants/research support, consulting fees and/or honoraria within the past 3 years from Angelini, AOP Orphan Pharmaceuticals AG, AstraZeneca, Celegne GmbH, Eli Lilly, Janssen-Cilag Pharma GmbH, KRKA-Pharma, Lundbeck A/S, Neuraxpharm, Pfizer, Pierre Fabre, Schwabe and Servier. R.L. received travel grants and/or conference speaker honoraria from Shire, AstraZeneca, Lundbeck A/S, Dr. Willmar Schwabe GmbH, Orphan Pharmaceuticals AG, Janssen-Cilag Pharma GmbH, and Roche Austria GmbH.

Spatial analysis and high resolution mapping of the human whole-brain transcriptome for integrative analysis in neuroimaging.

The quantification of big pools of diverse molecules provides important insights on brain function, but is often restricted to a limited number of observations, which impairs integration with other modalities. To resolve this issue, a method allowing for the prediction of mRNA expression in the entire brain based on microarray data provided in the Allen Human Brain Atlas was developed. Microarray data of 3702 samples from 6 brain donors was registered to MNI and cortical surface space using FreeSurfer. For each of 18,686 genes, spatial dependence of transcription was assessed using variogram modelling. Variogram models were employed in Gaussian process regression to calculate best linear unbiased predictions for gene expression at all locations represented in well-established imaging atlases for cortex, subcortical structures and cerebellum. For validation, predicted whole-brain transcription of the HTR1A gene was correlated with [carbonyl-11C]WAY-100635 positron emission tomography data collected from 30 healthy subjects. Prediction results showed minimal bias ranging within ±0.016 (cortical surface), ±0.12 (subcortical regions) and ±0.14 (cerebellum) in units of log2 expression intensity for all genes. Across genes, the correlation of predicted and observed mRNA expression in leave-one-out cross-validation correlated with the strength of spatial dependence (cortical surface: r = 0.91, subcortical regions: r = 0.85, cerebellum: r = 0.84). 816 out of 18,686 genes exhibited a high spatial dependence accounting for more than 50% of variance in the difference of gene expression on the cortical surface. In subcortical regions and cerebellum, different sets of genes were implicated by high spatially structured variability. For the serotonin 1A receptor, correlation between PET binding potentials and predicted comprehensive mRNA expression was markedly higher (Spearman ρ = 0.72 for cortical surface, ρ = 0.84 for subcortical regions) than correlation of PET and discrete samples only (ρ = 0.55 and ρ = 0.63, respectively). Prediction of mRNA expression in the entire human brain allows for intuitive visualization of gene transcription and seamless integration in multimodal analysis without bias arising from non-uniform distribution of available samples. Extension of this methodology promises to facilitate translation of omics research and enable investigation of human brain function at a systems level.