Tracing HIV-1 strains that imprint broadly neutralizing antibody responses.

Understanding the determinants of broadly neutralizing antibody (bNAb) evolution is crucial for the development of bNAb-based HIV vaccines1. Despite emerging information on cofactors that promote bNAb evolution in natural HIV-1 infections, in which the induction of bNAbs is genuinely rare2, information on the impact of the infecting virus strain on determining the breadth and specificity of the antibody responses to HIV-1 is lacking. Here we analyse the influence of viral antigens in shaping antibody responses in humans. We call the ability of a virus strain to induce similar antibody responses across different hosts its antibody-imprinting capacity, which from an evolutionary biology perspective corresponds to the viral heritability of the antibody responses. Analysis of 53 measured parameters of HIV-1-binding and neutralizing antibody responses in a cohort of 303 HIV-1 transmission pairs (individuals who harboured highly related HIV-1 strains and were putative direct transmission partners or members of an HIV-1 transmission chain) revealed that the effect of the infecting virus on the outcome of the bNAb response is moderate in magnitude but highly significant. We introduce the concept of bNAb-imprinting viruses and provide evidence for the existence of such viruses in a systematic screening of our cohort. The bNAb-imprinting capacity can be substantial, as indicated by a transmission pair with highly similar HIV-1 antibody responses and strong bNAb activity. Identification of viruses that have bNAb-imprinting capacities and their characterization may thus provide the potential to develop lead immunogens.

Absence of Proviral Human Immunodeficiency Virus (HIV) Type 1 Evolution in Early-Treated Individuals With HIV Switching to Dolutegravir Monotherapy During 48 Weeks.

Human immunodeficiency virus type 1 (HIV-1) infection is treated with antiretroviral therapy (ART), usually consisting of 2-3 different drugs, referred to as combination ART (cART). Our recent randomized clinical trial comparing a switch to dolutegravir monotherapy with continuation of cART in early-treated individuals demonstrated sustained virological suppression over 48 weeks. Here, we characterize the longitudinal landscape of the HIV-1 reservoir in these participants, with particular attention to potential differences between treatment groups regarding evidence of evolution as a proxy for low-level replication. Near full-length HIV-1 proviral polymerase chain reaction and next-generation sequencing was applied to longitudinal peripheral blood mononuclear cell samples to assess proviral evolution and the potential emergence of drug resistance mutations (DRMs). Neither an increase in genetic distance nor diversity over time was detected in participants of both treatment groups. Single proviral analysis showed high proportions of defective proviruses and low DRM numbers. No evidence for evolution during dolutegravir monotherapy was found in these early-treated individuals.

Prevalence of HIV-1 drug resistance mutations in proviral DNA in the Swiss HIV Cohort Study, a retrospective study from 1995 to 2018.

BACKGROUND: Genotypic resistance testing (GRT) is routinely performed upon diagnosis of HIV-1 infection or during virological failure using plasma viral RNA. An alternative source for GRT could be cellular HIV-1 DNA. OBJECTIVES: A substantial number of participants in the Swiss HIV Cohort Study (SHCS) never received GRT. We applied a method that enables access to the near full-length proviral HIV-1 genome without requiring detectable viraemia. METHODS: Nine hundred and sixty-two PBMC specimens were received. Our two-step nested PCR protocol was applied to generate two overlapping long-range amplicons of the HIV-1 genome, sequenced by next-generation sequencing (NGS) and analysed by MinVar, a pipeline to detect drug resistance mutations (DRMs). RESULTS: Six hundred and eighty-one (70.8%) of the samples were successfully amplified, sequenced and analysed by MinVar. Only partial information of the pol gene was contained in 82/681 (12%), probably due to naturally occurring deletions in the proviral sequence. All common HIV-1 subtypes were successfully sequenced. We detected at least one major DRM at high frequency (≥15%) in 331/599 (55.3%) individuals. Excluding APOBEC-signature (G-to-A mutation) DRMs, 145/599 (24.2%) individuals carried at least one major DRM. RT-inhibitor DRMs were most prevalent. The experienced time on ART was significantly longer in DRM carriers (P = 0.001) independent of inclusion or exclusion of APOBEC-signature DRMs. CONCLUSIONS: We successfully applied a reliable and efficient method to analyse near full-length HIV-1 proviral DNA and investigated DRMs in individuals with undetectable or low viraemia. Additionally, our data underscore the need for new computational tools to exclude APOBEC-related hypermutated NGS sequence reads for reporting DRMs.

Long-term experimental evolution of HIV-1 reveals effects of environment and mutational history.

An often-returning question for not only HIV-1, but also other organisms, is how predictable evolutionary paths are. The environment, mutational history, and random processes can all impact the exact evolutionary paths, but to which extent these factors contribute to the evolutionary dynamics of a particular system is an open question. Especially in a virus like HIV-1, with a large mutation rate and large population sizes, evolution is expected to be highly predictable if the impact of environment and history is low, and evolution is not neutral. We investigated the effect of environment and mutational history by analyzing sequences from a long-term evolution experiment, in which HIV-1 was passaged on 2 different cell types in 8 independent evolutionary lines and 8 derived lines, 4 of which involved a switch of the environment. The experiments lasted for 240-300 passages, corresponding to approximately 400-600 generations or almost 3 years. The sequences show signs of extensive parallel evolution-the majority of mutations that are shared between independent lines appear in both cell types, but we also find that both environment and mutational history significantly impact the evolutionary paths. We conclude that HIV-1 evolution is robust to small changes in the environment, similar to a transmission event in the absence of an immune response or drug pressure. We also find that the fitness landscape of HIV-1 is largely smooth, although we find some evidence for both positive and negative epistatic interactions between mutations.

The Interplay Between Replication Capacity of HIV-1 and Surrogate Markers of Disease.

BACKGROUND: HIV-1 replication capacity (RC) of transmitted/founder viruses may influence the further course of HIV-1 infection. METHODS: RCs of 355 whole-genome primary HIV-1 isolates derived from samples acquired during acute and recent primary HIV-1 infection (PHI) were determined using a novel high-throughput infection assay in primary cells. The RCs were used to elucidate potential factors that could be associated with RC during PHI. RESULTS: Increased RC was found to be associated with increased set point viral load (VL), and significant differences in RCs among 13 different HIV-1 subtypes were discerned. Notably, we observed an increase in RCs for primary HIV-1 isolates of HIV-1 subtype B over a 17-year period. Associations were not observed between RC and CD4 count at sample date of RC measurement, CD4 recovery after initiation of antiretroviral treatment, CD4 decline in untreated individuals, and acute retroviral syndrome severity scores. CONCLUSIONS: These findings highlight that RCs of primary HIV-1 isolates acquired during the acute and recent phase of infection are more associated with viral factors, that is set point VL, than with host factors. Furthermore, we observed a temporal increase in RC for HIV-1 subtype B viruses over a period of 17 years. CLINICAL TRIALS REGISTRATION: NCT00537966.

A Systematic Molecular Epidemiology Screen Reveals Numerous Human Immunodeficiency Virus (HIV) Type 1 Superinfections in the Swiss HIV Cohort Study.

BACKGROUND: Studying human immunodeficiency virus type 1 (HIV-1) superinfection is important to understand virus transmission, disease progression, and vaccine design. But detection remains challenging, with low sampling frequencies and insufficient longitudinal samples. METHODS: Using the Swiss HIV Cohort Study (SHCS), we developed a molecular epidemiology screening for superinfections. A phylogeny built from 22 243 HIV-1 partial polymerase sequences was used to identify potential superinfections among 4575 SHCS participants with longitudinal sequences. A subset of potential superinfections was tested by near-full-length viral genome sequencing (NFVGS) of biobanked plasma samples. RESULTS: Based on phylogenetic and distance criteria, 325 potential HIV-1 superinfections were identified and categorized by their likelihood of being detected as superinfections due to sample misidentification. NFVGS was performed for 128 potential superinfections; of these, 52 were confirmed by NFVGS, 15 were not confirmed, and for 61 sampling did not allow confirming or rejecting superinfection because the sequenced samples did not include the relevant time points causing the superinfection signal in the original screen. Thus, NFVGS could support 52 of 67 adequately sampled potential superinfections. CONCLUSIONS: This cohort-based molecular approach identified, to our knowledge, the largest population of confirmed superinfections, showing that, while rare with a prevalence of 1%-7%, superinfections are not negligible events.

Prevalence of integrase strand transfer inhibitor resistance mutations in antiretroviral-naive HIV-1-infected individuals in Cameroon.

