Linear regression equations to predict β-lactam, macrolide, lincosamide and fluoroquinolone minimum inhibitory concentrations from molecular antimicrobial resistance determinants in Streptococcus pneumoniae

Holistic understanding of trimethoprim resistance in Streptococcus pneumoniae using an integrative approach of genome-wide association study, resistance reconstruction, and machine learning

Antimicrobial resistance (AMR) is a public health threat worldwide. Next-generation sequencing (NGS) has opened unprecedented opportunities to accelerate AMR mechanism discovery and diagnostics. Here, we present an integrative approach to investigate trimethoprim (TMP) resistance in the key pathogen Streptococcus pneumoniae. We explored a collection of 662 S. pneumoniae genomes by conducting a genome-wide association study (GWAS), followed by functional validation using resistance reconstruction experiments, combined with machine learning (ML) approaches to predict TMP minimum inhibitory concentration (MIC). Our study showed that multiple additive mutations in the folA and sulA loci are responsible for TMP non-susceptibility in S. pneumoniae and can be used as key features to build ML models for digital MIC prediction, reaching an average accuracy within ±1 twofold dilution factor of 86.3%. Our roadmap of in silico analysis-wet-lab validation-diagnostic tool building could be adapted to explore AMR in other combinations of bacteria-antibiotic.

IMPORTANCE

In the age of next-generation sequencing (NGS), while data-driven methods such as genome-wide association study (GWAS) and machine learning (ML) excel at finding patterns, functional validation can be challenging due to the high numbers of candidate variants. We designed an integrative approach combining a GWAS on S. pneumoniae clinical isolates, followed by whole-genome transformation coupled with NGS to functionally characterize a large set of GWAS candidates. Our study validated several phenotypic folA mutations beyond the standard Ile100Leu mutation, and showed that the overexpression of the sulA locus produces trimethoprim (TMP) resistance in Streptococcus pneumoniae. These validated loci, when used to build ML models, were found to be the best inputs for predicting TMP minimal inhibitory concentrations. Integrative approaches can bridge the genotype-phenotype gap by biological insights that can be incorporated in ML models for accurate prediction of drug susceptibility.

Phenotypic and molecular characterization of Streptococcus pneumoniae serotype 3 isolates from blood and respiratory samples in Canada: CANWARD 2007-21

BACKGROUND

Lower respiratory infections and invasive disease caused by Streptococcus pneumoniae serotype 3 remain major clinical challenges around the world, despite widespread availability of updated vaccines.

METHODS

As part of CANWARD, antimicrobial susceptibility testing and serotyping were performed on all S. pneumoniae isolates from 2007 to 2021. A subset of 226/264 (85.6%) serotype 3 isolates were selected for WGS to determine sequence type (ST)/clonal cluster (CC) and correspondence of antimicrobial resistance determinants (erm, mefAE, tetM, cat, folA, folP) with resistance phenotype.

RESULTS

Of the 3,039 S. pneumoniae isolates obtained from 2007 to 2021, 8.7% (n = 264) were serotype 3, with 64.0% of respiratory origin and 36.0% from blood. Of 226 sequenced serotype 3 isolates, 184 (81.4%) were ST180 (GPSC12). The proportion of ST8561 (single locus variant of ST180) increased from 7.2% to 16.6% during the study period. An increasing proportion of serotype 3 isolates had phenotypic resistance (P = 0.0007) and genetic resistance determinants (P = 0.004), comparing 2017-21 to 2007-11, largely due to a recently expanded ST180 clade with cat, tetM and mef determinants.

CONCLUSIONS

S. pneumoniae serotype 3 from GPSC12 continues to dominate throughout Canada, with an increase in the proportion of ST8561. The proportion of serotype 3 isolates that are phenotypically resistant and with genetic resistance determinants is increasing over time, reflecting a global increase in GPSC12 genotypes with known resistance determinants. Phylogenomic characterization of isolates collected over time and from around the world may facilitate improved treatment and enhanced prevention strategies, including new vaccines with activity against S. pneumoniae serotype 3.

Current research update on group B streptococcal infection related to obstetrics and gynecology

Group B streptococcal (GBS) is a Gram-positive bacterium that is commonly found in the gastrointestinal tract and urogenital tract. GBS infestation during pregnancy is a significant contributor to maternal and neonatal morbidity and mortality globally. This article aims to discuss the infectious diseases caused by GBS in the field of obstetrics and gynecology, as well as the challenges associated with the detection, treatment, and prevention of GBS.

Predictive Modeling of Phenotypic Antimicrobial Susceptibility of Selected Beta-Lactam Antimicrobials from Beta-Lactamase Resistance Genes

The outcome of bacterial infection management relies on prompt diagnosis and effective treatment, but conventional antimicrobial susceptibility testing can be slow and labor-intensive. Therefore, this study aims to predict phenotypic antimicrobial susceptibility of selected beta-lactam antimicrobials in the bacteria of the family Enterobacteriaceae from different beta-lactamase resistance genotypes. Using human datasets extracted from the Antimicrobial Testing Leadership and Surveillance (ATLAS) program conducted by Pfizer and retail meat datasets from the National Antimicrobial Resistance Monitoring System for Enteric Bacteria (NARMS), we used a robust or weighted least square multivariable linear regression modeling framework to explore the relationship between antimicrobial susceptibility data of beta-lactam antimicrobials and different types of beta-lactamase resistance genes. In humans, in the presence of the blaCTX-M-1, blaCTX-M-2, blaCTX-M-8/25, and blaCTX-M-9 groups, MICs of cephalosporins significantly increased by values between 0.34-3.07 μg/mL, however, the MICs of carbapenem significantly decreased by values between 0.81-0.87 μg/mL. In the presence of carbapenemase genes (blaKPC, blaNDM, blaIMP, and blaVIM), the MICs of cephalosporin antimicrobials significantly increased by values between 1.06-5.77 μg/mL, while the MICs of carbapenem antimicrobials significantly increased by values between 5.39-67.38 μg/mL. In retail meat, MIC of ceftriaxone increased significantly in the presence of blaCMY-2, blaCTX-M-1, blaCTX-M-55, blaCTX-M-65, and blaSHV-2 by 55.16 μg/mL, 222.70 μg/mL, 250.81 μg/mL, 204.89 μg/mL, and 31.51 μg/mL respectively. MIC of cefoxitin increased significantly in the presence of blaCTX-M-65 and blaTEM-1 by 1.57 μg/mL and 1.04 μg/mL respectively. In the presence of blaCMY-2, MIC of cefoxitin increased by an average of 8.66 μg/mL over 17 years. Compared to E. coli isolates, MIC of cefoxitin in Salmonella enterica isolates decreased significantly by 0.67 μg/mL. On the other hand, MIC of ceftiofur increased in the presence of blaCTX-M-1, blaCTX-M-65, blaSHV-2, and blaTEM-1 by 8.82 μg/mL, 9.11 μg/mL, 8.18 μg/mL, and 1.04 μg/mL respectively. In the presence of blaCMY-2, MIC of ceftiofur increased by an average of 10.20 μg/mL over 14 years. The ability to predict antimicrobial susceptibility of beta-lactam antimicrobials directly from beta-lactamase resistance genes may help reduce the reliance on routine phenotypic testing with higher turnaround times in diagnostic, therapeutic, and surveillance of antimicrobial-resistant bacteria of the family Enterobacteriaceae.

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.

Comparison of contemporary invasive and non-invasive Streptococcus pneumoniae isolates reveals new insights into circulating anti-microbial resistance determinants

Streptococcus pneumoniae is a major human pathogen with a high burden of disease. Non-invasive isolates (those found in non-sterile sites) are thought to be a key source of invasive isolates (those found in sterile sites) and a reservoir of anti-microbial resistance (AMR) determinants. Despite this, pneumococcal surveillance has almost exclusively focused on invasive isolates. We aimed to compare contemporaneous invasive and non-invasive isolate populations to understand how they interact and identify differences in AMR gene distribution. We used a combination of whole-genome sequencing and phenotypic anti-microbial susceptibility testing and a data set of invasive (n = 1,288) and non-invasive (n = 186) pneumococcal isolates, collected in Victoria, Australia, between 2018 and 2022. The non-invasive population had increased levels of antibiotic resistance to multiple classes of antibiotics including beta-lactam antibiotics penicillin and ceftriaxone. We identified genomic intersections between the invasive and non-invasive populations and no distinct phylogenetic clustering of the two populations. However, this analysis revealed sub-populations overrepresented in each population. The sub-populations that had high levels of AMR were overrepresented in the non-invasive population. We determined that WamR-Pneumo was the most accurate in silico tool for predicting resistance to the antibiotics tested. This tool was then used to assess the allelic diversity of the penicillin-binding protein genes, which acquire mutations leading to beta-lactam antibiotic resistance, and found that they were highly conserved (≥80% shared) between the two populations. These findings show the potential of non-invasive isolates to serve as reservoirs of AMR determinants.

Clonal diversity of Streptococcus pneumoniae serotype 19A collected from children < 5 years old in Québec, Canada, 2016-2021

Streptococcus pneumoniae serotype 19A is a highly diverse, often antimicrobial-resistant Gram-positive bacterium which can cause invasive pneumococcal disease (IPD). In 2021, public health authorities in the Canadian province of Québec observed an increase of serotype 19A IPD in children <5 years. The purpose of this study was to determine the clonal composition of serotype 19A isolates collected from this age group in Québec, from 2016 to 2021. Forty-one and 37 IPD isolates from children <5 years from Québec and the remainder of Canada, respectively, were sequenced using the Illumina NextSeq platform. Phylogenetic analysis using SNVPhyl identified three clusters, corresponding to three common clones of serotype 19A: CC199, CC320 and ST695. CC199, predominantly represented by ST416, accounted for similar proportions of serotype 19A isolates collected from children in Québec (19.5 %) and other Canadian jurisdictions (OCJs, 21.6 %), with significant presence of ermB (62.5 % and 60 % of ST416 isolates, respectively). CC320 was more commonly identified from OCJs in comparison to Québec (18.9 % vs. 7.3 %, respectively), but were highly antimicrobial-resistant regardless of region. ST695 was the most common clone of serotype 19A collected in Québec from children <5 years, representing 65.9 % of isolates collected over the study period (40.5 % of isolates collected in OCJs). Phylogenetic analysis identified geographical differences in ST695 across Canada; including a large clade specific to Québec (with both susceptible and macrolide-resistant [ermB] subclades), and a separate macrolide-resistant (mefA) clade associated with OCJs. The Québec-specific ermB-ST695 clone represented 48.1 % of ST695 collected from the province. Continued genomic surveillance of S. pneumoniae serotype 19A is required to: i) track the prevalence and clonal composition of serotype 19A in Québec in future years; ii) characterize the clonal distribution of serotype 19A in adult populations; and iii) monitor whether the currently geographically restricted ermB-ST695 clone observed in Québec expands to OCJs.

An accurate and interpretable model for antimicrobial resistance in pathogenic Escherichia coli from livestock and companion animal species

Understanding the microbial genomic contributors to antimicrobial resistance (AMR) is essential for early detection of emerging AMR infections, a pressing global health threat in human and veterinary medicine. Here we used whole genome sequencing and antibiotic susceptibility test data from 980 disease causing Escherichia coli isolated from companion and farm animals to model AMR genotypes and phenotypes for 24 antibiotics. We determined the strength of genotype-to-phenotype relationships for 197 AMR genes with elastic net logistic regression. Model predictors were designed to evaluate different potential modes of AMR genotype translation into resistance phenotypes. Our results show a model that considers the presence of individual AMR genes and total number of AMR genes present from a set of genes known to confer resistance was able to accurately predict isolate resistance on average (mean F1 score = 98.0%, SD = 2.3%, mean accuracy = 98.2%, SD = 2.7%). However, fitted models sometimes varied for antibiotics in the same class and for the same antibiotic across animal hosts, suggesting heterogeneity in the genetic determinants of AMR resistance. We conclude that an interpretable AMR prediction model can be used to accurately predict resistance phenotypes across multiple host species and reveal testable hypotheses about how the mechanism of resistance may vary across antibiotics within the same class and across animal hosts for the same antibiotic.

Population structure, serotype distribution and antibiotic resistance of Streptococcus pneumoniae causing invasive disease in Victoria, Australia

Streptococcus pneumoniae is a major human pathogen and can cause a range of conditions from asymptomatic colonization to invasive pneumococcal disease (IPD). The epidemiology and distribution of IPD-causing serotypes in Australia has undergone large changes following the introduction of the 7-valent pneumococcal conjugate vaccine (PCV) in 2005 and the 13-valent PCV in 2011. In this study, to provide a contemporary understanding of the IPD causing population in Victoria, Australia, we aimed to examine the population structure and prevalence of antimicrobial resistance using whole-genome sequencing and comprehensive antimicrobial susceptibility data of 1288 isolates collected between 2018 and 2022. We observed high diversity among the isolates with 52 serotypes, 203 sequence types (STs) and 70 Global Pneumococcal Sequencing Project Clusters (GPSCs) identified. Serotypes contained in the 13v-PCV represented 35.3 % (n=405) of isolates. Antimicrobial resistance (AMR) to at least one antibiotic was identified in 23.8 % (n=358) of isolates with penicillin resistance the most prevalent (20.3 %, n=261 using meningitis breakpoints and 5.1 % n=65 using oral breakpoints). Of the AMR isolates, 28 % (n=101) were multidrug resistant (MDR) (resistant to three or more drug classes). Vaccination status of cases was determined for a subset of isolates with 34 cases classified as vaccine failure events (fully vaccinated IPD cases of vaccine serotype). However, no phylogenetic association with failure events was observed. Within the highly diverse IPD population, we identified six high-risk sub-populations of public health concern characterized by high prevalence, high rates of AMR and MDR, or serotype inclusion in vaccines. High-risk serotypes included serotypes 3, 19F, 19A, 14, 11A, 15A and serofamily 23. In addition, we present our data validating seroBA for in silico serotyping to facilitate ISO-accreditation of this test in routine use in a public health reference laboratory and have made this data set available. This study provides insights into the population dynamics, highlights non-vaccine serotypes of concern that are highly resistant, and provides a genomic framework for the ongoing surveillance of IPD in Australia which can inform next-generation IPD prevention strategies.

Genomic investigation of the most common Streptococcus pneumoniae serotypes causing invasive infections in Canada: the SAVE study, 2011-2020

OBJECTIVES

To investigate the lineages and genomic antimicrobial resistance (AMR) determinants of the 10 most common pneumococcal serotypes identified in Canada during the five most recent years of the SAVE study, in the context of the 10-year post-PCV13 period in Canada.

METHODS

The 10 most common invasive Streptococcus pneumoniae serotypes collected by the SAVE study from 2016 to 2020 were 3, 22F, 9N, 8, 4, 12F, 19A, 33F, 23A and 15A. A random sample comprising ∼5% of each of these serotypes collected during each year of the full SAVE study (2011-2020) were selected for whole-genome sequencing (WGS) using the Illumina NextSeq platform. Phylogenomic analysis was performed using the SNVPhyl pipeline. WGS data were used to identify virulence genes of interest, sequence types, global pneumococcal sequence clusters (GPSC) and AMR determinants.

RESULTS

Of the 10 serotypes analysed in this study, six increased significantly in prevalence from 2011 to 2020: 3, 4, 8, 9N, 23A and 33F (P ≤ 0.0201). Serotypes 12F and 15A remained stable in prevalence over time, while serotype 19A decreased in prevalence (P < 0.0001). The investigated serotypes represented four of the most prevalent international lineages causing non-vaccine serotype pneumococcal disease in the PCV13 era: GPSC3 (serotypes 8/33F), GPSC19 (22F), GPSC5 (23A) and GPSC26 (12F). Of these lineages, GPSC5 isolates were found to consistently possess the most AMR determinants. Commonly collected vaccine serotypes 3 and 4 were associated with GPSC12 and GPSC27, respectively. However, a more recently collected lineage of serotype 4 (GPSC192) was highly clonal and possessed AMR determinants.

CONCLUSIONS

Continued genomic surveillance of S. pneumoniae in Canada is essential to monitor for the appearance of new and evolving lineages, including antimicrobial-resistant GPSC5 and GPSC162.

Streptococcus pneumoniae serotyping and antimicrobial susceptibility: assessment for vaccine efficacy in Canada after the introduction of PCV13

BACKGROUND

Streptococcus pneumoniae continues to be an important bacterial pathogen associated with invasive (e.g. bacteraemia, meningitis) and non-invasive (e.g. community-acquired respiratory tract) infections worldwide. Surveillance studies conducted nationally and globally assist in determining trends over geographical areas and allow comparisons between countries.

OBJECTIVES

To characterize invasive isolates of S. pneumoniae in terms of their serotype, antimicrobial resistance, genotype and virulence and to use the serotype data to determine the level of coverage by different generations of pneumococcal vaccines.

METHODS

SAVE (Streptococcus pneumoniae Serotyping and Antimicrobial Susceptibility: Assessment for Vaccine Efficacy in Canada) is an ongoing, annual, national collaborative study between the Canadian Antimicrobial Resistance Alliance (CARE) and the National Microbiology Laboratory, focused on characterizing invasive isolates of S. pneumoniae obtained across Canada. Clinical isolates from normally sterile sites were forwarded by participating hospital public health laboratories to the Public Health Agency of Canada-National Microbiology Laboratory and CARE for centralized phenotypic and genotypic investigation.

RESULTS

The four articles in this Supplement provide a comprehensive examination of the changing patterns of antimicrobial resistance and MDR, serotype distribution, genotypic relatedness and virulence of invasive S. pneumoniae obtained across Canada over a 10 year period (2011-2020).

CONCLUSIONS

The data highlight the evolution of S. pneumoniae under pressure by vaccination and antimicrobial usage, as well as vaccine coverage, allowing both clinicians and researchers nationally and globally to view the current status of invasive pneumococcal infections in Canada.