Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Drivers of respiratory syncytial virus seasonal epidemics in children under 5 years in Kilifi, coastal Kenya

  • James Wambua ,

    Roles Conceptualization, Data curation, Formal analysis, Methodology, Software, Visualization, Writing – original draft, Writing – review & editing

    james.wambua@uhasselt.be

    Current address: Data Science Institute, I-BioStat, Hasselt University, Hasselt, Belgium

    Affiliation Kenya Medical Research Institute, KEMRI -Wellcome Trust Research Programme (KWTRP), Kilifi, Kenya

  • Patrick K. Munywoki,

    Roles Writing – review & editing

    Affiliation Kenya Medical Research Institute, KEMRI -Wellcome Trust Research Programme (KWTRP), Kilifi, Kenya

  • Pietro Coletti,

    Roles Writing – review & editing

    Affiliation Data Science Institute, I-BioStat, Hasselt University, Hasselt, Belgium

  • Bryan O. Nyawanda,

    Roles Writing – review & editing

    Affiliation Kenya Medical Research Institute, Center for Global Health Research, Kisumu, Kenya

  • Nickson Murunga,

    Roles Data curation, Writing – review & editing

    Affiliation Kenya Medical Research Institute, KEMRI -Wellcome Trust Research Programme (KWTRP), Kilifi, Kenya

  • D. James Nokes,

    Roles Conceptualization, Funding acquisition, Supervision, Writing – review & editing

    Affiliations Kenya Medical Research Institute, KEMRI -Wellcome Trust Research Programme (KWTRP), Kilifi, Kenya, School of Life Sciences, University of Warwick, Coventry, United Kingdom

  • Niel Hens

    Roles Conceptualization, Methodology, Supervision, Writing – review & editing

    Affiliations Data Science Institute, I-BioStat, Hasselt University, Hasselt, Belgium, Centre for Health Economics Research and Modelling Infectious Diseases, Vaccine & Infectious Disease Institute, University of Antwerp, Antwerp, Belgium

Abstract

Respiratory syncytial virus (RSV) causes significant childhood morbidity and mortality in the developing world. The determinants of RSV seasonality are of importance in designing interventions. They are poorly understood in tropical and sub-tropical regions in low- and middle-income countries. Our study utilized long-term surveillance data on cases of RSV associated with severe or very severe pneumonia in children aged 1 day to 59 months admitted to the Kilifi County Hospital. A generalized additive model was used to investigate the association between RSV admissions and meteorological variables (maximum temperature, rainfall, absolute humidity); weekly number of births within the catchment population; and school term dates. Furthermore, a time-series-susceptible-infected-recovered (TSIR) model was used to reconstruct an empirical transmission rate which was used as a dependent variable in linear regression and generalized additive models with meteorological variables and school term dates. Maximum temperature, absolute humidity, and weekly number of births were significantly associated with RSV activity in the generalized additive model. Results from the TSIR model indicated that maximum temperature and absolute humidity were significant factors. Rainfall and school term did not yield significant relationships. Our study indicates that meteorological parameters and weekly number of births potentially play a role in the RSV seasonality in this region. More research is required to explore the underlying mechanisms underpinning the observed relationships.

Introduction

Respiratory syncytial virus (RSV) is an ubiquitous RNA virus that is a significant cause of early childhood morbidity and mortality worldwide [1,2]. In low- and middle-income settings, RSV is recognized as an important cause of hospitalized severe pneumonia in children aged between 28 days and 5 years [1]. In 2015, it was estimated that there were 33.1 million (uncertainty range 21.6–50.3) RSV- Acute Lower Respiratory Infections (ALRI) episodes globally, leading to about 3.2 (2.7–3.8) million hospital admissions, and 59,600 (48,000–74,500) hospital deaths, in children < 5 years of age [2]. An approximate 1.4 million (1.2–1.7) in hospital admissions, and 27,300 (20,700–36,200) in-hospital deaths were in young children < 6 months of age. The vast majority of both morbidity and mortality arises from low- and middle-income countries [2].

There is no recognised prophylactic for RSV-ALRI for general use. After over 50 years of vaccine research and development, there is none, to date, that has been licensed [3]. Several RSV vaccines and monoclonal antibody products are at various stages of pre-clinical and clinical development targeting infants and pregnant women [4]. Most recently, a prophylactic high-potency long half-life anti-F protein monoclonal antibody met its primary endpoint in a Phase III trial [5]. Two maternal booster F protein vaccines are in Phase III trials [6]. The successful use of these products in early infants and pregnant women to prevent severe disease in children in the first few months of life will be affected by the timing of seasonal RSV transmission [7].

RSV is characterized by marked annual or biennial seasonal patterns [8]. In temperate climates, RSV infections have been found to happen in the winter months [9,10]. In tropical regions, a wide variation in the timing has been observed [10]. However, comparatively little attention has been directed to most tropical settings in Africa, hence the RSV seasonality pattern remains largely unexplored in these regions, despite the considerable ALRI burden and deaths [10,11]. To understand RSV seasonality, factors that drive these patterns should be well understood. In regions with temperate climates, RSV seasonal epidemics peak in winter months [11]. Previous studies in different tropical regions worldwide have yielded inconsistent results. Temperature, humidity and precipitation have been found to be strongly, weakly or not at all associated with RSV activity [1218]. Seasonality in births has also been shown to impact the seasonality in infectious diseases, especially with respect to first infections [19,20]. Malnutrition has been found to be associated with RSV seasonality [21,22]. Demographic factors such as household crowding indices have been found to influence the length of RSV seasons [23].

Seasonal patterns of RSV in Kenya, through which the equator passes have been compared in 3 distinct counties (Nairobi, Siaya and Kilifi) [24]. However, the potential factors driving these seasonal patterns have not been assessed. The current study was carried out using previously collected datasets from the Kenya Medical Research Institute (KEMRI)-Wellcome Trust Research Programme (KWTRP), located in Kilifi Town along the coastal region bordering the Indian Ocean. The study aimed to investigate whether meteorological parameters (maximum temperature, absolute humidity, rainfall), weekly number of births, and school terms could be determinants of seasonal RSV epidemics. Identifying the potential factors driving RSV seasonality would enhance the understanding of the virus transmission mechanisms in this setting.

Methods

Ethics statement

The ethical approval for the study was granted by the Kenya Medical Research Institute–Scientific Ethics Review Unit (reference number SERU 3178). Written informed consent was obtained from parents or guardians for all eligible participants before any specimen(s) was collected for RSV. The data was analysed anonymously.

Inclusivity in global research

Additional information regarding the ethical, cultural, and scientific considerations specific to inclusivity in global research is included in the Supporting Information (S1 Checklist).

Data

Data used were collected from routine longitudinal hospital case surveillance for RSV, and routine longitudinal surveillance for the number of births and meteorological variables. The datasets are described below.

RSV data

RSV data was obtained from Kilifi County Hospital located in Kilifi Town. Data was collected from children aged 1 day to 59 months who presented with a clinical syndrome for either severe or very severe pneumonia [25]. Written informed consent was obtained from parents or guardians for all eligible participants before any specimen(s) was collected. The data spanned the years 2002 to 2018. RSV was screened through immunofluorescence antibody test (IFAT) for all years and real time reverse transcription-polymerase chain reaction (RT-rtPCR) from 2007 [25]. Only data from locations covered by the Kilifi Health and Demographic Surveillance System (KHDSS) were included [26]. The KHDSS area extends roughly 40km north and south of Kilifi town and 30km in land, largely rural agricultural and the population size has approximately increased from 200,000 to 300,000 over the study period (Fig 1).

thumbnail
Fig 1. Geographical location of the study area.

Left: Location of Kilifi county in Kenya. Middle: Map of Kilifi county showing KHDSS boundaries. Right: KHDSS map showing Kilifi County Hospital location. Reprinted from [27] under a CC BY license, with permission from BioMed Central, Alice Kamau (original copyright owner) [2017].

https://doi.org/10.1371/journal.pone.0278066.g001

Births data

Weekly number of birth data was obtained from the KHDSS registers. The system records the number of births, deaths, pregnancies and migration events. The data spanned the years 2002 to 2018.

Meteorological data

Meteorological data was obtained from two weather stations. One station at Pwani University, Kilifi Town and the other at KWTRP. Data from the Pwani station was recorded manually and consisted of daily measurements of four main variables; maximum/minimum temperature (in degree Celsius), rainfall (in millimeters), and average relative humidity (expressed in %). The data comprised of measurements from January 2002 to July 2014. Data from KWTRP station for maximum/minimum temperature, rainfall, and average relative humidity was measured by an automatic weather station (Synoptic Automatic Weather Stations, Sutron, USA). Data from this station was collected from July 2010 to November 2018. Absolute humidity was obtained using equation 1, see supplementary material.

