Development of an amplicon-based sequencing approach in response to the global emergence of mpox

The 2022 multicountry mpox outbreak concurrent with the ongoing Coronavirus Disease 2019 (COVID-19) pandemic further highlighted the need for genomic surveillance and rapid pathogen whole-genome sequencing. While metagenomic sequencing approaches have been used to sequence many of the early mpox infections, these methods are resource intensive and require samples with high viral DNA concentrations. Given the atypical clinical presentation of cases associated with the outbreak and uncertainty regarding viral load across both the course of infection and anatomical body sites, there was an urgent need for a more sensitive and broadly applicable sequencing approach. Highly multiplexed amplicon-based sequencing (PrimalSeq) was initially developed for sequencing of Zika virus, and later adapted as the main sequencing approach for Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). Here, we used PrimalScheme to develop a primer scheme for human monkeypox virus that can be used with many sequencing and bioinformatics pipelines implemented in public health laboratories during the COVID-19 pandemic. We sequenced clinical specimens that tested presumptively positive for human monkeypox virus with amplicon-based and metagenomic sequencing approaches. We found notably higher genome coverage across the virus genome, with minimal amplicon drop-outs, in using the amplicon-based sequencing approach, particularly in higher PCR cycle threshold (Ct) (lower DNA titer) samples. Further testing demonstrated that Ct value correlated with the number of sequencing reads and influenced the percent genome coverage. To maximize genome coverage when resources are limited, we recommend selecting samples with a PCR Ct below 31 Ct and generating 1 million sequencing reads per sample. To support national and international public health genomic surveillance efforts, we sent out primer pool aliquots to 10 laboratories across the United States, United Kingdom, Brazil, and Portugal. These public health laboratories successfully implemented the human monkeypox virus primer scheme in various amplicon sequencing workflows and with different sample types across a range of Ct values. Thus, we show that amplicon-based sequencing can provide a rapidly deployable, cost-effective, and flexible approach to pathogen whole-genome sequencing in response to newly emerging pathogens. Importantly, through the implementation of our primer scheme into existing SARS-CoV-2 workflows and across a range of sample types and sequencing platforms, we further demonstrate the potential of this approach for rapid outbreak response.

[1]  N. K. Shah,et al.  An Mpox-Related Death in the United States. , 2023, The New England journal of medicine.

[2]  C. Carrington,et al.  Toward a global virus genomic surveillance network , 2023, Cell Host & Microbe.

[3]  P. Hemarajata,et al.  Identification of Human Monkeypox Virus Genome Deletions That Impact Diagnostic Assays , 2022, Journal of clinical microbiology.

[4]  G. Wallau,et al.  Possible Occupational Infection of Healthcare Workers with Monkeypox Virus, Brazil , 2022, Emerging infectious diseases.

[5]  H Zhao,et al.  Genomic deletions and rearrangements in monkeypox virus from the 2022 outbreak, USA , 2022, bioRxiv.

[6]  M. Kraemer,et al.  Urgent need for a non-discriminatory and non-stigmatizing nomenclature for monkeypox virus , 2022, PLoS biology.

[7]  A. Antinori,et al.  Monkeypox Virus Infection in Humans across 16 Countries - April-June 2022. , 2022, The New England journal of medicine.

[8]  N. Green,et al.  Rapid Diagnostic Testing for Response to the Monkeypox Outbreak — Laboratory Response Network, United States, May 17–June 30, 2022 , 2022, MMWR. Morbidity and mortality weekly report.

[9]  J. Blanco,et al.  Frequent detection of monkeypox virus DNA in saliva, semen, and other clinical samples from 12 patients, Barcelona, Spain, May to June 2022 , 2022, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[10]  D. Sobral,et al.  Phylogenomic characterization and signs of microevolution in the 2022 multi-country outbreak of monkeypox virus , 2022, Nature Medicine.

[11]  Todd G. Smith,et al.  Multiple lineages of Monkeypox virus detected in the United States, 2021- 2022 , 2022, bioRxiv.

[12]  J. Walley,et al.  Global human monkeypox outbreak: atypical presentation demanding urgent public health action , 2022, The Lancet. Microbe.

[13]  Jennifer L. Small,et al.  Clinical features and management of human monkeypox: a retrospective observational study in the UK , 2022, The Lancet. Infectious diseases.

[14]  A. Christoffels,et al.  The evolving SARS-CoV-2 epidemic in Africa: Insights from rapidly expanding genomic surveillance , 2022, Science.

[15]  A. Azman,et al.  Global landscape of SARS-CoV-2 genomic surveillance and data sharing , 2022, Nature Genetics.

[16]  P. Fariselli,et al.  SeqFu: A Suite of Utilities for the Robust and Reproducible Manipulation of Sequence Files , 2021, Bioengineering.

[17]  Andrew D Smith,et al.  Improvements to the ARTIC multiplex PCR method for SARS-CoV-2 genome sequencing using nanopore , 2020, bioRxiv.

[18]  Jaime Fern'andez del R'io,et al.  Array programming with NumPy , 2020, Nature.

[19]  C. Chiu,et al.  Clinical metagenomics , 2019, Nature Reviews Genetics.

[20]  P. Lemey,et al.  Tracking virus outbreaks in the twenty-first century , 2018, Nature Microbiology.

[21]  Karthik Gangavarapu,et al.  An amplicon-based sequencing framework for accurately measuring intrahost virus diversity using PrimalSeq and iVar , 2018, Genome Biology.

[22]  A. Lauring,et al.  Complexities of Viral Mutation Rates , 2018, Journal of Virology.

[23]  N. Segata,et al.  Shotgun metagenomics, from sampling to analysis , 2017, Nature Biotechnology.

[24]  Trevor Bedford,et al.  Multiplex PCR method for MinION and Illumina sequencing of Zika and other virus genomes directly from clinical samples , 2017, Nature Protocols.

[25]  C. Wilke Streamlined Plot Theme and Plot Annotations for 'ggplot2' , 2015 .

[26]  Andrew Rambaut,et al.  Real-time digital pathogen surveillance — the time is now , 2015, Genome Biology.

[27]  Yihui Xie,et al.  Create Interactive Web Maps with the JavaScript 'Leaflet'Library , 2015 .

[28]  Rachel S. G. Sealfon,et al.  Genomic surveillance elucidates Ebola virus origin and transmission during the 2014 outbreak , 2014, Science.

[29]  S. Shchelkunov,et al.  Species-specific identification of variola, monkeypox, cowpox, and vaccinia viruses by multiplex real-time PCR assay , 2011, Journal of virological methods.

[30]  I. Damon,et al.  Real-time PCR assays for the specific detection of monkeypox virus West African and Congo Basin strain DNA , 2010, Journal of virological methods.

[31]  Aaron R. Quinlan,et al.  BEDTools: a flexible suite of utilities for comparing genomic features , 2010, Bioinform..

[32]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[33]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[34]  Bartek Wilczynski,et al.  Biopython: freely available Python tools for computational molecular biology and bioinformatics , 2009, Bioinform..

[35]  G. Pauli,et al.  Orthopoxvirus Detection in Environmental Specimens during Suspected Bioterror Attacks: Inhibitory Influences of Common Household Products , 2007, Applied and Environmental Microbiology.

[36]  I. Damon,et al.  Detection of monkeypox virus with real-time PCR assays , 2006, Journal of clinical virology : the official publication of the Pan American Society for Clinical Virology.

[37]  R. Jordan,et al.  An Orally Bioavailable Antipoxvirus Compound (ST-246) Inhibits Extracellular Virus Formation and Protects Mice from Lethal Orthopoxvirus Challenge , 2005, Journal of Virology.

[38]  H. Ellerbrok,et al.  Detection of Orthopoxvirus DNA by Real-Time PCR and Identification of Variola Virus DNA by Melting Analysis , 2004, Journal of Clinical Microbiology.

[39]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .