ORF-Interrupting Mutations in Monkeypox Virus Genomes from Washington and Ohio, 2022

Monkeypox virus, the causative agent of the 2022 monkeypox outbreak, is a double-stranded DNA virus in the Orthopoxvirus genus of the Poxviridae family. Genes in terminal regions of Orthopoxvirus genomes mostly code for host-pathogen interaction proteins and are prone to selective pressure and modification events. Using viral whole genome sequencing, we identified twenty-five total clinical samples with ORF-disrupting mutations, including twenty samples encoding nonsense mutations in MPXVgp001/191 (OPG001), MPXVgp004/188 (OPG015), MPXVgp010 (OPG023), MPXVgp030 (OPG042), MPXVgp159 (OPG0178), or MPXVgp161 (OPG181). Additional mutations include a frameshift leading to an alternative C-terminus in MPXVgp010 (OPG023) and an insertion in an adenine homopolymer at the beginning of the annotated ORF for MPXVgp153 (OPG151), encoding a subunit of the RNA polymerase, suggesting the virus may instead use the start codon that encodes Met9 as annotated. Finally, we detected three samples with large (>900 bp) deletions. These included a 913 bp deletion that truncates the C-terminus of MPXVgp010 (OPG023); a 4205 bp deletion that eliminates MPXVgp012 (OPG025), MPXVgp013 (OPG027), and MPXVgp014 (OPG029) and truncates MPXVgp011 (OPG024; D8L) and MPXVgp015 (OPG030); and a 6881 bp deletion that truncates MPXVgp182 (OPG210) and eliminates putative ORFs MPXVgp184, MPXVgp185 (OPG005), and MPXVgp186, as well as MPXVgp187 (OPG016), and MPXVgp188 (OPG015) from the 3’ ITR only. MPXVgp182 encodes the monkeypox-specific, highly immunogenic surface glycoprotein B21R which has been proposed as a serological target. Overall, we find greater than one-tenth of our sequenced MPXV isolates have at least one gene inactivating mutation and these genes together comprised greater than one-tenth of annotated MPXV genes. Our findings highlight non-essential genes in monkeypox virus that may be evolving as a result of selective pressure in humans, as well as the limitations of targeting them for therapeutics and diagnostic testing.

[1]  A. Greninger,et al.  Genomic analysis of early spread of monkeypox virus in Washington State , 2022, medRxiv.

[2]  A. Greninger,et al.  Clinical Performance and Trends during the First Two Months of Monkeypox Virus PCR Testing at Two United States Reference Labs , 2022, medRxiv.

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

[4]  G. Cheng,et al.  Genomic annotation and molecular evolution of monkeypox virus outbreak in 2022 , 2022, Journal of medical virology.

[5]  V. Corman,et al.  Genetic variability, including gene duplication and deletion, in early sequences from the 2022 European monkeypox outbreak , 2022, bioRxiv.

[6]  P. Desingu,et al.  Genomic regions insertion and deletion in Monkeypox virus causing multi-country outbreak-2022 , 2022, bioRxiv.

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

[8]  David H. Philpott,et al.  Monkeypox Outbreak — Nine States, May 2022 , 2022, MMWR. Morbidity and mortality weekly report.

[9]  Konstantinos D. Tsirigos,et al.  SignalP 6.0 predicts all five types of signal peptides using protein language models , 2022, Nature Biotechnology.

[10]  Silvio C. E. Tosatto,et al.  Pfam: The protein families database in 2021 , 2020, Nucleic Acids Res..

[11]  U. Moens,et al.  Monkeypox Virus in Nigeria: Infection Biology, Epidemiology, and Evolution , 2020, Viruses.

[12]  Anna M. Mandra,et al.  Exportation of Monkeypox Virus From the African Continent , 2020, The Journal of infectious diseases.

[13]  Narmada Thanki,et al.  CDD/SPARCLE: the conserved domain database in 2020 , 2019, Nucleic Acids Res..

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

[15]  Trevor Bedford,et al.  Nextstrain: real-time tracking of pathogen evolution , 2017, bioRxiv.

[16]  S. Shchelkunov,et al.  Species-specific differentiation of variola, monkeypox, and varicella-zoster viruses by multiplex real-time PCR assay , 2016, Journal of virological methods.

[17]  B. Moss,et al.  Deciphering Poxvirus Gene Expression by RNA Sequencing and Ribosome Profiling , 2015, Journal of Virology.

[18]  K. Früh,et al.  T Cell Inactivation by Poxviral B22 Family Proteins Increases Viral Virulence , 2014, PLoS pathogens.

[19]  Björn Usadel,et al.  Trimmomatic: a flexible trimmer for Illumina sequence data , 2014, Bioinform..

[20]  Brian Bushnell,et al.  BBMap: A Fast, Accurate, Splice-Aware Aligner , 2014 .

[21]  Jeffrey R. Kugelman,et al.  Genomic Variability of Monkeypox Virus among Humans, Democratic Republic of the Congo , 2014, Emerging infectious diseases.

[22]  E. Kroon,et al.  Characterization of a New Vaccinia virus Isolate Reveals the C23L Gene as a Putative Genetic Marker for Autochthonous Group 1 Brazilian Vaccinia virus , 2012, PloS one.

[23]  Jay Shendure,et al.  Poxviruses Deploy Genomic Accordions to Adapt Rapidly against Host Antiviral Defenses , 2012, Cell.

[24]  S. Shchelkunov.  Orthopoxvirus Genes That Mediate Disease Virulence and Host Tropism , 2012, Advances in virology.

[25]  A. Legasse,et al.  Deletion of the Monkeypox Virus Inhibitor of Complement Enzymes Locus Impacts the Adaptive Immune Response to Monkeypox Virus in a Nonhuman Primate Model of Infection , 2011, Journal of Virology.

[26]  D. Esteban,et al.  Genes in the terminal regions of orthopoxvirus genomes experience adaptive molecular evolution , 2011, BMC Genomics.

[27]  Eneida L. Hatcher,et al.  Orthopoxvirus Genome Evolution: The Role of Gene Loss , 2010, Viruses.

[28]  Marc A Suchard,et al.  Using Time-Structured Data to Estimate Evolutionary Rates of Double-Stranded DNA Viruses , 2010, Molecular biology and evolution.

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

[30]  S. Pochet,et al.  Crystal structure of poxvirus thymidylate kinase: An unexpected dimerization has implications for antiviral therapy , 2008, Proceedings of the National Academy of Sciences.

[31]  R. Owens,et al.  Vaccinia Virus Proteins A52 and B14 Share a Bcl-2–Like Fold but Have Evolved to Inhibit NF-κB rather than Apoptosis , 2008, PLoS pathogens.

[32]  David S. Trossman,et al.  Zoonotic orthopoxviruses encode a high-affinity antagonist of NKG2D , 2007, The Journal of experimental medicine.

[33]  Jon M Hanifin,et al.  Multiple diagnostic techniques identify previously vaccinated individuals with protective immunity against monkeypox , 2005, Nature Medicine.

[34]  P. Jahrling,et al.  Virulence differences between monkeypox virus isolates from West Africa and the Congo basin , 2005, Virology.