Transmission of HIV-1 drug resistance mutations within partner-pairs: A cross-sectional study of a primary HIV infection cohort.

BACKGROUND: Transmission of human immunodeficiency virus type 1 (HIV-1) drug resistance mutations, particularly that of minority drug-resistant variants, remains poorly understood. Population-based studies suggest that drug-resistant HIV-1 is less transmissible than drug-susceptible viruses. We compared HIV-1 drug-resistant genotypes among partner-pairs in order to assess the likelihood of transmission of drug resistance mutations and investigate the role of minority variants in HIV transmission. METHODS AND FINDINGS: From 1992-2010, 340 persons with primary HIV-1 infection and their partners were enrolled into observational research studies at the University of Washington Primary Infection Clinic (UWPIC). Out of 50 partner-pairs enrolled, 36 (72%) transmission relationships were confirmed by phylogenetic distance analysis of HIV-1 envelope (env) sequences, and 31 partner-pairs enrolled after 1995 met criteria for this study. Drug resistance mutations in the region of the HIV-1 polymerase gene (pol) that encodes protease and reverse transcriptase were assessed by 454-pyrosequencing. In 25 partner-pairs where the transmission direction could be determined, 12 (48%) transmitters had 1-4 drug resistance mutations (23 total) detected in their HIV-1 populations at a median frequency of 6.0% (IQR 1.5%-98.7%, range 1.0%-99.6%). Of 10 major mutations detected in five transmitters at a frequency >95%, 100% (95% CI 69.2%-100%) were detected in recipients. All of these transmitters were antiretroviral (ARV)-naïve at the time of specimen collection. Fourteen mutations (eight major mutations and six accessory mutations) were detected in nine transmitters at low frequencies (1.0%-11.8%); four of these transmitters had previously received ARV therapy. Two (14% [95% CI 1.8%-42.8%]) G73S accessory mutations were detected in both transmitter and recipient. This number is not significantly different from the number expected based on the observed frequencies of drug-resistant viruses in transmitting partners. Limitations of this study include the small sample size and uncertainties in determining the timing of virus transmission and mutation history. CONCLUSIONS: Drug-resistant majority variants appeared to be commonly transmitted by ARV-naïve participants in our analysis and may contribute significantly to transmitted drug resistance on a population level. When present at low frequency, no major mutation was observed to be shared between partner-pairs; identification of accessory mutations shared within a pair could be due to transmission, laboratory artifact, or apolipoprotein B mRNA-editing enzyme, catalytic polypeptides (APOBECs), and warrants further study.

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.

Fitness-Balanced Escape Determines Resolution of Dynamic Founder Virus Escape Processes in HIV-1 Infection.

UNLABELLED: To understand the interplay between host cytotoxic T-lymphocyte (CTL) responses and the mechanisms by which HIV-1 evades them, we studied viral evolutionary patterns associated with host CTL responses in six linked transmission pairs. HIV-1 sequences corresponding to full-length p17 and p24 gag were generated by 454 pyrosequencing for all pairs near the time of transmission, and seroconverting partners were followed for a median of 847 days postinfection. T-cell responses were screened by gamma interferon/interleukin-2 (IFN-γ/IL-2) FluoroSpot using autologous peptide sets reflecting any Gag variant present in at least 5% of sequence reads in the individual's viral population. While we found little evidence for the occurrence of CTL reversions, CTL escape processes were found to be highly dynamic, with multiple epitope variants emerging simultaneously. We found a correlation between epitope entropy and the number of epitope variants per response (r = 0.43; P = 0.05). In cases in which multiple escape mutations developed within a targeted epitope, a variant with no fitness cost became fixed in the viral population. When multiple mutations within an epitope achieved fitness-balanced escape, these escape mutants were each maintained in the viral population. Additional mutations found to confer escape but undetected in viral populations incurred high fitness costs, suggesting that functional constraints limit the available sites tolerable to escape mutations. These results further our understanding of the impact of CTL escape and reversion from the founder virus in HIV infection and contribute to the identification of immunogenic Gag regions most vulnerable to a targeted T-cell attack. IMPORTANCE: Rapid diversification of the viral population is a hallmark of HIV-1 infection, and understanding the selective forces driving the emergence of viral variants can provide critical insight into the interplay between host immune responses and viral evolution. We used deep sequencing to comprehensively follow viral evolution over time in six linked HIV transmission pairs. We then mapped T-cell responses to explore if mutations arose due to adaption to the host and found that escape processes were often highly dynamic, with multiple mutations arising within targeted epitopes. When we explored the impact of these mutations on replicative capacity, we found that dynamic escape processes only resolve with the selection of mutations that conferred escape with no fitness cost to the virus. These results provide further understanding of the complicated viral-host interactions that occur during early HIV-1 infection and may help inform the design of future vaccine immunogens.

