ASU Scholarship Showcase

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: On 31 March 2013, the first human infections with the novel influenza A/H7N9 virus were reported in Eastern China. The outbreak expanded rapidly in geographic scope and size, with a total of 132 laboratory-confirmed cases reported by 3 June 2013, in 10 Chinese provinces and Taiwan. The incidence of A/H7N9 cases has stalled in recent weeks, presumably as a consequence of live bird market closures in the most heavily affected areas. Here we compare the transmission potential of influenza A/H7N9 with that of other emerging pathogens and evaluate the impact of intervention measures in an effort to guide pandemic preparedness.

Methods: We used a Bayesian approach combined with a SEIR (Susceptible-Exposed-Infectious-Removed) transmission model fitted to daily case data to assess the reproduction number (R) of A/H7N9 by province and to evaluate the impact of live bird market closures in April and May 2013. Simulation studies helped quantify the performance of our approach in the context of an emerging pathogen, where human-to-human transmission is limited and most cases arise from spillover events. We also used alternative approaches to estimate R based on individual-level information on prior exposure and compared the transmission potential of influenza A/H7N9 with that of other recent zoonoses.

Results: Estimates of R for the A/H7N9 outbreak were below the epidemic threshold required for sustained human-to-human transmission and remained near 0.1 throughout the study period, with broad 95% credible intervals by the Bayesian method (0.01 to 0.49). The Bayesian estimation approach was dominated by the prior distribution, however, due to relatively little information contained in the case data. We observe a statistically significant deceleration in growth rate after 6 April 2013, which is consistent with a reduction in A/H7N9 transmission associated with the preemptive closure of live bird markets. Although confidence intervals are broad, the estimated transmission potential of A/H7N9 appears lower than that of recent zoonotic threats, including avian influenza A/H5N1, swine influenza H3N2sw and Nipah virus.

Conclusion: Although uncertainty remains high in R estimates for H7N9 due to limited epidemiological information, all available evidence points to a low transmission potential. Continued monitoring of the transmission potential of A/H7N9 is critical in the coming months as intervention measures may be relaxed and seasonal factors could promote disease transmission in colder months.

Background: Immunosignaturing is a new peptide microarray based technology for profiling of humoral immune responses. Despite new challenges, immunosignaturing gives us the opportunity to explore new and fundamentally different research questions. In addition to classifying samples based on disease status, the complex patterns and latent factors underlying immunosignatures, which we attempt to model, may have a diverse range of applications.

Methods: We investigate the utility of a number of statistical methods to determine model performance and address challenges inherent in analyzing immunosignatures. Some of these methods include exploratory and confirmatory factor analyses, classical significance testing, structural equation and mixture modeling.

Results: We demonstrate an ability to classify samples based on disease status and show that immunosignaturing is a very promising technology for screening and presymptomatic screening of disease. In addition, we are able to model complex patterns and latent factors underlying immunosignatures. These latent factors may serve as biomarkers for disease and may play a key role in a bioinformatic method for antibody discovery.

Conclusion: Based on this research, we lay out an analytic framework illustrating how immunosignatures may be useful as a general method for screening and presymptomatic screening of disease as well as antibody discovery.

Background: The impact of socio-demographic factors and baseline health on the mortality burden of seasonal and pandemic influenza remains debated. Here we analyzed the spatial-temporal mortality patterns of the 1918 influenza pandemic in Spain, one of the countries of Europe that experienced the highest mortality burden.

Methods: We analyzed monthly death rates from respiratory diseases and all-causes across 49 provinces of Spain, including the Canary and Balearic Islands, during the period January-1915 to June-1919. We estimated the influenza-related excess death rates and risk of death relative to baseline mortality by pandemic wave and province. We then explored the association between pandemic excess mortality rates and health and socio-demographic factors, which included population size and age structure, population density, infant mortality rates, baseline death rates, and urbanization.

Results: Our analysis revealed high geographic heterogeneity in pandemic mortality impact. We identified 3 pandemic waves of varying timing and intensity covering the period from Jan-1918 to Jun-1919, with the highest pandemic-related excess mortality rates occurring during the months of October-November 1918 across all Spanish provinces. Cumulative excess mortality rates followed a south–north gradient after controlling for demographic factors, with the North experiencing highest excess mortality rates. A model that included latitude, population density, and the proportion of children living in provinces explained about 40% of the geographic variability in cumulative excess death rates during 1918–19, but different factors explained mortality variation in each wave.

Conclusions: A substantial fraction of the variability in excess mortality rates across Spanish provinces remained unexplained, which suggests that other unidentified factors such as comorbidities, climate and background immunity may have affected the 1918-19 pandemic mortality rates. Further archeo-epidemiological research should concentrate on identifying settings with combined availability of local historical mortality records and information on the prevalence of underlying risk factors, or patient-level clinical data, to further clarify the drivers of 1918 pandemic influenza mortality.

Background: Microarray image analysis processes scanned digital images of hybridized arrays to produce the input spot-level data for downstream analysis, so it can have a potentially large impact on those and subsequent analysis. Signal saturation is an optical effect that occurs when some pixel values for highly expressed genes or peptides exceed the upper detection threshold of the scanner software (2¹⁶ - 1 = 65, 535 for 16-bit images). In practice, spots with a sizable number of saturated pixels are often flagged and discarded. Alternatively, the saturated values are used without adjustments for estimating spot intensities. The resulting expression data tend to be biased downwards and can distort high-level analysis that relies on these data. Hence, it is crucial to effectively correct for signal saturation.

Results: We developed a flexible mixture model-based segmentation and spot intensity estimation procedure that accounts for saturated pixels by incorporating a censored component in the mixture model. As demonstrated with biological data and simulation, our method extends the dynamic range of expression data beyond the saturation threshold and is effective in correcting saturation-induced bias when the lost information is not tremendous. We further illustrate the impact of image processing on downstream classification, showing that the proposed method can increase diagnostic accuracy using data from a lymphoma cancer diagnosis study.

Conclusions: The presented method adjusts for signal saturation at the segmentation stage that identifies a pixel as part of the foreground, background or other. The cluster membership of a pixel can be altered versus treating saturated values as truly observed. Thus, the resulting spot intensity estimates may be more accurate than those obtained from existing methods that correct for saturation based on already segmented data. As a model-based segmentation method, our procedure is able to identify inner holes, fuzzy edges and blank spots that are common in microarray images. The approach is independent of microarray platform and applicable to both single- and dual-channel microarrays.

Background: High-throughput technologies such as DNA, RNA, protein, antibody and peptide microarrays are often used to examine differences across drug treatments, diseases, transgenic animals, and others. Typically one trains a classification system by gathering large amounts of probe-level data, selecting informative features, and classifies test samples using a small number of features. As new microarrays are invented, classification systems that worked well for other array types may not be ideal. Expression microarrays, arguably one of the most prevalent array types, have been used for years to help develop classification algorithms. Many biological assumptions are built into classifiers that were designed for these types of data. One of the more problematic is the assumption of independence, both at the probe level and again at the biological level. Probes for RNA transcripts are designed to bind single transcripts. At the biological level, many genes have dependencies across transcriptional pathways where co-regulation of transcriptional units may make many genes appear as being completely dependent. Thus, algorithms that perform well for gene expression data may not be suitable when other technologies with different binding characteristics exist. The immunosignaturing microarray is based on complex mixtures of antibodies binding to arrays of random sequence peptides. It relies on many-to-many binding of antibodies to the random sequence peptides. Each peptide can bind multiple antibodies and each antibody can bind multiple peptides. This technology has been shown to be highly reproducible and appears promising for diagnosing a variety of disease states. However, it is not clear what is the optimal classification algorithm for analyzing this new type of data.

Results: We characterized several classification algorithms to analyze immunosignaturing data. We selected several datasets that range from easy to difficult to classify, from simple monoclonal binding to complex binding patterns in asthma patients. We then classified the biological samples using 17 different classification algorithms. Using a wide variety of assessment criteria, we found ‘Naïve Bayes’ far more useful than other widely used methods due to its simplicity, robustness, speed and accuracy.

Conclusions: ‘Naïve Bayes’ algorithm appears to accommodate the complex patterns hidden within multilayered immunosignaturing microarray data due to its fundamental mathematical properties.

Background: Elucidating the role of the underlying risk factors for severe outcomes of the 2009 A/H1N1 influenza pandemic could be crucial to define priority risk groups in resource-limited settings in future pandemics.

Methods: We use individual-level clinical data on a large series of ARI (acute respiratory infection) hospitalizations from a prospective surveillance system of the Mexican Social Security medical system to analyze clinical features at presentation, admission delays, selected comorbidities and receipt of seasonal vaccine on the risk of A/H1N1-related death. We considered ARI hospitalizations and inpatient-deaths, and recorded demographic, geographic, and medical information on individual patients during August-December, 2009.

Results: Seasonal influenza vaccination was associated with a reduced risk of death among A/H1N1 inpatients (OR = 0.43 (95% CI: 0.25, 0.74)) after adjustment for age, gender, geography, antiviral treatment, admission delays, comorbidities and medical conditions. However, this result should be interpreted with caution as it could have been affected by factors not directly measured in our study. Moreover, the effect of antiviral treatment against A/H1N1 inpatient death did not reach statistical significance (OR = 0.56 (95% CI: 0.29, 1.10)) probably because only 8.9% of A/H1N1 inpatients received antiviral treatment. Moreover, diabetes (OR = 1.6) and immune suppression (OR = 2.3) were statistically significant risk factors for death whereas asthmatic persons (OR = 0.3) or pregnant women (OR = 0.4) experienced a reduced fatality rate among A/H1N1 inpatients. We also observed an increased risk of death among A/H1N1 inpatients with admission delays >2 days after symptom onset (OR = 2.7). Similar associations were also observed for A/H1N1-negative inpatients.

Conclusions: Geographical variation in identified medical risk factors including prevalence of diabetes and immune suppression may in part explain between-country differences in pandemic mortality burden. Furthermore, access to care including hospitalization without delay and antiviral treatment and are also important factors, as well as vaccination coverage with the 2008–09 trivalent inactivated influenza vaccine.

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: 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.

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.