Abstract
Human immunodeficiency virus, a primate lentivirus (PLV), causes AIDS in humans, whereas most PLVs are less or not pathogenic in monkeys. These notions suggest that the co-evolutionary process of PLVs and their hosts associates with viral pathogenicity and therefore, that elucidating the history of virus-host co-evolution is one of the most intriguing topics in the field of virology. To address this, recent studies have focused on the interplay between intrinsic anti-viral proteins, such as tetherin and viral antagonists. Through an experimental-phylogenetic approach, here we investigate the co-evolutionary interplay between tribe Cercopithecini tetherin and viral antagonists, Nef and Vpu. We reveal that tribe Cercopithecini tetherins are positively selected, possibly triggered by ancient Nef-like factor(s). We reconstruct the ancestral sequence of tribe Cercopithecini tetherin and demonstrate that all Nef proteins are capable of antagonizing ancestral Cercopithecini tetherin. Further, we consider the significance of evolutionary arms race between tribe Cercopithecini and their PLVs.
Similar content being viewed by others
Introduction
Based on the sequence similarity, the following two issues have been widely accepted: (i) human immunodeficiency virus type 1 (HIV-1), the causative agent of acquired immunodeficiency syndrome, emerged from zoonotic transmission of a simian immunodeficiency virus (SIV) in chimpanzee (SIVcpz) to humans around 100 years ago1,2,3; and (ii) SIVcpz appears to have emerged from the recombination of two lineages of SIVs from Old World monkeys (OWMs): SIVgsn/mon/mus lineage from greater-spot nosed monkey (Cercopithecus nictitans; GSN), mona monkey (Cercopithecus mona; MON) and mustached monkey (Cercopithecus cephus; MUS) and SIVrcm from red-capped mangabey (Cercocebus torquatus)4. Understanding the evolutionary history of primate lentiviruses (PLVs) including HIVs and SIVs is one of the most important and interesting topics in the field of retrovirology. However, because of their multiple cross-species transmissions and complicated recombination, it is difficult to elucidate how genetic conflicts between the ancient SIVs and their respective host species resulted in evolution and diversification.
OWMs, the family Cercopithecidae, are composed of 12 genera and a subfamily5. All SIVs identified so far encode 8 common genes: gag, pol, env, tat, rev, vpr, vif and nef 6. Among more than 40 SIVs, which have been identified in OWMs residing in Africa7, only 4 kinds of SIVs, SIVgsn in GSN, SIVmon in MON, SIVmus in MUS and SIVden in Dent’s mona monkey (Cercopithecus denti) encode an additional accessory gene, vpu8. Importantly, the SIVs encoding vpu have been identified only in the monkeys belonging to tribe Cercopithecini including the genus Cercopithecus, strongly suggesting that the vpu gene has emerged in the evolution and transmission of SIVs in this tribe5,7,8.
To elucidate the co-evolutionary relationship between SIVs and their hosts, recent investigations have experimentally addressed the evolutionary conflict between viral and host proteins8,9,10 that stems from the “Red Queen hypothesis”11 or “evolutionary arms race” concept. Such an approach can be the way to explain the co-evolutionary history of SIVs and their host species. For example, Vif, a common protein encoded by all PLVs, has a robust ability to counteract a cellular anti-PLV restriction factor, apolipoprotein B mRNA editing enzyme catalytic polypeptide-like 3G12. In addition, another anti-PLV restriction factor, SAM domain and HD domain 1 (SAMHD1), can be antagonized by the viral accessory proteins, Vpr or Vpx13,14. The vpx gene is encoded in certain SIV lineages and HIV type 2 and it has been assumed that the vpx gene evolved from gene duplication of its ancestral gene, vpr15. Moreover, based on an experimental-phylogenetic investigation, Lim et al. recently proposed that the evolutionary interaction between Vpr/Vpx and SAMHD1 has undergone the following four steps: (i) Vpr acquired anti-SAMHD1 activity; (ii) ancestral SIV(s) created vpx by the gene duplication of vpr; (iii) Vpr transferred its anti-SAMHD1 activity to Vpx16. Namely, anti-SAMHD1 ability has been transferred from an old gene (vpr) to a new gene (vpx) during the co-evolution of SIVs and their hosts.
One of the most complicated examples of the co-evolutionary relationship between OWMs and their SIVs is tetherin (also known as bone marrow stromal antigen 2, CD317 and HM1.24) and its viral antagonists8,17. Tetherin inhibits the release of nascent viral particles from virus-producing cells18,19. Similar to the case of SAMHD1 and Vpr/Vpx, SIVs encode two kinds of anti-tetherin antagonists: Nef and Vpu8. Nef is encoded in all SIVs and most SIVs antagonize tetherins of their hosts by Nef8. On the other hand, as described above, Vpu is encoded in certain SIVs in OWMs and the Vpu proteins of these SIVs potently antagonize tetherins of their hosts8,20. However, when, why and how the vpu gene was acquired in certain SIV lineages during their evolution is still unclear.
In this study, we particularly focus on the OWMs belonging to the tribe Cercopithecini and their SIVs and perform investigations based on molecular phylogenetics and evolution, experimental virology and structural biology. We reveal that the tetherins of the tribe Cercopithecini are under strong positive selection. In addition, we construct the ancestral sequences of tribe Cercopithecini tetherin and experimentally demonstrate that all Nef proteins of the SIVs isolated from the tribe Cercopithecini retain antagonistic ability to the ancestral tetherin of tribe Cercopithecini. Moreover, we estimate the time of vpu acquisition in certain SIV lineages and further, discuss the reason why vpu has been created and/or acquired from various scientific fields of view.
Results
Evolution of primate tetherin and CD4
Since lentiviral Nef and Vpu proteins have the common ability to down-regulate tetherin as well as CD417,19,20,21, we set out to perform molecular phylogenetic analyses on primate tetherin and CD4. In this study, we newly identified 11 tetherin sequences of 8 different OWMs belonging to the tribe Cercopithecini (1 Campbell’s mona monkey, 1 mustached monkey, 2 Sclater’s monkeys, 1 L’Hoest’s monkey, 2 Sykes’ monkeys, 2 red-eared monkeys, 1 red-tailed monkey and 1 sun-tailed monkey; listed in Table 1) and 3 CD4 sequences of 2 different OWMs (1 MUS and 2 sooty mangabeys; listed in Table 2). As shown in Fig. 1a,b, each family or infraorder (i.e., Hominoids, OWMs, or NWMs) respectively formed a monophyletic cluster on the reconstructed trees of both tetherin and CD4. On the other hand, within the cluster of Cercopithecini, the tetherins of certain Cercopithecini monkeys, particularly mustached monkey, red-eared monkey and Sclater’s monkey, did not form a monophyletic subcluster, respectively (Fig. 1a). This indicates that the nucleotide sequence of certain Cercopithecini tetherins, particularly mustached monkey, Red-eared monkey and Sclater’s monkey, are highly similar.
To detect positive selection in the evolution of primate tetherin and CD4, we estimated the nonsynonymous to synonymous (dN/dS) ratios. The two pairs of site models in PAML produced similar results and the results obtained from M7 (neutral model) versus M8 (selection model) comparisons are shown in Fig. 1c,d. Consistent with previous reports22,23,24 including ours25, the dN/dS ratio of primate tetherin was significantly greater than one for full-length (43.04), cytoplasmic tail (CT; 26.76) and transmembrane domain (TMD; 11.08) (Fig. 1c). Also, three codons, 9 (dN/dS = 5.2), 14 (dN/dS = 5.1) and 17 (dN/dS = 5.2) in primate tetherin, were identified to be positively selected sites with posterior probability greater than 0.95 (Fig. 1e). These findings on primate tetherin indicate that the functionally important regions of primate tetherin, particularly CT and TMD, have evolved under strong positive selection, which is in agreement with previous reports22,23,24,25.
On the other hand, the dN/dS ratio of primate CD4 was significantly greater than one for full-length (28.91) and the extracellular domain (ECD; 31.64) and six codons in the ECD, 48 (dN/dS = 3.6), 73 (dN/dS = 3.6), 77 (dN/dS = 3.6), 80 (dN/dS = 3.7), 113 (dN/dS = 3.7) and 265 (dN/dS = 3.7), were identified as positively selected codons by the site model analysis (Fig. 1f). These findings on primate CD4 suggest that the ECD of primate CD4 has evolved under strong positive selection. Because the ECD of CD4 molecule is homologous to immunoglobulins and plays a crucial role for immune recognition and immune responses26, the evolution and diversification of the immune system may closely associate with the positive selection observed in this study (Fig. 1d,f). In contrast to primate tetherin, positive selection was detected in neither TMD nor CT of primate CD4 (Fig. 1d,f). To down-regulate these cellular proteins, Nef targets the CTs of tetherin and CD421,27, whereas Vpu targets the TMD of tetherin and the CT of CD4, respectively8,28. Therefore, our findings suggest that primate tetherin but not CD4 has experienced the positive selection elicited by Nef and/or Vpu during evolution.
Positive selection detected in the evolution of tribe Cercopithecini tetherin
Among the SIVs in OWMs (identical to the family Cercopithecidae), vpu-positive SIVs have been identified only in the monkeys belonging to the tribe Cercopithecini (represented in pink in Fig. 1a,b). To elucidate the evolutionary interplay between Nef/Vpu and tetherin, we particularly focused on the tetherins of this tribe. The phylogenetic tree of 22 tetherins belonging to the tribe Cercopithecini showed that the tetherins of the hosts of vpu-positive SIVs intermingle with those of vpu-negative SIVs (Fig. 2a), suggesting that the presence of vpu did not result in the convergent evolution of Cercopithecini tetherin. Also, the site model analysis revealed that the dN/dS ratio of Cercopithecini tetherin was significantly greater than one for full-length (47.56) and CT (21.13) as well as ECD (12.15) (Fig. 2b). Four codons in the CT, 14 (dN/dS = 8.6), 16 (dN/dS = 8.5), 17 (dN/dS = 8.6) and 24 (dN/dS = 8.4) and two codons in the ECD, 67 (dN/dS = 8.3) and 99 (dN/dS = 8.4), were identified to be positively selected (Fig. 2c).
We then classified 22 Cercopithecini tetherins into two groups: the hosts of vpu-positive SIVs and those of vpu-negative SIVs. Because SIV has not been identified in Sclater’s monkeys29, we excluded the tetherin sequences of 2 Sclater’s monkeys from this classification. As shown in Fig. 3a, the branch-site tests in PAML revealed that the likelihood ratio test was significant with P < 0.01 in the analysis of the tetherins of the hosts of vpu-negative SIVs as well as 22 Cercopithecini tetherins, suggesting that positive selection has most likely operated on the tetherins of the monkeys infected with vpu-negative SIVs. Also, the site model revealed that the dN/dS ratio of the tetherins of the hosts of vpu-negative SIVs was significantly greater than one for full-length (42.55), CT (25.92) and ECD (12.15) (Fig. 3b, left) and nine codons positioned at 14, 16, 17, 24, 34, 67, 99, 100 and 159 were identified to be positively selected (Fig. 3c. left). In addition, thirteen codons including the nine codons detected by the site model (indicated by asterisks in Fig. 3d) were identified as positively selected sites by the random effects likelihood (REL) analysis implemented in the HyPhy package with Bayes factor greater than 50 (Fig. 3d, left). Furthermore, we constructed the ancestral sequence of the 22 tetherins of tribe Cercopithecini, using the codeml program in PAML (indicated by a red star in Fig. 2a. The sequence information is available in Supplemental dataset) and constructed the structure homology model of the ECD of ancestral Cercopithecini tetherin (Fig. 3e). By mapping the five positively selected sites in the ECD (positioned at 63, 67, 99, 100 and 159), we found that these amino acids were located on the same aspect of the alpha helix structure (Fig. 3e).
In contrast to the tetherins from the hosts of vpu-negative SIVs, it was notable that positive selection was not detected in the tetherins from the monkeys infected with vpu-positive SIVs by the branch-site model (Fig. 3a) and REL analyses (Fig. 3d, right). Although the site model showed that the dN/dS ratio of the tetherins of the hosts of vpu-positive SIVs was significantly greater than one for full-length (2Δl = 10.42, P < 0.01) (Fig. 3b, right), no positive selection was detected at the significant level (P < 0.05) for respective domains (Fig. 3b, right) and codons (Fig. 3c, right), which basically agreed with the result obtained from the branch-site test (Fig. 3a). To ask whether the difference in the positive selection between these two groups could be attributed to the difference in the number of tetherin sequences included, we performed genetic diversity analyses using MEGA630. As shown in Fig. 3f, the genetic diversity of tetherins of these two groups was comparable. Taken together, these results indicate that the much weaker selective pressure detected for the tetherins of the hosts of vpu-positive SIVs is not likely due to the smaller sample size of this group, but suggesting that Vpu did not exert a strong selective pressure on the tetherins of the hosts of vpu-positive SIVs.
Antagonism of Cercopithecini tetherin by SIV Nef
To directly evaluate the anti-viral activity of Cercopithecini tetherin and the antagonistic ability of SIV Nef proteins, we prepared an expression plasmid for the constructed ancestral Cercopithecini tetherin. Western blotting (Fig. 4a) and TZM-bl assay (Fig. 4b) revealed that increasing amounts of the ancestral Cercopithecini tetherin resulted in a dose-dependent decrease in the release of nascent virions. We confirmed that the infectious virus in the culture supernatant correlated strongly with the amount of supernatant viral p24 antigen, a physical measure of virion content (r = 0.947, P = 0.000031; Supplementary Fig. 2), which is consistent with previous reports20,22,31 and validates the tetherin impact on virus release.
We then prepared expression plasmids for 14 strains of SIV Nefs, which have been identified in the 12 species of Cercopithecini monkeys so far and investigated whether these Nef proteins have the ability to antagonize the ancestral Cercopithecini tetherin. As shown in Fig. 4a, none of the Nef proteins affected expression levels of tetherin and Gag, particularly Gag precursor (Pr55Gag; Fig. 4a), which is consistent with previous reports22,25,27. Moreover, we revealed that all Nef proteins enhanced viral release in the presence of the ancestral Cercopithecini tetherin (Fig. 4b). Importantly, the 6 Nef proteins of vpu-positive SIVs significantly augmented viral release (Fig. 4b). These findings directly demonstrate that all SIV Nef proteins isolated from the tribe Cercopithecini potently antagonize Cercopithecini tetherin regardless of whether or not they had acquired vpu.
Genetic and geographical consideration of the acquisition of vpu gene
To further assess the possibility that Vpu has exerted selective pressure on the tetherins of tribe Cercopithecini, we performed a Bayesian evolutionary analysis for dating the time of vpu gene acquisition. As shown in Fig. 5a, our analysis revealed that the vpu gene has been independently acquired twice, which were estimated to have occurred in SIVgsn/mon/mus lineage (nodes 8) around 19,418 years ago and in SIVden lineage (node 13) around 19,218 years ago, respectively (Table 3). In addition, it has been reported that the gain-of-function of Vpr to degrade SAMHD1 occurred after the divergence from SIVsun/lhoest lineages16. Our analyses revealed that the Vpr neofunction occurred after 35,766 years ago (Fig. 5a and Table 3). These findings suggest that vpu was acquired by the two lineages of SIVs around 20,000 years ago, which were relatively recent events in the evolutionary history of SIVs compared to the neofunction of Vpr.
We then considered the acquisition of vpu gene in geographic terms. As considered in the previous papers29,32,33, the habitats of the three species of OWMs, GSN, MON and MUS overlapped in West Africa including Cameroon, Gabon and Nigeria and Republic of the Congo, while the habitat of DEN was geographically separated (Fig. 5b). Therefore, it is plausible that vpu gene has been independently acquired by certain SIV(s) infecting the monkeys in these two separate geographic areas.
Furthermore, we considered the event of vpu acquisition in terms of viral genetics. Because of the restricted genome size of RNA viruses including lentiviruses, ‘genome compression’, which is caused by the use of overlapping genes, is a known characteristic of RNA viruses34. In fact, the 3′ end of all vpu genes of PLVs overlaps with the 5′ end of env (Fig. 5c, top). Although a previous study has suggested that there is no preference in the direction of frameshift in newly acquired viral genes34, we found that all vpu genes of OWM SIVs involved + 1 (forward) frameshift compared to env (data not shown). This suggests that the manner of vpu acquisition may be common in the two lineages of vpu-positive SIVs (Fig. 5a). Moreover, we measured the nucleotide length between the 3′ ends of 4 viral genes (tat1, rev1, vpr and vif) and the 5′ end of env, where vpu is encoded in certain SIVs. As expected, the nucleotide lengths of these 4 regions in vpu-positive SIVs were clearly longer than those in vpu-negative SIVs (Fig. 5c). However, it was of interest that the nucleotide lengths from the 3′ ends of tat1, rev1 and vpr, but not of vif, in the 3 strains of SIVdeb were significantly longer than those of the other vpu-negative SIVs (Fig. 5c). Although the mechanism of new gene acquisition by viruses is still unknown34,35, these findings suggest that the nucleotide length in this region can vary in vpu-negative SIVs and that this genomic region may be adequate for SIVs to create and/or acquire new gene(s).
Discussion
In this study, we newly determined 11 tetherin sequences of the 8 species of the tribe Cercopithecini as well as 3 CD4 sequences of OWMs. In addition, we performed the in-depth molecular phylogenetic analyses and revealed that the tetherins of the tribe Cercopithecini, particularly those of the host monkeys of vpu-negative SIVs, are under strong positive selection. Furthermore, we constructed the ancestral sequence of tribe Cercopithecini tetherin and demonstrated that the ancestral Cercopithecini tetherin has a robust ability to inhibit viral release. In this regard, since a previous paper has shown that the “artificial” tetherin, which artificially forms the same topology to tetherin, sufficiently confers anti-viral activity36, it might not be so surprising that the ancestral Cercopithecini tetherin estimated in this study exhibited anti-viral ability. Nevertheless, here we demonstrated that the anti-viral activity of the ancestral Cercopithecini tetherin is strongly antagonized by all SIV Nef proteins we used. Particularly noteworthy was that Nef’s antagonistic activity against Cercopithecini tetherin is not associated with the presence of vpu gene. Moreover, we estimated the time of vpu gene acquisition in certain SIVs and further considered its significance.
Among the order Primates, the domains of tetherin and CD4 under positive selection differed: primate tetherin has been under positive selection in the CT and TMD, whereas the ECD of primate CD4 was positively selected (Fig. 1). These findings suggest that the factors triggering selective pressure on these molecules differ from each other. In the case of CD4, the ECD plays critical roles in immune control such as the recognition of major histocompatibility complex class II26. Moreover, CD4 is utilized for the invasion of PLV through interacting its ECD with their envelope glycoprotein (Env)26. Therefore, it is plausible that the ECD of CD4 has been positively selected through the immune pressures outside of the cells and the interaction with PLV Env26,37. Moreover, positive selection was detected in neither TMD nor CT of primate CD4 (Fig. 1d,f), suggesting that these domains are evolutionary stable. It is known that both Nef and Vpu target CD4 CT for the down-regulation8,21,27. Therefore, these findings suggest that primate CD4 has not been under positive selection caused by Nef, Vpu and their ancestors. To the best of our knowledge, this is the first study evaluating the selective pressure on primate CD4 gene in-depth.
It was reported that the tribe Cercopithecini diversified approximately 8.2 million years ago (indicated in Fig. 1a)38. On the other hand, here we revealed that the vpu gene has been acquired in certain SIV lineages around 20,000 years ago (Fig. 5a and Table 3). These findings strongly suggest that the Vpu proteins encoded by the present-day SIVs cannot be the source of the selective pressure on Cercopithecini tetherin. In addition, the convergent evolution was not observed on the tetherins of vpu-positive SIV hosts (Fig. 2a). Furthermore, it was surprising that significant positive selection was not detected in any domains and codons of the tetherins of vpu-positive SIV hosts (Fig. 3a–d). Since there was no positive selection detected on the TMD of Cercopithecini tetherin (Fig. 2b), these findings suggest that ancestral Vpu-like factor(s) were not the selective pressure on the tetherins of tribe Cercopithecini.
In contrast to the tetherins of vpu-positive SIV hosts, the tetherins of the monkeys infected with vpu-negative SIVs exhibited strong positive selection in the CT and ECD (Fig. 3b–d). Interestingly, four out of the five positively selected sites in the ECD are located on the same aspect of its alpha-helix structure (Fig. 3e). This observation is reminiscent of the four amino acids, I34, L37, L41 and T45, in the TMD of human tetherin28. We have previously reported that these four amino acids in the TMD of human tetherin are responsible for HIV-1 Vpu-mediated antagonism and are located on the same helical face of the TMD28. Further study has revealed that the four amino acids in the TMD of human tetherin form the helix-helix intermolecular interaction with HIV-1 Vpu39. These findings raise the possibility that the four amino acids located on the same face of ECD may be positively selected from pressure caused by the other viral antagonists of tetherin. In fact, it has been revealed that the glycoprotein of Ebola virus, which sporadically causes epidemics in humans and primates residing in the central African countries such as the Democratic Republic of the Congo (formerly Zaire) and Sudan40, interacts with and antagonizes tetherin41. Moreover, Env of certain lentiviruses potently antagonize tetherin42,43,44,45. Therefore, these viral antagonists of tetherins and/or their ancestors could have exerted the selective pressure on the ECD of Cercopithecini tetherin. Moreover, although the main habitat of vpu-positive SIV hosts is West Africa (Fig. 5b), the monkeys infected with vpu-negative SIVs reside in a broad area of Africa including the central African countries46. This further suggests that the tetherins of vpu-negative SIV hosts have had the opportunity to be exposed to different pressures compared to those of vpu-positive SIVs.
As shown in Fig. 4b, we demonstrated that the infectivity of culture supernatant was significantly suppressed by the ancestral Cercopithecini tetherin. These findings suggest that the ancestral Cercopithecini tetherin possesses a robust activity to impair viral release and that the anti-viral ability of tetherin has been maintained in its evolution. On the other hand, it was surprising that all SIV Nef proteins including those of vpu-positive SIVs (e.g., SIVgsn, SIVmon, SIVmus and SIVden) were capable of antagonizing Cercopithecini tetherin-mediated anti-viral ability (Fig. 4). It has been demonstrated that the Vpu proteins of vpu-positive SIVs antagonize tetherins of their natural hosts20. Therefore, these observations suggest that the vpu gene was not necessarily been acquired by certain SIVs to gain a novel anti-tetherin antagonist. In the case of the evolutionary interplay between Vpr/Vpx and SAMHD1, three evolutionary steps have been proposed: (i) acquisition of anti-SAMHD1 activity by ancestral Vpr (i.e., Vpr neofunction); (ii) creation of vpx by gene duplication in certain SIV lineages; and (iii) transfer of anti-SAMHD1 activity from Vpr to Vpx16. In contrast to the scenario of Vpr/Vpx and SAMHD1, our findings suggest that Nef proteins of vpu-positive SIVs have not lost their anti-tetherin activity even though a new tetherin antagonist, Vpu, was acquired and that the transfer of anti-tetherin activity from Nef to Vpu has not occurred. Moreover, here we estimated that Vpr neofunction occurred around 36,000 years ago, which is relatively older than the acquisition of the vpu gene (Fig. 5a). Therefore, these findings imply that Nef still maintains anti-tetherin activity in vpu-positive SIVs because vpu is relatively younger than vpr/vpx and that Nefs of vpu-positive SIVs may transfer their anti-tetherin activity to Vpu in the future.
As shown in Fig. 5a, our results suggest that the vpu gene was independently acquired in two SIV lineages: SIVgsn/mon/mus and SIVden. This raises three possibilities. First, it might be possible that the vpu gene acquired in SIVgsn/mon/mus lineage has been horizontally transferred to SIVden lineage and vice versa. Because the two virus lineages, SIVgsn/mon/mus and SIVden, share strikingly similar genomic features (e.g., very similar vpu genes with nearly identical locations), the most parsimonious explanation is that the acquisition of vpu gene was not independent events, but rather descended from a single original event. However, SIVgsn/mon/mus is phylogenetically divergent from SIVden (Fig. 5a) and the habitats of the monkeys infected with these two SIV lineages are geographically separated (Fig. 5b)46. These two notions argue against this first hypothesis. Second, there is a possibility that the vpu gene was acquired in the common ancestor of these two SIV lineages (i.e., node 14 of Fig. 5a) and then certain SIVs (e.g., SIVtal, SIVsyk and SIVdeb) lost their own vpu. It might be possible for some viruses to lose the vpu gene because there are other routes to overcome OWM tetherin (e.g., Nef, Env)8. If something similar happened in an ancestral virus, it would relieve selective pressure to maintain a vpu gene. The third possibility is that the acquisition of vpu independently occurred twice in the two SIV lineages, perhaps the most feasible possibility. Further, it should be notified that the molecular clock analyses shown herein rely on the accessible information to date. It means that the information obtained in the future may affect the estimated age of vpu acquisition. Nevertheless, this is the first study inferring the time of accessory gene acquisition/generation by PLVs.
So, how was the vpu gene acquired and/or created ? Because the genome size of RNA viruses is strictly restricted, it has been assumed that RNA viruses evolutionary repeat trial-and-error to obtain new favorable genes34,35. Interestingly, we found that the nucleotide lengths between the 3′ end of tat1 and the 5′ end of env of prosimian endogenous lentiviruses (PSIVs), which share a common ancestor with modern SIVs47,48, were much shorter than those of SIVs (Fig. 5c). In this regard, PSIVs encode dUTPase in pol region, while PLVs including SIVs do not6. This raises a possibility that the loss of dUTPase in SIV ancestors has relaxed the restriction of viral genome space, which allowed viruses to acquire and/or create novel genes. In fact, PLVs encode multiple viral genes around the region encoding vpu more than the other lentiviruses, which encode dUTPase6,10. Therefore, it is plausible that PLVs have gained the chance to acquire and/or create new genes by losing dUTPase and that the acquisition of vpu might be one of the consequences.
Methods
Ethic statement
To determine the sequence of tetherin, blood was collected from wild-caught monkeys: 1 Campbell’s mona monkey (Cercopithecus campbelli), 1 mustached monkey (Cercopithecus cephus), 2 Sclater’s monkeys (Cercopithecus sclateri), 1 L’Hoest’s monkey (Cercopithecus lhoesti), 2 red-eared monkeys (Cercopithecus erythrotis), 1 red-tailed monkey (Cercopithecus ascanis) and 1 sun-tailed monkey (Cercopithecus solatus) according to the Guide for the Care and Use of Laboratory Animals49 under a NIAID Animal Care and Use Committee-approved protocol50,51. These procedures were approved by NIAID.
Sequencing PCR
Genomic DNA was extracted from cryopreserved peripheral blood mononuclear cells (PBMCs) of these 9 monkeys50,51 by using DNeasy kit (Qiagen). Also, genomic DNA was extracted from the body hair root of 2 Sykes’ monkeys (Cercopithecus albogularis), which are kept in the Japan Monkey Centre, Inuyama, Aichi, Japan) by using DNA Extractor FM kit (Wako). PCR was performed by using PfuUltra High Fidelity DNA polymerase (Agilent Technologies) and the following primers: 5′-CAG CTA GAG GGG AGA TCT GGA TG-3′; 5′-CTC ACT GAC CAG CTT CCT GGG-3′, which were used in our previous study25. The obtained PCR products were purified by gel extraction and directly sequenced by using BigDye Terminator v3.1 cycle sequencing kit (Applied Biosystems) with the two primers described above and the following 4 primers: 5′-GGA CTT CAC CAG ACC CTG AA-3′; 5′-TTC AGG GTC TGG TGA AGT CC-3′; 5′-TCT CTC CTT TGC TCC CAA AA-3′; 5′-TTT TGG GAG CAA AGG AGA GA-3′. To determine the sequence of CD4, RNA was extracted from the cryopreserved PBMCs by using RNeasy Mini Kit (Qiagen). Reverse transcription was performed by using Thermoscript RT-PCR System (Life Technologies) and RT-PCR was performed by using Platinum Taq DNA polymerase High Fidelity (Life Technologies) and the following primers: 5′-CAG CAA GGC CAC AAT GAA C-3′ and 5′-TGC CTC AAA TGG GGC TAC-3′. The obtained RT-PCR product was purified by gel extraction and then cloned by using TOPO TA Cloning Kit (Life Technologies). The sequencing PCR was performed by using ABI Prism 3130 xl genetic analyzer (Applied Biosystems) and the data was analyzed by Sequencher v5.1 software (Gene Codes Corporation).
Molecular phylogenetic analyses
The molecular phylogenetic analyses were performed as previously described25. Briefly, the 11 tetherin sequences newly identified in this study were aligned with 47 primate tetherin sequences (listed in Table 1) by using ClustalW implemented in MEGA630. Also, the 3 CD4 sequences newly identified in this study were aligned with 19 primate CD4 sequences (listed in Table 2) as described above. The alignments were verified manually at amino acid level. Then the phylogenetic trees were reconstructed using neighbor-joining (NJ) method52 with MEGA630 (Fig. 1a,b) and maximum-likelihood (ML) method with PhyML53 (Fig. S1A). Note that the phylogenetic trees of 58 primate tetherins reconstructed by these two methods yielded similar topology with partial difference in the relationships between certain tetherins; particularly mustached monkey, Red-eared monkey and Sclater’s monkey. We assume that these minor differences are due to higher similarity among these species, because each species did not form a monophyletic cluster on the tree and bootstrap support for most splits between the sequences were fairly low (< 75%, data not shown). Moreover, both NJ and ML trees were used for further PAML analyses and the small topological difference between these two trees did not affect the results (Fig. 1c and S1B). Furthermore, we confirmed that the nucleotide sequence of the ancestral Cercopithecini tetherin inferred by NJ tree was identical to that by ML tree (data not shown). We then conducted the analysis to detect positive selection along the tree. To infer positive selection across various primate lineages, two pairs of site models implemented in the PAML package v 4.754 were used to conduct the likelihood ratio tests for 58 tetherin genes (Fig. 1c,e) and 22 CD4 genes (Fig. 1d,f), respectively: M1 (neutral model) versus M2 (selection model) and M7 (neutral model) versus M8 (selection model). The REL method in HyPhy55 was also employed to detect positive selection (Fig. 3d). The ancestral Cercopithecini tetherin was inferred by using site model in the PAML analysis (Fig. 2a). Since we were particularly interested in whether the clades of tribe Cercopithecini, the hosts of vpu-negative SIVs and those of vpu-positive SIVs have evolved under positive selection, we further focused on the 22 tetherins of this clade. First, the branch-site model in PAML was employed for the analysis. This model allows dN/dS ratio to vary both among sites and branches, which is very useful for detecting positive selection along a particular lineage or clade (pre-specified as foreground branches)54. In our analysis, all the 22 tetherins of tribe Cercopithecini, 14 tetherins of the hosts of vpu-negative SIVs and 6 tetherins of the hosts of vpu-positive SIVs, were respectively specified as the foreground branches (Fig. 3a). Next, the site model in PAML and REL method in HyPhy were performed to these sequences (Fig. 3b–d). Moreover, the genetic distance (Fig. 3f) was calculated for the 14 tetherin genes of the hosts of vpu-negative SIV and the 6 genes of the hosts of vpu-positive SIV, respectively, by using MEGA630. We computed the overall mean distance was computed by using Tamura-Nei model56 with 100 bootstrap replications.
Protein homology modeling
The 3D structure of the ECD of ancestral Cercopithecini tetherin (Fig. 3e) was simulated by the Swiss-Model server (http://swissmodel.expasy.org/) using the crystal structure of the ECD of human tetherin (PDB code: 3MQB) as the template31.
BEAST analysis
The full-genome sequences of 34 SIV strains (listed in Table 4) were retrieved from the HIV Sequences Database (http://www.hiv.lanl.gov/content/sequence). Then, the gag, pol, env and vif genes were extracted from each viral sequence and were respectively aligned using MAFFT57. The resulting alignments were manually verified at the amino acid level. We then performed the single breakpoint analysis55 implemented in the HyPhy package to test for the recombination in each aligned dataset. The results showed that no evidence of recombination was detected in all four analyses. We further performed Gblocks (http://molevol.cmima.csic.es/castresana/Gblocks_server.html) v 9.1b to remove poorly aligned regions from these alignments. Finally, a concatenated dataset of all four genes was created by using BioEdit (http://www.mbio.ncsu.edu/bioedit/bioedit.html) v 7.2.5. To infer the time of most recent common ancestors of these SIVs, a Bayesian approach implemented in the BEAST package (http://beast.bio.ed.ac.uk) v 1.7.558 was employed. This analysis was conducted at the amino acid level, as the 34 SIV strains were highly divergent. We used a JTT substitution model with gamma-distributed rate variation among sites. The uncorrelated lognormal relaxed molecular clock model was employed to estimate substitution rates and the Yule process of speciation was used for the tree prior. We specified a uniform distributed prior (30,000–130,000 yr, initial = 70,000 yr) for the age of the root of the tree59. Five independent Markov Chain Monte Carlo (MCMC) analyses were run for 10–20 million generations with sampling every 1000 generations. We then used the program Tracer (tree.bio.ed.ac.uk/software/tracer) v 1.6 to check for the convergence and to confirm that the effective sample size (ESS) value was higher than 200 for all runs. The maximum clade credibility (MCC) tree was generated by summarizing the sample of trees produced by BEAST after a 10% burn-in using the TreeAnnotator program v 1.7.5 and the resulting MCC tree was viewed using FigTree (http://tree.bio.ed.ac.uk/) v 1.4.2.
Plasmid construction
The HA-tagged Nef expression plasmids of SIVgsn (strain 99CM166), SIVmon (strains 99CMCML1 and NG1), (SIVmus strains 01CM1085 and 01CM1239) were used in our previous study25. The Nef open reading frames (ORFs) of SIVden (strain CD1), SIVagmSab (strain SAB1), SIVsyk (strain KE44), SIVlhoest (strain LHO7), SIVsun (strain L14), SIVtal (strain 00CM266) and SIVdeb (strains CM5 and Kin1) were obtained from GeneArt Gene Synthesis service (Life Technologies). The Nef ORF of SIVagmVer strain TYO1 was obtained by PCR using pSA212 (an infectious molecular clone of SIVagmVer strain TYO1)60 as the template and the following primers: TYO1 Nef-Fwd, 5′-TTT TTT TTT CTA GAA TGG GCT CGC AGA ACT CA-3′; TYO1 Nef-Rev, 5′-TAT ATA TAT ATA GAT ATC CTT CCT CTT CAC CAG CC-3′. The resultant DNA fragment was digested with XbaI and EcoRV and was inserted into the XbaI-EcoRV site of pCGCG vector.
Cell culture and transfection
HEK293T cells and TZM-bl cells (obtained through NIH AIDS Research and Reference Reagent Program) were maintained in Dulbecco’s modified Eagle medium (Sigma) containing 10% heat-inactivated FCS and antibiotics. Transfection was performed by using PEI Max (GE Healthcare) according to the manufacturer’s protocol. Various amounts of KGC-tagged ancestral Cercopithecini tetherin expression plasmids (0, 10, 40 ng) and pNL4-3ΔvpuΔnef (1,200 ng) were cointroduced with or without respective SIV Nef expression plasmid (400 ng) into HEK293T cells. At 48 hours post-transfection, the culture supernatants and transfected cells were harvested and were respectively used for TZM-bl assay and Western blotting as described below.
Western blotting and TZM-bl assay
The culture supernatant harvested at 48 hours post-transfection was centrifuged to remove cells and produce virus suspensions. The infectivity of virus suspensions was measured by TZM-bl assay as previously described25. Briefly, 100 μl of the virus solution was inoculated into TZM-bl cells in 96-well plate (Nunc) and the β-galactosidase activity was measured by using the Galacto-Star mammalian reporter gene assay system (Applied Biosystems) and a 2030 ARVO X multilabel counter instrument (PerkinElmer) according to the manufacturers’ procedure. Western blotting was performed as previously described25 by using the following antibodies: anti-p24 polyclonal antibody (ViroStat), anti-KGC antibody (clone 21B10; Medical and Biological Laboratories, Inc.), anti-HA antibody (3F10; Roche) and anti-alpha-Tubulin (TUBA) antibody (DM1A; Sigma).
Statistical analyses
The data expressed as average with standard error (Fig. 3f) or standard deviation (Fig. 4b) and significant differences were determined by Student’s t test (Fig. 4b) or Welch’s t test (Fig. 5c).
Additional Information
How to cite this article: Takeuchi, J. S. et al. Coevolutionary dynamics between tribe Cercopithecini tetherins and their lentiviruses. Sci. Rep. 5, 16021; doi: 10.1038/srep16021 (2015).
References
Keele, B. F. et al. Chimpanzee reservoirs of pandemic and nonpandemic HIV-1. Science 313, 523–526 (2006).
Sharp, P. M. & Hahn, B. H. The evolution of HIV-1 and the origin of AIDS. Philos. Trans. R. Soc. Lond. B Biol. Sci. 365, 2487–2494 (2010).
Worobey, M. et al. Direct evidence of extensive diversity of HIV-1 in Kinshasa by 1960. Nature 455, 661–664 (2008).
Bailes, E. et al. Hybrid origin of SIV in chimpanzees. Science 300, 1713 (2003).
Groves, C. P. Order primates. Mammal species of the world: a taxonomic and geographic reference (eds Wilson D. E, Reeder D. M ), 3rd edn (Johns Hopkins University Press, 2005).
Desrosiers, R. C. Nonhuman lentiviruses. Fields Virology (eds Knipe DM, Howley PM ), 5th edn (Lippincott Williams & Wilkins, Philadelphia, 2007).
Klatt, N. R., Silvestri, G. & Hirsch, V. Nonpathogenic simian immunodeficiency virus infections. Cold Spring Harb. Perspect. Med. 2, a007153 (2012).
Kirchhoff, F. Immune evasion and counteraction of restriction factors by HIV-1 and other primate lentiviruses. Cell Host Microbe 8, 55–67 (2010).
Duggal, N. K. & Emerman, M. Evolutionary conflicts between viruses and restriction factors shape immunity. Nat. Rev. Immunol. 12, 687–695 (2012).
Gifford, R. J. Viral evolution in deep time: lentiviruses and mammals. Trends Genet. 28, 89–100 (2012).
Dawkins, R. & Krebs, J. R. Arms races between and within species. Proc. R. Soc. Lond. B. Biol. Sci. 205, 489–511 (1979).
Sheehy, A. M., Gaddis, N. C., Choi, J. D. & Malim, M. H. Isolation of a human gene that inhibits HIV-1 infection and is suppressed by the viral Vif protein. Nature 418, 646–650 (2002).
Hrecka, K. et al. Vpx relieves inhibition of HIV-1 infection of macrophages mediated by the SAMHD1 protein. Nature 474, 658–661 (2011).
Laguette, N. et al. SAMHD1 is the dendritic- and myeloid-cell-specific HIV-1 restriction factor counteracted by Vpx. Nature 474, 654–657 (2011).
Tristem, M., Marshall, C., Karpas, A. & Hill, F. Evolution of the primate lentiviruses: evidence from vpx and vpr. EMBO J. 11, 3405–3412 (1992).
Lim, E. S. et al. The ability of primate lentiviruses to degrade the monocyte restriction factor SAMHD1 preceded the birth of the viral accessory protein Vpx. Cell Host Microbe 11, 194–204 (2012).
Evans, D. T., Serra-Moreno, R., Singh, R. K. & Guatelli, J. C. BST-2/tetherin: a new component of the innate immune response to enveloped viruses. Trends Microbiol. 18, 388–396 (2010).
Neil, S. J., Zang, T. & Bieniasz, P. D. Tetherin inhibits retrovirus release and is antagonized by HIV-1 Vpu. Nature 451, 425–430 (2008).
Van Damme, N. et al. The interferon-induced protein BST-2 restricts HIV-1 release and is downregulated from the cell surface by the viral Vpu protein. Cell Host Microbe 3, 245–252 (2008).
Sauter, D. et al. Tetherin-driven adaptation of Vpu and Nef function and the evolution of pandemic and nonpandemic HIV-1 strains. Cell Host Microbe 6, 409–421 (2009).
Landi, A., Iannucci, V., Nuffel, A. V., Meuwissen, P. & Verhasselt, B. One protein to rule them all: modulation of cell surface receptors and molecules by HIV Nef. Curr. HIV Res. 9, 496–504 (2011).
Lim, E. S., Malik, H. S. & Emerman, M. Ancient adaptive evolution of tetherin shaped the functions of Vpu and Nef in human immunodeficiency virus and primate lentiviruses. J. Virol. 84, 7124–7134 (2010).
Liu, J., Chen, K., Wang, J. H. & Zhang, C. Molecular evolution of the primate antiviral restriction factor tetherin. PLoS One 5, e11904 (2010).
McNatt, M. W. et al. Species-specific activity of HIV-1 Vpu and positive selection of tetherin transmembrane domain variants. PLoS Pathog. 5, e1000300 (2009).
Kobayashi, T. et al. Characterization of red-capped mangabey tetherin: implication for the co-evolution of primates and their lentiviruses. Sci. Rep. 4, 5529 (2014).
Miceli, M. C. & Parnes, J. R. Role of CD4 and CD8 in T cell activation and differentiation. Adv. Immunol. 53, 59–122 (1993).
Jia, B. et al. Species-specific activity of SIV Nef and HIV-1 Vpu in overcoming restriction by tetherin/BST2. PLoS Pathog. 5, e1000429 (2009).
Kobayashi, T. et al. Identification of amino acids in the human tetherin transmembrane domain responsible for HIV-1 Vpu interaction and susceptibility. J. Virol. 85, 932–945 (2011).
Bibollet-Ruche, F. et al. New simian immunodeficiency virus infecting De Brazza’s monkeys (Cercopithecus neglectus): evidence for a cercopithecus monkey virus clade. J. Virol. 78, 7748–7762 (2004).
Tamura, K., Stecher, G., Peterson, D., Filipski, A. & Kumar, S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 30, 2725–2729 (2013).
Yang, H. et al. Structural insight into the mechanisms of enveloped virus tethering by tetherin. Proc. Natl. Acad. Sci. USA. 107, 18428–18432 (2010).
Courgnaud, V. et al. Identification of a new simian immunodeficiency virus lineage with a vpu gene present among different cercopithecus monkeys (C. mona, C. cephus and C. nictitans) from Cameroon. J. Virol. 77, 12523–12534 (2003).
Sharp, P. M., Shaw, G. M. & Hahn, B. H. Simian immunodeficiency virus infection of chimpanzees. J. Virol. 79, 3891–3902 (2005).
Belshaw, R., Pybus, O. G. & Rambaut, A. The evolution of genome compression and genomic novelty in RNA viruses. Genome Res. 17, 1496–1504 (2007).
Homles, E. C. The Evolution and Emergence of RNA Viruses. Oxford University Press (2009).
Perez-Caballero, D. et al. Tetherin inhibits HIV-1 release by directly tethering virions to cells. Cell 139, 499–511 (2009).
Demogines, A., Abraham, J., Choe, H., Farzan, M. & Sawyer, S. L. Dual host-virus arms races shape an essential housekeeping protein. PLoS Biol. 11, e1001571 (2013).
Perelman, P. et al. A molecular phylogeny of living primates. PLoS Genet. 7, e1001342 (2011).
Skasko, M. et al. HIV-1 Vpu protein antagonizes innate restriction factor BST-2 via lipid-embedded helix-helix interactions. J. Biol. Chem. 287, 58–67 (2012).
Muyembe-Tamfum, J. J., Kipasa, M., Kiyungu, C. & Colebunders, R. Ebola outbreak in Kikwit, democratic republic of the congo: discovery and control measures. J. Infect. Dis. 179 Suppl 1, S259–262 (1999).
Kaletsky, R. L., Francica, J. R., Agrawal-Gamse, C. & Bates, P. Tetherin-mediated restriction of filovirus budding is antagonized by the Ebola glycoprotein. Proc. Natl. Acad. Sci. USA. 106, 2886–2891 (2009).
Le Tortorec, A. & Neil, S. J. Antagonism to and intracellular sequestration of human tetherin by the human immunodeficiency virus type 2 envelope glycoprotein. J. Virol. 83, 11966–11978 (2009).
Gupta, R. K. et al. Simian immunodeficiency virus envelope glycoprotein counteracts tetherin/BST-2/CD317 by intracellular sequestration. Proc. Natl. Acad. Sci. USA. 106, 20889–20894 (2009).
Morrison, J. H. et al. Feline immunodeficiency virus envelope glycoproteins antagonize tetherin through a distinctive mechanism that requires virion incorporation. J. Virol. 88, 3255–3272 (2014).
Yin, X. et al. Equine tetherin blocks retrovirus release and its activity is antagonized by equine infectious anemia virus envelope protein. J. Virol. 88, 1259–1270 (2014).
IUCN 2014. The IUCN red list of threatened species. Version 2014.3. Downloaded on 17 November 2014.
Gifford, R. J. et al. A transitional endogenous lentivirus from the genome of a basal primate and implications for lentivirus evolution. Proc. Natl. Acad. Sci. USA. 105, 20362–20367 (2008).
Gilbert, C., Maxfield, D. G., Goodman, S. M. & Feschotte, C. Parallel germline infiltration of a lentivirus in two Malagasy lemurs. PLoS Genet. 5, e1000425 (2009).
National Research Council. Guide for the care and use of laboratory animals (National Academies Press, Washington, D.C., 1996).
Beer, B. E. et al. Characterization of novel simian immunodeficiency viruses from red-capped mangabeys from Nigeria (SIVrcmNG409 and -NG411). J. Virol. 75, 12014–12027 (2001).
Vinton, C. et al. CD4-like immunological function by CD4- T cells in multiple natural hosts of simian immunodeficiency virus. J. Virol. 85, 8702–8708 (2011).
Saitou, N. & Nei, M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4, 406–425 (1987).
Guindon, S. et al. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst. Biol. 59, 307–321 (2010).
Yang, Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol. Biol. Evol. 24, 1586–1591 (2007).
Kosakovsky Pond, S. L. & Frost, S. D. Not so different after all: a comparison of methods for detecting amino acid sites under selection. Mol. Biol. Evol. 22, 1208–1222 (2005).
Tamura, K. & Nei, M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol. Biol. Evol. 10, 512–526 (1993).
Katoh, K., Misawa, K., Kuma, K. & Miyata, T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 30, 3059–3066 (2002).
Drummond, A. J., Suchard, M. A., Xie, D. & Rambaut, A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol. Biol. Evol. 29, 1969–1973 (2012).
Worobey, M. et al. Island biogeography reveals the deep history of SIV. Science 329, 1487 (2010).
Shibata, R. et al. Construction and characterization of an infectious DNA clone and of mutants of simian immunodeficiency virus isolated from the African green monkey. J. Virol. 64, 307–312 (1990).
Acknowledgements
We would like to thank Dr. Tomoyuki Miura (Institute for Virus Research, Kyoto University, Japan), Dr. Akio Adachi (Institute of Health Biosciences, The University of Tokushima) and Dr. Frank Kirchhoff (University of Ulm, Germany) for providing pSA212 and pCGCG vector, respectively, Dr. Sarah Elliott (University of Pennsylvania) for technical advice, Dr. Akira Kato (Japan Monkey Centre) for sampling and providing the body hair of Sykes’ monkeys and Dr. So Nakagawa (Tokai University) for helpful suggestion for the analyses on the molecular phylogenetic and evolution. This study was supported in-part by the following grants: the Aihara Innovative Mathematical Modeling Project, the Japan Society for the Promotion of Science (JSPS) through the “Funding Program for World-Leading Innovative R&D on Science and Technology (FIRST Program),” initiated by the Council for Science and Technology Policy of Japan (to K.S.); CREST, Japan Science and Technology Agency (to K.S.); Health Labour Sciences Research Grant 26361601 from the Ministry of Health, Labor and Welfare of Japan (MHLW) (to K.S.); Takeda Science Foundation (to K.S.); Sumitomo Foundation Research Grant (to K.S.); Senshin Medical Research Foundation (to K.S.); Imai Memorial Trust for AIDS Research (to K.S.); Ichiro Kanehara Foundation (to K.S.); Kanae Foundation for the Promotion of Medical Science (to K.S.); Suzuken Memorial Foundation (to K.S.); Uehara Memorial Foundation (to K.S.); the intramural research program of NIAID, NIH (to K.M. and V.M.H.); NIH grants R56-AI091516 (to R.G.C.) and R01-MH061139 (to R.G.C.); JSPS Research Fellowship PD 15J06242 (to R.Y.); Grants-in-Aid for Research Activity Start-up (to T.I.) and Scientific Research C15K07166 (to K.S.), B24390112 (to Y.K.) and S22220007 (to Y.K.) from JSPS; a Grant-in-Aid for Scientific Research on Innovative Areas 24115008 (to Y.K.) from the Ministry of Education, Culture, Sports, Science and Technology of Japan; Research on HIV/AIDS from Japan Agency for Medical Research and Development (AMED) 15Afk0410013h0001 (to Y.K.); the Platform Project for Supporting in Drug Discovery and Life Science Research (Platform for Dynamic Approaches to Living System) from AMED (to E.Y.); and JSPS Core-to-Core Program, A. Advanced Research Networks (to K.M., V.M.H. and Y.K.).
Author information
Authors and Affiliations
Contributions
F.R., Y.K. and K.S. wrote the main manuscript text; J.S.T., F.R. and H.T. performed molecular phylogenetic analyses and prepared Figures 1, 2, 3d,f and 5a; T.I. performed homology modeling and prepared Figure 3e; J.S.T., R.Y., E.Y., Y.N., T.K. and N.M. performed the experiments and prepared Figure 4; K.S. performed in silico analyses and prepared Figures 5b,c; K.M., Y.S., K.S.W., R.G.C. and V.M.H. contributed to the monkey samples; K.S. conceived and designed the experiments. All authors reviewed the manuscript.
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Electronic supplementary material
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Takeuchi, J., Ren, F., Yoshikawa, R. et al. Coevolutionary dynamics between tribe Cercopithecini tetherins and their lentiviruses. Sci Rep 5, 16021 (2015). https://doi.org/10.1038/srep16021
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/srep16021
This article is cited by
-
Phylogenetic analysis of small ruminant lentiviruses in Germany and Iran suggests their expansion with domestic sheep
Scientific Reports (2020)
-
New World feline APOBEC3 potently controls inter-genus lentiviral transmission
Retrovirology (2018)
-
Species-specific host factors rather than virus-intrinsic virulence determine primate lentiviral pathogenicity
Nature Communications (2018)
-
A naturally occurring bovine APOBEC3 confers resistance to bovine lentiviruses: implication for the co-evolution of bovids and their lentiviruses
Scientific Reports (2016)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.