Impact of moderate and late preterm birth on neurodevelopment, brain development and respiratory health at school age: protocol for a longitudinal cohort study (LaPrem study).

INTRODUCTION: Children born moderate to late preterm (MLP, 32-36 weeks' gestation) account for approximately 85% of all preterm births globally. Compared with children born at term, children born MLP are at increased risk of poor neurodevelopmental outcomes. Despite making up the largest group of preterm children, developmental outcomes of children born MLP are less well studied than in other preterm groups. This study aimed to (1) compare neurodevelopmental, respiratory health and brain magnetic resonance imaging (MRI) outcomes between children born MLP and term at 9 years of age; (2) examine the differences in brain growth trajectory from infancy to 9 years between children born MLP and term; and in children born MLP; (3) examine the relationship between brain development and neurodevelopment at 9 years; and (4) identify risk factors for poorer outcomes at 9 years. METHODS AND ANALYSIS: The "LaPrem" (Late Preterm MRI Study) study is a longitudinal cohort study of children born MLP and term controls, born at the Royal Women's Hospital in Melbourne, Australia, between 2010 and 2013. Participants were recruited in the neonatal period and were previously followed up at 2 and 5 years. This 9-year school-age follow-up includes neuropsychology, motor and physical activities, and lung function assessments, as well as brain MRI. Outcomes at 9 years will be compared between birth groups using linear and logistic regressions. Trajectories of brain development will be compared between birth groups using mixed effects models. The relationships between MRI and neurodevelopmental outcomes, as well as other early predictors of poor 9-year outcomes, will be explored using linear and logistic regression. ETHICS AND DISSEMINATION: This study was approved by the human research ethics committee at the Royal Children's Hospital, Melbourne, Australia. Study outcomes will be disseminated through peer-reviewed publications, conference presentations and social media. © Author(s) (or their employer(s)) 2021. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.

This is the first cohort study of Australian children born moderate to late preterm to be assessed at school age with longitudinal MRI brain data.

Extensive data collected at term equivalent, 2, 5 and 9 years of age will enable us to compare a wide range of important outcomes between birth groups and to better understand developmental trajectories in this population.

We expect some attrition at the 9-year follow-up and recognise that differences between children who do and do not complete an assessment may introduce bias.

As this is a hospital-based cohort, the results may not apply to the general population.

We are unable to assess all important health outcomes due to resource and time constraints.

Introduction

Preterm birth (<37 weeks’ gestation) accounts for an estimated 10% of all live birth globally, of which almost 85% are born moderate to late preterm (MLP: 32–36 weeks’ gestation).1 In an Australian context, 22 000 children were born MLP in 2017, comprising approximately 7% of all live births.2 Despite making up the largest group of preterm births, developmental outcomes of children born MLP are less well understood than preterm groups born earlier (ie, <32 weeks). Emerging research from our team and others over the last decade has highlighted that MLP birth is not as benign as previously thought.3–9 Compared with children born at term, MLP children have (1) smaller and less mature brains at term equivalent age using brain MRI5; (2) increased rates of developmental delay and poorer social–emotional development at 2 years of age3; (3) poorer school readiness, worse cognitive performance and greater special education needs6 7; and (4) increased respiratory morbidity in infancy and early childhood.8 9

Such findings have increased awareness of the morbidities associated with MLP birth and are influencing obstetric decision-making about the timing of MLP births. However, permanent and meaningful effects on obstetric decision-making await reports of clinically important school-age outcomes in an Australian context, which we will obtain with this study. Empowering obstetricians with the knowledge of whether adverse outcomes in the MLP group persist into school age will enable better decision-making regarding the risks and benefits of delivering a baby ‘just a bit early’. Knowledge of developmental trajectories for children born MLP, as well as risk factors for poor outcomes at this age, will provide the evidence base for long-term MLP developmental and health surveillance guidelines. Furthermore, an improved understanding of school-age outcomes for children born MLP will have the potential to inform tailored interventions to suit children’s specific needs, information that is not available from assessments in the first few years of development.

Early developmental assessments are only moderately predictive of school-age cognitive functioning and academic achievement, with more accurate information obtained through school-age assessments.10 Likewise, lung function is better assessed at school age. Currently, it is recommended that very preterm (<32 weeks’ gestation) children have regular developmental and health follow-up throughout childhood.11 However, similar guidelines are lacking for MLP children, presumably because the few studies describing outcomes during critical periods of development, such as at school age, have not characterised which developmental domains may be affected compared with children born at term. Without this vital information, we cannot design and implement effective interventions to improve outcomes for this population. Moreover, little is known of the risk factors to alert parents and health practitioners to MLP children at greatest risk and who might benefit most from surveillance.

This 9-year follow-up of a pre-existing MLP cohort and term-born controls will allow us to identify specific domains affected by MLP birth, as well as enable us to identify risk factors for poor outcomes, which in turn will inform the design and implementation of effective early interventions to improve outcomes for these vulnerable children.

Aims and hypotheses

Aim 1

The first aim was to compare outcomes at 9 years between children born MLP and at term, including (1) neurodevelopment outcomes (neuropsychological, academic, behavioural, physical activity and motor function); (2) brain structure and function, using advanced brain MRI measures; and (3) respiratory health and lung function.

Hypothesis 1

At 9 years, compared with term controls, children born MLP will have (1) lower IQ, poorer attention, executive function, language, memory and academic achievement, more problems with behaviour, impaired motor function and lower physical activity levels; (2) smaller brain volumes, aberrant cortical morphology and less mature white matter microstructure; and (3) higher rates of wheezing illnesses and more airflow obstruction on lung function testing.

Aim 2

The second aim was to examine the differences in brain growth trajectory, from infancy to 9 years, between children born MLP and term using longitudinal MRI measures.

Hypothesis 2

Children born MLP will have slower trajectory of brain growth from infancy to 9 years compared with children born at term.

Aim 3

The third aim was to examine the relationships between neurodevelopment at 9 years and (1) brain structure at 9 years and (2) the changes in brain measures between term-equivalent age and 9 years in MLP children.

Hypothesis 3

Aberrant brain development at 9 years, and slower growth and maturation in brain volumes, cortical and white matter brain measures between term-equivalent age and 9 years will be related to poorer developmental outcomes at 9 years for children born MLP.

Aim 4

The fourth aim was to identify risk factors in pregnancy and early infancy for poorer school age outcomes

(neurodevelopment, brain structure and function and respiratory health) for children born MLP.

Hypothesis 4

Foetal growth restriction, newborn respiratory problems, smaller brain volumes and less mature white matter microstructure at term-equivalent age, suboptimal neurobehaviour at term equivalent age, and developmental delay and respiratory ill health at 2 years will be associated with poorer school age outcomes. Environmental risk factors (eg, passive smoke exposure and adverse social risk) during pregnancy and early infancy will also be associated with poorer outcomes.

Methods and analysis

Protocol details

Protocol details: V.4. Dated 6 April 2020. Updated to include physically distanced assessment on 26 August 2020.

Design

This study is a prospective longitudinal follow-up of an established cohort.

Study population

The LaPrem study is a cohort study of children born MLP (n=201) and a control group of children born at term (≥37 weeks’ gestation and >2499 g birth weight, n=201). Children were recruited in the neonatal period from the Royal Women’s Hospital, a large tertiary hospital in Melbourne, Australia, between 2010 and 2013. Children with congenital abnormalities or genetic syndromes known to affect development were excluded. In addition, term-born infants were excluded if they were unwell at birth, received resuscitation, were admitted to the neonatal nursery or were identified as having conditions affecting growth or development.

The LaPrem cohort has previously been assessed at term equivalent age and at 2 and 5 years. Results of these assessments have been published previously.3–5 12–18 This protocol paper details the 9-year follow-up of this cohort. All surviving children still enrolled in the LaPrem cohort will be eligible to participate in this study.

Assessment procedure

Study visits

Families will be contacted by a research nurse and invited to participate in the 9-year assessment. Information about the assessment will be provided through a plain language information and consent form, as well as through a handout designed specifically for the participating children. If parents agree for their child to participate in the 9-year follow-up, they will be asked to sign the participant information and consent form. This consent form has separate, optional consent options for (1) an MRI scan (which is preceded by a session in a mock MRI); (2) a lung function assessment; (3) permission to access the child’s National Assessment Programme—Literacy and Numeracy (NAPLAN) results; and (4) permission to access the child’s Medicare record, which includes data on outpatient (hospital and health practitioners) and inpatient visits, as well as prescription medications. Parents may consent to some or all of these study aspects on behalf of their child. More information about each of these components is detailed under the specific outcome measures.

Assessments will take place over 2 days, with regular breaks as required to minimise fatigue. All assessors will be blinded to birth group and results of previous assessments. Assessments commenced in June 2019, with data collection proposed to finish in December 2022.

Assessment procedure: COVID-19 adaptations

The COVID-19 pandemic, which escalated in March 2020 in Victoria, Australia, has necessitated changes to the assessment procedure for this study to ensure, first, that relevant institutional and government guidelines are adhered to and, second, that study participants and assessors are provided with a safe environment. While community transmission of COVID-19 in Victoria, Australia, remains a concern, participants will be offered one of two assessment formats:

A face-to-face assessment, with physical distancing measures and hygiene procedures in place to minimise risk of viral transmission.

In the case of restrictions prohibiting face-to-face assessments for research purposes, a remote assessment will be offered, in which parts of the cognitive and motor components will be conducted via the video conferencing platform, with the participants in their own home. When these restrictions are lifted, these participants may be offered the opportunity to complete the remaining components of the assessments, including brain MRI and lung function, if possible.

The option of remote assessment may also be offered outside the context of such restrictions to reduce the total time spent onsite or to participants who do not wish to attend the hospital at all due to concerns about viral transmission, or to regional, interstate or international participants who are unable to travel. Further details about adapted assessment procedures are provided in online supplemental material 1.

Outcome measures

All assessments will be conducted by trained professionals blinded to group allocation and clinical history.

General cognitive ability will be assessed using the Wechsler Intelligence Scale for Children—Fifth Edition, Australia and New Zealand.19 Ten core subtests will be administered, which will provide an estimate of full-scale IQ as well as composite scores on five domains of cognitive functioning: verbal comprehension, visual spatial, fluid reasoning, working memory and processing speed.

Attention will be assessed using subtests from the Test of Everyday Attention for Children, Second Edition.20 The five selected subtests include (1) Vigil, a measure of sustained attention where participants keep count of the number of sounds they hear on a recording; (2) Hector, a measure of focused visual attention where participants search for target objects in busy visual arrays; (3) Troy dual task, a more complex divided attention task that combines the previous two tasks, so that the participant searches for visual targets while counting sounds at the same time; (4) simple reaction time, a measure of reaction time where participants press a computer key in response to an on-screen stimulus; and (5) Reds, Blues, Bags, and Shoes, a measure of the participants’ ability to rapidly and fluently switch attention between different task requirements.

Executive functions will be examined using two tests: (1) the Contingency Naming Test, which measures aspects of executive function such as response inhibition, working memory and mental flexibility,21 and involves four trials that increase in task complexity, where the participant names shapes and colours according to learnt rules; (2) the Rey Complex Figure Test, which involves reproducing a complicated line drawing, which is assessed for accuracy and strategy.22

Language skills will be assessed using one subtest from the NEPSY-II,23 the Comprehension of Instructions Subtest, which assesses ability to process and execute oral instructions of increasing syntactic complexity. For each item, the child points to stimuli in response to verbal instructions.

Memory will be assessed using the California Verbal Learning Test—Children’s Version, testing participants’ ability to learn a verbally presented word list over multiple learning trials.24 The participant’s ability to remember the list of words over both short and long time delays will be assessed.

Academic achievement will be assessed using the Wechsler Individual Achievement Test—Third Edition, Australia and New Zealand.25 The selected subtests will assess reading (single word reading and pseudoword decoding), spelling and math computation. NAPLAN results will provide further data on academic achievement.

Motor function will be assessed using the Movement Assessment Battery for Children, Second Edition,26 a gold standard for assessing motor impairment in children 3–16 years of age. The assessment is divided into three subscales: manual dexterity, aiming and catching and balance.

Cerebral palsy will be diagnosed based on a neurological examination including loss of motor function, abnormal tone and tendon reflexes, or recorded if the child has previously received a diagnosis. The five-level Gross Motor Function Classification System (GMFCS) will be used to further classify motor function for children with cerebral palsy. The GMFCS is a reliable and widely used scale that classifies children’s gross motor function based on self-initiated movement.27

Physical activity levels will be assessed with a high-resolution (13-bit), water resistant to 1.5 m ankle mounted triaxial accelerometer (Axivity AX3, Newcastle on Tyne, UK) worn for seven consecutive days (including a full weekend) for 24 hours/day. The key outcome measure is mean steps per day, with additional measures including time spent inactive and the time spent at different exercise intensities per day. Ankle-mounted accelerometers have been used previously to assess activity in children with motor impairments,28 and the ankle is well established in research as the optimal position for stepping movements across a range of devices.29 The monitor will be given to children at the follow-up assessment and returned by post.

Physical activity participation will be measured using the Physical Activity Questionnaire for Older Children (PAQ-C).30 This questionnaire will be completed by the child during the assessment, with assistance from the assessor, if necessary, using an online survey linked to a Research Electronic Data Capture (REDCap) database.31 The PAQ-C is a self-report 7-day recall questionnaire which was developed for children aged 8–14 years and includes a checklist of common sporting and leisure activities, as well as questions about physical activity levels during the school day, after school and on weekends.30 The sporting and leisure activity list was updated for the current study to include common activities undertaken by Australian children, with permission from the authors of the outcome measure. The reliability and validity of the PAQ-C have been investigated in several different geographical regions, and studies have demonstrated good reliability and fair–moderate validity with accelerometers.32 In addition to the PAQ-C, two brief questions assessing participation in sedentary screen behaviours will be administered. The questions are adapted from the Centers for Disease Control and Prevention 2017 Middle School Youth Risk Behaviour Survey and are as follows:

‘On an average school day, how many hours do you watch TV?’

On an average school day, how many hours do you play video or computer games or use a computer for something that is not school work? (Count time spent on things such as Xbox, PlayStation, an iPad or other tablet, a smartphone, texting, YouTube, Instagram, Facebook or other social media.)33

Grip strength will be measured using a hand-held dynamometer. Strength (force in kilogram) will be measured with the child in a seated position, with their elbow flexed at 90°, shoulder adducted into their body and their forearm parallel to the floor. Each child will be given a practice trial, followed by three subsequent trials for each of (1) the child’s preferred hand; (2) the child’s non-preferred hand; and (3) bimanual, from which the maximum for each condition will be taken. Grip strength is an inexpensive and easily implemented assessment which, in adulthood, can be an indicator of health as poorer grip strength has been associated with functional impairment,34 cardiovascular disease and all-cause mortality.35

Anthropometric measures include height, weight, waist circumference, chest circumference, head circumference, mid-arm circumference and lower limb length, and will be measured according to standard guidelines. In addition, blood pressure, mean arterial pressure and heart rate will be taken at rest using an automated device.

Brain MRI

A trained research assistant will lead each child through a ‘mock’ MRI to familiarise the child with the process before undertaking the MRI scan. Only children who are comfortable and able to stay still during the mock MRI procedure will undergo MRI, ensuring high compliance and quality images. Imaging will be undertaken using a 3-Tesla Siemens MAGNETOM Prisma scanner at the Royal Children’s Hospital MRI Centre using sequences listed in table 1. Children will be scanned within 2 weeks of their neuropsychological assessment. We will quantify brain tissue volumes for over 100 brain regions and structures, as well as cortical morphology, including surface area, thickness, curvature, folding and sulcal depth for around 68 cortical regions, using the validated software FreeSurfer V.7.36 A mathematical ratio of the T1 to T2 images will be calculated to derive myelin maps,37 providing a sensitive measure of myelination throughout the entire brain. Diffusion modelling techniques will be used to measure microstructural complexity, including traditional diffusion tensor imaging,38 as well as more advanced models: (1) the novel Spherical Mean Technique, which can estimate diffusivity and neurite density without influence from crossing fibres39; and (2) fixel-based analyses, which can quantify the apparent fibre density of individual fibres within each diffusion image voxel.40 Seventy-two major white matter tracts will be virtually dissected using TractSeg, a novel convolutional neural network-based approach that directly segments tracts in the field of fibre orientation distribution function peaks with high accuracy, without using tractography, image registration or parcellation.41 We will use modern mathematical methods of graph theory and network-based statistics42 to determine how different parts of the brain interconnect both structurally (using constrained spherical deconvolution43) and functionally (using resting state functional MRI). The sensitive MRI analysis techniques will be useful for determining brain developmental maturity and hold great promise for uncovering the subtler brain abnormalities associated with MLP birth.

Table 1

MRI sequences

1	T1- weighted multiecho magnetisation prepared rapid gradient echo (MP-RAGE) images with echo planar image-navigated prospective motion compensation (TR, 2550 ms; TE, TE1 2.14 ms, TE2 3.94 ms, TE3 5.77 ms, TE4 7.5 ms; IV, 0.9 mm3; TA, 06:47)
2	3D T2-weighted turbo spin echo images (TR, 3200 ms; TE, 408 ms; IV, 0.9 mm3; TA, 04:43)
3	Multishell simultaneous multislice echo planar diffusion images (TR, 3500 ms; TE, 67 ms; IV, 2.0 mm3; b=750 s/mm2, 25 gradient directions; b=2000 s/mm2, 45 directions; and b=2800 s/mm2, 60 directions; 19 b=0 s/mm2; multiband acceleration factor 2; TA, 10:59) with matching reverse-phase encoding images (TA, 01:40)
4	Quantitative susceptibility mapping gradient echo images (TR, 30 ms; TE1 6.32 ms, TE2 11.03 ms, TE3 15.74 ms, TE4 20.45 ms, TE5 25.16 ms; IV, 1.0 mm3; TA, 05:38), acquired on a subset starting November 2019
5	Multiecho gradient recalled echo planar resting state functional MRI images with prospective acquisition correction (TR, 1500 ms; TE, 33 ms; IV, 2.5 mm3; multiband acceleration factor 3; total TA, 13:04, in two consecutive acquisitions of TA, 06:32, for compliance) with matching reverse-phase encoding images (TA, 0:56)

IV, isotropic voxel; TA, time of acquisition; TE, echo time; TR, repetition time.

Lung function

Lung function (detailed in table 2) will be measured in the Respiratory Function Laboratory. Spirometry will be measured according to standard guidelines of the American Thoracic Society44 and expressed primarily as z-scores.45 Plethysmography to assess lung volumes and measurement of gas transfer are standardised tests routinely used in clinical practice. The Lung Clearance Index is an innovative test that measures ventilation inhomogeneity; it is more sensitive than spirometry in detecting early airway dysfunction, and it has reference data available.46

Table 2

Lung function variables

Lung function	Variable measured
Airflow	FEV1, FEF25%–75%
Lung volumes	FVC, TLC, RV
Reversibility	Bronchodilator to determine how much airway obstruction is reversible.
Gas exchange	DLco: measures alveolar–capillary membrane pathology
Ventilation inhomogeneity	Lung Clearance Index, second (airway conductance)

DLco, diffusing capacity of the lung for carbon monoxide; FEF25%–75%, forced mid-expiratory flow; FEV1, forced expiratory volume in one second; FVC, forced vital capacity; RV, residual volume; TLC, total lung capacity.

Parent questionnaires

Parents/caregivers will be asked to complete several questionnaires (table 3) via online surveys which are linked to a REDCap database. Questionnaires, which are to be completed by the primary caregiver, include four measures of child behaviour, parent-reported assessments of child motor skills and respiratory symptoms, as well as demographic data to determine social risk. In addition, both parents will be asked to complete the Hospital Anxiety and Depression Scale.47 Parents who report a moderate–severe level of symptoms (score above 11) will be contacted, and relevant support services will be suggested to them. If they are unable to be contacted by phone within a reasonable time frame, a handout of contact information for available local supports will be provided by email.

Table 3

Parent questionnaires

Questionnaire	Description
Developmental Coordination Disorder Questionnaire49	Parent-reported measure designed to identify subtle motor problems in children; construct validity has been established, and cut-off scores to determine risk of motor impairment are available for this age group.
Strengths and Difficulties Questionnaire50	Parent-reported, well-validated measure that assesses children’s overall behaviour problems, emotional symptoms, hyperactivity/attention, peer relationship problems and prosocial behaviour
Behaviour Rating Inventory of Executive Functioning, Second Edition51	Parent-reported measure designed to assess children’s executive function behaviours in the home environment
Social Communication Questionnaire—Current Version52	Parent-reported measure used to evaluate children’s communication skills and social functioning
ADHD Rating Scale-5 for Children and Adolescents53	Parent-reported measure to assess potential symptoms of ADHD. The ADHD Rating Scale-5 is a valid and reliable instrument that is widely used by practitioners in screening, diagnosis and treatment evaluation.
The International Study of Asthma and Allergies in Childhood Core Questionnaire54	Parent-reported measure to assess for symptoms of asthma, allergic rhinitis and eczema
Social Risk Index55	Assesses six aspects of social risk, including family structure, education of primary caregiver, occupation of primary income earner, employment status of primary income earner, language spoken at home and maternal age at birth. This measure has been used successfully in previous studies considering children born preterm.56
The Hospital Anxiety and Depression Scale 47	To be completed by both caregivers, this questionnaire is a self-rating scale that measures anxiety and depression in both hospital and community settings

ADHD, attention deficit hyperactivity disorder.

Assessment feedback for participants

All participants will receive written results from all assessments in a formal report. Results will be articulated in terms of performance ranges (ie, below average, average and above average) for neuropsychology and motor assessments. The report will include contact details for the study team, and families will be encouraged to contact the team if they wish to discuss the information in the report. If there are assessment results that require clinical follow-up, a researcher will contact the family directly to discuss the findings and options for referral to appropriate services as indicated.

Sample size

The sample size is determined by the size of the existing cohort (n=402). Based on our past experience of assessing many similar age cohorts, we conservatively estimate an overall follow-up rate of 90% at the 9-year assessment; this will result in 180 MLP and 180 controls. There are 27 pairs of twins in the MLP group (27% of participants) and none in the term control group. Assuming an intraclass correlation coefficient of 0.6 between outcomes for twins and that 27% of those who attend the follow-up are twins (n=48), we will have an effective sample size of 155 MLP (based on a design effect of 1.16). Samples of 155 MLP and 180 controls will enable us to detect differences between the groups as small as 0.31 SD with 80% power (based on a two-sided test with α=0.05). For most outcomes that are measured continuously, a difference of 0.31 SD between groups would be clinically important and is a difference that we could potentially observe. For proportions, if the event rate is 50% in the controls, we will have 80% power to detect differences of ±16%. As the event rate moves away from 50%, we will be able to detect smaller absolute differences between the two groups with 80% power, which are again clinically important differences, and differences that we could potentially observe.

Data management

All neuropsychology, motor and lung function outcome data will be entered into an electronic REDCap database by a member of the research team. Parent-completed questionnaires, as well as the PAQ-C, which is completed by the child, will be completed as online surveys linked to the REDCap database. Data will be stored in a secure data storage facility based at the Murdoch Children’s Research Institute (MCRI) and backed up daily onto a central server. Data will be stored until the youngest participant reaches 25 years of age. Participant’s deidentified MRI data will be stored and backed up in the secure developmental imaging facility initially on-site at MCRI and later archived at a secure off-site facility. The Quantitative susceptibility mapping (QSM) MRI images will be analysed in collaboration with the Florey Institute of Neuroscience and Mental Health, which has previously applied these techniques in neurodegeneration research.48 A mutual transfer agreement will be set up prior to commencement of analysis. Any modifications to the approved study protocol will need a formal amendment to the research and ethics committee of the Royal Children’s Hospital.

Statistical analysis

Outcomes will be compared between MLP and term-born controls using linear (continuous) and logistic (binary) regression, both unadjusted and adjusted for potential confounding variables, including gestation at birth, sex and social risk. Models will be fitted using generalised estimating equations (GEEs) and reported with robust (sandwich) estimates of SEs to allow for clustering of siblings within a family.

Longitudinal differences in brain MRI measures from term equivalent age to 9 years will be compared between MLP and controls using mixed effects models, including a fixed effect of age, group and the interaction between age and group, and a random effect to allow for the repeated observations within an individual.

To determine the relationships between imaging metrics and neurodevelopmental outcomes, we will use linear and logistic regressions (fitted with GEEs), with adjustment for social and environmental factors, including gestation at birth, sex and social risk, to determine the relationships over and above these factors.

Linear and logistic regressions (fitted using GEEs) will be used to assess predictors (brain MRI, including quantitative and qualitative measures at term equivalent age, gestation at birth, sex, birth weight z-score and 2 year neurodevelopment) of outcomes at 9 years. Initially, potential predictors will be explored using univariable regression for each predictor–outcome combination, before combining predictors into a multivariable model for each outcome to assess independent predictors. Social factors (maternal education and occupation of primary caregiver) will be assessed as modifiers in these relationships by the inclusion of interaction terms between the predictors and the social factor to the multivariable regression models.

Patient and public involvement

The development of the research question and outcome measures are based on the findings of our previous research into the health and developmental outcomes of MLP from infancy to 2 years, and other research.3–9 The importance of long-term outcomes was discussed in a workshop sponsored by our National Health and Medical Research Council-funded Centre of Research Excellence (CRE) in Newborn Medicine, which included parent representatives.11 While there was no direct involvement of MLP individuals or their families in the design of the study, the CRE in the Newborn Medicine Consumer Advisory Group, which includes individuals born preterm and their families, will be involved in designing a knowledge translation plan to disseminate study findings.

Ethics and dissemination

This study has ethical approval from the human research ethics committee at the Royal Children’s Hospital, Melbourne, Australia. Study outcomes will be presented through peer-reviewed publications and conference presentations. In addition, our results will be disseminated to key stakeholders, including clinicians, parents and teachers, facilitated by our CRE in the Newborn Medicine Consumer Advisory Group and partnerships with parent support groups.