Reciprocal interactions between tumour cell populations enhance growth and reduce radiation sensitivity in prostate cancer.

Intratumoural heterogeneity (ITH) contributes to local recurrence following radiotherapy in prostate cancer. Recent studies also show that ecological interactions between heterogeneous tumour cell populations can lead to resistance in chemotherapy. Here, we evaluated whether interactions between heterogenous populations could impact growth and response to radiotherapy in prostate cancer. Using mixed 3D cultures of parental and radioresistant populations from two prostate cancer cell lines and a predator-prey mathematical model to investigate various types of ecological interactions, we show that reciprocal interactions between heterogeneous populations enhance overall growth and reduce radiation sensitivity. The type of interaction influences the time of regrowth after radiation, and, at the population level, alters the survival and cell cycle of each population without eliminating either one. These interactions can arise from oxygen constraints and from cellular cross-talk that alter the tumour microenvironment. These findings suggest that ecological-type interactions are important in radiation response and could be targeted to reduce local recurrence.

Systematic review and meta-analysis of small bowel dose-volume and acute toxicity in conventionally-fractionated rectal cancer radiotherapy.

INTRODUCTION: The limited radiation tolerance of the small-bowel causes toxicity for patients receiving conventionally-fractionated radiotherapy for rectal cancer. Safe radiotherapy dose-escalation will require a better understanding of such toxicity. We conducted a systematic review and meta-analysis using published datasets of small bowel dose-volume and outcomes to analyse the relationship with acute toxicity. MATERIALS AND METHODS: SCOPUS, EMBASE & MEDLINE were searched to identify twelve publications reporting small-bowel dose-volumes and toxicity data or analysis. Where suitable data were available (mean absolute volume with parametric error measures), fixed-effects inverse-variance meta-analysis was used to compare cohorts of patients according to Grade ≥3 toxicity. For other data, non-parametric examinations of irradiated small-bowel dose-volume and incidence of toxicity were conducted, and a univariate logistic regression model was fitted. RESULTS: On fixed-effects meta-analysis of three studies (203 patients), each of the dose-volume measures V5Gy-V40Gy were significantly greater (p < 0.00001) for patients with Grade ≥3 toxicity than for those without. Absolute difference was largest for the lowest dose-volume parameter; however relative difference increases with increasing dose. On logistic regression multiple small-bowel DVH parameters were predictive of toxicity risk (V5Gy, V10Gy, V30Gy - V45Gy), with V10Gy the strongest (p = 0.004). CONCLUSIONS: Analysis of published clinical cohort dose-volume data provides evidence for a significant dose-volume-toxicity response effect for a wide range of clinically-relevant doses in the treatment of rectal cancer. Both low dose and high dose are shown to predict toxicity risk, which has important implications for radiotherapy planning and consideration of dose escalation for these patients.

Functional Parameters Derived from Magnetic Resonance Imaging Reflect Vascular Morphology in Preclinical Tumors and in Human Liver Metastases.

Purpose: Tumor vessels influence the growth and response of tumors to therapy. Imaging vascular changes in vivo using dynamic contrast-enhanced MRI (DCE-MRI) has shown potential to guide clinical decision making for treatment. However, quantitative MR imaging biomarkers of vascular function have not been widely adopted, partly because their relationship to structural changes in vessels remains unclear. We aimed to elucidate the relationships between vessel function and morphology in vivo Experimental Design: Untreated preclinical tumors with different levels of vascularization were imaged sequentially using DCE-MRI and CT. Relationships between functional parameters from MR (iAUC, K trans, and BATfrac) and structural parameters from CT (vessel volume, radius, and tortuosity) were assessed using linear models. Tumors treated with anti-VEGFR2 antibody were then imaged to determine whether antiangiogenic therapy altered these relationships. Finally, functional-structural relationships were measured in 10 patients with liver metastases from colorectal cancer.Results: Functional parameters iAUC and K trans primarily reflected vessel volume in untreated preclinical tumors. The relationships varied spatially and with tumor vascularity, and were altered by antiangiogenic treatment. In human liver metastases, all three structural parameters were linearly correlated with iAUC and K trans For iAUC, structural parameters also modified each other's effect.Conclusions: Our findings suggest that MR imaging biomarkers of vascular function are linked to structural changes in tumor vessels and that antiangiogenic therapy can affect this link. Our work also demonstrates the feasibility of three-dimensional functional-structural validation of MR biomarkers in vivo to improve their biological interpretation and clinical utility. Clin Cancer Res; 24(19); 4694-704. ©2018 AACR.

Potential of Proton Therapy to Reduce Acute Hematologic Toxicity in Concurrent Chemoradiation Therapy for Esophageal Cancer.

PURPOSE: Radiation therapy dose escalation using a simultaneous integrated boost (SIB) is predicted to improve local tumor control in esophageal cancer; however, any increase in acute hematologic toxicity (HT) could limit the predicted improvement in patient outcomes. Proton therapy has been shown to significantly reduce HT in lung cancer patients receiving concurrent chemotherapy. Therefore, we investigated the potential of bone marrow sparing with protons for esophageal tumors. METHODS AND MATERIALS: Twenty-one patients with mid-esophageal cancer who had undergone conformal radiation therapy (3D50) were selected. Two surrogates for bone marrow were created by outlining the thoracic bones (bone) and only the body of the thoracic vertebrae (TV) in Eclipse. The percentage of overlap of the TV with the planning treatment volume was recorded for each patient. Additional plans were created retrospectively, including a volumetric modulated arc therapy (VMAT) plan with the same dose as for 3D50; a VMAT SIB plan with a dose prescription of 62.5 Gy to the high-risk subregion within the planning treatment volume; a reoptimized TV-sparing VMAT plan; and a proton therapy plan with the same SIB dose prescription. The bone and TV dose metrics were recorded and compared across all plans and variations with respect to PTV and percentage of overlap for each patient. RESULTS: The 3D50 plans showed the highest bone mean dose and TV percentage of volume receiving ≥30 Gy (V30Gy) for each patient. The VMAT plans irradiated a larger bone V10Gy than did the 3D50 plans. The reoptimized VMAT62.5 VT plans showed improved sparing of the TV volume, but only the proton plans showed significant sparing for bone V10Gy and bone mean dose, especially for patients with a larger PTV. CONCLUSIONS: The results of the present study have shown that proton therapy can reduced bone marrow toxicity.

[18F]Fluoromisonidazole PET in rectal cancer.

BACKGROUND: There is an increasing interest in developing predictive biomarkers of tissue hypoxia using functional imaging for personalised radiotherapy in patients with rectal cancer that are considered for neoadjuvant chemoradiotherapy (CRT). The study explores [18F]fluoromisonidazole ([18F]FMISO) positron emission tomography (PET) scans for predicting clinical response in rectal cancer patients receiving neoadjuvant CRT. METHODS: Patients with biopsy-proven rectal adenocarcinoma were imaged at 0-45 min, 2 and 4 h, at baseline and after 8-10 fractions of CRT (week 2). The first 6 patients did not receive an enema (the non-enema group) and the last 4 patients received an enema before PET-CT scan (the enema group). [18F]FMISO production failed on 2 occasions. Static PET images at 4 h were analysed using tumour-to-muscle (T:M) SUVmax and tumour-to-blood (T:B) SUVmax. The 0-45 min dynamic PET scans were analysed using Casciari model to report hypoxia and perfusion. Akaike information criteria (AIC) were used to compare data fittings for different pharmacokinetic models. Pathological tumour regression grade was scored using American Joint Committee on Cancer (AJCC) 7.0. Shapiro-Wilk test was used to evaluate the normality of the data. RESULTS: Five out of eleven (5/11) patients were classed as good responders (AJCC 0/1 or good clinical response) and 6/11 as poor responders (AJCC 2/3 or poor clinical response). The median T:M SUVmax was 2.14 (IQR 0.58) at baseline and 1.30 (IQR 0.19) at week 2, and the corresponding median tumour hypoxia volume was 1.08 (IQR 1.31) cm3 and 0 (IQR 0.15) cm3, respectively. The median T:B SUVmax was 2.46 (IQR 1.50) at baseline and 1.61 (IQR 0.14) at week 2, and the corresponding median tumour hypoxia volume was 5.68 (IQR 5.86) cm3 and 0.76 (IQR 0.78) cm3, respectively. For 0-45 min tumour modelling, the median hypoxia was 0.92 (IQR 0.41) min-1 at baseline and 0.70 (IQR 0.10) min-1 at week 2. The median perfusion was 4.10 (IQR 1.71) ml g-1 min-1 at baseline and 2.48 (IQR 3.62) ml g-1 min-1 at week 2. In 9/11 patients with both PET scans, tumour perfusion decreased in non-responders and increased in responders except in one patient. None of the changes in other PET parameters showed any clear trend with clinical outcome. CONCLUSIONS: This pilot study with small number of datasets revealed significant challenges in delivery and interpretation of [18F]FMISO PET scans of rectal cancer. There are two principal problems namely spill-in from non-tumour tracer activity from rectal and bladder contents. Emphasis should be made on reducing spill-in effects from the bladder to improve data quality. This preliminary study has shown fundamental difficulties in the interpretation of [18F]FMISO PET scans for rectal cancer, limiting its clinical applicability.

An approximate analytical solution of the Bethe equation for charged particles in the radiotherapeutic energy range.

Charged particles such as protons and carbon ions are an increasingly important tool in radiotherapy. There are however unresolved physics issues impeding optimal implementation, including estimation of dose deposition in non-homogeneous tissue, an essential aspect of treatment optimization. Monte Carlo (MC) methods can be employed to estimate radiation profile, and whilst powerful, these are computationally expensive, limiting practicality. In this work, we start from fundamental physics in the form of the Bethe equation to yield a novel approximate analytical solution for particle range, energy and linear energy transfer (LET). The solution is given in terms of the exponential integral function with relativistic co-ordinate transform, allowing application at radiotherapeutic energy levels (50-350 MeV protons, 100-600 Mev/a.m.u carbon ions). Model results agreed closely for protons and carbon-ions (mean error within ≈1%) of literature values. Agreement was high along particle track, with some discrepancy manifesting at track-end. The model presented has applications within a charged particle radiotherapy optimization framework as a rapid method for dose and LET estimation, capable of accounting for heterogeneity in electron density and ionization potential.

Modelling duodenum radiotherapy toxicity using cohort dose-volume-histogram data.

BACKGROUND AND PURPOSE: Gastro-intestinal toxicity is dose-limiting in abdominal radiotherapy and correlated with duodenum dose-volume parameters. We aimed to derive updated NTCP model parameters using published data and prospective radiotherapy quality-assured cohort data. MATERIAL AND METHODS: A systematic search identified publications providing duodenum dose-volume histogram (DVH) statistics for clinical studies of conventionally-fractionated radiotherapy. Values for the Lyman-Kutcher-Burman (LKB) NTCP model were derived through sum-squared-error minimisation and using leave-one-out cross-validation. Data were corrected for fraction size and weighted according to patient numbers, and the model refined using individual patient DVH data for two further cohorts from prospective clinical trials. RESULTS: Six studies with published DVH data were utilised, and with individual patient data included outcomes for 531 patients in total (median follow-up 16months). Observed gastro-intestinal toxicity rates ranged from 0% to 14% (median 8%). LKB parameter values for unconstrained fit to published data were: n=0.070, m=0.46, TD50(1) [Gy]=183.8, while the values for the model incorporating the individual patient data were n=0.193, m=0.51, TD50(1) [Gy]=299.1. CONCLUSIONS: LKB parameters derived using published data are shown to be consistent to those previously obtained using individual patient data, supporting a small volume-effect and dependence on exposure to high threshold dose.

