Performance assessment of the ALBIRA II pre-clinical SPECT S102 system for 99mTc imaging.

OBJECTIVE: The performance characteristics of the SPECT sub-system S102 of the ALBIRA II PET/SPECT/CT are analyzed for the 80 mm field of view (FOV) to evaluate the potential in-vivo imaging in rats, based on measurements of the system response for the commonly used Technetium-99 m (99mTc) in small animal imaging. METHODS: The ALBIRA II tri-modal µPET/SPECT/CT pre-clinical system (Bruker BioSpin, Ettlingen, Germany) was used. The SPECT modality is made up of two opposite gamma cameras (Version S102) with Sodium doped Cesium Iodide (CsI(Na)) single continuous crystal detectors coupled to position-sensitive photomultipliers (PSPMTs). Imaging was performed with the NEMA NU-4 image quality phantom (Data Spectrum Corporation, Durham, USA). Measurements were performed with a starting activity concentration of 4.76 MBq/mL 99mTc. An energy window of 20% at 140 keV was selected in this study. The system offers a 20 mm, 40 mm, 60 mm and an 80 mm field of view (FOV) and in this study the 80 mm FOV was used for all the acquisitions. The data were reconstructed with an ordered subset expectation maximization (OSEM) algorithm. Sensitivity, spatial resolution, count rate linearity, convergence of the algorithm and the recovery coefficients (RC) were analyzed. All analyses were performed with PMOD and MATLAB software. RESULTS: The sensitivities measured at the center of the 80 mm FOV with the point source were 23.1 ± 0.3 cps/MBq (single pinhole SPH) and 105.6 ± 5.5 cps/MBq (multi pinhole MPH). The values for the axial, tangential and radial full width at half maximum (FWHM) were 2.51, 2.54, and 2.55 mm with SPH and 2.35, 2.44 and 2.32 mm with MPH, respectively. The corresponding RC values for the 5 mm, 4 mm, 3 mm and 2 mm rods were 0.60 ± 0.28, 0.61 ± 0.24, 0.29 ± 0.11 and 0.20 ± 0.06 with SPH and 0.56 ± 0.20, 0.50 ± 0.18, 0.38 ± 0.09 and 0.23 ± 0.06 with MPH. To obtain quantitative imaging data, the image reconstructions should be performed with 12 iterations. CONCLUSION: The ALBIRA II preclinical SPECT sub-system S102 has a favorable sensitivity and spatial resolution for the 80 mm FOV setting for both the SPH and MPH configurations and is a valuable tool for small animal imaging.

Low-Dose Imaging in a New Preclinical Total-Body PET/CT Scanner.

Ionizing radiation constitutes a health risk to imaging scientists and study animals. Both PET and CT produce ionizing radiation. CT doses in pre-clinical in vivo imaging typically range from 50 to 1,000 mGy and biological effects in mice at this dose range have been previously described. [18F]FDG body doses in mice have been estimated to be in the range of 100 mGy for [18F]FDG. Yearly, the average whole body doses due to handling of activity by PET technologists are reported to be 3-8 mSv. A preclinical PET/CT system is presented with design features which make it suitable for small animal low-dose imaging. The CT subsystem uses a X-source power that is optimized for small animal imaging. The system design incorporates a spatial beam shaper coupled with a highly sensitive flat-panel detector and very fast acquisition (<10 s) which allows for whole body scans with doses as low as 3 mGy. The mouse total-body PET subsystem uses a detector architecture based on continuous crystals, coupled to SiPM arrays and a readout based in rows and columns. The PET field of view is 150 mm axial and 80 mm transaxial. The high solid-angle coverage of the sample and the use of continuous crystals achieve a sensitivity of 9% (NEMA) that can be leveraged for use of low tracer doses and/or performing rapid scans. The low-dose imaging capabilities of the total-body PET subsystem were tested with NEMA phantoms, in tumor models, a mouse bone metabolism scan and a rat heart dynamic scan. The CT imaging capabilities were tested in mice and in a low contrast phantom. The PET low-dose phantom and animal experiments provide evidence that image quality suitable for preclinical PET studies is achieved. Furthermore, CT image contrast using low dose scan settings was suitable as a reference for PET scans. Total-body mouse PET/CT studies could be completed with total doses of <10 mGy.

