Skip directly to site content Skip directly to page options Skip directly to A-Z link Skip directly to A-Z link Skip directly to A-Z link
Volume 21, Number 7—July 2015

Monitoring of Ebola Virus Makona Evolution through Establishment of Advanced Genomic Capability in Liberia

Jeffrey R. Kugelman1, Michael R. Wiley1, Suzanne Mate1, Jason T. Ladner1, Brett Beitzel, Lawrence Fakoli, Fahn Taweh, Karla Prieto, Joseph W. Diclaro, Timothy Minogue, Randal J. Schoepp, Kurt E. Schaecher, James Pettitt, Stacey Bateman, Joseph Fair, Jens H. Kuhn, Lisa Hensley, Daniel J. Park, Pardis C. Sabeti, Mariano Sanchez-Lockhart, Fatorma K. Bolay, Gustavo PalaciosComments to Author , on behalf of US Army Medical Research Institute of Infectious DiseasesNational Institutes of HealthIntegrated Research Facility–Frederick Ebola Response Team 2014–2015
Author affiliations: US Army Medical Research Institute of Infectious Diseases, Fort Detrick, Frederick, Maryland, USA (J.R. Kugelman, M.R. Wiley, S. Mate, J.T. Ladner, B. Beitzel, K. Prieto, T. Minogue, R.J. Schoepp, K.E. Schaecher, S. Bateman, M. Sanchez-Lockhart, G. Palacios); Liberian Institute for Biomedical Research, Charlesville, Liberia (L. Fakoli, F. Taweh, F.K. Bolay); Naval Medical Research Unit 3, Cairo, Egypt (J.W. Diclaro); Integrated Research Facility at Fort Detrick, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Frederick (J. Pettitt, J.H. Kuhn, L. Hensley); Foundation Merieux, Washington, DC, USA (J. Fair); Broad Institute, Cambridge, Massachusetts, USA (D.J. Park, P.C. Sabeti); Center for Systems Biology, Harvard University, Cambridge (P.C. Sabeti)

Cite This Article


To support Liberia’s response to the ongoing Ebola virus (EBOV) disease epidemic in Western Africa, we established in-country advanced genomic capabilities to monitor EBOV evolution. Twenty-five EBOV genomes were sequenced at the Liberian Institute for Biomedical Research, which provided an in-depth view of EBOV diversity in Liberia during September 2014–February 2015. These sequences were consistent with a single virus introduction to Liberia; however, shared ancestry with isolates from Mali indicated at least 1 additional instance of movement into or out of Liberia. The pace of change is generally consistent with previous estimates of mutation rate. We observed 23 nonsynonymous mutations and 1 nonsense mutation. Six of these changes are within known binding sites for sequence-based EBOV medical countermeasures; however, the diagnostic and therapeutic impact of EBOV evolution within Liberia appears to be low.

The outbreak of Ebola virus disease (EVD) in Western Africa that started in November 2013 (1) is the largest recorded filovirus disease outbreak. As the outbreak continues, public health and emerging infectious disease officials have declared a continuing need for real-time monitoring of Ebola virus (EBOV) evolution (2,3). As of March 11, 2015, a total of 41% of reported cases had been fatal (4). By the end of March 2015, the intensity of the outbreak, which throughout its course affected 6 Western Africa countries, appeared to be receding, with near 0 activity in Liberia and no cases in Mali, Nigeria, and Senegal. However, EBOV continues to spread in Guinea and Sierra Leone. The epidemic is still causing more infections per week than have been recorded in previous EVD outbreaks (5). Therefore, public health officials continue to use media to maintain public awareness, to advocate for diligent handwashing and use of other protective measures, and to avoid complacency that could lead to reemergence (5). Vigilance is of paramount importance because currently used assays for EVD diagnosis, and many medical countermeasures in development, were designed using EBOV reference genome variants from previous outbreaks (69). Therefore, monitoring EBOV genomic drift is crucial because genetic changes can affect the efficacy of sequence-based therapeutics and diagnostics.

