Assessing the role of multiple mechanisms increasing the age of dengue cases in Thailand.

The mean age of dengue hemorrhagic fever (DHF) cases increased considerably in Thailand from 8.1 to 24.3 y between 1981 and 2017 (mean annual increase of 0.45 y). Alternative proposed explanations for this trend, such as changes in surveillance practices, reduced mosquito­human contact, and shifts in population demographics, have different implications for global dengue epidemiology. To evaluate the contribution of each of these hypothesized mechanisms to the observed data, we developed 20 nested epidemiological models of dengue virus infection, allowing for variation over time in population demographics, infection hazards, and reporting rates. We also quantified the effect of removing or retaining each source of variation in simulations of the age trajectory. Shifts in the age structure of susceptibility explained 58% of the observed change in age. Adding heterogeneous reporting by age and reductions in per-serotype infection hazard to models with shifts in susceptibility explained an additional 42%. Reductions in infection hazards were mostly driven by changes in the number of infectious individuals at any time (another consequence of shifting age demographics) rather than changes in the transmissibility of individual infections. We conclude that the demographic transition drives the overwhelming majority of the observed change as it changes both the age structure of susceptibility and the number of infectious individuals. With the projected Thai population age structure, our results suggest a continuing increase in age of DHF cases, shifting the burden toward individuals with more comorbidity. These insights into dengue epidemiology may be relevant to many regions of the globe currently undergoing comparable changes in population demographics.

Effects of long-term exposure to outdoor air pollution on COVID-19 incidence: A population-based cohort study accounting for SARS-CoV-2 exposure levels in the Netherlands.

Several studies have linked air pollution to COVID-19 morbidity and severity. However, these studies do not account for exposure levels to SARS-CoV-2, nor for different sources of air pollution. We analyzed individual-level data for 8.3 million adults in the Netherlands to assess associations between long-term exposure to ambient air pollution and SARS-CoV-2 infection (i.e., positive test) and COVID-19 hospitalisation risks, accounting for spatiotemporal variation in SARS-CoV-2 exposure levels during the first two major epidemic waves (February 2020-February 2021). We estimated average annual concentrations of PM10, PM2.5 and NO2 at residential addresses, overall and by PM source (road traffic, industry, livestock, other agricultural sources, foreign sources, other Dutch sources), at 1 × 1 km resolution, and weekly SARS-CoV-2 exposure at municipal level. Using generalized additive models, we performed interval-censored survival analyses to assess associations between individuals' average exposure to PM10, PM2.5 and NO2 in the three years before the pandemic (2017-2019) and COVID-19-outcomes, adjusting for SARS-CoV-2 exposure, individual and area-specific confounders. In single-pollutant models, per interquartile (IQR) increase in exposure, PM10 was associated with 7% increased infection risk and 16% increased hospitalisation risk, PM2.5 with 8% increased infection risk and 18% increased hospitalisation risk, and NO2 with 3% increased infection risk and 11% increased hospitalisation risk. Bi-pollutant models suggested that effects were mainly driven by PM. Associations for PM were confirmed when stratifying by urbanization degree, epidemic wave and testing policy. All emission sources of PM, except industry, showed adverse effects on both outcomes. Livestock showed the most detrimental effects per unit exposure, whereas road traffic affected severity (hospitalisation) more than infection risk. This study shows that long-term exposure to air pollution increases both SARS-CoV-2 infection and COVID-19 hospitalisation risks, even after controlling for SARS-CoV-2 exposure levels, and that PM may have differential effects on these COVID-19 outcomes depending on the emission source.

Can incorporating genotyping data into efficacy estimators improve efficiency of early phase malaria vaccine trials?

