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.

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.

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 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.

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.

Invasive disease caused by Streptococcus pneumoniae (IPD) is one of the leading causes of morbidity and mortality in young children worldwide. In Argentina, PCV13 was introduced into the childhood immunization programme nationwide in 2012 and PCV7 was available from 2000, but only in the private market. Since 1993 the National IPD Surveillance Programme, consisting of 150 hospitals, has conducted nationwide pneumococcal surveillance in Argentina in children under 6 years of age, as part of the SIREVA II-OPS network. A total of 1713 pneumococcal isolates characterized by serotype (Quellung) and antimicrobial resistance (agar dilution) to ten antibiotics, belonging to three study periods: pre-PCV7 era 1998-1999 (pre-PCV), before the introduction of PCV13 2010-2011 (PCV7) and after the introduction of PCV13 2012-2013 (PCV13), were available for inclusion. Fifty-four serotypes were identified in the entire collection and serotypes 14, 5 and 1 represented 50 % of the isolates. Resistance to penicillin was 34.9 %, cefotaxime 10.6 %, meropenem 4.9 %, cotrimoxazole 45 %, erythromycin 21.5 %, tetracycline 15.4 % and chloramphenicol 0.4 %. All the isolates were susceptible to levofloxacin, rifampin and vancomycin. Of 1713 isolates, 1061 (61.9 %) were non-susceptible to at least one antibiotic and 235(13.7 %) were multidrug resistant. A subset of 413 isolates was randomly selected and whole-genome sequenced as part of Global Pneumococcal Sequencing Project (GPS). The genome data was used to investigate the population structure of S. pneumoniae defining pneumococcal lineages using Global Pneumococcal Sequence Clusters (GPSCs), sequence types (STs) and clonal complexes (CCs), prevalent serotypes and their associated pneumococcal lineages and genomic inference of antimicrobial resistance. The collection showed a great diversity of strains. Among the 413 isolates, 73 known and 36 new STs were identified belonging to 38 CCs and 25 singletons, grouped into 52 GPSCs. Important changes were observed among vaccine types when pre-PCV and PCV13 periods were compared; a significant decrease in serotypes 14, 6B and 19F and a significant increase in 7F and 3. Among non-PCV13 types, serogroup 24 increased from 0 % in pre-PCV to 3.2 % in the PCV13 period. Our analysis showed that 66.1 % (273/413) of the isolates were predicted to be non-susceptible to at least one antibiotic and 11.9 % (49/413) were multidrug resistant. We found an agreement of 100 % when comparing the serotype determined by Quellung and WGS-based serotyping and 98.4 % of agreement in antimicrobial resistance. Continued surveillance of the pneumococcal population is needed to reveal the dynamics of pneumococcal isolates in Argentina in post-PCV13. This article contains data hosted by Microreact.

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-1); 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.

Globally, India has a high burden of pneumococcal disease, and pneumococcal conjugate vaccine (PCV) has been rolled out in different phases across the country since May 2017 in the national infant immunization programme (NIP). To provide a baseline for assessing the impact of the vaccine on circulating pneumococci in India, genetic characterization of pneumococcal isolates detected prior to introduction of PCV would be helpful. Here we present a population genomic study of 480 Streptococcus pneumoniae isolates collected across India and from all age groups before vaccine introduction (2009-2017), including 294 isolates from pneumococcal disease and 186 collected through nasopharyngeal surveys. Population genetic structure, serotype and antimicrobial susceptibility profile were characterized and predicted from whole-genome sequencing data. Our findings revealed high levels of genetic diversity represented by 110 Global Pneumococcal Sequence Clusters (GPSCs) and 54 serotypes. Serotype 19F and GPSC1 (CC320) was the most common serotype and pneumococcal lineage, respectively. Coverage of PCV13 (Pfizer) and 10-valent Pneumosil (Serum Institute of India) serotypes in age groups of ≤2 and 3-5 years were 63-75 % and 60-69 %, respectively. Coverage of PPV23 (Merck) serotypes in age groups of ≥50 years was 62 % (98/158). Among the top five lineages causing disease, GPSC10 (CC230), which ranked second, is the only lineage that expressed both PCV13 (serotypes 3, 6A, 14, 19A and 19F) and non-PCV13 (7B, 13, 10A, 11A, 13, 15B/C, 22F, 24F) serotypes. It exhibited multidrug resistance and was the largest contributor (17 %, 18/103) of NVTs in the disease-causing population. Overall, 42 % (202/480) of isolates were penicillin-resistant (minimum inhibitory concentration ≥0.12 µg ml-1) and 45 % (217/480) were multidrug-resistant. Nine GPSCs (GPSC1, 6, 9, 10, 13, 16, 43, 91, 376) were penicillin-resistant and among them six were multidrug-resistant. Pneumococci expressing PCV13 serotypes had a higher prevalence of antibiotic resistance. Sequencing of pneumococcal genomes has significantly improved our understanding of the biology of these bacteria. This study, describing the pneumococcal disease and carriage epidemiology pre-PCV introduction, demonstrates that 60-75 % of pneumococcal serotypes in children ≤5 years are covered by PCV13 and Pneumosil. Vaccination against pneumococci is very likely to reduce antibiotic resistance. A multidrug-resistant pneumococcal lineage, GPSC10 (CC230), is a high-risk clone that could mediate serotype replacement.

