Molecular Characterization Based on Whole-Genome Sequencing of Streptococcus pneumoniae in Children Living in Southwest China During 2017-2019

Limited protection of pneumococcal vaccines against emergent Streptococcus pneumoniae serotype 14/ST876 strains

PURPOSE

Streptococcus pneumoniae (Spn) is a major cause of child death. We investigated the epidemiology of S. pneumoniae in a pediatric fever clinic and explored the genomics basis of the limited vaccine response of serotype 14 strains worldwide.

METHODS

Febrile disease and pneumonia were diagnosed following criteria from the WHO at the end of 2019 at a tertiary children's hospital. Spn was isolated by culture from nasopharyngeal (NP) swabs. The density was determined by lytA-base qPCR. Isolates were serotyped by Quellung and underwent antimicrobial susceptibility testing. Whole-genome sequencing was employed for molecular serotyping, MLST, antibiotic gene determination, SNP calling, recombination prediction, and phylogenetic analysis.

RESULTS

The presence of pneumococcus in the nasopharynx (87.5%, 7/8, p = 0.0227) and a high carriage (100%, 7/7, p = 0.0123) were significantly associated with pneumonia development. Living with siblings (73.7%, 14/19, p = 0.0125) and non-vaccination (56.0%, 28/50, p = 0.0377) contributed significantly to the Spn carriage. Serotype 14 was the most prevalent strain (16.67%, 5/30). The genome analysis of 1497 serotype 14 strains indicated S14/ST876 strains were only prevalent in China, presented limited vaccine responses with higher recombination activities within its cps locus, and unique variation patterns in the genes wzg and lrp.

CONCLUSION

With the lifting of the one-child policy, it will be crucial for families with multiple children to get PCV vaccinations in China. Due to the highly variant cps locus and distinctive variation patterns in capsule shedding and binding proteins genes, the prevalent S14/ST876 strains have shown poor response to current vaccines. It is necessary to continue monitoring the molecular epidemiology of this vaccine escape clone.

Transcriptional regulation of TacL-mediated lipoteichoic acids biosynthesis by ComE during competence impacts pneumococcal transformation

Competence development is essential for bacterial transformation since it enables bacteria to take up free DNA from the surrounding environment. The regulation of teichoic acid biosynthesis is tightly controlled during pneumococcal competence; however, the mechanism governing this regulation and its impact on transformation remains poorly understood. We demonstrated that a defect in lipoteichoic acid ligase (TacL)-mediated lipoteichoic acids (LTAs) biosynthesis was associated with impaired pneumococcal transformation. Using a fragment of tacL regulatory probe as bait in a DNA pulldown assay, we successfully identified several regulatory proteins, including ComE. Electrophoretic mobility shift assays revealed that phosphomimetic ComE, but not wild-type ComE, exhibited specific binding to the probe. DNase I footprinting assays revealed the specific binding sequences encompassing around 30 base pairs located 31 base pairs upstream from the start codon of tacL. Expression of tacL was found to be upregulated in the ΔcomE strain, and the addition of exogenous competence-stimulating peptide repressed the tacL transcription in the wild-type strain but not the ΔcomE mutant, indicating that ComE exerted a negative regulatory effect on the transcription of tacL. Mutation in the JH2 region of tacL upstream regulatory sequence led to increased LTAs abundance and displayed higher transformation efficiency. Collectively, our work identified the regulatory mechanisms that control LTAs biosynthesis during competence and thereby unveiled a repression mechanism underlying pneumococcal transformation.

Whole genome sequencing data of Streptococcus pneumoniae isolated from Indonesian population

Streptococcus pneumoniae is the leading cause of bacterial pneumonia, bacteremia, and meningitis. Indonesia introduced the pneumococcal conjugate vaccine (PCV) nationwide in 2022. In this study, we present whole genome sequence (WGS) data of 94 S. pneumoniae isolates that were obtained from hospitalized patients, healthy children, and adult groups from different regions prior to PCV program in Indonesia. DNA sequences of S. pneumoniae were obtained using the TruSeq Nano DNA kit (Illumina NovaSeq6000 Platform). The genome data of S. pneumoniae features a 1,969,562 bp to 2,741,371 bp circular chromosome with 39-40% G+C content. The genome includes 1935-3319 coding sequences (CDS), 2 to 5 rRNA genes, 43 to 49 tRNA genes, and 56 to 71 ncRNA. These data will be useful for analyzing the serotype, sequence type, virulence genes, antimicrobial resistance genes, and the impact of pneumococcal vaccination in Indonesia. The FASTQ raw files of these sequences are available under BioProject accession number PRJNA995903 and Sequence Read Archive accession numbers SRR25316461-SRR25316554.

Serotype, antibiotic susceptibility and whole-genome characterization of Streptococcus pneumoniae in all age groups living in Southwest China during 2018-2022

Background

Streptococcus pneumoniae is a common pathogen that colonizes the human upper respiratory tract, causing high morbidity and mortality worldwide. This study aimed to investigate the prevalence status of S. pneumoniae isolated from patients of all ages in Southwest China, including serotype, antibiotic susceptibility and other molecular characteristics, to provide a basis for clinical antibiotic usage and vaccine development.

Methods

This study was conducted from January 2018 to March 2022 at West China Hospital, West China Second University Hospital, First People's Hospital of Longquanyi District (West China Longquan Hospital), Meishan Women and Children's Hospital (Alliance Hospital of West China Second University Hospital) and Chengdu Jinjiang Hospital for Women and Children Health. Demographic and clinical characteristics of 263 pneumococcal disease (PD) all-age patients were collected and analyzed. The serotypes, sequence types (STs), and antibiotic resistance of the strains were determined by next-generation sequencing, sequence analysis and the microdilution broth method.

Results

The most common pneumococcal serotypes were 19F (17.87%), 19A (11.41%), 3 (8.75%), 23F (6.46%) and 6A (5.70%). Coverage rates for PCV10, PCV13, PCV15, PCV20 and PCV24 were 36.12, 61.98, 61.98, 63.12 and 64.26%, respectively. Prevalent STs were ST271 (12.55%), ST320 (11.79%), ST90 (4.18%), ST876 (4.18%) and ST11972 (3.42%). Penicillin-resistant S. pneumoniae (PRSP) accounted for 82.35 and 1.22% of meningitis and nonmeningitis PD cases, respectively. Resistance genes msrD (32.7%), mefA (32.7%), ermB (95.8%), tetM (97.3%) and catTC (7.6%) were found among 263 isolates. Most isolates showed high resistance to erythromycin (96.96%) and tetracycline (79.85%), with more than half being resistant to SXT (58.94%). A few isolates were resistant to AMX (9.89%), CTX (11.03%), MEN (9.13%), OFX (1.14%), LVX (1.14%) and MXF (0.38%). All isolates were susceptible to vancomycin and linezolid.

Conclusion

Our study provides reliable information, including the prevalence, molecular characterization and antimicrobial resistance of S. pneumoniae isolates causing pneumococcal diseases in Southwest China. The findings contribute to informed and clinical policy decisions for prevention and treatment.

Disease-Associated Streptococcus pneumoniae Genetic Variation

Streptococcus pneumoniae is an opportunistic pathogen that causes substantial illness and death among children worldwide. The genetic backgrounds of pneumococci that cause infection versus asymptomatic carriage vary substantially. To determine the evolutionary mechanisms of opportunistic pathogenicity, we conducted a genomic surveillance study in China. We collected 783 S. pneumoniae isolates from infected and asymptomatic children. By using a 2-stage genomewide association study process, we compared genomic differences between infection and carriage isolates to address genomic variation associated with pathogenicity. We identified 8 consensus k-mers associated with adherence, antimicrobial resistance, and immune modulation, which were unevenly distributed in the infection isolates. Classification accuracy of the best k-mer predictor for S. pneumoniae infection was good, giving a simple target for predicting pathogenic isolates. Our findings suggest that S. pneumoniae pathogenicity is complex and multifactorial, and we provide genetic evidence for precise targeted interventions.

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.

Integrated genomic analysis of antibiotic resistance and virulence determinants in invasive strains of Streptococcus pneumoniae

Introduction

Streptococcus pneumoniae is an important human pathogen that may cause severe invasive pneumococcal diseases (IPDs) in young children and the elderly. A comprehensive comparative whole-genome analysis of invasive and non-invasive serotype strains offers great insights that are applicable to vaccine development and disease control.

Methods

In this study, 58 invasive (strains isolated from sterile sites) and 71 non-invasive (serotypes that have not been identified as invasive in our study) pneumococcal isolates were identified among the 756 pneumococcal isolates obtained from seven hospitals in Zhejiang, China (2010-2022). Serotyping, antimicrobial resistance tests, and genomic analyses were conducted to characterize these strains.

Results and discussion

The three most invasive serotypes were 23F, 14, and 6B. The invasive pneumococcal isolates' respective resistance rates against penicillin, ceftriaxone, tetracycline, and erythromycin were 34.5%, 15.5%, 98.3%, and 94.7%. Whole-genome sequencing indicated that the predominant invasive clonal complexes were CC271, CC876, and CC81. The high rate of penicillin non-susceptible Streptococcus pneumoniae (PNSP) is related to the clonal distribution of resistance-conferring penicillin-binding proteins (PBP). Interestingly, we found a negative correlation between invasiveness and resistance in the invasive pneumococcal serotype strains, which might be due to the proclivity of certain serotypes to retain their β-lactam resistance. Moreover, the mutually exclusive nature of zmpC and rrgC+srtBCD suggests their intricate and potentially redundant roles in promoting the development of IPD. These findings reveal significant implications for pneumococcal vaccine development in China, potentially informing treatment strategies and measures to mitigate disease transmission.

[Dose-response relationship between age and Streptococcus pneumoniae vaccination coverage in kindergarten children]

OBJECTIVES

To investigate the potential relationship between age and Streptococcus pneumoniae vaccination coverage in kindergarten children, and to provide a basis for guiding vaccination and developing new protein vaccines.

METHODS

The stratified cluster random sampling method was used to select 1 830 healthy children from six kindergartens in Shunde District, Foshan City, China, and nasopharyngeal swabs were collected for the isolation and identification of Streptococcus pneumoniae. The logistic regression model based on restricted cubic spline was used to analyze the dose-response relationship between age and Streptococcus pneumoniae vaccination coverage.

RESULTS

The rate of nasal Streptococcus pneumoniae carriage was 22.46% (411/1 830) among the kindergarten children, with the predominant serotypes of 6B, 19F, 15A, 23A, 34, and 23F. The coverage rates of 10-valent pneumococcal conjugate vaccine (PCV10) and 13-valent pneumococcal conjugate vaccine (PCV13) were 53.0% and 57.9%, respectively, and there was a significant non-linear dose-response relationship between age and the coverage rates of PCV10 and PCV13 (P<0.05), with a higher coverage rate of PCV10 (88.0%) and PCV13 (91.1%) in the children aged 2 years. There was a significant non-linear dose-response relationship between age and the coverage rates of pilus islet 1 (PI-1) and pilus islet 2 (PI-2) (P<0.05), with a lower vaccination coverage rate for PI-1 (37.7%) and PI-2 (16.1%). The coverage rates of PI-1 (13.0%-58.5%) and PI-2 (6.0%-29.4%) were lower in all age groups. The virulence genes lytA (99.5%) and ply (99.0%) associated with candidate protein vaccines showed higher vaccination coverage rates.

CONCLUSIONS

There is a significant non-linear dose-response relationship between the age of kindergarten children and the coverage rates of PCV10 and PCV13 serotypes, and kindergarten children aged 2 years have a relatively high coverage rate of PCV. The high prevalence of the virulence genes lytA and ply shows that they are expected to become candidate virulence factors for the development of a new generation of recombinant protein vaccines.

Molecular characterization of invasive Streptococcus pneumoniae clinical isolates from a tertiary children's hospital in eastern China

Streptococcus pneumoniae is a common opportunistic pathogen that causes invasive pneumococcal disease (IPD), especially in children. This study aimed to determine the prevalence and molecular characteristics of S. pneumoniae isolated from children with IPD. A total of 78 S. pneumoniae isolates from aseptic body fluids of 70 IPD patients were collected at the Children's Hospital of Nanjing Medical University (Jiangsu Province, China) during 2017-2021. Whole-genome sequencing technology was used to analyze the serotype, sequence type (ST), virulence, and antibiotic resistance of the 78 invasive S. pneumoniae clinical isolates. Our results showed that the pneumococcal infection rate declined after the COVID-19 outbreak in 2019. Serotypes 19F, 14, 6A, 23F, 19A, and 6B were the most common strains. The pneumococcal conjugate vaccine (PCV) 13 serotype coverage rate was 87.1%. All isolates were classified by multi-locus sequence typing (MLST) analysis into 27 different STs, including 3 novel STs (ST17941, ST17942, and ST17944) and 1 novel allele [recP (558)]. The most predominant ST was ST271, followed by ST320 and ST876. All isolates carried the following virulence genes: cbpG, lytB, lytC, pce (cbpE), pavA, slrA, plr (gapA), hysA, nanA, eno, piuA, psaA, cppA, iga, htrA (degP), tig (ropA), zmpB, and ply. All isolates were multidrug resistant and had high levels of resistance to macrolides, tetracyclines, and sulfonamides. Taken together, this study revealed extensive genetic diversity among S. pneumoniae isolates from a single Chinese hospital. Wearing masks, universal infant vaccination with PCV13, and the launch of recombinant protein vaccine development programs could reduce the burden of IPD in children. IMPORTANCE Invasive pneumococcal disease (IPD) caused by Streptococcus pneumoniae in children remains a global burden and should be given more attention due to the fact that the pneumococcal vaccine is not fully covered globally. The molecular epidemiological characteristics of S. pneumoniae are not so clear, especially in these years of COVID-19. In this study, we collected S. pneumoniae isolates from the aseptic body fluid of children with IPD from 2017 to 2021 in a tertiary children's hospital in China and revealed the extensive genetic diversity of these isolates. Most importantly, we first found that the rate of pneumococcal infection has declined since the COVID-19 outbreak in 2019, which means that wearing masks could reduce the transmission of S. pneumoniae. In addition, it was shown that universal infant vaccination with PCV13 seems essential for reducing the burden of IPD in children.

