| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, December 2007, p. 4267-4275, Vol. 51, No. 12
0066-4804/07/$08.00+0 doi:10.1128/AAC.00962-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Department of Microbiology and Immunology, F. Edward Hébert School of Medicine, Uniformed Services University of the Health Sciences, 4301 Jones Bridge Road, Bethesda, Maryland 20814-4799
Received 25 July 2007/ Returned for modification 28 August 2007/ Accepted 25 September 2007
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The obligate intracellular nature of these pathogens and their unique biphasic lifestyle pose challenges for the treatment of associated infections (53). Indeed, an efficient antichlamydial drug not only must achieve adequate intracellular penetration and concentration but also must be able to inhibit the metabolically active form of the organism (RBs) in its ability to undergo DNA or protein synthesis, cell division, or differentiation into infectious particles (EBs). Currently, the recommended first-line therapeutic regimens for chlamydial infections are the tetracyclines and the macrolide azithromycin (AZM), which inhibit bacterial translation by binding to the 30S and 50S ribosomal subunits, respectively (53). In addition to its chemotherapeutic use, AZM is used for chemoprophylaxis of blinding trachoma, a chronic infection caused by ocular serovars of C. trachomatis, common in developing countries (15).
Although chlamydial infections are characterized by a high recurrence rate despite appropriate drug therapy (23), the majority of clinical failures has been attributed to reinfection or relapse following deviation of the organism to persistent, nonreplicating chlamydial forms that are phenotypically antibiotic resistant but can revert to typical RBs at the end of the treatment (22). On the other hand, clinical failures linked to true genotypic resistance due to chromosomal mutations have rarely been reported, suggesting that mutations which confer antibiotic resistance in chlamydiae are not selected for in vivo (23, 61).
Resistance to ribosome-targeting antibiotics due to mutations in rRNA genes has been observed mainly in pathogens possessing a low copy number of rRNA operons, because the selective advantage of a mutation in one rRNA copy is usually masked by the abundance of wild-type drug-sensitive rRNAs transcribed from the other unmutated gene copy (or copies) (49). Interestingly, the Chlamydiaceae harbor either one or two rRNA operons depending on the species, implying that resistance to the current antichlamydial drugs could emerge from ribosomal mutations. Previously, we showed that spontaneous resistant mutants of C. psittaci 6BC with distinct mutations in the unique 16S rRNA gene could be isolated in the plaque assay in the presence of the aminoglycoside spectinomycin (at a frequency of 5 x 10–5) but not in the presence of tetracycline, another 16S-targeting antibiotic (frequency, <3 x 10–9) (6, 7). Because the physiological burden of resistance mutations is an important factor that affects the appearance, stability, and maintenance of the phenotype (3, 27), we suggested that chromosomal point mutations conferring resistance to tetracycline incurred a high fitness cost for the bacteria that was reflected in the rarity of genotypic antibiotic resistance among chlamydial clinical isolates. Indeed, stable tetracycline resistance in chlamydiae has been reported only for Chlamydia suis, and this resistance probably arose from horizontal gene transfer rather than spontaneous genetic mutations affecting the bacterial ribosome (13).
In this study, we examined the contribution of spontaneous chromosomal mutation to emergence of AZM resistance in chlamydiae in vitro. AZM-resistant chlamydial isolates were isolated either in the plaque assay using a clonal population of C. psittaci 6BC (one rRNA operon) or following serial passages in AZM for C. trachomatis L2, which harbors two rRNA operons. We found that mutations in the unique C. psittaci 23S rRNA gene at position 2058 or 2059 (Escherichia coli numbering system) and in C. trachomatis ribosomal protein L4 lowered bacterial sensitivities to multiple antibiotics. However, these mutations were associated with a high physiological burden, as evidenced by a reduced production of infectious particles in tissue culture in the absence of antibiotic. The implications of these results for the use of macrolides in treatment of chlamydial infections are discussed.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Titration and antimicrobial susceptibility testing of C. trachomatis and C. psittaci. The susceptibilities of the chlamydial strains to AZM, erythromycin (ERM), josamycin (JOS), spiramycin (SPI), clindamycin (CLI), virginiamycin M1 (VIR), and chloramphenicol (CLM) were determined in a plaque assay. The MIC was defined as the drug concentration that inhibited the development of 105 chlamydial PFU in a confluent L2 monolayer in a 60-mm dish (7). All antibiotics were purchased from Sigma Chemical Co.
Isolation of chlamydial mutants. To isolate spontaneous AZM-resistant variants, 60-mm dishes were infected with 107 to 108 PFU, corresponding to multiplicities of infection (MOI) of 1 and 10, respectively, and AZM was added 2 h postinoculation (p.i.) at concentrations ranging from 200 ng/ml to 1 µg/ml. The frequency of spontaneous mutation to drug resistance was determined by dividing the number of PFU on selective medium by the number of PFU added to the monolayer (as measured by titration of PFU in the absence of antibiotic) (7).
AZM-resistant populations of C. trachomatis L2 were isolated following enrichment by serial passage in subinhibitory concentrations of AZM. Starting with 2 x 109 wild-type C. trachomatis L2 infectious particles, we ended up with a population that grew in the presence of 200 ng/ml of AZM after 30 passages.
To monitor the stability of the resistance phenotype, we compared the number of PFU obtained in the absence or presence of antibiotic following the growth of each variant for a minimum of 14 days in the plaque assay in the absence of AZM.
PCR and DNA sequencing of the macrolide resistance targets. PCR amplification and DNA sequencing were used to determine whether resistance to AZM was due to a mutation in the 23S rRNA gene or in rplD and rplV (encoding L4 protein and ribosomal L22 protein, respectively). PCR primers are listed in Table 1. Total genomic DNA was prepared from infected cells with DNeasy tissue kits (Qiagen). Alternatively, single plaques were collected in 10 µl sterile H2O and diluted with 230 µl SPG (250 mM sucrose, 10 mM sodium phosphate, 5 mM L-glutamic acid). PCR amplification was then performed on 10 µl of the plaque lysate which had been denatured for 10 min at 96°C beforehand, in a 50-µl reaction mixture using Platinum Taq high-fidelity DNA polymerase (Invitrogen) supplemented with 0.8 µg/ml bovine serum albumin. After 30 cycles, the PCR product was ligated in pGEM-T (Promega) and sequenced with primers annealing to the vector (i.e., PUC-R and PUC-F) or to the insert (Table 1).
|
DNA sequences for each antibiotic resistant isolate were aligned using Clone Manager 8 (Scientific & Educational Software, Durham, NC) and compared to the respective DNA sequence obtained for the wild-type parental strain.
Physiological cost associated with the mutations. (i) Pairwise competition experiment.
C. psittaci 6BC wild-type BCRB and one isogenic representative of each AZM-resistant variant were coinfected at a ratio of
1:1 to a MOI of 1 in confluent mouse fibroblast monolayers in 60-mm dishes and incubated at 37°C in 5% CO2. After 2 h of infection, the inoculum was removed. The cells were washed twice with Dulbecco's modified Eagle medium (DMEM; GIBCO) and incubated in DMEM supplemented with 10% fetal bovine serum, 1x MEM nonessential amino acids (Sigma-Aldrich), and 2 µg of cycloheximide per ml. EBs were harvested in triplicate after sonication of the infected cells at 19, 24, 29, 34, and 46 h p.i. and stored at –80°C in 400 µl SPG. Approximately 8 x 106 infectious particles from the mixed infection obtained at 46 h p.i. were passed a second time into fresh monolayers in 60-mm dishes, allowed to grow for another 46 h, and harvested as before. Titers of serial dilutions of each harvest were determined in duplicate in the plaque assay in both drug-free DMEM (total PFU) and drug-containing DMEM (AZM-resistant PFU). Plaques were counted after 10 days of incubation. The plating efficiency of each resistant mutant was unaffected by the presence of AZM. The twofold PFU increase rate was estimated by using Prism 3.0 software from a plot of ln(PFU) = f(time), where the slope is ln2/twofold PFU increase rate (hours). The competition index was defined as the ratio of the output mutant/wild-type ratio to the input mutant/wild-type ratio (6). Additionally, the sizes of a minimum of 60 plaques formed by each C. psittaci 6BC strain in the absence of selection at 10 days p.i. were determined and averaged.
(ii) Pure culture. C. trachomatis L2 infectivity was determined while the strains were growing in pure culture in the absence of selection. Confluent mouse fibroblast monolayers in 60-mm dishes were infected with each strain to an MOI of 0.1 and incubated at 37°C in 5% CO2. After 2 h of infection, the inoculum was removed and the cells were incubated in DMEM supplemented with 10% fetal bovine serum, 1x MEM nonessential amino acids, and 2 µg of cycloheximide per ml. EBs were harvested in triplicate after sonication of the infected cells at 46 h p.i. and stored at –80°C in 400 µl SPG. Titers of serial dilutions of each harvest were determined in duplicate in the plaque assay in the absence of selection. Plaques were counted after 14 days of incubation. Sizes of 173 C. trachomatis L2 wild-type plaques and 287 and 270 plaques for L2AZM#23 and L2AZM#31, respectively, were determined and averaged. The EB generation rate at 46 h p.i. was defined as the number of total PFU obtained at 46 h p.i. divided by the number of PFU used to infect initially.
Nucleotide sequence accession numbers. C. psittaci 6BC rplV and rplD sequences determined in the present study have been deposited in GenBank under accession numbers EU035809 and EU035810, respectively.
| RESULTS |
|---|
|
|
|---|
In our laboratory, sensitivity of Chlamydia spp. to antibiotics is measured in the plaque assay, and the MIC is defined as the concentration of drug that inhibits the development of 105 PFU. The number of input bacteria corresponds to an MOI of 0.01 in a confluent monolayer of L2 mouse fibroblasts in 60-mm dishes (7). The AZM MIC is 100 ng/ml for C. psittaci 6BC. When monolayers were infected with a minimum of 8 x 107 PFU in the presence of AZM ranging from 200 to 1,000 ng/ml, resistant plaques appeared at a frequency of 1.35 x 10–8 ± 0.15 x 10–8 for the laboratory stock of C. psittaci 6BC (population genetically heterogeneous) and 0.75 x 10–8 ± 0.5 x 10–8 for the clonal population BCRB. Fifteen independent AZM-resistant plaques were isolated from both the heterogeneous and the clonal populations of C. psittaci 6BC, expanded in the presence of AZM, and further analyzed.
Biochemical studies have shown that AZM reversibly binds to the large ribosomal subunit in the vicinity of the peptidyl transferase center and causes cell growth arrest due to inhibition of protein synthesis (11, 44, 52, 59). More precisely, AZM interacts with bacterial 23S rRNA by connecting hairpin 35 in domain II of the rRNA and the peptidyl transferase loop in domain V. The sequence of the AZM binding site in the 23S rRNA gene was determined by sequencing approximately 900 nucleotides of the 1,400-bp PCR fragment amplified from the parent strain C. psittaci 6BC and 30 independent AZM-resistant mutants (Table 2). Each mutant showed a single mutation at A2058 or A2059 (E. coli numbering system), both of which are known to confer the highest levels of macrolide resistance in other organisms (60). One AZM-resistant representative of each mutant class, i.e., BCRBAZ1, BCRBAZ2, and BCRBAZ5, was chosen from the C. psittaci 6BC clonal population and expanded for two developmental cycles to perform further phenotypic and physiological characterizations on bacterial populations which were as isogenic as possible (Table 2). Because mutations in ribosomal proteins L4 (rplD) and L22 (rplV) have also been reported to confer macrolide resistance in some organisms (12), we amplified, cloned, and sequenced both genes in the isogenic mutant strains. We did not find any changes in the DNA sequences of these two genes. This result indicates that spontaneous AZM resistance arose in C. psittaci 6BC from single mutations in the unique 23S rRNA gene, resulting in either an A2058C, A2059C, or A2059G mutation.
|
|
The stability of the acquired resistance was then tested by growth in the absence of selective pressure (i.e., no AZM). Plaques formed by each C. psittaci 6BC variant grown for 14 days in the absence of selection displayed the same number of infectious particles when titers were determined in the presence or absence of 200 ng/ml of AZM, indicating that the resistance phenotype of these three chlamydial mutants was stable (data not shown).
Ribosomal mutations observed in AZM-tolerant isolates of C. trachomatis L2. Although C. trachomatis L2 has the same sensitivity to AZM as C. psittaci 6BC (MIC of 100 ng/ml in the plaque assay), we were unable to obtain spontaneous AZM-resistant isolates of C. trachomatis L2 (frequency, <1.3 x 10–8 on 200 and 400 ng/ml AZM and <4.0 x 10–9 on 800 ng/ml and 1.5 µg/ml AZM). Because C. trachomatis L2 harbors two copies of the 23S rRNA gene, we hypothesized that a single mutation in only one ribosomal operon might not confer a level of AZM resistance high enough to be selected in one step in the plaque assay. Therefore, we reasoned that serial passage in subinhibitory concentrations of antibiotic would allow the enrichment of putative bacteria with low-level AZM resistance. Starting with 2 x 109 wild-type C. trachomatis L2 infectious particles, we isolated an AZM-tolerant population of C. trachomatis L2 growing in the presence of twice the MIC after 30 passages with AZM.
Ten individual plaques formed after expansion of the last harvest in the presence or absence of AZM for 2 weeks revealed no mutations in either of the two 23S rRNA genes of the mutants. We also sequenced the rplD and rplV genes, coding for ribosomal proteins L4 and L22, respectively. These genes represent two other known mutational targets for macrolide resistance (12). In all 10 isolates, one single mutation, C196 to A, was found in rplD, creating to a Gln-to-Lys alteration at position 66 in the C. trachomatis L2 ribosomal protein L4. This mutation arose in the population as early as the 23rd passage (in L2AZM#23).
Drug sensitivity phenotypes associated with the Q66K alteration in C. trachomatis L2 ribosomal protein L4. Gln66 in C. trachomatis L2 ribosomal protein L4 corresponds to Gln62 in the homologous E. coli protein and lies within the highly conserved region of L4 that is responsible for binding to 23S rRNA (28, 64). The alteration in ribosomal protein L4 conferred an eightfold decrease in sensitivity for AZM, a fourfold decrease in ERM, JOS, SPI, and CLI sensitivity, and a twofold decrease for VIR sensitivity but did not affect the binding of CLM to ribosomes (Table 4). However, the MICs for the macrolides, with the exception of SPI, and the lincosamide were still lower than 1 µg/ml, suggesting that these antibiotics should still show clinical efficacy against these chlamydial mutants in vivo (24, 53).
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
Although treatment failures (defined as Chlamydia persistence 1 month after treatment) following macrolide therapy have been reported, most reports do not address the role of genetic resistance in the recurrence of chlamydial infections (23, 61). For example, Golden et al. (18) reported a treatment failure rate of 8% among women who reported no new sexual exposures, but that study did not include phenotypic or genotypic susceptibility analysis of the bacteria. Similarly, no genetic analysis was initiated in a patient who failed to clear C. trachomatis infection with JOS (51). In veterinary practice, too, Owen et al. reported that infections with C. felis (C. psittaci feline pneumonitis agent) appeared to be insensitive to AZM in four of five cats studied but did not analyze the possible contribution of mutations to treatment failure (37).
In the present study, high-level AZM resistance arose in C. psittaci 6BC at a frequency of around 1 x 10–8. This frequency is similar to that reported for Mycobacterium avium, which, like C. psittaci, harbors a single rRNA operon (21, 31). We previously noted a similarity in the development of aminoglycoside resistance by both species (7). Here, we found that spontaneous macrolide resistance arose in C. psittaci 6BC from single mutations in the unique 23S rRNA gene resulting in an A2058C, A2059C, or A2059G substitution. This is the first report showing stable, high-level resistance to multiple clinically relevant antibiotics in a chlamydial strain following single point mutations in the 23S rRNA gene. Interestingly, although the A2058G substitution in 23S rRNA is the most common macrolide resistance mutation encountered in bacterial pathogens (12, 60), it was found in none of the 30 independent mutants of C. psittaci 6BC analyzed in this study. This suggests that the A2058G change is deleterious for C. psittaci 6BC. Pfister et al. (42) recently proposed that the fitness cost associated with this particular mutation is dependent on the nature of the adjacent 2057-2611 base pair, as bacteria with A2057-U2611, such as mycobacteria, are more tolerant to the resistance mutation than species with G2057-C2611, such as E. coli, Helicobacter pylori, and Streptococcus pneumoniae. However, C. psittaci 6BC and other Chlamydiaceae display the A2057-U2611 pair as in the mycobacteria. Therefore, one would expect that the biological cost of an A2058G mutation would be tolerated in C. psittaci as it is in mycobacteria. We isolated no such mutants. It is possible that, compared to mycobacteria, C. psittaci 6BC lacks additional genetic factors such as an intragenic or extragenic compensatory mutation to balance the cost of the A2058G mutation (8, 29). Moreover, no A-to-U transition was observed at position 2058 or 2059 in the AZM-resistant C. psittaci mutants, and we did not detect any mutation at positions 2057, 2452, and 2611, although those have been shown to confer low-level macrolide resistance in some other bacteria (60). Because bacteria growing at high density exhibit a low level of antibiotic resistance, we used high concentrations of antibiotic to select for resistant mutants (7). This strategy precluded isolating low-level-resistance isolates of Chlamydia by plaque assay. Furthermore, isolation of chlamydial variants by this technique depends on both growth and cell-to-cell transmission of the organism. In the absence of a genetic system for chlamydiae, we cannot determine at this time if any of these mutations would be associated with a level of AZM resistance too low to allow selection in the plaque assay or if the mutations would be deleterious for the organism.
In contrast to C. psittaci 6BC, we were unable to obtain highly AZM-resistant isolates of C. trachomatis L2 in the plaque assay. We showed previously that the frequencies of resistance to rifampin are on the same order for both chlamydial species (7), suggesting that mutations in the 23S rRNA gene should form at equivalent rates. However, because C. trachomatis harbors two 23S rRNA gene copies, selection of spontaneous AZM-resistant isolates in the plaque assay requires the mutation to be dominant over the wild-type (unmutated) copy. Although a resistant phenotype is apparently codominant in mycobacteria (49), Staphylococcus aureus (47), H. pylori (57), Mycoplasma hominis (38), and Ureaplasma parvum (40), homogenization of the mutation to all ribosomal copies by gene conversion has been linked to high-level macrolide resistance (32, 56) as well as better stability or homogeneity of the resistance phenotype (4, 57, 62). At the present time, we lack the genetic tools to examine the contribution of single mutations in one 23S rRNA gene copy to the level of macrolide resistance in C. trachomatis L2. Interestingly, a double mutation, A2058C T2611C, was recently characterized in clinical isolates of C. trachomatis after AZM treatment failure, suggesting that multiple mutations in the 23S rRNA gene are necessary for expression of high-level resistance (30). However, the percentages of heterozygous and homozygous populations are not clear in that report, as both wild-type and mutated copies of the 23S rRNA were detected in the clinical isolates.
In the present study, we isolated a population of C. trachomatis L2 with an eightfold decrease in AZM susceptibility due to a mutation in rplD coding for ribosomal protein L4. Gln66 (Gln62 in E. coli) lies in a phylogenetically conserved disordered loop that displays many basic residues proposed to be the central element in binding of the protein to the 23S rRNA (28, 64). Replacement of the uncharged Gln with the positively charged Lys is likely to affect the binding of the chlamydial ribosomal protein L4 to the cognate 23S rRNA molecules. Indeed, mutations at residues in this conserved L4 disordered loop were found to affect the overall folding of 23S rRNA in domains II, III, and V, perturbing both the translational activity of ribosomes and the action of antibiotics known to interact with nucleotide residues in the peptidyl transferase center (28, 33, 58). Although rplD mutations linked to low-level macrolide resistance have been selected in vitro in clinically relevant organisms (10, 39, 54), such mutations are frequently found to be paired with additional mutations in 23S rRNA or in rplV (ribosomal protein L22) in clinical isolates (41, 48). Wolter et al. elegantly showed that persistence of macrolide-resistant clinical isolates of Streptococcus pneumoniae resulting from alterations in ribosomal protein L4 is linked to acquisition of secondary mutations that compensate for the defect in bacterial growth as well as increasing the level of macrolide resistance (63). In the present work, we found that the physiological cost associated with the Q66K mutation in C. trachomatis L4 protein could be alleviated without reversion of the drug sensitivity characteristics of the strain (Table 4 and 6). Further experiments are planned to test whether the new "adapted" genetic background would now be more favorable for acquisition of high-level AZM-resistant mutations.
The rate of increase in frequency of resistance to an antibiotic is directly proportional to the efficacy of the drug and the extent of its use and is inversely proportional to the cost that resistance imposes on bacterial fitness, i.e., its rate of infectious transmission and its ability to compete with other strains (3). We analyzed the biological costs associated with point mutations in the 23S rRNA gene conferring macrolide resistance to C. psittaci 6BC by comparing growth of the susceptible parent to that of three isogenic macrolide-resistant variants in the absence of selection. Each mutant was delayed in the formation of infectious particles, as seen by the extended eclipse phase of the development cycle. Additionally, each mutant was severely outcompeted by the wild-type strain at the end of the cycle (CI less than 0.2). These results indicate that mutations in the 23S RNA gene associated with AZM resistance also impose a high physiological cost on C. psittaci. It is worth adding that chlamydial growth conditions in this study were optimized for cell culture, i.e., cycloheximide was added to the medium to create a more favorable environment for bacterial growth by inhibiting host protein synthesis, thus making more nutrients available to the bacteria (20). One might expect even stronger developmental differences between each mutant and the parent strain when they have to compete for survival in the natural host. Therefore, we predict that these mutations would not be maintained in vivo unless compensatory mutations are selected to adapt to the costs of chromosomal antibiotic resistance (8, 27, 29).
It is clear that many interacting factors influence the probability of macrolide resistance development in chlamydial infections. Antibiotic-resistant mutants selected in vivo tend to be ones with low or no fitness cost in vitro (16, 34, 50). Unfortunately, the lack of a small animal model for C. trachomatis L2 precludes testing the consequence of the change in ribosomal protein L4 on the pathogenicity of the strain. In addition, our laboratory strain of C. psittaci 6BC has been passaged in tissue culture for several years, and its current level of virulence in parrots has not been studied. Nevertheless, our results give rise to some general considerations regarding the likelihood of development and persistence of AZM-resistant chlamydial strains. In this study, mutants were readily obtained in vitro for C. psittaci 6BC following mutations in the unique 23S rRNA gene, while only rplD mutants of C. trachomatis L2 with increased tolerance to AZM were isolated. Each mutant strain was attenuated for growth in the absence of selection, suggesting that emergence of AZM resistance in chlamydiae in vivo depends on acquisition of compensatory mutations. Frequent intermittent antimicrobial treatments as a consequence of recurrent infections or repeated mass prophylaxis therapy currently in place to control trachoma in countries where it is endemic (15), as well as for experimental prevention of coronary artery diseases (19), might select for isolates adapted to the high cost of chromosomal resistance. Allowing the C. trachomatis rplD mutant to evolve during serial passages permitted selection of a more fit bacterial population in vitro without changing the antibiotic sensitivity characteristics of the strain. However, we still lack information on whether chlamydiae are evolving in vivo in response to antibiotic selection pressure. It will be important to examine the consequences of drug resistance mutations on the pathogenicity and survival of chlamydiae such as C. muridarum or C. caviae, which can be tested in their natural small mammal hosts (i.e., mice and guinea pigs, respectively). Additionally, follow-up of potential treatment failures with both phenotypic and genotypic susceptibility assays of C. trachomatis isolates will be important in order to understand whether and how antibiotic resistance occurs in the clinical setting. In the meantime, caution should be used in relying upon AZM for treatment of human and veterinary psittacosis.
| ACKNOWLEDGMENTS |
|---|
We thank Nancy E. Adams for technical assistance, Anne-Laure Prunier for her comments and critical reading of the manuscript, and Michael N. Flora and the USUHS Biomedical Instrumentation Center for DNA sequencing and oligonucleotide synthesis services.
The opinions or assertions contained herein are those of the authors and are not to be construed as official or reflecting the views of the Department of Defense or the Uniformed Services University of the Health Sciences.
| FOOTNOTES |
|---|
Published ahead of print on 1 October 2007. ![]()
| REFERENCES |
|---|
|
|
|---|
G alteration mediates ketolide resistance in combination with deletion in L22. Antimicrob. Agents Chemother. 50:3816-3823.
G. Proc. Natl. Acad. Sci. USA 102:5180-5185.
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Clin. Vaccine Immunol. | Clin. Microbiol. Rev. |
|---|---|