Next Article in Journal
Characterization of Two Pseudomonas aeruginosa Viruses vB_PaeM_SCUT-S1 and vB_PaeM_SCUT-S2
Previous Article in Journal
Infectious Virions of Bombyx Mori Latent Virus Are Incorporated into Bombyx Mori Nucleopolyhedrovirus Occlusion Bodies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of a Versatile, Near Full Genome Amplification and Sequencing Approach for a Broad Variety of HIV-1 Group M Variants

1
Department of Pathology, New York University School of Medicine, New York, NY 10016, USA
2
Faculty of Medicine and Biomedical Sciences, Department of Biochemistry, University of Yaoundé 1, BP 1364 Yaoundé, Cameroon
3
Faculty of Science, Department of Biochemistry, BP 1364 Yaoundé, Cameroon
4
Center of Research for Emerging and Re-Emerging Diseases (CREMER), Institute of Medical Research and Study of Medicinal Plants, BP 906 Yaoundé, Cameroon
5
Applied Bioinformatics Laboratories (ABL) and Genome Technology Center (GTC), Division of Advanced Research Technologies (DART), New York University Langone Medical Center, New York, NY 10016, USA
6
Medical Diagnostic Center, BP 15810 Yaoundé, Cameroon
7
Chantal Biya International Reference Center for Research on HIV/AIDS Prevention and Management, BP 3077 Messa Yaoundé, Cameroon
8
Manhattan Veterans Affairs New York Harbor Healthcare System, New York, NY 10010, USA
9
Faculty of Medicine and Biomedical Sciences, Department of Microbiology, Parasitology and Infectious Diseases, University of Yaoundé 1, BP 1364 Yaoundé, Cameroon
10
Faculty of Medicine and Biomedical Sciences, Department of Pharmacotoxicology & Pharmacokinetics, University of Yaoundé 1, BP 1364 Yaoundé, Cameroon
*
Author to whom correspondence should be addressed.
Deceased.
Viruses 2019, 11(4), 317; https://doi.org/10.3390/v11040317
Submission received: 5 March 2019 / Revised: 26 March 2019 / Accepted: 27 March 2019 / Published: 1 April 2019
(This article belongs to the Section Animal Viruses)

Abstract

:
Near full genome sequencing (NFGS) of HIV-1 is required to assess the genetic composition of HIV-1 strains comprehensively. Population-wide, it enables a determination of the heterogeneity of HIV-1 and the emergence of novel/recombinant strains, while for each individual it constitutes a diagnostic instrument to assist targeted therapeutic measures against viral components. There is still a lack of robust and adaptable techniques for efficient NFGS from miscellaneous HIV-1 subtypes. Using rational primer design, a broad primer set was developed for the amplification and sequencing of diverse HIV-1 group M variants from plasma. Using pure subtypes as well as diverse, unique recombinant forms (URF), variable amplicon approaches were developed for NFGS comprising all functional genes. Twenty-three different genomes composed of subtypes A (A1), B, F (F2), G, CRF01_AE, CRF02_AG, and CRF22_01A1 were successfully determined. The NFGS approach was robust irrespective of viral loads (≥306 copies/mL) and amplification method. Third-generation sequencing (TGS), single genome amplification (SGA), cloning, and bulk sequencing yielded similar outcomes concerning subtype composition and recombinant breakpoint patterns. The introduction of a simple and versatile near full genome amplification, sequencing, and cloning method enables broad application in phylogenetic studies of diverse HIV-1 subtypes and can contribute to personalized HIV therapy and diagnosis.

1. Introduction

HIV-1 full genome sequencing is a challenging task due to the broad degree of HIV-1 genomic diversity worldwide. Four major groups of HIV-1 exist, groups M, N, O, and P [1,2,3], which are highly divergent and trace back to separate cross-species transmission events from primates to humans [4]. HIV-1 group M has spread globally, causing more than 85% of global HIV infections and can be subdivided into nine subtypes (A–D, F–H, J, and K), at least six sub-subtypes of A and F (A1–A6, F1–F2), currently at least 98 circulating recombinant forms (CRFs) [5], and numerous unique recombinant forms (URFs) [6]. CRFs and URFs are composed of two or more (sub)-subtypes; while CRFs have been identified in at least three epidemiologically unlinked individuals, URFs are defined as recombinants without evidence of onward transmission [5]. The continuous emergence of novel HIV-1 M strains [6,7,8,9] including sequences that remain unclassified or can be attributed to ancient strains [10] poses a significant challenge for HIV-1 molecular surveillance, diagnosis, and therapy [11,12]. Consequently, robust and versatile amplification and sequencing protocols are needed that apply to such strains over the near full genome (NFG).
The vast majority of sequences in the Los Alamos National Library (LANL) HIV sequence database are derived from partial genome sequencing [13]. Partial genome sequencing obscures subtype determinations and limits the identification of epidemiologic signatures and recombination breakpoints, which can be distributed across the entire viral genome [14,15]. To overcome these limitations, PCR-based near full genome sequencing (NFGS) protocols have been introduced, which were mostly applied in cohorts restricted to one or a few subtypes [16,17,18,19,20]. Protocols have been optimized for next-generation sequencing (NGS) [21,22], or the use of oligonucleotide probes combined with real-time PCR [23,24]. The coverage of multiple subtypes has been addressed by universal NFGS approaches [22,25]. Universal approaches rely on primer binding to semi-conserved regions in the HIV genome, and limited coverage and efficiencies can still occur depending on viral load and subtype. The combination of NGS with HIV-specific probes increased detection sensitivity and sequencing depth [26]; however, the detection of emerging and rapidly evolving strains is expected to be challenging. NGS approaches provide large-scale sequencing with in-depth coverage and analyses of minority variants; however, their usage is limited in resource-constrained settings due to the required equipment and expertise, and do not allow downstream processing of the amplicons and cloning. There is still a need for robust and straightforward NFGS protocols suited for a broad diversity of HIV-1 M subtypes and the generation of NFG amplicons adaptable to downstream biological applications.
A flexible NFGS approach has therefore been established by our team using rationally designed primers that bind at highly conserved domains across major HIV-1 group M subtypes. The set of primers enables nested PCRs to span single or multiple functional gene regions. It provides NFG amplification and sequencing through one-amplicon, two-amplicon or multiple-amplicon approaches. Pure HIV-1 M subtypes and highly diverse HIV-1 M URFs from Cameroon, characterized in detail in Banin et al. [27], were used to validate the protocol, which was applied to bulk sequencing, cloning, single-genome amplification (SGA), and third-generation sequencing (TGS).

2. Materials and Methods

2.1. Ethical Clearance

This study was conducted in accordance with the Declaration of Helsinki and the protocol was approved by the Institutional Ethical Review Board of the Cameroon Ministry of Public Health and by the Institutional Review Board at New York University School of Medicine (NYUSoM), New York, USA on April 6, 2018 (protocol: i09-0431). Before sample collection, informed consent was obtained from the study participants, who were all part of a cohort of HIV positive individuals; this cohort is monitored at the Medical Diagnostic Center (MDC) in Yaoundé, Cameroon in collaboration with the New York University School of Medicine (NYUSoM), New York, USA.

2.2. Study Samples

Whole blood samples were collected at the MDC between 2000 and 2015. They were shipped under standard regulatory conditions for biological materials to NYUSoM for processing of plasma and peripheral blood mononuclear cells (PBMCs) using Ficoll gradient centrifugation (Histopaque, Sigma-Aldrich, St. Louis, MO, USA), and stored at −80 °C. Study samples were selected based on previously identified URF sequences through partial genome sequencing in our group and on the availability of plasma volumes ≥500 µL [28,29,30].

2.3. RNA Extraction

Viral RNA was extracted from 500 µL HIV-1 positive plasma using QIAamp Viral RNA Minikit (Qiagen, Valencia, CA, USA). Virions were concentrated by ultracentrifugation at 14,000× g for 2 h at 4 °C. All but 140 µL of the supernatant was removed before addition of lysis buffer and subsequent RNA extraction according to the manufacturer’s instructions. Extracted RNA (30 µL) was either immediately reverse transcribed or stored at −80 °C.

2.4. cDNA Synthesis

Reverse transcription of RNA was performed using SuperScript III (Life Technologies, Carlsbad, CA, USA) with primer OFM19 [31] corresponding to HXB2 position 9632–9604. If amplicons within the 5’ half of the HIV-1 genome could not be obtained using cDNA generated with OFM19 primer, HIV 6352 rev [25] was used as an alternative. For the generation of cDNA, 5.25 µL nuclease-free water (Promega) was mixed with 0.75 µL reverse primer OFM19 or HIV 6352 rev (20µM), 3 µL (10 mM) dNTP mix and 30 µL extracted template RNA. The mix was incubated at 65 °C for 5 min and removed on ice for 1 min. 12 µL 5× first strand buffer, 3 µL (100 mM) DTT solution, 3 µL (40 U/µL) RNaseOUT, and 3 µL (200 U/μL) SuperScript III Reverse Transcriptase were added to the mix to get a total volume of 60 µL. The mix was vortexed and incubated at 50 °C for 1 h followed by an increase to 55 °C for 1 h. SuperScript III Reverse Transcriptase was inactivated at 70 °C for 15 min. Remaining RNA was digested by adding 3 µL RNase H (2 U/µL) and incubation at 37 °C for 20 min.

2.5. Primer Design

