- Open Access
Discovery of prosimian and afrotherian foamy viruses and potential cross species transmissions amidst stable and ancient mammalian co-evolution
- Aris Katzourakis†1Email author,
- Pakorn Aiewsakun†1,
- Hongwei Jia2,
- Nathan D Wolfe3, 4, 5,
- Matthew LeBreton6,
- Anne D Yoder7 and
- William M Switzer2Email author
© Katzourakis et al.; licensee BioMed Central Ltd. 2014
Received: 19 February 2014
Accepted: 2 July 2014
Published: 4 August 2014
Foamy viruses (FVs) are a unique subfamily of retroviruses that are widely distributed in mammals. Owing to the availability of sequences from diverse mammals coupled with their pattern of codivergence with their hosts, FVs have one of the best-understood viral evolutionary histories ever documented, estimated to have an ancient origin. Nonetheless, our knowledge of some parts of FV evolution, notably that of prosimian and afrotherian FVs, is far from complete due to the lack of sequence data.
Here, we report the complete genome of the first extant prosimian FV (PSFV) isolated from a lorisiforme galago (PSFVgal), and a novel partial endogenous viral element with high sequence similarity to FVs, present in the afrotherian Cape golden mole genome (ChrEFV). We also further characterize a previously discovered endogenous PSFV present in the aye-aye genome (PSFVaye). Using phylogenetic methods and available FV sequence data, we show a deep divergence and stable co-evolution of FVs in eutherian mammals over 100 million years. Nonetheless, we found that the evolutionary histories of bat, aye-aye, and New World monkey FVs conflict with the evolutionary histories of their hosts. By combining sequence analysis and biogeographical knowledge, we propose explanations for these mismatches in FV-host evolutionary history.
Our discovery of ChrEFV has expanded the FV host range to cover the whole eutherian clade, and our evolutionary analyses suggest a stable mammalian FV-host co-speciation pattern which extends as deep as the exafroplacentalian basal diversification. Nonetheless, two possible cases of host switching were observed. One was among New World monkey FVs, and the other involves PSFVaye and a bat FV which may involve cross-species transmission at the level of mammalian orders. Our results highlight the value of integrating multiple sources of information to elucidate the evolutionary history of viruses, including continental and geographical histories, ancestral host locations, in addition to the natural history of host and virus.
Foamy viruses (FVs) are complex retroviruses in the Spumaretrovirinae subfamily . All known contemporary FVs cause persistent but non-virulent, asymptomatic infections exclusively among boreoeutherian mammals (Additional file 1: Table S1) . The first FV was discovered in a macaque in 1954  and was isolated in 1955 . Shortly thereafter, numerous FVs were isolated from other boreoeutherian mammals, including non-human primates (NHPs) [4–13], cats , cattle , and horses , for most of which complete genomes are available [16–25]. Recently, a novel FV was discovered by metagenomics in a bat (Rhinolophus affinis) from China (RhiFV) and was partially sequenced ; however, whether RhiFV is bat-specific will require confirmation with population and epidemiological studies. Although virtually every NHP sampled has a species-specific FV, a human-specific FV has not been identified. Moreover, all FVs isolated from humans have been demonstrated to originate from zoonotic infections with simian FVs (SFVs) [27–36].
Similar to other retroviruses, FVs occasionally leave ‘fossils’ in host genomes, which are the relics of past infections. These genomic fossils are called endogenous retroviruses (ERVs) and result from germ line infections that are followed by vertical inheritance and thus are present in every cell of an infected organism. In contrast, exogenous retroviruses and hence exogenous FVs are transmitted horizontally from organism-to-organism via infectious body fluids and may be confined to certain cell types defined by their host tropism. Only two endogenous FVs have been discovered in mammalian genomes to date. One is present in the two-toed sloth (Choloepus hoffmanni) genome (SloEFV) which was found to have a genome organization characteristic of all contemporary FVs . The second endogenous FV was found in the aye-aye (Daubentonia madagascariensis) genome (PSFVaye) which to date has only been partially characterized . Given a broad range of FV mammalian hosts and the wealth of available whole mammalian genomes, the finding of only two endogenous FVs so far suggests that FV endogenization is a rare event. While a number of defective ERVs found in fish genomes have been reported to exhibit some similarity to known exogenous FVs [39–41], more evidence is required to definitively determine if these ERVs are truly endogenous FVs, or alternatively, are distinct lineages that branched early on in FV or retroviral evolution.
Owing to the availability of sequences from diverse mammals (Additional file 1: Table S1), coupled with the stable pattern of codivergence with their hosts, the evolutionary history of mammalian FVs can be studied in great depth. For example, phylogenetic analysis of extant SFVs revealed strong evidence of SFV-host co-speciation over more than 30 million years (Myr) . The discovery and phylogenetic analysis of SloEFV  supported and extended the FV-host co-speciation hypothesis further to the basal radiation of eutherians which occurred more than 100 Myr ago (Ma)  and simultaneously expanded FV mammalian host range to cover xenarthrans, which are one of the four superorders of the Eutheria clade (Laurasiatheria, Euarchontoglires and Afrotheria are the other superorders). This is in contrast to lentiviruses, for which sequence availability is comparable to that of FVs, but their evolutionary history is not as well understood and much more difficult to investigate, as they do not always co-speciate with their hosts . Nevertheless, at the present day, evidence of FV infection in, and sequence data from prosimians and afrotherian mammals is still lacking, and this is required to better understand the FV radiation and evolutionary history in eutherians. Here, we report the discovery and describe the first extant exogenous PSFV isolated from a lorisiforme galago (Otolemur crassicaudatus panganiensis) (PSFVgal) and a novel ERV, that is potentially an endogenous FV, present in the genome of an afrotherian Cape golden mole (Chrysochloris asiatica) (ChrEFV). We also further describe PSFVaye that was previously only partially characterized . Together with the currently available FV sequences, we use detailed molecular analyses to re-examine the co-evolutionary history of mammalian hosts and their FVs.
Characterization of the complete genome of the extant PSFVgal from a galago
Sequence similarity analysis revealed that the LTR length is 1,267 nt. An asparagine-1,2 tRNA-utilizing primer binding site (PBS) was identified downstream of the 5′-LTR before gag (TGGCGTCCCTGGGTGGGC, nt 1,270-1,287), but which is tRNALys (TGGCGCCCAACGTGGGGG/C) in all other mammalian FVs. We confirmed the tRNAAsp PBS sequences by population-based sequencing of multiple PCR products amplified on different dates using primers surrounding the PBS with two different primer pairs and two different concentrations (0.1 and 1.0 ug) of day 52 tissue PSFVgal-infected HeLa cell DNA. By comparison with the SFVmac LTR , the cap site defining the 5′ boundary of the R region was determined to be 20 nt downstream from the TATA box (TATATAA, nt 928–934) at nt 955 within the 5′-LTR, and the polyadenylation site within the 3′-LTR, which defines the 3′ boundary of the R region, was located at the CA dinucleotide (nt 11,955) 20 nt downstream from the polyadenylation signal (AATAAA, nt 11,930-11,935). Thus, the PSFVgal LTR U3, R, and U5 regions are 954, 150, and 163 nt long, respectively. Two polypurine tracts (PPT: AAGGAGAGGG) for the dual initiation of plus-stand DNA synthesis  were located at the 5′ boundary of the 3′-LTR (nt 10,842-10,853) and at the center of the genome toward the 3′ terminus of pol (nt 6,068-6,077) as anticipated.
The gag gene was predicted to be 1,932 nt long (nt 1,363-3,294) to generate a 644 amino acid (aa) protein. The cytoplasmic retention signal (CTRS), important for particle assembly and viral budding , was located at aa 73–90 (GNWGPGDRFARIEVLLRD). The position of a late-assembly domain proline rich PSAP motif, which is essential for primate FV particle assembly and budding [47, 48], was located at aa positions 217–220. The PSAP motif is highly conserved in all SFVs but its position within Gag is variable . The three FV-specific glycine-arginine rich (GR I-III) motifs or boxes within the C-terminal region of Gag  were also determined (aa 479–493, 541–563, 567–597, respectively).
The Pol protein is 1,156 aa long (nt 3,248-6,718) and protease (PR), reverse transcriptase (RT), RNase H, and integrase (IN) were predicted to be located at aa positions 10–168, 154–355, 581–737, and 742–1,140, respectively. The catalytic centers of PR (DTG) and RT (YVDD) were located at aa 20–22 and 304–307, respectively. RNase H active sites were found to be at D589 and D659 and IN catalytic motif (DD35E) within the IN core domain (aa 862–972) was at D926-D938-E962.
The Env protein is 999 aa long (nt 6,603-9,602). We found a highly conserved WXXW motif (WLRW, aa 10–13) in the N-terminus which is essential for Gag-Env protein interaction during the budding process and is found in all known extant FVs [50, 51]. The internal promoter (TATAAAA) within the env gene, which serves for the initiation of bel1 transcription, is present at nt 9,299-9,305. We also identified a hydrophobic transmembrane region, which contains a putative signal-peptide peptidase cleavage site (TMGWCIGLFCLLLILLFS↓LVIVIL), at aa 81–104. The cleavage occurs within the endoplasmic reticulum to remove the signal peptide . A potential furin-cleavage site was also found (TRPNYTAARSRR↓SVE) at aa 131–145. It is this site at which the FV Env protein is cleaved, either by furin or furin-like proteases, to produce an Env leader protein – a component that interacts with the N-terminus of Gag via its WXXW motif and is absolutely required for FV particle budding .
The bel1 and bel2 genes were predicted to be present at nt 9,569-10,402 and nt 10,026-11,414, respectively. The putative Bel1 protein (277 aa) exhibits only a weak similarity (<15% identity) to extant FV Bel1 proteins, which could only be shown by comparison with Hidden Markov Model (HMM) of FV Bel1 proteins (HMM profile GyDB (http://www.gydb.org/): Bel1 spumaretroviridae; E = 3.7E-06). In contrast, the putative Bel2 protein (511 aa) exhibits a relatively high similarity with other FV Bel2 proteins with the best BLASTp score to the BFV Bel2 (E = 4E-08; maximum identity = 26%).
Discovery and characterization of PSFVaye and ChrEFV
Identification of foamy virus (FV)-like sequences in the aye-aye and Cape golden mole genome
FV genomic region
BLASTp best hit
Aye-aye ( Daubentonia madagascariensis )
(5′ truncated) 5′-LTR
Primer binding site
CAA70074 Gag FFV
AAF64414 Pol EFV
(3′ truncated) env
AFK85016 Env RhiFV
Cape golden mole ( Chrysochloris asiatica )
(5′ truncated) pol
<615.. > 3534
AFX98084 Pol SFV
<3464.. > 6371
AFR79245 Env BFV
Within PSFVaye, via sequence homology, we identified i) a PBS (TGACACCCAATGTGGGGC; nt 114–131) with a broad tRNA usage (tRNALys, Asn, Thr, Pro) predicted, ii) a complete putative gag gene (nt 198–1,971) which is interrupted by a SINE (AluJo , nt 1,449-1,748), iii) a complete putative pol gene (nt 1,941-5,727; PR: nt 1,968-2,459; RT: nt 2,415-3,306; RNase H: nt 3,985-4,467; IN: nt 4,480-5,676) which is also interrupted by SINEs (AluJ Mim  & ALU , nt 2,871-3,165), and iv) a 3′-truncated putative env gene (nt 5,612-7,867). We also found ~113 nt on the 5′ end of the PBS to be uniquely-mapping within the aye-aye genome. This might represent some portion of the 5′ end of the 5′ PSFVaye LTR (nt 1–111). We confirmed the truncation of the putative 5′-LTR and env gene and the absence of a 3′-LTR in the PSFVaye genome by using a genome walking of genomic DNA from two different aye-ayes. We found that PSFVaye exhibits several characteristic FV features. For example, its Gag protein contains a degenerate, but still recognizable, putative CTRS (GPR?VGD*WQRICLAFQY, aa 37–54) and three GR I-III boxes (box I: nt 1,248-1,358; box II: nt 1,422-1,781 interrupted by a SINE (nt 1,449-1,748); box III nt 1,851-1,871); the putative pol gene harbors a PPT (AAGGACAAAG, nt 5,092-5,101); and the N-terminus of the putative Env contains a highly conserved WXXW motif (WLAW, aa 10–13). A hydrophobic transmembrane region containing a signal-peptide peptidase cleavage site (ILIWIMLFLILFSAILVS↓TLIAVF) was also located in the N-terminus of the putative Env at aa 63–86. We could not identify a furin-cleavage site in the Env of PSFVaye, most likely due to its defective nature. However, by comparing the Env sequence of PSFVaye to those of other exogenous non-defective FVs, we speculate that it should be located between aa 113 and 127 (YFSQAHV*KSRAIHF). Moreover, its putative Gag and Env proteins only exhibit similarities to FV proteins (Table 1) and share no detectible similarities with other (non-FV) retroviral proteins. All these findings strongly suggest that the PSFVaye sequence is an endogenous FV, and not merely a distantly related FV-like ERV.
ChrEFV contains i) a 5′-truncated putative Pol coding gene (nt <1- < 2,937) which only contains full RT (nt 1–597), RNase H (nt 1,270-1,737), and IN coding domains (nt 1,750-2,931), ii) a full putative Env coding gene (nt <2,867- > 5,785) and iii) a full putative Bel1 coding gene (nt 5,752-6,393). Following the putative bel1 gene is a long non-repetitive sequence of ~1,700 nt containing a SINE (AFROSINE1B ). The bel2 gene and the 3′-LTR of ChrEFV may be located within this region but cannot be identified, most likely due to divergence from the probes we used. Nonetheless, the N-terminus of the ChrEFV Env protein harbors the highly conserved WXXW motif (WMRW, aa 10–13) and only shows significant similarity to FV Env proteins (Table 1). A hydrophobic transmembrane region harboring a putative signal-peptide peptidase cleavage site (VMTSYVS?LILLGIIITA↓SFI TIC, aa 74–97) was also found in the putative Env at the anticipated location, and a potential furin-cleavage site (GQIISNSSSNRRKI) could be located at aa 125–138 within the putative Env as well. Furthermore, the putative Bel1 protein also exhibits some similarity to FV Bel1 proteins which could be demonstrated by using the HMM searching technique (HMM profile GyDB (http://www.gydb.org/): Bel1 spumaretroviridae; E = 0.24; note that this weak support is rather typical and expected given that this is a part of an endogenous FV and that FV Bel1 protein is not conserved among FVs ). Although a number of FV-specific genomic features as outlined in the PSFVgal characterization section could not be located, possibly due to the lack of sequence data and/or divergence, this genomic information supports a FV progenitor of ChrEFV.
FV phylogenetic analyses
A retrovirus phylogeny (Additional file 3: Figure S2) was estimated using an alignment of RT protein sequences to establish the phylogenetic positions of PSFVgal, PSFVaye, and ChrEFV among retroviruses. All three new FV sequences are placed with sequences of extant FVs. Together with the genomic characterization above, the phylogeny strongly supports the classification of PSFVgal as a true exogenous FV and that both PSFVaye and ChrEFV are real endogenous FVs. To the best of our knowledge, no exogenous lemuriform and afrotherian FVs have been identified to date. Thus, our finding of PSFVaye, and ChrEFV raises the possibility of the presence of undetected FVs circulating among other lemurs and afrotherians.
Interestingly, while we can confirm that PSFVaye and ChrEFV are endogenous FVs, a number of ERVs found in fish genomes, including the genomes of coelacanth (CoeEFV) , zebra fish (DrFV-1) , platyfish (platyfishEFV) , and cod (CodEFV) , have been reported to show some detectible similarity to exogenous FVs, but more evidence is required to definitively determine that these ERVs truly have a FV origin. Since these fish ERVs are phylogenetically basal to all known extant mammalian FVs, it is difficult to say with certainty that the progenitors of these ERVs are “true” FVs, and not distinct lineages that branched earlier in retroviral evolution but which are “FV-like” enough to be detected by sequence similarity using the true FVs as probes. Moreover, in terms of genomic organization, these fish ERVs do not show all the features characteristic of mammalian FV genomes. For example, the two identified accessory genes of CoeEFV show no significant similarity to those of extant mammalian FVs , only two out of five determined DrFV-1 proteins, Gag and Pol, are FV-like , FV-like accessory genes in platyfishEFV have not been located , and only FV-like RT could be found for CodEFV . Although it is known that the accessory genes are the least conserved FV genes , and that this partial characterization of their genomes could be due to an ancient divergence of fish and mammalian FVs or simply the lack of sequence data, it raises a possibility that the progenitors of these ERVs may not be FVs but merely class III FV-like viruses. Resolution of this debate will require identification and analysis of FV genomes of other vertebrates like amphibians, reptiles, and birds, if they indeed exist . Phylogenetic analysis of extant fish FVs and ERVs would also help resolve this debate.
Despite ChrEFV having an inferred co-speciation event (posterior probability = 1), the ChrEFV branch lengths, including both the terminal branch and the internal branch leading to other FVs, are longer than would be expected based on the FV-host divergence linear relationship (Figure 3B; identified as outliers, green and blue dots, respectively, with Cook’s distance investigation). Given that ChrEFV is basal to all known extant FVs, this finding challenges the notion that ChrEFV is an ERV that has a FV origin. Nevertheless, ChrEFV is a long terminal branch in the inferred FV phylogeny and has accumulated numerous neutral genetic changes since its ancient endogenisation. Coupled with long-branch attraction, this ‘pseudogene effect’  could cause artificial inflation of the branch length with a concomitant deep placement of ChrEFV within the FV phylogeny. Thus, given the congruent ChrEFV/host phylogenies and the FV-like features of ChrEFV outlined above, it is still likely that ChrEFV is a genuine endogenous FV and not merely a class III FV-like ERV. Resolution of ChrEFV evolutionary history requires further analysis with other extant afrotherian, marsupial or monotreme FVs when and if they are identified. The results from these analyses could also elucidate the earlier events of mammalian FV diversification. To the best of our knowledge, this is the first (endogenous) afrotherian FV to be discovered, extending the FV host range to cover the whole eutherian clade.
PSFVgal was found to be a sister taxon of SFVs with robust support (posterior probability = 0.95), consistent with a FV-host co-speciation history. Additionally, we also identified a number of extant FVs from other lorises, including the silvery greater galago (Otolemur monteiri monteiri), southern lesser bush baby (Galago senegalensis moholi), and a potto (Perodicticus potto) (Additional file 1: Table S2), and partially sequenced their genomes (IN core domain: ~259 nt, 85 aa). When these new sequences were included in our phylogenetic analysis, they were all found to be distinct FVs forming a clade with PSFVgal with relatively high statistical support (posterior probability = 0.85) within which the branching orders generally mirror those of their lorisiforme hosts (Additional file 6: Figure S5). These results suggest that co-speciation between exogenous lorisiforme FVs and their hosts is very stable, extending to the species level. Together, our results indicate an ancient distribution and evolution of lorisiforme FVs in continental Africa dating to the divergence between strepsirrhine primates (lemurs and lorises) and the anthropoid primates about 70–95 Ma [43, 62, 63]. The absence of PSFV in some lorises (Additional file 1: Table S2) may be due to the testing of limited numbers of individuals and/or testing of captive animals that were captured as infants and thus were likely PSFV-negative at that time.
In contrast, PSFVaye was not a sister-taxon of PSFVgal as would be expected, showing a conflict in the co-evolutionary history of PSFVaye (Figure 3A). Instead, the sequence is robustly placed outside the boreoeutherian FV clade (posterior probability = 1), but still remains together with exafroplacentalian FVs (posterior probability = 1). The same evolutionary history was inferred for RhiFV (Figure 3A); instead of being grouped together with other fereungulata FVs (bovine, equine, and feline FVs) as would be expected if it were to co-speciate with its host, it is placed robustly outside the boreoeutherian FV clade (posterior probability = 1), but still remains together with exafroplacentalian FVs (posterior probability = 1). We compared our best estimated phylogeny against hypothetical alternative phylogenies where PSFVaye and/or RhiFV co-speciate with their hosts using approximately unbiased tests , and found that the co-speciation hypotheses of PSFVaye and RhiFV are both rejected (Additional file 4: Figure S3-A). Interestingly, our analysis also inferred a sister group relationship between PSFVaye and RhiFV although statistical support for this relationship is extremely weak (posterior probabilities = 0.56, Figure 3A).
Consistent with previous findings [37, 42], our analyses suggest an extremely stable FV-host co-speciation history across a timescale spanning millions of years. Interestingly, this is in contrast to what is suggested by various in vitro experiments. For example, an in vitro study has shown prototype FV (PFV) to be capable of infecting not only primate, but also rodent, laurasiatheria, avian and reptile cells, with an extremely broad range of susceptible cell types . Furthermore, in vitro investigations of FV cell-attachment and entry also suggested that FVs of different species likely utilize the same, perhaps promiscuous, receptor molecule(s) for cell-attachment and/or entry [69–71]. Although these findings might in part help to explain FV cross-species transmissions (FV speciation events by means of jumping from one host species to another, without the host speciating at the time of the jump) observed in nature (see below), one would expect to see more host-switches happening given these findings. Analyses of several anti-retroviral restriction factors have shown that they are specific to particular FVs. For example, while the tripartite motif protein 5αs (TRIM5αs) of most New World monkeys (NWMs) have been shown to be able to restrict some NWM FVs, PFV and a SFV from macaque, TRIM5αs from apes cannot, but instead can restrict feline FV [19, 72]. Inhibition of apolipoprotein B mRNA-editing, enzyme-catalytic, polypeptide-like 3C (APOBEC3C) by FVs has also been shown to be species-specific . Whether this species-specific FV-host antagonistic interaction is one of the factors underlining the stable co-speciation history or the result of it is still unclear. To date, the factors that determine this extremely stable pattern of FV-host co-speciation in nature are still very poorly understood.
A more extreme case of host switching across taxonomic orders may explain the phylogenetic placement of PSFVaye and RhiFV, where both are placed robustly outside the boreoeutherian FV clade (posterior probability = 1). Nevertheless, one possibility for the inferred placement of PSFVaye may be that it is an artefact due to neutral genetic changes that have accumulated since its endogenization . To examine whether or not neutral evolution alone could be responsible for the observation, Bayesian phylogenies of concatenated Pol-Env protein alignments containing artificially and uniformly mutated SFVagm, SFVsqu, and PSFVgal sequences were inferred. These sequences were mutated in silico with an overestimate of the mammalian neutral substitution rate of 10E-9 substitutions per site per year (s/n/y) over a mean estimate of PSFVaye endogenization date of 40 Myr ago (Figure 2). The substitution process was assumed to be neutral, homogenous, and independent across sites and base types. Re-aligning the sequences after the simulation using MUSCLE implemented in MEGA6  does not change the relative position of mutated sequences in other SFV sequences. Our results suggest that, even with this biased substitution rate setting of 10E-9 s/n/y, the neutral genetic change accumulation alone is insufficient to explain the inferred phylogenetic position of PSFVaye. Our in silico analysis only increased the branch length and decreased the branch support without changing the tree topology (Additional file 7: Figure S6). It is thus unlikely that the placement of PSFVaye is an artefact. Similarly, the placement of RhiFV is also likely to be genuine since it is an extant exogenous FV, supported by the absence of in-frame stop codons or frameshift mutations, and that it was found to be present in viral particles and not host genomic material . In addition, the RhiFV lineage remains in the same phylogenetic position even after removal of PSFVaye from the phylogenetic analysis.
We then asked whether there is a reasonable FV host switching model that can explain both (i) the deep phylogenetic placement of PSFVaye and RhiFV which is basal to all known exafroplacentalian FVs, and (ii) the dispersed geographical distribution of their hosts (Madagascar and China, respectively). We propose the following scenario to explain the phylogenetic relationships of PSFVaye and RhiFV within the context of the evolutionary and geographical history of FV hosts (Figure 4). Approximately 180–160 Ma, Gondwana was split into two separate landmasses: the Madagascar-India-Antarctica (MIA) landmass and the Africa-South America (ASA) landmass [78, 79]. The MIA landmass then split into the Madagascar-India and the Antarctica landmasses about 135 Ma [78, 79]. Thereafter, but still in the very early period of eutherian evolutionary history, eutherians arrived onto the Madagascar-India landmass upon which their FVs established stable populations (Figure 4A). Since then, these early eutherian FVs evolved independently from, and hence have distant evolutionary relationships with, FVs circulating on the African landmass, which later would give rise to all extant SFVs, lorisiforme FVs and fereungulata FVs. The India landmass was then split from the Madagascar landmass about 80–90 Ma [78, 79] (Figure 4B), separating the FV populations into two groups: one circulating in Madagascar which gave rise to the progenitor of PSFVaye and the second transferred to Laurasia via the India landmass continental drift which independently gave rise to RhiFV (Figure 4C). However, the model does not suggest which ancestral species introduced FV to the lemur and bat lineages and when these transmissions occurred. Nonetheless, both PSFVaye and RhiFV are placed robustly with exafroplacentalian FVs but outside the boreoeutherian FV clade (posterior probability = 1 and 1, respectively). This suggests that the original hosts were exafroplacentalians, and not boreoeutherians, and that the host switches occurred between species across deep clades long diverged from each other. Sequence data of FVs from additional mammals, especially bats and lemurs, are required to examine further the evolutionary histories of PSFVaye and RhiFV.
Our hypothesis helps to explain the deep phylogenetic placement and dispersed geographical distribution of PSFVaye and RhiFV. In addition, this hypothesis also gives an explicit prediction that PSFVaye and RhiFV should share a MRCA prior to 80–90 Ma based on the estimated Madagascar-India landmass split date [78, 79]. Indeed, although not well-supported (posterior probability = 0.56), our best phylogenetic estimate suggests that PSFVaye and RhiFV have a sister-group relationship (Figure 3A). Assuming that our inferred PSFVaye-RhiFV relationship is correct, we estimated the tMRCA of PSFVaye and RhiFV to be 93 (85–101) Myr based on the linear relationship of FV-host diversity (Figure 3B). Although slightly high, our estimate range largely overlaps with the date range derived from this hypothesis. The accumulated neutral changes in the PSFVaye sequences may explain both the weak branch support and the slightly overestimated divergence date.
Here, we report the characterization of the complete PSFVgal genome obtained from a galago, and describe in more detail the partial genome of PSFVaye from the aye-aye. The genomic organization of PSFVgal and PSFaye is characteristic of FVs, and they are phylogenetically placed within the FVs. The defective nature of the PSFVaye genome, coupled with its robust amplification from genomic DNA and an absence of antibodies to the genetically related PSFVgal in infected animals, indicates that it is endogenous. In contrast, PSFVgal has a complete viral transcriptome, elicits a FV-specific immune response, and is not present in all individuals, consistent with other exogenous FVs. We also report the discovery and describe a novel ERV present in the Cape golden mole genome, ChrEFV. Genomic inspection and analysis suggests that ChrEFV is also an endogenous FV, and its phylogenetic placement is consistent with a history of co-speciation between FVs and their mammalian hosts. However, ChrEFV has a high level of accumulated sequence divergence compared to other FVs, likely an artefact of relaxation of evolutionary constraints since becoming endogenous.
We found a general, stable pattern of mammalian FV-host co-divergence which extends as deep as the exafroplacentalian basal diversification, spanning more than 100 Myr. To date, it is still poorly understood why FVs stably co-speciate with their natural hosts. Furthermore, we also identified two possible cases of host switching in the evolutionary history of FVs. The first was observed in NWMs, as previously shown , which may have happened during captivity or in the wild and which has previously been reported to rarely occur in Old World monkeys and apes . However, others have shown congruent NWM and SFV phylogenies using a larger distribution of species-specific sequences . The second involves PSFVaye and RhiFV which may involve cross-species transmission at the level of mammalian orders. We propose a scenario based upon geographical knowledge of continental drift and hypothetical migration of ancient eutherians to explain this observation. Our results highlight the value of integrating multiple sources of information to elucidate the evolutionary history of mammals and their viruses, including continental and geographical histories, ancestral host locations, in addition to the natural history of host and virus.
Nonhuman primate (NHP) samples and nucleic acid preparation
NHP specimens (serum, frozen or fresh whole blood, tissues) were purchased from the Duke Lemur Center, were collected on an opportunistic basis from seven U.S. zoos following approved animal use protocols, or dried blood spots (DBS) via NHP hunters in Cameroon from freshly hunted NHP bushmeat in a study approved by Institutional Animal Care and Use Committees and the Cameroon government (Additional file 1: Table S2). Hunters were educated about the risks associated with direct contact with NHPs and about appropriate prevention measures. Prior to processing, all specimens were stored at -20°C or -80°C, except DBS which were stored in nytran ziplock bags with dessicant at room temperature. DNA lysates or nucleic acids were prepared from blood specimens and DBS as previously described . Nucleic acids were extracted from archived tissue specimens from the Duke Lemur Center (liver, spleen, kidney) using a Qiagen QIAmp DNA mini kit. DNA concentrations were determined using a Nanodrop spectrophotometer and DNA integrity was confirmed with ß-actin PCR.
Isolation of PSFVgal
HeLa cells were infected with PSFVgal (SFV-5; ATCC# VR-644), originally isolated from the throat swab of an Otolemur crassicaudatus panganiensis, and grown in complete DMEM medium supplemented with 10% FBS, 2 mM L-glutamine, 100 ug/ml streptomycin and 100 U/ml penicillin (Invitrogen). PSFVgal is the only viral isolate analyzed in our study. Cell cultures were incubated at 37°C and 5% CO2 and were split when 90% confluent. Cellular DNA was prepared using a Qiagen kit and quantified with a Nanodrop spectrophotometer when cytopathic effect was observed in >50% of the cells.
Western blot (WB) assay
PSFVgal-infected and uninfected cells were pelleted by centrifugation at 1,500 rpm for 10 min and washed 2X with phosphate-buffered saline. Antigen for WB testing was prepared by treating cell pellets with WB sample buffer (Invitrogen) at 100°C for 10 min. Protein concentrations were determined using the BioRad DC Protein Assay kit and 150 ug of infected and uninfected HeLa cell lysates were applied separately to 4-12% polyacrylamide gels (Invitrogen) for WB analysis. Serum specimens were tested using a 1:50 dilution as previously described .
PCR-amplification of 5′ and 3′ PSFVgal genome halves, plasmid cloning, and sequence analysis
To obtain the full-length genomic sequence of PSFVgal we first PCR-amplified small regions in the LTR and polymerase (pol) gene by using nested PCR with degenerate SFV primers and conditions provided elsewhere [80, 81]. PSFVgal-specific primers were then designed from the LTR and pol fragments to amplify two overlapping genomic halves using PCR and a Roche Expand 20 kb PCR System Kit following the manufacturer’s instructions. The Expand kit contains thermostable Taq DNA polymerase and a thermostable DNA polymerase with proofreading activity to minimize PCR-induced sequence artefacts. Briefly, the PCR reaction mixture was prepared in 50 ul reaction volumes containing 500 uM of each dNTP, 400 nM of each forward and reverse primer, and 5 units of Expand 20 kb enzyme mix using 0.5 ug of PSFVgal tissue culture DNA. The LTR-pol fragment was amplified using the primers PSFVgal-LF1 5′ GGC TTG GAT AAT TAA TTG TTA GAT GCT CTG 3′ and PSFVgalpolF2R 5′ GTT CCA AAC GTA TGC CCC TCT CCT T 3′. The PSFVgal pol-LTR fragment primers are PSFVgalpolF1 5′ GTC AGC ATT CAC CTC TTC CAC CTT G 3′ and PSFVGAL-LR1 5′ GAC TTA TTT ATT ACT GCA AGA CCC GAG AGG G 3′. Following denaturation at 92°C for 2 min, amplification consisted of 10 cycles at 92°C for 10 secs, 60°C for 30 secs, and 68°C for 6 mins, followed by 20 cycles at 92°C for 10 secs, 60°C for 30 secs, and 68°C for 6 mins with an additional 10 sec per cycle.
PCR products were visualized with 0.8% agarose gel electrophoresis and amplicons of the expected size were collected using a QIAGEN gel purification kit. Purified amplicons were cloned into the pCR-XL-TOPO vector (Invitrogen) for genome sequencing. Two plasmids were obtained; pCR-XL-TOPO- PSFVgal-LP containing the LTR-pol insert (5,148-bp) and pCR-XL-TOPO- PSFVgal-PL containing the pol-LTR insert (6,201-bp). The PSFVgal plasmid DNAs were purified using the Purelink Hipure Plasmid Filter Purification Kit (Invitrogen) and sequenced in both directions using a 3130XL Genetic Analyzer (Applied Biosystems). Complete genomes were assembled using the software program Vector NTI v11.1. 5′ and 3′ LTRs were determined based on overlapping sequences and comparison with available FV genome sequences.
Homology of PBS sequences to tRNAs was inferred using the transfer RNA database (http://trnadb.bioinf.uni-leipzig.de/). Confirmation of the tRNA PBS sequence was done by PCR amplification with primers flanking that region and using 0.1 and 1.0 ug of PSFVgal-infected HeLa tissue culture DNA. Primary and nested primers were used to independently amplify two products of different lengths and amplification was done separate from that to obtain the 5′ halve of the genome containing the PBS. Primary PCR primer sequences are F982 (5′-GAC CAG TGT GAG ATT GGT GTC TC-3′) and R1547 (5′-CCA GTT GCC TCC TGT GAT TCT AAC-3′) and the internal primers are F1133 (5′-CTC CTG GTT GAG GAC AAG GGA AC-3′) and R1475 (5′-CCA ATT TGT GCT CGT GGC ACT GG-3′) and standard PCR conditions were used. Numbers in the primer names reference their locations in the PSFVgal genome. PCR products from both primary and nested and 0.1 and 1.0 ug DNA inputs were sequenced.
Discovery and characterization of ChrEFV and PSFVaye
Publically available GenBank whole genome shotgun (WGS) sequences were screened using tBLASTn and various FV protein sequences (Additional file 1: Table S3). As previously reported by Han and Worobey , several matches were returned from the WGS assembled sequence of the aye-aye, Daubentonia madagascariensis, spanning four non-overlapping contigs. The adjacency of these four contigs and sequences flanking the gag and env regions, especially the LTR and host integration sequences, were determined using a Seegene genome walking kit (Seegene, MD, USA). Three rounds of PCR were performed using universal primers provided in the kit (DW2-ACP, DW2-ACPN, UNIP2) with PSFVaye-specific primers (1st PCR primers: upstream PSFVAYE-1174R/DW2-ACP; downstream, PSFVAYE-6122 F/DW2-ACP; 2nd PCR primers: upstream, PSAGF2REV/DW2-ACPN; downstream, PSAER2FOR/DW2-ACPN; 3rd PCR primers: upstream, PSAGF1REV/UniP2; downstream, PSAER1FOR/UniP2). Positive bands were cloned and Sanger sequenced using an ABI7700 instrument. Furthermore, this in silico screening process also returned one contig from the WGS assembled sequence of a Cape golden mole (Chrysochloris asiatica), designated ‘ChrEFV’ here. Both were annotated via sequence homology. The presence of simple repetitive elements was determined by CENSOR (http://www.girinst.org/censor/).
PCR testing of Prosimian genomic DNA specimens for FV
1 ug of prosimian DNA was PCR tested for FV pol sequences using a combination of nested primer sets or a single round of PCR to detect PSFVaye gag sequences depending on the species (Additional file 1: Table S2). The first PCR assay is PSFVgal-specific and is based on the PSFVgal genome and uses the first and second round PCR primers SIF1GAL 5′ CTT GCT GTG CAG AGC AGT CAC AAG GT 3′ and SIR1GAL 5′ GTT TTA TTT CAC TGT TTT TCC GTT CCA C 3′ and SIF3GAL 5′ CCA AGT CTG GAT GCA GAG CTT ATC CA 3′ and SIR3GAL 5′ ACT TTG GGG GTG ATA CGG AGT ACT 3′ to generate 712-bp and 635-bp fragments, respectively. The second assay was designed using an alignment of complete Old World and New World primate SFV genomes and uses the primary primers SIF5N 5′ TAC ATG GTT ATA CCC CAC KAA GGC TCC TCC 3′ and SIR5N 5′ AAT AAW GGA TAC CAC TTT GTA GGT CTT CC 3′ and nested primers SIP4N 5′ TGC ATT CCG ATC AAG GAT CAG CAT T 3′ and SIR1NN 5′ GTT TTA TYT CCY TGT TTT TCC TYT CCA CCA T 3′ to generate to generate 282-bp and 141-bp pol sequences, respectively. The third assay was designed using an alignment of PSFVgal and PSFVaye genomes and uses the primary primers 3′ FVPF05 5′ KKM TAY TGG TGR CCT AAT ATG 3′ and FVPFR5 5′ GGT CTW CCA ACY ART AGT TTA G 3′ and nested primers FVPF01 5′ CCT TTT GAT AAA ATY TAT ATG G 3′ and FVPR01 5′ CAS CTT TCC ACT ACT TTG G 3′ to generate 483-bp and 292-bp products, respectively. To detect PSFVaye gag sequences the primers PSAGF1 5′ AAG ACC CTT GCT GCC TAA TGT TGG 3′ and PSAGR1 5′ TAT TTG TAA CCA GGG CTT GAC CAG 3′ were used to amplify a 475-bp sequence. 5 ul of primary PCR product was used as template in the nested PCR reaction. Selected PCR products were visualized by agarose gel electrophoresis analysis, extracted using a QIAquick Gel Extraction kit, and sequenced in both directions using an ABI7700 instrument.
Estimating the integration dates of PSFVaye and ChrEFV with Monte Carlo simulation
To estimate the age of PSFVaye and ChrEFV, we analyzed the frequencies of in-frame stop-codons in their Pol protein sequences. Excluding a hypothetical last stop codon, PSFVaye (1,163 codons) and ChrEFV (975 codons) Pols contain 9 and 14 in-frame stop codon mutations, which are equivalent to stop codon frequencies of 0.00774 and 0.0144, respectively. Using a Monte Carlo simulation approach and pols of PSFVgal, SFVcpz, SFVagm, SFVsqu, FFV, and SloEFV as model sequences, we estimated two cumulative distribution functions (CDFs) of time duration for in-frame stop codons to be accumulated at these two frequencies for each model sequence. Each pol model sequence was mutated in silico over various hypothetical time durations from 5 Myr to 200 Myr with an increment of 5 Myr. The mutation (substitution) process was assumed to be neutral, homogenous, and independent across all sites and base types with the rate of 2.2E-9 s/n/y, an average substitution rate of mammalian genomes . 1,000 simulations were performed for each period of time and the frequency of in-frame stop codons was then calculated for each simulation. Two separate CDFs of time duration were constructed for each sequence model using stop-codon frequencies of 0.00774 and 0.0144. In total, we constructed six CDFs for each of the two stop-codon frequencies. Means, medians, modes and 95% confidence intervals of the time duration for in-frame stop codons to accumulate at the two frequencies were calculated.
To investigate the phylogenetic relationships of PSFVgal, PSFVaye, and ChrEFV with other retroviruses, a consensus unrooted phylogenetic tree was constructed using a manually curated alignment of 240 RT protein sequences (162 aa) from alpha-, beta-, gamma-, delta-, epsilon-, lenti-, and spuma-like retroviruses using RAxML 7.2.8-HPC2 on XSEDE . The alignment is available from the authors upon request. The LG + G + F model was determined to fit the data best by using ProtTest 2.4 . Bootstrap support values for the branching order were calculated using 5,000 pseudoreplicates.
To further investigate the phylogenetic relationships of PSFVgal, PSFVaye, and ChrEFV with other FVs, a Bayesian FV phylogeny was estimated using MrBayes 3.2.1  and a manually-curated alignment of concatenated Pol-Env sequences (893 and 603 amino acids, respectively). We used CoeEFV  to root our mammalian FV tree since it is the most immediate outgroup of mammalian FVs known to date although it is still debatable whether this ERV is a real endogenous FV or not. Gag sequences were not included in the alignment since (i) the Gag of RhiFV is not available and (ii) the alignable region of Gag for the rest of the FVs is relatively short (~180 aa). Potential recombination among the sequences within the alignment were checked using VisRD3  both at nucleotide and protein levels. For both analyses, the null distribution was built based on 1,000 sets of randomly-shuffled sequences and the extended statistical geometry (Hamming) weighting option was selected. In the nucleotide analysis, the window size and step size was 300 and 60 nt, respectively; while in the protein analysis, the window size and step size was 100 and 20 aa, respectively. The results showed no significant evidence for recombination. The alignment is available from the authors upon request. The protein phylogeny reconstruction was performed using best available protein-specific model partitions of rtREV + Γ(4) + F (Pol) and WAG + I + Γ(4) + F (Env) as determined by ProtTest 2.4  under the AIC criterion. The MCMC was run for 5,000,000 steps with the initial 25% discarded as burn-in. Trees and parameters were logged every 100 steps. Parameter value convergence and sampling independency were manually inspected and had effective sample sizes >1,000.
We evaluated the FV-host co-speciation hypothesis using two analyses: (i) FV-host phylogenetic reconciliation analysis and (ii) FV-host divergence correlation analysis. Phylogenetic reconciliation analysis of FV and host trees was performed using Jane v4.0  (Genetic algorithm: number of generations = 100, population size = 100). The vertex-based cost mode was used, and the costs were set as follows: co-speciation = -1, duplication = 0, duplication & host switch = 0, loss = 0, and failure to diverge = 0. This setting was adopted in order to maximize the number of co-speciation events. The probability of observing the inferred number of potential co-speciation events by chance was calculated by random tip mapping in Jane v4.0  (Genetic algorithm: number of generation = 100, population size = 100, sample size = 1,000). This method of analyzing phylogenetic incongruence cannot formally differentially weight transmission at different taxonomic levels however, and the derivation of p-values relies upon the interchangeability of branches. Thus, it does not discriminate between distant and close interspecies transmissions. Once the potential FV-host co-speciation events were located, we then identified FV-host corresponding branches that represent FV-host co-evolutionary histories for the FV-host divergence correlation analysis. In this test, we investigated whether or not the lengths of these corresponding branches were linearly correlated, i.e. whether or not the accumulation of FV sequence divergence occurred in proportion to host divergence, using linear regression implemented in MATLAB . Outliers were determined by using Cook’s distance inspection. Outliers were removed one at a time using a threshold value of 3x mean of the observed Cook’s distances and the model was re-fitted until no outliers were found. The model was then used to estimate how long the parental branch of RhiFV and PSFVaye branches (0.0564 amino-acid substitutions per site) is in units of time. Given that the radiation of exafroplacentalians was 101.1 Ma , we then re-calculated when PSFVaye and RhiFV shared a MRCA to be 93.342 (85.095-101.589) Ma.
Our inferred phylogeny suggests two possible cases of FV host-switches: one involves PSFVaye and RhiFV and the other involves NWM FVs. Approximate unbiased (AU) tests using multi-scale bootstrapping  were performed to investigate whether or not our Pol-Env/Pol/Env alignments reject the co-speciation model for PSFVaye and RhiFV, and NWM FVs (Additional file 4: Figure S3). Four completing alternative phylogenetic placements of PSFVaye and RhiFV (Additional file 4: Figure S3-A), and two alternative placements of NWM FVs (Additional file 4: Figure S3-B) are compared to each other. The site-wise log likelihoods used in the tests were computed using PAML 4.7a . The best available amino acid substitution model, as determined by ProtTest 2.4  under the AIC criterion, was used (Pol: rtREV + Γ(4) + F; Env: WAG + Γ(4) + F). A molecular clock was not imposed, and ambiguous sites were included. AU tests were performed in Consel  with default settings.
Additional Bayesian FV phylogenies were constructed using the same settings as described above. One included a number of short fragments of other lorisiforme FV Pol sequences (IN core domain sequences, 85 aa) to further investigate phylogenetic relationships of lorisiforme FVs with other FVs. The additional Bayesian phylogenies included several in silico mutated sequences to examine the effects of neutral genetic changes on the inferred phylogenetic relationships. The alignment was not re-aligned after the sequence simulation (see main text). Both alignments are available from the authors upon request.
We thank the seven U.S. zoos, including the San Diego Zoo, Zoo New England, the Gladys Porter Zoo, the Baltimore Zoo, the Detroit Zoo, and two other zoos that wish to remain anonymous, and the Duke Lemur Center for providing us with specimens for the study. We also thank David Sintasath and Michelle Sturgeon for help preparing nucleic acids from the potto DBS. Funding for this study was provided in part by the Faucett Family and Rawji Foundations and Global Viral Forecasting (ANRS12182/12255 and NIH RO1 AI50529), graciously supported by the U.S. Department of Defense Armed Forces Health Surveillance Center, Division of Global Emerging Infections, Surveillance Operations (AFHSC GEIS), the Defense Threat Reduction Agency Cooperative Biological Engagement Program (DTRA-CBEP), Department of Defense HIV/AIDS Prevention Program (DHAPP), Google.org, the Skoll Foundation, and the U.S. Agency for International Development (USAID) Emerging Pandemic Threats Program, PREDICT project, under the terms of Cooperative Agreement Number GHN-A-OO-09-00010-00. Use of trade names is for identification only and does not imply endorsement by the U.S. Department of Health and Human Services, the Public Health Service, or the Centers for Disease Control and Prevention. The findings and conclusions in this report are those of the authors and do not necessarily represent the views of the Centers for Disease Control and Prevention. P.A. is funded by the Royal Thai Government. AK is funded by the Royal Society. We thank R. Blakey for tectonic maps. This is DLC publication number 1273.
- Linial ML: Foamy viruses are unconventional retroviruses. J Virol. 1999, 73: 1747-1755.PubMed CentralPubMedGoogle Scholar
- Enders JF, Peebles TC: Propagation in tissue cultures of cytopathogenic agents from patients with measles. Proc Soc Exp Biol Med. 1954, 86: 277-286.View ArticlePubMedGoogle Scholar
- Rustigian R, Johnston P, Reihart H: Infection of monkey kidney tissue cultures with virus-like agents. Proc Soc Exp Biol Med. 1955, 88: 8-16.View ArticlePubMedGoogle Scholar
- Johnston PB: A second immunologic type of simian foamy virus: monkey throat infections and unmasking by both types. J Infect Dis. 1961, 109: 1-9.View ArticlePubMedGoogle Scholar
- Stiles GE, Bittle JL, Cabasso VJ: Comparison of simian foamy virus strains including a new serological type. Nature. 1964, 201: 1350-1351.View ArticlePubMedGoogle Scholar
- Johnston PB: Taxonomic features of seven serotypes of simian and ape foamy viruses. Infect Immun. 1971, 3: 793-799.PubMed CentralPubMedGoogle Scholar
- Rogers NG, Basnight M, Gibbs CJ, Gajdusek DC: Latent viruses in chimpanzees with experimental kuru. Nature. 1967, 216: 446-449.View ArticlePubMedGoogle Scholar
- Hooks JJ, Gibbs CJ, Chou S, Howk R, Lewis M, Gajdusek DC: Isolation of a new simian foamy virus from a spider monkey brain culture. Infect Immun. 1973, 8: 804-813.PubMed CentralPubMedGoogle Scholar
- Hooks JJ, Gibbs CJ: The foamy viruses. Bacteriol Rev. 1975, 39: 169-185.PubMed CentralPubMedGoogle Scholar
- Heberling RL, Kalter SS: Isolation of foamy viruses from baboon (Papio cynocephalus) tissues. Am J Epidemiol. 1975, 102: 25-29.PubMedGoogle Scholar
- McClure MO, Bieniasz PD, Schulz TF, Chrystie IL, Simpson G, Aguzzi A, Hoad JG, Cunningham A, Kirkwood J, Weiss RA: Isolation of a new foamy retrovirus from orangutans. J Virol. 1994, 68: 7124-7130.PubMed CentralPubMedGoogle Scholar
- Bieniasz PD, Rethwilm A, Pitman R, Daniel MD, Chrystie I, McClure MO: A comparative study of higher primate foamy viruses, including a new virus from a gorilla. Virology. 1995, 207: 217-228.View ArticlePubMedGoogle Scholar
- Marczynska B, Jones CJ, Wolfe LG: Syncytium-forming virus of common marmosets (Callithrix jacchus jacchus). Infect Immun. 1981, 31: 1261-1269.PubMed CentralPubMedGoogle Scholar
- Riggs JL, Oshirls , Taylor DO, Lennette EH: Syncytium-forming agent isolated from domestic cats. Nature. 1969, 222: 1190-1191.View ArticlePubMedGoogle Scholar
- Malmquist WA, Van der Maaten MJ, Boothe AD: Isolation, immunodiffusion, immunofluorescence, and electron microscopy of a syncytial virus of lymphosarcomatous and apparently normal cattle. Cancer Res. 1969, 29: 188-200.PubMedGoogle Scholar
- Tobaly-Tapiero J, Bittoun P, Neves M, Guillemin MC, Lecellier CH, Puvion-Dutilleul F, Gicquel B, Zientara S, Giron ML, de Thé H, Saïb A: Isolation and characterization of an equine foamy virus. J Virol. 2000, 74: 4064-4073.PubMed CentralView ArticlePubMedGoogle Scholar
- Kupiec JJ, Kay A, Hayat M, Ravier R, Périès J, Galibert F: Sequence analysis of the simian foamy virus type 1 genome. Gene. 1991, 101: 185-194.View ArticlePubMedGoogle Scholar
- Renne R, Friedl E, Schweizer M, Fleps U, Turek R, Neumann-Haefelin D: Genomic organization and expression of simian foamy virus type 3 (SFV-3). Virology. 1992, 186: 597-608.View ArticlePubMedGoogle Scholar
- Pacheco B, Finzi A, McGee-Estrada K, Sodroski J: Species-specific inhibition of foamy viruses from South American monkeys by New World Monkey TRIM5α proteins. J Virol. 2010, 84: 4095-4099.PubMed CentralView ArticlePubMedGoogle Scholar
- Herchenröder O, Renne R, Loncar D, Cobb EK, Murthy KK, Schneider J, Mergia A, Luciw PA: Isolation, cloning, and sequencing of simian foamy viruses from chimpanzees (SFVcpz): high homology to human foamy virus (HFV). Virology. 1994, 201: 187-199.View ArticlePubMedGoogle Scholar
- Thümer L, Rethwilm A, Holmes EC, Bodem J: The complete nucleotide sequence of a New World simian foamy virus. Virology. 2007, 369: 191-197.View ArticlePubMedGoogle Scholar
- Verschoor EJ, Langenhuijzen S, van den Engel S, Niphuis H, Warren KS, Heeney JL: Structural and evolutionary analysis of an orangutan foamy virus. J Virol. 2003, 77: 8584-8587.PubMed CentralView ArticlePubMedGoogle Scholar
- Schulze A, Lemey P, Schubert J, McClure MO, Rethwilm A, Bodem J: Complete nucleotide sequence and evolutionary analysis of a gorilla foamy virus. J Gen Virol. 2011, 92: 582-586.View ArticlePubMedGoogle Scholar
- Winkler I, Bodem J, Haas L, Zemba M, Delius H, Flower R, Flügel RM, Löchelt M: Characterization of the genome of feline foamy virus and its proteins shows distinct features different from those of primate spumaviruses. J Virol. 1997, 71: 6727-6741.PubMed CentralPubMedGoogle Scholar
- Renshaw RW, Casey JW: Transcriptional mapping of the 3′ end of the bovine syncytial virus genome. J Virol. 1994, 68: 1021-1028.PubMed CentralPubMedGoogle Scholar
- Wu Z, Ren X, Yang L, Hu Y, Yang J, He G, Zhang J, Dong J, Sun L, Du J, Liu L, Xue Y, Wang J, Yang F, Zhang S, Jin Q: Virome analysis for identification of novel mammalian viruses in bat species from Chinese provinces. J Virol. 2012, 86: 10999-11012.PubMed CentralView ArticlePubMedGoogle Scholar
- Achong BG, Mansell PW, Epstein MA, Clifford P: An unusual virus in cultures from a human nasopharyngeal carcinoma. J Natl Cancer Inst. 1971, 46: 299-307.PubMedGoogle Scholar
- Switzer WM, Bhullar V, Shanmugam V, Cong M-E, Parekh B, Lerche NW, Yee JL, Ely JJ, Boneva R, Chapman LE, Folks TM, Heneine W: Frequent simian foamy virus infection in persons occupationally exposed to nonhuman primates. J Virol. 2004, 78: 2780-2789.PubMed CentralView ArticlePubMedGoogle Scholar
- Meiering CD, Linial ML: Historical perspective of foamy virus epidemiology and infection. Clin Microbiol Rev. 2001, 14: 165-176.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu W, Worobey M, Li Y, Keele BF, Bibollet-Ruche F, Guo Y, Goepfert PA, Santiago ML, Ndjango J-BN, Neel C, Clifford SL, Sanz C, Kamenya S, Wilson ML, Pusey AE, Gross-Camp N, Boesch C, Smith V, Zamma K, Huffman MA, Mitani JC, Watts DP, Peeters M, Shaw GM, Switzer WM, Sharp PM, Hahn BH: Molecular ecology and natural history of simian foamy virus infection in wild-living chimpanzees. PLoS Pathog. 2008, 4: e1000097-PubMed CentralView ArticlePubMedGoogle Scholar
- Nemo GJ, Brown PW, Gibbs CJ, Gajdusek DC: Antigenic relationship of human foamy virus to the simian foamy viruses. Infect Immun. 1978, 20: 69-72.PubMed CentralPubMedGoogle Scholar
- Brown P, Nemo G, Gajdusek DC: Human foamy virus: further characterization, seroepidemiology, and relationship to chimpanzee foamy viruses. J Infect Dis. 1978, 137: 421-427.View ArticlePubMedGoogle Scholar
- Callahan ME, Switzer WM, Matthews AL, Roberts BD, Heneine W, Folks TM, Sandstrom PA: Persistent zoonotic infection of a human with simian foamy virus in the absence of an intact orf-2 accessory gene. J Virol. 1999, 73: 9619-9624.PubMed CentralPubMedGoogle Scholar
- CDC: Nonhuman primate spumavirus infections among persons with occupational exposure--United States, 1996. MMWR Morb Mortal Wkly Rep. 1997, 46: 129-131.Google Scholar
- Schweizer M, Falcone V, Gänge J, Turek R, Neumann-Haefelin D: Simian foamy virus isolated from an accidentally infected human individual. J Virol. 1997, 71: 4821-4824.PubMed CentralPubMedGoogle Scholar
- Heneine W, Switzer WM, Sandstrom P, Brown J, Vedapuri S, Schable CA, Khan AS, Lerche NW, Schweizer M, Neumann-Haefelin D, Chapman LE, Folks TM: Identification of a human population infected with simian foamy viruses. Nat Med. 1998, 4: 403-407.View ArticlePubMedGoogle Scholar
- Katzourakis A, Gifford RJ, Tristem M, Gilbert MTP, Pybus OG: Macroevolution of complex retroviruses. Science. 2009, 325: 1512-View ArticlePubMedGoogle Scholar
- Han G-Z, Worobey M: An endogenous foamy virus in the aye-aye (Daubentonia madagascariensis). J Virol. 2012, 86: 7696-7698.PubMed CentralView ArticlePubMedGoogle Scholar
- Han G-Z, Worobey M: An endogenous foamy-like viral element in the coelacanth genome. PLoS Pathog. 2012, 8: e1002790-PubMed CentralView ArticlePubMedGoogle Scholar
- Llorens C, Muñoz-Pomer A, Bernad L, Botella H, Moya A: Network dynamics of eukaryotic LTR retroelements beyond phylogenetic trees. Biol Direct. 2009, 4: 41-PubMed CentralView ArticlePubMedGoogle Scholar
- Schartl M, Walter RB, Shen Y, Garcia T, Catchen J, Amores A, Braasch I, Chalopin D, Volff J-N, Lesch K-P, Bisazza A, Minx P, Hillier L, Wilson RK, Fuerstenberg S, Boore J, Searle S, Postlethwait JH, Warren WC: The genome of the platyfish, Xiphophorus maculatus, provides insights into evolutionary adaptation and several complex traits. Nat Genet. 2013, 45: 567-572.View ArticlePubMedGoogle Scholar
- Switzer WM, Salemi M, Shanmugam V, Gao F, Cong M-E, Kuiken C, Bhullar V, Beer BE, Vallet D, Gautier-Hion A, Tooze Z, Villinger F, Holmes EC, Heneine W: Ancient co-speciation of simian foamy viruses and primates. Nature. 2005, 434: 376-380.View ArticlePubMedGoogle Scholar
- Bininda-Emonds ORP, Cardillo M, Jones KE, MacPhee RDE, Beck RMD, Grenyer R, Price SA, Vos RA, Gittleman JL, Purvis A: The delayed rise of present-day mammals. Nature. 2007, 446: 507-512.View ArticlePubMedGoogle Scholar
- Charleston MA, Robertson DL: Preferential host switching by primate lentiviruses can account for phylogenetic similarity with the primate phylogeny. Syst Biol. 2002, 51: 528-535.View ArticlePubMedGoogle Scholar
- Kupiec JJ, Sonigo P: Reverse transcriptase jumps and gaps. J Gen Virol. 1996, 77 (Pt 9): 1987-1991.View ArticlePubMedGoogle Scholar
- Eastman SW, Linial ML: Identification of a conserved residue of foamy virus Gag required for intracellular capsid assembly. J Virol. 2001, 75: 6857-6864.PubMed CentralView ArticlePubMedGoogle Scholar
- Patton GS, Morris SA, Chung W, Bieniasz PD, McClure MO: Identification of domains in gag important for prototypic foamy virus egress. J Virol. 2005, 79: 6392-6399.PubMed CentralView ArticlePubMedGoogle Scholar
- Stange A, Mannigel I, Peters K, Heinkelein M, Stanke N, Cartellieri M, Göttlinger H, Rethwilm A, Zentgraf H, Lindemann D: Characterization of prototype foamy virus gag late assembly domain motifs and their role in particle egress and infectivity. J Virol. 2005, 79: 5466-5476.PubMed CentralView ArticlePubMedGoogle Scholar
- Schliephake AW, Rethwilm A: Nuclear localization of foamy virus Gag precursor protein. J Virol. 1994, 68: 4946-4954.PubMed CentralPubMedGoogle Scholar
- Lindemann D, Pietschmann T, Picard-Maureau M, Berg A, Heinkelein M, Thurow J, Knaus P, Zentgraf H, Rethwilm A: A particle-associated glycoprotein signal peptide essential for virus maturation and infectivity. J Virol. 2001, 75: 5762-5771.PubMed CentralView ArticlePubMedGoogle Scholar
- Wilk T, Geiselhart V, Frech M, Fuller SD, Flügel RM, Löchelt M: Specific interaction of a novel foamy virus Env leader protein with the N-terminal Gag domain. J Virol. 2001, 75: 7995-8007.PubMed CentralView ArticlePubMedGoogle Scholar
- Geiselhart V, Bastone P, Kempf T, Schnölzer M, Löchelt M: Furin-mediated cleavage of the feline foamy virus Env leader protein. J Virol. 2004, 78: 13573-13581.PubMed CentralView ArticlePubMedGoogle Scholar
- Perry GH, Reeves D, Melsted P, Ratan A, Miller W, Michelini K, Louis EE, Pritchard JK, Mason CE, Gilad Y: A genome sequence resource for the aye-aye (Daubentonia madagascariensis), a nocturnal lemur from Madagascar. Genome Biol Evol. 2012, 4: 126-135.PubMed CentralView ArticlePubMedGoogle Scholar
- The draft genome of Chrysochloris asiatica. http://www.ncbi.nlm.nih.gov/nuccore/AMDV00000000.1,
- Jurka J: Origin and evolution of Alu repetitive elements. Mol Biol Intell Unit impact short interspersed Elem host genome. Edited by: Maraia RJ. 1995, Austin: R.G. Landes Company, 25-41.Google Scholar
- Jurka J: SINE elements from the mouse lemur. Repbase Reports. 2010, 10: 2171-Google Scholar
- Jurka J, Smith T: A fundamental division in the Alu family of repeated sequences. Proc Natl Acad Sci U S A. 1988, 85: 4775-4778.PubMed CentralView ArticlePubMedGoogle Scholar
- Jurka J: SINE elements from African elephants. Repbase Reports. 2010, 10: 491-Google Scholar
- Bodem J: Regulation of foamy viral transcription and RNA export. Adv Virus Res. 2011, 81: 1-31.View ArticlePubMedGoogle Scholar
- Perelman P, Johnson WE, Roos C, Seuánez HN, Horvath JE, Moreira MAM, Kessing B, Pontius J, Roelke M, Rumpler Y, Schneider MPC, Silva A, O’Brien SJ, Pecon-Slattery J: A molecular phylogeny of living primates. PLoS Genet. 2011, 7: e1001342-PubMed CentralView ArticlePubMedGoogle Scholar
- Horvath JE, Weisrock DW, Embry SL, Fiorentino I, Balhoff JP, Kappeler P, Wray GA, Willard HF, Yoder AD: Development and application of a phylogenomic toolkit: resolving the evolutionary history of Madagascar’s lemurs. Genome Res. 2008, 18: 489-499.PubMed CentralView ArticlePubMedGoogle Scholar
- Yoder AD, Yang Z: Divergence dates for Malagasy lemurs estimated from multiple gene loci: geological and evolutionary context. Mol Ecol. 2004, 13: 757-773.View ArticlePubMedGoogle Scholar
- Springer MS, Murphy WJ, Eizirik E, O’Brien SJ: Placental mammal diversification and the Cretaceous-Tertiary boundary. Proc Natl Acad Sci U S A. 2003, 100: 1056-1061.PubMed CentralView ArticlePubMedGoogle Scholar
- Lemey P, Lott M, Martin DP, Moulton V: Identifying recombinants in human and primate immunodeficiency virus sequence alignments using quartet scanning. BMC Bioinformatics. 2009, 10: 126-PubMed CentralView ArticlePubMedGoogle Scholar
- Conow C, Fielder D, Ovadia Y, Libeskind-Hadas R: Jane: a new tool for the cophylogeny reconstruction problem. Algorithms Mol Biol. 2010, 5: 16-PubMed CentralView ArticlePubMedGoogle Scholar
- McAllister BF, Werren JH: Phylogenetic analysis of a retrotransposon with implications for strong evolutionary constraints on reverse transcriptase. Mol Biol Evol. 1997, 14: 69-80.View ArticlePubMedGoogle Scholar
- Shimodaira H: An approximately unbiased test of phylogenetic tree selection. Syst Biol. 2002, 51: 492-508.View ArticlePubMedGoogle Scholar
- Hill CL, Bieniasz PD, McClure MO: Properties of human foamy virus relevant to its development as a vector for gene therapy. J Gen Virol. 1999, 80 (Pt 8): 2003-2009.View ArticlePubMedGoogle Scholar
- Berg A, Pietschmann T, Rethwilm A, Lindemann D: Determinants of foamy virus envelope glycoprotein mediated resistance to superinfection. Virology. 2003, 314: 243-252.View ArticlePubMedGoogle Scholar
- Plochmann K, Horn A, Gschmack E, Armbruster N, Krieg J, Wiktorowicz T, Weber C, Stirnnagel K, Lindemann D, Rethwilm A, Scheller C: Heparan sulfate is an attachment factor for foamy virus entry. J Virol. 2012, 86: 10028-10035.PubMed CentralView ArticlePubMedGoogle Scholar
- Herchenröder O, Moosmayer D, Bock M, Pietschmann T, Rethwilm A, Bieniasz PD, McClure MO, Weis R, Schneider J: Specific binding of recombinant foamy virus envelope protein to host cells correlates with susceptibility to infection. Virology. 1999, 255: 228-236.View ArticlePubMedGoogle Scholar
- Yap MW, Lindemann D, Stanke N, Reh J, Westphal D, Hanenberg H, Ohkura S, Stoye JP: Restriction of foamy viruses by primate Trim5alpha. J Virol. 2008, 82: 5429-5439.PubMed CentralView ArticlePubMedGoogle Scholar
- Perkovic M, Schmidt S, Marino D, Russell RA, Stauch B, Hofmann H, Kopietz F, Kloke B-P, Zielonka J, Ströver H, Hermle J, Lindemann D, Pathak VK, Schneider G, Löchelt M, Cichutek K, Münk C: Species-specific inhibition of APOBEC3C by the prototype foamy virus protein bet. J Biol Chem. 2009, 284: 5819-5826.PubMed CentralView ArticlePubMedGoogle Scholar
- Opazo JC, Wildman DE, Prychitko T, Johnson RM, Goodman M: Phylogenetic relationships and divergence times among New World monkeys (Platyrrhini, Primates). Mol Phylogenet Evol. 2006, 40: 274-280.View ArticlePubMedGoogle Scholar
- Schneider H: The current status of the New World monkey phylogeny. An Acad Bras Cienc. 2000, 72: 165-172.View ArticlePubMedGoogle Scholar
- Muniz CP, Troncoso LL, Moreira M, Soares E, Pissinatti A, Bonvicino CR, Seuánez HN, Sharma B, Jia H, Shankar A, Switzer WM, Santos AF, Soares M: Identification and characterization of highly divergent simian foamy viruses in a wide range of new world primates from Brazil. PLoS One. 2013, 8: e67568-PubMed CentralView ArticlePubMedGoogle Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S: MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013, 30: 2725-2729.PubMed CentralView ArticlePubMedGoogle Scholar
- Hawkesworth C, Kelley S, Turner S, Le Roex A, Storey B: Mantle processes during Gondwana break-up and dispersal. J Afr Earth Sci. 1999, 28: 239-261.View ArticleGoogle Scholar
- de Wit MJ: MADAGASCAR: heads it’s a continent, tails it’s an island. Annu Rev Earth Planet Sci. 2003, 31: 213-248.View ArticleGoogle Scholar
- Hussain AI, Shanmugam V, Bhullar VB, Beer BE, Vallet D, Gautier-Hion A, Wolfe ND, Karesh WB, Kilbourn AM, Tooze Z, Heneine W, Switzer WM: Screening for simian foamy virus infection by using a combined antigen Western blot assay: evidence for a wide distribution among Old World primates and identification of four new divergent viruses. Virology. 2003, 309: 248-257.View ArticlePubMedGoogle Scholar
- Wolfe ND, Switzer WM, Carr JK, Bhullar VB, Shanmugam V, Tamoufe U, Prosser AT, Torimiro JN, Wright A, Mpoudi-Ngole E, McCutchan FE, Birx DL, Folks TM, Burke DS, Heneine W: Naturally acquired simian retrovirus infections in central African hunters. Lancet. 2004, 363: 932-937.View ArticlePubMedGoogle Scholar
- Kumar S, Subramanian S: Mutation rates in mammalian genomes. Proc Natl Acad Sci U S A. 2002, 99: 803-808.PubMed CentralView ArticlePubMedGoogle Scholar
- Stamatakis A, Hoover P, Rougemont J: A rapid bootstrap algorithm for the RAxML Web servers. Syst Biol. 2008, 57: 758-771.View ArticlePubMedGoogle Scholar
- Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21: 2104-2105.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574.View ArticlePubMedGoogle Scholar
- MATLAB: MATLAB and Statistics Toolbox Release R2012a. 2010, Massachusetts: The MathWorks IncGoogle Scholar
- Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24: 1586-1591.View ArticlePubMedGoogle Scholar
- Shimodaira H, Hasegawa M: CONSEL: for assessing the confidence of phylogenetic tree selection. Bioinformatics. 2001, 17: 1246-1247.View 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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.