Feasibility and accuracy of UF/NCI phantoms and Monte Carlo retrospective dosimetry in children treated on National Wilms Tumor Study protocols

A Historical Survey of Key Epidemiological Studies of Ionizing Radiation Exposure

In this article we review the history of key epidemiological studies of populations exposed to ionizing radiation. We highlight historical and recent findings regarding radiation-associated risks for incidence and mortality of cancer and non-cancer outcomes with emphasis on study design and methods of exposure assessment and dose estimation along with brief consideration of sources of bias for a few of the more important studies. We examine the findings from the epidemiological studies of the Japanese atomic bomb survivors, persons exposed to radiation for diagnostic or therapeutic purposes, those exposed to environmental sources including Chornobyl and other reactor accidents, and occupationally exposed cohorts. We also summarize results of pooled studies. These summaries are necessarily brief, but we provide references to more detailed information. We discuss possible future directions of study, to include assessment of susceptible populations, and possible new populations, data sources, study designs and methods of analysis.

Harmonica consensus, controversies, and future directions in radiotherapy for pediatric Wilms tumors

Radiotherapy (RT) is essential for multimodality treatment of pediatric renal tumors, particularly in higher-risk and metastatic disease. Despite decades of use, particularly for Wilms tumor, there remain controversies regarding RT indications, timing, dose, and targets. To align global management, we address these issues in this international HARMONIsation and CollAboration (HARMONICA) project. There are multiple knowledge gaps and opportunities for future research including: (1) utilization of advanced RT technologies, including intensity-modulated RT, proton beam therapy, combined with image-guided RT to reduce target volumes; (2) impact of molecular biomarkers including loss of heterozygosity at 1p, 16q, and 1q gain on RT indications; (3) mitigation of reproductive toxicity following RT; (4) promotion of RT late effects research; and (5) support to overcome challenges in RT utilization in low- and middle-income countries where 90% of the world's children reside. Here, we outline current status and future directions for RT in pediatric renal tumors.

Body region-specific 3D age-scaling functions for scaling whole-body computed tomography anatomy for pediatric late effects studies

Purpose.Radiation epidemiology studies of childhood cancer survivors treated in the pre-computed tomography (CT) era reconstruct the patients' treatment fields on computational phantoms. For such studies, the phantoms are commonly scaled to age at the time of radiotherapy treatment because age is the generally available anthropometric parameter. Several reference size phantoms are used in such studies, but reference size phantoms are only available at discrete ages (e.g.: newborn, 1, 5, 10, 15, and Adult). When such phantoms are used for RT dose reconstructions, the nearest discrete-aged phantom is selected to represent a survivor of a specific age. In this work, we (1) conducted a feasibility study to scale reference size phantoms at discrete ages to various other ages, and (2) evaluated the dosimetric impact of using exact age-scaled phantoms as opposed to nearest age-matched phantoms at discrete ages.Methods.We have adopted the University of Florida/National Cancer Institute (UF/NCI) computational phantom library for our studies. For the feasibility study, eight male and female reference size UF/NCI phantoms (5, 10, 15, and 35 years) were downscaled to fourteen different ages which included next nearest available lower discrete ages (1, 5, 10 and 15 years) and the median ages at the time of RT for Wilms' tumor (3.9 years), craniospinal (8.0 years), and all survivors (9.1 years old) in the Childhood Cancer Survivor Study (CCSS) expansion cohort treated with RT. The downscaling was performed using our in-house age scaling functions (ASFs). To geometrically validate the scaling, Dice similarity coefficient (DSC), mean distance to agreement (MDA), and Euclidean distance (ED) were calculated between the scaled and ground-truth discrete-aged phantom (unscaled UF/NCI) for whole-body, brain, heart, liver, pancreas, and kidneys. Additionally, heights of the scaled phantoms were compared with ground-truth phantoms' height, and the Centers for Disease Control and Prevention (CDC) reported 50th percentile height. Scaled organ masses were compared with ground-truth organ masses. For the dosimetric assessment, one reference size phantom and seventeen body-size dependent 5-year-old phantoms (9 male and 8 female) of varying body mass indices (BMI) were downscaled to 3.9-year-old dimensions for two different radiation dose studies. For the first study, we simulated a 6 MV photon right-sided flank field RT plan on a reference size 5-year-old and 3.9-year-old (both of healthy BMI), keeping the field size the same in both cases. Percent of volume receiving dose ≥15 Gy (V₁₅) and the mean dose were calculated for the pancreas, liver, and stomach. For the second study, the same treatment plan, but with patient anatomy-dependent field sizes, was simulated on seventeen body-size dependent 5- and 3.9-year-old phantoms with varying BMIs. V₁₅, mean dose, and minimum dose received by 1% of the volume (D₁), and by 95% of the volume (D₉₅) were calculated for pancreas, liver, stomach, left kidney (contralateral), right kidney, right and left colons, gallbladder, thoracic vertebrae, and lumbar vertebrae. A non-parametric Wilcoxon rank-sum test was performed to determine if the dose to organs of exact age-scaled and nearest age-matched phantoms were significantly different (p < 0.05).Results.In the feasibility study, the best DSCs were obtained for the brain (median: 0.86) and whole-body (median: 0.91) while kidneys (median: 0.58) and pancreas (median: 0.32) showed poorer agreement. In the case of MDA and ED, whole-body, brain, and kidneys showed tighter distribution and lower median values as compared to other organs. For height comparison, the overall agreement was within 2.8% (3.9 cm) and 3.0% (3.2 cm) of ground-truth UF/NCI and CDC reported 50th percentile heights, respectively. For mass comparison, the maximum percent and absolute differences between the scaled and ground-truth organ masses were within 31.3% (29.8 g) and 211.8 g (16.4%), respectively (across all ages). In the first dosimetric study, absolute difference up to 6% and 1.3 Gy was found for V₁₅and mean dose, respectively. In the second dosimetric study, V₁₅and mean dose were significantly different (p < 0.05) for all studied organs except the fully in-beam organs. D₁and D₉₅were not significantly different for most organs (p > 0.05).Conclusion.We have successfully evaluated our ASFs by scaling UF/NCI computational phantoms from one age to another age, which demonstrates the feasibility of scaling any CT-based anatomy. We have found that dose to organs of exact age-scaled and nearest aged-matched phantoms are significantly different (p < 0.05) which indicates that using the exact age-scaled phantoms for retrospective dosimetric studies is a better approach.

Automatic generation of three-dimensional dose reconstruction data for two-dimensional radiotherapy plans for historically treated patients

Performing large-scale three-dimensional radiation dose reconstruction for patients requires a large amount of manual work. We present an image processing-based pipeline to automatically reconstruct radiation dose. The pipeline was designed for childhood cancer survivors that received abdominal radiotherapy with anterior-to-posterior and posterior-to-anterior field set-up. First, anatomical landmarks are automatically identified on two-dimensional radiographs. Second, these landmarks are used to derive parameters to emulate the geometry of the plan on a surrogate computed tomography. Finally, the plan is emulated and used as input for dose calculation. For qualitative evaluation, 100 cases of automatic and manual plan emulations were assessed by two experienced radiation dosimetrists in a blinded comparison. The two radiation dosimetrists approved 100%/100% and 92%/91% of the automatic/manual plan emulations, respectively. Similar approval rates of 100% and 94% hold when the automatic pipeline is applied on another 50 cases. Further, quantitative comparisons resulted in on average < 5    mm difference in plan isocenter/borders, and < 0.9    Gy in organ mean dose (prescribed dose: 14.4 Gy) calculated from the automatic and manual plan emulations. No statistically significant difference in terms of dose reconstruction accuracy was found for most organs at risk. Ultimately, our automatic pipeline results are of sufficient quality to enable effortless scaling of dose reconstruction data generation.

Conversion of computational human phantoms into DICOM-RT for normal tissue dose assessment in radiotherapy patients

Radiotherapy (RT) treatment planning systems (TPS) are designed for the fast calculation of dose to the tumor bed and nearby organs at risk using x-ray computed tomography (CT) images. However, CT images for a patient are typically available for only a small portion of the body, and in some cases, such as for retrospective epidemiological studies, no images may be available at all. When dose to organs that lie out-of-scan must be estimated, a convenient alternative for the unknown patient anatomy is to use a matching whole-body computational phantom as a surrogate. The purpose of the current work is to connect such computational phantoms to commercial RT TPS for retrospective organ dose estimation. A custom software with graphical user interface (GUI), called the DICOM-RT Generator, was developed in MATLAB to convert voxel computational phantoms into the digital imaging and communications in medicine radiotherapy (DICOM-RT) format, compatible with commercial TPS. DICOM CT image sets for the phantoms are created via a density-to-Hounsfield unit (HU) conversion curve. Accompanying structure sets containing the organ contours are automatically generated by tracing binary masks of user-specified organs on each phantom CT slice. The software was tested on a library of body size-dependent phantoms, the International Commission on Radiological Protection reference phantoms, and a canine voxel phantom, taking only a few minutes per conversion. The resulting DICOM-RT files were tested on several commercial TPS. As an example application, a library of converted phantoms was used to estimate organ doses for members of the National Wilms Tumor Study (NWTS) cohort. The converted phantom library, in DICOM format, and a standalone MATLAB-compiled executable of the DICOM-RT Generator are available for others to use for research purposes (http://ncidose.cancer.gov).