Identification of multiple independent horizontal gene transfers into poxviruses using a comparative genomics approach
© Bratke and McLysaght; licensee BioMed Central Ltd. 2008
Received: 10 October 2007
Accepted: 27 February 2008
Published: 27 February 2008
Poxviruses are important pathogens of humans, livestock and wild animals. These large dsDNA viruses have a set of core orthologs whose gene order is extremely well conserved throughout poxvirus genera. They also contain many genes with sequence and functional similarity to host genes which were probably acquired by horizontal gene transfer.
Although phylogenetic trees can indicate the occurrence of horizontal gene transfer and even uncover multiple events, their use may be hampered by uncertainties in both the topology and the rooting of the tree. We propose to use synteny conservation around the horizontally transferred gene (HTgene) to distinguish between single and multiple events.
Here we devise a method that incorporates comparative genomic information into the investigation of horizontal gene transfer, and we apply this method to poxvirus genomes. We examined the synteny conservation around twenty four pox genes that we identified, or which were reported in the literature, as candidate HTgenes. We found support for multiple independent transfers into poxviruses for five HTgenes. Three of these genes are known to be important for the survival of the virus in or out of the host cell and one of them increases susceptibility to some antiviral drugs.
In related genomes conserved synteny information can provide convincing evidence for multiple independent horizontal gene transfer events even in the absence of a robust phylogenetic tree for the HTgene.
Horizontal gene transfer (HGT, also: lateral gene transfer) is an important process in viral evolution. HGT describes an event that differs from the regular 'vertical' transfer of genes from parent to offspring: for example, host gene capture in viruses, where pieces of host DNA are incorporated into the newly formed virion. HGT is thought to play a major role in evolution, as there has apparently been frequent transfer of protein domains or entire genes within and between all divisions of life with important implications for molecular phylogeny studies [1–3]. Especially in pathogens these horizontally transferred genes (HTgenes) are of enormous interest as they often increase virulence by encoding homologs to the host's immune system and interfering with it [4–6]. Multiple independent horizontal gene transfer events are thus particularly significant; they highlight the advantage that is conferred on the virus genome by a gene that has been gained and maintained several times.
Poxviruses are large double-stranded DNA viruses (up to 360 kb) that infect a wide range of hosts from insects to mammals. They are classified into entomopox (EPV, insect-infecting) and chordopox viruses (ChPV, vertebrate-infecting), with ChPV sub-divided into at least nine genera. Several members of the Poxviridae are associated with human disease, including variola (smallpox), Molluscum contagiosum and monkeypox, which has been classifed as an emerging virus [7, 8] and was responsible for an outbreak of monkeypox in the US in 2003. The most well known poxviruses are smallpox, an orthopox virus (OPV) responsible for devastating human pandemics until it was declared eradicated by the WHO in 1980, and cowpox, used in the development of the first vaccine by Edward Jenner in 1796. Over 50 genomes have been fully sequenced and many attempts have been made to reconstruct a robust phylogeny of poxviruses [9–12].
Poxviruses remain in the cytoplasm of the host cell and to a large extent provide their own replication machinery. Many of these genes are not part of a viral core set and it is possible that they were acquired from an ancient host cell by a poxvirus ancestor . The fact that poxviruses do not enter the nucleus poses a practical problem when it comes to the mechanism of host gene capture: how HGT from host to virus takes place is not yet known, but as poxvirus genes lack introns, HGT may proceed through reverse transcription of host mRNAs and integration of the cDNA into the virus genome. . It has recently been suggested that certain poxviruses may act as vectors for the transfer of transposable elements between animal phyla .
There are several approaches to detecting HGT, all of which are best used in conjunction with each other . The first evidence usually comes from sequence similarity to a gene in a distantly related organism or in the viral host. Phylogenetic gene trees which include sequences from both viruses and cellular organisms can also provide an indication of host gene capture where virus genes group closely with homologs from the host or organisms related to the host, i.e. where gene history differs from species history . Atypical base composition of a particular gene, for example a GC content that differs significantly from that of most other genes in the genome, provides another piece of evidence that HGT may have recently taken place and this approach has been used to detect host-derived regions in poxviruses  and other small genomes . This signal, however, is expected to fade over time and therefore this approach is only useful for detecting recent HGT.
These difficulties are compounded when one wishes to test whether there have been multiple independent HGT events of the same gene into related viruses. This may be indicated by a polyphyletic relationship in the phylogenetic tree, but such a tree would require strong support (e.g., high bootstrap values) for several branches in order to be reliable – frequently this support is not there. A mixed pattern of presence and absence of a particular gene in a group of closely related organisms is another indication for horizontal transfer of the gene in question, but an ancestral origin with loss in some genomes can also explain this observation.
We use comparative gene order information to distinguish between a shared history of the HTgene in the virus genomes and several independent transfer events. Where a given HTgene is in an equivalent location in terms of its neighboring genes in all of the genomes only a single transfer event is supported. A HTgene that is located in a different genomic region in related viruses strengthens the case for several transfer events rather than a common origin. This is an especially strong tool in poxviruses where the gene order is highly conserved throughout the syntenic core of the genome [9, 20].
We systematically searched for evidence of HGT between poxviruses and other organisms using strict search criteria and identified nine candidate HTgenes, some of which had been previously reported. We examined both the tree topology and the genomic context for evidence of the number of HGT events. We also applied the comparative synteny method to several candidate horizontally transferred genes reported in the literature [11, 21].
Results and Discussion
Poxvirus gene families and species phylogeny
Core Conserved Protein Families
Transcription factor VLTF
Viral core protein p4b
Subunit of VETF (82 kDa)
Virion core protein p4a
Myristylated protein/membrane protein
DNA helicase/transcriptional elongation factor
Holliday junction endonuclease/palmitylprotein
Subunit of VITF-3 (45 kDa)
RNA polymerase subunit rpol 32
ATPase/DNA packaging/virion assembly
Subunit of VETF (70 kDa)
ATPase/nucleoside triphosphate phosphohydrolase/transcription termination factor
Subunit of mRNA capping enzyme/transcription termination factor
Poly A polymerase subunit
Sulfhydryl oxidase/redox protein/S-S bond formation/virus assembly
Membrane protein/S-S bond formation pathway protein
Structural virion core protein
Myristylated protein/membrane protein
Virion core cysteine protease
Poly A polymerase subunite/cap methyltransferase
RNA polymerase subunit
ImV envelope/membrane protein
Complete Poxvirus Genomes
Amsacta moorei entomopoxvirus
Melanoplus sanguinipes entomopoxvirus
Canarypox virus strain ATCC VR-111
Fowlpox virus isolate HP-438 (Munich)
Molluscum contagiosum virus subtype 1
Bovine popular stomatitis virus strain BV-AR02
Orf virus strain OV-SA00
Orf virus strain OV-IA82
Orf virus strain NZ2
Ectromelia virus strain Moscow
Cowpox virus strain Brighton Red
Monkeypox virus strain Zaire-96-I-16
Monkeypox virus strain Congo_2003_358
Monkeypox virus strain Zaire_1979-005
Monkeypox virus strain USA_2003_044
Monkeypox virus strain USA_2003_039
Monkeypox virus isolate Sierra Leone
Monkeypox virus strain Liberia_1970_184
Monkeypox virus strain WRAIR7-61
Monkeypox virus strain COP-58
Camelpox virus M-96 (Kazakhstan)
Camelpox virus CMS
Variolar minor virus strain Garcia-1966
Variola major virus strain Bangladesh-1975
Variola virus strain India-1967
Cowpox virus strain GRI-90
Vaccinia virus strain Acambis 3000 Mod. Virus Ankara
Vaccinia virus strain 3737
Vaccinia virus strain Acambis clone 2000
Vaccinia virus strain Acambis clone 3
Vaccinia virus Western Reserve
Vaccinia virus strain Tian Tan
Vaccinia virus Copenhagen
Vaccinia virus strain Lister
Vaccinia virus strain LC16m8
Vaccinia virus strain LC16mO
Yaba monkey tumor virus
Yaba-like disease virus
Myxoma virus strain Lausanne
Rabbit fibroma virus
Deerpox virus W-848-83
Deerpox virus W-1170-84
Sheeppox virus 10700-99 strain TU-V02127
Sheeppox virus strain
Sheeppox virus strain NISKHI
Goatpox virus strain Pellor
Goatpox virus strain G20-LKV 19
Lumpy skin disease virus isolate Neethling vaccine LW1959
Lumpy skin disease virus isolate Neethling Warmbaths LW
Lumpy skin disease virus isolate Neethling 2490
Within-genus relationships for genera with many representatives were inferred by the inclusion of as many gene sequences as possible. For parapox, capripox and orthopox, we aligned all ortholog families conserved in the genus of interest (i.e. genes connected by yellow lines in Additional File 1) and inferred trees for those genera using ML and NJ (see Methods). NJ and ML results were congruent and the bootstrapped ML results are shown in Fig. 1.
Evidence for horizontal gene transfer from sequence similarity
Gene Families With Potential HGT
Number of HGT events
Family Size ¶
Candidate HTgenes identified by PSI-BLAST search
entomo, avi, para, ortho, clade II
avi, para, yata, capri
Ribonucleotide reductase, large subunit
Ribonucleotide reductase, small subunit
croc, avi, or-tho, clade II
avi, ortho, clade II
entomo, avi, ortho, clade II
Candidate HTgenes from the literature
NlpC/P60 peptidase (G6R)
ortho, clade II
avi, mollusci, ortho, yata, deer, sui
yata, sui, deer
entomo, avi, lepori
α-soluble NSF attachment protein (SNAP)
DNA-directed RNA poly-merase, subunit A
entomo, ortho, lepori, sui, deer, capri
Uracil DNA glycolase
Interferon gamma receptor
ortho, lepori, sui, deer, capri
Evidence from gene order comparisons
Poxviruses have extremely well conserved gene order, particularly of the core families (Additional File 1 and ). We estimated the overall level of synteny conservation throughout poxvirus genomes by considering the conservation of neighboring genes between orthologs in pox. For these purposes we examined only genes and neighbors which had been grouped into families without paralogs and which were present in at least two complete genomes. We considered that if a given gene always has the same two immediate neighbors in every genome in which it is found, then it has perfect neighbor conservation. By contrast, if the neighbors are always different, then there is zero conservation. We measured this using the following statistic: (M-N)/(M-2); where M is the maximum possible number of unique neighbours (i.e., different neighbors in every case) and N is the observed number of neighbours. This will always return a value between 0 and 1. We found that 69% of ortholog families had a neighbor conservation value of at least 0.9 and that 89% had at least 0.7 neighbor conservation.
We made three kinds of illustrations of the genomic context of HTgenes, distinguished by the type of gene families that the HTgene is plotted alongside. In the first kind the HTgene location is plotted with respect to the "core families" only. There are only 33 of these genes, so the distances within the diagram may be quite large, however because these are shared by all pox genomes, this view permits a low-resolution comparison across all genera. The second type of diagram plots the HTgene with respect to ortholog families, which are much greater in number. This increases the "resolution" of the plot and still avoids the confounding effect of paralogous genes. Finally we plot the location of the HTgene with respect to any of the gene families and examine it with the caveat that paralogous genes may be misleading.
We applied this analysis to the nine HTgenes uncovered by our PSI-BLAST search and compared the results to the scenarios supported by ML and NJ tree of the pox gene family and the non-pox homologs. We also applied the comparative gene arrangement analysis to candidate HTgenes from the literature (Table 3).
In the case of dUTPase and Uracil DNA glycolase there was insufficient gene order conservation to allow comparisons – this reflects the fact that many of the transfers were into the poxvirus chromosome ends where the genome appears to be more labile. Also, the serpin gene family contains many paralogs within poxviruses which made comparison of gene order between viruses intractable. However, we found that in most cases gene order comparisons between the poxviruses can distinguish between single or multiple origins of a HTgene in poxviruses, even in the absence of a robust phylogenetic tree. The gene location comparison provided supporting evidence for at least two independent origins of: Glutathione Peroxidase; Interleukin-10; Ribonucleotide Reductase, small subunit; Thymidine Kinase; Deoxyribo-Dipyrimidine Photolyase (discussed below). In some cases, such as Interleukin-10, the synteny information is suggestive of more HGT events than were apparent from the phylogenetic tree alone. In other cases, such as Ribonucleotide Reductase, small subunit, the converse was true, and the synteny data indicated that there were fewer events than did the (poorly supported) phylogenetic tree.
Interleukin-10 (IL-10) is a multifunctional cytokine that suppresses inflammation, antiviral response and other host functions. It inhibits activation and maturation of dendritic cells and thus reduces their antigen-presenting ability and inhibits cytokine synthesis via T- and NK- cells. Homologs of IL-10 have been found in many different viruses [24–26]. Viral IL-10 has been shown to play an important role in delaying the development of acquired immunity to orf virus in humans  and to vaccinia in mice . The origin of some viral homologs to IL-10 has been investigated and it is suggested that orf virus (a parapoxvirus) captured the gene from sheep or goat [25, 29].
The NJ tree for IL-10 genes (based on their aa sequences) found in pox and other organisms indicates that there are two origins of this gene in poxviruses because the pox genes group into two clades, i.e., they are polyphyletic (Additional File 2; and ). However, the low bootstrap support for the tree topology means that this tree alone does not give strong support for multiple horizontal transfer events. In fact, the ML tree has a different topology indicating monophyly and also has low bootstraps (Additional File 2). In this case, the lack of clarity in the phylogenetic tree possibly stems from a lack of information within the alignment (Additional File 3), as indicated by the low bootstrap support.
Including more genes in the synteny comparison distinguishes the genomic locations of the yatapox and capripox IL-10 genes as each is present in a region of genome that is well-conserved between these pox genera (thus permitting comparison) yet different for the two sets of viruses (Fig. 3B), indicating distinct origins for the yatapox and capripox IL-10 genes. This same distinction was found by Hughes based on phylogenetic analysis . The IL-10 genomic locations in parapox and avipox are not sufficiently conserved to permit comparison. However, we note that if one or both of these originated in the same event as the yatapox or capripox IL-10 then that also requires many loss events in other pox genera as well as a loss and subsequent re-acquisition in either yatapox or capripox. The gene order with respect to all neighboring gene families (figure not shown) supports the distinction of yatapox and capripox IL-10 but is not otherwise informative of the history of IL-10 due to lack of significant conservation, so cannot distinguish these possibilities. We therefore infer at least two independent HGT events of the IL-10 gene from vertebrates to poxviruses.
Thymidine kinase (TK) is present in many poxviruses [30–32]. It phosphorylates a wide range of substrates and, surprisingly, has been found to increase susceptibility to antiviral drugs by selectively phosphorylating some thymidine analog drugs to a greater degree than the host TK1 gene and thus increasing the efficacy of the antiviral drug . Separate origins for avipox and orthopox TK have been suggested based on phylogenetic evidence [29, 34] and this is supported here by the different genomic environments in avipox and other chordopox (Additional File 4; and ).
The gene order is completely scrambled in entomopox, however, so these data do not distinguish between shared origins for entomopox and either the avipox or other chordopox TK genes and distinct origins (i.e., the data do not distinguish between two or three transfer events). The former scenario would require that the TK gene was present in the common ancestor of all pox and subsequently lost in either avipox or the other chordopox only to be regained independently. If the entomopox and avipox TK gene are descended from the same HGT event, then this scenario would require three events (HGT into pox ancestor, loss from chordopox excluding avipox, HGT into another chordopox lineage). In a simplistic sense this is equally parsimonious with three HGT events. A more biologically realistic analysis should include consideration of any selection to retain, lose or acquire this gene; properties that we do not currently know. The fact that this gene has been acquired and retained independently at least twice indicates that multiple HGT events are possible. Loss and reacquisition in the chordopox lineage may indicate significant changes in the environment have changed the nature of the selection acting on this gene.
In avipox, the TK gene and the adjacent HT-motif protein gene (HT) are in the location occupied by ribonucleotide reductase, large subunit (RRL) in orthopox and suipox. This, along with reported 15 bp direct repeats flanking TK and HT led to the suggestion that these genes transposed within avipox to this location in avipox and replaced RRL . This suggestion is at odds with the experimental strategy employed here because if intra-genomic translocations are a common occurrence it would mean that completely distinct genomic neighborhoods are not sufficient to indicate distinct origins. However, there is no real support for this transposition hypothesis from the many data now available: (i) The phylogenetic evidence supports an independent origin of avipox TK; (ii) There is no evidence that the RRL gene was ever present in avipox because it is not found outside orthopox and suipox (whose most recent common ancestor is not shared with avipox); (iii) Furthermore, we cannot identify the 15 bp repeat referred to in the 1992 paper  but we did identify a 12 bp direct repeat (starting nucleotide positions 87698 and 88647 in AF198100). One of the repeat copies is inside the open reading frame of the I5L gene, which adds the constraint that if it is the result of an insertion-site duplication then it must not have disrupted the function of this gene. The relationship of the entomopox TK gene to the other pox TK genes remains inconclusive, with a third independent transfer being one possibility.
We conclude from the gene order data that there were at least two independent horizontal transfers of TK into poxviruses, and possibly a third (into entomopox). The independent transfer into the same locus of the TK gene in avipox and the RRL gene into the ancestor of orthopox and suipox suggests the intriguing possibility that this genomic location has properties that make it particularly receptive to gene transfers.
Ribonucleotide Reductase, small subunit
Ribonucleotide reductase (RR) functions as a heterodimeric tetramer and leads to the reduction of all four ribonucleotides to deoxyribonucleotides in DNA biosynthesis. The small subunit (RRS) was sequenced in vaccinia and found to be highly similar to eukaryotic RSS genes, including the mouse RSS, suggesting host gene capture . RRS is present in most chordopox genera except parapox, molluscipox and some avipox. The avipox RRS gene is present in a different synteny block than the RRS gene in other chordopox (Additional File 5). These data support two independent transfers of RRS into pox: one into the avipox lineage, and one into the ancestor of orthopox and clade II.
The viral glutathione peroxidase (GP) protects virus and infected cells against oxidative damage resulting from interaction with the immune system . It is not present in orthopox and was first identified in molluscum contagiosum virus . GP is found in avipox and molluscipox and separate transfer events have previously been suggested . The synteny diagram with ortholog families clearly shows that the gene is in an equivalent location in all avipox genomes (Additional File 6), but is in a distinct location in molluscipox and thus supports independent transfers of this gene into these two lineages.
Photolyases are enzymes which catalyse the repair of UV-induced DNA damage. In poxviruses, photolyases are implicated in the survival of the virus outside host cells in between infections . Deoxyribodipyrimidine Photolyase (DP) is found in entomopox, avipox and leporipox. It was identified as a candidate HTgene through sequence similarity with genes in insect and vertebrate genomes . The avipox gene lies in a region of conserved synteny with other ChPV, but the DP gene is not at the syntenous location in leporipox (Additional File 7). This supports independent origins for the avipox and leporipox DP genes. It is not possible to compare the location of the entomopox genes with each other or with ChPV due to insufficient synteny conservation. We conclude that there were at least two independent transfers of DP into poxviruses.
Here we have gone beyond the comparison of the sequences of homologous genes and have introduced comparative synteny data into the analysis of horizontal gene transfer. These data are completely independent of phylogenetic tree and sequence composition data. For the 24 genes with evidence of horizontal transfer analysed here, all but four had sufficiently conserved synteny between poxviral genomes to permit comparison. In five cases, the presence of the HTgene in a well-conserved, yet distinct, region in different poxvirus genomes supported a conclusion of multiple independent transfers of a homologous gene into Poxviridae. In the other cases, this study confirmed that a single event gave rise to the HTgene in poxviruses.
This strategy is based on the assumption that there is no bias in the insertion site of a HTgene with respect to the neighboring genes in the receiving (pox) genome. The only known bias in terms of the location of new genes in poxviruses is that they tend to be located towards the ends of chromosomes  and this is not expected to influence the neighbors for this analysis. More importantly, we interpret the conservation of synteny across poxvirus genomes as evidence against genomic rearrangement events in that region. This means that under these circumstances we can exclude the possibility that the new genomic location is caused by genome rearrangement within the pox lineage rather than independent transfer into the lineage.
A major strength of this method is that it does not rely on phylogenetic trees (and thus multiple sequence alignments) in order to interpret the number of transfer events. In previous studies many HTgenes were refractory to phylogenetic analysis due to difficulties rooting the trees or insignificant support for critical branches (e.g.  in pox, and ). Phylogenetic trees are frequently plagued by problems which abolish their usefulness for analysis of horizontal gene transfer . Because the conserved synteny data are independent of the phylogenetic tree for the HTgene these results could be combined.
The HTgenes with multiple origins identified here have very interesting biological properties. The presence of viral TK increases susceptibilty to antiviral drugs, a surprising finding . The particular benefit to the virus, if any, of having this gene is not clear. Three other HTgenes, IL-10, GP and DP are known to improve the survival of the virus against the host immune system (IL-10) and environmental damage (GP and DP). In the case of these genes it is easy to speculate that they are particularly advantageous to the virus, and the fact that they have been successfully transferred into poxvirus genomes on more than one occasion lends support to this hypothesis.
We downloaded all available poxvirus proteins from GenBank in June 2006 and predicted proteins for unannotated complete genomes using the EMBOSS getorf application  (total 17,185 sequences). These proteins were compared in an all-against-all BLAST search using default parameters . We used The Markov Cluster Algorithm (MCL)  to group proteins into families with the Centering parameter set to 1.2 and the Inflation parameter set to 1.2. These parameters were chosen after experimentation in order to maximize the family size. Alignments of all of these families were performed using MUSCLE and are provided as supplementary material (Additional File 8).
Species Phylogeny Construction
Thirty-three families were found in single-copy in every completely sequenced poxvirus genome ("core families"; Table 1). We generated alignments of the protein sequences of each of these genes using MUSCLE  and concatenated the alignments. Based on this concatenated sequence alignment, phylogenetic trees of 54 completely sequenced poxviruses (Table 2) were reconstructed by different methods: NJ (PHYLIP package  and ClustalW ), MP (PHYLIP package) and ML (PhyML package ). The best models for ML analysis were estimated using PhyML  by experimenting with a range of settings for several parameters (gamma parameter, substitution model, number of substitution rate categories) and choosing the ones that gave the highest likelihood. The JTT and TN93 models were used for protein and nucleotide alignments, respectively. To allow for bootstrapping, a less computationally intensive, representative phylogeny including at most two members of each genus was constructed. Single gene trees were also constructed using PhyML from the 33 conserved families, and their congruency with the concatenated gene tree was noted. To infer within-genus relationships, ortholog families shared by all members of a genus were aligned at the protein level and the gaps were copied into the nucleotide sequence to generate nucleotide (nt) alignments. The nt alignments were concatenated and the tree of the genus was inferred using PhyML with the TN93 model. For comparison we also constructed the trees using NJ (with both the PHYLIP package and ClustalW).
Detection of non-poxviral homologs
We used PSI-BLAST  to search all 17,815 poxvirus proteins against the GenBank non-redundant (nr) database (May 2007; 4.9 million sequences). The e-value threshold for inclusion in the PSSM was e-20, with a maximum of five iterations. Only alignments of 70% identity or higher and at least 80% length of query sequence were kept as homologs. We selected these parameters after some experimentation with different values in order to retrieve only strong candidates for HGT and to keep detection of false positives to a minimum. We excluded hits to synthetic constructs (e.g., artificial vectors).
The protein sequences of members of the nine gene families with PSI-BLAST hits outside poxviruses were aligned with their non-poxviral homologs using MUSCLE . NJ trees were drawn with ClustalW  and bootstrapped 1000 times. ML trees were drawn with PhyML using model JTT, proportion of invariable sites = 0, number of substitution rate categories = 6 and estimating the gamma distribution parameter.
Three kinds of comparative genomic plots were produced. Only genes that have been classified into families are considered. (1) A low resolution plot illustrates the location of the gene of interest in each genome with respect to "core families". The actual distance between genes on this plot may be quite large. (2) A medium resolution compares the location of the HTgene with respect to "ortholog families". (3) The maximum resolution plot includes neighboring genes from any families.
In all of these diagrams when the gene of interest is present in a given genome, the genes plotted along side are immediate neighbors (of the given type) with no intervening genes. That is, in the ortholog families plot, the genes plotted are all the nearest-neighbor genes with identified orthologs in any genome.
When the gene of interest is absent from a given genome, we allowed some intervening genes and indicate this with the thickness of the line (as per figure legend). We chose ranges of 1–2 and 3–6 intervening genes to indicate different levels of conservation after experimenting with different values as we considered these to be the most clear and informative.
Horizontal gene transfer
Markov Cluster Algorithm
Ribonucleotide Reductase large subunit
Ribonucleotide Reductase small subunit
Horizontally transferred gene
Operational taxonomic unit
We would like to thank DG Knowles and D Kelleher for valuable comments and technical help. This work is supported by an Irish Research Council for Science and Engineering Technology (IRCSET) Postgraduate Fellowship to KAB and the Science Foundation Ireland (SFI) President of Ireland Young Researcher Award to AMcL.
- Doolittle WF: Lateral genomics. Trends Cell Biol. 1999, 9 (12): M5-8. 10.1016/S0962-8924(99)01664-5.View ArticlePubMedGoogle Scholar
- Brown JR: Ancient horizontal gene transfer. Nat Rev Genet. 2003, 4 (2): 121-32. 10.1038/nrg1000.View ArticlePubMedGoogle Scholar
- Choi IG, Kim SH: Global extent of horizontal gene transfer. Proc Natl Acad Sci USA. 2007, 104 (11): 4489-94. 10.1073/pnas.0611557104.PubMed CentralView ArticlePubMedGoogle Scholar
- Barry M, McFadden G: Virus encoded cytokines and cytokine receptors. Parasitology. 1997, 115 (): S89-100. 10.1017/S0031182097001820.View ArticlePubMedGoogle Scholar
- Lawrence JG, Ochman H: Molecular archaeology of the Escherichia coli genome. Proc Natl Acad Sci USA. 1998, 95 (16): 9413-7. 10.1073/pnas.95.16.9413.PubMed CentralView ArticlePubMedGoogle Scholar
- Ochman H, Lawrence JG, Groisman EA: Lateral gene transfer and the nature of bacterial innovation. Nature. 2000, 405 (6784): 299-304. 10.1038/35012500.View ArticlePubMedGoogle Scholar
- Mahalingam S, Damon IK, Lidbury BA: 25 years since the eradication of smallpox: why poxvirus research is still relevant. Trends Immunol. 2004, 25 (12): 636-9. 10.1016/j.it.2004.10.002.View ArticlePubMedGoogle Scholar
- Su JR: Emerging viral infections. Clin Lab Med. 2004, 24 (3): 773-95. 10.1016/j.cll.2004.05.002. viiView ArticlePubMedGoogle Scholar
- McLysaght A, Baldi PF, Gaut BS: Extensive gene gain associated with adaptive evolution of poxviruses. Proc Natl Acad Sci USA. 2003, 100 (26): 15655-60. 10.1073/pnas.2136653100.PubMed CentralView ArticlePubMedGoogle Scholar
- Gubser C, Hue S, Kellam P, Smith GL: Poxvirus genomes: a phylogenetic analysis. J Gen Virol. 2004, 85 (Pt 1): 105-17. 10.1099/vir.0.19565-0.View ArticlePubMedGoogle Scholar
- Hughes AL, Friedman R: Poxvirus genome evolution by gene gain and loss. Mol Phylogenet Evol. 2005, 35: 186-95. 10.1016/j.ympev.2004.12.008.View ArticlePubMedGoogle Scholar
- Xing K, Deng R, Wang J, Feng J, Huang M, Wang X: Genome-based phylogeny of poxvirus. Intervirology. 2006, 49 (4): 207-14. 10.1159/000090790.View ArticlePubMedGoogle Scholar
- Iyer LM, Balaji S, Koonin EV, Aravind L: Evolutionary genomics of nucleo-cytoplasmic large DNA viruses. Virus Res. 2006, 117: 156-84. 10.1016/j.virusres.2006.01.009.View ArticlePubMedGoogle Scholar
- Sekiguchi J, Shuman S: Ligation of RNA-containing duplexes by vaccinia DNA ligase. Biochemistry. 1997, 36 (29): 9073-9. 10.1021/bi970705m.View ArticlePubMedGoogle Scholar
- Piskurek O, Okada N: Poxviruses as possible vectors for horizontal transfer of retroposons from reptiles to mammals. Proc Natl Acad Sci USA. 2007, 104 (29): 12046-51. 10.1073/pnas.0700531104. Epub 2007 Jul 10.PubMed CentralView ArticlePubMedGoogle Scholar
- Lawrence JG, Ochman H: Reconciling the many faces of lateral gene transfer. Trends Microbiol. 2002, 10: 1-4. 10.1016/S0966-842X(01)02282-X.View ArticlePubMedGoogle Scholar
- Shackelton LA, Holmes EC: The evolution of large DNA viruses: combining genomic information of viruses and their hosts. Trends Microbiol. 2004, 12 (10): 458-65. 10.1016/j.tim.2004.08.005.View ArticlePubMedGoogle Scholar
- Da Silva M, Upton C: Host-derived pathogenicity islands in poxviruses. Virol J. 2005, 2 (30): 30-10.1186/1743-422X-2-30.PubMed CentralView ArticlePubMedGoogle Scholar
- Tsirigos A, Rigoutsos I: A sensitive, support-vector-machine method for the detection of horizontal gene transfers in viral, archaeal and bacterial genomes. Nucleic Acids Res. 2005, 33 (12): 3699-707. 10.1093/nar/gki660.PubMed CentralView ArticlePubMedGoogle Scholar
- Mackett M, Archard LC: Conservation and variation in Orthopoxvirus genome structure. J Gen Virol. 1979, 45 (3): 683-701.View ArticlePubMedGoogle Scholar
- Anantharaman V, Aravind L: Evolutionary history, structural features and biochemical diversity of the NlpC/P60 superfamily of enzymes. Genome Biol. 2003, 4 (2): R11-10.1186/gb-2003-4-2-r11.PubMed CentralView ArticlePubMedGoogle Scholar
- Hertig C, Coupar BE, Gould AR, Boyle DB: Field and vaccine strains of fowlpox virus carry integrated sequences from the avian retrovirus, reticuloendotheliosis virus. Virology. 1997, 235 (2): 367-76. 10.1006/viro.1997.8691.View ArticlePubMedGoogle Scholar
- Hayakawa T, Xu J, Hukuhara T: Cloning and sequencing of the gene for an enhancing factor from Pseudaletia separata entomopoxvirus. Gene. 1996, 177 (1–2): 269-70. 10.1016/0378-1119(96)00297-1.View ArticlePubMedGoogle Scholar
- Vieira P, de Waal-Malefyt R, Dang MN, Johnson KE, Kastelein R, Fiorentino DF, deVries JE, Roncarolo MG, Mosmann TR, Moore KW: Isolation and expression of human cytokine synthesis inhibitory factor cDNA clones: homology to Epstein-Barr virus open reading frame BCRFI. Proc Natl Acad Sci USA. 1991, 88 (4): 1172-6. 10.1073/pnas.88.4.1172.PubMed CentralView ArticlePubMedGoogle Scholar
- Fleming SB, McCaughan CA, Andrews AE, Nash AD, Mercer AA: A homolog of interleukin-10 is encoded by the poxvirus orf virus. J Virol. 1997, 71 (6): 4857-61.PubMed CentralPubMedGoogle Scholar
- Kotenko SV, Saccani S, Izotova LS, Mirochnitchenko OV, Pestka S: Human cytomegalovirus harbors its own unique IL-10 homolog (cmvIL-10). Proc Natl Acad Sci USA. 2000, 97 (4): 1695-700. 10.1073/pnas.97.4.1695.PubMed CentralView ArticlePubMedGoogle Scholar
- Chan A, Baird M, Mercer AA, Fleming SB: Maturation and function of human dendritic cells are inhibited by orf virus-encoded interleukin-10. J Gen Virol. 2006, 87 (Pt 11): 3177-81. 10.1099/vir.0.82238-0.View ArticlePubMedGoogle Scholar
- Kurilla MG, Swaminathan S, Welsh RM, Kieff E, Brutkiewicz RR: Effects of virally expressed interleukin-10 on vaccinia virus infection in mice. J Virol. 1993, 67 (12): 7623-8.PubMed CentralPubMedGoogle Scholar
- Hughes AL: Origin and evolution of viral interleukin-10 and other DNA virus genes with vertebrate homologues. J Mol Evol. 2002, 54: 90-101. 10.1007/s00239-001-0021-1.View ArticlePubMedGoogle Scholar
- Upton C, McFadden G: Identification and nucleotide sequence of the thymidine kinase gene of Shope fibroma virus. J Virol. 1986, 60 (3): 920-7.PubMed CentralPubMedGoogle Scholar
- Boyle DB, Coupar BE, Gibbs AJ, Seigman LJ, Both GW: Fowlpox virus thymidine kinase: nucleotide sequence and relationships to other thymidine kinases. Virology. 1987, 156 (2): 355-65. 10.1016/0042-6822(87)90415-6.View ArticlePubMedGoogle Scholar
- Amano H, Ueda Y, Miyamura T: Identification and characterization of the thymidine kinase gene of Yaba virus. J Gen Virol. 1995, 76 (Pt 5): 1109-15.View ArticlePubMedGoogle Scholar
- Prichard MN, Keith KA, Johnson MP, Harden EA, McBrayer A, Luo M, Qiu S, Chattopadhyay D, Fan X, Torrence PF, Kern ER: Selective phosphorylation of antiviral drugs by vaccinia virus thymidine kinase. Antimicrob Agents Chemother. 2007, 51 (5): 1795-803. 10.1128/AAC.01447-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Koonin EV, Senkevich TG: Evolution of thymidine and thymidylate kinases: the possibility of independent capture of TK genes by different groups of viruses. Virus Genes. 1992, 6 (2): 187-96. 10.1007/BF01703067.View ArticlePubMedGoogle Scholar
- Binns MM, Tomley FM, Campbell J, Boursnell ME: Comparison of a conserved region in fowlpox virus and vaccinia virus genomes and the translocation of the fowlpox virus thymidine kinase gene. J Gen Virol. 1988, 69 (Pt 6): 1275-83.View ArticlePubMedGoogle Scholar
- Binns MM, Boursnell ME, Skinner MA: Gene translocations in poxviruses: the fowlpox virus thymidine kinase gene is flanked by 15 bp direct repeats and occupies the locus which in vaccinia virus is occupied by the ribonucleotide reductase large subunit gene. Virus Res. 1992, 24 (2): 161-72. 10.1016/0168-1702(92)90004-S.View ArticlePubMedGoogle Scholar
- Slabaugh M, Roseman N, Davis R, Mathews C: Vaccinia virus-encoded ribonucleotide reductase: sequence conservation of the gene for the small subunit and its amplification in hydroxyurea-resistant mutants. J Virol. 1988, 62 (2): 519-27.PubMed CentralPubMedGoogle Scholar
- Moss B, Shisler JL, Xiang Y, Senkevich TG: Immune-defense molecules of molluscum contagiosum virus, a human poxvirus. Trends Microbiol. 2000, 8 (10): 473-7. 10.1016/S0966-842X(00)01838-2.View ArticlePubMedGoogle Scholar
- Senkevich TG, Bugert JJ, Sisler JR, Koonin EV, Darai G, Moss B: Genome sequence of a human tumorigenic poxvirus: prediction of specific host response-evasion genes. Science. 1996, 273 (5276): 813-6. 10.1126/science.273.5276.813.View ArticlePubMedGoogle Scholar
- Srinivasan V, Schnitzlein WM, Tripathy DN: Fowlpox virus encodes a novel DNA repair enzyme, CPD-photolyase, that restores infectivity of UV light-damaged virus. J Virol. 2001, 75 (4): 1681-8. 10.1128/JVI.75.4.1681-1688.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang J, Mullapudi N, Lancto CA, Scott M, Abrahamsen MS, Kissinger JC: Phylogenomic evidence supports past endosymbiosis, intracellular and horizontal gene transfer in Cryptosporidium parvum. Genome Biol. 2004, 5 (11): R88-10.1186/gb-2004-5-11-r88.PubMed CentralView ArticlePubMedGoogle Scholar
- Podell S, Gaasterland T: DarkHorse: a method for genome-wide prediction of horizontal gene transfer. Genome Biol. 2007, 8 (2): R16-10.1186/gb-2007-8-2-r16.PubMed CentralView ArticlePubMedGoogle Scholar
- Rice P, Longden I, Bleasby A: EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet. 2000, 16 (6): 276-7. 10.1016/S0168-9525(00)02024-2.View ArticlePubMedGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-10.View ArticlePubMedGoogle Scholar
- van Dongen S: Graph Clustering by Flow Simulation. PhD thesis. 2000, University of UtrechtGoogle Scholar
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32 (5): 1792-7. 10.1093/nar/gkh340.PubMed CentralView ArticlePubMedGoogle Scholar
- Felsenstein J: PHYLIP (Phylogeny Inference Package) version 3.6. 2005Google Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22 (22): 4673-80. 10.1093/nar/22.22.4673.PubMed CentralView ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52 (5): 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- Baldo AM, McClure MA: Evolution and horizontal transfer of dUTPase-encoding genes in viruses and their hosts. J Virol. 1999, 73 (9): 7710-21.PubMed CentralPubMedGoogle Scholar
- McGeoch DJ: Protein sequence comparisons show that the 'pseudoproteases' encoded by poxviruses and certain retroviruses belong to the deoxyuridine triphosphatase family. Nucleic Acids Res. 1990, 18 (14): 4105-10. 10.1093/nar/18.14.4105.PubMed CentralView ArticlePubMedGoogle Scholar
- Upton C, McFadden G: DNA sequence homology between the terminal inverted repeats of Shope fibroma virus and an endogenous cellular plasmid species. Mol Cell Biol. 1986, 6: 265-76.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.