# Evolution of Gut Microbiome and Metabolome in Suspected Necrotizing Enterocolitis: A Case-Control Study

Camille Brehin, Damien Dubois, Odile Dicky, Sophie Breinig, Eric Oswald, Matteo Serino

► To cite this version:

Camille Brehin, Damien Dubois, Odile Dicky, Sophie Breinig, Eric Oswald, et al.. Evolution of Gut Microbiome and Metabolome in Suspected Necrotizing Enterocolitis: A Case-Control Study. Journal of Clinical Medicine, 2020, 9 (7), pp.2278. 10.3390/jcm9072278 . hal-02968291

HAL Id: hal-02968291

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

Submitted on 15 Oct 2020

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>

Article

# Evolution of Gut Microbiome and Metabolome in Suspected Necrotizing Enterocolitis: A Case-Control Study

Camille Brehin <fnref n="1" />,<fnref n="2" />, Damien Dubois <fnref n="2" />,<fnref n="3" />, Odile Dicky <fnref n="4" />, Sophie Breinig <fnref n="5" /> ID, Eric Oswald <fnref n="2" />,<fnref n="3" /> and Matteo Serino <fnref n="2" />,* ID

Received: 2 June 2020; Accepted: 15 July 2020; Published: 17 July 2020

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

check for
updates

**Abstract:** Background: Necrotizing enterocolitis (NEC) is a devastating condition in preterm infants due to multiple factors, including gut microbiota dysbiosis. NEC development is poorly understood, due to the focus on severe NEC (NEC-2/3). Methods: We studied the gut microbiota, microbiome and metabolome of children with suspected NEC (NEC-1). Results: NEC-1 gut microbiota had a higher abundance of the Streptococcus (second 10-days of life) and Staphylococcus (third 10-days of life) species. NEC-1 children showed a microbiome evolution in the third 10-days of life being the most divergent, and were associated with a different metabolomic signature than in healthy children. The NEC-1 microbiome had increased glycosaminoglycan degradation and lysosome activity by the first 10-days of life, and was more sensitive to childbirth, low birth weight and gestational age, than healthy microbiome. NEC-1 fecal metabolome was more divergent by the second month of life. Conclusions: NEC-1 gut microbiota and microbiome modifications appear more distinguishable by the third 10-days of life, compared to healthy children. These data identify a precise window of time (i.e., the third 10-days of life) and provide microbial targets to fight/blunt NEC-1 progression.

**Keywords:** necrotizing enterocolitis; intestinal microbiology; microbiome; infant gut; metabolomics

# 1. Introduction

Necrotizing enterocolitis (NEC), defined by the Bell classification [1–3], is the most severe intestinal disease in preterm infants, with a mortality score of 25% and long-term neurological morbidity [4]. Yet, a precise initiating factor of this pathology is missing. In the last decade, gut microbiota was identified and recognized as a specific organ with functions widely beyond digestion [5]. Both its taxonomic (relative abundance) and functional (microbial pathway) alterations, named dysbiosis, were described in several pathologies, in particular, metabolic diseases such as type 2 diabetes and obesity [6–8], and intestinal inflammatory diseases [9,10]. Importantly, a dysbiotic gut microbiota associated with a very high inflammatory status of the gut [11,12] may trigger NEC development, since germ-free mice do not develop NEC [13].

<footnote n="1">¹ General Pediatrics Department, Children Hospital, Toulouse University Hospital, 31300 Toulouse, France;
brehin.c@chu-toulouse.fr</footnote>

<footnote n="2">² IRSD, Université de Toulouse, INSERM, INRAE, ENVT, UPS, 31024 Toulouse, France;
dubois.d@chu-toulouse.fr (D.D.); eric.oswald@inserm.fr (E.O.)</footnote>

<footnote n="3">³ Service of Bacteriology-Hygiene, Toulouse University Hospital, 31300 Toulouse, France</footnote>

<footnote n="4">⁴ Department of Neonatology, Children Hospital, University Hospital, UMR 1027, INSERM, Paul Sabatier
University, 31000 Toulouse, France; dicky.o@chu-toulouse.fr</footnote>

<footnote n="5">5 Neonatal and Pediatric Intensive Care Unit, Toulouse University Hospital, 31300 Toulouse, France;
breinig.s@chu-toulouse.fr</footnote>

<footnote>* Correspondence: matteo.serino@inserm.fr; Tel.: +33-(0)5-6274-4525</footnote>
From a clinical and microbiological point of view, studies of NEC were focused only on established and severe phenotypes, such as NEC-2 and NEC-3. Based on the French study Etude épidémiologique sur les petits âges gestationnels (EPIPAGE 2), the incidence of proved NEC-2 and NEC-3 is 1–5% in preterm infants born at less than 32 weeks of gestation [14].

By contrast, NEC suspicions such as lethargy, bradycardia, thermic instability associated with biliary gastric residues, vomiting, abdominal distension with or without rectal bleeding, with a normal abdominal x-ray image or a simple dilatation, which identifies suspected NEC (NEC-1), have not been studied yet. In fact, enteropathies are frequent in the first weeks of life in preterm infants, though no data are available about NEC-1 incidence. This induces the end of alimentation, a prolonged (sometime life-lasting) parenteral nutrition, with a delayed gut maturation and failure to thrive [15]. Therefore, to study the evolution of gut microbiota and microbiome during the early onset of NEC, we focused on NEC-1 children within the first two months of life. We studied the fecal metabolome, to understand how a change in gut microbiota may drive alterations in intestinal metabolites. To further understand which factor of mother and child may affect the evolution of gut microbiota, microbiome and fecal metabolome during NEC-1, we analyzed: the presence of neonatal antibiotherapy (ABx), ABx treatment on the mother, childbirth (Cesarean-section (C-sec) vs. vaginal birth (VB)), very low birth weight (VLBW), extreme low birth weight (ELBW) and gestational age (GA) > or ≤ 28 weeks.

