Ultrasound-based sensors for respiratory motion assessment in multimodality PET imaging.

Breathing motion can displace internal organs by up to several cm; as such, it is a primary factor limiting image quality in medical imaging. Motion can also complicate matters when trying to fuse images from different modalities, acquired at different locations and/or on different days. Currently available devices for monitoring breathing motion often do so indirectly, by detecting changes in the outline of the torso rather than the internal motion itself, and these devices are often fixed to floors, ceilings or walls, and thus cannot accompany patients from one location to another. We have developed small ultrasound-based sensors, referred to as 'organ configuration motion' (OCM) sensors, that attach to the skin and provide rich motion-sensitive information. In the present work we tested the ability of OCM sensors to enable respiratory gating duringin vivoPET imaging. A motion phantom involving an FDG solution was assembled, and two cancer patients scheduled for a clinical PET/CT exam were recruited for this study. OCM signals were used to help reconstruct phantom andin vivodata into time series of motion-resolved images. As expected, the motion-resolved images captured the underlying motion. In Patient #1, a single large lesion proved to be mostly stationary through the breathing cycle. However, in Patient #2, several small lesions were mobile during breathing, and our proposed new approach captured their breathing-related displacements. In summary, a relatively inexpensive hardware solution was developed here for respiration monitoring. Because the proposed sensors attach to the skin, as opposed to walls or ceilings, they can accompany patients from one procedure to the next, potentially allowing data gathered in different places and at different times to be combined and compared in ways that account for breathing motion.

Partial volume correction for improved PET quantification in 18F-NaF imaging of atherosclerotic plaques.

BACKGROUND: Accurate quantification of plaque imaging using 18F-NaF PET requires partial volume correction (PVC). METHODS: PVC of PET data was implemented by the use of a local projection (LP) method. LP-based PVC was evaluated with an image quality (NEMA) and with a thorax phantom with "plaque-type" lesions of 18-36 mL. The validated PVC method was then applied to a cohort of 17 patients, each with at least one plaque in the carotid or ascending aortic arteries. In total, 51 calcified (HU > 110) and 16 non-calcified plaque lesions (HU < 110) were analyzed. The lesion-to-background ratio (LBR) and the relative change of LBR (ΔLBR) were measured on PET. RESULTS: Following PVC, LBR of the spheres (NEMA phantom) was within 10% of the original values. LBR of the thoracic lesions increased by 155% to 440% when the LP-PVC method was applied to the PET images. In patients, PVC increased the LBR in both calcified [mean = 78% (-8% to 227%)] and non-calcified plaques [mean = 41%, (-9%-104%)]. CONCLUSIONS: PVC helps to improve LBR of plaque-type lesions in both phantom studies and clinical patients. Better results were obtained when the PVC method was applied to images reconstructed with point spread function modeling.

Dose reduction in half-time myocardial perfusion SPECT-CT with multifocal collimation.

BACKGROUND: Recent technological advances in myocardial perfusion imaging may warrant the use of lower injected activity. We evaluated whether quantitative measures of stress myocardial perfusion defects using Tc-99m sestamibi and low-energy high-resolution (LEHR) collimators are equivalent to lower dose SPECT-CT with cardiac multifocal collimators and software (IQ·SPECT). METHODS: 93 patients underwent one-day rest-stress gated SPECT-CT. Following conventional rest imaging, 925-1100 MBq (25-30 mCi) of Tc-99m sestamibi was injected during stress testing. Stress SPECT-CT images were acquired two ways: with LEHR (13 minutes) and IQ·SPECT (7 minutes). Low-dose IQ·SPECT stress was simulated by subsampling the full-dose data to half-, quarter-, and eighth-count levels. Abnormalities were quantified using the total perfusion deficit (TPD) score and dose-specific databases. RESULTS: The mean ± SD of the differences between LEHR and IQ·SPECT TPD scores were -1.01 ± 5.36%, -0.10 ± 5.81%, 1.78 ± 4.81%, and 1.75 ± 6.05% at full, half, quarter, and eighth doses, respectively. Differences were statistically significant for quarter and eighth doses. Correlation between LEHR and IQ·SPECT was excellent at all doses (R ≥ 0.93). Bland-Altman plots demonstrated minimal bias. CONCLUSIONS: With IQ·SPECT, quantitative stress SPECT-CT imaging is possible with half of the standard injected activity in half the time.

Mitochondrial iron chelation ameliorates cigarette smoke-induced bronchitis and emphysema in mice.

Chronic obstructive pulmonary disease (COPD) is linked to both cigarette smoking and genetic determinants. We have previously identified iron-responsive element-binding protein 2 (IRP2) as an important COPD susceptibility gene and have shown that IRP2 protein is increased in the lungs of individuals with COPD. Here we demonstrate that mice deficient in Irp2 were protected from cigarette smoke (CS)-induced experimental COPD. By integrating RNA immunoprecipitation followed by sequencing (RIP-seq), RNA sequencing (RNA-seq), and gene expression and functional enrichment clustering analysis, we identified Irp2 as a regulator of mitochondrial function in the lungs of mice. Irp2 increased mitochondrial iron loading and levels of cytochrome c oxidase (COX), which led to mitochondrial dysfunction and subsequent experimental COPD. Frataxin-deficient mice, which had higher mitochondrial iron loading, showed impaired airway mucociliary clearance (MCC) and higher pulmonary inflammation at baseline, whereas mice deficient in the synthesis of cytochrome c oxidase, which have reduced COX, were protected from CS-induced pulmonary inflammation and impairment of MCC. Mice treated with a mitochondrial iron chelator or mice fed a low-iron diet were protected from CS-induced COPD. Mitochondrial iron chelation also alleviated CS-induced impairment of MCC, CS-induced pulmonary inflammation and CS-associated lung injury in mice with established COPD, suggesting a critical functional role and potential therapeutic intervention for the mitochondrial-iron axis in COPD.

Use of radiopharmaceuticals in diagnostic nuclear medicine in the United States: 1960-2010.

To reconstruct reliable nuclear medicine-related occupational radiation doses or doses received as patients from radiopharmaceuticals over the last five decades, the authors assessed which radiopharmaceuticals were used in different time periods, their relative frequency of use, and typical values of the administered activity. This paper presents data on the changing patterns of clinical use of radiopharmaceuticals and documents the range of activity administered to adult patients undergoing diagnostic nuclear medicine procedures in the U.S. between 1960 and 2010. Data are presented for 15 diagnostic imaging procedures that include thyroid scan and thyroid uptake; brain scan; brain blood flow; lung perfusion and ventilation; bone, liver, hepatobiliary, bone marrow, pancreas, and kidney scans; cardiac imaging procedures; tumor localization studies; localization of gastrointestinal bleeding; and non-imaging studies of blood volume and iron metabolism. Data on the relative use of radiopharmaceuticals were collected using key informant interviews and comprehensive literature reviews of typical administered activities of these diagnostic nuclear medicine studies. Responses of key informants on relative use of radiopharmaceuticals are in agreement with published literature. Results of this study will be used for retrospective reconstruction of occupational and personal medical radiation doses from diagnostic radiopharmaceuticals to members of the U.S. radiologic technologists' cohort and in reconstructing radiation doses from occupational or patient radiation exposures to other U.S. workers or patient populations.

Approaches to reducing radiation dose from radionuclide myocardial perfusion imaging.

Radionuclide myocardial perfusion imaging (MPI) plays a vital role in the evaluation and management of patients with coronary artery disease. However, because of a steep growth in MPI in the mid 2000s, concerns about inappropriate use of MPI and imaging-related radiation exposure increased. In response, the professional societies developed appropriate-use criteria for MPI. Simultaneously, novel technology, image-reconstruction software for traditional scanners, and dedicated cardiac scanners emerged and facilitated the performance of MPI with low-dose and ultra-low-dose radiotracers. This paper provides a practical approach to performing low-radiation-dose MPI using traditional and novel technologies.

Imaging cardiac amyloidosis: a pilot study using ¹8F-florbetapir positron emission tomography.

