Impact of image-space resolution modeling for studies with the high-resolution research tomograph

Reconstruction of multi-animal PET acquisitions with anisotropically variant PSF

Among other factors such as random, attenuation and scatter corrections, uniform spatial resolution is key to performing accurate quantitative studies in Positron emission tomography (PET). Particularly in preclinical PET studies involving simultaneous acquisition of multiple animals, the degradation of image resolution due to the depth of interaction (DOI) effect far from the center of the Field of View (FOV) becomes a significant concern. In this work, we incorporated a spatially-variant resolution model into a real time iterative reconstruction code to obtain accurate images of multi-animal acquisition. We estimated the spatially variant point spread function (SV-PSF) across the FOV using measurements and Monte Carlo (MC) simulations. The SV-PSF obtained was implemented in a GPU-based Ordered subset expectation maximization (OSEM) reconstruction code, which includes scatter, attenuation and random corrections. The method was evaluated with acquisitions from two preclinical PET/CT scanners of the SEDECAL Argus family: a Derenzo phantom placed 2 cm off center in the 4R-SuperArgus, and a multi-animal study with 4 mice in the 6R-SuperArgus. The SV-PSF reconstructions showed uniform spatial resolution without significant increase in reconstruction time, with superior image quality compared to the uniform PSF model.

AI for PET image reconstruction

Image reconstruction for positron emission tomography (PET) has been developed over many decades, with advances coming from improved modelling of the data statistics and improved modelling of the imaging physics. However, high noise and limited spatial resolution have remained issues in PET imaging, and state-of-the-art PET reconstruction has started to exploit other medical imaging modalities (such as MRI) to assist in noise reduction and enhancement of PET's spatial resolution. Nonetheless, there is an ongoing drive towards not only improving image quality, but also reducing the injected radiation dose and reducing scanning times. While the arrival of new PET scanners (such as total body PET) is helping, there is always a need to improve reconstructed image quality due to the time and count limited imaging conditions. Artificial intelligence (AI) methods are now at the frontier of research for PET image reconstruction. While AI can learn the imaging physics as well as the noise in the data (when given sufficient examples), one of the most common uses of AI arises from exploiting databases of high-quality reference examples, to provide advanced noise compensation and resolution recovery. There are three main AI reconstruction approaches: (i) direct data-driven AI methods which rely on supervised learning from reference data, (ii) iterative (unrolled) methods which combine our physics and statistical models with AI learning from data, and (iii) methods which exploit AI with our known models, but crucially can offer benefits even in the absence of any example training data whatsoever. This article reviews these methods, considering opportunities and challenges of AI for PET reconstruction.

Effects of escitalopram on synaptic density in the healthy human brain: a randomized controlled trial

Selective serotonin reuptake inhibitors (SSRIs) are widely used for treating neuropsychiatric disorders. However, the exact mechanism of action and why effects can take several weeks to manifest is not clear. The hypothesis of neuroplasticity is supported by preclinical studies, but the evidence in humans is limited. Here, we investigate the effects of the SSRI escitalopram on presynaptic density as a proxy for synaptic plasticity. In a double-blind placebo-controlled study (NCT04239339), 32 healthy participants with no history of psychiatric or cognitive disorders were randomized to receive daily oral dosing of either 20 mg escitalopram (n = 17) or a placebo (n = 15). After an intervention period of 3-5 weeks, participants underwent a [11C]UCB-J PET scan (29 with full arterial input function) to quantify synaptic vesicle glycoprotein 2A (SV2A) density in the hippocampus and the neocortex. Whereas we find no statistically significant group difference in SV2A binding after an average of 29 (range: 24-38) days of intervention, our secondary analyses show a time-dependent effect of escitalopram on cerebral SV2A binding with positive associations between [11C]UCB-J binding and duration of escitalopram intervention. Our findings suggest that brain synaptic plasticity evolves over 3-5 weeks in healthy humans following daily intake of escitalopram. This is the first in vivo evidence to support the hypothesis of neuroplasticity as a mechanism of action for SSRIs in humans and it offers a plausible biological explanation for the delayed treatment response commonly observed in patients treated with SSRIs. While replication is warranted, these results have important implications for the design of future clinical studies investigating the neurobiological effects of SSRIs.

Genetic contributions to brain serotonin transporter levels in healthy adults

The serotonin transporter (5-HTT) critically shapes serotonin neurotransmission by regulating extracellular brain serotonin levels; it remains unclear to what extent 5-HTT levels in the human brain are genetically determined. Here we applied [11C]DASB positron emission tomography to image brain 5-HTT levels and evaluated associations with five common serotonin-related genetic variants that might indirectly regulate 5-HTT levels (BDNF rs6265, SLC6A4 5-HTTLPR, HTR1A rs6295, HTR2A rs7333412, and MAOA rs1137070) in 140 healthy volunteers. In addition, we explored whether these variants could predict in vivo 5-HTT levels using a five-fold cross-validation random forest framework. MAOA rs1137070 T-carriers showed significantly higher brain 5-HTT levels compared to C-homozygotes (2-11% across caudate, putamen, midbrain, thalamus, hippocampus, amygdala and neocortex). We did not observe significant associations for the HTR1A rs6295 and HTR2A rs7333412 genotypes. Our previously observed lower subcortical 5-HTT availability for rs6265 met-carriers remained in the presence of these additional variants. Despite this significant association, our prediction models showed that genotype moderately improved prediction of 5-HTT in caudate, but effects were not statistically significant after adjustment for multiple comparisons. Our observations provide additional evidence that serotonin-related genetic variants modulate adult human brain serotonin neurotransmission.

Kinetic models for estimating occupancy from single-scan PET displacement studies

The traditional design of PET target engagement studies is based on a baseline scan and one or more scans after drug administration. We here evaluate an alternative design in which the drug is administered during an on-going scan (i.e., a displacement study). This approach results both in lower radiation exposure and lower costs. Existing kinetic models assume steady state. This condition is not present during a drug displacement and consequently, our aim here was to develop kinetic models for analysing PET displacement data. We modified existing compartment models to accommodate a time-variant increase in occupancy following the pharmacological in-scan intervention. Since this implies the use of differential equations that cannot be solved analytically, we developed instead one approximate and one numerical solution. Through simulations, we show that if the occupancy is relatively high, it can be estimated without bias and with good accuracy. The models were applied to PET data from six pigs where [11C]UCB-J was displaced by intravenous brivaracetam. The dose-occupancy relationship estimated from these scans showed good agreement with occupancies calculated with Lassen plot applied to baseline-block scans of two pigs. In summary, the proposed models provide a framework to determine target occupancy from a single displacement scan.

Impact of list-mode reconstruction and image-space point spread function correction on PET image contrast and quantitative value using SiPM-based PET/CT system

We evaluate the effects of list-mode reconstruction and the image-space point spread function (iPSF) on the contrast and quantitative values of positron emission tomography (PET) images using a SiPM-PET/CT system. The evaluation is conducted on an NEMA body phantom and clinical images using a Cartesion Prime SiPM-PET/CT system. The signal-to-background ratio (SBR) of the phantom is set to 2, 4, 6, and 8, and all the PET image data are obtained and reconstructed using 3D-OSEM, time-of-flight, iPSF (-/ +), and a 4-mm Gaussian filter with several iterations. The evaluation criteria include % background variability (NB,10 mm), % contrast (QH,10 mm), iPSF change in QH,10 mm (ΔQH,10 mm) for edge artifact evaluation, profile curves, visual evaluation of edge artifacts, clinical imaging for the standardized uptake value (SUV) of lung nodules, and SNRliver. NB,10 mm demonstrates no significant difference in all SBRs with and without iPSF, whereas QH,10 mm is higher based on the SBR with and without iPSF. ΔQH,10 mm indicates increased iterations and a larger rate of change (> 5%) for small spheres of < 17 mm. The profile curves portrayed almost real concentrations, except for the 10-mm sphere of SBR2 without iPSF; however, with iPSF, an overshoot was observed in the 13-mm sphere of all SBRs. The degree of overshoot increased with increasing iteration and SBR. Edge artifacts were detected at values ≥ 17-22 mm in SBRs other than SBR2 with iPSF. Irrespective of the nodal size, SUV and SNRliver improved considerably after iPSF adjustment. Therefore, the effects of list-mode reconstruction and iPSF on PET image contrast were limited, and the overcorrection of the quantitative values was validated using iPSF.

Repeated low doses of psilocybin increase resilience to stress, lower compulsive actions, and strengthen cortical connections to the paraventricular thalamic nucleus in rats

