Identification of Nucleotide-Level Changes Impacting Gene Content and Genome Evolution in Orthopoxviruses

ABSTRACT Poxviruses are composed of large double-stranded DNA (dsDNA) genomes coding for several hundred genes whose variation has supported virus adaptation to a wide variety of hosts over their long evolutionary history. Comparative genomics has suggested that the Orthopoxvirus genus in particular has undergone reductive evolution, with the most recent common ancestor likely possessing a gene complement consisting of all genes present in any existing modern-day orthopoxvirus species, similar to the current Cowpox virus species. As orthopoxviruses adapt to new environments, the selection pressure on individual genes may be altered, driving sequence divergence and possible loss of function. This is evidenced by accumulation of mutations and loss of protein-coding open reading frames (ORFs) that progress from individual missense mutations to gene truncation through the introduction of early stop mutations (ESMs), gene fragmentation, and in some cases, a total loss of the ORF. In this study, we have constructed a whole-genome alignment for representative isolates from each Orthopoxvirus species and used it to identify the nucleotide-level changes that have led to gene content variation. By identifying the changes that have led to ESMs, we were able to determine that short indels were the major cause of gene truncations and that the genome length is inversely proportional to the number of ESMs present. We also identified the number and types of protein functional motifs still present in truncated genes to assess their functional significance. IMPORTANCE This work contributes to our understanding of reductive evolution in poxviruses by identifying genomic remnants such as single nucleotide polymorphisms (SNPs) and indels left behind by evolutionary processes. Our comprehensive analysis of the genomic changes leading to gene truncation and fragmentation was able to detect some of the remnants of these evolutionary processes still present in orthopoxvirus genomes and suggests that these viruses are under continual adaptation due to changes in their environment. These results further our understanding of the evolutionary mechanisms that drive virus variation, allowing orthopoxviruses to adapt to particular environmental niches. Understanding the evolutionary history of these virus pathogens may help predict their future evolutionary potential.

[1]  F. Fenner,et al.  The biological characters of several strains of vaccinia, cowpox and rabbitpox viruses. , 1958, Virology.

[2]  Rodrigo Lopez,et al.  Clustal W and Clustal X version 2.0 , 2007, Bioinform..

[3]  C. Upton,et al.  Characterization of indels in poxvirus genomes , 2011, Virus Genes.

[4]  C Upton,et al.  Poxviruses: past, present and future. , 2006, Virus research.

[5]  T. Traavik,et al.  Recombinant viruses obtained from co-infection in vitro with a live vaccinia-vectored influenza vaccine and a naturally occurring cowpox virus display different plaque phenotypes and loss of the transgene. , 2004, Vaccine.

[6]  Nicolas Lartillot,et al.  PhyloBayes 3: a Bayesian software package for phylogenetic reconstruction and molecular dating , 2009, Bioinform..

[7]  H. Meyer,et al.  Zoonotic poxviruses , 2009, Veterinary microbiology.

[8]  K. Katoh,et al.  MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. , 2002, Nucleic acids research.

[9]  D. Posada jModelTest: phylogenetic model averaging. , 2008, Molecular biology and evolution.

[10]  F. Fenner,et al.  Chapter 3 – Mousepox , 2007 .

[11]  John Healy,et al.  GapCoder automates the use of indel characters in phylogenetic analysis , 2003, BMC Bioinformatics.

[12]  Mark A. Miller,et al.  Creating the CIPRES Science Gateway for inference of large phylogenetic trees , 2010, 2010 Gateway Computing Environments Workshop (GCE).

[13]  Yun Zhang,et al.  ViPR: an open bioinformatics database and analysis resource for virology research , 2011, Nucleic Acids Res..

[14]  G. McFadden Poxvirus tropism , 2005, Nature Reviews Microbiology.

[15]  M. Rosenberg,et al.  How should gaps be treated in parsimony? A comparison of approaches using simulation. , 2007, Molecular phylogenetics and evolution.

[16]  N. Moran,et al.  Deletional bias and the evolution of bacterial genomes. , 2001, Trends in genetics : TIG.

[17]  G. McFadden,et al.  Technical knockout: understanding poxvirus pathogenesis by selectively deleting viral immunomodulatory genes , 2004, Cellular microbiology.

[18]  A. Mercer,et al.  Poxvirus Host Range Protein CP77 Contains an F-Box-Like Domain That Is Necessary To Suppress NF-κB Activation by Tumor Necrosis Factor Alpha but Is Independent of Its Host Range Function , 2009, Journal of Virology.

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

[20]  M. Frace,et al.  The Phylogenetics and Ecology of the Orthopoxviruses Endemic to North America , 2009, PloS one.

[21]  Kirsten A. Bratke,et al.  A survey of host range genes in poxvirus genomes. , 2013, Infection, genetics and evolution : journal of molecular epidemiology and evolutionary genetics in infectious diseases.

[22]  David H. L. Bishop,et al.  The International Committee on Taxonomy of Viruses , 1995 .

[23]  S. Shchelkunov,et al.  The genomic sequence analysis of the left and right species-specific terminal region of a cowpox virus strain reveals unique sequences and a cluster of intact ORFs for immunomodulatory and host range proteins. , 1998, Virology.

[24]  P. Kellam,et al.  Poxvirus genomes: a phylogenetic analysis. , 2004, The Journal of general virology.

[25]  F. Fenner,et al.  [Genetic studies with mammalian poxviruses. I. Demonstration of recombination between two strains of vaccina virus]. , 1960, Virology.

[26]  F. Fenner,et al.  Genetic studies with mammalian poxviruses. IV. Hybridization between several different poxviruses. , 1960, Virology.

[27]  Kirsten A. Bratke,et al.  Identification of multiple independent horizontal gene transfers into poxviruses using a comparative genomics approach , 2008, BMC Evolutionary Biology.

[28]  S. K. Rasmussen,et al.  Capturing the Natural Diversity of the Human Antibody Response against Vaccinia Virus , 2010, Journal of Virology.

[29]  E. Myers,et al.  Basic local alignment search tool. , 1990, Journal of molecular biology.

[30]  G. McFadden,et al.  Immunopathogenesis of poxvirus infections: forecasting the impending storm , 2007, Immunology and cell biology.

[31]  Pravech Ajawatanawong,et al.  Evolution of protein indels in plants, animals and fungi , 2013, BMC Evolutionary Biology.

[32]  Ramón Doallo,et al.  CircadiOmics: integrating circadian genomics, transcriptomics, proteomics and metabolomics , 2012, Nature Methods.

[33]  G. McFadden,et al.  When good vaccines go wild: Feral Orthopoxvirus in developing countries and beyond. , 2008, Journal of infection in developing countries.

[34]  Stephanie Irausquin,et al.  The evolutionary biology of poxviruses. , 2010, Infection, genetics and evolution : journal of molecular epidemiology and evolutionary genetics in infectious diseases.

[35]  C. Upton,et al.  Genomic Analysis of the Vaccinia Virus Strain Variants Found in Dryvax Vaccine , 2011, Journal of Virology.

[36]  N. Okada,et al.  Poxviruses as possible vectors for horizontal transfer of retroposons from reptiles to mammals , 2007, Proceedings of the National Academy of Sciences.

[37]  M. Reynolds,et al.  Factors affecting the likelihood of monkeypox's emergence and spread in the post-smallpox era , 2012, Current opinion in virology.

[38]  Yu Li,et al.  Genome Sequence Diversity and Clues to the Evolution of Variola (Smallpox) Virus , 2006, Science.

[39]  B. Moss,et al.  Simultaneous high-resolution analysis of vaccinia virus and host cell transcriptomes by deep RNA sequencing , 2010, Proceedings of the National Academy of Sciences.

[40]  John P. Huelsenbeck,et al.  MRBAYES: Bayesian inference of phylogenetic trees , 2001, Bioinform..

[41]  S. Shchelkunov,et al.  Orthopoxviruses Pathogenic for Humans , 2010 .

[42]  E. Holmes,et al.  Rates of evolutionary change in viruses: patterns and determinants , 2008, Nature Reviews Genetics.

[43]  M. Danila,et al.  Complete coding sequences of the rabbitpox virus genome. , 2005, The Journal of general virology.

[44]  A. Hughes,et al.  Poxvirus genome evolution by gene gain and loss. , 2005, Molecular phylogenetics and evolution.

[45]  Geoffrey L. Smith Genus Orthopoxvirus: Vaccinia virus , 2007 .

[46]  Elliot J. Lefkowitz,et al.  Poxvirus protein evolution: Family wide assessment of possible horizontal gene transfer events , 2009, Virus Research.

[47]  D. Rock,et al.  Genome of Horsepox Virus , 2006, Journal of Virology.

[48]  J. Gulley,et al.  Poxviral vectors for cancer immunotherapy , 2012, Expert opinion on biological therapy.

[49]  David H. Evans,et al.  Genome Scale Patterns of Recombination between Coinfecting Vaccinia Viruses , 2014, Journal of Virology.

[50]  H. Bedson,et al.  Hybrids derived from the viruses of alastrim and rabbit pox , 1964, Journal of Hygiene.

[51]  K. Xing,et al.  Genome-Based Phylogeny of Poxvirus , 2006, Intervirology.

[52]  B. Moss,et al.  Immunology 101 at poxvirus U: immune evasion genes. , 2001, Seminars in immunology.

[53]  John P. Huelsenbeck,et al.  MrBayes 3: Bayesian phylogenetic inference under mixed models , 2003, Bioinform..

[54]  K. Crandall,et al.  Incorporating gaps as phylogenetic characters across eight DNA regions: ramifications for North American Psoraleeae (Leguminosae). , 2008, Molecular phylogenetics and evolution.

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

[56]  Y. Nakazawa,et al.  Chasing Jenner's Vaccine: Revisiting Cowpox Virus Classification , 2011, PloS one.

[57]  C. Upton,et al.  Prediction of Steps in the Evolution of Variola Virus Host Range , 2014, PloS one.

[58]  G. McFadden,et al.  Poxvirus host range genes. , 2008, Advances in virus research.

[59]  Daniel H. Huson,et al.  Dendroscope: An interactive viewer for large phylogenetic trees , 2007, BMC Bioinformatics.

[60]  S. Salzberg,et al.  Microbial gene identification using interpolated Markov models. , 1998, Nucleic acids research.

[61]  P. Earl,et al.  Elimination of A-type inclusion formation enhances cowpox virus replication in mice: Implications for orthopoxvirus evolution , 2014, Virology.

[62]  P. Brown,et al.  Comparative Analysis of Viral Gene Expression Programs during Poxvirus Infection: A Transcriptional Map of the Vaccinia and Monkeypox Genomes , 2008, PloS one.

[63]  A. Mercer,et al.  A truncated two-alpha-helix F-box present in poxvirus ankyrin-repeat proteins is sufficient for binding the SCF1 ubiquitin ligase complex. , 2009, The Journal of general virology.

[64]  M. V. Regenmortel,et al.  Virus taxonomy: classification and nomenclature of viruses. Seventh report of the International Committee on Taxonomy of Viruses. , 2000 .

[65]  Chris Upton,et al.  Poxvirus Bioinformatics Resource Center: a comprehensive Poxviridae informational and analytical resource , 2004, Nucleic Acids Res..

[66]  Rolf Apweiler,et al.  InterProScan: protein domains identifier , 2005, Nucleic Acids Res..

[67]  Brandon S. Gaut,et al.  Extensive gene gain associated with adaptive evolution of poxviruses , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[68]  H. Bedson,et al.  Hybrids derived from the viruses of variola major and cowpox , 1964, Journal of Hygiene.