Fengrong Ren

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 models fit to data and its utility in phylogeny reconstruction and caution against use of overly complex substitution models.

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.

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.

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

ABSTRACT 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.

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.

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.

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

ABSTRACT 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.

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.

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.

A new algorithm for analysis of within-host HIV-1 evolution.

A new algorithm for inferring the evolution of within-host viral sequences is presented. A sequential-linking approach is developed so that a longitudinal phylogenetic tree can be reconstructed from sequential molecular data that are obtained at different time points from the same host. The algorithm employs a codon-based model, which uses a Markov process to describe substitutions between codons, to calculate nonsynonymous and synonymous substitution rates and to distinguish positive selection and neutral evolution. The algorithm is applied to a data set of the V3 region of the HIV-1 envelope genes sequenced at different years after the infection of a single patient. The results suggest that this algorithm may provide a more realistic description of viral evolution than traditional evolutionary models, because it accounts for both neutral and adaptive evolution, and reconstructs a longitudinal phylogenetic tree that describes the dynamic process of viral evolution.