Genomic epidemiology of penicillin-non-susceptible Streptococcus pneumoniae

Hybrid sequence-based analysis reveals the distribution of bacterial species and genes in the oral microbiome at a high resolution

Bacteria in the oral microbiome are poorly identified owing to the lack of established culture methods for them. Thus, this study aimed to use culture-free analysis techniques, including bacterial single-cell genome sequencing, to identify bacterial species and investigate gene distribution in saliva. Saliva samples from the same individual were classified as inactivated or viable and then analyzed using 16S rRNA sequencing, metagenomic shotgun sequencing, and bacterial single-cell sequencing. The results of 16S rRNA sequencing revealed similar microbiota structures in both samples, with Streptococcus being the predominant genus. Metagenomic shotgun sequencing showed that approximately 80 % of the DNA in the samples was of non-bacterial origin, whereas single-cell sequencing showed an average contamination rate of 10.4 % per genome. Single-cell sequencing also yielded genome sequences for 43 out of 48 wells for the inactivated samples and 45 out of 48 wells for the viable samples. With respect to resistance genes, four out of 88 isolates carried cfxA, which encodes a β-lactamase, and four isolates carried erythromycin resistance genes. Tetracycline resistance genes were found in nine bacteria. Metagenomic shotgun sequencing provided complete sequences of cfxA, ermF, and ermX, whereas other resistance genes, such as tetQ and tetM, were detected as fragments. In addition, virulence factors from Streptococcus pneumoniae were the most common, with 13 genes detected. Our average nucleotide identity analysis also suggested five single-cell-isolated bacteria as potential novel species. These data would contribute to expanding the oral microbiome data resource.

Changes in pneumococcal serotypes distribution and penicillin resistance in healthy children five years after generalization of PCV10

Objective

Streptococcus pneumoniae (S. pneumoniae) nasopharyngeal carriage has significantly decreased after the generalization of pneumococcal vaccination worldwide. This study sought to investigate changes in S. pneumoniae carriage rates, serotype distribution and penicillin non-susceptibility following the generalization of 10-valent pneumococcal conjugate vaccine.

Methods

A prospective study was conducted in Marrakesh, Morocco, between 2017 and 2018, among healthy children attending vaccination centers. We collected nasopharyngeal swabs and questionnaire data for each child. Using univariate logistic regression, we analyzed the association between S. pneumoniae carriage and various risk factors. Comparisons of serotype diversity and penicillin resistance between 2017 and 2018 and the period before introduction of vaccination (2008-2009, n = 660) were performed using Simpson index and the chi-squared test, respectively.

Results

During 2017-2018, 515 children aged between 6 and 36 months participated. The S. pneumoniae carriage rate was 43.3%. Looking at the distribution serotypes, the rate of PCV10 serotypes rate was only 9.6%. Among non-vaccine serotypes, an increase in serotypes 6C/6D (22; 14%), 19B/19C (17; 10.8%), and 15B/15C (11; 7%) was observed. A particular increase in serotype diversity was also observed after the generalization of PCV10 (p < 0.001). S. pneumoniae non-susceptible to penicillin decreased, reaching a rate of 26.6% in 2017-2018.

Conclusion

The significant change in S. pneumoniae carriage, serotype distribution, and penicillin resistance highlights the effectiveness of the pneumococcal conjugate vaccine among children in Marrakesh, Morocco.

Streptococcus suis serotype 9 in Italy: genomic insights into high-risk clones with emerging resistance to penicillin

BACKGROUND

Streptococcus suis is an important pig pathogen and an emerging zoonotic agent. In a previous study, we described a high proportion of penicillin-resistant serotype 9 S. suis (SS9) isolates on pig farms in Italy.

OBJECTIVES

We hypothesized that resistance to penicillin emerged in some SS9 lineages characterized by substitutions at the PBPs, contributing to the successful spread of these lineages in the last 20 years.

METHODS

Sixty-six SS9 isolates from cases of streptococcosis in pigs were investigated for susceptibility to penicillin, ceftiofur and ampicillin. The isolates were characterized for ST, virulence profile, and antimicrobial resistance genes through WGS. Multiple linear regression models were employed to investigate the associations between STs, year of isolation, substitutions at the PBPs and an increase in MIC values to β-lactams.

RESULTS

MIC values to penicillin increased by 4% each year in the study period. Higher MIC values for penicillin were also positively associated with ST123, ST1540 and ST1953 compared with ST16. The PBP sequences presented a mosaic organization of blocks. Within the same ST, substitutions at the PBPs were generally more frequent in recent isolates. Resistance to penicillin was driven by substitutions at PBP2b, including K479T, D512E and K513E, and PBP2x, including T551S, while reduced susceptibility to ceftiofur and ampicillin were largely dependent on substitutions at PBP2x.

CONCLUSIONS

Here, we identify the STs and substitutions at the PBPs responsible for increased resistance of SS9 to penicillin on Italian pig farms. Our data highlight the need for monitoring the evolution of S. suis in the coming years.

Genomic epidemiology of Streptococcus pneumoniae serotype 16F lineages

Due to the emergence of non-vaccine serotypes in vaccinated populations, Streptococcus pneumoniae remains a major global health challenge despite advances in vaccine development. Serotype 16F is among the predominant non-vaccine serotypes identified among vaccinated infants in South Africa (SA). To characterize lineages and antimicrobial resistance in 16F isolates obtained from South Africa and place the local findings in a global context, we analysed 10 923 S. pneumoniae carriage isolates obtained from infants recruited as part of a broader SA birth cohort. We inferred serotype, resistance profile for penicillin, chloramphenicol, cotrimoxazole, erythromycin and tetracycline, and global pneumococcal sequence clusters (GPSCs) from genomic data. To ensure global representation, we also included S. pneumoniae carriage and disease isolates from the Global Pneumococcal Sequencing (GPS) project database (n=19 607, collected from 49 countries across 5 continents, 1995-2018, accessed 17 March 2022). Nine per cent (934/10923) of isolates obtained from infants in the Drakenstein community in SA and 2 %(419/19607) of genomes in the GPS dataset were serotype 16F. Serotype 16F isolates were from 28 different lineages of S. pneumoniae, with GPSC33 and GPSC46 having the highest proportion of serotype 16F isolates at 26 % (346/1353) and 53 % (716/1353), respectively. Serotype 16F isolates were identified globally, but most isolates were collected from Africa. GPSC33 was associated with carriage [OR (95 % CI) 0.24 (0.09-0.66); P=0.003], while GPSC46 was associated with disease [OR (95 % CI) 19.9 (2.56-906.50); P=0.0004]. Ten per cent (37/346) and 15 % (53/346) of isolates within GPSC33 had genes associated with resistance to penicillin and co-trimoxazole, respectively, and 18 % (128/716) of isolates within GPSC46 had genes associated with resistance to co-trimoxazole. Resistant isolates formed genetic clusters, which may suggest emerging resistant lineages. Serotype 16F lineages were common in southern Africa. Some of these lineages were associated with disease and resistance to penicillin and cotrimoxazole. We recommend continuous genomic surveillance to determine the long-term impact of serotype 16F lineages on vaccine efficacy and antimicrobial therapy globally. Investing in vaccine strategies that offer protection over a wide range of serotypes/lineages remains essential. This paper contains data hosted by Microreact.

Evidence of virulence and antimicrobial resistance in Streptococcus pneumoniae serotype 16F lineages

Introduction

Due to the emergence of non-vaccine serotypes in vaccinated populations, Streptococcus pneumoniae remains a major global health challenge despite advances in vaccine development. Serotype 16F is among the predominant non-vaccine serotypes identified among vaccinated infants in South Africa (SA).

Aim

To characterise lineages and antimicrobial resistance in 16F isolates obtained from South Africa and placed the local findings in a global context.

Methodology

We analysed 10923 S. pneumoniae carriage isolates obtained from infants recruited as part of a broader SA birth cohort. We inferred serotype, resistance profile for penicillin, chloramphenicol, cotrimoxazole, erythromycin and tetracycline, and Global Pneumococcal Sequence Clusters (GPSCs) from genomic data. To ensure global representation, we also included S. pneumoniae carriage and disease isolates from the Global Pneumococcal Sequencing (GPS) project database (n=19,607, collected from 49 countries across five continents, years covered (1995 - 2018), accessed on 17 th March 2022).

Results

Nine percent (934/10923) of isolates obtained from infants in the Drakenstein community in SA and 2% (419/19607) of genomes in the GPS dataset were serotype 16F. Serotype 16F isolates were from 28 different lineages of S. pneumoniae, with GPSC33 and GPSC46 having the highest proportion of serotype 16F isolates at 26% (346/1353) and 53% (716/1353), respectively. Serotype 16F isolates were identified globally, however, most isolates were collected from Africa. GPSC33 was associated with carriage [OR (95% CI) 0.24 (0.09 - 0.66); p=0.003], while GPSC46 was associated with disease [OR (95% CI) 19.9 (2.56 - 906.50); p=0.0004]. 10% (37/346) and 15% (53/346) of isolates within GPSC33 had genes associated with resistance to penicillin and co-trimoxazole, respectively, and 18% (128/716) of isolates within GPSC46 had genes associated with resistance to co-trimoxazole. Resistant isolates formed genetic clusters which may suggest emerging resistant lineages.

Discussion

Serotype 16F lineages are common in Southern Africa. Some of these lineages are associated with disease, and resistance to penicillin and cotrimoxazole. We recommend continuous genomic surveillance to determine long term impact of serotype 16F lineages on vaccine efficacy and antimicrobial therapy globally. Investing in vaccine strategies that offer protection over a wide range of serotypes/lineages remains essential.

DATA SUMMARY

The sequencing reads for the genomes analysed have been deposited in the European Nucleotide Archive and the accession numbers for each isolate are listed in Supplementary Table1 . Phylogenetic tree of serotype 16F pneumococcal genomes and associated metadata are available for download and visualisation on the Microreact website: Phylogenies of seotype 16F, GPSC33 and GPSC46 are available on the Microreact serotype-16F , GPSC33 and GPSC46 , respectively.

IMPACT STATEMENT

This study shows that serotype 16F lineages are predominant in Southern Africa and are associated with disease and antimicrobial resistance. Although serotype 16F has been included in the newer formulation of the upcoming vaccine formulations of PCV21 and IVT-25, continuous surveillance to determine long term impact of serotype 16F lineages on vaccines and antimicrobial therapy remains essential.

Molecular Characterization and Antibiotic Susceptibility of Non-PCV13 Pneumococcal Serotypes among Vaccinated Children in Cape Coast, Ghana

Preventive strategies involving the use of pneumococcal conjugate vaccines (PCVs) are known to drastically reduce pneumococcal disease. However, PCV vaccination has been plagued with serotype replacement by non-PCV serotypes. In this study, we describe the prevalence and molecular characteristics of non-PCV13 serotypes (non-vaccine serotypes, NVTs) from pneumococcal carriage isolates obtained from children < 5 years old in Cape Coast, Ghana, after PCV introduction. The isolates were subjected to antibiotic susceptibility testing and multilocus sequence typing (MLST), and molecular techniques were used to detect the presence of virulence genes. Serotypes 11A, 13, 15B, 23B, and 34 formed the top five of the 93 NVT isolates. As such, 20 (21.5%), 49 (48.4%), and 70 (74.3%) isolates were non-susceptible to penicillin, tetracycline, and cotrimoxazole, respectively. Sixteen (17.2%) multidrug-resistant isolates were identified. However, non-susceptibility to ceftriaxone and erythromycin was low and all isolates were fully susceptible to levofloxacin, linezolid, and vancomycin. Whereas pcpA, pavB, lytA, and psrP genes were detected in nearly all serotypes, pilus islet genes were limited to serotypes 11A, 13, and 23B. MLST for predominant serotype 23B isolates revealed three known and seven novel sequence types (STs). ST172 and novel ST15111 were the most dominant and both STs were related to PMEN clone Columbia23F-26 (ST338). In conclusion, non-PCV13 serotype 23B was the most prevalent, with characteristics of rapid clonal expansion of ST172 and ST15111, which are related to international clones of the pneumococcus. Continuous monitoring of NVTs in Ghana is, therefore, essential, as they have the potential to cause invasive disease, show high antibiotic resistance, and attenuate the effects of PCV vaccination.

A Chimeric Penicillin Binding Protein 2X Significantly Decreases in Vitro Beta-Lactam Susceptibility and Increases in Vivo Fitness of Streptococcus pyogenes

All tested strains of Streptococcus pyogenes (group A streptococcus, GAS) remain susceptible to penicillin. However, GAS strains with amino acid substitutions in penicillin-binding proteins that confer decreased susceptibility to beta-lactam antibiotics have been identified recently. This discovery raises concerns about emergence of beta-lactam antibiotic resistance in GAS. Whole genome sequencing recently identified GAS strains with a chimeric penicillin-binding protein 2X (PBP2X) containing a recombinant segment from Streptococcus dysgalactiae subspecies equisimilis (SDSE). To directly test the hypothesis that the chimeric SDSE-like PBP2X alters beta-lactam susceptibility in vitro and fitness in vivo, an isogenic mutant strain was generated and virulence assessed in a mouse model of necrotizing myositis. Compared with naturally occurring and isogenic strains with a wild-type GAS-like PBP2X, strains with the chimeric SDSE-like PBP2X had reduced susceptibility in vitro to nine beta-lactam antibiotics. In a mouse model of necrotizing myositis, the strains had identical fitness in the absence of benzylpenicillin treatment. However, mice treated intermittently with a subtherapeutic dose of benzylpenicillin had significantly more colony-forming units recovered from limbs infected with strains with the chimeric SDSE-like PBP2X. These results show that mutations such as the PBP2X chimera may result in significantly decreased beta-lactam susceptibility and increased fitness and virulence. Expanded diagnostic laboratory surveillance, genome sequencing, and molecular pathogenesis study of potentially emergent beta-lactam antibiotic resistance among GAS are needed.

The role of interspecies recombination in the evolution of antibiotic-resistant pneumococci

Multidrug-resistant Streptococcus pneumoniae emerge through the modification of core genome loci by interspecies homologous recombinations, and acquisition of gene cassettes. Both occurred in the otherwise contrasting histories of the antibiotic-resistant S. pneumoniae lineages PMEN3 and PMEN9. A single PMEN3 clade spread globally, evading vaccine-induced immunity through frequent serotype switching, whereas locally circulating PMEN9 clades independently gained resistance. Both lineages repeatedly integrated Tn916-type and Tn1207.1-type elements, conferring tetracycline and macrolide resistance, respectively, through homologous recombination importing sequences originating in other species. A species-wide dataset found over 100 instances of such interspecific acquisitions of resistance cassettes and flanking homologous arms. Phylodynamic analysis of the most commonly sampled Tn1207.1-type insertion in PMEN9, originating from a commensal and disrupting a competence gene, suggested its expansion across Germany was driven by a high ratio of macrolide-to-β-lactam consumption. Hence, selection from antibiotic consumption was sufficient for these atypically large recombinations to overcome species boundaries across the pneumococcal chromosome.

Inferring genetic fitness from genomic data

The genetic composition of a naturally developing population is considered as due to mutation, selection, genetic drift, and recombination. Selection is modeled as single-locus terms (additive fitness) and two-loci terms (pairwise epistatic fitness). The problem is posed to infer epistatic fitness from population-wide whole-genome data from a time series of a developing population. We generate such data in silico and show that in the quasilinkage equilibrium phase of Kimura, Neher, and Shraiman, which pertains at high enough recombination rates and low enough mutation rates, epistatic fitness can be quantitatively correctly inferred using inverse Ising-Potts methods.

[Current microbiological aspects of community respiratory infection beyond COVID-19]

From a microbiological point of view, both empirical and targeted antimicrobial treatment in respiratory infection is based on the sensitivity profile of isolated microorganisms and the possible resistance mechanisms that they may present. The latter may vary in different geographic areas according to prescription profiles and vaccination programs. Beta-lactam antibiotics, fluoroquinolones, and macrolides are the most commonly used antimicrobials during the exacerbations of chronic obstructive pulmonary disease and community-acquired pneumonia. In their prescription, different aspects such as intrinsic activity, bactericidal effect or their ability to prevent the development of resistance must be taken into account. The latter is related to the PK/PD parameters, the mutant prevention concentration and the so-called selection window. More recently, the potential ecological impact has grown in importance, not only on the intestinal microbiota, but also on the respiratory one. Maintaining the state of eubiosis requires the use of antimicrobials with a low profile of action on anaerobic bacteria. With their use, the resilience of the bacterial populations belonging to the microbiota, the state of resistance of colonization and the collateral damage related to the emergence of resistance to the antimicrobials in pathogens causing the infections and in the bacterial populations integrating the microbiota.

The Emergence of Invasive Streptococcus pneumoniae Serotype 24F in Lebanon: Complete Genome Sequencing Reveals High Virulence and Antimicrobial Resistance Characteristics

Background

Invasive pneumococcal disease (IPD) remains a global health problem. IPD incidence has significantly decreased by the use of pneumococcal conjugate vaccines (PCV). Nevertheless, non-PCV serotypes remain a matter of concern. Eight Streptococcus pneumoniae serotype 24F isolates, belonging to a non-PCV serotype, were detected through the Lebanese Inter-Hospital Pneumococcal Surveillance Program. The aim of the study is to characterize phenotypic and genomic features of the 24F isolates in Lebanon.

Methods

WGS using long reads sequencing (PacBio) was performed to produce complete circular genomes and to determine clonality, antimicrobial resistance and virulence determinants.

Results

The sequencing results yielded eight closed circular genomes. Three multilocus sequence typing (MLST) types were identified (ST11618, ST14184, ST15253). Both MLST and WGS analyses revealed that these isolates from Lebanon were genetically homogenous belonging to clonal complex CC230 and clustered closely with isolates originating from Canada, United States of America, United Kingdom and Iceland. Their penicillin binding protein profiles correlated with both β-lactam susceptibility patterns and MLST types. Moreover, the isolates harbored the macrolide and tetracycline resistance genes and showed a similar virulence gene profile. To our knowledge, this study represents the first report of complete phenotypic and genomic characterization of the emerging Streptococcus pneumoniae, serotype 24F, in the Middle East and North Africa region.

Post-Vaccination Streptococcus pneumoniae Carriage and Virulence Gene Distribution among Children Less Than Five Years of Age, Cape Coast, Ghana

In 2012, Ghana introduced PCV13 into its childhood immunization program. To monitor the pneumococcus after PCV13 vaccination, we analyzed serotypes, antibiotic resistance, and virulence genes of pneumococcal carriage isolates among children under five years of age. We obtained nasopharyngeal swabs from 513 children from kindergartens and immunization centers in Cape Coast, Ghana. Pneumococcal serotypes were determined by multiplex-PCR and Quellung reaction. Antibiotic resistance and virulence genes prevalence were determined by disc diffusion and PCR respectively. Overall, carriage prevalence was 29.4% and PCV13 coverage was 38.4%. Over 60% of the isolates were non-PCV13 serotypes and serotype 23B was the most prevalent. One isolate showed full resistance to penicillin, while 35% showed intermediate resistance. Resistance to erythromycin and clindamycin remained low, while susceptibility to ceftriaxone, levofloxacin and vancomycin remained high. Penicillin resistance was associated with PCV13 serotypes. Forty-three (28.5%) strains were multidrug-resistant. Virulence genes pavB, pcpA, psrP, pilus-1, and pilus-2 were detected in 100%, 87%, 62.9%, 11.9%, and 6.6% of the strains, respectively. The pilus islets were associated with PCV13 and multidrug-resistant serotypes. PCV13 vaccination had impacted on pneumococcal carriage with a significant increase in non-PCV13 serotypes and lower penicillin resistance. Including PcpA and PsrP in pneumococcal protein-based vaccines could be beneficial to Ghanaian children.

A Choline-Recognizing Monomeric Lysin, ClyJ-3m, Shows Elevated Activity against Streptococcus pneumoniae

Streptococcus pneumoniae is a leading pathogen for bacterial pneumonia, which can be treated with bacteriophage lysins harboring a conserved choline binding module (CBM). Such lysins regularly function as choline-recognizing dimers. Previously, we reported a pneumococcus-specific lysin ClyJ comprising the binding domain from the putative endolysin gp20 from the Streptococcus phage SPSL1 and the CHAP (cysteine, histidine-dependent amidohydrolase/peptidase) catalytic domain from the PlyC lysin. A variant of ClyJ with a shortened linker, i.e., ClyJ-3, shows improved activity and reduced cytotoxicity. Resembling typical CBM-containing lysins, ClyJ-3 dimerized upon binding with choline. Herein, we further report a choline-recognizing variant of ClyJ-3, i.e., ClyJ-3m, constructed by deleting its C-terminal tail. Biochemical characterization showed that ClyJ-3m remains a monomer after it binds to choline yet exhibits improved bactericidal activity against multiple pneumococcal strains with different serotypes. In an S. pneumoniae-infected bacteremia model, a single intraperitoneal administration of 2.32 μg/mouse of ClyJ-3m showed 70% protection, while only 20% of mice survived in the group receiving an equal dose of ClyJ-3 (P < 0.05). A pharmacokinetic analysis following single intravenously doses of 0.29 and 1.16 mg/kg of ClyJ-3 or ClyJ-3m in BALB/c mice revealed that ClyJ-3m shows a similar half-life but less clearance and a greater area under curve than ClyJ-3. Taken together, the choline-recognizing monomer ClyJ-3m exhibited enhanced bactericidal activity and improved pharmacokinetic proprieties compared to those of its parental ClyJ-3 lysin. Our study also provides a new way for rational design and programmed engineering of lysins targeting S. pneumoniae.

Improved Prediction of Bacterial Genotype-Phenotype Associations Using Interpretable Pangenome-Spanning Regressions

Discovery of genetic variants underlying bacterial phenotypes and the prediction of phenotypes such as antibiotic resistance are fundamental tasks in bacterial genomics. Genome-wide association study (GWAS) methods have been applied to study these relations, but the plastic nature of bacterial genomes and the clonal structure of bacterial populations creates challenges. We introduce an alignment-free method which finds sets of loci associated with bacterial phenotypes, quantifies the total effect of genetics on the phenotype, and allows accurate phenotype prediction, all within a single computationally scalable joint modeling framework. Genetic variants covering the entire pangenome are compactly represented by extended DNA sequence words known as unitigs, and model fitting is achieved using elastic net penalization, an extension of standard multiple regression. Using an extensive set of state-of-the-art bacterial population genomic data sets, we demonstrate that our approach performs accurate phenotype prediction, comparable to popular machine learning methods, while retaining both interpretability and computational efficiency. Compared to those of previous approaches, which test each genotype-phenotype association separately for each variant and apply a significance threshold, the variants selected by our joint modeling approach overlap substantially.IMPORTANCE Being able to identify the genetic variants responsible for specific bacterial phenotypes has been the goal of bacterial genetics since its inception and is fundamental to our current level of understanding of bacteria. This identification has been based primarily on painstaking experimentation, but the availability of large data sets of whole genomes with associated phenotype metadata promises to revolutionize this approach, not least for important clinical phenotypes that are not amenable to laboratory analysis. These models of phenotype-genotype association can in the future be used for rapid prediction of clinically important phenotypes such as antibiotic resistance and virulence by rapid-turnaround or point-of-care tests. However, despite much effort being put into adapting genome-wide association study (GWAS) approaches to cope with bacterium-specific problems, such as strong population structure and horizontal gene exchange, current approaches are not yet optimal. We describe a method that advances methodology for both association and generation of portable prediction models.

Restricted Sequence Variation in Streptococcus pyogenes Penicillin Binding Proteins

β-Lactam antibiotics are the first-line therapeutic option for Streptococcus pyogenes infections. Despite the global high prevalence of S. pyogenes infections and widespread use of β-lactams worldwide, reports of resistance to β-lactam antibiotics, such as penicillin, have been incredibly rare. Recently, β-lactam resistance, as defined by clinical breakpoints, was detected in two clinical S. pyogenes isolates with accompanying mutations in the active site of the penicillin binding protein PBP2x, raising concerns that β-lactam resistance will become more widespread. We screened a global database of S. pyogenes genome sequences to investigate the frequency of PBP mutations, identifying that PBP mutations are uncommon relative to those of Streptococcus pneumoniae. These findings support clinical observations that β-lactam resistance is rare in S. pyogenes and suggest that there are considerable constraints on S. pyogenes PBP sequence variation.

A recent clinical report has linked Streptococcus pyogenes β-lactam antibiotic resistance to mutation in the penicillin binding protein (PBP) PBP2x. To determine whether this is an isolated case or reflects a broader prevalence of mutations that might confer reduced β-lactam susceptibility, we investigated the relative frequency of PBP sequence variation within a global database of 9,667 S. pyogenes isolates. We found that mutations in S. pyogenes PBPs (PBP2x, PBP1a, PBP1b, and PBP2a) occur infrequently across this global database, with fewer than 3 amino acid changes differing between >99% of the global population. Only 4 of the 9,667 strains contained mutations near transpeptidase active sites of PBP2x or PBP1a. The reported PBP2x T553K substitution was not identified. These findings are in contrast to those of 2,520 S. pneumococcus sequences where PBP mutations are relatively frequent and are often located in key β-lactam binding pockets. These data, combined with the general lack of penicillin resistance reported in S. pyogenes worldwide, suggests that extensive, unknown constraints restrict S. pyogenes PBP sequence plasticity. Our findings imply that while heavy antibiotic pressure may select for mutations in the PBPs, there is currently no evidence of such mutations becoming fixed in the S. pyogenes population or that mutations are being sequentially acquired in the PBPs.

IMPORTANCE β-Lactam antibiotics are the first-line therapeutic option for Streptococcus pyogenes infections. Despite the global high prevalence of S. pyogenes infections and widespread use of β-lactams worldwide, reports of resistance to β-lactam antibiotics, such as penicillin, have been incredibly rare. Recently, β-lactam resistance, as defined by clinical breakpoints, was detected in two clinical S. pyogenes isolates with accompanying mutations in the active site of the penicillin binding protein PBP2x, raising concerns that β-lactam resistance will become more widespread. We screened a global database of S. pyogenes genome sequences to investigate the frequency of PBP mutations, identifying that PBP mutations are uncommon relative to those of Streptococcus pneumoniae. These findings support clinical observations that β-lactam resistance is rare in S. pyogenes and suggest that there are considerable constraints on S. pyogenes PBP sequence variation.