Kidney dosimetry in [177Lu]Lu-DOTA-TATE therapy based on multiple small VOIs

PURPOSE

The aim was to investigate the use of multiple small VOIs for kidney dosimetry in [177Lu]Lu-DOTA-TATE therapy.

METHOD

The study was based on patient and simulated SPECT images in anthropomorphic geometries. Images were reconstructed using two reconstruction programs (local LundaDose and commercial Hermia) using OS-EM with and without resolution recovery (RR). Five small VOIs were placed to determine the average activity concentration (AC) in each kidney. The study consisted of three steps: (i) determination of the number of iterations for AC convergence based on simulated images; (ii) determination of recovery-coefficients (RCs) for 2 mL VOIs using a separate set of simulated images; (iii) assessment of operator variability in AC estimates for simulated and patient images. Five operators placed the VOIs, using for guidance: a) SPECT/CT with RR, b) SPECT/CT without RR, and c) CT only. For simulated images, time-integrated ACs (TIACs) were evaluated. For patient images, estimated ACs were compared with results of a previous method based on whole-kidney VOIs.

RESULTS

Eight iterations and ten subsets were sufficient for both programs and reconstruction settings. Mean RCs (mean ± SD) with RR were 1.03 ± 0.02 (LundaDose) and 1.10 ± 0.03 (Hermia), and without RR 0.91 ± 0.03 (LundaDose) and 0.94 ± 0.03 (Hermia). Most stable and accurate estimates of the AC were obtained using five 2-mL VOIs guided by SPECT/CT with RR, applying them to images without RR, and including an explicit RC for recovery correction.

CONCLUSION

The small VOI method based on five 2-mL VOIs was found efficient and sufficiently accurate for kidney dosimetry in [177Lu]Lu-DOTA-TATE therapy.

A PBPK model for PRRT with [177Lu]Lu-DOTA-TATE: Comparison of model implementations in SAAM II and MATLAB/SimBiology

Physiologically based pharmacokinetic (PBPK) models offer the ability to simulate and predict the biodistribution of radiopharmaceuticals and have the potential to enable individualised treatment planning in molecular radiotherapy. The objective of this study was to develop and implement a whole-body compartmental PBPK model for peptide receptor radionuclide therapy (PRRT) with [177Lu]Lu-DOTA-TATE in SimBiology to allow for more complex analyses. The correctness of the model implementation was ensured by comparing its outputs, such as the time-integrated activity (TIA), with those of a PBPK model implemented in SAAM II software.

METHODS

A combined PBPK model for [68Ga]Ga-DOTA-TATE and [177Lu]Lu-DOTA-TATE was developed and implemented in both SAAM II and SimBiology. A retrospective analysis of 12 patients with metastatic neuroendocrine tumours (NETs) was conducted. First, time-activity curves (TACs) and TIAs from the two software were calculated and compared for identical parameter values. Second, pharmacokinetic parameters were fitted to activity concentrations, analysed and compared.

RESULTS

The PBPK model implemented in SimBiology produced TIA results comparable to those generated by the model implemented in SAAM II, with a relative deviation of less than 0.5% when using the same input parameters. The relative deviation of the fitted TIAs was less than 5% when model parameter values were fitted to the measured activity concentrations.

CONCLUSION

The proposed PBPK model implemented in SimBiology can be used for dosimetry in radioligand therapy and TIA prediction. Its outputs are similar to those generated by the PBPK model implemented in SAAM II, confirming the correctness of the model implementation in SimBiology.

Feasibility of 177Lu activity quantification using a small portable CZT-based gamma-camera

BACKGROUND

In image processing for activity quantification, the end goal is to produce a metric that is independent of the measurement geometry. Photon attenuation needs to be accounted for and can be accomplished utilizing spectral information, avoiding the need of additional image acquisitions. The aim of this work is to investigate the feasibility of 177Lu activity quantification with a small CZT-based hand-held gamma-camera, using such an attenuation correction method.

METHODS

A previously presented dual photopeak method, based on the differential attenuation for two photon energies, is adapted for the three photopeaks at 55 keV, 113 keV, and 208 keV for 177Lu. The measurement model describes the count rates in each energy window as a function of source depth and activity, accounting for distance-dependent system sensitivity, attenuation, and build-up. Parameter values are estimated from characterizing measurements, and the source depth and activity are obtained by minimizing the difference between measured and modelled count rates. The method is applied and evaluated in phantom measurements, in a clinical setting for superficial lesions in two patients, and in a pre-clinical setting for one human tumour xenograft. Evaluation is made for a LEHR and an MEGP collimator.

RESULTS

For phantom measurements at clinically relevant depths, the average (and standard deviation) in activity errors are 17% ± 9.6% (LEHR) and 2.9% ± 3.6% (MEGP). For patient measurements, deviations from activity estimates from planar images from a full-sized gamma-camera are 0% ± 21% (LEHR) and 16% ± 18% (MEGP). For mouse measurements, average deviations of - 16% (LEHR) and - 6% (MEGP) are obtained when compared to a small-animal SPECT/CT system. The MEGP collimator appears to be better suited for activity quantification, yielding a smaller variability in activity estimates, whereas the LEHR results are more severely affected by septal penetration.

CONCLUSIONS

Activity quantification for 177Lu using the hand-held camera is found to be feasible. The readily available nature of the hand-held camera may enable more frequent activity quantification in e.g., superficial structures in patients or in the pre-clinical setting.

Implementation of dosimetry for molecular radiotherapy; results from a European survey

PURPOSE

The use of molecular radiotherapy (MRT) has been rapidly evolving over the last years. The aim of this study was to assess the current implementation of dosimetry for MRTs in Europe.

METHODS

A web-based questionnaire was open for treating centres between April and June 2022, and focused on 2020-2022. Questions addressed the application of 16 different MRTs, the availability and involvement of medical physicists, software used, quality assurance, as well as the target regions for dosimetry, whether treatment planning and/or verification were performed, and the dosimetric methods used.

RESULTS

A total of 173 responses suitable for analysis was received from centres performing MRT, geographically distributed over 27 European countries. Of these, 146 centres (84 %) indicated to perform some form of dosimetry, and 97 % of these centres had a medical physicist available and almost always involved in dosimetry. The most common MRTs were 131I-based treatments for thyroid diseases and thyroid cancer, and [223Ra]RaCl2 for bone metastases. The implementation of dosimetry varied widely between therapies, from almost all centres performing dosimetry-based planning for microsphere treatments to none for some of the less common treatments (like 32P sodium-phosphate for myeloproliferative disease and [89Sr]SrCl2 for bone metastases).

CONCLUSIONS

Over the last years, implementation of dosimetry, both for pre-therapeutic treatment planning and post-therapy absorbed dose verification, increased for several treatments, especially for microsphere treatments. For other treatments that have moved from research to clinical routine, the use of dosimetry decreased in recent years. However, there are still large differences both across and within countries.

The contest between internal and external-beam dosimetry: The Zeno's paradox of Achilles and the tortoise

Radionuclide therapy, also called molecular radiotherapy (MRT), has come of age, with several novel radiopharmaceuticals being approved for clinical use or under development in the last decade. External beam radiotherapy (EBRT) is a well-established treatment modality, with about half of all oncologic patients expected to receive at least one external radiation treatment over their disease course. The efficacy and the toxicity of both types of treatment rely on the interaction of radiation with biological tissues. Dosimetry played a fundamental role in the scientific and technological evolution of EBRT, and absorbed doses to the target and to the organs at risk are calculated on a routine basis. In contrast, in MRT the usefulness of internal dosimetry has long been questioned, and a structured path to include absorbed dose calculation is missing. However, following a similar route of development as EBRT, MRT treatments could probably be optimized in a significant proportion of patients, likely based on dosimetry and radiobiology. In the present paper we describe the differences and the similarities between internal and external-beam dosimetry in the context of radiation treatments, and we retrace the main stages of their development over the last decades.

Traceable calibration with 177Lu and comparison of activity meters at hospitals in Norway and Sweden

INTRODUCTION

The activity meter is used to determine the activity of delivered radiopharmaceuticals, administered activity to patients and reference activity when gamma-cameras are calibrated prior to imaged-based dosimetry. The aim is to describe a procedure for cross-calibration of activity meters at different clinical sites, and report on 177Lu activity results when using factory-set calibration factors compared to when calibration is performed with traceability to a primary standard.

METHODS

Thirty activity meters placed at seven hospitals in Norway and Sweden from four manufacturers: Capintec, Commecer, NuviaTech and Veenstra were included. A stock solution with 177Lu was prepared at the local sites and measured in each activity meter with factory settings. The solution was shipped to the reference site at Lund University for measurements in a secondary standard activity meter. Deviations between local and reference activity measurements were determined for three geometries: 25-mL vial, 10-mL syringe and 1-mL syringe.

RESULTS

The median of the deviations was 6.4 % for the 25 mL vial, 5.9 % for the 10 mL syringe and 6.8 % for the 1 mL syringe. The median of the deviations for the 25 mL vial, was 1.5 % for activity meters from Capintec, 7.0 % for Comecer, 11.0 % for NuviaTech and 2.4 % for Veenstra. The majority of the deviations were positive and the maximum deviation was 14.5 %.

CONCLUSION

The activity of 177Lu measured in an activity meter with factory-set dial settings may yield deviations up to 14.5%, compared to activities measured with traceability to a primary standard. This would imply an undertreatment of patients.

On the use of solid 133Ba sources as surrogate for liquid 131I in SPECT/CT calibration: a European multi-centre evaluation

INTRODUCTION

Commissioning, calibration, and quality control procedures for nuclear medicine imaging systems are typically performed using hollow containers filled with radionuclide solutions. This leads to multiple sources of uncertainty, many of which can be overcome by using traceable, sealed, long-lived surrogate sources containing a radionuclide of comparable energies and emission probabilities. This study presents the results of a quantitative SPECT/CT imaging comparison exercise performed within the MRTDosimetry consortium to assess the feasibility of using 133Ba as a surrogate for 131I imaging.

MATERIALS AND METHODS

Two sets of four traceable 133Ba sources were produced at two National Metrology Institutes and encapsulated in 3D-printed cylinders (volume range 1.68-107.4 mL). Corresponding hollow cylinders to be filled with liquid 131I and a mounting baseplate for repeatable positioning within a Jaszczak phantom were also produced. A quantitative SPECT/CT imaging comparison exercise was conducted between seven members of the consortium (eight SPECT/CT systems from two major vendors) based on a standardised protocol. Each site had to perform three measurements with the two sets of 133Ba sources and liquid 131I.

RESULTS

As anticipated, the 131I pseudo-image calibration factors (cps/MBq) were higher than those for 133Ba for all reconstructions and systems. A site-specific cross-calibration reduced the performance differences between both radionuclides with respect to a cross-calibration based on the ratio of emission probabilities from a median of 12-1.5%. The site-specific cross-calibration method also showed agreement between 133Ba and 131I for all cylinder volumes, which highlights the potential use of 133Ba sources to calculate recovery coefficients for partial volume correction.

CONCLUSION

This comparison exercise demonstrated that traceable solid 133Ba sources can be used as surrogate for liquid 131I imaging. The use of solid surrogate sources could solve the radiation protection problem inherent in the preparation of phantoms with 131I liquid activity solutions as well as reduce the measurement uncertainties in the activity. This is particularly relevant for stability measurements, which have to be carried out at regular intervals.

Averaging of absorbed doses: How matter matters

BACKGROUND

Dosimetry in radionuclide therapy often requires the calculation of average absorbed doses within and between spatial regions, for example, for voxel-based dosimetry methods, for paired organs, or across multiple tumors. Formation of such averages can be made in different ways, starting from different definitions.

PURPOSE

The aim of this study is to formally specify different averaging strategies for absorbed doses, and to compare their results when applied to absorbed dose distributions that are non-uniform within and between regions.

METHODS

For averaging within regions, two definitions of the average absorbed dose are considered: the simple average over the region (the region average) and the average when weighting by the mass density (density-weighted region average). The latter is shown to follow from the definition of mean absorbed dose according to the ICRU, and to be consistent with the MIRD formalism. For averaging between different spatial regions, three definitions follow: the volume-weighted, the mass-weighted, and the unweighted average. With respect to characterizing non-uniformity, the different average definitions lead to the use of dose-volume histograms (DVHs) (region average), dose-mass histograms (DMHs) (density-weighted region average), and unweighted histograms (unweighted average). Average absorbed doses are calculated for three worked examples, starting from the different definitions. The first, schematic, example concerns the calculation of the average absorbed dose between two regions with different volumes or mass densities. The second, stylized, example concerns voxel-based dosimetry, for which the average absorbed-dose rate within a region is calculated. The geometries studied include three 177 Lu-filled voxelized spheres, where the sphere masses are held constant while the material compositions, densities, and volumes are varied. For comparison, the mean absorbed-dose rates obtained using unit-density sphere S-values are also included. The third example concerns SPECT/CT-based tumor dosimetry for five patients undergoing therapy with 177 Lu-PSMA and six patients undergoing therapy with 177 Lu-DOTA-TATE, for which the average absorbed-dose rates across multiple tumors are calculated. For the second and third examples, analyses also include representations by histograms.

RESULTS

Example 1 shows that the average absorbed doses, calculated using different definitions, can differ considerably if the masses and absorbed doses for two regions are markedly different. From example 2 it is seen that the density-weighted region average is stable under different activity and density distributions and is also in line with results using S-values. In contrast, the region average varies as function of the activity distribution. In example 3, the absorbed dose rates for individual tumors differ by (1.1 ± 4.3)% and (-0.1 ± 0.4)% with maximum deviations of +34.4% and -1.4% for 177 Lu-PSMA and 177 Lu-DOTA-TATE, respectively, when calculated as region averages or density-weighted region averages, with largest deviations obtained when the density is non-uniform. The average absorbed doses calculated across all tumors are similar when comparing mass-weighted and volume-weighted averages but these differ substantially from unweighted averages.

CONCLUSION

Different strategies for averaging of absorbed doses within and between regions can lead to substantially different absorbed-dose estimates. At reporting of radionuclide therapy dosimetry, it is important to specify the averaging strategy applied.

Relationships between uptake of [68Ga]Ga-DOTA-TATE and absorbed dose in [177Lu]Lu-DOTA-TATE therapy

BACKGROUND

Somatostatin receptor ⁶⁸Ga PET imaging is standard for evaluation of a patient's suitability for ¹⁷⁷Lu peptide receptor radionuclide therapy of neuroendocrine tumours (NETs). The ⁶⁸Ga PET serves to ensure sufficient somatostatin receptor expression, commonly evaluated qualitatively. The aim of this study is to investigate the quantitative relationships between uptake in ⁶⁸Ga PET and absorbed doses in ¹⁷⁷Lu therapy.

METHOD

Eighteen patients underwent [⁶⁸Ga]Ga-DOTA-TATE PET imaging within 20 weeks prior to their first cycle of [¹⁷⁷Lu]Lu-DOTA-TATE. Absorbed doses for therapy were estimated for tumours, kidney, spleen, and normal liver parenchyma using a hybrid SPECT/CT-planar method. Gallium-68 activity concentrations were retrieved from PET images and also used to calculate SUVs and normalized SUVs, using blood and tissue for normalization. The ⁶⁸Ga activity concentrations per injected activity, SUVs, and normalized SUVs were compared with ¹⁷⁷Lu activity concentrations 1 d post-injection and ¹⁷⁷Lu absorbed doses. For tumours, for which there was a variable number per patient, both inter- and intra-patient correlations were analysed. Furthermore, the prediction of ¹⁷⁷Lu tumour absorbed doses based on a combination of tumour-specific ⁶⁸Ga activity concentrations and group-based estimates of the effective half-lives for grade 1 and 2 NETs was explored.

RESULTS

For normal organs, only spleen showed a significant correlation between the ⁶⁸Ga activity concentration and ¹⁷⁷Lu absorbed dose (r = 0.6). For tumours, significant, but moderate, correlations were obtained, with respect to both inter-patient (r = 0.7) and intra-patient (r = 0.45) analyses. The correlations to absorbed doses did not improve when using ⁶⁸Ga SUVs or normalized SUVs. The relationship between activity uptakes for ⁶⁸Ga PET and ¹⁷⁷Lu SPECT was stronger, with correlation coefficients r = 0.8 for both inter- and intra-patient analyses. The ¹⁷⁷Lu absorbed dose to tumour could be predicted from the ⁶⁸Ga activity concentrations with a 95% coverage interval of - 65% to 248%.

CONCLUSIONS

On a group level, a high uptake of [⁶⁸Ga]Ga-DOTA-TATE is associated with high absorbed doses at ¹⁷⁷Lu-DOTA-TATE therapy, but the relationship has a limited potential with respect to individual absorbed dose planning. Using SUV or SUV normalized to reference tissues do not improve correlations compared with using activity concentration per injected activity.

Monte Carlo modelling of a compact CZT-based gamma camera with application to 177Lu imaging

Background

Semiconductor gamma-camera systems based on cadmium zinc telluride (CZT) detectors present new challenges due to an energy-response that includes effects of low-energy tailing. In particular, such energy tails produce effects that need to be considered when imaging radionuclides with multiple emissions such as $$^{177}{\mathrm {Lu}}$$ 177 Lu . Monte Carlo simulation can be used to investigate the behaviour of such systems and optimise their use, provided that the detector model closely reflects the real physical detector. The aim of this work is to develop a CZT model applicable for simulation of CZT-based gamma cameras.

Methods

The equations describing the charge transport and signal induction are considered in three dimensions and are solved numerically, and the CZT model is then realised by coupling the detector-response to the photon-transport handled by the SIMIND Monte Carlo program. The CZT model is tuned to reproduce experimentally measured energy spectra of a hand-held gamma camera system for multiple radionuclides ($$^{99\mathrm {m}}{\mathrm {Tc}}$$ 99 m Tc , $$^{123}{\mathrm {I}}$$ 123 I and $$^{177}{\mathrm {Lu}}$$ 177 Lu ) and parallel-hole collimators (MEGP, LEHR) as well as an uncollimated system.

Results

Overall, the model results agree well with measurements across the range of experimental conditions. The applicability of the model is demonstrated by separating energy spectra into components to investigate the interference of high-energy photons on lower energy-windows, where pronounced effects of low-energy tailing for $$^{177}{\mathrm {Lu}}$$ 177 Lu are observed.

Conclusions

The developed model provides understanding of the specifics of the camera response and is expected to be helpful for future optimisation of gamma camera applications.

A multicentre and multi-national evaluation of the accuracy of quantitative Lu-177 SPECT/CT imaging performed within the MRTDosimetry project

PURPOSE

Patient-specific dosimetry is required to ensure the safety of molecular radiotherapy and to predict response. Dosimetry involves several steps, the first of which is the determination of the activity of the radiopharmaceutical taken up by an organ/lesion over time. As uncertainties propagate along each of the subsequent steps (integration of the time-activity curve, absorbed dose calculation), establishing a reliable activity quantification is essential. The MRTDosimetry project was a European initiative to bring together expertise in metrology and nuclear medicine research, with one main goal of standardizing quantitative 177Lu SPECT/CT imaging based on a calibration protocol developed and tested in a multicentre inter-comparison. This study presents the setup and results of this comparison exercise.

METHODS

The inter-comparison included nine SPECT/CT systems. Each site performed a set of three measurements with the same setup (system, acquisition and reconstruction): (1) Determination of an image calibration for conversion from counts to activity concentration (large cylinder phantom), (2) determination of recovery coefficients for partial volume correction (IEC NEMA PET body phantom with sphere inserts), (3) validation of the established quantitative imaging setup using a 3D printed two-organ phantom (ICRP110-based kidney and spleen). In contrast to previous efforts, traceability of the activity measurement was required for each participant, and all participants were asked to calculate uncertainties for their SPECT-based activities.

RESULTS

Similar combinations of imaging system and reconstruction lead to similar image calibration factors. The activity ratio results of the anthropomorphic phantom validation demonstrate significant harmonization of quantitative imaging performance between the sites with all sites falling within one standard deviation of the mean values for all inserts. Activity recovery was underestimated for total kidney, spleen, and kidney cortex, while it was overestimated for the medulla.

CONCLUSION

