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 next generation technologies and bioinformatics pipelines for capsular typing of Streptococcus pneumoniae

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.

Genomics for public health and international surveillance of antimicrobial resistance

Historically, epidemiological investigation and surveillance for bacterial antimicrobial resistance (AMR) has relied on low-resolution isolate-based phenotypic analyses undertaken at local and national reference laboratories. Genomic sequencing has the potential to provide a far more high-resolution picture of AMR evolution and transmission, and is already beginning to revolutionise how public health surveillance networks monitor and tackle bacterial AMR. However, the routine integration of genomics in surveillance pipelines still has considerable barriers to overcome. In 2022, a workshop series and online consultation brought together international experts in AMR and pathogen genomics to assess the status of genomic applications for AMR surveillance in a range of settings. Here we focus on discussions around the use of genomics for public health and international AMR surveillance, noting the potential advantages of, and barriers to, implementation, and proposing recommendations from the working group to help to drive the adoption of genomics in public health AMR surveillance. These recommendations include the need to build capacity for genome sequencing and analysis, harmonising and standardising surveillance systems, developing equitable data sharing and governance frameworks, and strengthening interactions and relationships among stakeholders at multiple levels.

Genomic epidemiology of Streptococcus pneumoniae serotype 16F lineages

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.

Novel pneumococcal capsule type 33E results from the inactivation of glycosyltransferase WciE in vaccine type 33F

The polysaccharide (PS) capsule is essential for immune evasion and virulence of Streptococcus pneumoniae. Existing pneumococcal vaccines are designed to elicit anticapsule antibodies; however, the effectiveness of these vaccines is being challenged by the emergence of new capsule types or variants. Herein, we characterize a newly discovered capsule type, 33E, that appears to have repeatedly emerged from vaccine type 33F via an inactivation mutation in the capsule glycosyltransferase gene, wciE. Structural analysis demonstrated that 33E and 33F share an identical repeat unit backbone [→5)-β-D-Galf2Ac-(1→3)-β-D-Galp-(1→3)-α-D-Galp-(1→3)-β-D-Galf-(1→3)-β-D-Glcp-(1→], except that a galactose (α-D-Galp) branch is present in 33F but not in 33E. Though the two capsule types were indistinguishable using conventional typing methods, the monoclonal antibody Hyp33FM1 selectively bound 33F but not 33E pneumococci. Further, we confirmed that wciE encodes a glycosyltransferase that catalyzes the addition of the branching α-D-Galp and that its inactivation in 33F strains results in the expression of the 33E capsule type. Though 33F and 33E share a structural and antigenic similarity, our pilot study suggested that immunization with a 23-valent pneumococcal PS vaccine containing 33F PS did not significantly elicit cross-opsonic antibodies to 33E. New conjugate vaccines that target capsule type 33F may not necessarily protect against 33E. Therefore, studies of new conjugate vaccines require knowledge of the newly identified capsule type 33E and reliable pneumococcal typing methods capable of distinguishing it from 33F.

Molecular characterization of Streptococcus pneumoniae causing disease among children in Nigeria during the introduction of PCV10 (GSK)

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.

Evidence of virulence and antimicrobial resistance in Streptococcus pneumoniae serotype 16F lineages

Introduction

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

Aim

To characterise lineages and antimicrobial resistance in 16F isolates obtained from South Africa and placed the local findings in a global context.

Methodology

We analysed 10923 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 five continents, years covered (1995 - 2018), accessed on 17 th March 2022).

Results

