Pulmonary magnetic resonance-guided online adaptive radiotherapy of locally advanced: the PUMA trial

BACKGROUND

Patients with locally-advanced non-small-cell lung cancer (LA-NSCLC) are often ineligible for surgery, so that definitive chemoradiotherapy (CRT) represents the treatment of choice. Nevertheless, long-term tumor control is often not achieved. Intensification of radiotherapy (RT) to improve locoregional tumor control is limited by the detrimental effect of higher radiation exposure of thoracic organs-at-risk (OAR). This narrow therapeutic ratio may be expanded by exploiting the advantages of magnetic resonance (MR) linear accelerators, mainly the online adaptation of the treatment plan to the current anatomy based on daily acquired MR images. However, MR-guidance is both labor-intensive and increases treatment times, which raises the question of its clinical feasibility to treat LA-NSCLC. Therefore, the PUMA trial was designed as a prospective, multicenter phase I trial to demonstrate the clinical feasibility of MR-guided online adaptive RT in LA-NSCLC.

METHODS

Thirty patients with LA-NSCLC in stage III A-C will be accrued at three German university hospitals to receive MR-guided online adaptive RT at two different MR-linac systems (MRIdian Linac®, View Ray Inc. and Elekta Unity®, Elekta AB) with concurrent chemotherapy. Conventionally fractioned RT with isotoxic dose escalation up to 70 Gy is applied. Online plan adaptation is performed once weekly or in case of major anatomical changes. Patients are followed-up by thoracic CT- and MR-imaging for 24 months after treatment. The primary endpoint is twofold: (1) successfully completed online adapted fractions, (2) on-table time. Main secondary endpoints include adaptation frequency, toxicity, local tumor control, progression-free and overall survival.

DISCUSSION

PUMA aims to demonstrate the clinical feasibility of MR-guided online adaptive RT of LA-NSCLC. If successful, PUMA will be followed by a clinical phase II trial that further investigates the clinical benefits of this approach. Moreover, PUMA is part of a large multidisciplinary project to develop MR-guidance techniques.

TRIAL REGISTRATION

ClinicalTrials.gov: NCT05237453 .

A deep image-to-image network organ segmentation algorithm for radiation treatment planning: principles and evaluation

Background

We describe and evaluate a deep network algorithm which automatically contours organs at risk in the thorax and pelvis on computed tomography (CT) images for radiation treatment planning.

Methods

The algorithm identifies the region of interest (ROI) automatically by detecting anatomical landmarks around the specific organs using a deep reinforcement learning technique. The segmentation is restricted to this ROI and performed by a deep image-to-image network (DI2IN) based on a convolutional encoder-decoder architecture combined with multi-level feature concatenation. The algorithm is commercially available in the medical products “syngo.via RT Image Suite VB50” and “AI-Rad Companion Organs RT VA20” (Siemens Healthineers). For evaluation, thoracic CT images of 237 patients and pelvic CT images of 102 patients were manually contoured following the Radiation Therapy Oncology Group (RTOG) guidelines and compared to the DI2IN results using metrics for volume, overlap and distance, e.g., Dice Similarity Coefficient (DSC) and Hausdorff Distance (HD₉₅). The contours were also compared visually slice by slice.

Results

We observed high correlations between automatic and manual contours. The best results were obtained for the lungs (DSC 0.97, HD₉₅ 2.7 mm/2.9 mm for left/right lung), followed by heart (DSC 0.92, HD₉₅ 4.4 mm), bladder (DSC 0.88, HD₉₅ 6.7 mm) and rectum (DSC 0.79, HD₉₅ 10.8 mm). Visual inspection showed excellent agreements with some exceptions for heart and rectum.

Conclusions

The DI2IN algorithm automatically generated contours for organs at risk close to those by a human expert, making the contouring step in radiation treatment planning simpler and faster. Few cases still required manual corrections, mainly for heart and rectum.

Statistical breathing curve sampling to quantify interplay effects of moving lung tumors in a 4D Monte Carlo dose calculation framework

PURPOSE

The interplay between respiratory tumor motion and dose application by intensity modulated radiotherapy (IMRT) techniques can potentially lead to undesirable and non-intuitive deviations from the planned dose distribution. We developed a 4D Monte Carlo (MC) dose recalculation framework featuring statistical breathing curve sampling, to precisely simulate the dose distribution for moving target volumes aiming at a comprehensive assessment of interplay effects.

METHODS

We implemented a dose accumulation tool that enables dose recalculations of arbitrary breathing curves including the actual breathing curve of the patient. This MC dose recalculation framework is based on linac log-files, facilitating a high temporal resolution up to 0.1 s. By statistical analysis of 128 different breathing curves, interplay susceptibility of different treatment parameters was evaluated for an exemplary patient case. To facilitate prospective clinical application in the treatment planning stage, in which patient breathing curves or linac log-files are not available, we derived a log-file free version with breathing curves generated by a random walk approach. Interplay was quantified by standard deviations σ in D_5%, D_50% and D_95%.

RESULTS

Interplay induced dose deviations for single fractions were observed and evaluated for IMRT and volumetric arc therapy (σ_D95% up to 1.3 %) showing a decrease with higher fraction doses and an increase with higher MU rates. Interplay effects for conformal treatment techniques were negligible (σ<0.1%). The log-file free version and the random walk generated breathing curves yielded similar results (deviations in σ< 0.1 %) and can be used as substitutes for interplay assessment.

CONCLUSION

It is feasible to combine statistically sampled breathing curves with MC dose calculations. The universality of the presented framework allows comprehensive assessment of interplay effects in retrospective and prospective clinically relevant scenarios.

Dosimetric Impact of the Positional Imaging Frequency for Hypofractionated Prostate Radiotherapy - A Voxel-by-Voxel Analysis

Background

To investigate deviations between planned and applied treatment doses for hypofractionated prostate radiotherapy and to quantify dosimetric accuracy in dependence of the image guidance frequency.

Methods

Daily diagnostic in-room CTs were carried out in 10 patients in treatment position as image guidance for hypofractionated prostate radiotherapy. Fraction doses were mapped to the planning CTs and recalculated, and applied doses were accumulated voxel-wise using deformable registration. Non-daily imaging schedules were simulated by deriving position correction vectors from individual scans and used to rigidly register the following scans until the next repositioning before dose recalculation and accumulation. Planned and applied doses were compared regarding dose-volume indices and TCP and NTCP values in dependence of the imaging and repositioning frequency.

Results

Daily image-guided repositioning was associated with only negligible deviations of analyzed dose-volume parameters and conformity/homogeneity indices for the prostate, bladder and rectum. Average CTV T did not significantly deviate from the plan values, and rectum NTCPs were highly comparable, while bladder NTCPs were reduced. For non-daily image-guided repositioning, there were significant deviations in the high-dose range from the planned values. Similarly, CTV dose conformity and homogeneity were reduced. While TCPs and rectal NTCPs did not significantly deteriorate for non-daily repositioning, bladder NTCPs appeared falsely diminished in dependence of the imaging frequency.

Conclusion

Using voxel-by-voxel dose accumulation, we showed for the first time that daily image-guided repositioning resulted in only negligible dosimetric deviations for hypofractionated prostate radiotherapy. Regarding dosimetric aberrations for non-daily imaging, daily imaging is required to adequately deliver treatment.

Dose-guided patient positioning in proton radiotherapy using multicriteria-optimization

Proton radiotherapy (PT) requires accurate target alignment before each treatment fraction, ideally utilizing 3D in-room X-ray computed tomography (CT) imaging. Typically, the optimal patient position is determined based on anatomical landmarks or implanted markers. In the presence of non-rigid anatomical changes, however, the planning scenario cannot be exactly reproduced and positioning should rather aim at finding the optimal position in terms of the actually applied dose. In this work, dose-guided patient alignment, implemented as multicriterial optimization (MCO) problem, was investigated in the scope of intensity-modulated and double-scattered PT (IMPT and DSPT) for the first time. A method for automatically determining the optimal patient position with respect to pre-defined clinical goals was implemented. Linear dose interpolation was used to access a continuous space of potential patient shifts. Fourteen head and neck (H&N) and eight prostate cancer patients with up to five repeated CTs were included. Dose interpolation accuracy was evaluated and the potential dosimetric advantages of dose-guided over bony-anatomy-based patient alignment investigated by comparison of clinically relevant target and organ-at-risk (OAR) dose-volume histogram (DVH) parameters. Dose interpolation was found sufficiently accurate with average pass-rates of 90% and 99% for an exemplary H&N and prostate patient, respectively, using a 2% dose-difference criterion. Compared to bony-anatomy-based alignment, the main impact of automated MCO-based dose-guided positioning was a reduced dose to the serial OARs (spinal cord and brain stem) for the H&N cohort. For the prostate cohort, under-dosage of the target structures could be efficiently diminished. Limitations of dose-guided positioning were mainly found in reducing target over-dosage due to weight loss for H&N patients, which might require adaptation of the treatment plan. Since labor-intense online quality-assurance is not required for dose-guided patient positioning, it might, nevertheless, be considered an interesting alternative to full online re-planning for initially mitigating the effects of anatomical changes.

Investigating CT to CBCT image registration for head and neck proton therapy as a tool for daily dose recalculation

PURPOSE

Intensity modulated proton therapy (IMPT) of head and neck (H&N) cancer patients may be improved by plan adaptation. The decision to adapt the treatment plan based on a dose recalculation on the current anatomy requires a diagnostic quality computed tomography (CT) scan of the patient. As gantry-mounted cone beam CT (CBCT) scanners are currently being offered by vendors, they may offer daily or weekly updates of patient anatomy. CBCT image quality may not be sufficient for accurate proton dose calculation and it is likely necessary to perform CBCT CT number correction. In this work, the authors investigated deformable image registration (DIR) of the planning CT (pCT) to the CBCT to generate a virtual CT (vCT) to be used for proton dose recalculation.

METHODS

Datasets of six H&N cancer patients undergoing photon intensity modulated radiation therapy were used in this study to validate the vCT approach. Each dataset contained a CBCT acquired within 3 days of a replanning CT (rpCT), in addition to a pCT. The pCT and rpCT were delineated by a physician. A Morphons algorithm was employed in this work to perform DIR of the pCT to CBCT following a rigid registration of the two images. The contours from the pCT were deformed using the vector field resulting from DIR to yield a contoured vCT. The DIR accuracy was evaluated with a scale invariant feature transform (SIFT) algorithm comparing automatically identified matching features between vCT and CBCT. The rpCT was used as reference for evaluation of the vCT. The vCT and rpCT CT numbers were converted to stopping power ratio and the water equivalent thickness (WET) was calculated. IMPT dose distributions from treatment plans optimized on the pCT were recalculated with a Monte Carlo algorithm on the rpCT and vCT for comparison in terms of gamma index, dose volume histogram (DVH) statistics as well as proton range. The DIR generated contours on the vCT were compared to physician-drawn contours on the rpCT.

RESULTS

The DIR accuracy was better than 1.4 mm according to the SIFT evaluation. The mean WET differences between vCT (pCT) and rpCT were below 1 mm (2.6 mm). The amount of voxels passing 3%/3 mm gamma criteria were above 95% for the vCT vs rpCT. When using the rpCT contour set to derive DVH statistics from dose distributions calculated on the rpCT and vCT the differences, expressed in terms of 30 fractions of 2 Gy, were within [-4, 2 Gy] for parotid glands (D(mean)), spinal cord (D(2%)), brainstem (D(2%)), and CTV (D(95%)). When using DIR generated contours for the vCT, those differences ranged within [-8, 11 Gy].

CONCLUSIONS

In this work, the authors generated CBCT based stopping power distributions using DIR of the pCT to a CBCT scan. DIR accuracy was below 1.4 mm as evaluated by the SIFT algorithm. Dose distributions calculated on the vCT agreed well to those calculated on the rpCT when using gamma index evaluation as well as DVH statistics based on the same contours. The use of DIR generated contours introduced variability in DVH statistics.

Feasibility of automated proton therapy plan adaptation for head and neck tumors using cone beam CT images

Background

Intensity modulated proton therapy (IMPT) of head and neck (H&N) tumors may benefit from plan adaptation to correct for the dose perturbations caused by weight loss and tumor volume changes observed in these patients. As cone beam CT (CBCT) is increasingly considered in proton therapy, it may be possible to use available CBCT images following intensity correction for plan adaptation. This is the first study exploring IMPT plan adaptation on CBCT images corrected and delineated by deformable image registration of the planning CT (pCT) to the CBCT, yielding a virtual CT (vCT).

Methods

A Morphons algorithm was used to deform the pCTs and corresponding delineations of 9 H&N cancer patients to a weekly CBCT acquired within ±3 days of a control replanning CT scan (rpCT). The IMPT treatment plans were adapted using the vCT and the adapted and original plans were recalculated on the rpCT for dose/volume parameter evaluation of the impact of adaptation.

Results

On the rpCT, the adapted plans were equivalent to the original plans in terms of target volumes D₉₅ and V₉₅, but showed a significant reduction of D₂ in these volumes. OAR doses were mostly equivalent or reduced. In particular, the adapted plans did not reduce parotid gland D_mean, but the dose to the optical system. For three patients the spinal cord or brain stem received higher, though well below tolerance, maximum dose. Subsequent tightening of the treatment planning constraints for these OARs on new vCT-adapted plans did not degrade target coverage and yielded pCT equivalent plans on the vCT.

Conclusions

An offline automated procedure to generate an adapted IMPT plan on CBCT images was developed and investigated. When evaluating the adapted plan on a control rpCT we observed reduced D₂ in target volumes as major improvement. OAR sparing was only partially improved by the procedure. Despite potential limitations in the accuracy of the vCT approach, an improved quality of the adapted plans could be achieved.

Electronic supplementary material

The online version of this article (doi:10.1186/s13014-016-0641-7) contains supplementary material, which is available to authorized users.

Long-term results in malignant pleural mesothelioma treated with neoadjuvant chemotherapy, extrapleural pneumonectomy and intensity-modulated radiotherapy

Introduction

We investigated the clinical outcome and the toxicity of trimodal therapy of malignant pleural mesothelioma (MPM) treated with neoadjuvant chemotherapy, extrapleural pneumonectomy (EPP) and adjuvant intensity-modulated radiotherapy (IMRT).

Methods

Chemotherapy regimens included Cisplatin/Pemetrexed, Carboplatin/Pemetrexed and Cisplatin/Gemcitabine, followed by EPP. 62 patients completed the adjuvant radiotherapy. IMRT was carried out in two techniques, either step&shoot or helical tomotherapy. Median target dose was 48 Gy to 54 Gy. Toxicity was scored with the Common Terminology Criteria (CTC) for Adverse Events. We used Kaplan-Meier method to estimate actuarial rate of locoregional control (LRC), distant control (DC) and overall survival (OS), measured from the date of surgery. Rates were compared using the logrank test. For multivariate analysis the Cox proportional hazard model was used.

Results

The median OS, LRC and DC times were 20.4, 31.4 and 21.4 months. The 1-, 2-, 3-year OS rates were 63, 42, 28 %, the LRC rates were 81, 60, 40 %, and the DC rates were 62, 48, 41 %. We observed no CTC grade 4 or grade 5 toxicity. Step&shoot and helical tomotherapy were equivalent both in dosimetric characteristics and clinical outcome. Biphasic tumor histology was associated with worse clinical outcome compared to epitheloid histology.

Conclusions

Mature clinical results of trimodal treatment for MPM were presented. They indicate that hemithoracic radiotherapy after EPP can be safely administered by either step&shoot IMRT and tomotherapy. However, the optimal prospective patient selection for this aggressive trimodal therapy approach remains unclear. This study can serve as a benchmark for current and future therapy concepts for MPM.

Electronic supplementary material

The online version of this article (doi:10.1186/s13014-015-0575-5) contains supplementary material, which is available to authorized users.

Comparing cone-beam CT intensity correction methods for dose recalculation in adaptive intensity-modulated photon and proton therapy for head and neck cancer

BACKGROUND

Adaptive intensity-modulated photon and proton radiotherapy (IMRT and IMPT) of head and neck (H&N) cancer requires frequent three-dimensional (3D) dose calculation. We compared two approaches for dose recalculation on the basis of intensity-corrected cone-beam (CB) x-ray computed tomography (CT) images.

MATERIAL AND METHODS

For nine H&N tumor patients, virtual CTs (vCT) were generated by deformable image registration of the planning CT (pCT) to the CBCT. The second intensity correction approach used population-based lookup tables for scaling CBCT intensities to the pCT HU range (CBCTLUT). IMRT and IMPT plans were generated with a commercial treatment planning system. Dose recalculations on vCT and CBCTLUT were analyzed using a (3%, 3 mm) gamma-index analysis and comparison of normal tissue and tumor dose/volume parameters. A replanning CT (rpCT) acquired within three days of the CBCT served as reference. Single field uniform dose (SFUD) proton plans were created and recalculated on vCT and CBCTLUT for proton range comparison.

RESULTS

Dose/volume parameters showed minor differences between rpCT, vCT and CBCTLUT in IMRT, but clinically relevant deviations between CBCTLUT and rpCT in the spinal cord for IMPT. Gamma-index pass-rates were found increased for vCT with respect to CBCTLUT in IMPT (by up to 21 percentage points) and IMRT (by up to 9 percentage points) for most cases. The SFUD-based proton range assessment showed improved agreement of vCT and rpCT, with 88-99% of the depth dose profiles in beam's eye view agreeing within 3 mm. For CBCTLUT, only 80-94% of the profiles fulfilled this criterion.

CONCLUSION

vCT and CBCTLUT are suitable options for dose recalculation in adaptive IMRT. In the scope of IMPT, the vCT approach is preferable.

In vivo measurement of dose distribution in patients' lymphocytes: helical tomotherapy versus step-and-shoot IMRT in prostate cancer

In radiotherapy, in vivo measurement of dose distribution within patients' lymphocytes can be performed by detecting gamma-H2AX foci in lymphocyte nuclei. This method can help in determining the whole-body dose. Options for risk estimations for toxicities in normal tissue and for the incidence of secondary malignancy are still under debate. In this investigation, helical tomotherapy (TOMO) is compared with step-and-shoot IMRT (SSIMRT) of the prostate gland by measuring the dose distribution within patients' lymphocytes. In this prospective study, blood was taken from 20 patients before and 10 min after their first irradiation fraction for each technique. The isolated leukocytes were fixed 2 h after radiation. DNA double-stranded breaks in lymphocyte nuclei were stained immunocytochemically using anti-gamma-H2AX antibodies. Gamma-H2AX foci distribution in lymphocytes was determined for each patient. Using a calibration line, dose distributions in patients' lymphocytes were determined by studying the gamma-H2AX foci distribution, and these data were used to generate a cumulative dose-lymphocyte histogram (DLH). Measured in vivo (DLH), significantly fewer lymphocytes indicated low-dose exposure (<40% of the applied dose) during TOMO compared with SSIMRT. The dose exposure range, between 45 and 100%, was equal with both radiation techniques. The mean number of gamma-H2AX foci per lymphocyte was significantly lower in the TOMO group compared with the SSIMRT group. In radiotherapy of the prostate gland, TOMO generates a smaller fraction of patients' lymphocytes with low-dose exposure relative to the whole body compared with SSIMRT. Differences in the constructional buildup of the different linear accelerator systems, e.g. the flattening filter, may be the cause thereof. The influence of these methods on the incidence of secondary malignancy should be investigated in further studies.

Combined use of (18)F-FDG and (18)F-FMISO in unresectable non-small cell lung cancer patients planned for radiotherapy: a dynamic PET/CT study

Aim of this study was to evaluate and compare, by means of dynamic and static PET/CT, the distribution patterns and pharmacokinetics of fluorine-18 fluorodeoxyglucose ((18)F-FDG) and of fluorine-18-fluoromisonidazole ((18)F-FMISO) in non-small cell lung cancer (NSCLC) patients scheduled for intensity modulated radiation therapy (IMRT). Thirteen patients suffering from inoperable stage III NSCLC underwent PET/CTs with (18)F-FDG and (18)F-FMISO for tumor metabolism and hypoxia assessment accordingly. Evaluation of PET/CT studies was based on visual analysis, semi-quantitative (SUV) calculations and absolute quantitative estimations, after application of a two-tissue compartment model and a non-compartmental approach. (18)F-FDG PET/CT revealed all thirteen primary lung tumors as sites of increased (18)F-FDG uptake. Six patients demonstrated also in total 43 (18)F-FDG avid metastases; these patients were excluded from radiotherapy. (18)F-MISO PET/CT demonstrated 12/13 primary lung tumors with faint tracer uptake. Only one tumor was clearly (18)F-FMISO avid, (SUVaverage = 3.4, SUVmax = 5.0). Mean values for (18)F-FDG, as derived from dPET/CT data, were SUVaverage = 8.9, SUVmax = 15.1, K1 = 0.23, k2 = 0.53, k3 = 0.17, k4 = 0.02, influx = 0.05 and fractal dimension (FD) = 1.25 for the primary tumors. The respective values for (18)F-FMISO were SUVaverage = 1.4, SUVmax = 2.2, K1 = 0.26, k2 = 0.56, k3 = 0.06, k4 = 0.06, influx = 0.02 and FD = 1.14. No statistically significant correlation was observed between the two tracers. (18)F-FDG PET/CT changed therapy management in six patients, by excluding them from planned IMRT. (18)F-FMISO PET/CT revealed absence of significant tracer uptake in the majority of the (18)F-FDG avid NSCLCs. Lack of correlation between the two tracers' kinetics indicates that they reflect different molecular mechanisms and implies the discordance between increased glycolysis and hypoxia in the malignancy.

Software for quantitative analysis of radiotherapy: overview, requirement analysis and design solutions

Radiotherapy is a fast-developing discipline which plays a major role in cancer care. Quantitative analysis of radiotherapy data can improve the success of the treatment and support the prediction of outcome. In this paper, we first identify functional, conceptional and general requirements on a software system for quantitative analysis of radiotherapy. Further we present an overview of existing radiotherapy analysis software tools and check them against the stated requirements. As none of them could meet all of the demands presented herein, we analyzed possible conceptional problems and present software design solutions and recommendations to meet the stated requirements (e.g. algorithmic decoupling via dose iterator pattern; analysis database design). As a proof of concept we developed a software library "RTToolbox" following the presented design principles. The RTToolbox is available as open source library and has already been tested in a larger-scale software system for different use cases. These examples demonstrate the benefit of the presented design principles.

Multimodal hypoxia imaging and intensity modulated radiation therapy for unresectable non-small-cell lung cancer: the HIL trial

Background

Radiotherapy, preferably combined with chemotherapy, is the treatment standard for locally advanced, unresectable non-small cell lung cancer (NSCLC). The tumor response to different therapy protocols is variable, with hypoxia known to be a major factor that negatively influences treatment effectiveness. Visualisation of tumor hypoxia prior to the use of modern radiation therapy strategies, such as intensity modulated radiation therapy (IMRT), might allow optimized dose applications to the target volume, leading to improvement of therapy outcome. ¹⁸ F-fluoromisonidazole dynamic positron emission tomography and computed tomography (¹⁸ F-FMISO dPET-CT) and functional magnetic resonance imaging (functional MRI) are attractive options for imaging tumor hypoxia.

Methods/design

The HIL trial is a single centre study combining multimodal hypoxia imaging with ¹⁸ F-FMISO dPET-CT and functional MRI, with intensity modulated radiation therapy (IMRT) in patients with inoperable stage III NSCLC. 15 patients will be recruited in the study. All patients undergo initial FDG PET-CT and serial ¹⁸ F-FMISO dPET-CT and functional MRI before treatment, at week 5 of radiotherapy and 6 weeks post treatment. Radiation therapy is performed as inversely planned IMRT based on 4D-CT.

Discussion

Primary objectives of the trial are to characterize the correlation of ¹⁸ F-FMISO dPET-CT and functional MRI for tumor hypoxia imaging in NSCLC and evaluate possible effects of radiation therapy on tumor re-oxygenation. Further objectives include the generation of data regarding the prognostic value of ¹⁸ F-FMISO dPET-CT and functional MRI for locoregional control, progression free survival and overall survival of NSCLC treated with IMRT, which will form the basis for larger clinical trials focusing on possible interactions between tumor oxygenation and radiotherapy outcome.

Trial registration

The ClinicalTrials.gov protocol ID is NCT01617980

A stochastic approach to estimate the uncertainty of dose mapping caused by uncertainties in b-spline registration

PURPOSE

In fractionated radiation therapy, image guidance with daily tomographic imaging becomes more and more clinical routine. In principle, this allows for daily computation of the delivered dose and for accumulation of these daily dose distributions to determine the actually delivered total dose to the patient. However, uncertainties in the mapping of the images can translate into errors of the accumulated total dose, depending on the dose gradient. In this work, an approach to estimate the uncertainty of mapping between medical images is proposed that identifies areas bearing a significant risk of inaccurate dose accumulation.

METHODS

This method accounts for the geometric uncertainty of image registration and the heterogeneity of the dose distribution, which is to be mapped. Its performance is demonstrated in context of dose mapping based on b-spline registration. It is based on evaluation of the sensitivity of dose mapping to variations of the b-spline coefficients combined with evaluation of the sensitivity of the registration metric with respect to the variations of the coefficients. It was evaluated based on patient data that was deformed based on a breathing model, where the ground truth of the deformation, and hence the actual true dose mapping error, is known.

RESULTS

The proposed approach has the potential to distinguish areas of the image where dose mapping is likely to be accurate from other areas of the same image, where a larger uncertainty must be expected.

CONCLUSIONS

An approach to identify areas where dose mapping is likely to be inaccurate was developed and implemented. This method was tested for dose mapping, but it may be applied in context of other mapping tasks as well.

Automated vs. manual pattern recognition of 3D (1)H MRSI data of patients with prostate cancer

RATIONALE AND OBJECTIVES

The aim of this study was to assess (1) automated analysis methods versus manual evaluation by human experts of three-dimensional proton magnetic resonance spectroscopic imaging (MRSI) data from patients with prostate cancer and (2) the contribution of spatial information to decision making.

MATERIALS AND METHODS

Three-dimensional proton MRSI was applied at 1.5 T. MRSI data from 10 patients with histologically proven prostate adenocarcinoma, scheduled either for prostatectomy or intensity-modulated radiation therapy, were evaluated. First, two readers manually labeled spectra using spatial information to identify the localization of spectra and neighborhood information, establishing the reference set of this study. Then, spectra were labeled again manually in a blinded and randomized manner and evaluated automatically using software that applied spectral line fitting as well as pattern recognition routines. Statistical analysis of the results of the different approaches was performed.

RESULTS

Altogether, 1018 spectra were evaluable by all methods. Numbers of evaluable spectra differed significantly depending on patient and evaluation method. Compared to automated analysis, the readers made rather binary decisions, using information from neighboring spectra in ambiguous cases, when evaluating MRSI data as a whole. Differences between anatomically blinded and unblinded evaluation were larger than differences between evaluations using blinded data and automated techniques.

CONCLUSIONS

An automated approach, which evaluates each spectrum individually, can be as good as an anatomy-blinded human reader. Spatial information is routinely used by human experts to support their final decisions. Automated procedures that consider anatomic information for spectral evaluation will enhance the diagnostic impact of MRSI of the human prostate.

Diffusion-weighted imaging-based probabilistic segmentation of high- and low-proliferative areas in high-grade gliomas

The apparent diffusion coefficient (ADC) derived from diffusion-weighted imaging (DWI) correlates inversely with tumor proliferation rates. High-grade gliomas are typically heterogeneous and the delineation of areas of high and low proliferation is impeded by partial volume effects and blurred borders. Commonly used manual delineation is further impeded by potential overlap with cerebrospinal fluid and necrosis. Here we present an algorithm to reproducibly delineate and probabilistically quantify the ADC in areas of high and low proliferation in heterogeneous gliomas, resulting in a reproducible quantification in regions of tissue inhomogeneity. We used an expectation maximization (EM) clustering algorithm, applied on a Gaussian mixture model, consisting of pure superpositions of Gaussian distributions. Soundness and reproducibility of this approach were evaluated in 10 patients with glioma. High- and low-proliferating areas found using the clustering correspond well with conservative regions of interest drawn using all available imaging data. Systematic placement of model initialization seeds shows good reproducibility of the method. Moreover, we illustrate an automatic initialization approach that completely removes user-induced variability. In conclusion, we present a rapid, reproducible and automatic method to separate and quantify heterogeneous regions in gliomas.

