The NCS code of practice for the quality assurance of treatment planning systems (NCS-35)

A subcommittee of the Netherlands Commission on Radiation Dosimetry (NCS) was initiated in 2018 with the task to update and extend a previous publication (NCS-15) on the quality assurance of treatment planning systems (TPS) (Bruinviset al2005). The field of treatment planning has changed considerably since 2005. Whereas the focus of the previous report was more on the technical aspects of the TPS, the scope of this report is broader with a focus on a department wide implementation of the TPS. New sections about education, automated planning, information technology (IT) and updates are therefore added. Although the scope is photon therapy, large parts of this report will also apply to all other treatment modalities. This paper is a condensed version of these guidelines; the full version of the report in English is freely available from the NCS website (http://radiationdosimetry.org/ncs/publications). The paper starts with the scope of this report in relation to earlier reports on this subject. Next, general aspects of the commissioning process are addressed, like e.g. project management, education, and safety. It then focusses more on technical aspects such as beam commissioning and patient modeling, dose representation, dose calculation and (automated) plan optimisation. The final chapters deal with IT-related subjects and scripting, and the process of updating or upgrading the TPS.

Validation of skeletal muscle mass assessment at the level of the third cervical vertebra in patients with head and neck cancer

BACKGROUND

Low skeletal muscle mass (SMM) is associated with adverse outcomes. SMM is often assessed at the third lumbar vertebra (L3) on abdominal imaging. Abdominal imaging is not routinely performed in patients with head and neck cancer (HNC). We aim to validate SMM measurement at the level of the third cervical vertebra (C3) on head and neck imaging.

MATERIAL AND METHODS

Patients with pre-treatment whole-body computed tomography (CT) between 2010 and 2018 were included. Cross-sectional muscle area (CSMA) was manually delineated at the level of C3 and L3. Correlation coefficients and intraclass correlation coefficients (ICCs) were calculated. Cohen's kappa was used to assess the reliability of identifying a patient with low SMM.

RESULTS

Two hundred patients were included. Correlation between CSMA at the level of C3 and L3 was good (r = 0.75, p < 0.01). Using a multivariate formula to estimate CSMA at L3, including gender, age, and weight, correlation improved (r = 0.82, p < 0.01). The agreement between estimated and actual CSMA at L3 was good (ICC 0.78, p < 0.01). There was moderate agreement in the identification of patients with low SMM based on the estimated lumbar skeletal muscle mass index (LSMI) and actual LSMI (Cohen's κ: 0.57, 95%CI 0.45-0.69).

CONCLUSIONS

CSMA at C3 correlates well with CSMA at L3. There is moderate agreement in the identification of patients with low SMM based on the estimated lumbar SMI (based on measurement at C3) and actual LSMI.

DeepDose: a robust deep learning-based dose engine for abdominal tumours in a 1.5 T MRI radiotherapy system

We present a robust deep learning-based framework for dose calculations of abdominal tumours in a 1.5 T MRI radiotherapy system. For a set of patient plans, a convolutional neural network is trained on the dose of individual multi-leaf-collimator segments following the DeepDose framework. It can then be used to predict the dose distribution per segment for a set of patient anatomies. The network was trained using data from three anatomical sites of the abdomen: prostate, rectal and oligometastatic tumours. A total of 216 patient fractions were used, previously treated in our clinic with fixed-beam IMRT using the Elekta MR-linac. For the purpose of training, 176 fractions were used with random gantry angles assigned to each segment, while 20 fractions were used for the validation of the network. The ground truth data were calculated with a Monte Carlo dose engine at 1% statistical uncertainty per segment. For a total of 20 independent abdominal test fractions with the clinical angles, the network was able to accurately predict the dose distributions, achieving 99.4% ± 0.6% for the whole plan prediction at the 3%/3 mm gamma test. The average dose difference and standard deviation per segment was 0.3% ± 0.7%. Additional dose prediction on one cervical and one pancreatic case yielded high dose agreement of 99.9% and 99.8% respectively for the 3%/3 mm criterion. Overall, we show that our deep learning-based dose engine calculates highly accurate dose distributions for a variety of abdominal tumour sites treated on the MR-linac, in terms of performance and generality.

Proof-of-concept delivery of intensity modulated arc therapy on the Elekta Unity 1.5 T MR-linac

In this work we present the first delivery of intensity modulated arc therapy on the Elekta Unity 1.5 T MR-linac. The machine's current intensity modulated radiation therapy based control system was modified suitably to enable dynamic delivery of radiation, for the purpose of exploring MRI-guided radiation therapy adaptation modes in a research setting. The proof-of-concept feasibility was demonstrated by planning and delivering two types of plans, each investigating the performance of different parts of a dynamic treatment. A series of fixed-speed arc plans was used to show the high-speed capabilities of the gantry during radiation, while several fully modulated prostate plans-optimised following the volumetric modulated arc therapy approach-were delivered in order to establish the performance of its multi-leaf collimator and diaphragms. These plans were delivered to Delta⁴ Phantom+ MR and film phantoms passing the clinical quality assurance criteria used in our clinic. In addition, we also performed some initial MR imaging experiments during dynamic therapy, demonstrating that the impact of radiation and moving gantry/collimator components on the image quality is negligible. These results show that arc therapy is feasible on the Elekta Unity system. The machine's high performance components enable dynamic delivery during fast gantry rotation and can be controlled in a stable fashion to deliver fully modulated plans.

Impact of body fat distribution and sarcopenia on the overall survival in patients with spinal metastases receiving radiotherapy treatment: a prospective cohort study

Introduction: An increasing number of patients is diagnosed with spinal metastases due to elevated cancer incidence and improved overall survival. Patients with symptomatic spinal bone metastases often receive radiotherapy with or without surgical stabilisation. Patients with a life expectancy of less than 3 months are generally deemed unfit for surgery, therefore adequate pre-treatment assessment of life expectancy is necessary. The aim of this study was to assess new factors associated with overall survival for this category of patients.Patients and methods: Patients who received radiotherapy for thoracic or lumbar spinal metastases from June 2013 to December 2016 were included in this study. The pre-treatment planning CT for radiotherapy treatment was used to assess the patient's visceral fat area, subcutaneous fat area, total muscle area and skeletal muscle density on a single transverse slice at the L3 level. The total muscle area was used to assess sarcopenia. Furthermore, data were collected on age, sex, primary tumour, Karnofsky performance score, medical history, number of bone metastases, non-bone metastases and neurological symptoms. Univariable and multivariable cox regressions were performed to determine the association between our variables of interest and the survival at 90 and 365 days.Results: A total of 310 patients was included. The median age was 67 years. Overall survival rates for 90 and 365 days were 71% and 36% respectively. For 90- and 365-day survival, the Karnofsky performance score, muscle density and primary tumour were independently significantly associated. The visceral or subcutaneous fat area and their ratio and sarcopenia were not independently associated with overall survival.Conclusions: Of the body morphology, only muscle density was statistically significant associated with overall survival after 90 and 365 days in patients with spinal bone metastases. Body fat distribution was not significantly associated with overall survival.

Interobserver agreement of skeletal muscle mass measurement on head and neck CT imaging at the level of the third cervical vertebra

OBJECTIVES

Skeletal muscle mass (SMM) is most often assessed in cancer patients on abdominal computed tomography (CT) imaging at the level of the third lumbar vertebra (L3). Abdominal CT imaging is not routinely performed in head and neck cancer (HNC) patients. Recently, a novel method to assess SMM on a single transversal CT slice at the level of the third cervical vertebra (C3) was published. The objective of this study was to assess the robustness of this novel C3 measurement method in terms of interobserver agreement.

PATIENTS AND METHODS

Patients diagnosed with locally advanced head and neck squamous cell carcinoma (LA-HNSCC) at our center between 2007 and 2011 were evaluated. Fifty-four patients with were randomly selected for analysis. Six observers independently measured the cross-sectional muscle area (CSMA) at the level of C3 using a predefined, written protocol as instruction. Interobserver agreement was assessed using intraclass correlation coefficients (ICCs), a Bland-Altman plot and Fleiss' kappa (κ).

RESULTS

The agreement in vertebra selection between all observers was excellent (Fleiss' κ: 0.96). There was a substantial agreement between all observers in single slice selection (Fleiss' κ: 0.61). For all CSMA measurements, ICCs were excellent (0.763-0.969; all p < 0.001). The Bland-Altman plot showed good agreement between measurements, with narrow limits of agreement.

CONCLUSION

Interobserver agreement for SMM measurement at the level of C3 was excellent. Assessment of SMM at the level of C3 is easy and robust and can performed on routinely available imaging in HNC patients.

Towards fast online intrafraction replanning for free-breathing stereotactic body radiation therapy with the MR-linac

The hybrid MRI-radiotherapy machines, like the MR-linac (Elekta AB, Stockholm, Sweden) installed at the UMC Utrecht (Utrecht, The Netherlands), will be able to provide real-time patient imaging during treatment. In order to take advantage of the system's capabilities and enable online adaptive treatments, a new generation of software should be developed, ranging from motion estimation to treatment plan adaptation. In this work we present a proof of principle adaptive pipeline designed for high precision stereotactic body radiation therapy (SBRT) suitable for sites affected by respiratory motion, like renal cell carcinoma (RCC). We utilized our research MRL treatment planning system (MRLTP) to simulate a single fraction 25 Gy free-breathing SBRT treatment for RCC by performing inter-beam replanning for two patients and one volunteer. The simulated pipeline included a combination of (pre-beam) 4D-MRI and (online) 2D cine-MR acquisitions. The 4DMRI was used to generate the mid-position reference volume, while the cine-MRI, via an in-house motion model, provided three-dimensional (3D) deformable vector fields (DVFs) describing the anatomical changes during treatment. During the treatment fraction, at an inter-beam interval, the mid-position volume of the patient was updated and the delivered dose was accurately reconstructed on the underlying motion calculated by the model. Fast online replanning, targeting the latest anatomy and incorporating the previously delivered dose was then simulated with MRLTP. The adaptive treatment was compared to a conventional mid-position SBRT plan with a 3 mm planning target volume margin reconstructed on the same motion trace. We demonstrate that our system produced tighter dose distributions and thus spared the healthy tissue, while delivering more dose to the target. The pipeline was able to account for baseline variations/drifts that occurred during treatment ensuring target coverage at the end of the treatment fraction.

Development and clinical introduction of automated radiotherapy treatment planning for prostate cancer

To develop an automated radiotherapy treatment planning and optimization workflow to efficiently create patient specifically optimized clinical grade treatment plans for prostate cancer and to implement it in clinical practice. A two-phased planning and optimization workflow was developed to automatically generate 77Gy 5-field simultaneously integrated boost intensity modulated radiation therapy (SIB-IMRT) plans for prostate cancer treatment. A retrospective planning study (n  =  100) was performed in which automatically and manually generated treatment plans were compared. A clinical pilot (n  =  21) was performed to investigate the usability of our method. Operator time for the planning process was reduced to  <5 min. The retrospective planning study showed that 98 plans met all clinical constraints. Significant improvements were made in the volume receiving 72Gy (V72Gy) for the bladder and rectum and the mean dose of the bladder and the body. A reduced plan variance was observed. During the clinical pilot 20 automatically generated plans met all constraints and 17 plans were selected for treatment. The automated radiotherapy treatment planning and optimization workflow is capable of efficiently generating patient specifically optimized and improved clinical grade plans. It has now been adopted as the current standard workflow in our clinic to generate treatment plans for prostate cancer.

A new methodology for inter- and intrafraction plan adaptation for the MR-linac

The new era of hybrid MRI and linear accelerator machines, including the MR-linac currently being installed in the University Medical Center Utrecht (Utrecht, The Netherlands), will be able to provide the actual anatomy and real-time anatomy changes of the patient's target(s) and organ(s) at risk (OARs) during radiation delivery. In order to be able to take advantage of this input, a new generation of treatment planning systems is needed, that will allow plan adaptation to the latest anatomy state in an online regime. In this paper, we present a treatment planning algorithm for intensity-modulated radiotherapy (IMRT), which is able to compensate for patient anatomy changes. The system consists of an iterative sequencing loop open to anatomy updates and an inter- and intrafraction adaptation scheme that enables convergence to the ideal dose distribution without the need of a final segment weight optimization (SWO). The ability of the system to take into account organ motion and adapt the plan to the latest anatomy state is illustrated using artificial baseline shifts created for three different kidney cases. Firstly, for two kidney cases of different target volumes, we show that the system can account for intrafraction motion, delivering the intended dose to the target with minimal dose deposition to the surroundings compared to conventional plans. Secondly, for a third kidney case we show that our algorithm combined with the interfraction scheme can be used to deliver the prescribed dose while adapting to the changing anatomy during multi-fraction treatments without performing a final SWO.

Patient position verification using small IMRT fields

A commonly used approach to quantify and minimize patient setup errors is by using electronic portal imaging devices (EPIDs). The position of the tumor can be verified indirectly by matching the bony anatomy to a reference image containing the same structures. In this paper we present two off-line methods for detecting the position of the bony anatomy automatically, even if every single portal image of each segment of an IMRT treatment beam contains insufficient matching information. Extra position verification fields will no longer be necessary, which reduces the total dose to the patient. The first method, the stack matching method (SMM), stacks the portal image of each segment of a beam to a three dimensional (3D) volume, and this volume is subsequently used during the matching phase. The second method [the averaged projection matching method (APMM)], is a simplification of the first one, since the initially created volume is reduced again to a 2D artificial image, which speeds up the matching procedure considerably, without a significant loss of accuracy. Matching is based on normalized mutual information. We demonstrate our methods by comparing them to existing matching routines, such as matching based on the largest segment. Both phantom and patient experiments show that our methods are comparable with the results obtained from standard position verification methods. The matches are verified by means of visual inspection. Furthermore, we show that when a distinct area of 40-60 cm2 of the EPID is exposed during one treatment beam, both SMM and APMM are able to deliver a good matching result.

Assessment of rigid multi-modality image registration consistency using the multiple sub-volume registration (MSR) method

Registration of different imaging modalities such as CT, MRI, functional MRI (fMRI), positron (PET) and single photon (SPECT) emission tomography is used in many clinical applications. Determining the quality of any automatic registration procedure has been a challenging part because no gold standard is available to evaluate the registration. In this note we present a method, called the 'multiple sub-volume registration' (MSR) method, for assessing the consistency of a rigid registration. This is done by registering sub-images of one data set on the other data set, performing a crude non-rigid registration. By analysing the deviations (local deformations) of the sub-volume registrations from the full registration we get a measure of the consistency of the rigid registration. Registration of 15 data sets which include CT, MR and PET images for brain, head and neck, cervix, prostate and lung was performed utilizing a rigid body registration with normalized mutual information as the similarity measure. The resulting registrations were classified as good or bad by visual inspection. The resulting registrations were also classified using our MSR method. The results of our MSR method agree with the classification obtained from visual inspection for all cases (p < 0.02 based on ANOVA of the good and bad groups). The proposed method is independent of the registration algorithm and similarity measure. It can be used for multi-modality image data sets and different anatomic sites of the patient.