Model predictions and evaluation of possible control strategies for the 2009 A/H1N1v influenza pandemic in Italy

Characterizing the transmission patterns of seasonal influenza in Italy: lessons from the last decade

BACKGROUND

Despite thousands of influenza cases annually recorded by surveillance systems around the globe, estimating the transmission patterns of seasonal influenza is challenging.

METHODS

We develop an age-structured mathematical model to influenza transmission to analyze ten consecutive seasons (from 2010 to 2011 to 2019-2020) of influenza epidemiological and virological data reported to the Italian surveillance system.

RESULTS

We estimate that 18.4-29.3% of influenza infections are detected by the surveillance system. Influenza infection attack rate varied between 12.7 and 30.5% and is generally larger for seasons characterized by the circulation of A/H3N2 and/or B types/subtypes. Individuals aged 14 years or less are the most affected age-segment of the population, with A viruses especially affecting children aged 0-4 years. For all influenza types/subtypes, the mean effective reproduction number is estimated to be generally in the range 1.09-1.33 (9 out of 10 seasons) and never exceeding 1.41. The age-specific susceptibility to infection appears to be a type/subtype-specific feature.

CONCLUSIONS

The results presented in this study provide insights on type/subtype-specific transmission patterns of seasonal influenza that could be instrumental to fine-tune immunization strategies and non-pharmaceutical interventions aimed at limiting seasonal influenza spread and burden.

The roles of medical examiners in the COVID-19 era: a comparison between the United States and Italy

Italy and the United States are two of the countries most affected by SARS-CoV-2 (COVID-19), with more than 240,760 confirmed cases in Italy and 2,699,658 in the United States (as of July 2, 2020). The current COVID-19 pandemic has led to substantial changes in many fields of medicine, specifically in the forensic discipline. Medicolegal activities related to conducting autopsies have been largely affected by the COVID-19 pandemic. Postmortem examinations are generally discouraged by government regulations due to the risk of spreading the disease further through the handling and dissection of bodies from patients who succumbed to COVID-19 infection. There is a paucity of data regarding the persistence of SARS-CoV-2 in bodies, as well as concerning the reliability of swabbing methods in human remains. On the other hand, the autopsy is an essential tool to provide necessary information about the pathophysiology of the disease that presents useful clinical and epidemiological insights. On this basis, we aim to address issues concerning general medical examiner/coroner organization, comparing the Italian and American systems. We also discuss the pivotal roles of forensic pathologists in informing infectious disease surveillance. Finally, we focus on the impact of COVID-19 emergency on medicolegal practices in Italy and the United States, as well as the responses of the forensic scientific community to the emerging concerns related to the pandemic. We believe that stronger efforts by authorities are necessary to facilitate completing postmortem examinations, as data derived from such assessments are expected to be paramount to improving patient management and disease prevention.

Multi-model forecasts of the ongoing Ebola epidemic in the Democratic Republic of Congo, March-October 2019

The 2018–2020 Ebola outbreak in the Democratic Republic of the Congo is the first to occur in an armed conflict zone. The resulting impact on population movement, treatment centres and surveillance has created an unprecedented challenge for real-time epidemic forecasting. Most standard mathematical models cannot capture the observed incidence trajectory when it deviates from a traditional epidemic logistic curve. We fit seven dynamic models of increasing complexity to the incidence data published in the World Health Organization Situation Reports, after adjusting for reporting delays. These models include a simple logistic model, a Richards model, an endemic Richards model, a double logistic growth model, a multi-model approach and two sub-epidemic models. We analyse model fit to the data and compare real-time forecasts throughout the ongoing epidemic across 29 weeks from 11 March to 23 September 2019. We observe that the modest extensions presented allow for capturing a wide range of epidemic behaviour. The multi-model approach yields the most reliable forecasts on average for this application, and the presented extensions improve model flexibility and forecasting accuracy, even in the context of limited epidemiological data.

Evaluation of the COVID-19 Pandemic Intervention Strategies with Hesitant F-AHP

In this study, a hesitant fuzzy AHP method is presented to help decision makers (DMs), especially policymakers, governors, and physicians, evaluate the importance of intervention strategy alternatives applied by various countries for the COVID-19 pandemic. In this research, a hesitant fuzzy multicriteria decision making (MCDM) method, hesitant fuzzy Analytic Hierarchy Process (hesitant F-AHP), is implemented to make pairwise comparison of COVID-19 country-level intervention strategies applied by various countries and determine relative importance scores. An illustrative study is presented where fifteen intervention strategies applied by various countries in the world during the COVID-19 pandemic are evaluated by seven physicians (a professor of infectious diseases and clinical microbiology, an infectious disease physician, a clinical microbiology physician, two internal medicine physicians, an anesthesiology and reanimation physician, and a family physician) in Turkey who act as DMs in the process.

Systematic selection between age and household structure for models aimed at emerging epidemic predictions

Numerous epidemic models have been developed to capture aspects of human contact patterns, making model selection challenging when they fit (often-scarce) early epidemic data equally well but differ in predictions. Here we consider the invasion of a novel directly transmissible infection and perform an extensive, systematic and transparent comparison of models with explicit age and/or household structure, to determine the accuracy loss in predictions in the absence of interventions when ignoring either or both social components. We conclude that, with heterogeneous and assortative contact patterns relevant to respiratory infections, the model’s age stratification is crucial for accurate predictions. Conversely, the household structure is only needed if transmission is highly concentrated in households, as suggested by an empirical but robust rule of thumb based on household secondary attack rate. This work serves as a template to guide the simplicity/accuracy trade-off in designing models aimed at initial, rapid assessment of potential epidemic severity.

Models of emerging epidemics can be exceedingly helpful in planning the response, but early on model selection is a difficult task. Here, the authors explore the joint contribution of age stratification and household structure on epidemic spread, and provides a rule of thumb to guide model choice.

Reactive school closure weakens the network of social interactions and reduces the spread of influenza

School-closure policies are considered one of the most promising nonpharmaceutical interventions for mitigating seasonal and pandemic influenza. However, their effectiveness is still debated, primarily due to the lack of empirical evidence about the behavior of the population during the implementation of the policy. Over the course of the 2015 to 2016 influenza season in Russia, we performed a diary-based contact survey to estimate the patterns of social interactions before and during the implementation of reactive school-closure strategies. We develop an innovative hybrid survey-modeling framework to estimate the time-varying network of human social interactions. By integrating this network with an infection transmission model, we reduce the uncertainty surrounding the impact of school-closure policies in mitigating the spread of influenza. When the school-closure policy is in place, we measure a significant reduction in the number of contacts made by students (14.2 vs. 6.5 contacts per day) and workers (11.2 vs. 8.7 contacts per day). This reduction is not offset by the measured increase in the number of contacts between students and nonhousehold relatives. Model simulations suggest that gradual reactive school-closure policies based on monitoring student absenteeism rates are capable of mitigating influenza spread. We estimate that without the implemented reactive strategies the attack rate of the 2015 to 2016 influenza season would have been 33% larger. Our study sheds light on the social mixing patterns of the population during the implementation of reactive school closures and provides key instruments for future cost-effectiveness analyses of school-closure policies.

Real-time predictive seasonal influenza model in Catalonia, Spain

Influenza surveillance is critical to monitoring the situation during epidemic seasons and predictive mathematic models may aid the early detection of epidemic patterns. The objective of this study was to design a real-time spatial predictive model of ILI (Influenza Like Illness) incidence rate in Catalonia using one- and two-week forecasts. The available data sources used to select explanatory variables to include in the model were the statutory reporting disease system and the sentinel surveillance system in Catalonia for influenza incidence rates, the official climate service in Catalonia for meteorological data, laboratory data and Google Flu Trend. Time series for every explanatory variable with data from the last 4 seasons (from 2010–2011 to 2013–2014) was created. A pilot test was conducted during the 2014–2015 season to select the explanatory variables to be included in the model and the type of model to be applied. During the 2015–2016 season a real-time model was applied weekly, obtaining the intensity level and predicted incidence rates with 95% confidence levels one and two weeks away for each health region. At the end of the season, the confidence interval success rate (CISR) and intensity level success rate (ILSR) were analysed. For the 2015–2016 season a CISR of 85.3% at one week and 87.1% at two weeks and an ILSR of 82.9% and 82% were observed, respectively. The model described is a useful tool although it is hard to evaluate due to uncertainty. The accuracy of prediction at one and two weeks was above 80% globally, but was lower during the peak epidemic period. In order to improve the predictive power, new explanatory variables should be included.

Detecting a Surprisingly Low Transmission Distance in the Early Phase of the 2009 Influenza Pandemic

The spread of the 2009 H1N1 influenza pandemic in England was characterized by two major waves of infections: the first one was highly spatially localized (mainly in the London area), while the second one spread homogeneously through the entire country. The reasons behind this complex spatiotemporal dynamics have yet to be clarified. In this study, we perform a Bayesian analysis of five models entailing different hypotheses on the possible determinants of the observed pattern. We find a consensus among all models in showing a surprisingly low transmission distance (defined as the geographic distance between the place of residence of the infectors and her/his infectees) during the first wave: about 1.5 km (2.2 km if infections linked to household and school transmission are excluded). The best-fitting model entails a change in human activity regarding contacts not related to household and school. By using this model we estimate that the transmission distance sharply increased to 5.3 km (10 km when excluding infections linked to household and school transmission) during the second wave. Our study reveals a possible explanation for the observed pattern and highlights the need of better understanding human mobility and activity patterns under the pressure posed by a pandemic threat.

Recalibrating disease parameters for increasing realism in modeling epidemics in closed settings

Background

The homogeneous mixing assumption is widely adopted in epidemic modelling for its parsimony and represents the building block of more complex approaches, including very detailed agent-based models. The latter assume homogeneous mixing within schools, workplaces and households, mostly for the lack of detailed information on human contact behaviour within these settings. The recent data availability on high-resolution face-to-face interactions makes it now possible to assess the goodness of this simplified scheme in reproducing relevant aspects of the infection dynamics.

Methods

We consider empirical contact networks gathered in different contexts, as well as synthetic data obtained through realistic models of contacts in structured populations. We perform stochastic spreading simulations on these contact networks and in populations of the same size under a homogeneous mixing hypothesis. We adjust the epidemiological parameters of the latter in order to fit the prevalence curve of the contact epidemic model. We quantify the agreement by comparing epidemic peak times, peak values, and epidemic sizes.

