Feline Immunodeficiency Virus Evolutionarily Acquires Two Proteins, Vif and Protease, Capable of Antagonizing Feline APOBEC3.

The interplay between viral and host proteins has been well studied to elucidate virus-host interactions and their relevance to virulence. Mammalian genes encode apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3) proteins, which act as intrinsic restriction factors against lentiviruses. To overcome APOBEC3-mediated antiviral actions, lentiviruses have evolutionarily acquired an accessory protein, viral infectivity factor (Vif), and Vif degrades host APOBEC3 proteins via a ubiquitin/proteasome-dependent pathway. Although the Vif-APOBEC3 interaction and its evolutionary significance, particularly those of primate lentiviruses (including HIV) and primates (including humans), have been well investigated, those of nonprimate lentiviruses and nonprimates are poorly understood. Moreover, the factors that determine lentiviral pathogenicity remain unclear. Here, we focus on feline immunodeficiency virus (FIV), a pathogenic lentivirus in domestic cats, and the interaction between FIV Vif and feline APOBEC3 in terms of viral virulence and evolution. We reveal the significantly reduced diversity of FIV subtype B compared to that of other subtypes, which may associate with the low pathogenicity of this subtype. We also demonstrate that FIV subtype B Vif is less active with regard to feline APOBEC3 degradation. More intriguingly, we further reveal that FIV protease cleaves feline APOBEC3 in released virions. Taken together, our findings provide evidence that a lentivirus encodes two types of anti-APOBEC3 factors, Vif and viral protease.IMPORTANCE During the history of mammalian evolution, mammals coevolved with retroviruses, including lentiviruses. All pathogenic lentiviruses, excluding equine infectious anemia virus, have acquired the vif gene via evolution to combat APOBEC3 proteins, which are intrinsic restriction factors against exogenous lentiviruses. Here we demonstrate that FIV, a pathogenic lentivirus in domestic cats, antagonizes feline APOBEC3 proteins by both Vif and a viral protease. Furthermore, the Vif proteins of an FIV subtype (subtype B) have attenuated their anti-APOBEC3 activity through evolution. Our findings can be a clue to elucidate the complicated evolutionary processes by which lentiviruses adapt to mammals.

A naturally occurring bovine APOBEC3 confers resistance to bovine lentiviruses: implication for the co-evolution of bovids and their lentiviruses.

Mammals have co-evolved with lentiviruses for a long time. As evidence, viral infectivity factor (Vif), encoded by lentiviruses, antagonizes the anti-viral action of cellular APOBEC3 of their hosts. Here, we address the co-evolutionary dynamics of bovine APOBEC3 and the following two bovine lentiviruses: bovine immunodeficiency virus (BIV) and Jembrana disease virus (JDV). We determined the sequences of three APOBEC3 genes of bovids belonging to the genera Bos and Bison and showed that bovine APOBEC3Z3 is under a strong positive selection. We found that APOBEC3Z3 of gaur, a bovid in the genus Bos, acquired resistance to JDV Vif-mediated degradation after diverging from the other bovids through conversion of the structural composition of the loop 1 domain. Interestingly, the resistance of gaur APOBEC3Z3 can be attributed to the positive selection of residue 62. This study provides the first evidence, suggesting that a co-evolutionary arms race between bovids and lentiviruses occurred in Asia.

Small ruminant lentiviral Vif proteins commonly utilize cyclophilin A, an evolutionarily and structurally conserved protein, to degrade ovine and caprine APOBEC3 proteins.