PURPOSE: Cardiac amyloidosis, a restrictive heart disease with high mortality and morbidity, is underdiagnosed due to limited targeted diagnostic imaging. The primary aim of this study was to evaluate the utility of (18)F-florbetapir for imaging cardiac amyloidosis. METHODS: We performed a pilot study of cardiac (18)F-florbetapir PET in 14 subjects: 5 control subjects without amyloidosis and 9 subjects with documented cardiac amyloidosis. Standardized uptake values (SUV) of (18)F-florbetapir in the left ventricular (LV) myocardium, blood pool, liver, and vertebral bone were determined. A (18)F-florbetapir retention index (RI) was computed. Mean LV myocardial SUVs, target-to-background ratio (TBR, myocardial/blood pool SUV ratio) and myocardial-to-liver SUV ratio between 0 and 30 min were calculated. RESULTS: Left and right ventricular myocardial uptake of (18)F-florbetapir were noted in all the amyloid subjects and in none of the control subjects. The RI, TBR, LV myocardial SUV and LV myocardial to liver SUV ratio were all significantly higher in the amyloidosis subjects than in the control subjects (RI median 0.043 min(-1), IQR 0.034 - 0.051 min(-1), vs. 0.023 min(-1), IQR 0.015 - 0.025 min(-1), P = 0.002; TBR 1.84, 1.64 - 2.50, vs. 1.26, IQR 0.91 - 1.36, P = 0.001; LV myocardial SUV 3.84, IQR 1.87 - 5.65, vs. 1.35, IQR 1.17 - 2.28, P = 0.029; ratio of LV myocardial to liver SUV 0.67, IQR 0.44 - 1.64, vs. 0.18, IQR 0.15 - 0.35, P = 0.004). The myocardial RI, TBR and myocardial to liver SUV ratio also distinguished the control subjects from subjects with transthyretin and those with light chain amyloid. CONCLUSION: (18)F-Florbetapir PET may be a promising technique to image light chain and transthyretin cardiac amyloidosis. Its role in diagnosing amyloid in other organ systems and in assessing response to therapy needs to be further studied.

Quantification of Myocardial Perfusion Reserve Using Dynamic SPECT Imaging in Humans: A Feasibility Study.

UNLABELLED: Myocardial perfusion imaging (MPI) is well established in the diagnosis and workup of patients with known or suspected coronary artery disease (CAD); however, it can underestimate the extent of obstructive CAD. Quantification of myocardial perfusion reserve with PET can assist in the diagnosis of multivessel CAD. We evaluated the feasibility of dynamic tomographic SPECT imaging and quantification of a retention index to describe global and regional myocardial perfusion reserve using a dedicated solid-state cardiac camera. METHODS: Ninety-five consecutive patients (64 men and 31 women; median age, 67 y) underwent dynamic SPECT imaging with (99m)Tc-sestamibi at rest and at peak vasodilator stress, followed by standard gated MPI. The dynamic images were reconstructed into 60-70 frames, 3-6 s/frame, using ordered-subsets expectation maximization with 4 iterations and 32 subsets. Factor analysis was used to estimate blood-pool time-activity curves, used as input functions in a 2-compartment kinetic model. K1 values ((99m)Tc-sestamibi uptake) were calculated for the stress and rest images, and K2 values ((99m)Tc-sestamibi washout) were set to zero. Myocardial perfusion reserve (MPR) index was calculated as the ratio of the stress and rest K1 values. Standard MPI was evaluated semiquantitatively, and total perfusion deficit (TPD) of at least 5% was defined as abnormal. RESULTS: Global MPR index was higher in patients with normal MPI (n = 51) than in patients with abnormal MPI (1.61 [interquartile range (IQR), 1.33-2.03] vs. 1.27 [IQR, 1.12-1.61], P = 0.0002). By multivariable regression analysis, global MPR index was associated with global stress TPD, age, and smoking. Regional MPR index was associated with the same variables and with regional stress TPD. Sixteen patients undergoing invasive coronary angiography had 20 vessels with stenosis of at least 50%. The MPR index was 1.11 (IQR, 1.01-1.21) versus 1.30 (IQR, 1.12-1.67) in territories supplied by obstructed and nonobstructed arteries, respectively (P = 0.02). MPR index showed a stepwise reduction with increasing extent of obstructive CAD (P = 0.02). CONCLUSION: Dynamic tomographic imaging and quantification of a retention index describing global and regional perfusion reserve are feasible using a solid-state camera. Preliminary results show that the MPR index is lower in patients with perfusion defects and in regions supplied by obstructed coronary arteries. Further studies are needed to establish the clinical role of this technique as an aid to semiquantitative analysis of MPI.

