skip to main content

Explore Research Products in the NSF-PAR

Title: Mathematical modelling of the dynamics of image-informed tumor habitats in a murine model of glioma
Abstract

Tumors exhibit high molecular, phenotypic, and physiological heterogeneity. In this effort, we employ quantitative magnetic resonance imaging (MRI) data to capture this heterogeneity through imaging-based subregions or “habitats” in a murine model of glioma. We then demonstrate the ability to model and predict the growth of the habitats using coupled ordinary differential equations (ODEs) in the presence and absence of radiotherapy. Female Wistar rats (N = 21) were inoculated intracranially with 106C6 glioma cells, a subset of which received 20 Gy (N = 5) or 40 Gy (N = 8) of radiation. All rats underwent diffusion-weighted and dynamic contrast-enhanced MRI at up to seven time points. All MRI data at each visit were subsequently clustered usingk-means to identify physiological tumor habitats. A family of four models consisting of three coupled ODEs were developed and calibrated to the habitat time series of control and treated rats and evaluated for predictive capability. The Akaike Information Criterion was used for model selection, and the normalized sum-of-square-error (SSE) was used to evaluate goodness-of-fit in model calibration and prediction. Three tumor habitats with significantly different imaging data characteristics (p < 0.05) were identified: high-vascularity high-cellularity, low-vascularity high-cellularity, and low-vascularity low-cellularity. Model selection resulted in a five-parameter model whose predictions of habitat dynamics yielded SSEs that were similar to the SSEs from the calibrated model. It is thus feasible to mathematically describe habitat dynamics in a preclinical model of glioma using biology-based ODEs, showing promise for forecasting heterogeneous tumor behavior.

  more » « less
NSF-PAR ID:
10397940
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Scientific Reports
Volume:
13
Issue:
1
ISSN:
2045-2322
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ; ; ; ; ; et al ( , Scientific Reports)
    null (Ed.)
    Abstract Radiogenomics uses machine-learning (ML) to directly connect the morphologic and physiological appearance of tumors on clinical imaging with underlying genomic features. Despite extensive growth in the area of radiogenomics across many cancers, and its potential role in advancing clinical decision making, no published studies have directly addressed uncertainty in these model predictions. We developed a radiogenomics ML model to quantify uncertainty using transductive Gaussian Processes (GP) and a unique dataset of 95 image-localized biopsies with spatially matched MRI from 25 untreated Glioblastoma (GBM) patients. The model generated predictions for regional EGFR amplification status (a common and important target in GBM) to resolve the intratumoral genetic heterogeneity across each individual tumor—a key factor for future personalized therapeutic paradigms. The model used probability distributions for each sample prediction to quantify uncertainty, and used transductive learning to reduce the overall uncertainty. We compared predictive accuracy and uncertainty of the transductive learning GP model against a standard GP model using leave-one-patient-out cross validation. Additionally, we used a separate dataset containing 24 image-localized biopsies from 7 high-grade glioma patients to validate the model. Predictive uncertainty informed the likelihood of achieving an accurate sample prediction. When stratifying predictions based on uncertainty, we observed substantially higher performance in the group cohort (75% accuracy, n = 95) and amongst sample predictions with the lowest uncertainty (83% accuracy, n = 72) compared to predictions with higher uncertainty (48% accuracy, n = 23), due largely to data interpolation (rather than extrapolation). On the separate validation set, our model achieved 78% accuracy amongst the sample predictions with lowest uncertainty. We present a novel approach to quantify radiogenomics uncertainty to enhance model performance and clinical interpretability. This should help integrate more reliable radiogenomics models for improved medical decision-making. 
    more » « less
  2. ; ; ; ; ; ; ( , Magnetic Resonance in Medicine)
    Purposes

    To develop and evaluate a boundary informed electrical properties tomography (BIEPT) technique for high‐resolution imaging of tumor electrical properties (EPs) heterogeneity on a rodent tumor xenograft model.

     Methods

    Tumor EP distributions were inferred from a reference area external to the tumor, as well as internal EP spatial variations derived from a plurality of relative transmit B1measurements at 7T. Edge sparsity constraint was enforced to enhance numerical stability. Phantom experiments were performed to determine the imaging accuracy and sensitivity for structures of various EP values, as well as geometrical sizes down to 1.5 mm. Numerical simulation of a realistic rodent model was used to quantify the algorithm performance in the presence of noise. Eleven athymic rats with human breast cancer xenograft were imaged in vivo, and representative pathological samples were acquired for comparison.

     Results

    Reconstructed EPs of the phantoms correspond well to the ground truth acquired from dielectric probe measurements, with the smallest structure reliably detectable being 3 mm. EPs heterogeneity inside a tumor is successfully retrieved in both simulated and experimental cases. In vivo tumor imaging results demonstrate similar local features and spatial patterns to anatomical MRI and pathological slides. The imaged conductivity of necrotic tissue is higher than that of viable tissues, which agrees with our expectation.

     Conclusion

    BIEPT enables robust detection of tumor EPs heterogeneity with high accuracy and sensitivity to small structures. The retrieved quantitative EPs reflect tumor pathological features (e.g., necrosis). These results provide strong rationale to further expand BIEPT studies toward pathological conditions where EPs may yield valuable, non‐invasive biomarkers.

     
    more » « less
  3. ; ; ; ; ( , Magnetic Resonance in Medicine)
    Purpose

    Diffusion MRI offers insight into ischemic stroke progression in both human and rodent models. However, diffusion MRI to evaluate therapeutic application of mesenchymal stem cells is limited. Robust analytical techniques are required to identify potential physiological changes as a function of cell therapy in stroke. Here, we seek to establish Neurite Orientation Dispersion and Density Imaging (NODDI) as a feasible method in evaluating stroke evolution in response to cell‐based therapeutics.

     Methods

    Diffusion MRI data at 21.1T were acquired from 16 male rats. Rats were grouped randomly: naïve (baseline, N = 5), stroke with injections of phosphate buffered saline (N = 6), stroke with injection of 2D human mesenchymal stem cells (hMSC, N = 5). Data were acquired on days 1, 3, 7, and 21 post‐surgery. DTI and NODDI maps were generated, with regions of interest placed in the ischemic hemisphere external capsule and striatum. Diffusion parameters were compared between groups each day, and within groups across hemispheres and longitudinally. Behavioral characterizations were on days 0 (pre‐surgery), 3, 7, 14, and 21.

     Results

    The 2D hMSC preserved diffusional restriction in the external capsule compared to saline (day 1: MD,P= .4060; AD,P= .0220). NODDI indicates that hMSC may have preserved intracellular volume fractions (ICVF: day 1,P= .0086; day 3,P= .0021; day 21,P= .0383). Diffusion metrics of hMSC treated animals were comparable to naïve for the external capsule.

     Conclusions

    NODDI compliments DTI metrics, enhances interpretation of tissue outcome in ischemic stroke following hMSC application, and may be useful in evaluating or predicting therapeutic response.

     
    more » « less
  4. ; ; ; ; ; ( , Ecology of Freshwater Fish)
    Abstract

    Landscape‐level geomorphic processes influence the spatial and temporal arrangement of fish habitats in freshwater ecosystems and fishes move across riverscapes, selecting a suite of habitats to maximise fitness. Here, we explore the influence of geomorphology on stream channel attributes and assess Broad Whitefish (Coregonus nasus) spawning habitat potential in the Colville River in Arctic Alaska. Using high‐resolution digital surface models (5 m2), we quantified the stream network extent and summarised channel habitat attributes continuously across the drainage network. Next, we developed an intrinsic potential (IP) model for Broad Whitefish by using geomorphic channel parameters previously understood to be associated with spawning habitats (channel width, median substrate size and channel braiding) to estimate the potential of streams across the Colville River watershed to provide spawning habitat. Our model results show the majority of habitat with high IP (≥0.6) was located within the braided sections of the main channel, which encompass >1548 km, but only 2% of the total channel network. The IP model was tested by tracking radio‐tagged Broad Whitefish using aerial surveys. Prespawn fish moved into the watershed starting mid‐July and mostly used habitat with moderate to very high IP in the middle and lower watershed. Several individuals were relocated in smaller multichannels with vegetated bars that contained very low IP (≤0.2), suggesting that other factors, such as hyporheic flow, may also influence spawning habitat selection. Our study demonstrates that IP modelling offers a useful method to quantify spawning habitat potential in data‐poor riverscapes, providing useful information for managers to assess potential anthropogenic impacts and develop conservation plans to protect essential Broad Whitefish habitat.

     
    more » « less
  5. ; ; ; ; ; ; ; ; ; ; et al ( , Ecological Monographs)
    Abstract

    Behavioral thermoregulation is an efficient mechanism to buffer the physiological effects of climate change. Thermal ecology studies have traditionally tested how thermal constraints shape thermoregulatory behaviors without accounting for the potential major effects of landscape structure and water availability. Thus, we lack a general understanding of the multifactorial determinants of thermoregulatory behaviors in natural populations. In this study, we quantified the relative contribution of elevation, thermal gradient, moisture gradient, and landscape structure in explaining geographic variation in thermoregulation strategies of a terrestrial ectotherm species. We measured field‐active body temperature, thermal preferences, and operative environmental temperatures to calculate thermoregulation indices, including thermal quality of the habitat and thermoregulation efficiency for a very large sample of common lizards (Zootoca vivipara) from 21 populations over 3 yr across the Massif Central mountain range in France. We used an information‐theoretic approach to compare eight a priori thermo‐hydroregulation hypotheses predicting how behavioral thermoregulation should respond to environmental conditions. Environmental characteristics exerted little influence on thermal preference with the exception that females from habitats with permanent access to water had lower thermal preferences. Field body temperatures and accuracy of thermoregulation were best predicted by the interaction between air temperature and a moisture index. In mesic environments, field body temperature and thermoregulation accuracy increased with air temperature, but they decreased in drier habitats. Thermoregulation efficiency (difference between thermoregulation inaccuracy and the thermal quality of the habitat) was maximized in cooler and more humid environments and was mostly influenced by the thermal quality of the habitat. Our study highlights complex patterns of variation in thermoregulation strategies, which are mostly explained by the interaction between temperature and water availability, independent of the elevation gradient or thermal heterogeneity. Although changes in landscape structure were expected to be the main driver of extinction rate of temperate zone ectotherms with ongoing global change, we conclude that changes in water availability coupled with rising temperatures might have a drastic impact on the population dynamics of some ectotherm species.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email