Three-dimensional maximum a posteriori (MAP) imaging with radiopharmaceuticals labeled with three Cu radionuclides

Fast and cost-effective cyclotron production of61Cu using anatZn liquid target: an opportunity for radiopharmaceutical production and R&D

Following our previous work on the production of radiometals, such as⁶⁴Cu and⁶⁸Ga, through the irradiation of liquid targets using a medical cyclotron, we describe in this paper a technique to produce⁶¹Cu through the irradiation of natural zinc using a liquid target.

Assessment of a fully 3D Monte Carlo reconstruction method for preclinical PET with iodine-124

Iodine-124 is a radionuclide well suited to the labeling of intact monoclonal antibodies. Yet, accurate quantification in preclinical imaging with I-124 is challenging due to the large positron range and a complex decay scheme including high-energy gammas. The aim of this work was to assess the quantitative performance of a fully 3D Monte Carlo (MC) reconstruction for preclinical I-124 PET. The high-resolution small animal PET Inveon (Siemens) was simulated using GATE 6.1. Three system matrices (SM) of different complexity were calculated in addition to a Siddon-based ray tracing approach for comparison purpose. Each system matrix accounted for a more or less complete description of the physics processes both in the scanned object and in the PET scanner. One homogeneous water phantom and three heterogeneous phantoms including water, lungs and bones were simulated, where hot and cold regions were used to assess activity recovery as well as the trade-off between contrast recovery and noise in different regions. The benefit of accounting for scatter, attenuation, positron range and spurious coincidences occurring in the object when calculating the system matrix used to reconstruct I-124 PET images was highlighted. We found that the use of an MC SM including a thorough modelling of the detector response and physical effects in a uniform water-equivalent phantom was efficient to get reasonable quantitative accuracy in homogeneous and heterogeneous phantoms. Modelling the phantom heterogeneities in the SM did not necessarily yield the most accurate estimate of the activity distribution, due to the high variance affecting many SM elements in the most sophisticated SM.

Isotope specific resolution recovery image reconstruction in high resolution PET imaging

PURPOSE

Measuring and incorporating a scanner-specific point spread function (PSF) within image reconstruction has been shown to improve spatial resolution in PET. However, due to the short half-life of clinically used isotopes, other long-lived isotopes not used in clinical practice are used to perform the PSF measurements. As such, non-optimal PSF models that do not correspond to those needed for the data to be reconstructed are used within resolution modeling (RM) image reconstruction, usually underestimating the true PSF owing to the difference in positron range. In high resolution brain and preclinical imaging, this effect is of particular importance since the PSFs become more positron range limited and isotope-specific PSFs can help maximize the performance benefit from using resolution recovery image reconstruction algorithms.

METHODS

In this work, the authors used a printing technique to simultaneously measure multiple point sources on the High Resolution Research Tomograph (HRRT), and the authors demonstrated the feasibility of deriving isotope-dependent system matrices from fluorine-18 and carbon-11 point sources. Furthermore, the authors evaluated the impact of incorporating them within RM image reconstruction, using carbon-11 phantom and clinical datasets on the HRRT.

RESULTS

The results obtained using these two isotopes illustrate that even small differences in positron range can result in different PSF maps, leading to further improvements in contrast recovery when used in image reconstruction. The difference is more pronounced in the centre of the field-of-view where the full width at half maximum (FWHM) from the positron range has a larger contribution to the overall FWHM compared to the edge where the parallax error dominates the overall FWHM.

CONCLUSIONS

Based on the proposed methodology, measured isotope-specific and spatially variant PSFs can be reliably derived and used for improved spatial resolution and variance performance in resolution recovery image reconstruction. The benefits are expected to be more substantial for more energetic positron emitting isotopes such as Oxygen-15 and Rubidium-82.

Noise propagation in resolution modeled PET imaging and its impact on detectability

Positron emission tomography imaging is affected by a number of resolution degrading phenomena, including positron range, photon non-collinearity and inter-crystal blurring. An approach to this issue is to model some or all of these effects within the image reconstruction task, referred to as resolution modeling (RM). This approach is commonly observed to yield images of higher resolution and subsequently contrast, and can be thought of as improving the modulation transfer function. Nonetheless, RM can substantially alter the noise distribution. In this work, we utilize noise propagation models in order to accurately characterize the noise texture of reconstructed images in the presence of RM. Furthermore we consider the task of lesion or defect detection, which is highly determined by the noise distribution as quantified using the noise power spectrum. Ultimately, we use this framework to demonstrate why conventional trade-off analyses (e.g. contrast versus noise, using simplistic noise metrics) do not provide a complete picture of the impact of RM and that improved performance of RM according to such analyses does not necessarily translate to the superiority of RM in detection task performance.

Resolution modeling in PET imaging: theory, practice, benefits, and pitfalls

In this paper, the authors review the field of resolution modeling in positron emission tomography (PET) image reconstruction, also referred to as point-spread-function modeling. The review includes theoretical analysis of the resolution modeling framework as well as an overview of various approaches in the literature. It also discusses potential advantages gained via this approach, as discussed with reference to various metrics and tasks, including lesion detection observer studies. Furthermore, attention is paid to issues arising from this approach including the pervasive problem of edge artifacts, as well as explanation and potential remedies for this phenomenon. Furthermore, the authors emphasize limitations encountered in the context of quantitative PET imaging, wherein increased intervoxel correlations due to resolution modeling can lead to significant loss of precision (reproducibility) for small regions of interest, which can be a considerable pitfall depending on the task of interest.

Development of copper based drugs, radiopharmaceuticals and medical materials

Copper is one of the most interesting elements for various biomedical applications. Copper compounds show vast array of biological actions, including anti-inflammatory, anti-proliferative, biocidal and other. It also offers a selection of radioisotopes, suitable for nuclear imaging and radiotherapy. Quick progress in nanotechnology opened new possibilities for design of copper based drugs and medical materials. To date, copper has not found many uses in medicine, but number of ongoing research, as well as preclinical and clinical studies, will most likely lead to many novel applications of copper in the near future.

The motivations and methodology for high-throughput PET imaging of small animals in cancer research

Over the last decade, small-animal PET imaging has become a vital platform technology in cancer research. With the development of molecularly targeted therapies and drug combinations requiring evaluation of different schedules, the number of animals to be imaged within a PET experiment has increased. This paper describes experimental design requirements to reach statistical significance, based on the expected change in tracer uptake in treated animals as compared to the control group, the number of groups that will be imaged, and the expected intra-animal variability for a given tracer. We also review how high-throughput studies can be performed in dedicated small-animal PET, high-resolution clinical PET systems and planar positron imaging systems by imaging more than one animal simultaneously. Customized beds designed to image more than one animal in large-bore small-animal PET scanners are described. Physics issues related to the presence of several rodents within the field of view (i.e. deterioration of spatial resolution and sensitivity as the radial and the axial offsets increase, respectively, as well as a larger effect of attenuation and the number of scatter events), which can be assessed by using the NEMA NU 4 image quality phantom, are detailed.

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

PURPOSE

Iterative reconstruction algorithms are becoming more commonly employed in positron emission tomography (PET) imaging; however, the quantitative accuracy of the reconstructed images still requires validation for various levels of contrast and counting statistics.

METHODS

The authors present an evaluation of the quantitative accuracy of the 3D maximum a posteriori (3D-MAP) image reconstruction algorithm for dynamic PET imaging with comparisons to two of the most widely used reconstruction algorithms: the 2D filtered-backprojection (2D-FBP) and 2D-ordered subsets expectation maximization (2D-OSEM) on the Siemens microPET scanners. The study was performed for various levels of count density encountered in typical dynamic scanning as well as the imaging of cardiac activity concentration in small animal studies on the Focus 120. Specially designed phantoms were used for evaluation of the spatial resolution, image quality, and quantitative accuracy. A normal mouse was employed to evaluate the accuracy of the blood time activity concentration extracted from left ventricle regions of interest (ROIs) within the images as compared to the actual blood activity concentration measured from arterial blood sampling.

RESULTS

For MAP reconstructions, the spatial resolution and contrast have been found to reach a stable value after 20 iterations independent of the β values (i.e., hyper parameter which controls the weight of the penalty term) and count density within the frame. The spatial resolution obtained with 3D-MAP reaches values of ∼1.0 mm with a β of 0.01 while the 2D-FBP has value of 1.8 mm and 2D-OSEM has a value of 1.6 mm. It has been observed that the lower the hyper parameter β used in MAP, more iterations are needed to reach the stable noise level (i.e., image roughness). The spatial resolution is improved by using a lower β value at the expense of higher image noise. However, with similar noise level the spatial resolution achieved by 3D-MAP was observed to be better than that by 2D-FBP or 2D-OSEM. Using an image quality phantom containing hot spheres, the estimated activity concentration in the largest sphere has the expected concentration relative to the background area for all the MAP images. The obtained recovery coefficients have been also shown to be almost independent of the count density. 2D-FBP and 2D-OSEM do not perform as well, yielding recovery coefficients lower than those observed with 3D-MAP (approximately 33% lower for the smallest sphere). However, a small positive bias was observed in MAP reconstructed images for frames of very low count density. This bias is present in the uniform area for count density of less than 0.05 × 10(6) counts/ml. For the dynamic mouse study, it was observed that 3D-MAP (even gated at diastole) cannot predict accurately the blood activity concentration due to residual spill-over activity from the myocardium into the left ventricle (approximately 15%). However, 3D-MAP predicts blood activity concentration closer to blood sampling than 2D-FBP.

CONCLUSIONS

The authors observed that 3D-MAP produces more accurate activity concentration estimates than 2D-FBP or 2D-OSEM at all practical levels of statistics and contrasts due to improved spatial resolution leading to lesser partial volume effect.

A practical guide to the construction of radiometallated bioconjugates for positron emission tomography

Positron emission tomography (PET) has become a vital imaging modality in the diagnosis and treatment of disease, most notably cancer. A wide array of small molecule PET radiotracers have been developed that employ the short half-life radionuclides (11)C, (13)N, (15)O, and (18)F. However, PET radiopharmaceuticals based on biomolecular targeting vectors have been the subject of dramatically increased research in both the laboratory and the clinic. Typically based on antibodies, oligopeptides, or oligonucleotides, these tracers have longer biological half-lives than their small molecule counterparts and thus require labeling with radionuclides with longer, complementary radioactive half-lives, such as the metallic isotopes (64)Cu, (68)Ga, (86)Y, and (89)Zr. Each bioconjugate radiopharmaceutical has four component parts: biomolecular vector, radiometal, chelator, and covalent link between chelator and biomolecule. With the exception of the radiometal, a tremendous variety of choices exists for each of these pieces, and a plethora of different chelation, conjugation, and radiometallation strategies have been utilized to create agents ranging from (68)Ga-labeled pentapeptides to (89)Zr-labeled monoclonal antibodies. Herein, the authors present a practical guide to the construction of radiometal-based PET bioconjugates, in which the design choices and synthetic details of a wide range of biomolecular tracers from the literature are collected in a single reference. In assembling this information, the authors hope both to illuminate the diverse methods employed in the synthesis of these agents and also to create a useful reference for molecular imaging researchers both experienced and new to the field.

Molecular SPECT Imaging: An Overview

Molecular imaging has witnessed a tremendous change over the last decade. Growing interest and emphasis are placed on this specialized technology represented by developing new scanners, pharmaceutical drugs, diagnostic agents, new therapeutic regimens, and ultimately, significant improvement of patient health care. Single photon emission computed tomography (SPECT) and positron emission tomography (PET) have their signature on paving the way to molecular diagnostics and personalized medicine. The former will be the topic of the current paper where the authors address the current position of the molecular SPECT imaging among other imaging techniques, describing strengths and weaknesses, differences between SPECT and PET, and focusing on different SPECT designs and detection systems. Radiopharmaceutical compounds of clinical as well-preclinical interest have also been reviewed. Moreover, the last section covers several application, of μSPECT imaging in many areas of disease detection and diagnosis.