Pan-Genome-Wide Association Study of Serotype 19A Pneumococci Identifies Disease-Associated Genes

Despite the widespread implementation of pneumococcal vaccines, hypervirulent Streptococcus pneumoniae serotype 19A is endemic worldwide. It is still unclear whether specific genetic elements contribute to complex pathogenicity of serotype 19A isolates. We performed a large-scale pan-genome-wide association study (pan-GWAS) of 1,292 serotype 19A isolates sampled from patients with invasive disease and asymptomatic carriers. To address the underlying disease-associated genotypes, a comprehensive analysis using three methods (Scoary, a linear mixed model, and random forest) was performed to compare disease and carriage isolates to identify genes consistently associated with disease phenotype. By using three pan-GWAS methods, we found consensus on statistically significant associations between genotypes and disease phenotypes (disease or carriage), with a subset of 30 consistently significant disease-associated genes. The results of functional annotation revealed that these disease-associated genes had diverse predicted functions, including those that participated in mobile genetic elements, antibiotic resistance, virulence, and cellular metabolism. Our findings suggest the multifactorial pathogenicity nature of this hypervirulent serotype and provide important evidence for the design of novel protein-based vaccines to prevent and control pneumococcal disease. IMPORTANCE It is important to understand the genetic and pathogenic characteristics of S. pneumoniae serotype 19A, which may provide important information for the prevention and treatment of pneumococcal disease. This global large-sample pan-GWAS study has identified a subset of 30 consistently significant disease-associated genes that are involved in mobile genetic elements, antibiotic resistance, virulence, and cellular metabolism. These findings suggest the multifactorial pathogenicity nature of hypervirulent S. pneumoniae serotype 19A isolates and provide implications for the design of novel protein-based vaccines.

Phylogenomic insights into evolutionary trajectories of multidrug resistant S. pneumoniae CC271 over a period of 14 years in China

BACKGROUND

Streptococcus pneumoniae is a gram-positive opportunistic pathogen, and infection risks of S. pneumoniae can be profoundly augmented by its acquired multidrug-resistance (MDR). The rapid development of MDR in S. pneumoniae was attributed to the international dissemination of a small number of multidrug-resistant "clones." Clonal complex (CC) 271 is a prevalent MDR CC in the world and the most prevalent CC in China. However, the evolutionary trajectories of multidrug-resistant S. pneumoniae CC271 in China still are largely unknown.

METHODS

We investigated a collection of 1312 S. pneumoniae isolates collected from 28 tertiary hospitals in China from 2007 to 2020. Recombination prediction and recombination-masked phylogenetic analysis were combined to determine the population structure and mode of evolution of CC271. Data from the Global Pneumococcal Sequencing program (GPS) were combined to understand the global distribution of clones identified in this study. Bayesian analysis were recruited to analysis the evolutionary dynamics of dominant clones within CC271 in China.

RESULTS

The phylogenomic analysis resulted in the discovery of two globally distributed clones, ST271-A and ST271-B. ST271-A was a derivative of ST236 and an ancestor of ST271-B and ST320, refining the internal phylogenetic relationship of CC271. ST271-B was the most dominant clone in China, with higher β-lactam resistance especially for cephalosporins comparing to other MDR clones. Bayesian skyline plot showed a rapid expansion of 19F ST271-B from 1995 to 2000, which correlates with the widespread use of cephalosporins in the 1990s in China. 19A ST320, a vaccine-escape clone, is the second largest population in China. The Bayesian skyline plot showed that the 19A ST320 began to expand rapidly around 2001, which appeared to coincide with the prevalence of 19A after application of PCV7 in 2000 in the USA. We also observed frequent transmission of 19A ST320 between countries. It suggests that mass vaccination in some countries could affect the prevalence of clones in unvaccinated countries in the context of high-frequency international transmission.

CONCLUSIONS

Our results refined the internal phylogenetic relationship of CC271, showing that the 19F ST271-B and 19A ST320 evolved independently from ST271-A, with different histories and driving forces for their evolution and dissemination in China.

Construction and protective efficacy of a novel Streptococcus pneumoniae fusion protein vaccine NanAT1-TufT1-PlyD4

During the past decades, with the implementation of pneumococcal polysaccharide vaccine (PPV) and pneumococcal conjugate vaccines (PCVs), a dramatic reduction in vaccine type diseases and transmissions has occurred. However, it is necessary to develop a less expensive, serotype-independent pneumococcal vaccine due to the emergence of nonvaccine-type pneumococcal diseases and the limited effect of vaccines on colonization. As next-generation vaccines, conserved proteins, such as neuraminidase A (NanA), elongation factor Tu (Tuf), and pneumolysin (Ply), are promising targets against pneumococcal infections. Here, we designed and constructed a novel fusion protein, NanAT1-TufT1-PlyD4, using the structural and functional domains of full-length NanA, Tuf and Ply proteins with suitable linkers based on bioinformatics analysis and molecular cloning technology. Then, we tested whether the protein protected against focal and lethal pneumococcal infections and examined its potential protective mechanisms. The fusion protein NanAT1-TufT1-PlyD4 consists of 627 amino acids, which exhibits a relatively high level of thermostability, high stability, solubility and a high antigenic index without allergenicity. The purified fusion protein was used to subcutaneously immunize C57BL/6 mice, and NanAT1-TufT1-PlyD4 induced a strong and significant humoral immune response. The anti-NanAT1-TufT1-PlyD4 specific IgG antibody assays increased after the first immunization and reached the highest value at the 35th day. The results from in vitro experiments showed that anti-NanAT1-TufT1-PlyD4 antisera could inhibit the adhesion of Streptococcus pneumoniae (S. pneumoniae) to A549 cells. In addition, immunization with NanAT1-TufT1-PlyD4 significantly reduced S. pneumoniae colonization in the lung and decreased the damage to the lung tissues induced by S. pneumoniae infection. After challenge with a lethal dose of serotype 3 (NC_WCSUH32403), a better protection effect was observed with NanAT1-TufT1-PlyD4-immunized mice than with the separate full-length proteins and the adjuvant control; the survival rate was 50%, which met the standard of the marketed vaccine. Moreover, we showed that the humoral immune response and the Th1, Th2 and Th17-cellular immune pathways are involved in the immune protection of NanAT1-TufT1-PlyD4 to the host. Collectively, our results support that the novel fusion protein NanAT1-TufT1-PlyD4 exhibits extensive immune stimulation and is effective against pneumococcal challenges, and these properties are partially attributed to humoral and cellular-mediated immune responses.

Molecular Characterization of Predominant Serotypes, Drug Resistance, and Virulence Genes of Streptococcus pneumoniae Isolates From East China

Streptococcus pneumoniae is a common diplococcus pathogen found worldwide. The characterization of predominant serotypes, drug resistance, and virulence genes of S. pneumoniae isolates prevailing in different areas and countries is clinically important for choice of antibiotics and improvement of vaccines. In this study, pneumonia (78.7%) and meningitis (37.0%) were the predominant diseases observed in the 282 (children) and 27 (adults) S. pneumoniae-infected patients (p < 0.05) from seven hospitals in different areas of East China. Of the 309 pneumococcal isolates, 90.3% were classified by PCR into 15 serotypes, with serotypes 19F (27.2%) and the 6A/B (19.1%) being most predominant (p < 0.05). Importantly, serotypes 15A and 15B/C combined for a total of 10.4% of the isolates, but these serotypes are not included in the 13-valent pneumococcal capsule conjugate vaccine used in China. Antimicrobial susceptibility analysis by the E-test showed that >95% of the 309 pneumococcal isolates were susceptible to moxifloxacin and levofloxacin, as well as 18.4, 85.8, and 81.6% of the isolates displayed susceptibility to penicillin, cefotaxime, and imipenem, respectively. A significant correlation between the prevalence of predominant serotypes and their penicillin resistance was observed (p < 0.05). In particular, >95% of all the pneumococcal isolates showed resistance to erythromycin and azithromycin. Of the nine detected virulence genes, the lytA, ply, hysA, and nanA were the most common with 95–100% positive rates in the 309 pneumococcal isolates, while the pavA and psaA genes displayed a significant correlation with pneumococcal bacteremia and meningitis (p < 0.05). Overall, our data suggested that the predominant serotypes, drug resistance, and virulence genes of the S. pneumoniae isolates prevailing in East China are distinct from those observed in other areas of China and adjacent countries.