|Home | About | Journals | Submit | Contact Us | Français|
Fluoroquinolones are an important class of antibiotics for the treatment of infections arising from the gram-positive respiratory pathogen Streptococcus pneumoniae. Although there is evidence supporting interspecific lateral DNA transfer of fluoroquinolone target loci, no studies have specifically been designed to assess the role of intraspecific lateral transfer of these genes in the spread of fluoroquinolone resistance. This study involves a comparative evolutionary perspective, in which the evolutionary history of a diverse set of S. pneumoniae clinical isolates is reconstructed from an expanded multilocus sequence typing data set, with putative recombinants excluded. This control history is then assessed against networks of each of the four fluoroquinolone target loci from the same isolates. The results indicate that although the majority of fluoroquinolone target loci from this set of 60 isolates are consistent with a clonal dissemination hypothesis, 3 to 10% of the sequences are consistent with an intraspecific lateral transfer hypothesis. Also evident were examples of interspecific transfer, with two isolates possessing a parE-parC gene region arising from viridans group streptococci. The Spain 23F-1 clone is the most dominant fluoroquinolone-nonsusceptible clone in this set of isolates, and the analysis suggests that its members act as frequent donors of fluoroquinolone-nonsusceptible loci. Although the majority of fluoroquinolone target gene sequences in this set of isolates can be explained on the basis of clonal dissemination, a significant number are more parsimoniously explained by intraspecific lateral DNA transfer, and in situations of high S. pneumoniae population density, such events could be an important means of resistance spread.
Fluoroquinolones were originally introduced in the mid-1980s principally for the treatment of gram-negative bacterial infections. Subsequent development of fluoroquinolones was targeted to also be effective against gram-positive cocci. Much of the reason for this development interest was related to the increasing resistance of the important respiratory tract pathogen Streptococcus pneumoniae to other antibiotics, particularly beta-lactams. Fluoroquinolone resistance in S. pneumoniae remains relatively low, but it is nonetheless increasing, particularly in some geographic areas, such as Hong Kong (16). Resistance to fluoroquinolones in gram-positive cocci, such as S. pneumoniae, is thought to be mainly due to spontaneous mutations that occur in target loci due to errors of polymerase during replication (17). The acquisition of new genetic material via lateral gene transfer (LGT; also referred to as recombination) is generally assumed to be of limited importance in fluoroquinolone resistance in gram-positive bacteria (17).
Several recent studies have specifically set out to examine the role of interspecific lateral DNA transfer of fluoroquinolone target loci from viridans group streptococci to S. pneumoniae. Ferrandiz et al. (12) and Balsalobre et al. (1) present evidence for interspecific lateral DNA transfer involving the parE-parC gene region from viridans group streptococci to S. pneumoniae. Balsalobre et al. (1) concluded that 11% of the 46 isolates examined were putative interspecific recombinants and suggest that viridans group streptococci act as donors of fluoroquinolone resistance genes to S. pneumoniae via horizontal transfer. Recently, de la Campa et al. (5) reported on the quinolone resistance-determining regions (QRDRs) of parE, parC, and gyrA from a set of 75 ciprofloxacin-resistant isolates from Spain, 5 of which are of viridans group streptococcus origin. In contrast, another study concluded that the interspecific lateral transfer of fluoroquinolone target loci was not occurring at all (30). Thus, although there is evidence for the lateral transfer of fluoroquinolone target loci from viridans group streptococci to S. pneumoniae, the frequency and importance of the phenomenon in the spread of resistance remain somewhat uncertain.
If interspecific lateral transfer of fluoroquinolone target loci is known to occur, then it is likely that intraspecific lateral transfer among S. pneumoniae isolates is also occurring. Indeed, it has been shown that the frequency of homologous recombination in S. pneumoniae decreases with the sequence divergence between donor and recipient (20, 25). Furthermore, it is clear that penicillin binding proteins (PBPs) are transferred both intraspecifically (4, 22) and interspecifically (7, 14). Nonetheless, we are not aware of any demonstration of intraspecific DNA exchange involving the fluoroquinolone target loci of S. pneumoniae clinical isolates. Undoubtedly, part of the reason for this is the fact that intraspecific events are much more difficult to detect than interspecific events, since differences in the donor and recipient sequences can be much more subtle and require an approach that concomitantly indexes both the evolutionary history of the isolates and that of the resistance genes. Several recent publications have commented on the importance of clonal dissemination in the spread of fluoroquinolone resistance in S. pneumoniae (15, 28, 31), while other studies from North America report a wider genetic diversity, with fewer isolates clearly associated with the international clones (21, 40). A complete understanding of the role of clonal spread is not possible without knowledge of the role of intraspecific lateral DNA transfer of fluoroquinolone target loci. At present, we are not aware of any published studies that attempt to address intraspecific lateral DNA transfer of S. pneumoniae fluoroquinolone target loci.
A classic means of detecting lateral DNA transfer is to assess branching congruence or incongruence in comparisons of a control evolutionary history (the best estimate of the true evolutionary history for the organisms in question) to that derived from a particular locus (gene history) obtained from the same biological entities (for a review, see reference 23). Incongruent histories support lateral DNA transfer, and congruent histories do not. This approach is analogous to that used in previous studies of intraspecific lateral transfer of PBPs in S. pneumoniae, except that the earlier studies did not compare reconstructed evolutionary histories but instead genotyped the isolates and assessed whether different clones had identical PBP gene sequences (horizontal transfer) or were in fact indistinguishable based on overall genetic relatedness and their PBP sequence (clonal spread) (4, 22). Evolutionary histories provide an additional level of information regarding relationships and, consequently, the prospect of a more thorough picture of clonal spread versus horizontal transfer. However, in the case of S. pneumoniae the reconstruction of a control evolutionary history is complicated by recombination and by the relatively low level of sequence variation present in housekeeping genes between different isolates (10). Although both these variables hamper phylogenetic reconstruction, they are much less of a problem with evolutionary networks, which are well suited to the representation of relationships between sequences of relatively low divergence, typical of population data, and allow the incorporation of the nonbifurcating genealogical information that can be associated with population-level divergences (3, 10). It is also possible to at least partially ameliorate the problems of recombination. A large number of methods for the detection of recombinant sequences are now available, and some of these have been applied to bacterial population sequence data (13). Once recombinant sequence regions are detected, they can be excised before the reconstruction of evolutionary histories. Thus, we believe that it is possible to reconstruct control evolutionary histories for S. pneumoniae isolates by employing network reconstruction of relatively large amounts of housekeeping sequence (with the multilocus sequence typing [MLST] data as a subset) from each isolate after the sequence has been carefully screened for putative recombinants. Such a control history can then serve as a meaningful framework for comparison against evolutionary networks of other genes of interest, in this case, fluoroquinolone target loci.
It is important to try and assess the existence and ultimately estimate the frequency of intraspecific lateral transfer of fluoroquinolone target loci because it could be a very important means of resistance spread in situations in which the population density of S. pneumoniae is high. Children represent such a situation; studies of day care and pediatric chronic care centers have found that children are often colonized with high densities of pneumococci in the nasopharynx (26, 27, 39). At present, the potential pediatric indications for fluoroquinolones include chronic ear infection, invasive gastrointestinal infection, complicated urinary tract infections, and pseudomonal bronchopulmonary complications in cystic fibrosis (26). There is, however, increasing pressure on regulatory authorities to approve fluoroquinolones for more general use in the pediatric population (26). Therefore, frequent intraspecific lateral transfer of fluoroquinolone target loci between S. pneumoniae isolates, combined with the approval of fluoroquinolones for more general use in the pediatric population, could represent an important set of circumstances for the development of fluoroquinolone resistance.
The purpose of this study is to use molecular evolution principles and techniques in an attempt to gain insight into the frequency with which S. pneumoniae fluoroquinolone target loci may be spread via intraspecific lateral DNA transfer. The evolutionary history of a diverse set of S. pneumoniae clinical isolates from the Alexander Project collection (11, 19) is derived from network reconstructions of an expanded MLST data set that includes a concatenated 8-kb sequence alignment of housekeeping genes, with putative recombinants excluded. This control history is then assessed against evolutionary networks of each of the four fluoroquinolone target loci from the same set of isolates. Congruent histories support a clonal dissemination hypothesis for a given fluoroquinolone target locus, and discordant histories suggest lateral DNA transfer.
The susceptibility data and descriptions of the tests for the Alexander Project have been published elsewhere (11, 19). Our choice of isolates from the Alexander Project collection involved a random subsample of 42 fluoroquinolone-resistant isolates (MIC ≥ 8 μg/ml for either ciprofloxacin or ofloxacin from the total of 170 fluoroquinolone-resistant isolates), 4 fluoroquinolone-intermediate isolates (ciprofloxacin and ofloxacin MICs = 4 μg/ml), and 14 susceptible isolates (ciprofloxacin and ofloxacin MICs < 4 μg/ml) and spanned the years 1992 to 2000. Two of the fluoroquinolone-resistant isolates were subsequently determined to have a parE-parC gene region that arose from an interspecific DNA exchange event involving viridans group streptococci and were not included in the evolutionary reconstructions. Table Table11 includes the details for the 60 isolates used in this study. The majority of isolates (n = 44) were nasopharyngeal in origin, the origins of 11 were unknown, and 5 were from blood. Colony PCR was performed with 45 μl of Invitrogen High Fidelity PCR Supermix per colony. The control histories were reconstructed from a concatenated set of housekeeping gene sequences; the resulting alignment of nine loci for the 58 S. pneumoniae clinical isolates, prior to any exclusions, was 9,473 bp long. This included all or nearly all of the sequences of the recP (transketolase), gki (glucose kinase), gdh (glucose-6-phosphate dehydrogenase), xpt (xanthine phosphoribosyltransferase), spi (signal peptidase I), and aroE (shikimate dehydrogenase) genes used for MLST (8); we did not sample ddl (the usual seventh MLST locus) because previous studies have shown that it is frequently recombinant (9). In addition to the sequencing of larger fragments of the customary MLST loci, we also included the genes neighboring aroE, spi, and xpt, which are aroB (3-dehydroquinate synthase), rnhb (RNase HII), and pbux (xanthine permease), respectively. The complete aroE and aroB genes were amplified with primers aroE+1f (5′-CACTGCGGATGTGACTGGTTCGA-3′) and aroE+2r (5′-CCCTCAATAATAGCTGTTAGACGGGG-3′). The complete spi and rnhb genes were amplified with primers spi+1f (5′-GATAGAAGAAGAGGCTGAGATTGGT-3′) and spi+1r (5′-CAATCTCACGGCTGAGCTGAGTT-3′). A fragment encompassing all of the xpt gene and most of the pbux gene was amplified with primers xpt+1f (5′-GTCTATGATACCACTACAACGGGA-3′) and xpt+1r (5′-CGGCATTGAGGACAATAGCGAGT-3′). We amplified gki (all but 33 bp of the 3′ end), gdh (all but 27 bp of the 3′ end and 40 bp of the 5′ end), and recP (all but 76 bp of the 3′ end) using the following respective primer pairs: gki1f (5′-CACGCAAACCTTTGCATAAGTGA-3′) and gki2r (5′-ACAAGTGATGCTGCTCCGATAAC-3′), gdh2f (5′-TGTTACAATTTTCGGTGCGAGTG-3′) and gdh2r (5′-TCTAAGCGACCATCTTGACGATA-3′), and recP1f (5′-CTTTGGAGGATTTCCGTATGTTG-3′) and recP2r (5′-CTGCCAATACTTTTGGTGCTGGG-3′).
In addition to the housekeeping genes, we sampled the four fluoroquinolone target genes from the same set of 58 isolates: parE (topoisomerase IV subunit B), parC (topoisomerase IV subunit A), gyrA (DNA gyrase subunit A), and gyrB (DNA gyrase subunit B). The entire gyrB, gyrA, and parE genes were amplified with the following respective primer pairs: gyrAp1 (5′-ACGTTTTAGTGGTTTAGAGG-3′) and gyrAp4 (5′-TGAGCATCCCATTGAGAAGTTAGA-3′), gyrBp1 (5′-ATTGGCACTGTATGGTATCAC-3′) and gyrBp4 (5′-AGCGCCTCTAATCTCCCCTCGTT-3′), and parEp1 (5′-GCCGATAACAGCCAAAAGTCC-3′) and parEp4 (5′-AATTCGTCCCGTGTTCAGTTAC-3′). A fragment encompassing all of parC except 6 bp of the 3′ end was amplified with primers parCp1 (5′-TTTATGGGCTTTGTATCTTATGTC-3′) and parCp4 (5′-GTTTTCCCAATCAATCGTTACTGT-3′).
The PCR products were purified with a QIAGEN QIAquick PCR purification kit, according to the manufacturer's instructions. The products were sequenced on both strands by using Applied Biosystems machines 3730, 3700, and 3100 and Applied Biosystems Big Dye.
Homologous DNA sequence alignments were generated by using ClustalX, version 1.83 (38), with any necessary additional manual editing in Genedoc, version 22.214.171.124. The nine housekeeping genes were concatenated into a single alignment, and then the program PLATO (13) was used to determine if any portion of the alignment was evolving in an anomalous fashion due to either recombination or selection. Recombination can be difficult to detect with the low levels of sequence divergence typical of S. pneumoniae housekeeping genes, and PLATO is likely to miss some of the more subtle instances of recombination (33). However, we have found that it performs very well in detecting sequences known to be frequently recombinant in S. pneumoniae (ddl, for example ), and it has a low propensity to detect recombination when it is not present (33). Thus, PLATO may miss some subtle instances of recombinant sequences, but at the same time, it has a high likelihood of detecting the majority of recombinant sequences in our alignments. Since PLATO does not tolerate polytomies, identical sequences were removed from the housekeeping data set prior to PLATO runs and were then added back to the alignment for further phylogenetic analyses. Multiple PLATO analyses were performed to ensure convergence by using the following parameters: 10 steps for the sliding window, 1,000 replications of Monte Carlo simulation, and an HKY model of sequence evolution. Any regions of the housekeeping alignment that PLATO found to be statistically anomalous were removed. The individual genes within this concatenated alignment encompassed the following regions: aroE = positions 1 to 855, aroB = positions 856 to 1923, gdh = positions 1924 to 3337, gki = positions 3338 to 4272, pbux = positions 4273 to 5489, recP = positions 5490 to 7390, rnhb = positions 7391 to 8276, spi = positions 8277 to 8891, and xpt = positions 8892 to 9473.
The program TCS (3) was used to implement the method of Templeton et al. (37) for the reconstruction of statistical parsimony networks for both the housekeeping alignment and each of the fluoroquinolone target gene alignments. To avoid ambiguities in haplotype assignment, all missing positions were removed from each alignment prior to TCS analysis (37). In order to provide a further means of verifying our clonal group designations, phylogenies were also reconstructed. For the final concatenated housekeeping alignment, as well as the four resistance gene alignments, the best-fitting model of sequence evolution and the corresponding values for the rate matrix, shape of the gamma distribution, and proportion of invariant sites were estimated by the program MODELTEST (32). Phylogenies were reconstructed by the maximum-likelihood (ML) and neighbor-joining (NJ) methods implemented in PAUP*4.0b (36). For the ML method, starting trees were obtained by neighbor joining, with nearest-neighbor interchange used as the branch-swapping algorithm. Bootstrap support values were obtained with 1,000 replicates for NJ analyses and 100 replicates for ML analyses. In addition, Bayesian phylogenetic reconstruction was performed using Mr. Bayes (18) with 3,000,000 generations, a sampling frequency of every 100 generations, GTR + gamma + invariants, four Markov chains, random starting trees, and a burn-in of 2,000,000 generations. Multiple Bayesian analyses were performed to ensure convergence. Phylogenetic statistical comparisons were made by the conservative Shimodaira-Hasegawa (SH) test (34, 35).
The sequences obtained in this study have been deposited in GenBank under the following accession numbers: DQ174780 to DQ175475 and DQ176454 to DQ176511.
Within the concatenated alignment of nine housekeeping genes encompassing 9,473 bp, PLATO identified several anomalous evolving regions, which were excluded before phylogenetic analysis. These regions, with their accompanying PLATO Z values (Z values greater than 4.2526 were judged to be significant; the value indicated is the lowest from the various independent runs) included the following: positions 238 to 244, Z = 11.7949; positions 1874 to 1881, Z = 7.0216; positions 1971 to 1977, Z = 13.71133; positions 2350 to 2354, Z = 4.4019; positions 2792 to 2796, Z = 4.5986; positions 3083 to 3092, Z = 6.4793; positions 3124 to 3128, Z = 5.7274; positions 3185 to 3189, Z = 19.5443; positions 3212 to 3233, Z = 11.5082; positions 3313 to 3317, Z = 4.6513; positions 3790 to 3796, Z = 5.8212; positions 3844 to 3935, Z = 19.8386; positions 4617 to 4638, Z = 4.8284; positions 4922 to 4927, Z = 5.3199; positions 5337 to 5343, Z = 6.4382; positions 5993 to 5999, Z = 9.1872; positions 6229 to 6233, Z = 5.8387; positions 6303 to 6307, Z = 13.2928; positions 6672 to 6676, Z = 8.4217; positions 7157 to 7220, Z = 4.5825; positions 7267 to 7271, Z = 6.8452; positions 7561 to 8555, Z = 18.5434; positions 9033 to 9038, Z = 17.4671; and positions 9206 to 9210, Z = 4.3690. After the exclusion of these 1,310 characters the resulting alignment was 8,163 bp in length. Nearly all of the rnhb locus, as well as the 5′ end of spi, were excluded (rnhb is immediately upstream of spi) due to their anomalous character. Individual gene trees involving rnhb provided clear evidence for several strongly supported conflicting nodes in comparison to trees constructed from the 8,163-bp housekeeping alignment. It seems likely, therefore, that this 993-bp stretch encompassing most of rnhb and the 5′ end of spi is frequently recombinant.
The control statistical parsimony network (3, 37) reconstructed from the 8,163 bp of concatenated housekeeping sequence data depicts two very distinct groups (Fig. (Fig.1).1). The main body of the network includes an assortment of isolates from around the world with various antibiotic resistance phenotypes, including isolates susceptible to all classes of antibiotics, as well as fluoroquinolone-nonsusceptible isolates (Table (Table1).1). A second group, disconnected from the rest of the network, includes only fluoroquinolone-nonsusceptible isolates, collected in 1999 and 2000 from China and the United States, that were so dissimilar from the other isolates that too many steps (>19) were required to accommodate their connection with the rest of the network. Based on our comparison of the necessary housekeeping sequence data to the MLST database (8), the latter group is quite clearly the globally distributed, multidrug-resistant clone known as Spain 23F-1 (MLST allelic profile 4-4-2-4-4-1-1; we amplified and sequenced the MLST portion of ddl from representative isolates of this clone to verify the ddl allele). The control phylogeny (Fig. (Fig.2)2) depicts a very similar arrangement, with two very distinct clades and a very long branch separating Spain 23F-1 from the rest of the isolates. To our knowledge this is the first indication that the Spain 23F-1 clone may be highly genetically distinct from much of the rest of S. pneumoniae. The reason that our data reveal this previously unrecognized difference appears to lie in the high level of sequence divergence of the aroE and aroB genes of Spain 23F-1 from the sequences of the aroE and aroB genes of the rest of the isolates. The customary MLST approach (8, 29) includes an approximately 400-bp piece of aroE, as well as similarly sized pieces of six other loci, whereas our approach involves sequence data derived from the complete aroE and aroB loci, concomitant with seven other complete (or very nearly complete) loci, of which the MLST sequence data are a subset. The fact that we have sequenced all of aroE, as well as aroB, results in the long branch leading to Spain 23F-1. The elimination of aroE and aroB from this data set eliminates the long branch, and reduction of the data set down to only the MLST sequence set does not reveal this long branch. Explanations for this difference in Spain 23F-1 aroE and aroB are not clear. Reconstruction of the phylogeny of various Streptococcus aroE and aroB sequences available in GenBank does not support an interspecific lateral transfer argument; the Spain 23F-1 aroE and aroB sequences tend to be monophyletic with other S. pneumoniae sequences and merely represent a divergent allele. However, it is possible that the putative donor sequence is simply not yet available in GenBank. Selection is a possibility, although many of the synapomorphies for Spain 23F-1 are synonymous substitutions. aroE codes for shikimate dehydrogenase and aroB codes for 3-dehydroquinate synthase, both of which are involved in the aromatic biosynthetic pathway; however, the selection cause-and-effect issues for such loci within a specific strain are far from clear.
We used the information from both the housekeeping network (Fig. (Fig.1)1) and phylogeny (Fig. (Fig.2)2) to demarcate several clonal groups (where a clonal group can be either a clone or isolates from a few closely related clones), the members of which are colored similarly in Fig. Fig.33 in order to track the evolution of their fluoroquinolone target genes. These colored clonal groups are based on a consensus of closely adjacent haplotypes apparent in the housekeeping network (Fig. (Fig.1)1) and strongly supported monophyletic groups with short (or no) internodes in the case of the phylogeny (Fig. (Fig.2).2). In general, our clonal groups correspond very closely to the clonal designations based on MLST criteria, which suggests that two isolates are different clones if they differ at three or more alleles of the seven MLST loci (29). An exception to this would include our “light blue” clonal group, where isolate 132.1US00 differs from the other two isolates by at least four alleles and would therefore be a different clone, based on MLST criteria. However, our network (Fig. (Fig.1)1) and phylogeny (Fig. (Fig.2)2) very clearly associate this isolate with the other two, and therefore, we assign isolate 132.1US00 to the same clonal group as isolates 68.1Wa00 and 94.1Nd00. Our purpose here is to follow the history of the fluoroquinolone resistance-associated genes of isolates which share a recent common ancestry; all of the similarly colored isolates in the network (Fig. (Fig.1)1) and phylogeny (Fig. (Fig.2)2) fit that criterion. Adoption of the MLST criteria for clonal designation (with the exclusion of ddl) results in the conservative estimate that the housekeeping network (Fig. (Fig.1)1) includes a minimum of 24 clones, 17 of which are fluoroquinolone nonsusceptible. These clones are scattered throughout the network and occasionally include fully susceptible isolates as part of the clone, suggesting that fluoroquinolone nonsusceptibility has evolved on numerous different occasions in this set of isolates. However, approximately 50% of the fluoroquinolone nonsusceptible isolates can be accounted for by only a few clones (three or four). Nonsusceptible clones occasionally come from the same geographic location (e.g., the “red” clone from Japan with isolates 319.1JA99S, 319.2JA99S, 294.1JA99S, and 319.1JA00S) but more often include isolates from different geographic locations (e.g., the “purple” [Taiwan 19F-14] clone with isolates 458.1JA00S, 526.1Tx99S, and 543.1Ca99S). In addition to Spain 23F-1, it is possible to attach several of our isolates to some of the Pneumococcal Molecular Epidemiology Network (PMEN) international clones. For example, nonsusceptible isolates from the Taiwan 23F-15 and the Taiwan 19F-14 clones are represented (six identical alleles; ddl unknown), with others at least closely related (vary at one locus; ddl unknown) to Poland 23F-16, England 14-9, Spain 9V-3, and Tennessee 23F-4 (Fig. (Fig.11 and and2).2). Fluoroquinolone resistance has previously been reported for all of these clones (2, 15, 28). In our set, Spain 23F-1 is the clone with the largest number of fluoroquinolone nonsusceptible isolates (n = 13). Ten nonsusceptible isolates are members of clones currently uncharacterized by the PMEN.
In the process of collecting the comparative sequence data from the Alexander Project collection isolates, we identified several interspecific lateral DNA transfer examples involving the parE and parC gene regions of viridans group streptococci similar to those described previously (1, 12). Sequencing of the parE gene of isolates 434.1FR99S and 591.1Ca99S yielded a gene sequence similar to that of S. mitis, and thus, we sequenced the intergenic region between parE and parC of these two isolates in order to compare to earlier reports involving interspecific transfer of this region (1). Comparative analysis indicated that both our isolates, similar to those described by Balsalobre et al. (1), have a complete viridans group streptococcus ant gene immediately downstream of parE. Between ant and parC, two of the isolates of Balsalobre et al. (specifically, isolates 4589 and 3870) possess an approximately 875-bp insert (homologous to the sequence of S. mutans GTP pyrophosphokinase) which is not shared with our isolates, although one of ours (isolate 434.1FR99S) does have a 97-bp homologous piece of the 3′ end of this insert. Balsalobre et al. (1) reported 5 cases (out of 46) of S. pneumoniae strains with interspecific DNA exchange involving the parE, parC gene region, and we have found an additional 2 such cases (out of 60). Of these seven examples, there are five different genomic arrangements for this region, suggesting that interspecific recombination in this gene region may not be that rare of an occurrence. Recently, de la Campa et al. (5) reported on the QRDRs of parE, parC, and gyrA from a set of 75 isolates, 5 of which were of viridans group streptococcus origin. These various examples serve to highlight the fact that foreign DNA does appear to be playing a role in shaping the parE- parC gene region within S. pneumoniae, although the degree to which this is the case appears to be variable between studies and may be related to the environmental situation from which the isolates were sampled (30).
It seems highly likely that if interspecific lateral gene transfer of fluoroquinolone target loci is occurring between S. pneumoniae and viridans group streptococci, intraspecific recombination is also occurring. This view stems from the fact that the frequency of homologous recombination in S. pneumoniae decreases with the sequence divergence between donor and recipient (25) and that high population densities of S. pneumoniae are common in certain environments and situations (26, 27, 39). Homologous recombination involving fluoroquinolone loci could involve recombination points within or outside the target locus. Recombinant points either outside the gene or at least proximal to the 5′ and 3′ ends involving a donor and a recipient that differ in sequence, result in a different branching position for that isolate, relative to the control history, in networks reconstructed from sequence data for that particular fluoroquinolone target locus. Recombinant points within the gene may not result in a different branching position in the resulting evolutionary history because the recombinant piece may be relatively small and may not carry sufficient evolutionary signal to alter its position compared to the control history. Thus, a convincing alternative placement for a fluoroquinolone target locus relative to the control history results in a conservative estimate of the extent of intraspecific lateral gene transfer for this set of isolates. Our criterion for lateral transfer is the convincing disruption of a clonal group (same color in the control network) in networks derived from the fluoroquinolone target genes. This disruption can take two forms: (i) dispersion of one or more control clonal group members across the fluoroquinolone target gene network (category 1) or (ii) inclusion of a novel (different colored) isolate(s) in a control clonal group on the fluoroquinolone target gene network (category 2). For category 1, convincing disruption refers to the isolates of the same color in highly disparate positions (separated by a large number of substitutions and no longer nearest neighbors) on the fluoroquinolone gene network. For category 2 the most convincing evidence is from isolates that have identical or nearly identical haplotypes on a fluoroquinolone gene network but that are not closely related on the control housekeeping network. This conservative approach means that category 1 transfers are evaluated by using only the colored clonal groups and not singleton clones. A focus only on clonal groups for category 1 transfers will underrepresent the absolute number of LGT events for the entire data set. A further point of significance with regard to these definitions is that for category 1 the putative donor and recipient involved in the LGT are not determinable, whereas in category 2 they generally are.
A comparison of the control histories to each of the gene histories derived for the fluoroquinolone target loci indicates that there are regions of both branching similarity and differences (Fig. (Fig.1;1; Fig. 3a to d). For example, there is a group of four Japanese fluoroquinolone nonsusceptible isolates (indicated in red) from 1999 and 2000 which are of the same haplotype for all fluoroquinolone target loci as well as the housekeeping genes (Fig. (Fig.1;1; Fig. 3a to d). This sort of arrangement for the control history (Fig. (Fig.1)1) clearly supports these isolates as members of the same clone, and the same arrangement for each of the fluoroquinolone target loci (Fig. 3a to d) indicates that the fluoroquinolone target genes of these isolates were acquired through common ancestry and descent. A similar example involves the susceptible isolates 0.49SP94 and 0.57SW99, which have identical haplotypes for all statistical parsimony networks (Fig. (Fig.1;1; Fig. 3a to d). Indeed, the majority of fluoroquinolone gene sequences from the isolates of the colored clones group together and thus support the view of clonal dissemination of the fluoroquinolone target loci. Other cases support lateral DNA transfer involving one or more isolates of a clonal group and one or more fluoroquinolone target genes but never all four. For example, there is an international fluoroquinolone nonsusceptible clone apparent in the control history which includes two isolates from Japan, one isolate from Texas, and another isolate from Italy (indicated in orange in Fig. Fig.1).1). These isolates have the identical haplotype for both gyrase genes, indicating no lateral transfer of these loci. However, the history underlying their parE genes is somewhat different (Fig. (Fig.3a).3a). All four of these “orange” isolates have the same parE sequence; however, unlike the isolates with the control history, four additional isolates (isolates O.47SP00, 215.1Oh00, 132.1NY98, and 68.1PR00) also share this same haplotype and are very widely separated from the “orange” clone on the housekeeping network (Fig. (Fig.1).1). This is an example of our category 2 (see above) for lateral transfer and indicates that members of the “orange” clone have acted as donors in lateral DNA transfer events involving these other lineages. An example of our category 1 lateral transfer would involve the “blue” clonal group, composed of isolates 136.1PL98, 95.1BR98, 0.34PL96, and 72.1In00 (Fig. (Fig.1);1); the gyrB network clearly dissociates 136.1PL98 from the rest of the clones, with placement at virtually opposite ends of the network, indicating that 136.1PL98 has received a recombinant gyrB locus (Fig. (Fig.3d).3d). Another example involves the “pink” clonal group (isolates 95.1Ma00 and 215.1Oh00) and the “light green” clonal group (isolates 0.19IS99 and 67.1FR92), which have identical gyrA sequences (Fig. (Fig.2b)2b) but which are very dissimilar on the control network (Fig. (Fig.1).1). This suggests lateral transfer of gyrA involving these two clones, but the directionality cannot be assessed.
One of the advantages of using phylogenetic trees is that there are statistical tests available that allow one to make comparisons between topologies (branching arrangement) and thus assess whether the highest likelihood topology is significantly different from a constrained topology that reflects an alternative hypothesis of evolutionary history. A similarly agreed upon set of tests are presently not available for network comparisons. Using the conservative SH test (34, 35), we compared the highest-likelihood trees reconstructed from the fluoroquinolone target sequence data against trees reconstructed from the same data, but with the constraint that they had to contain the 11 clonal groups present in the control phylogeny. If the fluoroquinolone target loci in this set of isolates were vertically inherited, then they should conform to the clonal groups depicted on the control histories and the SH test run for each gene would not be significant. If the fluoroquinolone target genes are laterally transferred, then their evolutionary history will not conform to the clonal groups on the control phylogeny and the SH test will be significant. The result indicates that for all four fluoroquinolone target genes for this set of 58 isolates the null hypothesis of no difference between the constrained or clonal phylogeny and the highest-likelihood phylogeny was rejected, and thus, we accept the alternative hypothesis of lateral DNA transfer: gyrA, difference in -ln L = 14.8934 (P = 0.024); gyrB, difference in -ln L = 106.8542 (P = 0.000); parC, difference in -ln L = 56.8522 (P = 0.009); and parE, difference in -ln L = 82.7417 (P = 0.000). It is important to realize that this test does not incorporate any information regarding interclonal group relationships but, instead, merely assesses clonal group composition. Furthermore, the resulting gene phylogenies and statistical tests, employing principles very different from those involved in the network reconstruction, are in complete agreement with the conclusions derived from the networks. Thus, this phylogenetic statistical approach serves to corroborate our network-based conclusions regarding the lateral transfer of fluoroquinolone target loci.
Networks, on the other hand, have the conceptual advantage that they do not force the data into a tree-like (bifurcating) structure; however, it can be much more difficult to make comparisons between them. This is arguably somewhat exacerbated in our situation by the difference in the length of the data set used to reconstruct the control history and that used for each of the fluoroquinolone target loci. More specifically, in a comparison of networks, the confidence that one can place in a particular LGT event is related to the number of steps (nucleotide substitutions) that distinguish the isolates in question on the different networks. In a further effort to ensure that our network-based conclusions on transfer are as conservative as possible, we reconstructed additional networks specifically designed to assess whether the larger number of sequence positions in our housekeeping data (8,163 bp) compared to the number of sequence positions of individual fluoroquinolone target loci (1.6 to 2.3 kb, depending on the gene) were in any way biasing our conclusions on transfer. This was accomplished by creating multiple (n = 20), randomly reduced (jackknifing) housekeeping data sets (using the 8,163-bp data set as input) that match the number of sequence positions for each of the respective fluoroquinolone target loci, repeating the network reconstruction for each of those smaller data sets, and then evaluating whether each individual transfer hypothesis was still supported with these data sets of comparable size. All apparent instances of either category 1 or category 2 transfer were also supported by using these more reduced data sets and in some instances arguably provided more convincing evidence. For example, the adjacent isolates comprising the turquoise clone (Taiwan 23F-15) (isolates 292.1AT96 and 588.1JA99) are separated by 10 steps in the 8,163-bp housekeeping network (Fig. (Fig.1)1) and by 17 steps on the gyrA network (Fig. (Fig.3b);3b); however, in housekeeping networks based on the same number of sequence positions as gyrA, these two isolates are separated by only two steps (n = 20; mean = 2.45; standard deviation = 1.39). Thus, random reduction of the data set simply reinforced the view that these two isolates have gyrA sequences that are very dissimilar (17 steps removed) from that expected from their evolutionary history (same clone; approximately 2 steps).
Another possible explanation for fluoroquinolone target loci from different clones with similar sequences is convergent molecular substitutions brought about by similar antibiotic selection pressure. If that were the case, then the majority of such substitutions should be amino acid changing in nature, and thus, evaluation of the lateral DNA transfer hypothesis by using only synonymous substitutions should yield a different result. Statistical parsimony networks reconstructed by using only synonymous substitutions tended to yield the same conclusions regarding lateral transfer events, indicating that convergent molecular selection due to antibiotic selective pressure is not the more parsimonious explanation. A conservative estimate of the number of lateral transfer events of fluoroquinolone target loci involving the lineages in this data set is approximately 17 (summarized in Table Table2).2). The number of examples of lateral transfer for each of the loci seems roughly similar except for parC, for which the current data set supports fewer instances (Table (Table22).
Many of the fluoroquinolone nonsusceptible isolates in this analysis carried substitutions in the QRDRs of the four loci, recognized as being important in conferring resistance to this class of antibiotics. This included the following substitutions and proportion of isolates, for each of the respective loci: parC, D78A/N (4.5%) and S79F/Y (54.5%), D83Y/G/N (11.4%); parE, D435N (13.6%), P454S (2.3%), and E474K (2.3%); gyrA, S81F/Y (52.3%) and E85K (2.3%); and gyrB, R379L (2.3%). Approximately 23% of the nonsusceptible isolates had wild-type parC; such proportions of fluoroquinolone nonsusceptible isolates without first-step parC mutations are in rough agreement with those from other studies (24). The vast majority of the isolates without parC mutations had mutations in gyrA and parE, with the most common parE mutations being L290F and A326V.
An important and interesting group to consider in the context of the lateral transfer of fluoroquinolone target genes is the Spain 23F-1 clone. We do not find any examples of these isolates “losing” their fluoroquinolone genes through recombination with another allele; note that “dark green” isolates are never dispersed in networks derived from their fluoroquinolone genes (Fig. (Fig.3).3). Instead, we find isolates from outside this clone with the same fluoroquinolone target genes as the Spain 23F-1 group (our category 2 indication of lateral transfer). This is true for all the fluoroquinolone target genes except parC (Fig. (Fig.3c).3c). These lateral transfer events involve isolates throughout Europe and the United States covering a range of years. The fact that we see Spain 23F-1 as a donor of fluoroquinolone resistance genes and never as a recipient suggests the possibility that selection is somehow inhibiting recombination in Spain 23F-1. This might seem a bit counterintuitive, since Spain 23F-1 may have acquired its fluoroquinolone resistance character through a recombinational history. However, now that it has developed resistance to the current crop of fluoroquinolones, the selective pressure is different, and it might well be more beneficial to retain these particular alleles than lose them through recombination. It has been argued elsewhere that bacterial mismatch repair, which affects rates of recombination and mutation, might be in a state of evolutionary flux in bacteria (6). This could result from some environmental situations in which it is more beneficial to incorporate new DNA more efficiently, whereas in other situations it might be a better strategy to keep the genome more stable, thereby selecting for efficient mismatch repair. Antibiotics represent such a scenario; new antibiotics represent an enormous selective pressure, but for resistant clones the pressure becomes one of not losing resistance determinants through an undesirable recombination event. One would, however, need a much larger isolate sample size to properly evaluate this recipient hypothesis. Whatever the cause-effect explanations, Spain 23F-1, one of the most widely dispersed, highly successful, multidrug-resistant S. pneumoniae clones (29), is not only a health menace because of these characteristics, but also because its fluoroquinolone resistance determining genes are being laterally transferred to other lineages. Undoubtedly, this is also true for other major and minor clones; but the ubiquitous nature of the Spain 23F-1 clone, concomitant with the lateral transfer of the fluoroquinolone resistance associated genes arising from members of this clone, means that Spain 23F-1 is a doubly important source of fluoroquinolone resistance spread.
Our data and analysis indicate that fluoroquinolone resistance associated genes are being laterally transferred between S. pneumoniae isolates (minimum of 17 LGT events in a set of 58 isolates). This in turn suggests a potentially important means of resistance spread; however, the degree to which it is actually responsible for the spread of resistance remains a question. The histories of the fluoroquinolone target genes for the isolates evaluated in this data set are more often consistent with clonal dissemination than with lateral DNA transfer. For example, if we confine ourselves to a consideration of just our designated clonal groups, in the case of parE, of the 38 nonsusceptible isolates within these groups, 32 of them possess parE sequences which are consistent with a clonal dissemination hypothesis for this locus, 4 are consistent with lateral DNA transfer, and 2 are ambiguous. Roughly similar results are apparent for the other loci (clonal, LGT, and ambiguous): gyrA, 33, 3, and 2, respectively; gyrB, 34, 2, and 2, respectively; and parC: 33, 1, and 4, respectively.
Although the majority of fluoroquinolone target gene sequences in this S. pneumoniae data set can be explained on the basis of clonal dissemination, there are still a significant number which are more parsimoniously explained by intraspecific lateral DNA transfer; and in situations of high S. pneumoniae population density, such LGT events could be a very important means of resistance spread. The accurate assessment of this phenomenon is not trivial; but in the process of attempting to do so, we have made a concerted effort to be highly conservative in our identification of LGT loci. The degree to which LGT has played and is continuing to play a role in the spread of resistance to fluoroquinolones and other classes of antibiotics deserves further attention. Studies specifically assessing the relative frequency of intraspecific lateral DNA transfer of resistance conferring loci in different environmental situations would be of particular interest.
We thank the GSK US-based sequencing facility, under the management of Ganesh Sathe, for all their efforts on this project.