SeroBA: rapid high-throughput serotyping of Streptococcus pneumoniae from whole genome sequence data

Changes in the incidence, viral coinfection pattern and outcomes of pneumococcal hospitalizations during and after the COVID-19 pandemic

BACKGROUND

The incidence of pneumococcal pneumonia in the context of the Coronavirus Disease 2019 (COVID-19) pandemic, along with the real-world data on the ratio of non-invasive to invasive pneumococcal pneumonia, is an area that has not been thoroughly studied. The outcomes associated with coinfection of influenza and COVID-19 remain unknown. This study examined the incidence, demographics, coinfection with influenza and/or COVID-19, and clinical outcomes of pneumococcal hospitalizations in Hong Kong during the baseline, pandemic, and post-pandemic periods.

METHODS

Hospitalization records of individuals aged 18 years and above with pneumococcal disease from January 2015 to August 2024 were extracted from the territory-wide electronic medical record database. Pneumococcal disease was categorized into invasive pneumococcal pneumonia (IPP), invasive pneumococcal disease without pneumonia (IPDWP), and non-invasive pneumococcal pneumonia (NIPP), followed by univariate and multivariate analyses. Effects of coinfection with influenza and COVID-19 were analyzed.

RESULTS

The incidence of pneumococcal disease decreased during the COVID-19 pandemic but rebounded in the post-pandemic period. There were no significant changes in the distribution of pneumococcal serotypes across the three periods. The study revealed a strong positive correlation between monthly pneumococcal hospitalizations and the indicator of influenza activity, while the correlation with the COVID-19 indicator was weak. Additionally, strong positive correlations were observed between the indicator of influenza activity and influenza coinfections, as well as between the indicator of COVID-19 activity and COVID-19 coinfections. Multivariate analyses identified male gender, a higher comorbidity index, older age, invasive pneumococcal disease (IPP and IPDWP), coinfection with influenza and COVID-19, and hospitalization during the pandemic period as factors associated with adverse outcomes.

CONCLUSIONS

The study showcases changes in the epidemiology of pneumococcal disease during and after the COVID-19 pandemic. It highlights the adverse effects of influenza and COVID-19 coinfections on the outcomes of patients with pneumococcal disease and emphasizes the need to vaccinate vulnerable populations against these infections.

An abundance of aliC and aliD genes were identified in saliva using a novel multiplex qPCR to characterize group II non-encapsulated pneumococci with improved specificity

Pneumococcal surveillance studies are reporting increasing prevalence of non-encapsulated pneumococci (NESp). NESp are an important reservoir for genetic exchange among streptococci, including for antimicrobial resistance, and are increasingly implicated in disease. Disease-associated NESp commonly carries the virulence genes pspK, or aliC and aliD in their cps locus instead of capsule genes. While molecular methods targeting the cps region are widely used for serotyping encapsulated strains, there are few assays available for the classification of NESp, meaning it is not widely undertaken. Therefore, we exploited these genes as targets for a novel qPCR assay for detecting and classifying NESp strains with improved efficiency and specificity. We conducted bioinformatic analysis on sequences from 402 NESp and 45 other mitis-group streptococci and developed a multiplex-qPCR, targeting pspK, aliD and two regions of aliC. The assay was validated using 16 previously identified NESp isolates. We then applied the assay to DNA extracted from culture-enriched saliva and isolated and characterized suspected NESp colonies, with confirmation by whole genome sequencing. Bioinformatic analyses demonstrated that previously published primers for aliC and aliD had low pneumococcal specificity but indicated that targeting two regions of aliC would improve species specificity, without compromising sensitivity. Our novel multiplex assay accurately typed all isolates. When screening saliva, we found a high prevalence of aliC and aliD, even in samples negative for pneumococcal genes lytA and piaB. Isolated colonies which were aliC and aliD positive could be differentiated as non-pneumococcal streptococci using our assay. Our multiplex-qPCR assay can be used to efficiently screen even highly polymicrobial samples, such as saliva, for NESp genes, to detect and differentiate potentially pathogenic NESp clades from closely related mitis-group streptococci. This will allow for a better understanding of the true prevalence of NESp and its impact on pneumococcal carriage and disease.

Comparative pangenomics of Streptococcus pneumoniae from Malawi: uncovering genetic variability and pathogenicity

Streptococcus pneumoniae is a significant cause of bacterial infections, including pneumonia, meningitis and septicemia, primarily affecting children, the elderly and immunocompromised individuals. This study aimed to elucidate the serotype and lineage distribution and molecular mechanisms underlying pneumococcal invasiveness through a comprehensive pangenomic analysis of 1416 isolates from Malawi. Our analysis comprised 810 isolates from asymptomatic carriers and 606 isolates from patients with bacteraemia or meningitis. We identified 58 serotypes, with serotypes 1, 5 and 12F exhibiting significantly higher prevalence among patients. These serotypes likely exhibit reduced nasopharyngeal colonization and demonstrate rapid dissemination to sterile sites. Notably, these serotypes form a distinct lineage with distinct genomic characteristics, including the absence of V-type ATP synthase. The pangenome analysis revealed two highly conserved surface protein complexes, F-type ATP synthase and SecA1-SecY, which deserve further investigation as potential targets for novel therapeutic interventions.

A global epidemic serotype 14 Streptococcus pneumoniae switching to non-vaccine types

Streptococcus pneumoniae infection is a leading cause of bacterial pneumonia worldwide, contributing to both invasive pneumococcal diseases and non-invasive diseases. Following the introduction of pneumococcal conjugate vaccines (PCV7, PCV10, and PCV13), a significant decrease in the prevalence of vaccine-type pneumococci has been observed, while the prevalence of non-vaccine-type (NVT) pneumococci has increased. This shift has been evident in both disease incidence and nasopharyngeal carriage. Serotype switching refers to the alteration or exchange of the serotype-related cps gene locus that has been flagged globally as serotypes of concern. In this study, we collected 166 S. pneumoniae clinical isolates and analyzed the serotypes of isolates to evaluate vaccine effectiveness. Additionally, we identified multidrug-resistant NVT S. pneumoniae isolates that recently emerged within the epidemic CC876 complex, a pandemic serotype 14 clone. More importantly, we identified a serotype 15B/ST4749 strain, which might have been generated through a serotype-switch event from serotype 14/ST4749. Mapping the genomic recombination sites within the serotype 15B/ST4749 progeny revealed a donated fragment encompassing cps, pbp1a, and additional key factors, including the DNA uptake-related gene comC. More importantly, as a strain that is already prevalent, the novel NVT strain possesses the potential to become a global pandemic strain.IMPORTANCEThe study employed high-throughput sequencing to analyze Streptococcus pneumoniae isolates from Xiamen Children's Hospital, China, to evaluate vaccine effectiveness in light of serotype changes and documented the occurrence of a serotype switch (from 14 to 15B) between vaccine-type and non-vaccine-type serotypes in the ST4749 strain. This observation indicates a genomic recombination and adaptive response of S. pneumoniae under selective pressure exerted by vaccination, offering novel insights into strategies for the prevention and control of S. pneumoniae and the optimization of vaccine deployment.

Update on the evolving landscape of pneumococcal capsule types: new discoveries and way forward

SUMMARYStreptococcus pneumoniae (the "pneumococcus") is a significant human pathogen. The key determinant of pneumococcal fitness and virulence is its ability to produce a protective polysaccharide (PS) capsule, and anti-capsule antibodies mediate serotype-specific opsonophagocytic killing of bacteria. Notably, immunization with pneumococcal conjugate vaccines (PCVs) has effectively reduced the burden of disease caused by serotypes included in vaccines but has also spurred a relative upsurge in the prevalence of non-vaccine serotypes. Recent advancements in serotyping and bioinformatics surveillance tools coupled with high-resolution analytical techniques have enabled the discovery of numerous new capsule types, thereby providing a fresh perspective on the dynamic pneumococcal landscape. This review offers insights into the current pneumococcal seroepidemiology highlighting important serotype shifts in different global regions in the PCV era. It also comprehensively summarizes newly discovered serotypes from 2007 to 2024, alongside updates on revised chemical structures and the de-novo determinations of structures for previously known serotypes. Furthermore, we spotlight emerging evidence on non-pneumococcal Mitis-group strains that express capsular PS that are serologically and biochemically related to the pneumococcal capsule types. We further discuss the implications of these recent findings on capsule nomenclature, pneumococcal carriage detection, and future PCV design. The review maps out the current status and also outlines the course for future research and vaccine strategies, ensuring a continued effective response to the evolving pneumococcal challenge.

Genomic characterization of Streptococcus pneumoniae isolates obtained from carriage and disease among paediatric patients in Addis Ababa, Ethiopia

