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.

Biological Soil Crusts (BSCs) are organosedimentary assemblages comprised of microbes and minerals in topsoil of terrestrial environments. BSCs strongly impact soil quality in dryland ecosystems (e.g., soil structure and nutrient yields) due to pioneer species such as Microcoleus vaginatus; phototrophs that produce filaments that bind the soil together, and support an array of heterotrophic microorganisms. These microorganisms in turn contribute to soil stability and biogeochemistry of BSCs. Non-cyanobacterial populations of BSCs are less well known than cyanobacterial populations. Therefore, we attempted to isolate a broad range of numerically significant and phylogenetically representative BSC aerobic heterotrophs. Combining simple pre-treatments (hydration of BSCs under dark and light) and isolation strategies (media with varying nutrient availability and protection from oxidative stress) we recovered 402 bacterial and one fungal isolate in axenic culture, which comprised 116 phylotypes (at 97% 16S rRNA gene sequence homology), 115 bacterial and one fungal. Each medium enriched a mostly distinct subset of phylotypes, and cultivated phylotypes varied due to the BSC pre-treatment. The fraction of the total phylotype diversity isolated, weighted by relative abundance in the community, was determined by the overlap between isolate sequences and OTUs reconstructed from metagenome or metatranscriptome reads. Together, more than 8% of relative abundance of OTUs in the metagenome was represented by our isolates, a cultivation efficiency much larger than typically expected from most soils. We conclude that simple cultivation procedures combined with specific pre-treatment of samples afford a significant reduction in the culturability gap, enabling physiological and metabolic assays that rely on ecologically relevant axenic cultures.

Background: Increasing our understanding of the factors affecting the severity of the 2009 A/H1N1 influenza pandemic in different regions of the world could lead to improved clinical practice and mitigation strategies for future influenza pandemics. Even though a number of studies have shed light into the risk factors associated with severe outcomes of 2009 A/H1N1 influenza infections in different populations (e.g., [1-5]), analyses of the determinants of mortality risk spanning multiple pandemic waves and geographic regions are scarce. Between-country differences in the mortality burden of the 2009 pandemic could be linked to differences in influenza case management, underlying population health, or intrinsic differences in disease transmission [6]. Additional studies elucidating the determinants of disease severity globally are warranted to guide prevention efforts in future influenza pandemics.

In Mexico, the 2009 A/H1N1 influenza pandemic was characterized by a three-wave pattern occurring in the spring, summer, and fall of 2009 with substantial geographical heterogeneity [7]. A recent study suggests that Mexico experienced high excess mortality burden during the 2009 A/H1N1 influenza pandemic relative to other countries [6]. However, an assessment of potential factors that contributed to the relatively high pandemic death toll in Mexico are lacking. Here, we fill this gap by analyzing a large series of laboratory-confirmed A/H1N1 influenza cases, hospitalizations, and deaths monitored by the Mexican Social Security medical system during April 1 through December 31, 2009 in Mexico. In particular, we quantify the association between disease severity, hospital admission delays, and neuraminidase inhibitor use by demographic characteristics, pandemic wave, and geographic regions of Mexico.

Methods: We analyzed a large series of laboratory-confirmed pandemic A/H1N1 influenza cases from a prospective surveillance system maintained by the Mexican Social Security system, April-December 2009. We considered a spectrum of disease severity encompassing outpatient visits, hospitalizations, and deaths, and recorded demographic and geographic information on individual patients. We assessed the impact of neuraminidase inhibitor treatment and hospital admission delay (≤ > 2 days after disease onset) on the risk of death by multivariate logistic regression.

Results: Approximately 50% of all A/H1N1-positive patients received antiviral medication during the Spring and Summer 2009 pandemic waves in Mexico while only 9% of A/H1N1 cases received antiviral medications during the fall wave (P < 0.0001). After adjustment for age, gender, and geography, antiviral treatment significantly reduced the risk of death (OR = 0.52 (95% CI: 0.30, 0.90)) while longer hospital admission delays increased the risk of death by 2.8-fold (95% CI: 2.25, 3.41).

Conclusions: Our findings underscore the potential impact of decreasing admission delays and increasing antiviral use to mitigate the mortality burden of future influenza pandemics.

The heterocyclic indole-alkaloid scytonemin is a sunscreen found exclusively among cyanobacteria. An 18-gene cluster is responsible for scytonemin production in Nostoc punctiforme ATCC 29133. The upstream genes scyABCDEF in the cluster are proposed to be responsible for scytonemin biosynthesis from aromatic amino acid substrates. In vitro studies of ScyA, ScyB, and ScyC proved that these enzymes indeed catalyze initial pathway reactions. Here we characterize the role of ScyD, ScyE, and ScyF, which were logically predicted to be responsible for late biosynthetic steps, in the biological context of N. punctiforme. In-frame deletion mutants of each were constructed (ΔscyD, ΔscyE, and ΔscyF) and their phenotypes studied. Expectedly, ΔscyE presents a scytoneminless phenotype, but no accumulation of the predicted intermediaries. Surprisingly, ΔscyD retains scytonemin production, implying that it is not required for biosynthesis. Indeed, scyD presents an interesting evolutionary paradox: it likely originated in a duplication event from scyE, and unlike other genes in the operon, it has not been subjected to purifying selection. This would suggest that it is a pseudogene, and yet scyD is highly conserved in the scytonemin operon of cyanobacteria. ΔscyF also retains scytonemin production, albeit exhibiting a reduction of the production yield compared with the wild-type. This indicates that ScyF is not essential but may play an adjuvant role for scytonemin synthesis. Altogether, our findings suggest that these downstream genes are not responsible, as expected, for the late steps of scytonemin synthesis and we must look for those functions elsewhere. These findings are particularly important for biotechnological production of this sunscreen through heterologous expression of its genes in more tractable organisms.

The 1918 influenza pandemic was a major epidemiological event of the twentieth century resulting in at least twenty million deaths worldwide; however, despite its historical, epidemiological, and biological relevance, it remains poorly understood. Here we examine the relationship between annual pneumonia and influenza death rates in the pre-pandemic (1910–17) and pandemic (1918–20) periods and the scaling of mortality with latitude, longitude and population size, using data from 66 large cities of the United States. The mean pre-pandemic pneumonia death rates were highly associated with pneumonia death rates during the pandemic period (Spearman ρ = 0.64–0.72; P<0.001). By contrast, there was a weak correlation between pre-pandemic and pandemic influenza mortality rates. Pneumonia mortality rates partially explained influenza mortality rates in 1918 (ρ = 0.34, P = 0.005) but not during any other year. Pneumonia death counts followed a linear relationship with population size in all study years, suggesting that pneumonia death rates were homogeneous across the range of population sizes studied. By contrast, influenza death counts followed a power law relationship with a scaling exponent of ∼0.81 (95%CI: 0.71, 0.91) in 1918, suggesting that smaller cities experienced worst outcomes during the pandemic. A linear relationship was observed for all other years. Our study suggests that mortality associated with the 1918–20 influenza pandemic was in part predetermined by pre-pandemic pneumonia death rates in 66 large US cities, perhaps through the impact of the physical and social structure of each city. Smaller cities suffered a disproportionately high per capita influenza mortality burden than larger ones in 1918, while city size did not affect pneumonia mortality rates in the pre-pandemic and pandemic periods.

Background: The historical Japanese influenza vaccination program targeted at schoolchildren provides a unique opportunity to evaluate the indirect benefits of vaccinating high-transmitter groups to mitigate disease burden among seniors. Here we characterize the indirect mortality benefits of vaccinating schoolchildren based on data from Japan and the US.

Methods: We compared age-specific influenza-related excess mortality rates in Japanese seniors aged ≥65 years during the schoolchildren vaccination program (1978–1994) and after the program was discontinued (1995–2006). Indirect vaccine benefits were adjusted for demographic changes, socioeconomics and dominant influenza subtype; US mortality data were used as a control.

Results: We estimate that the schoolchildren vaccination program conferred a 36% adjusted mortality reduction among Japanese seniors (95%CI: 17–51%), corresponding to ∼1,000 senior deaths averted by vaccination annually (95%CI: 400–1,800). In contrast, influenza-related mortality did not change among US seniors, despite increasing vaccine coverage in this population.

Conclusions: The Japanese schoolchildren vaccination program was associated with substantial indirect mortality benefits in seniors.

Background: The role of demographic factors, climatic conditions, school cycles, and connectivity patterns in shaping the spatio-temporal dynamics of pandemic influenza is not clearly understood. Here we analyzed the spatial, age and temporal evolution of the 2009 A/H1N1 influenza pandemic in Chile, a southern hemisphere country covering a long and narrow strip comprising latitudes 17°S to 56°S.

Methods: We analyzed the dissemination patterns of the 2009 A/H1N1 pandemic across 15 regions of Chile based on daily hospitalizations for severe acute respiratory disease and laboratory confirmed A/H1N1 influenza infection from 01-May to 31-December, 2009. We explored the association between timing of pandemic onset and peak pandemic activity and several geographical and demographic indicators, school vacations, climatic factors, and international passengers. We also estimated the reproduction number (R) based on the growth rate of the exponential pandemic phase by date of symptoms onset, estimated using maximum likelihood methods.

Results: While earlier pandemic onset was associated with larger population size, there was no association with connectivity, demographic, school or climatic factors. In contrast, there was a latitudinal gradient in peak pandemic timing, representing a 16-39-day lag in disease activity from the southern regions relative to the northernmost region (P < 0.001). Geographical differences in latitude of Chilean regions, maximum temperature and specific humidity explained 68.5% of the variability in peak timing (P = 0.01). In addition, there was a decreasing gradient in reproduction number from south to north Chile (P < 0.0001). The regional mean R estimates were 1.6-2.0, 1.3-1.5, and 1.2-1.3 for southern, central and northern regions, respectively, which were not affected by the winter vacation period.

Conclusions: There was a lag in the period of most intense 2009 pandemic influenza activity following a South to North traveling pattern across regions of Chile, significantly associated with geographical differences in minimum temperature and specific humidity. The latitudinal gradient in timing of pandemic activity was accompanied by a gradient in reproduction number (P < 0.0001). Intensified surveillance strategies in colder and drier southern regions could lead to earlier detection of pandemic influenza viruses and improved control outcomes.

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.

N₂ fixation and ammonia oxidation (AO) are the two most important processes in the nitrogen (N) cycle of biological soil crusts (BSCs). We studied the short-term response of acetylene reduction assay (ARA) rates, an indicator of potential N₂ fixation, and AO rates to temperature (T, -5°C to 35°C) in BSC of different successional stages along the BSC ecological succession and geographic origin (hot Chihuahuan and cooler Great Basin deserts). ARA in all BSCs increased with T until saturation occurred between 15 and 20°C, and declined at 30–35°C. Culture studies using cyanobacteria isolated from these crusts indicated that the saturating effect was traceable to their inability to grow well diazotrophically within the high temperature range. Below saturation, temperature response was exponential, with Q₁₀ significantly different in the two areas (~ 5 for Great Basin BSCs; 2–3 for Chihuahuan BSCs), but similar between the two successional stages. However, in contrast to ARA, AO showed a steady increase to 30–35°C in Great Basin, and Chihuhuan BSCs showed no inhibition at any tested temperature. The T response of AO also differed significantly between Great Basin (Q₁₀ of 4.5–4.8) and Chihuahuan (Q₁₀ of 2.4–2.6) BSCs, but not between successional stages. Response of ARA rates to T did not differ from that of AO in either desert. Thus, while both processes scaled to T in unison until 20°C, they separated to an increasing degree at higher temperature. As future warming is likely to occur in the regions where BSCs are often the dominant living cover, this predicted decoupling is expected to result in higher proportion of nitrates in soil relative to ammonium. As nitrate is more easily lost as leachate or to be reduced to gaseous forms, this could mean a depletion of soil N over large landscapes globally.