Prediction of time-integrated activity coefficients in PRRT using simulated dynamic PET and a pharmacokinetic model.

PURPOSE: To investigate the accuracy of predicted time-integrated activity coefficients (TIACs) in peptide-receptor radionuclide therapy (PRRT) using simulated dynamic PET data and a physiologically based pharmacokinetic (PBPK) model. METHODS: PBPK parameters were estimated using biokinetic data of 15 patients after injection of (152±15)MBq of 111In-DTPAOC (total peptide amount (5.78±0.25)nmol). True mathematical phantoms of patients (MPPs) were the PBPK model with the estimated parameters. Dynamic PET measurements were simulated as being done after bolus injection of 150MBq 68Ga-DOTATATE using the true MPPs. Dynamic PET scans around 35min p.i. (P1), 4h p.i. (P2) and the combination of P1 and P2 (P3) were simulated. Each measurement was simulated with four frames of 5min each and 2 bed positions. PBPK parameters were fitted to the PET data to derive the PET-predicted MPPs. Therapy was simulated assuming an infusion of 5.1GBq of 90Y-DOTATATE over 30min in both true and PET-predicted MPPs. TIACs of simulated therapy were calculated, true MPPs (true TIACs) and predicted MPPs (predicted TIACs) followed by the calculation of variabilities v. RESULTS: For P1 and P2 the population variabilities of kidneys, liver and spleen were acceptable (v<10%). For the tumours and the remainders, the values were large (up to 25%). For P3, population variabilities for all organs including the remainder further improved, except that of the tumour (v>10%). CONCLUSION: Treatment planning of PRRT based on dynamic PET data seems possible for the kidneys, liver and spleen using a PBPK model and patient specific information.

Investigation of the imaging characteristics of the ALBIRA II small animal PET system for 18F, 68Ga and 64Cu.

AIM: In this study the performance characteristics of the Albira II PET sub-system and the response of the system for the following radionuclides 18F, 68Ga and 64Cu was analyzed. MATERIALS AND METHODS: The Albira II tri-modal system (Bruker BioSpin MRI GmbH, Ettlingen, Germany) is a pre-clinical device for PET, SPECT and CT. The PET sub-system uses single continuous crystal detectors of lutetium yttrium orthosilicate (LYSO). The detector assembly consists of three rings of 8 detector modules. The transaxial field of view (FOV) has a diameter of 80mm and the axial FOV is 148mm. A NEMA NU-4 image quality phantom (Data Spectrum Corporation, Durham, USA) having five rods with diameters of 1, 2, 3, 4 and 5mm and a uniform central region was used. Measurements with 18F, 68Ga and 64Cu were performed in list mode acquisition over 10h. Data were reconstructed using a maximum-likelihood expectation-maximization (MLEM) algorithm with iteration numbers between 5 and 50. System sensitivity, count rate linearity, convergence and recovery coefficients were analyzed. RESULTS: The sensitivities for the entire FOV (non-NEMA method) for 18F, 68Ga and 64Cu were (3.78±0.05)%, (3.97±0.18)% and (3.79±0.37)%, respectively. The sensitivity based on the NEMA protocol using the 22Na point source yielded (5.53±0.06)%. Dead-time corrected true counts were linear for activities ≤7MBq (18F and 68Ga) and ≤17MBq (64Cu) in the phantom. The radial, tangential and axial full widths at half maximum (FWHMs) were 1.52, 1.47 and 1.48mm. Recovery coefficients for the uniform region with a total activity of 8MBq in the phantom were (0.97±0.05), (0.98±0.06), (0.98±0.06) for 18F, 68Ga and 64Cu, respectively. CONCLUSION: The Albira II pre-clinical PET system has an adequate sensitivity range and the system linearity is suitable for the range of activities used for pre-clinical imaging. Overall, the system showed a favorable image quality for pre-clinical applications.

Comparison of five cluster validity indices performance in brain [18 F]FET-PET image segmentation using k-means.

