# Health-related quality of life of children born very preterm: a multinational European cohort study

Sung Wook Kim, Lazaros Andronis, Anna-Veera Seppänen, Adrien Aubert, Henrique Barros, Elizabeth Draper, Mariane Sentenac, Jennifer Zeitlin, Stavros Petrou

<figure><img src="image_1.png" /></figure>

To cite this version:

Sung Wook Kim, Lazaros Andronis, Anna-Veera Seppänen, Adrien Aubert, Henrique Barros, et al. "Health-related quality of life of children born very preterm: a multinational European cohort study." Quality of Life Research, 2023, 32 (1), pp.47-58. 10.1007/s11136-022-03217-9. hal-04556213

HAL Id: hal-04556213

https://ut3-toulouseinp.hal.science/hal-04556213v1

Submitted on 7 Oct 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Health-related quality of life of children born very preterm: a multinational European cohort study

Sung Wook Kim<fnref n="1" />☐, Lazaros Andronis<fnref n="2" />, Anna-Veera Seppänen<fnref n="3" />, Adrien M. Aubert<fnref n="3" />, Henrique Barros<fnref n="4" />,

Elizabeth S. Draper<fnref n="5" />, Mariane Sentenac<fnref n="3" />, Jennifer Zeitlin<fnref n="3" />, Stavros Petrou<fnref n="1" />, On Behalf of the SHIPS Research Group

Accepted: 24 July 2022 / Published online: 17 August 2022

© The Author(s) 2022

## Abstract

**Purpose** This study aims to (1) describe the health-related quality of life (HRQoL) outcomes experienced by children born very preterm (28–31 weeks’ gestation) and extremely preterm (< 28 weeks’ gestation) at five years of age and (2) explore the mediation effects of bronchopulmonary dysplasia (BPD) and severe non-respiratory neonatal morbidity on those outcomes.

**Methods** This investigation was based on data for 3687 children born at < 32 weeks’ gestation that contributed to the EPICE and SHIPS studies conducted in 19 regions across 11 European countries. Descriptive statistics and multi-level ordinary linear squares (OLS) regression were used to explore the association between perinatal and sociodemographic characteristics and PedsQL™ GCS scores. A mediation analysis that applied generalised structural equation modelling explored the association between potential mediators and PedsQL™ GCS scores.

**Results** The multi-level OLS regression (fully adjusted model) revealed that birth at < 26 weeks’ gestation, BPD status and experience of severe non-respiratory morbidity were associated with mean decrements in the total PedsQL™ GCS score of 0.35, 3.71 and 5.87, respectively. The mediation analysis revealed that the indirect effects of BPD and severe non-respiratory morbidity on the total PedsQL™ GCS score translated into decrements of 1.73 and 17.56, respectively, at < 26 weeks’ gestation; 0.99 and 10.95, respectively, at 26–27 weeks’ gestation; and 0.34 and 4.80, respectively, at 28–29 weeks’ gestation (referent: birth at 30–31 weeks’ gestation).

**Conclusion** The findings suggest that HRQoL is particularly impaired by extremely preterm birth and the concomitant complications of preterm birth such as BPD and severe non-respiratory morbidity.

**Keywords** Preterm birth; health-related quality of life · Mediation analysis · Children · PedsQL

The members of SHIPS Research Group are mentioned in
Acknowledgements.

<figure><img src="image_3.png" /></figure>

Sung Wook Kim
sungwook.kim@phc.ox.ac.uk

## Plain English summary

Preterm birth is associated with an increased risk of mor-
tality in comparison to birth at term. Compared to chil-
dren born at term, preterm born children are more likely to
develop serious complications such as bronchopulmonary
dysplasia (BPD) and severe non-respiratory morbidity,
including intraventricular haemorrhage, periventricular leu-
komalacia, necrotising enterocolitis and retinopathy of pre-
maturity. Nevertheless, evidence on the mechanism through
which effects of neonatal morbidities are transmitted to the
health-related quality of life (HRQoL) of those born preterm
is scarce. Hence, this study aims to investigate the mecha-
nism. Furthermore, this study aims to describe the HRQoL
outcomes experienced by children born preterm at five years
of age. The study used data for 3687 children obtained from

<footnote n="1">¹ Nuffield Department of Primary Care Health Sciences,
University of Oxford, Radcliffe Observatory Quarter,
Woodstock Road, Oxford OX2 6GG, UK</footnote>

<footnote n="2">² Warwick Clinical Trials Unit, Warwick Medical School,
University of Warwick, Coventry, UK</footnote>

<footnote n="3">³ Université Paris Cité, Inserm, INRAE, Centre for Research
in Epidemiology and StatisticS (CRESS), Obstetrical
Perinatal and Pediatric Epidemiology Research Team,
EPOPÉ, 75004 Paris, France</footnote>

<footnote n="4">⁴ EPIUnit-Instituto de Saúde Pública da Universidade do Porto,
Porto, Portugal</footnote>

<footnote n="5">⁵ Department of Health Sciences, University of Leicester,
Leicester, UK</footnote>
the EPICE (Effective Perinatal Intensive care in Europe) and SHIPS (Screening to improve Health In very Preterm infantS) studies. The data were collected from 11 European countries. Data were analysed using a multi-level ordinary least squares and a separate mediation analysis. The new evidence provided by this study shows that neonatal morbidities and, in particular severe non-respiratory morbidity, are the main drivers of HRQoL outcomes in children born very preterm and extremely preterm. These results can be used to inform service targeting towards these higher risk groups.

## Introduction

Very preterm birth (VPT, birth at 28 to 31 weeks' gestation) and extremely preterm birth (EPT, birth at <28 weeks' gestation) are associated with an increased risk of mortality in comparison to birth at term [1]. Moreover, although survival rates following VPT and EPT have increased in recent years [2], surviving children are at increased risk of cerebral palsy, visual and auditory deficits, poor respiratory outcomes, impaired motor and cognitive ability and psychiatric disorders compared to children born at term [3–6]. Furthermore, up to one third of children born VPT or EPT and their parents face a life course with significant morbidity, dependency, low educational achievement and socioeconomic challenges [7–9].

There is increasing recognition of the need to measure the impact of preterm birth across multiple domains [10]. Assessments of health-related quality of life (HRQoL) outcomes provide valuable complementary data to biomedical assessments that have traditionally been reported [10]. Importantly, they also provide a framework for reflecting the multi-dimensional aspects of physical, mental and social well-being of those born preterm and their families. Several assessments of HRQoL following preterm birth have applied multi-dimensional health profile measures that disaggregate outcomes across dimensions or subscales [11–13]. A small number of recent assessments of the HRQoL of preterm-born children [14–16] have applied the Pediatric Quality of life Inventory™ (PedsQL™) [17], which was designed to provide an age-specific modular approach to measuring HRQoL in healthy children and adolescents, as well as those with acute and chronic health conditions, across the broadest empirically feasible age groups. However, these studies are based on case series or cohorts with relatively small samples, with each limited to application in one country.

It is well known that, compared to children born at term, preterm born children are more likely to develop serious complications such as bronchopulmonary dysplasia (BPD) [18, 19] and severe non-respiratory morbidity [20–22]. Those who develop BPD can suffer from long-term respiratory complications, such as impaired lung function, which

can have life-long consequences for those born preterm [18]. Likewise, the severe non-respiratory morbidity associated with preterm birth can also result in decrements in the HRQoL of those born preterm. For example, one study found that pre-school children born preterm who had been admitted to neonatal intensive care due to morbidity associated with preterm birth had significantly poorer HRQoL than children born at term across the measurement scales for lungs, stomach, eating disorders, motor functioning and communication of the TNO-AZL Preschool Quality Of Life questionnaire [23]. However, evidence on the mediation effects of BPD and severe non-respiratory morbidity on the HRQoL of those born preterm is scarce. More generally, studies are required to understand the mechanisms through which preterm birth affects the HRQoL of those born very or extremely preterm. The aims of this study are therefore two-fold: firstly, to describe the HRQoL experienced by children born VPT or EPT across multi-dimensional outcomes and, secondly, to explore the mediation effects of BPD and severe non-respiratory morbidity on those outcomes.

## Methods

### Study population

Data for this study were obtained from the EPICE (Effective Perinatal Intensive care in Europe) cohort study and SHIPS (Screening to improve Health In very Preterm infantS), a follow-up study for EPICE. A detailed description of the design, outcome measures and data collection processes in the EPICE and SHIPS studies can be found elsewhere [24]. In brief, the EPICE and SHIPS studies constituted and followed up a population-based prospective cohort of all children born between 22+0 weeks and 31+6 weeks of gestation over a 12-month period (except in France where it was conducted over a 6-month period) in 2011–2012 in 19 regions across 11 European countries: Belgium (Flanders); Denmark (Eastern Region); Estonia (entire country); France (Burgundy, Ile-de-France and the Northern region); Germany (Hesse and Saarland); Italy (Emilia-Romagna, Lazio and Marche); the Netherlands (Central and Eastern region), Poland (Wielkopolska); Portugal (Lisbon and Northern region); Sweden (Greater Stockholm) and the United Kingdom (UK) (East Midlands, Northern and Yorkshire & the Humber regions). Perinatal data were extracted from medical records whilst parental questionnaires at two and five years of age were used as the primary vehicle for assessments of longer term outcomes.
## Health-related quality of life

To measure HRQoL outcomes, the PedsQL™ 4.0 Generic Core Scales (hereafter PedsQL GCS for brevity) were completed via postal questionnaires by parents when their children were five years of age. The parent proxy-report [ages 5–7 (young child)] format was used. The PedsQL GCS have been developed to measure the core dimensions of health as well as school functioning [25]. The PedsQL GCS contain 23 items covering physical functioning (8 items), emotional functioning (5 items), social functioning (5 items) and school/day care functioning (5 items). A 5-point response scale is utilised across each item (0 = never a problem; 1 = almost never a problem; 2 = sometimes a problem; 3 = often a problem; 4 = almost always a problem). Raw responses were scored (0 to 4) and reversely transformed to a 0–100 scale. The mean of scale scores was computed as the sum across the relevant items over the number of items answered, thereby accounting for missing data [25]. The mean psychosocial health summary score was calculated as the sum of the item scores over the number of the items completed for the emotional, social and school functioning scales, whilst the physical health summary score was identical to the physical functioning scale score. The total score represents the sum of scores across all 23 items divided by the number of items answered. Higher scores indicate better HRQoL.

## Predictive factors

The key predictive factors incorporated into the analysis that explored mediation effects were gestational age at birth, BPD status and severe non-respiratory morbidity status, all of which were measured using data extracted from medical records. Gestational age at birth was classified into four categories of < 26 weeks, 26 to 27 weeks, 28 to 29 weeks and 30 to 31 weeks. BPD was defined as receipt of supplemental oxygen and/or ventilatory support (continuous positive airway pressure (CPAP) or mechanical ventilation) at 36 weeks of postmenstrual age and dichotomised by its presence or absence. Severe non-respiratory morbidity consisted of intraventricular haemorrhage (IVH) grade III/IV, cystic periventricular leukomalacia (cPVL), retinopathy of prematurity (ROP) stage III or more or necrotizing enterocolitis (NEC) requiring surgery or peritoneal drainage and dichotomised according to the presence or absence of at least one of these morbidities.

## Statistical analysis

Differences in the total PedsQL GCS score in relation to preterm status (VPT vs EPT), and presence of BPD (no, yes) or severe non-respiratory morbidity (no, yes), were tested using

Student's *t* test. These differences were estimated for each of the scales of the PedsQL GCS, namely physical functioning, emotional functioning, social functioning, school functioning, psychosocial functioning, as well as for the total score.

Multi-level regression analysis was initially carried out to explore the association between the total PedsQL GCS score and perinatal, neonatal and sociodemographic variables. As the overall PedsQL GCS score is a continuous variable, ordinary least square (OLS) regression was employed in the analysis. The following levels of the multi-level analysis were specified: 1) individual child, 2) parent/caregiver and 3) country. Independent variables in the model were chosen based on published evidence [26–28], including published empirical models using data from the EPICE and SHIPS studies [29]. These included the following: gestational age at birth (weeks) (< 26, 26–27, 28–29, 30–31) [30], small for gestational age (SGA) status (birth weight < 10th, ≥ 10th centile for gestational age) [31], congenital anomaly status (no or yes), sex (male or female), multiplicity (singletons, twins or more, based on known associations between multiplicity and increased risk of preterm birth) [30], BPD status (no or yes), severe non-respiratory morbidity status (no or yes), highest maternal education level (high school or less, more than high school) [27], maternal country of birth (native-born, non-native from other European country, non-native from non-European country), maternal age at child-birth (years) (< 25, 25 to 34, > 34) and parity (primiparous, multiparous).

Two regression models were constructed for the multi-level OLS regression: (1) a partially adjusted model that excluded BPD status and severe non-respiratory morbidity status and (2) a fully adjusted model that included all independent variables. Due to different education systems across the participating countries with children starting compulsory schooling at different ages, regions in Denmark and Germany provided a low number of responses to questions about school functioning within the PedsQL GCS. To address this issue, a sensitivity analysis excluding Denmark and Germany was performed.

BPD and severe non-respiratory morbidity can act as mediators that negatively affect health for preterm born children [28, 32]. Mediation analysis is a widely used tool for exploring whether the effect of an independent variable on a dependent variable is entirely or partially due to one or more intermediary (mediator) variables [33, 34]. In mediation analysis, the overall effect can be decomposed into direct and indirect effects. Total effects are calculated as the sum of the direct effects and indirect effects that are the product of path coefficients (e.g. preterm birth to BPD * BPD to PedsQL GCS score). This can be used to understand a relationship between an independent variable (X) and a dependent variable (Y) with the underlying mechanism that affects the dependent variable through a mediator variable (M). Simply
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Overview of mediation analysis</figcaption></figure>

adjusting for the intermediate variable can cause bias [35]. As a result, in a regression analysis, the coefficient of X on Y can be considered a direct effect, whereas the coefficient of M on Y can be considered an indirect effect. In cases where there are multiple mediators in the model, these can be estimated using either Structural Equation Modelling (SEM) or Generalised Structural Equation Modelling (GSEM) [36]. Direct effects indicate the effects of independent variables on a dependent variable whilst indirect effects indicate the effects that are mediated through the mediators on the dependent variable (Fig. 1). In general, if the indirect effects are statistically significant, it is believed that mediation effects exist [37]. In other words, the significance of the total effect is not necessary for the existence of mediation effects. In this study, the mediating variables were BPD and severe non-respiratory morbidity. Direct effects were estimated as the effects of gestational age per se on the total PedsQL GCS score whilst indirect effects were estimated as the effects of the mediating variables of BPD and severe non-respiratory morbidity on the total PedsQL GCS score. This is illustrated in Fig. 1. The GSEM framework was applied and standard errors were estimated using a bootstrap procedure run over 1000 iterations [38]. A sensitivity analysis that excluded data from Denmark and Germany was also conducted. Statistical analyses were carried out with using Stata 16.0 (StataCorp, College Station, TX).

## Results

Table 1 summarises the characteristics of the study children at birth, split by country. Estonia and the UK had the lowest (N=153) and highest numbers (N=1745) of preterm born study children, respectively. The number of children followed up at five years, using the denominator of all live and eligible children (N=6759), was 3687 (54.5%). Of the eligible children, singletons represented 5437 (80.4%) of the preterm born children, whereas twins and higher order births represented 1322 (19.6%) of the preterm born children.

