Previous Article | Next Article ![]()
Journal of Bacteriology, December 2007, p. 8685-8692, Vol. 189, No. 23
0021-9193/07/$08.00+0 doi:10.1128/JB.01150-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
,
Centre for Infectious Diseases, Royal (Dick) School of Veterinary Studies, The Chancellor's Building, New Royal Infirmary, University of Edinburgh, Edinburgh EH16 4SB, Scotland, United Kingdom,1 Centre for Preventive Medicine, Animal Health Trust, Newmarket, Suffolk CB8 7UU,United Kingdom,2 Department of Veterinary Pathobiology, Faculty of Life Science, University of Copenhagen, Frederiksberg C, Denmark,3 Dermatology Unit, Division of Veterinary Clinical Sciences, Royal (Dick) School of Veterinary Studies, Easter Bush Veterinary Centre, University of Edinburgh, Roslin, Midlothian, EH25 9RG, Scotland, United Kingdom4
Received 20 July 2007/ Accepted 13 September 2007
|
|
|---|
|
|
|---|
Despite its prevalence among many members of the animal kingdom, an understanding of the population genetics of S. intermedius is lacking. Until very recently, studies of the diversity of S. intermedius populations have been limited to phenotypic and molecular typing approaches (1, 3, 5, 12, 16-18, 40). A relatively high degree of phenotypic diversity exists within the species, leading to the identification of different biotypes associated with specific animal hosts (16, 28), and ribotyping studies have shown that strains are typically associated with only a single host species (1, 5). Canine strains were represented by a small number of distinct ribotypes, in contrast to strains from other animals, which generally had greater diversity (17). The level of phenotypic and genotypic diversity observed among S. intermedius isolates has led some investigators to speculate that an S. intermedius group (SIG) consisting of several species or subspecies may exist (5, 28). A closely related species, Staphylococcus pseudintermedius, was recently described (7) and has been isolated from several animal species, including healthy dogs in a veterinary clinic in Japan (32), and from a human infection (38). Very recently, a study of Japanese isolates based on DNA sequences of sodA and hsp60 genes identified three different species, S. intermedius, S. pseudintermedius, and Staphylococcus delphini, among isolates phenotypically identified as S. intermedius and suggested a reclassification of the species (33). However, the distribution of these different species among isolates from outside Japan was not examined.
Worryingly, methicillin-resistant strains of S. intermedius (MRSI) and S. pseudintermedius (MRSP) have recently emerged in veterinary clinics around the world (15, 22, 32, 44), but their evolutionary histories have not been examined. An understanding of the population genetics of S. intermedius is required to provide a much-needed framework for studies of the evolution, pathogenesis, and emerging methicillin resistance of S. intermedius. In order to investigate the population genetic structure of S. intermedius, we developed a multilocus sequencing approach that included five gene loci with a range of predicted nucleotide diversities to facilitate inter- and intraspecies differentiation. Broad new insights into the diversity of S. intermedius populations were obtained.
|
|
|---|
Bacterial growth conditions and genomic-DNA isolation. Staphylococcal strains were grown on tryptic soy agar (Oxoid) at 37°C overnight or in tryptic soy broth (Oxoid) at 37°C overnight with shaking at 200 rpm. Genomic-DNA extraction was carried out with a bacterial genomic-DNA purification kit (Edge Biosystems) according to the manufacturer's instructions. Prior to incubation at 37°C for 10 min, 125 µg/ml lysostaphin (Ambi products) was included.
Gene loci selected for DNA sequence analysis. Five loci, the 16S rRNA, cpn60 (hsp60), tuf, pta, and agrD genes, were selected for DNA sequence analysis. 16S rRNA, cpn60, and tuf have been used previously in single-locus approaches to differentiating staphylococcal species (7, 13, 24, 25), and cpn60, tuf, and pta have been incorporated into multilocus sequence-typing schemes for differentiating strains of several different bacterial species (2, 10, 43). The accessory gene regulator (agr) quorum-sensing system found in the majority of staphylococcal species examined (8) is thought to have diversified along lines of speciation, giving rise to a number of subspecies groups that have the capacity for intra- and interspecies inhibition of virulence (30, 42). We examined the allelic variation of agrD, which encodes the autoinducing peptide (AIP), to further investigate both inter- and intraspecies differentiation.
PCR amplification of gene fragments. Oligonucleotide primers (Sigma-Genosys or Invitrogen) were designed for the 16S rRNA, cpn60, and agrD genes based on gene sequences available for S. intermedius in GenBank (accession numbers Z26897, AF019773, AY965912, and AF346723), and pta primers were designed based on the S. intermedius pta nucleotide sequence (unpublished data). Oligonucleotide sequences specific for the tuf and mecA genes were designed in previous studies (21, 31). The oligonucleotide sequences and predicted PCR product sizes for specified gene fragments were as follows: 16S rRNA (370 bp), 5'-CCT CTT CGG AGG ACA AAG TGA-3' (forward) and 5'-GAC CCG GGA ACG TAT TCA CC-3' (reverse); tuf (500 bp), 5'-CAA TGC CAC AAA CTC G-3' (forward) and 5'-GCT TCA GCG TAG TCT A-3' (reverse); cpn60 (550 bp), 5'-GCG ACT GTA CTT GCA CAA GCA-3' (forward) and 5'-AAC TGC AAC CGC TGT AAA T G-3' (reverse); pta (570 bp), 5'-GTG CGT ATC GTA TTA CCA GAA GG-3' (forward) and 5'-GCA GAA CCT TTT GTT GAG AAG C-3' (reverse); agrD (300 bp), 5'-GGG GTA TTA TTA CAA TCA TTC-3' (forward) and 5'-CTG ATG CGA AAA TAA AGG ATT G-3' (reverse), and 5'-CTC ATG ACT ATT GCA TGC ATC G-3' (reverse; pigeon isolates only); and mecA (310 bp), 5'-GTA GAA ATG ACT GAA CGT CCG ATAA-3' (forward) and 5'-CCA ATT CCA CAT TGT TTC GGT CTA A-3' (reverse). PCR mixtures for the 16S RNA, tuf, cpn60, pta, and mecA genes included 0.2 µM each primer, 0.025 U/µl Taq polymerase (Promega), 1.5 mM MgCl2, 0.2 mM deoxynucleoside triphosphates (Promega), and 1 µl genomic DNA template in 96-well PCR microplates (Axygen/Thistle Scientific). The thermocycler program for the 16S RNA, tuf, cpn60, and pta genes consisted of an initial denaturation for 2 min at 95°C, followed by 30 cycles of denaturation for 1 min at 95°C, annealing for 1 min at 55°C, and an extension for 1 min at 72°C, and a final extension for 7 min at 72°C. The mecA thermocycler program consisted of an initial denaturation for 5 min at 95°C, followed by 30 cycles of denaturation for 45 s at 95°C, annealing for 45 s at 52°C, and an extension for 1 min at 72°C, and a final extension for 5 min at 72°C. For agrD amplification, the reaction mixtures included 0.5 µM each primer, 0.02 U/µl Vent DNA polymerase (New England Biolabs [NEB]), 1x ThermoPol Reaction Buffer (NEB), 0.2 mM deoxynucleoside triphosphates (Promega), and 1 µl genomic-DNA template. The thermocycler conditions were an initial denaturation for 2 min at 94°C, followed by 35 cycles of denaturation for 15 s at 94°C, annealing for 30 s at 45°C, and extension for 1 min at 72°C, with a final extension for 7 min at 72°C (36). PCR products were purified by incubation with exonuclease I and Antarctic phosphatase (NEB) at 37°C for 15 min, followed by inactivation at 80°C for 15 min.
DNA sequencing. Sequencing reactions were carried out using 5 µl purified PCR-amplified DNA (approximately 50 to 90 ng) plus 1 µl sequencing primer (3.2 pmol/µl) with the BigDye Terminator v3.0 Ready Reaction cycle sequencing kit (Applied Biosystems) according to the manufacturer's instructions. Each sequencing reaction included 25 cycles with denaturation for 30 s at 95°C, annealing for 20 s at 50°C, and extension for 4 min at 60°C. The nucleotide sequence was determined with a 3730 DNA analyzer (Applied Biosystems). For the 16S rRNA, cpn60, tuf, and pta genes, forward and reverse sequencing reactions were carried out with independently amplified PCR products to rule out the possibility of Taq polymerase-generated errors. For agrD, forward and reverse sequencing reactions were carried out from a single PCR product generated with Vent DNA polymerase (NEB), which contains proofreading activity and has a much lower predicted error rate than Taq polymerase.
DNA sequence and molecular evolutionary analyses.
DNA sequences were assembled using BioEdit Sequence Alignment Editor software (http://www.mbio.ncsu.edu/BioEdit/bioedit.html) and the Staden package (35). Nucleotide and amino acid sequences were aligned using AlignX in Vector NTI Advance10 (Invitrogen). Phylogenetic analyses were carried out with MEGA v.3.1 software (23). The neighbor-joining method was applied to construct phylogenetic trees using the Kimura two-parameter model, while the degree of statistical support for the nodes on the minimum evolution tree was evaluated by examining their percent recovery in 1,000 resample trees by the bootstrap test. For estimating events of recombination, RDP software, which includes the programs GENECONV, BOOTSCAN, MAXIMUM
2, CHIMAERA, and SISTER SCANNING, was used (27). The index of association standardized (IAS) between the different gene loci was calculated using the LIAN program (version 3.1; Department of Biotechnology and Bioinformatics, University of Applied Sciences, Weihenstephan, Germany) (http://adenine.biz.fh-weihenstephan.de/lian_3.1/). The value of the IAS would be expected to be zero for linkage equilibrium when recombination events occur frequently. If the IAS value differs significantly from zero (P < 0.05), recombination should be rare.
eBURST analysis. Predicted lines of evolutionary descent and clonal complexes in our collection of isolates were identified using the eBURST algorithm (http://eburst.mlst.net). Sequence types (STs) were included in the same group if they shared four of the five gene loci with at least one other ST within the group. Subgroups were defined by the existence of at least three single-locus variants.
Nucleotide sequence accession numbers. The DNA sequences generated in this study were deposited in GenBank, accession no. EU157195 to EU157715.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Summary of nucleotide sequence variation for the SIG and for the S. pseudintermedius, S. delphini, and S. intermedius phylotypes
|
![]() View larger version (25K): [in a new window] |
FIG. 1. Phylogenetic tree constructed by the neighbor-joining method with 1,000 bootstrap replicates of concatenated 16S rRNA gene, tuf, cpn60, and pta sequences. (A) Radiating tree indicating the three major phylotypes, with species type strains indicated by triangles. (B) Branch style tree indicating the host and geographic origins of all isolates. Methicillin-resistant, mecA-positive isolates (dots); human isolates (squares); previously identified S. pseudintermedius isolates (diamonds); and species-type isolates (triangles) are indicated, along with the agr type in roman numerals in parentheses. Bootstrap values over 50% are indicated.
|
S. pseudintermedius has a largely clonal population structure.
As mentioned above, trees constructed with sequences from each housekeeping locus appeared to be generally in good agreement, suggesting that recombination has not played a major role in the evolution of these genes and consistent with a clonal population structure (see Fig. S1 in the supplemental material). To further investigate the population structure, we calculated the degrees of linkage disequilibrium within the whole SIG and within the S. pseudintermedius and S. delphini species (Table 2). The value of the IAS for the SIG as a whole was 0.1510 (P < 0.001), and for the S. pseudintermedius and S. delphini species it was 0.0877 (P < 0.001) and 0.3699 (P < 0.001), respectively, indicating that there is linkage disequilibrium within each group tested (Table 2). In addition, the RDP2 program was used to search for evidence of recombination among the selected gene loci (27). For each locus, none of the recombination detection methods, GENECONV, BOOTSCAN, MAXIMUM
2, CHIMAERA, and SISTER SCANNING, detected any evidence for recombinant sequences among the isolates examined (data not shown). Taken together, these data indicate a largely clonal population structure for the SIG as a whole and for the S. pseudintermedius and S. delphini species.
|
View this table: [in a new window] |
TABLE 2. Analysis of linkage disequilibrium by calculation of IAS for the SIG and for the S. pseudintermedius and S. delphini phylotypes
|
Identification and distribution of agr types among the SIG of closely related species. DNA sequencing of the agrD locus revealed the presence of four predicted AIP variants among the strains examined, including a novel fourth variant that has not been described previously (designated type IV) which is specific to 27% of all strains (Table 3). The AIP variants may correspond to distinct agr interference groups (8, 30, 36, 42). All staphylococcal species examined to date have been shown to encode AIP peptides that are unique to each species. In contrast, the three closely related species examined in the current study had AIP peptide variants in common, suggesting the existence of a conserved agr quorum-sensing system. All four AIP variants (I to IV) were encoded among the S. pseudintermedius isolates; AIP variants I, II, and IV were identified among the S. delphini group of isolates; and the S. intermedius type strain and related pigeon isolates all encoded AIP type I (see Table S1 in the supplemental material).
|
View this table: [in a new window] |
TABLE 3. Amino acid sequences of the predicted agrD-encoded AIPs identified among the strains examined
|
MRSP strains have evolved by multiple mecA gene acquisitions by different clones. The presence of the mecA gene in 16 of 105 isolates was detected by PCR (see Table S1 in the supplemental material). Sequencing of the mecA gene in representative isolates revealed a high degree of homology with mecA of Staphylococcus aureus origin (95 to 100%; data not shown). These strains had previously been identified as MRSI in the centers where they were isolated and had oxacillin MIC levels of >256 µg/ml for all strains except one (strain 13), which had an MIC of 0.75 µg/ml. All strains previously identified as MRSI belonged to the S. pseudintermedius phylotype and are therefore reclassified as MRSP. MRSP genotypes are distributed widely across the diversity identified within the S. pseudintermedius phylogenetic tree, including one isolate, 3279, that is genetically indistinguishable from methicillin-sensitive isolates (ST29), suggesting horizontal transfer of the mecA gene (Fig. 1). The existence of multiple MRSP clones is supported by eBURST analysis, which shows that MRSP isolates belong to five distinct STs (ST29, ST68, ST69, ST70, and ST71) that differ at two to four of the five loci examined, strongly suggesting that they do not share a very recent ancestor (Fig. 2). Taken together, these data suggest that the mecA gene has been acquired multiple times by different S. pseudintermedius strains. Of note, 9 out of 10 MRSP isolates from five different centers in Sweden and Germany belong to ST71, indicating that a common clone may predominate in North and Central Europe (see Table S1 in the supplemental material). There was no sharing of STs among European and U.S. MRSP isolates (see Table S1 in the supplemental material).
![]() View larger version (9K): [in a new window] |
FIG. 2. Schematic diagram of the clonal relatedness of S. pseudintermedius STs predicted by eBURST analysis. Each ST is represented by a black dot, the size of which is proportional to the number of isolates of that ST. Blue and yellow dots denote the predicted group founder and subgroup founders, respectively. Single-locus variants are linked by lines. STs that include MRSP strains are indicated by red circles. STs that are not linked by branches to other STs do not share at least four identical alleles with any other ST.
|
|
|
|---|
The third major phylotype found among the isolates previously identified as S. intermedius was associated with several different host species, including horse, camel, and pigeon, and was phylogenetically aligned with the S. delphini type strain isolated from a dolphin (39). A previous DNA-DNA hybridization analysis confirmed that the S. delphini type strain represented a distinct species compared to the recently described S. pseudintermedius and the S. intermedius type strain (7). Taken together, these data suggest that S. delphini is commonly misidentified as S. intermedius and may be more clinically important than was previously thought. In fact, very few studies have reported the identification of S. delphini strains since the original description of the species, but one study in Norway isolated S. delphini from a case of bovine mastitis, extending the broad host range observed in the current study (4). Of note, Sasaki et al. indicated that more than one species may exist among isolates genetically allied with S. delphini (33).
The accessory gene regulator (agr) is conserved throughout the staphylococci but has diverged along lines that appear to parallel speciation within the genus (8, 42). This divergence has given rise to a novel type of interstrain and interspecies cross-inhibition that likely represents an essential aspect of staphylococcal biology and may be a predominant feature of the evolutionary forces that have driven it. In order to further investigate the diversity of SIG populations, we examined the diversity of the agrD locus, which encodes the AIP (8, 20, 36, 42). Sung et al. identified three agrD alleles corresponding to agr specificity groups encoding different AIP variants among 20 strains of S. intermedius isolated from dogs in a single veterinary hospital in the United Kingdom (36) and demonstrated biological activity for two of them. Here, we identified the same three predicted AIPs, in addition to a novel fourth AIP, which was encoded by approximately 27% of strains (Table 3). Of note, all of the predicted AIP variants of the SIG contain a serine in place of the conserved cysteine found in the AIPs of the other staphylococcal species analyzed to date, resulting in a lactone rather than a thiolactone ring (20). S. aureus populations are also divided into four distinct groups based on agr allelic variation. Interference in virulence gene expression caused by different S. aureus agr groups has been suggested to be a mechanism for isolating bacterial populations and a fundamental basis for subdividing the species (42). All previously examined staphylococcal species encode AIPs that are unique to each species (8). However, we found that agr alleles were shared among the three closely related species of the SIG, indicating that a common quorum-sensing capacity has been maintained despite species differentiation in largely distinct ecological niches.
A previous study of agr evolution in S. aureus reported that genotypes were not associated with more than one agr type, indicating that agr radiation preceded clonal diversification and that recombination has played a very limited role in the distribution of agr diversity (42). The study concluded that the species was phylogenetically structured according to agr group (42). More recently, Robinson et al. provided evidence that the S. aureus species could be divided into two subgroups that both contain multiple clonal complexes and agr groups (30). The authors proposed that recombination events had resulted in the sharing of agr groups between the two subspecies but concluded that recombination of agr has not occurred very frequently within S. aureus populations, as agr and clone tree topologies within the two subspecies were in agreement (30). In contrast, the association of different agr alleles with strains of S. pseudintermedius of identical genotype identified in the current study suggests that assortive recombination has frequently contributed to the distribution of agr alleles among S. pseudintermedius populations. The markedly different topology of the phylogenetic tree constructed with agrD sequences compared to trees based on the other four gene loci is consistent with this theory (Fig. 1). Overall, we have found that S. pseudintermedius has a largely clonal population structure, but recombination appears to have played an important role in the distribution of agr alleles within the S. pseudintermedius species and the SIG as a whole. The sharing of agr alleles between different species of staphylococci has not been previously observed. This discovery indicates that agr differentiation has not occurred strictly along lines that parallel speciation. The lack of an association between agr type and SIG species, host, disease, and clinical or geographic origin identified in the current study leads to the question of what selective pressure is driving agr diversification. The basis for agr diversity and the importance of its biological activity in the SIG remain to be elucidated.
Recently, an increasing number of episodes of S. intermedius infections that were refractory to treatment with methicillin have been reported (15, 22, 32, 44). Here, we found that the methicillin-resistant SIG strains examined are all classified as S. pseudintermedius and have evolved by multiple acquisitions of the mecA gene by different S. pseudintermedius clones. Of note, a common clone was identified among isolates from five different centers in Sweden and Germany, indicating the existence of a widespread successful clone in Northern and Central Europe. MRSP clones were not shared between Europe and North America, indicating geographic restriction and probably reflecting the very recent emergence of methicillin-resistant strains. The identification of the common methicillin-resistant clones in the current study means that surveillance can be carried out to monitor MRSP clonal dissemination and strategies for the control of MRSP infections can be targeted to the most widespread clones.
Taken together, these data have resolved the population genetic structure of the SIG, resulting in broad new insights into the ancient and recent evolution of this important group of animal pathogens.
Published ahead of print on 28 September 2007. ![]()
Supplemental material for this article may be found at http://jb.asm.org/. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»