Evaluating the Alimentary and Respiratory Tracts in Health and disease (EARTH) research programme: a protocol for prospective, longitudinal, controlled, observational studies in children with chronic disease at an Australian tertiary paediatric hospital.

INTRODUCTION: Chronic gastrointestinal and respiratory conditions of childhood can have long-lasting physical, psychosocial and economic effects on children and their families. Alterations in diet and intestinal and respiratory microbiomes may have important implications for physical and psychosocial health. Diet influences the intestinal microbiome and should be considered when exploring disease-specific alterations. The concepts of gut-brain and gut-lung axes provide novel perspectives for examining chronic childhood disease(s). We established the 'Evaluating the Alimentary and Respiratory Tracts in Health and disease' (EARTH) research programme to provide a structured, holistic evaluation of children with chronic gastrointestinal and/or respiratory conditions. METHODS AND ANALYSIS: The EARTH programme provides a framework for a series of prospective, longitudinal, controlled, observational studies (comprised of individual substudies), conducted at an Australian tertiary paediatric hospital (the methodology is applicable to other settings). Children with a chronic gastrointestinal and/or respiratory condition will be compared with age and gender matched healthy controls (HC) across a 12-month period. The following will be collected at baseline, 6 and 12 months: (i) stool, (ii) oropharyngeal swab/sputum, (iii) semi-quantitative food frequency questionnaire, (iv) details of disease symptomatology, (v) health-related quality of life and (vi) psychosocial factors. Data on the intestinal and respiratory microbiomes and diet will be compared between children with a condition and HC. Correlations between dietary intake (energy, macro-nutrients and micro-nutrients), intestinal and respiratory microbiomes within each group will be explored. Data on disease symptomatology, quality of life and psychosocial factors will be compared between condition and HC cohorts.Results will be hypothesis-generating and direct future focussed studies. There is future potential for direct translation into clinical care, as diet is a highly modifiable factor. ETHICS AND DISSEMINATION: Ethics approval: Sydney Children's Hospitals Network Human Research Ethics Committee (HREC/18/SCHN/26). Results will be presented at international conferences and published in peer-reviewed journals. TRIAL REGISTRATION NUMBER: NCT04071314. © Author(s) (or their employer(s)) 2020. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.

The prospective, longitudinal, controlled, observational design of this research programme provides a structured approach which can be simultaneously applied to multiple chronic gastrointestinal and/or respiratory conditions of childhood and utilises a universal control cohort (for age and gender matching).

This study will simultaneously evaluate dietary intake and the intestinal and respiratory microbiomes, which will tease out disease-causing alterations in the microbiomes, provide insights into the gut-lung axis and potentially identify modifiable dietary factors.

We will explore relationships between the primary outcomes (diet, intestinal and respiratory microbiomes) and health-related quality of life (including symptomatology), which may provide insights into the gut-brain axis and identify novel pathogenic mechanisms in these conditions.

A limitation of this research programme is that it currently includes a single centre, Sydney Children’s Hospital Randwick, Australia, however it is a tertiary referral centre for a diverse group of children across the state of New South Wales, Australia.

A further limitation is the arbitrary sample size targets given the exploratory nature of these studies.

Introduction

The primary disease burden in childhood has shifted over the last century from infectious to chronic diseases.1 Chronic childhood diseases, encompassing a wide spectrum of conditions with different pathogeneses, may have long-lasting physical, psychosocial and economic effects on children and their families.2 3 The human microbiome is a collection of all microorganisms (bacteria, viruses, archaea and eucaryotes) living in association with the human body.4 Our understanding of the human microbiomes in health and disease has begun to develop due to the advent of high-throughput sequencing and mass-spectrometry technologies, with the gut emerging as an ecosystem of particular interest. While the effects of an altered gut microbiome (dysbiosis) may not apply to all chronic diseases, there are conditions, disease-related complications and comorbidities linked to gut microbial dysbiosis. This is especially true in chronic gastrointestinal and respiratory conditions. Affected children are at risk of an imbalanced diet as well as mental health difficulties, which in turn can influence eating behaviours, attitudes and nutritional intake.5 6 Some of these conditions also require lifelong dietary modifications; for example, cystic fibrosis (CF).7 Additionally, the complex interaction between microbiota (ie, bacteria), available nutrients and the immune system is essential in maintaining homoeostasis and fighting against invading pathogens at mucosal sites.8 An important limitation common to most current publications on the human intestinal microbiome in chronic childhood disease(s) is the lack of quantifiable dietary data, as the diet has a marked influence on gut microbiota in health.9