PURPOSE: Dynamic [18 F]fluoro-ethyl-L-tyrosine positron emission tomography ([18 F]FET-PET) is used to identify tumor lesions for radiotherapy treatment planning, to differentiate glioma recurrence from radiation necrosis and to classify gliomas grading. To segment different regions in the brain k-means cluster analysis can be used. The main disadvantage of k-means is that the number of clusters must be pre-defined. In this study, we therefore compared different cluster validity indices for automated and reproducible determination of the optimal number of clusters based on the dynamic PET data. METHODS: The k-means algorithm was applied to dynamic [18 F]FET-PET images of 8 patients. Akaike information criterion (AIC), WB, I, modified Dunn's and Silhouette indices were compared on their ability to determine the optimal number of clusters based on requirements for an adequate cluster validity index. To check the reproducibility of k-means, the coefficients of variation CVs of the objective function values OFVs (sum of squared Euclidean distances within each cluster) were calculated using 100 random centroid initialization replications RCI100 for 2 to 50 clusters. k-means was performed independently on three neighboring slices containing tumor for each patient to investigate the stability of the optimal number of clusters within them. To check the independence of the validity indices on the number of voxels, cluster analysis was applied after duplication of a slice selected from each patient. CVs of index values were calculated at the optimal number of clusters using RCI100 to investigate the reproducibility of the validity indices. To check if the indices have a single extremum, visual inspection was performed on the replication with minimum OFV from RCI100 . RESULTS: The maximum CV of OFVs was 2.7 × 10-2 from all patients. The optimal number of clusters given by modified Dunn's and Silhouette indices was 2 or 3 leading to a very poor segmentation. WB and I indices suggested in median 5, [range 4-6] and 4, [range 3-6] clusters, respectively. For WB, I, modified Dunn's and Silhouette validity indices the suggested optimal number of clusters was not affected by the number of the voxels. The maximum coefficient of variation of WB, I, modified Dunn's, and Silhouette validity indices were 3 × 10-2 , 1, 2 × 10-1 and 3 × 10-3 , respectively. WB-index showed a single global maximum, whereas the other indices showed also local extrema. CONCLUSION: From the investigated cluster validity indices, the WB-index is best suited for automated determination of the optimal number of clusters for [18 F]FET-PET brain images for the investigated image reconstruction algorithm and the used scanner: it yields meaningful results allowing better differentiation of tissues with higher number of clusters, it is simple, reproducible and has an unique global minimum.

Treatment planning in PRRT based on simulated PET data and a PBPK model. Determination of accuracy using a PET noise model.

AIM: To investigate the accuracy of treatment planning in peptide-receptor radionuclide therapy (PRRT) based on simulated PET data (using a PET noise model) and a physiologically based pharmacokinetic (PBPK) model. METHODS: The parameters of a PBPK model were fitted to the biokinetic data of 15 patients. True mathematical phantoms of patients (MPPs) were the PBPK model with the fitted parameters. PET measurements after bolus injection of 150 MBq 68Ga-DOTATATE were simulated for the true MPPs. PET noise with typical noise levels was added to the data (i.e. c = 0.3 [low], 3, 30 and 300 [high]). Organ activity data in the kidneys, tumour, liver and spleen were simulated at 0.5, 1 and 4 h p.i. PBPK model parameters were fitted to the simulated noisy PET data to derive the PET-predicted MPPs. Therapy was simulated assuming an infusion of 3.3 GBq of 90Y-DOTATATE over 30 min. Time-integrated activity coefficients (TIACs) of simulated therapy in tumour, kidneys, liver, spleen and remainder were calculated from both, true MPPs (true TIACs) and predicted MPPs (predicted TIACs). Variability v between true TIACs and predicted TIACs were calculated and analysed. Variability ≤ 10 % was considered to be an accurate prediction. RESULTS: For all noise level, variabilities for the kidneys, liver, and spleen showed an accurate prediction for TIACs, e.g. c = 300: vkidney = (5 ± 2)%, vliver = (5 ± 2)%, vspleen = (4 ± 2)%. However, tumour TIAC predictions were not accurate for all noise levels, e.g. c = 0.3: vtumour = (8 ± 5)%. CONCLUSION: PET-based treatment planning with kidneys as the dose limiting organ seems possible for all reported noise levels using an adequate PBPK model and previous knowledge about the individual patient.

Physiologically based pharmacokinetic modeling of (18)F-SiFAlin-Asp3-PEG1-TATE in AR42J tumor bearing mice.

