Jazmin Schwartz

An iterative technique to segment PET lesions using a Monte Carlo based mathematical model

PURPOSE The need for an accurate lesion segmentation tool in 18FDG PET is a prerequisite for the estimation of lesion response to therapy, for radionuclide dosimetry, and for the application of 18FDG PET to radiotherapy planning. In this work, the authors have developed an iterative method based on a mathematical fit deduced from Monte Carlo simulations to estimate tumor segmentation thresholds. METHODS The GATE software, a GEANT4 based Monte Carlo tool, was used to model the GE Advance PET scanner geometry. Spheres ranging between 1 and 6 cm in diameters were simulated in a 10 cm high and 11 cm in diameter cylinder. The spheres were filled with water-equivalent density and simulated in both water and lung equivalent background. The simulations were performed with an infinite, 8/1, and 4/1 target-to-background ratio (T/B). A mathematical fit describing the correlation between the lesion volume and the corresponding optimum threshold value was then deduced through analysis of the reconstructed images. An iterative method, based on this mathematical fit, was developed to determine the optimum threshold value. The effects of the lesion volume and T/B on the threshold value were investigated. This method was evaluated experimentally using the NEMA NU2-2001 IEC phantom, the ACNP cardiac phantom, a randomly deformed aluminum can, and a spheroidal shape phantom implemented artificially in the lung, liver, and brain of patient PET images. Clinically, the algorithm was evaluated in six lesions from five patients. Clinical results were compared to CT volumes. RESULTS This mathematical fit predicts an existing relationship between the PET lesion size and the percent of maximum activity concentration within the target volume (or threshold). It also showed a dependence of the threshold value on the T/B, which could be eliminated by background subtraction. In the phantom studies, the volumes of the segmented PET targets in the PET images were within 10% of the nominal ones. Clinically, the PET target volumes were also within 10% of those measured from CT images. CONCLUSIONS This iterative algorithm enabled accurately segment PET lesions, independently of their contrast value.

Study of radiolabeled indium-111 and yttrium-90 ibritumomab tiuxetan in primary central nervous system lymphoma.

Therapeutic options for refractory or recurrent primary central nervous system lymphoma (PCNSL) are limited. The blood–brain barrier makes many agents used in systemic lymphomas ineffective in CNS lymphomas. The objective of this study was to determine whether intravenous radioimmunotherapy using anti‐CD20 antibody can be delivered to PCNSL.

Dynamic Small-Animal PET Imaging of Tumor Proliferation with 3′-Deoxy-3′-18F-Fluorothymidine in a Genetically Engineered Mouse Model of High-Grade Gliomas

3′-Deoxy-3′-18F-fluorothymidine (18F-FLT), a partially metabolized thymidine analog, has been used in preclinical and clinical settings for the diagnostic evaluation and therapeutic monitoring of tumor proliferation status. We investigated the use of 18F-FLT for detecting and characterizing genetically engineered mouse (GEM) high-grade gliomas and evaluating the pharmacokinetics in GEM gliomas and normal brain tissue. Our goal was to develop a robust and reproducible method of kinetic analysis for the quantitative evaluation of tumor proliferation. Methods: Dynamic 18F-FLT PET imaging was performed for 60 min in glioma-bearing mice (n = 10) and in non–tumor-bearing control mice (n = 4) by use of a dedicated small-animal PET scanner. A 3-compartment, 4-parameter model was used to characterize 18F-FLT kinetics in vivo. For compartmental analysis, the arterial input was measured by placing a region of interest over the left ventricular blood pool and was corrected for partial-volume averaging. The 18F-FLT “trapping” and tissue flux model parameters were correlated with measured uptake (percentage injected dose per gram [%ID/g]) values at 60 min. Results: 18F-FLT uptake values (%ID/g) at 1 h in brain tumors were significantly greater than those in control brains (mean ± SD: 4.33 ± 0.58 and 0.86 ± 0.22, respectively; P < 0.0004). Kinetic analyses of the measured time–activity curves yielded independent, robust estimates of tracer transport and metabolism, with compartmental model–derived time–activity data closely fitting the measured data. Except for tracer transport, statistically significant differences were found between the applicable model parameters for tumors and normal brains. The tracer retention rate constant strongly correlated with measured 18F-FLT uptake values (r = 0.85, P < 0.0025), whereas a more moderate correlation was found between net 18F-FLT flux and 18F-FLT uptake values (r = 0.61, P < 0.02). Conclusion: A clinically relevant mouse glioma model was characterized by both static and dynamic small-animal PET imaging of 18F-FLT uptake. Time–activity curves were kinetically modeled to distinguish early transport from a subsequent tracer retention phase. Estimated 18F-FLT rate constants correlated positively with %ID/g measurements. Dynamic evaluation of 18F-FLT uptake offers a promising approach for noninvasively assessing cellular proliferation in vivo and for quantitatively monitoring new antiproliferation therapies.

Renal uptake of bismuth-213 and its contribution to kidney radiation dose following administration of actinium-225-labeled antibody

Clinical therapeutic studies using (225)Ac-labeled antibodies have begun. Of major concern is renal toxicity that may result from the three alpha-emitting progeny generated following the decay of (225)Ac. The purpose of this study was to determine the amount of (225)Ac and non-equilibrium progeny in the mouse kidney after the injection of (225)Ac-huM195 antibody and examine the dosimetric consequences. Groups of mice were sacrificed at 24, 96 and 144 h after injection with (225)Ac-huM195 antibody and kidneys excised. One kidney was used for gamma ray spectroscopic measurements by a high-purity germanium (HPGe) detector. The second kidney was used to generate frozen tissue sections which were examined by digital autoradiography (DAR). Two measurements were performed on each kidney specimen: (1) immediately post-resection and (2) after sufficient time for any non-equilibrium excess (213)Bi to decay completely. Comparison of these measurements enabled estimation of the amount of excess (213)Bi reaching the kidney (γ-ray spectroscopy) and its sub-regional distribution (DAR). The average absorbed dose to whole kidney, determined by spectroscopy, was 0.77 (SD 0.21) Gy kBq(-1), of which 0.46 (SD 0.16) Gy kBq(-1) (i.e. 60%) was due to non-equilibrium excess (213)Bi. The relative contributions to renal cortex and medulla were determined by DAR. The estimated dose to the cortex from non-equilibrium excess (213)Bi (0.31 (SD 0.11) Gy kBq(-1)) represented ∼46% of the total. For the medulla the dose contribution from excess (213)Bi (0.81 (SD 0.28) Gy kBq(-1)) was ∼80% of the total. Based on these estimates, for human patients we project a kidney-absorbed dose of 0.28 Gy MBq(-1) following administration of (225)Ac-huM195 with non-equilibrium excess (213)Bi responsible for approximately 60% of the total. Methods to reduce renal accumulation of radioactive progeny appear to be necessary for the success of (225)Ac radioimmunotherapy.

Improved tumor imaging and therapy via i.v. IgG-mediated time-sequential modulation of neonatal Fc receptor.

The long plasma half-life of IgG, while allowing for enhanced tumor uptake of tumor-targeted IgG conjugates, also results in increased background activity and normal-tissue toxicity. Therefore, successful therapeutic uses of conjugated antibodies have been limited to the highly sensitive and readily accessible hematopoietic tumors. We report a therapeutic strategy to beneficially alter the pharmacokinetics of IgG antibodies via pharmacological inhibition of the neonatal Fc receptor (FcRn) using high-dose IgG therapy. IgG-treated mice displayed enhanced blood and whole-body clearance of radioactivity, resulting in better tumor-to-blood image contrast and protection of normal tissue from radiation. Tumor uptake and the resultant therapeutic response was unaltered. Furthermore, we demonstrated the use of this approach for imaging of tumors in humans and discuss its potential applications in cancer imaging and therapy. The ability to reduce the serum persistence of conjugated IgG antibodies after their infusion can enhance their therapeutic index, resulting in improved therapeutic and diagnostic efficacy.

Bone Marrow Dosimetry Using 124I-PET

