Faculty and Staff

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.

This paper studies the effect of targeted observations on state and parameter estimates determined with Kalman filter data assimilation (DA) techniques. We first provide an analytical result demonstrating that targeting observations within the Kalman filter for a linear model can significantly reduce state estimation error as opposed to fixed or randomly located observations. We next conduct observing system simulation experiments for a chaotic model of meteorological interest, where we demonstrate that the local ensemble transform Kalman filter (LETKF) with targeted observations based on largest ensemble variance is skillful in providing more accurate state estimates than the LETKF with randomly located observations. Additionally, we find that a hybrid ensemble Kalman filter parameter estimation method accurately updates model parameters within the targeted observation context to further improve state estimation.

The effects of urbanization on ozone levels have been widely investigated over cities primarily located in temperate and/or humid regions. In this study, nested WRF-Chem simulations with a finest grid resolution of 1 km are conducted to investigate ozone concentrations O₃ due to urbanization within cities in arid/semi-arid environments. First, a method based on a shape preserving Monotonic Cubic Interpolation (MCI) is developed and used to downscale anthropogenic emissions from the 4 km resolution 2005 National Emissions Inventory (NEI05) to the finest model resolution of 1 km. Using the rapidly expanding Phoenix metropolitan region as the area of focus, we demonstrate the proposed MCI method achieves ozone simulation results with appreciably improved correspondence to observations relative to the default interpolation method of the WRF-Chem system. Next, two additional sets of experiments are conducted, with the recommended MCI approach, to examine impacts of urbanization on ozone production: (1) the urban land cover is included (i.e., urbanization experiments) and, (2) the urban land cover is replaced with the region's native shrubland. Impacts due to the presence of the built environment on O₃ are highly heterogeneous across the metropolitan area. Increased near surface O₃ due to urbanization of 10–20 ppb is predominantly a nighttime phenomenon while simulated impacts during daytime are negligible. Urbanization narrows the daily O₃ range (by virtue of increasing nighttime minima), an impact largely due to the region's urban heat island. Our results demonstrate the importance of the MCI method for accurate representation of the diurnal profile of ozone, and highlight its utility for high-resolution air quality simulations for urban areas.

Introduction: Fluorescence-guided surgery is one of the rapidly emerging methods of surgical “theranostics.” In this review, we summarize current fluorescence techniques used in neurosurgical practice for brain tumor patients as well as future applications of recent laboratory and translational studies.

Methods: Review of the literature.

Results: A wide spectrum of fluorophores that have been tested for brain surgery is reviewed. Beginning with a fluorescein sodium application in 1948 by Moore, fluorescence-guided brain tumor surgery is either routinely applied in some centers or is under active study in clinical trials. Besides the trinity of commonly used drugs (fluorescein sodium, 5-aminolevulinic acid, and indocyanine green), less studied fluorescent stains, such as tetracyclines, cancer-selective alkylphosphocholine analogs, cresyl violet, acridine orange, and acriflavine, can be used for rapid tumor detection and pathological tissue examination. Other emerging agents, such as activity-based probes and targeted molecular probes that can provide biomolecular specificity for surgical visualization and treatment, are reviewed. Furthermore, we review available engineering and optical solutions for fluorescent surgical visualization. Instruments for fluorescent-guided surgery are divided into wide-field imaging systems and hand-held probes. Recent advancements in quantitative fluorescence-guided surgery are discussed.

Conclusion: We are standing on the threshold of the era of marker-assisted tumor management. Innovations in the fields of surgical optics, computer image analysis, and molecular bioengineering are advancing fluorescence-guided tumor resection paradigms, leading to cell-level approaches to visualization and resection of brain tumors.

