Assessing the effects of human mixing patterns on human immunodeficiency virus-1 interhost phylogenetics through social network simulation

Geneticists seeking to understand HIV-1 evolution among human hosts generally assume that hosts represent a panmictic population. Social science research demonstrates that the network patterns over which HIV-1 spreads are highly nonrandom, but the effect of these patterns on the genetic diversity of HIV-1 and other sexually transmitted pathogens has yet to be thoroughly examined. In addition, interhost phylogenetic models rarely account explicitly for genetic diversity arising from intrahost dynamics. This study outlines a graph-theoretic framework (exponential random graph modeling, ERGM) for the estimation, inference, and simulation of dynamic partnership networks. This approach is used to simulate HIV-1 transmission and evolution under eight mixing patterns resembling those observed in empirical human populations, while simultaneously incorporating intrahost viral diversity. Models of parametric growth fit panmictic populations well, yielding estimates of total viral effective population on the order of the product of infected host size and intrahost effective viral population size. Populations exhibiting patterns of nonrandom mixing differ more widely in estimates of effective population size they yield, however, and reconstructions of population dynamics can exhibit severe errors if panmixis is assumed. I discuss implications for HIV-1 phylogenetics and the potential for ERGM to provide a general framework for addressing these issues.     

Genome-wide patterns of intrahuman dengue virus diversity reveal associations with viral phylogenetic clade and interhost diversity

Analogous to observations in RNA viruses such as human immunodeficiency virus, genetic variation associated with intrahost dengue virus (DENV) populations has been postulated to influence viral fitness and disease pathogenesis. Previous attempts to investigate intrahost genetic variation in DENV characterized only a few viral genes or a limited number of full-length genomes. We developed a whole-genome amplification approach coupled with deep sequencing to capture intrahost diversity across the entire coding region of DENV-2. Using this approach, we sequenced DENV-2 genomes from the serum of 22 Nicaraguan individuals with secondary DENV infection and captured ～75% of the DENV genome in each sample (range, 40 to 98%). We identified and quantified variants using a highly sensitive and specific method and determined that the extent of diversity was considerably lower than previous estimates. Significant differences in intrahost diversity were detected between genes and also between antigenically distinct domains of the Envelope gene. Interestingly, a strong association was discerned between the extent of intrahost diversity in a few genes and viral clade identity. Additionally, the abundance of viral variants within a host, as well as the impact of viral mutations on amino acid encoding and predicted protein function, determined whether intrahost variants were observed at the interhost level in circulating Nicaraguan DENV-2 populations, strongly suggestive of purifying selection across transmission events. Our data illustrate the value of high-coverage genome-wide analysis of intrahost diversity for high-resolution mapping of the relationship between intrahost diversity and clinical, epidemiological, and virological parameters of viral infection.     

Baseline CD4+ T Cell Counts Correlates with HIV-1 Synonymous Rate in HLA-B*5701 Subjects with Different Risk of Disease Progression

HLA-B*5701 is the host factor most strongly associated with slow HIV-1 disease progression, although risk of progression may vary among patients carrying this allele. The interplay between HIV-1 evolutionary rate variation and risk of progression to AIDS in HLA-B*5701 subjects was studied using longitudinal viral sequences from high-risk progressors (HRPs) and low-risk progressors (LRPs). Posterior distributions of HIV-1 genealogies assuming a Bayesian relaxed molecular clock were used to estimate the absolute rates of nonsynonymous and synonymous substitutions for different set of branches. Rates of viral evolution, as well as in vitro viral replication capacity assessed using a novel phenotypic assay, were correlated with various clinical parameters. HIV-1 synonymous substitution rates were significantly lower in LRPs than HRPs, especially for sets of internal branches. The viral population infecting LRPs was also characterized by a slower increase in synonymous divergence over time. This pattern did not correlate to differences in viral fitness, as measured by in vitro replication capacity, nor could be explained by differences among subjects in T cell activation or selection pressure. Interestingly, a significant inverse correlation was found between baseline CD4⁺ T cell counts and mean HIV-1 synonymous rate (which is proportional to the viral replication rate) along branches representing viral lineages successfully propagating through time up to the last sampled time point. The observed lower replication rate in HLA-B*5701 subjects with higher baseline CD4⁺ T cell counts provides a potential model to explain differences in risk of disease progression among individuals carrying this allele.     

Generalized selection to overcome innate immunity selects for host breadth in an RNA virus

Virus‐host coevolution has selected for generalized host defense against viruses, exemplified by interferon production/signaling and other innate immune function in eukaryotes such as humans. Although cell‐surface binding primarily limits virus infection success, generalized adaptation to counteract innate immunity across disparate hosts may contribute to RNA virus emergence potential. We examined this idea using vesicular stomatitis virus (VSV) populations previously evolved on strictly immune‐deficient (HeLa) cells, strictly immune competent (MDCK) cells, or on alternating deficient/competent cells. By measuring viral fitness in unselected human cancer cells of differing innate immunity, we confirmed that HeLa‐adapted populations were specialized for innate immune‐deficient hosts, whereas MDCK‐adapted populations were relatively more generalized for fitness on hosts of differing innate immune capacity and of different species origin. We also confirmed that HeLa‐evolved populations maintained fitness in immune‐deficient nonhuman primate cells. These results suggest that innate immunity is more prominent than host species in determining viral fitness at the host‐cell level. Finally, our prediction was inexact that selection on alternating deficient/competent hosts should produce innate viral generalists. Rather, fitness differences among alternating host‐evolved VSV populations indicated variable capacities to evade innate immunity. Our results suggest that the evolutionary history of innate immune selection can affect whether RNA viruses evolve greater host‐breadth.     

Differential trends in the codon usage patterns in HIV-1 genes

Host-pathogen interactions underlie one of the most complex evolutionary phenomena resulting in continual adaptive genetic changes, where pathogens exploit the host's molecular resources for growth and survival, while hosts try to eliminate the pathogen. Deciphering the molecular basis of host-pathogen interactions is useful in understanding the factors governing pathogen evolution and disease propagation. In host-pathogen context, a balance between mutation, selection, and genetic drift is known to maintain codon bias in both organisms. Studies revealing determinants of the bias and its dynamics are central to the understanding of host-pathogen evolution. We considered the Human Immunodeficiency Virus (HIV) type 1 and its human host to search for evolutionary signatures in the viral genome. Positive selection is known to dominate intra-host evolution of HIV-1, whereas high genetic variability underlies the belief that neutral processes drive inter-host differences. In this study, we analyze the codon usage patterns of HIV-1 genomes across all subtypes and clades sequenced over a period of 23 years. We show presence of unique temporal correlations in the codon bias of three HIV-1 genes illustrating differential adaptation of the HIV-1 genes towards the host preferred codons. Our results point towards gene-specific translational selection to be an important force driving the evolution of HIV-1 at the population level.     

Intercompartmental recombination of HIV-1 contributes to env intrahost diversity and modulates viral tropism and sensitivity to entry inhibitors

HIV-1 circulates within an infected host as a genetically heterogeneous viral population. Viral intrahost diversity is shaped by substitutional evolution and recombination. Although many studies have speculated that recombination could have a significant impact on viral phenotype, this has never been definitively demonstrated. We report here phylogenetic and subsequent phenotypic analyses of envelope genes obtained from HIV-1 populations present in different anatomical compartments. Assessment of env compartmentalization from immunologically discrete tissues was assessed utilizing a single genome amplification approach, minimizing in vitro-generated artifacts. Genetic compartmentalization of variants was frequently observed. In addition, multiple incidences of intercompartment recombination, presumably facilitated by low-level migration of virus or infected cells between different anatomic sites and coinfection of susceptible cells by genetically divergent strains, were identified. These analyses demonstrate that intercompartment recombination is a fundamental evolutionary mechanism that helps to shape HIV-1 env intrahost diversity in natural infection. Analysis of the phenotypic consequences of these recombination events showed that genetic compartmentalization often correlates with phenotypic compartmentalization and that intercompartment recombination results in phenotype modulation. This represents definitive proof that recombination can generate novel combinations of phenotypic traits which differ subtly from those of parental strains, an important phenomenon that may have an impact on antiviral therapy and contribute to HIV-1 persistence in vivo.     

Distribution of HIV-1 resistance-conferring polymorphic allelesSDF-1-3′A, CCR2-64I andCCR5-Δ32 in diverse populations of Andhra Pradesh, South India

Polymorphic allelic variants of chemokine receptors CCR2 and CCR5, as well as of stromal-derived factor-1 SDF-1, the ligand for the chemokine receptor CXCR4, are known to have protective effects against HIV-1 infection and to be involved with delay in disease progression. We have studied the DNA polymorphisms at the loci that encode these proteins in 525 healthy individuals without any history of HIV-1 infection from 11 diverse populations of Andhra Pradesh, South India. The two protective allelesSDF-1-3′A andCCR2-64I at the SDF-1 and CCR2 loci, respectively, are present in all populations studied, although their frequencies differ considerably across populations (from 17% to 35% for theSDF-1-3′A allele, and from 3% to 17% forCCR2-64I). In contrast theCCR5-δ32 allele is observed only in three populations (Yamani, Pathan and Kamma), all in low frequencies (i.e. 1% to 3%). The mean number of mutant alleles (for the three loci together) carried by each individual varies from 0.475 (in Vizag Brahmins) to 0.959 (in Bohra Muslims). The estimated relative hazard values for the populations, computed from the three-locus genotype data, are comparable to those from Africa and Southeast Asia, where AIDS is known to be widespread.     

Dendritic cells preferentially transfer CXCR4-using human immunodeficiency virus type 1 variants to CD4+ T lymphocytes in trans

Human immunodeficiency virus type 1 (HIV-1) preferentially utilizes the CCR5 coreceptor for target cell entry in the acute phase of infection, while later in disease progression the virus switches to the CXCR4 coreceptor in approximately 50% of patients. In response to HIV-1 the adaptive immune response is triggered, and antibody (Ab) production is elicited to block HIV-1 entry. We recently determined that dendritic cells (DCs) can efficiently capture Ab-neutralized HIV-1, restore infectivity, and transmit infectious virus to target cells. Here, we tested the effect of Abs on trans transmission of CCR5 or CXCR4 HIV-1 variants. We observed that transmission of HIV-1 by immature as well as mature DCs was significantly higher for CXCR4- than CCR5-tropic viral strains. Additionally, neutralizing Abs directed against either the gp41 or gp120 region of the envelope such as 2F5, 4E10, and V3-directed Abs inhibited transmission of CCR5-tropic HIV-1, whereas Ab-treated CXCR4-tropic virus demonstrated unaltered or increased transmission. To further study the effects of coreceptor usage we tested molecularly cloned HIV-1 variants with modifications in the envelope that were based on longitudinal gp120 V1 and V3 variable loop sequences from a patient progressing to AIDS. We observed that DCs preferentially facilitated infection of CD4⁺ T lymphocytes of viral strains with an envelope phenotype found late in disease. Taken together, our results illustrate that DCs transmit CXCR4-tropic HIV-1 much more efficiently than CCR5 strains; we hypothesize that this discrimination could contribute to the in vivo coreceptor switch after seroconversion and could be responsible for the increase in viral load.     

Temporal dynamics of SARS-CoV-2 mutation accumulation within and across infected hosts

Analysis of SARS-CoV-2 genetic diversity within infected hosts can provide insight into the generation and spread of new viral variants and may enable high resolution inference of transmission chains. However, little is known about temporal aspects of SARS-CoV-2 intrahost diversity and the extent to which shared diversity reflects convergent evolution as opposed to transmission linkage. Here we use high depth of coverage sequencing to identify within-host genetic variants in 325 specimens from hospitalized COVID-19 patients and infected employees at a single medical center. We validated our variant calling by sequencing defined RNA mixtures and identified viral load as a critical factor in variant identification. By leveraging clinical metadata, we found that intrahost diversity is low and does not vary by time from symptom onset. This suggests that variants will only rarely rise to appreciable frequency prior to transmission. Although there was generally little shared variation across the sequenced cohort, we identified intrahost variants shared across individuals who were unlikely to be related by transmission. These variants did not precede a rise in frequency in global consensus genomes, suggesting that intrahost variants may have limited utility for predicting future lineages. These results provide important context for sequence-based inference in SARS-CoV-2 evolution and epidemiology.     

Highly sensitive and specific detection of rare variants in mixed viral populations from massively parallel sequence data

Viruses diversify over time within hosts, often undercutting the effectiveness of host defenses and therapeutic interventions. To design successful vaccines and therapeutics, it is critical to better understand viral diversification, including comprehensively characterizing the genetic variants in viral intra-host populations and modeling changes from transmission through the course of infection. Massively parallel sequencing technologies can overcome the cost constraints of older sequencing methods and obtain the high sequence coverage needed to detect rare genetic variants (< 1%) within an infected host, and to assay variants without prior knowledge. Critical to interpreting deep sequence data sets is the ability to distinguish biological variants from process errors with high sensitivity and specificity. To address this challenge, we describe V-Phaser, an algorithm able to recognize rare biological variants in mixed populations. V-Phaser uses covariation (i.e. phasing) between observed variants to increase sensitivity and an expectation maximization algorithm that iteratively recalibrates base quality scores to increase specificity. Overall, V-Phaser achieved > 97% sensitivity and > 97% specificity on control read sets. On data derived from a patient after four years of HIV-1 infection, V-Phaser detected 2,015 variants across the -10 kb genome, including 603 rare variants (< 1% frequency) detected only using phase information. V-Phaser identified variants at frequencies down to 0.2%, comparable to the detection threshold of allele-specific PCR, a method that requires prior knowledge of the variants. The high sensitivity and specificity of V-Phaser enables identifying and tracking changes in low frequency variants in mixed populations such as RNA viruses.     

Human immunodeficiency virus type 1 genetic evolution in children with different rates of development of disease.

The rate of development of disease varies considerably among human immunodeficiency virus type 1 (HIV-1)-infected children. The reasons for these observed differences are not clearly understood but most probably depend on the dynamic interplay between the HIV-1 quasispecies virus population and the immune constraints imposed by the host. To study the relationship between disease progression and genetic diversity, we analyzed the evolution of viral sequences within six perinatally infected children by examining proviral sequences spanning the C2 through V5 regions of the viral envelope gene by PCR of blood samples obtained at sequential visits. PCR product DNAs from four sample time points per child were cloned, and 10 to 13 clones from each sample were sequenced. Greater genetic distances relative to the time of infection were found for children with low virion-associated RNA burdens and slow progression to disease relative to those found for children with high virion-associated RNA burdens and rapid progression to disease. The greater branch lengths observed in the phylogenetic reconstructions correlated with a higher accumulation rate of nonsynonymous base substitutions per potential nonsynonymous site, consistent with positive selection for change rather than a difference in replication kinetics. Viral sequences from children with slow progression to disease also showed a tendency to form clusters that associated with different sampling times. These progressive shifts in the viral population were not found in viral sequences from children with rapid progression to disease. Therefore, despite the HIV-1 quasispecies being a diverse, rapidly evolving, and competing population of genetic variants, different rates of genetic evolution could be found under different selective constraints. These data suggest that the evolutionary dynamics exhibited by the HIV-1 quasispecies virus populations are compatible with a Darwinian system evolving under the constraints of natural selection.     

Evolution of CCR5 use before and during coreceptor switching

