Search

Using co-culture experiments and mathematical modelling, Paczkowski et al discover that prostate cancer spheroids comprising mixed tumour cell populations display enhanced growth and reduced radiation sensitivity due to competitive and antagonistic interactions between cell populations. This interdisciplinary approach reveals a role for ecological-type interactions in the radiation response and may be used to study other cancer types.

Abstract 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.

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.

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.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.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.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.

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.

Oxygen levels in cancerous tissue can have a significant effect on treatment response: hypoxic tissue is both more radioresistant and more chemoresistant than well-oxygenated tissue. While recent advances in medical imaging have facilitated real-time observation of macroscopic oxygenation, the underlying physics limits the resolution to the millimetre domain, whereas oxygen tension varies over a micrometre scale. If the distribution of oxygen in the tumour micro-environment can be accurately estimated, then the effect of potential dose escalation to these hypoxic regions could be better modelled, allowing more realistic simulation of biologically adaptive treatments. Reaction–diffusion models are commonly used for modelling oxygen dynamics, with a variety of functional forms assumed for the dependence of oxygen consumption rate (OCR) on cellular status and local oxygen availability. In this work, we examine reaction–diffusion models of oxygen consumption in spherically and cylindrically symmetric geometries. We consider two different descriptions of oxygen consumption: one in which the rate of consumption is constant and one in which it varies with oxygen tension in a hyperbolic manner. In each case, we derive analytic approximations to the steady-state oxygen distribution, which are shown to closely match the numerical solutions of the equations and accurately predict the extent to which oxygen can diffuse. The derived expressions relate the limit to which oxygen can diffuse into a tissue to the OCR of that tissue. We also demonstrate that differences between these functional forms are likely to be negligible within the range of literature estimates of the hyperbolic oxygen constant, suggesting that the constant consumption rate approximation suffices for modelling oxygen dynamics for most values of OCR. These approximations also allow the rapid identification of situations where hyperbolic consumption forms can result in significant differences from constant consumption rate models, and so can reduce the computational workload associated with numerical solutions, by estimating both the oxygen diffusion distances and resultant oxygen profile. Such analysis may be useful for parameter fitting in large imaging datasets and histological sections, and allows easy quantification of projected differences between functional forms of OCR.

Background and Purpose:Complex-shaped malignant pleural mesotheliomas (MPMs) with challenging volumes are extremely difficult to treat by conventional radiotherapy due to tolerance doses of the surrounding normal tissue. In a feasibility study, we evaluated if inversely planned stereotactic intensity-modulated radiation therapy (IMRT) could be applied in the treatment of MPM. Patients and Methods: Results: Conclusion: Patients and Methods: Eight patients with unresectable lesions were treated after failure of chemotherapy. All patients were positioned using noninvasive patient fixation techniques which can be attached to the applied extracranial stereotactic system. Due to craniocaudal extension of the tumor, it was necessary to develop a special software attached to the inverse planning program KonRad, which can connect two inverse treatment plans and consider the applied dose of the first treatment plan in the area of the matchline of the second treatment plan. Results: Except for one patient, in whom radiotherapy was canceled due to abdominal metastasis, treatment could be completed in all patients and was well tolerated. Median survival after diagnosis was 20 months and after IMRT 6.5 months. Therefore, both the 1-year actuarial overall survival from the start of radiotherapy and the 2-year actuarial overall survival since diagnosis were 28%. IMRT did not result in clinically significant acute side effects. By using the described inverse planning software, overor underdosage in the region of the field matchline could be prevented. Pure treatment time ranged between 10 and 21 min. Conclusion:This study showed that IMRT is feasible in advanced unresectable MPM. The presented possibilities of stereotactic IMRT in the treatment of MPM will justify the evaluation of IMRT in early-stage pleural mesothelioma combined with chemotherapy in a study protocol, in order to improve the outcome of these patients. Furthermore, dose escalation should be possible by using IMRT. [ABSTRACT FROM AUTHOR]

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 () 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 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 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 : 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 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The enormous advances in the understanding of human anatomy, physiology and pathology in recent decades have led to ever-improving methods of disease prevention, diagnosis and treatment. Many of these achievements have been enabled, at least in part, by advances in ionizing radiation detectors. Radiology has been transformed by the implementation of multi-slice CT and digital x-ray imaging systems, with silver halide films now largely obsolete for many applications. Nuclear medicine has benefited from more sensitive, faster and higher-resolution detectors delivering ever-higher SPECT and PET image quality. PET/MR systems have been enabled by the development of gamma ray detectors that can operate in high magnetic fields. These huge advances in imaging have enabled equally impressive steps forward in radiotherapy delivery accuracy, with 4DCT, PET and MRI routinely used in treatment planning and online image guidance provided by cone-beam CT.The challenge of ensuring safe, accurate and precise delivery of highly complex radiation fields has also both driven and benefited from advances in radiation detectors. Detector systems have been developed for the measurement of electron, intensity-modulated and modulated arc x-ray, proton and ion beams, and around brachytherapy sources based on a very wide range of technologies. The types of measurement performed are equally wide, encompassing commissioning and quality assurance, reference dosimetry, in vivo dosimetry and personal and environmental monitoring.In this article, we briefly introduce the general physical characteristics and properties that are commonly used to describe the behaviour and performance of both discrete and imaging detectors. The physical principles of operation of calorimeters; ionization and charge detectors; semiconductor, luminescent, scintillating and chemical detectors; and radiochromic and radiographic films are then reviewed and their principle applications discussed. Finally, a general discussion of the application of detectors for x-ray nuclear medicine and ion beam imaging and dosimetry is presented. [ABSTRACT FROM AUTHOR]