Background and aims. Despite the introduction of pneumococcal conjugate vaccines (PCVs), Streptococcus pneumoniae still remains an important cause of morbidity and mortality, especially among children under 5 years in sub-Saharan Africa. We sought to determine the distribution of serotypes, lineages and antimicrobial resistance of S. pneumoniae from carriage and disease among children presenting to health facilities, 5-6 years after the introduction of PCV10 in Ethiopia.Methods. Whole-genome sequencing (WGS) was performed on 103 S. pneumoniae (86 from nasopharyngeal swabs, 4 from blood and 13 from middle ear discharge) isolated from children aged <15 years at 3 healthcare facilities in Addis Ababa, Ethiopia, from September 2016 to August 2017. Using the WGS data, serotypes were predicted, isolates were assigned to clonal complexes, global pneumococcal sequence clusters (GPSCs) were inferred and screening for alleles and mutations that confer resistance to antibiotics was performed using multiple bioinformatic pipelines.Results. The 103 S. pneumoniae isolates were assigned to 38 serotypes (including nontypeable) and 46 different GPSCs. The most common serotype was serotype 19A. Common GPSCs were GPSC1 [14.6% (15/103), sequence type (ST) 320, serotype 19A], GPSC268 [8.7% (9/103), ST 6882 and novel STs; serotypes 16F, 11A and 35A] and GPSC10 [8.7% (9/103), STs 2013, 230 and 8804; serotype 19A]. The four invasive isolates were serotype 19A (n=2) and serotype 33C (n=2). Resistance to penicillin (>0.06 µg ml-1, CLSI meningitis cutoff) was predicted in 57% (59/103) of the isolates, and 43% (25/58) penicillin-binding protein allele combinations were predicted to be associated with penicillin resistance. Resistance mutations in folA (I100L) and/or folP (indel between fifty-sixth and sixty-seventh aa) were identified among 66% (68/103) of the isolates, whilst tetracycline (tetM) and macrolide (ermB and mefA) resistance genes were found in 46.6% (48/103), 20.4% (21/103) and 20.4% (21/103) of the isolates, respectively. Multidrug resistance (MDR) (≥3 antibiotic classes) was observed in 31.1% (32/103) of the isolates. GPSC1 and GPSC10 accounted for 46.8% (15/32) and 18.7% (6/32) of the overall MDR.Conclusion. Five to 6 years after the introduction of PCV10 in Ethiopia, the S. pneumoniae obtained from carriage and disease among paediatric patients showed diverse serotype and pneumococcal lineages. The most common serotype identified was 19A, expressed by the MDR lineages GPSC1 and GPSC10, which is not covered by PCV10 but is included in PCV13. Continued assessment of the impact of PCV on the population structure of S. pneumoniae in Ethiopia is warranted during and after PCV13 introduction.

Investigating two decades of Streptococcus pneumoniae bacteraemia in the Gelderland area, the Netherlands, using whole-genome sequencing

In the Netherlands, the 7-valent pneumococcal conjugate vaccine (PCV) was introduced to the childhood immunization programme in 2006 and replaced by the 10-valent PCV (PCV10, GSK) in 2011. To describe invasive pneumococcal disease in the era of childhood PCV vaccination on pneumococcal bacteraemia across all ages, we collected and sequenced 979 pneumococcal blood isolates from consecutive patients with pneumococcal bacteraemia in the Gelderland area, the Netherlands, between 2000 and 2020. In total, 58% of the bacteraemia cases (n=563/979) occurred in the elderly population. Compared to the pre-PCV period (2000-2005), the odds ratio for non-PCV10 bacteraemia was 17.5 (CI 9.9-31.6; P<0.001) in the late-PCV10 period, showing an overall increase in the proportion of bacteraemia cases being caused by non-vaccine serotype pneumococci (2016-2020). The increase in non-PCV10 serotypes is mainly driven by an expansion of lineage global pneumococcal sequencing cluster 3 (GPSC3) expressing serotype 8, alongside the emergence of serotype 12F that was mediated by multiple lineages (GPSC32/GPSC26/GPSC55). Both serotypes 8 and 12F were included in the latest PCV20 formulation that is licensed to be used in children and adults in Europe. Over 20 years, we observed a low prevalence of antimicrobial resistance (AMR) as predicted by genome data. There were no significant changes in AMR prevalence after vaccine introduction (P>0.05 for all comparisons). We saw a stably low prevalence of reduced penicillin susceptibility, which was observed in multiple pneumococcal lineages, with GPSC10 being the most common in the Gelderland collection, whilst GPSC1 and GPSC6 were common among the penicillin-resistant pneumococcal blood culture isolates provided by the Netherlands Reference Laboratory for Bacterial Meningitis. Comparison to global collections of GPSC10, GPSC1 and GPSC6 isolates favored the likelihood of separate introductions of penicillin-resistant isolates rather than cloncal expansion. Genomic surveillance of pneumococcal bacteraemia in this unbiased population sample in the Netherlands supports the use of higher valency PCVs, such as PCV20, especially in adults, to prevent future bacteraemia cases caused by Streptococcus pneumoniae in the Gelderland area, the Netherlands, while maintaining a low prevalence of AMR in the pneumococcal population.

PneusPage: A WEB-BASED TOOL for the analysis of Whole-Genome Sequencing Data of Streptococcus pneumonia

With the advent of whole-genome sequencing, opportunities to investigate the population structure, transmission patterns, antimicrobial resistance profiles, and virulence determinants of Streptococcus pneumoniae at high resolution have been increasingly expanding. Consequently, a user-friendly bioinformatics tool is needed to automate the analysis of Streptococcus pneumoniae whole-genome sequencing data, summarize clinically relevant genomic features, and further guide treatment options. Here, we developed PneusPage, a web-based tool that integrates functions for species prediction, molecular typing, drug resistance determination, and data visualization of Streptococcus pneumoniae. To evaluate the performance of PneusPage, we analyzed 80 pneumococcal genomes with different serotypes from the Global Pneumococcal Sequencing Project and compared the results with those from another platform, PathogenWatch. We observed a high concordance between the two platforms in terms of serotypes (100% concordance rate), multilocus sequence typing (100% concordance rate), penicillin-binding protein typing (88.8% concordance rate), and the Global Pneumococcal Sequencing Clusters (98.8% concordance rate). In addition, PneusPage offers integrated analysis functions for the detection of virulence and mobile genetic elements that are not provided by previous platforms. By automating the analysis pipeline, PneusPage makes whole-genome sequencing data more accessible to non-specialist users, including microbiologists, epidemiologists, and clinicians, thereby enhancing the utility of whole-genome sequencing in both research and clinical settings. PneusPage is available at https://pneuspage.minholee.net/.

A genome-based survey of invasive pneumococci in Norway over four decades reveals lineage-specific responses to vaccination

BACKGROUND

Streptococcus pneumoniae is a major cause of mortality globally. The introduction of pneumococcal conjugate vaccines (PCVs) has reduced the incidence of the targeted serotypes significantly, but expansion of non-targeted serotypes, serotype replacement, and incomplete vaccine-targeting contribute to pneumococcal disease in the vaccine era. Here, we characterize the changing population genetic landscape of S. pneumoniae in Norway over a 41-year period (1982-2022).

METHODS

Since 2018, all cases of invasive pneumococcal disease have undergone whole-genome sequencing (WGS) at the Norwegian Institute of Public Health. In order to characterize the changing population over time, historical isolates were re-cultured and sequenced, resulting in a historical WGS dataset. Isolates were assigned to global pneumococcal sequence clusters (GPSCs) using PathogenWatch and assigned to serotypes using in silico (SeroBA) and in vitro methods (Quellung reaction). Temporal phylogenetic analyses were performed on GPSCs of particular interest.

RESULTS

The availability of WGS data allowed us to study capsular variation at the level of individual lineages. We detect highly divergent fates for different GPSCs following the introduction of PCVs. For two out of eight major GPSCs, we identified multiple instances of serotype switching from vaccine types to non-vaccine types. Dating analyses suggest that most instances of serotype switching predated the introduction of PCVs, but expansion occurred after their introduction. Furthermore, selection for penicillin non-susceptibility was not a driving force for the changing serotype distribution within the GPSCs over time.

CONCLUSIONS

PCVs have been major shapers of the Norwegian disease-causing pneumococcal population, both at the level of serotype distributions and the underlying lineage dynamics. Overall, the introduction of PCVs has reduced the incidence of invasive disease. However, some GPSCs initially dominated by vaccine types escaped the effect of vaccination through expansion of non-vaccine serotypes. Close monitoring of circulating lineages and serotypes will be key for ensuring optimal vaccination coverage going forward.

Molecular epidemiology of Streptococcus pneumoniae isolates causing invasive and noninvasive infection in Ethiopia

