Osong Public Health Res PerspectOsong Public Health and Research Perspectives2210-90992233-6052241594433766904S2210-9099(11)00002-610.1016/j.phrp.2011.04.001Invited Original ArticleIncubation Period of Ebola Hemorrhagic Virus Subtype ZaireEichnerMartinmartin.eichner@uni-tuebingen.deaDowellScott F.aFireseNinabDepartment of Medical Biometry, University of Tübingen, Tübingen, GermanyCenters for Disease Control and Prevention, Atlanta, Georgia, USACorresponding author. martin.eichner@uni-tuebingen.de124201112420116201121372822011252011© 2011 Published by Elsevier B.V. on behalf of Korea Centers for Disease Control and Prevention.2011This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.Objectives

Ebola hemorrhagic fever has killed over 1300 people, mostly in equatorial Africa. There is still uncertainty about the natural reservoir of the virus and about some of the factors involved in disease transmission. Until now, a maximum incubation period of 21 days has been assumed.

Methods

We analyzed data collected during the Ebola outbreak (subtype Zaire) in Kikwit, Democratic Republic of the Congo, in 1995 using maximum likelihood inference and assuming a log-normally distributed incubation period.

Results

The mean incubation period was estimated to be 12.7 days (standard deviation 4.31 days), indicating that about 4.1% of patients may have incubation periods longer than 21 days.

Conclusion

If the risk of new cases is to be reduced to 1% then 25 days should be used when investigating the source of an outbreak, when determining the duration of surveillance for contacts, and when declaring the end of an outbreak.

Keywordsdisease outbreaks/prevention & controlEbola/epidemiologyEbola hemorrhagic feverEbola/prevention & controlstatistical models
Introduction

Ebola hemorrhagic fever, caused by the Ebola virus (EBOV) is a severe and often fatal disease. The first known outbreak occurred in Zaire (now the Democratic Republic of Congo) in 1976. Since then, single cases and large epidemics have recurred in equatorial Africa. Apart from Zaire, Sudan, Gabon, Ivory Coast, Uganda, Kenya, Angola and the Republic of the Congo have all been affected [1]. Until now, more than 1850 cases and about 1300 deaths have been reported [2]; the latest are from an outbreak in Uganda. Marburg virus, a close relative of EBOV [3], also causes severe outbreaks of hemorrhagic fever in Africa [4]. There are four known subtypes of EBOV named after the place of their first appearance: Zaire, Sudan, Ivory Coast, and Reston in Virginia, USA. The former three cause severe forms of the disease in humans. Subtype Reston is believed not to be dangerous in humans, but it can be fatal in non-human primates [5–7] and pigs [8]. A new strain of EBOV with the proposed name Bundibugyo EBOV was discovered in a Ugandan epidemic in November 2007 [9]. The subtypes have different mortality rates (about 90% for Zaire and 50% for Sudan) and may also have different incubation periods [10,11].

The natural reservoir of the virus remains unknown [12,13]; bats and fruit bats are suspected to be hosts [14,15], and monkeys could be a vector in transmission to humans [16–18]. The two most important routes for person-to-person transmission seem to be direct physical contact and contact with body fluids. It is very likely that droplets or aerosols play little or no role [6,19,20]. When no symptoms are visible, EBOV-infected persons do not seem to be contagious [21,22]. The incubation period is generally assumed to range from 2 to 21 days [10,23–25]; for EBOV most calculations put it at 6 to 10 days [23,26,27]. Knowledge of the incubation period will aid research on EBOV transmission and is fundamental in determining for how long people should be monitored after having had contact with cases. This knowledge is also essential when investigating the animal source that might have triggered an epidemic. In this paper, we investigate the incubation period of EBOV using data collected during an outbreak in Kikwit, Democratic Republic of the Congo, in 1995 [22]. This data set is especially useful because it contains precise information on the contact behavior of 173 persons. We apply an estimation method that was used to study smallpox [28], but has never been used for Ebola hemorrhagic fever.

Materials and MethodsDescription of the epidemic

In 1995, when the outbreak took place, Kikwit had about 200,000 inhabitants and only two hospitals that lacked electricity and running water [29]; protective equipment like gloves was in short supply [30]. On 6 January 1995, a forest worker fell ill with hemorrhagic symptoms and died some days later. In mid-April, a nosocomial outbreak began, leading to the spread of the disease within and between families in Kikwit and nearby areas. On 3 May, containment measures were implemented. One week later, EBOV was confirmed by the Centers for Disease Control and Prevention (Atlanta, USA) [31,32]. Having killed 81% of 315 infected people, the epidemic stopped on 16th July 1995 [5,19,33].

