A spatio-temporal assessment of simian/human immunodeficiency virus (SHIV) evolution reveals a highly dynamic process within the host

Reversions mask the contribution of adaptive evolution in microbiomes

When examining bacterial genomes for evidence of past selection, the results obtained depend heavily on the mutational distance between chosen genomes. Even within a bacterial species, genomes separated by larger mutational distances exhibit stronger evidence of purifying selection as assessed byd N / d S , the normalized ratio of nonsynonymous to synonymous mutations. Here, we show that the classical interpretation of this scale-dependence, weak purifying selection, leads to problematic mutation accumulation when applied to available gut microbiome data. We propose an alternative, adaptive reversion model with exactly opposite implications for dynamical intuition and applications ofd N / d S . Reversions that occur and sweep within-host populations are nearly guaranteed in microbiomes due to large population sizes, short generation times, and variable environments. Using analytical and simulation approaches, we show that adaptive reversion can explain thed N / d S decay given only dozens of locally-fluctuating selective pressures, which is realistic in the context of Bacteroides genomes. The success of the adaptive reversion model argues for interpreting low values ofd N / d S obtained from long-time scales with caution, as they may emerge even when adaptive sweeps are frequent. Our work thus inverts the interpretation of an old observation in bacterial evolution, illustrates the potential of mutational reversions to shape genomic landscapes over time, and highlights the importance of studying bacterial genomic evolution on short time scales.

Elevated HIV Viral Load is Associated with Higher Recombination Rate In Vivo

HIV's exceptionally high recombination rate drives its intrahost diversification, enabling immune escape and multidrug resistance within people living with HIV. While we know that HIV's recombination rate varies by genomic position, we have little understanding of how recombination varies throughout infection or between individuals as a function of the rate of cellular coinfection. We hypothesize that denser intrahost populations may have higher rates of coinfection and therefore recombination. To test this hypothesis, we develop a new approach (recombination analysis via time series linkage decay or RATS-LD) to quantify recombination using autocorrelation of linkage between mutations across time points. We validate RATS-LD on simulated data under short read sequencing conditions and then apply it to longitudinal, high-throughput intrahost viral sequencing data, stratifying populations by viral load (a proxy for density). Among sampled viral populations with the lowest viral loads (<26,800 copies/mL), we estimate a recombination rate of 1.5×10-5 events/bp/generation (95% CI: 7×10-6 to 2.9×10-5), similar to existing estimates. However, among samples with the highest viral loads (>82,000 copies/mL), our median estimate is approximately 6 times higher. In addition to co-varying across individuals, we also find that recombination rate and viral load are associated within single individuals across different time points. Our findings suggest that rather than acting as a constant, uniform force, recombination can vary dynamically and drastically across intrahost viral populations and within them over time. More broadly, we hypothesize that this phenomenon may affect other facultatively asexual populations where spatial co-localization varies.

Elevated HIV viral load is associated with higher recombination rate in vivo

HIV's exceptionally high recombination rate drives its intra-host diversification, enabling immune escape and multi-drug resistance within people living with HIV. While we know that HIV's recombination rate varies by genomic position, we have little understanding of how recombination varies throughout infection or between individuals as a function of the rate of cellular coinfection. We hypothesize that denser intra-host populations may have higher rates of coinfection and therefore recombination. To test this hypothesis, we develop a new approach (Recombination Analysis via Time Series Linkage Decay, or RATS-LD) to quantify recombination using autocorrelation of linkage between mutations across time points. We validate RATS-LD on simulated data under short read sequencing conditions and then apply it to longitudinal, high-throughput intra-host viral sequencing data, stratifying populations by viral load (a proxy for density). Among sampled viral populations with the lowest viral loads (< 26,800 copies/mL), we estimate a recombination rate of 1.5 × 10-5 events/bp/generation (95% CI: 7 × 10-6 - 2.9 × 10-5), similar to existing estimates. However, among samples with the highest viral loads (> 82,000 copies/mL), our median estimate is approximately 6 times higher. In addition to co-varying across individuals, we also find that recombination rate and viral load are associated within single individuals across different time points. Our findings suggest that rather than acting as a constant, uniform force, recombination can vary dynamically and drastically across intra-host viral populations and within them over time. More broadly, we hypothesize that this phenomenon may affect other facultatively asexual populations where spatial co-localization varies.

Design of an optimal combination therapy with broadly neutralizing antibodies to suppress HIV-1

Infusion of broadly neutralizing antibodies (bNAbs) has shown promise as an alternative to anti-retroviral therapy against HIV. A key challenge is to suppress viral escape, which is more effectively achieved with a combination of bNAbs. Here, we propose a computational approach to predict the efficacy of a bNAb therapy based on the population genetics of HIV escape, which we parametrize using high-throughput HIV sequence data from bNAb-naive patients. By quantifying the mutational target size and the fitness cost of HIV-1 escape from bNAbs, we predict the distribution of rebound times in three clinical trials. We show that a cocktail of three bNAbs is necessary to effectively suppress viral escape, and predict the optimal composition of such bNAb cocktail. Our results offer a rational therapy design for HIV, and show how genetic data can be used to predict treatment outcomes and design new approaches to pathogenic control.

HIV drug resistance in various body compartments

PURPOSE OF REVIEW

HIV drug resistance testing using blood plasma or dried blood spots forms part of international guidelines. However, as the clinical utility of assessing drug resistance in other body compartments is less well established, we review this for blood cells and samples from other body compartments.

RECENT EVIDENCE

Although clinical benefit is not clear, drug resistance testing in blood cells is often performed when patients with suppressed plasma viral loads require a treatment substitution. In patients with HIV neurocognitive disease, cerebral spinal fluid (CSF) drug resistance is rarely discordant with plasma but has nevertheless been used to guide antiretroviral drug substitutions. Cases with HIV drug resistance in genital fluids have been documented but this does not appear to indicate transmission risk when blood plasma viral loads are suppressed.

SUMMARY

Drug-resistant variants, which may be selected in tissues under conditions of variable adherence and drug penetration, appear to disseminate quickly, and become detectable in blood. This may explain why drug resistance discordance between plasma and these compartments is rarely found. Partial compartmentalization of HIV populations is well established for the CSF and the genital tract but other than blood plasma, evidence is lacking to support drug resistance testing in body compartments.

Repeated origins, widespread gene flow, and allelic interactions of target-site herbicide resistance mutations

Causal mutations and their frequency in agricultural fields are well-characterized for herbicide resistance. However, we still lack understanding of their evolutionary history: the extent of parallelism in the origins of target-site resistance (TSR), how long these mutations persist, how quickly they spread, and allelic interactions that mediate their selective advantage. We addressed these questions with genomic data from 19 agricultural populations of common waterhemp (Amaranthus tuberculatus), which we show to have undergone a massive expansion over the past century, with a contemporary effective population size estimate of 8 x 10⁷. We found variation at seven characterized TSR loci, two of which had multiple amino acid substitutions, and three of which were common. These three common resistance variants show extreme parallelism in their mutational origins, with gene flow having shaped their distribution across the landscape. Allele age estimates supported a strong role of adaptation from de novo mutations, with a median age of 30 suggesting that most resistance alleles arose soon after the onset of herbicide use. However, resistant lineages varied in both their age and evidence for selection over two different timescales, implying considerable heterogeneity in the forces that govern their persistence. Two such forces are intra- and inter-locus allelic interactions; we report a signal of extended haplotype competition between two common TSR alleles, and extreme linkage with genome-wide alleles with known functions in resistance adaptation. Together, this work reveals a remarkable example of spatial parallel evolution in a metapopulation, with important implications for the management of herbicide resistance.

SIV Evolutionary Dynamics in Cynomolgus Macaques during SIV-Mycobacterium tuberculosis Co-Infection

Co-infection with Mycobacterium tuberculosis (Mtb) and human immunodeficiency virus (HIV) is a worldwide public health concern, leading to worse clinical outcomes caused by both pathogens. We used a non-human primate model of simian immunodeficiency virus (SIV)-Mtb co-infection, in which latent Mtb infection was established prior to SIV_mac251 infection. The evolutionary dynamics of SIV env was evaluated from samples in plasma, lymph nodes, and lungs (including granulomas) of SIV-Mtb co-infected and SIV only control animals. While the diversity of the challenge virus was low and overall viral diversity remained relatively low over 6–9 weeks, changes in viral diversity and divergence were observed, including evidence for tissue compartmentalization. Overall, viral diversity was highest in SIV-Mtb animals that did not develop clinical Mtb reactivation compared to animals with Mtb reactivation. Among lung granulomas, viral diversity was positively correlated with the frequency of CD4+ T cells and negatively correlated with the frequency of CD8+ T cells. SIV diversity was highest in the thoracic lymph nodes compared to other sites, suggesting that lymphatic drainage from the lungs in co-infected animals provides an advantageous environment for SIV replication. This is the first assessment of SIV diversity across tissue compartments during SIV-Mtb co-infection after established Mtb latency.

