Paolo Bajardi

Seasonal transmission potential and activity peaks of the new influenza A(H1N1): a Monte Carlo likelihood analysis based on human mobility

BackgroundOn 11 June the World Health Organization officially raised the phase of pandemic alert (with regard to the new H1N1 influenza strain) to level 6. As of 19 July, 137,232 cases of the H1N1 influenza strain have been officially confirmed in 142 different countries, and the pandemic unfolding in the Southern hemisphere is now under scrutiny to gain insights about the next winter wave in the Northern hemisphere. A major challenge is pre-empted by the need to estimate the transmission potential of the virus and to assess its dependence on seasonality aspects in order to be able to use numerical models capable of projecting the spatiotemporal pattern of the pandemic.MethodsIn the present work, we use a global structured metapopulation model integrating mobility and transportation data worldwide. The model considers data on 3,362 subpopulations in 220 different countries and individual mobility across them. The model generates stochastic realizations of the epidemic evolution worldwide considering 6 billion individuals, from which we can gather information such as prevalence, morbidity, number of secondary cases and number and date of imported cases for each subpopulation, all with a time resolution of 1 day. In order to estimate the transmission potential and the relevant model parameters we used the data on the chronology of the 2009 novel influenza A(H1N1). The method is based on the maximum likelihood analysis of the arrival time distribution generated by the model in 12 countries seeded by Mexico by using 1 million computationally simulated epidemics. An extended chronology including 93 countries worldwide seeded before 18 June was used to ascertain the seasonality effects.ResultsWe found the best estimate R0 = 1.75 (95% confidence interval (CI) 1.64 to 1.88) for the basic reproductive number. Correlation analysis allows the selection of the most probable seasonal behavior based on the observed pattern, leading to the identification of plausible scenarios for the future unfolding of the pandemic and the estimate of pandemic activity peaks in the different hemispheres. We provide estimates for the number of hospitalizations and the attack rate for the next wave as well as an extensive sensitivity analysis on the disease parameter values. We also studied the effect of systematic therapeutic use of antiviral drugs on the epidemic timeline.ConclusionThe analysis shows the potential for an early epidemic peak occurring in October/November in the Northern hemisphere, likely before large-scale vaccination campaigns could be carried out. The baseline results refer to a worst-case scenario in which additional mitigation policies are not considered. We suggest that the planning of additional mitigation policies such as systematic antiviral treatments might be the key to delay the activity peak in order to restore the effectiveness of the vaccination programs.

Human Mobility Networks, Travel Restrictions, and the Global Spread of 2009 H1N1 Pandemic

After the emergence of the H1N1 influenza in 2009, some countries responded with travel-related controls during the early stage of the outbreak in an attempt to contain or slow down its international spread. These controls along with self-imposed travel limitations contributed to a decline of about 40% in international air traffic to/from Mexico following the international alert. However, no containment was achieved by such restrictions and the virus was able to reach pandemic proportions in a short time. When gauging the value and efficacy of mobility and travel restrictions it is crucial to rely on epidemic models that integrate the wide range of features characterizing human mobility and the many options available to public health organizations for responding to a pandemic. Here we present a comprehensive computational and theoretical study of the role of travel restrictions in halting and delaying pandemics by using a model that explicitly integrates air travel and short-range mobility data with high-resolution demographic data across the world and that is validated by the accumulation of data from the 2009 H1N1 pandemic. We explore alternative scenarios for the 2009 H1N1 pandemic by assessing the potential impact of mobility restrictions that vary with respect to their magnitude and their position in the pandemic timeline. We provide a quantitative discussion of the delay obtained by different mobility restrictions and the likelihood of containing outbreaks of infectious diseases at their source, confirming the limited value and feasibility of international travel restrictions. These results are rationalized in the theoretical framework characterizing the invasion dynamics of the epidemics at the metapopulation level.

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

BackgroundMathematical 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.MethodsWe 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.ResultsReal-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.ConclusionsOur 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.

Optimizing surveillance for livestock disease spreading through animal movements

