Phylogenetic relationship of prophages is affected by CRISPR selection in Group A Streptococcus

BACKGROUND: Group A Streptococcus (GAS) is a major human pathogen, which is associated with a wide spectrum of invasive diseases, such as pharyngitis, scarlet fever, rheumatic fever, and streptococcal toxic shock syndrome (STSS). It is hypothesized that differences in GAS pathogenicity are related to the acquisition of diverse bacteriophages (phages). Nevertheless, the GAS genome also harbors clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (cas) genes, which play an important role in eliminating foreign DNA, including those of phages. However, the structure of prophages in GAS strains is mosaic, and the phylogenetic relationship between prophages and CRISPR is not clear. In this study, we analyzed CRISPR and prophage structure using 118 complete genome sequences of GAS strains to elucidate the relationship between two genomic elements. Additionally, phylogenetic and M-type analyses were performed. RESULTS: Of the 118 GAS strains, 80 harbored type I-C and/or II-A CRISPR/cas loci. A total of 553 spacer sequences were identified from CRISPR/cas loci and sorted into 229 patterns. We identified and classified 373 prophages into 14 groups. Some prophage groups shared a common integration site, and were related to M-type. We further investigated the correlation between spacer sequences and prophages. Of the 229 spacer sequence patterns, 203 were similar to that of other GAS prophages. No spacer showed similarity with that of a specific prophage group with mutL integration site. Moreover, the average number of prophages in strains with type II-A CRISPR was significantly less than that in type I-C CRISPR and non-CRISPR strains. However, there was no statistical difference between the average number of prophages in type I-C strains and that in non-CRISPR strains. CONCLUSIONS: Our results indicated that type II-A CRISPR may play an important role in eliminating phages and that the prophage integration site may be an important criterion for the acceptance of foreign DNA by GAS. M type, spacer sequence, and prophage group data were correlated with the phylogenetic relationships of GAS. Therefore, we hypothesize that genetic characteristics and/or phylogenetic relationships of GAS may be estimated by analyzing its spacer sequences. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1186/s12866-019-1393-y) contains supplementary material, which is available to authorized users.