GATE based Monte Carlo simulation of planar scintigraphy to estimate the nodular dose in radioiodine therapy for autonomous thyroid adenoma

Abstract The recommended target dose in radioiodine therapy of solitary hyperfunctioning thyroid nodules is 300-400 Gy and therefore higher than in other radiotherapies. This is due to the fact that an unknown, yet significant portion of the activity is stored in extranodular areas but is neglected in the calculatory dosimetry. We investigate the feasibility of determining the ratio of nodular and extranodular activity concentrations (uptakes) from post-therapeutically acquired planar scintigrams with Monte Carlo simulations in GATE. The geometry of a gamma camera with a high energy collimator was emulated in GATE (Version 5). A geometrical thyroid-neck phantom (GP) and the ICRP reference voxel phantoms “Adult Female” (AF, 16 ml thyroid) and “Adult Male” (AM, 19 ml thyroid) were used as source regions. Nodules of 1 ml and 3 ml volume were placed in the phantoms. For each phantom and each nodule 200 scintigraphic acquisitions were simulated. Uptake ratios of nodule and rest of thyroid ranging from 1 to 20 could be created by summation. Quantitative image analysis was performed by investigating the number of simulated counts in regions of interest (ROIs). ROIs were created by perpendicular projection of the phantom onto the camera plane to avoid a user dependant bias. The ratio of count densities in ROIs over the nodule and over the contralateral lobe, which should be least affected by nodular activity, was taken to be the best available measure for the uptake ratios. However, the predefined uptake ratios are underestimated by these count density ratios: For an uptake ratio of 20 the count ratios range from 4.5 (AF, 1 ml nodule) to 15.3 (AM, 3 ml nodule). Furthermore, the contralateral ROI is more strongly affected by nodular activity than expected: For an uptake ratio of 20 between nodule and rest of thyroid up to 29% of total counts in the ROI over the contralateral lobe are caused by decays in the nodule (AF 3 ml). In the case of the 1 ml nodules this effect is smaller: 9-11% (AF) respectively 7-8% (AM). For each phantom, the dependency of count density ratios upon uptake ratios can be modeled well by both linear and quadratic regression (quadratic: r 2  > 0.99), yielding sets of parameters which in reverse allow the computation of uptake ratios (and thus dose) from count density ratios. A single regression model obtained by fitting the data of all simulations simultaneously did not provide satisfactory results except for GP, while underestimating the true uptake ratios in AF and overestimating them in AM. The scintigraphic count density ratios depend upon the uptake ratios between nodule and rest of thyroid, upon their volumes, and their respective position in a non-trivial way. Further investigations are required to derive a comprehensive rule to calculate the uptake or dose ratios based on post-therapeutic scintigraphy.