Article Text

Effect of antibiotics in the first week of life on faecal microbiota development
  1. Emmy Van Daele1,
  2. Kim Kamphorst2,3,
  3. Arine M Vlieger3,
  4. Gerben Hermes1,
  5. Christian Milani4,5,
  6. Marco Ventura4,5,
  7. Clara Belzer1,
  8. Hauke Smidt1,
  9. Ruurd M van Elburg6,
  10. Jan Knol1,7
  1. 1 Laboratory of Microbiology, Wageningen University & Research, Wageningen, The Netherlands
  2. 2 Pediatrics, Amsterdam Gastroenterology, Metabolism & Nutrition, Amsterdam Reproduction & Development, Amsterdam UMC Locatie AMC, Amsterdam, The Netherlands
  3. 3 Pediatrics, St. Antonius Hospital, Nieuwegein, The Netherlands
  4. 4 Laboratory of Probiogenomics, Department of Chemistry, Life Sciences and Environmental Sustainability, University of Parma Department of Chemical Life Sciences and Environmental Sustainability, Parma, Emilia-Romagna, Italy
  5. 5 Interdepartmental Research Centre "Microbiome Research Hub", University of Parma, Parma, Emilia-Romagna, Italy
  6. 6 Pediatrics, Amsterdam Gastroenterology, Metabolism & Nutrition, Amsterdam Reproduction & Development Amsterdam, Amsterdam UMC Locatie AMC, Amsterdam, North Holland, The Netherlands
  7. 7 Nutricia Research BV, Utrecht, The Netherlands
  1. Correspondence to Professor Jan Knol, Nutricia Research BV, 3584 CT Utrecht, Utrecht, The Netherlands; jan.knol{at}danone.com

Abstract

Background Infants are frequently exposed to antibiotics (AB) in the first week of life for suspected bacterial infections. Little is known about the effect of AB on the developing intestinal microbiota. Therefore, we studied intestinal microbiota development with and without AB exposure in the first week of life in term born infants.

Methods We analysed the faecal microbiota from birth until 2.5 years of age by 16S rRNA gene amplicon sequencing in a cohort with 56 term born infants, exposed to AB in the first week of life (AB+) (AB for 2–3 days (AB2, n=20), AB for 7 days (AB7, n=36)), compared with 126 healthy controls (AB-). The effects of AB and duration were examined in relation to delivery and feeding mode.

Results AB+ was associated with significantly increased relative abundance of Enterobacteriaceae at 3 weeks and 1 year and a decrease of Bifidobacteriaceae, from 1 week until 3 months of age only in vaginally delivered, but not in C-section born infants. Similar deviations were noted in AB7, but not in AB2. After AB, breastfed infants had lower relative abundance of potentially pathogenic Enterobacteriaceae compared with formula fed infants and recovered 2 weeks faster towards controls.

Conclusions AB exposure in the first week of life alters faecal microbiota development with deviations in the relative abundance of individual taxa until 1 year of age. These alterations can have long-term health consequences, which emphasises the need for future studies aiming at restoring intestinal microbiota after AB administration.

  • Microbiology
  • Infant Development
  • Gastroenterology
  • Epidemiology
  • Neonatology

Data availability statement

Data are available on reasonable request. The data are available on request with the corresponding author.

http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

What is already known on this topic

  • Up to 20% of neonates receive antibiotics because of (suspected) early-onset neonatal sepsis. In older infants, antibiotics have been shown to disturb the intestinal microbiota, but studies in newborns with a developing microbiota are limited.

  • Feeding type and delivery mode also affect the microbiota, but their effect in relation to antibiotic exposure in the first week of life has not yet been studied.

What this study adds

  • Antibiotic exposure in the first week of life affects microbiota development throughout the first year of life, with more profound deviations in infants born at term exposed for 7 days compared with only 2–3 days.

  • Antibiotic exposure in the first week of life resulted in deviations in the faecal microbiota development of vaginally delivered infants but not of C-section delivered infants compared with their respective controls. This could be attributed to the difference caused by delivery mode itself, with C-section delivered infants deviating from vaginal controls up to 1 month of age regardless of antibiotic exposure.

How this study might affect research, practice and/or policy

  • This study underlines the importance of early cessation of antibiotics started at birth because of the prolonged effect on the intestinal microbiota and possible impact on health.

Introduction

During the first 1000 days of life, the intestinal microbiota impacts health in later life through the interdependent development of microbiota, immune system, growth and cognitive function.1 2 After birth, the intestinal microbiota develops rapidly, driven by exposure to microbes from maternal, environmental and dietary sources.3 During this early development, the intestinal microbiota is unstable and susceptible to perturbations such as those caused by antibiotic (AB) exposure. These perturbations may have long-term consequences on the developing microbiota and also on the developing immune system,4 growth5 6 and have already been associated with increased prevalence of asthma, allergies, coeliac disease, eczema, eosinophilic esophagitis, infantile colic, inflammatory bowel disease and obesity.7–11

In mice-based studies, AB exposure in the first week of life altered microbiota composition and immune function, but gavage with mature untreated microbiota, restored the perturbation and reduced the negative health effects.12 To understand the full impact in humans and develop restoration strategies, more knowledge is needed.

Worldwide, up to 20% of all neonates are prescribed ABs because of (suspected) early-onset neonatal sepsis, although in most cases, sepsis is unconfirmed and ABs could be discontinued after 48–72 hours.13–15 More prolonged AB exposure can gradually reduce overall diversity16 and richness17 of the neonate’s intestinal ecosystem. Additionally, AB type also determines the microbial perturbation through specific mechanisms of action and host interactions. Studies in infants have been inconclusive,17–19 but indicated types with faster recovery towards the microbiota composition of controls.20 Therefore, more comparative studies are needed between durations and types in AB regimes to optimise AB administration.21

In this study, we investigated the microbiota development in a subset of the Intestinal Microbiota Composition after AB treatment in early life (INCA) cohort.22 The primary aim was to investigate the impact of AB exposure in the first days of life on microbiota development during the first 2.5 years of life. Secondary aims were to examine (1) short (2–3 days) versus long (7 days) AB administration, (2) different AB types and (3) the impact of feeding and delivery mode on AB perturbation.

Materials and methods

Study design

This prospective, observational study has been described previously.10 22 To study the impact of feeding mode on AB-induced deviations, we selected a subset of 182 infants with 1128 samples who were exclusively breast-fed (BF) or formula-fed (FF) in the first 3 months of life. An overview of collected samples and their selection for analyses reported here is provided in online supplemental figure 1). All AB+ infants received gentamicin, which was combined with penicillin (ABPen), amoxicillin (ABAMX) or amoxicillin with clavulanic acid (ABAMC).

Data and faeces collection

Baseline characteristics such as birth mode were assessed through hospital records. Feeding mode was reported monthly during the first year of life. Nine faecal samples were collected from infants (figure 1). Until discharge from the hospital, faeces were sampled from diapers by hospital staff and immediately frozen at −20°C. Sampling continued at home by the parents, using sampling instructions and freezer storage. After 1 year, parents delivered the samples to the clinic after transport on ice. A final sample was taken around 2 years, stored in the home freezer and collected by the study nurse. At the hospital, all samples were stored at −20°C.

Figure 1

Overview of the collected questionnaires, faecal sampling points and the age categories based on the age (days) at sampling. N, number of infants for which samples were available for a given time point. d(ays), m(onths) and y(ears) indicate the age category ranging from birth until around 2 years of age.

16S rRNA gene amplicon sequencing

DNA was extracted from a maximum of 200 mg of the homogenised faecal sample at GenProbio srl (Parma, Italy) with the QIAamp DNA Stool Mini kit according to manufacturer’s instructions (Qiagen Ltd, Strasse, Germany). Sequencing libraries were prepared according to the 16S Metagenomic Sequencing Library Preparation Protocol (Part No. 15 044 223 Rev. B—Illumina) at GenProbio srl (Parma, Italy). Minor adaptations to the published protocol23 are noted in online supplemental file 1). Sequencing resulted in ~44 934 (SD 17205) reads per sample. Data were processed using NG-Tax 2.0 on demultiplexed fastq files, using default settings.24 Taxonomy was assigned using SILVA reference database V.128.25 Amplicon sequence variants (ASVs) were defined as unique sequence variants. Two synthetic mock communities were sequenced as positive controls.26

Statistical analysis of the baseline characteristics

Baseline characteristics were calculated using R 3.6.127 and tableone28 (table 1). Differences between groups were examined by using the Fisher exact test for the categorical variables and analysis of variance (ANOVA) for continuous variables. Not normally distributed variables were tested with the Kruskal-Wallis test and indicated by their median and IQR. Bonferroni correction was performed to correct for multiple testing.

Table 1

Baseline characteristics of the INCA cohort subset included in this study

Bioinformatic and statistical analysis of the sequence data

All analyses were performed in R 3.6.1.27 Samples were stratified into 11 age-based and right-closed intervals for statistical analysis (figure 1). The Jenks Natural Breaks Classification (classInt package) was used to calculate the optimal ranges.29 Because of increasing increments between sampling, the x-axis (age) was log2 transformed for visualisation.

Alpha diversity (within sample diversity) was calculated at ASV level using Picante30 and Microbiome31 packages and following metrics: Faith’s phylogenetic diversity,32 ASV richness, Shannon and Inverse Simpson. All except Shannon diversity needed logarithmic transformation to obtain normal distribution for one-way analysis of covariance. Consecutive analyses were corrected for the baseline characteristic that differed significantly between AB+ and AB- and between AB2 and AB7.

Temporal trends in relative abundance were visualised using local regression with locally estimated scatterplot smoothing using ggplot2. These relative abundances did not meet normality requirements and were therefore compared using beta regression with BetaReg33 per age interval. The effects of AB on the relative abundance of phyla were modelled including and excluding delivery mode, feeding mode and additional AB exposure between one and 6 months. The optimal beta regression models, based on Akaike and Bayesian information criteria, only included AB exposure in the first week of life without additional terms.

Beta diversity (between sample diversity) was calculated using pairwise Weighted (WU)34 and Unweighted UniFrac (UU) distances.35 WU takes the relative abundance of each ASV into account, whereas UU uses presence or absence of ASVs. Pairwise UU and UW distance matrices were used to plot principle response curves (PRCs).36 PRC analysis is a redundancy analysis for interpreting longitudinal data. It visualises multivariate responses in a repeated observation design.36 37 The method was designed to analyse the treatment effect over time compared with controls, as it can disentangle the time-dependent effects from other possible determinants.38 39 In this study, time was displayed on the x-axis and the intestinal microbiota development was shown compared with AB- infants as the baseline reference group. Differences between AB groups were assessed per age interval using ANOVA in the vegan package.40

Results

Baseline characteristics of INCA cohort subset

The baseline characteristics differed between AB- and AB+ for gestational age, birth weight and additional AB exposure between one to 6 months (table 1). Birth weight z-score (birth weight corrected for gestational age), however, was comparable. AB2 differed from AB7 with regard to additional AB exposure between 1 and 6 months (5% vs 34%, respectively, p=0.019). Baseline characteristics were comparable between the AB type groups (online supplemental table S1).

Antibiotic-induced alterations to intestinal microbiota development

AB exposure during the first week of life did not alter microbial alpha diversity between birth and 2.5 years (online supplemental figure 1). The temporal patterns of the four major phyla (95% of the average relative abundance) were compared with univariate analyses (figure 2). Relative abundance of Proteobacteria was high overall, during the first months of life. AB exposure further increased this relative abundance at 1 month (mean 43.6% AB+, 31.5% AB-) and 1 year (mean 17.6% AB+, 6.5% AB-). Actinobacteria peaked around 3 months, but AB exposure decreased their relative abundance at 1 week (mean 6.3% AB+, 18.2% AB-), 2 weeks (mean 8.4% AB+, 24.4% AB-), 1 month (mean 17.5% AB+, 27.2% AB-) and 3 months (mean 26.5% AB+, 34.4% AB-). The average relative abundance of Bacteroidetes was stable at approximately 10% during the first year of life. Firmicutes drastically increased in relative abundance towards 2.5 years. Both were unaffected by AB.

Figure 2

Temporal trajectories of the relative abundance of the four main phyla in the developing intestinal microbiota from birth to 2.5 years of age. Thick lines represent the average with shading showing the 95% CI. Differences between AB+ and AB- were calculated using beta regression for each age category. * and vertical grey shading: difference in relative abundance (p<0.05), AB+, infants who received AB during their first week of life, AB-, infants who did not receive AB during their first week of life, LOESS, locally estimated scatterplot smoothing.

Effect of antibiotic duration on long-term microbiota development

Based on univariate statistics, the temporal trajectories of Actinobacteria and Proteobacteria were most affected by AB. Based on UU, ABs increased several members of the Enterobacteriaceae and decreased Bifidobacteriaceae at one and 2 weeks (online supplemental figure 2A). For WU, ABs impacted similar taxa at 1 year (p<0.05) (figure 3A). These AB deviations were similar in AB2 and AB7, but only the microbiota composition of AB7 differed from AB- baseline (figure 3B and online supplemental figure 2B).