Nine percent (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, however, 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]. 10% (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.

Discussion

Serotype 16F lineages are common in Southern Africa. Some of these lineages are associated with disease, and resistance to penicillin and cotrimoxazole. We recommend continuous genomic surveillance to determine 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.

DATA SUMMARY

The sequencing reads for the genomes analysed have been deposited in the European Nucleotide Archive and the accession numbers for each isolate are listed in Supplementary Table1 . Phylogenetic tree of serotype 16F pneumococcal genomes and associated metadata are available for download and visualisation on the Microreact website: Phylogenies of seotype 16F, GPSC33 and GPSC46 are available on the Microreact serotype-16F , GPSC33 and GPSC46 , respectively.

IMPACT STATEMENT

This study shows that serotype 16F lineages are predominant in Southern Africa and are associated with disease and antimicrobial resistance. Although serotype 16F has been included in the newer formulation of the upcoming vaccine formulations of PCV21 and IVT-25, continuous surveillance to determine long term impact of serotype 16F lineages on vaccines and antimicrobial therapy remains essential.

Discovery and Characterization of Pneumococcal Serogroup 36 Capsule Subtypes, Serotypes 36A and 36B

Streptococcus pneumoniae can produce a wide breadth of antigenically diverse capsule types, a fact that poses a looming threat to the success of vaccines that target pneumococcal polysaccharide (PS) capsule. Yet, many pneumococcal capsule types remain undiscovered and/or uncharacterized. Prior sequence analysis of pneumococcal capsule synthesis (cps) loci suggested the existence of capsule subtypes among isolates identified as "serotype 36" according to conventional capsule typing methods. We discovered these subtypes represent two antigenically similar but distinguishable pneumococcal capsule serotypes, 36A and 36B. Biochemical analysis of their capsule PS structure reveals that both have the shared repeat unit backbone [→5)-α-d-Galf-(1→1)-d-Rib-ol-(5→P→6)-β-d-ManpNAc-(1→4)-β-d-Glcp-(1→] with two branching structures. Both serotypes have a β-d-Galp branch to Ribitol. Serotypes 36A and 36B differ by the presence of a α-d-Glcp-(1→3)-β-d-ManpNAc or α-d-Galp-(1→3)-β-d-ManpNAc branch, respectively. Comparison of the phylogenetically distant serogroup 9 and 36 cps loci, which all encode this distinguishing glycosidic bond, revealed that the incorporation of Glcp (in types 9N and 36A) versus Galp (in types 9A, 9V, 9L, and 36B) is associated with the identity of four amino acids in the cps-encoded glycosyltransferase WcjA. Identifying functional determinants of cps-encoded enzymes and their impact on capsule PS structure is key to improving the resolution and reliability of sequencing-based capsule typing methods and discovering novel capsule variants indistinguishable by conventional serotyping methods.

A global genomic perspective on the multidrug-resistant Streptococcus pneumoniae 15A-CC63 sub-lineage following pneumococcal conjugate vaccine introduction

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.

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.

Emergence of a multidrug-resistant and virulent Streptococcus pneumoniae lineage mediates serotype replacement after PCV13: an international whole-genome sequencing study

BACKGROUND

Serotype 24F is one of the emerging pneumococcal serotypes after the introduction of pneumococcal conjugate vaccine (PCV). We aimed to identify lineages driving the increase of serotype 24F in France and place these findings into a global context.

METHODS

Whole-genome sequencing was performed on a collection of serotype 24F pneumococci from asymptomatic colonisation (n=229) and invasive disease (n=190) isolates among individuals younger than 18 years in France, from 2003 to 2018. To provide a global context, we included an additional collection of 24F isolates in the Global Pneumococcal Sequencing (GPS) project database for analysis. A Global Pneumococcal Sequence Cluster (GPSC) and a clonal complex (CC) were assigned to each genome. Phylogenetic, evolutionary, and spatiotemporal analysis were conducted using the same 24F collection and supplemented with a global collection of genomes belonging to the lineage of interest from the GPS project database (n=25 590).

FINDINGS

Serotype 24F was identified in numerous countries mainly due to the clonal spread of three lineages: GPSC10 (CC230), GPSC16 (CC156), and GPSC206 (CC7701). GPSC10 was the only multidrug-resistant lineage. GPSC10 drove the increase in 24F in France and had high invasive disease potential. The international dataset of GPSC10 (n=888) revealed that this lineage expressed 16 other serotypes, with only six included in 13-valent PCV (PCV13). All serotype 24F isolates were clustered in a single clade within the GPSC10 phylogeny and long-range transmissions were detected from Europe to other continents. Spatiotemporal analysis showed GPSC10-24F took 3-5 years to spread across France and a rapid change of serotype composition from PCV13 serotype 19A to 24F during the introduction of PCV13 was observed in neighbouring country Spain.

INTERPRETATION

Our work reveals that GPSC10 alone is a challenge for serotype-based vaccine strategy. More systematic investigation to identify lineages like GPSC10 will better inform and improve next-generation preventive strategies against pneumococcal diseases.

FUNDING

Bill & Melinda Gates Foundation, Wellcome Sanger Institute, and the US Centers for Disease Control and Prevention.

Key features of pneumococcal isolates recovered in Central and Northwestern Russia in 2011–2018 determined through whole-genome sequencing

Invasive pneumococcal disease remains one of the leading causes of morbidity and mortality worldwide. In Russia, 13- valent pneumococcal conjugate vaccine (PCV13) was introduced into the childhood immunization programme nationwide in 2014. As part of the Global Pneumococcal Sequencing Project (GPS), we used genome data to characterize 179 pneumococcal isolates collected from Russia in 2011–2018 to investigate the circulating pneumococcal strains using a standardized genomic definition of pneumococcal lineages (global pneumococcal sequence clusters, GPSCs), prevalent serotypes and antimicrobial resistance profiles.

We observed high serotype and lineage diversity among the 179 isolates recovered from cerebrospinal fluid (n=77), nasopharyngeal swabs (n=99) and other non-sterile site swabs (n=3). Overall, 60 GPSCs were identified, including 48 clonal complexes (CCs) and 14 singletons, and expressed 42 serotypes (including non-typable). Among PCV13 serotypes, 19F, 6B and 23F were the top three serotypes while 11A, 15B/C and 8 were the top three among non-PCV13 serotypes in the collection. Two lineages (GPSC6 and GPSC47) expressed both PCV13 and non-PCV13 serotypes that caused invasive disease, and were penicillin- and multidrug-resistant (MDR), highlighting their potential to adapt and continue to cause infections under vaccine and antibiotic selective pressure. PCV13 serotypes comprised 92 % (11/12) of the CSF isolates from the children aged below 5 years; however, the prevalence of PCV13 serotype isolates dropped to 53 % (31/58) among the nasopharyngeal isolates. Our analysis showed that 59 % (105/179) of the isolates were predicted to be non-susceptible to at least one class of antibiotics and 26 % (46/179) were MDR. Four MDR lineages (GPSC1, GPSC6, GPSC10 and GPSC47) accounted for 65 % (30/46) of the MDR isolates and expressed PCV13 serotypes (93 %, 28/30).

This study provides evidence of high genetic and serotype diversity contributed by a mix of globally spreading and regionally circulating lineages in Russia. The observations suggest that the PCV13 vaccine could be important in reducing both invasive disease and antimicrobial resistance. We also identify potential lineages (GPSC6 and GPSC47) that may evade the vaccine.

Variants of Streptococcus pneumoniae Serotype 14 from Papua New Guinea with the Potential to Be Mistyped and Escape Vaccine-Induced Protection

Streptococcus pneumoniae (the pneumococcus) is a human pathogen of global importance, classified into serotypes based on the type of capsular polysaccharide produced. Serotyping of pneumococci is essential for disease surveillance and vaccine impact measurement. However, the accuracy of serotyping methods can be affected by previously undiscovered variants. Previous studies have identified variants of serotype 14, a highly invasive serotype included in all licensed vaccine formulations. However, the potential of these variants to influence serotyping accuracy and evade vaccine-induced protection has not been investigated. In this study, we screened 1,386 nasopharyngeal swabs from children hospitalized with acute respiratory infection in Papua New Guinea for pneumococci. Swabs containing pneumococci (n = 1,226) were serotyped by microarray to identify pneumococci with a divergent serotype 14 capsule locus. Three serotype 14 variants ('14-like') were isolated and characterized further. The serotyping results of these isolates using molecular methods varied depending on the method, with 3/3 typing as nontypeable (PneumoCaT), 3/3 typing as serotype 14 (seroBA), and 2/3 typing as serotype 14 (SeroCall and quantitative PCR). All three isolates were nontypeable by phenotypic methods (Quellung and latex agglutination), indicating the absence of capsule. Illumina and nanopore sequencing were employed to examine their capsule loci and revealed unique mutations. Lastly, when incubated with sera from vaccinated individuals, the 14-like isolates evaded serotype-specific opsonophagocytic killing. Our study highlights the need for phenotypic testing to validate serotyping data derived from molecular methods. The convergent evolution of capsule loss underscores the importance of studying pneumococcal population biology to monitor the emergence of pneumococci capable of vaccine escape, globally. IMPORTANCE Pneumococcus is a pathogen of major public health importance. Current vaccines have limited valency, targeting a subset (up to 20) of the more than 100 capsule types (serotypes). Precise serotyping methods are therefore essential to avoid mistyping, which can reduce the accuracy of data used to inform decisions around vaccine introduction and/or maintenance of national vaccination programs. In this study, we examine a variant of serotype 14 (14-like), a virulent serotype present in all currently licensed vaccine formulations. Although these 14-like pneumococci no longer produce a serotype 14 capsule, widely used molecular methods can mistype them as serotype 14. Importantly, we show that 14-like pneumococci can evade opsonophagocytic killing mediated by vaccination. Despite the high accuracy of molecular methods for serotyping, our study reemphasizes their limitations. This is particularly relevant in situations where nonvaccine type pneumococci (e.g., the 14-likes in this study) could potentially be misidentified as a vaccine type (e.g., serotype 14).

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.

Comparative Genomics of Disease and Carriage Serotype 1 Pneumococci

The isolation of Streptococcus pneumoniae serotypes in systemic tissues of patients with invasive disease versus the nasopharynx of healthy individuals with asymptomatic carriage varies widely. Some serotypes are hyper-invasive, particularly serotype 1, but the underlying genetics remain poorly understood due to the rarity of carriage isolates, reducing the power of comparison with invasive isolates. Here, we use a well-controlled genome-wide association study to search for genetic variation associated with invasiveness of serotype 1 pneumococci from a serotype 1 endemic setting in Africa. We found no consensus evidence that certain genomic variation is overrepresented among isolates from patients with invasive disease than asymptomatic carriage. Overall, the genomic variation explained negligible phenotypic variability, suggesting a minimal effect on the disease status. Furthermore, changes in lineage distribution were seen with lineages replacing each other over time, highlighting the importance of continued pathogen surveillance. Our findings suggest that the hyper-invasiveness is an intrinsic property of the serotype 1 strains, not specific for a "disease-associated" subpopulation disproportionately harboring unique genomic variation.

Widespread sharing of pneumococcal strains in a rural African setting: proximate villages are more likely to share similar strains that are carried at multiple timepoints

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 (χ²=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.

Endocarditis caused by non-typeable Streptococcus pneumoniae

The Streptococcus pneumoniae capsule is regarded as indispensable in bacteremia. We report an infant with a ventricular septal defect and infective endocarditis caused by nontypeable S. pneumoniae. In-depth investigation confirmed a deficient capsule yet favored pneumococcal fitness for causing infective endocarditis, rather than a host immune disorder, as the cause of infective endocarditis in this case.

RCandy: an R package for visualizing homologous recombinations in bacterial genomes

SUMMARY

Homologous recombination is an important evolutionary process in bacteria and other prokaryotes, which increases genomic sequence diversity and can facilitate adaptation. Several methods and tools have been developed to detect genomic regions recently affected by recombination. Exploration and visualization of such recombination events can reveal valuable biological insights, but it remains challenging. Here, we present RCandy, a platform-independent R package for rapid, simple and flexible visualization of recombination events in bacterial genomes.

AVAILABILITY AND IMPLEMENTATION

RCandy is an R package freely available for use under the MIT license. It is platform-independent and has been tested on Windows, Linux and MacOSX. The source code comes together with a detailed vignette available on GitHub at https://github.com/ChrispinChaguza/RCandy.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Population genetic structure, serotype distribution and antibiotic resistance of Streptococcus pneumoniae causing invasive disease in children in Argentina

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.

Streptococcus pneumoniae genomic datasets from an Indian population describing pre-vaccine evolutionary epidemiology using a whole genome sequencing approach

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⁻¹) 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.

Phylogeography and resistome of pneumococcal meningitis in West Africa before and after vaccine introduction

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.

Using genomics to examine the persistence of Streptococcus pneumoniae serotype 19A in Ireland and the emergence of a sub-clade associated with vaccine failures

BACKGROUND

Streptococcus pneumoniae serotype 19A remains a significant cause of invasive pneumococcal disease (IPD) in Ireland despite the successful introduction of a 13-valent pneumococcal conjugate vaccine (PCV13) in 2010 which reduced the overall incidence of IPD in children.

METHODS

Invasive Streptococcus pneumoniae serotype 19A isolates from the Irish reference laboratory between 2007-08 and 2017-18 were analysed using whole genome sequencing (WGS) to investigate the persistence of this vaccine-preventable serotype. We compared the entire national 19A collection to other international collections using a standardised nomenclature of Global Pneumococcal Sequencing Clusters (GPSC).

RESULTS

Expansion of GPSCs and clonal complexes (CCs) may have been associated with vaccine introduction and antimicrobial prescribing policies. A sub-clade of GPSC1-CC320 (n = 25) unique to Ireland, included five of the ten vaccine failures/breakthrough cases identified (p = 0.0086). This sub-clade was not observed in a global GPSC1-CC320 collection. All isolates within the sub-clade (n = 25) contained a galE gene variant rarely observed in a global pneumococcal collection (n = 37/13454, p < 0.001) nor within GPSC1-CC320 (n = 19/227) (p < 0.001). The sub-clade was estimated to have emerged at the start of the PCV-vaccine era (ancestral origin 2000, range 1995-2004) and expanded in Ireland, with most isolated after PCV13 introduction (n = 24/25).

CONCLUSIONS

The identification of a sub-clade/variant of serotype 19A highlights the benefit of using WGS to analyse genotypes associated with persistence of a preventable serotype of S. pneumoniae. Particularly as this sub-clade identified was more likely to be associated with IPD in vaccinated children than other 19A genotypes. It is possible that changes to the galE gene, which is involved in capsule production but outside of the capsular polysaccharide biosynthesis locus, may affect bacterial persistence within the population. Discrete changes associated with vaccine-serotype persistence should be further investigated and may inform vaccine strategies.

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.

A mosaic tetracycline resistance gene tet(S/M) detected in an MDR pneumococcal CC230 lineage that underwent capsular switching in South Africa

OBJECTIVES

We reported tet(S/M) in Streptococcus pneumoniae and investigated its temporal spread in relation to nationwide clinical interventions.

METHODS

We whole-genome sequenced 12 254 pneumococcal isolates from 29 countries on an Illumina HiSeq sequencer. Serotype, multilocus ST and antibiotic resistance were inferred from genomes. An SNP tree was built using Gubbins. Temporal spread was reconstructed using a birth-death model.

RESULTS

We identified tet(S/M) in 131 pneumococcal isolates and none carried other known tet genes. Tetracycline susceptibility testing results were available for 121 tet(S/M)-positive isolates and all were resistant. A majority (74%) of tet(S/M)-positive isolates were from South Africa and caused invasive diseases among young children (59% HIV positive, where HIV status was available). All but two tet(S/M)-positive isolates belonged to clonal complex (CC) 230. A global phylogeny of CC230 (n=389) revealed that tet(S/M)-positive isolates formed a sublineage predicted to exhibit resistance to penicillin, co-trimoxazole, erythromycin and tetracycline. The birth-death model detected an unrecognized outbreak of this sublineage in South Africa between 2000 and 2004 with expected secondary infections (effective reproductive number, R) of ∼2.5. R declined to ∼1.0 in 2005 and <1.0 in 2012. The declining epidemic could be related to improved access to ART in 2004 and introduction of pneumococcal conjugate vaccine (PCV) in 2009. Capsular switching from vaccine serotype 14 to non-vaccine serotype 23A was observed within the sublineage.

CONCLUSIONS

The prevalence of tet(S/M) in pneumococci was low and its dissemination was due to an unrecognized outbreak of CC230 in South Africa. Capsular switching in this MDR sublineage highlighted its potential to continue to cause disease in the post-PCV13 era.

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.

A Nationwide Outbreak of Invasive Pneumococcal Disease in Israel Caused by Streptococcus Pneumoniae Serotype 2

BACKGROUND

Invasive pneumococcal disease (IPD) caused by Streptococcus pneumoniae serotype 2 (Sp2) is infrequent. Large scale outbreaks have not been reported following pneumococcal conjugate vaccine (PCV) implementation. We describe a Sp2 IPD outbreak in Israel, in the 13-valent PCV (PCV13) era, with focus on Sp2 population structure and evolutionary dynamics.

METHODS

The data derived from a population-based, nationwide active surveillance of IPD since 2009. 7-valent PCV (PCV7)/PCV13 vaccines were introduced in July 2009 and November 2010, respectively. Sp2 isolates were tested for antimicrobial susceptibility, Multilocus Sequence Typing (MLST) and Whole Genome Sequencing (WGS) analysis.

RESULTS

Overall, 170 Sp2 IPD cases were identified during 2009-2019; Sp2 increased in 2015 and caused 6% of IPD during 2015-2019, a 7-fold increase compared with 2009-2014.The outbreak was caused by a previously unreported molecular type (ST-13578), initially observed in Israel in 2014. This clone caused 88% of Sp2 during 2015-2019. ST-13578 is a single-locus variant of ST-1504, previously reported globally, including in Israel. WGS analysis confirmed clonality among the ST-13578 population. Single-nucleotide polymorphisms-dense regions support a hypothesis that the ST-13578 outbreak clone evolved from ST-1504 by recombination.All tested strains were penicillin-susceptible (MIC <0.06 μg/mL). The ST-13578 clone was identified almost exclusively (99%) in the Jewish population and was mainly distributed in 3/7 Israeli districts. The outbreak is still ongoing, although declining since 2017.Conclusions: To the best of our knowledge, this is the first widespread Sp2 outbreak since PCV13 introduction worldwide, caused by the emerging ST-13578 clone.

Characterization of Pneumococcal Colonization Dynamics and Antimicrobial Resistance Using Shotgun Metagenomic Sequencing in Intensively Sampled South African Infants

Background: There remains a significant proportion of deaths due to pneumococcal pneumonia in infants from low- and middle-income countries despite the marginal global declines recorded in the past decade. Monitoring changes in pneumococcal carriage is key to understanding vaccination-induced shifts in the ecology of carriage, patterns of antimicrobial resistance, and impact on health. We longitudinally investigated pneumococcal carriage dynamics in PCV-13 vaccinated infants by collecting nasopharyngeal (NP) samples at 2-weekly intervals from birth through the first year of life from 137 infants. As a proof of concept, 196 NP samples were retrieved from a subset of 23 infants to explore strain-level pneumococcal colonization patterns and associated antimicrobial-resistance determinants. These were selected on the basis of changes in serotype and antibiogram over time. NP samples underwent short-term enrichment for streptococci prior to total nucleic acid extraction and whole metagenome shotgun sequencing (WMGS). Reads were assembled and aligned to pneumococcal reference genomes for the extraction of pneumococcal and non-pneumococcal bacterial reads. Pneumococcal contigs were aligned to the Antibiotic Resistance Gene-ANNOTation database of acquired AMR genes. In silico pneumococcal capsular and multilocus sequence typing were performed.

Results: Of the 196 samples sequenced, 174 had corresponding positive cultures for pneumococci, of which, 152 were assigned an in silico serotype. Metagenomic sequencing detected a single pneumococcal serotype in 85% (129/152), and co-colonization in 15% (23/152) of the samples. Twenty-two different pneumococcal serotypes were identified, with 15B/15C and 16F being the most common non-PCV13 serotypes, while 23F and 19A were the most common PCV13 serotypes. Twenty-six different sequence types (STs), including four novel STs were identified in silico. Mutations in the folA and folP genes, associated with cotrimoxazole resistance, were detected in 89% (87/98) of cotrimoxazole-non-susceptible pneumococci, as well as in the pbp1a and pbp2x genes, in penicillin non-susceptible ST7052^15B/15C isolates.

Conclusions: Metagenomic sequencing of NP samples is a valuable culture-independent technique for a detailed evaluation of the pneumococcal component and resistome of the NP microbiome. This method allowed for the detection of novel STs, as well as co-colonization, with a predominance of non-PCV13 serotypes in this cohort. Forty-eight resistance genes, as well as mutations associated with resistance were detected, but the correlation with phenotypic non-susceptibility was lower than expected.

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.

Visualizing variation within Global Pneumococcal Sequence Clusters (GPSCs) and country population snapshots to contextualize pneumococcal isolates

Knowledge of pneumococcal lineages, their geographic distribution and antibiotic resistance patterns, can give insights into global pneumococcal disease. We provide interactive bioinformatic outputs to explore such topics, aiming to increase dissemination of genomic insights to the wider community, without the need for specialist training. We prepared 12 country-specific phylogenetic snapshots, and international phylogenetic snapshots of 73 common Global Pneumococcal Sequence Clusters (GPSCs) previously defined using PopPUNK, and present them in Microreact. Gene presence and absence defined using Roary, and recombination profiles derived from Gubbins are presented in Phandango for each GPSC. Temporal phylogenetic signal was assessed for each GPSC using BactDating. We provide examples of how such resources can be used. In our example use of a country-specific phylogenetic snapshot we determined that serotype 14 was observed in nine unrelated genetic backgrounds in South Africa. The international phylogenetic snapshot of GPSC9, in which most serotype 14 isolates from South Africa were observed, highlights that there were three independent sub-clusters represented by South African serotype 14 isolates. We estimated from the GPSC9-dated tree that the sub-clusters were each established in South Africa during the 1980s. We show how recombination plots allowed the identification of a 20 kb recombination spanning the capsular polysaccharide locus within GPSC97. This was consistent with a switch from serotype 6A to 19A estimated to have occured in the 1990s from the GPSC97-dated tree. Plots of gene presence/absence of resistance genes (tet, erm, cat) across the GPSC23 phylogeny were consistent with acquisition of a composite transposon. We estimated from the GPSC23-dated tree that the acquisition occurred between 1953 and 1975. Finally, we demonstrate the assignment of GPSC31 to 17 externally generated pneumococcal serotype 1 assemblies from Utah via Pathogenwatch. Most of the Utah isolates clustered within GPSC31 in a USA-specific clade with the most recent common ancestor estimated between 1958 and 1981. The resources we have provided can be used to explore to data, test hypothesis and generate new hypotheses. The accessible assignment of GPSCs allows others to contextualize their own collections beyond the data presented here.

Carriage Dynamics of Pneumococcal Serotypes in Naturally Colonized Infants in a Rural African Setting During the First Year of Life

Streptococcus pneumoniae (the pneumococcus) carriage precedes invasive disease and influences population-wide strain dynamics, but limited data exist on temporal carriage patterns of serotypes due to the prohibitive costs of longitudinal studies. Here, we report carriage prevalence, clearance and acquisition rates of pneumococcal serotypes sampled from newborn infants bi-weekly from weeks 1 to 27, and then bi-monthly from weeks 35 to 52 in the Gambia. We used sweep latex agglutination and whole genome sequencing to serotype the isolates. We show rapid pneumococcal acquisition with nearly 31% of the infants colonized by the end of first week after birth and quickly exceeding 95% after 2 months. Co-colonization with multiple serotypes was consistently observed in over 40% of the infants at each sampling point during the first year of life. Overall, the mean acquisition time and carriage duration regardless of serotype was 38 and 24 days, respectively, but varied considerably between serotypes comparable to observations from other regions. Our data will inform disease prevention and control measures including providing baseline data for parameterising infectious disease mathematical models including those assessing the impact of clinical interventions such as pneumococcal conjugate vaccines.

Pneumococcal lineages associated with serotype replacement and antibiotic resistance in childhood invasive pneumococcal disease in the post-PCV13 era: an international whole-genome sequencing study

BACKGROUND

Invasive pneumococcal disease remains an important health priority owing to increasing disease incidence caused by pneumococci expressing non-vaccine serotypes. We previously defined 621 Global Pneumococcal Sequence Clusters (GPSCs) by analysing 20 027 pneumococcal isolates collected worldwide and from previously published genomic data. In this study, we aimed to investigate the pneumococcal lineages behind the predominant serotypes, the mechanism of serotype replacement in disease, as well as the major pneumococcal lineages contributing to invasive pneumococcal disease in the post-vaccine era and their antibiotic resistant traits.

METHODS

We whole-genome sequenced 3233 invasive pneumococcal disease isolates from laboratory-based surveillance programmes in Hong Kong (n=78), Israel (n=701), Malawi (n=226), South Africa (n=1351), The Gambia (n=203), and the USA (n=674). The genomes represented pneumococci from before and after pneumococcal conjugate vaccine (PCV) introductions and were from children younger than 3 years. We identified predominant serotypes by prevalence and their major contributing lineages in each country, and assessed any serotype replacement by comparing the incidence rate between the pre-PCV and PCV periods for Israel, South Africa, and the USA. We defined the status of a lineage as vaccine-type GPSC (≥50% 13-valent PCV [PCV13] serotypes) or non-vaccine-type GPSC (>50% non-PCV13 serotypes) on the basis of its initial serotype composition detected in the earliest vaccine period to measure their individual contribution toward serotype replacement in each country. Major pneumococcal lineages in the PCV period were identified by pooled incidence rate using a random effects model.

FINDINGS

The five most prevalent serotypes in the PCV13 period varied between countries, with only serotypes 5, 12F, 15B/C, 19A, 33F, and 35B/D common to two or more countries. The five most prevalent serotypes in the PCV13 period varied between countries, with only serotypes 5, 12F, 15B/C, 19A, 33F, and 35B/D common to two or more countries. These serotypes were associated with more than one lineage, except for serotype 5 (GPSC8). Serotype replacement was mainly mediated by expansion of non-vaccine serotypes within vaccine-type GPSCs and, to a lesser extent, by increases in non-vaccine-type GPSCs. A globally spreading lineage, GPSC3, expressing invasive serotypes 8 in South Africa and 33F in the USA and Israel, was the most common lineage causing non-vaccine serotype invasive pneumococcal disease in the PCV13 period. We observed that same prevalent non-vaccine serotypes could be associated with distinctive lineages in different countries, which exhibited dissimilar antibiotic resistance profiles. In non-vaccine serotype isolates, we detected significant increases in the prevalence of resistance to penicillin (52 [21%] of 249 vs 169 [29%] of 575, p=0·0016) and erythromycin (three [1%] of 249 vs 65 [11%] of 575, p=0·0031) in the PCV13 period compared with the pre-PCV period.

INTERPRETATION

Globally spreading lineages expressing invasive serotypes have an important role in serotype replacement, and emerging non-vaccine serotypes associated with different pneumococcal lineages in different countries might be explained by local antibiotic-selective pressures. Continued genomic surveillance of the dynamics of the pneumococcal population with increased geographical representation in the post-vaccine period will generate further knowledge for optimising future vaccine design.

FUNDING

Bill & Melinda Gates Foundation, Wellcome Sanger Institute, and the US Centers for Disease Control.

Putative novel cps loci in a large global collection of pneumococci

The pneumococcus produces a polysaccharide capsule, encoded by the cps locus, that provides protection against phagocytosis and determines serotype. Nearly 100 serotypes have been identified with new serotypes still being discovered, especially in previously understudied regions. Here we present an analysis of the cps loci of more than 18  000 genomes from the Global Pneumococcal Sequencing (GPS) project with the aim of identifying novel cps loci with the potential to produce previously unrecognized capsule structures. Serotypes were assigned using whole genome sequence data and 66 of the approximately 100 known serotypes were included in the final dataset. Closer examination of each serotype’s sequences identified nine putative novel cps loci (9X, 11X, 16X, 18X1, 18X2, 18X3, 29X, 33X and 36X) found in ~2.6  % of the genomes. The large number and global distribution of GPS genomes provided an unprecedented opportunity to identify novel cps loci and consider their phylogenetic and geographical distribution. Nine putative novel cps loci were identified and examples of each will undergo subsequent structural and immunological analysis.

International genomic definition of pneumococcal lineages, to contextualise disease, antibiotic resistance and vaccine impact

BACKGROUND

Pneumococcal conjugate vaccines have reduced the incidence of invasive pneumococcal disease, caused by vaccine serotypes, but non-vaccine-serotypes remain a concern. We used whole genome sequencing to study pneumococcal serotype, antibiotic resistance and invasiveness, in the context of genetic background.

METHODS

Our dataset of 13,454 genomes, combined with four published genomic datasets, represented Africa (40%), Asia (25%), Europe (19%), North America (12%), and South America (5%). These 20,027 pneumococcal genomes were clustered into lineages using PopPUNK, and named Global Pneumococcal Sequence Clusters (GPSCs). From our dataset, we additionally derived serotype and sequence type, and predicted antibiotic sensitivity. We then measured invasiveness using odds ratios that relating prevalence in invasive pneumococcal disease to carriage.

FINDINGS

The combined collections (n = 20,027) were clustered into 621 GPSCs. Thirty-five GPSCs observed in our dataset were represented by >100 isolates, and subsequently classed as dominant-GPSCs. In 22/35 (63%) of dominant-GPSCs both non-vaccine serotypes and vaccine serotypes were observed in the years up until, and including, the first year of pneumococcal conjugate vaccine introduction. Penicillin and multidrug resistance were higher (p < .05) in a subset dominant-GPSCs (14/35, 9/35 respectively), and resistance to an increasing number of antibiotic classes was associated with increased recombination (R2 = 0.27 p < .0001). In 28/35 dominant-GPSCs, the country of isolation was a significant predictor (p < .05) of its antibiogram (mean misclassification error 0.28, SD ± 0.13). We detected increased invasiveness of six genetic backgrounds, when compared to other genetic backgrounds expressing the same serotype. Up to 1.6-fold changes in invasiveness odds ratio were observed.

INTERPRETATION

We define GPSCs that can be assigned to any pneumococcal genomic dataset, to aid international comparisons. Existing non-vaccine-serotypes in most GPSCs preclude the removal of these lineages by pneumococcal conjugate vaccines; leaving potential for serotype replacement. A subset of GPSCs have increased resistance, and/or serotype-independent invasiveness.

Fast and flexible bacterial genomic epidemiology with PopPUNK

The routine use of genomics for disease surveillance provides the opportunity for high-resolution bacterial epidemiology. Current whole-genome clustering and multilocus typing approaches do not fully exploit core and accessory genomic variation, and they cannot both automatically identify, and subsequently expand, clusters of significantly similar isolates in large data sets spanning entire species. Here, we describe PopPUNK (Population Partitioning Using Nucleotide K-mers), a software implementing scalable and expandable annotation- and alignment-free methods for population analysis and clustering. Variable-length k-mer comparisons are used to distinguish isolates’ divergence in shared sequence and gene content, which we demonstrate to be accurate over multiple orders of magnitude using data from both simulations and genomic collections representing 10 taxonomically widespread species. Connections between closely related isolates of the same strain are robustly identified, despite interspecies variation in the pairwise distance distributions that reflects species’ diverse evolutionary patterns. PopPUNK can process 10³–10⁴ genomes in a single batch, with minimal memory use and runtimes up to 200-fold faster than existing model-based methods. Clusters of strains remain consistent as new batches of genomes are added, which is achieved without needing to reanalyze all genomes de novo. This facilitates real-time surveillance with consistent cluster naming between studies and allows for outbreak detection using hundreds of genomes in minutes. Interactive visualization and online publication is streamlined through the automatic output of results to multiple platforms. PopPUNK has been designed as a flexible platform that addresses important issues with currently used whole-genome clustering and typing methods, and has potential uses across bacterial genetics and public health research.

Electricity production from municipal solid waste using microbial fuel cells

The organic content of municipal solid waste has long been an attractive source of renewable energy, mainly as a solid fuel in waste-to-energy plants. This study focuses on the potential to use microbial fuel cells to convert municipal solid waste organics into energy using various operational conditions. The results showed that two-chamber microbial fuel cells with carbon felt and carbon felt allocation had a higher maximal power density (20.12 and 30.47 mW m(-2) for 1.5 and 4 L, respectively) than those of other electrode plate allocations. Most two-chamber microbial fuel cells (1.5 and 4 L) had a higher maximal power density than single-chamber ones with corresponding electrode plate allocations. Municipal solid waste with alkali hydrolysis pre-treatment and K3Fe(CN)6 as an electron acceptor improved the maximal power density to 1817.88 mW m(-2) (~0.49% coulomb efficiency, from 0.05-0.49%). The maximal power density from experiments using individual 1.5 and 4 L two-chamber microbial fuel cells, and serial and parallel connections of 1.5 and 4 L two-chamber microbial fuel cells, was found to be in the order of individual 4 L (30.47 mW m(-2)) > serial connection of 1.5 and 4 L (27.75) > individual 1.5 L (20.12) > parallel connection of 1.5 and 4 L (17.04) two-chamber microbial fuel cells . The power density using municipal solid waste microbial fuel cells was compared with information in the literature and discussed.

Effects of different SRT on anaerobic digestion of MSW dosed with various MSWI ashes

This study investigated different solid retention time (SRT) on municipal solid waste (MSW) anaerobic digestion with various MSW incinerator fly ash (FA) and bottom ash (BA) addition. Results showed that biogas production rates (BPRs, ≈ 200 to ≈ 400 mL/gVS) with organic loading rate of ≈ 0.053 gVS/gVS(reactor) (Day 1-435, SRT 20 days, SRT20) at FA 1g/d (FA1), BA 12 g/d (BA12) and BA 24 g/d (BA24) dosed bioreactors increased after adaptation. BPRs with SRT10 and SRT5 decreased while BPRs with SRT40 showed to increase compared to initial BPRs (≈ 200 mL/gVS) with SRT20. SRT5 operation reduced the BPRs (≈ 10 - ≈ 90 mL/gVS) significantly and only BA12 and BA24 dosed bioreactors could recover the BPRs (≈ 100 - ≈ 200 mL/gVS) after SRT20 operation (Day 613-617) compared to FA1 and FA3 and control. Released levels of Co, Mo and W at BA12 and BA24 dosed bioreactors showed most potential to improve MSW anaerobic digestion.

Effects of micro-nano and non micro-nano MSWI ashes addition on MSW anaerobic digestion

This study aims at investigating the effects of micro-nano municipal solid waste (MSW) incinerator (MSWI) fly ash (FA) and bottom ash (BA) on the MSW anaerobic digestion. Results showed that suitable micro-nano and non micro-nano MSWI ashes addition (FA/MSW 3, 6, 18 and 30 g g(-1) VS and BA/MSW 12, 36, 60 and 120 g g(-1) VS) could enhance the biogas production compared to the control. It was particularly found to have the highest biogas production at the micro-nano MSWI BA/MSW ratio of 36 g g(-1) VS (∼193 mL g(-1) VS MSW, ∼3.5 times to the control). Micro-nano MSWI FA and BA added bioreactors had the higher biogas production than the corresponding non micro-nano MSWI FA and BA added ones. Suitable MSWI ashes addition could improve the biogas production due to the released metals levels suitable for the MSW anaerobic digestion particularly found in the micro-nano added bioreactors.

Genetic identity of aminoglycoside-resistance genes in Escherichia coli isolates from human and animal sources

A bacterial collection (n=249) obtained in Hong Kong from 2002 to 2004 was used to investigate the molecular epidemiology of aminoglycoside resistance among Escherichia coli isolates from humans and food-producing animals. Of these, 89 isolates were gentamicin-sensitive (human n=60, animal n=29) and 160 isolates were gentamicin-resistant (human n=107, animal n=53). Overall, 84.1% (90/107) and 75.5% (40/53) of the gentamicin-resistant isolates from human and animal sources, respectively, were found to possess the aacC2 gene. The aacC2 gene for 20 isolates (10 each for human and animal isolates) was sequenced. Two alleles were found that were equally distributed in human and animal isolates. PFGE showed that the gentamicin-resistant isolates exhibited diverse patterns with little clonality. In some isolates, the aacC2 gene was encoded on large transferable plasmids of multiple incompatibility groups (IncF, IncI1 and IncN). An IncFII plasmid of 140 kb in size was shared by one human and three animal isolates. In summary, this study showed that human and animal isolates share the same pool of resistance genes.

Presynaptic mechanisms underlying cannabinoid inhibition of excitatory synaptic transmission in rat striatal neurons

The striatum is a crucial site of action for the motor effects of cannabinoids (CBs). However, the electrophysiological consequences of activation of CB receptors on the striatal neurons have not been established. Here we report for the first time that the cannabimimetic aminoalkylindole WIN 55,212-2 and the endogenous cannabinoid anandamide substantially depress corticostriatal glutamatergic synaptic transmission onto striatal neurons in the brain slice preparation. The selective CB1 receptor antagonist SR 141716 effectively reversed this inhibition. WIN 55,212-2 significantly increased the paired-pulse facilitation of synaptically evoked EPSCs, while having no effect on the sensitivity of postsynaptic neurons to [alpha]-amino-3-hydroxy-5-methylisoxazole-4-propionic acid. WIN 55,212-2 also reduced the frequency of spontaneous, action potential-dependent EPSCs (sEPSCs) without altering their amplitude distribution. Superfusion of WIN 55,212-2 elicited a membrane hyperpolarization accompanied by a decrease in input resistance. Both effects were blocked by intracellular caesium. In contrast, intracellular caesium failed to affect WIN 55,212-2-mediated synaptic inhibition. The WIN 55,212-2-mediated synaptic inhibition was blocked by the Gi/o protein inhibitor pertussis toxin (PTX), but not by the GABA(A) receptor antagonist bicuculline or GABA(B) receptor antagonist SCH 50911. Pretreatment with the N-type Ca2+ channel antagonist [omega]-conotoxin GVIA selectively abolished the WIN-55,212-2-mediated synaptic inhibition. These results suggest that cannabinoids depress the corticostriatal glutamatergic synaptic transmission through the activation of presynaptic CB1 receptors to inhibit N-type Ca2+ channel activity, which in turn reduces glutamate release. The presynaptic action of cannabinoids is mediated by a PTX-sensitive Gi/o protein-coupled signalling pathway.

Cortical renal leiomyoma with extension to renal pelvis

We describe a case of renal leiomyoma in a 21-year-old woman who presented with flank pain and hematuria. Urographic and computed tomographic (CT) studies revealed a large right renal mass with polypoid outgrowth protruding into the renal pelvis. Cortical renal leiomyoma with this radiographic manifestation is extremely rare.