This international comparison exercise demonstrates that harmonization of quantitative SPECT/CT is feasible when following very specific instructions of a dedicated calibration protocol, as developed within the MRTDosimetry project. While quantitative imaging performance demonstrates significant harmonization, an over- and underestimation of the activity recovery highlights the limitations of any partial volume correction in the presence of spill-in and spill-out between two adjacent volumes of interests.

Dosimetric quantities of neuroendocrine tumors over treatment cycles with 177Lu-DOTA-TATE

Tumor dosimetry was performed for ¹⁷⁷Lu-DOTA-TATE with the aims of better understanding i) the range and variation of the tumor absorbed doses (ADs), ii) how different dosimetric quantities evolve over the treatment cycles, and iii) whether this evolution differs depending on the tumor grade. Such information is important for radiobiological interpretation and may inform the design of alternative administration schemes. Methods: Data come from 41 patients with neuroendocrine tumors (NETs) of grade 1 (n = 23) or 2 (n = 18), that had received between 2 and 9 treatment cycles. Dosimetry was performed for 182 individual lesions, giving in total 880 individual AD assessments across all cycles. Hybrid planar-SPECT/CT imaging was used, including quantitative SPECT reconstruction, voxel-based absorbed-dose-rate calculation, semi-automatic image segmentation, and partial-volume correction. Linear mixed-effect models were used to analyze changes over cycles in tumor ADs, absorbed-dose rates and activity concentrations at day-1, effective half-times, and tumor volumes. Tumors smaller than 8 ml were excluded from analyses. Results: Tumor ADs ranged between 2 and 77 Gy per cycle. On average the AD decreased over the cycles, with significantly different rates (P < 0.05) for grade 1 and 2 NETs of 6% and 14% per cycle, respectively. The absorbed-dose rates and activity concentrations at day-1 decreased by similar amounts. The effective half-times were less variable but shorter for grade 2 than grade 1 (P < 0.001). For grade 2 NETS the tumor volumes decreased, with a similar tendency in grade 1. Conclusion: The tumor AD, absorbed-dose rate and activity uptake decrease, in parallel with tumor volumes, between ¹⁷⁷Lu-DOTA-TATE treatment cycles, particularly for grade 2 NETs. The effective half-times vary less but are lower for grade 2 than grade 1 NETs. These results may indicate the development of radiation-induced fibrosis and could have implications for the design of future treatment and dosimetry protocols.

Characterisation of a hand-held CZT-based gamma camera for 177Lu imaging

BACKGROUND

Currently, hand-held gamma cameras are being developed for ^99mTc imaging, mainly for sentinel lymph node detection. These cameras offer advantages, such as mobility and ease of access, and may be useful also for other applications such as biokinetic studies in animals or for imaging of small, superficial structures in patients. In this work, the suitability of a CZT-based hand-held camera for ¹⁷⁷Lu imaging is investigated. The energy response of CZT-based detectors combined with the multiple photon emissions of ¹⁷⁷Lu poses new challenges compared to ^99mTc imaging, and a thorough camera characterisation is thus warranted.

METHODS

Three collimators (LEHR, LEHS, and MEGP) and three energy windows (55 keV, 113 keV, and 208 keV) are investigated. Characterised camera properties include the system spatial resolution, energy resolution, sensitivity, image uniformity, septal penetration, and temperature dependence. Characterisations are made starting from NEMA guidelines when applicable, with adjustments made when required. The applicability of the camera is demonstrated by imaging of a superficially located tumour in a patient undergoing [¹⁷⁷ Lu]Lu-DOTA-TATE therapy.

RESULTS

Overall, the results are encouraging. Compared to a conventional gamma camera, the hand-held camera generally has a higher sensitivity for a given collimator. For source-collimator distances below 3 cm, the spatial resolution FWHM is within 6 mm for the LEHR and MEGP collimators. Before uniformity correction, the central field-of-view integral uniformity shows best results for the 113-keV window, with values obtained between 11 and 14%. The corresponding values after uniformity correction are within 3%. Effects of septal penetration are observed but are manageable with a proper combination of collimator and energy window setting. Septal penetration and collimator scatter not only affect the 208-keV window but also contribute with counts in lower windows due to energy-tailing effects. The patient study revealed non-uniform uptake patterns in a region that appeared uniform in a conventional gamma camera image.

CONCLUSIONS

The results show that the hand-held camera can be used for ¹⁷⁷Lu imaging. A 113-keV energy window combined with LEHR or MEGP collimators provides the best image system characteristics.

Bone Marrow Absorbed Doses and Correlations with Hematologic Response During 177Lu-DOTATATE Treatments Are Influenced by Image-Based Dosimetry Method and Presence of Skeletal Metastases

This study aimed to compare different image-based methods for bone marrow dosimetry and study the dose–response relationship during treatment with ¹⁷⁷Lu-DOTATATE in patients with and without skeletal metastases. Methods: This study included 46 patients with advanced neuroendocrine tumors treated with at least 2 fractions of ¹⁷⁷Lu-DOTATATE at Sahlgrenska University Hospital. High- and low-uptake compartments were automatically outlined in planar images collected at 2, 24, 48, and 168 h after injection. The bone marrow absorbed doses were calculated from the cross doses of the high- and low-uptake compartments and the self-dose, using the time–activity concentration curve for the low-uptake compartment. This time–activity concentration curve was adjusted using a fixed constant of 1.8 for the planar dosimetry method and using the activity concentrations in vertebral bodies in SPECT images at 24 h after injection of ¹⁷⁷Lu-DOTATATE in 4 hybrid methods: L4-SPECT used the activity concentration in the L4 vertebra, whereas V-SPECT, L-SPECT, and T-SPECT used the median activity concentration in all visible vertebrae, lumbar vertebrae, and thoracic vertebrae, respectively. Results: Using the planar method, L4-SPECT, V-SPECT, L-SPECT, and T-SPECT, the estimated median bone marrow absorbed doses were 0.19, 0.36, 0.40, 0.39, and 0.46 Gy/7.4 GBq, respectively, with respective ranges of 0.12–0.33, 0.15–1.44, 0.19–1.71, 0.21–1.60, and 0.18–2.12 Gy/7.4 GBq. For all methods, the bone marrow absorbed dose significantly correlated with decreased platelet counts. This correlation increased after treatment fraction 2: the Spearman correlation (r_s) were −0.49 for the planar method, −0.61 for L4-SPECT, −0.63 for V-SPECT, −0.63 for L-SPECT, and −0.57 for T-SPECT. A separate analysis revealed an increased correlation for patients without skeletal metastases using the planar method (r_s = −0.67). In contrast, hybrid methods had poor correlations for patients without metastases and stronger correlations for patients with skeletal metastases (r_s = −0.61 to −0.74). The mean bone marrow absorbed doses were 3%–69% higher for patients with skeletal metastases than for patients without. Conclusion: The estimated bone marrow absorbed doses by image-based techniques and the correlation with platelets are influenced by the choice of measured vertebrae and the presence of skeletal metastases.

EANM practical guidance on uncertainty analysis for molecular radiotherapy absorbed dose calculations

A framework is proposed for modelling the uncertainty in the measurement processes constituting the dosimetry chain that are involved in internal absorbed dose calculations. The starting point is the basic model for absorbed dose in a site of interest as the product of the cumulated activity and a dose factor. In turn, the cumulated activity is given by the area under a time-activity curve derived from a time sequence of activity values. Each activity value is obtained in terms of a count rate, a calibration factor and a recovery coefficient (a correction for partial volume effects). The method to determine the recovery coefficient and the dose factor, both of which are dependent on the size of the volume of interest (VOI), are described. Consideration is given to propagating estimates of the quantities concerned and their associated uncertainties through the dosimetry chain to obtain an estimate of mean absorbed dose in the VOI and its associated uncertainty. This approach is demonstrated in a clinical example.

A method for tumor dosimetry based on hybrid planar-SPECT/CT images and semiautomatic segmentation

PURPOSE

A hybrid planar-SPECT/CT method for tumor dosimetry in 177 Lu-DOTATATE therapy, applicable to datasets consisting of multiple conjugate-view images and one SPECT/CT, is developed and evaluated.

METHODS

The imaging protocol includes conjugate-view imaging at 1, 24, 96, and 168 h post infusion (p.i.) and a SPECT/CT acquisition 24 h p.i. The dosimetry method uses the planar images to estimate the shape of the time-activity concentration curve, which is then rescaled to absolute units using the SPECT-derived activity concentration. The resulting time-integrated activity concentration coefficient (TIACC) is used to calculate the tumor-absorbed dose. Semiautomatic segmentation techniques are applied for tumor delineation in both planar and SPECT images, where the planar image segmentation is accomplished using an active-rays-based technique. The selection of tumors is done by visual inspection of planar and SPECT images and applying a set of criteria concerning the tumor visibility and possible interference from superimposed activity uptakes in the planar images. Five different strategies for determining values from planar regions of interest (ROIs), based on entire or partial ROIs, and with and without background correction, are evaluated. Evaluation is performed against a SPECT/CT-based method on data from six patients where sequential conjugate-view and SPECT/CT imaging have been performed in parallel and against ground truths in Monte Carlo simulated images. The patient data are also used to evaluate the interoperator variability and to assess the validity of the developed criteria for tumor selection.

RESULTS

For patient images, the hybrid method produces TIACCs that are on average 6% below those of the SPECT/CT only method, with standard deviations for the relative TIACC differences of 8%-11%. Simulations show that the hybrid and SPECT-based methods estimate the TIACCs to within approximately 10% for tumors larger than around 10 ml, while for smaller tumors, all methods underestimate the TIACCs due to underestimations of the activity concentrations in the SPECT images. The planar image segmentation has a low operator dependence, with a median Dice similarity coefficient of 0.97 between operators. The adopted criteria for tumor selection manage to discriminate the tumors for which the absorbed-dose deviations between the hybrid and SPECT methods are the highest.

CONCLUSIONS

The hybrid method is found suitable for studies of tumor-absorbed doses in radionuclide therapy, provided that selection criteria regarding the visibility and overlapping activities in the planar images are applied.

Feasibility of simplifying renal dosimetry in 177Lu peptide receptor radionuclide therapy

BACKGROUND

Recently, 177Lu-dotatate therapy for neuroendocrine tumours has received regulatory approval. Dosimetry can be used to optimize treatment on an individual basis, but there is no international consensus as to how it should be done. The aim of this study is to determine a feasible and accurate dosimetry method to guide individualized peptide receptor radionuclide therapy (PRRT) for patients with neuroendocrine tumours. As part of a clinical trial on 177Lu-dotatate therapy, renal dosimetry was performed for all patients in each treatment cycle, using a hybrid planar-SPECT/CT method. In the present study, we use the image data acquired from 22 patients and 119 cycles and define a set of alternative treatment planning strategies, each representing a simplification in terms of image acquisition and dosimetric calculations. The results from the simplified strategies are compared to the results from the protocol-prescribed hybrid planar-SPECT/CT-based method by analysing differences both in per-cycle and total cumulative absorbed dose (AD) analyses.

RESULTS

In general, the SPECT-based methods gave results that were largely consistent with the protocol-specified hybrid method, both in the per-cycle and cumulative AD analyses. Notably, performing one SPECT/CT per cycle at 96 h yielded ADs that were very similar to the protocol method. The methods using planar dosimetry resulted in larger variations, as expected, while giving 4 cycles to all patients resulted in the largest inter-individual differences in cumulative AD.

CONCLUSIONS

Performing one SPECT/CT at 96 h in every treatment cycle gives sufficiently reliable dosimetric results to base individualized treatment planning on, with a reasonable demand on resources.

Variations in the practice of molecular radiotherapy and implementation of dosimetry: results from a European survey

BACKGROUND

Currently, the implementation of dosimetry in molecular radiotherapy (MRT) is not well investigated, and in view of the Council Directive (2013/59/Euratom), there is a need to understand the current availability of dosimetry-based MRT in clinical practice and research studies. The aim of this study was to assess the current practice of MRT and dosimetry across European countries.

METHODS

An electronic questionnaire was distributed to European countries. This addressed 18 explicitly considered therapies, and for each therapy, a similar set of questions were included. Questions covered the number of patients and treatments during 2015, involvement of medical specialties and medical physicists, implementation of absorbed dose planning, post-therapy imaging and dosimetry, and the basis of therapy prescription.

RESULTS

Responses were obtained from 26 countries and 208 hospitals, administering in total 42,853 treatments. The most common therapies were 131I-NaI for benign thyroid diseases and thyroid ablation of adults. The involvement of a medical physicist (mean over all 18 therapies) was reported to be either minority or never by 32% of the responders. The percentage of responders that reported that dosimetry was included on an always/majority basis differed between the therapies and showed a median value of 36%. The highest percentages were obtained for 177Lu-PSMA therapy (100%), 90Y microspheres of glass (84%) and resin (82%), 131I-mIBG for neuroblastoma (59%), and 131I-NaI for benign thyroid diseases (54%). The majority of therapies were prescribed based on fixed-activity protocols. The highest number of absorbed-dose based prescriptions were reported for 90Y microsphere treatments in the liver (64% and 96% of responses for resin and glass, respectively), 131I-NaI treatment of benign thyroid diseases (38% of responses), and for 131I-mIBG treatment of neuroblastoma (18% of responses).

CONCLUSIONS

There is a wide variation in MRT practice across Europe and for different therapies, including the extent of medical-physicist involvement and the implementation of dosimetry-guided treatments.

Dosimetry-based treatment planning for molecular radiotherapy: a summary of the 2017 report from the Internal Dosimetry Task Force