Streptococcus pneumoniae, a medically important opportunistic bacterial pathogen of the upper respiratory tract, is a major public health concern, causing a wide range of pneumococcal illnesses, both invasive and noninvasive. It is associated with significant global morbidity and mortality, including pneumonia, meningitis, sepsis, and acute otitis media. The major purpose of this study was to determine the molecular epidemiology of Streptococcus pneumoniae strains that cause invasive and noninvasive infections in Ethiopia. A prospective study was undertaken in two regional hospitals between January 2018 and December 2019. Whole-genome sequencing was used to analyze all isolates. Serotypes and multilocus sequence types (MLST) were derived from genomic data. The E-test was used for antimicrobial susceptibility testing. Patient samples obtained 54 Streptococcus pneumoniae isolates, 33 from invasive and 21 from noninvasive specimens. Our findings identified 32 serotypes expressed by 25 Global Pneumococcal Sequence Clusters (GPSCs) and 42 sequence types (STs), including 21 new STs. The most common sequence types among the invasive isolates were ST3500, ST5368, ST11162, ST15425, ST15555, ST15559, and ST15561 (2/33, 6% each). These sequence types were linked to serotypes 8, 7 C, 15B/C, 16 F, 10 A, 15B, and 6 A, respectively. Among the noninvasive isolates, only ST15432, associated with serotype 23 A, had numerous isolates (4/21, 19%). Serotype 14 was revealed as the most resistant strain to penicillin G, whereas isolates from serotypes 3, 8, 7 C, and 10 A were resistant to erythromycin. Notably, all serotype 6 A isolates were resistant to both erythromycin and penicillin G. Our findings revealed an abnormally significant number of novel STs, as well as extremely diversified serotypes and sequence types, implying that Ethiopia may serve as a breeding ground for novel STs. Recombination can produce novel STs that cause capsular switching. This has the potential to influence how immunization campaigns affect the burden of invasive pneumococcal illness. The findings highlight the importance of continuous genetic surveillance of the pneumococcal population as a vital step toward enhancing future vaccine design.

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.

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.

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.

A narrative review of genomic characteristics, serotype, immunogenicity, and vaccine development of Streptococcus pneumoniae capsular polysaccharide

This narrative review describes genomic characteristic, serotyping, immunogenicity, and vaccine development of Streptococcus pneumoniae capsular polysaccharide (CPS). CPS is a primary virulence factor of S. pneumoniae. The genomic characteristics of S. pneumoniae CPS, including the role of biosynthetic gene and genetic variation within cps (capsule polysaccharide) locus which may lead to serotype replacement are still being investigated. One hundred unique serotypes of S. pneumoniae have been identified through various methods of serotyping using phenotypic and genotypic approach. The advantages and limitations of each method are various, emphasizing the need for accurate and comprehensive serotyping for effective disease surveillance and vaccine targeting. In addition, we elaborate the critical role of CPS in vaccine development by providing an overview of immunogenicity, ongoing research of pneumococcal vaccines, and the impact on disease burden.

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.

Prevalence and resistance characteristics of multidrug-resistant Streptococcus pneumoniae isolated from the respiratory tracts of hospitalized children in Shenzhen, China

Background

PCV13 introduction in China has led to a significant reduction of vaccine serotype Streptococcus pneumoniae. However, non-vaccine serotypes with highly resistance and invasiveness were often reported in the post-pneumococcal conjugate vaccine era and there was regional differences.

Methods

A total of 669 S. pneumoniae strains were collected from the respiratory tracts of hospitalized children at Shenzhen Children's Hospital in 2021 and 2022. Antimicrobial resistance (AMR) characteristics were assessed through antibiotic susceptibility testing performed with the VITEK 2 compact system. AMR genes and single nucleotide polymorphisms (SNPs) in pbp1a, pbp2b, and pbp2x were identified via analysis of whole genome sequencing data. Statistical examination of the data was conducted employing chi-square and Fisher's exact tests.

Results

We found that non-vaccine serotypes strains had accounted for 46.6% of all the pneumococcal isolated strains. The most common non-vaccine serotype is 23A, with a prevalence rate of 8.9%, followed by 15A (6.6%), 6E (5.7%), 34 (3.2%), and 15B (2.9%). The multidrug resistance rates (MDR) of vaccine serotypes were 19F (99.36%), 19A (100%), 23F (98.08%), 6B (100%), and 6C (100%). Meanwhile, the MDR of non-vaccine serotypes were 15B (100.00%), 6E (100%), 15C (100%), 34 (95.24%), and 23A (98.31%). Resistance rates of 6E to more than six antibiotic classes reached 89.47%, which is similar to 19F (83.33%) and 19A (90%). Unique resistance profiles were also identified for non-vaccine serotypes, including significantly higher resistance to chloramphenicol in 6E, 15B, and 15C than in 19F and 19A. Furthermore, through genome sequencing, we revealed strong correlation of cat-TC with chloramphenicol resistance, patA/patB with tetracycline resistance, ermB and pmrA with erythromycin resistance.

Conclusion

The introduction of PCV13 into China from 2017 has led to a shift in the dominant composition of pneumococcal strains. There has been a notable rise and spread of multidrug-resistant non-vaccine serotypes among children. Specifically, the non-vaccine serotype 6E, which was not widely reported in China previously, has emerged. To comprehend the resistance mechanisms, it is crucial to further investigate the molecular and genetic characteristics of these non-vaccine serotypes.

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.

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

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

Exploring clonality and virulence gene associations in bloodstream infections using whole-genome sequencing and clinical data

BACKGROUND

Bloodstream infections (BSIs) remain a significant cause of mortality worldwide. Causative pathogens are routinely identified and susceptibility tested but only very rarely investigated for their resistance genes, virulence factors, and clonality. Our aim was to gain insight into the clonality patterns of different species causing BSI and the clinical relevance of distinct virulence genes.

METHODS

For this study, we whole-genome-sequenced over 400 randomly selected important pathogens isolated from blood cultures in our diagnostic department between 2016 and 2021. Genomic data on virulence factors, resistance genes, and clonality were cross-linked with in-vitro data and demographic and clinical information.

RESULTS

The investigation yielded extensive and informative data on the distribution of genes implicated in BSI as well as on the clonality of isolates across various species.

CONCLUSION

Associations between survival outcomes and the presence of specific genes must be interpreted with caution, and conducting replication studies with larger sample sizes for each species appears mandatory. Likewise, a deeper knowledge of virulence and host factors will aid in the interpretation of results and might lead to more targeted therapeutic and preventive measures. Monitoring transmission dynamics more efficiently holds promise to serve as a valuable tool in preventing in particular BSI caused by nosocomial pathogens.

Integrated genomic analysis of antibiotic resistance and virulence determinants in invasive strains of Streptococcus pneumoniae

Introduction

Streptococcus pneumoniae is an important human pathogen that may cause severe invasive pneumococcal diseases (IPDs) in young children and the elderly. A comprehensive comparative whole-genome analysis of invasive and non-invasive serotype strains offers great insights that are applicable to vaccine development and disease control.

Methods

In this study, 58 invasive (strains isolated from sterile sites) and 71 non-invasive (serotypes that have not been identified as invasive in our study) pneumococcal isolates were identified among the 756 pneumococcal isolates obtained from seven hospitals in Zhejiang, China (2010-2022). Serotyping, antimicrobial resistance tests, and genomic analyses were conducted to characterize these strains.

Results and discussion

The three most invasive serotypes were 23F, 14, and 6B. The invasive pneumococcal isolates' respective resistance rates against penicillin, ceftriaxone, tetracycline, and erythromycin were 34.5%, 15.5%, 98.3%, and 94.7%. Whole-genome sequencing indicated that the predominant invasive clonal complexes were CC271, CC876, and CC81. The high rate of penicillin non-susceptible Streptococcus pneumoniae (PNSP) is related to the clonal distribution of resistance-conferring penicillin-binding proteins (PBP). Interestingly, we found a negative correlation between invasiveness and resistance in the invasive pneumococcal serotype strains, which might be due to the proclivity of certain serotypes to retain their β-lactam resistance. Moreover, the mutually exclusive nature of zmpC and rrgC+srtBCD suggests their intricate and potentially redundant roles in promoting the development of IPD. These findings reveal significant implications for pneumococcal vaccine development in China, potentially informing treatment strategies and measures to mitigate disease transmission.

Phylogenetic inference of pneumococcal transmission from cross-sectional data, a pilot study

Background: Inference on pneumococcal transmission has mostly relied on longitudinal studies which are costly and resource intensive. Therefore, we conducted a pilot study to test the ability to infer who infected whom from cross-sectional pneumococcal sequences using phylogenetic inference. Methods: Five suspected transmission pairs, for which there was epidemiological evidence of who infected whom, were selected from a household study. For each pair, Streptococcus pneumoniae full genomes were sequenced from nasopharyngeal swabs collected on the same day. The within-host genetic diversity of the pneumococcal population was used to infer the transmission direction and then cross-validated with the direction suggested by the epidemiological records. Results: The pneumococcal genomes clustered into the five households from which the samples were taken. The proportion of concordantly inferred transmission direction generally increased with increasing minimum genome fragment size and single nucleotide polymorphisms. We observed a larger proportion of unique polymorphic sites in the source bacterial population compared to that of the recipient in four of the five pairs, as expected in the case of a transmission bottleneck. The only pair that did not exhibit this effect was also the pair that had consistent discordant transmission direction compared to the epidemiological records suggesting potential misdirection as a result of false-negative sampling. Conclusions: This pilot provided support for further studies to test if the direction of pneumococcal transmission can be reliably inferred from cross-sectional samples if sequenced with sufficient depth and fragment length.

Genomic classification and antimicrobial resistance profiling of Streptococcus pneumoniae and Haemophilus influenza isolates associated with paediatric otitis media and upper respiratory infection