Software-guided standardization of manual landmark data in medical images

It is crucial to evaluate registration algorithms in order to make them available in clinical practice. Several evaluation strategies have been proposed in the past, and one approach is to evaluate these algorithms with intrinsic anatomical landmarks identified by a health professional. The acquisition and handling of large amounts of these landmark data is a time-consuming task for the health professional, and it is vulnerable to errors and inconsistencies. Additionally, limited access to appropriate tools makes dealing with landmark data considerably more difficult. We introduce a strategy for the acquisition of landmarks for the landmark-based evaluation of registration algorithms and we present an ontology-driven software tool that assists the different partners involved to act according to that strategy. This tool provides the user with intrinsic knowledge of the registration problems, the possibility to conveniently make the acquired data available to further processing, and an easy-to-use graphical interface.

Respiratory-induced prostate motion: characterization and quantification in dynamic MRI

BACKGROUND AND PURPOSE

To investigate prostate movement during deep breathing and contraction of abdominal musculature by means of dynamic MRI and analyze implications for image-guided radiotherapy of prostate cancer.

PATIENTS AND METHODS

A total of 43 patients and 8 healthy volunteers were examined with MRI. Images during deep respiration and during contraction of abdominal musculature (via a coughing maneuver) were obtained with dynamic two-dimensional (2D) balanced SSFP; 3 frames/s were obtained over an acquisition time of 15 s. Images were acquired in sagittal orientation to evaluate motion along both the craniocaudal (cc)-axis and anteroposterior (ap)-axis. Prostate motion was quantified semi-automatically using dedicated software tools.

RESULTS

Respiratory induced mean cc-axis displacement of the prostate was 2.7 ± 1.9 (SD) mm (range, 0.5-10.6 mm) and mean ap-axis displacement 1.8 ± 1.0 (SD) mm (range, 0.3-10 mm). In 69% of the subjects, breathing-related prostate movements were found to be negligible (< 3 mm). The prostate displacement for abdominal contraction was significantly higher: mean cc-axis displacement was max. 8.4 ± 6.7 (SD) mm (range, 0.6-27 mm); mean anteroposterior movement was 8.3 ± 7.7 (SD) mm (range, 0.7-26 mm).

CONCLUSION

Dynamic MRI is an excellent tool for noninvasive real-time imaging of prostate movement. Further investigations regarding possible applications in image-guided radiotherapy, e.g., for individualized planning and in integrated linac/MRI systems, are warranted.

Comparative analysis of Pareto surfaces in multi-criteria IMRT planning

In the multi-criteria optimization approach to IMRT planning, a given dose distribution is evaluated by a number of convex objective functions that measure tumor coverage and sparing of the different organs at risk. Within this context optimizing the intensity profiles for any fixed set of beams yields a convex Pareto set in the objective space. However, if the number of beam directions and irradiation angles are included as free parameters in the formulation of the optimization problem, the resulting Pareto set becomes more intricate. In this work, a method is presented that allows for the comparison of two convex Pareto sets emerging from two distinct beam configuration choices. For the two competing beam settings, the non-dominated and the dominated points of the corresponding Pareto sets are identified and the distance between the two sets in the objective space is calculated and subsequently plotted. The obtained information enables the planner to decide if, for a given compromise, the current beam setup is optimal. He may then re-adjust his choice accordingly during navigation. The method is applied to an artificial case and two clinical head neck cases. In all cases no configuration is dominating its competitor over the whole Pareto set. For example, in one of the head neck cases a seven-beam configuration turns out to be superior to a nine-beam configuration if the highest priority is the sparing of the spinal cord. The presented method of comparing Pareto sets is not restricted to comparing different beam angle configurations, but will allow for more comprehensive comparisons of competing treatment techniques (e.g., photons versus protons) than with the classical method of comparing single treatment plans.

Biological in-vivo measurement of dose distribution in patients' lymphocytes by gamma-H2AX immunofluorescence staining: 3D conformal- vs. step-and-shoot IMRT of the prostate gland

BACKGROUND

Different radiation-techniques in treating local staged prostate cancer differ in their dose- distribution. Physical phantom measurements indicate that for 3D, less healthy tissue is exposed to a relatively higher dose compared to SSIMRT. The purpose is to substantiate a dose distribution in lymphocytes in-vivo and to discuss the possibility of comparing it to the physical model of total body dose distribution.

METHODS

For each technique (3D and SSIMRT), blood was taken from 20 patients before and 10 min after their first fraction of radiotherapy. The isolated leukocytes were fixed 2 hours after radiation. DNA double-strand breaks (DSB) in lymphocytes' nuclei were stained immunocytochemically using the gamma-H2AX protein. Gamma-H2AX foci inside each nucleus were counted in 300 irradiated as well as 50 non-irradiated lymphocytes per patient. In addition, lymphocytes of 5 volunteer subjects were irradiated externally at different doses and processed under same conditions as the patients' lymphocytes in order to generate a calibration-line. This calibration-line assigns dose-value to mean number of gamma-H2AX foci/ nucleus. So the dose distributions in patients' lymphocytes were determined regarding to the gamma-H2AX foci distribution. With this information a cumulative dose-lymphocyte-histogram (DLH) was generated. Visualized distribution of gamma-H2AX foci, correspondingly dose per nucleus, was compared to the technical dose-volume-histogram (DVH), related to the whole body-volume.

