The role of recombination for the coevolutionary dynamics of HIV and the immune response

Social, Epidemiological, and Virological Characteristics from Peruvian Subjects Living with HIV-1/AIDS with Different Sexual Risk Behavior

HIV-1 genetic diversity and resistance profile might change according to the risky sexual behavior of the host. To show this, we recruited 134 individuals between the years 2015 and 2017 identified as transgender women sex workers (TWSW, n = 73) and Heterosexual Military Officers (HET-MO, n = 61). After obtaining informed consent, we collected a blood sample to perform the HIV genotyping, CD4 cell count, and viral load. We used bioinformatics approaches for detecting resistance mutations and recombination events. Epidemiological data showed that both groups reported sexually transmitted diseases and they were widespread among TWSW, especially syphilis and herpes virus (35.6%). Illegal drugs consumption was higher among TWSW (71.2%), whereas condom use was inconsistent for both HET-MO (57.4%) and TWSW (74.0%). TWSW showed the shortest time exposition to antiretroviral therapy (ART) (3.5 years) and the lowest access to ART (34.2%) that conducted treatment failure (>4 logs). HIV-1 sequences from TWSW and HET-MO were analyzed to determine the genetic diversity and antiretroviral drug resistance. Phylogeny analysis revealed 125 (93%) cases of subtype B, 01 subtype A (0.76%), 07 (5.30%) BF recombinants, and 01 (0.76%) AG recombinant. Also, TWSW showed a higher recombination index (9.5%, 7/73) than HET-MO (1.5%, 1/68). HET-MO only showed acquired resistance (26.23%, 16/61), whereas TWSW showed both acquired as transmitted resistance (9.59% for each). In conclusion, TWSW and HET-MO showed significant differences considering the epidemiological characteristics, genetic diversity, recombination events, and HIV resistance profile.

Genetic Diversity of HIV-1 in Krasnoyarsk Krai: Area with High Levels of HIV-1 Recombination in Russia

More than a quarter of HIV-infected individuals registered in Russia live in Siberia. Unlike Central Russia where HIV-1 subtype A6 is predominant, in most Siberian regions since 2012, a new HIV-1 CRF63_02A1 genetic variant has spread, with the share of this variant attaining 75-85% among newly identified HIV cases. Krasnoyarsk Krai is considered to be a high-risk territory according to morbidity rate and HIV infection incidence among the population. The current paper aims to study the molecular epidemiologic characteristics of HIV-1 spreading in Krasnoyarsk Krai. Phylogenetic and recombination analyses of pol (PR-RT, IN) and env regions of the virus were used for genotyping 159 HIV-1 isolated in Krasnoyarsk Krai. 57.2% of the isolates belonged to subtype A (A6) specific to Russia, 12.6% to CRF63_02A1, and 0.6% to CRF02_AGСА, and in 29.6% HIV-1 URFs were detected, including URF63/А (23.9%), URFА/В (4.4%), and URF02/А (1.3%). In 6 of 7, HIV-1 URFА/В identical recombination model was detected; the origin of 38 URF63/А was proven to be the result of individual recombination events. Since 2015, a share of the population with newly diagnosed HIV who were infected with HIV-1 URF reached an exceptionally high rate of 38.6%. As distinct from adjacent Siberian regions, the HIV-1 CRF63_02A1 prevalence rate in Krasnoyarsk Krai is within 16%; however, the increased contribution of new HIV-1 into the regional epidemic development was observed due to the recombination of viruses of subtypes А, В, and CRF63_02A1. The difference between the described molecular epidemiologic picture in Krasnoyarsk Krai and in adjacent areas is likely caused by differences in predominant routes of HIV transmission and by more recent HIV-1 CRF63_02A1 transmission in the PWID group, which had a high prevalence of HIV-1 subtype A by the time of the new virus transmission, resulting in increased possibility of coinfection with various HIV-1 genetic variants.

Effect of synaptic cell-to-cell transmission and recombination on the evolution of double mutants in HIV

