Normal Tissue Complication Probability Models Research Articles

Skin dose-volume predictors of moderate-severe late side effects after whole breast radiotherapy

BackgroundToxicity after whole breast Radiotherapy is a relevant issue, impacting the quality-of-life of a not negligible number of patients. We aimed to develop a Normal Tissue Complication Probability (NTCP) model predicting late toxicities by combining dosimetric parameters of the breast dermis and clinical factors. MethodsThe skin structure was defined as the outer CT body contour's 5 mm inner isotropic expansion. It was retrospectively segmented on a large mono-institutional cohort of early-stage breast cancer patients enrolled between 2009 and 2017 (n = 1066). Patients were treated with tangential-field RT, delivering 40 Gy in 15 fractions to the whole breast. Toxicity was reported during Follow-Up (FU) using SOMA/LENT scoring. The study endpoint was moderate-severe late side effects consisting of Fibrosis-Atrophy-Telangiectasia-Pain (FATP G ≥ 2) developed within 42 months after RT completion. A machine learning pipeline was designed with a logistic model combining clinical factors and absolute skin DVH (cc) parameters as output. ResultsThe FATP G2 + rate was 3.8 %, with 40/1066 patients experiencing side effects. After the preprocessing of variables, a cross-validation was applied to define the best-performing model. We selected a 4-variable model with Post-Surgery Cosmetic alterations (Odds Ratio, OR = 7.3), Aromatase Inhibitors (as a protective factor with OR = 0.45), V20 Gy (50 % of the prescribed dose, OR = 1.02), and V42 Gy (105 %, OR = 1.09). Factors were also converted into an adjusted V20Gy. ConclusionsThe association between late reactions and skin DVH when delivering 40 Gy/15 fr was quantified, suggesting an independent role of V20 and V42. Few clinical factors heavily modulate the risk.

Interobserver delineation variability of computed tomography-based radiomic features of the parotid gland.

To assess the interobserver delineation variability of radiomic features of the parotid gland from computed tomography (CT) images and evaluate the correlation of these features for head and neck cancer (HNC) radiotherapy patients. Contrast-enhanced CT images of 20 HNC patients were utilized. The parotid glands were delineated by treating radiation oncologists (ROs), a selected RO and AccuContour auto-segmentation software. Dice similarity coefficients (DSCs) between each pair of observers were calculated. A total of 107 radiomic features were extracted, whose robustness to interobserver delineation was assessed using the intraclass correlation coefficient (ICC). Pearson correlation coefficients (r) were calculated to determine the relationship between the features. The influence of excluding unrobust features from normal tissue complication probability (NTCP) modeling was investigated for severe oral mucositis (grade ≥3). The average DSC was 0.84 (95% confidence interval, 0.83-0.86). Most of the shape features demonstrated robustness (ICC ≥0.75), while the first-order and texture features were influenced by delineation variability. Among the three observers investigated, 42 features were sufficiently robust, out of which 36 features exhibited weak correlation (|r|<0.8). No significant difference in the robustness level was found when comparing manual segmentation by a single RO or automated segmentation with the actual clinical contour data made by treating ROs. Excluding unrobust features from the NTCP model for severe oral mucositis did not deteriorate the model performance. Interobserver delineation variability had substantial impact on radiomic features of the parotid gland. Both manual and automated segmentation methods contributed similarly to this variation.

Mechanistic model of radiotherapy-induced lung fibrosis using coupled 3D agent-based and Monte Carlo simulations

BackgroundMechanistic modelling of normal tissue toxicities is unfolding as an alternative to the phenomenological normal tissue complication probability models. The latter, currently used in the clinics, rely exclusively on limited patient data and neglect spatial dose distribution information. Among the various approaches, agent-based models are appealing as they provide the means to include patient-specific parameters and simulate long-term effects in complex systems. However, Monte Carlo tools remain the state-of-the-art for modelling radiation transport and provide measurements of the delivered dose with unmatched precision.MethodsIn this work, we develop and characterize a coupled 3D agent-based – Monte Carlo model that mechanistically simulates the onset of the radiation-induced lung fibrosis in an alveolar segment. To the best of our knowledge, this is the first such model.ResultsOur model replicates extracellular matrix patterns, radiation-induced lung fibrosis severity indexes and functional subunits survivals that show qualitative agreement with experimental studies and are consistent with our past results. Moreover, in accordance with experimental results, higher functional subunits survival and lower radiation-induced lung fibrosis severity indexes are achieved when a 5-fractions treatment is simulated. Finally, the model shows increased sensitivity to more uniform protons dose distributions with respect to more heterogeneous ones from photon irradiation.ConclusionsThis study lays thus the groundwork for further investigating the effects of different radiotherapeutic treatments on the onset of radiation-induced lung fibrosis via mechanistic modelling.