Psilocybin (a classic serotonergic psychedelic drug) has received appraisal for use in psychedelic-assisted therapy of several psychiatric disorders. A less explored topic concerns the use of repeated low doses of psychedelics, at a dose that is well below the psychedelic dose used in psychedelic-assisted therapy and often referred to as microdosing. Psilocybin microdose users frequently report increases in mental health, yet such reports are often highly biased and vulnerable to placebo effects. Here we establish and validate a psilocybin microdose-like regimen in rats with repeated low doses of psilocybin administration at a dose derived from occupancy at rat brain 5-HT2A receptors in vivo. The rats tolerated the repeated low doses of psilocybin well and did not manifest signs of anhedonia, anxiety, or altered locomotor activity. There were no deficits in pre-pulse inhibition of the startle reflex, nor did the treatment downregulate or desensitize the 5-HT2A receptors. However, the repeated low doses of psilocybin imparted resilience against the stress of multiple subcutaneous injections, and reduced the frequency of self-grooming, a proxy for human compulsive actions, while also increasing 5-HT7 receptor expression and synaptic density in the paraventricular nucleus of the thalamus. These results establish a well-validated regimen for further experiments probing the effects of repeated low doses of psilocybin. Results further substantiate anecdotal reports of the benefits of psilocybin microdosing as a therapeutic intervention, while pointing to a possible physiological mechanism.

Could tau-PET imaging contribute to a better understanding of the different patterns of clinical progression in Alzheimer's disease? A 2-year longitudinal study

BACKGROUND

Monitoring the progression of Tau pathology makes it possible to study the clinical diversity of Alzheimer's disease. In this 2-year longitudinal PET study, we aimed to determine the progression of [¹⁸F]-flortaucipir binding and of cortical atrophy, and their relationships with cognitive decline.

METHODS

Twenty-seven AD patients at the mild cognitive impairment/mild dementia stages and twelve amyloid-negative controls underwent a neuropsychological assessment, 3 T brain MRI, and [¹⁸F]-flortaucipir PET imaging (Tau1) and were monitored annually over 2 years with a second brain MRI and tau-PET imaging after 2 years (Tau2). We analyzed the progression of tau standardized uptake value ratio (SUVr) and grey matter atrophy both at the regional and voxelwise levels. We used mixed effects models to explore the relations between the progression of SUVr values, cortical atrophy, and cognitive decline.

RESULTS

We found an average longitudinal increase in tau SUVr values, except for the lateral temporoparietal cortex where the average SUVr values decreased. Individual analyses revealed distinct profiles of SUVr progression according to temporoparietal Tau1 uptake: high-Tau1 patients demonstrated an increase in SUVr values over time in the frontal lobe, but a decrease in the temporoparietal cortex and a rapid clinical decline, while low-Tau1 patients displayed an increase in SUVr values in all cortical regions and a slower clinical decline. Cognitive decline was strongly associated with the progression of regional cortical atrophy, but only weakly associated with SUVr progression.

CONCLUSIONS

Despite a relatively small sample size, our results suggest that tau-PET imaging could identify patients with a potentially "more aggressive" clinical course characterized by high temporoparietal Tau1 SUVr values and a rapid clinical progression. In these patients, the paradoxical decrease in temporoparietal SUVr values over time could be due to the rapid transition to ghost tangles, for which the affinity of the radiotracer is lower. They could particularly benefit from future therapeutic trials, the neuroimaging outcome measures of which deserve to be discussed.

Evaluation of the α-synuclein PET radiotracer (d3)-[11C]MODAG-001 in pigs

BACKGROUND

A positron emission tomography (PET) radiotracer to neuroimage α-synuclein aggregates would be a crucial addition for early diagnosis and treatment development in disorders such as Parkinson's disease, where elevated aggregate levels are a histopathological hallmark. The radiotracer (d₃)-[¹¹C]MODAG-001 has recently shown promise for visualization of α-synuclein pre-formed fibrils (α-PFF) in rodents. We here test the radiotracer in a pig model where proteins are intracerebrally injected immediately before scanning. Four pigs were injected in one hemisphere with 150 μg α-PFF, and in the other hemisphere, either 75 μg α-PFF or human brain homogenate from either dementia with Lewy bodies (DLB) or Alzheimer's disease (AD) was injected. All pigs underwent one or two (d₃)-[¹¹C]MODAG-001 PET scans, quantified with the non-invasive Logan graphical analysis using the occipital cortex as a reference region.

RESULTS

The α-PFF and AD homogenate injected brain regions had high uptake of (d₃)-[¹¹C]MODAG-001 compared to the occipital cortex or cerebellum. BP_ND values in 150 μg α-PFF injected regions was 0.78, and in the AD homogenate injected regions was 0.73. By contrast, the DLB homogenate injected region did not differ in uptake and clearance compared to the reference regions. The time-activity curves and BP_ND values in the 150 μg and 75 μg injected regions of α-PFFs show a dose-dependent effect, and the PET signal could be blocked by pretreatment with unlabeled MODAG-001.

CONCLUSION

We find that both α-PFF and AD brain homogenates give rise to increased binding of (d₃)-[¹¹C]MODAG-001 when injected into the pig brain. Despite its limited specificity for cerebral α-synuclein pathology, (d₃)-[¹¹C]MODAG-001 shows promise as a lead tracer for future radiotracer development.

Evaluation of in vivo staging of amyloid deposition in cognitively unimpaired elderly aged 78-94

Amyloid-beta (Aβ) deposition is common in cognitively unimpaired (CU) elderly >85 years. This study investigated amyloid distribution and evaluated three published in vivo amyloid-PET staging schemes from a cognitively unimpaired (CU) cohort aged 84.9 ± 4.3 years (n = 75). SUV-based principal component analysis (PCA) was applied to 18F-flutemetamol PET data to determine an unbiased regional covariance pattern of tracer uptake across grey matter regions. PET staging schemes were applied to the data and compared to the PCA output. Concentration of p-tau181 was measured in blood plasma. The PCA revealed three distinct components accounting for 91.2% of total SUV variance. PC1 driven by the large common variance of uptake in neocortical and striatal regions was significantly positively correlated with global SUVRs, APOE4 status and p-tau181 concentration. PC2 represented mainly non-specific uptake in typical amyloid-PET reference regions, and PC3 the occipital lobe. Application of the staging schemes demonstrated that the majority of the CU cohort (up to 93%) were classified as having pathological amount and distribution of Aβ. Good correspondence existed between binary (+/-) classification and later amyloid stages, however, substantial differences existed between schemes for low stages with 8-17% of individuals being unstageable, i.e., not following the sequential progression of Aβ deposition. In spite of the difference in staging outcomes there was broad spatial overlap between earlier stages and PC1, most prominently in default mode network regions. This study critically evaluated the utility of in vivo amyloid staging from a single PET scan in CU elderly and found that early amyloid stages could not be consistently classified. The majority of the cohort had pathological Aβ, thus, it remains an open topic what constitutes abnormal brain Aβ in the oldest-old and what is the best method to determine that.

Evaluation of applying space-variant resolution modeling to attenuation correction in PET

Attenuation correction aims to recover the underestimated tracer uptake and improve the image contrast recovery in positron emission tomography (PET). However, traditional ray-tracing-based projection of attenuation maps is inaccurate as some physical effects are not considered, such as finite crystal size, inter-crystal penetration and inter-crystal scatter. In this study, we evaluated the effects of applying resolution modeling (RM) to attenuation correction by implementing space-variant RM to complement physical effects which are usually omitted in the traditional projection model. We verified this method on a brain PET scanner developed by our group, in both Monte Carlo simulation and real-world data, in comparison with space-invariant Gaussian RM, average-depth-of-interaction, and multi-ray tracing methods. The results indicate that the space-variant RM is superior in terms of artifacts reduction and contrast recovery.

Tau-PET imaging predicts cognitive decline and brain atrophy progression in early Alzheimer's disease

OBJECTIVES

To explore whether regional tau binding measured at baseline is associated with the rapidity of Alzheimer's disease (AD) progression over 2 years, as assessed by the decline in specified cognitive domains, and the progression of regional brain atrophy, in comparison with amyloid-positron emission tomography (PET), MRI and cerebrospinal fluid (CSF) biomarkers.

METHODS

Thirty-six patients with AD (positive CSF biomarkers and amyloid-PET) and 15 controls underwent a complete neuropsychological assessment, 3T brain MRI, [¹¹C]-PiB and [¹⁸F]-flortaucipir PET imaging, and were monitored annually over 2 years, with a second brain MRI after 2 years. We used mixed effects models to explore the relations between tau-PET, amyloid-PET, CSF biomarkers and MRI at baseline and cognitive decline and the progression of brain atrophy over 2 years in patients with AD.

RESULTS

Baseline tau-PET was strongly associated with the subsequent cognitive decline in regions that are usually associated with each cognitive domain. No significant relationship was observed between the cognitive decline and initial amyloid load, regional cortical atrophy or CSF biomarkers. Baseline tau tracer binding in the superior temporal gyrus was associated with subsequent atrophy in an inferomedial temporal volume of interest, as was the voxelwise tau tracer binding with subsequent cortical atrophy in the superior temporal, parietal and frontal association cortices.

