Molecular epidemiology of Streptococcus Pneumoniae isolated from children with community-acquired pneumonia under 5 years in Chengdu, China

DOI: https://doi.org/10.21203/rs.3.rs-1464814/v1

Abstract

Background: Streptococcus pneumoniae (S. Pneumoniae) is one of the most common community-associated pathogens responsible for pneumonia in children. This retrospective study aimed to investigate the molecular characteristics of S. Pneumoniae isolated from children with community-acquired pneumonia (CAP) under 5 years in Chengdu, China.

Methods: Molecular characteristics of S. Pneumoniae included serotype and virulence factor performed by using PCR method and sequence types (STs) determined by sequencing seven housekeeping genes. In addition, the potential relationships between molecular characteristics were depicted by minimum spanning tree and correspondence analysis.

Results: The prevailing serotypes were 19F (18.52%), 6B (17.59%), 19A (13.89%), 6A (6.48%) and 23F (5.56%) among 108 isolates. The overall coverage rates of 7-valent, 10-valent, 13-valent pneumococcal conjugate vaccines (PCVs) were 47.32%, 48.1%, 75%, respectively. Meanwhile, the coverage rates of PCV13 among the isolates from CAP patients aged <1 year were high up to 84.2%. MLST analysis results showed that there were 56 different STs identified, of which the dominant STs were ST271 (22.22%) and ST320 (12.04%). Additionally, most of the isolates carried ply, psaA, nanA, pavA, piaA, and CC271 isolates expressed more of nanA than non-CC271 isolates. Moreover, there were strong relevant relationships among STs, serotypes, virulence factors.

Conclusion: Considering serotypes and virulence factors together can be used as the foundation for the formulation of vaccine strategy.

1. Introduction

Community-acquired pneumonia (CAP) is a common cause of respiratory leading to high mortality in children under 5 years old all over the world, especially in low-income and developing countries. According to the World Health Organization (WHO) report, despite the number of death declining from 920,000 in 2016 to 740,000 in 2019, CAP was still the deadliest infectious disease accounting for 14% of all deaths in children younger than 5 years[1]. The etiology of CAP is complex constantly caused by various microorganisms, such as bacteria, viruses. It was reported by Ning et al[2] on the etiology of nearly 100,000 CAP children, Streptococcus pneumoniae (S. pneumoniae) was ranked top 3 among bacteria isolated from children aged < 5 years hospitalized with CAP, indicating S. pneumoniae pneumonia should not be underestimated.

S. pneumoniae possesses a series of virulence factors that are required for immune evasion against the host complement system[3]. Capsular polysaccharide is one of the most vital virulence factors in S. pneumoniae, which currently around 100 different capsular serotypes, divided into 25 individual types and 21 serogroups, each composed of several serotypes with related capsulate determinants, have been identified[4]. Pneumococcal conjugate vaccines (PCVs) can enormously decline the incidence of pneumococcal diseases in terms of targeting certain specific capsular polysaccharide of S. pneumoniae. For example, after the introduction of PCVs into the Turkish National Immunization Program in children, the potential coverage rate approximately decreased 20% for PCV7, PCV10, and PCV13 in 2008–2014, respectively[5]. The reduction of PCV coverage rates was attributed to increasing vaccinated children.

Although it was proved to be effective by vaccination in preventing pneumococcal diseases, the prevalence of CAP remained up to 65.80% in Chinese children under 5 years old because of a low rate of vaccination[6]. In addition, currently available PCVs cannot well dispose of non-PCV serotypes or non-capsular serotypes in isolates[7]. Therefore, many researchers have turned their attention from PCVs to other virulence factors. Well-known virulence factors naturally existed in S. pneumoniae, including autolysin (lytA), pneumolysin (ply), pneumococcal surface adhesin A (psaA), pneumococcal surface protein A (pspA)[8]. These factors seemed to be highly conserved and widely distributed in human pneumococci, implying the virulence target might be a fantastic vaccine candidate in the future.

Even though currently many studies on the molecular epidemiology of S. pneumoniae have been published in China (especially serotype and resistance), little is focused on the underlying relationship among STs, serotypes and virulence factors. So we aimed to investigate the molecular epidemiology of S. pneumoniae isolated from children under 5 years with CAP in Chengdu, clarify the relationship between STs, serotypes and virulence factors in order to provide laboratory-related data for pneumococcal infection prevention strategies and vaccines.

2. Materials And Methods

2.1 Study design

This retrospective study was conducted from January 2019 to January 2021 at Chengdu Women's and Children's Central Hospital, affiliated University of Electronic Science and Technology of China, which has about 2.6 million outpatient visits each year and 1840 establishment beds. The suspected diagnose of CAP was based on the following symptoms: the onset of systemic infection such as elevated temperature (> 38.0℃) or chills, and suffering from acute lower respiratory tract infection such as cough, gasp, accelerated breathing and abnormal auscultation, which was confirmed by X-ray or computerized tomographic scan of the chest. The eligibility criteria for enrolled participants included the following: (1) diagnosed as CAP (fulfilled above conditions); (2) younger than 5 years old; (3) specimens from patients like sputum, blood or alveolar lavage fluid isolated S. pneumoniae strains; (4) not vaccinated any PCVs. Excluding multiple admissions for the same patient and strains consecutively isolated from the same patient. In addition, hospitalized patient information was retrospectively extracted from the clinical records system, including elementary demographics, clinical manifestation, examination results.

2.2 Specimen collection, culture and identification

Specimens were collected by professional physicians or nurses in accordance with standard operating procedures and then transported as soon as possible to clinical microbiology laboratory, generally which were streaked on 5% sheep blood Columbia agar plate and cultivated overnight at 37°C in the 5% CO2 incubator. Blood specimen was cultured using the BacT/Alert 3D system (BD, US), and then the positive alert blood specimen was transferred and cultured to 5% sheep blood plate. Preliminary identification was performed on basis of morphology and Gram staining of the isolates, subsequently confirmed by VITEK 2 Compact (BioMérieux, France).

2.3 DNA preparation

Firstly sterilized normal saline was used to adjust the suspension to 1 McFarland turbidity. The bacterial suspension was subsequently centrifuged at 10000 rpm for one minute at 4℃, and the supernatant was discarded. Ultimately DNA was extracted from S. pneumoniae strains using TIANamp Bacteria DNA Kits (TIANGEN, China) according to the manufacturer’s instructions.

2.4 Serotyping

Serotyping was determined by multiplex polymerase chain reaction (PCR) as described in the previous article[9]. A total of nineteen primers were assigned to five multiplex reaction systems. A specific primer pair to the cpsA was used as the internal positive control. Serotype 6 was divided into 6A, 6B and 6C using the method referred to the previous study[10]. The serotype of isolates belonged to untyped when the serotype was not detected using the method mentioned above. Amplified PCR products were Semi-quantitatively analyzed by agarose gel electrophoresis assay, positive bands sequencing was subsequently entrusted to HuaDa Gene Biological Technology Company, Limited, HuaDa, China. Afterward, the coverage rates of PCV were calculated by the proportion of isolates classified serotypes.

2.5 Multilocus Sequence typing

