The Contribution of Genetic Variation of Streptococcus pneumoniae to the Clinical Manifestation of Invasive Pneumococcal Disease

Serotype epidemiology and case-fatality risk of invasive pneumococcal disease: a nationwide population study from Switzerland, 2012-2022

In Switzerland, thirteen-valent pneumococcal conjugate vaccine (PCV13) has been introduced in 2011. During the COVID-19 pandemic, cases of invasive pneumococcal disease (IPD) have decreased but consequences on the serotype epidemiology are less clear. The objective of the study has been to analyse the impact of PCV13 introduction and the COVID-19 pandemic on the IPD epidemiology and investigate the changes in the case fatality risk (CFR). We analysed data from the Swiss nationwide surveillance for the period 2012-2022. Poisson and logistic regression analyses were performed allowing us to inspect trends over time and to define serotypes that are associated with case fatality. In total, 8747 IPD cases were included from 2012 to 2022. IPD incidence dropped in the years 2020 (6.0/100,000) and 2021 (5.5/100,000) but recovered in 2022 (9.1/100,000). While the incidence numbers of patients >65 years did not reach the pre-pandemic level, numbers significantly increased in infants <1 year in 2022 (IRR 1.08, 95%CI: 1.01-1.16). The incidence of PCV13 serotypes among all IPD cases decreased until 2019 before increasing again during the pandemic (in 2022). Logistic regression analyses revealed that the PCV20 serotype 11A (OR: 1.76, 95%CI: 1.14-2.64), and the PCV13 serotypes 3 (OR: 1.26, 95% CI: 1.04-1.53) and 19F (OR: 1.76, 95%CI: 1.14-2.65) were significantly associated with increased CFR. In conclusion, the COVID-19 pandemic has had only minor temporary effects on the serotype distribution. Continued use of vaccines with extended serotype coverage may further reduce IPD disease burden and mortality.

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

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

Pneumococcal genetic variability in age-dependent bacterial carriage

The characteristics of pneumococcal carriage vary between infants and adults. Host immune factors have been shown to contribute to these age-specific differences, but the role of pathogen sequence variation is currently less well-known. Identification of age-associated pathogen genetic factors could leadto improved vaccine formulations. We therefore performed genome sequencing in a large carriage cohort of children and adults and combined this with data from an existing age-stratified carriage study. We compiled a dictionary of pathogen genetic variation, including serotype, strain, sequence elements, single-nucleotide polymorphisms (SNPs), and clusters of orthologous genes (COGs) for each cohort - all of which were used in a genome-wide association with host age. Age-dependent colonization showed weak evidence of being heritable in the first cohort (h2 = 0.10, 95% CI 0.00-0.69) and stronger evidence in the second cohort (h2 = 0.56, 95% CI 0.23-0.87). We found that serotypes and genetic background (strain) explained a proportion of the heritability in the first cohort (h2serotype = 0.07, 95% CI 0.04-0.14 and h2GPSC = 0.06, 95% CI 0.03-0.13) and the second cohort (h2serotype = 0.11, 95% CI 0.05-0.21 and h2GPSC = 0.20, 95% CI 0.12-0.31). In a meta-analysis of these cohorts, we found one candidate association (p=1.2 × 10-9) upstream of an accessory Sec-dependent serine-rich glycoprotein adhesin. Overall, while we did find a small effect of pathogen genome variation on pneumococcal carriage between child and adult hosts, this was variable between populations and does not appear to be caused by strong effects of individual genes. This supports proposals for adaptive future vaccination strategies that are primarily targeted at dominant circulating serotypes and tailored to the composition of the pathogen populations.

PowerBacGWAS: a computational pipeline to perform power calculations for bacterial genome-wide association studies

Genome-wide association studies (GWAS) are increasingly being applied to investigate the genetic basis of bacterial traits. However, approaches to perform power calculations for bacterial GWAS are limited. Here we implemented two alternative approaches to conduct power calculations using existing collections of bacterial genomes. First, a sub-sampling approach was undertaken to reduce the allele frequency and effect size of a known and detectable genotype-phenotype relationship by modifying phenotype labels. Second, a phenotype-simulation approach was conducted to simulate phenotypes from existing genetic variants. We implemented both approaches into a computational pipeline (PowerBacGWAS) that supports power calculations for burden testing, pan-genome and variant GWAS; and applied it to collections of Enterococcus faecium, Klebsiella pneumoniae and Mycobacterium tuberculosis. We used this pipeline to determine sample sizes required to detect causal variants of different minor allele frequencies (MAF), effect sizes and phenotype heritability, and studied the effect of homoplasy and population diversity on the power to detect causal variants. Our pipeline and user documentation are made available and can be applied to other bacterial populations. PowerBacGWAS can be used to determine sample sizes required to find statistically significant associations, or the associations detectable with a given sample size. We recommend to perform power calculations using existing genomes of the bacterial species and population of study.

Streptococcus pneumoniae: a Plethora of Temperate Bacteriophages With a Role in Host Genome Rearrangement

Bacteriophages (phages) are viruses that infect bacteria. They are the most abundant biological entity on Earth (current estimates suggest there to be perhaps 10³¹ particles) and are found nearly everywhere. Temperate phages can integrate into the chromosome of their host, and prophages have been found in abundance in sequenced bacterial genomes. Prophages may modulate the virulence of their host in different ways, e.g., by the secretion of phage-encoded toxins or by mediating bacterial infectivity. Some 70% of Streptococcus pneumoniae (the pneumococcus)—a frequent cause of otitis media, pneumonia, bacteremia and meningitis—isolates harbor one or more prophages. In the present study, over 4000 S. pneumoniae genomes were examined for the presence of prophages, and nearly 90% were found to contain at least one prophage, either defective (47%) or present in full (43%). More than 7000 complete putative integrases, either of the tyrosine (6243) or serine (957) families, and 1210 full-sized endolysins (among them 1180 enzymes corresponding to 318 amino acid-long N-acetylmuramoyl-L-alanine amidases [LytA_PPH]) were found. Based on their integration site, 26 different pneumococcal prophage groups were documented. Prophages coding for tRNAs, putative virulence factors and different methyltransferases were also detected. The members of one group of diverse prophages (PPH090) were found to integrate into the 3’ end of the host lytA_Spn gene encoding the major S. pneumoniae autolysin without disrupting it. The great similarity of the lytA_Spnand lytA_PPH genes (85–92% identity) allowed them to recombine, via an apparent integrase-independent mechanism, to produce different DNA rearrangements within the pneumococcal chromosome. This study provides a complete dataset that can be used to further analyze pneumococcal prophages, their evolutionary relationships, and their role in the pathogenesis of pneumococcal disease.

Examining the Distribution and Impact of Single-Nucleotide Polymorphisms in the Capsular Locus of Streptococcus pneumoniae Serotype 19A

Streptococcus pneumoniae serotype 19A prevalence has increased after the implementation of the PCV7 and PCV10 vaccines. In this study, we have provided, with high accuracy, the genetic diversity of the 19A serotype in a cohort of Dutch invasive pneumococcal disease patients and asymptomatic carriers obtained in the period from 2004 to 2016. The whole genomes of the 338 pneumococcal isolates in this cohort were sequenced and their capsule (cps) loci compared to examine their diversity and determine the impact on the production of capsular polysaccharide (CPS) sugar precursors and CPS shedding. We discovered 79 types with a unique cps locus sequence. Most variation was observed in the rmlB and rmlD genes of the TDP-Rha synthesis pathway and in the wzg gene, which is of unknown function. Interestingly, gene variation in the cps locus was conserved in multiple alleles. Using RmlB and RmlD protein models, we predict that enzymatic function is not affected by the single-nucleotide polymorphisms as identified. To determine if RmlB and RmlD function was affected, we analyzed nucleotide sugar levels using ultrahigh-performance liquid chromatography-mass spectrometry (UHPLC-MS). CPS precursors differed between 19A cps locus subtypes, including TDP-Rha, but no clear correlation was observed. Also, significant differences in multiple nucleotide sugar levels were observed between phylogenetically branched groups. Because of indications of a role for Wzg in capsule shedding, we analyzed if this was affected. No clear indication of a direct role in shedding was found. We thus describe genotypic variety in rmlB, rmlD, and wzg in serotype 19A in the Netherlands, for which we have not discovered an associated phenotype.

