|
|
||||||||
| ABSTRACT |
|
|
|---|
| INTRODUCTION |
|
|
|---|
Ebola hemorrhagic fever was not reported again until the end of 1994, when three outbreaks started almost simultaneously. In November 1994, ethnologists studying chimpanzees (Pan troglodytes verus) in Tai National Park in Côte dIvoire found dead chimpanzees and noticed the absence of others. One dead chimpanzee was found to be infected with Ebola virus, and one researcher was infected with Ebola during a chimpanzee necropsy.6,7 The following month, Ebola outbreaks were reported in Gabon and the Democratic Republic of Congo. Multiple human cases were reported in northeastern Gabon in the gold panning camps of Mekouka, Andock, and Minkebe,8,9 and a large human Ebola outbreak began in the Kikwit District in the Democratic Republic of Congo.10 The Kikwit outbreak resembled the 1976 episodes in Yambuku, Democratic Republic of Congo where secondary transmission of the virus in Kikwit occurred through close personal contact in families and hospitals where infection control mechanisms were not in place because of economic constraints. Retrospective case analysis suggests that the index case may have been a charcoal maker that worked in the forest outside Kikwit and was presumably exposed through an unknown mechanism to the virus. In Kikwit, human-to-human transmission occurred without being recognized until the end of April 1995. Additional human outbreaks were reported in February 1996 in Mayibout II, Gabon, a village 40 km south of the original outbreak in the gold panning camps, and in July 1996 at a logging camp between Ovan and Koumameyong, near Booue.9
Ape carcasses were found near the sites of three Ebola hemorrhagic fever outbreaks in Gabon; however from 1994 to 1997, only one chimpanzee was found to be Ebola positive in Lopé Park near Booué in August 1996. Circumstantial reports attribute the February 1996 Gabon outbreak to infection after butchering a dead chimpanzee, although this has not been verified.9
The largest Ebola hemorrhagic fever epidemic to date occurred in the Gulu District in Uganda from August 2000 to January 2001.11,12 Epidemiologic investigation of this epidemic showed intrafamilial and nosocomial transmission similar to the Kikwit outbreak.
From October 2001 to July 2002, several Ebola outbreaks were reported in the Ogooue-Ivindo Province of Gabon and in the Mbomo and Kéllé districts in the Republic of the Congo.13 In August and September 2001, a die-off of chimpanzees, gorillas, and duikers was reported in Gabon.14 Then in December 2002, another outbreak of Ebola occurred in humans in Kéllé. This outbreak was again preceded by a die-off of chimpanzees, gorillas, and duikers in November 2002 that were laboratory-positive for Ebola infection.15
Ebola virus affects non-human primates, and recent work suggested a gradual spread of Ebola hemorrhagic fever virus within gorilla and chimpanzee communities in western equatorial Africa with a dramatic decrease noted in population numbers.16 Non-human primates have been implicated as the source of several, but not all, human outbreaks through contact with infected ape meat acquired through hunting. We consider the possibility that Ebola outbreaks in non-human primates may occur independently of human cases. This raises the potential problem of distinguishing primary Ebola outbreaks from undetected, residual outbreaks that reemerge from a non-human primate source, usually through the bush meat trade.
Despite extensive field investigations to define the ecology of the virus, the mechanism of transmission from its reservoir to non-human primates and humans remains a mystery. Filoviruses do not persist in experimentally infected non-human primates. Non-human primates are probably not the natural reservoir and, similar to humans, these species are probably infected when direct or indirect contact is made with the natural host. No Ebola viruses or antibodies against Ebola virus were identified in the more than 32,000 animals including arthropods, mammals, avians, and reptiles tested at several sites in Africa as potential candidate reservoir species.1215,1721
Multiple factors likely contribute to the appearance of Ebola hemorrhagic fever in non-human primates and humans. Of interest is the spatial specificity, seasonal context, and occasional temporal clustering of Ebola hemorrhagic fever outbreaks. Ebola epidemics in Nzara, Sudan and Yambuku, Democratic Republic of Congo in 1976 occurred within two months of each other in two geographic locations separated by hundreds of kilometers involving two separate viral strains (Sudan and Zaire). The outbreaks of Tai, Côte dIvoire, Mekouka, Gabon, and Kikwit, Democratic Republic of Congo in late 1994 also occurred within two months of each other in three different geographic regions involving two viral strains (Côte dIvoire and Zaire). Fifteen years passed between the 19761979 and 19941996 temporal clusters of Ebola cases without identification of additional outbreaks in equatorial Africa. The late 2000 Uganda outbreak and the late 2001 Gabon outbreaks occurred seven and eight years after the 1994 outbreaks.
Prior studies have suggested possible enviroclimatic influence of Ebola fever incidence due to the predominant appearance of epidemics at the end of the rainy season and/or the start of the dry season.6,22,23 Tucker and others analyzed Landsat data from all reported Ebola outbreak locations and found either gallery tropical forest or continuous tropical forest at every outbreak location, in addition to a significant negative climatological anomaly associated with the 1994 cluster of outbreaks.24 These initial studies suggest Ebola virus emerges from its cryptic reservoir in a specific geotemporal and enviroclimatic context. While we acknowledge the multi-factorial nature of epidemic triggering and propagation, these studies implied an ongoing need to evaluate remotely sensed data for use in predictive disease risk modeling, otherwise known as Remotely Sensed Epidemic Surveillance.25 Here we investigate the possibility that enviroclimatic trigger events preceded all initial Ebola outbreaks from 1994 to 2002, using time series satellite data from July 1981 to the present.
| MATERIALS AND METHODS |
|
|
|---|
|
High values of NDVI (~+0.6) are representative of dense green canopies with a high concentration of chlorophyll, while low values (~0.1) are indicative of sparse vegetation cover and stressed vegetation.31,32 The NDVI data provide a temporally and spatially consistent inventory of global gridded time series at 8 km resolution using bimonthly maximum composites.2833
A number of studies have shown that NDVI time series are a surrogate for vegetation response to rainfall and evapotranspiration in a wide range of environmental conditions. The NDVI is a surrogate for photosynthetic capacity since it is highly correlated to the absorbed fraction of photo-synthetically active radiation (FPAR) and thus gross photosynthesis.33 We use the NDVI to infer the FPAR, which is directly influenced by rainfall. The NDVI time series has been used to investigate ecologic dynamics and enviroclimatically coupled diseases.3438
The documented Ebola hemorrhagic fever outbreaks that occurred within the 19812003 satellite record were located by the closest geographic coordinates and their respective NDVI were extracted. The mean value, variance, and seasonal profiles of NDVI were used to identify regions (8 x 8 km grid points) of continental Africa with similar characteristics to those from the Ebola sites.
Figure 1
shows the NDVI seasonal patterns for the sites of reported Ebola outbreaks within the satellite record and the NDVI time series during that period. These signals present a characteristic bimodal/unimodal rainy season pattern proper of tropical forest and tropical gallery forest within a savanna matrix of equatorial Africa. Of interest is the unique unimodal seasonal context of the Democratic Republic of Congo and Uganda Ebola incidences compared with the bimodal rainy season patterns from the Gabon, Congo and Côte dIvoire signals. The Nzara, Sudan site was added for comparison with the Democratic Republic of Congo and Uganda sites because of noted similarities in the unimodal seasonal behavior; however, cases of Ebola were not documented in Sudan within the time period of the satellite record analyzed. The bimodal signals show a large deficit of the NDVI with respect to the seasonal pattern (negative anomaly) during the months prior to first appearance of the virus. The Kéllé, Republic of Congo NDVI time series does not show this negative anomaly; however, the Mekambo, Gabon time series does, which supports the observation that the cases in Kéllé, Republic of Congo represented an extension of the Mekambo, Gabon outbreak. Note also that a negative NDVI anomaly was found at the Kikwit site during the first months of 1994, but was not associated with the immediate outbreak. A major drought was observed in 1994 in the immediate areas of the bimodal signals (i.e., the Côte dIvoire and Gabon sites) but not the unimodal sites (i.e., the Democratic Republic of Congo and Uganda sites).
|
0.92) using a two-tailed t-test. High spatial coherence was found among the pixels at this confidence level, i.e., the distance to Ebola outbreaks were less than 220 km or 2 degrees.
|
The SVD is a method of finding linear combinations of unitary matrices U and V such that the individual covariances are maximized. The solution given by an SVD matrix decomposition is [U, S, V] = SVD(A), where A = USV.
The notation means the transpose of a vector or matrix. In this expression, U and V are matrices whose columns are orthonormal (unitary matrices) and S is a diagonal matrix whose diagonal elements are the positive singular values sort in descending order.
The CCA is a multivariate statistical technique used to examine and describe the strength of a linear association between two sets of random variables (left and right fields).3941 This analysis produces a set of orthonormal vectors for the satellite NDVI signals (U) and Ebola sites (V), which are referred to as the canonical factors and are sorted in descending order in terms of the percentage of the covariance explained by the canonical factors (S). Of great importance is that the leading paired-modes between annual NDVI images and annual NDVI time series of regions of Ebola incidences from the CCA are the largest signal involved between them in both the temporal and spatial domains. In this case, this largest signal can be used to identify for all Africa, regions that can be associated with Ebola incidences through the canonical factor leading modes, and thus develop a temporal-spatial risk index map for Africa.
To focus the analysis on the dominant modes of covariability, we performed CCA using standardized, seasonal-mean leading principal factors.39 Before the results are shown, a few preliminary comments regarding statistical significance are in order. As in many applications of CCA, we are dealing with fields in which the number of grid points used have a relatively small number of realizations in the time domain. The number of degree of freedom in the correlation coefficients derived from the normalized NDVI data set would be 22 in the absence of year-to-year correlation. In this report, we rely on simple two-tailed t-tests performed on 64 Monte Carlo simulations for assessing the statistical significance of the results in which 5% of the grid points were selected randomly. A 5% significance level (r
0.7) using a two-tailed t-test is reported.
Since the CCA isolates linear combinations of variables that optimize the correlations within two fields, a criterion to determine when a pixel is anomalous can be defined in terms of the integration of the first five canonical vectors for each pixel with a negative anomaly, hereafter called the Ebola Geo-Temporal Trigger (EGT2) index
![]() |
where U, S, and V are the canonical components given by the CCA, Nc represents the five most important components used in this study, and the notation represents the transpose of a vector or matrix. A pixel with an EGT2 index > 0.9 is counted as a pixel involved in each trigger event. For display purposes, this index is normalized by their maximum value (2.5). Adding all EGT2 pixel values for a particular year provides an annual overall Trigger Index (R).
| RESULTS |
|
|
|---|
|
We found several important criteria for definition of an Ebola trigger event. Table 2
shows the number of pixels (highlighted) associated with a high EGT2 index for all of Africa for each year (Np); the average of minimum distance to Ebola sites (D); the overall transmission trigger index (R); the percentage of Np associated with each individual Ebola site (% of Np); and the minimum/mode distances from a trigger event to an Ebola site (d). Notice the high concentration of anomalous grid points around Gabon sites (50% with minimum distance of 2 km and mode distance less than 100 km). The repeated occurrences of environmental stress conditions in this region are a warning of possible future outbreaks of this enigmatic disease. Figure 4
shows the transmission risk R values in Table 2
obtained from the results of the CCA in the spatial domain for each year. For display purposes the R values less than 2.75 were rescaled by one-third. Note that Figure 4
shows 1994 as a year of a significant trigger event that spanned Côte dIvoire, Gabon, and the Democratic Republic of Congo. Years highlighted with no associated appearance of Ebola, such as 1989, were associated with a low overall transmission risk index (< 3.5), indicating low environmental transmission competency potential. The year 1991 was the year with maximum number of pixels and trigger index highlighted. This was due to the general disruption in recorded NDVI over global tropical forest after the eruption of Mount Pinatubo.42,43
|
|
|
| DISCUSSION |
|
|
|---|
We recognized the possibility of introducing errors into our analysis due to 1) the hypothesis that Ebola can be transmitted among non-human primate species, 2) the difficulty in determining a primary Ebola outbreak from a secondary or tertiary one; and 3) undetected cases likely occurring. Lastly, the emergence of Ebola virus might not involve primates but other animals such as bats or rats. Nevertheless, we found two time periods during the past 23 years when trigger events were present in equatorial Africa and no Ebola outbreaks were reported (i.e., 1989, 1991; Figures 4
and 5
). We interpret this to mean that either no transmission of the Ebola virus occurred in non-human primates and humans or an unidentified epidemic may have occurred at that time.
It has been previously suggested the consistent appearance of Ebola during the rainy season may be related to increased numbers of insects and mammals.6 Insect biomass in Makokou, Gabon has been reported to be 3.9 times greater during the rainy season from mid-October to mid-December.44 Dramatic seasonal changes in fruiting and foraging behavior in primates, squirrels, and birds have been documented in Côte dIvoire and Gabon, where diet overlap and competition for the same food source has been observed during the period of maximum fruit production during the rainy season.6,4447 The trigger event described in this study, beyond the seasonal enviroclimatic influences on Ebola transmission, may provide an added influence that results in greater contact of candidate host insect and animals that become exposed to the cryptic Ebola reservoir through changes in foraging behavior.
The findings of this investigation suggest that highlighting conditions favorable for Ebola virus transmission is critical information for health care workers in Côte dIvoire, the Republic of Congo, the Democratic Republic of Congo, Uganda, and Gabon. We will continue our analysis and publicize our predictions to enable confirmation or refutation of our predictions. We have mapped the areas of potential Ebola outbreaks based on the time series analysis for the unimodal and bimodal rainy season locations shown in Figure 1
. If our analysis is correct, future Ebola outbreaks will only occur in these areas and not elsewhere in equatorial Africa (Figure 5
).
Received February 9, 2004. Accepted for publication April 24, 2004.
Acknowledgments: We thank Brian Rothman (University of Cincinnati College of Medicine), Neal Woolen (United States Army Medical Research Institute of Infectious Diseases), David Heymann, Guenael Rodier (World Health Organization Communicable Diseases Cluster), Richard Hatchett (United States Department of Health and Human Services), and Dr. Robert Vanessea for their valuable assistance.
Financial support: The appointments of James M. Wilson to the World Health Organization Ebola Tai Forest Project and the NASA-Goddard Space Flight Center were supported jointly by the Department of Communicable Diseases Surveillance and Response, the World Health Organization, and the Office of Applications, NASA.
Authors addresses: Jorge E. Pinzon, Biospheric Sciences Branch, Laboratory for Terrestrial Physics, Code 923, National Aeronautics and Space Administration-Goddard Space Flight Center, Greenbelt, MD 20771 and Science Systems & Applications, Inc., Lanham, MD 20706, Telephone: 301-614-6685, Fax: 301-614-6015, E-mail: pinzon{at}negev.gsfc.nasa.gov. James M. Wilson, Global Alert and Response Team, Department of Communicable Diseases Surveillance and Response, World Health Organization, Geneva, Switzerland, E-mail: wilson{at}isis.imac.georgetown.edu. Compton J. Tucker, Sciences Branch, Laboratory for Terrestrial Physics, Code 923, National Aeronautics and Space Administration-Goddard Space Flight Center, Greenbelt, MD 20771, E-mail: compton{at}kratmos.gsfc.nasa.gov. Ray Arthur and Pierre Formenty, Global Alert and Response Team, Department of Communicable Diseases Surveillance and Response, World Health Organization, Geneva, Switzerland, E-mails: rca8{at}cdc.gov and formenty{at}who.int. Peter B. Jahlring, U.S. Army Medical Research Institute for Infectious Disease, Fort Detrick, Frederick, MD 21702-5011, E-mail: Peter.Jahring{at}det.amedd.army.mil.
Reprint requests: Jorge E. Pinzon, Laboratory for Terrestrial Physics, Code 923, NASA/Goddard Space Flight Center, Greenbelt, MD 20771, Telephone: 301-614-6685, Fax: 301-614-6015, E-mail: pinzon{at}negev.gsfc.nasa.gov.
| REFERENCES |
|
|
|---|
This article has been cited by other articles:
![]() |
H. Frumkin, J. Hess, G. Luber, J. Malilay, and M. McGeehin Climate Change: The Public Health Response Am J Public Health, March 1, 2008; 98(3): 435 - 445. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. J. Wittmann, R. Biek, A. Hassanin, P. Rouquet, P. Reed, P. Yaba, X. Pourrut, L. A. Real, J.-P. Gonzalez, and E. M. Leroy Isolates of Zaire ebolavirus from wild apes reveal genetic lineage and recombinants PNAS, October 23, 2007; 104(43): 17123 - 17127. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. T. PETERSON, R. R. LASH, D. S. CARROLL, and K. M. JOHNSON GEOGRAPHIC POTENTIAL FOR OUTBREAKS OF MARBURG HEMORRHAGIC FEVER. Am J Trop Med Hyg, July 1, 2006; 75(1): 9 - 15. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |