2.1 Basic information of the subjects
Among the 22 enrolled preterm infants, there were 14 males and 8 females. Regarding the gestational age, there were 8 cases with 28 weeks <= gestational age < 30 weeks, 6 cases with 30 weeks <= gestational age < 31 weeks and 8 cases with 31 weeks <= gestational age < 32 weeks. There was no significant difference in the gender ratio, delivery mode, feeding mode or gestational age stratification between the EUGR and control groups (Table 1).
2.2 Sequencing results and species accumulation curve
DNA extraction and PCR were successfully performed in all the stool samples, and the MiSeq libraries were prepared and sequenced. The species accumulation curve showed that the number of OTUs rapidly increased with the increase of the stool samples at the beginning, then slowly increased and entered a stationary phase, which proved that the number of samples in this study was sufficient to reflect the species richness in the community.(Fig.1).
2.3 OTU analysis
On the 14th day after birth, 8 infants had growth retardation, and a total of 1288 OTUs were obtained. Among the obtained OTUs, 331 were specific to the control group, 159 were specific to the EUGR group, and 798 were shared by the two groups. On the 28th day after birth, 10 infants developed growth retardation, and a total of 1177 OTUs were obtained. These OTUs included 238 ones specific to the control group, 226 ones specific to the EUGR group, and 713 OTUs shared by the two groups (Fig.2).
2.4 Alpha diversity analysis
Both the species richness (observed_species) and Shannon index (shannon) showed a stable pattern of the curve, which indicated that the sequencing depth was sufficient to reflect the vast majority of the microbial diversity in the samples, and an increased amount of data could result in limited new species. At the 2nd week after birth, there was an obvious difference between the Shannon index of the EUGR group and that of the control group, but this difference was insignificant. Meanwhile, there was no significant difference in any of the other indexes such as chao1, observed_species and PD_whole_tree between the two groups at 2 weeks or 4 weeks after birth (Fig. 3).
2.5 Beta diversity analysis
The principal component analysis (PCA) and non-metric multidimensional scaling (NMDS) analysis based on the OTU abundance in the stool samples showed that after 2 or 4 weeks of birth, the two groups had similar principal components of the intestinal microflora, but they could be roughly separated. This finding suggested that there was a certain difference in the structure of intestinal microflora between the two groups (Fig. 4).
2.6 LeFSe analysis
We performed the linear discriminant analysis effect size (LefSe) analysis, which is mainly used to identify species with a significantly different abundance among different groups. Based on the obtained species, an inter-group difference analysis was carried out, followed by a linear discriminant analysis (LDA) to estimate the effect of each species on the difference. At the 2nd week after birth, the items with an absolute LDA higher than 4 in the EUGR group included Streptococcaceae, Streptococcus, Bacteroidetes, Bacteroidales and Stenotrophomonas, while those in the control group included Enterococcaceae and Enterococcus. On the other hand, at the 4th week after birth, items with an absolute LDA higher than 3 in the EUGR group included Clostriaceae, Eubacteriaceae and Eubacterium (Fig.5).
2.7 Comparison of the community composition
On the 14th day after birth, there was a significant difference between the community principal components of the infants in the EUGR group and those of the normal control group on the family and genus levels. On the family level, the principal components in the control group were Enterococcaceae (53.95%) and Enterobacteriaceae (24.96%), while the Enterococcaceae content in the EUGR group significantly decreased to 23.45%, and there was no significant difference in the Enterobacteriaceae content (24.05%). Streptococcaceae accounted for 11.33% in the EUGR group and 2.58% in the control group, with no significant difference. On the genus level, the principal component of the control group was Enterococcus (53.92%), which accounted only for 26.44% in the EUGR group. The proportion of Streptococcus significantly differed between the two groups, accounting for 2.57% in the control group and 11.33% in the EUGR group. On the 28th day after birth, there was no significant difference between the community composition of the infants in the EUGR and control groups.