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.

S-cysteinylated albumin and methionine-oxidized apolipoprotein A-I (apoA-I) have been posed as candidate markers of diseases associated with oxidative stress. Here, a dilute-and-shoot form of LC–electrospray ionization–MS requiring half a microliter of blood plasma was employed to simultaneously quantify the relative abundance of these oxidized proteoforms in samples stored at −80 °C, −20 °C, and room temperature and exposed to multiple freeze-thaw cycles and other adverse conditions in order to assess the possibility that protein oxidation may occur as a result of poor sample storage or handling. Samples from a healthy donor and a participant with poorly controlled type 2 diabetes started at the same low level of protein oxidation and behaved similarly; significant increases in albumin oxidation via S-cysteinylation were found to occur within hours at room temperature and days at −20 °C. Methionine oxidation of apoA-I took place on a longer time scale, setting in after albumin oxidation reached a plateau. Freeze–thaw cycles had a minimal effect on protein oxidation. In matched collections, protein oxidation in serum was the same as that in plasma. Albumin and apoA-I oxidation were not affected by sample headspace or the degree to which vials were sealed. ApoA-I, however, was unexpectedly found to oxidize faster in samples with lower surface-area-to-volume ratios. An initial survey of samples from patients with inflammatory conditions normally associated with elevated oxidative stress-including acute myocardial infarction and prostate cancer—demonstrated a lack of detectable apoA-I oxidation. Albumin S-cysteinylation in these samples was consistent with known but relatively brief exposures to temperatures above −30 °C (the freezing point of blood plasma). Given their properties and ease of analysis, these oxidized proteoforms, once fully validated, may represent the first markers of blood plasma specimen integrity based on direct measurement of oxidative molecular damage that can occur under suboptimal storage conditions.

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.

Human societies are unique in the level of cooperation among non-kin. Evolutionary models explaining this behavior typically assume pure strategies of cooperation and defection. Behavioral experiments, however, demonstrate that humans are typically conditional co-operators who have other-regarding preferences. Building on existing models on the evolution of cooperation and costly punishment, we use a utilitarian formulation of agent decision making to explore conditions that support the emergence of cooperative behavior. Our results indicate that cooperation levels are significantly lower for larger groups in contrast to the original pure strategy model. Here, defection behavior not only diminishes the public good, but also affects the expectations of group members leading conditional co-operators to change their strategies. Hence defection has a more damaging effect when decisions are based on expectations and not only pure strategies.

Collective behaviors in social insect societies often emerge from simple local rules. However, little is known about how these behaviors are dynamically regulated in response to environmental changes. Here, we use a compartmental modeling approach to identify factors that allow harvester ant colonies to regulate collective foraging activity in response to their environment. We propose a set of differential equations describing the dynamics of: (1) available foragers inside the nest, (2) active foragers outside the nest, and (3) successful returning foragers, to understand how colony-specific parameters, such as baseline number of foragers, interactions among foragers, food discovery rates, successful forager return rates, and foraging duration might influence collective foraging dynamics, while maintaining functional robustness to perturbations. Our analysis indicates that the model can undergo a forward (transcritical) bifurcation or a backward bifurcation depending on colony-specific parameters. In the former case, foraging activity persists when the average number of recruits per successful returning forager is larger than one. In the latter case, the backward bifurcation creates a region of bistability in which the size and fate of foraging activity depends on the distribution of the foraging workforce among the model׳s compartments. We validate the model with experimental data from harvester ants (Pogonomyrmex barbatus) and perform sensitivity analysis. Our model provides insights on how simple, local interactions can achieve an emergent and robust regulatory system of collective foraging activity in ant colonies.

Most studies on the response of socioeconomic systems to a sudden shift focus on long-term equilibria or end points. Such narrow focus forgoes many valuable insights. Here we examine the transient dynamics of regime shift on a divided population, exemplified by societies divided ideologically, politically, economically, or technologically. Replicator dynamics is used to investigate the complex transient dynamics of the population response. Though simple, our modeling approach exhibits a surprisingly rich and diverse array of dynamics. Our results highlight the critical roles played by diversity in strategies and the magnitude of the shift. Importantly, it allows for a variety of strategies to arise organically as an integral part of the transient dynamics-as opposed to an independent process-of population response to a regime shift, providing a link between the population's past and future diversity patterns. Several combinations of different populations' strategy distributions and shifts were systematically investigated. Such rich dynamics highlight the challenges of anticipating the response of a divided population to a change. The findings in this paper can potentially improve our understanding of a wide range of socio-ecological and technological transitions.

