Optimization of a Low Cost and Broadly Sensitive Genotyping Assay for HIV-1 Drug Resistance Surveillance and Monitoring in Resource-Limited Settings

Commercially available HIV-1 drug resistance (HIVDR) genotyping assays are expensive and have limitations in detecting non-B subtypes and circulating recombinant forms that are co-circulating in resource-limited settings (RLS). This study aimed to optimize a low cost and broadly sensitive in-house assay in detecting HIVDR mutations in the protease (PR) and reverse transcriptase (RT) regions of pol gene. The overall plasma genotyping sensitivity was 95.8% (N = 96). Compared to the original in-house assay and two commercially available genotyping systems, TRUGENE® and ViroSeq®, the optimized in-house assay showed a nucleotide sequence concordance of 99.3%, 99.6% and 99.1%, respectively. The optimized in-house assay was more sensitive in detecting mixture bases than the original in-house (N = 87, P<0.001) and TRUGENE® and ViroSeq® assays. When the optimized in-house assay was applied to genotype samples collected for HIVDR surveys (N = 230), all 72 (100%) plasma and 69 (95.8%) of the matched dried blood spots (DBS) in the Vietnam transmitted HIVDR survey were genotyped and nucleotide sequence concordance was 98.8%; Testing of treatment-experienced patient plasmas with viral load (VL) ≥ and <3 log10 copies/ml from the Nigeria and Malawi surveys yielded 100% (N = 46) and 78.6% (N = 14) genotyping rates, respectively. Furthermore, all 18 matched DBS stored at room temperature from the Nigeria survey were genotyped. Phylogenetic analysis of the 236 sequences revealed that 43.6% were CRF01_AE, 25.9% subtype C, 13.1% CRF02_AG, 5.1% subtype G, 4.2% subtype B, 2.5% subtype A, 2.1% each subtype F and unclassifiable, 0.4% each CRF06_CPX, CRF07_BC and CRF09_CPX. Conclusions The optimized in-house assay is broadly sensitive in genotyping HIV-1 group M viral strains and more sensitive than the original in-house, TRUGENE® and ViroSeq® in detecting mixed viral populations. The broad sensitivity and substantial reagent cost saving make this assay more accessible for RLS where HIVDR surveillance is recommended to minimize the development and transmission of HIVDR.


Introduction
Treatment of HIV-1 infection with highly active antiretroviral therapy (HAART) in the past decades has remarkably reduced HIV/ AIDS related mortality and morbidity. However, the emergence of drug resistance in persons on antiretroviral therapy (ART) and the transmission of drug-resistant HIV strains to newly infected persons are a major threat to the global success on HIV prevention and treatment effort [1,2,3]. Recent years, under multilateral supports for HIV treatment and prevention programs, especially the U.S. President Emergency Plan for AIDS Relief (PEPFAR) with the targets of treating two million HIV-infected people with ART, preventing five million new HIV infection and care for 10 million HIV-infected people and AIDS orphans, access to antiretroviral drugs (ARVs) has been scaled up rapidly in resource-limited countries where availability of laboratory monitoring is often limited or lacking [4,5]. This creates the potential for HIV drug resistance (HIVDR) emergence and transmission in these settings. Detection and monitoring of HIVDR by molecular genotyping is pivotal to ensure ongoing regimen efficacy. It is the standard of care in resource-rich countries [2,6]; however in resource-limited countries, HIVDR testing is not generally available or it is too costly to be used in routine monitoring of patients receiving ARVs. Therefore, the World Health Organization (WHO) recommends population-based surveillance and monitoring of HIVDR in resource-limited settings [2,4,7]. Pattern and rates of transmitted and acquired drug resistant HIV variants will collectively inform regional and global recommendations on which ARVs to maintain or change in first and second-line regimens [7].
Population sequencing-based genotyping methods including ViroSeqH, TRUGENEH and in-house assays are widely used, and the most informative and affordable genotyping methods for monitoring patients on ART in clinical practice [8,9,10,11]. However, ViroSeqH and TRUGENEH, the two FDA-approved genotyping assays were designed for HIV-1 group M subtype B viruses which are the predominant HIV-1 strains in resource-rich countries. In addition, these commercial kits are expensive and less sensitive to non-B subtypes, limiting their utility in resource-limited settings [12,13,14]. There have been no commercially available HIV-1 genotyping assays designed for non-B subtypes and circulating recombinant forms (CRFs) that are predominant viral strains in resource-limited countries. Moreover, the demand for low cost and sensitive genotyping methods is increasing with the establishment and expansion of laboratory molecular monitoring in these settings [15,16].
The most frequently used HIVDR genotypic assays are assays that detect resistance mutations in the reverse-transcriptase (RT) and protease (PR) genes [17,18,19]. The minimal genotyping requirements for these two regions are PR codons 10-99 and RT codons   [3,20,21]. Our original in-house assay [22] has its limitations: (1). It does not cover the entire PR gene region required for resistance testing; (2). For some HIV-1 subtypes or CRFs, some sequencing primers generate higher background noises which would affect the detection of mixture bases. We redesigned two primers for RT-PCR, one for nested PCR and four sequencing primers in the optimized assay.
In-house assays are relatively inexpensive and sensitive for multiple subtypes, but in-house assays should only be implemented after adequate validation, including evaluating assay's performance with various HIV-1 subtypes and CRFs [3,20,23]. Factors that could contribute to genotyping quality include type of assay/ kit used, specimen handling and storage, level of experience of technicians performing the analysis, heterozygosity of sequences, and viral subtypes in clinical samples [11,21]. In this study, therefore, we validated sensitivity, accuracy and specificity of the optimized in-house assay. We also evaluated the assay applicability to HIVDR surveillance and monitoring.

