Development of an age-scalable 3D computational phantom in DICOM standard for late effects studies of childhood cancer survivors

Physics Considerations for Evaluation of Dose for Dose-Response Models of Pediatric Late Effects From Radiation Therapy: A PENTEC Introductory Review

PURPOSE

We describe the methods used to estimate the accuracy of dosimetric data found in literature sources used to construct the Pediatric Normal Tissue Effects in the Clinic (PENTEC) dose-response models, summarize these findings of each organ-specific task force, describe some of the dosimetric challenges and the extent to which these efforts affected the final modeling results, and provide guidance on the interpretation of the dose-response results given the various dosimetric uncertainties.

METHODS AND MATERIALS

Each of the PENTEC task force medical physicists reviewed all the journal articles used for dose-response modeling to identify, categorize, and quantify dosimetric uncertainties. These uncertainties fell into 6 broad categories. A uniform nomenclature was developed for describing the "dosimetric quality" of the articles used in the PENTEC reviews. Among the multidisciplinary experts in the PENTEC effort, the medical physicists were charged with the dosimetric evaluation, as they are most expert in this subject.

RESULTS

The percentage dosimetric uncertainty was estimated for each late effect endpoint for all PENTEC organ reports. Twelve specific sources of dose uncertainty were identified related to the 6 broad categories. The most common reason for organ dose uncertainty was that prescribed dose rather than organ dose was reported. Percentage dose uncertainties ranged from 5% to 200%. Systematic uncertainties were used to correct the dose component of the models. Random uncertainties were also described in each report and in some cases used to modify dose axis error bars. In addition, the potential effects of dose binning were described.

CONCLUSIONS

PENTEC reports are designed to provide guidance to radiation oncologists and treatment planners for organ dose constraints. It is critical that these dose constraint recommendations are as accurate as possible, acknowledging the large error bars for many. Achieving this accuracy is important as it enables clinicians to better balance target dose coverage with risk of late effects. Evidence-based dose constraints for pediatric patients have been lacking and, in this regard, PENTEC fills an important unmet need. One must be aware of the limitations of our recommendations, and that for some organ systems, large uncertainties exist in the dose-response model because of clinical endpoint uncertainty, dosimetric uncertainty, or both.

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.

Late-onset kidney failure in survivors of childhood cancer: a report from the Childhood Cancer Survivor Study

BACKGROUND

The incidence of and risk factors for late-onset kidney failure among survivors over the very long term remains understudied.

MATERIALS AND METHODS

A total of 25,530 childhood cancer survivors (median follow-up 22.3 years, interquartile range 17.4-28.8) diagnosed between 1970 and 1999, and 5045 siblings from the Childhood Cancer Survivor Study were assessed for self-reported late-onset kidney failure, defined as dialysis, renal transplantation, or death attributable to kidney disease. Piecewise exponential models evaluated associations between risk factors and the rate of late-onset kidney failure.

RESULTS

A total of 206 survivors and 10 siblings developed late-onset kidney failure, a 35-year cumulative incidence of 1.7% (95% confidence interval [CI] = 1.4-1.9) and 0.2% (95% confidence interval [CI] = 0.1-0.4), respectively, corresponding to an adjusted rate ratio (RR) of 4.9 (95% CI = 2.6-9.2). High kidney dose from radiotherapy (≥15Gy; RR = 4.0, 95% CI = 2.1-7.4), exposure to high-dose anthracycline (≥250 mg/m2; RR = 1.6, 95% CI = 1.0-2.6) or any ifosfamide chemotherapy (RR = 2.6, 95% CI = 1.2-5.7), and nephrectomy (RR = 1.9, 95% CI = 1.0-3.4) were independently associated with elevated risk for late-onset kidney failure among survivors. Survivors who developed hypertension, particularly in the context of prior nephrectomy (RR = 14.4, 95% CI = 7.1-29.4 hypertension with prior nephrectomy; RR = 5.9, 95% CI = 3.3-10.5 hypertension without prior nephrectomy), or diabetes (RR = 2.2, 95%CI = 1.2-4.2) were also at elevated risk for late-onset kidney failure.

CONCLUSIONS

Survivors of childhood cancer are at increased risk for late-onset kidney failure. Kidney dose from radiotherapy ≥15 Gy, high-dose anthracycline, any ifosfamide, and nephrectomy were associated with increased risk of late-onset kidney failure among survivors. Successful diagnosis and management of modifiable risk factors such as diabetes and hypertension may mitigate the risk for late-onset kidney failure. The association of late-onset kidney failure with anthracycline chemotherapy represents a novel finding that warrants further study.

Radiation therapy related cardiac disease risk in childhood cancer survivors: Updated dosimetry analysis from the Childhood Cancer Survivor Study

BACKGROUND AND PURPOSE

We previously evaluated late cardiac disease in long-term survivors in the Childhood Cancer Survivor Study (CCSS) based on heart radiation therapy (RT) doses estimated from an age-scaled phantom with a simple atlas-based heart model (H_Atlas). We enhanced our phantom with a high-resolution CT-based anatomically realistic and validated age-scalable cardiac model (H_Hybrid). We aimed to evaluate how this update would impact our prior estimates of RT-related late cardiac disease risk in the CCSS cohort.

METHODS

We evaluated 24,214 survivors from the CCSS diagnosed from 1970 to 1999. RT fields were reconstructed on an age-scaled phantom with H_Hybrid and mean heart dose (D_m), percent volume receiving ≥ 20 Gy (V₂₀) and ≥ 5 Gy with V₂₀ = 0 ( [Formula: see text] ) were calculated. We reevaluated cumulative incidences and adjusted relative rates of grade 3-5 Common Terminology Criteria for Adverse Events outcomes for any cardiac disease, coronary artery disease (CAD), and heart failure (HF) in association with D_m, V₂₀, and [Formula: see text] (as categorical variables). Dose-response relationships were evaluated using piecewise-exponential models, adjusting for attained age, sex, cancer diagnosis age, race/ethnicity, time-dependent smoking history, diagnosis year, and chemotherapy exposure and doses. For relative rates, D_m was also considered as a continuous variable.

RESULTS

Consistent with previous findings with H_Atlas, reevaluation using H_Hybrid dosimetry found that, D_m ≥ 10 Gy, V₂₀ ≥ 0.1%, and [Formula: see text]  ≥ 50% were all associated with increased cumulative incidences and relative rates for any cardiac disease, CAD, and HF. While updated risk estimates were consistent with previous estimates overall without statistically significant changes, there were some important and significant (P < 0.05) increases in risk with updated dosimetry for D_m in the category of 20 to 29.9 Gy and V₂₀ in the category of 30% to 79.9%. When changes in the linear dose-response relationship for D_m were assessed, the slopes of the dose response were steeper (P < 0.001) with updated dosimetry. Changes were primarily observed among individuals with chest-directed RT with prescribed doses ≥ 20 Gy.

CONCLUSION

These findings present a methodological advancement in heart RT dosimetry with improved estimates of RT-related late cardiac disease risk. While results are broadly consistent with our prior study, we report that, with updated cardiac dosimetry, risks of cardiac disease are significantly higher in two dose and volume categories and slopes of the Dm-specific RT-response relationships are steeper. These data support the use of contemporary RT to achieve lower heart doses for pediatric patients, particularly those requiring chest-directed RT.

Development and validation of an age-scalable cardiac model with substructures for dosimetry in late-effects studies of childhood cancer survivors

BACKGROUND AND PURPOSE

Radiation therapy is a risk factor for late cardiac disease in childhood cancer survivors. Several pediatric cohort studies have established whole heart dose and dose-volume response models. Emerging data suggest that dose to cardiac substructures may be more predictive than whole heart metrics. In order to develop substructure dose-response models, the heart model previously used for pediatric cohort dosimetry needed enhancement and substructure delineation.

METHODS

To enhance our heart model, we combined the age-scalable capability of our computational phantom with the anatomically-delineated (with substructures) heart models from an international humanoid phantom series. We examined cardiac volume similarity/overlap between registered age-scaled phantoms (1, 5, 10, and 15 years) with the enhanced heart model and the reference phantoms of the same age; dice similarity coefficient (DSC) and overlap coefficient (OC) were calculated for each matched pair. To assess the accuracy of our enhanced heart model, we compared doses from computed tomography-based planning (ground truth) with reconstructed heart doses. We also compared doses calculated with the prior and enhanced heart models for a cohort of nearly 5000 childhood cancer survivors.

RESULTS

We developed a realistic cardiac model with 14-substructures, scalable across a broad age range (1-15 years); average DSC and OC were 0.84 ± 0.05 and 0.90 ± 0.05, respectively. The average percent difference between reconstructed and ground truth mean heart doses was 4.2%. In the cohort dosimetry analysis, dose and dose-volume metrics were approximately 10% lower on average when the enhanced heart model was used for dose reconstructions.

CONCLUSION

We successfully developed and validated an anatomically realistic age-scalable cardiac model that can be used to establish substructure dose-response models for late cardiac disease in childhood cancer survivor cohorts.