Recombination in HIV infection can impact virus evolution in vivo in complex ways, as has been shown both experimentally and mathematically. The effect of free virus versus synaptic, cell-to-cell transmission on the evolution of double mutants, however, has not been investigated. Here, we do so by using a stochastic agent-based model. Consistent with data, we assume spatial constraints for synaptic but not for free-virus transmission. Two important effects of the viral spread mode are observed: (i) for disadvantageous mutants, synaptic transmission protects against detrimental effects of recombination on double mutant persistence. Under free virus transmission, recombination increases double mutant levels for negative epistasis, but reduces them for positive epistasis. This reduction for positive epistasis is much diminished under predominantly synaptic transmission, and recombination can, in fact, lead to increased mutant levels. (ii) The mode of virus spread also directly influences the evolutionary fate of double mutants. For disadvantageous mutants, double mutant production is the predominant driving force, and hence synaptic transmission leads to highest double mutant levels due to increased transmission efficiency. For advantageous mutants, double mutant spread is the most important force, and hence free virus transmission leads to fastest invasion due to better mixing. For neutral mutants, both production and spread of double mutants are important, and hence an optimal mixture of free virus and synaptic transmission maximizes double mutant fractions. Therefore, both free virus and synaptic transmission can enhance or delay double mutant evolution. Implications for drug resistance in HIV are discussed.

Within-Epitope Interactions Can Bias CTL Escape Estimation in Early HIV Infection

As human immunodeficiency virus (HIV) begins to replicate within hosts, immune responses are elicited against it. Escape mutations in viral epitopes—immunogenic peptide parts presented on the surface of infected cells—allow HIV to partially evade these responses, and thus rapidly go to fixation. The faster they go to fixation, i.e., the higher their escape rate, the larger the selective pressure exerted by the immune system is assumed to be. This relation underpins the rationale for using escapes to assess the strength of immune responses. However, escape rate estimates are often obtained by employing an aggregation procedure, where several mutations that affect the same epitope are aggregated into a single, composite epitope mutation. The aggregation procedure thus rests upon the assumption that all within-epitope mutations have indistinguishable effects on immune recognition. In this study, we investigate how violation of this assumption affects escape rate estimates. To this end, we extend a previously developed simulation model of HIV that accounts for mutation, selection, and recombination to include different distributions of fitness effects (DFEs) and inter-mutational genomic distances. We use this discrete time Wright–Fisher based model to simulate early within-host evolution of HIV for DFEs and apply standard estimation methods to infer the escape rates. We then compare true with estimated escape rate values. We also compare escape rate values obtained by applying the aggregation procedure with values estimated without use of that procedure. We find that across the DFEs analyzed, the aggregation procedure alters the detectability of escape mutations: large-effect mutations are overrepresented while small-effect mutations are concealed. The effect of the aggregation procedure is similar to extracting the largest-effect mutation appearing within an epitope. Furthermore, the more pronounced the over-exponential decay of the DFEs, the more severely true escape rates are underestimated. We conclude that the aggregation procedure has two main consequences. On the one hand, it leads to a misrepresentation of the DFE of fixed mutations. On the other hand, it conceals within-epitope interactions that may generate irregularities in mutation frequency trajectories that are thus left unexplained.

JAK-STAT Signaling Pathways and Inhibitors Affect Reversion of Envelope-Mutated HIV-1

HIV can spread by both cell-free and cell-to-cell transmission. Here, we show that many of the amino acid changes in Env that are close to the CD4 binding pocket can affect HIV replication. We generated a number of mutant viruses that were unable to infect T cells as cell-free viruses but were nevertheless able to infect certain T cell lines as cell-associated viruses, which was followed by reversion to the wild type. However, the activation of JAK-STAT signaling pathways caused the inhibition of such cell-to-cell infection as well as the reversion of multiple HIV Env mutants that displayed differences in their abilities to bind to the CD4 receptor. Specifically, two T cell activators, interleukin-2 (IL-2) and phorbol 12-myristate 13-acetate (PMA), both capable of activation of JAK-STAT pathways, were able to inhibit cell-to-cell viral transmission. In contrast, but consistent with the above result, a number of JAK-STAT and mTOR inhibitors actually promoted HIV-1 transmission and reversion. Hence, JAK-STAT signaling pathways may differentially affect the replication of a variety of HIV Env mutants in ways that differ from the role that these pathways play in the replication of wild-type viruses.IMPORTANCE Specific alterations in HIV Env close to the CD4 binding site can differentially change the ability of HIV to mediate infection for cell-free and cell-associated viruses. However, such differences are dependent to some extent on the types of target cells used. JAK-STAT signaling pathways are able to play major roles in these processes. This work sheds new light on factors that can govern HIV infection of target cells.