Mammals have co-evolved with retroviruses, including lentiviruses, over a long period. Evidence supporting this contention is that viral infectivity factor (Vif) encoded by lentiviruses antagonizes the anti-viral action of cellular apolipoprotein B mRNA editing enzyme catalytic polypeptide-like 3 (APOBEC3) of the host. To orchestrate E3 ubiquitin ligase complex for APOBEC3 degradation, Vifs utilize mammalian proteins such as core-binding factor beta (CBFB; for primate lentiviruses) or cyclophilin A (CYPA; for Maedi-Visna virus [MVV]). However, the co-evolutionary relationship between lentiviral Vif and the mammalian proteins associated with Vif-mediated APOBEC3 degradation is poorly understood. Moreover, it is unclear whether Vif proteins of small ruminant lentiviruses (SRLVs), including MVV and caprine arthritis encephalitis virus (CAEV), commonly utilize CYPA to degrade the APOBEC3 of their hosts. In this study, molecular phylogenetic and protein homology modeling revealed that Vif co-factors are evolutionarily and structurally conserved. It was also found that not only MVV but also CAEV Vifs degrade APOBEC3 of both sheep and goats and that CAEV Vifs interact with CYPA. These findings suggest that lentiviral Vifs chose evolutionarily and structurally stable proteins as their partners (e.g., CBFB or CYPA) for APOBEC3 degradation and, particularly, that SRLV Vifs evolved to utilize CYPA as their co-factor in degradation of ovine and caprine APOBEC3.

A Naturally Occurring Domestic Cat APOBEC3 Variant Confers Resistance to Feline Immunodeficiency Virus Infection.

UNLABELLED: Apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3; A3) DNA cytosine deaminases can be incorporated into progeny virions and inhibit lentiviral replication. On the other hand, viral infectivity factor (Vif) of lentiviruses antagonizes A3-mediated antiviral activities by degrading A3 proteins. It is known that domestic cat (Felis catus) APOBEC3Z3 (A3Z3), the ortholog of human APOBEC3H, potently suppresses the infectivity of vif-defective feline immunodeficiency virus (FIV). Although a recent report has shown that domestic cat encodes 7 haplotypes (hap I to hap VII) of A3Z3, the relevance of A3Z3 polymorphism in domestic cats with FIV Vif has not yet been addressed. In this study, we demonstrated that these feline A3Z3 variants suppress vif-defective FIV infectivity. We also revealed that codon 65 of feline A3Z3 is a positively selected site and that A3Z3 hap V is subject to positive selection during evolution. It is particularly noteworthy that feline A3Z3 hap V is resistant to FIV Vif-mediated degradation and still inhibits vif-proficient viral infection. Moreover, the side chain size, but not the hydrophobicity, of the amino acid at position 65 determines the resistance to FIV Vif-mediated degradation. Furthermore, phylogenetic analyses have led to the inference that feline A3Z3 hap V emerged approximately 60,000 years ago. Taken together, these findings suggest that feline A3Z3 hap V may have been selected for escape from an ancestral FIV. This is the first evidence for an evolutionary "arms race" between the domestic cat and its cognate lentivirus. IMPORTANCE: Gene diversity and selective pressure are intriguing topics in the field of evolutionary biology. A direct interaction between a cellular protein and a viral protein can precipitate an evolutionary arms race between host and virus. One example is primate APOBEC3G, which potently restricts the replication of primate lentiviruses (e.g., human immunodeficiency virus type 1 [HIV-1] and simian immunodeficiency virus [SIV]) if its activity is not counteracted by the viral Vif protein. Here we investigate the ability of 7 naturally occurring variants of feline APOBEC3, APOBEC3Z3 (A3Z3), to inhibit FIV replication. Interestingly, one feline A3Z3 variant is dominant, restrictive, and naturally resistant to FIV Vif-mediated degradation. Phylogenetic analyses revealed that the ancestral change that generated this variant could have been caused by positive Darwinian selection, presumably due to an ancestral FIV infection. The experimental-phylogenetic investigation sheds light on the evolutionary history of the domestic cat, which was likely influenced by lentiviral infection.

Coevolutionary dynamics between tribe Cercopithecini tetherins and their lentiviruses.

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.

Pandemic HIV-1 Vpu overcomes intrinsic herd immunity mediated by tetherin.

Among the four groups of HIV-1 (M, N, O, and P), HIV-1M alone is pandemic and has rapidly expanded across the world. However, why HIV-1M has caused a devastating pandemic while the other groups remain contained is unclear. Interestingly, only HIV-1M Vpu, a viral protein, can robustly counteract human tetherin, which tethers budding virions. Therefore, we hypothesize that this property of HIV-1M Vpu facilitates human-to-human viral transmission. Adopting a multilayered experimental-mathematical approach, we demonstrate that HIV-1M Vpu confers a 2.38-fold increase in the prevalence of HIV-1 transmission. When Vpu activity is lost, protected human populations emerge (i.e., intrinsic herd immunity develops) through the anti-viral effect of tetherin. We also reveal that all Vpus of transmitted/founder HIV-1M viruses maintain anti-tetherin activity. These findings indicate that tetherin plays the role of a host restriction factor, providing 'intrinsic herd immunity', whereas Vpu has evolved in HIV-1M as a tetherin antagonist.

Vif determines the requirement for CBF-ß in APOBEC3 degradation.

APOBEC3 (apolipoprotein B mRNA editing enzyme catalytic polypeptide-like 3) proteins are cellular DNA deaminases that restrict a broad spectrum of lentiviruses. This process is counteracted by Vif (viral infectivity factor) of lentiviruses, which binds APOBEC3s and promotes their degradation. CBF-ß (core binding factor subunit ß) is an essential co-factor for the function of human immunodeficiency virus type 1 Vif to degrade human APOBEC3s. However, the requirement for CBF-ß in Vif-mediated degradation of other mammalian APOBEC3 proteins is less clear. Here, we determined the sequence of feline CBFB and performed phylogenetic analyses. These analyses revealed that mammalian CBFB is under purifying selection. Moreover, we demonstrated that CBF-ß is dispensable for feline immunodeficiency virus Vif-mediated degradation of APOBEC3s of its host. These findings suggested that primate lentiviruses have adapted to use CBF-ß, an evolutionary stable protein, to counteract APOBEC3 proteins of their hosts after diverging from other lentiviruses.

Characterization of red-capped mangabey tetherin: implication for the co-evolution of primates and their lentiviruses.

Primate lentiviruses including human immunodeficiency virus type 1 (HIV-1) and simian immunodeficiency viruses (SIVs) evolved through the acquisition of antagonists against intrinsic host restriction factors, such as tetherin. It is widely accepted that HIV-1 has emerged by zoonotic transmission of SIV in chimpanzee (SIVcpz), and that SIVcpz Nef protein antagonizes chimpanzee tetherin. Although Nef of SIVcpz shares a common ancestor with that of SIVrcm, an SIV in red-capped mangabey (Cercocebus torquatus), it remains unclear whether SIVrcm Nef can antagonize tetherin of its natural host. In this study, we determine the sequence of red-capped mangabey tetherin for the first time and directly demonstrate that SIVrcm Nef is the bona fide antagonist of red-capped mangabey tetherin. These findings suggest that SIVrcm Nef is the functional ancestor of SIVcpz Nef. Moreover, molecular phylogenetic analyses reveal that tetherins of the genus Cercocebus have experienced adaptive evolution, which is presumably promoted by primate lentiviruses.

Change of positive selection pressure on HIV-1 envelope gene inferred by early and recent samples.

HIV-1 infection has been on the rise in Japan recently, and the main transmission route has changed from blood transmission in the 1980s to homo- and/or hetero-sexual transmission in the 2000s. The lack of early viral samples with clinical information made it difficult to investigate the possible virological changes over time. In this study, we sequenced 142 full-length env genes collected from 16 Japanese subjects infected with HIV-1 in the 1980s and in the 2000s. We examined the diversity change in sequences and potential adaptive evolution of the virus to the host population. We used a codon-based likelihood method under the branch-site and clade models to detect positive selection operating on the virus. The clade model was extended to account for different positive selection pressures in different viral populations. The result showed that the selection pressure was weaker in the 2000s than in the 1980s, indicating that it might have become easier for the HIV to infect a new host and to develop into AIDS now than 20 years ago and that the HIV may be becoming more virulent in the Japanese population. The study provides useful information on the surveillance of HIV infection and highlights the utility of the extended clade models in analysis of virus populations which may be under different selection pressures.

Within-host co-evolution of Gag P453L and protease D30N/N88D demonstrates virological advantage in a highly protease inhibitor-exposed HIV-1 case.

