Research | Open | Published:
Molecular characterization and phylogenetics of Fennoscandian cowpox virus isolates based on the p4c and atip genes
Virology Journalvolume 11, Article number: 119 (2014)
Cowpox virus (CPXV), a rodent-borne Orthopoxvirus (OPV) that is indigenous to Eurasia can infect humans, cattle, felidae and other animals. Molecular characterization of CPXVs isolated from different geographic locations is important for the understanding of their biology, geographic distribution, classification and evolution. Our aim was to characterize CPXVs isolated from Fennoscandia on the basis of A-type inclusion (ATI) phenotype, restriction fragment length polymorphism (RFLP) profiles of atip gene fragment amplicon, and phylogenetic tree topology in conjunction with the patristic and genetic distances based on full length DNA sequence of the atip and p4c genes.
ATI phenotypes were determined by transmission electron microcopy and RFLP profiles were obtained by restriction enzyme digestion of the atip gene fragment PCR product. A 6.2 kbp region spanning the entire atip and p4c genes of Fennoscandian CPXV isolates was amplified and sequenced. The phylogenetic affinity of Fennoscandian CPXV isolates to OPVs isolated from other geographic regions was determined on the basis of the atip and p4c genes.
Fennoscandian CPXV isolates encoded full length atip and p4c genes. They produce wild type V+ ATI except for CPXV-No-H2. CPXVs were resolved into six and seven species clusters based on the phylogeny of the atip and p4c genes respectively. The CPXVs isolated from Fennoscandia were grouped into three distinct clusters that corresponded to isolates from Norway, Sweden and Finland.
CPXV is a polyphyletic assemblage of six or seven distinct clusters and the current classification in which CPXVs are united as one single species should be re-considered. Our results are of significance to the classification and evolution of OPVs.
The Poxviridae is a family of large DNA viruses that multiplies in the cytoplasm of infected cells. The Poxviridae family members are classified into two sub-families based on the hosts they infect; viruses belonging to Chordopoxvirinae infect vertebrates whereas members of Entomopoxvirinae infect insects. The Chordopoxvirinae consists of at least eight genera and the Orthopoxvirus (OPV) is the most studied genus of the chordopoxviruses. OPVs are broadly divided into “Old World” and “North American” species. Members of the Old World OPVs include Variola virus (VARV), the etiologic agent of smallpox; Vaccinia virus (VACV), the vaccine virus used to eradicate smallpox; Cowpox virus (CPXV), a rodent- borne zoonotic OPV that is indigenous to Eurasia; Monkeypox virus (MPXV), a zoonotic OPV that causes smallpox-like diseases in humans; Ectromelia virus (ECTV), the etiologic agent of mousepox (lab mice); Camelpox virus (CMLV) and Taterapox virus (TATV). OPV species presumed to be endemic to North America include Raccoonpox virus (RCNV), Volepox virus (VPXV), and Skunkpox virus (SKPV)[1, 3].
With the eradication and subsequent cessation of vaccination against smallpox, other OPVs especially MPXV in Central/West Africa, CPXV in Europe, VACV in Brazil and the Indian subcontinent are staging a comeback[4–11]. There is an increasing incidence of CPXV infections in humans, domestic cats, zoo animals and wild life[12, 13]. Although human cowpox virus infection is in general characterized by mild and self-limiting lesions, they can also be fatal[14, 15] especially in immune-compromised individuals. Molecular characterization of CPXVs isolated from different geographic locations is important in understanding their geographic distribution, variability and evolution, as well as monitoring the emergence of atypical CPXV strains with enhanced virulence. Such highly virulent CPXVs may arise due to (i) adaptive mutations in rodent or human host, (ii) recombination with other naturally occurring OPVs, (iii) recombination with genetically engineered OPVs or (iv) accidental escape of a lethal recombinant CPXV from laboratory containment. Recently, a CPXV cluster whose genomes are closer to CMLV, TATV and VARV has been identified. In addition, the molecular characterization of CPXVs from different geographic regions will provide a baseline for assessing the potential for recombination between poxvirus-vectored vaccines and naturally circulating OPVs in regions where the former will be used to vaccinate against human, pets, production animals and wild-life diseases.
CPXVs are genetically heterogeneous[17, 18]. They also show marked differences in phenotypic properties. A recent phylogenetic analysis of the genome of 12 CPXVs shows that CPXV is a polyphyletic assemblage with at least five clades. Although this study is the most comprehensive to date, it included only two isolates from Fennoscandia. The under-representation of the Fennoscandian isolates biased the assignment of the Norwegian human CPXV isolate to the same clade as strains isolated in the United Kingdom. The inclusion of more isolates from Scandinavia would help to clarify the phylogenetic position of Scandinavian CPXVs in relation to other OPVs. In order to study the evolutionary relationship between CPXVs and other OPVs, we sequenced the region spanning the entire p4c and atip open reading frames (ORFs) of CPXV isolated from specific Fennoscandian geographical locations and determined the phylogenetic relationships between these CPXVs and other OPVs on the basis of these sequences. The atip gene encodes the A-type inclusion protein (ATIP)[21, 22] while the p4c encodes a protein that is required but not sufficient[24, 25] for the formation of wild type V+ A-type inclusion (ATI). It has been demonstrated that functionally intact atip, p4c and VACV Copenhagen A27L homologue are required for the formation of wild type V+ ATI. The atip and p4c genes were selected because of the genetic stability being located in the central region of the genome, their high degree of conservation, their role in production of ATI, an inclusion body with a significant role in host to host transmission and viral survival outside the host, as well as history of use as a marker for differentiating OPVs. Both the p4c and atip genes are located in the central part of the genome and it is preferable to perform evolutionary analysis of OPV relatedness or divergence based on genes located in the central genome region, since genes located in the terminal region are highly variable due to adaptive selection and recombination[17, 27, 28]. Previously, amplification of an atip gene fragment with specific primers in conjunction with RFLP generated by restriction enzyme digestion of the atip gene fragment amplicon have been used to differentiate OPV species[29, 30].
In this study, we investigated the sequence diversity and phylogeny of Fennoscandian CPXVs based on the atip and p4c genes. CPXVs isolated from Fennoscandia encoded intact atip and p4c ORFs. They also produced wild type V+ ATI except for CPXV-No-H2. Phylogenetically, CPXVs were resolved into six and seven distinct clusters based on the atip and p4c genes respectively. Fennoscandian CPXVs were segregated into three of the distinct clusters and isolates from one Fennoscandian country have closer phylogenetic relationship to each other than isolates from another country or geographical region.
Fennoscandian CPXV isolates produce wild type V+ A-type inclusions
The production of ATI was one of the first biological properties used to differentiate CPXV from VACV. Among OPVs, only CPXV, ECTV and RCNV are known to produce ATI in infected cells[22, 32, 33]. Three strain specific ATI phenotypes are known to exist; the V+ ATI in which virions are within the ATI matrix, the V− phenotype in which the ATI lacks virions within or on the surface, and the V+/ phenotype in which the ATIs lack internalized virions but are encrusted with virions on their surface[24, 32, 33]. To determine the ATI phenotype produced in cells infected with Fennoscandian CPXVs, we quantified ATI phenotypes in 50 cell sections that were clearly infected. Absolute and relative amounts of ATI phenotypes produced by each virus strain at 36 hours post infection were determined in infected Vero and A549 cells. Our results show that all the Fennoscandian CPXVs except CPXV-No-H2 produced the wild type V+ ATI in both Vero and A549 cells (Table 1, Figure 1). In Vero cells, the V+ATI ranged from 76.5% in CPXV-Swe-H1 to 98.6% in CPXV-FIN/T2000, while in A549 cells, the V+ ATI accounted for 93.9% - 100% of all ATI phenotype produced for each Fennoscandian CPXV strain (Table 1). The reason why more ATIs are produced in A549 cells compared to Vero was not investigated but may be related to the severity of cytopathic effects (CPE) in the respective cell lines. Fennoscandian CPXVs (except CPXV-FIN/T2000) were more cytopathic in Vero than in A549 cells (data not shown). Presumably, the reduced CPE in A549 leaves cells virtually intact allowing virion morphogenesis (including ATI production) to proceed optimally. The V+ ATI is very stable and remains intact even after cell lysis (Figure 1F).
The atip gene of different geographic CPXV isolates displays distinct XbaI restriction enzyme profiles
The XbaI restriction enzyme digestion of PCR products amplified with ATI-2 primer pairs[30, 33] has proved to be robust in identifying and differentiating OPVs. To differentiate Fennoscandian CPXVs, we amplified a fragment of the atip gene with ATI-2 primer pairs (Figure 2A) and digested the obtained PCR products for 2 hours with XbaI (Figure 2B). Samples showing partial digests after two hours incubation with XbaI were further incubated for another two hours to obtain full digest (Figure 2C). Our results show that an approximate 1659 bp amplicon was generated for CPXV-BR (2) and CPXV-FIN/T2000 (9), 1650 bp product for CPXV-No-H1 (3), CPXV-No-F1 (5), CPXV-No-F2 (6), CPXV-Swe-H1 (7), CPXV-Swe-H2 (8), and approximately 1186 bp and 1279 bp products were amplified for CPXV-No-H2 (4) and ECTV-MOS (10) (Figure 2A). Six groups of XbaI restriction patterns were identified (Figure 2B and C). CPXV-BR belonged to group 1 while CPXV-No-H1, CPXV-No-F1, and CPXV-No-F2 belonged to group 2. Group 3 viruses contained the two Swedish isolates and the Finnish isolate belonged to group 4. Group 5 viruses were made up of CPXV-No-H2 and group 6 contained ECTV-MOS. DNA sequencing of the ATI-2 PCR products in conjunction with in silico XbaI digestion of the DNA sequences gave the exact sizes of the amplicons and the XbaI restriction fragments (Table 2). Group 1 (CPXV-BR) produced fragments of 645, 515, 346, 100, and 67 bp); Group 2 (CPXV-No-H1, CPXV-No-F1 and CPXV-No-F2) produced fragments of 645, 543, 346, and 67 bp; Group 3 (CPXV-Swe-H1, CPXV-Swe-H2) generated fragments of 645, 443, 346, 100, and 67 bp) (Table 2). The group 4 viruses (CPXV-FIN/T2000) produced fragments of 645, 345, 299, 144, 100, 72, and 67 bp. The group 5 (CPXV-No-H2) generated fragments of 482, 343, 151, and 151 bp while the group 6 (ECTV-MOS) produced fragments of 575, 343, 151, 151 (Table 2). The restriction pattern of Fennoscandian CPXVs is different from that of the reference strain, CPXV-BR. Thus with the exception of CPXV-No-H2, the restriction pattern of Fennoscandian CPXVs correlates with the geographic origin of the isolates. The CPXV-No-H2 restriction pattern is similar to that of ECTV-MOS. Characterization of CPXV-No-H2 as a naturally occurring recombinant between ECTV and CPXV has been reported elsewhere.
Fennoscandian CPXVs encode intact atip and p4c ORFs
All the Fennoscandian CPXVs (except CPXV-No-H2) produced wild type V+ ATI in infected cells and since functionally intact atip and p4c genes are essential for the formation of V+ ATI[23, 25], we wanted to find out whether these two genes are intact or disrupted in CPXVs isolated from Fennoscandia. We therefore sequenced the region spanning the atip and p4c genes as well as their flanking sequences and identified the atip and p4c ORFs. The Fennoscandian CPXV isolates contained intact atip ORFs encoding a polypeptide of 1258 amino acids (aa) in two Norwegian feline isolates, 1262 aa in the two Swedish isolates and 1280 aa in the Finish isolate [HQ680374-HQ680378]. Compared to the reference strain CPXV-BR, the atip gene of CPXV isolates from Norway and Sweden contained a 72 bp deletion, but this deletion was absent in the Finnish isolate (Additional file1). The 72 bp deletion is also present in the atip gene of other OPV species including strains of CMLV, MPXV, TATV, VACV and VARV (Additional file1). The p4c ORFs of the Swedish isolates encoded a polypeptide of 529 aa, CPXV-No-F1 and CPXV-No-F2 contained a polypeptide of 522 aa and 517 aa respectively, and that of the Finnish isolate has a p4c coding region corresponding to 518 aa [GenBank accession numbers HQ680374-HQ680378]. The P4c protein is highly conserved among strains of OPVs. Variability in the P4c aa sequence of Fennoscandian CPXVs and other OPVs was present only in the C-terminally located polyaspartate tract (Additional file2). The number of consecutive aspartic acid residues in the P4c polypeptide of CPXV-No-F1, CPXV-No-F2, CPXV-Swe-H1, and CPXV-FIN/T2000 was 24, 19, 31, and 20 respectively (Additional file2). In CPXV-Swe-H2 a valine residue was sandwiched between 10 and 20 consecutive aspartic acid residues (Additional file2). Compared to the Fennoscandian CPXVs, the reference strain CPXV-BR has a disrupted p4c ORF due a single nucleotide (Adenine) deletion at position 765 of the alignment (Additional file3). This deletion resulted in a frame shift mutation that introduced a premature stop codon (TAG) at positions 782 – 784 of the alignment (Additional file3) and this in turn truncated the CPXV-BR p4c gene into CPXV-BR 161 and CPXV-BR 159 ORFs (http://poxvirus.org/gene_detail.asp?gene_id=41865, Additional file3). A disrupted p4c gene is responsible for the V− ATI produced in cells infected with CPXV-BR The presence of full length atip and p4c genes in CPXVs isolated from Fennoscandia may account for their production of V+ ATIs in infected cells.
Phylogeny, genetic and patristic distances based on the atip gene
To examine the evolutionary relationships between CPXVs isolated from Fennoscandia with other CPXVs/OPVs, the individual alignments of the atip and p4c nucleotide sequences were subjected to phylogenetic analysis using Maximum - Likelihood (ML), Neighbor – Joining (NJ) and Bayesian – Inference (BI) methods. The atip and p4c sequences were not concatenated because the Kishino-Hasegawa (KH), Shimodaira-Hasegawa (SH) and Approximately Unbiased (AU) paired tests showed that the phylogenetic signals from the two genes were significantly incongruent (P < 0.05) as to warrant separate analysis. CPXV-No-H2 was excluded from evolutionary analysis of the atip gene because a recombination event was detected just downstream the start of the atip gene. In addition to tree topology, CPXVs were also separated into different clusters on the basis that the genetic and patristic distances between different clusters or groups were equal or more than the TATV-CMLV threshold. The TATV-CMLV threshold was chosen because it represents the lowest distance between distinct OPV species and thus can serve as the lowest reference distance for grouping OPVs into the same or different species. The ML tree based on the atip gene showed that CPXVs were divided into six monophyletic clusters (Figure 3). Group 1 is made up of Norwegian isolates, group 2 contains isolates from the United Kingdom, group 3 included Swedish isolates and some German isolates, group 4 and group 5 were made of one French and one German isolates, while group 6 was made up of CPXV isolates from Finland, Russia and Austria (Figure 3). Both the patristic and genetic distances for each of these groups or clusters were equal to or exceeded the TATV-CMLV threshold (Table 3). The TATV-CMLV patristic distance cutoff as well as the TATV-CMLV genetic distance threshold is 0.013 for the atip gene. CPXV groups 1–5 are described as CPXV-like because they are closer to reference CPXV isolate (CPXV-BR) while CPXV group 6 are VACV-like because they are closer to VACV than to any other CPXV including CPXV-BR. Indeed, the patristic and genetic distances between CPXV group 6 (CPXV-FIN/T2000, CPXV_FIN2000_MAN, CPXV_GRI_90, CPXV_AUS1999_867) and VACV (VACV-WR, VACV-Acambis, VACV-lister, VACV-CVA) were 0.026 and 0.022 but the patristic distance between CPXV 6 and any other cluster (CPXV1-5) ranged from 0.038 to 0.061 while the genetic distance varied from 0.032 to 0.049 (Table 3). The NJ tree topology was similar to trees generated by the ML method while the BI trees did not show good resolution like the ML and NJ trees (data not shown). In addition, phylogeny based on the ATIP amino acid sequences yielded tree topologies that were similar to those based on atip nucleotide sequences (data not shown). Overall the atip gene phylogeny showed that CPXVs resolved into two major monophyletic clades (CPXV-like and VACV-like) and these two major clades were further resolved into six monophyletic clusters or groups. The CPXVs isolated from Fennoscandia were resolved into three of the six monophyletic clusters.
Phylogeny, genetic and patristic distances based on the p4c gene
Taken into cognizance the results of the Kishino-Hasegawa (KH), Shimodaira-Hasegawa (SH) and Approximately Unbiased (AU) paired tests, it was expected that the p4c phylogenetic tree topology would be different from that of the atip in some nodes and branches. The p4c gene phylogenetic tree was constructed with ECTV (Figure 4, Additional file4) and CPXV-No-H2 (Figure 5) as outgroup taxa as well as an unrooted tree (without assumptions on common decent) in which both ECTV and CPXV-No-H2 were excluded (Additional file5). The BI tree (Figure 4) and the ML trees (Additional file4, Figure 5, and Additional file5) based on the p4c gene resolved CPXVs into seven clusters or groups. The patristic and genetic distances between any of the seven clusters were in general equal to or exceeded the TATV-CMLV threshold although there were a few exceptions (Table 4). The TATV-CMLV patristic and genetic distance thresholds were 0.011 and 0.012 respectively for the p4c gene (Table 4). The CPXV isolates from Norway and the United Kingdom resolved into CPXV groups 1 and 2 respectively, similar to what was obtained with the tree based on the atip gene. However, unlike the atip gene phylogeny, CPXV 3 contained the Swedish and French isolates while two and three German isolates belonged to CPXV 4 and CPXV 5 (Figure 4, Additional file4, Figure 5 and Additional file5). In addition, the CPXV 6 as obtained with the atip based phylogeny was split into groups 6 and 7 in the p4c based phylogeny. Group 6 contained a single isolate from Austria while group 7 was made up of isolates from Russia and Finland (Figure 4, Additional file4, Figure 5 and Additional file5). The group 7 CPXVs (CPXV-FIN/T2000, CPXV_FIN2000_MAN, and CPXV_GRI_90) was closer to VACV than any other CPXV clusters (Figure 5, Additional file5). Thus the patristic and genetic distances between CPXV 7 and VACV was 0.013 and 0.014, whereas the patristic and genetic distances between other CPXV clusters (groups 1–6) and VACV had a range of 0.023 – 0.029 and 0.021 – 0.029 (Table 4). Unlike the p4c phylogenetic tree rooted with CPXV-No-H2 (Figure 5) and the unrooted tree (Additional file5), the p4c gene phylogeny with ECTV as the outgroup taxa (Figure 4 and Additional file4) appear not to corroborate the observation that CPXV group 7 is closer to VACV than to any other CPXV cluster. Although the BI (Figure 4) and ML (Additional file4) phylogenetic trees have similar topology, the BI tree has very strong bootstrap support in almost all the nodes (posterior probability of 1.0), while the ML tree has low clade support (<70%) in some of the nodes. The NJ tree topology for the p4c gene was similar to that of the BI and ML trees (data not shown). Also, the tree topologies obtained for the P4c amino acids were similar to what was obtained with the p4c nucleotide sequences (data not shown). Overall the phylogeny based on the p4c gene as well as the atip gene has demonstrated the profound genetic heterogeneity among CPXVs, dividing CPXVs into six or seven distinct clusters. The Fennoscandian CPXVs were found to belong to three of the seven distinct clusters.
The objective of this study was to characterize CPXVs isolated from Fennoscandia and compare their biological and genetic characteristics to that of CPXVs isolated from other geographic regions, as well as to other OPV species. The Fennoscandian CPXVs were characterized on the basis of ATI phenotype, RFLP of atip gene fragment amplicon, sequence and phylogenetic analysis based on full length atip and p4c genes. We have demonstrated that CPXVs isolated from Fennoscandia produced wild type V+ ATI (except for CPXV-No-H2) and encode full length atip and p4c genes. The V+ ATI was produced in both Vero and A549 cells, suggesting that it is a strain specific trait. Functionally intact atip and p4c genes have been shown to be essential but not sufficient[24, 33] for the formation of V+ ATI. In addition to functionally intact atip and p4c genes, it has been demonstrated that VACV Copenhagen A27L homologue is required for the formation of V+ ATI. The production of V+ ATI in cells infected with Fennoscandian CPXVs may be due to the presence of functionally intact atip, p4c and A27L sequences. Although these three genes have been shown to be essential for the formation of V+ ATI, it cannot be excluded that other genes are involved in the formation of this phenotype. We are currently investigating whether the three genes (atip, p4c, A27L homologue) are sufficient for the production of wild type V+ ATI or if additional genes are required. The occlusion of virions within ATI may facilitate host to host transmission by protecting the virions from the harsh environment during transmission between hosts. Thus, Fennoscandian CPXVs with the exception of CPXV-No-H2 may have evolved the V+ATI phenotype to increase virus survival capability outside the host and aid host-to-host transmission under extreme climatic conditions in Fennoscandia, especially the seasonal freeze-thaw cycles.
The RFLP profiles generated by digesting PCR products amplified with ATI-2 primer pairs have enabled correct species assignment of 73 OPV isolates already known to belong to different OPV species. Our application of this method to CPXVs isolated in Fennoscandia yielded RFLP profiles that correlate with the geographic region of the isolates. Although the RFLP profiles 2–4 were unique to Fennoscandian CPXVs, they were more or less related to published profiles of other CPXVs. However the RFLP profile 5 for CPXV-No-H2 was similar to that of ECTV and was completely different from any known CPXV profile. This observation was the first clue of a suspected recombination event. We have reported elsewhere that CPXV-No-H2 is a novel recombinant between CPXV and ECTV. Thus, in addition to robust differentiation and classification of OPV strains, RFLP profiles generated by digesting ATI-2 primer amplicons can serve as a first indicator of atypical or recombinant CPXV.
The phylogenetic tree topology in tandem with genetic and patristic distances has been used for robust molecular taxonomy of OPVs. We employed the same method in classifying Fennoscandian CPXVs and other CPXVs/OPVs, with the modification that we chose to use the TATV-CMLV threshold rather than VARV-TATV that was reported by others. We chose the TATV-CMLV threshold because TATV is closer to CMLV than to VARV, thus the distances between TATV and CMLV represents the lowest distance between distinct OPV species, and CPXVs whose genetic/patristic distances equals or exceeds the TATV-CMLV threshold should not be grouped as the same species. The atip gene phylogeny showed that OPVs were divided into two major monophyletic clades that were further subdivided into six clusters. The CPXV-like clade was exclusively made of five different CPXV clusters (CPXV 1–5) that did not contain any other OPV species, while the VACV-like clade contained VACV-like CPXVs (CPXV 6) in addition to MPXV, VACV, TATV and CMLV. These results are in agreement with a recent phylogenetic analysis based on ten concatenated conserved genes. However, it has to be noted that the bootstrap support for the VACV-like clade is low. One possible explanation for the low VACV-clade support is that vaccine or laboratory strains of VACV were used for the construction of the phylogenetic tree. These vaccine/laboratory strains usually have high number of passages in cell cultures or laboratory animal, and these passages may have introduced mutations into the genome of these strains. Presumably, these mutations might have interfered with the phylogenetic signals. To test this hypothesis; we reconstructed the atip gene phylogeny using sequences of naturally isolated VACV strains. The most improved VACV clade support was obtained when only HSPV-MNR76 and VACV-3737 sequences were used. VACV clade support was improved from less than 50% (Figure 3) to 61% (Additional file6) while the CPXV-clade was improved from 73% (Figure 3) to 82% (Additional file6). These results suggest that mutations introduced in the atip gene of vaccine/laboratory strains of VACV following passages in cell culture or laboratory animals may in part account for the low VACV clade support. Alternatively, low clade support may indicate that members of the clade are not phylogenetically related. However this is unlikely in this case as phylogenetic analysis based on multiple genes as well as the entire conserved central region of OPV genome have yielded tree topology consistent with the atip gene phylogeny reported in this study with the exception that bootstrap support for the VACV-like clade was 95% or more. It has been estimated that CPXV-like and VACV-like clades of OPVs diverged from the common ancestor some ten thousand years ago (TYA), a period that corresponds to the last ice age over Fennoscandia. It has been speculated that the wide abundance of various rodents (natural host to CPXVs and other OPVs) may have fuelled OPV divergent evolution. However, full genome OPV phylogeny grouped CPXVs classified as CPXV-like as sister to CMLV, TATV and VARV. This is in disagreement with the atip gene phylogeny reported in this paper and also with phylogeny based on concatenated conserved genes published elsewhere. The reason(s) for this discrepancy is unclear but it may be that phylogeny based on whole genome provides better resolved trees than those obtained with single genes or concatenated multiple genes. Again, while the atip gene phylogeny reported in this study resolved CPXV_NOR_1994_MAN into the same species cluster as other CPXVs isolated from Norway (CPXV group 1) and distinct from isolates from the United Kingdom (CPXV group 2), the genome based phylogenetic tree grouped CPXV_NOR_1994_MAN into the same cluster as isolates from the United Kingdom. There are two likely reasons for this. First, only one isolate from Norway and no isolate from Sweden were used in the whole genome phylogenetic tree reconstruction. This might have under-represented phylogenetic signals for CPXVs isolated from Norway. Secondly, the higher TATV-VARV threshold was used in the whole genome study as opposed to the lower TATV-CMLV threshold that was used in this study. Indeed if the lower TATV-CMLV threshold was used in the whole genome study, CPXV_NOR_1994_MAN would be distinct from isolates from the United Kingdom. Thus, it can be concluded that CPXV_NOR_1994_MAN clustered with other CPXV isolates from Norway and is distinct from CPXVs isolated in the United Kingdom.
The phylogeny based on the p4c gene showed that CPXVs classified as CPXV-like (CPXV group 1–5) formed a major OPV monophyletic clade that included CMLV, TATV and VARV. Thus, as opposed to the atip gene phylogram, the p4c tree shows that CPXVs classified as CPXV-like (CPXV group 1–5) were sister to TATV, CMLV and VARV. This is in accordance with the results of whole genome phylogeny, but in contrast to the phylogram obtained from concatenating conserved genes located at the central region of the genome. In addition, the p4c genetic and patristic distance measures showed that CPXV group 7 was closer to VACV than to any other OPV. Surprisingly, the p4c gene phylogeny with ECTV as outgroup taxa (Figure 4, Additional file4) does not corroborate the findings obtained from patristic and genetic distance measures. We suspected that the reason for this is that the ECTV sequence is not an uncontroversial outgroup taxa as it is too close to VACV. We tested this hypothesis by reconstructing an unrooted p4c gene phylogenetic tree and our result confirmed that the ECTV p4c sequence was very close to the VACV (data not shown). When ECTV was excluded and an unrooted tree was reconstructed, it clearly showed that CPXV group 7 formed the same monophyletic clade as VACV (Additional file5). Previously, we have shown that the p4c gene of CPXV-No-H2 is diverged from the homologues in all OPV species and probably represent an ancestral sequence. A reconstruction of the p4c gene phylogeny with CPXV-No-H2 sequence as out group taxa demonstrated that CPXV group 7 formed the same clade with VACV (Figure 5). This is in agreement with the whole genome phylogeny. The evidence that CPXV-No-H2 has ECTV atip gene while its p4c gene was shown to be distant from other CPXVs (Figure 5) may raise the suspicion that it is not a CPXV strain. In addition to clinical history we have shown that CPXV-No-H2 is a CPXV strain based (i) presence of two copies of cytokine response modifier B (crmB) gene, (ii) sequence and phylogenetic analysis based on multiple genes including crmB, Chinese hamster ovary host range (CHOhr) gene, and the haemagglutinin (HA) gene[18, 33]. Our previous HA phylogenetic tree construction did not include all the CPXVs used in this study. Therefore we reconstructed the HA phylogenetic tree including all the CPXVs used in this study as well as representatives of both “Old World” and “North America” OPV species. Our result clearly showed that CPXV-No-H2 has phylogenetic affinity with other CPXVs belonging to CPXV-like clade and it is phylogenetically distinct from other “Old World” and “North American” species of OPVs (Additional file7). Unlike the p4c BI phylogenetic tree which has strong posterior probability of 1.0 (100%) in most nodes including the node that grouped CPXV 1–5 with CMLV, TATV and VARV, the p4c ML tree has low bootstrap support in some of nodes. The low clade support in some nodes in the p4c gene ML tree may be due to under-estimation of clade support by the ML algorithm.
The conflicting phylogenetic signals between the atip and p4c genes as shown by both the paired tests and some of the branching and nodes in the tree topologies were rather surprising since both genes are located in the same central part of the genome (the p4c gene is just upstream of the atip) and are required for the same function, that is, the formation of ATI phenotypes. The conflicting phylogenetic signals between atip and p4c genes may indicate that either or both genes are involved in functions other than the formation of ATI, and thus may have different rates of evolution. In spite of this conflict in phylogenetic signals between the atip and p4c, both phylograms have demonstrated that (i) CPXVs are genetically heterogeneous and can be subdivided into six or seven species clusters, (ii) CPXVs isolated from Fennoscandia belong to three of the distinct clusters and that each of the clusters contain isolates from one specific geographic region or country, (iii) CPXV_NOR_1994_MAN belonged to the same cluster as other Norwegian isolates and was distinct from isolates from the United Kingdom, (iv) VACV-like CPXVs (CPXV_GRI_90, CPXV-FIN/T2000, CPXV_FIN2000_MAN, CPXV_AUS1999_867) are closer to VACV than any of the other 14 CPXV isolates (CPXV group 1–5), (v) TATV, CMLV and VARV belong to the same monophyletic cluster, and TATV is closer to CMLV than to VARV, (vi) HSPV is grouped together with VACV isolates, and (vii) MPXVs were resolved into Congo basin and West Africa sub-clusters. Other investigators have demonstrated similar results with respect to OPV phylogeny and evolution[17, 35–37].
The fact that CPXV isolates from Norway and Sweden were classified as CPXV-like while isolates from Finland in conjunction with the isolate from Russia were VACV-like may be an indication that CPXV-like and VACV-like CPXVs have two distinct evolutionary histories in geographically or otherwise separated rodent lineages that re-colonized Fennoscandia after the post glacial retraction 12 – 8 TYA. A recent retrospective study on OPV molecular evolution showed that CPXV-like viruses may have separated from the common ancestor approximately 10 TYA and commenced individual evolution. This event may have triggered the divergent evolution of OPVs that lead to the emergence of VACV-like CPXVs and other OPV species. Intriguingly, the estimated 10 TYA age for OPV divergent evolution corresponds well to the time of post glacial retraction (12 – 8 TYA), thus supporting the post glacial retraction hypothesis for OPV evolution. Moreover, post-glacial re-colonization hypothesis has been used to explain the distribution of different Puumala Hantavirus (PUU) genotypes in Norway and Sweden[38, 39], and it is interesting that both PUU and CPXV have bank voles as reservoir species[38–40]. Detailed studies including a large number of CPXVs isolated from different geographic regions in Fennoscandia as well as their genome sequences will be required in order to determine whether or not (or to what degree) post-glacial re-colonization could explain the genetic diversity among CPXVs.
CPXVs isolated from Fennoscandia encode full length atip and p4c genes and produced wild type V+ ATI except for CPXV-No-H2. CPXVs were genetically heterogeneous and were resolved into six or seven distinct clusters. The Fennoscandian CPXV belonged to three of the distinct clusters and isolates from the same country or geographic region belonged to the same cluster. Our results show that VACV-like CPXVs and CPXV-like CPXVs may have distinct evolutionary histories, presumably in rodent host species. These results are of relevance to the evolution and classification of CPXVs in particular and OPVs in general.
Cells and viruses
Vero (ATCC CCL-81) and A549 (ATCC CCL-185) were obtained from the American Type Culture Collection (ATCC) (Rockville, MD, USA). Both cell lines were grown under conditions suggested by ATCC. The origin of Fennoscandian CPXV isolates used in this study has been described elsewhere[41–43]. The Swedish CPXV isolates designated Swe-H1 and Swe-H2 were isolated from human clinical cases in the district of Skåne in 1990. The Norwegian isolates No-H1 and No-FI were isolated in 1994 from a human and a cat in the district of Bergen. The Norwegian isolate, No-F2, was isolated from a cat from the same Bergen district in 1999, and another Norwegian human isolate No-H2 was isolated in 2001 from a male teenager in Nordland, Norway. The Finnish isolate, FIN/T2000, was isolated from a four year old girl in 2000. The reference strain CPXV Brighton Red (CPXV-BR) was included for comparison/control and ECTV strain Moscow (ECTV-MOS) was included for comparison of RFLP profiles. ECTV-MOS was purchased from ATCC. All isolates except No-F1 and No-F2 have been passaged in cell cultures prior to arrival in our laboratory. The isolates have been passaged once or twice prior to arrival in our laboratory. Stocks of virus isolates were prepared from infected Vero cells after three rounds of plaque purification/amplification and titer of stocks were determined by plaque assay.
Transmission electron microscopy
Confluent monolayers of Vero and A549 cells in six well tissue culture plates (NUNC, Sweden) were infected with Fennoscandian CPXVs and CPXV-BR at a multiplicity of infection (m.o.i.) of 5 plaque forming units (pfu) per cell. Viruses were adsorbed to monolayers for one hour at 4°C. Infected cells were washed twice with phosphate buffered saline (PBS) and incubated in a medium containing 2.5% fetal bovine serum (FBS) at 37°C in a 5% CO2 atmosphere. At 36 hours post infection (hpi), infected cells were washed with fresh medium and fixed in MacDowell’s solution, pH7.4, for one hour at room temperature. Fixed cells were processed for transmission electron microscopy as described elsewhere. Quantification of ATI phenotypes by electron microscopy in sections of infected cells was as we previously described. Collection of images from thin sections was done in JEOL electron microscopy operating at an accelerating voltage of 100 K.
Viral DNA was extracted from the lowest cell culture passage of viral stocks using DNA Minikit (QIAGEN GmbH Strasse 1, Hilden, Germany). Eight primer pairs at different molar concentrations were used to generate amplicons spanning the entire atip and p4c ORF as well as their flanking sequences. All the eight primer pairs; ATIN-1, ATIN-3, ATIN-4, ATI-2, ATI-5, p4c-1, p4c-2, p4c-3 have been published previously. The PCR reaction for each gene fragment was performed with 2 ng DNA in a total volume of 50 μl PCR mix. The PCR mixes contained appropriate concentration of primer pairs, 200 μM of each dNTP, 1 x reaction buffer (10 mM Tris–HCl, pH 8.3, 50 mM KCL, 1.5 mM MgCl2) and 1.25U Amplitaq DNA polymerase (Applied Biosystems, Foster City, USA). The PCR cycling conditions for all the primer pairs have been described elsewhere[24, 33]. The PCR products were resolved in 1% Seakem agarose gels (Cambrex Bioscience, Rockland, ME, USA) in 1 x TAE buffer. Pictures of ethidium bromide stained DNA fragments were photographed with Gel Doc 2000 (Biorad).
Restriction enzyme digestion of ATI-2 PCR product
The ATI-2 gene fragment of Fennoscandian CPXV, CPXV-BR, and ECTV-MOS was amplified from virus DNA as previously described[30, 33]. The ATI-2 PCR product was digested with 20 units of XbaI and incubated at 37°C for 2 hours. Amplicons that were not fully digested within the two hour period were subjected to additional two hours of digestion. The restriction digests were separated using 2.0% Metaphore agarose (Medprobe) gels. The RFLP profiles were photographed with Gel Doc 2000 (Biorad).
The PCR products of the eight overlapping fragments for each Fennoscandian CPXVs were purified using GFX PCR DNA and Gel Band Purification Kit (GE Health, Uppsala, Sweden), following the manufacturers instruction. Cycle sequencing reactions were performed using Big Dye 3.1 Sequencing Kit (Applied Biosystems, Foster City, CA, USA). Each purified DNA fragment was sequenced in both orientations and at least two independent DNA cycling reactions were performed for each fragment. The cycle sequencing extension products were electrophoresed using ABI Prism™ 377 DNA automatic sequencer (Applied Biosystems).
Contig assembly, sequence analysis and recombination detection
Raw sequence data was edited in TREV 1.9 (http://staden.sourceforge.net). Forward and reverse DNA sequences for fragments 1–8 for each of the Fennoscandian CPXVs were assembled into a single contig using Chromas Pro 1.5 software (http://technelysium.com.au/?page_id=27). The atip and p4c ORFs were identified from respective contigs using ORF finder (http://www.ncbi.nlm.nih.gov/gorf/gorf.html). The atip and p4c sequences (nucleotide and amino acids) of Fennoscandian CPXV isolates and those of other OPVs (including 12 other CPXVs) (Table 5) were retrieved from the GenBank and aligned with ClustalW version 1.8. Positions with gaps were removed from the atip and p4c alignments (nucleotide and amino-acids). The DNA and amino acid identity matrices were calculated from aligned sequences using Bioedit version 18.104.22.168 (http://www.mbio.ncsu.edu/BioEdit/page2.html). For recombination detection, the approximately 6200 nt region containing the atip and p4c ORFs of Fennoscandian CPXVs, other CPXVs and ECTV-MOS were aligned with ClustalW 1.8, and subjected to recombination analysis. Recombination detection was carried out with SimPlot version 3.5.1 and Recombination detection program 3 (RPD3). The sequences reported in this paper were deposited under GenBank accession numbers [HQ680372, HQ680374-HQ680378].
Phylogenetic analysis, patristic and genetic distance matrices
The best fit nucleotide substitution model was selected based on JmodelTest while the best fit model for protein evolution was selected using ProTest. Kishino-Hasegawa (KH), Shimodaira and Hasegawa (SH), and Approximately Ubiased (AU) paired tests implemented in Treefinder version march 2011 (http://www.treefinder.de/) were used to examine the possibility of concatenating the atip and p4c aligned sequences for phylogenetic analysis. A significant p-value (P < 0.05) rejects the null hypothesis (all the tree topologies are good explanations of the data) and precludes combining the atip and p4c molecular datasets. Phylogenetic trees were constructed for atip a nd p4c alignments using ML and NJ methods as implemented in MEGA version 5.05, BI as implemented in MrBayes version 3.1.2. ECTV-MOS and ECTV strain Naval (ECTV-NAV) were included as outgroup taxa in all phylogenetic analysis. CPXV-No-H2 p4c gene was included as an outgroup for the p4c phylogeny. Bayesian analysis was conducted with the following settings: nst = 6, rates = invgamma, mcmc (2 simultaenous runs), samplefreq = 1000, nchains = 4 (3 hot, 1 cold), 10 million generations, and burnin = 2500. A 50% majority rule consensus tree was generated and posterior probabilities of 95% or more was considered statistically significant for clade support. ML gene trees were constructed using GTR + G substitution model while the Maximum Composite Likelihood model was used for the construction of NJ gene trees. The robustness of ML and NJ trees was evaluated by bootstrap analysis of 1000 replicates. The pairwise genetic distance without reference to a tree was generated for the atip and p4c genes using the Maximum Composite Likelihood method while the patristic distance (tree branch lengths between taxa) between the isolates was extracted from the ML and BI trees using the program Patristic version 1.0 (http://www.bioinformatics.org/patristic/). For both genetic and patristic distance methods, the distances were averaged across taxa to produce a value for each node. Since TATV and CMLV represent the two closest but distinct OPV species, the distance (genetic and patristic) between these two species represents a threshold value. In addition to tree topology, CPXVs were further separated into different groups if their genetic/patristic distance is equal or exceeds the TATV-CMLV threshold.
Moss B: Poxviridae: The Viruses and Their Replication. In Fields Virology. 4th edition. Edited by: Roizman BHP, Straus S, Martin M, DE G, Lamb R, Knipe D. Philadelphia: Lippincott Williams and Wilkins; 2001:2849-2883.
Mercer AA, Schmidt A, Weber O: Poxviruses. In Advances in Infectious Diseases. Basel: Birkhäuser Verlag; 2007:113-116. 355–373
Emerson GL, Li Y, Frace MA, Olsen-Rasmussen MA, Khristova ML, Govil D, Sammons SA, Regnery RL, Karem KL, Damon IK, Carroll DS: The phylogenetics and ecology of the orthopoxviruses endemic to North America. PLoS One 2009, 4: e7666.
Shah S: New threat from poxviruses. Sci Am 2013, 308: 66-71.
Damaso CR, Esposito JJ, Condit RC, Moussatche N: An emergent poxvirus from humans and cattle in Rio de Janeiro State: cantagalo virus may derive from Brazilian smallpox vaccine. Virology 2000, 277: 439-449.
Quixabeira-Santos JC, Medaglia ML, Pescador CA, Damaso CR: Animal movement and establishment of vaccinia virus Cantagalo strain in Amazon biome, Brazil. Emerg Infect Dis 2011, 17: 726-729.
Rimoin AW, Mulembakani PM, Johnston SC, Lloyd Smith JO, Kisalu NK, Kinkela TL, Blumberg S, Thomassen HA, Pike BL, Fair JN, Wolfe ND, Shongo RL, Graham BS, Formenty P, Okitolonda E, Hensley LE, Meyer H, Wright LL, Muyembe JJ: Major increase in human monkeypox incidence 30 years after smallpox vaccination campaigns cease in the democratic republic of Congo. Proc Natl Acad Sci U S A 2010, 107: 16262-16267.
Kroon EG, Mota BE, Abrahao JS, da Fonseca FG, de Souza Trindade G: Zoonotic Brazilian Vaccinia virus: from field to therapy. Antiviral Res 2011, 92: 150-163.
Assis FL, Almeida GM, Oliveira DB, Franco-Luiz AP, Campos RK, Guedes MI, Fonseca FG, Trindade GS, Drumond BP, Kroon EG, Abrahao JS: Characterization of a new vaccinia virus isolate reveals the C23L gene as a putative genetic marker for autochthonous group 1 Brazilian vaccinia virus. PLoS One 2012, 7: e50413.
Singh RK, Hosamani M, Balamurugan V, Bhanuprakash V, Rasool TJ, Yadav MP: Buffalopox: an emerging and re-emerging zoonosis. Anim Health Res Rev 2007, 8: 105-114.
Goyal T, Varshney A, Bakshi SK, Barua S, Bera BC, Singh RK: Buffalo pox outbreak with atypical features: a word of caution and need for early intervention! Int J Dermatol 2013, 52: 1224-1230.
Campe H, Zimmermann P, Glos K, Bayer M, Bergemann H, Dreweck C, Graf P, Weber BK, Meyer H, Buttner M, Busch U, Sing A: Cowpox virus transmission from pet rats to humans, Germany. Emerg Infect Dis 2009, 15: 777-780.
Essbauer S, Pfeffer M, Meyer H: Zoonotic poxviruses. Vet Microbiol 2010, 140: 229-236.
Eis-Hubinger AM, Gerritzen A, Schneweis KE, Pfeiff B, Pullmann H, Mayr A, Czerny CP: Fatal cowpox-like virus infection transmitted by cat. Lancet 1990, 336: 880.
Czerny CP, Waldmann R, Scheubeck T: Identification of three distinct antigenic sites in parapoxviruses. Arch Virol 1997, 142: 807-821.
Dabrowski PW, Radonic A, Kurth A, Nitsche A: Genome-wide comparison of cowpox viruses reveals a new clade related to Variola virus. PLoS One 2013, 8: e79953.
Gubser C, Hue S, Kellam P, Smith GL: Poxvirus genomes: a phylogenetic analysis. J Gen Virol 2004, 85: 105-117.
Hansen H, Okeke MI, Nilssen O, Traavik T: Comparison and phylogenetic analysis of cowpox viruses isolated from cats and humans in Fennoscandia. Arch Virol 2009, 154: 1293-1302.
Duraffour S, Mertens B, Meyer H, van den Oord JJ, Mitera T, Matthys P, Snoeck R, Andrei G: Emergence of cowpox: study of the virulence of clinical strains and evaluation of antivirals. PLoS One 2013, 8: e55808.
Carroll DS, Emerson GL, Li Y, Sammons S, Olson V, Frace M, Nakazawa Y, Czerny CP, Tryland M, Kolodziejek J, Nowotny N, Olsen-Rasmussen M, Khristova M, Govil D, Karem K, Damon IK, Meyer H: Chasing Jenner’s vaccine: revisiting cowpox virus classification. PLoS One 2011, 6: e23086.
Patel DD, Pickup DJ: Messenger RNAs of a strongly-expressed late gene of cowpox virus contain 5′-terminal poly(A) sequences. EMBO J 1987, 6: 3787-3794.
Patel DD, Pickup DJ, Joklik WK: Isolation of cowpox virus A-type inclusions and characterization of their major protein component. Virology 1986, 149: 174-189.
McKelvey TA, Andrews SC, Miller SE, Ray CA, Pickup DJ: Identification of the orthopoxvirus p4c gene, which encodes a structural protein that directs intracellular mature virus particles into A-type inclusions. J Virol 2002, 76: 11216-11225.
Okeke MI, Adekoya OA, Moens U, Tryland M, Traavik T, Nilssen O: Comparative sequence analysis of A-type inclusion (ATI) and P4c proteins of orthopoxviruses that produce typical and atypical ATI phenotypes. Virus Genes 2009, 39: 200-209.
Howard AR, Weisberg AS, Moss B: Congregation of orthopoxvirus virions in cytoplasmic A-type inclusions is mediated by interactions of a bridging protein (A26p) with a matrix protein (ATIp) and a virion membrane-associated protein (A27p). J Virol 2010, 84: 7592-7602.
Smith GL, Murphy BJ, Law M: Vaccinia virus motility. Annu Rev Microbiol 2003, 57: 323-342.
Babkin IV, Babkina IN: A retrospective study of the orthopoxvirus molecular evolution. Infect Genet Evol 2012, 12: 1597-1604.
Coulson D, Upton C: Characterization of indels in poxvirus genomes. Virus Genes 2011, 42: 171-177.
Meyer H, Pfeffer M, Rziha HJ: Sequence alterations within and downstream of the A-type inclusion protein genes allow differentiation of Orthopoxvirus species by polymerase chain reaction. J Gen Virol 1994,75(Pt 8):1975-1981.
Meyer H, Ropp SL, Esposito JJ: Gene for A-type inclusion body protein is useful for a polymerase chain reaction assay to differentiate orthopoxviruses. J Virol Methods 1997, 64: 217-221.
Downie AW: A study of the lesions produced experimentally by cowpox virus. J Pathol Bacteriol 1939, 48: 361-379.
Shida H, Tanabe K, Matsumoto S: Mechanism of virus occlusion into A-type inclusion during poxvirus infection. Virology 1977, 76: 217-233.
Okeke MI, Hansen H, Traavik T: A naturally occurring cowpox virus with an ectromelia virus A-type inclusion protein gene displays atypical A-type inclusions. Infect Genet Evol 2012, 12: 160-168.
Meyer H, Damon IK, Esposito JJ: Orthopoxvirus diagnostics. Methods Mol Biol 2004, 269: 119-134.
Li Y, Carroll DS, Gardner SN, Walsh MC, Vitalis EA, Damon IK: On the origin of smallpox: correlating variola phylogenics with historical smallpox records. Proc Natl Acad Sci U S A 2007, 104: 15787-15792.
Tulman ER, Delhon G, Afonso CL, Lu Z, Zsak L, Sandybaev NT, Kerembekova UZ, Zaitsev VL, Kutish GF, Rock DL: Genome of horsepox virus. J Virol 2006, 80: 9244-9258.
Likos AM, Sammons SA, Olson VA, Frace AM, Li Y, Olsen-Rasmussen M, Davidson W, Galloway R, Khristova ML, Reynolds MG, Zhao H, Carroll DS, Gurns A, Formenty P, Esposito JJ, Regnery RL, Damon IK: A tale of two clades: monkeypox viruses. J Gen Virol 2005, 86: 2661-2672.
Horling J, Lundkvist A, Jaarola M, Plyusnin A, Tegelstrom H, Persson K, Lehvaslaiho H, Hornfeldt B, Vaheri A, Niklasson B: Distribution and genetic heterogeneity of Puumala virus in Sweden. J Gen Virol 1996,77(Pt 10):2555-2562.
Lundkvist A, Wiger D, Horling J, Sjolander KB, Plyusnina A, Mehl R, Vaheri A, Plyusnin A: Isolation and characterization of Puumala hantavirus from Norway: evidence for a distinct phylogenetic sublineage. J Gen Virol 1998,79(Pt 11):2603-2614.
Sandvik T, Tryland M, Hansen H, Mehl R, Moens U, Olsvik O, Traavik T: Naturally occurring orthopoxviruses: potential for recombination with vaccine vectors. J Clin Microbiol 1998, 36: 2542-2547.
Pelkonen PM, Tarvainen K, Hynninen A, Kallio ER, Henttonen K, Palva A, Vaheri A, Vapalahti O: Cowpox with severe generalized eruption, Finland. Emerg Infect Dis 2003, 9: 1458-1461.
Tryland M, Myrmel H, Holtet L, Haukenes G, Traavik T: Clinical cowpox cases in Norway. Scand J Infect Dis 1998, 30: 301-303.
Cronqvist J, Ekdahl K, Kjartansdottir A: Cowpox - a cat disease in man. Larkartidningen 2006, 88: 2605-2606.
Katz E, Wolffe E, Moss B: Identification of second-site mutations that enhance release and spread of vaccinia virus. J Virol 2002, 76: 11637-11644.
Thompson JD, Higgins DG, Gibson TJ: Improved sensitivity of profile searches through the use of sequence weights and gap excision. Comput Appl Biosci 1994, 10: 19-29.
Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, Ingersoll R, Sheppard HW, Ray SC: Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol 1999, 73: 152-160.
Martin D, Rybicki E: RDP: detection of recombination amongst aligned sequences. Bioinformatics 2000, 16: 562-563.
Posada D: jModelTest: phylogenetic model averaging. Mol Biol Evol 2008, 25: 1253-1256.
Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics 2005, 21: 2104-2105.
Jobb G, von Haeseler A, Strimmer K: TREEFINDER: a powerful graphical analysis environment for molecular phylogenetics. BMC Evol Biol 2004, 4: 18.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol 2011, 28: 2731-2739.
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 2003, 19: 1572-1574.
Huelsenbeck J, Rannala B: Frequentist properties of Bayesian posterior probabilities of phylogenetic trees under simple and complex substitution models. Syst Biol 2004, 53: 904-913.
Fourment M, Gibbs MJ: PATRISTIC: a program for calculating patristic distances and graphically comparing the components of genetic change. BMC Evol Biol 2006, 6: 1.
This study was funded by GenØk-Centre for Biosafety Tromsø Norway. We thank Randi Olsen and Helga-Marie Bye of the Department of Electron Microscopy, UiT The Arctic University of Norway for their technical assistance. We also thank Olli Vapalahti of Hartman Institute Finland for providing the Finnish cowpox virus isolate.
The authors declare that they have no competing interests.
MIO designed the study, carried out the experiments and wrote the draft of the manuscript. ASO participated in the design of the study, interpreted data and revised the manuscript. ØN contributed to the conception and design of the study, interpreted the data and revised the manuscript. UM interpreted the data, revised the manuscript and participated in the design of the study. MT contributed to the study design, interpreted data and revised the manuscript. TB interpreted the data and revised the manuscript. TT conceived the study, contributed to its design, interpreted data and revised the manuscript. All authors read and approved the final manuscript.