BACKGROUND: Early phase malaria vaccine field trials typically measure malaria infection by PCR or thick blood smear microscopy performed on serially sampled blood. Vaccine efficacy (VE) is the proportion reduction in an endpoint due to vaccination and is often calculated as VEHR = 1-hazard ratio or VERR = 1-risk ratio. Genotyping information can distinguish different clones and distinguish multiple infections over time, potentially increasing statistical power. This paper investigates two alternative VE endpoints incorporating genotyping information: VEmolFOI, the vaccine-induced proportion reduction in incidence of new clones acquired over time, and VEC, the vaccine-induced proportion reduction in mean number of infecting clones per exposure. METHODS: Power of VEmolFOI and VEC was compared to that of VEHR and VERR by simulations and analytic derivations, and the four VE methods were applied to three data sets: a Phase 3 trial of RTS,S malaria vaccine in 6912 African infants, a Phase 2 trial of PfSPZ Vaccine in 80 Burkina Faso adults, and a trial comparing Plasmodium vivax incidence in 466 Papua New Guinean children after receiving chloroquine + artemether lumefantrine with or without primaquine (as these VE methods can also quantify effects of other prevention measures). By destroying hibernating liver-stage P. vivax, primaquine reduces subsequent reactivations after treatment completion. RESULTS: In the trial of RTS,S vaccine, a significantly reduced number of clones at first infection was observed, but this was not the case in trials of PfSPZ Vaccine or primaquine, although the PfSPZ trial lacked power to show a reduction. Resampling smaller data sets from the large RTS,S trial to simulate phase 2 trials showed modest power gains from VEC compared to VEHR for data like those from RTS,S, but VEC is less powerful than VEHR for trials in which the number of clones at first infection is not reduced. VEmolFOI was most powerful in model-based simulations, but only the primaquine trial collected enough serial samples to precisely estimate VEmolFOI. The primaquine VEmolFOI estimate decreased after most control arm liver-stage infections reactivated (which mathematically resembles a waning vaccine), preventing VEmolFOI from improving power. CONCLUSIONS: The power gain from the genotyping methods depends on the context. Because input parameters for early phase power calculations are often uncertain, these estimators are not recommended as primary endpoints for small trials unless supported by targeted data analysis. TRIAL REGISTRATIONS: NCT00866619, NCT02663700, NCT02143934.

Irrigation-Induced Environmental Changes Sustain Malaria Transmission and Compromise Intervention Effectiveness.

BACKGROUND: Irrigated agriculture enhances food security, but it potentially promotes mosquito-borne disease transmission and affects vector intervention effectiveness. This study was conducted in the irrigated and nonirrigated areas of rural Homa Bay and Kisumu Counties, Kenya. METHODS: We performed cross-sectional and longitudinal surveys to determine Plasmodium infection prevalence, clinical malaria incidence, molecular force of infection (molFOI), and multiplicity of infection. We examined the impact of irrigation on the effectiveness of the new interventions. RESULTS: We found that irrigation was associated with >2-fold higher Plasmodium infection prevalence and 3-fold higher clinical malaria incidence compared to the nonirrigated area. Residents in the irrigated area experienced persistent, low-density parasite infections and higher molFOI. Addition of indoor residual spraying was effective in reducing malaria burden, but the reduction was more pronounced in the nonirrigated area than in the irrigated area. CONCLUSIONS: Our findings collectively suggest that irrigation may sustain and enhance Plasmodium transmission and affects intervention effectiveness.

Individual, Household, and Community Drivers of Dengue Virus Infection Risk in Kamphaeng Phet Province, Thailand.

BACKGROUND: Dengue virus (DENV) often circulates endemically. In such settings with high levels of transmission, it remains unclear whether there are risk factors that alter individual infection risk. METHODS: We tested blood taken from individuals living in multigenerational households in Kamphaeng Phet province, Thailand for DENV antibodies (N = 2364, mean age 31 years). Seropositivity ranged from 45.4% among those 1-5 years old to 99.5% for those >30 years. Using spatially explicit catalytic models, we estimated that 11.8% of the susceptible population gets infected annually. RESULTS: We found that 37.5% of the variance in seropositivity was explained by unmeasured household-level effects with only 4.2% explained by spatial differences between households. The serostatus of individuals from the same household remained significantly correlated even when separated by up to 15 years in age. CONCLUSIONS: These findings show that despite highly endemic transmission, persistent differences in infection risk exist across households, the reasons for which remain unclear.

Host preferences inhibit transmission from potential superspreader host species.

Host species that are particularly abundant, infectious and/or infected tend to contribute disproportionately to symbiont (parasite or mutualist) maintenance in multi-host systems. Therefore, in a facultative multi-host system where two host species had high densities, high symbiont infestation intensities and high infestation prevalence, we expected interspecific transmission rates to be high. Instead, we found that interspecific symbiont transmission rates to caged sentinel hosts were an order of magnitude lower than intraspecific transmission rates in the wild. Using laboratory experiments to decompose transmission rates, we found that opportunities for interspecific transmission were frequent, where interspecific and intraspecific contact rate functions were statistically indistinguishable. However, most interspecific contacts did not lead to transmission events owing to a previously unrecognized transmission barrier: strong host preferences. During laboratory choice experiments, the symbiont preferred staying on or dispersing to its current host species, even though the oligochaete symbiont is a globally distributed host generalist that can survive and reproduce on many snail host species. These surprising results suggest that when managing symbiont transmission, identifying key host species is still important, but it may be equally important to identify and manage transmission barriers that keep potential superspreader host species in check.

Spatiotemporal variations in exposure: Chagas disease in Colombia as a case study.

Age-stratified serosurvey data are often used to understand spatiotemporal trends in disease incidence and exposure through estimating the Force-of-Infection (FoI). Typically, median or mean FoI estimates are used as the response variable in predictive models, often overlooking the uncertainty in estimated FoI values when fitting models and evaluating their predictive ability. To assess how this uncertainty impact predictions, we compared three approaches with three levels of uncertainty integration. We propose a performance indicator to assess how predictions reflect initial uncertainty.In Colombia, 76 serosurveys (1980-2014) conducted at municipality level provided age-stratified Chagas disease prevalence data. The yearly FoI was estimated at the serosurvey level using a time-varying catalytic model. Environmental, demographic and entomological predictors were used to fit and predict the FoI at municipality level from 1980 to 2010 across Colombia.A stratified bootstrap method was used to fit the models without temporal autocorrelation at the serosurvey level. The predictive ability of each model was evaluated to select the best-fit models within urban, rural and (Amerindian) indigenous settings. Model averaging, with the 10 best-fit models identified, was used to generate predictions.Our analysis shows a risk of overconfidence in model predictions when median estimates of FoI alone are used to fit and evaluate models, failing to account for uncertainty in FoI estimates. Our proposed methodology fully propagates uncertainty in the estimated FoI onto the generated predictions, providing realistic assessments of both central tendency and current uncertainty surrounding exposure to Chagas disease.

Estimating the Force of Infection for Dengue Virus Using Repeated Serosurveys, Ouagadougou, Burkina Faso.

Because of limited data on dengue virus in Burkina Faso, we conducted 4 consecutive age-stratified longitudinal serologic surveys, ≈6 months apart, among persons 1-55 years of age, during June 2015-March 2017, which included a 2016 outbreak. The seroconversion rate before the serosurvey enrollment was estimated by binomial regression, taking age as the duration of exposure, and assuming constant force of infection (FOI) over age and calendar time. We calculated FOI between consecutive surveys and rate ratios for potentially associated characteristics based on seroconversion using the duration of intervals. Among 2,897 persons at enrollment, 66.3% were IgG-positive, and estimated annual FOI was 5.95%. Of 1,269 enrollees participating in all 4 serosurveys, 438 were IgG-negative at enrollment. The annualized FOI ranged from 10% to 20% (during the 2016 outbreak). Overall, we observed high FOI for dengue. These results could support decision-making about control and preventive measures for dengue.

Parasite evolution in an age-structured population.

Although mortality increases with age in most organisms, senescence is missing from models of parasite evolution. Since virulence evolves according to the host's mortality, and since virulence influences the intensity of transmission, which determines the average age at infection and thus the mortality rate of a senescing host, we expected that epi-evolutionary feedbacks would underlie the evolution of virulence in a population of senescing hosts. We tested this idea by extending an age-structured model of epidemiological dynamics with the parasite's evolution. A straightforward prediction of our model is that stronger senescence forces the evolution of higher virulence. However, the model also reveals that the evolved virulence depends on the average age at infection, giving an evolutionary feedback with the epidemiological situation, a prediction not found when assuming a constant mortality rate with age. Additionally, and in contrast to most models of parasite evolution, we found that the virulence at the evolutionary equilibrium is influenced by whether the force of infection depends on the density or on the frequency of infected hosts, due to changes in the average age at infection. Our findings suggest that ignoring age-specific effects, and in particular senescence, can give misleading predictions about parasite evolution.

Dynamics and determinants of the force of infection of dengue virus from 1994 to 2015 in Managua, Nicaragua.

Dengue virus (DENV) is the most prevalent human vector-borne viral disease. The force of infection (FoI), the rate at which susceptible individuals are infected in a population, is an important metric for infectious disease modeling. Understanding how and why the FoI of DENV changes over time is critical for developing immunization and vector control policies. We used age-stratified seroprevalence data from 12 years of the Pediatric Dengue Cohort Study in Nicaragua to estimate the annual FoI of DENV from 1994 to 2015. Seroprevalence data revealed a change in the rate at which children acquire DENV-specific immunity: in 2004, 50% of children age >4 years were seropositive, but by 2015, 50% seropositivity was reached only by age 11 years. We estimated a spike in the FoI in 1997-1998 and 1998-1999 and a gradual decline thereafter, and children age <4 years experienced a lower FoI. Two hypotheses to explain the change in the FoI were tested: (i) a transition from introduction of specific DENV serotypes to their endemic transmission and (ii) a population demographic transition due to declining birth rates and increasing life expectancy. We used mathematical models to simulate these hypotheses. We show that the initial high FoI can be explained by the introduction of DENV-3 in 1994-1998, and that the overall gradual decline in the FoI can be attributed to demographic shifts. Changes in immunity and demographics strongly impacted DENV transmission in Nicaragua. Population-level measures of transmission intensity are dynamic and thus challenging to use to guide vaccine implementation locally and globally.

Epidemic Dynamics and Adaptive Vaccination Strategy: Renewal Equation Approach.

We use analytical methods to investigate a continuous vaccination strategy's effects on the infectious disease dynamics in a closed population and a demographically open population. The methodology and key assumptions are based on Breda et al. (J Biol Dyn 6(Sup2):103-117, 2012). We show that the cumulative force of infection for the closed population and the endemic force of infection in the demographically open population can be reduced significantly by combining two factors: the vaccine effectiveness and the vaccination rate. The impact of these factors on the force of infection can transform an endemic steady state into a disease-free state.

Force of Infection and True Infection Rate of Dengue in Singapore: Implications for Dengue Control and Management.

National data on dengue notifications do not capture all dengue infections and do not reflect the true intensity of disease transmission. To assess the true dengue infection rate and disease control efforts in Singapore, we conducted age-stratified serosurveys among residents after a 2013 outbreak that was the largest dengue outbreak on record. The age-weighted prevalence of dengue immunoglobulin G among residents was 49.8% (95% confidence interval: 48.4, 51.1) in 2013 and 48.6% (95% confidence interval: 47.0, 50.0) in 2017; prevalence increased with age. Combining these data with those from previous serosurveys, the year-on-year estimates of the dengue force of infection from 1930 to 2017 revealed a significant decrease from the late 1960s to the mid-1990s, after which the force of infection remained stable at approximately 10 per 1,000 persons per year. The reproduction number (R0) had also declined since the 1960s. The reduction in dengue transmission may be attributed to the sustained national vector program and partly to a change in the age structure of the population. The improved estimated ratio of notified cases to true infections, from 1:14 in 2005-2009 to 1:6 in 2014-2017, signifies that the national notification system, which relies on diagnosed cases, has improved over time. The data also suggest that the magnitudes of dengue epidemics cannot be fairly compared across calendar years and that the current disease control program remains applicable.

Climate drives spatial variation in Zika epidemics in Latin America.

Between 2015 and 2017, Zika virus spread rapidly through populations in the Americas with no prior exposure to the disease. Although climate is a known determinant of many Aedes-transmitted diseases, it is currently unclear whether climate was a major driver of the Zika epidemic and how climate might have differentially impacted outbreak intensity across locations within Latin America. Here, we estimated force of infection for Zika over time and across provinces in Latin America using a time-varying susceptible-infectious-recovered model. Climate factors explained less than 5% of the variation in weekly transmission intensity in a spatio-temporal model of force of infection by province over time, suggesting that week to week transmission within provinces may be too stochastic to predict. By contrast, climate and population factors were highly predictive of spatial variation in the presence and intensity of Zika transmission among provinces, with pseudo-R2 values between 0.33 and 0.60. Temperature, temperature range, rainfall and population size were the most important predictors of where Zika transmission occurred, while rainfall, relative humidity and a nonlinear effect of temperature were the best predictors of Zika intensity and burden. Surprisingly, force of infection was greatest in locations with temperatures near 24°C, much lower than previous estimates from mechanistic models, potentially suggesting that existing vector control programmes and/or prior exposure to other mosquito-borne diseases may have limited transmission in locations most suitable for Aedes aegypti, the main vector of Zika, dengue and chikungunya viruses in Latin America.

Disentangling serology to elucidate henipa- and filovirus transmission in Madagascar fruit bats.

Bats are reservoirs for emerging human pathogens, including Hendra and Nipah henipaviruses and Ebola and Marburg filoviruses. These viruses demonstrate predictable patterns in seasonality and age structure across multiple systems; previous work suggests that they may circulate in Madagascar's endemic fruit bats, which are widely consumed as human food. We aimed to (a) document the extent of henipa- and filovirus exposure among Malagasy fruit bats, (b) explore seasonality in seroprevalence and serostatus in these bat populations and (c) compare mechanistic hypotheses for possible transmission dynamics underlying these data. To this end, we amassed and analysed a unique dataset documenting longitudinal serological henipa- and filovirus dynamics in three Madagascar fruit bat species. We uncovered serological evidence of exposure to Hendra-/Nipah-related henipaviruses in Eidolon dupreanum, Pteropus rufus and Rousettus madagascariensis, to Cedar-related henipaviruses in E. dupreanum and R. madagascariensis and to Ebola-related filoviruses in P. rufus and R. madagascariensis. We demonstrated significant seasonality in population-level seroprevalence and individual serostatus for multiple viruses across these species, linked to the female reproductive calendar. An age-structured subset of the data highlighted evidence of waning maternal antibodies in neonates, increasing seroprevalence in young and decreasing seroprevalence late in life. Comparison of mechanistic epidemiological models fit to these data offered support for transmission hypotheses permitting waning antibodies but retained immunity in adult-age bats. Our findings suggest that bats may seasonally modulate mechanisms of pathogen control, with consequences for population-level transmission. Additionally, we narrow the field of candidate transmission hypotheses by which bats are presumed to host and transmit potentially zoonotic viruses globally.

Force of infection of Middle East respiratory syndrome in dromedary camels in Kenya.

Middle East respiratory syndrome coronavirus (MERS-CoV) is a zoonotic disease transmitted from dromedary camels to people, which can result in outbreaks with human-to-human transmission. Because it is a subclinical infection in camels, epidemiological measures other than prevalence are challenging to assess. This study estimated the force of infection (FOI) of MERS-CoV in camel populations from age-stratified serological data. A cross-sectional study of MERS-CoV was conducted in Kenya from July 2016 to July 2017. Seroprevalence was stratified into four age groups: <1, 1-2, 2-3 and >3 years old. Age-independent and age-dependent linear and quadratic generalised linear models were used to estimate FOI in pastoral and ranching camel herds. Models were compared based on computed AIC values. Among pastoral herds, the age-dependent quadratic FOI was the best fit model, while the age-independent FOI was the best fit for the ranching herd data. FOI provides an indirect estimate of infection risk, which is especially valuable where direct estimates of incidence and other measures of infection are challenging to obtain. The FOIs estimated in this study provide important insight about MERS-CoV dynamics in the reservoir species, and contribute to our understanding of the zoonotic risks of this important public health threat.

Host size and proximity to diseased neighbours drive the spread of a coral disease outbreak in Hawai'i.

Understanding how disease risk varies over time and across heterogeneous populations is critical for managing disease outbreaks, but this information is rarely known for wildlife diseases. Here, we demonstrate that variation in host and pathogen factors drive the direction, duration and intensity of a coral disease outbreak. We collected longitudinal health data for 200 coral colonies, and found that disease risk increased with host size and severity of diseased neighbours, and disease spread was highest among individuals between 5 and 20 m apart. Disease risk increased by 2% with every 10 cm increase in host size. Healthy colonies with severely diseased neighbours (greater than 75% affected tissue) were 1.6 times more likely to develop disease signs compared with colonies with moderately diseased neighbours (25-75% affected tissue). Force of infection ranged from 7 to 20 disease cases per 1000 colonies (mean = 15 cases per 1000 colonies). The effective reproductive ratio, or average number of secondary infections per infectious individual, ranged from 0.16 to 1.22. Probability of transmission depended strongly on proximity to diseased neighbours, which demonstrates that marine disease spread can be highly constrained within patch reefs.

Force of infection of Helicobacter pylori in Mexico: evidence from a national survey using a hierarchical Bayesian model.

Helicobacter pylori (H. pylori) is present in the stomach of half of the world's population. The force of infection describes the rate at which susceptibles acquire infection. In this article, we estimated the age-specific force of infection of H. pylori in Mexico. Data came from a national H. pylori seroepidemiology survey collected in Mexico in 1987-88. We modelled the number of individuals with H. pylori at a given age as a binomial random variable. We assumed that the cumulative risk of infection by a given age follows a modified exponential catalytic model, allowing some fraction of the population to remain uninfected. The cumulative risk of infection was modelled for each state in Mexico and were shrunk towards the overall national cumulative risk curve using Bayesian hierarchical models. The proportion of the population that can be infected (i.e. susceptible population) is 85.9% (95% credible interval (CR) 84.3%-87.5%). The constant rate of infection per year of age among the susceptible population is 0.092 (95% CR 0.084-0.100). The estimated force of infection was highest at birth 0.079 (95% CR 0.071-0.087) decreasing to zero as age increases. This Bayesian hierarchical model allows stable estimation of state-specific force of infection by pooling information between the states, resulting in more realistic estimates.

Inferring infection hazard in wildlife populations by linking data across individual and population scales.

Our ability to infer unobservable disease-dynamic processes such as force of infection (infection hazard for susceptible hosts) has transformed our understanding of disease transmission mechanisms and capacity to predict disease dynamics. Conventional methods for inferring FOI estimate a time-averaged value and are based on population-level processes. Because many pathogens exhibit epidemic cycling and FOI is the result of processes acting across the scales of individuals and populations, a flexible framework that extends to epidemic dynamics and links within-host processes to FOI is needed. Specifically, within-host antibody kinetics in wildlife hosts can be short-lived and produce patterns that are repeatable across individuals, suggesting individual-level antibody concentrations could be used to infer time since infection and hence FOI. Using simulations and case studies (influenza A in lesser snow geese and Yersinia pestis in coyotes), we argue that with careful experimental and surveillance design, the population-level FOI signal can be recovered from individual-level antibody kinetics, despite substantial individual-level variation. In addition to improving inference, the cross-scale quantitative antibody approach we describe can reveal insights into drivers of individual-based variation in disease response, and the role of poorly understood processes such as secondary infections, in population-level dynamics of disease.

Contact and contagion: Probability of transmission given contact varies with demographic state in bighorn sheep.

Understanding both contact and probability of transmission given contact are key to managing wildlife disease. However, wildlife disease research tends to focus on contact heterogeneity, in part because the probability of transmission given contact is notoriously difficult to measure. Here, we present a first step towards empirically investigating the probability of transmission given contact in free-ranging wildlife. We used measured contact networks to test whether bighorn sheep demographic states vary systematically in infectiousness or susceptibility to Mycoplasma ovipneumoniae, an agent responsible for bighorn sheep pneumonia. We built covariates using contact network metrics, demographic information and infection status, and used logistic regression to relate those covariates to lamb survival. The covariate set contained degree, a classic network metric describing node centrality, but also included covariates breaking the network metrics into subsets that differentiated between contacts with yearlings, ewes with lambs, and ewes without lambs, and animals with and without active infections. Yearlings, ewes with lambs, and ewes without lambs showed similar group membership patterns, but direct interactions involving touch occurred at a rate two orders of magnitude higher between lambs and reproductive ewes than between any classes of adults or yearlings, and one order of magnitude higher than direct interactions between multiple lambs. Although yearlings and non-reproductive bighorn ewes regularly carried M. ovipneumoniae, our models suggest that a contact with an infected reproductive ewe had approximately five times the odds of producing a lamb mortality event of an identical contact with an infected dry ewe or yearling. Consequently, management actions targeting infected animals might lead to unnecessary removal of young animals that carry pathogens but rarely transmit. This analysis demonstrates a simple logistic regression approach for testing a priori hypotheses about variation in the odds of transmission given contact for free-ranging hosts, and may be broadly applicable for investigations in wildlife disease ecology.

On the origin and timing of Zika virus introduction in Brazil.

The timing and origin of Zika virus (ZIKV) introduction in Brazil has been the subject of controversy. Initially, it was assumed that the virus was introduced during the FIFA World Cup in June-July 2014. Then, it was speculated that ZIKV may have been introduced by athletes from French Polynesia (FP) who competed in a canoe race in Rio de Janeiro in August 2014. We attempted to apply mathematical models to determine the most likely time window of ZIKV introduction in Brazil. Given that the timing and origin of ZIKV introduction in Brazil may be a politically sensitive issue, its determination (or the provision of a plausible hypothesis) may help to prevent undeserved blame. We used a simple mathematical model to estimate the force of infection and the corresponding individual probability of being infected with ZIKV in FP. Taking into account the air travel volume from FP to Brazil between October 2013 and March 2014, we estimated the expected number of infected travellers arriving at Brazilian airports during that period. During the period between December 2013 and February 2014, 51 individuals travelled from FP airports to 11 Brazilian cities. Basing on the calculated force of ZIKV infection (the per capita rate of new infections per time unit) and risk of infection (probability of at least one new infection), we estimated that 18 (95% CI 12-22) individuals who arrived in seven of the evaluated cities were infected. When basic ZIKV reproduction numbers greater than one were assumed in the seven evaluated cities, ZIKV could have been introduced in any one of the cities. Based on the force of infection in FP, basic reproduction ZIKV number in selected Brazilian cities, and estimated travel volume, we concluded that ZIKV was most likely introduced and established in Brazil by infected travellers arriving from FP in the period between October 2013 and March 2014, which was prior to the two aforementioned sporting events.