ASU Scholarship Showcase

Five immunocompetent C57BL/6-cBrd/cBrd/Cr (albino C57BL/6) mice were injected with GL261-luc2 cells, a cell line sharing characteristics of human glioblastoma multiforme (GBM). The mice were imaged using magnetic resonance (MR) at five separate time points to characterize growth and development of the tumor. After 25 days, the final tumor volumes of the mice varied from 12 mm³ to 62 mm³, even though mice were inoculated from the same tumor cell line under carefully controlled conditions. We generated hypotheses to explore large variances in final tumor size and tested them with our simple reaction-diffusion model in both a 3-dimensional (3D) finite difference method and a 2-dimensional (2D) level set method. The parameters obtained from a best-fit procedure, designed to yield simulated tumors as close as possible to the observed ones, vary by an order of magnitude between the three mice analyzed in detail. These differences may reflect morphological and biological variability in tumor growth, as well as errors in the mathematical model, perhaps from an oversimplification of the tumor dynamics or nonidentifiability of parameters. Our results generate parameters that match other experimental in vitro and in vivo measurements. Additionally, we calculate wave speed, which matches with other rat and human measurements.

Background:
Data assimilation refers to methods for updating the state vector (initial condition) of a complex spatiotemporal model (such as a numerical weather model) by combining new observations with one or more prior forecasts. We consider the potential feasibility of this approach for making short-term (60-day) forecasts of the growth and spread of a malignant brain cancer (glioblastoma multiforme) in individual patient cases, where the observations are synthetic magnetic resonance images of a hypothetical tumor.

Results:
We apply a modern state estimation algorithm (the Local Ensemble Transform Kalman Filter), previously developed for numerical weather prediction, to two different mathematical models of glioblastoma, taking into account likely errors in model parameters and measurement uncertainties in magnetic resonance imaging. The filter can accurately shadow the growth of a representative synthetic tumor for 360 days (six 60-day forecast/update cycles) in the presence of a moderate degree of systematic model error and measurement noise.

Conclusions:
The mathematical methodology described here may prove useful for other modeling efforts in biology and oncology. An accurate forecast system for glioblastoma may prove useful in clinical settings for treatment planning and patient counseling.

The apolipoprotein E (APOE) e4 genotype is a powerful risk factor for late-onset Alzheimer’s disease (AD). In the Alzheimer’s Disease Neuroimaging Initiative (ADNI) cohort, we previously reported significant baseline structural differences in APOE e4 carriers relative to non-carriers, involving the left hippocampus more than the right—a difference more pronounced in e4 homozygotes than heterozygotes. We now examine the longitudinal effects of APOE genotype on hippocampal morphometry at 6-, 12- and 24-months, in the ADNI cohort. We employed a new automated surface registration system based on conformal geometry and tensor-based morphometry. Among different hippocampal surfaces, we computed high-order correspondences, using a novel inverse-consistent surface-based fluid registration method and multivariate statistics consisting of multivariate tensor-based morphometry (mTBM) and radial distance. At each time point, using Hotelling’s T² test, we found significant morphological deformation in APOE e4 carriers relative to non-carriers in the full cohort as well as in the non-demented (pooled MCI and control) subjects at each follow-up interval. In the complete ADNI cohort, we found greater atrophy of the left hippocampus than the right, and this asymmetry was more pronounced in e4 homozygotes than heterozygotes. These findings, combined with our earlier investigations, demonstrate an e4 dose effect on accelerated hippocampal atrophy, and support the enrichment of prevention trial cohorts with e4 carriers.

Many children born preterm exhibit frontal executive dysfunction, behavioral problems including attentional deficit/hyperactivity disorder and attention related learning disabilities. Anomalies in regional specificity of cortico-striato-thalamo-cortical circuits may underlie deficits in these disorders. Nonspecific volumetric deficits of striatal structures have been documented in these subjects, but little is known about surface deformation in these structures. For the first time, here we found regional surface morphological differences in the preterm neonatal ventral striatum. We performed regional group comparisons of the surface anatomy of the striatum (putamen and globus pallidus) between 17 preterm and 19 term-born neonates at term-equivalent age. We reconstructed striatal surfaces from manually segmented brain magnetic resonance images and analyzed them using our in-house conformal mapping program. All surfaces were registered to a template with a new surface fluid registration method. Vertex-based statistical comparisons between the two groups were performed via four methods: univariate and multivariate tensor-based morphometry, the commonly used medial axis distance, and a combination of the last two statistics. We found statistically significant differences in regional morphology between the two groups that are consistent across statistics, but more extensive for multivariate measures. Differences were localized to the ventral aspect of the striatum. In particular, we found abnormalities in the preterm anterior/inferior putamen, which is interconnected with the medial orbital/prefrontal cortex and the midline thalamic nuclei including the medial dorsal nucleus and pulvinar. These findings support the hypothesis that the ventral striatum is vulnerable, within the cortico-stiato-thalamo-cortical neural circuitry, which may underlie the risk for long-term development of frontal executive dysfunction, attention deficit hyperactivity disorder and attention-related learning disabilities in preterm neonates.

Land surface energy balance in a built environment is widely modelled using urban canopy models with representation of building arrays as big street canyons. Modification of this simplified geometric representation, however, leads to challenging numerical difficulties in improving physical parameterization schemes that are deterministic in nature. In this paper, we develop a stochastic algorithm to estimate view factors between canyon facets in the presence of shade trees based on Monte Carlo simulation, where an analytical formulation is inhibited by the complex geometry. The model is validated against analytical solutions of benchmark radiative problems as well as field measurements in real street canyons. In conjunction with the matrix method resolving infinite number of reflections, the proposed model is capable of predicting the radiative exchange inside the street canyon with good accuracy. Modeling of transient evolution of thermal filed inside the street canyon using the proposed method demonstrate the potential of shade trees in mitigating canyon surface temperatures as well as saving of building energy use. This new numerical framework also deepens our insight into the fundamental physics of radiative heat transfer and surface energy balance for urban climate modeling.

