Simulation of an SEIR infectious disease model on the dynamic contact network of conference attendees

The spread of infectious diseases crucially depends on the pattern of contacts among individuals. Knowledge of these patterns is thus essential to inform models and computational efforts. Few empirical studies are however available that provide estimates of the number and duration of contacts among social groups. Moreover, their space and time resolution are limited, so that data is not explicit at the person-to-person level, and the dynamical aspect of the contacts is disregarded. Here, we want to assess the role of data-driven dynamic contact patterns among individuals, and in particular of their temporal aspects, in shaping the spread of a simulated epidemic in the population. We consider high resolution data of face-to-face interactions between the attendees of a conference, obtained from the deployment of an infrastructure based on Radio Frequency Identification (RFID) devices that assess mutual face-to-face proximity. The spread of epidemics along these interactions is simulated through an SEIR model, using both the dynamical network of contacts defined by the collected data, and two aggregated versions of such network, in order to assess the role of the data temporal aspects. We show that, on the timescales considered, an aggregated network taking into account the daily duration of contacts is a good approximation to the full resolution network, whereas a homogeneous representation which retains only the topology of the contact network fails in reproducing the size of the epidemic. These results have important implications in understanding the level of detail needed to correctly inform computational models for the study and management of real epidemics

Spatially Resolved Distribution Function and the Medium-Range Order in Metallic Liquid and Glass

The structural description of disordered systems has been a longstanding challenge in physical science. We propose an atomic cluster alignment method to reveal the development of three-dimensional topological ordering in a metallic liquid as it undercools to form a glass. By analyzing molecular dynamic (MD) simulation trajectories of a Cu64.5Zr35.5 alloy, we show that medium-range order (MRO) develops in the liquid as it approaches the glass transition. Specifically, around Cu sites, we observe “Bergman triacontahedron” packing (icosahedron, dodecahedron and icosahedron) that extends out to the fourth shell, forming an interpenetrating backbone network in the glass. The discovery of Bergman-type MRO from our order-mining technique provides unique insights into the topological ordering near the glass transition and the relationship between metallic glasses and quasicrystals

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

