- Open Access
New World feline APOBEC3 potently controls inter-genus lentiviral transmission
Retrovirology volume 15, Article number: 31 (2018)
The apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3; A3) gene family appears only in mammalian genomes. Some A3 proteins can be incorporated into progeny virions and inhibit lentiviral replication. In turn, the lentiviral viral infectivity factor (Vif) counteracts the A3-mediated antiviral effect by degrading A3 proteins. Recent investigations have suggested that lentiviral vif genes evolved to combat mammalian APOBEC3 proteins, and have further proposed that the Vif-A3 interaction may help determine the co-evolutionary history of cross-species lentiviral transmission in mammals.
Here we address the co-evolutionary relationship between two New World felids, the puma (Puma concolor) and the bobcat (Lynx rufus), and their lentiviruses, which are designated puma lentiviruses (PLVs). We demonstrate that PLV-A Vif counteracts the antiviral action of APOBEC3Z3 (A3Z3) of both puma and bobcat, whereas PLV-B Vif counteracts only puma A3Z3. The species specificity of PLV-B Vif is irrespective of the phylogenic relationships of feline species in the genera Puma, Lynx and Acinonyx. We reveal that the amino acid at position 178 in the puma and bobcat A3Z3 is exposed on the protein surface and determines the sensitivity to PLV-B Vif-mediated degradation. Moreover, although both the puma and bobcat A3Z3 genes are polymorphic, their sensitivity/resistance to PLV Vif-mediated degradation is conserved.
To the best of our knowledge, this is the first study suggesting that the host A3 protein potently controls inter-genus lentiviral transmission. Our findings provide the first evidence suggesting that the co-evolutionary arms race between lentiviruses and mammals has occurred in the New World.
The apolipoprotein B mRNA editing enzyme catalytic polypeptide-like 3 (APOBEC3; A3) proteins are cellular cytidine deaminases that are specifically encoded in mammals but not in other vertebrates [1, 2]. Mammalian A3 proteins, particularly primate A3 proteins, are considered cellular intrinsic immune factors that potently restrict lentiviral replication. To exhibit antiviral activity, some A3 proteins are incorporated into the released progeny virions and enzymatically insert guanine-to-adenine hypermutations into the viral genome, thereby halting viral replication. In turn, the lentiviral protein, viral infectivity factor (Vif), antagonizes the A3-mediated antiviral action by degrading A3 proteins via the ubiquitin/proteasome-dependent pathway (reviewed in ).
Elucidating co-evolutionary relationships between hosts and their viruses is an intriguing topic in virology and is crucial to understanding how viruses influence their hosts’ evolution and vice versa. Cell-based virological experiments are essential for a better understanding of the evolutionary conflict between mammals and their viruses, including lentiviruses. For instance, since the interaction between host A3 and lentiviral Vif is species specific, various investigations focusing on the functional relationship between A3 and Vif have recently been conducted in combination with molecular phylogenetic approaches (reviewed in ). This strategy stems from the “Red Queen hypothesis ”, which proposes that host and viral proteins have competed with one another for survival over time [6,7,8]. However, since most previous observations are based the Old World evolutionary events , whether the evolutionary arms race between mammals and lentiviruses has occurred in the New World is unclear.
Feline immunodeficiency virus (FIV) is a lentivirus that was first isolated in 1987 from domestic cats (Felis catus) with chronic AIDS-like disorders . Domestic cats encode multiple A3 genes, and the feline A3 protein (designated A3Z3) potently impairs FIV replication by incorporating into nascent virions [10,11,12,13,14,15]. In response, FIV Vif antagonizes the antiviral activity of feline A3Z3 by degrading this protein [11,12,13,14,15]. To elucidate the evolutionary relationship between FIV and felids, we have recently reported that the domestic cat A3Z3 is polymorphic and that a haplotype of the domestic cat A3Z3 is resistant to Vif-mediated degradation by FIVfca, which is the FIV that infects domestic cats . Our findings suggest that the domestic cat A3Z3 is under positive selection due to evolutionary selective pressure caused by FIVfca or related ancestral viruses [12, 14].
In addition to FIVfca in domestic cats, various FIV types have been identified in other felids, such as FIVple in lions (Panthera leo), FIVpco in pumas (Puma concolor), and FIVlru in bobcats (Lynx rufus) (reviewed in ). Interestingly, although these felids become infected with a specific FIV species, Lee et al.  recently reported that a subcluster of FIVpco co-circulated in both pumas and bobcats in North America. The authors designated these viruses (FIVpco and FIVlru) puma lentiviruses (PLVs) and re-classified them as follows: PLV-A includes the FIVpco and FIVlru strains co-circulating in both pumas and bobcats, and PLV-B includes FIVpco strains circulating only in pumas (Fig. 1a) [17,18,19]. This scenario is the first known case indicating cross-species transmission (CST) of a lentivirus between hosts of different genera (Puma and Lynx) at present. Since A3 genes are under positive selection [12, 20,21,22], the A3 sequences are highly variable among hosts of different genera. Furthermore, we can reasonably hypothesize that the specificity in which Vif counteracts the host A3 also differs in each host genus. These findings and insights raise the hypothesis that Vif-A3 interplay between FIV and the two New World felids (puma and bobcat) is closely associated with the mode of CST among lentiviruses, which may illustrate the history of the co-evolutionary arms race between lentiviruses and felids in the New World.
In this study, we perform cell-based virological experiments and demonstrate that the species specificity of the Vif proteins of PLV-A and PLV-B is different between each other. Through the combinational investigations of experimental virology and molecular phylogenetics, we also provide evidence suggesting that bobcat A3Z3 plays a role as the species barrier specifically against PLV-B. Moreover, we determine the amino acid residue responsible for the sensitivity to PLV-B Vif-dependent degradation. This is the first report addressing the co-evolutionary interplay between antiviral A3 of the two New World felids and feline lentiviral Vif.
Bobcat A3Z3 is resistant to PLV-B Vif-mediated degradation
To elucidate the co-evolutionary relationship between PLVs and their host felids, particularly between viral Vif and host A3Z3, first we constructed a phylogenetic tree of FIV vif gene. As shown in Fig. 1a, FIVfca and FIVple formed clusters based on their hosts. In sharp contrast, the FIVlru sequences were co-mingled with some FIVpco sequences, whereas the other FIVpco sequences formed a cluster (Fig. 1a). These results are consistent with the phylogeny of the full-length FIV sequences shown in a previous study . According to this previous report, here we designated the FIVlru and co-mingled FIVpco sequences PLV-A and the other FIVpco sequences PLV-B. The sequence identities of PLV-A Vif and PLV-B Vif were 90.1 ± 4.8 and 85.7 ± 7.1%, respectively, whereas the Vif sequence identity between PLV-A and PLV-B was 40.6 ± 1.1%. Thus, the PLV-A Vif and PLV-B Vif sequences clearly differed.
Next, we focused on the host felid evolutionary relationships. Based on comprehensive genetic information, a previous study  described the evolutionary history of the felids of interest in this study (Fig. 1b). The common ancestor of the puma (Puma concolor) and the bobcat (Lynx rufus) is estimated to have crossed the Bering Isthmus from Eurasia to the New World approximately 8.0–8.5 million years ago (Mya) and to have diverged approximately 7.2 Mya . The puma diverged from the cheetah (Acinonyx jubatus) approximately 4.9 Mya, whereas the bobcat was branched from the lynx (Lynx lynx) approximately 3.2 Mya (Fig. 1b) . Puma and bobcat resided in the New World, whereas lynx and cheetah returned to Eurasia approximately 1.2–1.6 Mya (Fig. 1b) . To assess the phylogenetic relationships of the A3Z3 genes of these felids, we collected body hairs from a bobcat and a cheetah from Japanese zoos and determined their A3Z3 sequences (see Additional file 1: Fig. S1 and Additional file 2: Table S1). As shown in Fig. 1c, the topology of the phylogenetic tree of the A3Z3 genes form these four felids corresponded to their evolutionary relationships.
Based on these observations, two scenarios are assumed. The first scenario was that both PLV-A and PLV-B co-circulated in pumas and bobcats in the New World but PLV-B became extinct only in the bobcat population (“scenario 1” in Fig. 1d). The other possibility was that PLV-A and PLV-B were specific viruses in bobcats and pumas, respectively, in the past. After co-habitation of these two feline species in the New World, PLV-A was transferred from bobcats to pumas, whereas PLV-B CST from pumas to bobcats was hampered for unknown reasons (“scenario 2” in Fig. 1d).
To experimentally investigate the interplay between PLVs and their hosts, we analyzed the antiviral activity of the puma and bobcat A3Z3 proteins. In the absence of PLV Vif, both puma and bobcat A3Z3 were expressed at comparable levels and are incorporated into the released particles of vif-deficient FIV (strain Petaluma) in a dose-dependent manner (Fig. 2a). In addition, the infectivity of vif-deficient FIV was suppressed by both puma and bobcat A3Z3 proteins in a dose-dependent manner (Fig. 2b). Interestingly, although the protein expression (Fig. 2a, top) and incorporation levels in the released virions (Fig. 2a, bottom) were similar between the puma and bobcat A3Z3 s, the antiviral effect was significantly higher for the puma A3Z3 than for the bobcat A3Z3 (Fig. 2b).
Next, we analyzed the sensitivity of puma and bobcat A3Z3 to PLV Vif proteins. We used expression plasmids for three PLV-A Vif strains (Lru1, Lru20 and Pco5) and three PLV-B Vif strains (Pco14, Pco16 and Pco27). The sampling years and locations for each virus are summarized in Table 1. As shown in Fig. 2c, the puma A3Z3 was degraded by all of the PLV Vif proteins tested in this study, and the antiviral effect of puma A3Z3 was counteracted by these PLV Vif proteins (Fig. 2d). In sharp contrast, the bobcat A3Z3 was degraded by the PLV-A Vif proteins but was resistant to degradation mediated by the PLV-B Vif proteins (Fig. 2e, top). Bobcat A3Z3 was incorporated into nascent virions even in the presence of PLV-B Vif (Fig. 2e, bottom), which significantly decreased viral infectivity (Fig. 2f). Taken together, these findings suggest that bobcat A3Z3 is resistant to PLV-B Vif-mediated degradation.
Puma A3Z3 is uniquely sensitive to PLV-B Vif-mediated degradation
As shown in Fig. 1b, the puma and bobcat are evolutionarily similar to the cheetah (Acinonyx jubatus) and the lynx (Lynx lynx), respectively. To address whether the sensitivity of host A3Z3 to PLV-B Vif is associated with the host phylogeny, we constructed expression plasmids for cheetah and lynx A3Z3 proteins and performed similar experiments. Our results revealed that both cheetah and lynx A3Z3 were degraded by PLV-A Vif but were resistant to PLV-B Vif-mediated degradation (Fig. 3a, b, top). Similar to bobcat A3Z3, the cheetah and lynx A3Z3 proteins were incorporated into the released virions even in the presence of PLV-B Vif (Fig. 3a, b, bottom) resulting in an antiviral effect (Fig. 3c, d). These findings suggest that the PLV-A Vif antagonizes all four feline A3Z3 tested in this study, whereas PLV-B Vif antagonizes only the puma A3Z3.
The threonine residue at position 178 of bobcat A3Z3 confers the resistance to PLV-B Vif-mediated degradation
We compared the amino acid sequence of puma, bobcat, cheetah and lynx A3Z3s and found three amino acids at positions 90, 131 and 178 that consistently differed between the puma and the other felids (Fig. 4a). To assess the positions of these three residues in the tertiary structures of the puma and bobcat A3Z3 proteins, we constructed A3Z3 protein homology models for these felids (Fig. 4b). Residues 90 and 178 are located in alpha-helices 3 and 6, respectively, whereas residue 131 is positioned in the loop (Fig. 4c). Additionally, residues 90 and 178 are localized on the protein surface (Fig. 4d).
To determine the amino acid residue(s) responsible for the resistance to PLV-B Vif-mediated degradation, we constructed lines of bobcat A3Z3 mutants and performed cell-based loss-of-function experiments. Although the bobcat A3Z3 E90K and I131T mutants were resistant to PLV-B Vif-mediated degradation (strain Pco27), the bobcat A3Z3 T178M mutant was degraded by PLV-B Vif (Fig. 5a). We also used combination mutants and found those harboring the T178M substitution lost the ability to resist PLV-B Vif-mediated degradation (strain Pco27) (Fig. 5a). In the presence of PLV-B Vif, incorporation of bobcat A3Z3 derivatives possessing the T178M mutation into the released viruses was impaired (Fig. 5a, bottom), and the derivatives’ antiviral effects were abrogated (Fig. 5b). To determine whether this observation was strain-specific, we performed similar experiments using the other PLV-B strains (Pco14 and Pco16). Similar to strain Pco27 (Fig. 5a, b), the Vif proteins of the other PLV-B strains degraded the bobcat A3Z3 T178M mutant (Fig. 5c) and significantly recovered viral infectivity (Fig. 5d). These findings suggest that the resistance of bobcat A3Z3 to PLV-B Vif-mediated degradation is determined by the amino acid residue at position 178.
To validate the importance of the amino acid at position 178 to sensitivity to PLV-B Vif-mediated degradation, we performed gain-of-function experiments based on the puma A3Z3. By substituting the methionine residue at position 178 of the puma A3Z3 to threonine, the puma A3Z3 mutant became resistant to PLV-B Vif-mediated degradation (Fig. 6a, top). The puma A3Z3 M178T mutant was efficiently incorporated into the released virions (Fig. 6a, bottom) and exhibited significant antiviral activity (Fig. 6b). Taken together, these findings suggest that the amino acid at position 178 in the bobcat/puma A3Z3 plays a pivotal role in conferring resistance to PLV-B Vif-mediated degradation.
Puma and bobcat A3Z3 are polymorphic
Previous studies have revealed that mammalian A3 genes, including feline A3Z3 genes, are under positive selection and are highly diversified due to evolutionary selective pressures presumably caused by the ancestral FIV Vif . For instance, the domestic cat A3Z3 is polymorphic , and a haplotype of the domestic cat A3Z3 renders resistance to FIVfca Vif-mediated degradation . These previous findings raise the possibility that the puma and bobcat A3Z3 genes are also polymorphic and that the PLV Vif sensitivity may differ among haplotypes. To address this possibility, we additionally determined the A3Z3 sequences of five pumas, one bobcat as well as one cheetah. As shown in Fig. 7a, we detected additional puma and bobcat A3Z3 haplotypes and designated them haplotype II (hap II). One nonsynonymous mutation was detected in the puma A3Z3 hap II (a c392t mutation resulting in T131I amino acid substitution), and a heterozygous sequence was detected in the bobcat A3Z3 hap II (g347r; a g347a mutation resulting in a R116H amino acid substitution) (Fig. 7a). Next, we mapped these residues on the protein homology model. Both residue 131 of the puma A3Z3 (Fig. 4b) and residue 116 of the bobcat A3Z3 (Additional file 3: Fig. S2) are positioned in the loop structure. Then, we prepared expression plasmids for these haplotypes and assessed their antiviral activity and sensitivity to PLV Vif. As shown in Fig. 7b, puma and bobcat A3Z3 hap II were expressed at similar levels to their hap I derivatives. Notably, in the absence of PLV Vif, the antiviral effects of the puma and bobcat A3Z3 hap II were significantly higher than those of the hap I derivatives (Fig. 7c). In particular, the bobcat A3Z3 hap II exhibited stronger antiviral activity in a dose-dependent manner (Fig. 7c). Moreover, puma A3Z3 hap II was sensitive to both the PLV-A and PLV-B Vifs (Fig. 7d, e), whereas the bobcat A3Z3 hap II was sensitive to PLV-A Vif but was resistant to PLV-B Vif-mediated degradation (Fig. 7e, f). Taken together, these findings suggest that the puma and bobcat A3Z3 genes are polymorphic and that the bobcat A3Z3 hap II exhibits a higher antiviral effect than the hap I derivative, whereas the PLV Vif sensitivity phenotype is conserved within each species.
In this study, we identified the A3Z3 sequences of two felids, the bobcat and cheetah, and demonstrated that PLV-A Vif counteracted the A3Z3 proteins of both the puma and bobcat lineages, whereas PLV-B Vif counteracted only the puma A3Z3 (Figs. 1, 2, 3). Through loss-of-function (Fig. 5) and gain-of-function (Fig. 6) experiments, we also determined that the amino acid at position 178 of the puma and bobcat A3Z3 was responsible for the sensitivity to PLV-B Vif. In addition, structural modeling suggested that this residue was exposed to the protein surface (Fig. 4). Furthermore, although the puma and bobcat A3Z3 genes were polymorphic, the sensitivity to PLV-B Vif was conserved in each species (Fig. 7). Some previous studies have addressed the potential of OWM A3G to restrict HIV-1 infection [25, 26]. However, OWMs do not infect with HIV-1 in nature. Therefore, to the best of our knowledge, this study is the first to suggest that inter-species lentiviral transmission in nature is controlled by a host A3 protein.
As summarized in Fig. 1a, d, PLV-A is shared by both pumas and bobcats in North America, whereas PLV-B is detected only in pumas. Through cell-based virological experiments, here we demonstrated that PLV-A Vif degraded the A3Z3 protein of feline lineages including puma and bobcat, whereas PLV-B Vif counteracted only the puma A3Z3. Additionally, the bobcat A3Z3 exhibited a significant antiviral effect even in the presence of PLV-B Vif. Because the functionality of PLV Vif is independent of the sampling location and year (Fig. 2 and Table 1), our results suggest that the functional relationship between PLV Vif and feline A3Z3 is conserved. Moreover, we found a polymorphism in the puma A3Z3, but the sensitivity of the puma A3Z3 variants to PLV-B Vif-mediated degradation was also conserved (Fig. 6). These findings suggest that “scenario 2” reasonably explains the difference in species tropism between PLV-A and PLV-B and that bobcat A3Z3 acts as a species barrier to restrict PLV-B CST of PLV-B from pumas (Fig. 1d). To the best of our knowledge, this is the first report demonstrating that the host A3 protein is the factor that restricts non-primate lentiviral CST.
Consistent with our results, PLV-A and PLV-B form clusters by sampling location [17,18,19]. Moreover, the PLV-As detected in pumas (FIVpco) and bobcats (FIVlru), particularly those isolated in California, were highly co-mingled in a cluster (Fig. 1a). These observations suggest that PLV-As have been transmitted relatively recently from bobcats to pumas after geological sequestration, and we can assume that this bobcat-to-puma CST event occurs frequently at least in California.
Our findings in New World felids were reminiscent of the fact that simian immunodeficiency viruses (SIVs) in Old World monkeys (OWMs) were transferred to chimpanzees, which originated in the emergence of a novel SIV in chimpanzees (SIVcpz) (reviewed in [27, 28]). SIV CST from OWMs to chimpanzees in Africa is explained phylogenetically and is assumed to be due to the prey-predator relationship between OWMs and chimpanzees. Small OWMs are the prey of chimpanzees in the wild; thus, chimpanzees are frequently exposed to various SIVs that infect their OWM prey species [29, 30]. Similarly, pumas are a top carnivore in the New World  and bobcats are their prey [32, 33]. Therefore, pumas may be frequently exposed to PLV-A in bobcats via prey-predator interactions, leading to bobcat-to-puma CST; this possibility is consistent with a previous assumption based on phylogenetic analyses .
Although PLV-A Vif degraded all of the feline A3Z3 proteins tested in this study, PLV-B Vif counteracted only the puma A3Z3, suggesting the evolutionary convergence of the puma and PLV-B. Compton and Emerman reported that the Vif protein of an SIV infecting colobus monkey (Colobus guereza), which is an OWM in Africa, counteracted only colobus A3G . Based on findings for the A3G-Vif interaction in colobus monkeys, the authors proposed convergent evolution of colobus monkeys and SIVcor and further suggested that the convergent co-evolution of these monkeys and viruses occurred at least 12 Mya . Similarly, the puma diversified from the cheetah approximately 4.9 Mya (Fig. 1c) . Puma and cheetah habitats were separated approximately 1.2–1.6 Mya after cheetahs immigrated from the New World to Eurasia (Fig. 1b) . Since the cheetah A3Z3 was not counteracted by the PLV-B Vif (Fig. 3), our findings suggested that the co-evolutionary convergence of the puma and PLV-B arose at least 1.2–1.6 Mya. Moreover, OWMs exhibit higher viral RNA loads without any disorders . Since SIVs circumvent the immune pressure from host monkeys because they are naturally infected, SIVs are less diversified in African OWMs, which is assumed to be due to the co-existence of SIVs in OWMs over a long period . Lee et al.  recently showed that the PLV-B plasma viral RNA load in pumas was significantly higher than PLV-A loads in both pumas and bobcats, whereas PLV-B was less diversified than PLV-A. This is the first study suggesting the convergent co-evolution of a mammal (puma) and a lentivirus (PLV-B) in the New World, and PLV-B can be concluded to have co-evolved with pumas over a long period, similar to apathogenic SIV infections in OWMs.
Through mutagenesis experiments, we demonstrated that the amino acid at position 178 determined the sensitivity to PLV-B Vif-dependent degradation (Fig. 5). Since both the gain- and loss-of-function experiments indicated the importance of this residue, the species specificity of PLV-B Vif against puma and bobcat A3Z3 was determined by this single amino acid. In addition, the structural biology investigation revealed that this amino acid residue was exposed on the A3Z3 protein surface (Fig. 4). In our previous study , we revealed that there are at least seven haplotypes in the domestic cat A3Z3 and that a domestic cat A3Z3 haplotype (designated haplotype V in our previous study) was resistant to FIVfca Vif-mediated degradation. Furthermore, we demonstrated that amino acid residue 65 of the domestic cat A3Z3 determines the sensitivity to FIVfca Vif . Similar to this study, residue 65 in the domestic cat A3Z3 was also exposed on the protein surface, which suggested that this residue was associated with the interaction with FIV Vif . However, the amino acid residue at position 65 was positioned at alpha-helix 2 in the domestic cat A3Z3 , whereas the residue 178 was located at alpha-helix 6 in the puma and bobcat A3Z3 (Fig. 5). Therefore, our findings suggest that these residues render resistance to FIV Vif in a different manner.
In addition to A3, other cellular proteins such as tripartite motif containing protein 5 (TRIM5)  and tetherin [37, 38] are known as restriction factors that inhibit lentiviral replication in primates. Therefore, it might be possible that these restriction factors in pumas and/or bobcats may also play roles in preventing FIV CST. However, feline TRIM5 is not functional because of the insertion of premature stop codon , and feline tetherin is unable to restrict spreading FIV infection . Furthermore, although FIV Vif is a functional antagonist against feline A3, FIV lacks functional proteins that potently counteract TRIM5 and tetherin [39, 40]. Primates encode several restriction factors against lentiviruses; however, only A3 is convincing to be a restriction factor in felids. Therefore, feline A3 may play a pivotal role in controlling CST in felids.
In summary, we demonstrated that bobcat A3Z3 was resistant to PLV-B Vif-dependent degradation, whereas PLV-A Vif overcame the antiviral action mediated by both the puma and bobcat A3Z3. The co-evolutionary relationship between primate A3 proteins and their lentiviral Vifs has been rigorously investigated [20, 22, 34, 41,42,43,44]; however, few studies, including ours, have addressed the evolutionary dynamics of non-primate A3 and non-primate lentiviruses [4, 12, 14]. Here, we provided evidence suggesting that lentiviral CST between different genera (Puma and Lynx) was controlled by the Vif-A3 interaction. To the best of our knowledge, this is the first report providing evidence of the co-evolutionary arms race between mammals and lentiviruses in the New World.
To determine the feline A3Z3 sequences, blood, body hair or cryopreserved muscle tissue of pumas, bobcats and cheetahs were kindly provided by the following facilities: Tennoji zoo, Osaka, Japan; Kobe City Oji zoo, Hyogo, Japan; Omoriyama zoo, Akita, Japan; Tama Zoological Park, Tokyo, Japan; Izu Animal Kingdom, Shizuoka, Japan; and Shizuoka Municipal Nihondaira zoo, Shizuoka, Japan. Sampling was performed in accordance with the guideline of Tokyo University of Agriculture, Japan. All experimental protocols were approved by a committee at Tokyo University of Agriculture, Japan.
Sequencing PCR of the feline A3Z3 genes
Sequencing PCR of feline A3Z3 was performed as previously described [45, 46]. Briefly, genomic DNA was extracted from the samples described above using the DNA Extractor FM kit (Wako) or the DNeasy Blood & Tissue it (Qiagen). PCR was performed using the PrimeSTAR GXL DNA polymerase (Takara) and the primers are listed in Additional file 2: Table S1. The obtained PCR products were purified by gel extraction using the QIAquick gel extraction kit (Qiagen). The nucleotide sequences were determined by a DNA sequencing service (Fasmac, Kanagawa, Japan) and the data were analyzed using Sequencher v5.1 software (Gene Codes Corporation).
Molecular phylogenetic analysis of the FIV vif and feline A3Z3 genes
Molecular phylogenetic analyses were performed as previously described [12, 14, 45,46,47]. Briefly, the sequences of PLV Vif and feline A3Z3, some of which were newly identified in this study (Figs. 1a, c, 7a), were aligned using ClustalW implemented in MEGA7 . The alignments were verified manually at the amino acid level. Phylogenetic trees (Figs. 1a, c, 7a) were reconstructed using the maximum likelihood (ML) method with PhyML . To calculate the amino acid sequence diversity of PLV Vif, 26 and 34 amino acid sequences from PLV-A and PLV-B Vif (summarized in Fig. 1a) are used, respectively. A multiple alignment was generated using L-INS-i in MAFFT . The gapped regions were removed using trimAl with the nogaps option , and 227 amino acid sites were used for the analysis. Then, we performed pairwise comparisons of 60 amino acid sequences to calculate the amino acid sequence identity using MEGA7 .
The expressing plasmids for HA-tagged puma and lynx A3Z3 were kindly provided by Dr. Carsten Münk . The expressing plasmids for HA-tagged bobcat and cheetah A3Z3 were constructed by using the genomic DNA fragments and the primers listed in Additional file 2: Table S1. The point mutants of HA-tagged puma and bobcat A3Z3 were constructed by using a GeneArt site-directed mutagenesis system (Thermo Fisher Scientific). Each wild-type plasmid was used as the template, and the primers used are listed in Additional file 2: Table S1. The His-tagged PLV-A Vif (strains Lru1, Lru20 and Pco5) and PLV-B (strains Pco14, Pco16 and Pco27) were obtained from GeneArt gene synthesis service (Thermo Fisher Scientific). The obtained DNA fragments were digested with BamHI and SalI and inserted into the BamHI-SalI site of pDON-AI plasmid (Takara). The nucleotide sequences were determined by a DNA sequencing service (Fasmac, Kanagawa, Japan) and the data were analyzed by Sequencher v5.1 software (Gene Codes Corporation).
Cell culture and transfection
HEK293T cells (CRL-11268; ATCC) were cultured in Dulbecco’s modified Eagle medium (Sigma-Aldrich) supplemented with 10% heat-inactivated fetal calf serum and antibiotics (Thermo Fisher Scientific). Transfection was performed by using PEI Max (GE Healthcare) in accordance with the manufacturers procedures and described previously [12, 14, 45,46,47, 52]. To analyze the dose-dependent anti-FIV activity of feline A3Z3, pFP93 (pFIVgagpolΔvif; a replication incompetent vif-deficient FIV packaging construct derived from clone FIV-34TF10 [GenBank accession number M25381]; kindly provided by Dr. Eric M. Poeschla) (200 ng), pTiger-luc (pFIVΨ-luc) (150 ng), and pMD.G (pVSVg; a vesicular stomatitis virus G [VSVg] expression plasmid) (50 ng) were co-transfected into HEK293T cells (1 × 105 cells) with feline A3Z3 expression plasmid (50, 100, 200, or 400 ng). To analyze the functional relationship between feline A3Z3 and PLV-Vif, feline A3Z3 expression plasmid (200 ng), pFP93 (200 ng), pTiger-luc (150 ng), and pMD.G (50 ng) were co-transfected into HEK293T cells with or without His-tagged PLV Vif expression plasmid (400 ng). After 48 h post-transfection, the transfected cells and culture supernatants were harvested as previously described [12,13,14, 45,46,47, 52].
Western blotting and virus reporter assay
Western blotting and reporter assay were performed as previously described [12,13,14, 45,46,47, 52]. For the Western blotting of virus particles, 340 μl of the culture supernatant was ultracentrifuged at 100,000×g for 1 h at 4 °C using a TL-100 instrument (Beckman), and the pellet was lysed with 1 × SDS buffer. For the Western blotting of transfected cells, the cells were lysed with RIPA buffer (50 mM Tris–HCl buffer [pH 7.6], 150 mM NaCl, 1% Nonidet P-40, 0.5% sodium deoxycholate, 0.1% SDS) with protease inhibitor cocktail (Roche). The following antibodies for Western blotting: anti-His polyclonal antibody (OGHis; Medical and Biological Laboratories), anti-HA antibody (3F10; Roche), anti-FIV p24 Capsid antibody (PAK3-2C1; Santa Cruz Biotechnology); anti-alpha-tubulin (TUBA) antibody (DM1A; Sigma), and anti-VSVg antibody (P5DA; Roche). For FIV reporter assays, HEK293T cells were used for the target of infection. Ten microliter of the culture supernatant of transfected cells was inoculated into HEK293T cells in a 96-well plate (Nunc), and the firefly luciferase activity was measured by using the BrillianStar-LT assay system (Toyo-b-net) and the 2030 ARVO X multilabel counter instrument (PerkinElmer) according to the manufacturers’ procedures.
Protein homology modeling
Homology modeling was performed using the SWISS-MODEL server [53,54,55,56]. After BLAST searches  of the bobcat and puma A3Z3 amino acid sequences against protein data bank sequence entries (http://www.rcsb.org/pdb/), the crystal structure of human APOBEC3A (PDB: 5KEG)  was selected as the best template for homology modeling per the Global Model Quality Estimations, QMEAN statistical parameters, and modeled sequence length. Each generated model was minimized and refined using Discovery Studio (Dassault Systèmes BIOVIA, Discovery Studio Modeling Environment, Release 4.1, San Diego: Dassault Systèmes, 2007). Mutant models and 3D images were generated with PyMOL (The PyMOL Molecular Graphics System, version 1.8 Schrödinger, LLC).
The data are expressed as averages with the standard deviations (SDs), and statistically significant differences were determined using Student’s t test.
The A3Z3 sequences of five pumas, two bobcats and two cheetahs were submitted to DDBJ (accession numbers LC376039-LC376042).
apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3)
feline immunodeficiency virus
million years ago
Old World monkey
simian immunodeficiency virus
tripartite motif containing protein 5
viral infectivity factor
vesicular stomatitis virus G
Koito A, Ikeda T. Intrinsic restriction activity by AID/APOBEC family of enzymes against the mobility of retroelements. Mob Genet Elem. 2011;1:197–202.
LaRue RS, Andresdottir V, Blanchard Y, Conticello SG, Derse D, Emerman M, et al. Guidelines for naming nonprimate APOBEC3 genes and proteins. J Virol. 2009;83:494–7.
Harris RS, Dudley JP. APOBECs and virus restriction. Virology. 2015;479–480:131–45.
Nakano Y, Aso H, Soper A, Yamada E, Moriwaki M, Juarez-Fernandez G, et al. A conflict of interest: the evolutionary arms race between mammalian APOBEC3 and lentiviral Vif. Retrovirology. 2017;14:31.
Dawkins R, Krebs JR. Arms races between and within species. Proc R Soc Lond B Biol Sci. 1979;205:489–511.
Duggal NK, Emerman M. Evolutionary conflicts between viruses and restriction factors shape immunity. Nat Rev Immunol. 2012;12:687–95.
Gifford RJ. Viral evolution in deep time: lentiviruses and mammals. Trends Genet. 2012;28:89–100.
Kirchhoff F. Immune evasion and counteraction of restriction factors by HIV-1 and other primate lentiviruses. Cell Host Microbe. 2010;8:55–67.
Pedersen NC, Ho EW, Brown ML, Yamamoto JK. Isolation of a T-lymphotropic virus from domestic cats with an immunodeficiency-like syndrome. Science. 1987;235:790–3.
Münk C, Beck T, Zielonka J, Hotz-Wagenblatt A, Chareza S, Battenberg M, et al. Functions, structure, and read-through alternative splicing of feline APOBEC3 genes. Genome Biol. 2008;9:R48.
Wang J, Zhang W, Lv M, Zuo T, Kong W, Yu X. Identification of a Cullin5-ElonginB-ElonginC E3 complex in degradation of feline immunodeficiency virus Vif-mediated feline APOBEC3 proteins. J Virol. 2011;85:12482–91.
Yoshikawa R, Izumi T, Yamada E, Nakano Y, Misawa N, Ren F, et al. A naturally occurring domestic cat APOBEC3 variant confers resistance to FIV infection. J Virol. 2015;90:474–85.
Yoshikawa R, Nakano Y, Yamada E, Izumi T, Misawa N, Koyanagi Y, et al. Species-specific differences in the ability of feline lentiviral Vif to degrade feline APOBEC3 proteins. Microbiol Immunol. 2016;60:272–9.
Yoshikawa R, Takeuchi JS, Yamada E, Nakano Y, Misawa N, Kimura Y, et al. Feline immunodeficiency virus evolutionarily acquires two proteins, vif and protease, capable of antagonizing Feline APOBEC3. J Virol. 2017;91:e00250-17.
Zielonka J, Marino D, Hofmann H, Yuhki N, Löchelt M, Münk C. Vif of feline immunodeficiency virus from domestic cats protects against APOBEC3 restriction factors from many felids. J Virol. 2010;84:7312–24.
Pecon-Slattery J, Troyer JL, Johnson WE, O’Brien SJ. Evolution of feline immunodeficiency virus in Felidae: implications for human health and wildlife ecology. Vet Immunol Immunopathol. 2008;123:32–44.
Lee J, Malmberg JL, Wood BA, Hladky S, Troyer R, Roelke M, et al. Feline immunodeficiency virus cross-species transmission: implications for emergence of new lentiviral infections. J Virol. 2017;91:e02134-16.
Franklin SP, Troyer JL, Terwee JA, Lyren LM, Boyce WM, Riley SP, et al. Frequent transmission of immunodeficiency viruses among bobcats and pumas. J Virol. 2007;81:10961–9.
Lee JS, Bevins SN, Serieys LE, Vickers W, Logan KA, Aldredge M, et al. Evolution of puma lentivirus in bobcats (Lynx rufus) and mountain lions (Puma concolor) in North America. J Virol. 2014;88:7727–37.
Sawyer SL, Emerman M, Malik HS. Ancient adaptive evolution of the primate antiviral DNA-editing enzyme APOBEC3G. PLoS Biol. 2004;2:E275.
Vallender EJ, Lahn BT. Positive selection on the human genome. Hum Mol Genet. 2004;13(2):R245–54.
Krupp A, McCarthy KR, Ooms M, Letko M, Morgan JS, Simon V, et al. APOBEC3G polymorphism as a selective barrier to cross-species transmission and emergence of pathogenic SIV and AIDS in a primate host. PLoS Pathog. 2013;9:e1003641.
Johnson WE, Eizirik E, Pecon-Slattery J, Murphy WJ, Antunes A, Teeling E, et al. The late Miocene radiation of modern Felidae: a genetic assessment. Science. 2006;311:73–7.
de Castro FL, Junqueira DM, de Medeiros RM, da Silva TR, Costenaro JG, Knak MB, et al. Analysis of single-nucleotide polymorphisms in the APOBEC3H gene of domestic cats (Felis catus) and their association with the susceptibility to feline immunodeficiency virus and feline leukemia virus infections. Infect Genet Evol. 2014;27:389–94.
Hatziioannou T, Princiotta M, Piatak M Jr, Yuan F, Zhang F, Lifson JD, et al. Generation of simian-tropic HIV-1 by restriction factor evasion. Science. 2006;314:95.
Schrofelbauer B, Senger T, Manning G, Landau NR. Mutational alteration of human immunodeficiency virus type 1 Vif allows for functional interaction with nonhuman primate APOBEC3G. J Virol. 2006;80:5984–91.
Hahn BH, Shaw GM, De Cock KM, Sharp PM. AIDS as a zoonosis: scientific and public health implications. Science. 2000;287:607–14.
Sharp PM, Hahn BH. The evolution of HIV-1 and the origin of AIDS. Philos Trans R Soc Lond B Biol Sci. 2010;365:2487–94.
Gogarten JF, Akoua-Koffi C, Calvignac-Spencer S, Leendertz SA, Weiss S, Couacy-Hymann E, et al. The ecology of primate retroviruses: an assessment of 12 years of retroviral studies in the Tai national park area, Cote d Ivoire. Virology. 2014;460–461:147–53.
Leendertz SA, Locatelli S, Boesch C, Kucherer C, Formenty P, Liegeois F, et al. No evidence for transmission of SIVwrc from western red colobus monkeys (Piliocolobus badius badius) to wild West African chimpanzees (Pan troglodytes verus) despite high exposure through hunting. BMC Microbiol. 2011;11:24.
Lomolino MV, Riddle BR, Whittaker RJ, Brown JH: Distributions of Communities. In: Biogeography. 4th ed. Sunderland: Oxford University Press; 2010. pp. 121–164.
Wynne EM, Mathew WA, Jonathan NP. Quantifying risk and resource use for a large carnivore in an expanding urban-wildland interface. J Appl Ecol. 2015;53:371–8.
Smith JA, Wang Y, Wilmers CC. Top carnivores increase their kill rates on prey as a response to human-induced fear. Proc Biol Sci. 2015;282(1802):20142711.
Compton AA, Emerman M. Convergence and divergence in the evolution of the APOBEC3G-Vif interaction reveal ancient origins of simian immunodeficiency viruses. PLoS Pathog. 2013;9:e1003135.
Coffin JM. HIV population dynamics in vivo: implications for genetic variation, pathogenesis, and therapy. Science. 1995;267:483–9.
Stremlau M, Owens CM, Perron MJ, Kiessling M, Autissier P, Sodroski J. The cytoplasmic body component TRIM5alpha restricts HIV-1 infection in Old World monkeys. Nature. 2004;427:848–53.
Neil SJ, Zang T, Bieniasz PD. Tetherin inhibits retrovirus release and is antagonized by HIV-1 Vpu. Nature. 2008;451:425–30.
Van Damme N, Goff D, Katsura C, Jorgenson RL, Mitchell R, Johnson MC, 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. 2008;3:245–52.
McEwan WA, Schaller T, Ylinen LM, Hosie MJ, Towers GJ, Willett BJ. Truncation of TRIM5 in the Feliformia explains the absence of retroviral restriction in cells of the domestic cat. J Virol. 2009;83:8270–5.
Dietrich I, McMonagle EL, Petit SJ, Vijayakrishnan S, Logan N, Chan CN, et al. Feline tetherin efficiently restricts release of feline immunodeficiency virus but not spreading of infection. J Virol. 2011;85:5840–52.
Compton AA, Hirsch VM, Emerman M. The host restriction factor APOBEC3G and retroviral Vif protein coevolve due to ongoing genetic conflict. Cell Host Microbe. 2012;11:91–8.
Compton AA, Malik HS, Emerman M. Host gene evolution traces the evolutionary history of ancient primate lentiviruses. Philos Trans R Soc Lond B Biol Sci. 2013;368:20120496.
Jern P, Coffin JM. Host-retrovirus arms race: trimming the budget. Cell Host Microbe. 2008;4:196–7.
OhAinle M, Kerns JA, Li MM, Malik HS, Emerman M. Antiretroelement activity of APOBEC3H was lost twice in recent human evolution. Cell Host Microbe. 2008;4:249–59.
Kobayashi T, Takeuchi JS, Ren F, Matsuda K, Sato K, Kimura Y, et al. Characterization of red-capped mangabey tetherin: implication for the co-evolution of primates and their lentiviruses. Sci Rep. 2014;4:5529.
Yamada E, Yoshikawa R, Nakano Y, Misawa N, Kobayashi T, Ren F, et al. A naturally occurring bovine APOBEC3 confers resistance to bovine lentiviruses: implication for the co-evolution of bovids and their lentiviruses. Sci Rep. 2016;6:33988.
Takeuchi JS, Ren F, Yoshikawa R, Yamada E, Nakano Y, Kobayashi T, et al. Coevolutionary dynamics between tribe Cercopithecini tetherins and their lentiviruses. Sci Rep. 2015;5:16021.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.
Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59:307–21.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.
Capella-Gutierrez S, Silla-Martinez JM, Gabaldon T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009;25:1972–3.
Nakano Y, Misawa N, Juarez-Fernandez G, Moriwaki M, Nakaoka S, Funo T, et al. HIV-1 competition experiments in humanized mice show that APOBEC3H imposes selective pressure and promotes virus adaptation. PLoS Pathog. 2017;13:e1006348.
Biasini M, Bienert S, Waterhouse A, Arnold K, Studer G, Schmidt T, et al. SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information. Nucleic Acids Res. 2014;42:W252–8.
Guex N, Peitsch MC, Schwede T. Automated comparative protein structure modeling with SWISS-MODEL and Swiss-PdbViewer: a historical perspective. Electrophoresis. 2009;30(Suppl 1):S162–73.
Bordoli L, Kiefer F, Arnold K, Benkert P, Battey J, Schwede T. Protein structure homology modeling using SWISS-MODEL workspace. Nat Protoc. 2009;4:1–13.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–402.
Arnold K, Bordoli L, Kopp J, Schwede T. The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics. 2006;22:195–201.
Kouno T, Silvas TV, Hilbert BJ, Shandilya SMD, Bohn MF, Kelch BA, et al. Crystal structure of APOBEC3A bound to single-stranded DNA reveals structural basis for cytidine deamination and specificity. Nat Commun. 2017;8:15024.
Bruen TC, Poss M. Recombination in feline immunodeficiency virus genomes from naturally infected cougars. Virology. 2007;364:362–70.
KS conceived the study and designed experiments; Y.Konno, S.Nagaoka, IK, KY, Y.Kagawa, RK, HA and TK performed experiments; Y.Konno, S.Nagaoka, IK, KY and HA analyzed data; MTU and S.Nakagawa constructed protein homology models; S.Nakagawa and KS performed phylogenetic analyses; KS wrote the paper; and Y.Konno, S.Nagaoka, IK, MTU, S.Nakagawa, Y.Koyanagi and KS edited the manuscript. All authors read and approved the final manuscript.
We thank Carsten Münk (Heinrich-Heine-Universität, Düsseldorf, Germany) and Eric M. Poeschla (University of Colorado, USA) for providing experimental materials and Tennoji zoo (Osaka, Japan), Kobe City Oji zoo (Hyogo, Japan), Omoriyama zoo (Akita, Japan), Tama Zoological Park (Tokyo, Japan), Izu Animal Kingdom (Shizuoka, Japan), and Shizuoka Municipal Nihondaira zoo (Shizuoka, Japan) for kindly providing feline samples.
The authors declare that no competing interests exist.
Consent for publication
Ethics approval and consent to participate
This study was supported in part by AMED J-PRIDE 17fm0208006h0001 (to KS); JST CREST (to KS); JSPS KAKENHI Grants-in-Aid for Scientific Research C 15K07166 (to KS), Scientific Research B (Generative Research Fields) 16KT0111 (to KS), and Scientific Research on Innovative Areas 16H06429 (to S. Nakagawa and KS), 16K21723 (to S.Nakagawa and KS), 17H05823 (to S. Nakagawa) and 17H05813 (to KS); Takeda Science Foundation (to KS); Salt Science Research Foundation (to KS); Smoking Research Foundation (to KS); Chube Ito Foundation (to KS); Fordays Self-Reliance Support in Japan (to KS); Mishima Kaiun Memorial Foundation (to KS); Tobemaki Foundation (to KS); ONO Medical Research Foundation (to KS); MEXT-Supported Program for the Strategic Research Foundation at Private Universities (to S.Nakagawa and MTU); Joint Usage/Research Center program of Institute for Frontier Life and Medical Sciences Kyoto University (to S.Nakagawa); JSPS Core-to-Core program, A. Advanced Research Networks (to Y.Koyanagi); and AMED Research on HIV/AIDS 16fk0410203h002 (to Y.Koyanagi).
Springer Nature remains neutral with regard to jurisdictional claims in maps published and institutional affiliations.