The role of necrosis, acute hypoxia and chronic hypoxia in 18F-FMISO PET image contrast: a computational modelling study.

Positron emission tomography (PET) using 18F-fluoromisonidazole (FMISO) is a promising technique for imaging tumour hypoxia, and a potential target for radiotherapy dose-painting. However, the relationship between FMISO uptake and oxygen partial pressure ([Formula: see text]) is yet to be quantified fully. Tissue oxygenation varies over distances much smaller than clinical PET resolution (<100 µm versus ∼4 mm), and cyclic variations in tumour perfusion have been observed on timescales shorter than typical FMISO PET studies (∼20 min versus a few hours). Furthermore, tracer uptake may be decreased in voxels containing some degree of necrosis. This work develops a computational model of FMISO uptake in millimetre-scale tumour regions. Coupled partial differential equations govern the evolution of oxygen and FMISO distributions, and a dynamic vascular source map represents temporal variations in perfusion. Local FMISO binding capacity is modulated by the necrotic fraction. Outputs include spatiotemporal maps of [Formula: see text] and tracer accumulation, enabling calculation of tissue-to-blood ratios (TBRs) and time-activity curves (TACs) as a function of mean tissue oxygenation. The model is characterised using experimental data, finding half-maximal FMISO binding at local [Formula: see text] of 1.4 mmHg (95% CI: 0.3-2.6 mmHg) and half-maximal necrosis at 1.2 mmHg (0.1-4.9 mmHg). Simulations predict a non-linear non-monotonic relationship between FMISO activity (4 hr post-injection) and mean tissue [Formula: see text] : tracer uptake rises sharply from negligible levels in avascular tissue, peaking at ∼5 mmHg and declining towards blood activity in well-oxygenated conditions. Greater temporal variation in perfusion increases peak TBRs (range 2.20-5.27) as a result of smaller predicted necrotic fraction, rather than fundamental differences in FMISO accumulation under acute hypoxia. Identical late FMISO uptake can occur in regions with differing [Formula: see text] and necrotic fraction, but simulated TACs indicate that additional early-phase information may allow discrimination of hypoxic and necrotic signals. We conclude that a robust approach to FMISO interpretation (and dose-painting prescription) is likely to be based on dynamic PET analysis.

Dose and Fractionation in Radiation Therapy of Curative Intent for Non-Small Cell Lung Cancer: Meta-Analysis of Randomized Trials.

PURPOSE: The optimum dose and fractionation in radiation therapy of curative intent for non-small cell lung cancer remains uncertain. We undertook a published data meta-analysis of randomized trials to examine whether radiation therapy regimens with higher time-corrected biologically equivalent doses resulted in longer survival, either when given alone or when given with chemotherapy. METHODS AND MATERIALS: Eligible studies were randomized comparisons of 2 or more radiation therapy regimens, with other treatments identical. Median survival ratios were calculated for each comparison and pooled. RESULTS: 3795 patients in 25 randomized comparisons of radiation therapy dose were studied. The median survival ratio, higher versus lower corrected dose, was 1.13 (95% confidence interval [CI] 1.04-1.22) when radiation therapy was given alone and 0.83 (95% CI 0.71-0.97) when it was given with concurrent chemotherapy (P for difference=.001). In comparisons of radiation therapy given alone, the survival benefit increased with increasing dose difference between randomized treatment arms (P for trend=.004). The benefit increased with increasing dose in the lower-dose arm (P for trend=.01) without reaching a level beyond which no further survival benefit was achieved. The survival benefit did not differ significantly between randomized comparisons where the higher-dose arm was hyperfractionated and those where it was not. There was heterogeneity in the median survival ratio by geographic region (P<.001), average age at randomization (P<.001), and year trial started (P for trend=.004), but not for proportion of patients with squamous cell carcinoma (P=.2). CONCLUSIONS: In trials with concurrent chemotherapy, higher radiation therapy doses resulted in poorer survival, possibly caused, at least in part, by high levels of toxicity. Where radiation therapy was given without chemotherapy, progressively higher radiation therapy doses resulted in progressively longer survival, and no upper dose level was found above which there was no further benefit. These findings support the consideration of further radiation therapy dose escalation trials, making use of modern treatment methods to reduce toxicity.

Conformity analysis to demonstrate reproducibility of target volumes for Margin-Intense Stereotactic Radiotherapy for borderline-resectable pancreatic cancer.

BACKGROUND AND PURPOSE: Margin-directed neoadjuvant radiotherapy for borderline-resectable pancreatic cancer (BRPC) aims to facilitate clear surgical margins. A systematic method was developed for definition of a boost target volume prior to a formal phase-I study. MATERIAL AND METHODS: Reference structures were defined by two oncologists and one radiologist, target structures were submitted by eight oncologist investigators and compared using conformity indices. Resultant risk of duodenal bleed (NTCP) was modelled. RESULTS: For GTV, reference volume was 2.1cm3 and investigator mean was 6.03cm3 (95% CI 3.92-8.13cm3), for boost volume 1.1cm3 and 1.25cm3 (1.02-1.48cm3). Mean Dice conformity coefficient for GTV was 0.47 (0.38-0.56), and for boost volume was significantly higher at 0.61 (0.52-0.70, p=0.01). Discordance index (DI) for GTV was 0.65 (0.56-0.75) and for boost volume was significantly lower at 0.39 (0.28-0.49, p=0.001). NTCP using reference contours was 2.95%, with mean for investigator contour plans 3.93% (3.63-4.22%). Correlations were seen between NTCP and GTV volume (p=0.02) and NTCP and DI (correlation coefficient 0.83 (0.29-0.97), p=0.01). CONCLUSIONS: Better conformity with reference was shown for boost volume compared with GTV. Investigator GTV volumes were larger than reference, had higher DI scores and modelled toxicity risk. A consistent method of target structure definition for margin-directed pancreatic radiotherapy is demonstrated.

An Analysis of Plan Robustness for Esophageal Tumors: Comparing Volumetric Modulated Arc Therapy Plans and Spot Scanning Proton Planning.

PURPOSE: Planning studies to compare x-ray and proton techniques and to select the most suitable technique for each patient have been hampered by the nonequivalence of several aspects of treatment planning and delivery. A fair comparison should compare similarly advanced delivery techniques from current clinical practice and also assess the robustness of each technique. The present study therefore compared volumetric modulated arc therapy (VMAT) and single-field optimization (SFO) spot scanning proton therapy plans created using a simultaneous integrated boost (SIB) for dose escalation in midesophageal cancer and analyzed the effect of setup and range uncertainties on these plans. METHODS AND MATERIALS: For 21 patients, SIB plans with a physical dose prescription of 2 Gy or 2.5 Gy/fraction in 25 fractions to planning target volume (PTV)50Gy or PTV62.5Gy (primary tumor with 0.5 cm margins) were created and evaluated for robustness to random setup errors and proton range errors. Dose-volume metrics were compared for the optimal and uncertainty plans, with P<.05 (Wilcoxon) considered significant. RESULTS: SFO reduced the mean lung dose by 51.4% (range 35.1%-76.1%) and the mean heart dose by 40.9% (range 15.0%-57.4%) compared with VMAT. Proton plan robustness to a 3.5% range error was acceptable. For all patients, the clinical target volume D98 was 95.0% to 100.4% of the prescribed dose and gross tumor volume (GTV) D98 was 98.8% to 101%. Setup error robustness was patient anatomy dependent, and the potential minimum dose per fraction was always lower with SFO than with VMAT. The clinical target volume D98 was lower by 0.6% to 7.8% of the prescribed dose, and the GTV D98 was lower by 0.3% to 2.2% of the prescribed GTV dose. CONCLUSIONS: The SFO plans achieved significant sparing of normal tissue compared with the VMAT plans for midesophageal cancer. The target dose coverage in the SIB proton plans was less robust to random setup errors and might be unacceptable for certain patients. Robust optimization to ensure adequate target coverage of SIB proton plans might be beneficial.