Gold nanocages as photothermal transducers for cancer treatment

Gold nanocages represent a new class of nanomaterials with compact size and tunable optical properties for biomedical applications. They exhibit strong light absorption in the near-infrared region in which light can penetrate deeply into soft tissue. After PEGylation, the Au nanocages can be passively delivered to tumors in animals. Analysis of tissue distribution for the PEGylated Au nanocages showed that the tumor uptake was 5.7 %ID/g at 96 h post injection. The Au nanocages were found not only on the surface, but also in the core of the tumor. By exposing tumors to a near-infrared diode laser (0.7 W/cm2, CW, λ=808 nm) for 10 min, the photothermal effect of the Au nanocages could selectively destroy tumor tissue with minimum damage to the surrounding healthy tissue. Data from functional [18F]fluorodexoyglucose positron emission tomography revealed a decrease in tumor metabolic activity upon the photothermal treatment. Histological examination identified extensive damage to the nuclei of tumor cells and tumor interstitium.

Radioimmunoimaging with longer-lived positron-emitting radionuclides: potentials and challenges

Radioimmunoimaging and therapy has been an area of interest for several decades. Steady progress has been made toward clinical translation of radiolabeled monoclonal antibodies for diagnosis and treatment of diseases. Tremendous advances have been made in imaging technologies such as positron emission tomography (PET). However, these advances have so far eluded routine translation into clinical radioimmunoimaging applications due to the mismatch between the short half-lives of routinely used positron-emitting radionuclides such as (18)F versus the pharmacokinetics of most intact monoclonal antibodies of interest. The lack of suitable positron-emitting radionuclides that match the pharmacokinetics of intact antibodies has generated interest in exploring the use of longer-lived positron emitters that are more suitable for radioimmunoimaging and dosimetry applications with intact monoclonal antibodies. In this review, we examine the opportunities and challenges of radioimmunoimaging with select longer-lived positron-emitting radionuclides such as (124)I, (89)Zr, and (86)Y with respect to radionuclide production, ease of radiolabeling intact antibodies, imaging characteristics, radiation dosimetry, and clinical translation potential.

Quantitative small animal PET imaging with nonconventional nuclides

Positron Emission Tomography (PET) has gained a tremendous momentum recently for clinical applications notably with the availability of (18)F-fluorodeoxyglucose for staging and evaluation of therapy efficacy in various types of cancers. Nonconventional positron emitting nuclides are now being investigated for the development of novel imaging and therapeutic strategies. However, these nuclides have less than ideal imaging properties. This article compares the performance for imaging of nonconventional nuclides such as (61)Cu, (68)Ga, (86)Y and (94m)Tc with the standard imaging nuclide (18)F for high-resolution small animal PET imaging. Quantitative imaging performance was evaluated in terms of spatial resolution and hot spheres recovery coefficients from image resolution and image quality phantoms representing the mouse. The data were reconstructed using algorithms of 2D filtered-back-projection, 2D ordered-subsets expectation maximization and maximum-a-posteriori. It is shown that the spatial resolution point spread function can be well explained by a double-gaussian function due to the generally long range of the positron. We show that, with the knowledge of the measured point spread functions, the accurate activity concentration in small lesions can be recovered when imaging with long-range positron emitters.

Optical observation of energy loss distribution and practical range of positrons from a 18F water solution in a water-equivalent phantom

The energy loss distribution of beta+ particles is closely related to their maximum penetration depth distribution and annihilation point distribution. The latter is of practical importance for positron emission tomography. Experimental data related to the energy loss distribution are important for comprehensive validation of physics and simulation models of beta+ interactions. In this paper the authors present a new experimental approach that allows them to visually observe the beta+ energy loss distribution of a solution of nuclear medicine radioisotopes in a plastic scintillator using an optical camera. The authors also report a set of the first experimental results. A water solution of 18F was localized in a small hole in a plastic scintillator (BC430). Optical imaging of the scintillator yielded visual images of the energy loss distribution with a submillimeter resolution. The radial dependence in the energy distribution was quantitatively measured by analysis of the images, and exponential fitting parameters were obtained. The authors observed that the results of Monte Carlo simulation with EGS5 (version 1.0.2) and GEANT4 (version 4.9.01.p01) were consistent with those obtained experimentally. The results of the Monte Carlo simulation indicated that for a linear scale, the energy loss distribution in the scintillator was approximately the same as that in water, and the relative shape of the energy loss distribution was close to those of the maximum penetration depth distribution and annihilation point distribution. This paper also presents discussions about the further possibilities of this optical imaging approach. Thus, optical observation of the beta+ energy loss distribution in a scintillator is a promising technique for visual and quantitative experimental studies of beta+ emission from a solution of radioisotopes that are used in nuclear medicine.

