
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39223167
70796
10.1038/s41598-024-70796-7
Article
Respiratory rates among rural Gambian children: a community-based cohort study
Mogeni Polycarp pkambona11@gmail.com

12
Amima Sharon 3
Gunther Jennifer 1
Pinder Margaret 45
Tusting Lucy S. 6
D’Alessandro Umberto 5
Cousens Simon 1
Lindsay Steve W. s.w.lindsay@durham.ac.uk

4
Bradley John 1
1 https://ror.org/00a0jsq62 grid.8991.9 0000 0004 0425 469X Department of Infectious Disease Epidemiology, London School of Hygiene and Tropical Medicine, London, UK
2 https://ror.org/04r1cxt79 grid.33058.3d 0000 0001 0155 5938 Kenya Medical Research Institute (KEMRI), Nairobi, Kenya
3 https://ror.org/02y9nww90 grid.10604.33 0000 0001 2019 0495 Department of Food Science, Nutrition and Technology, University of Nairobi, Nairobi, Kenya
4 https://ror.org/01v29qb04 grid.8250.f 0000 0000 8700 0572 Department of Biosciences, Durham University, Durham, UK
5 grid.415063.5 0000 0004 0606 294X Medical Research Council’s (MRC) Unit The Gambia at the London School of Hygiene and Tropical Medicine, Banjul, The Gambia
6 https://ror.org/00a0jsq62 grid.8991.9 0000 0004 0425 469X Department of Disease Control, London School of Hygiene and Tropical Medicine, London, UK
2 9 2024
2 9 2024
2024
14 2035427 2 2024
21 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Normal respiratory rates (RR) for children under five in the tropics are well-documented, but data for older children are limited. This study tracked RR changes with age and examined associations with nutritional status and environmental factors. We monitored rural Gambian children aged 6 months to 14 years, recording RR during home visits twice weekly over two rainy seasons. Using a generalized additive model, we constructed RR reference curves, and a linear mixed-effect model identified factors influencing RR. A total of 830 children provided 67,512 RR measurements. Their median age was 6.07 years (interquartile range 4.21–8.55) and 400 (48.2%) were female. Age, stunting, ambient temperature, and time of RR measurement were independent predictors of respiratory rate. Strikingly, children showing signs of illness had greater variability in repeat RR measurements. We constructed a RR reference chart for children aged one to 13 years and proposed a cutoff of > 26 breaths/min for raised RR among children aged > 5 years bridging an important gap in this age group. Although the time of data collection, nutritional status, and ambient temperature were predictors of RR, their effect size is not clinically significant enough to warrant a change in the current WHO guidelines owing to the prevailing uncertainty in the measurement of RR. The finding that RRs between repeat measurements were more variable among children with signs of illness suggests that a single RR measurement may be inadequate to reliably assess the status of sick children—a population in which accurate diagnosis is essential to enable targeted interventions with lifesaving treatment.

Keywords

Respiratory rate
Centile charts
Nutrition
Ambient temperature
Subject terms

Health care
Medical research
Risk factors
Signs and symptoms
MRC-DfID-Wellcome TrustMR/M007383/1 Global Challenges Research Fund’s BOVA NetworkBB/R00532X/1 Lindsay Steve W. MRC/FCDOMR/R010161/1 Bradley John issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcBackground

Childhood respiratory illness constitutes the largest cause of post-neonatal mortality globally, a burden that could be considerably reduced by timely diagnosis and treatment1–3. Respiratory rate (RR) is an important vital sign monitored in clinical settings for clinical deterioration and respiratory distress and is used to support the diagnosis of severe respiratory infections4–6. As an indicator of respiratory function, healthcare professionals use RR to provide timely intervention, manage disease, and monitor treatment effect6.

The World Health Organization (WHO) defines tachypnoea (abnormally rapid breathing) as a RR ≥ 60 breaths per minute among children aged < 2 months, RR ≥ 50 among children aged between 2 and < 12 months, and RR ≥ 40 among children aged 12 months to < 60 months4. However, RR measurement in children is often challenging7,8 and although various pediatric early warning scores incorporating RR exist, such as the Pediatric Early Warning Score, there are no clear guidelines specifically tailored for children older than 5 years residing in hot, humid, and low-altitude areas in sub-Saharan Africa (sSA). RRs are typically counted manually in low-resource settings as breaths per minute using electronic timers and counting beads9. Though manual measurement is often the gold standard in rural sSA, it can be imprecise and subject to intra-observer variability compounded by whether the child is asleep or awake at the time when measurements are taken8,10. In response to calls for better pneumonia diagnosis11,12, automated RR counting aids are increasingly under development and evaluation13–16. However, there are concerns that some of these aids are out of reach or are unsuitable in resource-limited settings8.

RR measurements are subject to variations associated with both intrinsic factors (e.g. child age, gender, underlying health conditions, child state (asleep, awake or agitated), underlying health conditions) and extrinsic factors that include environmental and socioeconomic factors8,17. Younger children have higher and more variable RR measurements compared to older children whilst those residing at higher altitudes have higher RR measurements, likely due to lower oxygen concentration at higher altitudes17–19.

Here, we report the results of 92 consecutive rounds of community follow-up of children recruited to participate in a trial investigating the effect of improved housing on the incidence of clinical malaria in The Gambia20. The study was undertaken over two consecutive malaria transmission seasons between June and December of 2016 and 201720. We utilized a large longitudinal dataset of RR measurements taken twice weekly to construct the RR-for-age reference curves for children aged 1–13 years and provide a comprehensive description of the predictors of RR among children residing in rural Gambia; a country at low altitude. This is one of the few studies to report RR in children aged > 5 years and the first to do so among rural Gambian children.

Methods

Study setting

The study was conducted in villages in the Upper River Region of The Gambia, located in the east of the country. This is an area of Sudanian savanna at an altitude of less than 50 m above sea level. The mean annual rainfall is 876 mm, occurring mainly between May and October, followed by a long dry season. The enrolled villages are located on both the north and south banks of the River Gambia; an area that experiences an average maximum daily temperature of 36 °C and a minimum of 18 °C. Further details have been reported previously21.

Study design and participants

Data came from a household randomized controlled trial involving children aged 6 months to 14 years, who were recruited to a study investigating the effect of improved housing on the incidence of clinical malaria20,21. The study took place over a 2-year period in which clinical follow-ups were conducted over two consecutive malaria transmission seasons (June to December of 2016 and 2017)21. Details of the study design, data collection, consenting, supervision, and training procedures have been presented previously in a protocol21 and the main study results20. Briefly, alongside the primary trial endpoint of clinical malaria, RR data were recorded twice weekly on every participating child by trained healthcare providers equipped with electronic timers. The data were collected because of concerns that the housing intervention could cause respiratory illness due to restricting airflow through the house. The final trial results, however, showed no evidence of an increase in respiratory illness in the intervention group20. The assessment was done twice weekly21 to increase the case detection rate given the high burden of childhood pneumonia in West Africa22. The process involved the development of a standardized protocol containing guidelines for collecting all clinical measurements, including recording respiratory rate measurements when the child was calm. If any measurement was ≥ 50 breaths/min for children under 1 year, ≥ 40 breaths/min for children aged 1–5 years, or ≥ 25 breaths/min among children aged 5–14 years of age, a repeat measurement was taken at least 5 min after the initial measurement by the same staff member20. Where repeat measurements were taken, we used a second (repeat) measurement for our main analyses. Measurements were taken by observing the child’s naked chest movement, where one rise and fall of the chest was counted as one breath and recorded as the number of breaths per minute.

Statistical analysis

Construction of reference curves

The main analytical objective was to construct a RR reference chart that can be used to inform clinical management of respiratory illnesses in low-altitude areas. The outcome of interest was the RR measurement derived from children who were not reported to be sick during data collection in the household. Here, we excluded data from unwell children, defined by signs of respiratory illness (cough, chest indrawing, wheezing or difficulty in breathing), fever (axillary temperature ≥ 37.5 °C), signs of gastrointestinal illness (vomiting and diarrhea), child reported as being unwell by caregiver or on medication (Fig. 1). In addition, we excluded an implausible outlier (recorded as 181 breaths/min). This measurement was not followed by a second measurement as required by the protocol21.Fig. 1 Number of participants, observations contributed, and reasons for exclusion at each analysis stage.

We conducted an exploratory analysis to examine potential outliers in the dataset and visually inspected the association between respiratory rate and age using the lowess function and scatter plots. Smooth percentile curves of respiratory rate as a function of age were constructed using the generalized additive models for location, scale, and shape (GAMLSS)23. The Box–Cox-t model24 was the best fit for our data among the plausible GAMLSS models (the Box–Cox-t, Box–Cox-Power-Exponential, and Box–Cox–Cole–green) that were examined in a data-driven model selection procedure25,26. The selected model was used to construct the reference chart (at − 2 standard deviation (SD), − 1SD, median, 1SD, and 2SD) of RR as a function of age. The results were truncated at 1 and 13 years to reduce edge effects and the impact of the limited number of study participants < 1 (n = 8) and > 13 years (n = 32) of age.

Predictors of respiratory rate

The second analytical objective was to examine the predictors of respiratory rate to identify potential subgroups that may influence the accuracy of the constructed reference chart. We considered nutritional status, rainfall, minimum daily temperature as a proxy for ambient temperature, time of day when the measurement was taken, sex, and the randomization arm as potential predictors in the regression models. We linked the respiratory rate data to three cross-sectional anthropometric surveys conducted in December 2016, June 2017, and December 2017. The survey datasets were linked to the child’s corresponding respiratory rate observations within a month of each survey to reduce potential bias resulting from linking observations that are far apart in time. Anthropometric z-scores were calculated using the 2006 and 2007 WHO growth references for children aged < 5 years25 and children aged ≥ 5-years26 respectively. We defined stunting as a height-for-age z-score (HAZ) < − 2, underweight as weight-for-age z-score (WAZ) < − 2 and wasting as weight-for-height/length z-score (WHZ) < − 2. Daily rainfall and minimum temperatures were obtained from the meteorological department in Basse Santa Su town. We used univariable analyses (retaining covariates at p value < 0.1) and backward exclusion of non-significant covariates (p > 0.05) to arrive at a final parsimonious model. Age, sex, and randomization variables were included in multivariable models regardless of their significance in the univariable analyses. The variance inflation factor was used to examine potential collinearity.

Associations between respiratory rate and potential predictor variables were assessed using mixed-effects linear regression with random effects for each child. Age was categorized following the trend observed in the lowess curve to accommodate the non-linear trend. The likelihood ratio test statistic was used to compare a model with random intercepts only and a model containing both random intercepts and slopes. Effect estimates and their corresponding 95% confidence intervals (95% CI) were used to assess the magnitude and direction of effect for the various predictor variables included in the regression models.

In an ancillary analysis, we assessed the repeatability of RR measurements by examining the degree of agreement between repeat measurements that were taken within 5 min of each other using the Bland–Altman method of estimating the degree of agreement27. These repeat measurements were taken only when the first RR measurement was deemed raised for the child’s age20. We also examined the variations in the absolute difference between paired repeat measurements using the one-way ANOVA. Intra-observer variability was defined as the variation in repeated measurements by the same observer, and the intra-class correlation coefficient (ICC) was used to assess the reliability of these repeat measurements. The GAMLSS package23 in R statistical environment was used to estimate the centile curves whilst the remaining analyses were conducted using Stata 17 (Stata Corp, College Station, TX, USA).

Role of the funding source

The funders had no role in study design, data collection, data analysis, interpretation of results, manuscript writing, or the decision to submit the manuscript for publication. The corresponding author had full access to the data and made the final decision to submit the manuscript for publication following written approval from the co-authors.

Ethical approval

