Viral diversity and diversification of major non-structural genes vif, vpr, vpu, tat exon 1 and rev exon 1 during primary HIV-1 subtype C infection

Inferring the multiplicity of founder variants initiating HIV-1 infection: a systematic review and individual patient data meta-analysis

BACKGROUND

HIV-1 infections initiated by multiple founder variants are characterised by a higher viral load and a worse clinical prognosis than those initiated with single founder variants, yet little is known about the routes of exposure through which transmission of multiple founder variants is most probable. Here we used individual patient data to calculate the probability of multiple founders stratified by route of HIV exposure and study methodology.

METHODS

We conducted a systematic review and meta-analysis of studies that estimated founder variant multiplicity in HIV-1 infection, searching MEDLINE, Embase, and Global Health databases for papers published between Jan 1, 1990, and Sept 14, 2020. Eligible studies must have reported original estimates of founder variant multiplicity in people with acute or early HIV-1 infections, have clearly detailed the methods used, and reported the route of exposure. Studies were excluded if they reported data concerning people living with HIV-1 who had known or suspected superinfection, who were documented as having received pre-exposure prophylaxis, or if the transmitting partner was known to be receiving antiretroviral treatment. Individual patient data were collated from all studies, with authors contacted if these data were not publicly available. We applied logistic meta-regression to these data to estimate the probability that an HIV infection is initiated by multiple founder variants. We calculated a pooled estimate using a random effects model, subsequently stratifying this estimate across exposure routes in a univariable analysis. We then extended our model to adjust for different study methods in a multivariable analysis, recalculating estimates across the exposure routes. This study is registered with PROSPERO, CRD42020202672.

FINDINGS

We included 70 publications in our analysis, comprising 1657 individual patients. Our pooled estimate of the probability that an infection is initiated by multiple founder variants was 0·25 (95% CI 0·21-0·29), with moderate heterogeneity (Q=132·3, p<0·0001, I²=64·2%). Our multivariable analysis uncovered differences in the probability of multiple variant infection by exposure route. Relative to a baseline of male-to-female transmission, the predicted probability for female-to-male multiple variant transmission was significantly lower at 0·13 (95% CI 0·08-0·20), and the probabilities were significantly higher for transmissions in people who inject drugs (0·37 [0·24-0·53]) and men who have sex with men (0·30 [0·33-0·40]). There was no significant difference in the probability of multiple variant transmission between male-to-female transmission (0·21 [0·14-0·31]), post-partum transmission (0·18 [0·03-0·57]), pre-partum transmission (0·17 [0·08-0·33]), and intra-partum transmission (0·27 [0·14-0·45]).

INTERPRETATION

We identified that transmissions in people who inject drugs and men who have sex with men are significantly more likely to result in an infection initiated by multiple founder variants, and female-to-male infections are significantly less probable. Quantifying how the routes of HIV infection affect the transmission of multiple variants allows us to better understand how the evolution and epidemiology of HIV-1 determine clinical outcomes.

FUNDING

Medical Research Council Precision Medicine Doctoral Training Programme and a European Research Council Starting Grant.

Sequence and Functional Variation in the HIV-1 Rev Regulatory Axis

BACKGROUND

To complete its replication cycle, HIV-1 requires the nucleocytoplasmic export of intron-containing viral mRNAs. This process is ordinarily restricted by the cell, but HIV overcomes the block by means of a viral protein, Rev, and an RNA secondary structure found in all unspliced and incompletely spliced viral mRNAs called the Rev Response Element (RRE). In vivo activity of the Rev-RRE axis requires Rev binding to the RRE, oligomerization of Rev to form a competent ribonucleoprotein complex, and recruitment of cellular factors including Crm1 and RanGTP in order to export the targeted transcript. Sequence variability is observed among primary isolates in both Rev and the RRE, and the activity of both can be modulated through relatively small sequence changes. Primary isolates show differences in Rev-RRE activity and a few studies have found a correlation between lower Rev-RRE activity and slower progression of clinical disease. Lower Rev-RRE activity has also been associated with the evasion of cytotoxic T lymphocyte mediated killing.

CONCLUSION

The HIV-1 Rev-RRE regulatory axis is an understudied mechanism by which viral adaptation to diverse immune milieus may take place. There is evidence that this adaptation plays a role in HIV pathogenesis, particularly in immune evasion and latency, but further studies with larger sample sizes are warranted.

Combining Viral Genetics and Statistical Modeling to Improve HIV-1 Time-of-infection Estimation towards Enhanced Vaccine Efficacy Assessment

Knowledge of the time of HIV-1 infection and the multiplicity of viruses that establish HIV-1 infection is crucial for the in-depth analysis of clinical prevention efficacy trial outcomes. Better estimation methods would improve the ability to characterize immunological and genetic sequence correlates of efficacy within preventive efficacy trials of HIV-1 vaccines and monoclonal antibodies. We developed new methods for infection timing and multiplicity estimation using maximum likelihood estimators that shift and scale (calibrate) estimates by fitting true infection times and founder virus multiplicities to a linear regression model with independent variables defined by data on HIV-1 sequences, viral load, diagnostics, and sequence alignment statistics. Using Poisson models of measured mutation counts and phylogenetic trees, we analyzed longitudinal HIV-1 sequence data together with diagnostic and viral load data from the RV217 and CAPRISA 002 acute HIV-1 infection cohort studies. We used leave-one-out cross validation to evaluate the prediction error of these calibrated estimators versus that of existing estimators and found that both infection time and founder multiplicity can be estimated with improved accuracy and precision by calibration. Calibration considerably improved all estimators of time since HIV-1 infection, in terms of reducing bias to near zero and reducing root mean squared error (RMSE) to 5-10 days for sequences collected 1-2 months after infection. The calibration of multiplicity assessments yielded strong improvements with accurate predictions (ROC-AUC above 0.85) in all cases. These results have not yet been validated on external data, and the best-fitting models are likely to be less robust than simpler models to variation in sequencing conditions. For all evaluated models, these results demonstrate the value of calibration for improved estimation of founder multiplicity and of time since HIV-1 infection.

HIV-1 Tat interactions with cellular 7SK and viral TAR RNAs identifies dual structural mimicry

The HIV Tat protein competes with the 7SK:HEXIM interaction to hijack pTEFb from 7SK snRNP and recruit it to the TAR motif on stalled viral transcripts. Here we solve structures of 7SK stemloop-1 and TAR in complex with Tat’s RNA binding domain (RBD) to gain insights into this process. We find that 7SK is peppered with arginine sandwich motifs (ASM)—three classical and one with a pseudo configuration. Despite having similar RBDs, the presence of an additional arginine, R52, confers Tat the ability to remodel the pseudo configuration, required for HEXIM binding, into a classical sandwich, thus displacing HEXIM. Tat also uses R52 to remodel the TAR bulge into an ASM whose structure is identical to that of the remodeled ASM in 7SK. Together, our structures reveal a dual structural mimicry wherein viral Tat and TAR have co-opted structural motifs present in cellular HEXIM and 7SK for productive transcription of its genome.

The HIV Tat protein recruits a host elongation factor from the cellular 7SK complex to the viral TAR RNA to ensure transcriptional elongation. Here, Pham et al. solve the structures of both 7SK and TAR RNAs in complex with Tat’s RNA binding domain and gain mechanistic insights into the process.

HIV-1 Tat and Viral Latency: What We Can Learn from Naturally Occurring Sequence Variations

Despite the effective use of antiretroviral therapy, the remainder of a latently HIV-1-infected reservoir mainly in the resting memory CD4⁺ T lymphocyte subset has provided a great setback toward viral eradication. While host transcriptional silencing machinery is thought to play a dominant role in HIV-1 latency, HIV-1 protein such as Tat, may affect both the establishment and the reversal of latency. Indeed, mutational studies have demonstrated that insufficient Tat transactivation activity can result in impaired transcription of viral genes and the establishment of latency in cell culture experiments. Because Tat protein is one of highly variable proteins within HIV-1 proteome, it is conceivable that naturally occurring Tat mutations may differentially modulate Tat functions, thereby influencing the establishment and/or the reversal of viral latency in vivo. In this mini review, we summarize the recent findings of Tat naturally occurring polymorphisms associating with host immune responses and we highlight the implication of Tat sequence variations in relation to HIV latency.

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.

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.

Functional bottlenecks for generation of HIV-1 intersubtype Env recombinants

Background

Intersubtype recombination is a powerful driving force for HIV evolution, impacting both HIV-1 diversity within an infected individual and within the global epidemic. This study examines if viral protein function/fitness is the major constraint shaping selection of recombination hotspots in replication-competent HIV-1 progeny. A better understanding of the interplay between viral protein structure-function and recombination may provide insights into both vaccine design and drug development.

Results

In vitro HIV-1 dual infections were used to recombine subtypes A and D isolates and examine breakpoints in the Env glycoproteins. The entire env genes of 21 A/D recombinants with breakpoints in gp120 were non-functional when cloned into the laboratory strain, NL4-3. Likewise, cloning of A/D gp120 coding regions also produced dead viruses with non-functional Envs. 4/9 replication-competent viruses with functional Env’s were obtained when just the V1-V5 regions of these same A/D recombinants (i.e. same A/D breakpoints as above) were cloned into NL4-3.

Conclusion

These findings on functional A/D Env recombinants combined with structural models of Env suggest a conserved interplay between the C1 domain with C5 domain of gp120 and extracellular domain of gp41. Models also reveal a co-evolution within C1, C5, and ecto-gp41 domains which might explain the paucity of intersubtype recombination in the gp120 V1-V5 regions, despite their hypervariability. At least HIV-1 A/D intersubtype recombination in gp120 may result in a C1 from one subtype incompatible with a C5/gp41 from another subtype.

Transmitted/founder simian immunodeficiency virus envelope sequences in vesicular stomatitis and Semliki forest virus vector immunized rhesus macaques

Identification of transmitted/founder simian immunodeficiency virus (SIV) envelope sequences responsible for infection may prove critical for understanding HIV/SIV mucosal transmission. We used single genome amplification and phylogenetic analyses to characterize transmitted/founder SIVs both in the inoculum and in immunized-infected rhesus monkeys. Single genome amplification of the SIVsmE660 inoculum revealed a maximum diversity of 1.4%. We also noted that the consensus sequence of the challenge stock differed from the vaccine construct in 10 amino acids including 3 changes in the V4 loop. Viral env was prepared from rhesus plasma in 3 groups of 6 immunized with vesicular stomatitis virus (VSV) vectors and boosted with Semliki forest virus (SFV) replicons expressing (a) SIVsmE660 gag-env (b) SIVsmE660 gag-env plus rhesus GM-CSF and (c) control influenza hemagglutinin protein. Macaques were immunized twice with VSV-vectors and once with SFV vector and challenged intrarectally with 4000 TCID50. Single genome amplification characterized the infections of 2 unprotected animals in the gag-env immunized group, both of which had reduced acute plasma viral loads that ended as transient infections indicating partial immune control. Four of 6 rhesus were infected in the gag-env + GM-CSF group which demonstrated that GM-CSF abrogated protection. All 6 animals from the control group were infected having high plasma viral loads. We obtained 246 full-length envelope sequences from SIVsmE660 infected macaques at the peak of infection and determined the number of transmitted/founder variants per animal. Our analysis found that 2 of 2 gag-env vaccinated but infected macaques exhibited single but distinct virus envelope lineages whereas rhesus vaccinated with gag-env-GM-CSF or HA control exhibited both single and multiple env lineages. Because there were only 2 infected animals in the gag-env vaccinated rhesus compared to 10 infected rhesus in the other 2 groups, the significance of finding single env variants in the gag-env vaccinated group could not be established.

Role of host immune responses in sequence variability of HIV-1 Vpu

Viral protein U (Vpu) is an accessory protein associated with two main functions important in human immunodeficiency virus type 1 (HIV-1) replication and dissemination; these are down-regulation of CD4 receptor through mediating its proteasomal degradation and enhancement of virion release by antagonizing tetherin/BST2. It is also well established that Vpu is one of the most highly variable proteins in the HIV-1 proteome. However it is still unclear what drives Vpu sequence variability, whether Vpu acquires polymorphisms as a means of immune escape, functional advantage, or otherwise. It is assumed that the host-pathogen interaction is a cause of polymorphic phenotype of Vpu and that the resulting functional heterogeneity of Vpu may have critical significance in vivo. In order to comprehensively understand Vpu variability, it is important to integrate at the population level the genetic association approaches to identify specific amino acid residues and the immune escape kinetics which may impose Vpu functional constraints in vivo. This review will focus on HIV-1 accessory protein Vpu in the context of its sequence variability at population level and also bring forward evidence on the role of the host immune responses in driving Vpu sequence variability; we will also highlight the recent findings that illustrate Vpu functional implication in HIV-1 pathogenesis.

Can combination prevention strategies reduce HIV transmission in generalized epidemic settings in Africa? The HPTN 071 (PopART) study plan in South Africa and Zambia

The HIV Prevention Trials Network (HPTN) is conducting the HPTN 071 (PopART) study in 21 communities in Zambia and South Africa with support from a consortium of funders. HPTN 071 (PopART) is a community-randomized trial of a combination prevention strategy to reduce HIV incidence in the context of the generalized epidemic of southern Africa. The full PopART intervention strategy is anchored in home-based HIV testing and facilitated linkage of HIV-infected persons to care through community health workers and universal antiretroviral therapy for seropositive persons regardless of CD4+ cell count or HIV viral load. To further reduce the risk of HIV acquisition among uninfected individuals, the study aims to expand voluntary medical male circumcision, diagnosis and treatment of sexually transmitted infections, behavioral counseling, and condom distribution. The full PopART intervention strategy also incorporates promotion of other interventions designed to reduce HIV and tuberculosis transmission, including optimization of the prevention of mother-to-child HIV transmission and enhanced individual and public health tuberculosis services. Success for the PopART strategy depends on the ability to increase coverage for the study interventions whose uptake is a necessary antecedent to a prevention effect. Processes will be measured to assess the degree of penetration of the interventions into the communities. A randomly sampled population cohort from each community will be used to measure the impact of the PopART strategy on HIV incidence over 3 years. We describe the strategy being tested and progress to date in the HPTN 071 (PopART) study.

HIV populations are large and accumulate high genetic diversity in a nonlinear fashion

HIV infection is characterized by rapid and error-prone viral replication resulting in genetically diverse virus populations. The rate of accumulation of diversity and the mechanisms involved are under intense study to provide useful information to understand immune evasion and the development of drug resistance. To characterize the development of viral diversity after infection, we carried out an in-depth analysis of single genome sequences of HIV pro-pol to assess diversity and divergence and to estimate replicating population sizes in a group of treatment-naive HIV-infected individuals sampled at single (n = 22) or multiple, longitudinal (n = 11) time points. Analysis of single genome sequences revealed nonlinear accumulation of sequence diversity during the course of infection. Diversity accumulated in recently infected individuals at rates 30-fold higher than in patients with chronic infection. Accumulation of synonymous changes accounted for most of the diversity during chronic infection. Accumulation of diversity resulted in population shifts, but the rates of change were low relative to estimated replication cycle times, consistent with relatively large population sizes. Analysis of changes in allele frequencies revealed effective population sizes that are substantially higher than previous estimates of approximately 1,000 infectious particles/infected individual. Taken together, these observations indicate that HIV populations are large, diverse, and slow to change in chronic infection and that the emergence of new mutations, including drug resistance mutations, is governed by both selection forces and drift.

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.