Estimate of the reproduction number of the 2015 Zika virus outbreak in Barranquilla, Colombia, and estimation of the relative role of sexual transmission

130258-Thumbnail Image.png
Description

Background
In 2015, the Zika arbovirus (ZIKV) began circulating in the Americas, rapidly expanding its global geographic range in explosive outbreaks. Unusual among mosquito-borne diseases, ZIKV has been shown to also be sexually transmitted, although sustained autochthonous transmission due to sexual

Background
In 2015, the Zika arbovirus (ZIKV) began circulating in the Americas, rapidly expanding its global geographic range in explosive outbreaks. Unusual among mosquito-borne diseases, ZIKV has been shown to also be sexually transmitted, although sustained autochthonous transmission due to sexual transmission alone has not been observed, indicating the reproduction number (R0) for sexual transmission alone is less than 1. Critical to the assessment of outbreak risk, estimation of the potential attack rates, and assessment of control measures, are estimates of the basic reproduction number, R0.
Methods
We estimated the R0 of the 2015 ZIKV outbreak in Barranquilla, Colombia, through an analysis of the exponential rise in clinically identified ZIKV cases (n = 359 to the end of November, 2015).
Findings
The rate of exponential rise in cases was ρ = 0.076 days[superscript −1], with 95% CI [0.066,0.087] days[superscript −1]. We used a vector-borne disease model with additional direct transmission to estimate the R0; assuming the R0 of sexual transmission alone is less than 1, we estimated the total R0 = 3.8 [2.4,5.6], and that the fraction of cases due to sexual transmission was 0.23 [0.01,0.47] with 95% confidence.
Interpretation
This is among the first estimates of R0 for a ZIKV outbreak in the Americas, and also among the first quantifications of the relative impact of sexual transmission.

Date Created
2016-10-17

The Timing and Targeting of Treatment in Influenza Pandemics Influences the Emergence of Resistance in Structured Populations

127996-Thumbnail Image.png
Description

Antiviral resistance in influenza is rampant and has the possibility of causing major morbidity and mortality. Previous models have identified treatment regimes to minimize total infections and keep resistance low. However, the bulk of these studies have ignored stochasticity and

Antiviral resistance in influenza is rampant and has the possibility of causing major morbidity and mortality. Previous models have identified treatment regimes to minimize total infections and keep resistance low. However, the bulk of these studies have ignored stochasticity and heterogeneous contact structures. Here we develop a network model of influenza transmission with treatment and resistance, and present both standard mean-field approximations as well as simulated dynamics. We find differences in the final epidemic sizes for identical transmission parameters (bistability) leading to different optimal treatment timing depending on the number initially infected. We also find, contrary to previous results, that treatment targeted by number of contacts per individual (node degree) gives rise to more resistance at lower levels of treatment than non-targeted treatment. Finally we highlight important differences between the two methods of analysis (mean-field versus stochastic simulations), and show where traditional mean-field approximations fail. Our results have important implications not only for the timing and distribution of influenza chemotherapy, but also for mathematical epidemiological modeling in general. Antiviral resistance in influenza may carry large consequences for pandemic mitigation efforts, and models ignoring contact heterogeneity and stochasticity may provide misleading policy recommendations.

Date Created
2013-02-07

A Data-Driven Mathematical Model of CA-MRSA Transmission Among Age Groups: Evaluating the Effect of Control Interventions

128055-Thumbnail Image.png
Description

Community associated methicillin-resistant Staphylococcus aureus (CA-MRSA) has become a major cause of skin and soft tissue infections (SSTIs) in the US. We developed an age-structured compartmental model to study the spread of CA-MRSA at the population level and assess the

Community associated methicillin-resistant Staphylococcus aureus (CA-MRSA) has become a major cause of skin and soft tissue infections (SSTIs) in the US. We developed an age-structured compartmental model to study the spread of CA-MRSA at the population level and assess the effect of control intervention strategies. We used Monte-Carlo Markov Chain (MCMC) techniques to parameterize our model using monthly time series data on SSTIs incidence in children (≤19 years) during January 2004 -December 2006 in Maricopa County, Arizona. Our model-based forecast for the period January 2007–December 2008 also provided a good fit to data. We also carried out an uncertainty and sensitivity analysis on the control reproduction number, Rc which we estimated at 1.3 (95% CI [1.2,1.4]) based on the model fit to data. Using our calibrated model, we evaluated the effect of typical intervention strategies namely reducing the contact rate of infected individuals owing to awareness of infection and decolonization strategies targeting symptomatic infected individuals on both and the long-term disease dynamics. We also evaluated the impact of hypothetical decolonization strategies targeting asymptomatic colonized individuals. We found that strategies focused on infected individuals were not capable of achieving disease control when implemented alone or in combination. In contrast, our results suggest that decolonization strategies targeting the pediatric population colonized with CA-MRSA have the potential of achieving disease elimination.

Date Created
2013-11-21

On the dynamics of dengue virus type 2 with residence times and vertical transmission

130280-Thumbnail Image.png
Description
A two-patch mathematical model of Dengue virus type 2 (DENV-2) that accounts for vectors’ vertical transmission and between patches human dispersal is introduced. Dispersal is modelled via a Lagrangian approach. A host-patch residence-times basic reproduction number is derived and conditions

A two-patch mathematical model of Dengue virus type 2 (DENV-2) that accounts for vectors’ vertical transmission and between patches human dispersal is introduced. Dispersal is modelled via a Lagrangian approach. A host-patch residence-times basic reproduction number is derived and conditions under which the disease dies out or persists are established. Analytical and numerical results highlight the role of hosts’ dispersal in mitigating or exacerbating disease dynamics. The framework is used to explore dengue dynamics using, as a starting point, the 2002 outbreak in the state of Colima, Mexico.
Date Created
2016-08-05

Homo-Psychologicus: Reactionary Behavioural Aspects of Epidemics

128259-Thumbnail Image.png
Description

We formulate an in silico model of pathogen avoidance mechanism and investigate its impact on defensive behavioural measures (e.g., spontaneous social exclusions and distancing, crowd avoidance and voluntary vaccination adaptation). In particular, we use SIR(B)S (e.g., susceptible-infected-recovered with additional behavioural

We formulate an in silico model of pathogen avoidance mechanism and investigate its impact on defensive behavioural measures (e.g., spontaneous social exclusions and distancing, crowd avoidance and voluntary vaccination adaptation). In particular, we use SIR(B)S (e.g., susceptible-infected-recovered with additional behavioural component) model to investigate the impact of homo-psychologicus aspects of epidemics. We focus on reactionary behavioural changes, which apply to both social distancing and voluntary vaccination participations. Our analyses reveal complex relationships between spontaneous and uncoordinated behavioural changes, the emergence of its contagion properties, and mitigation of infectious diseases. We find that the presence of effective behavioural changes can impede the persistence of disease. Furthermore, it was found that under perfect effective behavioural change, there are three regions in the response factor (e.g., imitation and/or reactionary) and behavioural scale factor (e.g., global/local) factors ρ–α behavioural space. Mainly, (1) disease is always endemic even in the presence of behavioural change, (2) behavioural-prevalence plasticity is observed and disease can sometimes be eradication, and (3) elimination of endemic disease under permanence of permanent behavioural change is achieved. These results suggest that preventive behavioural changes (e.g., non-pharmaceutical prophylactic measures, social distancing and exclusion, crowd avoidance) are influenced by individual differences in perception of risks and are a salient feature of epidemics. Additionally, these findings indicates that care needs to be taken when considering the effect of adaptive behavioural change in predicting the course of epidemics, and as well as the interpretation and development of the public health measures that account for spontaneous behavioural changes.

Date Created
2015-10-14

A Graphical Approach to a Model of a Neuronal Tree With a Variable Diameter

128263-Thumbnail Image.png
Description

Tree-like structures are ubiquitous in nature. In particular, neuronal axons and dendrites have tree-like geometries that mediate electrical signaling within and between cells. Electrical activity in neuronal trees is typically modeled using coupled cable equations on multi-compartment representations, where each

Tree-like structures are ubiquitous in nature. In particular, neuronal axons and dendrites have tree-like geometries that mediate electrical signaling within and between cells. Electrical activity in neuronal trees is typically modeled using coupled cable equations on multi-compartment representations, where each compartment represents a small segment of the neuronal membrane. The geometry of each compartment is usually defined as a cylinder or, at best, a surface of revolution based on a linear approximation of the radial change in the neurite. The resulting geometry of the model neuron is coarse, with non-smooth or even discontinuous jumps at the boundaries between compartments. We propose a hyperbolic approximation to model the geometry of neurite compartments, a branched, multi-compartment extension, and a simple graphical approach to calculate steady-state solutions of an associated system of coupled cable equations. A simple case of transient solutions is also briefly discussed.

Date Created
2014-07-09

Apis Mellifera Octopamine Receptor 1 (AmOA1) Expression in Antennal Lobe Networks of the Honey Bee (Apis Mellifera) and Fruit Fly (Drosophila Melanogaster)

128321-Thumbnail Image.png
Description

Octopamine (OA) underlies reinforcement during appetitive conditioning in the honey bee and fruit fly, acting via different subtypes of receptors. Recently, antibodies raised against a peptide sequence of one honey bee OA receptor, AmOA1, were used to study the distribution

Octopamine (OA) underlies reinforcement during appetitive conditioning in the honey bee and fruit fly, acting via different subtypes of receptors. Recently, antibodies raised against a peptide sequence of one honey bee OA receptor, AmOA1, were used to study the distribution of these receptors in the honey bee brain (Sinakevitch et al., 2011). These antibodies also recognize an isoform of the AmOA1 ortholog in the fruit fly (OAMB, mushroom body OA receptor). Here we describe in detail the distribution of AmOA1 receptors in different types of neurons in the honey bee and fruit fly antennal lobes. We integrate this information into a detailed anatomical analysis of olfactory receptor neurons (ORNs), uni- and multi-glomerular projection neurons (uPNs, and mPNs) and local interneurons (LNs) in glomeruli of the antennal lobe. These neurons were revealed by dye injection into the antennal nerve, antennal lobe, medial and lateral antenno-protocerbral tracts (m-APT and l-APT), and lateral protocerebral lobe (LPL) by use of labeled cell lines in the fruit fly or by staining with anti-GABA. We found that ORN receptor terminals and uPNs largely do not show immunostaining for AmOA1. About seventeen GABAergic mPNs leave the antennal lobe through the ml-APT and branch into the LPL. Many, but not all, mPNs show staining for AmOA1. AmOA1 receptors are also in glomeruli on GABAergic processes associated with LNs. The data suggest that in both species one important action of OA in the antennal lobe involves modulation of different types of inhibitory neurons via AmOA1 receptors. We integrated this new information into a model of circuitry within glomeruli of the antennal lobes of these species.

Date Created
2013-10-25
Agent

Where Two Are Fighting, the Third Wins: Stronger Selection Facilitates Greater Polymorphism in Traits Conferring Competition-Dispersal Tradeoffs

128738-Thumbnail Image.png
Description

A major conundrum in evolution is that, despite natural selection, polymorphism is still omnipresent in nature: Numerous species exhibit multiple morphs, namely several abundant values of an important trait. Polymorphism is particularly prevalent in asymmetric traits, which are beneficial to

A major conundrum in evolution is that, despite natural selection, polymorphism is still omnipresent in nature: Numerous species exhibit multiple morphs, namely several abundant values of an important trait. Polymorphism is particularly prevalent in asymmetric traits, which are beneficial to their carrier in disruptive competitive interference but at the same time bear disadvantages in other aspects, such as greater mortality or lower fecundity. Here we focus on asymmetric traits in which a better competitor disperses fewer offspring in the absence of competition. We report a general pattern in which polymorphic populations emerge when disruptive selection increases: The stronger the selection, the greater the number of morphs that evolve. This pattern is general and is insensitive to the form of the fitness function. The pattern is somewhat counterintuitive since directional selection is excepted to sharpen the trait distribution and thereby reduce its diversity (but note that similar patterns were suggested in studies that demonstrated increased biodiversity as local selection increases in ecological communities). We explain the underlying mechanism in which stronger selection drives the population towards more competitive values of the trait, which in turn reduces the population density, thereby enabling lesser competitors to stably persist with reduced need to directly compete. Thus, we believe that the pattern is more general and may apply to asymmetric traits more broadly. This robust pattern suggests a comparative, unified explanation to a variety of polymorphic traits in nature.

Date Created
2016-02-04

Urban Economies and Occupation Space: Can They Get “There” From “Here”?

128749-Thumbnail Image.png
Description

Much of the socioeconomic life in the United States occurs in its urban areas. While an urban economy is defined to a large extent by its network of occupational specializations, an examination of this important network is absent from the

Much of the socioeconomic life in the United States occurs in its urban areas. While an urban economy is defined to a large extent by its network of occupational specializations, an examination of this important network is absent from the considerable body of work on the determinants of urban economic performance. Here we develop a structure-based analysis addressing how the network of interdependencies among occupational specializations affects the ease with which urban economies can transform themselves. While most occupational specializations exhibit positive relationships between one another, many exhibit negative ones, and the balance between the two partially explains the productivity of an urban economy. The current set of occupational specializations of an urban economy and its location in the occupation space constrain its future development paths. Important tradeoffs exist between different alternatives for altering an occupational specialization pattern, both at a single occupation and an entire occupational portfolio levels.

Date Created
2013-09-09

Time-Specific Ecologic Niche Models Forecast the Risk of Hemorrhagic Fever With Renal Syndrome in Dongting Lake District, China, 2005-2010

128753-Thumbnail Image.png
Description

Background: Hemorrhagic fever with renal syndrome (HFRS), a rodent-borne infectious disease, is one of the most serious public health threats in China. Increasing our understanding of the spatial and temporal patterns of HFRS infections could guide local prevention and control strategies.

Methodology/Principal

Background: Hemorrhagic fever with renal syndrome (HFRS), a rodent-borne infectious disease, is one of the most serious public health threats in China. Increasing our understanding of the spatial and temporal patterns of HFRS infections could guide local prevention and control strategies.

Methodology/Principal Findings: We employed statistical models to analyze HFRS case data together with environmental data from the Dongting Lake district during 2005–2010. Specifically, time-specific ecologic niche models (ENMs) were used to quantify and identify risk factors associated with HFRS transmission as well as forecast seasonal variation in risk across geographic areas. Results showed that the Maximum Entropy model provided the best predictive ability (AUC = 0.755). Time-specific Maximum Entropy models showed that the potential risk areas of HFRS significantly varied across seasons. High-risk areas were mainly found in the southeastern and southwestern areas of the Dongting Lake district. Our findings based on models focused on the spring and winter seasons showed particularly good performance. The potential risk areas were smaller in March, May and August compared with those identified for June, July and October to December. Both normalized difference vegetation index (NDVI) and land use types were found to be the dominant risk factors.

Conclusions/Significance: Our findings indicate that time-specific ENMs provide a useful tool to forecast the spatial and temporal risk of HFRS.

Date Created
2014-09-03
Agent