<table><caption>**Table 1** Descriptive statistics of participant characteristics (entire cohort of live births), by country</caption><thead><tr><th></th><th>Belgium</th><th>Denmark</th><th>Estonia</th><th>France</th><th>Germany</th><th>Italy</th><th>Netherlands</th><th>Poland</th><th>Portugal</th><th>UK</th><th>Sweden</th></tr></thead><tbody><tr><td>Infants (N)</td><td>752</td><td>351</td><td>153</td><td>1307</td><td>758</td><td>1134</td><td>393</td><td>316</td><td>724</td><td>1745</td><td>267</td></tr><tr><td>Gestational age, completed weeks, Mean (SD)</td><td>28.6 (2.3)</td><td>28.0 (2.4)</td><td>28.7 (2.2)</td><td>28.4 (2.3)</td><td>28.2 (2.6)</td><td>28.5 (2.4)</td><td>28.4 (2.5)</td><td>28.1 (2.8)</td><td>28.5 (2.2)</td><td>28.5 (2.4)</td><td>28.2 (2.5)</td></tr><tr><td>Male, N (%)</td><td>426 (56.6)</td><td>185 (52.7)</td><td>80 (52.3)</td><td>701 (53.6)</td><td>407 (53.7)</td><td>563 (49.6)</td><td>195 (49.6)</td><td>176 (55.7)</td><td>423 (58.4)</td><td>965 (55.3)</td><td>145 (54.3)</td></tr><tr><td>Singleton, N (%)</td><td>464 (61.7)</td><td>208 (59.3)</td><td>119 (77.8)</td><td>881 (67.4)</td><td>487 (64.2)</td><td>768 (67.7)</td><td>274 (69.7)</td><td>251 (79.4)</td><td>500 (69.1)</td><td>1300 (74.6)</td><td>185 (69.3)</td></tr><tr><td>Twins or more, N (%)</td><td>288 (38.3)</td><td>143 (40.7)</td><td>34 (22.2)</td><td>426 (32.6)</td><td>271 (35.8)</td><td>366 (32.3)</td><td>119 (30.3)</td><td>65 (20.6)</td><td>224 (30.9)</td><td>443 (25.4)</td><td>82 (30.7)</td></tr><tr><td>Birth weight (grams), Mean (SD)</td><td>1260.0 (420.8)</td><td>1141.6 (415.4)</td><td>1287.0 (390.1)</td><td>1162.3 (372.1)</td><td>1165.1 (418.0)</td><td>1195.4 (420.3)</td><td>1198.2 (412.1)</td><td>1224.7 (474.2)</td><td>1149.4 (363.5)</td><td>1210.4 (383.9)</td><td>1216.7 (436.4)</td></tr><tr><td>Follow up at five years of age, N (%)</td><td>280 (37.2)</td><td>152 (43.3)</td><td>134 (87.6)</td><td>779 (59.6)</td><td>280 (36.9)</td><td>693 (61.1)</td><td>155 (39.4)</td><td>189 (59.8)</td><td>433 (59.8)</td><td>448 (25.7)</td><td>144 (53.9)</td></tr></tbody></table>

UK denotes United Kingdom; SD denotes standard deviation

<figure><img src="image_5.png" /></figure>
Estonia had the highest follow-up rate (88%) whilst the UK had the lowest follow-up rate (26%).

Table 2 presents differences in mean PedsQL GCS scores by preterm status, BPD status and severe non-respiratory morbidity status. The results show that the total PedsQL GCS score differed significantly by preterm status (78.61 for VPT vs 75.70 for EPT; *p* < 0.001), BPD status (78.60 for no BPD vs 72.70 for BPD; *p* < 0.001) and severe non-respiratory morbidity status (78.59 for no morbidity vs 71.98 for morbidity; *p* < 0.001). Significant differences were also observed in PedsQL GCS physical, emotional, social, school and psychosocial functioning scores by preterm status, BPD status and severe non-respiratory morbidity status. A table that replicates these analyses for the four gestational age categories (< 26, 26–27, 28–29, 30–31 weeks) is presented in Online Appendix 1.

Table 3 shows the results of the OLS regression based on multi-level modelling. Both partially and fully adjusted models present similar patterns in terms of the direction of signs, statistical significance and magnitude of coefficients. In the fully adjusted model, the total PedsQL GCS score of the children who were born at < 26 gestational weeks, with BPD or with severe non-respiratory morbidity were lower than the respective values for the reference group of children born at 30–31 gestational weeks, without BPD or without

severe non-respiratory morbidity by 0.35 (*p* = 0.834), 3.71 (*p* < 0.001) and 5.87 (*p* < 0.001), respectively. Birth earlier than 26 weeks’ gestation was associated with a significant reduction in the total PedsQL GCS score in the partially adjusted model (*p* < 0.05), but not following further adjustment for BPD status and severe non-respiratory morbidity status. In addition, congenital anomaly status, sex, maternal education and maternal country of birth all had statistically significant effects in both the partially and fully adjusted models. For instance, in the fully adjusted model, the coefficients for presence of a congenital anomaly, male sex, highest maternal education of high school or less and non-European maternal country of birth were – 3.34 (*p* < 0.001), – 2.86 (*p* < 0.001), – 1.56 (*p* < 0.05) and – 6.45 (*p* < 0.001), respectively. The results of the sensitivity analysis excluding Denmark and Germany showed a similar pattern and are presented in Online Appendix 2. The results of a multi-level OLS regression that treated gestational age as a continuous variable is presented in Online Appendix 3.

Table 4 shows the results of the mediation analysis using BPD and severe non-respiratory morbidity as mediators based on the GSEM (Online Appendix 4). The indirect effects of BPD and severe non-respiratory morbidity on the total PedsQL GCS score were estimated at – 1.73 (bootstrap 95% confidence interval (CI) – 2.68 to – 0.77) and

<table><caption>**Table 2** Differences in PedsQL GCS scores by preterm, BPD and severe non-respiratory morbidity status</caption><thead><tr><th>Preterm status</th><th>VPT, mean (SD)</th><th>EPT, mean (SD)</th><th>Mean difference</th><th>95% Confidence interval</th><th>p value</th></tr></thead><tbody><tr><td>Physical functioning</td><td>81.34 (20.39)</td><td>77.82 (22.40)</td><td>3.41</td><td>(1.69, 5.13)</td><td>&lt;0.001</td></tr><tr><td>Emotional functioning</td><td>76.19 (16.70)</td><td>74.83 (17.77)</td><td>1.09</td><td>(-0.30, 2.48)</td><td>0.037</td></tr><tr><td>Social functioning</td><td>82.74 (18.73)</td><td>79.03 (20.20)</td><td>3.26</td><td>(1.71, 4.81)</td><td>&lt;0.001</td></tr><tr><td>School functioning</td><td>76.51 (19.27)</td><td>72.91 (19.66)</td><td>3.64</td><td>(2.06, 5.21)</td><td>&lt;0.001</td></tr><tr><td>Psychosocial functioning</td><td>79.33 (14.90)</td><td>76.27 (16.19)</td><td>2.66</td><td>(1.60, 4.10)</td><td>&lt;0.001</td></tr><tr><td>Total score</td><td>78.61 (14.93)</td><td>75.70 (16.06)</td><td>2.85</td><td>(1.42, 3.90)</td><td>&lt;0.001</td></tr><tr><td>BPD</td><td>No, mean (SD)</td><td>Yes, mean (SD)</td><td>Mean difference</td><td>95% Confidence interval</td><td>p value</td></tr><tr><td>Physical functioning</td><td>81.27 (20.33)</td><td>74.41 (24.40)</td><td>6.86</td><td>(4.86, 8.86)</td><td>&lt;0.001</td></tr><tr><td>Emotional functioning</td><td>76.15 (16.65)</td><td>73.23 (18.91)</td><td>2.92</td><td>(1.29, 4.55)</td><td>0.001</td></tr><tr><td>Social functioning</td><td>82.76 (18.51)</td><td>75.14 (21.87)</td><td>7.62</td><td>(5.79, 9.45)</td><td>&lt;0.001</td></tr><tr><td>School functioning</td><td>76.61 (18.85)</td><td>69.10 (21.34)</td><td>7.51</td><td>(5.54, 9.48)</td><td>&lt;0.001</td></tr><tr><td>Psychosocial functioning</td><td>79.31 (14.78)</td><td>72.97 (17.99)</td><td>6.35</td><td>(4.89, 7.80)</td><td>&lt;0.001</td></tr><tr><td>Total score</td><td>78.60 (14.79)</td><td>72.70 (17.26)</td><td>5.90</td><td>(4.45, 7.35)</td><td>&lt;0.001</td></tr><tr><td>Severe non-respiratory morbidity</td><td>No, mean (SD)</td><td>Yes, mean (SD)</td><td>Mean difference</td><td>95% Confidence interval</td><td>p value</td></tr><tr><td>Physical functioning</td><td>81.52 (19.84)</td><td>71.43 (26.85)</td><td>10.1</td><td>(7.83, 12.37)</td><td>&lt;0.001</td></tr><tr><td>Emotional functioning</td><td>76.24 (16.85)</td><td>73.32 (17.59)</td><td>2.93</td><td>(1.06, 4.79)</td><td>0.002</td></tr><tr><td>Social functioning</td><td>82.70 (18.64)</td><td>74.20 (21.74)</td><td>8.51</td><td>(6.41, 10.60)</td><td>&lt;0.001</td></tr><tr><td>School functioning</td><td>76.44 (18.89)</td><td>68.04 (21.95)</td><td>8.40</td><td>(6.14, 10.67)</td><td>&lt;0.001</td></tr><tr><td>Psychosocial functioning</td><td>79.37 (14.70)</td><td>71.83 (18.00)</td><td>7.54</td><td>(5.88, 9.19)</td><td>&lt;0.001</td></tr><tr><td>Total score</td><td>78.59 (14.84)</td><td>71.98 (17.38)</td><td>6.62</td><td>(4.96, 8.28)</td><td>&lt;0.001</td></tr></tbody></table>