CONCLUSIONS

These results suggest that tau tracer binding is predictive of cognitive decline in AD in domain-specific brain areas, which provides important insights into the interaction between tau burden and neurodegeneration, and is of the utmost importance to develop new prognostic markers that will help improve the design of therapeutic trials.

An in vivo Pig Model for Testing Novel Positron Emission Tomography Radioligands Targeting Cerebral Protein Aggregates

Positron emission tomography (PET) has become an essential clinical tool for diagnosing neurodegenerative diseases with abnormal accumulation of proteins like amyloid-β or tau. Despite many attempts, it has not been possible to develop an appropriate radioligand for imaging aggregated α-synuclein in the brain for diagnosing, e.g., Parkinson’s Disease. Access to a large animal model with α-synuclein pathology would critically enable a more translationally appropriate evaluation of novel radioligands. We here establish a pig model with cerebral injections of α-synuclein preformed fibrils or brain homogenate from postmortem human brain tissue from individuals with Alzheimer’s disease (AD) or dementia with Lewy body (DLB) into the pig’s brain, using minimally invasive surgery and validated against saline injections. In the absence of a suitable α-synuclein radioligand, we validated the model with the unselective amyloid-β tracer [¹¹C]PIB, which has a high affinity for β-sheet structures in aggregates. Gadolinium-enhanced MRI confirmed that the blood-brain barrier was intact. A few hours post-injection, pigs were PET scanned with [¹¹C]PIB. Quantification was done with Logan invasive graphical analysis and simplified reference tissue model 2 using the occipital cortex as a reference region. After the scan, we retrieved the brains to confirm successful injection using autoradiography and immunohistochemistry. We found four times higher [¹¹C]PIB uptake in AD-homogenate-injected regions and two times higher uptake in regions injected with α-synuclein-preformed-fibrils compared to saline. The [¹¹C]PIB uptake was the same in non-injected (occipital cortex, cerebellum) and injected (DLB-homogenate, saline) regions. With its large brain and ability to undergo repeated PET scans as well as neurosurgical procedures, the pig provides a robust, cost-effective, and good translational model for assessment of novel radioligands including, but not limited to, proteinopathies.

Brain partial volume correction with point spreading function reconstruction in high-resolution digital PET: comparison with an MR-based method in FDG imaging

OBJECTIVE

In quantitative positron emission tomography (PET) of the brain, partial volume effect due mainly to the finite spatial resolution of the PET scanner (> 3 mm full width at half maximum [FWHM]) is a primary source of error in the measurement of tracer uptake, especially in small structures such as the cerebral cortex (typically < 3 mm thickness). The aim of this study was to evaluate the partial volume correction (PVC) performance of point spread function-incorporated reconstruction (PSF reconstruction) in combination with the latest digital PET scanner. This evaluation was performed through direct comparisons with magnetic resonance imaging (MR)-based PVC (used as a reference method) in a human brain study.

METHODS

Ten healthy subjects underwent brain ¹⁸F-FDG PET (30-min acquisition) on a digital PET/CT system (Siemens Biograph Vision, 3.5-mm FWHM scanner resolution at the center of the field of view) and anatomical T1-weighted MR imaging for MR-based PVC. PSF reconstruction was applied with a wide range of iterations (4 to 256; 5 subsets). FDG uptake in the cerebral cortex was evaluated using the standardized uptake value ratio (SUVR) and compared between PSF reconstruction and MR-based PVC.

RESULTS

Cortical structures were visualized by PSF reconstruction with several tens of iterations and were anatomically well matched with the MR-derived cortical segments. Higher numbers of iterations resulted in higher cortical SUVRs, which approached those of MR-based PVC (1.76), although even with the maximum number of iterations they were still smaller by 16% (1.47), corresponding to approximately 1.5-mm FWHM of the effective spatial resolution.

CONCLUSION

With the latest digital PET scanner, PSF reconstruction can be used as a PVC technique in brain PET, albeit with suboptimal resolution recovery. A relative advantage of PSF reconstruction is that it can be applied not only to cerebral cortical regions, but also to various small structures such as small brain nuclei that are hardly visualized on anatomical T1-weighted imaging, and thus hardly recovered by MR-based PVC.

Distinct amyloid and tau PET signatures are associated with diverging clinical and imaging trajectories in patients with amnestic syndrome of the hippocampal type

We aimed to investigate the amyloid and tau PET imaging signatures of patients with amnestic syndrome of the hippocampal type (ASHT) and study their clinical and imaging progression according to their initial PET imaging status. Thirty-six patients with a progressive ASHT and 30 controls underwent a complete neuropsychological assessment, 3 T brain MRI, [¹¹C]-PiB and [¹⁸F]-Flortaucipir PET imaging. Subjects were clinically followed-up annually over 2 years, with a second 3 T MRI (n = 27 ASHT patients, n = 28 controls) and tau-PET (n = 20 ASHT patients) at the last visit. At baseline, in accordance with the recent biological definition of Alzheimer's disease (AD), the AD PET signature was defined as the combination of (i) positive cortical amyloid load, and (ii) increased tau tracer binding in the entorhinal cortices and at least one of the following regions: amygdala, parahippocampal gyri, fusiform gyri. Patients who did not meet these criteria were considered to have a non-AD pathology (SNAP). Twenty-one patients were classified as AD and 15 as SNAP. We found a circumscribed tau tracer retention in the entorhinal cortices and/or amygdala in 5 amyloid-negative SNAP patients. At baseline, the SNAP patients were older and had lower ApoE ε4 allele frequency than the AD patients, but both groups did not differ regarding the neuropsychological testing and medial temporal lobe atrophy. During the 2-year follow-up, the episodic memory and language decline, as well as the temporo-parietal atrophy progression, were more pronounced in the AD sub-group, while the SNAP patients had a more pronounced progression of atrophy in the frontal lobes. Longitudinal tau tracer binding increased in AD patients but remained stable in SNAP patients. At baseline, distinct amyloid and tau PET signatures differentiated early AD and SNAP patients despite identical cognitive profiles characterized by an isolated ASHT and a similar degree of medial temporal atrophy. During the longitudinal follow-up, AD and SNAP patients diverged regarding clinical and imaging progression. Among SNAP patients, tau PET imaging could detect a tauopathy restricted to the medial temporal lobes, which was possibly explained by primary age-related tauopathy.

False positive rates in positron emission tomography (PET) voxelwise analyses

Issues with inflated false positive rates (FPRs) in brain imaging have recently received significant attention. However, to what extent FPRs present a problem for voxelwise analyses of Positron Emission Tomography (PET) data remains unknown. In this work, we evaluate the FPR using real PET data under group assignments that should yield no significant results after correcting for multiple comparisons. We used data from 159 healthy participants, imaged with the serotonin transporter ([¹¹C]DASB; N = 100) or the 5-HT₄ receptor ([¹¹C]SB207145; N = 59). Using this null data, we estimated the FPR by performing 1,000 group analyses with randomly assigned groups of either 10 or 20, for each tracer, and corrected for multiple comparisons using parametric Monte Carlo simulations (MCZ) or non-parametric permutation testing. Our analyses show that for group sizes of 10 or 20, the FPR for both tracers was 5-99% using MCZ, much higher than the expected 5%. This was caused by a heavier-than-Gaussian spatial autocorrelation, violating the parametric assumptions. Permutation correctly controlled the FPR in all cases. In conclusion, either a conservative cluster forming threshold and high smoothing levels, or a non-parametric correction for multiple comparisons should be performed in voxelwise analyses of brain PET data.

Imaging HDACs In Vivo: Cross-Validation of the [11C]Martinostat Radioligand in the Pig Brain

PURPOSE

With the emerging knowledge about the impact of epigenetic alterations on behavior and brain disorders, the ability to measure epigenetic alterations in brain tissue in vivo has become critically important. We present the first in vivo/in vitro cross-validation of the novel positron emission tomography (PET) radioligand [¹¹C]Martinostat in the pig brain with regard to its ability to measure histone deacetylase 1-3 (HDAC1-3) levels in vivo.

PROCEDURES

Nine female Danish landrace pigs underwent 121-min dynamic PET scans with [¹¹C]Martinostat. We quantified [¹¹C]Martinostat uptake using both a simple ratio method and kinetic models with arterial input function. By the end of the scan, the animals were euthanized and the brains were extracted. We measured HDAC1-3 protein levels in frontal cortex, cerebellum vermis, and hippocampus and compared the protein levels and regional outcome values to the [¹¹C]Martinostat PET quantification.

RESULTS

[¹¹C]Martinostat distributed widely across brain regions, with the highest uptake in the cerebellum vermis and the lowest in the olfactory bulbs. Based on the Akaike information criterion, the quantification was most reliably performed by Ichise MA1 kinetic modeling, but since the radioligand displayed very slow kinetics, we also calculated standard uptake value (SUV) ratios which correlated well with V_T. The western blots revealed higher brain tissue protein levels of HDAC1/2 compared to HDAC3, and HDAC1 and HDAC2 levels were highly correlated in all three investigated brain regions. The in vivo SUV ratio measure correlated well with the in vitro HDAC1-3 levels, whereas no correlation was found between V_T values and HDAC levels.

