Quantitative accuracy of MAP reconstruction for dynamic PET imaging in small animals

Optimized Methodology for Reference Region and Image-Derived Input Function Kinetic Modeling in Preclinical PET

PET imaging of small animals is often used for assessing biodistribution of a novel radioligand and pharmacology in small animal models of disease. PET acquisition and processing settings may affect reference region or image-derived input function (IDIF) kinetic modeling estimates. We examined four different factors in comparing quantitative results: 1) effect of reconstruction algorithm, 2) number of MAP iterations, 3) strength of the MAP prior, and 4) Attenuation and scatter. The effect of these parameters has not been explored for small-animal reference region and IDIF kinetic modeling approaches. Dynamic PET/CT scans were performed in 3 species with 3 different tracers: house sparrows with [¹¹C]raclopride, rats with [¹⁸F]AS2471907 (11βHSD1) and mice with [¹¹C]UCB-J (SV2A). FBP yielded lower kinetic modeling estimates compared to 3D-OSEM-MAP reconstructions, in sparrow and rat studies. Target resolutions (MAP prior strength) of 1.5 and 3.0mm demonstrated reduced V_T in rats but only 3.0mm reduced BP _ND in sparrows. Therefore, use of the highest target resolution (0.8mm) is warranted. We demonstrated using kinetic modeling that forgoing CT-based attenuation and scatter correction may be appropriate to improve animal throughput when using short-lived radioisotopes in sparrows and mice. This work provides recommendations and a framework for future optimization of kinetic modeling for preclinical PET methodology with novel radioligands.

Harmonization of PET image reconstruction parameters in simultaneous PET/MRI

Objective

Simultaneous PET/MRIs vary in their quantitative PET performance due to inherent differences in the physical systems and differences in the image reconstruction implementation. This variability in quantitative accuracy confounds the ability to meaningfully combine and compare data across scanners. In this work, we define image reconstruction parameters that lead to comparable contrast recovery curves across simultaneous PET/MRI systems.

Method

The NEMA NU-2 image quality phantom was imaged on one GE Signa and on one Siemens mMR PET/MRI scanner. The phantom was imaged at 9.7:1 contrast with standard spheres (diameter 10, 13, 17, 22, 28, 37 mm) and with custom spheres (diameter: 8.5, 11.5, 15, 25, 32.5, 44 mm) using a standardized methodology. Analysis was performed on a 30 min listmode data acquisition and on 6 realizations of 5 min from the listmode data. Images were reconstructed with the manufacturer provided iterative image reconstruction algorithms with and without point spread function (PSF) modeling. For both scanners, a post-reconstruction Gaussian filter of 3–7 mm in steps of 1 mm was applied. Attenuation correction was provided from a scaled computed tomography (CT) image of the phantom registered to the MR-based attenuation images and verified to align on the non-attenuation corrected PET images. For each of these image reconstruction parameter sets, contrast recovery coefficients (CRCs) were determined for the SUV_mean, SUV_max and SUV_peak for each sphere. A hybrid metric combining the root-mean-squared discrepancy (RMSD) and the absolute CRC values was used to simultaneously optimize for best match in CRC between the two scanners while simultaneously weighting toward higher resolution reconstructions. The image reconstruction parameter set was identified as the best candidate reconstruction for each vendor for harmonized PET image reconstruction.

Results

The range of clinically relevant image reconstruction parameters demonstrated widely different quantitative performance across cameras. The best match of CRC curves was obtained at the lowest RMSD values with: for CRC_mean, 2 iterations-7 mm filter on the GE Signa and 4 iterations-6 mm filter on the Siemens mMR, for CRC_max, 4 iterations-6 mm filter on the GE Signa, 4 iterations-5 mm filter on the Siemens mMR and for CRC_peak, 4 iterations-7 mm filter with PSF on the GE Signa and 4 iterations-7 mm filter on the Siemens mMR. Over all reconstructions, the RMSD between CRCs was 1.8%, 3.6% and 2.9% for CRC mean, max and peak, respectively. The solution of 2 iterations-3 mm on the GE Signa and 4 iterations-3 mm on Siemens mMR, both with PSF, led to simultaneous harmonization and with high CRC and low RMSD for CRC mean, max and peak with RMSD values of 2.8%, 5.8% and 3.2%, respectively.

Conclusions

For two commercially available PET/MRI scanners, user-selectable parameters that control iterative updates, image smoothing and PSF modeling provide a range of contrast recovery curves that allow harmonization in harmonization strategies of optimal match in CRC or high CRC values. This work demonstrates that nearly identical CRC curves can be obtained on different commercially available scanners by selecting appropriate image reconstruction parameters.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40658-021-00416-0.

Influence of Co-57 and CT Transmission Measurements on the Quantification Accuracy and Partial Volume Effect of a Small Animal PET Scanner

PURPOSE

