Introduction
Pneumococcal conjugated vaccines (PCV) represent a global strategy to reduce the burden of diseases caused by S. pneumoniae (pneumococcus). However, with their widespread implementation, concerns have arisen about the potential increase in S. aureus infections [Reference Ebruke1]. S. aureus is an opportunistic organism that colonizes mucous membranes and the skin of 20%–80% of the human population without symptoms [Reference Brown, Leech, Rogers and McLoughlin2]. It is a common cause of bacterial infections both in the community and hospital settings. In 2024, the World Health Organization (WHO) updated the list of the main drug-resistant pathogens, including S. aureus methicillin-resistant and S. pneumoniae macrolide-resistant, as bacteria with high and medium priority, respectively [3].
Following the introduction of PCVs, there has been a replacement of pneumococcal serotypes with an increase in the prevalence of non-vaccine serotypes (not included in the PCVs) that are clinically relevant, such as serotype 19A. The first clinical trial evaluating PCV7 and its impact on S. aureus colonization reported an increase in the frequency of colonization after receiving two doses [Reference van Gils4]. Furthermore, an inverse association between S. pneumoniae vaccine serotypes (serotypes included in the PCVs, PCV-serotypes) and S. aureus is also reported [Reference van Gils4–Reference Quintero6], suggesting that colonization dynamics involve a complex interplay of bacterial and host factors.
Describing the patterns of antibiotic resistance in community-acquired S. aureus is important because antibiotic treatment with a drug to which S. aureus is resistant promotes colonization and, consequently, the development of potentially serious infections [Reference Tong7]. Additionally, there is no information on the rate of methicillin-resistant S. aureus (MRSA) colonizing healthy Peruvian children. The study aimed to assess the antibiotic resistance of S. aureus in nasopharyngeal carriers under 24 months of age between 2018 and 2019 in Lima, Peru, determine the prevalence of carriage and co-colonization with S. pneumoniae, and identify potential factors associated with colonization status.
Methods
Study design and population
We used samples from a cross-sectional multicentre study conducted to evaluate changes in pneumococcal serotypes after PCV13 introduction [Reference Gonzales8]. The previous study enrolled 1,000 healthy children between 2 and 24 months of age who attended the well-child outpatient or the vaccination services of five hospitals in Lima, Peru. Only 894 children were enrolled from the previous study since their parents or caregivers provided consent to use their data and samples for future studies.
The main study defined a ‘healthy child’ as any child who attended the growth, development, and/or immunization outpatient clinic at Cayetano Heredia Hospital, Instituto de Salud del Niño Sede Breña, Edgardo Rebagliati Martins Hospital, Hospital Docente Madre-Niño San Bartolome, or Daniel Alcides Carrion Hospital, and whose parent or caregiver reported that, at the time of enrolment, the child did not suffer from any major illness such as pneumonia, sepsis, bacteraemia, or meningitis/encephalitis. Children might have had mild respiratory infections such as rhinorrhea, mild cough, sneezing, or a temperature of 38.0 to 38.5 °C at enrolment [Reference Gonzales8].
Sample collection
Nasopharyngeal swab samples were taken and placed in STGG transport medium (skim milk, tryptic soy broth, glucose, glycerol) and transported in coolers within 8 h to the Paediatric Infectious Diseases Laboratory (LIP) at Cayetano Heredia University (the central laboratory of the main study), for subsequent microbiological analysis.
Laboratory studies
As part of the main study, the swabs were placed in 3 ml of THY medium (Todd-Hewitt Broth with 0.5% yeast extract and 5% rabbit serum) and incubated for 4 to 6 h at 37 °C with 5.0% CO2. The bacterial growth was cultured on blood agar (Trypto-Casein Soy Agar enriched with 5% sheep blood) under the same previously described conditions for 24 h. S. pneumoniae isolates were identified through colony morphology, the presence of alpha haemolysis, Gram staining, bile solubility, and optochin sensitivity [Reference Gonzales9]. S. aureus isolates were identified through the bacterial growth obtained in THY medium, plated on Mannitol Salt Agar, and with the Gram staining, catalase, and coagulase tests [Reference Chien10]. The ATCC 25923 and ATCC 43300 strains were used as controls for S. aureus and D39/NCTC 7466 for S. pneumoniae. Pneumococcal serotype was determined by whole genome sequencing (WGS).
Susceptibility testing
Antibiotic susceptibility of S. aureus to erythromycin, clindamycin, cefoxitin (to determine methicillin resistance) [11], penicillin, and trimethoprim-sulfamethoxazole (SXT) were determined by Kirby Bauer method using antibiotic disks (Oxoid Ltd., Basingstoke, Hans, UK), as well as D-Test to identify clindamycin resistance induced by erythromycin on Mueller Hinton agar and incubated for 24 h at 37 °C [11]. Antibiotic susceptibility of S. pneumoniae to azithromycin, clindamycin, ceftriaxone, penicillin, and SXT was determined by MIC using E-test® (AB Biodisk, Solna, Sweden) on Mueller Hinton agar with 5% sheep blood and incubated for 24 h at 37 °C in an atmosphere of 5% CO2. Interpretation was carried out according to the Clinical and Laboratory Standards Institute (CLSI) guidelines [11].
Statistical analysis
Bivariate analyses were performed to determine the associations between the dependent variable (colonization status), the independent variable, and covariates of interest using parametric and non-parametric tests such as χ2 and Fisher’s exact test (variables that meet or do not meet Cochran’s rules, respectively), ANOVA (variables with normality distribution), and Kruskal–Wallis (variables without normality distribution). Multinomial regressions were conducted to determine associations with the colonization status of S. pneumoniae, S. aureus, and co-colonized subjects. Additionally, Poisson regression with robust variances was performed to explore the association between children colonized by S. pneumoniae and S. aureus (this model only considered these two variables independently and not as a polytomous variable as used in the multinomial regression model). In the multivariate analysis, covariates were included as identified in the literature and by the clinical experience of previous studies that resulted in statistical significance in our bivariate analyses, as identified in the Directed Acyclic Graph (DAG) (Supplementary Figure S1). The category ‘not colonized’ was considered the comparator group. All measures of association were the prevalence ratio (PR) (this was approximated using risk ratios (RR) and adjusted using the estimation option ‘rrr’ (multinomial model) and robust standard errors (Poisson model)) due to the study design. The analyses were performed in Stata/SE version 18.0 (StataCorp LP, College Station, TX), considering a 95% confidence interval, and the statistical significance level was set a p < 0.05.
Ethical aspects
This study was approved by the Institutional Review Board of Cayetano Heredia University (Lima, Peru)
Results
A total of 894 samples from healthy children enrolled between 2018 and 2019 were analyzed; 19.7% were positive only for S. aureus, 20.3% were positive only for S. pneumoniae, and 2.9% had co-colonization with both bacteria (Table 1). The mean age of the children was 11.1 months; 20.9% were preterm, and 37.2% completed the recommended 6 months of exclusive breastfeeding. Mild respiratory illness, characterized by mild cough, rhinorrhea, sneezing, or a temperature between 38°C and 38.5°C, was presented in 36.3%. PCV vaccination coverage, with at least one dose, was 89.4% (Table 1).
Table 1. Characteristics of children included in the study (N = 894)a

