The role of vaccination coverage, individual behaviors, and the public health response in the control of measles epidemics: an agent-based simulation for California.

BACKGROUND: Measles cases continue to occur among susceptible individuals despite the elimination of endemic measles transmission in the United States. Clustering of disease susceptibility can threaten herd immunity and impact the likelihood of disease outbreaks in a highly vaccinated population. Previous studies have examined the role of contact tracing to control infectious diseases among clustered populations, but have not explicitly modeled the public health response using an agent-based model. METHODS: We developed an agent-based simulation model of measles transmission using the Framework for Reconstructing Epidemiological Dynamics (FRED) and the Synthetic Population Database maintained by RTI International. The simulation of measles transmission was based on interactions among individuals in different places: households, schools, daycares, workplaces, and neighborhoods. The model simulated different levels of immunity clustering, vaccination coverage, and contact investigations with delays caused by individuals' behaviors and/or the delay in a health department's response. We examined the effects of these characteristics on the probability of uncontrolled measles outbreaks and the outbreak size in 365 days after the introduction of one index case into a synthetic population. RESULTS: We found that large measles outbreaks can be prevented with contact investigations and moderate contact rates by having (1) a very high vaccination coverage (≥ 95%) with a moderate to low level of immunity clustering (≤ 0.5) for individuals aged less than or equal to 18 years, or (2) a moderate vaccination coverage (85% or 90%) with no immunity clustering for individuals (≤ 18 years of age), a short intervention delay, and a high probability that a contact can be traced. Without contact investigations, measles outbreaks may be prevented by the highest vaccination coverage with no immunity clustering for individuals (≤ 18 years of age) with moderate contact rates; but for the highest contact rates, even the highest coverage with no immunity clustering for individuals (≤ 18 years of age) cannot completely prevent measles outbreaks. CONCLUSIONS: The simulation results demonstrated the importance of vaccination coverage, clustering of immunity, and contact investigations in preventing uncontrolled measles outbreaks.

Quantifying the economic value and quality of life impact of earlier influenza vaccination.

BACKGROUND: Influenza vaccination is administered throughout the influenza disease season, even as late as March. Given such timing, what is the value of vaccinating the population earlier than currently being practiced? METHODS: We used real data on when individuals were vaccinated in Allegheny County, Pennsylvania, and the following 2 models to determine the value of vaccinating individuals earlier (by the end of September, October, and November): Framework for Reconstructing Epidemiological Dynamics (FRED), an agent-based model (ABM), and FluEcon, our influenza economic model that translates cases from the ABM to outcomes and costs [health care and lost productivity costs and quality-adjusted life-years (QALYs)]. We varied the reproductive number (R0) from 1.2 to 1.6. RESULTS: Applying the current timing of vaccinations averted 223,761 influenza cases, $16.3 million in direct health care costs, $50.0 million in productivity losses, and 804 in QALYs, compared with no vaccination (February peak, R0 1.2). When the population does not have preexisting immunity and the influenza season peaks in February (R0 1.2-1.6), moving individuals who currently received the vaccine after September to the end of September could avert an additional 9634-17,794 influenza cases, $0.6-$1.4 million in direct costs, $2.1-$4.0 million in productivity losses, and 35-64 QALYs. Moving the vaccination of just children to September (R0 1.2-1.6) averted 11,366-1660 influenza cases, $0.6-$0.03 million in direct costs, $2.3-$0.2 million in productivity losses, and 42-8 QALYs. Moving the season peak to December increased these benefits, whereas increasing preexisting immunity reduced these benefits. CONCLUSION: Even though many people are vaccinated well after September/October, they likely are still vaccinated early enough to provide substantial cost-savings.

FRED (a Framework for Reconstructing Epidemic Dynamics): an open-source software system for modeling infectious diseases and control strategies using census-based populations.

BACKGROUND: Mathematical and computational models provide valuable tools that help public health planners to evaluate competing health interventions, especially for novel circumstances that cannot be examined through observational or controlled studies, such as pandemic influenza. The spread of diseases like influenza depends on the mixing patterns within the population, and these mixing patterns depend in part on local factors including the spatial distribution and age structure of the population, the distribution of size and composition of households, employment status and commuting patterns of adults, and the size and age structure of schools. Finally, public health planners must take into account the health behavior patterns of the population, patterns that often vary according to socioeconomic factors such as race, household income, and education levels. RESULTS: FRED (a Framework for Reconstructing Epidemic Dynamics) is a freely available open-source agent-based modeling system based closely on models used in previously published studies of pandemic influenza. This version of FRED uses open-access census-based synthetic populations that capture the demographic and geographic heterogeneities of the population, including realistic household, school, and workplace social networks. FRED epidemic models are currently available for every state and county in the United States, and for selected international locations. CONCLUSIONS: State and county public health planners can use FRED to explore the effects of possible influenza epidemics in specific geographic regions of interest and to help evaluate the effect of interventions such as vaccination programs and school closure policies. FRED is available under a free open source license in order to contribute to the development of better modeling tools and to encourage open discussion of modeling tools being used to evaluate public health policies. We also welcome participation by other researchers in the further development of FRED.

Geospatial analytics to evaluate point-of-dispensing sites for mass immunizations in Allegheny County, Pennsylvania.

CONTEXT: Public health agencies use mass immunization locations to quickly administer vaccines to protect a population against an epidemic. The selection of such locations is frequently determined by available staffing levels and in some places, not all potential sites can be opened, often because of a lack of resources. Public health agencies need assistance in determining which n sites are the prime ones to open given available staff to minimize travel time and travel distance for those in the population who need to get to a site to receive treatment. OBJECTIVE: Employ geospatial analytical methods to identify the prime n locations from a predetermined set of potential locations (eg, schools) and determine which locations may not be able to achieve the throughput necessary to reach the herd immunity threshold based on varying R0 values. DESIGN: Spatial location-allocation algorithms were used to select the ideal n mass vaccination locations. SETTING: Allegheny County, Pennsylvania, served as the study area. MAIN OUTCOME MEASURES: The most favorable sites were selected and the number of individuals required to be vaccinated to achieve the herd immunity threshold for a given R0, ranging from 1.5 to 7, was determined. Locations that did not meet the Centers for Disease Control and Prevention throughput recommendation for smallpox were identified. RESULTS: At R0 = 1.5, all mass immunization locations met the required throughput to achieve the herd immunity threshold within 5 days. As R0s increased from 2 to 7, an increasing number of sites were inadequate to meet throughput requirements. CONCLUSIONS: Identifying the top n sites and categorizing those with throughput challenges allows health departments to adjust staffing, shift length, or the number of sites. This method has the potential to be expanded to select immunization locations under a number of additional scenarios.

Synthesized Population Databases: A Geospatial Database of US Poultry Farms.

The pervasive and potentially severe economic, social, and public health consequences of infectious disease in farmed animals require that plans be in place for a rapid response. Increasingly, agent-based models are being used to analyze the spread of animal-borne infectious disease outbreaks and derive policy alternatives to control future outbreaks. Although the locations, types, and sizes of animal farms are essential model inputs, no public domain nationwide geospatial database of actual farm locations and characteristics currently exists in the United States. This report describes a novel method to develop a synthetic dataset that replicates the spatial distribution of poultry farms, as well as the type and number of birds raised on them. It combines county-aggregated poultry farm counts, land use/land cover, transportation, business, and topographic data to generate locations in the conterminous United States where poultry farms are likely to be found. Simulation approaches used to evaluate the accuracy of this method when compared to that of a random placement alternative found this method to be superior. The results suggest the viability of adapting this method to simulate other livestock farms of interest to infectious disease researchers.

The role of subway travel in an influenza epidemic: a New York City simulation.

The interactions of people using public transportation in large metropolitan areas may help spread an influenza epidemic. An agent-based model computer simulation of New York City's (NYC's) five boroughs was developed that incorporated subway ridership into a Susceptible-Exposed-Infected-Recovered disease model framework. The model contains a total of 7,847,465 virtual people. Each person resides in one of the five boroughs of NYC and has a set of socio-demographic characteristics and daily behaviors that include age, sex, employment status, income, occupation, and household location and membership. The model simulates the interactions of subway riders with their workplaces, schools, households, and community activities. It was calibrated using historical data from the 1957-1958 influenza pandemics and from NYC travel surveys. The surveys were necessary to enable inclusion of subway riders into the model. The model results estimate that if influenza did occur in NYC with the characteristics of the 1957-1958 pandemic, 4% of transmissions would occur on the subway. This suggests that interventions targeted at subway riders would be relatively ineffective in containing the epidemic. A number of hypothetical examples demonstrate this feature. This information could prove useful to public health officials planning responses to epidemics.

Including the Group Quarters Population in the US Synthesized Population Database.

In 2005, RTI International researchers developed methods to generate synthesized population data on US households for the US Synthesized Population Database. These data are used in agent-based modeling, which simulates large-scale social networks to test how changes in the behaviors of individuals affect the overall network. Group quarters are residences where individuals live in close proximity and interact frequently. Although the Synthesized Population Database represents the population living in households, data for the nation's group quarters residents are not easily quantified because of US Census Bureau reporting methods designed to protect individuals' privacy.Including group quarters population data can be an important factor in agent-based modeling because the number of residents and the frequency of their interactions are variables that directly affect modeling results. Particularly with infectious disease modeling, the increased frequency of agent interaction may increase the probability of infectious disease transmission between individuals and the probability of disease outbreaks.This report reviews our methods to synthesize data on group quarters residents to match US Census Bureau data. Our goal in developing the Group Quarters Population Database was to enable its use with RTI's US Synthesized Population Database in the Modeling of Infectious Diseases Agent Study.

The benefits to all of ensuring equal and timely access to influenza vaccines in poor communities.

When influenza vaccines are in short supply, allocating vaccines equitably among different jurisdictions can be challenging. But justice is not the only reason to ensure that poorer counties have the same access to influenza vaccines as do wealthier ones. Using a detailed computer simulation model of the Washington, D.C., metropolitan region, we found that limiting or delaying vaccination of residents of poorer counties could raise the total number of influenza infections and the number of new infections per day at the peak of an epidemic throughout the region-even in the wealthier counties that had received more timely and abundant vaccine access. Among other underlying reasons, poorer counties tend to have high-density populations and more children and other higher-risk people per household, resulting in more interactions and both increased transmission of influenza and greater risk for worse influenza outcomes. Thus, policy makers across the country, in poor and wealthy areas alike, have an incentive to ensure that poorer residents have equal access to vaccines.

Would school closure for the 2009 H1N1 influenza epidemic have been worth the cost?: a computational simulation of Pennsylvania.

BACKGROUND: During the 2009 H1N1 influenza epidemic, policy makers debated over whether, when, and how long to close schools. While closing schools could have reduced influenza transmission thereby preventing cases, deaths, and health care costs, it may also have incurred substantial costs from increased childcare needs and lost productivity by teachers and other school employees. METHODS: A combination of agent-based and Monte Carlo economic simulation modeling was used to determine the cost-benefit of closing schools (vs. not closing schools) for different durations (range: 1 to 8 weeks) and symptomatic case incidence triggers (range: 1 to 30) for the state of Pennsylvania during the 2009 H1N1 epidemic. Different scenarios varied the basic reproductive rate (R(0)) from 1.2, 1.6, to 2.0 and used case-hospitalization and case-fatality rates from the 2009 epidemic. Additional analyses determined the cost per influenza case averted of implementing school closure. RESULTS: For all scenarios explored, closing schools resulted in substantially higher net costs than not closing schools. For R(0) = 1.2, 1.6, and 2.0 epidemics, closing schools for 8 weeks would have resulted in median net costs of $21.0 billion (95% Range: $8.0 - $45.3 billion). The median cost per influenza case averted would have been $14,185 ($5,423 - $30,565) for R(0) = 1.2, $25,253 ($9,501 - $53,461) for R(0) = 1.6, and $23,483 ($8,870 - $50,926) for R(0) = 2.0. CONCLUSIONS: Our study suggests that closing schools during the 2009 H1N1 epidemic could have resulted in substantial costs to society as the potential costs of lost productivity and childcare could have far outweighed the cost savings in preventing influenza cases.

Vaccination deep into a pandemic wave potential mechanisms for a "third wave" and the impact of vaccination.

BACKGROUND: In December 2009, when the H1N1 influenza pandemic appeared to be subsiding, public health officials and unvaccinated individuals faced the question of whether continued H1N1 immunization was still worthwhile. PURPOSE: To delineate what combinations of possible mechanisms could generate a third pandemic wave and then explore whether vaccinating the population at different rates and times would mitigate the wave. METHODS: As part of ongoing work with the Office of the Assistant Secretary for Preparedness and Response at the USDHHS during the H1N1 influenza pandemic, the University of Pittsburgh Models of Infectious Disease Agent Study team employed an agent-based computer simulation model of the Washington DC metropolitan region to delineate what mechanisms could generate a "third pandemic wave" and explored whether vaccinating the population at different rates and times would mitigate the wave. This model included explicit representations of the region's individuals, school systems, workplaces/commutes, households, and communities. RESULTS: Three mechanisms were identified that could cause a third pandemic wave; substantially increased viral transmissibility from seasonal forcing (changing influenza transmission with changing environmental conditions, i.e., seasons) and progressive viral adaptation; an immune escape variant; and changes in social mixing from holiday school closures. Implementing vaccination for these mechanisms, even during the down-slope of the fall epidemic wave, significantly mitigated the third wave. Scenarios showed the gains from initiating vaccination earlier, increasing the speed of vaccination, and prioritizing population subgroups based on Advisory Committee on Immunization Practices recommendations. CONCLUSIONS: Additional waves in an epidemic can be mitigated by vaccination even when an epidemic appears to be waning.

A computer simulation of vaccine prioritization, allocation, and rationing during the 2009 H1N1 influenza pandemic.

In the fall 2009, the University of Pittsburgh Models of Infectious Disease Agent Study (MIDAS) team employed an agent-based computer simulation model (ABM) of the greater Washington, DC, metropolitan region to assist the Office of the Assistant Secretary of Public Preparedness and Response, Department of Health and Human Services, to address several key questions regarding vaccine allocation during the 2009 H1N1 influenza pandemic, including comparing a vaccinating children (i.e., highest transmitters)-first policy versus the Advisory Committee on Immunization Practices (ACIP)-recommended vaccinating at-risk individuals-first policy. Our study supported adherence to the ACIP (instead of a children-first policy) prioritization recommendations for the H1N1 influenza vaccine when vaccine is in limited supply and that within the ACIP groups, children should receive highest priority.

Molecular complexity of successive bacterial epidemics deconvoluted by comparative pathogenomics.

Understanding the fine-structure molecular architecture of bacterial epidemics has been a long-sought goal of infectious disease research. We used short-read-length DNA sequencing coupled with mass spectroscopy analysis of SNPs to study the molecular pathogenomics of three successive epidemics of invasive infections involving 344 serotype M3 group A Streptococcus in Ontario, Canada. Sequencing the genome of 95 strains from the three epidemics, coupled with analysis of 280 biallelic SNPs in all 344 strains, revealed an unexpectedly complex population structure composed of a dynamic mixture of distinct clonally related complexes. We discovered that each epidemic is dominated by micro- and macrobursts of multiple emergent clones, some with distinct strain genotype-patient phenotype relationships. On average, strains were differentiated from one another by only 49 SNPs and 11 insertion-deletion events (indels) in the core genome. Ten percent of SNPs are strain specific; that is, each strain has a unique genome sequence. We identified nonrandom temporal-spatial patterns of strain distribution within and between the epidemic peaks. The extensive full-genome data permitted us to identify genes with significantly increased rates of nonsynonymous (amino acid-altering) nucleotide polymorphisms, thereby providing clues about selective forces operative in the host. Comparative expression microarray analysis revealed that closely related strains differentiated by seemingly modest genetic changes can have significantly divergent transcriptomes. We conclude that enhanced understanding of bacterial epidemics requires a deep-sequencing, geographically centric, comparative pathogenomics strategy.

Protecting health care workers: a pandemic simulation based on Allegheny County.

BACKGROUND AND OBJECTIVES: The Advisory Committee on Immunization Practices has identified health care workers (HCWs) as a priority group to receive influenza vaccine. Although the importance of HCW to the health care system is well understood, the potential role of HCW in transmission during an epidemic has not been clearly established. METHODS: Using a standard SIR (Susceptible-Infected-Recovered) framework similar to previously developed pandemic models, we developed an agent-based model (ABM) of Allegheny County, PA, that incorporates the key health care system features to simulate the spread of an influenza epidemic and its effect on hospital-based HCWs. FINDINGS: Our simulation runs found the secondary attack rate among unprotected HCWs to be approximately 60% higher (54.3%) as that of all adults (34.1%), which would result in substantial absenteeism and additional risk to HCW families. Understanding how a pandemic may affect HCWs, who must be available to treat infected patients as well as patients with other medical conditions, is crucial to policy makers' and hospital administrators' preparedness planning.

A computer simulation of employee vaccination to mitigate an influenza epidemic.

BACKGROUND: Better understanding the possible effects of vaccinating employees is important and can help policymakers and businesses plan vaccine distribution and administration logistics, especially with the current H1N1 influenza vaccine in short supply. PURPOSE: This article aims to determine the effects of varying vaccine coverage, compliance, administration rates, prioritization, and timing among employees during an influenza pandemic. METHODS: As part of the H1N1 influenza planning efforts of the Models of Infectious Disease Agent Study network, an agent-based computer simulation model was developed for the Washington DC metropolitan region, encompassing five metropolitan statistical areas. Each simulation run involved introducing 100 infectious individuals to initiate a 1.3 reproductive-rate (R(0)) epidemic, consistent with H1N1 parameters to date. Another set of scenarios represented a R(0)=1.6 epidemic. RESULTS: An unmitigated epidemic resulted in substantial productivity losses (a mean of $112.6 million for a serologic 15% attack rate and $193.8 million for a serologic 25% attack rate), even with the relatively low estimated mortality impact of H1N1. Although vaccinating Advisory Committee on Immunization Practices-defined priority groups resulted in the largest savings, vaccinating all remaining workers captured additional savings and, in fact, reduced healthcare workers' and critical infrastructure workers' chances of infection. Moreover, although employee vaccination compliance affected the epidemic, once 20% compliance was achieved, additional increases in compliance provided less incremental benefit. Even though a vast majority of the workplaces in the DC metropolitan region had fewer than 100 employees, focusing on vaccinating only those in larger firms (> or =100 employees) was just as effective in mitigating the epidemic as trying to vaccinate employees in all workplaces. CONCLUSIONS: Timely vaccination of at least 20% of the large-company workforce can play an important role in epidemic mitigation.

Attribute Assignment to a Synthetic Population in Support of Agent-Based Disease Modeling.

Communicable-disease transmission models are useful for the testing of prevention and intervention strategies. Agent-based models (ABMs) represent a new and important class of the many types of disease transmission models in use. Agent-based disease models benefit from their ability to assign disease transmission probabilities based on characteristics shared by individual agents. These shared characteristics allow ABMs to apply transmission probabilities when agents come together in geographic space. Modeling these types of social interactions requires data, and the results of the model largely depend on the quality of these input data. We initially generated a synthetic population for the United States, in support of the Models of Infectious Disease Agent Study. Subsequently, we created shared characteristics to use in ABMs. The specific goals for this task were to assign the appropriately aged populations to schools, workplaces, and public transit. Each goal presented its own challenges and problems; therefore, we used different techniques to create each type of shared characteristic. These shared characteristics have allowed disease models to more realistically predict the spread of disease, both spatially and temporally.

Simulating school closure strategies to mitigate an influenza epidemic.

BACKGROUND: There remains substantial debate over the impact of school closure as a mitigation strategy during an influenza pandemic. The ongoing 2009 H1N1 influenza pandemic has provided an unparalleled opportunity to test interventions with the most up-to-date simulations. METHODS: To assist the Allegheny County Health Department during the 2009 H1N1 influenza pandemic, the University of Pittsburgh Models of Infectious Disease Agents Study group employed an agent-based computer simulation model (ABM) of Allegheny County, Pennsylvania, to explore the effects of various school closure strategies on mitigating influenza epidemics of different reproductive rates (R0). RESULTS: Entire school system closures were not more effective than individual school closures. Any type of school closure may need to be maintained throughout most of the epidemic (ie, at least 8 weeks) to have any significant effect on the overall serologic attack rate. In fact, relatively short school closures (ie, 2 weeks or less) may actually slightly increase the overall attack rate by returning susceptible students back into schools in the middle of the epidemic. Varying the illness threshold at which school closures are triggered did not seem to have substantial impact on the effectiveness of school closures, suggesting that short delays in closing schools should not cause concern. CONCLUSIONS: School closures alone may not be able to quell an epidemic but, when maintained for at least 8 weeks, could delay the epidemic peak for up to a week, providing additional time to implement a second more effective intervention such as vaccination.

Mapping of Sample Collection Data: GIS Tools for the Natural Product Researcher.

Scientists engaged in the research of natural products often either conduct field collections themselves or collaborate with partners who do, such as botanists, mycologists, or SCUBA divers. The information gleaned from such collecting trips (e.g. longitude/latitude coordinates, geography, elevation, and a multitude of other field observations) have provided valuable data to the scientific community (e.g., biodiversity), even if it is tangential to the direct aims of the natural products research, which are often focused on drug discovery and/or chemical ecology. Geographic Information Systems (GIS) have been used to display, manage, and analyze geographic data, including collection sites for natural products. However, to the uninitiated, these tools are often beyond the financial and/or computational means of the natural product scientist. With new, free, and easy-to-use geospatial visualization tools, such as Google Earth, mapping and geographic imaging of sampling data are now within the reach of natural products scientists. The goals of the present study were to develop simple tools that are tailored for the natural products setting, thereby presenting a means to map such information, particularly via open source software like Google Earth.

Using Geographic Information Systems to Define and Map Commuting Patterns as Inputs to Agent-Based Models.

By understanding the movement patterns of people, mathematical modelers can develop models that can better analyze and predict the spread of infectious diseases. People can come into close contact in their workplaces. This report describes methods to develop georeferenced commuting patterns that can be used to characterize the work-related movement of US populations and help agent-based modelers predict workplace contacts that result in disease transmission. We used a census data product called "Census Spatial Tabulation: Census Track of Work by Census Tract of Residence (STP64)" as the data source to develop commuting pattern data for agent-based synthesized populations databases and to develop map products to visualize commuting patterns in the United States. The three primary maps we developed show inbound, outbound, and net change levels of inbound versus outbound commuters by census tract for the year 2000. Net change counts of commuters are visualized as elevations. The results can be used to quantify and assign commuting patterns of synthesized populations among different census tracts.

Synthesized Population Databases: A US Geospatial Database for Agent-Based Models.

Agent-based models simulate large-scale social systems. They assign behaviors and activities to "agents" (individuals) within the population being modeled and then allow the agents to interact with the environment and each other in complex simulations. Agent-based models are frequently used to simulate infectious disease outbreaks, among other uses.RTI used and extended an iterative proportional fitting method to generate a synthesized, geospatially explicit, human agent database that represents the US population in the 50 states and the District of Columbia in the year 2000. Each agent is assigned to a household; other agents make up the household occupants.For this database, RTI developed the methods for generating synthesized households and personsassigning agents to schools and workplaces so that complex interactions among agents as they go about their daily activities can be taken into accountgenerating synthesized human agents who occupy group quarters (military bases, college dormitories, prisons, nursing homes).In this report, we describe both the methods used to generate the synthesized population database and the final data structure and data content of the database. This information will provide researchers with the information they need to use the database in developing agent-based models.Portions of the synthesized agent database are available to any user upon request. RTI will extract a portion (a county, region, or state) of the database for users who wish to use this database in their own agent-based models.

Antioxidant activity and total phenolic content of aqueous and methanolic extracts of Jordanian plants: an ICBG project.

As part of an International Cooperative Biodiversity Groups (ICBG) program to study Jordan's biodiversity, the relative levels of antioxidant activity and the total phenolic content of aqueous and methanolic extracts of a total of 95 plant species, all of Jordanian origin and those collected at random, have been measured. The total phenolic content of aqueous and methanolic extracts of the investigated plant species ranged from 4.4 to 78.3 mg and from 2.1 to 52.8 mg gallic acid equivalents g(-1) dry weight, respectively, while the total antioxidant capacity ranged from 20.0 to 916.7 and from 15.1 to 915.6 micromol Trolox equivalents g(-1) dry weight, respectively. Based on this collection, approximately 5% of assayed plants showed high levels of antioxidant activity. There was a significant linear correlation between antioxidant activity and total phenolic content for aqueous and methanolic extracts, suggesting that phenolic compounds were the predominant antioxidant components in the investigated plant species. Interestingly, a few of the collected plants had high-antioxidant activity yet "low" phenolic content includes Ceratonia siliqua and Viscum cruciatum. These plants may serve as sources of antioxidants with new chemotypes.