Abstract
Free full text
![Logo of aac](/http/europepmc.org/corehtml/pmc/pmcgifs/logo-aac.png)
Sensitive Deep-Sequencing-Based HIV-1 Genotyping Assay To Simultaneously Determine Susceptibility to Protease, Reverse Transcriptase, Integrase, and Maturation Inhibitors, as Well as HIV-1 Coreceptor Tropism
Associated Data
Abstract
With 29 individual antiretroviral drugs available from six classes that are approved for the treatment of HIV-1 infection, a combination of different phenotypic and genotypic tests is currently needed to monitor HIV-infected individuals. In this study, we developed a novel HIV-1 genotypic assay based on deep sequencing (DeepGen HIV) to simultaneously assess HIV-1 susceptibilities to all drugs targeting the three viral enzymes and to predict HIV-1 coreceptor tropism. Patient-derived gag-p2/NCp7/p1/p6/pol-PR/RT/IN- and env-C2V3 PCR products were sequenced using the Ion Torrent Personal Genome Machine. Reads spanning the 3′ end of the Gag, protease (PR), reverse transcriptase (RT), integrase (IN), and V3 regions were extracted, truncated, translated, and assembled for genotype and HIV-1 coreceptor tropism determination. DeepGen HIV consistently detected both minority drug-resistant viruses and non-R5 HIV-1 variants from clinical specimens with viral loads of ≥1,000 copies/ml and from B and non-B subtypes. Additional mutations associated with resistance to PR, RT, and IN inhibitors, previously undetected by standard (Sanger) population sequencing, were reliably identified at frequencies as low as 1%. DeepGen HIV results correlated with phenotypic (original Trofile, 92%; enhanced-sensitivity Trofile assay [ESTA], 80%; TROCAI, 81%; and VeriTrop, 80%) and genotypic (population sequencing/Geno2Pheno with a 10% false-positive rate [FPR], 84%) HIV-1 tropism test results. DeepGen HIV (83%) and Trofile (85%) showed similar concordances with the clinical response following an 8-day course of maraviroc monotherapy (MCT). In summary, this novel all-inclusive HIV-1 genotypic and coreceptor tropism assay, based on deep sequencing of the PR, RT, IN, and V3 regions, permits simultaneous multiplex detection of low-level drug-resistant and/or non-R5 viruses in up to 96 clinical samples. This comprehensive test, the first of its class, will be instrumental in the development of new antiretroviral drugs and, more importantly, will aid in the treatment and management of HIV-infected individuals.
INTRODUCTION
According to the Joint United Nations Programme on HIV/AIDS (UNAIDS), at the end of 2012, close to 10 million of the approximately 34 million people living with HIV had access to antiretroviral therapy, with the goal being to reach 15 million people with HIV treatment by 2015 (1). Not only has this broader access to antiretroviral drugs led to considerable reductions in morbidity and mortality (2, 3), but unfortunately, it has increased the risk of virologic failure due to the emergence (3, 4) and potential transmission (5) of drug-resistant viruses. Therefore, detecting and quantifying drug resistance, and in the case of CCR5 receptor antagonists, determining HIV-1 coreceptor tropism, have become the standard of care prior to designing new antiretroviral regimens (3, 6,–8).
To date, 29 individual antiretroviral drugs from six drug classes have been approved by the U.S. Food and Drug Administration (FDA) to be used in the treatment of HIV-1 infection, including protease inhibitors (PI), nucleoside/nucleotide reverse transcriptase inhibitors (NRTI), nonnucleoside reverse transcriptase inhibitors (NNRTI), integrase inhibitors (INI), fusion inhibitors (FI), and entry inhibitors (EI). HIV-1 resistance to PI, NRTI, NNRTI, and INI can be determined using (i) indirect methods based on the detection of specific amino acid substitutions (due to underlying nucleotide mutations) in the respective coding regions previously associated with resistance to specific antiretroviral drugs (i.e., genotyping) (6, 9, 10), (ii) more direct methods that test the ability of a patient-derived virus to replicate in the presence of antiretroviral drugs in a cell-based assay (i.e., phenotyping) (11,–13), or (iii) a combination of the two approaches that takes advantage of a large database in order to infer the level of HIV-1 drug resistance based on genotyping and its relationship with matched phenotypic data (14). Similarly, since treatment with CCR5 antagonists requires prior knowledge of the HIV-1 coreceptor tropism in the patient, i.e., CCR5- or CXCR4-tropic viruses (R5 and X4, respectively), dual tropism (R5/X4), or a mixture of both R5 and X4 viruses (7, 15), a multitude of phenotypic and genotypic approaches to determine HIV-1 coreceptor tropism have been developed (8, 16, 17). Phenotypic assays to determine HIV-1 drug resistance or tropism usually involve the generation of patient-derived pol (11,–13) or env recombinant (18,–20) viruses, respectively, to quantify their ability to infect susceptible cell lines expressing the appropriate HIV-1 receptors and coreceptors; in the case of HIV-1 tropism, these measures may also be based on the quantification of cell-to-cell fusion events (21,–23). Conversely, genotypic HIV-1 tropism tests (8, 15, 16, 24,–26) take advantage of the properties of specific regions in the env gene as determinants of CCR5 or CXCR4 tropism, mainly the V3 region of gp120, and their interpretations are based on a series of bioinformatics methods to infer the ability of HIV-1 to use either or both coreceptors to enter host cells (27,–30).
As expected, phenotypic (experimental) and genotypic (computational) approaches to determining HIV-1 drug resistance or HIV-1 coreceptor tropism have some disadvantages, including the longer turnaround times and higher cost of the phenotypic assays or the intrinsic predictive nature of the genotypic tests. Particular emphasis has been made on the limited sensitivities of genotypic HIV-1 tropism assays to detect minor non-R5 variants (16, 31), and to a lesser extent on the ability of genotypic HIV-1 drug resistance tests to detect minority drug-resistant variants (32,–34). In the case of HIV-1 drug resistance, the vast amount of information accumulated during the last 2 decades by correlating mutations with phenotypic data has led to the almost exclusive use of genotypic antiretroviral testing based on population (Sanger) sequencing to manage patients infected with HIV-1 (2, 35). In contrast, although several studies have shown significant concordance and similar predictive values (36,–40), genotypic HIV-1 tropism assays based on population sequencing seem to be less sensitive and specific than phenotypic assays (8, 16, 17, 41). Thus, a cell-based assay (Trofile; Monogram Biosciences) (19, 42) is currently the standard method in the United States for determining HIV-1 coreceptor tropism, while genotypic HIV-1 tropism tests are largely used in Europe (16, 31).
To date, all current commercial genotypic HIV-1 drug resistance assays are based on population sequencing (10, 43, 44), which can detect only minority variants that are present in >20% of the viral population (44,–48). However, and although this is still uncertain, drug-resistant HIV-1 minority variants (i.e., those present in as low as 1% of the viral population) have been suggested to be clinically relevant, as they have a high chance of selection under antiretroviral drug pressure conditions (49,–57). For that reason, a series of ultrasensitive assays have been developed to detect drug-resistant HIV-1 minority variants, e.g., allele-specific PCR (49, 58), oligonucleotide ligation assays (33, 59), and deep (next-generation) sequencing (60,–62). On the other hand, as described above, the adoption of genotypic HIV-1 tropism assays in the clinical setting has been hampered by the limited sensitivities of the population-based sequencing assays to detect minor non-R5 variants. Therefore, more sensitive genotypic HIV-1 tropism assays based on deep sequencing have been developed to detect non-R5 variants present at frequencies of <20% of the population, and these have been shown to correlate well with both phenotypic assays (36, 63,–67) and the virological response to CCR5-receptor antagonists, such as maraviroc (Selzentry/Celsentri, Pfizer, NY) (36, 63, 66). Nevertheless, a combination of at least two different genotypic assays is still needed to assess the susceptibility of a patient-derived HIV-1 infection to all FDA-approved antiretroviral drugs, including CCR5 antagonists. Therefore, in this study, we have developed, characterized, and validated a novel HIV-1 genotyping assay based on deep sequencing to simplify the monitoring of patients infected with HIV-1. This all-inclusive sensitive methodology accurately provides drug resistance information for all protease, reverse transcriptase, integrase, and maturation inhibitors, as well as HIV-1 coreceptor tropism, in a single, more efficient, rapid, and affordable clinical assay.
(This research was presented in part at the International HIV & Hepatitis Virus Drug Resistance Workshop and Curative Strategies, Toronto, Ontario, Canada, 4 to 8 June 2013.)
MATERIALS AND METHODS
Viruses and plasmids.
The following viruses were obtained from the AIDS Research and Reference Reagent Program, Division of AIDS, NIAID, NIH: HIV-1A-92RW009, HIV-1A-93RW020, HIV-1A-92UG029, HIV-1B-92BR014, HIV-1B-92TH593, HIV-1B-US714, HIV-1B-92US727, HIV-1B-92US076, HIV-1C-92BR025, HIV-1D-94UG108, HIV-1D-92UG038, HIV-1D-93UG065, HIV-1F-93BR029, HIV-1F-93BR020, HIV-1G-RU570, HIV-1G-RU132, HIV-1AE-CMU02, HIV-1AE-CMU06, HIV-1AE-92TH021, HIV-1BF-93BR029, and HIV-2CBL-20. Other viruses were from Eric J. Arts' laboratory at Case Western Reserve University (CWRU), Cleveland, OH: HIV-1A-V115, HIV-1A-V120, HIV-1C-C18, HIV-1C-C20, HIV-1C-C21, HIV-1C-C22, HIV-1D-V89, HIV-1D-V122, HIV-1D-V126, HIV-1F-VI820, HIV-1F-V164, HIV-1F-CA16, and HIV-1F-CA20. Aliquots of additional RNA or DNA viruses were obtained from the Molecular Diagnostics or Medical Microbiology laboratories at University Hospitals Case Medical Center (UHCMC), Cleveland, OH (BK virus [BKV], cytomegalovirus [CMV], herpes simplex virus 1 and 2 [HSV-1 and HSV-2], and varicella zoster virus [VZV]) or the Division of Infectious Diseases, School of Medicine at CWRU (hepatitis B virus [HBV], hepatitis C virus [HCV], and Epstein-Barr virus [EBV]). Plasmids containing patient-derived HIV-1 gag-p2/NCp7/p1/p6/pol-PR/RT/IN-coding sequences from multidrug-resistant viruses, i.e., 08-180 and 08-194, have been described (68), along with a pNL4-3 plasmid containing the env gene from the R5 HIV-1YU2 virus (69). Five plasmid mixtures containing drug resistance mutations in the HIV-1 pol gene, i.e., K65R (5%) plus wild type (95%), K103N (5%) plus wild type (95%), K101E (5%) plus E138K (5%) plus wild type (90%), K101E plus E138K (10%) plus wild type (90%), and M184V (RT) plus E92Q (IN) (10%) plus wild type (90%), were obtained from Gilead Sciences, Inc. (Foster City, CA).
Clinical samples.
Blood plasma samples for the characterization and verification of the novel HIV-1 genotypic and coreceptor tropism assay were obtained during routine patient monitoring from a well-characterized cohort of HIV-infected individuals at the AIDS Clinical Trials Unit (ACTU) at CWRU/UHCMC, with the understanding and written consent of each participant. RNA specimens, derived from plasma samples collected from HIV-infected individuals enrolled in the (i) maraviroc expanded-access program in Europe or (ii) the ALLEGRO trial, were obtained from the Hospital Carlos III (Madrid, Spain) (40). Written informed consent was obtained from the patients before participation in the study, as previously described (40, 70). HIV-1 coreceptor tropism was determined at baseline using two phenotypic assays, i.e., the original version of the Trofile (19) and VeriTrop (23), and by population sequencing analyzed with Geno2Pheno (27), with false-positive rates (FPR) (predicted frequency of classifying an R5 sequence as a non-R5 virus) based on optimized cutoffs associated with the analysis of clinical data from MOTIVATE (2.5% and 5.75%) (63). Finally, plasma samples were obtained from HIV-infected individuals at the Infectious Diseases Unit Virgen del Rocio University Hospital (Seville, Spain) participating in a study to evaluate the use of an 8-day maraviroc monotherapy clinical test (MCT) (71, 72). Patients provided written informed consent, and the ethics committee of the hospital approved the study (72). HIV-1 coreceptor tropism in these samples was determined at baseline using two different phenotypic assays, i.e., the enhanced-sensitivity Trofile assay (ESTA) (42) and TROCAI (73), and by population sequencing analyzed with Geno2Pheno (27) with an FPR of 10%, according to the recommendations from the European Consensus Group on the clinical management of HIV-1 tropism testing, as described on the Geno2Pheno website (http://coreceptor.bioinf.mpi-inf.mpg.de/index.php).
Reverse transcription-PCR amplification of gag-p2/NCp7/p1/p6/pol-PR/RT/IN- and env-C2V3-coding regions.
Plasma viral RNA was purified from pelleted virus particles by centrifuging 1 ml of plasma at 18,000 × g for 60 min at 4°C, removing 860 μl of cell-free supernatant, resuspending the pellet in the remaining 140 μl, and finally extracting viral RNA using the QIAamp viral RNA minikit (Qiagen; Valencia, CA). Viral RNA was reverse transcribed using AccuScript high-fidelity reverse transcriptase (Stratagene Agilent; Santa Clara, CA) and the corresponding antisense external primers in a 20-μl reaction mixture containing 1 mM deoxynucleoside triphosphates (dNTPs), 10 mM dithiothreitol (DTT), and 10 units of RNase inhibitor. The HIV-1 genomic region encoding the Gag proteins p2, p7, p1, and p6 and the protease, reverse transcriptase, and integrase enzymes was amplified as two overlapping fragments (1,657 nucleotides [nt] and 2,002 nt corresponding to the p2-5′half RT and 3′half RT-INT, respectively) using a series of external and nested primers with defined cycling conditions (13). External PCRs were carried out in a 50-μl mixture containing 0.2 mM dNTPs, 1 mM MgCl2, and 2.5 units of Pfu Turbo DNA polymerase (Stratagene). Nested PCRs were carried out in a 50-μl mixture containing 0.2 mM dNTPs, 0.3 units of Pfu Turbo DNA polymerase, and 1.9 units of Taq polymerase (Denville Scientific; Metuchen, NJ). A fragment corresponding to the C2V3 region (480 nt) of the surface glycoprotein (gp120) in the envelope gene was amplified using a series of external and nested primers with defined cycling conditions, as previously described (74).
Population (Sanger) sequencing analysis.
PCR products corresponding to the gag-p2/NCp7/p1/p6/pol-PR/RT/IN- and env-C2V3-coding regions of HIV-1 were purified with the QIAquick PCR purification kit (Qiagen) and sequenced (Sanger, population, or global sequence) using AP Biotech DYEnamic ET Terminator cycle with Thermosequenase II (Davis Sequencing LLC, Davis, CA). Nucleotide sequences were analyzed using DNAStar Lasergene Software Suite version 10.0.1 (Madison, WI).
Deep sequencing of gag-p2/NCp7/p1/p6/pol-PR/RT/IN- and env-C2V3-coding regions.
The three PCR products corresponding to the gag-p2/NCp7/p1/p6/pol-PR/RT/IN- (1,657-nt and 2,002-nt fragments) and env-C2V3- (480-nt fragment) coding regions of HIV-1 were purified (Agencourt AMPure XP; Beckman Coulter) and quantified (2100 Bioanalyzer DNA 7500; Agilent Technologies) prior to using the Ion Xpress fragment library kit (Life Technologies, Carlsbad CA) to construct a multiplexed library for shotgun sequencing on the Ion Personal Genome Machine (PGM) (Life Technologies) (Fig. 1). Briefly, a mixture of all three purified DNA amplicons (33 ng each) was randomly fragmented and the blunt ends repaired using the Ion Shear Plus reagent (Life Technologies), followed by DNA purification (Agencourt AMPure XP). The P1 adapter (5′-CCA CTA CGC CTC CGC TTT CCT CTC TAT GGG CAG TCG GTG AT, 5′-ATC ACC GAC TGC CCA TAG AGA GGA AAG CGG AGG CGT AGT GG*T*T [asterisks indicate phosphorothioate bonds]) and one of 96 barcodes were ligated to the repaired fragment ends prior to DNA purification (Agencourt AMPure XP). DNA fragments were then selected by size (i.e., 300 bp; Pippin Prep; Life Technologies) and each bar-coded library, i.e., a mixture of all three amplicons per sample, was purified (Agencourt AMPure XP) and normalized using the Ion Library Equalizer kit (Life Technologies). All bar-coded DNA libraries corresponding to patient-derived amplicons plus the HIV-1NL4-3 control were pooled in equimolar concentrations and the templates prepared and enriched for sequencing on the Ion Sphere particles (ISPs) using the Ion OneTouch 200 template kit version 2 (Life Technologies) in the Ion OneTouch 2 system (Life Technologies). Templated ISPs were quantified (Qubit 2.0; Life Technologies) and loaded into an Ion 318 Chip (Life Technologies) to be sequenced on the Ion PGM using the Ion PGM sequencing 200 kit version 2 (Life Technologies). Following a 4-h-and-20-min sequencing run, signal processing and base calling were performed with the Torrent Analysis Suite version 3.4.2.
![An external file that holds a picture, illustration, etc.
Object name is zac0041427490001.jpg An external file that holds a picture, illustration, etc.
Object name is zac0041427490001.jpg](/http/europepmc.org/articles/PMC4023761/bin/zac0041427490001.jpg)
Overview of the protocol for the novel HIV-1 genotyping and coreceptor tropism assay (DeepGen HIV). Three PCR products corresponding to the gag-p2/NCp7/p1/p6/pol-PR/RT/IN- (1,657 bp and 2,002 bp) and env-C2V3- (480 bp) coding regions of HIV-1 were used to construct a multiplexed library for shotgun sequencing on the Ion PGM. Signal processing and base calling were performed with the Torrent analysis suite version 3.4.2 and sequences analyzed using DEEPGEN Software Tool Suite. The HIVdb Program Genotypic Resistance Interpretation Algorithm from the Stanford University HIV Drug Resistance Database (see http://hivdb.stanford.edu) and Geno2Pheno (27) were used to infer the levels of susceptibility to PI, RTI, and INI and for HIV-1 coreceptor tropism determination, respectively.
Read mapping, variant calling, and phylogenetic analysis.
As part of the novel HIV-1 genotypic and coreceptor tropism assay, we developed the DeepGen Software Tool Suite for the processing of HIV-1 deep sequencing data and HIV-1 drug resistance determination. DeepGen uses two main tools: Viral Read Mapper and Variant Caller.
(i) Viral Read Mapper.
To minimize the amount of data loss during mapping due to the high HIV-1 sequence variability and to allow for interpatient indel variation across the gag-p2/NCp7/p1/p6/pol-PR/RT/IN- and env-C2V3-coding regions, sample-specific reference sequences were constructed for each of these two genomic regions, i.e., positions 1807 to 5096 and 6900 to 7400 in the HXB2 reference strain (GenBank accession no. K03455), respectively. The mapping of reads from each sample/region occurred in three stages. First, a guide template for mapping was selected from the Los Alamos HIV Sequence Database (see http://www.hiv.lanl.gov/content/sequence/HIV/mainpage.html) by comparing 100 randomly selected reads to the corresponding region within all full-length sequences present within the HIV Sequence Database. This comparison was performed using a k-mer approach that rapidly identifies regions of similarity between any two sequences to select a guide sequence for mapping with minimal divergence from the read data, as such divergence is the primary cause of biased data loss (75). Following the selection of a guide sequence, the reads were mapped and aligned using the mapping algorithm described previously (76). During mapping, site indexes in relation to HXB2 were also maintained. Next, to reduce the diversity between the reads and reference sequence, a consensus was generated across each site of the guide sequence, and the reads were remapped to this final consensus template. Reads spanning the 3′ end of Gag, PR, RT, and INT were then translated and assembled for genotyping.
(ii) Variant Caller.
Variant calling, which is the identification and calculation of the frequency of each amino acid present in each genomic position, was calculated using a table generated by the Viral Read Mapper as input, which includes the nucleotide frequencies at each position relative to the reference sequence and numbering relative to the HIV-1B-HXB2 reference strain. The coverage, indel, codon, and residue frequencies at each position were also listed. Variant Caller summarized the results in a graphical interface, with a particular focus on the sites of known drug resistance based on the latest edition of the International AIDS Society-USA (IAS-USA) HIV drug resistance mutations list (77). A list of the amino acids at these positions and their frequencies was exported as a tabulated text file and used with the HIVdb Program Genotypic Resistance Interpretation Algorithm from the Stanford University HIV Drug Resistance Database (see http://hivdb.stanford.edu) to infer the levels of susceptibility to protease, reverse transcriptase, and integrase inhibitors.
In addition, for each data set, reads spanning amino acid positions (i) 50 to 85 in the protease (HXB2 2400 to 2508), (ii) 180 to 215 in the RT (HXB2 3087 to 3195), (iii) 130 to 165 in the integrase (HXB2 4617 to 4725), and (iv) 1 to 35 in the V3 region (HXB2 7110 to 7217) were extracted, truncated, and translated for phylogenetic analysis and HIV-1 coreceptor tropism prediction as described below. Within each data set, only one representative of any identical variant was maintained, but the overall frequency was stored. All variants with a frequency of >1 within the population were aligned using ClustalW (78) and the phylogeny reconstructed using the neighbor-joining statistical method as implemented within MEGA 5.05 (79). In this study, a minority variant was defined as a variation detected at >1% (based on the intrinsic error rate of the system as described below) and <20% of the virus population, corresponding to those mutations that cannot be determined using population sequencing (44,–48).
Genotypic HIV-1 coreceptor tropism determination.
HIV-1 coreceptor tropism was predicted from population and deep sequencing V3 sequences using Geno2Pheno (27). Regarding the population V3 sequences, nucleotide mixtures were considered when the second highest peak in the electropherogram was >25% and the nucleotide mixtures translated into all possible permutations. Geno2Pheno with FPR of 2.5% and 5.75%, based on optimized cutoffs associated with the analysis of clinical data from MOTIVATE (2.5% and 5.75%) (63), or an FPR of 10%, according to the recommendations from the European Consensus Group on the clinical management of HIV-1 tropism testing as described on the Geno2Pheno website (http://coreceptor.bioinf.mpi-inf.mpg.de/index.php), was used for the clinical samples obtained from the Madrid and Seville cohorts, respectively. In the case of deep sequencing V3 sequences, reads spanning amino acid positions 1 to 35 in the V3 region (HXB2 7110 to 7217) were extracted and truncated for HIV-1 coreceptor tropism determination using Geno2Pheno (27) with an FPR of 3.5%, based on optimized cutoffs for determining HIV-1 coreceptor usage, as previously described (36, 63, 80). Deep sequencing V3 sequences usually spanned 105 nucleotides (35 amino acids), with some minor discrepancies associated with natural HIV-1 variation (81, 82), which led to V3 sequences with an open reading frame of 96, 99, 102, 108, or 111 nucleotides, all starting and ending with a cysteine codon, i.e., TG(T/C). V3 reads with stop codons (TGA, TAA, or TAG) and/or where the nucleotide length was not a multiple of 3 (e.g., 101, 103, 104, and 106), mostly associated with naturally or methodology (PCR or sequencing)-induced insertions and/or deletions, were not included in the analysis. Deep sequencing of the V3 region was considered unsuccessful if reads from the majority variants had to be omitted from the analysis. Finally, blood plasma samples were classified as containing non-R5 viruses if ≥2% of the individual sequences, as determined by deep sequencing, were predicted to be non-R5 (36, 63).
Statistical analyses.
The descriptive results are expressed as median values, interquartile ranges, standard deviations, and confidence intervals. Pearson's correlation coefficient was used to determine the strength of the associations between categorical variables. A paired t test was used to compare the number of drug resistance mutations detected by population and deep sequencing in the same sample. All differences with a P value of <0.05 were considered statistically significant. The kappa coefficient, which assesses a chance-adjusted measure of the agreement between any number of categories, was calculated using ComKappa3 version 3.0.1 (83) to quantify the concordance among the different HIV-1 tropism determinations. All statistical analyses were performed using GraphPad Prism version 6.0b (GraphPad Software, La Jolla, CA), unless otherwise specified. The gag-p2/NCp7/p1/p6/pol-PR/RT/IN and/or env-V3 nucleotide sequences obtained in this study by deep sequencing have been submitted to the Los Alamos National Laboratory HIV-DB Next-Generation Sequence Archive.
RESULTS
Characterization of the reverse transcription-PCR amplification step.
As described in Materials and Methods and as shown in Fig. 1, the novel HIV-1 genotyping and coreceptor tropism assay requires reverse transcription-PCR amplification of three amplicons covering the HIV-1 gag-p2/NCp7/p1/p6/pol-PR/RT/IN- and env-C2V3-coding regions. The characterization of the two overlapping PCR products (1,657 nt and 2,002 nt) comprising the subgenomic HIV-1 region spanning the Gag proteins p2, p7, p1, and p6 and the protease, reverse transcriptase, and integrase coding regions was detailed previously as part of a study describing ViralARTS HIV, an HIV-1 phenotyping assay (13). Thus, here, we focused in the characterization of the 480-nt fragment of the HIV-1 env gene, corresponding to the C2V3 region of the surface glycoprotein (gp120). The sensitivity of the reverse transcription-PCR amplification step was tested by analyzing 79 plasma samples obtained from the ACTU (Cleveland, OH). Blood samples from HIV-infected individuals with plasma viral loads ranging from 1,000 to >10,000 copies of viral RNA/ml were used to PCR amplify the C2V3 fragment. Similar to the results observed with the two overlapping fragments spanning the gag-p2/NCp7/p1/p6/pol-PR/RT/IN-coding region (13), reverse transcription-PCR products of the correct size were consistently obtained (92% [73/79]) in plasma samples with ≥1,000 copies/ml of HIV RNA (Table 1).
TABLE 1
Sensitivity of reverse transcription-PCR amplification of the C2V3 region of the HIV-1 env gene (480 nt)
Viral load (copies/ml) | % positive samples by RT-PCR (no. of positive samples/total no. of samples tested)a |
---|---|
1,001–5,000 | 95 (19/20) |
5,001–10,000 | 85 (17/20) |
>10,000 | 95 (37/39) |
Highly reproducible success in reverse transcription-PCR amplification of the specific HIV-1 gag-p2/NCp7/p1/p6/pol-PR/RT/IN- and env-C2V3 products was obtained when testing 15 plasma samples with different viral loads. The details of these tests using two different operators with different lots of critical reagents and over a 7-day period are described in Weber et al. (13) for the gag-p2/NCp7/p1/p6/pol-PR/RT/IN fragments, and in Table S1 in the supplemental material for the env-C2V3 fragment. The specificities of the reverse transcription-PCR primers and reactions for the env-C2V3 fragment were analyzed using nucleic acids from a series of RNA and DNA viruses (i.e., BKV, CMV, HSV-1, HSV-2, VZV, HBV, HCV, and EBV). As expected, no cross-reactivity was observed with any of these viruses, as all reverse transcription-PCRs failed to generate any detectable amplicons (see Fig. S1 in the supplemental material). Similar results were obtained for the gag-p2/NCp7/p1/p6/pol-PR/RT/IN fragments, as previously described (13).
Finally, although most of the HIV-1 genotyping and coreceptor tropism determinations are performed in North America, Europe, and Australia, where subtype HIV-1 strains are predominant (see http://www.who.int/hiv/pub/global_report2010/en/index.html), it was important to test the ability of the assay to work with more worldwide prevalent non-B HIV-1 variants. For that, the env-C2V3 fragment was reverse transcription-PCR amplified from 33 diverse HIV-1 isolates, including five subtype A (HIV-1A-92RW009, HIV-1A-93RW020, HIV-1A-92UG029, HIV-1A-V115, and HIV-1A-V120,), five subtype B (HIV-1B-92BR014, HIV-1B-92TH593, HIV-1B-US714, HIV-1B-92US727, and HIV-1B-92US076), five subtype C (HIV-1C-92BR025, HIV-1C-C18, HIV-1C-C20, HIV-1C-C21, and HIV-1C-C22,), six subtype D (HIV-1D-94UG108, HIV-1D-92UG038, HIV-1D-93UG065, HIV-1D-V89, HIV-1D-V122, and HIV-1D-V126,), six subtype F (HIV-1F-93BR029, HIV-1F-93BR020, HIV-1F-VI820, HIV-1F-V164, HIV-1F-CA16, and HIV-1F-CA20), two subtype G (HIV-1G-RU570 and HIV-1G-RU132,), and four circulating recombinant forms (HIV-1AE-CMU02, HIV-1AE-CMU06, HIV-1AE-92TH021, and HIV-1BF-93BR029). Amplicons of the correct sizes were obtained for the env-C2V3 (see Table S2 in the supplemental material) and gag-p2/NCp7/p1/p6/pol-PR/RT/IN fragments (13) from all HIV-1 group M isolates analyzed, while negative or inconclusive results were obtained with the HIV-2CBL-20 strain (data not shown).
Estimation of the intrinsic error rate of the assay.
Point mutations, insertions, and deletions (indels) can be introduced in the PCR amplification and sequencing steps of any deep-sequencing-based assay (84). Therefore, it was important to calculate the intrinsic (combined) error rate of our novel HIV-1 genotyping and coreceptor tropism assay, since this value might affect the practical limit of detection of the assay. For that, the pNL4-3-hRluc plasmid containing the entire genome of the wild-type HIV-1NL4-3 strain (69) was transformed into Electrocomp TOP10 bacteria (Invitrogen). One bacterial colony was grown overnight in 10 ml of bacterial culture, and the plasmid DNA was purified and transformed again into bacteria. Ten individual and theoretically identical bacterial colonies were used for direct PCR amplification of the gag-p2/NCp7/p1/p6/pol-PR/RT/IN and env-C2V3 fragments (Fig. 2A), and these were sequenced using the same protocol utilized with the clinical samples. The quality of the DNA sequences was analyzed and the reads were filtered in the Ion Torrent server using a Phred quality score of 20 (Q20), which provides a base call accuracy of 99% (i.e., a 1 in 100 probability of an incorrect base call). The average coverages (sequencing depth) per nucleotide position for the 10 clones were 5,750 (range, 681 to 15,614) and 3,797 (range, 942 to 9,981) for the gag-p2/NCp7/p1/p6/pol-PR/RT/IN and env-V3 regions, respectively (Fig. 2B). For each individual NL4-3 clone, the reads were independently mapped to the pNL4-3-hRluc reference sequence (13, 69), and all point mutation and indel information in relation to the reference was analyzed using Segminator II (76).
![An external file that holds a picture, illustration, etc.
Object name is zac0041427490002.jpg An external file that holds a picture, illustration, etc.
Object name is zac0041427490002.jpg](/http/europepmc.org/articles/PMC4023761/bin/zac0041427490002.jpg)
Error rate determination. (A) The pNL4-3-hRluc plasmid containing the entire genome of the wild-type HIV-1NL4-3 strain (69) was transformed into bacteria and the gag-p2/NCp7/p1/p6/pol-PR/RT/IN and env-C2V3 fragments were PCR amplified and deeply sequenced from 10 individual colonies. The reads from each individual NL4-3 clone were independently mapped to the pNL4-3-hRluc reference sequence (13, 69) using Segminator II (76). (B) Coverage, i.e., number of reads per nucleotide position, for the 10 NL4-3 clones. (C) Overall (point mutation, insertions, and deletions) error rate per nucleotide position calculated using a Phred quality score of 20. The mean error rates for the gag-p2/NCp7/p1/p6/pol-PR/RT/IN and env-V3 regions are indicated in green, while the minimum thresholds to detect mutations in the minority HIV-1 variants with this novel assay (1%) are indicated by dashed red lines. (D) Overall error rate in positions associated with drug resistance. Only codon changes with error rates >1% are indicated, i.e., L10, M46, and F53 in the protease region, F77, K101, V179, and G190 in the RT region, G193 in the integrase region, and amino acid 11 in the V3 region. Homopolymeric regions, defined as four or more identical consecutive nucleotides, are indicated as vertical gray bars.
Although all 10 NL4-3 clones were expected to have no mutations (i.e., point mutation and/or indels) relative to the pNL4-3-hRluc reference sequence, a number of errors were observed throughout the p2/NCp7/p1/p6/pol-PR/RT/IN and V3 regions, with error rates ranging from 0% to 29% (mean, 0.39%) and 0% to 9.5% (mean, 0.37%), respectively (Fig. 2C). The average error frequencies due to point mutations were 0.17% (range, 0% to 2.5%) and 0.12% (0% to 0.3%) for the p2/NCp7/p1/p6/pol-PR/RT/IN and V3 regions, respectively, whereas the average error rates associated with indels were 0.22% (range, 0% to 28%) and 0.25% (0% to 9.2%), respectively. Most of the positions with a total (point mutation plus indels) error rate of >1% corresponded to the last nucleotide of a homopolymeric region, defined as four or more identical consecutive nucleotides (data not shown). Some of these nucleotide positions corresponded to codons that have been associated with resistance to antiretroviral drugs, e.g., L10 in the protease (3.5%), K101 in the RT (3%), and G193 in the integrase (10.5%), or with coreceptor tropism, e.g., position 11 in the V3 region (2.6%) (Fig. 2D, Table 2). Interestingly, most of the errors in these (homopolymeric) positions corresponded to indels, with a limited number of point mutation errors, e.g., L10 (3.3% versus 0.22%), K101 (2.7% versus 0.24%), G193 (9.6% versus 0.89%), and position 11 in the V3 (2.3% versus 0.25%), respectively (Table 2). Therefore, considering that (i) the overall error rates for the gag-p2/NCp7/p1/p6/pol-PR/RT/IN and env-V3 regions were 0.39% and 0.37%, respectively, (ii) the point mutation error rates were <1% for all the codons associated with drug resistance, and (iii) the Variant Caller in the DeepGen Software Tool Suite identifies and filters out the indels, it was reasonable to define a frequency of 1% as the minimum threshold to detect mutations in minority HIV-1 variants with this novel assay.
TABLE 2
Error rate distribution
Genomic regiona | Codonb | Sequencec | Error rate (mean ± SD)d | ||
---|---|---|---|---|---|
Total | Point mutations | Indels | |||
gag-pol | NAe | NA | 0.0039 ± 0.0001 | 0.0017 ± 0.0001 | 0.0022 ± 0.0001 |
env-V3 | NA | NA | 0.0037 ± 0.0002 | 0.0012 ± 0.0001 | 0.0025 ± 0.0002 |
Protease | L10 | cgacccCTCgtc | 0.0353 ± 0.0086 | 0.0022 ± 0.0002 | 0.0331 ± 0.0051 |
M46 | accaaaaATGata | 0.0920 ± 0.0056 | 0.0013 ± 0.0001 | 0.0907 ± 0.0043 | |
F53 | aggtTTTatc | 0.0436 ± 0.0034 | 0.0021 ± 0.0002 | 0.0415 ± 0.0022 | |
RT | F77 | agatTTCag | 0.0205 ± 0.0028 | 0.0018 ± 0.0001 | 0.0187 ± 0.0011 |
K101 | gttaAAAcag | 0.0295 ± 0.0044 | 0.0024 ± 0.0002 | 0.0271 ± 0.0032 | |
V179 | ataGTCatc | 0.0756 ± 0.0122 | 0.0036 ± 0.0002 | 0.0720 ± 0.0091 | |
G190 | gtaGGAtct | 0.0221 ± 0.0036 | 0.0012 ± 0.0001 | 0.0209 ± 0.0021 | |
Integrase | G193 | attgggGGGtac | 0.1048 ± 0.0277 | 0.0089 ± 0.0003 | 0.0959 ± 0.0145 |
V3 | R11 | aagaaaaAGTatc | 0.0256 ± 0.0037 | 0.0025 ± 0.0002 | 0.0231 ± 0.0020 |
Performance of the novel deep-sequencing-based HIV-1 genotypic and coreceptor tropism assay.
The measure of success for any deep-sequencing-based assay depends on its ability to generate the maximum number of reads per sequencing run (individual sequences), which then allows for the detection of minority variants within the HIV-1 population. This inherent quality is the sum of a series of metrics, including, but not limited to (i) the number of samples multiplexed and sequenced per run, (ii) chip-loading efficiency, (iii) the total number of quality reads, (iv) mean read length, and (v) the sequencing coverage at each nucleotide position. Most of the deep sequencing runs described in this study involved multiplexing up to 96 individual samples per sequencing reaction, a number that ensured a minimum coverage of 1,000 per nucleotide position sequenced that was required to secure the detection of a minor variant present in at least 1% of the population (85). Efficient loading of Ion Sphere particles into the Ion 318 Chip proved to be user dependent (mean, 72%; range, 60% to 84%). The total number of quality reads was proportional to the chip-loading efficiency (empty wells), with other parameters, such as enrichment (no template), polyclonality (ISPs with excess DNA library), test fragments, and primer dimers, potentially affecting the final number of total reads in this study (mean, 3,827,323 reads; range, 3,051,463 to 4,936,375 reads). We used the Ion PGM sequencing 200 kit version 2 in all sequencing runs, generating an average read length of 147 bp (range, 119 bp to 178 bp). As expected, the average coverage varied with each sequencing run, correlating mostly with the number of multiplexed samples per sequencing reaction, e.g., 20 samples (mean, 9,008; range, 3,776 to 15,458 in the gag-p2/NCp7/p1/p6/pol-PR/RT/IN region and mean, 6,494; range, 2,322 to 8,599 in the env-V3 region) or 96 samples (mean, 4,485; range, 1,612 to 7,274 in the gag-p2/NCp7/p1/p6/pol-PR/RT/IN region and mean, 1,017; range, 966 to 1,070 in the env-V3 region).
As described above, we calculated the error rate of the HIV-1 genotyping and coreceptor tropism assay to be <1%, and we incorporated this error-defined cutoff into the evaluation of the analytical sensitivity and the limit of detection for minority HIV-1 variants. For that, we evaluated extensively the analytical sensitivity of the test to detect and quantify drug resistance mutations (gag-p2/NCp7/p1/p6/pol-PR/RT/IN) and non-R5 variants (env-V3) within mixtures of viral populations. First, we sequenced five plasmid mixtures that contained one or two drug resistance mutations in the RT- and/or IN-coding regions at a frequency of 5% or 10% (i.e., mixture of plasmids containing the respective mutations with a plasmid comprising the wild-type HIV-1HXB2 sequence). As shown in Fig. S2 in the supplemental material, all mutations were detected and quantified at the expected proportions, including those at a frequency of 5% of the total population. Next, in order to quantify more accurately the analytical sensitivity of the assay, we mixed DNA from a plasmid containing a patient-derived multidrug-resistant gag-p2/NCp7/p1/p6/pol-PR/RT/IN fragment in the X4 HIV-1NL4-3 backbone (08-180) (68) with DNA from a plasmid containing the genome of the wild-type HIV-1NL4-3 virus carrying the env gene from the R5 HIV-1YU2 virus (69). Plasmid DNA was quantified and dilutions were used to prepare eight mixtures containing the X4 multidrug-resistant 08-180 plasmid at 0%, 0.1%, 1%, 2%, 3%, 5%, 10%, and 100% concentrations at a final concentration of 0.1 ng/ml. This total plasmid concentration (0.1 ng/ml or 100,000 fg/ml) theoretically allowed the detection of 100 fg of the plasmid when diluted to 0.1% of the population using nested PCR (86). Plasmid 08-180pol/NL43-(X4)env was generated by the yeast cloning method, which allows for a better representation of the in vivo HIV-1 quasispecies (13). It contained numerous drug resistance mutations in the protease, RT, and integrase, most of them as majority members of the quasispecies (>99% of the population); however, two amino acid substitutions in the protease were present as minority variants, i.e., L33F at 21.9% and F53Y at 1.7%. Interestingly, most drug resistance mutations were detected in the plasmid mixtures containing approximately 1% of the 08-180pol/NL43-(X4) plasmid, with the exception of substitution T215Y in the RT, which was identified at 0.95% (Fig. 3A; see also Fig. S3 in the supplemental material). As expected, the detection of minority mutations leading to two amino acid substitutions (L33F and F53Y) faded quickly and in proportion to their frequencies in the original population. Similar results were observed during the detection of X4 (NL4-3) V3 sequences, which were detected in a plasmid mixture containing approximately 1% of the 08-180pol/NL43-(X4) env plasmid (Fig. 3A; see also Fig. S3). Unfortunately, and most likely due to the need to remove V3 sequences with odd open reading frames (as described in Materials and Methods), the quantification of X4 sequences in the mixtures containing the 08-180pol/NL43-(X4) env plasmid at 2% and 3% of the population failed or was not accurate, respectively (Fig. 3A; see also Fig. S3).
![An external file that holds a picture, illustration, etc.
Object name is zac0041427490003.jpg An external file that holds a picture, illustration, etc.
Object name is zac0041427490003.jpg](/http/europepmc.org/articles/PMC4023761/bin/zac0041427490003.jpg)
The ability of the novel HIV-1 genotyping and coreceptor tropism assay to detect drug resistance mutations (in the gag-p2/NCp7/p1/p6/pol-PR/RT/IN fragment) and non-R5 variants (in the env-V3 region) within mixtures of viral populations was used to calculate the analytical sensitivity of the test. (A) A gag-p2/NCp7/p1/p6/pol-PR/RT/IN PCR product was obtained from an antiretroviral-experienced patient (08-180) and used to construct p2-INT recombinant viruses based on the yeast cloning method to maintain the HIV-1 quasispecies (13, 68). This plasmid preparation contained the pol gene from the patient and the env gene from the CXCR4-tropic HIV-1NL4-3 strain, i.e., 08-180 pol/NL43-(X4)env. Plasmid NL4-3 pol/YU2(R5)env contains the genome of the wild-type HIV-1NL4-3 virus carrying the env gene from the R5 HIV-1YU2 virus (69). A series of plasmid mixtures were created by mixing 0.1%, 1%, 2%, 3%, 5%, and 10% of the 08-180pol/NL43-(X4)env plasmid with the corresponding amount of the NL4-3pol/YU2(R5)env plasmid at a final concentration of 0.1 ng/ml. DNA from the entire plasmid mixtures, together with the two individual plasmids as controls (100%), was purified and deeply sequenced as described in Materials and Methods. (B) Plasma containing a patient-derived multidrug-resistant gag-p2/NCp7/p1/p6/pol-PR/RT/IN recombinant virus constructed using the X4 HIV-1NL4-3 backbone (08-194) and a wild-type HIV-1NL4-3 virus carrying the env gene from the R5 HIV-1YU2 virus was mixed at 0%, 1%, 5%, and 100% of the 08-194 virus at a final concentration of 100,000 copies/ml. The frequency of each mutation detected in the original population at ≥1% of the population, the threshold calculated based on the intrinsic error rate of the assay, is indicated (mean ± standard deviation from quadruplicate experiments in the case of the experiment using mixtures of HIV-infected plasma). Amino acid substitutions detected at a frequency of <90% of the population are indicated in red.
Finally, in order to mimic the first steps of the assay (i.e., RNA purification and reverse transcription-PCR) under controlled conditions, HIV-1-seronegative plasma samples were spiked with two viruses, the first being a patient-derived multidrug-resistant gag-p2/NCp7/p1/p6/pol-PR/RT/IN recombinant virus constructed using the X4 HIV-1NL4-3 backbone (08-194) (68) and the second being a wild-type HIV-1NL4-3 virus carrying the env gene from the R5 HIV-1YU2 virus (69). The plasma HIV-1 RNA (viral) load was determined (Cobas AmpliPrep/Cobas TaqMan HIV-1 test version 2.0; Roche) and dilutions were used to prepare four mixtures containing the X4 multidrug-resistant 08-194 virus at 0%, 1%, 5%, and 100% concentrations in a final viral load of 100,000 copies/ml. This was the average viral load in plasma samples obtained from highly antiretroviral-experienced patients, usually carrying multidrug-resistant viruses, from recent previous studies in our laboratory (13, 68). The viral RNA was purified, reverse transcription-PCR amplified, bar coded in quadruplicate, and deeply sequenced as described in Materials and Methods. As expected, all drug resistance mutations from 08-194 were detected when the mixture contained 100% of this virus (Fig. 3B). Interestingly, a few amino acid substitutions were identified by deep sequencing that were not detected in the original study using Sanger sequencing (68), e.g., F53L (2.2%), V77I (17.4%), and I93V (19.7%) in the protease and L100I (3.5%) in the RT coding regions (Fig. 3B). All mutations present at a frequency of >50% in the original 08-194 virus were detected in the 5%:95% (08-194-to-wild type) mixture; however, none of these mutations were identified when the mixture included 1% 08-194 virus (Fig. 3B). Similar results were observed in the env gene, i.e., X4 sequences corresponding to the V3 region of the HIV-1NL4-3 were detected when present at a frequency of 100% and 5% in the viral mixture (Fig. 3B).
The reproducibility of the HIV-1 genotyping and coreceptor tropism assay was evaluated by testing samples from the wild-type HIV-1 control strain (NL4-3), an antiretroviral-naive individual (12-596), and two antiretroviral-experienced (08-198 and 12-069) individuals. The four samples were reverse transcription-PCR amplified in triplicate (3×), each amplicon was bar coded in quadruplicate (4×), and DNA libraries were prepared in duplicate (2×) and then sequenced twice (2×), for a total of 48 sequences per virus (Fig. 4A). First, reads with a frequency of >1 corresponding to 105-bp fragments from the protease, RT, integrase, and V3 regions were used to construct neighbor-joining phylogenetic trees to quantify intra- and interpatient genetic distances and rule out any potential cross-contamination. Figure 4B shows a clear virus-dependent clustering of sequences in all four HIV-1 regions. As expected, the interpatient genetic distances were larger than the range of intrapatient genetic diversity in the four HIV-1 regions, i.e., 0.0495 (range, 0.0023 to 0.0113), 0.0698 (range, 0.0019 to 0.0150), 0.0554 (range, 0.0001 to 0.0145), and 0.3671 (range, 0.0046 to 0.1067) substitutions per site in the protease, RT, integrase, and V3 regions, respectively. Next, the frequency of each nucleotide at each position was compared among the 16 sequences obtained for each one of the triplicate amplicons (n = 48) for all four viruses in the gag-p2/NCp7/p1/p6/pol-PR/RT/IN and env-V3 regions. Statistically significant correlations were observed when the three sets of 16 sequences were compared for each virus, with r values ranging from 0.9857 to 0.9996 (P < 0.0001, Pearson's coefficient correlation) (Fig. 4C). More important, all 48 sequences detected the same amino acids with similar frequencies in each position in all four viruses. This was evident when only positions associated with drug resistance in the protease, RT, and integrase regions were evaluated (Fig. 4D). Wild-type amino acids were basically the only ones identified in these positions (range, 99.6% to 100%) in the NL4-3 reference virus, while various mutations and different frequencies were detected in the patient-derived samples (Fig. 4D and andE).E). Finally, a series of minor amino acid substitutions were identified repeatedly in the patient-derived samples (08-198, 12-069, and 12-596) at frequencies below the limit of detection of Sanger sequencing, e.g., A98G (mean ± standard deviation, 3.3% ± 0.8%) and K103N (8.5% ± 2.6%) in virus 12-069 (Fig. 4E).
![An external file that holds a picture, illustration, etc.
Object name is zac0041427490004.jpg An external file that holds a picture, illustration, etc.
Object name is zac0041427490004.jpg](/http/europepmc.org/articles/PMC4023761/bin/zac0041427490004.jpg)
Assay reproducibility. (A) Samples from two antiretroviral-naive (NL4-3 and 12-596) and two antiretroviral-experienced (08-198 and 12-069) individuals were reverse transcription-PCR amplified in triplicate, each amplicon was bar coded four times, and two DNA libraries were prepared and sequenced in duplicate for a total of 48 sequences per sample. (B) Neighbor-joining phylogenetic trees constructed using reads with a frequency of >1 corresponding to 105-bp fragments from the protease, RT, integrase, and V3 regions. Each color-coded dot represents a unique variant; frequency is not depicted. Bootstrap resampling (1,000 data sets) of the multiple alignments tested the statistical robustness of the trees, with percentage values of >75% indicated by an asterisk. s/nt, substitutions per nucleotide. (C) Pearson's correlation coefficient was used to determine the strength of association between the frequency of each nucleotide at each position among the 16 sequences obtained for each one of triplicate amplicons (n = 48) for all four viruses in the gag-p2/NCp7/p1/p6/pol-PR/RT/IN and env-V3 regions. Over 135,000 and 5,000 points are included in each one of the gag-p2/NCp7/p1/p6/pol-PR/RT/IN and env-V3 plots, respectively. r, correlation coefficient; p, two-tailed P value. (D) Amino acids detected in codons associated with drug resistance in the protease, RT, and integrase regions according to the IAS-USA (77). Drug resistance mutations with a frequency of ≥20% (blue, green, or purple), <20% (yellow), or any other amino acid changes (gray) are indicated. Only amino acid substitutions with a frequency >1% are depicted. (E) Frequency of amino acids in positions associated with drug resistance (gag-p2/NCp7/p1/p6/pol-PR/RT/IN) or X4 variants (env-V3) found in any of the four samples. Each dot represents the mean and 95% confidence interval, with the exception of the inset (sample 12-069), in which each dot indicates the frequency of amino acids detected in each of the 48 replicates, including their mean ± standard deviation.
Comparison of drug susceptibility determination using deep sequencing (DeepGen HIV) with the current standard HIV-1 genotypic assays based on population sequencing.
As described above, a multitude of HIV-1 drug resistance methods have been developed, but only a few have been deployed in the clinical setting, including several genotypic tests based on population sequencing (10, 43, 44). Here, blood plasma samples from 166 treatment-experienced HIV-infected individuals from two cohorts of patients (Seville and Madrid) were analyzed using standard population-based HIV-1 genotyping and the novel deep-sequencing assay. The mean CD4+ T-cell count in these patients was 353 cells/μl (interquartile range [IQR], 190 to 488) and their mean plasma viral load was 69,459 copies/ml (IQR, 5,218 to 87,000). The length of the antiretroviral treatment varied among individuals, averaging 8.2 years (between years 1989 and 2012), and included a diversity of treatment regimens using a multitude of PI, NRTI, NNRTI, raltegravir, and/or maraviroc. Altogether, a total of 1,701 mutations (379 and 1,322 in the Seville and Madrid cohorts, respectively) in positions associated with drug resistance were detected by both methodologies (i.e., 954 in the protease, 613 in the RT, and 134 in the integrase region) (Fig. 5A). As expected, all drug resistance mutations identified by population sequencing were also detected by deep sequencing, while an additional 1,073 drug resistance mutations (337 and 736 in the Seville and Madrid cohorts, respectively) were detected by deep sequencing only (i.e., 511 in the protease, 1,015 in the RT, and 97 in the integrase region) (Fig. 5A). Overall, the difference in the numbers of drug resistance mutations detected by the two methods was significant, even when the mutations were quantified by drug class; i.e., an average of 3.1, 2.8, and 0.6 additional mutations associated with PI, RTI, and INI, respectively, were detected by deep sequencing compared to population sequencing (paired t test, P < 0.0001) (Fig. 5A). Interestingly, additional PI (mean, 5.3 versus 3.1) and RTI (3.1 versus 1.1) but not INI (0.9 versus 0.8) resistance mutations were identified by deep sequencing in patients from the Seville cohort (Fig. 5A). Unlike some of the HIV-infected individuals from the Madrid cohort, these patients from Seville were not treated with raltegravir. The slight difference in the number of INI mutations detected by deep sequencing in the Seville patients corresponded to the identification of the L101I polymorphism, which has been associated with decreased susceptibility to the INI dolutegravir (87). Figure S4 in the supplemental material shows a comparison of the frequency of amino acids detected by population and deep sequencing in a virus carrying multiple mutations associated with resistance to PI, RTI, and INI.
![An external file that holds a picture, illustration, etc.
Object name is zac0041427490005.jpg An external file that holds a picture, illustration, etc.
Object name is zac0041427490005.jpg](/http/europepmc.org/articles/PMC4023761/bin/zac0041427490005.jpg)
Comparison of the novel HIV-1 genotyping and coreceptor tropism assay with other HIV-1 genotypic phenotypic tests. Plasma samples from 166 treatment-experienced HIV-infected individuals from two cohorts of patients (Seville and Madrid) were analyzed as described in Materials and Methods. (A) Top two plots compare the number of drug resistance mutations detected by standard population (Sanger) and deep sequencing in each patient. The total numbers of drug resistance mutations identified by each sequencing method are indicated. The difference in the numbers of drug resistance mutations (mean ± standard deviation) detected by population and deep sequencing in the protease (PR), reverse transcriptase (RT), and integrase (INT) regions is indicated in the bottom two plots. The mean values are indicated. Statistically significant differences are marked with an asterisk (paired t test, P < 0.0001). Deep, deep sequencing; Pop, population sequencing. The dashed gray line reflects the exact number of mutations detected by deep and population sequencing. (B) Hierarchical clustering analysis was used to group the different HIV-1 coreceptor tropism determinations by similarity. Dendrograms (bootstrap values are indicated) were calculated using the Euclidean distance and complete cluster methods with 100 bootstrap iterations (as described at http://www.hiv.lanl.gov/content/sequence/HEATMAP/heatmap.html). Bootstrap values of >60% are indicated with an asterisk. Green and red blocks indicate the absence or presence of non-R5 (X4) viruses, respectively, as determined by each assay. Concordance between DeepGen HIV and the other HIV-1 coreceptor tropism assays is indicated in bubbles below the graphs. G2Pclin, Geno2Pheno with an FPR of 10%; MCT, 8-day maraviroc monotherapy clinical test (71, 72); ESTA, enhanced sensitivity Trofile assay (42); Trofile, the original version of the Trofile assay (19); VeriTrop, phenotypic HIV-1 tropism assay (23); G2Pmot, Geno2Pheno with a FPR of 2.5% and 5.75% based on optimized cutoffs associated with the analysis of clinical data from MOTIVATE (63).
Comparison of HIV-1 coreceptor tropism data obtained with DeepGen HIV to other phenotypic or genotypic HIV-1 coreceptor tropism assays.
Plasma samples from 114 HIV-infected individuals screened to be treated with a maraviroc-containing regimen, a subset of samples from the same Seville and Madrid cohorts of patients, were analyzed using the novel deep-sequencing-based HIV-1 coreceptor tropism assay. These results were compared with those from a series of genotypic (population sequencing) and phenotypic (ESTA [42], TROCAI [73], and VeriTrop [23]) HIV-1 tropism tests. Only samples with results from all the different tests (i.e., 38 and 76 from Seville and Madrid, respectively) were included. Hierarchical clustering analysis grouped the different HIV-1 coreceptor tropism determinations based on their ability to detect R5 and non-R5 (X4, dual tropic, and/or dual mixed [D/M]) sequences (Fig. 5B). The plasma samples from the Seville cohort were from patients participating in a study evaluating the use of an 8-day maraviroc monotherapy clinical test (MCT), the rationale being that the plasma viral load in patients carrying non-R5 viruses will fail to decrease at least one log10 or will be undetectable in subjects with <1,000 copies/ml after receiving maraviroc in this period of time (71, 72). Overall, in this cohort of patients, the concordance and agreement were high among the different HIV-1 tropism methods, with DeepGen HIV (Geno2Pheno with an FPR of 3.5%) showing good agreement with population sequencing analyzed using Geno2Pheno with an FPR of 10% (84.4%, kappa = 0.37), MCT (82.9%, kappa = 0.44), and ESTA (80%, kappa = 0.47) (Fig. 5B). Similar concordance was observed between ESTA and MCT (85%, kappa = 0.64) and between population sequencing/Geno2Pheno with an FPR of 10% and MCT (82.9%, kappa = 0.51). Interestingly, a perfect agreement (100%, kappa = 1) was observed between MCT and TROCAI, the phenotypic HIV-1 tropism assay performed in Seville (data not shown). Using slightly older samples from the Madrid cohort, DeepGen HIV showed excellent agreement with the original Trofile assay (91.7%, kappa = 0.79) and good concordance with VeriTrop (79.8%, kappa = 0.58) or population sequencing/Geno2Pheno with an FPR of 2.5% to 5.75% (74.4%, kappa = 0.37) (Fig. 5B). Concordance between the original Trofile assay and population sequencing/Geno2Pheno with an FPR of 2.5% to 5.75% was comparable (80.3%, kappa = 0.54), while a 73.7% (kappa = 0.5) agreement was observed between the two phenotypic tests (VeriTrop and Trofile) in baseline samples from these patients (Fig. 5B).
DISCUSSION
Almost 25 years since drug-resistant HIV-1 strains were first reported following initial treatments with zidovudine (AZT) (88), and 28 additional antiretroviral drugs later, the detection and characterization of drug-resistant variants have become essential in the clinical care of HIV-infected patients (3, 6,–8). A multitude of genotypic (sequencing-based) and phenotypic (cell-based) assays have been developed to identify and quantify HIV-1 drug resistance, and in the case of CCR5 receptor antagonists, to determine HIV-1 coreceptor tropism (8, 34, 35, 89). Most of these HIV-1 drug resistance or tropism tests correspond to methods developed and used in research settings; however, a few selected assays are commercially available. All current commercial genotypic HIV-1 drug resistance assays are based on population (Sanger) sequencing (10, 43, 44), limiting the detection of minority variants to those in >20% of the viral population (44,–48). In the case of HIV-1 tropism determination, only a couple of ultrasensitive deep-sequencing-based genotypic tests are currently used in the clinical setting (63, 66). Here, we have developed, characterized, and validated the first all-inclusive HIV-1 genotyping and tropism test (DeepGen HIV) based on deep sequencing, which is able to detect minority drug-resistant and non-R5 HIV-1 variants in a single assay.
Deep sequencing is here to stay. The efficiency, sensitivity, and cost-effectiveness of deep-sequencing-based methodologies have allowed their use in a multitude of biological fields, including the sequencing of whole genomes of animals and plants (90, 91), targeted studies on polymorphisms related to various genetic disorders (92) and cancer (93, 94), and research aimed at decoding the genomes of bacteria (95, 96), fungi (97), and viruses (98, 99). We have also observed an exponential increase during the last 5 years in the number of HIV studies based on deep sequencing. Many of them are related to HIV-1 pathogenesis, transmission, and evolution (98, 100), but the majority are aimed at detecting minority non-R5 HIV-1 variants (24, 26, 36, 63, 65,–67, 75, 101,–104) or low-frequency drug-resistant variants that might lead to treatment failure (51, 60,–62, 84, 105,–108). However, none of the methods used in these studies have been able to simultaneously evaluate resistance to antiretroviral drugs targeting Gag, protease, reverse transcriptase, and integrase and the ability of the virus to use the CCR5 coreceptor to enter the host cell all in a single assay. DeepGen HIV was designed to use the same two overlapping amplicons covering the 3′ end of Gag and the entire pol gene used in the phenotypic HIV-1 assay ViralARTS HIV (13), with the addition of the V3 region of the gp120 in the env gene, which allows the possibility to reflex from the ultrasensitive genotypic test (DeepGen HIV) to the comprehensive cell-based ViralARTS HIV assay using the same PCR products. The 97% amplification success of the two overlapping fragments covering the gag-p2/NCp7/p1/p6/pol-PR/RT/IN region from plasma samples with ≥1,000 copies/ml of HIV RNA (13) was matched here by the successful amplification (92%) of the env-V3 fragment from plasma samples with a similar viral load, which is comparable to the success rates observed in other studies (23, 62, 109). Moreover, we have been able to PCR amplify intermittently these long fragments from plasma samples with viral loads between 50 and 1,000 copies/ml, depending on the quality of the clinical specimen (data not shown). The reverse transcription-PCR step was not only reproducible but ensured successful amplification of the samples from diverse HIV-1 subtypes while avoiding amplification of nonspecific products from endogenous or any of the related viruses tested.
Any methodology based on PCR amplification and deep sequencing endures the same fundamental problem, that is, errors are introduced during the process (67, 84, 85, 110, 111). In fact, a limited number of errors are introduced during the PCR step, with most of the errors produced during deep sequencing, mainly insertions and deletions (84). Thus, it was important here to calculate the intrinsic error rate of the entire system since it undoubtedly affects the limit of detection of the assay. The combined error rates, i.e., those for point mutations, insertions, and deletions, of our new HIV-1 genotyping and coreceptor tropism assay were 0.39% and 0.37% for the gag-p2/NCp7/p1/p6/pol-PR/RT/IN and env-V3 regions, respectively. These values were similar to the average error rates previously reported in the HIV-1 pol or env genes using other deep-sequencing platforms, ranging from 0.3% to 0.98% (62, 66, 67, 84, 85, 110, 112,–115). Approximately a 10-fold-higher combined error rate was observed in nucleotide positions associated with homopolymeric regions, resembling findings from other studies describing similar difficulties when sequencing regions with identical consecutive nucleotides (61, 62, 84, 85, 113). More importantly, some of these homopolymeric regions encompass positions associated with resistance to antiretroviral drugs, which might represent a challenge during the interpretation of minority mutations detected at these positions. Using computational methods, we were able to discern between genuine genetic variation and errors introduced during the sequencing process. For example, codon 193 in the integrase region showed the highest error rate in the entire HIV-1 genomic region analyzed, i.e., 9.6% for indels and 10.5% overall. In HIV-1NL4-3, the 25 nucleotides upstream of this position correspond to a series of identical consecutive nucleotides (Fig. 6), which might explain the elevated error rate in this particular region. However, the Ion Torrent software was able to filter most of the reads with sequencing errors at the expense of reducing approximately 5-fold the coverage around this region (Fig. 2B and and6)6) but still above the 1,000 reads required to guarantee the detection of a minor variant present in ≥1% of the population (85). Moreover, the Variant Caller in the DeepGen Software Tool Suite filtered out all the indels at position 193, and as a consequence, we were able to accurately and repeatedly determine the correct amino acids in all wild-type HIV-1NL4-3 sequences and the different variants from the patient-derived samples (Fig. 6).
![An external file that holds a picture, illustration, etc.
Object name is zac0041427490006.jpg An external file that holds a picture, illustration, etc.
Object name is zac0041427490006.jpg](/http/europepmc.org/articles/PMC4023761/bin/zac0041427490006.jpg)
Coverage around codon 193 of the HIV-1 integrase region. Top two panels summarize the coverage of the 10 HIV-1NL4-3 sequences (top) obtained from the clones used to calculate the intrinsic error rate of the assay (gray lines) and 6 HIV-1NL4-3 sequences obtained from the positive controls in six regular sequencing runs (blue lines). The middle panel includes the coverage of 8 patient-derived HIV-1 sequences (clinical samples labeled 13-xxx). Gray bar indicates the homopolymeric region at and upstream of codon 193, i.e., 25 nucleotides in bold in the nucleotide alignment. Red vertical line depicts codon 193. The bottom panel shows the nucleotide sequences of the 8 patient-derived HIV-1 sequences.
The approach of using a deep-sequencing-based assay to detect minority HIV-1 variants relies on the capacity of the methodology to truthfully detect low-frequency sequences from usually complex virus populations. The limits of detection of these assays are a combination of the intrinsic error rate of the method and the analytical sensitivity of the test. As described above, the ability of our HIV-1 genotyping and coreceptor tropism assay to discern between genetic variation and sequencing error was calculated to be 1%. The analytical sensitivity of the assay, determined using mixtures of plasmid DNA, seems to suggest that the deep-sequencing-based test is able to detect minority variants when present in as low as 1% of the population; however, mutations were detected in the HIV-infected plasma mixtures at a 5% frequency but not at 1%. Five percent is the limit of detection previously calculated when studying HIV using 454 pyrosequencing (62, 66). It is possible that although sequences at a frequency of 1% theoretically can be identified, the bias for the majority members of the population during reverse transcription-PCR amplification, library preparation, and/or sequencing might limit the ability of the assay to efficiently detect minority variants of <5%. Additional studies using a larger number of clinical samples will be needed to verify this limit of detection. Nevertheless, our novel assay was able to reproducibly detect HIV-1 minority variants at frequencies at least 4-fold lower than the 20% reported for the standard Sanger-based sequencing (16, 44,–48).
A multitude of studies have compared the efficacies of phenotypic and genotypic assays to detect and quantify HIV-1 drug resistance (116,–120) or coreceptor tropism (38, 121,–124). While there is usually high concordance between drug resistance methodologies (13, 119, 120), in general, population-based sequencing HIV-1 tropism tests are less sensitive and less specific than phenotypic assays (8, 41). Here, DeepGen HIV identified all the drug resistance mutations detected by the standard population sequencing-based HIV-1 genotyping tests used to monitor 166 individuals in two well-established cohorts of patients. More importantly, >1,000 additional minor drug resistance mutations at frequencies between 1% and 20% were detected using the novel deep-sequencing-based assay. In the cases of detection of minority non-R5 variants and determination of HIV-1 coreceptor tropism, DeepGen HIV showed good concordance with both ESTA and the original Trofile assay (80% and 92%, respectively), similar to the agreement of these phenotypic tests with 454-based sequencing systems used to determine HIV-1 tropism (80% to 87%) (24, 36, 63,–66, 104, 125, 126). Deep-sequencing-based HIV-1 coreceptor tropism assays have shown good concordance with virological responses to CCR5 antagonists, such as maraviroc (36, 63); thus, more sensitive HIV-1 genotyping methods based on deep sequencing may also help detect minority drug-resistant variants earlier, possibly improving the design of antiretroviral treatment schemes.
As described above, DeepGen HIV detected additional (minority) HIV-1 drug resistance mutations and non-R5 variants compared to population sequencing; however, the clinical relevance of these minority members of the viral population is still under debate (33, 51, 55, 57, 127,–129). It is reasonable to assume that under appropriate selection (i.e., drug pressure), minority variants carrying drug resistance mutations will eventually outcompete other members of the viral population, and the earlier these mutations are detected, the sooner the proper strategy can be defined to control the growth of these viruses. However, the threshold for the identification of significant low-level variants may be mutation and/or antiretroviral drug dependent (32); therefore, early detection of drug-resistant HIV-1 minority variants, or non-R5 variants in the case of CCR5 antagonists, may contribute to the prediction of clinical or treatment outcome (32, 36, 57, 63, 127). Additional effort, most likely in the form of longitudinal studies, is needed to better understand the clinical relevance of these minority HIV-1 variants. Here, we have developed, characterized, and validated a novel all-inclusive and ultrasensitive HIV-1 genotyping assay based on deep sequencing to simplify the monitoring of patients infected with HIV-1. This first-in-class methodology accurately provides drug resistance information for all protease, reverse transcriptase, integrase, and maturation inhibitors, as well as HIV-1 coreceptor tropism, in a single, more efficient, rapid, and affordable test. This unique assay may provide the platform not only for predicting clinical outcome in HIV-infected individuals—who otherwise might fail treatment regimens—but in the development and clinical assessment of novel antiretroviral drugs.
ACKNOWLEDGMENTS
We thank Eric J. Arts (HIV-1 isolates), Michael R. Jacobs (BKV, CMV, HSV-1, HSV-2, and VZV), David Canaday (EBV), and Donald Anthony (HBV and HCV) from Case Western Reserve University (Cleveland, OH) for providing access to the respective RNA or DNA viruses. We also thank Benigno Rodriguez and the Case Western Reserve University/University Hospitals Center for AIDS Research (grant no. NIH P30 AI036219) and Vicente Soriano (Hospital Carlos III, Madrid, Spain) for providing clinical samples for the characterization and validation studies. We are grateful to Michael D. Miller and Kirsten L. White (Gilead Sciences, Inc., Foster City, CA) for providing some of the plasmid mixtures used in the analytical sensitivity experiments.
E.R.-M. and M.L. were supported by Redes Telemáticas de Investigación Cooperativa en Salud (RETICS), Red de Investigación en SIDA (no. RIS RD12/0017/0029).
M.E.Q.-M. designed the study, collected and assembled the data, and wrote and drafted the manuscript. R.M.G., A.M.M., and D.W. performed all cell culture, molecular, and sequencing experiments. J.A., F.F., M.E.Q.-M. and D.L.R. designed the DEEPGEN Software Tool Suite, while J.A. and F.F. developed the suite. E.R.-M. and M.L. performed the studies associated with the Seville cohort, including population sequencing, MCT, and TROCAI analyses. R.M.G., C.L.S., and M.E.Q.-M. contributed to the overall analysis of the data. All the authors read and approved the final manuscript.
Footnotes
Published ahead of print 27 January 2014
Supplemental material for this article may be found at http://dx.doi.org/10.1128/AAC.02710-13.
REFERENCES
Articles from Antimicrobial Agents and Chemotherapy are provided here courtesy of American Society for Microbiology (ASM)
Full text links
Read article at publisher's site: https://doi.org/10.1128/aac.02710-13
Read article for free, from open access legal sources, via Unpaywall:
https://aac.asm.org/content/aac/58/4/2167.full.pdf
Citations & impact
Impact metrics
Citations of article over time
Alternative metrics
![Altmetric item for https://www.altmetric.com/details/2256084](https://api.altmetric.com/v1/donut/2256084_64.png)
Discover the attention surrounding your research
https://www.altmetric.com/details/2256084
Smart citations by scite.ai
Explore citation contexts and check if this article has been
supported or disputed.
https://scite.ai/reports/10.1128/aac.02710-13
Article citations
Rate of response to initial antiretroviral therapy according to level of pre-existing HIV-1 drug resistance detected by next-generation sequencing in the strategic timing of antiretroviral treatment (START) study.
HIV Med, 25(2):212-222, 29 Sep 2023
Cited by: 0 articles | PMID: 37775947
Sensitive HIV-1 DNA Pol Next-Generation Sequencing for the Characterisation of Archived Antiretroviral Drug Resistance.
Viruses, 15(9):1811, 25 Aug 2023
Cited by: 0 articles | PMID: 37766218 | PMCID: PMC10536450
Next Generation Sequencing Approaches to Characterize the Respiratory Tract Virome.
Microorganisms, 10(12):2327, 24 Nov 2022
Cited by: 3 articles | PMID: 36557580 | PMCID: PMC9785614
Review Free full text in Europe PMC
Towards Next-Generation Sequencing for HIV-1 Drug Resistance Testing in a Clinical Setting.
Viruses, 14(10):2208, 07 Oct 2022
Cited by: 3 articles | PMID: 36298763 | PMCID: PMC9608942
Protease and gag diversity and drug resistance mutations among treatment-naive Mexican people living with HIV.
BMC Infect Dis, 22(1):447, 10 May 2022
Cited by: 3 articles | PMID: 35538426 | PMCID: PMC9088029
Go to all (49) article citations
Data
Data behind the article
This data has been text mined from the article, or deposited into data resources.
BioStudies: supplemental material and supporting data
Nucleotide Sequences
- (2 citations) ENA - K03455
Similar Articles
To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.
Contribution of human immunodeficiency virus type 1 minority variants to reduced drug susceptibility in patients on an integrase strand transfer inhibitor-based therapy.
PLoS One, 9(8):e104512, 11 Aug 2014
Cited by: 8 articles | PMID: 25110880 | PMCID: PMC4128663
Sensitive cell-based assay for determination of human immunodeficiency virus type 1 coreceptor tropism.
J Clin Microbiol, 51(5):1517-1527, 13 Mar 2013
Cited by: 16 articles | PMID: 23486708 | PMCID: PMC3647936
Novel method for simultaneous quantification of phenotypic resistance to maturation, protease, reverse transcriptase, and integrase HIV inhibitors based on 3'Gag(p2/p7/p1/p6)/PR/RT/INT-recombinant viruses: a useful tool in the multitarget era of antiretroviral therapy.
Antimicrob Agents Chemother, 55(8):3729-3742, 31 May 2011
Cited by: 18 articles | PMID: 21628544 | PMCID: PMC3147656
Emerging drug targets for antiretroviral therapy.
Drugs, 65(13):1747-1766, 01 Jan 2005
Cited by: 66 articles | PMID: 16114975
Review
Phenotypic assays for the determination of coreceptor tropism in HIV-1 infected individuals.
Eur J Med Res, 12(9):463-472, 01 Oct 2007
Cited by: 17 articles | PMID: 17933728
Review
Funding
Funders who supported this work.
Biotechnology and Biological Sciences Research Council (1)
Understanding the evolution and diversity of viral pathogens using next generation sequencing technologies
Professor David Robertson, The University of Manchester
Grant ID: BB/H012419/1
NIAID NIH HHS (1)
Grant ID: P30 AI036219