The introduction of pneumococcal conjugate vaccines (PCV7, PCV10, PCV13) around the world has proved successful in preventing invasive pneumococcal disease. However, immunization against Streptococcus pneumoniae has led to serotype replacement by non-vaccine serotypes, including serotype 15A. Clonal complex 63 (CC63) is associated with many serotypes and has been reported in association with 15A after introduction of PCVs. A total of 865 CC63 isolates were included in this study, from the USA (n=391) and a global collection (n=474) from 1998-2019 and 1995-2018, respectively. We analysed the genomic sequences to identify serotypes and penicillin-binding protein (PBP) genes 1A, 2B and 2X, and other resistance determinants, to predict minimum inhibitory concentrations (MICs) against penicillin, erythromycin, clindamycin, co-trimoxazole and tetracycline. We conducted phylogenetic and spatiotemporal analyses to understand the evolutionary history of the 15A-CC63 sub-lineage. Overall, most (89.5 %, n=247) pre-PCV isolates in the CC63 cluster belonged to serotype 14, with 15A representing 6.5 % of isolates. Conversely, serotype 14 isolates represented 28.2 % of post-PCV CC63 isolates (n=618), whilst serotype 15A isolates represented 65.4 %. Dating of the CC63 lineage determined the most recent common ancestor emerged in the 1980s, suggesting the 15A-CC63 sub-lineage emerged from its closest serotype 14 ancestor prior to the development of pneumococcal vaccines. This sub-lineage was predominant in the USA, Israel and China. Multidrug resistance (to three or more drug classes) was widespread among isolates in this sub-lineage. We show that the CC63 lineage is globally distributed and most of the isolates are penicillin non-susceptible, and thus should be monitored.

Bacteriocins are antimicrobial peptides produced by bacteria to inhibit other bacteria in the surrounding environment. Streptococcus pneumoniae is a leading cause of disease worldwide and colonises the healthy human nasopharynx, where it competes for space and nutrients. Pneumococcal conjugate vaccines have reduced the incidence of disease, but they also restructure the bacterial population, and this restructuring likely alters the nasopharyngeal competition dynamics. Here, the distribution of bacteriocins was examined in over 5000 carriage and disease-causing pneumococci from Iceland and Kenya, recovered before and after the introduction of pneumococcal vaccination. Overall, up to eleven different bacteriocin gene clusters were identified per pneumococcus. Significant differences in the prevalence of bacteriocins were observed before and after vaccine introduction, and among carriage and disease-causing pneumococci, which were largely explained by the bacterial population structure. Genetically similar pneumococci generally harboured the same bacteriocins although sometimes different repertoires of bacteriocins were observed, which suggested that horizontal transfer of bacteriocin clusters had occurred. These findings demonstrated that vaccine-mediated changes in the pneumococcal population altered the prevalence and distribution of bacteriocins. The consequences of this for pneumococcal colonisation and disease remain to be determined.

Mitis group Streptococcus are human obligate bacteria residing in the nasopharynx and oral cavity. They comprise both commensal and pathogenic species with the most well-known being Streptococcus pneumoniae - a leading cause of meningitis and pneumonia. A primary difference between the commensal and pathogenic species is the presence of the polysaccharide capsule - a major virulence factor in S. pneumoniae, also present in other commensal species. Our current understanding of the evolutionary divergence of the pathogenic and commensal species has been inferred from extant strains. Ancient genomes can further elucidate streptococcal evolutionary history. We extracted streptococcal genome reads from a 5700-year-old ancient metagenome and worked towards characterizing them. Due to excessive within- and between-species recombination common among streptococci we were unable to parse individual species. Further, the composite reads of the ancient metagenome do not fit within the diversity of any specific extant species. Using a capsular gene database and AT-content analysis we determined that this ancient metagenome is missing polysaccharide synthesis genes integral to streptococcal capsule formation. The presence of multiple zinc metalloproteases suggests that adaptation to host IgA1 had begun and the presence of other virulence factors further implies development of close host-microbe interactions, though the absence of a capsule suggests an inability to cause invasive disease. The presence of specific virulence factors such as pneumolysin implies stable maintenance of such genes through streptococcal evolution that may strengthen their value as anti-pneumococcal vaccine antigens, while maintaining awareness of their potential presence in commensal species. Following from Jensen et al.'s initial analysis we provide historical context for this long time human nasopharyngeal resident, the Mitis group Streptococcus.

Recent studies have provided evidence for rapid pathogen genome diversification, some of which could potentially affect the course of disease. We have previously described such variation seen between isolates infecting the blood and cerebrospinal fluid (CSF) of a single patient during a case of bacterial meningitis. Here, we performed whole-genome sequencing of paired isolates from the blood and CSF of 869 meningitis patients to determine whether such variation frequently occurs between these two niches in cases of bacterial meningitis. Using a combination of reference-free variant calling approaches, we show that no genetic adaptation occurs in either invaded niche during bacterial meningitis for two major pathogen species, Streptococcus pneumoniae and Neisseria meningitidis. This study therefore shows that the bacteria capable of causing meningitis are already able to do this upon entering the blood, and no further sequence change is necessary to cross the blood-brain barrier. Our findings place the focus back on bacterial evolution between nasopharyngeal carriage and invasion, or diversity of the host, as likely mechanisms for determining invasiveness.

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.

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.

The potential for genome-wide modelling of epistasis has recently surfaced given the possibility of sequencing densely sampled populations and the emerging families of statistical interaction models. Direct coupling analysis (DCA) has previously been shown to yield valuable predictions for single protein structures, and has recently been extended to genome-wide analysis of bacteria, identifying novel interactions in the co-evolution between resistance, virulence and core genome elements. However, earlier computational DCA methods have not been scalable to enable model fitting simultaneously to 104-105 polymorphisms, representing the amount of core genomic variation observed in analyses of many bacterial species. Here, we introduce a novel inference method (SuperDCA) that employs a new scoring principle, efficient parallelization, optimization and filtering on phylogenetic information to achieve scalability for up to 105 polymorphisms. Using two large population samples of Streptococcus pneumoniae, we demonstrate the ability of SuperDCA to make additional significant biological findings about this major human pathogen. We also show that our method can uncover signals of selection that are not detectable by genome-wide association analysis, even though our analysis does not require phenotypic measurements. SuperDCA, thus, holds considerable potential in building understanding about numerous organisms at a systems biological level.

The implementation of pneumococcal conjugate vaccines (PCVs) has led to a decline in vaccine-type disease. However, there is evidence that the epidemiology of non-typeable Haemophilus influenzae (NTHi) carriage and disease can be altered as a consequence of PCV introduction. We explored the epidemiological shifts in NTHi carriage using whole genome sequencing over a 5-year period that included PCV13 replacement of PCV7 in the UK's National Immunization Programme in 2010. Between 2008/09 and 2012/13 (October to March), nasopharyngeal swabs were taken from children <5 years of age. Significantly increased carriage post-PCV13 was observed and lineage-specific associations with Streptococcus pneumoniae were seen before but not after PCV13 introduction. NTHi were characterized into 11 discrete, temporally stable lineages, congruent with current knowledge regarding the clonality of NTHi. The increased carriage could not be linked to the expansion of a particular clone and different co-carriage dynamics were seen before PCV13 implementation when NTHi co-carried with vaccine serotype pneumococci. In summary, PCV13 introduction has been shown to have an indirect effect on NTHi epidemiology and there exists both negative and positive, distinct associations between pneumococci and NTHi. This should be considered when evaluating the impacts of pneumococcal vaccine design and policy.

