ASU Scholarship Showcase

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 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: Advancements in geographic information systems over the past two decades have increased the specificity by which an individual’s neighborhood environment may be spatially defined for physical activity and health research. This study investigated how different types of street network buffering methods compared in measuring a set of commonly used built environment measures (BEMs) and tested their performance on associations with physical activity outcomes.

Methods: An internationally-developed set of objective BEMs using three different spatial buffering techniques were used to evaluate the relative differences in resulting explanatory power on self-reported physical activity outcomes. BEMs were developed in five countries using ‘sausage,’ ‘detailed-trimmed,’ and ‘detailed,’ network buffers at a distance of 1 km around participant household addresses (n = 5883).

Results: BEM values were significantly different (p < 0.05) for 96% of sausage versus detailed-trimmed buffer comparisons and 89% of sausage versus detailed network buffer comparisons. Results showed that BEM coefficients in physical activity models did not differ significantly across buffering methods, and in most cases BEM associations with physical activity outcomes had the same level of statistical significance across buffer types. However, BEM coefficients differed in significance for 9% of the sausage versus detailed models, which may warrant further investigation.

Conclusions: Results of this study inform the selection of spatial buffering methods to estimate physical activity outcomes using an internationally consistent set of BEMs. Using three different network-based buffering methods, the findings indicate significant variation among BEM values, however associations with physical activity outcomes were similar across each buffering technique. The study advances knowledge by presenting consistently assessed relationships between three different network buffer types and utilitarian travel, sedentary behavior, and leisure-oriented physical activity outcomes.

Background: The World Health Organization recommends strategies to improve urban design, public transportation, and recreation facilities to facilitate physical activity for non-communicable disease prevention for an increasingly urbanized global population. Most evidence supporting environmental associations with physical activity comes from single countries or regions with limited variation in urban form. This paper documents variation in comparable built environment features across countries from diverse regions.

Methods: The International Physical Activity and the Environment Network (IPEN) study of adults aimed to measure the full range of variation in the built environment using geographic information systems (GIS) across 12 countries on 5 continents. Investigators in Australia, Belgium, Brazil, Colombia, the Czech Republic, Denmark, China, Mexico, New Zealand, Spain, the United Kingdom, and the United States followed a common research protocol to develop internationally comparable measures. Using detailed instructions, GIS-based measures included features such as walkability (i.e., residential density, street connectivity, mix of land uses), and access to public transit, parks, and private recreation facilities around each participant’s residential address using 1-km and 500-m street network buffers.

Results: Eleven of 12 countries and 15 cities had objective GIS data on built environment features. We observed a 38-fold difference in median residential densities, a 5-fold difference in median intersection densities and an 18-fold difference in median park densities. Hong Kong had the highest and North Shore, New Zealand had the lowest median walkability index values, representing a difference of 9 standard deviations in GIS-measured walkability.

Conclusions: Results show that comparable measures can be created across a range of cultural settings revealing profound global differences in urban form relevant to physical activity. These measures allow cities to be ranked more precisely than previously possible. The highly variable measures of urban form will be used to explain individuals’ physical activity, sedentary behaviors, body mass index, and other health outcomes on an international basis. Present measures provide the ability to estimate dose–response relationships from projected changes to the built environment that would otherwise be impossible.

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: 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: To more accurately quantify the potential impact of the neighbourhood environment on adults’ physical activity (PA), it is important to compare environment-PA associations between periods of the day or week when adults are more versus less likely to be in their neighbourhood and utilise its PA resources. We examined whether, among adults from 10 countries, associations between objectively-assessed neighbourhood environment attributes and moderate-to-vigorous physical activity (MVPA) varied by time of the day and day of the week. The secondary aim was to examine whether such associations varied by employment status, gender and city.

Methods: This cross-sectional study included 6,712 adults from 14 cities across 10 countries with ≥1 day of valid accelerometer-assessed MVPA and complete information on socio-demographic and objectively-assessed environmental characteristics within 0.5 and 1 km street-network buffers around the home. Accelerometer measures (MVPA min/h) were created for six time periods from early morning until late evening/night, for weekdays and weekend days separately. Associations were estimated using generalized additive mixed models.

Results: Time of the day, day of week, gender and employment status were significant moderators of environment-MVPA associations. Land use mix was positively associated with MVPA in women who were employed and in men irrespective of their employment status. The positive associations between MVPA and net residential density, intersection density and land use mix were stronger in the mornings of weekdays and the afternoon/evening periods of both weekdays and weekend days. Associations between number of parks and MVPA were stronger in the mornings and afternoon/evenings irrespective of day of the week. Public transport density showed consistent positive associations with MVPA during weekends, while stronger effects on weekdays were observed in the morning and early evenings.

Conclusions: This study suggests that space and time constraints in adults’ daily activities are important factors that determine the impact of neighbourhood attributes on PA. Consideration of time-specific associations is important to better characterise the magnitude of the effects of the neighbourhood environment on PA. Future research will need to examine the contribution of built environment characteristics of areas surrounding other types of daily life centres (e.g., workplaces) to explaining adults’ PA at specific times of the day.

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.