Data set

The data that we used were collected between 17 May and 3 June 1995 [22]. All members of every household from which a primary case had died or was discharged from hospital between 1 January and 7 May 1995 were interviewed. The resultant data set contains information on 27 households with a total of 23 primary cases and 173 contacts persons. It also contains detailed exposure information for each household member. We used this data to classify the contacts into strong contacts (contact with body fluids), weak contacts (direct physical contact but no contact with body fluids) and non-relevant contacts (only non-physical contact).

Model description and methods

Because the exact time of infection is not known, the duration of the incubation period could only be derived from the duration and intensity of contact with the index case. We were also forced to estimate the rates at which contacts were infected during the different periods of the disease in the index cases. Because the available information was not enough to allow us to estimate separate rates for the case’s stay in the hospital and for the funeral, we used a common rate of infection for both periods. To distinguish between the effects of strong and weak contact and between contact at home and contact in hospital and during the funeral, and to keep the number of parameters to a minimum, we proposed three different contact models (Table 1). In Model 1, the infection rates caused by weak contact are the same at home and in hospital, whereas the infection rates produced by strong contact differ in the two places. In Model 2, the infection rates for weak contact differ at home and in hospital, and the infection rates for strong contact are obtained by adding a constant to the infection rate of the corresponding weak contact. In Model 3, the infection rates produced by strong contact differ in the two settings, and the infection rates produced by weak contact are obtained by multiplying the infection rate of the corresponding strong contact with a constant.

The incubation period is assumed to be lognormally distributed with mean μ and standard deviation (SD) σ. This type of distribution is frequently used to describe incubation periods of acute infectious diseases [34,35]. Having parameterized the distribution with mean μ and SD σ, the density function of the lognormal density is given byf(x)=12πσxexp((lnxμ)22σ2),where x is the duration of the incubation period.

We used maximum likelihood inference to estimate the parameters. This method is similar to the procedure proposed by Eichner and Dietz for smallpox incubation periods [25]. The likelihood contribution of each household member who escaped infection, iset0λ(τ)dτand the likelihood contribution of each household member who developed the disease, ist0t1λ(t)et0tλ(τ)dτf(t1t)t,where t0 is the time when the index case first developed symptoms and t1 is the time when the household contact became ill. λ(t) is the infection rate at time t (which depends on the type of contact and the stage of the disease in the index case; see Table 1). f(t1 − t) is the density of the incubation period (see f(x) above) for the delay t1 − t between the household member becoming infected (unknown time t) and the onset of the disease (time t1). Because of lack of information, we assumed that all household members who contracted the infection were infected by the index case in that household.

Results

Of the three models examined, the best was Model 1 with c = 0 (Table 1). The results of the maximum likelihood estimates and the supported ranges are given in Table 2.

The mean incubation period μ was estimated to be 12.7 days with a SD σ of 4.31 days. The force of infection a for all types of contact during hospitalization and the funeral, and for weak contact at home, was 0.0254 per day. This indicates that the probability of escaping infection for one week was 0.84 in spite of having weak contact with the index case at home or at the hospital. The force of infection for strong contacts at home is a + b = 0.161 per day, indicating that the probability of escaping infection for 1 week is, in this case, only 0.32.

Discussion

In this study, the mean incubation period for EBOV hemorrhagic fever was estimated to be 12.7 days (SD 4.31 days). Earlier estimates for EBOV mostly indicated shorter incubation periods ranging from 6 days to 10 days [23,26,27]. Bwaka et al reported a longer incubation period for human-to-human transmission compared to the incubation period for infection caused by a needle prick [26]. In our study, infections were transmitted directly and not by needle prick [36], which may partly explain the longer incubation period estimate. A similar longer mean incubation period of 11.7 days was detected in EBOV infected gorillas [37]. For the Ebola-Sudan strain epidemic in Uganda in 2001, the biggest Ebola epidemic ever reported, the mean incubation period was calculated to range from 6.3 days to 12 days [38]. Bwaka et al and Breman et al estimated a mean incubation period for EBOV of 6.2 and 6.3 days, respectively, but both estimates were based on small data sets [23,26]. Lekone et al calculated a mean incubation period of 10.1 days for the same Kikwit 1995 epidemic [27] for which our data were collected; however, their estimate is not based on the same data set as ours and the authors assumed an exponential distribution for the incubation period. If our parameter estimations were also based on an exponentially distributed incubation period instead of a lognormal one, the estimated mean incubation period would drop to 12.0 days. A more empirical computation that distinguished between a minimum and a maximum incubation period was carried out by Dowell et al [22]. The minimum incubation period was calculated from the death of the index case to the onset of fever in the secondary case. The mean came out to be 7 days, and the range was from 1 day to 15 days. The maximum incubation period was calculated from the onset of fever in the index case to the onset of fever in the secondary case and the mean came out to be 17 days with the range going from 9 days to 25 days. Because infections occur between the onset of fever and the day of death, the mean incubation period should lie between 7 and 17 days.

Our study had a number of limitations: (i) when estimating the parameters, 31 of the 200 household members had to be excluded from the data set because important contact information was missing (two households with four and 19 members, respectively, that included seven secondary cases) or because the possibility that tertiary cases had occurred could not be ruled out (one household with eight members that included three secondary cases); (ii) because we had to estimate the time of infection, there may be more variability in our estimate than if the time of infection had been known precisely; and (iii) different strains of the EBOV behave differently and the number of cases exposed via needles and via direct contact varies in different epidemics, making it difficult to extrapolate from historical data to future outbreaks.

However, our results suggest that a longer incubation period than previously assumed should be taken into consideration when exploring Ebola transmission possibilities and when searching for the original source of the infection that triggered an epidemic. In the past, EBOV transmission occurred when hunters [39] or scientists [40] came into contact with monkeys that were retrospectively assumed to be infected with EBOV. EBOV was propagated in human populations when needles and syringes were reused [19] or when health care workers or family members cared for Ebola cases [22,41]. Important tools in Ebola hemorrhagic fever containment are the surveillance of persons who have been in contact with Ebola-infected patients and the monitoring of convalescent patients before and after they have been discharged from hospital. The length of time these two measures should be carried out and the announcement of the end of an epidemic depend on the length of the incubation period. For all three scenarios that we modeled, a maximum incubation period of 21 days is normally assumed [19,33,36,38,42–44], meaning that in practice surveillance usually stops three weeks after the last contact. Epidemiologists have suggested that two consecutive incubation periods (a total of 42 days) must elapse before an outbreak is declared to be controlled [45,46]. By combining a rigorously collected data set with a new approach for calculating infectious disease incubation periods, our results suggest that the standard assumption of a maximum incubation period of 21 days may be incorrect. Our density function supports durations longer than 21 days. The risk of developing the disease after the 21st day is 4.1 %. To reduce the risk to 1%, 25 days should be taken as the incubation period. A surveillance duration of this length could be an acceptable compromise because: (i) the direct costs (surveillance) and indirect costs (economic cost of quarantines and lost livelihoods) for each day of observation are relatively small; (ii) there is minimal risk of further cases after the 25th day, and (iii) vigilance for Ebola infections by the population and the public health system should be encouraged so that new cases are detected early enough.

ReferencesFeldmannH.GeisbertT.W.Ebola haemorrhagic feverLancet37797682011 Mar 584986215 Nov 2010 [Epub ahead of print]21084112WHOEbola haemorrhagic fever — fact sheet revised in May 2004 [monograph on the Internet]Available at:2004World Health OrganizationGenevahttp://www.who.int/wer/2004/wer7949.pdf[Date accessed: 17 December 2008]ChlibekR.SmetanaJ.VackovaM.Ebola and Marburg fever — outbreaks of viral haemorrhagic feverKlin Mikrobiol Infekc Lek1262006 Dec21722317230375Marburg hemorrhagic fever — Uganda (03). August 10, 2007. Available at: http://www.promedmail.org archive 20070810.2609 [Date accessed: 08 December 2008].CDCEbola Hemorrhagic Fever information packet [monograph on the Internet]Available at:2002CDC, Special Pathogens Branch — Division of Viral and Rickettsial DiseasesAtlantahttp://www.cdc.gov/ncidod/dvrd/spb/mnpages/dispages/Fact_Sheets/Ebola_Fact_Booklet.pdf[cited 2008 Dec 8, Date accessed: 08 December 2008]KlenkH.D.FeldmannH.Ebola and Marburg viruses: molecular and cellular biology2004Horizon BioscienceNorfolk, UKMcCormickJ.B.Ebola virus ecologyJ Infect Dis190112004 Dec 11893189415529250NormileD.Emerging infectious diseases. Scientists puzzle over Ebola-Reston virus in pigsScience23359132009 Jan 2345119164717TownerJ.S.SealyT.K.KhristovaM.L.Newly discovered ebola virus associated with hemorrhagic fever outbreak in UgandaPLoS Pathog4112008 Nove100021219023410MacNeilA.FarnonE.C.WamalaJ.Proportion of deaths and clinical features in Bundibugyo Ebola virus infection, UgandaEmerg Inf Dis16122010 Dec19691972ChowellG.HengartnerN.W.Castillo-ChavezC.The basic reproductive number of Ebola and the effects of public health measures: the cases of Congo and UgandaJ Theor Biol22912004 Jul 711912615178190FeldmannH.Wahl-JensenV.JonesS.M.Ebola virus ecology: a continuing mysteryTrends Microbiol129102004 Oct43343715381189LeirsH.MillsJ.N.KrebsJ.W.Search for the Ebola virus reservoir in Kikwit, Democratic Republic of the Congo: reflections on a vertebrate collectionJ Infect Dis179Suppl. 11999 FebS155S1639988179HaymanD.T.S.EmmerichP.YuM.Long-term survival of an urban fruit bat seropositive for Ebola and Lagos bat virusesPLoS One582010 Aug 4e1197820694141LeroyE.M.EpelboinA.MondongeV.Human Ebola outbreak resulting from direct exposure to fruit bats in Luebo, Democratic Republic of CongoVector-borne and Zoonotic Diseases962009 Dec72372819323614LeroyE.M.KumulunguiB.PourrutX.Fruit bats as reservoirs of Ebola virusNature43870682005 Dec 157557616319873LeroyE.M.RouquetP.FormentyP.Multiple Ebola virus transmission events and rapid decline of central African wildlifeScience30356562004 Jan 1538739014726594BiekR.WalshP.D.LeroyE.M.Recent common ancestry of Ebola Zaire virus found in a bat reservoirPLoS Pathog2102006 Octe9017069458KhanA.S.TshiokoF.K.HeymannD.L.The reemergence of Ebola hemorrhagic fever, Democratic Republic of the Congo, 1995. Commission de Lutte contre les Epidemies a KikwitJ Infect Dis179Suppl. 11999 FebS76S869988168BaronR.C.McCormickJ.B.ZubeirO.A.Ebola virus disease in southern Sudan: hospital dissemination and intrafamilial spreadBull World Health Organ616198399710036370486BennettD.BrownD.Ebola virusBMJ31069911995 May134413457787519DowellS.F.MukunuR.KsiazekT.G.Transmission of Ebola hemorrhagic fever: a study of risk factors in family members, Kikwit, Democratic Republic of the Congo, 1995. Commission de Lutte contre les Epidemies a KikwitJ Infect Dis179Suppl. 11999 FebS87S919988169BremanJ.G.PiotP.JohnsonK.M.The epidemiology of Ebola hemorrhagic fever in Zaire, 1976PattynS.Ebola Virus Haemorrhagic Fever. Proceedings of the international colloquium on Ebola virus infections and other haemorrhagic fevers, 1977 Dec 6–8, Antwerpen, Belgien1978Elsevier/North Holland Biomedical PressAmsterdamGeisbertT.W.JahrlingP.B.Exotic emerging viral diseases: progress and challengesNat Med10Suppl. 122004 DecS110S12115577929NyamathiA.M.FaheyJ.L.SandsH.Ebola virus: immune mechanisms of protection and vaccine developmentBiol Res Nurs442003 Apr27628112698920BwakaM.A.BonnetM.J.CalainP.Ebola hemorrhagic fever in Kikwit, Democratic Republic of the Congo: clinical observations in 103 patientsJ Infect Dis179Suppl. 11999 FebS1S79988155LekoneP.E.FinkenstädtB.F.Statistical inference in a stochastic epidemic SEIR model with control intervention: Ebola as a case studyBiometrics6242006 Dec1170117717156292EichnerM.DietzK.Transmission potential of smallpox: estimates based on detailed data from an outbreakAm J Epidemiol15822003 Jul 1511011712851223HewlettB.L.HewlettB.S.Providing care and facing death: nursing during Ebola outbreaks in central AfricaJ Transcult Nurs16402005 Oct28929716160191HeymannD.L.BarakamfitiyeD.SzczeniowskiM.Ebola hemorrhagic fever: lessons from Kikwit, Democratic Republic of the CongoJ Infect Dis179Suppl. 11999 FebS283S2869988197GuimardY.BwakaM.A.ColebundersR.Organization of patient care during the Ebola hemorrhagic fever epidemic in Kikwit, Democratic Republic of the Congo, 1995J Infect Dis179Suppl. 11999 FebS268S2739988194Muyembe-TamfumJ.J.KipasaM.KiyunguC.Ebola outbreak in Kikwit, Democratic Republic of the Congo: discovery and control measuresJ Infect Dis179Suppl. 11999 FebS259S2629988192IVSViral haemorrhagic feversEuro surveill72002315412631941HorowitzM.L.CohenN.D.TakaiS.Application of Sartwell’s model (lognormal distribution of incubation periods) to age at onset and age at death of foals with Rhodococcus equi pneumonia as evidence of perinatal infectionJ Vet Intern Med15932001 May–Jun17117511380023SartwellP.E.The distribution of incubation periods of infectious diseaseAm J Hyg5131950 May31031815413610Centres for Diseases Control and PreventionUpdate: Outbreak of Ebola hemorrhagic fever — ZaireMMWR Morb Mortal Wkly Rep44191995 May 193813827739512BermejoM.Rodriguez-TeijeiroJ.D.IlleraG.Ebola outbreak killed 5000 gorillasScience31458052006 Dec 8156417158318OkwareS.I.OmaswaF.G.ZarambaS.An outbreak of Ebola in UgandaTrop Med Int Health7122002 Dec1068107512460399WolfeN.D.ProsserT.A.CarrJ.K.Exposure to nonhuman primates in rural CameroonEmerg Infect Dis2004 DecAvailable at:http://www.cdc.gov/ncidod/EID/vol10no12/04–0062.htm[serial online, cited 2008 Dec 19, Date accessed: 08 December 2008]WHOWHO recommended Guidelines for Epidemic Preparedness and Response: Ebola Haemorrhagic Fever (EHF) [monograph on the Internet]Available at:1997World Health OrganizationGenevahttp://www.who.int/csr/resources/publications/ebola/whoemcdis977E.pdf[cited 2008 Nov 25, Date accessed: 08 December 2008]SadekR.F.KhanA.S.StevensG.Ebola hemorrhagic fever, Democratic Republic of the Congo, 1995: determinants of survivalJ Infect Dis179Suppl. 11999 FebS24S279988161FrancesconiP.YotiZ.DeclichS.Ebola hemorrhagic fever transmission and risk factors of contacts, UgandaEmerg Infect Dis2003 NovAvailable at:http://www.cdc.gov/ncidod/EID/vol9no11/03–0339.htm[serial online, cited 2008 Nov 25, Date accessed: 08 December 2008]LloydE.S.ZakiS.R.RollinP.E.Long-term disease surveillance in Bandundu region, Democratic Republic of the Congo: a model for early detection and prevention of Ebola hemorrhagic feverJ Infect Dis179Suppl. 11999 FebS274S2809988195WHOOutbreak(s) of Ebola haemorrhagic fever in the Republic of the Congo, January–April 2003Wkly Epidemiol Rec78200328529614509121HoppM.Ebola Hemorrhagic Fever — Congo DR (11): WHO. ProMedAvailable at:October 3, 2007http://www.promedmail.orgarchive 20071003.3270 [Date accessed: 08 December 2008]MasonC.The strains of EbolaCMAJ178102008 May 61266126718458253Acknowledgements

This work was supported by the EU project INFTRANS (FP6 STREP; contract no. 513715), by the MODELREL project funded by DG SANCO (no. 2003206-SI 2378802) and by the German Ministry of Health. We thank to Mrs. Mariana Nold for her expert help with data handling and parameter inference.

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Infection rates λ(t) for different contact types and disease periods for the Ebola outbreak in Kikwit 1995. a, b and c are parameters which are used to construct the infection rate for each day on which a contact between a case and a household member is reported

At home
During hospitalization and funeral
Weak contactaStrong contactbWeak contactaStrong contactb
Model 1aa + baa + c
Model 2aa + bcc + b
Model 3a × bba × cc

Weak contact refers to direct physical contact.

Strong contact additionally involves contact with body fluids. For a detailed explanation, see the text.

Maximum likelihood estimates and supported ranges for mean μ and standard deviation σ of the incubation period and for transmission possibilities of Ebola hemorrhagic fever in Kikwit in 1995

ParameterMLELower supported rangeUpper supported range
μ12.7 d10.1 d16.1 d
σ4.31 d2.56 d8.60 d
a0.0254/d0.0118/d0.0466/d
b0.136/d0.00550/d0.370/d