The transmission dynamics of Streptococcus pneumoniae in sub-Saharan Africa are poorly understood due to a lack of adequate epidemiological and genomic data. Here we leverage a longitudinal cohort from 21 neighbouring villages in rural Africa to study how closely related strains of S. pneumoniae are shared among infants. We analysed 1074 pneumococcal genomes isolated from 102 infants from 21 villages. Strains were designated for unique serotype and sequence-type combinations, and we arbitrarily defined strain sharing where the pairwise genetic distance between strains could be accounted for by the mean within host intra-strain diversity. We used non-parametric statistical tests to assess the role of spatial distance and prolonged carriage on strain sharing using a logistic regression model. We recorded 458 carriage episodes including 318 (69.4 %) where the carried strain was shared with at least one other infant. The odds of strain sharing varied significantly across villages (χ2=47.5, df=21, P-value <0.001). Infants in close proximity to each other were more likely to be involved in strain sharing, but we also show a considerable amount of strain sharing across longer distances. Close geographic proximity (<5 km) between shared strains was associated with a significantly lower pairwise SNP distance compared to strains shared over longer distances (P-value <0.005). Sustained carriage of a shared strain among the infants was significantly more likely to occur if they resided in villages within a 5 km radius of each other (P-value <0.005, OR 3.7). Conversely, where both infants were transiently colonized by the shared strain, they were more likely to reside in villages separated by over 15 km (P-value <0.05, OR 1.5). PCV7 serotypes were rare (13.5 %) and were significantly less likely to be shared (P-value <0.001, OR -1.07). Strain sharing was more likely to occur over short geographical distances, especially where accompanied by sustained colonization. Our results show that strain sharing is a useful proxy for studying transmission dynamics in an under-sampled population with limited genomic data. This article contains data hosted by Microreact.

Despite contributing to the large disease burden in West Africa, little is known about the genomic epidemiology of Streptococcus pneumoniae which cause meningitis among children under 5 years old in the region. We analysed whole-genome sequencing data from 185 S. pneumoniae isolates recovered from suspected paediatric meningitis cases as part of the World Health Organization (WHO) invasive bacterial diseases surveillance from 2010 to 2016. The phylogeny was reconstructed, accessory genome similarity was computed and antimicrobial-resistance patterns were inferred from the genome data and compared to phenotypic resistance from disc diffusion. We studied the changes in the distribution of serotypes pre- and post-pneumococcal conjugate vaccine (PCV) introduction in the Central and Western sub-regions separately. The overall distribution of non-vaccine, PCV7 (4, 6B, 9V, 14, 18C, 19F and 23F) and additional PCV13 serotypes (1, 3, 5, 6A, 19A and 7F) did not change significantly before and after PCV introduction in the Central region (Fisher's test P value 0.27) despite an increase in the proportion of non-vaccine serotypes to 40 % (n=6) in the post-PCV introduction period compared to 21.9 % (n=14). In the Western sub-region, PCV13 serotypes were more dominant among isolates from The Gambia following the introduction of PCV7, 81 % (n=17), compared to the pre-PCV period in neighbouring Senegal, 51 % (n=27). The phylogeny illustrated the diversity of strains associated with paediatric meningitis in West Africa and highlighted the existence of phylogeographical clustering, with isolates from the same sub-region clustering and sharing similar accessory genome content. Antibiotic-resistance genotypes known to confer resistance to penicillin, chloramphenicol, co-trimoxazole and tetracycline were detected across all sub-regions. However, there was no discernible trend linking the presence of resistance genotypes with the vaccine introduction period or whether the strain was a vaccine or non-vaccine serotype. Resistance genotypes appeared to be conserved within selected sub-clades of the phylogenetic tree, suggesting clonal inheritance. Our data underscore the need for continued surveillance on the emergence of non-vaccine serotypes as well as chloramphenicol and penicillin resistance, as these antibiotics are likely still being used for empirical treatment in low-resource settings. This article contains data hosted by Microreact.