Analysis of HLA A*02 association with vaccine efficacy in the RV144 HIV-1 vaccine trial.

UNLABELLED: The RV144 HIV-1 vaccine trial demonstrated partial efficacy of 31% against HIV-1 infection. Studies into possible correlates of protection found that antibodies specific to the V1 and V2 (V1/V2) region of envelope correlated inversely with infection risk and that viruses isolated from trial participants contained genetic signatures of vaccine-induced pressure in the V1/V2 region. We explored the hypothesis that the genetic signatures in V1 and V2 could be partly attributed to selection by vaccine-primed T cells. We performed a T-cell-based sieve analysis of breakthrough viruses in the RV144 trial and found evidence of predicted HLA binding escape that was greater in vaccine versus placebo recipients. The predicted escape depended on class I HLA A*02- and A*11-restricted epitopes in the MN strain rgp120 vaccine immunogen. Though we hypothesized that this was indicative of postacquisition selection pressure, we also found that vaccine efficacy (VE) was greater in A*02-positive (A*02(+)) participants than in A*02(-) participants (VE = 54% versus 3%, P = 0.05). Vaccine efficacy against viruses with a lysine residue at site 169, important to antibody binding and implicated in vaccine-induced immune pressure, was also greater in A*02(+) participants (VE = 74% versus 15%, P = 0.02). Additionally, a reanalysis of vaccine-induced immune responses that focused on those that were shown to correlate with infection risk suggested that the humoral responses may have differed in A*02(+) participants. These exploratory and hypothesis-generating analyses indicate there may be an association between a class I HLA allele and vaccine efficacy, highlighting the importance of considering HLA alleles and host immune genetics in HIV vaccine trials. IMPORTANCE: The RV144 trial was the first to show efficacy against HIV-1 infection. Subsequently, much effort has been directed toward understanding the mechanisms of protection. Here, we conducted a T-cell-based sieve analysis, which compared the genetic sequences of viruses isolated from infected vaccine and placebo recipients. Though we hypothesized that the observed sieve effect indicated postacquisition T-cell selection, we also found that vaccine efficacy was greater for participants who expressed HLA A*02, an allele implicated in the sieve analysis. Though HLA alleles have been associated with disease progression and viral load in HIV-1 infection, these data are the first to suggest the association of a class I HLA allele and vaccine efficacy. While these statistical analyses do not provide mechanistic evidence of protection in RV144, they generate testable hypotheses for the HIV vaccine community and they highlight the importance of assessing the impact of host immune genetics in vaccine-induced immunity and protection. (This study has been registered at ClinicalTrials.gov under registration no. NCT00223080.).

Indel and Carryforward Correction (ICC): a new analysis approach for processing 454 pyrosequencing data.

MOTIVATION: Pyrosequencing technology provides an important new approach to more extensively characterize diverse sequence populations and detect low frequency variants. However, the promise of this technology has been difficult to realize, as careful correction of sequencing errors is crucial to distinguish rare variants (∼1%) in an infected host with high sensitivity and specificity. RESULTS: We developed a new approach, referred to as Indel and Carryforward Correction (ICC), to cluster sequences without substitutions and locally correct only indel and carryforward sequencing errors within clusters to ensure that no rare variants are lost. ICC performs sequence clustering in the order of (i) homopolymer indel patterns only, (ii) indel patterns only and (iii) carryforward errors only, without the requirement of a distance cutoff value. Overall, ICC removed 93-95% of sequencing errors found in control datasets. On pyrosequencing data from a PCR fragment derived from 15 HIV-1 plasmid clones mixed at various frequencies as low as 0.1%, ICC achieved the highest sensitivity and similar specificity compared with other commonly used error correction and variant calling algorithms. AVAILABILITY AND IMPLEMENTATION: Source code is freely available for download at http://indra.mullins.microbiol.washington.edu/ICC. It is implemented in Perl and supported on Linux, Mac OS X and MS Windows.

Cervicovaginal Human Papillomavirus Genomes, Microbiota Composition and Cytokine Concentrations in South African Adolescents.

The interaction between cervicovaginal virome, bacteriome and genital inflammation has not been extensively investigated. We assessed the vaginal DNA virome from 33 South African adolescents (15-19 years old) using shotgun DNA sequencing of purified virions. We present analyses of eukaryote-infecting DNA viruses, with a focus on human papillomavirus (HPV) genomes and relate these to the vaginal bacterial microbiota (assessed by 16S rRNA gene sequencing) and cytokines (assessed by Luminex). The DNA virome included single-stranded (Anelloviridae, Genomoviridae) and double-stranded DNA viruses (Adenoviridae, Alloherpesviridae, Herpesviridae, Marseilleviridae, Mimiviridae, Polyomaviridae, Poxviridae). We identified 110 unique, complete HPV genomes within two genera (Alphapapillomavirus and Gammapapillomavirus) representing 40 HPV types and 12 species. Of the 40 HPV types identified, 35 showed positive co-infection patterns with at least one other type, mainly HPV-16. HPV-35, a high-risk genotype currently not targeted by available vaccines, was the most prevalent HPV type identified in this cohort. Bacterial taxa commonly associated with bacterial vaginosis also correlated with the presence of HPV. Bacterial vaginosis, rather than HPV, was associated with increased genital inflammation. This study lays the foundation for future work characterizing the vaginal virome and its role in women's health.

Bacterial microbiome and host inflammatory gene expression in foreskin tissue.

The penile epithelial microbiome remains underexplored. We sequenced human RNA and a segment of the bacterial 16S rRNA gene from the foreskin tissue of 144 adolescents from South Africa and Uganda collected during penile circumcision after receipt of 1-2 doses of placebo, emtricitabine + tenofovir disoproxil fumarate, or emtricitabine + tenofovir alafenamide to investigate the microbiome of foreskin tissue and its potential changes with antiretroviral use. We identified a large number of anaerobic species, including Corynebacterium acnes, which was detected more frequently in participants from South Africa than Uganda. Bacterial populations did not differ by treatment received, and no differentially abundant taxa were identified between placebo versus active drug recipients. The relative abundance of specific bacterial taxa was negatively correlated with expression of genes downstream of the innate immune response to bacteria and regulation of inflammation. Our results show no difference in the tissue microbiome of the foreskin with short-course antiretroviral use but that bacterial taxa were largely inversely correlated with inflammatory gene expression, consistent with commensal colonization.

Demographic processes affect HIV-1 evolution in primary infection before the onset of selective processes.

HIV-1 transmission and viral evolution in the first year of infection were studied in 11 individuals representing four transmitter-recipient pairs and three independent seroconverters. Nine of these individuals were enrolled during acute infection; all were men who have sex with men (MSM) infected with HIV-1 subtype B. A total of 475 nearly full-length HIV-1 genome sequences were generated, representing on average 10 genomes per specimen at 2 to 12 visits over the first year of infection. Single founding variants with nearly homogeneous viral populations were detected in eight of the nine individuals who were enrolled during acute HIV-1 infection. Restriction to a single founder variant was not due to a lack of diversity in the transmitter as homogeneous populations were found in recipients from transmitters with chronic infection. Mutational patterns indicative of rapid viral population growth dominated during the first 5 weeks of infection and included a slight contraction of viral genetic diversity over the first 20 to 40 days. Subsequently, selection dominated, most markedly in env and nef. Mutants were detected in the first week and became consensus as early as day 21 after the onset of symptoms of primary HIV infection. We found multiple indications of cytotoxic T lymphocyte (CTL) escape mutations while reversions appeared limited. Putative escape mutations were often rapidly replaced with mutually exclusive mutations nearby, indicating the existence of a maturational escape process, possibly in adaptation to viral fitness constraints or to immune responses against new variants. We showed that establishment of HIV-1 infection is likely due to a biological mechanism that restricts transmission rather than to early adaptive evolution during acute infection. Furthermore, the diversity of HIV strains coupled with complex and individual-specific patterns of CTL escape did not reveal shared sequence characteristics of acute infection that could be harnessed for vaccine design.

Multistage genomewide association study identifies a locus at 1q41 associated with rate of HIV-1 disease progression to clinical AIDS.

BACKGROUND: A mean of 9-10 years of human immunodeficiency virus type 1 (HIV-1) infection elapse before clinical AIDS develops in untreated persons, but this rate of disease progression varies substantially among individuals. To investigate host genetic determinants of the rate of progression to clinical AIDS, we performed a multistage genomewide association study. METHODS: The discovery stage comprised 156 individuals from the Multicenter AIDS Cohort Study, enriched with rapid and long-term nonprogressors to increase statistical power. This was followed by replication tests of putatively associated genotypes in an independent population of 590 HIV-1-infected seroconverters. RESULTS: Significant associations with delayed AIDS progression were observed in a haplotype located at 1q41, 36 kb upstream of PROX1 on chromosome 1 (relative hazard ratio, 0.69; Fisher's combined P = 6.23 X 10(-7)). This association was replicated further in an analysis stratified by transmission mode, with the effect consistent in sexual or mucosal and parenteral transmission (relative hazard ratios, 0.72 and 0.63, respectively; combined P = 1.63 X 10(-6)). CONCLUSIONS: This study identified and replicated a locus upstream of PROX1 that is associated with delayed progression to clinical AIDS. PROX1 is a negative regulator of interferon-gamma expression in T cells and also mitigates the advancement of vascular neoplasms, such as Kaposi sarcoma, a common AIDS-defining malignancy. This study adds to the cumulative polygenic host component that effectively regulates the progression to clinical AIDS among HIV-1-infected individuals, raising prospects for potential new avenues for therapy and improvements in AIDS prognosis.

Presence and Persistence of Putative Lytic and Temperate Bacteriophages in Vaginal Metagenomes from South African Adolescents.

The interaction between gut bacterial and viral microbiota is thought to be important in human health. While fluctuations in female genital tract (FGT) bacterial microbiota similarly determine sexual health, little is known about the presence, persistence, and function of vaginal bacteriophages. We conducted shotgun metagenome sequencing of cervicovaginal samples from South African adolescents collected longitudinally, who received no antibiotics. We annotated viral reads and circular bacteriophages, identified CRISPR loci and putative prophages, and assessed their diversity, persistence, and associations with bacterial microbiota composition. Siphoviridae was the most prevalent bacteriophage family, followed by Myoviridae, Podoviridae, Herelleviridae, and Inoviridae. Full-length siphoviruses targeting bacterial vaginosis (BV)-associated bacteria were identified, suggesting their presence in vivo. CRISPR loci and prophage-like elements were common, and genomic analysis suggested higher diversity among Gardnerella than Lactobacillus prophages. We found that some prophages were highly persistent within participants, and identical prophages were present in cervicovaginal secretions of multiple participants, suggesting that prophages, and thus bacterial strains, are shared between adolescents. The number of CRISPR loci and prophages were associated with vaginal microbiota stability and absence of BV. Our analysis suggests that (pro)phages are common in the FGT and vaginal bacteria and (pro)phages may interact.

HLA class I-driven evolution of human immunodeficiency virus type 1 subtype c proteome: immune escape and viral load.

Human immunodeficiency virus type 1 (HIV-1) mutations that confer escape from cytotoxic T-lymphocyte (CTL) recognition can sometimes result in lower viral fitness. These mutations can then revert upon transmission to a new host in the absence of CTL-mediated immune selection pressure restricted by the HLA alleles of the prior host. To identify these potentially critical recognition points on the virus, we assessed HLA-driven viral evolution using three phylogenetic correction methods across full HIV-1 subtype C proteomes from a cohort of 261 South Africans and identified amino acids conferring either susceptibility or resistance to CTLs. A total of 558 CTL-susceptible and -resistant HLA-amino acid associations were identified and organized into 310 immunological sets (groups of individual associations related to a single HLA/epitope combination). Mutations away from seven susceptible residues, including four in Gag, were associated with lower plasma viral-RNA loads (q < 0.2 [where q is the expected false-discovery rate]) in individuals with the corresponding HLA alleles. The ratio of susceptible to resistant residues among those without the corresponding HLA alleles varied in the order Vpr > Gag > Rev > Pol > Nef > Vif > Tat > Env > Vpu (Fisher's exact test; P < or = 0.0009 for each comparison), suggesting the same ranking of fitness costs by genes associated with CTL escape. Significantly more HLA-B (chi(2); P = 3.59 x 10(-5)) and HLA-C (chi(2); P = 4.71 x 10(-6)) alleles were associated with amino acid changes than HLA-A, highlighting their importance in driving viral evolution. In conclusion, specific HIV-1 residues (enriched in Vpr, Gag, and Rev) and HLA alleles (particularly B and C) confer susceptibility to the CTL response and are likely to be important in the development of vaccines targeted to decrease the viral load.

ViroBLAST: a stand-alone BLAST web server for flexible queries of multiple databases and user's datasets.

UNLABELLED: ViroBLAST is a stand-alone BLAST web interface for nucleotide and amino acid sequence similarity searches. It extends the utility of BLAST to query against multiple sequence databases and user sequence datasets, and provides a friendly output to easily parse and navigate BLAST results. ViroBLAST is readily useful for all research areas that require BLAST functions and is available online and as a downloadable archive for independent installation. AVAILABILITY: http://indra.mullins.microbiol.washington.edu/blast/viroblast.php.

Pairwise growth competition assay for determining the replication fitness of human immunodeficiency viruses.

In vitro fitness assays are essential tools for determining viral replication fitness for viruses such as HIV-1. Various measurements have been used to extrapolate viral replication fitness, ranging from the number of viral particles per infectious unit, growth rate in cell culture, and relative fitness derived from multiple-cycle growth competition assays. Growth competition assays provide a particularly sensitive measurement of fitness since the viruses are competing for cellular targets under identical growth conditions. There are several experimental factors to consider when conducting growth competition assays, including the multiplicity of infection (MOI), sampling times, and viral detection and fitness calculation methods. Each factor can affect the end result and hence must be considered carefully during the experimental design. The protocol presented here includes steps from constructing a new recombinant HIV-1 clone to performing growth competition assays and analyzing the experimental results. This protocol utilizes experimental parameter values previously shown to yield consistent and robust results. Alternatives are discussed, as some parameters need to be adjusted according to the cell type and viruses being studied. The protocol contains two alternative viral detection methods to provide flexibility as the availability of instruments, reagents and expertise varies between laboratories.

CD8 and CD4 epitope predictions in RV144: no strong evidence of a T-cell driven sieve effect in HIV-1 breakthrough sequences from trial participants.

The modest protection afforded by the RV144 vaccine offers an opportunity to evaluate its mechanisms of protection. Differences between HIV-1 breakthrough viruses from vaccine and placebo recipients can be attributed to the RV144 vaccine as this was a randomized and double-blinded trial. CD8 and CD4 T cell epitope repertoires were predicted in HIV-1 proteomes from 110 RV144 participants. Predicted Gag epitope repertoires were smaller in vaccine than in placebo recipients (p = 0.019). After comparing participant-derived epitopes to corresponding epitopes in the RV144 vaccine, the proportion of epitopes that could be matched differed depending on the protein conservation (only 36% of epitopes in Env vs 84-91% in Gag/Pol/Nef for CD8 predicted epitopes) or on vaccine insert subtype (55% against CRF01_AE vs 7% against subtype B). To compare predicted epitopes to the vaccine, we analyzed predicted binding affinity and evolutionary distance measurements. Comparisons between the vaccine and placebo arm did not reveal robust evidence for a T cell driven sieve effect, although some differences were noted in Env-V2 (0.022≤p-value≤0.231). The paucity of CD8 T cell responses identified following RV144 vaccination, with no evidence for V2 specificity, considered together both with the association of decreased infection risk in RV 144 participants with V-specific antibody responses and a V2 sieve effect, lead us to hypothesize that this sieve effect was not T cell specific. Overall, our results did not reveal a strong differential impact of vaccine-induced T cell responses among breakthrough infections in RV144 participants.

Factors affecting relative fitness measurements in pairwise competition assays of human immunodeficiency viruses.

Cell culture growth competition assays of human immunodeficiency virus type 1 (HIV-1) are used to estimate viral fitness and quantify the impact of mutations conferring drug resistance and immunological escape. A comprehensive study of growth competition assays was conducted and identified experimental parameters that can impact measurements of relative fitness including multiplicity of infection, viral input ratio, number, timing and interval of time points used to evaluate selective outgrowth, and the algorithm for calculating fitness values. An optimized protocol is developed here that is a multi-point growth competition assay that resolves reproducibly small differences in viral fitness. The optimized protocol uses an MOI of 0.005, a consistent ratio of mutant: parental viruses (70:30), and a multipoint [1+s 4,7] algorithm that uses data points within the logarithmic phase of viral growth for assessing fitness differences.

Improved detection of rare HIV-1 variants using 454 pyrosequencing.

454 pyrosequencing, a massively parallel sequencing (MPS) technology, is often used to study HIV genetic variation. However, the substantial mismatch error rate of the PCR required to prepare HIV-containing samples for pyrosequencing has limited the detection of rare variants within viral populations to those present above ~1%. To improve detection of rare variants, we varied PCR enzymes and conditions to identify those that combined high sensitivity with a low error rate. Substitution errors were found to vary up to 3-fold between the different enzymes tested. The sensitivity of each enzyme, which impacts the number of templates amplified for pyrosequencing, was shown to vary, although not consistently across genes and different samples. We also describe an amplicon-based method to improve the consistency of read coverage over stretches of the HIV-1 genome. Twenty-two primers were designed to amplify 11 overlapping amplicons in the HIV-1 clade B gag-pol and env gp120 coding regions to encompass 4.7 kb of the viral genome per sample at sensitivities as low as 0.01-0.2%.

Is the virulence of HIV changing? A meta-analysis of trends in prognostic markers of HIV disease progression and transmission.

OBJECTIVE: The potential for changing HIV-1 virulence has significant implications for the AIDS epidemic, including changing HIV transmission rates, rapidity of disease progression, and timing of ART. Published data to date have provided conflicting results. DESIGN: We conducted a meta-analysis of changes in baseline CD4(+) T-cell counts and set point plasma viral RNA load over time in order to establish whether summary trends are consistent with changing HIV-1 virulence. METHODS: We searched PubMed for studies of trends in HIV-1 prognostic markers of disease progression and supplemented findings with publications referenced in epidemiological or virulence studies. We identified 12 studies of trends in baseline CD4(+) T-cell counts (21, 052 total individuals), and eight studies of trends in set point viral loads (10 ,785 total individuals), spanning the years 1984-2010. Using random-effects meta-analysis, we estimated summary effect sizes for trends in HIV-1 plasma viral loads and CD4(+) T-cell counts. RESULTS: Baseline CD4(+) T-cell counts showed a summary trend of decreasing cell counts [effect = -4.93  cells/µl per year, 95% confidence interval (CI) -6.53 to -3.3]. Set point viral loads showed a summary trend of increasing plasma viral RNA loads (effect = 0.013  log(10)  copies/ml per year, 95% CI -0.001 to 0.03). The trend rates decelerated in recent years for both prognostic markers. CONCLUSION: Our results are consistent with increased virulence of HIV-1 over the course of the epidemic. Extrapolating over the 30 years since the first description of AIDS, this represents a CD4(+) T cells loss of approximately 148  cells/µl and a gain of 0.39  log(10)  copies/ml of viral RNA measured during early infection. These effect sizes would predict increasing rates of disease progression, and need for ART as well as increasing transmission risk.

Genetic impact of vaccination on breakthrough HIV-1 sequences from the STEP trial.

We analyzed HIV-1 genome sequences from 68 newly infected volunteers in the STEP HIV-1 vaccine trial. To determine whether the vaccine exerted selective T cell pressure on breakthrough viruses, we identified potential T cell epitopes in the founder sequences and compared them to epitopes in the vaccine. We found greater distances to the vaccine sequence for sequences from vaccine recipients than from placebo recipients. The most significant signature site distinguishing vaccine from placebo recipients was Gag amino acid 84, a site encompassed by several epitopes contained in the vaccine and restricted by human leukocyte antigen (HLA) alleles common in the study cohort. Moreover, the extended divergence was confined to the vaccine components of the virus (HIV-1 Gag, Pol and Nef) and not found in other HIV-1 proteins. These results represent what is to our knowledge the first evidence of selective pressure from vaccine-induced T cell responses on HIV-1 infection in humans.

DIVEIN: a web server to analyze phylogenies, sequence divergence, diversity, and informative sites.

DIVEIN is a web interface that performs automated phylogenetic and other analyses of nucleotide and amino acid sequences. Starting with a set of aligned sequences, DIVEIN estimates evolutionary parameters and phylogenetic trees while allowing the user to choose from a variety of evolutionary models; it then reconstructs the consensus (CON), most recent common ancestor (MRCA), and center of tree (COT) sequences. DIVEIN also provides tools for further analyses, including condensing sequence alignments to show only informative sites or private mutations; computing phylogenetic or pairwise divergence from any user-specified sequence (CON, MRCA, COT, or existing sequence from the alignment); computing and outputting all genetic distances in column format; calculating summary statistics of diversity and divergence from pairwise distances; and graphically representing the inferred tree and plots of divergence, diversity, and distance distribution histograms. DIVEIN is available at http://indra.mullins.microbiol.washington.edu/DIVEIN.

Amino-acid co-variation in HIV-1 Gag subtype C: HLA-mediated selection pressure and compensatory dynamics.

BACKGROUND: Despite high potential for HIV-1 genetic variation, the emergence of some mutations is constrained by fitness costs, and may be associated with compensatory amino acid (AA) co-variation. To characterize the interplay between Cytotoxic T Lymphocyte (CTL)-mediated pressure and HIV-1 evolutionary pathways, we investigated AA co-variation in Gag sequences obtained from 449 South African individuals chronically infected with HIV-1 subtype C. METHODOLOGY/PRINCIPAL FINDINGS: Individuals with CTL responses biased toward Gag presented lower viral loads than individuals with under-represented Gag-specific CTL responses. Using methods that account for founder effects and HLA linkage disequilibrium, we identified 35 AA sites under Human Leukocyte Antigen (HLA)-restricted CTL selection pressure and 534 AA-to-AA interactions. Analysis of two-dimensional distances between co-varying residues revealed local stabilization mechanisms since 40% of associations involved neighboring residues. Key features of our co-variation analysis included sites with a high number of co-varying partners, such as HLA-associated sites, which had on average 55% more connections than other co-varying sites. CONCLUSIONS/SIGNIFICANCE: Clusters of co-varying AA around HLA-associated sites (especially at typically conserved sites) suggested that cooperative interactions act to preserve the local structural stability and protein function when CTL escape mutations occur. These results expose HLA-imprinted HIV-1 polymorphisms and their interlinked mutational paths in Gag that are likely due to opposite selective pressures from host CTL-mediated responses and viral fitness constraints.