Evaluation of a method for projection-based tissue-activity estimation within small volumes of interest.

A new method of compensating for tissue-fraction and count-spillover effects, which require tissue segmentation only within a small volume surrounding the primary lesion of interest, was evaluated for SPECT imaging. Tissue-activity concentration estimates are obtained by fitting the measured projection data to a statistical model of the segmented tissue projections. Multiple realizations of two simulated human-torso phantoms, each containing 20 spherical 'tumours', 1.6 cm in diameter, with tumour-to-background ratios of 8:1 and 4:1, were simulated. Estimates of tumour- and background-activity concentration values for homogeneous as well as inhomogeneous tissue activities were compared to the standard uptake value (SUV) metrics on the basis of accuracy and precision. For perfectly registered, high-contrast, superficial lesions in a homogeneous background without scatter, the method yielded accurate (<0.4% bias) and precise (<6.1%) recovery of the simulated activity values, significantly outperforming the SUV metrics. Tissue inhomogeneities, greater tumour depth and lower contrast ratios degraded precision (up to 11.7%), but the estimates remained almost unbiased. The method was comparable in accuracy but more precise than a well-established matrix inversion approach, even when errors in tumour size and position were introduced to simulate moderate inaccuracies in segmentation and image registration. Photon scatter in the object did not significantly affect the accuracy or precision of the estimates.

Improved regional activity quantitation in nuclear medicine using a new approach to correct for tissue partial volume and spillover effects.

We have developed a new method of compensating for effects of partial volume and spillover in dual-modality imaging. The approach requires segmentation of just a few tissue types within a small volume-of-interest (VOI) surrounding a lesion; the algorithm estimates simultaneously, from projection data, the activity concentration within each segmented tissue inside the VOI. Measured emission projections were fitted to the sum of resolution-blurred projections of each such tissue, scaled by its unknown activity concentration, plus a global background contribution obtained by reprojection through the reconstructed image volume outside the VOI. The method was evaluated using multiple-pinhole µSPECT data simulated for the MOBY mouse phantom containing two spherical lung tumors and one liver tumor, as well as using multiple-bead phantom data acquired on µSPECT and µCT scanners. Each VOI in the simulation study was 4.8 mm (12 voxels) cubed and, depending on location, contained up to four tissues (tumor, liver, heart, lung) with different values of relative (99m)Tc concentration. All tumor activity estimates achieved bias after ∼ 15 ordered-subsets expectation maximization (OSEM) iterations (×10 subsets) , with better than 8% precision ( ≤ 25% greater than the Cramer-Rao lower bound). The projection-based fitting approach also outperformed three standardized uptake value (SUV)-like metrics, one of which was corrected for count spillover. In the bead phantom experiment, the mean ± standard deviation of the bias of VOI estimates of bead concentration were 0.9±9.5%, comparable to those of a perturbation geometric transfer matrix (pGTM) approach (-5.4±8.6%); however, VOI estimates were more stable with increasing iteration number than pGTM estimates, even in the presence of substantial axial misalignment between µCT and µSPECT image volumes.

Joint optimization of collimator and reconstruction parameters in SPECT imaging for lesion quantification.

Obtaining the best possible task performance using reconstructed SPECT images requires optimization of both the collimator and reconstruction parameters. The goal of this study is to determine how to perform this optimization, namely whether the collimator parameters can be optimized solely from projection data, or whether reconstruction parameters should also be considered. In order to answer this question, and to determine the optimal collimation, a digital phantom representing a human torso with 16 mm diameter hot lesions (activity ratio 8:1) was generated and used to simulate clinical SPECT studies with parallel-hole collimation. Two approaches to optimizing the SPECT system were then compared in a lesion quantification task: sequential optimization, where collimation was optimized on projection data using the Cramer­Rao bound, and joint optimization, which simultaneously optimized collimator and reconstruction parameters. For every condition, quantification performance in reconstructed images was evaluated using the root-mean-squared-error of 400 estimates of lesion activity. Compared to the joint-optimization approach, the sequential-optimization approach favoured a poorer resolution collimator, which, under some conditions, resulted in sub-optimal estimation performance. This implies that inclusion of the reconstruction parameters in the optimization procedure is important in obtaining the best possible task performance; in this study, this was achieved with a collimator resolution similar to that of a general-purpose (LEGP) collimator. This collimator was found to outperform the more commonly used high-resolution (LEHR) collimator, in agreement with other task-based studies, using both quantification and detection tasks.

Improved activity estimation with MC-JOSEM versus TEW-JOSEM in 111In SPECT.

We have previously developed a fast Monte Carlo (MC)-based joint ordered-subset expectation maximization (JOSEM) iterative reconstruction algorithm, MC-JOSEM. A phantom study was performed to compare quantitative imaging performance of MC-JOSEM with that of a triple-energy-window approach (TEW) in which estimated scatter was also included additively within JOSEM, TEW-JOSEM. We acquired high-count projections of a 5.5 cm3 sphere of 111In at different locations in the water-filled torso phantom; high-count projections were then obtained with 111In only in the liver or only in the soft-tissue background compartment, so that we could generate synthetic projections for spheres surrounded by various activity distributions. MC scatter estimates used by MC-JOSEM were computed once after five iterations of TEW-JOSEM. Images of different combinations of liver/background and sphere/background activity concentration ratios were reconstructed by both TEW-JOSEM and MC-JOSEM for 40 iterations. For activity estimation in the sphere, MC-JOSEM always produced better relative bias and relative standard deviation than TEW-JOSEM for each sphere location, iteration number, and activity combination. The average relative bias of activity estimates in the sphere for MC-JOSEM after 40 iterations was -6.9%, versus -15.8% for TEW-JOSEM, while the average relative standard deviation of the sphere activity estimates was 16.1% for MC-JOSEM, versus 27.4% for TEW-JOSEM. Additionally, the average relative bias of activity concentration estimates in the liver and the background for MC-JOSEM after 40 iterations was -3.9%, versus -12.2% for TEW-JOSEM, while the average relative standard deviation of these estimates was 2.5% for MC-JOSEM, versus 3.4% for TEW-JOSEM. MC-JOSEM is a promising approach for quantitative activity estimation in 111In SPECT.

Adsorption of metallic radionuclides on plastic phantom walls.

Metal radionuclide solutions at neutral pH adhere to plastic containers. Adsorption of radionuclides on the walls of phantoms leads to a nonuniform activity distribution, which could adversely affect imaging studies, as well as phantom-based validations of absorbed dose calculations used in radioimmunotherapy, requiring accurate knowledge of the underlying activity distribution. In the work reported here, the authors determined the degree of metal chelation required to minimize metallic radionuclide oxide formation and adsorption on phantom walls in order to yield more reliable experimental data for validating image-based dosimetry. Using hollow spherical plastic phantoms, the authors evaluated three different radionuclides, I-131, In-111, and Y-90, in solutions containing three different concentrations of the chelator, ethylenediaminetetraacetate (EDTA). Adsorption to plastic walls was determined using microSPECT imaging and/or by counting aliquots of solutions. Reconstructed images and measurements of I-131 activity showed that it was uniformly distributed within all spheres; however, images of In-111 in 0.25-microM EDTA indicated that the activity concentration near the wall was much higher than that in the middle of the sphere. The decrease in activity concentration near the center of the spheres was approximately 47%. Y-90 in 0.25-microM EDTA behaved similarly; the activity concentration of Y-90 decreased by 46%. For an In-111 or Y-90 radioactivity concentration of 0.74 MBq/mL, a 2.5-microM EDTA solution was required to achieve a uniform distribution, suggesting that, under our experimental conditions, approximately 700 EDTA molecules were required for each radiometal atom to prevent precipitation and adsorption on poly(methylmethacrylate). For certain radiometals, e.g., In-111 or Y-90, adequate chelation is essential to achieve uniform activity concentration values and homogeneous distribution within the phantom compartments.

Quantitative simultaneous In-111/Tc-99m planar imaging in a long-bone infection phantom.