Analytic system matrix resolution modeling in PET: an application to Rb-82 cardiac imaging

This work explores application of a novel resolution modeling technique based on analytic physical models which individually models the various resolution degrading effects in PET (positron range, photon non-collinearity, inter-crystal scattering and inter-crystal penetration) followed by their combination and incorporation within the image reconstruction task. In addition to phantom studies, the proposed technique was particularly applied to and studied in the task of clinical Rb-82 myocardial perfusion imaging, which presently suffers from poor statistics and resolution properties in the reconstructed images. Overall, the approach is able to produce considerable enhancements in image quality. The reconstructed FWHM for a Discovery RX PET/CT scanner was seen to improve from 5.1 mm to 7.7 mm across the field-of-view (FoV) to approximately 3.5 mm nearly uniformly across the FoV. Furthermore, extended-source phantom studies indicated clearly improved images in terms of contrast versus noise performance. Using Monte Carlo simulations of clinical Rb-82 imaging, the resolution modeling technique was seen to significantly outperform standard reconstructions qualitatively, and also quantitatively in terms of contrast versus noise (contrast between the myocardium and other organs, as well as between myocardial defects and the left ventricle).

Quantitative serial imaging of an 124I anti-CEA monoclonal antibody in tumor-bearing mice

OBJECTIVE

The 4.2-day half-life of (124)I favors its use for positron emission tomography (PET) of monoclonal antibodies (mAbs). However, high positron energy and beta(+) -associated cascade gamma rays pose image resolution and background noise problems for (124)I. This study evaluated quantitative PET of an (124)I mAb in tumor-bearing mice.

METHODS

An R4 microPETtrade mark (Siemens/CTIMI, Knoxville, TN) was used with standard energy and coincidence timing windows (350-750 keV and 6 ns, respectively), delayed random coincidence subtraction, iterative image reconstruction, and no attenuation or scatter correction. Image resolution, contrast, and response linearity were compared for (124)I and (18)F, using phantoms. Nude mice bearing human colon tumors (LS-174T) were injected intravenously with a chimeric (124)I anti-CEA mAb (cT84.66) and imaged serially 1 hour to 7 days postinjection. Venous blood was sampled to validate image-derived blood curves. Mice were sacrificed after the final scan, and the biodistribution of (124)I was measured by direct tissue assay. Images were converted to units of kBq/g for each tissue of interest by comparing the final scans with the direct assays.

RESULTS

Measured resolution (FWHM) 0-16 mm from the scanner axis was 2.3-2.7 mm for (124)I versus 1.9-2.0 mm for (18)F. Due to true coincidence events between annihilation photons and cascade gamma rays, background was greater for (124)I than (18)F, but the signal-to-background ratio was still more than 20, and (124)I image intensities varied linearly with activity concentration. Tissue-based calibration worked well (i.e., PET blood curves agreed with direct measurements within 12% at all time points), while calibration, based on a cylindrical phantom approximating the mouse body, yielded tumor quantitation that was 46%-66% low, compared with direct assay.

CONCLUSIONS

Images of quantitative accuracy sufficient for biodistribution measurements can be obtained from tumor-bearing mice by using (124)I anti-CEA mAbs with standard microPET acquisition and processing techniques, provided the calibration is based on the direct assay of excised tissue samples.

Annihilation photon acollinearity in PET: volunteer and phantom FDG studies

Annihilation photon acollinearity is a fundamental but little investigated problem in positron emission tomography (PET). In this paper, the cause of the angular deviation from 180.00 degrees is described as well as how to evaluate it under conditions of a spatially distributed radiation source and a limited acquisition time for the human body. A relationship between the shape of the photopeak spectrum and the angular distribution is formulated using conservation laws of momentum and energy over the pair annihilation. Then the formula is used to evaluate the acollinearity for a pool phantom and the human body with FDG injected. The angular distribution for the pool phantom agrees well with that for pure water which had been directly measured by Colombino et al in 1965 (Nuovo Cimento 38 707-23), and also with that for the human body determined in this study. Pure water can be considered as a good approximation of the human body regarding the angular deviation. The blurring coefficient to be multiplied by the ring diameter in calculations of the PET spatial resolution is experimentally determined for the first time as 0.00243 +/- 0.00014; this is 10% larger than the value widely used by investigators.

Iterative image reconstruction using inverse Fourier rebinning for fully 3-D PET

We describe a fast forward and back projector pair based on inverse Fourier rebinning for use in iterative image reconstruction for fully 3-D positron emission tomography (PET). The projector pair is used as part of a factored system matrix that takes into account detector-pair response by using shift-variant sinogram blur kernels, thereby combining the computational advantages of Fourier rebinning with iterative reconstruction using accurate system models. The forward projector consists of a 2-D projector, which maps 3-D images into 2-D direct sinograms, followed by exact inverse rebinning which maps the 2-D into fully 3-D sinograms. The back projector is implemented as the transpose of the forward projector and differs from the true exact rebinning operator in the sense that it does not require reprojection to compute missing lines of response (LORs). We compensate for two types of inaccuracies that arise in a cylindrical PET scanner when using inverse Fourier rebinning: 1) nonuniform radial sampling and 2) nonconstant oblique angles in the radial direction in a single oblique sinogram. We examine the effects of these corrections on sinogram accuracy and reconstructed image quality. We evaluate performance of the new projector pair for maximum a posteriori (MAP) reconstruction of simulated and in vivo data. The new projector results in only a small loss in resolution towards the edge of the field-of-view when compared to the fully 3-D geometric projector and requires an order of magnitude less computation.

Advances in the production, processing and microPET image quality of technetium-94m

This work involves the production, processing and imaging of the short-lived, rarely used positron emission tomography (PET) radionuclide technetium-94m (94mTc). Our procedures are an extension of methods reported in the literature and are detailed within. A key modification was the development of a single step that combines purification and concentration of an aqueous 94mTc-pertechnetate solution, which both reduces processing time and increases the final concentration of the solution. Additionally, a convenient method for the direct recovery of 94mTc into an organic solvent was developed, eliminating the solvent transfer step needed for organic syntheses using 94mTc. Each of these advances potentially extends the scope of syntheses possible with this short-lived radionuclide. To explore the imaging potential of 94mTc, we carried out phantom imaging studies on small-scale high-resolution PET scanners to estimate the limitations of detection associated with 94mTc and PET. Preliminary studies demonstrate that useful images can be obtained with modern image reconstruction algorithms when using a correction for the cascade gamma ray contamination.

Production of non-standard PET radionuclides and the application of radiopharmaceuticals labeled with these nuclides

The field of positron emission tomography (PET) has expanded dramatically over recent years. In spite of this expansion the large majority of clinical studies are carried out utilizing one radiopharmaceutical-2-fluoro-2-deoxyglucose. Many research groups are developing novel radiopharmaceuticals. A major emphasis is on other agents labeled with 18F. Several other positron emitting radionuclides can be prepared in high yields in small biomedical cyclotrons. Some of these have half-lives that make delivery significantly easier than the delivery of 18F compounds. These radionuclides include: 64Cu (half life 12.7 h), 76Br (half life 16.2 h), 86Y (half life 14.74 h) and 124I (half life 4.2 days). The method of production of these and other 'non-standard' PET radionuclides will be discussed and the method of labeling radiopharmaceuticals with these radionuclides described. Several of these radiopharmaceuticals have been studied in animal models as well and a limited number translated to the human situation.