Microbial Composition in ASD and HC
In our study, a total of 15 fecal samples from children were collected by us, 9 ASDs and 6 HCs. The basic features of the 15 children included in our study are shown in table 1. The microbial composition at the Phylum level (relative abundance > 5% in the 15 samples) is shown in Figure 1A-1. Firmicutes was the most abundant genus in all samples. The microbial composition at the Class level (relative abundance > 1% in the 15 samples) is shown in Figure 1A-2. The level of Bacteroidales in the ASD group (29.291±8.689) was significantly less than that in the HC group (49.310±17.384) (t=2.963, p = 0.0110). The abundance of Selenomonadales in the ASD group (0.069±0.031) was significantly less than that in the HC group (7.423±3.186) (t= 3.157, p = 0.0076) (Figure 1B-1). The microbiota at the Order level is shown in Figure 1A-3. The average abundance was 1.0%. The abundance of 31 species was above the average in all 15 samples (Figure 1A-4). The average abundance of Ruminococcaceae was higher in the ASD group (25.378±9.562) compared to that in the HC group (13.561±8.313), (t=2.463, p = 0.0285). The average abundance of Prevotellaceae (5.540±8.836) in the ASD group was lower than that in the HC group (29.313±21.913), (t=2.957, p = 0.0111), as show in Figure 1B-2. The microbial composition at the Family level (relative abundance > 1% in the 15 samples) is shown in Figure 1A-4. Bacteroides and Prevotella_9 were the most abundant genus across all samples.
Microbiota and Clinical Features of ASD children
Table 2 shows the baseline data associated to the ASD and HC children enrolled in this study. The correlation between ASD and HC and gender, cesarean section and breastfeeding was further analyzed. A significant difference in gender, cesarean section and breastfeeding was not obtained between the two groups because of the small sample size. However, based on the differential analysis of gut microbiota of ASD and HC group, striking results were obtained. As regard the principal component analysis, ASD and HC children are distinguished in different quadrants, and they have distinct gut microbiota composition (Figure 2-A). The same result was obtained with the Tags Sampled analysis. The Chao1 index of ASD and HC group also showed the characteristics of two different biomes (Figure 2-B). In other alpha diversity analyses, the analysis of the number of OTUs and the Shannon Index also reflected the differences between ASD and HC children (Figure 2-C and Figure 2-D, respectively).
Predictive Microbiota Functional Profiling
The functional contribution of the bacteria in the ASD and HC group was predicted based on OTUs using the Tax4Fun package in R software. A total of 53 KEGG orthology were found across all samples, mainly belonging to the pathway Organismal Systems, Environmental Information Processing, and Cellular processes. The most abundant functional pathways are presented in Figure 3. Cluster analysis of the KEGG pathway revealed that the gut microbes differentially expressed in the ASD group are mostly representing the Organismal Systems pathway and the difference between the two groups was mainly in the nervous system pathway. Similarly, the Environmental Information pathway was mostly represented in the ADS group and was also significantly lower in the ADS group than in the HC group.