PURPOSE: Peptide receptor radionuclide therapy (PRRT) is commonly performed in the treatment of neuroendocrine tumors (NET), where somatostatin analogs (DOTATATE) are radiolabeled with (90)Y, (68)Ga or (111)In for pre-therapeutic and therapeutic purposes. Quantitative evaluation of the biokinetic data can be performed by using physiologically based pharmacokinetic (PBPK) models. Knowledge about the biodistribution in a pre-clinical setting would allow optimizing the translation from bench to bedside. The aim of this study was to develop a PBPK model to describe the biodistribution of a novel sst2-targeting radiotracer. METHODS: Biokinetic data of six mice after injection of (18)F-SiFAlin-Asp3-PEG1-TATE were investigated using two PBPK models. The PBPK models describe the biodistribution of the tracer in the tumor, kidneys, liver, remainder and whole body via blood flow to these organs via absorption, distribution, metabolism and excretion. A recently published sst2 PBPK model for humans (model 1) was used to describe the data. Physiological information in this model was adapted to that of a mouse. Model 1 was further modified by implementing receptor-mediated endocytosis (model 2). Model parameters were fitted to the biokinetic data of each mouse. Model selection was performed by calculating Akaike weights wi using the corrected Akaike Information Criterion (AICc). RESULTS: The implementation of receptor-mediated endocytosis considerably improved the description of the biodistribution (Akaike weights w1=0% and w2=100% for model 1 and 2, respectively). The resulting time-integrated activity coefficients determined by model 2 were for tumor (0.05 ± 0.02) h, kidneys (0.11 ± 0.01) h and liver (0.02 ± 0.01) h. CONCLUSION: Simply downscaling a human PBPK model does not allow for an accurate description of (18)F-SiFAlin-Asp3-PEG1-TATE in mice. Biokinetics of this tracer can be accurately and adequately described using a physiologically based pharmacokinetic model including receptor-mediated endocytosis. Thus, an optimized translation from bench to bedside is possible.

The role of patient-based treatment planning in peptide receptor radionuclide therapy.

BACKGROUND: Accurate treatment planning is recommended in peptide-receptor radionuclide therapy (PRRT) to minimize the toxicity to organs at risk while maximizing tumor cell sterilization. The aim of this study was to quantify the effect of different degrees of individualization on the prediction accuracy of individual therapeutic biodistributions in patients with neuroendocrine tumors (NETs). METHODS: A recently developed physiologically based pharmacokinetic (PBPK) model was fitted to the biokinetic data of 15 patients with NETs after pre-therapeutic injection of (111)In-DTPAOC. Mathematical phantom patients (MPP) were defined using the assumed true (true MPP), mean (MPP 1A) and median (MPP 1B) parameter values of the patient group. Alterations of the degree of individualization were introduced to both mean and median patients by including patient-specific information as a priori knowledge: physical parameters and hematocrit (MPP 2A/2B). Successively, measurable individual biokinetic parameters were added: tumor volume V tu (MPP 3A/3B), glomerular filtration rate GFR (MPP 4A/4B), and tumor perfusion f tu (MPP 5A/5B). Furthermore, parameters of MPP 5A/5B and a simulated (68)Ga-DOTATATE PET measurement 60 min p.i. were used together with the population values used as Bayesian parameters (MPP 6A/6B). Therapeutic biodistributions were simulated assuming an infusion of (90)Y-DOTATATE (3.3 GBq) over 30 min to all MPPs. Time-integrated activity coefficients were predicted for all MPPs and compared to the true MPPs for each patient in tumor, kidneys, spleen, liver, remainder, and whole body to obtain the relative differences RD. RESULTS: The large RD values of MPP 1A [RDtumor = (625 ± 1266)%, RDkidneys = (11 ± 38)%], and MPP 1B [RDtumor = (197 ± 505)%, RDkidneys = (11 ± 39)%] demonstrate that individual treatment planning is needed due to large physiological differences between patients. Although addition of individual patient parameters reduced the deviations considerably [MPP 5A: RDtumor = (-2 ± 27)% and RDkidneys = (16 ± 43)%; MPP 5B: RDtumor = (2 ± 28)% and RDkidneys = (7 ± 40)%] errors were still large. For the kidneys, prediction accuracy was considerably improved by including the PET measurement [MPP 6A/MPP 6B: RDtumor = (-2 ± 22)% and RDkidneys = (-0.1 ± 0.5)%]. CONCLUSION: Individualized treatment planning is needed in the investigated patient group. The use of a PBPK model and the inclusion of patient specific data, e.g., weight, tumor volume, and glomerular filtration rate, do not suffice to predict the therapeutic biodistribution. Integrating all available a priori information in the PBPK model and using additionally PET data measured at one time point for tumor, kidneys, spleen, and liver could possibly be sufficient to perform an individualized treatment planning.

