Michel Koole

GATE: a simulation toolkit for PET and SPECT.

Monte Carlo simulation is an essential tool in emission tomography that can assist in the design of new medical imaging devices, the optimization of acquisition protocols and the development or assessment of image reconstruction algorithms and correction techniques. GATE, the Geant4 Application for Tomographic Emission, encapsulates the Geant4 libraries to achieve a modular, versatile, scripted simulation toolkit adapted to the field of nuclear medicine. In particular, GATE allows the description of time-dependent phenomena such as source or detector movement, and source decay kinetics. This feature makes it possible to simulate time curves under realistic acquisition conditions and to test dynamic reconstruction algorithms. This paper gives a detailed description of the design and development of GATE by the OpenGATE collaboration, whose continuing objective is to improve, document and validate GATE by simulating commercially available imaging systems for PET and SPECT. Large effort is also invested in the ability and the flexibility to model novel detection systems or systems still under design. A public release of GATE licensed under the GNU Lesser General Public License can be downloaded at http:/www-lphe.epfl.ch/GATE/. Two benchmarks developed for PET and SPECT to test the installation of GATE and to serve as a tutorial for the users are presented. Extensive validation of the GATE simulation platform has been started, comparing simulations and measurements on commercially available acquisition systems. References to those results are listed. The future prospects towards the gridification of GATE and its extension to other domains such as dosimetry are also discussed.

Phase 1 Study of the Pittsburgh Compound B Derivative 18F-Flutemetamol in Healthy Volunteers and Patients with Probable Alzheimer Disease

11C-Pittsburgh compound B (PiB) marks Aβ amyloidosis, a key pathogenetic process in Alzheimer disease (AD). The use of 11C-PiB is limited to centers with a cyclotron. Development of the 18F-labeled thioflavin derivative of PiB, 18F-flutemetamol, could hugely increase the availability of this new technology. The aims of this phase 1 study were to perform brain kinetic modeling of 18F-flutemetamol, optimize the image acquisition procedure, and compare methods of analysis (step 1) and to compare 18F-flutemetamol brain retention in AD patients versus healthy controls in a proof-of-concept study (steps 1 and 2). Methods: In step 1, 3 AD patients (Mini-Mental State Examination, 22–24) and 3 elderly healthy controls were scanned dynamically during windows of 0–90, 150–180, and 220–250 min after injection of approximately 180 MBq of 18F-flutemetamol, with arterial sampling. We compared different analysis methods (compartmental modeling, Logan graphical analysis, and standardized uptake value ratios) and determined the optimal acquisition window for step 2. In step 2, 5 AD patients (Mini-Mental State Examination, 20–26) and 5 elderly healthy controls were scanned from 80 to 170 min after injection. To determine overall efficacy, steps 1 and 2 were pooled and standardized uptake value ratios were calculated using cerebellar cortex as a reference region. Results: No adverse events were reported. There was a strong correlation between uptake values obtained with the different analysis methods. From 80 min after injection onward, the ratio of neocortical to cerebellar uptake was maximal and only marginally affected by scan start time or duration. AD patients showed significantly increased standardized uptake value ratios in neocortical association zones and striatum, compared with healthy controls, whereas uptake in white matter, cerebellum, and pons did not differ between groups. Two AD patients were 18F-flutemetamol–negative and 1 healthy control was 18F-flutemetamol–positive. Conclusion: 18F-flutemetamol uptake can be readily quantified. This phase 1 study warrants further studies to validate this 18F-labeled derivative of PiB as a biomarker for Aβ amyloidosis.

Whole-Body Biodistribution and Radiation Dosimetry of 18F-GE067: A Radioligand for In Vivo Brain Amyloid Imaging

We have characterized the biodistribution and dosimetry of 18F-3′-F-6-OH-BTA1 (18F-GE067), a newly developed radioligand to visualize and quantify amyloid burden, in healthy elderly human subjects. Methods: Six subjects (5 men and 1 woman; age range, 51–74 y) underwent dynamic whole-body PET/CT for 6 h after a bolus injection of 18F-GE067. Source organs were delineated on PET/CT. Individual organ doses and effective doses were determined. Results: No adverse events or clinically significant changes were observed. 18F-GE067 is excreted predominantly through the hepatobiliary system. The gallbladder, upper large intestine, and small intestine are the organs with the highest absorbed dose (average, 287, 173, and 155 μGy/MBq, respectively). The mean effective dose was 33.8 ± 3.4 μSv/MBq, a dose comparable to that of many other 18F-labeled radiopharmaceuticals. Conclusion: The estimated effective dose of 18F-GE067 for PET amyloid imaging was acceptable (class II-b defined by the World Health Organization), and relatively low variability between subjects was observed.

