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.

To address the need to study frozen clinical specimens using next-generation RNA, DNA, chromatin immunoprecipitation (ChIP) sequencing and protein analyses, we developed a biobank work flow to prospectively collect biospecimens from patients with renal cell carcinoma (RCC). We describe our standard operating procedures and work flow to annotate pathologic results and clinical outcomes. We report quality control outcomes and nucleic acid yields of our RCC submissions (N=16) to The Cancer Genome Atlas (TCGA) project, as well as newer discovery platforms, by describing mass spectrometry analysis of albumin oxidation in plasma and 6 ChIP sequencing libraries generated from nephrectomy specimens after histone H3 lysine 36 trimethylation (H3K36me3) immunoprecipitation. From June 1, 2010, through January 1, 2013, we enrolled 328 patients with RCC. Our mean (SD) TCGA RNA integrity numbers (RINs) were 8.1 (0.8) for papillary RCC, with a 12.5% overall rate of sample disqualification for RIN <7. Banked plasma had significantly less albumin oxidation (by mass spectrometry analysis) than plasma kept at 25°C (P<.001). For ChIP sequencing, the FastQC score for average read quality was at least 30 for 91% to 95% of paired-end reads. In parallel, we analyzed frozen tissue by RNA sequencing; after genome alignment, only 0.2% to 0.4% of total reads failed the default quality check steps of Bowtie2, which was comparable to the disqualification ratio (0.1%) of the 786-O RCC cell line that was prepared under optimal RNA isolation conditions. The overall correlation coefficients for gene expression between Mayo Clinic vs TCGA tissues ranged from 0.75 to 0.82. These data support the generation of high-quality nucleic acids for genomic analyses from banked RCC. Importantly, the protocol does not interfere with routine clinical care. Collections over defined time points during disease treatment further enhance collaborative efforts to integrate genomic information with outcomes.

Cancer therapy selects for cancer cells resistant to treatment, a process that is fundamentally evolutionary. To what extent, however, is the evolutionary perspective employed in research on therapeutic resistance and relapse? We analyzed 6,228 papers on therapeutic resistance and/or relapse in cancers and found that the use of evolution terms in abstracts has remained at about 1% since the 1980s. However, detailed coding of 22 recent papers revealed a higher proportion of papers using evolutionary methods or evolutionary theory, although this number is still less than 10%. Despite the fact that relapse and therapeutic resistance is essentially an evolutionary process, it appears that this framework has not permeated research. This represents an unrealized opportunity for advances in research on therapeutic resistance.

Insulin-like growth factor 1 (IGF1) is an important biomarker for the management of growth hormone disorders. Recently there has been rising interest in deploying mass spectrometric (MS) methods of detection for measuring IGF1. However, widespread clinical adoption of any MS-based IGF1 assay will require increased throughput and speed to justify the costs of analyses, and robust industrial platforms that are reproducible across laboratories. Presented here is an MS-based quantitative IGF1 assay with performance rating of >1,000 samples/day, and a capability of quantifying IGF1 point mutations and posttranslational modifications. The throughput of the IGF1 mass spectrometric immunoassay (MSIA) benefited from a simplified sample preparation step, IGF1 immunocapture in a tip format, and high-throughput MALDI-TOF MS analysis. The Limit of Detection and Limit of Quantification of the resulting assay were 1.5 μg/L and 5 μg/L, respectively, with intra- and inter-assay precision CVs of less than 10%, and good linearity and recovery characteristics. The IGF1 MSIA was benchmarked against commercially available IGF1 ELISA via Bland-Altman method comparison test, resulting in a slight positive bias of 16%. The IGF1 MSIA was employed in an optimized parallel workflow utilizing two pipetting robots and MALDI-TOF-MS instruments synced into one-hour phases of sample preparation, extraction and MSIA pipette tip elution, MS data collection, and data processing. Using this workflow, high-throughput IGF1 quantification of 1,054 human samples was achieved in approximately 9 hours. This rate of assaying is a significant improvement over existing MS-based IGF1 assays, and is on par with that of the enzyme-based immunoassays. Furthermore, a mutation was detected in ∼1% of the samples (SNP: rs17884626, creating an A→T substitution at position 67 of the IGF1), demonstrating the capability of IGF1 MSIA to detect point mutations and posttranslational modifications.

Serum Amyloid A (SAA) is an acute phase protein complex consisting of several abundant isoforms. The N- terminus of SAA is critical to its function in amyloid formation. SAA is frequently truncated, either missing an arginine or an arginine-serine dipeptide, resulting in isoforms that may influence the capacity to form amyloid. However, the relative abundance of truncated SAA in diabetes and chronic kidney disease is not known.