Differences in characteristics between preterm status, BPD and severe non-respiratory morbidity status were tested using the Student t test. BPD denotes bronchopulmonary dysplasia; VPT denotes very preterm; EPT denotes extremely preterm
<table><caption>**Table 3** Multi-level regression analysis of the association between the total PedsQL GCS score and clinical and sociodemographic factors</caption><thead><tr><th></th><th></th><th>Partially adjusted<sup>+</sup></th><th>Fully adjusted</th></tr></thead><tbody><tr><td rowspan="3">Gestational age (weeks) (ref: 30–31)</td><td>&lt; 26</td><td>− 4.46 (2.22)*</td><td>− 0.35 (1.64)</td></tr><tr><td>26–27</td><td>− 1.46 (0.83)</td><td>0.42 (0.71)</td></tr><tr><td>28–29</td><td>− 0.68 (0.89)</td><td>− 0.23 (0.82)</td></tr><tr><td>SGA (ref: ≥ 10th centile)</td><td>&lt; 10th centile</td><td>− 0.74 (0.60)</td><td>− 0.39 (0.57)</td></tr><tr><td>Congenital anomalies (ref: no)</td><td>Yes</td><td>− 3.95 (0.87)***</td><td>− 3.34 (0.90) ***</td></tr><tr><td>Sex (ref: female)</td><td>Male</td><td>− 2.93 (0.69)***</td><td>− 2.86 (0.69) ***</td></tr><tr><td>Multiplicity (ref: singleton)</td><td>Twins and higher</td><td>2.02 (0.83)*</td><td>1.83 (0.87) *</td></tr><tr><td rowspan="2">BPD (ref: no)</td><td>Yes</td><td></td><td>− 3.71 (0.92) ***</td></tr><tr><td>Yes</td><td></td><td>− 5.87 (1.13) ***</td></tr><tr><td>Mother's education (ref: higher education)</td><td>High school or less</td><td>− 1.71 (0.62)**</td><td>− 1.56 (0.66)*</td></tr><tr><td rowspan="2">Country of birth for mothers (ref: native)</td><td>Non-native, European born</td><td>− 0.18 (1.24)</td><td>0.59 (1.29)</td></tr><tr><td>Non-native, non-European born</td><td>− 5.81 (0.99)***</td><td>− 6.45 (1.02)***</td></tr><tr><td rowspan="2">Mother's age at childbirth (years) (ref: 25–34)</td><td>&lt; 25</td><td>− 1.07 (0.99)</td><td>− 0.78 (0.95)</td></tr><tr><td>&gt; 34</td><td>1.31 (1.03)</td><td>1.33 (1.09)</td></tr><tr><td>Parity (ref: multiparous)</td><td>Primiparous</td><td>− 0.26 (0.65)</td><td>− 0.20 (0.55)</td></tr></tbody></table>

SGA small for gestational age, BPD bronchopulmonary dysplasia. This is defined as receipt of supplemental oxygen and/or ventilatory support (CPAP or mechanical ventilation) at 36 weeks of postmenstrual age

Severe non-respiratory morbidity: Intraventricular haemorrhage grades III–IV (IVH), periventricular leukomalacia (PVL), retinopathy of prematurity stages III-V (ROP) or necrotising enterocolitis needing surgery (NEC)

Likelihood Ratio (LR) test for linear regression = *p* < 0.001

– 17.56 (bootstrap CI – 24.15 to – 10.96), respectively, at < 26 weeks’ gestation (referent: birth at 30–31 weeks’ gestation). The indirect effects of BPD and severe non-respiratory morbidity on the total PedsQL GCS score were estimated at – 0.99 (bootstrap 95% CI – 1.53 to – 0.46) and – 10.95 (bootstrap 95% CI – 15.19 to – 6.70), respectively, at 26–27 weeks’ gestation, and – 0.34 (bootstrap 95% CI – 0.54 to – 0.15) and – 4.80 (bootstrap 95% CI – 6.96 to – 2.64), respectively, at 28–29 weeks’ gestation (referent: birth at 30–31 weeks’ gestation). In contrast, the direct effects of birth at < 26 weeks, 26–27 weeks and 28–29 weeks were – 1.06 (bootstrap 95% CI – 3.38 to 1.26), 0.10 (bootstrap 95% CI – 1.37 to 1.57) and – 0.45 (bootstrap 95% CI – 1.73 to 0.82), respectively (referent: birth at 30–31 weeks’ gestation), which were not statistically significant. In short, Table 4 illustrates the importance of the indirect effects of BPD and severe non-respiratory morbidity on the total PedsQL GCS score, with diminishing indirect effects with increasing gestational age. This result is graphically presented in Fig. 2. Online Appendix 5 presents the sensitivity analysis that excluded Denmark and Germany from the mediation analysis and reveal similar magnitude of coefficients observed when the data for the entire study population were analysed.

Discussion

This study provides new evidence on HRQoL outcomes in mid-childhood of children born very or extremely preterm in 11 European countries. To the best of our knowledge, this is the first study that provides information with respect to HRQoL outcomes, as measured by the PedsQL GCS, of children born preterm across several European countries. The estimated PedsQL GCS scores are based on 3687 VPT born children in Europe. In addition, this study explores the mediation effects of BPD and severe non-respiratory morbidity on PedsQL GCS scores.

The study confirms that the HRQoL of VPT born children is significantly better in those who did not experience BPD or severe non-respiratory morbidity. This result is consistent with evidence from existing studies [28, 39, 40]. A study conducted in the USA [41] reported that in children who survived to 18–36 months of corrected age with severe BPD, the mean physical health summary score of the PedsQL GCS was 78.0 and the mean psychosocial health summary score of the PedsQL GCS was 75.3. This compared with physical health and psychosocial health summary scores of 74.4 and 73.0, respectively, observed

<footnote>*Partially adjusted model excludes BPD and severe non-respiratory morbidity</footnote>

<footnote>* *p* < 0.05, ** *p* < 0.01, *** *p* < 0.001</footnote>
<table><caption>**Table 4** Mediation analysis of the indirect effects and total effects of preterm status on the total PedsQL GCS score (four categories of gestational age: <26 weeks, 26–27 weeks; 28–29 weeks; 30–31 weeks (referent))</caption><thead><tr><th></th><th></th><th>Coef (SE)</th><th>Bootstrap 95% confidence interval</th></tr></thead><tbody><tr><td rowspan="6">&lt;26</td><td>BPD: indirect effects<sup>a</sup></td><td>-1.73 (0.49)</td><td>(-2.68, -0.77)</td></tr><tr><td>Direct effects<sup>b</sup></td><td>-1.06 (1.18)</td><td>(-3.38, 1.26)</td></tr><tr><td>BPD: total effects<sup>c</sup></td><td>-2.79 (1.16)</td><td>(-5.06, -0.51)</td></tr><tr><td>Severe non-respiratory morbidity: indirect effects<sup>a</sup></td><td>-17.56 (3.37)</td><td>(-24.15, -10.96)</td></tr><tr><td>Direct effects<sup>b</sup></td><td>-1.06 (1.18)</td><td>(-3.38, 1.26)</td></tr><tr><td>Severe non-respiratory morbidity: total effects<sup>c</sup></td><td>-18.62 (3.29)</td><td>(-25.06, -12.17)</td></tr><tr><td rowspan="6">26-27</td><td>BPD: indirect effects<sup>a</sup></td><td>-0.99 (0.27)</td><td>(-1.53, -0.46)</td></tr><tr><td>Direct effects<sup>b</sup></td><td>0.10 (0.75)</td><td>(-1.37, 1.57)</td></tr><tr><td>BPD: total effects<sup>c</sup></td><td>-0.90 (0.75)</td><td>(-2.37, 0.58)</td></tr><tr><td>Severe non-respiratory morbidity: indirect effects<sup>a</sup></td><td>-10.95 (2.16)</td><td>(-15.19, -6.70)</td></tr><tr><td>Direct effects<sup>b</sup></td><td>0.10 (0.75)</td><td>(-1.37, 1.57)</td></tr><tr><td>Severe non-respiratory morbidity: total effects<sup>c</sup></td><td>-10.85 (2.17)</td><td>(-15.09, -6.60)</td></tr><tr><td rowspan="6">28-29</td><td>BPD: indirect effects<sup>a</sup></td><td>-0.34 (0.10)</td><td>(-0.54, -0.15)</td></tr><tr><td>Direct effects<sup>b</sup></td><td>-0.45 (0.65)</td><td>(-1.73, 0.82)</td></tr><tr><td>BPD: total effects<sup>c</sup></td><td>-0.80 (0.65)</td><td>(-2.08, 0.49)</td></tr><tr><td>Severe non-respiratory morbidity: indirect effects<sup>a</sup></td><td>-4.80 (1.10)</td><td>(-6.96, -2.64)</td></tr><tr><td>Direct effects<sup>b</sup></td><td>-0.45 (0.65)</td><td>(-1.73, 0.82)</td></tr><tr><td>Severe non-respiratory morbidity: total effects<sup>c</sup></td><td>-5.25 (1.24)</td><td>(-7.68, -2.82)</td></tr></tbody></table>

This was performed based on 1000 bootstrap simulations

<figure><img src="image_6.png" /><figcaption>**Fig. 2** Mediation analysis of the indirect effects and total effects of preterm status on the total PedsQL GCS score*</figcaption></figure>

<footnote><sup>a</sup>This presents the indirect effect of each mediator (BPD and Severe non-respiratory morbidity) on the PedsQL GCS score. See Fig. 1</footnote>

<footnote><sup>b</sup>This represents the direct effect of gestational age (GA) on the PedsQL GCS score. See Fig. 1</footnote>

<footnote><sup>c</sup>Total effects are calculated as the sum of the direct effects and indirect effects that are the product of path coefficients (e.g. preterm birth to BPD * BPD to PedsQL GCS score)</footnote>
amongst children with BPD in our study (descriptive analyses). Although it is not straightforward to compare and contrast our PedsQL GCS score findings with those in other studies and with ‘population’ norms, a study conducted in Canada [42] for liver transplant children aged under 6 found that healthy children had a mean total PedsQL score of 82, compared to 73 and 75 for chronically ill children and liver transplant children, respectively.

It is unclear whether the approximate 3 point difference in mean total PedsQL GCS that we observed between children born VPT and EPT can be considered clinically meaningful to stakeholders as there is, to our knowledge, no generally accepted minimum clinically important difference in the total PedsQL GCS score. However, there have been a few studies that have generated estimates of minimally important clinical differences for the PedsQL GCS score [42–44]. Amongst these studies, a US study [43] that estimated PedsQL GCS scores of 20,031 families with children aged 2–16 years showed that healthy children had a mean total PedsQL GCS score of 84, compared to 74 for chronically ill children. Furthermore, as noted above, a Canadian study [42] of liver transplant children aged under 6 years found that healthy children had a mean total PedsQL GCS score of 82, compared to 73 and 75 for chronically ill children and liver transplant children, respectively.

In our partially adjusted multi-level OLS model, birth at <26 weeks' gestation was associated with a statistically significant decrement in HRQoL, consistent with evidence from other research studies [40, 45, 46], but this effect dissipated after adjustment for BPD status and severe non-respiratory morbidity status (fully adjusted model). Our multi-level OLS model also found that presence of a congenital anomaly, male sex, low maternal education status and non-European maternal country of birth are independent predictors of the HRQoL of VPT born children. Comparative evidence of the association between these factors and the HRQoL of VPT born children is sparse. One study conducted in Finland found that presence of a congenital anomaly was negatively associated with the HRQoL of preterm born children born between 2000 and 2003 [47]. A study of early-school age children born preterm, conducted in France, revealed that poorer HRQoL was reported for boys than for girls when the ‘Vécu et Santé Perçue de l’Adolescent et de l’Enfant’ (VSP-A) questionnaire (parent version) was used, but the difference was not statistically significant [48]. Low parental education status was also negatively associated with the HRQoL of preterm born children in a Taiwanese study, but the association was not statistically significant [49]. We are not, however, aware of any other studies that have explored the association between maternal country of birth and the HRQoL of preterm born children, although the association between migrant status and preterm birth rates is well established [50–52].

This study provides clear evidence around the mediation effects of BPD and severe non-respiratory morbidity on the HRQoL outcomes of VPT born children in European countries. Although there have been numerous studies describing the negative effects of perinatal complications on HRQoL outcomes, the mechanisms through which preterm birth translates into poorer HRQoL have not been elucidated. Existing studies [32, 41, 53] focus on either the association between complications and HRQoL outcomes or the association between preterm birth per se and HRQoL outcomes [40] rather than exploring or focussing on the mediating mechanisms through which those outcomes are experienced.

The mediation analysis performed in this study shows that an association exists between the complications associated with preterm birth including BPD and, in particular, severe non-respiratory morbidity which explains a high proportion of the total effect of gestational age and total PedsQL GCS scores. It was found that gestational age at birth no longer had a continued direct effect on outcomes after these indirect effects were taken into consideration as the direct effects of gestational age were not statistically significant. This implies that poorer HRQoL as measured by the total PedsQL GCS score is primarily driven by the complications of preterm birth rather than preterm status per se. The negative effects of BPD and severe non-respiratory morbidity for the health and development of preterm born children have been well reported, with severe non-respiratory morbidities, namely brain lesions, being strongly associated with some of the most severe impairments, such as cerebral palsy [54]. Further, one study reported that, within a preterm group, infants who had neonatal steroid treatment for BPD presented lower HRQoL scores as measured by the Child Health Questionnaire (CHQ-Pf 50) than infants who did not have this treatment [40].

Mediation analysis is, in principle, a causal inference model [55]. However, a mediation analysis proposes the causal relationship between the independent variable and the dependent variable by mediating variables rather than directly showing the causal relationship between the independent variable and the dependent variable. It has been argued that statistical significance only confirms the effect that the presumed model provides [56]. The strength of a mediation analysis comes from the appropriate study design based on the clinical and epidemiological context. Given this, the results of this study should be carefully interpreted because mediation is not defined statistically.

The main strength of this study comes from its large sample size of approximately 3600 born very preterm across Europe. The population-based sample was drawn from defined geographic areas rather than clinic-based populations; consequently, selection biases are unlikely to represent a major problem. It also used a validated generic measure of children’s HRQoL that facilitates comparisons with the

<figure><img src="image_7.png" /></figure>
HRQoL outcomes experienced by general childhood and clinical populations [48, 57]. Notably, a comprehensive analytical strategy was adopted to disentangle the direct and indirect effects of gestational age at birth on HRQoL outcomes.

Inevitably, there are limitations to this study. Firstly, the methods for collecting HRQoL data were heterogeneous across countries. For example, Denmark and Germany provided substantially lower numbers of observations due to the different schooling systems in those countries as registration for schooling starts later. To overcome this, a sensitivity analysis excluding these two countries was performed, the results of which showed a similar pattern. Secondly, the study is based on the untestable assumption of no omitted variables for a causal mediation analysis [58]. In order to minimise this potential bias, we adjusted for a broad range of covariates that can act as potential confounders in the multilevel OLS regression. Thirdly, this study was performed using parental reports of children's HRQoL and there may be biases associated with parental reporting. There are a number of studies with respect to cross reporting of children's HRQoL from both parents and children [59, 60], but they are mostly concentrated in children who can accurately express their thoughts and feelings (e.g. 8–16 years old) as pre-school aged children's communicative expressions are generally constrained. In addition, a recent study revealed that five-year-old children could not independently provide PedsQL GCS data of sufficient psychometric quality [61]. As such, this is not a problem unique to this study.

