Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Real-time, portable genome sequencing for Ebola surveillance

Abstract

The Ebola virus disease epidemic in West Africa is the largest on record, responsible for over 28,599 cases and more than 11,299 deaths1. Genome sequencing in viral outbreaks is desirable to characterize the infectious agent and determine its evolutionary rate. Genome sequencing also allows the identification of signatures of host adaptation, identification and monitoring of diagnostic targets, and characterization of responses to vaccines and treatments. The Ebola virus (EBOV) genome substitution rate in the Makona strain has been estimated at between 0.87 × 10−3 and 1.42 × 10−3 mutations per site per year. This is equivalent to 16–27 mutations in each genome, meaning that sequences diverge rapidly enough to identify distinct sub-lineages during a prolonged epidemic2,3,4,5,6,7. Genome sequencing provides a high-resolution view of pathogen evolution and is increasingly sought after for outbreak surveillance. Sequence data may be used to guide control measures, but only if the results are generated quickly enough to inform interventions8. Genomic surveillance during the epidemic has been sporadic owing to a lack of local sequencing capacity coupled with practical difficulties transporting samples to remote sequencing facilities9. To address this problem, here we devise a genomic surveillance system that utilizes a novel nanopore DNA sequencing instrument. In April 2015 this system was transported in standard airline luggage to Guinea and used for real-time genomic surveillance of the ongoing epidemic. We present sequence data and analysis of 142 EBOV samples collected during the period March to October 2015. We were able to generate results less than 24 h after receiving an Ebola-positive sample, with the sequencing process taking as little as 15–60 min. We show that real-time genomic surveillance is possible in resource-limited settings and can be established rapidly to monitor outbreaks.

This is a preview of subscription content, access via your institution

Access options

Rent or buy this article

Prices vary by article type

from$1.95

to$39.95

Prices may be subject to local taxes which are calculated during checkout

Figure 1: Deployment of the portable genome surveillance system in Guinea.
Figure 2: Real-time genomics surveillance in context of the Guinea Ebola virus disease epidemic.
Figure 3: Evolution of EBOV over the course of the Ebola virus disease epidemic.

Similar content being viewed by others

Accession codes

Primary accessions

European Nucleotide Archive

References

  1. World Health Organisation. Ebola Situation Report - 11 November 2015. http://apps.who.int/ebola/current-situation/ebola-situation-report-11-november-2015 (2015)

  2. Gire, S. K. et al. Genomic surveillance elucidates Ebola virus origin and transmission during the 2014 outbreak. Science 345, 1369–1372 (2014)

    Article  CAS  ADS  Google Scholar 

  3. Carroll, M. W. et al. Temporal and spatial analysis of the 2014–2015 Ebola virus outbreak in West Africa. Nature 524, 97–101 (2015)

    Article  CAS  ADS  Google Scholar 

  4. Simon-Loriere, E. et al. Distinct lineages of Ebola virus in Guinea during the 2014 West African epidemic. Nature 524, 102–104 (2015)

    Article  CAS  ADS  Google Scholar 

  5. Park, D. J. et al. Ebola virus epidemiology, transmission, and evolution during seven months in Sierra Leone. Cell 161, 1516–1526 (2015)

    Article  CAS  Google Scholar 

  6. Tong, Y.-G. et al. Genetic diversity and evolutionary dynamics of Ebola virus in Sierra Leone. Nature 524, 93–96 (2015)

    Article  CAS  Google Scholar 

  7. Kugelman, J. R. et al. Monitoring of Ebola Virus Makona evolution through establishment of advanced genomic capability in Liberia. Emerg. Infect. Dis. 21, 1135–1143 (2015)

    Article  CAS  Google Scholar 

  8. Gardy, J., Loman, N. J. & Rambaut, A. Real-time digital pathogen surveillance—the time is now. Genome Biol. 16, 155 (2015)

    Article  Google Scholar 

  9. Yozwiak, N. L., Schaffner, S. F. & Sabeti, P. C. Data sharing: make outbreak research open access. Nature 518, 477–479 (2015)

    Article  CAS  ADS  Google Scholar 

  10. Heger, M. Liberia’s LIBR Genome Center monitors Ebola outbreak, emerging pathogens. https://www.genomeweb.com/sequencing-technology/liberias-libr-genome-center-monitors-ebola-outbreak-emerging-pathogens (2015)

  11. Quick, J., Quinlan, A. R. & Loman, N. J. A reference bacterial genome dataset generated on the MinION™ portable single-molecule nanopore sequencer. Gigascience 3, 22 (2014)

    Article  Google Scholar 

  12. Jain, M. et al. Improved data analysis for the MinION nanopore sequencer. Nature Methods 12, 351–356 (2015)

    Article  CAS  Google Scholar 

  13. Urban, J. M., Bliss, J., Lawrence, C. E. & Gerbi, S. A. Sequencing ultra-long DNA molecules with the Oxford Nanopore MinION. Preprint at http://biorxiv.org/content/early/2015/05/13/019281 (2015)

  14. Quick, J. et al. Rapid draft sequencing and real-time nanopore sequencing in a hospital outbreak of Salmonella. Genome Biol. 16, 114 (2015)

    Article  Google Scholar 

  15. Greninger, A. L. et al. Rapid metagenomic identification of viral pathogens in clinical samples by real-time nanopore sequencing analysis. Genome Med . 7, 99 (2015)

    Article  Google Scholar 

  16. Loman, N. J., Quick, J. & Simpson, J. T. A complete bacterial genome assembled de novo using only nanopore sequencing data. Nature Methods 12, 733–735 (2015)

    Article  CAS  Google Scholar 

  17. Goodfellow, I. et al. Recent evolution patterns of Ebola virus obtained by direct sequencing in Sierra Leone. http://virological.org/t/recent-evolution-patterns-of-ebola-virus-obtained-by-direct-sequencing-in-sierra-leone/150 (2015)

  18. Herold, K. E., Sergeev, N., Matviyenko, A. & Rasooly, A. in Biosensors and Biodetection 504, 441–458 (Humana, 2009)

    Article  CAS  Google Scholar 

  19. Mate, S. E. et al. Molecular evidence of sexual transmission of Ebola virus. N. Engl. J. Med. 373, 2448–2454 (2015)

    Article  CAS  Google Scholar 

  20. Smither, S. J. et al. Buffer AVL alone does not inactivate Ebola virus in a representative clinical sample type. J. Clin. Microbiol. 53, 3148–3154 (2015)

    Article  Google Scholar 

  21. Loman, N. J. & Quinlan, A. R. Poretools: a toolkit for analyzing nanopore sequence data. Bioinformatics 30, 3399–3401 (2014)

    Article  CAS  Google Scholar 

  22. Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. Preprint at http://arxiv.org/abs/1303.3997 (2013)

  23. Garrison, E. & Marth, G. Haplotype-based variant detection from short-read sequencing. Preprint at http://arxiv.org/abs/1207.3907 (2012)

  24. Kurtz, S. et al. Versatile and open software for comparing large genomes. Genome Biol. 5, R12 (2004)

    Article  Google Scholar 

  25. Stamatakis, A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312–1313 (2014)

    Article  CAS  Google Scholar 

  26. Edgar, R. C. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797 (2004)

    Article  CAS  Google Scholar 

  27. Price, M. N., Dehal, P. S. & Arkin, A. P. FastTree 2–approximately maximum-likelihood trees for large alignments. PLoS One 5, e9490 (2010)

    Article  ADS  Google Scholar 

  28. Drummond, A. J. & Rambaut, A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol. Biol. 7, 214 (2007)

    Article  Google Scholar 

  29. Hasegawa, M., Kishino, H. & Yano, T. Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J. Mol. Evol. 22, 160–174 (1985)

    Article  CAS  ADS  Google Scholar 

  30. Yang, Z. Maximum likelihood phylogenetic estimation from DNA sequences with variable rates over sites: approximate methods. J. Mol. Evol. 39, 306–314 (1994)

    Article  CAS  ADS  Google Scholar 

  31. Gill, M. S. et al. Improving Bayesian population dynamics inference: a coalescent-based model for multiple loci. Mol. Biol. Evol. 30, 713–724 (2013)

    Article  CAS  Google Scholar 

  32. Drummond, A. J., Ho, S. Y. W., Phillips, M. J. & Rambaut, A. Relaxed phylogenetics and dating with confidence. PLoS Biol. 4, e88 (2006)

    Article  Google Scholar 

  33. Cock, P. J. et al. Biopython: freely available Python tools for computational molecular biology and bioinformatics. Bioinformatics 25, 1422–1423 (2009)

    Article  CAS  Google Scholar 

Download references

Acknowledgements

The EMLab is a technical partner in the WHO Emerging and Dangerous Pathogens Laboratory Network (EDPLN), and the Global Outbreak Alert and Response Network (GOARN) and the deployments in West Africa have been coordinated and supported by the GOARN Operational Support Team at WHO/HQ and the African Union. This work was carried out in the context of the project EVIDENT (Ebola virus disease: correlates of protection, determinants of outcome, and clinical management) that received funding from the European Union’s Horizon 2020 research and innovation program under grant agreement No 666100 and in the context of service contract IFS/2011/272-372 funded by Directorate-General for International Cooperation and Development. J.Q. is funded by the NIHR Surgical Reconstruction and Microbiology Research Centre (SRMRC). N.J.L. is funded by a Medical Research Council Special Training Fellowship in Biomedical Informatics (to September 2015) and a Medical Research Council Bioinformatics Fellowship. J.T.S. is supported by the Ontario Institute for Cancer Research through funding provided by the Government of Ontario. Dstl support was funded by the UK Ministry of Defence (MOD). Dstl authors thank S. Lonsdale, C. Lonsdale and C. Mayers for supply of RNA, previous assistance, and review of the manuscript. The views expressed in this paper are not necessarily endorsed by the UK MOD. A.R. was supported by EU Seventh Framework Programme [FP7/2007-2013] under Grant Agreement no. 278433-PREDEMICS and ERC Grant agreement no. 260864. We are grateful for the generous support of University of Birmingham alumni for donations in support of the pilot work. The MRC Cloud Infrastructure for Microbial Bioinformatics (CLIMB) cyberinfrastructure was used to conduct bioinformatics analysis. The authors would like to thank B. Oppenheim and C. Wardius for help with logistics and the staff of Alta Biosciences, University of Birmingham and Sigma-Aldrich for generating PCR primers especially rapidly for this project. The authors would like to thank scientists deployed from the Special Pathogens Program from the National Microbiology Laboratory, Public Health Agency of Canada, who worked on EBOV diagnostics in Guinea. We are grateful to I. Goodfellow, M. Cotten and P. Kellam for permission to include sequences from Sierra Leone in this analysis. We thank R. Vipond for assistance with validation experiments. We thank H. Eno and B. Myers for help with proofreading. We are thankful for the generous support of reagents and technical support from Oxford Nanopore. We thank the staff at Oxford Nanopore for technical and logistical support during this project with special thanks to S. Brooking, O. Hartwell, R. Pettett, C. Brown, G. Sanghera and R. Ronan. We thank T. Bedford and R. Neher for developing the Nextstrain website.

Author information

Authors and Affiliations

Authors

Contributions

N.J.L., J.Q., M.K.O’S., D.W., S.G., M.W.C. conceived the study. N.J.L., J.Q., M.K.O’S., S.A.W., J.T., P.R., D.T. designed the lab in a suitcase and laboratory protocol and initial validation. J.Q., S.D., L.C., J.A.B., R.K., L.E.K., and A.Ma. performed MinION sequencing. N.J.L., J.Q. and J.T.S. performed bioinformatics analysis and wrote software. J.T.S. added variant calling support to the nanopolish software. N.J.L., J.Q., S.D., E.S., P.F., L.C., A.Mi., N.M. and I.R. analysed the data. G.D., A.R., N.J.L., J.Q. and G.P. performed phylogenetic analysis. J.A.H., D.A.M., G.P., K.L., B.A. assisted further validation experiments. M.W.C., M.Ga., S.G., A.D.C., K.S., E.F. and R.W. coordinated activities for the European Mobile Laboratories. N.J.L., J.Q., S.D., M.W.C., S.G., M.K.O’S., A.R., E.S., P.F., I.R., A.Mi., and L.C. wrote the manuscript. All other authors were involved either in sample collection, and/or logistical support and strategic oversight for the work.

Corresponding author

Correspondence to Nicholas J. Loman.

Ethics declarations

Competing interests

J.Q., N.J.L. and J.T.S. have all received travel expenses and accommodation from Oxford Nanopore to speak at organised symposia. J.Q. and N.J.L. have received an honorarium payment to speak at an Oxford Nanopore meeting. N.J.L. is a member of the Oxford Nanopore MinION Access Programme and has received reagents free of charge as part of the MinION Access Programme and in support of this project but does not receive other financial compensation or hold shares. D.T. is an employee of Oxford Nanopore.

Additional information

MinION and Illumina raw sequence files have been deposited into the European Nucleotide Archive under project code PRJEB10571.

Extended data figures and tables

Extended Data Figure 1 Primer schemes employed during the study.

We designed PCR primers to generate amplicons that would span the EBOV genome. a, We initially designed 38 primer pairs which were used in the initial validation study and which cover >97% of the EBOV genome. During in-field sequencing we used a 19-reaction scheme or 11-reaction scheme, which generated longer products. The predicted amplicon products are shown with forward primers and reverse primers indicated by green bars on the forward and reverse strand, respectively, scaled according to the EBOV virus coordinates. b, c, The amplicon product sizes expected are shown for the 19-reaction scheme (b) and the 11-reaction scheme (c). No amplicon covers the extreme 3′ region of the genome. The last primer pair, 38_R, ends at position 18578, 381 bases away from the end of the virus genome. The primer diagram was created with Biopython33.