OBJECTIVES: In Cameroon, the integrase (IN) strand transfer inhibitor (INSTI) dolutegravir was recently introduced for the treatment of HIV-1 infection. Since pretreatment HIV-1 drug resistance can jeopardize the success of ART, and considering the high heterogeneity of circulating HIV-1 subtypes in Cameroon, we investigated the prevalence of pretreatment HIV-1 resistance to INSTIs. METHODS: Fingerprick dried blood spot samples were collected from 339 newly diagnosed HIV-1-infected individuals between 2015 and 2016 in four hospitals in Cameroon. Universal primers were designed to amplify the HIV-1 IN region from amino acid 1 to 276. Amplicons were sequenced with Illumina next-generation sequencing and analysed with the Polymorphism Analysis Sequencing (PASeq) platform, using the Stanford HIV Drug Resistance Database to interpret HIV-1 drug resistance mutations (DRMs). RESULTS: The amplification/sequencing success rate was 75.2% with 255/339 sequences obtained. Applying a cut-off of 1%, major DRMs to INSTIs were detected in 13 (5.1%) individuals, but only 1 individual harboured an INSTI DRM (E92G) at a nucleotide frequency ≥15%. However, 140/255 (54.9%) individuals harboured polymorphic accessory INSTI DRMs, mainly at high frequencies. In line with that observation, HIV-1 subtype diversity among individuals was high. CONCLUSIONS: Pretreatment HIV-1 resistance to INSTIs was low in the study sites, which supports the use of INSTIs in Cameroon. Nevertheless, further studies are necessary to assess the impact of polymorphic accessory INSTI DRMs on INSTI-based ART regimens.

Parallel Evolution of HIV-1 in a Long-Term Experiment.

One of the most intriguing puzzles in biology is the degree to which evolution is repeatable. The repeatability of evolution, or parallel evolution, has been studied in a variety of model systems, but has rarely been investigated with clinically relevant viruses. To investigate parallel evolution of HIV-1, we passaged two replicate HIV-1 populations for almost 1 year in each of two human T-cell lines. For each of the four evolution lines, we determined the genetic composition of the viral population at nine time points by deep sequencing the entire genome. Mutations that were carried by the majority of the viral population accumulated continuously over 1 year in each evolution line. Many majority mutations appeared in more than one evolution line, that is, our experiments showed an extreme degree of parallel evolution. In one of the evolution lines, 62% of the majority mutations also occur in another line. The parallelism impairs our ability to reconstruct the evolutionary history by phylogenetic methods. We show that one can infer the correct phylogenetic topology by including minority mutations in our analysis. We also find that mutation diversity at the beginning of the experiment is predictive of the frequency of majority mutations at the end of the experiment.

Viral Diversity Based on Next-Generation Sequencing of HIV-1 Provides Precise Estimates of Infection Recency and Time Since Infection.

BACKGROUND: Human immunodeficiency virus type 1 (HIV-1) genetic diversity increases over the course of infection and can be used to infer the time since infection and, consequently, infection recency, which are crucial for HIV-1 surveillance and the understanding of viral pathogenesis. METHODS: We considered 313 HIV-infected individuals for whom reliable estimates of infection dates and next-generation sequencing (NGS)-derived nucleotide frequency data were available. Fractions of ambiguous nucleotides, obtained by population sequencing, were available for 207 samples. We assessed whether the average pairwise diversity calculated using NGS sequences provided a more exact prediction of the time since infection and classification of infection recency (<1 year after infection), compared with the fraction of ambiguous nucleotides. RESULTS: NGS-derived average pairwise diversity classified an infection as recent with a sensitivity of 88% and a specificity of 85%. When considering only the 207 samples for which fractions of ambiguous nucleotides were available, the NGS-derived average pairwise diversity exhibited a higher sensitivity (90% vs 78%) and specificity (95% vs 67%) than the fraction of ambiguous nucleotides. Additionally, the average pairwise diversity could be used to estimate the time since infection with a mean absolute error of 0.84 years, compared with 1.03 years for the fraction of ambiguous nucleotides. CONCLUSIONS: Viral diversity based on NGS data is more precise than that based on population sequencing in its ability to predict infection recency and provides an estimated time since infection that has a mean absolute error of <1 year.

Full-length haplotype reconstruction to infer the structure of heterogeneous virus populations.

