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.

Although insulin resistance in skeletal muscle is well-characterized, the role of circulating whole blood in the metabolic syndrome phenotype is not well understood. We set out to test the hypothesis that genes involved in inflammation, insulin signaling and mitochondrial function would be altered in expression in the whole blood of individuals with metabolic syndrome. We further wanted to examine whether similar relationships that we have found previously in skeletal muscle exist in peripheral whole blood cells. All subjects (n=184) were Latino descent from the Arizona Insulin Resistance registry. Subjects were classified based on the metabolic syndrome phenotype according to the National Cholesterol Education Program’s Adult Treatment Panel III. Of the 184 Latino subjects in the study, 74 were classified with the metabolic syndrome and 110 were without. Whole blood gene expression profiling was performed using the Agilent 4x44K Whole Human Genome Microarray. Whole blood microarray analysis identified 1,432 probes that were altered in expression ≥1.2 fold and P<0.05 after Benjamini-Hochberg in the metabolic syndrome subjects. KEGG pathway analysis revealed significant enrichment for pathways including ribosome, oxidative phosphorylation and MAPK signaling (all Benjamini-Hochberg P<0.05). Whole blood mRNA expression changes observed in the microarray data were confirmed by quantitative RT-PCR. Transcription factor binding motif enrichment analysis revealed E2F1, ELK1, NF-kappaB, STAT1 and STAT3 significantly enriched after Bonferroni correction (all P<0.05). The results of the present study demonstrate that whole blood is a useful tissue for studying the metabolic syndrome and its underlying insulin resistance although the relationship between blood and skeletal muscle differs.

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.

Type 2 diabetes (T2D) is a complex metabolic disease that is more prevalent in ethnic groups such as Mexican Americans, and is strongly associated with the risk factors obesity and insulin resistance. The goal of this study was to perform whole genome gene expression profiling in adipose tissue to detect common patterns of gene regulation associated with obesity and insulin resistance. We used phenotypic and genotypic data from 308 Mexican American participants from the Veterans Administration Genetic Epidemiology Study (VAGES). Basal fasting RNA was extracted from adipose tissue biopsies from a subset of 75 unrelated individuals, and gene expression data generated on the Illumina BeadArray platform. The number of gene probes with significant expression above baseline was approximately 31,000. We performed multiple regression analysis of all probes with 15 metabolic traits. Adipose tissue had 3,012 genes significantly associated with the traits of interest (false discovery rate, FDR ≤ 0.05). The significance of gene expression changes was used to select 52 genes with significant (FDR ≤ 10^-4) gene expression changes across multiple traits. Gene sets/Pathways analysis identified one gene, alcohol dehydrogenase 1B (ADH1B) that was significantly enriched (P < 10^-60) as a prime candidate for involvement in multiple relevant metabolic pathways. Illumina BeadChip derived ADH1B expression data was consistent with quantitative real time PCR data. We observed significant inverse correlations with waist circumference (2.8 x 10[superscript -9]), BMI (5.4 x 10^-6), and fasting plasma insulin (P < 0.001). These findings are consistent with a central role for ADH1B in obesity and insulin resistance and provide evidence for a novel genetic regulatory mechanism for human metabolic diseases related to these traits.

