Previous Article | Next Article ![]()
Journal of Clinical Microbiology, June 2004, p. 2711-2717, Vol. 42, No. 6
0095-1137/04/$08.00+0 DOI: 10.1128/JCM.42.6.2711-2717.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Department of Pathology, The Johns Hopkins Medical Institutions, Baltimore, Maryland,1 AIDS Research and Retrovirus Discovery, Abbott Laboratories, Abbott Park, Illinois,2 Laboratoire de Santé Hygiène Mobile,4 Université de Yaoundé, Yaoundé, Cameroon,5 Celera Diagnostics, Alameda, California3
Received 24 September 2003/ Returned for modification 8 January 2004/ Accepted 1 March 2004
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
HIV-1 strains are genetically diverse. Nine phylogenetically pure subtypes (A, B, C, D, F, G, H, J, and K) and 15 circulating recombinant forms (CRFs) are currently recognized for HIV-1 group M (HIV-1 nomenclature proposal, Los Alamos HIV Sequence Data Base [http://hiv-web.lanl.gov]). HIV-1 drug resistance mutations and their effects on drug susceptibility have been studied predominantly in patients infected with subtype B strains. Although HIV-1 subtype B causes the majority of infections in the United States and Europe, other subtypes and CRFs cause the overwhelming majority of HIV-1 infections worldwide. Analysis of drug resistance in non-subtype B HIV-1 is becoming increasingly important given the expanding availability of antiretroviral drugs throughout the world and the increasing prevalence of non-subtype B HIV-1 in regions where antiretroviral drugs are widely used (19). HIV-1 subtype can influence viral susceptibility to antiretroviral drugs (2, 6), as well as the pattern of resistance mutations that emerge after drug exposure (4, 15). HIV-1 subtype may also influence the rate at which specific resistance mutations emerge after drug exposure (7, 8, 20) and the amino acids selected at specific positions under drug pressure (14, 21). Novel, subtype-specific drug resistance mutations may also emerge in non-subtype B strains at positions not associated with drug resistance in subtype B (1). Further studies are needed to examine the emergence and fading of drug resistance mutations in individuals with non-subtype B HIV-1 and to examine genotypic correlates of drug resistance in different HIV-1 subtypes.
The ViroSeq HIV-1 Genotyping System (Celera Diagnostics, Alameda, Calif.) is an integrated system for genotyping HIV-1 protease and RT that is cleared by the FDA for clinical use. This system utilizes 9 or 10 DNA primers for analysis: 1 for reverse transcription, 2 for PCR amplification, and 6 or 7 for sequence analysis. The performance of the ViroSeq system was intensively evaluated and validated for the genotyping of HIV-1 subtype B. Subtype-specific sequence differences at the sites of primer binding may complicate the genotyping of non-subtype B HIV-1. In a previous report, a high level of performance of the ViroSeq system for analysis of 130 plasma samples with HIV-1 non-subtype B strains from Uganda (mostly subtypes A and D) was demonstrated (16). Two recent reports further evaluated the performance of the ViroSeq system for analysis of non-subtype B HIV-1 (10, 12). However, those studies used predominantly cultured isolates rather than plasma and included relatively small sample sets. In this report, we evaluated the performance of the ViroSeq system for analysis of genetically diverse strains of HIV-1 by using plasma samples obtained from 126 HIV-1-infected individuals from Brazil, Cameroon, Uganda, Thailand, Argentina, and South Africa. This analysis focused on the ability of the ViroSeq system to successfully amplify DNA for sequencing and provide complete DNA sequences of the regions of interest. The analysis included an interlaboratory comparison of assay performance.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Viral load determination. The LCx HIV RNA Quantitative Assay (Abbott Laboratories, Abbott Park, Ill.) (not licensed for use in the United States) was performed according to the manufacturer's specifications by using a sample volume of 1.0 ml. This competitive RT PCR assay targets the pol integrase region of HIV-1 (13). The upper and lower limits of quantification for the 1.0-ml sample preparation volume are 1,000,000 (6.0 log10) and 50 (1.7 log10) HIV RNA copies/ml, respectively.
HIV-1 genotyping. HIV-1 genotyping was performed by using the Celera Diagnostics ViroSeq HIV-1 Genotyping System (version 2.0) according to the manufacturer's instructions with 0.5 ml of plasma. Analysis with this system begins with HIV-1 RNA isolation and involves high-speed centrifugation of plasma to pellet virus particles, lysis of virus particles with a chaotropic agent to release viral RNA, and isopropanol-ethanol precipitation for RNA recovery. Ten microliters of the resuspended RNA was used for the reverse transcription (65°C for 30 s, 42°C for 65 min, 99°C for 5 min) with Moloney murine leukemia virus RT followed by a single 40-cycle PCR (50°C for 10 min, 93°C for 12 min, 93°C for 20 s, 64°C for 45 s, 66°C for 3 min, 72°C for 10 min) with AmpliTaq Gold DNA polymerase (Applied Biosystems, Foster City, Calif.). The PCR yielded a 1.8-kb DNA product. PCRs were performed with the Amperase uracil N-glycosylase contamination control system to reduce the risk of contaminating PCRs with products from previous amplification reactions. PCR products were purified by using spin columns and analyzed by agarose gel electrophoresis prior to sequencing. DNA sequencing was performed by using premixed BigDye terminator (Applied Biosystems) sequencing reagents with seven different primers. Conditions for cycle sequencing (25 cycles) were 96°C for 10 s, 50°C for 5 s, and 60°C for 4 min. BigDye terminator chemistry provides 99.91% accuracy at 580 bases for the ABI PRISM 3100 Genetic Analyzer (N. Marlowe, unpublished observations). Sequence data were analyzed by using Celera Diagnostics ViroSeq HIV-1 Genotyping System software (version 2.5), which assembles sequence data from the primers into a contiguous sequence that can be inspected for the identification of drug resistance mutations. Processed sequences include the region coding for protease amino acids 1 to 99 and RT amino acids 1 to 335 (full-length sequences). Corresponding genomic positions within the HXB2 reference isolate are nucleotides 2253 to 2549 for protease and 2550 to 3554 for RT. Bidirectional sequence data were obtained for protease and for RT amino acids 1 to 321. Operators at both sites (Johns Hopkins University and Celera Diagnostics) were trained and certified by the manufacturer and were proficient in running the ViroSeq assay. Each laboratory was blinded to the results obtained by the other laboratory for assay comparison.
Phylogenetic analysis. HIV-1 sequences trimmed to 1,260 nucleotides (nt), corresponding to protease amino acids 1 to 99 and RT amino acids 1 to 321 (297 and 963 nt, respectively), were used to determine subtypes. Nucleotide alignments were obtained with the biosequence editor Seqapp v1.9a169 (obtained from http://iubio.bio.indiana.edu/soft/molbio/seqapp/). Alignments included 47 reference sequences from the Los Alamos National Laboratory (http://hiv-web.lanl.gov) recommended for HIV-1 subtyping, as well as 4 additional reference sequences, for a total of 51 sequences: A1-U455 (M62320), A1-SE7253 (AF069670), A1-92UG037 (U51190), A1-Q23-17 (AF004885), A2-94CY017-41 (AF286237), A2-97CDKTB48 (AF286238), B-HXB2 (K03455), B-JRFL (U63632), B-WEAU160 (U21135), B-RF (M17451), C-ETH2220 (U46016), C-96BW0502 (AF110967), C-95IN21068 (AF067155), C-92BR025 (U52953), D-ELI (K03454), D-NDK (M27323), D-84ZR085 (U88822), D-94UG114 (U88824), F1-VI850 (AF077336), F1-FIN9363 (AF075703), F1-93BR020 (AF005494), F2-MP255 (AJ249236), F2-MP257 (AJ249237), G-HH8793-12.1 (AF061641), G-SE6165 (AF061642), H-VI997 (AF190128), H-90CF056 (AF005496), H-VI991 (AF190127), J-SE91733 (AF082395), J-SE9280 (AF082394), K-EQTB11C (AJ249235), K-MP535 (AJ249239), CRF01_AE-93TH253 (U51189), CRF01_AE-CM240 (U54771), CRF01_AE-90CF402 (U51188), CRF02_AG-DJ263 (AF063223), CRF02_AG-DJ264 (AF063224), CRF02_AG-IBNG (L39106), CRF03_AB-KAL153-2 (AF193276), CRF03_AB-RU98001 (AF193277), CRF04_cpx-97PVCH (AF119820), CRF04_cpx-97PVMY (AF119819), CRF04_cpx-CY032 (AF049337), CRF05_DF-VI1310 (AF193253), CRF05_DF-VI961 (AF076998), CRF12_BF-A32989 (AF408630), CRF12_BF-ARMA159 (AF385936), N-YBF30 (AJ006022), SIV-CPZGAB (X52154), O-ANT70 (L20587), and O-MVP5180 (L20571). Alignment and gap stripping were performed manually. Distances between the sequences were calculated with DNADist by using the Kimura two-parameter model as an optimal substitution model with a transition-to-transversion ratio of 1.5 (PHYLIP, version 3.572, obtained from http://evolution.genetics.washington.edu/phylip.html). Neighbor-joining and consense were used to create phylogenetic trees with 500 bootstrap replications (SeqBoot). Consensus trees were displayed with TreeView. Bootstrap values of >80 were considered acceptable for subtype assignment.
Sequencing primer mismatch analysis. Nucleotide sequences generated by the ViroSeq HIV-1 genotyping system software (version 2.5) were aligned with ViroSeq sequencing primers B, C, F, G, and H by the Clustal W method (MegAlign module, Lasergene version 5.0.1; DNASTAR Inc., Madison, Wis.) to identify nucleotide mismatches between oligonucleotide primer and target sequences. Since primers A and D are not encompassed within the processed sequence, raw sequence files generated from primer F were aligned with sequencing primers A and D.
Nucleotide sequence accession number. GenBank accession numbers for nucleotide sequences of protease and RT are AY444180 to AY444303.
| RESULTS |
|---|
|
|
|---|
|
ViroSeq assay performance: sequencing. The locations and positions of sequencing primers in the ViroSeq system are shown in Fig. 1. Primers A, B, and C sequence the sense strand of the PCR product, and primers F, G, and H sequence the antisense strand of the PCR product. An analysis of the performance of individual sequencing primers was performed at Johns Hopkins University (Table 2). Full-length sequences were obtained for 124 (99%) of the 125 samples sequenced. For 95 (76%) of those samples, bidirectional sequencing data were obtained along the entire region analyzed. Sequencing was successful on the first attempt for 671 (89%) of 750 reactions with these six primers and was successful for 41 (52%) of the 79 reactions that were repeated after an initial primer failure. Sequencing was successful for 108 (86%) of 125 primer A reactions. The success rate for primer A was the lowest success rate among the six primers. Since primer A binds to the highly heterogeneous gag region of the HIV-1 genome, primer D may be used as an alternate sense primer if sequencing is unsuccessful with primer A. Sequencing with primer D was successful for 6 (30%) of 20 samples that failed to sequence with primer A; for these samples, only unidirectional sequence data were obtained for the region encoding the first 17 to 18 amino acids of protease. Sequencing was successful with either of the two alternate primers (primer A or primer D) for 114 (91%) of the samples tested and for 100, 99, 92, 98, and 94% of reactions with primers B, C, F, G, and H, respectively (Table 2). One sample failed sequencing with primers A, D, F, and G and was considered a sequencing failure. This sample also failed sequencing at Celera Diagnostics. To confirm the absence of sample cross-contamination, sequences from each sample were trimmed to 1,260 nt, aligned, and compared phylogenetically. No evidence of sample cross-contamination was observed for either site.
|
|
Performance of sequencing primers for different subtypes and intersubtype recombinant forms. The performance of the ViroSeq sequencing primers at Johns Hopkins University was analyzed for samples of each subtype or CRF (Table 2). The majority of samples of all subtype categories yielded full bidirectional sequence data with no primer failures. All six primers (primers A or D, B, C, F, G, and H) were successful for over 90% of samples in all subtype categories with the following exceptions: two primer A/D failures and two primer F failures among 16 subtype A/A2 samples, two primer G failures among 13 subtype C samples, three primer A/D failures among 9 CRF01_AE samples, one primer F and one primer H failure among 9 subtype F/F2 samples, four primer H failures among 32 CRF02_AG samples, and two primer F failures among 14 intersubtype recombinant samples.
To examine nucleotide conservation at the sequencing primer binding sites, ViroSeq-generated sequences were aligned with the oligonucleotide primer sequences. The mean number and range of mismatches for the primary sequencing primers (primers A, B, C, F, G, and H) were tabulated for each subtype and CRF (Table 3). Overall, primer sites were well conserved. In most cases, the mean number of mismatches was 1 or less. The mean and range of mismatches were slightly higher for primer A with subtype F and CRF02_AG, primer B with CRF01_AE, and primer G with subtype C. The least overall conservation was observed for primer H. This sequencing primer had a mean of one mismatch or fewer with subtypes A, B, D, and H but had means of 1.1 to 2.6 mismatches with subtypes C, F, and G and with CRF01_AE, CRF02_AG, and other intersubtype recombinant strains. For primer H, the mean number of mismatches was highest (2.6) for CRF02_AG strains, ranging from 0 to 5.
|
The nucleotide sequence identity levels of the 124 sequence pairs ranged from 98 to 100%, with a median identity of 99.8%. Identity was >99% for 118 (95%) of the 124 paired sequences. The nucleotides identified (base calls) differed at 488 (0.3%) of the 156,364 nucleotide positions analyzed. Of these base-calling discrepancies, 479 represented the identification of a nucleotide mixture in one laboratory and the identification of one of the two corresponding nucleotides in the other laboratory (e.g., A+G = R in one laboratory and A in the other laboratory). Five of the 460 different base calls influenced the amino acids detected at positions associated with antiretroviral drug resistance: L10L/V (two base calls), L63L/P (two base calls), and Y181Y/C (one base call). Review of electropherogram data from both laboratories for these positions revealed the following: (i) both of the L10L/V mixtures were clearly detected in data from one laboratory but were not detected in data from the other laboratory, (ii) both of the L63L/P mixtures were clearly detected in data from one laboratory and were suggested by data from the second laboratory (background signal interfered with the identification of the mixtures), and (iii) the Y181C substitution was detected in one laboratory in a region where sequences from three primers overlapped. The mixture was seen in only one direction but met the criteria for mixture identification described in the ViroSeq system software manual; it was not detected in the second laboratory. Clear differences in the sequences (e.g., A in one laboratory and G in the other) were detected at only 9 (<0.01%) of the 156,364 nucleotide positions analyzed. None of these differences were at positions associated with antiretroviral drug resistance. A retrospective review of electropherograms revealed that all nine of the differences reflected the misidentification of bases at one laboratory. In six cases, misidentification may have occurred due to mobility shifts that were missed by the operator during manual inspection of the automated calls made by ViroSeq software; four of these nucleotides were in a single sequence encoding three amino acids in RT. In two cases, there was high background signal in the region.
Detection of amino acid changes at positions associated with antiretroviral drug resistance. With the ViroSeq software, a contiguous consensus sequence derived for each sample is compared to a reference sequence, HXB-2 (GenBank accession number K03455), and differences from the reference sequence are tabulated. Non-subtype B sequences often differ from subtype B reference sequences at several positions, reflecting naturally occurring amino acid polymorphisms. The following amino acid changes were detected (the percentage of sequences with each amino acid change is shown in parentheses): in protease, M36I (90%), L63P (16%), K20R (14%), V77I (10%), L10V (9%), L10I (4%), A71T (2.4%), and M46I (one sample), and in RT, S68G (1.6%) and V118I (1.6%), plus the substitutions E44D, D67N, K70R, Y181C, T215F, K219Q, and G333E, each detected in a single sample. Four of these substitutions (D67N, K70R, T215F, and K219Q) are associated with resistance to zidovudine and were detected in a single individual from Brazil with subtype B HIV-1. The Y181C substitution is associated with high-level resistance to nevirapine; this amino acid substitution was detected in a single Cameroonian sample infected with a CRF02_AG strain.
| DISCUSSION |
|---|
|
|
|---|
The most consistent difference in performance observed in the present study compared to previously reported results for subtype B strains (5) was in the performance of the A and D alternate primers, which bind to the highly heterogeneous gag region. For subtype B strains, sequencing was successful with the A and D primers for 92 and 88% of reactions, respectively, and at least one of the primers (either primer A or primer D) was successful with 98.4% of the samples. In contrast, for the genetically diverse strains in this study, the A and D primers were successful for 86 and 30% of the reactions, respectively. Both primers (A and D) failed for 9% of the genetically diverse samples, resulting in some regions of unidirectional sequence data.
The genetic heterogeneity of HIV-1 viruses presents a challenge for molecular assays that depend on oligonucleotide primer hybridization. Nucleotide mismatches between the primer(s) and template have the potential to reduce the efficiency of hybridization, leading to inefficient amplification or failed sequencing runs. The genetically diverse panel of samples used in the present study provided an opportunity to examine the level of nucleotide conservation within primer sites used for sequence analysis across multiple subtypes and CRFs. Consistent with the performance (>95% success rate) of the primary set of ViroSeq sequencing primers (primers A to C and F to H), a relatively high degree of nucleotide conservation was observed at the primer binding sites (Table 3). One or no mismatches were present in primer A, B, C, F, and G sites in 90% or more of the samples. For primer A, B, C, F, and G sites, two or fewer mismatches were present in 97, 99, 100, 99, and 99% of the panel members, respectively. The primer H site was less conserved, with two or fewer mismatches for 77% of this genetically diverse sample panel. Interestingly, although the level of conservation was somewhat lower in primer H, it accounted for only 20% of the failed sequencing reactions in our study. Primers A and F accounted for, respectively, 46 and 26% of the failed sequencing reactions, although their overall nucleotide conservation was higher than that for the primer H binding site. The overall level of sequence conservation within the target binding site of a primer is just one of many factors that can influence the efficiency of primer hybridization. The position and nature of the mismatch, primer length, and stringency (temperature and buffer conditions) are additional factors that contribute to mismatch tolerance. For primer F, all eight sequencing failures (reproducible at both laboratories) could be attributed to a mismatch at the 3' terminus of the primer. The samples with primer failures included a subset (6 to 11%) of subtype A, C, D, and F and CRF02_AG specimens. In all cases, primer A provided sequence coverage on the opposite strand. Two primer G sequencing failures, both due to 3'-terminal mismatches on subtype C strains, were observed. Thus, the position of the mismatch was critical for all sequencing failures with primers F and G. Sequencing failures with primers A, C, and H were generally associated with internal mismatches near the 3' ends of the primers or due to a higher total numbers of mismatches (three to five mismatches for primer H).
Overall, the sequencing success rates for the primers were highly concordant between laboratories. However, in five of six cases (primer F, one subtype A specimen; primer A, one subtype A, one subtype D, and two CRF01_AE specimens) in which sequencing failures were observed at Johns Hopkins University (Table 2), no mismatches were present in the primer binding sites. In the sixth case, only a single mismatch, near the 5' end of primer A, was identified. Thus, no obvious molecular basis for failure was evident. For all six specimens, sequences were successfully generated in the second laboratory, suggesting that the failures may have resulted from technical issues.
Primer D is included in the ViroSeq kit as an accessory primer in the event of sequencing failure with primer A. Consistent with the relatively low success rate (30%) for non-subtype B strains, the primer D binding site is the least conserved of the seven ViroSeq sequencing primer binding sites. Nearly all cases of failure with primer D can be attributed to a relatively high number of total mismatches. Our data are consistent with a previous report of a low success rate for primer D with CRF01_AE strains (17). Nevertheless, primer D was successful in five cases in which primer A failed (one subtype B specimen, one subtype D specimen, two CRF02_AG specimens, and one CRF01_AE specimen) and thus is of value, even with non-subtype B strains.
The performance of the ViroSeq system (version 2.0) with non-subtype B HIV-1 was evaluated in three previous reports. Those reports also examined the performance of another commercial genotyping system, the TRUGENE HIV-1 genotyping kit (Bayer Diagnostics, Tarrytown, N.Y.), which is also cleared by the FDA for clinical use. In the first report (10), the performance of the ViroSeq system was evaluated with 15 cultured isolates of non-subtype B HIV-1 (3 subtype A, 2 subtype C, 3 subtype D, 2 subtype F, 1 subtype G, 1 subtype J, and 3 recombinant [1 A/D and 2 CRF02_A/G] isolates) and three non-subtype B plasma samples (one each of subtypes C, D, and H). All of the samples were successfully amplified and sequenced with the ViroSeq system (version 2.0). Full bidirectional sequence data were obtained for 11 (73%) of the 18 non-subtype B samples. Primer failures included both the A and D primers for three samples (subtypes F, J, and H), the F and G primers for one subtype A sample, the F primer for one subtype C sample, and the G primer for one subtype D and one CRF02_AG sample. Genotypes were also obtained for all of the non-B samples using the TRUGENE system and supplemental primers provided by the manufacturer. Samples analyzed with the TRUGENE system were noted to have more frequent sequence ambiguities.
In the second report (12), 34 cultured isolates, 27 of which were non-subtype B isolates (2 subtype A, 2 subtype A/G, 6 subtype C, 2 subtype D, 8 CRF01_AE, 4 subtype F, 2 subtype G, and 1 subtype H), were analyzed. Genotyping was successful for all of the non-subtype B samples with the ViroSeq system. Sequencing with primer A was successful for 24 (89%) of the non-subtype B samples, while sequencing with the alternate primer D was successful for only 8 (30%) of the samples. Both primers A and D failed for three isolates (of subtypes C, E, and H) and one sample for which primer H failed (a subtype C sample). Analysis with the TRUGENE system required the use of the version 1.5 primers for three (11%) of the 27 non-B samples (1 subtype C and 2 subtype G samples). Sequencing primer failures, primarily in the protease region, were described for 22 (81%) of the 27 non-subtype B samples. Samples with subtypes B and D performed better with the PR set of protease sequencing primers, whereas the P2 protease primer set was preferable for the other subtypes. The P2 primer set was successful for all samples tested but provides a shorter sequence in the protease region. Both sets of protease primers are provided with the FDA-cleared version of the TRUGENE system.
In the third report (3), genotyping was performed by using viral isolates from cell culture. Genotyping with the version 2.0 ViroSeq system was successful for 20 (80%) of 25 non-subtype B isolates. In contrast, genotyping with the TRUGENE system (version 1, without the use of supplemental primers) was successful for only 12 (50%) of 24 non-subtype B isolates. The five non-subtype B isolates that failed with the ViroSeq system included three subtype A isolates, one subtype C isolate, and one subtype J isolate. In all five cases, all seven sequencing primers failed. This type of result was not observed in the present study or in the other two studies cited. It seems unlikely that all five samples diverged genetically at the binding sites of all seven sequencing primers. Instead, these results suggest that a contaminant (e.g., from cell culture) might have interfered with the sequencing reactions. In contrast to the three previous studies described, the present study used only plasma samples for analysis.
The present report confirms the value of the ViroSeq system for the analysis of plasma samples from countries harboring diverse HIV-1 subtypes. In the present study, a high level of interlaboratory assay reproducibility was observed with this system. Further studies are needed to extend this analysis to include non-subtype B samples with low viral loads and non-subtype B samples from patients failing antiretroviral therapy. The FDA-cleared ViroSeq system offers an advantage for drug resistance monitoring in that a single set of PCR and sequencing primers can be used for the analysis of genetically diverse strains of HIV-1. The inclusion of the uracil N-glycosylase contamination control system reduces the likelihood of PCR carryover. The ViroSeq system provides a single, uninterrupted pol region sequence, and sequence data are provided in the standard FASTA format, simplifying phylogenetic analysis for subtype determination and other studies. Moreover, the ViroSeq system has been FDA-cleared for use on multiple platforms, including the high-throughput ABI Prism 3100 and 3700 DNA analyzers, making it an attractive tool for patient monitoring.
| ACKNOWLEDGMENTS |
|---|
This work was supported by the HIV Network for Prevention Trials (HIVNET), sponsored by the U.S. National Institute of Allergy and Infectious Diseases (NIAID), National Institutes of Health (NIH), DHHS, through contract N01-AI-35173 with Family Health International, contract N01-AI-45200 with Fred Hutchinson Cancer Research Center, and subcontract NOI-AI-35173-417 with Makerere Univ.; by the HIV Prevention Trials Network (HPTN), sponsored by the NIAID, National Institutes of Child Health and Human Development (NICH/HD), National Institute on Drug Abuse, National Institute of Mental Health, and Office of AIDS Research of the NIH, DHHS (contracts U01-AI-46745 and U01-AI-48054); by the Adult AIDS Clinical Trials Groups (NIH, Division of AIDS, NIAID); and by grant R01-HD-42965-01.
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Antimicrob. Agents Chemother. | Clin. Microbiol. Rev. |
|---|---|
| Clin. Vaccine Immunol. | ALL ASM JOURNALS |
|---|