In-111-white-blood-cell and Tc-99m-sulfur-colloid dual-radionuclide imaging are frequently utilized in the evaluation of patients with suspected osteomyelitis. We have developed a quantitative planar imaging method in which Tc-99m and In-111 scans are acquired simultaneously in accurate spatial registration. Long, thin tubes containing only In-111 or Tc-99m were first imaged in a list mode within a water bath inclined with respect to the water surface; from these, 12 energy spectra corresponding to different Tc/In ratios were synthesized. Triple-energy-window (TEW) parameters for scatter and radionuclide crosstalk correction, including scatter windows and weights, were optimized using 100 noise realizations of each of the spectra (1200 total). A long-bone phantom containing a simulated infection site was then imaged in water with five In/Tc ratios; 100 noise realizations of two conjugate-view images were generated from each acquisition (500 total). Two regions of interest (ROIs) were defined, and the ratio of In/Tc count ratios in these two ROIs was evaluated with and without the TEW scatter correction and geometric mean attenuation compensation. The average bias improved from 17.2% to 5.3%, with comparable precision. TEW corrections with non-optimized but practical energy windows also improved the bias to 6.4%. Compared with subjective visual assessment, quantitation of In-111/Tc-99m ratios may improve diagnostic accuracy and could eventually permit grading of osteomyelitis.

Collimator optimization for detection and quantitation tasks: application to gallium-67 imaging.

We describe a new approach to the problem of collimator optimization in nuclear medicine; our methodology is illustrated for the challenging case of gallium-67 imaging. Collimator-design methods based on empirical rules, such as specification of an allowable level of single-septal penetration (SSP) at a fixed energy, are especially inappropriate for radionuclides characterized by an abundance of high-energy contaminant photons that scatter in the patient, collimator, and/or detector before detection within one of a few photopeak energy windows. Lead X-rays produced in the collimator are an additional source of contamination. We designed optimal collimation for 67Ga based on relevant clinical imaging tasks and a realistic simulation of photon transport in a phantom, collimator, and detector. Collimator designs were compared on the basis of performance in lesion detection, as predicted by a three-channel Hotelling observer (CHO), as well as in tumor and background activity estimation (EST), quantified by task-specific signal-to-noise ratios (SNRs). The optimal values of collimator lead content were 22.0 and 23.8 g/cm2, respectively, for CHO and EST, while the optimal geometric resolution values were 1.8 and 1.6 cm full-width at half-maximum (FWHM), respectively, at a distance of 23.5 cm. The resolution of a commercially available medium-energy low-penetration collimator (MELP) is 1.9 cm FWHM at this distance. The optimal values for SSP at 300 keV were 7.3% and 5.8% based on CHO and EST, respectively, compared to 5.2% for the MELP collimator. Compared with the commercial MELP collimator, the 67Ga collimator optimized for tumor detection or activity estimation tasks provided improved geometric spatial resolution with reduced geometric efficiency and, surprisingly, allowed an increased level of single-septal penetration.

Optimization of Ga-67 imaging for detection and estimation tasks: dependence of imaging performance on spectral acquisition parameters.

UNLABELLED: We have compared the use of two (93 and 185 keV) and three (93, 185, and 300 keV) photopeaks for Ga-67 tumor imaging and optimized the placement of each energy window. METHODS: The bases for optimization and evaluation were ideal and Bayesian signal-to-noise ratios (SNR) for the detection of spheres embedded in a realistic anthropomorphic digital torso phantom and ideal SNR for the estimation of their size and activity concentration. Seven spheres of radii ranging from 1 to 3 cm, located at several sites in the torso, were simulated using a realistic Monte Carlo program. We also calculated the ideal SNR for the detection from simple phantom acquisitions. RESULTS: For detection and estimation tasks, the optimum windows were identical for all sphere sizes and locations. For the 93 keV photopeak, the optimal window was 84-102 keV for the detection and 87-102 keV for estimation; these windows are narrower than the 20% window often used in the clinic (83-101 keV). For the 185 keV photopeak, the optimal window was 170-220 keV for the detection and 170-215 keV for estimation; these are substantially different than the 15% window used in our clinic (171-199 keV). For the 300 keV photopeak, the optimal window for detection was 270-320 keV, and for estimation, 280-320 keV. Using the three optimized, rather than only the two lower-energy, windows yielded a 9% increase in the SNR for the detection of the 3 cm diam sphere (a 12% increase for a 2 cm diam sphere) and a 7% increase in the SNR for estimation of its size. For the acquired phantom data, detection also increased by 9%-12% when using three, rather than two, energy windows.