Incorporating Cardiac Substructures Into Radiation Therapy For Improved Cardiac Sparing

Growing evidence suggests that radiation therapy (RT) doses to the heart and cardiac substructures (CS) are strongly linked to cardiac toxicities, though only the heart is considered clinically. This work aimed to utilize the superior soft-tissue contrast of magnetic resonance (MR) to segment CS, quantify uncertainties in their position, assess their effect on treatment planning and an MR-guided environment. Automatic substructure segmentation of 12 CS was completed using a novel hybrid MR/computed tomography (CT) atlas method and was improved upon using a 3-dimensional neural network (U-Net) from deep learning. Intra-fraction motion due to respiration was then quantified. The inter-fraction setup uncertainties utilizing a novel MR-linear accelerator were also quantified. Treatment planning comparisons were performed with and without substructure inclusions and methods to reduce radiation dose to sensitive CS were evaluated. Lastly, these described technologies (deep learning U-Net) were translated to an MR-linear accelerator and a segmentation pipeline was created. Automatic segmentations from the hybrid MR/CT atlas was able to generate accurate segmentations for the chambers and great vessels (Dice similarity coefficient (DSC) \u3e 0.75) but coronary artery segmentations were unsuccessful (DSC\u3c0.3). After implementing deep learning, DSC for the chambers and great vessels was ≥0.85 along with an improvement in the coronary arteries (DSC\u3e0.5). Similar accuracy was achieved when implementing deep learning for MR-guided RT. On average, automatic segmentations required ~10 minutes to generate per patient and deep learning only required 14 seconds. The inclusion of CS in the treatment planning process did not yield statistically significant changes in plan complexity, PTV, or OAR dose. Automatic segmentation results from deep learning pose major efficiency and accuracy gains for CS segmentation offering high potential for rapid implementation into radiation therapy planning for improved cardiac sparing. Introducing CS into RT planning for MR-guided RT presented an opportunity for more effective sparing with limited increase in plan complexity

Incorporating Cardiac Substructures Into Radiation Therapy For Improved Cardiac Sparing

Growing evidence suggests that radiation therapy (RT) doses to the heart and cardiac substructures (CS) are strongly linked to cardiac toxicities, though only the heart is considered clinically. This work aimed to utilize the superior soft-tissue contrast of magnetic resonance (MR) to segment CS, quantify uncertainties in their position, assess their effect on treatment planning and an MR-guided environment. Automatic substructure segmentation of 12 CS was completed using a novel hybrid MR/computed tomography (CT) atlas method and was improved upon using a 3-dimensional neural network (U-Net) from deep learning. Intra-fraction motion due to respiration was then quantified. The inter-fraction setup uncertainties utilizing a novel MR-linear accelerator were also quantified. Treatment planning comparisons were performed with and without substructure inclusions and methods to reduce radiation dose to sensitive CS were evaluated. Lastly, these described technologies (deep learning U-Net) were translated to an MR-linear accelerator and a segmentation pipeline was created. Automatic segmentations from the hybrid MR/CT atlas was able to generate accurate segmentations for the chambers and great vessels (Dice similarity coefficient (DSC) \u3e 0.75) but coronary artery segmentations were unsuccessful (DSC\u3c0.3). After implementing deep learning, DSC for the chambers and great vessels was ≥0.85 along with an improvement in the coronary arteries (DSC\u3e0.5). Similar accuracy was achieved when implementing deep learning for MR-guided RT. On average, automatic segmentations required ~10 minutes to generate per patient and deep learning only required 14 seconds. The inclusion of CS in the treatment planning process did not yield statistically significant changes in plan complexity, PTV, or OAR dose. Automatic segmentation results from deep learning pose major efficiency and accuracy gains for CS segmentation offering high potential for rapid implementation into radiation therapy planning for improved cardiac sparing. Introducing CS into RT planning for MR-guided RT presented an opportunity for more effective sparing with limited increase in plan complexity

Quantitative analysis with machine learning models for multi-parametric brain imaging data