Predicting radiation pneumonitis in lung cancer: a EUD-based machine learning approach for volumetric modulated arc therapy patients.

This study aims to develop an optimal machine learning model that uses lung equivalent uniform dose (lung EUD to predict radiation pneumonitis (RP) occurrence in lung cancer patients treated with volumetric modulated arc therapy (VMAT). We analyzed a cohort of 77 patients diagnosed with locally advanced squamous cell lung cancer (LASCLC) receiving concurrent chemoradiotherapy with VMAT. Patients were categorized based on the onset of grade II or higher radiation pneumonitis (RP 2+). Dose volume histogram data, extracted from the treatment planning system, were used to compute the lung EUD values for both groups using a specialized numerical analysis code. We identified the parameter α, representing the most significant relative difference in lung EUD between the two groups. The predictive potential of variables for RP2+, including physical dose metrics, lung EUD, normal tissue complication probability (NTCP) from the Lyman-Kutcher-Burman (LKB) model, and lung EUD-calibrated NTCP for affected and whole lung, underwent both univariate and multivariate analyses. Relevant variables were then employed as inputs for machine learning models: multiple logistic regression (MLR), support vector machine (SVM), decision tree (DT), and K-nearest neighbor (KNN). Each model's performance was gauged using the area under the curve (AUC), determining the best-performing model. The optimal α-value for lung EUD was 0.3, maximizing the relative lung EUD difference between the RP 2+ and non-RP 2+ groups. A strong correlation coefficient of 0.929 (P< 0.01) was observed between lung EUD (α = 0.3) and physical dose metrics. When examining predictive capabilities, lung EUD-based NTCP for the affected lung (AUC: 0.862) and whole lung (AUC: 0.815) surpassed LKB-based NTCP for the respective lungs. The decision tree (DT) model using lung EUD-based predictors emerged as the superior model, achieving an AUC of 0.98 in both training and validation datasets. The likelihood of developing RP 2+ has shown a significant correlation with the advancements in RT technology. From traditional 3-D conformal RT, lung cancer treatment methodologies have transitioned to sophisticated techniques like static IMRT. Accurately deriving such a dose-effect relationship through NTCP modeling of RP incidence is statistically challenging due to the increased number of degrees-of-freedom. To the best of our knowledge, many studies have not clarified the rationale behind setting the α-value to 0.99 or 1, despite the closely aligned calculated lung EUD and lung mean dose MLD. Perfect independence among variables is rarely achievable in real-world scenarios. Four prominent machine learning algorithms were used to devise our prediction models. The inclusion of lung EUD-based factors substantially enhanced their predictive performance for RP 2+. Our results advocate for the decision tree model with lung EUD-based predictors as the optimal prediction tool for VMAT-treated lung cancer patients. Which could replace conventional dosimetric parameters, potentially simplifying complex neural network structures in prediction models.

MR-guided stereotactic radiotherapy of infra-diaphragmatic oligometastases: Evaluation of toxicity and dosimetric parameters

Background and purposeThe SOFT trial is a prospective, multicenter, phase 2 trial investigating magnetic resonance (MR)-guided stereotactic ablative radiotherapy (SABR) for abdominal, soft tissue metastases in patients with oligometastatic disease (OMD) (clinicaltrials.gov ID NCT04407897). We present the primary endpoint analysis of 1-year treatment-related toxicity (TRAE). Materials and methodsPatients with up to five oligometastases from non-hematological cancers were eligible for inclusion. A risk-adapted strategy prioritized fixed organs at risk (OAR) constraints over target coverage. Fractionation schemes were 45–67.5 Gy in 3–8 fractions. The primary endpoint was grade ≥ 4 TRAE within 12 months post-SABR. The association between the risk of gastrointestinal (GI) toxicity and clinical and dosimetric parameters was tested using a normal tissue complication probability model. ResultsWe included 121 patients with 147 oligometastatic targets, mainly located in the liver (41 %), lymph nodes (35 %), or adrenal glands (14 %). Nearly half of all targets (48 %, n = 71) were within 10 mm of a radiosensitive OAR. No grade 4 or 5 TRAEs, 3.5 % grade 3 TRAEs, and 43.7 % grade 2 TRAEs were reported within the first year of follow-up. We found a significant association between grade ≥ 2 GI toxicity and the parameters GI OAR D0.1cc, D1cc, and D20cc. ConclusionIn this phase II study of MR-guided SABR of oligometastases in the infra-diaphragmatic region, we found a low incidence of toxicity despite half of the lesions being within 10 mm of a radiosensitive OAR. GI OAR D0.1cc, D1cc, and D20cc were associated with grade ≥ 2 GI toxicity.

Using meta-analysis and CNN-NLP to review and classify the medical literature for normal tissue complication probability in head and neck cancer

PurposeThe study aims to enhance the efficiency and accuracy of literature reviews on normal tissue complication probability (NTCP) in head and neck cancer patients using radiation therapy. It employs meta-analysis (MA) and natural language processing (NLP).Material and methodsThe study consists of two parts. First, it employs MA to assess NTCP models for xerostomia, dysphagia, and mucositis after radiation therapy, using Python 3.10.5 for statistical analysis. Second, it integrates NLP with convolutional neural networks (CNN) to optimize literature search, reducing 3256 articles to 12. CNN settings include a batch size of 50, 50–200 epoch range and a 0.001 learning rate.ResultsThe study's CNN-NLP model achieved a notable accuracy of 0.94 after 200 epochs with Adamax optimization. MA showed an AUC of 0.67 for early-effect xerostomia and 0.74 for late-effect, indicating moderate to high predictive accuracy but with high variability across studies. Initial CNN accuracy of 66.70% improved to 94.87% post-tuning by optimizer and hyperparameters.ConclusionThe study successfully merges MA and NLP, confirming high predictive accuracy for specific model-feature combinations. It introduces a time-based metric, words per minute (WPM), for efficiency and highlights the utility of MA and NLP in clinical research.

Normal tissue complication probability model of temporal lobe injury following re-irradiation of IMRT for local recurrent nasopharyngeal carcinoma.

We tried to establish the normal tissue complication probability (NTCP) model of temporal lobe injury of recurrent nasopharyngeal carcinoma (NPC) patients after two courses of intensity modulated radiotherapy (IMRT) to provide more reliable dose-volume data reference to set the temporal lobe tolerance dose for recurrent NPC patients in the future. Recurrent NPC patients were randomly divided into training data set and validation data set in a ratio of 2:1, All the temporal lobes (TLs) were re-contoured as R/L structures and named separately in the MIM system. The dose distribution of the initial IMRT plan was deformed into the second course planning CT via MIM software to get the deformed dose. Equivalent dose of TLs in 2Gy fractions was calculated via linear quadratic model, using an α/β=3 for temporal lobes. NTCP model that correlated the irradiated volume of the temporal lobe and? the clinical variables were evaluated in a multivariate prediction model using AUC analysis. From Jan. 2010 to Dec. 2020, 78 patients were enrolled into our study. Among which 26 (33.3%) developed TLI. The most important factors affecting TLI was the sum-dose d1.5cc of TL, while the possible clinical factors did not reach statistically significant differences in multivariate analysis. According to NTCP model, the TD5 and TD50 EQD2 dose of sum-dose d1.5cc were 65.26Gy (46.72-80.69Gy) and 125.25Gy (89.51-152.18Gy), respectively. For the accumulated EQD2 dose, the area under ROC shadow was 0.8702 (0.7577-0.9828) in model validation, p<0.001. In this study, a NTCP model of temporal lobe injury after a second course of IMRT for recurrent nasopharyngeal carcinoma was established. TD5 and TD50 doses of temporal lobe injury after re-RT were obtained according to the model, and the model was verified by validation set data.

Estimating how contouring differences affect normal tissue complication probability modelling

Background and purposeNormal tissue complication probability (NTCP) models are developed from large retrospective datasets where automatic contouring is often used to contour the organs at risk. This study proposes a methodology to estimate how discrepancies between two sets of contours are reflected on NTCP model performance. We apply this methodology to heart contours within a dataset of non-small cell lung cancer (NSCLC) patients. Materials and methodsOne of the contour sets is designated the ground truth and a dosimetric parameter derived from it is used to simulate outcomes via a predefined NTCP relationship. For each simulated outcome, the selected dosimetric parameters associated with each contour set are individually used to fit a toxicity model and their performance is compared. Our dataset comprised 605 stage IIA-IIIB NSCLC patients. Manual, deep learning, and atlas-based heart contours were available. ResultsHow contour differences were reflected in NTCP model performance depended on the slope of the predefined model, the dosimetric parameter utilized, and the size of the cohort. The impact of contour differences on NTCP model performance increased with steeper NTCP curves. In our dataset, parameters on the low range of the dose-volume histogram were more robust to contour differences. ConclusionsOur methodology can be used to estimate whether a given contouring model is fit for NTCP model development. For the heart in comparable datasets, average Dice should be at least as high as between our manual and deep learning contours for shallow NTCP relationships (88.5 ± 4.5 %) and higher for steep relationships.

Proton Therapy in Breast Cancer: A Review of Potential Approaches for Patient Selection.

Proton radiotherapy may be a compelling technical option for the treatment of breast cancer due to its unique physical property known as the "Bragg peak." This feature offers distinct advantages, promising superior dose conformity within the tumor area and reduced radiation exposure to surrounding healthy tissues, enhancing the potential for better treatment outcomes. However, proton therapy is accompanied by inherent challenges, primarily higher costs and limited accessibility when compared to well-developed photon irradiation. Thus, in clinical practice, it is important for radiation oncologists to carefully select patients before recommendation of proton therapy to ensure the transformation of dosimetric benefits into tangible clinical benefits. Yet, the optimal indications for proton therapy in breast cancer patients remain uncertain. While there is no widely recognized methodology for patient selection, numerous attempts have been made in this direction. In this review, we intended to present an inspiring summarization and discussion about the current practices and exploration on the approaches of this treatment decision-making process in terms of treatment-related side-effects, tumor control, and cost-efficiency, including the normal tissue complication probability (NTCP) model, the tumor control probability (TCP) model, genomic biomarkers, cost-effectiveness analyses (CEAs), and so on. Additionally, we conducted an evaluation of the eligibility criteria in ongoing randomized controlled trials and analyzed their reference value in patient selection. We evaluated the pros and cons of various potential patient selection approaches and proposed possible directions for further optimization and exploration. In summary, while proton therapy holds significant promise in breast cancer treatment, its integration into clinical practice calls for a thoughtful, evidence-driven strategy. By continuously refining the patient selection criteria, we can harness the full potential of proton radiotherapy while ensuring maximum benefit for breast cancer patients.

Dose-response modeling and treatment plan assessment with a python software toolkit

This software assistant aims at calculating the dose-response relations of tumors and normal tissues, or clinically assessing already determined values by other researchers. It can also indicate the optimal dose prescription by optimizing the expected treatment outcome. The software is developed solely in python programming language, and it employs PSFL license for its Graphical User Interface (GUI), NUMPY, MATPLOTLIB, and SCIPY libraries. It comprises of two components. The first is the Dose-response relations derivation component, which takes as input the dose volume histograms (DVHs) of patients and their recorded responses regarding a given clinical endpoint to determine the parameters of different tumor control probability (TCP) or normal tissue complication probability (NTCP) models. The second is the Treatment Plan Assessment component, which uses the DVHs of a plan and the dose-response parameters values of the involved tumors and organs at risk (OARs) to calculate their expected responses. Additionally, the overall probabilities of benefit (PB), injury (PI) and complication-free tumor control (P+) are calculated. The software calculates rapidly the corresponding generalized equivalent uniform doses (gEUD) and biologically effective uniform doses (D‾‾) for the Lyman-Kutcher-Burman (LKB), parallel volume (PV) and relative seriality (RS) models respectively, determining the model parameters. In the Dose-Response Relations Derivation component, the software plots the dose-response curves of the irradiated organ with the relevant confidence internals along with the data of the patients with and without toxicity. It also calculates the odds ratio (OR) and the area under the curve (AUC) of different dose metrics or model parameter values against the individual patient outcomes to determine their discrimination capacity. It also performs a goodness-of-fit evaluation of any model parameter set. The user has the option of viewing plots like Scatter, 3D surfaces, and Bootstrap plots. In the Treatment Plan Assessment part, the software calculates the TCP and NTCP values of the involved tumors and OARs, respectively. Furthermore, it plots the dose-response curves of the TCPs, NTCPs, PB, PI, and P+ for a range of prescription doses for different treatment plans. The presented software is ideal for efficiently conducting studies of radiobiological modeling. Furthermore, it is ideal for performing treatment plan assessment, comparison, and optimization studies.

Proton beam therapy in paediatric cancer: Anticipating the opening of the Australian Bragg Centre for Proton Therapy and Research.

Proton Beam Therapy (PBT) has the potential to improve paediatric cancer care by reducing radiation exposure and thus long-term toxicities. Ethical concerns and debates surrounding the treatment, such as eligibility and accessibility, are ongoing in Australia. The Australian Bragg Centre for Proton Therapy and Research (ABCPTR) (named after Sir William Henry Bragg who described the Bragg peak in his laboratory at the University of Adelaide in 1903) aims to increase access to PBT in Australasia and offer a patient-centred care approach. Research is underway to assess PBT's safety and cost-effectiveness, using tools including Normal Tissue Complication Probability (NTCP) models. Collaborative efforts are focused on developing tailored survivorship clinics to enhance patient follow-up and quality of life. With the anticipated opening of the ABCPTR, Australia is preparing to take a significant step in radiation oncology, offering new research opportunities and creating a publicly funded treatment centre. The initiative aims to balance treatment efficacy with patient care, setting the stage for a future in which radiation therapy will reduce long-term side effects compared to the current standard of care. The implementation of PBT in Australia represents a complex and promising approach to paediatric oncology. This article provides an overview of the current landscape, highlighting the potential benefits and challenges of a treatment that could redefine the quality of survivorship and contribute to global research and best clinical practice.

Normal Tissue Complication Probability Modeling of Severe Radiation-Induced Lymphopenia Using Blood Dose for Patients With Hepatocellular Carcinoma

PurposeThis study aims to develop a normal tissue complication probability (NTCP) model to estimate the risk of severe radiation-induced lymphopenia (SRIL, absolute lymphocyte count (ALC) < 500/μL) by using the blood dose of hepatocellular carcinoma (HCC) patients. Methods and MaterialsWe retrospectively collected data from 75 HCC patients who received radiotherapy (RT) between 2015 and 2018. The hematological dose framework calculated blood dose-volume histograms (DVHs) using a predefined blood flow model, organ DVHs, the number of treatment fractions, and beam delivery time. A Lyman-Kutcher-Burman (LKB) model with a generalized equivalent dose was used to establish the NTCP model, reflecting the whole blood DVHs. The optimization of the LKB parameters was conducted by minimizing a negative log-likelihood function. ResultsThere were 6, 4, 18, 33, and 14 patients in the groups with radiation-induced lymphopenia grades 0, 1, 2, 3, and 4, respectively. The median pre/post-RT ALC were 1410/μL (range, 520–3710/μL) and 470/μL (range, 60–1760/μL), respectively. There was a correlation between blood mean dose and ALC depletion (Pearson r = -0.664; p < 0.001). The average blood mean doses in each radiation-induced lymphopenia group were 2.90 Gy (95% CI, 1.96–3.85 Gy) for grade 0–1, 5.29 Gy (95% CI, 4.12–6.45 Gy) for grade 2, 8.81 Gy (95% CI, 7.55–10.07 Gy) for grade 3, and 11.69 Gy (95% CI, 9.82–17.57 Gy) for grade 4. When applying developed an NTCP model to predict SRIL, the area under the receiver operating characteristic curve and Brier score were 0.89 and 0.12, respectively. ConclusionWe developed the first NTCP model based on whole blood DVHs for estimating SRIL after abdominal RT in HCC patients. Our results showed a strong correlation between blood dose and ALC depletion, suggesting the potential to predict the risk of SRIL occurrence for using blood dose.

The Role of Thoracic Vertebral Body Dosimetry in Minimizing Acute Hematologic Toxicities of Patients With Non-Small Cell Lung Cancer Receiving Lung Radiation Therapy and Immunotherapy

IntroductionHematologic toxicities (HTs) are among the most common toxicities of combined immunotherapy and radiotherapy (RT). It remains essential to prevent RT-induced HTs because they can cause treatment discontinuation (influencing antitumoral effects), and because lymphopenia might dampen the effects of immunotherapy. To date, there are no studies examining the impact of thoracic vertebral body (TVB) RT dose on HTs in non-small cell lung cancer (NSCLC) patients receiving combined lung RT and PD(L)-1 immunotherapy. MethodsFor standardization, all doses were reported as 2-Gy equivalents (EQD2). Mirroring publications prior to the immunotherapy era, TVB volumes referred to T1-T10 and specific dosimetric parameters (DmeanEQD2, V5EQD2‒V60EQD2) were analyzed. Logistic regression estimated associations between grade ≥3 HTs (HT3+) and dosimetric/clinical parameters. Normal tissue complication probability (NTCP) models were constructed by logistic regression analysis modeling for HT3+. Receiver operating characteristic (ROC) analysis delineated TVB dosimetric thresholds, stratification of which was able to evaluate post-RT absolute lymphocyte count (ALC) and immunotherapy responses. Areas under the curve (AUCs) for NTCP models were corroborated by bootstrapping (optimism-corrected) methodology. ResultsIn 132 patients, there were 26 (19.7%) instances of HT3+. On multivariate analysis, DmeanEQD2 and V5EQD2‒V20EQD2 associated with HT3+ (P<0.05 for all). The NTCP models illustrated a 50% probability of HT3+ at a DmeanEQD2=39.8 Gy, V5EQD2=87.4%, V10EQD2=77.0%, and V20EQD2=68.4%. ROC analysis delineated optimal thresholds of HT3+ with DmeanEQD2 ±30.2 Gy, V5EQD2 ±69.1%, V10EQD2 ±64.6%, and V20EQD2 ±53.5%. Patients treated with values above those cutoffs had over double the risk of HT3+, with significant differences in post-RT ALC and immunotherapy responses (P<0.05 for all). AUCs for each individual parameter ranged from 0.743 to 0.798, and combining all four aforementioned cutoffs into a ROC curve resulted in a qualitatively higher AUC (0.836). ConclusionsThis hypothesis-generating work suggests that TVB dosimetry may equate with HT3+ in NSCLC patients undergoing combined lung RT/immunotherapy. Applying TVB dose constraints in this population could reduce HT3+ and avoid dampening of immunotherapy responses, but prospective validation is required.

A multivariable normal tissue complication probability model for predicting radiation-induced hypothyroidism in nasopharyngeal carcinoma patients in the modern radiotherapy era.

Radiation-induced hypothyroidism (RHT) is a common long-term complication for nasopharyngeal carcinoma (NPC) survivors. A model using clinical and dosimetric factors for predicting risk of RHT could suggest a proper dose-volume parameters for the treatment planning in an individual level. We aim to develop a multivariable normal tissue complication probability (NTCP) model for RHT in NPC patients after intensity-modulated radiotherapy or volumetric modulated arc therapy. The model was developed using retrospective clinical data and dose-volume data of the thyroid and pituitary gland based on a standard backward stepwise multivariable logistic regression analysis and was then internally validated using 10-fold cross-validation. The final NTCP model consisted of age, pretreatment thyroid-stimulating hormone and mean thyroid dose. The model performance was good with an area under the receiver operating characteristic curve of 0.749 on an internal (200 patients) and 0.812 on an external (25 patients) validation. The mean thyroid dose at ≤45Gy was suggested for treatment plan, owing to an RHT incidence of 2% versus 61% in the >45Gy group.

NTCP modelling of xerostomia after radiotherapy for oropharyngeal cancer using the PRO-CTCAE and CTCAE scoring systems at different time-points post-RT

PurposeThis study aims at determining the parameter values of three normal tissue complication probability (NTCP) models for the contralateral parotid gland, contralateral submandibular gland (SMG) and contralateral salivary glands regarding the endpoint of xerostomia 6–24 months after radiotherapy for oropharynx cancer. MethodsThe treatment and outcome data of 231 patients with favorable risk, HPV-associated oropharyngeal squamous cell carcinoma are analyzed. 60 Gy intensity modulated radiotherapy was delivered to all the patients. The presence and severity of xerostomia was recorded (pre- and post- radiotherapy) by the PRO-CTCAE and the CTCAE scoring systems. In both scoring systems, patients with a change in symptom severity (from baseline) of ≥ 2 were considered responders. ResultsXerostomia was observed in 61.3 %, 39.2 %, 28.6 % and 27.0 % of the patients based on the PRO-CTCAE scoring system at 6-, 12-, 18- and 24-months post-RT, respectively. The AUCs of the contralateral salivary glands ranged between 0.58–0.64 in the LKB model with the gEUD ranging between 20.3 Gy and 24.7 Gy. ConclusionsBased on the PRO-CTCAE scores, mean dose < 22 Gy, V50 < 10 % for the contralateral salivary glands and mean dose < 18 Gy, V45 < 10 % for the contralateral parotid were found to significantly reduce by a factor of 2–3 the risk for radiation induced xerostomia that is observed at 6–24 months post-RT, respectively. Also, gEUD < 22 Gy to the contralateral salivary glands and < 18 Gy to the contralateral parotid was found to significantly reduce the risk for radiation induced xerostomia that is observed at 6–24 months post-RT by 2.0–2.3 times.

Radiomics and dosiomics for predicting radiation-induced hypothyroidism and guiding intensity-modulated radiotherapy

To guide individualized intensity-modulated radiotherapy (IMRT), we developed and prospectively validated a multiview radiomics risk model for predicting radiation-induced hypothyroidism in patients with nasopharyngeal carcinoma. And simulated radiotherapy plans with same dose-volume-histogram (DVH) but different dose distributions were redesigned to explore the clinical application of the multiview radiomics risk model. The radiomics and dosiomics were built based on selected radiomics and dosiomics features from planning computed tomography and dose distribution, respectively. The multiview radiomics risk model that integrated radiomics, dosiomics, DVH parameters, and clinical factors had better performance than traditional normal tissue complication probability models. And multiview radiomics risk model could identify differences of patient hypothyroidism-free survival that cannot be stratified by traditional models. Besides, two redesigned simulated plans further verified the clinical application and advantage of the multiview radiomics risk model. The multiview radiomics risk model was a promising method to predict radiation-induced hypothyroidism and guide individualized IMRT.

Dosimetric and radiobiological analyses of a de-escalation strategy for elective nodal regions in human papillomavirus-associated oropharyngeal cancer

IntroductionIn this simulation study, we examined the effects of a de-escalation strategy with a reduced dose to subclinical nodal regions in patients with human papillomavirus (HPV)-associated oropharyngeal carcinoma (OPC). MethodsWe created two patterns of intensity-modulated radiotherapy for 16 patients with HPV-associated OPC. In the standard and de-escalation plans, the initial field including elective nodal regions received 46 and 30 Gy, followed by 20 and 36 Gy to the cutdown field, respectively. Comparison metrics were set for each organ at risk (OAR). We compared these metric values and the probability of adverse effects based on the normal tissue complication probability (NTCP) model between the two plans. ResultsBoth plans generally met the dose constraints for the targets and all OAR. Among the comparison metrics, the mean doses to the brain, pharyngeal constrictor muscle, thyroid, and skin and the dose to a 1 % volume of the skin were higher in the standard plan than in the de-escalation plan (P = 0.031, 0.007, < 0.001, < 0.001, and 0.006, respectively). NTCP analyses revealed that the probability of adverse effects in the ipsilateral parotid gland and thyroid was higher in the standard plan than in the de-escalation plan (standard vs. de-escalation plans: ipsilateral parotid gland, 6.4 % vs. 5.0 %, P = 0.016; thyroid, 3.3 % vs. 0.5 %, P < 0.001). ConclusionsA de-escalation strategy with elective nodal regions is a promising treatment to prevent a decline in the quality of life in patients with HPV-associated OPC, particularly xerostomia, dysphagia, and hypothyroidism.

Model development of dose and volume predictors for esophagitis induced during chemoradiotherapy for lung cancer as a step towards radiobiological treatment planning

BackgroundCurrently, radiation therapy treatment planning system intends biological optimization that relies heavily upon plan metrics from tumor control probability (TCP) and normal tissue complication probability (NTCP) modeling. Implementation and expansion of TCP and NTCP models with alternative data is an important step towards reliable radiobiological treatment planning. In this retrospective single institution study, the treatment charts of 139 lung cancer patients treated with chemo-radiotherapy were reviewed and correlated dosimetric predictors with the incidence of esophagitis and established NTCP model of esophagitis grade 1 and 2 for lung cancer patients.MethodsEsophagus is an organ at risk (OAR) in lung cancer radiotherapy (RT). Esophagitis is a common toxicity induced by RT. In this study, dose volume parameters Vx (Vx: percentage esophageal volume receiving ≥ x Gy) and mean esophagus dose (MED) as quantitative dose-volume metrics, the esophagitis grade 1 and 2 as endpoints, were reviewed and derived from the treatment planning system and the electronic medical record system. Statistical analysis of binary logistic regression and probit were performed to have correlated the probability of grade 1 and 2 esophagitis to MED and Vx. IBM SPSS software version 24 at 5% significant level (α = 0.05) was used in the statistical analysis.ResultsThe probabilities of incidence of grade 1 and 2 esophagitis proportionally increased with increasing the values of Vx and MED. V20, V30, V40, V50 and MED are statistically significant good dosimetric predictors of esophagitis grade 1. 50% incidence probability (TD50) of MED for grade 1 and 2 esophagitis were determined. Lyman Kutcher Burman model parameters, such as, n, m and TD50, were fitted and compared with other published findings. Furthermore, the sigmoid shaped dose responding curve between probability of esophagitis grade 1 and MED were generated respecting to races, gender, age and smoking status.ConclusionsV20, V30, V40 and V50 were added onto Quantitative Analysis of Normal Tissue Effects in the clinic, or QUANTEC group’s dose constrains of V35, V50, V70 and MED. Our findings may be useful as both validation of 3-Dimensional planning era models and also additional clinical guidelines in treatment planning and plan evaluation using radiobiology optimization.

Effects of Radiation Therapy on the Female Reproductive Tract in Childhood Cancer Survivors: A PENTEC Comprehensive Review

The PENTEC (Pediatric Normal Tissue Effects in the Clinic) task force aimed to quantify effects of radiation therapy (RT) dose to the female reproductive organs after treatment for childhood cancer. Relevant studies published 1970 to 2017 were identified systematically through PubMed, Medline, and Cochrane databases with additional articles before 2021 identified by the group. Two large studies reported sufficient data to allow modeling of acute ovarian failure (AOF; loss of function ≤5 year from diagnosis) and premature ovarian insufficiency (POI; loss of function at attained age <40 years) based on maximum dose to least affected ovary. Although normal tissue complication probability modeling was not feasible for the uterus due to limited data, the relationship between ultrasound-measured uterine volume and estimated amount of RT was plotted. Limited data regarding vaginal toxicity were available. The risk of AOF increases with RT dose to least affected ovary, alkylating agent cumulative dose (cyclophosphamide equivalent dose [CED] in g/m2), age at RT, and stem cell transplantation: Two Gy to the least affected ovary resulted in AOF risk of 1% to 5% (CED = 0, risk increasing with age), 4% to 7% (CED = 10 g/m2, risk increasing with age), and 6% to 13% (CED = 30 g/m2, risk increasing with age). For patients aged 1 and 20 years at time of RT, AOF risk was ≥50% at doses of 24 Gy and 20 Gy with no alkylating chemotherapy, 22.5 Gy and 17 Gy with intermediate alkylator dose (10 g/m2), and 17 Gy and 13 Gy with high alkylator dose (30 g/m2). Risk of POI increases with survivor (attained) age (rather than age at time of RT), radiation dose to least affected ovary, and alkylator dose. Data review suggested that higher radiation doses to the uterus are associated with uterine toxicity, with uterine size considerably restricted after 12 Gy. Vaginal radiation in children is associated with high toxicity risk, although dose-volume data are not available for quantification. Risk of AOF increases with age at RT, CED exposure, and RT dose; risk of POI likewise increases with RT dose, CED exposure, and survivor age. Both AOF and POI are expected to affect fertility and estrogen production. Data suggest that RT uterine dose >12 Gy may be associated with uterine size restriction. Adult literature suggests that maintaining vaginal dose <5 Gy may limit toxicity. Treatment of life-threatening malignancy remains a priority over reproductive preservation; however, when possible, radiation and surgical techniques should be considered to minimize dose to least affected ovary, uterus, and vagina. Survivors should receive endocrine and gynecologic support; those desiring pregnancy should be counseled early to maximize reproductive options.

External Validation of 8 Normal Tissue Complication Probability Models for Radiation-Associated Primary Hypothyroidism Using Long-Term Thyroid Function Outcomes

Primary hypothyroidism is a common late complication following neck irradiation. The reported incidences were as high as 20-50% among head and neck cancer survivors. Several normal tissue complication probability models for radiation-associated primary hypothyroidism have been proposed, but none has been externally validated using long-term thyroid function outcome. We undertook a systematic review of the literature and identified 8 distinct normal tissue complication probability (NTCP) models for radiation-associated primary hypothyroidism (6 logistic regression models and 2 nomograms). We retrospectively contoured and reviewed the thyroid gland dosimetry of all consecutive head and neck cancer patients who underwent definitive radiotherapy between January 2013 and December 2015 in two tertiary oncology centers. Patients with abnormal pre-radiotherapy thyroid function were excluded. Thyroid function was monitored at least yearly after radiotherapy. Primary hypothyroidism was defined as serum thyroid-stimulating hormone levels above the upper limit of normal. Euthyroid controls were defined as patients who had normal thyroid function at least 5 years beyond the completion of radiotherapy. Diagnostic properties (sensitivity, specificity and accuracy) of the 8 NTCP models were determined. Discrimination and calibration of the models were assessed using c-statistics, discrimination slopes, calibration curves and Brier score. With a median follow-up of 7.0 years, we identified 216 cases of post-radiation primary hypothyroidism and 120 euthyroid controls within the study period. All four logistic regression models (Boomsma 2012, Rønjom 2013, Shen 2021 and Huang 2021) which were based on mean thyroid dose and thyroid volume demonstrated satisfactory diagnostic performances. The accuracy, c-statistic, discrimination slope and Brier score of these four models were > 0.70 (0.711 - 0.732), > 0.75 (0.763 - 0.765), > 0.20 (0.203 - 0.295) and < 0.20 (0.187 - 0.198), respectively. Calibrations of prediction were excellent at ranges of high predicted probabilities. At ranges of low predicted probabilities, these models underestimated the risk of primary hypothyroidism by 10-20%. Other models or nomograms which utilized alternative radiation dose-volume parameters (e.g., V30, V35 or V50), or without factoring in the variation in thyroid volume, had significantly worse performance in terms of both discrimination and calibration. This large external validation study showed that NTCP models that included thyroid volume and mean thyroid dose had satisfactory performances in predicting long-term thyroid function outcomes following neck irradiation. Refinement is required to optimize calibration in patients deemed at low risk of post-radiation primary hypothyroidism.