- Study protocol
- Open Access
- Open Peer Review
A comparison of intestinal microbiota in a population of low-risk infants exposed and not exposed to intrapartum antibiotics: The Baby & Microbiota of the Intestine cohort study protocol
BMC Pediatricsvolume 16, Article number: 183 (2016)
The intestinal microbiota influences metabolic, nutritional, and immunologic processes and has been associated with a broad range of adverse health outcomes including asthma, obesity and Type 2 diabetes. Early life exposures may alter the course of gut microbial colonization leading to differences in metabolic and immune regulation throughout life. Although approximately 50 % of low-risk full-term infants born in Canada are exposed to intrapartum antibiotics, little is known about the influence of this common prophylactic treatment on the developing neonatal intestinal microbiota. The purpose of this study is to describe the intestinal microbiome over the first 3 years of life among healthy, breastfed infants born to women with low-risk pregnancies at full term gestation and to determine if at 1 year of age, the intestinal microbiome of infants exposed to intrapartum antibiotics differs in type and quantity from the infants that are not exposed.
A prospectively followed cohort of 240 mother-infant pairs will be formed by enrolling eligible pregnant women from midwifery practices in the City of Hamilton and surrounding area in Ontario, Canada. Participants will be followed until the age of 3 years. Women are eligible to participate in the study if they are considered to be low-risk, planning a vaginal birth and able to communicate in English. Women are excluded if they have a multiple pregnancy or a preterm birth. Study questionnaires are completed, anthropometric measures are taken and biological samples are acquired including eight infant stool samples between 3 days and 3 years of age.
Our experience to date indicates that midwifery practices and clients are keen to participate in this research. The midwifery client population is likely to have high rates of breastfeeding and low rates of intervention, allowing us to examine the comparative development of the microbiome in a relatively healthy and homogenous population. Results from this study will make an important contribution to the growing understanding of the patterns of intestinal microbiome colonization in the early years of life and may have implications for best practices to support the establishment of the microbiome at birth.
The intestinal microbiota is essential to metabolic, nutritional and immunologic processes. The microbiome of healthy adults varies significantly between individuals. Differences in diversity or relative abundance of microbes has been associated with a broad range of adverse health outcomes including: obesity, inflammatory bowel disease, atopic disease, Type 2 diabetes,  multiple sclerosis,  mental health or behaviour disorders, including autism  and alcoholism . The adult microbiota can be transiently altered by diet, infection and antibiotic use, but in most individuals it reverts to the composition observed before the perturbation [5–8].
In the newborn infant, following birth, colonisation of the gut typically proceeds rapidly in the first hours and days following initial exposure to maternal vaginal and rectal microbiota and subsequently to maternal skin and breast milk [9, 10]. Compared to adults, the infant microbiome is quite unstable, has greater inter- and intra-individual variation, has a lower number of species and has a higher proportion of Bifidobacterium . Microbial communities in the young child stabilize over time and begin to resemble those in adults between 1 and 3 years of age . Factors such as the timing of weaning from breast milk as well as the type and timing of the introduction of solid foods may contribute to this process, [13, 14] but require further exploration. Membership of microbial communities and colonization patterns has been shown to differ by mode of birth. Caesarean section and exposure of the newborn to antibiotics following birth cause deviation from normal colonization patterns with disturbances in and poor diversity of intestinal microbiota and these differences have been shown to persist for up to 1 year [15–18]. Caesarean section has been associated with a number of long-term health outcomes and the gut microbiome has been implicated in the biological pathway [19, 20]. Thus, early exposure to the appropriate colonizing organisms is potentially critical to long-term health.
As many as 50 % of low-risk, full term infants born in Canada are exposed to intrapartum antibiotic prophylaxis (IAP); about 30 % because of prophylactic management for Group B Streptococcus (GBS)  and the remainder due to antibiotic coverage used for Caesarean birth. IAP has the potential to impact early development of neonatal intestinal microbiota, first because of the changed maternal vaginal and rectal microbiota to which the neonate is exposed during birth, and secondly due to the direct exposure of the fetus to the antibiotic. Therefore, the primary objective of our study is to determine if infants born to women who receive IAP for GBS or a similar regimen have intestinal microbiota at 1 year that differs significantly in type and quantity from those not exposed to IAP (Fig. 1).
Study design and setting
A prospectively followed cohort of 240 mother-infant pairs will be formed by enrolling eligible pregnant women from midwifery practices in the City of Hamilton and surrounding area in Ontario, Canada. A pilot cohort (83 mother-infant pairs) was enrolled between July 1st 2012 and December 31st 2013. Recruitment for the second phase of the study began October 1st 2014 and is ongoing. Participants will be followed until the age of 3 years. The study was approved by the joint Hamilton Health Sciences - McMaster University Research Ethics Board and by Research Ethics Boards at all participating healthcare organizations.
Women are eligible to participate in the study if they are considered to be low-risk (defined as being under the care of a midwife), planning a vaginal birth and able to communicate in English in order to provide signed informed consent. Women are excluded from the study if they have a known multiple pregnancy or a preterm birth (before 37 completed weeks gestation). Women who are enrolled prior to 37 completed weeks gestation are contacted by telephone after 37 completed weeks to reassess their eligibility.
The Baby & Mi study is advertised through pamphlets and posters in the waiting areas at eight participating midwifery practices and on social media. The midwives and administrative staff are asked to inform their pregnant clients about the study. ‘Consent to Contact’ forms are completed by mothers interested in hearing more about the study and research personnel then contact them by telephone or email. Study visits are arranged with women who agree to participate.
After providing informed consent, women are asked to complete a baseline questionnaire that collects information about their pregnancy and family medical history. Women enrolled at a gestational age less than 37 completed weeks are contacted by telephone within 3 weeks of their estimated due date to reassess time-dependent characteristics such as smoking status and medication use during pregnancy. Case report forms are available in Additional file 1.
Data about the pregnancy and birth are recorded onto paper case report forms from antenatal forms, the birth record and mother and infant charts by the midwives. On day 3, day 10 and 6 weeks postpartum participants are instructed to complete follow up questionnaires and to place a diaper containing stool in a labeled study bag and store the sample frozen in their household freezer until their study visit when the baby is 12 weeks of age. The women are reminded to complete these tasks on the appropriate days by research personnel through email or text message. In order to facilitate recall throughout the follow up period, mothers are given a study diary that provides them with a place to record events such as infections, medications and sleep disturbances.
Participants are asked to attend follow up visits at McMaster University Medical Centre at 12 weeks, 5 months, 1 year, 2 years and 3 years of age for stool sample delivery, follow up data collection and measurements of infant growth. Participants of the pilot study were also asked to attend a follow up visit at 6 weeks postpartum. Women are contacted between study visits to complete questionnaires online or by telephone. A summary of the data collection time points can be found in Table 1.
Infants born to women who are administered prophylactic antibiotic treatment for GBS or a similar regimen in labour will be classified as IAP exposed. The unexposed group includes infants born to women who are not exposed to IAP. Women received IAP under the recommendations of their care providers according to standard care practices. Information about intrapartum antibiotic use is collected from the birth record including indication, dose, frequency and time to delivery from first dose.
The primary outcome is the type and abundance of species present in the intestinal microbiota at 1 year of age, which will be determined through bacterial diversity analysis and microbiota profiling of stool specimens This comparison will also be made using stool samples collected at day 3 and 10, 6 and 12 weeks, 5 months, 2 and 3 years of age.
Gut microbiome measurement
DNA extraction from stool is carried out using a previously described protocol that enhances DNA recovery from microbial communities [22, 23] with modifications to increase quantitative recovery of bacteria across different taxa . Upon delivery to the laboratory, stool samples are thawed and DNA is extracted. Approximately 100 mg of stool is sampled from each diaper. Where the stool is fully absorbed into the diaper, a 1 cm by 1 cm square of diaper is cut using sterile scissors and used for DNA extraction. Bacterial community profiling of the 16S rRNA gene is carried out using paired end reads of the V3 region using barcoded Illumina sequencing as described previously . Sample preparation and 250 paired-end sequencing is carried out on a MiSeq Illumina sequencer as per manufacturer’s instructions. This provides approximately 50,000–100,000 reads per sample that is processed by an in-house bioinformatics pipeline  and the output includes clustered sequences in operational taxonomic units (OTUs) and taxonomic assignments as described previously [26, 27].
Growth is monitored at each follow-up visit using standard approaches for the measurements of length, weight (including weight to length ratio and peak weight velocity), skinfold thickness (tricep, subscapular, bicep and suprailiac), head and hip circumference. The rate of fat accretion will be determined using measurements of body composition at 6 weeks (pilot participants only), 12 weeks, 5 months and 3 years. At 6 weeks, 12 weeks and 5 months of age an air displacement plethysmography system (PEA POD) is used to measure percent fat. At 3 years of age body composition is measured using dual-energy X-ray absorptiometry (DXA).
Atopic disease is assessed during study visits and telephone contacts by asking mothers if they suspect or have been told by a physician that their child has eczema, asthma, reactive airway disease or allergies. Symptoms of eczema are also assessed using the Sampson Oranje criteria for eczema  and the core questionnaire for eczema that was developed by the International Study of Asthma and Allergies in Childhood . Sleep patterns are reported at each contact using questions from the Brief Index Sleep Questionnaire . Symptomatic gastroesophageal reflux is evaluated at 6 weeks, 12 weeks, 5 months and 1 year using the Infant Gastroesophageal Reflux Questionnaire Revised . Gastrointestinal problems are measured at 30 months using the PedsQL-Gastrointestinal Symptoms Scale . Measures of cardiometabolic health and systemic inflammation will be made using fasting blood collected from each participant at 3 years of age. This will include lipids, glucose and inflammatory markers. Fecal calprotectin will be measured at 3 years of age as a marker for gut inflammation. Child behaviour and temperament will be measured using the Strengths and Difficulties  and the Children’s Behaviour Very Short Form questionnaires .
Additional information is collected to describe the study population. Several variables will be assessed as confounders or effect modifiers of the relationship between IAP exposure and intestinal microbiota colonisation and secondary outcomes. Covariates that are collected at baseline include: maternal age, maternal pre-pregnancy weight, maternal and household smoking status, medications taken during pregnancy and first degree relative history of atopic disease, obesity, allergies or heart disease. Information about the pregnancy and birth are collected from birth records and charts including: GBS screening results, mode of delivery, place of birth (home or hospital), birth weight, length at birth, sex of the newborn, Apgar scores, maternal and infant antibiotic use prior to hospital discharge, intensive care unit admission and time from birth until discharge or the midwife leaving (homebirth). At follow up visits and through telephone contacts, parents are asked about time-varying covariates such as maternal and child exposure to antibiotics or antifungals, other medication or health product use, vaccine history, travel history, pets, daycare, exposure to tobacco smoke and changes in feeding (breast, formula or combination) including the addition of solid food.
Developing a sample size for studies in which the microbiota is the primary outcome is challenging. To date, most studies investigating the infant microbiome have included less than 100 subjects. To assess the feasibility of our study design, a sample size of 80 mother-infant pairs was selected. Recruitment and the collection of stool and data was found to be feasible, and attrition was low (6 %). Therefore we determined that enrolling 160 mother-infant pairs in Phase 2, for a total sample size of 240, was achievable.
We will use a diagram to summarize the patient flow in the study. Demographic and prognostic baseline characteristics will be reported as mean (standard deviation) for continuous variables and count (percent) for categorical variables. We will use regression analysis to analyze data for all clinical outcomes with intrapartum exposure to antibiotics as an independent variable after adjusting for confounders. Clinically important covariates will be investigated as potential confounders where appropriate for each outcome. The results will be expressed as odds ratio [OR] (for logistic regression for binary outcomes) or coefficient (for linear regression for continuous outcomes), corresponding standard error, 95 % confidence intervals and associated p-values. P-values will be reported to three decimal places with p-values less than 0.001 reported as p < 0.001. For all tests, we will use alpha = 0.05 level of significance. Assessment of model assumptions for regression analyses will be done by examining the residuals.
There is likely to be missing data that will likely increase with duration of the follow-up. We will use multiple imputation  to handle missing data. We will use generalized estimating equations (GEE)  to account for possible serial correlation of measurements within a participant overtime. Unlike ordinary regression analysis, GEE allows accounting for the possible correlation of outcomes for participants over time. Lastly, we will use propensity score methods to address the differential propensity for infant antibiotic exposure. All analyses will be performed using R or SPSS statistical software.
To understand the influence of IAP on the infant gut microbiota, we will analyze the microbial communities using culture-independent microbiota profiling methods. Microbiome analysis will include α-diversity metrics for each sample and β-diversity measures (Bray-Curtis) and other statistical analysis in R . Association of microbial community differences with sample groups will be assessed with permutational multivariate analysis of variance using Bray-Curtis dissimlarities (with vegan package in R) . Association of taxon abundance with sample groups will be assessed with a generalized linear mixed model with age and intrapartum exposure as fixed effects and individual child as a random effect using the lme4 package in R . Functional properties of the microbiota will be inferred using PICRUSt (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States) .
We enrolled 32 % of women who completed a ‘Consent to Contact’ form during the recruitment phase of our pilot study. The rates of follow-up and of data completion have been excellent. Of the first 83 participants recruited, one participant became ineligible prior to delivery (high risk pregnancy identified). Two participants have withdrawn; one prior to delivery and one at 12 weeks postpartum. Two participants were lost to follow up, resulting in 78 (94 %) participants remaining in the pilot cohort. In the pilot cohort, 94 % of follow up visits up to and including the 1-year visit were attended.
Strengths & limitations
We are developing a cohort of mother-infant pairs derived from the midwifery client population, providing us with the opportunity to examine a healthy pregnant and newborn population with low intervention rates and high breast-feeding rates. The homogeneity of this population may allow us to better examine relationships between the microbiome and exposures of interest. The midwifery population also gives us the opportunity for novel observations including place of birth (home and hospital) as a potential covariate. However, our selection of a low-risk cohort will preclude evaluation of nutritional intake that does not include breastfeeding or of exposures that occur with higher risk pregnancies or following preterm birth. Thus, our findings will need to be verified in these populations.
The investigative team is inter-professional in nature, has complementary areas of expertise, and includes national leaders in their respective fields. Our team provides expertise in pediatrics, obstetrics, midwifery, biostatistics, epidemiology, basic science including animal modeling, and microbiology, and is well suited both to undertake this research and to bring results back to the clinical practice arenas where the findings will be most relevant.
Our research will provide documentation of the development of the gut microbiome in a population of healthy full term infants born to low-risk women at full term gestation. By comparing the outcomes of exposed and unexposed infants we provide a first step in evaluating potential unforeseen negative consequences later in life of the current GBS clinical guidelines advocating IAP for GBS colonized women. Our team will explore this complex issue with examination of outcomes at both the microbiological level by examining intestinal colonisation in the first 3 years of life and in exploring the influence on systemic inflammation, adiposity, atopic disease and other outcomes of interest. Results from this study will make an important contribution to the growing body of knowledge and understanding of the patterns of intestinal microbiome establishment in the early years of life and the impact of IAP on the development of intestinal microbiota in newborns.
Dual-energy X-ray absorptiometry
Group B Streptococcus
Generalized estimating equations
Intrapartum antibiotic prophylaxis
Operational taxonomic units
Phylogenetic investigation of communities by reconstruction of unobserved states
Larsen N, Vogensen FK, van den Berg FW, Nielsen DS, Andreasen AS, Pedersen BK, Al-Soud WA, Sorensen SJ, Hansen LH, Jakobsen M. Gut microbiota in human adults with type 2 diabetes differs from non-diabetic adults. PLoS One. 2010;5:e9085.
Berer K, Mues M, Koutrolos M, Rasbi ZA, Boziki M, Johner C, Wekerle H, Krishnamoorthy G. Commensal microbiota and myelin autoantigen cooperate to trigger autoimmune demyelination. Nature. 2011;479:538–41.
Williams BL, Hornig M, Parekh T, Lipkin WI. Application of novel PCR-based methods for detection, quantitation, and phylogenetic characterization of Sutterella species in intestinal biopsy samples from children with autism and gastrointestinal disturbances. MBio. 2012;3. doi:10.1128/mBio.00261–11. Print 2012.
Mutlu EA, Gillevet PM, Rangwala H, Sikaroodi M, Naqvi A, Engen PA, Kwasny M, Lau CK, Keshavarzian A. Colonic microbiome is altered in alcoholism. Am J Physiol liver Physiol. 2012;302:G966–78.
David LA, Maurice CF, Carmody RN, Gootenberg DB, Button JE, Wolfe BE, Ling AV, Devlin AS, Varma Y, Fischbach MA, Biddinger SB, Dutton RJ, Turnbaugh PJ. Diet rapidly and reproducibly alters the human gut microbiome. Nature. 2014;505:559–63.
Stein RR, Bucci V, Toussaint NC, Buffie CG, Ratsch G, Pamer EG, Sander C, Xavier JB. Ecological modeling from time-series inference: insight into dynamics and stability of intestinal microbiota. PLoS Comput Biol. 2013;9:e1003388.
Costello EK, Stagaman K, Dethlefsen L, Bohannan BJ, Relman DA. The application of ecological theory toward an understanding of the human microbiome. Science. 2012;336:1255–62.
Caporaso JG, Lauber CL, Costello EK, Berg-Lyons D, Gonzalez A, Stombaugh J, Knights D, Gajer P, Ravel J, Fierer N, Gordon JI, Knight R. Moving pictures of the human microbiome. Genome Biol. 2011;12:R50. doi:10.1186/gb-2011-12-5-r50.
Gosalbes MJ, Llop S, Valles Y, Moya A, Ballester F, Francino MP. Meconium microbiota types dominated by lactic acid or enteric bacteria are differentially associated with maternal eczema and respiratory problems in infants. Clin Exp Allergy. 2013;43:198–211.
Tsuji H, Oozeer R, Matsuda K, Matsuki T, Ohta T, Nomoto K, Tanaka R, Kawashima M, Kawashima K, Nagata S, Yamashiro Y. Molecular monitoring of the development of intestinal microbiota in Japanese infants. Benef Microbes. 2012;3:113–25.
Matamoros S, Gras-Leguen C, Le Vacon F, Potel G, de La Cochetiere MF. Development of intestinal microbiota in infants and its impact on health. Trends Microbiol. 2013;21:167–73.
Yatsunenko T, Rey FE, Manary MJ, Trehan I, Dominguez-Bello MG, Contreras M, Magris M, Hidalgo G, Baldassano RN, Anokhin AP, Heath AC, Warner B, Reeder J, Kuczynski J, Caporaso JG, Lozupone CA, Lauber C, Clemente JC, Knights D, Knight R, Gordon JI. Human gut microbiome viewed across age and geography. Nature. 2012;486:222–7.
Koenig JE, Spor A, Scalfone N, Fricker AD, Stombaugh J, Knight R, Angenent LT, Ley RE. Succession of microbial consortia in the developing infant gut microbiome. Proc Natl Acad Sci U S A. 2011;108 Suppl 1:4578–85.
Bäckhed F, Roswall J, Peng Y, Feng Q, Jia H, Kovatcheva-Datchary P, Li Y, Xia Y, Xie H, Zhong H, Khan MT, Zhang J, Li J, Xiao L, Al-Aama J, Zhang D, Lee YS, Kotowska D, Colding C, Tremaroli V, Yin Y, Bergman S, Xu X, Madsen L, Kristiansen K, Dahlgren J, Jun W. Dynamics and Stabilization of the Human Gut Microbiome during the First Year of Life. Cell Host Microbe. 2015;17:690–703.
Fouhy F, Guinane CM, Hussey S, Wall R, Ryan CA, Dempsey EM, Murphy B, Ross RP, Fitzgerald GF, Stanton C, Cotter PD. High-throughput sequencing reveals the incomplete, short-term recovery of infant gut microbiota following parenteral antibiotic treatment with ampicillin and gentamicin. Antimicrob Agents Chemother. 2012;56:5811–20.
Penders J, Thijs C, Vink C, Stelma FF, Snijders B, Kummeling I, van den Brandt PA, Stobberingh EE. Factors influencing the composition of the intestinal microbiota in early infancy. Pediatrics. 2006;118:511–21.
Jakobsson HE, Abrahamsson TR, Jenmalm MC, Harris K, Quince C, Jernberg C, Bjorksten B, Engstrand L, Andersson AF. Decreased gut microbiota diversity, delayed Bacteroidetes colonisation and reduced Th1 responses in infants delivered by Caesarean section. Gut. 2013;63:559–66.
Fujimura KE, Sitarik AR, Havstad S, Lin DL, Levan S, Fadrosh D, Panzer AR, LaMere B, Rackaityte E, Lukacs NW, Wegienka G, Boushey HA, Ownby DR, Zoratti EM, Levin AM, Johnson CC, Lynch SV. Neonatal gut microbiota associates with childhood multisensitized atopy and T cell differentiation. Nature Medicine. 2016;22(10):1187–1191.
Black M, Bhattacharya S, Philip S, Norman JE, McLernon DJ. Planned Cesarean Delivery at Term and Adverse Outcomes in Childhood Health. JAMA. 2015;314:2271–9.
Cho CE, Norman M. Cesarean section and development of the immune system in the offspring. Am J Obstet Gynecol. 2013;208:249–54.
Money D, Allen VM, Canada S of O and G of. The prevention of early-onset neonatal group B streptococcal disease. J Obstet Gynaecol Can. 2013;35:939–51.
Sibley CD, Grinwis ME, Field TR, Eshaghurshan CS, Faria MM, Dowd SE, Parkins MD, Rabin HR, Surette MG. Culture enriched molecular profiling of the cystic fibrosis airway microbiome. PLoS One. 2011;6:e22702.
Sibley CD, Parkins MD, Rabin HR, Duan K, Norgaard JC, Surette MG. A polymicrobial perspective of pulmonary infections exposes an enigmatic pathogen in cystic fibrosis patients. Proc Natl Acad Sci U S A. 2008;105:15070–5.
Whelan FJ, Verschoor CP, Stearns JC, Rossi L, Luinstra K, Loeb M, Smieja M, Johnstone J, Surette MG, Bowdish DME. The loss of topography in the microbial communities of the upper respiratory tract in the elderly. Ann Am Thorac Soc. 2014;11:513–21.
Bartram AK, Lynch MD, Stearns JC, Moreno-Hagelsieb G, Neufeld JD. Generation of multimillion-sequence 16S rRNA gene libraries from complex microbial communities by assembling paired-end illumina reads. Appl Environ Microbiol. 2011;77:3846–52.
McDonald D, Price MN, Goodrich J, Nawrocki EP, DeSantis TZ, Probst A, Andersen GL, Knight R, Hugenholtz P. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012;6:610–8.
DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, Huber T, Dalevi D, Hu P, Andersen GL. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72:5069–72.
Muraro A, Dreborg S, Halken S, Host A, Niggemann B, Aalberse R, Arshad SH, von Berg A, Carlsen KH, Duschen K, Eigenmann P, Hill D, Jones C, Mellon M, Oldeus G, Oranje A, Pascual C, Prescott S, Sampson H, Svartengren M, Vandenplas Y, Wahn U, Warner JA, Warner JO, Wickman M, Zeiger RS. Dietary prevention of allergic diseases in infants and small children. Part II. Evaluation of methods in allergy prevention studies and sensitization markers. Definitions and diagnostic criteria of allergic diseases. Pediatr Allergy Immunol. 2004;15:196–205.
Asher MI, Keil U, Anderson HR, Beasley R, Crane J, Martinez F, Mitchell EA, Pearce N, Sibbald B, Stewart AW. International Study of Asthma and Allergies in Childhood (ISAAC): rationale and methods. Eur Respir J Off J Eur Soc Clin Respir Physiol. 1995;8:483–91.
Sadeh A. A brief screening questionnaire for infant sleep problems: validation and findings for an Internet sample. Pediatrics. 2004;113:e570–7.
Kleinman L, Rothman M, Strauss R, Orenstein SR, Nelson S, Vandenplas Y, Cucchiara S, Revicki DA. The infant gastroesophageal reflux questionnaire revised: development and validation as an evaluative instrument. Clin Gastroenterol Hepatol. 2006;4:588–96.
Varni JW, Bendo CB, Denham J, Shulman RJ, Self MM, Neigut DA, Nurko S, Patel AS, Franciosi JP, Saps M, Verga B, Smith A, Yeckes A, Heinz N, Langseder A, Saeed S, Zacur GM, Pohl JF. PedsQL gastrointestinal symptoms module: feasibility, reliability, and validity. J Pediatr Gastroenterol Nutr. 2014;59:347–55.
Croft S, Stride C, Maughan B, Rowe R. Validity of the strengths and difficulties questionnaire in preschool-aged children. Pediatrics. 2015;135:e1210–9.
de la Osa N, Granero R, Penelo E, Domènech JM, Ezpeleta L. The short and very short forms of the Children’s Behavior Questionnaire in a community sample of preschoolers. Assessment. 2014;21:463–76.
Little RJA, Rubin DB. Statistical analysis with missing data. New York: John Wiley & Sons; 1987.
Everitt BS, Howell DC. Encyclopedia of Statistics in Behavioral Science. Chichester: John Wiley & Sons, Ltd; 2005.
R Development Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2011. http://www.r-project.org/. ISBN 3-900051-07-0.
Oksanen J, Guillaume Blanchet F, Kindt R, Legendre P, Minchin P, O’hara R, Simpson G, Solymos P, Stevens M, Wagner H. Vegan: Community Ecology Package, R Package Version 2.3-0. 2015.
Bates D, Maechler M, Bolker B, Walker S. _lme4: Linear Mixed-Effects Models Using Eigen and S4_, R Package Version 1.1–8. 2015.
Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, Clemente JC, Burkepile DE, Vega Thurber RL, Knight R, Beiko RG, Huttenhower C. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31:814–21.
We would like to thank the women and their babies who are participating in this study for their valuable time. We gratefully acknowledge the midwives and administrators at The Hamilton Midwives, Community Midwives of Hamilton, Access Midwives, Burlington and Area Midwives, Community Midwives of Brantford, Lincoln Community Midwives, Niagara Midwives and Renaissance Midwifery for the time they committed to recruiting participants and collecting data. We would like to thank Laura Rossi and Michelle Shah for their work in processing the stool samples for analysis. We would like to acknowledge the assistance that has been provided by Sarah Kanji, Glara Rhee, Kimberly Farrell and Samantha Krueger in developing study protocols and conducting study visits.
This study is funded by grants from the Hamilton Academic Health Sciences Organization and the Canadian Institutes of Health Research (MOP-136811). JDS holds CDA Scholar (SC-5-12-3891-JS) CIHR New Investigator awards (MSH-136665).
Availability of data and materials
Study concept and design: EKH, ACH, MGS, KMM. Drafting of the manuscript: JS, EKH, EG, JCS, KMM. Critical revision of the manuscript for important intellectual content: ACH, HM, AM, JDS, EMR, LT. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The study was approved by the Research Ethics boards at Hamilton Health Sciences (12–201), St. Joseph’s Healthcare Hamilton (12–3721), Joseph Brant Hospital (000–022–14), Niagara Health System (2014–12–001) and Brant Community Healthcare System. Written informed consent was provided by all participants prior to enrollment in the study.
Case Report Forms. This file contains the case report forms to be completed by participants and research staff from the time of eligibility assessment and enrollment to the time of the 3 year study visit. (PDF 1007 kb)