Gliomas are considered to be the most common primary adult malignant brain tumor. With the dramatic increases in computational power and improvements in image analysis algorithms, computer-aided medical image analysis has been introduced into clinical applications. Precision tumor grading and genotyping play an indispensable role in clinical diagnosis, treatment and prognosis. Gliomas diagnostic procedures include histopathological imaging tests, molecular imaging scans and tumor grading. Pathologic review of tumor morphology in histologic sections is the traditional method for cancer classification and grading, yet human study has limitations that can result in low reproducibility and inter-observer agreement. Compared with histopathological images, Magnetic resonance (MR) imaging present the different structure and functional features, which might serve as noninvasive surrogates for tumor genotypes. Therefore, computer-aided image analysis has been adopted in clinical application, which might partially overcome these shortcomings due to its capacity to quantitatively and reproducibly measure multilevel features on multi-parametric medical information. Imaging features obtained from a single modal image do not fully represent the disease, so quantitative imaging features, including morphological, structural, cellular and molecular level features, derived from multi-modality medical images should be integrated into computer-aided medical image analysis. The image quality differentiation between multi-modality images is a challenge in the field of computer-aided medical image analysis. In this thesis, we aim to integrate the quantitative imaging data obtained from multiple modalities into mathematical models of tumor prediction response to achieve additional insights into practical predictive value. Our major contributions in this thesis are: 1. Firstly, to resolve the imaging quality difference and observer-dependent in histological image diagnosis, we proposed an automated machine-learning brain tumor-grading platform to investigate contributions of multi-parameters from multimodal data including imaging parameters or features from Whole Slide Images (WSI) and the proliferation marker KI-67. For each WSI, we extract both visual parameters such as morphology parameters and sub-visual parameters including first-order and second-order features. A quantitative interpretable machine learning approach (Local Interpretable Model-Agnostic Explanations) was followed to measure the contribution of features for single case. Most grading systems based on machine learning models are considered “black boxes,” whereas with this system the clinically trusted reasoning could be revealed. The quantitative analysis and explanation may assist clinicians to better understand the disease and accordingly to choose optimal treatments for improving clinical outcomes. 2. Based on the automated brain tumor-grading platform we propose, multimodal Magnetic Resonance Images (MRIs) have been introduced in our research. A new imaging–tissue correlation based approach called RA-PA-Thomics was proposed to predict the IDH genotype. Inspired by the concept of image fusion, we integrate multimodal MRIs and the scans of histopathological images for indirect, fast, and cost saving IDH genotyping. The proposed model has been verified by multiple evaluation criteria for the integrated data set and compared to the results in the prior art. The experimental data set includes public data sets and image information from two hospitals. Experimental results indicate that the model provided improves the accuracy of glioma grading and genotyping

A Semi-Automated Approach to Medical Image Segmentation using Conditional Random Field Inference

Medical image segmentation plays a crucial role in delivering effective patient care in various diagnostic and treatment modalities. Manual delineation of target volumes and all critical structures is a very tedious and highly time-consuming process and introduce uncertainties of treatment outcomes of patients. Fully automatic methods holds great promise for reducing cost and time, while at the same time improving accuracy and eliminating expert variability, yet there are still great challenges. Legally and ethically, human oversight must be integrated with ”smart tools” favoring a semi-automatic technique which can leverage the best aspects of both human and computer. In this work we show that we can formulate a semi-automatic framework for the segmentation problem by formulating it as an energy minimization problem in Conditional Random Field (CRF). We show that human input can be used as adaptive training data to condition a probabilistic boundary term modeled for the heterogeneous boundary characteristics of anatomical structures. We demonstrated that our method can effortlessly adapt to multiple structures and image modalities using a single CRF framework and tools to learn probabilistic terms interactively. To tackle a more difficult multi-class segmentation problem, we developed a new ensemble one-vs-rest graph cut algorithm. Each graph in the ensemble performs a simple and efficient bi-class (a target class vs the rest of the classes) segmentation. The final segmentation is obtained by majority vote. Our algorithm is both faster and more accurate when compared with the prior multi-class method which iteratively swaps classes. In this Thesis, we also include novel volumetric segmentation algorithms which employ deep learning and indicate how to synthesize our CRF framework with convolutional neural networks (CNN). This would allow incorporating user guidance into CNN based deep learning for this task. We think a deep learning based method interactively guided by human expert is the ideal solution for medical image segmentation

