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.

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: Immunosignaturing is a new peptide microarray based technology for profiling of humoral immune responses. Despite new challenges, immunosignaturing gives us the opportunity to explore new and fundamentally different research questions. In addition to classifying samples based on disease status, the complex patterns and latent factors underlying immunosignatures, which we attempt to model, may have a diverse range of applications.

Methods: We investigate the utility of a number of statistical methods to determine model performance and address challenges inherent in analyzing immunosignatures. Some of these methods include exploratory and confirmatory factor analyses, classical significance testing, structural equation and mixture modeling.

Results: We demonstrate an ability to classify samples based on disease status and show that immunosignaturing is a very promising technology for screening and presymptomatic screening of disease. In addition, we are able to model complex patterns and latent factors underlying immunosignatures. These latent factors may serve as biomarkers for disease and may play a key role in a bioinformatic method for antibody discovery.

Conclusion: Based on this research, we lay out an analytic framework illustrating how immunosignatures may be useful as a general method for screening and presymptomatic screening of disease as well as antibody discovery.

Background: Microarray image analysis processes scanned digital images of hybridized arrays to produce the input spot-level data for downstream analysis, so it can have a potentially large impact on those and subsequent analysis. Signal saturation is an optical effect that occurs when some pixel values for highly expressed genes or peptides exceed the upper detection threshold of the scanner software (2¹⁶ - 1 = 65, 535 for 16-bit images). In practice, spots with a sizable number of saturated pixels are often flagged and discarded. Alternatively, the saturated values are used without adjustments for estimating spot intensities. The resulting expression data tend to be biased downwards and can distort high-level analysis that relies on these data. Hence, it is crucial to effectively correct for signal saturation.

Results: We developed a flexible mixture model-based segmentation and spot intensity estimation procedure that accounts for saturated pixels by incorporating a censored component in the mixture model. As demonstrated with biological data and simulation, our method extends the dynamic range of expression data beyond the saturation threshold and is effective in correcting saturation-induced bias when the lost information is not tremendous. We further illustrate the impact of image processing on downstream classification, showing that the proposed method can increase diagnostic accuracy using data from a lymphoma cancer diagnosis study.

Conclusions: The presented method adjusts for signal saturation at the segmentation stage that identifies a pixel as part of the foreground, background or other. The cluster membership of a pixel can be altered versus treating saturated values as truly observed. Thus, the resulting spot intensity estimates may be more accurate than those obtained from existing methods that correct for saturation based on already segmented data. As a model-based segmentation method, our procedure is able to identify inner holes, fuzzy edges and blank spots that are common in microarray images. The approach is independent of microarray platform and applicable to both single- and dual-channel microarrays.

Background: High-throughput technologies such as DNA, RNA, protein, antibody and peptide microarrays are often used to examine differences across drug treatments, diseases, transgenic animals, and others. Typically one trains a classification system by gathering large amounts of probe-level data, selecting informative features, and classifies test samples using a small number of features. As new microarrays are invented, classification systems that worked well for other array types may not be ideal. Expression microarrays, arguably one of the most prevalent array types, have been used for years to help develop classification algorithms. Many biological assumptions are built into classifiers that were designed for these types of data. One of the more problematic is the assumption of independence, both at the probe level and again at the biological level. Probes for RNA transcripts are designed to bind single transcripts. At the biological level, many genes have dependencies across transcriptional pathways where co-regulation of transcriptional units may make many genes appear as being completely dependent. Thus, algorithms that perform well for gene expression data may not be suitable when other technologies with different binding characteristics exist. The immunosignaturing microarray is based on complex mixtures of antibodies binding to arrays of random sequence peptides. It relies on many-to-many binding of antibodies to the random sequence peptides. Each peptide can bind multiple antibodies and each antibody can bind multiple peptides. This technology has been shown to be highly reproducible and appears promising for diagnosing a variety of disease states. However, it is not clear what is the optimal classification algorithm for analyzing this new type of data.

Results: We characterized several classification algorithms to analyze immunosignaturing data. We selected several datasets that range from easy to difficult to classify, from simple monoclonal binding to complex binding patterns in asthma patients. We then classified the biological samples using 17 different classification algorithms. Using a wide variety of assessment criteria, we found ‘Naïve Bayes’ far more useful than other widely used methods due to its simplicity, robustness, speed and accuracy.

Conclusions: ‘Naïve Bayes’ algorithm appears to accommodate the complex patterns hidden within multilayered immunosignaturing microarray data due to its fundamental mathematical properties.

Introduction: The ketogenic diet (KD) is a high-fat, low-carbohydrate diet that alters metabolism by increasing the level of ketone bodies in the blood. KetoCal® (KC) is a nutritionally complete, commercially available 4∶1 (fat∶ carbohydrate+protein) ketogenic formula that is an effective non-pharmacologic treatment for the management of refractory pediatric epilepsy. Diet-induced ketosis causes changes to brain homeostasis that have potential for the treatment of other neurological diseases such as malignant gliomas.

Methods: We used an intracranial bioluminescent mouse model of malignant glioma. Following implantation animals were maintained on standard diet (SD) or KC. The mice received 2×4 Gy of whole brain radiation and tumor growth was followed by in vivo imaging.

Results: Animals fed KC had elevated levels of β-hydroxybutyrate (p = 0.0173) and an increased median survival of approximately 5 days relative to animals maintained on SD. KC plus radiation treatment were more than additive, and in 9 of 11 irradiated animals maintained on KC the bioluminescent signal from the tumor cells diminished below the level of detection (p<0.0001). Animals were switched to SD 101 days after implantation and no signs of tumor recurrence were seen for over 200 days.

Conclusions: KC significantly enhances the anti-tumor effect of radiation. This suggests that cellular metabolic alterations induced through KC may be useful as an adjuvant to the current standard of care for the treatment of human malignant gliomas.

Immunosignaturing shows promise as a general approach to diagnosis. It has been shown to detect immunological signs of infection early during the course of disease and to distinguish Alzheimer’s disease from healthy controls. Here we test whether immunosignatures correspond to clinical classifications of disease using samples from people with brain tumors. Blood samples from patients undergoing craniotomies for therapeutically naïve brain tumors with diagnoses of astrocytoma (23 samples), Glioblastoma multiforme (22 samples), mixed oligodendroglioma/astrocytoma (16 samples), oligodendroglioma (18 samples), and 34 otherwise healthy controls were tested by immunosignature. Because samples were taken prior to adjuvant therapy, they are unlikely to be perturbed by non-cancer related affects. The immunosignaturing platform distinguished not only brain cancer from controls, but also pathologically important features about the tumor including type, grade, and the presence or absence of O⁶-methyl-guanine-DNA methyltransferase methylation promoter (MGMT), an important biomarker that predicts response to temozolomide in Glioblastoma multiformae patients.

Background: Malignant brain tumors affect people of all ages and are the second leading cause of cancer deaths in children. While current treatments are effective and improve survival, there remains a substantial need for more efficacious therapeutic modalities. The ketogenic diet (KD) - a high-fat, low-carbohydrate treatment for medically refractory epilepsy - has been suggested as an alternative strategy to inhibit tumor growth by altering intrinsic metabolism, especially by inducing glycopenia.

Methods: Here, we examined the effects of an experimental KD on a mouse model of glioma, and compared patterns of gene expression in tumors vs. normal brain from animals fed either a KD or a standard diet.

Results: Animals received intracranial injections of bioluminescent GL261-luc cells and tumor growth was followed in vivo. KD treatment significantly reduced the rate of tumor growth and prolonged survival. Further, the KD reduced reactive oxygen species (ROS) production in tumor cells. Gene expression profiling demonstrated that the KD induces an overall reversion to expression patterns seen in non-tumor specimens. Notably, genes involved in modulating ROS levels and oxidative stress were altered, including those encoding cyclooxygenase 2, glutathione peroxidases 3 and 7, and periredoxin 4.

Conclusions: Our data demonstrate that the KD improves survivability in our mouse model of glioma, and suggests that the mechanisms accounting for this protective effect likely involve complex alterations in cellular metabolism beyond simply a reduction in glucose.

Astronauts are exposed to a unique combination of stressors during spaceflight, which leads to alterations in their physiology and potentially increases their susceptibility to disease, including infectious diseases. To evaluate the potential impact of the spaceflight environment on the regulation of molecular pathways mediating cellular stress responses, we performed a first-of-its-kind pilot study to assess spaceflight-related gene-expression changes in the whole blood of astronauts. Using an array comprised of 234 well-characterized stress-response genes, we profiled transcriptomic changes in six astronauts (four men and two women) from blood preserved before and immediately following the spaceflight. Differentially regulated transcripts included those important for DNA repair, oxidative stress, and protein folding/degradation, including HSP90AB1, HSP27, GPX1, XRCC1, BAG-1, HHR23A, FAP48, and C-FOS. No gender-specific differences or relationship to number of missions flown was observed. This study provides a first assessment of transcriptomic changes occurring in the whole blood of astronauts in response to spaceflight.

There is an increasing awareness that health care must move from post-symptomatic treatment to presymptomatic intervention. An ideal system would allow regular inexpensive monitoring of health status using circulating antibodies to report on health fluctuations. Recently, we demonstrated that peptide microarrays can do this through antibody signatures (immunosignatures). Unfortunately, printed microarrays are not scalable. Here we demonstrate a platform based on fabricating microarrays (~10 M peptides per slide, 330,000 peptides per assay) on silicon wafers using equipment common to semiconductor manufacturing. The potential of these microarrays for comprehensive health monitoring is verified through the simultaneous detection and classification of six different infectious diseases and six different cancers. Besides diagnostics, these high-density peptide chips have numerous other applications both in health care and elsewhere.