The rational design of universal HIV-1 group M amplification and sequencing primers was based on full genome multiple sequence alignments with reference strains downloaded from the LANL database. Our approach is related to PrimerDesign-M, an alignment-based primer design tool from the LANL database, which also served as the first guide to identify appropriate genomic regions for some of our primers. Generally, most primers were designed manually based on sequence alignments and primer criteria as follows. An HIV-1 full genome alignment comprised of 480 reference sequences from all major group M (sub)-subtypes and CRFs [32] was shortened to a representative 40 strain-panel. The latter included four references per eight (sub)-subtypes A1, A2, B, C, D, F1, F2, G, and two CRFs, 01_AE, and 02_AG (Results 3.2), which make up to the most dominant (sub-)subtypes/CRFs prevailing the five inhabited continents of the world.
The procedure was similar for the design of amplification and sequencing primers with the only difference of targeted regions. Amplification primers mostly targeted conserved domains adjacent or at the ends of the genes of interest. Sequencing primers had no spatial constraints; however, appropriate distances (~500 bp) between consecutive primer binding sites should provide sufficient sequence overlap between the multiple complementary sequence fragments (≥100 bp overlap) while aiming for the lowest possible number of primers to cover NFGS. Preferentially, primers were designed to have a length of 18–30 base pairs with 40–60% G/C content and a melting temperature (Tm) of 50–65 °C. Primer pairs should be within 5 °C difference in Tm and exempt of complementary regions to avoid heterodimers. Primers with long stretches of identical base pairs (especially G and C), hairpin structures, or self-complementary regions (homodimers) were avoided, as can be analyzed with a primer design software such as IDT OligoAnalyzer (https://www.idtdna.com/calc/analyzer). If possible, primers should start and end with one or two G/C base pairs. The 3’ end of the primer was most critical to accurately match the homologous regions of the majority of reference sequences. Few ambiguous base pairs were tolerated in the central region of the primer sequence and up to three ambiguity characters were introduced for ≤70% conserved sites per primer. Since accurate primer design and positive prediction do not guarantee efficient, practical implementation, trial and error have to succeed. Amplification primers were screened with at least three study samples (cDNA templates), and if possible, in conjunction with already established pairing primers. Primer selection was based upon efficient PCR amplification yielding strong and clean single bands without smear or evidence of unspecific amplification products in at least one NFG PCR or at least 1/3 of PCR reactions amplifying smaller regions. Out of >70 newly tested amplification primers, 26 primers passed the selection criteria (Results 3.2). Sequencing primers were screened with at least two study samples, i.e., bulk-amplified nested PCR products with single strong bands that enabled positive sequencing reactions using established sequencing primers. Out of >90 newly tested sequencing primers, 64 primers worked sufficiently to determine sequences ≥700 bp in length (quality value ≥16, Macrogen sequencing service, New York, NY, USA) from half of the test samples (Results 3.2) [33].

2.6. PCR Amplification

PCR reactions were performed with PrimeSTAR GXL DNA polymerase, (Clontech, Mountain View, CA, USA) or High-Fidelity Taq polymerase (Life Technologies, Carlsbad, CA, USA). For Prime polymerase, nested PCRs were performed in a total reaction volume of 25 µL per sample containing 15.5 µL nuclease-free water, 5 µL 5X PrimeSTAR GXL Buffer (Mg2+ plus), 2.5 µL dNTP mix (2.5 mM each), 0.5 µL Prime polymerase, 0.5 µL forward and reverse primers, and 1 µL cDNA. For the Taq polymerase, a total reaction volume of 25 µL per sample constituted of 18.85 µL nuclease-free water, 2.5 µL buffer 10X, 1 µL MgSO4, 0.5 µL dNTPs, 0.15 µL Taq polymerase, 0.5 µL forward and reverse primers, and 1 µL cDNA.
The PCR reactions were performed with the following cycling conditions: initial denaturation step of 30 s at 98 °C and 30 cycles of 10 s at 98 °C (denaturation), 15 s at 53 °C (annealing), and elongation at 68 °C for 1 min per 1 kb of the amplification product. After 30 cycles, a final elongation step was performed at 68 °C, and the mix stored at 4 °C. The PCR products were analyzed by gel electrophoresis at 120 V for 30 min on a 0.8% agarose gel. To exclude/minimize unspecific amplification and cross-contamination, samples from different participants were processed separately or sequentially, and nucleotide removal treatments were applied frequently. Multiple sequence alignments (MEGA5.2) and highlighter plots (Los Alamos HIV Database) were done with current and historical PCR products from our lab to exclude cross-contamination.

2.7. Sanger Sequencing and Sequence Editing

Sanger Sequencing of amplified PCR products was done by Macrogen (Macrogen Corp., New York, NY, USA). Assembly and editing of the sequences were done with DNAStar package programs SeqMan Pro and EditSeq (Lasergene, Madison, WI, USA).

2.8. Simplot and Recombinant Drawing Tool

The Simplot software package version 3.2 (Window size 200, step size 20 and 250 bootstrap replicate) was used for breakpoint determination. Schematic illustrations of the various full genomes were done using the recombinant drawing tool provided by the Los Alamos HIV sequence database. A bootstrap support of 70% between a subtype reference strain and a query sequence was used as criteria to assign a subtype to a breakpoint region. If bootstrap values remained below 70% using several pure and CRF reference subtypes, the region was considered unidentified.

2.9. Third-Generation Sequencing

Third-generation sequencing (TGS) was performed using Oxford Nanopore MinIon flow cells, capable of high-throughput sequencing in real-time and generating long DNA reads. TGS was carried out at the NYU sequencing and genomic unit. Unsheared HIV DNA from 12 samples was quantitated with the qubit HS DNA kit (Thermo Fisher Scientific). DNA was normalized to 1 µg in 45 uL of buffer and individually input into the Oxford Nanopore PCR Barcoding kit (EXP-PBC001, R9 Version) according to the protocol outlined by Oxford Nanopore (PBGE_9007_V7_revA_16May2016\ PBGE_9007_v7_revJ_16May2016). Each sample was barcoded with a unique adapter index provided in the kit (using the NEB LongAmp Taq master mix, NEB # M0287S) and run with the following PCR conditions: Initial Denaturation for 3 min at 95 °C, Denaturation for 15 s at 95 °C, Annealing for 15 s at 62 °C, Extension for 12 min at 65 °C, Denaturation to Extension was repeated for an additional 14 cycles to compensate for the large fragment sizes (15 cycles total), and one final extension at the end for 14 min at 65 °C. After PCR, each sample was cleaned with Ampure XP (Beckman Coulter) and barcoded samples were quantitated with the Qubit HS DNA kit. Twelve samples were normalized to 110 ng each and combined into one pool for input into the Nanopore Sequencing kit from Oxford Nanopore (SQK-NSK007). The library prep was followed, as described in the protocol above. The pooled library was eluted for 30 min at 37 °C. The final library was quantitated with the qubit HS DNA kit and added onto a primed SpotON Flow Cell (Oxford Nanopore, #FLO-MIN105). The SpotON Flow cell was run on the MinION Mk 1 (#MIN-MAP002) with the script NC_48Hr_Sequencing_Run_FLO_MIN105. FASTQ files were generated using poretools (version 0.6.0) (PMID: 25143291) and general quality controls were performed. Barcodes were trimmed using cutadapt (version 1.9.1) (http://journal.embnet.org/index.php/embnetjournal/article/view/200) and all the reads were mapped to HIV-I (HXB2) reference genome (K03455.1) for an initial quality assessment using the Burrows-Wheeler Aligner (BWA) (version 0.7.13) mem algorithm [34] and visualized using Integrative Genomics Viewer (IGV) software packages [35]. Phylogenetic analyses were done using MEGA5.2 software [36]. Maximum likelihood phylogenetic trees were generated using 1000 bootstrap replicates [37]. Any HIV unspecific reads or minor HIV variants that clustered with a viral population from an unrelated sample were removed. TGS (sub)populations were averaged to consensus (con) sequences using Consensus Maker (Los Alamos National Library (LANL) Database) (www.hiv.lanl.gov) or SeqMan Pro.

2.10. Single-Gnome Amplification

Single genome amplifications (SGA) are PCRs performed on endpoint-diluted cDNA to preclude amplification artifacts such as template switching that can occur when mixtures of genetically diverse templates are used. Template cDNA dilutions that resulted in no more than 30% positive nested PCR reactions were assumed to contain one amplifiable cDNA template per positive PCR (>80% likelihood) [31]. Specifically, ten nested PCR reactions were performed per cDNA dilution. Multiple cDNA dilutions in the range between 1:2 and 1:1,000,000 were tested for each sample. Dilutions with no more than three positive PCR reactions were considered for SGA analyses and the positive amplicons subjected to direct sequencing. SGAs were performed using High-Fidelity PrimeSTAR GXL DNA polymerase (Clontech, Mountain View, CA, USA) using the same PCR and cycling conditions as described under 2.5.

2.11. Cloning

PCR products were ligated into pcDNA3.1 (blunt-end cloning, PrimeSTAR PCR products) or pCR4 TOPO (3’-A overhang cloning, Taq PCR products) cloning vectors (Life Technologies, Carlsbad, CA, USA). The plasmids were subsequently transformed into One Shot Top10 Chemically Competent E. coli, (Life Technologies, Carlsbad, CA, USA) and cultured overnight using LB-ampicillin agar plates [38].

2.12. Data Storage and Documentation

Near full genome sequences are available from GenBank with accession numbers MK086109-MK086132. Whole sets of TGS sequences are available upon request. Correspondence, data and material requests should be addressed to Ralf Duerr ([email protected]).

3. Results

3.1. Design of an Adaptable NFG Amplification Strategy

The goal was to establish a flexible and straightforward methodology for NFG amplification, cloning, and sequencing applicable to the broad variety of group M subtypes and recombinant forms. The protocol should be suitable for several downstream methods including high-throughput sequencing (NGS or TGS), SGA, cloning, and bulk sequencing, and applicable to modern practice in well-equipped labs as well as to basic procedures in resource-limited countries. The NFGS approaches aimed to cover >8.5 kb including all encoding HIV-1 M genes, i.e., from gag to nef. The primary goal was the reliable and consistent determination of NFGS of highly diverse group M strains irrespective of the number of amplicons needed (Figure 1). The first option was the amplification of HIV-1 NFG in a single PCR using primers targeting the 5’ UTR and 3’ UTR regions. Also, we established a two-amplicon and multiple-amplicon approach. The two-amplicon approach targeted two overlapping “half genome” (HG) sequences. The 5’ HG region (HG1) included gag and pol, while the 3’ HG region (HG2) included env and nef. The accessory genes (vif, vpr, vpu) were included in both half genomes in their overlapping part. The multiple amplicon strategies aimed to cover near full genomes using three to five overlapping amplicons from a selection of eight PCR constructs in total. Each PCR construct encompassed complete genomic units of encoding genes (env also separated into gp120 and gp41).

3.2. Rational Primer Design for NFG Amplification and Sequencing of Diverse Group M Viruses

NFGS approaches on our highly diverse Cameroonian URF samples had only limited success using published NFG primers [16,25] (see below). Therefore, rational primer design was performed using full genome reference sequence alignments with diverse HIV-1 group M strains downloaded from the LANL database (Table 1). Initially, an alignment composed of 480 reference sequences was used belonging to several known pure and CRF clades covering the major branches of each subtype [32].
Subsequently, the alignment was reduced and simplified to a representative 40 strain-panel comprising four reference sequences per eight (sub)-subtypes A1, A2, B, C, D, F1, F2, G, and two CRFs, 01_AE, and 02_AG (Figure 2). The designed amplification primers mostly targeted conserved domains adjacent or at the ends of the HIV-1 M genes of interest under consideration of primer characteristics such as melting temperature, GC content and secondary structure (using IDT web application and/or PrimerDesign-M, LANL database). The newly designed primers were combined with published primers [25,31] that exhibited significant amplification efficiencies in our Cameroonian cohort (Table 1). The newly developed primers targeted regions in the HIV-1 M genome, which were highly conserved across several reference sequences of diverse subtypes, as shown for the NFG and two HG amplicon approaches (Figure 2; Figures S1 and S2). Of significant importance, the terminal nucleotides of the primer binding sites were highly conserved, whereas a few internal positions exhibited moderate variance for some primers/reference strains. The binding sites of three primers for half genome amplifications were fully conserved across the 40 reference strains panel (Figures S1 and S2). Overall, we designed/selected 26 new amplification primers, complemented by ten primers from previous studies [25,31] that proved to be efficient. They were used for the amplification of 15 different constructs including four NFGs, three HGs, and eight smaller regions within the full genome (Figure 1 and Figure 3, Table 1).
Sequencing primers were designed based on the same criteria as described above for amplification primers. They were designed to generate fragments with ≥100 bp overlap when sequences ≥700 bp are obtained (Figure S3). Seventy-one sequencing primers were applied including 43 newly designed forward primers and 21 newly designed reverse primers, complemented by seven reverse primers published by other groups [25,31] (Table 2 and Table 3).
Fourteen sequencing primers were most frequently used based on their broad reactivity with our diverse study samples (Figure S3). These 14 primers provided complete sequence coverage from HxB2 positions 596 to 9604 (~9.1 kb) in NFG one-amplicon approaches. For two-amplicon approaches, eight primers were used to sequence the 1st HG (HxB2 position 776 to 5956; ~5.2 kb) and six primers for the 2nd HG (HxB2 position 5037 to 9533; ~4.6 kb).

3.3. Subtype-Independent Amplification and Sequencing of HIV-1 Near Full Genomes

Initially, we screened published NFGS nested primers for suitability with our cohort samples from Cameroon. Nested PCRs with a primer set previously optimized for NFGS of clade C viruses [16] did not result in any NFG PCR products using three study samples including #6541-6, which worked for most other NFGS approaches tested in this study. However, using a nested primer set recently optimized for multiple HIV-1 subtypes [25], we obtained positive reactions for #6541-6 and two more study samples out of 23 tested specimens (Table 4). Based on these results, we screened the newly designed primers as well as combinations of our and published primers to optimize NFGS across multiple HIV-1 group M subtypes (Table 1).
To establish and optimize the method, we first used four study samples from infections with pure (sub)-subtype viruses (B and F2; #Strain0526 and #LB002-1), a CRF (CRF02_AG; #6506-8), and a URF (#6541-6) (based on preliminary and published partial genome sequencing data of Duerr and Nyambi labs [29,40]) (Figure 4). Subtypes F2 and CRF02_AG were chosen, because of their high prevalence in Cameroon [41]. Two of those constructs could be amplified with a one-amplicon approach, the other two with a two-amplicon approach (Table 4). Subsequently, we analyzed 19 more samples, including another subtype F2 virus (#LB022), another CRF02_AG virus (#LB006-1), and different highly diverse URF samples from Cameroon (Table 4). The 18 URF genomes were composed of subtypes A (A1), F (F2), G, and recombinants CRF01_AE, CRF02_AG, and CRF22_01A1, and revealed non-uniform mosaic NFG sequences [27].
Out of 23 samples, 11 NFGS were obtained with a one-amplicon approach. Another 11 samples needed two amplicons and one sample at least three amplicons for NFGS determination (Table 4). The one-amplicon approaches were based on nested PCRs with primers targeting the 5’ UTR and 3’ UTR ends to obtain approximately 9.1 kb products in a single PCR reaction (Figure 1). The four optimized NFGS primer pairs (NFG1 to NFG 4) obtained yields between 13% and 21.7 % (Table 1; Table 4); the best yield was obtained with the NFG1 approach based on four newly designed primers (21.7% yield). Only one sample (#6541-6) could be amplified using all four NFGS approaches (Table 4). The four one-amplicon approaches complemented each other to an overall 43% yield. In the two-amplicon approach, NFGS were obtained from two overlapping half genomes. Half genome 1 (HG1) consisted of gag, pol and the accessory genes vif, vpr, and vpu (~5.2 kb), whereas the second half genome (HG2) included vif, vpr, vpu, env, and nef (~4.6 kb) (Figure 1 and Table 1). HG1a and b had 45% and 80% yield, respectively, and complemented each other except for one sample, to achieve a 95% yield in combination. Of interest, HG2, which was obtained using a mix of newly designed primers and published primers [25,31] achieved 100% yield and proved to be the most efficient nested primer combination in our set of study samples (Table 4). The amplification of smaller regions for the multiple-amplicon approach included gag (1.5 kb), gagpol (4.5 kb), pol (3.2 kb), vifvpu (1.0 kb), vif/gp120 (3.0 kb), env (3.4 kb), gp120 (1,7 kb), and gp41/nef (1.8 kb) (Figure 1; Table 1 and Table 4) with yields >55%. Since plasma viral load (PVL) has been described as a critical parameter for successful NFG sequencing [16,17,21], we analyzed whether the viral load of a sample predicted the number of amplicons needed for characterizing the NFG. Using the diverse study set of Cameroonian URF viruses, no association was observed between the number of amplicons needed for NFGS and PVL (Figure S4). The lowest PVL in our set of study samples with determined NFGS was 306 cps/mL.

3.4. Comparative NFGS Analysis Using Different Amplification, Cloning, and Sequencing Strategies

To validate our NFGS approach, different amplification, cloning, and sequencing procedures including bulk amplification, high-throughput/deep sequencing, and SGA were applied. The NFGS outcome was compared across platforms using sample NYU6541-6, which gave positive results for most amplification approaches. Bulk sequencing is the traditional method for the genetic determination of viral variants that make up ≥20% of a viral population in a sample, e.g., it is still the standard for genotypic drug resistance testing [41,42]. High-throughput sequencing enables more in-depth analyses and can assess the diversity of viral populations in a sample. As proof of principle, we used the portable MinION TGS technology (Oxford Nanopore Technologies) with the capacity to determine long reads in real-time [43]. Despite limitations in accuracy and depth [44,45], it gave an estimate of whether our diverse amplicon approaches are suited for high-throughput/deep sequencing (Figure 5 and Table 5) [27].
Since amplification artifacts such as recombination events can occur in vitro during bulk amplification, single genome amplification (SGA) [31] was used for validation of our results. Six samples were tested for SGA including NFG and HG constructs. Our results confirmed the results reported by other groups [46,47] that major viral populations detected by deep or bulk sequencing show a good match with SGA sequences with comparable genetic coverage (Figure 5 and Table 5) [27]. Specifically, in vitro recombination rates <1.6% were reported for bulk amplification [46,47]. As a precaution, TGS minority variants with <2% prevalence that were not confirmed by SGA were consequently excluded in our analyses. In contrast to high-throughput sequencing such as NGS or TGS, bulk sequencing and SGA allow the cloning of NFG constructs, which enables further functional characterization in downstream biological assays, such as replicative fitness or phylogenetic drug resistance tests [48,49,50]. An NFGS plasmid generated by Topo cloning of an NFG bulk amplicon was included in the comparative analysis, which revealed comparable characteristics to bulk sequencing as evident by <0.1% genetic distance and similar recombination breakpoint patterns (Figure 5). A comparative Simplot Bootscan analysis was performed to underline the significance of the introduced NFGS approach for accurate subtype/recombinant form determination across platforms. Bulk sequencing, cloning, SGA, and TGS (consensus sequence) yielded similar NFGS results with genetic distances <2% across platforms and overall comparable recombinant breakpoints patterns (Figure 5).

4. Discussion

Here, we report a simple and efficient method for the amplification and sequencing of HIV-1 near full genomes from plasma to cover a broad range of recombinant (URF and CRF) viruses and pure group M (sub)-subtypes. Due to the flexible approach based on one, two, or more amplicons, 100% of the study samples were genotyped over the near full genome. Recent studies performing NFGS out of plasma achieved different efficiencies in the range of 45–100% [16,17,18,20,21,22,25,46], from which most approaches remained below 100% efficiency and/or were restricted to specific subtypes. A one-amplicon approach obtained 67% NFGS efficiency in a subtype C cohort [16], whereas a two-amplicon approach applied to samples of diverse clades yielded 90–92% success rates [25,46]. A study in a US clade B cohort used an alternative two- or three-amplicon strategy to obtain 76% yield for samples with viral loads >10,000 cps/mL. A four-amplicon strategy achieved 93% yield when applied to diverse HIV-1 subtypes and groups [21], however when used in an Indian cohort with clade C and AC recombinant samples, the yield dropped to 45% despite adjustments to an optimized six amplicon strategy [18]. The results of these studies imply and confirm our findings in achieving higher percentages of NFGS when diverse amplicon approaches are complementarily used. In our study, the 43% efficiency of the one-amplicon approach, based on four different NFG protocols, increased to 95% when a two-amplicon strategy based on two different half genome protocols was used, and eventually achieved 100% with the help of multiple amplicons. Two NGS-based methods were recently introduced, which also obtained 100% NFGS rates. They were applied in a clade B cohort using a four amplicon strategy [20] or in a highly diverse study population using mixed four- and six-amplicon NGS libraries [22]. The latter study used an elegant method based on degenerate primers fused to common adapter sequences for efficient library generation, which enabled broad coverage of HIV subtypes and groups; however, the need of NGS equipment and expert knowledge limits its application in resource-constrained countries. Also, the method required considerable optimizations with regards to RNA extraction and cDNA synthesis for suitability with clinical samples.
For NFGS, a viral load threshold of 3,000 cps/mL [16,21] or even 10,000 cps/mL [16,17,21] has been reported. NGS can generate full genome coverage at viral loads >log3.5 cps/mL, while even most sensitive NGS/probes approach gradually decline to partial genome coverage for <log3.5 cps/mL [22,26]. Of interest, we obtained NFGS over a broad range of PVL down to ~300 cps/mL, although it is worth mentioning that in our study, we only had one sample below 3,000 cps/mL PVL included, and plasma was routinely concentrated ~3.5 fold (from 500 µL to 140 µL) before RNA extraction. Other than that, our method needed no additional optimizations. Other groups reported that a substantially higher yield could be obtained using Benzonase pretreatment of plasma before RNA extraction to reduce the background of human RNA and DNA and by increasing the reverse transcription temperature from 42 °C to 47 °C [22]. Doubling the amount of RT enzyme from 200 U to 400 U significantly increased the critical RT step for NFGS in another study [17].
Despite the limitation in sample size (23 subjects), the analysis of a range of pure (sub)-subtypes and recombinant forms including A (A1), B, F (F2), G, CRF01_AE, CRF02_AG, and CRF22_01A1 implies a broad applicability of the presented NFGS approach across diverse (sub)-subtypes/recombinants and geographic regions of the world. Universal usage of our approach is further corroborated by the set of primers, designed to bind several most prevalent pure (sub)-subtypes and recombinant forms on the globe. Nevertheless, future studies are needed to show that the primers and the NFGS strategy are robust on larger sample sets from diverse regions of the world. The introduced procedure of rational primer design coupled with a flexible amplicon strategy allows researchers to improve the protocol further according to their needs.
Our approach provides a universal toolbox that can be used and adapted for different amplification and sequencing applications and is easy to handle for laboratories in resource-limited countries. Besides the flexible methodology, it allows both NFGS as well as focused molecular studies for distinct genomic regions. Since our method does not necessarily depend on high-throughput sequencing (NGS/TGS), which usually determines a one-way road towards genetic read-out, the SGA or bulk amplification can be used for the generation of molecular clones and subsequent downstream applications in biological assays. Notably, the introduced combination of primers from different studies [25,31] including our newly designed primers brought up highly efficient nested PCR primer pairs that may prove useful for purposes beyond NFGS and subtype monitoring. For example, we achieved 100% yield in our HG2 approach, which includes the amplification of the complete env region and accessory genes. It could be a useful alternative for env SGA, pseudovirus generation, screening for neutralizing antibody resistant strains, and genotypic tropism analysis. Another construct encompassing whole pol, gag, and accessory genes (HG1b) was obtained with 80% amplification yield when using four newly designed nested primers. The latter represents a promising approach for genotypic drug resistance testing in research and clinical care covering several essential target sites. Although not tested in the present study, we expect our NFG amplification approaches to performing well with state-of-the-art NGS systems to capture minority variants in quasispecies populations, as we showed recently for the sensitive Illumina MiSeq platform using a related amplification approach restricted to the reverse transcriptase region of pol [41].
Hence, the proposed protocols have broad applicability and can form the basis for downstream applications in phylogenetic and molecular studies, and further project-specific optimizations of primer combinations. While the studied in-depth sequencing approach based on MinION technology needs further optimization regarding accuracy and depth, it indicates future implications in real-time NFGS on-site.

5. Conclusions

In summary, we designed and validated a robust and adaptable protocol for NFGS from plasma, capable of amplifying and sequencing a broad diversity of HIV-1 group M strains including pure (sub)-subtypes and mosaic viral strains (URF and CRF). Despite the limited number of 23 samples, our approach suggests broad applicability based on the semi-conserved primer design and the flexible amplicon approach. Future studies using larger sample numbers from different geographic regions and epidemics will need to confirm its suitability in global HIV surveillance. The introduced NFGS approach may prove valuable to complement existing strategies in structural and functional HIV research including regions with rapid evolution and emergence of novel recombinant viral strains.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/11/4/317/s1, Figure S1: Conserved primer binding sites for HIV-1 half genome 1 amplification. Figure S2: Conserved primer binding sites for HIV-1 half genome 2 amplification. Figure S3: Stepwise sequencing approach. Figure S4: Correlation analysis between plasma viral load and the number of amplicons needed for near full genome sequencing.

Author Contributions

Conceptualization, R.D.; Formal analysis, A.N.B., M.T., M.T., A.K.-J., A.J.N., J.M., D.M., J.N., A.H., P.N.N., C.F. and R.D.; Investigation, A.N.B., M.T., J.S.B., P.Z. and X.W.; Project administration, A.J.N., D.M., A.H., P.N.N. and C.F.; Writing—original draft, A.N.B. and R.D.; Writing—review & editing, A.N.B. and R.D.

Funding

This research was funded by grant number AI083142 from the National Institute of Allergy and Infectious Diseases and grant number TW009604 from the Fogarty International Center, National Institutes of Health. This research was also supported by the Sub-Saharan African Network for TB/HIV Research Excellence (SANTHE), a DELTAS Africa Initiative (grant # DEL-15-006). The DELTAS Africa Initiative is an independent funding scheme of the African Academy of Sciences (AAS)’s Alliance for Accelerating Excellence in Science in Africa (AESA) and supported by the New Partnership for Africa’s Development Planning and Coordinating Agency (NEPAD Agency) with funding from the Welcome Trust (grant # 107752/Z/15/Z) and the UK government.

Acknowledgments

The authors thank the individuals who consented to participate and donated their blood samples for this study. Our appreciation is extended to the Cameroon Ministry of Public Health for supporting this study. Thanks further go to Caroline Kakam, Bladine Asaah, and Flavia Camacho for their assistance in sample collection and manuscript preparation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gurtler, L.G.; Hauser, P.H.; Eberle, J.; von Brunn, A.; Knapp, S.; Zekeng, L.; Tsague, J.M.; Kaptue, L. A new subtype of human immunodeficiency virus type 1 (MVP-5180) from Cameroon. J. Virol. 1994, 68, 1581–1585. [Google Scholar]
  2. Plantier, J.C.; Leoz, M.; Dickerson, J.E.; de Oliveira, F.; Cordonnier, F.; Lemee, V.; Damond, F.; Robertson, D.L.; Simon, F. A new human immunodeficiency virus derived from gorillas. Nat. Med. 2009, 15, 871–872. [Google Scholar] [CrossRef]
  3. Vallari, A.; Holzmayer, V.; Harris, B.; Yamaguchi, J.; Ngansop, C.; Makamche, F.; Mbanya, D.; Kaptue, L.; Ndembi, N.; Gurtler, L.; et al. Confirmation of putative HIV-1 group P in Cameroon. J. Virol. 2011, 85, 1403–1407. [Google Scholar] [CrossRef] [PubMed]
  4. Keele, B.F.; Van Heuverswyn, F.; Li, Y.; Bailes, E.; Takehisa, J.; Santiago, M.L.; Bibollet-Ruche, F.; Chen, Y.; Wain, L.V.; Liegeois, F.; et al. Chimpanzee reservoirs of pandemic and nonpandemic HIV-1. Science 2006, 313, 523–526. [Google Scholar] [CrossRef]
  5. HIV Circulating Recombinant Forms (CRFs). Available online: https://www.hiv.lanl.gov/content/sequence/HIV/CRFs/CRFs.html (accessed on 27 January 2019).
  6. Hemelaar, J.; Elangovan, R.; Yun, J.; Dickson-Tetteh, L.; Fleminger, I.; Kirtley, S.; Williams, B.; Gouws-Williams, E.; Ghys, P.D.; Characterisation, W.-U.N.f.H.I. Global and regional molecular epidemiology of HIV-1, 1990-2015: A systematic review, global survey, and trend analysis. Lancet Infect. Dis. 2018. [Google Scholar] [CrossRef]
  7. Reis, M.; Bello, G.; Guimaraes, M.L.; Stefani, M.M.A. Characterization of HIV-1 CRF90_BF1 and putative novel CRFs_BF1 in Central West, North and Northeast Brazilian regions. PLoS ONE 2017, 12, e0178578. [Google Scholar] [CrossRef]
  8. Villabona Arenas, C.J.; Vidal, N.; Ahuka Mundeke, S.; Muwonga, J.; Serrano, L.; Muyembe, J.J.; Boillot, F.; Delaporte, E.; Peeters, M. Divergent HIV-1 strains (CRF92_C2U and CRF93_cpx) co-circulating in the Democratic Republic of the Congo: Phylogenetic insights on the early evolutionary history of subtype C. Virus Evol. 2017, 3, vex032. [Google Scholar] [CrossRef] [PubMed]
  9. Miao, J.; Ran, J.; Song, Y.; Liu, Y.; Gao, L.; Miao, Z.; Zhang, C.; Feng, Y.; Xia, X. Characterization of a Novel HIV-1 Circulating Recombinant Form, CRF01_AE/B’/C (CRF96_cpx), in Yunnan, China. AIDS Res. Hum. Retroviruses 2017. [Google Scholar] [CrossRef] [PubMed]
  10. Tongo, M.; Dorfman, J.R.; Martin, D.P. High Degree of HIV-1 Group M (HIV-1M) Genetic Diversity within Circulating Recombinant Forms: Insight into the Early Events of HIV-1M Evolution. J. Virol. 2015, 90, 2221–2229. [Google Scholar] [CrossRef]
  11. Hemelaar, J. Implications of HIV diversity for the HIV-1 pandemic. J. Infect. 2013, 66, 391–400. [Google Scholar] [CrossRef] [PubMed]
  12. Zhang, M.; Foley, B.; Schultz, A.K.; Macke, J.P.; Bulla, I.; Stanke, M.; Morgenstern, B.; Korber, B.; Leitner, T. The role of recombination in the emergence of a complex and dynamic HIV epidemic. Retrovirology 2010, 7, 25. [Google Scholar] [CrossRef]
  13. HIV Sequence Database. Available online: https://www.hiv.lanl.gov/content/sequence/HIV/mainpage.html (accessed on 27 January 2019).
  14. Smyth, R.P.; Schlub, T.E.; Grimm, A.J.; Waugh, C.; Ellenberg, P.; Chopra, A.; Mallal, S.; Cromer, D.; Mak, J.; Davenport, M.P. Identifying recombination hot spots in the HIV-1 genome. J. Virol 2014, 88, 2891–2902. [Google Scholar] [CrossRef]
  15. Fan, J.; Negroni, M.; Robertson, D.L. The distribution of HIV-1 recombination breakpoints. Infect. Genet. Evol. 2007, 7, 717–723. [Google Scholar] [CrossRef] [PubMed]
  16. Rousseau, C.M.; Birditt, B.A.; McKay, A.R.; Stoddard, J.N.; Lee, T.C.; McLaughlin, S.; Moore, S.W.; Shindo, N.; Learn, G.H.; Korber, B.T.; et al. Large-scale amplification, cloning and sequencing of near full-length HIV-1 subtype C genomes. J. Virol. Methods 2006, 136, 118–125. [Google Scholar] [CrossRef] [PubMed]
  17. Nadai, Y.; Eyzaguirre, L.M.; Constantine, N.T.; Sill, A.M.; Cleghorn, F.; Blattner, W.A.; Carr, J.K. Protocol for nearly full-length sequencing of HIV-1 RNA from plasma. PLoS ONE 2008, 3, e1420. [Google Scholar] [CrossRef] [PubMed]
  18. Alampalli, S.V.; Thomson, M.M.; Sampathkumar, R.; Sivaraman, K.; U K J, A.J.; Dhar, C.; D Souza, G.; Berry, N.; Vyakarnam, A. Deep sequencing of near full-length HIV-1 genomes from plasma identifies circulating subtype C and infrequent occurrence of AC recombinant form in Southern India. PLoS ONE 2017, 12, e0188603. [Google Scholar] [CrossRef]
  19. Alves, B.M.; Siqueira, J.D.; Garrido, M.M.; Botelho, O.M.; Prellwitz, I.M.; Ribeiro, S.R.; Soares, E.A.; Soares, M.A. Characterization of HIV-1 Near Full-Length Proviral Genome Quasispecies from Patients with Undetectable Viral Load Undergoing First-Line HAART Therapy. Viruses 2017, 9, 392. [Google Scholar] [CrossRef] [PubMed]
  20. Henn, M.R.; Boutwell, C.L.; Charlebois, P.; Lennon, N.J.; Power, K.A.; Macalalad, A.R.; Berlin, A.M.; Malboeuf, C.M.; Ryan, E.M.; Gnerre, S.; et al. Whole genome deep sequencing of HIV-1 reveals the impact of early minor variants upon immune recognition during acute infection. PLoS Pathog. 2012, 8, e1002529. [Google Scholar] [CrossRef] [PubMed]
  21. Gall, A.; Ferns, B.; Morris, C.; Watson, S.; Cotten, M.; Robinson, M.; Berry, N.; Pillay, D.; Kellam, P. Universal amplification, next-generation sequencing, and assembly of HIV-1 genomes. J. Clin. Microbiol. 2012, 50, 3838–3844. [Google Scholar] [CrossRef] [PubMed]
  22. Berg, M.G.; Yamaguchi, J.; Alessandri-Gradt, E.; Tell, R.W.; Plantier, J.C.; Brennan, C.A. A Pan-HIV Strategy for Complete Genome Sequencing. J. Clin. Microbiol. 2016, 54, 868–882. [Google Scholar] [CrossRef]
  23. Heipertz, R.A., Jr.; Ayemoba, O.; Sanders-Buell, E.; Poltavee, K.; Pham, P.; Kijak, G.H.; Lei, E.; Bose, M.; Howell, S.; O’Sullivan, A.M.; et al. Significant contribution of subtype G to HIV-1 genetic complexity in Nigeria identified by a newly developed subtyping assay specific for subtype G and CRF02_AG. Medicine (Baltimore) 2016, 95, e4346. [Google Scholar] [CrossRef]
  24. Kijak, G.H.; Tovanabutra, S.; Sanders-Buell, E.; Watanaveeradej, V.; de Souza, M.S.; Nelson, K.E.; Ketsararat, V.; Gulgolgarn, V.; Wera-arpachai, M.; Sriplienchan, S.; et al. Distinguishing molecular forms of HIV-1 in Asia with a high-throughput, fluorescent genotyping assay, MHAbce v.2. Virology 2007, 358, 178–191. [Google Scholar] [CrossRef]
  25. Grossmann, S.; Nowak, P.; Neogi, U. Subtype-independent near full-length HIV-1 genome sequencing and assembly to be used in large molecular epidemiological studies and clinical management. J. Int. AIDS Soc. 2015, 18, 20035. [Google Scholar] [CrossRef]
  26. Yamaguchi, J.; Olivo, A.; Laeyendecker, O.; Forberg, K.; Ndembi, N.; Mbanya, D.; Kaptue, L.; Quinn, T.C.; Cloherty, G.A.; Rodgers, M.A.; et al. Universal Target Capture of HIV Sequences from NGS Libraries. Front. Microbiol. 2018, 9, 2150. [Google Scholar] [CrossRef] [PubMed]
  27. Banin, A.N.; Tuen, M.; Tongo, M.; Bimela, J.S.; Nanfack, A.; Courtney, C.; Zappile, P.; Heguy, A.; Fokunang, C.; Mbanya, D.; Ngogang, J.; Nyambi, P.N.; Duerr, R. Near Full Genome Sequencing of Novel HIV-1 Unique Recombinant Forms Circulating in Cameroon. J. Int AIDS Soc. 2019. Submitted. [Google Scholar]
  28. Courtney, C.R.; Agyingi, L.; Fokou, A.; Christie, S.; Asaah, B.; Meli, J.; Ngai, J.; Hewlett, I.; Nyambi, P.N. Monitoring HIV-1 Group M Subtypes in Yaounde, Cameroon Reveals Broad Genetic Diversity and a Novel CRF02_AG/F2 Infection. AIDS Res. Hum. Retroviruses 2016, 32, 381–385. [Google Scholar] [CrossRef]
  29. Agyingi, L.; Mayr, L.M.; Kinge, T.; Orock, G.E.; Ngai, J.; Asaah, B.; Mpoame, M.; Hewlett, I.; Nyambi, P. The evolution of HIV-1 group M genetic variability in Southern Cameroon is characterized by several emerging recombinant forms of CRF02_AG and viruses with drug resistance mutations. J. Med. Virol. 2014, 86, 385–393. [Google Scholar] [CrossRef]
  30. Ragupathy, V.; Zhao, J.; Wood, O.; Tang, S.; Lee, S.; Nyambi, P.; Hewlett, I. Identification of new, emerging HIV-1 unique recombinant forms and drug resistant viruses circulating in Cameroon. Virol. J. 2011, 8, 185. [Google Scholar] [CrossRef]
  31. Salazar-Gonzalez, J.F.; Bailes, E.; Pham, K.T.; Salazar, M.G.; Guffey, M.B.; Keele, B.F.; Derdeyn, C.A.; Farmer, P.; Hunter, E.; Allen, S.; et al. Deciphering human immunodeficiency virus type 1 transmission and early envelope diversification by single-genome amplification and sequencing. J. Virol 2008, 82, 3952–3970. [Google Scholar] [CrossRef] [PubMed]
  32. Tongo, M.; Dorfman, J.R.; Abrahams, M.R.; Mpoudi-Ngole, E.; Burgers, W.A.; Martin, D.P. Near full-length HIV type 1M genomic sequences from Cameroon: Evidence of early diverging under-sampled lineages in the country. Evol. Med. Public Health 2015, 2015, 254–265. [Google Scholar] [CrossRef] [PubMed]
  33. Ewing, B.; Hillier, L.; Wendl, M.C.; Green, P. Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998, 8, 175–185. [Google Scholar] [CrossRef]
  34. Li, H.; Durbin, R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 2010, 26, 589–595. [Google Scholar] [CrossRef] [PubMed]
  35. Robinson, J.T.; Thorvaldsdottir, H.; Winckler, W.; Guttman, M.; Lander, E.S.; Getz, G.; Mesirov, J.P. Integrative genomics viewer. Nat. Biotechnol. 2011, 29, 24–26. [Google Scholar] [CrossRef] [PubMed]
  36. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef]
  37. Tamura, K.; Dudley, J.; Nei, M.; Kumar, S. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol. Biol. Evol. 2007, 24, 1596–1599. [Google Scholar] [CrossRef] [PubMed]
  38. Courtney, C.R.; Mayr, L.; Nanfack, A.J.; Banin, A.N.; Tuen, M.; Pan, R.; Jiang, X.; Kong, X.P.; Kirkpatrick, A.R.; Bruno, D.; et al. Contrasting antibody responses to intrasubtype superinfection with CRF02_AG. PLoS ONE 2017, 12, e0173705. [Google Scholar] [CrossRef]
  39. Mayr, L.M.; Powell, R.L.; Ngai, J.N.; Takang, W.A.; Nadas, A.; Nyambi, P.N. Superinfection by discordant subtypes of HIV-1 does not enhance the neutralizing antibody response against autologous virus. PLoS ONE 2012, 7, e38989. [Google Scholar] [CrossRef]
  40. Powell, R.L.; Urbanski, M.M.; Burda, S.; Nanfack, A.; Kinge, T.; Nyambi, P.N. Utility of the heteroduplex assay (HDA) as a simple and cost-effective tool for the identification of HIV type 1 dual infections in resource-limited settings. AIDS Res. Hum. Retroviruses 2008, 24, 100–105. [Google Scholar] [CrossRef] [PubMed]
  41. Nanfack, A.J.; Redd, A.D.; Bimela, J.S.; Ncham, G.; Achem, E.; Banin, A.N.; Kirkpatrick, A.R.; Porcella, S.F.; Agyingi, L.A.; Meli, J.; et al. Multimethod Longitudinal HIV Drug Resistance Analysis in Antiretroviral-Therapy-Naive Patients. J. Clin. Microbiol 2017, 55, 2785–2800. [Google Scholar] [CrossRef] [PubMed]
  42. Alidjinou, E.K.; Deldalle, J.; Hallaert, C.; Robineau, O.; Ajana, F.; Choisy, P.; Hober, D.; Bocket, L. RNA and DNA Sanger sequencing versus next-generation sequencing for HIV-1 drug resistance testing in treatment-naive patients. J. Antimicrob. Chemother. 2017, 72, 2823–2830. [Google Scholar] [CrossRef] [PubMed]
  43. Jain, M.; Olsen, H.E.; Paten, B.; Akeson, M. The Oxford Nanopore MinION: Delivery of nanopore sequencing to the genomics community. Genome Biol. 2016, 17, 239. [Google Scholar] [CrossRef] [PubMed]
  44. Laver, T.; Harrison, J.; O’Neill, P.A.; Moore, K.; Farbos, A.; Paszkiewicz, K.; Studholme, D.J. Assessing the performance of the Oxford Nanopore Technologies MinION. Biomol. Detect. Quantif. 2015, 3, 1–8. [Google Scholar] [CrossRef] [PubMed]
  45. Tyler, A.D.; Mataseje, L.; Urfano, C.J.; Schmidt, L.; Antonation, K.S.; Mulvey, M.R.; Corbett, C.R. Evaluation of Oxford Nanopore’s MinION Sequencing Device for Microbial Whole Genome Sequencing Applications. Sci. Rep. 2018, 8, 10931. [Google Scholar] [CrossRef] [PubMed]
  46. Aralaguppe, S.G.; Siddik, A.B.; Manickam, A.; Ambikan, A.T.; Kumar, M.M.; Fernandes, S.J.; Amogne, W.; Bangaruswamy, D.K.; Hanna, L.E.; Sonnerborg, A.; et al. Multiplexed next-generation sequencing and de novo assembly to obtain near full-length HIV-1 genome from plasma virus. J. Virol. Methods 2016, 236, 98–104. [Google Scholar] [CrossRef]
  47. Jordan, M.R.; Kearney, M.; Palmer, S.; Shao, W.; Maldarelli, F.; Coakley, E.P.; Chappey, C.; Wanke, C.; Coffin, J.M. Comparison of standard PCR/cloning to single genome sequencing for analysis of HIV-1 populations. J. Virol. Methods 2010, 168, 114–120. [Google Scholar] [CrossRef] [PubMed]
  48. Lu, J.; Sista, P.; Giguel, F.; Greenberg, M.; Kuritzkes, D.R. Relative replicative fitness of human immunodeficiency virus type 1 mutants resistant to enfuvirtide (T-20). J. Virol. 2004, 78, 4628–4637. [Google Scholar] [CrossRef] [PubMed]
  49. Parrish, N.F.; Wilen, C.B.; Banks, L.B.; Iyer, S.S.; Pfaff, J.M.; Salazar-Gonzalez, J.F.; Salazar, M.G.; Decker, J.M.; Parrish, E.H.; Berg, A.; et al. Transmitted/founder and chronic subtype C HIV-1 use CD4 and CCR5 receptors with equal efficiency and are not inhibited by blocking the integrin alpha4beta7. PLoS Pathog. 2012, 8, e1002686. [Google Scholar] [CrossRef]
  50. Weber, J.; Vazquez, A.C.; Winner, D.; Rose, J.D.; Wylie, D.; Rhea, A.M.; Henry, K.; Pappas, J.; Wright, A.; Mohamed, N.; et al. 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. 2011, 55, 3729–3742. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Near full genome sequence amplification approaches. (A) Genome map of HIV-1 reference strain HxB2 (GenBank: K03455). (B) One-amplicon approach to obtain a single 9.1 kb PCR fragment (near full genome) shown in green. (C) Two-amplicon approach with two overlapping PCR fragments (shown in red), with half genome 1 (HG1) 5.2 kb and half genome 2 (HG2) 4.6 kb in size. (D) Multiple-amplicon approach: Selection of eight overlapping PCR fragments (shown in black), i.e. gag (1.5 kb), gagpol (4.5 kb), pol (3.2 kb), vifvpu (1.5 kb), vif/gp120 (3 kb), env (3.4 kb), gp120 (1.7 kb) and gp41/nef (1.8 kb). Positions of primers are schematically shown with gray arrows in (A) and (B).
Figure 1. Near full genome sequence amplification approaches. (A) Genome map of HIV-1 reference strain HxB2 (GenBank: K03455). (B) One-amplicon approach to obtain a single 9.1 kb PCR fragment (near full genome) shown in green. (C) Two-amplicon approach with two overlapping PCR fragments (shown in red), with half genome 1 (HG1) 5.2 kb and half genome 2 (HG2) 4.6 kb in size. (D) Multiple-amplicon approach: Selection of eight overlapping PCR fragments (shown in black), i.e. gag (1.5 kb), gagpol (4.5 kb), pol (3.2 kb), vifvpu (1.5 kb), vif/gp120 (3 kb), env (3.4 kb), gp120 (1.7 kb) and gp41/nef (1.8 kb). Positions of primers are schematically shown with gray arrows in (A) and (B).
Viruses 11 00317 g001
Figure 2. Conserved primer binding sites for HIV-1 near full genome amplification. Multiple sequence alignments of primer binding sites across ten different pure (sub)-subtypes and CRFs, as relevant for the near full genome amplification approach 1. For each (sub)-subtype/CRF, four reference sequences were selected that broadly cover the respective clade; for the binding site of primer 3’UTR rev 1, only two A2 and F2 reference sequences were available from the LANL database. The numbering of primer binding regions is based on the HxB2 reference strain (GenBank: K03455). Positions (pos) of primer binding sites are shown in brackets (HxB2 pos. 596–9542). Orange and green arrows along the HxB2 genome map indicate first round and second round primer positions, respectively. Black dots and colored letters indicate matches and mismatches with primer nucleotides shown on top, respectively.
Figure 2. Conserved primer binding sites for HIV-1 near full genome amplification. Multiple sequence alignments of primer binding sites across ten different pure (sub)-subtypes and CRFs, as relevant for the near full genome amplification approach 1. For each (sub)-subtype/CRF, four reference sequences were selected that broadly cover the respective clade; for the binding site of primer 3’UTR rev 1, only two A2 and F2 reference sequences were available from the LANL database. The numbering of primer binding regions is based on the HxB2 reference strain (GenBank: K03455). Positions (pos) of primer binding sites are shown in brackets (HxB2 pos. 596–9542). Orange and green arrows along the HxB2 genome map indicate first round and second round primer positions, respectively. Black dots and colored letters indicate matches and mismatches with primer nucleotides shown on top, respectively.
Viruses 11 00317 g002
Figure 3. Gel electrophoresis of PCR products of different lengths for composite near full genome analysis. (A) HIV-1 near full genome (NFG; green) and half genome (HG; red) amplicons obtained with different primer sets. (B) Diverse amplicons within half genome 1 (HG1). (C) Diverse amplicons within half genome 2 (HG2). Amplicons were separated on an ethidium bromide-stained 0.8% agarose gel. MassRuler DNA ladder was used for amplicon size determination. The amplicons shown in the gels were obtained from the following samples: NYU6541-6 (NFG1), NYU129-5 (NFG2), LB089-1 (NFG3 and env), NYU1122-1 (NFG4, HG1, and HG1a), NYU119-3 (HG2, pol, and gag), NYU1999 (gagpol), NYU6556-3 (HG1b), BDHS33 (vif/gp120), NYU124-2 (gp120), and MDC131-1 (gp41/nef and vifvpu).
Figure 3. Gel electrophoresis of PCR products of different lengths for composite near full genome analysis. (A) HIV-1 near full genome (NFG; green) and half genome (HG; red) amplicons obtained with different primer sets. (B) Diverse amplicons within half genome 1 (HG1). (C) Diverse amplicons within half genome 2 (HG2). Amplicons were separated on an ethidium bromide-stained 0.8% agarose gel. MassRuler DNA ladder was used for amplicon size determination. The amplicons shown in the gels were obtained from the following samples: NYU6541-6 (NFG1), NYU129-5 (NFG2), LB089-1 (NFG3 and env), NYU1122-1 (NFG4, HG1, and HG1a), NYU119-3 (HG2, pol, and gag), NYU1999 (gagpol), NYU6556-3 (HG1b), BDHS33 (vif/gp120), NYU124-2 (gp120), and MDC131-1 (gp41/nef and vifvpu).
Viruses 11 00317 g003
Figure 4. Development of a near full genome sequencing approach using recombinant and pure subtype strains. Genome maps of four HIV-1 strains, which were used to develop the near full genome sequencing approach. The determined subtypes are indicated using the recombinant drawing tool of the LANL database. NYU6541_6, NYU6506_8, LB002_1 and strain 0526 correspond to a unique recombinant form (URF), the circulating recombinant form CRF02_AG, sub-subtype F2, and subtype B, respectively.
Figure 4. Development of a near full genome sequencing approach using recombinant and pure subtype strains. Genome maps of four HIV-1 strains, which were used to develop the near full genome sequencing approach. The determined subtypes are indicated using the recombinant drawing tool of the LANL database. NYU6541_6, NYU6506_8, LB002_1 and strain 0526 correspond to a unique recombinant form (URF), the circulating recombinant form CRF02_AG, sub-subtype F2, and subtype B, respectively.
Viruses 11 00317 g004
Figure 5. Comparative recombinant breakpoint analysis of near full genome sequences obtained with four different amplification, sequencing, and cloning strategies. Using sample NYU6541_6, different amplification, sequencing, and cloning techniques were studied: (A) third-generation sequencing (TGS), (B) bulk amplification and Sanger sequencing, (C) cloning after bulk amplification, followed by Sanger sequencing, and (D) single genome amplification (SGA) and Sanger sequencing. BootScan plots (Simplot software) are shown for near full genome sequences of NYU6541_6 with window size 200 and step size 20. Breakpoint patterns were determined by analyzing the query sequence against CRF02_AG (red), F2 (cyan), C (dark gray), and B (light gray) reference strains. The genomic regions gag, pol, vif/vpr/vpu and env are schematically shown along the x-axis in green, red, yellow, and purple bars, respectively. The y-axis indicates bootstrap values obtained after 250 replicative measurements.
Figure 5. Comparative recombinant breakpoint analysis of near full genome sequences obtained with four different amplification, sequencing, and cloning strategies. Using sample NYU6541_6, different amplification, sequencing, and cloning techniques were studied: (A) third-generation sequencing (TGS), (B) bulk amplification and Sanger sequencing, (C) cloning after bulk amplification, followed by Sanger sequencing, and (D) single genome amplification (SGA) and Sanger sequencing. BootScan plots (Simplot software) are shown for near full genome sequences of NYU6541_6 with window size 200 and step size 20. Breakpoint patterns were determined by analyzing the query sequence against CRF02_AG (red), F2 (cyan), C (dark gray), and B (light gray) reference strains. The genomic regions gag, pol, vif/vpr/vpu and env are schematically shown along the x-axis in green, red, yellow, and purple bars, respectively. The y-axis indicates bootstrap values obtained after 250 replicative measurements.
Viruses 11 00317 g005
Table 1. Amplification primers.
Table 1. Amplification primers.
Portion AnalyzedNested PCR RoundPrimer IDSequencePosition in HxB2Temp.Reference
Near full genome 1 (NFG1)1st Round5’ UTR_For1TGCCTTGAGTGCTTCAAGTAGTGTG536–56058.4 °CThis study
3’ UTR_Rev1TATTGAGGCTTAAGCAGTGGGTTC9615–959257 °CThis study
2nd Round5’ UTR_For2TCTGTTGTGTGACTCTGGTAACTAGAGATC567–59658.7 °CThis study
3’ UTR_Rev4GCTCAAATCTGGTCTAGCAAGAGAGA9567–954258 °CThis study
Near full genome 2 (NFG2)1st RoundHIV_682_ForTCTCTCGACGCAGGACTCGGCTTGCTG682–70867 °C[25]
HIV_9555_RevTCTACCTAGAGAGACCCAGTACA9555–953355.5 °C[25]
2nd RoundHIV_776_ForCTAGAAGGAGAGAGAGATGGGTGCGAG776–80061 °C[25]
HIV_9555_RevTCTACCTAGAGAGACCCAGTACA9555–953355.5 °C[25]
Near full genome 3 (NFG3)1st RoundHIV_682_ForTCTCTCGACGCAGGACTCGGCTTGCTG682–70867 °C[25]
3’ UTR_Rev3AGAGCTCCCAGGCTCAAATCTGGTCTA9578–955262.5 °CThis study
2nd RoundHIV_776_ForCTAGAAGGAGAGAGAGATGGGTGCGAG776–80061 °C[25]
3’ UTR_Rev4GCTCAAATCTGGTCTAGCAAGAGAGA9567–954258 °CThis study
Near full genome 4 (NFG4)1st RoundHIV_682_ForTCTCTCGACGCAGGACTCGGCTTGCTG682–70867 °C[25]
3’ UTR_Rev5CAGTACAGGCGAGAAGCAGCTGCT9539–951662 °CThis Study
2nd RoundHIV_776_ForCTAGAAGGAGAGAGAGATGGGTGCGAG776–80061 °C[25]
3’ UTR_Rev6AGCAGCTGCTTATATGCAGCATCTGAG9525–949960.7 °CThis Study
Half genome 1a (HG1a)1st RoundHIV_682_ForTCTCTCGACGCAGGACTCGGCTTGCTG682–70867 °C[25]
Vpu2_RevCCGCTTCTTCCTGCCATAGGA5985–596659.1 °CThis study
2nd RoundHIV_776_ForCTAGAAGGAGAGAGAGATGGGTGCGAG776–80061 °C[25]
Vpu3_RevTCCTGCCATAGGAGATGCCTAAG5978–595658.1 °CThis study
Half genome 1b (HG1b)1st RoundGag3_ForGAGAGATGGGTGCGAGAGC785–80357 °CThis study
Vpu2_RevCCGCTTCTTCCTGCCATAGGA5985–596659.1 °CThis study
2nd RoundGag4_ForTAGTATGGGCAAGCAGGGA890–90856 °CThis study
Vpu3_RevTCCTGCCATAGGAGATGCCTAAG5978–595658.1 °CThis study
Half genome 2 (HG2)1st RoundVif2_ForTGGAAAGGTGAAGGGGCAGTA4956–497658.4 °CThis study
OFM19GCACTCAAGGCAAGCTTTATTGAGGCTTA9632–960460.9 °C[31]
2nd RoundVif3_ForGATTATGGAAAACAGATGGCAGGT5037–506055.2 °CThis study
HIV_9555_RevTCTACCTAGAGAGACCCAGTACA9555–953355.5 °C[25]
gagpol1st RoundGag2_ForGACTAGCGGAGGCTAGAAG764–78254 °CThis study
Pol2_RevCCATGTTCTAATCCTCATCCTGTC5103–508055 °CThis study
2nd RoundGag3_ForGAGAGATGGGTGCGAGAGC785–80357 °CThis study
Pol3_RevCTGTCTACCTGCCACACA5084–506754 °CThis study
gag1st RoundGag2_ForGACTAGCGGAGGCTAGAAG764–78254 °CThis study
Gag1_RevCCAATTCCCCCTATCAT2404–238848 °CThis study
2nd RoundGag4_ForTAGTATGGGCAAGCAGGGA890–90856 °CThis study
Gag3_RevGGTCGTTGCCAAAGAGTGA2278–226055 °CThis study
pol1st RoundPol1_ForGAAGAAATGATGACAGC1819–183545 °CThis Study
Pol1_RevTGCCAGTCTCTTTCTCCTG5279–516154 °CThis study
2nd RoundPol2_ForAAGTGTTTCAACTGTGG1960–197647 °CThis study
Pol2_RevCCATGTTCTAATCCTCATCCTGTC5103–508055 °CThis study
vifvpu1st RoundVif 1_ForGGGTTTATTACAGGGACAGCAGAG4900–492357.3 °C[31]
Vpu1_RevTTGCCACTYTCTTCTGCTCTTTC6225–620356.4 °CThis study
2nd RoundVif 2_ForTGGAAAGGTGAAGGGGCAGTA4956–497658.4 °CThis study
Vpu2_RevCCGCTTCTTCCTGCCATAGGA5985–596659.1 °CThis study
vif/gp1201st RoundVif1_ForGGGTTTATTACAGGGACAGCAGAG4900–492357.2 °CThis study
Gp120outGCARCCCCAAAKYCCTAGG8018–800057.2 °CThis study
2nd RoundVif2_ForTGGAAAGGTGAAGGGGCAGTA4956–497658.4 °CThis study
Gp120inCGTCAGCGTYATTGACGCYGC7838–781861.4 °CThis study
env1st RoundVif1_ForGGGTTTATTACAGGGACAGCAGAG4900–492357.2 °C[31]
OFM19GCACTCAAGGCAAGCTTTATTGAGGCTTA9632–960460.9 °C[31]
2nd RoundEnv A(CACC)GGCTTAGGCATCTCCTATGGCAGGAAGAA5954–598263 °C[31]
02AG-Env NGTTCTGCCAATCTGGGAAGAATCCTTGTGTG9174–914462.3 °CThis study
gp1201st RoundEnv AGGCTTAGGCATCTCCTATGGCAGGAAGAA5954–598262.8 °C[31]
Gp120outGCARCCCCAAAKYCCTAGG8018–800057.2 °C[39]
2nd RoundEnv BAGAAAGAGCAGAAGACAGTGGCA6202–622458.2 °C[31]
Gp120inCGTCAGCGTYATTGACGCYGC7838–781861.4 °C[39]
gp41/nef1st RoundGp120in_For1CAGCAGGAAGCACTATGGGCG7798–781860.7 °CThis study
3’UTR_Rev1TATTGAGGCTTAAGCAGTGGGTTC9615–959257 °CThis study
2nd RoundGp120in_For2GCRGCGTCAATRACGCTGACG7818–783861.4 °CThis study
3’UTR_Rev4GCTCAAATCTGGTCTAGCAAGAGAGA9567–954258 °CThis study
Table 2. Forward sequencing primers.
Table 2. Forward sequencing primers.
Primer IDSequencePosition in HxB2Temp.Reference
HIV_580_ForTCTGGTAACTAGAGATCC580–59746.6 °CThis study
HIV_625_ForATCTCTAGCAGTGGCGCCCGA625–64163 °CThis study
HIV_788_ForAGATGGGTGCGAGAGCGT788–80959.2 °CThis study
HIV_1250_ForCATGGGTAAAGGTAATAGAAG1250–127047.9 °CThis study
HIV_1250b_ForCATGGGTAAAAGTAATAGAA1250–126944.5 °CThis study
HIV_1400_ForCCATCAATGAGGAAGCTGCA1400–141955.4 °CThis study
HIV_1830_ForTGATGACAGCATGCCAGG1830–184755.3 °CThis study
HIV_1970_ForTTCAACTGTGGCAAAGAAGG1970–198953.4 °CThis study
HIV_2075_ForGACAGGCTAATTTTTTAGGGA2075–209549.9 °CThis study
HIV_2590_ForCAGGAATGGATGGCCCAA2590–260755.2 °CThis study
HIV_2700_ForGGGCCTGAAAATCCATACAATACT2700–272355.1 °CThis study
HIV_2756_ForGTACTAAATGGAGAAAATTAG2756–277643.8 °CThis study
HIV_3300_ForAGCTGGACTGTCAATGA3300–331650.4 °CThis study
HIV_3350_ForGGGCAAGTCAAATTTATCCAG3350–337051.8 °CThis study
HIV_3355_ForAGCCAGATTTATCCAGG3355–337148.2 °CThis study
HIV_4000_ForTAGCCTTGCAGGATTCAGGAT4000–402056.0 °CThis study
HIV_4175_ForTGGAGGAAATGAACAAGTAGA4175–419551 °CThis study
HIV_4542_ForGCAGGAAGATGGCCAGT4542–455855.0 °CThis study
HIV_4646_ForTGGAATTCCCTACAATCC4646–466348.6 °CThis study
HIV_4747_ForAGACAGCAGTACAGATGGCAG4747–476756.7 °CThis study
HIV_4900_ForGGGTTTATTACAGGGACAGCA4900–492054.6 °CThis study
HIV_5388_ForTTTCAGAATCTGCCATAAG5388–540647.1 °CThis study
HIV_5769_ForCATTTCAGAATYGGGTG5769–578646.5 °CThis study
HIV_5840_ForGTAGATCCTARCCTAGA5840–585644.1 °CThis study
HIV_5970_ForATGGCAGGAAGAAGCGGAGAC5970–599059.5 °CThis study
HIV_6100_ForAGTAGCATTCATAGCAGCCAT6100–612054.0 °CThis study
HIV_6125_ForGTGTGGACTATAGTATATATAG6125–614644.1 °CThis study
HIV_6380_ForATTTTGTGCATCAGATGC6380–639744.7 °CThis study
HIV_6543_ForGATATAATTAGTCTATGGG6543–656140.8 °CThis study
HIV_6745_ForCACTTTTTTATAGACTTGAT6745–676442.7 °CThis study
HIV_6826_ForTTAMACAGGCTTGTCC6826–684147.2 °CThis study
HIV_6840_ForCCAAAGGTATCCTTTGAGCCA6840–686054.9 °CThis study
HIV_6855_ForGAGCCAATTCCCATACAT6855–687249.3 °CThis study
HIV_7500_ForATGTGGCAGAAAGTAGGACAAGC7500–752257.1 °CThis study
HIV_7520_ForAGCAATGTATGCCCCTC7520–753652 °CThis study
HIV_7633_ForCTGGAGGAGGAGATATGAG7633–765150.7 °CThis study
HIV_7660_ForGGAGAAGTGAATTATATAA7660–767840.6 °CThis study
HIV_7803_ForGGAAGCACTATGGGCGC7803–781956.7 °CThis study
HIV_8800_ForGGTGGCAAGTGGTCAAA8800–881653.0 °CThis study
HIV_8180_ForGCAGGAAAAGAATGAACAAG8180–819949.7 °CThis study
HIV_8970_ForGGTTAGAAGCACAAGAG8970–898647.0 °CThis study
HIV_9015_ForGGTACCTTTAAGACCAATGA9015–903449.3 °CThis study
HIV_9030_ForGACCAATGACTTATAAGG9030–904744.1 °CThis study
Table 3. Reverse Sequencing Primer.
Table 3. Reverse Sequencing Primer.
Primer IDSequencePosition in HxB2Temp.Reference
HIV_809_RevACGCTCTCGCACCCATCT809–79258.9 °CThis Study
HIV_912_RevTCCCTGCTTGCCCATACTA912–89455.6 °CThis study
HIV_1270_RevCTTCTATTACTTTTACCCATG1270–125045.9 °CThis study
HIV_1505_RevGTTCCTGCTATRTCACTTCC1505–148651.5 °CThis study
HIV_2095_RevTCCCTAAAAAATTAGCCTGTC2095–207550.2 °CThis study
HIV_3577_RevTGATAAATTTGATATGTCCA3577–355843.9 °CThis study
HIV_5060_RevACCTGCCATCTGTTTTCCATA5060–504054.2 °CThis study
HIV_5537_RevACACTAGGCAAAGGYGG5537–552153.5 °CThis study
HIV_6352_RevGGTACCCCATAATAGACTGTRACCCACAA6352–632459.9 °C[25]
HIV_6764_RevATCAAGTCTATAAAAAAGTG6764–674542.7 °CThis study
HIV_6841_RevGGACAAGCCTGTKTAA6841–682647.2 °CThis study
HIV_6867_RevTGGGAATTGGCTCAAA6867–685248.1 °CThis study
HIV_6904_RevTTTAGAATCGCAAAACCAGC6904–688551.4 °CThis study
HIV_7025_RevTTCTGCTAGRCTGCCATT7025–700852.4 °CThis study
HIV_7085_RevCTGTACTATTATGGTTT7085–706939.1 °CThis study
HIV_7351_RevAAACTATGTGTTGTAATTTC7351–733243.5 °CThis study
HIV_7819_RevGCGCCCATAGTGCTTCC7819–780356.7 °CThis study
HIV_7845_RevCCYGTACCGTCAGCGT7845–783056 °CThis study
HIV_8075_RevTTTTTACACCGTCTAG8075–805942 °CThis study
HIV_8448_RevTGTCTTGCTCKCCACCT8448–843242 °CThis study
HIV_8528_RevGTAGCTGAAGAGGCACAG8528-851154.9 °CThis study
HIV_8529_RevGGTAGCTGAAGAGGCACAGG8529-851052.6 °CThis study
REV14ACCATGTTATTTTTCCACATGTTAAA6526–650154.6 °C[31]
REV15CTGCCATTTAACAGCAGTTGAGTTGA7015–699057.8 °C[31]
REV 16ATGGGAGGGGCATACATTGCT7540–752059 °C[31]
REV 17CCTGGAGCTGTTTAATGCCCCAGAC7956–793261.7 °C[31]
REV 18GGTGAGTATCCCTGCCTAACTCTAT8365–834157.2 °C[31]
REV 19ACTTTTTGACCACTTGCCACCCAT8820–879759.6 °C[31]
Table 4. Summary of amplification approaches and yield using 23 pure or recombinant HIV-1 samples.
Table 4. Summary of amplification approaches and yield using 23 pure or recombinant HIV-1 samples.
Near Full
Genome 1
Near Full
Genome 2
Near Full
Genome 3
Near Full
Genome 4
Half
Genome 1a
Half
Genome 1b
Half
Genome 2
gagpolvif/gp120envgp120
1st Rd Primers5’UTRfor1/
3’UTRrev1
HIV682for/
HIV9555rev
HIV682for/
3’UTRrev3
HIV682for/
3’UTRrev5
HIV682for/
Vpu2 rev
Gag3for/
Vpu2 rev
Vif2for/
OFM19
Gag2for/
pol2 rev
Vif1for/
Gp120out
Vif1/
OFM19
EnvA/
Gp120out
Sample ID2nd Rd Primers5’UTRfor2/
3’UTR rev4
HIV776for/
HIV9555rev
HIV776for/
3’UTR rev4
HIV776for/
3’UTR rev6
HIV776for/
Vpu3rev
Gag4for/
Vpu3 rev
Vif3for/
HIV9555rev
Gag3for/
pol3rev
Vif2for/
Gp120in
EnvA/
02AG-EnvN
EnvB/
Gp120in
LB016-1positivepositivepositivepositivepositivepositivepositive
LB069-1positivepositivepositivepositivepositivepositive
LB082-1positivepositivepositivepositivepositivepositive
LB089-1positivepositivepositivepositivepositivepositivepositivepositive
LB095-1positivepositivepositivepositivepositivepositivepositive
LB104-1positivepositivepositivepositive
MDC131-1positivepositivepositive
MDC179-2positivepositivepositivepositive
BDHS24-2positivepositivepositivepositivepositivepositive
BDHS33positivepositivepositivepositive
NYU119-3positivepositivepositivepositivepositive
NYU124-2positivepositivepositivepositivepositivepositive
NYU129-5positivepositivepositivepositivepositivepositive
NYU1122-1positivepositivepositivepositivepositivepositivepositive
NYU1999positivepositivepositivepositive
NYU2140-1positivepositivepositivepositivepositive
NYU6556-3positivepositivepositivepositivepositivepositive
NYU6541-6positivepositivepositivepositivepositivepositivepositivepositivepositivepositive
NYU6506-8positivepositive
Strain0526positivepositive
LB002-1positiveNANANANANANANA
LB022-1positiveNANANANANANANA
LB006-1positiveNANANANANANANA
Miscellaneous HIV-1 genome regions amplified
gagpolvifvpugp120gp41/nef
1st Rd primersGag2for/
Gag1 rev
Pol1for/
Pol1rev
Vif1for/
Vpu1rev
EnvA/
Gp120out
Gp120in for1/
3’UTRrev1
Sample ID2nd Rd primersGag4for/
Gag3rev
Pol2for/
Pol2rev
Vif2for/
Vpu2rev
EnvB/
Gp120in
Gp120in for2/
3’UTRrev4
MDC131-1positivepositivepositiveNApositive
NYU119-3positivepositivepositivepositivepositive
NYU124-2positivepositivepositivepositivepositive
NA: not applicable; −: negative.
Table 5. Summary of third-generation sequencing (TGS) and single-genome amplification (SGA).
Table 5. Summary of third-generation sequencing (TGS) and single-genome amplification (SGA).
#Sample IDTGS
Total Number of HIV Reads
TGS
Average Read Length
TGS
Number of Long HIV Reads
TGS
Genomic Region
SGA
Genomic Region
1LB016-1713173117 (>5000 bp)Near full genome 1NA
2LB069-117013929488 (>5000 bp)Near full genome 1NA
3LB082-16104232179 (>4000 bp)Half genome 2Half genome 2
4LB089-1NANANANANA
5LB095-1NANANANAHalf genome 2
6LB104-11120383061 (>4000 bp)Half genome 2NA
7MDC131-1NANANANANA
8MDC179-219562807504 (>2000 bp)vif/gp120NA
9BDHS024-234634047799 (>4000 bp)Half genome 2NA
10BDHS3331844122813 (>4000 bp)Half genome 2NA
11NYU119-3NANANANANA
12NYU124-220762837964 (>2000 bp)vif/gp120NA
13NYU129-5NANANANAHalf genome 2
14NYU1122-126482478752 (>2000 bp)vif/gp120NA
15NYU1999-1NANANANANA
16NYU2140-121084044448 (>4000 bp)Half genome 2NA
17NYU6556-3NANANANAHalf genome 2
18NYU6541-6262141541 (>5000 bp)Near full genome 1Near full genome 1
19NYU6506-8NANANANANA
20Strain0526NANANANAHalf genome 2
21LB002-1NANANANANA
22LB022-1NANANANANA
23LB006-1NANANANANA
NA: not applicable.

Share and Cite

MDPI and ACS Style

Banin, A.N.; Tuen, M.; Bimela, J.S.; Tongo, M.; Zappile, P.; Khodadadi-Jamayran, A.; Nanfack, A.J.; Meli, J.; Wang, X.; Mbanya, D.; et al. Development of a Versatile, Near Full Genome Amplification and Sequencing Approach for a Broad Variety of HIV-1 Group M Variants. Viruses 2019, 11, 317. https://doi.org/10.3390/v11040317

AMA Style

Banin AN, Tuen M, Bimela JS, Tongo M, Zappile P, Khodadadi-Jamayran A, Nanfack AJ, Meli J, Wang X, Mbanya D, et al. Development of a Versatile, Near Full Genome Amplification and Sequencing Approach for a Broad Variety of HIV-1 Group M Variants. Viruses. 2019; 11(4):317. https://doi.org/10.3390/v11040317

Chicago/Turabian Style

Banin, Andrew N., Michael Tuen, Jude S. Bimela, Marcel Tongo, Paul Zappile, Alireza Khodadadi-Jamayran, Aubin J. Nanfack, Josephine Meli, Xiaohong Wang, Dora Mbanya, and et al. 2019. "Development of a Versatile, Near Full Genome Amplification and Sequencing Approach for a Broad Variety of HIV-1 Group M Variants" Viruses 11, no. 4: 317. https://doi.org/10.3390/v11040317

APA Style

Banin, A. N., Tuen, M., Bimela, J. S., Tongo, M., Zappile, P., Khodadadi-Jamayran, A., Nanfack, A. J., Meli, J., Wang, X., Mbanya, D., Ngogang, J., Heguy, A., Nyambi, P. N., Fokunang, C., & Duerr, R. (2019). Development of a Versatile, Near Full Genome Amplification and Sequencing Approach for a Broad Variety of HIV-1 Group M Variants. Viruses, 11(4), 317. https://doi.org/10.3390/v11040317

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop