Prediction of steps in the evolution of variola virus host range

Whole genome sequencing of recombinant viruses obtained from co-infection and superinfection of Vero cells with modified vaccinia virus ankara vectored influenza vaccine and a naturally occurring cowpox virus

Modified vaccinia virus Ankara (MVA) has been widely tested in clinical trials as recombinant vector vaccine against infectious diseases and cancers in humans and animals. However, one biosafety concern about the use of MVA vectored vaccine is the potential for MVA to recombine with naturally occurring orthopoxviruses in cells and hosts in which it multiplies poorly and, therefore, producing viruses with mosaic genomes with altered genetic and phenotypic properties. We previously conducted co-infection and superinfection experiments with MVA vectored influenza vaccine (MVA-HANP) and a feline Cowpox virus (CPXV-No-F1) in Vero cells (that were semi-permissive to MVA infection) and showed that recombination occurred in both co-infected and superinfected cells. In this study, we selected the putative recombinant viruses and performed genomic characterization of these viruses. Some putative recombinant viruses displayed plaque morphology distinct of that of the parental viruses. Our analysis demonstrated that they had mosaic genomes of different lengths. The recombinant viruses, with a genome more similar to MVA-HANP (>50%), rescued deleted and/or fragmented genes in MVA and gained new host ranges genes. Our analysis also revealed that some MVA-HANP contained a partially deleted transgene expression cassette and one recombinant virus contained part of the transgene expression cassette similar to that incomplete MVA-HANP. The recombination in co-infected and superinfected Vero cells resulted in recombinant viruses with unpredictable biological and genetic properties as well as recovery of delete/fragmented genes in MVA and transfer of the transgene into replication competent CPXV. These results are relevant to hazard characterization and risk assessment of MVA vectored biologicals.

Lumpy skin disease: history, current understanding and research gaps in the context of recent geographic expansion

Lumpy skin disease is recognized as a transboundary and emerging disease of cattle, buffaloes and other wild ruminants. Being initially restricted to Africa, and since 1989 the Middle East, the unprecedented recent spread across Eurasia demonstrates how underestimated and neglected this disease is. The initial identification of the causative agent of LSD as a poxvirus called LSD virus, was well as findings on LSDV transmission and epidemiology were pioneered at Onderstepoort, South Africa, from as early as the 1940s by researchers such as Weiss, Haig and Alexander. As more data emerges from an ever-increasing number of epidemiological studies, previously emphasized research gaps are being revisited and discussed. The currently available knowledge is in agreement with the previously described South African research experience that LSDV transmission can occur by multiple routes, including indirect contact, shared water sources and arthropods. The virus population is prone to molecular evolution, generating novel phylogenetically distinct variants resulting from a diverse range of selective pressures, including recombination between field and homologous vaccine strains in cell culture that produce virulent recombinants which pose diagnostic challenges. Host restriction is not limited to livestock, with certain wild ruminants being susceptible, with unknown consequences for the epidemiology of the disease.

Capripoxviruses, leporipoxviruses, and orthopoxviruses: Occurrences of recombination

Poxviruses are double-stranded DNA viruses with several members displaying restricted host ranges. They are genetically stable with low nucleotide mutation rates compared to other viruses, due to the poxviral high-fidelity DNA polymerase. Despite the low accumulation of mutations per replication cycle, poxvirus genomes can recombine with each other to generate genetically rearranged viruses through recombination, a process directly associated with replication and the aforementioned DNA polymerase. Orthopoxvirus replication is intimately tethered to high frequencies of homologous recombination between co-infecting viruses, duplicated sequences of the same virus, and plasmid DNA transfected into poxvirus-infected cells. Unfortunately, the effect of these genomic alterations on the cellular context for all poxviruses across the family Poxviridae remains elusive. However, emerging sequence data on currently circulating and archived poxviruses, such as the genera orthopoxviruses and capripoxviruses, display a wide degree of divergence. This genetic variability cannot be explained by clonality or genetic drift alone, but are probably a result of significant genomic alterations, such as homologous recombination, gene loss and gain, or gene duplications as the major selection forces acting on viral progeny. The objective of this review is to cross-sectionally overview the currently available findings on natural and laboratory observations of recombination in orthopoxviruses, capripoxviruses, and leporipoxviruses, as well as the possible mechanisms involved. Overall, the reviewed available evidence allows us to conclude that the current state of knowledge is limited in terms of the relevance of genetic variations across even a genus of poxviruses as well as fundamental features governing and precipitating intrinsic gene flow and recombination events.