Results

Good approximations of the peak times and peak values are obtained with the homogeneous mixing approach, with a median relative difference smaller than 20 % in all cases investigated. Accuracy in reproducing the peak time depends on the setting under study, while for the peak value it is independent of the setting. Recalibration is found to be linear in the epidemic parameters used in the contact data simulations, showing changes across empirical settings but robustness across groups and population sizes.

Conclusions

An adequate rescaling of the epidemiological parameters can yield a good agreement between the epidemic curves obtained with a real contact network and a homogeneous mixing approach in a population of the same size. The use of such recalibrated homogeneous mixing approximations would enhance the accuracy and realism of agent-based simulations and limit the intrinsic biases of the homogeneous mixing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12879-016-2003-3) contains supplementary material, which is available to authorized users.

School closure policies at municipality level for mitigating influenza spread: a model-based evaluation

Background

Nearly every year Influenza affects most countries worldwide and the risk of a new pandemic is always present. Therefore, influenza is a major concern for public health. School-age individuals are often the most affected group, suggesting that the inclusion in preparedness plans of school closure policies may represent an option for influenza mitigation. However, their applicability remains uncertain and their implementation should carefully be weighed on the basis of cost-benefit considerations.

Methods

We developed an individual-based model for influenza transmission integrating data on sociodemography and time use of the Italian population, face-to-face contacts in schools, and influenza natural history. The model was calibrated on the basis of epidemiological data from the 2009 influenza pandemic and was used to evaluate the effectiveness of three reactive school closure strategies, all based on school absenteeism.

Results

In the case of a new influenza pandemic sharing similar features with the 2009 H1N1 pandemic, gradual school closure strategies (i.e., strategies closing classes first, then grades or the entire school) could lead to attack rate reduction up to 20–25 % and to peak weekly incidence reduction up to 50–55 %, at the cost of about three school weeks lost per student. Gradual strategies are quite stable to variations in the start of policy application and to the threshold on student absenteeism triggering class (and school) closures. In the case of a new influenza pandemic showing different characteristics with respect to the 2009 H1N1 pandemic, we found that the most critical features determining the effectiveness of school closure policies are the reproduction number and the age-specific susceptibility to infection, suggesting that these two epidemiological quantities should be estimated early on in the spread of a new pandemic for properly informing response planners.

Conclusions

Our results highlight a potential beneficial effect of reactive gradual school closure policies in mitigating influenza spread, conditioned on the effort that decision makers are willing to afford. Moreover, the suggested strategies are solely based on routinely collected and easily accessible data (such as student absenteeism irrespective of the cause and ILI incidence) and thus they appear to be applicable in real world situations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12879-016-1918-z) contains supplementary material, which is available to authorized users.

Estimating transmission probability in schools for the 2009 H1N1 influenza pandemic in Italy

BACKGROUND

Epidemic models are being extensively used to understand the main pathways of spread of infectious diseases, and thus to assess control methods. Schools are well known to represent hot spots for epidemic spread; hence, understanding typical patterns of infection transmission within schools is crucial for designing adequate control strategies. The attention that was given to the 2009 A/H1N1pdm09 flu pandemic has made it possible to collect detailed data on the occurrence of influenza-like illness (ILI) symptoms in two primary schools of Trento, Italy.

RESULTS

The data collected in the two schools were used to calibrate a discrete-time SIR model, which was designed to estimate the probabilities of influenza transmission within the classes, grades and schools using Markov Chain Monte Carlo (MCMC) methods. We found that the virus was mainly transmitted within class, with lower levels of transmission between students in the same grade and even lower, though not significantly so, among different grades within the schools. We estimated median values of R 0 from the epidemic curves in the two schools of 1.16 and 1.40; on the other hand, we estimated the average number of students infected by the first school case to be 0.85 and 1.09 in the two schools.

CONCLUSIONS

The discrepancy between the values of R 0 estimated from the epidemic curve or from the within-school transmission probabilities suggests that household and community transmission played an important role in sustaining the school epidemics. The high probability of infection between students in the same class confirms that targeting within-class transmission is key to controlling the spread of influenza in school settings and, as a consequence, in the general population.

Measuring the potential of individual airports for pandemic spread over the world airline network

Background

Massive growth in human mobility has dramatically increased the risk and rate of pandemic spread. Macro-level descriptors of the topology of the World Airline Network (WAN) explains middle and late stage dynamics of pandemic spread mediated by this network, but necessarily regard early stage variation as stochastic. We propose that much of this early stage variation can be explained by appropriately characterizing the local network topology surrounding an outbreak’s debut location.

Methods

Based on a model of the WAN derived from public data, we measure for each airport the expected force of infection (AEF) which a pandemic originating at that airport would generate, assuming an epidemic process which transmits from airport to airport via scheduled commercial flights. We observe, for a subset of world airports, the minimum transmission rate at which a disease becomes pandemically competent at each airport. We also observe, for a larger subset, the time until a pandemically competent outbreak achieves pandemic status given its debut location. Observations are generated using a highly sophisticated metapopulation reaction-diffusion simulator under a disease model known to well replicate the 2009 influenza pandemic. The robustness of the AEF measure to model misspecification is examined by degrading the underlying model WAN.

