The 7-valent pneumococcal conjugated vaccine (PCV7) has affected the genetic population of Streptococcus pneumoniae in pediatric carriage. Little is known however about pneumococcal population genomics in adult invasive pneumococcal disease (IPD) under vaccine pressure. We sequenced and serotyped 349 strains of S. pneumoniae isolated from IPD patients in Nijmegen between 2001 and 2011. Introduction of PCV7 in the Dutch National Immunization Program in 2006 preluded substantial alterations in the IPD population structure caused by serotype replacement. No evidence could be found for vaccine induced capsular switches. We observed that after a temporary bottleneck in gene diversity after the introduction of PCV7, the accessory gene pool re-expanded mainly by genes already circulating pre-PCV7. In the post-vaccine genomic population a number of genes changed frequency, certain genes became overrepresented in vaccine serotypes, while others shifted towards non-vaccine serotypes. Whether these dynamics in the invasive pneumococcal population have truly contributed to invasiveness and manifestations of disease remains to be further elucidated. We suggest the use of whole genome sequencing for surveillance of pneumococcal population dynamics that could give a prospect on the course of disease, facilitating effective prevention and management of IPD.

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

Predicting how pathogen populations will change over time is challenging. Such has been the case with Streptococcus pneumoniae, an important human pathogen, and the pneumococcal conjugate vaccines (PCVs), which target only a fraction of the strains in the population. Here, we use the frequencies of accessory genes to predict changes in the pneumococcal population after vaccination, hypothesizing that these frequencies reflect negative frequency-dependent selection (NFDS) on the gene products. We find that the standardized predicted fitness of a strain, estimated by an NFDS-based model at the time the vaccine is introduced, enables us to predict whether the strain increases or decreases in prevalence following vaccination. Further, we are able to forecast the equilibrium post-vaccine population composition and assess the invasion capacity of emerging lineages. Overall, we provide a method for predicting the impact of an intervention on pneumococcal populations with potential application to other bacterial pathogens in which NFDS is a driving force.

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.

A diverse set of mobile genetic elements (MGEs) transmit between Streptococcus pneumoniae cells, but many isolates remain uninfected. The best-characterised defences against horizontal transmission of MGEs are restriction-modification systems (RMSs), of which there are two phase-variable examples in S. pneumoniae. Additionally, the transformation machinery has been proposed to limit vertical transmission of chromosomally integrated MGEs. This work describes how these mechanisms can act in concert. Experimental data demonstrate RMS phase variation occurs at a sub-maximal rate. Simulations suggest this may be optimal if MGEs are sometimes vertically inherited, as it reduces the probability that an infected cell will switch between RMS variants while the MGE is invading the population, and thereby undermine the restriction barrier. Such vertically inherited MGEs can be deleted by transformation. The lack of between-strain transformation hotspots at known prophage att sites suggests transformation cannot remove an MGE from a strain in which it is fixed. However, simulations confirmed that transformation was nevertheless effective at preventing the spread of MGEs into a previously uninfected cell population, if a recombination barrier existed between co-colonising strains. Further simulations combining these effects of phase variable RMSs and transformation found they synergistically inhibited MGEs spreading, through limiting both vertical and horizontal transmission.

Serotype 1 Streptococcus pneumoniae is a leading cause of invasive pneumococcal disease (IPD) worldwide, with the highest burden in developing countries. We report the whole-genome sequencing analysis of 448 serotype 1 isolates from 27 countries worldwide (including 11 in Africa). The global serotype 1 population shows a strong phylogeographic structure at the continental level, and within Africa there is further region-specific structure. Our results demonstrate that region-specific diversification within Africa has been driven by limited cross-region transfer events, genetic recombination and antimicrobial selective pressure. Clonal replacement of the dominant serotype 1 clones circulating within regions is uncommon; however, here we report on the accessory gene content that has contributed to a rare clonal replacement event of ST3081 with ST618 as the dominant cause of IPD in the Gambia.

Streptococcus pneumoniae is common nasopharyngeal commensal bacterium and important human pathogen. Vaccines against a subset of pneumococcal antigenic diversity have reduced rates of disease, without changing the frequency of asymptomatic carriage, through altering the bacterial population structure. These changes can be studied in detail through using genome sequencing to characterise systematically-sampled collections of carried S. pneumoniae. This dataset consists of 616 annotated draft genomes of isolates collected from children during routine visits to primary care physicians in Massachusetts between 2001, shortly after the seven valent polysaccharide conjugate vaccine was introduced, and 2007. Also made available are a core genome alignment and phylogeny describing the overall population structure, clusters of orthologous protein sequences, software for inferring serotype from Illumina reads, and whole genome alignments for the analysis of closely-related sets of pneumococci. These data can be used to study both bacterial evolution and the epidemiology of a pathogen population under selection from vaccine-induced immunity.

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.

