Previous Article | Next Article ![]()
Journal of Clinical Microbiology, July 2003, p. 3306-3311, Vol. 41, No. 7
0095-1137/03/$08.00+0 DOI: 10.1128/JCM.41.7.3306-3311.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Department of Pediatrics and Laboratory of Biochemical Pharmacology, Emory University School of Medicine and Veterans Affairs Medical Center, Decatur, Georgia 30033
Received 9 December 2002/ Returned for modification 18 March 2003/ Accepted 11 April 2003
|
|
|---|
|
|
|---|
The introduction of automated DNA sequencing technology revolutionized the rapid sequencing of novel genes, identification of drug-selected changes in the HIV-1 genome, and the appreciation of mixtures of viral populations in clinical samples (11). Further, hybridization strategies, such as line probe assay (LiPA), that can discriminate between wild-type (WT) and mutated codon were also introduced (27). Recent studies have shown that HIV-LiPA is more sensitive than sequencing for the detection of mixed virus populations but can result in a high level of hybridization failures (22). In addition, LiPA in general is limited, since probes for certain important mutations, such as K65R are unavailable.
During the past few years, clinical trials have shown that multidrug resistance mutants selected with highly active antiretroviral therapy are generally attenuated for replication. As soon as the role of drug-selected reverse transcriptase (RT) and protease mutations affecting HIV-1 replication became evident (2, 3, 5, 6, 16, 24, 25), automated sequencing strategies were used for quantification of viral growth fitness by comparing the peak heights of two nucleotides at a single locus in a chromatogram (6-10, 12, 24, 25).
In a previous study, investigators compared the Sequenase fluorescent dye-labeled dideoxynucleotide system (dye-terminator system) with the dye-labeled primer system. The dye-labeled primer system is more appropriate for the quantitative detection of mixtures of HIV-1 drug resistance mutations due to a better correlation of peak heights with the amount of relative DNA concentration (15). We have previously compared the growth fitness of WT virus (pNL4-3) and didanosine-selected mutant L74V by using the Taq fluorescent dye-labeled dideoxynucleotide system. Growth competition assays were performed by infecting human PBM cells with equivalent amounts of infectious virus, and the infected cells were collected at various time points. Sequencing of PCR-amplified genomic DNA and analysis of peak heights in chromatograms reveal that both WT T and mutated G nucleotide peaks were virtually equal on days 4 and 6. However, after day 8, the peak height of G was significantly lower than the WT nucleotide T (24). Furthermore, analysis of replication fitness by utilizing a mathematical model (7, 18) revealed a loss of 11% fitness for L74V viruses (24). This suggested that the dye-labeled terminator sequencing may be useful for quantitative analysis of mixtures containing two viruses. Thus, the utility of a dye-labeled dideoxynucleotide system for quantitative analysis remains controversial and requires further validation.
The dye-labeled dideoxynucleotide system is fast, cost-effective, and available in most sequencing core laboratories. We performed a rigorous analysis to validate the use of dye-terminator system for quantitation of two species of viruses whose genomes differ by a single nucleotide at the same position. Since several antiretroviral drugs, including highly active antiretroviral therapy, select RT mutations L74V and M184V (21, 26) and since these mutations confer a loss of replication fitness to the virus (2, 3, 24, 25), we used these markers for the validation of quantitative assay. Various ratios of WT and mutated DNA fragments were mixed, and sequencing reactions were conducted at the actual melting temperatures of sequencing primer. A strong linear relationship between ratios of peak height versus DNA concentration ratios was demonstrated.
(This work was presented in part at the Third HIV DRP Symposium on Antiviral Drug Resistance of the HIV Drug Resistance Program, National Cancer Institute, at Chantilly, Va., 8 to 11 December 2002.)
|
|
|---|
RT-PCR. Total RNA was isolated from culture supernatants by using the commercial QIAmp Viral RNA mini kit (Qiagen, Valencia, Calif.). A 784-bp RT fragment was amplified by single step RT PCR of RNA with primers RT21 and RT16 (Table 1). Both WT and mutated fragments were purified by using QIAquick PCR purification kit (Qiagen). DNA was quantified by measuring the optical density at wavelengths of 260 and 280 nm in a Beckman spectrophotometer. To avoid variability among independent sequencing reactions due to repeated freeze and thaw, aliquots of purified DNA were stored at -20°C. Each aliquot of DNA was used only once during the sequencing reaction.
|
View this table: [in a new window] |
TABLE 1. Oligonucleotides used in this study
|
Statistical analysis. At least three independent sequencing reactions with five different mixtures of WT and DNA containing L74V or M184V mutations were performed for each of the four primers (3 x 5 variables = 15). The peak heights were measured, and the ratios of WT (T) and mutated (G) nucleotides with forward primer 74F and WT (A) and mutated (C) nucleotides with reverse primer 74R were calculated for codon 74. Similarly, fractions of peak heights for nucleotides A/G and T/C were calculated at RT codon 184 by using the forward primer 184F and the reverse primer 184R, respectively. A regression analysis was performed on peak height ratios and DNA concentration ratios.
|
|
|---|
![]() View larger version (19K): [in a new window] |
FIG. 1. Location of sequencing primers in pol gene of proviral clone pNL4-3. All four forward (F) and reverse (R) primers, 74F, 74R, 184F, and 184R, were selected at equivalent distances with respect to RT codons 74 and 184. Primer sequences are presented in Table 1.
|
The mixtures of WT T (TTA) and mutated G (GTA) nucleotides with forward primer 74F and of WT A (AAT) and mutated C (CAT) nucleotides with reverse primer 74R at RT codon 74 were subjected to PCR, followed by sequencing. Similarly, WT A (ATG) and mutated G (GTG) nucleotides with forward primer 184F and WT T (TAC) and mutated C (CAC) nucleotides with reverse primer 184R at codon 184 were analyzed. Different ratios of the WT and mutated DNA were used in the sequencing reaction to correct for variability. A total of 100 ng of mixtures of DNA was used for the PCRs. The ratios of WT to mutated DNA were 50:50, 60:40, 40:60, 80:20, and 20:80. Both plus and minus DNA strands were sequenced by using forward and reverse primers. The primer locations are presented in Fig. 1 and Table 1.
It is a common practice to use an annealing temperature of 50°C in DNA sequencing core laboratories. However, how this affects the optimal hybridization of primer to target molecules during PCR of sequencing reactions is not clear. We performed all our sequencing reactions at actual melting temperatures (Tm) of specific primers (Table 1). Visual inspection of chromatograms (Fig. 2) revealed that the peak heights of nucleotides were not identical because of differences in base incorporation for nucleotides due to four different fluorescent dye-tagged terminators. However, the ratios of peak heights in relation to DNA concentrations were not affected. For example, we always observed relatively smaller peak heights with primer 74F in comparison to other primers, but the relative ratios of the peak heights of two nucleotides were always in accordance with DNA concentrations (Fig. 2 and Table 2).
![]() View larger version (37K): [in a new window] |
FIG. 2. Quantitation of DNA concentrations in relation to peak heights. Sequencing reactions were performed with mixtures of purified PCR fragments of WT DNA and the DNA containing the L74V and M184V mutations. To analyze nucleotide peaks at RT codon 74, both coding and noncoding DNA strands were sequenced with the forward and reverse primers 74F and 74R, respectively. The peak heights of nucleotides T and G (TTA/GTA) obtained with primer 74F and nucleotides A and C with primer 74R were measured, and the percentage of each nucleotide in the mixture was calculated. Similarly, peak heights of A and G with primer 184F and T and C with primer 184R were determined at RT codon 184. N, equivalent mixture of two nucleotides resulting in equal peak heights (preset values in software of ABI 377). This chromatogram is representative of at least three independent sequencing reactions. The average values for three sequencing runs are presented in Table 2.
|
|
View this table: [in a new window] |
TABLE 2. Fractions of peak heights in relation to DNA concentration
|
Correlation between ratios of DNA concentration and ratios of peak heights. The data obtained from three independent sequencing runs with each of the four primers is presented in Fig. 3. Since five different ratios of WT and mutated DNA were used, each plot (Fig. 3) shows 15 datum points. The peak heights of individual nucleotides were measured at RT codons 74 and 184, and the fractions (%) of the peak heights were calculated for both nucleotides. A scatter plot was used to depict 15 datum points (peak height ratios), and a regression analysis was performed to establish the relation between theDNA concentration of two variables and the peak heights obtained in the chromatogram (Fig. 2). The regression analysis between ratios of peak heights and ratios of related DNA concentrations showed a remarkable linear relation. This linear correlation was observed among independent experiments and also when the values of three experiments were considered together (Fig. 2 and 3). Regression analysis was also performed on the data obtained from individual sequencing runs. For three sequencing runs, the R2 values were as follows: primer 74F, 0.998, 0.998, and 0.999 (combined R2 = 0.994); primer 74R, 0.999, 0.996, and 0.997 (combined R2 = 0.997); primer 184F, 0.996, 0.998, and 0.999 (combined R2 = 0.976); and primer 184R, 0.999, 0.999, and 0.994 (combined R2 = 0.998).
![]() View larger version (17K): [in a new window] |
FIG. 3. Demonstration of a linear relation between DNA concentration and peak heights. Peak-height ratios (mutant/wild type) were plotted versus the corresponding DNA concentration ratios for the four sequencing primers: 74F (A), 74R (B), 184F (C), and 184R (D). Each plot shows 15 datum points (ratios of peak heights) obtained from three independent sequencing runs. A linear regression analysis was performed for each of the three independent runs (see Results). The combined linear regression correlation coefficients (R2) for primers 74F, 74R, 184F, and 184R were 0.994, 0.997, 0.976, and 0.998, respectively.
|
|
|
|---|
In an earlier study, Larder et al. (15) demonstrated that the quantitative detection of HIV-1 drug resistance mutations is best achieved by the fluorescence-labeled sequencing primers. These authors also indicated that Sequenase dye-terminator chemistry failed to produce equal peaks in the presence of 50% of mixture. In effect, visualizing the Sequenase terminator data in the chromatogram (Fig. 2b in reference 15) and the table in that study would indicate a lack of correlation between the DNA concentration and peak heights at several ratios of WT to mutant nucleotide. Moreover, due to the presence of background peaks in the chromatogram obtained from dye-terminator chemistry, the determination of the nucleotide with a minor population was impossible. Our data argue against those observations, and in 12 independent sequencing reactions using dye-terminator chemistry and four different primers, the peak size was always equal in a 50% mixture of WT and mutant DNA. Indeed, various ratios of WT to mutant were found to have a high degree of correlation with peak heights. In at least three independent sequencing runs with four primers (Fig. 2 and 3) and the same sequencer, the variability between runs was minimal. Unlike previous reports, we performed sequencing reactions at actual annealing temperatures (Tm) of primers, allowing an optimal hybridization of primer to target molecules. Additionally, the location of primer with respect to the nucleotide to be analyzed may play a role in determining the peak heights of nucleotides in a mixture. The evaluation of nucleotide peak heights in the same region of the chromatograms allowed a better comparison between the data obtained with four different primers.
Our values of peak heights obtained at mixed-base positions by using dye-labeled terminator chemistry were in fact similar to those reported previously with dye-labeled primer chemistry (15). However, the dye-labeled primer system is cumbersome and expensive due to the requirement of four dye-labeled primers for the detection of four nucleotides and incorporation of 5' universal-tailed M13 sequences in the PCR products. In comparison, dye-terminator chemistry is a single-tube method and easy to assay in a short period of time.
There has been a considerable improvement in sequencing chemistry during the last decade. Specifically, the discovery of Taq polymerase allowed us to perform sequencing reactions at considerably higher temperatures than Sequenase polymerase. This resulted in increased specificity and eliminated problems due to DNA secondary structure. Additionally, substitute nucleotides ITP and UTP can be used to replace dGTP and dTTP to overcome the problem in sequencing G/C-rich regions. As shown in Fig. 2, most of our chromatograms had very few to no background peaks, allowing us to determine 20% of the minor population in the mixture. This is in contrast to a previous report, in which others were unable to detect reliably a 20 to 25% minor population in the mixture with a primer-labeled sequencing chemistry (12). In summary, dye-labeled terminator chemistry could be used reliably for measuring relative fitness of cloned viruses in a mixture, and the strategy described here could be extended to determine replication efficiencies of HIV-1 variants selected during therapy.
We thank Rathna Giri Polavarapu, Emory University DNA Core Facility, and Alisa Morales of Applied Biosystems for valuable suggestions during this project. We also thank Cheryl Strauss for proofreading the manuscript.
|
|
|---|
Gly) at codon 69. J. Virol. 74:10958-10964.
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»