Results

AEF powerfully explains pandemic risk, showing correlation of 0.90 to the transmission level needed to give a disease pandemic competence, and correlation of 0.85 to the delay until an outbreak becomes a pandemic. The AEF is robust to model misspecification. For 97 % of airports, removing 15 % of airports from the model changes their AEF metric by less than 1 %.

Conclusions

Appropriately summarizing the size, shape, and diversity of an airport’s local neighborhood in the WAN accurately explains much of the macro-level stochasticity in pandemic outcomes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12879-016-1350-4) contains supplementary material, which is available to authorized users.

Model-Based Comprehensive Analysis of School Closure Policies for Mitigating Influenza Epidemics and Pandemics

School closure policies are among the non-pharmaceutical measures taken into consideration to mitigate influenza epidemics and pandemics spread. However, a systematic review of the effectiveness of alternative closure policies has yet to emerge. Here we perform a model-based analysis of four types of school closure, ranging from the nationwide closure of all schools at the same time to reactive gradual closure, starting from class-by-class, then grades and finally the whole school. We consider policies based on triggers that are feasible to monitor, such as school absenteeism and national ILI surveillance system. We found that, under specific constraints on the average number of weeks lost per student, reactive school-by-school, gradual, and county-wide closure give comparable outcomes in terms of optimal infection attack rate reduction, peak incidence reduction or peak delay. Optimal implementations generally require short closures of one week each; this duration is long enough to break the transmission chain without leading to unnecessarily long periods of class interruption. Moreover, we found that gradual and county closures may be slightly more easily applicable in practice as they are less sensitive to the value of the excess absenteeism threshold triggering the start of the intervention. These findings suggest that policy makers could consider school closure policies more diffusely as response strategy to influenza epidemics and pandemics, and the fact that some countries already have some experience of gradual or regional closures for seasonal influenza outbreaks demonstrates that logistic and feasibility challenges of school closure strategies can be to some extent overcome.

Connecting within and between-hosts dynamics in the influenza infection-staged epidemiological models with behavior change

Influenza viruses are a major public health problem worldwide. Although influenza has been extensively researched, there are still many aspects that are not fully understood such as the effects of within and between-hosts dynamics and their impact on behavior change. Here, we develop mathematical models with multiple infection stages and estimate parameters based on within-host data to investigate the impact of behavior change on influenza dynamics. We divide the infected population into three and four groups based on the age of the infection, which corresponds to viral load shedding. We consider within-host data on viral shedding to estimate the length and force of infection of the different infectivity stages. Our results show that behavior changes, due to exogenous events (e.g., media coverage) and disease symptoms, are effective in delaying and lowering an epidemic peak. We show that the dynamics of viral shedding and symptoms, during the infection, are key features when considering epidemic prevention strategies. This study improves our understanding of the spread of influenza virus infection in the population and provides information about the impact of emergent behavior and its connection to the within and between-hosts dynamics.

The role of different social contexts in shaping influenza transmission during the 2009 pandemic

Evaluating the relative importance of different social contexts in which infection transmission occurs is critical for identifying optimal intervention strategies. Nonetheless, an overall picture of influenza transmission in different social contexts has yet to emerge. Here we provide estimates of the fraction of infections generated in different social contexts during the 2009 H1N1 pandemic in Italy by making use of a highly detailed individual-based model accounting for time use data and parametrized on the basis of observed age-specific seroprevalence. We found that 41.6% (95%CI: 39-43.7%) of infections occurred in households, 26.7% (95%CI: 21-33.2) in schools, 3.3% (95%CI: 1.7-5%) in workplaces, and 28.4% (95%CI: 24.6-31.9%) in the general community. The above estimates strongly depend on the lower susceptibility to infection of individuals 19+ years old compared to younger ones, estimated to be 0.2 (95%CI 0.12-0.28). We also found that school closure over the weekends contributed to decrease the effective reproduction number of about 8% and significantly affected the pattern of transmission. These results highlight the pivotal role played by schools in the transmission of the 2009 H1N1 influenza. They may be relevant in the evaluation of intervention options and, hence, for informing policy decisions.

Network epidemiology and plant trade networks

Models of epidemics in complex networks are improving our predictive understanding of infectious disease outbreaks. Nonetheless, applying network theory to plant pathology is still a challenge. This overview summarizes some key developments in network epidemiology that are likely to facilitate its application in the study and management of plant diseases. Recent surveys have provided much-needed datasets on contact patterns and human mobility in social networks, but plant trade networks are still understudied. Human (and plant) mobility levels across the planet are unprecedented-there is thus much potential in the use of network theory by plant health authorities and researchers. Given the directed and hierarchical nature of plant trade networks, there is a need for plant epidemiologists to further develop models based on undirected and homogeneous networks. More realistic plant health scenarios would also be obtained by developing epidemic models in dynamic, rather than static, networks. For plant diseases spread by the horticultural and ornamental trade, there is the challenge of developing spatio-temporal epidemic simulations integrating network data. The use of network theory in plant epidemiology is a promising avenue and could contribute to anticipating and preventing plant health emergencies such as European ash dieback.