The envelope gene (env) of human immunodeficiency virus type 1 (HIV-1) undergoes rapid divergence from the transmitted sequence and increasing diversification during the prolonged course of chronic infection in humans. In about half of infected individuals or more, env evolution leads to expansion of the use of entry coreceptor from CCR5 alone to CCR5 and CXCR4. The stochastic nature of this coreceptor switch is not well explained by host selective forces that should be relatively constant between infected individuals. Moreover, differences in the incidence of coreceptor switching among different HIV-1 subtypes suggest that properties of the evolving virus population drive the switch. We evaluated the functional properties of sequential env clones from a patient with evidence of coreceptor switching at 5.67 years of infection. We found an abrupt decline in the ability of viruses to use CCR5 for entry at this time, manifested by a 1- to 2-log increase in susceptibility to CCR5 inhibitors and a reduced ability to infect cell lines with low CCR5 expression. There was an abnormally rapid 5.4% divergence in env sequences from 4.10 to 5.76 years of infection, with the V3 and V4/V5 regions showing the greatest divergence and evidence of positive selection. These observations suggest that a decline in the fitness of R5 virus populations may be one driving force that permits the emergence of R5X4 variants.     

CCR2-64I allele is associated with the progression of AIDS in a Han Chinese population

It has been reported that the CCR5-Delta32, CCR2-64I and SDF1-3′A polymorphisms have protective effects against HIV-1 infection and can delay the progression of AIDS in European populations. The aim of this study was to investigate the associations of the three loci with HIV-1 infection and the progression of AIDS in the Han population of northern China. We recruited 78 HIV-1-seropositive individuals and 70 population-matched healthy controls from north China. PCR-RFLP was used to genotype these three polymorphisms in all samples. There were no significant differences in the frequencies of any allele between HIV-1-seropositive individuals and healthy controls. In cases, the associations between the three loci and CD4⁺ T-lymphocyte counts were not significant; however, there was a significant association between the genotype frequency of CCR2-64I and clinical category (P < 0.05). We have identified a potentially important role of CCR2-64I allele in AIDS progression in the northern Chinese population. Further studies are necessary to clarify the contribution of the CCR2-64I allele to the progression of AIDS for the relatively small sample size.     

Chemokine receptor CCR5 genotype influences the kinetics of human immunodeficiency virus type 1 infection in human PBL-SCID mice.

Individuals homozygous for a 32-bp deletion (delta 32) in the CCR5 gene encoding the coreceptor for macrophage-tropic human immunodeficiency virus type 1 (HIV-1) are resistant to virus infection, and heterozygous individuals show some slowing of disease progression. The impact of the CCR5 genotype on HIV-1 infection was assessed in vitro and in the human PBL-SCID (hu-PBL-SCID) model. Cells and hu-PBL-SCID mice from CCR5 delta 32/delta 32 donors were resistant to infection with macrophage-tropic HIV-1 and showed slower replication of dual-tropic HIV-1. hu-PBL-SCID mice derived from CCR5 delta 32/+ heterozygotes showed delayed replication of macrophage-tropic HIV-1 despite a small and variable effect of heterozygosity on viral replication in vitro. The level of CCR5 expression appears to limit replication of macrophage-tropic and dual-tropic HIV-1 strains in vivo.     

Faster Adaptation in Smaller Populations: Counterintuitive Evolution of HIV during Childhood Infection

Analysis of HIV-1 gene sequences sampled longitudinally from infected individuals can reveal the evolutionary dynamics that underlie associations between disease outcome and viral genetic diversity and divergence. Here we extend a statistical framework to estimate rates of viral molecular adaptation by considering sampling error when computing nucleotide site-frequencies. This is particularly beneficial when analyzing viral sequences from within-host viral infections if the number of sequences per time point is limited. To demonstrate the utility of this approach, we apply our method to a cohort of 24 patients infected with HIV-1 at birth. Our approach finds that viral adaptation arising from recurrent positive natural selection is associated with the rate of HIV-1 disease progression, in contrast to previous analyses of these data that found no significant association. Most surprisingly, we discover a strong negative correlation between viral population size and the rate of viral adaptation, the opposite of that predicted by standard molecular evolutionary theory. We argue that this observation is most likely due to the existence of a confounding third variable, namely variation in selective pressure among hosts. A conceptual non-linear model of virus adaptation that incorporates the two opposing effects of host immunity on the virus population can explain this counterintuitive result.     

Intra-Host Diversity and Emergence of Unique GBV-C Viral Lineages in HIV Infected Subjects in Central China

GB virus C (GBV-C), which is highly prevalent among HIV/AIDS, seemed to slow the HIV disease progression. The HIV/GBV-C co-infected individuals may represent an interesting model for the investigation of the role played by HIV infection and/or the immune system in driving the evolution of the GBV-C viral populations. The present study investigated the prevalence and population dynamics of GB virus C in HIV infected individuals representing 13 geographic regions of Hubei Province of China. Approximately 37% of HIV-1 infected individuals were infected with GBV-C and genotype 3 is appeared to be predominant. Utilizing the 196 complete E2 nucleotide sequence data from 10 HIV/GBV-C infected individuals and employing coalescence based phylogenetic approaches; the present study has investigated the intra-host dynamics of GBV-C. The results revealed patient-specific unique GBV-C viral lineages and each viral lineage showed the evidence of rapid population expansion in respective HIV-1 infected patients, thus suggesting HIV-1 was unlikely to have been inhibiting effect on the GBV-C viral replication. GBV-C in all patients has experienced intense purifying selection, suggesting the GBV-C viral invasion and subsequent expansion within the HIV-1 infected hosts without any modification of the functional epitopes at their membrane protein. The finding of within host GBV-C recombinant sequences indicated recombination was one of the significant forces in the evolution and divergence of GBV-C.     

Dynamic Correlation between Intrahost HIV-1 Quasispecies Evolution and Disease Progression

Quantifying the dynamics of intrahost HIV-1 sequence evolution is one means of uncovering information about the interaction between HIV-1 and the host immune system. In the chronic phase of infection, common dynamics of sequence divergence and diversity have been reported. We developed an HIV-1 sequence evolution model that simulated the effects of mutation and fitness of sequence variants. The amount of evolution was described by the distance from the founder strain, and fitness was described by the number of offspring a parent sequence produces. Analysis of the model suggested that the previously observed saturation of divergence and decrease of diversity in later stages of infection can be explained by a decrease in the proportion of offspring that are mutants as the distance from the founder strain increases rather than due to an increase of viral fitness. The prediction of the model was examined by performing phylogenetic analysis to estimate the change in the rate of evolution during infection. In agreement with our modeling, in 13 out of 15 patients (followed for 3–12 years) we found that the rate of intrahost HIV-1 evolution was not constant but rather slowed down at a rate correlated with the rate of CD4+ T-cell decline. The correlation between the dynamics of the evolutionary rate and the rate of CD4+ T-cell decline, coupled with our HIV-1 sequence evolution model, explains previously conflicting observations of the relationships between the rate of HIV-1 quasispecies evolution and disease progression.     

Potential Contributions of Viral Envelope and Host Genetic Factors in a Human Immunodeficiency Virus Type 1-Infected Long-Term Survivor

The lack of clinical progression in some individuals despite prolonged human immunodeficiency virus type 1 (HIV-1) infection may result from infection with less-pathogenic viral strains. To address this question, we examined the HIV-1 envelope protein from a donor with a low viral burden, stable CD4⁺ T-lymphocyte counts, and little evidence of CD8⁺ T-cell expansion, activation, or immune activity. To avoid potential changes in envelope function resulting from selection in vitro, envelope clones were constructed by using viral RNA isolated from uncultured peripheral blood mononuclear cells (PBMC). The data showed that recombinant viruses containing envelope sequences derived from RNA isolated from patient PBMC replicated poorly in primary CD4⁺ T cells but demonstrated efficient growth in macrophages. The unusual phenotype of these viruses could not be explained solely by differential utilization of coreceptors since the chimeric viruses, as well as an uncloned isolate obtained from the same visit date, can utilize CCR5. In addition, the donor’s own cells appeared resistant to infection with chimeric viruses containing autologous envelope sequences. Genotype analysis revealed that the donor was heterozygous for the previously described 32-bp deletion in CCR5 which may be linked with prolonged survival in HIV-1-infected individuals. These data suggest that the changes in envelope sequences confer properties of viral attenuation, which together with the CCR5 +/Δ32 genotype could account for the long-term survival of this patient.