a Some variables may add <894 due to missing data.
b Median (IQR).
c Rhinorrhea, mild cough, sneezing, temperature ≤ 38.5 °C.
d Within 3 months preceding sampling.
e Season of swab collection.
f χ2 test.
g ANOVA.
h Kruskal–Wallis test.
i Fisher’s exact test.
Characteristics and factors associated with colonization
According to the multivariate analysis, the frequency of colonization by S. aureus alone was 0.59 times lower (95% CI: 0.37–0.94; p = 0.028) among those who completed exclusive breastfeeding and 0.48 times lower (95% CI: 0.26–0.87; p = 0.017) among those who received at least one dose of PCV13 (Table 2). The frequency of colonization by S. pneumoniae alone was 3.00 times higher (95% CI: 1.17–7.69; p = 0.022) among those who received at least one dose of PCV13, 1.95 times higher (CI: 1.32–2.90; p = 0.001) among those with mild respiratory illness, and 1.99 times higher (95% CI: 1.35–2.90; p = 0.001) among children living with other children. In contrast, the prevalence was 0.29 times lower (95% CI: 0.09–0.99; p = 0.048) among those who had been previously hospitalized. No significant association was found between co-colonization (S. pneumoniae + S. aureus) and PCV13 (p = 0.521), except for sex, where the prevalence of co-colonization prevalence was 2.50 times higher in males than in females (95% CI: 1.00–6.21; p = 0.049, (Table 2)).
Table 2. Factors associated with nasopharyngeal colonization in Peruvian children (multivariate analysis)

Note: Multivariate analysis using robust (modified) multinomial logistic regression. Non-colonized was the comparator group.
aPR: Prevalence ratio. 95% CI: 95% confidence interval. Ref: Reference category. -: Without observations.
Adjusting for sex, age (months), exclusive breastfeeding (6 months), day-care centre attendance, living with children <6 years old, mild respiratory illness, prior hospitalization, antibiotics use, and immunized with at least one PCV13 dose.
Additionally, children colonized by S. pneumoniae exhibited a reduced prevalence ratio (RP: 0.68, 95% CI: 0.46–0.99; p = 0.047) for colonization by S. aureus, suggesting an inverse association between the two potentially pathogenic species. A similar trend was observed with non-vaccine serotypes (RP: 0.69; 95% CI: 0.44–1.06; p = 0.087), although this finding did not reach statistical significance. In contrast, no association was identified with vaccine serotypes (RP: 0.70; 95% CI: 0.28–1.79; p = 0.461) (Supplementary Table S1). Information about pneumococcus serotypes of this series is described in a previous publication of our group [Reference Gonzales8].
Antibiotic susceptibility
We isolated S. aureus strains from the nasopharynx of 176 children (Table 1), with only three strains (1.7%) classified as MRSA based on the CLSI-recommended test of resistance to cefoxitin. In contrast, the majority of these strains (89.2%) were resistant to penicillin. Regarding macrolide resistance, 59 (33.5%) S. aureus strains were non-susceptible to erythromycin, while 36 strains (20.5%) were non-susceptible to clindamycin. The resistance to clindamycin was also evaluated through the D-test, which was positive in 34/59 erythromycin-resistant isolates (57.6%). Regarding this specific antibiotic, out of 36 strains non-susceptible to clindamycin, 34 showed inducible resistance, with 33 identified as MSSA (methicillin-sensitive S. aureus). Among 182 pneumococcal isolates evaluated, all were susceptible to penicillin using the non-meningeal breakpoint. Additionally, 48.9% of the isolates were resistant to the macrolide azithromycin (Table 3).
Table 3. Antibiotic resistance of S. aureus and Streptococcus pneumoniae strains

a S. aureus strains isolated from children colonized by only S. aures and co-colonized.
b S. pneumoniae strains isolated from children colonized by only S. pneumoniae and co-colonized.
c Non-meningitis breakpoint for S. pneumoniae strains.
d Erythromycin tested by S. aureus and azithromycin by S. pneumoniae.
e 34 isolates exhibited inducible resistance to clindamycin (D-test positive).
f SXT: trimethoprim-sulfamethoxazole.
Regarding the resistance pattern among S. aureus isolates, penicillin resistance was the most frequent (56.3%), while 14 strains (7.9%) were classified as pan-susceptibility. In S. pneumoniae, the most common resistance pattern involved azithromycin, clindamycin, and SXT (23.1%), followed by resistance to azithromycin and SXT (16.5%). Notably, 12.1% of pneumococcal isolates demonstrated pan-susceptibility to the five antibiotics evaluated (Table 4).
Table 4. Antibiotic resistance patterns S. aureus and Streptococcus pneumoniae strains isolated from healthy children

Abbreviations: FOX = cefoxitin, ERI/AZM = erythromycin tested by S. aureus or azithromycin by S. pneumoniae, PEN = penicillin, CLD = clindamycin, SXT = trimethoprim-sulfamethoxazole.
a 7 pneumococcal isolates were intermediate to azithromycin and/or trimethoprim-sulfamethoxazole.
b Pansensitive: Susceptible to cefoxitin, erythromycin, trimethoprim-sulfamethoxazole, penicillin and clindamycin.
Discussion
This study in Peruvian children observed nasopharyngeal colonization rates of 19.7% for S. aureus, 20.3% for S. pneumoniae of and 2.9% for co-colonization with both species. Importantly, despite concerns about antibiotic resistance, carriage strains of S. aureus and S. pneumoniae exhibited high susceptibility to methicillin and penicillin, respectively, in this population. While exclusive breastfeeding and PCV13 vaccination were associated with lower S. aureus colonization, PCV13 vaccination was associated with increased S. pneumoniae colonization, and an increased risk of colonization was observed in children experiencing mild respiratory illness and living with other children.
The observed S. aureus colonization rate of our study is lower than reported in other studies. This discrepancy may be attributed to the age range of our population, as S. aureus colonization rates typically exhibit age-related variations. Colonization rates are higher in infants (<6 months), decrease during early childhood, increase again between 6 and 12 years of age, and subsequently decline in adulthood [Reference Vieira12, Reference Sollid, Furberg, Hanssen and Johannessen13]. Previous studies have reported varying S. aureus colonization rates in children. For instance, Paraguay reported a higher colonization frequency (30.8%) in 2016 [Reference Benítez14], with similar rates in Brazil (31.1%), Colombia (33%), and Japan (28.2%). In contrast, Mexico (59.8%) and Venezuela (56%) showed higher rates [Reference Quintero6, Reference Hisata15–Reference Castro-Orozco18]. A study in children from Cajamarca, Peru, in 2009, reported a lower colonization rate of 12.7% using microbiological methods [Reference Chien10]. These variations in colonization rates across different regions can be attributed to a multitude of factors, including environmental and sociodemographic factors specific to each population [Reference Benítez14].
Exclusive breastfeeding and PCV13 vaccination were associated with decreased S. aureus colonization, aligning with the known benefits of breastfeeding on gut microbiota development [Reference Ho19–Reference Heisel21]. While some studies have suggested an association between PCV13 vaccination and decreased S. aureus colonization, this relationship is often attributed to the impact of the vaccine on specific pneumococcal serotypes [Reference van Gils4]. However, this mechanism was not observed in our study.
It is well-established that individuals colonized with MRSA have a higher risk of developing MRSA-related infections compared to those colonized with MSSA [Reference Kliegman22]. Risk factors for MRSA colonization include recent hospitalization, recent antibiotic treatment, and the use of beta-lactam antibiotics within the past 3 months [Reference Castro-Orozco18]. In our study, the prevalence of MRSA among S. aureus carriers was low (1.7%), consistent with the range reported in Latin America and globally (0.2–10%) [Reference Neyra23, Reference Moyo24]. In Peru, the reported MRSA prevalence was 0.9% in the general population and 0.3% in the military [Reference Neyra23, Reference Carmona, Sandoval and García25]. In Brazil, MRSA prevalence varied, with higher rates observed in day-care children (6.2%) compared to older than 5 years (1.2%) [Reference Lamaro-Cardoso17, Reference Braga26]. In contrast, significantly higher MRSA colonization rates have been observed in children older than 5 years in Tanzania (10.5%) and Ethiopia (9.8%) [Reference Moyo24, Reference Tigabu, Tiruneh and Mekonnen27].
Resistance to erythromycin and clindamycin among S. aureus isolates was observed at 33.5% and 20.5%, respectively, suggesting the presence of diverse macrolide resistance mechanisms, a trend consistent with observations in other Latin American countries [Reference Hamdan-Partida, Sainz-Espuñes and Bustos-Martínez16]. While clindamycin remains an option for treating S. aureus skin and soft tissue infections, its clinical utility is limited by emerging resistance, particularly inducible clindamycin resistance [Reference Khashei, Malekzadegan, Sedigh Ebrahim-Saraie and Razavi28, Reference DeLeo, Otto, Kreiswirth and Chambers29]. Notably, the D-test revealed a high prevalence of inducible clindamycin resistance (20.5%) in our study, emphasizing the importance of routine D-test performance in clinical microbiology [Reference Siberry, Tekle, Carroll and Dick30]. This finding raises concerns about the empirical use of clindamycin [Reference Moyo24, Reference Sedighi, Moez and Alikhani31]. In contrast, no resistance to trimethoprim-sulfamethoxazole was observed in our study, highlighting the continued importance of ongoing local antibiotic susceptibility surveillance to guide appropriate clinical management, especially considering the high resistance rates to trimethoprim-sulfamethoxazole (30%–60%) reported in Africa [Reference Moyo24].
The observed pneumococcal colonization frequency (20.2%) aligns with previously reported rates in low- and middle-income countries [Reference Gonzales8, Reference Adegbola32]. This frequency is slightly lower than observed in pre-PCV7 era studies in the region [Reference Chien10, Reference Mercado33], potentially reflecting the impact of PCV vaccination on reducing carriage of vaccine serotypes [Reference Gonzales8]. However, our study also found an association between PCV13 vaccination and increased pneumococcal colonization, primarily by non-PCV13 serotypes, suggesting potential selective pressure exerted by the vaccine [Reference Gonzales9]. Consistent with previous findings [Reference Mercado33–Reference Narváez35], factors associated with increased pneumococcal colonization included mild respiratory illness, living with other children, and younger age. Of note, hospitalization within 3 months prior to sample collection was associated with a decreased prevalence of pneumococcal colonization (aPR: 0.29; 95% CI: 0.09–0.99; p = 0.048), potentially due to recent antibiotic treatment in hospitalized children [Reference Gonzales9].
Macrolide resistance in S. pneumoniae isolates was observed at 48.9%, a significant increase compared to pre-PCV7 era findings in the same population [Reference Mercado33]. Molecular analysis from our previous studies [Reference Gonzales9] revealed that the macrolide resistance mechanisms were primarily mediated by the erm(B) gene (ribosomal methylation) and the MEGA element (mefA/E) (macrolide efflux pump).
The observed co-colonization frequency of 2.9% was comparable to that reported in a previous study among Peruvian children in Cajamarca (4.8%) using microbiological assays [Reference Chien10]. However, this rate is lower than observed in studies using molecular methods (28.2%), which generally exhibit higher sensitivity in detecting potential nasopharyngeal pathogens [Reference Chien10]. Care should be taken with prior antibiotic use, as DNA detection from non-viable bacteria may not indicate current colonization [Reference da Gloria S Carvalho36]. An inverse association was found between S. pneumoniae and S. aureus colonization. These two species, while both colonizing the nasopharynx, may exhibit competitive interactions [Reference van Gils4, Reference Bogaert5, Reference McNally37]. For instance, some S. pneumoniae strains can produce competence-stimulating peptides that inhibit S. aureus biofilm formation [Reference Khan38], while others produce H2O2, which can eradicate S. aureus in vitro [Reference Regev-Yochay39]. However, the impact of this interaction may vary depending on the age of the child, with competition potentially being less pronounced in older children [Reference Esposito40].
Our study had some limitations. Firstly, while S. pneumoniae strains were isolated promptly, S. aureus isolation was delayed (8–12 months). Although nasopharyngeal swabs were stored at −80 °C in STGG medium to maintain viability, this delay may have introduced potential biases. Secondly, the study relied on convenience sampling, potentially limiting the generalizability of the findings beyond the specific population of children attending public hospitals in Lima and Callao. While this population represents a significant proportion (27.2%) of children under 2 years old in Peru, it may not fully reflect the characteristics of children seeking care in private healthcare facilities or residing in other regions. Moreover, there is insufficient evidence to conclude an association with the co-colonization status, as it was an infrequent event, which may have affected the statistical power for this category; therefore, we recommend interpreting the results for this group as exploratory. Finally, limitations in immunization data collection, including incomplete records and the inability to assess adherence to the national immunization schedule, may have influenced the study findings.
Despite the introduction of three different PCVs into Peru’s National Immunization Program, the country lacks population-level surveillance of emerging bacterial pathogens, a crucial step following vaccine introduction, as exemplified by post-Hib vaccine surveillance [41]. This study highlights the need for continued surveillance of nasopharyngeal pathogens, such as S. pneumoniae and S. aureus, particularly given the increasing prevalence of antimicrobial resistance. Our findings suggest that promoting exclusive breastfeeding may contribute to a reduction in S. aureus colonization. Furthermore, while PCV13 vaccination showed an association with decreased S. aureus colonization, it was also associated with increased S. pneumoniae colonization, potentially due to selective pressure on non-vaccine serotypes. These findings emphasize the importance of ongoing surveillance to monitor the impact of PCVs on the nasopharyngeal microbiota and the emergence of potential vaccine-related consequences.
Supplementary material
The supplementary material for this article can be found at http://doi.org/10.1017/S0950268825100277.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Acknowledgements
We thank all members and study coordinators of Peruvian Pneumococcus Research Group for providing us with the samples and databases of the healthy carrier children study.
Author contribution
All authors initiated and designed this research study. Brayan Gonzales, Franco Castillo-Tokumori and Theresa Jean Ochoa were in charge of patient enrolment, sample collection and data acquisition. Valeria Mariana Li Valverde and Paulo Cesar Aguirre Castañeda performed the experiments. Valeria Mariana Li Valverde, Paulo Cesar Aguirre Castañeda and Brayan Gonzales analyzed and interpreted the results. Valeria Mariana Li Valverde, Paulo Cesar Aguirre Castañeda and Brayan Gonzales wrote the first draft of the manuscript. All authors contributed to the article and approved the submitted version.
Competing interests
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.
Ethical standard
The study involving human participants was reviewed and approved by Institutional Review Board of Cayetano Heredia University with unique code 206782.
Funding statement
The study was funded by internal funds of the Paediatric Infectious Diseases Laboratory at Cayetano Heredia University. JEV is in part supported by grants from the National Institutes of Health (NIH; 5R21AI144571-03, and 1R01AI175461-01A1). The content is solely the responsibility of the authors and does not necessarily represent the official view of the NIH.