Non-invasive in vivo positron emission tomography (PET) provides high detection sensitivity in the nano- to picomolar range and in addition to other advantages, the possibility to absolutely quantify the acquired data. The present study focuses on the comparison of transmission data acquired with an X-ray computed tomography (CT) scanner or a Co-57 source for the Inveon small animal PET scanner (Siemens Healthcare, Knoxville, TN, USA), as well as determines their influences on the quantification accuracy and partial volume effect (PVE). A special focus included the impact of the performed calibration on the quantification accuracy.

PROCEDURES

Phantom measurements were carried out to determine the quantification accuracy, the influence of the object size on the quantification, and the PVE for different sphere sizes, along the field of view and for different contrast ratios.

RESULTS

An influence of the emission activity on the Co-57 transmission measurements was discovered (deviations up to 24.06 % measured to true activity), whereas no influence of the emission activity on the CT attenuation correction was identified (deviations <3 % for measured to true activity). The quantification accuracy was substantially influenced by the applied calibration factor and by the object size. The PVE demonstrated a dependency on the sphere size, the position within the field of view, the reconstruction and correction algorithms and the count statistics. Depending on the reconstruction algorithm, only ∼30-40 % of the true activity within a small sphere could be resolved. The iterative 3D reconstruction algorithms uncovered substantially increased recovery values compared to the analytical and 2D iterative reconstruction algorithms (up to 70.46 % and 80.82 % recovery for the smallest and largest sphere using iterative 3D reconstruction algorithms). The transmission measurement (CT or Co-57 source) to correct for attenuation did not severely influence the PVE.

CONCLUSIONS

The analysis of the quantification accuracy and the PVE revealed an influence of the object size, the reconstruction algorithm and the applied corrections. Particularly, the influence of the emission activity during the transmission measurement performed with a Co-57 source must be considered. To receive comparable results, also among different scanner configurations, standardization of the acquisition (imaging parameters, as well as applied reconstruction and correction protocols) is necessary.

Oxygen-15 labeled CO2, O2, and CO PET in small animals: evaluation using a 3D-mode microPET scanner and impact of reconstruction algorithms

Background

Positron emission tomography (PET) studies using ¹⁵O-labeled CO₂, O₂, and CO have been used in humans to evaluate cerebral blood flow (CBF), the cerebral oxygen extraction fraction (OEF), and the cerebral metabolic rate of oxygen (CMRO₂) and cerebral blood volume (CBV), respectively. In preclinical studies, however, PET studies using ¹⁵O-labeled gases are not widely performed because of the technical difficulties associated with handling labeled gases with a short half-life. The aims of the present study were to evaluate the scatter fraction using 3D-mode micro-PET for ¹⁵O-labeled gas studies and the influence of reconstruction algorithms on quantitative values.

Nine male SD rats were studied using the steady state inhalation method for ¹⁵O-labeled gases with arterial blood sampling. The resulting PET images were reconstructed using filtered back projection (FBP), ordered-subset expectation maximization (OSEM) 2D, or OSEM 3D followed by maximum a posteriori (OSEM3D-MAP). The quantitative values for each brain region and each reconstruction method were calculated by applying different reconstruction methods.

Results

The quantitative values for the whole brain as calculated using FBP were 46.6 ± 12.5 mL/100 mL/min (CBF), 63.7 ± 7.2% (OEF), 5.72 ± 0.34 mL/100 mL/min (CMRO₂), and 5.66 ± 0.34 mL/100 mL (CBV), respectively. The CBF and CMRO₂ values were significantly higher when the OSEM2D and OSEM3D-MAP reconstruction methods were used, compared with FBP, whereas the OEF values were significantly lower when reconstructed using OSEM3D-MAP.

Conclusions

We evaluated the difference in quantitative values among the reconstruction algorithms using 3D-mode micro-PET. The iterative reconstruction method resulted in significantly higher quantitative values for CBF and CMRO₂, compared with the values calculated using the FBP reconstruction method.

Electronic supplementary material

The online version of this article (10.1186/s13550-017-0335-7) contains supplementary material, which is available to authorized users.

Image quality of Zr-89 PET imaging in the Siemens microPET Focus 220 preclinical scanner

PURPOSE

Zr-89 positron emission tomography (PET) is a valuable tool for understanding the biodistribution and pharmacokinetics of antibody-based therapeutics. We compared the image quality of Zr-89 PET and F-18 PET in the Siemens microPET Focus 220 preclinical scanner using different reconstruction methods.

PROCEDURES

Image quality metrics were measured in various Zr-89 and F-18 PET phantoms, including the NEMA NU 4-2008 image quality phantom. Images were reconstructed using various algorithms.

RESULTS

Zr-89 PET had greater image noise, inferior spatial resolution, and greater spillover than F-18 PET, but comparable recovery coefficients for cylinders of various diameters. Of the reconstruction methods, OSEM3D resulted in the lowest noise, highest recovery coefficients, best spatial resolution, but also the greatest spillover. Scatter correction results were found to be sensitive to varying object sizes.

CONCLUSIONS

Zr-89 PET image quality was inferior to that of F-18, and no single reconstruction method was superior in all aspects of image quality.

Iterative Image Processing for Early Diagnostic of Beta-Amyloid Plaque Deposition in Pre-Clinical Alzheimer's Disease Studies

PURPOSE

To test and evaluate an efficient iterative image processing strategy to improve the quality of sub-optimal pre-clinical PET images. A novel iterative resolution subsets-based method to reduce noise and enhance resolution (RSEMD) has been demonstrated on examples of PET imaging studies of Alzheimer's disease (AD) plaques deposition in mice brains.

MATERIALS AND METHODS

The RSEMD method was applied to imaging studies of non-invasive detection of beta-amyloid plaque in transgenic mouse models of AD. Data acquisition utilized a Siemens Inveon® micro PET/CT device. Quantitative uptake of the tracer in control and AD mice brains was determined by counting the extent of plaque deposition by histological staining. The pre-clinical imaging software inviCRO® was used for fitting the recovery PET images to the mouse brain atlas and obtaining the time activity curves (TAC) from different brain areas.

RESULTS

In all of the AD studies the post-processed images proved to have higher resolution and lower noise as compared with images reconstructed by conventional OSEM method. In general, the values of SNR reached a plateau at around 10 iterations with an improvement factor of about 2 over sub-optimal PET brain images.

CONCLUSIONS

A rapidly converging, iterative deconvolution image processing algorithm with a resolution subsets-based approach RSEMD has been used for quantitative studies of changes in Alzheimer's pathology over time. The RSEMD method can be applied to sub-optimal clinical PET brain images to improve image quality to diagnostically acceptable levels and will be crucial in order to facilitate diagnosis of AD progression at the earliest stages.

Detective quantum efficiency (DQE) in PET scanners: A simulation study

The aim of the present study is to introduce the detective quantum efficiency (DQE) for the image quality assessment of positron emission tomography (PET) scanners. For this purpose, a thin layer chromatography (TLC) plane source was simulated using a previously validated, scanner and source geometry, Monte Carlo (MC) model. The model was developed with the Geant4 application for tomographic emission (GATE) MC package and reconstructed images obtained with the software for tomographic image reconstruction (STIR), with cluster computing. The GE Discovery ST PET scanner was simulated by using a previously validated code. A plane source consisting of a TLC plate, was simulated by a layer of silica gel on aluminum (Al) foil substrate, immersed in 18F-FDG bath solution (1MBq). Image quality was assessed in terms of the modulation transfer function (MTF) and the normalized noise power spectrum (NNPS) in order to obtain the detective quantum efficiency (DQE). MTF curves were estimated from transverse reconstructed images of the plane source, whereas the NNPS data were estimated from the corresponding coronal images. Images were reconstructed by the maximum likelihood estimation ordered subsets maximum a posteriori one step late (MLE)-OS-MAP-OSL algorithm, by using various subsets 1-21) and iterations 1-20). MTF values were found to increase up to the 12th iteration whereas remain almost constant thereafter. However, the range of the increase in the MTF is limited as the number of subsets increases. The noise levels were found to increase with the corresponding increase of both the number of iterations and subsets. The maximum NNPS value (0.517mm²) was observed for the 420 MLEM-equivalent iterations reconstructed image at 0cycles/mm. Finally DQE values were found to increase for spatial frequencies up to 0.038cycles/mm and to decrease thereafter with the corresponding increase in both number of iterations and subsets. The maximum DQE value (0.48 at 0.038cycles/mm) was obtained for the 8 MLEM-equivalent iterations image. The simulated PET evaluation method based on the TLC plane source can be useful in the quality control and in the further development of PET and SPECT scanners though GATE simulations.

In-vivo quantitative assessment of the therapeutic response in a mouse model of collagen-induced arthritis using 18 F-fluorodeoxyglucose positron emission tomography

Mouse collagen-induced arthritis (CIA) is the most commonly used animal model to investigate underlying pathogenesis of autoimmune arthritis and to demonstrate the therapeutic efficacy of novel drugs in autoimmune arthritis. The conventional read-outs of CIA are clinical score and histopathology, which have several limitations, including (i) subjected to observer bias; and (ii) longitudinal therapeutic efficacy of a new drug cannot be determined. Thus, a robust, non-invasive, in-vivo drug screening tool is currently an unmet need. Here we have assessed the utility of ¹⁸ F-fluorodeoxyglucose positron emission tomography (¹⁸ F-FDG) as an in-vivo screening tool for anti-inflammatory drugs using the mouse CIA model. The radiotracer ¹⁸ F-FDG and a PET scanner were employed to monitor CIA disease activity before and after murine anti-tumour necrosis factor (TNF)-α antibody (CNTO5048) therapy in the mouse CIA model. Radiotracer concentration was derived from PET images for individual limb joints and on a per-limb basis, and Spearman's correlation coefficient (ρ) was determined with clinical score and histology of the affected limbs. CNTO5048 improved arthritis efficiently, as evidenced by clinical score and histopathology. PET showed an increased uptake of ¹⁸ F-FDG with the progression of the disease and a significant decrease in the post-treatment group. ¹⁸ F-FDG uptake patterns showed a strong correlation with clinical score (ρ = 0·71, P < 0·05) and histopathology (ρ = 0·76, P < 0·05). This study demonstrates the potential of ¹⁸ F-FDG PET as a tool for in-vivo drug screening for inflammatory arthritis and to monitor the therapeutic effects in a longitudinal setting.

