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.

In vitro models that mimic in vivo host-pathogen interactions are needed to evaluate candidate drugs that inhibit bacterial virulence traits. We established a new approach to study Pseudomonas aeruginosa biofilm susceptibility on biotic surfaces, using a three-dimensional (3-D) lung epithelial cell model. P. aeruginosa formed antibiotic resistant biofilms on 3-D cells without affecting cell viability. The biofilm-inhibitory activity of antibiotics and/or the anti-biofilm peptide DJK-5 were evaluated on 3-D cells compared to a plastic surface, in medium with and without fetal bovine serum (FBS). In both media, aminoglycosides were more efficacious in the 3-D cell model. In serum-free medium, most antibiotics (except polymyxins) showed enhanced efficacy when 3-D cells were present. In medium with FBS, colistin was less efficacious in the 3-D cell model. DJK-5 exerted potent inhibition of P. aeruginosa association with both substrates, only in serum-free medium. DJK-5 showed stronger inhibitory activity against P. aeruginosa associated with plastic compared to 3-D cells. The combined addition of tobramycin and DJK-5 exhibited more potent ability to inhibit P. aeruginosa association with both substrates. In conclusion, lung epithelial cells influence the efficacy of most antimicrobials against P. aeruginosa biofilm formation, which in turn depends on the presence or absence of FBS.

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.

This study presents the first global transcriptional profiling and phenotypic characterization of the major human opportunistic fungal pathogen, Candida albicans, grown in spaceflight conditions. Microarray analysis revealed that C. albicans subjected to short-term spaceflight culture differentially regulated 452 genes compared to synchronous ground controls, which represented 8.3% of the analyzed ORFs. Spaceflight-cultured C. albicans–induced genes involved in cell aggregation (similar to flocculation), which was validated by microscopic and flow cytometry analysis. We also observed enhanced random budding of spaceflight-cultured cells as opposed to bipolar budding patterns for ground samples, in accordance with the gene expression data. Furthermore, genes involved in antifungal agent and stress resistance were differentially regulated in spaceflight, including induction of ABC transporters and members of the major facilitator family, downregulation of ergosterol-encoding genes, and upregulation of genes involved in oxidative stress resistance.

Finally, downregulation of genes involved in actin cytoskeleton was observed. Interestingly, the transcriptional regulator Cap1 and over 30% of the Cap1 regulon was differentially expressed in spaceflight-cultured C. albicans. A potential role for Cap1 in the spaceflight response of C. albicans is suggested, as this regulator is involved in random budding, cell aggregation, and oxidative stress resistance; all related to observed spaceflight-associated changes of C. albicans. While culture of C. albicans in microgravity potentiates a global change in gene expression that could induce a virulence-related phenotype, no increased virulence in a murine intraperitoneal (i.p.) infection model was observed under the conditions of this study. Collectively, our data represent an important basis for the assessment of the risk that commensal flora could play during human spaceflight missions. Furthermore, since the low fluid-shear environment of microgravity is relevant to physical forces encountered by pathogens during the infection process, insights gained from this study could identify novel infectious disease mechanisms, with downstream benefits for the general public.

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.

The probiotic effects of Lactobacillus reuteri have been speculated to partly depend on its capacity to produce the antimicrobial substance reuterin during the reduction of glycerol in the gut. In this study, the potential of this process to protect human intestinal epithelial cells against infection with Salmonella enterica serovar Typhimurium was investigated. We used a three-dimensional (3-D) organotypic model of human colonic epithelium that was previously validated and applied to study interactions between S. Typhimurium and the intestinal epithelium that lead to enteric salmonellosis. Using this model system, we show that L. reuteri protects the intestinal cells against the early stages of Salmonella infection and that this effect is significantly increased when L. reuteri is stimulated to produce reuterin from glycerol. More specifically, the reuterin-containing ferment of L. reuteri caused a reduction in Salmonella adherence and invasion (1 log unit), and intracellular survival (2 log units). In contrast, the L. reuteri ferment without reuterin stimulated growth of the intracellular Salmonella population with 1 log unit. The short-term exposure to reuterin or the reuterin-containing ferment had no observed negative impact on intestinal epithelial cell health. However, long-term exposure (24 h) induced a complete loss of cell-cell contact within the epithelial aggregates and compromised cell viability. Collectively, these results shed light on a potential role for reuterin in inhibiting Salmonella-induced intestinal infections and may support the combined application of glycerol and L. reuteri. While future in vitro and in vivo studies of reuterin on intestinal health should fine-tune our understanding of the mechanistic effects, in particular in the presence of a complex gut microbiota, this the first report of a reuterin effect on the enteric infection process in any mammalian cell type.