The drop in reported invasive pneumococcal disease among adults during the first COVID-19 wave in the Netherlands explained

Objectives: Streptococcus pneumoniae is the leading bacterial pathogen causing respiratory infections. Since the COVID-19 pandemic emerged, less invasive pneumococcal disease (IPD) was identified by surveillance systems worldwide. Measures to prevent transmission of SARS-CoV-2 also reduce transmission of pneumococci, but this would gradually lead to lower disease rates. Design: Here, we explore additional factors contributing to the instant drop in pneumococcal disease cases captured in surveillance. Results: Our observations on referral practices and other impediments to diagnostic testing indicate that residual IPD has likely occurred but remained undetected by conventional hospital-based surveillance. Conclusions: Depending on the setting, we discuss alternative monitoring strategies that could improve understanding of pneumococcal disease dynamics.

Forest and Trees: Exploring Bacterial Virulence with Genome-wide Association Studies and Machine Learning

The advent of inexpensive and rapid sequencing technologies has allowed bacterial whole-genome sequences to be generated at an unprecedented pace. This wealth of information has revealed an unanticipated degree of strain-to-strain genetic diversity within many bacterial species. Awareness of this genetic heterogeneity has corresponded with a greater appreciation of intraspecies variation in virulence. A number of comparative genomic strategies have been developed to link these genotypic and pathogenic differences with the aim of discovering novel virulence factors. Here, we review recent advances in comparative genomic approaches to identify bacterial virulence determinants, with a focus on genome-wide association studies and machine learning.

A Narrative Review of the Molecular Epidemiology and Laboratory Surveillance of Vaccine Preventable Bacterial Meningitis Agents: Streptococcus pneumoniae, Neisseria meningitidis, Haemophilus influenzae and Streptococcus agalactiae

This narrative review describes the public health importance of four most common bacterial meningitis agents, Streptococcus pneumoniae, Neisseria meningitidis, Haemophilus influenzae, and S. agalactiae (group B Streptococcus). Three of them are strict human pathogens that normally colonize the nasopharynx and may invade the blood stream to cause systemic infections and meningitis. S. agalactiae colonizes the genito-gastrointestinal tract and is an important meningitis agent in newborns, but also causes invasive infections in infants or adults. These four bacteria have polysaccharide capsules that protect them against the host complement defense. Currently licensed conjugate vaccines (against S. pneumoniae, H. influenza, and N. meningitidis only but not S. agalactiae) can induce protective serum antibodies in infants as young as two months old offering protection to the most vulnerable groups, and the ability to eliminate carriage of homologous serotype strains in vaccinated subjects lending further protection to those not vaccinated through herd immunity. However, the serotype-specific nature of these vaccines have driven the bacteria to adapt by mechanisms that affect the capsule antigens through either capsule switching or capsule replacement in addition to the possibility of unmasking of strains or serotypes not covered by the vaccines. The post-vaccine molecular epidemiology of vaccine-preventable bacterial meningitis is discussed based on findings obtained with newer genomic laboratory surveillance methods.

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.

Different Types of Coagulase Are Associated With 28-Day Mortality in Patients With Staphylococcus aureus Bloodstream Infections

Background:Staphylococcus aureus (S. aureus), a leading cause of bacteremia and infective endocarditis, exploits the human coagulation system by using a wide range of specific virulence factors. However, the impact of these host-pathogen interactions on the outcome of patients with Staphylococcus aureus bacteremia (SAB) remains unclear.

Methods: A total of 178 patients with S. aureus bacteremia were included and analyzed regarding bacterial factors (coa gene size, vWbp, clfA, clfB, fnbA, fnbB, fib) and clinical parameters. A stepwise multivariate Cox regression model and a Partitioning Around Medoids (PAM) cluster algorithm were used for statistical analysis.

Results: Patients' risk factors for 28-day mortality were creatinine (OR 1.49, p < 0.001), age (OR 1.9, p < 0.002), fibrinogen (OR 0.44, p < 0.004), albumin (OR 0.63, p < 0.02), hemoglobin (OR 0.59, p < 0.03), and CRP (OR 1.72, p < 0.04). Five distinct bacterial clusters with different mortality rates were unveiled, whereof two showed a 2-fold increased mortality and an accumulation of specific coagulase gene sizes, 547-base pairs and 660-base pairs.

Conclusions: Based on the data obtained in the present study an association of coagulase gene size and fib regarding 28-day mortality was observed in patients with S. aureus bloodstream infections. Further animal and prospective clinical studies are needed to confirm our preliminary findings.

Genome-wide association studies of Shigella spp. and Enteroinvasive Escherichia coli isolates demonstrate an absence of genetic markers for prediction of disease severity

BACKGROUND

We investigated the association of symptoms and disease severity of shigellosis patients with genetic determinants of infecting Shigella and entero-invasive Escherichia coli (EIEC), because determinants that predict disease outcome per individual patient could be used to prioritize control measures. For this purpose, genome wide association studies (GWAS) were performed using presence or absence of single genes, combinations of genes, and k-mers. All genetic variants were derived from draft genome sequences of isolates from a multicenter cross-sectional study conducted in the Netherlands during 2016 and 2017. Clinical data of patients consisting of binary/dichotomous representation of symptoms and their calculated severity scores were also available from this study. To verify the suitability of the methods used, the genetic differences between the genera Shigella and Escherichia were used as control.

RESULTS

The isolates obtained were representative of the population structure encountered in other Western European countries. No association was found between single genes or combinations of genes and separate symptoms or disease severity scores. Our benchmark characteristic, genus, resulted in eight associated genes and > 3,000,000 k-mers, indicating adequate performance of the algorithms used.

CONCLUSIONS

To conclude, using several microbial GWAS methods, genetic variants in Shigella spp. and EIEC that can predict specific symptoms or a more severe course of disease were not identified, suggesting that disease severity of shigellosis is dependent on other factors than the genetic variation of the infecting bacteria. Specific genes or gene fragments of isolates from patients are unsuitable to predict outcomes and cannot be used for development, prioritization and optimization of guidelines for control measures of shigellosis or infections with EIEC.

Joint sequencing of human and pathogen genomes reveals the genetics of pneumococcal meningitis

Streptococcus pneumoniae is a common nasopharyngeal colonizer, but can also cause life-threatening invasive diseases such as empyema, bacteremia and meningitis. Genetic variation of host and pathogen is known to play a role in invasive pneumococcal disease, though to what extent is unknown. In a genome-wide association study of human and pathogen we show that human variation explains almost half of variation in susceptibility to pneumococcal meningitis and one-third of variation in severity, identifying variants in CCDC33 associated with susceptibility. Pneumococcal genetic variation explains a large amount of invasive potential (70%), but has no effect on severity. Serotype alone is insufficient to explain invasiveness, suggesting other pneumococcal factors are involved in progression to invasive disease. We identify pneumococcal genes involved in invasiveness including pspC and zmpD, and perform a human-bacteria interaction analysis. These genes are potential candidates for the development of more broadly-acting pneumococcal vaccines.

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.

Blood‒Brain Barrier Pathology and CNS Outcomes in Streptococcus pneumoniae Meningitis

Streptococcus pneumoniae is a major meningitis-causing pathogen globally, bringing about significant morbidity and mortality, as well as long-term neurological sequelae in almost half of the survivors. Subsequent to nasopharyngeal colonisation and systemic invasion, translocation across the blood‒brain barrier (BBB) by S. pneumoniae is a crucial early step in the pathogenesis of meningitis. The BBB, which normally protects the central nervous system (CNS) from deleterious molecules within the circulation, becomes dysfunctional in S. pneumoniae invasion due to the effects of pneumococcal toxins and a heightened host inflammatory environment of cytokines, chemokines and reactive oxygen species intracranially. The bacteria‒host interplay within the CNS likely determines not only the degree of BBB pathological changes, but also host survival and the extent of neurological damage. This review explores the relationship between S. pneumoniae bacteria and the host inflammatory response, with an emphasis on the BBB and its roles in CNS protection, as well as both the acute and long-term pathogenesis of meningitis.

A molecular warning system for invasive pneumococcus

A genome-wide association study of Streptococcus pneumoniae identifies and validates bacterial loci associated with invasive infection.