Our previous studies show reduced abundance of the β-subunit of mitochondrial H+-ATP synthase (β-F1-ATPase) in skeletal muscle of obese individuals. The β-F1-ATPase forms the catalytic core of the ATP synthase, and it is critical for ATP production in muscle. The mechanism(s) impairing β-F1-ATPase metabolism in obesity, however, are not completely understood. First, we studied total muscle protein synthesis and the translation efficiency of β-F1-ATPase in obese (BMI, 36±1 kg/m²) and lean (BMI, 22±1 kg/m²) subjects. Both total protein synthesis (0.044±0.006 vs 0.066±0.006%·h^-1) and translation efficiency of β-F1-ATPase (0.0031±0.0007 vs 0.0073±0.0004) were lower in muscle from the obese subjects when compared to the lean controls (P<0.05). We then evaluated these same responses in a primary cell culture model, and tested the specific hypothesis that circulating non-esterified fatty acids (NEFA) in obesity play a role in the responses observed in humans. The findings on total protein synthesis and translation efficiency of β-F1-ATPase in primary myotubes cultured from a lean subject, and after exposure to NEFA extracted from serum of an obese subject, were similar to those obtained in humans. Among candidate microRNAs (i.e., non-coding RNAs regulating gene expression), we identified miR-127-5p in preventing the production of β-F1-ATPase. Muscle expression of miR-127-5p negatively correlated with β-F1-ATPase protein translation efficiency in humans (r = – 0.6744; P<0.01), and could be modeled in vitro by prolonged exposure of primary myotubes derived from the lean subject to NEFA extracted from the obese subject. On the other hand, locked nucleic acid inhibitor synthesized to target miR-127-5p significantly increased β-F1-ATPase translation efficiency in myotubes (0.6±0.1 vs 1.3±0.3, in control vs exposure to 50 nM inhibitor; P<0.05). Our experiments implicate circulating NEFA in obesity in suppressing muscle protein metabolism, and establish impaired β-F1-ATPase translation as an important consequence of obesity.

Background: Obesity is a metabolic disease caused by environmental and genetic factors. However, the epigenetic mechanisms of obesity are incompletely understood. The aim of our study was to investigate the role of skeletal muscle DNA methylation in combination with transcriptomic changes in obesity.

Results: Muscle biopsies were obtained basally from lean (n = 12; BMI = 23.4 ± 0.7 kg/m[superscript 2]) and obese (n = 10; BMI = 32.9 ± 0.7 kg/m[superscript 2]) participants in combination with euglycemic-hyperinsulinemic clamps to assess insulin sensitivity. We performed reduced representation bisulfite sequencing (RRBS) next-generation methylation and microarray analyses on DNA and RNA isolated from vastus lateralis muscle biopsies. There were 13,130 differentially methylated cytosines (DMC; uncorrected P < 0.05) that were altered in the promoter and untranslated (5' and 3'UTR) regions in the obese versus lean analysis. Microarray analysis revealed 99 probes that were significantly (corrected P < 0.05) altered. Of these, 12 genes (encompassing 22 methylation sites) demonstrated a negative relationship between gene expression and DNA methylation. Specifically, sorbin and SH3 domain containing 3 (SORBS3) which codes for the adapter protein vinexin was significantly decreased in gene expression (fold change −1.9) and had nine DMCs that were significantly increased in methylation in obesity (methylation differences ranged from 5.0 to 24.4 %). Moreover, differentially methylated region (DMR) analysis identified a region in the 5'UTR (Chr.8:22,423,530–22,423,569) of SORBS3 that was increased in methylation by 11.2 % in the obese group. The negative relationship observed between DNA methylation and gene expression for SORBS3 was validated by a site-specific sequencing approach, pyrosequencing, and qRT-PCR. Additionally, we performed transcription factor binding analysis and identified a number of transcription factors whose binding to the differentially methylated sites or region may contribute to obesity.

Conclusions: These results demonstrate that obesity alters the epigenome through DNA methylation and highlights novel transcriptomic changes in SORBS3 in skeletal muscle.

Background: Although the effect of the fat mass and obesity-associated (FTO) gene on adiposity is well established, there is a lack of evidence whether physical activity (PA) modifies the effect of FTO variants on obesity in Latino populations. Therefore, the purpose of this study was to examine PA influences and interactive effects between FTO variants and PA on measures of adiposity in Latinos.

Results: After controlling for age and sex, participants who did not engage in regular PA exhibited higher BMI, fat mass, HC, and WC with statistical significance (P < 0.001). Although significant associations between the three FTO genotypes and adiposity measures were found, none of the FTO genotype by PA interaction assessments revealed nominally significant associations. However, several of such interactive influences exhibited considerable trend towards association.

Conclusions: These data suggest that adiposity measures are associated with PA and FTO variants in Latinos, but the impact of their interactive influences on these obesity measures appear to be minimal. Future studies with large sample sizes may help to determine whether individuals with specific FTO variants exhibit differential responses to PA interventions.