RESULTS

Measured in-vivo (DLH) and according to the physical treatment-planning (DVH), more lymphocytes resulted with low-dose exposure (< 20% of the applied dose) and significantly fewer lymphocytes with middle-dose exposure (30%-60%) during Step-and-Shoot-IMRT, compared to conventional 3D conformal radiotherapy. The high-dose exposure (> 80%) was equal in both radiation techniques. The mean number of gamma-H2AX foci per lymphocyte was 0.49 (3D) and 0.47 (SSIMRT) without significant difference.

CONCLUSIONS

In-vivo measurement of the dose distribution within patients' lymphocytes can be performed by detecting gamma-H2AX foci. In case of 3D and SSIMRT, the results of this method correlate with the physical calculated total body dose-distribution, but cannot be interpreted unrestrictedly due to the blood circulation. One possible application of the present method could be in radiation-protection for in-vivo dose estimation after accidental exposure to radiation.

Intensity modulated or fractionated stereotactic reirradiation in patients with recurrent nasopharyngeal cancer

PURPOSE

To report our experience with intensity-modulated or stereotactic reirradiation in patients suffering from recurrent nasopharyngeal carcinoma

PATIENTS AND METHODS

The records of 17 patients with recurrent nasopharygeal carcinoma treated by intensity-modulated (n = 14) or stereotactic (n = 3) reirradiation in our institution were reviewed. Median age was 53 years and most patients (n = 14) were male. The majority of tumors showed undifferentiated histology (n = 14) and infiltration of intracranial structures (n = 12). Simultaneous systemic therapy was applied in 8 patients. Initial treatment covered the gross tumor volume with a median dose of 66 Gy (50-72 Gy) and the cervical nodal regions with a median dose of 56 Gy (50-60 Gy). Reirradiation was confined to the local relapse region with a median dose of 50.4 Gy (36-64Gy), resulting in a median cumulative dose of 112 Gy (91-134 Gy). The median time interval between initial and subsequent treatment was 52 months (6-132).

RESULTS

The median follow up for the entire cohort was 20 months and 31 months for survivors (10-84). Five patients (29%) developed isolated local recurrences and three patients (18%) suffered from isolated nodal recurrences. The actuarial 1- and 2-year rates of local/locoregional control were 76%/59% and 69%/52%, respectively. Six patients developed distant metastasis during the follow up period. The median actuarial overall survival for the entire cohort was 23 months, transferring into 1-, 2-, and 3-year overall survival rates of 82%, 44% and 37%. Univariate subset analyses showed significantly increased overall survival and local control for patients with less advanced rT stage, retreatment doses > 50 Gy, concurrent systemic treatment and complete response. Severe late toxicity (Grad III) attributable to reirradiation occurred in five patients (29%), particularly as hearing loss, alterations of taste/smell, cranial neuropathy, trismus and xerostomia.

CONCLUSION

Reirradiation with intensity-modulated or stereotactic techniques in recurrent nasopharyngeal carcinoma is feasible and yields encouraging results in terms of local control and overall survival in patients with acceptable toxicity in patients with less advanced recurrences. However, the achievable outcome is limited in patients with involvement of intracranial structures, emphasising the need for close monitoring after primary therapy.

IMRT reirradiation with concurrent cetuximab immunotherapy in recurrent head and neck cancer

PURPOSE

In this retrospective investigation, the outcome and toxicity after reirradiation with concurrent cetuximab immunotherapy of recurrent head and neck cancer (HNC) in patients who had contraindications to platinum-based chemotherapy were analyzed.

MATERIALS AND METHODS

Ten patients with locally advanced recurrent HNC were retrospectively evaluated. In 9 cases, histology was squamous cell carcinoma, in one case adenoid cystic carcinoma. External beam radiotherapy was part of the initial treatment in all cases. Reirradiation was carried out using step-and-shoot intensity-modulated radiotherapy (IMRT) with a median dose of 50.4 Gy. Cetuximab was applied as loading dose (400 mg/m(2)) 1 week prior to reirradiation and then weekly concurrently with radiotherapy (250 mg/m(2)).

RESULTS

The median overall survival time after initiation of reirradiation was 7 months; the 1-year overall survival (OS) rate was 40%. Local failure was found in 3 patients, resulting in a 1-year local control (LC) rate of 61%. The 1-year locoregional control (LRC) rate was 44%, while the 1-year distant metastasis-free survival (DMFS) was 75%. Acute hematological toxicity was not observed in the group. Severe acute toxicity included one fatal infield arterial bleeding and one flap necrosis. Severe late toxicities were noted in 2 patients: fibrosis of the temporomandibular joint in 1 patient and stenosis of the cervical esophagus in another.

CONCLUSIONS

IMRT reirradiation with concurrent cetuximab immunotherapy in recurrent HNC is feasible with acceptable acute toxicity. Further investigations are necessary to determine the clinical role of this therapy concept.

Comparison of peripheral zone and central gland volume in patients undergoing intensity-modulated radiotherapy

BACKGROUND AND PURPOSE

The shrinking effect of androgen deprivation therapy (ADT) and radiotherapy (RT) on prostate gland volume is a known clinical finding. Until now, it is not known which part of the prostate shrinks more. We examined patients with and without ADT undergoing intensity-modulated RT (IMRT) and performed 3-dimensional measurements of the peripheral zone (PZ) and central gland (CG) with magnetic resonance imaging (MRI).

METHODS AND MATERIALS

Prostate gland volumes of PZ and CG between planning MRI and first available follow-up MRI were retrospectively determined in 44 patients with localized prostate carcinoma. A total of 24 patients had ADT with a median time interval of 5 months (range, 1.5-24 months). Median time interval between both MRI time points was 132 days (range, 104-224 days). Two observers performed PZ and CG delineation in consensus using planimetry. Volume changes over time were determined and compared.