Studies on urban heat island (UHI) have been more than a century after the phenomenon was first discovered in the early 1800s. UHI emerges as the source of many urban environmental problems and exacerbates the living environment in cities. Under the challenges of increasing urbanization and future climate changes, there is a pressing need for sustainable adaptation/mitigation strategies for UHI effects, one popular option being the use of reflective materials. While it is introduced as an effective method to reduce temperature and energy consumption in cities, its impacts on environmental sustainability and large-scale non-local effect are inadequately explored. This paper provides a synthetic overview of potential environmental impacts of reflective materials at a variety of scales, ranging from energy load on a single building to regional hydroclimate. The review shows that mitigation potential of reflective materials depends on a set of factors, including building characteristics, urban environment, meteorological and geographical conditions, to name a few. Precaution needs to be exercised by city planners and policy makers for large-scale deployment of reflective materials before their environmental impacts, especially on regional hydroclimates, are better understood. In general, it is recommended that optimal strategy for UHI needs to be determined on a city-by-city basis, rather than adopting a “one-solution-fits-all” strategy.

In this paper, we develop a new automated surface registration system based on surface conformal parameterization by holomorphic 1-forms, inverse consistent surface fluid registration, and multivariate tensor-based morphometty (mTBM). First, we conformally map a surface onto a planar rectangle space with holomorphic 1-forms. Second, we compute surface conformal representation by combining its local conformal factor and mean curvature and linearly scale the dynamic range of the conformal representation to form the feature image of the surface. Third, we align the feature image with a chosen template image via the fluid image registration algorithm, which has been extended into the curvilinear coordinates to adjust for the distortion introduced by surface parameterization. The inverse consistent image registration algorithm is also incorporated in the system to jointly estimate the forward and inverse transformations between the study and template images. This alignment induces a corresponding deformation on the surface. We tested the system on Alzheimer's Disease Neuroimaging Initiative (ADNI) baseline dataset to study AD symptoms on hippocampus. In our system, by modeling a hippocampus as a 3D parametric surface, we nonlinearly registered each surface with a selected template surface. Then we used mTBM to analyze the morphometry difference between diagnostic groups. Experimental results show that the new system has better performance than two publicly available subcortical surface registration tools: FIRST and SPHARM. We also analyzed the genetic influence of the Apolipoprotein E(is an element of)4 allele (ApoE4), which is considered as the most prevalent risk factor for AD. Our work successfully detected statistically significant difference between ApoE4 carriers and non-carriers in both patients of mild cognitive impairment (MCI) and healthy control subjects. The results show evidence that the ApoE genotype may be associated with accelerated brain atrophy so that our work provides a new MRI analysis tool that may help presymptomatic AD research.

The hysteresis effect in diurnal cycles of net radiation R-n and ground heat flux G(0) has been observed in many studies, while the governing mechanism remains vague. In this study, we link the phenomenology of hysteresis loops to the wave phase difference between the diurnal evolutions of various terms in the surface energy balance. R-n and G(0) are parameterized with the incoming solar radiation and the surface temperature as two control parameters of the surface energy partitioning. The theoretical analysis shows that the vertical water flux W and the scaled ratio A(s)*/A(T)* (net shortwave radiation to outgoing longwave radiation) play crucial roles in shaping hysteresis loops of R-n and G(0). Comparisons to field measurements indicate that hysteresis loops for different land covers can be well captured by the theoretical model, which is also consistent with Camuffo-Bernadi formula. This study provides insight into the surface partitioning and temporal evolution of the energy budget at the land surface.

This study examines the spatial and temporal patterns of the surface urban heat island (SUHI) intensity in the Phoenix metropolitan area and the relationship with land use land cover (LULC) change between 2000 and 2014. The objective is to identify specific regions in Phoenix that have been increasingly heated and cooled to further understand how LULC change influences the SUHI intensity. The data employed include MODerate-resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST) 8-day composite June imagery, and classified LULC maps generated using 2000 and 2014 Landsat imagery. Results show that the regions that experienced the most significant LST changes during the study period are primarily on the outskirts of the Phoenix metropolitan area for both daytime and nighttime. The conversion to urban, residential, and impervious surfaces from all other LULC types has been identified as the primary cause of the UHI effect in Phoenix. Vegetation cover has been shown to significantly lower LST for both daytime and nighttime due to its strong cooling effect by producing more latent heat flux and less sensible heat flux. We suggest that urban planners, decision-makers, and city managers formulate new policies and regulations that encourage residential, commercial, and industrial developers to include more vegetation when planning new construction.

Predicting the timing of a castrate resistant prostate cancer is critical to lowering medical costs and improving the quality of life of advanced prostate cancer patients. We formulate, compare and analyze two mathematical models that aim to forecast future levels of prostate-specific antigen (PSA). We accomplish these tasks by employing clinical data of locally advanced prostate cancer patients undergoing androgen deprivation therapy (ADT). While these models are simplifications of a previously published model, they fit data with similar accuracy and improve forecasting results. Both models describe the progression of androgen resistance. Although Model 1 is simpler than the more realistic Model 2, it can fit clinical data to a greater precision. However, we found that Model 2 can forecast future PSA levels more accurately. These findings suggest that including more realistic mechanisms of androgen dynamics in a two population model may help androgen resistance timing prediction.