Quantitative accuracy of 177Lu SPECT imaging for molecular radiotherapy

Predicting the spatio-temporal response of recurrent glioblastoma treated with rhenium-186 labelled nanoliposomes

Rhenium-186 (186Re) labeled nanoliposome (RNL) therapy for recurrent glioblastoma patients has shown promise to improve outcomes by locally delivering radiation to affected areas. To optimize the delivery of RNL, we have developed a framework to predict patient-specific response to RNL using image-guided mathematical models.

Methods

We calibrated a family of reaction-diffusion type models with multi-modality imaging data from ten patients (NCR01906385) to predict the spatio-temporal dynamics of each patient's tumor. The data consisted of longitudinal magnetic resonance imaging (MRI) and single photon emission computed tomography (SPECT) to estimate tumor burden and local RNL activity, respectively. The optimal model from the family was selected and used to predict future growth. A simplified version of the model was used in a leave-one-out analysis to predict the development of an individual patient's tumor, based on cohort parameters.

Results

Across the cohort, predictions using patient-specific parameters with the selected model were able to achieve Spearman correlation coefficients (SCC) of 0.98 and 0.93 for tumor volume and total cell number, respectively, when compared to the measured data. Predictions utilizing the leave-one-out method achieved SCCs of 0.89 and 0.88 for volume and total cell number across the population, respectively.

Conclusion

We have shown that patient-specific calibrations of a biology-based mathematical model can be used to make early predictions of response to RNL therapy. Furthermore, the leave-one-out framework indicates that radiation doses determined by SPECT can be used to assign model parameters to make predictions directly following the conclusion of RNL treatment.

Statement of Significance

This manuscript explores the application of computational models to predict response to radionuclide therapy in glioblastoma. There are few, to our knowledge, examples of mathematical models used in clinical radionuclide therapy. We have tested a family of models to determine the applicability of different radiation coupling terms for response to the localized radiation delivery. We show that with patient-specific parameter estimation, we can make accurate predictions of future glioblastoma response to the treatment. As a comparison, we have shown that population trends in response can be used to forecast growth from the moment the treatment has been delivered.In addition to the high simulation and prediction accuracy our modeling methods have achieved, the evaluation of a family of models has given insight into the response dynamics of radionuclide therapy. These dynamics, while different than we had initially hypothesized, should encourage future imaging studies involving high dosage radiation treatments, with specific emphasis on the local immune and vascular response.

Quantitative 177Lu SPECT/CT imaging for personalized dosimetry using a ring-shaped CZT-based camera

BACKGROUND

Dosimetry after radiopharmaceutical therapy with 177Lu (177Lu-RPT) relies on quantitative SPECT/CT imaging, for which suitable reconstruction protocols are required. In this study, we characterized for the first time the quantitative performance of a ring-shaped CZT-based camera using two different reconstruction algorithms: an ordered subset expectation maximization (OSEM) and a block sequential regularized expectation maximization (BSREM) combined with noise reduction regularization. This study lays the foundations for the definition of a reconstruction protocol enabling accurate dosimetry for patients treated with 177Lu-RPT.

METHODS

A series of 177Lu-filled phantoms were acquired on a StarGuide™ (GE HealthCare), with energy and scatter windows centred at 208 (± 6%) keV and 185 (± 5%) keV, respectively. Images were reconstructed with the manufacturer implementations of OSEM (GE-OSEM) and BSREM (Q.Clear) algorithms, and various combinations of iterations and subsets. Additionally, the manufacturer-recommended Q.Clear-based reconstruction protocol was evaluated. Quantification accuracy, measured as the difference between the SPECT-based and the radionuclide calibrator-based activity, and noise were evaluated in a large cylinder. Recovery coefficients (RCs) and spatial resolution were assessed in a NEMA IEC phantom with sphere inserts. The reconstruction protocols considered suitable for clinical applications were tested on a cohort of patients treated with [177Lu]Lu-PSMA-I&T.

RESULTS

The accuracy of the activity from the cylinder, although affected by septal penetration, was < 10% for all reconstructions. Both algorithms featured improved spatial resolution and higher RCs with increasing updates at the cost of noise build-up, but Q.Clear outperformed GE-OSEM in reducing noise accumulation. When the reconstruction parameters were carefully selected, similar values for noise (~0.15), spatial resolution (~1 cm) and RCs were found, irrespective of the reconstruction algorithm. Analogue results were found in patients.

CONCLUSIONS

Accurate activity quantification is possible when imaging 177Lu with StarGuide™. However, the impact of septal penetration requires further investigations. GE-OSEM is a valid alternative to the recommended Q.Clear reconstruction algorithm, featuring comparable performances assessed on phantoms and patients.

Reconstructed SPECT images of 177Lu homogeneous cylindrical phantom used for calibration and texture analysis

In a clinical contest, it is common to use dedicated phantoms to perform quality assurance test to check the performance of a SPECT system. Some of these phantoms are also used to calibrate the system for dosimetric evaluation of patients undergoing radiometabolic cancer therapy. In this work, a 3D-OSEM reconstructed ¹⁷⁷Lu SPECT dataset of a homogeneous cylindrical phantom is described. This dataset was acquired to investigate the variation of the SPECT calibration factor, counts convergence, noise and uniformity by varying the number of subsets and iterations. In particular, the dataset is composed of images reconstructed using five different numbers of subsets and sixteen different numbers of iterations, for a total of 80 different configurations. The dataset is suitable for comparison with other reconstruction algorithms (e.g. FBP, MLEM, etc.) and radionuclides (e.g. technetium, yttrium). In regards to the uniformity issue, the same dataset allows the user to perform radiomic investigations on the influence of the border effect on the reconstructed images.

Measurement(s)	Calibration Factor • coefficient of variation • Texture Descriptor
Technology Type(s)	Single Photon Emission Computed Tomography • Image Feature

EANM dosimetry committee recommendations for dosimetry of 177Lu-labelled somatostatin-receptor- and PSMA-targeting ligands

The purpose of the EANM Dosimetry Committee is to provide recommendations and guidance to scientists and clinicians on patient-specific dosimetry. Radiopharmaceuticals labelled with lutetium-177 (177Lu) are increasingly used for therapeutic applications, in particular for the treatment of metastatic neuroendocrine tumours using ligands for somatostatin receptors and prostate adenocarcinoma with small-molecule PSMA-targeting ligands. This paper provides an overview of reported dosimetry data for these therapies and summarises current knowledge about radiation-induced side effects on normal tissues and dose-effect relationships for tumours. Dosimetry methods and data are summarised for kidneys, bone marrow, salivary glands, lacrimal glands, pituitary glands, tumours, and the skin in case of radiopharmaceutical extravasation. Where applicable, taking into account the present status of the field and recent evidence in the literature, guidance is provided. The purpose of these recommendations is to encourage the practice of patient-specific dosimetry in therapy with 177Lu-labelled compounds. The proposed methods should be within the scope of centres offering therapy with 177Lu-labelled ligands for somatostatin receptors or small-molecule PSMA.

Quantitative SPECT/CT for dosimetry of peptide receptor radionuclide therapy

Neuroendocrine tumors (NETs) are uncommon malignancies of increasing incidence and prevalence. As these slow growing tumors usually overexpress somatostatin receptors (SSTRs), the use of ⁶⁸Ga-DOTA-peptides (gallium-68 chelated with dodecane tetra-acetic acid to somatostatin), which bind to the SSTRs, allows for PET based imaging and selection of patients for peptide receptor radionuclide therapy (PRRT). PRRT with radiolabeled somatostatin analogues such as ¹⁷⁷Lu-DOTATATE (lutetium-177-[DOTA,Tyr³]-octreotate), is mainly used for the treatment of metastatic or inoperable NETs. However, PRRT is generally administered at a fixed injected activity in order not to exceed dose limits in critical organs, which is suboptimal given the variability in radiopharmaceutical uptake among patients. Advances in SPECT (single photon emission computed tomography) imaging enable the absolute quantitative measure of the true radiopharmaceutical distribution providing for PRRT dosimetry in each patient. Personalized PRRT based on patient-specific dosimetry could improve therapeutic efficacy by optimizing effective tumor absorbed dose while limiting treatment related radiotoxicity.

Imaging and dosimetry for alpha-particle emitter radiopharmaceutical therapy: improving radiopharmaceutical therapy by looking into the black box

Radiopharmaceutical therapy using α-particle emitting radionuclides (αRPT) is a novel treatment modality that delivers highly potent alpha-particles to cancer cells or their environment. We review the advantages and challenges of imaging and dosimetry in implementing αRPT for cancer patients.

Description of the methodology for dosimetric quantification in treatments with 177Lu-DOTATATE

Implementation of dosimetry calculations in the daily practice of Nuclear Medicine Departments is, at this time, a controversial issue, partly due to the lack of a standardized methodology that is accepted by all interested parties (patients, nuclear medicine physicians and medical physicists). However, since the publication of RD 601/2019 there is a legal obligation to implement it, despite the fact that it is a complex and high resource consumption procedure. The aim of this article is to review the theoretical bases of in vivo dosimetry in treatments with 177Lu-DOTATATE. The exposed methodology is the one proposed by the MIRD Committee (Medical Internal Radiation Dose) of the SNMMI (Society of Nuclear Medicine & Molecular Imaging). According to this method, the absorbed dose is obtained as the product of 2factors: the time-integrated activity of the radiopharmaceutical present in a source region and a geometrical factor S. This approach, which a priori seems simple, in practice requires several SPECT/CT acquisitions, several measurements of the whole body activity and taking several blood samples, as well as hours of image processing and computation. The systematic implementation of these calculations, in all the patients we treat, will allow us to obtain homogeneous data to correlate the absorbed doses in the lesions with the biological effect of the treatment. The final purpose of the dosimetry calculations is to be able to maximize the therapeutic effect in the lesions, controlling the radiotoxicity in the organs at risk.

The effect of calibration factors and recovery coefficients on 177Lu SPECT activity quantification accuracy: a Monte Carlo study

BACKGROUND

Different gamma camera calibration factor (CF) geometries have been proposed to convert SPECT data into units of activity concentration. However, no consensus has been reached on a standardised geometry. The CF is dependent on the selected geometry and is further affected by partial volume effects. This study investigated the effect of two CF geometries and their corresponding recovery coefficients (RCs) on the quantification accuracy of 177Lu SPECT images using Monte Carlo simulations.

METHODS

The CF geometries investigated were (i) a radioactive-sphere surrounded by non-radioactive water (sphere-CF) and (ii) a cylindrical phantom uniformly filled with radioactive water (cylinder-CF). Recovery coefficients were obtained using the sphere-CF and cylinder-CF, yielding the sphere-RC and cylinder-RC values, respectively, for partial volume correction (PVC). The quantification accuracy was evaluated using four different-sized spheres (15.6-65.4 ml) and a kidney model with known activity concentrations inside a cylindrical, torso and patient phantom. Images were reconstructed with the 3D OS-EM algorithm incorporating attenuation, scatter and detector-response corrections. Segmentation was performed using the physical size and a small cylindrical volume inside the cylinder for the sphere-CF and cylinder-CF, respectively.

RESULTS

The sphere quantification error (without PVC) was better for the sphere-CF (≤ - 5.54%) compared to the cylinder-CF (≤ - 20.90%), attributed to the similar geometry of the quantified and CF spheres. Partial volume correction yielded comparable results for the sphere-CF-RC (≤ 3.47%) and cylinder-CF-RC (≤ 3.53%). The accuracy of the kidney quantification was poorer (≤ 22.34%) for the sphere-CF without PVC compared to the cylinder-CF (≤ 2.44%). With PVC, the kidney quantification results improved and compared well for the sphere-CF-RC (≤ 3.50%) and the cylinder-CF-RC (≤ 3.45%).

CONCLUSION

The study demonstrated that upon careful selection of CF-RC combinations, comparable quantification errors (≤ 3.53%) were obtained between the sphere-CF-RC and cylinder-CF-RC, when all corrections were applied.

177 Lutetium SPECT/CT: Evaluation of collimator, photopeak and scatter correction

Purpose

The goal of this study was to find the optimal combination of collimator, photopeak and scatter correction for ¹⁷⁷Lutetium (¹⁷⁷Lu) SPECT/CT imaging.

Methods

Three experiments [sphere‐to‐background ratios (SBR) 50:1, 10:1, and 2:1] were performed with the NEMA Image Quality phantom filled with ¹⁷⁷Lu‐trichloride. SPECT/CT acquisitions were performed with the medium‐energy low‐penetration (MELP) collimator and ^99mTc/Krypton collimator. For each acquisition six reconstructions, all with attenuation correction (AC), were made: the 113‐keV photopeak only, the 208‐keV photopeak only and both photopeaks combined, each with or without scatter correction (SC). Image quality was assessed using contrast‐to‐noise ratios (CNR), quantification accuracy by means of recovery coefficients (RCs) and the spatial resolution using line profiles.

Results

With SBR 50:1 and 10:1, both collimators met the Rose criterion (CNR > 5), whereas the MELP collimator showed a higher CNR for the 2:1 ratio. The RC_mean was higher with the MELP collimator, most explicit after the 208‐keV AC/SC reconstruction for all acquisitions. The line profiles showed a better spatial resolution for the MELP collimator and the 208‐keV AC/SC reconstructions.

Conclusion

¹⁷⁷Lu SPECT/CT image quality and quantification was most optimal when acquired with the MELP collimator and reconstructed using the 208‐keV photopeak, with AC and SC.

Dosimetry and safety of 177Lu PSMA-617 along with polyglutamate parotid gland protector: preliminary results in metastatic castration-resistant prostate cancer patients

PURPOSE

Radioligand therapy (RLT) with ¹⁷⁷Lu-PSMA-617 is a promising option for patients with metastatic castration-resistant prostate cancer (mCRPC). The present study was designed to define the safety and initial response to a minimal effective injected activity/cycle of ¹⁷⁷Lu-PSMA-617 in mCRPC patients. New protective agents for salivary glands and kidney were co-administered and dosimetry was carried out.

PATIENTS AND METHODS

A prospective single-arm, open label phase II study on mCRPC was activated at our institute in April 2017. Patients with histologically confirmed advanced mCRPC previously treated with standard life-prolonging agents were enrolled. Folic polyglutamate tablets were orally administered as parotid gland protectors and 500 mL of a 10% mannitol solution was intravenously infused to reduce kidney uptake before the injection of 3.7-5.5 GBq of ¹⁷⁷Lu-PSMA-617 repeated four times at interval of 8 weeks. The adsorbed dose calculation was performed with MIRD formalism (OLINDA/EXM software). The Bryant and Day design was used to estimate the sample size taking account of both activity and toxicity.

RESULTS

Forty-three eligible patients were evaluated for toxicity and initial response. Dosimetry was carried out in 13 patients. Two (4.8%) patients had G3 and 8 (19.5%) had G2 hematological toxicity. Only 3 (6.9%) patients had mild G1 salivary gland toxicity and 8 (19.5%) had G1 renal toxicity. A decrease of ≥ 30% in prostate-specific antigen (PSA) was achieved after the first cycle in 17 (40.5%) patients, of whom 13 had a PSA decline of >50% after the second cycle. The median adsorbed doses were 0.65 mGy/MBq (range 0.33-2.63) for parotid glands, 0.42 mGy/MBq (0.14-0.81) for kidneys, 0.036 mGy/MBq (0.023-0.067) for red marrow, and 0.038 mGy/MBq (0.018-0.135) for the whole body.

CONCLUSION

In advanced, heavily pre-treated mCRPC patients, 3.7 GBq/cycle of ¹⁷⁷Lu-PSMA-617 was safe and produced early biochemical and imaging responses at PSMA whole-body scan post injection. Dosimetry of salivary glands suggests that the co-administration of polyglutamate tablets may reduce salivary gland uptake.

CLINICAL TRIAL REGISTRATION

EU Clinical Trials Register No.: 2016-002732-32; NCT03454750. Collection and assembly of data: April 2017 and February 2019.

Preclinical Voxel-Based Dosimetry in Theranostics: a Review

Due to the increasing use of preclinical targeted radionuclide therapy (TRT) studies for the development of novel theranostic agents, several studies have been performed to accurately estimate absorbed doses to mice at the voxel level using reference mouse phantoms and Monte Carlo (MC) simulations. Accurate dosimetry is important in preclinical theranostics to interpret radiobiological dose-response relationships and to translate results for clinical use. Direct MC (DMC) simulation is believed to produce more realistic voxel-level dose distribution with high precision because tissue heterogeneities and nonuniform source distributions in patients or animals are considered. Although MC simulation is considered to be an accurate method for voxel-based absorbed dose calculations, it is time-consuming, computationally demanding, and often impractical in daily practice. In this review, we focus on the current status of voxel-based dosimetry methods applied in preclinical theranostics and discuss the need for accurate and fast voxel-based dosimetry methods for pretherapy absorbed dose calculations to optimize the dose computation time in preclinical TRT.