A major conundrum in evolution is that, despite natural selection, polymorphism is still omnipresent in nature: Numerous species exhibit multiple morphs, namely several abundant values of an important trait. Polymorphism is particularly prevalent in asymmetric traits, which are beneficial to their carrier in disruptive competitive interference but at the same time bear disadvantages in other aspects, such as greater mortality or lower fecundity. Here we focus on asymmetric traits in which a better competitor disperses fewer offspring in the absence of competition. We report a general pattern in which polymorphic populations emerge when disruptive selection increases: The stronger the selection, the greater the number of morphs that evolve. This pattern is general and is insensitive to the form of the fitness function. The pattern is somewhat counterintuitive since directional selection is excepted to sharpen the trait distribution and thereby reduce its diversity (but note that similar patterns were suggested in studies that demonstrated increased biodiversity as local selection increases in ecological communities). We explain the underlying mechanism in which stronger selection drives the population towards more competitive values of the trait, which in turn reduces the population density, thereby enabling lesser competitors to stably persist with reduced need to directly compete. Thus, we believe that the pattern is more general and may apply to asymmetric traits more broadly. This robust pattern suggests a comparative, unified explanation to a variety of polymorphic traits in nature.

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.

Background: While prior studies have quantified the mortality burden of the 1957 H2N2 influenza pandemic at broad geographic regions in the United States, little is known about the pandemic impact at a local level. Here we focus on analyzing the transmissibility and mortality burden of this pandemic in Arizona, a setting where the dry climate was promoted as reducing respiratory illness transmission yet tuberculosis prevalence was high.

Methods: Using archival death certificates from 1954 to 1961, we quantified the age-specific seasonal patterns, excess-mortality rates, and transmissibility patterns of the 1957 H2N2 pandemic in Maricopa County, Arizona. By applying cyclical Serfling linear regression models to weekly mortality rates, the excess-mortality rates due to respiratory and all-causes were estimated for each age group during the pandemic period. The reproduction number was quantified from weekly data using a simple growth rate method and assumed generation intervals of 3 and 4 days. Local newspaper articles published during 1957–1958 were also examined.

Results: Excess-mortality rates varied between waves, age groups, and causes of death, but overall remained low. From October 1959-June 1960, the most severe wave of the pandemic, the absolute excess-mortality rate based on respiratory deaths per 10,000 population was 16.59 in the elderly (≥65 years). All other age groups exhibit very low excess-mortality and the typical U-shaped age-pattern was absent. However, the standardized mortality ratio was greatest (4.06) among children and young adolescents (5–14 years) from October 1957-March 1958, based on mortality rates of respiratory deaths. Transmissibility was greatest during the same 1957–1958 period, when the mean reproduction number was estimated at 1.08–1.11, assuming 3- or 4-day generation intervals with exponential or fixed distributions.

Conclusions: Maricopa County exhibited very low mortality impact associated with the 1957 influenza pandemic. Understanding the relatively low excess-mortality rates and transmissibility in Maricopa County during this historic pandemic may help public health officials prepare for and mitigate future outbreaks of influenza.

Background: Influenza viruses are a major cause of morbidity and mortality worldwide. Vaccination remains a powerful tool for preventing or mitigating influenza outbreaks. Yet, vaccine supplies and daily administration capacities are limited, even in developed countries. Understanding how such constraints can alter the mitigating effects of vaccination is a crucial part of influenza preparedness plans. Mathematical models provide tools for government and medical officials to assess the impact of different vaccination strategies and plan accordingly. However, many existing models of vaccination employ several questionable assumptions, including a rate of vaccination proportional to the population at each point in time.

Methods: We present a SIR-like model that explicitly takes into account vaccine supply and the number of vaccines administered per day and places data-informed limits on these parameters. We refer to this as the non-proportional model of vaccination and compare it to the proportional scheme typically found in the literature.

Results: The proportional and non-proportional models behave similarly for a few different vaccination scenarios. However, there are parameter regimes involving the vaccination campaign duration and daily supply limit for which the non-proportional model predicts smaller epidemics that peak later, but may last longer, than those of the proportional model. We also use the non-proportional model to predict the mitigating effects of variably timed vaccination campaigns for different levels of vaccination coverage, using specific constraints on daily administration capacity.

Conclusions: The non-proportional model of vaccination is a theoretical improvement that provides more accurate predictions of the mitigating effects of vaccination on influenza outbreaks than the proportional model. In addition, parameters such as vaccine supply and daily administration limit can be easily adjusted to simulate conditions in developed and developing nations with a wide variety of financial and medical resources. Finally, the model can be used by government and medical officials to create customized pandemic preparedness plans based on the supply and administration constraints of specific communities.