Samples
A total of 381 samples were used in this study: 151 samples were used for validation and 230 for application. For validation, we included 111 HIV-1 positive plasma and 10 dried blood spot (DBS) samples, as well as 30 HIV-1 negative DBS samples. For application, we tested 132 plasma and 98 DBS samples. The detailed information on these samples is described in Table 1.

Dried blood spot sample preparation and storage
Dried blood spot samples were prepared by spotting 100 ml of whole blood onto each of the five preprinted circles on a Whatman 903 filter paper (Whatman Inc, Piscataway, NJ) and were then dried overnight at ambient temperature. The next day, a piece of glassine paper was folded and a DBS card was placed in the folded paper, and 10-20 wrapped DBS cards were then packaged in a Bitran bag containing desiccant sacks and a humidity indicator card and sealed. The packaged DBS cards were stored either at 270uC for Vietnam samples or at room temperature for an average of 85.31642.66 days (median 83.5 days) for Nigeria samples. The DBS cards were shipped to the WHO Specialized Drug Resistance Laboratory at the Centers for Disease Control and Prevention (CDC) (Atlanta, GA, USA) either on dry ice for the samples from Vietnam or at ambient temperature for the DBS samples from Nigeria. All samples were stored at 270uC upon the arrival at CDC.

Viral RNA and total nucleic acid extraction
The QIAamp mini-viral RNA kit (Qiagen, Valencia, CA) was used to extract RNA from all plasma samples for validation purpose. Details for viral load (VL) measurement on samples from Cameroon, Thailand and Zambia were described previously  [22,24,25]. For Malawi samples, Abbott m2000 automatic sample preparation system (0.2 ml extraction protocol) was used to extract the plasma RNA. For Nigeria and Vietnam plasma and DBS samples, the NucliSensH EasyMAG TM automatic sample preparation system (BioMérieux, Durbam, NC) was used to extract the plasma RNA and blood total nucleic acid (TNA). To extract TNA, one DBS spot was cut out per specimen and placed in a 2ml of NucliSENSH lysis buffer (Biomeriuex, Durham, NC) for 30 min at room temperature with gentle rotation. Nucleic acid was then extracted from DBS samples using the NucliSENSH EasyMagH automated extraction system following the manufacturer's instructions. Nucleic acid was eluted in 25 ml of NucliSENSH Extraction Buffer 3 and stored at 280uC until use.

