Timing constraints of in vivo gag mutations during primary HIV-1 subtype C infection

No Difference in the Prevalence of HIV-1 gag Cytotoxic T-Lymphocyte-Associated Escape Mutations in Viral Sequences from Early and Late Parts of the HIV-1 Subtype C Pandemic in Botswana

HIV is known to accumulate escape mutations in the gag gene in response to the immune response from cytotoxic T lymphocytes (CTLs). These mutations can occur within an individual as well as at a population level. The population of Botswana exhibits a high prevalence of HLA*B57 and HLA*B58, which are associated with effective immune control of HIV. In this retrospective cross-sectional investigation, HIV-1 gag gene sequences were analyzed from recently infected participants across two time periods which were 10 years apart: the early time point (ETP) and late time point (LTP). The prevalence of CTL escape mutations was relatively similar between the two time points-ETP (10.6%) and LTP (9.7%). The P17 protein had the most mutations (9.4%) out of the 36 mutations that were identified. Three mutations (A83T, K18R, Y79H) in P17 and T190A in P24 were unique to the ETP sequences at a prevalence of 2.4%, 4.9%, 7.3%, and 5%, respectively. Mutations unique to the LTP sequences were all in the P24 protein, including T190V (3%), E177D (6%), R264K (3%), G248D (1%), and M228L (11%). Mutation K331R was statistically higher in the ETP (10%) compared to the LTP (1%) sequences (p < 0.01), while H219Q was higher in the LTP (21%) compared to the ETP (5%) (p < 0.01). Phylogenetically, the gag sequences clustered dependently on the time points. We observed a slower adaptation of HIV-1C to CTL immune pressure at a population level in Botswana. These insights into the genetic diversity and sequence clustering of HIV-1C can aid in the design of future vaccine strategies.

HIV-1 subtype C Nef-mediated SERINC5 down-regulation significantly contributes to overall Nef activity

BACKGROUND

Nef performs multiple cellular activities that enhance HIV-1 pathogenesis. The role of Nef-mediated down-regulation of the host restriction factor SERINC5 in HIV-1 pathogenesis is not well-defined. We aimed to investigate if SERINC5 down-regulation activity contributes to HIV-1 subtype C disease progression, to assess the relative contribution of this activity to overall Nef function, and to identify amino acids required for optimal activity. We measured the SERINC5 down-regulation activity of 106 subtype C Nef clones, isolated from individuals in early infection, for which the Nef activities of CD4 and HLA-I down-regulation as well as alteration of TCR signalling were previously measured. The relationship between SERINC5 down-regulation and markers of disease progression, and the relative contribution of SERINC5 down-regulation to a Nef fitness model-derived E value (a proxy for overall Nef fitness in vivo), were assessed.

RESULTS

No overall relationship was found between SERINC5 down-regulation and viral load set point (p = 0.28) or rate of CD4⁺ T cell decline (p = 0.45). CD4 down-regulation (p = 0.02) and SERINC5 down-regulation (p = 0.003) were significant determinants of E values in univariate analyses, with the greatest relative contribution for SERINC5 down-regulation, and only SERINC5 down-regulation remained significant in the multivariate analysis (p = 0.003). Using a codon-by-codon analysis, several amino acids were significantly associated with increased (10I, 11V, 38D, 51T, 65D, 101V, 188H and, 191H) or decreased (10K, 38E, 65E, 135F, 173T, 176T and, 191R) SERINC5 down-regulation activity. Site-directed mutagenesis experiments of selected mutants confirmed a substantial reduction in SERINC5 down-regulation activity associated with the mutation 173T, while mutations 10K, 135F, and 176T were associated with more modest reductions in activity that were not statistically significant.

CONCLUSIONS

These results suggest that SERINC5 down-regulation is a significant contributor to overall Nef function and identify potential genetic determinants of this Nef function that may have relevance for vaccines or therapeutics.

Genome-wide association study reveals genetic variants associated with HIV-1C infection in a Botswana study population

Although there have been many studies of gene variant association with different stages of HIV/AIDS progression in United States and European cohorts, few gene-association studies have assessed genic determinants in sub-Saharan African populations, which have the highest density of HIV infections worldwide. We carried out genome-wide association studies on 766 study participants at risk for HIV-1 subtype C (HIV-1C) infection in Botswana. Three gene associations (AP3B1, PTPRA, and NEO1) were shown to have significant association with HIV-1C acquisition. Each gene association was replicated within Botswana or in the United States-African American or United States-European American AIDS cohorts or in both. Each associated gene has a prior reported influence on HIV/AIDS pathogenesis. Thirteen previously discovered AIDS restriction genes were further replicated in the Botswana cohorts, extending our confidence in these prior AIDS restriction gene reports. This work presents an early step toward the identification of genetic variants associated with and affecting HIV acquisition or AIDS progression in the understudied HIV-1C afflicted Botswana population.

Phylogenetic approach to recover integration dates of latent HIV sequences within-host

Given that HIV evolution and latent reservoir establishment occur continually within-host, and that latently infected cells can persist long-term, the HIV reservoir should comprise a genetically heterogeneous archive recapitulating within-host HIV evolution. However, this has yet to be conclusively demonstrated, in part due to the challenges of reconstructing within-host reservoir establishment dynamics over long timescales. We developed a phylogenetic framework to reconstruct the integration dates of individual latent HIV lineages. The framework first involves inference and rooting of a maximum-likelihood phylogeny relating plasma HIV RNA sequences serially sampled before the initiation of suppressive antiretroviral therapy, along with putative latent sequences sampled thereafter. A linear model relating root-to-tip distances of plasma HIV RNA sequences to their sampling dates is used to convert root-to-tip distances of putative latent lineages to their establishment (integration) dates. Reconstruction of the ages of putative latent sequences sampled from chronically HIV-infected individuals up to 10 y following initiation of suppressive therapy revealed a genetically heterogeneous reservoir that recapitulated HIV's within-host evolutionary history. Reservoir sequences were interspersed throughout multiple within-host lineages, with the oldest dating to >20 y before sampling; historic genetic bottleneck events were also recorded therein. Notably, plasma HIV RNA sequences isolated from a viremia blip in an individual receiving otherwise suppressive therapy were highly genetically diverse and spanned a 20-y age range, suggestive of spontaneous in vivo HIV reactivation from a large latently infected cell pool. Our framework for reservoir dating provides a potentially powerful addition to the HIV persistence research toolkit.

Genome-Wide Analyses Reveal Gene Influence on HIV Disease Progression and HIV-1C Acquisition in Southern Africa

Sub-Saharan Africans infected with HIV-1C make up the largest AIDS patient population in the world and exhibit large heterogeneity in disease progression before initiating antiretroviral therapy. To identify host variants associated with HIV disease progression, we performed genome-wide association studies on a total of 556 treatment-naive HIV-infected individuals in Botswana. We characterized the pattern of HIV disease progression using a novel functional principal component analysis, which can better capture longitudinal CD4 and viral load (VL) trajectories. Two single-nucleotide polymorphisms (SNPs) near HCG22 (chr6, peak variant rs2535307, combined p = 3.72 × 10^-7, minor allele as risky allele) and CCNG1 (chr5, peak variant kgp22385164, combined p = 1.88 × 10^-6, minor allele as risky allele) were significantly associated with CD4 and VL dynamics. Inspection of SNPs in these gene regions in a third Botswana cohort (using GWATCH) also revealed a strong association of HCG22 with HIV-1C acquisition, suggesting that this region is associated with infection as well as disease progression. Our study uncovered two genetic regions that are significant and have specific effects on HIV-1C acquisition or progression in sub-Saharan Africans, and the result suggested new potential targets for AIDS prevention and treatment. In addition, our results also indicate the possibility of using genetic markers as HIV disease progression indicators in sub-Saharan Africans to prioritize fast progressors for antiretroviral treatment.

Pairwise diversity and tMRCA as potential markers for HIV infection recency

Intrahost human immunodeficiency virus (HIV)-1 diversity increases linearly over time. We assessed the extent to which mean pairwise distances and the time to the most recent common ancestor (tMRCA) inferred from intrahost HIV-1C env sequences were associated with the estimated time of HIV infection. Data from a primary HIV-1C infection study in Botswana were used for this analysis (N = 42). A total of 2540 HIV-1C env gp120 variable loop region 1 to conserved region 5 (V1C5) of the HIV-1 envelope gp120 viral sequences were generated by single genome amplification and sequencing, with an average of 61 viral sequences per participant and 11 sequences per time point per participant. Raw pairwise distances were calculated for each time point and participant using the ape package in R software. The tMRCA was estimated using phylogenetic inference implemented in Bayesian Evolutionary Analysis by Sampling Trees v1.8.2. Pairwise distances and tMRCA were significantly associated with the estimated time since HIV infection (both P < 0.001). Taking into account multiplicity of HIV infection strengthened these associations. HIV-1C env-based pairwise distances and tMRCA can be used as potential markers for HIV recency. However, the tMRCA estimates demonstrated no advantage over the pairwise distances estimates.

Deciphering Multiplicity of HIV-1C Infection: Transmission of Closely Related Multiple Viral Lineages

BACKGROUND

A single viral variant is transmitted in the majority of HIV infections. However, about 20% of heterosexually transmitted HIV infections are caused by multiple viral variants. Detection of transmitted HIV variants is not trivial, as it involves analysis of multiple viral sequences representing intra-host HIV-1 quasispecies.

METHODOLOGY

We distinguish two types of multiple virus transmission in HIV infection: (1) HIV transmission from the same source, and (2) transmission from different sources. Viral sequences representing intra-host quasispecies in a longitudinally sampled cohort of 42 individuals with primary HIV-1C infection in Botswana were generated by single-genome amplification and sequencing and spanned the V1C5 region of HIV-1C env gp120. The Maximum Likelihood phylogeny and distribution of pairwise raw distances were assessed at each sampling time point (n = 217; 42 patients; median 5 (IQR: 4-6) time points per patient, range 2-12 time points per patient).

RESULTS

Transmission of multiple viral variants from the same source (likely from the partner with established HIV infection) was found in 9 out of 42 individuals (21%; 95 CI 10-37%). HIV super-infection was identified in 2 patients (5%; 95% CI 1-17%) with an estimated rate of 3.9 per 100 person-years. Transmission of multiple viruses combined with HIV super-infection at a later time point was observed in one individual.

CONCLUSIONS

Multiple HIV lineages transmitted from the same source produce a monophyletic clade in the inferred phylogenetic tree. Such a clade has transiently distinct sub-clusters in the early stage of HIV infection, and follows a predictable evolutionary pathway. Over time, the gap between initially distinct viral lineages fills in and initially distinct sub-clusters converge. Identification of cases with transmission of multiple viral lineages from the same source needs to be taken into account in cross-sectional estimation of HIV recency in epidemiological and population studies.

Transmitted/Founder HIV-1 Subtype C Viruses Show Distinctive Signature Patterns in Vif, Vpr, and Vpu That Are Under Subsequent Immune Pressure During Early Infection

Viral variants that predominate during early infection may exhibit constrained diversity compared with those found during chronic infection and could contain amino acid signature patterns that may enhance transmission, establish productive infection, and influence early events that modulate the infection course. We compared amino acid distributions in 17 patients recently infected with HIV-1C with patients with chronic infection. We found significantly lower entropy in inferred transmitted/founder (t/f) compared with chronic viruses and identified signature patterns in Vif and Vpr from inferred t/f viruses. We investigated sequence evolution longitudinally up to 500 days postseroconversion and compared the impact of selected substitutions on predicted human leukocyte antigen (HLA) binding affinities of published and predicted cytotoxic T-lymphocyte epitopes. Polymorphisms in Vif and Vpr during early infection occurred more frequently at epitope-HLA anchor residues and significantly decreased predicted epitope-HLA binding. Transmission-associated sequence signatures may have implications for novel strategies to prevent HIV-1 transmission.

Assessing Biases in the Evaluation of Classification Assays for HIV Infection Recency

Identifying recent HIV infection cases has important public health and clinical implications. It is essential for estimating incidence rates to monitor epidemic trends and evaluate the effectiveness of interventions. Detecting recent cases is also important for HIV prevention given the crucial role that recently infected individuals play in disease transmission, and because early treatment onset can improve the clinical outlook of patients while reducing transmission risk. Critical to this enterprise is the development and proper assessment of accurate classification assays that, based on cross-sectional samples of viral sequences, help determine infection recency status. In this work we assess some of the biases present in the evaluation of HIV recency classification algorithms that rely on measures of within-host viral diversity. Particularly, we examine how the time since infection (TSI) distribution of the infected subjects from which viral samples are drawn affect performance metrics (e.g., area under the ROC curve, sensitivity, specificity, accuracy and precision), potentially leading to misguided conclusions about the efficacy of classification assays. By comparing the performance of a given HIV recency assay using six different TSI distributions (four simulated TSI distributions representing different epidemic scenarios, and two empirical TSI distributions), we show that conclusions about the overall efficacy of the assay depend critically on properties of the TSI distribution. Moreover, we demonstrate that an assay with high overall classification accuracy, mainly due to properly sorting members of the well-represented groups in the validation dataset, can still perform notoriously poorly when sorting members of the less represented groups. This is an inherent issue of classification and diagnostics procedures that is often underappreciated. Thus, this work underscores the importance of acknowledging and properly addressing evaluation biases when proposing new HIV recency assays.

A Generalized Entropy Measure of Within-Host Viral Diversity for Identifying Recent HIV-1 Infections

There is a need for incidence assays that accurately estimate HIV incidence based on cross-sectional specimens. Viral diversity-based assays have shown promises but are not particularly accurate. We hypothesize that certain viral genetic regions are more predictive of recent infection than others and aim to improve assay accuracy by using classification algorithms that focus on highly informative regions (HIRs).We analyzed HIV gag sequences from a cohort in Botswana. Forty-two subjects newly infected by HIV-1 Subtype C were followed through 500 days post-seroconversion. Using sliding window analysis, we screened for genetic regions within gag that best differentiate recent versus chronic infections. We used both nonparametric and parametric approaches to evaluate the discriminatory abilities of sequence regions. Segmented Shannon Entropy measures of HIRs were aggregated to develop generalized entropy measures to improve prediction of recency. Using logistic regression as the basis for our classification algorithm, we evaluated the predictive power of these novel biomarkers and compared them with recently reported viral diversity measures using area under the curve (AUC) analysis.Change of diversity over time varied across different sequence regions within gag. We identified the top 50% of the most informative regions by both nonparametric and parametric approaches. In both cases, HIRs were in more variable regions of gag and less likely in the p24 coding region. Entropy measures based on HIRs outperformed previously reported viral-diversity-based biomarkers. These methods are better suited for population-level estimation of HIV recency.The patterns of diversification of certain regions within the gag gene are more predictive of recency of infection than others. We expect this result to apply in other HIV genetic regions as well. Focusing on these informative regions, our generalized entropy measure of viral diversity demonstrates the potential for improving accuracy when identifying recent HIV-1 infections.

Exploring the benefits of antibody immune response in HIV-1 infection using a discrete model

The role of antibodies in HIV-1 infection is investigated using a discrete-time mathematical model that considers cell-free and cell-associated transmission of the virus. Model analysis shows that the effect of each type of antibody is dependent on the stage of the infection. Neutralizing antibodies are efficient in controlling the viral levels in the early days after seroconversion and antibodies that coat HIV-1-infected cells and recruit effector cells to either kill the HIV-1-infected cells or inhibit viral replication are efficient when the infection becomes established. Model simulations show that antibodies that inhibit viral replication are more effective in controlling the infection than those that recruit Natural Killer T cells after infection establishment. The model was fitted to subjects of the Tsedimoso study conducted in Botswana and conclusions similar to elasticity analysis results were obtained. Model fitting results predicted that neutralizing antibodies are more efficient in controlling the viral levels than antibodies that coat HIV-1-infected cells and recruit effector cells to either kill the HIV-1-infected cells or inhibit viral replication in the early days after seroconversion.

High frequency of transmitted HIV-1 Gag HLA class I-driven immune escape variants but minimal immune selection over the first year of clade C infection

In chronic HIV infection, CD8+ T cell responses to Gag are associated with lower viral loads, but longitudinal studies of HLA-restricted CD8+ T cell-driven selection pressure in Gag from the time of acute infection are limited. In this study we examined Gag sequence evolution over the first year of infection in 22 patients identified prior to seroconversion. A total of 310 and 337 full-length Gag sequences from the earliest available samples (median = 14 days after infection [Fiebig stage I/II]) and at one-year post infection respectively were generated. Six of 22 (27%) individuals were infected with multiple variants. There was a trend towards early intra-patient viral sequence diversity correlating with viral load set point (p = 0.07, r = 0.39). At 14 days post infection, 59.7% of Gag CTL epitopes contained non-consensus polymorphisms and over half of these (35.3%) comprised of previously described CTL escape variants. Consensus and variant CTL epitope proportions were equally distributed irrespective of the selecting host HLA allele and most epitopes remained unchanged over 12 months post infection. These data suggest that intrapatient diversity during acute infection is an indicator of disease outcome. In this setting, there is a high rate of transmitted CTL escape variants and limited immune selection in Gag during the first year of infection. These data have relevance for vaccine strategies designed to elicit effective CD8+ T cell immune responses.

Detection of HIV-1 matrix protein p17 quasispecies variants in plasma of chronic HIV-1-infected patients by ultra-deep pyrosequencing

BACKGROUND

The HIV-1 matrix protein p17 (p17MA) is a pleiotropic protein that plays a key role in the HIV-1 life cycle. It has been long believed to have a highly conserved primary amino acid sequence and a well-preserved structural integrity to avoid severe fitness consequences. However, recent data revealed that the carboxy (COOH)-terminus of p17MA possesses high levels of predicted intrinsic disorder, which would subtend to at least partially unfolded status of this region. This finding pointed to the need of investigating p17MA heterogeneity.

METHODS

The degree of intrapatient variations in the p17MA primary sequence was assessed on plasma viral RNA by using ultra-deep pyrosequencing.

RESULTS

Data obtained support a complex nature of p17MA quasispecies, with variants present at variable frequency virtually in all patients. Clusters of mutations were scattered along the entire sequence of the viral protein, but they were more frequently detected within the COOH-terminal region of p17MA. Moreover, deletions and insertions also occurred in a restricted area of the COOH-terminal region.

CONCLUSIONS

On the whole, our data show that the intrapatient level of sequence diversity in the p17MA is much higher than predicted before. Our results pave the way for further studies aimed at unraveling possible correlations between the presence of distinct p17MA variants and peculiar clinical evolutions of HIV-1 disease. The presence of p17MA quasispecies diversity may offer new tools to improve our understanding of the viral adaptation during the natural history of HIV-1 infection.

Nef-mediated down-regulation of CD4 and HLA class I in HIV-1 subtype C infection: association with disease progression and influence of immune pressure

