BackgroundUltra-deep pyrosequencing (UDPS) allows identification of rare HIV-1 variants and minority drug resistance mutations, which are not detectable by standard sequencing.Principal FindingsHere, UDPS was used to analyze the dynamics of HIV-1 genetic variation in reverse transcriptase (RT) (amino acids 180–220) in six individuals consecutively sampled before, during and after failing 3TC and AZT containing antiretroviral treatment. Optimized UDPS protocols and bioinformatic software were developed to generate, clean and analyze the data. The data cleaning strategy reduced the error rate of UDPS to an average of 0.05%, which is lower than previously reported. Consequently, the cut-off for detection of resistance mutations was very low. A median of 16,016 (range 2,406–35,401) sequence reads were obtained per sample, which allowed detection and quantification of minority resistance mutations at amino acid position 181, 184, 188, 190, 210, 215 and 219 in RT. In four of five pre-treatment samples low levels (0.07–0.09%) of the M184I mutation were observed. Other resistance mutations, except T215A and T215I were below the detection limit. During treatment failure, M184V replaced M184I and dominated the population in combination with T215Y, while wild-type variants were rarely detected. Resistant virus disappeared rapidly after treatment interruption and was undetectable as early as after 3 months. In most patients, drug resistant variants were replaced by wild-type variants identical to those present before treatment, suggesting rebound from latent reservoirs.ConclusionsWith this highly sensitive UDPS protocol preexisting drug resistance was infrequently observed; only M184I, T215A and T215I were detected at very low levels. Similarly, drug resistant variants in plasma quickly decreased to undetectable levels after treatment interruption. The study gives important insights into the dynamics of the HIV-1 quasispecies and is of relevance for future research and clinical use of the UDPS technology.

No abstract

Increased knowledge about HIV-1 transmission dynamics in different transmission groups and geographical regions is fundamental for assessing and designing prevention efforts against HIV-1 spread. Since the first reported cases of HIV infection during the early 1980s, the HIV-1 epidemic in the Nordic countries has been dominated by HIV-1 subtype B and MSM transmission. HIV-1 pol sequences and clinical data of 51 per cent of all newly diagnosed HIV-1 infections in Sweden, Denmark, and Finland in the period 2000–2012 (N = 3,802) were analysed together with a large reference sequence dataset (N = 4,537) by trend analysis and phylogenetics. Analysis of the eight dominating subtypes and CRFs in the Nordic countries (A, B, C, D, G, CRF01_AE, CRF02_AG, and CRF06_cpx) showed that the subtype B proportion decreased while the CRF proportion increased over the study period. A majority (57 per cent) of the Nordic sequences formed transmission clusters, with evidence of mixing both geographically and between transmission groups. Detailed analyses showed multiple occasions of transmissions from MSM to heterosexuals and that active transmission clusters more often involved single than multiple Nordic countries. The strongest geographical link was between Denmark and Sweden. Finally, Denmark had a larger proportion of heterosexual domestic spread of HIV-1 subtype B (75 per cent) compared with Sweden (49 per cent) and Finland (57 per cent). We describe different HIV-1 transmission patterns between countries and transmission groups in a large geographical region. Our results may have implications for public health interventions in targeting HIV-1 transmission networks and identifying where to introduce such interventions.

BackgroundUltra-deep pyrosequencing (UDPS) is used to identify rare sequence variants. The sequence depth is influenced by several factors including the error frequency of PCR and UDPS. This study investigated the characteristics and source of errors in raw and cleaned UDPS data.ResultsUDPS of a 167-nucleotide fragment of the HIV-1 SG3Δenv plasmid was performed on the Roche/454 platform. The plasmid was diluted to one copy, PCR amplified and subjected to bidirectional UDPS on three occasions. The dataset consisted of 47,693 UDPS reads. Raw UDPS data had an average error frequency of 0.30% per nucleotide site. Most errors were insertions and deletions in homopolymeric regions. We used a cleaning strategy that removed almost all indel errors, but had little effect on substitution errors, which reduced the error frequency to 0.056% per nucleotide. In cleaned data the error frequency was similar in homopolymeric and non-homopolymeric regions, but varied considerably across sites. These site-specific error frequencies were moderately, but still significantly, correlated between runs (r = 0.15–0.65) and between forward and reverse sequencing directions within runs (r = 0.33–0.65). Furthermore, transition errors were 48-times more common than transversion errors (0.052% vs. 0.001%; p<0.0001). Collectively the results indicate that a considerable proportion of the sequencing errors that remained after data cleaning were generated during the PCR that preceded UDPS.ConclusionsA majority of the sequencing errors that remained after data cleaning were introduced by PCR prior to sequencing, which means that they will be independent of platform used for next-generation sequencing. The transition vs. transversion error bias in cleaned UDPS data will influence the detection limits of rare mutations and sequence variants.

Transmitted drug resistance (TDR) is a clinical and epidemiological problem because it may contribute to failure of antiretroviral treatment. The prevalence of TDR varies geographically, and its prevalence in Sweden during the last decade has not been reported. Plasma samples from 1,463 patients newly diagnosed with HIV-1 infection between 2003 and 2010, representing 44% of all patients diagnosed in Sweden during this period, were analyzed using the WHO 2009 list of mutations for surveillance of TDR. Maximum likelihood phylogenetic analyses were used to determine genetic subtype and to investigate the relatedness of the sequences. Eighty-two patients showed evidence of TDR, representing a prevalence of 5.6% (95% CI: 4.5%–6.9%) without any significant time trends or differences between patients infected in Sweden or abroad. Multivariable logistic regression showed that TDR was positively associated with men who have sex with men (MSM) and subtype B infection and negatively associated with CD4 cell counts. Among patients with TDR, 54 (68%) had single resistance mutations, whereas five patients had multi-drug resistant HIV-1. Phylogenetic analyses identified nine significantly supported clusters involving 29 of the patients with TDR, including 23 of 42 (55%) of the patients with TDR acquired in Sweden. One cluster contained 18 viruses with a M41L resistance mutation, which had spread among MSM in Stockholm over a period of at least 16 years (1994–2010). Another cluster, which contained the five multidrug resistant viruses, also involved MSM from Stockholm. The prevalence of TDR in Sweden 2003–2010 was lower than in many other European countries. TDR was concentrated among MSM, where clustering of TDR strains was observed, which highlights the need for continued and improved measures for targeted interventions.

Emergence of human immunodeficiency virus type 1 (HIV-1) populations that switch or broaden coreceptor usage from CCR5 to CXCR4 is intimately coupled to CD4؉ cell depletion and disease progression toward AIDS. To better understand the molecular mechanisms involved in the coreceptor switch, we determined the nucleotide sequences of 253 V1 to V3 env clones from 27 sequential HIV-1 subtype B isolates from four patients with virus populations that switch coreceptor usage. Coreceptor usage of clones from dualtropic R5X4 isolates was characterized experimentally. Sequence analysis revealed that 9% of the clones from CXCR4-using isolates had originated by recombination events between R5 and X4 viruses. The majority (73%) of the recombinants used CXCR4. Furthermore, coreceptor usage of the recombinants was determined by a small region of the envelope, including V3. This is the first report demonstrating that intrapatient recombination between viruses with distinct coreceptor usage occurs frequently. It has been proposed that X4 viruses are more easily suppressed by the immune system than R5 viruses. We hypothesize that recombination between circulating R5 viruses and X4 viruses can result in chimeric viruses with the potential to both evade the immune system and infect CXCR4-expressing cells. The broadening in cell tropism of the viral population to include CXCR4-expressing cells would gradually impair the immune system and eventually allow the X4 population to expand. In conclusion, intrapatient recombination between viruses with distinct coreceptor usage may contribute to the emergence of X4 viruses in later stages of infection.

To better understand the evolution of the viral envelope glycoproteins (Env) in HIV-1 infected individuals who progress to AIDS maintaining an exclusive CCR5-using (R5) virus population, we cloned and sequenced the env gene of longitudinally obtained primary isolates. A shift in the electrostatic potential towards an increased net positive charge was revealed in gp120 of end-stage viruses. Residues with increased positive charge were primarily localized in the gp120 variable regions, with the exception of the V3 loop. Molecular modeling indicated that the modifications clustered on the gp120 surface. Furthermore, correlations between increased Env net charge and lowered CD4(+) T cell counts, enhanced viral fitness, reduced sensitivity to entry inhibitors and augmented cell attachment were disclosed. In summary, this study suggests that R5 HIV-1 variants with increased gp120 net charge emerge in an opportunistic manner during severe immunodeficiency. Thus, we here propose a new mechanism by which HIV-1 may gain fitness.

No abstract