# Very Preterm Children Gut Microbiota Comparison at the Neonatal Period of 1 Month and 3.5 Years of Life

Gaël Toubon, Marie-José Butel, Jean-Christophe Rozé, Patricia Lepage, Johanne Delannoy, Pierre-Yves Ancel, Marie-Aline Charles, Julio Aires

► To cite this version:

Gaël Toubon, Marie-José Butel, Jean-Christophe Rozé, Patricia Lepage, Johanne Delannoy, et al. Very Preterm Children Gut Microbiota Comparison at the Neonatal Period of 1 Month and 3.5 Years of Life. Frontiers in Microbiology, 2022, 13, 10.3389/fmicb.2022.919317. [https://doi.org/10.3389/fmicb.2022.919317](https://doi.org/10.3389/fmicb.2022.919317). This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International License.

HAL Id: hal-03901708

https://hal.inrae.fr/hal-03901708

Submitted on 15 Dec 2022

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Very Preterm Children Gut Microbiota Comparison at the Neonatal Period of 1 Month and 3.5 Years of Life

Gaël Toubon<sup>1,2,3</sup>, Marie-José Butel<sup>2,3</sup>, Jean-Christophe Rozé<sup>4</sup>, Patricia Lepage<sup>5</sup>, Johanne Delannoy<sup>2,3</sup>, Pierre-Yves Ancel<sup>1,3</sup>, Marie-Aline Charles<sup>1</sup> and Julio Aires<sup>2,3*</sup> for the EPIFLORE Study Group

Prematurity is a risk factor for dysbiosis of the gut microbiota due to particular birth conditions and frequent prolonged hospitalization of neonates. Although gut microbiota colonization after birth and its establishment during the hospitalization period have been studied in preterm infants, data on gut microbiota following discharge, particularly during early childhood, are scarce. The present study investigated the relationship between gut microbiota at 1 month after birth (hospitalization period) and 3.5 years of age in 159 preterm children belonging to the French EPIFLORE prospective observational cohort study. Analysis using bacterial 16S rRNA gene sequencing showed that the gut microbiota of preterm neonates at 1 month was highly variable and characterized by six distinct enterotypes. In contrast, the gut microbiota of the same children at 3.5 years of age showed less variability, with only two discrete enterotypes. An absence of association between enterotypes at 1 month and 3.5 years of age was observed. While the alpha diversity of gut microbiota significantly increased between 1 month and 3.5 years of age, for both alpha and beta diversities, there was no correlation between the 1-month and 3.5-years time points. Comparison at 3.5 years between children born either preterm (n = 159) or full-term (n = 200) showed no differences in terms of enterotypes, but preterm children harbored a lower Shannon diversity index and a different overall composition of microbiota than full-term children. This study suggests that the characteristics of the early gut microbiota of preterm children are not predictive of the microbial community composition at 3.5 years of age. However, the impact of gestational age is still noticeable on the gut microbiota up to 3.5 years of age.

**Keywords:** prematurity, gut microbiota, DOHaD, children, enterotypes

## INTRODUCTION

The human microbiota has become a research topic of utmost interest in recent years, with several studies highlighting its potential impact on physiology and, as a result, its possible effects on health and disease (Fan and Pedersen, 2021; Sarkar et al., 2021). Recently, substantial efforts have focused on studying infant gut microbiota with the emergence of the concept of Developmental

**OPEN ACCESS**

**Edited by:**
Huixia Yang,
Peking University, China

**Reviewed by:**
Bernard La Scola,
Aix-Marseille Université, France
Gregory R. Young,
Northumbria University,
United Kingdom

*Correspondence:
Julio Aires
julio.aires@u-paris.fr

**Specialty section:**
This article was submitted to Microorganisms in Vertebrate Digestive Systems, a section of the journal Frontiers in Microbiology

Received: 13 April 2022

Accepted: 22 June 2022

Published: 22 July 2022

**Citation:**
Toubon G, Butel M-J, Rozé J-C, Lepage P, Delannoy J, Ancel P-Y, Charles M-A and Aires J (2022) Very Preterm Children Gut Microbiota Comparison at the Neonatal Period of 1 Month and 3.5 Years of Life. *Front. Microbiol.* 13:919317. doi: 10.3389/fmicb.2022.919317

<footnote n="1"><sup>1</sup> INSERM, UMR1153 Centre de Recherche Épidémiologie et Statistiques (CRESS), Université Paris Cité, Paris, France,</footnote>

<footnote n="2"><sup>2</sup> Université Paris Cité, INSERM, UMR-S 1139, Physiopathologie et Pharmacotoxicologie Placentaire Humaine Microbiote Pré & Postnatal (3PHM), Paris, France, <sup>3</sup> FHU PREMA, Fighting Prematurity, Paris, France, <sup>4</sup> INRAE, UMR 1280, Physiologie des Adaptations Nutritionnelles (PhAN), Université hospitalière de Nantes, Nantes, France, <sup>5</sup> INRAE, UMR 1319, AgrosParisTech, Institut Micalis, Université Paris-Saclay, Paris, France</footnote>
Origins of Health and Disease (DOHaD) (Butel et al., 2018) and the critical first 1.000 days of life window hypothesis (Wopereis et al., 2014; Robertson et al., 2019). Prematurity, accounting for approximately 10% of births worldwide (Blencowe et al., 2012; Chawanpaiboon et al., 2019), is a risk factor for early gut dysbiosis, which is suggested to be associated with necrotizing enterocolitis (Rozé et al., 2017), growth delay (Tirone et al., 2019), cognitive impairment (Rozé et al., 2020), and brain damage (Seki et al., 2021). In very preterm neonates, gut microbiota is characterized by high inter-individual variability (Yee et al., 2019; Rozé et al., 2020) that can be attributed to varying birth conditions and extended exposure to the neonatal intensive care unit (NICU) environment and practices. Indeed, preterm infants are more often born via cesarean section and are frequently administered broad-spectrum antibiotics. Furthermore, gestational age influences the dynamics of the intestinal bacterial establishment (La Rosa et al., 2014). The early establishment of preterm gut microbiota is characterized by a dominance of opportunistic and potentially pathogenic bacteria, such as *Enterobacter*, *Enterococcus*, and *Staphylococcus*, and a delay in colonization by *Bifidobacterium* and *Bacteroides* (Arboleya et al., 2012; Barrett et al., 2013; La Rosa et al., 2014; Hill et al., 2017; Korpela et al., 2018). Several studies have highlighted the influence of the NICU environment on the establishment of preterm neonate gut microbiota (Hartz et al., 2015; Moles et al., 2015; D'Agata et al., 2019; Tauchi et al., 2019; Yap et al., 2021), notably during hospitalization (Rozé et al., 2020), and the existence of overlap of microbial taxa between NICU surfaces and the gut of preterm neonates (Brooks et al., 2014). Thus, while the infant gut microbiota establishment during the first weeks and months of life is documented, only a few studies have reported on the late development of gut microbiota in preterm children. Additionally, the relationship between gut microbiota and early-life neonatal factors, such as birth mode, breastfeeding, and antibiotic treatments after NICU discharge and during early childhood of preterm children, has been rarely studied, although it has been emphasized as a research priority (Groer et al., 2014). Indeed, to date, only a few studies with small cohorts have explored the development of gut microbiota in preterm children after 1 year of age (Gómez et al., 2017; Fouhy et al., 2019; Yee et al., 2019). These studies report an increase in alpha diversity over time and lower inter-individual variability in the gut microbiota between NICU discharge and early childhood. Furthermore, it has been suggested that gestational age at birth results in distinct microbial profiles extending up to 4 years of age (Fouhy et al., 2019). To the best of our knowledge, these data constitute the only available longitudinal evidence describing the gut microbiota of preterm children during early childhood. Therefore, given the paucity of data regarding this subject, there is a need to add substantial evidence related to prematurity and gut microbiota during early childhood, particularly in large cohorts.

This study aims to compare and characterize the gut microbiota of 159 preterm children from the EPIPAGE2 French nationwide cohort at two time points: at 1 month of life (NICU hospitalization) and at 3.5 years of age, when the child's gut microbiota supposedly stabilizes toward a mature adult-like gut microbiota (Stewart et al., 2018). Additionally, we investigated

factors previously reported to influence early gut microbiota development, such as gestational age (Groer et al., 2014; Gritz and Bhandari, 2015; Chernikova et al., 2018), birth mode (Dominguez-Bello et al., 2010; Portela et al., 2015; Stewart et al., 2018), birth weight (Groer et al., 2014; Unger et al., 2015), sex (Martin et al., 2016), antibiotics (Gibson et al., 2015; Tamburini et al., 2016), breastfeeding (Portela et al., 2015; Lemas et al., 2016; Pannaraj et al., 2017; Stewart et al., 2018), skin-to-skin contact (Hartz et al., 2015), and NICU strategies (Rozé et al., 2020). This study provides informative data concerning prematurity and the gut microbiota in early childhood based on a large and well-documented French multicentric cohort of preterm children.

# MATERIALS AND METHODS

## EPIFLORE Study and ELFE Cohort

The preterm children included in this study were part of the EPIFLORE study (Ancel and Goffinet, 2014), an ancillary study of the EPIPAGE2 French national birth cohort that was launched in 2011 and included preterm neonates born between 24 and 31 weeks of gestation (Lorthe et al., 2021). The aim of EPIFLORE was to investigate the development of the intestinal microbiota and its relationship with early-life factors that influence its establishment and subsequent health outcomes in early childhood among very preterm children (<32 weeks of gestational age). The EPIFLORE study included clinical data and stool samples collected from 24 French NICUs (Rozé et al., 2020). In the present study, full-term children who were born in 2011 from the ELFE French national birth cohort (Charles et al., 2020) were included as controls to investigate the differences between the gut microbiota of preterm and full-term children at 3.5 years of age. Full-term singletons born after 37 weeks of gestation with fecal samples available at 3.5 years of age and at least 2 years of subsequent follow-up were randomly selected from the ELFE cohort.

## Participants and Sample Collection

The EPIFLORE study included 729 neonates, and stool samples were collected from 574 preterm neonates during hospitalization at a median age of 23 days (interquartile range [IQR], 22-26 days), which is referred to as 1-month sampling. An additional stool sample was obtained for 208 preterm children at a median age of 42.6 months (IQR, 42.2-43.4 months), referred to as 3.5 years of age sampling. In the current study, 159 preterm children with fecal samples collected at both time points were included. Additionally, stool samples of 200 full-term children from the ELFE cohort were collected at 3.5 years of age (median age, 42.1 months; IQR, 41.4-43.6 months) during the same period and using the same methodology as for preterm children (Supplementary Figure 1).

## DNA Extraction, Sequencing, and Data Processing

The total fecal DNA was extracted according to the International Human Microbiome Standards operating procedure (Dore
et al., 2015), as previously reported (Rozé et al., 2020). We included negative controls that went through the same process as samples, from DNA extraction to sequencing. Sequencing was performed on the GeT-PlaGe platform of Génopole (Toulouse Midi-Pyrénées, France) using Illumina Miseq technology (V3, 2 × 250 bp) targeting the V3-V4 primers (V3fwd: ACGGRAGGCAGCAG, V4rev: AGGATTAGATACCCTGGTA) regions of the 16S bacterial rRNA gene. Raw sequences were analyzed using the pipeline "Find Rapidly OTU with Galaxy Solution" (FROGS) version 3.2 (September 2021) from the Galaxy software framework (Escudié et al., 2018). Sequences were checked for quality using FastQC (threshold 20-38). After trimming barcodes, the sequences were filtered (length of the reads ranged from 435 to 459 bp) and clustered into operational taxonomic units (OTUs) using the swarm clustering method implemented in the FROGS version 3.2 pipeline. OTUs representing less than 0.005% (788,210 reads) of all the sequences, the majority corresponding to singleton OTUs, were discarded (Bokulich et al., 2013). Taxonomic affiliation was assigned to the OTUs using the Silva 138.1 pintail100 database (September 2021) yielding 99% of the sequences affiliated with ≥99% identity and 100% coverage. A total of 11,526,400 reads (median, 22,619 reads per sample) were obtained from 16S rRNA gene sequencing. To adjust for the influence of uneven sampling depth, each sample was rarefied to the minimum sampling depth of the dataset (5,201 reads, **Supplementary Figure 2**), and rarefied data were used for all downstream analyses unless stated otherwise.

## Microbiota Analysis: Diversity and Composition

The Chao1 richness estimate (richness) and Shannon index (diversity) were calculated to investigate the patterns of microbial community diversity in the gut microbiota of preterm and full-term children. For beta diversity analysis, we computed dissimilarity matrices using Bray-Curtis and Unifrac distances. The Bray-Curtis dissimilarity distance reflects community composition considering the abundance of taxa, while the Unifrac distance considers the phylogenetic relationships among members of the bacterial communities. Gut microbiota inter-individual variability was visualized using principal coordinate analysis (PCoA/MDS) based on dissimilarity matrices. Stratification of the cohort based on gut microbial composition was performed according to the previously reported enterotyping guidelines (Arumugam et al., 2011). Briefly, the relative abundances of the classified genera were used to produce a Jensen-Shannon divergence (JSD) distance matrix between samples. Based on the obtained distance matrix, samples were clustered using Partitioning Around Medoids (PAM) algorithm. The optimal number of clusters was chosen by maximizing the Calinski-Harabasz index and was cross-validated with the silhouette index and prediction strength. The results of PAM clustering were visualized on PCoA biplots, and driver genera vectors were assessed using the "envfit" function. The vectors show the direction in the ordination space toward which genera change most rapidly and to which they have maximal

correlations within the ordination configuration. The genus with the highest relative abundance in each group was considered the main contributor to each enterotype.

## Neonatal Characteristics and Neonatal Intensive Care Unit Strategies

In each NICU, data on children and mothers during the perinatal and neonatal periods were collected prospectively during hospitalization until discharge. Available clinical data included gestational age, birth mode, birth weight, sex, antibiotics, breastfeeding, skin-to-skin contact, and NICU strategies. The NICU strategies were previously characterized in the EPIPAGE2 study regarding early extubation or no intubation, use of sedation, direct breastfeeding, skin-to-skin practice, speed of progression of enteral feeding, and duration of primary and secondary antibiotic therapy (Rozé et al., 2020). All these practices concern the early period of hospitalization before the 1-month stool sampling period. The methods to characterize these strategies were applied as previously described (Rozé et al., 2020).

## Statistical Analysis

For the comparative analysis of alpha diversity and pairwise beta diversity between enterotypes and time points, Wilcoxon rank-sum tests and Kruskal-Wallis tests were performed. Hypothesis testing was carried out using permutational multivariate analysis of variance (PERMANOVA) based on 999 permutations to analyze the influence of neonatal factors on the gut microbiota and to compare the gut microbiota between preterm and full-term children. For all PERMANOVA analyses, Benjamini-Hochberg FDR (false discovery rate) correction was used to adjust for multiple testing. We implemented a differential abundance testing analysis in order to assess significant changes in abundance at the genus level between samples at 1 month and 3.5 years in the preterm population. As these methods can produce heterogeneous results (Nearing et al., 2021; Wallen, 2021), we used two different methods: Analysis of Compositions of Microbiomes with Bias Correction (ANCOM-BC) (Lin and Peddada, 2020) and ANOVA-Like Differential Expression tool for compositional data (ALDEX2) (Fernandes et al., 2014). Differences in abundance were based on non-rarefied counts for both ANCOM-BC and ALDEX2, as both methods used their own normalization methods, and only genera with an abundance of 0.1 in at least 1% of the samples were used. We considered only concordant results between the two methods. P-values were corrected using the Benjamini-Hochberg method. Correlations of alpha- and beta-diversity community structures between 1 month and 3.5 years of age samples were assessed as follows: pairwise Spearman rank correlation test for alpha-diversity indexes; Mantel test for beta-diversity matrices using the Spearman rank correlation method based on 999 permutations. The association between enterotypes at both ages was tested using Fisher's exact test. To analyze the association between microbiota enterotypes found at 3.5 years of age and both neonatal characteristics and NICU strategies (neonates hospitalized at day 7), we used the chi-squared tests for categorical variables and student tests for continuous variables in univariate analyses.
Mixed-effects logistic regression with a random hospital intercept to consider the correlation between newborns in the same NICU was performed in multivariate analyses. All analyses were performed using the R software version 4.0.4 (R Foundation) and the packages phyloseq (v1.34.0), vegan (v2.6-0), cluster (v2.1.2), clusterSim (v0.49-2), fpc (v 2.2-9), ade4 (v1.7-16), ggplot2 (v 3.3.5), and lme4 (v1.1-27.1).

# RESULTS

## Enrolled Children Population

Our analyzed sample of preterm children (n = 159) was compared with the total population enrolled in the EPIFLORE study (n = 570), without stool samples available at the two time points, that is, 1 month and 3.5 years of age (**Supplementary Table 1**). We found that the characteristics of our study population did not significantly differ from the total population enrolled in the EPIFLORE study except in terms of the mother's age, the rate of mothers born outside France, and maternal level of education (**Supplementary Table 1**).

## Preterm Gut Microbiota at 1 Month of Age

Among the 159 stool samples analyzed at 1 month, 18 were non-amplifiable. The 141 neonates for whom sequencing data were obtained were partitioned into five bacterial enterotypes similar to those previously described by Rozé et al. (2020). In the present study, enterotypes 1–5 were characterized by a dominance of *Enterobacter* (n = 69), *Clostridium sensu stricto* 1 (n = 18), *Escherichia* (n = 25), *Enterococcus* (n = 18), and *Staphylococcus* (n = 11), respectively (**Figures 1A,B**). At the OTU level, enterotypes 1–5 were characterized by the dominance of OTU_1_E. hormaechei, OTU_13_C. neonatale, OTU_2_E. coli, OTU_4_E. faecalis, and OTU_19_S. caprae, respectively (**Supplementary Figure 3**). Enterotype 2 (dominance of *Clostridium sensu stricto* 1) and enterotype 3 (dominance of *E. coli*) were associated with greater gestational age (**Supplementary Figure 4**). Enterotype 3 was considered to be the more mature microbiota and defined as the reference enterotype at 1 month, as previously reported (Rozé et al., 2020). Concerning the 18 non-amplifiable samples, they represented the enterotype 6 gathering infants with a low bacterial load (Rozé et al., 2020).

Community complexities of the five different enterotypes were compared at the genus level according to alpha (Chao 1 richness and Shannon index) and beta (Bray-Curtis and Unifrac distances) diversities (**Figure 1C**). Significant differences were observed in alpha diversity across all enterotypes (Shannon diversity index, *p* = 0.003). Compared to reference enterotype 3, enterotypes 4 and 5 showed a significantly lower alpha diversity (Shannon index, *p* < 0.05 and *p* < 0.001, respectively). However, pairwise beta-diversity comparisons between enterotypes showed significant overall differences across groups (Bray-Curtis and Unifrac distance dissimilarities, *p* < 10⁻⁴ for both). Higher distance dissimilarities among infants belonging to enterotypes 1 and 2 were observed when compared to the reference enterotype

3 (Bray-Curtis, *p* < 0.0001 for both enterotypes), suggesting higher inter-individual variation. Enterotypes 4 and 5 were characterized by lower pairwise distance dissimilarities within groups (Bray-Curtis, [*p* < 0.0001], Unifrac [*p* = 0.007] for Ent4, and Bray-Curtis [*p* < 0.0001] for Ent5). At the OTU level, significant differences were observed in alpha diversity across all enterotypes (Chao1 and Shannon diversity index, *p* = 0.005 and *p* = 0.008, respectively). Differences at the OTU level are further presented in the **Supplementary Figure 3**.

## Preterm Gut Microbiota at 3.5 Years of Age

At 3.5 years of age, the gut microbiota of the 159 preterm children was characterized by the presence of six different phyla as follows: *Actinobacteriota*, *Bacteroidota*, *Desulfobacterota*, *Firmicuteota*, *Fusobacteriota*, and *Proteobacteriota*, with *Bacteroidota* and *Firmicuteota* being the most abundant phyla (**Supplementary Figure 5**). Concerning the stratification of the gut microbiota, it is optimally separated into two distinct enterotypes (**Figure 2A**). Clustering based on the genus taxonomic level showed that *Prevotella* had the strongest positive correlation and *Bacteroides* had the strongest negative correlation with the PCoA 1 dimension (**Figure 2A**). The enterotypes were enriched in either *Prevotella* (*P_type*) (n = 31) or *Bacteroides* (*B_type*) (n = 128) (**Figure 2B**). Interestingly, *Prevotella* was depleted in most children with *B_type* stools, with a median abundance of 0.02% (IQR, 0.02–0.06%), whereas *Bacteroides* remained the second most abundant genus in children with *P_type* stools, with a median abundance of 14% (IQR, 6–19%). At the OTU level, the *B_type* enterotype was characterized by the dominance of OTU_3_B. vulgatus, and the *P_type* enterotype was characterized by the dominance of OTU_30_P. copri (**Supplementary Figure 6**).

When comparing the microbiota of both *P_type* and *B_type* enterotypes according to alpha and beta diversities, the *P_type* enterotype showed greater richness (Chao1 index, *p* = 0.031, genus level) (**Figure 2C**). At the OTU level, no difference in richness was observed between the two enterotypes (**Supplementary Figures 6A,C**).

The Bray-Curtis distance pairwise dissimilarities between *P_type* and *B_type* enterotypes based on the genus taxonomic profile showed that *B_type* and *P_type* children shared a similar bacterial community within each enterotype (**Figure 2D**). The dissimilarity, indicated by the heatmap color and boxplot of pairwise dissimilarities, was significantly higher in the *P_type* enterotype than in the *B_type* enterotype, suggesting that the *P_type* community was more heterogeneous among children (**Figure 2D**). In contrast, at the OTU level, the *P_type* community was less heterogeneous and was characterized by lower pairwise dissimilarity distances (**Supplementary Figure 6**).

## Preterm and Full-Term Children Gut Microbiota Comparison at 3.5 Years of Age

To investigate the differences between the gut microbiota of preterm and full-term children at 3.5 years of age, we included
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Gut microbiota composition and diversity of the 141 preterm infants at 1 month after birth. (A) Clustering based on the genus taxonomic profiles. Biplot arrows indicate the top six genera that drive the sample to different locations on the plot. (B) The boxplots represent the relative abundance of the top 10 genera distributed among the five enterotypes. The sixth enterotype is not represented because it is constituted by infants without sample amplification (low bacterial load). The boxplots show the smallest and largest values, 25 and 75% quartiles, the median, and outliers. (C) Boxplots of alpha diversity (top) assessed by Chao1 and Shannon index, and pairwise beta diversity dissimilarities (bottom) assessed by Bray-Curtis and Unifrac distances within enterotypes based on genus taxonomic profiles (*p* < 0.05, **p** < 0.001, ns = *p* > 0.05, Reference group = Ent3); Ent, Enterotype.</figcaption></figure>

200 full-term children (ELFE cohort). Full-term children were significantly different in terms of gestational age, delivery mode, birth weight, the rate of mothers born outside of France, and the level of mother's education (Supplementary Table 2). *Firmicuteota* and *Bacteroidota* were the main phyla identified in the gut microbiota of full-term children. Moreover, the top 10 genera found in preterm children (except for the genus *Blautia*) also characterized the gut microbiota of full-term children (Figure 3A). Enterotyping at the genus taxonomic level showed that the gut microbiota in full-term and preterm children was optimally clustered into two discrete B_type (*n* = 293) and P_type (*n* = 66) enterotypes, with the two largest predictors being *Bacteroides* and *Prevotella*, respectively (Figure 3B). Interestingly, full-term children showed the same patterns as those seen in preterm children, where *Prevotella* was depleted from the majority of B_type children [median abundance of 0.04% (IQR, 0.02–0.08%)], whereas *Bacteroides* remained the second most abundant genus in P_type children (median abundance of 14% [IQR, 9–19%]) (Figure 3A). No differences in enterotype proportion distribution were observed between preterm and full-term children (*p* = 1.000; Figure 3C).

Analysis of alpha diversity showed that the gut microbiota of preterm children was characterized by a significantly lower Shannon diversity index than that of full-term children at the genus (*p* = 0.012) (Figure 3D) and OTU levels (*p* = 0.037). Overall, there were no differences in richness assessed using the Chao1 index at both the genus and OTU levels.

Using both Bray-Curtis and Unifrac distances at the genus level, the PERMANOVA suggested differences in the microbial community between preterm and full-term children (*p* = 0.012 and *p* = 0.006, FDR-corrected, respectively, Supplementary Figure 7).

## Preterm Gut Microbiota Comparison at 1 Month and 3.5 Years of Age

Figures 1B, 2B present the top 10 abundant genera at 1 month and 3.5 years of age, respectively, showing that the major genera present in the gut microbiota of preterm children changed between the two ages. In addition to *Bifidobacterium* and *Bacteroides*, the major genera initially found at 1 month (*Enterobacter*, *Clostridium sensu stricto*,
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Gut microbiota composition and diversity of the 159 preterm children at 3.5 years of age. (A) Clustering based on the genus taxonomic profile. Biplot arrows indicate the top four genera that drive samples to different locations on the plot. (B) The boxplots represent the relative abundance of the top 10 genera distributed among the P_ and B_type enterotypes. The boxplots show the smallest and largest values, 25 and 75% quartiles, the median, and outliers. (C) Boxplot of alpha diversity assessed by Chao1 estimate based on the genus taxonomic profile. (D) Inter-individual dissimilarity of the genus community based on the Bray-Curtis distance and represented by the heat map and boxplot. Each cell in the heat map, ordered according to B_ and P_type enterotypes, represents the dissimilarity level as per the color scale beside the box plot (*p* < 0.05, **p** < 0.0001).</figcaption></figure>

Escherichia, Enterococcus, and Staphylococcus) were replaced by other genera (Agathobacter, Faecalibacterium, Roseburia, Prevotella, and Ruminococcus). Except for Bifidobacterium and Prevotella, all genera cited above were significantly differentially abundant between the two time points (Supplementary Table 3). Based on the genus composition, the alpha diversity (Chao1 richness and Shannon index) of the gut microbiota significantly increased between both ages (Figure 4A). Pairwise dissimilarities based on genus composition within each age group showed that the gut microbiota of children at 1 month of age harbored higher beta diversity (Bray-Curtis and Unifrac) than at 3.5 years of age (*p* < 0.0001 for both; Figure 4A). Investigation of the relationship between gut microbiota at 1 month and 3.5 years of age in preterm children showed the absence of correlations of diversities between the two age groups either according to alpha diversity (*rho* = -0.054, *p* = 0.521 [Chao1 index]; *rho* = 0.063, *p* = 0.458 [Shannon index]) or beta diversity (*r* = 0.008, *p* = 0.410 [Bray-Curtis]; *r* = 0.030, *p* = 0.176 [Unifrac]). Comparison of the proportions between the two 3.5-year and the six 1-month enterotypes showed that enterotypes 3, 4, and 5 tended to harbor

the smallest proportion of P_type, but overall, the association was not significant (Fisher's exact test, *p* = 0.490) (Figure 4B). Compared to the reference enterotype 3, the gut microbiota of children originally clustered into enterotype 2 presented a higher alpha diversity at 3.5 years of age, but no overall differences were observed (Shannon index, *p* = 0.355; Figure 4C). No overall differences were noted at the OTU level (*p* = 0.173).

## Preterm Gut Microbiota at 3.5 Years of Age and Neonatal Factor Relations

Since we found that enterotypes at 3.5 years of age explained most of the variance in the gut microbiota (PERMANOVA analysis; Supplementary Table 4), we further explored the association between neonatal characteristics and enterotypes. Higher birth weight was significantly associated with the P_type enterotype in the univariate analysis (Table 1). In the multivariate model adjusted for maternal age, education level, and country of birth, we included covariates known to influence early gut microbiota development, such as gestational age, birth
<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Clustering of the gut microbiota of preterm (n = 159) and full-term (n = 200) children at 3.5 years of age. (A) The boxplots represent the relative abundance of the top 10 genera distributed between the B_ and P_types among the full-term children. The boxplots show the smallest and largest values, 25 and 75% quartiles, the median, and outliers. (B) Gut microbiota clustering of the 159 preterm and 200 full-term children at 3.5 years of age based on the genus taxonomic profile. Biplot arrows indicate the top four genera that drive samples to different locations on the plot. (C) Distribution of enterotypes among full-term and preterm children. The distribution is expressed as the proportion of each enterotype among each group of children. (D) Boxplot of alpha diversity assessed by Shannon index based on the genus taxonomic profile between preterm and full-term children (\*p < 0.05). The boxplots show the smallest and largest values, 25 and 75% quartiles, the median, and outliers.</figcaption></figure>

mode, birth weight, sex, antibiotics, breastfeeding, skin-to-skin contact, and NICU strategies. The results showed that birth weight was no longer associated with enterotypes, and no other neonatal factors were associated with enterotypes (Supplementary Table 5). Concerning the NICU strategies, we observed a tendency between the practice of a low volume of enteral nutrition and an increased probability of belonging to the P_type enterotype (OR: 5.3, 95% confidence interval [0.9–30.8]) (Supplementary Table 5).

DISCUSSION

This study showed that the gut microbiota of preterm and full-term children at 3.5 years of age is characterized by two enterotypes dominated by either *Bacteroides* or *Prevotella*. Interestingly, we found that prematurity still imprints the gut

microbiota of children, as the microbiota of preterm children showed lower diversity and different community composition than that seen in full-term children's microbiota. Additionally, the gut microbiota at 3.5 years of age was not related to that at 1 month in preterm children.

We confirmed that the gut microbiota of preterm neonates at 1 month of age is highly variable and clustered into six distinct enterotypes. Preterm neonate enterotypes were driven by a specific genus (i.e., *Enterobacter*, *Clostridium sensu stricto* 1, *Escherichia*, *Enterococcus*, and *Staphylococcus*), except for the sixth enterotype, which was characterized by a low bacterial load, leading to the absence of DNA amplification. These findings are in agreement with those of previous studies indicating that the early intestinal microbiota composition in preterm neonates is characterized by a strong dominance of only a few taxa (Stewart et al., 2016; Korpela et al., 2018). In a previous study, one of four genera, *Bifidobacterium*, *Enterobacter*, *Staphylococcus*, or
A

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

B

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

C

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

**FIGURE 4** | Comparison of the gut microbiota of the same preterm children at 1 month and 3.5 years of age. (A) Boxplots of alpha (top) and beta diversities (bottom) based on both Chao1 and Shannon index and both Bray–Curtis and Unifrac dissimilarity distances, respectively, at genus taxonomic level in the 141 preterm children with sequencing data at two time points. (B) Distribution of B_ and P_type enterotypes among the six enterotypes at 1 month in 159 preterm children with a fecal sample at two time points. The values are expressed as proportions of each enterotype at 3.5 years of age among each enterotype at 1 month. (C) Boxplots of Shannon diversity based on genus level at 3.5 years of age according to enterotype at 1 month (*p* < 0.05, ns = *p* > 0.05; Reference group = Ent3, Ent, Enterotype). The boxplots show the smallest and largest values, 25 and 75% quartiles, the median, and outliers. *****p* < 0.0001.

*Enterococcus*, represented > 50% of the reads in a given sample (Korpela et al., 2018). The authors associated *Staphylococcus* and *Enterococcus* dominance with younger postmenstrual age. In our cohort, *Staphylococcus*- and *Enterococcus*-driven enterotypes were also characterized by a lower gestational age (**Supplementary Figure 4**). Additionally, we found a dominance of *Clostridium sensu stricto* 1 and *Escherichia*, and an absence of a *Bifidobacterium*-driven enterotype, even though *Bifidobacterium* was among the top 10 dominant genera at 1 month. Korpela et al. (2018) reported an association between postnatal age and *Bifidobacterium* colonization. As their study was conducted with stool samples collected up to 60 days after birth, the lack of a *Bifidobacterium*-driven enterotype in our study could be explained by the fact that in our cohort, stool samples were collected at a median age of 23 days (IQR, 22–26 days), which

is too early for this genus to reach a state of dominance. In a study by Stewart et al. (2016), the gut microbiota of preterm neonates, examined using stool samples collected from birth to 100 days of life, was grouped into six clusters characterized by the dominance of *Klebsiella*, *Staphylococcus*, *Enterococcus*, *Escherichia*, and *Bifidobacterium*, or dominance of both *Klebsiella* and *Enterococcus*.

The present study showed that alpha diversity of the gut microbiota in preterm children increases with age and that the gut microbiota community is significantly different between 1 month and 3.5 years of age. The latter is characterized by two discrete enterotypes driven by two dominant genera, *Bacteroides* and *Prevotella*. Previous studies monitoring the gut microbiota of preterm children between the NICU stay and 2 years of age also showed a significant increase in alpha diversity with age
<table><caption>**TABLE 1** | Univariate association between preterm children enterotypes at 3.5 years of age, neonatal characteristics, and NICU strategies.</caption><thead><tr><th rowspan="2">Variables</th><th colspan="2">Enterotype</th><th rowspan="2">P-value</th></tr><tr><th>B_type (n = 128)</th><th>P_type (n = 31)</th></tr></thead><tbody><tr><td><strong>Gestational age (weeks)</strong></td><td></td><td></td><td></td></tr><tr><td>Mean (SD)</td><td>28.7 (1.96)</td><td>29.3 (1.89)</td><td>0.117</td></tr><tr><td><strong>Maternal age (years)</strong></td><td></td><td></td><td></td></tr><tr><td>&lt;25</td><td>6 (4.7%)</td><td>3 (9.7%)</td><td>0.354</td></tr><tr><td>25-35</td><td>85 (66.4%)</td><td>23 (74.2%)</td><td></td></tr><tr><td>≥35</td><td>32 (25.0%)</td><td>5 (16.1%)</td><td></td></tr><tr><td>Missing</td><td>5 (3.9%)</td><td>0 (0%)</td><td></td></tr><tr><td><strong>Country of birth of the mother</strong></td><td></td><td></td><td></td></tr><tr><td>France</td><td>112 (87.5%)</td><td>26 (83.9%)</td><td>0.810</td></tr><tr><td>Other</td><td>16 (12.5%)</td><td>5 (16.1%)</td><td></td></tr><tr><td><strong>Maternal level of education</strong></td><td></td><td></td><td></td></tr><tr><td>&lt;High school</td><td>24 (18.8%)</td><td>11 (35.5%)</td><td>0.011</td></tr><tr><td>High school</td><td>19 (14.8%)</td><td>8 (25.8%)</td><td></td></tr><tr><td>High school diploma +1 +2</td><td>30 (23.4%)</td><td>7 (22.6%)</td><td></td></tr><tr><td>&gt;High school diploma +3</td><td>55 (43.0%)</td><td>4 (12.9%)</td><td></td></tr><tr><td>Missing</td><td>0 (0%)</td><td>1 (3.2%)</td><td></td></tr><tr><td><strong>Neonatal factors</strong></td><td></td><td></td><td></td></tr><tr><td>Male sex</td><td>66 (51.6%)</td><td>20 (64.5%)</td><td>0.272</td></tr><tr><td>Birth weight (g), Mean (SD)</td><td>1130 (335)</td><td>1280 (317)</td><td>0.021</td></tr><tr><td>Birth weight (Z-score), Mean (SD)*</td><td>-0.913 (1.38)</td><td>-0.483 (1.36)</td><td>0.123</td></tr><tr><td>C-section delivery mode</td><td>82 (64.1%)</td><td>16 (51.6%)</td><td>0.283</td></tr><tr><td><strong>Practice of skin-to-skin contact during the first week of life</strong></td><td></td><td></td><td></td></tr><tr><td>Started between 0 and 3 days</td><td>29 (22.7%)</td><td>8 (25.8%)</td><td>0.359</td></tr><tr><td>Started between 4 and 7 days</td><td>39 (30.5%)</td><td>13 (41.9%)</td><td></td></tr><tr><td>Not practiced</td><td>53 (41.4%)</td><td>9 (29.0%)</td><td></td></tr><tr><td>Missing</td><td>7 (5.5%)</td><td>1 (3.2%)</td><td></td></tr><tr><td>Antibiotic therapy during neonatal period**</td><td>104 (81.3%)</td><td>28 (90.3%)</td><td>1.000</td></tr><tr><td>Missing</td><td>11 (8.6%)</td><td>0 (0%)</td><td></td></tr><tr><td>Did the child receive human milk during the neonatal period**</td><td>109 (85.2%)</td><td>25 (80.6%)</td><td>0.796</td></tr><tr><td>Missing</td><td>3 (2.3%)</td><td>1 (3.2%)</td><td></td></tr><tr><td><strong>NICU strategy***</strong></td><td></td><td></td><td></td></tr><tr><td>Direct breastfeeding during the first week</td><td>11 (8.6%)</td><td>2 (6.5%)</td><td>0.980</td></tr><tr><td>Skin-to-skin contact during the first week</td><td>74 (57.8%)</td><td>21 (67.7%)</td><td>0.419</td></tr><tr><td>Longer duration of primary antibiotic therapy</td><td>38 (29.7%)</td><td>13 (41.9%)</td><td>0.273</td></tr><tr><td>Longer duration of secondary antibiotic therapy</td><td>60 (46.9%)</td><td>17 (54.8%)</td><td>0.551</td></tr><tr><td>Sedation during the first week</td><td>81 (63.3%)</td><td>20 (64.5%)</td><td>1.000</td></tr><tr><td>No intubation or extubation at day 1</td><td>46 (35.9%)</td><td>14 (45.2%)</td><td>0.457</td></tr><tr><td>Low volume of enteral nutrition at day 7</td><td>35 (27.3%)</td><td>14 (45.2%)</td><td>0.087</td></tr></tbody></table>

Complete case analysis. Data are presented as number of events (percentages). Missing data are noted if any (B_type, enterotype enriched in Bacteroides; P_type, enterotype enriched in Prevotella).

NICU, neonatal intensive care unit.

*Score based on Olsen curves.

**Neonatal period is defined as the first 28 days of life.

***Favorable strategy: the observed percentage was greater than the expected percentage of newborns receiving the treatment or practice.

P-values marked with bold indicate statistically significant p-values (p < 0.05).

(Gómez et al., 2017; Yee et al., 2019) and reported differential beta-diversity structures between NICU stay and 4 years of age (Yee et al., 2019). We observed greater variability in the diversity of preterm children's gut microbiota at 1 month than at 3.5 years; this finding is in accordance with previous reports (Barrett et al., 2013; Stewart et al., 2016; Korpela et al., 2018). This indicates that preterm neonates' early gut microbiota has

greater inter-individual variability and lower stability but evolves and stabilizes over time. Furthermore, between 1 month and 3.5 years, 68 genera were significantly differentially abundant (25 negatively and 43 positively associated with the gut microbiota at 3.5 years). Indeed, in addition to *Bifidobacterium* and *Bacteroides*, the major genera shifted from organisms previously found to be abundant in preterm neonates during the NICU stay to common
organisms found in adults, such as Faecalibacterium, Roseburia, and Ruminococcus (Roswall et al., 2021). Altogether, these data support the fact that the intestinal microbiota of children born preterm evolves toward a richer and more diverse adult-like microbiota over time.

Analyses of correlations between preterm children's gut microbiota diversity during the NICU stay and early childhood showed no significant relationship. Children with higher alpha diversities at 1 month of age did not harbor higher alpha diversities at 3.5 years of age. Based on beta diversity, we found no correlation when examining community composition among close-together (i.e., lower dissimilarity distance) preterm children at 1 month and at 3.5 years of age, demonstrating that close-together newborns do not conserve their proximity over time. Similarly, there was no association between preterm neonate enterotypes at 1 month and 3.5 years of age, suggesting that belonging to one of the six enterotypes at 1 month did not predispose to harbor a certain type of enterotype at 3.5 years. If clustering methods are sensitive to sample size (Koren et al., 2013; Henry et al., 2015; Preud'homme et al., 2021), in the present study, we demonstrated good reproducibility of the clustering at 1 month based on our population of 141 EPIFLORE preterm neonates (i.e., 159 neonates excluding the 18 preterm neonates belonging to enterotype 6) compared to the previous clustering performed on 484 EPIFLORE preterm neonates (Rozé et al., 2020). Because we identified the same clustering patterns based on smaller sample size, we can assume that the shift from six to two enterotypes that we observed is not an artifact of the smaller sample size of the present study.

In the present study, the analysis between preterm and full-
term children at 3.5 years of age showed that both populations
shared similar microbiota profiles characterized by Bacteroides-
and Prevotella-driven enterotypes. This clustering is consistent
with the gut microbiota studies among full-term children
(Bergström et al., 2014; Nakayama et al., 2015; Méndez-Salazar
et al., 2018; Xiao et al., 2021) and adults (Roager et al., 2014;
Levy et al., 2020). A comparison of the microbiota community
complexities at 3.5 years of age indicated that preterm and
full-term children shared a similar richness composition, but
preterm children had a lower alpha diversity and a different beta
diversity (Supplementary Figure 7). It is known that the early
gut microbiota differs between full-term and preterm children
because of differences in exposure to prenatal and postnatal
determinants (Groer et al., 2014; Tauchi et al., 2019; Aguilar-
Lopez et al., 2021). Nonetheless, little is known about how
early differences in gut colonization and establishment evolve
as preterm neonates grow older. Differences in gut microbiota
diversity between preterm and full-term children are observed
within the first 6 months of life and are either still present
(Vandenplas et al., 2020) or disappear (Dahl et al., 2018; Yap et al.,
2021) at 12 months of age. Recently, the transition of enterotypes
shaping the gut microbiota of children at an early age revealed
that gut microbiota maturity of the preterm children was delayed
up to 12 months compared to that of full-term children (Xiao
et al., 2021). The present study indicates that at 3.5 years of age,
the gut microbiota of the preterm children differs in terms of
composition and is less diverse than that of full-term children's

gut microbiota, suggesting that the transition recovery of the
initially delayed gut microbiota colonization in preterm children
still continues at 3.5 years of age. Our results are in accordance
with those of Fouhy et al. (2019), who reported that gestational
age imprints gut microbiota up to 4 years of age. Over the course
of 3.5 years, we observed the evolution of preterm children's
gut microbiota toward an enterotype-like gut microbiota similar
to those found in full-term children. Because there is currently
a lack of data on the development of preterm children's gut
microbiota during early childhood and its impact on later health
outcomes, the present study provides another insight into our
understanding of the impact of gestational age on gut microbiota
development up to 3.5 years of age.

The neonatal factors, such as gestational age, delivery mode,
breastfeeding, early antibiotic therapy medication, skin-to-skin
practice, and NICU practices, are known to influence early gut
microbiota colonization (Rodríguez et al., 2015; Martin et al.,
2016; Tamburini et al., 2016; Kumbhare et al., 2019; Rozé et al.,
2020; Vandenplas et al., 2020). However, little is known about
the long-term impact of these factors on the gut microbiota.
Therefore, we investigated the potential effects of these factors on
the gut microbiota of preterm children aged 3.5 years. We found
no effects on the gut microbiota composition or any significant
associations with the enterotype prevalence in preterm children
at 3.5 years of age, indicating that these factors may not persist
up to this age. Nonetheless, we observed a tendency between
a low enteral volume at 7 days and an increased probability
of belonging to the *Prevotella*-driven enterotype; the odds ratio
indicates a strong association, but the wide 95% confidence
interval indicates poor precision, which may be due to the limited
sample size. Previously, the EPIFLORE study reported that
gestational age and birth mode were associated with enterotypes
at 1 month of age. Furthermore, no assisted ventilation on day 1
was associated with a decreased risk of belonging to enterotypes
5 or 6, while sedation and low volume of enteral nutrition were
associated with an increased risk (Rozé et al., 2020), [enterotypes
5 and 6 were the less mature enterotypes characterized by a
smaller gestational age (Supplementary Figure 4)]. If NICU
practices are associated with preterm neonates' gut microbiota
during the NICU stay, some of these practices may remain weakly
associated with the microbiota 3.5 years later.

**Study Limitations and Strengths**

The use of the 16S rRNA gene sequencing approach limits this
study to the description of bacterial composition. Analysis using
shotgun metagenomics would allow greater sequencing depth,
thus providing additional information regarding microbiome
functional profiles and meaning at the species level of the
microbiome that can differ between preterm and full-term
populations, despite similar bacterial composition. Microbiota
analysis was only performed for two times, that is, at 1 month
and 3.5 years of age, which does not take into account possible
dynamic changes in microbiota between these two time points.
Moreover, we studied the influence of prenatal factors on the
gut microbiota at 3.5 years in a preterm population, but we
did not take into account more contemporary factors occurring
during childhood, such as the use of antibiotics or dietary
habits that might influence the gut microbiota at 3.5 years. Finally, the number of children included may represent a limitation to definitely conclude about the association between the gut microbiota at 3.5 years of age and neonatal factors, but a tendency was observed, indicating a potential lack of statistical power.

In terms of major strengths, the present study is rare in the number of preterm children included within a multi-year follow-up and was performed on a large nationwide multicentric cohort, providing an accurate description of the neonatal characteristics of preterm children and characterization of the NICU strategies during hospitalization. Therefore, it allows a confident description of the preterm gut microbiota at both 1 month and 3.5 years of age, when the children's gut microbiota is supposed to mature toward an adult-like microbiota.

## CONCLUSION

This study demonstrates that during the NICU stay, the gut microbiota of preterm neonates is characterized by high inter-individual variation and low stability shaped by the NICU environment. However, NICU-shaped gut microbiota was not predictive of microbial community composition in the same child at 3.5 years of age. At 3.5 years of age, preterm children share gut microbiota similar to that of children born full-term and correspond to a mature adult-like microbiota in terms of enterotypes. However, the gut microbiota of preterm children was characterized by a lower Shannon diversity and different overall microbial composition compared to that seen in the gut microbiota of full-term children, indicating that prematurity still imprints the gut microbiota up to 3.5 years of age. These results must be replicated in other large-scale longitudinal cohorts to confirm the relationship between gestational age and other perinatal factors and the gut microbiota during early childhood.

## DATA AVAILABILITY STATEMENT

Personal data of children from the ELFE and EPIPAGE2 cohorts cannot be made publicly available for ethical reasons. They are available upon reasonable request from the authors under data-security conditions. The 16S rRNA gene reads are publicly available from the National Center for Biotechnology Information (nih.gov) Sequence Read Archive (SRA) under the Bioproject accession number: PRJNA798897, https://www.ncbi.nlm.nih.gov/bioproject/PRJNA798897.

## ETHICS STATEMENT

This study was approved by the National Data Protection Authority (Commission Nationale de l'Informatique et des Libertes [CNIL]), by the national advisory committee on information processing in health research (Comité Consultatif sur le Traitement de l'Information en matière de Recherche dans

le domaine de la Santé [CCTIRS]) and by the Committee for the Protection of People Participating in Biomedical Research (Comité de Protection des Personnes [CPP]). Recruitment and data collection occurred only after the families had received information and agreed to participate in these cohorts with informed consent.

## THE EPIFLORE STUDY GROUP

The members of the group include: M-JB, JA, Clotilde Rousseau, and JD (Federation of University Hospital, PREMA, UMR-S 1139, Faculty of Pharmacy, INSERM and Paris Descartes University); PL, Joel Dore, PhD, Ziad Al Nabhani, PhD, Karine Le Roux, and Celine Monot (INRA, UMR 1319 MICALIS); J-CR (Department of Neonatal Medicine, Nantes University Hospital); P-YA, Laetitia MartinMarchand, and Melanie Durox (INSERM, U1153, Obstetrical, Perinatal and Pediatric Epidemiology Team, Epidemiology and Biostatistics Sorbonne); Alexandre Lapillonne, MD, PhD (Department of Neonatal Medicine, Assistance Publique Hopitaux de Paris, Necker Enfants Malades Hospital); Jean-Charles Picaud, MD, PhD (Department of Neonatal Medicine, Hopital de la Croix-Rousse, Hospices Civils de Lyon); Farid Boudred, MD, PhD (Department of Neonatology, Faculte de Medecine, Aix-Marseille Université); Delphine Mitanchez, MD, PhD (Division of Neonatology, Department of Perinatology, Armand Trousseau Hospital); Valerie Biran (Department of Neonatology, Université Paris 7, Robert-Debre Hospital, Assistance Publique Hopitaux de Paris); Laurent Storme, MD, PhD (Department of Neonatal Medicine, Lille University Hospital); Olivier Claris, MD, PhD (Mothers and Children Hospital, Hospices Civils de Lyon); Gilles Cambonie, MD, PhD (Department of Neonatal Medicine, Montpellier University Hospital); Cyril Flamant, MD, PhD (Department of Neonatal Medicine, Nantes University Hospital); Anne Sauret, MD (Department of Neonatal Medicine, Rennes University Hospital); Odile Dicky, MD (Department of Neonatal Medicine, Toulouse University Hospital); Geraldine Favrais, MD, PhD (Department of Neonatology, Tours University Hospital); Jean-Michel Hascoet, MD, PhD (Department of Neonatal Medicine, Nancy University Hospital); Geraldine Gascoin (Department of Neonatal Medicine, Angers University Hospital); Gerard Thiriez, MD, PhD (Department of Pediatrics, Besançon University Hospital); Luc Desfrere, MD, PhD (Department of Neonatal Medicine, Louis Mourier Hospital, Assistance Publique Hopitaux de Paris); Xavier Durrmeyer, MD, PhD (Department of Neonatal Medicine, Centre Hospitalier Intercommunal); and Clement Chollat, MD (Department of Neonatal Medicine, Cochin University Hospital).

## AUTHOR CONTRIBUTIONS

J-CR, M-JB, M-AC, P-YA, and JA conceived and designed the study. M-JB, PL, GT, JD, and J-CR generated the data. JD and GT analyzed the data. GT drafted the manuscript. M-JB, J-CR,
PL, JD, P-YA, M-AC, and JA critically revised the manuscript. JA, M-JB, and M-AC supervised the study. All authors approved the submitted version.

## FUNDING

This work was supported by l'Agence Nationale pour la Recherche under Grants (EPIFLORE: ANR-12-BSV3-0025 and ELFE: ANR-11-EQPX-0038), the Nestec Research Center (Constitution of the collection of stool samples during EPIPAGE2), and the Biostime Institute for Nutrition and Care—Geneva (BINC-Geneva) (16S rRNA gene sequencing at 3.5 years).

## REFERENCES

Aguilar-Lopez, M., Dinsmoor, A. M., Ho, T. T. B., and Donovan, S. M. (2021). A systematic review of the factors influencing microbial colonization of the preterm infant gut. *Gut Microbes* 13, 1–33. doi: 10.1080/19490976.2021.1884514

Ancel, P. Y., and Goffinet, F. (2014). EPIPAGE 2: a preterm birth cohort in France in 2011. *BMC Pediatr.* 14:97. doi: 10.1186/1471-2431-14-97

Arboleya, S., Binetti, A., Salazar, N., Fernández, N., Solís, G., Hernández-Barranco, A., et al. (2012). Establishment and development of intestinal microbiota in preterm neonates. *FEMS Microbiol. Ecol.* 79, 763–772. doi: 10.1111/j.1574-6941.2011.01261.x

Arumugam, M., Raes, J., Pelletier, E., Le Paslier, D., Yamada, T., Mende, D. R., et al. (2011). Enterotypes of the human gut microbiome. *Nature* 473, 174–180. doi: 10.1038/nature09944

Barrett, E., Kerr, C., Murphy, K., O'Sullivan, O., Ryan, C. A., Dempsey, E. M., et al. (2013). The individual-specific and diverse nature of the preterm infant microbiota. *Arch. Dis. Childhood Fetal Neonatal Ed.* 98, F334–F340. doi: 10.1136/archdischild-2012-303035

Bergström, A., Skov, T. H., Bahl, M. I., Roager, H. M., Christensen, L. B., Ejlerskov, K. T., et al. (2014). Establishment of intestinal microbiota during early life: a longitudinal, explorative study of a large cohort of danish infants. *Appl. Environ. Microbiol.* 80, 2889–2900. doi: 10.1128/AEM.00342-14

Blencowe, H., Cousens, S., Oestergaard, M. Z., Chou, D., Moller, A. B., Narwal, R., et al. (2012). National, regional, and worldwide estimates of preterm birth rates in the year 2010 with time trends since 1990 for selected countries: a systematic analysis and implications. *Lancet* 379, 2162–2172. doi: 10.1016/S0140-6736(12)60820-4

Bokulich, N. A., Subramanian, S., Faith, J. J., Gevers, D., Gordon, J. I., Knight, R., et al. (2013). Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. *Nat. Methods* 10, 57–59. doi: 10.1038/nmeth.2s276

Brooks, B., Firek, B. A., Miller, C. S., Sharon, I., Thomas, B. C., Baker, R., et al. (2014). Microbes in the neonatal intensive care unit resemble those found in the gut of premature infants. *Microbiome* 2:1. doi: 10.1186/2049-2618-2-1

Butel, M. J., Waligora-Dupriet, A. J., and Wydaud-Dematteis, S. (2018). The developing gut microbiota and its consequences for health. *J. Dev. Origins Health Dis.* 9, 590–597. doi: 10.1017/S2040174418000119

Charles, M. A., Thierry, X., Lanoe, J. L., Bois, C., Dufourg, M. N., Popa, R., et al. (2020). Cohort Profile: The French national cohort of children (ELFE): birth to 5 years. *Int. J. Epidemiol.* 49, 368–369. doi: 10.1093/ije/dyz227

Chawanpaiboon, S., Vogel, J. P., Moller, A. B., Lumbiganon, P., Petzold, M., Hogan, D., et al. (2019). Global, regional, and national estimates of levels of preterm birth in 2014: a systematic review and modelling analysis. *Lancet Glob. Health* 7, e37–e46. doi: 10.1016/S2214-109X(18)30451-0

Chernikova, D. A., Madan, J. C., Housman, M. L., Zain-ul-abideen, M., Lundgren, S. N., Morrison, H. G., et al. (2018). The premature infant gut microbiome during the first 6 weeks of life differs based on gestational maturity at birth. *Pediatr. Res.* 84, 71–79. doi: 10.1038/s41390-018-0022-z

D'Agata, A. L., Wu, J., Welandawe, M. K. V., Dutra, S. V. O., Kane, B., and Groer, M. W. (2019). Effects of early life NICU stress on the developing gut microbiome. *Dev. Psychobiol.* 61, 650–660. doi: 10.1002/dev.21826

## ACKNOWLEDGMENTS

We thank the Genotoul Get-PlaGe sequencing platform and GenoToul bioinformatics facility (http://bioinfo.genotoul.fr) for providing computational resources.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.919317/full#supplementary-material

Dahl, C., Stigum, H., Valeur, J., Iszatt, N., Lenters, V., Peddada, S., et al. (2018). Preterm infants have distinct microbiomes not explained by mode of delivery, breastfeeding duration or antibiotic exposure. *Int. J. Epidemiol.* 47, 1658–1669. doi: 10.1093/ije/dyy064

Dominguez-Bello, M. G., Costello, E. K., Contreras, M., Magris, M., Hidalgo, G., Fierer, N., et al. (2010). Delivery mode shapes the acquisition and structure of the initial microbiota across multiple body habitats in newborns. *Proc. Natl. Acad. Sci. U.S.A.* 107, 11971–11975. doi: 10.1073/pnas.1002601107

Dore, J., Ehrlich, S. D., Levenez, F., Pelletier, E., Alberti, A., Bertrand, L., et al. (2015). *HMS_SOP_07_V1: standard operating procedure for fecal samples DNA extraction, protocol H.* *International Human Microbiome Standards*. London: International SOS.

Escudié, F., Auer, L., Bernard, M., Mariadassou, M., Cauquil, L., Vidal, K., et al. (2018). Frogs: find, rapidly, otus with galaxy solution. *Bioinformatics* 34, 1287–1294. doi: 10.1093/bioinformatics/btx791

Fan, Y., and Pedersen, O. (2021). Gut microbiota in human metabolic health and disease. *Nat. Rev. Microbiol.* 19, 55–71. doi: 10.1038/s41579-020-0433-9

Fernandes, A. D., Reid, J. N., Macklaim, J. M., McMurrough, T. A., Edgell, D. R., and Gloor, G. B. (2014). Unifying the analysis of high-throughput sequencing datasets: characterizing RNA-seq, 16S rRNA gene sequencing and selective growth experiments by compositional data analysis. *Microbiome* 2:15. doi: 10.1186/2049-2618-2-15

Fouhy, F., Watkins, C., Hill, C. J., O'Shea, C. A., Nagle, B., Dempsey, E. M., et al. (2019). Perinatal factors affect the gut microbiota up to four years after birth. *Nat. Commun.* 10:1517. doi: 10.1038/s41467-019-09252-4

Gibson, M. K., Crofts, T. S., and Dantas, G. (2015). Antibiotics and the developing infant gut microbiota and resistome. *Curr. Opin. Microbiol.* 27, 51–56. doi: 10.1016/j.mib.2015.07.007

Gómez, M., Moles, L., Espinosa-Martos, I., Bustos, G., de Vos, W., Fernández, L., et al. (2017). Bacteriological and immunological profiling of meconium and fecal samples from preterm infants: a two-year follow-up study. *Nutrients* 9:1293. doi: 10.3390/nu9121293

Gritz, E. C., and Bhandari, V. (2015). The human neonatal gut microbiome: a brief review. *Front. Pediatr.* 3:17. doi: 10.3389/fped.2015.00017

Groer, M. W., Luciano, A. A., Dishaw, L. J., Ashmeade, T. L., Miller, E., and Gilbert, J. A. (2014). Development of the preterm infant gut microbiome: a research priority. *Microbiome* 2:38. doi: 10.1186/2049-2618-2-38

Hartz, L. E., Bradshaw, W., and Brandon, D. H. (2015). Potential NICU environmental influences on the neonate's microbiome: a systematic review. *Adv. Neonatal Care* 15, 324–335. doi: 10.1097/ANC.0000000000000220

Henry, D., Dymnicki, A. B., Mohatt, N., Allen, J., and Kelly, J. G. (2015). Clustering methods with qualitative data: a mixed-methods approach for prevention research with small samples. *Prev. Sci.* 16, 1007–1016. doi: 10.1007/s11121-015-0561-z

Hill, C. J., Lynch, D. B., Murphy, K., Ulaszewska, M., Jeffery, I. B., O'Shea, C. A., et al. (2017). Evolution of gut microbiota composition from birth to 24 weeks in the INFANTMET Cohort. *Microbiome* 5:4. doi: 10.1186/s40168-016-0213-y

Koren, O., Knights, D., Gonzalez, A., Waldron, L., Segata, N., Knight, R., et al. (2013). A guide to enterotypes across the human body: meta-analysis of
microbial community structures in human microbiome datasets. *PLoS Comput. Biol.* 9:e1002863. doi: 10.1371/journal.pcbi.1002863

Korpela, K., Blakstad, F. W., Moltu, S. J., Strömmen, K., Nakstad, B., Rønnestad, A. E., et al. (2018). Intestinal microbiota development and gestational age in preterm neonates. *Sci. Rep.* 8:2453. doi: 10.1038/s41598-018-20827-x

Kumbhare, S. V., Patangia, D. V., Patil, R. H., Shouche, Y. S., and Patil, N. P. (2019). Factors influencing the gut microbiome in children: from infancy to childhood. *J. Biosci.* 44:49. doi: 10.1007/s12038-019-9860-z

La Rosa, P. S., Warner, B. B., Zhou, Y., Weinstock, G. M., Sodergren, E., Hall-Moore, C. M., et al. (2014). Patterned progression of bacterial populations in the premature infant gut. *Proc. Natl. Acad. Sci. U.S.A.* 111, 12522–12527. doi: 10.1073/pnas.1409497111

Lemas, D. J., Yee, S., Cacho, N., Miller, D., Cardel, M., Gurka, M., et al. (2016). Exploring the contribution of maternal antibiotics and breastfeeding to development of the infant microbiome and pediatric obesity. *Semin. Fetal Neonatal Med.* 21, 406–409. doi: 10.1016/j.siny.2016.04.013

Levy, R., Magis, A. T., Earls, J. C., Manor, O., Wilmanski, T., Lovejoy, J., et al. (2020). Longitudinal analysis reveals transition barriers between dominant ecological states in the gut microbiome. *Proc. Natl. Acad. Sci. U.S.A.* 117, 13839–13845. doi: 10.1073/pnas.1922498117

Lin, H., and Peddada, S. D. (2020). Analysis of compositions of microbiomes with bias correction. *Nat. Commun.* 11:3514. doi: 10.1038/s41467-020-17041-7

Lorthe, E., Benhammou, V., Marchand-Martin, L., Pierrat, V., Lebeaux, C., Durox, M., et al. (2021). Cohort Profile: The Etude Epidémiologique sur les Petits Ages Gestationnels-2 (EPIPAGE-2) preterm birth cohort. *Int. J. Epidemiol.* 50, 1428–1429m. doi: 10.1093/ije/dyaa282

Martin, R., Makino, H., Cetinyurek Yavuz, A., Ben-Amor, K., Roelofs, M., Ishikawa, E., et al. (2016). Early-life events, including mode of delivery and type of feeding, siblings and gender, shape the developing gut microbiota Suchodolski. *PLoS One* 11:e0158498. doi: 10.1371/journal.pone.0158498

Méndez-Salazar, E. O., Ortiz-López, M. G., Granados-Silvestre, M., de los Á, Palacios-González, B., and Menjivar, M. (2018). Altered gut microbiota and compositional changes in firmicutes and *proteobacteria* in mexican undernourished and obese children. *Front. Microbiol.* 9:2494. doi: 10.3389/fmicb.2018.02494

Moles, L., Gómez, M., Jiménez, E., Fernández, L., Bustos, G., Chaves, F., et al. (2015). Preterm infant gut colonization in the neonatal ICU and complete restoration 2 years later. *Clin. Microbiol. Infect.* 21:936.e1-10. doi: 10.1016/j.cmi.2015.06.003

Nakayama, J., Watanabe, K., Jiang, J., Matsuda, K., Chao, S.-H., Haryono, P., et al. (2015). Diversity in gut bacterial community of school-age children in Asia. *Sci. Rep.* 5:8397. doi: 10.1038/srep08397

Nearing, J. T., Douglas, G. M., Hayes, M., MacDonald, J., Desai, D., Allward, N., et al. (2021). Microbiome differential abundance methods produce disturbingly different results across 38 datasets. *Nat. Commun.* [Preprint]. doi: 10.1101/2021.05.10.443486

Pannaraj, P. S., Li, F., Cerini, C., Bender, J. M., Yang, S., Rollie, A., et al. (2017). Association between breast milk bacterial communities and establishment and development of the infant gut microbiome. *JAMA Pediatr.* 171:647. doi: 10.1001/jamapediatrics.2017.0378

Portela, D. S., Vieira, T. O., Matos, S. M., de Oliveira, N. F., and Vieira, G. O. (2015). Maternal obesity, environmental factors, cesarean delivery and breastfeeding as determinants of overweight and obesity in children: results from a cohort. *BMC Pregnancy Childbirth* 15:94. doi: 10.1186/s12884-015-0518-z

Preud'homme, G., Duarte, K., Dalleau, K., Lacomblez, C., Bresso, E., Smaïl-Tabbone, M., et al. (2021). Head-to-head comparison of clustering methods for heterogeneous data: a simulation-driven benchmark. *Sci. Rep.* 11:4202. doi: 10.1038/s41598-021-83340-8

Roager, H. M., Licht, T. R., Poulsen, S. K., Larsen, T. M., and Bahl, M. I. (2014). Microbial Enterotypes, Inferred by the Prevotella-to-*Bacteroides* Ratio, Remained Stable during a 6-month Randomized Controlled Diet Intervention with the New Nordic Diet. *Appl. Environ. Microbiol.* 80, 1142–1149. doi: 10.1128/AEM.03549-13

Robertson, R. C., Manges, A. R., Finlay, B. B., and Prendergast, A. J. (2019). The human microbiome and child growth – first 1000 days

and beyond. *Trends Microbiol.* 27, 131–147. doi: 10.1016/j.tim.2018.09.008

Rodríguez, J. M., Murphy, K., Stanton, C., Ross, R. P., Kober, O. I., Juge, N., et al. (2015). The composition of the gut microbiota throughout life, with an emphasis on early life. *Microb. Ecol. Health Dis.* 26:26050. doi: 10.3402/mehd.v26.26050

Roswall, J., Olsson, L. M., Kovatcheva-Datchary, P., Nilsson, S., Tremaroli, V., Simon, M.-C., et al. (2021). Developmental trajectory of the healthy human gut microbiota during the first 5 years of life. *Cell Host Microbe* 29, 765–776.e3. doi: 10.1016/j.chom.2021.02.021

Rozé, J.-C., Ancel, P.-Y., Lepage, P., Martin-Marchand, L., Al Nabhani, Z., Delannoy, J., et al. (2017). Nutritional strategies and gut microbiota composition as risk factors for necrotizing enterocolitis in very-preterm infants. *Am. J. Clin. Nutr.* 106, 821–830. doi: 10.3945/ajcn.117.152967

Rozé, J. C., Ancel, P. Y., Marchand-Martin, L., Rousseau, C., Montassier, E., Monot, C., et al. (2020). Assessment of neonatal intensive care unit practices and preterm newborn gut microbiota and 2-year neurodevelopmental outcomes. *JAMA Netw. Open* 3:e2018119. doi: 10.1001/jamanetworkopen.2020.18119

Sarkar, A., Yoo, J. Y., Valeria Ozorio Dutra, S., Morgan, K. H., and Groer, M. (2021). The Association between Early-Life Gut Microbiota and Long-Term Health and Diseases. *J. Clin. Med.* 10:459. doi: 10.3390/jcm10030459

Seki, D., Mayer, M., Hausmann, B., Pjevac, P., Giordano, V., Goeral, K., et al. (2021). Aberrant gut-microbiota-immune-brain axis development in premature neonates with brain damage. *Cell Host Microbe* 29, 1558–1572.e6. doi: 10.1016/j.chom.2021.08.004

Stewart, C. J., Ajami, N. J., O'Brien, J. L., Hutchinson, D. S., Smith, D. P., Wong, M. C., et al. (2018). Temporal development of the gut microbiome in early childhood from the TEDDY study. *Nature* 562, 583–588. doi: 10.1038/s41586-018-0617-x

Stewart, C. J., Embleton, N. D., Marrs, E. C. L., Smith, D. P., Nelson, A., Abdulkadir, B., et al. (2016). Temporal bacterial and metabolic development of the preterm gut reveals specific signatures in health and disease. *Microbiome* 4:67. doi: 10.1186/s40168-016-0216-8

Tamburini, S., Shen, N., Wu, H. C., and Clemente, J. C. (2016). The microbiome in early life: implications for health outcomes. *Nat. Med.* 22, 713–722. doi: 10.1038/nm.4142

Tauchi, H., Yahagi, K., Yamauchi, T., Hara, T., Yamaoka, R., Tsukuda, N., et al. (2019). Gut microbiota development of preterm infants hospitalised in intensive care units. *Benef. Microbes* 10, 641–651. doi: 10.3920/BM2019.0003

Tirone, C., Pezza, L., Paladini, A., Tana, M., Aurilia, C., Lio, A., et al. (2019). Gut and lung microbiota in preterm infants: immunological modulation and implication in neonatal outcomes. *Front. Immunol.* 10:2910. doi: 10.3389/fimmu.2019.02910

Unger, S., Spintzi, A., Shah, P., Mack, D., and O'Connor, D. L. (2015). Gut microbiota of the very-low-birth-weight infant. *Pediatr. Res.* 77, 205–213. doi: 10.1038/pr.2014.162

Vandenplas, Y., Carnielli, V. P., Ksiazyk, J., Luna, M. S., Migacheva, N., Mosselmans, J. M., et al. (2020). Factors affecting early-life intestinal microbiota development. *Nutrition* 78:110812. doi: 10.1016/j.nut.2020.110812

Wallen, Z. D. (2021). Comparison study of sixteen differential abundance testing methods using two large Parkinson disease gut microbiome datasets. *Genomics* [Preprint]. doi: 10.1101/2021.02.24.432717

Wopereis, H., Oozeer, R., Knipping, K., Belzer, C., and Knol, J. (2014). The first thousand days - intestinal microbiology of early life: establishing a symbiosis. *Pediatric Aller. Immunol.* 25, 428–438. doi: 10.1111/pai.12232

Xiao, L., Wang, J., Zheng, J., Li, X., and Zhao, F. (2021). Deterministic transition of enterotypes shapes the infant gut microbiome at an early age. *Genome Biol.* 22:243. doi: 10.1186/s13059-021-02463-3

Yap, P. S. X., Chong, C. W., Ahmad Kamar, A., Yap, I. K. S., Choo, Y. M., Lai, N. M., et al. (2021). Neonatal Intensive Care Unit (NICU) exposures exert a sustained influence on the progression of gut microbiota and metabolome in the first year of life. *Sci. Rep.* 11:1353. doi: 10.1038/s41598-020-80278-1

Yee, A. L., Miller, E., Dishaw, L. J., Gordon, J. M., Ji, M., Dutra, S., et al. (2019). Longitudinal microbiome composition and stability correlate with increased
weight and length of very-low-birth-weight infants. *mSystems* 4:e00229-18. doi:
10.1128/mSystems.00229-18

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

**Publisher's Note:** All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in

this article, or claim that may be made by its manufacturer, is not guaranteed or
endorsed by the publisher.

Copyright © 2022 Toubon, Butel, Rozé, Lepage, Delannoy, Ancel, Charles and Aires. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.