Acute otitis media (AOM) is the most common childhood bacterial infectious disease requiring antimicrobial therapy. Most cases of AOM are caused by translocation of Streptococcus pneumoniae or Haemophilus influenzae from the nasopharynx to the middle ear during an upper respiratory tract infection (URI). Ongoing genomic surveillance of these pathogens is important for vaccine design and tracking of emerging variants, as well as for monitoring patterns of antibiotic resistance to inform treatment strategies and stewardship.In this work, we examined the ability of a genomics-based workflow to determine microbiological and clinically relevant information from cultured bacterial isolates obtained from patients with AOM or an URI. We performed whole genome sequencing (WGS) and analysis of 148 bacterial isolates cultured from the nasopharynx (N = 124, 94 AOM and 30 URI) and ear (N = 24, all AOM) of 101 children aged 6-35 months presenting with AOM or an URI. We then performed WGS-based sequence typing and antimicrobial resistance profiling of each strain and compared results to those obtained from traditional microbiological phenotyping.WGS of clinical isolates resulted in 71 S. pneumoniae genomes and 76 H. influenzae genomes. Multilocus sequencing typing (MSLT) identified 33 sequence types for S. pneumoniae and 19 predicted serotypes including the most frequent serotypes 35B and 3. Genome analysis predicted 30% of S. pneumoniae isolates to have complete or intermediate penicillin resistance. AMR predictions for S. pneumoniae isolates had strong agreement with clinical susceptibility testing results for beta-lactam and non beta-lactam antibiotics, with a mean sensitivity of 93% (86-100%) and a mean specificity of 98% (94-100%). MLST identified 29 H. influenzae sequence types. Genome analysis identified beta-lactamase genes in 30% of H. influenzae strains, which was 100% in agreement with clinical beta-lactamase testing. We also identified a divergent highly antibiotic-resistant strain of S. pneumoniae, and found its closest sequenced strains, also isolated from nasopharyngeal samples from over 15 years ago.Ultimately, our work provides the groundwork for clinical WGS-based workflows to aid in detection and analysis of H. influenzae and S. pneumoniae isolates.

Pan-Genome-Wide Association Study of Serotype 19A Pneumococci Identifies Disease-Associated Genes

Despite the widespread implementation of pneumococcal vaccines, hypervirulent Streptococcus pneumoniae serotype 19A is endemic worldwide. It is still unclear whether specific genetic elements contribute to complex pathogenicity of serotype 19A isolates. We performed a large-scale pan-genome-wide association study (pan-GWAS) of 1,292 serotype 19A isolates sampled from patients with invasive disease and asymptomatic carriers. To address the underlying disease-associated genotypes, a comprehensive analysis using three methods (Scoary, a linear mixed model, and random forest) was performed to compare disease and carriage isolates to identify genes consistently associated with disease phenotype. By using three pan-GWAS methods, we found consensus on statistically significant associations between genotypes and disease phenotypes (disease or carriage), with a subset of 30 consistently significant disease-associated genes. The results of functional annotation revealed that these disease-associated genes had diverse predicted functions, including those that participated in mobile genetic elements, antibiotic resistance, virulence, and cellular metabolism. Our findings suggest the multifactorial pathogenicity nature of this hypervirulent serotype and provide important evidence for the design of novel protein-based vaccines to prevent and control pneumococcal disease. IMPORTANCE It is important to understand the genetic and pathogenic characteristics of S. pneumoniae serotype 19A, which may provide important information for the prevention and treatment of pneumococcal disease. This global large-sample pan-GWAS study has identified a subset of 30 consistently significant disease-associated genes that are involved in mobile genetic elements, antibiotic resistance, virulence, and cellular metabolism. These findings suggest the multifactorial pathogenicity nature of hypervirulent S. pneumoniae serotype 19A isolates and provide implications for the design of novel protein-based vaccines.

Accelerating bioinformatics implementation in public health

We have adopted an open bioinformatics ecosystem to address the challenges of bioinformatics implementation in public health laboratories (PHLs). Bioinformatics implementation for public health requires practitioners to undertake standardized bioinformatic analyses and generate reproducible, validated and auditable results. It is essential that data storage and analysis are scalable, portable and secure, and that implementation of bioinformatics fits within the operational constraints of the laboratory. We address these requirements using Terra, a web-based data analysis platform with a graphical user interface connecting users to bioinformatics analyses without the use of code. We have developed bioinformatics workflows for use with Terra that specifically meet the needs of public health practitioners. These Theiagen workflows perform genome assembly, quality control, and characterization, as well as construction of phylogeny for insights into genomic epidemiology. Additonally, these workflows use open-source containerized software and the WDL workflow language to ensure standardization and interoperability with other bioinformatics solutions, whilst being adaptable by the user. They are all open source and publicly available in Dockstore with the version-controlled code available in public GitHub repositories. They have been written to generate outputs in standardized file formats to allow for further downstream analysis and visualization with separate genomic epidemiology software. Testament to this solution meeting the requirements for bioinformatic implementation in public health, Theiagen workflows have collectively been used for over 5 million sample analyses in the last 2 years by over 90 public health laboratories in at least 40 different countries. Continued adoption of technological innovations and development of further workflows will ensure that this ecosystem continues to benefit PHLs.

The Molecular Approaches and Challenges of Streptococcus pneumoniae Serotyping for Epidemiological Surveillance in the Vaccine Era

Streptococcus pneumoniae (pneumococcus) belongs to the Gram-positive cocci. This bacterium typically colonizes the nasopharyngeal region of healthy individuals. It has a distinct polysaccharide capsule - a virulence factor allowing the bacteria to elude the immune defense mechanisms. Consequently, it might trigger aggressive conditions like septicemia and meningitis in immunocompromised or older individuals. Moreover, children below five years of age are at risk of morbidity and mortality. Studies have found 101 S. pneumoniae capsular serotypes, of which several correlate with clinical and carriage isolates with distinct disease aggressiveness. Introducing pneumococcal conjugate vaccines (PCV) targets the most common disease-associated serotypes. Nevertheless, vaccine selection pressure leads to replacing the formerly dominant vaccine serotypes (VTs) by non-vaccine types (NVTs). Therefore, serotyping must be conducted for epidemiological surveillance and vaccine assessment. Serotyping can be performed using numerous techniques, either by the conventional antisera-based (Quellung and latex agglutination) or molecular-based approaches (sequetyping, multiplex PCR, real-time PCR, and PCR-RFLP). A cost-effective and practical approach must be used to enhance serotyping accuracy to monitor the prevalence of VTs and NVTs. Therefore, dependable pneumococcal serotyping techniques are essential to precisely monitor virulent lineages, NVT emergence, and genetic associations of isolates. This review discusses the principles, associated benefits, and drawbacks of the respective available conventional and molecular approaches, and potentially the whole genome sequencing (WGS) to be directed for future exploration.

Comparative genomic epidemiology of serotype 3 IPD and carriage isolates from Southampton, UK between 2005 and 2017

Serotype 3 pneumococci remains a significant cause of disease despite its inclusion in PCV13. Whilst clonal complex 180 (CC180) represents the major clone, recent studies have refined the population structure into three clades: Iα, Iβ and II, with the last being a recent divergent and more antibiotic-resistant. We present a genomic analysis of serotype 3 isolates from paediatric carriage and all-age invasive disease, collected between 2005 and 2017 in Southampton, UK. Forty-one isolates were available for analysis. Eighteen were isolated during the annual cross-sectional surveillance of paediatric pneumococcal carriage. The remaining 23 were isolated from blood/cerebrospinal fluid specimens at the University Hospital Southampton NHS Foundation Trust laboratory. All carriage isolates were CC180 GPSC12. Greater diversity was seen with invasive pneumococcal disease (IPD) with three GPSC83 (ST1377: n=2, ST260: n=1) and one GPSC3 (ST1716). For both carriage and IPD, Clade Iα was dominant (94.4 and 73.9 % respectively). Two isolates were Clade II with one from carriage (a 34-month-old, October 2017) and one invasive isolate (49-year-old, August 2015). Four IPD isolates were outside the CC180 clade. All isolates were genotypically susceptible to penicillin, erythromycin, tetracycline, co-trimoxazole and chloramphenicol. Two isolates (one each from carriage and IPD; both CC180 GPSC12) were phenotypically resistant to erythromycin and tetracycline; the IPD isolate was also resistant to oxacillin.In the Southampton area, carriage and invasive disease associated with serotype 3 is predominantly caused by Clade Iα CC180 GPSC12.

Meningitis-associated pneumococcal serotype 8, ST 53, strain is hypervirulent in a rat model and has non-haemolytic pneumolysin which can be attenuated by liposomes

Introduction

Streptococcus pneumoniae bacteria cause life-threatening invasive pneumococcal disease (IPD), including meningitis. Pneumococci are classified into serotypes, determined by differences in capsular polysaccharide and both serotype and pneumolysin toxin are associated with disease severity. Strains of serotype 8, ST 53, are increasing in prevalence in IPD in several countries.

Methods

Here we tested the virulence of such an isolate in a rat model of meningitis in comparison with a serotype 15B and a serotype 14 isolate. All three were isolated from meningitis patients in South Africa in 2019, where serotype 8 is currently the most common serotype in IPD.

Results and Discussion

Only the serotype 8 isolate was hypervirulent causing brain injury and a high mortality rate. It induced a greater inflammatory cytokine response than either the serotype 15B or 14 strain in the rat model and from primary mixed-glia cells isolated from mouse brains. It had the thickest capsule of the three strains and produced non-haemolytic pneumolysin. Pneumolysin-sequestering liposomes reduced the neuroinflammatory cytokine response in vitro indicating that liposomes have the potential to be an effective adjuvant therapy even for hypervirulent pneumococcal strains with non-haemolytic pneumolysin.

Predicting β-lactam susceptibility from the genome of Streptococcus pneumoniae and other mitis group streptococci

Introduction

For Streptococcus pneumoniae, β-lactam susceptibility can be predicted from the amino acid sequence of the penicillin-binding proteins PBP1a, PBP2b, and PBP2x. The combination of PBP-subtypes provides a PBP-profile, which correlates to a phenotypic minimal inhibitory concentration (MIC). The non-S. pneumoniae Mitis-group streptococci (MGS) have similar PBPs and exchange pbp-alleles with S. pneumoniae. We studied whether a simple BLAST analysis could be used to predict phenotypic susceptibility in Danish S. pneumoniae isolates and in internationally collected MGS.

Method

Isolates with available WGS and phenotypic susceptibility data were included. For each isolate, the best matching PBP-profile was identified by BLAST analysis. The corresponding MICs for penicillin and ceftriaxone was retrieved. Category agreement (CA), minor-, major-, and very major discrepancy was calculated. Genotypic-phenotypic accuracy was examined with Deming regression.

Results

Among 88 S. pneumoniae isolates, 55 isolates had a recognized PBP-profile, and CA was 100% for penicillin and 98.2% for ceftriaxone. In 33 S. pneumoniae isolates with a new PBP-profile, CA was 90.9% (penicillin) and 93.8% (ceftriaxone) using the nearest recognized PBP-profile. Applying the S. pneumoniae database to non-S. pneumoniae MGS revealed that none had a recognized PBP-profile. For Streptococcus pseudopneumoniae, CA was 100% for penicillin and ceftriaxone in 19 susceptible isolates. In 33 Streptococcus mitis isolates, CA was 75.8% (penicillin) and 86.2% (ceftriaxone) and in 25 Streptococcus oralis isolates CA was 8% (penicillin) and 100% (ceftriaxone).

Conclusion

Using a simple BLAST analysis, genotypic susceptibility prediction was accurate in Danish S. pneumoniae isolates, particularly in isolates with recognized PBP-profiles. Susceptibility was poorly predicted in other MGS using the current database.

PneumoKITy: A fast, flexible, specific, and sensitive tool for Streptococcus pneumoniae serotype screening and mixed serotype detection from genome sequence data

Determination of serotypes of Streptococcus pneumoniae is essential for monitoring current vaccine programmes. Since October 2017, pneumococcal serotypes in England have been derived from whole genome sequencing (WGS) data using our bioinformatic tool PneumoCaT. That tool was designed for serotype determination from pure cultures in a reference laboratory. To help determine multiple serotypes in pneumococcal carriage samples, we developed a new software tool named PneumoKITy (Pneumococcal K-mer Integrated Typing) that uses the powerful Mash k-mer screening method for pneumococcal serotyping. Mash k-mer screening is more sequence specific and much faster than the mapping method used in PneumoCaT and can determine 54 (58.1  %) of the 93 serotypes in the SSI Diagnostica phenotypical serotyping scheme to type level with the remainder called to serogroup or subgroup level (e.g., 11A/D). PneumoKITy can be run on both FastQ and assembly input, requiring up to 11× less memory and running up to 29× faster than the current version of PneumoCaT (1.2.1) on FastQ files. PneumoKITy can be used as a rapid, flexible serotype screening method which adds sensitive detection of mixed serotypes, e.g., for nasopharyngeal carriage studies where the presence of multiple serotypes is common. PneumoKITy's ability to function from assembly file, for pure culture serotype detection, increases its speed. This speed potentially enables the software to be run using low infrastructure overhead via web-based platforms. PneumoKITy could be used as a fast initial screening method with other tools used for those serotypes that could not be fully determined to type level if necessary. PneumoKITy was found to be highly accurate and sensitive when run on a panel of FastQ files derived from mixed cultures with all serotypes in 47/51 (92.2  %) of samples being accurately detected. PneumoKITy was also able to accurately estimate the relative abundance of serotypes in the same sample. Estimates being within a mean relative abundance of 1.5 % of the expected abundance in mixtures with known concentrations. PneumoKITy was able to detect minor serotypes with expected abundance of 1 % in the known mixture serotypes. PneumoKITy is a rapid, flexible tool with wide-ranging applications outside of the pure-culture, reference laboratory serotyping remit of PneumoCaT.

Characterization of Emerging Serotype 19A Pneumococcal Strains in Invasive Disease and Carriage, Belgium

After switching from 13-valent to 10-valent pneumococcal conjugate vaccine (PCV10) (2015-2016) for children in Belgium, we observed rapid reemergence of serotype 19A invasive pneumococcal disease (IPD). Whole-genome sequencing of 166 serotype 19A IPD isolates from children (n = 54) and older adults (n = 56) and carriage isolates from healthy children (n = 56) collected after the vaccine switch (2017-2018) showed 24 sequence types (STs). ST416 (global pneumococcal sequence cluster [GPSC] 4) and ST994 (GPSC146) accounted for 75.9% of IPD strains from children and 65.7% of IPD (children and older adults) and carriage isolates in the PCV10 period (2017-2018). These STs differed from predominant 19A IPD STs after introduction of PCV7 (2011) in Belgium (ST193 [GPSC11] and ST276 [GPSC10]), which indicates that prediction of emerging strains cannot be based solely on historical emerging strains. Despite their susceptible antimicrobial drug profiles, these clones spread in carriage and IPD during PCV10 use.

Distinct Streptococcus pneumoniae cause invasive disease in Papua New Guinea

Streptococcus pneumoniae is a key contributor to childhood morbidity and mortality in Papua New Guinea (PNG). For the first time, whole genome sequencing of 174 isolates has enabled detailed characterisation of diverse S. pneumoniae causing invasive disease in young children in PNG, 1989-2014. This study captures the baseline S. pneumoniae population prior to the introduction of 13-valent pneumococcal conjugate vaccine (PCV13) into the national childhood immunisation programme in 2014. Relationships amongst lineages, serotypes and antimicrobial resistance traits were characterised, and the population was viewed in the context of a global collection of isolates. The analyses highlighted adiverse S. pneumoniae population associated with invasive disease in PNG, with 45 unique Global Pneumococcal Sequence Clusters (GPSCs) observed amongst the 174 isolates reflecting multiple lineages observed in PNG that have not been identified in other geographic locations. The majority of isolates were from children with meningitis, of which 52% (n=72) expressed non-PCV13 serotypes. Over a third of isolates were predicted to be resistant to at least one antimicrobial. PCV13 serotype isolates had 10.1 times the odds of being multidrug-resistant (MDR) compared to non-vaccine serotype isolates, and no isolates with GPSCs unique to PNG were MDR. Serotype 2 was the most commonly identified serotype; we identified a highly clonal cluster of serotype 2 isolates unique to PNG, and a distinct second cluster indicative of long-distance transmission. Ongoing surveillance, including whole-genome sequencing, is needed to ascertain the impact of the national PCV13 programme upon the S. pneumoniae population, including serotype replacement and antimicrobial resistance traits.

Validation of Fourier Transform Infrared Spectroscopy for Serotyping of Streptococcus pneumoniae

Fourier transform infrared (FT-IR) spectroscopy (IR Biotyper; Bruker) allows highly discriminatory fingerprinting of closely related bacterial strains. In this study, FT-IR spectroscopy-based capsular typing of Streptococcus pneumoniae was validated as a rapid, cost-effective, and medium-throughput alternative to the classical phenotypic techniques. A training set of 233 strains was defined, comprising 34 different serotypes and including all 24 vaccine types (VTs) and 10 non-vaccine types (NVTs). The acquired spectra were used to (i) create a dendrogram where strains clustered together according to their serotypes and (ii) train an artificial neural network (ANN) model to predict unknown pneumococcal serotypes. During validation using 153 additional strains, we reached 98.0% accuracy for determining serotypes represented in the training set. Next, the performance of the IR Biotyper was assessed using 124 strains representing 59 non-training set serotypes. In this setting, 42 of 59 serotypes (71.1%) could be accurately categorized as being non-training set serotypes. Furthermore, it was observed that comparability of spectra was affected by the source of the Columbia medium used to grow the pneumococci and that this complicated the robustness and standardization potential of FT-IR spectroscopy. A rigorous laboratory workflow in combination with specific ANN models that account for environmental noise parameters can be applied to overcome this issue in the near future. The IR Biotyper has the potential to be used as a fast, cost-effective, and accurate phenotypic serotyping tool for S. pneumoniae.

Genetic background of Cambodian pneumococcal carriage isolates following pneumococcal conjugate vaccine 13

Streptococcus pneumoniae (the pneumococcus) is a leading cause of childhood mortality globally and in Cambodia. It is commensal in the human nasopharynx, occasionally resulting in invasive disease. Monitoring population genetic shifts, characterized by lineage and serotype expansions, as well as antimicrobial-resistance (AMR) patterns is crucial for assessing and predicting the impact of vaccination campaigns. We sought to elucidate the genetic background (global pneumococcal sequence clusters; GPSCs) of pneumococci carried by Cambodian children following perturbation by pneumococcal conjugate vaccine (PCV) 13. We sequenced pre-PCV13 (01/2013-12/2015, N=258) and post-PCV13 carriage isolates (01/2016-02/2017, N=428) and used PopPUNK and SeroBA to determine lineage prevalence and serotype composition. Following PCV13 implementation in Cambodia, we saw expansions of non-vaccine type (NVT) serotypes 23A (GPSC626), 34 (GPSC45) and 6D (GPSC16). We predicted antimicrobial susceptibility using the CDC-AMR pipeline and determined concordance with phenotypic data. The CDC-AMR pipeline had >90 % concordance with the phenotypic antimicrobial-susceptibility testing. We detected a high prevalence of AMR in both expanding non-vaccine serotypes and residual vaccine serotype 6B. Persistently high levels of AMR, specifically persisting multidrug-resistant lineages, warrant concern. The implementation of PCV13 in Cambodia has resulted in NVT serotype expansion reflected in the carriage population and driven by specific genetic backgrounds. Continued monitoring of these GPSCs during the ongoing collection of additional carriage isolates in this population is necessary.

Population genomics of pneumococcal carriage in South Africa following the introduction of the 13-valent pneumococcal conjugate vaccine (PCV13) immunization

Streptococcus pneumoniae is a major human pathogen responsible for over 317000 deaths in children <5 years of age with the burden of the disease being highest in low- and middle-income countries including South Africa. Following the introduction of the 7-valent and 13-valent pneumococcal conjugate vaccine (PCV) in South Africa in 2009 and 2011, respectively, a decrease in both invasive pneumococcal infections and asymptomatic carriage of vaccine-type pneumococci were reported. In this study, we described the changing epidemiology of the pneumococcal carriage population in South Africa, by sequencing the genomes of 1825 isolates collected between 2009 and 2013. Using these genomic data, we reported the changes in serotypes, Global Pneumococcal Sequence Clusters (GPSCs), and antibiotic resistance before and after the introduction of PCV13. The pneumococcal carriage population in South Africa has a high level of diversity, comprising of 126 GPSCs and 49 serotypes. Of the ten most prevalent GPSCs detected, six were predominantly found in Africa (GPSC22, GPSC21, GPSC17, GPSC33, GPSC34 and GPSC52). We found a significant decrease in PCV7 serotypes (19F, 6B, 23F and 14) and an increase in non-vaccine serotypes (NVT) (16F, 34, 35B and 11A) among children <2 years of age. The increase in NVTs was driven by pneumococcal lineages GPSC33, GPSC34, GPSC5 and GPSC22. Overall, a decrease in antibiotic resistance for 11 antimicrobials was detected in the PCV13 era. Further, we reported a higher resistance prevalence among vaccine types (VTs), as compared to NVTs; however, an increase in penicillin resistance among NVT was observed between the PCV7 and PCV13 eras. The carriage isolates from South Africa predominantly belonged to pneumococcal lineages, which are endemic to Africa. While the introduction of PCV resulted in an overall reduction of resistance in pneumococcal carriage isolates, an increase in penicillin resistance among NVTs was detected in children aged between 3 and 5 years, driven by the expansion of penicillin-resistant clones associated with NVTs in the PCV13 era.

Effect of pneumococcal conjugate vaccine availability on Streptococcus pneumoniae infections and genetic recombination in Zhejiang, China from 2009 to 2019

Pneumococcal pneumonia is one of the main reasons for child death worldwide. Pneumococcal conjugate vaccines (PCVs) are considered the most effective strategy for pneumococcal disease (PD) prevention, but how a pause in PCV vaccination affects the prevalence of PD or the genetic evolution of Streptococcus pneumoniae genetic evolution is unknown. Based on the unique PCV introduction timeline (vaccine unavailable during April 2015-April 2017) in China, we aimed to evaluate the effect of interrupted PCV availability on PD and pneumococcal genome variation. Pneumococcal isolates (n = 386) were collected retrospectively from eight sites in Zhejiang, China from 2009 to 2019 in which 184 pathogenic (isolates from sterile and infection sites) strains were identified. An interrupted time series analysis was conducted to estimate changes in PD and the recombination frequency of whole genome-sequenced strains was estimated via SNP calling. We found that both PD and pneumococcal genome variation were affected by interrupted PCV availability. The proportion (∼70%) of vaccine-type pneumococcal LRTI (VT-LRTI) in all LRTI cases decreased to ∼30% in the later PCV7 period and rebounded to ∼70% in children once PCV7 became unavailable in April 2015 (p = 0.0007). The major clone CC271 strains showed slowed (p = 0.0293) recombination frequency (decreased from 2.82 ± 1.16–0.72 ± 0.21) upon PCV removal. Our study illustrated for the first time that VT-LRTI fluctuated upon interrupted vaccine availability in China and causing a decreased of recombination frequency of vaccine types. Promoting a nationwide continuous vaccination programme and strengthening S. pneumoniae molecular epidemiology surveillance are essential for PD prevention.

GBS-SBG - GBS Serotyping by Genome Sequencing

Group B Streptococcus (GBS; Streptococcus agalactiae) is the most common cause of neonatal meningitis and a rising cause of sepsis in adults. Recently, it has also been shown to cause foodborne disease. As with many other bacteria, the polysaccharide capsule of GBS is antigenic, enabling its use for strain serotyping. Recent advances in DNA sequencing have made sequence-based typing attractive (as has been implemented for several other bacteria, including Escherichia coli, Klebsiella pneumoniae species complex, Streptococcus pyogenes, and others). For GBS, existing WGS-based serotyping systems do not provide complete coverage of all known GBS serotypes (specifically including subtypes of serotype III), and none are simultaneously compatible with the two most common data types, raw short reads and assembled sequences. Here, we create a serotyping database (GBS-SBG, GBS Serotyping by Genome Sequencing), with associated scripts and running instructions, that can be used to call all currently described GBS serotypes, including subtypes of serotype III, using both direct short-read- and assembly-based typing. We achieved higher concordance using GBS-SBG on a previously reported data set of 790 strains. We further validated GBS-SBG on a new set of 572 strains, achieving 99.8% concordance with PCR-based molecular serotyping using either short-read- or assembly-based typing. The GBS-SBG package is publicly available and will hopefully accelerate and simplify serotyping by sequencing for GBS.

Genomic surveillance of invasive Streptococcus pneumoniae isolates in the period pre-PCV10 and post-PCV10 introduction in Brazil

In 2010, Brazil introduced the 10-valent pneumococcal conjugate vaccine (PCV10) into the national children’s immunization programme. This study describes the genetic characteristics of invasive Streptococcus pneumoniae isolates before and after PCV10 introduction. A subset of 466 [pre-PCV10 (2008–2009): n=232, post-PCV10 (2012–2013): n=234;<5 years old: n=310, ≥5 years old: n=156] pneumococcal isolates, collected through national laboratory surveillance, were whole-genome sequenced (WGS) to determine serotype, pilus locus, antimicrobial resistance and genetic lineages. Following PCV10 introduction, in the <5 years age group, non-vaccine serotypes (NVT) serotype 3 and serotype 19A were the most frequent, and serotypes 12F, 8 and 9 N in the ≥5 years old group. The study identified 65 Global Pneumococcal Sequence Clusters (GPSCs): 49 (88 %) were GPSCs previously described and 16 (12 %) were Brazilian clusters. In total, 36 GPSCs (55 %) were NVT lineages, 18 (28 %) vaccine serotypes (VT) and 11 (17 %) were both VT and NVT lineages. In both sampling periods, the most frequent lineage was GPSC6 (CC156, serotypes 14/9V). In the <5 years old group, a decrease in penicillin (P=0.0123) and cotrimoxazole (P<0.0001) resistance and an increase in tetracycline (P=0.019) were observed. Penicillin nonsusceptibility was predicted in 40 % of the isolates; 127 PBP combinations were identified (51 predicted MIC≥0.125 mg l⁻¹); cotrimoxazole (folA and/or folP alterations), macrolide (mef and/or ermB) and tetracycline (tetM, tetO or tetS/M) resistance were predicted in 63, 13 and 21.6 % of pneumococci studied, respectively. The main lineages associated with multidrug resistance in the post-PCV10 period were composed of NVT, GPSC1 (CC320, serotype 19A), and GPSC47 (ST386, serotype 6C). The study provides a baseline for future comparisons and identified important NVT lineages in the post-PCV10 period in Brazil.

Streptococcus pneumoniae serotypes that frequently colonise the human nasopharynx are common recipients of penicillin-binding protein gene fragments from Streptococcus mitis

Streptococcus pneumoniae is an important global pathogen that causes bacterial pneumonia, sepsis and meningitis. Beta-lactam antibiotics are the first-line treatment for pneumococcal disease, however, their effectiveness is hampered by beta-lactam resistance facilitated by horizontal genetic transfer (HGT) with closely related species. Although interspecies HGT is known to occur among the species of the genus Streptococcus, the rates and effects of HGT between Streptococcus pneumoniae and its close relatives involving the penicillin binding protein (pbp) genes remain poorly understood. Here we applied the fastGEAR tool to investigate interspecies HGT in pbp genes using a global collection of whole-genome sequences of Streptococcus mitis, Streptococcus oralis and S. pneumoniae. With these data, we established that pneumococcal serotypes 6A, 13, 14, 16F, 19A, 19F, 23F and 35B were the highest-ranking serotypes with acquired pbp fragments. S. mitis was a more frequent pneumococcal donor of pbp fragments and a source of higher pbp nucleotide diversity when compared with S. oralis. Pneumococci that acquired pbp fragments were associated with a higher minimum inhibitory concentration (MIC) for penicillin compared with pneumococci without acquired fragments. Together these data indicate that S. mitis contributes to reduced β-lactam susceptibility among commonly carried pneumococcal serotypes that are associated with long carriage duration and high recombination frequencies. As pneumococcal vaccine programmes mature, placing increasing pressure on the pneumococcal population structure, it will be important to monitor the influence of antimicrobial resistance HGT from commensal streptococci such as S. mitis.

VPsero: Rapid Serotyping of Vibrio parahaemolyticus Using Serogroup-Specific Genes Based on Whole-Genome Sequencing Data

Vibrio parahaemolyticus has emerged as a significant enteropathogen in human and marine habitats worldwide, notably in regions where aquaculture products constitute a major nutritional source. It is a growing cause of diseases including gastroenteritis, wound infections, and septicemia. Serotyping assays use commercially available antisera to identify V. parahaemolyticus strains, but this approach is limited by high costs, complicated procedures, cross-immunoreactivity, and often subjective interpretation. By leveraging high-throughput sequencing technologies, we developed an in silico method based on comparison of gene clusters for lipopolysaccharide (LPSgc) and capsular polysaccharide (CPSgc) by firstly using the unique-gene strategy. The algorithm, VPsero, which exploits serogroup-specific genes as markers, covers 43 K and all 12 O serogroups in serotyping assays. VPsero is capable of predicting serotypes from assembled draft genomes, outputting LPSgc/CPSgc sequences, and recognizing possible novel serogroups or populations. Our tool displays high specificity and sensitivity in prediction toward V. parahaemolyticus strains, with an average sensitivity in serogroup prediction of 0.910 for O and 0.961 for K serogroups and a corresponding average specificity of 0.990 for O and 0.998 for K serogroups.

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

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

The effect of recombination on the evolution of a population of Neisseria meningitidis

Neisseria meningitidis (the meningococcus) is a major human pathogen with a history of high invasive disease burden, particularly in sub-Saharan Africa. Our current understanding of the evolution of meningococcal genomes is limited by the rarity of large-scale genomic population studies and lack of in-depth investigation of the genomic events associated with routine pathogen transmission. Here, we fill this knowledge gap by a detailed analysis of 2839 meningococcal genomes obtained through a carriage study of over 50,000 samples collected systematically in Burkina Faso, West Africa, before, during, and after the serogroup A vaccine rollout, 2009-2012. Our findings indicate that the meningococcal genome is highly dynamic, with highly recombinant loci and frequent gene sharing across deeply separated lineages in a structured population. Furthermore, our findings illustrate how population structure can correlate with genome flexibility, as some lineages in Burkina Faso are orders of magnitude more recombinant than others. We also examine the effect of selection on the population, in particular how it is correlated with recombination. We find that recombination principally acts to prevent the accumulation of deleterious mutations, although we do also find an example of recombination acting to speed the adaptation of a gene. In general, we show the importance of recombination in the evolution of a geographically expansive population with deep population structure in a short timescale. This has important consequences for our ability to both foresee the outcomes of vaccination programs and, using surveillance data, predict when lineages of the meningococcus are likely to become a public health concern.

Global genomic pathogen surveillance to inform vaccine strategies: a decade-long expedition in pneumococcal genomics

Vaccines are powerful agents in infectious disease prevention but often designed to protect against some strains that are most likely to spread and cause diseases. Most vaccines do not succeed in eradicating the pathogen and thus allow the potential emergence of vaccine evading strains. As with most evolutionary processes, being able to capture all variations across the entire genome gives us the best chance of monitoring and understanding the processes of vaccine evasion. Genomics is being widely adopted as the optimum approach for pathogen surveillance with the potential for early and precise identification of high-risk strains. Given sufficient longitudinal data, genomics also has the potential to forecast the emergence of such strains enabling immediate or pre-emptive intervention. In this review, we consider the strengths and challenges for pathogen genomic surveillance using the experience of the Global Pneumococcal Sequencing (GPS) project as an early example. We highlight the multifaceted nature of genome data and recent advances in genome-based tools to extract useful information relevant to inform vaccine strategies and treatment options. We conclude with future perspectives for genomic pathogen surveillance.

Serotype Distribution of Streptococcus pneumoniae Isolates Causing Invasive and Non-Invasive Infections Using Whole-Genome Sequencing in Ethiopia

Background

In Ethiopia, pneumococcal conjugate vaccine 10 (PCV10) was introduced in 2011 in the national vaccination program. This study was aimed to assess serotype distribution of invasive and non-invasive Streptococcus pneumoniae isolates using whole-genome sequencing.

Methods

A hospital-based prospective study was conducted from 2018 to 2019 at Addis Ababa and Amhara region referral hospitals, from all patients. Clinical Samples were collected and initially cultured onto 5% sheep blood agar at 37°C in a 5% CO₂ atmosphere. Sequencing was done using the Illumina NextSeq 500 and SeroBA was used to predict serotypes from whole-genome sequencing raw data.

Results

Of the 57 S. pneumoniae isolates, there were 32 circulating serotypes. The most common serotypes were 15A/B/C (n=5, 8.8%), 6A (n=4, 7.0%), 10A/F (n=4, 7.0%), 23A (n=4, 7.0%) and 7C (n=3, 5.3%). The serotype coverage of PCV10 and PCV13 were 12.3% and 26.3% respectively. The most common invasive serotypes were 15A/B/C (n=5, 8.8%) and 6A (n=4, 7.0%), and non-invasive serotypes were 23A (n=4, 7.0%) and 10A/F (n=3, 5.3%). The most prevalent serotype obtained from PCV10 eligible children was 3 (n=2, 3.5%). The prevalent serotype obtained from PCV10 non-eligible patients were type 23A (n=4, 7%) and type 6A (n=3, 5.2%). The most common serotypes among children ≤18 years old were 10A/F, 7C, 35A/B, 16F, 19A, 3 and 38. However, the proportions of some non-vaccine serotypes (11A/B and 15A/B/C) were higher in adult patients.

Conclusion

In this study a shift in the distribution of non-vaccinated S. pneumoniae serotypes increases in the population, and PCV10 serotype coverage was reduced as compared to PCV13. Therefore, it is important to continue monitoring serotype changes among all patients in addition to assessing the impact and effectiveness brought by vaccines and provides a foundation for prevention strategies and vaccine policies.

Early Signals of Vaccine-driven Perturbation Seen in Pneumococcal Carriage Population Genomic Data

Background

Pneumococcal conjugate vaccines (PCVs) have reduced pneumococcal diseases globally. Pneumococcal genomic surveys elucidate PCV effects on population structure but are rarely conducted in low-income settings despite the high disease burden.

Methods

We undertook whole-genome sequencing (WGS) of 660 pneumococcal isolates collected through surveys from healthy carriers 2 years from 13-valent PCV (PCV13) introduction and 1 year after rollout in northern Malawi. We investigated changes in population structure, within-lineage serotype dynamics, serotype diversity, and frequency of antibiotic resistance (ABR) and accessory genes.

Results

In children <5 years of age, frequency and diversity of vaccine serotypes (VTs) decreased significantly post-PCV, but no significant changes occurred in persons ≥5 years of age. Clearance of VT serotypes was consistent across different genetic backgrounds (lineages). There was an increase of nonvaccine serotypes (NVTs)—namely 7C, 15B/C, and 23A—in children <5 years of age, but 28F increased in both age groups. While carriage rates have been recently shown to remain stable post-PCV due to replacement serotypes, there was no change in diversity of NVTs. Additionally, frequency of intermediate-penicillin-resistant lineages decreased post-PCV. Although frequency of ABR genes remained stable, other accessory genes, especially those associated with mobile genetic element and bacteriocins, showed changes in frequency post-PCV.

Conclusions

We demonstrate evidence of significant population restructuring post-PCV driven by decreasing frequency of vaccine serotypes and increasing frequency of few NVTs mainly in children under 5. Continued surveillance with WGS remains crucial to fully understand dynamics of the residual VTs and replacement NVT serotypes post-PCV.

Evaluation of Pneumococcal Serotyping of Nasopharyngeal-Carriage Isolates by Latex Agglutination, Whole-Genome Sequencing (PneumoCaT), and DNA Microarray in a High-Pneumococcal-Carriage-Prevalence Population in Malawi

Accurate assessment of the serotype distribution associated with pneumococcal colonization and disease is essential for evaluating and formulating pneumococcal vaccines and for informing vaccine policy. For this reason, we evaluated the concordance between pneumococcal serotyping results by latex agglutination, whole-genome sequencing (WGS) with PneumoCaT, and DNA microarray for samples from community carriage surveillance in Blantyre, Malawi. Nasopharyngeal swabs were collected according to WHO recommendations between 2015 and 2017 by using stratified random sampling among study populations. Participants included healthy children 3 to 6 years old (vaccinated with the 13-valent pneumococcal conjugate vaccine [PCV13] as part of the Expanded Program on Immunization [EPI]), healthy children 5 to 10 years old (age-ineligible for PCV13), and HIV-infected adults (18 to 40 years old) on antiretroviral therapy (ART). For phenotypic serotyping, we used a 13-valent latex kit (Statens Serum Institut [SSI], Denmark). For genomic serotyping, we applied the PneumoCaT pipeline to whole-genome sequence libraries. For molecular serotyping by microarray, we used the BUGS Bioscience Senti-SP microarray. A total of 1,347 samples were analyzed. Concordance was 90.7% (95% confidence interval [CI], 89.0 to 92.2%) between latex agglutination and PneumoCaT, 95.2% (95% CI, 93.9 to 96.3%) between latex agglutination and the microarray, and 96.6% (95% CI, 95.5 to 97.5%) between the microarray and PneumoCaT. By detecting additional vaccine serotype (VT) pneumococci carried at low relative abundances (median, 8%), the microarray increased VT detection by 31.5% over that by latex serotyping. To conclude, all three serotyping methods were highly concordant in identifying dominant serotypes. Latex serotyping is accurate in identifying vaccine serotypes and requires the least expertise and resources for field implementation and analysis. However, WGS, which adds population structure, and microarray, which adds multiple-serotype carriage, should be considered at regional reference laboratories for investigating the importance of vaccine serotypes at low relative abundances in transmission and disease.

Two multi-fragment recombination events resulted in the β-lactam-resistant serotype 11A-ST6521 related to Spain9V-ST156 pneumococcal clone spreading in south-western Europe, 2008 to 2016

BackgroundThe successful pneumococcal clone Spain9V-ST156 (PMEN3) is usually associated with vaccine serotypes 9V and 14.AimOur objective was to analyse the increase of a serotype 11A variant of PMEN3 as cause of invasive pneumococcal disease (IPD) in Spain and its spread in south-western Europe.MethodsWe conducted a prospective multicentre study of adult IPD in Spain (2008-16). Furthermore, a subset of 61 penicillin-resistant serotype 11A isolates from France, Italy, Portugal and Spain were subjected to whole genome sequencing (WGS) and compared with 238 genomes from the European Nucleotide Archive (ENA).ResultsAlthough the incidence of serotype 11A in IPD was stable, a clonal shift was detected from CC62 (penicillin-susceptible) to CC156 (penicillin-resistant). By WGS, three major 11A-CC156 lineages were identified, linked to ST156 (n = 5 isolates; France, Italy and Portugal), ST166 (n = 4 isolates; France and Portugal) and ST838/6521 (n = 52 isolates; France, Portugal and Spain). Acquisition of the 11A capsule allowed to escape vaccine effect. AP200 (11A-ST62) was the donor for ST156 and ST838/6521 but not for ST166. In-depth analysis of ST838/6521 lineage showed two multi-fragment recombination events including four and seven fragments from an 11A-ST62 and an NT-ST344 representative, respectively.ConclusionThe increase in penicillin-resistant serotype 11A IPD in Spain was linked to the spread of a vaccine escape PMEN3 recombinant clone. Several recombination events were observed in PMEN3 acquiring an 11A capsule. The most successful 11A-PMEN3 lineage spreading in south-western Europe appeared after two multi-fragment recombination events with representatives of two major pneumococcal clones (11A-ST62 and NT-ST344).

Bacterial genome-wide association study of hyper-virulent pneumococcal serotype 1 identifies genetic variation associated with neurotropism

Hyper-virulent Streptococcus pneumoniae serotype 1 strains are endemic in Sub-Saharan Africa and frequently cause lethal meningitis outbreaks. It remains unknown whether genetic variation in serotype 1 strains modulates tropism into cerebrospinal fluid to cause central nervous system (CNS) infections, particularly meningitis. Here, we address this question through a large-scale linear mixed model genome-wide association study of 909 African pneumococcal serotype 1 isolates collected from CNS and non-CNS human samples. By controlling for host age, geography, and strain population structure, we identify genome-wide statistically significant genotype-phenotype associations in surface-exposed choline-binding (P = 5.00 × 10-08) and helicase proteins (P = 1.32 × 10-06) important for invasion, immune evasion and pneumococcal tropism to CNS. The small effect sizes and negligible heritability indicated that causation of CNS infection requires multiple genetic and other factors reflecting a complex and polygenic aetiology. Our findings suggest that certain pathogen genetic variation modulate pneumococcal survival and tropism to CNS tissue, and therefore, virulence for meningitis.

Within-host microevolution of Streptococcus pneumoniae is rapid and adaptive during natural colonisation

Genomic evolution, transmission and pathogenesis of Streptococcus pneumoniae, an opportunistic human-adapted pathogen, is driven principally by nasopharyngeal carriage. However, little is known about genomic changes during natural colonisation. Here, we use whole-genome sequencing to investigate within-host microevolution of naturally carried pneumococci in ninety-eight infants intensively sampled sequentially from birth until twelve months in a high-carriage African setting. We show that neutral evolution and nucleotide substitution rates up to forty-fold faster than observed over longer timescales in S. pneumoniae and other bacteria drives high within-host pneumococcal genetic diversity. Highly divergent co-existing strain variants emerge during colonisation episodes through real-time intra-host homologous recombination while the rest are co-transmitted or acquired independently during multiple colonisation episodes. Genic and intergenic parallel evolution occur particularly in antibiotic resistance, immune evasion and epithelial adhesion genes. Our findings suggest that within-host microevolution is rapid and adaptive during natural colonisation.

Characterization of Streptococcus pneumoniae strains causing invasive infections using whole-genome sequencing

Purpose: antigenic and genetic characterization of Streptococcus pneumoniae strains isolated from patients with invasive forms of pneumococcal infection using whole-genome sequencing.Materials and Methods. The study was performed on 46 S. pneumoniae strains isolated during the PEHASus multicenter studies in 2015-2018. Sequencing was performed using Illumina protocols and equipment. The SPAdes, SeroBA, PneumoCaT software were used for data processing, as well as BIGSdb software (PubMLST.org).Results and Discussion. Whole-genome sequences of strains were obtained; the information was entered into the PubMLST database (id: 51080-51125). Ten (21%) strains were found to have serotype 3. Five (11%) strains belonged to serotype 19F and five to serogroup 6; two of them belonged to serotype 6A; one strain had 6B and 1 had 6BE serotype; 1 strain showed discordant result (6A or 6BE). Serotype 15B was identified in 3 (6.5%) strains. Serotypes 7F, 8, 9V, 14, 22F, 23F and 28A were identified in two strains each; serotypes 1, 4, 9N, 10C, 12F, 18C, 35F, 37 and 38 were found once. The proportion of strains with serotypes included in PCV13 and PPV23 vaccines was 65% and 80%, respectively. 36 sequence types were found in strains; out of them, 6 sequence types were found for the first time. A dominant sequence type or clone complexes could not be identified using multilocus sequence typing except for serotype 3 strains. The inability to identify clonal complexes is in congruence with the previously obtained data on the absence of S. pneumoniae clones associated with pneumococcal meningitis in Russia.Conclusion. The information about serotypes of S. pneumoniae causing invasive infections together with epidemiologic data about strain sources and vaccination allows us to evaluate the effectiveness of pneumococcal vaccines and provide information for improving the PCR-based routine serotyping.

Streptococcus pneumoniae Nasopharyngeal Carriage among PCV-10-Vaccinated HIV-1-Infected Children with Maintained Serological Memory in Ethiopia

Streptococcus pneumoniae (S. pneumoniae) vaccines have substantially reduced the burden of invasive pneumococcal diseases (IPDs) worldwide. Despite high coverage with S. pneumoniae vaccination, upper-respiratory-tract colonization by S. pneumoniae is still common. We assessed maintenance of serological responses to S. pneumoniae serotypes included in PCV-10 by ELISA in HIV-1-infected children (n = 50) and age-matched controls (n = 50) in Ethiopia. We isolated S. pneumoniae in nasopharyngeal swabs and determined S. pneumoniae serotype by whole genome sequencing (WGS). Comparable levels of S. pneumoniae serotype-specific IgG concentrations were detected in plasma of HIV-1-infected children and matched controls, with geometric mean concentrations (GMCs) consistently higher than the protective threshold for PCV-10 serotypes of 0.35 μg/mL. We isolated S. pneumoniae from 38 (out of 97) nasopharyngeal swabs, 25 from HIV-1-infected children and 13 from controls. WGS based serotyping revealed 22 known S. pneumoniae serotypes and 2 nontypeable (NT) isolates. Non-PCV-10 serotypes represented >90% of isolates. We showed that HIV-1-infected children and matched controls in Ethiopia carry a level of maintained serological memory to PCV-10 considered protective for IPDs. We identified a higher proportion of nasopharyngeal carriage with highly pathogenic S. pneumoniae non-PCV strains among HIV-1-infected children compared to controls.