Introduction
Extrauterine growth restriction (EUGR) is defined as growth values less than the 10th percentile of intrauterine growth expectation in preterm neonates at the time of discharge.Reference Clark, Thomas and Peabody1 Although enteral nutrient supply for preterm or low birth weight neonates might support their growth at a rate similar to those of normal intrauterine growth within the first 2–3 months after birth, EUGR frequently prevails in very low birth weight (VLBW) (<1500 g) and extremely low birth weight (<1000 g) infants,Reference Radmacher, Looney, Rafail and Adamkin2, Reference De Curtis and Rigo3 indicating that low birth weight is a major risk factor for EUGR.Reference Lee, Kim and Namgung4, Reference Goulart, Morais and Kopelman5 In the follow-up, preterm infants with EUGR without interventions for 6–12 months after birth were found to suffer from slower cognitive and communicative development, limitations of movement (fetal head, sitting and walking), and 5.5–6.5% of them presented with transient disorders in spatial orientation and body position transformation.Reference Wright, Dawson and Fallis6 Since the physical and mental development of preterm infants is significantly affected by EUGR in both the short-term and long-term, it is important to provide adequate and effective nutritional support for preterm infants with EUGR at three stages, including the first few weeks after birth (acute stage), the development stage and the post-discharge stage.Reference Bashir, Louie, Shi and Nagler-Anderson7–Reference Sakurai, Itabashi, Sato, Hibino and Mizuno11
While recent studies have mainly focused on the short- and long-term effects of metabolic functions and neurodevelopment of preterm infants with EUGR, previous studies have shown that the incidence of EUGR is related to factors such as birth weight, gestational age (GA) and the time required for catch-up growth.Reference Lima, Carvalho, Costa and Moreira12 Although it remains unclear whether there is a relationship between the incidence of EUGR and neonatal necrotizing enterocolitis, both are common in very preterm infants.13–Reference Stoll, Hansen and Bell18 Moreover, the intake of adequate nutrients in neonates plays a vital role in the development of cell signal transduction of the gastrointestinal system, nervous system and the immune system. Defects in energy accumulation and deficiency of microelements greatly affect normal growth and development, at both the cellular and histological levels. Several factors affect preterm neonates’ intake of nutrients, including gastrointestinal dysfunction, endocrine abnormalities and other complications, such as respiratory problems, abnormal immunological functions, central nervous system damage and difficulties in coordinating sucking and swallowing.Reference De Curtis and Rigo3, Reference Embleton, Pang and Cooke19, Reference Panigrahi20 Enteral support to meet the nutritional needs of premature neonates at the critical postnatal period should decrease the risk of contracting various conditions.Reference Jefferis, Power and Hertzman21, Reference Jie22 However, although previous studies have suggested that the intestinal flora affects the incidence of ulcerative colitis, asthma and inflammatory bowel disease,Reference Eriksson and Forsen23 few studies have been carried out on the intestinal flora of preterm infants. Understanding alterations in the intestinal flora of preterm infants may yield appropriate treatment options, and thereby promote the growth and development of preterm infants with EUGR.
There are approximately 100 trillion live microorganisms present in the human intestine, of which 99% are bacteria of 500–1000 species. The tremendous diversity of bacteria is roughly divided into three categories, including probiotics, pathogenic bacteria and neutral bacteria. Significant differences in the intestinal flora have been found between premature and term birth infants, along with different postnatal growth and development.24 The purpose of the current study was to investigate the classification of intestinal bacteria in preterm neonates with EUGR and compare it with that of preterm neonates without EUGR.
Patients and methods
Subjects
Our study was performed in accordance with the Declaration of Helsinki with regard to ethical principles for research involving human subjects. The study protocols were approved by the Research Ethics Committee of the East Campus of Shanghai Sixth People’s Hospital Affiliated to Shanghai Jiaotong University. Written, informed consent was obtained. A single-center cohort of 59 premature babies (<1000 g and/or <30 weeks GA) born in the Children’s Hospital of Fudan University between January 1, 2016 and December 1, 2017 was prospectively included in the current study. Infants were excluded if they had severe apnea, a congenital malformation and/or required surgery.
EUGR has been diagnosed as a weight 10% lower than that of infants in the same GA when they were discharged from hospital. The discharge standard has been that the infant is >36 weeks old, weighting more than 2 kg without any underlying disease needing hospital care and can be feed orally. As a result, a total of 23 preterm infants with definitely diagnosed EUGR and 36 preterm infants without EUGR were recruited as the study group (EUGR group) and the control group (non-EUGR), respectively.
Birth weight recovery time has been defined as the period of physiological weight loss after birth (generally no more than 10%) and the return to birth weight (approximately 7–10 days). Weight increasing time refers to the growth phase after the end of weight loss.
Stool sample preparation
Fresh stool samples were collected at the point when infants were discharged from hospital and immediately preserved in liquid nitrogen, after which they were transported to our laboratory in a mobile refrigerator and frozen at −80°C until DNA extraction within 12 weeks. A total of 200 mg of stool sample was used for DNA extraction each time.
DNA extraction
Total community genomic DNA extraction was performed using an E.Z.N.A. Soil DNA Kit (Omega Bio-tek, USA) following the manufacturer’s instruction. The level of DNA was measured using a Qubit 2.0 fluorometer (Life Technologies, USA) to ensure that adequate amounts of high-quality genomic DNA had been extracted. The extracted total DNA was preserved at −80°C.
16S rRNA gene amplification by PCR
Our target was the V3–V4 hypervariable region of the bacterial 16S rRNA gene. Polymerase chain reaction (PCR) was started immediately after the DNA was extracted. The 16S rRNA V3–V4 amplicon was amplified using KAPA HiFi HotStart ReadyMix (2×) (TaKaRa Bio Inc., Japan). Two universal bacterial 16S rRNA gene amplicon PCR primers (PAGE purified) were used: the amplicon PCR forward primer (CCTACGGGNGGCWGCAG) and the amplicon PCR reverse primer (GACTACHVGGGTATCTAATCC).
The reaction was set up as follows: Microbial DNA (10 ng/μl) 2 μl; amplicon PCR forward primer (10 μM) 1 μl; amplicon PCR reverse primer (10 μM) 1 μl; 2×KAPA HiFi HotStart ReadyMix 15 μl (total 30 μl). The plate was sealed and PCR was performed in a thermal instrument (Applied Biosystems 9700, USA) using the following program: one cycle of denaturing at 95°C for 3 min, five cycles of denaturing at 95°C for 30 s, annealing at 45°C for 30 s, elongation at 72°C for 30 s, then 20 cycles of denaturing at 95°C for 30 s, annealing at 55°C for 30 s, elongation at 72°C for 30 s with a final extension at 72°C for 5 min. The PCR products were checked using electrophoresis in 1% (w/v) agarose gel in Tris, boric acid, ethylenediaminetetraacetic acid (TBE) buffer stained with ethidium bromide and visualized under UV light.
16S gene library construction, quantification and high-throughput sequencing
AMPure XP beads were used to purify the free primers and primer dimer species in the amplicon products. Samples were delivered to Sangon BioTech (Shanghai) for library construction using a universal Illumina adapter and index. Before sequencing, the DNA concentration of each PCR product was determined using a Qubit® 2.0 Green double-stranded DNA assay and the quality was confirmed using a bioanalyzer (Agilent 2100, USA). Depending on coverage needs, all libraries were pooled for one run. The amplicons from each reaction mixture were pooled in equimolar ratios based on their concentrations. High-throughput sequencing was performed at Sangon Biotech in Shanghai using the MiSeq PE300 Sequencing System (Illumina, San Diego, CA, USA).
Statistical and bioinformatics analysis
After sequencing, data were collected as follows: The two short Illumina readings were assembled by PEAR software (version 0.9.6) according to the overlap, and the FASTQ files were processed to generate individual FASTA and QUAL files, which could then be analyzed by standard methods.
Sequences containing ambiguous bases or sequences longer than 480 base pairs (bp) were dislodged, with an allowance of a maximum homopolymer length if 6 bp.Reference Kochling, Sanz, Gavazza and Florencio 25 Sequences shorter than 200 bp were removed and any ambiguous base calls in the adapter and barcode sequences were deleted from the dataset.
The pre-cluster and chimeras.uchime commands of MothurReference Schloss, Westcott and Ryabin 26 were used to detect and remove chimera sequences. Distance matrices were established using the dist.seqs command, with the operational taxonomic units (OTUs) defined by using the furthest neighbor clustering algorithm at a phylogenetic similarity of 93–97%. Good’s coverage estimator was used to confirm the completeness of sampling. The rarefaction curves of OTUs, Good’s coverage and other richness and diversity indices of the bacterial community (i.e., ACE, Chao1, Shannon and Simpson) were estimated using the Mothur software.Reference Schloss, Westcott and Ryabin 26
Taxonomic classifications were conducted using the online Ribosomal Database Project (RDP) classifier, with a confidence threshold of 80%.Reference Wang, Garrity, Tiedje and Cole 27 A one-sample Kolmogorov–Smirnov (K–S) test was used to test the normality of the data. A general linear model (for normally distributed data) and a generalized linear model (for non-normally distributed data) were used to quantify the effects of host age and host species on the relative abundance of the top 5 phyla.
An independent sample t-test (for normally distributed data) or a Mann–Whitney U-test (for non-normally distributed data) was used to compare data between groups with the same classification.
A multivariate logistic model was employed to compare the positive rates of bacterial genera between the EUGR and control groups after adjusting for potential confounding factors, including birth weight, birth weight recovery days, start time of breast feeding, weight increase, parenteral nutrition days, antibiotics administration days and human milk or human milk plus thickener. A false-discovery rate (FDR) correction for multiple testing was applied for comparisons of bacterial genera to minimize type I errors, and only a corrected P-value of <0.05 was considered to be significant.
The Bray–Curtis similarity index was used as a metric of similarity between the bacterial communities based on the abundance of OTUs between samples. A heatmap analysis was conducted to compare the overall bacterial composition associated with the species and the age of the hosts. Venn diagrams and statistical clustering were used to determine the shared OTUs by all group members that were defined as core microbiome. The heatmap figures and Venn diagrams were produced using R1.
Results
Characteristics of preterm infants
A total of 59 premature infants (EUGR group, n=23; non-EUGR group, n=36) were eligible for our study. As shown in Table 1, no significant difference was found between the EUGR and non-EUGR control group in terms of GA, birth weight, birth weight recovery time, feeding age and weight increase time. However, preterm infants in the EUGR group had a significantly longer time of parenteral nutrition absorption and antibiotics administration (22±14 v. 17±10 days, P<0.05; 21±10 v. 14±9.5 days, P<0.05, respectively).
EUGR, extrauterine growth restriction.
In addition, the nurturing form and food composition for preterm infants in the two groups were also significantly different, as shown in Table 2. Significantly increased percentages of preterm infants with EUGR were nurtured with preterm formula and starch-based thickener, compared with the non-EUGR (13.0 v. 8.3%, P<0.05; 21.7 v. 0%, P<0.05, respectively), in contrast to significantly decreased percentages of human milk and carob-based thickener (60.9 v. 83.3%, P<0.05; 60.9 v. 83.3%, P<0.05, respectively).
EUGR, extrauterine growth restriction.
Table 3 summarizes and compares the abundance of 344 genera and 456 species of bacteria identified using high-throughput sequencing of all stool samples in the two groups. The total abundance of bacteria and the frequency of bacteria in terms of genus (7568 v. 13,760, P<0.000), and species (10,032 v. 18,240, P<0.000) were significantly decreased in the EUGR group compared with the non-EUGR group, while positive percentages of bacterial strains in the different groups were not different (10.0 v. 9.8%, P=0.746 and 8.2 v. 7.7%, P=0.187).
EUGR, extrauterine growth restriction.
* P<0.001, the abundance number of EUGR compared to number of control
Comparison of gastrointestinal bacterial abundance in terms of genus
Figure 1 shows heatmaps of bacterial genera (left side) and species (right side) abundances in stool samples of EUGR and control infant groups.
As shown in Fig. 2, there were 20 and 11 bacteria genera with significantly increased (upper panel) and decreased (lower panel) abundance in the EUGR group compared with the control group, respectively.
Comparison of gastrointestinal bacterial abundance in terms of species
As shown in Fig. 3, preterm infants with EUGR had significantly increased species of Streptococcus, Bacterium_RB5FF6 and Staphylococcus compared with the non-EUGR group. The abundance of Enterococcus was different due to the presence of different subgroups. Lactis was significantly decreased, while Faecium significantly increased in the EUGR group compared with the non-EUGR control group.
In addition, as shown in Fig. 4, some species of Pseudomonas (Fluorescens and Putida) significantly decreased, while others (Moraviensis, sp_Cra38 and sp_REmamp_189) significantly increased in the EUGR group compared with the non-EUGR group. The baumannii species of Acinetobacter significantly decreased, in contrast to the significantly increased Sp_Endosymbiont of Acinetobacter in the EUGR group. Similarly, preterm infants with EUGR showed significantly increased Aerogenes, Ludwigii and sp_FF3, in contrast to significantly decreased sp_CCBAU_15567 species of Enterobacter, compared with non-EUGR preterm infants.
Moreover, as shown in Fig. 5, there were significant changes in eight Klebsiella species, with seven species (Granulomatis, Michiganenis, Pneumoniae, sp_enrichment_culture, sp_ok1_1_9_814, sp_SCAUS856 and Variicola) significantly increased and one species (sp_enrichment_culture1) significantly decreased in the EUGR group. Furthermore, there were 11 other species that showed either significantly increased (Eubacterium_budayi, Marine_bacterium_BPYW9, Microbacterium_sp_TSWCW12, Propionibacterium_avidum_44067, Rothia_ mucilaginosa_ M508 and Unidentified marine_bacterioplan) or decreased (Bacillus_sp_3B1332, Bifidobacterium_animalis, Carnobacterium_maltaromaticum_ LM, Endosymbiont_of_Sphenophorus_lev and Pantoea_agglomerans) abundance in the EUGR group compared with the control group.
Comparison of gastrointestinal bacterial abundance in terms of genus and species with multivariate logistic analysis and FDR correction
After using a multivariate logistic model to compare the positive rates of bacterial genera and species between the EUGR and control groups, after adjusting for potential confounding factors including birth weight, birth weight recovery days, start times of breast feeding, weight increase, parenteral nutrition days, antibiotics administration days and human milk or human milk plus thickener and FDR correction, we found five genera with significantly increased or decreased abundance in the EUGR compared with the control group. In addition, the abundance of nine species differed significantly between the two groups, including three Klebsiella and one Enterobacter species, Streptococcus parasanguinis, Bacterium RB5FF6 and Endosymbiont_of_Sphenophorus_lev, which occurred more abundantly and one Microbacterium and one Acinetobacter species, which occurred less abundantly in the EGFR group samples (Table 4, Figs 2–5).
EUGR, extrauterine growth restriction.
Discussion
The incidence rate of EUGR is around 63% in our hospital. The main risk factors include premature birth and low birth weight, long time use of antibiotics, not breast-fed and contracting hospital infections more than once. In the current study, we found significant differences in the intestinal bacterial species between premature infants with and without EUGR. First, we observed that although having similar baseline characteristics, different feeding patterns, longer times of parenteral nutrition and antibiotic administration occurred in EUGR preterm infants, which resulted in significant differences in the total abundance of bacterial genus and species compared with the non-EUGR control group. There were significant differences in five genera and nine species (Table 4) after adjusting for confounding factors and FDR correction.
Interestingly, especially in terms of bacterial species, preterm infants with EUGR had significantly increased Staphylococcus compared with the non-EUGR group. In addition, significantly increased Enterococcus faecium was found in the EUGR infants without adjusting analysis, which is in line with a previous study that found four phases of microbiota development in newborns and the Enterococcus phase was a characteristic of prematurity.Reference Korpela, Blakstad and Moltu 28 However, after adjustments for cofounding factors and FDR corrections these differences disappeared.
Previous studies have also demonstrated that the growth of microflora could be affected by the birth state of the neonate, infections and exposure to antibiotics, which may result in irreversible alteration of bacterial species.Reference Antonopoulos, Huse and Morrison 29 , Reference Penders, Thijs and Vink 30 In addition, preterm infants might develop a unique neonatal intensive care unit flora, due to prolonged exposure.Reference Groer, Luciano and Dishaw 31 In the present study, before applying the adjustment analysis, Lactobacillus, Pseudomonas fluorescens and Pseudomonas putida were significantly decreased in the EUGR group. Especially delayed development of a Lactobacilli-rich microbiota has been associated with poor health in preterm infants.Reference Butel, Suau and Campeotto 32 The decreased Lactobacillus could have been caused by antibiotic administration to premature infants with EUGR or by EUGR itself, suggesting that the immunity and resistance of EUGR infants could be strengthened by supplementary Lactobacillus, and the low birth weight of EUGR neonates could be increased by promoting the absorption of protein, monosaccharides and other mineral nutrients, such as calcium and magnesium.Reference Sukanyaa, Vinoth and Ramesh 33 However, after correction analysis, we identified that these genera changes were all caused by confounding factors such antibiotic usage, human milk or lower birth weight.
Previous studies have shown that autoimmune and allergic diseases in the fetus and newborn were caused by immune system disorders, which resulted from intestinal flora maladjustment.Reference Huffnagle 34 This finding suggests that intestinal microflora disorders or abnormalities could affect the immunological functions of preterm infants, since the disappearance of some bacterial species is associated with a shortage of some immunocytes involved in the regulation of immune responses. Therefore, the alteration of intestinal bacterial species at the initial stage of the neonate is very likely to affect the development of the cellular and immune system by influencing the recognition of various antigens, which first appear at the mucosal barrier and then throughout the whole body. In the current study, we report for the first time, that preterm infants with EUGR had a significantly different abundance of intestinal bacterial genus and species compared with the non-EUGR premature infants, by analyzing high-throughput sequencing of the bacterial rRNA library from the total community genomic stool DNA. As suggested by other studies, the maladjustment of bacterial species in EUGR preterm infants is likely to result in susceptibility to food allergies through Toll-like receptor 4 signaling and promote T(H)1/T(H)2 immune deviation.Reference Bashir, Louie, Shi and Nagler-Anderson 7 , Reference Oyama, Sudo, Sogawa and Kubo 35
In addition, previous studies have demonstrated an increased presence of bio indicators such as C-reactive protein (CRP), human growth factor (HGF), interleukin-8 (IL-8), monocyte chemotactic protein 1 (MCP-1) and tumor necrosis factor α (TNF-α) in EUGR neonates, suggesting early changes in intestinal metabolism, nutrition absorption and immunological functions.Reference Ortiz-Espejo, Perez-Navero and Olza-Meneses 36 A concerning finding in the present study is the increasing prevalence over time of Streptococcus, Bacterium_RB5FF6, as well as Enterobacter and Klebsiella species in EUGR infants compared to healthy infants. It has been reported that Enterococcaceae were also composed of potential pathogens that have been frequently associated with late-onset sepsis in VLBW infants, and infrequently associated with necrotizing enterocolitis (Clostridiaceae).Reference de la Cochetiere, Piloquet and des Robert 37 However, it remains to be established what causes these changes in intestinal microflora in EUGR neonates. It is not known whether they are due to decreased immune functions and/or decreased metabolic functions of the EUGR neonates themselves. A prospective study with completely identical baseline characteristics in the affected and the control groups will be needed to investigate the possible causes of intestinal microflora changes in premature infants.
In summary, the current study has shown obvious changes in intestinal microflora between premature neonates with EUGR and without EUGR, with a significantly increased abundance of Klebsiella and Streptococcus strains as well as an Enterobacter species, while abundances of Microbacterium and Acinetobacter were significantly decreased in EUGR compared with non-EUGR preterm infants. However, we cannot unequivocally specify whether the abundance difference of these bacterial genera and species directly affects the growth of EUGR infants.
In conclusion, there was a significant difference in gut microbiota between EUGR and control children at discharge from the hospital, which might be attributed to different feeding or the use of antibiotics during hospitalization. Whether these different microorganisms in EUGR patients will be maintained and have an effect on their childhood growth and development will be the focus of a future study.
Acknowledgments
This study was supported by the National Natural Scientific Foundation [grant number: 81200456] as well as the Project of Chinese and Canadian Newborn Collaboration Network and the Shanghai University of Medicine & Health Science hundred teacher plan and Science and Technology Commission livelihood projects in Pudong District (PKJ2017-Y09).
Ethical Standards. Ethics approval and consent to participate
Our study was performed in accordance with the Declaration of Helsinki with regard to ethical principles for research involving human subjects. The protocols for our study were approved by the Research Ethics Committee of Children’s Hospital of Fudan University. Written, informed consent was obtained.
Data availability
The datasets supporting the conclusions of this article are included in the article.
Conflicts of interest
The authors declare no conflicts of interest.