Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, January 2005, p. 366-379, Vol. 49, No. 1
0066-4804/05/$08.00+0 doi:10.1128/AAC.49.1.366-379.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Laboratory of Molecular Genetics,1 Biomathematics Group, Instituto de Tecnologia Química e Biológica, Universidade Nova de Lisboa, Oeiras, Portugal,2 Laboratory of Computational Genomics,3 Laboratory of Microbiology, The Rockefeller University, New York, New York4
Received 26 November 2003/ Returned for modification 6 February 2004/ Accepted 11 September 2004
|
|
|---|
|
|
|---|
Several microbial surface proteins (adhesins) mediate the adherence of S. aureus to host proteins, such as fibrinogen and fibronectin. These plasma proteins coat indwelling medical devices, and the ability of the bacteria to adhere to the deposited proteins is believed to be an important factor in the pathogenesis of wound and foreign body infections (25). In S. aureus, the adhesin genes which have been characterized in detail include clf (45, 48) and fnb (36, 69), which encode the fibrinogen- and the fibronectin-binding proteins, respectively. Although the clumping factor genes, clfA and clfB, encode proteins whose structural organizations resemble those of other cell wall-anchored surface proteins, one distinctive feature is the presence of a repeat region, the R domain, that connects the cell wall-spanning domain to the A domain, which contains the ligand-binding site. It has been proposed that the R domain may function as a "stalk" that allows the exposure of the A domain at the bacterial surface for ligand interaction or that acts in the attachment of the protein to the cell wall (44). The R domain, both in ClfA and in ClfB, is mainly composed of repeats of the dipeptide serine-aspartate encoded by an 18-bp variable repeat, GAYTCNGAYTCNGAYAGY (45, 48). Owing to its potential for genetic variability, the R domain of the clumping factor genes was one of the target regions for sequencing in this study.
The fibronectin-binding protein genes, fnbA and fnbB, encode proteins which, in addition to fibronectin, also bind to fibrinogen (78). Besides their role in the pathogenesis of medical device-related infections, the fibronectin-binding proteins of S. aureus have been shown to be sufficient for the invasion of human cell lines (43, 71) and factors that contribute to the colonization of the mammary gland in a mouse model of mastitis (10). The fnb genes are in tandem in the S. aureus chromosome and share sequence homology, in particular at the 3' end: in the region containing the D repeats, in the wall-spanning WR repeats, and in hydrophobic membrane-spanning domain M (36). Originally, the fibronectin-binding activity was placed in the D region of FnbA (69), but recently it has been shown that this protein has multiple fibronectin-binding regions which include the B and C domains (35, 43). In FnbB, an additional fibronectin-binding region has been observed upstream to the D domain (36). In this study the DNA region containing the D, W, and M domains in both fnbA and fnbB has been sequenced.
Expression of these and other virulence factors is, at least in part, differentially regulated by the agr operon (58) and several other global regulator loci, which is likely to be important for the adaptation and survival of the microorganism in the host. In our study, in addition to regions of the clf and fnb genes, the following regions were also sequenced in the agr operon: (i) the 5' region of agrC; (ii) the DNA region between the P2 and P3 promoters containing the binding site for SarA (11), another transcription regulator related to virulence factor expression; and (iii) the complete DNA sequence coding for the effector molecule, RNAIII.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Molecular data for the representative collection of MRSA and MSSA strains used in this study
|
|
View this table: [in a new window] |
TABLE 2. Oligonucleotide primers used for amplification and sequencing of virulence-related loci in this study
|
(ii) fnb genes. A 3' fragment of the fnb genes that included the repeated D region, the wall-spanning W region, membrane-spanning domain M, and the translation stop codon was amplified with one forward primer common to both fnbA and fnbB (primer fnbF1) and one reverse primer specific for each gene located in the intergenic region (primer fnbAR1) or downstream of fnbB (primer fnbBR1).
(iii) agr operon. The 5' variable region of the agrC gene was amplified with primers agrCF1 and agrCR1. For strains HDE288, BK2464, BM18, COB3, JP1, CN1, N315, E2104, and E3001, however, for which no amplifications were obtained, primers agrCF2 and agrCR2 were used instead, with an extension time of 3' at 72°C and sequencing with primer agrCR2 only. (After sequencing, these strains were all shown to belong to agr group II.) The DNA fragment containing the complete RNAIII sequence and the region between the P2 and P3 promoters was amplified with primers RNAIIIF1 and RNAIIIR1; however, for the strains mentioned above, as well as strains HAR24 and E1410, the same forward primer and primer RNAIIIR2 were used under the same amplification conditions.
DNA sequencing. PCR products were purified with a Wizard PCR Preps purification system (Promega) and used as templates for automated sequencing with a BigDye Terminator cycle sequencing kit (Applied Biosystems) run on an ABI Prism 3700 DNA analyzer at the DNA Sequencing Resource Center, The Rockefeller University. Sequence data analysis was performed with the DNAStar package (Lasergene).
Representation of clf sequences. The following encoding was devised in order to represent the unexpectedly high degree of variation in the units composing the R domains of the clf sequences. We represented each of the 4 nucleotides by a three-dimensional floating-point vector in a way that all pairwise distances between these points were equal; we thus obtained for the 18-nucleotide repeat units of the form GAYTCNGAYTCNGAYAGY a set of 54-dimensional floating-point vectors, with 1 vector corresponding to one occurrence of the pattern in a clf sequence. We included in our analysis all occurrences of a repeat which differed by, at the most, four positions from the consensus sequence presented above. While this encoding can be used to specify any sequence of nucleotides precisely, it is impossible to visualize. In order to try and achieve a graphic representation of the R-domain sequences, we used the method of singular value decomposition to represent the points in a space with fewer dimensions and color coded the coordinates of the lower-dimensional space. We projected the 54-dimensional space onto a 3-dimensional space and used these three coordinates to color code each occurrence of the repeat in terms of hue, brightness, and saturation. Since hue is the most striking feature of color, it was used for the first coordinate, which reflects the highest variance of the scatter of points. We mapped the second coordinate to brightness and the third coordinate to saturation. As a result of this representation, repeats which were close to each other in the 54-dimensional space remained close to each other in the 3-dimensional space and displayed similar shades of color, and, thus, differences in color reflect differences in the nucleotide sequence of the repeat units. The sequences of the clf genes were displayed as strings of colored segments representing the repeat units of the R domain and black fragments representing sequences which did not follow the repeat unit consensus pattern defined above.
Construction of dendrograms from sequence data. Dendrograms were generated from the amino acid sequence data for the fnb genes on the basis of the percent identity of optimal pairwise alignments obtained with the program lalign (29) and with a BLOSUM-50 score matrix and default gap opening and extension penalties of 14 and 4. Sequences were clustered by the unweighted pair group method with arithmetic means (UPGMA). The dendrogram for the DNA sequences of agrC was generated from a distance matrix obtained from a multiple alignment by use of the CLUSTALX program (75), and sequences were also clustered by UPGMA.
PFGE. PFGE of SmaI digests of chromosomal DNA from the strains shown in Table 1 was performed as described previously (12). Relatedness among the PFGE profiles was evaluated with Bionumerics software (version 3.0; Applied Maths, Ghent, Belgium). The dendrogram was generated from a similarity matrix calculated with the Jaccard coefficient, and patterns were clustered by UPGMA.
MLST and spa typing. MLST was performed as described previously (22), with the exception that primer arcCF2 (5'-CCT TTA TTT GAT TCA CCA GCG-3') (15) was used. MLST alleles and STs were identified by using the MLST database, available at http://www.mlst.net. Molecular typing based on the sequence of the polymorphic region of protein A (spa typing) was performed as described previously (67).
SCCmec typing.
SCCmec types were determined by a multiplex PCR strategy which establishes a specific amplification pattern for each structural type (54). The exceptions were strain CA04, whose SCCmec type was previously assigned by PCR analysis of the ccrAB genes and the mec complex (16, 31, 42), and strain MW2, whose SCCmec element has been fully sequenced (6). In strains of SCCmec type IV, as determined by the multiplex strategy, the ccrAB allele type was also determined by PCR with primer sets specific for alleles 1, 2, and 3 and a control set of primers specific for the ccrAB locus, as described elsewhere (31, 52). According to the criteria defined by these previous investigators, SCCmec types are defined by the combination of the ccrAB allele with the genetic organization of the mecA regulon: the presence (class A mec, mecI-mecR1-mecA) or the absence (class B mec, IS1272-
mecR1-mecA) of the mecA transcription repressor (the mecI gene), information also provided by the multiplex strategy. SCCmec types I and IV are negative for mecI and have ccrAB alleles 1 and 2, respectively; and SCCmec types II and III are mecI positive and have ccrAB alleles 2 and 3, respectively. Since SCCmec type IV is defined by the multiplex strategy by the presence of only two bands (one corresponding to the dcs region and the other corresponding to the mecA gene), all SCCmec type IV strains were further characterized by ccrAB allele determination for confirmation purposes. Most SCCmec type IV strains in this collection characterized by the multiplex strategy had ccrAB allele 2. The exceptions were strains HDE288 and HDE1, which had ccrAB allele 4, as described previously (55); and although they were originally assigned to SCCmec type IV, they are now designated IV-like. In strains CPS22, CPS68, and ICP5011, which are variants of SCCmec type III (mecI positive), ccrAB alleles were also characterized: strains CPS22 and ICP5011 have ccrAB allele 3, while strain CPS68 was nontypeable for the ccrAB allele.
Comparison of virulence loci sequencing with typing methods. Discriminatory power was measured with Simpson's index of diversity, which calculates the probability that two unrelated strains sampled from the test population will be placed into different typing groups (30). Cross-classification concordance levels between clf sequence-based clustering and the results of known typing methods for this collection of isolates were determined by comparing pairwise matches and calculating the percentage of classification agreement (37).
|
|
|---|
![]() View larger version (59K): [in a new window] |
FIG. 1. Clumping factor sequences represented by color-coded repeats of the R domain. (A) clfA R-domain sequences 1 to 13 are unique sequences representative of the following strains: sequence 1, strains E2125, E2453, 10395, and E213 (Archaic clone); strains HSJ216, HU25, and PLN104 (Brazilian clone); and strains CPS22, CPS68, and ICP5011 (Portuguese clone); sequence 2, strain TAW9 (Hungarian clone); sequence 3, strains HU101 and HUSA304 (Hungarian clone); sequence 4, strains HPV107, PER34, and BK1953 (Iberian clone); sequence 5, strain MW2; sequence 6, strains BK2464 and CN1 (NY/Japan clone); sequence 7, strains HDE1 and HDE288 (Pediatric clone); sequence 8, strain N315; strains BM18 and COB3 (Pediatric clone); and strains JP1 and E3001 (NY/Japan clone); sequence 9, strain E2104 (NY/Japan clone); sequence 10, strain HAR22 (EMRSA-15 clone); sequence 11, clone EMRSA-16 (www.sanger.ac.uk/Projects/S_aureus/); sequence 12, strains PLN49 and E3812 (Berlin clone); sequence 13, strain CA04 (Berlin clone). (B) clfB R-domain sequences 1 to 21 are unique sequences representative of the following strains: sequence 1, strains E2125 and E2453 (Archaic clone); sequence 2, strains 10395 and E213 (Archaic clone); sequence 3, strains HPV107, PER34, and BK1953 (Iberian clone); sequence 4, strain N315 and strain HDE288 (Pediatric clone); sequence 5, COB3 (Pediatric clone); sequence 6, strain CN1 (NY/Japan clone); sequence 7, strain BK2464 (NY/Japan clone); sequence 8, strain JP1 (NY/Japan clone); sequence 9, strain E3001 (NY/Japan clone); sequence 10, strain BM18 (Pediatric clone); sequence 11, strain E2104 (NY/Japan clone); sequence 12, strain MW2; sequence 13, strains HSJ216 and PLN104 (Brazilian clone); sequence 14, strains HU101 and HUSA304 (Hungarian clone); sequence 15, strain HU25 (Brazilian clone); strain TAW9 (Hungarian clone); and strains CPS22, CPS68, and ICP5011 (Portuguese clone); sequence 16, strain HAR24 (EMRSA-16 clone); sequence 17, strain E1410 (EMRSA-16 clone); sequence 18, strain HAR22 (EMRSA-15 clone); sequence 19, strain PLN49 (Berlin clone); sequence 20, strain CA04 (Berlin clone); sequence 21, strain E3812 (Berlin clone). Lines perpendicular to the sequences delimit the 3'-end eight-repeat region (navy blue lines), which groups the sequences into classes circumscribed by rectangles and the minimum amount of additional sequence which is necessary to differentiate sequences within each class (turquoise lines). Regions of homology are underlined in instances in which the 3'-end eight-repeat region is interrupted.
|
|
View this table: [in a new window] |
TABLE 3. Discriminatory abilities of typing methods and sequencing results for virulence-related loci
|
Pairwise comparisons of translated fnb sequences resulted in identities between 85.2 and 96% for fnbA and 86.1 and 94.6% for fnbB when MRSA strains with distinct genetic backgrounds were analyzed. For both genes the most similar sequences were those of strain MW2 and those of the Archaic, Iberian, Brazilian, Hungarian, and Portuguese clones; and the most divergent ones were those of the Berlin clone and the NY/Japan clone (in fnbA) or strain MW2 (fnbB) (Fig. 2A and 3A). As expected, the fnb sequences of MSSA strains with genetic backgrounds similar to those of MRSA strains were nearly identical to (identity, 98.7 to 99%) or the same as the corresponding sequences of MRSA strains. Similar to what was observed for the clf genes, seven unique fnbB sequences defined six classes of clonal types: Archaic-Iberian, Pediatric-NY/Japan, EMRSA-15, EMRSA-16, MW2, and Berlin (Fig. 3A). These sequence types have an equivalent unique sequence in fnbA; EMRSA-15 and EMRSA-16 clones, however, do not appear to have this gene (Fig. 2A).
![]() View larger version (41K): [in a new window] |
FIG. 2. fnbA amino acid sequences derived from nucleotide sequences encoding the D, W, and M domains. Unique sequences 1 to 5 represent the following strains: sequence 1, strains E2125, E2453, 10395, and E213 (Archaic clone); strains HPV107, PER34, and BK1953 (Iberian clone); strains HSJ216, HU25, and PLN104 (Brazilian clone); strains HU101, HUSA304, and TAW9 (Hungarian clone); and strains CPS22, CPS68, and ICP5011 (Portuguese clone); sequence 2, strain N315; strains HDE288, BM18, and COB3 (Pediatric clone); and strains BK2464, JP1, CN1, and E3001 (NY/Japan clone); sequence 3, strain E2104 (NY/Japan clone); sequence 4, strains PLN49, CA04, and E3812 (Berlin clone); sequence 5, strain MW2. (A) Dendrogram based on percent identities of fnbA amino acid sequences; (B) fnbA amino acid sequences aligned with the corresponding published sequence for strain 8325-4 (GenBank accession no. J04151) with the CLUSTALX program (75); D1 to D5, repeats of the D domain; WR1 to WR7, repeats of the repetitive region of the cell wall-spanning (W) domain; WC, nonrepetitive region of the W domain; M, membrane-spanning domain (69). The GIDFVED and PIVP motifs are underlined.
|
![]() View larger version (45K): [in a new window] |
FIG. 3. fnbB amino acid sequences derived from nucleotide sequences encoding the D, W, and M domains. Unique sequences 1 to 7 represent the following strains: sequence 1, strains E2125, E2453, 10395, and E213 (Archaic clone); strains HPV107, PER34, and BK1953 (Iberian clone); strains HSJ216, HU25, and PLN104 (Brazilian clone); strains HU101, HUSA304, and TAW9 (Hungarian clone); and strains CPS22, CPS68, and ICP5011 (Portuguese clone); sequence 2, strains N315; strains HDE288, BM18, and COB3 (Pediatric clone); and strains BK2464, JP1, CN1, E2104, and E3001 (NY/Japan clone); sequence 3, strain HAR22 (EMRSA-15 clone); sequence 4, strain HAR24 (EMRSA-16 clone); sequence 5, strain E1410 (EMRSA-16 clone); sequence 6, strains PLN49, CA04, and E3812 (Berlin clone); sequence 7, strain MW2. (A) Dendrogram based on percent identities of fnbB amino acid sequences; (B) fnbB amino acid sequences aligned with the corresponding published sequence of strain 8325-4 (GenBank accession no. X62992) with the CLUSTALX program (75); D1 to D4, repeats of the D domain; WR1 to WR5, repeats of the repetitive region of the cell wall-spanning (W) domain; WC, nonrepetitive region of the W domain; M, membrane-spanning domain (36). The PIVP motif is underlined.
|
agr. Three DNA regions in the agr operon were sequenced: the 5' region of agrC, the region between the P2 and P3 promoters, and the region coding for RNAIII. On the basis of the agrC 5' sequence, all strains in this study could be grouped into previously identified agr groups. Most clones belonged to agr group I: the Archaic, Iberian, Brazilian, Portuguese, Hungarian, EMRSA-15, and Berlin clones had agrC sequences with 98.3 to 99.8% identities to the prototype sequence for agr group I. Interestingly, in strain HUSA304 (Hungarian clone) the tnp gene from IS256 was inserted at nucleotide 472 from the agrC translation start site, which raises issues about the functionality of the agr receptor protein in this strain. The remaining clones belonged to groups II (NY/Japan and Pediatric clones and strain N315) and III (EMRSA-16 clone and strain MW2) (Fig. 4).
![]() View larger version (10K): [in a new window] |
FIG. 4. Dendrogram based on a multiple alignment of agrC 5' DNA sequences performed with the CLUSTALX program (75). Unique sequences 1 to 13 represent the following strains: sequence 1, strain HUSA304 (Hungarian clone); sequence 2, prototype agrC 5' variable sequence for agr group I (50); sequence 3, strains E2125, E2453, 10395, and E213 (Archaic clone), strains HPV107, PER34, and BK1953 (Iberian clone); strains HSJ216, HU25, and PLN104 (Brazilian clone): strains CPS22, CPS68, and ICP5011 (Portuguese clone); and strain HU101 (Hungarian clone); sequence 4, strain TAW9 (Hungarian clone); sequence 5, strain HAR22 (EMRSA-15 clone); sequence 6, strains PLN49, CA04, and E3812 (Berlin clone); sequence 7, prototype agrC 5' variable sequence for agr group IV (GenBank accession no. AF288215); sequence 8, strain HAR24 (EMRSA-16 clone), prototype agrC 5' variable sequence for agr group III (GenBank accession no. AF001783); sequence 9, strain E1410 (EMRSA-16 clone); sequence 10, strain MW2; sequence 11, strains BK2464, JP1, CN1, E2104, and E3001 (NY/Japan clone) and strains HDE288, BM18, and COB3 (Pediatric clone); sequence 12, prototype agrC 5' variable sequence for agr group II (GenBank accession no. AF001782); sequence 13, strain N315.
|
![]() View larger version (31K): [in a new window] |
FIG. 5. Nucleotide sequences of the agr interpromoter region were aligned with the published sequence for a derivative of strain 8325 (GenBank accession no. X52543). Unique sequences 1 to 5 represent the following strains: sequence 1, all strains of the Archaic, Iberian, Brazilian, Hungarian, NY/Japan, and Pediatric clonal types (Table 1) and strains N315 and MW2; sequence 2, strains CPS22, CPS68, and ICP5011 (Portuguese clone); sequence 3, strain HAR22 (EMRSA-15 clone); sequence 4, strains HAR24 and E1410 (EMRSA-16 clone); sequence 5, strains PLN49, CA04, and E3812 (Berlin clone). The annotation is according to Novick et al. (50), Morfeldt et al. (47), and Chien and Cheung (11). Arrows underline direct repeats of the P2 and P3 promoters. The putative 35 and 10 boxes and the SarA-binding site are also underlined.
|
![]() View larger version (71K): [in a new window] |
FIG. 6. Comparison of PFGE patterns of the representative epidemic MRSA strains and the other MRSA and MSSA strains used in this study (Table 1). The dendrogram was generated from a similarity matrix calculated with the Jaccard coefficient, and patterns were clustered by UPGMA. Band comparisons were performed with a tolerance of 1.06. The scale on the top of the dendrogram represents similarity. Shaded branches are below the cluster cutoff value calculated with Bionumerics software.
|
|
|
|---|
Our studies documenting the clone- and strain-specific sequence variations in S. aureus genes, which are important for interactions with the human host, may also provide useful information for future studies on the mechanisms responsible for the superior epidemicity and geographic dominance of globally spread MRSA clones.
Sequence variations in clf genes. The genetic relatedness of the MRSA clones evaluated in this study was classified from the sequences for the R domains of the clumping factor genes. The R-domain sequences themselves provided two levels of analysis: clustering into six classes of clones on the basis of the last eight 3' repeats of clfA (Fig. 1) and identification of individual clones within each cluster or even specific strains of each clone when the sequence was read farther upstream (i.e., to the left from the navy blue perpendicular line markers dissecting the sequences in Fig. 1). The results of the same analysis performed with clfB were concordant with the data for clfA albeit more discriminating; a seventh class was identified, and a larger number of specific sequence types was observed for individual strains from each clone. Interestingly, the six classes of clones which were identified on the basis of the 3' clf sequences, the Archaic-Iberian, MW2, Pediatric-NY/Japan, EMRSA-15, EMRSA-16, and Berlin clones, are in agreement with the clonal complex classifications of the same strains defined by MLST: CC8, CC1, CC5, CC22, CC30, and CC45, respectively (23, 24). They also coincide with the sequence data for the fnb genes. Therefore, by sequencing of a single locus it was possible to identify not only the clonal complex corresponding to, for instance, the Pediatric and NY/Japan clones but also to distinguish between the NY/Japan strains isolated in the United States from one isolated in Japan. Owing to this dual capability, clf sequencing may have potential as a complementary typing method. Sequencing of the R domains of the clf genes needs to be performed with a larger number of isolates, and the in vitro and in vivo stabilities of these DNA regions need to be assayed in order to establish meaningful comparisons with other sequence typing data.
Some clfB sequences seem to register a genetic event already documented by spa typing and other sequencing data: the recombination of a large fragment of 557 kb from ST30 (EMRSA-16 background) into the branch of ST8 (Archaic-Iberian background) which gave rise to ST239 (Brazilian, Portuguese, and Hungarian lineage) (63). According to our data, the recombined fragment common to the EMRSA-16 clone and the Brazilian, Portuguese, and Hungarian cluster includes clfB, in addition to the spa and arcC genes, but did not include the fnb genes, which classified these clones into separate sequence types (Fig. 7). Further evidence for the relatedness between EMRSA-16 and the Brazilian, Portuguese, and Hungarian cluster in this stretch of the genome is the striking similarity between the clfB sequence in strain E1410, the MSSA strain from the 1960s which has the same MLST background as EMRSA-16, and the clfB sequences in strains belonging to the Brazilian, Portuguese, and Hungarian cluster. Yet, nucleotide mutations within the repeats allowed differentiation of strains with the EMRSA-16 clonal background from those belonging to the Brazilian, Portuguese, and Hungarian cluster. The importance of the background rate of nucleotide mutation within repeats, along with the extent of repeat number variation, has been documented for spa typing, which was recently been proposed as a means to address both long- and short-term epidemiological issues (37).
![]() View larger version (14K): [in a new window] |
FIG. 7. Map locations of virulence-related loci sequenced in this study (in boldface), as well as those of the spa and MLST genes. Sequence data for S. aureus strain N315 were downloaded from www.ncbi.nlm.nhi.gov/genomes.
|
Sequence variations in fnb genes. Sequence data for the fnb genes provided another source of information concerning the relationships between the S. aureus strains evaluated in this study. One fnbA sequence and one fnbB sequence corresponded to each class of clones defined by clfA sequencing, and these were the same for all clones within the class, thus confirming the classification obtained by clf sequencing and MLST. Most clones possessed two fnb genes, the only exceptions being clones EMRSA-16 and EMRSA-15, which had fnbB only. This is in agreement with previous studies in which not all but only between 77 and 91% of S. aureus clinical isolates possessed both fnbA and fnbB (57, 59). It was reported that isolates which differed in the number of fnb genes did not exhibit significant differences in fibronectin binding. On the other hand, fibronectin binding was negatively correlated with protease activity, which, in turn, is regulated by agr expression (59). This was documented in some detail for Canadian strain CMRSA-4, which had only one fnb gene (fnbB), was indistinguishable from EMRSA-16 by PFGE (70), and exhibited a high level of protease activity and a low-level capacity to bind to fibronectin.
Analysis of the fnb sequences in the strains that we tested resulted in the observation of a closer sequence similarity between the sequences of strain MW2 and those of the group comprising the Archaic, Iberian, Brazilian, Hungarian, and Portuguese clones (the Archaic-Iberian group) and complete identity between the sequences of this group and those of strain 8325-4. Previously reported motifs, such as the PIVP motif at the end of the fourth D repeat described for strain 8325-4 (36, 69), which is unique to the fnb sequences of this strain and the Archaic-Iberian group, or the GIDFVED motif described for Canadian epidemic CMRSA-1 (59), were also observed. In CMRSA-1 the GIDFVED motif has been reported to replace the SVDFEED epitope, which is essential for fibronectin binding, in one of the D repeats of FnbA as a possible strategy of immune evasion, but with a cost in terms of fibronectin binding. In our study, the GIDFVED motif replaced the SVDFEED motif in the fourth D repeat of the fnbA sequence for the Berlin clone, which, like CMRSA-1, possesses an additional D repeat in this gene. Are CMRSA-1 and the Berlin clone from the same lineage? Indeed, it appears to be so, although the Berlin clone was not included in the survey which compared Canadian MRSA isolates with international epidemic clones (70). However, the investigators referred to a previous report (7) in which CMRSA-1, under the designation of OE-MRSA, clustered with epidemic strains from Belgium, Switzerland, and Germany after randomly amplified polymorphic DNA analysis. In yet two other studies, the Canadian OE clone (8), also referred to as the Ontario epidemic clone (9), was shown by PFGE to belong to the same clonal type as the Berlin clone. The PIVP motif of fnb and the GIDFVED motif of fnbA may be useful as specific markers for clonal classification.
Sequence variations in some regions of agr. All strains in this study could be assigned to one of three major agr specificity groups. Most clones belonged to agr group I; but this reflected the fact that the Archaic, Iberian, Brazilian, Hungarian, and Portuguese clones (the Archaic-Iberian cluster of clones) have similar genetic backgrounds. If agr grouping is looked at as a genomic classification, the clones in this study were distributed more evenly: the Archaic-Iberian cluster, EMRSA-15, and Berlin clones in group I; the Pediatric and NY/Japan clones and strain N315 in group II; and the EMRSA-16 clone and strain MW2 in group III. Novick (49) has proposed that agr groups may represent ancient evolutionary divisions in terms of the organism's fundamental biology, and subsequent studies have linked the agr type to the genetic background of S. aureus from both disease (32, 33) and colonization (68) isolates. As observed previously (26), agr group I was the interference group which showed more genetic variation at the agr locus and which displayed the largest number of agrC sequence variants (which may again reflect the fact that more strains belonging to group I were analyzed). In particular, in one strain of the Hungarian clone, the tnp gene from IS256 was inserted into the 5' variable region of agrC. Instability in the agrC-coding region that leads to a truncated or a mutated protein has been documented (46, 72, 76). The mutations which occurred in AgrC during in vitro serial passage of S. aureus decreased the levels of production of secreted virulence factors and increased the growth yields of the bacteria, suggesting that the fitness of agr variants may be increased in certain ecological niches (72). On the other hand, the sequences of RNAIII and the region between the P2 and the P3 promoters were highly conserved in the set of clones analyzed in the present study. Yet, it would be interesting to clarify whether the point mutations in the P3 promoter in the Berlin and EMRSA-16 clones and in the P2 promoter in the Portuguese and EMRSA-15 clones have an effect on the overall regulation of agr in these clones.
Geographic dominance of clonal types and agr type. The form of bacterial interference mediated by inhibition of the synthesis of virulence factors and other extracellular proteins in S. aureus strains of different agr groups has been well documented in vitro (32, 34). Although different agr groups have been observed for several years in strains isolated from S. aureus carriers (76), only one agr type was detected at a single time in each healthy individual's nasal flora (27, 41). The observation that the clones with different agr types in our study correspond to distinct geographic areas in which these clones are dominant (i.e., are most frequently recovered) again raises the question of whether some type of agr-associated interference exists in vivo: strains of agr group I, represented by the Iberian, Brazilian, Portuguese, Hungarian, Berlin, and EMRSA-15 clones, are predominant in Europe and some South American countries; strains of group II, represented by the Pediatric and NY/Japan clones, have mainly been isolated in Japan and North America (but also in some European countries); and strains of group III, which were represented only by the EMRSA-16 clone, are also mainly isolated in Europe (see references 23 and 56 and the references therein). Although the agr types showed geographic overlap, some paradigmatic instances are noteworthy, such as the one in which two different agr types have coexisted for years in the same hospital but in different wards and in which one type does not overtake the other, as described for the Iberian (group I) and Pediatric (group II) clones in a Portuguese hospital (65). A similar case was recorded in Colombia, in which the Brazilian clone (group I) widely disseminated in other Southern American countries was totally absent and instead the dominant clone was the Pediatric clone of MRSA (group II) (28). On the other hand, the displacement of one major local clone with another was also observed in Portuguese hospitals: the Portuguese clone which was most frequently recovered during the 1985 surveillance study (18) was replaced by the Iberian clone (66) in 1992-1993, followed most recently by further replacement of the Iberian clone by the Brazilian clone (4). Each of these three clones belongs to the same agr type (group I). Similarly, clones belonging to agr group I were dominant in German hospitals in the 1990s (the Northern German [ST247] and Hannover [ST254] clones), and agr group I was still predominant in 2002, even though it was represented by different clones (the Berlin [ST45] and Barnim [ST22] clones). The rise of the Southern German (ST228) clone from agr group II observed in 2000 was not sustained in the following years, and it would be interesting to monitor the evolution of the recently emerged Rhine-Hesse (ST5) clone, also from group II (73). We hypothesize that, due to differences in genomic characteristics associated with a given agr type, MRSA epidemic clones belonging to three agr types may be competing for dominance in the hospital setting throughout the world.
In conclusion, the sequence polymorphisms observed in virulence-related loci may be associated with differential regulation by a global regulator of virulence genes (agr interpromoter region) or mechanisms that interfere with epidemiological dynamics (agr receptor variable region) in a collection of isolates representative of MRSA epidemic clones. Further experiments are warranted in order to investigate these issues. The sequences of the fnb genes also presented polymorphisms at the amino acid level in a region which is important for fibronectin binding. However, recent findings that the fibronectin-binding proteins have multiple, substituting fibronectin-binding regions (35, 43) suggest that the polymorphisms in the D region alone may not reflect a functional difference in fibronectin-binding capability and, consequently, may not have an effect on the capacity to initiate infection. Nevertheless, the fnb sequences together with the clf sequences have provided useful tools for genotypic characterization of MRSA isolates at a resolution higher than that provided by MLST. A specific motif in fnbA (GIDFVED) allowed the identification of a clonal type (Berlin) previously reported under several unrelated designations. In broad evolutionary terms, the types detected by fnb and clf sequencing were in agreement with those obtained by MLST and allowed the recognition of six lineages among the collection of MRSA isolates evaluated. In particular, the results of clfB sequencing, which had a discriminatory capacity greater than that of spa typing (Table 3), were also highly congruent with those of PFGE and spa typing, and if the last eight repeats (144 bp) of the R domain of this single locus are considered, the results of clfB sequencing are also congruent with those of MLST (Table 4). Thus, sequences from the R domain of clfB have strong potential for use in the typing of S. aureus strains.
|
View this table: [in a new window] |
TABLE 4. Cross-classification concordance levels for the collection of S. aureus isolates
|
This work was partially supported by Project POCTI/1999/ESP/34872, Fundação para a Ciência e a Tecnologia, Lisbon, Portugal. A. R. Gomes was supported by grants PRAXIS XXI/BPD/16375/98 and SFRH/BPD/9373/2002; and S. Vinga was supported by grant SFRH/BD/3134/2000, Fundação para a Ciência e a Tecnologia, Lisbon, Portugal.
|
|
|---|
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»