CONCLUSIONS

We found good correlation between in vivo measured SUV ratios and in vitro measures of HDAC 1-3 proteins, supporting that [¹¹C]Martinostat provides a good in vivo measure of the cerebral HDAC1-3 protein levels.

Improved Detection of Postoperative Residual Meningioma with [68Ga]Ga-DOTA-TOC PET Imaging Using a High-resolution Research Tomograph PET Scanner

PURPOSE

PET with somatostatin receptor ligand [⁶⁸Ga]Ga-DOTA-D-Phe¹-Tyr³-octreotide ([⁶⁸Ga]Ga-DOTA-TOC) is an established method in radiotherapy planning because of the improved detection and delineation of meningioma tissue. We investigated the diagnostic accuracy of supplementary [⁶⁸Ga]Ga-DOTA-TOC PET in patients with a 3-month postoperative MRI reporting gross-total resection (GTR).

EXPERIMENTAL DESIGN

Thirty-seven patients with a histologically proven meningioma and GTR on postoperative MRI were prospectively referred to [⁶⁸Ga]Ga-DOTA-TOC PET. Detection and volume measurements of [⁶⁸Ga]Ga-DOTA-TOC-avid lesions in relation to the primary tumor site were recorded. Residual tumor in suspicious lesions suggested by [⁶⁸Ga]Ga-DOTA-TOC PET was verified by (i) tumor recurrence/progression on subsequent MRI scans according to the Response Assessment of Neuro-Oncology criteria, (ii) subsequent histology, and (iii) follow-up [⁶⁸Ga]Ga-DOTA-TOC PET scan.

RESULTS

Twenty-three PET scans demonstrated [⁶⁸Ga]Ga-DOTA-TOC-avid lesions suspicious of residual meningioma, where 18 could be verified by (i) tumor progression on subsequent MRI scans (n = 6), (ii) histologic confirmation (n = 3), and (iii) follow-up [⁶⁸Ga]Ga-DOTA-TOC PET scans confirming the initial PET findings (n = 9) after an overall median follow-up time of 17 months (range, 9-35 months). In contrast, disease recurrence was seen in only 2 of 14 patients without [⁶⁸Ga]Ga-DOTA-TOC-avid lesions (P < 0.0001). The sensitivity, specificity, and diagnostic accuracy of [⁶⁸Ga]Ga-DOTA-TOC PET in detecting meningioma residue was 90% [95% confidence interval (CI), 67-99], 92% (95% CI, 62-100), and 90% (95% CI, 74-98; P < 0.0001), respectively.

CONCLUSIONS

The majority of patients with GTR on 3-month postoperative MRI may have small unrecognized meningioma residues that can be detected using [⁶⁸Ga]Ga-DOTA-TOC PET.

A high-resolution in vivo atlas of the human brain's benzodiazepine binding site of GABAA receptors

Gamma-aminobutyric acid (GABA) is the main inhibitory neurotransmitter in the human brain and plays a key role in several brain functions and neuropsychiatric disorders such as anxiety, epilepsy, and depression. For decades, several in vivo and ex vivo techniques have been used to highlight the mechanisms of the GABA system, however, no studies have currently combined the techniques to create a high-resolution multimodal view of the GABA system.

Here, we present a quantitative high-resolution in vivo atlas of the human brain benzodiazepine receptor sites (BZR) located on postsynaptic ionotropic GABA_A receptors (GABA_A Rs), generated on the basis of in vivo [¹¹C]flumazenil Positron Emission Tomography (PET) data. Next, based on ex vivo autoradiography data, we transform the PET-generated atlas from binding values into BZR protein density. Finally, we examine the brain regional association between BZR protein density and ex vivo mRNA expression for the 19 subunits in the GABA_A R, including an estimation of the minimally required expression of mRNA levels for each subunit to translate into BZR protein.

This represents the first publicly available quantitative high-resolution in vivo atlas of the spatial distribution of BZR densities in the healthy human brain. The atlas provides a unique neuroscientific tool as well as novel insights into the association between mRNA expression for individual subunits in the GABA_A R and the BZR density at each location in the brain.

Effects of a single dose of psilocybin on behaviour, brain 5-HT2A receptor occupancy and gene expression in the pig

Psilocybin has in some studies shown promise as treatment of major depressive disorder and psilocybin therapy was in 2019 twice designated as breakthrough therapy by the U.S. Food and Drug Administration (FDA). A very particular feature is that ingestion of just a single dose of psilocybin is associated with lasting changes in personality and mood. The underlying molecular mechanism behind its effect is, however, unknown. In a translational pig model, we here present the effects of a single dose of psilocybin on pig behaviour, receptor occupancy and gene expression in the brain. An acute i.v. injection of 0.08 mg/kg psilocybin to awake female pigs induced characteristic behavioural changes in terms of headshakes, scratching and rubbing, lasting around 20 min. A similar dose was associated with a cerebral 5-HT_2A receptor occupancy of 67%, as determined by positron emission tomography, and plasma psilocin levels were comparable to what in humans is associated with an intense psychedelic experience. We found that 19 genes were differentially expressed in prefrontal cortex one day after psilocybin injection, and 3 genes after 1 week. Gene Set Enrichment Analysis demonstrated that multiple immunological pathways were regulated 1 week after psilocybin exposure. This provides a framework for future investigations of the lasting molecular mechanisms induced by a single dose of psilocybin. In the light of an ongoing debate as to whether psilocybin is a safe treatment for depression and other mental illnesses, it is reassuring that our data suggest that any effects on gene expression are very modest.

Advances in multimodal data fusion in neuroimaging: Overview, challenges, and novel orientation

Multimodal fusion in neuroimaging combines data from multiple imaging modalities to overcome the fundamental limitations of individual modalities. Neuroimaging fusion can achieve higher temporal and spatial resolution, enhance contrast, correct imaging distortions, and bridge physiological and cognitive information. In this study, we analyzed over 450 references from PubMed, Google Scholar, IEEE, ScienceDirect, Web of Science, and various sources published from 1978 to 2020. We provide a review that encompasses (1) an overview of current challenges in multimodal fusion (2) the current medical applications of fusion for specific neurological diseases, (3) strengths and limitations of available imaging modalities, (4) fundamental fusion rules, (5) fusion quality assessment methods, and (6) the applications of fusion for atlas-based segmentation and quantification. Overall, multimodal fusion shows significant benefits in clinical diagnosis and neuroscience research. Widespread education and further research amongst engineers, researchers and clinicians will benefit the field of multimodal neuroimaging.

Increased microglial activation in patients with Parkinson disease using [18F]-DPA714 TSPO PET imaging

INTRODUCTION

Increasing evidence suggests that neuroinflammation is active in Parkinson disease (PD) and contributes to neurodegeneration. This process can be studied in vivo with PET and radioligands targeting TSPO, upregulated in activated microglia. Initial PET studies investigating microglial activation in PD with the [¹¹C]-PK11195 have provided inconclusive results. Here we assess the presence and distribution of neuroinflammatory response in PD patients using [¹⁸F]-DPA714 and to correlate imaging biomarkers to dopamine transporter imaging and clinical status.

METHODS

PD patients (n = 24, Hoehn and Yahr I-III) and 28 healthy controls were scanned with [¹⁸F]-DPA714 and [¹¹C]-PE2I and analyzed. They were all genotyped for TSPO polymorphism. Regional binding parameters were estimated (reference Logan graphical approach with supervised cluster analysis). Impact of TSPO genotype was analyzed using Wilcoxon signed-rank test. Differences between groups were investigated using a two-way ANOVA and Tukey post hoc tests.

RESULTS

PD patients showed significantly higher [¹⁸F]-DPA714 binding compared to healthy controls bilaterally in the midbrain (p < 0.001), the frontal cortex (p = 0.001), and the putamen contralateral to the more clinically affected hemibody (p = 0.038). Microglial activation in these regions did not correlate with the severity of motor symptoms, disease duration nor putaminal [¹¹C]-PE2I uptake. However, there was a trend toward a correlation between cortical TSPO binding and disease duration (p = 0.015 uncorrected, p = 0.07 after Bonferroni correction).

CONCLUSION

[¹⁸F]-DPA714 binding confirmed that there is a specific topographic pattern of microglial activation in the nigro-striatal pathway and the frontal cortex of PD patients.

TRIAL REGISTRATION

Trial registration: INFLAPARK, NCT02319382. Registered 18 December 2014- Retrospectively registered, https://clinicaltrials.gov/ct2/show/NCT02319382.

