Next Article in Journal
BST2/Tetherin Inhibition of Alphavirus Exit
Next Article in Special Issue
A Phylogeographic Investigation of African Monkeypox
Previous Article in Journal
Baculovirus-mediated Gene Delivery and RNAi Applications
Previous Article in Special Issue
Gene Acquisition Convergence between Entomopoxviruses and Baculoviruses
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome Variability and Gene Content in Chordopoxviruses: Dependence on Microsatellites

1
Department of Microbiology, University of Alabama at Birmingham, BBRB 276/11, 845 19th St S, Birmingham, AL 35222, USA
2
Stanford Genome Technology Center, Stanford University, 855 California Ave, Palo Alto, CA 94304, USA
*
Author to whom correspondence should be addressed.
Submission received: 27 February 2015 / Revised: 24 March 2015 / Accepted: 17 April 2015 / Published: 22 April 2015
(This article belongs to the Special Issue Poxvirus Evolution)

Abstract

:
To investigate gene loss in poxviruses belonging to the Chordopoxvirinae subfamily, we assessed the gene content of representative members of the subfamily, and determined whether individual genes present in each genome were intact, truncated, or fragmented. When nonintact genes were identified, the early stop mutations (ESMs) leading to gene truncation or fragmentation were analyzed. Of all the ESMs present in these poxvirus genomes, over 65% co-localized with microsatellites—simple sequence nucleotide repeats. On average, microsatellites comprise 24% of the nucleotide sequence of these poxvirus genomes. These simple repeats have been shown to exhibit high rates of variation, and represent a target for poxvirus protein variation, gene truncation, and reductive evolution.

1. Introduction

The Poxviridae family consists of viruses with a wide host range and potential for causing disease. The family is divided into subfamilies based on host range, consisting of the Entomopoxvirinae which infect insects, and the Chordopoxvirinae (ChPV) which infect vertebrates [1]. The chordopoxviruses are further divided into several genera, including the Orthopoxvirus genus, which includes the most well-known and the best characterized poxviruses, variola virus and vaccinia virus. Smallpox disease is caused by the viruses belonging to the Variola virus (VARV) species, and is well known both for its severity and its eradication through a worldwide vaccination effort. The vaccine used was based on viruses belonging to the species Vaccinia virus (VACV), which at some point replaced cowpox virus that was originally used in Jenner’s smallpox vaccine, although the natural host of vaccinia viruses is unknown [2]. In addition to the Orthopoxvirus genus, chordopoxviruses are categorized into nine other genera, namely Avipoxvirus, Capripoxvirus, Cervidpoxvirus, Crocodylidpoxvirus, Leporipoxvirus, Molluscipoxvirus, Parapoxvirus, Suipoxvirus, and Yatapoxvirus, as well as the unassigned species Squirrelpox virus, and the yet to be classified, yoka poxvirus (YKPV) [1,3].
Poxvirus genomes consist of linear double-stranded DNA of 133,000–360,000 bp bounded by inverted terminal repeats [4]. Poxviruses code for 133–328 genes, and at least 90 gene families exhibit significant homology across the ChPV subfamily [4]. Viruses assigned to the same genus exhibit higher degrees of gene conservation. For example, all orthopoxviruses share a core set of approximately 174 genes, however, the remaining genes outside of this core set are always present among the set of 214 genes present in cowpox viruses, in addition to the other orthopoxviruses that may also contain the gene [4]. Poxviruses are unique among the DNA viruses because DNA replication occurs in the cytoplasm as opposed to the nucleus [5,6]. The large coding capacity of their genomes is one of the ways poxviruses support replication in the absence of the cellular DNA replication machinery, by coding for many of their own transcription and replication proteins. Most poxvirus genomes share similar genetic organization, with genes involved in essential functions such as transcription, DNA replication, and virion assembly located in a highly conserved central region [7]. Genes involved in immunomodulation and host range determination tend to cluster at the ends of the genomes, which are much more variable in terms of sequence homology and gene content.
Viruses within the species Cowpox virus (CPXV) have the largest genomes in comparison to other orthopoxviruses, and are believed to be most genetically similar to the ancestral orthopoxvirus in terms of gene content [4,8]. Several species including Cowpox virus and Monkeypox virus (MPXV) are capable of causing zoonotic infections in humans, and both species have a wide host range [9]. The orthopoxviruses are antigenically similar, and infection with one can induce long-lasting immunity against the others [10]. Only limited genomic data is available for orthopoxviruses endemic to North America, however they appear to be phylogenetically distinct from other orthopoxviruses [11].
The Avipoxvirus genus infects birds, can be transmitted via arthropod vectors as well as aerosols, and members are being investigated as vaccine vectors [12]. Although most chordopoxviruses share similar genome organization, many avipoxviruses contain an extensive genomic rearrangement relative to other ChPV [13]. Chordopoxviruses are often named for the animal from which they were isolated, although in many cases it is not the reservoir host. Ruminants such as cattle, goats, and sheep are the hosts of viruses in the Capripoxvirus genus [14]. Members of the Cervidpoxvirus genus are normally found in mule deer, but have recently been isolated from a gazelle, viruses of the Leporipoxvirus species infect rabbits, members of the Suipoxvirus species infect only swine, and primates are the hosts for members of the Yatapoxvirus genus [15,16].
Some members of ChPV are distinct from the others due to their high genomic GC content. The GC content in chordopoxviruses ranges from 25% in the capripoxviruses to 67% in squirrelpox virus [17], although a biological role for variation in GC content has not yet been identified. The high GC poxviruses include viruses in the Molluscipoxvirus genus, which is specific to humans, Nile crocodilepox virus, a virus in the Crocodylidpoxvirus genus, and members of the Parapoxvirus genus, which have been identified in hoofed mammals, shrews, weasels, and other members of the mammalian Laurasiatheria superorder [18,19,20,21].
The large host range of some poxviruses may be a result of host jumping, where a virus will occasionally infect an animal other than its typical host, and adaptive mutations may allow successful colonization of the new species [22,23]. DNA viruses have a comparatively low genetic mutation rate of approximately 10−7–10−9 mutations per site per round of replication which is closer to that of their hosts than to RNA viruses [24,25]. It is not well understood how the potentially large number of mutations necessary to cross host species barriers are able to accumulate quickly enough to allow the virus to adapt to a new host. Poxviruses are subject to several mechanisms which introduce genomic variability, including horizontal gene transfer (HGT), single nucleotide polymorphisms (SNPs), small insertions and deletions (indels), larger deletions of sequence which may include entire genes, and recombination [26]. The mechanisms utilized by viruses within any particular poxvirus genus to support their variation and evolution may differ from one genus to another. For example, the orthopoxviruses do not seem to have acquired any new genes through HGT since their most recent common ancestor, however, evidence supports progressive gene fragmentation and loss as one of the major determinants of speciation within the Orthopoxvirus genus [4]. Short indels are the most common cause for variation in orthopoxvirus gene content resulting in gene truncation, fragmentation, and loss [27].
Microsatellites, also called Simple Sequence Repeats (SSRs), or Short Tandem Repeats (STRs), are motifs of 1–6 nucleotides arranged in tandem repeats [28,29]. They are present in eukaryotes, prokaryotes, and viruses. The presence and characteristics of microsatellites have been shown to vary according to genome length and GC content, and the number and length of repeats varies by species and their location in coding and noncoding regions [30,31]. Microsatellites can function as hypervariable regions compared to surrounding non-repetitive genomic sequence. The motifs are not always replicated perfectly during DNA replication, which can lead to changes in the number of repetitions or single nucleotide polymorphisms (SNPs). Slipped-strand mispairing is one of the most likely causes of mutation in microsatellites, and occurs when DNA polymerase “slips” backwards or forwards along the template, and is facilitated by the short repeats [30,32]. When insertions and deletions (indels) occur in gene coding regions, they result in frameshift mutations and frequently therefore early stop mutations and truncation, unless the indel length is an in-frame multiple of 3, in which case the protein sequence will vary by the deletion or addition of an amino acid. Microsatellites can also serve as “hot spots” for recombination, which can lead to duplication or deletion of larger stretches of DNA, sometimes including whole genes [33,34].
Repetitive genetic elements have been shown to mediate different aspects of virus biology such as latency in some herpesviruses [35]. Poxviruses contain clusters of tandem repeats in their telomeres near the exterior ends of inverted terminal repeats bookending the genome, although the repeat units are longer than microsatellites [36]. The high variability often identified at short sequence repeats has led to great interest in characterizing microsatellites in many different species. Microsatellite distribution has been investigated in ssDNA viruses such as members of the family Geminiviridae that infect plants and the families Circoviridae, Parvoviridae, and Anelloviridae that infect vertebrates [37,38]. Short repeats were also identified throughout viruses in the Herpesviridae family, six genotypes of hepatitis C viruses, adenoviruses, influenza viruses, sin nombre virus, and human immunodeficiency virus type 1 [39,40,41,42,43].
Despite widespread identification of microsatellites and their functional significance in eukaryotic and prokaryotic genomes, the biological role of microsatellites in viral genomes remains largely unknown. We previously identified an inverse association between the length of orthopoxvirus genomes and the number of mutations they contain that could lead to truncated and fragmented genes [27]. While analyzing the data, we observed that many of the ESMs occurred at microsatellites (data was not reported). In this report, we have assessed the microsatellite content of viruses throughout the Chordopoxvirinae subfamily, characterized the microsatellites present in the viral genomes, and analyze the relationship between ESMs and microsatellites.

2. Materials and Methods

2.1. Genome Sequences and Assessment of Gene Content

We chose representative isolates from each genus within the Chordopoxvirinae subfamily, as well as each species within the Orthopoxvirus genus with available complete genome sequences (Table 1). Nucleotide sequences were downloaded from the Viral Bioinformatics Resource Center (vbrc.org) [44], with the exception of SQPV, which was downloaded from the National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov).
Table 1. Selected chordopoxvirus genomes used in these analyses.
Table 1. Selected chordopoxvirus genomes used in these analyses.
GenusAbbreviationNCBI strainAccession no.Intact genesTruncated genesFragmented genesLength (bp)
Species
Avipoxvirus
Fowlpox virusFWPVFowlpox virus strain IowaNC_002188231513288,539
Canarypox virusCNPVCanarypox virus strain ATCC VR111NC_00530931624359,853
Capripoxvirus
Sheeppox virusSPPVSheeppox virus strain TU-V02127NC_00400215018149,955
Goatpox virusGTPVGoatpox virus strain PellorNC_00400315027149,599
Lumpy skin disease virusLSDVLumpy skin disease virus strain Neethling 2490NC_00302715900150,773
Cervidpoxvirus
Mule deerpox virusDPVDeerpox virus W-848-83NC_00696617201166,259
Crocodylidpoxvirus
Nile crocodilepox virusCRVCrocodilepox virus strain ZimbabweNC_00803017500190,054
Leporipoxvirus
Myxoma virusMYXVMyxoma virus strain LausanneNC_00113216601161,773
Rabbit fibroma virusRFVRabbit fibroma virus strain KaszaNC_00126616104159,857
Molluscipoxvirus
Molluscum contagiosum virusMOCVMolluscum contagiosum virus strain subtype 1NC_00173115900190,289
Orthopoxvirus
Cowpox virusCPXV-GerCowpox virus strain Germany 91-3DQ43759320833228,250
Cowpox virusCPXV-BRCowpox virus strain Brighton RedAF48275820631224,499
Cowpox virusCPXV-GriCowpox virus strain GRI-90X9435520930223,666
Vaccinia virusHSPVHorsepox virus strain MNR-76DQ792504181224212,633
Vaccinia virusRPXVRabbitpox virus strain UtrechtAY484669179134197,731
Vaccinia virusVACV-WRVaccinia virus strain WR (Western Reserve)AY243312178123194,711
Monkeypox virusMPXV-WRMonkeypox virus strain MPXV-WRAIR7-61; Walter Reed 267AY603973175714199,195
Monkeypox virusMPXV-ZAIMonkeypox virus strain Zaire-96-I-16AF380138176714196,858
Ectromelia virusECTVEctromelia virus strain MoscowAF012825172217209,771
Taterapox virusTATVTaterapox virus strain Dahomey 1968DQ437594163267198,050
Camelpox virusCMLVCamelpox virus strain M-96 from KazakhstanAF438165174148205,719
Variola virusVARV-BRZVariola virus strain Brazil 1966 (v66-39 Sao Paulo)DQ4414191621812188,062
Variola virusVARV-SLEVariola virus strain Sierra Leone 1969 (V68-258)DQ4414371621711187,014
Variola virusVARV-SAFVariola virus strain South Africa 1965 (103 T'vaal, Nelspruit)DQ4414361621711185,881
Variola virusVARV-KWTVariola virus strain Kuwait 1967 (K1629)DQ4414331621711185,853
Parapoxvirus
Orf virusORFVOrf virus strain OV-SA00NC_00533613400139,962
Bovine papular stomatitis virusBPSVBovine papular stomatitis virus strain BV-AR02NC_00533713600134,431
Suipoxvirus
Swinepox virusSWPVSwinepox virus strain Nebraska 17077-99NC_00338915100146,454
Yatapoxvirus
Tanapox virusTANVTanapox virus strain KenyaNC_00988815120144,565
Yaba monkey tumor virusYMTVYaba monkey tumor virus strain AmanoNC_00517913920134,721
Unassigned
Yoka poxvirus (not classified)YKPVYoka poxvirusNC_015960.1164125175,699
Squirrelpox virusSQPVSquirrel poxvirus strain Red squirrel UKNC_022563.114200148,803
As described in Hendrickson et al. [4] and Hatcher et al. [27], we used the Poxvirus Genome Annotation System (PGAS) to determine the coding potential of open reading frames (ORFs), to annotate predicted genes according to their coding state (i.e., intact, truncated, fragmented, or missing genes), and to compare syntenic regions in closely related strains or species. MAFFT alignments using the FFT-NS-2 method were generated for the avipoxviruses, the high GC isolates, and all isolates except the avipoxviruses, and were used to identify genes that have similar location but low homology [45]. An ORF was considered intact if it was similar in length to the longest ORF for all orthologs of that gene. Truncated genes were defined as genes greater than 30 amino acids long and retaining a predicted promoter sequence, but were less than 80% of the amino acid length of intact orthologs. This cutoff is based on a comparison of gene length conservation (described in Hendrickson et al. [4]). Fragmented genes were identified as genomic regions containing identifiable homologous sequence to an orthologous, intact ORF, but with the fragmented ORF coding for less than 30 amino acids, missing a predicted promoter region, or missing the 5' end of the gene, including a start codon.
ORFs that were flagged as potential genes through PGAS automated prediction but did not have an identified promoter, had an absent or weak Kozak consensus sequence, had a low Glimmer score, and showed no orthology with intact poxvirus genes were not labeled as genes or gene remnants unless they were identified as genes when the genome was originally published. To differentiate orthologous from paralogous genes, we assessed sequence homology as well as gene synteny, the conservation of genomic location and gene neighbors [46]. We utilized the term “syntelog” to represent orthologous genes shared across isolates with a common syntenic genome location [4].

2.2. Detection of Early Stop Mutations (ESMs)

ESMs are defined as mutations that give rise to a stop codon in the sequence of a gene that either interrupts the start codon or truncates the ORF to a length 80% or less of the length of the intact orthologous gene. ESMs were identified as in Hatcher et al. [27] by visualizing the nucleotide sequence of each open reading frame and the corresponding amino acid translation of the coding frame and the two alternate frames, and annotating mutations that introduced either a stop codon in the normal coding frame, introduced an indel resulting in a frameshift mutation and subsequent early stop codon in the new reading frame, or altered the start codon of a gene. In the case of an altered reading frame, the mutation that caused the change in reading frame is coded as the ESM and not the newly introduced (now in-frame) stop codon.

2.3. Microsatellite Identification

Each of the viruses was assessed individually for the presence of microsatellites using the IMEx program available at http://imex.cdfd.org.in [42]. Repeats with motifs of 1–6 bp were identified if present with a minimum repeat number of 4 for motifs of 1 base, 3 for motifs of 2 or 3 bases, and 2 for motifs of 4, 5, or 6 bases. Perfect and imperfect repeats were allowed, and the imperfections were limited to a variation of 10% within the repeat tract. Microsatellite content was determined by dividing the total number of nucleotides present within the repeats by the length of the genome. When counting ESMs that overlapped microsatellites, if an ESM only partially overlapped a microsatellite, it was not included in the number of ESMs overlapping repeats.

2.4. Statistical Analysis

The total number of microsatellites was normalized according to genome size. Relative abundance was calculated as the number of microsatellites per kilobase of genomic sequence. Relative density was calculated as the sum total of microsatellite nucleotides per kilobase of genomic sequence.
Linear regression was used to determine the relationship between various genomic features and microsatellites, and was performed in Microsoft Excel 2010.

3. Results

3.1. Assessment of Chordopoxvirus Gene Content

We annotated the genes of representative viruses from each chordopoxvirus genus, and characterized those genes as intact, truncated, or fragmented (Table 1). Using whole genome multiple sequence alignments, we were able to use homology and synteny to assist with gene annotation, including identifying truncated and fragmented genes by comparing the genomic sequence in these regions with the sequences of intact, syntenic genes present in other isolates. When identifying non-intact genes, we define truncated genes as less than 80% the length of intact genes but with the 5' end intact, as opposed to fragmented genes which maintain nucleotide sequence homology to intact genes but either do not have an intact 5' end of the gene, or the remaining orthologous ORF is less than 30 amino acids long. We were able to identify at least one non-intact gene in most of the viruses, although the majority were in orthopoxviruses (reported in Hatcher et al. [27]), yoka poxvirus, capripoxviruses, and the avipoxviruses. Yoka poxvirus is the closest phylogenetic relative to the orthopoxviruses which has been identified at this time [3]. We were unable to identify any truncated or fragmented genes in the high GC viruses. Figure 1 provides a comparison of the similarities of the isolates as pairwise nucleotide identities for the most conserved central portion of the genomes.
Figure 1. Pairwise nucleotide % identities for the conserved central region of the poxvirus isolates used in this study. Values below the diagonal included gaps when calculating the pairwise identity, whereas values above the diagonal had all gapped columns removed before analysis. The conserved central region spans the genomic sequences from Cop-F9L (CPXV-Ger position 52,398) to Cop-A34R (CPXV-Ger position 160,094). Darker shading reflects higher levels of identity. *FWPV and CNPV, contain a genome inversion in this region, and therefore their sequences were rearranged in order to align them with the other viruses.
Figure 1. Pairwise nucleotide % identities for the conserved central region of the poxvirus isolates used in this study. Values below the diagonal included gaps when calculating the pairwise identity, whereas values above the diagonal had all gapped columns removed before analysis. The conserved central region spans the genomic sequences from Cop-F9L (CPXV-Ger position 52,398) to Cop-A34R (CPXV-Ger position 160,094). Darker shading reflects higher levels of identity. *FWPV and CNPV, contain a genome inversion in this region, and therefore their sequences were rearranged in order to align them with the other viruses.
Viruses 07 02126 g001

3.2. Assessment of Truncated and Fragmented Genes, and Early Stop Mutation Content

Genome reduction through gene fragmentation and loss is one of the primary drivers of speciation in the orthopoxviruses, and previously we described an analysis of early stop mutations (ESMs) leading to gene truncation and fragmentation in orthopoxviruses [27]. We expanded that analysis here, characterizing the ESMs in viruses in the other genera of the Chordopoxvirinae subfamily, and were able to identify the causative ESMs for all but two of the non-intact genes (Table 2). The two instances were due to high variability between the non-intact genes and their homologs, which prevented us from being able to reliably identify the ESMs.
Orthopoxvirus genes range from highly intact in cowpox viruses (0.013–0.026 non-intact genes per 1000 bp) to highly fragmented in the variola viruses and TATV (0.150–0.167 non-intact genes per 1000 bp). Orthopoxvirus truncated and fragmented genes are interrupted by an average of 2.8 ESMs per gene, similar to the average for the entire subfamily, although some isolates average as high as five ESMs per non-intact gene. YKPV is not only genetically similar to the orthopoxviruses, but has a comparable number of truncated or fragmented genes (17 in YKPV and an average of 20.7 in orthopoxviruses), and ESMs per nonintact gene (2.8 in both YKPV and orthopoxviruses).
Table 2. Comparative genomic content of early stop mutationss and microsatellites.
Table 2. Comparative genomic content of early stop mutationss and microsatellites.
Virus GenomeESM CountNon-intact genesNon-intact genes/1000 bpESMs/non-intact geneESMs overlapping microsatellites% of ESMs overlapping microsatellitesGenome microsatellite content (%)GC% of genomeGC% of microsatellites
FWPV53180.0622.91528.324.03115
CNPV2260.0173.7940.923.53022
SPPV10790.06011.86661.731.7259
GTPV8290.0609.15668.331.72510
LSDV000-0-31.12610
DPV210.0060.0150.029.72610
CRV000-0-22.26266
MYXV210.0061.0150.020.84435
RFV4440.02511.01125.021.14028
MOCV000-0-22.56370
CPXV-Ger1160.0261.8763.622.53420
CPXV-BR340.0180.8133.322.93320
CPXV-Gri230.0130.7150.022.53420
HSPV49260.1221.92959.222.63319
RPXV44170.0862.63681.822.53419
VACV-WR41150.0772.73175.622.83319
MPXV-WR104210.1055.06663.522.73319
MPXV-ZAI95210.1074.55658.922.63319
ECTV73280.1332.63750.722.53319
TATV87330.1672.65563.222.73319
CMLV65220.1073.03756.922.73319
VARV-BRZ92300.1603.16065.222.63318
VARV-SLE100280.1503.66767.022.63318
VARV-SAF93280.1513.35862.422.63318
VARV-KWT93280.1513.35963.422.73318
ORFV000-0-22.86468
BPSV000-0-12.26576
SWPV000-0-27.12812
TANV420.0141.0125.030.22710
YMTV420.0151.5375.028.03012
YKPV47170.0972.82144.728.8 *2611
SQPV000-00.023.56773
Total11403790.060*3.4 *78468.8 *24.137 *26 *
* These figures are averages, not totals.
In the avipoxviruses, FWPV has many more truncated and fragmented genes than CNPV (18 versus 6), a greater number of ESMs (53 versus 22), but the frequency of ESMs per nonintact gene in FWPV is lower than in CNPV (2.9 versus 3.7). Two members of the capripoxviruses have nine non-intact genes each, with an average of 11.8 ESMs per gene in SPPV, and an average of 9.1 in GTPV. The third capripoxvirus isolate included, LSDV, does not appear to have any truncated or fragmented genes. Because we were unable to identify truncated or fragmented genes in the high GC content viruses, we were also unable to identify ESMs in those viruses. It is difficult to say whether the lack of non-intact genes may be due to the low number of sequenced isolates available for high GC viruses compared to high AT viruses, or due to different evolutionary mechanisms that prevent the accumulation of ESMs. The availability of additional sequenced high GC poxvirus isolates may provide the data necessary to answer this question.
Short repeats in DNA sequences are associated with higher rates of mutation in many species, and we wanted to understand if microsatellites may have a role in the introduction of early stop mutations in poxviruses. To determine if there is an association between ESMs and short repeats, we identified the microsatellites in the genomes, and identified the ESMs that occurred at the microsatellites (Figure 2). The short repeats are present in all of the isolates and account for an average 24% of the genome in chordopoxviruses, ranging from 12.2% in BPSV to 31.7% in the capripoxviruses (Table 2). The orthopoxviruses are within a very narrow range of 22.4%–23.0% of the genome consisting of repeats. In most viruses where ESMs were identified, a greater number of ESMs are located at microsatellites than can be explained by the microsatellite content of the genomes (Table 2 and Figure 3A,B). In some examples, the incidence of colocalization is striking, for example in RPXV, 81.1% of the ESMs occur at microsatellites, but the microsatellites comprise only 22.5% of the genome. This is not a universal observation however since in FWPV, RFV, and TANV, the distributions of ESMs and the genome microsatellite content are close to equal. We did not, however, find any correlation between the frequency of ESMs overlapping microsatellites, and the overall genome microsatellite content (Figure 3C).
Figure 2. Example of an early stop mutation (ESM) overlapping a microsatellite. For each isolate, a portion of the coding sequence for the thymidylate kinase gene (Cop-A48R) is shown. The top row is the nucleotide sequence, the amino acid translation is the middle row, and a bar is shown to indicate the gene ORF. The green bar represents the intact gene in CPXV-Ger which continues off the right of the figure, and the grey bar represents the homologous truncated gene in TATV. The pink bars highlight microsatellites consisting of a monomeric repeat. An ESM in TATV in the form of a 2 bp deletion is shown by the red bar.
Figure 2. Example of an early stop mutation (ESM) overlapping a microsatellite. For each isolate, a portion of the coding sequence for the thymidylate kinase gene (Cop-A48R) is shown. The top row is the nucleotide sequence, the amino acid translation is the middle row, and a bar is shown to indicate the gene ORF. The green bar represents the intact gene in CPXV-Ger which continues off the right of the figure, and the grey bar represents the homologous truncated gene in TATV. The pink bars highlight microsatellites consisting of a monomeric repeat. An ESM in TATV in the form of a 2 bp deletion is shown by the red bar.
Viruses 07 02126 g002
Figure 3. Presence of early stop mutations at microsatellites is independent of repeat content of the genome. (a) The number of ESMs which overlap or do not overlap microsatellites. Only isolates with non-intact genes and ESMs are shown; (b) The percentage of ESMs which overlap microsatellites is shown, along with the percentage of each genome which consists of microsatellites; (c) No relationship could be identified between the percentages of ESMs which overlap microsatellites and the percentage of each genome which consists of microsatellites.
Figure 3. Presence of early stop mutations at microsatellites is independent of repeat content of the genome. (a) The number of ESMs which overlap or do not overlap microsatellites. Only isolates with non-intact genes and ESMs are shown; (b) The percentage of ESMs which overlap microsatellites is shown, along with the percentage of each genome which consists of microsatellites; (c) No relationship could be identified between the percentages of ESMs which overlap microsatellites and the percentage of each genome which consists of microsatellites.
Viruses 07 02126 g003aViruses 07 02126 g003b

3.3. Chordopoxvirus Microsatellite Content

Counts of microsatellites of oligonucleotides from 1–6 bp motifs are described in Table 3 and Figure 4. The relative abundance (RA) of microsatellites, or the average number of microsatellites in every kilobase of genomic sequence varies from 18 in BSPV to 55 in SPPV and GTPV. The orthopoxviruses are very similar to each other, and each member has a RA between 34 and 39 microsatellites/kb. As the length of repeat motifs increases, the number of repeats decreases. Although the reported hexamer RA for several of the viruses is close to 0, repeats of 6 bp motifs were observed in all of the viruses. The portion of the genome made up of microsatellites is the relative density (RD), and ranges from 122 microsatellite nucleotides per kilobase of genomic nucleotides in BSPV to 317 nt/kb in SPPV and GTPV. The most numerous repeats are monomers and trimers, with the greatest number of repeats being monomers in high AT viruses, with high GC viruses having greater numbers of trimers. This is true whether comparing actual counts or RA for each virus, however because RD takes into account the size of each repeat, there is no clear pattern. Microsatellites appear to be distributed relatively evenly across genomes, as shown in Figure 5; however, there is a significant difference between the RA of microsatellites in noncoding and coding regions (T test, p < 0.0001, 95% CI of the difference between the means is 10.146–26.034). The mean RA in noncoding regions is 57 microsatellites/1000 bp, and in coding regions is 39 microsatellites/1000 bp. This is best explained by more relaxed selection pressure on noncoding regions, where mutations due to the increased variability in microsatellites are less likely to affect virus viability.
Table 3. Description of the microsatellite sequences identified in Poxviridae genomes. Comparative abundance of repeats is shown. The high GC viruses are shown at the bottom of the table.
Table 3. Description of the microsatellite sequences identified in Poxviridae genomes. Comparative abundance of repeats is shown. The high GC viruses are shown at the bottom of the table.
IsolateTotal numberTotal number of 1, 2, 3, 4, 5, 6 bp
microsatellites
RARA of 1, 2, 3, 4, 5, 6 bp
Microsatellites a
Total length of microsatellitesRDRD of 1, 2, 3, 4, 5, 6 bp
microsatellites
RA in noncoding regionsRA in coding regions
FWPV120635896, 894, 4110, 838, 217, 1084220, 3, 14, 3, 1, 0a6935924093, 20, 89, 26, 8, 54940
CNPV144656687, 1100, 5214, 999, 310, 1474019, 3, 14, 3, 1, 0a8450823585, 20, 91, 25, 9, 54739
SPPV82834964, 435, 2207, 461, 143, 735533, 3, 15, 3, 1, 0a47574317162, 19, 93, 28, 10, 68553
GTPV82144863, 431, 2209, 481, 154, 765533, 3, 15, 3, 1, 147406317159, 19, 92, 29, 11, 79152
LSDV81814888, 430, 2173, 463, 149, 785432, 3, 14, 3, 1, 146886311157, 19, 91, 28, 10, 712951
DPV86204780, 583, 2490, 510, 154, 1035229, 4, 15, 3, 1, 149369297134, 23, 94, 28, 10, 88149
MYXV57802599, 373, 2174, 464, 123, 473616, 2, 13, 3, 1, 0a3365220873, 15, 83, 26, 8, 48833
RFV57812683, 423, 2018, 469, 130, 583617, 3, 13, 3, 1, 0a3380021176, 17, 78, 25, 9, 66234
CPXV-Ger86483650, 704, 3307, 698, 173, 1163816, 3, 14, 3, 1, 15141522573, 20, 91, 27, 8, 76235
CPXV-BR86563758, 716, 3178, 685, 172, 1473917, 3, 14, 3, 1, 15144322976, 20, 88, 27, 8, 95736
CPXV-Gri84663613, 691, 3202, 668, 178, 1143816, 3, 14, 3, 1, 15025122573, 20, 89, 27, 8, 75336
HSPV82053662, 632, 3012, 643, 158, 983917, 3, 14, 3, 1, 0a4810922679, 19, 88, 27, 8, 64537
RPXV75703370, 584, 2766, 596, 154, 1003817, 3, 14, 3, 1, 14444722578, 19, 87, 27, 8, 74437
VACV-WR75863506, 588, 2664, 574, 162, 923918, 3, 14, 3, 1, 0a4432322883, 19, 85, 26, 9, 64537
MPXV-WR76453281, 643, 2854, 626, 155, 863816, 3, 14, 3, 1, 0a4518722775, 20, 89, 28, 8, 64437
MPXV-ZAI75403262, 617, 2830, 589, 157, 853817, 3, 14, 3, 1, 0a4455922675, 20, 90, 26, 9, 64637
ECTV80183519, 667, 2942, 651, 153, 863817, 3, 14, 3, 1, 0a4722022576, 20, 87, 28, 8, 64536
TATV75593273, 582, 2880, 575, 157, 923817, 3, 15, 3, 1, 0a4494322775, 19, 91, 25, 9, 84237
CMLV79003446, 644, 2945, 607, 159, 993817, 3, 14, 3, 1, 0a4668822777, 20, 89, 26, 8, 84437
VARV-BRZ73053294, 558, 2665, 564, 146, 783918, 3, 14, 3, 1, 0a4251322680, 19, 88, 26, 8, 54338
VARV-SLE72413237, 554, 2667, 561, 145, 773917, 3, 14, 3, 1, 0a4225122679, 19, 89, 26, 8, 54337
VARV-SAF72143234, 553, 2650, 556, 143, 783917, 3, 14, 3, 1, 0a4207822679, 19, 89, 26, 8, 54337
VARV-KWT72283238, 553, 2658, 557, 144, 783917, 3, 14, 3, 1, 0a4214022779, 19, 89, 26, 8, 54338
SWPV66423001, 672, 2296, 467, 205, 744520, 5, 16, 3, 1, 13973727196, 31, 98, 29, 16, 76844
TANV77464835, 371, 1929, 415, 117, 795433, 3, 13, 3, 1, 143730302162, 17, 83, 26, 8, 76752
YMTV67034184, 322, 1668, 379, 93, 575031, 2, 12, 3, 1, 0a37689280150, 15, 77, 25, 7, 56348
YKPV85173902, 821, 2973, 554, 180, 874822, 5, 17, 3, 1, 0a50540288103, 31, 107, 29, 11, 76047
CRV66261711, 640, 3533, 460, 139, 145359, 3, 19, 2, 1, 14220922240, 24, 119, 22, 8, 106234
MOCV64561475, 1274, 2960, 454, 130, 163348, 7, 16, 2, 1, 14285222534, 52, 98, 21, 7, 124132
ORFV49591123, 748, 2540, 335, 94, 119358, 5, 18, 2, 1, 13194022835, 38, 116, 22, 7, 113036
BSPV24041123, 522, 240, 318, 80, 121188, 4, 2, 2, 1, 11643712236, 27, 20, 22, 7, 113417
SQPV55651423, 393, 3159, 337, 85, 1683710, 3, 21, 2, 1, 13491123541, 18, 135, 20, 6, 145636
a The RA of 6bp microsatelites is closer to 0 than to 1, but is not equal to 0.
Figure 4. Relative abundance and density of microsatellites. RA for (a) all microsatellites and for (b) each motif length measured (1 bp to 6 bp) for each isolate. RD for (c) all microsatellites and for (d) each motif length measured for each isolate. The high GC viruses are shown at the right end of each figure.
Figure 4. Relative abundance and density of microsatellites. RA for (a) all microsatellites and for (b) each motif length measured (1 bp to 6 bp) for each isolate. RD for (c) all microsatellites and for (d) each motif length measured for each isolate. The high GC viruses are shown at the right end of each figure.
Viruses 07 02126 g004
Figure 5. Distribution of microsatellites across genomes. The number of microsatellites is shown for bins of 2000 bp. The number of bins reflects the length of each genome. The isolates shown are (a) CNPV; (b) CPXV-Ger; (c) GTPV; (d) ORFV; and (e) VARV-BRZ.
Figure 5. Distribution of microsatellites across genomes. The number of microsatellites is shown for bins of 2000 bp. The number of bins reflects the length of each genome. The isolates shown are (a) CNPV; (b) CPXV-Ger; (c) GTPV; (d) ORFV; and (e) VARV-BRZ.
Viruses 07 02126 g005
We next attempted to identify genomic characteristics that may be associated with microsatellite content. The genome length for viruses included in this study ranges from approximately 134,000–360,000 bp. Neither RA nor RD show significant correlation with genome size in the chordopoxvirus isolates we analyzed (Figure 6). We also analyzed the microsatellites for their GC content, since it has previously been shown that the GC content of the genome is associated with the microsatellite composition in other organisms [31,47]. All of the microsatellites in a virus were assessed for the proportion of G and C nucleotides they contained, and this number was compared to the total GC content of each genome. In these poxviruses, the GC content of the microsatellites was strongly associated with the GC content of the genome (Figure 7).
Figure 6. Genome length is not correlated with microsatellite abundance or density. (a) RA and (b) RD are shown for each isolate, plotted versus genomic length.
Figure 6. Genome length is not correlated with microsatellite abundance or density. (a) RA and (b) RD are shown for each isolate, plotted versus genomic length.
Viruses 07 02126 g006
Figure 7. Correlation of microsatellite GC content and genome GC content.
Figure 7. Correlation of microsatellite GC content and genome GC content.
Viruses 07 02126 g007

4. Discussion

Poxviruses have been shown to evolve through a combination of vertical gene transfer, horizontal gene transfer from eukaryotes and other viruses, genetic drift, and, at least in the orthopoxviruses, gene loss. In this report, we cataloged all of the truncated and fragmented genes present in completely sequenced genomes of representative poxviruses from the Chordopoxvirinae subfamily, and we identified the early stop mutations leading to these gene truncations and fragmentations. We also showed that these ESMs were associated with microsatellites, and we characterized the microsatellite content in the genomes.
We observed a wide variation in the number of nonintact genes both in the subfamily as a whole, and often within genera. While reductive evolution has commonly been accepted as a major form of adaptation in the orthopoxviruses, this finding may signify that it is also active in other poxvirus genera such as the avipoxviruses. Although FWPV and CNPV share a similar genomic arrangement, CNPV contains over 75 kbp of additional sequence compared to FWPV [48]. FWPV also contains three times as many nonintact genes compared to CNPV (Table 2). Therefore, genome reduction through a process of gene truncation and fragmentation may have played a significant role in the evolution of FWPV. Similar features can be seen in the genus Capripoxvirus, where LSDV appears to only code for intact genes. In comparison, SPPV and GTPV each have nine non-intact genes.
One of the limitations with this study is that gene truncations and fragments cannot be identified without an example of an intact gene. When analyzing genomes for genes that are truncated or fragmented, we were unable to identify any in the high GC viruses. There are at least two possible interpretations for this. The high GC viruses contain a large number of unique genes and a low number of viruses sequenced for each genus, making it difficult to identify a gene as truncated or fragmented since there is no intact gene in a closely related virus with which to compare. Because of this, several genes may be classified as intact, when they are in fact truncated or fragmented. The second possible reason is that the high GC viruses may be subject to different evolutionary mechanisms, including a lack of gene reduction as compared to viruses with higher AT content and higher numbers of nonintact genes. However, a number of non-high GC viruses also contain low numbers of nonintact genes including LSDV, DPV, SWPV, the leporipoxviruses, and the yatapoxviruses. In addition, CPXV isolates only contain 3–6 nonintact genes in contrast to the other orthopoxviruses that contain between 15 and 33 nonintact genes. This emphasizes that GC content alone does not directly determine the number of nonintact genes present in any single genome, and that other factors must play a role in determining the extent to which gene truncation and fragmentation play a role in the evolution of poxvirus species.
While most viruses show an association between ESMs and microsatellites, this does not appear to be the case for FWPV, RFV, and TANV, (Table 2, Figure 3B). This may be due to different mechanisms being used for the introduction of ESMs compared to other poxviruses, but it may also be due to the accumulation of mutations in microsatellites that originally colocalized with ESMs, but whose sequence has now diverged significantly making it difficult to identify the former microsatellites.
Another interesting observation is the relatively low number of microsatellites in BPSV. BPSV and ORFV are both in the genus Parapoxvirus, and have genomes that are 134,431 and 139,952 bp long, respectively. We were unable to identify any non-intact genes in either ORFV or BPSV (both are high GC viruses). BPSV has a genomic microsatellite content of 12.2%, compared to ORFV, which is 22.8%, and 20.8%–31.7% for all other chordopoxviruses.
DNA viruses, and especially dsDNA viruses, have been shown to have relatively low mutation rates, more similar to eukaryotes than to RNA viruses [25]. Therefore, mechanisms that would support an increase in the number of mutations or that would target mutations to specific regions of the genome may provide an opportunity to adapt to alterations in the environment, selecting for new genome variants that increase the fitness of the virus. Recently, Elde et al. identified just such a mechanism for vaccinia virus that was able to expand the copy number of an immunomodulatory gene when exposed to selective conditions favoring expression of that gene [49]. The genomic regions containing the gene were most likely duplicated by recombination, however the authors also witnessed loss of the duplications when selection pressure was removed. Gene duplication and reduction may be one method of rapidly altering the poxvirus genome. Increased numbers of mutations due to the presence of microsatellites may be another. The two methods are not mutually exclusive, since non-homologous recombination preferentially occurs at repeat regions [33,34], and this may be a mechanism by which the accordion-like expansion and contraction of genes like the instance observed by Elde et al. may occur.

5. Conclusions

This study highlights the role of microsatellites in poxvirus sequence diversity, and how they may affect viral host range and pathogenesis by playing a role in gene variation and inactivation. We found a high incidence of statistically significant co-localization of ESMs with microsatellites. Genes which are no longer impacted by selection pressure to remain functional may be more likely to accumulate mutations at microsatellites, since these short repeats often show higher variability than the surrounding sequences [50,51]. Microsatellites may therefore serve as a major source of genomic variability in chordopoxviruses, contributing to the introduction of ESMs and therefore impacting gene content, virus biology, and evolution. Microsatellite hypervariability may allow the virus to adapt more quickly to environmental changes, such as the phenotypic changes needed when adapting to newly encountered hosts.

Acknowledgments

This work was supported by NIH/NIAID Grant No. U01 AI48706 and NIH/NIAID Contract No. HHSN266200400036C to E.J.L. We thank Robert Curtis Hendrickson, Don Dempsey, and John Osborne for their contributions to various aspects of this work.

Author Contributions

E. L. H. performed or verified poxvirus genome annotations, developed the methods used to analyze the data, performed the analysis, and wrote the manuscript. C. W. developed the approach and applications used to annotate poxvirus genomes, annotated many of the isolates evaluated in this study, and edited the manuscript. E. J. L. supervised the work and edited the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. International Committee on Taxonomy of Viruses. Virus Taxonomy: Classification and Nomenclature of Viruses: Ninth Report of the International Committee on Taxonomy of Viruses; Elsevier Academic Press: London, UK, 2011. [Google Scholar]
  2. Carroll, D.S.; Emerson, G.L.; Li, Y.; Sammons, S.; Olson, V.; Frace, M.; Nakazawa, Y.; Czerny, C.P.; Tryland, M.; Kolodziejek, J.; et al. Chasing jenner's vaccine: Revisiting cowpox virus classification. PLOS ONE 2011, 6, e23086. [Google Scholar]
  3. Zhao, G.; Droit, L.; Tesh, R.B.; Popov, V.L.; Little, N.S.; Upton, C.; Virgin, H.W.; Wang, D. The genome of yoka poxvirus. J. Virol. 2011, 85, 10230–10238. [Google Scholar] [CrossRef] [PubMed]
  4. Hendrickson, R.C.; Wang, C.; Hatcher, E.L.; Lefkowitz, E.J. Orthopoxvirus genome evolution: The role of gene loss. Viruses 2010, 2, 1933–1967. [Google Scholar] [CrossRef] [PubMed]
  5. Schramm, B.; Locker, J.K. Cytoplasmic organization of poxvirus DNA replication. Traffic 2005, 6, 839–846. [Google Scholar] [CrossRef] [PubMed]
  6. Moss, B. Poxvirus DNA replication. Cold Spring Harb. Perspect. Biol. 2013, 5, a010199. [Google Scholar] [CrossRef]
  7. Lefkowitz, E.J.; Wang, C.; Upton, C. Poxviruses: Past, present and future. Virus Res. 2006, 117, 105–118. [Google Scholar] [CrossRef] [PubMed]
  8. McLysaght, A.; Baldi, P.F.; Gaut, B.S. Extensive gene gain associated with adaptive evolution of poxviruses. Proc. Natl. Acad. Sci. USA 2003, 100, 15655–15660. [Google Scholar] [CrossRef] [PubMed]
  9. Essbauer, S.; Pfeffer, M.; Meyer, H. Zoonotic poxviruses. Vet. Microbiol. 2010, 140, 229–236. [Google Scholar] [CrossRef] [PubMed]
  10. Smith, G.L. Genus Orthopoxvirus: Vaccinia virus. In Poxviruses; Mercer, A.A., Weber, O., Schmidt, A., Eds.; Birkhäuser Verlag: Basel, Switzerland, 2007; pp. 1–45. [Google Scholar]
  11. Emerson, G.L.; Li, Y.; Frace, M.A.; Olsen-Rasmussen, M.A.; Khristova, M.L.; Govil, D.; Sammons, S.A.; Regnery, R.L.; Karem, K.L.; Damon, I.K.; et al. The phylogenetics and ecology of the orthopoxviruses endemic to north america. PLOS ONE 2009, 4, e7666. [Google Scholar]
  12. Weli, S.C.; Tryland, M. Avipoxviruses: Infection biology and their use as vaccine vectors. Virol. J. 2011, 8, e49. [Google Scholar] [CrossRef]
  13. Mockett, B.; Binns, M.M.; Boursnell, M.E.; Skinner, M.A. Comparison of the locations of homologous fowlpox and vaccinia virus genes reveals major genome reorganization. J. Gen. Virol. 1992, 73, 2661–2668. [Google Scholar] [CrossRef] [PubMed]
  14. Tulman, E.R.; Afonso, C.L.; Lu, Z.; Zsak, L.; Kutish, G.F.; Rock, D.L. Genome of lumpy skin disease virus. J. Virol. 2001, 75, 7122–7130. [Google Scholar] [CrossRef] [PubMed]
  15. Afonso, C.L.; Tulman, E.R.; Lu, Z.; Zsak, L.; Osorio, F.A.; Balinsky, C.; Kutish, G.F.; Rock, D.L. The genome of swinepox virus. J. Virol. 2002, 76, 783–790. [Google Scholar] [CrossRef] [PubMed]
  16. Bracht, A.J.; Armien, A.G.; Carrillo, C.; O'Hearn, E.S.; Fabian, A.W.; Moran, K.E.; Lu, Z.; Ariyakumer, D.S.; Rasmussen, J.M.; Metwally, S.A. Isolation and characterization of a cervidpoxvirus from a goitered gazelle (gazella subgutturosa) from a zoologic park in minnesota. J. Zoo Wildl. Med. 2013, 44, 589–595. [Google Scholar] [CrossRef] [PubMed]
  17. Darby, A.C.; McInnes, C.J.; Kjaer, K.H.; Wood, A.R.; Hughes, M.; Martensen, P.M.; Radford, A.D.; Hall, N.; Chantrey, J. Novel host-related virulence factors are encoded by squirrelpox virus, the main causative agent of epidemic disease in red squirrels in the UK. PLOS ONE 2014, 9, e96439. [Google Scholar]
  18. MacNeil, A.; Lederman, E.; Reynolds, M.G.; Ragade, N.J.; Talken, R.; Friedman, D.; Hall, W.; Shwe, T.; Li, Y.; Zhao, H.; et al. Diagnosis of bovine-associated parapoxvirus infections in humans: Molecular and epidemiological evidence. Zoonoses Public Health 2010, 57, e161–e164. [Google Scholar]
  19. Afonso, C.L.; Delhon, G.; Tulman, E.R.; Lu, Z.; Zsak, A.; Becerra, V.M.; Zsak, L.; Kutish, G.F.; Rock, D.L. Genome of deerpox virus. J. Virol. 2005, 79, 966–977. [Google Scholar] [CrossRef] [PubMed]
  20. Hanson, D.; Diven, D.G. Molluscum contagiosum. Dermatol. Online J. 2003, 9, e2. [Google Scholar]
  21. Afonso, C.L.; Tulman, E.R.; Delhon, G.; Lu, Z.; Viljoen, G.J.; Wallace, D.B.; Kutish, G.F.; Rock, D.L. Genome of crocodilepox virus. J. Virol. 2006, 80, 4978–4991. [Google Scholar] [CrossRef] [PubMed]
  22. Li, Y.; Carroll, D.S.; Gardner, S.N.; Walsh, M.C.; Vitalis, E.A.; Damon, I.K. On the origin of smallpox: Correlating variola phylogenics with historical smallpox records. Proc. Natl. Acad. Sci. USA 2007, 104, 15787–15792. [Google Scholar] [CrossRef] [PubMed]
  23. Likos, A.M.; Sammons, S.A.; Olson, V.A.; Frace, A.M.; Li, Y.; Olsen-Rasmussen, M.; Davidson, W.; Galloway, R.; Khristova, M.L.; Reynolds, M.G.; et al. A tale of two clades: Monkeypox viruses. J. Gen. Virol. 2005, 86, 2661–2672. [Google Scholar]
  24. Drake, J.W.; Hwang, C.B. On the mutation rate of herpes simplex virus type 1. Genetics 2005, 170, 969–970. [Google Scholar] [CrossRef] [PubMed]
  25. Duffy, S.; Shackelton, L.A.; Holmes, E.C. Rates of evolutionary change in viruses: Patterns and determinants. Nat. Rev. Genet. 2008, 9, 267–276. [Google Scholar] [CrossRef] [PubMed]
  26. Odom, M.R.; Hendrickson, R.C.; Lefkowitz, E.J. Poxvirus protein evolution: Family wide assessment of possible horizontal gene transfer events. Virus Res. 2009, 144, 233–249. [Google Scholar] [CrossRef] [PubMed]
  27. Hatcher, E.L.; Hendrickson, R.C.; Lefkowitz, E.J. Identification of nucleotide-level changes impacting gene content and genome evolution in orthopoxviruses. J. Virol. 2014, 88, 13651–13668. [Google Scholar] [CrossRef] [PubMed]
  28. Schlotterer, C. Evolutionary dynamics of microsatellite DNA. Chromosoma 2000, 109, 365–371. [Google Scholar] [CrossRef] [PubMed]
  29. Behe, M.J. Tracts of adenosine and cytidine residues in the genomes of prokaryotes and eukaryotes. DNA Seq. 1998, 8, 375–383. [Google Scholar] [PubMed]
  30. Dieringer, D.; Schlotterer, C. Two distinct modes of microsatellite mutation processes: Evidence from the complete genomic sequences of nine species. Genome Res. 2003, 13, 2242–2251. [Google Scholar] [CrossRef] [PubMed]
  31. Coenye, T.; Vandamme, P. Characterization of mononucleotide repeats in sequenced prokaryotic genomes. DNA Res. 2005, 12, 221–233. [Google Scholar] [CrossRef] [PubMed]
  32. Ellegren, H. Microsatellites: Simple sequences with complex evolution. Nat. Rev. Genet. 2004, 5, 435–445. [Google Scholar] [CrossRef] [PubMed]
  33. Moxon, E.R.; Rainey, P.B.; Nowak, M.A.; Lenski, R.E. Adaptive evolution of highly mutable loci in pathogenic bacteria. Curr. Biol. 1994, 4, 24–33. [Google Scholar] [CrossRef] [PubMed]
  34. Toth, G.; Gaspari, Z.; Jurka, J. Microsatellites in different eukaryotic genomes: Survey and analysis. Genome Res. 2000, 10, 967–981. [Google Scholar] [CrossRef] [PubMed]
  35. Kaufer, B.B.; Jarosinski, K.W.; Osterrieder, N. Herpesvirus telomeric repeats facilitate genomic integration into host telomeres and mobilization of viral DNA during reactivation. J. Exp. Med. 2011, 208, 605–615. [Google Scholar] [CrossRef] [PubMed]
  36. Baroudy, B.M.; Moss, B. Sequence homologies of diverse length tandem repetitions near ends of vaccinia virus genome suggest unequal crossing over. Nucleic Acids Res. 1982, 10, 5673–5679. [Google Scholar] [CrossRef] [PubMed]
  37. Jain, A.; Mittal, N.; Sharma, P.C. Genome wide survey of microsatellites in ssdna viruses infecting vertebrates. Gene 2014, 552, 209–218. [Google Scholar] [CrossRef] [PubMed]
  38. George, B.; Mashhood Alam, C.; Jain, S.K.; Sharfuddin, C.; Chakraborty, S. Differential distribution and occurrence of simple sequence repeats in diverse geminivirus genomes. Virus Genes 2012, 45, 556–566. [Google Scholar] [CrossRef] [PubMed]
  39. Chen, M.; Tan, Z.; Jiang, J.; Li, M.; Chen, H.; Shen, G.; Yu, R. Similar distribution of simple sequence repeats in diverse completed human immunodeficiency virus type 1 genomes. FEBS Lett. 2009, 583, 2959–2963. [Google Scholar] [CrossRef] [PubMed]
  40. Chen, M.; Tan, Z.; Zeng, G. Microsatellite is an important component of complete hepatitis c virus genomes. Infect. Genet. Evol. 2011, 11, 1646–1654. [Google Scholar] [CrossRef] [PubMed]
  41. Houng, H.S.; Lott, L.; Gong, H.; Kuschner, R.A.; Lynch, J.A.; Metzgar, D. Adenovirus microsatellite reveals dynamics of transmission during a recent epidemic of human adenovirus serotype 14 infection. J. Clin. Microbiol. 2009, 47, 2243–2248. [Google Scholar] [CrossRef] [PubMed]
  42. Mudunuri, S.B.; Nagarajaram, H.A. Imex: Imperfect microsatellite extractor. Bioinformatics 2007, 23, 1181–1187. [Google Scholar] [CrossRef] [PubMed]
  43. Szpara, M.L.; Tafuri, Y.R.; Parsons, L.; Shamim, S.R.; Verstrepen, K.J.; Legendre, M.; Enquist, L.W. A wide extent of inter-strain diversity in virulent and vaccine strains of alphaherpesviruses. PLOS Pathog. 2011, 7, e1002282. [Google Scholar] [CrossRef]
  44. Viral Bioinformatics Resource Center. Available online: www.vbrc.org (accessed on 5 June 2008).
  45. Katoh, K.; Misawa, K.; Kuma, K.; Miyata, T. Mafft: A novel method for rapid multiple sequence alignment based on fast fourier transform. Nucleic Acids Res. 2002, 30, 3059–3066. [Google Scholar] [CrossRef] [PubMed]
  46. Bratke, K.A.; McLysaght, A. Identification of multiple independent horizontal gene transfers into poxviruses using a comparative genomics approach. BMC Evol. Biol. 2008, 8, e67. [Google Scholar] [CrossRef]
  47. Ouyang, Q.; Zhao, X.; Feng, H.; Tian, Y.; Li, D.; Li, M.; Tan, Z. High gc content of simple sequence repeats in herpes simplex virus type 1 genome. Gene 2012, 499, 37–40. [Google Scholar] [CrossRef] [PubMed]
  48. Tulman, E.R.; Afonso, C.L.; Lu, Z.; Zsak, L.; Kutish, G.F.; Rock, D.L. The genome of canarypox virus. J. Virol. 2004, 78, 353–366. [Google Scholar] [CrossRef] [PubMed]
  49. Elde, N.C.; Child, S.J.; Eickbush, M.T.; Kitzman, J.O.; Rogers, K.S.; Shendure, J.; Geballe, A.P.; Malik, H.S. Poxviruses deploy genomic accordions to adapt rapidly against host antiviral defenses. Cell 2012, 150, 831–841. [Google Scholar] [CrossRef] [PubMed]
  50. Martin, D.P.; Biagini, P.; Lefeuvre, P.; Golden, M.; Roumagnac, P.; Varsani, A. Recombination in eukaryotic single stranded DNA viruses. Viruses 2011, 3, 1699–1738. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Krupovic, M.; Koonin, E.V. Evolution of eukaryotic single-stranded DNA viruses of the bidnaviridae family from genes of four other groups of widely different viruses. Sci. Rep. 2014, 4, 5347. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Hatcher, E.L.; Wang, C.; Lefkowitz, E.J. Genome Variability and Gene Content in Chordopoxviruses: Dependence on Microsatellites. Viruses 2015, 7, 2126-2146. https://0-doi-org.brum.beds.ac.uk/10.3390/v7042126

AMA Style

Hatcher EL, Wang C, Lefkowitz EJ. Genome Variability and Gene Content in Chordopoxviruses: Dependence on Microsatellites. Viruses. 2015; 7(4):2126-2146. https://0-doi-org.brum.beds.ac.uk/10.3390/v7042126

Chicago/Turabian Style

Hatcher, Eneida L., Chunlin Wang, and Elliot J. Lefkowitz. 2015. "Genome Variability and Gene Content in Chordopoxviruses: Dependence on Microsatellites" Viruses 7, no. 4: 2126-2146. https://0-doi-org.brum.beds.ac.uk/10.3390/v7042126

Article Metrics

Back to TopTop