BACKGROUND

The European directive on basic safety standards (Council directive 2013/59 Euratom) mandates dosimetry-based treatment planning for radiopharmaceutical therapies. The directive comes into operation February 2018, and the aim of a report produced by the Internal Dosimetry Task Force of the European Association of Nuclear Medicine is to address this aspect of the directive. A summary of the report is presented.

RESULTS

A brief review of five of the most common therapy procedures is included in the current text, focused on the potential to perform patient-specific dosimetry. In the full report, 11 different therapeutic procedures are included, allowing additional considerations of effectiveness, references to specific literature on quantitative imaging and dosimetry, and existing evidence for absorbed dose-effect correlations for each treatment. Individualized treatment planning with tracer diagnostics and verification of the absorbed doses delivered following therapy is found to be scientifically feasible for almost all procedures investigated, using quantitative imaging and/or external monitoring. Translation of this directive into clinical practice will have significant implications for resource requirements.

CONCLUSIONS

Molecular radiotherapy is undergoing a significant expansion, and the groundwork for dosimetry-based treatment planning is already in place. The mandated individualization is likely to improve the effectiveness of the treatments, although must be adequately resourced.

SPECT image segmentation for estimation of tumour volume and activity concentration in 177Lu-DOTATATE radionuclide therapy

Background

Dosimetry in radionuclide therapy has the potential to allow for a treatment tailored to the individual patient. One therapeutic radiopharmaceutical where patient-specific dosimetry is feasible is ¹⁷⁷Lu-DOTATATE, used for the treatment of neuroendocrine tumours. The emission of gamma photons by ¹⁷⁷Lu allows for imaging with SPECT (single photon emission computed tomography). One important step for dosimetry using this imaging technique is the SPECT image segmentation, which needs to be robust and accurate for the estimated quantities to be reliable. This work investigates different methods for automatic tumour delineation in ¹⁷⁷Lu-DOTATATE SPECT images.

Three segmentation methods are considered: a fixed 42% threshold (FT), the Otsu method (OM) and a method based on Fourier surfaces (FS). Effects of including resolution compensation in the iterative SPECT image reconstruction are also studied. Evaluation is performed based on Monte Carlo-simulated SPECT images from 24 h and 336 h post injection (p.i.), for determination of the volume, activity concentration and dice similarity coefficient. In addition, patient data are used to investigate the correspondence of tumour volumes when delineated in SPECT or morphological CT or MR images. Patient data are also used to examine the sensitivity to the operator-dependent initialization.

Results

For simulated images from 24 h p.i. reconstructed without resolution compensation, a volume and activity-concentration root-mean-square error below 15% is typically obtained for tumours above approximately 10 cm³ when using OM or FS, while FT performs considerably worse. When including resolution compensation, the tumour volume becomes underestimated and the activity concentration overestimated. The FS method appears to be robust to noise, as seen for the 336 h images. The differences between the tumour volumes estimated from the SPECT images and the volumes estimated from morphological images are generally larger than the discrepancies seen for the simulated data sets.

Conclusions

Segmentation results are encouraging for future dosimetry of tumours with volumes above approximately 10 cm³. Using resolution compensation in the reconstruction may have a negative effect on volume estimation.

Whole-remnant and maximum-voxel SPECT/CT dosimetry in 131 I-NaI treatments of differentiated thyroid cancer

PURPOSE

To investigate the possible differences between SPECT/CT based whole-remnant and maximum-voxel dosimetry in patients receiving radio-iodine ablation treatment of differentiated thyroid cancer (DTC).

METHODS

Eighteen DTC patients were administered 1.11 GBq of 131 I-NaI after near-total thyroidectomy and rhTSH stimulation. Two patients had two remnants, so in total dosimetry was performed for 20 sites. Three SPECT/CT scans were performed for each patient at 1, 2, and 3-7 days after administration. The activity, the remnant mass, and the maximum-voxel activity were determined from these images and from a recovery-coefficient curve derived from experimental phantom measurements. The cumulated activity was estimated using trapezoidal-exponential integration. Finally, the absorbed dose was calculated using S-values for unit-density spheres in whole-remnant dosimetry and S-values for voxels in maximum-voxel dosimetry.

RESULTS

The mean absorbed dose obtained from whole-remnant dosimetry was 40 Gy (range 2-176 Gy) and from maximum-voxel dosimetry 34 Gy (range 2-145 Gy). For any given patient, the activity concentrations for each of the three time-points were approximately the same for the two methods. The effective half-lives varied (R = 0.865), mainly due to discrepancies in estimation of the longer effective half-lives. On average, absorbed doses obtained from whole-remnant dosimetry were 1.2 ± 0.2 (1 SD) higher than for maximum-voxel dosimetry, mainly due to differences in theS-values. The method-related differences were however small in comparison to the wide range of absorbed doses obtained in patients.

CONCLUSIONS

Simple and consistent procedures for SPECT/CT based whole-volume and maximum-voxel dosimetry have been described, both based on experimentally determined recovery coefficients. Generally the results from the two approaches are consistent, although there is a small, systematic difference in the absorbed dose due to differences in the S-values, and some variability due to differences in the estimated effective half-lives, especially when the effective half-life is long. Irrespective of the method used, the patient absorbed doses obtained span over two orders of magnitude.

Biologically effective dose in fractionated molecular radiotherapy--application to treatment of neuroblastoma with (131)I-mIBG

In this work, the biologically effective dose (BED) is investigated for fractionated molecular radiotherapy (MRT). A formula for the Lea-Catcheside G-factor is derived which takes the possibility of combinations of sub-lethal damage due to radiation from different administrations of activity into account. In contrast to the previous formula, the new G-factor has an explicit dependence on the time interval between administrations. The BED of tumour and liver is analysed in MRT of neuroblastoma with (131)I-mIBG, following a common two-administration protocol with a mass-based activity prescription. A BED analysis is also made for modified schedules, when due to local regulations there is a maximum permitted activity for each administration. Modifications include both the simplistic approach of delivering this maximum permitted activity in each of the two administrations, and also the introduction of additional administrations while maintaining the protocol-prescribed total activity. For the cases studied with additional (i.e. more than two) administrations, BED of tumour and liver decreases at most 12% and 29%, respectively. The decrease in BED of the tumour is however modest compared to the two-administration schedule using the maximum permitted activity, where the decrease compared to the original schedule is 47%.

Pharmacokinetic digital phantoms for accuracy assessment of image-based dosimetry in (177)Lu-DOTATATE peptide receptor radionuclide therapy