Methods: Using mass spectrometric immunoassay, the abundance of SAA truncations relative to the native variants was examined in plasma of 91 participants with type 2 diabetes and chronic kidney disease and 69 participants without diabetes.

Results: The ratio of SAA 1.1 (missing N-terminal arginine) to native SAA 1.1 was lower in diabetics compared to non-diabetics (p = 0.004), and in males compared to females (p<0.001). This ratio was negatively correlated with glycated hemoglobin (r = −0.32, p<0.001) and triglyceride concentrations (r = −0.37, p<0.001), and positively correlated with HDL cholesterol concentrations (r = 0.32, p<0.001).

Conclusion: The relative abundance of the N-terminal arginine truncation of SAA1.1 is significantly decreased in diabetes and negatively correlates with measures of glycemic and lipid control.

Background: Cysteine sulfenic acid (Cys-SOH) plays important roles in the redox regulation of numerous proteins. As a relatively unstable posttranslational protein modification it is difficult to quantify the degree to which any particular protein is modified by Cys-SOH within a complex biological environment. The goal of these studies was to move a step beyond detection and into the relative quantification of Cys-SOH within specific proteins found in a complex biological setting--namely, human plasma.

Results: This report describes the possibilities and limitations of performing such analyses based on the use of thionitrobenzoic acid and dimedone-based probes which are commonly employed to trap Cys-SOH. Results obtained by electrospray ionization-based mass spectrometric immunoassay reveal the optimal type of probe for such analyses as well as the reproducible relative quantification of Cys-SOH within albumin and transthyretin extracted from human plasma--the latter as a protein previously unknown to be modified by Cys-SOH.

Conclusions: The relative quantification of Cys-SOH within specific proteins in a complex biological setting can be accomplished, but several analytical precautions related to trapping, detecting, and quantifying Cys-SOH must be taken into account prior to pursuing its study in such matrices.

Land surface energy balance in a built environment is widely modelled using urban canopy models with representation of building arrays as big street canyons. Modification of this simplified geometric representation, however, leads to challenging numerical difficulties in improving physical parameterization schemes that are deterministic in nature. In this paper, we develop a stochastic algorithm to estimate view factors between canyon facets in the presence of shade trees based on Monte Carlo simulation, where an analytical formulation is inhibited by the complex geometry. The model is validated against analytical solutions of benchmark radiative problems as well as field measurements in real street canyons. In conjunction with the matrix method resolving infinite number of reflections, the proposed model is capable of predicting the radiative exchange inside the street canyon with good accuracy. Modeling of transient evolution of thermal filed inside the street canyon using the proposed method demonstrate the potential of shade trees in mitigating canyon surface temperatures as well as saving of building energy use. This new numerical framework also deepens our insight into the fundamental physics of radiative heat transfer and surface energy balance for urban climate modeling.

Studies on urban heat island (UHI) have been more than a century after the phenomenon was first discovered in the early 1800s. UHI emerges as the source of many urban environmental problems and exacerbates the living environment in cities. Under the challenges of increasing urbanization and future climate changes, there is a pressing need for sustainable adaptation/mitigation strategies for UHI effects, one popular option being the use of reflective materials. While it is introduced as an effective method to reduce temperature and energy consumption in cities, its impacts on environmental sustainability and large-scale non-local effect are inadequately explored. This paper provides a synthetic overview of potential environmental impacts of reflective materials at a variety of scales, ranging from energy load on a single building to regional hydroclimate. The review shows that mitigation potential of reflective materials depends on a set of factors, including building characteristics, urban environment, meteorological and geographical conditions, to name a few. Precaution needs to be exercised by city planners and policy makers for large-scale deployment of reflective materials before their environmental impacts, especially on regional hydroclimates, are better understood. In general, it is recommended that optimal strategy for UHI needs to be determined on a city-by-city basis, rather than adopting a “one-solution-fits-all” strategy.

The hysteresis effect in diurnal cycles of net radiation R-n and ground heat flux G(0) has been observed in many studies, while the governing mechanism remains vague. In this study, we link the phenomenology of hysteresis loops to the wave phase difference between the diurnal evolutions of various terms in the surface energy balance. R-n and G(0) are parameterized with the incoming solar radiation and the surface temperature as two control parameters of the surface energy partitioning. The theoretical analysis shows that the vertical water flux W and the scaled ratio A(s)*/A(T)* (net shortwave radiation to outgoing longwave radiation) play crucial roles in shaping hysteresis loops of R-n and G(0). Comparisons to field measurements indicate that hysteresis loops for different land covers can be well captured by the theoretical model, which is also consistent with Camuffo-Bernadi formula. This study provides insight into the surface partitioning and temporal evolution of the energy budget at the land surface.