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.

Quantitative three-dimensional (3D) computed tomography (CT) imaging of living single cells enables orientation-independent morphometric analysis of the intricacies of cellular physiology. Since its invention, x-ray CT has become indispensable in the clinic for diagnostic and prognostic purposes due to its quantitative absorption-based imaging in true 3D that allows objects of interest to be viewed and measured from any orientation. However, x-ray CT has not been useful at the level of single cells because there is insufficient contrast to form an image. Recently, optical CT has been developed successfully for fixed cells, but this technology called Cell-CT is incompatible with live-cell imaging due to the use of stains, such as hematoxylin, that are not compatible with cell viability. We present a novel development of optical CT for quantitative, multispectral functional 4D (three spatial + one spectral dimension) imaging of living single cells. The method applied to immune system cells offers truly isotropic 3D spatial resolution and enables time-resolved imaging studies of cells suspended in aqueous medium. Using live-cell optical CT, we found a heterogeneous response to mitochondrial fission inhibition in mouse macrophages and differential basal remodeling of small (0.1 to 1 fl) and large (1 to 20 fl) nuclear and mitochondrial structures on a 20- to 30-s time scale in human myelogenous leukemia cells. Because of its robust 3D measurement capabilities, live-cell optical CT represents a powerful new tool in the biomedical research field.

Driven by an increasing number of studies demonstrating its relevance to a broad variety of disease states, the bioenergy production phenotype has been widely characterized at the bulk sample level. Its cell-to-cell variability, a key player associated with cancer cell survival and recurrence, however, remains poorly understood due to ensemble averaging of the current approaches. We present a technology platform for performing oxygen consumption and extracellular acidification measurements of several hundreds to 1,000 individual cells per assay, while offering simultaneous analysis of cellular communication effects on the energy production phenotype. The platform comprises two major components: a tandem optical sensor for combined oxygen and pH detection, and a microwell device for isolation and analysis of single and few cells in hermetically sealed sub-nanoliter chambers. Our approach revealed subpopulations of cells with aberrant energy production profiles and enables determination of cellular response variability to electron transfer chain inhibitors and ion uncouplers.

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.

The histone deacetylase (HDAC) inhibitor vorinostat has received significant attention in recent years as an ‘epigenetic’ drug used to treat solid tumors. However, its mechanisms of action are not entirely understood, particularly with regard to its interaction with the aberrations in 3D nuclear structure that accompany neoplastic progression. We investigated the impact of vorinostat on human esophageal epithelial cell lines derived from normal, metaplastic (pre-cancerous), and malignant tissue. Using a combination of novel optical computed tomography (CT)-based quantitative 3D absorption microscopy and conventional confocal fluorescence microscopy, we show that subjecting malignant cells to vorinostat preferentially alters their 3D nuclear architecture relative to non-cancerous cells. Optical CT (cell CT) imaging of fixed single cells showed that drug-treated cancer cells exhibit significant alterations in nuclear morphometry. Confocal microscopy revealed that vorinostat caused changes in the distribution of H3K9ac-marked euchromatin and H3K9me3-marked constitutive heterochromatin. Additionally, 3D immuno-FISH showed that drug-induced expression of the DNA repair gene MGMT was accompanied by spatial relocation toward the center of the nucleus in the nuclei of metaplastic but not in non-neoplastic cells. Our data suggest that vorinostat’s differential modulation of 3D nuclear architecture in normal and abnormal cells could play a functional role in its anti-cancer action.

Chronic manganese (Mn) exposure is associated with neuromotor and neurocognitive deficits, but the exact mechanism of Mn neurotoxicity is still unclear. With the advent of magnetic resonance imaging (MRI), in-vivo analysis of brain structures has become possible. Among different sub-cortical structures, the basal ganglia (BG) has been investigated as a putative anatomical biomarker in MR-based studies of Mn toxicity. However, previous investigations have yielded inconsistent results in terms of regional MR signal intensity changes. These discrepancies may be due to the subtlety of brain alterations caused by Mn toxicity, coupled to analysis techniques that lack the requisite detection power. Here, based on brain MRI, we apply a 3D surface-based morphometry method on 3 bilateral basal ganglia structures in school-age children chronically exposed to Mn through drinking water to investigate the effect of Mn exposure on brain anatomy. Our method successfully pinpointed significant enlargement of many areas of the basal ganglia structures, preferentially affecting the putamen. Moreover, these areas showed significant correlations with fine motor performance, indicating a possible link between altered basal ganglia neurodevelopment and declined motor performance in high Mn exposed children.

The apolipoprotein E (APOE) e4 allele is the most prevalent genetic risk factor for Alzheimer's disease (AD). Hippocampal volumes are generally smaller in AD patients carrying the e4 allele compared to e4 noncarriers. Here we examined the effect of APOE e4 on hippocampal morphometry in a large imaging database—the Alzheimer's Disease Neuroimaging Initiative (ADNI). We automatically segmented and constructed hippocampal surfaces from the baseline MR images of 725 subjects with known APOE genotype information including 167 with AD, 354 with mild cognitive impairment (MCI), and 204 normal controls. High-order correspondences between hippocampal surfaces were enforced across subjects with a novel inverse consistent surface fluid registration method. Multivariate statistics consisting of multivariate tensor-based morphometry (mTBM) and radial distance were computed for surface deformation analysis. Using Hotelling's T2 test, we found significant morphological deformation in APOE e4 carriers relative to noncarriers in the entire cohort as well as in the nondemented (pooled MCI and control) subjects, affecting the left hippocampus more than the right, and this effect was more pronounced in e4 homozygotes than heterozygotes. Our findings are consistent with previous studies that showed e4 carriers exhibit accelerated hippocampal atrophy; we extend these findings to a novel measure of hippocampal morphometry. Hippocampal morphometry has significant potential as an imaging biomarker of early stage AD.

Gompertz’s empirical equation remains the most popular one in describing cancer cell population growth in a wide spectrum of bio-medical situations due to its good fit to data and simplicity. Many efforts were documented in the literature aimed at understanding the mechanisms that may support Gompertz’s elegant model equation. One of the most convincing efforts was carried out by Gyllenberg and Webb. They divide the cancer cell population into the proliferative cells and the quiescent cells. In their two dimensional model, the dead cells are assumed to be removed from the tumor instantly. In this paper, we modify their model by keeping track of the dead cells remaining in the tumor. We perform mathematical and computational studies on this three dimensional model and compare the model dynamics to that of the model of Gyllenberg and Webb. Our mathematical findings suggest that if an avascular tumor grows according to our three-compartment model, then as the death rate of quiescent cells decreases to zero, the percentage of proliferative cells also approaches to zero. Moreover, a slow dying quiescent population will increase the size of the tumor. On the other hand, while the tumor size does not depend on the dead cell removal rate, its early and intermediate growth stages are very sensitive to it.

Background: Androgens bind to the androgen receptor (AR) in prostate cells and are essential survival factors for healthy prostate epithelium. Most untreated prostate cancers retain some dependence upon the AR and respond, at least transiently, to androgen ablation therapy. However, the relationship between endogenous androgen levels and cancer etiology is unclear. High levels of androgens have traditionally been viewed as driving abnormal proliferation leading to cancer, but it has also been suggested that low levels of androgen could induce selective pressure for abnormal cells. We formulate a mathematical model of androgen regulated prostate growth to study the effects of abnormal androgen levels on selection for pre-malignant phenotypes in early prostate cancer development.

Results: We find that cell turnover rate increases with decreasing androgen levels, which may increase the rate of mutation and malignant evolution. We model the evolution of a heterogeneous prostate cell population using a continuous state-transition model. Using this model we study selection for AR expression under different androgen levels and find that low androgen environments, caused either by low serum testosterone or by reduced 5α-reductase activity, select more strongly for elevated AR expression than do normal environments. High androgen actually slightly reduces selective pressure for AR upregulation. Moreover, our results suggest that an aberrant androgen environment may delay progression to a malignant phenotype, but result in a more dangerous cancer should one arise.

Conclusions: The model represents a useful initial framework for understanding the role of androgens in prostate cancer etiology, and it suggests that low androgen levels can increase selection for phenotypes resistant to hormonal therapy that may also be more aggressive. Moreover, clinical treatment with 5α-reductase inhibitors such as finasteride may increase the incidence of therapy resistant cancers.