Tumor invasion margin on the Riemannian space of brain fibers

Glioma is one of the most challenging types of brain tumors to treat or control locally. One of the main problems is to determine which areas of the apparently normal brain contain glioma cells, as gliomas are known to infiltrate several centimeters beyond the clinically apparent lesion that is visualized on standard Computed Tomography scans (CT) or Magnetic Resonance Images (MRIs). To ensure that radiation treatment encompasses the whole tumor, including the cancerous cells not revealed by MRI, doctors treat the volume of brain that extends 2cm out from the margin of the visible tumor. This approach does not consider varying tumor-growth dynamics in different brain tissues, thus it may result in killing some healthy cells while leaving cancerous cells alive in the other areas. These cells may cause recurrence of the tumor later in time, which limits the effectiveness of the therapy. Knowing that glioma cells preferentially spread along nerve fibers, we propose the use of a geodesic distance on the Riemannian manifold of brain diffusion tensors to replace the Euclidean distance used in the clinical practice and to correctly identify the tumor invasion margin. This mathematical model results in a first-order Partial Differential Equation (PDE) that can be numerically solved in a stable and consistent way. To compute the geodesic distance, we use actual Diffusion Weighted Imaging (DWI) data from 11 patients with glioma and compare our predicted infiltration distance map with actual grwoth in follow-up MRI scans. Results show improvement in predicting the invasion margin when using the geodesic distance as opposed to the 2cm conventional Euclidean distance.