## 2. Research Design and Methods

### 2.1. Cohort Constitution

We conducted a prospective monocentric case-control cohort study. This study was approved (number of the approval: DC 2016-2804) by Neonatal and Pediatric Intensive Care Unit and Neonatology Department of Purpan Hospital in Toulouse, France. The parents of the children involved in this study gave their approval by written consent. The inclusion criteria regarding all of the children hospitalized into the Neonatal and Pediatric Intensive Care Unit or Neonatology Departments of the Purpan Hospital, were:

- Newborn of gestational age under 34 weeks of gestation;

- Diagnosis of suspected necrotizing enterocolitis (NEC-1) made by a neonatologist;

- Obtainment of the non-opposition from parents or their legal representative.

Following the inclusion of every case, we conducted in parallel a search for two controls, according to the following matching criteria, listed in decreasing priority:

- Gestational age (±1 week of gestation, priority to matched age);

- Body weight;

- Neonatal antibiotherapy;

- Childbirth (C-section vs. vaginal);

- Maternal antibiotherapy.

Inclusion criteria for controls were:

- Newborn of gestational age under 34 weeks of gestation;

- Respect of the matching according to the priority order of the established criteria;

- Obtainment of the oral non-opposition from parents or their legal representative.

Children with complex congenital cardiopathy or with spontaneous intestinal perforation without a radiological evidence of NEC were excluded from the study.

Based on these criteria, we included 11 NEC-1 children, with 27 feces collection (4 fecal samples for time-point 1–10 days (d); 10 fecal samples for time-point 11–20 d; 7 fecal samples for time-point 21–30 d; 6 fecal samples for time-point > 30 d) and 21 healthy children, with 53 feces collection (15 fecal samples for time-point 1–10 d; 14 fecal samples for time-point 11–20 d; 13 fecal samples for time-point
21–30 d; 11 fecal samples for time-point > 30 d). Hence, a total of 80 fecal samples was analyzed in our study. The period of collection was day 1 to day 68 of life of the newborn.

## 2.2. Taxonomic and Functional Analysis of Gut Microbiota

Feces analyzed in this study were collected by nurses in the related department in the first week of life and once a week till the end of the hospitalization. Feces were firstly kept at 4 °C in a 5 mL Eppendorf tube with 20% glycerol/Lysogeny Broth and then stored at –80 °C. Total DNA was extracted from feces as previously described [16], with a modification: a thermic shock of 30 seconds was performed between each bead-shaking step (3 bead-shaking steps of 30 seconds each at maximum speed). The 16S bacterial DNA V3–V4 regions were targeted by 357wf-785R primers and analyzed by MiSeq (RTLGenomics, Texas, USA). An average of 68,669 sequences was generated per sample. Bioinformatic filters applied as already described [17]. Cladogram and LDA scores were drawn via the LEfSe algorithm [18]. Diversity indices were calculated using the software Past 4.02 (Hammer, Ø., Harper, D.A.T., and P. D. Ryan, 2001. PAST: Paleontological Statistics Software Package for Education and Data Analysis. Palaeontologia Electronica 4(1): 9pp). The predictive functional analysis of the gut microbiota was performed via PICRUSt [19]. Diseases and host genetic variation linked to NEC-1_21–30 d-associated gut microbiota were identified via MicrobiomeAnalyst [20], with the Taxon Set Enrichment Analysis module.

## 2.3. Fecal Metabolome Analysis

The metabolome (total metabolites) analysis of the feces was performed as previously described [17]. Fecal samples have been prepared as it follows: 50 mg of feces were homogenized for 30 seconds in 500 µL of a pH 7.0 phosphate buffer, prepared in D₂O. Then, the homogenate was chilled into ice for 1 minute and centrifuged at 12,000 RPM for 10 minutes at 4 °C. The supernatant was then recovered, and the pellet was re-homogenized again at the same conditions of 12,000 RPM for 10 minutes at 4 °C. All the supernatants were then pooled and centrifuged at 18,000 RPM for 30 minutes at 4 °C. The supernatant was recovered and centrifuged again at the same conditions of 18,000 RPM for 30 minutes at 4 °C. A total of 600 µL of the final supernatant was then analyzed into nuclear magnetic resonance (NMR) tubes of 5 mm of diameter. The conformity criterium to validate the final sample was the aspect, to be crystal clear.

Pathway-associated metabolite sets and SNP-associated metabolite sets (Supplementary Figure S3 C,D, Supplementary Figure S4 B,C and Supplementary Figure S6 G) were analyzed via MetaboAnalyst 4.0 [21], with the enrichment analysis module.

## 2.4. Statistical Analysis

The results are presented as mean ± SEM for histograms and box and whiskers graphs. Statistical analyses were performed by two-way analysis of variance (ANOVA) followed by a two-stage linear step-up procedure of Benjamini, Krieger and Yekutieli to correct for multiple comparisons by controlling the false discovery rate (<0.05) (for histograms) or the Mann-Whitney test (for box and whiskers), as indicated in the figure legend, by using GraphPad Prism version 7.05 for Windows Vista (GraphPad Software, San Diego, CA, USA). For Table 1, results are presented as median or as indicated and *P* value was calculated using Fisher's exact test. Significant values were considered starting at *P* < 0.05. For the taxonomical and predictive functional analysis of gut microbiota, significant values were considered, starting at *P* < 0.05 or *P* < 0.01 when indicated. Principal component analysis (PCA) graphs were drawn by using Past 4.02.
<table><caption>**Table 1. Cohorts characteristics.**</caption><thead><tr><th>Variables, Description</th><th>NEC-1<br>n = 11</th><th>Healthy<br>n = 21</th><th>P (Fisher's Exact Test)</th></tr></thead><tbody><tr><td>Birth weight, median, inter-quartile (IQ)</td><td>1150 (845–1815)</td><td>1360 (700–2105)</td><td>Not significant (ns) (&gt; 0.05)</td></tr><tr><td>Gestational age, median (weeks) (IQ)</td><td>28.4 (26–31)</td><td>30 (26, 4–32)</td><td>ns</td></tr><tr><td>Gender</td><td></td><td></td><td></td></tr><tr><td>Girls, number (%)</td><td>2 (18)</td><td>8 (38)</td><td>ns</td></tr><tr><td>Boys, number</td><td>9</td><td>13</td><td></td></tr><tr><td>Patent Ductus arteriosus, number (%)</td><td>3 (27)</td><td>6 (28)</td><td>ns</td></tr><tr><td>Parity, number (%)</td><td>3 (27)</td><td>3 (10)</td><td>ns</td></tr><tr><td>Antenatal corticosteroids, number (%)</td><td>11 (100)</td><td>19 (90)</td><td>ns</td></tr><tr><td>Hypertension, eclampsia, number (%)</td><td>2 (18)</td><td>3 (14)</td><td>ns</td></tr><tr><td>Multiple births, number (%)</td><td>2 (18)</td><td>4 (19)</td><td>ns</td></tr><tr><td>Antenatal antibiotics, number (%)</td><td>4 (36)</td><td>5 (24)</td><td>ns</td></tr><tr><td>Chorioamniotitis, number (%)</td><td>2 (18)</td><td>1 (5)</td><td>ns</td></tr><tr><td>Apgar Score, 1 min (IQ)</td><td>8 (1–10)</td><td>7 (1–10)</td><td>ns</td></tr><tr><td>5 min (IQ)</td><td>10 (4–10)</td><td>8 (1–10)</td><td>ns</td></tr><tr><td>Cordon pH (IQ)</td><td>7.23 (6.8–7.4)</td><td>7.31 (7.04–7.43)</td><td>ns</td></tr><tr><td>Cordon lactates (IQ)</td><td>5.7 (2.9–9.6)</td><td>3 (1–7.2)</td><td>0.04 *<br>(&lt;0.05)</td></tr><tr><td>Mean arterial pressure at hospital<br>admission (IQ)</td><td>29 (20–47)</td><td>29.5 (17–43)</td><td>ns</td></tr><tr><td>Hospital Admission T (°C) (IQ)</td><td>36.5 (35.1–37.8)</td><td>36.8 (35.8–37.5)</td><td>ns</td></tr><tr><td>Antibiotics in the first week of life (%)</td><td>10 (90)</td><td>18 (85)</td><td>ns</td></tr><tr><td>Days under antibiotics (IQ)</td><td>7.5 (2–17)</td><td>3 (0–18)</td><td>ns</td></tr><tr><td>Days under antibiotics<br>(Third-Generation-Cephalosporin (3GC) ±<br>Penicillin A, ±aminoglycoside) in the first<br>week of life (IQ)</td><td>3 (0–7)</td><td>3 (0–7)</td><td>ns</td></tr><tr><td>Children under glycopeptides number (%)</td><td>8 (72)</td><td>3 (14)</td><td>0.0018 **<br>(&lt;0.01)</td></tr><tr><td>Bacteremia, number (%)</td><td>5 (45)</td><td>1 (5)</td><td>0.01 *<br>(&lt;0.05)</td></tr><tr><td>Exposition to mother milk, number (%)</td><td>11 (100)</td><td>21 (100)</td><td>ns</td></tr><tr><td>Age of enteropathy (days) (IQ)</td><td>12 (4–60)</td><td>-</td><td>-</td></tr><tr><td>Exposition to inotropes, number (%)</td><td>1 (9)</td><td>0 (0)</td><td>ns</td></tr><tr><td>Blood transfusion, number (%)</td><td>2 (18)</td><td>4 (19)</td><td>ns</td></tr><tr><td>Full enteral feeding (days) (IQ)</td><td>23 (14–39)</td><td>11 (3–29)</td><td>0.0002 ***<br>(&lt;0.001)</td></tr><tr><td>Median fasting time (median in days) (IQ)</td><td>1 (1–4)</td><td>-</td><td>-</td></tr></tbody></table>

## 3. Results

### 3.1. *Analysis of Gut Microbiota, Microbiome and Fecal Metabolome During NEC-1*

To understand the microbial and metabolomic evolution during the early onset of necrotizing enterocolitis (NEC), we studied clinical profile suspected NEC (NEC-1) preterm infants. NEC-1 children underwent more glycopeptides treatment, showed significantly higher cordon lactates, bacteremia and a longer full enteral feeding, when compared to age-matched healthy children (Table 1). NEC-1 children also displayed a lower plasma pH and enteral milk volume at day 7 (Supplementary Figure S1A,B) and a higher abundance of *Streptococcus* species (Supplementary Figure S2A) compared to healthy children. Both populations of children showed a high intragroup variance in terms of gut microbiota (Supplementary Figure S2B) and overall microbial diversity (Supplementary Figure S2C). NEC-1 microbiome showed increased activity for pathway related to transcription, glycosaminoglycan degradation and lysosome, compared to healthy children (Supplementary Figure S2D). Then, we analyzed the fecal metabolome to appreciate NEC-1-induced changes in gut microbial metabolic activity. NEC-1 children displayed a reduced intragroup variation and significantly lower levels of ethanol (Supplementary Figure S2E). Overall, these data show that NEC-1 is characterized by a precise gut microbiota, microbiome and gut microbial metabolites profile.
## 3.2. Analysis of Gut Microbiota, Microbiome and Fecal Metabolome During the Evolution of NEC-1 over the First Two Months of Life

Given the presence of a NEC-1-specific gut microbiota and microbiome profile, we aimed to identify at what time these profiles establish. We divided both NEC-1 and healthy children populations in subgroups according to periods of ten days of life as it follows: 1–10 d (d stands for “days”), 11–20 d, 21–30 d for the first month of life and > 30 d for the second one. In the first 10 days, NEC-1 children displayed a divergent and more homogenous gut microbiota compared to healthy children, with the latter characterized by a higher abundance of *Klebsiella* species (Figure 1A,B). At this stage of life, gut microbiota in NEC-1 had a lower diversity based on Chao-1 index (Figure 1C) and a different microbial activity related to replication, recombination and repair proteins, lysosome and glycosaminoglycan degradation (Figure 1D). No significant changes were observed in fecal metabolites (Figure 1E). Overall, these data show that gut microbiome starts to diverge at the early onset of NEC-1.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Analysis of gut microbiota, microbiome and metabolome in the first 10-days of life in healthy vs. necrotizing enterocolitis (NEC)-1 children. (A) Gut microbiota analysis via linear discriminant analysis (LDA) score between healthy (H) vs. NEC-1 children, in the first 10-days of life, 1 to 10 days; (B) principal component analysis (PCA) of the gut microbiota; (C) indices of gut microbiota diversity; (D) LDA score for microbial pathways; (E) histogram of the overall fecal metabolites and PCA as inset.

\*\*P < 0.01. Two-way ANOVA, followed by a two-stage linear step-up procedure of Benjamini, Krieger and Yekutieli to correct for multiple comparisons, by controlling the false discovery rate (<0.05); N = 15 for H and N = 4 for NEC-1.</figcaption></figure>

In the second 10-days of life, NEC-1 gut microbiota was characterized again by a higher abundance of *Streptococcus* species and bacteria from the Micrococcales order (Figure 2A), with a high intragroup variance (Figure 2B). At this stage of life, NEC-1 gut microbiota also showed a higher diversity based on Chao-1 index (Figure 2C), but no microbial pathway differently regulated (Figure 2D). As for the
fecal metabolome, NEC-1 children displayed significant lower levels of serine (Figure 2E). Overall, these data show a stronger evolution of gut microbiota than gut microbiome in the second 10-days of life, between NEC-1 and healthy children.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Analysis of gut microbiota, microbiome and metabolome in the second 10-days of life in healthy vs. NEC-1 children. (A) Gut microbiota analysis via LDA score between healthy (H) vs. NEC-1 children, in the second 10-days of life, 11 to 20 days (d) (the LDA score is only shown for NEC-1 children meaning that no bacteria are significantly higher in the H group vs. NEC-1); (B) PCA of the gut microbiota; (C) indices of gut microbiota diversity; (D) null cladogram for microbial pathways; (E) histogram of the overall fecal metabolites and PCA as inset. \*\*P < 0.01. \*\*\*P < 0.001. Two-way ANOVA, followed by a two-stage linear step-up procedure of Benjamini, Krieger and Yekutieli to correct for multiple comparisons by controlling the false discovery rate (<0.05); N = 14 for H and N = 10 for NEC-1.</figcaption></figure>

In the third 10-days of life, changes in NEC-1 gut microbiota compared to healthy children occurred to a bigger extent and were related to increased *Staphylococcus* and *Streptococcus* species (Figure 3A,B), together with a high intragroup variance (Figure 3C) and no change in the overall diversity indices (Figure 3D). We also observed a NEC-1 microbiome profile, mainly based on thiamine and seleno-compound metabolism (Figure 3E). The NEC-1 gut microbiota profile of the third 10-days of life was associated with: (i) multiple diseases and found significantly increased in ulcerative colitis (Figure 4A); (ii) host genetic variation and significantly related to ANP32E, a gene involved in ulcerative colitis [22], in line with previous reports. In terms of fecal metabolome, we observed no significant changes in NEC-1 vs. healthy children (Figure 4C). Then, we studied feces collected in the second month of life. In this period of life, the taxonomical differences in the gut microbiota of NEC-1 vs. healthy children were related to the increase in *Raoultella* species in NEC-1 gut microbiota (Figure 5A), with a still high intragroup variance (Figure 5B), and no change in the overall microbial diversity indices (Figure 5C). We also observed microbial functions related to DNA repair increased in the NEC-1
gut microbiome (Figure 5D). This period of life was characterized by the highest separation in terms of fecal metabolome, with significant lower levels of ethanol and leucine in NEC-1 children (Figure 5E).

<figure><img src="image_6.png" /><figcaption>**Figure 3.** A specific gut microbiota and microbiome exist in the third 10-days of life in healthy vs. NEC-1 children. (A) Comparative analysis of the gut microbiota by LDA effect size (LEfSe): the cladogram shows bacterial taxa significantly higher in the group of children of the same color, in the fecal microbiota between healthy (H) vs. NEC-1 children, in the third 10-days of life, 21 to 30 days (D) (the cladogram shows the taxonomic levels represented by rings with phyla at the innermost and genera at the outermost ring and each circle is a bacterial member within that level); (B) LDA score used to build the cladogram in (A); (C) PCA of the gut microbiota; (D) indices of gut microbiota diversity; (E) LDA score for microbial pathways. *N* = 13 for H and *N* = 7 for NEC-1.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Diseases, host genetic variation and metabolome analysis in the third 10-days of life during NEC-1. (A) Diseases and (B) host genetic variation linked to NEC-1_21-30d associated gut microbiota; (C) histogram of the overall fecal metabolites and PCA, as inset. *N* = 13 for H and *N* = 7 for NEC-1.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Analysis of gut microbiota, microbiome and metabolome in the second month of life in healthy vs. NEC-1 children. (A) Gut microbiota analysis via LDA score between healthy (H) vs. NEC-1 children, in the second month of life > 30 days (d); (B) principal component analysis (PCA) of the gut microbiota; (C) indices of gut microbiota diversity; (D) LDA score for predictive microbial pathways (*P* < 0.01); (E) histogram of the overall fecal metabolites and PCA, as inset. ***P* < 0.001. Two-way ANOVA, followed by a two-stage linear step-up procedure of Benjamini, Krieger and Yekutieli, to correct for multiple comparisons by controlling the false discovery rate (<0.05); *N* = 11 for H and *N* = 6 for NEC-1.</figcaption></figure>

### 3.3. Specific Impact of NEC-1 on the Evolution of Gut Microbiota. Microbiome and Fecal Metabolome over the First Two Months of Life, Compared to Healthy Children

To investigate the evolution of gut microbiota, microbiome and fecal metabolome over the first two months of life, we conducted an intra-group study in both NEC-1 and healthy children, according to the four groups reported above: 1–10 d, 11–20 d, 21–30 d and > 30 d. We did not observe any taxonomic significant change in the gut microbiota of NEC-1 children. However, the group NEC-1_21–30 d had a specific gut microbiome with an increased restriction enzyme activity, among others (Supplementary Figure S3A). The four NEC-1 groups also differed in terms of fecal metabolome, with regard to leucine, ethanol and serine amounts (Supplementary Figure S3B). Based on these results, we performed a metabolomic enrichment analysis on two levels: (i) pathway-associated metabolite sets (Supplementary Figure S3C) and (ii) single nucleotide polymorphism (SNP)-associated metabolite sets (Supplementary Figure S3D). NEC-1 metabolomic profile (increased ethanol and serine) was significantly associated with both homocysteine degradation and phosphatidylethanolamine biosynthesis (Supplementary Figure S3C), with serine being the metabolite the most linked to NEC-1-associated SNP (Supplementary Figure S3D). By contrast, in healthy children, the four groups reported above did not differ in terms of both gut microbiota and microbiome, but only with regard to fecal metabolome (Supplementary Figure S4A). Healthy metabolomic profile (increased leucine, ethanol and dihydroxyacetone) was significantly associated with valine, leucine and isoleucine degradation and to ketone body metabolism (Supplementary Figure S4B), with leucine being the metabolite the most linked to healthy-associated SNP (Supplementary Figure S4C). Overall, these data suggest that: (i) a different intragroup evolution
exists between NEC-1 and healthy children, with regard to gut microbiota and microbiome; and (ii) the NEC-1 microbiome appears to be more sensitive to mother-related factors.

### 3.4. Maternal and Child Factors Influencing the Gut Microbiota, Microbiome and Fecal Metabolome During NEC-1

Next, we asked which factor related to both mother and child may affect the most the above reported parameters. We analyzed six conditions: neonatal antibiotherapy (ABx), ABx treatment on mother, childbirth (C-section (C-sec) vs. vaginal birth (VB)), very low birth weight (VLBW), extreme low birth weight (ELBW) and gestational age (GA) > or ≤ 28 weeks.

Only neonatal ABx treatment affected the gut microbiota in both NEC-1 and healthy children (Supplementary Figure S5A). By contrast, all the above factors, except the VLBW, affected the gut microbiome (Supplementary Figure S5B–F). Note that childbirth modality, ELBW and GA affected the gut microbiome only in NEC-1 children (Supplementary Figure S5D–F). Moreover, all the above factors, except the neonatal ABx treatment and ELBW, affected the fecal metabolome between NEC-1 and healthy children (Supplementary Figure S6A–F). Then, we performed a metabolomic enrichment analysis on the pathway-associated metabolite sets, based on Supplementary Figure S6F, in which there is an increase in ethanol and succinate within in the NEC-1_GA ≤ 28 w. Ketone body and butyrate metabolism were the most significantly associated with this metabolomic set (Supplementary Figure S6G).

## 4. Discussion

In this prospective study, we focused on suspected necrotizing enterocolitis NEC-1 preterm infants. NEC-1 phenotype has been poorly clinically investigated, with no data available on gut microbiota, microbiome and fecal metabolome. By contrast, NEC-2 and NEC-3, more severe and established phenotypes, have been more characterized. As for clinical parameters, the increased cordon lactate levels we found in NEC-1 has been recently positively correlated to the development of enteropathy [23]. Hence, the hypothesis of hypoxic lesions in utero or during birth may not be excluded, and could even be predictive of neonatal morbidity. Importantly, the observed reduced enteral nutrition volume in NEC-1 is not a protective factor during the development of severe NEC (NEC-2 and NEC-3), but rather, it may lengthen hospitalization and infections risk [24]. Moreover, a recent multicentric study showed that a slow rate enteral feeding is associated with an increased risk of developing NEC-2 and NEC-3 [14].

NEC-1 children showed a high general variance for gut microbiota and fecal metabolome, which is in line with a personalized microbiota and fecal metabolome profiles of preterm infant [25]. This is evident from PCA analyses in Figure 4 C and Figure 5 E, where the metabolome profiles of the healthy vs. NEC-1 children appear to diverge, but are still presenting some overlap, due to high intragroup variance. Both this evidence and the delayed intestinal colonization of preterm infants [26,27] may explain the lack of NEC-1-specific microbial group in the first 10-days of life. The analysis by periods of ten days of life revealed a divergence for both gut microbiota and microbiome in NEC-1 by the third 10-days of life. In particular, the higher abundance of *Staphylococcus* in NEC-1 is in accordance with the early colonization by *Staphylococcus* bacteria of the intestine of preterm infants [28]. This datum suggests the third 10-days as an optimal time window to be targeted by antibiotics directed against bacterial species higher in NEC-1, such as *Staphylococcus*. However, in our study NEC-1 children who underwent glycopeptide and aminoglycoside therapy were more numerous than healthy children. Therefore, this evidence suggests that NEC-1 may be associated with glycopeptide and/or aminoglycoside-resistance, since NEC-1 gut microbiota was characterized by an increase, and not a decrease, of *Staphylococcus*. Since aminoglycosides are active antibiotics against enterobacteria, their administration could delay intestinal colonization by Proteobacteria, and thus promote the implantation of resistant genera, such as *Staphylococcus* and *Streptococcus*. Based on this evidence, our data suggest not to prolong antibiotic therapy beyond the first week of life in preterm infants, if inflammatory parameters have
been normalized. Furthermore, NEC-1 gut microbiota profile was associated with ulcerative colitis and host genetic variation in the ANP32E gene, encoding a protein implicated in cortico-resistance during ulcerative colitis [22]. NEC-1 children showed increased exposition to antenatal corticosteroids compared to healthy children, even though a study has not identified antenatal corticosteroids as a NEC-inducing factor [29]. Despite Anp32e-deficient mice display no sign of disease [30], it has not to be excluded the role of Anp32e in a model of gut inflammation mimicking ulcerative colitis. Hence, further studies are warranted on genetic factors of NEC. In terms of microbial functions, the intragroup analysis showed in the third 10-days of life a higher restriction enzyme activity in the NEC-1 gut microbiome. This bacterial activity, directed against bacteriophages and enriched in the newborn intestine [31], suggests an increased virus activity and, hence, a virome dysbiosis, beyond a microbiota dysbiosis, during NEC-1 evolution. Notably, with regard to the time-point > 30 days, microbiota and microbiome changes are still present and, most importantly, different from those observed at earlier time-points. This suggests that NEC-1 onset may have influenced the evolution of gut microbiota. All these microbial data are associated with our observation about a change in fecal amino-acids, such as leucine and serine, confirming the association between gut microbiota dysbiosis and a change in amino-acids metabolism [32].

Finally, several factors may affect the gut microbiota of preterm infants, e.g., hospital regimens, even with regard to the putative use of probiotics to blunt NEC. In a very recent publication, Kurath-Koller et al. showed the efficacy of some probiotics in improving the gut microbiome of very low birth weight infants during the first two weeks of age, in a triple-center cohort study [33]. In our mono-centric study, the very low birth weight had an impact on the fecal metabolome of both healthy and NEC-1 children, but not on their gut microbiota nor on their microbiome. This evidence underlines the importance of the hospital environment for the evolution of both an eubiotic and a dysbiotic, e.g., during NEC-1, gut microbiota.

## 5. Conclusions

The small cohort represents a limitation of our work, due to the difficulty to investigate NEC at its early onset. However, our data may provide neonatal departments with immediate indications to blunt NEC-1 evolution such as: i) increasing the enteral volume of nutrition, especially in the first days of life; ii) revising and reducing antibiotic therapy up to the first week of life in preterm infants.