RESULTS

Patients who had ADT showed smaller prostate volume in the first MRI (mean [SD], 32 [16.7] mL), which was still present after IMRT (28.1 [16.7] mL). Patients who had no ADT started with 44.6 (16.9) mL and showed 37.5 (13.9) mL after IMRT. Shrinking effect in PZ was significantly larger than in CG for all patients (-18.3% vs -6.3%, P < 0.05).

CONCLUSIONS

Because, typically, most tumors are located in PZ and this area also shows the largest shrinkage effect after IMRT, this should be taken into account for planning purposes. Notably, there are only minor differences in the relative shrinking effects between patients with and without ADT, although they start with different volumes.

New multileaf collimator with a leaf width of 5 mm improves plan quality compared to 10 mm in step-and-shoot IMRT of HNC using integrated boost procedure

PURPOSE

To investigate whether a new multileaf collimator with a leaf width of 5 mm (MLC-5) over the entire field size of 40 x 40 cm(2) improves plan quality compared to a leaf width of 10 mm (MLC-10) in intensity-modulated radiotherapy (IMRT) with integrated boost for head and neck cancer.

PATIENTS AND METHODS

A plan comparison was performed for ten patients with head and neck cancer. For each patient, seven plans were calculated: one plan with MLC-10 and nine beams, four plans with MLC-5 and nine beams (with different intensity levels and two-dimensional median filter sizes [2D-MFS]), and one seven-beam plan with MLC-5 and MLC-10, respectively. Isocenter, beam angles and planning constraints were not changed. Mean values of common plan parameters over all ten patients were estimated, and plan groups of MLC-5 and MLC-10 with nine and seven beams were compared.

RESULTS

The use of MLC-5 led to a significantly higher conformity index and an improvement of the 90% coverage of PTV1 (planning target volume) and PTV2 compared with MLC-10. This was noted in the nine- and seven-beam plans. Within the nine-beam group with MLC-5, a reduction of the segment number by up to 25% at reduced intensity levels and for increased 2D-MFS did not markedly worsen plan quality. Interestingly, a seven-beam IMRT with MLC-5 was inferior to a nine-beam IMRT with MLC-5, but superior to a nine-beam IMRT with MLC-10.

CONCLUSION

The use of an MLC-5 has significant advantages over an MLC-10 with respect to target coverage and protection of normal tissues in step-and-shoot IMRT of head and neck cancer.

Long-term survival of cancer patients compared to heart failure and stroke: a systematic review

BACKGROUND

Cancer, heart failure and stroke are among the most common causes of death worldwide. Investigation of the prognostic impact of each disease is important, especially for a better understanding of competing risks. Aim of this study is to provide an overview of long term survival of cancer, heart failure and stroke patients based on the results of large population- and hospital-based studies.

METHODS

Records for our study were identified by searches of Medline via Pubmed. We focused on observed and relative age- and sex-adjusted 5-year survival rates for cancer in general and for the four most common malignancies in developed countries, i.e. lung, breast, prostate and colorectal cancer, as well as for heart failure and stroke.

RESULTS

Twenty studies were identified and included for analysis. Five-year observed survival was about 43% for all cancer entities, 40-68% for stroke and 26-52% for heart failure. Five-year age and sex adjusted relative survival was 50-57% for all cancer entities, about 50% for stroke and about 62% for heart failure. In regard to the four most common malignancies in developed countries 5-year relative survival was 12-18% for lung cancer, 73-89% for breast cancer, 50-99% for prostate cancer and about 43-63% for colorectal cancer. Trend analysis revealed a survival improvement over the last decades.

CONCLUSIONS

The results indicate that long term survival and prognosis of cancer is not necessarily worse than that of heart failure and stroke. However, a comparison of the prognostic impact of the different diseases is limited, corroborating the necessity for further systematic investigation of competing risks.

Intensity modulated radiotherapy (IMRT) in benign giant cell tumors--a single institution case series and a short review of the literature

BACKGROUND

Giant cell tumors are rare neoplasms, representing less than 5% of all bone tumors. The vast majority of giant cell tumors occurs in extremity sites and is treated by surgery alone. However, a small percentage occurs in pelvis, spine or skull bones, where complete resection is challenging. Radiation therapy seems to be an option in these patients, despite the lack of a generally accepted dose or fractionation concept. Here we present a series of five cases treated with high dose IMRT.

PATIENTS AND METHODS

From 2000 and 2006 a total of five patients with histologically proven benign giant cell tumors have been treated with IMRT in our institution. Two patients were male, three female, and median age was 30 years (range 20-60). The tumor was located in the sacral region in four and in the sphenoid sinus in one patient. All patients had measurable gross disease prior to radiotherapy with a median size of 9 cm. All patients were treated with IMRT to a median total dose of 64 Gy (range 57.6 Gy to 66 Gy) in conventional fractionation.

RESULTS

Median follow up was 46 months ranging from 30 to 107 months. Overall survival was 100%. One patient developed local disease progression three months after radiotherapy and needed extensive surgical salvage. The remaining four patients have been locally controlled, resulting in a local control rate of 80%. We found no substantial tumor shrinkage after radiotherapy but in two patients morphological signs of extensive tumor necrosis were present on MRI scans. Decline of pain and/or neurological symptoms were seen in all four locally controlled patients. The patient who needed surgical salvage showed markedly reduced pain but developed functional deficits of bladder, rectum and lower extremity due to surgery. No severe acute or late toxicities attributable to radiation therapy were observed so far.

CONCLUSION

IMRT is a feasible option in giant cells tumors not amendable to complete surgical removal. In our case series local control was achieved in four out of five patients with marked symptom relief in the majority of cases. No severe toxicity was observed.