Figures
Abstract
Background
The surveillance of emm types and macrolide susceptibility of group A streptococcus (GAS) in various areas and time periods enhances the understanding of the epidemiology of GAS infections and may guide treatment strategies and the formulation of type-specific vaccines. Greece has emerged as a country with high macrolide use. However, studies suggest a gradual reduction in macrolide consumption after 2007.
Methods
During a 7-year period (2011–2017), 604 GAS isolates were recovered from consecutive children presenting with pharyngeal or nonpharyngeal infections in Central Greece; 517 viable isolates underwent molecular analysis, including emm typing.
Results
Isolates belonged to 20 different emm types (in decreasing order of prevalence: 1, 89, 4, 12, 28, 3, 75 and 6, accounting for 88.2% of total isolates). The emm types comprised 10 emm clusters (five most common clusters: E4, A-C3, E1, A-C4 and A-C5). The emm89 isolates were acapsular (‘new clade‘). Overall macrolide resistance rate was 15.4%, and cMLSB emerged as the predominant resistance phenotype (56.4%). The lowest annual resistance rates occurred in 2014 (13.1%), 2016 (5.5%) and 2017(8.0%) (P for trend = 0.002). Consumption of macrolide/lincosamide/streptogramin B declined by 22.6% during 2011–2017. Macrolide resistance and emm28 and emm77 types were associated (both P<0.001). The most frequently identified genetic lineages of macrolide-resistant GAS included emm28/ST52, emm77/ST63, emm12/ST36, emm89/ST101 and emm4/ST39. We estimated that 98.8% of the isolates belonged to emm types incorporated into a novel 30-valent M protein vaccine.
Conclusions
In Central Greece during 2011–2017, the acapsular emm89 isolates comprised the second most prevalent type. Susceptibility testing and molecular analyses revealed decreasing GAS macrolide resistance rates, which may be attributed to the reduction in the consumption of macrolides and/or the reduced circulation of macrolide-resistant clones in recent years. Such data may provide valuable baseline information in targeting therapeutic intervention and the formulation of type-specific GAS vaccines.
Citation: Grivea IN, Syrogiannopoulos GA, Michoula AN, Gazeti G, Malli E, Tsilipounidaki K, et al. (2020) emm Types and clusters and macrolide resistance of pediatric group A streptococcal isolates in Central Greece during 2011-2017. PLoS ONE 15(5): e0232777. https://doi.org/10.1371/journal.pone.0232777
Editor: Jose Melo-Cristino, Universidade de Lisboa Faculdade de Medicina, PORTUGAL
Received: December 1, 2019; Accepted: April 21, 2020; Published: May 7, 2020
Copyright: © 2020 Grivea et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: The author(s) received no specific funding for this work.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Streptococcus pyogenes [group A streptococcus (GAS)] infections are a major cause of morbidity and mortality worldwide and are responsible for a diverse array of noninvasive, invasive and immune-related diseases [1]. The M protein of GAS is a surface protein, which constitutes a major virulence factor as well as the major immunologic epitope of GAS [2]; it possesses a hypervariable region of the amino-terminal, 40 to 50 amino acid residues. The gene encoding M protein, emm, is the basis for sequence typing used to differentiate among strains of GAS; this differentiation is founded on relatively minor sequence differences in the 5ʹ regions of the gene [3–6]. More than 220 emm types have been recognized to date [7].
A protein-based multivalent type-specific GAS vaccine containing amino-terminal M peptides from 30 different emm types has been developed [8, 9]. N-terminal M protein peptides evoke protective antibodies against epidemiologically important GAS types with the greatest bactericidal activity.
Recently, an emm cluster system based on strong phylogenetic support has been described, which serves as a functional classification scheme for GAS M proteins [10]. It adds valuable information regarding tissue tropism and immune response to GAS infections in various settings on a worldwide level. The emm cluster system is expected to guide vaccine design when attempting to predict vaccine efficacy [10–12]. This system has already been used to analyze the GAS epidemiology in the Pacific region, which is characterized by a great variety of circulating emm types [11].
The most important known mechanisms implicated in GAS macrolide resistance are a 14- and 15-membered ring macrolide-specific efflux (M phenotype) [13], encoded by the mef(A) gene [14], as well as the modification of the ribosomal target by a methylase encoded by the erm(B) [15] or the erm(TR) gene [16]; the latter is currently referred to as erm(A) or erm(A) subclass erm(TR) gene [17]. Methylation results in reduced binding of and co-resistance to 14-, 15-, and 16-membered ring macrolide, lincosamide, and streptogramin B (MLSB) antibiotics. Methylase can be expressed either constitutively (cMLSB phenotype) or be induced (iMLSB phenotype).
In several European countries increased rates of macrolide-resistant GAS isolates have been recorded in the late 1990s and early 2000s [18–24]. In these countries various dynamics have been observed −from 2005 onwards− regarding the prevalence of macrolide resistance, the phenotypes involved, the emm type distribution and the circulating resistant clones [25–28]. Therefore, continuing surveillance is in order as it will enhance the understanding of the epidemiology of GAS disease and may guide the formulation of multivalent type-specific vaccines.
High rates of macrolide-resistant GAS isolates have repeatedly been reported in Greece as compared to other European countries [29, 30]. The University General Hospital of Larissa (UGHL) serves as the academic tertiary care referral center for the broader area of Central Greece. The aim of this 7-year survey was to investigate three clinically important epidemiological characteristics of pharyngeal and nonpharyngeal GAS isolates: (i) emm type and cluster distribution, (ii) phenotypes, genotypes and clones of macrolide-resistant isolates, and (iii) proportion of isolates covered by the currently proposed 30-valent M protein-based GAS vaccine.
Materials and methods
Ethics statement
The research protocol was approved by the Ethics Committee of UGHL. The data was analyzed anonymously.
Patients and specimens
During the period January 1, 2011 to December 31, 2017, 604 GAS isolates were recovered from consecutive pediatric patients, who presented with pharyngeal or nonpharyngeal infections to the outpatient pediatric clinic and/or were hospitalized in the pediatric ward of UGHL. The population pool was derived from regions of Central Greece. In UGHL, it has been routine practice to perform swab culture on all children presenting with signs and symptoms of streptococcal pharyngitis. A single GAS isolate was identified in each child; each isolate was derived from a single pharyngeal or nonpharyngeal location.
Antimicrobial susceptibility testing and detection of macrolide resistance genes
Isolates were identified as GAS by colony morphology, β-hemolysis on sheep blood agar and Lancefield grouping with the use of a commercially available agglutination technique (Slidex, Streptokit; BioMérieux, Marcy l' Etoile, France). Susceptibility testing of isolates to benzylpenicillin, erythromycin and clindamycin was performed according to EUCAST guidelines (www.eucast.breakpoints v10.0). D test was performed in order to differentiate the iMLSB phenotype from the M phenotype. The D test was performed by disk diffusion, i.e. by placing a 15 μg erythromycin disk in 12–16 mm proximity to a 2 μg clindamycin disk on a Mueller Hinton agar plate which was inoculated with a streptococcal isolate; detection of the iMLSB phenotype was indicated by blunting of the clindamycin inhibition zone (www.eucast.breakpoints v10.0).
Streptococcal DNA was extracted according to CDC laboratory protocols for S. pyogenes (http://www.cdc.gov/streplab/protocols.html). All viable erythromycin-resistant isolates were screened for the presence of genes associated with macrolide resistance by Polymerase Chain Reaction (PCR), using specific primers for erm(B), erm(TR) and mef(A) [31, 32]. PCR products were separated in 2% gel electrophoresis for 30–40 minutes and samples were identified as positive if they exhibited a clearly defined band at 639bp, 371bp and 348bp, respectively.
Molecular typing
A total of 517 isolates were initially classified by emm typing [6]. The emm types were determined by amplification of the M protein gene (emm). All protocols and assignment of emm types and subtypes were as described in the S. pyogenes emm sequence database (http://www.cdc.gov/streplab/protocol-emm-type.html). The assignment of emm clusters was based on the CDC database (http://www.cdc.gov/streplab/downloads/distribution-emm-types.pdf). The isolates which belonged to type emm89 were further characterized for the presence of the hasABC locus, which is responsible for the synthesis of the hyaluronic acid capsule, using two different pairs of primers. In ‘new clade’ emm89 isolates, which lack the presence of hasABC locus, the first primer pair amplified across the hasABC locus from surrounding genes and generated a product of 632bp; the second primer pair amplified within a 157bp region, which is present only in clade-associated emm89 isolates −in place of the hasABC locus− and generated a product of 127bp [33].
All macrolide-resistant isolates and seven randomly selected macrolide-susceptible ones which belonged to emm89 −each study year represented by one isolate− were, in addition, characterized by Multilocus Sequence Typing (MLST) of the seven gene loci gki, gtr, murI, mutS, recP, xpt and yqiL [34]. The data for S. pyogenes alleles and sequence types (STs) was obtained through the MLST database at https://pubmlst.org/spyogenes/. Electrophoresis, purification and sequencing of PCR products for emm typing and MLST were performed as described above. New emm subtypes and new MLST alleles and STs were assigned to https://www2.cdc.gov/vaccines/biotech/strepblast.asp and https://pubmlst.org/spyogenes/, respectively.
Consumption of macrolides
The IMS Health Xponent Greece database was used to obtain data on the annual outpatient consumption of total and pediatric suspensions packages of macrolide, lincosamide and streptogramin B, expressed in packages per 1,000 inhabitants per day [35].
Statistical analysis
For comparison between groups categorical parameters were compared using the Fisher’s exact test. We assessed the difference between macrolide-susceptible and macrolide-resistant rates for each emm type using the standardized residuals post hoc method for contingency tables [36], with Bonferroni adjustment for multiple comparisons (critical P value = 0.0025). Yearly fluctuation was tested using a log-linear model, in which the exploratory variable was year effect and the response was measured as frequency of a certain emm type or rate of macrolide-resistant isolates. Difference in trend over time was analyzed by the χ2 test for trend, and the association between parameters was analyzed by Pearson’s correlation test. Two-sided tests were used. All analyses were performed with the IBM SPSS software version 25.0 (IBM Corp., Armonk, NY). An effect was considered significant when P<0.05.
Results
Sample and patient characteristics
During the 7-year study period, 604 GAS isolates were obtained from children with pharyngeal and nonpharyngeal GAS infections. Isolates were recovered in 596 noninvasive disease cases from cultures of throat swab (n = 541, 89.6%), skin and soft tissue exudates (pus) (n = 30, 5.0%) or purulent middle ear fluid (n = 25, 4.1%), and in eight cases of invasive disease, from blood (n = 7, 1.1%) or synovial fluid (septic arthritis: n = 1, 0.2%) culture. The diagnoses of seven bacteremic children were meningitis (n = 1), pneumonia with empyema (n = 1), septic arthritis (n = 1), osteomyelitis (n = 1), and bacteremia without an identified focus of infection (n = 3). Children were aged 0.8–15.5 years (median 6.5 years; interquartile range 4.5–9.0 years); 57.9% were boys.
emm types, subtypes and clusters of GAS isolates
Of the total 604 pharyngeal and nonpharyngeal GAS isolates obtained, 517 were viable and underwent emm typing; they belonged to 20 distinct emm types (Table 1). We were able to analyze all eight isolates that were recovered from children with invasive GAS disease in the 7-year study period.
The seven most prominent emm types among pharyngeal isolates were emm1, emm89, emm4, emm12, emm3, emm28 and emm75, which accounted for 84.1% of the isolates. Among nonpharyngeal isolates, the most prominent types were emm1, emm89, emm28, emm12, emm4, emm6, emm2 and emm75, which accounted for 88.5% of isolates in this group.
The 71 emm89 isolates (Table 1) were negative for the presence of the hasABC locus and are therefore presumed to belong to the recently emerged acapsular clade [33].
In addition to subtype emm.0, the following 23 emm subtypes were detected: emm1.16 (n = 2), emm1.24 (n = 1), emm1.119 (n = 1), emm1.120 (n = 2), emm3.1(n = 45), emm3.24 (n = 1), emm3.77 (n = 9), emm3.150 (n = 1), emm4.19 (n = 4), emm6.4 (n = 15), emm6.54 (n = 8), emm6.102 (n = 1), emm11.1(n = 1), emm12.1 (n = 1), emm12.19 (n = 2), emm12.37 (n = 4), emm12.4 (n = 5), emm12.102 (n = 2), emm28.21 (n = 6), emm48.1 (n = 2), emm58.7 (n = 1), emm65.6 (n = 1) and emm118.12 (n = 1).
Analysis of emm types according to emm cluster classification indicated that 20 emm types belonged to 10 emm clusters (Table 1). The most prevalent emm clusters were E4, A-C3, E1, A-C4 and A-C5, which, in total, accounted for 85.1% of all (pharyngeal and nonpharyngeal) isolates.
Yearly fluctuation of emm types and clusters
There was yearly fluctuation in the predominance of emm types. The variation in those belonging to emm4, emm89, emm12 and emm28 reached statistical significance (P<0.05).
Yearly variation of the emm clusters and level of significance in the fluctuation of pharyngeal and nonpharyngeal isolates within these clusters is presented in Fig 1. Significant yearly variation was also observed in clusters E4 (P<0.001), E1 (P<0.001) and A-C4 (P = 0.028) among pharyngeal isolates and in cluster E1 (P = 0.035) among nonpharyngeal noninvasive ones.
Only the five most common (>10%) emm clusters are represented with the P values of their fluctuation, while the five least common ones are grouped as ‘Other’. Numbers inside bars represent emm types in each cluster. Bar colors are explained in the Figure.
Macrolide resistance and MLST typing
All GAS isolates were susceptible to benzylpenicillin. The number of GAS isolates recovered each year, the frequency of macrolide resistance and the resistance phenotypes appear in Fig 2. Over the 7-year study period, resistance to macrolides was found in 93 (15.4%) of the collected 604 GAS isolates. The great majority of these macrolide-resistant isolates (92.9%) exhibited either the iMLSB or the cMLSB phenotype. Specifically, 56.4% had the cMLSB phenotype, 36.5% the iMLSB and 7.1% the M phenotype (Fig 2). There was no significant difference in macrolide resistance between pharyngeal isolates and those from other noninvasive infections (15.5% versus 16.4%, P = 0.85). All eight isolates recovered from children with invasive disease (n = 8) were macrolide susceptible.
The continuous line indicates the national outpatient consumption of total macrolide, lincosamide and streptogramin B packages (tablets and suspensions); IMS data are expressed as packages consumed per 1,000 inhabitants per day. Data on national outpatient consumption in the period preceding the current survey (2007–2010) are given in the text. The resistance phenotype is not available in eight isolates, as follows: 2012 (n = 1), 2014 (n = 3), 2015 (n = 3), and 2016 (n = 1).
During the study period, we noted a tendency towards significant decline in resistance rate (χ2 for trend P = 0.002). The lowest rates occurred in years 2014 (13.1%), 2016 (5.5%) and 2017 (8.0%). However, in 2015, an increased circulation of the cMLSB and iMLSB phenotype isolates was noted (Fig 2).
Among the 517 emm typed isolates, 85 were macrolide-resistant and were genotyped. They belonged to a limited number of emm types (Table 2); emm28, emm77 and emm12 accounted for 81.2% of the macrolide-resistant GAS isolates. Specifically, emm77 and emm4 accounted for 90.3% of the isolates carrying erm(TR), either alone or in combination with mef(A); emm28 and emm12 accounted for 87.5% of those carrying erm(B), either alone or in combination with mef(A) (Table 2).
In 2015, we noted an increase of the four most prominent types of macrolide-resistant isolates (emm77, emm28, emm12 and emm89).
All macrolide-resistant isolates that possessed erm(B), alone or in combination with mef(A), exhibited the cMLSB phenotype. All those possessing erm(TR), alone or in combination with mef(A), exhibited the iMLSB phenotype. Of those harboring solely mef(A) all exhibited the M phenotype.
Macrolide-resistant isolates were studied by MLST (Table 2). One emm28 isolate with a new allele combination, which produced a novel ST −ST1117−, was detected in a pharyngeal sample. In this new ST1117, the gki locus was 99% identical to allele 11. The difference was a T to C alteration at position 458bp and was assigned in MLST database as 172 allele number. All other isolates of each emm type tested belonged to a single ST.
Of the 71 emm89 isolates, 66 (93%) were macrolide susceptible. All seven typed ‘new clade’ macrolide-susceptible isolates were ST101, as follows: one from blood culture (bacteremia without an identified focus of infection), two from skin and soft tissue exudate, one from middle ear exudate and three pharyngeal swab specimens.
Associations between emm types and macrolide resistance
Of the 517 emm typed isolates, those that were macrolide-susceptible belonged to 18 emm types and the macrolide-resistant ones to nine emm types. Statistical associations between emm types and macrolide resistance are depicted in Table 3.
Consumption of macrolides
In addition to the data obtained during the study period on outpatient macrolide, lincosamide and streptogramin B consumption (Fig 2), we also obtained information on their use in previous years. At a national level, a substantial decrease (26.6%) in macrolide, lincosamide and streptogramin B consumption was noted during 2007–2010, from 1.24 packages/1,000 inhabitants/day in 2007 to 0.91 packages/1,000 inhabitants/day in 2010 (authors’ unpublished data based on IMS). The decline in consumption between 2011 and 2017 was 22.6%; the greatest decline (20.4%) was due to the decrease during the 2011–2012 period (Fig 2).
Specifically, the total macrolide, lincosamide and streptogramin B package consumption (tablets and suspensions) in Greece showed a reduction from 3,761,780 in 2011 to 2,835,653 packages in 2017 (24.6% decrease, IMS Health data). During the same time period, consumption of pediatric suspension packages also decreased from 840,981 packages in 2011 to 541,893 packages in 2017 (35.6% decrease, IMS Health data). During the 7-year study period, there was no association between annual macrolide resistance and macrolide, lincosamide and streptogramin B consumption regardless of whether use of both tablets and suspensions (Pearson’s correlation r = 0.404, P = 0.369) or suspensions alone (Pearson’s correlation r = 0.473, P = 0.283) were taken into account.
emm types of the isolates in relation to an experimental 30-valent GAS vaccine
Fifteen of the 20 emm types found in our population −both pharyngeal and nonpharyngeal− are included in the 30-valent M protein vaccine under development. The estimated overall coverage rate by this vaccine is 98.8% of the isolates identified in our study. The cumulative distribution of emm types of GAS isolates is presented in Fig 3. Four hundred and sixty-one (99.1%) of the 465 pharyngeal and 50 (96.2%) of the 52 nonpharyngeal isolates belonged to emm types that are included in the proposed vaccine.
Discussion
We analyzed 604 GAS isolates recovered from pediatric patients in the area of Central Greece from 2011 through 2017. Twenty distinct emm types comprised 10 emm clusters, among which −in descending order− E4, A-C3, E1, A-C4 and A-C5 prevailed. We found that the macrolide resistance rate was 15.4% in this 7-year period. The cMLSB phenotype predominated and persisted throughout the study, whereas the prevalence of M phenotype was 8-fold lower. The present survey updates available information on GAS in this area of Central Greece obtained by the same laboratory in two previous surveys, in 2001 [37] and in 2007–2009 [38]. Our recent findings show lower macrolide resistance rates than those found in the two previous surveys and an inversion of the dominant resistance phenotypes (S1 Table). The resistance rate increased from 19.3% in 2001 [37] to 24% in 2007–2009 [38] and dropped in the recent 2011–2017 survey. Moreover, the macrolide resistance rates in the two most recent years of the present study remained below 10%. Nevertheless, it should not be assumed that Greece has entered a sustained decline of macrolide resistance, and surveillance of GAS macrolide resistance should be maintained in the years to follow before safe conclusions can be drawn.
Among pharyngeal isolates, the seven most frequent emm clusters accounted for 98.5% of the emm typed isolates in this study. Other reports have shown similar results: 96.5% of those collected in a prospective study conducted in Western Greece between 1999 and 2005 [29], authors’ unpublished data], 92.8% of 7,040 isolates in the USA [12] and 91.1% in Southern Taiwan [39]. High diversity of emm types has been observed in the Pacific region and other low-income settings [11, 40].
A clade of emm89 isolates lacking the hyaluronic acid capsule biosynthesis locus has emerged in the 2000s and spread throughout several countries in North America and Europe [33, 41]. Seventy-one emm89 isolates were tested and were all inferred to belong to the new clade based on the absence of capsule locus. Seven percent of these was macrolide resistant and exhibited various phenotypes and resistance genotypes. The 12 tested isolates −seven macrolide-susceptible and five macrolide-resistant− belonged to ST101, which is the most common type among the acapsular clade [41]. Given the high prevalence of emm89 in this study, it appears that this clade has also become dominant in central Greece and macrolide resistance has emerged in selected isolates.
A significant association was found between macrolide susceptibility and emm1. We also established significant associations between macrolide resistance and the emm28 and emm77 types. Other studies report a similar correlation for emm28 [25–27, 30, 41] and emm77 [25, 28, 30, 42].
The distribution of our most common macrolide-resistant clones, −emm28/ST52, emm77/ST63, emm12/ST36, emm89/ST101 and emm4/ST39− varied widely in different European regions over time, while certain macrolide-resistant lineages frequently encountered in some European countries −emm22/ST46, emm11/ST403 and emm1/ST28− were either absent or uncommon in our study [25–28, 43–46]. Of note, in our survey, macrolide-resistant emm4/ST39 GAS lineage exhibited iMLSB phenotype, whereas in other European countries, this lineage has been most frequently found among resistant isolates which expressed the M phenotype [25, 26, 43–45].
Temporal changes in clonal composition have been noted in Central Greece. In the previous survey conducted in the same region between January 2007 and June 2009, no isolates harboring erm(B) gene and expressing cMLSB phenotype occurred; emm28/ST52 and emm58/ST176 were absent, whereas emm89/ST101, emm12/ST36 and emm75/ST150 were identified among macrolide-susceptible isolates [38]. The lineages emm94/ST89 and emm77/ST550 identified in the previous survey were not detected in this recent one (S1 Table). High clonal instability has been reported in Portugal from 2007 to 2013 [27].
Greece has been, and remains, an outstanding macrolide consumer in Europe [47]. At a national level, a gradual decline in the outpatient macrolide, lincosamide and streptogramin B consumption was noted between 2007 and 2010 and was further enhanced in 2012. From 2013 onwards, there has been limited change. A 0.67 value of packages/1,000 inhabitants/day of macrolide, lincosamide, streptogramin B and tetracycline has been proposed as representing a critical threshold. It was shown that when the use of these classes of antimicrobial agents was maintained below this critical threshold, a low level of macrolide-resistant GAS isolates and an increased number of erm(A)-harboring emm77 GAS with low fitness costs occurred [35]. The national consumption has not been lowered to this threshold in Greece. This may provide an explanation for the dominance of the erm(B)-associated cMLSB phenotype in Greece despite the reduction in the macrolide resistance rate. The small (2%) difference in macrolide, lincosamide and streptogramin B consumption between the two calculation methods, i.e., crude number of packages versus packages/1,000 inhabitants/day index, is due to small population changes during the study years and the logistics of the calculations.
In the surveillance period of our study the lowest rates of macrolide resistance were observed in 2014 (13.1%), 2016 (5.5%) and 2017 (8.0%). This decrease may be attributed to the gradual reduction in the consumption of macrolides, which commenced in the pre-study period, i.e. 2007–2010, in association with the dynamics of specific major macrolide-resistant clones [26, 27]. We speculate that the niche of macrolide-resistant clones that remains in a community after such reduction may occasionally lead to a temporary increased spread of these clones and cause a transient surge of macrolide resistance rate, similar to the one observed in our study during 2015. The subsequent nadir rates observed in 2016 and 2017 are echoed by the results of other independent studies including an abrupt drop from a high rate of approximately 22% in 2012 to a low one of 10.95% in 2013 reported by a major Children’s Hospital in Athens [30], and a similar reduction from a high rate of 16.3% in 2008 to a low one of 5.1% in 2009 in a region of Central Italy [28].
Finally, according to our study data, a 98.8% coverage of the isolates would be achieved by the 30-valent GAS vaccine under development. A 93.7% coverage by the same vaccine has been estimated for northern Spain [48]. It is evident that the potential efficacy of vaccines based on the M protein is limited by their suitability to local epidemiological characteristics.
There are certain limitations to this study. First, it is a single-center study. However, our tertiary care Center serves a substantial proportion of the pediatric population of Central Greece. Second, in our analysis we compared the regional macrolide resistance rate with the total national macrolide consumption, as it was not feasible to obtain an accurate estimation of the consumption −of pediatric packaging (suspension), in particular− at the local level. On the other hand, there are important strengths to our study. It is comparable to two previous surveys and all three studies span a period of 17 years (2001–2017). Indeed, these surveys have shared the same region of surveillance and method of sampling, and all samples have been analyzed by the same laboratory.
In conclusion, in Central Greece during 2011–2017, macrolide resistance rate of GAS isolates decreased to 15.4% as compared to 19.3% in 2001 and 24% in 2007–2009; it remained below 10% in the last two years of the study. During the study period the cMLSB phenotype predominated. A relatively limited number of emm types prevailed among GAS isolates recovered from children with pharyngeal and nonpharyngeal infections. A yearly fluctuation in the frequency of emm types and clusters was evident. The most frequently found genetic lineages of macrolide-resistant GAS isolates include those defined by emm28/ST52, emm77/ST63, emm12/ST36, emm89/ST101 and emm4/ST39. Interestingly, 98.8% of the GAS isolates recovered in our population belonged to emm types that are included in the 30-valent M protein vaccine under development. We propose that molecular surveillance of GAS isolates at different time-periods and in various areas, especially in the event of increased prevalence of macrolide-resistant GAS, may guide treatment strategies and the formulation and targeting of multivalent type-specific vaccines.
Supporting information
S1 Table. Comparison of GAS isolates recovered in Central Greece in 2001, 2007–2009 and 2011–2017.
https://doi.org/10.1371/journal.pone.0232777.s001
(PDF)
References
- 1.
Stevens DL. Group A beta-hemolytic streptococci: virulence factors, pathogenesis, and spectrum of clinical infections. In: Stevens DL, Kaplan EL, editors. Streptococcal infections: clinical aspects, microbiology, and molecular pathogenesis. New York: Oxford University Press. 2000. pp. 19–36.
- 2. Lancefield RC. Current knowledge of type-specific M antigens of group A streptococci. J Immunol. 1962;89: 307–313. pmid:14461914
- 3. Beachey EH, Seyer JM, Dale JB, Simpson WA, Kang AH. Type-specific protective immunity evoked by synthetic peptide of Streptococcus pyogenes M protein. Nature 1981;292: 457–459. pmid:6166868
- 4. Jones KF, Manjula BN, Johnston KH, Hollingshead SK, Scott JR, Fischetti VA. Location of variable and conserved epitopes among the multiple serotypes of streptococcal M protein. J Exp Med. 1985;161: 623–628. pmid:2579187
- 5. Dale JB, Beachey EH. Localization of protective epitopes of the amino terminus of type 5 streptococcal M protein. J Exp Med. 1986;163: 1191–1202. pmid:2422314
- 6. Beall B, Facklam R, Thompson T. Sequencing emm-specific PCR products for routine and accurate typing of group A streptococci. J Clin Microbiol. 1996;34: 953–958. pmid:8815115
- 7. McMillan DJ, Drèze PA, Vu T, Bessen DE, Guglielmini J, Steer AC, et al. Updated model of group A Streptococcus M proteins based on a comprehensive worldwide study. Clin Microbiol Infect. 2013;19: E222–E229. pmid:23464795
- 8. Dale JB, Penfound TA, Chiang EY, Walton WJ. New 30-valent M protein-based vaccine evokes cross-opsonic antibodies against non-vaccine serotypes of group A streptococci. Vaccine. 2011;29: 8175–8178. pmid:21920403
- 9. Dale JB, Penfound TA, Tamboura B, Sow SO, Nataro JP, Tapia M, et al. Potential coverage of a multivalent M protein-based group A streptococcal vaccine. Vaccine. 2013;31: 1576–1581. pmid:23375817
- 10. Sanderson-Smith M, De Oliveira DM, Guglielmini J, McMillan DJ, Vu T, Holien JK, et al. A systematic and functional classification of Streptococcus pyogenes that serves as a new tool for molecular typing and vaccine development. J Infect Dis. 2014;210: 1325–1338. pmid:24799598
- 11. Baroux G, D’Ortenzio E, Amédéo N, Baker C, Ali Alsuwayyid B, Dupont-Rouzeyrol M, et al. The emm-cluster typing system for Group A Streptococcus identifies epidemiologic similarities across the Pacific region. Clin Infect Dis. 2014;59, e84–e92. pmid:24965347
- 12. Shulman ST, Tanz RR, Dale JB, Steer AC, Smeesters PR. Added value of the emm-cluster typing system to analyze group A Streptococcus epidemiology in high-income settings. Clin Infect Dis. 2014;59: 1651–1652. pmid:25115872
- 13. Sutcliffe J, Tait-Kamradt A, Wondrack L. Streptococcus pneumoniae and Streptococcus pyogenes resistant to macrolides but sensitive to clindamycin: a common resistance pattern mediated by an efflux system. Antimicrob Agents Chemother. 1996;40: 1817–1824. pmid:8843287
- 14. Clancy J, Petitpas J, Dib-Hajj F, Yuan W, Cronan M, Kamath AV et al. Molecular cloning and functional analysis of a novel macrolide-resistance determinant, mefA, from Streptococcus pyogenes. Mol Microbiol. 1996;22: 867–879. pmid:8971709
- 15. Weisblum B. Erythromycin resistance by ribosome modification. Antimicrob Agents Chemother. 1995;39: 577–585. pmid:7793855
- 16. Seppälä H, Skurnik M, Soini H, Roberts MC, Huovinen P. A novel erythromycin resistance methylase gene (ermTR) in Streptococcus pyogenes. Antimicrob Agents Chemother. 1998;42: 257–262. pmid:9527769
- 17. Roberts MC, Sutcliffe J, Courvalin P, Jensen LB, Rood J, Seppälä H. Nomenclature for macrolide and macrolide-lincosamide-streptogramin B resistance determinants. Antimicrob Agents Chemother. 1999;43: 2823–2830. pmid:10582867
- 18. Bingen E, Bidet P, Mihaila-Amrouche L, Doit C, Forcet S, Brahimi N, et al. Emergence of macrolide-rresistant Streptococcus pyogenes strains in French children. Antimicrob Agents Chemother. 2004;48: 3559–3562. pmid:15328126
- 19. Szczypa K, Sadowy E, Izdebski R, Hryniewicz W. A rapid increase in macrolide resistance in Streptococcus pyogenes isolated in Poland during 1996–2002. J Antimicrob Chemother. 2004;54: 828–831. pmid:15329367
- 20. Silva-Costa C, Ramirez M, Melo-Cristino J, the Portuguese Surveillance Group for the Study of Respiratory Pathogens. Rapid inversion of the prevalences of macrolide resistance phenotypes paralleled by a diversification of T and emm types among Streptococcus pyogenes in Portugal. Antimicrob Agents Chemother. 2005;49: 2109–2111. pmid:15855540
- 21. Malhotra-Kumar S, Lammens C, Chapelle S, Wijdooghe M, Piessens J, Van Herck K, et al. Macrolide- and telithromycin-resistant Streptococcus pyogenes, Belgium, 1999–2003. Emerg Infect Dis. 2005;11: 939–942. pmid:15963292
- 22. Čižman M, Beović B, Seme K, Paragi M, Štrumbelj I, Müller-Premru M, et al. Macrolide resistance rates in respiratory pathogens in Slovenia following reduced macrolide use. Int J Antimicrob Agents 2006;28: 537–542. pmid:17101264
- 23. Mazzariol A, Koncan R, Bahar G, Cornaglia G. Susceptibilities of Streptococcus pyogenes and Streptococcus pneumoniae to macrolides and telithromycin: data from an Italian multicenter study. J Chemother. 2007;19: 500–507. pmid:18073148
- 24. Rubio-López V, Valdezate S, Álvarez D, Villalón P, Medina MJ, Salcedo C, et al. Molecular epidemiology, antimicrobial susceptibilities and resistance mechanisms of Streptococcus pyogenes isolates resistant to erythromycin and tetracycline in Spain (1994–2006). BMC Microbiol. 2012;12: 215. pmid:22998619
- 25. d’Humières C, Cohen R, Levy C, Bidet P, Thollot F, Wollner A, et al. Decline in macrolide-resistant Streptococcus pyogenes isolates from French children. Int J Med Microbiol. 2012;302: 300–303. pmid:23103047
- 26. Montes M, Tamayo E, Mojica C, García-Arenzana JM, Esnal O, Pérez-Trallero E. What causes decreased erythromycin resistance in Streptococcus pyogenes? Dynamics of four clones in a southern European region from 2005 to 2012. J Antimicrob Chemother. 2014;69: 1474–1482. pmid:24562616
- 27. Silva-Costa C, Ramirez M, Melo-Cristino J; the Portuguese Group for the Study of Streptococcal Infections. Declining macrolide resistance in Streptococcus pyogenes in Portugal (2007–13) was accompanied by continuous clonal changes. J Antimicrob Chemother. 2015;70: 2729–2733. pmid:26142474
- 28. Olivieri R, Morandi M, Zanchi A, Tordini G, Pozzi G, De Luca A, et al. Evolution of macrolide resistance in Streptococcus pyogenes over 14 years in an area of central Italy. J Med Microbiol. 2015;64: 1186–1195. pmid:26224594
- 29. Syrogiannopoulos GA, Grivea IN, Al-Lahham A, Panagiotou M, Tsantouli AG, Michoula AN, et al. Seven-year surveillance of emm types of pediatric group A streptococcal pharyngitis isolates in western Greece. PLoS One 2013;8: e71558. pmid:23977078
- 30. Michos A, Koutouzi FI, Tsakris A, Chatzichristou P, Koutouzis EI, Daikos GL, et al. Molecular analysis of Streptococcus pyogenes macrolide resistance of paediatric isolates during a 7 year period (2007–13). J Antimicrob Chemother. 2016;71: 2113–2117. pmid:27118782
- 31. Sutcliffe J, Grebe T, Tait-Kamradt A, Wondrack L. Detection of erythromycin-resistant determinants by PCR. Antimicrob Agents Chemother. 1996;40: 2562–2566. pmid:8913465
- 32. De Azevedo JCS, Yeung RH, Bast DJ, Duncan CL, Borgia SB, Low DE. Prevalence and mechanisms of macrolide resistance in clinical isolates of group A streptococci from Ontario, Canada. Antimicrob Agents Chemother. 1999;43: 2144–2147. pmid:10471555
- 33. Turner CE, Abbott J, Lamagni T, Holden MTG, David S, Jones MD, et al. Emergence of a new highly successful acapsular group A Streptococcus clade of genotype emm89 in the United Kingdom. mBio 2015;6: e00622–15. pmid:26173696
- 34. Enright MC, Spratt BG, Kalia A, Cross JS, Bessen DE. Multilocus sequence typing of Streptococcus pyogenes and the relationships between emm type and clone. Infect Immun. 2001;69: 2416–2427. pmid:11254602
- 35. Van Heirstraeten L, Coenen S, Lammens C, Hens N, Goossens H, Malhotra-Kumar S. Antimicrobial drug use and macrolide-resistant Streptococcus pyogenes, Belgium. Emerg Infect Dis. 2012;18: 1515–1518. pmid:22932671
- 36. Beasley TM, Schumacker RE. Multiple regression approach to analyzing contingency tables: post hoc and planned comparison procedures. J Exp Educ. 1995;64: 79–93.
- 37. Petinaki E, Kontos F, Pratti A, Skulakis C, Maniatis AN. Clinical isolates of macrolide-resistant Streptococcus pyogenes in Central Greece. Int J Antimicrob Agents. 2003;21: 67–70. pmid:12507840
- 38. Malli E, Tatsidou E, Damani A, Pantelidi K, Petinaki E, Skoulakis C, et al. Macrolide-resistant Streptococcus pyogenes in Central Greece: prevalence; mechanism and molecular identification. Int J Antimicrob Agents. 2010;35: 614–615. pmid:20206477
- 39. Chiang-Ni C, Zheng PX, Wang SY, Tsai PJ, Chuang WJ, Lin YS, et al. Epidemiology analysis of Streptococcus pyogenes in a hospital in southern Taiwan by use of the updated emm cluster typing system. J Clin Microbiol. 2016;54: 157–162. pmid:26560544
- 40. Abdissa A, Asrat D, Kronvall G, Shittu B, Achiko D, Zeidan M, et al. High diversity of group A streptococcal emm types among healthy schoolchildren in Ethiopia. Clin Infect Dis. 2006;42: 1362–1367. pmid:16619146
- 41. Friães A, Machado MP, Pato C, Carriço J, Melo-Cristino J, Ramirez M. Emergence of the same successful clade among distinct populations of emm89 Streptococcus pyogenes in multiple geographic regions. mBio 2015;6: e01780–15. pmid:26628724
- 42. Koutouzi F, Tsakris A, Chatzichristou P, Koutouzis E, Daikos GL, Kirikou E, et al. Streptococcus pyogenes emm types and clusters during a 7-year period (2007 to 2013) in pharyngeal and nonpharyngeal pediatric isolates. J Clin Microbiol. 2015; 53: 2015–2021. pmid:25878351
- 43. Reinert RR, Lütticken R, Sutcliffe JA, Tait-Kamradt A, Cil MY, Schorn HM, et al. Clonal relatedness of erythromycin-resistant Streptococcus pyogenes isolates in Germany. Antimicrob Agents Chemother. 2004;48: 1369–1373. pmid:15047546
- 44. Silva-Costa C, Ramirez M, Melo-Cristino J. Identification of macrolide-resistant clones of Streptococcus pyogenes in Portugal. Clin Microbiol Infect. 2006;12: 513–518. pmid:16700698
- 45. Silva-Costa C, Pinto FR, Ramirez M, Melo-Cristino J; the Portuguese Surveillance Group for the Study of Respiratory Pathogens. Decrease in macrolide resistance and clonal instability among Streptococcus pyogenes in Portugal. Clin Microbiol Infect. 2008;14: 1152–1159. pmid:19046174
- 46. Littauer P, Caugant DA, Sangvik M, Høiby EA, Sundsfjord A, Simonsen GS. Macrolide-resistant Streptococcus pyogenes in Norway: population structure and resistance determinants. Antimicrob Agents Chemother. 2006;50: 1896–1899. pmid:16641473
- 47. Adriaenssens N, Coenen S, Versporten A, Muller A, Minalu G, Faes C, et al. European Surveillance of Antimicrobial Consumption (ESAC): outpatient macrolide, lincosamide and streptogramin (MLS) use in Europe (1997–2009). J Antimicrob Chemother. 2011;66 (Suppl 6): vi37–vi45.
- 48. Tamayo E, Montes M, Garcia-Arenzana JM, Pérez-Trallero E. Streptococcus pyogenes emm-types in northern Spain; population dynamics over a 7-year period. J Infect. 2014;68: 50–57. pmid:23999149