Multilocus sequence typing (MLST) was performed by PCR amplification, purification, and sequencing of seven housekeeping genes (aroE, gdh, gki, recP, spi, xpt and ddl) according to PubMLST database (https://pubmlst.org/spneumoniae/)[11]. Alleles and sequence types (STs) were determined by sequence alignment using MEGA 7.0 software (https://www.megasoftware.net/), next querying MLST database mentioned above. Alleles or STs that were distinct from the current database were submitted to the curator of the database in order to achieve a new number assignment. Clonal complex (CC) was defined as a group of STs that shared five or more of seven alleles. Furthermore, the phylogenetic tree based on spliced genetic sequences was constructed by the UPGMA method with 1000-fold bootstrap resamplings and Tamura-Nei substitution model using MEGA 7.0 software. Interactive Tree Of Life (iTOL) v5 (https://itol.embl.de/) was used to display, annotate and embellish the phylogenetic tree[12]. Eventually, we used the goeBURST algorithm that a refinement of the eBURST algorithm by Feil et al[13] implemented in PHYLOViZ 2.0 software (http://www.phyloviz.net/) to display the relationship between the STs and serotypes.

2.6 Virulence genes detection

A total of seven virulence genes were detected using PCR method as previous study described[14], including autolysin (lytA), pneumolysin (ply), pneumococcal surface adhesin A (psaA), pneumococcal surface protein A (pspA), neuraminidase (nanA), pneumococcal adherence and virulence factor A (pavA), ion transporters (piaA).

2.7 Statistical analysis

All statistics in this study were performed using SPSS 22.0 software (https://www.ibm.com/analytics/spss-statistics-software). Continuous variables were represented by medians and quartile intervals, and categorical variables were represented by percentage. Categorical data were analyzed using a Chi-square test or Fisher’s exact test for difference. Correspondence analysis was used to further determine the relationships between STs and serotypes. Two-sided P-values < 0.05 were considered statistically significant.

3. Results

3.1 Demographic and clinical characteristics

A total of 108 children suffering from S. pneumoniae CAP were included in this study. Demographic and clinical characteristics from participants were as shown in Table 1. The ratio of male to female was 1.77 (69/39) and their age ranged from 0 to 5 years old, with median age (P25, P75) being 1 (0.75, 2) year. The number of groups based on age was 3 groups with 38 (35.19%) aged < 1 year, 39 (36.11%) aged 1–2 years and 31 (28.70%) aged 2–5 years. The majority of specimen types isolated from the patients was sputum accounting for 97.22% and the rest of specimen types included alveolar lavage fluid (2, 1.85%) and blood (1, 0.93%).

3.2 Serotypes distribution and vaccine coverage

There were 15 serotypes detected among 100 S. pneumoniae isolates in this study, while 8 of all strains were untyped by PCR method. The prevalent serotype was 19F accounting for 18.52%, followed by 6B (17.59%), 19A (13.89%), 6A (6.48%), 23F (5.56%). Classified by age, the predominant serotypes were 19F (21.05%), 6B (20.51%) and 19F (22.58%) among the isolates from CAP patients aged < 1 year, aged 1–2 years and aged 2–5 years, respectively.

The overall coverage rates of PCV7, 10, 13 among all isolates were 47.32%, 48.1%, 75%, respectively (Fig. 1). The coverage rates of PCV7 among the isolates from CAP patients aged < 1 year, aged 1–2 years and aged 2–5 years were 44.7%, 43.6%, 54.8%, respectively, while the coverage rates of PCV13 among those were 84.2%, 66.7%, 74.2%, respectively. Though the coverage rates of each PCV among all isolates between different groups sorted by age had no significant difference, generally, the coverage rate of PCV13 elicited significantly higher than that of PCV7 or PCV10.

3.3 MLST

A total of 56 different STs were identified among 108 S. pneumoniae isolates, of which there were 6 novel alleles (AroE552-553, gdh667, xpt937, xpt985, ddl1015) and 25 novel STs (ST15235, ST15272, ST15665, ST15667-15668, ST15672-15673, ST16093-16094, ST16445-16456, ST16467-16470) found in this study. The dominant STs were ST271 (22.22%) and ST320 (12.04%). As shown in Fig. 2, According to the substitution of housekeeping genes nucleotide among these STs, 56 STs were divided into 9 groups defined by bootstrap values of ≥ 70%. Every color strip indicated a CC and the grey strip represented singletons. The prevailing CC was CC271 (41.67%) with 8 different STs, including ST271, ST320, ST236, ST1937, ST16212, ST16093, ST16452 and ST16468. The correlation analysis between STs and serotypes showed that CC271 was closely associated with serotypes 19F and 19A, especially ST271 19F and ST320 19A, accounting for about 25% of all strains (Fig. 3).

3.4 Virulence genes distribution and pattern

Virulence genes distribution was listed in Table 2. Generally, all isolates carried lytA and most of the isolates carried ply (99.07%), psaA (98.15%), nanA (88.89%), pavA (96.30%) and piaA (97.22%). Approximately 45.37% of the isolates harbored pspA. It was worth noting that parts of virulence factors were significantly related to CC. In brief, the carriage rate of nanA in the CC271 group was higher than that in the non-CC271 group (P = 0.013), inversely the carriage rate of pspA in the CC271 group was lower than that in the non-CC271 group (P < 0.001).

There were 11 virulence patterns found in this study. The most common patterns were lytA-ply-psaA-nanA-pavA-piaA profile among all isolates accounting for 46.30%, with the relevant serotype being 19F and 19A and the lytA-ply-psaA-pspA-nanA-pavA-piaA profile among the isolates accounting for 35.19%, with correlated serotype being 6B (Table 3).

Table 1

Demographic and clinical characteristics from participants. 

Characteristics

No. of patients

%

Total

108

100

Gender

    

Male

69

63.89

female

39

36.11

Age (years)

    

< 1

38

35.19

1–2

39

36.11

2–5

31

28.70

Peoples

    

Han

103

95.37

Minority

5

4.63

Specimen types from

    

sputum

105

97.22

alveolar lavage fluid

2

1.85

blood

1

0.93

Clinical manifestations

    

Fever

76

70.37

Cough

103

95.37

Moist rales

58

53.70

Complications

    

Sepsis

14

12.96

Otitis media

5

4.63

Diarrhea

4

3.70

Endobronchitis

2

1.85

Respiratory failure

2

1.85

Length of stay (days; median and IQR)

7 (6–8)

-

Antibiotic use before admission

    

Yes

79

73.15

No

29

26.85

Prognosis

    

Cured

95

87.96

Improved

13

12.04

Dead

0

0


Table 2

Virulence genes distribution of S. pneumoniaeisolated from children with CAP. 

Virulence genes

Total n(%) (n = 108)

CC271 n(%) (n = 45)

Non-CC271 n(%) (n = 63)

χ2

P
 

positive

negative

positive

negative

positive

negative

    

lytA

108 (100)

0 (0)

45 (100)

0 (0)

63 (100)

0 (0)

-

-

ply

107 (99.07)

1 (0.93)

44 (97.78)

1 (2.22)

63 (100)

0 (0)

1.413

0.235

psaA

106 (98.15)

2 (1.85)

45 (100)

0 (0)

61 (96.83)

2 (3.17)

1.456

0.228

pspA

49 (45.37)

59 (54.63)

1 (2.22)

44 (97.78)

48 (76.19)

15 (23.81)

57.945

< 0.001

nanA

96 (88.89)

12 (11.11)

44 (97.78)

1 (2.22)

52 (82.54)

11 (17.46)

6.171

0.013

pavA

104 (96.30)

4 (3.70)

42 (93.33)

3 (6.67)

62 (98.41)

1 (1.59)

1.899

0.168

piaA

105 (97.22)

3 (2.78)

42 (93.33)

3 (6.67)

63 (100)

0 (0)

2.204

0.138
Bold indicated statistical significance among different CC groups (P < 0.05)
CAP: community acquired pneumonia


Table 3

The relationship between virulence pattern and serotypes of S. pneumoniae isolated from children with CAP. 

Virulence pattern

Isolates (No.)

Proportion (%)

Related serotypes (No.)

lytA-ply-psaA-pspA-nanA-pavA-piaA

38

35.19

6B (16), 3 (6), 23F (1), 15A (2), 6C (3), 7 (1), 23A(2), 15B (1), 18 (2), untyped (4)

lytA-ply-psaA-pspA-nanA-piaA

1

0.93

19F (1)

lytA-ply-psaA-nanA-piaA

2

1.85

19F (1), 23F (1)

lytA-ply-psaA-piaA

1

0.93

1 (1)

lytA-ply-psaA-nanA-pavA

2

1.85

19F (1), 19A (1)

lytA-psaA -nanA-pavA

1

0.93

19F (1)

lytA-ply-psaA-pspA-pavA-piaA

9

8.33

14 (1), 6C (2), 23F (3), 6B (2), untyped (1)

lytA-ply-psaA-pavA-piaA

2

1.85

14 (1), 7 (1)

lytA-ply-nanA-pavA-piaA

1

0.93

23F (1)

lytA-ply-pspA-nanA-pavA-piaA

1

0.93

untyped (1)

lytA-ply-psaA-nanA-pavA-piaA

50

46.30

19F (16), 19A (14), 15 A(1), 15B (3), 18B/C (1), 23A (2), 6A (7), 6B (1), 6C (1), 9V (1),3 (1), untyped (2)
CAP: community acquired pneumonia


4 Discussion

S. Pneumoniae is one of the common pathogens in CAP that leading cause of mortality in children under 5 years old. To prevent pneumococcal-related diseases and reduce the economic burden of that, many developed countries have introduced PCVs targeting serotypes of S. Pneumoniae into National Immunization Program. However, PCVs targeting capsular polysaccharide locus cannot overcome replacement serotype or non-capsular serotypes among S. Pneumoniae. Virulence factors, as good antigen targets, can effectively supplement the failure of PCVs. Therefore, we conducted this retrospective study to investigate molecular characteristics of STs, serotypes and virulence factors of S. Pneumoniae isolated from children under 5 years with CAP in Chengdu, so as to grasp epidemiological surveillance data and guide for vaccine vaccination and development.

Large, accurate, long-term serotype distribution monitoring of S. Pneumoniae isolates contributes to better multivalent vaccine development and guides for vaccination. In this study, The prevalent serotype was 19F, followed by 6B, 19A, 6A, 23F, which resembled other cities in China or Asian countries, but their rank order varied by region[15–17]. However, the prevailing serotypes were 19A, 7F, 3 and 1 in Europe and that was 14, 1 and 5 in America[18, 19]. Exception for region factor, S. Pneumoniae serotype distribution was limited by vaccine usage, age and so on. For example, the prevalent serotypes were 19F, 4, 6, 18 in the pre-PCV7 period in Russian children, and that serotypes were 4, 6, 18 after the introduction of PCV7[20], indicating a substantial reduction of 19F by vaccination. It also can be seen that different multivalent vaccines can only prevent infection with certain serotypes of S. pneumoniae. In this study, The overall coverage rates of PCV7, 10 among all isolates were 47.32%, 48.1%, respectively. PCV10 seemed not better than PCV7 because 1, 5, 7F were rarely detected in our study. When taking considering into age factor, we found that the coverage rate of PCV13 among isolates from children aged < 1 year was as high as 84.2%, which surpassed that of PCV7 and PCV10 among isolates from all children with no vaccination. We strongly suggest that it would be best to vaccinate with PCV13 for infants in the high coverage areas of serotypes, which can release the pressure caused by most vaccine-covered pneumococcal diseases.

Compared to pulsed-field gel electrophoresis (PFGE) technology, MLST is a simple, effective, comparable method to describe the molecular characterization of microorganisms. It has been applied in many aspects, such as epidemiological investigation, molecular evolution of bacteria[21]. In this study, CC271 was the most common group of STs, including the prevalent ST271 and ST320. A few previous research have suggested that there is a good correlation between STs and serotypes. For example, ST320 19A was reported as the prevalent clone in Latin American[22], while ST271 19F was prevalent in China[23]. By using minimum spanning tree analysis, we found that ST271 19F and ST320 19A originated from Taiwan19F-14 were common clones in this study. Additionally, there may exist also consistency between some STs and serotypes due to the limited sample size in this study, such as ST505 with serotype 3, ST3173 with serotype 6B, ST876 with serotype 14, which was similar to the previous studies reported in China[8, 24]. Interestingly, one ST may contain diverse serotypes sometimes. In this study, ST271 included 6A and 1 besides the absolutely prevailing 19F. The results by Tocheva et al[25] also showed that ST199 contained serotypes 19A, 14, 15 and ST162 contained 19F, 9V. Taken together, these results of studies indicated there was the evolution of serotypes switching among STs in order to adapt vaccine or antimicrobial pressure[26].

It is established that virulence factors, as important components of bacteria, can be alternative vaccine antigens. Recent report proved that conserved proteins that have been proposed for candidate vaccines should possess favorable immunogenicity and immunoreactivity of antigen, such as pspA, inhibiting the activation of the C3 complement protein and elevating the colonization ability of S. Pneumoniae [27]. Besides the above-mentioned conditions, an additional requirement for selecting virulence factors as candidate vaccines is widespread distribution in the population. In this study, all isolates carried lytA and most of the isolates carried ply, psaA, nanA, pavA and piaA, which was similar to other cities of China[14]. Based on our results, we suggested that these virulence factors including lytA, ply, pasA, nanA, pavA, piaA might be candidate vaccines in the future. However, the carrier rate of pspA was significantly lower than that of other virulence factors, which might be not suitable for vaccine candidate. Currently, little is known about the relationship between ST/CC and virulence factor. A few previous studies reported that CC271 carried more virulence factors than non-CC271[28, 29]. In this study, the positive rate of nanA in the CC271 was higher than that in the non-CC271. Recent research has revealed that pneumococcal neuraminidase A encoded by nanA gene can synergize with that of influenza A virus strengthening cytotoxicity[30]. These results together indicated CC271 harboring more nanA might exhibit much pathogenicity. Additionally, the most common virulence pattern was lytA-ply-psaA-nanA-pavA-piaA mainly distributed in serotype 19F and 19A in this study, while lytA-ply-psaA-pspA-nanA-pavA-piaA pattern in serotype 6B. It can be seen that the distribution of virulence factors might be associated with serotypes. Virulence factor pspA consisted of two families (family 1 and family 2), whose distribution was different from region. It was reported that the carriage rate of family 2 was higher than that of family 1 in Japan[31]. Although pspA of 19F and 6B belong to family 1, their distribution is also distinct. Our results indicated that serotype 6B isolates were more aggressive than 19F isolates. Furthermore, purified protein-vaccine from the homologous family 2 cannot against infection caused by S. pneumoniae carrying family 1[32, 33]. Therefore, we considered that supplemental vaccine strategy that PCVs coupled with certain virulence factors maybe reduce the burden of pneumococcal disease based on the distribution of STs, serotype and virulence factors[34].

To the best of our knowledge, this is the first study on potential relationships between molecular characteristic of S. pneumoniae from children with CAP under 5 years in Chengdu via minimum spanning tree and correspondence analysis, while several limitations warrant attention in this study. Firstly, the sample size of pneumococcus needs to be further expanded to a more comprehensive assessment of molecular epidemiology of S. pneumoniae. Secondly, some serotypes have not been identified so far, which might need other methods to identify. Ultimately, The results only represent the data of the last two years when the time span of specimens collection was between January 2019 to January 2021 in this article. More high-quality and long-term surveillance on STs, serotypes, virulence factors needs to be carried out in the future.

In summary, our findings indicated that CC271 was the predominant epidemic clone, especially ST271 19F and ST320 19A. Moreover, there were strong relationships between STs, serotypes, virulence factors, which indicated the formulation of vaccine strategies needs to carefully consider molecular epidemiological factors.

Declarations

Ethics statement

This study was approved by Chengdu Women's and Children's Central Hospital Ethics Committee, School of Medicine, University of Electronic Science and Technology of China (Approval No. 202035), confirming that all methods were performed in accordance with relevant guidelines and regulations. Meanwhile, informed consent was obtained from their legal guardians.

Consent for publication 

Not applicable

Availability of data and materials

All datasets generated for this study are included in the manuscript.

Competing interests

The authors declare no conflict of interest.

Funding

Not applicable

Author contributions

HJC conceived and designed the experiments. All authors contributed reagents/materials/analysis tools. HJC wrote the first draft of the manuscript and CGL reviewed the paper. All authors together agree to submit the manuscript.

Acknowledgements

We thank the patients who provided clinical records and specimens, and professor Zhao from Chengdu Medical College for his help on experiment.

References

  1. World Health Organization. Fact Sheet-Pneumonia.2019. https://www.who.int/en/news-room/fact-sheets/detail/pneumonia. Accessed 8 November 2021.
  2. Ning G, Wang X, Wu D, Yin Z, Li Y, Wang H, et al. The etiology of community-acquired pneumonia among children under 5years of age in mainland China, 2001-2015: A systematic review. Hum Vaccin Immunother 2017;13:2742-2750.
  3. Kadioglu A, Weiser J N, Paton J C, Andrew P W. The role of Streptococcus pneumoniae virulence factors in host respiratory colonization and disease. Nat Rev Microbiol 2008;6:288-301.
  4. Campanero-Rhodes M A, Lacoma A, Prat C, García E, Solís D. Development and Evaluation of a Microarray Platform for Detection of Serum Antibodies Against Streptococcus pneumoniae Capsular Polysaccharides. Anal Chem 2020;92:7437-7443.
  5. Ceyhan M, Ozsurekci Y, Gürler N, Öksüz L, Aydemir S, Ozkan S, et al. Serotype distribution of Streptococcus pneumoniae in children with invasive diseases in Turkey: 2008-2014. Hum Vaccin Immunother 2016;12:308-313.
  6. Sun Y, Li H, Pei Z, Wang S, Feng J, Xu L, et al. Incidence of community-acquired pneumonia in urban China: A national population-based study. Vaccine 2020;38:8362-8370.
  7. Keller L E, Robinson D A, McDaniel L S. Nonencapsulated Streptococcus pneumoniae: Emergence and Pathogenesis. MBio 2016;7:e01792.
  8. Li L, Fu J, Li S, Guo D, Chen Z, Chen S, et al. Phenotypic and molecular characterization of Streptococcus pneumoniae in pre-conjugate vaccine era: A Chinese hospital-based retrospective study. Vaccine 2018;36:599-605.
  9. Li X X, Xiao S Z, Gu F F, Zhao S Y, Xie Q, Sheng Z K, et al. Serotype Distribution, Antimicrobial Susceptibility, and Multilocus Sequencing Type (MLST) of Streptococcus pneumoniae From Adults of Three Hospitals in Shanghai, China. Front Cell Infect Microbiol 2019;9:407.
  10. Jin P, Xiao M, Kong F, Oftadeh S, Zhou F, Liu C, et al. Simple, accurate, serotype-specific PCR assay to differentiate Streptococcus pneumoniae serotypes 6A, 6B, and 6C. J Clin Microbiol 2009;47:2470-2474.
  11. Jolley K A, Bray J E, Maiden M C J. Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications. Wellcome Open Res 2018;3:124.
  12. Letunic I, Bork P. Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res 2021;49:W293-W296.
  13. Feil E J, Li B C, Aanensen D M, Hanage W P, Spratt B G. eBURST: inferring patterns of evolutionary descent among clusters of related bacterial genotypes from multilocus sequence typing data. J Bacteriol 2004;186:1518-1530.
  14. Zhao W, Pan F, Wang B, Wang C, Sun Y, Zhang T, et al. Epidemiology Characteristics of Streptococcus pneumoniae From Children With Pneumonia in Shanghai: A Retrospective Study. Front Cell Infect Microbiol 2019;9:258.
  15. Wang Q, Shi W, Li Y, Gao W, Yuan L, Dong F, et al. Serotype distribution of Streptococcus pneumoniae isolated from children hospitalized in Beijing children's hospital (2013-2019). Vaccine 2020;38:7858-7864.
  16. Lyu S, Hu H L, Yang Y H, Yao K H. A systematic review about Streptococcus Pneumoniae serotype distribution in children in mainland of China before the PCV13 was licensed. Expert Rev Vaccines 2017;16:997-1006.
  17. Tai S S. Streptococcus pneumoniae Serotype Distribution and Pneumococcal Conjugate Vaccine Serotype Coverage among Pediatric Patients in East and Southeast Asia, 2000-2014: a Pooled Data Analysis. Vaccines (Basel) 2016;4:4.
  18. Tin Tin Htar M, Christopoulou D, Schmitt H J. Pneumococcal serotype evolution in Western Europe. BMC Infect Dis 2015;15:419.
  19. Gentile A, Bardach A, Ciapponi A, Garcia-Marti S, Aruj P, Glujovsky D, et al. Epidemiology of community-acquired pneumonia in children of Latin America and the Caribbean: a systematic review and meta-analysis. Int J Infect Dis 2012;16:e5-15.
  20. Sidorenko S, Rennert W, Lobzin Y, Briko N, Kozlov R, Namazova-Baranova L, et al. Multicenter study of serotype distribution of Streptococcus pneumoniae nasopharyngeal isolates from healthy children in the Russian Federation after introduction of PCV13 into the National Vaccination Calendar. Diagn Microbiol Infect Dis 2020;96:114914.
  21. Maiden M C. Multilocus sequence typing of bacteria. Annu Rev Microbiol 2006;60:561-588.
  22. Moreno J, Duarte C, Cassiolato A P, Chacón G C, Alarcon P, Sánchez J, et al. Molecular characterization of Latin American invasive Streptococcus pneumoniae serotype 19A isolates. Vaccine 2020;38:3524-3530.
  23. Wang X, Cong Z, Huang W, Li C. Molecular characterization of Streptococcus pneumoniae isolated from pediatric patients in Shanghai, China. Pediatr Pulmonol 2020;55:2135-2141.
  24. Yang Y, Hua C Z, Fang C, Xie Y P, Li W, Fu Y, et al. Properties of Mucoid Serotype 3 Streptococcus pneumoniae From Children in China. Front Cell Infect Microbiol 2021;11:648040.
  25. Tocheva A S, Jefferies J M, Christodoulides M, Faust S N, Clarke S C. Distribution of carried pneumococcal clones in UK children following the introduction of the 7-valent pneumococcal conjugate vaccine: a 3-year cross-sectional population based analysis. Vaccine 2013;31:3187-3190.
  26. Croucher N J, Kagedan L, Thompson C M, Parkhill J, Bentley S D, Finkelstein J A, et al. Selective and genetic constraints on pneumococcal serotype switching. PLoS Genet 2015;11:e1005095.
  27. Morais V, Texeira E, Suarez N. Next-Generation Whole-Cell Pneumococcal Vaccine. Vaccines (Basel) 2019;7:151.
  28. Liang Z, Fu J, Li L, Yi R, Xu S, Chen J, et al. Molecular epidemiology of Streptococcus pneumoniae isolated from pediatric community-acquired pneumonia in pre-conjugate vaccine era in Western China. Ann Clin Microbiol Antimicrob 2021;20:4.
  29. Fu J, Li L, Liang Z, Xu S, Lin N, Qin P, et al. Etiology of acute otitis media and phenotypic-molecular characterization of Streptococcus pneumoniae isolated from children in Liuzhou, China. BMC Infect Dis 2019;19:168.
  30. Wren J T, Blevins L K, Pang B, Basu Roy A, Oliver M B, Reimche J L, et al. Pneumococcal Neuraminidase A (NanA) Promotes Biofilm Formation and Synergizes with Influenza A Virus in Nasal Colonization and Middle Ear Infection. Infect Immun 2017;85: e01044-16.
  31. Hotomi M, Togawa A, Kono M, Ikeda Y, Takei S, Hollingshead S K, et al. PspA family distribution, antimicrobial resistance and serotype of Streptococcus pneumoniae isolated from upper respiratory tract infections in Japan. PLoS One 2013;8:e58124.
  32. Rodrigues T C, Oliveira M L S, Soares-Schanoski A, Chavez-Rico S L, Figueiredo D B, Gonçalves V M, et al. Mucosal immunization with PspA (Pneumococcal surface protein A)-adsorbed nanoparticles targeting the lungs for protection against pneumococcal infection. PLoS One 2018;13:e0191692.
  33. Lin H, Peng Y, Lin Z, Zhang S, Guo Y. Development of a conjugate vaccine against invasive pneumococcal disease based on capsular polysaccharides coupled with PspA/family 1 protein of Streptococcus pneumoniae. Microb Pathog 2015;83-84:35-40.
  34. Colichio G B, Oliveira G S, Rodrigues T C, Oliveira M L S, Miyaji E N. Efficacy of a Protein Vaccine and a Conjugate Vaccine Against Co-colonization with Vaccine-type and Non-vaccine Type Pneumococci in Mice. Pathogens 2020;9:278.