The COVID-19 Pandemic as an Opportunity for Unravelling the Causative Association between Respiratory Viruses and Pneumococcus-Associated Disease in Young Children: A Prospective Study ======================================================================================================================================================================================== * Ron Dagan * Bart Adriaan van der Beek * Shalom Ben-Shimol * David Greenberg * Yonat Shemer-Avni * Daniel M. Weinberger * Dana Danino ## Summary **BACKGROUND** In young children, rates of lower respiratory infections (LRI) and invasive pneumococcal disease (IPD) have been associated with respiratory syncytial virus (RSV), human metapneumovirus (hMPV), influenza (flu), and parainfluenza (PIV) (collectively termed here as pneumonia and pneumococcal disease-associated viruses [PDA-viruses]). However, their contribution to the pathogenesis of these disease endpoints has not yet been elucidated. The COVID-19 pandemic provided a unique opportunity to examine the question. **METHODS** This prospective study comprised all children <5 years, living in southern Israel, during 2016 through 2021. The data were previously collected in multiple ongoing prospective surveillance programs and include: hospital visits for community-acquired alveolar pneumonia (CAAP), non-CAAP LRI; nasopharyngeal pneumococcal carriage (<3 years of age); respiratory virus activity; and nationwide, all-ages COVID-19 episodes and IPD in children <5 years. A hierarchical statistical model was developed to estimate the proportion of the different clinical endpoints attributable to each virus from monthly time series data, stratified by age and ethnicity. A separate model was fit for each endpoint, with covariates that included a linear time trend, 12-month harmonic variables to capture unexplained seasonal variations, and the proportion of tests positive for each virus in that month. **FINDINGS** During 2016 through 2021, 3,204, 26,695, 257, and 619 episodes of CAAP, non-CAAP LRI, pneumococcal bacteremic pneumonia and non-pneumonia IPD, respectively, were reported. Compared to 2016-2019, broad declines in the disease endpoints were observed shortly after the pandemic surge, coincident with a complete disappearance of all PDA-viruses and continued circulation of rhinovirus (RhV) and adenovirus (AdV). From April 2021, off-season and abrupt surges of all disease endpoints occurred, associated with similar dynamics among the PDA-viruses, which re-emerged sequentially. Using our model fit to the entire 2016-2021 period, 82% (95% CI, 75-88%) of CAAP episodes in 2021 were attributable to the common respiratory viruses, as were 22%-31% of the other disease endpoints. Virus-specific contributions to CAAP were: RSV, 49% (95% CI, 43-55%); hMPV, 13% (10-17%); PIV, 11% (7-15%); flu, 7% (1-13%). RhV and AdV did not contribute. RSV was the main contributor in all endpoints, especially in infants. Pneumococcal carriage prevalence remained largely stable throughout the study. **INTERPRETATION** RSV and hMPV play a critical role in the burden of CAAP and pneumococcal disease in children. Interventions targeting these viruses could have a secondary effect on the burden of disease typically attributed to bacteria. **FUNDING** There was no funding for this study. **Evidence before this study** Lower respiratory infections (LRI) and invasive pneumococcal disease (IPD) in young children, have often been associated with specific respiratory viruses, namely respiratory syncytial virus (RSV), human metapneumovirus (hMPV) influenza viruses (flu), and parainfluenza viruses (PIV) (termed in the current article pneumonia and pneumococcal disease-associated viruses [PDA-viruses]). However, their causative role as co-pathogens has not yet been fully elucidated. There is already ample evidence that bacteria and viruses interact to cause severe disease. This could be seen after the introduction of pneumococcal conjugate vaccines (PCVs), when there was a significant reduction in hospitalisation for viral lower respiratory infections (LRIs). This suggests that viral-pneumococcal coinfections are common and play a role in the pathogenesis of pneumococcal respiratory infections. To demonstrate the contribution of viruses to the burden of pneumococcal disease specifically, and pneumonia in general, it would be necessary to eliminate one or more of the respiratory viruses. Shortly after the start of the COVID-19 pandemic, multiple reports demonstrated reduced IPD and LRI rates among young children, coincident with dramatically reduced rates of the PDA-viruses globally. Initially, the reduced rates of pneumococcal disease were attributed to non-pharmaceutical interventions that might reduce pneumococcal transmission in the community. However, continuous, virtually unchanged pneumococcal carriage rates were reported in multiple studies, strongly suggesting the reduced circulation of *S. pneumoniae* was not significantly contributing to disease reduction. Surprisingly, pneumococcus-associated diseases and PDA-viruses simultaneously re-emerged in 2021 during the off-season. In contrast to PDA-viruses, other viruses, such as adenovirus and rhinovirus did not show any of the patterns discussed above. We searched PubMed on June 1st, 2022, for studies since 2020 using the following terms: (“COVID-19” or “SARS-Cov-2”) and (“*S. pneumoniae*” or “pneumococcus” or “IPD” or “respiratory virus” or respiratory syncytial virus” or “hMPV” or “influenza” or “parainfluenza” or “adenovirus” or “rhinovirus” or “lower respiratory infection”). The search was for English literature and unrestricted by date. **Added value of this study** Three unique characteristics of the COVID-19 pandemic-induced abnormal dynamics, coupled with multiple ongoing cohort studies in young children, contributed to the historic opportunity to model and quantify the attributable role of the various common respiratory viruses to four pneumococcus-associated disease endpoints (in particular community-acquired alveolar pneumonia (CAAP), non-CAAP LRIs, pneumococcal bacteremic pneumonia and non-pneumonia IPD): First, the full seasonal disappearance of all PDA-viruses shortly after the start of the pandemic, in the presence of continuous, uninterrupted pneumococcal carriage and continuous unchanged rhinovirus and adenovirus activity. Second, the off-season resurgence of the PDA-viruses in 2021. Third, the sequential, rather than simultaneous, re-emergence of the PDA-viruses. The analysis in this study suggests that several of the respiratory viruses, particularly RSV and hMPV, play an important causative role in the pathogenesis of pneumococcal diseases and other respiratory infections. Furthermore, the proportion attributable to each of the PDA-viruses for each of the four studied disease endpoints, and each of the age groups (<1, 1, and 2-4 years of age) could be estimated. **Implication of all the available findings** Our findings add evidence about the absolute and relative contribution of common respiratory viruses to the burden of pneumonia and pneumococcal diseases and related conditions in young children. The strong contribution of RSV to disease burden compared to other viruses in all studied disease endpoints suggests that interventions that target viruses could have secondary effects on the burden of diseases typically attributed to bacteria. ## Introduction Respiratory viruses are major contributors to the burden of pneumonia in general and disease caused by pneumococcus specifically. Respiratory syncytial virus (RSV), human metapneumovirus (hMPV), influenza viruses (flu), and parainfluenza viruses (PIV) have been associated with both community-acquired alveolar pneumonia (CAAP; considered mostly bacterial, in particular pneumococcal),1, 2 and invasive pneumococcal disease (IPD).1-12 However, for other common viruses, such as adenoviruses (AdV) and rhinoviruses (RhV), there is little evidence of a causal association with pneumonia or pneumococcal disease in young children.13 There is compelling evidence that pneumococcus, a bacterial pathogen, contributes to the burden of viral pneumonia and non-pneumonia lower respiratory tract infections (LRIs) in young children. This was seen in declining rates of hospitalisation for viral pneumonia following use of pneumococcal conjugate vaccines (PCVs), both in clinical trials and observational studies10, 14, 15 This supports the notion that pneumococcus may increase the severity of viral infections. However, there is limited evidence about whether reducing or eliminating the activity of specific viruses could reduce severity or rates of pneumonia or pneumococcal disease. The COVID-19 pandemic provided a unique opportunity to examine this question. Shortly after the emergence of the COVID-19 pandemic, temporary elimination of the circulation of RSV, hMPV, flu, and PIV occurred globally, while RhV and AdV continued to circulate.13, 16 At the same time, a drop in pneumococcal disease rates was seen.17-21 This was initially believed to result from reduced pneumococcal circulation, caused by the mitigation measures adopted during the COVID-19 pandemic. However, nasopharyngeal carriage rates of pneumococcus among young children were largely unchanged during the pandemic, suggesting that this was not the main factor behind the decline in rates of pneumococcal disease.12, 13, 16, 22 An alternative possibility, based on previous observations of bacterial-viral coinfections, is that the dramatic suppression of the respiratory viruses led to the decline in rates of IPD. Starting in 2021, an off-season sequential re-emergence of the suppressed viruses occurred in Israel, starting with PIV, and followed by hMPV, RSV, and flu. This provided a unique opportunity to examine the association of these viruses with rates of community acquired alveolar pneumonia (CAAP), IPD and non-CAAP LRIs. We hypothesized that RSV, hMPV, flu, and PIV (collectively termed in the current study pneumonia and pneumococcus disease-associated viruses [PDA-viruses]) are important determinants of these disease endpoints in young children. In a previous study we described events that occurred from the start of the pandemic in 2020 through February 2021, mainly describing the period when the PDA-viruses were absent. The current study extends the period to the end of 2021, which includes the sequential re-emergence of the PDA-viruses and pneumococcus-associated disease and adds a model to quantify the respective role of individual viruses in the dynamics of each of the four studied endpoints by age group. This study provides a comprehensive picture of the unusual dynamics of the respiratory viruses, pneumococcal disease, and other respiratory disease endpoints. ## Methods ### Setting The Soroka University Medical Center (SUMC) is the only hospital in the Negev, southern Israel, providing primary health care to >95% of the children living in the region, enabling incidence calculations. Two distinct paediatric ethnic populations reside side-by-side in southern Israel: The Bedouin population, with infectious disease epidemiology, including pneumococcal serotype distribution, similar to lower-middle income countries; and the Jewish population, whose epidemiology resembles higher-income Western countries. Higher rates of respiratory and invasive disease, and pneumococcal carriage were reported among Bedouin than among Jewish children.23 In 2020, there were ∼96,500 children <5 years in the Negev; ∼50% were Bedouin.24 Implementation of PCVs in the National Immunization Plan occurred in 2009 (PCV7) and 2010 (PCV13); since 2013 ∼90% of both Jewish and Bedouin children have received all three required doses.25 ### The COVID-19 pandemic in 2019-2021 The details of the COVID-19 pandemic and associated mitigation measures are described in **Supplementary Figure 1**. ### Study design The study population for respiratory disease and virus-testing comprised all children living in the region from January 2016 through December 2021. For IPD, we included all children nationwide. All disease databases included children <5 years old. However, for pneumococcal carriage, we included children <3 years since the complete database of children without respiratory symptoms has been created only for children in this age group. The data were derived from multiple ongoing, prospective long-standing cohort surveillance programs. Monthly incidence rates of the various disease endpoints and prevalence rates of pneumococcal nasopharyngeal carriage and detection of respiratory viruses during 2020-2021 were compared to the corresponding values during 2016-2019. The total numbers for each database and the comparison of demographic characteristics are presented in **Table 1**. View this table: [Table 1:](http://medrxiv.org/content/early/2022/12/01/2022.09.06.22279606/T1) Table 1: Demographics of the study population #### 1) Surveillance of community-acquired alveolar pneumonia (CAAP) and non-CAAP LRI requiring chest-radiography in children <5 years of age, southern Israel This is an ongoing, prospective, population-based, active surveillance system, initiated in 2002. The current report deals with data from January 2016 through December 2021. All children <5 years residing in the Negev district, seen at the SUMC paediatric emergency room (PER) (either ambulatory or hospitalised) who had chest radiological examination (CXR) within <48 hours from admission were included in the analysis. Clinical practice in the SUMC PER regarding LRI referral, evaluation, and the need for CXR examination did not change throughout the study period. Chest radiographs were analysed according to the World Health Organization (WHO) Standardization of Interpretation of Chest Radiograph Working Group.26, 27 All chest radiographs were recorded daily and were evaluated separately by the 2 paediatric infectious disease specialists who each read all the chest radiographs, independently. Further evaluation was done by a paediatric radiologist who was unaware of the clinical data and of the paediatricians’ evaluation. CAAP diagnosis was confirmed by agreement on the presence of alveolar infiltrates or pleural effusion of at least one of the study paediatric infectious disease specialists and the paediatric radiologist.28 Since November 2020 the paediatric radiologist was replaced by a third paediatric infectious disease specialist. For this period the diagnosis of CAAP needed agreement of at least two of the three readers. The diagnosis of non-CAAP LRI included all visits of LRI with CXR examination that did not show alveolar infiltrates. This group included both: 1) CXR with readings of infiltrates that were not alveolar (<5% of the group), including linear and patchy densities (interstitial infiltrates) in a lacy pattern involving both lungs, featuring peri bronchial thickening and multiple areas of atelectasis; or minor patchy infiltrate that are not of sufficient magnitude to constitute alveolar pneumonia, and small areas of atelectasis;26 and 2) absence of consolidation, other infiltrates or pleural effusion (but including CXR with hyperinflation, >95% of this group). All non-CAAP-LRIs were grouped together because previous studies showed a very low inter-reader consistency.29 #### 2) Nationwide Invasive Pneumococcal Diseases (IPD) Surveillance This has been an ongoing, nationwide, prospective, population-based, active surveillance system, initiated in 1989. Surveillance is conducted by the Israeli Pediatric Bacteremia and Meningitis Group (IPBMG) in 27 medical centres routinely obtaining blood cultures from children: All 26 hospitals admitting children and one major outpatient health maintenance organisation (HMO, Maccabi Healthcare Services) central laboratory. Less than 1% of blood cultures in all ages and no cerebrospinal fluid (CSF) cultures are obtained outside these centres. IPD episodes were defined as illness episodes during which *S. pneumoniae* was isolated from blood or CSF. Study population: All children <5 years in Israel. Non-culture diagnoses (polymerase chain reaction [PCR], antigen testing, gram stain or clinical diagnosis only) were excluded. Positive cultures from non-blood or CSF sterile sites (i.e., joint/pleural fluid) were excluded. Local investigators in each centre responded to a questionnaire sent monthly by the principal investigator, located at the Pediatric Infectious Disease Unit (PIDU), of the SUMC, which served as the study headquarters. Several measures were used to ensure consistent data collection; (1) weekly contact with all 27 laboratories; (2) assuring transportation of all *S. pneumoniae* isolates to the study headquarters’ laboratory; and (3) monthly contact with the local investigators. Completed reports included the following data: Isolate source (blood/CSF), culture date, birthdate, sex, ethnicity (Jewish/non-Jewish), main diagnoses, endpoint (mortality) and hospitalisation duration. IPD cases reported from each site have been constantly compared to the list of isolates obtained by the reference Ministry of Health laboratory. Since July 2009, >95% of the cases have been retrieved.30, 31 #### 3) Pneumococcal nasopharyngeal (NP) carriage surveillance The detailed nasopharyngeal testing methodology was recently described.25 Nasopharyngeal samples were obtained using a flexible dacron-tipped swab, introduced through the nostrils and advanced until resistance was found. These swabs were inoculated into modified Stewart transport medium (Medical Wire and Equipment Co., Ltd., Corsham, England) and processed within 16 hours at the SUMC Clinical Microbiology Laboratory. Material from swabs was plated on Columbia agar with 5% sheep blood and 5.0 μg/mL gentamicin and incubated aerobically at 35°C in a CO2-enriched atmosphere for 48h. The presumptive identification of *S. pneumoniae* was based on the presence of α-hemolysis and the inhibition by optochin, and the identity of the bacteria present was confirmed by a positive slide agglutination test result (Phadebact; Pharmacia Diagnostics). One *S. pneumoniae* colony per plate was then subcultured, harvested, and kept frozen at -70°C for further testing. Density analysis was performed at the time of initial sample processing. A semiquantitative plating method was used and graded 0 (negative) to 4 (growth in all plate quadrants). The NP swabs were plated by rolling the swab over one quarter of the plate and streaking the sample onto 4 quadrants using a sterile loop. Growth was termed 1+, 2+, 3+, or 4+ when colonies were seen in 1, 2, 3, or 4 quadrants, respectively.32 ##### a) Children <3 years without respiratory infections visiting the SUMC PER This is an ongoing, prospective, population-based, active surveillance system, initiated in November 2009. Each working day, nasopharyngeal cultures were obtained from the first 4 Jewish and 4 Bedouin children <3 years old, resident of the Negev region, presenting at the PER, after a written informed consent was signed by the parents. Children not residing in the Negev region were excluded. In cases where multiple swabs were obtained from a child during the same month, only the first culture was included25 If <8 children were available, nasopharyngeal carriage evaluation was obtained from randomly selected hospitalized children. For the current analysis, we excluded those with any diagnosis of respiratory infections. These children constituted 3114/7921 (34·3%) of tested children during 2016-2019. In 2020 and 2021 these children constituted 331/1190 (27·8%) and 457/1323 (34·5%), respectively, of all children tested for carriage. ##### b) Healthy children <3 years presenting for vaccinations Since 2011, nasopharyngeal cultures have been obtained from healthy children <3 years old, presenting to the Maternal and Child Health Care Centres in southern Israel for vaccination. A nasopharyngeal swab was obtained after a written informed consent was signed by the parents.33 These children constituted 4,807/7,921 (60·7%) of all tested children in 2016-2019. In 2020 and 2021 they constituted 859/1,190 (72·2%) and 866/1,323 (65·5%), respectively. #### 4) Nasopharyngeal detection of respiratory viruses Nasopharyngeal specimens for respiratory viruses were obtained from hospitalised children following the request of the treating physician for clinical indication. Until February 2019 the specimens were nasopharyngeal washes and the procedures, including the laboratory method for PCR testing, were previously published.34 The viruses tested were RSV, influenza A and influenza B viruses, parainfluenza viruses, adenovirus and human metapneumovirus. Since February 2019, nasal swabs were obtained and the detection method was switched to a commercial kit (Seeplex RV7 detection kit, Seengene Inc., Korea). ## Ethics The study and all its components were approved by the SUMC Ethics Committee with the following approval numbers: SOR-0405-16; SOR-0220-20; SOR-0062-11; 10374; 4908, 3075. ## Statistical Analysis Since respiratory infection rates, relative importance of viruses and clinical presentations differ by age and ethnic group, monthly incidence rates for all endpoints were stratified by age (<1, 1, 2-4 years) and ethnicity (Bedouin, Jewish). Viruses were divided into 2 groups: PDA-viruses (RSV, hMPV, flu, and PIV) and non-PDA-viruses (RhV and AdV). For the modelling we analysed each virus separately, even if they were detected as coinfections. The monthly proportion of positive pneumococcal nasopharyngeal cultures were calculated out of all collected samples for children <1 and 1-2 years, except during April-May 2020, due to disruptions in sample collection caused by a strict lockdown. No difference in the age distribution of the population <5 was observed during the study years. Ratios of incidence rates and proportions were calculated to compare the COVID-19 and pre-COVID-19 periods. All rate ratios were adjusted for ethnicity and age using the Cochran-Mantel-Haenszel method. The statistical significance threshold was p <0·05. Data were analysed using R 4.0.2.35 ### Modelling The goal was to estimate the proportion of each disease endpoint attributable to the different viruses. To do this, we developed hierarchical Bayesian regression models with a negative binomial likelihood and identity link.36 The clinical outcomes evaluated were CAAP, non-CAAP LRI, bacteremic pneumococcal pneumonia and non-pneumonia IPD. For each endpoint, we have a time series of the number of clinical cases, by month, stratified by age category (<1, 1, 2-4 years) and ethnicity (Jewish *vs*. Bedouin for CAAP and non-CAAP; Jewish *vs*. non-Jewish for IPD). The covariates are a linear trend for time, 12-month harmonic variables to capture unexplained seasonal variations, and the proportion of tests positive for each virus in that month for the corresponding ethnicity (combined across all three age groups). We use an identity link, rather than a log-link so that the effects of each virus would be additive. This model generally follows the structure used by Zheng et al..36 The magnitude of the effects of each virus and the seasonal and trends terms were allowed to vary by age and ethnicity. For age group a, ethnicity e, time point t, the model structure was ![Formula][1] Yt,a,e denotes the number of clinical cases for each time point *t*, age group *a*, and ethnicity *e. λ* t,a,e represents the expected value of Yt,a,e with a variance of *λ**ijk*(1 + *λ**ijk*/*r*) and an overdispersion parameter *r*. b0a,e represents the intercept for the group, b1a,e the time trend for the group, b3a,e and b4a,e capture the effect of baseline seasonality, and b4a,e – b8a,e represent the effect of each virus by group. After fitting the models, the contribution of each virus was estimated from the fitted parameters. *λ**t,a,e* was recalculated by setting the bK*a,e* term corresponding to the viruses to 0 and comparing these predicted values with the predicted values base on the fitted values of the parameters. Predicted values were summed by age group or overall to obtain aggregate estimates of the contribution of each virus. The viral time series were scaled to vary between 0 and 1. Each group-level effect was centred around a global effect, which allowed the coefficient to vary by group ![Formula][2] for the trend and harmonic coefficients (b1-b3), ![Formula][3] while for the intercept and for the effects of viruses, the effects were restricted to be non-negative: ![Formula][4] The model was fit in a Bayesian framework. The variance components (***τ***k) were given a uniform[0,100] prior on the standard deviation (e.g., sqrt(1/ ***τ***k2). ***β***K were assigned minimally informative priors of Normal(0, 1). All model fitting was performed using rjags37 with three separate Markov chains and a burn-in period of 10,000 iterations (i.e., before convergence of the models) in each chain. Posterior inference was based on 30,000 samples (10,000 from each chain). The analysis code can be found at [https://github.com/DanWeinberger/Israel\_pneumo_covid](https://github.com/DanWeinberger/Israel_pneumo_covid). Since RhV testing was initiated only in February 2019, we did not include this virus in the model. However, a separate regression analysis, adjusting for the same set of covariates as above, for the period of February 2019 through December 2021 showed no correlation of RhV to any of the studied clinical endpoints. ## Results ### I. Dynamics of Disease During 2016 through December 2021, a total of 3,204, 26,695, 257, and 619 episodes of CAAP, non-CAAP LRI, pneumococcal bacteremic pneumonia, and non-pneumonia IPD, respectively, were recorded. The highest incidence occurred in children <1 year. (**Table 1, Supplementary Tables 1, 2**) The dynamics of the four disease endpoints during the early pandemic period (through February 2021) were described in detail in a previous publication,13 and are included as part of **Figures 1 and 2; and Supplementary Tables 3 and 4**. In brief, broad declines in all four endpoints were observed. The largest decline in incidence was seen in CAAP and pneumococcal bacteremic pneumonia, while for non-CAAP-LRIs, although statistically significant, the reduction was less pronounced. For non-pneumomia IPD, the decline was not statistically significant. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/12/01/2022.09.06.22279606/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2022/12/01/2022.09.06.22279606/F1) Figure 1: Monthly dynamics of community-acquired alveolar pneumonia (CAAP) **(a)**, pneumococcal bacteremic pneumonia **(b)**, non-community-acquired alveolar pneumonia lower respiratory infections (non-CAAP LRIs) **(c)**, and non-pneumonia invasive pneumococcal disease (non-pneumonia IPD) **(d)**. The red line represents incidences in 2020 and 2021, the blue line shows the mean incidence and the shaded blue area the minimum and maximum incidence during 2016 through 2019. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/12/01/2022.09.06.22279606/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2022/12/01/2022.09.06.22279606/F2) Figure 2: Incidence rates ratio (±95% CI) by quarters comparing incidence of 2020 and 2021 to the mean incidence of 2016 through 2019; community-acquired alveolar pneumonia (CAAP) **(a)**, pneumococcal bacteremic pneumonia **(b)**, non-community-acquired alveolar pneumonia lower respiratory infections (non-CAAP LRIs) **(c)**, and non-pneumonia invasive pneumococcal disease (non-pneumonia IPD) **(d)**. Black dots are quarterly IRRs; the vertical lines show 95% confidence intervals. The numerical values are presented in **Supplementary Table 3 and 4**. Starting in spring 2021, an off-season, abrupt surge was observed with unexpectedly high rates, particularly for the 3 respiratory endpoints (CAAP, non-CAAP LRIs, and pneumococcal bacteremic pneumonia). In all endpoints, the off-season rates reached a magnitude similar to those observed during 2016-2019. (**Figures 1, 2, Supplementary Tables 3, 4**) The dynamic patterns of all disease endpoints were similar when examined by age (<1, 1, and 2-4 years). (**Supplementary Figures 2-4**) ### II. Dynamics of pneumococcal nasopharyngeal carriage in children <3 year A total of 7,921, 1,190, and 1,323 samples were obtained in 2016-2019, 2020, and 2021, respectively. (**Table 1**) During April and May 2020, we interrupted surveillance activities because of the severe lockdown but resumed from June 2020. (**Figure 3**) ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/12/01/2022.09.06.22279606/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2022/12/01/2022.09.06.22279606/F3) Figure 3: **(a)** Monthly dynamics of the proportion of all tested nasopharyngeal swabs that were positive for *S. pneumoniae* in children <3 years of age, in 2020 and 2021 compared to 2016 through 2019. The red line represents the proportion of 2020 and 2021, the blue line shows the mean proportion and the blue shaded area the minimum and maximum proportion during 2016 through 2019; **(b)** Dynamics of pneumococcus positive proportions of swab samples by semi quantitatively determined density: low density (1-2) represented by the blue line *vs*. high density (3-4), represented by the red line. During 2016-2019, the mean monthly proportion of carriers ranged from 34·3% to 56·7% in 2016-2019, 29·0% to 52·1% in 2020, and 31·3% to 45·6% in 2021. Although from September 2020 to February 2021, carriage was somewhat lower than during 2016-2019 (*P*=0·067) (**Figure 3a, Supplementary Tables 5, 6**), no overall trend was observed, and changes in carriage prevalence did not correlate with changes in any of the disease endpoints. Furthermore, the semi-quantitative analysis showed a decline in lower-density colonization. but the prevalence of high-density colonization remained stable. (**Figure 3b, Supplementary Table 7**) No significant changes in serotype distribution were found during the pandemic. (**Supplementary Table 8, Supplementary Figure 5**) ### III. Dynamics of viral activity in children <5 years A total of 27,147 nasopharyngeal samples were obtained for virus detection, 17,555; 3,444; and 6,143 in 2016-2019, 2020, and 2021, respectively. (**Table 1**) Of these, 41·7%, 58·0%, and 63·4%, respectively, were positive for ≥1 virus: RSV 4,024, AdV 3,049, PIV 1,143, flu A 1,026, hMPV 883, and flu B 308. Testing for RhV was only initiated in February 2019 with 4,845 positive samples. (**Supplementary Table 9**) Co-detection of ≥2 viruses was common, especially starting in February 2019, when RhV testing was added. However, a co-detection of any PDA-virus with another PDA-virus was rare (1·3%). **(Supplementary Table 10)** Strikingly, from April 2020 through October 2020, none of the PDA-viruses were detected. (**Figure 4, Supplementary Figure 6**) However, since November 2020, an off-season sequential re-emergence of the PDA-viruses was observed, each with a peak reaching an equal or higher magnitude compared to the pre-pandemic seasonal peaks: PIV (November 2020), hMPV (March 2021), RSV (May 2021), flu A (H3N2 only; September 2021). However, flu B was not detected from March 2020 through the end of the study. In contrast, although RhV and AdV had a short nadir during the time of the first lockdown (April-May 2020), they were frequently detected throughout the study. From March 2021, AdV had a higher-than-usual activity. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/12/01/2022.09.06.22279606/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2022/12/01/2022.09.06.22279606/F4) Figure 4: Dynamics of monthly numbers of all virus-positive nasal samples in children <5 years in southern Israel for specific viruses in 2020 and 2021 compared to 2016 through 2021. The red line represents the number of positive samples in 2020 and 2021, the blue line shows the mean number of positive samples and the shaded blue area the minimum and maximum number of positive samples during 2016 through 2019. ### IV. RSV and hMPV are major contributors to the burden of LRI and IPD Due to the striking similarity of the dynamics of all 4 disease endpoints and the dynamics of the PDA-viruses, we sought to quantify the contribution of each virus. We first fit a model to CAAP data from the pre-pandemic period only. The main contributor was RSV (32%), followed by hMPV (5%), flu (4%), PIV (3%), and AdV (1%). (**Supplementary Table 11**) Projecting the model forward based on observed virus activity, this model accurately captured the surge in CAAP cases in 2021 among children <1 year but underestimated the number of CAAP cases among children aged 1 and 2-4 years (**Supplementary Figure 7**), suggesting that the model under-estimated the importance of the viruses in the older age groups. Since the activity of these viruses typically peaks around the same time in winter, it can be difficult to accurately attribute changes to specific viruses. To disentangle the contribution of the viruses to the burden of pneumonia, we leveraged the different timing of the re-emergence of the different viruses in 2020-2021 and refit the models to the entire dataset. (**Figure 5, Table 2**) With these models, we estimated that 82% of CAAP cases could be attributed to the common respiratory viruses. The main contributor to this burden was RSV (49%), followed by hMPV (13%), PIV (11%), and flu (7%). AdV activity was not associated with variations in CAAP. The importance of the different viruses varied somewhat by age with smaller contributions overall from the viruses in older age groups. (**Figure 6, Table 2**) RSV was a particularly large contributor (62%) to CAAP among children <1 year. ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/12/01/2022.09.06.22279606/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2022/12/01/2022.09.06.22279606/F5) Figure 5: Observed (dark grey lines) and model-fitted values (black line) for the number of cases of community-acquired alveolar pneumonia (CAAP), non-community-acquired alveolar pneumonia lower respiratory infections (non-CAAP LRIs), pneumococcal bacteremic pneumonia, and non-pneumonia invasive pneumococcal disease (non-pneumonia IPD), based on a model fit to the entire data set. The red dashed line and shaded area represent the estimate for how many cases there would have been in the absence of viruses (±95% CI). ![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/12/01/2022.09.06.22279606/F6.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2022/12/01/2022.09.06.22279606/F6) Figure 6: Estimated contribution of different viruses to the burden of community-acquired alveolar pneumonia (CAAP) among children <5 years of age in southern Israel. Parainfluenza – yellow; hMPV – green; influenza – blue; RSV – orange; other viruses/non – grey. View this table: [Table 2:](http://medrxiv.org/content/early/2022/12/01/2022.09.06.22279606/T2) Table 2: Percent of CAAP attributable to RSV, hMPV, Flu, AdV, PIV (adjusted for ethnicity), stratified by age The viruses also made important contributions to the other endpoints with 31% of non-CAAP cases, 29% of IPD pneumonia and 22% of non-pneumonia IPD associated with viruses. (**Table 2, Supplementary Figure 9**) For both non-CAAP and IPD pneumonia, RSV was the most important contributor, 21% and 18%, respectively. For non-pneumonia IPD, none of the PDA-viruses, individually, contributed significantly. AdV did not contribute to any of the endpoints at any age. ## Discussion The COVID-19 pandemic has temporarily but strikingly modified the epidemiology of both PDA-virus activity and the studied clinical endpoints in which pneumococcus plays an important causative role. These clinical endpoints declined sharply in 2020 and returned in 2021, together with the re-emergence of the PDA-viruses. These dramatic changes occurred in the presence of virtually unchanged nasopharyngeal carriage of pneumococcus, strongly suggesting its continuous circulation in the community. While PDA-virus activity was associated with the studied disease endpoints, that of RhV and AdV was not. A high proportion of the burden of CAAP and non-CAAP LRIs was associated with the activity of PDA-viruses, especially during the first year of life. We also showed that the modified dynamics of both disease and virus activity during the COVID-19 pandemic helped to improve estimates of the PDA-virus contribution to the studied disease endpoints. Finally, we showed that RSV had the highest attributable contribution to all studied endpoints, except for non-pneumonia IPD. Multiple sites have described the abrupt reduction of pneumococcal diseases (mostly IPD) and all-cause respiratory diseases shortly following the start of the COVID-19 pandemic. The initial speculation was that the widely undertaken non-pharmaceutical interventions largely resulted in a reduction of pneumococcal transmission.17, 18 However, continuous active surveillance of pneumococcal carriage in young children showed normal or close to normal carriage rates, and no association with IPD or pneumococcal pneumonia could be demonstrated.8, 13, 19, 22 A recent report from Vietnam, while confirming no association of overall nasopharyngeal carriage during NPI, noted a marked decrease of the density of pneumococcus in carriage.38 In contrast, we found no reduction of carriage density, although measured by semi-quantitative methods. In addition, no change in serotype distribution was found, similar to findings from France.12 Therefore, it is unlikely that changes in carriage played an important role in the drastic disease incidence dynamics during COVID-19. No previous data on the carriage dynamics during the resurgence of pneumococcus-associated diseases have been reported. It is not clear why the PDA-viruses largely disappeared during the initial phase of the pandemic, while other viruses and pneumococci continued to circulate in the presence of the various mitigation measures. It is possible that pneumococcus, AdV and RhV rely on different age groups or routes for transmission, they could be more transmissible, or additional mechanism such as viral-viral interactions could play a role. During the study period, the dynamic patterns of the 4 studied disease endpoints were similar but not identical. The most affected rates were those of CAAP, with similar dynamics seen for pneumococcal bacteremic pneumonia. However, the smaller number of endpoint events of the latter resulted in a higher uncertainty range. Non-CAAP LRIs were also reduced during the first year of the pandemic, but the reduction was only 50% compared to that of CAAP. The non-pneumonia IPD incidence reductions, although exhibiting similar trends, did not reach significance in any of the time points, when compared to those in 2016-2019, consistent with the absence of demonstrable association with the PDA-viruses. The contribution to the various disease endpoints differed between individual viruses, and between age groups. By far, the most important contributor to both CAAP, non-CAAP LRIs, and pneumococcal bacteremic pneumonia was RSV. In children <2 years, over half of the cases were attributable to RSV. Its contribution to non-CAAP LRIs was only around half compared to that in CAAP. This was somewhat unexpected, due to the common association with bronchiolitis. However, this difference suggests a potentially significant pathogenic role of RSV as a coinfecting agent with bacterial pathogens causing CAAP. The extent of pneumococcal involvement in CAAP and non-CAAP-LRIs is still not fully clarified. This cannot be directly answered since detection of pneumococcus in the respiratory tract does not necessarily imply its causative role in CAAP. However, the profound impact of PCV on CAAP and even non-CAAP LRIs provides a powerful tool for inference on the likely causative role of *S. pneumoniae* in these disease endpoints.14, 39 Studies from Israel and elsewhere, have demonstrated a ∼50% reduction in CAAP in young children following PCV implementation,14, 28 thus strongly suggesting an extensive causative role of *S. pneumoniae* in CAAP during the pre-PCV era. Our non-CAAP LRI cases were severe enough to warrant a chest radiograph, but per definition, those with CAAP were excluded. Studies have shown that such cases have a more diverse aetiology. Even so, a rate reduction of ≥25% was observed in young children post PCV implementation, including a 34% reduction in our region.14 Using a similar logic as per CAAP, it is plausible that vaccine-serotype pneumococci also played an important role in non-CAAP LRIs, although of a lesser magnitude than in CAAP. The similarity in dynamics between CAAP and non-CAAP LRIs is therefore not surprising. Second to RSV, the estimated contribution of hMPV in CAAP is notable, contributing to approximately one in seven CAAP episodes in all age groups. The role of hMPV in CAAP has been well recognized.11, 40-42 However, in the current study we could compare its role in CAAP to that in other studied endpoints. Surprisingly, no significant contribution to any of the non-CAAP endpoints could be demonstrated. The specific mechanisms by which hMPV contributed mainly to CAAP are not elucidated. The attributable proportions of flu and PIV to both CAAP and non-CAAP, although significant, were modest, and the relative contribution of PIV increased with age. The lower relative contribution of flu to hospitalisation for LRIs compared to RSV and hMPV was described previously41 and is in contrast to reports from adults.43 We examined separately the dynamics of pneumococcal bacteremic pneumonia and non-pneumonia IPD. We have previously shown that after PCV7/PCV13 implementation in Israel, the dynamics of the two endpoints differed by specific pneumococcal serotype proportions and seasonality. The seasonality of bacteremic pneumonia was typical of respiratory viruses, while that of non-pneumonia IPD was not.2 In the current study we estimate that RSV played a significant role in pneumococcal bacteremic pneumonia, but its role in non-pneumonia IPD was negligible. This suggests some differences in the pathogenesis, further supporting previous observations.44 This might also explain, at least in part, the lack of significant reduction of the non-pneumonia IPD in the presence of significant reduction of the other three analysed endpoints. The rates of these viruses were completely unrelated to those of the studied disease endpoints. Detection of RhV and AdV is common, both in healthy children and children with LRI, and their attributable role to LRI is often questionable.45-47 Furthermore, the co-detection of RhV and AdV with other viruses is common, and the co-presence with any of the PDA-viruses usually bears the clinical presentation typical of the specific coinfecting PDA-virus. In contrast, we found that co-detection of ≥2 PDA-viruses was rare (1·3%). This study has a number of important strengths. We leverage unique and high-quality datasets that cover a long time period including pre-pandemic years. Together with the natural disruptions caused by the pandemic, these data provide a comprehensive understanding of the contributions of respiratory viruses to the aetiology of pneumococcal disease and other related respiratory diseases. These data are integrated using regressions models to estimate the contribution of each virus to the observed variations and estimates from models fit exclusively to the pre-pandemic period and those fit to the entire study period providing a range of estimates for the contribution of each virus. The analyses also have important limitations. We did not account for the possible effects of coinfections or interactions between viruses. The data on viral activity are also limited by being conducted in the hospital setting. However, we expect that viral activity in the hospital likely reflects viral activity in the community. Finally, the causal relationship of the PDA-viruses with the studied endpoints could be confounded by a third causative factor. However, given the high degree of association, the multiple endpoints studied, and the results of the model, this possibility seems unlikely. In conclusion, this study leverages the unique disruptions of pathogen circulation that followed COVID-19 mitigation measures to gain a deeper understanding of the role of respiratory viruses in driving the burden of pneumococcal disease and related conditions. These analyses suggest that several of the respiratory viruses, particularly RSV and hMPV, play a critical role in pathogenesis in children. This suggests that interventions that target these viruses could have secondary effects on the burden of disease typically attributed to bacteria. ## Supporting information Supplemenatry Tables [[supplements/279606_file03.pdf]](pending:yes) Supplementary Figures [[supplements/279606_file04.pptx]](pending:yes) ## Data Availability All data produced in the present study are available upon reasonable request to the authors ## Contributors RD was responsible for concept, design, planning of the study; acquisition, analysis, and interpretation of the data; drafting and critically revising the manuscript. BAVDB contributed to the concept of the study, analysis, and interpretation of the data, as well as reviewing the manuscript for important intellectual content. SB-S was involved in the planning of the study and critically reviewed the manuscript for important intellectual content. DG reviewed the manuscript for important intellectual content. YS-A contributed to the acquisition of the data and reviewed the manuscript for important intellectual content. DMW participated in the design of the study, analysis and interpretation of the data (modelling component), and reviewed the manuscript for important intellectual content. DD contributed to the design of the study, analysis, and interpretation of the data, as well as reviewing the manuscript for important intellectual content. ## Data sharing statement All relevant data will be shared by the lead contact upon request. ## Declaration of Interest RD has received grants from Pfizer, MSD, MedImmune/AstraZenaca. He serves as a scientific consultant and on the advisory board of Pfizer and MSD. He also is part of the speakers’ bureau of Pfizer, MSD, Sanofi Pasteur, and GSK. SB-S has received grants from Pfizer and serves as a scientific consultant and on the advisory board of Pfizer and MSD. He also is part of the speakers’ bureau of Pfizer, MSD, and GSK. DG has received grants from MSD, and he serves as scientific consultant for Pfizer, MSD, and GSK where he is also part of the speakers’ bureau. DMW has received consulting fees from Pfizer, Merck, Affinivax, Matrivax, and GSK and is principal investigator on grants from Pfizer and Merck to Yale University. DD has received grants from Pfizer. ## Footnotes * Text from the Appendix was moved to the main text. * Received September 6, 2022. * Revision received December 1, 2022. * Accepted December 1, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. 1.Ampofo K, Bender J, Sheng X, Korgenski K, Daly J, Pavia AT, et al. Seasonal invasive pneumococcal disease in children: role of preceding respiratory viral infection. Pediatrics. 2008;122:229–37. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1542/peds.2007-3192&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18676537&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F01%2F2022.09.06.22279606.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000258142500001&link_type=ISI) 2. 2.Ben-Shimol S, Greenberg D, Hazan G, Shemer-Avni Y, Givon-Lavi N, Dagan R. Seasonality of both bacteremic and nonbacteremic pneumonia coincides with viral lower respiratory tract infections in early childhood, in contrast to nonpneumonia invasive pneumococcal disease, in the pre-pneumococcal conjugate vaccine era. Clin Infect Dis. 2015;60:1384–7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/civ023&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25595749&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F01%2F2022.09.06.22279606.atom) 3. 3.Chiavenna C, Presanis AM, Charlett A, de Lusignan S, Ladhani S, Pebody RG, et al. Estimating age-stratified influenza-associated invasive pneumococcal disease in England: A time-series model based on population surveillance data. PLoS Med. 2019;16:e1002829. 4. 4.Greenberg D, Givon-Lavi N, Faingelernt Y, Ben-Shimol S, Avni YS, Bar-Ziv J, et al. Nasopharyngeal Pneumococcal Carriage During Childhood Community-Acquired Alveolar Pneumonia: Relationship Between Specific Serotypes and Coinfecting Viruses. J Infect Dis. 2017;215:1111–6. 5. 5.Hernandez S, Munoz-Almagro C, Ciruela P, Soldevila N, Izquierdo C, Codina MG, et al. Invasive Pneumococcal Disease and Influenza Activity in a Pediatric Population: Impact of PCV13 Vaccination in Pandemic and Nonpandemic Influenza Periods. J Clin Microbiol. 2019;57. 6. 6.Murdoch DR, Jennings LC. Association of respiratory virus activity and environmental factors with the incidence of invasive pneumococcal disease. J Infect. 2009;58:37–46. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jinf.2008.10.011&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19042025&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F01%2F2022.09.06.22279606.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000263624100006&link_type=ISI) 7. 7.O’Brien KL, Walters MI, Sellman J, Quinlisk P, Regnery H, Schwartz B, et al. Severe pneumococcal pneumonia in previously healthy children: the role of preceding influenza infection. Clin Infect Dis. 2000;30:784–9. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/313772&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10816149&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F01%2F2022.09.06.22279606.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000087557900011&link_type=ISI) 8. 8.Watson M, Gilmour R, Menzies R, Ferson M, McIntyre P, New South Wales Pneumococcal N. The association of respiratory viruses, temperature, and other climatic parameters with the incidence of invasive pneumococcal disease in Sydney, Australia. Clin Infect Dis. 2006;42:211–5. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/498897&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16355331&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F01%2F2022.09.06.22279606.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000233981800008&link_type=ISI) 9. 9.Weinberger DM, Grant LR, Steiner CA, Weatherholtz R, Santosham M, Viboud C, et al. Seasonal drivers of pneumococcal disease incidence: impact of bacterial carriage and viral activity. Clin Infect Dis. 2014;58:188–94. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/cit721&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24190895&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F01%2F2022.09.06.22279606.atom) 10. 10.Weinberger DM, Klugman KP, Steiner CA, Simonsen L, Viboud C. Association between respiratory syncytial virus activity and pneumococcal disease in infants: a time series analysis of US hospitalization data. PLoS Med. 2015;12:e1001776. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1001776&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25562317&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F01%2F2022.09.06.22279606.atom) 11. 11.Wolf DG, Greenberg D, Shemer-Avni Y, Givon-Lavi N, Bar-Ziv J, Dagan R. Association of human metapneumovirus with radiologically diagnosed community-acquired alveolar pneumonia in young children. J Pediatr. 2010;156:115–20. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jpeds.2009.07.014&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19782998&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F01%2F2022.09.06.22279606.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000273200400025&link_type=ISI) 12. 12.Rybak A, Levy C, Angoulvant F, Auvrignon A, Gembara P, Danis K, et al. Association of Nonpharmaceutical Interventions During the COVID-19 Pandemic With Invasive Pneumococcal Disease, Pneumococcal Carriage, and Respiratory Viral Infections Among Children in France. JAMA Netw Open. 2022;5:e2218959. 13. 13.Danino D, Ben-Shimol S, Van Der Beek BA, Givon-Lavi N, Avni YS, Greenberg D, et al. Decline in Pneumococcal Disease in Young Children during the COVID-19 Pandemic in Israel Associated with Suppression of seasonal Respiratory Viruses, despite Persistent Pneumococcal Carriage: A Prospective Cohort Study. Clin Infect Dis. 2021. 14. 14.Ben-Shimol S, Dagan R, Givon-Lavi N, Avital D, Bar-Ziv J, Greenberg D. Use of Chest Radiography Examination as a Probe for Pneumococcal Conjugate Vaccine Impact on Lower Respiratory Tract Infections in Young Children. Clin Infect Dis. 2020;71:177–87. 15. 15.Madhi SA, Klugman KP, Vaccine Trialist G. A role for Streptococcus pneumoniae in virus-associated pneumonia. Nat Med. 2004;10:811–3. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nm1077&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15247911&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F01%2F2022.09.06.22279606.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000223055700033&link_type=ISI) 16. 16.Dagan R, Danino D, Weinberger DM. The Pneumococcus-Respiratory Virus Connection-Unexpected Lessons From the COVID-19 Pandemic. JAMA Netw Open. 2022;5:e2218966. 17. 17.Amin-Chowdhury Z, Collins S, Sheppard C, Litt D, Fry NK, Andrews N, et al. Characteristics of Invasive Pneumococcal Disease Caused by Emerging Serotypes After the Introduction of the 13-Valent Pneumococcal Conjugate Vaccine in England: A Prospective Observational Cohort Study, 2014-2018. Clin Infect Dis. 2020;71:e235–e43. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F01%2F2022.09.06.22279606.atom) 18. 18.Brueggemann AB, Jansen van Rensburg MJ, Shaw D, McCarthy ND, Jolley KA, Maiden MCJ. Changes in the incidence of invasive disease due to Streptococcus pneumoniae, Haemophilus influenzae, and Neisseria meningitidis during the COVID-19 pandemic in 26 countries and territories in the Invasive Respiratory Infection Surveillance Initiative: a prospective analysis of surveillance data. The Lancet Digital Health. 2021;3:e360–e70. 19. 19.Danino D, van der Beek BA, Givon-Lavi N, Greenberg D, Ben-Shimol S, Dagan R. Unraveling the Impact of Pneumococcal Conjugate Vaccines on Ambulatory Antibiotic Drug Consumption in Young Children: An Interrupted Time-Series Analysis. Clin Infect Dis. 2021;73:1268–78. 20. 20.Rybak A, Levy C, Bonacorsi S, Bechet S, Vie le Sage F, Elbez A, et al. Antibiotic Resistance of Potential Otopathogens Isolated From Nasopharyngeal Flora of Children With Acute Otitis Media Before, During and After Pneumococcal Conjugate Vaccines Implementation. Pediatr Infect Dis J. 2018;37:e72–e8. 21. 21.Teng JLL, Fok KMN, Lin KPK, Chan E, Ma Y, Lau SKP, et al. Substantial Decline in Invasive Pneumococcal Disease During Coronavirus Disease 2019 Pandemic in Hong Kong Clin Infect Dis. 2022;74:335–8. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciab382&link_type=DOI) 22. 22.Willen L, Ekinci E, Cuypers L, Theeten H, Desmet S. Infant Pneumococcal Carriage in Belgium Not Affected by COVID-19 Containment Measures. Front Cell Infect Microbiol. 2021;11:825427. 23. 23.Shouval DS, Greenberg D, Givon-Lavi N, Porat N, Dagan R. Serotype coverage of invasive and mucosal pneumococcal disease in Israeli children younger than 3 years by various pneumococcal conjugate vaccines. Pediatr Infect Dis J. 2009;28:277–82. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/INF.0b013e31818e0e2e&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19258924&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F01%2F2022.09.06.22279606.atom) 24. 24.Israel CSAo. Table 3.11 Live births, deaths and infant deaths by district and sub-district, population groups and religion [https://wwwcbsgovil/he/publications/doclib/2018/3%20shnatonvitalstatistics/st03\_11xpdf](https://wwwcbsgovil/he/publications/doclib/2018/3%20shnatonvitalstatistics/st03_11xpdf), accessed November 2022. 2018. 25. 25.Ben-Shimol S, Givon-Lavi N, Greenberg D, Dagan R. Pneumococcal nasopharyngeal carriage in children <5 years of age visiting the pediatric emergency room in relation to PCV7 and PCV13 introduction in southern Israel. Hum Vaccin Immunother. 2016;12:268–76. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/21645515.2015.1095414&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26430921&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F01%2F2022.09.06.22279606.atom) 26. 26.Cherian T, Mulholland EK, Carlin JB, Ostensen H, Amin R, de Campo M, et al. Standardized interpretation of paediatric chest radiographs for the diagnosis of pneumonia in epidemiological studies. Bull World Health Organ. 2005;83:353–9. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=S0042-96862005000500011&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15976876&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F01%2F2022.09.06.22279606.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000229056700009&link_type=ISI) 27. 27.WHO. Pneumonia Vaccine Trial Investigators’ Group & World Health Organization. Standardization of interpretation of chest radiographs for the diagnosis of pneumonia in children World Health Organization. 2001:[https://apps.who.int/iris/handle/10665/66956](https://apps.who.int/iris/handle/10665/66956). 28. 28.Greenberg D, Givon-Lavi N, Ben-Shimol S, Ziv JB, Dagan R. Impact of PCV7/PCV13 introduction on community-acquired alveolar pneumonia in children <5 years. Vaccine. 2015;33:4623–9. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.vaccine.2015.06.062&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F01%2F2022.09.06.22279606.atom) 29. 29.Ben Shimol S, Dagan R, Givon-Lavi N, Tal A, Aviram M, Bar-Ziv J, et al. Evaluation of the World Health Organization criteria for chest radiographs for pneumonia diagnosis in children. Eur J Pediatr. 2012;171:369–74. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00431-011-1543-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21870077&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F01%2F2022.09.06.22279606.atom) 30. 30.Ben-Shimol S, Greenberg D, Givon-Lavi N, Elias N, Glikman D, Rubinstein U, et al. Rapid reduction in invasive pneumococcal disease after introduction of PCV7 into the National Immunization Plan in Israel. Vaccine. 2012;30:6600–7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.vaccine.2012.08.012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22939907&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F01%2F2022.09.06.22279606.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000310419700014&link_type=ISI) 31. 31.Ben-Shimol S, Greenberg D, Givon-Lavi N, Schlesinger Y, Miron D, Aviner S, et al. Impact of PCV7/PCV13 introduction on invasive pneumococcal disease (IPD) in young children: Comparison between meningitis and non-meningitis IPD. Vaccine. 2016;34:4543–50. 32. 32.Dagan R, Juergens C, Trammel J, Patterson S, Greenberg D, Givon-Lavi N, et al. PCV13-vaccinated children still carrying PCV13 additional serotypes show similar carriage density to a control group of PCV7-vaccinated children. Vaccine. 2017;35:945–50. 33. 33.Ben-Shimol S, Givon-Lavi N, Kotler L, Adriaan van der Beek B, Greenberg D, Dagan R. Post-13-Valent Pneumococcal Conjugate Vaccine Dynamics in Young Children of Serotypes Included in Candidate Extended-Spectrum Conjugate Vaccines. Emerg Infect Dis. 2021;27:150–60. 34. 34.Lieberman D, Lieberman D, Shimoni A, Keren-Naus A, Steinberg R, Shemer-Avni Y. Identification of respiratory viruses in adults: nasopharyngeal versus oropharyngeal sampling. J Clin Microbiol. 2009;47:3439–43. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNtIjtzOjU6InJlc2lkIjtzOjEwOiI0Ny8xMS8zNDM5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMTIvMDEvMjAyMi4wOS4wNi4yMjI3OTYwNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 35. 35.Team RC. A language and environment for statistical computing. R Foundation for Statistical Computing. 2020:[http://www.r-project.org/index.html](http://www.r-project.org/index.html). 36. 36.Zheng Z, Warren JL, Shapiro ED, Pitzer VE, Weinberger DM. Estimated incidence of respiratory hospitalizations attributable to RSV infections across age and socioeconomic groups. Pneumonia (Nathan). 2022;14:6. 37. 37.Plummer M. rjags: Bayesian graphical models using MCMC. R package version 3-14. Google Scholar. 2014. 38. 38.Nation ML, Manna S, Tran HP, Nguyen CD, Vy LTT, Uyen DY, et al. Impact of COVID-19 non-pharmaceutical interventions on pneumococcal carriage prevalence and density in Vietnam. medRxiv. 2022:2022.05.05.22274646. 39. 39.Feikin DR, Scott JA, Gessner BD. Use of vaccines as probes to define disease burden. Lancet. 2014;383:1762–70. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(13)61682-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24553294&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F01%2F2022.09.06.22279606.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000336267300030&link_type=ISI) 40. 40.Howard LM, Edwards KM, Zhu Y, Grijalva CG, Self WH, Jain S, et al. Clinical Features of Human Metapneumovirus-Associated Community-acquired Pneumonia Hospitalizations. Clin Infect Dis. 2021;72:108–17. 41. 41.Jain S, Williams DJ, Arnold SR, Ampofo K, Bramley AM, Reed C, et al. Community-acquired pneumonia requiring hospitalization among U.S. children. N Engl J Med. 2015;372:835–45. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1405870&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25714161&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F01%2F2022.09.06.22279606.atom) 42. 42.Wolf DG, Greenberg D, Kalkstein D, Shemer-Avni Y, Givon-Lavi N, Saleh N, et al. Comparison of human metapneumovirus, respiratory syncytial virus and influenza A virus lower respiratory tract infections in hospitalized young children. Pediatr Infect Dis J. 2006;25:320–4. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/01.inf.0000207395.80657.cf&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16567983&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F01%2F2022.09.06.22279606.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000236920400007&link_type=ISI) 43. 43.Hu Z, Lin J, Chen J, Cai T, Xia L, Liu Y, et al. Overview of Viral Pneumonia Associated With Influenza Virus, Respiratory Syncytial Virus, and Coronavirus, and Therapeutics Based on Natural Products of Medicinal Plants. Front Pharmacol. 2021;12:630834. 44. 44.Ben-Shimol S, Greenberg D, Hazan G, Givon-Lavi N, Gottesman G, Grisaru-Soen G, et al. Differential impact of pneumococcal conjugate vaccines on bacteremic pneumonia versus other invasive pneumococcal disease. Pediatr Infect Dis J. 2015;34:409–16. 45. 45.Singleton RJ, Bulkow LR, Miernyk K, DeByle C, Pruitt L, Hummel KB, et al. Viral respiratory infections in hospitalized and community control children in Alaska. J Med Virol. 2010;82:1282–90. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jmv.21790&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20513097&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F01%2F2022.09.06.22279606.atom) 46. 46.Tapia MD, Sylla M, Driscoll AJ, Toure A, Kourouma N, Sissoko S, et al. The Etiology of Childhood Pneumonia in Mali: Findings From the Pneumonia Etiology Research for Child Health (PERCH) Study. Pediatr Infect Dis J. 2021;40:S18–S28. 47. 47.Yun KW, Wallihan R, Desai A, Alter S, Ambroggio L, Cohen DM, et al. Clinical Characteristics and Etiology of Community-acquired Pneumonia in US Children, 2015-2018. Pediatr Infect Dis J. 2022;41:381–7. [1]: /embed/graphic-2.gif [2]: /embed/graphic-3.gif [3]: /embed/graphic-4.gif [4]: /embed/graphic-5.gif