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 number and variety of connectivity estimation methods is likely to continue to grow over the coming decade. Comparisons between methods are necessary to prune this growth to only the most accurate and robust methods. However, the nature of connectivity is elusive with different methods potentially attempting to identify different aspects of connectivity. Commonalities of connectivity definitions across methods upon which base direct comparisons can be difficult to derive. Here, we explicitly define “effective connectivity” using a common set of observation and state equations that are appropriate for three connectivity methods: dynamic causal modeling (DCM), multivariate autoregressive modeling (MAR), and switching linear dynamic systems for fMRI (sLDSf). In addition while deriving this set, we show how many other popular functional and effective connectivity methods are actually simplifications of these equations. We discuss implications of these connections for the practice of using one method to simulate data for another method. After mathematically connecting the three effective connectivity methods, simulated fMRI data with varying numbers of regions and task conditions is generated from the common equation. This simulated data explicitly contains the type of the connectivity that the three models were intended to identify. Each method is applied to the simulated data sets and the accuracy of parameter identification is analyzed. All methods perform above chance levels at identifying correct connectivity parameters. The sLDSf method was superior in parameter estimation accuracy to both DCM and MAR for all types of comparisons.

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.

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.

Purpose: PET (positron emission tomography) imaging researches of functional metabolism using fluorodeoxyglucose ([superscript 18]F-FDG) of animal brain are important in neuroscience studies. FDG-PET imaging studies are often performed on groups of rats, so it is desirable to establish an objective voxel-based statistical methodology for group data analysis.

Material and Methods: This study establishes a statistical parametric mapping (SPM) toolbox (plug-ins) named spmratIHEP for voxel-wise analysis of FDG-PET images of rat brain, in which an FDG-PET template and an intracranial mask image of rat brain in Paxinos & Watson space were constructed, and the default settings were modified according to features of rat brain. Compared to previous studies, our constructed rat brain template comprises not only the cerebrum and cerebellum, but also the whole olfactory bulb which made the later cognitive studies much more exhaustive. And with an intracranial mask image in the template space, the brain tissues of individuals could be extracted automatically. Moreover, an atlas space is used for anatomically labeling the functional findings in the Paxinos & Watson space. In order to standardize the template image with the atlas accurately, a synthetic FDG-PET image with six main anatomy structures is constructed from the atlas, which performs as a target image in the co-registration.

Results: The spatial normalization procedure is evaluated, by which the individual rat brain images could be standardized into the Paxinos & Watson space successfully and the intracranial tissues could also be extracted accurately. The practical usability of this toolbox is evaluated using FDG-PET functional images from rats with left side middle cerebral artery occlusion (MCAO) in comparison to normal control rats. And the two-sample t-test statistical result is almost related to the left side MCA.

Conclusion: We established a toolbox of SPM8 named spmratIHEP for voxel-wise analysis of FDG-PET images of rat brain.

Background: Carriers of the APOE ε4 allele are at increased risk of developing Alzheimer’s disease (AD), and have been shown to have reduced cerebral metabolic rate of glucose (CMRgl) in the same brain areas frequently affected in AD. These individuals also exhibit reduced plasma levels of apolipoprotein E (apoE) attributed to a specific decrease in the apoE4 isoform as determined by quantification of individual apoE isoforms in APOE ε4 heterozygotes. Whether low plasma apoE levels are associated with structural and functional brain measurements and cognitive performance remains to be investigated.

Methods: Using quantitative mass spectrometry we quantified the plasma levels of total apoE and the individual apoE3 and apoE4 isoforms in 128 cognitively normal APOE ε3/ε4 individuals included in the Arizona APOE cohort. All included individuals had undergone extensive neuropsychological testing and 25 had in addition undergone FDG-PET and MRI to determine CMRgl and regional gray matter volume (GMV).

Results: Our results demonstrated higher apoE4 levels in females versus males and an age-dependent increase in the apoE3 isoform levels in females only. Importantly, a higher relative ratio of apoE4 over apoE3 was associated with GMV loss in the right posterior cingulate and with reduced CMRgl bilaterally in the anterior cingulate and in the right hippocampal area. Additional exploratory analysis revealed several negative associations between total plasma apoE, individual apoE isoform levels, GMV and CMRgl predominantly in the frontal, occipital and temporal areas. Finally, our results indicated only weak associations between apoE plasma levels and cognitive performance which further appear to be affected by sex.

Conclusions: Our study proposes a sex-dependent and age-dependent variation in plasma apoE isoform levels and concludes that peripheral apoE levels are associated with GMV, CMRgl and possibly cognitive performance in cognitively healthy individuals with a genetic predisposition to AD.

Background: We introduced a hypometabolic convergence index (HCI) to characterize in a single measurement the extent to which a person’s fluorodeoxyglucose positron emission tomogram (FDG PET) corresponds to that in Alzheimer’s disease (AD). Apolipoprotein E ε4 (APOE ε4) gene dose is associated with three levels of risk for late-onset AD. We explored the association between gene dose and HCI in cognitively normal ε4 homozygotes, heterozygotes, and non-carriers.

Methods: An algorithm was used to characterize and compare AD-related HCIs in cognitively normal individuals, including 36 ε4 homozygotes, 46 heterozygotes, and 78 non-carriers.

Results: These three groups differed significantly in their HCIs (ANOVA, p = 0.004), and there was a significant association between HCIs and gene dose (linear trend, p = 0.001).

Conclusions: The HCI is associated with three levels of genetic risk for late-onset AD. This supports the possibility of using a single FDG PET measurement to help in the preclinical detection and tracking of AD.

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.