An interesting occurrence of a Rossby wave breaking event observed during the VORCORE experiment is presented and explained. Twenty-seven balloons were launched inside the Antarctic polar vortex. Almost all of these balloons evolved in the stratosphere around 500K within the vortex, except the one launched on 28 October 2005. In this case, the balloon was caught within a tongue of high potential vorticity (PV), and was ejected from the polar vortex. The evolution of this event is studied for the period between 19 and 25 November 2005. It is found that at the beginning of this period, the polar vortex experienced distortions due to the presence of Rossby waves. Then, these waves break and a tongue of high PV develops. On 25 November, the tongue became separated from the vortex and the balloon was ejected into the surf zone. Lagrangian simulations demonstrate that the air masses surrounding the balloon after its ejection were originating from the vortex edge. The wave breaking and the development of the tongue are confined within a region where a planetary Quasi-Stationary Wave 1 (QSW1) induces wind speeds with weaker values. The QSW1 causes asymmetry in the wind speed and the horizontal PV gradient along the edge of the polar vortex, resulting in a localized jet. Rossby waves with smaller scales propagating on top of this jet amplify as they enter the jet exit region and then break. The role of the QSW1 on the formation of the weak flow conditions that caused the non-linear wave breaking observed near the vortex edge is confirmed by three-dimensional numerical simulations using forcing with and without the contribution of the QSW1.

Background: The successful treatment of malignant gliomas remains a challenge despite the current standard of care, which consists of surgery, radiation and temozolomide. Advances in the survival of brain cancer patients require the design of new therapeutic approaches that take advantage of common phenotypes such as the altered metabolism found in cancer cells. It has therefore been postulated that the high-fat, low-carbohydrate, adequate protein ketogenic diet (KD) may be useful in the treatment of brain tumors. We have demonstrated that the KD enhances survival and potentiates standard therapy in a mouse model of malignant glioma, yet the mechanisms are not fully understood.

Methods: To explore the effects of the KD on various aspects of tumor growth and progression, we used the immunocompetent, syngeneic GL261-Luc2 mouse model of malignant glioma.

Results: Tumors from animals maintained on KD showed reduced expression of the hypoxia marker carbonic anhydrase 9, hypoxia inducible factor 1-alpha, and decreased activation of nuclear factor kappa B. Additionally, tumors from animals maintained on KD had reduced tumor microvasculature and decreased expression of vascular endothelial growth factor receptor 2, matrix metalloproteinase-2 and vimentin. Peritumoral edema was significantly reduced in animals fed the KD and protein analyses showed altered expression of zona occludens-1 and aquaporin-4.

Conclusions: The KD directly or indirectly alters the expression of several proteins involved in malignant progression and may be a useful tool for the treatment of gliomas.

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.

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.

Physical mechanisms of incongruency between observations and Weather Research and Forecasting (WRF) Model predictions are examined. Limitations of evaluation are constrained by (i) parameterizations of model physics, (ii) parameterizations of input data, (iii) model resolution, and (iv) flux observation resolution. Observations from a new 22.1-m flux tower situated within a residential neighborhood in Phoenix, Arizona, are utilized to evaluate the ability of the urbanized WRF to resolve finescale surface energy balance (SEB) when using the urban classes derived from the 30-m-resolution National Land Cover Database. Modeled SEB response to a large seasonal variation of net radiation forcing was tested during synoptically quiescent periods of high pressure in winter 2011 and premonsoon summer 2012. Results are presented from simulations employing five nested domains down to 333-m horizontal resolution. A comparative analysis of model cases testing parameterization of physical processes was done using four configurations of urban parameterization for the bulk urban scheme versus three representations with the Urban Canopy Model (UCM) scheme, and also for two types of planetary boundary layer parameterization: the local Mellor–Yamada–Janjić scheme and the nonlocal Yonsei University scheme. Diurnal variation in SEB constituent fluxes is examined in relation to surface-layer stability and modeled diagnostic variables. Improvement is found when adapting UCM for Phoenix with reduced errors in the SEB components. Finer model resolution is seen to have insignificant (<1 standard deviation) influence on mean absolute percent difference of 30-min diurnal mean SEB terms.