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.

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.

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.

The electronic band structure of MoS₂, MoSe₂, WS₂, and WSe₂, crystals has been studied at various hydrostatic pressures experimentally by photoreflectance (PR) spectroscopy and theoretically within the density functional theory (DFT). In the PR spectra direct optical transitions (A and B) have been clearly observed and pressure coefficients have been determined for these transitions to be: α_A = 2.0 ± 0.1 and α_B = 3.6 ± 0.1 meV/kbar for MoS₂, α_A = 2.3 ± 0.1 and α_B = 4.0 ± 0.1 meV/kbar for MoSe₂, α_A = 2.6 ± 0.1 and α_B = 4.1 ± 0.1 meV/kbar for WS₂, α_A = 3.4 ± 0.1 and α_B = 5.0 ± 0.5 meV/kbar for WSe₂. It has been found that these coefficients are in an excellent agreement with theoretical predictions. In addition, a comparative study of different computational DFT approaches has been performed and analyzed. For indirect gap the pressure coefficient have been determined theoretically to be −7.9, −5.51, −6.11, and −3.79, meV/kbar for MoS₂, MoSe₂, WS₂, and WSe₂, respectively. The negative values of this coefficients imply a narrowing of the fundamental band gap with the increase in hydrostatic pressure and a semiconductor to metal transition for MoS₂, MoSe₂, WS₂, and WSe₂, crystals at around 140, 180, 190, and 240 kbar, respectively.

The ongoing Zika virus (ZIKV) epidemic in the Americas poses a major global public health emergency. While ZIKV is transmitted from human to human by bites of Aedes mosquitoes, recent evidence indicates that ZIKV can also be transmitted via sexual contact with cases of sexually transmitted ZIKV reported in Argentina, Canada, Chile, France, Italy, New Zealand, Peru, Portugal, and the USA. Yet, the role of sexual transmission on the spread and control of ZIKV infection is not well-understood. We introduce a mathematical model to investigate the impact of mosquito-borne and sexual transmission on the spread and control of ZIKV and calibrate the model to ZIKV epidemic data from Brazil, Colombia, and El Salvador. Parameter estimates yielded a basic reproduction number R_0 = 2.055 (95% CI: 0.523–6.300), in which the percentage contribution of sexual transmission is 3.044% (95% CI: 0.123–45.73). Our sensitivity analyses indicate that R_0 is most sensitive to the biting rate and mortality rate of mosquitoes while sexual transmission increases the risk of infection and epidemic size and prolongs the outbreak. Prevention and control efforts against ZIKV should target both the mosquito-borne and sexual transmission routes.

Binary transition metal dichalcogenide monolayers share common properties such as a direct optical bandgap, spin-orbit splittings of hundreds of meV, light–matter interaction dominated by robust excitons and coupled spin-valley states. Here we demonstrate spin-orbit-engineering in Mo[_(1-x)]W_xSe₂ alloy monolayers for optoelectronics and applications based on spin- and valley-control. We probe the impact of the tuning of the conduction band spin-orbit spin-splitting on the bright versus dark exciton population. For MoSe₂ monolayers, the photoluminescence intensity decreases as a function of temperature by an order of magnitude (4–300 K), whereas for WSe₂ we measure surprisingly an order of magnitude increase. The ternary material shows a trend between these two extreme behaviors. We also show a non-linear increase of the valley polarization as a function of tungsten concentration, where 40% tungsten incorporation is sufficient to achieve valley polarization as high as in binary WSe₂.

We present two-dimensional Mg(OH)₂ sheets and their vertical heterojunctions with CVD-MoS₂ for the first time as flexible 2D insulators with anomalous lattice vibration and chemical and physical properties. New hydrothermal crystal growth technique enabled isolation of environmentally stable monolayer Mg(OH)₂ sheets. Raman spectroscopy and vibrational calculations reveal that the lattice vibrations of Mg(OH)₂ have fundamentally different signature peaks and dimensionality effects compared to other 2D material systems known to date. Sub-wavelength electron energy-loss spectroscopy measurements and theoretical calculations show that Mg(OH)₂ is a 6 eV direct-gap insulator in 2D, and its optical band gap displays strong band renormalization effects from monolayer to bulk, marking the first experimental confirmation of confinement effects in 2D insulators. Interestingly, 2D-Mg(OH)₂ sheets possess rather strong surface polarization (charge) effects which is in contrast to electrically neutral h-BN materials. Using 2D-Mg(OH)₂ sheets together with CVD-MoS₂ in the vertical stacking shows that a strong change transfer occurs from n-doped CVD-MoS₂ sheets to Mg(OH)₂, naturally depleting the semiconductor, pushing towards intrinsic doping limit and enhancing overall optical performance of 2D semiconductors. Results not only establish unusual confinement effects in 2D-Mg(OH)₂, but also offer novel 2D-insulating material with unique physical, vibrational, and chemical properties for potential applications in flexible optoelectronics.