99mTc-ECD brain perfusion SPET: variability, asymmetry and effects of age and gender in healthy adults.

Abstract. Reliable and high-resolution reference data for regional cerebral blood flow measured with single-photon emission tomography (SPET) are necessary for optimal clinical and research use. Therefore, a large dataset of normal technetium-99m labelled ethylene cysteine dimer (ECD) perfusion SPET in carefully screened healthy volunteers with an age range spanning six decades was created, with correction for non-uniform attenuation and scatter and based on an anatomically standardised analysis. Eighty-nine healthy volunteers, stratified for gender (46 females, 43 males; age 20–81 years), were included. Twelve volunteers underwent repeated 99mTc-ECD SPET after 2.5±2.3 weeks. An automated whole-brain volume of interest analysis with MANOVA as well as voxelwise analysis using SPM99 was conducted. Average intersubject variability was 4.8% while intrasubject reproducibility was 3.0%. An age-related decline in tracer uptake was found in the anterior cingulate gyrus, bilateral basal ganglia, left prefrontal, left lateral frontal and left superior temporal and insular cortex (all P=0.001–0.02). There was an overall increase in right/left asymmetry with age, which was most pronounced in the frontal and temporal neocortex. The most significant correlations between AI and age decade were found in the prefrontal (R=0.35, P=0.001) and superior temporal neocortex (R=0.43, P<0.001). Women had significantly higher uptake in the right parietal cortex (P<0.001), while men showed higher uptake in the cerebellum and the left anterior temporal and orbitofrontal cortex (all P<0.01). This normative dataset allows age- and gender-specific patient and group assessment of 99mTc-ECD perfusion SPET under a wide variety of clinical circumstances in relation to normal variations and highlights the importance of both age- and gender-specific normal datasets for optimal analysis sensitivity.

Monte Carlo simulations of a scintillation camera using GATE: validation and application modelling

Geant4 application for tomographic emission (GATE) is a recently developed simulation platform based on Geant4, specifically designed for PET and SPECT studies. In this paper we present validation results of GATE based on the comparison of simulations against experimental data, acquired with a standard SPECT camera. The most important components of the scintillation camera were modelled. The photoelectric effect. Compton and Rayleigh scatter are included in the gamma transport process. Special attention was paid to the processes involved in the collimator: scatter, penetration and lead fluorescence. A LEHR and a MEGP collimator were modelled as closely as possible to their shape and dimensions. In the validation study, we compared the simulated and measured energy spectra of different isotopes: 99mTc, 22Na, 57Co and 67Ga. The sensitivity was evaluated by using sources at varying distances from the detector surface. Scatter component analysis was performed in different energy windows at different distances from the detector and for different attenuation geometries. Spatial resolution was evaluated using a 99mTc source at various distances. Overall results showed very good agreement between the acquisitions and the simulations. The clinical usefulness of GATE depends on its ability to use voxelized datasets. Therefore, a clinical extension was written so that digital patient data can be read in by the simulator as a source distribution or as an attenuating geometry. Following this validation we modelled two additional camera designs: the Beacon transmission device for attenuation correction and the Solstice scanner prototype with a rotating collimator. For the first setup a scatter analysis was performed and for the latter design. the simulated sensitivity results were compared against theoretical predictions. Both case studies demonstrated the flexibility and accuracy of GATE and exemplified its potential benefits in protocol optimization and in system design.

Preclinical Evaluation of 18F-JNJ41510417 as a Radioligand for PET Imaging of Phosphodiesterase-10A in the Brain

Phosphodiesterases are enzymes that inactivate the intracellular second messengers 3′,5′-cyclic adenosine-monophosphate and/or cyclic guanosine-monophosphate. Of all 11 known phosphodiesterase families, phosphodiesterase-10A (PDE10A) has the most restricted distribution, with high expression in the striatum. PDE10A inhibitors are pursued as drugs for treatment of neuropsychiatric disorders. We have synthesized and evaluated 18F-JNJ41510417 as a selective and high-affinity radioligand for in vivo brain imaging of PDE10A using PET. Methods: The biodistribution of 18F-JNJ41510417 was evaluated in rats. Rat plasma and perfused brain homogenates were analyzed by high-performance liquid chromatography to quantify radiometabolites. Dynamic small-animal PET was performed in rats and in wild-type and PDE10A knock-out mice and compared with ex vivo autoradiography. Blocking and displacement experiments were performed using the nonradioactive analog and other selective PDE10A inhibitors. Results: Tissue distribution studies showed predominant hepatobiliary excretion, sufficient brain uptake (0.56 ± 0.00 percentage injected dose at 2 min after tracer injection), and continuous accumulation of the tracer in the striatum over time; rapid washout of nonspecific binding from other brain regions was observed. Polar radiometabolites were detected in plasma and brain tissue. Dynamic small-animal PET showed continuous tracer accumulation in the striatum, with rapid decline in the cortex and cerebellum. Pretreatment and chase experiments with PDE10A inhibitors showed that the tracer binding to PDE10A was specific and reversible. Imaging in PDE10A knock-out and wild-type mice further confirmed that binding in the striatum was specific for PDE10A. Conclusion: Experiments in rats and PDE10A knock-out mice indicate that 18F-JNJ41510417 binds specifically and reversibly to PDE10A in the striatum, suggesting that this new fluorinated quinoline derivative is a promising candidate for in vivo imaging of PDE10A using PET.

Experimental Performance Assessment of SPM for SPECT Neuroactivation Studies Using a Subresolution Sandwich Phantom Design

The validity domain of voxel-based statistical analysis of SPECT neuroactivation studies with statistical parametric mapping (SPM) has been investigated by a limited number of theoretical and simulation studies. In this work, an experimental setup is described with an assessment of its activation detection performance together with the influence of acquisition and processing parameters. A subresolution sandwich phantom was constructed using a printed high-resolution digital Hoffman phantom with a (99m)TcO(4)-ink mixture. Activations of 8, 16, and 24 mm diameter, with an intensity ranging from 5 to 60%, were constructed in the right frontal cortex, anterior and posterior cingulate, and left striatum. Triple-headed SPECT acquisitions were carried out using both fan-beam and parallel beam geometry. The impact of activation characteristics (size, intensity and location), study size, physical degradation factors, and reconstruction technique were studied using SPM99 in a group comparison design with correction for multiple comparisons. For a 15 x 15 design, all 24-mm activations of 5% intensity were detected for the fan-beam data. Decreased focus or study size, lower SPECT resolution, absence of scatter, and attenuation correction resulted in an increase in minimally detectable activation. For a single study referred to 15 control studies, only 24-mm activation foci with a minimal intensity of 10% were detected in the optimal configuration. This approach allows experimental parameter optimization of SPM-based group or single-subject SPECT activation studies compared to normal data, as used in clinical applications. In principle, these findings can be extended to SPECT receptor studies or PET data.

Whole-Body Biodistribution and Radiation Dosimetry of the Human Cannabinoid Type-1 Receptor Ligand 18F-MK-9470 in Healthy Subjects

The cannabinoid type-1 (CB1) receptor is one of the most abundant G-coupled protein receptors in the human body and is responsible for signal transduction of both endogenous and exogenous cannabinoids. The endocannabinoid system is strongly implicated in regulation of homeostasis and several neuropsychiatric disorders, obesity, and associated comorbidities, such as dyslipidemia and metabolic syndrome. We have used whole-body PET/CT to characterize the biodistribution and dosimetry of a novel high-affinity, subtype-selective radioligand, 18F-MK-9470, in healthy male and female subjects. Methods: Eight nonobese subjects (5 men, 3 women; age, 22–54 y) underwent serial whole-body PET/CT for 6 h after a bolus injection of 251 ± 25 MBq 18F-MK-9470 (N-[2-(3-cyano-phenyl)-3-(4-(2-18F-fluorethoxy)phenyl)-1-methylpropyl]-2-(5-methyl-2-pyridyloxy)-2-methylproponamide). Source organs were delineated 3-dimensionally using the combined morphologic and functional data. Residence times were derived from time–activity profiles using both the trapezoid rule and curve fitting. Individual organ doses and effective doses were determined using the OLINDA software package, with different approaches for gastrointestinal and urinary excretion modeling. Results: 18F-MK-9470 is taken up slowly in the brain, reaching a plateau at approximately 90–120 min after bolus injection and is excreted predominantly through the hepatobiliary system. The gallbladder, upper large intestine, small intestine, and liver are the organs with the highest absorbed dose (average: 159, 98, 87, and 86 μGy/MBq, respectively). The mean effective dose (ED) was 22.8 ± 4.3 μSv/MBq, indicating relatively low intersubject variability and a mean value in the range of many commercially available 18F-labeled radiopharmaceuticals. Brain uptake was relatively high compared with that of existing central nervous system ligands for other receptors, between 3.2% and 4.9% of the injected dose. Conclusion: The estimated radiation burden of 18F-MK-9470 for PET CB1 receptor imaging shows relatively low variability between subjects and has an acceptable ED, which allows multiple serial cerebral scans of good image quality, while remaining within the risk category class II-b defined by the World Health Organization and the International Commission for Radiation Protection for a standard injected activity (185–370 MBq).

Micro-Positron Emission Tomography Imaging of Rat Brain Metabolism during Expression of Contextual Conditioning

Using 18F-fluorodeoxygluose microPET imaging, we investigated the neurocircuitry of contextual anxiety versus control in awake, conditioned rats (n = 7–10 per group). In addition, we imaged a group expressing cued fear. Simultaneous measurements of startle amplitude and freezing time were used to assess conditioning. To the best of our knowledge, no neuroimaging studies in conditioned rats have been conducted thus far, although visualizing and quantifying the metabolism of the intact brain in behaving animals is clearly of interest. In addition, more insight into the neurocircuitry involved in contextual anxiety may stimulate the development of new treatments for anxiety disorders. Our main finding was hypermetabolism in a cluster comprising the bed nucleus of the stria terminalis (BST) in rats expressing contextual anxiety compared with controls. Analysis of a subset of rats showing the best behavioral results (n = 5 per subgroup) confirmed this finding. We also observed hypermetabolism in the same cluster in rats expressing contextual anxiety compared with rats expressing cued fear. Our results provide novel evidence for a role of the BST in the expression of contextual anxiety.

A three-dimensional theoretical model incorporating spatial detection uncertainty in continuous detector PET.

In this paper, we will describe a theoretical model of the spatial uncertainty for a line of response, due to the imperfect localization of events on the detector heads of a positron emission tomography (PET) camera. The forward acquisition problem is modelled by a Gaussian distribution of the position of interaction on a detector head, centred at the measured position. The a posteriori probability that an event originates from a certain point in the field of view (FOV) is calculated by integrating all the possible lines of response (LORs) through this point, weighted with the Gaussian detection likelihood at the LORs end points. We have calculated these a posteriori probabilities both for perpendicular and oblique coincidences. For the oblique coincidence case it was necessary to incorporate the effect of the crystal thickness in the calculations. We found in the perpendicular incidence case as well as in the oblique incidence case that the probability density function cannot be analytically expressed in a closed form, and it was thus calculated by means of numerical integration. A Gaussian was fit to the transversal profiles of this function for a given distance to the detectors. From these fits, we can conclude that the profiles can be accurately approximated by a Gaussian, both for perpendicular and oblique coincidences. The FWHM reaches a maximum at the detector heads, and decreases towards the centre of the FOV, as was expected. Afterwards we extended this two-dimensional model to three dimensions, thus incorporating the spatial uncertainty in both transversal directions. This theoretical model was then evaluated and a very good agreement was found with theoretical calculations and with geometric Monte Carlo simulations. Possible improvements for the above-described incorporation of crystal thickness are discussed. Therefore a detailed Monte Carlo study has been performed in order to investigate the interaction probability of photons of different energies along their path in several detector materials dedicated to PET. Finally two approaches for the incorporation of this theoretical model in reconstruction algorithms are outlined.