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

BACKGROUND: Pneumococcal conjugate vaccines (PCVs) have reduced pneumococcal diseases globally. Pneumococcal genomic surveys elucidate PCV effects on population structure but are rarely conducted in low-income settings despite the high disease burden. METHODS: We undertook whole-genome sequencing (WGS) of 660 pneumococcal isolates collected through surveys from healthy carriers 2 years from 13-valent PCV (PCV13) introduction and 1 year after rollout in northern Malawi. We investigated changes in population structure, within-lineage serotype dynamics, serotype diversity, and frequency of antibiotic resistance (ABR) and accessory genes. RESULTS: In children <5 years of age, frequency and diversity of vaccine serotypes (VTs) decreased significantly post-PCV, but no significant changes occurred in persons ≥5 years of age. Clearance of VT serotypes was consistent across different genetic backgrounds (lineages). There was an increase of nonvaccine serotypes (NVTs)-namely 7C, 15B/C, and 23A-in children <5 years of age, but 28F increased in both age groups. While carriage rates have been recently shown to remain stable post-PCV due to replacement serotypes, there was no change in diversity of NVTs. Additionally, frequency of intermediate-penicillin-resistant lineages decreased post-PCV. Although frequency of ABR genes remained stable, other accessory genes, especially those associated with mobile genetic element and bacteriocins, showed changes in frequency post-PCV. CONCLUSIONS: We demonstrate evidence of significant population restructuring post-PCV driven by decreasing frequency of vaccine serotypes and increasing frequency of few NVTs mainly in children under 5. Continued surveillance with WGS remains crucial to fully understand dynamics of the residual VTs and replacement NVT serotypes post-PCV.

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 Distribution of Invasive Serotype 35D Streptococcus pneumoniae Isolates following Introduction of 13-Valent Pneumococcal Conjugate Vaccine.

A newly recognized pneumococcal serotype, 35D, which differs from the 35B polysaccharide in structure and serology by not binding to factor serum 35a, was recently reported. The genetic basis for this distinctive serology is due to the presence of an inactivating mutation in wciG, which encodes an O-acetyltransferase responsible for O-acetylation of a galactofuranose. Here, we assessed the genomic data of a worldwide pneumococcal collection to identify serotype 35D isolates and understand their geographical distribution, genetic background, and invasiveness potential. Of 21,980 pneumococcal isolates, 444 were originally typed as serotype 35B by PneumoCaT. Analysis of the wciG gene revealed 23 isolates from carriage (n = 4) and disease (n = 19) with partial or complete loss-of-function mutations, including mutations resulting in premature stop codons (n = 22) and an in-frame mutation (n = 1). These were selected for further analysis. The putative 35D isolates were geographically widespread, and 65.2% (15/23) of them was recovered after the introduction of pneumococcal conjugate vaccine 13 (PCV13). Compared with serotype 35B isolates, putative serotype 35D isolates have higher invasive disease potentials based on odds ratios (OR) (11.58; 95% confidence interval[CI], 1.42 to 94.19 versus 0.61; 95% CI, 0.40 to 0.92) and a higher prevalence of macrolide resistance mediated by mefA (26.1% versus 7.6%; P = 0.009). Using the Quellung reaction, 50% (10/20) of viable isolates were identified as serotype 35D, 25% (5/20) as serotype 35B, and 25% (5/20) as a mixture of 35B/35D. The discrepancy between phenotype and genotype requires further investigation. These findings illustrated a global distribution of an invasive serotype, 35D, among young children post-PCV13 introduction and underlined the invasive potential conferred by the loss of O-acetylation in the pneumococcal capsule.

Epidemiological and Molecular Characterization of an Invasive Group A Streptococcus emm32.2 Outbreak.

An emm32.2 invasive group A streptococcus (iGAS) outbreak occurred in Liverpool from January 2010 to September 2012. This genotype had not previously been identified in Liverpool, but was responsible for 32% (14/44) of all iGAS cases reported during this time period. We performed a case-case comparison of emm32.2 iGAS cases with non-emm32.2 control iGAS cases identified in the Liverpool population over the same time period to assess patient risk factors for emm32.2 iGAS infection. The emm32.2 iGAS cases were confined to the adult population. We show that homelessness, intravenous drug use, and alcohol abuse predisposed patients to emm32.2 iGAS disease; however, no obvious epidemiological linkage between the patients with emm32.2 iGAS could be identified. Comparative whole-genome sequencing analysis of emm32.2 iGAS and non-emm32.2 control isolates was also performed to identify pathogen factors which might have driven the outbreak. We identified 19 genes, five of which had previously been implicated in virulence, which were present in all of the emm32.2 iGAS isolates but not present in any of the non-emm32.2 control isolates. We report that a novel emm32.2 genotype emerged in Liverpool in 2010 and identified a specific subset of genes, which could have allowed this novel emm32.2 genotype to persist in a disadvantaged population in the region over a 3-year period.

Phylogenetic Analysis of Invasive Serotype 1 Pneumococcus in South Africa, 1989 to 2013.

Serotype 1 is an important cause of invasive pneumococcal disease in South Africa and has declined following the introduction of the 13-valent pneumococcal conjugate vaccine in 2011. We genetically characterized 912 invasive serotype 1 isolates from 1989 to 2013. Simpson's diversity index (D) and recombination ratios were calculated. Factors associated with sequence types (STs) were assessed. Clonal complex 217 represented 96% (872/912) of the sampled isolates. Following the introduction of the 13-valent pneumococcal conjugate vaccine (PCV13), ST diversity increased in children <5 years (D, 0.39 to 0.63, P = 0.002) and individuals >14 years (D, 0.35 to 0.54, P < 0.001): ST-217 declined proportionately in children <5 years (153/203 [75%] versus 21/37 [57%], P = 0.027) and individuals >14 years (242/305 [79%] versus 96/148 [65%], P = 0.001), whereas ST-9067 increased (4/684 [0.6%] versus 24/228 [11%], P < 0.001). Three subclades were identified within ST-217: ST-217C1 (353/382 [92%]), ST-217C2 (15/382 [4%]), and ST-217C3 (14/382 [4%]). ST-217C2, ST-217C3, and single-locus variant (SLV) ST-8314 (20/912 [2%]) were associated with nonsusceptibility to chloramphenicol, tetracycline, and co-trimoxazole. ST-8314 (20/912 [2%]) was also associated with increased nonsusceptibility to penicillin (P < 0.001). ST-217C3 and newly reported ST-9067 had higher recombination ratios than those of ST-217C1 (4.344 versus 0.091, P < 0.001; and 0.086 versus 0.013, P < 0.001, respectively). Increases in genetic diversity were noted post-PCV13, and lineages associated with antimicrobial nonsusceptibility were identified.

Understanding pneumococcal serotype 1 biology through population genomic analysis.

BACKGROUND: Pneumococcus kills over one million children annually and over 90 % of these deaths occur in low-income countries especially in Sub-Saharan Africa (SSA) where HIV exacerbates the disease burden. In SSA, serotype 1 pneumococci particularly the endemic ST217 clone, causes majority of the pneumococcal disease burden. To understand the evolution of the virulent ST217 clone, we analysed ST217 whole genomes from isolates sampled from African and Asian countries. METHODS: We analysed 226 whole genome sequences from the ST217 lineage sampled from 9 African and 4 Asian countries. We constructed a whole genome alignment and used it for phylogenetic and coalescent analyses. We also screened the genomes to determine presence of antibiotic resistance conferring genes. RESULTS: Population structure analysis grouped the ST217 isolates into five sequence clusters (SCs), which were highly associated with different geographical regions and showed limited intracontinental and intercontinental spread. The SCs showed lower than expected genomic sequence, which suggested strong purifying selection and small population sizes caused by bottlenecks. Recombination rates varied between the SCs but were lower than in other successful clones such as PMEN1. African isolates showed higher prevalence of antibiotic resistance genes than Asian isolates. Interestingly, certain West African isolates harbored a defective chloramphenicol and tetracycline resistance-conferring element (Tn5253) with a deletion in the loci encoding the chloramphenicol resistance gene (cat pC194), which caused lower chloramphenicol than tetracycline resistance. Furthermore, certain genes that promote colonisation were absent in the isolates, which may contribute to serotype 1's rarity in carriage and consequently its lower recombination rates. CONCLUSIONS: The high phylogeographic diversity of the ST217 clone shows that this clone has been in circulation globally for a long time, which allowed its diversification and adaptation in different geographical regions. Such geographic adaptation reflects local variations in selection pressures in different locales. Further studies will be required to fully understand the biological mechanisms which makes the ST217 clone highly invasive but unable to successfully colonise the human nasopharynx for long durations which results in lower recombination rates.

Comparative Genomic Analysis of Meningitis- and Bacteremia-Causing Pneumococci Identifies a Common Core Genome.

Streptococcus pneumoniae is a nasopharyngeal commensal that occasionally invades normally sterile sites to cause bloodstream infection and meningitis. Although the pneumococcal population structure and evolutionary genetics are well defined, it is not clear whether pneumococci that cause meningitis are genetically distinct from those that do not. Here, we used whole-genome sequencing of 140 isolates of S. pneumoniae recovered from bloodstream infection (n = 70) and meningitis (n = 70) to compare their genetic contents. By fitting a double-exponential decaying-function model, we show that these isolates share a core of 1,427 genes (95% confidence interval [CI], 1,425 to 1,435 genes) and that there is no difference in the core genome or accessory gene content from these disease manifestations. Gene presence/absence alone therefore does not explain the virulence behavior of pneumococci that reach the meninges. Our analysis, however, supports the requirement of a range of previously described virulence factors and vaccine candidates for both meningitis- and bacteremia-causing pneumococci. This high-resolution view suggests that, despite considerable competency for genetic exchange, all pneumococci are under considerable pressure to retain key components advantageous for colonization and transmission and that these components are essential for access to and survival in sterile sites.

High multiple carriage and emergence of Streptococcus pneumoniae vaccine serotype variants in Malawian children.

BACKGROUND: Carriage of either single or multiple pneumococcal serotypes (multiple carriage) is a prerequisite for developing invasive pneumococcal disease. However, despite the reported high rates of pneumococcal carriage in Malawi, no data on carriage of multiple serotypes has been reported previously. Our study provides the first description of the prevalence of multiple pneumococcal carriage in Malawi. METHODS: The study was conducted in Blantyre and Karonga districts in Malawi, from 2008 to 2012. We recruited 116 children aged 0-13 years. These children were either HIV-infected (N = 44) or uninfected (N = 72). Nasopharyngeal samples were collected using sterile swabs. Pneumococcal serotypes in the samples were identified by microarray. Strains that could not be typed by microarray were sequenced to characterise possible genetic alterations within the capsular polysaccharide (CPS) locus. RESULTS: The microarray identified 179 pneumococcal strains (from 116 subjects), encompassing 43 distinct serotypes and non-typeable (NT) strains. Forty per cent (46/116) of children carried multiple serotypes. Carriage of vaccine type (VT) strains was higher (p = 0.028) in younger (0-2 years) children (71 %, 40/56) compared to older (3-13 years) children (50 %, 30/60). Genetic variations within the CPS locus of known serotypes were observed in 19 % (34/179) of the strains identified. The variants included 13-valent pneumococcal conjugate vaccine (PCV13) serotypes 6B and 19A, and the polysaccharide vaccine serotype 20. Serotype 6B variants were the most frequently isolated (47 %, 16/34). Unlike the wild type, the CPS locus of the 6B variants contained an insertion of the licD-family phosphotransferase gene. The CPS locus of 19A- and 20-variants contained an inversion in the sugar-biosynthesis (rmlD) gene and a 717 bp deletion within the transferase (whaF) gene, respectively. CONCLUSIONS: The high multiple carriage in Malawian children provides opportunities for genetic exchange through horizontal gene transfer. This may potentially lead to CPS locus variants and vaccine escape. Variants reported here occurred naturally, however, PCV13 introduction could exacerbate the CPS genetic variations. Further studies are therefore recommended to assess the invasive potential of these variants and establish whether PCV13 would offer cross-protection. We have shown that younger children (0-2 years) are a reservoir of VT serotypes, which makes them an ideal target for vaccination.

Genomic identification of a novel co-trimoxazole resistance genotype and its prevalence amongst Streptococcus pneumoniae in Malawi.

OBJECTIVES: This study aimed to define the molecular basis of co-trimoxazole resistance in Malawian pneumococci under the dual selective pressure of widespread co-trimoxazole and sulfadoxine/pyrimethamine use. METHODS: We measured the trimethoprim and sulfamethoxazole MICs and analysed folA and folP nucleotide and translated amino acid sequences for 143 pneumococci isolated from carriage and invasive disease in Malawi (2002-08). RESULTS: Pneumococci were highly resistant to both trimethoprim and sulfamethoxazole (96%, 137/143). Sulfamethoxazole-resistant isolates showed a 3 or 6 bp insertion in the sulphonamide-binding site of folP. The trimethoprim-resistant isolates fell into three genotypic groups based on dihydrofolate reductase (encoded by folA) mutations: Ile-100-Leu (10%), the Ile-100-Leu substitution together with a residue 92 substitution (56%) and those with a novel uncharacterized resistance genotype (34%). The nucleotide sequence divergence and dN/dS of folA and folP remained stable from 2004 onwards. CONCLUSIONS: S. pneumoniae exhibit almost universal co-trimoxazole resistance in vitro and in silico that we believe is driven by extensive co-trimoxazole and sulfadoxine/pyrimethamine use. More than one-third of pneumococci employ a novel mechanism of co-trimoxazole resistance. Resistance has now reached a point of stabilizing evolution. The use of co-trimoxazole to prevent pneumococcal infection in HIV/AIDS patients in sub-Saharan Africa should be re-evaluated.

Shigella Detection and Molecular Serotyping With a Customized TaqMan Array Card in the Enterics for Global Health (EFGH): Shigella Surveillance Study.

Background: Quantitative polymerase chain reaction (qPCR) targeting ipaH has been proven to be highly efficient in detecting Shigella in clinical samples compared to culture-based methods, which underestimate Shigella burden by 2- to 3-fold. qPCR assays have also been developed for Shigella speciation and serotyping, which is critical for both vaccine development and evaluation. Methods: The Enterics for Global Health (EFGH) Shigella surveillance study will utilize a customized real-time PCR-based TaqMan Array Card (TAC) interrogating 82 targets, for the detection and differentiation of Shigella spp, Shigella sonnei, Shigella flexneri serotypes, other diarrhea-associated enteropathogens, and antimicrobial resistance (AMR) genes. Total nucleic acid will be extracted from rectal swabs or stool samples, and assayed on TAC. Quantitative analysis will be performed to determine the likely attribution of Shigella and other particular etiologies of diarrhea using the quantification cycle cutoffs derived from previous studies. The qPCR results will be compared to conventional culture, serotyping, and phenotypic susceptibility approaches in EFGH. Conclusions: TAC enables simultaneous detection of diarrheal etiologies, the principal pathogen subtypes, and AMR genes. The high sensitivity of the assay enables more accurate estimation of Shigella-attributed disease burden, which is critical to informing policy and in the design of future clinical trials.

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.

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.

Invasive Streptococcus pneumoniae in children, Malawi, 2004-2006.

Of 176 invasive Streptococcus pneumoniae isolates from children in Malawi, common serotypes were 1 (23%), 6A/B (18%), 14 (6%), and 23F (6%). Coverage with the 7-valent pneumococcal conjugate vaccine (PCV) was 39%; PCV10 and PCV13 increased coverage to 66% and 88%, respectively. We found chloramphenicol resistance in 27% of isolates and penicillin nonsusceptibility in 10% (by using meningitis breakpoints); all were ceftriaxone susceptible.

Changes in Invasive Pneumococcal Disease Caused by Streptococcus pneumoniae Serotype 1 Following Introduction of PCV10 and PCV13: Findings from the PSERENADE Project.

Streptococcus pneumoniae serotype 1 (ST1) was an important cause of invasive pneumococcal disease (IPD) globally before the introduction of pneumococcal conjugate vaccines (PCVs) containing ST1 antigen. The Pneumococcal Serotype Replacement and Distribution Estimation (PSERENADE) project gathered ST1 IPD surveillance data from sites globally and aimed to estimate PCV10/13 impact on ST1 IPD incidence. We estimated ST1 IPD incidence rate ratios (IRRs) comparing the pre-PCV10/13 period to each post-PCV10/13 year by site using a Bayesian multi-level, mixed-effects Poisson regression and all-site IRRs using a linear mixed-effects regression (N = 45 sites). Following PCV10/13 introduction, the incidence rate (IR) of ST1 IPD declined among all ages. After six years of PCV10/13 use, the all-site IRR was 0.05 (95% credibility interval 0.04-0.06) for all ages, 0.05 (0.04-0.05) for <5 years of age, 0.08 (0.06-0.09) for 5-17 years, 0.06 (0.05-0.08) for 18-49 years, 0.06 (0.05-0.07) for 50-64 years, and 0.05 (0.04-0.06) for ≥65 years. PCV10/13 use in infant immunization programs was followed by a 95% reduction in ST1 IPD in all ages after approximately 6 years. Limited data availability from the highest ST1 disease burden countries using a 3+0 schedule constrains generalizability and data from these settings are needed.

Lower Density and Shorter Duration of Nasopharyngeal Carriage by Pneumococcal Serotype 1 (ST217) May Explain Its Increased Invasiveness over Other Serotypes.

Streptococcus pneumoniae is a frequent colonizer of the human nasopharynx and a major cause of life-threating invasive infections such as pneumonia, meningitis and sepsis. Over 1 million people die every year due to invasive pneumococcal disease (IPD), mainly in developing countries. Serotype 1 is a common cause of IPD; however, unlike other serotypes, it is rarely found in the carrier state in the nasopharynx, which is often considered a prerequisite for disease. The aim of this study was to understand this dichotomy. We used murine models of carriage and IPD to characterize the pathogenesis of African serotype 1 (sequence type 217) pneumococcal strains obtained from the Queen Elizabeth Central Hospital in Blantyre, Malawi. We found that ST217 pneumococcal strains were highly virulent in a mouse model of invasive pneumonia, but in contrast to the generally accepted assumption, can also successfully establish nasopharyngeal carriage. Interestingly, we found that cocolonizing serotypes may proliferate in the presence of serotype 1, suggesting that acquisition of serotype 1 carriage could increase the risk of developing IPD by other serotypes. RNA sequencing analysis confirmed that key virulence genes associated with inflammation and tissue invasiveness were upregulated in serotype 1. These data reveal important new insights into serotype 1 pathogenesis, with implications for carriage potential and risk of invasive disease through interactions with other cocolonizing serotypes, an often-overlooked factor in transmission and disease progression.IMPORTANCE The pneumococcus causes serious diseases such as pneumonia, sepsis, and meningitis and is a major cause of morbidity and mortality worldwide. Serotype 1 accounts for the majority of invasive pneumococcal disease cases in sub-Saharan Africa but is rarely found during nasopharyngeal carriage. Understanding the mechanisms leading to nasopharyngeal carriage and invasive disease by this serotype can help reduce its burden on health care systems worldwide. In this study, we also uncovered the potential impact of serotype 1 on disease progression of other coinfecting serotypes, which can have important implications for vaccine efficacy. Understanding the interactions between different serotypes during nasopharyngeal carriage may lead to improved intervention methods and therapies to reduce pneumococcal invasive disease levels.

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.

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.

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.