Predator-Prey Dynamics of Intra-Host Simian Immunodeficiency Virus Evolution Within the Untreated Host

The dynamic nature of the SIV population during disease progression in the SIV/macaque model of AIDS and the factors responsible for its behavior have not been documented, largely owing to the lack of sufficient spatial and temporal sampling of both viral and host data from SIV-infected animals. In this study, we detail Bayesian coalescent inference of the changing collective intra-host viral effective population size (N_e ) from various tissues over the course of infection and its relationship with what we demonstrate is a continuously changing immune cell repertoire within the blood. Although the relative contribution of these factors varied among hosts and time points, the adaptive immune response best explained the overall periodic dynamic behavior of the effective virus population. Data exposing the nature of the relationship between the virus and immune cell populations revealed the plausibility of an eco-evolutionary mathematical model, which was able to mimic the large-scale oscillations in N_e through virus escape from relatively few, early immunodominant responses, followed by slower escape from several subdominant and weakened immune populations. The results of this study suggest that SIV diversity within the untreated host is governed by a predator-prey relationship, wherein differing phases of infection are the result of adaptation in response to varying immune responses. Previous investigations into viral population dynamics using sequence data have focused on single estimates of the effective viral population size (N_e ) or point estimates over sparse sampling data to provide insight into the precise impact of immune selection on virus adaptive behavior. Herein, we describe the use of the coalescent phylogenetic frame- work to estimate the relative changes in N_e over time in order to quantify the relationship with empirical data on the dynamic immune composition of the host. This relationship has allowed us to expand on earlier simulations to build a predator-prey model that explains the deterministic behavior of the virus over the course of disease progression. We show that sequential viral adaptation can occur in response to phases of varying immune pressure, providing a broader picture of the viral response throughout the entire course of progression to AIDS.

Single-genome sequencing reveals within-host evolution of human malaria parasites

Population genomics of bulk malaria infections is unable to examine intrahost evolution; therefore, most work has focused on the role of recombination in generating genetic variation. We used single-cell sequencing protocol for low-parasitaemia infections to generate 406 near-complete single Plasmodium vivax genomes from 11 patients sampled during sequential febrile episodes. Parasite genomes contain hundreds of de novo mutations, showing strong signatures of selection, which are enriched in the ApiAP2 family of transcription factors, known targets of adaptation. Comparing 315 P. falciparum single-cell genomes from 15 patients with our P. vivax data, we find broad complementary patterns of de novo mutation at the gene and pathway level, revealing the importance of within-host evolution during malaria infections.

Quo vadis? Central Rules of Pathogen and Disease Tropism

Understanding why certain people get sick and die while others recover or never become ill is a fundamental question in biomedical research. A key determinant of this process is pathogen and disease tropism: the locations that become infected (pathogen tropism), and the locations that become damaged (disease tropism). Identifying the factors that regulate tropism is essential to understand disease processes, but also to drive the development of new interventions. This review intersects research from across infectious diseases to define the central mediators of disease and pathogen tropism. This review also highlights methods of study, and translational implications. Overall, tropism is a central but under-appreciated aspect of infection pathogenesis which should be at the forefront when considering the development of new methods of intervention.

The clarifying role of time series data in the population genetics of HIV

HIV can evolve remarkably quickly in response to antiretroviral therapies and the immune system. This evolution stymies treatment effectiveness and prevents the development of an HIV vaccine. Consequently, there has been a great interest in using population genetics to disentangle the forces that govern the HIV adaptive landscape (selection, drift, mutation, and recombination). Traditional population genetics approaches look at the current state of genetic variation and infer the processes that can generate it. However, because HIV evolves rapidly, we can also sample populations repeatedly over time and watch evolution in action. In this paper, we demonstrate how time series data can bound evolutionary parameters in a way that complements and informs traditional population genetic approaches. Specifically, we focus on our recent paper (Feder et al., 2016, eLife), in which we show that, as improved HIV drugs have led to fewer patients failing therapy due to resistance evolution, less genetic diversity has been maintained following the fixation of drug resistance mutations. Because soft sweeps of multiple drug resistance mutations spreading simultaneously have been previously documented in response to the less effective HIV therapies used early in the epidemic, we interpret the maintenance of post-sweep diversity in response to poor therapies as further evidence of soft sweeps and therefore a high population mutation rate (θ) in these intra-patient HIV populations. Because improved drugs resulted in rarer resistance evolution accompanied by lower post-sweep diversity, we suggest that both observations can be explained by decreased population mutation rates and a resultant transition to hard selective sweeps. A recent paper (Harris et al., 2018, PLOS Genetics) proposed an alternative interpretation: Diversity maintenance following drug resistance evolution in response to poor therapies may have been driven by recombination during slow, hard selective sweeps of single mutations. Then, if better drugs have led to faster hard selective sweeps of resistance, recombination will have less time to rescue diversity during the sweep, recapitulating the decrease in post-sweep diversity as drugs have improved. In this paper, we use time series data to show that drug resistance evolution during ineffective treatment is very fast, providing new evidence that soft sweeps drove early HIV treatment failure.

Long-term experimental evolution of HIV-1 reveals effects of environment and mutational history

An often-returning question for not only HIV-1, but also other organisms, is how predictable evolutionary paths are. The environment, mutational history, and random processes can all impact the exact evolutionary paths, but to which extent these factors contribute to the evolutionary dynamics of a particular system is an open question. Especially in a virus like HIV-1, with a large mutation rate and large population sizes, evolution is expected to be highly predictable if the impact of environment and history is low, and evolution is not neutral. We investigated the effect of environment and mutational history by analyzing sequences from a long-term evolution experiment, in which HIV-1 was passaged on 2 different cell types in 8 independent evolutionary lines and 8 derived lines, 4 of which involved a switch of the environment. The experiments lasted for 240–300 passages, corresponding to approximately 400–600 generations or almost 3 years. The sequences show signs of extensive parallel evolution—the majority of mutations that are shared between independent lines appear in both cell types, but we also find that both environment and mutational history significantly impact the evolutionary paths. We conclude that HIV-1 evolution is robust to small changes in the environment, similar to a transmission event in the absence of an immune response or drug pressure. We also find that the fitness landscape of HIV-1 is largely smooth, although we find some evidence for both positive and negative epistatic interactions between mutations.

Analysis of the longest evolutionary experiment with HIV-1 to-date reveals continuous viral adaptation over several years. The authors quantify the environment-specific mutations that arise and determine the fraction of mutations that co-occur with significantly different frequencies than expected by chance.

Robust Estimation of Recent Effective Population Size from Number of Independent Origins in Soft Sweeps

Estimating recent effective population size is of great importance in characterizing and predicting the evolution of natural populations. Methods based on nucleotide diversity may underestimate current day effective population sizes due to historical bottlenecks, whereas methods that reconstruct demographic history typically only detect long-term variations. However, soft selective sweeps, which leave a fingerprint of mutational history by recurrent mutations on independent haplotype backgrounds, holds promise of an estimate more representative of recent population history. Here, we present a simple and robust method of estimation based only on knowledge of the number of independent recurrent origins and the current frequency of the beneficial allele in a population sample, independent of the strength of selection and age of the mutation. Using a forward-time theoretical framework, we show the mean number of origins is a function of θ=2Nμ and current allele frequency, through a simple equation, and the distribution is approximately Poisson. This estimate is robust to whether mutants preexisted before selection arose and is equally accurate for diploid populations with incomplete dominance. For fast (e.g., seasonal) demographic changes compared with time scale for fixation of the mutant allele, and for moderate peak-to-trough ratios, we show our constant population size estimate can be used to bound the maximum and minimum population size. Applied to the Vgsc gene of Anopheles gambiae, we estimate an effective population size of roughly 6×107, and including seasonal demographic oscillations, a minimum effective population size >3×107, and a maximum <6×109, suggesting a mean ∼109.

Evolutionary Dynamics in Structured Populations Under Strong Population Genetic Forces

In the long-term neutral equilibrium, high rates of migration between subpopulations result in little population differentiation. However, in the short-term, even very abundant migration may not be enough for subpopulations to equilibrate immediately. In this study, we investigate dynamical patterns of short-term population differentiation in adapting populations via stochastic and analytical modeling through time. We characterize a regime in which selection and migration interact to create non-monotonic patterns of population differentiation over time when migration is weaker than selection, but stronger than drift. We demonstrate how these patterns can be leveraged to estimate high migration rates using approximate Bayesian computation. We apply this approach to estimate fast migration in a rapidly adapting intra-host Simian-HIV population sampled from different anatomical locations. We find differences in estimated migration rates between different compartments, even though all are above [Formula: see text] = 1. This work demonstrates how studying demographic processes on the timescale of selective sweeps illuminates processes too fast to leave signatures on neutral timescales.

Virus dynamics and phyloanatomy: Merging population dynamic and phylogenetic approaches

In evolutionary biology and epidemiology, phylodynamic methods are widely used to infer population biological characteristics, such as the rates of replication, death, migration, or, in the epidemiological context, pathogen spread. More recently, these methods have been used to elucidate the dynamics of viruses within their hosts. Especially the application of phylogeographic approaches has the potential to shed light on anatomical colonization pathways and the exchange of viruses between distinct anatomical compartments. We and others have termed this phyloanatomy. Here, we review the promise and challenges of phyloanatomy, and compare them to more classical virus dynamics and population genetic approaches. We argue that the extremely strong selection pressures that exist within the host may represent the main obstacle to reliable phyloanatomic analysis.

Intrahost Selection Pressure Drives Equine Arteritis Virus Evolution during Persistent Infection in the Stallion Reproductive Tract

Equine arteritis virus (EAV) is the causative agent of equine viral arteritis (EVA), a reproductive and respiratory disease of horses. Following natural infection, 10 to 70% of infected stallions can become carriers of EAV and continue to shed virus in the semen. In this study, sequential viruses isolated from nasal secretions, buffy coat cells, and semen of seven experimentally infected and two naturally infected EAV carrier stallions were deep sequenced to elucidate the intrahost microevolutionary process after a single transmission event. Analysis of variants from nasal secretions and buffy coat cells lacked extensive positive selection; however, characteristics of the mutant spectra were different in the two sample types. In contrast, the initial semen virus populations during acute infection have undergone a selective bottleneck, as reflected by the reduction in population size and diversifying selection at multiple sites in the viral genome. Furthermore, during persistent infection, extensive genome-wide purifying selection shaped variant diversity in the stallion reproductive tract. Overall, the nonstochastic nature of EAV evolution during persistent infection was driven by active intrahost selection pressure. Among the open reading frames within the viral genome, ORF3, ORF5, and the nsp2-coding region of ORF1a accumulated the majority of nucleotide substitutions during persistence, with ORF3 and ORF5 having the highest intrahost evolutionary rates. The findings presented here provide a novel insight into the evolutionary mechanisms of EAV and identified critical regions of the viral genome likely associated with the establishment and maintenance of persistent infection in the stallion reproductive tract.IMPORTANCE EAV can persist in the reproductive tract of infected stallions, and consequently, long-term carrier stallions constitute its sole natural reservoir. Previous studies demonstrated that the ampullae of the vas deferens are the primary site of viral persistence in the stallion reproductive tract and the persistence is associated with a significant inflammatory response that is unable to clear the infection. This is the first study that describes EAV full-length genomic evolution during acute and long-term persistent infection in the stallion reproductive tract using next-generation sequencing and contemporary sequence analysis techniques. The data provide novel insight into the intrahost evolution of EAV during acute and persistent infection and demonstrate that persistent infection is characterized by extensive genome-wide purifying selection and a nonstochastic evolutionary pattern mediated by intrahost selective pressure, with important nucleotide substitutions occurring in ORF1a (region encoding nsp2), ORF3, and ORF5.

Genetic diversity in the env V1-V2 region of proviral quasispecies from long-term controller MHC-typed cynomolgus macaques infected with SHIVSF162P4cy

Intra-host evolution of human immunodeficiency virus (HIV) and simian immunodeficiency virus (SIV) has been shown by viral RNA analysis in subjects who naturally suppress plasma viremia to low levels, known as controllers. However, little is known about the variability of proviral DNA and the inter-relationships among contained systemic viremia, rate of reservoir reseeding and specific major histocompatibility complex (MHC) genotypes, in controllers. Here, we analysed the proviral DNA quasispecies of the env V1-V2 region, in PBMCs and in anatomical compartments of 13 long-term controller monkeys after 3.2 years of infection with simian/human immunodeficiency virus (SHIV)SF162P4cy. A considerable variation in the genetic diversity of proviral quasispecies was present among animals. Seven monkeys exhibited env V1-V2 proviral populations composed of both clusters of identical ancestral sequences and new variants, whereas the other six monkeys displayed relatively high env V1-V2 genetic diversity with a large proportion of diverse novel sequences. Our results demonstrate that in SHIVSF162P4cy-infected monkeys there exists a disparate pattern of intra-host viral diversity and that reseeding of the proviral reservoir occurs in some animals. Moreover, even though no particular association has been observed between MHC haplotypes and the long-term control of infection, a remarkably similar pattern of intra-host viral diversity and divergence was found within animals carrying the M3 haplotype. This suggests that in animals bearing the same MHC haplotype and infected with the same virus, viral diversity follows a similar pattern with similar outcomes and control of infection.

Mechanisms and Concepts in RNA Virus Population Dynamics and Evolution

RNA viruses are unique in their evolutionary capacity, exhibiting high mutation rates and frequent recombination. They rapidly adapt to environmental changes, such as shifts in immune pressure or pharmacological challenge. The evolution of RNA viruses has been brought into new focus with the recent developments of genetic and experimental tools to explore and manipulate the evolutionary dynamics of viral populations. These studies have uncovered new mechanisms that enable viruses to overcome evolutionary challenges in the environment and have emphasized the intimate relationship of viral populations with evolution. Here, we review some of the emerging viral and host mechanisms that underlie the evolution of RNA viruses. We also discuss new studies that demonstrate that the relationship between evolutionary dynamics and virus biology spans many spatial and temporal scales, affecting transmission dynamics within and between hosts as well as pathogenesis.

On the effective depth of viral sequence data

Genome sequence data are of great value in describing evolutionary processes in viral populations. However, in such studies, the extent to which data accurately describes the viral population is a matter of importance. Multiple factors may influence the accuracy of a dataset, including the quantity and nature of the sample collected, and the subsequent steps in viral processing. To investigate this phenomenon, we sequenced replica datasets spanning a range of viruses, and in which the point at which samples were split was different in each case, from a dataset in which independent samples were collected from a single patient to another in which all processing steps up to sequencing were applied to a single sample before splitting the sample and sequencing each replicate. We conclude that neither a high read depth nor a high template number in a sample guarantee the precision of a dataset. Measures of consistency calculated from within a single biological sample may also be insufficient; distortion of the composition of a population by the experimental procedure or genuine within-host diversity between samples may each affect the results. Where it is possible, data from replicate samples should be collected to validate the consistency of short-read sequence data.

Allele Age Under Non-Classical Assumptions is Clarified by an Exact Computational Markov Chain Approach

Determination of the age of an allele based on its population frequency is a well-studied problem in population genetics, for which a variety of approximations have been proposed. We present a new result that, surprisingly, allows the expectation and variance of allele age to be computed exactly (within machine precision) for any finite absorbing Markov chain model in a matter of seconds. This approach makes none of the classical assumptions (e.g., weak selection, reversibility, infinite sites), exploits modern sparse linear algebra techniques, integrates over all sample paths, and is rapidly computable for Wright-Fisher populations up to N_e = 100,000. With this approach, we study the joint effect of recurrent mutation, dominance, and selection, and demonstrate new examples of “selective strolls” where the classical symmetry of allele age with respect to selection is violated by weakly selected alleles that are older than neutral alleles at the same frequency. We also show evidence for a strong age imbalance, where rare deleterious alleles are expected to be substantially older than advantageous alleles observed at the same frequency when population-scaled mutation rates are large. These results highlight the under-appreciated utility of computational methods for the direct analysis of Markov chain models in population genetics.