Figure 3

WU-based PRC analysis is a special case of RDA for multivariate responses in a repeated observation design, which is applied here on the longitudinal infant intestinal microbiota dataset from birth till 2.5 years of age. (A) The infants who did not receive antibiotics during the first week of life (AB-), were compared as a baseline to the antibiotic exposed infants (AB+). Bacterial genera shown are the main drivers of the differences between AB+ and AB-: taxa on the same side of baseline as the curve are linked to an increased relative abundance at that time point, opposite sides indicate a decrease. (B) AB- was also compared as a baseline with the different antibiotic durations of 2–3 (AB2) or 7 days (AB7). Significance was tested at the different time points using an ANOVA like permutation test (*p<0.05 compared with baseline AB-). Covariates that were controlled for included additional AB exposure between the age of one and 6 months. AB+, infants who received AB during their first week of life; AB-, infants who did not receive AB during their first week of life; AB2, AB exposure for 2–3 days in the first week of life indicated with light grey shading; AB7, AB exposure for 7 days indicated with grey shading; ANOVA, analysis of variance; ASV, amplicon sequence variants; PRC, principal response curve; RDA, redundancy analysis; WU, weighted UniFrac.

AB types had a different impact on the microbiota (online supplemental figure 3A,B) with ABAMX not deviating from AB- baseline and ABPEN deviating at 1 week and 1 year with increased relative abundance of Enterobacteriaceae members in WU and UU. UU-based deviations between ABAMC and AB- baseline were limited to week one and involved different ASVs compared with ABPEN. ABAMC affected WU in the long-term with increased Enterobacteriaceae and Enterococcaceae at 2 weeks and 1 month and also Bifidobacterium at 1 week and 2.5 years.

Impact of delivery and feeding mode on AB-associated deviations in the faecal microbiota development

Due to the relatively low number of faecal samples in the first week per feeding and delivery type (figure 1), effects were only reported in samples collected between 1 week and 2.5 years. Within AB-, microbiota deviated based on delivery mode from 1 week until 1 month (figure 4D). In vaginally delivered infants, the AB effect on the microbiota was still significant at 1 year with an increase of several Enterobacteriaceae, Enterococcaceae and Streptococcaceae and decrease in Bifidobacterium and Escherichia-Shigella. In contrast, no AB-mediated deviations were noted between C-section born infants.

Figure 4

WU-based PRC analyses in the different delivery and feeding mode groups. Bacterial genera shown are the main drivers for differences between the groups and the baseline: positive effect on the curve is linked to increased ASVs in the positive spectrum and decrease of those in the negative spectrum. (A) Breastfed controls (AB-BF) were compared as a baseline to antibiotic exposed (AB+) and FF infants. (B) Specifically featuring the AB effect in breastfed children, with breastfed control children as a baseline (AB-BF) and (C) formula-fed children, with formula-fed control children as a baseline (AB-). (D) Vaginally delivered control infants (AB-Vag) were compared as a baseline to antibiotic exposed (AB+) and C-section delivered (Sectio) infants. Significance was tested at the different time points using an ANOVA like permutation test (*p<0.05). Delivery mode analyses were controlled for feeding mode and vice versa. AB-, infants who did not receive AB during their first week of life; AB+, infants who received AB during their first week of life also indicated within grey shading; ASV, amplicon sequence variants; BF, infants exclusively breast-fed in the first 3 months of life; FF, infants exclusively formula-fed in the first 3 months of life; PRC, principal response curve; Sectio, infants delivered through C-section, Vag, vaginally delivered infants; WU, weighted UniFrac.

Compared with AB- BF baseline, AB+BF infants only deviated at 2 weeks, whereas AB+FF infants showed longer deviations from the first week up to 1 month (figure 4A). Because there was also a feeding effect during the first 6 months, the AB effect was also analysed within the separate feeding groups. AB+ was associated with decreased Bifidobacterium relative abundance in FF (1 month) (figure 4C), which occurred later than in BF infants (2 weeks) (figure 4B). In turn, AB+ was associated with increased relative abundance of Parabacteroides in BF, and with increased relative abundance of Enterobacteriaceae and Enterococcus in FF infants.

Discussion