RT-PCR and nested PCR
Sequences of RT-PCR and sequencing primers that were redesigned or modified based on the original assay [22] and HIV-1 sequences available at the Los Alamos HIV Database (www.hiv. lanl.gov) are shown in Table 2. Three RT-PCR and four sequencing primers were replaced during the optimization process. All primers used were synthesized at CDC Biotechnology Core Facility. Two oligonucleotides that mixed at a ratio of 1:1 (w/w) were used as the forward primer for one step RT-PCR. One-step RT-PCR was performed in a 50 ml reaction, which consisted of 10 ml of RNA or TNA extracts, 0.16 mM each of primers PRTM-F1 and RT-R1, and 0.5 ml SuperScript TM III one step RT/ PlatinumH Taq high Fidelity Enzyme Mix and 1x reaction buffer mixture containing Mg2+ and eoxyribonucleotide triphosphates (dNTPs) (Invitrogen, Carlsbad, CA). RT-PCR condition was an initial cycle RT step at 50uC for 45 min and 94uC for 2 min, and followed by 40 cycles of PCR at 94uC for 15 sec, 50uC for 20 sec, 72uC for 2 min and an extension at 72uC for 10 min. For nested PCR, 2 ml of RT-PCR product was added to a 50 ml reaction containing 0.12 mM of each of the inner primers PRT-F2 and RT-R2, 1x GeneAmp Gold Buffer II, 2 mM MgCl 2 , 400 mM each dNTP and 2.5 U of AmpliTaq Gold LD DNA polymerase (Applied Biosystems, Foster City, CA). After initial denaturation at 94uC for 4 min, 40 cycles of PCR were performed in GeneAmp 9700 thermocycler with the PCR conditions as 94uC for 15 sec, 55uC for 20 sec and 72 for 2 min and following an extension at 72uC for 10 min. In the case of the failed first RT-PCR attempt, PRTM2-F1 was used as rescue primer to replace PRTM-F1 to account for mutations occurring within the primer binding site.
The nested PCR product was confirmed by 1% agarose gel electrophoresis with a product size of 1,084 base pairs. The confirmed PCR products were purified using Exo-SAP IT PCR purification kits and used for cycle sequencing reaction with BigDye terminator cycle sequencing kit 3.1 (Applied Biosystems, CA).