To better understand the mechanism of HIV group-specific antigen (Gag) and protease (PR) co-evolution in drug-resistance acquisition, we analyzed a drug-resistance case by both bioinformatics and virological methods. We especially considered the quality of sequence data and analytical accuracy by introducing single-genome sequencing (SGS) and Spidermonkey/Bayesian graphical models (BGM) analysis, respectively. We analyzed 129 HIV-1 Gag-PR linkage sequences obtained from 8 time points, and the resulting sequences were applied to the Spidermonkey co-evolution analysis program, which identified ten mutation pairs as significantly co-evolving. Among these, we focused on associations between Gag-P453L, the P5' position of the p1/p6 cleavage-site mutation, and PR-D30N/N88D nelfinavir-resistant mutations, and attempted to clarify their virological significance in vitro by constructing recombinant clones. The results showed that P453L(Gag) has the potential to improve replication capacity and the Gag processing efficiency of viruses with D30N(PR)/N88D(PR) but has little effect on nelfinavir susceptibility. Homology modeling analysis suggested that hydrogen bonds between the 30th PR residue and the R452Gag are disturbed by the D30N(PR) mutation, but the impaired interaction is compensated by P453L(Gag) generating new hydrophobic interactions. Furthermore, database analysis indicated that the P453L(Gag)/D30N(PR)/N88D(PR) association was not specific only to our clinical case, but was common among AIDS patients.

Inferring within-patient HIV-1 evolutionary dynamics under anti-HIV therapy using serial virus samples with vSPA.

BACKGROUND: Analysis of within-patient HIV evolution under anti-HIV therapy is crucial to a better understanding the possible mechanisms of HIV drug-resistance acquisition. The high evolutionary rate of HIV allows us to trace its evolutionary process in real time by analyzing virus samples serially collected from the same patient. However, such studies are still uncommon due to the lack of powerful computational methods designed for serial virus samples. In this study, we develop a computational method, vSPA (viral Sequential Pathway Analysis), which groups viral sequences from the same sampling time into clusters and traces the evolution between clusters over sampling times. The method makes use of information of different sampling times and traces the evolution of important amino acid mutations. Second, a permutation test at the codon level is conducted to determine the threshold of the correlation coefficient for clustering viral quasispecies. We applied vSPA to four large data sets of HIV-1 protease and reverse transcriptase genes serially collected from two AIDS patients undergoing anti-HIV therapy over several years. RESULTS: The results show that vSPA can trace within-patient HIV evolution by detecting many amino acid changes, including important drug-resistant mutations, and by classifying different viral quasispecies coexisting during different periods of the therapy. CONCLUSION: Given that many new anti-HIV drugs will be available in the near future, vSPA may be useful for quickly providing information on the acquisition of HIV drug-resistant mutations by monitoring the within-patient HIV evolution under anti-HIV therapy as a computational approach.

A likelihood look at the supermatrix-supertree controversy.

Supermatrix and supertree methods are two strategies advocated for phylogenetic analysis of sequence data from multiple gene loci, especially when some species are missing at some loci. The supermatrix method concatenates sequences from multiple genes into a data supermatrix for phylogenetic analysis, and ignores differences in evolutionary dynamics among the genes. The supertree method analyzes each gene separately and assembles the subtrees estimated from individual genes into a supertree for all species. Most algorithms suggested for supertree construction lack statistical justifications and ignore uncertainties in the subtrees. Instead of supermatrix or supertree, we advocate the use of likelihood function to combine data from multiple genes while accommodating their differences in the evolutionary process. This combines the strengths of the supermatrix and supertree methods while avoiding their drawbacks. We conduct computer simulation to evaluate the performance of the supermatrix, supertree, and maximum likelihood methods applied to two phylogenetic problems: molecular-clock dating of species divergences and reconstruction of species phylogenies. The results confirm the theoretical superiority of the likelihood method. Supertree or separate analyses of data of multiple genes may be useful in revealing the characteristics of the evolutionary process of multiple gene loci, and the information may be used to formulate realistic models for combined analysis of all genes by likelihood.

