Evolution of endogenous retroviruses in the Suidae: evidence for different viral subpopulations in African and Eurasian host species
© Nascimento et al; licensee BioMed Central Ltd. 2011
Received: 23 November 2010
Accepted: 24 May 2011
Published: 24 May 2011
Porcine endogenous retroviruses (PERVs) represent remnants of an exogenous form that have become integrated in the domestic pig (Sus scrofa) genome. Although they are usually inactive, the capacity of γ1 ERVs to infect human cells in vitro has raised concerns about xenotransplantation because the viruses could cross the species barrier to humans. Here we have analyzed the evolution of γ1 ERVs in ten species of Suidae (suids, pigs and hogs) from Eurasia and Africa using DNA sequences for their coding domains (gag, pro/pol and env genes). For comparison with γ1 PERVs, we have also analysed γ2 ERVs which in domestic pigs are known to be inactive and do not pose a risk to xenotransplantation.
Phylogenetic analysis using Bayesian inference showed that γ1 and γ2 ERVs have distinctive evolutionary histories. Firstly, two different viral lineages of γ1 ERVs were found and a coevolutionary analysis demonstrated that they correspond broadly to their host phylogeny, one of Eurasian and another of African species, and show no evidence of horizontal transmission. γ2 ERVs, however, show a bush-like evolution, suggesting a rapid viral radiation from a single common ancestor with no correspondence between host and viral evolutionary trees. Furthermore, though γ1 ERV env genes do not possess frequent stop codons, γ2 env genes do. To understand whether γ1 suid ERVs may be still replicating, we have also evaluated their likely mechanism of proliferation by statistically testing internal to terminal branches using nonsynonymous versus synonymous substitution ratios. Our results suggest that γ1 ERVs are increasing in copy number by reinfection, which requires the translocation of the virus from one cell to another.
Evidence of at least two viral subpopulations was observed in γ1 ERVs from Eurasian and African host species. These results should be taken into account in xenotransplantation since γ1 ERVs appear to be codiverging with their host and maintaining ongoing capacity to infect somatic and germ cells.
Endogenous retroviruses (ERVs) are the remnants of an exogenous viral form that infected and became integrated in germ cell genomes . ERVs have been found in a variety of vertebrates, including reptiles, amphibian, birds and mammals . However, the occurrence and evolution of ERVs in different organisms are still not well understood and they are sometimes lumped with "junk" DNA of unknown function [3, 4].
In the domestic pig and wild boar (Sus scrofa) genome, ERVs are referred to as PERVs (porcine endogenous retroviruses) and belong to the β and γ genera comprising several viral groups [5, 6]. The γ1 genus is divided in PERV-A, B and C classes on the basis of env sequences, while the γ2 genus comprises only PERV-E. Although most PERVs carry mutations that render them functionally inactive, some γ1 PERVs retain the ability to express and can infect human cells in vitro [7, 8]. The capacity of γ1 ERVs to infect human cells in tissue culture raised serious concerns on the safety of xenotransplantations [9, 10] because pigs are the first choice of xenograft donors [11–13]; consequently γ1 ERVs have been extensively studied in Sus scrofa.
In contrast, few DNA sequences from γ2 ERVs have been reported to date [5, 14, 15], although a full length provirus has been described and named PERV-E because of its similarity to the human HERV-E. Furthermore, PERV-E forms a distinct clade with viruses from the HERV-E family and does not cluster with γ1 PERVs . Genetic and expression analyses also have not indicated an obvious infectious risk from γ2 retroviruses for xenotransplantation .
At present, little is known about either γ1 or γ2 ERVs in other suid species than Sus scrofa [5, 16]. Patience et al.  screened by PCR for the presence of ERVs (both γ and β genera) in five different suid species, while Niebert and Tönjes  sequenced γ1 DNA fragments in seven different suid species. However, they did not perform any evolutionary analyses of ERVs . Clearly both of these studies have left many questions unanswered.
The mechanisms by which γ1 ERVs increase their copy number in suid genomes are also poorly understood and three main mechanisms have been reported for human ERVs (HERVs) [17, 18]; retrotransposition in cis, complementation in trans, and reinfection . The first mechanism involves replication by co-packaging mRNAs in virus-like particles built up from viral proteins encoded by their own gag and pol genes. As retrotransposition in cis does not result in cell to cell translocation, the env gene is not required to be functional. Conversely, complementation in trans involves replication of two retroviral elements of the same family, with one defective provirus using the functional protein produced by another provirus. Finally, ERVs might form infectious exogenous particles and reinfect either a germ cell or a somatic cell. Because reinfection requires the translocation of viruses from one cell to another, the presence of functional genes is required, although this may also result from complementation in trans [17, 20]. New ERVs will likely acquire new mutations during reinfection or retrotransposition due to their error prone reverse transcriptase in virus replication . Thus proviruses might occasionally give rise to infectious exogenous forms which can reinfect the animal itself or animals of the same species and sometimes cross to other species, and subsequently integrate in a germ cell .
Recently a study of the host phylogeny [22, 23] clarified the relationship between Eurasian and African suids based on nuclear and mitochondrial DNA sequences, enabling the study of ERV evolution in these suid species and the evolutionary congruence of host and viral lineages. Here, we analyzed the γ1 ERV coding domains (gag, pro/pol and env) sampled from 10 different suid host species from Eurasia and Africa. We have analyzed their evolutionary relationships as well as their likely mechanism of replication. We have also compared the evolution of γ1 and γ2 ERVs using their env genes and implication for xenotransplantation.
Screening of ERV sequences
List of species of Suidae showing their origin and number of sequences generated.
Species and subspecies
Red river hog
Duisburg Zoo, Germany
Yorkshire Farm, UK
Sus barbatus oi
Western bearded pig
Singapore Zoo, Singapore
Sus barbatus barbatus
Bornean bearded pig
Singapore Zoo, Singapore
Sulawesi warty pig
Javan warty pig
Poznan Zoo, Poland
Surabaya Zoo, Indonesia
Blastn confirmed that amplicons corresponded to ERVs. Primers targeting env classes A and C cross amplified. Sequences generated with these primers were initially analyzed together and then separated into A and C according to blastn results. After reorganization, a total of 18 and 22 sequences for class A and C env were analyzed. env class C was not amplified from any African specimens, obviating the need for the non-parametric test for panmixia for this gene.
Three gag and one env class E sequences showed indels longer than 50 bp and were removed from alignments to increase accuracy of evolutionary analysis. This allowed a higher number of base pairs to be analyzed because these indels were observed in different positions. After removing these sequences, a total of 141 novel γ1 sequences were analyzed along with 28 sequences from the draft pig genome and 110 sequences from GenBank (Additional file 1). These sequences showed a pairwise genetic distances in expected number of substitutions per site ranging from 0 - 0.10 for pol, 0 - 0.11 for gag, 0 - 0.17 for env class A, 0 - 0.27 for env class B, and 0 - 0.15 for env class C. Identical sequences were rarely observed, and they mainly reflect identical sequences from GenBank and between Genbank and Sus scrofa genome sequences. Non-parametric tests revealed highly significant differentiation between ERVs from African versus Eurasian species (P < 0.001). On the other hand, no sequence was removed from γ2 env dataset for which pairwise genetic distances ranged from 0.001 - 0.11. The non-parametric test of population subdivision suggested that Eurasian, African and babirusa sequences came from different subpopulations (P < 0.02).
In addition, most DNA sequences of γ1 ERV genes lacked stop codons or, when present, were limited to one stop codon per sequence. This contrasted with γ2 env, which possessed several stop codons (up to 5 per sequence) and the presence of several indels, some of which were in unique positions and others at conserved positions across Eurasian, African and babirusa sequences, suggesting that they are unlikely to be replicating as env is necessary for viral reinfection. Possibly, they might retrotranspose in cis if other coding domains were devoid of stop codons, although this is apparently unlikely because Sus scrofa gag and pro/pol also showed stop codons and frame shift mutations. Interestingly, there were three γ2 sequences from African species (two from common warthog and one from desert warthog) and two from babirusa that did not show stop codons.
Test for detecting recombination
From 17% to 40% of sequences were identified as recombinants for γ1 genes by the PHI-NNet algorithm , and a higher frequency was mainly observed in sequences from African suids but less frequently than with RDP 3 (Figure 1). Apparent inter-group recombinants between Eurasian and African sequences were identified only for env, albeit with a low frequency (Figures 1C and 1D). It is interesting to note that approximately 50% of env class A sequences in African suids appeared to be recombinants (Figure 1C) while all env class C sequences, found only in Eurasian pigs, were likely to be recombinants and thus could not be used for evolutionary analyses. Because this algorithm identified a smaller number of recombinants, we used the observed non-recombinant dataset in our phylogenetic reconstructions.
The clades shown for African class A and B env phylogenies (Figures 3 and 4) were generally coincident with African host clades, except for the internal arrangement of Potamochoerus porcus and Hylochoerus meinertzhageni in the env class B phylogeny (Figure 4), which differed from the host tree [22, 23] where Hylochoerus grouped with Phacochoerus instead of Potamochoerus. On the other hand, Eurasian γ1 ERVs showed a less conserved arrangement in all phylogenies, with internal groups not corresponding to the host phylogeny [22, 23] except for the env class B internal arrangement, which was very similar in the host phylogeny (highlighted in Figure 4).
Mechanism of ERV proliferation
Most sequences of γ1 ERVs did not possess stop codons and/or indels. Furthermore, when present, they were usually located at unique positions in any given sequence. Low frequency of stop codons is evidence that complementation in trans has been a rare event .
d N /d S (ω) ratios and likelihood ratio tests (LRT) for ERV genes in the family Suidae to test whether ω for internal branches are significantly less than 1.
env (class A)
env (class B)
For the γ2 env gene, we found statistical support for the "one-ratio" model (P = 0.053), which made testing the "two-ratio" model unnecessary.
We found strong evidence that γ1 ERVs from Eurasian and African species based on the pol and env phylogenies belong to two different lineages. This pattern was also observed for LTR phylogenies sampled from Eurasian and African host species . Host species from Eurasia are more similar to one another, showing shorter branch lengths compared to the more differentiated African host species with longer branch lengths, consistent with traditional taxonomy which puts all Eurasian species in the genus Sus, whereas there are several genera recognized in Africa [32–34]. A similar pattern was generally observed in γ1 ERV class A and B env phylogenies and for the LTRs, suggesting that proviruses might be codiverging with their respective host species.
The resolution of the pol phylogeny into distinct Eurasian and African lineages, which was not found with gag, is not a consequence of analysis of longer alignments. Bayesian analysis based on 885 bp of pol sequence (the same size as the gag alignment) showed an almost identical topology to that observed in Figure 5, showing a single African clade of ERVs (including the one ERV sequence of the draft pig genome). Similarly, when ca. 400 bp of env is analysed, an almost identical topology of that of Figure 3 is also observed .
The gag phylogeny is different from the other viral gene phylogenies and to the host phylogeny [22, 23]. This may be explained by the gene tree versus species tree problem in which gene phylogenies sometimes conflict with a species tree because individual gene sequences can generate different topologies and consequently conflicting results [35–37], but is also likely to be due simply to the lack of resolution in the gag phylogeny which cannot be attributed to the length of base pairs analysed compared to the other genes. The lack of resolution of this gene was also observed by other authors when analysing other retroviruses [for an example see [38, 39]]. Alternatively, the DNA sequence conservation of gag gene suggests that Gag proteins are important for viral replication. Experimental analysis of Moloney murine leukemia virus showed that mutations of several portions of gag were incompatible with viral replication , implying strong selection for conservation of the gag gene sequence. In this case, while other viral domains might have evolved distinct African and Eurasian sequences, gag might have been constrained from doing so.
Class A and B env phylogenies showed two likely subpopulations in African suids, one in Phacochoerus species, and another in Potamochoerus and Hylochoerus species more related to ERVs from Sus species, (Figures 3 and 4). Amplification of orthologous sequences as well as full length proviral sequences would provide a better basis for comparison, but for reasons of practicality and convenience paralogous sequences amplified by PCR have been routinely used to determine the evolutionary relationship of ERVs in different host species [41–43].
Recombination was apparent in several γ1 and γ2 sequences, presumably as a result of co-packaging of different RNAs in the same viral particle. In Sus scrofa, recombination between highly variable env sequences is easier to detect than for gag and pol. Moreover, selection would be also more likely to be operating on env than on the other genes because the viral envelope must evolve to escape immune detection, thus making recombination a likely source of adaptative phenotypes. ERVs from African species showed more evidence of recombination than ERVs from Eurasian species, although this apparent difference may be an artifact of the higher differentiation (longer branch lengths) for pol and for class A and B env in African species. Furthermore, the RDP 3 software identified a higher number of recombinants than the PHI-NNet algorithm for γ1 genes, a fact also observed in HIV by Lamers et al. (2009) and attributed to the molecular models implemented in some RDP 3 methods, to the number of sequences in each analysis and to variation within each subpopulation altering the number of identified recombinants . This might be why the γ1 pol sequence from the draft pig genome was identified as recombinant by RDP 3 and non-recombinant by the PHI-NNet algorithm. These discordances obviously result from the different methods of analyses, pointing to the need of improving their consistency and power especially in detecting recombinants between ERV sequences. This improvement is also important to better characterize the inter subpopulation recombinants observed in this study and confirm that they represent genuine recombinant sequences.
We also found evidence for past purifying selection on all γ1 genes. This suggested that retroviruses that remained functional were initially selected and that transposition by reinfection was relevant, while loss of function would have likely resulted in neutral evolution within the host genome ω ≈ 1. Furthermore, the rarity of stop codons and/or indels also suggested that complementation in trans and retrotransposition in cis were not relevant mechanisms of transposition [17, 18].
Evolutionary comparisons of γ1 and γ2 envgenes
In contrast to the γ1 phylogeny, γ2 ERV env gene showed a bush-like unresolved tree, (Figure 7) confirmed by the likelihood-mapping, either resulting from inadequate sampling or data (soft polytomy) or reflecting the actual evolutionary history of γ2 ERVs in the Suidae (hard polytomy) [45–47]. The bush-like pattern might reflect an incomplete sampling because env from ten other suid species is missing in the phylogeny. However, the γ1 ERV env phylogenies (Figures 3 and 4) using fewer species showed a better resolved tree with species from Eurasia grouping apart from African suids, suggesting that the number of viral sequences samples from each suid host would not be an explanation for the poor resolution of γ2 ERV env phylogeny.
An unresolved phylogeny, like that observed for γ2 ERVs, may also result from data saturation, but the env sequences showed very low saturation. Although we cannot rule out the possibility that the γ2 star-like phylogeny resulted from inadequate sampling, this seems unlikely as samples used here were available from ERVs from three different suid lineages observed by Nascimento [22, 23]. In this case, a hard polytomy seems more plausible. Although Vandamme  emphasised the difficulty of proving the existence of hard polytomies in real life, Poe and Chubb , Barth et al. , Willerslev et al.  and several others [47, 50] have observed some apparent examples in different taxa, including birds, mammals, plants and protozoa, and this pattern was also observed for human ERVs .
This hard polytomy may then represent a rapid radiation from a single common ancestor giving rise to multiple distinct retroviral lineages almost at the same time [19, 45]. These ERVs have subsequently become inactivated resulting in the loss of the last active lineage (also suggested by the presence of several stop codons), following a long period of inactivity leading to the star shaped phylogeny .
The reason why distinct evolutionary histories are observed for γ1 and γ2 ERVs is not known. It would be reasonable to attribute it to the evolution of the host defensive mechanisms. However, Katzourakis et al.  by modelling ERV evolutionary dynamics demonstrated that such a star shaped phylogeny can be generated by a "null model" in which all parameters are set constant through out time.
The different lineages of γ1 ERVs in the Suidae are generally congruent with the host phylogeny. This implies that horizontal transmission across these very different host lineages has not occurred. Moreover, the frequent evolutionary occurrence of viral recombination may be relevant to risk analysis of xenotransplantation, especially since viral recombination will favour adaptation and confer capacity for maintaining ongoing infection of somatic and germ cells. This highlights the importance of improving methods to prevent PERVs from crossing the species barriers. Finally, the evolution of γ2 ERVs confirms their lack of risk in xenotransplantation because the star phylogeny suggests loss of the last active element which is also confirmed by the presence of numerous stop codons.
Sample collection, primer design and PCR assays
Oligonucleotide primer used to amplify and to sequence partial ERV genes.
Primer sequences (5' → 3')
Expected size (bp)
Cloning and sequencing
Purified PCR products were cloned using the TOPO TA cloning® kit (Invitrogen, Australia) and DNA from plasmids containing inserts was subsequently extracted with UltraClean™ Mini Plasmid Prep kit (Mo Bio). For each specimen, approximately five purified plasmids were sent to the Australian Genome Research Facility Ltd (AGRF; Brisbane, Australia) for DNA sequencing with primers listed in Table 3. Electropherograms were checked using BIOEDIT (version 184.108.40.206) .
DNA sequences were blasted using the blastn option (http://blast.ncbi.nlm.nih.gov/Blast.cgi) to confirm that ERV products had been amplified and to identify all previously reported PERV sequences. The pre-ensembl Sus scrofa genome database (http://pre.ensembl.org/index.html; version Sscrofa8) was searched using the blat option for gag, pro/pol and env, with best matching sequences being also included in alignments. Only sequences with 100% coverage were used in subsequent analyses. GenBank sequences were from genomic PERV sequences from pigs of various breeds, porcine cell lines, infected human primary cells and transcripts from virus particles released from porcine cell lines. The majority of GenBank sequences do not contain information on genomic position and we could not take this information into account. Sequences from the Sus scrofa genome were also not representative of the numerous PERV loci already described. This was because the Sscrofa8 version of the pig genome was not finished and because we have narrowed our search to sequences with 100% coverage.
Alignment and genetic distance estimations
All novel and available PERV sequences from GenBank and from the draft pig genome were aligned with MUSCLE (version 3.6) , with default options and were manually checked. Columns showing gaps in the majority of aligned sequences and nucleotides at the minority of sequences were manually removed. Conversely, columns showing gaps at sites where nucleotides were present in the majority of sequences were not removed, and these gaps were treated as missing data . pol and gag sequences from murine leukemia virus (AY277737 and EU075329 respectively), gibbon leukemia virus (U60065 and M26927 respectively) and koala retrovirus (AF151794) were used as outgroups. We initially aligned amino acid sequences of outgroups and suid ERVs, and subsequently back-translated to nucleotide using SEAVIEW (version 4.0) . We did not include outgroups when analyzing the env genes because of the very low conservation of ERV env sequences. Pairwise genetic distances were estimated using the modified Log-Det implemented in MEGA (version 4)  in view that it can estimate reliable distances for closely related taxa . Identical sequences, very similar (pairwise genetic distances equal to 0.001) and those with large indels were removed from alignments for a better phylogeny estimation, recombination detection and calculation of ω (see below).
Test of subpopulation structure and detection of recombination
To test whether viral sequences isolated from Eurasian, African and the babirusa (when available) suids belong to statistically different subpopulations, we explored the datasets with a test for population sub-division developed by Hudson et al.  and adapted to study HIV populations by Achaz et al. . This non-parametric test of population subdivision was performed using a web-based interface (http://wwwabi.snv.jussieu.fr/achaz/hudsontest.html) with a nominal significance threshold of 0.05 for accepting or rejecting a null hypothesis for structure. Because recombinants may misplace sequences in a phylogeny [59, 60], they were detected using the PHI-NNet algorithm . This could efficiently detect recombination among closely related sequences  from Eurasian host species, African host species and the babirusa (when available) and could even suggest Eurasian-African (-babirusa) recombinants. Recombinants were also independently detected using the software RDP 3 which implements seven different recombination programs, 1) the original RDP method , 2) the GENECONV method [62, 63], 3) the MaxChi method , 4) the Chimaera method , 5) the SiScan method  and 7) 3SEQ method . Default settings and the option "auto mask for optimal recombination detection" were used in all analyses.
Bayes factor comparisons between model M0 and M1 for each ERV gene showing also the number of generations for each Bayesian inference.
Number of generations
env (class A)
env (class B)
env (class E)
Test of substitution saturation
Because the γ1 gag and the γ2 env sequences showed unresolved phylogenies, we have calculated an index of DNA substitution saturation (Iss)  using the program DAMBE (version 5.0.85)  to evaluate whether sequences contained phylogenetic signals. Furthermore, because of the star-like γ2 env phylogeny, a likelihood-mapping analysis  was performed using TREE-PUZZLE which shows the percentage of tree-like, star-like and net-like signals in the dataset [for more information, see ].
Estimations of ω ratios
To evaluate whether γ1 ERVs may be still replicating  we have calculated the nonsynonymous (d N ) versus synonymous (d S ) substitution ratios (ω = d N /d S ) for γ1 genes using PAML (version 4.2)  according to Belshaw et al. . Internal branches were compared to terminal branches using the "two-ratio" model  where the null hypothesis of ω ≥ 1 (neutral or positive selection) was compared with the alternative hypothesis of ω < 1 (purifying selection). Significance was tested by fixing ω at 1 for internal branches and performing a likelihood ratio test (LRT), where twice the difference of the log-likelihoods was compared to critical values of the χ2 distribution for 1 degree of freedom .
For γ2 env gene, we have compared a single ω for the entire tree (the "one-ratio" model) to a "free-ratio" model (a separate ω ratio for each branch of the phylogeny) . The significance of differences between "one-ratio" and "free-ratio" models was assessed by a LRT for 34 degrees of freedom .
Outgroups were removed from the trees before performing these calculations.
GenBank accession numbers
env class E: [GQ906159-GQ906198]; env class A: [GQ906199-GQ906216]; env class B: [GQ906217-GQ906249]; env class C: [GQ906250-GQ906271]; gag: [GQ906272-GQ906309]; pol: [GQ906310-GQ906342] for sequences generated in this study. For sequences retrieved from GenBank see additional file 1.
Acknowledgements and Funding
F. F. Nascimento was a grantee of EIPRS and Stock and Meat Industries Grant-in Aid (Australia), and Boehringer Ingelheim Fonds (Germany). M. Charleston gratefully acknowledges funding from the Australian Research Council (grant number DP1094891). We are grateful for H. Klingel (Ugandan Institute of Ecology, Uganda, and University of Braunschweig, Germany) and to E. Randi (Istituto Superiore per la Protezione e la Ricerca Ambientale, Italy) for providing the DNA samples of the forest hog and the desert warthog respectively. We are thankful to H. Seuánez (Instituto Nacional de Câncer, Brazil) and G. Larson (Durham University, United Kingdom) for reviewing earlier drafts of this manuscript.
- Gifford R, Tristem M: The evolution, distribution and diversity of endogenous retroviruses. Virus Genes. 2003, 26: 291-315. 10.1023/A:1024455415443.View ArticlePubMedGoogle Scholar
- Herniou E, Martin J, Miller K, Cook J, Wilkinson M, Tristem M: Retroviral diversity and distribution in vertebrates. J Virol. 1998, 72: 5955-5966.PubMedPubMed CentralGoogle Scholar
- Weiss RA: The discovery of endogenous retroviruses. Retrovirology. 2006, 3: 67-10.1186/1742-4690-3-67.View ArticlePubMedPubMed CentralGoogle Scholar
- Biemont C, Vieira C: Genetics: junk DNA as an evolutionary force. Nature. 2006, 443: 521-524. 10.1038/443521a.View ArticlePubMedGoogle Scholar
- Patience C, Switzer WM, Takeuchi Y, Griffiths DJ, Goward ME, Heneine W, Stoye JP, Weiss RA: Multiple groups of novel retroviral genomes in pigs and related species. J Virol. 2001, 75: 2771-2775. 10.1128/JVI.75.6.2771-2775.2001.View ArticlePubMedPubMed CentralGoogle Scholar
- Klymiuk N, Müller M, Brem G, Aigner B: Characterization of porcine endogenous retrovirus gamma pro-pol nucleotide sequences. J Virol. 2002, 76: 11738-11743. 10.1128/JVI.76.22.11738-11743.2002.View ArticlePubMedPubMed CentralGoogle Scholar
- Martin U, Kiessig V, Blusch JH, Haverich A, von der Helm K, Herden T, Steinhoff G: Expression of pig endogenous retrovirus by primary porcine endothelial cells and infection of human cells. Lancet. 1998, 352: 692-694. 10.1016/S0140-6736(98)07144-X.View ArticlePubMedGoogle Scholar
- Takeuchi Y, Patience C, Magre S, Weiss RA, Banerjee PT, Le Tissier P, Stoye JP: Host range and interference studies of three classes of pig endogenous retrovirus. J Virol. 1998, 72: 9986-9991.PubMedPubMed CentralGoogle Scholar
- Patience C, Takeuchi Y, Weiss RA: Infection of human cells by an endogenous retrovirus of pigs. Nat Med. 1997, 3: 282-286. 10.1038/nm0397-282.View ArticlePubMedGoogle Scholar
- Stoye JP, Coffin JM: The dangers of xenotransplantation. Nat Med. 1995, 1: 1100-View ArticlePubMedGoogle Scholar
- Platt JL: Xenotransplantation. New risks, new gains. Nature. 2000, 407: 29-30. 27View ArticleGoogle Scholar
- Boneva RS, Folks TM: Xenotransplantation and risks of zoonotic infections. Ann Med. 2004, 36: 504-517. 10.1080/07853890410018826.View ArticlePubMedGoogle Scholar
- Smetanka C, Cooper DK: The ethics debate in relation to xenotransplantation. Rev Sci Tech. 2005, 24: 335-342.PubMedGoogle Scholar
- Mang R, Maas J, Chen X, Goudsmit J, van Der Kuyl AC: Identification of a novel type C porcine endogenous retrovirus: evidence that copy number of endogenous retroviruses increases during host inbreeding. J Gen Virol. 2001, 82: 1829-1834.View ArticlePubMedGoogle Scholar
- Klymiuk N, Muller M, Brem G, Aigner B: Phylogeny, recombination and expression of porcine endogenous retrovirus gamma2 nucleotide sequences. J Gen Virol. 2006, 87: 977-986. 10.1099/vir.0.81552-0.View ArticlePubMedGoogle Scholar
- Niebert M, Tönjes RR: Evolutionary spread and recombination of porcine endogenous retroviruses in the Suiformes. J Virol. 2005, 79: 649-654. 10.1128/JVI.79.1.649-654.2005.View ArticlePubMedPubMed CentralGoogle Scholar
- Belshaw R, Katzourakis A, Paces J, Burt A, Tristem M: High copy number in human endogenous retrovirus families is associated with copying mechanisms in addition to reinfection. Mol Biol Evol. 2005, 22: 814-817. 10.1093/molbev/msi088.View ArticlePubMedGoogle Scholar
- Belshaw R, Pereira V, Katzourakis A, Talbot G, Paces J, Burt A, Tristem M: Long-term reinfection of the human genome by endogenous retroviruses. Proc Natl Acad Sci USA. 2004, 101: 4894-4899. 10.1073/pnas.0307800101.View ArticlePubMedPubMed CentralGoogle Scholar
- Katzourakis A, Rambaut A, Pybus OG: The evolutionary dynamics of endogenous retroviruses. Trends Microbiol. 2005, 13: 463-468. 10.1016/j.tim.2005.08.004.View ArticlePubMedGoogle Scholar
- Finnegan DJ: Transposable elements. Curr Opin Genet Dev. 1992, 2: 861-867. 10.1016/S0959-437X(05)80108-X.View ArticlePubMedGoogle Scholar
- Jern P, Coffin JM: Effects of retroviruses on host genome function. Annu Rev Genet. 2008, 42: 709-732. 10.1146/annurev.genet.42.110807.091501.View ArticlePubMedGoogle Scholar
- Nascimento FF: Evolution of the family Suidae and their endogenous retroviruses. PhD. 2009, The University of SydneyGoogle Scholar
- Gongora J, Cuddahee RE, Nascimento FF, Palgrave CJ, Lowden S, Ho SYW, Simond D, Damayanti CS, White DJ, Tay WT, Randi E, Klingel H, Rodrigues-Zarate CJ, Allen K, Moran C, Larson G: Rethinking the evolution of sub-Saharan African suids (Suidae, Artiodactyla). Zool Scr. 2011Google Scholar
- Martin DP, Williamson C, Posada D: RDP2: recombination detection and analysis from sequence alignments. Bioinformatics. 2005, 21: 260-262. 10.1093/bioinformatics/bth490.View ArticlePubMedGoogle Scholar
- Salemi M, Gray RR, Goodenow MM: An exploratory algorithm to identify intra-host recombinant viral sequences. Mol Phylogenet Evol. 2008, 49: 618-628. 10.1016/j.ympev.2008.08.017.View ArticlePubMedPubMed CentralGoogle Scholar
- Xia X, Xie Z, Salemi M, Chen L, Wang Y: An index of substitution saturation and its application. Mol Phylogenet Evol. 2003, 26: 1-7. 10.1016/S1055-7903(02)00326-3.View ArticlePubMedGoogle Scholar
- Xia X: Assessing substitution saturation with DAMBE. The phylogenetic handbook : a practical approach to phylogenetic analysis and hypothesis testing. Edited by: Salemi M, Vandamme A-M, Lemey P. 2009, Cambridge: Cambridge University Press, 615-623. 2View ArticleGoogle Scholar
- Strimmer K, von Haeseler A: Likelihood-mapping: a simple method to visualize phylogenetic content of a sequence alignment. Proc Natl Acad Sci USA. 1997, 94: 6815-6819. 10.1073/pnas.94.13.6815.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang Z, Bielawski JP: Statistical methods for detecting molecular adaptation. Trends Ecol Evol. 2000, 15: 496-503. 10.1016/S0169-5347(00)01994-7.View ArticlePubMedGoogle Scholar
- Yang Z: Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. Mol Biol Evol. 1998, 15: 568-573.View ArticlePubMedGoogle Scholar
- Nascimento FF, Gongora J, Tristem M, Lowden S, Moran C: Distinctive differences in long terminal repeat sequences between gamma1 endogenous retroviruses of African and Eurasian suid species. Infect Genet Evol. 2011, 11: 686-693. 10.1016/j.meegid.2011.01.010.View ArticlePubMedGoogle Scholar
- Harris JM, White TD: Evolution of the Plio-Pleistocene African Suidae. Trans Am Phil Soc. 1979, 69: 1-128.View ArticleGoogle Scholar
- White TD, Harris JM: Suid evolution and correlation of African hominid localities. Science. 1977, 198: 13-21. 10.1126/science.331477.View ArticlePubMedGoogle Scholar
- Pickford M: Old world suoid systematics, phylogeny, biogeography and biostratigraphy. Paleontologia i Evolució. 1993, 26-27: 237-269.Google Scholar
- Jeffroy O, Brinkmann H, Delsuc F, Philippe H: Phylogenomics: the beginning of incongruence?. Trends Genet. 2006, 22: 225-231. 10.1016/j.tig.2006.02.003.View ArticlePubMedGoogle Scholar
- Rannala B, Yang Z: Phylogenetic inference using whole genomes. Annu Rev Genomics Hum Genet. 2008, 9: 217-231. 10.1146/annurev.genom.9.081307.164407.View ArticlePubMedGoogle Scholar
- Degnan JH, Rosenberg NA: Gene tree discordance, phylogenetic inference and the multispecies coalescent. Trends Ecol Evol. 2009, 24: 332-340. 10.1016/j.tree.2009.01.009.View ArticlePubMedGoogle Scholar
- Llorens C, Fares MA, Moya A: Relationships of gag-pol diversity between Ty3/Gypsy and Retroviridae LTR retroelements and the three kings hypothesis. BMC Evol Biol. 2008, 8: 276-10.1186/1471-2148-8-276.View ArticlePubMedPubMed CentralGoogle Scholar
- Howard TM, Sheng Z, Wang M, Wu Y, Rasheed S: Molecular and phylogenetic analyses of a new amphotropic murine leukemia virus (MuLV-1313). Virol J. 2006, 3: 101-10.1186/1743-422X-3-101.View ArticlePubMedPubMed CentralGoogle Scholar
- Auerbach MR, Shu C, Kaplan A, Singh IR: Functional characterization of a portion of the Moloney murine leukemia virus gag gene by genetic footprinting. Proc Natl Acad Sci USA. 2003, 100: 11678-11683. 10.1073/pnas.2034020100.View ArticlePubMedPubMed CentralGoogle Scholar
- Jaratlerdsiri W, Rodriguez-Zarate CJ, Isberg SR, Damayanti CS, Miles LG, Chansue N, Moran C, Melville L, Gongora J: Distribution of endogenous retroviruses in crocodilians. J Virol. 2009, 83: 10305-10308. 10.1128/JVI.00668-09.View ArticlePubMedPubMed CentralGoogle Scholar
- Xiao R, Park K, Lee H, Kim J, Park C: Identification and classification of endogenous retroviruses in cattle. J Virol. 2008, 82: 582-587. 10.1128/JVI.01451-07.View ArticlePubMedGoogle Scholar
- Tristem M, Kabat P, Lieberman L, Linde S, Karpas A, Hill F: Characterization of a novel murine leukemia virus-related subgroup within mammals. J Virol. 1996, 70: 8241-8246.PubMedPubMed CentralGoogle Scholar
- Lamers SL, Salemi M, Galligan DC, de Oliveira T, Fogel GB, Granier SC, Zhao L, Brown JN, Morris A, Masliah E, McGrath MS: Extensive HIV-1 intra-host recombination is common in tissues with abnormal histopathology. PLoS One. 2009, 4: e5065-10.1371/journal.pone.0005065.View ArticlePubMedPubMed CentralGoogle Scholar
- Vandamme AM: Basic concepts of molecular evolution. The phylogenetic handbook: a practical approach to phylogenetic analysis and hypothesis testing. Edited by: Salemi M, Vandamme A-M, Lemey P. 2009, Cambridge: Cambridge University Press, 3-29. 2View ArticleGoogle Scholar
- Barth D, Przybos E, Fokin SI, Schlegel M, Berendonk TU: Cytochrome b sequence data suggest rapid speciation within the Paramecium aurelia species complex. Mol Phylogenet Evol. 2008, 49: 669-673. 10.1016/j.ympev.2008.08.007.View ArticlePubMedGoogle Scholar
- Whitfield JB, Lockhart PJ: Deciphering ancient rapid radiations. Trends Ecol Evol. 2007, 22: 258-265. 10.1016/j.tree.2007.01.012.View ArticlePubMedGoogle Scholar
- Poe S, Chubb AL: Birds in a bush: five genes indicate explosive evolution of avian orders. Evolution. 2004, 58: 404-415.View ArticlePubMedGoogle Scholar
- Willerslev E, Gilbert MT, Binladen J, Ho SY, Campos PF, Ratan A, Tomsho LP, da Fonseca RR, Sher A, Kuznetsova TV, Nowak-Kemp M, Roth TL, Miller W, Schuster SC: Analysis of complete mitochondrial genomes from extinct and extant rhinoceroses reveals lack of phylogenetic resolution. BMC Evol Biol. 2009, 9: 95-10.1186/1471-2148-9-95.View ArticlePubMedPubMed CentralGoogle Scholar
- Rokas A, Carroll SB: Bushes in the tree of life. PLoS Biol. 2006, 4: e352-10.1371/journal.pbio.0040352.View ArticlePubMedPubMed CentralGoogle Scholar
- Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999, 41: 95-98.Google Scholar
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang Z: Computational molecular evolution. 2006, Oxford: Oxford University PressView ArticleGoogle Scholar
- Galtier N, Gouy M, Gautier C: SEAVIEW and PHYLO_WIN: two graphic tools for sequence alignment and molecular phylogeny. Comput Appl Biosci. 1996, 12: 543-548.PubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.View ArticlePubMedGoogle Scholar
- Tamura K, Kumar S: Evolutionary distance estimation under heterogeneous substitution pattern among lineages. Mol Biol Evol. 2002, 19: 1727-1736.View ArticlePubMedGoogle Scholar
- Hudson RR, Boos DD, Kaplan NL: A statistical test for detecting geographic subdivision. Mol Biol Evol. 1992, 9: 138-151.PubMedGoogle Scholar
- Achaz G, Palmer S, Kearney M, Maldarelli F, Mellors JW, Coffin JM, Wakeley J: A robust measure of HIV-1 population turnover within chronically infected individuals. Mol Biol Evol. 2004, 21: 1902-1912. 10.1093/molbev/msh196.View ArticlePubMedGoogle Scholar
- Kosiol C, Bofkin L, Whelan S: Phylogenetics by likelihood: evolutionary modeling as a tool for understanding the genome. J Biomed Inform. 2006, 39: 51-61. 10.1016/j.jbi.2005.08.003.View ArticlePubMedGoogle Scholar
- Jermiin LS, Jayaswal V, Ababneh F, Robinson J: Phylogenetic model evaluation. Methods Mol Biol. 2008, 452: 331-364. 10.1007/978-1-60327-159-2_16.View ArticlePubMedGoogle Scholar
- Martin D, Rybicki E: RDP: detection of recombination amongst aligned sequences. Bioinformatics. 2000, 16: 562-563. 10.1093/bioinformatics/16.6.562.View ArticlePubMedGoogle Scholar
- Sawyer S: Statistical tests for detecting gene conversion. Mol Biol Evol. 1989, 6: 526-538.PubMedGoogle Scholar
- Padidam M, Sawyer S, Fauquet CM: Possible emergence of new geminiviruses by frequent recombination. Virology. 1999, 265: 218-225. 10.1006/viro.1999.0056.View ArticlePubMedGoogle Scholar
- Smith JM: Analyzing the mosaic structure of genes. J Mol Evol. 1992, 34: 126-129.PubMedGoogle Scholar
- Posada D, Crandall KA: Evaluation of methods for detecting recombination from DNA sequences: computer simulations. Proc Natl Acad Sci USA. 2001, 98: 13757-13762. 10.1073/pnas.241370698.View ArticlePubMedPubMed CentralGoogle Scholar
- Gibbs MJ, Armstrong JS, Gibbs AJ: Sister-scanning: a Monte Carlo procedure for assessing signals in recombinant sequences. Bioinformatics. 2000, 16: 573-582. 10.1093/bioinformatics/16.7.573.View ArticlePubMedGoogle Scholar
- Boni MF, Posada D, Feldman MW: An exact nonparametric method for inferring mosaic structure in sequence triplets. Genetics. 2007, 176: 1035-1047.View ArticlePubMedPubMed CentralGoogle Scholar
- Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- Keane TM, Creevey CJ, Pentony MM, Naughton TJ, McLnerney JO: Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified. BMC Evol Biol. 2006, 6: 29-10.1186/1471-2148-6-29.View ArticlePubMedPubMed CentralGoogle Scholar
- Nylander JA, Ronquist F, Huelsenbeck JP, Nieves-Aldrey JL: Bayesian phylogenetic analysis of combined data. Syst Biol. 2004, 53: 47-67. 10.1080/10635150490264699.View ArticlePubMedGoogle Scholar
- Tracer v1.4. [http://beast.bio.ed.ac.uk/Tracer]
- Xia X, Xie Z, Salemi M, Chen L, Wang Y: An index of substitution saturation and its application. Mol Phylogenet Evol. 2003, 26: 1-7. 10.1016/S1055-7903(02)00326-3.View ArticlePubMedGoogle Scholar
- Xia X, Xie Z: DAMBE: software package for data analysis in molecular biology and evolution. J Hered. 2001, 92: 371-373. 10.1093/jhered/92.4.371.View ArticlePubMedGoogle Scholar
- Schmidt HA, Strimmer K, Vingron M, von Haeseler A: TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002, 18: 502-504. 10.1093/bioinformatics/18.3.502.View ArticlePubMedGoogle Scholar
- Schmidt HA, von Haeseler A: Phylogenetic inference using maximum likelihood methods. The phylogenetic handbook : a practical approach to phylogenetic analysis and hypothesis testing. Edited by: Salemi M, Vandamme A-M, Lemey P. 2009, Cambridge: Cambridge University Press, 181-198. 2View ArticleGoogle Scholar
- Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24: 1586-1591. 10.1093/molbev/msm088.View ArticlePubMedGoogle Scholar
- Hasegawa M, Kishino H, Yano T: Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol. 1985, 22: 160-174. 10.1007/BF02101694.View ArticlePubMedGoogle Scholar
- Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980, 16: 111-120. 10.1007/BF01731581.View ArticlePubMedGoogle Scholar
- Kass RE, Raftery AE: Bayes factors. J Am Stat Assoc. 1995, 90: 773-795. 10.2307/2291091.View ArticleGoogle 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.