Quantitative SPECT/CT-Technique and Clinical Applications

The continuous development of SPECT over the past 50 years has led to improved image quality and increased diagnostic confidence. The most influential developments include the realization of hybrid SPECT/CT devices, as well as the implementation of attenuation correction and iterative image reconstruction techniques. These developments have led to a preference for SPECT/CT devices over SPECT-only systems and to the widespread adoption of the former, strengthening the role of SPECT/CT as the workhorse of Nuclear Medicine imaging. New trends in the ongoing development of SPECT/CT are diverse. For example, whole-body SPECT/CT images, consisting of acquisitions from multiple consecutive bed positions in the manner of PET/CT, are increasingly performed. Additionally, in recent years, some interesting approaches in detector technology have found their way into commercial products. For example, some SPECT cameras dedicated to specific organs employ semiconductor detectors made of cadmium telluride or cadmium zinc telluride, which have been shown to increase the obtainable image quality by offering a higher sensitivity and energy resolution. However, the advent of quantitative SPECT/CT which, like PET, can quantify the amount of tracer in terms of Bq/mL or as a standardized uptake value could be regarded as most important development. It is a major innovation that will lead to increased diagnostic accuracy and confidence, especially in longitudinal studies and in the monitoring of treatment response. The current work comprises two main aspects. At first, physical and technical fundamentals of SPECT image formation are described and necessary prerequisites of quantitative SPECT/CT are reviewed. Additionally, the typically achievable quantitative accuracy based on reports from the literature is given. Second, an extensive list of studies reporting on clinical applications of quantitative SPECT/CT is provided and reviewed.

Energy Window and Collimator Optimization in Lutetium-177 Single-photon Emission Computed Tomography Imaging using Monte Carlo Simulation

Introduction:

In lutetium-177 (Lu-177) single-photon emission computed tomography (SPECT) imaging, the accuracy of activity quantification is degraded by penetrated and scattered photons. We assessed the scattered photon fractions in order to determine the optimal situation and development of correction method. This study proposes to compare the image quality that can be achieved by three collimators.

Materials and Methods:

Siemens Medical System Symbia fitted with high-energy (HE), medium-energy (ME), and low-energy high-resolution collimators was simulated using the SIMIND Monte Carlo code simulation code. Counts were collected in three different main-energy window widths (20%, 15%, and 10%) for Lu-177 point source. Primary and scattered point spread functions and also geometric, penetration, scattering were drawn and analyzed.

Results:

In Lu-177 imaging, a 20% of main-energy window and ME collimator were found to be optimal. HE collimator can be used when the resolution is not required.

Conclusion:

These results provide the optimal energy window and collimator in Lu-177 SPECT imaging and will help the quantification of Lu-177.

Texture analysis in 177Lu SPECT phantom images: Statistical assessment of uniformity requirements using texture features

The purpose of this study was to apply texture analysis (TA) to evaluate the uniformity of SPECT images reconstructed with the 3D Ordered Subsets Expectation Maximization (3D-OSEM) algorithm. For this purpose, a cylindrical homogeneous phantom filled with ¹⁷⁷Lu was used and a total of 24 spherical volumes of interest (VOIs) were considered inside the phantom. The location of the VOIs was chosen in order to define two different configurations, i.e. gravity and radial configuration. The former configuration was used to investigate the uniformity of distribution of ¹⁷⁷Lu inside the phantom, while the latter configuration was used to investigate the lack of uniformity from center towards edge of the images. For each VOI, the trend of different texture features considered as a function of 3D-OSEM updates was investigated in order to evaluate the influence of reconstruction parameters. TA was performed using CGITA software. The equality of the average texture feature trends in both spatial configurations was assumed as the null hypothesis and was tested by functional analysis of variance (fANOVA). With regard to the gravity configuration, no texture feature rejected the null hypothesis when the number of subsets increased. For the radial configuration, the statistical analysis revealed that, depending on the 3D-OSEM parameters used, a few texture features were capable of detecting the non-uniformity of ¹⁷⁷Lu distribution inside the phantom moving from the center of the image towards its edge. Finally, cross-correlation coefficients were calculated to better identify the features that could play an important role in assessing quality assurance procedures performed on SPECT systems.