High-Sequence Diversity and Rapid Virus Turnover Contribute to Higher Rates of Coreceptor Switching in Treatment-Experienced Subjects with HIV-1 Viremia

Coreceptor switching from CCR5 to CXCR4 is common during chronic HIV-1 infection, but is even more common in individuals who have failed antiretroviral therapy (ART). Prior studies have suggested rapid mutation and/or recombination of HIV-1 envelope (env) genes during coreceptor switching. We compared the functional and genotypic changes in env of viruses from viremic subjects who had failed ART just before and after coreceptor switching and compared those to viruses from matched subjects without coreceptor switching. Analysis of multiple unique functional env clones from each subject revealed extensive diversity at both sample time points and rapid diversification of sequences during the 4-month interval in viruses from both 9 subjects with coreceptor switching and 15 control subjects. Only two subjects had envs with evidence of recombination. Three findings distinguished env clones from subjects with coreceptor switching from controls: (1) lower entry efficiency via CCR5; (2) longer V1/V2 regions; and (3), lower nadir CD4 T cell counts during prior years of infection. Most of these subjects harbored virus with lower replicative capacity associated with protease (PR) and/or reverse transcriptase inhibitor resistance mutations, and the extensive diversification tended to lead either to improved entry efficiency via CCR5 or the gain of entry function via CXCR4. These results suggest that R5X4 or X4 variants emerge from a diverse, low-fitness landscape shaped by chronic infection, multiple ART resistance mutations, the availability of target cells, and reduced entry efficiency via CCR5.

Models of Viral Population Dynamics

Models of viral population dynamics have contributed enormously to our understanding of the pathogenesis and transmission of several infectious diseases, the coevolutionary dynamics of viruses and their hosts, the mechanisms of action of drugs, and the effectiveness of interventions. In this chapter, we review major advances in the modeling of the population dynamics of the human immunodeficiency virus (HIV) and briefly discuss adaptations to other viruses.

Investigating the Consequences of Interference between Multiple CD8+ T Cell Escape Mutations in Early HIV Infection

During early human immunodeficiency virus (HIV) infection multiple CD8⁺ T cell responses are elicited almost simultaneously. These responses exert strong selective pressures on different parts of HIV’s genome, and select for mutations that escape recognition and are thus beneficial to the virus. Some studies reveal that the later these escape mutations emerge, the more slowly they go to fixation. This pattern of escape rate decrease(ERD) can arise by distinct mechanisms. In particular, in large populations with high beneficial mutation rates interference among different escape strains –an effect that can emerge in evolution with asexual reproduction and results in delayed fixation times of beneficial mutations compared to sexual reproduction– could significantly impact the escape rates of mutations. In this paper, we investigated how interference between these concurrent escape mutations affects their escape rates in systems with multiple epitopes, and whether it could be a source of the ERD pattern. To address these issues, we developed a multilocus Wright-Fisher model of HIV dynamics with selection, mutation and recombination, serving as a null-model for interference. We also derived an interference-free null model assuming initial neutral evolution before immune response elicitation. We found that interference between several equally selectively advantageous mutations can generate the observed ERD pattern. We also found that the number of loci, as well as recombination rates substantially affect ERD. These effects can be explained by the underexponential decline of escape rates over time. Lastly, we found that the observed ERD pattern in HIV infected individuals is consistent with both independent, interference-free mutations as well as interference effects. Our results confirm that interference effects should be considered when analyzing HIV escape mutations. The challenge in estimating escape rates and mutation-associated selective coefficients posed by interference effects cannot simply be overcome by improved sampling frequencies or sizes. This problem is a consequence of the fundamental shortcomings of current estimation techniques under interference regimes. Hence, accounting for the stochastic nature of competition between mutations demands novel estimation methodologies based on the analysis of HIV strains, rather than mutation frequencies.