**Availability of data and materials:** All data are available in the main text or the supplementary materials and via the following repositories: Sequence Read Archive (SRA) database with the assigned identifier PRJNA579480.

**Ethics approval and consent to participate:** This study was approved (number of the approval: DC 2016-2804) by Neonatal and Pediatric Intensive Care Unit and Neonatology Department of Purpan Hospital in Toulouse, France. The parents of the children involved in this study gave their approval by written consensus.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/2077-0383/9/7/2278/s1, Figure S1: Baseline plasma characteristics in healthy vs. NEC-1 children, Figure S2: Analysis of gut microbiota, microbiome and metabolome during NEC-1 over the first two months of life, Figure S3: A specific microbiome and metabolome exist in healthy vs. NEC-1 children over the first two months of life, Figure S4: Fecal metabolome progression over the first two months of life in healthy children, Figure S5: Maternal and child factors shaping gut microbiota and microbiome in healthy vs. NEC-1 children, Figure S6: Maternal and child factors shaping fecal metabolome in healthy vs. NEC-1 children.

**Author Contributions:** Conceptualization, C.B. and M.S.; Data curation, C.B. and M.S.; Formal analysis, C.B. and M.S.; Funding acquisition, M.S.; Investigation, M.S.; Methodology, C.B., D.D. and M.S.; Project administration, M.S.; Resources, C.B., O.D. and S.B.; Supervision, M.S. and E.O.; Validation, M.S.; Writing—original draft, C.B. and M.S.; Writing—review & editing, C.B. and M.S. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by LONGEVIE SA, grant to M.S.
**Acknowledgments:** We thank the departments of Neonatal and Pediatric Intensive Care Unit and Neonatology of the Purpan Hospital, Toulouse, France. We also thank Cecile Canlet and the Platform MetaToul-AXIOM for the metabolomic NMR analysis.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

## References

1. Bell, R.S. Neonatal necrotizing enterocolitis. *N. Engl. J. Med.* **1970**, *283*, 153–154. [CrossRef] [PubMed]

2. Bell, M.J. Neonatal necrotizing enterocolitis. *N. Engl. J. Med.* **1978**, *298*, 281–282. [CrossRef]

3. Kliegman, R.M.; Walsh, M.C. Neonatal necrotizing enterocolitis: Pathogenesis, classification, and spectrum of illness. *Curr. Probl. pediatr.* **1987**, *17*, 213–288. [CrossRef]

4. Ancel, P.Y.; Goffinet, F.; Group, E.-W.; Kuhn, P.; Langer, B.; Matis, J.; Hernandorena, X.; Chabanier, P.; Joly-Pedespan, L.; Lecomte, B.; et al. Survival and morbidity of preterm children born at 22 through 34 weeks' gestation in France in 2011: Results of the EPIPAGE-2 cohort study. *JAMA Pediatr.* **2015**, *169*, 230–238. [CrossRef] [PubMed]

5. Durazzo, M.; Ferro, A.; Gruden, G. Gastrointestinal Microbiota and Type 1 Diabetes Mellitus: The State of Art. *J. Clin. Med.* **2019**, *8*, 1843. [CrossRef] [PubMed]

6. Serino, M.; Blasco-Baque, V.; Nicolas, S.; Burcelin, R. Managing the manager: Gut microbes, stem cells and metabolism. *Diabetes Metab.* **2014**, *40*, 186–190. [CrossRef]

7. Serino, M.; Chabo, C.; Burcelin, R. Intestinal MicrobiOMICS to define health and disease in human and mice. *Curr. Pharm. Biotechnol.* **2012**, *13*, 746–758. [CrossRef]

8. Serino, M.; Luche, E.; Chabo, C.; Amar, J.; Burcelin, R. Intestinal microflora and metabolic diseases. *Diabetes Metab.* **2009**, *35*, 262–272. [CrossRef]

9. Bertani, L.; Caviglia, G.P.; Antonioli, L.; Pellicano, R.; Fagoonee, S.; Astegiano, M.; Saracco, G.M.; Bugianesi, E.; Blandizzi, C.; Costa, F.; et al. Serum Interleukin-6 and -8 as Predictors of Response to Vedolizumab in Inflammatory Bowel Diseases. *J. Clin. Med.* **2020**, *9*, 1323. [CrossRef]

10. Ni, J.; Wu, G.D.; Albenberg, L.; Tomov, V.T. Gut microbiota and IBD: Causation or correlation? *Nat. Rev. Gastroenterol. Hepatol.* **2017**, *14*, 573–584. [CrossRef]

11. Claud, E.C.; Walker, W.A. Bacterial colonization, probiotics, and necrotizing enterocolitis. *J. Clin. Gastroenterol.* **2008**, *42* (Suppl. S2), S46–S52. [CrossRef]

12. Pammi, M.; Cope, J.; Tarr, P.I.; Warner, B.B.; Morrow, A.L.; Mai, V.; Gregory, K.E.; Kroll, J.S.; McMurtry, V.; Ferris, M.J.; et al. Intestinal dysbiosis in preterm infants preceding necrotizing enterocolitis: A systematic review and meta-analysis. *Microbiome* **2017**, *5*, 31. [CrossRef] [PubMed]

13. Musemeche, C.A.; Kosloske, A.M.; Bartow, S.A.; Umland, E.T. Comparative effects of ischemia, bacteria, and substrate on the pathogenesis of intestinal necrosis. *J. Pediatr. Surg.* **1986**, *21*, 536–538. [CrossRef]

14. Roze, J.C.; Ancel, P.Y.; Lepage, P.; Martin-Marchand, L.; Al Nabhani, Z.; Delannoy, J.; Picaud, J.C.; Lapillonne, A.; Aires, J.; Durox, M.; et al. Nutritional strategies and gut microbiota composition as risk factors for necrotizing enterocolitis in very-preterm infants. *Am. J. Clin. Nutr.* **2017**, *106*, 821–830. [CrossRef]

15. Bazacliu, C.; Neu, J. Necrotizing Enterocolitis: Long Term Complications. *Curr. Pediatr. Rev.* **2019**, *15*, 115–124. [CrossRef] [PubMed]

16. Serino, M.; Luche, E.; Gres, S.; Baylac, A.; Berge, M.; Cenac, C.; Waget, A.; Klopp, P.; Iacovoni, J.; Klopp, C.; et al. Metabolic adaptation to a high-fat diet is associated with a change in the gut microbiota. *Gut* **2012**, *61*, 543–553. [CrossRef] [PubMed]

17. Nicolas, S.; Blasco-Baque, V.; Fournel, A.; Gilleron, J.; Klopp, P.; Waget, A.; Ceppo, F.; Marlin, A.; Padmanabhan, R.; Iacovoni, J.S.; et al. Transfer of dysbiotic gut microbiota has beneficial effects on host liver metabolism. *Mol. Syst. Biol.* **2017**, *13*, 921. [CrossRef]