Extended Data Figure 2 List of equipment and consumables to establish the genome surveillance system.

ac, We show the list of equipment (a), disposable consumables (b) and reagents (c) to establish in-field genomic surveillance. Sufficient reagents were shipped for 20 samples. MinION sequencing requires a mix of chilled and frozen reagents. Recommended shipping conditions are specified. The picture underneath depicts MinION flowcells ready for shipping with insulating material (left) and frozen reagents (right).

Extended Data Figure 3 Bioinformatics workflow.

This figure summarizes the steps performed during bioinformatics analysis (ordered from top to bottom), in order to generate consensus sequences. The right column shows the example software command executed at each step.

Extended Data Figure 4 Results of MinION validation.

a, The results of comparing four MinION sequences with Illumina sequences generated as part of a previous study3 are shown. Each row in the table demonstrates the number of true positives, false positives and false negatives for a sample. False negatives may result in masked sequences, owing to being outside of regions covered by the amplicon scheme, having low coverage or falling within a primer binding site. Results before and after quality filtering (log likelihood ratio of >200) are shown. After quality filtering, no false positive calls were detected. All detected false negatives were masked with Ns in the final consensus sequence. No positions were called incorrectly. b, The four consensus sequences, plus an additional sample that had missing coverage in one amplicon are shown as part of a phylogenetic reconstruction with genomes from Carroll et al.3. Sample labels in red, blue, pink, yellow and blue represent pairs of sequences generated on MinION and llumina. These fall into identical clusters.

Extended Data Figure 5 Relationship between coverqage and log-likelihood ratio for sample 076769.

Line-plot showing the relationship between sequence depth of coverage (x axis) and the log likelihood ratio for detected SNPs derived by subsampling reads from a single sequencing run to simulate the effect of low coverage. The horizontal and vertical line indicates the cut-offs (quality and coverage respectively) for consensus calling. Therefore, all variants are detected below 25× coverage, and the vast majority meet the threshold quality at 25× coverage or slightly above. Any combination of log likelihood ratio or coverage that placed variants in the grey box would be represented as a masked position in the final consensus sequence.

Extended Data Figure 6 Duration of MinION sequencing runs.

For each sequence run the sequencing duration, measured as the difference between timestamp of the first read seen and the last read transferred for analysis. 127 runs are shown, with 15 outliers with duration greater than 200 min excluded.

Extended Data Figure 7 Histogram of Ct values for study samples.

Ct values for samples in the study (where information was available) ranged between 13.8 and 35.7, with a mean of 22.

Extended Data Figure 8 Sequence accuracy for samples.

a, b, Accuracy measurements for the entire set of two-direction reads were made for the validation samples, sequenced in the United Kingdom (a) and each of the 142 samples from real-time genomic surveillance (b). Accuracy is defined according to the definition from Quick et al.11. Vertical dashed lines indicate the mean accuracy for the sample.

Extended Data Figure 9 Maximum likelihood phylogenetic inference of 125 Ebola virus samples from this study with 603 previously published sequences.

Coloured nodes are from this study. Node shape reflects country of origin. ac, the entire data set is shown (a), with zoomed regions focusing on lineages GN1 (b) and SL3 (c) identified during real-time sequencing. Map figure adapted from SimpleMaps website (http://simplemaps.com/resources/svg-gn).

Extended Data Figure 10 Root-to-tip divergence plot and mean evolutionary rate estimate.

a, Root-to-tip divergence plot for the 728 Ebola samples generated through maximum likelihood analysis. Samples from real-time genomic surveillance are coloured as per Fig. 3 and Extended Data Fig. 9. b, Mean evolutionary rate estimate (in substitutions per site per year) across the EBOV phylogeny recovered using BEAST under a relaxed lognormal molecular clock. Blue area corresponds to the 95% highest posterior density (HPD) (mean of the distribution is 1.19 × 10−3, 95% HPDs: 1.09–1.29 × 10−3 substitutions per site per year). Hatched regions in red are outside the 95% HPD intervals.

Supplementary information

Supplementary Information

This file contains a Field Guide to Nanopore Sequencing - a detailed discussion of logistical issues that arose during this project and Supplementary Tables 1-4. (PDF 428 kb)

PowerPoint slides

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Quick, J., Loman, N., Duraffour, S. et al. Real-time, portable genome sequencing for Ebola surveillance. Nature 530, 228–232 (2016). https://doi.org/10.1038/nature16996

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/nature16996

This article is cited by

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing