Covariability of selected amino acid positions for HIV type 1 subtypes C and B

Quantifying how single dose Ad26.COV2.S vaccine efficacy depends on Spike sequence features

In the ENSEMBLE randomized, placebo-controlled phase 3 trial (NCT04505722), estimated single-dose Ad26.COV2.S vaccine efficacy (VE) was 56% against moderate to severe-critical COVID-19. SARS-CoV-2 Spike sequences were determined from 484 vaccine and 1,067 placebo recipients who acquired COVID-19. In this set of prespecified analyses, we show that in Latin America, VE was significantly lower against Lambda vs. Reference and against Lambda vs. non-Lambda [family-wise error rate (FWER) p < 0.05]. VE differed by residue match vs. mismatch to the vaccine-insert at 16 amino acid positions (4 FWER p < 0.05; 12 q-value ≤ 0.20); significantly decreased with physicochemical-weighted Hamming distance to the vaccine-strain sequence for Spike, receptor-binding domain, N-terminal domain, and S1 (FWER p < 0.001); differed (FWER ≤ 0.05) by distance to the vaccine strain measured by 9 antibody-epitope escape scores and 4 NTD neutralization-impacting features; and decreased (p = 0.011) with neutralization resistance level to vaccinee sera. VE against severe-critical COVID-19 was stable across most sequence features but lower against the most distant viruses.

Quantifying how single dose Ad26.COV2.S vaccine efficacy depends on Spike sequence features

It is of interest to pinpoint SARS-CoV-2 sequence features defining vaccine resistance. In the ENSEMBLE randomized, placebo-controlled phase 3 trial, estimated single-dose Ad26.COV2.S vaccine efficacy (VE) was 56% against moderate to severe-critical COVID-19. SARS-CoV-2 Spike sequences were measured from 484 vaccine and 1,067 placebo recipients who acquired COVID-19 during the trial. In Latin America, where Spike diversity was greatest, VE was significantly lower against Lambda than against Reference and against all non-Lambda variants [family-wise error rate (FWER) p < 0.05]. VE also differed by residue match vs. mismatch to the vaccine-strain residue at 16 amino acid positions (4 FWER p < 0.05; 12 q-value ≤ 0.20). VE significantly decreased with physicochemical-weighted Hamming distance to the vaccine-strain sequence for Spike, receptor-binding domain, N-terminal domain, and S1 (FWER p < 0.001); differed (FWER ≤ 0.05) by distance to the vaccine strain measured by 9 different antibody-epitope escape scores and by 4 NTD neutralization-impacting features; and decreased (p = 0.011) with neutralization resistance level to vaccine recipient sera. VE against severe-critical COVID-19 was stable across most sequence features but lower against viruses with greatest distances. These results help map antigenic specificity of in vivo vaccine protection.

Prediction of VRC01 neutralization sensitivity by HIV-1 gp160 sequence features

The broadly neutralizing antibody (bnAb) VRC01 is being evaluated for its efficacy to prevent HIV-1 infection in the Antibody Mediated Prevention (AMP) trials. A secondary objective of AMP utilizes sieve analysis to investigate how VRC01 prevention efficacy (PE) varies with HIV-1 envelope (Env) amino acid (AA) sequence features. An exhaustive analysis that tests how PE depends on every AA feature with sufficient variation would have low statistical power. To design an adequately powered primary sieve analysis for AMP, we modeled VRC01 neutralization as a function of Env AA sequence features of 611 HIV-1 gp160 pseudoviruses from the CATNAP database, with objectives: (1) to develop models that best predict the neutralization readouts; and (2) to rank AA features by their predictive importance with classification and regression methods. The dataset was split in half, and machine learning algorithms were applied to each half, each analyzed separately using cross-validation and hold-out validation. We selected Super Learner, a nonparametric ensemble-based cross-validated learning method, for advancement to the primary sieve analysis. This method predicted the dichotomous resistance outcome of whether the IC50 neutralization titer of VRC01 for a given Env pseudovirus is right-censored (indicating resistance) with an average validated AUC of 0.868 across the two hold-out datasets. Quantitative log IC50 was predicted with an average validated R2 of 0.355. Features predicting neutralization sensitivity or resistance included 26 surface-accessible residues in the VRC01 and CD4 binding footprints, the length of gp120, the length of Env, the number of cysteines in gp120, the number of cysteines in Env, and 4 potential N-linked glycosylation sites; the top features will be advanced to the primary sieve analysis. This modeling framework may also inform the study of VRC01 in the treatment of HIV-infected persons.

Viral genetic diversity and protective efficacy of a tetravalent dengue vaccine in two phase 3 trials

Two phase 3 placebo-controlled trials of the CYD-TDV vaccine, evaluated in children aged 2-14 y (CYD14) and 9-16 y (CYD15), demonstrated vaccine efficacy (VE) of 56.5% and 60.8%, respectively, against symptomatic virologically confirmed dengue (VCD). Sieve analyses were conducted to evaluate whether and how VE varied with amino acid sequence features of dengue viruses (DENVs). DENV premembrane/envelope amino acid sequences from VCD endpoint cases were aligned with the vaccine insert sequences, and extensions of the proportional hazards model were applied to assess variation in VE with amino acid mismatch proportion distances from vaccine strains, individual amino acid residues, and phylogenetic genotypes. In CYD14, VE against VCD of any serotype (DENV-Any) decreased significantly with increasing amino acid distance from the vaccine, whereas in CYD15, VE against DENV-Any was distance-invariant. Restricting to the common age range and amino acid distance range between the trials and accounting for differential VE by serotype, however, showed no evidence of VE variation with distance in either trial. In serotype-specific analyses, VE against DENV4 decreased significantly with increasing amino acid distance from the DENV4 vaccine insert and was significantly greater against residue-matched DENV4 at eight signature positions. These effects were restricted to 2- to 8-y-olds, potentially because greater seropositivity of older children at baseline might facilitate a broader protective immune response. The relevance of an antigenic match between vaccine strains and circulating DENVs was also supported by greater estimated VE against serotypes and genotypes for which the circulating DENVs had shorter amino acid sequence distances from the vaccine.

HCV evolutionary genetics of SVR versus virologic failure assessed from the vaniprevir phase III registration trials

In the phase III registration studies conducted in Japan, Japanese HCV gt1 patients administered vaniprevir 300 mg twice daily plus pegylated interferon/ribavirin for 12 or 24 weeks achieved SVR24 rates of 83.7-84.5% among treatment-naïve patients, and 92.0-96.2% and 61.9% among breakthrough/relapsers or null-responders to prior interferon based therapy. As evidenced by direct sequencing, patients who did not achieve SVR24 principally failed due to treatment-emerging mutations at D168 or in a few cases R155. In this work, additional sequence analysis was conducted to address whether there were baseline polymorphisms associated with failure, evaluate the persistence of resistant virus among treatment failures, and assess for evidence of second site co-evolution with R155 or D168 mutations. To accomplish this, clonal sequencing (up to 40 clones per sample) was conducted on baseline, failure, and follow-up samples from all 38 patients among the vaniprevir treatment arms who met virologic failure criteria (37 gt1b, 1 gt1a, herein referred to as virologic failures) and baseline samples from 41 vaniprevir-treated SVR24 patients (all gt1b) selected among the three studies. SVR24 and virologic failure patients showed similar distributions of baseline polymorphisms previously associated with failure to one or more protease inhibitors. Furthermore, there was no evidence for baseline polymorphisms or a genetic signature across the NS3 protease domain specific to virologic failure patients, and which distinguishes them from baseline SVR24 sequences beyond a chance distribution. 24 of 32 virologic failures for whom baseline, failure, and follow-up samples were available showed reduced prevalence of the resistant virus first observed at the time of failure during the protocol-defined follow-up period of 24 weeks. Finally, pairwise analysis using either alignment or phylogenetic based methodologies provided no evidence for second site evolution with either the R155 or D168 mutations attributed to failure. This work supports and extends earlier findings based upon direct sequencing that attributed virologic failure to vaniprevir in the Phase III studies solely to the emergence of R155 or D168 mutations, with no apparent influence by other residues within the NS3 protease domain on treatment outcome. CLINICALTRIALS.

GOVIDENTIFIERS

NCT01370642, NCT01405937, NCT01405560.

Coevolution Analysis of HIV-1 Envelope Glycoprotein Complex

The HIV-1 Env spike is the main protein complex that facilitates HIV-1 entry into CD4⁺ host cells. HIV-1 entry is a multistep process that is not yet completely understood. This process involves several protein-protein interactions between HIV-1 Env and a variety of host cell receptors along with many conformational changes within the spike. HIV-1 Env developed due to high mutation rates and plasticity escape strategies from immense immune pressure and entry inhibitors. We applied a coevolution and residue-residue contact detecting method to identify coevolution patterns within HIV-1 Env protein sequences representing all group M subtypes. We identified 424 coevolving residue pairs within HIV-1 Env. The majority of predicted pairs are residue-residue contacts and are proximal in 3D structure. Furthermore, many of the detected pairs have functional implications due to contributions in either CD4 or coreceptor binding, or variable loop, gp120-gp41, and interdomain interactions. This study provides a new dimension of information in HIV research. The identified residue couplings may not only be important in assisting gp120 and gp41 coordinate structure prediction, but also in designing new and effective entry inhibitors that incorporate mutation patterns of HIV-1 Env.

Deep sequencing of protease inhibitor resistant HIV patient isolates reveals patterns of correlated mutations in Gag and protease

While the role of drug resistance mutations in HIV protease has been studied comprehensively, mutations in its substrate, Gag, have not been extensively cataloged. Using deep sequencing, we analyzed a unique collection of longitudinal viral samples from 93 patients who have been treated with therapies containing protease inhibitors (PIs). Due to the high sequence coverage within each sample, the frequencies of mutations at individual positions were calculated with high precision. We used this information to characterize the variability in the Gag polyprotein and its effects on PI-therapy outcomes. To examine covariation of mutations between two different sites using deep sequencing data, we developed an approach to estimate the tight bounds on the two-site bivariate probabilities in each viral sample, and the mutual information between pairs of positions based on all the bounds. Utilizing the new methodology we found that mutations in the matrix and p6 proteins contribute to continued therapy failure and have a major role in the network of strongly correlated mutations in the Gag polyprotein, as well as between Gag and protease. Although covariation is not direct evidence of structural propensities, we found the strongest correlations between residues on capsid and matrix of the same Gag protein were often due to structural proximity. This suggests that some of the strongest inter-protein Gag correlations are the result of structural proximity. Moreover, the strong covariation between residues in matrix and capsid at the N-terminus with p1 and p6 at the C-terminus is consistent with residue-residue contacts between these proteins at some point in the viral life cycle.

Understanding the structure of HIV proteins and the function of drug-resistant mutations of these proteins is critical for the development of effective HIV treatments. Selected gag mutations have been shown to provide compensatory functions for protease resistance mutations and may directly contribute to the development of drug resistance. To determine associations between protease inhibitor mutations and gag, we utilized deep sequencing of HIV gag and protease from a collection of viral isolates from patients treated with highly active retroviral protease inhibitors. Deep sequencing allows for accurate measurement of mutation frequencies at each position, allowing estimation, using a novel method we developed, of the covariation between any two residues on gag. Using this information, we characterize the variation within gag and protease and identify the most strongly correlated pairs of inter- and intra-protein residues. Our results suggest that matrix and p1/p6 mutations form the core of a network of strongly correlated gag mutations and contribute to recurrent treatment failure. Extracting gag residue covariation information from the deep sequencing of patient viral samples may provide insight into structural aspects of the Gag polyprotein as well new areas for small molecule targeting to disrupt Gag function.

Comprehensive sieve analysis of breakthrough HIV-1 sequences in the RV144 vaccine efficacy trial

The RV144 clinical trial showed the partial efficacy of a vaccine regimen with an estimated vaccine efficacy (VE) of 31% for protecting low-risk Thai volunteers against acquisition of HIV-1. The impact of vaccine-induced immune responses can be investigated through sieve analysis of HIV-1 breakthrough infections (infected vaccine and placebo recipients). A V1/V2-targeted comparison of the genomes of HIV-1 breakthrough viruses identified two V2 amino acid sites that differed between the vaccine and placebo groups. Here we extended the V1/V2 analysis to the entire HIV-1 genome using an array of methods based on individual sites, k-mers and genes/proteins. We identified 56 amino acid sites or "signatures" and 119 k-mers that differed between the vaccine and placebo groups. Of those, 19 sites and 38 k-mers were located in the regions comprising the RV144 vaccine (Env-gp120, Gag, and Pro). The nine signature sites in Env-gp120 were significantly enriched for known antibody-associated sites (p = 0.0021). In particular, site 317 in the third variable loop (V3) overlapped with a hotspot of antibody recognition, and sites 369 and 424 were linked to CD4 binding site neutralization. The identified signature sites significantly covaried with other sites across the genome (mean = 32.1) more than did non-signature sites (mean = 0.9) (p < 0.0001), suggesting functional and/or structural relevance of the signature sites. Since signature sites were not preferentially restricted to the vaccine immunogens and because most of the associations were insignificant following correction for multiple testing, we predict that few of the genetic differences are strongly linked to the RV144 vaccine-induced immune pressure. In addition to presenting results of the first complete-genome analysis of the breakthrough infections in the RV144 trial, this work describes a set of statistical methods and tools applicable to analysis of breakthrough infection genomes in general vaccine efficacy trials for diverse pathogens.

Vaccine-induced Human Antibodies Specific for the Third Variable Region of HIV-1 gp120 Impose Immune Pressure on Infecting Viruses

To evaluate the role of V3-specific IgG antibodies (Abs) in the RV144 clinical HIV vaccine trial, which reduced HIV-1 infection by 31.2%, the anti-V3 Ab response was assessed. Vaccinees' V3 Abs were highly cross-reactive with cyclic V3 peptides (cV3s) from diverse virus subtypes. Sieve analysis of CRF01_AE breakthrough viruses from 43 vaccine- and 66 placebo-recipients demonstrated an estimated vaccine efficacy of 85% against viruses with amino acids mismatching the vaccine at V3 site 317 (p = 0.004) and 52% against viruses matching the vaccine at V3 site 307 (p = 0.004). This analysis was supported by data showing that vaccinees' plasma Abs were less reactive with I³⁰⁷ when replaced with residues found more often in vaccinees' breakthrough viruses. Simultaneously, viruses with mutations at F³¹⁷ were less infectious, possibly due to the contribution of F³¹⁷ to optimal formation of the V3 hydrophobic core. These data suggest that RV144-induced V3-specific Abs imposed immune pressure on infecting viruses and inform efforts to design an HIV vaccine.

•

The RV144 vaccine reduced infection by viruses with isoleucine in V3 position 307.

•

Many vaccine-induced antibodies are cross-reactive and target an epitope including I³⁰⁷.

•

There was selection for breakthrough viruses carrying F³¹⁷ in V3 (p = 0.004).

•

F³¹⁷ is needed to maintain optimal infectivity.

•

F³¹⁷ is a poor or non-contact residue for vaccine induced V3 antibodies.

Increased HIV-1 vaccine efficacy against viruses with genetic signatures in Env V2

The RV144 trial demonstrated 31% vaccine efficacy at preventing human immunodeficiency virus (HIV)-1 infection. Antibodies against the HIV-1 envelope variable loops 1 and 2 (Env V1 and V2) correlated inversely with infection risk. We proposed that vaccine-induced immune responses against V1/V2 would have a selective effect against, or sieve, HIV-1 breakthrough viruses. A total of 936 HIV-1 genome sequences from 44 vaccine and 66 placebo recipients were examined. We show that vaccine-induced immune responses were associated with two signatures in V2 at amino acid positions 169 and 181. Vaccine efficacy against viruses matching the vaccine at position 169 was 48% (confidence interval 18% to 66%; P = 0.0036), whereas vaccine efficacy against viruses mismatching the vaccine at position 181 was 78% (confidence interval 35% to 93%; P = 0.0028). Residue 169 is in a cationic glycosylated region recognized by broadly neutralizing and RV144-derived antibodies. The predicted distance between the two signature sites (21 ± 7 Å) and their match/mismatch dichotomy indicate that multiple factors may be involved in the protection observed in RV144. Genetic signatures of RV144 vaccination in V2 complement the finding of an association between high V1/V2-binding antibodies and reduced risk of HIV-1 acquisition, and provide evidence that vaccine-induced V2 responses plausibly had a role in the partial protection conferred by the RV144 regimen.

Mutual information analysis reveals coevolving residues in Tat that compensate for two distinct functions in HIV-1 gene expression

Viral genomes are continually subjected to mutations, and functionally deleterious ones can be rescued by reversion or additional mutations that restore fitness. The error prone nature of HIV-1 replication has resulted in highly diverse viral sequences, and it is not clear how viral proteins such as Tat, which plays a critical role in viral gene expression and replication, retain their complex functions. Although several important amino acid positions in Tat are conserved, we hypothesized that it may also harbor functionally important residues that may not be individually conserved yet appear as correlated pairs, whose analysis could yield new mechanistic insights into Tat function and evolution. To identify such sites, we combined mutual information analysis and experimentation to identify coevolving positions and found that residues 35 and 39 are strongly correlated. Mutation of either residue of this pair into amino acids that appear in numerous viral isolates yields a defective virus; however, simultaneous introduction of both mutations into the heterologous Tat sequence restores gene expression close to wild-type Tat. Furthermore, in contrast to most coevolving protein residues that contribute to the same function, structural modeling and biochemical studies showed that these two residues contribute to two mechanistically distinct steps in gene expression: binding P-TEFb and promoting P-TEFb phosphorylation of the C-terminal domain in RNAPII. Moreover, Tat variants that mimic HIV-1 subtypes B or C at sites 35 and 39 have evolved orthogonal strengths of P-TEFb binding versus RNAPII phosphorylation, suggesting that subtypes have evolved alternate transcriptional strategies to achieve similar gene expression levels.

The hypervariable HIV-1 capsid protein residues comprise HLA-driven CD8+ T-cell escape mutations and covarying HLA-independent polymorphisms

One proposed HIV vaccine strategy is to induce Gag-specific CD8(+) T-cell responses that can corner the virus, through fitness cost of viral escape and unavailability of compensatory mutations. We show here that the most variable capsid residues principally comprise escape mutants driven by protective alleles HLA-B*57, -5801, and -8101 and covarying HLA-independent polymorphisms that arise in conjunction with these escape mutations. These covarying polymorphisms are potentially compensatory and are concentrated around three tropism-determining loops of p24, suggesting structural interdependencies. Our results demonstrate complex patterns of adaptation of HIV under immune selection pressure, the understanding of which should aid vaccine design.

Subtype-specific conservation of isoleucine 309 in the envelope V3 domain is linked to immune evasion in subtype C HIV-1 infection

The V3 region of the HIV-1 envelope (Env) glycoprotein gp120 is a key functional domain yet it exhibits distinct mutational patterns across subtypes. Here an invariant residue (Ile 309) was replaced with Leu in 7 subtype C patient-derived Envs from recent infection and 4 related neutralizing antibody escape variants that emerged later. For these 11 Envs, I309L did not alter replication in primary CD4 T cells; however, replication in monocyte-derived macrophages was enhanced. Infection of cell lines with low CD4 or CCR5 revealed that I309L enhanced utilization of CD4 but did not affect the ability to use CCR5. This CD4-enhanced phenotype tracked with sensitivity to sCD4, indicating increased exposure of the CD4 binding site. The results suggest that Ile 309 preserves a V3-mediated masking function that occludes the CD4 binding site. The findings point to an immune evasion strategy in subtype C Env to protect this vulnerable immune target.

Phylogenetic analysis of population-based and deep sequencing data to identify coevolving sites in the nef gene of HIV-1

Rapidly evolving viruses such as HIV-1 display extensive sequence variation in response to host-specific selection, while simultaneously maintaining functions that are critical to replication and infectivity. This apparent conflict between diversifying and purifying selection may be resolved by an abundance of epistatic interactions such that the same functional requirements can be met by highly divergent sequences. We investigate this hypothesis by conducting an extensive characterization of sequence variation in the HIV-1 nef gene that encodes a highly variable multifunctional protein. Population-based sequences were obtained from 686 patients enrolled in the HOMER cohort in British Columbia, Canada, from which the distribution of nonsynonymous substitutions in the phylogeny was reconstructed by maximum likelihood. We used a phylogenetic comparative method on these data to identify putative epistatic interactions between residues. Two interactions (Y120/Q125 and N157/S169) were chosen to further investigate within-host evolution using HIV-1 RNA extractions from plasma samples from eight patients. Clonal sequencing confirmed strong linkage between polymorphisms at these sites in every case. We used massively parallel pyrosequencing (MPP) to reconstruct within-host evolution in these patients. Experimental error associated with MPP was quantified by performing replicates at two different stages of the protocol, which were pooled prior to analysis to reduce this source of variation. Phylogenetic reconstruction from these data revealed correlated substitutions at Y120/Q125 or N157/S169 repeated across multiple lineages in every host, indicating convergent within-host evolution shaped by epistatic interactions.

Impact of amino acid variations in Gag and protease of HIV type 1 CRF01_AE strains on drug susceptibility of virus to protease inhibitors

BACKGROUND

Protease (PR) inhibitors (PIs) were designed against subtype B virus of human immunodeficiency virus type 1 (HIV-1), but believed to retain its activity against most of the other subtypes. CRF01_AE PR (AE-PR) contains background mutations that are presumed to alter the drug susceptibility of PR. In addition, amino acid variations found in HIV-1 Gag potentially affect the drug susceptibility or catalytic efficiency of PR.

METHODS

We studied the impact of naturally occurring amino acid substitutions found in AE-PR and CRF01_AE Gag (AE-Gag) on the drug susceptibility of PR to 9 currently available PIs, using the pNL4-3-derived luciferase reporter virus containing AE-Gag and/or AE-PR genes derived from drug treatment-naïve, HIV-1-infected Thai patients.

RESULTS

Sequencing analysis revealed that several mutations were detected in deduced amino acid sequences of AE-PR and AE-Gag genes, as compared to these genes of pNL4-3. Drug susceptibility tests revealed that AE-PR showed a variety of susceptibilities to 9 PIs compared with pNL4-3 PR. In addition, AE-Gag significantly reduced the drug susceptibility of AE-PR and pNL4-3 PR.

CONCLUSION

Our results suggest that amino acid variations in AE-PR and AE-Gag play roles in determining the drug susceptibility of CRF01_AE viruses to PIs.

Pairwise and higher-order correlations among drug-resistance mutations in HIV-1 subtype B protease

Background

The reaction of HIV protease to inhibitor therapy is characterized by the emergence of complex mutational patterns which confer drug resistance. The response of HIV protease to drugs often involves both primary mutations that directly inhibit the action of the drug, and a host of accessory resistance mutations that may occur far from the active site but may contribute to restoring the fitness or stability of the enzyme. Here we develop a probabilistic approach based on connected information that allows us to study residue, pair level and higher-order correlations within the same framework.

Results

We apply our methodology to a database of approximately 13,000 sequences which have been annotated by the treatment history of the patients from which the samples were obtained. We show that including pair interactions is essential for agreement with the mutational data, since neglect of these interactions results in order-of-magnitude errors in the probabilities of the simultaneous occurence of many mutations. The magnitude of these pair correlations changes dramatically between sequences obtained from patients that were or were not exposed to drugs. Higher-order effects make a contribution of as much as 10% for residues taken three at a time, but increase to more than twice that for 10 to 15-residue groups. The sequence data is insufficient to determine the higher-order effects for larger groups. We find that higher-order interactions have a significant effect on the predicted frequencies of sequences with large numbers of mutations. While relatively rare, such sequences are more prevalent after multi-drug therapy. The relative importance of these higher-order interactions increases with the number of drugs the patient had been exposed to.

Conclusion

Correlations are critical for the understanding of mutation patterns in HIV protease. Pair interactions have substantial qualitative effects, while higher-order interactions are individually smaller but may have a collective effect. Together they lead to correlations which could have an important impact on the dynamics of the evolution of cross-resistance, by allowing the virus to pass through otherwise unlikely mutational states. These findings also indicate that pairwise and possibly higher-order effects should be included in the models of protein evolution, instead of assuming that all residues mutate independently of one another.

Amino acid covariation in a functionally important human immunodeficiency virus type 1 protein region is associated with population subdivision

The frequently reported amino acid covariation of the highly polymorphic human immunodeficiency virus type 1 (HIV-1) exterior envelope glycoprotein V3 region has been assumed to reflect fitness epistasis between residues. However, nonrandom association of amino acids, or linkage disequilibrium, has many possible causes, including population subdivision. If the amino acids at a set of sequence sites differ in frequencies between subpopulations, then analysis of the whole population may reveal linkage disequilibrium even if it does not exist in any subpopulation. HIV-1 has a complex population structure, and the effects of this structure on linkage disequilibrium were investigated by estimating within- and among-subpopulation components of variance in linkage disequilibrium. The amino acid covariation previously reported is explained by differences in amino acid frequencies among virus subpopulations in different patients and by nonsystematic disequilibrium among patients. Disequilibrium within patients appears to be entirely due to differences in amino acid frequencies among sampling time points and among chemokine coreceptor usage phenotypes of virus particles, but not source tissues. Positive selection explains differences in allele frequencies among time points and phenotypes, indicating that these differences are adaptive rather than due to genetic drift. However, the absence of a correlation between linkage disequilibrium and phenotype suggests that fitness epistasis is an unlikely cause of disequilibrium. Indeed, when population structure is removed by analyzing sequences from a single time point and phenotype, no disequilibrium is detectable within patients. These results caution against interpreting amino acid covariation and coevolution as evidence for fitness epistasis.

Appreciating HIV type 1 diversity: subtype differences in Env

Human immunodeficiency virus type 1 (HIV-1) group M is responsible for the current AIDS pandemic and exhibits exceedingly high levels of viral genetic diversity around the world, necessitating categorization of viruses into distinct lineages, or subtypes. These subtypes can differ by around 35% in the envelope (Env) glycoproteins of the virus, which are displayed on the surface of the virion and are targets for both neutralizing antibody and cell-mediated immune responses. This diversity reflects the remarkable ability of the virus to adapt to selective pressures, the bulk of which is applied by the host immune response, and represents a serious obstacle for developing an effective vaccine with broad coverage. Thus, it is important to understand the underlying biological consequences of intersubtype diversity. Recent studies have revealed that some of the HIV-1 subtypes exhibit phenotypic differences stemming from subtle changes in Env structure, particularly within the highly immunogenic V3 domain, which participates directly in viral entry. This review will therefore explore current research that describes subtype differences in Env at the genetic and phenotypic level, focusing in particular on V3, and highlighting recent discoveries about the unique features of subtype C Env, which is the most globally prevalent subtype.

Predictive tests for linked changes

Mutations may confer a survival advantage to an organism and they can also reduce their fitness. In particular, we are interested in identifying correlated changes in genomic sequences. We consider the general situation where the observed characters at two genomic positions are summarized by an r x c contingency table. The test statistic focusses on double departures from the consensus configuration. When the original data are aggregated into two possible categories at each position (consensus vs non-consensus character), we obtain a 2 x 2 table to derive a test statistic that deals with the total number of double changes. Expected values and variances are predicted, under the assumption of independence, from table entries corresponding to single-mutation events. In some situations, the resulting tests are more powerful than those previously proposed.

Analysis of natural sequence variation and covariation in human immunodeficiency virus type 1 integrase

Human immunodeficiency virus type 1 (HIV-1) integrase inhibitors are in clinical trials, and raltegravir and elvitegravir are likely to be the first licensed drugs of this novel class of HIV antivirals. Understanding resistance to these inhibitors is important to maximize their efficacy. It has been shown that natural variation and covariation provide valuable insights into the development of resistance for established HIV inhibitors. Therefore, we have undertaken a study to fully characterize natural polymorphisms and amino acid covariation within an inhibitor-naïve sequence set spanning all defined HIV-1 subtypes. Inter- and intrasubtype variation was greatest in a 50-amino-acid segment of HIV-1 integrase incorporating the catalytic aspartic acid codon 116, suggesting that polymorphisms affect inhibitor binding and pathways to resistance. The critical mutations that determine the resistance pathways to raltegravir and elvitegravir (N155H, Q148K/R/H, and E92Q) were either rare or absent from the 1,165-sequence data set. However, 25 out of 41 mutations associated with integrase inhibitor resistance were present. These mutations were not subtype associated and were more prevalent in the subtypes that had been sampled frequently within the database. A novel modification of the Jaccard index was used to analyze amino acid covariation within HIV-1 integrase. A network of 10 covarying resistance-associated mutations was elucidated, along with a further 15 previously undescribed mutations that covaried with at least two of the resistance positions. The validation of covariation as a predictive tool will be dependent on monitoring the evolution of HIV-1 integrase under drug selection pressure.

Subtype-specific conformational differences within the V3 region of subtype B and subtype C human immunodeficiency virus type 1 Env proteins

The V3 region of the human immunodeficiency virus type 1 gp120 Env protein is a key domain in Env due to its role in interacting with the coreceptors CCR5 and CXCR4. We examined potential subtype-specific V3 region differences by comparing patterns of amino acid variability and probing for subtype-specific structures using 11 anti-V3 monoclonal antibodies (V3 MAbs). Differences between the subtypes in patterns of variability were most evident in the stem and turn regions of V3 (positions 9 to 24), with the two subtypes being very similar in the base region. The characteristics of the binding of V3 MAbs to Env proteins of the subtype B virus JR-FL and the subtype C virus BR025 suggested three patterns, as each group of MAbs recognized a specific conformation- or sequence-based epitope. Viruses pseudotyped with Env from JR-FL and BR025 were resistant to neutralization by the V3 MAbs, although the replacement of the Env V3 region of the SF162 virus with the JR-FL V3 created a pseudotyped virus that was hypersensitive to neutralization. A single mutation in V3 (H13R) made this chimeric Env selectively resistant to one group of V3 MAbs, consistent with the mAb binding properties. We hypothesize that there are intrinsic differences in V3 conformation between subtype B and subtype C that are localized to the stem and turn regions and that these differences have two important biological consequences: first, subtype B and subtype C V3 regions can have subtype-specific epitopes that will inherently limit antibody cross-reactivity, and second, V3 conformational differences may potentiate the frequent evolution of R5- into X4-tropic variants of subtype B but limit subtype C virus from using the same mechanism to evolve X4-tropic variants as efficiently.

An evolutionary-network model reveals stratified interactions in the V3 loop of the HIV-1 envelope

The third variable loop (V3) of the human immunodeficiency virus type 1 (HIV-1) envelope is a principal determinant of antibody neutralization and progression to AIDS. Although it is undoubtedly an important target for vaccine research, extensive genetic variation in V3 remains an obstacle to the development of an effective vaccine. Comparative methods that exploit the abundance of sequence data can detect interactions between residues of rapidly evolving proteins such as the HIV-1 envelope, revealing biological constraints on their variability. However, previous studies have relied implicitly on two biologically unrealistic assumptions: (1) that founder effects in the evolutionary history of the sequences can be ignored, and; (2) that statistical associations between residues occur exclusively in pairs. We show that comparative methods that neglect the evolutionary history of extant sequences are susceptible to a high rate of false positives (20%-40%). Therefore, we propose a new method to detect interactions that relaxes both of these assumptions. First, we reconstruct the evolutionary history of extant sequences by maximum likelihood, shifting focus from extant sequence variation to the underlying substitution events. Second, we analyze the joint distribution of substitution events among positions in the sequence as a Bayesian graphical model, in which each branch in the phylogeny is a unit of observation. We perform extensive validation of our models using both simulations and a control case of known interactions in HIV-1 protease, and apply this method to detect interactions within V3 from a sample of 1,154 HIV-1 envelope sequences. Our method greatly reduces the number of false positives due to founder effects, while capturing several higher-order interactions among V3 residues. By mapping these interactions to a structural model of the V3 loop, we find that the loop is stratified into distinct evolutionary clusters. We extend our model to detect interactions between the V3 and C4 domains of the HIV-1 envelope, and account for the uncertainty in mapping substitutions to the tree with a parametric bootstrap.

N-linked glycan modifications in gp120 of human immunodeficiency virus type 1 subtype C render partial sensitivity to 2G12 antibody neutralization

The monoclonal antibody (MAb) 2G12 recognizes a cluster of high-mannose oligosaccharides on the human immunodeficiency virus type 1 (HIV-1) envelope glycoprotein gp120 and is one of a select group of MAbs with broad neutralizing activity. However, subtype C viruses are generally resistant to 2G12 neutralization. This has been attributed to the absence of a glycosylation site at position 295 in most subtype C gp120s, which instead is typically occupied by a Val residue. Here we show that N-linked glycans in addition to the one at position 295 are important in the formation of the 2G12 epitope in subtype C gp120. Introduction of the glycosylation site at position 295 into three subtype C molecular clones, Du151.2, COT9.6, and COT6.15, did increase 2G12 binding to all three mutagenized gp120s, but at various levels. The COT9-V295N mutant showed the strongest 2G12 binding and was the only mutant to become sensitive to 2G12 neutralization, although very high antibody concentrations were required. Introduction of a glycosylation site at position 448 into mutant COT6-V295N, which occurs naturally in COT9, resulted in a virus that was partially sensitive to 2G12. Interestingly, a glycosylation site at position 442, which is common among subtype C viruses, also contributed to the 2G12 epitope. The addition of this glycan increased virus neutralization sensitivity to 2G12, whereas its deletion conferred resistance. Collectively, our results indicate that the 2G12 binding site cannot readily be reconstituted on the envelopes of subtype C viruses, suggesting structural differences from other HIV subtypes in which the 2G12 epitope is naturally expressed.

Distinguishing functional amino acid covariation from background linkage disequilibrium in HIV protease and reverse transcriptase

Correlated amino acid mutation analysis has been widely used to infer functional interactions between different sites in a protein. However, this analysis can be confounded by important phylogenetic effects broadly classifiable as background linkage disequilibrium (BLD). We have systematically separated the covariation induced by selective interactions between amino acids from background LD, using synonymous (S) vs. amino acid (A) mutations. Covariation between two amino acid mutations, (A,A), can be affected by selective interactions between amino acids, whereas covariation within (A,S) pairs or (S,S) pairs cannot. Our analysis of the pol gene — including the protease and the reverse transcriptase genes — in HIV reveals that (A,A) covariation levels are enormously higher than for either (A,S) or (S,S), and thus cannot be attributed to phylogenetic effects. The magnitude of these effects suggests that a large portion of (A,A) covariation in the HIV pol gene results from selective interactions. Inspection of the most prominent (A,A) interactions in the HIV pol gene showed that they are known sites of independently identified drug resistance mutations, and physically cluster around the drug binding site. Moreover, the specific set of (A,A) interaction pairs was reproducible in different drug treatment studies, and vanished in untreated HIV samples. The (S,S) covariation curves measured a low but detectable level of background LD in HIV.

CTL response to HIV type 1 subtype C is poorly predicted by known epitope motifs

Cytotoxic T lymphocyte (CTL) responses are thought to be essential for the control of HIV-1 replication in vivo and immunogens that elicit CTL responses are currently a major focus of HIV vaccine research. Here we investigated two aspects of the CTL response to HIV-1 subtype C that are important for vaccine design and efficacy monitoring. First, we assessed the relationship between the CTL response and sequence diversity, using a robust statistical method. While peptides that were most frequently recognized by the CTL response in Nef and p24 tended to be conserved, this was not the case for p17 where epitope recognition coincided with highly variable regions. Second, we investigated the relationship between observed and predicted CTL responses, given the HLA genotype of infected individuals. Only 52% of the Nef peptides and 64% of the Gag peptides that elicited a CTL response contained sequence motifs thought to be required for binding by the HLA-A or -B alleles found in the corresponding patient. In a comparable subtype B dataset a much higher proportion of the peptides that elicited a CTL response were consistent with the patient HLA genotype (96% and 83% for Nef and Gag, respectively). We demonstrate that this difference between subtypes C and B is likely to result from a combination of a tendency for HLA alleles common in Southern African populations to be poorly characterized, as well as a tendency for sequence motifs associated with HLA recognition to be overspecified for sequence variation found in the B clade. Our results suggest that knowledge of HLA binding motifs is likely to be biased toward certain populations and subtypes. This can have important implications for understanding immune escape and predicting vaccine efficacy in the context of populations primarily infected with non-B subtype HIV-1.

The reverse transcriptase 67N 70R 215Y genotype is the predominant TAM pathway associated with virologic failure among HIV type 1C-infected adults treated with ZDV/ddI-containing HAART in southern Africa

HIV-1C has become the dominant HIV-1 subtype in the global AIDS epidemic. Historically, the evolution of drug-resistant mutations was characterized primarily among antiretroviral (ARV)-treated HIV-1B infections. Whereas the non-B viruses are susceptible to the currently used ARVs, some differences between HIV-1 subtypes in response to ARV regimens have been reported. We analyzed the profile of ARV-associated mutations in HIV-1C infection treated with ZDV/ddI-containing regimens in an open-label, randomized 3 x 2 x 2 factorial study comparing ZDV/3TC vs. ZDV/ddI vs. d4T/3TC and EFV vs. NVP regimens in drug-naive adults in Botswana. The overall rate of virologic failure in the ZDV/ddI-containing arms was 14%. We addressed the development of NRTI-associated mutations in 23 virologically failed patients in the ZDV/ddI-containing arms. The 67N 70R 215Y genotype with wild-type amino acids at codon positions 41 and 210 was a dominant pattern of NRTI-associated mutations at the time of virologic failure. The mutation T215Y was the first step in the evolution of the 67N 70R 215Y genotype and was followed by mutations K70R and D67N. Representing a mixture of TAM-1 (41L/210W/215Y) and TAM-2 (67N/70R/215F /219Q) pathways, the 67N 70R 215Y genotype with wild-type amino acids at codon positions 41, 210, and 219 is a unique TAM pathway that is rarely seen in HIV-1B infection. Although limited by relatively small numbers, our data suggest that the 67N 70R 215Y genotype may be the HIV-1C-specific response to the first-line ZDV/ddI-containing regimen at the time of virologic failure. The presence of the 67N 70R 215Y genotype with wild-type amino acids at codon positions 41, 210, and 219 in HIV-1C infection suggests that the evolution of ARV-associated mutations and TAM pathways might be unique in non-B HIV-1 subtypes treated with particular ARV regimens.