Streptococcus pneumoniae is responsible for 240 000-460 000 deaths in children under 5 years of age each year. Accurate identification of pneumococcal serotypes is important for tracking the distribution and evolution of serotypes following the introduction of effective vaccines. Recent efforts have been made to infer serotypes directly from genomic data but current software approaches are limited and do not scale well. Here, we introduce a novel method, SeroBA, which uses a k-mer approach. We compare SeroBA against real and simulated data and present results on the concordance and computational performance against a validation dataset, the robustness and scalability when analysing a large dataset, and the impact of varying the depth of coverage on sequence-based serotyping. SeroBA can predict serotypes, by identifying the cps locus, directly from raw whole genome sequencing read data with 98 % concordance using a k-mer-based method, can process 10 000 samples in just over 1 day using a standard server and can call serotypes at a coverage as low as 15-21×. SeroBA is implemented in Python3 and is freely available under an open source GPLv3 licence from: https://github.com/sanger-pathogens/seroba.

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.

Streptococcus pneumoniae is one of the most important causes of microbial diseases in humans. The genomes of 44 diverse strains of S. pneumoniae were analyzed and compared with strains of non-pathogenic streptococci of the Mitis group.

Streptococcus pneumoniae is a leading cause of meningitis in countries where pneumococcal conjugate vaccines (PCV) targeting commonly occurring serotypes are not routinely used. However, effectiveness of PCV would be jeopardized by emergence of invasive pneumococcal diseases (IPD) caused by serotypes which are not included in PCV. Systematic hospital based surveillance in Bangladesh was established and progressively improved to determine the pathogens causing childhood sepsis and meningitis. This also provided the foundation for determining the spectrum of serotypes causing IPD. This article reports an unprecedented upsurge of serotype 2, an uncommon pneumococcal serotype, without any known intervention.

Carriage of either single or multiple pneumococcal serotypes (multiple carriage) is a prerequisite for developing invasive pneumococcal disease. However, despite the reported high rates of pneumococcal carriage in Malawi, no data on carriage of multiple serotypes has been reported previously. Our study provides the first description of the prevalence of multiple pneumococcal carriage in Malawi.

Whole genome sequencing (WGS)-based approaches for pneumococcal capsular typing have become an alternative to serological methods. In silico serotyping from WGS has not yet been applied to long-read sequences produced by third-generation technologies. The objective of the study was to determine the capsular types of pneumococci causing invasive disease in Catalonia (Spain) using serological typing and WGS and to compare the performance of different bioinformatics pipelines using short- and long-read data from WGS. All invasive pneumococcal pediatric isolates collected in Hospital Sant Joan de Déu (Barcelona) from 2013 to 2019 were included. Isolates were assigned a capsular type by serological testing based on anticapsular antisera and by different WGS-based pipelines: Illumina sequencing followed by serotyping with PneumoCaT, SeroBA, and Pathogenwatch vs MinION-ONT sequencing coupled with serotyping by Pathogenwatch from pneumococcal assembled genomes. A total of 119 out of 121 pneumococcal isolates were available for sequencing. Twenty-nine different serotypes were identified by serological typing, with 24F (n = 17; 14.3%), 14 (n = 10; 8.4%), and 15B/C (n = 8; 6.7%) being the most common serotypes. WGS-based pipelines showed initial concordance with serological typing (>91% of accuracy). The main discrepant results were found at the serotype level within a serogroup: 6A/B, 6C/D, 9A/V, 11A/D, and 18B/C. Only one discrepancy at the serogroup level was observed: serotype 29 by serological testing and serotype 35B/D by all WGS-based pipelines. Thus, bioinformatics WGS-based pipelines, including those using third-generation sequencing, are useful for pneumococcal capsular assignment. Possible discrepancies between serological typing and WGS-based approaches should be considered in pneumococcal capsular-type surveillance studies.

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

Streptococcus pneumoniae causes a high burden of invasive pneumococcal disease (IPD) globally, especially in children from resource-poor settings. Like many bacteria, the pneumococcus can import DNA from other strains or even species by transformation and homologous recombination, which has allowed the pneumococcus to evade clinical interventions such as antibiotics and pneumococcal conjugate vaccines (PCVs). Pneumococci are enclosed in a complex polysaccharide capsule that determines the serotype; the capsule varies in size and is associated with properties including carriage prevalence and virulence. We determined and quantified the association between capsule and recombination events using genomic data from a diverse collection of serotypes sampled in Malawi. We determined both the amount of variation introduced by recombination relative to mutation (the relative rate) and how many individual recombination events occur per isolate (the frequency). Using univariate analyses, we found an association between both recombination measures and multiple factors associated with the capsule, including duration and prevalence of carriage. Because many capsular factors are correlated, we used multivariate analysis to correct for collinearity. Capsule size and carriage duration remained positively associated with recombination, although with a reduced P value, and this effect may be mediated through some unassayed additional property associated with larger capsules. This work describes an important impact of serotype on recombination that has been previously overlooked. While the details of how this effect is achieved remain to be determined, it may have important consequences for the serotype-specific response to vaccines and other interventions.

The pneumococcal pilus has been shown to be an important determinant of adhesion and virulence in mouse models of colonization, pneumonia, and bacteremia. A pilus is capable of inducing protective immunity, supporting its inclusion in next-generation pneumococcal protein vaccine formulations. Whether this vaccine target is common among pneumococci in sub-Saharan Africa is uncertain. To define the prevalence and genetic diversity of type I and II pili among invasive pneumococci in Malawi prior to the introduction of the 13-valent pneumococcal conjugate vaccine (PCV13) into routine childhood immunization, we examined 188 Streptococcus pneumoniae isolates collected between 2002 and 2008 (17% serotype 1). In this region of high disease burden, we found a low frequency of invasive piliated pneumococci (14%) and pilus gene sequence diversity similar to that seen previously in multiple global pneumococcal lineages. All common serotypes with pilus were covered by PCV13 and so we predict that pilus prevalence will be reduced in the Malawian pneumococcal population after PCV13 introduction.

Despite enabling Streptococcus pneumoniae to acquire antibiotic resistance and evade vaccine-induced immunity, transformation occurs at variable rates across pneumococci. Phase variants of isolate RMV7, distinguished by altered methylation patterns driven by the translocating variable restriction-modification (tvr) locus, differed significantly in their transformation efficiencies and biofilm thicknesses. These differences were replicated when the corresponding tvr alleles were introduced into an RMV7 derivative lacking the locus. RNA-seq identified differential expression of the type 1 pilus, causing the variation in biofilm formation, and inhibition of competence induction in the less transformable variant, RMV7domi. This was partly attributable to RMV7domi's lower expression of ManLMN, which promoted competence induction through importing N-acetylglucosamine. This effect was potentiated by analogues of some proteobacterial competence regulatory machinery. Additionally, one of RMV7domi's phage-related chromosomal island was relatively active, which inhibited transformation by increasing expression of the stress response proteins ClpP and HrcA. However, HrcA increased competence induction in the other variant, with its effects depending on Ca2+ supplementation and heat shock. Hence the heterogeneity in transformation efficiency likely reflects the diverse signalling pathways by which it is affected. This regulatory complexity will modulate population-wide responses to synchronising quorum sensing signals to produce co-ordinated yet stochastic bet hedging behaviour.

This study aimed to define the molecular basis of co-trimoxazole resistance in Malawian pneumococci under the dual selective pressure of widespread co-trimoxazole and sulfadoxine/pyrimethamine use.

Streptococcus pneumoniae (pneumococcus) is a leading vaccine-preventable cause of childhood invasive disease. Nigeria has the second highest pneumococcal disease burden globally, with an estimated ~49 000 child deaths caused by pneumococcal infections each year. Ten-valent pneumococcal conjugate vaccine (GSK; PCV10) was introduced in December 2014 in a phased approach. However, few studies have characterized the disease-causing pneumococci from Nigeria. This study assessed the prevalence of serotypes, antibiotic susceptibility and genomic lineages using whole genome sequencing and identified lineages that could potentially escape PCV10 (GSK). We also investigated the potential differences in pneumococcal lineage features between children with and without sickle cell disease. A collection of 192 disease-causing pneumococcal isolates was obtained from Kano (n=189) and Abuja (n=3) states, Nigeria, between 1 January 2014 and 31 May 2018. The majority (99 %, 190/192) of specimens were recovered from children aged 5 years or under. Among them, 37 children had confirmed or traits of sickle cell disease. Our findings identified 25 serotypes expressed by 43 Global Pneumococcal Sequence Clusters (GPSCs) and 85 sequence types (STs). The most common serotypes were 14 (18 %, n=35), 6B (16 %, n=31), 1 (9 %, n=17), 5 (9 %, n=17) and 6A (9 %, n=17); all except serotype 6A are included in PCV10 (GSK). PCV10 (SII; PNEUMOSIL) and PCV13 formulations include serotypes 6A and 19A which would increase the overall coverage from 67 % by PCV10 (GSK) to 78 and 82 %, respectively. The pneumococcal lineages were a mix of globally spreading and unique local lineages. Following the use of PCV10 (GSK), GPSC5 expressing serotype 6A, GPSC10 (19A), GPSC26 (12F and 46) and GPSC627 (9L) are non-vaccine type lineages that could persist and potentially expand under vaccine-selective pressure. Approximately half (52 %, 99/192) of the pneumococcal isolates were resistant to the first-line antibiotic penicillin and 44 % (85/192) were multidrug-resistant. Erythromycin resistance was very low (2 %, 3/192). There was no significant difference in clinical manifestation, serotype prevalence or antibiotic resistance between children with and without traits of or confirmed sickle cell disease. In summary, our findings show that a high percentage of the pneumococcal disease were caused by the serotypes that are covered by currently available vaccines. Given the low prevalence of resistance, macrolide antibiotics, such as erythromycin, should be considered as an option to treat pneumococcal disease in Nigeria. However, appropriate use of macrolide antibiotics should be vigilantly monitored to prevent the potential increase in macrolide resistance.