Next-generation sequencing (NGS) technologies enable new insights into the diversity of virus populations within their hosts. Diversity estimation is currently restricted to single-nucleotide variants or to local fragments of no more than a few hundred nucleotides defined by the length of sequence reads. To study complex heterogeneous virus populations comprehensively, novel methods are required that allow for complete reconstruction of the individual viral haplotypes. Here, we show that assembly of whole viral genomes of ∼8600 nucleotides length is feasible from mixtures of heterogeneous HIV-1 strains derived from defined combinations of cloned virus strains and from clinical samples of an HIV-1 superinfected individual. Haplotype reconstruction was achieved using optimized experimental protocols and computational methods for amplification, sequencing and assembly. We comparatively assessed the performance of the three NGS platforms 454 Life Sciences/Roche, Illumina and Pacific Biosciences for this task. Our results prove and delineate the feasibility of NGS-based full-length viral haplotype reconstruction and provide new tools for studying evolution and pathogenesis of viruses.

Origin of minority drug-resistant HIV-1 variants in primary HIV-1 infection.

BACKGROUND: Drug-resistant human immunodeficiency virus type 1 (HIV-1) minority variants (MVs) are present in some antiretroviral therapy (ART)-naive patients. They may result from de novo mutagenesis or transmission. To date, the latter has not been proven. METHODS: MVs were quantified by allele-specific polymerase chain reaction in 204 acute or recent seroconverters from the Zurich Primary HIV Infection study and 382 ART-naive, chronically infected patients. Phylogenetic analyses identified transmission clusters. RESULTS: Three lines of evidence were observed in support of transmission of MVs. First, potential transmitters were identified for 12 of 16 acute or recent seroconverters harboring M184V MVs. These variants were also detected in plasma and/or peripheral blood mononuclear cells at the estimated time of transmission in 3 of 4 potential transmitters who experienced virological failure accompanied by the selection of the M184V mutation before transmission. Second, prevalence between MVs harboring the frequent mutation M184V and the particularly uncommon integrase mutation N155H differed highly significantly in acute or recent seroconverters (8.2% vs 0.5%; P < .001). Third, the prevalence of less-fit M184V MVs is significantly higher in acutely or recently than in chronically HIV-1-infected patients (8.2% vs 2.5%; P = .004). CONCLUSIONS: Drug-resistant HIV-1 MVs can be transmitted. To what extent the origin-transmission vs sporadic appearance-of these variants determines their impact on ART needs to be further explored.

Delay of HIV-1 rebound after cessation of antiretroviral therapy through passive transfer of human neutralizing antibodies.

To determine the protective potential of the humoral immune response against HIV-1 in vivo we evaluated the potency of three neutralizing antibodies (2G12, 2F5 and 4E10) in suppressing viral rebound in six acutely and eight chronically HIV-1-infected individuals undergoing interruption of antiretroviral treatment (ART). Only two of eight chronically infected individuals showed evidence of a delay in viral rebound during the passive immunization. Rebound in antibody-treated acutely infected individuals upon cessation of ART was substantially later than in a control group of 12 individuals with acute infection. Escape mutant analysis showed that the activity of 2G12 was crucial for the in vivo effect of the neutralizing antibody cocktail. By providing further direct evidence of the potency, breadth and titers of neutralizing antibodies that are required for in vivo activity, these data underline both the potential and the limits of humoral immunity in controlling HIV-1 infection.

Characterization of human immunodeficiency virus type 1 (HIV-1) diversity and tropism in 145 patients with primary HIV-1 infection.

BACKGROUND: In the context of sexual transmission of human immunodeficiency virus type 1 (HIV-1), current findings suggest that the mucosal barrier is the major site of viral selection, transforming the complex inoculum to a small, homogeneous founder virus population. We analyzed HIV-1 transmission in relation to viral and host characteristics within the Zurich primary HIV-1 infection study. METHODS: Clonal HIV-1 envelope sequences (on average 16 clones/patient) were isolated from the first available plasma samples during the early phase of infection from 145 patients with primary HIV-1 infection. Phylogenetic and tropism analyses were performed. Differences of viral diversities were investigated in association with several parameters potentially influencing HIV-1 transmission, eg, concomitant sexually transmitted infections (STIs) and mode of transmission. RESULTS: Median viral diversity within env C2-V3-C3 region was 0.39% (range 0.04%-3.23%). Viral diversity did not correlate with viral load, but it was slightly correlated with the duration of infection. Neither transmission mode, gender, nor STI predicted transmission of more heterogeneous founder virus populations that were found in 16 of 145 patients (11%; diversity >1%). Only 2 patients (1.4%) were assuredly infected with CXCR4-tropic HIV-1 within a R5/X4-tropic--mixed population, as revealed and confirmed using several genotypic prediction algorithms and phenotypic assays. CONCLUSIONS: Our findings suggest that transmission of multiple HIV-1 variants might be a complex process that is not dependent on mucosal factors alone. CXCR4-tropic viruses can be sexually transmitted in rare instances, but their clinical relevance remains to be determined.

Efficient suppression of minority drug-resistant HIV type 1 (HIV-1) variants present at primary HIV-1 infection by ritonavir-boosted protease inhibitor-containing antiretroviral therapy.

BACKGROUND: Selection of preexisting minority variants of drug-resistant human immunodeficiency virus type 1 (HIV-1) can lead to virological failure in patients who receive antiretroviral therapy (ART) with low genetic resistance barriers. We studied treatment response and dynamics of minority variants during the first weeks of ART containing a ritonavir-boosted protease inhibitor (PI) and 2 nucleoside reverse-transcriptase inhibitors (NRTIs), which is a regimen with a high genetic resistance barrier. METHODS: Plasma samples obtained prior to initiation of ART from 109 patients with primary HIV infection and samples obtained during viral decay during early ART from 17 of these 109 patients were tested by allele-specific polymerase chain reaction for K103N and M184V variants. RESULTS: K103N and/or M184V mutations were detected in 15 (13.8%) of 109 patients prior to ART as minority variants. No selection of these variants was observed within the first weeks of ART in 7 of 15 patients with preexisting drug resistance mutations, nor was any selection observed in 10 patients without preexisting drug resistance mutations. Most patients received ART immediately after diagnosis of HIV-1 infection, showed a rapid decrease in viral load, and experienced sufficient suppression of viremia for 48 months. CONCLUSIONS: Minority variants, in particular viruses harboring the M184V mutation, were efficiently suppressed in patients with acute infection who received a ritonavir-boosted PI and 2 NRTIs (most regimens included lamivudine). Under this high genetic resistance barrier regimen, the M184V was not further selected.

A Novel High Throughput, Parallel Infection Assay for Determining the Replication Capacities of 346 Primary HIV-1 Isolates of the Zurich Primary HIV-1 Infection Study in Primary Cells.

HIV-1 replication capacity is an important characteristic to understand the replication competence of single variants or virus populations. It can further aid in the understanding of HIV-1 pathogenicity, disease progression, and drug resistance mutations. To effectively study RC, many assays have been established. However, there is still demand for a high throughput replication capacity assay using primary cells which is robust and reproducible. In this study, we established such an assay and validated it using 346 primary HIV-1 isolates from patients enrolled in the Zurich Primary HIV Infection study (ZPHI) and two control viruses, HIV-1 JR-CSFWT and HIV-1 JR-CSFK65R_M184V. Replication capacity was determined by measuring the viral growth on PBMCs over 10 days by longitudinally transferring cell culture supernatant to TZM-bl reporter cells. By utilizing the TZM-bl luciferase reporter assay, we determined replication capacity by measuring viral infectivity. The simplicity of the experimental setup allowed for all 346 primary HIV-1 isolates to be replicated at one time. Although the infectious input dose for each virus was normalized, a broad range of replication capacity values over 4 logs was observed. The approach was confirmed by two repeated experiments and we demonstrated that the reproducibility of the replication capacity values is statistically comparable between the two separate experiments. In summary, these results endorse our high throughput replication capacity assay as reproducible and robust and can be utilized for large scale HIV-1 replication capacity experiments in primary cells.

HIV-1 integration sites in CD4+ T cells during primary, chronic, and late presentation of HIV-1 infection.

HIV-1 is capable of integrating its genome into that of its host cell. We examined the influence of the activation state of CD4+ T cells, the effect of antiretroviral therapy (ART), and the clinical stage of HIV-1 infection on HIV-1 integration site features and selection. HIV-1 integration sites were sequenced from longitudinally sampled resting and activated CD4+ T cells from 12 HIV-1-infected individuals. In total, 589 unique HIV-1 integration sites were analyzed: 147, 391, and 51 during primary, chronic, and late presentation of HIV-1 infection, respectively. As early as during primary HIV-1 infection and independent of the activation state of CD4+ T cells collected on and off ART, HIV-1 integration sites were preferentially detected in recurrent integration genes, genes associated with clonal expansion of latently HIV-1-infected CD4+ T cells, cancer-related genes, and highly expressed genes. The preference for cancer-related genes was more pronounced at late stages of HIV-1 infection. Host genomic features of HIV-1 integration site selection remained stable during HIV-1 infection in both resting and activated CD4+ T cells. In summary, characteristic HIV-1 integration site features are preestablished as early as during primary HIV-1 infection and are found in both resting and activated CD4+ T cells.