18. Segata, N.; Izard, J.; Waldron, L.; Gevers, D.; Miropolsky, L.; Garrett, W.S.; Huttenhower, C. Metagenomic biomarker discovery and explanation. *Genome. Biol.* **2011**, *12*, R60. [CrossRef]

19. Langille, M.G.; Zaneveld, J.; Caporaso, J.G.; McDonald, D.; Knights, D.; Reyes, J.A.; Clemente, J.C.; Burkepile, D.E.; Vega Thurber, R.L.; Knight, R.; et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. *Nat. Biotechnol.* **2013**, *31*, 814–821. [CrossRef] [PubMed]
20. Dhariwal, A.; Chong, J.; Habib, S.; King, I.L.; Agellon, L.B.; Xia, J. *Microbiome Analyst*: A web-based tool for comprehensive statistical, visual and meta-analysis of microbiome data. *Nucleic Acids Res.* **2017**, *45*, W180–W188. [CrossRef]

21. Xia, J.; Psychogios, N.; Young, N.; Wishart, D.S. *MetaboAnalyst*: A web server for metabolomic data analysis and interpretation. *Nucleic Acids Res.* **2009**, *37*, W652–W660. [CrossRef] [PubMed]

22. Loren, V.; Garcia-Jaraquemada, A.; Naves, J.E.; Carmona, X.; Manosa, M.; Aransay, A.M.; Lavin, J.L.; Sanchez, I.; Cabre, E.; Manye, J.; et al. ANP32E, a Protein Involved in Steroid-Refractoriness in Ulcerative Colitis, Identified by a Systems Biology Approach. *J. Crohns Colitis* **2019**, *13*, 351–361. [CrossRef] [PubMed]

23. Tuuli, M.G.; Stout, M.J.; Shanks, A.; Odibo, A.O.; Macones, G.A.; Cahill, A.G. Umbilical cord arterial lactate compared with pH for predicting neonatal morbidity at term. *Obstet. Gynecol.* **2014**, *124*, 756–761. [CrossRef] [PubMed]

24. Oddie, S.J.; Young, L.; McGuire, W. Slow advancement of enteral feed volumes to prevent necrotising enterocolitis in very low birth weight infants. *Cochrane Database Syst. Rev.* **2017**, *8*, CD001241. [CrossRef]

25. Wandro, S.; Osborne, S.; Enriquez, C.; Bixby, C.; Arrieta, A.; Whiteson, K. The Microbiome and Metabolome of Preterm Infant Stool Are Personalized and Not Driven by Health Outcomes, Including Necrotizing Enterocolitis and Late-Onset Sepsis. *Msphere* **2018**, *3*. [CrossRef]

26. Groer, M.W.; Luciano, A.A.; Dishaw, L.J.; Ashmeade, T.L.; Miller, E.; Gilbert, J.A. Development of the preterm infant gut microbiome: A research priority. *Microbiome* **2014**, *2*, 38. [CrossRef]

27. Ho, T.T.B.; Groer, M.W.; Kane, B.; Yee, A.L.; Torres, B.A.; Gilbert, J.A.; Maheshwari, A. Dichotomous development of the gut microbiome in preterm infants. *Microbiome* **2018**, *6*, 157. [CrossRef]

28. Gibson, M.K.; Wang, B.; Ahmadi, S.; Burnham, C.A.; Tarr, P.I.; Warner, B.B.; Dantas, G. Developmental dynamics of the preterm infant gut microbiota and antibiotic resistome. *Nat. Microbiol.* **2016**, *1*, 16024. [CrossRef]

29. Travers, C.P.; Clark, R.H.; Spitzer, A.R.; Das, A.; Garite, T.J.; Carlo, W.A. Exposure to any antenatal corticosteroids and outcomes in preterm infants by gestational age: Prospective cohort study. *BMJ* **2017**, *356*, j1039. [CrossRef]

30. Reilly, P.T.; Afzal, S.; Wakeham, A.; Haight, J.; You-Ten, A.; Zaugg, K.; Dembowy, J.; Young, A.; Mak, T.W. Generation and characterization of the Anp32e-deficient mouse. *PLoS ONE* **2010**, *5*, e13597. [CrossRef]

31. Lim, E.S.; Zhou, Y.; Zhao, G.; Bauer, I.K.; Droit, L.; Ndao, I.M.; Warner, B.B.; Tarr, P.I.; Wang, D.; Holtz, L.R. Early life dynamics of the human gut virome and bacterial microbiome in infants. *Nat. Med.* **2015**, *21*, 1228–1234. [CrossRef] [PubMed]

32. Mardinoglu, A.; Shoaie, S.; Bergentall, M.; Ghaffari, P.; Zhang, C.; Larsson, E.; Backhed, F.; Nielsen, J. The gut microbiota modulates host amino acid and glutathione metabolism in mice. *Mol. Syst. Biol.* **2015**, *11*, 834. [CrossRef] [PubMed]

33. Kurath-Koller, S.; Neumann, C.; Moissl-Eichinger, C.; Kraschl, R.; Kanduth, C.; Hopfer, B.; Pausan, M.R.; Urlesberger, B.; Resch, B. Hospital Regimens Including Probiotics Guide the Individual Development of the Gut Microbiome of Very Low Birth Weight Infants in the First Two Weeks of Life. *Nutrients* **2020**, *12*, 1256. [CrossRef] [PubMed]

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