Genomic Sequencing and Phylogenomics of Cowpox Virus

Cowpox virus (CPXV; genus Orthopoxvirus; family Poxviridae) is the causative agent of cowpox, a self-limiting zoonotic infection. CPXV is endemic in Eurasia, and human CPXV infections are associated with exposure to infected animals. In the Fennoscandian region, five CPXVs isolated from cats and humans were collected and used in this study. We report the complete sequence of their genomes, which ranged in size from 220–222 kbp, containing between 215 and 219 open reading frames. The phylogenetic analysis of 87 orthopoxvirus strains, including the Fennoscandian CPXV isolates, confirmed the division of CPXV strains into at least five distinct major clusters (CPXV-like 1, CPXV-like 2, VACV-like, VARV-like and ECTV-Abatino-like) and can be further divided into eighteen sub-species based on the genetic and patristic distances. Bayesian time-scaled evolutionary history of CPXV was reconstructed employing concatenated 62 non-recombinant conserved genes of 55 CPXV. The CPXV evolution rate was calculated to be 1.65 × 10⁻⁵ substitution/site/year. Our findings confirmed that CPXV is not a single species but a polyphyletic assemblage of several species and thus, a reclassification is warranted.

Recombination drives the emergence of orf virus diversity: evidence from the first complete genome sequence of an Indian orf virus isolate and comparative genomic analysis

Contagious pustular dermatitis is a disease that primarily infects small ruminants and possesses zoonotic potential. It is caused by orf virus (ORFV), a member of the genus Parapoxvirus. In this study, we evaluated an ORFV outbreak in goats in Madhya Pradesh, a state in central India, during 2017. The transboundary potential of this virus was evaluated by constructing phylogenetic trees. The complete genome sequence of an ORFV isolate named Ind/MP/17 was found to be 139,807 bp in length with 63.7% GC content and 132 open reading frames (ORFs) flanked by 3,910-bp inverted terminal repeats (ITRs). An investigation into evolutionary parameters such as selection pressure (θ = dN/dS) and nucleotide diversity (π) demonstrated that ORFV has undergone purifying selection. A total of 40 recombination events were identified, 21 of which were evident in the Ind/MP/17 genome, indicating its ability to generate new variants.

Genomic Sequencing and Analysis of a Novel Human Cowpox Virus With Mosaic Sequences From North America and Old World Orthopoxvirus

Orthopoxviruses (OPXVs) not only infect their natural hosts, but some OPXVs can also cause disease in humans. Previously, we partially characterized an OPXV isolated from an 18-year-old male living in Northern Norway. Restriction enzyme analysis and partial genome sequencing characterized this virus as an atypical cowpox virus (CPXV), which we named CPXV-No-H2. In this study, we determined the complete genome sequence of CPXV-No-H2 using Illumina and Nanopore sequencing. Our results showed that the whole CPXV-No-H2 genome is 220,276 base pairs (bp) in length, with inverted terminal repeat regions of approximately 7 kbp, containing 217 predicted genes. Seventeen predicted CPXV-No-H2 proteins were most similar to OPXV proteins from the Old World, including Ectromelia virus (ECTV) and Vaccinia virus, and North America, Alaskapox virus (AKPV). CPXV-No-H2 has a mosaic genome with genes most similar to other OPXV genes, and seven potential recombination events were identified. The phylogenetic analysis showed that CPXV-No-H2 formed a separate clade with the German CPXV isolates CPXV_GerMygEK938_17 and CPXV_Ger2010_MKY, sharing 96.4 and 96.3% nucleotide identity, respectively, and this clade clustered closely with the ECTV-OPXV Abatino clade. CPXV-No-H2 is a mosaic virus that may have arisen out of several recombination events between OPXVs, and its phylogenetic clustering suggests that ECTV-Abatino-like cowpox viruses form a distinct, new clade of cowpox viruses.

Teaching a new mouse old tricks: Humanized mice as an infection model for Variola virus

Smallpox, caused by the solely human pathogen Variola virus (VARV), was declared eradicated in 1980. While known VARV stocks are secure, smallpox remains a bioterrorist threat agent. Recent U.S. Food and Drug Administration approval of the first smallpox anti-viral (tecovirimat) therapeutic was a successful step forward in smallpox preparedness; however, orthopoxviruses can become resistant to treatment, suggesting a multi-therapeutic approach is necessary. Animal models are required for testing medical countermeasures (MCMs) and ideally MCMs are tested directly against the pathogen of interest. Since VARV only infects humans, a representative animal model for testing therapeutics directly against VARV remains a challenge. Here we show that three different humanized mice strains are highly susceptible to VARV infection, establishing the first small animal model using VARV. In comparison, the non-humanized, immunosuppressed background mouse was not susceptible to systemic VARV infection. Following an intranasal VARV challenge that mimics the natural route for human smallpox transmission, the virus spread systemically within the humanized mouse before mortality (~ 13 days post infection), similar to the time from exposure to symptom onset for ordinary human smallpox. Our identification of a permissive/representative VARV animal model can facilitate testing of MCMs in a manner consistent with their intended use.

Monkeypox virus emergence in wild chimpanzees reveals distinct clinical outcomes and viral diversity

Monkeypox is a viral zoonotic disease on the rise across endemic habitats. Despite the growing importance of monkeypox virus, our knowledge on its host spectrum and sylvatic maintenance is limited. Here, we describe the recent repeated emergence of monkeypox virus in a wild, human-habituated western chimpanzee (Pan troglodytes verus, hereafter chimpanzee) population from Taï National Park, Ivory Coast. Through daily monitoring, we show that further to causing its typical exanthematous syndrome, monkeypox can present itself as a severe respiratory disease without a diffuse rash. By analysing 949 non-invasively collected samples, we identify the circulation of at least two distinct monkeypox virus lineages and document the shedding of infectious particles in faeces and flies, suggesting that they could mediate indirect transmission. We also show that the carnivorous component of the Taï chimpanzees' diet, mainly consisting of the sympatric monkeys they regularly hunt, did not change nor shift towards rodent consumption (the presumed reservoir) before the outbreaks, suggesting that the sudden emergence of monkeypox virus in this population is probably due to changes in the ecology of the virus itself. Using long-term mortality surveillance data from Taï National Park, we provide evidence of little to no prior viral activity over at least two decades. We conclude that great ape sentinel systems devoted to the longitudinal collection of behavioural and health data can help clarify the epidemiology and clinical presentation of zoonotic pathogens.

Isolation and Characterization of Akhmeta Virus from Wild-Caught Rodents (Apodemus spp.) in Georgia

Akhmeta virus is a unique Orthopoxvirus that was described in 2013 from the country of Georgia. This paper presents the first isolation of this virus from small mammal (Rodentia; Apodemus spp.) samples and the molecular characterization of those isolates. The identification of the virus in small mammals is an essential component to understanding the natural history of this virus and its transmission to human populations and could guide public health interventions in Georgia. Akhmeta virus genomes harbor evidence suggestive of recombination with a variety of other orthopoxviruses; this has implications for the evolution of orthopoxviruses, their ability to infect mammalian hosts, and their ability to adapt to novel host species.

In 2013, a novel orthopoxvirus was detected in skin lesions of two cattle herders from the Kakheti region of Georgia (country); this virus was named Akhmeta virus. Subsequent investigation of these cases revealed that small mammals in the area had serological evidence of orthopoxvirus infections, suggesting their involvement in the maintenance of these viruses in nature. In October 2015, we began a longitudinal study assessing the natural history of orthopoxviruses in Georgia. As part of this effort, we trapped small mammals near Akhmeta (n = 176) and Gudauri (n = 110). Here, we describe the isolation and molecular characterization of Akhmeta virus from lesion material and pooled heart and lung samples collected from five wood mice (Apodemus uralensis and Apodemus flavicollis) in these two locations. The genomes of Akhmeta virus obtained from rodents group into 2 clades: one clade represented by viruses isolated from A. uralensis samples, and one clade represented by viruses isolated from A. flavicollis samples. These genomes also display several presumptive recombination events for which gene truncation and identity have been examined.

IMPORTANCE Akhmeta virus is a unique Orthopoxvirus that was described in 2013 from the country of Georgia. This paper presents the first isolation of this virus from small mammal (Rodentia; Apodemus spp.) samples and the molecular characterization of those isolates. The identification of the virus in small mammals is an essential component to understanding the natural history of this virus and its transmission to human populations and could guide public health interventions in Georgia. Akhmeta virus genomes harbor evidence suggestive of recombination with a variety of other orthopoxviruses; this has implications for the evolution of orthopoxviruses, their ability to infect mammalian hosts, and their ability to adapt to novel host species.

Bioinformatics for Analysis of Poxvirus Genomes

In recent years, there have been numerous technological advances in the field of molecular biology; these include next- and third-generation sequencing of DNA genomes and mRNA transcripts and mass spectrometry of proteins. Perhaps, however, it is genome sequencing that impacts a virologist the most. In 2017, more than 480 complete genome sequences of poxviruses have been generated, and are constantly used in many different ways by almost all molecular virologists. Matching this growth in data acquisition is an explosion of the relatively new field of bioinformatics, providing databases to store and organize this valuable/expensive data and algorithms to analyze it. For the bench virologist, access to intuitive, easy-to-use, software is often critical for performing bioinformatics-based experiments. Three common hurdles for the researcher are (1) selection, retrieval, and reformatting genomics data from large databases; (2) use of tools to compare/analyze the genomics data; and (3) display and interpretation of complex sets of results. This chapter is directed at the bench virologist and describes the software that helps overcome these obstacles, with a focus on the comparison and analysis of poxvirus genomes. Although poxvirus genomes are stored in public databases such as GenBank, this resource can be cumbersome and tedious to use if large amounts of data must to be collected. Therefore, we also highlight our Viral Orthologous Clusters database system and integrated tools that we developed specifically for the management and analysis of complete viral genomes.

Bayesian Phylogeography and Pathogenic Characterization of Smallpox Based on HA, ATI, and CrmB Genes

Variola virus is at risk of re-emergence either through accidental release, bioterrorism, or synthetic biology. The use of phylogenetics and phylogeography to support epidemic field response is expected to grow as sequencing technology becomes miniaturized, cheap, and ubiquitous. In this study, we aimed to explore the use of common VARV diagnostic targets hemagglutinin (HA), cytokine response modifier B (CrmB), and A-type inclusion protein (ATI) for phylogenetic characterization as well as the representativeness of modelling strategies in phylogeography to support epidemic response should smallpox re-emerge. We used Bayesian discrete-trait phylogeography using the most complete data set currently available of whole genome (n = 51) and partially sequenced (n = 20) VARV isolates. We show that multilocus models combining HA, ATI, and CrmB genes may represent a useful heuristic to differentiate between VARV Major and subclades of VARV Minor which have been associated with variable case-fatality rates. Where whole genome sequencing is unavailable, phylogeography models of HA, ATI, and CrmB may provide preliminary but uncertain estimates of transmission, while supplementing whole genome models with additional isolates sequenced only for HA can improve sample representativeness, maintaining similar support for transmission relative to whole genome models. We have also provided empirical evidence delineating historic international VARV transmission using phylogeography. Due to the persistent threat of re-emergence, our results provide important research for smallpox epidemic preparedness in the posteradication era as recommended by the World Health Organisation.

Molecular and microscopic characterization of a novel Eastern grey kangaroopox virus genome directly from a clinical sample

Poxviruses are large DNA viruses with varying zoonotic potential, and are recognised in a broad range of wildlife. Although poxviruses have been detected in kangaroos, their genetic relationships to poxviruses in other animals and humans is not well understood. Here, we present a novel genome sequence of a marsupial poxvirus, the Eastern grey kangaroopox virus (EKPV-NSW), isolated from a wild eastern grey kangaroo. In the present study, histopathologically confirmed epidermal pox lesions were used to recover the full-length viral genome and perform electron microscopic analysis, with both immature virions and intracellular mature virions detected. Subsequent analysis of the EKPV-NSW genome demonstrated the highest degree of sequence similarity with EKPV-SC strain (91.51%), followed by WKPV-WA (87.93%), and MOCV1 (44.05%). The novel EKPV-NSW complete genome encompasses most of the chordopoxviruses protein coding genes (138) that are required for genome replication and expression, with only three essential protein coding genes being absent. The novel EKPV-NSW is missing 28 predicted genes compared to the recently isolated EKPV-SC, and carries 21 additional unique genes, encoding unknown proteins. Phylogenetic and recombination analyses showed EKPV-NSW to be the distinct available candidate genome of chordopoxviruses.

Re-Assembly and Analysis of an Ancient Variola Virus Genome

We report a major improvement to the assembly of published short read sequencing data from an ancient variola virus (VARV) genome by the removal of contig-capping sequencing tags and manual searches for gap-spanning reads. The new assembly, together with camelpox and taterapox genomes, permitted new dates to be calculated for the last common ancestor of all VARV genomes. The analysis of recently sequenced VARV-like cowpox virus genomes showed that single nucleotide polymorphisms (SNPs) and amino acid changes in the vaccinia virus (VACV)-Cop-O1L ortholog, predicted to be associated with VARV host specificity and virulence, were introduced into the lineage before the divergence of these viruses. A comparison of the ancient and modern VARV genome sequences also revealed a measurable drift towards adenine + thymine (A + T) richness.

Two novel poxviruses with unusual genome rearrangements: NY_014 and Murmansk

The genome sequence and annotation of two novel poxviruses, NY_014 and Murmansk, are presented. Despite being isolated on different continents and from different hosts, the viruses are relatively similar, albeit distinct species. The closest known relative of the novel viruses is Yoka poxvirus. Five novel genes were found in these genomes, two of which were MHC class I homologs. Although the core of these genomes was well conserved, the terminal regions showed significant variability with large deletions and surprising evidence of recombination with orthopoxviruses.

Characterization of Two Historic Smallpox Specimens from a Czech Museum

Although smallpox has been known for centuries, the oldest available variola virus strains were isolated in the early 1940s. At that time, large regions of the world were already smallpox-free. Therefore, genetic information of these strains can represent only the very last fraction of a long evolutionary process. Based on the genomes of 48 strains, two clades are differentiated: Clade 1 includes variants of variola major, and clade 2 includes West African and variola minor (Alastrim) strains. Recently, the genome of an almost 400-year-old Lithuanian mummy was determined, which fell basal to all currently sequenced strains of variola virus on phylogenetic trees. Here, we determined two complete variola virus genomes from human tissues kept in a museum in Prague dating back 60 and 160 years, respectively. Moreover, mass spectrometry-based proteomic, chemical, and microscopic examinations were performed. The 60-year-old specimen was most likely an importation from India, a country with endemic smallpox at that time. The genome of the 160-year-old specimen is related to clade 2 West African and variola minor strains. This sequence likely represents a new endemic European variant of variola virus circulating in the midst of the 19th century in Europe.

Characterization of Eptesipoxvirus, a novel poxvirus from a microchiropteran bat

The genome of Eptesipoxvirus (EPTV) is the first poxvirus genome isolated from a microbat. The 176,688 nt sequence, which is believed to encompass the complete coding region of the virus, is 67% A+T and is predicted to encode 191 genes. 11 of these genes have no counterpart in GenBank and are therefore unique to EPTV. The presence of a distantly related ortholog of Vaccinia virus F5L in EPTV uncovered a link with fragmented F5L orthologs in Molluscum contagiosum virus/squirrelpox and clade II viruses. Consistent with the unique position of EPTV approximately mid-point between the orthopoxviruses and the clade II viruses, EPTV has 11 genes that are specific to the orthopoxviruses and 13 genes that are typical, if not exclusive, to the clade II poxviruses. This mosaic nature of EPTV blurs the distinction between the old description of the orthopoxvirus and clade II groups. Genome annotation and characterization failed to find any common virulence genes shared with the other poxvirus isolated from bat (pteropoxvirus); however, EPTV encodes 3 genes that may have been transferred to or from deerpox and squirrelpox viruses; 2 of these, a putative endothelin-like protein and a MHC class I-like protein are likely to have immunomodulatory roles.