The Role of Oxygen in Avascular Tumor Growth.

The oxygen status of a tumor has significant clinical implications for treatment prognosis, with well-oxygenated subvolumes responding markedly better to radiotherapy than poorly supplied regions. Oxygen is essential for tumor growth, yet estimation of local oxygen distribution can be difficult to ascertain in situ, due to chaotic patterns of vasculature. It is possible to avoid this confounding influence by using avascular tumor models, such as tumor spheroids, a much better approximation of realistic tumor dynamics than monolayers, where oxygen supply can be described by diffusion alone. Similar to in situ tumours, spheroids exhibit an approximately sigmoidal growth curve, often approximated and fitted by logistic and Gompertzian sigmoid functions. These describe the basic rate of growth well, but do not offer an explicitly mechanistic explanation. This work examines the oxygen dynamics of spheroids and demonstrates that this growth can be derived mechanistically with cellular doubling time and oxygen consumption rate (OCR) being key parameters. The model is fitted to growth curves for a range of cell lines and derived values of OCR are validated using clinical measurement. Finally, we illustrate how changes in OCR due to gemcitabine treatment can be directly inferred using this model.

ARCII: A phase II trial of the HIV protease inhibitor Nelfinavir in combination with chemoradiation for locally advanced inoperable pancreatic cancer.

BACKGROUND AND PURPOSE: Nelfinavir can enhance intrinsic radiosensitivity, reduce hypoxia and improve vascularity. We conducted a phase II trial combining nelfinavir with chemoradiotherapy (CRT) for locally advanced inoperable pancreatic cancer (LAPC). MATERIALS AND METHODS: Radiotherapy (50.4Gy/28 fractions; boost to 59.4Gy/33 fractions) was administered with weekly gemcitabine and cisplatin. Nelfinavir started 3-10days before and was continued during CRT. The primary end-point was 1-year overall survival (OS). Secondary end-points included histological downstaging, radiological response, 1-year progression free survival (PFS), overall survival (OS) and treatment toxicity. An imaging sub-study (n=6) evaluated hypoxia ((18)F-Fluoromisonidazole-PET) and perfusion (perfusion CT) during induction nelfinavir. RESULTS: The study closed after recruiting 23 patients, due to non-availability of Nelfinavir in Europe. The 1-year OS was 73.4% (90% CI: 54.5-85.5%) and median OS was 17.4months (90% CI: 12.8-18.8). The 1-year PFS was 21.8% (90% CI: 8.9-38.3%) and median PFS was 5.5months (90% CI: 4.1-8.3). All patients experienced Grade 3/4 toxicity, but many were asymptomatic laboratory abnormalities. Four of 6 patients on the imaging sub-study demonstrated reduced hypoxia and increased perfusion post-nelfinavir. CONCLUSIONS: CRT combined with nelfinavir showed acceptable toxicity and promising survival in pancreatic cancer.

Estimating oxygen distribution from vasculature in three-dimensional tumour tissue.

Regions of tissue which are well oxygenated respond better to radiotherapy than hypoxic regions by up to a factor of three. If these volumes could be accurately estimated, then it might be possible to selectively boost dose to radio-resistant regions, a concept known as dose-painting. While imaging modalities such as 18F-fluoromisonidazole positron emission tomography (PET) allow identification of hypoxic regions, they are intrinsically limited by the physics of such systems to the millimetre domain, whereas tumour oxygenation is known to vary over a micrometre scale. Mathematical modelling of microscopic tumour oxygen distribution therefore has the potential to complement and enhance macroscopic information derived from PET. In this work, we develop a general method of estimating oxygen distribution in three dimensions from a source vessel map. The method is applied analytically to line sources and quasi-linear idealized line source maps, and also applied to full three-dimensional vessel distributions through a kernel method and compared with oxygen distribution in tumour sections. The model outlined is flexible and stable, and can readily be applied to estimating likely microscopic oxygen distribution from any source geometry. We also investigate the problem of reconstructing three-dimensional oxygen maps from histological and confocal two-dimensional sections, concluding that two-dimensional histological sections are generally inadequate representations of the three-dimensional oxygen distribution.

The effect of dose escalation on gastric toxicity when treating lower oesophageal tumours: a radiobiological investigation.

PURPOSE: Using radiobiological modelling to estimate normal tissue toxicity, this study investigates the effects of dose escalation for concurrent chemoradiation therapy (CRT) in lower third oesophageal tumours on the stomach. METHODS AND MATERIALS: 10 patients with lower third oesophageal cancer were selected from the SCOPE 1 database (ISCRT47718479) with a mean planning target volume (PTV) of 348 cm(3). The original 3D conformal plans (50 Gy3D) were compared to newly created RapidArc plans of 50 GyRA and 60 GyRA, the latter using a simultaneous integrated boost (SIB) technique using a boost volume, PTV2. Dose-volume metrics and estimates of normal tissue complication probability (NTCP) were compared. RESULTS: There was a significant increase in NTCP of the stomach wall when moving from the 50 GyRA to the 60 GyRA plans (11-17 %, Wilcoxon signed rank test, p = 0.01). There was a strong correlation between the NTCP values of the stomach wall and the volume of the stomach wall/PTV 1 and stomach wall/PTV2 overlap structures (R = 0.80 and R = 0.82 respectively) for the 60 GyRA plans. CONCLUSION: Radiobiological modelling suggests that increasing the prescribed dose to 60 Gy may be associated with a significantly increased risk of toxicity to the stomach. It is recommended that stomach toxicity be closely monitored when treating patients with lower third oesophageal tumours with 60 Gy.

Characterizing Heterogeneity within Head and Neck Lesions Using Cluster Analysis of Multi-Parametric MRI Data.

PURPOSE: To describe a methodology, based on cluster analysis, to partition multi-parametric functional imaging data into groups (or clusters) of similar functional characteristics, with the aim of characterizing functional heterogeneity within head and neck tumour volumes. To evaluate the performance of the proposed approach on a set of longitudinal MRI data, analysing the evolution of the obtained sub-sets with treatment. MATERIAL AND METHODS: The cluster analysis workflow was applied to a combination of dynamic contrast-enhanced and diffusion-weighted imaging MRI data from a cohort of squamous cell carcinoma of the head and neck patients. Cumulative distributions of voxels, containing pre and post-treatment data and including both primary tumours and lymph nodes, were partitioned into k clusters (k = 2, 3 or 4). Principal component analysis and cluster validation were employed to investigate data composition and to independently determine the optimal number of clusters. The evolution of the resulting sub-regions with induction chemotherapy treatment was assessed relative to the number of clusters. RESULTS: The clustering algorithm was able to separate clusters which significantly reduced in voxel number following induction chemotherapy from clusters with a non-significant reduction. Partitioning with the optimal number of clusters (k = 4), determined with cluster validation, produced the best separation between reducing and non-reducing clusters. CONCLUSION: The proposed methodology was able to identify tumour sub-regions with distinct functional properties, independently separating clusters which were affected differently by treatment. This work demonstrates that unsupervised cluster analysis, with no prior knowledge of the data, can be employed to provide a multi-parametric characterization of functional heterogeneity within tumour volumes.

A tumor control probability model for anal squamous cell carcinoma.

BACKGROUND AND PURPOSE: A recent update of the RTOG 9811, reported differing relapse rates for early and late anal squamous cell carcinoma following chemoradiotherapy (CRT). There may be a role for dose-individualization, however the dose-response relationship for anal cancer is not currently known. Intensity-modulated radiotherapy (IMRT) has been widely adopted with multiple series published. The aim is to fit a tumor control probability (TCP) model to the published IMRT data. MATERIALS AND METHODS: We performed a systematic review of PubMed and Embase databases to identify thirteen appropriate papers, including 625 patients. Predefined data fields were collected. A standard linear quadratic TCP model, which included repopulation, was fit by least squares minimization. RESULTS: The fitted TCP curve demonstrated a dose-response relationship with α=0.196 Gy(-1). The curve suggests: in early stage tumours, a dose reduction from 50 Gy to 45 Gy reduces 2 year local control from 98% to 95%; in late stage tumours, a dose escalation from 50 Gy to 55 Gy improves the 2 year local control rate from approximately 50% to 80%. CONCLUSIONS: The published data are broadly consistent with a linear quadratic dose-response model. Dose-individualization in anal cancer should be further investigated in the context of clinical trials.

Improving In Vivo High-Resolution CT Imaging of the Tumour Vasculature in Xenograft Mouse Models through Reduction of Motion and Bone-Streak Artefacts.

INTRODUCTION: Preclinical in vivo CT is commonly used to visualise vessels at a macroscopic scale. However, it is prone to many artefacts which can degrade the quality of CT images significantly. Although some artefacts can be partially corrected for during image processing, they are best avoided during acquisition. Here, a novel imaging cradle and tumour holder was designed to maximise CT resolution. This approach was used to improve preclinical in vivo imaging of the tumour vasculature. PROCEDURES: A custom built cradle containing a tumour holder was developed and fix-mounted to the CT system gantry to avoid artefacts arising from scanner vibrations and out-of-field sample positioning. The tumour holder separated the tumour from bones along the axis of rotation of the CT scanner to avoid bone-streaking. It also kept the tumour stationary and insensitive to respiratory motion. System performance was evaluated in terms of tumour immobilisation and reduction of motion and bone artefacts. Pre- and post-contrast CT followed by sequential DCE-MRI of the tumour vasculature in xenograft transplanted mice was performed to confirm vessel patency and demonstrate the multimodal capacity of the new cradle. Vessel characteristics such as diameter, and branching were quantified. RESULTS: Image artefacts originating from bones and out-of-field sample positioning were avoided whilst those resulting from motions were reduced significantly, thereby maximising the resolution that can be achieved with CT imaging in vivo. Tumour vessels ≥ 77 µm could be resolved and blood flow to the tumour remained functional. The diameter of each tumour vessel was determined and plotted as histograms and vessel branching maps were created. Multimodal imaging using this cradle assembly was preserved and demonstrated. CONCLUSIONS: The presented imaging workflow minimised image artefacts arising from scanner induced vibrations, respiratory motion and radiopaque structures and enabled in vivo CT imaging and quantitative analysis of the tumour vasculature at higher resolution than was possible before. Moreover, it can be applied in a multimodal setting, therefore combining anatomical and dynamic information.

Improved calibration of mass stopping power in low density tissue for a proton pencil beam algorithm.

Dose distributions for proton therapy treatments are almost exclusively calculated using pencil beam algorithms. An essential input to these algorithms is the patient model, derived from x-ray computed tomography (CT), which is used to estimate proton stopping power along the pencil beam paths. This study highlights a potential inaccuracy in the mapping between mass density and proton stopping power used by a clinical pencil beam algorithm in materials less dense than water. It proposes an alternative physically-motivated function (the mass average, or MA, formula) for use in this region. Comparisons are made between dose-depth curves calculated by the pencil beam method and those calculated by the Monte Carlo particle transport code MCNPX in a one-dimensional lung model. Proton range differences of up to 3% are observed between the methods, reduced to <1% when using the MA function. The impact of these range errors on clinical dose distributions is demonstrated using treatment plans for a non-small cell lung cancer patient. The change in stopping power calculation methodology results in relatively minor differences in dose when plans use three fields, but differences are observed at the 2%-2 mm level when a single field uniform dose technique is adopted. It is therefore suggested that the MA formula is adopted by users of the pencil beam algorithm for optimal dose calculation in lung, and that a similar approach is considered when beams traverse other low density regions such as the paranasal sinuses and mastoid process.