Ethical approval for the study was granted by The Gambia Government/Medical Research Council Joint Ethics Committee (reference: SCC 1390v3) and the School of Biological and Biomedical Sciences Ethics Committee, Durham University, Durham, UK (reference: SBBS/EC/1401/RooPfs). Written informed consent was obtained from the head of the household and the parent or guardian of all participating children. The study was conducted following the International Conference on Harmonization Tripartite Guideline for Good Clinical Practice and the Declaration of Helsinki.

Results

Reference curves

A total of 830 children recruited to the study contributed 67,512 observations in the analyses. The median age was 6.07 years (interquartile range (IQR), 4.21–8.55) and 400 (48.2%) were female. Details of study participants and the number of observations contributed at each stage of the analysis are presented in Fig. 1. The RR-for-age reference chart was constructed using 60,127 observations from 823 children (Fig. 2). There was a marked nonlinear decline in RR with increasing age between 1 and 6 years of age, after which the decline was minimal (the predicted median RR was 31 breaths/min [IQR 29–34] among 1-year-olds, 22 breaths/min [IQR 21–23] among 6-year-olds and 21 breaths/min [IQR 21–22] among 13-year-olds). The pre-specified cutoff for raised RR (> + 2SD) declined from 46 breaths/min among 1-year-olds to 27 breaths/min among 6-year-olds and to 25 breaths/min among 13-year-olds (Table 1, Fig. 2).Fig. 2 Respiratory rate (breaths/min) for age reference curves. The black dashed line represents the median RR measurements, the blue dashed lines represent the 1SD from the median, the red dashed lines represent 2SD from the median and the green continuous line represents the WHO cutoff for tachypnoea (40 breaths/min) among children aged > 12 months and < 5 years. Note that the WHO does not provide guidelines on tachypnoea among children aged above 5 years.

Table 1 Summary of children without signs of illness and their corresponding crude and model estimate of 2SD cutoff for raised respiratory rate (RR).

Age in years	Crude estimate of respiratory rate cutoff	GAMLSS estimate of RR cutoff	
Observations	Number of children	Median (IQR)	2SD (cutoff for raised RR)	Age in Years	Median (IQR)	2SD (model cutoff for raised RR)	
1 to < 1.5	410	20	32 (28–35)	41	1	31 (29–34)	46	
1.5 to < 2.5	2609	89	32 (28–35)	40	2	29 (27–31)	42	
2.5 to < 3.5	4308	142	30 (26–34)	38	3	26 (25–28)	37	
3.5 to < 4.5	5273	177	28 (24–33)	38	4	24 (23–25)	33	
4.5 to < 5.5	7180	218	25 (21–30)	38	5	22 (22–23)	29	
5.5 to < 6.5	7875	237	21 (19–25)	30	6	22 (21–23)	27	
6.5 to < 7.5	6753	209	20 (19–24)	29	7	22 (21–22)	26	
7.5 to < 8.5	6593	193	20 (18–23)	28	8	22 (21–22)	26	
8.5 to < 9.5	5821	167	20 (18–24)	28	9	22 (21–22)	26	
9.5 to < 10.5	4600	142	19 (18–22)	29	10	22 (21–22)	26	
10.5 to < 11.5	3881	118	20 (18–22)	28	11	22 (21–22)	26	
11.5 to < 12.5	3092	92	19 (18–22)	28	12	21 (21–22)	25	
12.5 to < 13.5	1478	49	19 (17–22)	26	13	21 (21–22)	25	

Respiratory rate among sick children

A total of 4119 (6.4%) observations from sick children were excluded from the construction of the reference chart. The median age of sick children was 5.5 years (IQR 3.6–8.0) and 2016 (48.9%) observations were from female subjects. The distribution of observations from children with and without signs of illness by age and RR are presented in Supplementary Table 1 and Supplementary Fig. 1.

Predictors of respiratory rate among community children without signs of illness

The prevalence of stunting (1–14 years of age) was 21.1% (95% CI 19.4–22.9%), the prevalence of underweight (1–10 years of age) was 22.6% (95% CI 20.6–24.6%) and prevalence of wasting (1 to < 5 years of age) was 8.0% (95% CI 5.9–10.6%). Females contributed 1037 (47.7%) observations in the 3 cross-sectional surveys. The distribution of observations over age and RR measurements are presented in Supplementary Fig. 1A and C, and Table 1.

In a pooled adjusted mixed-effects linear regression (n = 800 children, contributing 12,215 observations); age (non-linear effect), stunting, minimum daily temperature (a proxy for ambient temperature), and time of day of data collection were associated with respiratory rate (Table 2). On average, stunting was associated with a 0.84 breaths/min (95% CI 0.40–1.28, p < 0.001) increase in RR. A 1 °C increase in minimum daily temperature was associated with a 0.38 breaths/min (95% CI 0.33–0.42, p < 0.001) increase in RR. In addition, there was a nonlinear relationship between RR rate measurements and the time of the day when the measurements were taken (Table 2). On average RR measurements increased from 0700 to ~ 1500 h after which the measurements declined towards 2000 h (p < 0.001). While rainfall was significantly associated with RR in the univariable model, it was not significant in multivariable models that included minimum temperature. No evidence of being female or enrolled in the intervention group was associated with RR in the adjusted regression models (Table 2). As observed in the GAMLSS model and consistent with previous studies18,19, older children had lower RR measurements in both the unadjusted and the adjusted regression models (Table 2). That is, compared to children aged 1-to-2-year-old; children aged 2-to-3-years, 3-to-4 years, 4-to-6 years, 6-to-8 years, and 8-to-13 years old had RRs that were 0.38, 2.15, 5.14, 9.56, and 10.64 breaths/min lower respectively (Table 2). Table 2 Predictors of respiratory rate.

Covariate	Model 1: Univariable Analysis (n = 12,215)	Model 2: Multivariable Analysis (n = 12,215)	
Change in RR	95% CI	p value	Change in RR	95% CI	p value	
Age (ref, 1 to < 2 years)	ref			ref			
 2 to < 3 years	− 0.18	− 1.87 to 1.50	0.831	− 0.38	− 2.22 to 1.48	0.698	
 3 to < 4 years	− 1.77	− 3.45 to − 0.08	0.041	− 2.15	− 4.04 to − 0.27	0.025	
 4 to < 6 years	− 4.93	− 6.49 to − 3.37	 < 0.001	− 5.14	− 6.93 to − 3.35	 < 0.001	
 6 to < 8 years	− 9.15	− 10.65 to − 7.65	 < 0.001	− 9.56	− 11.32 to − 7.81	 < 0.001	
 8 to < 13 years	− 10.12	− 11.61 to − 8.62	 < 0.001	− 10.64	− 12.39 to − 8.90	 < 0.001	
Sex (ref, Male)	ref			ref			
 Female	0.11	− 0.50 to 0.72	0.720	0.21	− 0.12 to 0.54	0.213	
Randomization (ref, Control)	ref			ref			
 Treatment	0.45	− 0.15 to 1.06	0.142	0.19	− 0.14 to 0.52	0.258	
Time when RR measurements were taken (ref, 07:00–09:00 h)i	ref			ref			
 09:00–11:00 h	0.36	− 0.07 to 0.78	0.099	0.36	0.06– 0.65	0.017	
 11:00–13:00 h	0.68	0.15–1.20	0.011	0.78	0.46–1.10	 < 0.001	
 13:00–15:00 h	1.09	0.49–1.69	 < 0.001	1.02	0.61–1.42	 < 0.001	
 15:00–17:00 h	0.34	− 0.19 to 0.87	0.208	0.26	− 0.12 to 0.64	0.183	
 17:00–20:00 h	− 0.15	− 0.76 to 0.46	0.629	0.51	0.07 to 0.94	0.023	
Rainfall (mm) on the day of data collectionii	0.03	0.02–0.04	 < 0.001	0.001	− 0.01 to 0.01	0.935	
Minimum daily temp (oC)	0.31	0.26–0.36	 < 0.001	0.38	0.33 to 0.42	 < 0.001	
Stunting (ref, normal)iii	Ref						
 Stunted	2.46	1.66–3.27	 < 0.001	0.84	0.40 to 1.28	 < 0.001	
Wasting (ref, normal)iii	ref						
 Wasted	− 0.74	− 2.27 to 0.80	0.347				
Underweight (ref, normal)iii	ref						
 Underweight	0.47	− 0.24 to 1.18	0.191				
Model 1 presents results of the univariable mixed-effect linear regression model and Model 2 presents results from the final multivariable mixed-effect linear regression model. Sex and the randomization group have been included as potential confounders in the multivariable model.

iTime is defined as the time when the respiratory rate measurement was taken and was restricted to between 07:00 and 20:00 h respectively.

iiRainfall (mm) and ambient temperature on the day when RR measurements were taken were retrospectively extracted from the meteorological department in Basse Santa Su town and linked to the RR dataset.

iiiStunting was defined as the height/length for age z-score < − 2 (HAZ/LAZ < − 2), wasting was defined as the weight for height z-score < − 2 (WHZ < − 2), and underweight was defined as the weight for age < − 2 (WAZ < − 2). The 2006 and 2007 WHO algorithm estimates HAZ/LAZ for children 6 months to 14 years, WAZ for children under 10 years, and WHZ for children under 5 years.

Intra-observer variability in respiratory rate measurement

A total of 2311 paired measurements were used in the reliability assessment. 1,159 (50.2%) of the paired measurements were from female subjects and 460 (19.9%) were from children with signs of illness. Overall, the ICC for RR among all children with paired measurements (age 1–13 years) was 0.93 (95% CI 0.93–0.94) whilst the ICC for children without signs of illness was 0.84 (95% CI 0.83–0.86). Among children aged less than 5 years, the ICC was 0.94 (95% CI 0.92–0.96) while the estimate for children above 5 years and < 14 years was 0.83 (95% CI 0.82–0.84). The Bland–Altman plots of intra-observer reliability in RR measurement are presented in Fig. 3A and B. The mean difference between paired RR measurements was small in both children with (1.49 breaths/min [limits of agreement: − 0.74 to 3.71]) and without (1.04 breaths/min [limits of agreement: − 3.94 to 6.02]) signs of illness (Fig. 3A and B). However, the variability between the differences in paired RR measurements among sick children (SD = 2.54) was higher than among children without (SD = 1.13) signs of illness (one-way ANOVA, p < 0.001).Fig. 3 Bland–Altman plots for RR measurements. Panel A displays data from children without signs of illness whilst panel B displays data from children with signs of illness. The black dashed line represents the level of bias while the red dashed lines display the limits of agreement. RR1 is the first RRs measurement while RR2 is the second RRs measurement taken in that order. The mean difference in RR measurements was 1.49 breaths/min (limits of agreement: -0.74 to 3.71) for children without signs of illness and 1.04 breaths/min (limits of agreement: -3.94 to 6.02) for children with signs of illness. Data points are jittered for clarity.

Discussion

The WHO definition of tachypnoea does not address in any detail children above 5 years of age28. To construct reference charts for children aged 1 to 13 years we used a longitudinal dataset of RR measurements derived from a large cohort of rural Gambian children enrolled in a trial investigating the impact of improved housing on the incidence of malaria20,21. These data were linked to potential environmental exposures, demographic factors, and three cross-sectional surveys on participant’s nutritional status. RR measurements were done to monitor pneumonia, a potential serious adverse effect, since there was concern that the housing modifications in the intervention arm would increase respiratory disease—however, this was shown not to be the case20,21. The constructed RR-for-age reference chart among children without signs of illness decreased with increasing age and is consistent with previous findings17–19.

On average, stunting was associated with a small increase in RR consistent with a recent multicounty pneumonia trial29 in which stunting was associated with a RR of 1 breath/min higher at hospital admission. Whereas the multicounty trial assessed children presenting with signs of pneumonia29, our study examined RR measurements among mainly healthy children in the community. Stunted children have been reported to have lower ventilatory reserve than children with normal height for age30. Stunting has been associated with reduced lung function resulting from impaired lung development31–33 while chronic undernutrition predisposes children to severe illnesses such as pneumonia that can potentially lead to long-term effects on lung development32–35. Therefore, the relatively higher RRs observed among stunted children may be compensatory for inadequate ventilation. Our results, however, contrast with a previous study conducted in The Gambia that documented relatively lower RR measurements among stunted children with signs of pneumonia36 thus highlighting the need for further research to confirm these findings.

Higher ambient temperature was associated with increased RR measurements consistent with previous studies in both human and animal models37. The higher RR observed at higher ambient temperatures may be a homeostatic mechanism by which the body maintains its normal temperature37. On average, RR measurements taken in the afternoon were higher compared to those taken in the morning or in the evening. Higher RR in the afternoon may be a result of the cumulative increase in the child’s activity during the day. Increased activity (exercise) has been associated with an increase in oxygen requirements leading to an increase in RR38. Another potential explanation is that ambient temperature increases towards the afternoon and may thus confound the time when measurements were taken. However, the precise ambient temperature when RR measurements were taken was not available for these analyses.

Our study has important strengths. First, few studies provide a comprehensive assessment of the factors associated with RR using a large cohort of children who were actively followed up to document their clinical disease status, nutritional status, and other potential environmental risk factors. Second, the RR measurements were taken from children drawn from randomly selected households in the community and therefore the findings are generalizable to the community. In contrast, most RR studies have relied on cross-sectional studies where data are collected in hospital or daycare centers and therefore difficult to generalize to the community. Third, although the repeat RR measurements were only taken among children whose initial measurements were considered raised for age, the intra-observer reliability index (ICC) in both the overall and subgroup analyses was > 0.8 suggesting good to excellent reliability39. In addition, we did not find any systematic differences between the repeat measurements as shown in the Bland Altman plot (Fig. 3A and B). The variability between the RR repeat measurements was higher among children with signs of illness than among children without signs of illness (p < 0.001), highlighting a potential limitation of utilizing a single RR measurement for the classification of pneumonia and other respiratory illnesses.

Our study has some limitations. We did not record whether the child was asleep or awake when the RR measurements were taken. Previous studies have shown that RR measurements taken when children are asleep are lower than those taken when children are awake18. However, measurements were only taken when the child was calm. In addition, our data were collected from rural Gambia, a country situated less than 60m above sea level with high seasonal temperatures. Therefore, the interpretation of our results may be limited to areas with similar geographical and climatic characteristics17.

In conclusion, the current WHO tachypnea guidelines do not address in any detail children above 5 years of age28. We constructed centile charts using data collected from a low altitude area in The Gambia and propose a cutoff of > 26 breaths/min for raised RR among children > 5 years and below 13 years of age. We demonstrate that age, time of data collection, nutritional status, and ambient temperature were associated with RRs. However, the small change in RR associated with nutritional status, time of day when RR measurements were taken and ambient temperature does not warrant a change in the current WHO tachypnea guidelines for children > 1 year and < 5 years of age. Consistent with a previous study17, we note that the WHO definitions for fast breathing in children under 5 years may misclassify vulnerable children with acute respiratory signs, leading to inappropriate management and misallocation of scarce resources like antibiotics and oxygen. The finding that RR differences between paired repeat measurements were more variable among children with signs of illness raises concerns over the reliability of a single RR measurement among children with signs of illness—a population needing accurate diagnosis to enable targeted interventions with lifesaving treatment.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-70796-7.

Acknowledgements

The study was funded by the MRC-DfID-Wellcome Trust (Global Health Trials; MR/M007383/1). SWL is supported by the Global Challenges Research Fund’s BOVA Network (Global Challenges Research Fund (BB/R00532X/1) supported by the Biotechnology and Biological Sciences Research Council and Medical Research Council). JB received support from the UK Medical Research Council (MRC) and the UK Foreign, Commonwealth, and Development Office (FCDO) under the MRC/FCDO Concordat agreement and is also part of the EDCTP2 programme supported by the European Union (Grant Ref: MR/R010161/1). We are grateful for the support of the villagers, the village health workers, the regional health team, and the staff in the health clinics. We also thank the MRCG at the London School of Hygiene and Tropical Medicine for hosting this study in Basse Field Station and their technical and logistical support, especially the nurse field assistants and their supervisors; and the data supervisor, manager, and data entry staff.

Author contributions

Conceptualization and design: PM, JG, MP, UD’, SL, JB Performed the experiments: MP, UD’, SL, JB Data curation: PM, SA, JG, JB Funding acquisition: MP, SL Project administration and investigation: MP, SL Formal data analysis: PM, SA Writing- original draft: PM, SA All authors reviewed the manuscript.

Data availability

Access to clinical data requires a formal application to the Scientific Coordinating Committee of the Medical Research Council Unit The Gambia (MRCG) and the Joint Gambian Government’s and MRCG Ethics Committee in The Gambia at https://www.mrc.gm/scientific-coordinating-committee/

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Polycarp Mogeni and Sharon Amima.
==== Refs
References

1. Liu L Oza S Hogan D Chu Y Perin J Zhu J Global, regional, and national causes of under-5 mortality in 2000–15: An updated systematic analysis with implications for the Sustainable Development Goals Lancet Lond. Engl. 2016 388 10063 3027 3035 10.1016/S0140-6736(16)31593-8
Liu, L. et al. Global, regional, and national causes of under-5 mortality in 2000–15: An updated systematic analysis with implications for the Sustainable Development Goals. Lancet Lond. Engl. 388(10063), 3027–3035 (2016).10.1016/S0140-6736(16)31593-8
2. Bryce J Boschi-Pinto C Shibuya K Black RE WHO Child Health Epidemiology Reference Group WHO estimates of the causes of death in children Lancet Lond. Engl. 2005 365 9465 1147 1152 10.1016/S0140-6736(05)71877-8
Bryce, J., Boschi-Pinto, C., Shibuya, K., Black, R. E., WHO Child Health Epidemiology Reference Group. WHO estimates of the causes of death in children. Lancet Lond. Engl. 365(9465), 1147–1152 (2005).10.1016/S0140-6736(05)71877-8
3. GBD 2015 LRI Collaborators Estimates of the global, regional, and national morbidity, mortality, and aetiologies of lower respiratory tract infections in 195 countries: A systematic analysis for the Global Burden of Disease Study 2015 Lancet Infect. Dis. 2017 17 11 1133 1161 10.1016/S1473-3099(17)30396-1 28843578
GBD 2015 LRI Collaborators. Estimates of the global, regional, and national morbidity, mortality, and aetiologies of lower respiratory tract infections in 195 countries: A systematic analysis for the Global Burden of Disease Study 2015. Lancet Infect. Dis. 17(11), 1133–1161 (2017).28843578 10.1016/S1473-3099(17)30396-1
4. WHO. Revised WHO Classification and Treatment of Pneumonia in Children at Health Facilities: Evidence Summaries [Internet]. Geneva: World Health Organization; 2014 [cited 2023 Feb 24]. (WHO Guidelines Approved by the Guidelines Review Committee). Available from: http://www.ncbi.nlm.nih.gov/books/NBK264162/
5. Lanata CF Rudan I Boschi-Pinto C Tomaskovic L Cherian T Weber M Methodological and quality issues in epidemiological studies of acute lower respiratory infections in children in developing countries Int. J. Epidemiol. 2004 33 6 1362 1372 10.1093/ije/dyh229 15166188
Lanata, C. F. et al. Methodological and quality issues in epidemiological studies of acute lower respiratory infections in children in developing countries. Int. J. Epidemiol. 33(6), 1362–1372 (2004).15166188 10.1093/ije/dyh229
6. Nicolò A Massaroni C Schena E Sacchetti M The importance of respiratory rate monitoring: From healthcare to sport and exercise Sensors 2020 20 21 6396 10.3390/s20216396 33182463
Nicolò, A., Massaroni, C., Schena, E. & Sacchetti, M. The importance of respiratory rate monitoring: From healthcare to sport and exercise. Sensors 20(21), 6396 (2020).33182463 10.3390/s20216396
7. Goodman D Crocker ME Pervaiz F McCollum ED Steenland K Simkovich SM Challenges in the diagnosis of paediatric pneumonia in intervention field trials: Recommendations from a pneumonia field trial working group Lancet Respir. Med. 2019 7 12 1068 1083 10.1016/S2213-2600(19)30249-8 31591066
Goodman, D. et al. Challenges in the diagnosis of paediatric pneumonia in intervention field trials: Recommendations from a pneumonia field trial working group. Lancet Respir. Med. 7(12), 1068–1083 (2019).31591066 10.1016/S2213-2600(19)30249-8
8. Ginsburg AS Lenahan JL Izadnegahdar R Ansermino JM A Systematic review of tools to measure respiratory rate in order to identify childhood pneumonia Am. J. Respir. Crit. Care Med. 2018 197 9 1116 1127 10.1164/rccm.201711-2233CI 29474107
Ginsburg, A. S., Lenahan, J. L., Izadnegahdar, R. & Ansermino, J. M. A Systematic review of tools to measure respiratory rate in order to identify childhood pneumonia. Am. J. Respir. Crit. Care Med. 197(9), 1116–1127 (2018).29474107 10.1164/rccm.201711-2233CI
9. Noordam AC Barberá Laínez Y Sadruddin S van Heck PM Chono AO Acaye GL The use of counting beads to improve the classification of fast breathing in low-resource settings: A multi-country review Health Policy Plan. 2015 30 6 696 704 10.1093/heapol/czu047 24974104
Noordam, A. C. et al. The use of counting beads to improve the classification of fast breathing in low-resource settings: A multi-country review. Health Policy Plan. 30(6), 696–704 (2015).24974104 10.1093/heapol/czu047
10. Spence H Baker K Wharton-Smith A Mucunguzi A Matata L Habte T Childhood pneumonia diagnostics: Community health workers’ and national stakeholders’ differing perspectives of new and existing aids Glob. Health Action 2017 10 1 1290340 10.1080/16549716.2017.1290340 28485694
Spence, H. et al. Childhood pneumonia diagnostics: Community health workers’ and national stakeholders’ differing perspectives of new and existing aids. Glob. Health Action 10(1), 1290340 (2017).28485694 10.1080/16549716.2017.1290340
11. Ginsburg AS Sadruddin S Klugman KP Innovations in pneumonia diagnosis and treatment: A call to action on World Pneumonia Day, 2013 Lancet Glob Health 2013 1 6 e326 e327 10.1016/S2214-109X(13)70117-7 25104591
Ginsburg, A. S., Sadruddin, S. & Klugman, K. P. Innovations in pneumonia diagnosis and treatment: A call to action on World Pneumonia Day, 2013. Lancet Glob Health 1(6), e326–e327 (2013).25104591 10.1016/S2214-109X(13)70117-7
12. Källander K Young M Qazi S Universal access to pneumonia prevention and care: A call for action Lancet Respir. Med. 2014 2 12 950 952 10.1016/S2213-2600(14)70248-6 25466341
Källander, K., Young, M. & Qazi, S. Universal access to pneumonia prevention and care: A call for action. Lancet Respir. Med. 2(12), 950–952 (2014).25466341 10.1016/S2213-2600(14)70248-6
13. Baker K Maurel A Ward C Getachew D Habte T McWhorter C Automated respiratory rate counter to assess children for symptoms of pneumonia: Protocol for cross-sectional usability and acceptability studies in Ethiopia and Nepal JMIR Res. Protoc. 2020 9 3 e14405 10.2196/14405 32224491
Baker, K. et al. Automated respiratory rate counter to assess children for symptoms of pneumonia: Protocol for cross-sectional usability and acceptability studies in Ethiopia and Nepal. JMIR Res. Protoc. 9(3), e14405 (2020).32224491 10.2196/14405
14. Rawat, A. et al. Health System Considerations for Community-Based Implementation of Automated Respiratory Counters to Identify Childhood Pneumonia in 5 Regions of Ethiopia: A Qualitative Study. Int J Health Policy Manag [Internet]. 2023 Oct 2 [cited 2023 Oct 21]; Available from: https://www.ijhpm.com/article_4512.html
15. Tack B Vita D Mbaki TN Lunguya O Toelen J Jacobs J Performance of automated point-of-care respiratory rate counting versus manual counting in children under five admitted with severe febrile illness to Kisantu hospital, DR Congo Diagnostics 2021 11 11 2078 10.3390/diagnostics11112078 34829427
Tack, B. et al. Performance of automated point-of-care respiratory rate counting versus manual counting in children under five admitted with severe febrile illness to Kisantu hospital, DR Congo. Diagnostics 11(11), 2078 (2021).34829427 10.3390/diagnostics11112078
16. Baker K Ward C Maurel A de Cola MA Smith H Getachew D Usability and acceptability of a multimodal respiratory rate and pulse oximeter device in case management of children with symptoms of pneumonia: A cross-sectional study in Ethiopia Acta Paediatr. 2021 110 5 1620 1632 10.1111/apa.15682 33220086
Baker, K. et al. Usability and acceptability of a multimodal respiratory rate and pulse oximeter device in case management of children with symptoms of pneumonia: A cross-sectional study in Ethiopia. Acta Paediatr. 110(5), 1620–1632 (2021).33220086 10.1111/apa.15682
17. Crocker ME Hossen S Goodman D Simkovich SM Kirby M Thompson LM Effects of high altitude on respiratory rate and oxygen saturation reference values in healthy infants and children younger than 2 years in four countries: A cross-sectional study Lancet Glob. Health 2020 8 3 e362 e373 10.1016/S2214-109X(19)30543-1 32087173
Crocker, M. E. et al. Effects of high altitude on respiratory rate and oxygen saturation reference values in healthy infants and children younger than 2 years in four countries: A cross-sectional study. Lancet Glob. Health 8(3), e362–e373 (2020).32087173 10.1016/S2214-109X(19)30543-1
18. Herbert A Pearn J Wilson S Normal percentiles for respiratory rate in children—Reference ranges determined from an optical sensor Children 2020 7 10 160 10.3390/children7100160 33023258
Herbert, A., Pearn, J. & Wilson, S. Normal percentiles for respiratory rate in children—Reference ranges determined from an optical sensor. Children 7(10), 160 (2020).33023258 10.3390/children7100160
19. Fleming S Thompson M Stevens R Heneghan C Plüddemann A Maconochie I Normal ranges of heart rate and respiratory rate in children from birth to 18 years of age: A systematic review of observational studies The Lancet 2011 377 9770 1011 1018 10.1016/S0140-6736(10)62226-X
Fleming, S. et al. Normal ranges of heart rate and respiratory rate in children from birth to 18 years of age: A systematic review of observational studies. The Lancet 377(9770), 1011–1018 (2011).10.1016/S0140-6736(10)62226-X
20. Pinder M Bradley J Jawara M Affara M Conteh L Correa S Improved housing versus usual practice for additional protection against clinical malaria in The Gambia (RooPfs): A household-randomised controlled trial Lancet Planet Health 2021 5 4 e220 e229 10.1016/S2542-5196(21)00002-4 33838737
Pinder, M. et al. Improved housing versus usual practice for additional protection against clinical malaria in The Gambia (RooPfs): A household-randomised controlled trial. Lancet Planet Health 5(4), e220–e229 (2021).33838737 10.1016/S2542-5196(21)00002-4
21. Pinder M Conteh L Jeffries D Jones C Knudsen J Kandeh B The RooPfs study to assess whether improved housing provides additional protection against clinical malaria over current best practice in The Gambia: Study protocol for a randomized controlled study and ancillary studies Trials 2016 17 1 275 10.1186/s13063-016-1400-7 27255167
Pinder, M. et al. The RooPfs study to assess whether improved housing provides additional protection against clinical malaria over current best practice in The Gambia: Study protocol for a randomized controlled study and ancillary studies. Trials 17(1), 275 (2016).27255167 10.1186/s13063-016-1400-7
22. Mackenzie GA Hill PC Jeffries DJ Hossain I Uchendu U Ameh D Effect of the introduction of pneumococcal conjugate vaccination on invasive pneumococcal disease in The Gambia: A population-based surveillance study Lancet Infect. Dis. 2016 16 6 703 711 10.1016/S1473-3099(16)00054-2 26897105
Mackenzie, G. A. et al. Effect of the introduction of pneumococcal conjugate vaccination on invasive pneumococcal disease in The Gambia: A population-based surveillance study. Lancet Infect. Dis. 16(6), 703–711 (2016).26897105 10.1016/S1473-3099(16)00054-2
23. Stasinopoulos DM Rigby RA Generalized additive models for location scale and shape (GAMLSS) in R J. Stat. Softw. 2008 23 1 46
Stasinopoulos, D. M. & Rigby, R. A. Generalized additive models for location scale and shape (GAMLSS) in R. J. Stat. Softw. 23, 1–46 (2008).
24. Rigby RA Stasinopoulos DM Using the Box–Cox t distribution in GAMLSS to model skewness and kurtosis Stat. Model. 2006 6 3 209 229 10.1191/1471082X06st122oa
Rigby, R. A. & Stasinopoulos, D. M. Using the Box–Cox t distribution in GAMLSS to model skewness and kurtosis. Stat. Model. 6(3), 209–229 (2006).10.1191/1471082X06st122oa
25. Borghi E de Onis M Garza C Van den Broeck J Frongillo EA Grummer-Strawn L Construction of the World Health Organization child growth standards: Selection of methods for attained growth curves Stat. Med. 2006 25 2 247 265 10.1002/sim.2227 16143968
Borghi, E. et al. Construction of the World Health Organization child growth standards: Selection of methods for attained growth curves. Stat. Med. 25(2), 247–265 (2006).16143968 10.1002/sim.2227
26. de Onis M Onyango AW Borghi E Siyam A Nishida C Siekmann J Development of a WHO growth reference for school-aged children and adolescents Bull. World Health Organ. 2007 85 9 660 667 10.2471/BLT.07.043497 18026621
de Onis, M. et al. Development of a WHO growth reference for school-aged children and adolescents. Bull. World Health Organ. 85(9), 660–667 (2007Sep).18026621 10.2471/BLT.07.043497
27. Bland JM Altman DG Statistical methods for assessing agreement between two methods of clinical measurement Lancet Lond. Engl. 1986 1 8476 307 310 10.1016/S0140-6736(86)90837-8
Bland, J. M. & Altman, D. G. Statistical methods for assessing agreement between two methods of clinical measurement. Lancet Lond. Engl. 1(8476), 307–310 (1986).10.1016/S0140-6736(86)90837-8
28. WHO. Pocket book of primary health care for children and adolescents: guidelines for health promotion, disease prevention and management from the newborn period to adolescence [Internet] (World Health Organization. Regional Office for Europe, accessed 21 October 2023); https://apps.who.int/iris/handle/10665/352485 (2022)
29. Moschovis PP Addo-Yobo EOD Banajeh S Chisaka N Christiani DC Hayden D Stunting is associated with poor outcomes in childhood pneumonia Trop. Med. Int. Health TMIH 2015 20 10 1320 1328 10.1111/tmi.12557
Moschovis, P. P. et al. Stunting is associated with poor outcomes in childhood pneumonia. Trop. Med. Int. Health TMIH 20(10), 1320–1328 (2015).10.1111/tmi.12557
30. Nguyen TTD Hoo AF Lum S Wade A Thia LP Stocks J New reference equations to improve interpretation of infant lung function Pediatric Pulmonol. 2013 48 4 370 380 10.1002/ppul.22656
Nguyen, T. T. D. et al. New reference equations to improve interpretation of infant lung function. Pediatric Pulmonol. 48(4), 370–380 (2013).10.1002/ppul.22656
31. Narayanan M Owers-Bradley J Beardsmore CS Mada M Ball I Garipov R Alveolarization continues during childhood and adolescence Am. J. Respir. Crit. Care Med. 2012 185 2 186 191 10.1164/rccm.201107-1348OC 22071328
Narayanan, M. et al. Alveolarization continues during childhood and adolescence. Am. J. Respir. Crit. Care Med. 185(2), 186–191 (2012).22071328 10.1164/rccm.201107-1348OC
32. Tennant PWG Gibson GJ Pearce MS Lifecourse predictors of adult respiratory function: results from the Newcastle Thousand Families Study Thorax 2008 63 9 823 830 10.1136/thx.2008.096388 18408051
Tennant, P. W. G., Gibson, G. J. & Pearce, M. S. Lifecourse predictors of adult respiratory function: results from the Newcastle Thousand Families Study. Thorax 63(9), 823–830 (2008).18408051 10.1136/thx.2008.096388
33. Mukhopadhvay S Macleod KA Ong TJ Ogston SA “Ethnic” variation in childhood lung function may relate to preventable nutritional deficiency Acta Paediatr. 2001 90 11 1299 1303 10.1111/j.1651-2227.2001.tb01579.x 11808903
Mukhopadhvay, S., Macleod, K. A., Ong, T. J. & Ogston, S. A. “Ethnic” variation in childhood lung function may relate to preventable nutritional deficiency. Acta Paediatr. 90(11), 1299–1303 (2001).11808903 10.1111/j.1651-2227.2001.tb01579.x
34. Kaali S Jack DW Prah RKD Chillrud SN Mujtaba MN Kinney PL Poor early childhood growth is associated with impaired lung function: Evidence from a Ghanaian pregnancy cohort Pediatric Pulmonol. 2022 57 9 2136 2146 10.1002/ppul.26015
Kaali, S. et al. Poor early childhood growth is associated with impaired lung function: Evidence from a Ghanaian pregnancy cohort. Pediatric Pulmonol. 57(9), 2136–2146 (2022).10.1002/ppul.26015
35. Barker DJ Godfrey KM Fall C Osmond C Winter PD Shaheen SO Relation of birth weight and childhood respiratory infection to adult lung function and death from chronic obstructive airways disease BMJ 1991 303 6804 671 675 10.1136/bmj.303.6804.671 1912913
Barker, D. J. et al. Relation of birth weight and childhood respiratory infection to adult lung function and death from chronic obstructive airways disease. BMJ 303(6804), 671–675 (1991).1912913 10.1136/bmj.303.6804.671
36. Falade AG Tschappeler H Greenwood BM Mulholland EK Use of simple clinical signs to predict pneumonia in young Gambian children: The influence of malnutrition Bull. World Health Organ. 1995 73 6
Falade, A. G., Tschappeler, H., Greenwood, B. M. & Mulholland, E. K. Use of simple clinical signs to predict pneumonia in young Gambian children: The influence of malnutrition. Bull. World Health Organ. 73, 6 (1995).
37. Zila I Calkovska A Effects of elevated body temperature on control of breathing Acta Med. Martiniana 2011 11 Supplement-1 24 30 10.2478/v10201-011-0008-6
Zila, I. & Calkovska, A. Effects of elevated body temperature on control of breathing. Acta Med. Martiniana 11(Supplement-1), 24–30 (2011).10.2478/v10201-011-0008-6
38. Nicolò A Massaroni C Passfield L Respiratory frequency during exercise: The neglected physiological measure Front. Physiol. 2017 10.3389/fphys.2017.00922 29321742
Nicolò, A., Massaroni, C. & Passfield, L. Respiratory frequency during exercise: The neglected physiological measure. Front. Physiol.10.3389/fphys.2017.00922 (2017).29321742 10.3389/fphys.2017.00922
39. Koo TK Li MY A guideline of selecting and reporting intraclass correlation coefficients for reliability research J. Chiropr. Med. 2016 15 2 155 163 10.1016/j.jcm.2016.02.012 27330520
Koo, T. K. & Li, M. Y. A guideline of selecting and reporting intraclass correlation coefficients for reliability research. J. Chiropr. Med. 15(2), 155–163 (2016).27330520 10.1016/j.jcm.2016.02.012