On 11 June the World Health Organization officially raised the phase of pandemic alert (with regard to the new H1N1 influenza strain) to level 6. We use a global structured metapopulation model integrating mobility and transportation data worldwide 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 1M computationally simulated epidemics. An extended chronology including 93 countries worldwide seeded before 18 June was used to ascertain the seasonality effects. We found the best estimate R0 = 1.75 (95% 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. The 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. We suggest that the planning of additional mitigation policies such as systematic antiviral treatments might be the key to delay the activity peak inorder to restore the effectiveness of the vaccination programs.Comment: Paper: 29 Pages, 3 Figures and 5 Tables. Supplementary Information: 29 Pages, 5 Figures and 7 Tables. Print version: http://www.biomedcentral.com/1741-7015/7/4

The GLEaMviz computational tool, a publicly available software to explore realistic epidemic spreading scenarios at the global scale

<p>Abstract</p> <p>Background</p> <p>Computational models play an increasingly important role in the assessment and control of public health crises, as demonstrated during the 2009 H1N1 influenza pandemic. Much research has been done in recent years in the development of sophisticated data-driven models for realistic computer-based simulations of infectious disease spreading. However, only a few computational tools are presently available for assessing scenarios, predicting epidemic evolutions, and managing health emergencies that can benefit a broad audience of users including policy makers and health institutions.</p> <p>Results</p> <p>We present "GLEaMviz", a publicly available software system that simulates the spread of emerging human-to-human infectious diseases across the world. The GLEaMviz tool comprises three components: the client application, the proxy middleware, and the simulation engine. The latter two components constitute the GLEaMviz server. The simulation engine leverages on the Global Epidemic and Mobility (GLEaM) framework, a stochastic computational scheme that integrates worldwide high-resolution demographic and mobility data to simulate disease spread on the global scale. The GLEaMviz design aims at maximizing flexibility in defining the disease compartmental model and configuring the simulation scenario; it allows the user to set a variety of parameters including: compartment-specific features, transition values, and environmental effects. The output is a dynamic map and a corresponding set of charts that quantitatively describe the geo-temporal evolution of the disease. The software is designed as a client-server system. The multi-platform client, which can be installed on the user's local machine, is used to set up simulations that will be executed on the server, thus avoiding specific requirements for large computational capabilities on the user side.</p> <p>Conclusions</p> <p>The user-friendly graphical interface of the GLEaMviz tool, along with its high level of detail and the realism of its embedded modeling approach, opens up the platform to simulate realistic epidemic scenarios. These features make the GLEaMviz computational tool a convenient teaching/training tool as well as a first step toward the development of a computational tool aimed at facilitating the use and exploitation of computational models for the policy making and scenario analysis of infectious disease outbreaks.</p

Resilience management during large-scale epidemic outbreaks

Assessing and managing the impact of large-scale epidemics considering only the individual risk and severity of the disease is exceedingly difficult and could be extremely expensive. Economic consequences, infrastructure and service disruption, as well as the recovery speed, are just a few of the many dimensions along which to quantify the effect of an epidemic on society's fabric. Here, we extend the concept of resilience to characterize epidemics in structured populations, by defining the system-wide critical functionality that combines an individual’s risk of getting the disease (disease attack rate) and the disruption to the system’s functionality (human mobility deterioration). By studying both conceptual and data-driven models, we show that the integrated consideration of individual risks and societal disruptions under resilience assessment framework provides an insightful picture of how an epidemic might impact society. In particular, containment interventions intended for a straightforward reduction of the risk may have net negative impact on the system by slowing down the recovery of basic societal functions. The presented study operationalizes the resilience framework, providing a more nuanced and comprehensive approach for optimizing containment schemes and mitigation policies in the case of epidemic outbreaks

Analysis of the effectiveness of interventions used during the 2009 A/H1N1 influenza pandemic

<p>Abstract</p> <p>Background</p> <p>Following the emergence of the A/H1N1 2009 influenza pandemic, public health interventions were activated to lessen its potential impact. Computer modelling and simulation can be used to determine the potential effectiveness of the social distancing and antiviral drug therapy interventions that were used at the early stages of the pandemic, providing guidance to public health policy makers as to intervention strategies in future pandemics involving a highly pathogenic influenza strain.</p> <p>Methods</p> <p>An individual-based model of a real community with a population of approximately 30,000 was used to determine the impact of alternative interventions strategies, including those used in the initial stages of the 2009 pandemic. Different interventions, namely school closure and antiviral strategies, were simulated in isolation and in combination to form different plausible scenarios. We simulated epidemics with reproduction numbers R₀of 1.5, which aligns with estimates in the range 1.4-1.6 determined from the initial outbreak in Mexico.</p> <p>Results</p> <p>School closure of 1 week was determined to have minimal effect on reducing overall illness attack rate. Antiviral drug treatment of 50% of symptomatic cases reduced the attack rate by 6.5%, from an unmitigated rate of 32.5% to 26%. Treatment of diagnosed individuals combined with additional household prophylaxis reduced the final attack rate to 19%. Further extension of prophylaxis to close contacts (in schools and workplaces) further reduced the overall attack rate to 13% and reduced the peak daily illness rate from 120 to 22 per 10,000 individuals. We determined the size of antiviral stockpile required; the ratio of the required number of antiviral courses to population was 13% for the treatment-only strategy, 25% for treatment and household prophylaxis and 40% for treatment, household and extended prophylaxis. Additional simulations suggest that coupling school closure with the antiviral strategies further reduces epidemic impact.</p> <p>Conclusions</p> <p>These results suggest that the aggressive use of antiviral drugs together with extended school closure may substantially slow the rate of influenza epidemic development. These strategies are more rigorous than those actually used during the early stages of the relatively mild 2009 pandemic, and are appropriate for future pandemics that have high morbidity and mortality rates.</p

Search for rare quark-annihilation decays, B --> Ds(*) Phi

We report on searches for B- --> Ds- Phi and B- --> Ds*- Phi. In the context of the Standard Model, these decays are expected to be highly suppressed since they proceed through annihilation of the b and u-bar quarks in the B- meson. Our results are based on 234 million Upsilon(4S) --> B Bbar decays collected with the BABAR detector at SLAC. We find no evidence for these decays, and we set Bayesian 90% confidence level upper limits on the branching fractions BF(B- --> Ds- Phi) Ds*- Phi)<1.2x10^(-5). These results are consistent with Standard Model expectations.Comment: 8 pages, 3 postscript figues, submitted to Phys. Rev. D (Rapid Communications

Rejuvenation of metallic glasses by non-affine thermal strain.

When a spatially uniform temperature change is imposed on a solid with more than one phase, or on a polycrystal of a single, non-cubic phase (showing anisotropic expansion-contraction), the resulting thermal strain is inhomogeneous (non-affine). Thermal cycling induces internal stresses, leading to structural and property changes that are usually deleterious. Glasses are the solids that form on cooling a liquid if crystallization is avoided--they might be considered the ultimate, uniform solids, without the microstructural features and defects associated with polycrystals. Here we explore the effects of cryogenic thermal cycling on glasses, specifically metallic glasses. We show that, contrary to the null effect expected from uniformity, thermal cycling induces rejuvenation, reaching less relaxed states of higher energy. We interpret these findings in the context that the dynamics in liquids become heterogeneous on cooling towards the glass transition, and that there may be consequent heterogeneities in the resulting glasses. For example, the vibrational dynamics of glassy silica at long wavelengths are those of an elastic continuum, but at wavelengths less than approximately three nanometres the vibrational dynamics are similar to those of a polycrystal with anisotropic grains. Thermal cycling of metallic glasses is easily applied, and gives improvements in compressive plasticity. The fact that such effects can be achieved is attributed to intrinsic non-uniformity of the glass structure, giving a non-uniform coefficient of thermal expansion. While metallic glasses may be particularly suitable for thermal cycling, the non-affine nature of strains in glasses in general deserves further study, whether they are induced by applied stresses or by temperature change.This research was supported by the World Premier International Research Center Initiative (WPI), MEXT, Japan, by NSF China and MOST 973 China, and by the Engineering and the Engineering and Physical Sciences Research Council, UK (Materials World Network project). Y.H.S. acknowledges support from a China Scholarship Council (CSC) scholarship.This is the author accepted manuscript. The final version is available from Nature Publishing Group via http://dx.doi.org/10.1038/nature1467

QTL meta-analysis of root traits in Brassica napus under contrasting phosphorus supply in two growth systems

A high-density SNP-based genetic linkage map was constructed and integrated with a previous map in the Tapidor x Ningyou7 (TNDH) Brassica napus population, giving a new map with a total of 2041 molecular markers and an average marker density which increased from 0.39 to 0.97 (0.82 SNP bin) per cM. Root and shoot traits were screened under low and ‘normal’ phosphate (Pi) supply using a ‘pouch and wick’ system, and had been screened previously in an agar based system. The P-efficient parent Ningyou7 had a shorter primary root length (PRL), greater lateral root density (LRD) and a greater shoot biomass than the P-inefficient parent Tapidor under both treatments and growth systems. Quantitative trait loci (QTL) analysis identified a total of 131 QTL, and QTL meta-analysis found four integrated QTL across the growth systems. Integration reduced the confidence interval by ~41%. QTL for root and shoot biomass were co-located on chromosome A3 and for lateral root emergence were co-located on chromosomes A4/C4 and C8/C9. There was a major QTL for LRD on chromosome C9 explaining ~18% of the phenotypic variation. QTL underlying an increased LRD may be a useful breeding target for P uptake efficiency in Brassica

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

<p>Abstract</p> <p>Background</p> <p>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.</p> <p>Methods</p> <p>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.</p> <p>Results</p> <p>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, <it>R</it>₀, 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.</p> <p>Conclusions</p> <p>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.</p