Patient-specific image-based dosimetry is considered to be a useful tool to limit toxicity associated with peptide receptor radionuclide therapy (PRRT). To facilitate the establishment and reliability of absorbed-dose response relationships, it is essential to assess the accuracy of dosimetry in clinically realistic scenarios. To this end, we developed pharmacokinetic digital phantoms corresponding to patients treated with (177)Lu-DOTATATE. Three individual voxel phantoms from the XCAT population were generated and assigned a dynamic activity distribution based on a compartment model for (177)Lu-DOTATATE, designed specifically for this purpose. The compartment model was fitted to time-activity data from 10 patients, primarily acquired using quantitative scintillation camera imaging. S values for all phantom source-target combinations were calculated based on Monte-Carlo simulations. Combining the S values and time-activity curves, reference values of the absorbed dose to the phantom kidneys, liver, spleen, tumours and whole-body were calculated. The phantoms were used in a virtual dosimetry study, using Monte-Carlo simulated gamma-camera images and conventional methods for absorbed-dose calculations. The characteristics of the SPECT and WB planar images were found to well represent those of real patient images, capturing the difficulties present in image-based dosimetry. The phantoms are expected to be useful for further studies and optimisation of clinical dosimetry in (177)Lu PRRT.

Hybrid Imaging for Patient-Specific Dosimetry in Radionuclide Therapy

Radionuclide therapy aims to treat malignant diseases by systemic administration of radiopharmaceuticals, often using carrier molecules such as peptides and antibodies. The radionuclides used emit electrons or alpha particles as a consequence of radioactive decay, thus leading to local energy deposition. Administration to individual patients can be tailored with regards to the risk of toxicity in normal organs by using absorbed dose planning. The scintillation camera, employed in planar imaging or single-photon emission computed tomography (SPECT), generates images of the spatially and temporally varying activity distribution. Recent commercially available combined SPECT and computed tomography (CT) systems have dramatically increased the possibility of performing accurate dose planning by using the CT information in several steps of the dose-planning calculation chain. This paper discusses the dosimetry chain used for individual absorbed-dose planning and highlights the areas where hybrid imaging makes significant contributions.

Dosimetric results in treatments of neuroblastoma and neuroendocrine tumors with (131)I-metaiodobenzylguanidine with implications for the activity to administer

PURPOSE

The aim was to investigate whole-body and red marrow absorbed doses in treatments of neuroblastoma (NB) and adult neuroendocrine tumors (NETs) with (131)I-metaiodobenzylguanidine and to propose a simple method for determining the activity to administer when dosimetric data for the individual patient are not available.

METHODS

Nine NB patients and six NET patients were included, giving in total 19 treatments as four patients were treated twice. Whole-body absorbed doses were determined from dose-rate measurements and planar gamma-camera imaging. For six NB and five NET treatments, red marrow absorbed doses were also determined using the blood-based method.

RESULTS

Dosimetric data from repeated administrations in the same patient were consistent. In groups of NB and NET patients, similar whole-body residence times were obtained, implying that whole-body absorbed dose per unit of administered activity could be reasonably well described as a power function of the patient mass. For NB, this functional form was found to be consistent with dosimetric data from previously published studies. The whole-body to red marrow absorbed dose ratio was similar among patients, with values of 1.4 ± 0.6-1.7 ± 0.7 (1 standard deviation) in NB treatments and between 1.5 ± 0.6 and 1.7 ± 0.7 (1 standard deviation) in NET treatments.

CONCLUSIONS

The consistency of dosimetric results between administrations for the same patient supports prescription of the activity based on dosimetry performed in pretreatment studies, or during the first administration in a fractionated schedule. The expressions obtained for whole-body absorbed doses per unit of administered activity as a function of patient mass for NB and NET treatments are believed to be a useful tool to estimate the activity to administer at the stage when the individual patient biokinetics has not yet been measured.

The evidence base for the use of internal dosimetry in the clinical practice of molecular radiotherapy

Molecular radiotherapy (MRT) has demonstrated unique therapeutic advantages in the treatment of an increasing number of cancers. As with other treatment modalities, there is related toxicity to a number of organs at risk. Despite the large number of clinical trials over the past several decades, considerable uncertainties still remain regarding the optimization of this therapeutic approach and one of the vital issues to be answered is whether an absorbed radiation dose-response exists that could be used to guide personalized treatment. There are only limited and sporadic data investigating MRT dosimetry. The determination of dose-effect relationships for MRT has yet to be the explicit aim of a clinical trial. The aim of this article was to collate and discuss the available evidence for an absorbed radiation dose-effect relationships in MRT through a review of published data. Based on a PubMed search, 92 papers were found. Out of 79 studies investigating dosimetry, an absorbed dose-effect correlation was found in 48. The application of radiobiological modelling to clinical data is of increasing importance and the limited published data on absorbed dose-effect relationships based on these models are also reviewed. Based on National Cancer Institute guideline definition, the studies had a moderate or low rate of clinical relevance due to the limited number of studies investigating overall survival and absorbed dose. Nevertheless, the evidence strongly implies a correlation between the absorbed doses delivered and the response and toxicity, indicating that dosimetry-based personalized treatments would improve outcome and increase survival.

Patient-specific whole-body attenuation correction maps from a CT system for conjugate-view-based activity quantification: method development and evaluation

For activity quantification based on planar scintillation camera measurements, photon attenuation is an important factor that needs to be corrected for in a patient- and organ-specific manner. One possibility for obtaining attenuation correction maps is to use X-ray CT scout images. Since the intensity of scout images is in relative numbers, their image values need to be multiplied by a factor to become quantitative and thus useful for attenuation correction. The calibration factor can for our current imaging system be obtained from a scanner system file, but is generally not available. For this purpose, a method based on the patient weight has been developed. Results based on 79 patient scout images show that the calibration factor thus determined correlates well with values that, in this case, are independently specified by the system. The accuracy of attenuation correction factors (ACFs) derived from the scout-based attenuation correction maps is evaluated by comparison to ACFs derived from three-dimensional CT studies. For photon energies of 208, 245, and 364 keV, scout-based ACFs are on average 1.2% and 0.5% from the CT-derived values, using the system-based and the weight-based values of the scout-image calibration factor, respectively. The imprecision is somewhat higher for the weight-based method, due to variability in the delineation of the patient contour used as a part of this method. In conclusion, X-ray scouts are found useful for attenuation correction with a satisfactory accuracy obtained, both using the new, weight-based method, and using the previous, system-based method, for determining the required calibration factor.

Dynamic (99m)Tc-MAG3 renography: images for quality control obtained by combining pharmacokinetic modelling, an anthropomorphic computer phantom and Monte Carlo simulated scintillation camera imaging

In dynamic renal scintigraphy, the main interest is the radiopharmaceutical redistribution as a function of time. Quality control (QC) of renal procedures often relies on phantom experiments to compare image-based results with the measurement setup. A phantom with a realistic anatomy and time-varying activity distribution is therefore desirable. This work describes a pharmacokinetic (PK) compartment model for (99m)Tc-MAG3, used for defining a dynamic whole-body activity distribution within a digital phantom (XCAT) for accurate Monte Carlo (MC)-based images for QC. Each phantom structure is assigned a time-activity curve provided by the PK model, employing parameter values consistent with MAG3 pharmacokinetics. This approach ensures that the total amount of tracer in the phantom is preserved between time points, and it allows for modifications of the pharmacokinetics in a controlled fashion. By adjusting parameter values in the PK model, different clinically realistic scenarios can be mimicked, regarding, e.g., the relative renal uptake and renal transit time. Using the MC code SIMIND, a complete set of renography images including effects of photon attenuation, scattering, limited spatial resolution and noise, are simulated. The obtained image data can be used to evaluate quantitative techniques and computer software in clinical renography.

On the biologically effective dose (BED)-using convolution for calculating the effects of repair: II. Numerical considerations

We have previously shown analytically that the biologically effective dose (BED), including effects of repair during irradiation and of incomplete repair between fractions, can be formulated using a convolution between the absorbed dose rate function and the function describing repair. In this work, a discrete formalism is derived along with its implementation via the fast Fourier transform. The implementation takes the intrinsic periodicity of the discrete Fourier transform into consideration, as well as possible inconsistencies that may arise due to discretization and truncation of the functions describing the absorbed dose rate and repair. Numerically and analytically calculated BED values are compared for various situations in external beam radiotherapy, brachytherapy and radionuclide therapy, including the use of different repair models. The numerical method is shown to be accurate and versatile since it can be applied to any kind of absorbed dose rate function and allows for the incorporation of different repair models. Typical accuracies for clinically realistic examples are in the order of 10(-3)% to 10(-5)%. The method has thus the potential of being a useful tool for the calculation of BED, also in situations with complicated irradiation patterns or repair functions.

High-dose iodine-131-metaiodobenzylguanidine with haploidentical stem cell transplantation and posttransplant immunotherapy in children with relapsed/refractory neuroblastoma

We evaluated the feasibility and efficacy of using high-dose iodine-131-metaiodobenzylguanidine ((131)I-MIBG) followed by reduced-intensity conditioning (RIC) and transplantation of T cell-depleted haploidentical peripheral blood stem cells (designated haplo-SCT) to treat relapsing/refractory neuroblastoma (RRNB). Five RRNB patients were enrolled: 4 with relapse (3 after autologous SCT) and 1 with induction therapy failure. The preparative regimen included high-dose (131)I-MIBG on day -20, followed by fludarabine (Flu), thiotepa, and melphalan (Mel) from day -8 to -1. Granulocyte-colony stimulating factor (G-CSF)-mobilized, T cell-depleted haploidentical paternal stem cells were infused on day 0 together with cultured donor mesenchymal stem cells. A single dose of rituximab was given on day +1. After cessation of short immunosuppression (mycophenolate, OKT3), 4 children received donor lymphocyte infusion (DLI). (131)I-MIBG infusion and RIC were well tolerated. All patients engrafted. No primary acute graft-versus-host disease (aGVHD) was observed. Four children developed aGVHD after DLI and were successfully treated. Analysis of immunologic recovery showed fast reappearance of potentially immunocompetent natural killer (NK) and T cells, which might have acted as effector cells responsible for the graft-versus-tumor (GVT) effect. Two children are alive and well, with no evidence of disease 40 and 42 months after transplantation. One patient experienced late progression with new bone lesions (sternum) 38 months after haplo-SCT, and is being treated with local irradiation and reinstituted DLI. One patient rejected the graft, was rescued with autologous backup, and died of progressive disease 5 months after transplantation. Another child relapsed 7 months after transplantation and died 5 months later. High-dose (131)I-MIBG followed by RIC and haplo-SCT for RRNB is feasible and promising, because 2 of 5 children on that regimen achieved long-lasting remission. Further studies are needed to evaluate targeted therapy and immune-mediated tumor control in high-risk neuroblastoma.

Consequences of inadvertent radioiodine treatment of Graves' disease and thyroid cancer in undiagnosed pregnancy. Can we rely on routine pregnancy testing?

INTRODUCTION

Radioiodine and most cytostatic treatments are contraindicated in pregnancy. Still, inadvertent therapy does occur. Radioiodine was given to two pregnant women with Graves' disease and thyroid cancer respectively, both in their 20th gestational week. Routine pregnancy tests based on urinary beta-hCG had failed to indicate pregnancy in both cases.

METHODS

Estimation of doses to the foetuses and foetal thyroids. Scrutiny of pregnancy testing.

RESULTS AND CONCLUSIONS

Doses to foetal thyroids were ablative (250-600 Gy). Total foetal dose in the Graves' patient was 100 mGy and compatible with survival, whereas a foetal dose of approximately 700 mGy together with induced hypothyroidism was fatal for the foetus of the cancer patient. Routine pregnancy tests may fail early and late in pregnancy. The possibility of pregnancy should be considered in all fertile women before therapy with radionuclides or cytostatic regimens, and a clinical investigation undertaken on wide indications with determination of serum beta-hCG, preferably together with an ultrasound examination.

Parametric images of antibody pharmacokinetics based on serial quantitative whole-body imaging and blood sampling

We present a method for pharmacokinetic modeling of distributions of (111)In-labeled monoclonal antibodies (mAbs) on individual pixels of planar scintillation-camera images.

METHODS

The method is applied to 2 sets of clinical whole-body images, each consisting of 6 consecutive images acquired over a week. Quantification is performed on a pixel basis, yielding images in units of Bq/pixel. The images acquired on the different occasions are registered using a nonrigid method, and for each pixel location a time-activity curve is obtained for which kinetic modeling is performed. The (111)In-mAb is assumed to be located in either the vascular or the extravascular space. The vascular content is assumed to follow the global blood kinetics as determined from blood samples, together with a model parameter alpha that describes the fraction of the whole-body blood volume present in the particular pixel. The rate of change of the extravascular compartment is described by a linear 1-tissue-compartment model with 2 rate constants, K'1 and k2, reflecting extravasation and washout, respectively. The model is optimized for each pixel position with regard to the values of the 3 parameters (alpha, K'1, and k2), resulting in 3 parametric images. From these, images of the cumulated activity in vascular and extravascular spaces are calculated, as is an image of the rate-constants ratio, which is closely related to the volume of distribution.

RESULTS

The resulting parametric images are analyzed in terms of the appearance of the time-activity curves at various locations. Results also include interpretation of the parametric images in their clinical context, and the location of regions that exhibit high extravasation and a low washout rate is compared with confirmed malignant sites.

CONCLUSION

Parametric imaging allows the study and analysis of the spatial and temporal distributions of mAbs simultaneously. Parametric imaging enhances regions where the pharmacokinetics differ from the surrounding tissue and provides a tool to detect and locate unexpected kinetic behavior, which is sometimes characteristic of malignant tissue. For dosimetry in radionuclide therapy, parametric imaging offers a less biased means of analyzing serial mAb images than traditional region-of-interest-based analysis.