Lin Liu1†
Yanfei Wang2,3,4†Lihong Ge5†
Dongping Hu6Xi Xiang7
Ying Fu8,9Jun Lu10
Xi Li1
Yunsong Yu2,3,4Yuexing Tu1,11*‡
Xueqing Wu2,3,4*‡- 1Laboratory Medicine Center, Department of Clinical Laboratory, Zhejiang Provincial People’s Hospital, Affiliated People’s Hospital, Hangzhou Medical College, Hangzhou, Zhejiang, China
- 2Department of Infectious Diseases, Sir Run Run Shaw Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang, China
- 3Sir Run Run Shaw Hospital, Zhejiang University School of Medicine, Key Laboratory of Microbial Technology and Bioinformatics of Zhejiang Province, Hangzhou, Zhejiang, China
- 4Regional Medical Center for National Institute of Respiratory Diseases, Sir Run Run Shaw Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang, China
- 5Department of Clinical Laboratory, The Children’s Hospital, Zhejiang University School of Medicine, National Clinical Research Center for Child Health, Hangzhou, Zhejiang, China
- 6Department of Infectious Disease, Affiliated Dongyang Hospital of Wenzhou Medical University, Dongyang, Zhejiang, China
- 7Department of Clinical Laboratory, Affiliated Jinhua Hospital, Zhejiang University School of Medicine, Jinhua, Zhejiang, China
- 8Department of Clinical Laboratory, Sir Run Run Shaw Hospital, School of Medicine, Zhejiang University, Hangzhou, Zhejiang, China
- 9Key Laboratory of Precision Medicine in Diagnosis and Monitoring Research of Zhejiang Province, Hangzhou, Zhejiang, China
- 10The Quzhou Affiliated Hospital of Wenzhou Medical University, Quzhou People’s Hospital, Quzhou, Zhejiang, China
- 11Department of Critical Care Medicine, Tongde Hospital of Zhejiang Province, Hangzhou, China
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.
Introduction
Streptococcus pneumoniae is an important pathogen that causes invasive pneumococcal diseases (IPDs) such as invasive pneumonia, sepsis, bacteremia, and meningitis (Weiser et al., 2018). The most infected population comprises young children, older adults, and immunocompromised patients (O’Brien et al., 2009). IPDs can be effectively controlled by introducing pneumococcal conjugate vaccines (PCVs) which are designed by target the most prevalent and virulent pneumococcal serotypes (Whitney et al., 2003; O’Brien et al., 2009; Pneumococcal Disease: Prevention | CDC, 2023). Epidemiological studies of IPD have consistently reported specific serotypes that cause invasive diseases, such as the recently reported serotype 4 and 24F IPD emergence in Israel and France, respectively (Kellner et al., 2021; Lo et al., 2022). To identify important virulent strains, a comparative investigation between isolates from different origins is widely accepted (Sharew et al., 2021; Yan et al., 2021). However, there is a potential risk of overlooking important differences between invasive and non-invasive serotype strains. Investigations of this type appear to be relatively scarce in the current literature, which are crucial for understanding the disease-causing potential of invasive pneumococcal serotypes and for contributing to vaccination strategies and vaccine development.
The selection of anti-infective treatment regimens should consider the drug resistance status of the strain. The increasing resistance of S. pneumoniae to β-lactams, macrolides, and tetracyclines has led to increasingly limited options for treating IPD (Thummeepak et al., 2015; Suaya et al., 2020). The high resistance of S. pneumoniae to macrolides and tetracyclines in China means that these drugs are seldom prescribed, and β-lactams are the first-line treatment for IPDs (Zhou et al., 2022b). However, macrolides are still recommended for regions with resistance rates below 25% (Gregory and Davis, 2020), and the global increase in the non-susceptibility S. pneumoniae to β-lactams highlights the critical need to monitor the antibiotic resistance status of invasive pneumococci.
The development of IPD usually starts with bacterial upper respiratory colonization, in which pneumococcus may asymptomatically colonize the host (Bogaert et al., 2004). If pneumococcal strains evade the host immune defense and reach the lower respiratory tract, they can cause inflammation and fluid accumulation in the lungs, which is normally confirmed as pneumonia (Coonrod, 1989; Bogaert et al., 2004). Once the pneumococcus spreads to the bloodstream or cerebrospinal fluid, it causes sepsis or meningitis, which can lead to organ failure and death (Mook-Kanamori et al., 2011; Weiser et al., 2018). During this process, various pneumococcal virulence factors contribute to immune evasion (capsule polysaccharides, Cps), epithelial cell adhesion (pneumococcal rrg pathogenic island), and host tissue invasion (zinc metalloproteinase, ZmpC) (Jonsson et al., 1985; Camilli et al., 2006; El Mortaji et al., 2012; Cremers et al., 2014). Comparative whole-genome sequencing (WGS) analysis of virulence factors among invasive and non-invasive serotype strains, for example the presence of critical virulence genes in invasive isolates rather than non-invasive serotype strains, would be a valuable contribution to our understanding of IPD development.
In this study, we identified 58 invasive pneumococcal strains which were isolated form sterile site from a pool of 756 isolates collected from seven hospitals in Zhejiang, China, during the period 2010–2022. Those 58 invasive pneumococcal cover 16 serotypes. A genomic comparison analysis was conducted between these strains and a set of 71 strains from serotypes that had not been previously identified as invasive in our collected cases. This study aimed to investigate the serotype distribution, antimicrobial susceptibility, molecular epidemiology, and virulence factors of invasive pneumococcal serotype strains in Zhejiang Province, China.
Materials and methods
Streptococcus pneumoniae isolation and serotyping
S. pneumoniae isolates (n=129) were collected from seven tertiary hospitals in Zhejiang, China, from July 2010 to January 2022, and included 58 invasive and 71 non-invasive serotype strains. Invasive isolates were collected from patients’ sterile site, for instance, blood, bronchoalveolar lavage fluid (BALF), and cerebrospinal fluid (CSF) specimens. Non-invasive serotype isolates were those strains obtained from sputum, nasopharynx, and oropharynx specimens and belong to serotypes that have never been identified in IPDs in all of our pneumococcal-positive cases. All isolates were obtained by culturing the clinical samples on blood agar plates at 37°C with 5% CO2 and were identified by optochin, bile solubility, and lytA PCR tests. Thereafter, all isolates were subjected to serotyping by the latex agglutination test and Quellung reaction (SSI Diagnostica, Denmark). We also conducted in silico serotyping after WGS (detailed below) using SeroBA software (https://github.com/sanger-pathogens/seroba) (Epping et al., 2018).
Clinical information collection
The clinical information of all patients in each IPD case was retrospectively extracted from their medical records with the approval of the Sir Run Shaw Hospital Ethics Review Committee (Zhejiang University School of Medicine, 20201112-32), which included data on sex, age, and primary diagnosis (Table 1). The specimen types of all cases are summarized in Table 2 (IPD cases) and Supplementary Table 1 (non-invasive pneumococcal infection cases).
Table 1 Demographic and clinical characteristics of IPD patients.
Table 2 Invasive specimen type and related pneumococcal diseases.
Antimicrobial susceptibility test
Broth microdilution assays were used to determine the minimal inhibitory concentrations (MICs) of the tested antimicrobial agents according to the Clinical and Laboratory Standard Institute (CLSI) protocols, as described previously (Lewis, 2023). The antimicrobial agents tested were penicillin (PEN), ceftriaxone (CRO), erythromycin (ERY), and tetracycline (TET). S. pneumoniae strain ATCC 49619 was used as the quality control strain. The results were defined according to the 2023 Clinical and CLSI Guidelines M100-Ed33 (Lewis, 2023).
WGS and analysis
Genomic DNA was extracted from all isolates using a QIAamp® DNA Mini Kit (Qiagen, Valencia, CA, USA) and sequenced using an Illumina HiSeq X 10 platform (Illumina, San Diego, CA, USA). The Illumina reads were assembled by end pairing using Shovill (Seemann T, https://github.com/tseemann/shovill), with a minimum splicing length of 200 bp and a minimum coverage of 10-fold. The final assemblies have the N50 not less than 60K and the minimum sequencing depth is 300X. In silico sequence type (ST) of each strain was then obtained by blasting our genome against the pubMLST database (https://pubmlst.org) via mlst (Seemann T, https://github.com/tseemann/mlst) (Jolley and Maiden, 2010). To confirm the serotyping results of our Quellung reactions, we also conducted in silico serotyping by comparing sequenced reads to a database containing key genes that determine serotypes in cps gene clusters using SeroBA. Virulence genes and antimicrobial resistance genes were screened using ABRicate software (Seemann T, Abricate, GitHub https://github.com/tseemann/abricate). Prominent amino acid substitutions in penicillin (PEN)-binding proteins (PBP1a, PBP2b, and PBP2x) were analyzed by BLST+ (2.13.0) (https://github.com/ncbi/blast_plus_docs) against the database of CDC, USA (https://www.cdc.gov/streplab/pneumococcus/mic.html). Phylogenetic trees were constructed using popPUNK (Lees et al., 2019) (https://github.com/bacpop/PopPUNK) and visualized using iTOL (v6, https://itol.embl.de).
Statistical analysis
To assess the difference in antibiotic resistance between invasive and non-invasive isolates, we conducted a Mann–Whitney U test on each column, where a p<0.0001 was considered statistically significant. The correlation between invasive and PEN insensitivity ratios for each serotype was calculated using a two-tailed method. A Pearson correlation coefficient (r)<-0.7 was considered a strong negative correlation, r<-0.5 was considered a moderate correlation, and r<-0.3 was considered a weak correlation. A two-tailed p<0.05 was considered statistically significant. All analyses were performed using GraphPad Prism v9.5.0.
Results
Clinical characteristics of patients with IPD
The demographic and clinical characteristics of patients with IPD are summarized in Tables 1, 2. Among 58 patients with IPD, males accounted for the majority (69.0%). In the different age groups, 44.8% (26/58) of IPD cases were identified in young children (0–5 years), and 19.0% (11/58) occurred in patients ≥65 years old. The three most common primary diagnoses were bronchitis, pneumonia, and sepsis. The specimen types used were blood (49/58, 84.5%), BALF (5/58, 8.6%), and cerebrospinal fluid (CSF, 4/58, 6.9%). Bronchitis, fever, and sepsis were the primary diagnoses of cases that later been confirmed as IPD due to pneumococcal blood culture positive. In cases of non-invasive pneumococcal infection, the most commonly primary diagnosed disease was pneumonia (Supplementary Table 1).
Serotype distribution
Among all invasive strains, 16 serotypes were identified, including 10 vaccine serotypes and six non-vaccine serotypes, which accounted for 84.5% (49/58) and 15.5% (9/58) of all invasive isolates, respectively (Figures 1A, B). The top three invasive serotypes were 23F (17.2%, 10/58), 14 (17.2%, 10/58), and 6B (10.3%, 6/58), all of which were covered by 7-valent pneumococcal polysaccharide conjugate vaccine (PCV7) (Figure 1A). The most frequently isolated non-vaccine-invasive serotypes were 34 and 15C (Figure 1C). Young children were the most frequently infected population with IPD, wherein especially the serotypes 6B, 14, and 23F caused IPD cases. We also included 71 non-invasive pneumococcal serotype strains in our study. The serotype distribution data showed that 6C, 15A, 15B, and 16F were the most frequently identified serotypes in cases of non-invasive pneumococcal infection (Figure 1D).
Figure 1 Serotype distribution and vaccine serotype coverage of pneumococcal isolates All tested S. pneumoniae isolates were serotyped by quellung reactions and in silico by whole genome sequencing via SeroAB. (A) The distribution of invasive pneumococcal vaccine serotypes in different age groups: 0-5 (blank), 6-65 (grey), and, >65 (light grey); (B) The proportion of pneumococcal conjugate vaccine (PCV) 7, PCV13-add, PCV20-add, and non-vaccine serotype (NVT) invasive pneumococcal isolates; (C) Serotype distribution of NVT invasive pneumococcal isolates; (D) Serotype distribution of pneumococcal isolates belong to the serotypes that are never been detected for invasive cases (non-invasive serotypes).
Antimicrobial susceptibility
Antimicrobial susceptibility results for the 58 invasive S. pneumoniae isolates are presented in Table 3. According to the non-meningitis breakpoint, the non-susceptibility rates of the isolates to PEN, CRO, TET, and ERY were 34.5%, 15.5%, 98.3%, and 94.7%, respectively. Regarding the meningitis breakpoint, the insensitivity rates of the isolates against PEN and CRO were 82.8% and 43.1%, respectively. Most PEN- and CRO-non-susceptible strains were PCV-covered serotypes, whereas serotypes 19F and 14 accounted for the majority. The MIC values of each isolate are presented in Supplementary Tables 3 and 4, showing that the MIC90 values of invasive and non-invasive pneumococcal isolates against PEN were 8 and 4 μg/mL, respectively.
Table 3 Antibiotic susceptibility of invasive S. pneumoniae.
Phylogenetic analysis and antimicrobial resistance gene determination
To further understand the genetic characteristics of the invasive S. pneumoniae strains, we conducted WGS analyses of all invasive pneumococcal isolates (n=58) and compared them with the non-invasive isolates (n=71). According to our phylogenetic analysis (Figure 2A), the most prevalent invasive clone clusters (CCs) were CC876 (15.5%, 9/58, serotype 14), CC271 (13.8%,8/58, serotype 19F/19A), and CC81 (12.1%, 7/58, serotype 23F). The major non-invasive S. pneumoniae clones were CC8250 (9.9%, 7/71, serotype 16F), CC99 (8.5%, 6/71, serotype 11A/35A), and CC5893 (8.5%, 6/71, serotype 29). The beta-lactam non-sensitive strains were mainly distributed in the invasive branches and were concentrated in the prevalent clones CC271, CC876, and CC81, which carry PBP1a-2b-2x types of 13-11-33, 60-16-312, and 15-12-18, respectively. Correspondingly, a comparison of MIC values between invasive and non-invasive strains showed significantly higher MICs in invasive pneumococcus against PEN and CRO (Figure 2B). We also identified a novel PBP2x type in strain dy19012 (ST342, serotype 23F) and a novel PBP2b type in strain sy19004 (ST81, serotype 23F) that may mediate PEN resistance (MIC=8 μg/mL). Regardless of their invasiveness, all strains presented high-level resistance against ERY and TET; most carried genes of ermB and tetM, and the carry rate for having both ermB and tetM is 90.6% (116/128). Another macrolide resistance-encoding gene, mef(A) was prevalent mainly in CC27, and the carry rate for having both mefA and tetM is 7.8% (10/128). The genome of hz21002 was obtained directly by sequencing the blood sample and did not pass quality control, resulting in unavailability of data for this isolate.
Figure 2 The phylogenetic tree and antibiotic resistance determinants of invasive and non-invasive serotype pneumococcal strains. (A) The phylogenetic tree of all tested pneumococcal isolates (n=129) was constructed in PopPUNK, where three major clone complexes (CC) for invasive and non-invasive serotype isolates were shaded in light red and blue, respectively. The metadata including specimen types, sequence type (ST), serotype, and antibiotic susceptibility test results (AST) was aligned for all sequenced pneumococcal isolates, which was followed by the detection of antibiotic resistance determinants of PBPs, tetM, ermB, and mefA; (B) Minimal inhibition concentration (MIC) comparison between invasive and non-invasive serotype strains. “****” indicate a significant difference with a p-value less than 0.0001, “ns” indicate no significance was detected, red dash line is the cut-off of resistant MIC.
Correlation of beta-lactam resistance and invasion
Careful analysis of the above datasets and taking into account our total pneumococcal strain storage (n=756) revealed that the serotype 19F strains presented the highest PEN insensitivity rate of 62.8% (108/172) and the lowest invasive ratio of 2.9% (5/172). In contrast, 80% (4/5) of the serotype 4 strains in our study were invasive isolates and none were PEN-insensitive (Figures 3A, B). This finding prompted us to perform correlation analysis for all invasive serotypes. As shown in Figure 3C, the invasive and PEN insensitivity ratios of the invasive serotypes were negatively and moderately correlated, respectively (r=-0.5444, p=0.292). Among all invasive serotypes, 19F and 19A were the most resistant and least invasive, while serotypes 4 and 8 were the most invasive and least resistant, respectively. Because 19F is the most prevalent strain in our study resulting in its high isolation rate in the invasive strains, we conducted a separate analysis that included all available genome data of serotype 19F strains (n=124) in our laboratory. As shown in Supplementary Figure 1, most serotype 19F isolates belonged to serotype CC271. The PEN insensitivity of this serotype is very high (62.8%), and is mediated by the same PBP1a-2b-2x type (13-11-33). Among more than hundred 19F strains, only five were invasive.
Figure 3 The correlation between invasiveness and penicillin non-susceptibility of invasive serotypes. The proportion of invasive (A) and PEN non-susceptible (B) isolate in each invasive serotype (included all our strain bank isolates, n=756); (C) The correlation analysis between the invasiveness and PEN non-susceptibility in all invasive serotypes. A Pearson correlation coefficient (r) <-0.7 was considered a strong negative correlation, r<-0.5 was considered a moderate correlation, and r<-0.3 was considered a weak correlation. A two-tailed p<0.05 was considered statistically significant.
Virulence factors analysis
To determine the virulence of invasive pneumococci, we conducted virulence gene screening for both invasive and non-invasive strains. As shown in Figure 4, the choline-binding protein gene cbpA was present in only one strain. The virulence factor-encoding genes pitAB (iron uptake) and srtG1,2 (cognate sortase) were only present in the CC271 strains. The genes cpsA (capsule synthesis), hysA (hyaluronidase), lytABC (autolysin), nanAB (neuraminidase), pavA (fibronectin-binding protein), pce (phosphorylcholine esterase), ply (pneumolysin), and psaA (pneumococcal surface adhesin A) were present in both invasive and non-invasive strains with no clone specificity. The non-invasive strains carried slightly more cbpGD, pfbA (plasmin- and fibronectin-binding protein A), and pspAC (pneumococcal surface protein A) genes. Moreover, the genes rrgABC (pilus adhesin), srtBCD (sortase), and zmpC (zinc metalloproteinase) tended to be carried by invasive strains. Interestingly, the co-occurrence of zmpC and rrgC+srtBCD was mutually exclusive, and zmpC was mainly carried by serotypes 4 and 8, whereas rrgC+srtBCD appeared more frequently in the 19F and 19A strains. Again, the analysis targeting only serotype 19F showed rrgABC+srtBCD was clonally distributed in these isolates (Supplementary Figure 2).
Figure 4 The virulence factor analysis of invasive and non-invasive serotype pneumococcal strains. (A) The detection of virulence factors was attached for each strain in the phylogenetic tree of invasive and non-invasive serotype strains. (B) Carrying proportion of virulence factors for invasive (black) and non-invasive (grey) serotype strains. Several virulence factors are clonally distributed (cbpA and pitAB.) or carried by all strains (ply and psaA), which were excluded from the analysis. (C) Distribution and proportion of zmpC in the carried serotypes of 3, 4, 8, and 14; D. Distribution and proportion of rrgC+srtBCD in the carried serotypes of 6A, 6B, 15C, 19A, 19F, and 34.
Discussion
The identification of prevalent invasive serotypes is crucial for effective prevention and management of patients with IPD. The data presented in this study indicate that the most prevalent IPD serotypes in Zhejiang, China, are 23F, 14, and 6B. However, the CSF-isolated pneumococcal serotypes were 34, 8, and 15C, of which 34 and 15C were the most commonly identified non-vaccine serotypes. Owing to the low PCV coverage in China, the prevalence of IPD serotypes is different from that in developed countries (Pick et al., 2020; Suaya et al., 2020; Yanagihara et al., 2021; Wu et al., 2022). A recent multicenter study of 300 invasive S. pneumoniae isolates predominantly from northern China indicated that the most prevalent invasive serotypes are 23F, 19F, and 19A, and the top CSF-isolated serotypes are 23F, 19F, and 14 (Zhou et al., 2022a). A similar study from western China showed that the most prevalent pneumococcal serotypes are 19F, 19A, and 6B (Yan et al., 2021). The inconsistent results found in the literature as well as in our study highlight the variation in serotype distribution across different regions of China. Therefore, it is crucial to conduct a national surveillance study on invasive S. pneumoniae to accurately assess this situation. The region-specific dominance of certain serotypes may direct different vaccine strategies. While PCV13 and PCV20 cover the majority of IPD serotypes in Zhejiang Province, non-vaccine serotypes 34 and 15C have caused highly virulent cases of meningitis that require continued vigilance.
Our AST results showed that the nonsusceptibility to penicillin and ceftriaxone in IPD isolates was significantly higher than that in non-invasive serotype strains. But further analysis indicates that β-lactams resistance was mainly mediated by the clonal distribution of CC271, CC876, and CC81 strains carrying certain PBP1a-2b-2x combinations which were confirmed as PEN-non-susceptible PBP types from a global pneumococcal database (Li et al., 2016). For this reason, according to the non-meningitis breakpoint, we reported a very high PNSP rate of over 30% for invasive S. pneumoniae tested in the current study. A study conducted in China that collected 993 strains of S. pneumoniae up to 2017 showed that the non-susceptibility rate of IPD strains to penicillin reached 22.4%, which was significantly higher than that of the non-invasive strains (Yan et al., 2021). However, another pneumococcal epidemiological study collected 300 invasive pneumococci from 2010 to 2015 and reported a very low PNSP rate (4.3%) (Zhou et al., 2022b). Later, the same group performed a deeper analysis of PBPs using conventional PCR and showed that the PBP site substitutions in PNSP strains matched our findings of non-susceptible PBP, which were distributed mostly in prevalent clones. Although in silico predicted MIC of strains carrying non-susceptible PBPs against penicillin was 4 μg/mL (Li et al., 2016), most of the isolates in their study had an MIC of 1–2 μg/mL. The reason for these inconsistent results might be due to MIC value interpretation, sample collection time periods, and regional differences. Macrolide and tetracycline resistance was maintained at high levels in both invasive and non-invasive pneumococci, which is mediated by the national distribution of ermB and tetM in China (Zhou et al., 2022b), indicating the limited clinical value of such drugs in China.
Notably, some serotype strains had very high invasive rates, but none of them presented a penicillin-resistant phenotype; therefore, we performed further analysis to examine the relationship between invasiveness and penicillin resistance. Among all the collected strains of invasive serotypes, we found a negative correlation between invasive and penicillin non-susceptibility rates. To date, no relevant reports have been published regarding S. pneumoniae. However, in the gram-negative bacterium Klebsiella pneumoniae, hypervirulent strains have been reported to exhibit a relatively low ability to acquire antibiotic-resistant plasmids and hardly simultaneously exhibit virulence and resistance (Tian et al., 2022). An ongoing study in our laboratory has shown that serotype 19F strains display greater proclivity to maintain their β-lactam resistance phenotype than that of other strains. In the top three invasive serotypes (4, 8, 33B), none of the strains to been found resistant to PEN in our study. Recent reports demonstrated the same PSSP for 190 invasive serotype 4 strains and 90 invasive serotype 8 strains (Hansen et al., 2021; Kellner et al., 2021). No report of PEN resistance was found in pneumococcal serotype 33B strains. The mechanism behind such phenotype need further identification. The negative correlation observed between the invasiveness and non-susceptibility of S. pneumoniae may be attributed to the presence and retention of specific virulence factors and resistance determinants in the highly invasive and resistant serotype strains, respectively. A more comprehensive understanding of this newly discovered epidemiological phenomenon requires further investigation into its underlying mechanisms.
During the invasive infection process, various virulence factors participate at different stages. We did not find a carrying difference between invasive and non-invasive pneumococcal serotype strains in several classical virulence factor-encoding genes, such as capsule synthesis, pneumococcal surface adhesin, autolysin, fibronectin-binding protein, and pneumolysin, which are similar to those reported previously (Yan et al., 2021). However, rrgABC (pilus adhesin), srtBCD (sortase), and zmpC (zinc metalloproteinase) are mostly carried by invasive pneumococci, and the rrgABC+srtBCD locus is clonally distributed in CC271 strains. There are two intriguing points in this section regarding these results. First, except for CC271, one of the pilus type 1 (P1) encoding genes, rrgC, would be independently expressed in the invasive strains together with the three P1 specific sortase encoding genes, srtBCD. RrgC has been reported as a lectin targeting different host glycosylations; however, it is the least understood pilus protein in S. pneumoniae (Day et al., 2017). Our findings demonstrate the importance of sortase for pilus expression and the small pili RrgC for S. pneumoniae causing invasive disease. Future pneumococci-host interaction studies would be valuable by focusing on these virulence factors. Furthermore, rrgABC+srtBCD was initially identified in the highly invasive serotype 4 strain TIGR4 (Barocchi et al., 2006). However, none of our serotype 4 invasive isolates was found to carry P1-related genes; instead, they encoded a zinc metalloproteinase. The co-occurrence of zmpC and rrgC+srtBCD was mutually exclusive, and zmpC was mainly carried by serotypes 4 and 8, which are highly correlated with IPDs. A study conducted in Italy (Camilli et al., 2006) demonstrated that serotypes 8 and 11A carried zmpC, and another study from the Netherlands (Cremers et al., 2014) revealed that zmpC was predominantly associated with serotypes 8, 4, 33A/F, and 11A/D. In contrast, our laboratory isolates of serotype 11A were collected from patients without IPD and no zmpC was detected in these strains. Pneumococcal ZmpC is involved in the breakdown of host tissues, whereas the pilus biogenesis proteins RrgC and SrtBCD play a role in the adhesion and colonization of the bacterium. There is no clear explanation for the mutual exclusion of these two different virulence factors; they may interfere with the expression of other functions or they are carried on different mobile genetic elements that are not transferable between different clones. Nevertheless, our findings demonstrate that both are important for the invasiveness of different invasive pneumococcal serotypes. Research on the serotype-specific virulence mechanism of S. pneumoniae would be significantly more meaningful than studying its overall pneumococcal virulence.
In conclusion, this study provides valuable insights into invasive S. pneumoniae serotypes and their management. A national surveillance study is necessary to understand the variation in invasive serotype distribution across different regions of China, which is crucial for PCV vaccination strategies. A high PNSP rate was related to the clonal distribution of non-susceptible PBP types, and we found a negative correlation between invasiveness and resistance in invasive pneumococcal strains. The mutually exclusive nature of zmpC and rrgC+srtBCD suggests their intricate and potentially redundant roles in promoting the development of IPD. Further mechanistic studies will contribute to the development of pneumococcal vaccines.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/genbank/, PRJNA977823 https://www.ncbi.nlm.nih.gov/genbank/, PRJNA795524 https://www.ncbi.nlm.nih.gov/genbank/,PRJNA924107.
Author contributions
XW, YT and YY conceptualized and designed the study, drafted the initial manuscript, and reviewed and revised the manuscript. LL and YW contributed equally to the pneumococcal isolates collection, carried out the initial analyses, and reviewed and revised the manuscript. LG, DH, XX, JL, YF, and XL, contributed to clinical data collection, pneumococcal serotyping, and genomic sequencing, and reviewed and revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
The National Natural Science Foundation of China (No.32000092); Jinhua Science and Technology Research Key Project 2021-3-070; Specific Fund (2022ZYC-B21) for Clinical Medical Research of the Zhejiang Medical Association; Quzhou Science and Technology Bureau, 2022K60.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2023.1238693/full#supplementary-material
Supplementary Figure 1 | The phylogenetic tree and antibiotic resistance determinants of pneumococcal serotype 19F strains. The phylogenetic tree of all sequenced pneumococcal serotype 19F isolates (n=124) was constructed in PopPUNK, where the major clone complex (CC) 271 was shaded in light red. The metadata including specimen types, sequence type (ST), antibiotic susceptibility test results (AST), and PBPs types were aligned for all isolates.
Supplementary Figure 2 | Virulence factor detection in serotype 19F strains. The detection of virulence factors was attached for all isolates in the phylogenetic tree along with specimen data and sequence type of pneumococcal serotype 19F strains (n=124). A similar virulence factor carrying pattern was noticed in this clone.
References
Barocchi, M. A., Ries, J., Zogaj, X., Hemsley, C., Albiger, B., Kanth, A., et al. (2006). A pneumococcal pilus influences virulence and host inflammatory responses. Proc. Natl. Acad. Sci. U.S.A. 103, 2857–2862. doi: 10.1073/pnas.0511017103
Bogaert, D., De Groot, R., Hermans, P. W. M. (2004). Streptococcus pneumoniae colonisation: the key to pneumococcal disease. Lancet Infect. Dis. 4, 144–154. doi: 10.1016/S1473-3099(04)00938-7
Camilli, R., Pettini, E., Grosso, M. D., Pozzi, G., Pantosti, A., Oggioni, M. R. (2006). Zinc metalloproteinase genes in clinical isolates of Streptococcus pneumoniae: association of the full array with a clonal cluster comprising serotypes 8 and 11A. Microbiol. (Reading) 152, 313–321. doi: 10.1099/mic.0.28417-0
Cremers, A. J. H., Kokmeijer, I., Groh, L., de Jonge, M. I., Ferwerda, G. (2014). The role of ZmpC in the clinical manifestation of invasive pneumococcal disease. Int. J. Med. Microbiol. 304, 984–989. doi: 10.1016/j.ijmm.2014.06.005
Day, C. J., Paton, A. W., Harvey, R. M., Hartley-Tassell, L. E., Seib, K. L., Tiralongo, J., et al. (2017). Lectin activity of the pneumococcal pilin proteins. Sci. Rep. 7, 17784. doi: 10.1038/s41598-017-17850-9
El Mortaji, L., Fenel, D., Vernet, T., Di Guilmi, A. M. (2012). Association of RrgA and RrgC into the Streptococcus pneumoniae pilus by sortases C-2 and C-3. Biochemistry 51, 342–352. doi: 10.1021/bi201591n
Epping, L., van Tonder, A. J., Gladstone, R. A., Consortium The Global Pneumococcal Sequencing, Bentley, S. D., Page, A. J., et al. (2018). SeroBA: rapid high-throughput serotyping of Streptococcus pneumoniae from whole genome sequence data. Microbial Genomics 4. doi: 10.1099/mgen.0.000186
Gregory, O., Davis, A. M. (2020). Diagnosis and treatment of adults with community-acquired pneumonia. JAMA 323, 885–886. doi: 10.1001/jama.2019.21118
Hansen, C. B., Fuursted, K., Valentiner-Branth, P., Dalby, T., Jørgensen, C. S., Slotved, H.-C. (2021). Molecular characterization and epidemiology of Streptococcus pneumoniae serotype 8 in Denmark. BMC Infect. Dis. 21, 421. doi: 10.1186/s12879-021-06103-w
Jolley, K. A., Maiden, M. C. (2010). BIGSdb: Scalable analysis of bacterial genome variation at the population level. BMC Bioinf. 11, 595. doi: 10.1186/1471-2105-11-595
Jonsson, S., Musher, D. M., Chapman, A., Goree, A., Lawrence, E. C. (1985). Phagocytosis and killing of common bacterial pathogens of the lung by human alveolar macrophages. J. Infect. Dis. 152, 4–13. doi: 10.1093/infdis/152.1.4
Kellner, J. D., Ricketson, L. J., Demczuk, W. H. B., Martin, I., Tyrrell, G. J., Vanderkooi, O. G., et al. (2021). Whole-genome analysis of Streptococcus pneumoniae serotype 4 causing outbreak of invasive pneumococcal disease, Alberta, Canada. Emerg. Infect. Dis. 27, 1867–1875. doi: 10.3201/eid2707.204403
Lees, J. A., Harris, S. R., Tonkin-Hill, G., Gladstone, R. A., Lo, S. W., Weiser, J. N., et al. (2019). Fast and flexible bacterial genomic epidemiology with PopPUNK. Genome Res. 29, 304–316. doi: 10.1101/gr.241455.118
Lewis, J. S. (2023). M100:Performance Standards for Antimicrobial Susceptibility Testing, 33rd Edition. Available at: https://clsi.org/standards/products/microbiology/documents/m100/ (Accessed April 1, 2023).
Li, Y., Metcalf, B. J., Chochua, S., Li, Z., Gertz, R. E., Walker, H., et al. (2016). Penicillin-binding protein transpeptidase signatures for tracking and predicting β-lactam resistance levels in Streptococcus pneumoniae. mBio 7, e00756–e00716. doi: 10.1128/mBio.00756-16
Lo, S. W., Mellor, K., Cohen, R., Alonso, A. R., Belman, S., Kumar, N., et al. (2022). Emergence of a multidrug-resistant and virulent Streptococcus pneumoniae lineage mediates serotype replacement after PCV13: an international whole-genome sequencing study. Lancet Microbe 3, e735–e743. doi: 10.1016/S2666-5247(22)00158-6
Mook-Kanamori, B. B., Geldhoff, M., van der Poll, T., van de Beek, D. (2011). Pathogenesis and pathophysiology of pneumococcal meningitis. Clin. Microbiol. Rev. 24, 557–591. doi: 10.1128/CMR.00008-11
O’Brien, K. L., Wolfson, L. J., Watt, J. P., Henkle, E., Deloria-Knoll, M., McCall, N., et al. (2009). Burden of disease caused by Streptococcus pneumoniae in children younger than 5 years: global estimates. Lancet 374, 893–902. doi: 10.1016/S0140-6736(09)61204-6
Pick, H., Daniel, P., Rodrigo, C., Bewick, T., Ashton, D., Lawrence, H., et al. (2020). Pneumococcal serotype trends, surveillance and risk factors in UK adult pneumonia 2013-18. Thorax 75, 38–49. doi: 10.1136/thoraxjnl-2019-213725
Pneumococcal Disease: Prevention | CDC. (2023). Available at: https://www.cdc.gov/pneumococcal/about/prevention.html (Accessed March 25, 2023).
Sharew, B., Moges, F., Yismaw, G., Mihret, A., Abebe, W., Fentaw, S., et al. (2021). Serotype distribution of Streptococcus pneumoniae isolates causing invasive and non-invasive infections using whole-genome sequencing in Ethiopia. Infect. Drug Resist. 14, 787–794. doi: 10.2147/IDR.S293578
Suaya, J. A., Mendes, R. E., Sings, H. L., Arguedas, A., Reinert, R.-R., Jodar, L., et al. (2020). Streptococcus pneumoniae serotype distribution and antimicrobial nonsusceptibility trends among adults with pneumonia in the United States 2009-2017. J. Infect. 81, 557–566. doi: 10.1016/j.jinf.2020.07.035
Thummeepak, R., Leerach, N., Kunthalert, D., Tangchaisuriya, U., Thanwisai, A., Sitthisak, S. (2015). High prevalence of multi-drug resistant Streptococcus pneumoniae among healthy children in Thailand. J. Infect. Public Health 8, 274–281. doi: 10.1016/j.jiph.2014.11.002
Tian, D., Liu, X., Chen, W., Zhou, Y., Hu, D., Wang, W., et al. (2022). Prevalence of hypervirulent and carbapenem-resistant Klebsiella pneumoniae under divergent evolutionary patterns. Emerg. Microbes Infect. 11, 1936–1949. doi: 10.1080/22221751.2022.2103454
Weiser, J. N., Ferreira, D. M., Paton, J. C. (2018). Streptococcus pneumoniae: transmission, colonization and invasion. Nat. Rev. Microbiol. 16, 355–367. doi: 10.1038/s41579-018-0001-8
Whitney, C. G., Farley, M. M., Hadler, J., Harrison, L. H., Bennett, N. M., Lynfield, R., et al. (2003). Decline in invasive pneumococcal disease after the introduction of protein-polysaccharide conjugate vaccine. N Engl. J. Med. 348, 1737–1746. doi: 10.1056/NEJMoa022823
Wu, X., Zhao, S., Jiang, Y., Xiang, X., Ge, L., Chen, Q., et al. (2022). Effect of pneumococcal conjugate vaccine availability on Streptococcus pneumoniae infections and genetic recombination in Zhejiang, China from 2009 to 2019. Emerging Microbes Infections 11, 606–615. doi: 10.1080/22221751.2022.2040921
Yan, Z., Cui, Y., Huang, X., Lei, S., Zhou, W., Tong, W., et al. (2021). Molecular characterization based on whole-genome sequencing of Streptococcus pneumoniae in children living in southwest China during 2017-2019. Front. Cell Infect. Microbiol. 11. doi: 10.3389/fcimb.2021.726740
Yanagihara, K., Kosai, K., Mikamo, H., Mukae, H., Takesue, Y., Abe, M., et al. (2021). Serotype distribution and antimicrobial susceptibility of Streptococcus pneumoniae associated with invasive pneumococcal disease among adults in Japan. Int. J. Infect. Dis. 102, 260–268. doi: 10.1016/j.ijid.2020.10.017
Zhou, M., Wang, L., Wang, Z., Kudinha, T., Wang, Y., Xu, Y., et al. (2022a). Molecular characterization of penicillin-binding protein2x, 2b and 1a of Streptococcus pneumoniae causing invasive pneumococcal diseases in China: A multicenter study. Front. Microbiol. 13. doi: 10.3389/fmicb.2022.838790
Keywords: invasive pneumococcal serotype, whole-genome sequencing, antibiotic resistance, rrgC, ZmpC
Citation: Liu L, Wang Y, Ge L, Hu D, Xiang X, Fu Y, Lu J, Li X, Yu Y, Tu Y and Wu X (2023) Integrated genomic analysis of antibiotic resistance and virulence determinants in invasive strains of Streptococcus pneumoniae. Front. Cell. Infect. Microbiol. 13:1238693. doi: 10.3389/fcimb.2023.1238693
Received: 12 June 2023; Accepted: 03 October 2023;
Published: 19 October 2023.
Edited by:
Yi-Wei Tang, Cepheid, United StatesReviewed by:
Amanda J. Chase, Nova Southeastern University, United StatesKaifeng Wu, Zunyi Medical University Third Affiliated Hospital, China
Jorge Eugenio Vidal, University of Mississippi Medical Center, United States
Copyright © 2023 Liu, Wang, Ge, Hu, Xiang, Fu, Lu, Li, Yu, Tu and Wu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xueqing Wu, eHVlcWluZy53dUB6anUuZWR1LmNu; Yuexing Tu, dHV5dWV4aW5nMTk4OEAxNjMuY29t
†These authors have contributed equally to this work
‡These authors have contributed equally to this work and share senior authorship