In this prospective, observational INCA study, we examined the microbiota development after AB exposure during the first week of life and found perturbations in the faecal microbiota development from 1 week until 1 year of age. These perturbations included decreased relative abundance of Bifidobacteriaceae while potentially pathogenic Enterobacteriaceae increased. This study adds new insights into long-term compositional shifts after neonatal AB exposure.16 41 Our results corroborate findings in older infants with increased Enterobacteriaceae and decreased Bifidobacteriaceae after AB administration.16 20 42–44 Importantly, the severity and duration of AB-mediated microbiota perturbations increased with longer AB administration (5–7 vs 2–3 days). The results also align with a small study in preterm infants, where >5 days AB exposure intensified perturbations compared with 2–3 days.45

Bifidobacteriaceae form a cornerstone in the early development of the immune system. They were shown to promote B-cell maturation and associations with decreased inflammatory responses and T-regulatory cell acquisition.46 Enterobacteriaceae, on the other hand, produce toxins and have lipopolysaccharides on their outer membranes, which causes inflammation.47–49 Therefore, it is not surprising that reduced Bifidobacteriaceae, often combined with an increase in potentially pathogenic Enterobacteriaceae, like Shigella, Klebsiella and Enterobacter, have been associated with immune-mediated disorders like asthma.50 Similar deviations were also found in functional disorders like infantile colic51 and irritable bowel syndrome.52

The long-term microbiota effect of ABs and its associated negative health outcomes reinforce the need for implementing AB stewardship programs53–55 to avoid AB overuse.21 56 57 The microbiota perturbations were only significant after 5–7 days compared with 2–3 days AB which could explain previous findings from the INCA cohort: namely higher incidence of infantile colic, wheezing and food allergies in infants exposed for 5–7 days, but not for 2–3 days.10 58 If this observed difference between 2–3 days and 5–7 days exposure is the result of longer AB exposure or the result of a concomitant infection or inflammatory response is yet unclear. The AB7 infants were treated because of suspected early onset sepsis (EOS). EOS is rare in term infants,59–61 but it is difficult to distinguish from normal neonatal physiology after birth, and laboratory tests cannot always reliably detect or rule out EOS.62 Because the consequences of delaying treatment are significant, on average 82 newborns without EOS are treated for each case.15 63–65 In our study population, only two of the 36 AB7 infants had a positive blood culture. The others were also treated for 5–7 days because of elevated inflammatory markers or clinical symptoms. Uzan-Yulzari et al showed that the association between neonatal AB exposure and growth was independent of the neonatal infection state.6 This suggests that the differences in microbiota development after AB treatment in our study are more likely the result of the AB treatment duration itself than caused by a possible EOS. Our new findings emphasise the need for microbiota restoration to minimise aberrant immune development. Suggested strategies include prebiotics, probiotics and synbiotics66 but also faecal transfers, which partially restored the microbiota of mice exposed to AB for 7 days.12

In vaginally delivered infants, the AB effect was most pronounced with microbiota deviations in the second week of life. C-section born infants, however, showed similar perturbations regardless of AB exposure. After C-section, microbiota perturbations occurred due to reduced vertical mother-infant transmission of important intestinal microbes such as Bacteroides and Bifidobacterium, while transmission of other microbes like skin and mouth bacteria increased,67 as well as due to maternal AB administration prior to cord clamping.3 C-section delivery already showed decreases in Bifidobacterium spp and increases in opportunistic pathogens from hospital environments like Enterobacter, Enterococcus and Klebsiella.68 This resembled the AB effect, which might explain the lack of an additional AB effect in C-section infants as these infants already lack the affected microbial groups from birth.

Feeding also has a major impact on early life microbiota development.69 In our study, ABs in the first week of life perturbed the microbiota of both BF and FF infants, but potentially pathogenic Enterobacteriaceae only increased in FF infants. Moreover, AB perturbations were still notable at 1 month in FF infants but only until 2 weeks in BF infants. Breastmilk probably aids restoration through components like human milk oligosaccharides and live bacteria, which stimulate the growth of bifidobacteria and reduce (potential) pathogens.70

The strengths of this study are the quantity of samples and long-term follow-up. This enabled the investigation of the interplay of AB with delivery and feeding modes. Moreover, the quantity of sampling points allowed for detailed and long-term detection of AB-induced perturbations within individuals. This was relevant as AB impact was not uniform over time, suggesting that limited sampling points could lead to misinterpretation. Finally, the number of infants receiving additional courses of AB in the first year was low in this cohort, thereby reducing an important confounding factor.

The methodology for profiling the intestinal microbiota, which targeted the V3 hypervariable regions of the bacterial 16S rRNA gene, provides a cost-effective overview of bacterial community composition, however, resolution at species level is limited, and amplification bias cannot be unequivocally ruled out.71 72 The applied Probio_Uni and Probio_Rev primers were validated and, compared with other primers, they seemed to represent relevant members of the intestinal microbiota such as bifidobacteria more accurately, which makes them especially fit for analysing the intestinal microbiota of infants.23 For future research, whole genome shotgun sequencing could be used to increase the accuracy of species and strain detection.73 74 Another limitation may be that environmental factors and maternal-infant interactions could have been confounders, because AB+infants were admitted to neonatal wards, whereas AB- infants stayed with their mothers on the maternity ward and were discharged earlier. Last, we did not have sufficient data on perinatal AB exposure and were therefore unable to correct for it, although it is questionable to what extent this confounder is important to take into account.75 76 Additionally, the study was not primarily designed (and thus underpowered) to conclude on AB types. Nevertheless, our results suggest that ABAMX induced less perturbations as it did not result in any differences from AB- (online supplemental figure 3B). The addition of the β-lactamase inhibitor clavulanic acid (ABAMC) was associated with higher levels of bifidobacteria compared with other AB types, which supports an earlier finding in a single subject.20 Dedicated studies are, however, needed to further elucidate the optimal regime with the least microbial perturbations.

In conclusion, AB exposure in the first week of life in term born infants disturbed the microbiota up to 1 year, with more significant deviations after longer AB exposure (5–7 days). Both C-section delivery and AB administration in the first week of life are associated with deviant intestinal microbiota, but the two combined are not associated with further deviation. Breastfeeding was associated with reduced severity and duration of perturbations compared with formula feeding. Our observations may help to elucidate why AB-exposed infants have more health problems. It may also support the development of preventive and curative strategies after AB exposure to stimulate the growth of beneficial microbiota in order to prevent future health problems.

Data availability statement

Data are available on reasonable request. The data are available on request with the corresponding author.

Ethics statements

Patient consent for publication

Ethics approval

The INCA study was approved by the ethical board of the St. Antonius Hospital in Nieuwegein, the Netherlands (registered as NCT02536560), and informed consent was obtained from both parents. Participants gave informed consent to participate in the study before taking part.

Acknowledgments

First of all, we want to thank the mothers and infants who participated in the INCA cohort for their time and effort. Second, we would like to acknowledge the help of all the hospital staff in the recruitment and collection of data, and especially Carin Bunkers and Nicole Rutten for their support. Last, we thank Marta Magnifesta at GenProbio srl (Parma, Italy) for performing the DNA extraction and sequencing of the faecal samples.

References

Supplementary materials

  • Supplementary Data

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

Footnotes

  • RMvE and JK contributed equally.

  • Contributors JK and RMvE contributed equally to this paper. EVD conceptualized the research questions, wrote the manuscript, performed the analyses and data interpretation, wrote and revised the manuscript. KK provided data and revised the manuscript. AMV conceptualised the original cohort study and provided a clinical outlook for data interpretation and writing and revised the manuscript. GH conceptualized the research questions, provided support for the data analysis and revised the manuscript. CM performed data analysis and revised the manuscript. MV performed data analysis and revised the manuscript. CB acquired funding and revised the manuscript. HS supervised the research activity and revised the manuscript. RMvE conceptualised the original cohort study and provided a clinical outlook for data interpretation and writing and revised the manuscript. JK is guarantor, supervised the research activity and revised the manuscript.

  • Funding This work was supported by JPI HDHL in conjunction with ZonMW and Danone Nutricia Research and grant number IM2015.

  • Competing interests This work and the PhD research by EVD was financed by a EU Joint Programming Initiative namely A Healthy Diet for a Healthy Life (JPIHDHL, http://www.healthydietforhealthylife.eu/) inconjunction with ZonMW and Danone Nutricia Research. GH is a full full-time employee of Chr Hansen A/S since January 2020. CB received a grant financed by JPI HDHL inconjunction with ZonMW and Danone Nutricia Research. RMvE was an employee at Danone Nutricia Research till 2020. JK is a full-time employee of Danone Nutricia.

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

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

Request Permissions

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