Low Dose PET Image Reconstruction with Total Variation Using Alternating Direction Method

In this paper, a total variation (TV) minimization strategy is proposed to overcome the problem of sparse spatial resolution and large amounts of noise in low dose positron emission tomography (PET) imaging reconstruction. Two types of objective function were established based on two statistical models of measured PET data, least-square (LS) TV for the Gaussian distribution and Poisson-TV for the Poisson distribution. To efficiently obtain high quality reconstructed images, the alternating direction method (ADM) is used to solve these objective functions. As compared with the iterative shrinkage/thresholding (IST) based algorithms, the proposed ADM can make full use of the TV constraint and its convergence rate is faster. The performance of the proposed approach is validated through comparisons with the expectation-maximization (EM) method using synthetic and experimental biological data. In the comparisons, the results of both LS-TV and Poisson-TV are taken into consideration to find which models are more suitable for PET imaging, in particular low-dose PET. To evaluate the results quantitatively, we computed bias, variance, and the contrast recovery coefficient (CRC) and drew profiles of the reconstructed images produced by the different methods. The results show that both Poisson-TV and LS-TV can provide a high visual quality at a low dose level. The bias and variance of the proposed LS-TV and Poisson-TV methods are 20% to 74% less at all counting levels than those of the EM method. Poisson-TV gives the best performance in terms of high-accuracy reconstruction with the lowest bias and variance as compared to the ground truth (14.3% less bias and 21.9% less variance). In contrast, LS-TV gives the best performance in terms of the high contrast of the reconstruction with the highest CRC.

OSSI-PET: Open-Access Database of Simulated [(11)C]Raclopride Scans for the Inveon Preclinical PET Scanner: Application to the Optimization of Reconstruction Methods for Dynamic Studies

A wide range of medical imaging applications benefits from the availability of realistic ground truth data. In the case of positron emission tomography (PET), ground truth data is crucial to validate processing algorithms and assessing their performances. The design of such ground truth data often relies on Monte-Carlo simulation techniques. Since the creation of a large dataset is not trivial both in terms of computing time and realism, we propose the OSSI-PET database containing 350 simulated [(11)C]Raclopride dynamic scans for rats, created specifically for the Inveon pre-clinical PET scanner. The originality of this database lies on the availability of several groups of scans with controlled biological variations in the striata. Besides, each group consists of a large number of realizations (i.e., noise replicates). We present the construction methodology of this database using rat pharmacokinetic and anatomical models. A first application using the OSSI-PET database is presented. Several commonly used reconstruction techniques were compared in terms of image quality, accuracy and variability of the activity estimates and of the computed kinetic parameters. The results showed that OP-OSEM3D iterative reconstruction method outperformed the other tested methods. Analytical methods such as FBP2D and 3DRP also produced satisfactory results. However, FORE followed by OSEM2D reconstructions should be avoided. Beyond the illustration of the potential of the database, this application will help scientists to understand the different sources of noise and bias that can occur at the different steps in the processing and will be very useful for choosing appropriate reconstruction methods and parameters.

Preclinical Evaluation and Quantification of 18F-FPEB as a Radioligand for PET Imaging of the Metabotropic Glutamate Receptor 5

The metabotropic glutamate receptor 5 (mGluR5) is a high-interest target for PET imaging because it plays a role in several pathologies, including addiction, schizophrenia, and fragile X syndrome.

METHODS

We studied the pharmacokinetics of (18)F-FPEB (3-(18)F-fluoro-5-(2-pyridinylethynyl)benzonitrile), a selective PET radioligand for mGluR5, and used it to quantify mGluR5 in rat brain. Quantification was performed using both arterial sampling in combination with compartment models and simplified reference methods. The simplified reference tissue model (SRTM), Ichise's original multi-linear reference tissue model (MRTMO), and Logan noninvasive were tested as reference models with nondisplaceable binding (BPND) as outcome parameter. Additionally, test-retest scans were obtained in 6 animals.

RESULTS

(18)F-FPEB uptake in rat brain was consistent with its known distribution. No radiometabolites were present in the brain, and binding was specific as shown in blocking experiments, which also confirmed the cerebellum as a viable reference region. A 2-tissue-compartment model was used to determine BPND for the striatum (11.7 ± 1.5), nucleus accumbens (10.6 ± 2.0), hippocampus (9.0 ± 1.2), cortex (7.2 ± 1.0), and thalamus (4.0 ± 0.9). Reference methods were able to estimate these values with small bias (<2%). Test-retest analysis showed high repeatability between scans below 6%, also for shorter scan durations of 30 and 60 min.

CONCLUSION

