Victoria A. Johnson,4,5 John W. Mellors,6 John L. Sullivan,7 Ronald Swanstrom,8 Richard T. D'Aquila,9,
Mark Van Gorder,10 Mark Holodniy,11 Robert M. Lloyd, Jr.,12,
Caroline Reid,12 Gillian F. Morgan,12 and Dean L. Winslow12
Gladstone Institute of Virology and Immunology,1 University of California, San Francisco,2 Consolidated Laboratories, Van Nuys,10 AIDS Research Center, VA Palo Alto Health Care System and Division of Infectious Diseases and Geographic Medicine, Stanford University, Palo Alto, California,11 Harvard University Medical School,9 University of Massachusetts, Boston, Massachusetts,7 University of North Carolina at Chapel Hill, Chapel Hill, North Carolina,8 University of Alabama at Birmingham,4 Veteran's Administration Medical Center, Birmingham, Alabama,5 University of Pittsburgh, Pittsburgh, Pennsylvania,6 Division of Infectious Diseases, University of Colorado Health Sciences Center, Denver, Colorado,3 and Visible Genetics, Inc., Toronto, Ontario, Canada,12
Received 11 February 2002/ Returned for modification 13 October 2002/ Accepted 21 December 2002
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The clinical trials demonstrating utility of drug resistance genotyping for the selection of antiretroviral therapy were supported by centralized reference laboratories with extensive experience with viral genotypic testing. In clinical practice, drug resistance genotypic testing is performed using a wide variety of methods whose performance characteristics have not been thoroughly established. Significant variations exist between laboratories in the extent of technical training, the reagents used, the amplification and sequencing protocols, and datum review and reporting procedures. This variation in methods may contribute to the variable quality of resistance assay results that have been observed in international evaluations (24).
We now describe a multicenter blinded study of the accuracy and reproducibility of results obtained with the TRUGENE HIV-1 Genotyping Kit and OpenGene DNA Sequencing System (Visible Genetics, Inc.). This is an integrated system that includes a training and certification program, positive and negative control specimens, reverse transcription (RT)-PCR amplification, dye primer labeled bi-directional sequencing, polyacrylamide gel electrophoresis, and sequence analysis including sequence contiguous element assembly, guided manual editing, and an interpretative HIV-1 resistance report based on algorithms developed and periodically revised by an international committee of experts. The study aimed to evaluate the accuracy of the assay system in the evaluation of plasma samples from HIV-1-infected persons and virus stocks derived from infectious molecular clones and to determine the variability in accuracy by laboratory site, by technician, by day, and by kit lot.
| MATERIALS AND METHODS |
|---|
|
|
|---|
"Gold standard" sequence determination of plasma from HIV-1-infected persons. Two aliquots of 1 ml each of plasma from each subject were centrifuged at 45,000 x g for 1 h at 4°C to concentrate viral particles, which were then resuspended in 150 µl of plasma. Viral RNA from the entire volume was mixed in lysis buffer then isolated using columns containing activated silica. The RNA was eluted in 50 µl according to the manufacturer's recommendations (QIAamp viral RNA kit; Qiagen Inc.). To allow broader sampling of the virus population, two aliquots of the extracted RNA (17 µl each) were independently amplified by RT-PCR using two distinct and different primer sets, neither of which are used in the TRUGENE HIV-1 Genotyping Kit. These extraction and amplification conditions are sufficient to consistently amplify specimens with viral loads as low as 60 copies per ml (R Lloyd, R. Schuurman, H. Stang, DeGroot, L. Hough, D. Burns, R. Mathis, and P. Feorino, 3rd Int. Workshop HIV Drug Resist. Treat. Strat., abstr. 52, 1999). PCR products were purified by agarose gel electrophoresis and extraction from the gel (QIAQUICK gel extraction kit; Qiagen Inc.) and then ligated into pGEM-T vectors (Promega). The ligation products were introduced into Escherichia coli by transformation followed by drug selection. The presence of inserted patient-derived gene inserts was verified by restriction digestion of color-selected colonies. Purified plasmid DNA was sequenced using Cy5.0 and Cy5.5 labeled T7 and SP6 primers (5'-GTAATACGACTCACTATAGGG-3' and 5'-ATTTAGGTGACACTATAGAATAC-3, respectively). Twenty cloned sequences from each specimen were derived by using the first 10 sequences from each of two independent PCR amplification reactions. The consensus sequence from aligned nucleotide sequences was used as the gold standard sequence from each specimen. In order to represent mixtures greater than or equal to 30% of the plasma HIV-1 population, bases reported in at least 6 of the 20 clones (30%) were included in the gold standard consensus sequence using standard ambiguity codes. Viral sequences in one of the nine plasma specimens (VA-MH-001) were difficult to clone using standard cloning techniques; therefore, only 16 sequences were used to establish the gold standard consensus in this instance.
Generation of virus stocks from infectious molecular clones. Mutant HIV-1 viruses containing 72 different drug resistance mutations were prepared from infectious molecular clones of HIV-1. Mutant variants of HIV-1 were prepared by site-directed mutagenesis of an infectious molecular clone of HIV-1LAI (pLAI.2) (kindly provided by L. Montagnier through the NIH AIDS Research and Reference Reagent Program). Mutations were introduced into the protease (PR)- or RT-coding region of pol using the Altered Sites II in vitro mutagenesis system (Promega, Madison, Wis.). The resulting plasmids were linearized by restriction enzyme digestion and electroporated into CEM-SS cells to generate infectious viruses, which were propagated to yield high-titer stocks. HIV-1 viral load was assayed (Roche Amplicor HIV-1 Monitor), and the virus stocks were diluted to the desired concentration in blood plasma from anti-HIV-1 seronegative donors and stored at -70°C with no temperature fluctuation greater than 10%.
Performance sites. The laboratory analyses of specimens were performed using the TRUGENE HIV-1 Genotyping Kit assay in eight laboratories in the United States; six laboratories were involved in the analysis of plasma samples from infected persons and three laboratories were involved in the analysis of virus samples prepared from infectious molecular clones. Samples were distributed among the sites in a balanced manner to maximize the power to detect differences in assay performance between laboratory sites. All technologists participating in the study were trained and certified to perform the TRUGENE HIV-1 Genotyping Kit assay using standardized protocols. The training involved 1 week of practical and theoretical training at a reference laboratory (Visible Genetics, Inc., Atlanta, Ga.) followed by proficiency testing of a blinded practice panel and a blinded certification panel of four unknown specimens plus two controls. Technologists were required to achieve greater than 95% sequencing accuracy on the blinded panels.
TRUGENE HIV-1 Genotyping Kit performance. All assay procedures were performed according to the manufacturer's recommendations (Visible Genetics, Inc.). Briefly, plasma viral RNA was extracted using activated silica in a column format (QIAamp viral RNA kit; Qiagen Inc.). For this study, 140 µl of unconcentrated plasma was mixed in a lysis buffer, placed over the column, washed, and then eluted in 60 µl. An aliquot of 17 µl of extracted RNA eluate was mixed with RT-PCR reagents containing RNase inhibitor, Taq polymerase, sequence specific primers, dATP, dCTP, dGTP, dTTP, magnesium, dithiothreitol, and pH buffers. The RT-PCR primers used in the kit amplify the entire protease and the first 250 codons of reverse transcriptase as a single amplicon. The reverse transcriptase reaction was performed at 90°C for 2 min, followed by incubation at 50°C for 5 min, at which point reverse transcriptase enzyme was added to the reaction mixture, which was then incubated for an additional 55 min before the reaction was terminated at 94°C for 2 min. DNA amplification was performed in the same tube using cycling conditions of 94°C for 30 s, 57°C for 30 s, and 68°C for 2 min for 20 cycles; this was followed by cycling conditions of 94°C for 30 s, 60°C for 30 s, and 68°C for 2.5 min for 17 cycles; and the cycling finished with one cycle of 68°C for 7 min and then a 4°C hold. Unpurified RT-PCR product (5 µl) from each specimen was used in each of 16 sequencing reactions, based on the CLIP principle (25), and employing four pairs of primers that sequence the protease reading frame (two pairs) and the beginning and middle of the reverse transcriptase reading frame (one pair each). For each segment, four sequencing reactions were used with four dideoxy terminators (ddATP, ddCTP, ddGTP, and ddTTP), Cy5-labeled forward primers, and Cy5.5-labeled reverse primers to generate CLIP sequencing products over 30 reaction cycles (94°C for 20 s, 56°C for 20 s, 70°C for 1.5 min) followed by a 5-min terminal extension at 70°C. Unpurified CLIP sequencing products were analyzed using 6% polyacrylamide-urea-Tris-borate-EDTA gel electrophoresis for 50 min at 2,000 V, with premanufactured Microcel cassettes, which were cast using Surefill 6% sequencing polyacrylamide gel cartridges (Visible Genetics, Inc.), containing premeasured amounts of acrylamide, urea, and Tris-borate-EDTA, which were polymerized using UV light.
Sequence data analysis. Sequence data acquisition and analysis were performed using the Gene Objects software (Visible Genetics, Inc.). The software allows automated and real-time acquisition of sequence data from the DNA sequencers (Long Read Tower Sequencers). Sequence fragments are assembled automatically and electropherogram patterns are displayed for semiautomated review after alignment with a clade B HIV-1 reference sequence (HIV-1LAI) and a mutated reference sequence containing codon changes typically observed in drug-resistant HIV-1. The use of both mutant and wild-type reference sequences serves to minimize biases toward the wild type during manual editing. The software can be set to skip automatically to positions where there are differences in base identification between the forward and reverse sequencing reactions or between the determined sequence and either of the reference strains or if results from only one sequencing direction is available. For this study, manual editing was limited to bases that were discordant between the two sequencing reactions. The edited sequences were exported from the software for analysis.
Statistical analysis. All statistical analysis was performed by an independent contractor (HHI, Hunt Valley, Md.), who had no communication with any of the laboratory testing sites. The sequences derived from each assay were compared to the gold standard sequence to assess the extent of agreement. Agreement was compared for all nucleotide positions, for all codons, for codons that were mutated in the gold standard (relative to HIV-1LAI), for codons that were wild type in the gold standard, and for codons that have been associated with decreased susceptibility or poor virological response to antiretroviral therapies. Only exact agreement was scored as an accurate result, defined as detection of all of the bases present in the gold standard sequences and no other bases. For example, if the gold standard sequence indicated that A was the predominant base present at a given site and G was present in 6 of 20 clones, the result would be scored as accurate only if A and G were detected and C and T were not detected. Detection of C would be considered inaccurate in the a priori data analysis plan even if C were detected in at least one of the gold standard sequences, but at a frequency of less than 30% of the clones. The mutation recognition rate was defined as the portion of mutated codons (relative to HIV-1LAI) that were exactly detected by the assay. The wild-type recognition rate was defined as the proportion of nonmutant codons that were exactly detected by the assay. Codons associated with antiretroviral resistance were defined as codons that have been associated with virological drug failure in vivo, decreased drug susceptibility in vitro, or drug selection in tissue culture, as compiled in published tables (6, 14). Differences in accuracy between laboratories, technicians, kit lots, and days were evaluated for statistical significance using analysis of variance across all groups, followed by pairwise comparisons if there was evidence of differences. A secondary analysis, designed after the primary analysis, allowed results to be considered accurate if bases or codons present in the gold standard in fewer than 30% of clones were detected, provided that all bases present in more than 30% of the gold standard clones were also detected. The accuracy of replicate testing of virus stocks derived from infectious molecular clones was calculated for each of the 72 drug resistance mutations present in the panel.
Nucleotide sequence accession numbers. The consensus sequences from the gold standard sequencing procedures have been submitted to GenBank with accession numbers AF472535 to AF472552.
| RESULTS |
|---|
|
|
|---|
|
|
Overall accuracy in clinical plasma samples. The overall accuracy of the base detection was 98.7% counting only exact matches, defined as the detection of all of the bases present in at least 30% of the gold standard clones and no other bases (Table 3) (30% rule). Detection at sequence positions with pure wild-type codons in the gold standard clones was more accurate than detection at positions with some evidence of mutation, which frequently represented mixtures of bases (0.99 versus 0.82 for wild-type and mutant detection, respectively). A large portion of results deemed inaccurate using the 30% rule were due to detection of sequence variants that were present in less than 30% of the gold standard clones, which was the cutoff for detection established a priori for the data analysis plan. A secondary analysis was performed which allowed detection of bases or codons to be considered accurate provided that the base or codon was present in at least one of the gold standard clones (e.g., 1 of 20 or 5%), and all of the bases or codons present in at least 30% of the gold standard clones were detected in the assay (Table 3, 5% rule). This secondary analysis revealed substantially higher accuracy using the 5% rule to define the gold standard. For example, accuracy of mutation detection improved from 0.82 to 0.94 using the 5% rule. Wild-type mutation recognition agreement remained higher than the mutation recognition agreement rate, as expected, because wild-type codons by definition had no variation in the sequences used to define the gold standard, while mutant codons could represent any mixture of mutant and wild-type codons. There was no difference in assay accuracy at sequence positions that are known to be important for drug resistance and all other positions.
|
Accuracy in virus samples prepared from infectious molecular clones. An analysis was performed to assess accuracy of detection of the entire spectrum of mutations associated with drug resistance in virus stocks derived from infectious molecular clones of HIV-1. Ten HIV-1 stocks derived from infectious molecular clones were spiked into plasma from HIV-1-seronegative donors and diluted to a concentration that just exceeded 10,000 HIV-1 RNA copies/ml. The infectious molecular clones were constructed to carry multiple antiretroviral resistance-associated mutations (Tables 1 and 2) and were sequenced 10 times each using the TRUGENE HIV-1 Genotyping Kit. Taken together, the panel contained 10 virus stocks containing 72 different drug resistance mutations. In a blinded fashion, the panel was tested six times in one laboratory and two times in each of two other laboratories. Fifty-four of 72 resistance mutations were correctly identified in all 10 replicate tests, and the remaining 18 mutated codons were correctly identified in >90% of test results. Fingerprint analysis of nucleotide sequences revealed that results from two aliquots of two specimens from one of the laboratory sites appeared to have been contaminated with sequences from another HIV-1 isolate, unrelated to the virus in the panel (data not shown). When results of these two contaminated sequences were excluded from the analysis, 31 of 31 (100%) resistance mutations in PR and 40 of 41 (97.6%) resistance mutations in RT were correctly identified in all of the remaining assay results; one RT mutation (K101E) was correctly identified in 9 of 10 (90%) of the test results.
Analysis of codon misidentifications in plasma specimens from infected persons. All assay results (n = 90) that were used to compare assay performance between six laboratory sites were also used for a detailed analysis of the clinical implications of codon misidentifications that occurred in the course of this study. In this analysis, we considered 54 codon positions that have been associated with drug resistance, virological drug failure, or drug selection (6). Of the 4,860 drug resistance codon positions evaluated (i.e., 90 assays x 54 codon positions per assay), assay results did not exactly match the gold standard at 70 (1.4%) positions, using the 5% rule to define the gold standard (Table 4). Of the 70 misidentified codon positions, 3 (4.3%) proved to be due to mislabeling of one specimen aliquot prior to distribution to one of the laboratories. Among the remaining 67 misidentified codon positions, the assay correctly identified at least one (but not all) of the codons present in the gold standard. In 61% (41 of 67) of the codon positions, the assay identified a codon that was not present in the gold standard, in addition to the gold standard codon. In the other 39% (26 of 67) of codon positions, the assay detected only one of several codons that were present in the gold standard. The majority (n = 55 of 67; 82%) of the codon misidentifications occurred at five codon positions (PR 10 and 63 and RT 184, 210, and 215) in three specimens (Table 4). Test results for these codon positions were identical in the majority of independent assays, indicating a high degree of concordance between laboratories despite discordance with the gold standard. The remaining 12 codon sites that were not exactly matched with the gold standard were distributed across eight sites in three specimens.
|
| DISCUSSION |
|---|
|
|
|---|
Several factors are likely to have contributed to the accuracy and reproducibility of this assay. First, the establishment of explicit training procedures and proficiency evaluation is important, especially for clinical laboratories where medical technologists may not have training or prior experience in molecular biology. Second, the detailed and explicit protocol for the assay contributes to consistent performance, which probably accounts for the improved performance of kit-based assays versus home brew genotyping assays in a recent international proficiency evaluation (W. Keulen, D. Brambilla, M. Buimer, J. Tijnagel, J. Bremer, S. Land, L. de Graaff, H. Versteeg, C. Boucher, and R. Schuurman, 5th Int. Workshop Drug Resist. Treat. Strat., abstr. 166, 2001; L. Shaker-Irwin, A. Scarsella, E. Rogolsky, J. Stryker, and S. Day, 5th Int. Workshop Drug Resist. Treat. Strat., abstr. 160, 2001). Third, the use of dye-labeled primers has been associated with more consistent detection of mixtures of sequences compared with dye-labeled terminators, which may have biased base incorporation rates depending on adjacent sequences. Fourth, the automated data acquisition, analysis, and reporting helps eliminate transcription errors which have accounted for a large portion of sequencing discrepancies in other studies (M. Hoover, D. Wentworth, J. Neaton, et al., 4th Int. Workshop HIV Drug Resist. Treat. Strat., abstr. 78, 2000). Fifth, strict adherence to good manufacturing practices during production of the of the kit, as well as hardware associated with the assay, may explain the lack of variation in accuracy from lot to lot. Finally, the TRUGENE HIV-1 Genotyping Kit has been simplified by elimination of DNA purification steps after the RT-PCR and the CLIP sequencing reactions, by use of the novel CLIP reaction and by use of preformed gel cassettes. These simplifications serve to make the assay more robust through minimization of assay steps.
The results of two assays of infectious molecular clones of HIV-1 were invalidated due to presence of a contaminating HIV-1 sequence. Following completion of data analysis, sequence contamination was investigated using the Genetic Fingerprint function of the OpenGene software. This software function compares the pattern of polymorphisms and mutations of a sequence to all other sequences in the software library's database. Sequences showing identity with another sequence in the database (other than those from the same patient) suggest the possibility of contamination and should be investigated prior to reporting the result. For this study of coded specimens, participating laboratories were instructed not to use the Fingerprint function, which could have biased accuracy results by revealing virus identities in the panel, thereby compromising blinding procedures. Use of the Fingerprint function is, however, recommended during routine laboratory use and was found to detect the contamination events that occurred in this study. Stringent precautions against contamination when performing PCR-based tests are needed (16), including architectural separation of laboratory areas used for PCR setup and DNA analysis, use of negative template controls, and sequence comparison with sequences determined previously in the same laboratory, as with fingerprinting. Uracil N-glycosylase is also used to decrease carry-over contamination in PCR-based assays, but was not included in the formulation of the TRUGENE assay which minimizes and detects contamination events by other methods. Further, uracil N-glycosylase treatment would not eliminate the risk of carryover contamination during the setup of CLIP reactions utilized in the TRUGENE assay.
Our study posed several real-life challenges to a drug resistance genotyping assay, including use of plasma from HIV-1-infected humans. Virus populations in vivo are diverse, which poses significant challenges to molecular diagnostic assays that involve oligonucleotide primers for amplification and sequence detection. These challenges are addressed in the TRUGENE HIV-1 Genotyping Kit through use of degenerate primers and ambiguity codes that allow information about sequence mixtures to be retained in the report. Nevertheless, all of the codon identification errors analyzed in Table 4 (other than those due to mislabeling of aliquots prior to distribution to testing sites) involved incomplete, but partially correct, identification of mixtures of codons present in either the gold standard consensus or the assay result. Even analysis of multiple viral sequence clones from two highly sensitive PCRs, as performed here to establish the gold standard sequence consensus, does not guarantee full representation of the diversity of clinically occurring HIV-1 populations. For example, the consistent detection of the RT L210W mutation (as a mixture with wild type) in 11 of 12 independent assays of specimen VA-MH-005 suggests that the mutation was present in the virus population at low frequency, although it was not detected in any of the 20 clones used to establish the gold standard. Failure to identify the mutation in the clones used for the gold standard could occur due to random sampling error of diverse populations. Ambiguity inherent in establishing gold standard measurements of naturally occurring virus populations is also indicated by the improvement in assay accuracy (Table 3) observed when using the 5% cutoff to define the gold standard, which allowed detection of codons present in less than 30% of the virus population to be considered correct. Mutation detection in homogeneous virus stocks derived from infectious molecular clones was also found to be highly accurate and reproducible, representing important measures of the assay performance that is not confounded by virus population diversity.
The challenges for training and quality control were also increased in our study because of the involvement of multiple laboratories (both academic and clinical), multiple testing personnel, and multiple kit lots. Use of large panels of specimens, unique labels for all aliquots, and an external data management and analysis agency eliminated any opportunities for partial unblinding of analysis or reporting in our study. Hence, we believe that the assay performance described here will reflect assay performance in a variety of clinical and research laboratories. Our estimates of assay performance cannot be directly compared with estimates based on less challenging evaluations that involve analysis of molecular clones of HIV-1, extensive use of reference laboratories, or proficiency panels that involve limited numbers of specimens and labeling that may allow laboratories to cross-validate results prior to reporting.
After the study was completed, the TRUGENE assay that is available to clinical laboratories was modified to include a centrifugation step to concentrate the virus from 1 ml of plasma prior to RNA extraction and provision of alternative PCR primers (version 1.5). The centrifugation step was found to decrease the minimum viral load required for genotypic evaluation to less than 100 copies/ml (R. Lloyd, R. Schuuman, H. Stang, DeGroot, L. Hough, D. Burns, R. Mathis, and P. Feorino, 3rd Int. Workshop HIV Drug Resist. Treat. Strat., abstr. 52, 1999). The alternative primers are designed to bind to diverse HIV-1 subtypes to permit testing of non-B HIV-1 subtypes using the TRUGENE system (L. Jagodzinksi, J. Cooley, S. Kelly, and N. Michael, 9th Conf. Retrovir. Opportun. Infect., abstr. 594-T, 2002). The version 1.5 primers are one of the primer sets used in this study to amplify the gold standard sequences. Both assay modifications aim to improve the assay success rate and did not affect the accuracy of sequence determinations in the clinical laboratory that participated in this study (data not shown).
The clinical utility of drug resistance genotyping assays has been consistently demonstrated in clinical trials where laboratory testing is performed in reference laboratories with extensive experience with molecular biology assays (2, 10). Translation of this experience from clinical trials to clinical practice will require the availability of assay technology that can be mastered and consistently controlled in clinical diagnostic laboratories. The results of this study indicate that the TRUGENE HIV-1 Genotyping Kit is an accurate and robust system suitable for clinical laboratory testing.
| ACKNOWLEDGMENTS |
|---|
We thank Jacqueline Javier, Birgit Drews, Darren Hazelwood, Jennifer Koel, Randy Huelsman, Pamela Beatty, Terri Watkins, Lorraine Sutton, Lijia Shi, Russell Young, Shing-Yi Lee, and Manjeet Sandhu for outstanding technical support. Statistical analysis and data management were performed by Clifford Meng (HHI). Anila Choudhry (VGI Clinical Research) provided clinical research support and monitoring of the study.
| FOOTNOTES |
|---|
Present address: Partners AIDS Research Center, Brigham and Womens Hospital, Cambridge, MA 02139. ![]()
Present address: Division of Infectious Disease, Department of Medicine, Vanderbilt University School of Medicine, Nashville, TN 37232. ![]()
Present address: Research Think Tank, Alpharetta, GA 30004. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Antimicrob. Agents Chemother. | Clin. Microbiol. Rev. |
|---|---|
| Clin. Vaccine Immunol. | ALL ASM JOURNALS |
|---|