Multi-Atlas Segmentation of Biomedical Images: A Survey

Abstract Multi-atlas segmentation (MAS), first introduced and popularized by the pioneering work of Rohlfing

INTEGRATION OF BIOMEDICAL IMAGING AND TRANSLATIONAL APPROACHES FOR MANAGEMENT OF HEAD AND NECK CANCER

The aim of the clinical component of this work was to determine whether the currently available clinical imaging tools can be integrated with radiotherapy (RT) platforms for monitoring and adaptation of radiation dose, prediction of tumor response and disease outcomes, and characterization of patterns of failure and normal tissue toxicity in head and neck cancer (HNC) patients with potentially curable tumors. In Aim 1, we showed that the currently available clinical imaging modalities can be successfully used to adapt RT dose based-on dynamic tumor response, predict oncologic disease outcomes, characterize RT-induced toxicity, and identify the patterns of disease failure. We used anatomical MRIs for the RT dose adaptation purpose. Our findings showed that after proper standardization of the immobilization and image acquisition techniques, we can achieve high geometric accuracy. These images can then be used to monitor the shrinkage of tumors during RT and optimize the clinical target volumes accordingly. Our results also showed that this MR-guided dose adaptation technique has a dosimetric advantage over the standard of care and was associated with a reduction in normal tissue doses that translated into a reduction of the odds of long-term RT-induced toxicity. In the second aim, we used quantitative MRIs to determine its benefit for prediction of oncologic outcomes and characterization of RT-induced normal tissue toxicity. Our findings showed that delta changes of apparent diffusion coefficient parameters derived from diffusion-weighted images at mid-RT can be used to predict local recurrence and recurrence free-survival. We also showed that Ktrans and Ve vascular parameters derived from dynamic contrast-enhanced MRIs can characterize the mandibular areas of osteoradionecrosis. In the final clinical aim, we used CT images of recurrence and baseline CT planning images to develop a methodology and workflow that involves the application of deformable image registration software as a tool to standardize image co-registration in addition to granular combined geometric- and dosimetric-based failure characterization to correctly attribute sites and causes of locoregional failure. We then successfully applied this methodology to identify the patterns of failure following postoperative and definitive IMRT in HNC patients. Using this methodology, we showed that most recurrences occurred in the central high dose regions for patients treated with definitive IMRT compared with mainly non-central high dose recurrences after postoperative IMRT. We also correlated recurrences with pretreatment FDG-PET and identified that most of the central high dose recurrences originated in an area that would be covered by a 10-mm margin on the volume of 50% of the maximum FDG uptake. In the translational component of this work, we integrated radiomic features derived from pre-RT CT images with whole-genome measurements using TCGA and TCIA data. Our results demonstrated a statistically significant associations between radiomic features characterizing different tumor phenotypes and different genomic features. These findings represent a promising potential towards non-invasively tract genomic changes in the tumor during treatment and use this information to adapt treatment accordingly. In the final project of this dissertation, we developed a high-throughput approach to identify effective systemic agents against aggressive head and neck tumors with poor prognosis like anaplastic thyroid cancer. We successfully identified three candidate drugs and performed extensive in vitro and in vivo validation using orthotopic and PDX models. Among these drugs, HDAC inhibitor and LBH-589 showed the most effective tumor growth inhibition that can be used in future clinical trials

Quantifying, Understanding and Predicting Differences Between Planned and Delivered Dose to Organs at Risk in Head & Neck Cancer Patients Undergoing Radical Radiotherapy to Promote Intelligently Targeted Adaptive Radiotherapy

