Previous Article | Next Article ![]()
Infection and Immunity, July 2007, p. 3673-3679, Vol. 75, No. 7
0019-9567/07/$08.00+0 doi:10.1128/IAI.01695-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

University Medical Center Utrecht, Heidelberglaan 100, 3584 CX Utrecht, The Netherlands,1 TNO Quality of Life, Business Unit Microbiology, Utrechtseweg 48, 3704HE Zeist, The Netherlands,2 National Institute of Public Health and the Environment (RIVM), Antonie van Leeuwenhoeklaan 9, 3721 MA Bilthoven, The Netherlands3
Received 24 October 2006/ Returned for modification 29 November 2006/ Accepted 28 March 2007
|
|
|---|
|
|
|---|
We compared the genetic compositions of M1 and M28 GAS strains from invasive and noninvasive diseases that were isolated long before and after the mid-1980s using a mixed-whole genome microarray to study whether genomic changes underlie the resurgence of invasive GAS disease in The Netherlands (5, 21). M1 and M28 are both predominant M types in The Netherlands. M1 is highly virulent and represents the predominant M type in The Netherlands (20); serotype M28 is particularly associated with puerperal sepsis. This is the first microarray study that compared genetic alterations in invasive and noninvasive GAS strains of different M types, spanning a period of four decades.
All M1 and M28 strains isolated after the 1980s were enriched in prophages encoding superantigens or streptodornases, irrespective of their source of isolation. As these phage-encoded virulence factors are associated with invasive disease (12, 18), these findings suggest that there was an overall increase in virulence among these strains. Thus, population-wide alterations in virulence rather than emergence of virulent subclones within given M types may have contributed to the resurgence of invasive GAS disease in The Netherlands.
|
|
|---|
Microarray construction. Random DNA fragments obtained from eight historical and more recent GAS strains were used to produce a mixed-genome DNA microarray (21). In short, random DNA fragments (1 to 1.9 kb) from these strains were cloned into Escherichia coli, PCR amplified, and spotted on the microarray. The probability P [0;1] to cover the entire GAS genome with an average length G (1.9 x 106 bp) and an average open reading frame length t (1,500 bp) can be estimated using Poisson distribution statistics with the following formula: P = 1 [1 (t + I 2r)G]N, where N is the number of significant biomarkers (2,704) (see Results) with average insert size I (1,450 bp). With an estimated minimal overlap size for successful hybridization between spotted biomarkers and labeled DNA of about 80 nucleotides (r), the probability for complete genome coverage was estimated to be 98%.
PCR-amplified probes specific for 34 known virulence genes (superantigens, adhesins, hemolysins/proteolytic enzymes, immunoreactive antigens, and regulatory elements) were also spotted on the microarray. Main genomic differences between the 41 GAS strains were identified by differentially hybridizing fragments on the microarray. DNA of differentiating biomarkers, defined as spots that had a positive Cy3 signal (reference DNA) for all isolates tested and no Cy5 signal (tester DNA) for at least one of the analyzed strains, was sequenced.
This microarray approach has several advantages. As the microarray is spotted with random DNA fragments obtained from relevant GAS strains, it does not require prior genome sequence information and allows identification of new genes. Furthermore, most GAS genomes published thus far are from North America, whereas the current approach provides a dedicated microarray to study well-documented GAS strains from The Netherlands.
Data analysis. Hierarchical clustering of differentiating biomarkers from all strains was done with TIGR software (14) (available at http://www.tigr.org/software/tm4). Clusters of biomarkers that emerged or disappeared in one of the M types over time were sequenced. To determine the possible function of a given sequence, ERGO bioinformatics was used (http://ergo.integratedgenomics.com/ERGO/), as were BLAST searches in the GenBank database. For the sake of simplicity we use the terms "acquisition" and "loss" to describe genetic differences between historical and recent strains.
Bacteriophage induction and lysis assay. GAS strains were grown overnight at 37°C in Todd-Hewitt medium plus yeast extract (THY). Each overnight GAS culture was diluted 1:100 with prewarmed THY and grown to an optical density at 660 nm of 0.2. Mitomycin C was added to the cultures to a final concentration of 0.2 µg/ml. Cultures were incubated for an additional 3 h at 37°C (2).
Mitomycin-treated bacteria were centrifuged at 4,000 x g for 15 min, and the supernatant was sterilized with a 0.22-µm-pore-size filter (Millipore). GAS target strains were grown to an optical density at 660 nm of 1.0, taken up in THY soft agar, and plated on THY agar plates. A phage mixture (10 µl) was spotted onto air-dried lawns. Plates were incubated overnight at 37°C, and lysis was defined as a visible clear area under the point where phage was applied. Experiments were performed three times.
Phage restriction analysis. Mitomycin-treated bacteria (see above) were centrifuged at 4,000 x g for 15 min. The supernatant was centrifuged at 141,000 x g for 4 h at 10°C, and the pellet was suspended in 400 µl of phage suspension buffer (0.15 M NaCl, 10 mM Tris HCl [pH 7.5], 5 mM MgCl2, 1 mM CaCl2). The phage particles were lysed with 0.5% sodium dodecyl sulfate, 10 mM EDTA, and 400 µg of proteinase K/ml for 1 h at 37°C. Phage DNA was extracted with phenol-chloroform-isoamyl alcohol (25:24:1), followed by chloroform-isoamyl alcohol (24:1), and precipitated with 0.2 volume of 3 M sodium acetate (pH 4.2) and 2 volumes of ethanol at 70°C overnight. Finally, DNA was washed with 70% ethanol and suspended in 40 µl distilled water (2). Digestion of the phage DNA was performed with restriction endonucleases EcoRI (Roche) and PinAI (Gibco BRL) according to the guidelines of the manufacturer. Restriction endonucleases were selected on the basis of their ability to distinguish between the different GAS prophages (8) using Webcutter 2.0 software (available at http://rna.lundberg.gu.se/cutter2/). Phage restriction fragments were analyzed on 1.0% agarose gels.
|
|
|---|
6180.1 (8),
6180.2 (8),
315.6 (3),
5005.1 (19), or
5005.3 (19) (Table 1). One cluster (cluster 4 in Fig. 1) matched both
315.6 and
5005.3 and consisted of biomarkers that span a 33-kb genetic fragment which is identical in these two prophages (1). All these phages carry at least one superantigen or streptodornase (Table 2).
![]() View larger version (38K): [in a new window] |
FIG. 1. Two-dimensional hierarchical clustering, illustrating differential biomarkers in recent (isolated after 1980) and old (isolated before 1980) M1 and M28 strains. The GAS strains are arranged along the x axis by M type and year of isolation (except isolates from the 1990s, which are [given the short time span] clustered in the biomarker profile). Isolates associated with invasive GAS disease are indicated by asterisks. Isolates associated with noninvasive GAS disease are indicated by dashes. The dendrogram on the y axis shows the clustering of 203 differentiating biomarkers. Red represents the presence of a biomarker, and green represents the absence of a biomarker. The bars at the top show the decade of isolation for each M type, as follows: pink indicates isolates obtained in the 1950s, brown indicates isolates obtained in the 1960s, yellow indicates isolates obtained in the 1970s, green indicates isolates obtained in the 1980s and blue indicates isolates obtained in the 1990s. The arrows indicate strains with the emm28.0'' allele. Braces indicate clusters of biomarkers, representing prophages. The numbering of the clusters corresponds to the numbering in Table 1.
|
|
View this table: [in a new window] |
TABLE 1. Differential biomarkers for old (isolated before 1980) and new (isolated after 1980) M1 and M28 strains and their putative functions
|
|
View this table: [in a new window] |
TABLE 2. Oligonucleotide primers used for detection of phage-specific virulence factors
|
5005.1 and
5005.3. In contrast, two strains from 1959 and 1960 lacked cluster
5005.3, whereas one of these strains also lacked cluster
5005.1. M28 strains showed a higher degree of variability. The 19 M28 isolates all belonged to emm28.0. However, within this cluster two different emm alleles were found. Fifteen of 19 isolates had sequences that were identical to the emm28.0 reference sequence (http://www.cdc.gov/ncidod/biotech/strep/emmdata.htm#emm28). The remaining four isolates showed an identical pattern consisting of seven base pair mutations in the emm sequence outside the emm typing region. These seven base pair mutations corresponded to three altered amino acids. These four M28 isolates were therefore designated emm28.0'' and are indicated in Fig. 1.
M28 strains from the 1990s possessed the
6180.2 cluster, the
315.6 cluster, or the
5005.3 and
5005.1 clusters. These clusters were all absent in strains isolated between the 1960s and the 1980s. All three GAS emm28.0'' isolates from the 1990s contained the
6180.1 cluster, in contrast to a GAS emm28.0'' isolate from the 1950s. Interestingly, the
6180.2 cluster was present in recent emm28.0 and emm28.0'' strains, whereas both of these emm variants lacked this cluster in the 1960s (Fig. 1).
To confirm that the different clusters containing prophage genes correspond to the prophages described in the GenBank database, additional prophage-specific PCRs, phage lytic assays, and phage restriction analyses were performed. For all 41 strains, the results of PCRs and sequencing of prophage-specific virulence factors (Table 2) were completely in line with the microarray-predicted prophage content (data not shown). All strains were challenged with mitomycin C to examine the phage release and phage sensitivity within the M1 and M28 strain collection using a bacterial lysis assay. Only from induced M1 and M28 strains that possessed the
5005.3 cluster were phages that demonstrated lytic activity collected. These phages lysed only M1 and M28 strains that did not contain either the
5005.3 or
315.6 cluster (data not shown). Restriction profiling of DNA isolated from mitomycin-induced phages essentially confirmed their identification based on microarray, PCR, and lysis experiments (Fig. 2).
![]() View larger version (39K): [in a new window] |
FIG. 2. EcoRI and PinAI DNA restriction profiles of prophages isolated from M1 and M28 strains. For each strain-prophage combination a representative restriction profile is shown. The observed restriction patterns corresponded to the predicted restriction patterns, except for a few fragments indicated with an asterisk. These aberrant fragments may be explained by mutations in the corresponding restriction sites. The numbers in parentheses identify isolates used for phage DNA restriction profiling in Fig. 1. The numbers next to the bands indicate the predicted DNA fragment sizes for the restriction profiles of prophage 5005.1, 5005.3, and 6108.2 sequences deposited in the GenBank database, as determined using Webcutter 2.0.
|
|
|
|---|
Previous findings in the United States showed that M1 isolates have acquired prophages 5005.1 and 5005.3 over time (19). We show that in The Netherlands, acquisition of phages 5005.1 and 5005.3 may already have occurred in the 1950s. The oldest M1 strain in our collection, isolated in 1959, appears to represent an "intermediate" M1 isolate that had already obtained
5005.1 but was still lacking
5005.3. Since all subsequent M1 isolates possessed both
5005.1 and
5005.3, these phages, or the M1 clone harboring them, may have rapidly gained predominance among all M1 isolates. The phage induction and lysis experiments suggest that
5005.3 remained highly mobile in all strains harboring this phage. The predominance of this phage among M1 strains, despite its mobility, suggests that it confers a selective advantage to the bacterium. Indeed,
5005.3 is likely to enhance the virulence of its host; Aziz and colleagues showed that phages 5005.1 (sphinx) and 5005.3 (phyramid) constituted the main difference between SF370, a sequenced strain that is infrequently associated with invasive infections, and a epidemic virulent M1 strain (1).
Prophages 5005.1 and 5005.3 were also present in an emm28.0 isolate from the 1990s. To our knowledge, this is the first report of the acquisition of M1 phages by another M type.
Analogous to M1, acquisition of these phages may have increased the virulence of M28 strains in The Netherlands. Given the early acquisition and predominance of prophages 5005.1 and 5005.3 in M1 isolates, these phages may have first been introduced into M1 and subsequently been transferred to M28. Indeed, lysis experiments indicated that phage 5005.3 could be induced from M1 strains and was able to infect and lyse M28 isolates that lacked this phage. Other emm28.0 isolates are likely to have obtained an M3 prophage (
315.6) containing the streptodornase-encoding sdn gene (2). Both microarray hybridization patterns and lysis assays suggest that prophages 315.6 and 5005.3 are mutually exclusive, probably since these phages are highly homologous and have the same phage repressor module. Phage 6180.2 may have been acquired by strains with both allelic variants of emm28.0 (emm28.0 and emm28.0'') on separate occasions (Fig. 3). Alternatively, phage 6180.2 may have been obtained by a single emm28.0 strain, which diverged into two emm28.0 variants after acquisition of the prophages. Since both emm28.0 subtypes (without this phage) were already present in 1960s isolates and all emm28.0'' strains had identical emm gene sequences, the latter explanation is unlikely. Phages 6180.1 and 6180.2 have recently been described as part of a completed M28 genome sequence (8). However, this is the first report describing the acquisition of these phages in M28 over time, based on the notion that the oldest M28 isolate in this study did not contain either one of these phages and
6180.1 was identified only in M28 strains isolated in the 1990s.
![]() View larger version (21K): [in a new window] |
FIG. 3. Hypothetical model of evolutionary events accommodating the data presented in this paper. M1 isolates may have acquired prophage 5005.1, encoding SpeA2, and 5005.3, encoding SdaD2, on separate occasions around 1960. A similar event may have taken place in emm28.0 in the 1990s. In these years, other emm28.0 isolates may have obtained an M3 prophage ( 315.6) containing the sdn gene. M28 isolates that did not acquire these M1 and M3 prophages may have obtained the M28 phages 6180.1 and 6180.2 with corresponding virulence factors. These events involved both types of the allelic variants of emm28.0. The presence of 6180.1 in historical emm28 isolates and its absence in recent emm28 isolates harboring prophage 315.6 or prophages 5005.1 and 5005.3 suggest that acquisition of the latter phages could be associated with loss of 6180.1 despite the fact that they do not have the same chromosomal insertion site (8) or repressor module (GenBank).
|
Taken together, our results suggest that over a period of four decades, M1 and especially M28 strains have acquired prophages on at least eight different occasions (Fig. 3). All the prophages carried superantigens (speA2, speC, speK), a phospholipase (sla), or streptodornases (sdaD2, sdn). These superantigens, streptodornases, and possibly the phospholipase are all associated with invasive GAS disease (12, 16, 18). Since prophage enrichment was similar in GAS strains associated with invasive disease and GAS strains associated with noninvasive disease, our results suggest that there was an overall increase in the virulence of M1 and M28 strains over the last four decades rather than emergence of hypervirulent subclones. This increased overall virulence potential may have enhanced the frequency of invasive GAS and hence contributed to the reemergence of invasive GAS disease in The Netherlands.
Published ahead of print on 23 April 2007. ![]()
|
|
|---|
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»