The size and spread of the current EVD outbreak reinforces the need to build public health infrastructure, including state-of-the-art diagnostic and surveillance capabilities, to implement and maintain effective EVD monitoring, treatment, and prevention platforms. The Liberian Institute for Biomedical Research (LIBR), established in 1975, is located in Charlesville, 50 km southeast of Liberia’s capital, Monrovia. As of April 2, 2015, it is one of the few local facilities within Liberia processing clinical samples from persons suspected to have EVD. A consortium comprising US government and nongovernment agencies has been working with the Liberian government to equip LIBR with advanced genomic sequencing capabilities. These capabilities are dedicated primarily to EVD surveillance activities, including genome sequencing of EBOV-positive samples. The new LIBR Genome Center has a Miseq sequencer (Illumina, San Diego, CA, USA) and ancillary supporting capabilities, including electrophoresis for qualification, fluorometry for quantitation, PCR for amplification, and fully functional computational analysis capabilities to perform pathogen discovery and microbial genome characterization. The US Army Medical Research Institute of Infectious Diseases (USAMRIID) Center for Genome Sciences supports LIBR operation and development. Sample preparation procedures under biosafety containment are provided within the same building complex by the Liberian National Reference Laboratories, operated by USAMRIID and the National Institutes of Health Integrated Research Facility Ebola Response Team (Fort Detrick, Frederick, MD, USA). Throughput at the LIBR Genome Center is 10–20 samples (≈10 billion bases of sequence data) per week, with a target turnaround time of 7 days from sample receipt for high-priority samples. To ensure long-term sustainment of surveillance-based sequencing capabilities, local biomedical scientists have been trained and can proficiently perform all daily activities.

Here we demonstrate the utility and capabilities of the LIBR Genome Center. With the immediate goal of continuing the natural history characterization of the EBOV Makona variant (EBOV/Mak [10]) currently circulating in Western Africa and to support ongoing clinical trials to evaluate candidate medical countermeasures, we describe 25 EBOV genome sequences from the first 5 sequencing runs conducted at the LIBR Genome Center. We chose these samples for full-genome characterization from ≈1,700 available samples on the basis of high viral load (cycle threshold [Ct]) value <24) and date of collection to ensure up-to-date temporal coverage.

Materials and Methods


Figure 1

Thumbnail of A) Median-joining haplotype network constructed from a full-genome alignment of 122 clinical Ebola virus Makona (EBOV/Mak) isolates (list of isolates in Technical Appendix 3). Each colored vertex represents a sampled viral haplotype, with the numbered vertices representing the centers of the 3 clusters described in (12). All sampled isolates from Liberia originated from cluster 2. The size of each vertex is relative to the number of sampled isolates, and the colors indicate country

Figure 1. A) Median-joining haplotype network constructed from a full-genome alignment of 122 clinical Ebola virus Makona (EBOV/Mak) isolates (list of isolates in Technical Appendix 3). Each colored vertex...

We chose samples from 25 patients from the larger collection (≈1,700 positive cases) on the basis of diagnostic Ct values that indicated a high enough viral load to provide a full genome (Ct<24), beginning with the most recent available at the time of preparation in February 2015. Sampling continued with progressively older samples to describe the lineages most likely to still be circulating at the time. These patients were treated in 7 different Ebola treatment units and had resided in 7 of the 15 counties in Liberia (Table 1; Technical Appendix 1 Figure 1). Plasma or oral swab samples from which viral RNA was recovered and sequenced were tested at LIBR during September 23, 2014–February 14, 2015. Patients’ ages were as follows: 1 infant (1 year), 6 children (2–15 years), 8 young adults (18–35 years), and 10 middle-aged adults (42–67 years). The male:female ratio was 2:1. However, among ≈1,700 samples at LIBR from persons with EVD, the ratio was close to 1:1 (48%/52%), and viral load did not differ by patient sex, which demonstrates that our higher ratio is a sampling artifact.

Sample Processing

RNA was converted to cDNA and amplified by using sequence-independent single-primer amplification (11). Amplified cDNA was quantified with a Qubit 3.0 fluorometer (Life Technologies, Carlsbad, CA, USA) and used as the starting material for the Illumina Nextera XT DNA library preparation kit (Illumina). Sequencing was performed on an Illumina Miseq by using either V2 or V3 reagent kits (Illumina) with a minimum of 2 × 151 cycles per run.