Because of its favorable reversible kinetics, high specificity, and absence of brain radiometabolites (18)F-FPEB proves a highly useful tracer for in vivo visualization of the mGluR5 in rat brain. Moreover, reference tissue models allow noninvasive, rapid scanning with good test-retest.

Test-retest repeatability of myocardial blood flow and infarct size using ¹¹C-acetate micro-PET imaging in mice

PURPOSE

Global and regional responses of absolute myocardial blood flow index (iMBF) are used as surrogate markers to assess response to therapies in coronary artery disease. In this study, we assessed the test-retest repeatability of iMBF imaging, and the accuracy of infarct sizing in mice using (11)C-acetate PET.

METHODS

(11)C-Acetate cardiac PET images were acquired in healthy controls, endothelial nitric oxide synthase (eNOS) knockout transgenic mice, and mice after myocardial infarction (MI) to estimate global and regional iMBF, and myocardial infarct size compared to (18)F-FDG PET and ex-vivo histology results.

RESULTS

Global test-retest iMBF values had good coefficients of repeatability (CR) in healthy mice, eNOS knockout mice and normally perfused regions in MI mice (CR = 1.6, 2.0 and 1.5 mL/min/g, respectively). Infarct size measured on (11)C-acetate iMBF images was also repeatable (CR = 17 %) and showed a good correlation with the infarct sizes found on (18)F-FDG PET and histopathology (r (2) > 0.77; p < 0.05).

CONCLUSION

(11)C-Acetate micro-PET assessment of iMBF and infarct size is repeatable and suitable for serial investigation of coronary artery disease progression and therapy.

In vivo quantification of mouse autoimmune arthritis by PET/CT

AIM

To quantify the progression and severity of mouse collagen-induced arthritis (CIA) using an in vivo imaging tool, (18) F-fluorodeoxyglucose ((18) F-FDG) PET/CT and validate it against gold standard 'histopathological' evaluation.

METHOD

The PET radiotracer (18) F-FDG, a marker for glucose metabolism, was injected in mice at different stages of CIA and the radiotracer distribution was imaged using a PET scanner. A sequential CT scan provided correlated anatomy. Radiotracer concentration was derived from PET/CT images for individual limb joints and on a per-limb basis at different stages of the disease. The imaging outcomes were subjected to correlation analysis with concurrently measured clinical and histological score.

RESULTS

Clinical and histological score, and hence disease severity, showed a strong linear correlation (r(2)  = 0.71, P = 0.001 and r(2)  = 0.87, P < 0.001, respectively) with radiotracer concentration measured from PET/CT during the progression of CIA.

CONCLUSIONS

The strong positive correlation of the (18) F-FDG PET/CT findings with the histopathological evaluation at different stages of the disease suggest the potential of this imaging tool for the non-invasive assessment of progression and severity in mouse autoimmune arthritis. Thus, in preclinical studies, (18) F-FDG PET/CT can be considered as a non-invasive tool to develop novel therapies of inflammatory arthritis.

Performance assessment of a preclinical PET scanner with pinhole collimation by comparison to a coincidence-based small-animal PET scanner

PET imaging of rodents is increasingly used in preclinical research, but its utility is limited by spatial resolution and signal-to-noise ratio of the images. A recently developed preclinical PET system uses a clustered-pinhole collimator, enabling high-resolution, simultaneous imaging of PET and SPECT tracers. Pinhole collimation strongly departs from traditional electronic collimation achieved via coincidence detection in PET. We investigated the potential of such a design by direct comparison to a traditional PET scanner.

METHODS

Two small-animal PET scanners, 1 with electronic collimation and 1 with physical collimation using clustered pinholes, were used to acquire data from Jaszczak (hot rod) and uniform phantoms. Mouse brain imaging using (18)F-FDG PET was performed on each system and compared with quantitative ex vivo autoradiography as a gold standard. Bone imaging using (18)F-NaF allowed comparison of imaging in the mouse body. Images were visually and quantitatively compared using measures of contrast and noise.

RESULTS

Pinhole PET resolved the smallest rods (diameter, 0.85 mm) in the Jaszczak phantom, whereas the coincidence system resolved 1.1-mm-diameter rods. Contrast-to-noise ratios were better for pinhole PET when imaging small rods (<1.1 mm) for a wide range of activity levels, but this reversed for larger rods. Image uniformity on the coincidence system (<3%) was superior to that on the pinhole system (5%). The high (18)F-FDG uptake in the striatum of the mouse brain was fully resolved using the pinhole system, with contrast to nearby regions equaling that from autoradiography; a lower contrast was found using the coincidence PET system. For short-duration images (low-count), the coincidence system was superior.

CONCLUSION

In the cases for which small regions need to be resolved in scans with reasonably high activity or reasonably long scan times, a first-generation clustered-pinhole system can provide image quality in terms of resolution, contrast, and the contrast-to-noise ratio superior to a traditional PET system.

In Vivo Labeling of Serum Albumin for PET

The purpose of this study was to develop a novel in vivo albumin-labeling method to allow PET of cardiac function after myocardial infarction and vascular leakage and increased permeability in inflammatory diseases and malignant tumors.

