Arkadiusz Komorowski1, Ana Weidenauer1, Matej Murgaš1, Ulrich Sauerzopf1, Wolfgang Wadsak2, Markus Mitterhauser3, Martin Bauer4, Marcus Hacker5, Nicole Praschak-Rieder1, Siegfried Kasper6, Rupert Lanzenberger7, Matthäus Willeit1. 1. Department of Psychiatry and Psychotherapy, Division of General Psychiatry, Medical University of Vienna, Währinger Gürtel 18-20, 1090 Vienna, Austria. 2. Department of Biomedical Imaging and Image-guided Therapy, Division of Nuclear Medicine, Medical University of Vienna, Vienna, Austria; Center for Biomarker Research in Medicine (CBmed), Graz, Austria. 3. Department of Biomedical Imaging and Image-guided Therapy, Division of Nuclear Medicine, Medical University of Vienna, Vienna, Austria; Ludwig Boltzmann Institute for Applied Diagnostics, Vienna, Austria. 4. Department of Clinical Pharmacology, Medical University of Vienna, Vienna, Austria. 5. Department of Biomedical Imaging and Image-guided Therapy, Division of Nuclear Medicine, Medical University of Vienna, Vienna, Austria. 6. Center for Brain Research, Medical University of Vienna, Vienna, Austria. 7. Department of Psychiatry and Psychotherapy, Division of General Psychiatry, Medical University of Vienna, Währinger Gürtel 18-20, 1090 Vienna, Austria. Electronic address: rupert.lanzenberger@meduniwien.ac.at.
Abstract
Open access post-mortem transcriptome atlases such as the Allen Human Brain Atlas (AHBA) can inform us about mRNA expression of numerous proteins of interest across the whole brain, while in vivo protein binding in the human brain can be quantified by means of neuroreceptor positron emission tomography (PET). By combining both modalities, the association between regional gene expression and receptor distribution in the living brain can be approximated. Here, we compare the characteristics of D2 and D3 dopamine receptor distribution by applying the dopamine D2/3 receptor agonist radioligand [11C]-(+)-PHNO and human gene expression data. Since [11C]-(+)-PHNO has a higher affinity for D3 compared to D2 receptors, we hypothesized that there is a stronger relationship between D2/3 non-displaceable binding potentials (BPND) and D3 mRNA expression. To investigate the relationship between D2/3 BPND and mRNA expression of DRD2 and DRD3 we performed [11C]-(+)-PHNO PET scans in 27 healthy subjects (12 females) and extracted gene expression data from the AHBA. We also calculated D2/D3 mRNA expression ratios to imitate the mixed D2/3 signal of [11C]-(+)-PHNO. In accordance with our a priori hypothesis, a strong correlation between [11C]-(+)-PHNO and DRD3 expression was found. However, there was no significant correlation with DRD2 expression. Calculated D2/D3 mRNA expression ratios also showed a positive correlation with [11C]-(+)-PHNO binding, reflecting the mixed D2/3 signal of the radioligand. Our study supports the usefulness of combining gene expression data from open access brain atlases with in vivo imaging data in order to gain more detailed knowledge on neurotransmitter signaling.
Open access post-mortem transcriptome atlases such as the Allen Human Brain Atlas (AHBA) can inform us about mRNA expression of numerous proteins of interest across the whole brain, while in vivo protein binding in the human brain can be quantified by means of neuroreceptor positron emission tomography (PET). By combining both modalities, the association between regional gene expression and receptor distribution in the living brain can be approximated. Here, we compare the characteristics of D2 and D3 dopamine receptor distribution by applying the dopamine D2/3 receptor agonist radioligand [11C]-(+)-PHNO and human gene expression data. Since [11C]-(+)-PHNO has a higher affinity for D3 compared to D2 receptors, we hypothesized that there is a stronger relationship between D2/3 non-displaceable binding potentials (BPND) and D3 mRNA expression. To investigate the relationship between D2/3 BPND and mRNA expression of DRD2 and DRD3 we performed [11C]-(+)-PHNO PET scans in 27 healthy subjects (12 females) and extracted gene expression data from the AHBA. We also calculated D2/D3 mRNA expression ratios to imitate the mixed D2/3 signal of [11C]-(+)-PHNO. In accordance with our a priori hypothesis, a strong correlation between [11C]-(+)-PHNO and DRD3 expression was found. However, there was no significant correlation with DRD2 expression. Calculated D2/D3 mRNA expression ratios also showed a positive correlation with [11C]-(+)-PHNO binding, reflecting the mixed D2/3 signal of the radioligand. Our study supports the usefulness of combining gene expression data from open access brain atlases with in vivo imaging data in order to gain more detailed knowledge on neurotransmitter signaling.
Increasingly sophisticated neuroimaging methods and data-sharing initiatives have furthered the understanding of human brain function. Using the Allen Human Brain Atlas (AHBA; http://human.brainmap.org) (Hawrylycz et al., 2012), numerous studies investigating the relationship between structural and functional brain parameters and the transcriptome were published recently (Burt, 2018; Goel et al., 2014; Richiardi et al., 2015; Shin et al., 2017). Several authors showed close associations of gene expression levels with in vivo protein densities, mainly focusing on serotonergic neurotransmitter receptors (Komorowski et al., 2017; Rizzo et al., 2016; Unterholzner et al., 2020). Similar to the serotonergic system, dopamine signaling is mediated through G-protein-coupled receptors and plays a crucial role in physiologic and pathologic processes in humans, including neuropsychiatric conditions (Beaulieu and Gainetdinov, 2011). To identify pathophysiological mechanisms within the dopamine system, different neuroimaging modalities including high-resolution whole-brain transcriptome maps (Gryglewski et al., 2018) and positron emission tomography (PET) imaging data of dopamine receptors can be integrated, extending previous studies on dopamine receptors. Labelled with carbon-11, (+)-4-propyl-3,4,4a,5,6,10b-hexahydro-2H-naphtho[1,2-b][1,4]oxazin-9-ol [11C]-(+)-PHNO) (Wilson et al., 2005), a radioligand with full agonistic properties at dopamine D2/3 receptors (Brown et al., 1997) offers an excellent signal-to-noise (SNR) ratio and favourable kinetics for PET imaging in humans (Ginovart et al., 2007; Jensen et al., 2007). Many studies using [11C]-(+)-PHNO have focused on schizophrenia (Graff-Guerrero et al., 2009; Mizrahi et al., 2012; Mizrahi et al., 2011; Mizrahi et al., 2014; Suridjan et al., 2013; Weidenauer et al., 2020), or substance use disorders (Boileau et al., 2015; Boileau et al., 2012; Payer et al., 2014; Worhunsky et al., 2017). Applying PET, several approaches have been described to determine the relative affinity of [11C]-(+)-PHNO as well as the regionally distinctive topology of D2 compared to D3 receptors, including pharmacological and animal studies (Graff-Guerrero et al., 2009). Depending on G-protein binding, the dopamine D2 receptor occurs in an interconvertible high- and low-affinity state (D2
high vs. D2
low) towards its natural agonist dopamine, with the high-affinity state being the functionally active one (George et al., 1985). In contrast, G-protein binding has comparably little influence on the affinity of D3 receptors for dopamine (Sokoloff et al., 1992). Thus, D3 receptors are assumed to be “locked” in the high-affinity state, affecting binding of D2/3 radiotracers with a yet unknown fraction. It is well supported that [11C]-(+)-PHNO binds more to D2 high than to D 2 low (Nobrega and Seeman, 1994; Seeman et al., 1993), and PET imaging studies in non-human primates indicate that in vivo [11C]-(+)-PHNO has a 25-48 fold higher affinity for D3 over D2 receptors (Gallezot et al., 2012). Being easily displaced by endogenous dopamine, [11C]-(+)-PHNO is a useful tool for challenge studies (Ginovart et al., 2006; Weidenauer et al., 2020; Willeit et al., 2008; Willeit et al., 2006). Still, although several studies addressed this issue in animals and humans (Graff-Guerrero et al., 2010; Rabiner et al., 2009; Searle et al., 2010; Seeman et al., 2007; Tziortzi et al., 2011), the relative contribution of extracellular dopamine, D2
high, D2
low and D3 binding to the [11C]-(+)-PHNO overall signal is insufficiently understood to this day.In this study, we aimed to contribute to the disentangling of the [11C]-(+)-PHNO PET signal into the respective fractions of D2 vs. D3 receptors by relating subcortical non-displaceable binding potential (BPND) values of healthy human subjects to gene expression data obtained from the AHBA and interpolated transcriptome maps according to Gryglewski et al. (Gryglewski et al., 2018).
Methods
Participants
Imaging data were acquired in a larger study investigating dopamine release in healthy subjects and patients with schizophrenia (Weidenauer et al., 2020) (Clinical Trial Registry: EUDRACT 2010-019586-29). Here, analysis was restricted to healthy participants. All study procedures were approved by the Ethics Committee of the Medical University of Vienna and pertinent federal regulatory authorities. After being informed on study procedures and study-related risks, participants gave written informed consent. Subjects were only included if they were free of medical or neurological disorders and if laboratory or urine tests did not show clinically relevant abnormalities. All participants refrained from alcohol consumption as well as other medication before PET imaging. After careful screening, 34 participants were enrolled into this study from 2012-12-04 to 2017-05-31 at the Department of Psychiatry and Psychotherapy, Medical University of Vienna. Seven subjects did not enter analysis due to poor data quality, leaving 27 healthy participants (mean age: 26.22, SD: 2.36, ranging from 22 to 31; 12 females) for final data analysis.
Radiosynthesis of [11C]-(+)-PHNO has been described in detail elsewhere (Rami-Mark et al., 2013). PET imaging was performed using a GE Advance scanner (General Electric Medical Systems, Milwaukee, WI, USA). After injection of 315.7 MBq ± 76.2 [11C]-(+)-PHNO, emission data were obtained over 90 min. Reconstruction of raw data was performed by filtered-back projection to return dynamic images with 15 one min frames and 15 five min frames. The software AFNI (AFNI_17.3.06, NIMH, Bethesda, MD, USA) was used for frame-wise motion correction and co-registration of average PET images to T1-weighted magnetic resonance imaging (MRI) data. MRI images were normalized to the Montreal Neurological Institute (MNI) space and transformation matrices were used for normalization of PET images. Decay-corrected time activity curves were extracted from dynamic PET images using regions of interest (ROIs) from the high-resolution probabilistic in vivo atlas of subcortical nuclei by Pauli et al. (Pauli et al., 2018), which includes the majority of relevant subcortical structures. [11C]-(+)-PHNO BPnd values were calculated using the simplified reference tissue model 2 (SRTM2) (Ginovart et al., 2007; Wu and Carson, 2002) as implemented in PMOD software (Version 3.6; PMOD Technologies Ltd, Zurich, Switzerland). The cerebellum was used as a reference region, since it is virtually devoid of dopamine D2/3 receptors (Camps et al., 1989; Hall et al., 1996; Levant, 1998).To test the reliability of local imaging procedures, five male healthy subjects (mean age: 26, SD: 1.1) underwent two [11C]-(+)-PHNO PET scans at least one weak apart. Data were preprocessed and analyzed as described above. Subsequently, test-retest reliability was assessed using Pearson Interclass Correlation Coefficients (Liljequist et al., 2019).
Magnetic resonance imaging
MRI measurements were conducted on an Achieva 3.0 T (Philips Medical Systems, Best, The Netherlands) or a Magnetom Skyra 3.0 T scanner (Siemens Healthineers, Erlangen, Germany) with a spatial resolution of 0.9 × 0.9 × 1.1 mm using comparable acquisition protocols.
Whole-brain gene expression
The AHBA comprises up to 3702 microarray samples with a voxel size of 2 mm across the whole brain from six healthy subjects (mean age: 42.5, SD: 13.4; one female). A detailed summary of mRNA data acquisition and processing was published upon creation of the atlas (Hawrylycz et al., 2012; Shen et al., 2012). Within the AHBA, gene expression values (log2) of the D2 and D 3 receptor are represented in MNI space by 5 and 7 probes, respectively. Out of all available probes for both receptors, the median mRNA values were included for further analyses. DRD2 and DRD3 gene expression data were obtained only for subcortical regions in order to perform region-wise analyses in ROIs that show sufficiently reliable [11C]-(+)-PHNO BPND and to avoid bias due to disparate gene expression levels between cortical and subcortical regions (Chen et al., 2016; Kirsch and Chechik, 2016). Due to the unbalanced sex distribution of the AHBA, differences in gene expression were not analyzed separately for male and female brain donors.Original gene expression data were averaged within ROIs and subsequently combined into a brain template, including mRNA expression from all AHBA donor brains. In case of identical sample coordinates within the resulting brain template, overlapping mRNA expression values were averaged beforehand within these voxels. Despite widespread anatomical coverage, data prediction from the AHBA is potentially hampered in small subcortical regions due to sparse anatomical sampling. For compensating this sampling bias, we additionally obtained interpolated predictions of gene expression for D2 and D3 receptors from corresponding transcriptome maps recently published by Gryglewski et al. (Gryglewski et al., 2018). Except for the receptors investigated in this study, further whole-brain transcriptome maps based on Gaussian process regression are available online at www.meduniwien.ac.at/neuroimaging/mRNA.html for a total of 18,686 genes.
Statistical analysis
Statistical analysis was performed using R (Version 3.1.1, http://www.R-project.org/). Group-averaged t11C]-(+)-PHNO BPND values were used for correlation analyses with gene expression data. Original mRNA data from the AHBA as well as interpolated expression values from whole-brain transcriptome maps were extracted within corresponding ROIs to allow further correlation with imaging parameters. Variability of PET and mRNA data within each brain region was evaluated in terms of mean and standard deviation values. The brain atlas by Pauli et al. (Pauli et al., 2018) was used both for brain parcellation of imaging data as well as for region-wise extraction of gene expression values. To evaluate possible influences of sex, [11C]-(+)-PHNO BPND values were also extracted for female and male subjects separately. Additionally, the relationship between radiotracer binding and age was examined to rule out a potentially confounding effect, since there is conflicting literature on this topic (Nakajima et al., 2015, Matuskey et al., 2016).To unravel information about binding to D2 and D3 receptors from [11C]-(+)-PHNO data, radioligand BPND were compared with gene expression values by means of Pearson Product Moment correlation coefficients after testing for normality using Shapiro-Wilk tests. Additionally, associations were assessed by means of Spearman’s correlation coefficients. Binding of [11C]-(+)-PHNO to D2
low receptors was considered negligible, leaving D2
high and D3 fractions of [11C]-(+)-PHNO BPND for statistical analysis. Given the mixed D2/3 signal of [11C]-(+)-PHNO with a stronger affinity to D3 receptors, we alternatively estimated the D3 ratio (D3 % mRNA) based on gene expression data. D3% mRNA was calculated from the ratio between DRD3 and the sum of DRD2 and DRD3 expression values (D3 % mRNA: D3/(D2 + D3)* 100).
Results
[11C]-(+)-PHNO BPND values showed the expected binding pattern described earlier (Willeit et al., 2008). Briefly, high radiotracer enrichment was found in the ventral pallidum (BPND : 3.14 ± 0.45), nucleus accumbens (BPND : 2.48 ± 0.31), as well as globus pallidus externus (BPND : 2.45 ± 0.35) and internus (BPND : 2.12 ± 0.41). Likewise, mRNA levels varied across these brain regions with comparable expression values between original AHBA and interpolated transcriptome data (see Table 1 and Fig. 1). There was no significant relationship between [11C]-(+)-PHNO BPnd and age in any of the analyzed regions (p > 0.05).
Table 1
[11C]-(+)-PHNO binding potentials obtained from 27 healthy subjects and mRNA expression from the Allen Human Brain Atlas are provided within regions-of-interest labeled according to Pauli et al. (Pauli et al., 2018). Gene expression values were extracted from interpolated transcriptome maps provided by Gryglewski et al. (Gryglewski et al., 2018).
ROI (Pauli)
[11C]-(+)-PHNO BPND
D2 rec. mRNA (log2)
D3 rec. mRNA (log2)
CAU
1.00 ± 0.16
5.19 ± 0.50
3.85 ± 1.06
EA
1.61 ± 0.25
4.49 ± 0.20
3.02 ± 0.69
GP ext.
2.45 ± 0.35
4.45 ± 0.22
2.55 ± 0.59
GP int.
2.12 ± 0.41
4.23 ± 0.10
2.19 ± 0.33
HN
0.11 ± 0.05
4.30 ± 0.12
1.83 ± 0.12
Hypo
0.91 ± 0.16
4.49 ± 0.08
1.97 ± 0.28
NAc
2.48 ± 0.31
4.40 ± 0.67
3.77 ± 1.46
PPN
0.49 ± 0.10
4.86 ± 0.12
2.09 ± 0.06
PUT
1.71 ± 0.22
4.85 ± 0.50
3.13 ± 0.95
RN
0.38 ± 0.08
4.75 ± 0.09
2.17 ± 0.08
SNC
0.59 ± 0.11
4.99 ± 0.40
2.12 ± 0.06
SNR
0.62 ± 0.11
5.00 ± 0.05
2.12 ± 0.07
VP
3.14 ± 0.45
4.49 ± 0.12
3.95 ± 0.74
VT
0.56 ± 0.12
4.83 ± 0.13
2.07 ± 0.04
Legend: ROI: region-of-interest, Pauli: Pauli atlas, CAU: caudate nucleus, EA: extended amygdala, GP ext.: globus pallidus externus, GP int.: globus pallidus internus, HN: habenular nuclei, Hypo: hypothalamus, NAc: nucleus accumbens, PPN: parabrachial pigmented nucleus, PUT: putamen, RN: red nucleus, SNC: substantia nigra pars compacta, SNR: substantia nigra pars reticulata, VP: ventral pallidum, VT: ventral tegmentum.
Fig. 1
Statistical parametric maps displaying distributions of dopamine D2 and D3 receptors. Interpolated gene expression patterns (log2) according to Gryglewski et al. (Gryglewski et al., 2018) are depicted in MNI space (x = −12, y = 10, z = −10) for both dopamine receptors in the upper two panels, while [11C]-(+)-PHNO binding (BPND) is provided in the lowest panel.
Data reliability in the test-retest sample according to Pearson interclass correlation coefficients was excellent (> 0.9) for the caudate and the putamen, good (> 0.75) for the extended amygdala, globus pallidus, habenular nuclei, hypothalamus, mammillary nucleus, nucleus accumbens, and moderate (> 0.5) for the parabrachial pigmented nucleus, red nucleus, substantia nigra compacta et reticulata, subthalamic nucleus, ventral pallidum and ventral tegmentum. Remarkably, no regions showed poor reliability (< 0.5).Associations of gene expression data with [11C]-(+)-PHNO BPND were primarily evaluated by means of Pearson Product Moment correlation coefficients. Correlation analyses using interpolated mRNA data yielded similar results for all analyses compared to the original AHBA data. As hypothesized, DRD3 expression data showed a strong correlation with [11C]-(+)-PHNO BPND across the whole brain (interpolated: r = 0.69, p = 0.006; original data: r = 0.71, p < 0.01). Brain areas rich in D3 receptors showed both high radioligand binding and high mRNA expression (see Fig. 2). In contrast, DRD2 expression showed no significant correlation with [11C]-(+)-PHNO BPND values (interpolated: r = -0.46, p = 0.1; original data: r = 0.11, p = 0.68). A significant correlation between [11C]-(+)-PHNO BPND and D3 % mRNA was also observed (interpolated: r = 0.8, p < 0.01; original data: r = 0.51, p = 0.04), supporting the D3 binding preference of [11C]-(+)-PHNO (see Fig. 3). The linear associations were subsequently confirmed by Spearman’s rank correlation coefficients, which yielded similar results (interpolated: D3 rho = 0.78, p = 0.001; D2 rho = -0.31, p = 0.28; original data: D3 rho = 0.68, p = 0.004; D 2 rho = 0.14, p = 0.60). Likewise, separate analyses for female and male subjects yielded comparable results, both for the D3 (female interpolated: r = 0.7, p = 0.005; male interpolated: r = 0.71, p = 0.004) and D2 receptor (female interpolated: r = -0.47, p = 0.09; male interpolated: r = -0.43, p = 0.13).
Fig. 2
Relationship between [11C]-(+2-PHNO BPND and interpolated mRNA data (log2) for A) the D2 receptor (t = -1.79, df = 12, r = -0.46, p = 0.1) and B) the D3 receptor (t = 3.30, df = 12, r = 0.69, p = 0.006). Regions-of-interest were labeled according to Pauli et al. (Pauli et al., 2018).
Legend: ROI: region-of-interest, Pauli: Pauli atlas, CAU: caudate nucleus, EA: extended amygdala, GP ext.: globus pallidus externus, GP int.: globus pallidus internus, HN: habenular nuclei, Hypo: hypothalamus, NAc: nucleus accumbens, PPN: parabrachial pigmented nucleus, PUT: putamen, RN: red nucleus, SNC: substantia nigra pars compacta, SNR: substantia nigra pars reticulata, VP: ventral pallidum, VT: ventral tegmentum.
Fig. 3
Relationship between [11C]-(+2-PHNO BPND and interpolated mRNA data (log2) for A) the D2 receptor (t = -1.79, df = 12, r = -0.46, p = 0.1) and B) the D3 receptor (t = 3.30, df = 12, r = 0.69, p = 0.006). Regions-of-interest were labeled according to Pauli et al. (Pauli et al., 2018).
Legend: ROI: region-of-interest, Pauli: Pauli atlas, CAU: caudate nucleus, EA: extended amygdala, GP ext.: globus pallidus externus, GP int.: globus pallidus internus, HN: habenular nuclei, Hypo: hypothalamus, NAc: nucleus accumbens, PPN: parabrachial pigmented nucleus, PUT: putamen, RN: red nucleus, SNC: substantia nigra pars compacta, SNR: substantia nigra pars reticulata, VP: ventral pallidum, VT: ventral tegmentum.
Discussion
In this study, we integrated high-resolution gene expression maps from the AHBA open-source database with specific molecular information on dopamine receptor binding acquired with PET imaging. Our main finding is that [11C]-(+)-PHNO BPND values were strongly correlated with D3 receptor mRNA expression levels. In contrast, we failed to observe a significant correlation for the D2 receptor. In a next step, we intended to mimic the mixed D2/3 receptor signal of [11C]-(+)-PHNO by calculating the D3 mRNA ratio from the sum of DRD3 and DRD2 expression values. Likewise, the significant correlations reflect a linear relationship between mRNA expression and agonist radioligand binding that exists with D3 but not D2 receptors. Displaying comparable invitro affinity for D2
high and D3 (Seeman et al., 1993), [11C]-(+)-PHNO, in vivo, provides a mixed signal composed of binding to D2 receptors in high affinity states and D3 receptors (Graff-Guerrero et al., 2010; Rabiner et al., 2009; Searle et al., 2010; Tziortzi et al., 2011). While strong correlations with D3 mRNA levels indicate a direct relationship between D3 protein expression and [11C]-(+)-PHNO BPND values, the failure to observe such a relationship with D 2 mRNA levels supports the notion that G-protein-dependent variations in affinity states are more relevant for D2 than for D3 receptors (Sokoloff et al., 1992). Although alternative explanations such as posttranslational modifications, trafficking processes or a simple lack of power need to be considered, our data highlight the relevance of the two-state model for dopamine receptor functioning in the living human brain. Thereby, the affinity states of D2 receptors depend on G-protein coupling (George et al., 1985), which explains the weak correlation for [11C]-(+)-PHNO BPND and overall D2 gene expression including a large proportion of D2
low receptors.In contrast to previous studies evaluating the role of mRNA expression patterns on PET binding potentials for serotonergic targets (Komorowski et al., 2017; Veronese et al., 2016), we now analyzed binding of a radiotracer that is well established in imaging of the dopamine system. Thereby, the role of D2 and D3 receptors in neuropsychiatric disorders and alterations of dopamine signaling, which is modified by cortical and other upstream regions, can be assessed by PET imaging using [11C]-(+)-PHNO. Our study shows how open-source mRNA expression data can be utilized to gain further information about specific protein targets in humans, surrogating the in vivo topological distribution. However, a closer investigation in regard to radiotracer binding and the molecular effects of gene variants is necessary to fully understand regional protein biosynthesis of dopamine receptors.In general, differing correlation coefficients between both dopamine receptors might arise from epitranscriptomic modifications or alternating in vivo conformational states (Mauer and Jaffrey, 2018). E.g. analyses for D2 receptors, which exist in high and low affinity states, were complicated by their varying affinity to [11C]-(+)-PHNO. To overcome this methodological limitation, human proteomic data comprising functional genetic variants as well as information on transcriptional, translational and posttranslational mechanisms are needed. Yet, when considering potential epitranscriptomic processes of dopaminergic receptors (Araki et al., 2007; Flagel et al., 2016), the noncongruent findings between both analyzed dopamine receptors indeed appear to originate from receptor-specific binding preferences of [11C]-(+)-PHNO rather than differing mRNA modifications. The results of this study are, however, hampered by the fact that post-mortem mRNA expression data from the AHBA were obtained in a small number of subjects including limited transcriptome information from only one female donor and one hemisphere in four out of six sampled brains (Hawrylycz et al., 2012). Upon creation of the atlas, ubiquitous coverage of gene expression samples across the entire brain was prioritized over a balanced donor distribution. Although separate association analyses with regard to sex, age, or ancestry were not possible for mRNA data, these disadvantages appear negligible. Notably, the relationship between male and female [11C]-(+)PHNO BPND and mRNA expression data didn’t show significant differences when analyzed separately. Regarding other potential confounders, Matuskey et al. have described a decline of the densities of dopamine receptors with age (Matuskey et al., 2016). While this predominantly affected D2 receptor availability rather than D3 receptors, other authors did not oberve age-associated changes in the dopamine system (Nakajima et al., 2015). Even though only participants aged from 22 to 31 were included in this study, the fact that we didn’t find a relationship between age and [11C]-(+)PHNO BPND indicated a rather insignificant role of age for subsequent correlation analyses.Gene expression values based on microarray data often suffer from background and cross-hybridization disturbances (Maier et al., 2009), complicating the assumption of comparable absolute mRNA values between both dopamine receptors. Still, possible differences in absolute gene expression values between both receptors didn’t impact on the respective correlation analyses with [11C]-(+)-PHNO BPND that were performed for each receptor separately. Beside these limitations, we still regard the AHBA as an extremely valuable data source for comprehensive neuroimaging studies. To compensate for the incomplete gene expression coverage of the brain due to sparse sampling of the AHBA within several brain regions, we have analyzed interpolated gene expression data provided by Gryglewski et al. (Gryglewski et al., 2018) in addition to the original data. Interpolated expression values of dopaminergic receptors were thereby approximated with a high spatial resolution across the whole brain. Aside from that, there is still ongoing controversy regarding the selection of probes included within the AHBA. Each method, e.g. performing separate analyses for each probe, or averaging all available probes, has its advantages and disadvantages (Arnatkevicĭūė et al., 2019; Komorowski et al., 2017; Veronese et al., 2016). By selecting the median gene expression value within each brain region, we minimized possible bias on our analysis introduced by insensitive probes from the AHBA, yielding representative gene expression values within all ROIs. Considering the interpolated data by Gryglewski et al. (Gryglewski et al., 2018), probes were selected by means of a step-wise variogram modeling approach, eliminating the least sensitive probes before calculating a mean value for each mRNA sample.In sum, we elucidated the relationship between gene expression and the mixed imaging signal of [11C]-(+)-PHNO. The results of this study also indicate unique transcriptional processes for different brain regions, leading to characteristic distribution patterns of the D2 and D3 dopamine receptor with a high variability between brain regions. Notably, neuroimaging studies significantly benefit from gene expression information, especially for targets that lack sufficiently specific radiotracers in humans. Furthering this knowledge might enhance the definition of new molecular targets for the treatment of neuropsychiatric disorders.
Conclusion
This study supports evidence of differing binding properties of [11C]-(+)-PHNO to D2 and D3 dopamine receptors in different brain regions. According to our data, regional dopamine D3 receptor mRNA expression is a more relevant determinant of [11C]-(+)-PHNO binding than expression of the D2 receptor, whose affinity states depend to a greater degree on G-protein coupling. Overall, combining in vivo imaging and post-mortem gene expression data contributes to the understanding of the [11C]-(+)-PHNO PET signal from a new perspective.
Authors: David Matuskey; Patrick Worhunksy; Elizabeth Correa; Brian Pittman; Jean-Dominique Gallezot; Nabeel Nabulsi; Jim Ropchan; Venkatesh Sreeram; Rohit Gudepu; Edward Gaiser; Kelly Cosgrove; Yu-Shin Ding; Marc N Potenza; Yiyun Huang; Robert T Malison; Richard E Carson Journal: Neuroimage Date: 2016-02-11 Impact factor: 6.556
Authors: Romina Mizrahi; Ofer Agid; Carol Borlido; Ivonne Suridjan; Pablo Rusjan; Sylvain Houle; Gary Remington; Alan A Wilson; Shitij Kapur Journal: Schizophr Res Date: 2011-09 Impact factor: 4.939
Authors: Andri C Tziortzi; Graham E Searle; Sofia Tzimopoulou; Cristian Salinas; John D Beaver; Mark Jenkinson; Marc Laruelle; Eugenii A Rabiner; Roger N Gunn Journal: Neuroimage Date: 2010-06-30 Impact factor: 6.556
Authors: P Sokoloff; M Andrieux; R Besançon; C Pilon; M P Martres; B Giros; J C Schwartz Journal: Eur J Pharmacol Date: 1992-04-10 Impact factor: 4.432
Authors: Joshua B Burt; Murat Demirtaş; William J Eckner; Natasha M Navejar; Jie Lisa Ji; William J Martin; Alberto Bernacchia; Alan Anticevic; John D Murray Journal: Nat Neurosci Date: 2018-08-06 Impact factor: 24.884
Authors: Béla Kiss; István Laszlovszky; Balázs Krámos; András Visegrády; Amrita Bobok; György Lévay; Balázs Lendvai; Viktor Román Journal: Biomolecules Date: 2021-01-14