The principles and framework of this research programme were developed to be applicable to many chronic gastrointestinal and/or respiratory conditions of childhood. Due to the clinical and/or research expertise of the authors and for the purposes of this manuscript, we will describe this programme based on three relevant chronic diseases: (i) CF, (ii) obstructive sleep apnoea (OSA) and Hirschsprung’s disease (HSCR). These conditions all have reported or expected changes in their intestinal and/or respiratory microbiomes.

CF is the most common life-shortening recessive disease in Caucasians.10 It is characterised by intestinal malabsorption, impaired growth and nutrition and lung disease.11 In CF, a high calorie, high fat diet (110% to 200% of recommended daily energy intake) is advised to prevent malnutrition and optimise growth.7 Recent reports suggest that children tend to achieve the recommended CF diet primarily by overconsumption of energy-dense, nutrient-poor foods rather than nutrient-dense foods.12 We have previously reported that children with CF, from as early as infancy, have alterations in their gut microbiota, impaired innate immunity and intestinal inflammation.13–17 We have also observed that poor growth in children with CF is significantly correlated with the degree of intestinal inflammation.14 The aetiology of gut microbial dysbiosis and inflammation in CF remains unclear. It is plausible that dietary intake plays a role, as enteric fat abundance (from a high-fat diet) may select for a pro-inflammatory microbiota.18 19 Alterations in intestinal metabolomic20 and proteomic21 profiles have also been reported. As the life expectancy of CF patients improves, age-related diseases such as gastrointestinal malignancies and cardiovascular disease (eg, myocardial infarcts in adults with CF) are a growing concern.22 Thus, optimal strategies to optimise health and reduce disease risk factors need to be determined.

In children, OSA can have cardiovascular, neurocognitive and behavioural consequences.23 Murine studies suggest intermittent hypoxia, hypercapnia and sleep fragmentation promote intestinal dysbiosis, increased visceral fat mass, systemic inflammation and atherosclerosis.24–27Additionally, the inhibition of gut microbial metabolites attenuating atherosclerosis26 and replication of hypertension after faecal transplant from hypertensive to normotensive OSA rats28 suggest the possibility of influencing clinical outcomes through affecting the gut microbiome. In adult studies, OSA is associated with gut epithelial damage,29and nasal dysbiosis and inflammation.30

HSCR is a congenital disorder where the distal intestine is aganglionic for a variable length. This results in a functional bowel obstruction that usually presents in newborns. Following corrective surgery, children often have ongoing intestinal symptoms, and Hirschsprung-associated enterocolitis (HAEC) remains the most frequent complication. This may result in frequent hospitalisations and even mortality. Children with and without HAEC often have an altered intestinal microbiome31 and altered composition of short chain fatty acids.32

To the best of our knowledge, there are no publications on the intestinal virome (ie, viruses) in children with CF, OSA or HSCR. Bacteriophages (viruses which infect bacteria) can influence bacterial populations via host lysis and horizontal gene transfer, as well as indirectly regulate immune function and inflammation.33–36

Despite accumulating evidence linking health, diet and the microbiomes, there is a paucity of research exploring this simultaneously in the context of chronic paediatric disease.37 Furthermore, potential gut-brain38 and gut-lung39 axes have yet to be well characterised in these conditions. The gut-brain axis refers to the bidirectional communication between the central nervous systems and gut microbiome, and is mediated by neural, endocrine and immune pathways.40 The gut-lung axis refers to the bidirectional relationship between the gut and lungs, as there appears to be an immunological relationship between them.39 41 Simultaneous, longitudinal studies using an integrated ‘omics’ approach will help to identify the functional consequences and pathogenic mechanisms that occur within the altered intestinal and respiratory milieu in chronic conditions. By exploring disease mechanisms and environmental interactions (eg, diet) we may in turn develop insights into potential therapeutic strategies. Additionally, we may be able to identify whether diet may be amenable to specific modifications which may in turn benefit the intestinal microbiome.

The EARTH programme has been established to provide a structured approach to analysing the gastrointestinal and respiratory microbiomes and diet in children with a chronic gastrointestinal and/or respiratory condition. The design improves efficiency by recruiting and assessing a healthy control (HC) group which can be used for comparison against each of the conditions (as opposed to recruiting a new HC group for each condition). Although our initial design is focussed on CF, OSA and HSCR, the programme framework is applicable to other chronic gastrointestinal and/or respiratory conditions of childhood.

Objectives

The objective of this research programme is to evaluate and compare children with a chronic gastrointestinal and/or respiratory condition and age and gender matched HC. The primary objectives include analysing the intestinal and respiratory microbiomes (using an integrated ‘omics’ approach) and dietary intake using validated, parent-report tools (table 1). The secondary objectives are also presented in table 1 and include evaluating:

Table 1

Primary and secondary objectives with related outcome measures

Domain	Data source	Technique	Outcome measures	Between group analyses*	Within group analyses†
Primary objectives
(1) Intestinal microbiome	(1) Stool sample	Bacterial communities (16S rRNA (V4)13or MSS)	Alpha diversity (richness and Shannon index)	Student’s t-test or Wilcoxon signed-rank test	Pearson or Spearman correlations with: Gastrointestinal microbiome Respiratory microbiome Diet Secondary objectives Descriptive
			Beta diversity (UniFrac) distances46	PERMANOVA73
			Relative abundances of bacteria	ANCOM74
		Viral communities (metagenomic sequencing49	Alpha diversity (richness and Shannon index)	Student’s t-test or Wilcoxon signed-rank test
			Beta diversity (Bray-Curtis dissimilarities)	PERMANOVA73
(2) Respiratory microbiome	(2) Oropharyngeal swab or sputum sample		Relative abundances of viruses	ANCOM74
		Proteomics (LC-MS)21	Protein z-score normalised LFQ intensities	Student’s t-test
			Pathway/network upregulation or downregulation	Condition/HC ratio
		Metabolomics (UHPLC-MS/MS)20	Metabolite normalised abundance	Student’s t-test
			Pathway/network upregulation or downregulation	Condition/HC ratio
(3) Diet	(i) ACAES (ages 2 to 18 year)(ii) 24 hours food recall(ages 0 up to 2 year)		Energy intake	Student’s t-test, Wilcoxon signed-rank test or Fisher’s exact test
			Percent energy from core foods
			Macronutrient intake
			Micronutrient intake
			Diet quality score‡
Secondary objectives
(1) Biomarkers	Stool, oropharyngeal swab or sputum sample	ELISA	Inflammation (calprotectin, M2-PK, CRP and interleukins)	Student’s t-test, Wilcoxon signed-rank test or Fisher’s exact test	Descriptive
(2) Symptomatology and HRQOL	PedsQL Infant Scales (0–2 year) and gastrointestinal symptoms module (2–18 year)54–56		HRQOL and gastrointestinal symptoms	Student’s t-test, Wilcoxon signed-rank test or Fisher’s exact test	Descriptive
	Rome IV Questionnaire57 58		Gastrointestinal symptoms
	Spence Children’s Anxiety Scale59 60		Anxiety symptoms
	Short Mood and Feelings Questionnaires61 62		Depressive symptoms
(3) Phenotypic and clinical information	Anthropometrics		Z-scores; weight, length/height, weight-for-length (ages 0 to 2 year) and BMI (ages 2 to 20 year)	Student’s t-test, Wilcoxon signed-rank test or Fisher’s exact test	Descriptive
	Clinical presentations		Number and length of hospitalisations, emergency department presentations, medications, vaccinations
	Results		Biochemistry, microbiology and imaging results
	Perinatal factors		Mode of delivery, feeding during infancy
(4) Socio-demographic factors			Ethnicity	Descriptive	Descriptive
			SEIFA code63	Descriptive	Descriptive

All samples, questionnaires and data will be collected from all participants at each time point.

*Between group analyses describe comparisons between a condition and healthy control groups.

†Within group analyses describe analyses of two outcome measures within subjects of the same condition group.

‡ACAES only.

ACAES, Australian Child and Adolescent Eating Survey; ANCOM, analysis of composition of microbiomes; BMI, body mass index; CRP, C-reactive protein; HC, healthy control; HRQOL, health-related quality of life; LC-MS, liquid chromatography-mass spectrometry; LFQ, label-free quantification; M2-PK, M2 pyruvate kinase; MSS, metagenomic shotgun sequencing; PERMANOVA, permutational multivariate analysis of variance; rRNA, Ribosomal nucleic acid; SEIFA, socio-economic indexes for areas (a measure of relative socio-economic advantage and disadvantage in Australia); UHPLC-MS/MS, ultra-high performance liquid chromatography-tandem mass spectrometry; V4, Hypervariable region V4.

Known inflammatory biomarkers.

Symptomatology and health-related quality of life (HRQOL) using validated measures.

Phenotypic and clinical information.

Sociodemographic factors.

Additional secondary objectives include correlating within children with the same condition: (i) dietary intake with the intestinal microbiome, (ii) dietary intake with the respiratory microbiome and (iii) the intestinal and respiratory microbiomes.

We hypothesise that:

Children with chronic gastrointestinal and/or respiratory conditions will have altered intestinal and respiratory microbiomes compared with healthy children, and

Diet plays a key role in influencing the intestinal and respiratory microbiomes and this may impact on clinical outcomes, biomarkers of disease and health-related quality of life.

To our knowledge, this programme will enable the first series of studies comparing the intestinal and respiratory microbiomes and diet in children with chronic gastrointestinal and/or respiratory conditions. Initial results will be hypothesis-generating and used to direct future studies tailored to a specific focus or line of inquiry. Additionally, studies from this research programme have potential for direct translation into clinical care as diet is a highly modifiable factor.

Methods and analyses

Study design

The EARTH programme provides a framework for a series of prospective, longitudinal, controlled, observational studies, with each individual study comparing children with a chronic gastrointestinal and/or respiratory condition to HC. A single healthy control group will be used for comparison against all conditions and healthy controls are defined as children who are free of any chronic disease. The standardised methodological approach will also allow for comparisons between different health conditions. The Standard Protocol Items: Recommendations for Interventional Trials reporting guidelines were used for this protocol.42

Setting

Studies will be carried out at a single centre; the Sydney Children’s Hospital (SCH) in Randwick, Australia. SCH is a tertiary paediatric hospital.

Participants

Children are eligible if they:

Are aged between 0 and 18 years;

Have been diagnosed with a chronic gastrointestinal and/or respiratory condition defined by consensus diagnostic criteria; or

Are free of any chronic health condition (healthy control group); and

Have a parent(s)/carer(s) who provides informed consent, or are at least 16 years old and provide informed consent.

Ineligibility criteria include:

Children with more than one concurrent or unrelated chronic disease;

Inability to comply with study requirements;

Parent(s)/guardian(s) are unable to speak English or do not have a reading level age of at least 12 years.

Participants with a chronic gastrointestinal and/or respiratory condition will be matched to a HC for gender and age (as closely as possible).

Recruitment strategy

Participants with chronic gastrointestinal and/or respiratory conditions will be approached at their routine clinical appointments in the outpatient department. Flyers will be placed in the hospital for recruitment of HC. Prior to study participation, detailed written and verbal information will be provided about the content and extent of the study. Written informed consent from the parent/legal guardian of each participant will be required. If the child is deemed Gillick competent,43 they will be encouraged to sign a specific child assent form. Parents/legal guardians and participants may withdraw consent at any time.

Outcome measures

The outcomes measures are presented in table 1. All samples, questionnaires and data will be collected from all participants at each time point. Presented below is a simplified explanation of each outcome/variable included in the research programme.

Primary outcomes/variables

1. Intestinal microbiome assessed from a stool sample using one or more of:

Bacterial community analysis (16S ribosomal RNA (V4)13 or metagenomic shotgun sequencing):

Alpha diversity indices:

Richness: the total number of unique species.44

Shannon index: a measure of both species abundance and evenness.45

Beta diversity indices:

UniFrac: a distance metric used to compare biological communities that incorporates phylogenetic distances between observed organisms.46

Bray-Curtis dissimilarity: a count metric used to quantify the compositional dissimilarity between two different sites.47

Relative abundance: the per cent composition of an organism relative to the total number of organisms in the area.48

Viral community analysis (metagenomic sequencing),49 as above for bacterial community analysis.

Proteomics (liquid chromatography-mass spectrometry (LC-MS)):21

Protein z-score normalised label-free quantification intensities.

Pathway/network upregulation or downregulation based on the ratio of condition/HC.

Metabolomics (ultra-high performance liquid chromatography-tandem mass spectrometry (UHPLC-MS/MS)):20

Metabolite normalised abundance.

Pathway/network upregulation or downregulation based on the ratio of condition/HC.

2. Respiratory microbiome assessed from an oropharyngeal swab or sputum sample, using one or more of the techniques listed above (1a–d).

A sputum sample will be obtained in children able to expectorate and an oropharyngeal swab will be collected in children unable to expectorate.

Associations with the intestinal microbiome (1) will be used to explore the gut-lung axis.

3. Dietary intake assessed using a validated semi-quantitative food frequency questionnaire:

For participates aged 0 up to 2 years, a 24 hours food recall:

Energy intake, percentage energy from core foods, macronutrient intake (total intake and proportion of energy intake) and micronutrient intake (total intake and proportion of energy intake).

For participants aged 2 to 18 years, the Australian Child and Adolescent Eating Survey (ACAES).50–53 The ACAES is a validated food frequency questionnaire (120 items, semi-quantitative) used to quantify food and nutrient intake over the preceding 6 months, developed and validated for use in Australian children:

Energy intake, percentage energy from core foods, macronutrient intake (total intake and proportion of energy intake), micronutrient intake (total intake and proportion of energy intake) and overall diet quality score.

Faecal and respiratory inflammatory biomarkers, such as calprotectin, M2 pyruvate kinase, C-reactive protein and interleukins.

Symptomatology and health-related quality of life (HRQOL) will be collected directly from children where age-appropriate measures exist and/or parents using age-appropriate measures:

PedsQL Infant Scales (0 to 2 year) and gastrointestinal symptoms module (2 to 18 year)54–56 (HRQOL and gastrointestinal symptoms).

Rome IV Questionnaire57 58 (gastrointestinal symptoms). Designed to diagnose functional gastrointestinal disorders, which are defined as disorders of the gut-brain interaction in children aged 0 to 18 years. These criteria capture gastrointestinal symptoms which are relevant to motility disturbance, visceral hypersensitivity, altered mucosal and immune function, altered gut bacteria and altered central nervous system processing;

Associations with the intestinal microbiome (1) will be used to explore the gut-brain axis.

Spence Children’s Anxiety Scale59 60 (anxiety symptoms);

Associations with the intestinal microbiome (1) will be used to explore the gut-brain axis.

Short Mood and Feelings Questionnaires61 62 (depressive symptoms).

Associations with the intestinal microbiome (1) will be used to explore the gut-brain axis.

Anthropometrics, including z-scores for weight, length/height, weight-for-length (ages 0 to 2 years) and body mass index (BMI) (ages 2 to 20 years).

Z-scores; weight, length/height, weight-for-length (ages 0 to 2 year) and BMI (ages 2 to 20 year).

Clinical information and biochemical results obtained through routine care, such as number and length of hospitalisations, emergency department presentations, perinatal factors (mode of delivery, feeding type(s) in infancy), medications, vaccination status (including timing of most recent vaccination), biochemistry, microbiology and imaging results.

Sociodemographic factors such as ethnicity and socio-economic indexes for areas code.63

Procedures

Each participant will be assessed on three occasions over a 12-month period; at study entry, 6 and 12 month follow-up. At each time point, the following will be collected:

A stool sample;

An oropharyngeal swab or sputum sample (a sputum sample will be obtained in children able to expectorate and an oropharyngeal swab will be collected in children unable to expectorate);

Dietary intake measured using the ACAES (2 to 18 years)or 24 hours food recall (0 up to 2 years);

A secure, password-protected online survey comprising:

PedsQL Infant Scales (0 to 2 year) and gastrointestinal symptoms module (2 to 18 year),54–56 tailored to age;

Rome IV Questionnaire57 58 (0 to 18 years);

Spence Children’s Anxiety Scale59 60 (3 to 18 years);

Short Mood and Feelings Questionnaires61 62 (6 to 18 years);

Clinical and biochemical results obtained through routine care and hospitalisations (if available);

Sociodemographic factors (baseline survey only);

Anthropometrics: height, weight and BMI z-scores.

Details regarding sample and data collection

Stool samples will be collected in a sterile specimen jar using a Feces Catcher (Abbexa Ltd, Cambridge, UK). Sputum samples will be collected from children who are able to expectorate and oropharyngeal swabs will be collected in those children who cannot expectorate sputum. All samples will be transported in a cooler bag (with a −18°C ice pack) to the hospital laboratory within 24 hours of collection. Stool, oropharyngeal swab and sputum samples will be aliquoted and stored at −80°C. Where practical, samples will be collected at least 4 weeks after completion of oral/intravenous antibiotic therapy (excluding prophylactic antibiotics) or an acute infectious illness.

For children aged 0 up to 2 years, a 24 hours dietary recall will be conducted by a study dietician. For children aged 2 to 18 years, the ACAES will be completed online by parents/guardians or by the child themselves if aged over 14 years.

Participant clinical information and questionnaires will be collected using the Qualtrics online survey software tool (www.qualtrics.com), which is distributed via a secure email link. Qualtrics is a secure, password-protected platform, which allows for the distribution of electronic surveys and collection of data. The surveys are programmed to be condition/control-specific and facilitate a personalised flow depending on the time point of the study and the age of the participant (to facilitate administration of age-appropriate questionnaires). Further details on each of the measures used for evaluating symptomatology and HRQOL are presented in table 2. Clinical and biochemical results will be obtained via the SCH electronic medical record system and recorded in a Qualtrics survey.

Table 2

Measures for symptomatology and health-related quality of life (HRQOL)

Measure	Domains (Items)	Scoring	Interpretation
(i) PedsQL
Infant Scales - parent report for infants (ages 1–12 months)56	Total (36): 5-point LS.Physical functioning (6), physical symptoms (10), emotional functioning (12), social functioning (4), cognitive functioning (4).	Items are reverse scored and linearly transformed on a scale from 0 to 100.	Higher scores indicate better HRQOL
Infant Scales - parent report for infants (ages 13–24 months)56	Total (45): 5-point LS.Physical functioning (9), physical symptoms (10), emotional functioning (12), social functioning (5), cognitive functioning (9).	Items are reverse scored and linearly transformed on a scale from 0 to 100.	Higher scores indicate better HRQOL
3.0 Gastrointestinal symptoms module – parent report for toddlers (ages 2–4)54 55	Total (74): 5-point LS.Stomach pain and hurt (6), stomach discomfort when eating (5), food and drink limits (6), trouble swallowing (3), heartburn and reflux (4), nausea and vomiting (4), gas and bloating (7), constipation (14), blood in poop (2), diarrhoea (7), worry about going poop (5), worry about stomach aches (2), medicines (4), communication (5).	Items are reverse scored and linearly transformed on a scale from 0 to 100.	Higher scores indicate lower problems.
Gastrointestinal symptoms module (acute V.3.0) – parent report for young children (ages 5–7)54 55	Total (74): 3-point and 5-point LS.Stomach pain and hurt (6), stomach discomfort when eating (5), food and drink limits (6), trouble swallowing (3), heartburn and reflux (4), nausea and vomiting (4), gas and bloating (7), constipation (14), blood in poop (2), diarrhoea (7), worry about going poop (5), worry about stomach aches (2), medicines (4), communication (5).	Items are reverse scored and linearly transformed on a scale from 0 to 100.	Higher scores indicate lower problems.
Gastrointestinal symptoms module (acute V.3.0) – young child report (ages 5–7)54 55
Gastrointestinal symptoms module (acute V.3.0) – parent report for children (ages 8–12)54 55	Total (74): 3-point and 5-point LS.Stomach pain and hurt (6), stomach discomfort when eating (5), food and drink limits (6), trouble swallowing (3), heartburn and reflux (4), nausea and vomiting (4), gas and bloating (7), constipation (14), blood in poop (bowel movement) (2), diarrhoea (7), worry about going poop (bowel movements) (5), worry about stomach aches (2), medicines (4), communication (5).	Items are reverse scored and linearly transformed on a scale from 0 to 100.	Higher scores indicate lower problems.
Gastrointestinal symptoms module (acute V.3.0) – child report (ages 8–12)54 55
Gastrointestinal symptoms module (acute V.3.0) – parent report for teens (ages 13–18)54 55
Gastrointestinal symptoms module (acute V.3.0) – teens report (ages 13–18)54 55
(ii) Rome IV
Rome IV – parent-report form for infants and toddlers (ages 0–3) (R49QG-toddler)57	Total (29 for ages 0–12 months; 18 for ages 1–3 years):Infant gastrointestinal problems (11), vomiting (9), bowel movements (9)	Defined diagnostic criteria for functional gastrointestinal disorders in neonates and toddlers:57Infant regurgitation, infant rumination syndrome, cyclic vomiting syndrome, infant colic, functional diarrhoea, infant dyschezia, functional constipation.
Parent-report form for children and adolescents (4 years of age and older) (R4PDQ-child)58	Total (42):Belly ache and uncomfortable feelings above the belly button (12), belly aches and abdominal pain around and below the belly button (10), bowel movements (7), nausea and vomiting (9), other symptoms (4).	Defined diagnostic criteria for functional gastrointestinal disorders in children and adolescents:58Cyclic vomiting syndrome, functional nausea and functional vomiting, rumination syndrome, aerophagia, functional dyspepsia, irritable bowel syndrome, abdominal migraine, functional abdominal pain – not otherwise specified, functional constipation, non-retentive faecal incontinence.
Self-report form for children and adolescents (10 years of age and older) (R4PDQ-child)58
(iii) Spence Children’s Anxiety Scale
Spence – Preschool Anxiety Scale (parent report) (ages 0 to 4)59 60	Total (34): 5-point LS.Generalised anxiety (5), social anxiety (6), obsessive compulsive disorder (5), physical injury fears (7), separation anxiety (5).	Responses are scored 0 (not true at all) to 4 (very often true). A maximum possible score of 112.	A score 1 SD above mean for a subscale or total score warrants further clinical investigation. A score of 0.5 SD above the mean on total score is indicative of an elevated, but not clinical level of anxiety.
Spence Children’s Anxiety Scale (parent report) (5 years and older)59 60	Total (38 scored, 39 total): 4-point LS.Panic attack and agoraphobia (9), separation anxiety (6), physical injury fears (5), social phobia (6), obsessive compulsive (6), generalised anxiety disorder/overanxious disorder (6).	Responses are scored 0 (never) to 3 (always). A maximum possible score of 114. T-score calculation.	A score 1 SD above mean (T-score of ≥60) for a subscale or total score is indicative of subclinical or elevated levels of anxiety warranting further clinical investigation.
Spence Children’s Anxiety Scale (8 years and older)59 60	Total (38 scores, 45 total): 4-point LS.Separation anxiety (6), social phobia (6), obsessive compulsive (6), panic attack and agoraphobia (9), physical injury fears (5), generalised anxiety (6).
(iv) Short Mood and Feelings Questionnaire
Mood and Feelings Questionnaire: Short Version (parent report on Child) (ages 6–18)61 62 75	Total (13): 3-point LS.Depressive symptoms (13).	Responses are scored 0 (not true) to 2 (true). A maximum possible score of 26.	Higher scores suggest more severe depressive symptoms. A score of ≥12 may indicate the presence of depression in the respondent.
Mood and Feelings Questionnaire: Short Version (child self-report) (ages 6–18)61 62 75	Total (13): 3-point LS.Depressive symptoms (13).

All questionnaires will be collected from all participants at each time point.

N.B. Disease-specific questionnaires can be added into the Qualtrics data collection form, that is, the Paediatric Sleep Questionnaire: Sleep-Disordered Breathing Subscale,76 for children with OSA.

HRQOL, health-related quality of life; LS, Likert scale; OSA, obstructive sleep apnoea.

Participants will be measured for their height and weight using standardised methods.

Sample and data processing techniques

Processing of stool, oropharyngeal swab and sputum samples is almost identical (sparing a few initial sample preparation steps). For bacterial community analysis, DNA will be extracted using QIAamp DNA kits (QIAGEN, Hilden, Germany) according to manufacturer’s instructions. For 16S rRNA gene analysis specifically, amplification will be performed with primers 515F and 806R spanning the V4 region and sequencing data will be processed using USEARCH.64

In the instance where species resolution of bacterial communities is thought to be beneficial, metagenomic shotgun sequencing (MSS) will be performed as an alternative to 16S rRNA gene sequencing. For MSS, no amplification step will be performed prior to sequencing. Sequencing data will be processed using a custom in-house pipeline.

For viral community analysis specifically, sample preparation will follow an adjusted NetoVIR (Novel Enrichment Technique Of VIRomes) protocol.49 All sequencing will be performed using the Illumina MiSeq platform at the Ramaciotti Centre for Genomics at the University of New South Wales (UNSW). Briefly, sequencing data will be processed using the Vipie platform65 for taxonomic assignment and VirSorter pipeline66 for functional annotation.

For untargeted proteomics, samples will undergo an adjusted Debyser et al protocol for protein extraction, gel electrophoresis and analysed using LC-MS/MS at the Bioanalytical Mass Spectrometry Facility (BMSF), UNSW.21 Briefly, proteomics data will be analysed using MaxQuant67 and Ingenuity Pathway Analysis (QIAGEN).

For untargeted metabolomics, metabolites will be extracted in 1:1 (v:v) acetonitrile:H20 and analysed using a U3000 UHPLC system coupled to a Q-Exactive mass spectrometer (MS; Thermo Fisher Scientific) at the BMSF, UNSW. Briefly, metabolomics data will be analysed using Progenesis CoMet (Waters/NonLinear Dynamics).

Faecal and respiratory biomarkers (listed above) will be measured using ELISA.

Nutrient intake data from the ACAES and 24 hours recall is computed using FoodWorks (V.3.02.581) and the following databases: Australian AusNut 1999 database (All Foods) Revision 14 and AusFoods (Brands) Revision 5 (Xyris Software (Australia) Pty Ltd, FoodWorks Professional V.3.02.581. 2004: Brisbane Australia). Outputs include a quantified estimate and the percentage of energy from a wide of range of macro-nutrients (protein, fat, carbohydrate) and micro-nutrients (vitamins A, B, C and minerals such as iron, zinc and calcium). In addition, overall diet quality score and the percentage of energy derived from nutrient rich core foods and energy-dense, nutrient-poor discretionary foods is calculated.

Administration of patient records and data

At the time of consent and enrolment, participants will be assigned a unique study ID number (nine alphanumeric characters). All patient records, samples and data are de-identified using the unique study ID. Data will be stored securely as per ethics review board guidelines.

Handling of abnormal outcomes or distress

The well-being of participants is of utmost importance. Participants and their parents/guardians will be advised to contact any of the study investigators if they have concerns regarding any aspects of their participation. It is possible that thinking about one’s health or the health of one’s child may elicit emotional distress in some participants. Depending on the nature of the concern or level of distress communicated, a relevant study investigator will contact the participant and/or his or her primary caregiver by telephone or in-person to assess any concerns and arrange appropriate follow-up or referral as soon as possible. Additionally, all Participant Information Sheets will provide the details for several, free, age-appropriate 24 hours telephone-based support services. All individuals will be clearly informed that choosing not to take part in the study, or withdrawing from the study at any stage, will not adversely affect their or their child’s healthcare or relationship with hospital staff in any way.

Bias, confounding factors and handling of missing data

The single-centre nature of this study is a limitation due to the restricted recruitment pool available and potential for selection bias; however, SCH is a tertiary referral centre for a diverse group of children across the state of New South Wales, which is the most populous state in Australia. Age and gender are known confounding factors for microbiome analyses and are controlled for with matching. There are rapid changes in the intestinal microbiota during the first 3 years of life, after which it becomes relatively stable.68–70Although we aim to match participants as closely as possible, our criteria for acceptable matching is as follows:

For children less than 1 year, age to be matched within 3 months,

For children aged 1 to 3 years, age to be matched within 6 months,

For children aged 4 years and older, age to be matched within 2 years.

Additional confounding factors for microbiome analyses include perinatal factors and ethnicity, for which sensitivity analyses will be performed. Condition specific medications (eg, pancreatic enzyme replacement therapy or antibiotic therapy in CF) are potential confounders and attempts to control for these factors will be made at the analysis stage. Vaccines may influence short-term cytokine production and analyses as a confounder will be considered.71 72 Missing data will be treated as missing and accounted for using linear mixed models (see statistical methods below).

Study size

In an exploratory research programme of this nature, with multiple conditions of interest, sample size calculations for the primary outcomes are difficult. As an initial, arbitrary target, three males and three females in each of the following age ranges (0 to 5, >5 to 10, >10 to 18 years) will be recruited to account for age-related and gender-related changes in microbiomes and diet. This calculation assumes that six participants will be required for most statistical tests of interest and an analysis can be performed on the smallest subgroup (eg, six CF vs six HC children aged 0 to 5 years). Therefore 18 participants for each condition and 18 HC (which can be used for comparison against multiple conditions) are an initial target sample size. Initial data from this sample size can then be used for subsequent power-focussed study designs.

Statistical methods

Statistical analyses will be performed in R V.3.4.4. All outcome measures will be analysed cross-sectionally and temporally. Descriptive statistics will be calculated for all outcome parameters for each cohort according to normality of distribution. Given the sample size, categorical variables will be compared using Fisher's exact test. Continuous variables will be analysed according to distribution with a Student's t-test or Wilcoxon signed-rank test for parametric and non-parametric data, respectively. A linear random-effects mixed model or variant of generalised linear-mixed model will be used to evaluate cross-sectional and temporal differences in outcome measures. This technique will allow for control of confounders and treatment of missing data as missing. Correlations between two continuous variables will be performed using Pearson or Spearman correlations according to distribution. Alpha diversity indices will be measured by richness (number of taxa) and Shannon index. Phylogeny-based and taxonomy-based beta diversity will be calculated using UniFrac distances46 and Bray-Curtis dissimilarities, respectively, and used to generate non-metric multidimensional scaling plots. Permutational multivariate analysis of variance tests (permutations=1000) will be used to test if beta diversity significantly differs between groups and age using the vegan function adonis.73 A significant difference in abundance of taxa, proteins or metabolites between groups will be assessed using the analysis of composition of microbiomes package V.1.1–3.74 For all analyses, p<0.05 (two-tailed) is considered significant except in the instance of multiple comparisons, in which case a Benjamini & Hochberg correction will be applied and q<0.05 will be considered significant.

Ethics and dissemination

The EARTH research programme received ethics approval from the Sydney Children’s Hospitals Network Human Research Ethics Committee (HREC/18/SCHN/26). Any amendment to the protocol which may impact the conduct of the study will be approved by the ethics committee before implementation.

The results of studies from this research programme will be presented in international conferences and will be published in peer-reviewed journals. Findings may also be presented as: (i) easy-to-read summaries for participants and the community; (ii) educational lectures and seminars for patients, families and the community; (iii) website and social media postings; (iv) newsletter updates for study participants; (v) reports for relevant advocacy groups and funding partners.

Expected outcomes and significance of the research project

To our knowledge the EARTH research programme will be the first in children with a chronic gastrointestinal and/or respiratory condition to simultaneously evaluate dietary intake and the intestinal and respiratory microbiomes. By exploring disease mechanisms and environmental interactions (ie, diet) we may in turn develop insights into potential therapeutic strategies. Studies from this programme have the potential for direct translation into clinical care as diet is a highly modifiable factor. This programme also provides a structured approach for performing prospective, longitudinal, controlled, observational studies which can be simultaneously applied to multiple health conditions, and utilised a universal control cohort.