Validation of diffuse correlation spectroscopy against 15O-water PET for regional cerebral blood flow measurement in neonatal piglets

Diffuse correlation spectroscopy (DCS) can non-invasively and continuously asses regional cerebral blood flow (rCBF) at the cot-side by measuring a blood flow index (BFI) in non-traditional units of cm2/s. We have validated DCS against positron emission tomography using 15O-labeled water (15O-water PET) in a piglet model allowing us to derive a conversion formula for BFI to rCBF in conventional units (ml/100g/min). Neonatal piglets were continuously monitored by the BabyLux device integrating DCS and time resolved near infrared spectroscopy (TRS) while acquiring 15O-water PET scans at baseline, after injection of acetazolamide and during induced hypoxic episodes. BFI by DCS was highly correlated with rCBF (R = 0.94, p < 0.001) by PET. A scaling factor of 0.89 (limits of agreement for individual measurement: 0.56, 1.39)×109× (ml/100g/min)/(cm2/s) was used to derive baseline rCBF from baseline BFI measurements of another group of piglets and of healthy newborn infants showing an agreement with expected values. These results pave the way towards non-invasive, cot-side absolute CBF measurements by DCS on neonates.

Extension of the ML-EM algorithm for dose estimation using PET in proton therapy: application to an inhomogeneous target

Positron emission tomography (PET) has been used for in vivo treatment verification, mainly for range verification, in proton therapy. Evaluating the direct dose from PET measurements remains challenging; however, it is highly desirable from a clinical perspective. In this study, a method for estimating the dose distribution from the positron emitter distributions was developed using the maximum likelihood expectation maximization algorithm. The 1D spatial relationship between positron emitter distributions and a dose distribution in an inhomogeneous target was inputted into the system matrix based on a filter framework. In contrast, spatial resolution of the PET system and total variation regularization (as prior knowledge for dose distribution) were considered in the 3D image-space. The dose estimation was demonstrated using Monte Carlo simulated PET activity distributions with substantial noise in a head and neck phantom. This mimicked the single field irradiation of the spread-out Bragg peak beams at clinical dose levels. Besides the simple implementation of the algorithm, this strategy achieved a high-speed calculation (30 s for a 3D dose estimation) and accurate dose and range estimations (less than 10% and 2 mm errors at 1-σ values, respectively). The proposed method could be key for using PET for in vivo dose monitoring.

Different preprocessing strategies lead to different conclusions: A [11C]DASB-PET reproducibility study

Positron emission tomography (PET) neuroimaging provides unique possibilities to study biological processes in vivo under basal and interventional conditions. For quantification of PET data, researchers commonly apply different arrays of sequential data analytic methods ("preprocessing pipeline"), but it is often unknown how the choice of preprocessing affects the final outcome. Here, we use an available data set from a double-blind, randomized, placebo-controlled [¹¹C]DASB-PET study as a case to evaluate how the choice of preprocessing affects the outcome of the study. We tested the impact of 384 commonly used preprocessing strategies on a previously reported positive association between the change from baseline in neocortical serotonin transporter binding determined with [¹¹C]DASB-PET, and change in depressive symptoms, following a pharmacological sex hormone manipulation intervention in 30 women. The two preprocessing steps that were most critical for the outcome were motion correction and kinetic modeling of the dynamic PET data. We found that 36% of the applied preprocessing strategies replicated the originally reported finding (p < 0.05). For preprocessing strategies with motion correction, the replication percentage was 72%, whereas it was 0% for strategies without motion correction. In conclusion, the choice of preprocessing strategy can have a major impact on a study outcome.

Common HTR2A variants and 5-HTTLPR are not associated with human in vivo serotonin 2A receptor levels

The serotonin 2A receptor (5‐HT2AR) is implicated in the pathophysiology and treatment of various psychiatric disorders. [¹⁸F]altanserin and [¹¹C]Cimbi‐36 positron emission tomography (PET) allow for high‐resolution imaging of 5‐HT2AR in the living human brain. Cerebral 5‐HT2AR binding is strongly genetically determined, though the impact of specific variants is poorly understood. Candidate gene studies suggest that HTR2A single nucleotide polymorphisms including rs6311/rs6313, rs6314, and rs7997012 may influence risk for psychiatric disorders and mediate treatment response. Although known to impact in vitro expression of 5‐HT2AR or other serotonin (5‐HT) proteins, their effect on human in vivo brain 5‐HT2AR binding has as of yet been insufficiently studied. We thus assessed the extent to which these variants and the commonly studied 5‐HTTLPR predict neocortex in vivo 5‐HT2AR binding in healthy adult humans. We used linear regression analyses and likelihood ratio tests in 197 subjects scanned with [¹⁸F]altanserin or [¹¹C]Cimbi‐36 PET. Although we observed genotype group differences in 5‐HT2AR binding of up to ~10%, no genetic variants were statistically significantly predictive of 5‐HT2AR binding in what is the largest human in vivo 5‐HT2AR imaging genetics study to date. Thus, in vitro and post mortem results suggesting effects on 5‐HT2AR expression did not carry over to the in vivo setting. To any extent these variants might affect clinical risk, our findings do not support that 5‐HT2AR binding mediates such effects. Our observations indicate that these individual variants do not significantly contribute to genetic load on human in vivo 5‐HT2AR binding.

Imaging the aging brain: study design and baseline findings of the SENIOR cohort

Background

Current demographic trends point towards an aging society entailing increasing occurrence and burden of neurodegenerative diseases. In this context, understanding physiological aging and its turning point into neurodegeneration is essential for the development of possible biomarkers and future therapeutics of brain disease.

Methods

The SENIOR study represents a longitudinal, observational study including cognitively healthy elderlies aged between 50 and 70 years old at the time of inclusion, being followed annually over 10 years. Our multimodal protocol includes structural, diffusion, functional, and sodium magnetic resonance imaging (MRI) at 3 T and 7 T, positron emission tomography (PET), blood samples, genetics, audiometry, and neuropsychological and neurological examinations as well as assessment of neuronal risk factors.

Results

One hundred forty-two participants (50% females) were enrolled in the SENIOR cohort with a mean age of 60 (SD 6.3) years at baseline. Baseline results with multiple regression analyses reveal that cerebral white matter lesions can be predicted by cardiovascular and cognitive risk factors and age. Cardiovascular risk factors were strongly associated with juxtacortical and periventricular lesions. Intra-subject across-test variability as a measure of neuropsychological test performance and possible cognitive marker predicts white matter volume and is significantly associated with risk profile. Division of the cohort into subjects with a higher and lower risk profile shows significant differences in intra-subject across-test variability and volumes as well as cortical thickness of brain regions of the temporal lobe. There is no difference between the lower- and higher-risk groups in amyloid load using PET data from a subset of 81 subjects.

Conclusions

We here describe the study protocol and baseline findings of the SENIOR observational study which aim is the establishment of integrated, multiparametric maps of normal aging and the identification of early biomarkers for neurodegeneration. We show that intra-subject across-test variability as a marker of neuropsychological test performance as well as age, gender, and combined risk factors influence neuronal decline as represented by decrease in brain volume, cortical thickness, and increase in white matter lesions. Baseline findings will be used as underlying basis for the further implications of aging and neuronal degeneration as well as examination of brain aging under different aspects of brain pathology versus physiological aging.

Association Between Sumatriptan Treatment During a Migraine Attack and Central 5-HT1B Receptor Binding

Importance

Triptans, the most efficient acute treatment for migraine attacks, are 5-HT1B/1D receptor agonists, but their precise mechanism of action is not completely understood. The extent to which triptans enter the central nervous system and bind to 5-HT1B receptors in the brain is unknown.

Objectives

To determine the occupancy of sumatriptan to central 5-HT1B receptors, and to investigate changes in brain serotonin levels during migraine attacks.

Design, Setting, and Participants

This study of 8 patients in Denmark used a within-participant design and was conducted from April 20, 2015, to December 5, 2016. Participants were otherwise healthy patients with untreated episodic migraine without aura, aged between 18 and 65 years, and recruited from the general community. Data analysis was performed from January 2017 to April 2018.

Interventions

All participants underwent positron emission tomographic scans after injection of [11C]AZ10419369, a specific 5-HT1B receptor radiotracer. All participants were scanned 3 times: (1) during an experimentally induced migraine attack, (2) after a subcutaneous injection of 6-mg subcutaneous sumatriptan, and (3) on a migraine attack-free day. Scans 1 and 2 were conducted on the same study day. Each scan lasted for 90 minutes.

Main Outcome and Measure

The primary outcome was the nondisplaceable binding potential of [11C]AZ10419369 across 7 brain regions involved in pain modulation. The binding potential reflects receptor density, and changes in binding potential reflects displacement of the radiotracer. The occupancy of sumatriptan was estimated from the 2 scans before and after sumatriptan administration.

Results