Quantitative and qualitative assessment of Yttrium-90 PET/CT imaging.

Yttrium-90 is known to have a low positron emission decay of 32 ppm that may allow for personalized dosimetry of liver cancer therapy with (90)Y labeled microspheres. The aim of this work was to image and quantify (90)Y so that accurate predictions of the absorbed dose can be made. The measurements were performed within the QUEST study (University of Sydney, and Sirtex Medical, Australia). A NEMA IEC body phantom containing 6 fillable spheres (10-37 mm ∅) was used to measure the 90Y distribution with a Biograph mCT PET/CT (Siemens, Erlangen, Germany) with time-of-flight (TOF) acquisition. A sphere to background ratio of 8:1, with a total (90)Y activity of 3 GBq was used. Measurements were performed for one week (0, 3, 5 and 7 d). he acquisition protocol consisted of 30 min-2 bed positions and 120 min-single bed position. Images were reconstructed with 3D ordered subset expectation maximization (OSEM) and point spread function (PSF) for iteration numbers of 1-12 with 21 (TOF) and 24 (non-TOF) subsets and CT based attenuation and scatter correction. Convergence of algorithms and activity recovery was assessed based on regions-of-interest (ROI) analysis of the background (100 voxels), spheres (4 voxels) and the central low density insert (25 voxels). For the largest sphere, the recovery coefficient (RC) values for the 30 min -2-bed position, 30 min-single bed and 120 min-single bed were 1.12 ± 0.20, 1.14 ± 0.13, 0.97 ± 0.07 respectively. For the smaller diameter spheres, the PSF algorithm with TOF and single bed acquisition provided a comparatively better activity recovery. Quantification of Y-90 using Biograph mCT PET/CT is possible with a reasonable accuracy, the limitations being the size of the lesion and the activity concentration present. At this stage, based on our study, it seems advantageous to use different protocols depending on the size of the lesion.

Disruption of the mouse Jhy gene causes abnormal ciliary microtubule patterning and juvenile hydrocephalus.

Congenital hydrocephalus, the accumulation of excess cerebrospinal fluid (CSF) in the ventricles of the brain, affects one of every 1000 children born today, making it one of the most common human developmental disorders. Genetic causes of hydrocephalus are poorly understood in humans, but animal models suggest a broad genetic program underlying the regulation of CSF balance. In this study, the random integration of a transgene into the mouse genome led to the development of an early onset and rapidly progressive hydrocephalus. Juvenile hydrocephalus transgenic mice (Jhy(lacZ)) inherit communicating hydrocephalus in an autosomal recessive fashion with dilation of the lateral ventricles observed as early as postnatal day 1.5. Ventricular dilation increases in severity over time, becoming fatal at 4-8 weeks of age. The ependymal cilia lining the lateral ventricles are morphologically abnormal and reduced in number in Jhy(lacZ/lacZ) brains, and ultrastructural analysis revealed disorganization of the expected 9+2 microtubule pattern. Rather, the majority of Jhy(lacZ/lacZ) cilia develop axonemes with 9+0 or 8+2 microtubule structures. Disruption of an unstudied gene, 4931429I11Rik (now named Jhy) appears to underlie the hydrocephalus of Jhy(lacZ/lacZ) mice, and the Jhy transcript and protein are decreased in Jhy(lacZ/lacZ) mice. Partial phenotypic rescue was achieved in Jhy(lacZ/lacZ) mice by the introduction of a bacterial artificial chromosome (BAC) carrying 60-70% of the JHY protein coding sequence. Jhy is evolutionarily conserved from humans to basal vertebrates, but the predicted JHY protein lacks identifiable functional domains. Ongoing studies are directed at uncovering the physiological function of JHY and its role in CSF homeostasis.