ASU Scholarship Showcase

Methicillin resistant Staphylococcus aureus (MRSA) is currently a major cause of skin and soft tissue infections (SSTI) in the United States. Seasonal variation of MRSA infections in hospital settings has been widely observed. However, systematic time-series analysis of incidence data is desirable to understand the seasonality of community acquired (CA)-MRSA infections at the population level. In this paper, using data on monthly SSTI incidence in children aged 0–19 years and enrolled in Medicaid in Maricopa County, Arizona, from January 2005 to December 2008, we carried out time-series and nonlinear regression analysis to determine the periodicity, trend, and peak timing in SSTI incidence in children at different age: 0-4 years, 5-9 years, 10-14 years, and 15-19 years. We also assessed the temporal correlation between SSTI incidence and meteorological variables including average temperature and humidity. Our analysis revealed a strong annual seasonal pattern of SSTI incidence with peak occurring in early September. This pattern was consistent across age groups. Moreover, SSTIs followed a significantly increasing trend over the 4-year study period with annual incidence increasing from 3.36% to 5.55% in our pediatric population of approximately 290,000. We also found a significant correlation between the temporal variation in SSTI incidence and mean temperature and specific humidity. Our findings could have potential implications on prevention and control efforts against CA-MRSA.

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: Dengue fever is a mosquito-borne disease that affects between 50 and 100 million people each year. Increasing our understanding of the heterogeneous transmission patterns of dengue at different spatial scales could have considerable public health value by guiding intervention strategies.

Methods: Based on the weekly number of dengue cases in Perú by province, we investigated the association between dengue incidence during the period 1994-2008 and demographic and climate factors across geographic regions of the country.

Results: Our findings support the presence of significant differences in the timing of dengue epidemics between jungle and coastal regions, with differences significantly associated with the timing of the seasonal cycle of mean temperature.

Conclusions: Dengue is highly persistent in jungle areas of Perú where epidemics peak most frequently around March when rainfall is abundant. Differences in the timing of dengue epidemics in jungle and coastal regions are significantly associated with the seasonal temperature cycle. Our results suggest that dengue is frequently imported into coastal regions through infective sparks from endemic jungle areas and/or cities of other neighboring endemic countries, where propitious environmental conditions promote year-round mosquito breeding sites. If jungle endemic areas are responsible for multiple dengue introductions into coastal areas, our findings suggest that curtailing the transmission of dengue in these most persistent areas could lead to significant reductions in dengue incidence in coastal areas where dengue incidence typically reaches low levels during the dry season.

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.

Most previous works on complete synchronization of chaotic oscillators focused on the one-channel interaction scheme where the oscillators are coupled through only one variable or a symmetric set of variables. Using the standard framework of master-stability function (MSF), we investigate the emergence of complex synchronization behaviors under all possible configurations of two-channel coupling, which include, for example, all possible cross coupling schemes among the dynamical variables. Utilizing the classic Rössler and Lorenz oscillators, we find a rich variety of synchronization phenomena not present in any previously extensively studied, single-channel coupling configurations. For example, in many cases two coupling channels can enhance or even generate synchronization where there is only weak or no synchronization under only one coupling channel, which has been verified in a coupled neuron system. There are also cases where the oscillators are originally synchronized under one coupling channel, but an additional synchronizable coupling channel can, however, destroy synchronization. Direct numerical simulations of actual synchronization dynamics verify the MSF-based predictions. Our extensive computation and heuristic analysis provide an atlas for synchronization of chaotic oscillators coupled through two channels, which can be used as a systematic reference to facilitate further research in this area.

Nonhyperbolicity, as characterized by the coexistence of Kolmogorov-Arnold-Moser (KAM) tori and chaos in the phase space, is generic in classical Hamiltonian systems. An open but fundamental question in physics concerns the relativistic quantum manifestations of nonhyperbolic dynamics. We choose the mushroom billiard that has been mathematically proven to be nonhyperbolic, and study the resonant tunneling dynamics of a massless Dirac fermion. We find that the tunneling rate as a function of the energy exhibits a striking "clustering" phenomenon, where the majority of the values of the rate concentrate on a narrow region, as a result of the chaos component in the classical phase space. Relatively few values of the tunneling rate, however, spread outside the clustering region due to the integrable component. Resonant tunneling of electrons in nonhyperbolic chaotic graphene systems exhibits a similar behavior. To understand these numerical results, we develop a theoretical framework by combining analytic solutions of the Dirac equation in certain integrable domains and physical intuitions gained from current understanding of the quantum manifestations of chaos. In particular, we employ a theoretical formalism based on the concept of self-energies to calculate the tunneling rate and analytically solve the Dirac equation in one dimension as well as in two dimensions for a circular-ring-type of tunneling systems exhibiting integrable dynamics in the classical limit. Because relatively few and distinct classical periodic orbits are present in the integrable component, the corresponding relativistic quantum states can have drastically different behaviors, leading to a wide spread in the values of the tunneling rate in the energy-rate plane. In contrast, the chaotic component has embedded within itself an infinite number of unstable periodic orbits, which provide far more quantum states for tunneling. Due to the nature of chaos, these states are characteristically similar, leading to clustering of the values of the tunneling rate in a narrow band. The appealing characteristic of our work is a demonstration and physical understanding of the "mixed" role played by chaos and regular dynamics in shaping relativistic quantum tunneling dynamics.

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.

Background: Ebola virus disease (EVD) has generated a large epidemic in West Africa since December 2013. This mini-review is aimed to clarify and illustrate different theoretical concepts of infectiousness in order to compare the infectiousness across different communicable diseases including EVD.

Methods: We employed a transmission model that rests on the renewal process in order to clarify theoretical concepts on infectiousness, namely the basic reproduction number, R₀, which measures the infectiousness per generation of cases, the force of infection (i.e. the hazard rate of infection), the intrinsic growth rate (i.e. infectiousness per unit time) and the per-contact probability of infection (i.e. infectiousness per effective contact).

Results: Whereas R₀ of EVD is similar to that of influenza, the growth rate (i.e. the measure of infectiousness per unit time) for EVD was shown to be comparatively lower than that for influenza. Moreover, EVD and influenza differ in mode of transmission whereby the probability of transmission per contact is lower for EVD compared to that of influenza.

Conclusions: The slow spread of EVD associated with the need for physical contact with body fluids supports social distancing measures including contact tracing and case isolation. Descriptions and interpretations of different variables quantifying infectiousness need to be used clearly and objectively in the scientific community and for risk communication.

Carefully calibrated transmission models have the potential to guide public health officials on the nature and scale of the interventions required to control epidemics. In the context of the ongoing Ebola virus disease (EVD) epidemic in Liberia, Drake and colleagues, in this issue of PLOS Biology, employed an elegant modeling approach to capture the distributions of the number of secondary cases that arise in the community and health care settings in the context of changing population behaviors and increasing hospital capacity. Their findings underscore the role of increasing the rate of safe burials and the fractions of infectious individuals who seek hospitalization together with hospital capacity to achieve epidemic control. However, further modeling efforts of EVD transmission and control in West Africa should utilize the spatial-temporal patterns of spread in the region by incorporating spatial heterogeneity in the transmission process. Detailed datasets are urgently needed to characterize temporal changes in population behaviors, contact networks at different spatial scales, population mobility patterns, adherence to infection control measures in hospital settings, and hospitalization and reporting rates.

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.