Previous Article | Next Article ![]()
Journal of Clinical Microbiology, September 2004, p. 3998-4006, Vol. 42, No. 9
0095-1137/04/$08.00+0 DOI: 10.1128/JCM.42.9.3998-4006.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
The Central Branch Office, Center for Disease Control, Taichung City, Taiwan
Received 11 February 2004/ Returned for modification 30 March 2004/ Accepted 9 May 2004
|
|
|---|
|
|
|---|
Analysis of the clonal relationships between clinical isolates from patients by various typing methods is a practical approach to elucidation of the epidemiology of a disease. To date, a number of phenotyping and genotyping methods have been described for S. pyogenes (1, 3, 10, 13, 18-20, 22, 26, 27, 30, 33). Among these methods, M serotyping has been taken as the "gold standard" for the characterization of S. pyogenes strains, in light of its importance to streptococcal virulence. However, application of M serotyping is restricted due to the lack of a comprehensive set of antisera in most laboratories and the high proportion of nontypeable isolates (21). In recent years, M serotyping has been replaced by sequencing of the 5' emm-coding region. Because the emm sequence can be used reliably to predict the M serotype and the sequences can be compared online, this method has been used worldwide as an epidemiologic tool to characterize S. pyogenes isolates recovered from patients with various diseases. This method has helped to identify many new emm and emm-like genes (12, 32). Vir typing and pulsed-field gel electrophoresis (PFGE) are the other two powerful subtyping methods for epidemiologic investigations of S. pyogenes infections. Vir typing, based on restriction length fragment polymorphisms of the 5- to 7-kb Vir regulon region of S. pyogenes, has been shown to be highly discriminatory and applicable to all S. pyogenes strains and provides a good correlation between the Vir type and the M serotype (13, 14). PFGE is based on restriction length fragment polymorphisms of the whole microbial genome and allows detection of variations among strains. Despite its usefulness, the results from Vir typing and PFGE generated between laboratories are not comparable unless a standardized protocol is followed.
In the study described here, 179 S. pyogenes isolates were characterized by the emm, Vir, and PFGE typing methods to investigate the epidemiologic aspects of scarlet fever in Taiwan. The protocols for the Vir and PFGE methods were standardized, and the typing patterns were analyzed with computer software that permits them to be compared online.
|
|
|---|
Prediction of suitable restriction enzymes for PFGE analysis of S. pyogenes. The Restriction Digest Tool provided on The Institute for Genome Research website (http://www.tigr.org/) was used to search for the restriction enzyme sites on the genome sequences of S. pyogenes strains MGAS315, MGAS8232, SF370, and SSI-1 to generate restriction fragment length profiles. Each of the four genomic sequences contained 8 to 57 restriction sites for the restriction enzymes, and the genomic sequences were further evaluated for their fragment size distributions. In the standard PFGE protocol, XbaI-digested genomic DNA fragments of S. enterica serovar Braenderup H9812 were used as reference size markers in PFGE gels. Therefore, only those fragments with sizes ranging from 21 to 1,135 kb were counted as analyzable by the pattern analysis software.
PFGE analysis. The standard PFGE protocol for S. pyogenes was developed on the basis of PulseNet's Listeria monocytogenes PFGE protocol, with minor modifications (15). Briefly, the S. pyogenes isolates were grown on blood agar plates incubated in 5% CO2 at 35°C for 16 to 24 h. The standard PFGE protocol for L. monocytogenes was followed to prepare the bacterial cell suspension, determine the bacterial concentration, make agarose plugs, lyse the bacterial cells, and wash the agarose plugs. Plug slices (width, 2 mm) were digested with 10 U of SmaI or 20 U of SgrAI. The DNA fragments were then separated in 1% Seakem Gold agarose gels (FMC BioProducts) at 14°C with a Bio-Rad CHEF Mapper apparatus (Bio-Rad Laboratories) in 0.5x Tris-borate-EDTA (TBE; pH 8) at a 120° fixed angle, at a fixed voltage (6 V/cm), and with pulse time intervals from 4 to 40 s for 20 h. The XbaI-digested genomic DNA fragments of S. enterica ser. Braenderup H9812 were used as reference size markers. The gel was stained with 1 mg of ethidium bromide per liter for 30 min and destained for 60 to 90 min with water that had been subjected to reverse osmosis, and the water was changed every 20 to 30 min. The gel was exposed on a UV transilluminator, and the image was captured digitally with a gel documentation system (AlphaImager 2000; Alpha Innotech Corporation, San Leandro, Calif.).
emm typing. The procedure developed by Beall et al. (1) was used to prepare the emm DNA fragments from the S. pyogenes isolates for nucleotide sequence determination. The amplified DNA amplicons and primer 1 were sent to a biotechnology company (Mission Biotech Corp., Taipei, Taiwan) for DNA sequencing. The first 160 bp of each of the 5' emm sequences was compared with the sequences in the emm database for determination of the degrees of homology(http://www.cdc.gov/ncidod/biotech/strep/strepindex.html) to determine emmtypes.
Vir typing. The standard Vir typing protocol was developed by modifying the procedures described by Gardiner et al. (13). Briefly, S. pyogenes was grown on 5% sheep blood agar and then incubated in 5% CO2 environment at 35°C for 16 to 24 h. A loop of bacterial growth was transferred into an Eppendorf tube containing 200 µl of lysozyme reaction solution (20 mM Tris-HCl [pH 8.0], 2 mM EDTA, 20 µg of lysozyme per ml) and then suspended by vigorously vortexing and incubated at 37°C for 30 min. The cells were then lysed by freezing and thawing the cell suspension five times in liquid nitrogen for 2 min and boiling water for 2 min. The lysate was then subjected to DNA extraction with a commercial kit (Blood & Tissue Genomic DNA Extraction Miniprep System; Viogene, Taipei, Taiwan). DNA was finally eluted with 100 µl of TE (Tris-EDTA) buffer and stored at 75°C until use.
The Vir regulon region was amplified with a commercial long PCR kit (TaKaRa Ex Taq; TaKaRa Shuzo Co., Kyoto, Japan). Briefly, a 50-µl PCR mixture containing 2 µl of the genomic DNA solution, 0.4 µM primers VUF and SBR (13), 250 µM deoxynucleoside triphosphates, 1x Ex Taq buffer, and 1.25 U of TaKaRa Ex Taq was subjected to 1 cycle at 94°C for 2 min and then 25 cycles at 94°C for 10 s, 60°C for 2 min, and 68°C for 6 min. Five microliters of the PCR product was digested with 5 U of HaeIII at 37°C for 1 h. The DNA fragments were then separated by electrophoresis in a 2% agarose gel in 0.5% TBE buffer at 6 V/cm for 2 h. The 100-bp ladder DNA (GeneTeks BioScience Inc., Taipei, Taiwan) was loaded at five-well intervals and used as a reference size marker. The gel was stained with 1 mg of ethidium bromide per liter for 30 min and destained for 60 min with water that had been subjected to reverse osmosis, and the water was changed every 20 to 30 min. The gel was exposed on a UV transilluminator, and the image was captured digitally with a gel documentation system (AlphaImager 2000, Alpha Innotech Corporation).
Computer analysis. The digital PFGE pattern and Vir pattern images were analyzed with BioNumerics software (Applied Maths, Kortrijk, Belgium). The computer software provided functions to set up a database of DNA patterns for pattern identification and phylogenetic analysis. With the aid of the computer software, dendrograms derived from the PFGE patterns and the Vir patterns were constructed by use of the unweighted pair group method with arithmetic averages (UPGMA) algorithm, based on Dice similarity coefficients.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Distribution of emm type by year
|
![]() View larger version (80K): [in a new window] |
FIG. 1. Representative Vir patterns of genomes of 16 S. pyogenes isolates digested with HaeIII. Vir typing was performed by the standard protocol described in Materials and Methods. Lanes M, 100-bp ladders used as reference size markers; lanes 1 to 16, S. pyogenes isolates Sp07233 (VT1, emm1), Sp15088 (VT1.2, emm1), Sp11029 (VT2, emm1), Sp10837 (VT3, emm4), Sp15066 (VT4, emm4), Sp15580 (VT4.2, emm4), Sp15581 (VT4.4, emm4), Sp15735 (VT14, emm6), Sp15459 (VT6, emm12), Sp18217 (VT6.2, emm12), Sp09459 (VT7, emm12), Sp04348 (VT8, emm12), Sp08850 (VT9, emm22), Sp11433 (VT10, emm33), Sp11792 (VT11, emm74), and Sp11014 (VT12, st11014), respectively.
|
![]() View larger version (39K): [in a new window] |
FIG. 2. Dendrogram and Vir patterns of S. pyogenes isolates and the association with emm types and the number of isolates belonging to the genotype. The dendrogram was constructed with BioNumerics software, with 3% optimization and 0.85% position tolerance, by using the UPGMA algorithm with Dice similarity coefficients.
|
![]() View larger version (21K): [in a new window] |
FIG. 3. Restriction profiles of 10 rarely cutting restriction enzymes for S. pyogenes strain MGAS315 and the number of restriction sites in the genomes of strains MGAS315, MGAS8232, SF370, and SSI-1. The prediction was done with the whole genome sequences of the indicated strains by use of the Restriction Digest Tool provided on The Institute for Genome Research website (http://www.tigr.org/). Fragments above the dashed line (>21 kb) are effective for analysis with the reference size markers, XbaI-digested DNA of S. enterica serovar Braenderup strain H9812. The numbers at the bottom indicate the numbers of restriction sites on the genome of the indicated strains.
|
![]() View larger version (87K): [in a new window] |
FIG. 4. Representative PFGE patterns of 11 S. pyogenes isolates obtained by the standard PFGE protocol with SmaI digestion described in Materials and Methods. Lanes M, chromosomal DNA of S. enterica serovar Braenderup H9812 digested with XbaI as reference size markers (size range, 20.5 to 1,135 kb). Lane 1, S. pyogenes isolate Sp15601 (SPYS16.0006, emm4); lane 2, Sp16314 (SPYS16.0026, emm4); lane 3, Sp15090 (SPYS16.0006, emm4); lane 4, Sp15581 (SPYS16.0006, emm4); lane 5, Sp15735 (SPYS16.0020, emm6); lane 6, Sp18297 (SPYS16.0021, emm6); lane 7, Sp04853 (SPYS16.0019, emm12); lane 8, Sp05265 (SPYS16.0019, emm12); lane 9, Sp7286 (SPYS16.0019, emm12); lane 10, Sp08504 (SPYS16.0013, emm12); lane 11, Sp8536 (SPYS16.0015, emm12). The chromosomal DNA of strain Sp16314 was resistant to SmaI digestion.
|
![]() View larger version (67K): [in a new window] |
FIG. 5. Dendrogram and PFGE patterns of SmaI-digested chromosomal DNA of S. pyogenes isolates and association with emm types, Vir patterns, and the number of isolates belonging to these genotypes. The dendrogram was constructed with BioNumerics software, with 4% optimization and 1% position tolerance, by using the UPGMA algorithm and Dice similarity coefficients. The clusters (clusters A, B, C, D, E, and F) contained isolates with similarity coefficients greater than 70% (indicated by the dashed line).
|
![]() View larger version (106K): [in a new window] |
FIG. 6. Representative PFGE patterns of genomes of 10 S. pyogenes isolates digested with SgrAI. The chromosomal DNA of the S. pyogenes isolates was resistant to SmaI digestion. Lanes M, chromosomal DNA of S. enterica serovar Braenderup H9812 digested with XbaI as reference size markers (size range, 20.5 to 1,135 kb); lane 1, S. pyogenes isolate Sp11433 (SPYS15.0006, emm33); lane 2, Sp15006 (SPYS15.0016, emm12); lane 3, Sp15066 (SPYS15.0001, emm4); lane 4, Sp15504 (SPYS15.0016, emm12); lane 5, Sp15685 (SPYS15.0014, emm12); lane 6, Sp15782 (SPYS15.0016, emm12); lane 7, Sp18278 (SPYS15.0016, emm12); lane 8, Sp18471 (SPYS15.0014, emm12); lanes 9 and 10, Sp16314 (SPYS15.0003, emm4).
|
![]() View larger version (87K): [in a new window] |
FIG. 7. Dendrogram and PFGE patterns of chromosomal DNA of S. pyogenes isolates digested with SgrAI and the association with SmaI PFGE patterns, emm types, and the numbers of isolates belonging to these genotypes. The dendrogram was constructed with BioNumerics software, with 3% optimization and 1% position tolerance, by using the UPGMA algorithm and Dice similarity coefficients. A 220-kb DNA fragment (shown by the arrow) appeared with the six strains (Sp15685, Sp18471, Sp15006, Sp15504, Sp15782, and Sp18278) resistant to digestion with SmaI.
|
|
View this table: [in a new window] |
TABLE 2. Distribution of major emm-Vir-PFGE genotypes by year
|
|
|
|---|
The M protein, which is encoded by the emm gene, is an important virulence factor for S. pyogenes (2). It has been reported that strains of certain M serotypes are epidemiologically associated with particular clinical syndromes. For example, both M1 and M3 strains were associated particularly with invasive diseases and fatal infections in Britain from 1980 to 1990 (8) and in the United States from 1995 to 1999 (23). M18 strains were associated with acute rheumatic fever (28). However, in a comparison study with control isolates from patients with noninvasive infections, Johnson et al. (17) indicated that the hypothesis that the increase in invasive infections in the United States during the late 1980s and 1990s was associated with particular M-protein types was not supported by the statistical evidence. Instead, the increase was a result of the increased numbers of strains of those particular M serotypes circulating in the general population at that time. Those investigators suggested that host factors, including individual and population-based immunity, must also be significant in influencing the infection potential. Among the M serotypes, M1, M2, M3, M4, M6, and M22 were detected among the isolates associated with scarlet fever (11, 24, 25, 34). In this study, we showed that the prevalent emm types for the S. pyogenes isolates collected in central Taiwan between 1996 and 1999 were emm4 (45%), emm12 (36%), emm1 (8%), and emm22 (7%). The epidemiologic pattern, however, is inconsistent with that in southern Taiwan, 180 km away. Yan et al. (34) showed that only the emm1 (72.7%), emm4 (26.0%), and emm25 (1.3%) types were found among the isolates recovered in southern Taiwan from 1993 to 2002. The emm1 isolates appeared in all years studied. In contrast, emm4 emerged only from 1995 to 2000, with a peak in 1997. Our data showed that in central Taiwan, emm12 is a major type among the isolates and the incidence rate was increasing in 1999. However, it was not identified in the isolates in southern Taiwan during the same period. Although the number of emm types identified in central and southern Taiwan was limited to a few types, it cannot be concluded that scarlet fever is associated only with those M serotypes. The high prevalence of those serotypes among the isolates causing scarlet fever during this period could result from an increase in the particular M serotypes circulating in the general population in that area.
Usually, only a small number of infected people develop symptoms in a scarlet fever outbreak (24). This also occurs with other invasive streptococcal diseases (6). As a streptococcal epidemic occurs, it is not easy to provide an outbreak alert by identification of only a few sporadic cases. Therefore, a sensitive surveillance system that uses modern molecular subtyping techniques is needed for the early detection of a streptococcal epidemic. PulseNet, a national molecular subtyping network for surveillance of food-borne diseases developed at the Centers for Disease Control and Prevention, is a successful example (31). This network has been proven to be effective for the early detection of food-borne disease outbreaks by PFGE analysis of the sporadic cases that occur in different geographical areas as a result of extensive food distribution (4). Molecular subtyping of S. pyogenes isolates by standard PFGE procedures and the development of a database of DNA fingerprints for these isolates will help to identify streptococcal epidemics at the early stages and allow monitoring of the changes in epidemiologic patterns.
PFGE patterns are generated in image form, but these are difficult to compare with those generated from different gels or those generated by different laboratories because the migration of the DNA bands in the gel can be affected by many factors, such as the strength of the agarose, the electrophoretic apparatus, and the electrophoretic conditions. To make PFGE image patterns comparable, the operating procedures that must be standardized include the DNA preparation methods, the restriction enzyme(s), the reference size markers, and the electrophoretic conditions used. In this study, we standardized a PFGE protocol for S. pyogenes by modifying the standard PFGE protocol for L. monocytogenes (15). Our results indicated that the DNA preparation and electrophoretic conditions used for Listeria work well for the S. pyogenes isolates and that the size range of the reference size markers made from S. enterica serovar Braenderup strain H9812 covered most of the DNA fragments generated by SmaI and SgrAI digestion of the genomic DNA. During the development of the PFGE standard protocol for S. pyogenes, the restriction enzymes were not selected from trial-and-error tests. Instead, they were searched for by using computer prediction of the whole genomic sequences of S. pyogenes strains and were chosen in view of enzyme cost and the size range of the reference size markers used. Among the 10 restriction enzymes found in the computer searches, RsrII, SanDI, SmaI, and SgrAI generated good PFGE patterns when the electrophoretic conditions and the reference size markers made from S. enterica serovar Braenderup strain H9812 were used. However, RsrII and SanDI are quite expensive and are not a good choice for use in routine PFGE analysis. Another five restriction enzymes (ApaI, BsePI, NaeI, PmeI, and SacII) generated many small fragments. They could be good for PFGE analysis under electrophoretic conditions that favor the separation of small fragments and in which reference size markers with a size range covering most of the fragments are used. Since many microbial genomic sequences have been deciphered, the restriction enzymes suitable for PFGE analysis of these microbes can be predicted by computer programs, which can save on the tremendous amounts of time and money usually spent on trial-and-error tests.
Our data showed that PFGE exhibited a higher discriminatory power than Vir typing for the local S. pyogenes isolates tested. However, the data showed that some isolates with identical PFGE patterns could be further discriminated by Vir typing (Fig. 5). Therefore, PFGE cannot completely replace Vir typing, which detects variations in the Vir regulon, a small region of 5 to 7 kb. Size variations in the small region cannot usually be detected by PFGE analysis. The dendrograms generated from the PFGE and Vir patterns showed that isolates with identical emm types are usually located in a distinct cluster (Fig. 2 and 5), implicating a common evolutionary origin for the emm isolates. The dendrograms generated from the SmaI PFGE patterns (Fig. 5) showed that emm12 isolates are located in two distinct clusters, indicating that the emm12 isolates could have been derived from two evolutionary origins or from a common origin that diverged into two groups at an earlier time. Compared to PFGE and Vir typing, emm typing exhibited much less discriminatory power. However, it was a good molecular marker that could be used to infer the evolutionary relationships among the local S. pyogenes isolates analyzed by PFGE and Vir typing.
The clustering analysis based on the SmaI PFGE patterns showed six distinct clusters. Each of the clusters contained isolates harboring the same emm type (Fig. 5). The different Vir-PFGE genotypes in a cluster could be derived from a common emm clone. A major emm-Vir-PFGE genotype usually existed in each cluster. For example, genotypes emm4-VT4-SPYS16.0006 and emm12-VT6-SPYS16.0013 were the major emm-Vir-PFGE genotypes in clusters B and C, respectively (Fig. 5). The clone of the major genotype usually emerged over several years and was distributed over a wide area. However, some emerged and then soon disappeared. For instance, emm1-VT1-SPYS16.0022 of cluster F, emm4-VT4-SPYS16.0006 of cluster B, emm12-VT6-SPYS16.0013 of cluster C, and emm12-VT6-SPYS16.0017 of cluster D were detected from 1997 to 1999; but emm4-VT3-SPYS16.0009 of cluster B, emm12-VT6-SPYS16.0019 of cluster D, and emm22-VT9-SPYS16.0001 of cluster A emerged over a short period of time (Table 2). A strain of emm12-VT6-SPYS16.0019 caused an outbreak in a kindergarten in December 1998 but did not spread into the community (5). Whether a major genotype can maintain its predominance or can be replaced by a genotype derived from it deserves to be watched. The database of fingerprints for the S. pyogenes isolates set up in this study will provide a basis of information for the long-term surveillance of the evolution of the genotypes.
PFGE identified nine isolates with DNA resistant to SmaI digestion. These isolates have emerged since 1998, and the number increased in 1999. Because DNA resistant to SmaI digestion was detected in isolates with different emm and Vir types, it is highly probable that this characteristic is carried on a mobile genetic element, which was supported by comparison of the PFGE patterns of the resistant and susceptible strains. By comparison of the SgrAI PFGE patterns, a 220-kb fragment was found in the SmaI digestion-resistant isolates but not in the SmaI digestion-susceptible emm12 isolates (Fig. 7), suggesting that this DNA had been inserted into an indigenous emm12-VT6-SPYS16.0013 clone that had been circulating in the area for a long time. However, there were no differences in the SgrAI PFGE patterns of the SmaI digestion-resistant and -susceptible emm4 isolates (Fig. 7). The genetic material was undetectable in the isolates by PFGE analysis. The SmaI digestion resistance phenomenon in an S. pyogenes strain was reported in 1997 by Cocuzza et al. (7). However, the genetic material responsible for this characteristic has not been studied yet. In addition to the modification on the SmaI site, whether the genetic material also carries a virulence factor(s) that leads to a dramatic increase in scarlet fever cases caused by the digestion-resistant strains in 2000 must be investigated further.
We have observed an increasing number of scarlet fever cases since 2000. In central Taiwan, the number of cases in 2002 was four times greater than that in 1999. This increase could result from the emergence of a virulent strain, a higher reporting rate, restrictions on antibiotic usage since 2000, or other environmental factors. The emm sequence, PFGE, and Vir fingerprint database for the S. pyogenes isolates made in this study will serve as a basis for information for epidemiologic studies of this disease.
We thank Bala Swaminathan and Lewis Graves of the Centers for Disease Control and Prevention for providing PFGE and BioNumerics training and Bernard Beall for identifying the new emm-like gene. We also thank our colleagues Shu-Ying Li and Hwa-Jen Teng for critical review of the manuscript.
|
|
|---|
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»