Article Text

Original research
Characterising airway obstructive, dysanaptic and PRISm phenotypes of prematurity-associated lung disease
  1. Michael Cousins1,
  2. Kylie Hart1,
  3. Sarah J Kotecha1,
  4. A John Henderson2,
  5. W John Watkins1,
  6. Andrew Bush3,
  7. Sailesh Kotecha1
  1. 1 Department of Child Health, Cardiff University School of Medicine, Cardiff, UK
  2. 2 MRC Integrative Epidemiology Unit, Population Health Sciences, Bristol Medical School, University of Bristol, Bristol, UK
  3. 3 Centre for Paediatrics and Child Health, Imperial College of Medicine, London, UK
  1. Correspondence to Professor Sailesh Kotecha, Department of Child Health, Cardiff University School of Medicine, Cardiff CF14 4XN, UK; kotechas{at}cardiff.ac.uk

Abstract

Introduction Although obstructive airway disease has been shown to be associated with prematurity, other spirometry phenotypes are less well described.

Objectives We characterised abnormal spirometry phenotypes in preterm-born children, including prematurity-associated obstructive lung disease (POLD, forced expiratory volume in 1 s (FEV1)<lower limit of normal (LLN), FEV1/forced vital capacity (FVC)<LLN), prematurity-associated preserved ratio of impaired spirometry (pPRISm, FEV1<LLN, FEV1/FVC≥LLN) and prematurity-associated dysanapsis (pDysanapsis, FEV1≥LLN, FEV1/FVC<LLN), and associated them with early life factors, bronchodilator responsiveness and fractional exhaled nitric oxide (FENO).

Methods 768 children, aged 7–12 years, underwent FENO measurements and spirometry before and after salbutamol. Groups were compared using parametric tests; multinomial regression was used.

Results 22.6% of 544 preterm-born (mean gestation: 31 weeks) and 9.2% of 195 term-born children, with satisfactory data available, were classified into one of four abnormal spirometry groups. Each phenotype was generally more prevalent in preterm-born children than in the term-born children. For the preterm group, POLD-reversible (4.4%) was associated with increased FENO, bronchopulmonary dysplasia (BPD) and intrauterine growth restriction. POLD-fixed group (3.3%) did not have increased FENO but was associated with BPD. 41% of the pDysanapsis group (5.9%) had bronchodilator response, 31% had increased FENO and was associated with postnatal weight gain. In the pPRISm group (9%), 13% responded to bronchodilators, FENO was not increased and was non-significantly associated with body mass index (p=0.064).

Conclusions Further to airway obstruction, we describe airway dysanapsis and pPRISm spirometry phenotypes in survivors of prematurity, both of which have poor outlook in other disease groups. By identifying specific phenotypes, targeted therapy can be developed to improve long-term outcomes.

  • paediatric lung disaese

Data availability statement

Data from the RHiNO study are available to research collaborators subject to confidentiality and non-disclosure agreements. Contact Professor Sailesh Kotecha (kotechas{at}cardiff.ac.uk) for any data requests.

https://creativecommons.org/licenses/by/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution 4.0 Unported (CC BY 4.0) license, which permits others to copy, redistribute, remix, transform and build upon this work for any purpose, provided the original work is properly cited, a link to the licence is given, and indication of whether changes were made. See: https://creativecommons.org/licenses/by/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

WHAT IS ALREADY KNOWN ON THIS TOPIC

  • Obstructive lung disease following preterm birth has been previously described; however, little is known about other potential patterns or phenotypes of lung disease in such children.

WHAT THIS STUDY ADDS

  • This study characterises four groups of preterm-born children with evidence of lung dysfunction, each with differential associations with spirometry, exhaled nitric oxide and reversibility, as well as early and current life factors.

HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY

  • Using the clinical phenotypes we have identified, spirometry measurements can be used to identify children who can benefit from treatment to improve lung dysfunction, and also identify children who need long-term surveillance due to later comorbidities associated with their pattern of lung disease.

Introduction

Low lung function observed in prematurity-associated lung disease (PLD),1–3 including those who develop bronchopulmonary dysplasia (BPD, also known as chronic lung disease of prematurity), is now considered a precursor for early onset chronic obstructive pulmonary disease (COPD).4 Although historical focus has been on preterm-born survivors who had BPD in infancy,5 6 it is increasingly recognised that late-preterm-born children (especially those born at 33–36 weeks’ gestation) also have lung function decrements.1 7 We have recently reported that gestational age at birth and intrauterine growth restriction (IUGR) are better predictors of low lung function in childhood than BPD in multivariable models.8 Our systematic review showed improved spirometry after acute administration of bronchodilators, but the specific group of preterm-born children with PLD that would benefit the most is unclear.9 Furthermore, fractional exhaled nitric oxide (FENO) has not shown to be increased in children who had BPD in infancy.10 In contrast, our recent data showed that FENO was decreased by inhaled corticosteroids (ICS) either alone or in combination with long-acting beta2 receptor agonists (LABA) but not after placebo treatment in a randomised controlled trial suggesting that some children with PLD may have eosinophilic airway disease.11

However, while airway obstruction has been described after preterm birth, little is known about other patterns of abnormal spirometry of PLD, including preserved ratio of impaired spirometry (PRISm) and dysanaptic airway growth. PRISm is an important phenotype as it is associated with increased risk of developing COPD and cardiovascular disease and with all-cause mortality in middle-aged and old-aged adult populations,12 13 but has yet to be evaluated in children including in the preterm population; and dysanaptic airway growth is associated with worse asthma.14

We hypothesised that there are multiple spirometry phenotypes in PLD, including:

  • prematurity-associated obstructive lung disease that is bronchodilator responsive (POLD-reversible),

  • persistent airflow limitation of prematurity that does not respond to bronchodilators (POLD-fixed),

  • prematurity-associated PRISm (pPRISm),

  • dysanapsis of prematurity (pDysanapsis).

We also examined early and current life factors which were associated with each phenotype, in order to begin to understand how to personalise preventative strategies.

Methods

Study population

The Respiratory Health Outcomes in Neonates (RHiNO, EudraCT: 2015-003712-20) study has been described previously.8 11 Briefly, children from a previous questionnaire study1 15 were supplemented with additional preterm-born children, and were mailed a respiratory questionnaire for their parents to complete. Responders to the questionnaire were invited to join RHiNO for a home or hospital visit to obtain anthropometric details, perinatal and respiratory history (supplemented by examination of the child’s medical records), measurement of FENO (NIOX VERO, Circassia, Oxford, UK) and spirometry (Microloop, CareFusion, Wokingham, UK) before and after administration of salbutamol (400 μg administered via a spacer device) as described in the online supplemental file. Spirometry was performed and quality controlled according to the American Thoracic Society/European Respiratory Society guidelines,16 and per cent predicted and z-scores were calculated using Global Lung Initiative reference equations.17

Supplemental material

Children aged 7–12 years, born at ≤34 weeks’ gestation for the preterm group or at ≥37 weeks’ gestation for the term control group were invited. Children born at 35–36 weeks’ gestation were not included as they have previously been shown to have lesser degrees of lung function deficits.7 Recruitment took place prospectively between November 2016 and September 2019. A subgroup underwent extended assessments including skin prick testing (see online supplemental file). Children with significant congenital, cardiac or neurodevelopmental abnormalities were excluded. Z-scores for birth weight, current weight, height and BMI and IUGR (defined as <10th centile for birth weight adjusted for sex and gestation) were calculated using the LMS growth program.18 Deprivation scores and quintiles were estimated from the participants’ postcodes using the Welsh Index of Multiple Deprivation score, which is a validated measure of deprivation based on eight domains including wealth, schooling and home ownership.19 BPD was defined as supplemental oxygen-dependency of 28 days of age or greater for those born <32 weeks’ gestation and at 56 days of age for those ≥32 weeks’ gestation).20

Definitions of spirometry phenotypes

POLD in the preterm group was defined as FEV1 and FEV1/FVC ratio less than the lower limit of normal (LLN) and further subdivided into reversible (POLD-reversible) and fixed airway (POLD-fixed) disease if their FEV1 improved by 10% or not, respectively, after bronchodilator administration. Dysanapsis was defined according to the American Thoracic Society/European Respiratory Society (ATS/ERS) definition as FEV1≥LLN and FEV1/FVC<LLN21 and pPRISm as FEV1<LLN and FEV1/FVC ratio≥LLN.11 12 Equivalent phenotypes were defined similarly for the term group. Preterm-born and term-born children with FEV1≥LLN and FEV1/FVC ratio≥LLN were considered as control groups. These definitions and other abbreviations are described in table 1.

Table 1

Abbreviations and definitions for grouping participants based on lung function

Statistical analyses

Data are given as means and SD, or medians and ranges, as appropriate. One-way analysis of variance with Bonferroni correction were used to compare means as appropriate. The χ2 test was used to compare independence of categorical variables. ORs (and 95% CIs) were calculated to compare each phenotype in the preterm group against the phenotype in the term group. Multinomial regression was used to identify early and current life factors that were associated with spirometry phenotypes. Data were analysed using SPSS V.26 (IBM, New York, USA). P value <0.05 was considered statistically significant.

Results

From 1426 returned questionnaires, 768, including 565 preterm-born and 203 term-born, children participated in a home or hospital visit. As previously reported, lower gestational age and birth weight; and fewer of the most-deprived population were noted for responders (31.0 weeks, 1703 g) and non-responders (31.6 weeks, 1828 g) to the questionnaires.8 Twenty-one preterm-born and 8 term-born children’s spirometry did not pass quality control criteria thus spirometry data were available from 544 preterm-born and 195 term-born children (table 2). Two hundred seventy-six children were born at <32 weeks’ gestation including 99 born at 28 weeks or less. One hundred eight children had mild/moderate/severe BPD diagnosed in infancy.

Table 2

Participant demographics by phenotypes

Classification of preterm-born and term-born children into spirometry phenotypes

Overall, 123/544 (22.6%) preterm-born children compared with 18/195 (9.2%) of term-born children were classified into one of the four spirometry phenotypes giving OR and 95% CIs of 2.87 (1.70 to 4.86, p=0.0001) when the preterm and term groups were compared (figure 1). For the preterm group, 91 (16.7%) with FEV1<LLN were further classified into 42 (7.7%) with FEV1/FVC ratio<LLN into who had an obstructive phenotype (POLD) and into 49 (9.0%) with FEV1/FVC ratio≥LLN into pPRISm phenotype. The POLD group was further classified into those who had (POLD-reversible, n=24, 4.4%) or did not have (POLD-fixed, n=18, 3.3%) bronchodilator responsiveness. Preterm-born children with FEV1≥LLN and FEV1/FVC ratio<LLN were classified into the pDysanapsis group (n=32, 5.9%) and the remaining 421 (77.4%) with FEV1 and FEV1/FVC ratio≥LLN were classified as the preterm control group. For the term group, 2 (1%, OR; 95% CI 5.05; 1.18 to 21.58; p=0.029), 1 (0.5%, 7.57; 1.00 to 57.12, p=0.049), 9 (4.6%, 2.29; 1.09 to 4.76; p=0.026) and 6 (4.6%; 2.24; 0.92 to 5.56; p=0.075) were classified into reversible and fixed obstructive, pPRISm and dysanapsis groups, respectively. The OR and 95% CI compared each phenotype between the preterm and term groups.

Figure 1

Classification of preterm and term groups based on their spirometry measures and bronchodilators responses. FEV1, forced expiratory volume in 1 s; FVC, forced vital capacity; LLN, lower limit of normal.

Participant demographics

When the six groups were compared (table 2), the pPRISm group were marginally older, but weight, height and BMI were similar between the groups. The POLD-reversible group had more IUGR. BPD was more prevalent in the two POLD groups regardless of bronchodilator response. The two POLD groups and the pDysanapsis group were associated with more respiratory symptoms, doctor-diagnosed asthma and inhaler use. Spirometry was lower in all phenotypes identified when compared with the term control group with the lowest measures noted in the two POLD groups. Preterm controls had very similar spirometry to the term control group, hence were used as reference population in multinomial regression analyses.

Bronchodilator response

Figure 2 shows the bronchodilator responses for the six groups. The largest and least increases were observed in the POLD-reversible group and POLD-fixed group (expected as defined on basis of bronchodilator response). FEV1 increased to ≥LLN after bronchodilator administration in 16/24 (67%) in the POLD-reversible group suggesting potential to reach spirometry within normal limits. Fewer of pPRISm (13%); and preterm (6%) or term (5%) control groups had a response but 41% of the pDysanapsis group responded to bronchodilators.

Figure 2

(A) Mean change in bronchodilator plotted against the different groups and (B) table showing the proportion of subjects with a positive bronchodilator response defined as 10% improvement in postbronchodilator %FEV1. Significant versus †POLD-fixed, ‡pPRISm, ¥pDysanapsis, ¶PTc, ₴Tc (single symbol denotes significance level <0.05, double symbol <0.01, triple symbol <0.001). FEV1, forced expiratory volume in 1 s; FVC, forced vital capacity; pDysanapsis, dysanapsis of prematurity; POLD, prematurity-associated obstructive lung disease; pPRISm, prematurity-associated PRISm.

FENO levels

FENO levels were increased the most in the POLD-reversible group (mean: 36.1, 95% CI 23.7 to 48.6 ppb) and the pDysanapsis group (27.9; 16.1 to 39.7) but were similar in the pPRISm, POLD-fixed and control groups (figure 3); 45% and 31%, respectively of the POLD-reversible and pDysanapsis groups had FENO greater than the National Institute for Health and Care Excellence and ATS recommended 35 ppb, but this proportion was 13% or less in the remaining groups.22 23

Figure 3

(A) Fractional exhaled nitric oxide (FENO) levels plotted against the different groups and (B) table showing the proportion of subjects with FENO>35 ppb. Significant versus *POLD-reversible, †PLD-fixed, ‡pPRISm, ¥pDysanapsis, ¶PTc, ₴Tc (single symbol denotes significance level <0.05, double symbol <0.01, triple symbol <0.001). pDysanapsis, dysanapsis of prematurity; POLD, prematurity-associated obstructive lung disease; pPRISm, prematurity-associated PRISm.

Skin prick testing

A subgroup of 237/739 (32%) children underwent skin prick testing (table 3). Eleven per cent of the term group and 20%–23% of the POLD-fixed, pPRISm and preterm controls had positive tests. Although the POLD-reversible group had greater positivity than the term groups (33% vs 11%, p<0.05), this difference was not significantly different when compared with preterm-born children (33% vs 23%).

Table 3

Skin prick testing for the different groups shown

Association of early and current life factors with the four phenotypes

Finally, we used multinomial regression to assess which early and current life factors, including sex, BPD, IUGR, antenatal maternal smoking, deprivation, current BMI and postnatal weight change in z-scores from birth to current weight, were associated with each of abnormal spirometry group with preterm controls as the reference population (table 4). Univariable analyses showed that BPD (OR 3.35; 95% CI 1.43 to 7.83) and IUGR (4.56; 1.89 to 10.96) were associated with the subsequent development of POLD-reversible; and BPD (3.75; 1.43 to 9.83) with the development of POLD-fixed. Weight gain change of z-scores from birth to the current weight (1.37; 1.07 to 1.76) was associated with development of pDysanapsis. The association of BMI z-scores with pPRISm (0.79; 0.62 to 1.01, p=0.064) and IUGR with pDysanapsis (2.13; 0.87 to 5.17, p=0.097) had p values of <0.1. Thus, BPD was associated with the POLD-fixed group, and possibly BMI z-scores with pPRISm. Since the POLD-reversible group was associated with both BPD and IUGR; and the pDysanapsis group with both weight gain and marginally with IUGR (table 4), these factors were included in a multivariable multinomial regression model with outcomes of POLD-reversible and pDysanapsis groups assessed against the preterm controls as the reference population. The results confirmed that POLD-reversible was associated with BPD and IUGR and pDysanapsis group with postnatal weight gain (table 4).

Table 4

(A) Univariable and (B) multivariable multinomial regression analyses using early and current life factors to associate with different spirometry phenotypes in the preterm population only

Discussion

As expected, we identified groups of children with fixed and variable airflow obstruction including, and, contrary to previous literature,10 a subgroup of children with increased FENO. We have also demonstrated that dysanaptic airway growth and PRISm are part of the spectrum of lung abnormalities following preterm birth. Thus, spirometry identified four phenotypes of PLD, representing 22.6% of the preterm population compared with 9.2% in the term group. We have also reported associations of each phenotype with presence of early or current life factors with BPD, IUGR, BMI and weight gain between birth and current weight being most relevant.

The study of respiratory outcomes after preterm birth has largely focused on those who survived the diagnosis of BPD in infancy. Increasingly, it is recognised that many preterm-born survivors of BPD do not develop respiratory disease in later life.8 In contrast, many born at later gestations especially at 33–34 weeks’ gestation are also at risk of later lung function decrements.8 Furthermore, it is unclear which groups are associated with bronchodilator responses and with FENO, which has not shown to be increased in those who had BPD in infancy.10

Several studies show that survivors of prematurity especially those who had BPD in infancy respond to single doses of bronchodilators but long-term evaluation of bronchodilators is limited.9 We identified two obstructive airway groups, one with low FEV1 (POLD groups) and a second in which FEV1 was ≥LLN and FEV1/FVC ratio<LLN (pDysanapsis) as suggested by the ATS/ERS guidelines.15 We subclassified the POLD group into those who did and did not respond to bronchodilators. In addition, a subgroup of the PRISm patients (13%) was also bronchodilator responsive. Bronchodilator responsiveness is an important group as it is reasonable to speculate that these children, especially those whose FEV1 improved to greater than LLN, should benefit from inhaler treatment with combined ICS and LABA as we have recently reported.11 It should be noted that the evidence for corticosteroid responsive airway eosinophilia is not present in most POLD children, but there are insufficient safety data for using LABA alone in this group. In our previous meta-analyses, FENO was not different between term and preterm groups including those who had BPD in infancy.10 Only 45% of the POLD-reversible group had an increased FENO thus mechanisms other than atopy and airway eosinophilia are important in this group. In support of this, only 33% had a positive skin test. Inevitably, a number of children born preterm will have coincident atopic asthma.1 Preterm-born infants who die from BPD have been shown to have smooth muscle extending much further distally in the airways at autopsy than non-BPD and term-born controls.24 25 It is likely that the smooth muscle remains responsive to bronchodilators at least in a proportion of the population with POLD. On multivariable modelling, both BPD and IUGR were strongly associated with the development of POLD-reversible disease. Historically, focus for long-term respiratory outcomes has been on preterm-born children with BPD.2 26 The BPD group with low lung function were more likely to develop POLD rather than pPRISm, when compared with non-BPD preterm-born children with low lung function. This may be due to the greater number of interventions, including respiratory support and prolonged supplemental oxygen, that these children will have been exposed to in the neonatal period. While this is speculative, some support is provided by the above mentioned increased distal airway smooth muscle extension that has been observed in infants who die from BPD when compared with non-BPD controls.24 25 Additional factors such as maternal antenatal smoking and postnatal exposure to smoking and pollution as well as antenatal and postnatal use of corticosteroids in infancy and childhood could potentially lead to, and/or potentially modify, the eventual phenotype. However, prospective studies will be required, especially to evaluate if treatment can modulate the lung disease and potentially modulate the phenotype.

The POLD-fixed group may potentially have progressed from a reversible to a fixed phenotype, but this is speculative. This lack of reversibility may be explained by the distal extension of smooth muscle hypertrophy remodelling into myofibroblasts which is unlikely to respond to inhaled therapy. Alternatively, failure of normal airway development may be a factor. This group did not have increased FENO and were no more likely to be skin prick test positive than the preterm group (22% vs 23%). Interestingly, the POLD-fixed group was associated with BPD in multinomial regression modelling but not with IUGR thus again suggesting that early intervention with mechanical ventilation, supplemental oxygen therapy, etc in the neonatal period is likely to result in obstructive airway disease possible due to lung injury to the parenchyma and to the airways.

Only 41% of the pDysanapsis group responded to bronchodilators, and 31% had FENO over 35 ppb. Positive skin prick tests were noted in 30% against 11% and 23% in the term and preterm groups. Early weight gain in infancy in both term and preterm born children has been shown to result in dysanapsis.27 28 On univariable analyses, weight gain between birth and current weight was significantly (p=0.013) and IUGR was weakly (p=0.097) associated with the pDysanapsis group with only the former remaining significantly associated on multivariable regression analyses. This group was not associated with BPD. The pDysanapsis group had increased respiratory symptoms, doctor diagnosed asthma and inhaler use. Of concern is that this group of children is likely to have similar worse outcomes as observed for asthma in school age children and adults.14

Although it is tempting to dismiss the group with %FEV1<LLN and FEV1/FVC ratio≥LLN as children with small lungs, longitudinal studies into late adulthood are increasingly showing that low lung function trajectories lead to premature development of COPD even without accelerated decline.29–31 Furthermore, PRISm phenotype is increasingly being recognised and associated with increased cardiorespiratory morbidity and all-cause mortality.12 13 Most studies have been of adults at middle age or beyond with little described in children. Since prematurity is associated with early life decreases in spirometry, it is very likely that this group of subjects will progress to develop a prematurity-associated PRISm phenotype. As with the observations in adults, we noted no increase in FENO and a subgroup of 13% responded to bronchodilators. On univariable multinomial regression, there was no association between weight gain between birth and current age but there was a near statistically significant negative association with BMI (beta=−0.23, p=0.064). There was no effect of sex or BPD.

We have shown that there are three different phenotypes that may respond to bronchodilators, some of which are associated with increased FENO and atopy. We have also recently showed that preterm-born children with POLD have decreased exercise capacity as well as increased functional residual capacity (FRC) and residual volume (RV) when compared with term-born controls.32 In contrast, those with pPRISm had intermediate decreases in exercise capacity measures but had decreased FRC and RV. Some but not all small studies thus far have shown an increase in neutrophils33 or oxidant injury34 but evidence for airway eosinophilic inflammation in this group is lacking. Credence is given to an active role for the increased FENO in the obstructive airway phenotypes as we noted decreased FENO in a recent randomised controlled trial in both arms of ICS alone and when used in combination with LABA. Clinically inhaled corticosteroids non-significantly improved %FEV1 by ~7% but the combination of ICS and LABA resulted in >14% significant improvement. The current data suggest that such a trial of therapy is not unreasonable for the POLD-obstructive and pDysanapsis groups and a supervised prolonged trial may be of benefit to those who do not respond to single doses of bronchodilators, especially if any inflammatory component contributes to the underlying airway disease. Whether the pPRISm group and the POLD-fixed group respond to prolonged therapy is speculative but a similar trial of supervised prolonged combined inhaled treatment is not unreasonable.

There are a number of strengths and weaknesses in this study. The main strength is that we have studied the largest group of preterm-born children with standardised objective assessments that has permitted the classification into different phenotypes. Although we have recently reported association of lung volumes with POLD and pPRISm,32 the majority of the children, especially those preterm-born children in the pPRISm group, did not have their lung volumes assessed to enable classification into restrictive or mixed lung disease. While we did not have total lung volumes for all the children, our primary aim was to identify specific phenotypes that can be identified in the general outpatient clinic where spirometry and FENO measurements are often routine. Although some of the data may be subject to recall bias, we have based the classification into phenotypes by using objective measures which are not influenced by recall bias. In addition, although we used the definition for dysanapsis suggested by the ATS/ERS guidelines, we are aware that there are a number of related definitions14 35 and agree further work is required to identify an ideal definition which is applicable to both children and adults.

In summary, we have identified four phenotypes of PLD which have differential associations with bronchodilators, FENO and with early life factors. These phenotypes especially the POLD group were recently shown to be associated with imaging abnormalities.36 The burden of respiratory disease is significant as 23% of the preterm population, compared with 9.2% of term-born children, was associated with at least one of these four phenotypes. By studying the mechanisms that underly these endotypes, we hope, specifically targeted therapies to prevent or improve respiratory outcomes of preterm-born children can be developed.

Data availability statement

Data from the RHiNO study are available to research collaborators subject to confidentiality and non-disclosure agreements. Contact Professor Sailesh Kotecha (kotechas{at}cardiff.ac.uk) for any data requests.

Ethics statements

Patient consent for publication

Ethics approval

Ethical approval was obtained from the South-West Bristol Research Ethics Committee (15/SW/0289). Written informed consent/assent were obtained from the parents/children, respectively.

Acknowledgments

We would like to thank the RHiNO team and to Louise Richards at NWIS for the information provided to contact the children. Most importantly, we would like to acknowledge and thank all the children and their parents who participated in the study for their great support and enthusiasm throughout the study.

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • MC and KH are joint first authors.

  • Contributors SK and JH conceived and designed the study. MC, KH and SK were involved in identifying and assessing the children and in data collection and its interpretation. SK, MC, KH, SJK, WJW and AB were involved in the data analysis and interpretation. MC, AB and SK drafted the manuscript. All authors were involved in revising the manuscript and approved the final submitted version.

  • Funding This study was funded by the Medical Research Council (reference: MR/M022552/1).

  • Competing interests SK reports securing a research grant from the Medical Research Council for this work. SK reports funding from HTA/NIHR, Moulton Foundation, GSK, Nutricia Foundation and Aspire Pharma outside this work. SJK and WJW reports funding from Moulton Foundation outside this work.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.