Genome Assembly

We assembled EBOV genomes by aligning reads to the genome of Ebola virus/H.sapiens-wt/SLE/2014/Makona-G3686.1 (GenBank accession no. KM034562.1) (12). Amplification primers were removed from the sequencing reads by using Cutadapt version 1.21 (13), and low-quality reads/bases were filtered by using Prinseq-lite version 0.20.4 (-min_qual_mean 25 -trim_left 20 -min_len 50) (14). Reads were aligned to the reference genome by using DNAStar Lasergene nGen (DNAStar, Madison, WI, USA), and a new consensus was generated by using a combination of Samtools v0.1.18 (15) and custom scripts. Only bases with Phred quality score >20 were used in consensus calling, and a minimum of 3× read-depth coverage, in support of the consensus, was required to make a call; positions lacking this depth of coverage were treated as missing (i.e., called as “N”).

Genetic Analysis

Consensus sequences generated here were aligned with additional publically available EBOV genomes by using Sequencher version 5.2.3 (Gene Codes, Ann Arbor, MI, USA). SnpEff version 4.1b (build 2015-02-13) was used to annotate all single-nucleotide polymorphisms (SNPs) by using the genome of Ebola virus/H.sapiens-wt/GIN/2014/Makona-C15 (GenBank accession no. KJ660346.2) as a reference (16). All 25 genomes from Liberia were used to identify variable sites. For the rest of the genetic analysis, we used only the 14 sequences with >90% genome coverage. A median-joining haplotype network was constructed in PopART version 1.7.2 ( Path-O-Gen version 1.4 (17) was used to calculate the root-to-tip distances by using a maximum-likelihood phylogeny (PhyML version 3.0 (18); general time reversible model) with rooting based on the EBOV phylogeny published by Gire et al. (12). BEAST version 1.8.2 (17) was used to estimate the mutation rate and the time to the most recent common ancestor for several evolutionary lineages that included Liberia EBOV isolates. For analysis, we divided the alignment into 3 partitions (i.e., first + second codon sites, third codon site, and noncoding sites). The substitution process was modeled independently for each by using the Hasegawa, Kishino, and Yano model with 4 gamma categories. An exponential growth coalescent model was used with a strict clock. The XML input file is available on request from the authors.


From the first 5 sequencing runs, we obtained 25 EBOV genomes with >50% coverage; 6 of these were coding complete (Table 2) (19). These genomes contained 97 new sequence variants: 47 synonymous, 23 nonsynonymous, 1 nonsense, and 26 noncoding mutations (Technical Appendix 2). Multiple distinct evolutionary lineages were detected, but all were consistent with a single introduction of a cluster 2–type (12) virus into Liberia followed by within-country diversification (Figure 1, panel A). Because 19 of the 25 genomes had calls at all 5 positions that discriminate clusters 1, 2, and 3, we have high confidence in cluster attribution.

Molecular dating places the common ancestor to all of the sampled isolates from Liberia during May 2–July 9, 2014 (95% highest posterior density [HPD] interval), which corresponds with the early days of the outbreak in Monrovia (3). However, we cannot rule out ongoing EBOV exchange among EVD-infected countries. In fact, shared ancestry among 3 isolates from Liberia and the 4 available sequences from Mali suggests some level of international movement. We estimated dates associated with 2 nodes along the shared Liberia/Mali EBOV lineage (labeled * and ** in Figure 1, panel A); these estimates ranged from July 6 through September 15, 2014, and from July 26 through September 27, 2014, respectively (95% HPD). Overall, collection dates correlated well with root-to-tip distances within the Western Africa EVD outbreak (Figure 1, panel B). Linear regression analysis (using the lm function in R version 3.1.1; estimated an overall rate of change of 9.17 × 10−4 substitutions/site/year (± 5.23 × 10−5). Bayesian analysis estimated a similar rate of change of 9.44–15.67 × 10−4 substitutions/site/year (95% HPD).

We reviewed all publicly available genomic information for EBOV/Mak (122 genome sequences [1,12]) to evaluate the effect of genomic drift on biomedical countermeasures (drugs and diagnostic assays). We assessed the potential impact of intra-outbreak genetic divergence on 13 drugs and 2 diagnostic assays (known to be used in Liberia) with the same approach previously used (6). Two sequence-binding treatment modalities are available for postexposure treatment of EVD: small interfering RNAs (siRNAs) (20) and phosphorodiamidate morpholino oligomers (21) targeting L, VP24, and/or VP35 gene transcripts, and passive immunotherapy based on antibodies or antibody cocktails targeting EBOV glycoprotein (2226). These treatments inhibit viral replication by targeting viral transcripts for degradation (siRNA) or by blocking translation (phosphorodiamidate morpholino oligomers), or they acutely neutralize the virus to enable the host to mount an effective immune response (passive immunotherapy). These countermeasures were originally designed specifically against sequences obtained during previous outbreaks (20,27) or were generated against their glycoproteins (e.g., the monoclonal antibodies [mAbs] were obtained after immunization with Ebola virus/H.sapiens-tc/COD/1995/Kikwit-9510621 [EBOV/Kik-9510621] [28]).

Figure 2

Thumbnail of Mutation analysis of candidate therapeutic drug and diagnostic binding sites used in outbreak of Ebola virus (EBOV) disease, Western Africa. A single-nucleotide polymorphism (SNP) table is combined with a heat map based on 2 categories: 1) mutations tolerated by the therapeutic drug or diagnostic target (highlighted in green); 2) mutations within the binding region of a therapeutic drug or diagnostic assay that have not yet been tested (highlighted in yellow/orange) (20–24,27,30,31)

Figure 2. Mutation analysis of candidate therapeutic drug and diagnostic binding sites used in outbreak of Ebola virus (EBOV) disease, Western Africa. A single-nucleotide polymorphism (SNP) table is combined with a heat map...

Since the Western Africa outbreak began, at least 33 viral mutations have occurred that could affect countermeasures. We previously reported 27 of these mutations (6). Twenty-six (79%) mutations induced nonsynonymous changes to epitopes recognized by mAbs included in passive immunotherapy cocktails. Another 5 (15%) were located in published binding regions of siRNA-based therapeutic drugs. Tekmira has adjusted its siRNAs to account for 4 of these 5 changes since its initial publication (29; E.P. Thi et al., unpub. data). The final 2 mutations were located in the published binding region of primers or probes for quantitative PCR diagnostic tests that have been used during outbreak control activities in Liberia: 1 change each in the binding sites of the Kulesh-TM assay and the Kulesh-MGB assay (9). Nevertheless, reassessment of the assays at USAMRIID has suggested that the changes will be tolerated without loss in sensitivity (data not shown). Changes in all EBOV/Mak sequences are considered “interoutbreak” (n = 23); changes observed only in some sequences from Western Africa are considered “intraoutbreak” sites (n = 10, EBOV-WA <100%). We also examined the binding sites of an additional 18 publicly available EBOV quantitative PCRs, which might (or might not) also be used in Western Africa (Technical Appendix 1 Figure 2, Technical Appendix 1 Table). We observed 25 changes, of which 6 were reported previously (12). Each SNP has the potential to affect the efficacy of available therapeutic drugs (original and updated versions) or diagnostic assays (Table 3; Figure 2; Technical Appendix 1 Figure 2, Technical Appendix 1 Table; nucleotide positions are reported relative to EBOV/Kik-9510621, for consistency [6]).

Several of the 27 previously identified changes (green in Figure 2) already have been demonstrated to be tolerated while maintaining efficacy (24,30,3234), thus minimizing their potential effect (6). Six of these 33 SNPs (EBOV-LIB <100%; orange in Figure 2) appeared during the surveillance period of this study (September 23, 2014–February 14, 2015) in samples obtained in Liberia (12). None of these changes have been previously associated with EBOV resistance to any therapeutic drug. Five of the new changes might affect 1 of the components of the ZMapp antibody cocktail (mAb 13C6). However, the conformational target site for this antibody (positions 1–295, soluble glycoprotein) is broader in length and more poorly defined than the other sequence-based countermeasure targets considered in our risk assessment. The sixth mutation might affect the binding site of the siRNA viral protein (VP) 35 target (for that particular sample, the mutation appears in an area of low sequencing coverage depth). Thus, when these new changes are combined with the changes observed previously (yellow in Figure 2), we can conclude that retesting several therapeutic drugs against isolates currently circulating might be necessary to determine whether any of these mutations impact their efficacy. In particular, it is important to reevaluate drugs that include mAb 13C6 (part of the ZMapp, ZMAb, and MB-003 antibody cocktails), mAb 13F6 (part of MB-003), mAb 1H3 (part of ZMAb), and the siRNA VP35 targets (Table 3, Figure 2) (6).


Our study details the establishment of a genomic sequencing and analysis center within Liberia for real-time monitoring of viral evolution. The initial sequences generated at this facility have provided a first glimpse into EBOV/Mak evolution from the end of 2014 to the beginning of 2015. Although genetically diverse, the viruses circulating in Liberia during this period are consistent with a single introduction event followed by diversification within Liberia. The cluster 2 haplotype from which all the sampled Liberia sequences radiate is thought to have been circulating in Guinea and Sierra Leone during late May 2014 (12). Moreover, it was the second most common sequence detected in Sierra Leone during late May through mid-June (12). Introduction of this haplotype from either of these neighboring countries could have resulted in the sampled diversity; however, we cannot rule out the possibility of multiple introductions. Additional spatial and temporal sampling within Liberia, Guinea, and Sierra Leone will help to differentiate these 2 scenarios.

The 25 Liberia EBOV/Mak genomes included 23 nonsynonymous mutations and 1 nonsense mutation that have not previously been seen in Western Africa (although some of these mutations have been observed in EBOV isolates from previous EVD outbreaks). A nonsense mutation, which is present within 2 of the 25 sequences, is predicted to result in premature truncation (6 aa) of VP30. VP30 is an essential protein for viral transcription; it is needed for the RNA-dependent RNA polymerase (L) to read beyond a cis-RNA element in the nucleoprotein mRNA 5′ untranslated region (35) and is required to reinitiate transcription at gene junctions (36). Moreover, VP30 phosphorylation modulates the composition and function of the RNA synthesis machinery (37). To our knowledge, no functional domains have been described in the truncated region. Further characterization is needed to determine whether this or any of the other detected mutations impacted the relative fitness of the affected EBOV isolates. Within Liberia, geography showed little correlation with phylogeny; most EBOV lineages within Liberia appear to be geographically widespread within the sampled regions.

Previous analysis of EBOV/Mak genomes from Sierra Leone and Guinea suggests that the evolutionary rate within the current EVD outbreak might be higher than the rate between outbreaks (12). After incorporation of sequences from Liberia, which were collected later in the outbreak, our estimates of substitution rate fell between the previous estimates for EBOV/Mak only and for all EBOV (12,38). As more sequence data become available, it will be interesting to see whether a significant change in the evolutionary rate can be detected within the current EVD outbreak.

Our ability to quantify international EBOV exchange is limited because few isolates from other countries were available during the sampled timeframe. However, shared ancestry between isolates from Mali and 3 isolates from Liberia suggests at least 1 transmission event across national borders (3). All EVD cases in Mali have been attributed to movement of infected persons into Mali from Guinea (39). With the current dataset, it is impossible to say whether the shared Liberia/Mali lineage originated in Liberia and was then transported to Mali through Guinea or whether the lineage emerged in Guinea and later moved independently to Liberia and Mali. Active EBOV outbreaks were occurring in both Liberia and Guinea during the period estimated for the emergence of this shared lineage (July–September 2014).

The genomic changes observed for EBOV/Mak during its circulation in Liberia append 5 additional mutations to the list of changes that might affect the binding of the 13C6 mAb, a component of ZMapp. All of these changes, however, were present at relatively low frequency (<12%) in our current sample, and none of the sampled lineages have accumulated >1 change per therapeutic drug type. We observed no significant changes (i.e., likely to affect efficacy) in the binding sites for the 2 diagnostic assays known to be used in Liberia. Overall, no dramatic changes were observed in the samples evaluated; the risk assessment for the impact of genomic drift during the outbreak should remain low. As previously stated (6), our analysis is not without caveats. Our current analysis covers only the late period of the outbreak in Liberia; no analysis has yet been published with data for similar time points from Guinea or Sierra Leone. In addition, to complete our assessment of the evolution of EBOV in Liberia, an earlier period of time from the introduction of the virus in March 2014 to early September 2014 needs to be investigated.

Our findings offer a concise evaluation of the potential impact of the evolution of EBOV/Mak based on genome reconstruction of 25 isolates from Liberia obtained during September 2014–February 2015. This work would not have been possible without the establishment of a genomic surveillance capability in Liberia, which emphasizes the need for global sequencing capabilities to be part of the first response during future virus outbreaks.

Dr. Kugelman is a biodefense research scientist, computational biologist, and head of Bioinformatics at the Center for Genome Sciences at USAMRIID. His research interests include the genomic study of filovirus and orthopoxvirus infections.



We thank Matthew A. Voorhees, Brian J. Kearney, Gary W. Carter, Jason Kindrachuk, Anthony R. Jones, Alex Hail, Alissa L. Byrne, Melissa A. Dugan, Mark M. Bailey, Krisztina Janosko, Joanna Fishback, Joshua Johnson, Richard Bennett, Susana Padilla, Michelle A. Jefferson, Jonathan Marchand, and Michael R. Holbrook for providing the infrastructure, isolation of viral genomic material, and Ebola diagnostic screening on which this study is built. We thank Aaron T. Momolu, Josiah S. George, Charlesetta Kanneh, Vera Yatta Walker, and Augustine Fallah for their help, support, and enthusiasm. We thank Neal Woollen, Dave Norwood, and Peter Jahrling for their full support on this endeavor. We also thank Deborah R. Malac for her support of the USAMRIID and NIH/IRF Ebola Response Team; Tabitha Austill for arranging the shipment of equipment and supplies to LIBR; and Sean Lovett and Brad Pfeffer for the expedited deposition of sequence data in public repositories. We also thank Laura Bollinger for providing technical writing services.

This work was supported by Defense Threat Reduction Agency, the Global Biosurveillance Technology Initiative, the US Agency for International Development, and Illumina.

J.H.K. performed this work as an employee of Tunnell Government Services, Inc., a subcontractor to Battelle Memorial Institute under its prime contract with the National Institute of Allergy and Infectious Diseases, under contract no. HHSN272200700016I.



  1. Baize  S, Pannetier  D, Oestereich  L, Rieger  T, Koivogui  L, Magassouba  N, Emergence of Zaire Ebola virus disease in Guinea—preliminary report. N Engl J Med. 2014;371:141825 . DOIPubMedGoogle Scholar
  2. Hawkes  N. Ebola outbreak is a public health emergency of international concern, WHO warns. BMJ. 2014;349:g5089. DOIPubMedGoogle Scholar
  3. Yozwiak  NL, Schaffner  SF, Sabeti  PC. Data sharing: make outbreak research open access. Nature. 2015;518:4779. DOIPubMedGoogle Scholar
  4. World Health Organization. Ebola situation report—11 March 2015 [cited 2015 Mar 30].
  5. Chan  M. WHO report by the Director-General to the special session of the Executive Board on Ebola [cited 2015 Mar 30].
  6. Kugelman  JR, Sanchez-Lockhart  M, Andersen  KG, Gire  S, Park  DJ, Sealfon  R, Evaluation of the potential impact of Ebola virus genomic drift on the efficacy of sequence-based candidate therapeutics. MBio. 2015;6:e0222714. DOIPubMedGoogle Scholar
  7. Shuchman  M. WHO enters new terrain in Ebola research. CMAJ. 2014;186:E5278. DOIPubMedGoogle Scholar
  8. Sayburn  A. WHO gives go ahead for experimental treatments to be used in Ebola outbreak. BMJ. 2014;349:g5161. DOIPubMedGoogle Scholar
  9. Trombley  AR, Wachter  L, Garrison  J, Buckley-Beason  VA, Jahrling  J, Hensley  LE, Comprehensive panel of real-time TaqMan polymerase chain reaction assays for detection and absolute quantification of filoviruses, arenaviruses, and New World hantaviruses. Am J Trop Med Hyg. 2010;82:95460. DOIPubMedGoogle Scholar
  10. Kuhn  JH, Andersen  KG, Baize  S, Bao  Y, Bavari  S, Berthet  N, Nomenclature- and database-compatible names for the two Ebola virus variants that emerged in Guinea and the Democratic Republic of the Congo in 2014. Viruses. 2014;6:476099. DOIPubMedGoogle Scholar
  11. Djikeng  A, Halpin  R, Kuzmickas  R, Depasse  J, Feldblyum  J, Sengamalay  N, Viral genome sequencing by random priming methods. BMC Genomics. 2008;9:5. DOIPubMedGoogle Scholar
  12. Gire  SK, Goba  A, Andersen  KG, Sealfon  RS, Park  DJ, Kanneh  L, Genomic surveillance elucidates Ebola virus origin and transmission during the 2014 outbreak. Science. 2014;345:136972. DOIPubMedGoogle Scholar
  13. Martin  M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnetjournal. 2011;17:102.
  14. Schmieder  R, Edwards  R. Quality control and preprocessing of metagenomic datasets. Bioinformatics. 2011;27:8634. DOIPubMedGoogle Scholar
  15. Li  H, Handsaker  B, Wysoker  A, Fennell  T, Ruan  J, Homer  N, The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:20789. DOIPubMedGoogle Scholar
  16. Cingolani  P, Platts  A. Wang le L, Coon M, Nguyen T, Wang L, et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin). 2012;6:80–92.
  17. Drummond  AJ, Suchard  MA, Xie  D, Rambaut  A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:196973. DOIPubMedGoogle Scholar
  18. Guindon  S, Dufayard  JF, Lefort  V, Anisimova  M, Hordijk  W, Gascuel  O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59:30721. DOIPubMedGoogle Scholar
  19. Ladner  JT, Beitzel  B, Chain  PS, Davenport  MG, Donaldson  EF, Frieman  M, Standards for sequencing viral genomes in the era of high-throughput sequencing. MBio. 2014;5:e0136014. DOIPubMedGoogle Scholar
  20. Geisbert  TW, Lee  AC, Robbins  M, Geisbert  JB, Honko  AN, Sood  V, Postexposure protection of non-human primates against a lethal Ebola virus challenge with RNA interference: a proof-of-concept study. Lancet. 2010;375:1896905. DOIPubMedGoogle Scholar
  21. Warren  TK, Warfield  KL, Wells  J, Swenson  DL, Donner  KS, Van Tongeren  SA, Advanced antisense therapies for postexposure protection against lethal filovirus infections. Nat Med. 2010;16:9914. DOIPubMedGoogle Scholar
  22. Qiu  X, Audet  J, Wong  G, Fernando  L, Bello  A, Pillet  S, Sustained protection against Ebola virus infection following treatment of infected nonhuman primates with ZMAb. Sci Rep. 2013;3:3365.
  23. Pettitt  J, Zeitlin  L, Kim  DH, Working  C, Johnson  JC, Bohorov  O, Therapeutic intervention of Ebola virus infection in rhesus macaques with the MB-003 monoclonal antibody cocktail. Sci Transl Med. 2013;5:199ra113. PMID: 23966302
  24. Qiu  X, Wong  G, Audet  J, Bello  A, Fernando  L, Alimonti  JB, Reversion of advanced Ebola virus disease in nonhuman primates with ZMapp. Nature. 2014;514:4753 .PubMedGoogle Scholar
  25. Olinger  GG Jr, Pettitt  J, Kim  D, Working  C, Bohorov  O, Bratcher  B, Delayed treatment of Ebola virus infection with plant-derived monoclonal antibodies provides protection in rhesus macaques. Proc Natl Acad Sci U S A. 2012;109:180305. DOIPubMedGoogle Scholar
  26. Dye  JM, Herbert  AS, Kuehne  AI, Barth  JF, Muhammad  MA, Zak  SE, Postexposure antibody prophylaxis protects nonhuman primates from filovirus disease. Proc Natl Acad Sci U S A. 2012;109:50349. DOIPubMedGoogle Scholar
  27. Heald  AE, Iversen  PL, Saoud  JB, Sazani  P, Charleston  JS, Axtelle  T, Safety and pharmacokinetic profiles of phosphorodiamidate morpholino oligomers with activity against Ebola virus and Marburg virus: results of two single-ascending-dose studies. Antimicrob Agents Chemother. 2014;58:663947. DOIPubMedGoogle Scholar
  28. Kuhn  JH, Bao  Y, Bavari  S, Becker  S, Bradfute  S, Brister  JR, Virus nomenclature below the species level: a standardized nomenclature for natural variants of viruses assigned to the family Filoviridae. Arch Virol. 2013;158:30111. DOIPubMedGoogle Scholar
  29. Rezler  JP. Tekmira provides periodic update on TKM-Ebola program [cited 2015 Mar 30].
  30. Wilson  JA, Hevey  M, Bakken  R, Guest  S, Bray  M, Schmaljohn  AL, Epitopes involved in antibody-mediated protection from Ebola virus. Science. 2000;287:16646. DOIPubMedGoogle Scholar
  31. Qiu  X, Alimonti  JB, Melito  PL, Fernando  L, Stroher  U, Jones  SM. Characterization of Zaire ebolavirus glycoprotein-specific monoclonal antibodies. Clin Immunol. 2011;141:21827 . DOIPubMedGoogle Scholar
  32. Towner  JS, Paragas  J, Dover  JE, Gupta  M, Goldsmith  CS, Huggins  JW, Generation of eGFP expressing recombinant Zaire ebolavirus for analysis of early pathogenesis events and high-throughput antiviral drug screening. Virology. 2005;332:207. DOIPubMedGoogle Scholar
  33. Swenson  DL, Warfield  KL, Warren  TK, Lovejoy  C, Hassinger  JN, Ruthel  G, Chemical modifications of antisense morpholino oligomers enhance their efficacy against Ebola virus infection. Antimicrob Agents Chemother. 2009;53:208999. DOIPubMedGoogle Scholar
  34. Geisbert  TW, Hensley  LE, Kagan  E, Yu  EZ, Geisbert  JB, Daddario-DiCaprio  K, Postexposure protection of guinea pigs against a lethal Ebola virus challenge is conferred by RNA interference. J Infect Dis. 2006;193:16507. DOIPubMedGoogle Scholar
  35. Theriault  S, Groseth  A, Neumann  G, Kawaoka  Y, Feldmann  H. Rescue of Ebola virus from cDNA using heterologous support proteins. Virus Res. 2004;106:4350. DOIPubMedGoogle Scholar
  36. Martínez  MJ, Biedenkopf  N, Volchkova  V, Hartlieb  B, Alazard-Dany  N, Reynard  O, Role of Ebola virus VP30 in transcription reinitiation. J Virol. 2008;82:1256973. DOIPubMedGoogle Scholar
  37. Biedenkopf  N, Hartlieb  B, Hoenen  T, Becker  S. Phosphorylation of Ebola virus VP30 influences the composition of the viral nucleocapsid complex: impact on viral transcription and replication. J Biol Chem. 2013;288:1116574. DOIPubMedGoogle Scholar
  38. Carroll  SA, Towner  JS, Sealy  TK, McMullan  LK, Khristova  ML, Burt  FJ, Molecular evolution of viruses of the family Filoviridae based on 97 whole-genome sequences. J Virol. 2013;87:260816. DOIPubMedGoogle Scholar
  39. World Health Organization. Mali: details of the additional cases of Ebola virus disease [cited 2015 Mar 30].




Cite This Article

DOI: 10.3201/eid2107.150522

1These authors contributed equally to this article.

Table of Contents – Volume 21, Number 7—July 2015

EID Search Options
presentation_01 Advanced Article Search – Search articles by author and/or keyword.
presentation_01 Articles by Country Search – Search articles by the topic country.
presentation_01 Article Type Search – Search articles by article type and issue.



Please use the form below to submit correspondence to the authors or contact them at the following address:

Gustavo Palacios, Center for Genome Sciences, USAMRIID, 1425 Porter St. Fort Detrick, Frederick, MD 21702, USA

Send To

10000 character(s) remaining.


Page created: June 12, 2015
Page updated: June 12, 2015
Page reviewed: June 12, 2015
The conclusions, findings, and opinions expressed by authors contributing to this journal do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions. Use of trade names is for identification only and does not imply endorsement by any of the groups named above.