Genomic characterization of two novel pathogenic avipoxviruses isolated from pacific shearwaters (Ardenna spp.)

Background

Over the past 20 years, many marine seabird populations have been gradually declining and the factors driving this ongoing deterioration are not always well understood. Avipoxvirus infections have been found in a wide range of bird species worldwide, however, very little is known about the disease ecology of avian poxviruses in seabirds. Here we present two novel avipoxviruses from pacific shearwaters (Ardenna spp), one from a Flesh-footed Shearwater (A. carneipes) (SWPV-1) and the other from a Wedge-tailed Shearwater (A. pacificus) (SWPV-2).

Results

Epidermal pox lesions, liver, and blood samples were examined from A. carneipes and A. pacificus of breeding colonies in eastern Australia. After histopathological confirmation of the disease, PCR screening was conducted for avipoxvirus, circovirus, reticuloendotheliosis virus, and fungal agents. Two samples that were PCR positive for poxvirus were further assessed by next generation sequencing, which yielded complete Shearwaterpox virus (SWPV) genomes from A. pacificus and A. carneipes, both showing the highest degree of similarity with Canarypox virus (98% and 67%, respectively). The novel SWPV-1 complete genome from A. carneipes is missing 43 genes compared to CNPV and contains 4 predicted genes which are not found in any other poxvirus, whilst, SWPV-2 complete genome was deemed to be missing 18 genes compared to CNPV and a further 15 genes significantly fragmented as to probably cause them to be non-functional.

Conclusion

These are the first avipoxvirus complete genome sequences that infect marine seabirds. In the comparison of SWPV-1 and −2 to existing avipoxvirus sequences, our results indicate that the SWPV complete genome from A. carneipes (SWPV-1) described here is not closely related to any other avipoxvirus genome isolated from avian or other natural host species, and that it likely should be considered a separate species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3680-z) contains supplementary material, which is available to authorized users.

The genomes of three North American orthopoxviruses

The complete genomes of a skunkpox, volepox, and raccoonpox virus were sequenced and annotated. Phylogenetic analysis of these genomes indicates that although these viruses are all orthopoxviruses, they form a distinct clade to the other known species. This supports the ancient divergence of the North American orthopoxviruses from other members of the orthopoxviruses. Only two open reading frames appear to be unique to this group of viruses, but a relatively small number of insertions/deletions contribute to the varied gene content of this clade. The availability of these genomes will help determine whether skunkpox and volepox viruses share the characteristics that make raccoonpox a useful vaccine vector.

The origin of the variola virus

The question of the origin of smallpox, one of the major menaces to humankind, is a constant concern for the scientific community. Smallpox is caused by the agent referred to as the variola virus (VARV), which belongs to the genus Orthopoxvirus. In the last century, smallpox was declared eradicated from the human community; however, the mechanisms responsible for the emergence of new dangerous pathogens have yet to be unraveled. Evolutionary analyses of the molecular biological genomic data of various orthopoxviruses, involving a wide range of epidemiological and historical information about smallpox, have made it possible to date the emergence of VARV. Comparisons of the VARV genome to the genomes of the most closely related orthopoxviruses and the examination of the distribution their natural hosts' ranges suggest that VARV emerged 3000 to 4000 years ago in the east of the African continent. The VARV evolution rate has been estimated to be approximately 2 × 10-6 substitutions/site/year for the central conserved genomic region and 4 × 10-6 substitutions/site/year for the synonymous substitutions in the genome. Presumably, the introduction of camels to Africa and the concurrent changes to the climate were the particular factors that triggered the divergent evolution of a cowpox-like ancestral virus and thereby led to the emergence of VARV.

Identification of nucleotide-level changes impacting gene content and genome evolution in orthopoxviruses

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.