This is only an attractor.
Please read the original page athttp://www.plospathogens.org/article/info%3Adoi%2F10.1371%2Fjournal.ppat.1000212
Newly Discovered Ebola Virus Associated with Hemorrhagic Fever Outbreak in Uganda
Jonathan S. Towner1, Tara K. Sealy1, Marina L. Khristova2, César G. Albariño1, Sean Conlan3, Serena A. Reeder1, Phenix-Lan Quan3, W. Ian Lipkin3, Robert Downing4, Jordan W. Tappero4, Samuel Okware5, Julius Lutwama6, Barnabas Bakamutumaho6, John Kayiwa6, James A. Comer1, Pierre E. Rollin1, Thomas G. Ksiazek1, Stuart T. Nichol1*
1 Special Pathogens Branch, Centers for Disease Control and Prevention, Atlanta, Georgia, United States of America, 2 Scientific Resources Program, Centers for Disease Control and Prevention, Atlanta, Georgia, United States of America, 3 Center for Infection and Immunity, Mailman School of Public Health, Columbia University, New York, New York, United States of America, 4 Global AIDS Program, Centers for Disease Control and Prevention, Entebbe, Uganda, 5 Ministry of Health, Republic of Uganda, Kampala, Uganda, 6 Uganda Virus Research Institute, Entebbe, Uganda
Over the past 30 years, Zaire and Sudan ebolaviruses have been responsible for large hemorrhagic fever (HF) outbreaks with case fatalities ranging from 53% to 90%, while a third species, Côte d'Ivoire ebolavirus, caused a single non-fatal HF case. In November 2007, HF cases were reported in Bundibugyo District, Western Uganda. Laboratory investigation of the initial 29 suspect-case blood specimens by classic methods (antigen capture, IgM and IgG ELISA) and a recently developed random-primed pyrosequencing approach quickly identified this to be an Ebola HF outbreak associated with a newly discovered ebolavirus species (Bundibugyo ebolavirus) distantly related to the Côte d'Ivoire ebolavirus found in western Africa. Due to the sequence divergence of this new virus relative to all previously recognized ebolaviruses, these findings have important implications for design of future diagnostic assays to monitor Ebola HF disease in humans and animals, and ongoing efforts to develop effective antivirals and vaccines.
In this report we describe a newly discovered ebolavirus species which caused a large hemorrhagic fever outbreak in western Uganda. The virus is genetically distinct, differing by more than 30% at the genome level from all other known ebolavirus species. The unique nature of this virus created challenges for traditional filovirus molecular based diagnostic assays and genome sequencing approaches. Instead, we quickly determined over 70% of the virus genome using a recently developed random-primed pyrosequencing approach that allowed the rapid development of a molecular detection assay that was deployed in the disease outbreak response. This draft sequence allowed easy completion of the whole genome sequence using a traditional primer walking approach and prompt confirmation that this virus represented a new ebolavirus species. Current efforts to design effective diagnostics, antivirals and vaccines will need to take into account the distinct nature of this important new member of the filovirus family.
Citation: Towner JS, Sealy TK, Khristova ML, Albariño CG, Conlan S, et al. (2008) Newly Discovered Ebola Virus Associated with Hemorrhagic Fever Outbreak in Uganda. PLoS Pathog 4(11): e1000212. doi:10.1371/journal.ppat.1000212
Editor: Christopher F. Basler, Mount Sinai School of Medicine, United States of America
Received: August 6, 2008; Accepted: October 20, 2008; Published: November 21, 2008
This is an open-access article distributed under the terms of the Creative Commons Public Domain declaration which stipulates that, once placed in the public domain, this work may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose.
Funding: All funding for this work was provided by the Centers for Disease Control and Prevention.
Competing interests: The authors have declared that no competing interests exist.
* E-mail: firstname.lastname@example.org
The family Filoviridae consists of two genera, Marburgvirus and Ebolavirus, which have likely evolved from a common ancestor . The genus Ebolavirus is comprised of four species, Zaire, Sudan, Reston and Côte d'Ivoire (Ivory Coast) ebolaviruses, which have, with the exception of Reston and Côte d'Ivoire ebolaviruses, been associated with large hemorrhagic fever (HF) outbreaks in Africa with high case fatality (53–90%) . Viruses of each species have genomes that are at least 30–40% divergent from one another, a level of diversity that presumably reflects differences in the ecologic niche they occupy and in their evolutionary history. Identification of the natural reservoir of ebolaviruses remains somewhat elusive, although recent PCR and antibody data suggest that three species of arboreal fruit bats may be carriers of Zaire ebolavirus. No data has yet been published to suggest reservoirs for the Sudan, Reston and Côte d'Ivoire ebolavirus species. However, a cave-dwelling fruit bat has been recently implicated as a natural host for marburgvirus ,, supporting the hypothesis that different bat species may be the reservoir hosts for the various filoviruses.
Filovirus outbreaks are sporadic, sometimes interspersed by years or even decades of no apparent disease activity. The last new species of ebolavirus was discovered 14 years ago (1994), in Cote d'Ivoire (Ivory Coast), and involved a single non-fatal case, a veterinarian who performed an autopsy on an infected chimpanzee found in the Tai Forest . No further disease reports have been associated with Côte d'Ivoire ebolavirus, in contrast to Zaire and Sudan ebolaviruses which have each caused multiple large outbreaks over the same time period. Here, we report the isolation and characterization of a new ebolavirus that was responsible for a large hemorrhagic fever outbreak in western Uganda. This new virus, proposed name Bundibugyo ebolavirus, is most closely related, albeit distantly, to Côte d'Ivoire ebolavirus, which was somewhat unexpected, given the large geographic distance between the two countries of origin.
In late November 2007 HF cases were reported in the townships of Bundibugyo and Kikyo in Bundibugyo District, Western Uganda (Figure 1A). A total of 29 blood samples were initially collected from suspect cases and sent in two air-transport shipments to the Centers for Disease Control and Prevention (CDC) for immediate testing. Evidence of acute ebolavirus infection was detected in eight specimens using a broadly reactive ebolavirus antigen capture assay known to cross-react with the different ebolavirus species  and an IgM capture assay based on Zaire ebolavirus reagents (Table 1). The Ugandan Ministry of Health was notified on November 28, 2007. These specimens were negative when initially tested with highly sensitive real-time RT-PCR assays specific for all known Zaire and Sudan ebolaviruses and marburgviruses. However, further evidence of acute ebolavirus infection was obtained using a traditionally less sensitive (relative to the real-time RT-PCR assays) but more broadly reactive filovirus L gene-specific RT-PCR assay (1 specimen) (Table 1). Sequence analysis of the PCR fragment (400 bp of the virus L gene) revealed the reason for the initial failure of the real-time RT-PCR assays was that the sequence was distinct from the four known species of ebolavirus, although it was distantly related to Cote d'Ivoire ebolavirus. In total, 9 of 29 specimens showed evidence of ebolavirus infection, and all tests were negative for marburgvirus (data not shown).
Figure 1. Geographic locations of Ebola HF outbreaks and phylogenetic relationships of representative filoviruses.
(A) Map of Africa showing the sites of all known ebolavirus outbreaks denoted by colored circles for Zaire ebolavirus (yellow), Sudan ebolavirus (green), and Côte d'Ivoire ebolavirus (red). The expanded map of Uganda shows the location of the communities of Bundibugyo and Kikyo (black circles) in western Uganda, the site of the recent outbreak of Bundibugyo ebolavirus. Also shown on the Uganda map are the cities of Kampala (capital), Entebbe (international airport) and Gulu (the site of an outbreak of Sudan ebolavirus in 2000, the largest known Ebola HF outbreak on record). (B) Phylogenetic tree comparing full-length genomes of ebolavirus and marburgvirus by Bayesian analysis. Posterior probabilities greater than 0.5 and maximum likelihood bootstrap values greater than 50 are indicated at the nodes.
Table 1. Ebolavirus diagnostic results of initial 29 specimens obtained from Bundibugyo District with numerical specimen numbers assigned.
Approximately 70% of the virus genome was rapidly (in less than 10 days) sequenced from total RNA extracted from a patient serum (#200706291) using a recently established metagenomics pyro-sequencing method developed at 454 Life Sciences which involves successive rounds of random DNA amplification . Using the newly derived draft sequence, a real-time RT-PCR assay specific for the NP gene of this virus was quickly developed and evaluated. The assay was shown to have excellent sensitivity (Table 1), finding positive all the initial six samples that tested positive by either virus antigen capture (five specimens) or virus isolation assays (four specimens). The antigen-capture, IgM, IgG and newly designed real-time PCR assays were quickly transferred to the Uganda Virus Research Institute during the course of the outbreak to facilitate rapid identification and isolation of Ebola cases in the affected area for efficient control of the outbreak. The outbreak continued through late December, 2007, and resulted in 149 suspected cases and 37 deaths .
The entire genome sequence of this virus was completed using a classic primer walking sequencing approach on RNA from the reference virus isolate (#811250). The complete genome of the Côte d'Ivoire ebolavirus was not available, so it too was derived by a similar combination of random primed pyrosequencing and primer walking approaches. Acquisition of these sequences allowed for the first time the phylogenetic analysis of the complete genomes of representatives of all known species of Ebola and Marburg viruses. The analysis revealed that the newly discovered virus differed from the four existing ebolavirus species (Figure 1B), with approximately 32% nucleotide difference from even the closest relative, Côte d'Ivoire ebolavirus (Table 2). Similar complete genome divergence (35–45%) is seen between the previously characterized ebolavirus species. Bundibugyo ebolavirus is the proposed name for this newly identified species. This high level of genetic diversity translates to considerable amino acid differences in the encoded virus proteins. The virus surface glycoprotein, which is an important determinant of virus tropism and pathogenicity, differs between Bundibugyo and Côte d'Ivoire and Zaire ebolaviruses by over 27 and 35%, respectively at the amino acid level. Similarly, important virulence factors such as the immunosuppressive VP35 protein differ between Bundibugyo and Côte d'Ivoire and Zaire ebolaviruses by 23 and 25%, respectively. This extent of divergence will likely be reflected in significant antigenic and pathogenicity differences among these viruses.
Table 2. Identity (percent) matrix based on comparisons of full-length genome sequences.
Current human prototype ebolavirus vaccines include Zaire and Sudan ebolaviruses–. Cross-protection studies will need to be done to assess whether vaccine designs will need to incorporate the Bundibugyo ebolavirus. The unique nature of this virus has other implications too, including screening of potential antivirals and pathogenicity studies. Retrospective analysis of case description, epidemiologic and laboratory data from the Bundibugyo outbreak are still ongoing, but it is clear that the case fatality (~36%) associated with Bundibugyo ebolavirus infection is lower than that observed for Zaire ebolavirus (approx. 80–90%) and Sudan ebolavirus (approximately 50–55%) . Studies in non-human primates need to be performed to compare the pathogenicity of these viruses. This investigation also highlights the power of molecular detection and characterization tools to quickly identify new pathogens, while providing a cautionary note regarding sole dependence on molecular techniques such as real-time PCR assays for detection of novel agents in biodefense or emerging disease surveillance programs.
Materials and Methods
Ebolavirus detection and virus isolation
Several diagnostic techniques were used for each sample: (i) antigen capture, IgG, and IgM assays were performed as previously described , (ii) virus isolation attempts were performed on Vero E6 cells  and monitored for 14 days; (iii) RNA was extracted and tested for Zaire and Sudan ebolavirus and marburgvirus  using real-time quantitative RT-PCR assays designed to detect all known strains of each respective virus species [the primers/probe for the Sudan ebolavirus assay were EboSudBMG 1(+) 5′-GCC ATG GIT TCA GGT TTG AG-3′, EboSudBMG 1(−) 5′-GGT IAC ATT GGG CAA CAA TTC A and EboSudBMG Probe 5′FAM-AC GGT GCA CAT TCT CCT TTT CTC GGA-BHQ1]; (iv) the conventional RT-PCR was performed with the filo A/B primer set as previously described  using Superscript III (Invitrogen) according to the manufacturer's instructions. The specimen 200706291 was selected as the reference sample for further sequence analysis.
Pyro-sequencing was carried out utilizing the approach developed by 454 Life Sciences, and the method described by Cox-Foster et al., . Subsequent virus whole genome primer walking was performed as previously described  but using the primers specific for Bundibugyo ebolavirus RT-PCR amplification. In total, the entire virus genome was amplified in six overlapping RT-PCR fragments (all primers listed 5′ to 3′): fragment A (predicted size 2.7 kb) was amplified using forward-GTGAGACAAAGAATCATTCCTG with reverse-CATCAATTGCTCAGAGATCCACC; fragment B (predicted size 3.0 kb) was amplified using forward-CCAACAACACTGCATGTAAGT with reverse-AGGTCGCGTTAATCTTCATC; fragment C (predicted size 3.5 kb) was amplified using forward-GATGGTTGAGTTACTTTCCGG with reverse-GTCTTGAGTCATCAATGCCC; fragment D (predicted size 3.1 kb) was amplified using forward-CCACCAGCACCAAAGGAC with reverse-CTATCGGCAATGTAACTATTGG; fragment E (predicted size 3.4 kb) was amplified using forward-GCCGTTGTAGAGGACACAC with reverse-CACATTAAATTGTTCTAACATGCAAG and fragment F (predicted size 3.5 kb) was amplified using forward-CCTAGGTTATTTAGAAGGGACTA with reverse-GGT AGA TGT ATT GAC AGC AAT ATC.
The exact 5′ and 3′ ends of Bundibugyo ebolavirus were determined by 3′ RACE from virus RNA extracted from virus infected Vero E6 cell monolayers using TriPure isolation reagent. RNAs were then polyadenylated in vitro using A-Plus poly(A) polymerase tailing kit (Epicenter Biotechnologies) following the manufacturer's instructions and then purified using an RNeasy kit (Qiagen) following standard protocols. Ten microliters of in vitro polyadenylated RNA were added as template in RT-PCR reactions, using SuperScript III One-Step RT-PCR system with Platinum Taq High Fidelity (Invitrogen) following the manufacturer's protocol. Two parallel RT-PCR reactions using the oligo(dT)-containing 3′RACE-AP primer (Invitrogen) mixed with 1 of 2 viral specific primers, Ebo-U 692(−) ACAAAAAGCTATCTGCACTAT and Ebo-U18269(+) CTCAGAAGCAAAATTAATGG, generated ~700 nt long fragments containing the 3′ ends of either genomic and antigenomic RNAs. The resulting RT-PCR products were analyzed by agarose electrophoresis, and DNA bands of the correct sizes were purified using QIAquick Gel Extraction Kit (Qiagen) and sequenced using standard protocols (ABI).
The nucleotide sequence of the Côte d'Ivoire ebolavirus isolate RNA was initially determined using the exact same pyro-sequencing strategy as that used for Bundibugyo ebolavirus described above. This method generated sequence for approximately 70% of the entire genome. This draft sequence was then used to design a whole genome primer walking strategy for filling any gaps and confirming the initial sequence. The following Côte d'Ivoire ebolavirus-specific primers were used to generate RT-PCR fragments, designated A–F, as follows: Fragment A (predicted size 3.0 kb) was amplified using forward-GTGTGCGAATAACTATGAGGAAG and reverse-GTCTGTGCAATGTTGATGAAGG; Fragment B (predicted size 3.2 kb) was amplified using forward-CATGAAAACCACACTCAACAAC and reverse-GTTGCCTTAATCTTCATCAAGTTC; Fragment C (predicted size 3.0 kb) was amplified using forward-GGCTATAATGAATTTCCTCCAG and reverse-CAAGTGTATTTGTGGTCCTAGC; fragment D (predicted size 3.5 kb) was amplified using forward-GCTGGAATAGGAATCACAGG and reverse-CGGTAGTCTACAGTTCTTTAG; fragment E (predicted size 4.0 kb) was amplified using forward-GACAAAGAGATTAGATTAGCTATAG and reverse-GTAATGAGAAGGTGTCATTTGG; fragment F (predicted size 2.9 kb) was amplified using forward-CACGACTTAGTTGGACAATTGG and reverse-CAGACACTAATTAGATCTGGAAG; fragment G (predicted size 1.3 kb) was amplified using forward-CGGACACACAAAAAGAAWRAA and reverse-CGTTCTTGACCTTAGCAGTTC; and fragment H (predicted size 2.5 kb) was amplified using forward-GCACTATAAGCTCGATGAAGTC and reverse-TGGACACACAAAAARGARAA. A gap in the sequence contig was located between fragments C and D and this was resolved using the following primers to generate a predicted fragment of 1.5 kb: forward-CTGAGAGGATCCAGAAGAAAG and reverse-GTGTAAGCGTTGATATACCTCC. The terminal ~20 nucleotides of the sequence were not experimentally determined but were inferred by comparing with the other known Ebola genome sequences.
Bundibugyo ebolavirus real-time RT-PCR assay
The primers and probe used in the Bundibugyo ebolavirus specific Q-RT-PCR assay were as follows: EboU965(+): 5′-GAGAAAAGGCCTGTCTGGAGAA-3′, EboU1039(−): 5′-TCGGGTATTGAATCAGACCTTGTT-3′ and EboU989 Prb: 5′Fam-TTCAACGACAAATCCAAGTGCACGCA-3′BHQ1. Q-RT-PCR reactions were set up using Superscript III One-Step Q-RT-PCR (Invitrogen) according to the manufacturer's instructions and run for 40 cycles with a 58°C annealing temperature.
Modeltest 3.7  was used to examine 56 models of nucleotide substitution to determine the model most appropriate for the data. The General Time Reversible model incorporating invariant sites and a gamma distribution (GTR+I+G) was selected using the Akaike Information Criterion (AIC). Nucleotide frequencies were A = 0.3278, C = 0.2101, G = 0.1832, T = 0.2789, the proportion of invariant sites = 0.1412, and the gamma shape parameter = 1.0593. A maximum likelihood analysis was subsequently performed in PAUP*4.0b10  using the GTR+I+G model parameters. Bootstrap support values were used to assess topological support and were calculated based on 1,000 pseudoreplicates .
In addition, a Bayesian phylogenetic analysis was conducted in MrBayes 3.2  using the GTR+I+G model of nucleotide substitution. Two simultaneous analyses, each with four Markov chains, were run for 5,000,000 generations sampling every 100 generations. Prior to termination of the run, the AWTY module was used to assess Markov Chain Monte Carlo convergence to ensure that the length of the analysis was sufficient . Trees generated before the stabilization of the likelihood scores were discarded (burn in = 40), and the remaining trees were used to construct a consensus tree. Nodal support was assessed by posterior probability values (≥95 = statistical support).
We thank the initial investigation team which included Joseph F. Wamala, Omoruto D. Atai, Zirus William Lali, Kagirita Atek, Enoc Bahati, James Ndyezuka and Monica M. Musenero, whose efforts resulted in collection of the specimens tested here. We would also like to acknowledge Drs. Michael Egholm and Stephen K. Hutchison at 454 Life Sciences for their critical expertise with pyro-sequencing and their willingness to quickly make available essential resources in response to a public health crisis. In addition, we thank Debi Cannon and Zachary Reed for expert technical assistance as well as Craig Manning for help with computer graphics. We are grateful to Dr. Edward Katongole-Mbidde, the Director of UVRI, for providing a safe environment for specimen processing, and Beryl West for logistic and administrative support. The findings and conclusions in this report are those of the authors and do not necessarily represent the views of the funding agencies.
Conceived and designed the experiments: JST TKS MLK CGA STN. Performed the experiments: JST TKS MLK CGA SAR BB JK JAC. Analyzed the data: JST TKS MLK CGA SC SAR BB JK JAC PER TGK STN. Contributed reagents/materials/analysis tools: JST TKS MLK CGA SC SAR PLQ WIL RD JWT SO JL BB JK JAC PER TGK STN. Wrote the paper: JST TKS CGA SAR STN.
Suzuki Y, Gojobori T (1997) The origin and evolution of Ebola and Marburg viruses. Mol Bio Evol 14(
: 800–806. Find this article online
Sanchez A, Geisbert TW, Feldmann H (2007) Filoviridae: Marburg and Ebola Viruses. In: Knipe DM, Howley PM, editors. Fields Virology. Philidelphia: Lippincott Williams and Williams. pp. 1409–1448.
Leroy EM, Kumulungui B, Pourrut X, Rouquet P, Hassanin A, et al. (2005) Fruit bats as reservoirs of Ebola virus. Nature 438: 575–6. Find this article online
Towner JS, Pourrut X, Albariño CG, Nze Nkogue C, Bird BH, et al. (2007) Marburg virus infection detected in a common African bat. PLoS ONE 2(
: e764. doi:10.1371/journal.pone.0000764. Find this article online
Swanepoel R, Smit SB, Rollin PE, Formenty P, Leman PA, et al. (2007) Studies of reservoir hosts for Marburg virus. Emerg Infect Dis 13(12): 1847–51. Find this article online
Le Guenno B, Formenty P, Wyers M, Gounon P, Walker F, et al. (1995) Isolation and partial characterization of a new strain of Ebola virus. Lancet 345: 1271–4. Find this article online
Ksiazek TG, Rollin PE, Williams AJ, Bressler DS, Martin ML, et al. (1999) Clinical virology of Ebola hemorrhagic fever (EHF): virus, virus antigen, IgG and IgM antibody findings among EHF patients in Kikwit, Democratic Republic of Congo, 1995. J Infect Dis 179: (suppl 1)S177–S187. Find this article online
Cox-Foster DL, Conlan S, Holmes EC, Palacios G, Evans JD, et al. (2007) A metagenomic survey of microbes in honey bee colony collapse disorder. Science 318: 283–7. Find this article online
World Health Organization (2008) Ebola outbreak contained in Uganda. Features. Available: http://www.who.int/features/2008/ebola_outbreak/en/
. Accessed 22 February 2008.
Jones SM, Feldmann H, Ströher U, Geisbert JB, Fernando L, et al. (2005) Live attenuated recombinant vaccine protects nonhuman primates against Ebola and Marburg viruses. Nat Med 11: 786–90. Find this article online
Geisbert TW, Daddario-DiCaprio KM, Williams KJ, Geisbert JB, Leung A, et al. (2008) Recombinant vesicular stomatitis virus vector mediates postexposure protection against Sudan Ebola hemorrhagic fever in nonhuman primates. J Virol 82: 5664–8. Find this article online
Sullivan NJ, Sanchez A, Rollin PE, Yang ZY, Nabel GJ (2000) Development of a preventive vaccine for Ebola virus infection in primates. Nature 408: 605–609. Find this article online
Ksiazek TG, West CP, Rollin PE, Jahrling PB, Peters C J (1999) ELISA for the detection of antibodies to Ebola viruses. J Infect Dis 179: (suppl 1)S192–S198. Find this article online
Rodriguez L, De Roo A, Guimard Y, Trappier SG, Sanchez A, et al. (1999) Persistence and genetic stability of Ebola virus during the outbreak in Kikwit, Democratic Republic of Congo, 1995. J Infect Dis 179: (suppl 1)S170–S176. Find this article online
Towner JS, Sealy TK, Ksiazek TG, Nichol ST (2007) High-throughput molecular detection of hemorrhagic fever virus threats with applications for outbreak settings. J Infect Dis 196: (suppl 2)S205–212. Find this article online
Sanchez A, Ksiazek TG, Rollin PE, Miranda MEG, Trappier SG, et al. (1999) Detection and molecular characterization of Ebola viruses causing disease in human and nonhuman primates. J Infect Dis 179: (suppl 1)S164–S169. Find this article online
Towner JS, Khristova ML, Sealy TK, Vincent MJ, Erickson BR, et al. (2006) Marburgvirus genomics and association with a large hemorrhagic fever outbreak in Angola. J Virol 80: 6497–516. Find this article online
Posada D, Crandall KA (1998) MODELTEST: testing the model of DNA substitution. Bioinformatics 14: 817–818. Find this article online
Swofford DL (2002) PAUP*: phylogenetic analysis using parsimony (*and other methods) version 4.0b10. Sunderland, Mass: Sinauer Associates.
Felsenstein J (1985) Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39: 783–791. Find this article online
Ronquist F, Huelsenbeck JP (2003) MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574. Find this article online
Nylander JAA, Wilgenbusch JC, Warren DL, Swofford DL (2008) AWTY (are we there yet?): a system for graphical exploration of MCMC convergence in Bayesian phylogenetics. Bioinformatics 24: 581–583. Find this article online
All site content, except where otherwise noted, is licensed under a Creative Commons Attribution License.
RC 0.9 beta Managed Colocation provided by UnitedLayer.
TG-1Strong supporters of the PLoS.