School term data

The data consisted of school start dates and end dates of half terms and holidays for pre-primary and primary [28]. Children attend pre-primary schooling for 2–3 years and encompasses the age group 3–5 years. Primary school is attended for 8 years and includes children 6 years and above.

Data processing

Weekly averages of meteorological variables were used in the analysis. Due to missing data in the two meteorological datasets, multiple imputation using multivariate imputation by chained equations (MICE package) implemented in the R Statistical software (version 4.1.1) was performed on the daily observations [29]. We performed five imputations resulting in five datasets for each meteorological dataset. The number of imputations (5) was chosen based on imputation guidelines by Rubin [30] on the percentages of missing data. The resulting datasets were merged with the RSV data, birth data, and data on school term by year and calendar week. We performed separate analyses for the two meteorological datasets. See supplementary material for more details on the rationale for separate analyses of the 2002–2014 and 2010–2018 datasets.

Statistical analysis

To explore the plausible drivers of RSV seasonal epidemics, we applied two analytical approaches. The first involved using a Generalized Additive Model (GAM), whilst the second involved a time-series-susceptible-infected-recovered model (TSIR). GAMs are generalized linear models with predictors specified as smooth functions of the covariates [31]. The cubic spline smoothing function was adopted as it has been shown to yield smoother interpolations when compared to other methods [31,32]. The model employed a negative binomial with a log-link function due to over-dispersion in the weekly number of RSV cases. The model included meteorological variables, weekly number of births, school term (in or out of term time), and RSV activity in the previous 3 weeks (due to the serial correlation resulting from the infection dynamics of RSV). We explored models for the lagged variables up to the previous one week (i.e. values observed in the previous week for the covariates) to explore delay effects. A cubic smoothing spline function of time was included to control for seasonal patterns and long-term trend. Statistical significance was determined by an alpha level of 5%. All analyses were done using the mgcv package in the R statistical software [33].

The TSIR method [3436] was applied to enhance the exploration of how changes in birth rates affect RSV dynamics as it has been shown for measles epidemics in England and Wales [36]. Furthermore, it also provided an alternative approach to explore the plausible drivers of RSV transmission dynamics in this region [34]. This discrete time epidemic model comprises two state variables, namely: the infected and the susceptible. Since the susceptible is unobserved, we reconstructed it using a set of recursive difference equations utilizing the weekly number of RSV cases, weekly number of births and interpolated total weekly population over the years using the tsiR package implemented in R statistical software [35]. Next, using the reconstructed susceptible population over time, we reconstructed an empirical estimate of the transmission rate which was subsequently used as a dependent variable in both linear regression and generalized additive negative binomial models with meteorological variables and school term data as independent variables. More information on the TSIR and GAM methodologies is contained in the supplementary material.

Both analyses were performed on each of the five imputed 2002–2014 datasets and the five imputed 2010–2018 datasets obtained as described above in the data pre-processing sub-section. Our choice of independent analysis was motivated by the fact that imputation of the missing data in the meteorological variables was performed at daily resolutions and the weekly averages were obtained to merge with the other considered datasets. Robust imputation methodologies for data collected in different resolutions are still lacking in the literature and thus future work in this area is warranted.

Results

In all our analyses, we observed that the weekly number of births was only significant at lag two and three (number of births observed in the previous two and three weeks, respectively). Hence in the GAM models, we considered a lag of two weeks of the weekly number of births to coincide with the meteorological variables and school term, and a lag of three weeks of the weekly births to coincide with a lag of one week of the meteorological variables and school term. The two approaches (GAM and TSIR) yielded similar results, hence, we report the GAM analyses below, the TSIR results are in the supplementary material.

In the 2002–2014 datasets, we observed that maximum temperature both at lags zero and one week was a significant predictor of the weekly number of RSV cases (p-value <0.05, Table 1). Absolute humidity was a significant predictor at lag one week (p-value <0.05, Table 1). The weekly number of births was a significant predictor of the weekly number of RSV cases at both lags two and three weeks (p-value <0.05, Table 1). Rainfall and school terms did not yield significant relationships. The plots of the smooth terms of maximum temperature, absolute humidity and weekly number of births indicated non-linear relationships with the weekly number of RSV cases (S1 Fig).

thumbnail
Table 1. Generalized additive negative binomial models with lags zero and one week of meteorological variables and lags two and three weeks of weekly number of births, respectively, in the 2002–2014 datasets.

https://doi.org/10.1371/journal.pone.0278066.t001

Analyses using the 2010–2018 datasets yielded similar results where maximum temperature and absolute humidity were significant predictors of the weekly RSV cases at both lags zero and one week (p-value <0.05, S1 Table in S1 File). Similarly, weekly number of births at both lags two and three weeks was a significant predictor of the RSV cases (S1 Table in S1 File). The plots of the smooth terms of maximum temperature, absolute humidity, and weekly number of births showed non-linear relationships with the weekly number of RSV cases (Fig 2).

thumbnail
Fig 2. Plots of the meteorological smooth terms in the first imputed 2010–2018 dataset.

The dashed lines are the 95% confidence intervals.

https://doi.org/10.1371/journal.pone.0278066.g002

Generally, the correlation coefficients between the observed and the predicted number of weekly RSV cases was high for both the 2002–2014 datasets (Table 1), and the 2010–2018 datasets (S1 Table in S1 File). The curves of the predicted and the observed weekly RSV cases closely followed each other for both the 2002–2014 datasets (Fig 3), and in the 2010–2018 datasets (Fig 4). The percentage of deviance explained was relatively higher for models including the lagged values of meteorological variables, weekly number of births and school terms (Table 1 & S1 Table in S1 File). The 2010–2018 datasets yielded higher values of the percentages of the deviance explained, the adjusted R-squared, and the correlation coefficients between the observed and predicted weekly number of RSV cases (Table 1 & S1 Table in S1 File).

thumbnail
Fig 3. Plot of observed and predicted weekly RSV cases in the first imputed 2002–2014 dataset.

The black line connects the observed weekly RSV cases while the red line connects the predicted weekly RSV cases.

https://doi.org/10.1371/journal.pone.0278066.g003

thumbnail
Fig 4. Plot of observed and predicted weekly RSV cases in the first imputed 2010–2018 dataset.

The black line connects the observed weekly RSV cases while the red line connects the predicted weekly RSV cases.

https://doi.org/10.1371/journal.pone.0278066.g004

Discussion

In this study, we found that maximum temperature, absolute humidity and weekly number of births are significant predictors of the weekly number of RSV cases from the independent analysis of ten imputed datasets.

Although humans can reproduce throughout the year, some human populations exhibit seasonal variation in reproduction leading to seasonal birth rates/patterns [37]. These variations can be utilized to characterize the seasonal variations in the incidence of infectious diseases in early childhood [20,37]. In our study, the weekly number of births was a significant predictor of RSV. Seasonality of births leads to pulses of new hosts in the population leading to an expansion of the susceptible population, a significant factor modulating the transmission dynamics of childhood infectious diseases [19,20]. A study in the US [38] found out that year-to-year and state-to-state variations in birth rates determined the timing of rotavirus epidemics. Another study in the US [39] showed variations in RSV epidemic timing where counties with larger population sizes had earlier epidemic peaks and recorded higher RSV incidence. A mathematical modelling study [40] found that changes in birth rates led to either annual or biennial RSV seasonal epidemics. Another modelling study [16] in the US showed that the transition to annual seasonal RSV epidemics in 2000 from biennial patterns in 1990s in California could have been modulated by the changes in the birth rates. Similarly, another study [20] in Sub-Saharan Africa showed that changes in seasonal birth rates can either lead to annual or biennial dynamics of infectious diseases occurring early in life such as measles. Thus, the observed associations between the weekly number of births and weekly RSV cases in this region could suggest that the seasonal changes in births could be a potential driver of the observed RSV seasonal epidemics. Further research is needed in this area.

School term did not yield significant relationships in our study. This might suggest that school terms do not play a major role in determining the RSV seasonality in this region. However, aggregation of children in schools has been found to play an important role in increasing the transmission rates of measles in England and Wales [41], Influenza in England and varicella zoster virus in several settings [42]. For RSV, it has been hypothesized that increased contact rates during school terms could play a crucial role in influencing its seasonality [43]. In addition, several studies have demonstrated that RSV is mostly introduced into the households by school going siblings consequently leading to the infections to the infants [44,45]. However, results from 4 recent studies do not appear to support this earlier hypothesis [1416,34] consistent with what we observed in our study. Insights into why ‘school terms’ is not identified as a significant contributor to RSV seasonality would be gained from further studies of RSV transmission patterns in the school setting, including linking RSV individual-level infections data and social contact patterns for school going children at school-level, community and household levels to gain deeper understanding of ‘who infects who’ in the community.

Rainfall did not yield significant relationships in our study implying that it might not be a plausible driver of RSV seasonal epidemics in this region. This is in contrast with several other studies that have found associations between rainfall and RSV activity in the following regions: Baguio city in the Philippines, Okinawa in Japan [46], tropical Kolkata in India [47], and some states in the US [16], Lombok island in Indonesia [12], Kuala Lumpur in Malaysia [47], the Philippines [15], Gambia in West Africa [17], and an Equatorial city in Brazil [48]. Maximum temperature showed non-linear relationship with RSV activity in this region. Interestingly, absolute humidity also exhibited a non-linear relationship with the RSV activity. A previous laboratory study [49] conducted to understand the stability of RSV in varying conditions found out that at 37°C, RSV inactivation took place over a period of 3 days. Varying the temperature to 4°C, indicated that 99% inactivation required a duration of 6 days or more. Another study [50] at a constant temperature of 20.5°C, the maximal stability of RSV was recorded at 80–90% humidity, and maximal inactivation occurred at 20–30% humidity. Thus, we postulate that stability of RSV in large particle aerosols in this tropical setting could be dependent on variations in maximum temperature and absolute humidity. This could presumably modulate aerosol transmission through altering the stability of RSV in the environment.

While meteorological factors and weekly number of births could likely influence RSV seasonality, other plausible drivers of RSV epidemics such as malnutrition [22], household crowding and population density [23] could be explored. In our study, we used outdoor measurements for the meteorological variables which could differ with indoor conditions where the study domain (children under 5 years) spend much of their time. School term data was inferred from the 2019 academic calendar, there might be year to year variations over the study period. Our study assumed that the meteorological data had the same spatial coverage in all locations where the study domain lived and there was no spatial heterogeneity based on the distribution of these RSV cases. Furthermore, although Kenya sits on the equator, it has quite different climatic regions. Thus, the chosen study location of coastal Kenya may differ from other regions and hence our results might not be generalizable. More studies, using data from geographically distinct locations with variation in the seasonal patterns and potential drivers are needed in this regard.

Conclusion

Our study has identified maximum temperature, absolute humidity and weekly number of births as potential drivers of RSV epidemics in this region. These findings contribute to the current understanding of the potential drivers of RSV seasonality in a tropical low-income setting where little is known about the plausible drivers of the RSV seasonality. Identifying these factors is important in improving the understanding of the transmission dynamics of the virus in this setting. Furthermore, these factors can be utilized in the formulation of prediction models of RSV activity in this setting. This potentially could provide early warnings of RSV seasons and subsequently inform the timings of RSV immunizations using vaccine candidates currently in the pipeline. However, more research is needed in other regions with different climatic classifications to shed more light on the observed relationships.

Supporting information

S1 Checklist. Completed questionnaire on inclusivity in global research.

https://doi.org/10.1371/journal.pone.0278066.s001

(DOCX)

S1 Fig. Plots of the meteorological smooth terms in the first imputed 2002–2014 dataset.

The dashed lines are the 95% confidence intervals.

https://doi.org/10.1371/journal.pone.0278066.s002

(TIF)

S2 Fig. The estimated transmission rate from the time-series-susceptible-infected-recovered (TSIR) model with the 95% confidence interval.

https://doi.org/10.1371/journal.pone.0278066.s003

(TIF)

S3 Fig. Plot of the mean of 100 stochastic simulations from the TSIR model fit (red) plotted alongside observed data of weekly RSV cases (blue).

https://doi.org/10.1371/journal.pone.0278066.s004

(TIF)

S4 Fig. Plot for the weekly number of RSV cases and the weekly number of births from the year 2002 to 2018.

https://doi.org/10.1371/journal.pone.0278066.s005

(TIF)

S5 Fig. Weekly number of RSV cases compared to mean weekly maximum temperature, absolute humidity, rainfall, and potential evapotranspiration between the year 2002 and 2014.

https://doi.org/10.1371/journal.pone.0278066.s006

(TIF)

S6 Fig. Weekly number of RSV cases compared to mean weekly maximum temperature, absolute humidity, rainfall, and potential evapotranspiration between the year 2010–2018.

https://doi.org/10.1371/journal.pone.0278066.s007

(TIF)

S1 File. Details on data pre-processing and Time Series Susceptible-Infected-Recovered (TSIR) methodology, and supplementary results.

https://doi.org/10.1371/journal.pone.0278066.s008

(DOCX)

Acknowledgments

This study is published with permission from the Director of KEMRI.

References

  1. 1. O’Brien KL, Baggett HC, Brooks WA, Feikin DR, Hammitt LL, Higdon MM, et al. Causes of severe pneumonia requiring hospital admission in children without HIV infection from Africa and Asia: the PERCH multi-country case-control study. The Lancet. 2019;394(10200):757–79. pmid:31257127
  2. 2. Shi T, McAllister DA, O’Brien KL, Simoes EAF, Madhi SA, Gessner BD, et al. Global, regional, and national disease burden estimates of acute lower respiratory infections due to respiratory syncytial virus in young children in 2015: a systematic review and modelling study. The Lancet. 2017;390(10098):946–58.
  3. 3. Simões EA, Bont L, Manzoni P, Fauroux B, Paes B, Figueras-Aloy J, et al. Past, present and future approaches to the prevention and treatment of respiratory syncytial virus infection in children. Infectious diseases and therapy. 2018;7(1):87–120. pmid:29470837
  4. 4. Mejias A, Rodríguez-Fernández R, Oliva S, Peeples ME, Ramilo O. The journey to a respiratory syncytial virus vaccine. Ann Allergy Asthma Immunol. 2020;125(1):36–46. pmid:32217187
  5. 5. Kemp A. Nirsevimab MELODY Phase III trial met primary endpoint of reducing RSV lower respiratory tract infections in healthy infants.
  6. 6. Gunatilaka A, Giles ML. Maternal RSV vaccine development. Where to from here? Human Vaccines & Immunotherapeutics. 2021;17(11):4542–8.
  7. 7. Li Y, Hodgson D, Wang X, Atkins KE, Feikin DR, Nair H. Respiratory syncytial virus seasonality and prevention strategy planning for passive immunisation of infants in low-income and middle-income countries: a modelling study. The Lancet Infectious Diseases. 2021;21(9):1303–12. pmid:33965062
  8. 8. White LJ, Waris M, Cane PA, Nokes DJ, Medley GF. The transmission dynamics of groups A and B human respiratory syncytial virus (hRSV) in England & Wales and Finland: seasonality and cross-protection. Epidemiol Infect. 2005;133(2):279–89.
  9. 9. Bloom-Feshbach K, Alonso WJ, Charu V, Tamerius J, Simonsen L, Miller MA, et al. Latitudinal variations in seasonal activity of influenza and respiratory syncytial virus (RSV): a global comparative review. PLoS One. 2013;8(2):e54445. pmid:23457451
  10. 10. Li Y, Reeves RM, Wang X, Bassat Q, Brooks WA, Cohen C, et al. Global patterns in monthly activity of influenza virus, respiratory syncytial virus, parainfluenza virus, and metapneumovirus: a systematic analysis. The Lancet Global Health. 2019;7(8):e1031–e45. pmid:31303294
  11. 11. Obando-Pacheco P, Justicia-Grande AJ, Rivero-Calle I, Rodríguez-Tenreiro C, Sly P, Ramilo O, et al. Respiratory syncytial virus seasonality: a global overview. The Journal of infectious diseases. 2018;217(9):1356–64. pmid:29390105
  12. 12. Omer SB, Sutanto A, Sarwo H, Linehan M, Djelantik IG, Mercer D, et al. Climatic, temporal, and geographic characteristics of respiratory syncytial virus disease in a tropical island population. Epidemiol Infect. 2008;136(10):1319–27. pmid:18177515
  13. 13. Paynter S. Humidity and respiratory virus transmission in tropical and temperate settings. Epidemiol Infect. 2015;143(6):1110–8. pmid:25307020
  14. 14. Paynter S, Ware RS, Sly PD, Weinstein P, Williams G. Respiratory syncytial virus seasonality in tropical Australia. Aust N Z J Public Health. 2015;39(1):8–10. pmid:25648729
  15. 15. Paynter S, Yakob L, Simões EAF, Lucero MG, Tallo V, Nohynek H, et al. Using Mathematical Transmission Modelling to Investigate Drivers of Respiratory Syncytial Virus Seasonality in Children in the Philippines. PLOS ONE. 2014;9(2):e90094. pmid:24587222
  16. 16. Pitzer VE, Viboud C, Alonso WJ, Wilcox T, Metcalf CJ, Steiner CA, et al. Environmental Drivers of the Spatiotemporal Dynamics of Respiratory Syncytial Virus in the United States. PLOS Pathogens. 2015;11(1):e1004591. pmid:25569275
  17. 17. Weber MW, Mulholland EK, Greenwood BM. Respiratory syncytial virus infection in tropical and developing countries. Trop Med Int Health. 1998;3(4):268–80. pmid:9623927
  18. 18. Yusuf S, Piedimonte G, Auais A, Demmler G, Krishnan S, Van Caeseele P, et al. The relationship of meteorological conditions to the epidemic activity of respiratory syncytial virus. Epidemiol Infect. 2007;135(7):1077–90. pmid:17346359
  19. 19. Altizer S, Dobson A, Hosseini P, Hudson P, Pascual M, Rohani P. Seasonality and the dynamics of infectious diseases. Ecol Lett. 2006;9(4):467–84. pmid:16623732
  20. 20. Dorélien AM, Ballesteros S, Grenfell BT. Impact of Birth Seasonality on Dynamics of Acute Immunizing Infections in Sub-Saharan Africa. PLOS ONE. 2013;8(10):e75806. pmid:24204580
  21. 21. Paynter S, Ware RS, Lucero MG, Tallo V, Nohynek H, Simões EA, et al. Poor growth and pneumonia seasonality in infants in the Philippines: cohort and time series studies. PLoS One. 2013;8(6):e67528. pmid:23840731
  22. 22. Paynter S, Ware RS, Lucero MG, Tallo V, Nohynek H, Weinstein P, et al. Malnutrition: a risk factor for severe respiratory syncytial virus infection and hospitalization. The Pediatric infectious disease journal. 2014;33(3):267–71. pmid:24168980
  23. 23. Zachariah P, Shah S, Gao D, Simões EA. Predictors of the duration of the respiratory syncytial virus season. Pediatr Infect Dis J. 2009;28(9):772–6. pmid:19710584
  24. 24. Rose EB, Nyawanda BO, Munywoki PK, Murunga N, Bigogo GM, Otieno NA, et al. Respiratory syncytial virus seasonality in three epidemiological zones of Kenya. Influenza and Other Respiratory Viruses. 2021;15(2):195–201. pmid:33305543
  25. 25. Nokes DJ, Ngama M, Bett A, Abwao J, Munywoki P, English M, et al. Incidence and severity of respiratory syncytial virus pneumonia in rural Kenyan children identified through hospital surveillance. Clinical Infectious Diseases. 2009;49(9):1341–9. pmid:19788358
  26. 26. Scott JA, Bauni E, Moisi JC, Ojal J, Gatakaa H, Nyundo C, et al. Profile: The Kilifi Health and Demographic Surveillance System (KHDSS). Int J Epidemiol. 2012;41(3):650–7. pmid:22544844
  27. 27. Kamau A, Nyaga V, Bauni E, Tsofa B, Noor AM, Bejon P, et al. Trends in bednet ownership and usage, and the effect of bednets on malaria hospitalization in the Kilifi Health and Demographic Surveillance System (KHDSS): 2008–2015. BMC infectious diseases. 2017;17(1):1–9.
  28. 28. Kenya Term Dates and School Holidays. 2019.
  29. 29. Van Buuren S, Groothuis-Oudshoorn K. mice: Multivariate imputation by chained equations in R. Journal of statistical software. 2011;45:1–67.
  30. 30. Rubin DB, editor An overview of multiple imputation. Proceedings of the survey research methods section of the American statistical association; 1988: Citeseer.
  31. 31. Wood SN. Generalized additive models: an introduction with R: CRC press; 2017.
  32. 32. Azizan I, Karim SABA, Raju SSK, editors. Fitting rainfall data by using cubic spline interpolation. MATEC Web of Conferences; 2018: EDP Sciences.
  33. 33. Wood S. Mixed GAM computation vehicle with automatic smoothness estimation. R package version 1.8–12. 2018.
  34. 34. Baker RE, Mahmud AS, Wagner CE, Yang W, Pitzer VE, Viboud C, et al. Epidemic dynamics of respiratory syncytial virus in current and future climates. Nature communications. 2019;10(1):1–8.
  35. 35. Becker AD, Grenfell BT. tsiR: An R package for time-series Susceptible-Infected-Recovered models of epidemics. PloS one. 2017;12(9):e0185528. pmid:28957408
  36. 36. Finkenstädt BF, Grenfell BT. Time series modelling of childhood diseases: a dynamical systems approach. Journal of the Royal Statistical Society: Series C (Applied Statistics). 2000;49(2):187–205.
  37. 37. Dorélien A. Birth seasonality in Sub-Saharan Africa. Demographic Research. 2016;34(27):761–96.
  38. 38. Pitzer VE, Viboud C, Simonsen L, Steiner C, Panozzo CA, Alonso WJ, et al. Demographic variability, vaccination, and the spatiotemporal dynamics of rotavirus epidemics. Science. 2009;325(5938):290–4. pmid:19608910
  39. 39. Noveroske DB, Warren JL, Pitzer VE, Weinberger DM. Local variations in the timing of RSV epidemics. BMC Infect Dis. 2016;16(1):674. pmid:27835988
  40. 40. Hogan AB, Glass K, Moore HC, Anderssen RS. Exploring the dynamics of respiratory syncytial virus (RSV) transmission in children. Theor Popul Biol. 2016;110:78–85. pmid:27155294
  41. 41. Fine PE, Clarkson JA. Measles in England and Wales—I: An analysis of factors underlying seasonal patterns. Int J Epidemiol. 1982;11(1):5–14. pmid:7085179
  42. 42. Jackson C, Mangtani P, Fine P, Vynnycky E. The Effects of School Holidays on Transmission of Varicella Zoster Virus, England and Wales, 1967–2008. PLOS ONE. 2014;9(6):e99762. pmid:24932994
  43. 43. Waris M, White LJ. Seasonality of respiratory syncytial virus infection. Clinical infectious diseases: an official publication of the Infectious Diseases Society of America. 2006;43(4):541. pmid:16838247
  44. 44. Hall CB, Geiman JM, Biggar R, Kotok DI, Hogan PM, Douglas GR Jr. Respiratory syncytial virus infections within families. N Engl J Med. 1976;294(8):414–9. pmid:173995
  45. 45. Munywoki PK, Koech DC, Agoti CN, Lewa C, Cane PA, Medley GF, et al. The source of respiratory syncytial virus infection in infants: a household cohort study in rural Kenya. The Journal of infectious diseases. 2014;209(11):1685–92. pmid:24367040
  46. 46. Kamigaki T, Chaw L, Tan AG, Tamaki R, Alday PP, Javier JB, et al. Seasonality of Influenza and Respiratory Syncytial Viruses and the Effect of Climate Factors in Subtropical–Tropical Asia Using Influenza-Like Illness Surveillance Data, 2010–2012. PLOS ONE. 2016;11(12):e0167712. pmid:28002419
  47. 47. Khor C-S, Sam I-C, Hooi P-S, Quek K-F, Chan Y-F. Epidemiology and seasonality of respiratory viral infections in hospitalized children in Kuala Lumpur, Malaysia: a retrospective study of 27 years. BMC pediatrics. 2012;12(1):1–9.
  48. 48. Alonso WJ, Laranjeira BJ, Pereira SA, Florencio CM, Moreno EC, Miller MA, et al. Comparative dynamics, morbidity and mortality burden of pediatric viral respiratory infections in an equatorial city. Pediatr Infect Dis J. 2012;31(1):e9–14. pmid:22001966
  49. 49. Hambling MH. SURVIVAL OF THE RESPIRATORY SYNCYTIAL VIRUS DURING STORAGE UNDER VARIOUS CONDITIONS. Br J Exp Pathol. 1964;45(6):647–55. pmid:14245166
  50. 50. Rechsteiner J. Inactivation of respiratory syncytial virus in air. Antonie Van Leeuwenhoek. 1969;35(2):238. pmid:5310454