Research article | Open | Open Peer Review | Published:
Estimation of Zika virus prevalence by appearance of microcephaly
BMC Infectious Diseasesvolume 16, Article number: 754 (2016)
There currently is a severe Zika Virus (ZIKV) epidemic in Brazil and other South American countries. Due to international travel, this poses severe public health risk of ZIKV importation to other countries. We estimate the prevalence of ZIKV in an import region by the time a microcephaly case is detected, since microcephaly is presently the most significant indication of ZIKV presence.
We establish a mathematical model to describe ZIKV spread from a source region to an import region. This model incorporates both vector transmission (between humans and mosquitoes) and sexual transmission (from males to females). We take account of population structure through a contact network for sexually active individuals. Parameter values of our model are either taken from the literature or estimated from travel data.
This model gives us the probability distribution of time until detection of the first microcephaly case. Based on current field observations, our results also indicate that the percentage of infected pregnant women that results in fetal abnormalities is more likely to be on the smaller end of the 1%–30% spectrum that is currently hypothesized. Our model predicts that for import regions with at least 250,000 people, on average 1,000–12,000 will have been infected by the time of the first detection of microcephaly, and on average 200–1,500 will be infectious at this time. Larger population sizes do not significantly change our predictions.
By the first detection of a microcephaly case, a sizable fraction of the population will have been infected by ZIKV. It is thus clear that adequate surveillance, isolation, and quarantine are needed in susceptible import regions to stop the dissemination of a Zika epidemic.
In March 2015, an outbreak of Zika virus (hereafter referred to as ZIKV) was discovered in Bahia, Brazil . Recent findings  attribute the introduction of ZIKV in Brazil to a single imported infection, estimated to be during the latter half of 2013. This outbreak is currently ongoing with public health efforts to mitigate the dispersal of ZIKV throughout other populations. Recent ZIKV outbreaks before the current outbreak in Brazil include outbreaks in Pacific Islands, such as in 2007 in Yap Island, Micronesia  and in 2013 in French Polynesia . Currently, there are two main strains of ZIKV, namely an African lineage and an Asian lineage .
Acute ZIKV infections in humans may result in fever, macopapular rash, conjunctivitis, muscle pain and joint aches . These symptoms of ZIKV infections are usually mild and difficult to differentiate from other viral infections thus resulting in vast underreporting of ZIKV cases. However, epidemiological studies of recent outbreaks of ZIKV indicate associations between ZIKV and neurological disorders, especially microcephaly and significant brain anomalies in fetuses of infected pregnant females . This has been confirmed during both the French Polynesia outbreak  and the current outbreak in Brazil ; see  and  for specific cases of fetal abnormalities during infected pregnancies. In addition, both the infection by the Brazilian strain  and the infection by the SZ01 Asian strain  have been demonstrated to cause microcephaly in mice fetuses.
For a population that is not under active surveillance, an abnormal fetus is very likely the first indication of a ZIKV epidemic in the population. However, starting from the importation of ZIKV into a population, it may take some time for ZIKV to infect a pregnant female due to the low fraction of such females in the population (approximately 1% in the United States ). Moreover, only infections in their first trimester are more likely to result in an increased risk for fetal abnormality such as microcephaly , and such an abnormal fetus can be detected at roughly 20 weeks . Thus, by the time a microcephaly case is detected, ZIKV could have spread in the population for a long time, and caused an epidemic. Estimates of the final size of ZIKV epidemics on Yap Island and in French Polynesia are 73% and 66%, respectively [3, 8], showing that a significant proportion of individuals in a susceptible population become infected with ZIKV during an outbreak. In conjunction with the association of ZIKV with neurological disorders, these percentages indicate potential repercussions on a large scale.
In this paper, we quantify the prevalence and total infections of ZIKV in a region (e.g., a city, a state or a country) at the time of the first detection of microcephaly, estimates that are crucial for public health planning at the start of a ZIKV outbreak.
We develop a mathematical model for the spread of ZIKV, which is transmitted to humans by mosquitoes of the Aedes genus; for e.g., A. aegypti, A. albopictus and A. africanus mosquitoes . These mosquitoes also transmit yellow fever, dengue and chikungunya. Mathematical models for these vector transmitted diseases have been extensively studied [16, 17]. However, ZIKV can also be transmitted through sexual contact from an infected male to a susceptible female [15, 18–20]. In fact, ZIKV persists in semen much longer than in blood and other bodily fluids . Thus vector transmission may have a shorter infectious period than sexual transmission. Recovery after infection seems to confer long-term immunity for individuals . There currently exist no cure, no treatment, and no vaccine for ZIKV. This effectively renders every individual that has not contracted ZIKV a susceptible one, and indicates a large potential for the dissemination of ZIKV through human populations.
Our mathematical model considers two regions to incorporate the importation of cases from a “source” region (Region 1) experiencing an epidemic into a completely susceptible “import” region (Region 2). Figure 1 illustrates our model for the transmission of ZIKV in two regions, and Table 1 summarizes parameter values and their ranges; see Additional file 1: for justification of the model and of parameter values.
The time to the first infection T p of a pregnant female in Region 2 can be estimated from this model. New cases per day in pregnant females that result in microcephaly can be modeled as a non-homogeneous Poisson process with rate
λ(t) = − p × z × (rate of change of susceptible sexually active females in Region 2),
where p is the proportion of pregnant females in the sexually active female population and z is the proportion of pregnancies that are affected by ZIKV and result in fetal abnormalities. The probability density function for T p is λ(t)exp(∫ o t λ(τ)dτ). For example, 1% of the US population are pregnant females , which corresponds to p = 3% since 1/3 of the population is sexually active females. It is estimated that z = 29% in Brazil , whereas it is estimated that z = 1% based on data from the ZIKV outbreak in French Polynesia .
The time to the first detection of a microcephaly case T d can then be computed by adding the delay from the infection of a pregnant female to the detection of abnormality in the fetus. Fetal abnormality most likely results from ZIKV infection in the first trimester, as shown in the 2013-2014 outbreak of ZIKV in French Polynesia . Because the infection of the mother can occur any time during her pregnancy, by the detection at 20 weeks, i.e., 140 days , T d = T p + 140 − T f where T f is a uniformly distributed random infection time between 0 and 84 days (i.e., during the first trimester).
Numerically simulating our model to the time of the first microcephaly detection T d yields the prevalence of ZIKV in the import region at that time. The simulations are done using baseline parameter values as in Table 1, and with the random network distributions for sexually active males and females in the same region being Poisson with mean 2.25, and for males and females of different regions being Poisson with mean 0.1. These values are chosen so that the resulting basic reproduction number R0 and final size agree with literature values; see the Results section. We assume that the population of Region 1 is 2 million and that initially one person in Region 1 is infectious, while all other people and all mosquitoes in Regions 1 and 2 are susceptible.
The uncertainty in parameter values in Table 1 is incorporated in the simulations using the Latin Hypercube Sampling (LHS) maximin criteria , assuming all parameters except the death rate of mosquitoes d have uniform distributions with ranges given in Table 1, and 1/d is normally distributed with mean μ = 50 days and standard deviation σ = 3 days. For estimating T p and T d distributions, we sample 2,000 points from the LH parameter space. For prevalence and total case estimates in the import region, we sample 100 points in the LH parameter space, and for each point selected, we compute the probability density function of the time until detection of ZIKV and sample a further 100 points from this distribution. We select LH data that gives rise to ZIKV epidemics with a least 1 pregnant female in a population of 2 million.
The basic reproduction number R0 is computed for the baseline parameters in Table 1: R0 ≈ 1.4 (TA, Additional file 1), agreeing with other estimated values . For these parameters, the final sizes of the outbreaks in Regions 1 and 2 are numerically calculated to be 63.95% and 64.66%, respectively. These values are near the reported final sizes of outbreaks on Yap Island in 2007 and French Polynesia in 2013-2014.
For a range of z values, Fig. 2a shows box plots of the length of time until a pregnant female contracts ZIKV and passes it on to her fetus for population sizes in Region 2 ranging from 0.25 to 2 million. This median time lies between 525 and 800 days, depending on population size and z values. Figure 2b shows the box plots of the time to the first microcephaly detection for the same settings as used in Fig. 2a, giving a median time lying between 625 and 900 days. For every doubling of population size, both median times go down by at most 30 days.
Based on our sensitivity analysis procedure, we present in Fig. 3a boxplots of the 10,000 data points of prevalence (number of infectious individuals) at the time of detection T d in the import region (Region 2) for fixed z values, and in Fig. 3b boxplots of 10,000 data points of the total cases by T d in the import region for fixed z values. In Fig. 3a and b the populations of the import region are fixed at 0.25, 0.5, 1, and 2 million, whereas that of the source region remains at 2 million. Figure 3 shows that by the time a microcephaly case is detected in an import region with a population of at least 250,000, the median prevalence of ZIKV ranges from 200-1,500, and the median total infected ranges from 1,000-12,000. Larger population sizes do not significantly change these predictions.
In order to quantify ZIKV transmission from a source region to an import region, we formulated a mathematical model for the transmission of ZIKV in two regions. As in other ZIKV models [24, 25], we incorporate in our model vector transmission (between mosquitoes and humans), and sexual transmission (from sexually active males to females). We take account of population structure through a contact network for sexually active individuals, which is more realistic than the homogeneous mixing assumption of other models [24, 25]. Using data from the literature and estimated from travel, we obtain the probability distributions of the time to the detection of microcephaly T d in the import region for various population sizes. Based on this, we compute the prevalence at detection of microcephaly in an import region. Our method for estimating T d can be adapted to other models of diseases that cause microcephaly or other neurological disorders in a fetus, for example, rubella.
Our values of T d presented in Fig. 3b agree with the recent identification of microcephaly in a region in Colombia around March 2016  as it is likely that a single infectious individual was responsible for the spread of ZIKV from Asia to Brazil in the second half of 2013 , corresponding to a time delay of about 800 days from the arrival of ZIKV in the source region to the detection of fetal abnormalities in a fetus in the import region in Colombia. Based on the time until microcephaly was observed in Columbia, our results on T d in Fig. 2a indicate it is more likely that the percentage of infected pregnancies is between 1% and 10%.
From our model estimates, upon detection of microcephaly in a region with a population of at least 250,000 people, on average 1,000–12,000 will have been infected, and on average 200–1,500 will be infectious at this time. Thus, it is imperative that potential import regions monitor for ZIKV infections at all times. If indeed 1% ≤ z ≤ 10%, total cases and prevalences are more likely to be on the higher end of the reported ranges, indicating even greater local dissemination of ZIKV, and a greater public health issue.
ZIKV may be also linked to Guillain-Barré Syndrome (GBS) and other neurological syndromes . However, GBS may arise from many other causes . These syndromes may have higher prevalences and shorter delays of onset than microcephaly. Thus our model suggests that monitoring for these syndromes seems to be beneficial for early detection of ZIKV, provided the link between ZIKV and GBS is established and quantified. Currently, detecting ZIKV by GBS cases may be less reliable than by microcephaly.
Through simulations of our model, we have shown that by the first detection of a microcephaly case, or similar fetal abnormalities, a sizable fraction (between 0.4% and 4.8%) of the population will have been infected by ZIKV. This makes controlling of local ZIKV spread very difficult. It is thus important for public health agencies to adopt surveillance (and possibly quarantine and isolation) measures to travelers from a region with an ongoing epidemic to avoid severe threats of ZIKV importation.
Technical Appendix (Additional file 1)
Campos G, Bandeira A, Sardi S. Zika virus outbreak, Bahia, Brazil. Emerg Infect Dis. 2015;21:1885–6.
Faria N, Azevedo R, Kraemer M, Souza R, Cunha M, Hill S, et al. Zika virus in the Americas: early epidemiological and genetic findings. Science. 2016;352:345–9.
Duffy M, Chen TH, Hancock W, Powers A, Kool J, Lanciotti R, et al. Zika virus outbreak on Yap Island, Federated States of Micronesia. N Engl J Med. 2009;360:2536–43.
Cao-Lormeau VM, Roche C, Teissier A, Robin E, Berry AL, Mallet HP, et al. Zika virus, French Polynesia, South Pacific, 2013. Emerg Infect Dis. 2014;20:1085–6.
Haddow A, Schuh A, Yasuda C, Kasper M, Heang V, Huy R, et al. Genetic characterization of Zika virus strains: geographic expansions of the Asian lineage. PLoS Negl Trop Dis. 2012;6:e1477.
MacFadden D, Bogoch I. Zika virus infection. Can Med Assoc J. 2016;188:367.
Rasmussen SA, Jamieson DJ, Honein MA, Petersen LR. Zika virus and birth defects --- reviewing the evidence for causality. N Engl J Med. 2007;374:1981–7.
Cauchemez S, Besnard M, Bompard P, Dub T, Guillemette-Artur P, Eyrolle-Guignot D, et al. Association between Zika virus and microcephaly in French Polynesia, 2013-2015: a retrospective study. Lancet. 2016;387:2125–32.
Brasil P, Pereira J, Gabaglia C, Damasceno L, Wakimoto M, Nogueira R, et al. Zika virus infection in pregnant women in Rio de Janeiro-preliminary report. N Engl J Med 2016; 10.1056/NEJMoa1602412:Published online before print.
Mlakar J, Korva M, Tul N, Popovic M, Poljsak-Prijatelj M, Mraz J, et al. Zika virus associated with microcephaly. N Engl J Med. 2016;374:951–8.
Driggers R, Ho CY, Korhonen E, Kuivanen S, Jskelinen A, Smura T, et al. Zika virus infection with prolonged maternal viremia and fetal brain abnormalities. N Engl J Med. 2016;374:2142–51.
Cugola F, Fernandes I, Russo F, Freitas B, Dias J, Guimarães KP, et al. The Brazilian Zika virus strain causes birth defects in experimental models. Nature. 2016;534:267–71.
Li C, Xu D, Ye Q, Hong S, Jiang Y, Liu X, et al. Zika virus disrupts neural progenitor development and leads to microcephaly in mice. Cell Stem Cell. 2016;19:120–6.
Jamieson D, Honein M, Rasmussen S, Williams J, DL Swerdlow MB, Lindstrom S, et al. H1N1 2009 influenza virus infection during pregnancy in the USA. Lancet. 2009;374:451–8.
D’Ortenzio E, Matheron S, de Lamballerie X, Hubert B, Piorkowski G, Maquart M, et al. Evidence of sexual transmission of Zika virus. N Engl J Med. 2016;374:2195–8.
Reiner R, Perkins T, Barker C, Niu T, Chaves L, Ellis A, et al. A systematic review of mathematical models of mosquito-borne pathogen transmission: 1970–2010. J R Soc Interface. 2013;10:20120921.
Manore C, Hickmann K, Xu S, Hyman H. Comparing dengue and chikungunya emergence and endemic transmission in A. aegypti and A. albopictus. J Theor Biol. 2014;356:174–91.
Foy B, Kobylinski K, Foy J, Blitvich B, Rosa A, Haddow A, et al. Probable non-vector-borne transmission of Zika virus, Colorado, USA. Emerg Infect Dis. 2011;17:880–2.
Musso D, Roche C, Robin E, Nhan T, Teissier A, Cao-Lormeau V. Potential sexual transmission of Zika virus. Emerg Infect Dis. 2015;21:359–61.
McCarthy M. Zika virus was transmitted by sexual contact in Texas, health officials report. BMJ. 2016;352:i720.
Atkinson B, Hearn P, Afrough B, Lumley S, Carter D, Aarons E, et al. Detection of Zika virus in semen. Emerg Infect Dis. 2016;22:940.
Christofferson R. Zika virus emergence and expansion: lessons learned from dengue and chikungunya may not provide all the answers. Am J Trop Med Hyg. 2016;95:15–8.
Blower S, Dowlatabadi H. Sensitivity and uncertainty analysis for complex models of disease transmission: an HIV model, as an example. Int Stat Rev. 1994;62:229–43.
Gao D, Lou Y, He D, Porco TC, Kuang Y, Chowell G, et al. Prevention and control of Zika as a mosquito-borne and sexually transmitted disease: A mathematical modeling analysis. Sci Rep 2016;10.1038/srep28070:Published online before print.
Towers S, Brauer F, Castillo-Chavez C, Falconar A, Mubayi A, Romero-Vivas C. Estimation of the reproduction number of the 2015 Zika virus outbreak in Barranquilla, Colombia, and a first estimate of the relative role of sexual transmission. Epidemics. 2016;17:50–5.
WHO. Zika virus microcephaly and Guillain-Barré syndrome situation report 14 April 2016. Technical report, World Health Organization, 2016.
Ferguson N, Cucunubá Z, Dorigatti I, Nedjati-Gilani G, Donnelly C, Basáñez MG, et al. Countering Zika in Latin America. Science 2016;10.1126/science.aag0219:Published online before print.
PAHO. Zika-Epidemiological Update 14 April 2016. Technical report, Panamerican Health Organization, 2016.
Riker D, Vila-Goulding J. The boom in Brazilians traveling to the United States. J Int Commerce Econ. 2013;5:1–15.
Petersen E, Polen K, Meaney-Delman D, Ellington S, Oduyebo T, Cohn A, et al. Update: interim guidance for health care providers caring for women of reproductive age with possible Zika virus exposure-United States, 2016. MMWR Morb Mortal Wkly Rep. 2016;65:315–22.
McMeniman C, Lane R, Cass B, Fong A, Sidhu M, Wang YF, et al. Stable introduction of a life-shortening Wolbachia infection into the mosquito Aedes aegypti. Science. 2009;323:141–4.
The authors thank two reviewers for their comments and suggestions, which have improved our exposition.
This research is supported by NSERC, through NSERC Discovery Grants (J.M. and P.vdD.) and NSERC-USRA (C.M.S.-R.).
Availability of data and materials
Technical Appendix online.
All three authors conceived and discussed the study as well as drafted the manuscript. CSR carried out the numerical simulations under the guidance of JM. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent to publish
Ethics approval and consent to participate
We predict that by the detection of the first case of microcephaly in a region, a significant number of people will have been infected with Zika virus, posing a significant public health threat.
This file contains the main methods that are involved in the model formulation and the estimation of parameters. This file also contains the technical details relating to the computation of the basic reproduction number R 0 . (DOCX 139 kb)