Within-host evolution of human immunodeficiency virus (HIV) is shaped by strong immune responses mounted against the virus. Multiple CD8⁺ T cell populations, each recognizing a specific part of an HIV protein, simultaneously suppress HIV growth. Escape mutations that arise in HIV genome regions coding for these virion protein parts, impair CD8⁺ T cell recognition and are consequently strongly selected. The emergence and rise of these escape mutations exhibits an intriguing temporal pattern: the earlier an escape mutation arises, the faster it goes to fixation. This pattern is termed escape rate decrease (ERD). In this paper, we tested computationally whether interference, i.e. the coexistence of multiple genetically distinct HIV strains engaged in competitive interaction within the host, could be a possible source of ERD. As an alternative, we also mathematically derived the temporal pattern of escapes under interference-free conditions, and compared this with data. We found that interference between multiple beneficial mutations could generate ERD. However, ERD does not imply the presence of interference. Thus, more detailed data is required to unambiguously determine whether interference effects influence ERD generation. Nevertheless, interference should be considered when studying the within-host evolution of HIV. Ignoring its effects on population dynamics can severely underestimate the protective capacity of CD8⁺ T cells.

Population genomics of intrapatient HIV-1 evolution

Many microbial populations rapidly adapt to changing environments with multiple variants competing for survival. To quantify such complex evolutionary dynamics in vivo, time resolved and genome wide data including rare variants are essential. We performed whole-genome deep sequencing of HIV-1 populations in 9 untreated patients, with 6-12 longitudinal samples per patient spanning 5-8 years of infection. The data can be accessed and explored via an interactive web application. We show that patterns of minor diversity are reproducible between patients and mirror global HIV-1 diversity, suggesting a universal landscape of fitness costs that control diversity. Reversions towards the ancestral HIV-1 sequence are observed throughout infection and account for almost one third of all sequence changes. Reversion rates depend strongly on conservation. Frequent recombination limits linkage disequilibrium to about 100 bp in most of the genome, but strong hitch-hiking due to short range linkage limits diversity.

The Effect of Interference on the CD8(+) T Cell Escape Rates in HIV

In early human immunodeficiency virus (HIV) infection, the virus population escapes from multiple CD8⁺ cell responses. The later an escape mutation emerges, the slower it outgrows its competition, i.e., the escape rate is lower. This pattern could indicate that the strength of the CD8⁺ cell responses is waning, or that later viral escape mutants carry a larger fitness cost. In this paper, we investigate whether the pattern of decreasing escape rates could also be caused by genetic interference among different escape strains. To this end, we developed a mathematical multi-epitope model of HIV dynamics, which incorporates stochastic effects, recombination, and mutation. We used cumulative linkage disequilibrium measures to quantify the amount of interference. We found that nearly synchronous, similarly strong immune responses in two-locus systems enhance the generation of genetic interference. This effect, combined with a scheme of densely spaced sampling times at the beginning of infection and sparse sampling times later, leads to decreasing successive escape rate estimates, even when there were no selection differences among alleles. These predictions are supported by empirical data from one HIV-infected patient. Thus, interference could explain why later escapes are slower. Considering escape mutations in isolation, neglecting their genetic linkage, conceals the underlying haplotype dynamics and can affect the estimation of the selective pressure exerted by CD8⁺ cells. In systems in which multiple escape mutations appear, the occurrence of interference dynamics should be assessed by measuring the linkage between different escape mutations.

The route of HIV escape from immune response targeting multiple sites is determined by the cost-benefit tradeoff of escape mutations

Cytotoxic T lymphocytes (CTL) are a major factor in the control of HIV replication. CTL arise in acute infection, causing escape mutations to spread rapidly through the population of infected cells. As a result, the virus develops partial resistance to the immune response. The factors controlling the order of mutating epitope sites are currently unknown and would provide a valuable tool for predicting conserved epitopes. In this work, we adapt a well-established mathematical model of HIV evolution under dynamical selection pressure from multiple CTL clones to include partial impairment of CTL recognition, , as well as cost to viral replication, . The process of escape is described in terms of the cost-benefit tradeoff of escape mutations and predicts a trajectory in the cost-benefit plane connecting sequentially escaped sites, which moves from high recognition loss/low fitness cost to low recognition loss/high fitness cost and has a larger slope for early escapes than for late escapes. The slope of the trajectory offers an interpretation of positive correlation between fitness costs and HLA binding impairment to HLA-A molecules and a protective subset of HLA-B molecules that was observed for clinically relevant escape mutations in the Pol gene. We estimate the value of from published experimental studies to be in the range (0.01–0.86) and show that the assumption of complete recognition loss () leads to an overestimate of mutation cost. Our analysis offers a consistent interpretation of the commonly observed pattern of escape, in which several escape mutations are observed transiently in an epitope. This non-nested pattern is a combined effect of temporal changes in selection pressure and partial recognition loss. We conclude that partial recognition loss is as important as fitness loss for predicting the order of escapes and, ultimately, for predicting conserved epitopes that can be targeted by vaccines.

Like many viruses, HIV has evolved mechanisms to evade the host immune response. As early as a few weeks after infection is initiated, mutations appear in the viral genome that reduce the ability of cytotoxic T lymphocytes (CTL) to control virus replication. However, of the many mutations in the viral genome that could potentially mediate viral escape from the CTL response, a specific subset are typically observed. This suggests that some mutations either entail too high a fitness cost for the virus, or are relatively inefficient escape mutations. A successful vaccine would target the CTL response to these regions in such a way that escape would not be possible. We use a computational model of HIV infection in order to study the factors that determine whether a given escape mutation will occur, how long it will be maintained in the population, and how these changes in the viral genome will affect the CTL response. Our analysis highlights the important role of partial recognition loss conferred by a mutation in producing the complex dynamics of escape that are observed during the course of infection.

Identification of an env-defective HIV-1 mutant capable of spontaneous reversion to a wild-type phenotype in certain T-cell lines

Background

Attempts to eradicate HIV from cellular reservoirs are vital but depend on a clear understanding of how viral variants are transmitted and survive in the different cell types that constitute such reservoirs. Mutations in the env gene of HIV may be able to exert a differential influence on viral transmission ability in regard to cell-free and cell-associated viral forms.

Methods

The ability of HIV containing an env G367R mutation in cell-free and cell-associated viruses to cause infection and to revert to wild-type was measured using several T cell lines. To determine factors that might potentially influence the reversion of G367R, we studied each of entry inhibitors, inhibitors of cellular endocytosis, and modulators of cell growth and activation.

Results

We demonstrate that an HIV-1 variant containing a G367R substitution within the CD4 binding site of gp120 was non-infectious as free virus in culture but was infectious when infected cells were co-cultured with certain T cell lines or when cells were transfected by a relevant proviral plasmid. Differences in viral infectivity by cell-associated G367R viruses were determined by the type of target cell employed, regardless which type of donor cell was used. Reversion was slowed or inhibited by entry inhibitors and by inhibitors of cellular endocytosis. Interleukin 2 was able to block G367R reversion in only one of the T cell lines studied but not in the other, while phorbol 12-myristate 13-acetate (PMA) inhibited G367R reversion in all the T cell lines.

Conclusions

Env-defective HIV may have a different phenotype as cell-free versus cell-associated virus. The persistence of defective forms can potentially lead to the emergence of virulent forms. The heterogeneity of cell types that constitute the HIV reservoir can contribute to viral variability, even among similar types of cells. This is the first demonstration of a mutation in the HIV envelope, i.e. G367R, that can compromise infection by cell-free virus but less severely by cell-associated virus and that does so in a cell type-dependent manner.

Recombination-mediated escape from primary CD8+ T cells in acute HIV-1 infection

Background

A major immune evasion mechanism of HIV-1 is the accumulation of non-synonymous mutations in and around T cell epitopes, resulting in loss of T cell recognition and virus escape.

Results

Here we analyze primary CD8+ T cell responses and virus escape in a HLA B*81 expressing subject who was infected with two T/F viruses from a single donor. In addition to classic escape through non-synonymous mutation/s, we also observed rapid selection of multiple recombinant viruses that conferred escape from T cells specific for two epitopes in Nef.

Conclusions

Our study shows that recombination between multiple T/F viruses provide greater options for acute escape from CD8+ T cell responses than seen in cases of single T/F virus infection. This process may contribute to the rapid disease progression in patients infected by multiple T/F viruses.

Electronic supplementary material

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

Identifying recombination hot spots in the HIV-1 genome

HIV-1 infection is characterized by the rapid generation of genetic diversity that facilitates viral escape from immune selection and antiretroviral therapy. Despite recombination's crucial role in viral diversity and evolution, little is known about the genomic factors that influence recombination between highly similar genomes. In this study, we use a minimally modified full-length HIV-1 genome and high-throughput sequence analysis to study recombination in gag and pol in T cells. We find that recombination is favored at a number of recombination hot spots, where recombination occurs six times more frequently than at corresponding cold spots. Interestingly, these hot spots occur near important features of the HIV-1 genome but do not occur at sites immediately around protease inhibitor or reverse transcriptase inhibitor drug resistance mutations. We show that the recombination hot and cold spots are consistent across five blood donors and are independent of coreceptor-mediated entry. Finally, we check common experimental confounders and find that these are not driving the location of recombination hot spots. This is the first study to identify the location of recombination hot spots between two similar viral genomes with great statistical power and under conditions that closely reflect natural recombination events among HIV-1 quasispecies.

IMPORTANCE

The ability of HIV-1 to evade the immune system and antiretroviral therapy depends on genetic diversity within the viral quasispecies. Retroviral recombination is an important mechanism that helps to generate and maintain this genetic diversity, but little is known about how recombination rates vary within the HIV-1 genome. We measured recombination rates in gag and pol and identified recombination hot and cold spots, demonstrating that recombination is not random but depends on the underlying gene sequence. The strength and location of these recombination hot and cold spots can be used to improve models of viral dynamics and evolution, which will be useful for the design of robust antiretroviral therapies.

VIRAPOPS: a forward simulator dedicated to rapidly evolved viral populations

SUMMARY

Daily, mutability and recombination of RNA viruses result in the production of million variants. All these rapid genomic changes directly influence the functional sites of the protein, its 3D structure or its drug resistances. Therefore, it is important to simulate these drastic switches to determine their effects on virus populations. Many computer programs are able to simulate specific variations in DNA genomes, but are generally non-adapted to RNA viruses. They simulate site-specific selection pressures, but rarely pressures on covariant or on higher order correlated sites and no at all on synthetic lethal groups. That is why we felt it important to create VIRAPOPS, a forward simulator that models specific RNA virus functions. It was designed for computational biologists, biologists and virologists.

AVAILABILITY AND IMPLEMENTATION

Free binaries are available through a software repository at http://petitjeanmichel.free.fr/itoweb.petitjean.freeware.html.

Superinfection with drug-resistant HIV is rare and does not contribute substantially to therapy failure in a large European cohort

Background

Superinfection with drug resistant HIV strains could potentially contribute to compromised therapy in patients initially infected with drug-sensitive virus and receiving antiretroviral therapy. To investigate the importance of this potential route to drug resistance, we developed a bioinformatics pipeline to detect superinfection from routinely collected genotyping data, and assessed whether superinfection contributed to increased drug resistance in a large European cohort of viremic, drug treated patients.

Methods

We used sequence data from routine genotypic tests spanning the protease and partial reverse transcriptase regions in the Virolab and EuResist databases that collated data from five European countries. Superinfection was indicated when sequences of a patient failed to cluster together in phylogenetic trees constructed with selected sets of control sequences. A subset of the indicated cases was validated by re-sequencing pol and env regions from the original samples.

Results

4425 patients had at least two sequences in the database, with a total of 13816 distinct sequence entries (of which 86% belonged to subtype B). We identified 107 patients with phylogenetic evidence for superinfection. In 14 of these cases, we analyzed newly amplified sequences from the original samples for validation purposes: only 2 cases were verified as superinfections in the repeated analyses, the other 12 cases turned out to involve sample or sequence misidentification. Resistance to drugs used at the time of strain replacement did not change in these two patients. A third case could not be validated by re-sequencing, but was supported as superinfection by an intermediate sequence with high degenerate base pair count within the time frame of strain switching. Drug resistance increased in this single patient.

Conclusions

Routine genotyping data are informative for the detection of HIV superinfection; however, most cases of non-monophyletic clustering in patient phylogenies arise from sample or sequence mix-up rather than from superinfection, which emphasizes the importance of validation. Non-transient superinfection was rare in our mainly treatment experienced cohort, and we found a single case of possible transmitted drug resistance by this route. We therefore conclude that in our large cohort, superinfection with drug resistant HIV did not compromise the efficiency of antiretroviral treatment.

Modeling sequence evolution in HIV-1 infection with recombination

Previously we proposed two simplified models of early HIV-1 evolution. Both showed that under a model of neutral evolution and exponential growth, the mean Hamming distance (HD) between genetic sequences grows linearly with time. In this paper we describe a more realistic continuous-time, age-dependent mathematical model of infection and viral replication, and show through simulations that even in this more complex description, the mean Hamming distance grows linearly with time. This remains unchanged when we introduce recombination, though the confidence intervals of the mean HD obtained ignoring recombination are overly conservative.

Accelerated crossing of fitness valleys through division of labor and cheating in asexual populations

Complex traits can require the accumulation of multiple mutations that are individually deleterious. Their evolution requires a fitness valley to be crossed, which can take relatively long time spans. A new evolutionary mechanism is described that accelerates the emergence of complex phenotypes, based on a “division of labor” game and the occurrence of cheaters. If each intermediate mutation leads to a product that can be shared with others, the complex type can arise relatively quickly as an emergent property among cooperating individuals, without any given individual having to accumulate all mutations. Moreover, the emergence of cheaters that destroy cooperative interactions can lead to the emergence of individuals that have accumulated all necessary mutations on a time scale that is significantly faster than observed in the absence of cooperation and cheating. Application of this mechanism to somatic and microbial evolution is discussed, including evolutionary processes in tumors, biofilms, and viral infections.

Estimating the fitness cost of escape from HLA presentation in HIV-1 protease and reverse transcriptase

Human immunodeficiency virus (HIV-1) is, like most pathogens, under selective pressure to escape the immune system of its host. In particular, HIV-1 can avoid recognition by cytotoxic T lymphocytes (CTLs) by altering the binding affinity of viral peptides to human leukocyte antigen (HLA) molecules, the role of which is to present those peptides to the immune system. It is generally assumed that HLA escape mutations carry a replicative fitness cost, but these costs have not been quantified. In this study, we assess the replicative cost of mutations which are likely to escape presentation by HLA molecules in the region of HIV-1 protease and reverse transcriptase. Specifically, we combine computational approaches for prediction of in vitro replicative fitness and peptide binding affinity to HLA molecules. We find that mutations which impair binding to HLA-A molecules tend to have lower in vitro replicative fitness than mutations which do not impair binding to HLA-A molecules, suggesting that HLA-A escape mutations carry higher fitness costs than non-escape mutations. We argue that the association between fitness and HLA-A binding impairment is probably due to an intrinsic cost of escape from HLA-A molecules, and these costs are particularly strong for HLA-A alleles associated with efficient virus control. Counter-intuitively, we do not observe a significant effect in the case of HLA-B, but, as discussed, this does not argue against the relevance of HLA-B in virus control. Overall, this article points to the intriguing possibility that HLA-A molecules preferentially target more conserved regions of HIV-1, emphasizing the importance of HLA-A genes in the evolution of HIV-1 and RNA viruses in general.

Our immune system can recognize and kill virus-infected cells by distinguishing between self and virus-derived protein fragments, called peptides, displayed on the surface of each cell. One requirement for a successful recognition is that those peptides bind to the human leukocyte antigen (HLA) class I molecules, which present them to the immune system. As a counter-strategy, human immunodeficiency virus type 1 (HIV-1) can acquire mutations that prevent this binding, thereby helping the virus to escape the surveillance of T-lymphocytes. It is likely that the virus pays a replicative cost for such escape mutations, but the magnitude of this cost has remained elusive. Here, we quantified this fitness cost in HIV-1 protease and reverse transcriptase by combining two computational systems biology approaches: one for prediction of in vitro replicative fitness, and one for the prediction of the efficiency of peptide binding to HLA. We found that in viral proteins targeted by HLA-A molecules, mutations which disrupt binding to those molecules carry a lower replicative fitness than mutations which do not have such an effect. We argue that these results are consistent with the hypothesis that our immune systems might have evolved to target genetic regions of RNA viruses which are costly for the pathogen to alter.

Evolutionary dynamics of HIV at multiple spatial and temporal scales

Infectious diseases remain a formidable challenge to human health, and understanding pathogen evolution is crucial to designing effective therapeutics and control strategies. Here, we review important evolutionary aspects of HIV infection, highlighting the concept of selection at multiple spatial and temporal scales. At the smallest scale, a single cell may be infected by multiple virions competing for intracellular resources. Recombination and phenotypic mixing introduce novel evolutionary dynamics. As the virus spreads between cells in an infected individual, it continually evolves to circumvent the immune system. We discuss evolutionary mechanisms of HIV pathogenesis and progression to AIDS. Viral spread throughout the human population can lead to changes in virulence and the transmission of immune-evading variation. HIV emerged as a human pathogen due to selection occurring between different species, adapting from related viruses of primates. HIV also evolves resistance to antiretroviral drugs within a single infected host, and we explore the possibility for the spread of these strains between hosts, leading to a drug-resistant epidemic. We investigate the role of latency, drug-protected compartments, and direct cell-to-cell transmission on viral evolution. The introduction of an HIV vaccine may select for viral variants that escape vaccine control, both within an individual and throughout the population. Due to the strong selective pressure exerted by HIV-induced morbidity and mortality in many parts of the world, the human population itself may be co-evolving in response to the HIV pandemic. Throughout the paper, we focus on trade-offs between costs and benefits that constrain viral evolution and accentuate how selection pressures differ at different levels of selection.

Dominant role of the 5' TAR bulge in dimerization of HIV-1 genomic RNA, but no evidence of TAR-TAR kissing during in vivo virus assembly

The 5' untranslated region of HIV-1 genomic RNA (gRNA) contains two stem-loop structures that appear to be equally important for gRNA dimerization: the 57-nucleotide 5' TAR, at the very 5' end, and the 35-nucleotide SL1 (nucleotides 243-277). SL1 is well-known for containing the dimerization initiation site (DIS) in its apical loop. The DIS is a six-nucleotide palindrome. Here, we investigated the mechanism of TAR-directed gRNA dimerization. We found that the trinucleotide bulge (UCU24) of the 5' TAR has dominant impacts on both formation of HIV-1 RNA dimers and maturation of the formed dimers. The ΔUCU trinucleotide deletion strongly inhibited the first process and blocked the other, thus impairing gRNA dimerization as severely as deletion of the entire 5' TAR, and more severely than deletion of the DIS, inactivation of the viral protease, or most severe mutations in the nucleocapsid protein. The apical loop of TAR contains a 10-nucleotide palindrome that has been postulated to stimulate gRNA dimerization by a TAR-TAR kissing mechanism analogous to the one used by SL1 to stimulate dimerization. Using mutations that strongly destabilize formation of the TAR palindrome duplex, as well as compensatory mutations that restore duplex formation to a wild-type-like level, we found no evidence of TAR-TAR kissing, even though mutations nullifying the kissing potential of the TAR palindrome could impair dimerization by a mechanism other than hindering of SL1. However, nullifying the kissing potential of TAR had much less severe effects than ΔUCU. By not uncovering a dimerization mechanism intrinsic to TAR, our data suggest that TAR mutations exert their effect 3' of TAR, yet not on SL1, because TAR and SL1 mutations have synergistic effects on gRNA dimerization.