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.

Background: Carriers of the APOE ε4 allele are at increased risk of developing Alzheimer’s disease (AD), and have been shown to have reduced cerebral metabolic rate of glucose (CMRgl) in the same brain areas frequently affected in AD. These individuals also exhibit reduced plasma levels of apolipoprotein E (apoE) attributed to a specific decrease in the apoE4 isoform as determined by quantification of individual apoE isoforms in APOE ε4 heterozygotes. Whether low plasma apoE levels are associated with structural and functional brain measurements and cognitive performance remains to be investigated.

Methods: Using quantitative mass spectrometry we quantified the plasma levels of total apoE and the individual apoE3 and apoE4 isoforms in 128 cognitively normal APOE ε3/ε4 individuals included in the Arizona APOE cohort. All included individuals had undergone extensive neuropsychological testing and 25 had in addition undergone FDG-PET and MRI to determine CMRgl and regional gray matter volume (GMV).

Results: Our results demonstrated higher apoE4 levels in females versus males and an age-dependent increase in the apoE3 isoform levels in females only. Importantly, a higher relative ratio of apoE4 over apoE3 was associated with GMV loss in the right posterior cingulate and with reduced CMRgl bilaterally in the anterior cingulate and in the right hippocampal area. Additional exploratory analysis revealed several negative associations between total plasma apoE, individual apoE isoform levels, GMV and CMRgl predominantly in the frontal, occipital and temporal areas. Finally, our results indicated only weak associations between apoE plasma levels and cognitive performance which further appear to be affected by sex.

Conclusions: Our study proposes a sex-dependent and age-dependent variation in plasma apoE isoform levels and concludes that peripheral apoE levels are associated with GMV, CMRgl and possibly cognitive performance in cognitively healthy individuals with a genetic predisposition to AD.

Background: We introduced a hypometabolic convergence index (HCI) to characterize in a single measurement the extent to which a person’s fluorodeoxyglucose positron emission tomogram (FDG PET) corresponds to that in Alzheimer’s disease (AD). Apolipoprotein E ε4 (APOE ε4) gene dose is associated with three levels of risk for late-onset AD. We explored the association between gene dose and HCI in cognitively normal ε4 homozygotes, heterozygotes, and non-carriers.

Methods: An algorithm was used to characterize and compare AD-related HCIs in cognitively normal individuals, including 36 ε4 homozygotes, 46 heterozygotes, and 78 non-carriers.

Results: These three groups differed significantly in their HCIs (ANOVA, p = 0.004), and there was a significant association between HCIs and gene dose (linear trend, p = 0.001).

Conclusions: The HCI is associated with three levels of genetic risk for late-onset AD. This supports the possibility of using a single FDG PET measurement to help in the preclinical detection and tracking of AD.

Predicting the timing of a castrate resistant prostate cancer is critical to lowering medical costs and improving the quality of life of advanced prostate cancer patients. We formulate, compare and analyze two mathematical models that aim to forecast future levels of prostate-specific antigen (PSA). We accomplish these tasks by employing clinical data of locally advanced prostate cancer patients undergoing androgen deprivation therapy (ADT). While these models are simplifications of a previously published model, they fit data with similar accuracy and improve forecasting results. Both models describe the progression of androgen resistance. Although Model 1 is simpler than the more realistic Model 2, it can fit clinical data to a greater precision. However, we found that Model 2 can forecast future PSA levels more accurately. These findings suggest that including more realistic mechanisms of androgen dynamics in a two population model may help androgen resistance timing prediction.

