Towards patient-specific modeling of brain tumor growth and formation of secondary nodes guided by DTI-MRI

Using mathematical modelling and AI to improve delivery and efficacy of therapies in cancer

Mathematical modelling has proven to be a valuable tool in predicting the delivery and efficacy of molecular, antibody-based, nano and cellular therapy in solid tumours. Mathematical models based on our understanding of the biological processes at subcellular, cellular and tissue level are known as mechanistic models that, in turn, are divided into continuous and discrete models. Continuous models are further divided into lumped parameter models - for describing the temporal distribution of medicine in tumours and normal organs - and distributed parameter models - for studying the spatiotemporal distribution of therapy in tumours. Discrete models capture interactions at the cellular and subcellular levels. Collectively, these models are useful for optimizing the delivery and efficacy of molecular, nanoscale and cellular therapy in tumours by incorporating the biological characteristics of tumours, the physicochemical properties of drugs, the interactions among drugs, cancer cells and various components of the tumour microenvironment, and for enabling patient-specific predictions when combined with medical imaging. Artificial intelligence-based methods, such as machine learning, have ushered in a new era in oncology. These data-driven approaches complement mechanistic models and have immense potential for improving cancer detection, treatment and drug discovery. Here we review these diverse approaches and suggest ways to combine mechanistic and artificial intelligence-based models to further improve patient treatment outcomes.

Chemotaxis effects on the vascular tumor growth: Phase-field model and simulations

In this paper, we propose a vascular tumor growth model that combines a phase-field tumor model with a phase-field angiogenesis model. By incorporating various tumor cell species, we capture the instabilities of the tumor in the presence of evolving neovasculature. The model not only considers different dynamics of tumor cell phase conversions, movement, and pressure effects but also provides a comprehensive representation of angiogenesis, encompassing chemotaxis of endothelial cells, sprouting, anastomoses, and blood flow in capillaries. This study evaluates the impact of chemotaxis on tumor cell movement in both avascular and vascular tumor growth scenarios. The results highlight the acceleration of tumor growth when angiogenesis is stimulated. Additionally, the investigation explores various initial distances of the tumor from neighboring vessels, revealing a critical threshold distance beyond which the angiogenesis factor fails to stimulate angiogenesis, resulting in the tumor maintaining a stable state. The integration of chemotaxis into the growth model induces instabilities, leading to increased nutrient availability and faster growth for the tumor. Furthermore, the study considers anti-angiogenesis therapy as an ideal approach, assuming complete inhibition of angiogenesis from the early stages. In this scenario, the tumor persists in a steady state, adhering to the avascular size limit in the absence of neovasculature. Conversely, when considering chemotaxis, anti-angiogenesis therapy loses efficiency, enabling unrestrained tumor growth towards neighboring vessels. This work sheds light on the intricate interplay among chemotaxis, angiogenesis, and anti-angiogenesis therapy in the context of vascular tumor growth, providing valuable insights for the development of targeted treatment strategies.

Heterogeneous Mechanical Stress and Interstitial Fluid Flow Predictions Derived from DCE-MRI for Rat U251N Orthotopic Gliomas

Mechanical stress and fluid flow influence glioma cell phenotype in vitro, but measuring these quantities in vivo continues to be challenging. The purpose of this study was to predict these quantities in vivo, thus providing insight into glioma physiology and potential mechanical biomarkers that may improve glioma detection, diagnosis, and treatment. Image-based finite element models of human U251N orthotopic glioma in athymic rats were developed to predict structural stress and interstitial flow in and around each animal's tumor. In addition to accounting for structural stress caused by tumor growth, our approach has the advantage of capturing fluid pressure-induced structural stress, which was informed by in vivo interstitial fluid pressure (IFP) measurements. Because gliomas and the brain are soft, elevated IFP contributed substantially to tumor structural stress, even inverting this stress from compressive to tensile in the most compliant cases. The combination of tumor growth and elevated IFP resulted in a concentration of structural stress near the tumor boundary where it has the greatest potential to influence cell proliferation and invasion. MRI-derived anatomical geometries and tissue property distributions resulted in heterogeneous interstitial fluid flow with local maxima near cerebrospinal fluid spaces, which may promote tumor invasion and hinder drug delivery. In addition, predicted structural stress and interstitial flow varied markedly between irradiated and radiation-naïve animals. Our modeling suggests that relative to tumors in stiffer tissues, gliomas experience unusual mechanical conditions with potentially important biological (e.g., proliferation and invasion) and clinical consequences (e.g., drug delivery and treatment monitoring).

Patient-Specific, Mechanistic Models of Tumor Growth Incorporating Artificial Intelligence and Big Data

Despite the remarkable advances in cancer diagnosis, treatment, and management over the past decade, malignant tumors remain a major public health problem. Further progress in combating cancer may be enabled by personalizing the delivery of therapies according to the predicted response for each individual patient. The design of personalized therapies requires the integration of patient-specific information with an appropriate mathematical model of tumor response. A fundamental barrier to realizing this paradigm is the current lack of a rigorous yet practical mathematical theory of tumor initiation, development, invasion, and response to therapy. We begin this review with an overview of different approaches to modeling tumor growth and treatment, including mechanistic as well as data-driven models based on big data and artificial intelligence. We then present illustrative examples of mathematical models manifesting their utility and discuss the limitations of stand-alone mechanistic and data-driven models. We then discuss the potential of mechanistic models for not only predicting but also optimizing response to therapy on a patient-specific basis. We describe current efforts and future possibilities to integrate mechanistic and data-driven models. We conclude by proposing five fundamental challenges that must be addressed to fully realize personalized care for cancer patients driven by computational models.

Investigating the synergistic effects of immunotherapy and normalization treatment in modulating tumor microenvironment and enhancing treatment efficacy

We developed a comprehensive mathematical model of cancer immunotherapy that takes into account: i) Immune checkpoint blockers (ICBs) and the interactions between cancer cells and the immune system, ii) characteristics of the tumor microenvironment, such as the tumor hydraulic conductivity, interstitial fluid pressure, and vascular permeability, iii) spatial and temporal variations of the modeled components within the tumor and the surrounding host tissue, iv) the transport of modeled components through the vasculature and between the tumor-host tissue with convection and diffusion, and v) modeling of the tumor draining lymph nodes were the antigen presentation and the development of cytotoxic immune cells take place. Our model successfully reproduced experimental data from various murine tumor types and predicted immune system profiling, which is challenging to achieve experimentally. It showed that combination of ICB therapy and normalization treatments, that aim to improve tumor perfusion, decreases interstitial fluid pressure and increases the concentration of both innate and adaptive immune cells at the tumor center rather than the periphery. Furthermore, using the model, we investigated the impact of modeled components on treatment outcomes. The analysis found that the number of functional vessels inside the tumor region and the ICB dose administered have the largest impact on treatment outcomes.

Applications of generative adversarial networks in neuroimaging and clinical neuroscience

Generative adversarial networks (GANs) are one powerful type of deep learning models that have been successfully utilized in numerous fields. They belong to the broader family of generative methods, which learn to generate realistic data with a probabilistic model by learning distributions from real samples. In the clinical context, GANs have shown enhanced capabilities in capturing spatially complex, nonlinear, and potentially subtle disease effects compared to traditional generative methods. This review critically appraises the existing literature on the applications of GANs in imaging studies of various neurological conditions, including Alzheimer's disease, brain tumors, brain aging, and multiple sclerosis. We provide an intuitive explanation of various GAN methods for each application and further discuss the main challenges, open questions, and promising future directions of leveraging GANs in neuroimaging. We aim to bridge the gap between advanced deep learning methods and neurology research by highlighting how GANs can be leveraged to support clinical decision making and contribute to a better understanding of the structural and functional patterns of brain diseases.

Data-driven spatio-temporal modelling of glioblastoma

Mathematical oncology provides unique and invaluable insights into tumour growth on both the microscopic and macroscopic levels. This review presents state-of-the-art modelling techniques and focuses on their role in understanding glioblastoma, a malignant form of brain cancer. For each approach, we summarize the scope, drawbacks and assets. We highlight the potential clinical applications of each modelling technique and discuss the connections between the mathematical models and the molecular and imaging data used to inform them. By doing so, we aim to prime cancer researchers with current and emerging computational tools for understanding tumour progression. By providing an in-depth picture of the different modelling techniques, we also aim to assist researchers who seek to build and develop their own models and the associated inference frameworks. Our article thus strikes a unique balance. On the one hand, we provide a comprehensive overview of the available modelling techniques and their applications, including key mathematical expressions. On the other hand, the content is accessible to mathematicians and biomedical scientists alike to accommodate the interdisciplinary nature of cancer research.

Coupling solid and fluid stresses with brain tumour growth and white matter tract deformations in a neuroimaging-informed model

Brain tumours are among the deadliest types of cancer, since they display a strong ability to invade the surrounding tissues and an extensive resistance to common therapeutic treatments. It is therefore important to reproduce the heterogeneity of brain microstructure through mathematical and computational models, that can provide powerful instruments to investigate cancer progression. However, only a few models include a proper mechanical and constitutive description of brain tissue, which instead may be relevant to predict the progression of the pathology and to analyse the reorganization of healthy tissues occurring during tumour growth and, possibly, after surgical resection. Motivated by the need to enrich the description of brain cancer growth through mechanics, in this paper we present a mathematical multiphase model that explicitly includes brain hyperelasticity. We find that our mechanical description allows to evaluate the impact of the growing tumour mass on the surrounding healthy tissue, quantifying the displacements, deformations, and stresses induced by its proliferation. At the same time, the knowledge of the mechanical variables may be used to model the stress-induced inhibition of growth, as well as to properly modify the preferential directions of white matter tracts as a consequence of deformations caused by the tumour. Finally, the simulations of our model are implemented in a personalized framework, which allows to incorporate the realistic brain geometry, the patient-specific diffusion and permeability tensors reconstructed from imaging data and to modify them as a consequence of the mechanical deformation due to cancer growth.

Network analysis shows decreased ipsilesional structural connectivity in glioma patients

Gliomas that infiltrate networks and systems, such as the motor system, often lead to substantial functional impairment in multiple systems. Network-based statistics (NBS) allow to assess local network differences and graph theoretical analyses enable investigation of global and local network properties. Here, we used network measures to characterize glioma-related decreases in structural connectivity by comparing the ipsi- with the contralesional hemispheres of patients and correlated findings with neurological assessment. We found that lesion location resulted in differential impairment of both short and long connectivity patterns. Network analysis showed reduced global and local efficiency in the ipsilesional hemisphere compared to the contralesional hemispheric networks, which reflect the impairment of information transfer across different regions of a network.

Tumors and their location distinctly alter both local and global brain connectivity within the ipsilesional hemisphere of glioma patients.

Inducing Biomechanical Heterogeneity in Brain Tumor Modeling by MR Elastography: Effects on Tumor Growth, Vascular Density and Delivery of Therapeutics

Simple Summary

Biomechanical forces aggravate brain tumor progression. In this study, magnetic resonance elastography (MRE) is employed to extract tissue biomechanical properties from five glioblastoma patients and a healthy subject, and data are incorporated in a mathematical model that simulates tumor growth. Mathematical modeling enables further understanding of glioblastoma development and allows patient-specific predictions for tumor vascularity and delivery of drugs. Incorporating MRE data results in a more realistic intratumoral distribution of mechanical stress and anisotropic tumor growth and a better description of subsequent events that are closely related to the development of stresses, including heterogeneity of the tumor vasculature and intrapatient variations in tumor perfusion and delivery of drugs.

Abstract

The purpose of this study is to develop a methodology that incorporates a more accurate assessment of tissue mechanical properties compared to current mathematical modeling by use of biomechanical data from magnetic resonance elastography. The elastography data were derived from five glioblastoma patients and a healthy subject and used in a model that simulates tumor growth, vascular changes due to mechanical stresses and delivery of therapeutic agents. The model investigates the effect of tumor-specific biomechanical properties on tumor anisotropic growth, vascular density heterogeneity and chemotherapy delivery. The results showed that including elastography data provides a more realistic distribution of the mechanical stresses in the tumor and induces anisotropic tumor growth. Solid stress distribution differs among patients, which, in turn, induces a distinct functional vascular density distribution—owing to the compression of tumor vessels—and intratumoral drug distribution for each patient. In conclusion, incorporating elastography data results in a more accurate calculation of intratumoral mechanical stresses and enables a better mathematical description of subsequent events, such as the heterogeneous development of the tumor vasculature and intrapatient variations in tumor perfusion and delivery of drugs.

A computational model of glioma reveals opposing, stiffness-sensitive effects of leaky vasculature and tumor growth on tissue mechanical stress and porosity

A biphasic computational model of a growing, vascularized glioma within brain tissue was developed to account for unique features of gliomas, including soft surrounding brain tissue, their low stiffness relative to brain tissue, and a lack of draining lymphatics. This model is the first to couple nonlinear tissue deformation with porosity and tissue hydraulic conductivity to study the mechanical interaction of leaky vasculature and solid growth in an embedded glioma. The present model showed that leaky vasculature and elevated interstitial fluid pressure produce tensile stress within the tumor in opposition to the compressive stress produced by tumor growth. This tensile effect was more pronounced in softer tissue and resulted in a compressive stress concentration at the tumor rim that increased when tumor was softer than host. Aside from generating solid stress, fluid pressure-driven tissue deformation decreased the effective stiffness of the tumor while growth increased it, potentially leading to elevated stiffness in the tumor rim. A novel prediction of reduced porosity at the tumor rim was corroborated by direct comparison with estimates from our in vivo imaging studies. Antiangiogenic and radiation therapy were simulated by varying vascular leakiness and tissue hydraulic conductivity. These led to greater solid compression and interstitial pressure in the tumor, respectively, the former of which may promote tumor infiltration of the host. Our findings suggest that vascular leakiness has an important influence on in vivo solid stress, stiffness, and porosity fields in gliomas given their unique mechanical microenvironment.

Imaging acute effects of bevacizumab on tumor vascular kinetics in a preclinical orthotopic model of U251 glioma

The effect of a human vascular endothelial growth factor antibody on the vasculature of human tumor grown in rat brain was studied. Using dynamic contrast-enhanced magnetic resonance imaging, the effects of intravenous bevacizumab (Avastin; 10 mg/kg) were examined before and at postadministration times of 1, 2, 4, 8, 12 and 24 h (N = 26; 4-5 per time point) in a rat model of orthotopic, U251 glioblastoma (GBM). The commonly estimated vascular parameters for an MR contrast agent were: (i) plasma distribution volume (v_p ), (ii) forward volumetric transfer constant (K^trans ) and (iii) reverse transfer constant (k_ep ). In addition, extracellular distribution volume (V_D ) was estimated in the tumor (V_D-tumor ), tumor edge (V_D-edge ) and the mostly normal tumor periphery (V_D-peri ), along with tumor blood flow (TBF), peri-tumoral hydraulic conductivity (K) and interstitial flow (Flux) and tumor interstitial fluid pressure (TIFP). Studied as % changes from baseline, the 2-h post-treatment time point began showing significant decreases in v_p , V_D-tumor, V_D-edge and V_D-peri , as well as K, with these changes persisting at 4 and 8 h in v_p , K, V_{D-tumor, -edge} and _-peri (t-tests; p < 0.05-0.01). Decreases in K^trans were observed at the 2- and 4-h time points (p < 0.05), while interstitial volume fraction (v_e ; = K^trans /k_ep ) showed a significant decrease only at the 2-h time point (p < 0.05). Sustained decreases in Flux were observed from 2 to 24 h (p < 0.01) while TBF and TIFP showed delayed responses, increases in the former at 12 and 24 h and a decrease in the latter only at 12 h. These imaging biomarkers of tumor vascular kinetics describe the short-term temporal changes in physical spaces and fluid flows in a model of GBM after Avastin administration.

In Silico Mathematical Modelling for Glioblastoma: A Critical Review and a Patient-Specific Case

Glioblastoma extensively infiltrates the brain; despite surgery and aggressive therapies, the prognosis is poor. A multidisciplinary approach combining mathematical, clinical and radiological data has the potential to foster our understanding of glioblastoma evolution in every single patient, with the aim of tailoring therapeutic weapons. In particular, the ultimate goal of biomathematics for cancer is the identification of the most suitable theoretical models and simulation tools, both to describe the biological complexity of carcinogenesis and to predict tumor evolution. In this report, we describe the results of a critical review about different mathematical models in neuro-oncology with their clinical implications. A comprehensive literature search and review for English-language articles concerning mathematical modelling in glioblastoma has been conducted. The review explored the different proposed models, classifying them and indicating the significative advances of each one. Furthermore, we present a specific case of a glioblastoma patient in which our recently proposed innovative mechanical model has been applied. The results of the mathematical models have the potential to provide a relevant benefit for clinicians and, more importantly, they might drive progress towards improving tumor control and patient's prognosis. Further prospective comparative trials, however, are still necessary to prove the impact of mathematical neuro-oncology in clinical practice.

Feasibility and clinical usefulness of modelling glioblastoma migration in adjuvant radiotherapy

Glioblastoma (GBM) is one of the most common primary brain tumours in adults, with a dismal prognosis despite aggressive multimodality treatment by a combination of surgery and adjuvant radiochemotherapy. A detailed knowledge of the spreading of glioma cells in the brain might allow for more targeted escalated radiotherapy, aiming to reduce locoregional relapse. Recent years have seen the development of a large variety of mathematical modelling approaches to predict glioma migration. The aim of this study is hence to evaluate the clinical applicability of a detailed micro- and meso-scale mathematical model in radiotherapy. First and foremost, a clinical workflow is established, in which the tumour is automatically segmented as input data and then followed in time mathematically based on the diffusion tensor imaging data. The influence of several free model parameters is individually evaluated, then the full model is retrospectively validated for a collective of 3 GBM patients treated at our institution by varying the most important model parameters to achieve optimum agreement with the tumour development during follow-up. Agreement of the model predictions with the real tumour growth as defined by manual contouring based on the follow-up MRI images is analyzed using the dice coefficient. The tumour evolution over 103-212 days follow-up could be predicted by the model with a dice coefficient better than 60% for all three patients. In all cases, the final tumour volume was overestimated by the model by a factor between 1.05 and 1.47. To evaluate the quality of the agreement between the model predictions and the ground truth, we must keep in mind that our gold standard relies on a single observer's (CB) manually-delineated tumour contours. We therefore decided to add a short validation of the stability and reliability of these contours by an inter-observer analysis including three other experienced radiation oncologists from our department. In total, a dice coefficient between 63% and 89% is achieved between the four different observers. Compared with this value, the model predictions (62-66%) perform reasonably well, given the fact that these tumour volumes were created based on the pre-operative segmentation and DTI.

Retrospective analysis of hemispheric structural network change as a function of location and size of glioma

Gliomas are neoplasms that arise from glial cell origin and represent the largest fraction of primary malignant brain tumours (77%). These highly infiltrative malignant cell clusters modify brain structure and function through expansion, invasion and intratumoral modification. Depending on the growth rate of the tumour, location and degree of expansion, functional reorganization may not lead to overt changes in behaviour despite significant cerebral adaptation. Studies in simulated lesion models and in patients with stroke reveal both local and distal functional disturbances, using measures of anatomical brain networks. Investigations over the last two decades have sought to use diffusion tensor imaging tractography data in the context of intracranial tumours to improve surgical planning, intraoperative functional localization, and post-operative interpretation of functional change. In this study, we used diffusion tensor imaging tractography to assess the impact of tumour location on the white matter structural network. To better understand how various lobe localized gliomas impact the topology underlying efficiency of information transfer between brain regions, we identified the major alterations in brain network connectivity patterns between the ipsilesional versus contralesional hemispheres in patients with gliomas localized to the frontal, parietal or temporal lobe. Results were indicative of altered network efficiency and the role of specific brain regions unique to different lobe localized gliomas. This work draws attention to connections and brain regions which have shared structural susceptibility in frontal, parietal and temporal lobe glioma cases. This study also provides a preliminary anatomical basis for understanding which affected white matter pathways may contribute to preoperative patient symptomology.

Tumor ablation due to inhomogeneous anisotropic diffusion in generic three-dimensional topologies

In recent decades computer-aided technologies have become prevalent in medicine, however, cancer drugs are often only tested on in vitro cell lines from biopsies. We derive a full three-dimensional model of inhomogeneous -anisotropic diffusion in a tumor region coupled to a binary population model, which simulates in vivo scenarios faster than traditional cell-line tests. The diffusion tensors are acquired using diffusion tensor magnetic resonance imaging from a patient diagnosed with glioblastoma multiform. Then we numerically simulate the full model with finite element methods and produce drug concentration heat maps, apoptosis hotspots, and dose-response curves. Finally, predictions are made about optimal injection locations and volumes, which are presented in a form that can be employed by doctors and oncologists.

Decompressive craniectomy of post-traumatic brain injury: an in silico modelling approach for intracranial hypertension management

Traumatic brain injury (TBI) causes brain edema that induces increased intracranial pressure and decreased cerebral perfusion. Decompressive craniectomy has been recommended as a surgical procedure for the management of swollen brain and intracranial hypertension. Proper location and size of a decompressive craniectomy, however, remain controversial and no clinical guidelines are available. Mathematical and computational (in silico) models can predict the optimum geometric conditions and provide insights for the brain mechanical response following a decompressive craniectomy. In this work, we present a finite element model of post-traumatic brain injury and decompressive craniectomy that incorporates a biphasic, nonlinear biomechanical model of the brain. A homogenous pressure is applied in the brain to represent the intracranial pressure loading caused by the tissue swelling and the models calculate the deformations and stresses in the brain as well as the herniated volume of the brain tissue that exits the skull following craniectomy. Simulations for different craniectomy geometries (unilateral, bifrontal and bifrontal with midline bar) and sizes are employed to identify optimal clinical conditions of decompressive craniectomy. The reported results for the herniated volume of the brain tissue as a function of the intracranial pressure loading under a specific geometry and size of craniectomy are exceptionally relevant for decompressive craniectomy planning.

Glioma Growth Prediction via Generative Adversarial Learning from Multi-Time Points Magnetic Resonance Images

Gliomas are the most dominant and lethal type of brain tumors. Growth prediction is significant to quantify tumor aggressiveness, improve therapy planning, and estimate patients' survival time. This is commonly addressed in literature using mathematical models guided by multi-time point scans of multi/single-modal data for the same subject. However, these models are mechanism-based and heavily rely on complicated mathematical formulations of partial differential equations with few parameters that are insufficient to capture different patterns and other characteristics of gliomas. In this paper, we propose a 3D generative adversarial networks (GANs) for glioma growth prediction. Specifically, we stack 2 GANs with conditional initialization of segmented feature maps. Furthermore, we employ Dice loss in our objective function and devised 3D U-Net architecture for better image generation. The proposed method is trained and validated using 3D patch-based strategy on real magnetic resonance images of 9 subjects with 3 time points. Experimental results show that the proposed method can be successfully used for glioma growth prediction with satisfactory performance.

From tumour perfusion to drug delivery and clinical translation of in silico cancer models

In silico cancer models have demonstrated great potential as a tool to improve drug design, optimise the delivery of drugs to target sites in the host tissue and, hence, improve therapeutic efficacy and patient outcome. However, there are significant barriers to the successful translation of in silico technology from bench to bedside. More precisely, the specification of unknown model parameters, the necessity for models to adequately reflect in vivo conditions, and the limited amount of pertinent validation data to evaluate models' accuracy and assess their reliability, pose major obstacles in the path towards their clinical translation. This review aims to capture the state-of-the-art in in silico cancer modelling of vascularised solid tumour growth, and identify the important advances and barriers to success of these models in clinical oncology. Particular emphasis has been put on continuum-based models of cancer since they - amongst the class of mechanistic spatio-temporal modelling approaches - are well-established in simulating transport phenomena and the biomechanics of tissues, and have demonstrated potential for clinical translation. Three important avenues in in silico modelling are considered in this contribution: first, since systemic therapy is a major cancer treatment approach, we start with an overview of the tumour perfusion and angiogenesis in silico models. Next, we present the state-of-the-art in silico work encompassing the delivery of chemotherapeutic agents to cancer nanomedicines through the bloodstream, and then review continuum-based modelling approaches that demonstrate great promise for successful clinical translation. We conclude with a discussion of what we view to be the key challenges and opportunities for in silico modelling in personalised and precision medicine.

On PTV definition for glioblastoma based on fiber tracking of diffusion tensor imaging data

Radiotherapy (RT) is commonly applied for the treatment of glioblastoma multiforme (GBM). Following the planning target volume (PTV) definition procedure standardized in guidelines, a 20% risk of missing non-local recurrences is present. Purpose of this study was to evaluate whether diffusion tensor imaging (DTI)-based fiber tracking may be beneficial for PTV definition taking into account the prediction of distant recurrences. 56 GBM patients were examined with magnetic resonance imaging (MRI) including DTI performed before RT after resection of the primary tumor. Follow-up MRIs were acquired in three month intervals. For the seven patients with a distant recurrence, fiber tracking was performed with three algorithms and it was evaluated whether connections existed from the primary tumor region to the distant recurrence. It depended strongly on the used tracking algorithm and the used tracking parameters whether a connection was observed. Most of the connections were weak and thus not usable for PTV definition. Only in one of the seven patients with a recurring tumor, a clear connection was present. It seems unlikely that DTI-based fiber tracking can be beneficial for predicting distant recurrences in the planning of PTVs for glioblastoma multiforme.

Fiber-tract localized diffusion coefficients highlight patterns of white matter disruption induced by proximity to glioma

Gliomas account for 26.5% of all primary central nervous system tumors. Recent studies have used diffusion tensor imaging (DTI) to extract white matter fibers and the diffusion coefficients derived from MR processing to provide useful, non-invasive insights into the extent of tumor invasion, axonal integrity, and gross differentiation of glioma from metastasis. Here, we extend this work by examining whether a tract-based analysis can improve non-invasive localization of tumor impact on white matter integrity. This study retrospectively analyzed preoperative magnetic resonance sequences highlighting contrast enhancement and DTI scans of 13 subjects that were biopsy-confirmed to have either high or low-grade glioma. We reconstructed the corticospinal tract and superior longitudinal fasciculus by applying atlas-based regions of interest to fibers derived from whole-brain deterministic streamline tractography. Within-subject comparison of hemispheric diffusion coefficients (e.g., fractional anisotropy and mean diffusivity) indicated higher levels of white matter degradation in the ipsilesional hemisphere. Novel application of along-tract analyses revealed that tracts traversing the tumor region showed significant white matter degradation compared to the contralesional hemisphere and ipsilesional tracts displaced by the tumor.