Estimating the subcritical transmissibility of the Zika outbreak in the State of Florida, USA, 2016
© The Author(s). 2016
Received: 30 September 2016
Accepted: 23 October 2016
Published: 9 November 2016
Florida State has reported autochthonous transmission of Zika virus since late July 2016. Here we assessed the transmissibility associated with the outbreak and generated a short-term forecast.
Time-dependent dynamics of imported cases reported in the state of Florida was approximated by a logistic growth equation. We estimated the reproduction number using the renewal equation in order to predict the incidence of local cases arising from both local and imported primary cases. Using a bootstrap method together with the logistic and renewal equations, a short-term forecast of local and imported cases was carried out.
The reproduction number was estimated at 0.16 (95 % Confidence Interval: 0.13, 0.19). Employing the logistic equation to capture a drastic decline in the number of imported cases expected through the course of 2016, together with the low estimate of the local reproduction number in Florida, the expected number of local reported cases was demonstrated to show an evident declining trend for the remainder of 2016.
The risk of local transmission in the state of Florida is predicted to dramatically decline by the end of 2016.
KeywordsPrediction Zika virus Epidemic Mathematical model Basic reproduction number
The state of Florida has reported over 735 travel-related Zika cases since February 2016, becoming the first state in the continental USA to report multiple laboratory-confirmed autochthonous cases of Zika . The state of Florida has not only experienced local transmission of Zika but also dengue and chikungunya viruses, which are transmitted via the common vector Aedes species.
To assess the risk of infection for local residents and travelers, it is critical to understand whether local transmission is sustained as well as forecast the duration and size of the outbreak. In this study, we carried out a risk assessment of Zika transmission in Florida aimed to estimate the extent of local transmission potential of Zika, e.g., are on-going local chains of transmission sustained in Florida? We also generated a short-term forecast of the expected burden of Zika for the remainder of 2016.
In order to keep Florida residents and visitors abreast of the presence of Zika cases in the state, the Florida Department of Health has maintained up-to-date counts of the number of Zika cases diagnosed in this state  while strengthening a robust mosquito-borne illness surveillance system. Our study relies on confirmed cases of Zika virus infection. Before 29 June 2016, except for pregnant women, only those cases exhibiting at least two of the following symptoms: fever, rash, joint pain and red eyes and having an epidemiological link (travel history or sexual contact with travelers from Zika-affected areas or suspected contact with cases) underwent laboratory diagnosis by serology or rRT-PCR. From 29 June, all laboratory-confirmed asymptomatic cases were counted. Non-pregnant cases with recent travel history to an area with widespread Zika virus transmission were classified as either travel-related or non-travel related cases . Hereafter, travel-related cases are referred to as imported cases, while non-travel related cases are referred to as local cases. We analysed the temporal evolution of confirmed Zika cases from 1 May to 23 September 2016. Since epidemiological and diagnostic procedures typically required 7 days , we analysed weekly case counts with week 0 starting on 1 May 2016.
as discussed elsewhere . Assuming that variation in case counts in week t is sufficiently captured by the Poisson distribution, in agreement with the underlying mechanism of the infection process in deterministic models, the maximum likelihood estimate of R is obtained by minimizing the Poisson-distributed likelihood function that uses the right-hand side of (2) for the expected value. A constant R is supported by the negligible impact of herd immunity on the transmission dynamics due to the limited scope of the epidemic in this area thus far.
where K is referred to as the carrying capacity (i.e., the expected total number of cases during the outbreak), γ, the steepness of the curve, and t 0 the time of the sigmoid’s midpoint. In addition to the recent interest in using phenomenological models that generalize the logistic equation , the SIR (susceptible-infectious-recovered) model can be approximated by the logistic curve . Hence, it may not be surprising that eq. (3) can be useful to capture single-epidemic transmission dynamics. Using (3), the expected weekly counts in week t was obtained by E(i t) = I(7 t)-I(7(t −1)) for t > 0. Assuming again that variation in imported case counts in week t is sufficiently captured by the Poisson distribution, maximum likelihood estimates of K, γ and t 0 are obtained by minimizing the likelihood function that uses the right-hand side of (3) for the expected value. Profile likelihood was employed to compute the 95 % confidence intervals (CI).
Lastly, we employ a parametric bootstrap method to resample R, K, γ and t 0 to generate short-term forecasts for case counts from week 21 (week starting with 25 September 2016) to 34 (week starting with 25 December 2016) . The number of imported cases can be simulated using (3) and the number of local cases can be obtained by additionally using (2), accounting for the dependence among the estimated parameters. Posterior 95 % prediction intervals were derived by taking 2.5th and 97.5th percentile points from 1000 bootstrap simulations.
Figure 2b compares observed and predicted local cases of Zika virus infection in Florida. Due to the expected decline in the number of imported cases together with our low and subcritical estimate of the reproduction number that lies substantially below 1.0, the number of local cases are also expected to wane over time, with only a few additional cases expected towards the end of 2016. Our predicted number of additional imported cases from week 21 to 34 was estimated at 41 (95 % PrI: 29, 54) cases. It is worth noting that at the last time of writing, case counts for weeks 22–23, indicated an additional 36 imported and 33 local reported cases, which is in line with our forecast and the predicted decreasing trend in case counts.
In this study we analysed case series of reported cases of Zika virus infection in Florida State, May-September 2016, in order to estimate the local transmission potential and generate forecasts of the number of expected reported cases for the remainder of 2016. The reproduction number R of local Zika virus transmission in Florida was estimated to lie in the range of 0.13–0.19, likely reflecting primarily limited transmission potential and perhaps the effects of intensive vector control efforts including the strong advice of draining standing water, covering clothing and bare skin with repellent, and covering windows with screens. This finding indicates that local transmission chains in this area have not been sustained in the absence of continued importation of infected individuals. Moreover, the number of imported Zika cases has been declining since mid-August. If the declining trend continues for the remainder of 2016, our forecast based on the logistic growth model indicates a drastic decline in the number of imported cases expected by the end of 2016. This finding and our subcritical estimate of the local reproduction number in Florida point to a declining trend in the expected number of local reported cases for the remainder of 2016. Although the risk of local transmission within 2016 will not necessarily decline to zero, as transmission via the sexual route cannot be ruled out, the infection risk will be likely very small as the epidemic approaches the end of 2016. Although we employed a phenomenological model (e.g., did not explicitly account for seasonal mosquito population dynamics), our forecast may even be an overestimate as a significantly reduced risk of mosquito transmission can be expected during the winter season.
Our study is the first to have demonstrated that the effective transmission potential of Zika virus is substantially low in the state of Florida even in an area with a history of sustained transmission of dengue and chikungunya during the summer months. Our estimate of the local reproduction number at 0.16 is far below those estimates previously reported from areas with widespread Zika virus transmissions [9–12]. Reflecting the decline in the number of infections in Central America (e.g. Mexico) and South America (e.g. Colombia), the frequency of imported cases is also expected to decrease over time. At the same time, we have shown that local cases are also expected to decline . With substantially decreased local transmission risk through the remainder of 2016, Florida could provide a basis for model-based transmission analyses and risk assessments of Zika across the world .
This study was conducted for real-time risk assessment of Zika virus infection, involving a number of important limitations. First, the dataset inherently involved reporting delays. It is likely that some symptomatic cases were retrospectively diagnosed by serology dating back to their date of illness onset. Rather than explicitly assessing the precise temporal transmission dynamics, we deemed sensible to generate an estimate of R based on a series of weekly case counts of limited precision. Second, asymptomatic transmission of Zika virus is known to be common , and the growing awareness of Zika virus among physicians is likely affecting differential diagnosis. Thus, observed confirmed cases must have experienced ascertainment bias. Third, the diagnostic criteria for reported cases confirmed with Zika was revised on 29 June 2016 to include asymptomatic individuals. This change was caused by the approval of a revised interim case definition in June 2016, which was originally issued in February 2016, along with the revision in their reporting system that partly relies on commercial laboratory testing. The ascertainment bias is likely to have been greatly reduced afterwards. Rather than attempting to provide a detailed analysis of the transmission dynamics of this outbreak using mechanistic models, we have intended to capture the overall temporal course of imported cases using a phenomenological approach. Fourth, our model did not explicitly account for the seasonality of local transmission. For this reason, our forecast for expected cases through the remainder of 2016 is likely an overestimate for a conservative assessment of the transmission risk associated with the outbreak.
While the revised version of this study was prepared in mid-October, a confirmed case was reported from Miami-Dade County. Despite possible transient increase in the number of cases, the expected decreasing trend of importation and limited transmission potential in Florida remains unchanged. Despite these limitations, we strongly believe that the present study sheds light on two critical issues, i.e. (i) transmission of Zika virus has not been locally sustained in Florida without continued external forcing and (ii) given a near-future reduction in imported cases, it is not farfetched to expect to dramatic decline in the number of local cases in Florida.
A case series of reported cases of Zika virus infection in Florida State was analysed, to estimate the local transmission potential and generate forecasts of the number of expected reported cases for the remainder of 2016. The reproduction number R of local Zika virus transmission in Florida was estimated to lie in the range of 0.13-0.19, indicating that local transmission chains cannot be sustained in the absence of continued importation of infected individuals. Moreover, the number of imported Zika cases has been declining since mid-August. The expected number of local reported cases is very likely to show an evident declining trend for the remainder of 2016.
Real-time reverse transcriptase polymerase chain reaction
HN received funding support from the Japan Agency for Medical Research and Development, the Japan Society for the Promotion of Science (JSPS) KAKENHI Grant Numbers 16KT0130, 16 K15356 and 26700028, Program for Advancing Strategic International Networks to Accelerate the Circulation of Talented Researchers, the Japan Science and Technology Agency (JST) CREST program and RISTEX program for Science of Science, Technology and Innovation Policy. GC acknowledges financial support from the NSF grant 1414374 as part of the joint NSF-NIH-USDA Ecology and Evolution of Infectious Diseases program; UK Biotechnology and Biological Sciences Research Council grant BB/M008894/1, NSF-IIS RAPID award #1518939, and NSF grant 1318788 III: Small: Data Management for Real-Time Data Driven Epidemic simulation, and the Division of International Epidemiology and Population Studies, The Fogarty International Center, US National Institutes of Health. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
Collected datasheet is available from the authors (LD) upon request.
LD collected and arranged epidemiological data. All authors jointly formulated mathematical modelling approach. HN and LD analysed the data. HN and GC drafted the manuscript and all authors revised the manuscript. All other authors gave comments on the revised manuscript and approved the final version of the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Centers for Disease Control and Prevention (USA). Florida investigation links four recent Zika cases to local mosquito-borne virus transmission [Internet]. Accessed 7 Nov 2016. Available from: http://www.cdc.gov/media/releases/2016/p0729-florida-zika-cases.html
- Florida Department of Health (USA). Department of Health Daily Zika Update, 26 September 2016 [Internet]. Accessed 7 Nov 2016. Available from: http://www.floridahealth.gov/newsroom/2016/09/092616-zika-update.html
- Florida Department of Health (USA). Department of Health’s Zika testing, investigation and notification process [Internet]. Accessed 7 Nov 2016. Available from: http://www.floridahealth.gov/diseases-and-conditions/zika-virus/_documents/zika-testing-process-final.pdf
- Majumder MS, Cohn E, Fish D, Brownstein JS. Estimating a feasible serial interval range for Zika fever. Bull WHO. 2016;in press. doi:10.2471/BLT.16.170860
- Roberts MG, Nishiura H. Early estimation of the reproduction number in the presence of imported cases: pandemic influenza H1N1-2009 in New Zealand. PLoS One. 2011;6:e17835.View ArticlePubMedPubMed CentralGoogle Scholar
- Chowell G, Viboud C, Simonsen L, Moghadas SM. Characterizing the reproduction number of epidemics with early subexponential growth dynamics. J R Soc Interface. 2016;13(123):20160659.
- Omori R, Nishiura H. Theoretical basis to measure the impact of short-lasting control of an infectious disease on the epidemic peak. Theor Biol Med Model. 2011;8:2.View ArticlePubMedPubMed CentralGoogle Scholar
- Efron B, Tibshirani RJ. An Introduction to the bootstrap. New York: Chapman & Hall; 1993.View ArticleGoogle Scholar
- Lessler J, Chaisson LH, Kucirka LM, Bi Q, Grantz K, Salje H, et al. Assessing the global threat from Zika virus. Science. 2016;353(6300):aaf8160.View ArticlePubMedGoogle Scholar
- Rojas DP, Dean NE, Yang Y, Kenah E, Quintero J, Tomasi S, et al. The epidemiology and transmissibility of Zika virus in Girardot and San Andres island, Colombia, September 2015 to January 2016. Euro Surveill. 2016;21(28). doi:10.2807/1560-7917.ES.2016.21.28.30283.
- Nishiura H, Kinoshita R, Mizumoto K, Yasuda Y, Nah K. Transmission potential of Zika virus infection in the South Pacific. Int J Infect Dis. 2016;45:95–7.View ArticlePubMedGoogle Scholar
- Chowell G, Hincapie-Palacio D, Ospina J, Pell B, Tariq A, Dahal S, Moghadas S, Smirnova A, Simonsen L, Viboud C. Using Phenomenological Models to Characterize Transmissibility and Forecast Patterns and Final Burden of Zika Epidemics. PLoS Curr. 2016;8. doi:10.1371/currents.outbreaks.f14b2217c902f453d9320a43a35b9583.
- Ferguson NM, Cucunubá ZM, Dorigatti I, Nedjati-Gilani GL, Donnelly CA, Basáñez MG, et al. Countering the Zika epidemic in Latin America. Science. 2016;353(6297):353–4. doi:10.1126/science.aag0219.View ArticlePubMedGoogle Scholar
- Robert MA, Christofferson RC, Silva NJ, Vasquez C, Mores CN, Wearing HJ. Modeling Mosquito-Borne Disease Spread in U.S. Urbanized Areas: The Case of Dengue in Miami. PLoS One. 2016;11(8):e0161365. doi:10.1371/journal.pone.0161365.View ArticlePubMedPubMed CentralGoogle Scholar
- Duffy MR, Chen TH, Hancock WT, Powers AM, Kool JL, Lanciotti RS, et al. Zika virus outbreak on Yap Island, Federated States of Micronesia. The New England Journal of Medicine. 2009;360(24):2536–43.View ArticlePubMedGoogle Scholar