Introduction: Radical radiotherapy (RT) is an effective but toxic treatment for head and neck cancer (HNC). Contemporary radiotherapy techniques sculpt dose to target disease and avoid organs at risk (OARs), but anatomical change during treatment mean that the radiation dose delivered to the patient – delivered dose (DA), is different to that anticipated at planning – planned dose (DP). Modifying the RT plan during treatment – Adaptive Radiotherapy (ART) – could mitigate these risks by reducing dose to OARs. However, clinical data to guide patient selection for, and timing of ART, are for lacking. Methods: 337 patients with HNC were recruited to the Cancer Research UK VoxTox study. Demographic, disease and treatment data were collated, and both DP and DA to organs at risk (OARs) were computed from daily megavoltage CT image guidance scans, using an open-source deformable image registration package (Elastix). Toxicity data were prospectively collected. Relationships between DP, DA and late toxicities were investigated with univariate, and logistic regression normal tissue complication probability (NTCP) modelling approaches. A sub-study of VoxTox recruited 18 patients who had MRI scans before RT fractions 1, 6, 16, and 26. Changes in salivary gland volumes and relative apparent diffusion coefficient (ADC) values were measured and related to toxicity events. Results: Spinal cord dose differences were small, and not predicted by weight loss or shape change. Mean DA to all other OARs was higher than DP; factors predicting higher DA included primary disease site, concomitant therapy, shape change and advanced neck disease. Nine patients (3.7%) saw DA>DP by 2Gy to more than half of the OARs assessed. These patients all had received bilateral neck RT for N-stage 2b oropharyngeal cancer. Strong uni- and multivariate relationships between OAR dose and toxicity were seen. Differences between DA and DP-based dose-toxicity models were minimal, and not statistically significant. On MRI, both parotid and submandibular glands shrank during treatment, whilst relative ADC rose. Relationships with toxicity were inconclusive. Conclusions: Small differences between OAR DP and DA mean that DA-based toxicity prediction models confer negligible additional benefit at the population level. Factors such as primary disease sub-site, concomitant systemic therapy, staging and shape change may help to select the patients that do develop clinically significant dose differences, and would benefit most from ART for toxicity reduction

Development of computer-based algorithms for unsupervised assessment of radiotherapy contouring

INTRODUCTION: Despite the advances in radiotherapy treatment delivery, target volume delineation remains one of the greatest sources of error in the radiotherapy delivery process, which can lead to poor tumour control probability and impact clinical outcome. Contouring assessments are performed to ensure high quality of target volume definition in clinical trials but this can be subjective and labour-intensive. This project addresses the hypothesis that computational segmentation techniques, with a given prior, can be used to develop an image-based tumour delineation process for contour assessments. This thesis focuses on the exploration of the segmentation techniques to develop an automated method for generating reference delineations in the setting of advanced lung cancer. The novelty of this project is in the use of the initial clinician outline as a prior for image segmentation. METHODS: Automated segmentation processes were developed for stage II and III non-small cell lung cancer using the IDEAL-CRT clinical trial dataset. Marker-controlled watershed segmentation, two active contour approaches (edge- and region-based) and graph-cut applied on superpixels were explored. k-nearest neighbour (k-NN) classification of tumour from normal tissues based on texture features was also investigated. RESULTS: 63 cases were used for development and training. Segmentation and classification performance were evaluated on an independent test set of 16 cases. Edge-based active contour segmentation achieved highest Dice similarity coefficient of 0.80 ± 0.06, followed by graphcut at 0.76 ± 0.06, watershed at 0.72 ± 0.08 and region-based active contour at 0.71 ± 0.07, with mean computational times of 192 ± 102 sec, 834 ± 438 sec, 21 ± 5 sec and 45 ± 18 sec per case respectively. Errors in accuracy of irregularly shaped lesions and segmentation leakages at the mediastinum were observed. In the distinction of tumour and non-tumour regions, misclassification errors of 14.5% and 15.5% were achieved using 16- and 8-pixel regions of interest (ROIs) respectively. Higher misclassification errors of 24.7% and 26.9% for 16- and 8-pixel ROIs were obtained in the analysis of the tumour boundary. CONCLUSIONS: Conventional image-based segmentation techniques with the application of priors are useful in automatic segmentation of tumours, although further developments are required to improve their performance. Texture classification can be useful in distinguishing tumour from non-tumour tissue, but the segmentation task at the tumour boundary is more difficult. Future work with deep-learning segmentation approaches need to be explored.Funded by National Radiotherapy Trials Quality Assurance (RTTQA) grou