The spatial propagation of many livestock infectious diseases critically depends on the animal movements among premises; so the knowledge of movement data may help us to detect, manage and control an outbreak. The identification of robust spreading features of the system is however hampered by the temporal dimension characterizing population interactions through movements. Traditional centrality measures do not provide relevant information as results strongly fluctuate in time and outbreak properties heavily depend on geotemporal initial conditions. By focusing on the case study of cattle displacements in Italy, we aim at characterizing livestock epidemics in terms of robust features useful for planning and control, to deal with temporal fluctuations, sensitivity to initial conditions and missing information during an outbreak. Through spatial disease simulations, we detect spreading paths that are stable across different initial conditions, allowing the clustering of the seeds and reducing the epidemic variability. Paths also allow us to identify premises, called sentinels, having a large probability of being infected and providing critical information on the outbreak origin, as encoded in the clusters. This novel procedure provides a general framework that can be applied to specific diseases, for aiding risk assessment analysis and informing the design of optimal surveillance systems.

Association between Recruitment Methods and Attrition in Internet-Based Studies

Internet-based systems for epidemiological studies have advantages over traditional approaches as they can potentially recruit and monitor a wider range of individuals in a relatively inexpensive fashion. We studied the association between communication strategies used for recruitment (offline, online, face-to-face) and follow-up participation in nine Internet-based cohorts: the Influenzanet network of platforms for influenza surveillance which includes seven cohorts in seven different European countries, the Italian birth cohort Ninfea and the New Zealand birth cohort ELF. Follow-up participation varied from 43% to 89% depending on the cohort. Although there were heterogeneities among studies, participants who became aware of the study through an online communication campaign compared with those through traditional offline media seemed to have a lower follow-up participation in 8 out of 9 cohorts. There were no clear differences in participation between participants enrolled face-to-face and those enrolled through other offline strategies. An Internet-based campaign for Internet-based epidemiological studies seems to be less effective than an offline one in enrolling volunteers who keep participating in follow-up questionnaires. This suggests that even for Internet-based epidemiological studies an offline enrollment campaign would be helpful in order to achieve a higher participation proportion and limit the cohort attrition.

Determinants of follow-up participation in the Internet-based European influenza surveillance platform Influenzanet.

Background “Influenzanet” is a network of Internet-based platforms aimed at collecting real-time data for influenza surveillance in several European countries. More than 30,000 European volunteers participate every year in the study, representing one of the largest existing Internet-based multicenter cohorts. Each week during the influenza season, participants are asked to report their symptoms (if any) along with a set of additional questions. Objective Focusing on the first influenza season of 2011-12, when the Influenzanet system was completely harmonized within a common framework in Sweden, the United Kingdom, the Netherlands, Belgium, France, Italy, and Portugal, we investigated the propensity of users to regularly come back to the platform to provide information about their health status. Our purpose was to investigate demographic and behavioral factors associated with participation in follow-up. Methods By means of a multilevel analysis, we evaluated the association between regular participation during the season and sociodemographic and behavioral characteristics as measured by a background questionnaire completed by participants on registration. Results We found that lower participation in follow-up was associated with lower educational status (odds ratio [OR] 0.80, 95% CI 0.75-0.85), smoking (OR 0.64, 95% CI 0.59-0.70), younger age (OR ranging from 0.30, 95% CI 0.26-0.33 to 0.70, 95% CI 0.64-0.77), not being vaccinated against seasonal influenza (OR 0.77, 95% CI 0.72-0.84), and living in a household with children (OR 0.69, 95% CI 0.65-0.74). Most of these results hold when single countries are analyzed separately. Conclusions Given the opportunistic enrollment of self-selected volunteers in the Influenzanet study, we have investigated how sociodemographic and behavioral characteristics may be associated with follow-up participation in the Influenzanet cohort. The study described in this paper shows that, overall, the most important determinants of participation are related to education and lifestyle: smoking, lower education level, younger age, people living with children, and people who have not been vaccinated against seasonal influenza tend to have a lower participation in follow-up. Despite the cross-country variation, the main findings are similar in the different national cohorts, and indeed the results are found to be valid also when performing a single-country analysis. Differences between countries do not seem to play a crucial role in determining the factors associated with participation in follow-up.

Modeling vaccination campaigns and the Fall/Winter 2009 activity of the new A(H1N1) influenza in the Northern Hemisphere

The unfolding of pandemic influenza A(H1N1) for Fall 2009 in the Northern Hemisphere is still uncertain. Plans for vaccination campaigns and vaccine trials are underway, with the first batches expected to be available early October. Several studies point to the possibility of an anticipated pandemic peak that could undermine the effectiveness of vaccination strategies. Here, we use a structured global epidemic and mobility metapopulation model to assess the effectiveness of massive vaccination campaigns for the Fall/Winter 2009. Mitigation effects are explored depending on the interplay between the predicted pandemic evolution and the expected delivery of vaccines. The model is calibrated using recent estimates on the transmissibility of the new A(H1N1) influenza. Results show that if additional intervention strategies were not used to delay the time of pandemic peak, vaccination may not be able to considerably reduce the cumulative number of cases, even when the mass vaccination campaign is started as early as mid-October. Prioritized vaccination would be crucial in slowing down the pandemic evolution and reducing its burden.

Unveiling patterns of international communities in a global city using mobile phone data

We analyse a large mobile phone activity dataset provided by Telecom Italia for the Telecom Big Data Challenge contest. The dataset reports the international country codes of every call/SMS made and received by mobile phone users in Milan, Italy, between November and December 2013, with a spatial resolution of about 200 meters. We first show that the observed spatial distribution of international codes well matches the distribution of international communities reported by official statistics, confirming the value of mobile phone data for demographic research. Next, we define an entropy function to measure the heterogeneity of the international phone activity in space and time. By comparing the entropy function to empirical data, we show that it can be used to identify the city’s hotspots, defined by the presence of points of interests. Eventually, we use the entropy function to characterize the spatial distribution of international communities in the city. Adopting a topological data analysis approach, we find that international mobile phone users exhibit some robust clustering patterns that correlate with basic socio-economic variables. Our results suggest that mobile phone records can be used in conjunction with topological data analysis tools to study the geography of migrant communities in a global city.

Interplay of network dynamics and heterogeneity of ties on spreading dynamics.

The structure of a network dramatically affects the spreading phenomena unfolding upon it. The contact distribution of the nodes has long been recognized as the key ingredient in influencing the outbreak events. However, limited knowledge is currently available on the role of the weight of the edges on the persistence of a pathogen. At the same time, recent works showed a strong influence of temporal network dynamics on disease spreading. In this work we provide an analytical understanding, corroborated by numerical simulations, about the conditions for infected stable state in weighted networks. In particular, we reveal the role of heterogeneity of edge weights and of the dynamic assignment of weights on the ties in the network in driving the spread of the epidemic. In this context we show that when weights are dynamically assigned to ties in the network, a heterogeneous distribution is able to hamper the diffusion of the disease, contrary to what happens when weights are fixed in time.

Pattern of Tick Aggregation on Mice: Larger Than Expected Distribution Tail Enhances the Spread of Tick-Borne Pathogens

The spread of tick-borne pathogens represents an important threat to human and animal health in many parts of Eurasia. Here, we analysed a 9-year time series of Ixodes ricinus ticks feeding on Apodemus flavicollis mice (main reservoir-competent host for tick-borne encephalitis, TBE) sampled in Trentino (Northern Italy). The tail of the distribution of the number of ticks per host was fitted by three theoretical distributions: Negative Binomial (NB), Poisson-LogNormal (PoiLN), and Power-Law (PL). The fit with theoretical distributions indicated that the tail of the tick infestation pattern on mice is better described by the PL distribution. Moreover, we found that the tail of the distribution significantly changes with seasonal variations in host abundance. In order to investigate the effect of different tails of tick distribution on the invasion of a non-systemically transmitted pathogen, we simulated the transmission of a TBE-like virus between susceptible and infective ticks using a stochastic model. Model simulations indicated different outcomes of disease spreading when considering different distribution laws of ticks among hosts. Specifically, we found that the epidemic threshold and the prevalence equilibria obtained in epidemiological simulations with PL distribution are a good approximation of those observed in simulations feed by the empirical distribution. Moreover, we also found that the epidemic threshold for disease invasion was lower when considering the seasonal variation of tick aggregation.