Influenza forecasting in human populations: a scoping review

Forecasts of influenza activity in human populations could help guide key preparedness tasks. We conducted a scoping review to characterize these methodological approaches and identify research gaps. Adapting the PRISMA methodology for systematic reviews, we searched PubMed, CINAHL, Project Euclid, and Cochrane Database of Systematic Reviews for publications in English since January 1, 2000 using the terms "influenza AND (forecast* OR predict*)", excluding studies that did not validate forecasts against independent data or incorporate influenza-related surveillance data from the season or pandemic for which the forecasts were applied. We included 35 publications describing population-based (N = 27), medical facility-based (N = 4), and regional or global pandemic spread (N = 4) forecasts. They included areas of North America (N = 15), Europe (N = 14), and/or Asia-Pacific region (N = 4), or had global scope (N = 3). Forecasting models were statistical (N = 18) or epidemiological (N = 17). Five studies used data assimilation methods to update forecasts with new surveillance data. Models used virological (N = 14), syndromic (N = 13), meteorological (N = 6), internet search query (N = 4), and/or other surveillance data as inputs. Forecasting outcomes and validation metrics varied widely. Two studies compared distinct modeling approaches using common data, 2 assessed model calibration, and 1 systematically incorporated expert input. Of the 17 studies using epidemiological models, 8 included sensitivity analysis. This review suggests need for use of good practices in influenza forecasting (e.g., sensitivity analysis); direct comparisons of diverse approaches; assessment of model calibration; integration of subjective expert input; operational research in pilot, real-world applications; and improved mutual understanding among modelers and public health officials.

Developing model-based public health policy through knowledge translation: the need for a 'Communities of Practice'

OBJECTIVES

The 2009 influenza A (H1N1) pandemic prompted public health agencies worldwide to respond in a context of substantial uncertainty. While many lessons around successful management strategies were learned during the influenza A (H1N1) pandemic, the usefulness and impact of mathematical models to optimize policy decisions in protecting public health were poorly realized. The authors explored the experiences of modellers and public health practitioners in trying to develop model-based public health policies in the management of the 2009 influenza A (H1N1) pandemic in Canada.

STUDY DESIGN

A qualitative case study design based on interviews and other textual data was used.

METHODS

Individual interviews were conducted with mathematical modellers and public health professionals from academia and government health departments during the second wave of the 2009 influenza A (H1N1) pandemic (both prior to and following the vaccine roll-out), using a convergent interviewing process. Interviews were supplemented with discussions held during three separate workshops involving representatives from these groups on the role of modelling in pandemic preparedness and responses. NVivo9™ was used to analyse interview data and associated notes.

RESULTS

Mathematical models were underutilized during the response phase of the 2009 influenza A (H1N1) pandemic, largely because many public health professionals were unaware of modelling infrastructure in Canada. Challenges were reflected in three ways: 1) the relevance of models to public health priorities; 2) the need for clear communication and plain language around modelling and its contributions and limitations; and 3) the need for increased trust and collaboration to develop strong working relationships.

CONCLUSIONS

Developing a 'Communities of Practice' between public health professionals and mathematical modellers during inter-pandemic periods based on common targeted goals, using plain language, and where relationships between individuals and organizations are developed early, could be an effective strategy to assist the process of public health policy decision-making, particularly when characterized by high levels of uncertainty.

Pandemic influenza A/H1N1pdm in Italy: age, risk and population susceptibility

Background

A common pattern emerging from several studies evaluating the impact of the 2009 A/H1N1 pandemic influenza (A/H1N1pdm) conducted in countries worldwide is the low attack rate observed in elderly compared to that observed in children and young adults. The biological or social mechanisms responsible for the observed age-specific risk of infection are still to be deeply investigated.

Methods

The level of immunity against the A/H1N1pdm in pre and post pandemic sera was determined using left over sera taken for diagnostic purposes or routine ascertainment obtained from clinical laboratories. The antibody titres were measured by the haemagglutination inhibition (HI) assay. To investigate whether certain age groups had higher risk of infection the presence of protective antibody (≥1∶40), was calculated using exact binomial 95% CI on both pre- and post- pandemic serological data in the age groups considered. To estimate age-specific susceptibility to infection we used an age-structured SEIR model.

Results

By comparing pre- and post-pandemic serological data in Italy we found age- specific attack rates similar to those observed in other countries. Cumulative attack rate at the end of the first A/H1N1pdm season in Italy was estimated to be 16.3% (95% CI 9.4%-23.1%). Modeling results allow ruling out the hypothesis that only age-specific characteristics of the contact network and levels of pre-pandemic immunity are responsible for the observed age-specific risk of infection. This means that age-specific susceptibility to infection, suspected to play an important role in the pandemic, was not only determined by pre-pandemic levels of H1N1pdm antibody measured by HI.

Conclusions

Our results claim for new studies to better identify the biological mechanisms, which might have determined the observed pattern of susceptibility with age. Moreover, our results highlight the need to obtain early estimates of differential susceptibility with age in any future pandemics to obtain more reliable real time estimates of critical epidemiological parameters.

Real-time numerical forecast of global epidemic spreading: case study of 2009 A/H1N1pdm

BACKGROUND