Bone marrow is usually dose-limiting for radioimmunotherapy. In this study, we directly estimated red marrow activity concentration and the self-dose component of absorbed radiation dose to red marrow based on PET/CT of 2 different 124I-labeled antibodies (cG250 and huA33) and compared the results with plasma activity concentration and plasma-based dose estimates. Methods: Two groups of patients injected with 124I-labeled monoclonal antibodies (11 patients with renal cancer receiving 124I-cG250 and 5 patients with colorectal cancer receiving 124I- huA33) were imaged by PET or PET/CT on 2 or 3 occasions after infusion. Regions of interest were drawn over several lumbar vertebrae, and red marrow activity concentration was quantified. Plasma activity concentration was also quantified using multiple patient blood samples. The red marrow–to–plasma activity concentration ratio (RMPR) was calculated at the times of imaging. The self-dose component of the absorbed radiation dose to the red marrow was estimated from the images, from the plasma measurements, and using a combination of both sets of measurements. Results: RMPR was observed to increase with time for both groups of patients. Mean (±SD) time-dependent RMPR (RMPR(t)) for the cG250 group increased from 0.13 ± 0.06 immediately after infusion to 0.23 ± 0.09 at approximately 6 d after infusion. For the huA33 group, mean RMPR(t) was 0.10 ± 0.04 immediately after infusion, 0.13 ± 0.05 approximately 2 d after infusion, and 0.20 ± 0.09 approximately 7 d after infusion. Plasma-based estimates of red marrow self-dose tended to be greater than image-based values by, on average, 11% and 47% for cG250 and huA33, respectively, but by as much as −73% to 62% for individual patients. The hybrid method combining RMPR(t) and plasma activity concentration provided a closer match to the image-based dose estimates (average discrepancies, −2% and 18% for cG250 and huA33, respectively). Conclusion: These results suggest that the assumption of time-independent proportionality between red marrow and plasma activity concentration may be too simplistic. Individualized imaged-based dosimetry is probably required for the optimal therapeutic delivery of radiolabeled antibodies, which does not compromise red marrow and may allow, for some patients, a substantial increase in administered activity and thus tumor dose.

Repeatability of SUV measurements in serial PET

PURPOSE The standardized uptake value (SUV) is a quantitative measure of FDG tumor uptake frequently used as a tool to monitor therapeutic response. This study aims to (i) assess the reproducibility and uncertainty of SUV max and SUV mean, due to purely statistical, i.e., nonbiological, effects and (ii) to establish the minimum uncertainty below which changes in SUV cannot be expected to be an indicator of physiological changes. METHODS Three sets of measurements were made using a GE Discovery STE PET/CT Scanner in 3D mode: (1) A uniform 68Ge 20 cm diameter cylindrical phantom was imaged. Thirty serial frames were acquired for durations of 3, 6, 10, 15, and 30 min. (2) Esser flangeless phantom (Data Spectrum, approximately 6.1 L) with fillable thin-walled cylinders inserts (diameters: 8, 12, 16, and 25 mm; height: approximately 3.8 mm) was scanned for five consecutive 3 min runs. The cylinders were filled with 18FDG with a 37 kBq/cc concentration, and with a target-to-background ratio (T/BKG) of 3/1. (3) Eight cancer patients with healthy livers were scanned approximately 1.5 h post injection. Three sequential 3 min scans were performed for one bed position covering the liver, with the patient and bed remaining at the same position for the entire length of the scan. Volumes of interest were drawn on all images using the corresponding CT and then transferred to the PET images. For each study (1-3), the average percent change in SUV mean and SUV max were determined for each run pair. Moreover, the repeatability coefficient was calculated for both the SUV mean and SUV max for each pair of runs. Finally, the overall ROI repeatability coefficient was determined for each pair of runs. RESULTS For the 68Ge phantom the average percent change in SUV max and SUV mean decrease as a function of increasing acquisition time from 4.7 +/- 3.1 to 1.1 +/- 0.6%, and from 0.14 +/- 0.09 to 0.04 +/- 0.03%, respectively. Similarly, the coefficients of repeatability also decrease between the 3 and 30 min acquisition scans, in the range of 10.9 +/- 3.9% - 2.6 +/- 0.9%, and 0.3 +/- 0.1% - 0.10 +/- 0.04%, for the SUV max and SUV mean, respectively. The overall ROI repeatability decreased from 18.9 +/- 0.2 to 6.0 +/- 0.1% between the 3 and 30 min acquisition scans. For the l8FDG phantom, the average percent change in SUV max and SUV mean decreases with target diameter from 3.6 +/- 2.0 to 1.5 +/- 0.8% and 1.5 +/- 1.3 to 0.26 +/- 0.15%, respectively, for targets from 8-25 mm in diameter and for a region in the background (BKG). The coefficients of repeatability for SUV max and SUV mean also decrease as a function of target diameter from 7.1 +/- 2.5 to 2.4 +/- 0.9 and 4.2 +/- 1.5 to 0.6 +/- 0.2, respectively, for targets from 8 mm to BKG in diameter. Finally, overall ROI repeatability decreased from 12.0 +/- 4.1 to 13.4 +/- 0.5 targets from 8 mm to BKG in diameter. Finally, for the measurements in healthy livers the average percent change in SUVmax and SUV mean were in the range of 0.5 +/- 0.2% - 6.2 +/- 3.9% and 0.4 +/- 0.1 and 1.6 +/- 1%, respectively. The coefficients of repeatability for SUV max and SUV men are in the range of 0.6 +/- 0.7% - 9.5 +/- 12% and 0.6 +/- 0.7% - 2.9 +/- 3.6%, respectively. The overall target repeatability varied between 27.9 +/- 0.5% and 41.1 +/- 1.0%. CONCLUSIONS The statistical fluctuations of the SUV mean are half as large as those of the SUV max in the absence of biological or physiological effects. In addition, for clinically applicable scan durations (i.e., approximately 3 min) and FDG concentrations, the SUV max and SUV mean have similar amounts of statistical fluctuation for small regions. However, the statistical fluctuations of the SUVmean rapidly decrease with respect tothe SUVmax as the statistical power of the data grows either due to longer scanning times or as the target regions encompass a larger volume.

Pharmacokinetic Analysis of Dynamic 18 F-Fluoromisonidazole PET Data in Non–Small Cell Lung Cancer

Hypoxic tumors exhibit increased resistance to radiation, chemical, and immune therapies. 18F-fluoromisonidazole (18F-FMISO) PET is a noninvasive, quantitative imaging technique used to evaluate the magnitude and spatial distribution of tumor hypoxia. In this study, pharmacokinetic analysis (PKA) of 18F-FMISO dynamic PET extended to 3 h after injection is reported for the first time, to our knowledge, in stage III–IV non–small cell lung cancer (NSCLC) patients. Methods: Sixteen patients diagnosed with NSCLC underwent 2 PET/CT scans (1–3 d apart) before radiation therapy: a 3-min static 18F-FDG and a dynamic 18F-FMISO scan lasting 168 ± 15 min. The latter data were acquired in 3 serial PET/CT dynamic imaging sessions, registered with each other and analyzed using pharmacokinetic modeling software. PKA was performed using a 2-tissue, 3-compartment irreversible model, and kinetic parameters were estimated for the volumes of interest determined using coregistered 18F-FDG images for both the volume of interest–averaged and the voxelwise time–activity curves for each patient’s lesions, normal lung, and muscle. Results: We derived average values of 18F-FMISO kinetic parameters for NSCLC lesions as well as for normal lung and muscle. We also investigated the correlation between the trapping rate (k3) and delivery rate (K1), influx rate (Ki) constants, and tissue-to-blood activity concentration ratios (TBRs) for all tissues. Lesions had trapping rates 1.6 times larger, on average, than those of normal lung and 4.4 times larger than those in muscle. Additionally, for almost all cases, k3 and Ki had a significant strong correlation for all tissue types. The TBR–k3 correlation was less straightforward, showing a moderate to strong correlation for only 41% of lesions. Finally, K1–k3 voxelwise correlations for tumors were varied, but negative for 76% of lesions, globally exhibiting a weak inverse relationship (average R = −0.23 ± 0.39). However, both normal tissue types exhibited significant positive correlations for more than 60% of patients, with 41% having moderate to strong correlations (R > 0.5). Conclusion: All lesions showed distinct 18F-FMISO uptake. Variable 18F-FMISO delivery was observed across lesions, as indicated by the variable values of the kinetic rate constant K1. Except for 3 cases, some degree of hypoxia was apparent in all lesions based on their nonzero k3 values.