In vivo efficacy of human immunodeficiency virus neutralizing antibodies: estimates for protective titers.

The definition of plasma neutralizing antibody titers capable of controlling human immunodeficiency virus (HIV) infection in vivo is considered a critical step in vaccine development. Here we provide estimates for effective neutralization titers by assessing samples from a recent passive immunization trial with the neutralizing monoclonal antibodies (MAbs) 2G12, 2F5, and 4E10 using an analytic strategy that dissects the contributions of these MAbs to the total neutralization activity in patient plasma. Assessment of neutralization activities for six responding patients with partial or complete control of viremia during the MAb treatment and for the eight nonresponding patients revealed a significant difference between these groups: Among responders, MAb-mediated activity exceeded the autologous neutralization response by 1 to 2 log units (median difference, 43.3-fold), while in the nonresponder group, the autologous activity prevailed (median difference, 0.63-fold). In order to reach a 50% proportion of the responders in our study cohort, MAb neutralizing titers higher than 1:200 were required based on this analysis. The disease stage appears to have a significant impact on the quantities needed, since titers above 1:1,000 were needed to reach the same effect in chronic infection. Although our analysis is based on very small sample numbers and thus cannot be conclusive, our data provide a first estimate on how in vitro-measured neutralizing antibody activity can relate to in vivo efficacy in controlling HIV infection and may therefore provide valuable information for vaccine development. Interestingly, lower neutralizing antibody levels showed an effect in acute compared to chronic infection, suggesting that in early disease stages, therapeutic vaccination may show promise. Equally, this raises hopes that a preventive vaccine could become effective at comparatively lower neutralizing antibody titers.

MinVar: A rapid and versatile tool for HIV-1 drug resistance genotyping by deep sequencing.

Genotypic monitoring of drug-resistance mutations (DRMs) in HIV-1 infected individuals is strongly recommended to guide selection of the initial antiretroviral therapy (ART) and changes of drug regimens. Traditionally, mutations conferring drug resistance are detected by population sequencing of the reverse transcribed viral RNA encoding the HIV-1 enzymes target by ART, followed by manual analysis and interpretation of Sanger sequencing traces. This process is labor intensive, relies on subjective interpretation from the operator, and offers limited sensitivity as only mutations above 20% frequency can be reliably detected. Here we present MinVar, a pipeline for the analysis of deep sequencing data, which allows reliable and automated detection of DRMs down to 5%. We evaluated MinVar with data from amplicon sequencing of defined mixtures of molecular virus clones with known DRM and plasma samples of viremic HIV-1 infected individuals and we compared it to VirVarSeq, another virus variant detection tool exclusively working on Illumina deep sequencing data. MinVar was designed to be compatible with a diverse range of sequencing platforms and allows the detection of DRMs and insertions/deletions from deep sequencing data without the need to perform additional bioinformatics analysis, a prerequisite to a widespread implementation of HIV-1 genotyping using deep sequencing in routine diagnostic settings.

Next-generation sequencing of HIV-1 RNA genomes: determination of error rates and minimizing artificial recombination.

Next-generation sequencing (NGS) is a valuable tool for the detection and quantification of HIV-1 variants in vivo. However, these technologies require detailed characterization and control of artificially induced errors to be applicable for accurate haplotype reconstruction. To investigate the occurrence of substitutions, insertions, and deletions at the individual steps of RT-PCR and NGS, 454 pyrosequencing was performed on amplified and non-amplified HIV-1 genomes. Artificial recombination was explored by mixing five different HIV-1 clonal strains (5-virus-mix) and applying different RT-PCR conditions followed by 454 pyrosequencing. Error rates ranged from 0.04-0.66% and were similar in amplified and non-amplified samples. Discrepancies were observed between forward and reverse reads, indicating that most errors were introduced during the pyrosequencing step. Using the 5-virus-mix, non-optimized, standard RT-PCR conditions introduced artificial recombinants in a fraction of at least 30% of the reads that subsequently led to an underestimation of true haplotype frequencies. We minimized the fraction of recombinants down to 0.9-2.6% by optimized, artifact-reducing RT-PCR conditions. This approach enabled correct haplotype reconstruction and frequency estimations consistent with reference data obtained by single genome amplification. RT-PCR conditions are crucial for correct frequency estimation and analysis of haplotypes in heterogeneous virus populations. We developed an RT-PCR procedure to generate NGS data useful for reliable haplotype reconstruction and quantification.