Mathematical and computational models for infectious diseases are increasingly used to support public-health decisions; however, their reliability is currently under debate. Real-time forecasts of epidemic spread using data-driven models have been hindered by the technical challenges posed by parameter estimation and validation. Data gathered for the 2009 H1N1 influenza crisis represent an unprecedented opportunity to validate real-time model predictions and define the main success criteria for different approaches.

METHODS

We used the Global Epidemic and Mobility Model to generate stochastic simulations of epidemic spread worldwide, yielding (among other measures) the incidence and seeding events at a daily resolution for 3,362 subpopulations in 220 countries. Using a Monte Carlo Maximum Likelihood analysis, the model provided an estimate of the seasonal transmission potential during the early phase of the H1N1 pandemic and generated ensemble forecasts for the activity peaks in the northern hemisphere in the fall/winter wave. These results were validated against the real-life surveillance data collected in 48 countries, and their robustness assessed by focusing on 1) the peak timing of the pandemic; 2) the level of spatial resolution allowed by the model; and 3) the clinical attack rate and the effectiveness of the vaccine. In addition, we studied the effect of data incompleteness on the prediction reliability.

RESULTS

Real-time predictions of the peak timing are found to be in good agreement with the empirical data, showing strong robustness to data that may not be accessible in real time (such as pre-exposure immunity and adherence to vaccination campaigns), but that affect the predictions for the attack rates. The timing and spatial unfolding of the pandemic are critically sensitive to the level of mobility data integrated into the model.

CONCLUSIONS

Our results show that large-scale models can be used to provide valuable real-time forecasts of influenza spreading, but they require high-performance computing. The quality of the forecast depends on the level of data integration, thus stressing the need for high-quality data in population-based models, and of progressive updates of validated available empirical knowledge to inform these models.

Transmission potential and design of adequate control measures for Marburg hemorrhagic fever

Marburg hemorrhagic fever is rare yet among the most severe diseases affecting humans, with case fatality ratio even higher than 80%. By analyzing the largest documented Marburg hemorrhagic fever epidemic, which occurred in Angola in 2005 and caused 329 deaths, and data on viral load over time in non-human primates, we make an assessment of transmissibility and severity of the disease. We also give insight into the control of new Marburg hemorrhagic fever epidemics to inform appropriate health responses. We estimated the distribution of the generation time to have mean 9 days (95%CI: 8.2–10 days) and standard deviation 5.4 days (95%CI: 3.9–8.6 days), and the basic reproduction number to be  = 1.59 (95%CI: 1.53–1.66). Model simulations suggest that a timely isolation of cases, starting no later than 2–3 days after symptoms onset, is sufficient to contain an outbreak. Our analysis reveals that Marburg hemorrhagic fever is characterized by a relatively small reproduction number and by a relatively long generation time. Such factors, along with the extremely high severity and fatality, support the rare occurrence of large epidemics in human populations. Our results also support the effectiveness of social distancing measures - case isolation in particular - to contain or at least to mitigate an emerging outbreak. This work represents an advance in the knowledge required to manage a potential Marburg hemorrhagic fever epidemic.

Inferring the structure of social contacts from demographic data in the analysis of infectious diseases spread

Social contact patterns among individuals encode the transmission route of infectious diseases and are a key ingredient in the realistic characterization and modeling of epidemics. Unfortunately, the gathering of high quality experimental data on contact patterns in human populations is a very difficult task even at the coarse level of mixing patterns among age groups. Here we propose an alternative route to the estimation of mixing patterns that relies on the construction of virtual populations parametrized with highly detailed census and demographic data. We present the modeling of the population of 26 European countries and the generation of the corresponding synthetic contact matrices among the population age groups. The method is validated by a detailed comparison with the matrices obtained in six European countries by the most extensive survey study on mixing patterns. The methodology presented here allows a large scale comparison of mixing patterns in Europe, highlighting general common features as well as country-specific differences. We find clear relations between epidemiologically relevant quantities (reproduction number and attack rate) and socio-demographic characteristics of the populations, such as the average age of the population and the duration of primary school cycle. This study provides a numerical approach for the generation of human mixing patterns that can be used to improve the accuracy of mathematical models in the absence of specific experimental data.

The dynamics of infectious diseases caused by pathogens transmissible from human to human strongly depends on contact patterns between individuals. High quality observational data on contact patterns, usually presented in the form of age-specific contact matrices, are difficult to gather and are currently available only for few countries worldwide. Here we propose a computational approach, based on the simulation of a virtual society of agents, allowing the estimation of contact patterns by age for 26 European countries. We validate the estimated contact matrices against those obtained by the most extensive field study on contact patterns, with data collected in eight European countries. We show that our contact matrices share some common features, e.g. individuals tend to mix preferentially with individuals their own age, and country-specific differences, which can be partly explained by differences in population structures due to different demographic trajectories followed after WWII. Our analysis highlights well defined correlations between epidemiological parameters and socio-demographic features of the populations. This study provides the first estimates of contact matrices for many European countries where specific experimental data are still not available.

Epidemiology and transmission dynamics of the 1918-19 pandemic influenza in Florence, Italy

To investigate the 1918/19 influenza pandemic daily number of new hospitalizations in the only hospital in Florence (Central Italy) were analyzed. In order to describe the transmission dynamics of the 1918/1919 pandemic influenza a compartmental epidemic model was used. Model simulations show a high level of agreement with the observed epidemic data. By assuming both latent and infectious period equal to 1.5 days, the estimated basic reproduction number was R(0)(1) = 1.03 (95% CI: 1.00-1.08) during the summer wave and R(0)(2) = 1.38 (95% CI: 1.32-1.48) during the fall wave. Varying the length of the generation time or the estimation method, R(0)(2) ranges from 1.32 to 1.71. The hospitalization rate was found significantly different between summer and fall waves. Notably, the estimated basic reproductive numbers are lower compared to those observed in other countries, while the age distribution of deaths resulted to be consistent with the patterns generally observed during of the 1918-1919 pandemic. Our knowledge on past pandemics, as for the 1918-19 Spanish influenza, would help improving mathematical modeling accuracy and understanding the mechanisms underlying the dynamics of future pandemics.

Variability in transmissibility of the 2009 H1N1 pandemic in Canadian communities

Background

The prevalence and severity of the 2009 H1N1 pandemic appeared to vary significantly across populations and geographic regions. We sought to investigate the variability in transmissibility of H1N1 pandemic in different health regions (including urban centres and remote, isolated communities) in the province of Manitoba, Canada.

Methods

The Richards model was used to fit to the daily number of laboratory-confirmed cases and estimate transmissibility (referred to as the basic reproduction number, R₀), doubling times, and turning points of outbreaks in both spring and fall waves of the H1N1 pandemic in several health regions.

Results

We observed considerable variation in R₀ estimates ranging from 1.55 to 2.24, with confidence intervals ranging from 1.45 to 2.88, for an average generation time of 2.9 days, and shorter doubling times in some remote and isolated communities compared to urban centres, suggesting a more rapid spread of disease in these communities during the first wave. For the second wave, R_e, the effective reproduction number, is estimated to be lower for remote and isolated communities; however, outbreaks appear to have been driven by somewhat higher transmissibility in urban centres.

Conclusions

There was considerable geographic variation in transmissibility of the 2009 pandemic outbreaks. While highlighting the importance of estimating R₀ for informing health responses, the findings indicate that projecting the transmissibility for large-scale epidemics may not faithfully characterize the early spread of disease in remote and isolated communities.

A new approach to characterising infectious disease transmission dynamics from sentinel surveillance: application to the Italian 2009-2010 A/H1N1 influenza pandemic

Syndromic and virological data are routinely collected by many countries and are often the only information available in real time. The analysis of surveillance data poses many statistical challenges that have not yet been addressed. For instance, the fraction of cases that seek healthcare and are thus detected is often unknown. Here, we propose a general statistical framework that explicitly takes into account the way the surveillance data are generated. Our approach couples a deterministic mathematical model with a statistical description of the reporting process and is applied to surveillance data collected in Italy during the 2009-2010 A/H1N1 influenza pandemic. We estimate that the reproduction number R was initially into the range 1.2-1.4 and that case detection in children was significantly higher than in adults. According to the best fit models, we estimate that school-age children experienced the highest infection rate overall. In terms of both estimated peak-incidence and overall attack rate, according to the Susceptibility and Immunity models the 5-14 years age-class was about 5 times more infected than the 65+ years old age-group and about twice more than the 15-64 years age-class. The multiplying factors are doubled using the Baseline model. Overall, the estimated attack rate was about 16% according to the Baseline model and 30% according to the Susceptibility and Immunity models.

Determinants of the spatiotemporal dynamics of the 2009 H1N1 pandemic in Europe: implications for real-time modelling

Influenza pandemics in the last century were characterized by successive waves and differences in impact and timing between different regions, for reasons not clearly understood. The 2009 H1N1 pandemic showed rapid global spread, but with substantial heterogeneity in timing within each hemisphere. Even within Europe substantial variation was observed, with the UK being unique in experiencing a major first wave of transmission in early summer and all other countries having a single major epidemic in the autumn/winter, with a West to East pattern of spread. Here we show that a microsimulation model, parameterised using data about H1N1pdm collected by the beginning of June 2009, explains the occurrence of two waves in UK and a single wave in the rest of Europe as a consequence of timing of H1N1pdm spread, fluxes of travels from US and Mexico, and timing of school vacations. The model provides a description of pandemic spread through Europe, depending on intra-European mobility patterns and socio-demographic structure of the European populations, which is in broad agreement with observed timing of the pandemic in different countries. Attack rates are predicted to depend on the socio-demographic structure, with age dependent attack rates broadly agreeing with available serological data. Results suggest that the observed heterogeneity can be partly explained by the between country differences in Europe: marked differences in school calendars, mobility patterns and sociodemographic structures. Moreover, higher susceptibility of children to infection played a key role in determining the epidemiology of the 2009 pandemic. Our work shows that it would have been possible to obtain a broad-brush prediction of timing of the European pandemic well before the autumn of 2009, much more difficult to achieve with simpler models or pre-pandemic parameterisation. This supports the use of models accounting for the structure of complex modern societies for giving insight to policy makers.

The 2009 H1N1pdm influenza pandemic spread rapidly but heterogeneously. A notable pattern occurred in Europe, with the UK exhibiting a first wave in early summer and a second wave in autumn, while all other European countries experienced a single wave in autumn/winter, resulting in a clear West to East pattern of spread. Our study asks which factors were most responsible for this variation, and to what extent the pattern of spread was predictable from data available in the first two months of the pandemic. Providing reliable answers to these questions would reduce uncertainty and improve situational awareness for policy-makers in the future, giving clearer expectations as to the likely impact and timing of a future pandemic and the potential effectiveness of mitigation measures. We found that that heterogeneity seen in 2009 can largely be explained by marked differences in school calendars, human mobility and demography across Europe. We also conclude that much of the variation in timing of the pandemic in Europe would have been predictable on the basis of data available in early June 2009. Our work supports the use of models accounting for the structure of complex modern societies for giving insight to policy makers in future pandemics.

Modeling for Estimating Influenza Patients from ILI Surveillance Data in Korea

Objective

Prediction of influenza incidence among outpatients from an influenza surveillance system is important for public influenza strategy.

Methods

We developed two influenza prediction models through influenza surveillance data of the Korea Centers for Disease Control and Prevention (each year, each province and metropolitan city; total reported patients with influenza-like illness stratified by age) for 6 years from 2005 to 2010 and disease-specific data (influenza code J09-J11, monthly number of influenza patients, total number of outpatients and hospital visits) from the Health Insurance Review and Assessment service.

Results

Incidence of influenza in each area, year, and month was estimated from our prediction models, which were validated by simulation processes. For example, in November 2009, Seoul and Joenbuk, the final number of influenza patients calculated by prediction models A and B underestimated actual reported cases by 64 and 833 patients, respectively, in Seoul and 6 and 9 patients, respectively, in Joenbuk. R-square demonstrated that prediction model A was more suitable than model B for estimating the number of influenza patients.

Conclusion

Our prediction models from the influenza surveillance system could estimate the nationwide incidence of influenza. This prediction will provide important basic data for national quarantine activities and distributing medical resources in future pandemics.

The effect of risk perception on the 2009 H1N1 pandemic influenza dynamics

Background

The 2009 H1N1 pandemic influenza dynamics in Italy was characterized by a notable pattern: as it emerged from the analysis of influenza-like illness data, after an initial period (September–mid-October 2009) characterized by a slow exponential increase in the weekly incidence, a sudden and sharp increase of the growth rate was observed by mid-October. The aim here is to understand whether spontaneous behavioral changes in the population could be responsible for such a pattern of epidemic spread.

Methodology/Principal Findings

In order to face this issue, a mathematical model of influenza transmission, accounting for spontaneous behavioral changes driven by cost/benefit considerations on the perceived risk of infection, is proposed and validated against empirical epidemiological data. The performed investigation revealed that an initial overestimation of the risk of infection in the general population, possibly induced by the high concern for the emergence of a new influenza pandemic, results in a pattern of spread compliant with the observed one. This finding is also supported by the analysis of antiviral drugs purchase over the epidemic period. Moreover, by assuming a generation time of 2.5 days, the initially diffuse misperception of the risk of infection led to a relatively low value of the reproductive number , which increased to in the subsequent phase of the pandemic.

Conclusions/Significance

This study highlights that spontaneous behavioral changes in the population, not accounted by the large majority of influenza transmission models, can not be neglected to correctly inform public health decisions. In fact, individual choices can drastically affect the epidemic spread, by altering timing, dynamics and overall number of cases.

Comparing large-scale computational approaches to epidemic modeling: agent-based versus structured metapopulation models

Background

In recent years large-scale computational models for the realistic simulation of epidemic outbreaks have been used with increased frequency. Methodologies adapt to the scale of interest and range from very detailed agent-based models to spatially-structured metapopulation models. One major issue thus concerns to what extent the geotemporal spreading pattern found by different modeling approaches may differ and depend on the different approximations and assumptions used.

Methods

We provide for the first time a side-by-side comparison of the results obtained with a stochastic agent-based model and a structured metapopulation stochastic model for the progression of a baseline pandemic event in Italy, a large and geographically heterogeneous European country. The agent-based model is based on the explicit representation of the Italian population through highly detailed data on the socio-demographic structure. The metapopulation simulations use the GLobal Epidemic and Mobility (GLEaM) model, based on high-resolution census data worldwide, and integrating airline travel flow data with short-range human mobility patterns at the global scale. The model also considers age structure data for Italy. GLEaM and the agent-based models are synchronized in their initial conditions by using the same disease parameterization, and by defining the same importation of infected cases from international travels.

Results

The results obtained show that both models provide epidemic patterns that are in very good agreement at the granularity levels accessible by both approaches, with differences in peak timing on the order of a few days. The relative difference of the epidemic size depends on the basic reproductive ratio, R₀, and on the fact that the metapopulation model consistently yields a larger incidence than the agent-based model, as expected due to the differences in the structure in the intra-population contact pattern of the approaches. The age breakdown analysis shows that similar attack rates are obtained for the younger age classes.

Conclusions

The good agreement between the two modeling approaches is very important for defining the tradeoff between data availability and the information provided by the models. The results we present define the possibility of hybrid models combining the agent-based and the metapopulation approaches according to the available data and computational resources.