The Dawn Framing Camera (FC) has imaged the northern hemisphere of the Asteroid (4) Vesta at high spatial resolution and coverage. This study represents the first investigation of the overall geology of the northern hemisphere (22–90°N, quadrangles Av-1, 2, 3, 4 and 5) using these unique Dawn mission observations. We have compiled a morphologic map and performed crater size–frequency distribution (CSFD) measurements to date the geologic units. The hemisphere is characterized by a heavily cratered surface with a few highly subdued basins up to ∼200 km in diameter. The most widespread unit is a plateau (cratered highland unit), similar to, although of lower elevation than the equatorial Vestalia Terra plateau. Large-scale troughs and ridges have regionally affected the surface. Between ∼180°E and ∼270°E, these tectonic features are well developed and related to the south pole Veneneia impact (Saturnalia Fossae trough unit), elsewhere on the hemisphere they are rare and subdued (Saturnalia Fossae cratered unit). In these pre-Rheasilvia units we observed an unexpectedly high frequency of impact craters up to ∼10 km in diameter, whose formation could in part be related to the Rheasilvia basin-forming event. The Rheasilvia impact has potentially affected the northern hemisphere also with S–N small-scale lineations, but without covering it with an ejecta blanket. Post-Rheasilvia impact craters are small (<60 km in diameter) and show a wide range of degradation states due to impact gardening and mass wasting processes. Where fresh, they display an ejecta blanket, bright rays and slope movements on walls. In places, crater rims have dark material ejecta and some crater floors are covered by ponded material interpreted as impact melt.

Oppia Quadrangle Av-10 (288–360°E, ±22°) is a junction of key geologic features that preserve a rough history of Asteroid (4) Vesta and serves as a case study of using geologic mapping to define a relative geologic timescale. Clear filter images, stereo-derived topography, slope maps, and multispectral color-ratio images from the Framing Camera on NASA’s Dawn spacecraft served as basemaps to create a geologic map and investigate the spatial and temporal relationships of the local stratigraphy. Geologic mapping reveals the oldest map unit within Av-10 is the cratered highlands terrain which possibly represents original crustal material on Vesta that was then excavated by one or more impacts to form the basin Feralia Planitia. Saturnalia Fossae and Divalia Fossae ridge and trough terrains intersect the wall of Feralia Planitia indicating that this impact basin is older than both the Veneneia and Rheasilvia impact structures, representing Pre-Veneneian crustal material. Two of the youngest geologic features in Av-10 are Lepida (∼45 km diameter) and Oppia (∼40 km diameter) impact craters that formed on the northern and southern wall of Feralia Planitia and each cross-cuts a trough terrain. The ejecta blanket of Oppia is mapped as ‘dark mantle’ material because it appears dark orange in the Framing Camera ‘Clementine-type’ color-ratio image and has a diffuse, gradational contact distributed to the south across the rim of Rheasilvia. Mapping of surface material that appears light orange in color in the Framing Camera ‘Clementine-type’ color-ratio image as ‘light mantle material’ supports previous interpretations of an impact ejecta origin. Some light mantle deposits are easily traced to nearby source craters, but other deposits may represent distal ejecta deposits (emplaced >5 crater radii away) in a microgravity environment.

The Quadrangles Av-11 and Av-12 on Vesta are located at the northern rim of the giant Rheasilvia south polar impact basin. The primary geologic units in Av-11 and Av-12 include material from the Rheasilvia impact basin formation, smooth material and different types of impact crater structures (such as bimodal craters, dark and bright crater ray material and dark ejecta material). Av-11 and Av-12 exhibit almost the full range of mass wasting features observed on Vesta, such as slump blocks, spur-and-gully morphologies and landslides within craters. Processes of collapse, slope instability and seismically triggered events force material to slump down crater walls or scarps and produce landslides or rotational slump blocks. The spur-and-gully morphology that is known to form on Mars is also observed on Vesta; however, on Vesta this morphology formed under dry conditions.

A variety of geologic landforms and features are observed within quadrangle Av-13 Tuccia in the southern hemisphere of Vesta. The quadrangle covers parts of the highland Vestalia Terra as well as the floors of the large Rheasilvia and Veneneia impact basins, which results in a substantial elevation difference of more than 40 km between the northern and the southern portions of the quadrangle. Measurements of crater size–frequency distributions within and surrounding the Rheasilvia basin indicate that gravity-driven mass wasting in the interior of the basin has been important, and that the basin has a more ancient formation age than would be expected from the crater density on the basin floor alone. Subsequent to its formation, Rheasilvia was superimposed by several mid-sized impact craters. The most prominent craters are Tuccia, Eusebia, Vibidia, Galeria, and Antonia, whose geology and formation ages are investigated in detail in this work. These impact structures provide a variety of morphologies indicating different sorts of subsequent impact-related or gravity-driven mass wasting processes. Understanding the geologic history of the relatively young craters in the Rheasilvia basin is important in order to understand the even more degraded craters in other regions of Vesta.