A unique amino acid substitution, T126I, in human genotype C of hepatitis B virus S gene and its possible influence on antigenic structural change.

Amino acid substitutions in the S gene of hepatitis B virus (HBV), especially in the 'a' determinant region, have been suggested to affect the antigenicity of the virus and the clinical outcome of the infected patient. However, no convincing evidence has been presented for this hypothesis, partly because the 3D structure of the S protein has not been determined. Comparative analysis of viral genes offers an approach to testing this hypothesis, as it may reveal signals of natural selection and provide insights into the functional significance of the observed amino acid substitutions. In this study, we analyze HBV S gene sequences obtained from 24 patients infected with HBV genotypes B or C, together with 16 representative viral strains of HBV genotypes A-F retrieved from GenBank. We use phylogenetic methods to infer evolutionary changes among HBV genotypes and to identify amino acid residues potentially under positive selective pressure. Furthermore, we employ the fragment assembly method to predict structural changes in the S protein. The results showed that an amino acid substitution within the 'a' determinant, T126I, was unique to genotype C, may affect the antigenicity of the HBsAg, and may result in poorer clinical outcomes of patients infected with genotype C viral strains. We suggest that an integrated approach of evolutionary comparison and structural prediction is useful in generating hypotheses for further laboratory validation.

An empirical examination of the utility of codon-substitution models in phylogeny reconstruction.

Models of codon substitution have been commonly used to compare protein-coding DNA sequences and are particularly effective in detecting signals of natural selection acting on the protein. Their utility in reconstructing molecular phylogenies and in dating species divergences has not been explored. Codon models naturally accommodate synonymous and nonsynonymous substitutions, which occur at very different rates and may be informative for recent and ancient divergences, respectively. Thus codon models may be expected to make an efficient use of phylogenetic information in protein-coding DNA sequences. Here we applied codon models to 106 protein-coding genes from eight yeast species to reconstruct phylogenies using the maximum likelihood method, in comparison with nucleotide- and amino acid-based analyses. The results appeared to confirm that expectation. Nucleotide-based analysis, under simplistic substitution models, were efficient in recovering recent divergences whereas amino acid-based analysis performed better at recovering deep divergences. Codon models appeared to combine the advantages of amino acid and nucleotide data and had good performance at recovering both recent and deep divergences. Estimation of relative species divergence times using amino acid and codon models suggested that translation of gene sequences into proteins led to information loss of from 30% for deep nodes to 66% for recent nodes. Although computational burden makes codon models unfeasible for tree search in large data sets, we suggest that they may be useful for comparing candidate trees. Nucleotide models that accommodate the differences in evolutionary dynamics at the three codon positions also performed well, at much less computational cost. We discuss the relationship between a model's fit to data and its utility in phylogeny reconstruction and caution against use of overly complex substitution models.

Longitudinal phylogenetic tree of within-host viral evolution from noncontemporaneous samples: a distance-based sequential-linking method.

A new method for reconstructing phylogenetic relationships of within-host (patient) viral evolution from noncontemporaneous samples is presented. This method has two important features: noncontemporaneous viral samples can be dealt with by a simple computing algorithm, and both neutral and adaptive evolution patterns occurring during the process of viral evolution can be estimated. In our previous study, we proposed a preliminary formulation of this algorithm that was based on the maximum likelihood method. However, that preliminary formulation was difficult to use because the calculation of the likelihood required an extremely large amount of time and the number of possible tree topologies increased exponentially according to the increase in the number of viral variants. In this paper, we propose another new algorithm, referred to as a distance-based sequential-linking algorithm, in which the neighbor-joining method is employed for reconstruction of the longitudinal phylogenetic tree from serial viral samples. This algorithm is applied to a longitudinal data set of the env gene (V3 region) of human immunodeficiency virus type 1 (HIV-1) obtained over 7 years after the infection of a single patient. The results suggest that this method can successfully reconstruct a longitudinal phylogenetic tree from noncontemporaneous viral samples within a reasonable calculation time. This revised method proved to be a useful tool for estimating the dynamic process of within-host viral evolution.