Sequence analysis
DNA sequencing of HIV-1 pol was performed in 3730 DNA genetic analyzer (Applied Biosystems). Six sequencing primers overlapping the entire amplicon were used ( Table 2). Sequencing raw data were edited with ChromasPro, v1.5 (Technelysium Pty Ltd, Australia) and confirmed by a second technician. To double check for all mixed bases, we also used a web-based sequence analysis tool, ReCall [26], in which minor peak calling was set at 15% of the main peak. To rule out PCR contamination, phylogenetic analyses were performed on all newly obtained sequences by MEGA 4 [27]. Sequence quality was also checked by Stanford HIVdb program. Sequences with frame shifts or stop codons were excluded from analysis. For transmitted HIVDR surveillance, WHO surveillance drug resistance mutation (SDRM) list was used [28]. For HIVDR monitoring surveys, drug resistance-associated mutations in PR and RT were interpreted using the Stanford Genotypic Resistance Interpretation Algorithm (http://hivdb.stanford.edu/pages/algs/HIVdb.html). Pairwise nucleo tide sequence identity and discrepancy were analyzed using BioEdit [29].
Sensitivity, accuracy and specificity of the assay As for HIV-1 drug resistance genotyping, there are no standardized or reference method (gold standard) to evaluate analytic and clinical performance in molecular genotyping for HIV-1 group M viruses. We validated the new method according to WHO/HIVResNet drug resistance guidelines [21], including participation in an external quality assessment (EQA) program, proficiency testing (PT) panels and comparing the results between new method and the original method already established in our laboratory [22].
Because the validation criteria were difficult to define based on the complexity of samples tested in this multi-subtype evaluation and all currently available assays (commercial or in-house) were unable to genotype 100% of the samples tested [16], for this study we used the genotyping sensitivity intervals as $95% for samples with VL $3 log copies/ml; accuracy was defined as detection of 99% of known DR mutation codons, and reproducibility/ precision was defined as $98% nucleotide identities in $ 90% of pairwise comparisons. In this study, sensitivity and reproducibility of the assay were assessed by comparing the current genotyping results from 96 field collected samples of known VL with those of the original assay [22] from Cameroon, Zambia and Thailand. The assay was also evaluated using TRUGENEH system GL12 with plasma from DigitalPT (N = 5), a HIVDR PT program offered by AccuTest at Boston, MA, USA and using ViroSeqH system v 2.8 with plasma from Virology Quality Assurance program (VQA, N = 5), a WHO-sponsored HIIVDR PT program and offered by VQA at Chicago, IL, USA. Additionally, the precision of the assay was evaluated using 4 replicates of a second VQA plasma PT panel (N = 5) and 3 of 4 replicates were tested with the optimized in-house assay and the remaining one was tested with TRUGENEH. The precision test was performed by 3 technicians. In addition, 10 DBS panels shipped at two different temperature conditions from VQA were also tested by the optimized assay. Specificity was determined by testing 30 HIV-negative DBS specimens collected from pregnant women in Tanzania.
Applying the assay for surveillance of transmitted HIVDR and HIVDR prevention monitoring surveys in resourcelimited countries Seventy two matched plasma and DBS samples from newly HIVdiagnosed persons in Vietnam were tested. For HIVDR monitoring surveys, we applied the optimized assay for resistance testing in samples collected from patients 12-15 months after the commencement of ART in two monitoring surveys. For the Malawi monitoring survey, 34 plasma samples from patients with VL ranged from 2.25 to 5.89 log10 copies/ml were tested. In Nigeria monitoring survey, 26 matched plasma and DBS samples with plasma VL ranged from 2.18 to 6.41 log10 copies/ml were analyzed.

HIV-1 subtyping
HIV-1 subtyping for the newly obtained sequences was performed using the REGA 4 HIV-1 Genotyping Tool [30]. Phylogenetic analyses were further conducted using neighborjoining method included in the MEGA 4 for sequences with unclassifiable subtypes. Reference sequences were obtained from the Los Alamos HIV Database (www.hiv.lanl.gov). Sequences obtained in the study were submitted to GenBank and their accession numbers are JN885633 to JN885719.

Reagent cost comparison
To estimate reagent cost savings by using the optimized broad sensitive genotyping assay, we calculated reagent cost per test of the assay and compared it with the reagent costs of commercially available genotyping systems, TRUGENEH and ViroSeqH. We used current U.S. market values in dollars for all the reagents we used in the in-house assay including RNA/TNA extraction, RT-PCR, nested PCR, PCR amplification confirmation, PCR purification and sequencing reactions. These reagent cost estimates did not include the cost for running test controls and any repetitions when needed.

Statistical analysis
Wilcoxon Signed-Rank test was used to analyze the difference in number of nucleotide mixtures detected between the optimized and original in-house assays. The statistical significance was considered when P value was ,0.05.

Ethical consideration
In accordance with United States regulations and international guidelines, the CDC human subjects review process determined this activity to be non-research and the protocol was approved by the Associate Director for Science (ADS), National Center for HIV, Hepatitis, STD and TB Prevention, CDC, Atlanta. All the study protocols were approved by local institutional review boards:

Validation of the optimized in-house assay
Sensitivity. The sensitivity of the optimized in-house assay was evaluated with 96 HIV-1 positive plasma samples collected from Cameroon, Thailand and Zambia. Of these, all 5 samples with VL ,3 log10 copies/ml and 87 (95.6%) of 91 samples with VL $3 log10 copies/ml were genotyped, resulting in an overall genotyping rate of 95.8% (92/96) comparing to 96.8% (93/96) by the original assay.
Accuracy. The accuracy of the optimized assay was first assessed by comparing 87 paired nucleotide sequences generated by the original in-house assay [22] and the optimized assay using ReCall and BioEdit programs. The mean nucleotide identity was 99.360.50% (mean 6 SD) among paired nucleotide sequences. Wilcoxon signed-rank test was used to compare original and optimized in-house assays in basecalling for mixed bases and revealed that the optimized assay detected significantly more mixed bases than the original one (P,0.001). However, this difference did not translate into differences in HIVDR mutations. Among 144 DR mutations detected in paired samples, we did not detect any complete discordant mutations at DR mutation sites and only 11 partially discordant DR mutation sites including 3 in PR and 8 in RT (Table 3) were seen. The overall DR codon agreement was 99.8% between the 87 paired samples.
Testing 10 plasma PT panel samples (five from DigitalPT and five from VQA) using the optimized in-house and TRUGENEH or ViroSeqH assays also indicated that the optimized in-house assay appeared to detect more mixed bases than commercial kits. However, DR mutation site differences only occurred in 2 of 76 DR mutations between the optimized in-house and TRUGENEH and 2 of 44 between the in-house and ViroSeqH in mixed bases. The overall sequence identity was 99.660.41% between the inhouse and TRUGENEH (Table 4), and 99.160.65% between the in-house and ViroSeqH (Table 5). Further examination to see whether more sensitive detection of base mixtures in the optimized in-house assay is a reproducible event, we analyzed 4 replicates of 5 samples that were tested by three optimized in-house assay runs and one TRUGENEH run under different operators. We found highly concordant sequence identities ranging from 98.22% to 99.65%. The minor differences observed in sequence identity were caused by base mixtures (Table 6). For example, at codons 37 and 41 of RT in sample 3 (Fig 1), one replicate detected mixture RY at the 2 nd and 3 rd positions of codon 37, and the second replicate did not find any mixtures, while the third replicate showed lower, yet visible, second peaks comparing to the first replicate. Similarly at the 2 nd position of codon 41, the first replicate revealed a Y (C/T), the second replicate showed a W (A/T), and the third replicate revealed an H (A/C/T) while TRUGENEH replicate detected a Y (C/T). Nucleotide mixtures also caused some mismatched DR mutations between the replicates. For instance, one in-house replicate missed 3 mixed codons (K65KR, D67DN and T69IT) in sample 2 among the 4 replicates. However, the minor peaks of nucleotide bases at these three codons could be seen, but were below the mixture cutoff (15%) on the chromatogram by ReCall. Thus these mixtures were not counted and resulted in the codon discrepancy. Another partial discordant example was the DR mutation M184MV, which was detected in sample 5 by all 3 inhouse replicates but not found in TRUGENEH replicate. These results indicated that the sequences generated by population-based sequencing were highly reproducible but the sensitivity at detecting low frequency of drug resistant HIV variants was very challenging. Table 3. Discordant drug resistance-associated amino acid positions in protease and reverse transcriptase from 87 plasma samples a genotyped by the original and optimized in-house assays.   Sensitivity on DBS samples. DBS samples are recommended by WHO for HIVDR surveillance in resource-limited settings in treatment-naïve populations [21]. The mission of our laboratory at CDC is to support HIVDR surveillance in PEPFARsupported countries; thus, we evaluated the assay sensitivity with two matched DBS experimental panels shipped under different temperature conditions from VQA. The optimized in-house assay was able to genotype all 5 DBS samples shipped with dry ice and 4 of 5 DBS samples shipped at ambient temperature with DBS VL ranging from 3.17 to 3.98 log10 copies/ml. The failed sample GEN001BS.04C was the sample with the lowest VL of 3.17 log10 copies/ml (f 7).
To evaluate the specificity of the optimized in-house assay, we tested HIV negative DBS samples (N = 30) collected from women attending ANC clinics in Tanzania and they were all found to be negative, resulting in the assay specificity of 100%.
EQA assessment results: Based on WHO/HIVResNet requirement to pass the PT panels, a drug resistance mutation (DRM) site score and nucleotide (nt) alignment score with consensus sequence of at least 99% (considering all 5 samples) must be achieved. The optimized in-house assay passed two sets of plasma PT panels with 100% DRM, 99.98% nt and 100% DRM, 99.88% nt scores, respectively. The assay showed high sequence concordance with the laboratories participating in the EQA program.

Application of the assay in the surveillance of HIVDR in resource-poor countries
With the satisfactory validation results of the optimized in-house assay, we applied this assay in the surveillance of HIVDR in three PEPFAR-supported countries.

Threshold survey of transmitted HIVDR in recently HIVinfected population in Vietnam.
For this survey, we tested 72 plasma and matched DBS specimens collected from individuals attending voluntary counseling and testing (VCT) in Ho Chi Minh City. We were able to genotype all 72 plasma and 69 (95.8%) DBS samples and sequence identity analysis (N = 69) indicated that overall nucleotide identity was 98.9%60.62% between matched plasma and DBS samples. The sequence differences were caused by partially discordant mixture bases located in the three HIVDR codons in the RT region.
Detection of HIVDR development in ART-experienced patients from Malawi and Nigeria. We next utilized the optimized in-house assay in detecting HIVDR development in patients treated with first-line ARVs for 12-15 months in two monitoring surveys conducted in Malawi (N = 34) and Nigeria (N = 26). Genotyping was successful for all 46 plasma samples collected from virologically failed patients defined as VL $3 log10 copies/ml according to the WHO definition [4]. For patients with VL between 2.18 and ,3 log10 copies/ml, 78.6% (11/14) plasma samples (7/8 from Nigeria and 4/6 from Malawi) were also successful. Furthermore, all 18 matched DBS samples from virologically failing patients and 4 of the 8 DBS samples with VL between 2.18 and ,3 log10 copies/ml from Nigeria were genotyped. The nucleotide sequence identity between the 22 plasma and DBS pairs was 98.860.80%. For DR mutations, 90.4% DR mutations identified in plasma were also found in DBS. For the 9.6% discordant DR mutations identified, majority of them was partially discordant (7.5%) and only 2.1% were completely discordant.

Reagent cost comparison
Using the current U.S. market values in dollars for all the reagents used in the optimized in-house assay, we estimated that the reagent cost per test for the optimized in-house assay was $40.00, comparing to $213.20 for TRUGENEH and $172.86 for ViroSeqH. In the reagent cost calculations, we did not include the cost for assay controls and any need for repetition of tests, which would increase the cost of reagents for all the assays compared here.

Discussion
The newly optimized in-house assay was broadly sensitive in genotyping multiple HIV-1 group M subtypes and CRFs from plasma and DBS collected from 6 resource-limited countries. The original in-house assay, although a success from the broad sensitivity perspective, was in need for improvement due to the concern of incomplete genotyping of PR gene and suboptimal Table 6. Genotyping reproducibility of replicate PCR products generated from independent RT-PCR amplification process by 3 different operators in a 5-member proficiency testing panel received from VQA. sequence quality due to primer design [22]. The validation of the re-designed RT-PCR and some of the sequencing primers in the current study confirmed that the newly optimized in-house assay is comparable to the original in-house assay in assay sensitivity and specificity and it is also broadly sensitive to all group M subtypes and CRFs circulating in PEPFAR-supported countries. Pairwise nucleotide sequence identity analyses from sequences generated by the optimized in-house assay and the ones obtained from the original in-house assay and two commercially available genotyping systems indicated the optimized in-house assay produced comparable genotyping results. More importantly, the optimized inhouse assay expanded genotyping codon coverage to include all PR mutations and improved sequence quality by reducing background noises to minimal, resulting in more sensitive mixture calling. The ability to detect nucleotide mixtures (low frequency viral strains) is important as recent studies have demonstrated that low frequency variants can grow rapidly and become predominant viral population under the selective drug pressure and lead to treatment failure [31,32]. Given that the newly optimized in-house assay was broadly sensitive in genotyping B and non-B subtype viral strains of HIV-1 group M viruses, it would be expected that the assay would efficiently genotype plasma and DBS samples from various geographical areas, and this was verified by applying the assay in genotyping samples collected from patients enrolled in two perspective HIVDR monitoring surveys. The optimized in-house assay was able to genotype 100% of plasma samples collected from virological failure patients, defined as VL $3 log10 copies/ml [4] in the present study and over three-fourths of the patients with VL between 2.18 and ,3 log 10 copies/ml. More importantly, genotyping was successful for all the 18 matched DBS samples collected from Nigerian patients with virological failure. It is worthy to note that these DBS samples had been stored at room temperature for an average of 85 days before shipping to our laboratory for testing. In addition, testing of two DBS PT panels shipped frozen or at ambient temperature revealed that genotyping was successful for all DBS samples except one with the lowest VL and shipped at ambient temperature. These results indicate that the optimized in-house assay is highly sensitive in genotyping both plasma and DBS samples. It is important to point out that interpretation of the genotyping results from the DBS samples stored at room temperature for a long period of time however, needs to be cautious since these DBS cards were packaged correctly and stored in an air-conditioning room with low Table 7. Genotyping efficiency and drug resistance-associated mutations identified in protease (PR) and reverse transcriptase (RT) by the optimized in-house assay from dried blood spots (DBS) PT panels.  humidity. These package and storage condition might have limited the true impact of suboptimal storage conditions existed in resource-limited settings on the quality of DBS cards and resulted in better genotyping efficiency. Studies have shown that correct packaging and storage of DBS are critical elements in ensuring successful genotyping results [33,34,35,36]. In fact, our own data here generated from two identical DBS PT panels and shipped at frozen or ambient temperature also indicate that even overnight exposure of DBS samples with low VL to ambient temperature in domestic shipment could have some detrimental impact on DBS quality for genotyping. Comparing the performance of the optimized in-house assay with TRUGENEH, ViroSeqH and the original in-house assay, high nucleotide sequence identity was revealed; however, minor differences existed in mixture base callings. The optimized in-house assay detected more mixed bases than the commercial kits and our original in-house assay. Many factors could contribute to the sequence discordances at the mixed nucleotide sites in HIV genotyping including viral quasispecies, primer binding preference and location, Taq polymerase misincorporation, sequence quality, basecalling criteria or technical errors [11,37,38]. Because HIV-1 viruses are rapidly evolving quasispecies [39], there are multiple HIV-1 variants in one patient [37,40]. Sequence identity and codon concordance are challenging when mixed bases are present [41]. It has been reported that ViroSeqH detected more mixtures (78%) than an in-house assay (22%) [10]. In contrast, our optimized in-house assay detected more mixture bases than other assays. This may be due to the fact that the optimized assay produced sequence chromatographs containing minimal background noise. To confirm this, we performed sequence editing for all validation samples (N = 102) including PT panels with ReCall program [26] using minor peak default mixture calling setting at .15% of the major peak in bidirectional sequences. We also independently tested one set of the PT panels by 3 different operators. These analyses showed that the optimized in-house assay gave more sensitive mixture calling. The variability in detecting nucleotide mixtures was likely due to the first-round RT-PCR [42] in sampling of quasispecies strains rather than by technical errors in the sequencing process [43]. The use of wide-spectrum degenerate primers and a mixture of two forward primers at slightly different binding sites in the initial run of RT-PCR are likely contributed to more mixtures calling in the optimized in-house assay. To what extent and by what factors the mixture variants could be affected and detected in HIV genotyping is a matter of speculation, which needs further studies. Phylogenetic analyses indicated that the optimized in-house assay could genotype HIV-1 group M subtypes A (A1, A2), B, C, F (F1, F2), G and CRFs including CRF01_AE, CRF02_AG, CRF06_CPX, CRF07_BC, CRF09_CPX and UC with an overall sensitivity of 96% using specimens from different geographical regions around the world. Due to the limited availability of HIVpositive samples, we only genotyped a small number of subtype A and F viral strains and further studies are needed to confirm our findings with larger sample sizes on these viral strains. It has been reported that genotyping sensitivity with two FDA-approved systems using non-B subtypes varies [44]. Some studies indicated that these two systems performed well for B and non-B subtypes [9,45,46,47] while others demonstrated that they were less sensitive to non-B subtypes and CRFs [12,13,14,48]. For instance, only 52% of serum samples were genotyped in an Ethiopian threshold survey using ViroSeqH and TRUGENEH methods sequentially [49]. For genotyping DBS samples collected from subtype B infected persons, one study reported 78.8% genotyping rate by TRUGENEH [50] while another study reported 57.5% genotyping rate using ViroSeqH for DBS samples stored for one year at 4uC [51]. An additional study reported an even lower DBS genotyping rate of 38.6% by ViroSeqH system [40]. However, a study using ViroSeqH documented 100% genotyping rate for DBS samples collected from subtype B-infected patients with VL great than 2,000 copies/ml and 54.5% genotyping rate with DBS samples from patients with VL less than 2,000 copies/ml when DBS samples were stored at optimal conditions [52]. Thus, commercial genotyping systems might work well with DBS samples collected from subtype-B-infected patients. Genotyping of DBS samples using these genotyping systems in non-B subtypes needs to be further studied. In comparison to our original in-house assay, our own experience using these two commercial assays with non-B subtypes was also not satisfactory. They often required repetition of RT-PCR or sequencing due to failure to amplify or sequence in TRUGENEH and ViroSeqH assays [12,14]. Compared to these commercial assays, the optimized in-house assay was not only sensitive, but also inexpensive. The assay could reduce the cost for genotyping reagents by 75%. The availability of low cost and broadly sensitive genotyping assay for plasma and DBS would make HIVDR surveillance and monitoring in resource-limited settings more accessible.
In conclusion, we have validated and improved a broadly sensitive and less expensive in-house genotyping assay for HIVDR surveillance and monitoring in resource-limited countries. Validation analyses indicate that the optimized in-house assay detected more mixed HIV-1 population than our original in-house assay and commercial genotyping kits. Given the high efficiency in genotyping diverse HIV-1 group M viral strains from plasma and DBS samples and substantial reagent cost saving, the optimized inhouse assay could be applicable to DR genotyping in both ARTnaive and -experienced populations according to current WHO recommendations for surveillance purpose [4].