Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userWhere did the tumor start? An inverse solver with sparse localization for tumor growth models
392
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a numerical scheme for solving an inverse problem for parameter estimation in tumor growth models for glioblastomas, a form of aggressive primary brain tumor. The growth model is a reaction-diffusion partial differential equation (PDE) for the tumor concentration. We use a PDE-constrained optimization formulation for the inverse problem. The unknown parameters are the reaction coefficient (proliferation), the diffusion coefficient (infiltration), and the initial condition field for the tumor PDE. Segmentation of Magnetic Resonance Imaging (MRI) scans from a single time snapshot drive the inverse problem where segmented tumor regions serve as partial observations of the tumor concentration. The precise time relative to tumor initiation is unknown, which poses an additional difficulty for inversion. We perform a frozen-coefficient spectral analysis and show that the inverse problem is severely ill-posed. We introduce a biophysically motivated regularization on the tumor initial condition. In particular, we assume that the tumor starts at a few locations (enforced with a sparsity constraint) and that the initial condition magnitude in the maximum norm equals one. We solve the resulting optimization problem using an inexact quasi-Newton method combined with a compressive sampling algorithm for the sparsity constraint. Our implementation uses PETSc and AccFFT libraries. We conduct numerical experiments on synthetic and clinical images to highlight the improved performance of our solver over an existing solver that uses a two-norm regularization for the calibration parameters. The existing solver is unable to localize the initial condition. Our new solver can localize the initial condition and recover infiltration and proliferation. In clinical datasets (for which the ground truth is unknown), our solver results in qualitatively different solutions compared to the existing solver.
rate research
Read More
In this article, we present a multispecies reaction-advection-diffusion partial differential equation (PDE) coupled with linear elasticity for modeling tumor growth. The model aims to capture the phenomenological features of glioblastoma multiforme observed in magnetic resonance imaging (MRI) scans. These include enhancing and necrotic tumor structures, brain edema and the so called mass effect, that is, the deformation of brain tissue due to the presence of the tumor. The multispecies model accounts for proliferating, invasive and necrotic tumor cells as well as a simple model for nutrition consumption and tumor-induced brain edema. The coupling of the model with linear elasticity equations with variable coefficients allows us to capture the mechanical deformations due to the tumor growth on surrounding tissues. We present the overall formulation along with a novel operator-splitting scheme with components that include linearly-implicit preconditioned elliptic solvers, and semi-Lagrangian method for advection. Also, we present results showing simulated MRI images which highlight the capability of our method to capture the overall structure of glioblastomas in MRIs.
We propose a method for extracting physics-based biomarkers from a single multiparametric Magnetic Resonance Imaging (mpMRI) scan bearing a glioma tumor. We account for mass effect, the deformation of brain parenchyma due to the growing tumor, which on its own is an important radiographic feature but its automatic quantification remains an open problem. In particular, we calibrate a partial differential equation (PDE) tumor growth model that captures mass effect, parameterized by a single scalar parameter, tumor proliferation, migration, while localizing the tumor initiation site. The single-scan calibration problem is severely ill-posed because the precancerous, healthy, brain anatomy is unknown. To address the ill-posedness, we introduce an ensemble inversion scheme that uses a number of normal subject brain templates as proxies for the healthy precancer subject anatomy. We verify our solver on a synthetic dataset and perform a retrospective analysis on a clinical dataset of 216 glioblastoma (GBM) patients. We analyze the reconstructions using our calibrated biophysical model and demonstrate that our solver provides both global and local quantitative measures of tumor biophysics and mass effect. We further highlight the improved performance in model calibration through the inclusion of mass effect in tumor growth models -- including mass effect in the model leads to 10% increase in average dice coefficients for patients with significant mass effect. We further evaluate our model by introducing novel biophysics-based features and using them for survival analysis. Our preliminary analysis suggests that including such features can improve patient stratification and survival prediction.
The comprehension of tumor growth is a intriguing subject for scientists. New researches has been constantly required to better understand the complexity of this phenomenon. In this paper, we pursue a physical description that account for some experimental facts involving avascular tumor growth. We have proposed an explanation of some phenomenological (macroscopic) aspects of tumor, as the spatial form and the way it growths, from a individual-level (microscopic) formulation. The model proposed here is based on a simple principle: competitive interaction between the cells dependent on their mutual distances. As a result, we reproduce many empirical evidences observed in real tumors, as exponential growth in their early stages followed by a power law growth. The model also reproduces the fractal space distribution of tumor cells and the universal behavior presented in animals and tumor growth, conform reported by West, Guiot {it et. al.}cite{West2001,Guiot2003}. The results suggest that the universal similarity between tumor and animal growth comes from the fact that both are described by the same growth equation - the Bertalanffy-Richards model - even they does not necessarily share the same biophysical properties.
We derive a full 3-dimensional (3-D) model of inhomogeneous -- anisotropic diffusion in a tumor region coupled to a binary population model. The diffusion tensors are acquired using Diffusion Tensor Magnetic Resonance Imaging (DTI) from a patient diagnosed with glioblastoma multiform (GBM). Then we numerically simulate the full model with Finite Element Method (FEM) and produce drug concentration heat maps, apoptosis regions, 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.
Presently 4T-1 luc cells were irradiated with proton under ultra-high dose rate FLASH or with gamma-ray with conventional dose rate, and then subcutaneous vaccination with or without Mn immuno-enhancing adjuvant into the mice for three times. One week later, we injected untreated 4T-1 luc cells on the other side of the vaccinated mice, and found that the untreated 4T-1 luc cells injected later nearly totally did not grow tumor (1/17) while controls without previous vaccination all grow tumors (18/18). The result is very interesting and the findings may help to explore in situ tumor vaccination as well as new combined radiotherapy strategies to effectively ablate primary and disseminated tumors. To our limited knowledge, this is the first paper reporting the high efficiency induction of systemic vaccination suppressing the metastasized/disseminated tumor progression.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