Eight patients with migraine were included in the study; of these participants, 7 (87%) were women. The mean (SD) age of participants on study day 1 was 29.5 (9.2) years and on study day 2 was 30.0 (8.9) years. Sumatriptan was associated with statistically significantly reduced 5-HT1B receptor binding across pain-modulating regions (mean [SD] binding potential, 1.20 [0.20] vs 1.02 [0.22]; P = .001), corresponding to a mean (SD) drug occupancy rate of 16.0% (5.3%). Furthermore, during migraine attacks, as compared with outside of attacks, 5-HT1B receptor binding was statistically significantly associated with reduced in pain-modulating regions (mean [SD] binding potential, 1.36 [0.22] vs 1.20 [0.20]; P = .02).

Conclusions and Relevance

Treatment with sumatriptan during migraine attacks appeared to be associated with a decrease in 5-HT1B receptor binding, a finding that is most likely associated with the binding of sumatriptan to central 5-HT1B receptors, but the contribution of ongoing cerebral serotonin release to the lower binding cannot be excluded; the migraine attack-associated decrease in binding could indicate that migraine attacks are associated with increases in endogenous serotonin.

Effects of TOF Resolution Models on Edge Artifacts in PET Reconstruction from Limited-Angle Data

Limited-angle data, such as data obtained from a dual-panel Breast-PET scanner, result in substantial image blur in directions coinciding with the missing cone of the image spectrum. On systems with time-of-flight (TOF) capabilities, this blur is reduced as given by the TOF uncertainty, with the image spectrum being correspondingly expanded into the missing spectral cone. Modeling of the TOF uncertainty in the reconstruction is expected to deconvolve this residual TOF blurring. We have however observed that, as a tradeoff, this TOF de-blurring process also introduces ringing artifacts at the edges, analogous to the edge effects observed with line-of-response (LOR) resolution modeling, which attempts to deconvolve the blur due to detector resolution effects. However, in the former case, the ringing artifacts are much wider due to the spatial extent of the TOF uncertainty as compared to the width of typical LOR resolution blur. We illustrate and investigate the effects of using matched, as well as under-modeled and over-modeled, TOF kernels on edge artifacts in reconstruction from limited-angle data, and compare them with TOF reconstructions of complete data. Although for the conventional data with full angular coverage the reconstruction is fairly insensitive to the exact size of the TOF kernel and TOF modeling does not produce ringing artifacts, it is not the case for the limited-angle data. We show that it is important to use some form of regularization of the TOF uncertainty deconvolution process within reconstruction of the limited-angle data, such as decreasing the TOF kernel size.

Evaluation of quantitative, efficient image reconstruction for VersaPET, a compact PET system

PURPOSE

Previously we developed a high-resolution positron emission tomography (PET) system-VersaPET-characterized by a block geometry with relatively large axial and transaxial interblock gaps and a compact geometry susceptible to parallax blurring effects. In this work, we report the qualitative and quantitative evaluation of a graphic processing unit (GPU)-accelerated maximum-likelihood by expectation-maximization (MLEM) image reconstruction framework for VersaPET which features accurate system geometry and projection space point-spread-function (PSF) modeling.

METHODS

We combined the ray-tracing module from software for tomographic image reconstruction (STIR), an open-source PET image reconstruction package, with VersaPET's exact block geometry for the geometric system matrix. Point-spread-function modeling of crystal penetration and scattering was achieved by a custom Monte-Carlo simulation for projection space blurring in all dimensions. We also parallelized the reconstruction in GPU taking advantage of the system's symmetry for PSF computation. To investigate the effects of PSF width, we generated and studied multiple kernels between one that reflects the true LYSO density in the MC simulation and another that reflects geometry only (no PSF). GATE simulations of hot and cold-sphere phantoms with spheres of different sizes, real microDerenzo phantom, and human blood vessel data were used to characterize the quantitative and qualitative performances of the reconstruction.

RESULTS

Reconstruction with an accurate system geometry effectively improved image quality compared to STIR (version 3.0) which assumes an idealized system geometry. Reconstructions of GATE-simulated hot-sphere phantom data showed that all PSF kernels achieved superior performance in contrast recovery and bias reduction compared to using no PSF, but may introduce edge artifact and lumped background noise pattern depending on the width of PSF kernels. Cold-sphere phantom simulation results also indicated improvement in contrast recovery and quantification with PSF modeling (compared to no PSF) for 5 and 10 mm cold spheres. Real microDerenzo phantom images with the PSF kernel that reflects the true LYSO density showed degraded resolving power of small sectors that could be resolved more clearly by underestimated PSF kernels, which is consistent with recent literature despite differences in scanner geometries and in approaches to system model estimation. The human vessel results resemble those of the hot-sphere phantom simulation with the PSF kernel that reflects the true LYSO density achieving the highest peak in the time activity curve (TAC) and similar lumped noise pattern.

CONCLUSIONS

We fully evaluated a practical MLEM reconstruction framework that we developed for VersaPET in terms of qualitative and quantitative performance. Different PSF kernels may be adopted for improving the results of specific imaging tasks but the underlying reasons for the variation in optimal kernel for the real and simulation studies requires further study.

Predicting Treatment Outcome in Major Depressive Disorder Using Serotonin 4 Receptor PET Brain Imaging, Functional MRI, Cognitive-, EEG-Based, and Peripheral Biomarkers: A NeuroPharm Open Label Clinical Trial Protocol

BACKGROUND

Between 30 and 50% of patients with major depressive disorder (MDD) do not respond sufficiently to antidepressant regimens. The conventional pharmacological treatments predominantly target serotonergic brain signaling but better tools to predict treatment response and identify relevant subgroups of MDD are needed to support individualized and mechanistically targeted treatment strategies. The aim of this study is to investigate antidepressant-free patients with MDD using neuroimaging, electrophysiological, molecular, cognitive, and clinical examinations and evaluate their ability to predict clinical response to SSRI treatment as individual or combined predictors.

METHODS

We will include 100 untreated patients with moderate to severe depression (>17 on the Hamilton Depression Rating Scale 17) in a non-randomized open clinical trial. We will collect data from serotonin 4 receptor positron emission tomography (PET) brain scans, functional magnetic resonance imaging (fMRI), electroencephalogram (EEG), cognitive tests, psychometry, and peripheral biomarkers, before (at baseline), during, and after 12 weeks of standard antidepressant treatment. Patients will be treated with escitalopram, and in case of non-response at week 4 or intolerable side effects, offered to switch to a second line treatment with duloxetine. Our primary outcome (treatment response) is assessed using the Hamilton depression rating subscale 6-item scores at week 8, compared to baseline. In a subset of the patients (n = ~40), we will re-assess the neurobiological response (using PET, fMRI, and EEG) 8 weeks after initiated pharmacological antidepressant treatment, to map neurobiological signatures of treatment responses. Data from matched controls will either be collected or is already available from other cohorts.

DISCUSSION

The extensive investigational program with follow-up in this large cohort of participants provides a unique possibility to (a) uncover potential biomarkers for antidepressant treatment response, (b) apply the findings for future stratification of MDD, (c) advance the understanding of pathophysiological underpinnings of MDD, and (d) uncover how putative biomarkers change in response to 8 weeks of pharmacological antidepressant treatment. Our data can pave the way for a precision medicine approach for optimized treatment of MDD and also provides a resource for future research and data sharing.

CLINICAL TRIAL REGISTRATION

