Large-scale meta-genome-wide association study reveals common genetic factors linked to radiation-induced acute toxicities across cancer types.

BACKGROUND: This study was designed to identify common genetic susceptibility and shared genetic variants associated with acute radiation-induced toxicity across 4 cancer types (prostate, head and neck, breast, and lung). METHODS: A genome-wide association study meta-analysis was performed using 19 cohorts totaling 12 042 patients. Acute standardized total average toxicity (STATacute) was modelled using a generalized linear regression model for additive effect of genetic variants, adjusted for demographic and clinical covariates (rSTATacute). Linkage disequilibrium score regression estimated shared single-nucleotide variation (SNV-formerly SNP)-based heritability of rSTATacute in all patients and for each cancer type. RESULTS: Shared SNV-based heritability of STATacute among all cancer types was estimated at 10% (SE = 0.02) and was higher for prostate (17%, SE = 0.07), head and neck (27%, SE = 0.09), and breast (16%, SE = 0.09) cancers. We identified 130 suggestive associated SNVs with rSTATacute (5.0 × 10‒8 < P < 1.0 × 10‒5) across 25 genomic regions. rs142667902 showed the strongest association (effect allele A; effect size ‒0.17; P = 1.7 × 10‒7), which is located near DPPA4, encoding a protein involved in pluripotency in stem cells, which are essential for repair of radiation-induced tissue injury. Gene-set enrichment analysis identified 'RNA splicing via endonucleolytic cleavage and ligation' (P = 5.1 × 10‒6, P = .079 corrected) as the top gene set associated with rSTATacute among all patients. In silico gene expression analysis showed that the genes associated with rSTATacute were statistically significantly up-regulated in skin (not sun exposed P = .004 corrected; sun exposed P = .026 corrected). CONCLUSIONS: There is shared SNV-based heritability for acute radiation-induced toxicity across and within individual cancer sites. Future meta-genome-wide association studies among large radiation therapy patient cohorts are worthwhile to identify the common causal variants for acute radiotoxicity across cancer types.

Sub-regional analysis of the parotid glands: model development for predicting late xerostomia with radiomics features in head and neck cancer patients.

BACKGROUND: The irradiation of sub-regions of the parotid has been linked to xerostomia development in patients with head and neck cancer (HNC). In this study, we compared the xerostomia classification performance of radiomics features calculated on clinically relevant and de novo sub-regions of the parotid glands of HNC patients. MATERIAL AND METHODS: All patients (N = 117) were treated with TomoTherapy in 30-35 fractions of 2-2.167 Gy per fraction with daily mega-voltage-CT (MVCT) acquisition for image-guidance purposes. Radiomics features (N = 123) were extracted from daily MVCTs for the whole parotid gland and nine sub-regions. The changes in feature values after each complete week of treatment were considered as predictors of xerostomia (CTCAEv4.03, grade ≥ 2) at 6 and 12 months. Combinations of predictors were generated following the removal of statistically redundant information and stepwise selection. The classification performance of the logistic regression models was evaluated on train and test sets of patients using the Area Under the Curve (AUC) associated with the different sub-regions at each week of treatment and benchmarked with the performance of models solely using dose and toxicity at baseline. RESULTS: In this study, radiomics-based models predicted xerostomia better than standard clinical predictors. Models combining dose to the parotid and xerostomia scores at baseline yielded an AUCtest of 0.63 and 0.61 for xerostomia prediction at 6 and 12 months after radiotherapy while models based on radiomics features extracted from the whole parotid yielded a maximum AUCtest of 0.67 and 0.75, respectively. Overall, across sub-regions, maximum AUCtest was 0.76 and 0.80 for xerostomia prediction at 6 and 12 months. Within the first two weeks of treatment, the cranial part of the parotid systematically yielded the highest AUCtest. CONCLUSION: Our results indicate that variations of radiomics features calculated on sub-regions of the parotid glands can lead to earlier and improved prediction of xerostomia in HNC patients.

A Cross-Sectional Study of Commercial Ewe Management Practices for Different Sheep Breeds across Southern Australia.

The management of ewes across southern Australia may vary with breed and can change over time and, as such, a greater understanding of producer management practices and the motivations that influence these practices is required. A cross-sectional study was performed by telephone interview with sheep producers managing Composite, Maternal, Merino or shedding ewe breeds mated in either spring, summer, or autumn. The surveyed producers were a unique subset of southern Australian producers. A large proportion of the surveyed producers followed current best practice guidelines for ewe mating and lambing nutritional management; however, some producers did not align with these targets. Further, some producers did not see the value in attaining the current recommendations. Pregnancy scanning was widely practiced, likely an artefact of the recruitment process; however, a few producers did not utilize this information for nutritional management at lambing time. Finally, most producers were active in their search for new information, seeking information regularly from a wide range of sources and reported making management changes within the last five years. Further work is required to understand why some producers are not adopting best practice where possible and to understand current barriers for adoption. Management guidelines for all sheep breeds are required to best manage sheep across southern Australia.

Assessing the generalisability of radiomics features previously identified as predictive of radiation-induced sticky saliva and xerostomia.

Background and purpose: While core to the scientific approach, reproducibility of experimental results is challenging in radiomics studies. A recent publication identified radiomics features that are predictive of late irradiation-induced toxicity in head and neck cancer (HNC) patients. In this study, we assessed the generalisability of these findings. Materials and Methods: The procedure described in the publication in question was applied to a cohort of 109 HNC patients treated with 50-70 Gy in 20-35 fractions using helical radiotherapy although there were inherent differences between the two patient populations and methodologies. On each slice of the planning CT with delineated parotid and submandibular glands, the imaging features that were previously identified as predictive of moderate-to-severe xerostomia and sticky saliva 12 months post radiotherapy (Xer12m and SS12m) were calculated. Specifically, Short Run Emphasis (SRE) and maximum CT intensity (maxHU) were evaluated for improvement in prediction of Xer12m and SS12m respectively, compared to models solely using baseline toxicity and mean dose to the salivary glands. Results: None of the associations previously identified as statistically significant and involving radiomics features in univariate or multivariate models could be reproduced on our cohort. Conclusion: The discrepancies observed between the results of the two studies delineate limits to the generalisability of the previously reported findings. This may be explained by the differences in the approaches, in particular the imaging characteristics and subsequent methodological implementation. This highlights the importance of external validation, high quality reporting guidelines and standardisation protocols to ensure generalisability, replication and ultimately clinical implementation.

Predicting radiotherapy-induced xerostomia in head and neck cancer patients using day-to-day kinetics of radiomics features.

Background and purpose: The images acquired during radiotherapy for image-guidance purposes could be used to monitor patient-specific response to irradiation and improve treatment personalisation. We investigated whether the kinetics of radiomics features from daily mega-voltage CT image-guidance scans (MVCT) improve prediction of moderate-to-severe xerostomia compared to dose/volume parameters in radiotherapy of head-and-neck cancer (HNC). Materials and Methods: All included HNC patients (N = 117) received 30 or more fractions of radiotherapy with daily MVCTs. Radiomics features were calculated on the contra-lateral parotid glands of daily MVCTs. Their variations over time after each complete week of treatment were used to predict moderate-to-severe xerostomia (CTCAEv4.03 grade ≥ 2) at 6, 12 and 24 months post-radiotherapy. After dimensionality reduction, backward/forward selection was used to generate combinations of predictors.Three types of logistic regression model were generated for each follow-up time: 1) a pre-treatment reference model using dose/volume parameters, 2) a combination of dose/volume and radiomics-based predictors, and 3) radiomics-based predictors. The models were internally validated by cross-validation and bootstrapping and their performance evaluated using Area Under the Curve (AUC) on separate training and testing sets. Results: Moderate-to-severe xerostomia was reported by 46 %, 33 % and 26 % of the patients at 6, 12 and 24 months respectively. The selected models using radiomics-based features extracted at or before mid-treatment outperformed the dose-based models with an AUCtrain/AUCtest of 0.70/0.69, 0.76/0.74, 0.86/0.86 at 6, 12 and 24 months, respectively. Conclusion: Our results suggest that radiomics features calculated on MVCTs from the first half of the radiotherapy course improve prediction of moderate-to-severe xerostomia in HNC patients compared to a dose-based pre-treatment model.

Mating Conditions and Management Practices Influence Pregnancy Scanning Outcomes Differently between Ewe Breeds.

Sheep production in southern Australia may vary by breed, time of year, production output (wool, meat, or both), region and seasonal influence. Sheep producers with flocks of approximately 300-500 ewes (n = 58) were recruited across southern Australia to take part in a survey and mating variables were collected from over 30,000 ewes between October 2020 and August 2021. A Bayesian Network (BN) was developed to identify the interrelatedness and most influential variable on pregnancy and fetal number (of pregnant ewes) outcomes under different scenarios. The BN analysis indicated a low association between the variables explored, however, were breed dependent. In wool-based breeds a mating liveweight of 60-69.5 kg predicted the lowest non-pregnant and greatest number of fetuses, and in shedding ewes 70-79.5 kg predicted the lowest non-pregnant rate and 90-99.5 kg the greatest number of fetuses. Pregnancy rate and fetuses per ewe were optimized at ram percentages of 1.5% for Composite and Merino ewes and 2% for Maternal ewes. A mating BCS 4 resulted in greatest pregnancy rate and number of fetuses across all breeds. Curvilinear relationships between mating liveweight, BCS and ram percentage were observed with pregnancy rate and fetal number. Practically, reproductive potential is best managed on a breed basis and with consideration of all variables explored.

Meta-GWAS identifies the heritability of acute radiation-induced toxicities in head and neck cancer.

BACKGROUND AND PURPOSE: We aimed to the genetic components and susceptibility variants associated with acute radiation-induced toxicities (RITs) in patients with head and neck cancer (HNC). MATERIALS AND METHODS: We performed the largest meta-GWAS of seven European cohorts (n = 4,042). Patients were scored weekly during radiotherapy for acute RITs including dysphagia, mucositis, and xerostomia. We analyzed the effect of variants on the average burden (measured as area under curve, AUC) per each RIT, and standardized total average acute toxicity (STATacute) score using a multivariate linear regression. We tested suggestive variants (p < 1.0x10-5) in discovery set (three cohorts; n = 2,640) in a replication set (four cohorts; n = 1,402). We meta-analysed all cohorts to calculate RITs specific SNP-based heritability, and effect of polygenic risk scores (PRSs), and genetic correlations among RITS. RESULTS: From 393 suggestive SNPs identified in discovery set; 37 were nominally significant (preplication < 0.05) in replication set, but none reached genome-wide significance (pcombined < 5 × 10-8). In-silico functional analyses identified "3'-5'-exoribonuclease activity" (FDR = 1.6e-10) for dysphagia, "inositol phosphate-mediated signalling" for mucositis (FDR = 2.20e-09), and "drug catabolic process" for STATacute (FDR = 3.57e-12) as the most enriched pathways by the RIT specific suggestive genes. The SNP-based heritability (±standard error) was 29 ± 0.08 % for dysphagia, 9 ± 0.12 % (mucositis) and 27 ± 0.09 % (STATacute). Positive genetic correlation was rg = 0.65 (p = 0.048) between dysphagia and STATacute. PRSs explained limited variation of dysphagia (3 %), mucositis (2.5 %), and STATacute (0.4 %). CONCLUSION: In HNC patients, acute RITs are modestly heritable, sharing 10 % genetic susceptibility, when PRS explains < 3 % of their variance. We identified numerus suggestive SNPs, which remain to be replicated in larger studies.

Associations between voxel-level accumulated dose and rectal toxicity in prostate radiotherapy.

BACKGROUND AND PURPOSE: Associations between dose and rectal toxicity in prostate radiotherapy are generally poorly understood. Evaluating spatial dose distributions to the rectal wall (RW) may lead to improvements in dose-toxicity modelling by incorporating geometric information, masked by dose-volume histograms. Furthermore, predictive power may be strengthened by incorporating the effects of interfraction motion into delivered dose calculations.Here we interrogate 3D dose distributions for patients with and without toxicity to identify rectal subregions at risk (SRR), and compare the discriminatory ability of planned and delivered dose. MATERIAL AND METHODS: Daily delivered dose to the rectum was calculated using image guidance scans, and accumulated at the voxel level using biomechanical finite element modelling. SRRs were statistically determined for rectal bleeding, proctitis, faecal incontinence and stool frequency from a training set (n = 139), and tested on a validation set (n = 47). RESULTS: SRR patterns differed per endpoint. Analysing dose to SRRs improved discriminative ability with respect to the full RW for three of four endpoints. Training set AUC and OR analysis produced stronger toxicity associations from accumulated dose than planned dose. For rectal bleeding in particular, accumulated dose to the SRR (AUC 0.76) improved upon dose-toxicity associations derived from planned dose to the RW (AUC 0.63). However, validation results could not be considered significant. CONCLUSIONS: Voxel-level analysis of dose to the RW revealed SRRs associated with rectal toxicity, suggesting non-homogeneous intra-organ radiosensitivity. Incorporating spatial features of accumulated delivered dose improved dose-toxicity associations. This may be an important tool for adaptive radiotherapy in the future.

Anatomical change during radiotherapy for head and neck cancer, and its effect on delivered dose to the spinal cord.

BACKGROUND AND PURPOSE: The impact of weight loss and anatomical change during head and neck (H&N) radiotherapy on spinal cord dosimetry is poorly understood, limiting evidence-based adaptive management strategies. MATERIALS AND METHODS: 133 H&N patients treated with daily mega-voltage CT image-guidance (MVCT-IG) on TomoTherapy, were selected. Elastix software was used to deform planning scan SC contours to MVCT-IG scans, and accumulate dose. Planned (DP) and delivered (DA) spinal cord D2% (SCD2%) were compared. Univariate relationships between neck irradiation strategy (unilateral vs bilateral), T-stage, N-stage, weight loss, and changes in lateral separation (LND) and CT slice surface area (SSA) at C1 and the superior thyroid notch (TN), and ΔSCD2% [(DA - DP) D2%] were examined. RESULTS: The mean value for (DA - DP) D2% was -0.07 Gy (95%CI -0.28 to 0.14, range -5.7 Gy to 3.8 Gy), and the mean absolute difference between DP and DA (independent of difference direction) was 0.9 Gy (95%CI 0.76-1.04 Gy). Neck treatment strategy (p = 0.39) and T-stage (p = 0.56) did not affect ΔSCD2%. Borderline significance (p = 0.09) was seen for higher N-stage (N2-3) and higher ΔSCD2%. Mean reductions in anatomical metrics were substantial: weight loss 6.8 kg; C1LND 12.9 mm; C1SSA 12.1 cm2; TNLND 5.3 mm; TNSSA 11.2 cm2, but no relationship between weight loss or anatomical change and ΔSCD2% was observed (all r2 < 0.1). CONCLUSIONS: Differences between delivered and planned spinal cord D2% are small in patients treated with daily IG. Even patients experiencing substantial weight loss or anatomical change during treatment do not require adaptive replanning for spinal cord safety.

Identifying risk factors for L'Hermitte's sign after IMRT for head and neck cancer.

BACKGROUND: L'Hermitte's sign (LS) after chemoradiotherapy for head and neck cancer appears related to higher spinal cord doses. IMRT plans limit spinal cord dose, but the incidence of LS remains high. METHODS: One hundred seventeen patients treated with TomoTherapy™ between 2008 and 2015 prospectively completed a side-effect questionnaire (VoxTox Trial Registration: UK CRN ID 13716). Baseline patient and treatment data were collected. Radiotherapy plans were analysed; mean and maximum spinal cord dose and volumes receiving 10, 20, 30 and 40 Gy were recorded. Dose variation across the cord was examined. These data were included in a logistic regression model. RESULTS: Forty two patients (35.9%) reported LS symptoms. Concurrent weekly cisplatin did not increase LS risk (p = 0.70, OR = 1.23 {95% CI 0.51-2.34}). Of 13 diabetic participants (9 taking metformin), only 1 developed LS (p = 0.025, OR = 0.13 {95% CI 0.051-3.27}). A refined binary logistic regression model showed that patients receiving unilateral radiation (p = 0.019, OR = 2.06 {95% CI 0.15-0.84}) were more likely to develop LS. Higher V40Gy (p = 0.047, OR = 1.06 {95% CI 1.00-1.12}), and younger age (mean age 56.6 vs 59.7, p = 0.060, OR = 0.96 {95% CI 0.92-1.00}) were associated with elevated risk of LS, with borderline significance. CONCLUSIONS: In this cohort, concomitant cisplatin did not increase risk, and LS incidence was lower in diabetic patients. Patient age and dose gradients across the spinal cord may be important factors.

Applying physical science techniques and CERN technology to an unsolved problem in radiation treatment for cancer: the multidisciplinary 'VoxTox' research programme.

The VoxTox research programme has applied expertise from the physical sciences to the problem of radiotherapy toxicity, bringing together expertise from engineering, mathematics, high energy physics (including the Large Hadron Collider), medical physics and radiation oncology. In our initial cohort of 109 men treated with curative radiotherapy for prostate cancer, daily image guidance computed tomography (CT) scans have been used to calculate delivered dose to the rectum, as distinct from planned dose, using an automated approach. Clinical toxicity data have been collected, allowing us to address the hypothesis that delivered dose provides a better predictor of toxicity than planned dose.

Recalculation of dose for each fraction of treatment on TomoTherapy.

OBJECTIVE: The VoxTox study, linking delivered dose to toxicity requires recalculation of typically 20-37 fractions per patient, for nearly 2000 patients. This requires a non-interactive interface permitting batch calculation with multiple computers. METHODS: Data are extracted from the TomoTherapy(®) archive and processed using the computational task-management system GANGA. Doses are calculated for each fraction of radiotherapy using the daily megavoltage (MV) CT images. The calculated dose cube is saved as a digital imaging and communications in medicine RTDOSE object, which can then be read by utilities that calculate dose-volume histograms or dose surface maps. The rectum is delineated on daily MV images using an implementation of the Chan-Vese algorithm. RESULTS: On a cluster of up to 117 central processing units, dose cubes for all fractions of 151 patients took 12 days to calculate. Outlining the rectum on all slices and fractions on 151 patients took 7 h. We also present results of the Hounsfield unit (HU) calibration of TomoTherapy MV images, measured over an 8-year period, showing that the HU calibration has become less variable over time, with no large changes observed after 2011. CONCLUSION: We have developed a system for automatic dose recalculation of TomoTherapy dose distributions. This does not tie up the clinically needed planning system but can be run on a cluster of independent machines, enabling recalculation of delivered dose without user intervention. ADVANCES IN KNOWLEDGE: The use of a task management system for automation of dose calculation and outlining enables work to be scaled up to the level required for large studies.

Accumulated dose to the rectum, measured using dose-volume histograms and dose-surface maps, is different from planned dose in all patients treated with radiotherapy for prostate cancer.

OBJECTIVE: We sought to calculate accumulated dose (DA) to the rectum in patients treated with radiotherapy for prostate cancer. We were particularly interested in whether dose-surface maps (DSMs) provide additional information to dose-volume histograms (DVHs). METHODS: Manual rectal contours were obtained for kilovoltage and daily megavoltage CT scans for 10 participants from the VoxTox study (380 scans). Daily delivered dose recalculation was performed using a ray-tracing algorithm. Delivered DVHs were summated to create accumulated DVHs. The rectum was considered as a cylinder, cut and unfolded to produce daily delivered DSMs; these were summated to produce accumulated DSMs. RESULTS: Accumulated dose-volumes were different from planned in all participants. For one participant, all DA levels were higher and all volumes were larger than planned. For four participants, all DA levels were lower and all volumes were smaller than planned. For each of these four participants, ≥1% of pixels on the accumulated DSM received ≥5 Gy more than had been planned. CONCLUSION: Differences between accumulated and planned dose-volumes were seen in all participants. DSMs were able to identify differences between DA and planned dose that could not be appreciated from the DVHs. Further work is needed to extract the dose data embedded in the DSMs. These will be correlated with toxicity as part of the VoxTox Programme. ADVANCES IN KNOWLEDGE: DSMs are able to identify differences between DA and planned dose that cannot be appreciated from DVHs alone and should be incorporated into future studies investigating links between DA and toxicity.