Our findings offer insights into factors that affect the HRQoL of children born very or extremely preterm. We anticipate that these insights will be useful input for policy making. By understanding which physical or mental ailments influence the HRQoL that VPT born children experience, researchers, clinicians and decision-makers should be able to concentrate their efforts on proactively targeting children with those particular morbidities. In terms of the health research agenda, the data presented in this study can be used in applied work. For example, our findings should be of help to researchers aiming to reflect the differences in HRQoL associated with different characteristics of preterm born children in prognostic models. Our findings can further be used to target services towards higher risk groups with severe morbidities rather than all children born very preterm. The data reported in this study can also act as a significant new resource that can inform HRQoL estimation within the context of complications associated with preterm birth.

In conclusion, this study describes the magnitude and mechanisms of the HRQoL outcomes of children born very preterm based on a large cohort recruited from 11 European countries. Equally importantly, the study's findings give rise to further questions that could trigger subsequent research on the topic. Based on our findings, further research is needed

to further clarify the mechanisms through which different morbidities affect the HRQoL of children born very preterm.

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1007/s11136-022-03217-9.

**Acknowledgements** We thank all the participants that completed the parental questionnaires. All members of the SHIPS study group reviewed and approved the final version for publication. SHIPS Research Group: Belgium: (J Lebeer¹, P Van Reempts², E Bruneel³, E Cloet⁴, A Oostra⁵, E Ortibus⁶, I Sarrechia⁷); Denmark (K Boerch⁸, P Pedersen⁹); Estonia (L Toome¹⁰, H Varendi¹¹, M Männamaa¹²); France (PY Ancel¹³, A Burguet¹⁴, PH Jarreau¹⁵, V Pierrat¹⁶, P Truffert¹⁷); Germany (RF Maier¹⁸, M Zemlin¹⁹, B Misselwitz²⁰, L Wohlers¹⁸, Italy (M Cuttini²¹, I Croci²¹, V Carnielli²², G Ancora²³, G Faldella²⁴, F Ferrari²⁵); The Netherlands (A van Heijst²⁶, C Koopman-Esseboom²⁷); Poland (J Gadzinowski²⁸, J Mazela²⁹, A Montgomery²⁸, T Pikuła²⁹) Portugal (H Barros³⁰, R Costa³⁰, C Rodrigues³⁰); Sweden (U Aden³¹); United Kingdom (ES Draper³², A Fenton³³, SJ Johnson³²); EFCNI (S Mader³⁴, N Thiele³⁴, JM Pfeil³⁴); Health Economics team (S Petrou³⁵, SW Kim³⁵, L Andronis³⁶); Inserm Coordination (J Zeitlin³⁷, AM Aubert³⁷, C Bonnet³⁷, R El Rafei³⁷, AV Seppänen³⁷). 1. Department of Family Medicine & Population Health (FAMPOP), Disability Studies, Faculty of Medicine & Health Sciences, University of Antwerp, Antwerp, Belgium. 2. Laboratory of Experimental Medicine and Pediatrics, Division of Neonatology, University of Antwerp, Antwerp, Belgium; Study Centre for Perinatal Epidemiology Flanders, Brussels, Belgium. 3. Centre for Developmental Disabilities, Neonatal Intensive Care, Oost Limburg Hospital, Genk, Belgium. 4. Vrije Universiteit Brussel Faculteit Geneeskunde en Farmacie; Paediatric Neurology, Universitair Ziekenhuis Brussels, Belgium. 5. Centre for Developmental Disabilities, Ghent University Hospital, Belgium. 6. Centre for Developmental Disabilities, Leuven University Hospital, Belgium; Department of Neuropediatrics, University of Leuven, Leuven, Belgium. 7 Centre for Developmental Disabilities, Antwerp University Hospital, Antwerp, Belgium; and Department of Primary & Interdisciplinary Care, Disability Studies, Faculty of Medicine, University of Antwerp, Antwerp, Belgium. 8. Department of Paediatrics, Hvidovre Hospital, Copenhagen University Hospital, Hvidovre, Denmark. 9. Department of Neonatology, Hvidovre Hospital, Hvidovre, Denmark. 10. Tallinn Children's Hospital, Tallinn; University of Tartu, Tartu, Estonia. 11. University of Tartu, Tartu University Hospital, Tartu, Estonia. 12. University of Tartu, Department of Paediatrics, Institute of Clinical Medicine, Tartu, Estonia. 13. Université Paris Cité, Inserm, INRAE, Centre for Research in Epidemiology and StatisticS (CRESS), Obstetrical Perinatal and Pediatric Epidemiology Research Team, EPOPé, F-75004 Paris, France. 14. Division of Pediatrics 2, Hôpital du Bocage; INSERM CIE1, CHRU Dijon, Université de Dijon, Dijon, France. 15. Université Paris Descartes and Assistance Publique Hôpitaux de Paris, Hôpitaux Universitaire Paris Centre Site Cochin, DHU Risks in pregnancy, Service de Médecine et Réanimation néonatales de Port-Royal, Paris, France. 16. Université Paris Cité, Inserm, INRAE, Centre for Research in Epidemiology and StatisticS (CRESS), Obstetrical Perinatal and Pediatric Epidemiology Research Team, EPOPé, F-75004 Paris, France; CHU Lille, Department of Neonatal Medicine, Jeanne de Flandre Hospital, 59000, Lille, France. 17. Department of Neonatology, Jeanne de Flandre Hospital, Lille CHRU Lille France. 18. Children's Hospital, University Hospital, Philipps University Marburg, Marburg, Germany Appendix; online. 19. University Medical Center, Homburg, Germany. 20. Institute of Quality Assurance Hesse, Eschborn, Germany. 21. Clinical Care and Management Innovation Research Area, Bambino Gesù Children's Hospital, IRCCS, Rome, Italy. 22. Maternal and Child Health Institute, Marche University and Salesi Hospital, Ancona, Italy. 23. Neonatal Intensive Care Unit, Ospedale degli Infermi, Rimini, Itay. 24. Neonatal Intensive Care Unit,
University Hospital S. Orsola-Malpighi, Bologna, Italy. 25. Department of Pediatrics and Neonatology, Modena University Hospital, Modena, Italy. 26. Department of Neonatology, Radboud University Medical Center, Nijmegen, the Netherlands. 27. Department of Neonatology, Wilhelmina Children's Hospital, Utrecht, The Netherlands. 28. Department of Neonatology, Poznan University of Medical Sciences, Poznan, Poland. 29. Department of Neonatology and Neonatal Infectious Diseases, Poznan University of Medical Sciences, Poznan, Poland. 30. EPIUnit-Institute of Public Health, University of Porto, Porto, Portugal. 31. Department of Women's and Children's Health, Karolinska Institutet, Stockholm, Sweden. 32. Department of Health Sciences, University of Leicester, Leicester, UK. 33. Newcastle University, Newcastle upon Tyne. 34. European Foundation for the Care of Newborn Infants (EFCNI), Munich, Germany. 35. Nuffield Department of Primary Care Health Sciences, University of Oxford, Oxford, UK. 36. Warwick Medical School, University of Warwick, Coventry, UK. 37. Université Paris Cité, Inserm, INRAE, Centre for Research in Epidemiology and StatisticS (CRESS), Obstetrical Perinatal and Pediatric Epidemiology Research Team, EPOPé, F-75004 Paris, France.

**Author contributions** SWK carried out the bulk of the analyses, drafted and revised the article, and approved the final article as submitted. SP designed the economic study, coordinated and supervised the analysis, critically reviewed the article and approved the final version for publication. LA drafted and revised the article critically and approved the final version for publication. AVS drafted and revised the article critically and approved the final version for publication. AMA drafted and revised the article critically and approved the final version for publication. MS drafted and revised the article critically and approved the final version for publication. JZ drafted and revised the article critically and approved the final version for publication. HB drafted and revised the article critically and approved the final version for publication. ESD drafted and revised the article critically and approved the final version for publication. All authors read the final manuscript and approved the final version for publication.

**Funding** This study was supported by a grant from the European Union's Horizon 2020 Research and Innovation Programme under grant agreement No 633724. Prof Petrou receives support as a National Institute for Health Research (NIHR) Senior Investigator (NF-SI-0616-10103) and from the NIHR Applied Research Collaboration Oxford and Thames Valley.

## Declarations

**Competing interests** The authors have no relevant financial or non-financial interests to disclose.

**Patient consent** Ethics approvals and parental consent were obtained according to national legislation in each country before data collection.

**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/.

## References

1. D'Onofrio, B. M., Class, Q. A., Rickert, M. E., Larsson, H., Långström, N., Lichtenstein, P. (2013). Preterm birth and mortality and morbidity: A population-based quasi-experimental study. *JAMA Psychiatry*, *70*(11), 1231–1240.

2. Arpino, C., Compagnone, E., Montanaro, M. L., Cacciatore, D., De Luca, A., Cerulli, A., Di Girolamo, S., Curatolo, P. (2010). Preterm birth and neurodevelopmental outcome: A review. *Child's nervous system*, *26*(9), 1139–1149.

3. CDC. *Preterm birth*. 2021 [cited 21 Oct 2021]; Available from: https://www.cdc.gov/reproductivehealth/maternalinfanthealth/pretermbirth.htm.

4. Moore, T., Hennessy, E. M., Myles, J., Johnson, S. J. Draper, E. S., Costeloe, K. L., Marlow, N. (2012) Neurological and developmental outcome in extremely preterm children born in England in 1995 and 2006: the EPICure studies. *BMJ* 345, e7961.

5. Larroque, B., Ancel, P-Y., Marret, S., Marchand, L., André, M., Arnaud, C., Pierrat, V., Rozé, J-C., Messer, J., Thiriez, G. (2008). Neurodevelopmental disabilities and special care of 5-year-old children born before 33 weeks of gestation (the EPIPAGE study): A longitudinal cohort study. *The Lancet*, *371*(9615), 813–820.

6. Serenius, F., Källén, K., Blennow, M., Ewald, U., Fellman, V., Holmström, G., Lindberg, E., Lundqvist, P., Maršál, K., Norman, M. (2013). Neurodevelopmental outcome in extremely preterm infants at 2.5 years after active perinatal care in Sweden. *JAMA*, *309*(17), 1810–1820.

7. Taylor, H. G., Klein, N., Anselmo, M. G., Minich, N., Espy, K. A., Hack, M. (2011). Learning problems in kindergarten students with extremely preterm birth. *Archives of pediatrics & adolescent medicine*, *165*(9), 819–825.

8. Khan, K. A., Petrou, S., Dritsaki, M., Johnson, S. J., Manktelow, B., Draper, E. S., Smith, L. K., Seaton, S. E., Marlow, N., Dorling, J. (2015). Economic costs associated with moderate and late preterm birth: a prospective population-based study. *BJOG: An International Journal of Obstetrics and Gynaecology*, *122*(11):1495-1505

9. Johnson, S., & Wolke, D. (2013). Behavioural outcomes and psychopathology during adolescence. *Early Human Development*, *89*(4), 199–207.

10. Saigal, S., & Doyle, L. W. (2008). An overview of mortality and sequelae of preterm birth from infancy to adulthood. *The Lancet*, *371*(9608), 261–269.

11. Petrou, S., Krabuanrat, N., & Khan, K. (2020). Preference-based health-related quality of life outcomes associated with preterm birth: A systematic review and meta-analysis. *PharmacoEconomics*, *38*(4), 357–373.

12. Zwicker, J. G., & Harris, S. R. (2008). Quality of life of formerly preterm and very low birth weight infants from preschool age to adulthood: A systematic review. *Pediatrics*, *121*(2), e366–e376.

13. van der Pal, S., Steinhof, M., Grevinga, M., Wolke, D., Verrips, G. (2020). Quality of life of adults born very preterm or very low birth weight: A systematic review. *Acta Paediatrica*, *109*(10), 1974–1988.

14. Olsen, J. E., Cheong, J. L. Y., Eeles, A. L., FitzGerald, T. L., Cameron, K. L. Albesher, R. A., Anderson, P. J., Doyle, L. W., Spittle, A. J. (2020). Early general movements are associated with developmental outcomes at 4.5–5 years. *Early Human Development*, *148*, 105115.

15. Dvir, Y., Frazier, J. A., Joseph, R. M., Mokrova, I., Moore, P. S. O'Shea, T. M., Hooper, S. R. Santos Jr, H. P. Kuban, K. (2019). Psychiatric symptoms: Prevalence, co-occurrence and
functioning among extremely low gestational age newborns at age ten years. *Journal of Developmental and Behavioral Pediatrics, 40*(9), 725.

16. Kelly, M. M. (2014). Assessment of life after prematurity in 9-to 10-year-old children. *MCN: The American Journal of Maternal/Child Nursing, 39*(1), 26–32.

17. Varni, J. W., Seid, M., & Kurtin, P. S. (2001). PedsQL™4.0: reliability and validity of the pediatric quality of life inventory™version 4.0 generic core scales in healthy and patient populations. *Medical Care, 39*, 800–812.

18. Thébaud, B., Goss, K. N., Laughon, M., Whitsett, J. A., Abman, S. H., Steinhorn, R. H., Aschner, J. L., Davis, P. G., McGrath-Morrow, S. A., Soll, R. F. (2019). Bronchopulmonary dysplasia. *Nature Reviews Disease Primers, 5*(1), 1–23.

19. Gou, X., Yang, L., Pan, L., Xiao, D. (2018). Association between bronchopulmonary dysplasia and cerebral palsy in children: A meta-analysis. *British Medical Journal Open, 8*(9), e020735.

20. Darlow, B. A., Hutchinson, J. L., Henderson-Smart, D. J., Donoghue, D. A., Simpson, J. M., Evans, N. J. (2005). Prenatal risk factors for severe retinopathy of prematurity among very preterm infants of the Australian and New Zealand Neonatal Network. *Pediatrics, 115*(4), 990–996.

21. Rogers, B., Msall, M., Owens, T., Guernsey, K., Brody, A., Buck, G., Hudak, M. (1994). Cystic periventricular leukomalacia and type of cerebral palsy in preterm infants. *The Journal of pediatrics, 125*(1), S1–S8.

22. Patole, S. (2007). Prevention and treatment of necrotising enterocolitis in preterm neonates. *Early Human Development, 83*(10), 635–642.

23. Theunissen, N. C. M., Veen, S., Fekkes, M., Koopman, H. M., Zwinderman, K. A. H., Brugman, E., Wit, J-M. (2001). Quality of life in preschool children born preterm. *Developmental Medicine and Child Neurology, 43*(7), 460–465.

24. Zeitlin, J., Maier, R. F., Cuttini, M., Aden, U., Boerch, K., Gadzinowski, J., Jarreau, P-H., Lebeer, J., Norman, M., Pedersen, P., Petrou, S., Pfeil, J. M., Toome, L., van Heijst, A., Van Reempts, P., Varendi, H., Barros, H., Draper, E. S., EPICE and SHIPS Research Group. (2020). Cohort profile: Effective Perinatal Intensive Care in Europe (EPICE) very preterm birth cohort. *International Journal of Epidemiology, 49*(2), 372–386.

25. PedsQL.org. PedsQL™. 2021 [cited 15 Aug 2021]; Available from: http://pedsql.org/index.html.

26. Khan, K. A., Petrou, S., Dritsaki, M., Johnson, S. J., Manktelow, B., Draper, E. S., Smith, L. K., Seaton, S. E., Marlow, N., Dorling, J., Field, D. J., Boyle, E. M. (2015). Economic costs associated with moderate and late preterm birth: a prospective population-based study. *BJOG: An International Journal of Obstetrics and Gynaecology, 122*(11), 1495–1505.

27. Stylianou-Riga, P., Kouis, P., Kinni, P., Rigas, A., Papadouri, T., Yiallouros, P. K., Theodorou, M. (2018). Maternal socioeconomic factors and the risk of premature birth and low birth weight in Cyprus: A case-control study. *Reproductive Health, 15*(1), 157.

28. Kwinta, P., & Pietrzyk, J. J. (2010). Preterm birth and respiratory disease in later life. *Expert Review of Respiratory Medicine, 4*(5), 593–604.

29. Kim, S.W., Andronis L., Seppänen, A-V., Aubert, A. M., Zeitlin, J., Barros, H., Draper, E. S., Petrou, S., on behalf of the SHIPS Research Group. (2021). Economic costs at age five associated with very preterm birth: Multinational European cohort study. *Pediatric Research 122*(11), 1495–1505.

30. Zeitlin, J., Szamotulska, K., Drewniak, N., Mohangoo, A. D., Chalmers, J., Sakkeus, L., Irgens, L., Gatt, M., Gissler, M., Blondel, B. (2013). Preterm birth time trends in Europe: a study of 19 countries. *BJOG: An International Journal of Obstetrics & Gynaecology, 120*(11), 1356–1365.

31. Zeitlin, J., Bonamy, A-K. E., Piedvache, A., Cuttini, M., Barros, H., Van Reempts, P., Mazela, J., Jarreau, P-H., Gortner, L., Draper, E. S., Maier, R. F., the EPICE Research Group. (2017). Variation in term birthweight across European countries affects the prevalence of small for gestational age among very preterm infants. *Acta Paediatrica, 106*(9), 1447–1455.

32. O'Reilly, M., Sozo, F., Harding, R. (2013). Impact of preterm birth and bronchopulmonary dysplasia on the developing lung: Long-term consequences for respiratory health. *Clinical and Experimental Pharmacology and Physiology, 40*(11), 765–773.

33. MacKinnon, D. P., Fairchild, A. J., & Fritz, M. S. (2007). Media-tion analysis. *J Annual Review of Psychology, 58*, 593–614.

34. Baron, R. M., & Kenny, D. A. (1986). The moderator–mediator variable distinction in social psychological research: Conceptual, strategic, and statistical considerations. *J Journal of personality social psychology, 51*(6), 1173.

35. VanderWeele, T. (2015). *Explanation in causal inference: methods for mediation and interaction*. Oxford University Press.

36. StataCorp. Stata Structural Equation Modeling Reference Manual release 17. 2021 [cited 15 Aug 2021]; Available from: https://www.stata.com/manuals/sem.pdf.

37. Rucker, D. D., Preacher, K. J., Tormala, Z. L., Petty, R. E. (2011). Mediation analysis in social psychology: Current practices and new recommendations. *Social and Personality Psychology Compass, 5*(6), 359–371.

38. Lockwood, C.M., & MacKinnon, D.P. (1998). Bootstrapping the standard error of the mediated effect. In: *Proceedings of the 23rd annual meeting of SAS Users Group International*. Citeeseer.

39. Sriram, S., Schreiber, M. D., Msall, M. E., Kuban, K. C. K., Joseph, R. M., O' Shea, T. M., Allred, E. N., Leviton, A., for the ELGAN Study Investigators. (2018). Cognitive development and quality of life associated with BPD in 10-year-olds born preterm. *Pediatrics 141*(6), e20172719.

40. Vederhus, B. J., Markestad, T., Eide, G. E., Graue M., Halvorsen, T. (2010). Health related quality of life after extremely preterm birth: A matched controlled cohort study. *Health and Quality of Life Outcomes, 8*(1), 1–8.

41. Brady, J. M., Zhang, H., Kirpalani, H., DeMauro, S. B. (2019). Living with severe bronchopulmonary dysplasia—Parental views of their child's quality of life. *The Journal of Pediatrics, 207*, 117–122.

42. Joffe, A. R., Wong, K., Bond, G. Y., Moez, E. K., Acton, B. V., Dinu, I. A., Yap, J. Y. K., Robertson, C. M. T., for the Western Canadian Complex Pediatric Therapies Follow-up Program. (2020). Kindergarten-age neurocognitive, functional, and quality-of-life outcomes after liver transplantation at under 6 years of age. *Pediatric Transplantation, 24*(2), e13624.

43. Varni, J. W., Burwinkle, T. M., Seid, M., Skarr, D. (2003). The PedsQL™* 4.0 as a pediatric population health measure: feasibility, reliability, and validity. *Ambulatory Pediatrics, 3*(6), 329–341.

44. Hilliard, M. E., Lawrence, J. M., Modi, A. C., Anderson, A., Crume, T., Dolan, L.M., Merchant, A. T., Yi-Frazier, J. P., Hood, K. K., for the SEARCH for Diabetes in Youth Study Group (2013). Identification of minimal clinically important difference scores of the PedsQL in children, adolescents, and young adults with type 1 and type 2 diabetes. *Diabetes Care, 36*(7), 1891–1897.

45. Gire, C., Resseguier, N., Brévaut-Malaty, V., Marret, S., Cambonie, G., Souksi-Medioni, I., Müller, J-B., Garcia, P., Berbis, J., Tosello, B., Pascal Auquier on behalf of the GPQoL study Group. (2019). Quality of life of extremely preterm school-age children without major handicap: A cross-sectional observational study. *Archives of Disease in Childhood, 104*(4), 333–339.

46. Ni, Y., Johnson, S., Marlow, N., Wolke, D. (2021). Reduced health-related quality of life in children born extremely preterm in 2006 compared with 1995: the EPICure Studies. Archives of Disease in Childhood-Fetal and Neonatal Edition
47. Rautava, L., Häkkinen, U., Korvenranta, E., Andersson, S., Gissler, M., Hallman, M., Korvenranta, H., Leipälä, J., Linna, M., Peltola, M., Tammela, O., Lehtonen, L. (2009). Health-related quality of life in 5-year-old very low birth weight infants. *The Journal of Pediatrics*, **155**(3), 338-343.e3.

48. Berbis, J., Einaudi, M. A., Simeoni, M. C., Brévaut-Malaty, V., Auquier, P., d' Ercole, C., Gire, C. (2012). Quality of life of early school-age French children born preterm: A cohort study. *European Journal of Obstetrics & Gynecology and Reproductive Biology*, **162**(1), 38–44.

49. Chien, L. Y., Chou, Y-H., Ko, Y-L., Lee, C-F. (2006). Health-related quality of life among 3–4-year-old children born with very low birthweight. *Journal of advanced Nursing*, **56**(1), 9–16.

50. Krieger, N., Huynh, M., Li, W., Waterman, P. D., Van Wye, G. (2018). Severe sociopolitical stressors and preterm births in New York City: 1 September 2015 to 31 August 2017. *Journal of Epidemiology and Community Health*, **72**(12), 1147–1152.

51. Khanolkar, A. R., Wedrén, S., Essén, B., Sparén, P., Koupil, I. (2015). Preterm and postterm birth in immigrant-and Swedish-born parents: A population register-based study. *European Journal of Epidemiology*, **30**(5), 435–447.

52. Vieira, M. E. B., & Linhares, M. B. M. (2016). Quality of life of individuals born preterm: A systematic review of assessment approaches. *Quality of Life Research*, **25**(9), 2123–2139.

53. Beaudoin, S., Tremblay, G. M., Croitoru, D., Benedetti, A., Landry, J. S. (2013). Healthcare utilization and health-related quality of life of adult survivors of preterm birth complicated by bronchopulmonary dysplasia. *J Acta Paediatrica*, **102**(6), 607–612.

54. Linsell, L., Malouf, R., Morris, J., Kurinczuk, J. J., Marlow, N. (2016). Prognostic factors for cerebral palsy and motor impairment in children born very preterm or very low birthweight: A systematic review. *Developmental Medicine & ChildNeurology*, **58**(6):554-569

55. Kenny, D. (2008). Reflections on mediation. *J Organizational Research Methods*, **11**(2), 353–358.

56. Kenny, D. *Mediation*. 2021 [cited 5 Nov 2021]; Available from: http://davidakenny.net/cm/mediate.htm.

57. Palta, M., & Sadek-Badawi, M. (2008). PedsQL relates to function and behavior in very low and normal birth weight 2- and 3-year-olds from a regional cohort. *Quality of Life Research*, **17**(5), 691–700.

58. Yu, Y., Liew, Z., Wang, A., Arah, O. A., Li, J., Olsen, J., Cnattin-gius, S., Qin, G., Obel, C., Fu, B. (2019). Mediating roles of pre-term birth and restricted fetal growth in the relationship between maternal education and infant mortality: A Danish population-based cohort study. *PLoS medicine*, **16**(6), e1002831.

59. Jozefiak, T., Larsson, B., Wichström, L., Mattejat, F., Ravens-Sieberer, U. (2008). Quality of Life as reported by school children and their parents: A cross-sectional survey. *Health Quality of Life Outcomes*, **6**(1), 1–11.

60. Theunissen, N. C., Vogels, T. G., Koopman, H. M., Verrips, G. H., Zwinderman, K. A., Verloove-Vanhorick, S. P., Wit, J. M. (1998). The proxy problem: Child report versus parent report in health-related quality of life research. *Quality of Life Research*, **7**(5), 387–397.

61. Conijn, J. M., Smits, N., & Hartman, E. E. (2020). Determining at what age children provide sound self-reports: An illustration of the validity-index approach. *Assessment*, **27**(7), 1604–1618.

**Publisher's Note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>