Fast and flexible bacterial genomic epidemiology with PopPUNK

Streptococcus pneumoniae carriage, serotypes, genotypes, and antimicrobial resistance trends among children in Portugal, after introduction of PCV13 in National Immunization Program: A cross-sectional study

Streptococcus pneumoniae carriage studies are crucial to monitor changes induced by use of pneumococcal conjugate vaccines and inform vaccination policies. In this cross-sectional study, we examined changes within the pneumococcal population following introduction of PCV13 in 2015 in the National Immunization Program (NIP), in Portugal. In 2018-2020 (NIP-PCV13), we obtained 1450 nasopharyngeal samples from children ≤6 years attending day-care. We assessed serotypes, antimicrobial resistance, and genotypes (MLST and GPSC) and compared findings with earlier periods: 2009-2010 (pre-PCV13), 2011-2012 (early-PCV13), and 2015-2016 (late-PCV13). Pneumococcal carriage prevalence remained stable at 60.2 %. Carriage of PCV13 serotypes was 10.7 %, markedly reduced compared to pre-PCV13 period (47.6 %). The most prevalent PCV13 serotypes were 19F, 3, and 19A all showing a significant decreasing trend compared to the pre-PCV13 period (from 7.1 % to 4.7 %, 10.1 % to 1.8 %, and 14.1 % to 1.8 %, respectively), a notable observation given the described limited effectiveness of PCV13 against serotype 3. Non-vaccinated children and children aged 4-6 years were more likely to carry PCV13 serotypes (2.5-fold, 95 %CI [1.1-5.6], and 2.9-fold, 95 %CI [1.3-6.8], respectively). The most prevalent non-PCV13 serotypes were 15B/C, 11A, 23B, 23A, and NT, collectively accounting for 51.9 % of all isolates. In total, 30.5 % of all pneumococci were potentially covered by PCV20. Resistance to penicillin (low-level) and macrolides increased significantly, from 9.3 % and 13.4 %, respectively, in the late-PCV13 period, to approximately 20 % each, mostly due to lineages expressing non-PCV13 serotypes, nearing pre-PCV13 levels. An expansion of lineages traditionally associated with PCV13 serotypes, like CC156-GPSC6 (serotype 14) and CC193-GPSC11 (serotype 19F), but now predominantly expressing non-PCV13 serotypes (11A, 15B/C, and 24F for GPSC6; and 15A and 21 for GPSC11) was noted. These findings indicate that the pneumococcal population is adapting to the pressures conferred by PCV13 and antimicrobial use and indicate the need to maintain close surveillance.

Population genomics of Streptococcus mitis in UK and Ireland bloodstream infection and infective endocarditis cases

Streptococcus mitis is a leading cause of infective endocarditis (IE). However, our understanding of the genomic epidemiology and pathogenicity of IE-associated S. mitis is hampered by low IE incidence. Here we use whole genome sequencing of 129 S. mitis bloodstream infection (BSI) isolates collected between 2001-2016 from clinically diagnosed IE cases in the UK to investigate genetic diversity, antimicrobial resistance, and pathogenicity. We show high genetic diversity of IE-associated S. mitis with virtually all isolates belonging to distinct lineages indicating no predominance of specific lineages. Additionally, we find a highly variable distribution of known pneumococcal virulence genes among the isolates, some of which are overrepresented in disease when compared to carriage strains. Our findings suggest that S. mitis in patients with clinically diagnosed IE is not primarily caused by specific hypervirulent or antimicrobial resistant lineages, highlighting the accidental pathogenic nature of S. mitis in patients with clinically diagnosed IE.

Prediction of post-PCV13 pneumococcal evolution using invasive disease data enhanced by inverse-invasiveness weighting

After introducing pneumococcal conjugate vaccines (PCVs), serotype replacement occurred in Streptococcus pneumoniae. Predicting which pneumococcal strains will become common in carriage after vaccination can enhance vaccine design, public health interventions, and understanding of pneumococcal evolution. Invasive pneumococcal isolates were collected during 1998-2018 by the Active Bacterial Core surveillance (ABCs). Carriage data from Massachusetts (MA) and Southwest United States were used to calculate weights. Using pre-vaccine data, serotype-specific inverse-invasiveness weights were defined as the ratio of the proportion of the serotype in carriage to the proportion in invasive data. Genomic data were processed under bioinformatic pipelines to define genetically similar sequence clusters (i.e., strains), and accessory genes (COGs) present in 5-95% of isolates. Weights were applied to adjust observed strain proportions and COG frequencies. The negative frequency-dependent selection (NFDS) model predicted strain proportions by calculating the post-vaccine strain composition in the weighted invasive disease population that would best match pre-vaccine COG frequencies. Inverse-invasiveness weighting increased the correlation of COG frequencies between invasive and carriage data in linear or logit scale for pre-vaccine, post-PCV7, and post-PCV13; and between different epochs in the invasive data. Weighting the invasive data significantly improved the NFDS model's accuracy in predicting strain proportions in the carriage population in the post-PCV13 epoch, with the adjusted R2 increasing from 0.254 before weighting to 0.545 after weighting. The weighting system adjusted invasive disease data to better represent the pneumococcal carriage population, allowing the NFDS mechanism to predict strain proportions in carriage in the post-PCV13 epoch. Our methods enrich the value of genomic sequences from invasive disease surveillance.IMPORTANCEStreptococcus pneumoniae, a common colonizer in the human nasopharynx, can cause invasive diseases including pneumonia, bacteremia, and meningitis mostly in children under 5 years or older adults. The PCV7 was introduced in 2000 in the United States within the pediatric population to prevent disease and reduce deaths, followed by PCV13 in 2010, PCV15 in 2022, and PCV20 in 2023. After the removal of vaccine serotypes, the prevalence of carriage remained stable as the vacated pediatric ecological niche was filled with certain non-vaccine serotypes. Predicting which pneumococcal clones, and which serotypes, will be most successful in colonization after vaccination can enhance vaccine design and public health interventions, while also improving our understanding of pneumococcal evolution. While carriage data, which are collected from the pneumococcal population that is competing to colonize and transmit, are most directly relevant to evolutionary studies, invasive disease data are often more plentiful. Previously, evolutionary models based on negative frequency-dependent selection (NFDS) on the accessory genome were shown to predict which non-vaccine strains and serotypes were most successful in colonization following the introduction of PCV7. Here, we show that an inverse-invasiveness weighting system applied to invasive disease surveillance data allows the NFDS model to predict strain proportions in the projected carriage population in the post-PCV13/pre-PCV15 and pre-PCV20 epoch. The significance of our research lies in using a sample of invasive disease surveillance data to extend the use of NFDS as an evolutionary mechanism to predict post-PCV13 population dynamics. This has shown that we can correct for biased sampling that arises from differences in virulence and can enrich the value of genomic data from disease surveillance and advance our understanding of how NFDS impacts carriage population dynamics after both PCV7 and PCV13 vaccination.

Mobile genetic elements define the non-random structure of the Salmonella enterica serovar Typhi pangenome

Bacterial relatedness measured using select chromosomal loci forms the basis of public health genomic surveillance. While approximating vertical evolution through this approach has proven exceptionally valuable for understanding pathogen dynamics, it excludes a fundamental dimension of bacterial evolution-horizontal gene transfer. Incorporating the accessory genome is the logical remediation and has recently shown promise in expanding epidemiological resolution for enteric pathogens. Employing k-mer-based Jaccard index analysis, and a novel genome length distance metric, we computed pangenome (i.e., core and accessory) relatedness for the globally important pathogen Salmonella enterica serotype Typhi (Typhi), and graphically express both vertical (homology-by-descent) and horizontal (homology-by-admixture) evolutionary relationships in a reticulate network of over 2,200 U.S. Typhi genomes. This analysis revealed non-random structure in the Typhi pangenome that is driven predominantly by the gain and loss of mobile genetic elements, confirming and expanding upon known epidemiological patterns, revealing novel plasmid dynamics, and identifying avenues for further genomic epidemiological exploration. With an eye to public health application, this work adds important biological context to the rapidly improving ways of analyzing bacterial genetic data and demonstrates the value of the accessory genome to infer pathogen epidemiology and evolution.IMPORTANCEGiven bacterial evolution occurs in both vertical and horizontal dimensions, inclusion of both core and accessory genetic material (i.e., the pangenome) is a logical step toward a more thorough understanding of pathogen dynamics. With an eye to public, and indeed, global health relevance, we couple contemporary tools for genomic analysis with decades of research on mobile genetic elements to demonstrate the value of the pangenome, known and unknown, annotated, and hypothetical, for stratification of Salmonella enterica serovar Typhi (Typhi) populations. We confirm and expand upon what is known about Typhi epidemiology, plasmids, and antimicrobial resistance dynamics, and offer new avenues of exploration to further deduce Typhi ecology and evolution, and ultimately to reduce the incidence of human disease.

Pangenome-spanning epistasis and coselection analysis via de Bruijn graphs

Studies of bacterial adaptation and evolution are hampered by the difficulty of measuring traits such as virulence, drug resistance, and transmissibility in large populations. In contrast, it is now feasible to obtain high-quality complete assemblies of many bacterial genomes thanks to scalable high-accuracy long-read sequencing technologies. To exploit this opportunity, we introduce a phenotype- and alignment-free method for discovering coselected and epistatically interacting genomic variation from genome assemblies covering both core and accessory parts of genomes. Our approach uses a compact colored de Bruijn graph to approximate the intragenome distances between pairs of loci for a collection of bacterial genomes to account for the impacts of linkage disequilibrium (LD). We demonstrate the versatility of our approach to efficiently identify associations between loci linked with drug resistance and adaptation to the hospital niche in the major human bacterial pathogens Streptococcus pneumoniae and Enterococcus faecalis.

Pan-pathogen deep sequencing of nosocomial bacterial pathogens in Italy in spring 2020: a prospective cohort study

BACKGROUND

Nosocomial infections pose a considerable risk to patients who are susceptible, and this is particularly acute in intensive care units when hospital-associated bacteria are endemic. During the first wave of the COVID-19 pandemic, the surge of patients presented a significant obstacle to the effectiveness of infection control measures. We aimed to assess the risks and extent of nosocomial pathogen transmission under a high patient burden by designing a novel bacterial pan-pathogen deep-sequencing approach that could be integrated with standard clinical surveillance and diagnostics workflows.

METHODS

We did a prospective cohort study in a region of northern Italy that was severely affected by the first wave of the COVID-19 pandemic. Inpatients on both ordinary and intensive care unit (ICU) wards at the San Matteo hospital, Pavia were sampled on multiple occasions to identify bacterial pathogens from respiratory, nasal, and rectal samples. Diagnostic samples collected between April 7 and May 10, 2020 were cultured on six different selective media designed to enrich for Acinetobacter baumannii, Escherichia coli, Enterococcus faecium, Enterococcus faecalis, Klebsiella pneumoniae, Pseudomonas aeruginosa, Staphylococcus aureus, and Streptococcus pneumoniae, and DNA from each plate with positive growth was deep sequenced en masse. We used mSWEEP and mGEMS to bin sequencing reads by sequence cluster for each species, followed by mapping with snippy to generate high quality alignments. Antimicrobial resistance genes were detected by use of ARIBA and CARD. Estimates of hospital transmission were obtained from pairwise bacterial single nucleotide polymorphism distances, partitioned by within-patient and between-patient samples. Finally, we compared the accuracy of our binned Acinetobacter baumannii genomes with those obtained by single colony whole-genome sequencing of isolates from the same hospital.

FINDINGS

We recruited patients from March 1 to May 7, 2020. The pathogen population among the patients was large and diverse, with 2148 species detections overall among the 2418 sequenced samples from the 256 patients. In total, 55 sequence clusters from key pathogen species were detected at least five times. The antimicrobial resistance gene prevalence was correspondingly high, with key carbapenemase and extended spectrum ß-lactamase genes detected in at least 50 (40%) of 125 patients in ICUs. Using high-resolution mapping to infer transmission, we established that hospital transmission was likely to be a significant mode of acquisition for each of the pathogen species. Finally, comparison with single colony Acinetobacter baumannii genomes showed that the resolution offered by deep sequencing was equivalent to single-colony sequencing, with the additional benefit of detection of co-colonisation of highly similar strains.

INTERPRETATION

Our study shows that a culture-based deep-sequencing approach is a possible route towards improving future pathogen surveillance and infection control at hospitals. Future studies should be designed to directly compare the accuracy, cost, and feasibility of culture-based deep sequencing with single colony whole-genome sequencing on a range of bacterial species.

FUNDING

Wellcome Trust, European Research Council, Academy of Finland Flagship program, Trond Mohn Foundation, and Research Council of Norway.

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.

Comprehensive blueprint of Salmonella genomic plasticity identifies hotspots for pathogenicity genes

Understanding the dynamic evolution of Salmonella is vital for effective bacterial infection management. This study explores the role of the flexible genome, organised in regions of genomic plasticity (RGP), in shaping the pathogenicity of Salmonella lineages. Through comprehensive genomic analysis of 12,244 Salmonella spp. genomes covering 2 species, 6 subspecies, and 46 serovars, we uncover distinct integration patterns of pathogenicity-related gene clusters into RGP, challenging traditional views of gene distribution. These RGP exhibit distinct preferences for specific genomic spots, and the presence or absence of such spots across Salmonella lineages profoundly shapes strain pathogenicity. RGP preferences are guided by conserved flanking genes surrounding integration spots, implicating their involvement in regulatory networks and functional synergies with integrated gene clusters. Additionally, we emphasise the multifaceted contributions of plasmids and prophages to the pathogenicity of diverse Salmonella lineages. Overall, this study provides a comprehensive blueprint of the pathogenicity potential of Salmonella. This unique insight identifies genomic spots in nonpathogenic lineages that hold the potential for harbouring pathogenicity genes, providing a foundation for predicting future adaptations and developing targeted strategies against emerging human pathogenic strains.

Comparison of gene-by-gene and genome-wide short nucleotide sequence-based approaches to define the global population structure of Streptococcus pneumoniae

Defining the population structure of a pathogen is a key part of epidemiology, as genomically related isolates are likely to share key clinical features such as antimicrobial resistance profiles and invasiveness. Multiple different methods are currently used to cluster together closely related genomes, potentially leading to inconsistency between studies. Here, we use a global dataset of 26 306 Streptococcus pneumoniae genomes to compare four clustering methods: gene-by-gene seven-locus MLST, core genome MLST (cgMLST)-based hierarchical clustering (HierCC) assignments, life identification number (LIN) barcoding and k-mer-based PopPUNK clustering (known as GPSCs in this species). We compare the clustering results with phylogenetic and pan-genome analyses to assess their relationship with genome diversity and evolution, as we would expect a good clustering method to form a single monophyletic cluster that has high within-cluster similarity of genomic content. We show that the four methods are generally able to accurately reflect the population structure based on these metrics and that the methods were broadly consistent with each other. We investigated further to study the discrepancies in clusters. The greatest concordance was seen between LIN barcoding and HierCC (adjusted mutual information score=0.950), which was expected given that both methods utilize cgMLST, but have different methods for defining an individual cluster and different core genome schema. However, the existence of differences between the two methods shows that the selection of a core genome schema can introduce inconsistencies between studies. GPSC and HierCC assignments were also highly concordant (AMI=0.946), showing that k-mer-based methods which use the whole genome and do not require the careful selection of a core genome schema are just as effective at representing the population structure. Additionally, where there were differences in clustering between these methods, this could be explained by differences in the accessory genome that were not identified in cgMLST. We conclude that for S. pneumoniae, standardized and stable nomenclature is important as the number of genomes available expands. Furthermore, the research community should transition away from seven-locus MLST, whilst cgMLST, GPSC and LIN assignments should be used more widely. However, to allow for easy comparison between studies and to make previous literature relevant, the reporting of multiple clustering names should be standardized within the research.

Development of the Pneumococcal Genome Library, a core genome multilocus sequence typing scheme, and a taxonomic life identification number barcoding system to investigate and define pneumococcal population structure

Investigating the genomic epidemiology of major bacterial pathogens is integral to understanding transmission, evolution, colonization, disease, antimicrobial resistance and vaccine impact. Furthermore, the recent accumulation of large numbers of whole genome sequences for many bacterial species enhances the development of robust genome-wide typing schemes to define the overall bacterial population structure and lineages within it. Using the previously published data, we developed the Pneumococcal Genome Library (PGL), a curated dataset of 30 976 genomes and contextual data for carriage and disease pneumococci recovered between 1916 and 2018 in 82 countries. We leveraged the size and diversity of the PGL to develop a core genome multilocus sequence typing (cgMLST) scheme comprised of 1222 loci. Finally, using multilevel single-linkage clustering, we stratified pneumococci into hierarchical clusters based on allelic similarity thresholds and defined these with a taxonomic life identification number (LIN) barcoding system. The PGL, cgMLST scheme and LIN barcodes represent a high-quality genomic resource and fine-scale clustering approaches for the analysis of pneumococcal populations, which support the genomic epidemiology and surveillance of this leading global pathogen.

Phylogenomic Analyses of Three Distinct Lineages Uniting Staphylococcus cohnii and Staphylococcus urealyticus from Diverse Hosts

We sequenced and assembled genomes for 17 isolates of Staphylococcus cohnii isolated from osteomyelitis lesions in young broilers from two separate experiments where we induced lameness using a hybrid wire-litter flooring system. Whole genome comparisons using three different methods support a close relationship of genomes from both S. cohnii and Staphylococcus urealyticus. The data support three different lineages, which we designated as Lineage 1, Lineage 2, and Lineage 3, uniting these two species within an evolving complex. We present evidence for horizontal transfer between lineages of genomic regions from 50-440 kbp. The transfer of a 186 kbp region from Lineage 1 to Lineage 2 appears to have generated Lineage 3. Human-associated isolates appear to be limited to Lineages 2 and 3 but Lineage 2 appears to contain a higher number of human pathogenic isolates. The chicken isolates from our lameness trials included genomically diverse isolates from both Lineage 1 and 2, and isolates from both lineages were obtained from osteomyelitis lesions of individual birds. Our results expand the diversity of Staphylococci associated with osteomyelitis in poultry and suggest a high diversity in the microbiome of day-old chicks. Our data also support a reevaluation and unification of the taxonomic classifications of S. cohnii and S. urealyticus.

Clonal Expansion of a Streptococcus pneumoniae Serotype 3 Capsule Variant Sequence Type 700 With Enhanced Vaccine Escape Potential After 13-Valent Pneumococcal Conjugate Vaccine Introduction

BACKGROUND

Streptococcus pneumoniae serotype 3 remains a problem globally. Malawi introduced 13-valent pneumococcal conjugate vaccine (PCV13) in 2011, but there has been no direct protection against serotype 3 carriage. We explored whether vaccine escape by serotype 3 is due to clonal expansion of a lineage with a competitive advantage.

METHODS

The distribution of serotype 3 Global Pneumococcal Sequence Clusters (GPSCs) and sequence types (STs) globally was assessed using sequences from the Global Pneumococcal Sequencing Project. Whole-genome sequences of 135 serotype 3 carriage isolates from Blantyre, Malawi (2015-2019) were analyzed. Comparative analysis of the capsule locus, entire genomes, antimicrobial resistance, and phylogenetic reconstructions were undertaken. Opsonophagocytosis was evaluated using serum samples from vaccinated adults and children.

RESULTS

Serotype 3 GPSC10-ST700 isolates were most prominent in Malawi. Compared with the prototypical serotype 3 capsular polysaccharide locus sequence, 6 genes are absent, with retention of capsule polysaccharide biosynthesis. This lineage is characterized by increased antimicrobial resistance and lower susceptibility to opsonophagocytic killing.

CONCLUSIONS

A serotype 3 variant in Malawi has genotypic and phenotypic characteristics that could enhance vaccine escape and clonal expansion after post-PCV13 introduction. Genomic surveillance among high-burden populations is essential to improve the effectiveness of next-generation pneumococcal vaccines.

Evidence of Reduced Virulence and Increased Colonization Among Pneumococcal Isolates of Serotype 3 Clade II Lineage in Mice

Recent phylogenetic profiling of pneumococcal serotype 3 (Pn3) isolates revealed a dynamic interplay among major lineages with the emergence and global spread of a variant termed clade II. The cause of Pn3 clade II dissemination along with epidemiological and clinical ramifications are currently unknown. Here, we sought to explore biological characteristics of dominant Pn3 clades in a mouse model of pneumococcal invasive disease and carriage. Carriage and virulence potential were strain dependent with marked differences among clades. We found that clinical isolates from Pn3 clade II are less virulent and less invasive in mice compared to clade I isolates. We also observed that clade II isolates are carried for longer and at higher bacterial densities in mice compared to clade I isolates. Taken together, our data suggest that the epidemiological success of Pn3 clade II could be related to alterations in the pathogen's ability to cause invasive disease and to establish a robust carriage episode.

Diversity and phage sensitivity to phages of porcine enterotoxigenic Escherichia coli

Enterotoxigenic Escherichia coli (ETEC) is a diverse and poorly characterized E. coli pathotype that causes diarrhea in humans and animals. Phages have been proposed for the veterinary biocontrol of ETEC, but effective solutions require understanding of porcine ETEC diversity that affects phage infection. Here, we sequenced and analyzed the genomes of the PHAGEBio ETEC collection, gathering 79 diverse ETEC strains isolated from European pigs with post-weaning diarrhea (PWD). We identified the virulence factors characterizing the pathotype and several antibiotic resistance genes on plasmids, while phage resistance genes and other virulence factors were mostly chromosome encoded. We experienced that ETEC strains were highly resistant to Enterobacteriaceae phage infection. It was only by enrichment of numerous diverse samples with different media and conditions, using the 41 ETEC strains of our collection as hosts, that we could isolate two lytic phages that could infect a large part of our diverse ETEC collection: vB_EcoP_ETEP21B and vB_EcoS_ETEP102. Based on genome and host range analyses, we discussed the infection strategies of the two phages and identified components of lipopolysaccharides ( LPS) as receptors for the two phages. Our detailed computational structural analysis highlights several loops and pockets in the tail fibers that may allow recognition and binding of ETEC strains, also in the presence of O-antigens. Despite the importance of receptor recognition, the diversity of the ETEC strains remains a significant challenge for isolating ETEC phages and developing sustainable phage-based products to address ETEC-induced PWD.IMPORTANCEEnterotoxigenic Escherichia coli (ETEC)-induced post-weaning diarrhea is a severe disease in piglets that leads to weight loss and potentially death, with high economic and animal welfare costs worldwide. Phage-based approaches have been proposed, but available data are insufficient to ensure efficacy. Genome analysis of an extensive collection of ETEC strains revealed that phage defense mechanisms were mostly chromosome encoded, suggesting a lower chance of spread and selection by phage exposure. The difficulty in isolating lytic phages and the molecular and structural analyses of two ETEC phages point toward a multifactorial resistance of ETEC to phage infection and the importance of extensive phage screenings specifically against clinically relevant strains. The PHAGEBio ETEC collection and these two phages are valuable tools for the scientific community to expand our knowledge on the most studied, but still enigmatic, bacterial species-E. coli.

Genetic diversity, determinants, and dissemination of Burkholderia pseudomallei lineages implicated in melioidosis in Northeast Thailand

Melioidosis is an often-fatal neglected tropical disease caused by an environmental bacterium Burkholderia pseudomallei. However, our understanding of the disease-causing bacterial lineages, their dissemination, and adaptive mechanisms remains limited. To address this, we conduct a comprehensive genomic analysis of 1,391 B. pseudomallei isolates collected from nine hospitals in northeast Thailand between 2015 and 2018, and contemporaneous isolates from neighbouring countries, representing the most densely sampled collection to date. Our study identifies three dominant lineages, each with unique gene sets potentially enhancing bacterial fitness in the environment. We find that recombination drives lineage-specific gene flow. Transcriptome analyses of representative clinical isolates from each dominant lineage reveal increased expression of lineage-specific genes under environmental conditions in two out of three lineages. This underscores the potential importance of environmental persistence for these dominant lineages. The study also highlights the influence of environmental factors such as terrain slope, altitude, and river direction on the geographical dispersal of B. pseudomallei. Collectively, our findings suggest that environmental persistence may play a role in facilitating the spread of B. pseudomallei, and as a prerequisite for exposure and infection, thereby providing useful insights for informing melioidosis prevention and control strategies.

DODGE: automated point source bacterial outbreak detection using cumulative long term genomic surveillance

SUMMARY

The reliable and timely recognition of outbreaks is a key component of public health surveillance for foodborne diseases. Whole genome sequencing (WGS) offers high resolution typing of foodborne bacterial pathogens and facilitates the accurate detection of outbreaks. This detection relies on grouping WGS data into clusters at an appropriate genetic threshold. However, methods and tools for selecting and adjusting such thresholds according to the required resolution of surveillance and epidemiological context are lacking. Here we present DODGE (Dynamic Outbreak Detection for Genomic Epidemiology), an algorithm to dynamically select and compare these genetic thresholds. DODGE can analyse expanding datasets over time and clusters that are predicted to correspond to outbreaks (or "investigation clusters") can be named with established genomic nomenclature systems to facilitate integrated analysis across jurisdictions. DODGE was tested in two real-world Salmonella genomic surveillance datasets of different duration, 2 months from Australia and 9 years from the United Kingdom. In both cases only a minority of isolates were identified as investigation clusters. Two known outbreaks in the United Kingdom dataset were detected by DODGE and were recognized at an earlier timepoint than the outbreaks were reported. These findings demonstrated the potential of the DODGE approach to improve the effectiveness and timeliness of genomic surveillance for foodborne diseases and the effectiveness of the algorithm developed.

AVAILABILITY AND IMPLEMENTATION

DODGE is freely available at https://github.com/LanLab/dodge and can easily be installed using Conda.

Geographical migration and fitness dynamics of Streptococcus pneumoniae

Streptococcus pneumoniae is a leading cause of pneumonia and meningitis worldwide. Many different serotypes co-circulate endemically in any one location1,2. The extent and mechanisms of spread and vaccine-driven changes in fitness and antimicrobial resistance remain largely unquantified. Here using geolocated genome sequences from South Africa (n = 6,910, collected from 2000 to 2014), we developed models to reconstruct spread, pairing detailed human mobility data and genomic data. Separately, we estimated the population-level changes in fitness of strains that are included (vaccine type (VT)) and not included (non-vaccine type (NVT)) in pneumococcal conjugate vaccines, first implemented in South Africa in 2009. Differences in strain fitness between those that are and are not resistant to penicillin were also evaluated. We found that pneumococci only become homogenously mixed across South Africa after 50 years of transmission, with the slow spread driven by the focal nature of human mobility. Furthermore, in the years following vaccine implementation, the relative fitness of NVT compared with VT strains increased (relative risk of 1.68; 95% confidence interval of 1.59-1.77), with an increasing proportion of these NVT strains becoming resistant to penicillin. Our findings point to highly entrenched, slow transmission and indicate that initial vaccine-linked decreases in antimicrobial resistance may be transient.

RATA: A novel class A carbapenemase with broad geographic distribution and potential for global spread

Carbapenem resistance's global proliferation poses a significant public health challenge. The primary resistance mechanism is carbapenemase production. In this study, we discovered a novel carbapenemase, RATA, located on the chromosome of Riemerella anatipestifer isolates. This enzyme shares ≤52 % amino acid sequence identity with other known β-lactamases. Antimicrobial susceptibility tests and kinetic assays demonstrated that RATA could hydrolyze not only penicillins and extended-spectrum cephalosporins but also monobactams, cephamycins, and carbapenems. Furthermore, its activity was readily inhibited by β-lactamase inhibitors. Bioinformatic analysis revealed 46 blaRATA-like genes encoding 27 variants in the NCBI database, involving 21 different species, including pathogens, host-associated bacteria, and environmental isolates. Notably, blaRATA-positive strains were globally distributed and primarily collected from marine environments. Concurrently, taxonomic analysis and GC content analysis indicated that blaRATA orthologue genes were predominantly located on the chromosomes of Flavobacteriaceae and shared a similar GC content as Flavobacteriaceae. Although no explicit mobile genetic elements were identified by genetic environment analysis, blaRATA-2 possessed the ability of horizontal transfer in R. anatipestifer via natural transformation. This work's data suggest that RATA is a new chromosome-encoded class A carbapenemase, and Flavobacteriaceae from marine environments could be the primary reservoir of the blaRATA gene.

Deep sequencing of Escherichia coli exposes colonisation diversity and impact of antibiotics in Punjab, Pakistan

Multi-drug resistant (MDR) E. coli constitute a major public health burden globally, reaching the highest prevalence in the global south yet frequently flowing with travellers to other regions. However, our comprehension of the entire genetic diversity of E. coli colonising local populations remains limited. We quantified this diversity, its associated antimicrobial resistance (AMR), and assessed the impact of antibiotic use by recruiting 494 outpatients and 423 community dwellers in the Punjab province, Pakistan. Rectal swab and stool samples were cultured on CLED agar and DNA extracted from plate sweeps was sequenced en masse to capture both the genetic and AMR diversity of E. coli. We assembled 5,247 E. coli genomes from 1,411 samples, displaying marked genetic diversity in gut colonisation. Compared with high income countries, the Punjabi population generally showed a markedly different distribution of genetic lineages and AMR determinants, while use of antibiotics elevated the prevalence of well-known globally circulating MDR clinical strains. These findings implicate that longitudinal multi-regional genomics-based surveillance of both colonisation and infections is a prerequisite for developing mechanistic understanding of the interplay between ecology and evolution in the maintenance and dissemination of (MDR) E. coli.

Genomic investigation of Lactococcus formosensis, Lactococcus garvieae, and Lactococcus petauri reveals differences in species distribution by human and animal sources

Lactococcus garvieae is a fish pathogen that can cause diseases in humans and cows. Two genetically related species, Lactococcus formosensis and Lactococcus petauri, may be misidentified as L. garvieae. It is unclear if these species differ in host specificity and virulence genes. This study analyzed the genomes of 120 L. petauri, 53 L. formosensis, and 39 L. garvieae isolates from various sources. The genetic diversity and virulence gene content of these isolates were compared. The results showed that 77 isolates previously reported as L. garvieae were actually L. formosensis or L. petauri. The distribution of the three species varied across different collection sources, with L. petauri being predominant in human infections, human fecal sources, and rainbow trout, while L. formosensis was more common in bovine isolates. The genetic diversity of isolates within each species was high and similar. Using a genomic clustering method, L. petauri, L. formosensis, and L. garvieae were divided into 45, 22, and 13 clusters, respectively. Most rainbow trout and human isolates of L. petauri belonged to different clusters, while L. formosensis isolates from bovine and human sources were also segregated into separate clusters. In L. garvieae, most human isolates were grouped into three clusters that also included isolates from food or other sources. Non-metric multidimensional scaling ordination revealed the differential association of 15 virulence genes, including 14 adherence genes and a bile salt hydrolase gene, with bacterial species and certain collection sources. In conclusion, this work provides evidence of host specificity among the three species.

IMPORTANCE

Lactococcus formosensis and Lactococcus petauri are two newly discovered bacteria, which are closely related to Lactococcus garvieae, a pathogen that affects farmed rainbow trout, as well as causes cow mastitis and human infections. It is unclear whether the three bacteria differ in their host preference and the presence of genes that contribute to the development of disease. This study shows that L. formosensis and L. petauri were commonly misidentified as L. garvieae. The three bacteria showed different distribution patterns across various sources. L. petauri was predominantly found in human infections and rainbow trout, while L. formosensis was more commonly detected in cow mastitis. Fifteen genes displayed a differential distribution among the three bacteria from certain sources, indicating a genetic basis for the observed host preference. This work indicates the importance of differentiating the three bacteria in diagnostic laboratories for surveillance and outbreak investigation purposes.

Estimating geographical spread of Streptococcus pneumoniae within Israel using genomic data

Understanding how pathogens spread across geographical space is fundamental for control measures such as vaccination. Streptococcus pneumoniae (the pneumococcus) is a respiratory bacterium responsible for a large proportion of infectious disease morbidity and mortality globally. Even in the post-vaccination era, the rates of invasive pneumococcal disease (IPD) remain stable in most countries, including Israel. To understand the geographical spread of the pneumococcus in Israel, we analysed 1174 pneumococcal genomes from patients with IPD across multiple regions. We included the evolutionary distance between pairs of isolates inferred using whole-genome data within a relative risk (RR) ratio framework to capture the geographical structure of S. pneumoniae. While we could not find geographical structure at the overall lineage level, the extra granularity provided by whole-genome sequence data showed that it takes approximately 5 years for invasive pneumococcal isolates to become fully mixed across the country.This article contains data hosted by Microreact.

Expansion of pneumococcal serotype 23F and 14 lineages with genotypic changes in capsule polysaccharide locus and virulence gene profiles post introduction of pneumococcal conjugate vaccine in Blantyre, Malawi

Since the introduction of the 13-valent pneumococcal conjugate vaccine (PCV13) in Malawi in 2011, there has been persistent carriage of vaccine serotype (VT) Streptococcus pneumoniae, despite high vaccine coverage. To determine if there has been a genetic change within the VT capsule polysaccharide (cps) loci since the vaccine's introduction, we compared 1022 whole-genome-sequenced VT isolates from 1998 to 2019. We identified the clonal expansion of a multidrug-resistant, penicillin non-susceptible serotype 23F GPSC14-ST2059 lineage, a serotype 14 GPSC9-ST782 lineage and a novel serotype 14 sequence type GPSC9-ST18728 lineage. Serotype 23F GPSC14-ST2059 had an I253T mutation within the capsule oligosaccharide repeat unit polymerase Wzy protein, which is predicted in silico to alter the protein pocket cavity. Moreover, serotype 23F GPSC14-ST2059 had SNPs in the DNA binding sites for the cps transcriptional repressors CspR and SpxR. Serotype 14 GPSC9-ST782 harbours a non-truncated version of the large repetitive protein (Lrp), containing a Cna protein B-type domain which is also present in proteins associated with infection and colonisation. These emergent lineages also harboured genes associated with antibiotic resistance, and the promotion of colonisation and infection which were absent in other lineages of the same serotype. Together these data suggest that in addition to serotype replacement, modifications of the capsule locus associated with changes in virulence factor expression and antibiotic resistance may promote vaccine escape. In summary, the study highlights that the persistence of vaccine serotype carriage despite high vaccine coverage in Malawi may be partly caused by expansion of VT lineages post-PCV13 rollout.

Limited protection of pneumococcal vaccines against emergent Streptococcus pneumoniae serotype 14/ST876 strains

PURPOSE

Streptococcus pneumoniae (Spn) is a major cause of child death. We investigated the epidemiology of S. pneumoniae in a pediatric fever clinic and explored the genomics basis of the limited vaccine response of serotype 14 strains worldwide.

METHODS

Febrile disease and pneumonia were diagnosed following criteria from the WHO at the end of 2019 at a tertiary children's hospital. Spn was isolated by culture from nasopharyngeal (NP) swabs. The density was determined by lytA-base qPCR. Isolates were serotyped by Quellung and underwent antimicrobial susceptibility testing. Whole-genome sequencing was employed for molecular serotyping, MLST, antibiotic gene determination, SNP calling, recombination prediction, and phylogenetic analysis.

RESULTS

The presence of pneumococcus in the nasopharynx (87.5%, 7/8, p = 0.0227) and a high carriage (100%, 7/7, p = 0.0123) were significantly associated with pneumonia development. Living with siblings (73.7%, 14/19, p = 0.0125) and non-vaccination (56.0%, 28/50, p = 0.0377) contributed significantly to the Spn carriage. Serotype 14 was the most prevalent strain (16.67%, 5/30). The genome analysis of 1497 serotype 14 strains indicated S14/ST876 strains were only prevalent in China, presented limited vaccine responses with higher recombination activities within its cps locus, and unique variation patterns in the genes wzg and lrp.

CONCLUSION

With the lifting of the one-child policy, it will be crucial for families with multiple children to get PCV vaccinations in China. Due to the highly variant cps locus and distinctive variation patterns in capsule shedding and binding proteins genes, the prevalent S14/ST876 strains have shown poor response to current vaccines. It is necessary to continue monitoring the molecular epidemiology of this vaccine escape clone.

Complete genome sequence of invasive Streptococcus pneumoniae serotype 11A causing meningitis in an adult patient, Italy 2022

Streptococcus pneumoniae is a major global health concern, being a common cause of meningitis in both children and adults. Here, we report the complete genome sequence of P10_PNE_LCR, a S. pneumoniae 11A strain isolated in Northern Italy from an adult patient diagnosed with meningitis.

Recombinational exchange of M-fibril and T-pilus genes generates extensive cell surface diversity in the global group A Streptococcus population

Among genes present in all group A streptococci (GAS), those encoding M-fibril and T-pilus proteins display the highest levels of sequence diversity, giving rise to the two primary serological typing schemes historically used to define strain. A new genotyping scheme for the pilin adhesin and backbone genes is developed and, when combined with emm typing, provides an account of the global GAS strain population. Cluster analysis based on nucleotide sequence similarity assigns most T-serotypes to discrete pilin backbone sequence clusters, yet the established T-types correspond to only half the clusters. The major pilin adhesin and backbone sequence clusters yield 98 unique combinations, defined as "pilin types." Numerous horizontal transfer events that involve pilin or emm genes generate extensive antigenic and functional diversity on the bacterial cell surface and lead to the emergence of new strains. Inferred pilin genotypes applied to a meta-analysis of global population-based collections of pharyngitis and impetigo isolates reveal highly significant associations between pilin genotypes and GAS infection at distinct ecological niches, consistent with a role for pilin gene products in adaptive evolution. Integration of emm and pilin typing into open-access online tools (pubmlst.org) ensures broad utility for end-users wanting to determine the architecture of M-fibril and T-pilus genes from genome assemblies.IMPORTANCEPrecision in defining the variant forms of infectious agents is critical to understanding their population biology and the epidemiology of associated diseases. Group A Streptococcus (GAS) is a global pathogen that causes a wide range of diseases and displays a highly diverse cell surface due to the antigenic heterogeneity of M-fibril and T-pilus proteins which also act as virulence factors of varied functions. emm genotyping is well-established and highly utilized, but there is no counterpart for pilin genes. A global GAS collection provides the basis for a comprehensive pilin typing scheme, and online tools for determining emm and pilin genotypes are developed. Application of these tools reveals the expansion of structural-functional diversity among GAS via horizontal gene transfer, as evidenced by unique combinations of surface protein genes. Pilin and emm genotype correlations with superficial throat vs skin infection provide new insights on the molecular determinants underlying key ecological and epidemiological trends.

Genome plasticity shapes the ecology and evolution of Phocaeicola dorei and Phocaeicola vulgatus

Phocaeicola dorei and Phocaeicola vulgatus are very common and abundant members of the human gut microbiome and play an important role in the infant gut microbiome. These species are closely related and often confused for one another; yet, their genome comparison, interspecific diversity, and evolutionary relationships have not been studied in detail so far. Here, we perform phylogenetic analysis and comparative genomic analyses of these two Phocaeicola species. We report that P. dorei has a larger genome yet a smaller pan-genome than P. vulgatus. We found that this is likely because P. vulgatus is more plastic than P. dorei, with a larger repertoire of genetic mobile elements and fewer anti-phage defense systems. We also found that P. dorei directly descends from a clade of P. vulgatus¸ and experienced genome expansion through genetic drift and horizontal gene transfer. Overall, P. dorei and P. vulgatus have very different functional and carbohydrate utilisation profiles, hinting at different ecological strategies, yet they present similar antimicrobial resistance profiles.

SaLTy: a novel Staphylococcus aureus Lineage Typer

Staphylococcus aureus asymptomatically colonises 30 % of humans but can also cause a range of diseases, which can be fatal. In 2017 S. aureus was associated with 20 000 deaths in the USA alone. Dividing S. aureus isolates into smaller sub-groups can reveal the emergence of distinct sub-populations with varying potential to cause infections. Despite multiple molecular typing methods categorising such sub-groups, they do not take full advantage of S. aureus genome sequences when describing the fundamental population structure of the species. In this study, we developed Staphylococcus aureus Lineage Typing (SaLTy), which rapidly divides the species into 61 phylogenetically congruent lineages. Alleles of three core genes were identified that uniquely define the 61 lineages and were used for SaLTy typing. SaLTy was validated on 5000 genomes and 99.12 % (4956/5000) of isolates were assigned the correct lineage. We compared SaLTy lineages to previously calculated clonal complexes (CCs) from BIGSdb (n=21 173). SALTy improves on CCs by grouping isolates congruently with phylogenetic structure. SaLTy lineages were further used to describe the carriage of Staphylococcal chromosomal cassette containing mecA (SCCmec) which is carried by methicillin-resistant S. aureus (MRSA). Most lineages had isolates lacking SCCmec and the four largest lineages varied in SCCmec over time. Classifying isolates into SaLTy lineages, which were further SCCmec typed, allowed SaLTy to describe high-level MRSA epidemiology. We provide SaLTy as a simple typing method that defines phylogenetic lineages (https://github.com/LanLab/SaLTy). SaLTy is highly accurate and can quickly analyse large amounts of S. aureus genome data. SaLTy will aid the characterisation of S. aureus populations and ongoing surveillance of sub-groups that threaten human health.

Ethanolamine metabolism through two genetically distinct loci enables Klebsiella pneumoniae to bypass nutritional competition in the gut

Successful microbial colonization of the gastrointestinal (GI) tract hinges on an organism's ability to overcome the intense competition for nutrients in the gut between the host and the resident gut microbiome. Enteric pathogens can exploit ethanolamine (EA) in the gut to bypass nutrient competition. However, Klebsiella pneumoniae (K. pneumoniae) is an asymptomatic gut colonizer and, unlike well-studied enteric pathogens, harbors two genetically distinct ethanolamine utilization (eut) loci. Our investigation uncovered unique roles for each eut locus depending on EA utilization as a carbon or nitrogen source. Murine gut colonization studies demonstrated the necessity of both eut loci in the presence of intact gut microbiota for robust GI colonization by K. pneumoniae. Additionally, while some Escherichia coli gut isolates could metabolize EA, other commensals were incapable, suggesting that EA metabolism likely provides K. pneumoniae a selective advantage in gut colonization. Molecular and bioinformatic analyses unveiled the conservation of two eut loci among K. pneumoniae and a subset of the related taxa in the K. pneumoniae species complex, with the NtrC-RpoN regulatory cascade playing a pivotal role in regulation. These findings identify EA metabolism as a critical driver of K. pneumoniae niche establishment in the gut and propose microbial metabolism as a potential therapeutic avenue to combat K. pneumoniae infections.

Invasive pneumococcal disease 3 years after introduction of a reduced 1 + 1 infant 13-valent pneumococcal conjugate vaccine immunisation schedule in England: a prospective national observational surveillance study

BACKGROUND

The UK transition from a 2 + 1 to a 1 + 1 infant immunisation schedule with the 13-valent pneumococcal conjugate vaccine (PCV13) on Jan 1, 2020, coincided with the start of the COVID-19 pandemic. We describe the epidemiology of invasive pneumococcal disease (IPD) in England over 6 financial years (April 1 to March 31) between 2017-18 and 2022-23.

METHODS

We used prospective national surveillance data, including serotyping and whole-genome sequencing of invasive isolates, to analyse IPD trends in England by age and financial year. We compared breakthrough infections and vaccine failure rates in 2022-23 among children eligible for the 1 + 1 schedule with rates in cohorts of children eligible for the 2 + 1 schedule between 2017-18 and 2019-20. We assessed genomic changes over time by comparing Global Pneumococcal Sequencing Clusters and multilocus sequence types among PCV13 serotypes causing IPD.

FINDINGS

There were 4598 laboratory-confirmed IPD cases in 2022-23, 3025 in 2021-22, 1240 in 2020-21, and 5316 in 2019-20. IPD incidence in 2022-23 was 14% lower than in 2019-20 (incidence rate ratio [IRR] 0·86, 95% CI 0·81-0·91; p<0·001). IPD incidence in 2022-23 compared with 2019-20 was 34% higher in children (aged <15 years) (378 cases vs 292 cases; IRR 1·34, 95% CI 1·08-1·68; p=0·009) and 17% lower in adults (aged 15 years and older; 4220 vs 5024; 0·83, 0·78-0·88; p<0·001). The proportion of PCV13-type IPD increased from 19·4% (95% CI 18·2-20·4; 957 of 4947) in 2019-20 to 29·7% (28·3-31·0; 1283 of 4326) in 2022-23, mainly due to serotype 3, but also serotypes 19F, 19A, and 4, alongside a decrease in non-PCV13 serotypes 8, 12F, and 9N. The increase in IPD incidence due to serotypes 3, 19A, and 19F was driven by clonal expansion of previously circulating strains, whereas serotype 4 expansion was driven by newer strains (ie, sequence types 801 and 15603). Breakthrough infections and vaccine failure rates were similar in children eligible for the 1 + 1 (1·08 per 100 000 person-years) and 2 + 1 (0·76 per 100 000 person-years; IRR 1·42, 95% CI 0·78-2·49; p=0·20) PCV13 schedules.

INTERPRETATION

Overall, IPD incidence in England was lower in 2022-23, 2 years after removal of pandemic restrictions, than in 2019-20. Breakthrough and vaccine failure rates were not significantly different between children who received the 1 + 1 compared with the 2 + 1 PCV13 immunisation schedule. The post-pandemic increase in childhood IPD incidence and especially PCV13-type IPD will require close monitoring.

FUNDING

None.

High-resolution genomics identifies pneumococcal diversity and persistence of vaccine types in children with community-acquired pneumonia in the UK and Ireland

BACKGROUND

Streptococcus pneumoniae is a global cause of community-acquired pneumonia (CAP) and invasive disease in children. The CAP-IT trial (grant No. 13/88/11; https://www.capitstudy.org.uk/ ) collected nasopharyngeal swabs from children discharged from hospitals with clinically diagnosed CAP, and found no differences in pneumococci susceptibility between higher and lower antibiotic doses and shorter and longer durations of oral amoxicillin treatment. Here, we studied in-depth the genomic epidemiology of pneumococcal (vaccine) serotypes and their antibiotic resistance profiles.

METHODS

Three-hundred and ninety pneumococci cultured from 1132 nasopharyngeal swabs from 718 children were whole-genome sequenced (Illumina) and tested for susceptibility to penicillin and amoxicillin. Genome heterogeneity analysis was performed using long-read sequenced isolates (PacBio, n = 10) and publicly available sequences.

RESULTS

Among 390 unique pneumococcal isolates, serotypes 15B/C, 11 A, 15 A and 23B1 were most prevalent (n = 145, 37.2%). PCV13 serotypes 3, 19A, and 19F were also identified (n = 25, 6.4%). STs associated with 19A and 19F demonstrated high genome variability, in contrast to serotype 3 (n = 13, 3.3%) that remained highly stable over a 20-year period. Non-susceptibility to penicillin (n = 61, 15.6%) and amoxicillin (n = 10, 2.6%) was low among the pneumococci analysed here and was independent of treatment dosage and duration. However, all 23B1 isolates (n = 27, 6.9%) were penicillin non-susceptible. This serotype was also identified in ST177, which is historically associated with the PCV13 serotype 19F and penicillin susceptibility, indicating a potential capsule-switch event.

CONCLUSIONS

Our data suggest that amoxicillin use does not drive pneumococcal serotype prevalence among children in the UK, and prompts consideration of PCVs with additional serotype coverage that are likely to further decrease CAP in this target population. Genotype 23B1 represents the convergence of a non-vaccine genotype with penicillin non-susceptibility and might provide a persistence strategy for ST types historically associated with vaccine serotypes. This highlights the need for continued genomic surveillance.

Overlapping Streptococcus pyogenes and Streptococcus dysgalactiae subspecies equisimilis household transmission and mobile genetic element exchange

Streptococcus dysgalactiae subspecies equisimilis (SDSE) and Streptococcus pyogenes share skin and throat niches with extensive genomic homology and horizontal gene transfer (HGT) possibly underlying shared disease phenotypes. It is unknown if cross-species transmission interaction occurs. Here, we conduct a genomic analysis of a longitudinal household survey in remote Australian First Nations communities for patterns of cross-species transmission interaction and HGT. Collected from 4547 person-consultations, we analyse 294 SDSE and 315 S. pyogenes genomes. We find SDSE and S. pyogenes transmission intersects extensively among households and show that patterns of co-occurrence and transmission links are consistent with independent transmission without inter-species interference. We identify at least one of three near-identical cross-species mobile genetic elements (MGEs) carrying antimicrobial resistance or streptodornase virulence genes in 55 (19%) SDSE and 23 (7%) S. pyogenes isolates. These findings demonstrate co-circulation of both pathogens and HGT in communities with a high burden of streptococcal disease, supporting a need to integrate SDSE and S. pyogenes surveillance and control efforts.

Genetic comparison of the Mycoplasma gallisepticum 6/85 vaccine strain and 6/85-like field isolates

Mycoplasma gallisepticum infection in poultry leads to disease and pathology that can reduce producer profits. Live attenuated vaccines are available that can limit or completely prevent the effects of infection. Field isolates that are genetically related to the attenuated vaccine strains have been isolated, raising the question of whether the attenuation of the vaccine strains is limited and can lead the strains to revert to more virulent forms. The 6/85 live attenuated vaccine is derived from a field isolate collected in the United States. Analysis of the genome of sequenced M. gallisepticum strains revealed a cluster of 10 6/85-like strains that group with the 6/85 vaccine strain. Four genomic regions were identified that allowed for strain differentiation. The genetic differences between strains points toward nine of the ten strains most likely being sister strains to the 6/85 vaccine strain. Insufficient differences are present in the tenth strain to make a definitive conclusion. These results suggest that most if not all strains similar to the live attenuated vaccine strain are field isolates of the parent strain used to derive the live attenuated vaccine.

Carbapenem-resistant Klebsiella pneumoniae among hospitalized patients in Cape Town, South Africa: molecular epidemiology and characterization

Background

The molecular epidemiology of carbapenem-resistant Enterobacterales in Cape Town remains largely unknown.

Objectives

This study aimed to describe the molecular epidemiology, resistome, virulome and mobilome of carbapenem-resistant Klebsiella pneumoniae (CRKP) within Cape Town to guide therapy, antimicrobial stewardship and infection prevention and control practices.

Methods

Eighty-five CRKP isolates from hospitalized patients underwent WGS as part of a prospective, multicentre, cross-sectional study, conducted between 1 November 2020 and 30 November 2022, across public-sector and private-sector hospitals in Cape Town, South Africa.

Results

MLST revealed three novel types, ST6785, ST6786 and ST6787, while the most common were ST219, ST307, ST17, ST13 and ST2497. Different predominant clones were noted in each hospital. The most common carbapenemase gene was blaOXA-48-like, detected in 71% of isolates, with blaNDM detected in 5%. Notably, co-detection of two carbapenemase genes (blaOXA-48-like and blaNDM) occurred in 13% of isolates. The yersiniabactin siderophore was detected in 73% of isolates, and was most commonly associated with the ICEKp5 mobile element. All carbapenemases were located on plasmids. The genes blaOXA-181 and blaOXA-232 colocalized with a ColKP3 replicon type on assembled contigs in 83% and 100% of cases, respectively.

Conclusions

CRKP epidemiology in Cape Town reflects institutionally dominant, rather than regional, clones. The most prevalent carbapenemase gene was blaOXA-48-like, in keeping with CRKP epidemiology in South Africa in general. Emerging clones harbouring both blaOXA-48-like and blaNDM, such as ST17, ST2497 and the novel ST6787, are a concern due to the limited availability of appropriate antimicrobial agents in South Africa.

Inter-species gene flow drives ongoing evolution of Streptococcus pyogenes and Streptococcus dysgalactiae subsp. equisimilis

Streptococcus dysgalactiae subsp. equisimilis (SDSE) is an emerging cause of human infection with invasive disease incidence and clinical manifestations comparable to the closely related species, Streptococcus pyogenes. Through systematic genomic analyses of 501 disseminated SDSE strains, we demonstrate extensive overlap between the genomes of SDSE and S. pyogenes. More than 75% of core genes are shared between the two species with one third demonstrating evidence of cross-species recombination. Twenty-five percent of mobile genetic element (MGE) clusters and 16 of 55 SDSE MGE insertion regions were shared across species. Assessing potential cross-protection from leading S. pyogenes vaccine candidates on SDSE, 12/34 preclinical vaccine antigen genes were shown to be present in >99% of isolates of both species. Relevant to possible vaccine evasion, six vaccine candidate genes demonstrated evidence of inter-species recombination. These findings demonstrate previously unappreciated levels of genomic overlap between these closely related pathogens with implications for streptococcal pathobiology, disease surveillance and prevention.

The first tigecycline resistant Enterococcus faecium in Norway was related to tigecycline exposure

OBJECTIVES

We describe the first tigecycline resistant enterococcal isolate in Norway and the mechanisms involved.

MATERIAL AND METHODS

The Norwegian National Advisory Unit on Detection of Antimicrobial Resistance (K-res). received in 2022 an Enterococcus faecium blood culture isolate with decreased susceptibility to tigecycline from a hospitalized patient in the South-Eastern Norway Health region for confirmatory testing. K-res verified a tigecycline-resistant E. faecium (TigR) with broth microdilution MIC of 0.5 mg/L. The patient had received treatment with tigecycline because of an infection with a linezolid- and vancomycin-resistant but tigecycline susceptible E. faecium (TigS) 47 days prior to the detection of the corresponding tigecycline-resistant isolate. Whole-genome comparisons, cgMLST and SNP analyses revealed that the two ST117 strains were closely related.

RESULTS

The TigR isolate showed a novel deletion of 2 amino acids (K57Y58) in a polymorphic region of ribosomal protein S10 previously associated with tigecycline resistance and a deletion of the tet(M) leader peptide previously related to increased expression of tet(M) and tigecycline resistance in enterococci.

CONCLUSIONS

Genomic and epidemiological analyses confirm that the two E. faecium (TigR and TigS) are closely related isolates of the same strain and that the two deletions (in rpsJ and of tet(M) leader peptide) account for the tigecycline resistance in TigR.

Optimising machine learning prediction of minimum inhibitory concentrations in Klebsiella pneumoniae

Minimum Inhibitory Concentrations (MICs) are the gold standard for quantitatively measuring antibiotic resistance. However, lab-based MIC determination can be time-consuming and suffers from low reproducibility, and interpretation as sensitive or resistant relies on guidelines which change over time. Genome sequencing and machine learning promise to allow in silico MIC prediction as an alternative approach which overcomes some of these difficulties, albeit the interpretation of MIC is still needed. Nevertheless, precisely how we should handle MIC data when dealing with predictive models remains unclear, since they are measured semi-quantitatively, with varying resolution, and are typically also left- and right-censored within varying ranges. We therefore investigated genome-based prediction of MICs in the pathogen Klebsiella pneumoniae using 4367 genomes with both simulated semi-quantitative traits and real MICs. As we were focused on clinical interpretation, we used interpretable rather than black-box machine learning models, namely, Elastic Net, Random Forests, and linear mixed models. Simulated traits were generated accounting for oligogenic, polygenic, and homoplastic genetic effects with different levels of heritability. Then we assessed how model prediction accuracy was affected when MICs were framed as regression and classification. Our results showed that treating the MICs differently depending on the number of concentration levels of antibiotic available was the most promising learning strategy. Specifically, to optimise both prediction accuracy and inference of the correct causal variants, we recommend considering the MICs as continuous and framing the learning problem as a regression when the number of observed antibiotic concentration levels is large, whereas with a smaller number of concentration levels they should be treated as a categorical variable and the learning problem should be framed as a classification. Our findings also underline how predictive models can be improved when prior biological knowledge is taken into account, due to the varying genetic architecture of each antibiotic resistance trait. Finally, we emphasise that incrementing the population database is pivotal for the future clinical implementation of these models to support routine machine-learning based diagnostics.

A convergent evolutionary pathway attenuating cellulose production drives enhanced virulence of some bacteria

Bacteria adapt to selective pressure in their immediate environment in multiple ways. One mechanism involves the acquisition of independent mutations that disable or modify a key pathway, providing a signature of adaptation via convergent evolution. Extra-intestinal pathogenic Escherichia coli (ExPEC) belonging to sequence type 95 (ST95) represent a global clone frequently associated with severe human infections including acute pyelonephritis, sepsis, and neonatal meningitis. Here, we analysed a publicly available dataset of 613 ST95 genomes and identified a series of loss-of-function mutations that disrupt cellulose production or its modification in 55.3% of strains. We show the inability to produce cellulose significantly enhances ST95 invasive infection in a rat model of neonatal meningitis, leading to the disruption of intestinal barrier integrity in newborn pups and enhanced dissemination to the liver, spleen and brain. Consistent with these observations, disruption of cellulose production in ST95 augmented innate immune signalling and tissue neutrophil infiltration in a mouse model of urinary tract infection. Mutations that disrupt cellulose production were also identified in other virulent ExPEC STs, Shigella and Salmonella, suggesting a correlative association with many Enterobacteriaceae that cause severe human infection. Together, our findings provide an explanation for the emergence of hypervirulent Enterobacteriaceae clones.

Genome-based assessment of antimicrobial resistance reveals the lineage specificity of resistance and resistance gene profiles in Riemerella anatipestifer from China

Riemerella anatipestifer (R. anatipestifer) is an important pathogen that causes severe systemic infections in domestic ducks, resulting in substantial economic losses for China's waterfowl industry. Controlling R. anatipestifer with antibiotics is extremely challenging due to its multidrug resistance. Notably, large-scale studies on antimicrobial resistance (AMR) and the corresponding genetic determinants in R. anatipestifer remain scarce. To solve this dilemma, more than 400 nonredundant R. anatipestifer isolates collected from 22 provinces in China between 1994 and 2021 were subjected to broth dilution antibiotic susceptibility assays, and their resistance-associated genetic determinants were characterized by whole-genome sequencing. While over 90% of the isolates was resistant to sulfamethoxazole, kanamycin, gentamicin, ofloxacin, norfloxacin, and trimethoprim, 88.48% of the isolates was resistant to the last-resort drug (tigecycline). Notably, R. anatipestifer resistance to oxacillin, norfloxacin, ofloxacin, and tetracycline was found to increase relatively over time. Genome-wide analysis revealed the alarmingly high prevalence of blaOXA-like (93.05%) and tet(X) (90.64%) genes and the uneven distribution of resistance genes among lineages. Overall, this study reveals a serious AMR situation regarding R. anatipestifer in China, with a high prevalence and high diversity of antimicrobial resistance genes, providing important data for the rational use of antibiotics in veterinary practice.IMPORTANCERiemerella anatipestifer (R. anatipestifer), an important waterfowl pathogen, has caused substantial economic losses worldwide, especially in China. Antimicrobial resistance (AMR) is a major challenge in controlling this pathogen. Although a few studies have reported antimicrobial resistance in R. anatipestifer, comprehensive data remain a gap. This study aims to address the lack of information on R. anatipestifer AMR and its genetic basis. By analyzing more than 400 isolates collected over two decades, this study reveals alarming levels of resistance to several antibiotics, including drugs of last resort. The study also revealed the lineage-specificity of resistance profiles and resistance gene profiles. Overall, this study provides new insights and updated data support for understanding AMR and its genetic determinants in R. anatipestifer.

Pneumococcal population genomics changes during the early time period of conjugate vaccine uptake in southern India

Streptococcus pneumoniae is a major cause of invasive disease of young children in low- and middle-income countries. In southern India, pneumococcal conjugate vaccines (PCVs) that can prevent invasive pneumococcal disease began to be used more frequently after 2015. To characterize pneumococcal evolution during the early time period of PCV uptake in southern India, genomes were sequenced and selected characteristics were determined for 402 invasive isolates collected from children <5 years of age during routine surveillance from 1991 to 2020. Overall, the prevalence and diversity of vaccine type (VT) and non-vaccine type (NVT) isolates did not significantly change post-uptake of PCV. Individually, serotype 1 and global pneumococcal sequence cluster (GPSC or strain lineage) 2 significantly decreased, whereas serotypes 6B, 9V and 19A and GPSCs 1, 6, 10 and 23 significantly increased in proportion post-uptake of PCV. Resistance determinants to penicillin, erythromycin, co-trimoxazole, fluoroquinolones and tetracycline, and multidrug resistance significantly increased in proportion post-uptake of PCV and especially among VT isolates. Co-trimoxazole resistance determinants were common pre- and post-uptake of PCV (85 and 93 %, respectively) and experienced the highest rates of recombination in the genome. Accessory gene frequencies were seen to be changing by small amounts across the frequency spectrum specifically among VT isolates, with the largest changes linked to antimicrobial resistance determinants. In summary, these results indicate that as of 2020 this pneumococcal population was not yet approaching a PCV-induced equilibrium and they highlight changes related to antimicrobial resistance. Augmenting PCV coverage and prudent use of antimicrobials are needed to counter invasive pneumococcal disease in this region.

Identification of a novel CG307 sub-clade in third-generation-cephalosporin-resistant Klebsiella pneumoniae causing invasive infections in the USA

Despite the notable clinical impact, recent molecular epidemiology regarding third-generation-cephalosporin-resistant (3GC-R) Klebsiella pneumoniae in the USA remains limited. We performed whole-genome sequencing of 3GC-R K. pneumoniae bacteraemia isolates collected from March 2016 to May 2022 at a tertiary care cancer centre in Houston, TX, USA, using Illumina and Oxford Nanopore Technologies platforms. A comprehensive comparative genomic analysis was performed to dissect population structure, transmission dynamics and pan-genomic signatures of our 3GC-R K. pneumoniae population. Of the 178 3GC-R K. pneumoniae bacteraemias that occurred during our study time frame, we were able to analyse 153 (86 %) bacteraemia isolates, 126 initial and 27 recurrent isolates. While isolates belonging to the widely prevalent clonal group (CG) 258 were rarely observed, the predominant CG, 307, accounted for 37 (29 %) index isolates and displayed a significant correlation (Pearson correlation test P value=0.03) with the annual frequency of 3GC-R K. pneumoniae bacteraemia. Interestingly, only 11 % (4/37) of CG307 isolates belonged to the commonly detected 'Texas-specific' clade that has been observed in previous Texas-based K. pneumoniae antimicrobial-resistance surveillance studies. We identified nearly half of our CG307 isolates (n=18) belonged to a novel, monophyletic CG307 sub-clade characterized by the chromosomally encoded bla SHV-205 and unique accessory genome content. This CG307 sub-clade was detected in various regions of the USA, with genome sequences from 24 additional strains becoming recently available in the National Center for Biotechnology Information (NCBI) SRA database. Collectively, this study underscores the emergence and dissemination of a distinct CG307 sub-clade that is a prevalent cause of 3GC-R K. pneumoniae bacteraemia among cancer patients seen in Houston, TX, and has recently been isolated throughout the USA.

Molecular Genomic Analyses of Enterococcus cecorum from Sepsis Outbreaks in Broilers

Extensive genomic analyses of Enterococcus cecorum isolates from sepsis outbreaks in broilers suggest a polyphyletic origin, likely arising from core genome mutations rather than gene acquisition. This species is a normal intestinal flora of avian species with particular isolates associated with osteomyelitis. More recently, this species has been associated with sepsis outbreaks affecting broilers during the first 3 weeks post-hatch. Understanding the genetic and management basis of this new phenotype is critical for developing strategies to mitigate this emerging problem. Phylogenomic analyses of 227 genomes suggest that sepsis isolates are polyphyletic and closely related to both commensal and osteomyelitis isolate genomes. Pangenome analyses detect no gene acquisitions that distinguish all the sepsis isolates. Core genome single nucleotide polymorphism analyses have identified a number of mutations, affecting the protein-coding sequences, that are enriched in sepsis isolates. The analysis of the protein substitutions supports the mutational origins of sepsis isolates.

Streptococcus pneumoniae serotype 33G: genetic, serological, and structural analysis of a new capsule type

IMPORTANCE

Streptococcus pneumoniae (the pneumococcus) is a bacterial pathogen with the greatest burden of disease in Asia and Africa. The pneumococcal capsular polysaccharide has biological relevance as a major virulence factor as well as public health importance as it is the target for currently licensed vaccines. These vaccines have limited valency, covering up to 23 of the >100 known capsular types (serotypes) with higher valency vaccines in development. Here, we have characterized a new pneumococcal serotype, which we have named 33G. We detected serotype 33G in nasopharyngeal swabs (n = 20) from children and adults hospitalized with pneumonia, as well as healthy children in Mongolia. We show that the genetic, serological, and biochemical properties of 33G differ from existing serotypes, satisfying the criteria to be designated as a new serotype. Future studies should focus on the geographical distribution of 33G and any changes in prevalence following vaccine introduction.

WGS of a cluster of MDR Shigella sonnei utilizing Oxford Nanopore R10.4.1 long-read sequencing

OBJECTIVES

To utilize long-read nanopore sequencing (R10.4.1 flowcells) for WGS of a cluster of MDR Shigella sonnei, specifically characterizing genetic predictors of antimicrobial resistance (AMR).

METHODS

WGS was performed on S. sonnei isolates identified from stool and blood between September 2021 and October 2022. Bacterial DNA from clinical isolates was extracted on the MagNA Pure 24 and sequenced on the GridION utilizing R10.4.1 flowcells. Phenotypic antimicrobial susceptibility testing was interpreted based on CLSI breakpoints. Sequencing data were processed with BugSeq, and AMR was assessed with BugSplit and ResFinder.

RESULTS

Fifty-six isolates were sequenced, including 53 related to the cluster of cases. All cluster isolates were identified as S. sonnei by sequencing, with global genotype 3.6.1.1.2 (CipR.MSM5), MLST 152 and PopPUNK cluster 3. Core genome MLST (cgMLST, examining 2513 loci) and reference-based MLST (refMLST, examining 4091 loci) both confirmed the clonality of the isolates. Cluster isolates were resistant to ampicillin (blaTEM-1), trimethoprim/sulfamethoxazole (dfA1, dfrA17; sul1, sul2), azithromycin (ermB, mphA) and ciprofloxacin (gyrA S83L, gyrA D87G, parC S80I). No genomic predictors of resistance to carbapenems were identified.

CONCLUSIONS

WGS with R10.4.1 enabled rapid sequencing and identification of an MDR S. sonnei community cluster. Genetic predictors of AMR were concordant with phenotypic antimicrobial susceptibility testing.

Strain features of pneumococcal isolates in the pre- and post-PCV10 era in Pakistan

Pakistan is amongst the four countries with the highest number of pneumococcal deaths. While the PCV10 vaccine was introduced in Pakistan in October 2012, data regarding the impact of the vaccine on the population dynamics of Streptococcus pneumoniae in Pakistan remain obscure. Using whole genome sequencing of 190 isolates (nasopharyngeal carriage=75, disease=113, unknown sites=2) collected between 2002 and 2020, this study presents characteristics of pneumococcal strains in Pakistan in the pre- and post-vaccine era. The isolates were characterized on the basis of serotype distribution, genetic lineages (or Global Pneumococcal Sequence Cluster, GPSC) and antibiotic resistance. A high level of diversity in serotype and genetic lineages of pneumococci was observed in Pakistan. Among 190 isolates, we identified 54 serotypes, 67 GPSCs and 116 sequence types (STs) including 23 new STs. The most prevalent GPSCs and their associated serotypes in nasopharyngeal carriage were GPSC54 (expressing serotype 9V), GPSC5 (15A and 7B, and serogroup 24), GPSC25 (15B/15C), GPSC67 (18C) and GPSC376 (6A and 6D). Similarly, among 113 disease-causing isolates, the most prevalent GPSC/serotype combinations were GPSC2 (serotype 1), GPSC10 (serotypes 14, 10A, 19A and 19F), GPSC43 (serotypes 13, 11A, 23B, 35A and 9V), GPSC67 (serotypes 18A and 18C) and GPSC642 (serotype 11A). Of the 190 isolates, the highest levels of resistance were observed against penicillin (58.9 %, n=122), erythromycin (29.5 %, n=56), clindamycin (13.2 %, n=25), co-trimoxazole (94.2 %, n=179) and tetracycline/doxycycline (53.2 %, n=101). A higher proportion of disease-causing isolates were multidrug resistant as compared to carriage isolates (54 % vs 25 %). Our data suggest limited coverage of PCV10 in nasopharyngeal (21.6 %, 16/74) as well as disease-causing (38.1 %, 16/42) isolates among children ≤5 years old; however, higher valent vaccine PCV13 would increase the coverage rates to 33.8 % in nasopharyngeal and 54.8 % in disease-causing isolates, whereas PCV24/25 would offer the highest coverage rates. Owing to the diversity of serotypes observed during the post-vaccine period, the suggested inclusion of serotype in future vaccine formulations will require investigations with larger data sets with an extended temporal window. This article contains data hosted by Microreact.

Emergence and mobilization of a novel lincosamide resistance gene lnu(I): From environmental reservoirs to pathogenic bacteria

Antimicrobial resistance remains an utmost concern in human and veterinary medicine, impacting humans, animals, and the environment while significantly influencing the principles of One Health. While Riemerella anatipestifer (R. anatipestifer) is recognized as a waterfowl pathogen with multidrug-resistant properties, the specifics of its lincosamide resistance mechanism are inadequately understood. In this study, we identified a novel lincosamide resistance gene, lnu(I), in R. anatipestifer RCAD0121, and investigated its potential origin, transfer mechanisms, and dissemination status through genomic epidemiology. This exhibited 74.80 % amino acid identity with a previously reported gene, lnu(H). PCR analysis revealed lnu(I) prevalence in at least 44 R. anatipestifer isolates collected from multiple provinces in China. Furthermore, genomic mining unveiled 56 lnu(I) sequences within publicly available databases, primarily originating from environmental sources. In addition, members of the family Flavobacteriaceae were the dominant (16/56, 28.57 %) bacteria carrying the lnu(I) gene, with Flavobacterium exhibiting a similar GC content as lnu(I). Notably, specific instances of the lnu(I) gene were linked to mobile genetic elements within human and animal pathogenic bacteria. These findings suggest that Flavobacterium species within the environment could serve as potential ancestral sources of the novel lnu(I) gene, which has undergone mobilization events toward pathogenic bacteria.

The population-level impact of Enterococcus faecalis genetics on intestinal colonization and extraintestinal infection

IMPORTANCE

Enterococcus faecalis causes life-threatening invasive hospital- and community-associated infections that are usually associated with multidrug resistance globally. Although E. faecalis infections cause opportunistic infections typically associated with antibiotic use, immunocompromised immune status, and other factors, they also possess an arsenal of virulence factors crucial for their pathogenicity. Despite this, the relative contribution of these virulence factors and other genetic changes to the pathogenicity of E. faecalis strains remain poorly understood. Here, we investigated whether specific genomic changes in the genome of E. faecalis isolates influence its pathogenicity-infection of hospitalized and nonhospitalized individuals and the propensity to cause extraintestinal infection and intestinal colonization. Our findings indicate that E. faecalis genetics partially influence the infection of hospitalized and nonhospitalized individuals and the propensity to cause extraintestinal infection, possibly due to gut-to-bloodstream translocation, highlighting the potential substantial role of host and environmental factors, including gut microbiota, on the opportunistic pathogenic lifestyle of this bacterium.

Municipal wastewater monitoring revealed the predominance of bla GES genes with diverse variants among carbapenemase-producing organisms: high occurrence and persistence of Aeromonas caviae harboring the new bla GES variant bla GES-48

IMPORTANCE

The emergence and spread of carbapenemase-producing organisms (CPOs) represent a global health threat because they are associated with limited treatment options and poor clinical outcomes. Wastewater is considered a hotspot for the evolution and dissemination of antimicrobial resistance. Thus, analyses of municipal wastewater are critical for understanding the circulation of these CPOs and carbapenemase genes in local communities, which remains scarcely known in Japan. This study resulted in several key observations: (i) the vast majority of bla GES genes, including six new bla GES variants, and less frequent bla IMP genes were carbapenemase genes encountered exclusively in wastewater influent; (ii) the most dominant CPO species were Aeromonas spp., in which a remarkable diversity of new sequence types was observed; and (iii) CPOs were detected from combined sewer wastewater, but not from separate sewer wastewater, suggesting that the load of CPOs from unrecognized environmental sources could greatly contribute to their detection in influent wastewater.

Streptococcus ruminantium-associated sheep mastitis outbreak detected in Italy is distinct from bovine isolates

Streptococcus ruminantium is the causative agent of several bovine and ovine diseases, however reports are uncommon and application of whole genome sequencing to identify is rare. We report for the first time, a severe ovine mastitis outbreak caused by S. ruminantium in Italy, 2022. S. ruminantium was isolated from 12 adult lactating ewes with diffuse nodules in the mammary parenchyma and predominantly serous and clotted milk. All outbreak isolates, along with five additional historical Italian isolates (between 2011 and 2017), were genomically characterised and then analysed in the context of all publicly available S. ruminantium genomes. Antimicrobial susceptibility testing was performed to determine the MICs of 16 antibiotics. The results showed that all isolates were susceptible to all antimicrobials tested except kanamycin. Single Nucleotide Variant analysis confirmed this as a clonal outbreak across 10 sheep (≤ 15 SNVs), while the two others were colonised by more distantly related clones (≤ 53 pairwise SNVs), indicating the presence of multiple infecting lineages. The five historical S. ruminantium isolates were comprised of genetically-distant singletons (between 1259 and 5430 pairwise SNVs to 2022 outbreak isolates). Ovine isolates were found to be genetically distinct to bovine isolates, forming monophyletic groups. Bovine isolates were similarly made up of singleton clones in all but two isolates. Taken together, our genomic analysis using all globally available genomes is consistent with general opportunistic pathogenesis of S. ruminantium. We encourage future genomic surveillance efforts to facilitate outbreak detection, as well as improve our understanding of this poorly-understood, multi-host, zoonotic pathogen.

Prediction of post-PCV13 pneumococcal evolution using invasive disease data enhanced by inverse-invasiveness weighting

Background

After introduction of pneumococcal conjugate vaccines (PCVs), serotype replacement occurred in the population of Streptococcus pneumoniae. Predicting which pneumococcal clones and serotypes will become more common in carriage after vaccination can enhance vaccine design and public health interventions, while also improving our understanding of pneumococcal evolution. We sought to use invasive disease data to assess how well negative frequency-dependent selection (NFDS) models could explain pneumococcal carriage population evolution in the post-PCV13 epoch by weighting invasive data to approximate strain proportions in the carriage population.

Methods

Invasive pneumococcal isolates were collected and sequenced during 1998-2018 by the Active Bacterial Core surveillance (ABCs) from the Centers for Disease Control and Prevention (CDC). To predict the post-PCV13 population dynamics in the carriage population using a NFDS model, all genomic data were processed under a bioinformatic pipeline of assembly, annotation, and pangenome analysis to define genetically similar sequence clusters (i.e., strains) and a set of accessory genes present in 5% to 95% of the isolates. The NFDS model predicted the strain proportion by calculating the post-vaccine strain composition in the weighted invasive disease population that would best match pre-vaccine accessory gene frequencies. To overcome the biases of invasive disease data, serotype-specific inverse-invasiveness weights were defined as the ratio of the proportion of the serotype in the carriage data to the proportion in the invasive data, using data from 1998-2001 in the United States, before conjugate vaccine introduction. The weights were applied to adjust both the observed strain proportion and the accessory gene frequencies.

Results

Inverse-invasiveness weighting increased the correlation of accessory gene frequencies between invasive and carriage data with reduced residuals in linear or logit scale for pre-vaccine, post-PCV7, and post-PCV13. Similarly, weighting increased the correlation of accessory gene frequencies between different time periods in the invasive data. By weighting the invasive data, we were able to use the NFDS model to predict strain proportions in the carriage population in the post-PCV13 epoch, with the adjusted R-squared between predicted and observed strain proportions increasing from 0.176 to 0.544 after weighting.

Conclusions

The weighting system adjusted the invasive disease surveillance data to better represent the carriage population of S. pneumoniae. The NFDS mechanism predicted the strain proportions in the projected carriage population as estimated from the weighted invasive disease frequencies in the post-PCV13 epoch. Our methods enrich the value of genomic sequences from invasive disease surveillance, which is readily available, easy to collect, and of direct interest to public health.