Synonymous substitution rates predict HIV disease progression as a result of underlying replication dynamics

Dimensionality reduction distills complex evolutionary relationships in seasonal influenza and SARS-CoV-2

Public health researchers and practitioners commonly infer phylogenies from viral genome sequences to understand transmission dynamics and identify clusters of genetically-related samples. However, viruses that reassort or recombine violate phylogenetic assumptions and require more sophisticated methods. Even when phylogenies are appropriate, they can be unnecessary or difficult to interpret without specialty knowledge. For example, pairwise distances between sequences can be enough to identify clusters of related samples or assign new samples to existing phylogenetic clusters. In this work, we tested whether dimensionality reduction methods could capture known genetic groups within two human pathogenic viruses that cause substantial human morbidity and mortality and frequently reassort or recombine, respectively: seasonal influenza A/H3N2 and SARS-CoV-2. We applied principal component analysis (PCA), multidimensional scaling (MDS), t-distributed stochastic neighbor embedding (t-SNE), and uniform manifold approximation and projection (UMAP) to sequences with well-defined phylogenetic clades and either reassortment (H3N2) or recombination (SARS-CoV-2). For each low-dimensional embedding of sequences, we calculated the correlation between pairwise genetic and Euclidean distances in the embedding and applied a hierarchical clustering method to identify clusters in the embedding. We measured the accuracy of clusters compared to previously defined phylogenetic clades, reassortment clusters, or recombinant lineages. We found that MDS embeddings accurately represented pairwise genetic distances including the intermediate placement of recombinant SARS-CoV-2 lineages between parental lineages. Clusters from t-SNE embeddings accurately recapitulated known phylogenetic clades, H3N2 reassortment groups, and SARS-CoV-2 recombinant lineages. We show that simple statistical methods without a biological model can accurately represent known genetic relationships for relevant human pathogenic viruses. Our open source implementation of these methods for analysis of viral genome sequences can be easily applied when phylogenetic methods are either unnecessary or inappropriate.

Higher HIV-1 evolutionary rate is associated with cytotoxic T lymphocyte escape mutations in infants

Escape from cytotoxic T lymphocyte (CTL) responses toward HIV-1 Gag and Nef has been associated with reduced control of HIV-1 replication in adults. However, less is known about CTL-driven immune selection in infants as longitudinal studies of infants are limited. Here, 1,210 gag and 1,264 nef sequences longitudinally collected within 15 months after birth from 14 HIV-1 perinatally infected infants and their mothers were analyzed. The number of transmitted founder (T/F) viruses and associations between virus evolution, selection, CTL escape, and disease progression were determined. The analyses indicated that a paraphyletic-monophyletic relationship between the mother-infant sequences was common (80%), and that the HIV-1 infection was established by a single T/F virus in 10 of the 12 analyzed infants (83%). Furthermore, most HIV-1 CTL escape mutations among infants were transmitted from the mothers and did not revert during the first year of infection. Still, immune-driven selection was observed at approximately 3 months after HIV-1 infection in infants. Moreover, virus populations with CTL escape mutations in gag evolved faster than those without, independently of disease progression rate. These findings expand the current knowledge of HIV-1 transmission, evolution, and CTL escape in infant HIV-1 infection and are relevant for the development of immune-directed interventions in infants.IMPORTANCEDespite increased coverage in antiretroviral therapy for the prevention of perinatal transmission, paediatric HIV-1 infection remains a significant public health concern, especially in areas of high HIV-1 prevalence. Understanding HIV-1 transmission and the subsequent virus adaptation from the mother to the infant's host environment, as well as the viral factors that affect disease outcome, is important for the development of early immune-directed interventions for infants. This study advances our understanding of vertical HIV-1 transmission, and how infant immune selection pressure is shaping the intra-host evolutionary dynamics of HIV-1.

The genomic evolutionary dynamics and global circulation patterns of respiratory syncytial virus

Respiratory syncytial virus (RSV) is a leading cause of acute lower respiratory tract infection in young children and the second leading cause of infant death worldwide. While global circulation has been extensively studied for respiratory viruses such as seasonal influenza, and more recently also in great detail for SARS-CoV-2, a lack of global multi-annual sampling of complete RSV genomes limits our understanding of RSV molecular epidemiology. Here, we capitalise on the genomic surveillance by the INFORM-RSV study and apply phylodynamic approaches to uncover how selection and neutral epidemiological processes shape RSV diversity. Using complete viral genome sequences, we show similar patterns of site-specific diversifying selection among RSVA and RSVB and recover the imprint of non-neutral epidemic processes on their genealogies. Using a phylogeographic approach, we provide evidence for air travel governing the global patterns of RSVA and RSVB spread, which results in a considerable degree of phylogenetic mixing across countries. Our findings highlight the potential of systematic global RSV genomic surveillance for transforming our understanding of global RSV spread.

PCR-Based Detection and Genetic Characterization of Parainfluenza Virus 5 Detected in Pigs in Korea from 2016 to 2018

This study applied a molecular-based method to detect parainfluenza virus 5 (PIV5) collected from 2016 to 2018 in nine provinces of Republic of Korea. We demonstrated that PIV5 was detectable in both serum and pooled organs at an average positive rate of 1.78% (99/5566). Among these, the complete genome sequence of 15,246 nucleotides was obtained for 12 field strains. Three out of the 12 strains had the lowest genetic identity (96.20-96.68%) among the 21 porcine PIV5 genomes collected in Germany, China, India, and Republic of Korea from 1998 to 2017. By analyzing a large collection of complete genome sequences of the structural protein-coding F and HN genes, this study proposed a classification of PIV5 into two lineages, 1 and 2, and identified that group 2.2.2 within sub-lineage 2.2 was substantially divergent. The evolution of two structural protein-coding genes was largely under purifying selection. A few codons (6/9 for the F gene, 7/8 for the HN gene) had elevated dN/dS values, which were loaded on internal branches and were predicted to be related to beneficial trait(s) of the virus.

Molecular Basis of the Toxigenic Vibrio cholerae O1 Serotype Switch from Ogawa to Inaba in Haiti

Toxigenic Vibrio cholerae O1 serotype Ogawa was introduced involuntarily into Haiti in October 2010, and virtually all of the clinical strains isolated during the first 5 years of the epidemic were Ogawa. Inaba strains were identified intermittently prior to 2015, with diverse mutations resulting in a common phenotype. In 2015, the percentage of clinical infections due to the Inaba serotype began to rapidly increase, with Inaba supplanting Ogawa as the dominant serotype during the subsequent 4 years. We investigated the molecular basis of the serotype switch and confirmed that all Inaba strains had the same level of mRNA expression of the wbeT genes, as well as the same translation levels for the truncated WbeT proteins in the V. cholerae Inaba isolates. Neither wbeT gene expression levels, differential mutations, or truncation size of the WbeT proteins appeared to be responsible for the successful Inaba switch in 2015. Our phylodynamic analysis demonstrated that the V. cholerae Inaba strains in Haiti evolved directly from Ogawa strains and that a significant increase of diversifying selection at the population level occurred at the time of the Ogawa-Inaba switch. We conclude that the emergence of the Inaba serotype was driven by diversifying selection, independent of the mutational pattern in the wbeT gene. IMPORTANCE Our phylodynamic analysis demonstrated that Vibrio cholerae Inaba strains in Haiti evolved directly from Ogawa strains. Our results support the hypothesis that after an initial Ogawa-dominated epidemic wave, V. cholerae Inaba was able to become the dominant strain thanks to a selective advantage driven by ongoing diversifying selection, independently from the mutational pattern in the wbeT gene.

Deep-sequencing of viral genomes from a large and diverse cohort of treatment-naive HIV-infected persons shows associations between intrahost genetic diversity and viral load

BACKGROUND

Infection with human immunodeficiency virus type 1 (HIV) typically results from transmission of a small and genetically uniform viral population. Following transmission, the virus population becomes more diverse because of recombination and acquired mutations through genetic drift and selection. Viral intrahost genetic diversity remains a major obstacle to the cure of HIV; however, the association between intrahost diversity and disease progression markers has not been investigated in large and diverse cohorts for which the majority of the genome has been deep-sequenced. Viral load (VL) is a key progression marker and understanding of its relationship to viral intrahost genetic diversity could help design future strategies for HIV monitoring and treatment.

METHODS

We analysed deep-sequenced viral genomes from 2,650 treatment-naive HIV-infected persons to measure the intrahost genetic diversity of 2,447 genomic codon positions as calculated by Shannon entropy. We tested for associations between VL and amino acid (AA) entropy accounting for sex, age, race, duration of infection, and HIV population structure.

RESULTS

We confirmed that the intrahost genetic diversity is highest in the env gene. Furthermore, we showed that mean Shannon entropy is significantly associated with VL, especially in infections of >24 months duration. We identified 16 significant associations between VL (p-value<2.0x10-5) and Shannon entropy at AA positions which in our association analysis explained 13% of the variance in VL. Finally, equivalent analysis based on variation in HIV consensus sequences explained only 2% of VL variance.

CONCLUSIONS

Our results elucidate that viral intrahost genetic diversity is associated with VL and could be used as a better disease progression marker than HIV consensus sequence variants, especially in infections of longer duration. We emphasize that viral intrahost diversity should be considered when studying viral genomes and infection outcomes.

TRIAL REGISTRATION

Samples included in this study were derived from participants who consented in the clinical trial, START (NCT00867048) (23), run by the International Network for Strategic Initiatives in Global HIV Trials (INSIGHT). All the participant sites are listed here: http://www.insight-trials.org/start/my_phpscript/participating.php?by=site.

Emergence and Evolutionary Response of Vibrio cholerae to Novel Bacteriophage, Democratic Republic of the Congo1

Cholera causes substantial illness and death in Africa. We analyzed 24 toxigenic Vibrio cholerae O1 strains isolated in 2015-2017 from patients in the Great Lakes region of the Democratic Republic of the Congo. Strains originating in southern Asia appeared to be part of the T10 introduction event in eastern Africa. We identified 2 main strain lineages, most recently a lineage corresponding to sequence type 515, a V. cholerae cluster previously reported in the Lake Kivu region. In 41% of fecal samples from cholera patients, we also identified a novel ICP1 (Bangladesh cholera phage 1) bacteriophage, genetically distinct from ICP1 isolates previously detected in Asia. Bacteriophage resistance occurred in distinct clades along both internal and external branches of the cholera phylogeny. This bacteriophage appears to have served as a major driver for cholera evolution and spread, and its appearance highlights the complex evolutionary dynamic that occurs between predatory phage and bacterial host.

Intrahost evolution of the HIV-2 capsid correlates with progression to AIDS

HIV-2 infection will progress to AIDS in most patients without treatment, albeit at approximately half the rate of HIV-1 infection. HIV-2 capsid (p26) amino acid polymorphisms are associated with lower viral loads and enhanced processing of T cell epitopes, which may lead to protective Gag-specific T cell responses common in slower progressors. Lower virus evolutionary rates, and positive selection on conserved residues in HIV-2 env have been associated with slower progression to AIDS. In this study we analysed 369 heterochronous HIV-2 p26 sequences from 12 participants with a median age of 30 years at enrolment. CD4% change over time was used to stratify participants into relative faster and slower progressor groups. We analysed p26 sequence diversity evolution, measured site-specific selection pressures and evolutionary rates, and determined if these evolutionary parameters were associated with progression status. Faster progressors had lower CD4% and faster CD4% decline rates. Median pairwise sequence diversity was higher in faster progressors (5.7x10-3 versus 1.4x10-3 base substitutions per site, P<0.001). p26 evolved under negative selection in both groups (dN/dS=0.12). Median virus evolutionary rates were higher in faster than slower progressors - synonymous rates: 4.6x10-3 vs. 2.3x10-3; and nonsynonymous rates: 6.9x10-4 vs. 2.7x10-4 substitutions/site/year, respectively. Virus evolutionary rates correlated negatively with CD4% change rates (ρ = -0.8, P=0.02), but not CD4% level. The signature amino acid at p26 positions 6, 12 and 119 differed between faster (6A, 12I, 119A) and slower (6G, 12V, 119P) progressors. These amino acid positions clustered near to the TRIM5α/p26 hexamer interface surface. p26 evolutionary rates were associated with progression to AIDS and were mostly driven by synonymous substitutions. Nonsynonymous evolutionary rates were an order of magnitude lower than synonymous rates, with limited amino acid sequence evolution over time within hosts. These results indicate HIV-2 p26 may be an attractive therapeutic target.

Plasma-Derived HIV-1 Virions Contain Considerable Levels of Defective Genomes

Genetically-characterizing full-length HIV-1 RNA is critical for identifying genetically-intact genomes and for comparing these RNA genomes to proviral DNA. We have developed a method for sequencing plasma-derived RNA using long-range sequencing (PRLS assay; ∼8.3 kb from gag to the 3′ end or ∼5 kb from integrase to the 3′ end). We employed the gag-3′ PRLS assay to sequence HIV-1 RNA genomes from ART-naive participants during acute/early infection (n = 6) or chronic infection (n = 2). On average, only 65% of plasma-derived genomes were genetically-intact. Defects were found in all genomic regions but were concentrated in env and pol. We compared these genomes to near-full-length proviral sequences from paired peripheral blood mononuclear cell (PBMC) samples for the acute/early group and found that near-identical (>99.98% identical) sequences were identified only during acute infection. For three participants who initiated therapy during acute infection, we used the int-3′ PRLS assay to sequence plasma-derived genomes from an analytical treatment interruption and identified 100% identical genomes between pretherapy and rebound time points. The PRLS assay provides a new level of sensitivity for understanding the genetic composition of plasma-derived HIV-1 RNA from viremic individuals either pretherapy or after treatment interruption, which will be invaluable in assessing possible HIV-1 curative strategies.

IMPORTANCE We developed novel plasma-derived RNA using long-range sequencing assays (PRLS assay; 8.3 kb, gag-3′, and 5.0 kb, int-3′). Employing the gag-3′ PRLS assay, we found that 26% to 51% of plasma-derived genomes are genetically-defective, largely as a result of frameshift mutations and deletions. These genetic defects were concentrated in the env region compared to gag and pol, likely a reflection of viral immune escape in env during untreated HIV-1 infection. Employing the int-3′ PRLS assay, we found that analytical treatment interruption (ATI) plasma-derived sequences were identical and genetically-intact. Several sequences from the ATI plasma samples were identical to viral sequences from pretherapy plasma and PBMC samples, indicating that HIV-1 reservoirs established prior to therapy contribute to viral rebound during an ATI. Therefore, near-full-length sequencing of HIV-1 particles is required to gain an accurate picture of the genetic landscape of plasma HIV-1 virions in studies of HIV-1 replication and persistence.

HIV/AIDS: Current Updates on the Disease, Treatment and Prevention

CCR5-delta 32 homozygous stem cell transplantation for HIV-infected individuals is being treated as a milestone in the global AIDS epidemic. Since 2008, when the second Berlin patient was cured from HIV after undergoing transplantation from a donor with delta-32 mutation, scientists are aiming for a long-term cure for the wider population. In 2019, a London patient became the second person to be free of HIV and came off the antiretroviral drugs completely. CCR5 gene is now being treated as a viable target for HIV treatment. It can be used in the treatment of HIV either through administration of drugs that bind to CCR5 and stop the receptor from working or through gene therapy to alter the CCR5 gene using CRISPR/Cas9 and prevent protein production. This review article aims to identify the obstacles and the need to overcome them in order to bridge the gap between current research and future potential cures for HIV.

Therapy Implications of Hepatitis C Virus Genetic Diversity

Hepatitis C virus (HCV) is an important human pathogen with a high chronicity rate. An estimated 71 million people worldwide are living with chronic hepatitis C (CHC) infection, which carries the risk of progression to hepatic fibrosis, cirrhosis, and hepatocellular carcinoma (HCC). Similar to other RNA viruses, HCV has a high rate of genetic variability generated by its high mutation rate and the actions of evolutionary forces over time. There are two levels of HCV genetic variability: intra-host variability, characterized by the distribution of HCV mutant genomes present in an infected individual, and inter-host variability, represented by the globally circulating viruses that give rise to different HCV genotypes and subtypes. HCV genetic diversity has important implications for virus persistence, pathogenesis, immune responses, transmission, and the development of successful vaccines and antiviral strategies. Here we will discuss how HCV genetic heterogeneity impacts viral spread and therapeutic control.

Deep Profiling Identifies Selection of Nonsynonymous Amino Acid Substitutions in HIV-1 Envelope During Early Infection

Understanding the evolutionary dynamics of the viruses within an individual at or near the moment of transmission can provide critical inputs for the design of an effective vaccine for HIV infection. In this study, high-throughput sequencing technology was employed to analyze the evolutionary rate in viruses obtained at a single time point from drug-naive recently infected infants and adults in the chronic stage of disease. Gene-wise nonsynonymous (pN) and synonymous (pS) mutation rates were estimated and compared between the two groups. Significant differences were observed in the evolutionary rates between viruses in the early and late stages of infection. Higher rates of adaptive mutations in the HIV-1 envelope gene (env) were found in the chronic viruses as compared with those in the early stages of HIV infection. Conversely, percentage of nonsynonymous substitutions in env was found to be higher in recently transmitted viruses. In addition, a positive correlation was found between mutation and the evolutionary rate, and infectivity titer in recent infection. Despite the small sample size, the study identified useful information about viral evolution on transmission-associated bottlenecks. The effect of intraindividual HIV-1 evolution at the population level was highly contemporary, and the higher percentage of nonsynonymous substitutions seen in env during recent HIV-1 infection has suggested a pattern of convergent evolution leading to a positive selection for survival fitness and disease progression.

Purifying Selection in Human Immunodeficiency Virus-1 pol Gene in Perinatally Human Immunodeficiency Virus-1-Infected Children Harboring Discordant Immunological Response and Virological Nonresponse to Long-Term Antiretroviral Therapy

Background

Biological monitoring of antiretroviral treatment (ART) in human immunodeficiency virus (HIV)-infected pediatric population remains challenging. The aim of the present study was to assess the long-term HIV-1 genetic diversity in pol gene in HIV-1-infected children in virological failure under antiretroviral regimen adapted according to the successive World Health Organization (WHO) guidelines for resource-constrained settings.

Methods

HIV-1 diversity in pol gene was assessed in HIV-1-infected children and adolescents born from HIV-infected mothers (median age at follow-up: 13.8 years) in virological failure (VF⁺) despite long-term regimen recommended by the WHO. The numbers of nonsynonymous substitutions per potential nonsynonymous site (dN) and of synonymous substitutions at potential synonymous sites (dS) in HIV-1 pol gene and the dN/dS ratios were used to estimate the selective pressure on circulating HIV-1.

Results

The immunological responses to ART basically corresponded to: 1) Full therapeutic failure with immunological (I^-) and virological nonresponses in one-quarter (24.6%) of study children ((I^-, VF⁺) subgroup); 2) Discordant immunovirological responses with paradoxical high CD4 T cell counts (I⁺) and high HIV-1 RNA load in the remaining cohort patients (75.4%) ((I⁺, VF⁺) subgroup). The mean dS was 1.8-fold higher in (I⁺, VF⁺) than (I^-, VF⁺) subgroup (25.9 ± 18.4 vs. 14.3 ± 10.8). In the (I⁺, VF⁺) subgroup, the mean dS was 1.6-fold higher than the mean dN. Finally, the mean dN/dS ratio was 2.1-fold lower in (I⁺, VF⁺) than (I^-, VF⁺) subgroup (0.6 ± 0.3 vs. 1.3 ± 0.7), indicating purifying selection in the immunovirological discordant (I⁺, VF⁺) subgroup and positive selection in the immunovirological failure (I^-, VF⁺) subgroup.

Conclusions

Children and adolescents in immunovirological therapeutic failure harbor positive selection of HIV-1 strains favoring diversifying in pol-encoded amino acids. In contrast, children with persistent discordant immunovirological responses show accumulation of mutations and purifying selection in pol gene sequences, indicating limited genetic evolution and likely suggesting genetic adaptation of viruses to host functional constraints.

Toxigenic Vibrio cholerae evolution and establishment of reservoirs in aquatic ecosystems

The spread of cholera in the midst of an epidemic is largely driven by direct transmission from person to person, although it is well-recognized that Vibrio cholerae is also capable of growth and long-term survival in aquatic ecosystems. While prior studies have shown that aquatic reservoirs are important in the persistence of the disease on the Indian subcontinent, an epidemiological view postulating that locally evolving environmental V. cholerae contributes to outbreaks outside Asia remains debated. The single-source introduction of toxigenic V. cholerae O1 in Haiti, one of the largest outbreaks occurring this century, with 812,586 suspected cases and 9,606 deaths reported through July 2018, provided a unique opportunity to evaluate the role of aquatic reservoirs and assess bacterial transmission dynamics across environmental boundaries. To this end, we investigated the phylogeography of both clinical and aquatic toxigenic V. cholerae O1 isolates and show robust evidence of the establishment of aquatic reservoirs as well as ongoing evolution of V. cholerae isolates from aquatic sites. Novel environmental lineages emerged from sequential population bottlenecks, carrying mutations potentially involved in adaptation to the aquatic ecosystem. Based on such empirical data, we developed a mixed-transmission dynamic model of V. cholerae, where aquatic reservoirs actively contribute to genetic diversification and epidemic emergence, which underscores the complexity of transmission pathways in epidemics and endemic settings and the need for long-term investments in cholera control at both human and environmental levels.

Source identification of HIV-1 transmission in three lawsuits Using Ultra-Deep pyrosequencing and phylogenetic analysis

BACKGROUND/PURPOSE

Intentional transmission of HIV-1 is a crime. Identifying the source of transmission between HIV-1 infected cases using phylogenetic analysis has limitations, including delayed examinations after the initiation of infection and ambiguity of phyletic relationships. This study was the first to introduce phylogenetic tree Results as forensic evidence in a trial in Taiwan.

METHODS

Three lawsuit cases from different district courts in Taiwan were chosen for this study. We identified the source of transmission between individuals in each lawsuit based on the maximum likelihood and Bayesian phylogenetic tree analyses using the HIV-1 sequences from molecular cloning and ultra-deep pyrosequencing (UDPS). Two gene regions of the HIV genome, env and gag, were involved.

RESULTS

The results of phylogenetic analysis using sequences from molecular cloning were clear and evidential enough in lawsuits 1 and 3. Due to the delayed sampling time, the result of sequences from molecular cloning in lawsuit 2 was ambiguous. Combined with the analyzed result of sequences from UDPS and epidemiological information, the source of transmission in lawsuit 2 was further identified.

CONCLUSION

Hence phylogenetic analyses cannot exclude the possibility of unsampled intermediaries, the data interpretation should be more careful and conservative, and it should not be considered as the only evidence for the source identification in a trial without epidemiological or serological information. The evaluation of the introduced UDPS method in the identification of transmission source has shown that the validity and evidential effects were still limited and need further optimization.

A non-parametric analytic framework for within-host viral phylogenies and a test for HIV-1 founder multiplicity

Phylogenetics is a powerful tool for understanding the diversification dynamics of viral pathogens. Here we present an extension of the spectral density profile of the modified graph Laplacian, which facilitates the characterization of within-host molecular evolution of viruses and the direct comparison of diversification dynamics between hosts. This approach is non-parametric and therefore fast and model-free. We used simulations of within-host evolutionary scenarios to evaluate the efficiency of our approach and to demonstrate the significance of interpreting a viral phylogeny by its spectral density profile in terms of diversification dynamics. The key features that are captured by the profile are positive selection on the viral gene (or genome), temporal changes in substitution rates, mutational fitness, and time between sampling. Using sequences from individuals infected with HIV-1, we showed the utility of this approach for characterizing within-host diversification dynamics, for comparing dynamics between hosts, and for charting disease progression in infected individuals sampled over multiple years. We furthermore propose a heuristic test for assessing founder heterogeneity, which allows us to classify infections with single and multiple HIV-1 founder viruses. This non-parametric approach can be a valuable complement to existing parametric approaches.

Variation in Intra-individual Lentiviral Evolution Rates: a Systematic Review of Human, Nonhuman Primate, and Felid Species

Lentiviral replication mediated by reverse transcriptase is considered to be highly error prone, leading to a high intra-individual evolution rate that promotes evasion of neutralization and persistent infection. Understanding lentiviral intra-individual evolutionary dynamics on a comparative basis can therefore inform research strategies to aid in studies of pathogenesis, vaccine design, and therapeutic intervention. We conducted a systematic review of intra-individual evolution rates for three species groups of lentiviruses-feline immunodeficiency virus (FIV), simian immunodeficiency virus (SIV), and human immunodeficiency virus (HIV). Overall, intra-individual rate estimates differed by virus but not by host, gene, or viral strain. Lentiviral infections in spillover (nonadapted) hosts approximated infections in primary (adapted) hosts. Our review consistently documents that FIV evolution rates within individuals are significantly lower than the rates recorded for HIV and SIV. FIV intra-individual evolution rates were noted to be equivalent to FIV interindividual rates. These findings document inherent differences in the evolution of FIV relative to that of primate lentiviruses, which may signal intrinsic difference of reverse transcriptase between these viral species or different host-viral interactions. Analysis of lentiviral evolutionary selection pressures at the individual versus population level is valuable for understanding transmission dynamics and the emergence of virulent and avirulent strains and provides novel insight for approaches to interrupt lentiviral infections.IMPORTANCE To the best of our knowledge, this is the first study that compares intra-individual evolution rates for FIV, SIV, and HIV following systematic review of the literature. Our findings have important implications for informing research strategies in the field of intra-individual virus dynamics for lentiviruses. We observed that FIV evolves more slowly than HIV and SIV at the intra-individual level and found that mutation rates may differ by gene sequence length but not by host, gene, strain, an experimental setting relative to a natural setting, or spillover host infection relative to primary host infection.

Divergence dating using mixed effects clock modelling: An application to HIV-1

The need to estimate divergence times in evolutionary histories in the presence of various sources of substitution rate variation has stimulated a rich development of relaxed molecular clock models. Viral evolutionary studies frequently adopt an uncorrelated clock model as a generic relaxed molecular clock process, but this may impose considerable estimation bias if discrete rate variation exists among clades or lineages. For HIV-1 group M, rate variation among subtypes has been shown to result in inconsistencies in time to the most recent common ancestor estimation. Although this calls into question the adequacy of available molecular dating methods, no solution to this problem has been offered so far. Here, we investigate the use of mixed effects molecular clock models, which combine both fixed and random effects in the evolutionary rate, to estimate divergence times. Using simulation, we demonstrate that this model outperforms existing molecular clock models in a Bayesian framework for estimating time-measured phylogenies in the presence of mixed sources of rate variation, while also maintaining good performance in simpler scenarios. By analysing a comprehensive HIV-1 group M complete genome data set we confirm considerable rate variation among subtypes that is not adequately modelled by uncorrelated relaxed clock models. The mixed effects clock model can accommodate this rate variation and results in a time to the most recent common ancestor of HIV-1 group M of 1920 (1915-25), which is only slightly earlier than the uncorrelated relaxed clock estimate for the same data set. The use of complete genome data appears to have a more profound impact than the molecular clock model because it reduces the credible intervals by 50 per cent relative to similar estimates based on short envelope gene sequences.

Molecular Signatures of HIV-1 Envelope Associated with HIV-Associated Neurocognitive Disorders

PURPOSE OF REVIEW

The HIV-1 envelope gene (env) has been an intense focus of investigation in the search for genetic determinants of viral entry and persistence in the central nervous system (CNS).

RECENT FINDINGS

Molecular signatures of CNS-derived HIV-1 env reflect the immune characteristics and cellular constraints of the CNS compartment. Although more readily found in those with advanced HIV-1 and HIV-associated neurocognitive disorders (HAND), molecular signatures distinguishing CNS-derived quasispecies can be identified early in HIV-1 infection, in the presence or absence of combination antiretroviral therapy (cART), and are dynamic. Amino acid signatures of CNS-compartmentalization and HAND have been identified across populations. While some significant overlap exists, none are universal. Detailed analyses of CNS-derived HIV-1 env have allowed researchers to identify a number of molecular determinants associated with neuroadaptation. Future investigations using comprehensive cohorts and longitudinal databases have the greatest potential for the identification of robust, validated signatures of HAND in the cART era.

Delayed disease progression in HIV-2: the importance of TRIM5α and the retroviral capsid

HIV-2 is thought to have entered the human population in the 1930s through cross-species transmission of SIV from sooty mangabeys in West Africa. Unlike HIV-1, HIV-2 has not led to a global pandemic, and recent data suggest that HIV-2 prevalence is declining in some West African states where it was formerly endemic. Although many early isolates of HIV-2 were derived from patients presenting with AIDS-defining illnesses, it was noted that a much larger proportion of HIV-2-infected subjects behaved as long-term non-progressors (LTNP) than their HIV-1-infected counterparts. Many HIV-2-infected adults are asymptomatic, maintaining an undetectable viral load for over a decade. However, despite lower viral loads, HIV-2 progresses to clinical AIDS without therapeutic intervention in most patients. In addition, successful treatment with anti-retroviral therapy (ART) is more challenging than for HIV-1. HIV-2 is significantly more sensitive to restriction by host restriction factor tripartite motif TRIM5α than HIV-1, and this difference in sensitivity is linked to differences in capsid structure. In this review we discuss the determinants of HIV-2 disease progression and focus on the important interactions between TRIM5α and HIV-2 capsid in long-term viral control.

Selection and Neutral Mutations Drive Pervasive Mutability Losses in Long-Lived Anti-HIV B-Cell Lineages

High-affinity antibodies arise within weeks of infection from the evolution of B-cell receptors under selection to improve antigen recognition. This rapid adaptation is enabled by the distribution of highly mutable "hotspot" motifs in B-cell receptor genes. High mutability in antigen-binding regions (complementarity determining regions [CDRs]) creates variation in binding affinity, whereas low mutability in structurally important regions (framework regions [FRs]) may reduce the frequency of destabilizing mutations. During the response, loss of mutational hotspots and changes in their distribution across CDRs and FRs are predicted to compromise the adaptability of B-cell receptors, yet the contributions of different mechanisms to gains and losses of hotspots remain unclear. We reconstructed changes in anti-HIV B-cell receptor sequences and show that mutability losses were ∼56% more frequent than gains in both CDRs and FRs, with the higher relative mutability of CDRs maintained throughout the response. At least 21% of the total mutability loss was caused by synonymous mutations. However, nonsynonymous substitutions caused most (79%) of the mutability loss in CDRs. Because CDRs also show strong positive selection, this result suggests that selection for mutations that increase binding affinity contributed to loss of mutability in antigen-binding regions. Although recurrent adaptation to evolving viruses could indirectly select for high mutation rates, we found no evidence of indirect selection to increase or retain hotspots. Our results suggest mutability losses are intrinsic to both the neutral and adaptive evolution of B-cell populations and might constrain their adaptation to rapidly evolving pathogens such as HIV and influenza.

Positive Selection at Key Residues in the HIV Envelope Distinguishes Broad and Strain-Specific Plasma Neutralizing Antibodies

Millions of people are still being infected with HIV decades after the first recognition of the virus. Currently, no vaccine is able to elicit bNAbs that will prevent infection by global HIV strains. Several studies have implicated HIV Env diversity in the development of breadth. However, Env evolution in individuals who fail to develop breadth despite mounting potent strain-specific neutralizing responses has not been well defined. Using longitudinal neutralization, epitope mapping, and sequence data from 14 participants, we found that overall measures of viral diversity were similar in all donors. However, the number of positively selected sites within Env epitopes was higher in bNAb participants than in strain-specific donors. We further identified common sites that were positively selected as bNAbs developed. These data indicate that while viral diversity is required for breadth, this should be highly targeted to specific residues to shape the elicitation of bNAbs by vaccination.

The development of HIV broadly neutralizing antibodies (bNAbs) has previously been shown to be associated with viral evolution and high levels of genetic diversity in the HIV envelope (Env) glycoprotein. However, few studies have examined Env evolution in those who fail to develop neutralization breadth in order to assess whether bNAbs result from distinct evolutionary pathways. We compared Env evolution in eight HIV-1-infected participants who developed bNAbs to six donors with similar viral loads who did not develop bNAbs over three years of infection. We focused on Env V1V2 and C3V4, as these are major targets for both strain-specific neutralizing antibodies (nAbs) and bNAbs. Overall evolutionary rates (ranging from 9.92 × 10⁻³ to 4.1 × 10⁻² substitutions/site/year) and viral diversity (from 1.1% to 6.5%) across Env, and within targeted epitopes, did not distinguish bNAb donors from non-bNAb donors. However, bNAb participants had more positively selected residues within epitopes than those without bNAbs, and several of these were common among bNAb donors. A comparison of the kinetics of strain-specific nAbs and bNAbs indicated that selection pressure at these residues increased with the onset of breadth. These data suggest that highly targeted viral evolution rather than overall envelope diversity is associated with neutralization breadth. The association of shared positively selected sites with the onset of breadth highlights the importance of diversity at specific positions in these epitopes for bNAb development, with implications for the development of sequential and cocktail immunization strategies.

IMPORTANCE Millions of people are still being infected with HIV decades after the first recognition of the virus. Currently, no vaccine is able to elicit bNAbs that will prevent infection by global HIV strains. Several studies have implicated HIV Env diversity in the development of breadth. However, Env evolution in individuals who fail to develop breadth despite mounting potent strain-specific neutralizing responses has not been well defined. Using longitudinal neutralization, epitope mapping, and sequence data from 14 participants, we found that overall measures of viral diversity were similar in all donors. However, the number of positively selected sites within Env epitopes was higher in bNAb participants than in strain-specific donors. We further identified common sites that were positively selected as bNAbs developed. These data indicate that while viral diversity is required for breadth, this should be highly targeted to specific residues to shape the elicitation of bNAbs by vaccination.

Low Postseroconversion CD4+ T-cell Level Is Associated with Faster Disease Progression and Higher Viral Evolutionary Rate in HIV-2 Infection

The relationship between HIV evolution and disease progression is fundamental to our understanding of HIV immune control and vaccine design. There are no clear definitions for faster and slower HIV-2 disease progression and for the relationship of the rate of progression with HIV-2 evolution. To address the hypothesis that viral evolution is correlated with disease progression in HIV-2 infection, we determined faster and slower disease progression based on follow-up data from a prospective cohort of police officers in Guinea-Bissau. The analysis showed that although the CD4⁺ T-cell level and the decline in the level were independently associated with progression to AIDS, only the CD4⁺ T-cell level or a combined CD4⁺ T-cell level/decline stratification was associated with the rate of HIV-2 evolution. The HIV-2 evolutionary rate was almost twice as high among the faster progressors as among the slower progressors. Importantly, this report defines previously unknown characteristics linking HIV-2 disease progression with virus evolution.

A positive correlation between virus evolutionary rate and disease progression has been shown for human immunodeficiency virus type 1 (HIV-1) infection. Much less is known about HIV-2, the second causative agent of AIDS. We analyzed 528 HIV-2 env V1-C3 sequences generated from longitudinal plasma samples that were collected from 16 study participants during a median observation time of 7.9 years (interquartile range [IQR], 5.2 to 14.0 years). Study participants were classified as faster or slower disease progressors based on longitudinal CD4⁺ T-cell data. The HIV-2 evolutionary rate was significantly associated with CD4⁺ T-cell levels and was almost twice as high among the faster progressors as among the slower progressors. Higher evolutionary rates were accounted for by both synonymous and nonsynonymous nucleotide substitutions. Moreover, slow disease progression was associated with stronger positive selection on HIV-2/SIVsm (simian immunodeficiency virus infecting sooty mangabey) surface-exposed conserved residues. This study demonstrated a number of previously unknown characteristics linking HIV-2 disease progression with virus evolution. Some of these findings distinguish HIV-2 from HIV-1 and may contribute to the understanding of differences in pathogenesis.

A majority of HIV persistence during antiretroviral therapy is due to infected cell proliferation

Antiretroviral therapy (ART) suppresses viral replication in people living with HIV. Yet, infected cells persist for decades on ART and viremia returns if ART is stopped. Persistence has been attributed to viral replication in an ART sanctuary and long-lived and/or proliferating latently infected cells. Using ecological methods and existing data, we infer that >99% of infected cells are members of clonal populations after one year of ART. We reconcile our results with observations from the first months of ART, demonstrating mathematically how a fossil record of historic HIV replication permits observed viral evolution even while most new infected cells arise from proliferation. Together, our results imply cellular proliferation generates a majority of infected cells during ART. Therefore, reducing proliferation could decrease the size of the HIV reservoir and help achieve a functional cure.

HIV infected cells persist for decades in patients under ART, but the mechanisms responsible remain unclear. Here, Reeves et al. use modeling approaches adapted from ecology to show that cellular proliferation, rather than viral replication, generates a majority of infected cells during ART.

Evolution of HIV-1 within untreated individuals and at the population scale in Uganda

HIV-1 undergoes multiple rounds of error-prone replication between transmission events, resulting in diverse viral populations within and among individuals. In addition, the virus experiences different selective pressures at multiple levels: during the course of infection, at transmission, and among individuals. Disentangling how these evolutionary forces shape the evolution of the virus at the population scale is important for understanding pathogenesis, how drug- and immune-escape variants are likely to spread in populations, and the development of preventive vaccines. To address this, we deep-sequenced two regions of the HIV-1 genome (p24 and gp41) from 34 longitudinally-sampled untreated individuals from Rakai District in Uganda, infected with subtypes A, D, and inter-subtype recombinants. This dataset substantially increases the availability of HIV-1 sequence data that spans multiple years of untreated infection, in particular for different geographical regions and viral subtypes. In line with previous studies, we estimated an approximately five-fold faster rate of evolution at the within-host compared to the population scale for both synonymous and nonsynonymous substitutions, and for all subtypes. We determined the extent to which this mismatch in evolutionary rates can be explained by the evolution of the virus towards population-level consensus, or the transmission of viruses similar to those that establish infection within individuals. Our findings indicate that both processes are likely to be important.

Neutral Theory and Rapidly Evolving Viral Pathogens

The evolution of viral pathogens is shaped by strong selective forces that are exerted during jumps to new hosts, confrontations with host immune responses and antiviral drugs, and numerous other processes. However, while undeniably strong and frequent, adaptive evolution is largely confined to small parts of information-packed viral genomes, and the majority of observed variation is effectively neutral. The predictions and implications of the neutral theory have proven immensely useful in this context, with applications spanning understanding within-host population structure, tracing the origins and spread of viral pathogens, predicting evolutionary dynamics, and modeling the emergence of drug resistance. We highlight the multiple ways in which the neutral theory has had an impact, which has been accelerated in the age of high-throughput, high-resolution genomics.

Subtle differences in selective pressures applied on the envelope gene of HIV-1 in pregnant versus non-pregnant women

Pregnancy is associated with modulations of maternal immunity that contribute to foeto-maternal tolerance. To understand whether and how these alterations impact antiviral immunity, a detailed cross-sectional analysis of selective pressures exerted on HIV-1 envelope amino-acid sequences was performed in a group of pregnant (n = 32) and non-pregnant (n = 44) HIV-infected women in absence of treatment with antiretroviral therapy (ART). Independent of HIV-1 subtype, p-distance, dN and dS were all strongly correlated with one another but were not significantly different in pregnant as compared to non-pregnant patients. Differential levels of selective pressure applied on different Env subdomains displayed similar yet non-identical patterns between the two groups, with pressure applied on C1 being significantly lower in constant regions C1 and C2 than in V1, V2, V3 and C3. To draw a general picture of the selection applied on the envelope and compensate for inter-individual variations, we performed a binomial test on selection frequency data pooled from pregnant and non-pregnant women. This analysis uncovered 42 positions, present in both groups, exhibiting statistically-significant frequency of selection that invariably mapped to the surface of the Env protein, with the great majority located within epitopes recognized by Env-specific antibodies or sites associated with the development of cross-reactive neutralizing activity. The median frequency of occurrence of positive selection per site was significantly lower in pregnant versus non-pregnant women. Furthermore, examination of the distribution of positively selected sites using a hypergeometric test revealed that only 2 positions (D137 and S142) significantly differed between the 2 groups. Taken together, these result indicate that pregnancy is associated with subtle yet distinctive changes in selective pressures exerted on the HIV-1 Env protein that are compatible with transient modulations of maternal immunity.

Impact of Clinical Parameters in the Intrahost Evolution of HIV-1 Subtype B in Pediatric Patients: A Machine Learning Approach

Determining the factors modulating the genetic diversity of HIV-1 populations is essential to understand viral evolution. This study analyzes the relative importance of clinical factors in the intrahost HIV-1 subtype B (HIV-1B) evolution and in the fixation of drug resistance mutations (DRM) during longitudinal pediatric HIV-1 infection. We recovered 162 partial HIV-1B pol sequences (from 3 to 24 per patient) from 24 perinatally infected patients from the Madrid Cohort of HIV-1 infected children and adolescents in a time interval ranging from 2.2 to 20.3 years. We applied machine learning classification methods to analyze the relative importance of 28 clinical/epidemiological/virological factors in the HIV-1B evolution to predict HIV-1B genetic diversity (d), nonsynonymous and synonymous mutations (dN, dS) and DRM presence. Most of the 24 HIV-1B infected pediatric patients were Spanish (91.7%), diagnosed before 2000 (83.3%), and all were antiretroviral therapy experienced. They had from 0.3 to 18.8 years of HIV-1 exposure at sampling time. Most sequences presented DRM. The best-predictor variables for HIV-1B evolutionary parameters were the age of HIV-1 diagnosis for d, the age at first antiretroviral treatment for dN and the year of HIV-1 diagnosis for ds. The year of infection (birth year) and year of sampling seemed to be relevant for fixation of both DRM at large and, considering drug families, to protease inhibitors (PI). This study identifies, for the first time using machine learning, the factors affecting more HIV-1B pol evolution and those affecting DRM fixation in HIV-1B infected pediatric patients.

Analysis of evolutionary rate of HIV-1 subtype B using blood donor samples in Japan

There are few reports on HIV-1 intra-host evolutionary rate in asymptomatic treatment-naïve patients. Here, the HIV-1 intra-host evolutionary rate was estimated based on HIV-1 RNA sequences from plasma samples of blood donors in Japan. Blood donors were assumed to have received no treatment for and have no symptoms of HIV-1 infection because they were healthy, and declared no risky behaviors of HIV-1 infection on a self-reported questionnaire or interview followed by donation. HIV-1 RNA was obtained from 85 plasma samples from 36 blood donors who donated blood multiple times and were HIV-1-positive. The C2V3C3 region which encodes for a part of the envelope protein, and the V3 loop in the C2V3C3 region were analyzed by RT-PCR and direct sequencing, and the sequences were compared. The nucleotide substitution rate was calculated by linear regression. All HIV-1 samples analyzed were classified as subtype B. The mean nucleotide substitution rate in C2V3C3 was calculated to be 6.2 × 10^-3-1.8 × 10^-2/site/year (V3: 4.5 × 10^-3-2.3 × 10^-2/site/year). The mean non-synonymous substitution rate in C2V3C3 was calculated to be 5.2 × 10^-3-1.7 × 10^-2/site/year (V3: 4.5 × 10^-3-2.1 × 10^-2/site/year). The mean synonymous substitution rate in C2V3C3 was calculated to be 1.1 × 10^-4-2.3 × 10^-3/site/year (V3: 2.9 × 10^-3/site/year). Among HIV-1 subtype B RNA-positive blood donors in Japan, the nucleotide substitution rate in C2V3C3 was estimated to be higher than that of reported cases using HIV-1 samples mainly obtained from AIDS patients. Compared to AIDS patients, immune responses against HIV-1 are probably more effective in HIV-1 RNA-positive blood donors. Consequently, immune pressure presumably promotes mutation of the virus genome.

The impact of HIV-1 within-host evolution on transmission dynamics

The adaptive potential of HIV-1 is a vital mechanism to evade host immune responses and antiviral treatment. However, high evolutionary rates during persistent infection can impair transmission efficiency and alter disease progression in the new host, resulting in a delicate trade-off between within-host virulence and between-host infectiousness. This trade-off is visible in the disparity in evolutionary rates at within-host and between-host levels, and preferential transmission of ancestral donor viruses. Understanding the impact of within-host evolution for epidemiological studies is essential for the design of preventive and therapeutic measures. Herein, we review recent theoretical and experimental work that generated new insights into the complex link between within-host evolution and between-host fitness, revealing temporal and selective processes underlying the structure and dynamics of HIV-1 transmission.

Patterns and rates of viral evolution in HIV-1 subtype B infected females and males

Biological sex differences affect the course of HIV infection, with untreated women having lower viral loads compared to their male counterparts but, for a given viral load, women have a higher rate of progression to AIDS. However, the vast majority of data on viral evolution, a process that is clearly impacted by host immunity and could be impacted by sex differences, has been derived from men. We conducted an intensive analysis of HIV-1 gag and env-gp120 evolution taken over the first 6–11 years of infection from 8 Women’s Interagency HIV Study (WIHS) participants who had not received combination antiretroviral therapy (ART). This was compared to similar data previously collected from men, with both groups infected with HIV-1 subtype B. Early virus populations in men and women were generally homogenous with no differences in diversity between sexes. No differences in ensuing nucleotide substitution rates were found between the female and male cohorts studied herein. As previously reported for men, time to peak diversity in env-gp120 in women was positively associated with time to CD4+ cell count below 200 (P = 0.017), and the number of predicted N-linked glycosylation sites generally increased over time, followed by a plateau or decline, with the majority of changes localized to the V1-V2 region. These findings strongly suggest that the sex differences in HIV-1 disease progression attributed to immune system composition and sensitivities are not revealed by, nor do they impact, global patterns of viral evolution, the latter of which proceeds similarly in women and men.

Mutation and recombination in pathogen evolution: Relevance, methods and controversies

Mutation and recombination drive the evolution of most pathogens by generating the genetic variants upon which selection operates. Those variants can, for example, confer resistance to host immune systems and drug therapies or lead to epidemic outbreaks. Given their importance, diverse evolutionary studies have investigated the abundance and consequences of mutation and recombination in pathogen populations. However, some controversies persist regarding the contribution of each evolutionary force to the development of particular phenotypic observations (e.g., drug resistance). In this study, we revise the importance of mutation and recombination in the evolution of pathogens at both intra-host and inter-host levels. We also describe state-of-the-art analytical methodologies to detect and quantify these two evolutionary forces, including biases that are often ignored in evolutionary studies. Finally, we present some of our former studies involving pathogenic taxa where mutation and recombination played crucial roles in the recovery of pathogenic fitness, the generation of interspecific genetic diversity, or the design of centralized vaccines. This review also illustrates several common controversies and pitfalls in the analysis and in the evaluation and interpretation of mutation and recombination outcomes.

HIV Progression Depends on Codon and Amino Acid Usage Profile of Envelope Protein and Associated Host-Genetic Influence

Acquired immune deficiency syndrome (AIDS) is a spectrum of conditions caused by infection with the human immunodeficiency virus (HIV). Two types of HIV have been characterized: HIV-1 and HIV-2. The present study investigated whether evolutionary selection pressure differs between rapid progressor (RP), slow progressor (SP), and long-term non-progressor (LTNP) of HIV-I infected individuals. An unexpected association between the evolutionary rate of substitution in envelope (env) gene and disease progression is observed. Our present study suggests that env genes of LTNP are subject to unusually strong functional constraint with respect to RP. We also observed that the three categories of env genes i.e., RP, SP, and LTNP, had their own characteristic pattern of amino acid usage and SP and LTNP sequences shared similar patterns of amino acid usage different from RP sequences and evolutionary rate significantly influenced the amino acid usage pattern of the three different types of env gene sequences. It was also noted that the evolutionary rate for the glycosylation sites of LTNP and SP sequences were even significantly less than the RP sequences. Comparative analysis on the influence of human host on the three categories of env genes are well correlated with the rates of disease progression suggesting the adaptive strategies of the viruses for successful residence and infection. Host associated selective constraints appeared most relaxed on the RP sequences and strongest in LTNP sequences. The present study clearly portrays how evolutionary selection pressure differs between three categories of env genes i.e., RP, SP, and LTNP. The env genes, coding for the env glycoproteins, experience severe selection constraints from the host due to their constant exposure to the host immune system. In this perspective it might be suggested that env gene evolution occurs mainly by negative selection with the occurrence of mutation that might not reach fixation in the viral population. This work also confers a deeper insight into the crucial effects of host factors that govern the overall progression of HIV infection.

Association of HIV diversity and virologic outcomes in early antiretroviral treatment: HPTN 052

Higher HIV diversity has been associated with virologic outcomes in children on antiretroviral treatment (ART). We examined the association of HIV diversity with virologic outcomes in adults from the HPTN 052 trial who initiated ART at CD4 cell counts of 350-550 cells/mm3. A high resolution melting (HRM) assay was used to analyze baseline (pre-treatment) HIV diversity in six regions in the HIV genome (two in gag, one in pol, and three in env) from 95 participants who failed ART. We analyzed the association of HIV diversity in each genomic region with baseline (pre-treatment) factors and three clinical outcomes: time to virologic suppression after ART initiation, time to ART failure, and emergence of HIV drug resistance at ART failure. After correcting for multiple comparisons, we did not find any association of baseline HIV diversity with demographic, laboratory, or clinical characteristics. For the 18 analyses performed for clinical outcomes evaluated, there was only one significant association: higher baseline HIV diversity in one of the three HIV env regions was associated with longer time to ART failure (p = 0.008). The HRM diversity assay may be useful in future studies exploring the relationship between HIV diversity and clinical outcomes in individuals with HIV infection.

Short-Sighted Virus Evolution and a Germline Hypothesis for Chronic Viral Infections

With extremely short generation times and high mutability, many viruses can rapidly evolve and adapt to changing environments. This ability is generally beneficial to viruses as it allows them to evade host immune responses, evolve new behaviours, and exploit ecological niches. However, natural selection typically generates adaptation in response to the immediate selection pressures that a virus experiences in its current host. Consequently, we argue that some viruses, particularly those characterised by long durations of infection and ongoing replication, may be susceptible to short-sighted evolution, whereby a virus' adaptation to its current host will be detrimental to its onward transmission within the host population. Here we outline the concept of short-sighted viral evolution and provide examples of how it may negatively impact viral transmission among hosts. We also propose that viruses that are vulnerable to short-sighted evolution may exhibit strategies that minimise its effects. We speculate on the various mechanisms by which this may be achieved, including viral life history strategies that result in low rates of within-host evolution, or the establishment of a 'germline' lineage of viruses that avoids short-sighted evolution. These concepts provide a new perspective on the way in which some viruses have been able to establish and maintain global pandemics.

Recent advances in understanding HIV evolution

The human immunodeficiency virus (HIV) evolves rapidly owing to the combined activity of error-prone reverse transcriptase, recombination, and short generation times, leading to extensive viral diversity both within and between hosts. This diversity is a major contributing factor in the failure of the immune system to eradicate the virus and has important implications for the development of suitable drugs and vaccines to combat infection. This review will discuss the recent technological advances that have shed light on HIV evolution and will summarise emerging concepts in this field.

Genotype-Specific Evolution of Hepatitis E Virus

Hepatitis E virus (HEV) is the most common cause of acute viral hepatitis globally. HEV comprises four genotypes with different geographic distributions and host ranges. We utilize this natural case-control study for investigating the evolution of zoonotic viruses compared to single-host viruses, using 244 near-full-length HEV genomes. Genome-wide estimates of the ratio of nonsynonymous to synonymous evolutionary changes (dN/dS ratio) located a region of overlapping reading frames, which is subject to positive selection in genotypes 3 and 4. The open reading frames (ORFs) involved have functions related to host-pathogen interaction, so genotype-specific evolution of these regions may reflect their fitness. Bayesian inference of evolutionary rates shows that genotypes 3 and 4 have significantly higher rates than genotype 1 across all ORFs. Reconstruction of the phylogenies of zoonotic genotypes demonstrates significant intermingling of isolates between hosts. We speculate that the genotype-specific differences may result from cyclical adaptation to different hosts in genotypes 3 and 4.IMPORTANCE Hepatitis E virus (HEV) is increasingly recognized as a pathogen that affects both the developing and the developed world. While most often clinically mild, HEV can be severe or fatal in certain demographics, such as expectant mothers. Like many other viral pathogens, HEV has been classified into several distinct genotypes. We show that most of the HEV genome is evolutionarily constrained. One locus of positive selection is unusual in that it encodes two distinct protein products. We are the first to detect positive selection in this overlap region. Genotype 1, which infects humans only, appears to be evolving differently from genotypes 3 and 4, which infect multiple species, possibly because genotypes 3 and 4 are unable to achieve the same fitness due to repeated host jumps.

Accurate predictions of population-level changes in sequence and structural properties of HIV-1 Env using a volatility-controlled diffusion model

The envelope glycoproteins (Envs) of HIV-1 continuously evolve in the host by random mutations and recombination events. The resulting diversity of Env variants circulating in the population and their continuing diversification process limit the efficacy of AIDS vaccines. We examined the historic changes in Env sequence and structural features (measured by integrity of epitopes on the Env trimer) in a geographically defined population in the United States. As expected, many Env features were relatively conserved during the 1980s. From this state, some features diversified whereas others remained conserved across the years. We sought to identify “clues” to predict the observed historic diversification patterns. Comparison of viruses that cocirculate in patients at any given time revealed that each feature of Env (sequence or structural) exists at a defined level of variance. The in-host variance of each feature is highly conserved among individuals but can vary between different HIV-1 clades. We designate this property “volatility” and apply it to model evolution of features as a linear diffusion process that progresses with increasing genetic distance. Volatilities of different features are highly correlated with their divergence in longitudinally monitored patients. Volatilities of features also correlate highly with their population-level diversification. Using volatility indices measured from a small number of patient samples, we accurately predict the population diversity that developed for each feature over the course of 30 years. Amino acid variants that evolved at key antigenic sites are also predicted well. Therefore, small “fluctuations” in feature values measured in isolated patient samples accurately describe their potential for population-level diversification. These tools will likely contribute to the design of population-targeted AIDS vaccines by effectively capturing the diversity of currently circulating strains and addressing properties of variants expected to appear in the future.

HIV-1 is the causative agent of the global AIDS pandemic. The envelope glycoproteins (Envs) of HIV-1 constitute a primary target for antibody-based vaccines. However, the diversity of Envs in the population limits the potential efficacy of this approach. Accurate estimates of the range of variants that currently infect patients and those expected to appear in the future will likely contribute to the design of population-targeted immunogens. We found that different properties (features) of Env have different propensities for small “fluctuations” in their values among viruses that infect patients at any given time point. This propensity of each feature for in-host variance, which we designate “volatility”, is conserved among patients. We apply this parameter to model the evolution of features (in patients and population) as a diffusion process driven by their “diffusion coefficients” (volatilities). Using volatilities measured from a few patient samples from the 1980s, we accurately predict properties of viruses that evolved in the population over the course of 30 years. The diffusion-based model described here efficiently captures evolution of phenotypes in biological systems controlled by a dominant random component.

Cyclophilins and nucleoporins are required for infection mediated by capsids from circulating HIV-2 primary isolates

HIV-2 groups have emerged from sooty mangabey SIV and entered the human population in Africa on several separate occasions. Compared to world pandemic HIV-1 that arose from the chimpanzee SIVcpz virus, the SIVsm-derived HIV-2, largely confined to West Africa, is less replicative, less transmissible and less pathogenic. Here, we evaluated the interactions between host cellular factors, which control HIV-1 infection and target the capsid, and HIV-2 capsids obtained from primary isolates from patients with different disease progression status. We showed that, like HIV-1, all HIV-2 CA we tested exhibited a dependence on cyclophilin A. However, we observed no correlation between HIV-2 viremia and susceptibility to hu-TRIM5alpha or dependence to CypA. Finally, we found that all CA from HIV-2 primary isolates exploit Nup358 and Nup153 for nucleus transposition. Altogether, these findings indicate that the ability to use the two latter nucleoporins is essential to infection of human cells for both HIV-1 and HIV-2. This dependence provides another molecular target that could be used for antiviral strategies against both HIV-1 and 2, based on both nucleoporins.

Fidelity of classwide-resistant HIV-2 reverse transcriptase and differential contribution of K65R to the accuracy of HIV-1 and HIV-2 reverse transcriptases

Nucleoside reverse transcriptase (RT) inhibitors constitute the backbone of current therapies against human immunodeficiency virus type 1 and type 2 (HIV-1 and HIV-2, respectively). However, mutational pathways leading to the development of nucleoside analogue resistance are different in both types of HIV. In HIV-2, resistance to all approved nucleoside analogues is conferred by the combination of RT substitutions K65R, Q151M and M184V. Nucleotide incorporation kinetic analyses of mutant and wild-type (WT) HIV-2 RTs show that the triple-mutant has decreased catalytic efficiency due to the presence of M184V. Although similar effects were previously reported for equivalent mutations in HIV-1 RT, the HIV-2 enzymes were catalytically less efficient. Interestingly, in highly divergent HIV-1 RTs, K65R confers several-fold increased accuracy of DNA synthesis. We have determined the intrinsic fidelity of DNA synthesis of WT HIV-2 RT and mutants K65R and K65R/Q151M/M184V. Our results show that those changes in HIV-2 RT have a relatively small impact on nucleotide selectivity. Furthermore, we found that there were less than two-fold differences in error rates obtained with forward mutation assays using mutant and WT HIV-2 RTs. A different conformation of the β3-β4 hairpin loop in HIV-1 and HIV-2 RTs could probably explain the differential effects of K65R.

Effect of the Latent Reservoir on the Evolution of HIV at the Within- and Between-Host Levels

The existence of long-lived reservoirs of latently infected CD4+ T cells is the major barrier to curing HIV, and has been extensively studied in this light. However, the effect of these reservoirs on the evolutionary dynamics of the virus has received little attention. Here, we present a within-host quasispecies model that incorporates a long-lived reservoir, which we then nest into an epidemiological model of HIV dynamics. For biologically plausible parameter values, we find that the presence of a latent reservoir can severely delay evolutionary dynamics within a single host, with longer delays associated with larger relative reservoir sizes and/or homeostatic proliferation of cells within the reservoir. These delays can fundamentally change the dynamics of the virus at the epidemiological scale. In particular, the delay in within-host evolutionary dynamics can be sufficient for the virus to evolve intermediate viral loads consistent with maximising transmission, as is observed, and not the very high viral loads that previous models have predicted, an effect that can be further enhanced if viruses similar to those that initiate infection are preferentially transmitted. These results depend strongly on within-host characteristics such as the relative reservoir size, with the evolution of intermediate viral loads observed only when the within-host dynamics are sufficiently delayed. In conclusion, we argue that the latent reservoir has important, and hitherto under-appreciated, roles in both within- and between-host viral evolution.

Differential Infection Patterns and Recent Evolutionary Origins of Equine Hepaciviruses in Donkeys

The hepatitis C virus (HCV) is a major human pathogen. Genetically related viruses in animals suggest a zoonotic origin of HCV. The closest relative of HCV is found in horses (termed equine hepacivirus [EqHV]). However, low EqHV genetic diversity implies relatively recent acquisition of EqHV by horses, making a derivation of HCV from EqHV unlikely. To unravel the EqHV evolutionary history within equid sister species, we analyzed 829 donkeys and 53 mules sampled in nine European, Asian, African, and American countries by molecular and serologic tools for EqHV infection. Antibodies were found in 278 animals (31.5%), and viral RNA was found in 3 animals (0.3%), all of which were simultaneously seropositive. A low RNA prevalence in spite of high seroprevalence suggests a predominance of acute infection, a possible difference from the mostly chronic hepacivirus infection pattern seen in horses and humans. Limitation of transmission due to short courses of infection may explain the existence of entirely seronegative groups of animals. Donkey and horse EqHV strains were paraphyletic and 97.5 to 98.2% identical in their translated polyprotein sequences, making virus/host cospeciation unlikely. Evolutionary reconstructions supported host switches of EqHV between horses and donkeys without the involvement of adaptive evolution. Global admixture of donkey and horse hepaciviruses was compatible with anthropogenic alterations of EqHV ecology. In summary, our findings do not support EqHV as the origin of the significantly more diversified HCV. Identification of a host system with predominantly acute hepacivirus infection may enable new insights into the chronic infection pattern associated with HCV.

IMPORTANCE

The evolutionary origins of the human hepatitis C virus (HCV) are unclear. The closest animal-associated relative of HCV occurs in horses (equine hepacivirus [EqHV]). The low EqHV genetic diversity implies a relatively recent acquisition of EqHV by horses, limiting the time span for potential horse-to-human infections in the past. Horses are genetically related to donkeys, and EqHV may have cospeciated with these host species. Here, we investigated a large panel of donkeys from various countries using serologic and molecular tools. We found EqHV to be globally widespread in donkeys and identify potential differences in EqHV infection patterns, with donkeys potentially showing enhanced EqHV clearance compared to horses. We provide strong evidence against EqHV cospeciation and for its capability to switch hosts among equines. Differential hepacivirus infection patterns in horses and donkeys may enable new insights into the chronic infection pattern associated with HCV.

Analysis of HIV Diversity in HIV-Infected Black Men Who Have Sex with Men (HPTN 061)

Background

HIV populations often diversify in response to selective pressures, such as the immune response and antiretroviral drug use. We analyzed HIV diversity in Black men who have sex with men who were enrolled in the HIV Prevention Trials Network 061 study.

Methods

A high resolution melting (HRM) diversity assay was used to measure diversity in six regions of the HIV genome: two in gag, one in pol, and three in env. HIV diversity was analyzed for 146 men who were HIV infected at study enrollment, including three with acute infection and 13 with recent infection (identified using a multi-assay algorithm), and for 21 men who seroconverted during the study. HIV diversification was analyzed in a paired analysis for 62 HIV-infected men using plasma samples from the enrollment and 12-month (end of study) visits.

Results

Men with acute or recent infection at enrollment and seroconverters had lower median HRM scores (lower HIV diversity) than men with non-recent infection in all six regions analyzed. In univariate analyses, younger age, higher CD4 cell count, and HIV drug resistance were associated with lower median HRM scores in multiple regions; ARV drug detection was marginally associated with lower diversity in the pol region. In multivariate analysis, acute or recent infection (all six regions) and HIV drug resistance (both gag regions) were associated with lower median HRM scores. Diversification in the pol region over 12 months was greater for men with acute or recent infection, higher CD4 cell count, and lower HIV viral load at study enrollment.

Conclusions

HIV diversity was significantly associated with duration of HIV infection, and lower gag diversity was observed in men who had HIV drug resistance. HIV pol diversification was more pronounced in men with acute or recent infection, higher CD4 cell count, and lower HIV viral load.

Clinical Determinants of HIV-1B Between-Host Evolution and their Association with Drug Resistance in Pediatric Patients

Understanding the factors that modulate the evolution of virus populations is essential to design efficient control strategies. Mathematical models predict that factors affecting viral within-host evolution may also determine that at the between-host level. Although HIV-1 within-host evolution has been associated with clinical factors used to monitor AIDS progression, such as patient age, CD4 cells count, viral load, and antiretroviral experience, little is known about the role of these clinical factors in determining between-host HIV-1 evolution. Moreover, whether the relative importance of such factors in HIV-1 evolution vary in adult and children patients, in which the course of infection is different, has seldom been analysed. To address these questions, HIV-1 subtype B (HIV-1B) pol sequences of 163 infected children and 450 adults of Madrid, Spain, were used to estimate genetic diversity, rates of synonymous and non-synonymous mutations, selection pressures and frequency of drug-resistance mutations (DRMs). The role and relative importance of patient age, %CD4, CD4/mm3, viral load, and antiretroviral experience in HIV-1B evolution was analysed. In the pediatric HIV-1B population, three clinical factors were primary predictors of virus evolution: Higher HIV-1B genetic diversity was observed with increasing children age, decreasing CD4/mm3 and upon antiretroviral experience. This was mostly due to higher rates of non-synonymous mutations, which were associated with higher frequency of DRMs. Using this data, we have also constructed a simple multivariate model explaining between 55% and 66% of the variance in HIV-1B evolutionary parameters in pediatric populations. On the other hand, the analysed clinical factors had little effect in adult-infecting HIV-1B evolution. These findings highlight the different evolutionary dynamics of HIV-1B in children and adults, and contribute to understand the factors shaping HIV-1B evolution and the appearance of drug-resistance mutation in pediatric patients.

Exceptional Heterogeneity in Viral Evolutionary Dynamics Characterises Chronic Hepatitis C Virus Infection

The treatment of HCV infection has seen significant progress, particularly since the approval of new direct-acting antiviral drugs. However these clinical achievements have been made despite an incomplete understanding of HCV replication and within-host evolution, especially compared with HIV-1. Here, we undertake a comprehensive analysis of HCV within-host evolution during chronic infection by investigating over 4000 viral sequences sampled longitudinally from 15 HCV-infected patients. We compare our HCV results to those from a well-studied HIV-1 cohort, revealing key differences in the evolutionary behaviour of these two chronic-infecting pathogens. Notably, we find an exceptional level of heterogeneity in the molecular evolution of HCV, both within and among infected individuals. Furthermore, these patterns are associated with the long-term maintenance of viral lineages within patients, which fluctuate in relative frequency in peripheral blood. Together, our findings demonstrate that HCV replication behavior is complex and likely comprises multiple viral subpopulations with distinct evolutionary dynamics. The presence of a structured viral population can explain apparent paradoxes in chronic HCV infection, such as rapid fluctuations in viral diversity and the reappearance of viral strains years after their initial detection.

Our knowledge of HCV within-host evolution is substantially limited, which is surprising given that highly successful therapies against the virus have been developed. Key aspects of HCV infection, such as rapid fluctuations in viral diversity and the reappearance of viral strains years after their initial detection, remain unexplained. To better understand this problem, we analyse viral sequences from HCV-infected patients sampled over several years. Our findings suggest that the replication dynamics during chronic HCV infection are distinct from those of HIV-1, and dominated by the co-circulation of multiple viral strains. Although a major difference between the two chronic-infecting viruses is the level of recombination, our results indicate that HCV within-host evolution is most likely to be shaped by a structured viral population. Crucially, our study shows that HCV sampled from blood does not fully represent the within-host viral population at that time. This may have important implications for HCV treatment, especially in patients that have seemingly cleared the virus, as well as for molecular epidemiology studies investigating HCV transmission.

Identifying predictors of time-inhomogeneous viral evolutionary processes

Various factors determine the rate at which mutations are generated and fixed in viral genomes. Viral evolutionary rates may vary over the course of a single persistent infection and can reflect changes in replication rates and selective dynamics. Dedicated statistical inference approaches are required to understand how the complex interplay of these processes shapes the genetic diversity and divergence in viral populations. Although evolutionary models accommodating a high degree of complexity can now be formalized, adequately informing these models by potentially sparse data, and assessing the association of the resulting estimates with external predictors, remains a major challenge. In this article, we present a novel Bayesian evolutionary inference method, which integrates multiple potential predictors and tests their association with variation in the absolute rates of synonymous and non-synonymous substitutions along the evolutionary history. We consider clinical and virological measures as predictors, but also changes in population size trajectories that are simultaneously inferred using coalescent modelling. We demonstrate the potential of our method in an application to within-host HIV-1 sequence data sampled throughout the infection of multiple patients. While analyses of individual patient populations lack statistical power, we detect significant evidence for an abrupt drop in non-synonymous rates in late stage infection and a more gradual increase in synonymous rates over the course of infection in a joint analysis across all patients. The former is predicted by the immune relaxation hypothesis while the latter may be in line with increasing replicative fitness during the asymptomatic stage.

Viral Evolution and Cytotoxic T Cell Restricted Selection in Acute Infant HIV-1 Infection

Antiretroviral therapy-naive HIV-1 infected infants experience poor viral containment and rapid disease progression compared to adults. Viral factors (e.g. transmitted cytotoxic T- lymphocyte (CTL) escape mutations) or infant factors (e.g. reduced CTL functional capacity) may explain this observation. We assessed CTL functionality by analysing selection in CTL-targeted HIV-1 epitopes following perinatal infection. HIV-1 gag, pol and nef sequences were generated from a historical repository of longitudinal specimens from 19 vertically infected infants. Evolutionary rate and selection were estimated for each gene and in CTL-restricted and non-restricted epitopes. Evolutionary rate was higher in nef and gag vs. pol, and lower in infants with non-severe immunosuppression vs. severe immunosuppression across gag and nef. Selection pressure was stronger in infants with non-severe immunosuppression vs. severe immunosuppression across gag. The analysis also showed that infants with non-severe immunosuppression had stronger selection in CTL-restricted vs. non-restricted epitopes in gag and nef. Evidence of stronger CTL selection was absent in infants with severe immunosuppression. These data indicate that infant CTLs can exert selection pressure on gag and nef epitopes in early infection and that stronger selection across CTL epitopes is associated with favourable clinical outcomes. These results have implications for the development of paediatric HIV-1 vaccines.

Bayesian codon substitution modelling to identify sources of pathogen evolutionary rate variation

Phylodynamic reconstructions rely on a measurable molecular footprint of epidemic processes in pathogen genomes. Identifying the factors that govern the tempo and mode by which these processes leave a footprint in pathogen genomes represents an important goal towards understanding infectious disease evolution. Discriminating between synonymous and non-synonymous substitution rates is crucial for testing hypotheses about the sources of evolutionary rate variation. Here, we implement a codon substitution model in a Bayesian statistical framework to estimate absolute rates of synonymous and non-synonymous substitution in unknown evolutionary histories. To demonstrate how this model can provide critical insights into pathogen evolutionary dynamics, we adopt hierarchical phylogenetic modelling with fixed effects and apply it to two viral examples. Using within-host HIV-1 data from patients with different host genetic background and different disease progression rates, we show that viral populations undergo faster absolute synonymous substitution rates in patients with faster disease progression, probably reflecting faster replication rates. We also re-analyse rabies data from different bat species in the Americas to demonstrate that climate predicts absolute synonymous substitution rates, which can be attributed to climate-associated bat activity and viral transmission dynamics. In conclusion, our model to estimate absolute rates of synonymous and non-synonymous substitution can provide a powerful approach to investigate how host ecology can shape the tempo of pathogen evolution.

Molecular clock of HIV-1 envelope genes under early immune selection

Background

The molecular clock hypothesis that genes or proteins evolve at a constant rate is a key tool to reveal phylogenetic relationships among species. Using the molecular clock, we can trace an infection back to transmission using HIV-1 sequences from a single time point. Whether or not a strict molecular clock applies to HIV-1’s early evolution in the presence of immune selection has not yet been fully examined.

Results

We identified molecular clock signatures from 1587 previously published HIV-1 full envelope gene sequences obtained since acute infection in 15 subjects. Each subject’s sequence diversity linearly increased during the first 150 days post infection, with rates ranging from \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.54 \times 10^{ - 5}$$\end{document}1.54×10-5 to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3.91 \times 10^{ - 5}$$\end{document}3.91×10-5 with a mean of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2.69 \times 10^{ - 5}$$\end{document}2.69×10-5 per base per day. The rate of diversification for 12 out of the 15 subjects was comparable to the neutral evolution rate. While temporal diversification was consistent with evolution patterns in the absence of selection, mutations from the founder virus were highly clustered on statistically identified selection sites, which diversified more than 65 times faster than non-selection sites. By mathematically quantifying deviations from the molecular clock under various selection scenarios, we demonstrate that the deviation from a constant clock becomes negligible as multiple escape lineages emerge. The most recent common ancestor of a virus pair from distinct escape lineages is most likely the transmitted founder virus, indicating that HIV-1 molecular dating is feasible even after the founder viruses are no longer detectable.

Conclusions

The ability of HIV-1 to escape from immune surveillance in many different directions is the driving force of molecular clock persistence. This finding advances our understanding of the robustness of HIV-1’s molecular clock under immune selection, implying the potential for molecular dating.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-016-0269-6) contains supplementary material, which is available to authorized users.

Intra-epidemic evolutionary dynamics of a Dengue virus type 1 population reveal mutant spectra that correlate with disease transmission

Dengue virus (DENV) is currently the most prevalent mosquito-borne viral pathogen. DENVs naturally exist as highly heterogeneous populations. Even though the descriptions on DENV diversity are plentiful, only a few studies have narrated the dynamics of intra-epidemic virus diversity at a fine scale. Such accounts are important to decipher the reciprocal relationship between viral evolutionary dynamics and disease transmission that shape dengue epidemiology. In the current study, we present a micro-scale genetic analysis of a monophyletic lineage of DENV-1 genotype III (epidemic lineage) detected from November 2012 to May 2014. The lineage was involved in an unprecedented dengue epidemic in Singapore during 2013–2014. Our findings showed that the epidemic lineage was an ensemble of mutants (variants) originated from an initial mixed viral population. The composition of mutant spectrum was dynamic and positively correlated with case load. The close interaction between viral evolution and transmission intensity indicated that tracking genetic diversity through time is potentially a useful tool to infer DENV transmission dynamics and thereby, to assess the epidemic risk in a disease control perspective. Moreover, such information is salient to understand the viral basis of clinical outcome and immune response variations that is imperative to effective vaccine design.