Evolutionary paths of streptococcal and staphylococcal superantigens
© Okumura et al.; licensee BioMed Central Ltd. 2012
Received: 3 April 2012
Accepted: 30 June 2012
Published: 17 August 2012
Streptococcus pyogenes (GAS) harbors several superantigens (SAgs) in the prophage region of its genome, although speG and smez are not located in this region. The diversity of SAgs is thought to arise during horizontal transfer, but their evolutionary pathways have not yet been determined. We recently completed sequencing the entire genome of S. dysgalactiae subsp. equisimilis (SDSE), the closest relative of GAS. Although speG is the only SAg gene of SDSE, speG was present in only 50% of clinical SDSE strains and smez in none. In this study, we analyzed the evolutionary paths of streptococcal and staphylococcal SAgs.
We compared the sequences of the 12–60 kb speG regions of nine SDSE strains, five speG+ and four speG – . We found that the synteny of this region was highly conserved, whether or not the speG gene was present. Synteny analyses based on genome-wide comparisons of GAS and SDSE indicated that speG is the direct descendant of a common ancestor of streptococcal SAgs, whereas smez was deleted from SDSE after SDSE and GAS split from a common ancestor. Cumulative nucleotide skew analysis of SDSE genomes suggested that speG was located outside segments of steeper slopes than the stable region in the genome, whereas the region flanking smez was unstable, as expected from the results of GAS. We also detected a previously undescribed staphylococcal SAg gene, selW, and a staphylococcal SAg -like gene, ssl, in the core genomes of all Staphylococcus aureus strains sequenced. Amino acid substitution analyses, based on dN/dS window analysis of the products encoded by speG, selW and ssl suggested that all three genes have been subjected to strong positive selection. Evolutionary analysis based on the Bayesian Markov chain Monte Carlo method showed that each clade included at least one direct descendant.
Our findings reveal a plausible model for the comprehensive evolutionary pathway of streptococcal and staphylococcal SAgs.
KeywordsSuperantigen Streptococcus Staphylococcus Fenome comparison Bayes MCMC
Bacterial superantigens (SAgs) have been shown to cause the massive activation of host T cells, strongly influencing immunological disorders. To date, nearly 50 bacterial SAgs and related molecules have been described, primarily from Gram-positive bacteria[1–3]. Streptococcus pyogenes (GAS) is one species of bacteria that harbors SAg genes. Analyses of the entire genomes of 13 GAS isolates have shown that each contains two to seven SAg genes (Additional file1), almost all located in the prophage regions of the genome. In contrast, genes encoding the SAgs speG and smez in GAS strains are not located on these mobile genetic elements, although some are surrounded by transposons. Thus, speG and smez in GAS may have been inherited from an ancestor by horizontal gene transfer. Although speJ in M1 GAS is not located on these mobile genetic elements, speJ is not conserved in the genome sequence of other GAS isolates, except for MGAS6180 (data not shown); in some strains, an SAg similar to speC is called speJ. We recently sequenced the entire genome of Streptococcus dysgalactiae subsp. equisimilis (SDSE) [DDBJ: AP010935], a bacterium that causes life-threatening infectious diseases, including sepsis and streptococcal toxic shock syndrome, similar to GAS[5–7]. Analyses of its sequence showed that SDSE is the closest relative of GAS sequenced to date, with 65% identity (Additional file2). Streptococcal bacteria other than GAS, such as S. dysgalactiae subsp. dysgalactiae and S. equi[9–11], have been reported to harbor more than one gene encoding proteins similar to SAgs. In contrast, targeted microarray analyses of 216 GAS virulence genes including SAgs in 58 SDSE strains isolated from human infections showed that the only SAg gene present in SDSE was speG, with about 50% of SDSE strains not harboring this gene[13–15].
Other representative bacterial SAgs and their related products have been identified in Staphylococcus aureus. At least 20 distinct staphylococcal SAgs have been described, including toxic shock syndrome toxin-1 (TSST-1), staphylococcal enterotoxins (SEs), and staphylococcal superantigen-like proteins (SSL), also called staphylococcal enterotoxin-like proteins (SEls)[1–3]. Almost all staphylococcal SAg genes are located in mobile genetic elements, such as prophages, transposons, plasmids, and pathogenicity islands (PIs). The distribution of these mobile elements among S. aureus isolates varies considerably. PIs that harbor the gene encoding TSST-1 can be excised and transduced with high efficiency by a staphylococcal phage.
In addition to these staphylococcal SAgs, recent studies have identified staphylococcal superantigen-like proteins (SSLs, also known as SETs), which have structural features similar to those of SAgs but do not possess SAg activity. All of the SSLs described to date are located in mobile genetic elements. Interestingly TSST-1, a functional SAg, shows higher sequence and structural similarity to SSL than to staphylococcal SAgs.
Structural analysis of SAgs has suggested that they evolved through the recombination of two smaller β-strand motifs, similar to the immunoglobulin binding motifs of streptococcal proteins G and L and the oligosaccharide/oligonucleotide binding family, such as the B subunits of AB(5) heat-labile enterotoxins, including cholera toxin, pertussis toxin, and verotoxin[19, 20]. However, the origin and evolutionary pathways of streptococcal and staphylococcal SAgs have not been well described.
To elucidate the origin of streptococcal SAgs based on de novo sequencing of SDSE strains and whole genome sequences, we have analyzed the synteny of the regions surrounding speG and smez in 13 GAS and 9 SDSE genomes. We also analyzed the genomic structures of all S. aureus strains for which whole genome data are available. We detected a previously undescribed gene that encodes a SEA-like protein (designated selW) and genes encoding SSL-like proteins, all of which are conserved in all S. aureus strains sequenced to date and are located in the core chromosome, not in any mobile elements. These findings, in addition to amino acid substitution analyses based on window analysis, cumulative TA-skew analysis and evolutionary analysis according to the Bayesian Markov chain Monte Carlo method, which allows the evolutionary path of SAg to be determined in chronological order, we were able to trace the origin and molecular evolution of streptococcal and staphylococcal SAgs.
Results and discussion
Comparison of sequences of speG regions in the GAS and SDSE genomes
To exclude the possibility that speG was acquired from a streptococcal phage, we compared the 50 kb sequences surrounding speG, a size sufficient to detect sequences derived from prophages. Synteny maps of the respective speG regions were essentially conserved in GAS strains and GGS_124 (Figure1A), except for MGAS10750, which did not harbor the speG sequence present in the corresponding speG regions of GAS and GGS_124. We found that the speG region of each GAS genome contains two to ten genes, which encode factors similar to mobile elements and phage-related genes, such as transposase, IS and co-activator of prophage gene expression. In contrast, these mobile elements could not be detected in the corresponding speG region of GGS_124 (Figure1A). The synteny of the regions surrounding the speG gene was highly conserved in eight GAS genomes (i.e. SF370, MGAS5005, MGAS2096, MGAS9429, MGAS10270, SSI-1, NZ131, and MGAS6180), each of which contains seven to eight transposase- and phage-related genes. These regions were 94% to 100% identical with each other. In the Manfredo genome, we found that IS1239, which is widely distributed in various isolates of GAS, had been inserted into the speG coding sequence, resulting in speG being a pseudogene in this strain (Figure1A). In the MGAS8232 genome, IS1239 flanked speG.
Although a previous study suggested that speG transferred from SDSE to GAS, our results clearly indicate that the synteny surrounding speG in the GAS and GGS_124 genomes has been essentially conserved and that modifications of this context, by insertion of mobile elements, occurred only in GAS strains. These results strongly suggest that speG in GAS and SDSE is an orthologous, not a xenologous, gene, the latter defined as a gene displaced by horizontal transfer from another lineage. Moreover, speG in GAS and SDSE is a descendant of an ancestral streptococcal SAg and has been conserved in evolution.
Comparison of the speG region among SDSE strains
Although the synteny of speG regions of GAS strains and GGS_124 has been highly conserved, about 50% of SDSE strains do not harbor speG[13–15]. We therefore selected nine SDSE isolates, five with (GGS_124, 163, 164, 168 and 170) and four without (RE378, SDSE_118, 160, and 165) speG (Additional file1). Following direct genome sequencing (GGS_124 and RE378) or PCR amplification using speG specific primers (Additional file3), we compared the sequences of these nine strains. Each of these isolates harbored a different emm type (Additional file1), widely used to type GAS and SDSE strains. We also included the full sequence of the ATCC 12394 genome [GenBank: CP002215], an SDSE that does not harbor speG.
When we analyzed the genetic structures surrounding speG (12 to 60 kb) in these SDSE isolates (Figure1B), we found that, in general, these structures were highly conserved, especially in the 12 kb regions between pgi (pink) and perR (blue), but that speG itself and its corresponding regions were not. Outside these 12 kb regions, we found that most of these strains contained 1 or 2 coding sequences similar to transposase or IS elements, including several that appeared to be common to the sequenced GAS genomes.
Remarkably, all four speG-negative strains (e.g. RE378, SDSE_118, 160, and 165) showed the insertion of an approximately 20 kb fragment between the hypothetical protein gene (locus_tag: SDEG_1990) and the gene similar to peptidoglycan endo-beta-N-acetylglucosaminidase (locus_tag: SDEG_1992) present in the GGS_124 genome, replacing speG (locus_tag: SDEG_1991) at the exact same site. These 20 kb fragments were composed of 19 or 22 coding sequences, which were similar to genes derived from evolutionally distant species such as Clostridium botulinum and C. tetani. However, the arrangements of these genes did not exactly match those of the clostridial genomes (data not shown), with most coding sequences sharing <60% similarity (e.g. Additional file4). In contrast, genetic structures other than these 20 kb fragments were highly conserved among the speG-positive and -negative strains (Figure1B). These findings indicated that synteny had been conserved in the regions surrounding speG, or the inserted 20 kb fragments, of these SDSE strains.
Analysis of the expression of speG region genes
Analysis of the smez region
We also analyzed the genome context surrounding smez, another chromosomally encoded SAg, in GAS and SDSE strains. We found that the synteny of approximately 20 kb regions containing smez (see details for MGAS10270 and NZ131, below) were highly conserved in all GAS strains sequenced (Figure2). Of the 13 completely sequenced GAS strains, 11 harbored smez genes, primarily at approximately 1.7 Mb. Although smez was present at the same site in SSI-1 as in the other GAS genomes, the former was functionally inactive due to a frame-shift mutation. In contrast, the MGAS10270 and NZ131 genomes did not contain smez fragments, even at other locations, despite their corresponding surrounding genome structures being highly conserved when compared with the other GAS strains. All GAS genomes contained highly similar dpp operons (dppA, dppB, dppC, dppD, and dppE) immediately downstream of smez, and all contained flaR and trpG, located upstream of smez coding sequences (Additional file5).
Analysis of the GGS_124, RE378, and ATCC 12394 strains revealed that none contained fragments similar to the 702 bp smez coding sequence derived from the SF370 genome. We therefore searched for flaR and the dpp operon, which were highly conserved in the smez flanking regions of GAS genomes (Figure2). In these three SDSE genomes, flaR was located at about 0.2 Mb, whereas the dpp operon was located at about 0.9 Mb, far from the position of flaR (Figure2 and Additional file5). Furthermore, synteny of the regions surrounding flaR and the dpp operon was not well conserved in these three SDSE genomes, suggesting rearrangement of the genome context. The flaR gene and the dpp operon show high similarities in GAS and SDSE (Additional file5), with concomitant sequences observed only in GAS and SDSE but not in other streptococci (data not shown).
We next plotted cumulative TA-skew diagrams of the three sequenced SDSE chromosomes (GGS_124, RE378, and ATCC 12394). Use of a similar method on 12 sequenced GAS genomes showed that all cumulative TA-skew curves of GAS genomes displayed a V-shape, interrupted by segments of steeper slopes, called steep-slope regions (SSRs). Diagram distortions including SSRs are thought to correspond to positions in which foreign genetic elements are integrated, including prophage-related genes, horizontally acquired elements, and pathogenicity islands, and in which genome rearrangements occur. The SSR was conserved among GAS strains, with smez at the border of the SSR, suggesting that this region is predisposed to be unstable.
Identification of conserved enterotoxin like and staphylococcal superantigen like genes in all S. aureus genomes
S. aureus also harbors many types of SAgs, such as TSST-1, SEs, and SSLs. We identified a relatively unknown staphylococcal SAg, selW, and an ssl gene cluster, both of which are conserved in all S. aureus genomes examined to date. Moreover, we found that each of these genes was located in the same chromosomal region of the S. aureus genomes, not within any mobile elements. The highly syntenic conservation of selW an d the ssl gene cluster among S. aureus genomes and their similarity to SEs and SSLs, respectively, suggest that they are likely the direct descendants of common ancestral SEs and SSLs, respectively.
Evolutionary analysis of streptococcal and staphylococcal SAgs and SSL
The physiological activities and three-dimensional structures of SAgs are quite similar in streptococci and staphylococci. Although many studies have focused on staphylococcal SAgs in mobile elements, little is known about staphylococcal SAg-related gene(s) located on the core chromosome. To analyze the relationship among SAgs, we employed the Bayesian MCMC method. Although the phylogenetic tree we obtained was similar to that observed previously report, the method we used makes possible the determination of the temporal evolution of SAgs. Evolutionary analysis of the streptococcal and staphylococcal SAgs, and their related products, SSLs, showed that those molecules could be divided into three clades, each of which contains at least one direct descendant of an ancestor. SAgs of clades I and III consist of streptococcal and staphylococcal SAgs, respectively. In contrast, clade III is a mixture of streptococcal and staphylococcal SAgs, containing only SELW of S. aureus.
Streptococcal SAgs are one of the important virulence factors involved in life-threatening diseases such as streptococcal toxic shock syndrome (STSS) and scarlet fever. At present, a total of 11 SAgs have been identified by GAS genome sequencing, with most GAS isolates possessing several SAg genes in their genomes. Although the diversity of SAgs is thought to arise during horizontal transfer, their evolutionary pathway has not been determined. To better understand SAg evolution, we sequenced the entire genome of SDSE, the closest relative of GAS, which harbors speG as its only SAg gene. Genome-wide comparisons of GAS and SDSE provided evidence that speG is the direct descendant of a common ancestor of the streptococcal SAg. Furthermore, we also detected previously undescribed inter-species horizontal SAg gene transfer events among three pathogens, S. pyogenes, S. dysgalactiae subsp. equisimili s and S. aureus. This study is the first time to describe the origin and evolution of SAgs in pathogenic streptococci and staphylococci. These findings suggest that horizontal gene transfer is a more ubiquitous genetic exchange system than previously known, and that it sometimes crosses interspecies barriers.
Bacterial strains and media
All S. dysgalactiae subsp. equisimilis (SDSE) strains used in this study were isolated from patients with invasive infections in different hospitals throughout Japan (Additional file1). Each SDSE isolate was cultured in 5% sheep blood agar or Brain Heart Infusion medium at 37°C under 5% CO2 as described.
Preparation of genomic DNA and sequencing
Streptococci were lysed as described, and genomic DNA was purified using Wizard® Genomic DNA Purification Kits (Promega).
PCR reactions were performed in volumes of 50 μl containing TaKaRa ExTaq DNA polymerase (TaKaRa), with amplification on a GeneAmp PCR System 9700 (Applied Biosystems). Primer sets for direct sequencing were based on GGS_124 and RE378 genome sequence data, with each set designed to amplify 5 kbp PCR products with 500 bp overlapping regions. The PCR primer set for the speG-specific region has been described (see also Additional file3). PCR products were electrophoresed on 1.0% agarose gels and purified using QIAquick PCR Purification kits (QIAGEN). All DNA fragments were sequenced on an ABI3100 DNA sequencer with a redundancy of 4.
All isolates were grown overnight in 10 ml of Brain Heart Infusion medium at 37°C under 5% CO2 in 15 ml conical tubes. The cells were harvested, total RNA was purified using RNeasy Mini Kits (QIAGEN), and RNA concentrations were measured using a NanoDrop™ 1000 spectrophotometer (Thermo Scientific).
RT-PCR analysis of target genes
Total RNA was reversed transcribed into cDNA using Superscript III reverse transcriptase kits (Invitrogen) and oligo dT primers. PCR amplifications were performed using the primer sequences in Additional file3. The PCR products were electrophoresed on 1.0% agarose gels and detected by UV-fluorescence after ethidium bromide staining.
Bioinformatics and evolutionary analyses
Homology searches and IS searches were performed using BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi) and IS finder (http://www-is.biotoul.fr/is.html), respectively. Cumulative TA skew analysis was performed using GenSkew (http://genskew.csb.univie.ac.at/). Both window size and stepsize sequence length were set at 1000 bp. Rates of evolution were estimated by a Window Analysis of dN and dS, using the online interface of WINA 0.34, in a sliding window size of 60 bp (20 codons) at 3 bp intervals. A phylogenetic tree was constructed with CLUSTALW (http://clustalw.ddbj.nig.ac.jp/top-j.html), MrBayes 3.1.2 (http://mrbayes.csit.fsu.edu/index.php) and TreeView X (http://darwin.zoology.gla.ac.uk/~rpage/treeviewx/index.html) software. Sequences were manually corrected using GENETYX-Mac (GENETYX Co.) and gene analysis was performed by in silico molecular cloning (in silico Biology Co.).
Nucleotide and amino acid sequence accession number
The DNA sequences of the region surrounding speG (10–60 kb) in each strain have been deposited in the DDBJ under the accession numbers listed in Additional file1. Accession numbers for SAgs used in this study are listed in Additional file9.
The authors thank Dr. Omoe at Iwate University for valuable discussions about the evolution of SEs. This work was partly supported by a grant for Research on Emerging and Reemerging Infectious Diseases (H22 Shinkouh-013). T. M. A was supported by JSPS KAKENHI Grant Number 21590503 and 24390109.
- Alouf JE, Müller H: What are superantigens?. The comprehensive sourcebook of bacterial protein toxins. Edited by: Alouf JE, Popoff MR. 2006, New York: Academic Press, 821-829. 3View Article
- Fraser JD, Proft T: The bacterial superantigen and superantigen-like proteins. Immunol Rev. 2008, 225: 226-243. 10.1111/j.1600-065X.2008.00681.x.View ArticlePubMed
- Brosnahan AJ, Schlievert PM: Gram-positive bacterial superantigen outside-in signaling causes toxic shock syndrome. FEBS J. 2011, 278 (23): 4649-4667. 10.1111/j.1742-4658.2011.08151.x.PubMed CentralView ArticlePubMed
- Shimomura Y, Okumura K, Murayama SY, Yagi J, Ubukata K, Kirikae T, Miyoshi-Akiyama T: Complete genome sequencing and analysis of a Lancefield group G Streptococcus dysgalactiae subsp. equisimilis strain causing streptococcal toxic shock syndrome (STSS). BMC Genomics. 2011, 12 (1): 17-10.1186/1471-2164-12-17.PubMed CentralView ArticlePubMed
- McDonald M, Towers RJ, Andrews RM, Carapetis JR, Currie BJ: Epidemiology of Streptococcus dysgalactiae subsp. equisimilis in tropical communities, Northern Australia. Emerg Infect Dis. 2007, 13 (11): 1694-1700. 10.3201/eid1311.061258.PubMed CentralView ArticlePubMed
- Tanaka D, Isobe J, Watahiki M, Nagai Y, Katsukawa C, Kawahara R, Endoh M, Okuno R, Kumagai N, Matsumoto M, Working Group for Group A Streptococci in Japan, et al: Genetic features of clinical isolates of Streptococcus dysgalactiae subsp. equisimilis possessing Lancefield's group A antigen. J Clin Microbiol. 2008, 46 (4): 1526-1529. 10.1128/JCM.02188-07.PubMed CentralView ArticlePubMed
- Broyles LN, Van Beneden C, Beall B, Facklam R, Shewmaker PL, Malpiedi P, Daily P, Reingold A, Farley MM: Population-based study of invasive disease due to beta-hemolytic streptococci of groups other than A and B. Clin Infect Dis. 2009, 48 (6): 706-712. 10.1086/597035.View ArticlePubMed
- Miyoshi-Akiyama T, Zhao J, Kato H, Kikuchi K, Totsuka K, Kataoka Y, Katsumi M, Uchiyama T: Streptococcus dysgalactiae-derived mitogen (SDM), a novel bacterial superantigen: characterization of its biological activity and predicted tertiary structure. Mol Microbiol. 2003, 47 (6): 1589-1599. 10.1046/j.1365-2958.2003.03411.x.View ArticlePubMed
- Proft T, Webb PD, Handley V, Fraser JD: Two novel superantigens found in both group A and group C Streptococcus. Infect Immun. 2003, 71 (3): 1361-1369. 10.1128/IAI.71.3.1361-1369.2003.PubMed CentralView ArticlePubMed
- Thomas P, Webb PD, Handley V, Fraser JD: Identification & characterisation of the two novel streptococcal pyrogenic exotoxins SPE-L & SPE-M. Indian J Med Res. 2004, 119 (Suppl): 37-43.PubMed
- Alber J, El-Sayed A, Estoepangestie S, Lammler C, Zschock M: Dissemination of the superantigen encoding genes seeL, seeM, szeL and szeM in Streptococcus equi subsp. equi and Streptococcus equi subsp. zooepidemicus. Vet Microbiol. 2005, 109 (1–2): 135-141.View ArticlePubMed
- Igwe EI, Shewmaker PL, Facklam RR, Farley MM, van Beneden C, Beall B: Identification of superantigen genes speM, ssa, and smeZ in invasive strains of beta-hemolytic group C and G streptococci recovered from humans. FEMS Microbiol Lett. 2003, 229 (2): 259-264. 10.1016/S0378-1097(03)00842-5.View ArticlePubMed
- Davies MR, McMillan DJ, Beiko RG, Barroso V, Geffers R, Sriprakash KS, Chhatwal GS: Virulence profiling of Streptococcus dysgalactiae subspecies equisimilis isolated from infected humans reveals 2 distinct genetic lineages that do not segregate with their phenotypes or propensity to cause diseases. Clin Infect Dis. 2007, 44 (11): 1442-1454. 10.1086/516780.View ArticlePubMed
- Hashikawa S, Iinuma Y, Furushita M, Ohkura T, Nada T, Torii K, Hasegawa T, Ohta M: Characterization of group C and G streptococcal strains that cause streptococcal toxic shock syndrome. J Clin Microbiol. 2004, 42 (1): 186-192. 10.1128/JCM.42.1.186-192.2004.PubMed CentralView ArticlePubMed
- Brandt CM, Schweizer KG, Holland R, Lutticken R, Freyaldenhoven BS: Lack of mitogenic activity ofspeG- andspeGdys-positiveStreptococcus dysgalactiaesubspeciesequisimilisisolates from patients with invasive infections. Int J Med Microbiol. 2005, 295 (8): 539-546. 10.1016/j.ijmm.2005.07.013.View ArticlePubMed
- Baba T, Takeuchi F, Kuroda M, Yuzawa H, Aoki K, Oguchi A, Nagai Y, Iwama N, Asano K, Naimi T, et al: Genome and virulence determinants of high virulence community-acquired MRSA. Lancet. 2002, 359 (9320): 1819-1827. 10.1016/S0140-6736(02)08713-5.View ArticlePubMed
- Ruzin A, Lindsay J, Novick RP: Molecular genetics of SaPI1–a mobile pathogenicity island in Staphylococcus aureus. Mol Microbiol. 2001, 41 (2): 365-377. 10.1046/j.1365-2958.2001.02488.x.View ArticlePubMed
- Williams RJ, Ward JM, Henderson B, Poole S, O'Hara BP, Wilson M, Nair SP: Identification of a novel gene cluster encoding staphylococcal exotoxin-like proteins: characterization of the prototypic gene and its protein product, SET1. Infect Immun. 2000, 68 (8): 4407-4415. 10.1128/IAI.68.8.4407-4415.2000.PubMed CentralView ArticlePubMed
- Arcus VL, Proft T, Sigrell JA, Baker HM, Fraser JD, Baker EN: Conservation and variation in superantigen structure and activity highlighted by the three-dimensional structures of two new superantigens from Streptococcus pyogenes. J Mol Biol. 2000, 299 (1): 157-168. 10.1006/jmbi.2000.3725.View ArticlePubMed
- Mitchell DT, Levitt DG, Schlievert PM, Ohlendorf DH: Structural evidence for the evolution of pyrogenic toxin superantigens. J Mol Evol. 2000, 51 (6): 520-531.View ArticlePubMed
- Lina G, Bohach GA, Nair SP, Hiramatsu K, Jouvin-Marche E, Mariuzza R, International Nomenclature Committee for Staphylococcal Superantigens: Standard nomenclature for the superantigens expressed by Staphylococcus. J Infect Dis. 2004, 189 (12): 2334-2336. 10.1086/420852.View ArticlePubMed
- Kapur V, Reda KB, Li LL, Ho LJ, Rich RR, Musser JM: Characterization and distribution of insertion sequence IS1239 in Streptococcus pyogenes. Gene. 1994, 150 (1): 135-140. 10.1016/0378-1119(94)90872-9.View ArticlePubMed
- Sachse S, Seidel P, Gerlach D, Gunther E, Rodel J, Straube E, Schmidt KH: Superantigen-like gene(s) in human pathogenicStreptococcus dysgalactiaesubsp.equisimilis: genomic localisation of the gene encoding streptococcal pyrogenic exotoxin G (speGdys). FEMS Immunol Med Microbiol. 2002, 34 (2): 159-167. 10.1111/j.1574-695X.2002.tb00618.x.View ArticlePubMed
- Koonin EV, Makarova KS, Aravind L: Horizontal gene transfer in prokaryotes: quantification and classification. Annu Rev Microbiol. 2001, 55: 709-742. 10.1146/annurev.micro.55.1.709.View ArticlePubMed
- Zhao J, Hayashi T, Saarinen S, Papageorgiou AC, Kato H, Imanishi K, Kirikae T, Abe R, Uchiyama T, Miyoshi-Akiyama T: Cloning, expression, and characterization of the superantigen streptococcal pyrogenic exotoxin G from Streptococcus dysgalactiae. Infect Immun. 2007, 75 (4): 1721-1729. 10.1128/IAI.01183-06.PubMed CentralView ArticlePubMed
- Steer AC, Law I, Matatolu L, Beall BW, Carapetis JR: Global emm type distribution of group A streptococci: systematic review and implications for vaccine development. Lancet Infect Dis. 2009, 9 (10): 611-616. 10.1016/S1473-3099(09)70178-1.View ArticlePubMed
- Suzuki H, Lefebure T, Hubisz MJ, Pavinski Bitar P, Lang P, Siepel A, Stanhope MJ: Comparative genomic analysis of the Streptococcus dysgalactiae species group: gene content, molecular adaptation, and promoter evolution. Genome Biol Evol. 2011, 3: 168-185. 10.1093/gbe/evr006.PubMed CentralView ArticlePubMed
- Mulligan ME, Hawley DK, Entriken R, McClure WR: Escherichia coli promoter sequences predict in vitro RNA polymerase selectivity. Nucleic Acids Res. 1984, 12 (1 Pt 2): 789-800.PubMed CentralView ArticlePubMed
- Panchaud A, Guy L, Collyn F, Haenni M, Nakata M, Podbielski A, Moreillon P, Roten CA: M-protein and other intrinsic virulence factors of Streptococcus pyogenes are encoded on an ancient pathogenicity island. BMC Genomics. 2009, 10: 198-10.1186/1471-2164-10-198.PubMed CentralView ArticlePubMed
- Roten CA, Gamba P, Barblan JL, Karamata D: Comparative Genometrics (CG): a database dedicated to biometric comparisons of whole genomes. Nucleic Acids Res. 2002, 30 (1): 142-144. 10.1093/nar/30.1.142.PubMed CentralView ArticlePubMed
- Grigoriev A: Analyzing genomes with cumulative skew diagrams. Nucleic Acids Res. 1998, 26 (10): 2286-2290. 10.1093/nar/26.10.2286.PubMed CentralView ArticlePubMed
- Highlander SK, Hulten KG, Qin X, Jiang H, Yerrapragada S, Mason EO, Shang Y, Williams TM, Fortunov RM, Liu Y, et al: Subtle genetic changes enhance virulence of methicillin resistant and sensitive Staphylococcus aureus. BMC Microbiol. 2007, 7: 99-10.1186/1471-2180-7-99.PubMed CentralView ArticlePubMed
- Rosenstein R, Nerz C, Biswas L, Resch A, Raddatz G, Schuster SC, Gotz F: Genome analysis of the meat starter culture bacterium Staphylococcus carnosus TM300. Appl Environ Microbiol. 2009, 75 (3): 811-822. 10.1128/AEM.01982-08.PubMed CentralView ArticlePubMed
- Davies MR, McMillan DJ, Van Domselaar GH, Jones MK, Sriprakash KS: Phage 3396 from a Streptococcus dysgalactiae subsp. equisimilis pathovar may have its origins in Streptococcus pyogenes. J Bacteriol. 2007, 189 (7): 2646-2652. 10.1128/JB.01590-06.PubMed CentralView ArticlePubMed
- Vojtek I, Pirzada ZA, Henriques-Normark B, Mastny M, Janapatla RP, Charpentier E: Lysogenic transfer of group A Streptococcus superantigen gene among Streptococci. J Infect Dis. 2008, 197 (2): 225-234. 10.1086/524687.PubMed CentralView ArticlePubMed
- Holden MT, Heather Z, Paillot R, Steward KF, Webb K, Ainslie F, Jourdan T, Bason NC, Holroyd NE, Mungall K: Genomic evidence for the evolution of Streptococcus equi: host restriction, increased virulence, and genetic exchange with human pathogens. PLoS Pathog. 2009, 5 (3): e1000346-10.1371/journal.ppat.1000346.PubMed CentralView ArticlePubMed
- Aziz RK, Edwards RA, Taylor WW, Low DE, McGeer A, Kotb M: Mosaic prophages with horizontally acquired genes account for the emergence and diversification of the globally disseminated M1T1 clone of Streptococcus pyogenes. J Bacteriol. 2005, 187 (10): 3311-3318. 10.1128/JB.187.10.3311-3318.2005.PubMed CentralView ArticlePubMed
- Synnott AJ, Kuang Y, Kurimoto M, Yamamichi K, Iwano H, Tanji Y: Isolation from sewage influent and characterization of novel Staphylococcus aureus bacteriophages with wide host ranges and potent lytic capabilities. Appl Environ Microbiol. 2009, 75 (13): 4483-4490. 10.1128/AEM.02641-08.PubMed CentralView ArticlePubMed
- Chen J, Novick RP: Phage-mediated intergeneric transfer of toxin genes. Science. 2009, 323 (5910): 139-141. 10.1126/science.1164783.View ArticlePubMed
- Endo T, Ikeo K, Gojobori T: Large-scale search for genes on which positive selection may operate. Mol Biol Evol. 1996, 13 (5): 685-690. 10.1093/oxfordjournals.molbev.a025629.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.