The study was registered at clinicaltrials.gov prior to initiation (NCT02869035; 08.16.2016, URL: https://clinicaltrials.gov/ct2/results?cond=&term=NCT02869035&cntry=&state=&city=&dist=).

Spatially-variant image-based modeling of PSF deformations with application to a limited angle geometry from a dual-panel breast-PET imager

Dual-panel PET system configuration can lead to spatially variable point-spread functions (PSF) of considerable deformations due to depth-of-interaction effects and limited angular coverage. If not modelled properly, these effects result in decreased and inconsistent recovery of lesion activity across the field-of-view (FOV), as well as mispositioning of lesions in the reconstructed image caused by strong PSF asymmetries. We implemented and evaluated models of such PSF deformations with spatially-variant image-based resolution modeling (IRM) within reconstruction (varRM) using the Direct Image REConstruction for Time-of-flight (DIRECT) method and within post-reconstruction deconvolution methods. In addition, DIRECT reconstruction was performed with a spatially-invariant IRM (invRM) and without resolution modeling (noRM) for comparison. The methods were evaluated using simulated data for a realistic breast model with a set of 5 mm lesions located throughout the FOV of a dual-panel Breast-PET scanner. We simulated high-count data to focus on the ability of each method to correctly recover the PSF deformations, and a clinically realistic count level to assess the impact of low count data on the quantitative performance of the evaluated techniques. Performance of the methods evaluated herein was assessed by comparing lesion activity recovery (%BIAS), consistency (%SD) across the FOV, overall error (%RMSE), and recovery of each lesion location. As expected, all techniques using IRM provide considerable improvement over the noRM reconstruction. For the high-count cases, the overall quantitative performance of all IRM techniques, whether within reconstruction or within post-reconstruction, is similar if the lesion location misplacements are ignored. However, invRM provides less consistent performance on activity across lesions and is not able to recover accurate lesion locations. For a clinically realistic count level, varRM reconstruction consistently outperforms all compared approaches, while the post-reconstruction IRM approaches exhibit higher %SD and %RMSE values due to being more affected by the data noise than the within-reconstruction IRM approaches.

Improved image resolution on thoracic carcinomas by quantitative 18F-FDG coincidence SPECT/CT in comparison to 18F-FDG PET/CT

Currently, ¹⁸F-FDG coincidence SPECT (Co-SPECT)/CT scan still serves as an important tool for diagnosis, staging, and evaluation of cancer treatment in developing countries. We implemented full physical corrections (FPC) to Co-SPECT (quantitative Co-SPECT) to improve the image resolution and contrast along with the capability for image quantitation. FPC included attenuation, scatter, resolution recovery, and noise reduction. A standard NEMA phantom filled with 10:1 F-18 activity concentration ratio in spheres and background was utilized to evaluate image performance. Subsequently, 15 patients with histologically confirmed thoracic carcinomas were included to undergo a ¹⁸F-FDG Co-SPECT/CT scan followed by a ¹⁸F-FDG PET/CT scan. Functional parameters as SUVmax, SUVmean, SULpeak, and MTV from both quantitative Co-SPECT and PET were analyzed. Image resolution of Co-SPECT for NEMA phantom was improved to reveal the smallest sphere from a diameter of 28 mm to 22 mm (17 mm for PET). The image contrast was enhanced from 1.7 to 6.32 (6.69 for PET) with slightly degraded uniformity in background (3.1% vs. 6.7%) (5.6% for PET). Patients' SUVmax, SUVmean, SULpeak, and MTV measured from quantitative Co-SPECT were overall highly correlated with those from PET (r=0.82–0.88). Adjustment of the threshold of SUVmax and SUV to determine SUVmean and MTV did not further change the correlations with PET (r=0.81–0.88). Adding full physical corrections to Co-SPECT images can significantly improve image resolution and contrast to reveal smaller tumor lesions along with the capability to quantify functional parameters like PET/CT.

Implementation of cylindrical PET scanners with block detector geometry in STIR

BACKGROUND

Software for Tomographic Image Reconstruction (STIR) is an open-source library for PET and SPECT image reconstruction, implementing iterative reconstruction as well as 2D- and 3D-filtered back projection. Quantitative reconstruction of PET data requires the knowledge of the scanner geometry. Typical scanners, clinical as well as pre-clinical ones, use a block-type geometry. Several rectangular blocks of crystals are arranged into regular polygons. Multiple of such polygons are arranged along the scanner axis. However, the geometrical representation of a scanner provided by STIR is a cylinder made of rings of individual crystals equally distributed in axial and transaxial directions. The data of realistic scanners are projected onto such virtual scanners prior to image reconstruction. This results in reduced quality of the reconstructed image. In this study, we implemented the above-described block geometry into the STIR library, permitting the image reconstruction without the interpolation step. In order to evaluate the difference in image quality, we performed Monte Carlo simulation studies of three different scanner designs: two scanners with multiple crystals per block and one with a single crystal per block. Simulated data were reconstructed using the standard STIR method and the newly implemented block geometry.

RESULTS

Visual comparison between the images reconstructed by the two models for the block-type scanners shows that the new implementation enhances the image quality to the extent that the results before normalization correction are comparable with those after normalization correction. The simulation result of a uniform cylinder shows that the coefficient of variation decreases from 25.8% to 20.9% by using the new implementation in STIR. Spatial resolution is enhanced resulting in a lower partial loss of intensity in sources of small size, e.g., the spill-over ratio for spherical sources of 1.8 mm diameter is 0.19 in the block and 0.34 in the cylindrical model.

CONCLUSIONS

Results indicate a significant improvement for the new model in comparison with the old one which mapped the polygonal geometry into a cylinder. The new implementation was tested and is available for use via the library of Swiss Federal Institute of Technology in Zurich (ETH).

Resolution modeling in projection space using a factorized multi-block detector response function for PET image reconstruction

Positron emission tomography (PET) images usually suffer from limited resolution and statistical uncertainties. However, a technique known as resolution modeling (RM) can be used to improve image quality by accurately modeling the system's detection process within the iterative reconstruction. In this study, we present an accurate RM method in projection space based on a simulated multi-block detector response function (DRF) and evaluate it on the Siemens hybrid MR-BrainPET system. The DRF is obtained using GATE simulations that consider nearly all the possible annihilation photons from the field-of-view (FOV). Intrinsically, the multi-block DRF allows the block crosstalk to be modeled. The RM blurring kernel is further generated by factorizing the blurring matrix of one line-of-response (LOR) into two independent detector responses, which can then be addressed with the DRF. Such a kernel is shift-variant in 4D projection space without any distance or angle compression, and is integrated into the image reconstruction for the BrainPET insert with single instruction multiple data (SIMD) and multi-thread support. Evaluation of simulations and measured data demonstrate that the reconstruction with RM yields significantly improved resolutions and reduced mean squared error (MSE) values at different locations of the FOV, compared with reconstruction without RM. Furthermore, the shift-variant RM kernel models the varying blurring intensity for different LORs due to the depth-of-interaction (DOI) dependencies, thus avoiding severe edge artifacts in the images. Additionally, compared to RM in single-block mode, the multi-block mode shows significantly improved resolution and edge recovery at locations beyond 10 cm from the center of BrainPET insert in the transverse plane. However, the differences have been observed to be low for patient data between single-block and multi-block mode RM, due to the brain size and location as well as the geometry of the BrainPET insert. In conclusion, the RM method proposed in this study can yield better reconstructed images in terms of resolution and MSE value, compared to conventional reconstruction without RM.

Optimization of preprocessing strategies in Positron Emission Tomography (PET) neuroimaging: A [11C]DASB PET study

Positron Emission Tomography (PET) is an important neuroimaging tool to quantify the distribution of specific molecules in the brain. The quantification is based on a series of individually designed data preprocessing steps (pipeline) and an optimal preprocessing strategy is per definition associated with less noise and improved statistical power, potentially allowing for more valid neurobiological interpretations. In spite of this, it is currently unclear how to design the best preprocessing pipeline and to what extent the choice of each preprocessing step in the pipeline minimizes subject-specific errors. To evaluate the impact of various preprocessing strategies, we systematically examined 384 different pipeline strategies in data from 30 healthy participants scanned twice with the serotonin transporter (5-HTT) radioligand [¹¹C]DASB. Five commonly used preprocessing steps with two to four options were investigated: (1) motion correction (MC) (2) co-registration (3) delineation of volumes of interest (VOI's) (4) partial volume correction (PVC), and (5) kinetic modeling. To quantitatively compare and evaluate the impact of various preprocessing strategies, we used the performance metrics: test-retest bias, within- and between-subject variability, the intraclass-correlation coefficient, and global signal-to-noise ratio. We also performed a power analysis to estimate the required sample size to detect either a 5% or 10% difference in 5-HTT binding as a function of preprocessing pipeline. The results showed a complex downstream dependency between the various preprocessing steps on the performance metrics. The choice of MC had the most profound effect on 5-HTT binding, prior to the effects caused by PVC and kinetic modeling, and the effects differed across VOI's. Notably, we observed a negative bias in 5-HTT binding across test and retest in 98% of pipelines, ranging from 0 to 6% depending on the pipeline. Optimization of the performance metrics revealed a trade-off in within- and between-subject variability at the group-level with opposite effects (i.e. minimization of within-subject variability increased between-subject variability and vice versa). The sample size required to detect a given effect size was also compromised by the preprocessing strategy, resulting in up to 80% increases in sample size needed to detect a 5% difference in 5-HTT binding. This is the first study to systematically investigate and demonstrate the effect of choosing different preprocessing strategies on the outcome of dynamic PET studies. We provide a framework to show how optimal and maximally powered neuroimaging results can be obtained by choosing appropriate preprocessing strategies and we provide recommendations depending on the study design. In addition, the results contribute to a better understanding of methodological uncertainty and variability in preprocessing decisions for future group- and/or longitudinal PET studies.

Metabolic regional and network changes in Alzheimer's disease subtypes

Clinical variants of Alzheimer's disease (AD) include the common amnestic subtype as well as subtypes characterised by leading visual processing impairments or by multimodal neurocognitive deficits. We investigated regional metabolic patterns and networks between AD subtypes. The study comprised 9 age-matched controls and 25 patients with mild to moderate AD. Methods included clinical and neuropsychological assessment, high-resolution FDG PET and T1-weighted 3D MR imaging with PET-MR coregistration, grey matter segmentation, atlas-based regions-of-interest, linear mixed effects and regional correlation analysis. Regional metabolic patterns differed significantly between groups, but significant hypometabolism in the posterior cingulate cortex (PCC) was common to all subtypes. The most distinctive regional abnormality was occipital hypometabolism in the visual subtype. In controls, two large clusters of positive regional metabolic correlations were observed. The most pronounced breakdown of the normal correlation pattern was found in amnestic patients who, in contrast, showed the least regional focal metabolic deficits. The normal positive correlation between PCC and hippocampus was lost in all subtypes. In conclusion, PCC hypometabolism and metabolic correlation breakdown between PCC and hippocampus are the common functional core of all AD subtypes. Network alterations exceed focal regional impairment and are most prominent in the amnestic subtype.

Whole-Brain Multimodal Neuroimaging Model Using Serotonin Receptor Maps Explains Non-linear Functional Effects of LSD

Understanding the underlying mechanisms of the human brain in health and disease will require models with necessary and sufficient details to explain how function emerges from the underlying anatomy and is shaped by neuromodulation. Here, we provide such a detailed causal explanation using a whole-brain model integrating multimodal imaging in healthy human participants undergoing manipulation of the serotonin system. Specifically, we combined anatomical data from diffusion magnetic resonance imaging (dMRI) and functional magnetic resonance imaging (fMRI) with neurotransmitter data obtained with positron emission tomography (PET) of the detailed serotonin 2A receptor (5-HT_2AR) density map. This allowed us to model the resting state (with and without concurrent music listening) and mechanistically explain the functional effects of 5-HT_2AR stimulation with lysergic acid diethylamide (LSD) on healthy participants. The whole-brain model used a dynamical mean-field quantitative description of populations of excitatory and inhibitory neurons as well as the associated synaptic dynamics, where the neuronal gain function of the model is modulated by the 5-HT_2AR density. The model identified the causative mechanisms for the non-linear interactions between the neuronal and neurotransmitter system, which are uniquely linked to (1) the underlying anatomical connectivity, (2) the modulation by the specific brainwide distribution of neurotransmitter receptor density, and (3) the non-linear interactions between the two. Taking neuromodulatory activity into account when modeling global brain dynamics will lead to novel insights into human brain function in health and disease and opens exciting possibilities for drug discovery and design in neuropsychiatric disorders.

Migraine is associated with high brain 5-HT levels as indexed by 5-HT4 receptor binding

INTRODUCTION

Serotonin (5-HT) plays a role in migraine pathophysiology, but whether brain 5-HT is involved in the conversion from episodic to chronic migraine is unknown. Here, we investigated brain 5-HT levels, as indexed by 5-HT₄ receptor binding, in chronic migraine patients and evaluated whether these were associated with migraine frequency.

METHODS

Sixteen chronic migraine patients underwent a dynamic PET scan after injection of [¹¹C]SB207145, a specific 5-HT₄ receptor radioligand. Data from 15 episodic migraine patients and 16 controls were included for comparison. Quantification of 5-HT₄ receptor binding was used as a proxy for brain 5-HT levels, since 5-HT₄ receptor binding is inversely related to brain 5-HT levels.

RESULTS

Chronic migraine patients had 9.1% (95% CI: [-17%; -1.0%]) lower 5-HT₄ receptor binding compared to controls ( p = 0.039). There was no difference in 5-HT₄ receptor binding between chronic and episodic migraine patients ( p = 0.48) and no association between number of monthly migraine days and 5-HT₄ receptor binding (slope estimate 0.003, 95% CI: [-0.004; 0.715], p = 0.39).

CONCLUSION

The finding of low 5-HT₄ receptor binding suggests that cerebral levels of 5-HT are elevated in chronic migraine patients. This is in line with observations made in patients with episodic migraine. Elevated brain 5-HT levels may thus be an inherent trait of the migraine brain rather than a risk factor for conversion from episodic to chronic migraine.

The importance of small polar radiometabolites in molecular neuroimaging: A PET study with [11C]Cimbi-36 labeled in two positions

[11C]Cimbi-36, a 5-HT2A receptor agonist PET radioligand, contains three methoxy groups amenable to [11C]-labeling. In pigs, [11C]Cimbi-36 yields a polar (M1) and a less polar (M2) radiometabolite fraction, while changing the labeling to [11C]Cimbi-36_5 yields only the M1 fraction. We investigate whether changing the labeling position of [11C]Cimbi-36 eliminates M2 in humans, and if this changes the signal-to-background ratio. Six healthy volunteers each underwent two dynamic PET scans; after injection of [11C]Cimbi-36, both the M1 and M2 fraction appeared in plasma, whereas only the M1 appeared after [11C]Cimbi-36_5 injection. [11C]Cimbi-36_5 generated higher uptake than [11C]Cimbi-36 in both neocortex and cerebellum. With the simplified reference tissue model mean neocortical non-displaceable binding potential for [11C]Cimbi-36 was 1.38 ± 0.07, whereas for [11C]Cimbi-36_5, it was 1.18 ± 0.14. This significant difference can be explained by higher non-displaceable binding caused by demethylation products in the M1 fraction such as [11C]formaldehyde and/or [11C]carbon dioxide/bicarbonate. Although often considered without any impact on binding measures, we show that small polar radiometabolites can substantially decrease the signal-to-background ratio of PET radioligands for neuroimaging. Further, we find that [11C]Cimbi-36 has a better signal-to-background ratio than [11C]Cimbi-36_5, and thus will be more sensitive to changes in 5-HT2A receptor levels in the brain.

High brain serotonin levels in migraine between attacks: A 5-HT4 receptor binding PET study

Migraine has been hypothesized to be a syndrome of chronic low serotonin (5-HT) levels, but investigations of brain 5-HT levels have given equivocal results. Here, we used positron emission tomography (PET) imaging of the 5-HT4 receptor as a proxy for brain 5-HT levels. Given that the 5-HT4 receptor is inversely related to brain 5-HT levels, we hypothesized that between attacks migraine patients would have higher 5-HT4 receptor binding compared to controls. Eighteen migraine patients without aura (migraine free >48 h), and 16 age- and sex-matched controls underwent PET scans after injection of [11C]SB207145, a specific 5-HT4 receptor radioligand. An investigator blinded to group calculated a neocortical mean [11C]SB207145 binding potential (BPND). Three migraine patients reported a migraine attack within 48 h after the scan and were excluded from the primary analysis. Comparing 15 migraine patients and 16 controls, we found that migraine patients have significantly lower neocortical 5-HT4 receptor binding than controls (0.60 ± 0.09 vs. 0.67 ± 0.05, p = .024), corrected for 5-HTTLPR genotype, sex and age. We found no association between 5-HT4 receptor binding and attack frequency, years with migraine or time since last migraine attack. Our finding of lower 5-HT4 receptor binding in migraine patients is suggestive of higher brain 5-HT levels. This is in contrast with the current belief that migraine is associated with low brain 5-HT levels. High brain 5-HT levels may represent a trait of the migraine brain or it could be a consequence of migraine attacks.

No evidence for direct effects of recombinant human erythropoietin on cerebral blood flow and metabolism in healthy humans

Erythropoietin (EPO) is expressed in human brain tissue, but its exact role is unknown. EPO may improve the efficiency of oxidative metabolism and has neuroprotective properties against hypoxic injuries in animal models. We aimed to investigate the effect of recombinant human EPO (rHuEPO) administration on healthy cerebral metabolism in humans during normoxia and during metabolic stress by inhalation of 10% O₂ hypoxic air. Twenty-four healthy men participated in a two-arm double-blind placebo-controlled trial. rHuEPO was administered as a low dose (5,000 IU) over 4 wk ( n = 12) or as a high dose (500 IU·kg body wt^-1·day^-1) for three consecutive days ( n = 12). Global cerebral blood flow (CBF) and metabolic rate of glucose (CMR_glc) were measured with positron emission tomography. CBF, metabolic rate of oxygen ([Formula: see text]), and cerebral lactate concentration were measured by magnetic resonance imaging and spectroscopy. Low-dose treatment increased hemoglobin and was associated with a near-significant decrease in CBF during baseline normoxia. High-dose treatment caused no change in CBF. Neither treatment had an effect on normoxia CMR_glc, [Formula: see text], or lactate concentration or an effect on the cerebral metabolic response to inhalation of hypoxic air. In conclusion, the study found no evidence for a direct effect of rHuEPO on cerebral metabolism. NEW & NOTEWORTHY We demonstrate with magnetic resonance imaging and positron emission tomography that administration of erythropoietin does not have a substantial direct effect on healthy human resting cerebral blood flow or effect on cerebral glucose and oxygen metabolism. Also, administration of erythropoietin did not have a direct effect on the metabolic response to acute hypoxic stress in healthy humans, and a suggested neuroprotective effect from erythropoietin is therefore likely not a direct effect of erythropoietin on cerebral metabolism.