Quality Assurance Considerations in Radiopharmaceutical Therapy Dosimetry Using PLANETDose: An International Atomic Energy Agency Study.

Implementation of radiopharmaceutical therapy dosimetry varies depending on the clinical application, dosimetry protocol, software, and ultimately the operator. Assessing clinical dosimetry accuracy and precision is therefore a challenging task. This work emphasizes some pitfalls encountered during a structured analysis, performed on a single-patient dataset consisting of SPECT/CT images by various participants using a standard protocol and clinically approved commercial software. Methods: The clinical dataset consisted of the dosimetric study of a patient administered with [177Lu]Lu-DOTATATE at Tygerberg Hospital, South Africa, as a part of International Atomic Energy Agency-coordinated research project E23005. SPECT/CT images were acquired at 5 time points postinjection. Patient and calibration images were reconstructed on a workstation, and a calibration factor of 122.6 Bq/count was derived independently and provided to the participants. A standard dosimetric protocol was defined, and PLANETDose (version 3.1.1) software was installed at 9 centers to perform the dosimetry of 3 treatment cycles. The protocol included rigid image registration, segmentation (semimanual for organs, activity threshold for tumors), and dose voxel kernel convolution of activity followed by absorbed dose (AD) rate integration to obtain the ADs. Iterations of the protocol were performed by participants individually and within collective training, the results of which were analyzed for dosimetric variability, as well as for quality assurance and error analysis. Intermediary checkpoints were developed to understand possible sources of variation and to differentiate user error from legitimate user variability. Results: Initial dosimetric results for organs (liver and kidneys) and lesions showed considerable interoperator variability. Not only was the generation of intermediate checkpoints such as total counts, volumes, and activity required, but also activity-to-count ratio, activity concentration, and AD rate-to-activity concentration ratio to determine the source of variability. Conclusion: When the same patient dataset was analyzed using the same dosimetry procedure and software, significant disparities were observed in the results despite multiple sessions of training and feedback. Variations due to human error could be minimized or avoided by performing intensive training sessions, establishing intermediate checkpoints, conducting sanity checks, and cross-validating results across physicists or with standardized datasets. This finding promotes the development of quality assurance in clinical dosimetry.

Laryngeal EBV-positive Inflammatory Follicular Dendritic cell/fibroblastic Reticular cell Tumour.

Epstein-Barr virus-positive Inflammatory follicular dendritic cell/fibroblastic reticular cell tumour (EBV-IFDC/FRCT) is a rare neoplasm that occurs almost exclusively in the liver or spleen. Extra-hepatosplenic presentation is infrequent and exceptional cases have been described arising in the gastrointestinal tract or in the pharynx. However, EBV-IFDC/FRCT cases have not been previously reported in the larynx. This report describes a case of a 32-year-old woman who arrived to the emergency department due to progressive dyspnea with associated inspiratory stridor and non-productive cough. Direct laryngoscopy showed a nodular tumour arising on the left posterior subglottic mucosa obstructing 90% of the airway. A preoperative dual energy contrast enhanced computed tomography (CECT) was performed demonstrating a low attenuation lesion on virtual non-contrast (VNC) images and vivid iodine uptake on the iodine map. The tumour was excised and the histopathological analysis led to the diagnosis of an EBV-IFDC/FRCT. A fibre-optic laryngoscopy six months after the surgery did not show any abnormalities. Although the vast majority of EBV-IFDC/FRCT occur in the liver or spleen, some extra hepatosplenic tumours have been reported affecting the head and neck region. We describe here the first case arising in the larynx, as well as the usefulness of preoperative dual energy imaging techniques to assess these lesions, thus providing information that could have management implications.

Normal organ dosimetry for thyroid cancer patients treated with radioiodine as part of the multi-centre multi-national Horizon 2020 MEDIRAD project.

PURPOSE: Dosimetry is rarely performed for the treatment of differentiated thyroid cancer patients with Na[131I]I (radioiodine), and information regarding absorbed doses delivered is limited. Collection of dosimetry data in a multi-centre setting requires standardised quantitative imaging and dosimetry. A multi-national, multi-centre clinical study was performed to assess absorbed doses delivered to normal organs for differentiated thyroid cancer patients treated with Na[131I]I. METHODS: Patients were enrolled in four centres and administered fixed activities of 1.1 or 3.7 GBq of Na[131I]I using rhTSH stimulation or under thyroid hormone withdrawal according to local protocols. Patients were imaged using SPECT(/CT) at variable imaging time-points following standardised acquisition and reconstruction protocols. Whole-body retention data were collected. Dosimetry for normal organs was performed at two dosimetry centres and results collated. RESULTS: One hundred and five patients were recruited. Median absorbed doses per unit administered activity of 0.44, 0.14, 0.05 and 0.16 mGy/MBq were determined for the salivary glands of patients treated at centre 1, 2, 3 and 4, respectively. Median whole-body absorbed doses for 1.1 and 3.7 GBq were 0.05 Gy and 0.16 Gy, respectively. Median whole-body absorbed doses per unit administered activity of 0.04, 0.05, 0.04 and 0.04 mGy/MBq were calculated for centre 1, 2, 3 and 4, respectively. CONCLUSIONS: A wide range of normal organ doses were observed for differentiated thyroid cancer patients treated with Na[131I]I, highlighting the necessity for individualised dosimetry. The results show that data may be collated from multiple centres if minimum standards for the acquisition and dosimetry protocols can be achieved.

Modelling SPECT auto-contouring acquisitions for 177Lu & 131I molecular radiotherapy using new developments in Geant4/GATE.

PURPOSE: Monte Carlo modelling of SPECT imaging in Molecular Radiotherapy can improve activity quantification. Until now, SPECT modelling with GATE only considered circular orbit (CO) acquisitions. This cannot reproduce auto-contour acquisitions, where the detector head moves close to the patient to improve image resolution. The aim of this work is to develop and validate an auto-contouring step-and-shoot acquisition mode for GATE SPECT modelling. METHODS: 177Lu and 131I SPECT experimental acquisitions performed on a Siemens Symbia T2 and GE Discovery 670 gamma camera, respectively, were modelled. SPECT projections were obtained for a cylindrical Jaszczak phantom and a lung and spine phantom. Detector head parameters (radial positions and acquisition angles) were extracted from the experimental projections to model the non-circular orbit (NCO) detector motion. The gamma camera model was validated against the experimental projections obtained with the cylindrical Jaszczak (177Lu) and lung and spine phantom (131I). Then, 177Lu and 131I CO and NCO SPECT projections were simulated to validate the impact of explicit NCO modelling on simulated projections. RESULTS: Experimental and simulated SPECT images were compared using the gamma index, and were in good agreement with gamma index passing rate (GIPR) and gammaavg of 96.27%, 0.242 (177Lu) and 92.89%, 0.36 (131I). Then, simulated 177Lu and 131I CO and NCO SPECT projections were compared. The GIPR, gammaavg between the two gamma camera motions was 99.85%, 0.108 for 177Lu and 75.58%, 0.6 for 131I. CONCLUSION: This work thereby justifies the need for auto-contouring modelling for isotopes with high septal penetration.

Generation of clinical 177Lu SPECT/CT images based on Monte Carlo simulation with GATE.

PURPOSE: Patient-specific dosimetry in MRT relies on quantitative imaging, pharmacokinetic assessment and absorbed dose calculation. The DosiTest project was initiated to evaluate the uncertainties associated with each step of the clinical dosimetry workflow through a virtual multicentric clinical trial. This work presents the generation of simulated clinical SPECT datasets based on GATE Monte Carlo modelling with its corresponding experimental CT image, which can subsequently be processed by commercial image workstations. METHODS: This study considers a therapy cycle of 6.85 GBq 177Lu-labelled DOTATATE derived from an IAEA-Coordinated Research Project (E23005) on "Dosimetry in Radiopharmaceutical therapy for personalised patient treatment". Patient images were acquired on a GE Infinia-Hawkeye 4 gamma camera using a medium energy (ME) collimator. Simulated SPECT projections were generated based on experimental time points and validated against experimental SPECT projections using flattened profiles and gamma index. The simulated projections were then incorporated into the patient SPECT/CT DICOM envelopes for processing and their reconstruction within a commercial image workstation. RESULTS: Gamma index passing rate (2% - 1 pixel criteria) between 95 and 98% and average gamma between 0.28 and 0.35 among different time points revealed high similarity between simulated and experimental images. Image reconstruction of the simulated projections was successful on HERMES and Xeleris workstations, a major step forward for the initiation of a multicentric virtual clinical dosimetry trial based on simulated SPECT/CT images. CONCLUSIONS: Realistic 177Lu patient SPECT projections were generated in GATE. These modelled datasets will be circulated to different clinical departments to perform dosimetry in order to assess the uncertainties in the entire dosimetric chain.

Setting up a quantitative SPECT imaging network for a European multi-centre dosimetry study of radioiodine treatment for thyroid cancer as part of the MEDIRAD project.

BACKGROUND: Differentiated thyroid cancer has been treated with radioiodine for almost 80 years, although controversial questions regarding radiation-related risks and the optimisation of treatment regimens remain unresolved. Multi-centre clinical studies are required to ensure recruitment of sufficient patients to achieve the statistical significance required to address these issues. Optimisation and standardisation of data acquisition and processing are necessary to ensure quantitative imaging and patient-specific dosimetry. MATERIAL AND METHODS: A European network of centres able to perform standardised quantitative imaging of radioiodine therapy of thyroid cancer patients was set-up within the EU consortium MEDIRAD. This network will support a concurrent series of clinical studies to determine accurately absorbed doses for thyroid cancer patients treated with radioiodine. Five SPECT(/CT) systems at four European centres were characterised with respect to their system volume sensitivity, recovery coefficients and dead time. RESULTS: System volume sensitivities of the Siemens Intevo systems (crystal thickness 3/8″) ranged from 62.1 to 73.5 cps/MBq. For a GE Discovery 670 (crystal thickness 5/8″) a system volume sensitivity of 92.2 cps/MBq was measured. Recovery coefficients measured on three Siemens Intevo systems show good agreement. For volumes larger than 10 ml, the maximum observed difference between recovery coefficients was found to be ± 0.02. Furthermore, dead-time coefficients measured on two Siemens Intevo systems agreed well with previously published dead-time values. CONCLUSIONS: Results presented here provide additional support for the proposal to use global calibration parameters for cameras of the same make and model. This could potentially facilitate the extension of the imaging network for further dosimetry-based studies.

OpenDose: Open-Access Resource for Nuclear Medicine Dosimetry.

Radiopharmaceutical dosimetry depends on the localization in space and time of radioactive sources and requires the estimation of the amount of energy emitted by the sources deposited within targets. In particular, when computing resources are not accessible, this task can be performed using precomputed tables of specific absorbed fractions (SAFs) or S values based on dosimetric models. The aim of the OpenDose collaboration is to generate and make freely available a range of dosimetric data and tools. Methods: OpenDose brings together resources and expertise from 18 international teams to produce and compare traceable dosimetric data using 6 of the most popular Monte Carlo codes in radiation transport (EGSnrc/EGS++, FLUKA, GATE, Geant4, MCNP/MCNPX, and PENELOPE). SAFs are uploaded, together with their associated statistical uncertainties, in a relational database. S values are then calculated from monoenergetic SAFs on the basis of the radioisotope decay data presented in International Commission on Radiological Protection Publication 107. Results: The OpenDose collaboration produced SAFs for all source region and target combinations of the 2 International Commission on Radiological Protection Publication 110 adult reference models. SAFs computed from the different Monte Carlo codes were in good agreement at all energies, with SDs below individual statistical uncertainties. Calculated S values were in good agreement with OLINDA/EXM 2.0 (commercial) and IDAC-Dose 2.1 (free) software. A dedicated website (www.opendose.org) has been developed to provide easy and open access to all data. Conclusion: The OpenDose website allows the display and downloading of SAFs and the corresponding S values for 1,252 radionuclides. The OpenDose collaboration, open to new research teams, will extend data production to other dosimetric models and implement new free features, such as online dosimetric tools and patient-specific absorbed dose calculation software, together with educational resources.

A semantic database for integrated management of image and dosimetric data in low radiation dose research in medical imaging.

Medical ionizing radiation procedures and especially medical imaging are a non negligible source of exposure to patients. Whereas the biological effects of high absorbed doses are relatively well known, the effects of low absorbed doses are still debated. This work presents the development of a computer platform called Image and Radiation Dose BioBank (IRDBB) to manage research data produced in the context of the MEDIRAD project, a European project focusing on research on low doses in the context of medical procedures. More precisely, the paper describes a semantic database linking dosimetric data (such as absorbed doses to organs) to the images corresponding to X-rays exposure (such as CT images) or scintigraphic images (such as SPECT or PET images) that allow measuring the distribution of a radiopharmaceutical. The main contributions of this work are: 1) the implementation of the semantic database of the IRDBB system and 2) an ontology called OntoMEDIRAD covering the domain of discourse involved in MEDIRAD research data, especially many concepts from the DICOM standard modelled according to a realist approach.

SPECT/CT: principales aplicaciones en la medicina nuclear / SPECT/CT: main applications in nuclear medicine

El SPECT/CT ha significado no solo la posibilidad de adquirir las imágenes anatómicas y funcionales en un único estudio, sino una verdadera revolución en el manejo de muchas enfermedades a partir de lo mejor que es capaz de aportar cada una de estas modalidades de imagen. El objetivo del presente trabajo es brindar una panorámica de las principales aplicaciones del SPECT/CT en el campo de la oncología, la cardiología y la neurología. Nuevos avances tecnológicos en el diseño de novedosos detectores de estado sólido y equipos dedicados han mejorado mucho las prestaciones de este tipo de modalidad dual. Esta técnica híbrida ha mostrado que mejora la sensibilidad y la especificidad de los estudios gammagráficos, a la vez que acorta los tiempos de adquisición y brinda imágenes corregistradas corregidas por atenuación que facilitan el mejor análisis de las mismas. Entre las principales aplicaciones en el estudio de las enfermedades oncológicas se pueden mencionar las siguientes: localización y seguimiento de diferentes tipos de tumores y sus metástasis, así como para la optimización de las dosis en los estudios de radioterapia. Esta técnica ha resultado muy útil en la valoración de la enfermedad de arterias coronarias al permitir una adecuada corrección de la atenuación, determinar el score de calcio y realizar estudios de angio-CT, en dependencia de la calidad de la CT. En el estudio de algunas enfermedades neurológicas también ha ido ganando terreno el SPECT-CT. La incorporación de nuevos avances tecnológicos, así como de novedosos radiofármacos augura un peso cada vez más relevante del SPECT/CT en la práctica clínica

SPECT/CT has represented not only the possibility of acquiring anatomical and functional images in one single study, but also a revolution for the clinical management of several diseases, taking the better of each one of these imaging modalities. The present work is aimed at presenting an overview of the most important applications of the SPECT/CT in the field of oncology, cardiology and neurology. New technological advances in the design of innovative solid state detectors and related equipment have had a positive effect on the performance of this kind of dual modality. This hybrid technique improves the sensitivity and the specificity of gammagraphic studies, as well as shortens the acquisition times and gives attenuation correction of co-registered images, which, in turn, makes their analysis easier. Some of the main applications for the study of oncological diseases are the following: localization and follow-up of different kinds of tumors, their metastasis and relapses, as well as the optimization of radiotherapy doses. This technique has been useful to evaluate the coronary artery disease allowing an adequate attenuation correction of images, the determination of calcium score, and performing angio-CT studies, according to the CT quality. SPECT/CT has also gained ground in the assessment of some neurological diseases. Conclusions: The introduction of new technological advances and radiopharmaceuticals thus predicting a more relevant place for SPECT/CT in clinical practice.

Impact of SPECT corrections on 3D-dosimetry for liver transarterial radioembolization using the patient relative calibration methodology.

PURPOSE: Many centers aim to plan liver transarterial radioembolization (TARE) with dosimetry, even without CT-based attenuation correction (AC), or with unoptimized scatter correction (SC) methods. This work investigates the impact of presence vs absence of such corrections, and limited spatial resolution, on 3D dosimetry for TARE. METHODS: Three voxelized phantoms were derived from CT images of real patients with different body sizes. Simulations of (99m)Tc-SPECT projections were performed with the SIMIND code, assuming three activity distributions in the liver: uniform, inside a "liver's segment," or distributing multiple uptaking nodules ("nonuniform liver"), with a tumoral liver/healthy parenchyma ratio of 5:1. Projection data were reconstructed by a commercial workstation, with OSEM protocol not specifically optimized for dosimetry (spatial resolution of 12.6 mm), with/without SC (optimized, or with parameters predefined by the manufacturer; dual energy window), and with/without AC. Activity in voxels was calculated by a relative calibration, assuming identical microspheres and (99m)Tc-SPECT counts spatial distribution. 3D dose distributions were calculated by convolution with (90)Y voxel S-values, assuming permanent trapping of microspheres. Cumulative dose-volume histograms in lesions and healthy parenchyma from different reconstructions were compared with those obtained from the reference biodistribution (the "gold standard," GS), assessing differences for D95%, D70%, and D50% (i.e., minimum value of the absorbed dose to a percentage of the irradiated volume). γ tool analysis with tolerance of 3%/13 mm was used to evaluate the agreement between GS and simulated cases. The influence of deep-breathing was studied, blurring the reference biodistributions with a 3D anisotropic gaussian kernel, and performing the simulations once again. RESULTS: Differences of the dosimetric indicators were noticeable in some cases, always negative for lesions and distributed around zero for parenchyma. Application of AC and SC reduced systematically the differences for lesions by 5%-14% for a liver segment, and by 7%-12% for a nonuniform liver. For parenchyma, the data trend was less clear, but the overall range of variability passed from -10%/40% for a liver segment, and -10%/20% for a nonuniform liver, to -13%/6% in both cases. Applying AC, SC with preset parameters gave similar results to optimized SC, as confirmed by γ tool analysis. Moreover, γ analysis confirmed that solely AC and SC are not sufficient to obtain accurate 3D dose distribution. With breathing, the accuracy worsened severely for all dosimetric indicators, above all for lesions: with AC and optimized SC, -38%/-13% in liver's segment, -61%/-40% in the nonuniform liver. For parenchyma, D50% resulted always less sensitive to breathing and sub-optimal correction methods (difference overall range: -7%/13%). CONCLUSIONS: Reconstruction protocol optimization, AC, SC, PVE and respiratory motion corrections should be implemented to obtain the best possible dosimetric accuracy. On the other side, thanks to the relative calibration, D50% inaccuracy for the healthy parenchyma from absence of AC was less than expected, while the optimization of SC was scarcely influent. The relative calibration therefore allows to perform TARE planning, basing on D50% for the healthy parenchyma, even without AC or with suboptimal corrections, rather than rely on nondosimetric methods.

Differences in 3D dose distributions due to calculation method of voxel S-values and the influence of image blurring in SPECT.

This study compares 3D dose distributions obtained with voxel S values (VSVs) for soft tissue, calculated by several methods at their current state-of-the-art, varying the degree of image blurring. The methods were: 1) convolution of Dose Point Kernel (DPK) for water, using a scaling factor method; 2) an analytical model (AM), fitting the deposited energy as a function of the source-target distance; 3) a rescaling method (RSM) based on a set of high-resolution VSVs for each isotope; 4) local energy deposition (LED). VSVs calculated by direct Monte Carlo simulations were assumed as reference. Dose distributions were calculated considering spheroidal clusters with various sizes (251, 1237 and 4139 voxels of 3 mm size), uniformly filled with (131)I, (177)Lu, (188)Re or (90)Y. The activity distributions were blurred with Gaussian filters of various widths (6, 8 and 12 mm). Moreover, 3D-dosimetry was performed for 10 treatments with (90)Y derivatives. Cumulative Dose Volume Histograms (cDVHs) were compared, studying the differences in D95%, D50% or Dmax (ΔD95%, ΔD50% and ΔDmax) and dose profiles.For unblurred spheroidal clusters, ΔD95%, ΔD50% and ΔDmax were mostly within some percents, slightly higher for (177)Lu with DPK (8%) and RSM (12%) and considerably higher for LED (ΔD95% up to 59%). Increasing the blurring, differences decreased and also LED yielded very similar results, but D95% and D50% underestimations between 30-60% and 15-50%, respectively (with respect to 3D-dosimetry with unblurred distributions), were evidenced. Also for clinical images (affected by blurring as well), cDVHs differences for most methods were within few percents, except for slightly higher differences with LED, and almost systematic for dose profiles with DPK (-1.2%), AM (-3.0%) and RSM (4.5%), whereas showed an oscillating trend with LED.The major concern for 3D-dosimetry on clinical SPECT images is more strongly represented by image blurring than by differences among the VSVs calculation methods. For volume sizes about 2-fold the spatial resolution, D95% and D50% underestimations up to about 60 and 50% could result, so the usefulness of 3D-dosimetry is highly questionable for small tumors, unless adequate corrections for partial volume effects are adopted.