METHODS

To label albumin in vivo, we synthesized a NOTA (1,4,7-triazacyclononane-N,N',N″-triacetic acid)-conjugated truncated form of Evans blue (NEB). (18)F labeling was achieved by the formation of an (18)F-aluminum fluoride ((18)F-AlF) complex, and (64)Cu labeling was obtained by a standard chelation method. Sixty-minute dynamic PET imaging was performed on normal mice to evaluate the distribution of (18)F-AlF-NEB, which was compared with in vitro-labeled mouse serum albumin ((18)F-fluorobenzyl-MSA). Electrocardiography-gated PET imaging was performed in a mouse model of myocardial infarction. Both dynamic and static PET scans were obtained in a mouse inflammation model induced by local injection of turpentine to evaluate vascular leakage. Tumor permeability was studied by dynamic and late-point static PET using (64)Cu-NEB in a UM-22B xenograft model.

RESULTS

NEB was successfully synthesized, and (18)F labeling including work-up took about 20-30 min, with a radiochemical purity greater than 95% without the need for high-performance liquid chromatography purification. Most of the radioactivity was retained in the circulation system at 60 min after injection (26.35 ± 1.52 percentage injected dose per gram [%ID/g]). With electrocardiography-gated PET, ventricles of the heart and major arteries were clearly visualized. The myocardial infarction mice showed much lower left ventricular ejection fraction than the control mice. Inflammatory muscles showed significantly higher tracer accumulation than the contralateral healthy ones. UM-22B tumor uptake of (64)Cu-NEB gradually increased with time (5.73 ± 1.11 %ID/g at 1 h and 8.03 ± 0.77 %ID/g at 2 h after injection).

CONCLUSION

The distribution and local accumulation of serum albumin can be noninvasively visualized and quantified by (18)F-AlF-NEB and (64)Cu-NEB PET. The simple labeling and broad applications make these imaging probes attractive for clinical translation.

[68Ga]-albumin-PET in the monitoring of left ventricular function in murine models of ischemic and dilated cardiomyopathy: comparison with cardiac MRI

PURPOSE

The purpose of this study is to evaluate left ventricular functional parameters in healthy mice and in different murine models of cardiomyopathy with the novel blood pool (BP) positron emission tomography (PET) tracer [68Ga]-albumin.

PROCEDURES

ECG-gated microPET examinations were obtained in healthy mice, and mice with dilative (DCM) and ischemic cardiomyopathy (ICM) using the novel BP tracer [68Ga]-albumin (AlbBP), as well as [18F]-FDG microPET. Cine-magnetic resonance imaging (MRI) examination performed on a clinical 1.5-T MRI provided the reference standard measurements.

RESULTS

When considering the combined group of healthy controls, DCM and ICM AlbBP-PET significantly overestimated the magnitudes of EDV (AlbBP, 181±86 μl; cine-MRI, 125±80 μl; P<0.001) and ESV (AlbBP, 136±92 μl; cine-MRI, 96±77 μl; P<0.001), whereas the EF (AlbBP, 31±16%; cine-MRI, 33±21%; P=0.910) matched closely to cine-MRI results, as did findings with [18F]-FDG. High correlations were found between the measured cardiac parameters (EDV: R=0.978, ESV: R=0.989, and LVEF: R=0.992).

CONCLUSIONS

Measuring left ventricular function in mice with [68Ga]-albumin BP PET is feasible and showed a high correlation compared to cine-MRI, which was used as a reference standard.

Optimization of a Model Corrected Blood Input Function from Dynamic FDG-PET Images of Small Animal Heart In Vivo

Quantitative evaluation of dynamic Positron Emission Tomography (PET) of mouse heart in vivo is challenging due to the small size of the heart and limited intrinsic spatial resolution of the PET scanner. Here, we optimized a compartment model which can simultaneously correct for spill over and partial volume effects for both blood pool and the myocardium, compute kinetic rate parameters and generate model corrected blood input function (MCBIF) from ordered subset expectation maximization - maximum a posteriori (OSEM-MAP) cardiac and respiratory gated ¹⁸F-FDG PET images of mouse heart with attenuation correction in vivo, without any invasive blood sampling. Arterial blood samples were collected from a single mouse to indicate the feasibility of the proposed method. In order to establish statistical significance, venous blood samples from n=6 mice were obtained at 2 late time points, when SP contamination from the tissue to the blood is maximum. We observed that correct bounds and initial guesses for the PV and SP coefficients accurately model the wash-in and wash-out dynamics of the tracer from mouse blood. The residual plot indicated an average difference of about 1.7% between the blood samples and MCBIF. The downstream rate of myocardial FDG influx constant, Ki (0.15±0.03 min^-1), compared well with Ki obtained from arterial blood samples (P=0.716). In conclusion, the proposed methodology is not only quantitative but also reproducible.

Impact of partial volume effect correction on cerebral β-amyloid imaging in APP-Swe mice using [(18)F]-florbetaben PET

We previously investigated the progression of β-amyloid deposition in brain of mice over-expressing amyloid-precursor protein (APP-Swe), a model of Alzheimer's disease (AD), in a longitudinal PET study with the novel β-amyloid tracer [(18)F]-florbetaben. There were certain discrepancies between PET and autoradiographic findings, which seemed to arise from partial volume effects (PVE). Since this phenomenon can lead to bias, most especially in the quantitation of brain microPET studies of mice, we aimed in the present study to investigate the magnitude of PVE on [(18)F]-florbetaben quantitation in murine brain, and to establish and validate a useful correction method (PVEC). Phantom studies with solutions of known radioactivity concentration were performed to measure the full-width-at-half-maximum (FWHM) resolution of the Siemens Inveon DPET and to validate a volume-of-interest (VOI)-based PVEC algorithm. Several VOI-brain-masks were applied to perform in vivo PVEC on [(18)F]-florbetaben data from C57BL/6(N=6) mice, while uncorrected and PVE-corrected data were cross-validated with gamma counting and autoradiography. Next, PVEC was performed on longitudinal PET data set consisting of 43 PET scans in APP-Swe (13-20months) and age-matched wild-type (WT) mice using the previously defined masks. VOI-based cortex-to-cerebellum ratios (SUVR) were compared for uncorrected and PVE-corrected results. Brains from a subset of transgenic mice were ultimately examined by autoradiography ex vivo and histochemistry in vitro as gold standard assessments, and compared to VOI-based PET results. The phantom study indicated a FWHM of 1.72mm. Applying a VOI-brain-mask including extracerebral regions gave robust PVEC, with increased precision of the SUVR results. Cortical SUVR increased with age in APP-Swe mice compared to baseline measurements (16months: +5.5%, p<0.005; 20months: +15.5%, p<0.05) with uncorrected data, and to a substantially greater extent with PVEC (16months: +12.2% p<0.005; 20months: +36.4% p<0.05). WT animals showed no binding changes, irrespective of PVEC. Relative to autoradiographic results, the error [%] for uncorrected cortical SUVR was 18.9% for native PET data, and declined to 4.8% upon PVEC, in high correlation with histochemistry results. We calculate that PVEC increases by 10% statistical power for detecting altered [(18)F]-florbetaben uptake in aging APP-Swe mice in planned studies of disease modifying treatments on amyloidogenesis.

Accurate and efficient modeling of the detector response in small animal multi-head PET systems

In fully three-dimensional PET imaging, iterative image reconstruction techniques usually outperform analytical algorithms in terms of image quality provided that an appropriate system model is used. In this study we concentrate on the calculation of an accurate system model for the YAP-(S)PET II small animal scanner, with the aim to obtain fully resolution- and contrast-recovered images at low levels of image roughness. For this purpose we calculate the system model by decomposing it into a product of five matrices: (1) a detector response component obtained via Monte Carlo simulations, (2) a geometric component which describes the scanner geometry and which is calculated via a multi-ray method, (3) a detector normalization component derived from the acquisition of a planar source, (4) a photon attenuation component calculated from x-ray computed tomography data, and finally, (5) a positron range component is formally included. This system model factorization allows the optimization of each component in terms of computation time, storage requirements and accuracy. The main contribution of this work is a new, efficient way to calculate the detector response component for rotating, planar detectors, that consists of a GEANT4 based simulation of a subset of lines of flight (LOFs) for a single detector head whereas the missing LOFs are obtained by using intrinsic detector symmetries. Additionally, we introduce and analyze a probability threshold for matrix elements of the detector component to optimize the trade-off between the matrix size in terms of non-zero elements and the resulting quality of the reconstructed images. In order to evaluate our proposed system model we reconstructed various images of objects, acquired according to the NEMA NU 4-2008 standard, and we compared them to the images reconstructed with two other system models: a model that does not include any detector response component and a model that approximates analytically the depth of interaction as detector response component. The comparisons confirm previous research results, showing that the usage of an accurate system model with a realistic detector response leads to reconstructed images with better resolution and contrast recovery at low levels of image roughness.

MicroPET-CT evaluation of interstitial brachytherapy in pancreatic carcinoma xenografts

BACKGROUND

Interstitial brachytherapy (IBT) has been introduced as treatment for unresectable pancreatic cancers to maximize local dose and minimize irradiation of surrounding normal tissue. MicroPET-CT systems provide excellent anatomic and molecular information.

PURPOSE

To use 18F-FDG micro-positron emission tomography (PET)/computed tomography (CT) to evaluate the therapeutic effect of (125)I interstitial brachytherapy on transplantation tumor of human pancreatic carcinoma in Balb/c-nu mice.

MATERIAL AND METHODS

Xenograft models were created by subcutaneous injection of Swl990 human pancreatic cancer cell suspensions into the immunodeficient Balb/c-nu mice. The study was randomly divided into three groups: control group (n = 6), empty seed implant group (n = 6), and (125)I seed implant group (n = 6), respectively. Before and 1 week after treatment, 18F-FDG microPET-CT scan was performed. In-vivo cell proliferation and apoptosis were monitored by thymidine kinase 1 (TK1) immunostaining and Dutp-biotin nick end labeling (TUNEL) assay.

RESULTS

Our results showed that before treatment the SUVmax and SUVmean values among three groups had no significant statistical difference. One week after treatment the SUVmax and SUVmean in (125)I seed implant group were significantly lower than before, while for the empty seed group and control group there were no significant difference compared with before treatment. Immunohistochemical analysis of tumor tissue revealed significantly less TK1 positive cells in (125)I seed implant group than in empty seed group and control group. The index of apoptosis was slightly higher in (125)I seed implant group than in empty seed group and control group as evaluated by TUNEL assay.

CONCLUSION

These results suggest that 18F-FDG microPET-CT may be useful as a non-invasive imaging modality to assess early response to (125)I seed brachytherapy in a pancreatic carcinoma Xenograft.

Repeatable noninvasive measurement of mouse myocardial glucose uptake with 18F-FDG: evaluation of tracer kinetics in a type 1 diabetes model

A noninvasive and repeatable method for assessing mouse myocardial glucose uptake with (18)F-FDG PET and Patlak kinetic analysis was systematically assessed using the vena cava image-derived blood input function (IDIF).

METHODS

Contrast CT and computer modeling was used to determine the vena cava recovery coefficient. Vena cava IDIF (n = 7) was compared with the left ventricular cavity IDIF, with blood and liver activity measured ex vivo at 60 min. The test-retest repeatability (n = 9) of Patlak influx constant K(i) at 10-40 min was assessed quantitatively using Bland-Altman analysis. Myocardial glucose uptake rates (rMGU) using the vena cava IDIF were calculated at baseline (n = 8), after induction of type 1 diabetes (streptozotocin [50 mg/kg] intraperitoneally, 5 d), and after acute insulin stimulation (0.08 mU/kg of body weight intraperitoneally). These changes were analyzed with a standardized uptake value calculation at 20 and 40 min after injection to correlate to the Patlak time interval.

RESULTS

The proximal mouse vena cava diameter was 2.54 ± 0.30 mm. The estimated recovery coefficient, calculated using nonlinear image reconstruction, decreased from 0.76 initially (time 0 to peak activity) to 0.61 for the duration of the scan. There was a 17% difference in the image-derived vena cava blood activity at 60 min, compared with the ex vivo blood activity measured in the γ-counter. The coefficient of variability for Patlak K(i) values between mice was found to be 23% with the proposed method, compared with 51% when using the left ventricular cavity IDIF (P < 0.05). No significant bias in K(i) was found between repeated scans with a coefficient of repeatability of 0.16 mL/min/g. Calculated rMGU values were reduced by 60% in type 1 diabetic mice from baseline scans (P < 0.03, ANOVA), with a subsequent increase of 40% to a level not significantly different from baseline after acute insulin treatment. These results were confirmed with a standardized uptake value measured at 20 and 40 min.

CONCLUSION

The mouse vena cava IDIF provides repeatable assessment of the blood time-activity curve for Patlak kinetic modeling of rMGU. An expected significant reduction in myocardial glucose uptake was demonstrated in a type 1 diabetic mouse model, with significant recovery after acute insulin treatment, using a mouse vena cava IDIF approach.

Noninvasive image derived heart input function for CMRglc measurements in small animal slow infusion FDG PET studies

Absolute quantitation of the cerebral metabolic rate for glucose (CMRglc) can be obtained in positron emission tomography (PET) studies when serial measurements of the arterial [(18)F]-fluoro-deoxyglucose (FDG) input are available. Since this is not always practical in PET studies of rodents, there has been considerable interest in defining an image-derived input function (IDIF) by placing a volume of interest (VOI) within the left ventricle of the heart. However, spill-in arising from trapping of FDG in the myocardium often leads to progressive contamination of the IDIF, which propagates to underestimation of the magnitude of CMRglc. We therefore developed a novel, non-invasive method for correcting the IDIF without scaling to a blood sample. To this end, we first obtained serial arterial samples and dynamic FDG-PET data of the head and heart in a group of eight anaesthetized rats. We fitted a bi-exponential function to the serial measurements of the IDIF, and then used the linear graphical Gjedde-Patlak method to describe the accumulation in myocardium. We next estimated the magnitude of myocardial spill-in reaching the left ventricle VOI by assuming a Gaussian point-spread function, and corrected the measured IDIF for this estimated spill-in. Finally, we calculated parametric maps of CMRglc using the corrected IDIF, and for the sake of comparison, relative to serial blood sampling from the femoral artery. The uncorrected IDIF resulted in 20% underestimation of the magnitude of CMRglc relative to the gold standard arterial input method. However, there was no bias with the corrected IDIF, which was robust to the variable extent of myocardial tracer uptake, such that there was a very high correlation between individual CMRglc measurements using the corrected IDIF with gold-standard arterial input results. Based on simulation, we furthermore find that electrocardiogram-gating, i.e. ECG-gating is not necessary for IDIF quantitation using our approach.