Nef plays a major role in HIV-1 pathogenicity. We studied HIV-1 subtype C infected individuals in acute/early (n = 120) or chronic (n = 207) infection to investigate the relationship between Nef-mediated CD4/HLA-I down-regulation activities and disease progression, and the influence of immune-driven sequence variation on these Nef functions. A single Nef sequence per individual was cloned into an expression plasmid, followed by transfection of a T cell line and measurement of CD4 and HLA-I expression. In early infection, a trend of higher CD4 down-regulation ability correlating with higher viral load set point was observed (r = 0.19, p = 0.05), and higher HLA-I down-regulation activity was significantly associated with faster rate of CD4 decline (p = 0.02). HLA-I down-regulation function correlated inversely with the number HLA-associated polymorphisms previously associated with reversion in the absence of the selecting HLA allele (r = -0.21, p = 0.0002). These data support consideration of certain Nef regions in HIV-1 vaccine strategies designed to attenuate the infection course.

Reversion and T cell escape mutations compensate the fitness loss of a CD8+ T cell escape mutant in their cognate transmitted/founder virus

Immune escape mutations that revert back to the consensus sequence frequently occur in newly HIV-1-infected individuals and have been thought to render the viruses more fit. However, their impact on viral fitness and their interaction with other immune escape mutations have not been evaluated in the background of their cognate transmitted/founder (T/F) viral genomes. To precisely determine the role of reversion mutations, we introduced reversion mutations alone or together with CD8⁺ T cell escape mutations in their unmodified cognate T/F viral genome and determined their impact on viral fitness in primary CD4⁺ T cells. Two reversion mutations, V247I and I64T, were identified in Gag and Tat, respectively, but neither had measurable effect on the fitness of their cognate T/F virus. The V247I and G248A mutations that were detected before and concurrently with the potent T cell escape mutation T242N, respectively, were selected by early T cell responses. The V247I or the G248A mutation alone partially restored the fitness loss caused by the T242N mutation. Together they could fully restore the fitness of the T242N mutant to the T/F level. These results demonstrate that the fitness loss caused by a T cell escape mutation could be compensated by preexisting or concurrent reversion and other T cell escape mutations. Our findings indicate that the overall viral fitness is modulated by the complex interplay among T cell escape, compensatory and reversion mutations to maintain the balance between immune escape and viral replication capacity.

Intra-host evolutionary rates in HIV-1C env and gag during primary infection

BACKGROUND

HIV-1 nucleotide substitution rates are central for understanding the evolution of HIV-1. Their accurate estimation is critical for analysis of viral dynamics, identification of divergence time of HIV variants, inference of HIV transmission clusters, and modeling of viral evolution.

METHODS

Intra-patient nucleotide substitution rates in HIV-1C gag and env gp120 V1C5 were analyzed in a longitudinal cohort of 32 individuals infected with a single viral variant. Viral quasispecies were derived by single genome amplification/sequencing from serially sampled blood specimens collected at median (IQR) of 5 (4-6) times per subject from enrollment (during Fiebig stages II to V) over a median (IQR) of 417 (351-471) days post-seroconversion (p/s). HIV-1C evolutionary rates were estimated by BEAST v.1.7 using a relaxed lognormal molecular clock model. The effect of antiretroviral therapy (ART) on substitution rates in gag and env was assessed in a subset of six individuals who started ARV therapy during the follow-up period.

RESULTS

During primary HIV-1C infection, the intra-patient substitution rates were estimated at a median (IQR) of 5.22E-03 (3.28E-03-7.55E-03) substitutions per site per year of infection within gag, and 1.58E-02 (9.99E-03-2.04E-02) substitutions per site per year within env gp120 V1C5. The substitution rates in env gp120 V1C5 were higher than in gag (p<0.001, Wilcoxon signed rank test). The median (IQR) relative rates of evolution at codon positions 1, 2, and 3 were 0.73 (0.48-0.84), 0.67 (0.52-0.86), and 1.54 (1.21-1.71) in gag, and 1.01 (0.86-1.15), 1.05 (0.99-1.21), and 0.86 (0.67-0.94) in env gp120 V1C5, respectively. A first to the third position codon rate ratio >1.0 within env was found in 25 (78.1%) cases, but only in 4 (12.5%) cases in gag, while a second to the third position codon rate ratio >1.0 in env was observed in 26 (81.3%) cases, but in gag only in 2 (6.3%) cases (p<0.001 for both comparisons, Fisher's exact test). No ART effect on substitution rates in gag and env was found, at least within the first 3-4 months after ART initiation. Individuals with early viral set point ⩾4.0 log10 copies/ml had higher substitution rates in env gp120 V1C5 (median (IQR) 1.88E-02 (1.54E-02-2.46E-02) vs. 1.04E (7.24E-03-1.55E-02) substitutions per site per year; p=0.017, Mann-Whitney sum rank test), while individuals with early viral set point ⩾3.0 log10 copies/ml had higher substitution rates in gag (median (IQR) 5.66E-03 (3.45E-03-7.94E-03) vs. 1.78E-03 (4.57E-04-5.15E-03); p=0.028; Mann-Whitney sum rank test).

CONCLUSIONS

The results suggest that in primary HIV-1C infection, (1) intra-host evolutionary rates in env gp120 V1C5 are about 3-fold higher than in gag; (2) selection pressure in env is more frequent than in gag; (3) initiation of ART does not change substitution rates in HIV-1C env or gag, at least within the first 3-4 months after starting ART; and (4) intra-host evolutionary rates in gag and env gp120 V1C5 are higher in individuals with elevated levels of early viral set point.

tat Exon 1 exhibits functional diversity during HIV-1 subtype C primary infection

Human immunodeficiency virus type 1 (HIV-1) Tat is a mediator of viral transcription and is involved in the control of virus replication. However, associations between HIV-1 Tat diversity and functional effects during primary HIV-1 infection are still unclear. We estimated selection pressures in tat exon 1 using the mixed-effects model of evolution with 672 viral sequences generated from 20 patients infected with HIV-1 subtype C (HIV-1C) over 500 days postseroconversion. tat exon 1 residues 3, 4, 21, 24, 29, 39, and 68 were under positive selection, and we established that specific amino acid signature patterns were apparent in primary HIV-1C infection compared with chronic infection. We assessed the impact of these mutations on long terminal repeat (LTR) activity and found that Tat activity was negatively affected by the Ala(21) substitution identified in 13/20 (65%) of patients, which reduced LTR activity by 88% (± 1%) (P < 0.001). The greatest increase in Tat activity was seen with the Gln(35)/Lys(39) double mutant that resulted in an additional 49% (± 14%) production of LTR-driven luciferase (P = 0.012). There was a moderate positive correlation between Tat-mediated LTR activity and HIV-1 RNA in plasma (P = 0.026; r = 0.400) after 180 days postseroconversion that was reduced by 500 days postseroconversion (P = 0.043; r = 0.266). Although Tat activation of the LTR is not a strong predictor of these clinical variables, there are significant linear relationships between Tat transactivation and patients' plasma viral loads and CD4 counts, highlighting the complex interplay between Tat mutations in early HIV-1C infection.

The emerging role of HLA-C in HIV-1 infection

Recently, genome-wide association studies have identified the major histocompatibility complex class I protein HLA-C as an important molecule that affects HIV disease progression. The association between HLA-C and HIV disease outcome was originally determined through a single nucleotide polymorphism (SNP) 35 kb upstream of the HLA-C locus. More recent work has focused on elucidating the functional significance of the -35 SNP, and several groups now have demonstrated HLA-C surface expression to be a key element in control of HIV viral load, with higher surface expression associating with slower disease progression. Most recently, control of HLA-C surface expression has been correlated with the presence of microRNA binding sites that affect HLA-C expression and control of HIV disease. This review highlights these results and explores the ways in which HLA-C surface expression could affect immune system function in the setting of HIV disease.

Evolutionary gamut of in vivo Gag substitutions during early HIV-1 subtype C infection

Two analyses of HIV-1 subtype C Gag quasispecies were performed in a prospective cohort of 42 acutely and recently infected individuals by SGA on viral RNA/proviral DNA templates. First, in vivo Gag substitutions were assessed in relation to the HIV-1C consensus sequence, which revealed that 29.3% of detected amino acid substitutions can be classified as reversions to subtype consensus, 61.3% as forward substitutions from subtype consensus, and 9.3% as polymorphisms not associated with the subtype consensus sequence. Second, the proportion, dynamics, and relationships within individual pools of viral quasispecies were analyzed. Among reverse substitutions, 16.1% were minor, 11.0% transient, 13.6% dominant, and 59.2% fixed. In contrast, 31.6% of forward substitutions were minor, 59.3% transient, 3.8% dominant, and 5.3% fixed. The distinct patterns in the spectrum and dynamics of reverse and forward Gag substitutions suggest that these differences should be considered in HIV-1 evolutionary studies and analyses of viral mutational pathways.

Human leukocyte antigen variants B*44 and B*57 are consistently favorable during two distinct phases of primary HIV-1 infection in sub-Saharan Africans with several viral subtypes

As part of an ongoing study of early human immunodeficiency virus type 1 (HIV-1) infection in sub-Saharan African countries, we have identified 134 seroconverters (SCs) with distinct acute-phase (peak) and early chronic-phase (set-point) viremias. SCs with class I human leukocyte antigen (HLA) variants B*44 and B*57 had much lower peak viral loads (VLs) than SCs without these variants (adjusted linear regression beta values of -1.08 ± 0.26 log(10) [mean ± standard error] and -0.83 ± 0.27 log(10), respectively; P < 0.005 for both), after accounting for several nongenetic factors, including gender, age at estimated date of infection, duration of infection, and country of origin. These findings were confirmed by alternative models in which major viral subtypes (A1, C, and others) in the same SCs replaced country of origin as a covariate (P ≤ 0.03). Both B*44 and B*57 were also highly favorable (P ≤ 0.03) in analyses of set-point VLs. Moreover, B*44 was associated with relatively high CD4(+) T-cell counts during early chronic infection (P = 0.02). Thus, at least two common HLA-B variants showed strong influences on acute-phase as well as early chronic-phase VL, regardless of the infecting viral subtype. If confirmed, the identification of B*44 as another favorable marker in primary HIV-1 infection should help dissect mechanisms of early immune protection against HIV-1 infection.

Transmission of single and multiple viral variants in primary HIV-1 subtype C infection

To address whether sequences of viral gag and env quasispecies collected during the early post-acute period can be utilized to determine multiplicity of transmitted HIV's, recently developed approaches for analysis of viral evolution in acute HIV-1 infection [1], [2] were applied. Specifically, phylogenetic reconstruction, inter- and intra-patient distribution of maximum and mean genetic distances, analysis of Poisson fitness, shape of highlighter plots, recombination analysis, and estimation of time to the most recent common ancestor (tMRCA) were utilized for resolving multiplicity of HIV-1 transmission in a set of viral quasispecies collected within 50 days post-seroconversion (p/s) in 25 HIV-infected individuals with estimated time of seroconversion. The decision on multiplicity of HIV infection was made based on the model's fit with, or failure to explain, the observed extent of viral sequence heterogeneity. The initial analysis was based on phylogeny, inter-patient distribution of maximum and mean distances, and Poisson fitness, and was able to resolve multiplicity of HIV transmission in 20 of 25 (80%) cases. Additional analysis involved distribution of individual viral distances, highlighter plots, recombination analysis, and estimation of tMRCA, and resolved 4 of the 5 remaining cases. Overall, transmission of a single viral variant was identified in 16 of 25 (64%) cases, and transmission of multiple variants was evident in 8 of 25 (32%) cases. In one case multiplicity of HIV-1 transmission could not be determined. In primary HIV-1 subtype C infection, samples collected within 50 days p/s and analyzed by a single-genome amplification/sequencing technique can provide reliable identification of transmission multiplicity in 24 of 25 (96%) cases. Observed transmission frequency of a single viral variant and multiple viral variants were within the ranges of 64% to 68%, and 32% to 36%, respectively.

Influence of Gag-protease-mediated replication capacity on disease progression in individuals recently infected with HIV-1 subtype C

HLA class I-mediated selection of immune escape mutations in functionally important Gag epitopes may partly explain slower disease progression in HIV-1-infected individuals with protective HLA alleles. To investigate the impact of Gag function on disease progression, the replication capacities of viruses encoding Gag-protease from 60 individuals in early HIV-1 subtype C infection were assayed in an HIV-1-inducible green fluorescent protein reporter cell line and were correlated with subsequent disease progression. Replication capacities did not correlate with viral load set points (P = 0.37) but were significantly lower in individuals with below-median viral load set points (P = 0.03), and there was a trend of correlation between lower replication capacities and lower rates of CD4 decline (P = 0.09). Overall, the proportion of host HLA-specific Gag polymorphisms in or adjacent to epitopes was negatively associated with replication capacities (P = 0.04), but host HLA-B-specific polymorphisms were associated with higher viral load set points (P = 0.01). Further, polymorphisms associated with host-specific protective HLA alleles were linked with higher viral load set points (P = 0.03). These data suggest that transmission or early HLA-driven selection of Gag polymorphisms results in reduced early cytotoxic T-lymphocyte (CTL) responses and higher viral load set points. In support of the former, 46% of individuals with nonprotective alleles harbored a Gag polymorphism exclusively associated with a protective HLA allele, indicating a high rate of their transmission in sub-Saharan Africa. Overall, HIV disease progression is likely to be affected by the ability to mount effective Gag CTL responses as well as the replication capacity of the transmitted virus.

Identifying and characterizing recently transmitted viruses

PURPOSE OF REVIEW

Improvements in sequencing approaches and robust mathematical modeling have dramatically increased information on viral genetics during acute infection with HIV and simian immunodeficiency virus, providing unprecedented insight into viral transmission and viral/immune interactions.

RECENT FINDINGS

Overall viral genetic diversity is reduced significantly during mucosal transmission. Remarkably, in the vast majority of sexual transmissions, this diversity is reduced to a single viral variant that establishes the initial productive clinical infection. By identifying and enumerating transmitted/founder viruses, researchers can begin to define the characteristics that are necessary and sufficient for successful viral replication within a new host.

SUMMARY

Acute HIV infection is a critical window of opportunity for vaccine and therapeutic intervention. New sequencing technologies and mathematical modeling of transmission and early evolution have provided a clearer understanding of the number of founder viruses that establish infection, the rapid generation of diversity in these viruses and the subsequent evasion of host immunity. The information gained by identifying transmitted viruses, monitoring the initial host responses to these viruses and then identifying mechanisms of viral escape could provide better strategies for vaccine development, preexposure prophylaxis, microbicides, or other therapeutic interventions.

Dynamics and timing of in vivo mutations at Gag residue 242 during primary HIV-1 subtype C infection

Viral mutations at Gag residue 242 and relevant viral polymorphisms were analyzed in a cohort of 42 individuals with primary HIV-1 subtype C infection using single-genome amplification/sequencing. In HLA-B*57/5801-negative subjects infected with 242N escape variant, reversion to Asn appeared at median (IQR) 103 days (97-213 days) post-seroconversion (p/s) and became dominant at 193 days (170-215 days) p/s. In subjects expressing HLA-B*57/5801 and infected with the wild-type virus, the T242N escape appeared at 203 days (196-231) p/s, reached dominance at 277 days (265-315 days) p/s, and became complete at 323 days (289-373 days) p/s. HLA-B*57/5801-negative subjects infected with 242N escape variant did not show reduced viral load or increased CD4 count. The study highlights the differential selection of T242N escape by HLA-B*57 and B*5801 and suggests that the presence of HLA-B*57/5801-mediated immune pressure is able to control replication of the wild-type virus encoding Thr at Gag residue 242 but fails to suppress the T242N escape variant.

HIV-1 subtype C-infected individuals maintaining high viral load as potential targets for the "test-and-treat" approach to reduce HIV transmission

The first aim of the study is to assess the distribution of HIV-1 RNA levels in subtype C infection. Among 4,348 drug-naïve HIV-positive individuals participating in clinical studies in Botswana, the median baseline plasma HIV-1 RNA levels differed between the general population cohorts (4.1–4.2 log₁₀) and cART-initiating cohorts (5.1–5.3 log₁₀) by about one log₁₀. The proportion of individuals with high (≥50,000 (4.7 log₁₀) copies/ml) HIV-1 RNA levels ranged from 24%–28% in the general HIV-positive population cohorts to 65%–83% in cART-initiating cohorts. The second aim is to estimate the proportion of individuals who maintain high HIV-1 RNA levels for an extended time and the duration of this period. For this analysis, we estimate the proportion of individuals who could be identified by repeated 6- vs. 12-month-interval HIV testing, as well as the potential reduction of HIV transmission time that can be achieved by testing and ARV treating. Longitudinal analysis of 42 seroconverters revealed that 33% (95% CI: 20%–50%) of individuals maintain high HIV-1 RNA levels for at least 180 days post seroconversion (p/s) and the median duration of high viral load period was 350 (269; 428) days p/s. We found that it would be possible to identify all HIV-infected individuals with viral load ≥50,000 (4.7 log₁₀) copies/ml using repeated six-month-interval HIV testing. Assuming individuals with high viral load initiate cART after being identified, the period of high transmissibility due to high viral load can potentially be reduced by 77% (95% CI: 71%–82%). Therefore, if HIV-infected individuals maintaining high levels of plasma HIV-1 RNA for extended period of time contribute disproportionally to HIV transmission, a modified “test-and-treat” strategy targeting such individuals by repeated HIV testing (followed by initiation of cART) might be a useful public health strategy for mitigating the HIV epidemic in some communities.

Human leukocyte antigens and HIV type 1 viral load in early and chronic infection: predominance of evolving relationships

Background

During untreated, chronic HIV-1 infection, plasma viral load (VL) is a relatively stable quantitative trait that has clinical and epidemiological implications. Immunogenetic research has established various human genetic factors, especially human leukocyte antigen (HLA) variants, as independent determinants of VL set-point.

Methodology/Principal Findings

To identify and clarify HLA alleles that are associated with either transient or durable immune control of HIV-1 infection, we evaluated the relationships of HLA class I and class II alleles with VL among 563 seroprevalent Zambians (SPs) who were seropositive at enrollment and 221 seroconverters (SCs) who became seropositive during quarterly follow-up visits. After statistical adjustments for non-genetic factors (sex and age), two unfavorable alleles (A*3601 and DRB1*0102) were independently associated with high VL in SPs (p<0.01) but not in SCs. In contrast, favorable HLA variants, mainly A*74, B*13, B*57 (or Cw*18), and one HLA-A and HLA-C combination (A*30+Cw*03), dominated in SCs; their independent associations with low VL were reflected in regression beta estimates that ranged from −0.47±0.23 to −0.92±0.32 log₁₀ in SCs (p<0.05). Except for Cw*18, all favorable variants had diminishing or vanishing association with VL in SPs (p≤0.86).

Conclusions/Significance

Overall, each of the three HLA class I genes had at least one allele that might contribute to effective immune control, especially during the early course of HIV-1 infection. These observations can provide a useful framework for ongoing analyses of viral mutations induced by protective immune responses.