Feasibility of 18F-Fluoromisonidazole Kinetic Modeling in Head and Neck Cancer Using Shortened Acquisition Times

18F-fluoromisonidazole dynamic PET (dPET) is used to identify tumor hypoxia noninvasively. Its routine clinical implementation, however, has been hampered by the long acquisition times required. We investigated the feasibility of kinetic modeling using shortened acquisition times in 18F-fluoromisonidazole dPET, with the goal of expediting the clinical implementation of 18F-fluoromisonidazole dPET protocols. Methods: Six patients with squamous cell carcinoma of the head and neck and 10 HT29 colorectal carcinoma–bearing nude rats were studied. In addition to an 18F-FDG PET scan, each patient underwent a 45-min 18F-fluoromisonidazole dPET scan, followed by 10-min acquisitions at 96 ± 4 and 163 ± 17 min after injection. Ninety-minute 18F-fluoromisonidazole dPET scans were acquired in animals. Intratumor voxels were classified into 4 clusters based on their kinetic behavior using k-means clustering. Kinetic modeling was performed using the foregoing full datasets (FD) and repeated for each of 2 shortened datasets corresponding to the first approximately 100 min (SD1; patients only) or the first 45 min (SD2) of dPET data. The kinetic rate constants (KRCs) as calculated with a 2-compartment model for both SD1 and SD2 were compared with those derived from FD by correlation (Pearson), regression (Passing–Bablok), deviation (Bland–Altman), and classification (area-under-the-receiver-operating characteristic curve) analyses. Simulations were performed to assess uncertainties due to statistical noise. Results: Strong correlation (r ≥ 0.75, P < 0.001) existed between all KRCs deduced from both SD1 and SD2, and from FD. Significant differences between KRCs were found only for FD-SD2 correlations in patient studies. K1 and k3 were reproducible to within approximately 6% and approximately 30% (FD-SD1; patients) and approximately 4% and approximately 75% (FD-SD2; animals). Area-under-the-receiver-operating characteristic curve values for classification of patient clusters as hypoxic, using a tumor-to-blood ratio greater than 1.2, were 0.91 (SD1) and 0.86 (SD2). The percentage SD in estimating K1 and k3 from 45-min shortened datasets due to noise was less than 1% and between 2% and 12%, respectively. Conclusion: Using single-session 45-min shortened 18F-fluoromisonidazole dPET datasets appears to be adequate for the identification of intratumor regions of hypoxia. However, k3 was significantly overestimated in the clinical cohort. Further studies are necessary to evaluate the clinical significance of differences between the results as calculated from full and shortened datasets.

Inter-operator variability in compartmental kinetic analysis of 18 F-fluoromisonidazole dynamic PET

PURPOSE To assess the inter-operator variability in compartment analysis (CA) of dynamic-FMISO (dyn-FMISO) PET. METHODS Study-I: Five investigators conducted CA for 23 NSCLC dyn-FMISO tumor time-activity-curves. Study-II: Four operators performed CA for four NSCLC dyn-FMISO datasets. Repeatability of Kinetic-Rate-Constants (KRCs) was assessed. RESULTS Study-I: Strong correlation (ICC > 0.9) and interchangeable results among operators existed for all KRCs. Study-II: Up to 103% variability in tumor segmentation, and weaker ICC in KRCs (ICC-VB = 0.53; ICC-K1 = 0.91; ICC-K1/k2 = 0.25; ICC-k3 = 0.32; ICC-Ki = 0.54) existed. All KRCs were repeatable among the different operators. CONCLUSIONS Inter-operator variability in CA of dyn-FMISO was shown to be within statistical errors.