Hui Li1†Wei Guo2†Sijie Li3Bishao Sun4Ningshan Li1Dongjing Xie5Zongming Dong4Dan Luo6
Wei Chen4Weihua Fu4Ji Zheng4*
Jingzhen Zhu4*- 1Department of Ultrasound, Xinqiao Hospital, Army Medical University, Chongqing, China
- 2Stroke Center, Puyang People's Hospital, Puyang, China
- 3Department of Pediatrics, Xinqiao Hospital, Army Medical University, Chongqing, China
- 4Department of Urology, Urologic Surgery Center, Xinqiao Hospital, Army Medical University, Chongqing, China
- 5Department of Neurology, Xinqiao Hospital, Army Medical University, Chongqing, China
- 6Department of Neurology, Yunyang People's Hospital, Yunyang, China
Background: Autism spectrum disorder (ASD) is associated with alterations in the gut microbiome. However, there are few studies on gut microbiota of children with ASD in China, and there is a lack of consensus on the changes of bacterial species.
Purpose: Autism spectrum disorder (ASD) is associated with alterations in the gut microbiome. However, there are few studies on gut microbiota of children with ASD in China, and there is a lack of consensus on the changes of bacterial species.
Methods: We used 16S rRNA sequencing to analyze ASD children (2 to 12 years), HC (2 to 12 years).
Results: Our findings showed that the α-diversity, composition, and relative abundance of gut microbiota in the ASD group were significantly different from those in the HC groups. Compared with the HC group, the α-diversity in the ASD group was significantly decreased. At the genus level, the relative abundance of g_Faecalibacterium, g_Blautia, g_Eubacterium_eligens_group, g_Parasutterella, g_Lachnospiraceae_NK4A136_group and g_Veillonella in ASD group was significantly increased than that in HC groups, while the relative abundance of g_Prevotella 9 and g_Agathobacter was significantly decreased than that in HC groups. In addition, KEGG pathway analysis showed that the microbial functional abnormalities in ASD patients were mainly concentrated in metabolic pathways related to fatty acid, amino acid metabolism and aromatic compound metabolism, and were partially involved in neurotransmitter metabolism.
Conclusion: This study revealed the characteristics of gut microbiota of Chinese children with ASD and provided further evidence of gut microbial dysbiosis in ASD.
Introduction
Autism spectrum disorder (ASD) refers to a group of complex neurodevelopmental disorders that typically present in early childhood and are characterized by stereotyped behavior, language, and social interaction disorders (Valentino et al., 2021). The prevalence of ASD has increased significantly in recent decades, affecting 1 in 59 children in the United States (Maenner et al., 2020) and approximately 39.23/10000 Chinese children aged 1.6–8 years (Wang et al., 2018). Patients with ASD require lifelong treatment and intervention, making it a serious public health concern (Joon et al., 2021). Both genetic and environmental factors are believed to be potential triggers for ASD (Wang et al., 2016; Modabbernia et al., 2017; Lord et al., 2018), with abnormal gut microbiota being an important environmental factor. Gastrointestinal problems are common in children with ASD, affecting 9–91% of cases (Ferguson et al., 2019) and often correlating with the severity of ASD (Adams et al., 2011).
Research has shown that the human gut microbiota plays a crucial role in influencing the brain through a communication pathway, called the gut-brain axis (Fowlie et al., 2018; Iannone et al., 2019; Lacorte et al., 2019; Sharon et al., 2019; Srikantha and Mohajeri, 2019). Changes in the gut microbiota can regulate gastrointestinal physiology, immune function, and behavior through the gut-microbiota-brain axis (Li et al., 2019; Sharon et al., 2019; Tomova et al., 2020; Nogay and Nahikian-Nelms, 2021). Germ-free animal studies have demonstrated that the gut microbiota shapes the brain by influencing synaptic signaling and gene transmission (Sanlier and Kocabas, 2021). Lactobacillus reuteri treatment has been shown to selectively reverse social deficits in mouse models of ASD, further highlighting the involvement of gut microbiota in the pathogenesis of ASD (Sgritta et al., 2019). While several human studies have shown alterations in the gut microflora in ASD, there is little consensus on the specific bacterial species involved. Conflicting results have been reported in different population studies, with some showing increased bacterial diversity (Wan et al., 2022), others showing no change (Kandeel et al., 2020), and still others showing a decrease (Dan et al., 2020) compared with healthy controls. Some studies have found that potentially harmful microbiota, such as Clostridium (Wang et al., 2019), Sutterellaceae (de Angelis et al., 2013), and Enterobacteriaceae (Ding et al., 2020), were more abundant in children with ASD, while Firmicutes (Liu et al., 2019) and Prevotellaceae (Pulikkan et al., 2018) were decreased in the ASD group. However, other studies have shown the opposite results (Rose et al., 2018; Zou et al., 2020; Ye et al., 2021). Studies have also found that Clostridium and Bacteroidetes/Firmicutes ratios were increased in ASD children with functional gastrointestinal disease (Luna et al., 2017), while a study in Slovakia showed increased numbers of Lactobacillus and a significant decrease in the Bacteroidetes/Firmicutes ratio (Tomova et al., 2015). Another study suggested that decreased Bifdobacteria abundance might lead to reduced folate production in individuals with ASD, so abnormal folate metabolism might be associated with ASD (Frye et al., 2017). In conclusion, while the results on gut microbiota and ASD are still inconsistent and mainly focused on Western populations, there is evidence to suggest that altered gut microbiota is associated with ASD. Further research is needed to better understand the role of gut microbiota in ASD and to identify non-invasive biomarkers for early diagnosis.
Given the inconsistent findings on the gut microbiota in ASD and the limited research on Chinese populations, it is important to investigate the changes in gut microbiota in Chinese children with ASD. Based on fecal 16S ribosomal DNA (rDNA) sequencing data from 957 children with ASD and 161 healthy controls in China, we aimed to determine the taxonomic composition of gut microbiota of children with ASD and the changes in the gut microbiota compared with healthy controls. Additionally, we identified potential biomarkers of bacterial species that could serve as non-invasive tools for early diagnosis of ASD.
Materials and methods
Ethics statement
This study was approved by the Ethics Committee of Second Affiliated Hospital, Army Medical University (Ethics NO. 2023–001-01). Written informed consents were signed by all the children’s legal guardians prior to the study. The clinical trial number: ChiCTR2300074832.
Study subject recruitment and fecal sample collection
In total, 1,118 participants from 12 provinces of China including 957 participants with clinical definition as ASD (aged between 2 and 12 years, average age 4.6) and 161 health children (HC) (aged between 2 and 12 years, average age 4.8) were recruited (Supplementary Table S1).
All children with ASD participating in this study were enrolled from social groups of unrelated families with autism (Participants had been diagnosed with ASD according to DSM-5 (Diagnostic and Statistical Manual of Mental Disorders - 5th Edition) before enrollment). The exclusion criteria included diseases such as depressive disorder, cerebral palsy, schizophrenia, bipolar disorder, significant sensory impairment, and clinically significant inflammatory conditions.
The control group of children were recruited from kindergartens and primary school. Children with psychiatric conditions (such as depressive disorder, schizophrenia, and bipolar disorder) were excluded on the basis of entrance examinations and parent interviews.
All subjects had not taken antibiotics, antipsychotics, probiotics, and prebiotics in the month prior to the sample collection. Feces samples were used to analyze the genes associated with microbiota. Supplementary Figure S1 presents the study flowchart.
Feces were collected by parents in hospitals or at home. The feces were collected according to the instructions and delivered immediately at a low temperature. The frozen feces were transported overnight with dry ice to CNNC Yili (Tianjin) Medical Laboratory Co., Ltd., where they were frozen at −80°C until DNA was extracted.
16S rRNA sequencing and analysis
DNA from intestinal contents was prepared from subjects in the ASD and HC groups using a fecal DNA extraction kit (DP712, Tiangen Company, Beijing, China). The analysis was repeated with three biological replicates in each group. DNA quality was monitored on 1% agarose gels. The hypervariable V3-V4 region of the 16S-rDNA gene was amplified by PCR using the primers 338F (5’-ACTCCTACGGGAGGCAGCAG-3′) and 806R (5’-GGACTACHVGGGTWTCTAAT-3′), where the barcode was an eight-base sequence unique to each sample. PCR was performed in 30-μL reactions composed of 15 μL of Phusion® High-Fidelity PCR Master Mix (New England Biolabs), 0.2 μM forward and reverse primers, and 10 ng of template DNA. The thermal cycle consisted of an initial denaturation step at 98°C for 1 min, followed by 30 cycles of denaturation at 98°C for 10 s, annealing at 50°C for 30 s, and elongation at 72°C for 30 s. Finally, the reactions were incubated at 72°C for 5 min. PCR products were detected on 2% agarose gels by electrophoresis and then purified using the GeneJET Gel Extraction Kit (Thermo Scientific). Sequencing libraries were generated using the Illumina TruSeq DNA PCR-Free Library Preparation Kit (Illumina, United States) according to the manufacturer’s recommendations, and index codes were added. Sequencing was performed using the Illumina NovaSeq MiSeq PE300 platform (Genecloud Co. Ltd., Chongqing, China). The sequence analysis was performed using the QIIME software package 2 (version 2020.2). The raw data of 16S rRNA sequencing underwent decomposition and quality filtering using version 0.20.0 of fastp, and merged using FLASH version 1.2.7 with the following merging conditions: (i) Any reads with an average quality score < 20 in a 50 bp sliding window were truncated at any position. Reads shorter than 50 bp were discarded, as well as reads containing ambiguous characters. (ii) Assembly was performed only based on overlapping sequence lengths greater than 10 bp, with a maximum mismatch ratio of 0.2 in the overlapping region. Reads that could not be assembled were discarded. (iii) Samples were differentiated based on barcodes and primers, and the sequence orientation was adjusted. Barcodes were matched exactly, while primer matching allowed for 2 nucleotide mismatches. Using the MOTHUR workflow, cluster effective tags with a similarity of ≥97% into operational taxonomic units (OTUs). The representative sequence within each cluster was selected as the tag sequence with the highest abundance. Bacterial alpha diversity was determined based on OTU analysis of samples and expressed as Chao index, Shannon index and Simpson index, which were calculated using the R program package “vegan.” Principal coordinate analysis (pCoA) was performed using the R package1 to represent the β diversity of the microbiome. Using the phyloseq package to calculate Bray-Curtis metric distances, UniFrac distances, and weighted UniFrac distances. The wilcoxon rank sum test was used to compare the differences of bacterial phyla, classes, orders, families and genera between the two groups. The linear discriminant analysis (LDA) effect size (LEfSe) method was used to analyze the bacterial community dominance between groups.2 LEfSe and Kruskal−Wallis rank sum test (p < 0.05) to identify features with significant differences in abundance levels between specified taxa, and to evaluate the effect size of each feature using LDA (LDA score (log10) =3.5 as the cutoff value). The R language mixOmics package was used to calculate Pearson correlation coefficient, and the correlation coefficient r2 and p value of differential gut microbiota and differential metabolites were calculated. The metabolic pathways of intestinal microorganisms were analyzed based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database.3 We evaluated the metabolic pathway enrichment of the two groups. The metabolic pathway was considered enriched, while when the p value of the metabolic pathway was <0.05. Correlation analysis was performed by calculating Pearson’s correlation coefficient using the R language mixOmics package to calculate the correlation coefficient r2 and the p value of differentially bacterial and differentially metabolites.
The raw data of 16S rRNA sequence of this paper have been deposited in the Genome sequence Archive (Genomics, Proteomics & Bioinformatics 2021) of the National Genomics Data Center (Nucleic Acids Res 2022) and the National Center for Biological Information of China/Beijing Institute of Genomics, Chinese Academy of Sciences (GSA - human: HRA004410), and can be in open access at https://ngdc.cncb.ac.cn/gsa-human.
Co-occurrence network analysis
To understand the correlations among different genera, we constructed co-occurrence network based on the 16S rRNA data (Frye et al., 2017; Wang et al., 2018). The bacterial correlations in the ASD and HC samples were analyzed, respectively, according to the relative abundance of each genus using Spearman’s correlation coefficient to construct the co-occurrence network. The significant correlated genus (p < 0.05, r ≥ 0.7) were visualized by Gephi 0.9.7 (www.netbeans.org). Then, the similarity between the two network structures was measured by node closeness and shared correlations. Closeness of the nodes was analyzed by Gephi to predicate node centralities in each network. The shared correlations between two groups were defined the edges with the same nodes in two co-occurrence networks. Only genera existed in at least 10% sample were included in the network analysis.
Statistical analysis
The data was presented as mean ± standard error (SEM). SPSS 26.0 software (SPSS, Inc., Chicago, IL, United States) was used to perform independent sample t test for the two groups of data. A two-tailed Wilcoxon rank-sum test was performed for the analysis of gut microbiota sequencing data by R Project. p < 0.05 was considered statistically significant.
Results
Differences of microbial diversity between ASD and healthy controls
The alpha diversity of gut microbiota was measured using various indices, including the Shannon index, Chao1 index, and Simpson index. The alpha diversity of gut microbiota in the ASD group showed a significant difference when compared with those of the HC groups (Figure 1). Specifically, these diversity indices showed a similar trend, that is, the alpha diversity of gut microbiota was higher in the HC group than in the ASD group (Figures 1A−C, online Supplementary Figure S2). The Anosim-analysis also showed that the differences between HC and ASD groups were significantly greater than those within the group, indicating meaningful grouping (Figure 1D). Figure 1E illustrated that the samples obtained from the ASD group had 1,152 features, while only 490 features were detected in the samples from the HC group, both groups shared a total of 432 features. These findings suggested that a significant number of features were present exclusively in the ASD group, indicating a higher relative abundance of certain species in the gut microbiota of children with ASD compared to healthy children. The pCoA principal component analysis based on bray-curtis clearly showed the microbiota differences between the ASD and HC groups. There were significant differences along the PCO1 axis, accounting for 12.56% of the total variation (Figure 1F). Sample cluster evolutionary tree analysis also effectively clustered the ASD and HC group (Figure 1G). In conclusion, the species richness and diversity of gut microbiota in children with ASD were significantly different from those in HC group.
Figure 1. Diversity of gut microbiota are different across ASD and HC. (A−C) Index of Alpha diversity (the Shannon, simpson, and Chao1). (D) Anosim-analysis results, between represents the difference between groups; others are within groups; the greater the distance is, the greater the difference is; and the thickness is the sample size. (E) UpSet plot of intersections and unique among different group. (F) pCoA of beta diversity. (G) Cluster evolutionary tree analysis. *p < 0.05, **p < 0.01 and ***p < 0.001.
Differences of microbial comparison between ASD and healthy controls
Next, the composition of gut microbiota in ASD, and HC group was further analyzed. At the phylum level, Firmicutes, Bacteroidetes, Proteobacteria and Actinobacteria were found to be dominant in the different groupings (Figure 2A). In comparison to the HC group, the ASD group exhibited higher levels of Firmicutes and a significantly lower Bacteroidetes/Firmicutes ratio (ASD = 0.54, HC = 0.89). At the genus level, g_Bacteroides, g_Faecalibacterium, g_Prevotella_9, g_Bifidobacterium, g_Phascolarctobacterium and g_Blautia were found to be dominant in the different groupings (Figure 2B). At the phylum level, the ASD group showed a significant decrease in p_Lentisphaerae, p_Bacteroidetes, p_Euryarchaeota, p_Patescibacteria, while p_Fusobacteria, p_Firmicutes, and p_Verrucomicrobia were significantly increased compared with the HC group (Figure 3A). At the genus level, we observed 316 significantly different species between ASD and HC groups (Figure 3B). Among them, the top 20 genera that showed a significant increase in the ASD group compared with the HC group included g_Blautia, g_Lachnoclostridium, g_Lachnospira, g_Lachnospiraceae_NK4A136_group, g_Eubacterium_eligens_group, g_Faecalibacterium, g_Dialister, g_Parasutterella, the significantly reduced bacterial flora included g_Prevotella_9, g_Alistipes, g_Agathobacter (Figure 3C). Then analysis at the species level was carried out to find out the bacteria with significant differences, among which the s_Bacteroides_fragilis in the ASD group was significantly increased compared with the HC group, while s_Bacteroides_coprocola_DSM_17136, s_Bacteroides_plebeius, s_Bacteroides_stercoris_ATCC_43183, and the s_Parabacteroides_merdae was significantly reduced (Figure 3D). It showed that indicating that these differential flora may be involved in the gastrointestinal symptomsof ASD.
Figure 2. Phylum and Genera of gut microbiota are strikingly different across ASD and HC. (A) Relative abundance of the indicated phylum. (B) Relative abundance of the indicated genus. (C) Genera co-occurrence network between ASD and HC group based on the Spearman correlation algorithms. Each node presents a bacterial genus. The node size indicates the weighted value of each genus per group, and the thickness of the line represents the Spearman coefficient (p < 0.05, r ≥ 0.7).
Figure 3. The shift of gut microbiota in ASD and HC. (A)Heatmap of selected differentially abundant phylum in ASD and HC (p < 0.05). (B) Heatmap of selected differentially abundant genus in ASD and HC (p < 0.05). (C,D) The top 20 of differentially abundant analysis between ASD and HC at the level of genus (C) and species (D). (E) LDA scores for the bacterial taxa showed different abundant between ASD and HC (LAD score > 3.5) (p < 0.05). *p < 0.05, **p < 0.01, and ***p < 0.001.
According to the former study, it demonstrated that the complex microbial ecosystem of the human intestinal tract is unevenly influenced by individualtaxa within different microbial communities (Swiss IBD Cohort Investigators et al., 2019). To characterize potential relationships among bacteria in gut microbial communities, we further constructed co-occurrence networks for each taxon based on significant Spearman correlations. The HC taxa were mainly composed of six co-occurrence networks consisting of dispersed genera from four primary phyla: Proteobacteria, Firmicutes, Actinobacteria and Bacteroidetes, while the ASD taxa were mainly composed of 15 co-occurrence networks consisting of dispersed genera from seven primary phyla: Proteobacteria, Firmicutes, Actinobacteria, Bacteroidetes, Chloroflexi, Acidobacteria, and Patescibacteria (Figure 2C). Two group displayed a co-occurrence network with strong positive correlation among genera. As shown in Figure 2C, the microbial community of the ASD group featured more complicated network. The correlation between the microbiota in the ASD group was distinctly increased compared to that of HC group. To quantify such differences, the average degree of edges (connections) and nodes (genera) were counted in the three microbial networks [ASD (49.78), HC (26.71)]. Taken together, the above analyses further indicating that there may be a disorder in the intestinal microecology of children with ASD.
The linear discriminant analysis (LDA) distribution diagram analysis (LAD score > 3.5) showed a clear alteration of the microbiota characterized by higher p_Frimicutes and p_Actinobacteria levels in ASD individuals (Figure 3E). However, p_Bacteroidetes levels were significantly decreased in ASD group (Figure 3E). The g_Faecalibacterium, g_Blautia, g_Eubacterium_eligens_group, g_Parasutterella, g_Lachnospiraceae_NK4A136_group and g_Veillonella were more abundant in ASD group, while the g_Prevotella 9, g_Agathobacter were more abundant in HC group (Figure 3E). The f_Ruminoccaceae, f_Lachnospiraceae was more abundant in ASD group, while f_Prevotellaceae and f_Enterobacteriaeae was more abundant in HC group (online Supplementary Figure S3). These apparently distinct species could be used as non-invasive biomarkers for the early diagnosis of autism. Taken together, compared with the HC group, the microbial composition of the ASD group showed alterations, further indicating that gut microecology dysbiosis existed in children with ASD.
Significant differences in metabolic pathways between the ASD group and the HC group
Ninety-seven Kegg pathways were different between the ASD and HC groups By KEGG pathway analysis (Figure 4A). In particular, among the Top 20 KEGG pathways, One carbon pool By folate pathways, carbon fixation In photosynthetic organisms pathways were significantly less enriched In The ASD group compared with The HC group, while The biosynthesis of ansamycins pathways, valine, leucine and isoleucine biosynthesis pathways, C5 − branched dibasic acid metabolism pathways, pentose phosphate pathway pathways, fatty acid biosynthesis pathways, lysine biosynthesis pathways were significantly more enriched (Figure 4B).
Figure 4. Different metabolic patterns in ASD and HC. (A) The average abundance of KEGG pathway differentially enriched in ASD and HC (p < 0.05). (B) The top 20 of KEGG pathway differentially enriched in ASD and HC (p < 0.05). (C) The heatmap of correlation between differential gut microbiota and differential metabolites (p < 0.05). Red color represents positive correlation, blue color represents negative correlation, and darker color corresponds to a stronger correlation. *p < 0.05, **p < 0.01, and ***p < 0.001.
To explore the potential relationships between the gut microbiome changes and metabolic products, a correlation matrix was generated using Spearman correlation (Figure 4C). The abundance of most species such as g_Blautia, g_Prevotella_9, g_Agathobacter, and g_Faecalibacterium, was positively correlated with the level of Penicillin_and_cephalosporin_biosynthesis, Carbon_fixation_in_photosynthetic_organisms, Arachidonic_acid_metabolism, Staphylococcus_aureus_infection. In summary, the abnormal gut microbiota in children with ASD leads to abnormal gut metabolites, which may be the potential reason for the effect of gut microbiota on ASD, although further research was needed to validate the results.
Discussion
Based on our study, we found significant differences in the diversity, structure, and relative abundance of intestinal microflora between ASD children and healthy children. Therefore, it is crucial to detect the gut microbiota in ASD children for intervention treatment.
Our results showed a significant decrease in gut microbiota diversity in ASD children compared with the HC group (Dan et al., 2020; Wan et al., 2022). These results were consistent with some previous studies (Carissimi et al., 2019; Wang et al., 2019; Dan et al., 2020; Ye et al., 2021). In addition, the Bacteroidetes/Firmicutes ratio of gut microbiota in ASD children was significantly reduced, which was consistent with the observation of gut microbiota in other autistic patients (Tomova et al., 2015; Strati et al., 2017). We found that different bacterial genera, such as g_Blautia, g_Faecalibacterium, g_Lachnospiraceae_NK4A136_group, g_Parasutterella, g_Eubacterium_eligens_group, g_Parasutterella, g_Akkermansia, and g_Anaerostipes significantly increased in the ASD group, while g_Prevotella_9, g_Alistipes, and g_Agathobacter were significantly reduced. These differential gut microbiotas can be used as biomarkers for early screening of ASD. In fact, the research results about ASD gut microbiota are often contradictory (Ho et al., 2020; Iglesias-Vázquez et al., 2020). At the genus level, more controversial findings have been reported, with g_Akkermansia reported to be more abundant in children with ASD (de Angelis et al., 2013; Li et al., 2021), but some studies have shown that the relative abundance of g_Akkermansia is significantly reduced in the gut microbiota of children with ASD (Wang et al., 2011; Zhang et al., 2018). Our data found that the relative abundance of g_Akkermansia in autistic patients was significantly increased in the ASD group. In addition, studies have shown that the relative abundance of g_faecalibacterium in the ASD sleep disorder group was significantly lower compared to healthy controls (Hua et al., 2020), while a meta-analysis reported that g_Faecalibacterium was more abundant in ASD children (Iglesias-Vázquez et al., 2020), consistent with our current research results. Some studies have shown that the Ruminococcaceae family, including Faecalibacterium, Subdoligranulum varia, Clostridium, and Eubacterium, is positively correlated with caproic acid levels (Dan et al., 2020), and ASD is associated with higher levels of caproic acid in the blood (Wang et al., 2012). The acetic acid produced by these intestinal bacteria is a ligand for GPR84, which can enhance the production of pro-inflammatory cytokine IL-12 and its subunit IL-12 P40 induced by lipopolysaccharide (Zhao et al., 2012). It has been reported that IL-12 or IL-12 P40 significantly increased in ASD children (Ashwood et al., 2011; Saghazadeh et al., 2019).
In addition, it was noted that g_Prevotella_9 was depleted in the gut microbiota of ASD children, while g_Blautia was enriched. The results of our g_Blautiade study are contrary to those of previous studies (Xie et al., 2022). However, the decline in Prevotella levels in ASD children was consistent with the results of previous studies (Dan et al., 2020; Zhang et al., 2021). Prevotella was the main bacterial species in the human gut, which could use polysaccharide to produce succinic acid (Kovatcheva-Datchary et al., 2015; Strati et al., 2017). Succinic acid could bind to the succinic acid receptor GPR91 on the surface of dendritic cells to enhance the immune response of antigen-specific T cells and protect the health of the host (Rubic et al., 2008). There is growing evidence that many ASD children have immune dysfunction (Ashwood et al., 2011). Therefore, it was speculated that the immune dysfunction of ASD might be related to the decrease of Prevatella bacteria in the gut, which further affected the physiological and behavioral characteristics of children with ASD. Why do different studies reach different conclusions, which we suspect that it may be due to differences in genetic background, varying manifestations of autism symptoms, as well as differences in age and living environment. Anyhow, these results again indicated that the gut microbiota of children with ASD was significantly different from that of non-ASD individuals.
According to the former study, it showed that Ecological networks of gut microbiota are considered critical for host health and well-being, because it shows that beneficial symbionts and their associated functions are maintained over time (Lozupone et al., 2012; Relman, 2012). In our study, the gut microbial community of children with ASD showed a more complex ecological network where more interactions between bacterial species were observed. It has been shown that highly cooperative microbial communities are perceived as less stable (Coyte and Rakoff-Nahoum, 2019). As previous studies have shown, only gradually increasing the proportion of cooperative interactions within a community almost always reduces the overall rate of return and the likelihood of stability (Coyte et al., 2015). As for key taxa in the network in ASD, Bacteroidetes produce propionic acid and other short-chain fatty acids (SCFA). Rat models administered with propionic acid show increased restrictive/repetitive behavior as well as impaired social behavior (MacFabe et al., 2011). In our study, Bacteroidetes showed the most interactions and identified as key species in the network in children with ASD. The associations between these bacteria and ASD requires further exploration. Besides this, Clostridium genus associated with neurological disorders (Averina et al., 2020) were connected closely with each other under the same genera. Interactions between harmful bacteria may promote and accelerate disease progression (Singer, 2010). Taken together, we observed an unstable and unfavorable ecosystem in the gut microbiome of ASD.
Numerous studies have shown that gut microbiota plays a key role in the biological and physiological characteristics of neurodevelopment, and the establishment of early community relationships among microorganisms may affect the neurodevelopment of children (Dinan and Cryan, 2017; Srikantha and Mohajeri, 2019; Liu et al., 2022). The results showed that, compared with HC group, abnormal microbial function in ASD patients was mainly concentrated in metabolic pathways related to fatty acid, amino acid and aromatic compound metabolism, and partly involved in the metabolism of neurotransmitters. Previous studies have shown that these functions were related to individual nervous system development, neurotransmitter biosynthesis and neuronal response regulation (Needham et al., 2020). For example, abnormalities in the bacterial pathway of tryptophan metabolism have been observed in patients with ASD, and its production of neuroprotective canine and neurotoxic quinolinic acid has been shown to be significantly associated with the severity of ASD (de Angelis et al., 2013; Lou et al., 2022). Our functional analysis of the gut microbiome provided some clues to these potential mechanisms that need to be further investigated in the future.
The differential microbiota and abnormal metabolic pathways identified in our study may have important clinical implications in the early diagnosis and treatment of autism. A large number of previous studies have shown that gut microbiota can be used as a new method for early screening and treatment of autism. Xinyan Xie et al. conducted a case–control study on 101 ASD children and 103 healthy controls in China, and found that Actinobacteria, Proteobacteria, Enterobacteriaceae, and Escherichia were significantly higher in ASD children. However, the numbers of monospheraceae, Blautia and unclassified lachnospiraceae were significantly reduced in ASD children. The gut microbiota of ASD children may have disorders in functional pathways such as amino acid metabolism, coenzyme and vitamin metabolism, and AMPK signaling pathway (Xie et al., 2022). Zhou Dan et al. performed 16S sequencing and metagenomic sequencing analysis on 143 children with clinically diagnosed ASD and 143 age- and sex-matched typically developing (TD) individuals. They found that ASD was associated with alterations in gut microbial profiles and abnormal metabolic activity. ASD patients showed gut ecological imbalance at the phylum, genus and species levels. The findings provide the possibility for future ASD interventions targeting specific microbiota related to neurotransmitter metabolism (Dan et al., 2020). Yingxin Zhao et al. analyzed and compared the gut microbiota of 36 ASD children with gastrointestinal symptoms and 40 TD children. Compared with TD, the gut microbiota of ASD patients with GI symptoms showed a decrease in alpha diversity and butyrate-producing bacteria (such as Faecalibacterium and Coprococcus). The results support the use of gut microbiota as a potential biomarker for early identification of ASD and intervention targeting specific gut microbiota (Zhao et al., 2023).
A strength of this study was the relatively large sample size compared with previous studies on the gut microbiota of ASD children (Supplementary Table S1). However, there were several limitations that need to be addressed in future studies. First, this study was observational only, and longitudinal studies were needed to explore how a disturbed gut microbiome contributes to the development of ASD symptoms. Secondly, due to the availability of data, we were only able to obtain and analyze gut microbiota data, and lacked other data (such as diet, clinical indications, physiology and biochemistry, metabolomics, etc.), which limited the ability of in-depth causal inference. These problems were expected to be solved in the following research.
In summary, our study showed that children with ASD exhibited intestinal microecological disorders at the phylum, genus, and species levels compared with healthy children. These findings supported previous studies on gut microbiota changes in children with ASD and contributed to the development of new approaches for early screening and treatment of ASD based on gut microbiota. However, due to the heterogeneity of ASD phenotypes, further studies in independent populations with larger sample sizes are still needed.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Ethics statement
The studies involving humans were approved by his study was approved by the Ethics Committee of Second Affiliated Hospital, Army Medical University (Ethics No. 2023–001-01). Written informed consents were signed by all the children’s legal guardians prior to the study. The clinical trial number: ChiCTR2300074832. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin.
Author contributions
HL: Conceptualization, Formal analysis, Investigation, Methodology, Project administration, Software, Validation, Writing – original draft, Writing – review & editing. WG: Investigation, Resources, Writing – review & editing. SL: Investigation, Resources, Writing – review & editing. BS: Writing – review & editing. NL: Writing – review & editing. DX: Investigation, Writing – review & editing. ZD: Writing – review & editing. DL: Investigation, Writing – review & editing. WC: Writing – review & editing. WF: Writing – review & editing. JZhe: Writing – review & editing. JZhu: Conceptualization, Funding acquisition, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The National Natural Science Foundation of China (no. 81900690) and the Natural Science Foundation of Chongqing (no. cstc2020jcyj-msxmX0065).
Acknowledgments
The authors thank Li Tang for technical assistance.
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.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1326870/full#supplementary-material
SUPPLEMENTARY FIGURE S1 | Study fowchart.
Supplementary FIGURE S2 | Alpha diversity across microbiome samples of ASD and HC.
Supplementary FIGURE S3 | Cladograms generated by LEfSe indicating differences in the bacterial taxa between ASD and HC.
Supplementary Table S1 | Characteristics of study participants.
Footnotes
References
Adams, J. B., Johansen, L. J., Powell, L. D., Quig, D., and Rubin, R. A. (2011). Gastrointestinal flora and gastrointestinal status in children with autism – comparisons to typical children and correlation with autism severity. BMC Gastroenterol. 11:22. doi: 10.1186/1471-230X-11-22
Ashwood, P., Krakowiak, P., Hertz-Picciotto, I., Hansen, R., Pessah, I., and van de Water, J. (2011). Elevated plasma cytokines in autism spectrum disorders provide evidence of immune dysfunction and are associated with impaired behavioral outcome. Brain Behav. Immun. 25, 40–45. doi: 10.1016/j.bbi.2010.08.003
Averina, O. V., Kovtun, A. S., Polyakova, S. I., Savilova, A. M., Rebrikov, D. V., and Danilenko, V. N. (2020). The bacterial neurometabolic signature of the gut microbiota of young children with autism spectrum disorders. J. Med. Microbiol. 69, 558–571. doi: 10.1099/jmm.0.001178
Carissimi, C., Laudadio, I., Palone, F., Fulci, V., Cesi, V., Cardona, F., et al. (2019). Functional analysis of gut microbiota and immunoinflammation in children with autism spectrum disorders. Dig. Liver Dis. 51, 1366–1374. doi: 10.1016/j.dld.2019.06.006
Coyte, K. Z., and Rakoff-Nahoum, S. (2019). Understanding competition and cooperation within the mammalian gut microbiome. Curr. Biol. 29, R538–R544. doi: 10.1016/j.cub.2019.04.017
Coyte, K. Z., Schluter, J., and Foster, K. R. (2015). The ecology of the microbiome: networks, competition, and stability. Science 350, 663–666. doi: 10.1126/science.aad2602
Dan, Z., Mao, X., Liu, Q., Guo, M., Zhuang, Y., Liu, Z., et al. (2020). Altered gut microbial profile is associated with abnormal metabolism activity of autism Spectrum disorder. Gut Microbes 11, 1246–1267. doi: 10.1080/19490976.2020.1747329
de Angelis, M., Piccolo, M., Vannini, L., Siragusa, S., de Giacomo, A., Serrazzanetti, D. I., et al. (2013). Fecal microbiota and metabolome of children with autism and pervasive developmental disorder not otherwise specified. PLoS One 8:e76993. doi: 10.1371/journal.pone.0076993
Dinan, T. G., and Cryan, J. F. (2017). Gut instincts: microbiota as a key regulator of brain development, ageing and neurodegeneration. J. Physiol. 595, 489–503. doi: 10.1113/JP273106
Ding, X., Xu, Y., Zhang, X., Zhang, L., Duan, G., Song, C., et al. (2020). Gut microbiota changes in patients with autism spectrum disorders. J. Psychiatr. Res. 129, 149–159. doi: 10.1016/j.jpsychires.2020.06.032
Ferguson, B. J., Dovgan, K., Takahashi, N., and Beversdorf, D. Q. (2019). The relationship among gastrointestinal symptoms, problem Behaviors, and internalizing symptoms in children and adolescents with autism Spectrum disorder. Front. Psych. 10:194. doi: 10.3389/fpsyt.2019.00194
Fowlie, G., Cohen, N., and Ming, X. (2018). The perturbance of microbiome and gut-brain Axis in autism Spectrum disorders. Int. J. Mol. Sci. 19:2251. doi: 10.3390/ijms19082251
Frye, R. E., Slattery, J. C., and Quadros, E. V. (2017). Folate metabolism abnormalities in autism: potential biomarkers. Biomark. Med 11, 687–699. doi: 10.2217/bmm-2017-0109
Ho, L., Tong, V., Syn, N., Nagarajan, N., Tham, E. H., Tay, S. K., et al. (2020). Gut microbiota changes in children with autism spectrum disorder: a systematic review. Gut Pathog. 12:6. doi: 10.1186/s13099-020-0346-1
Hua, X., Zhu, J., Yang, T., Guo, M., Li, Q., Chen, J., et al. (2020). The gut microbiota and associated metabolites are altered in sleep disorder of children with autism Spectrum disorders. Front. Psych. 11:855. doi: 10.3389/fpsyt.2020.00855
Iannone, L. F., Preda, A., Blottière, H. M., Clarke, G., Albani, D., Belcastro, V., et al. (2019). Microbiota-gut brain axis involvement in neuropsychiatric disorders. Expert. Rev. Neurother. 19, 1037–1050. doi: 10.1080/14737175.2019.1638763
Iglesias-Vázquez, L., Van Ginkel, R. G., Arija, V., and Canals, J. (2020). Composition of gut microbiota in children with autism Spectrum disorder: a systematic review and Meta-analysis. Nutrients 12:792. doi: 10.3390/nu12030792
Joon, P., Kumar, A., and Parle, M. (2021). What is autism? Pharmacol. Rep. 73, 1255–1264. doi: 10.1007/s43440-021-00244-0
Kandeel, W. A., Meguid, N. A., Bjørklund, G., Eid, E. M., Farid, M., Mohamed, S. K., et al. (2020). Impact of Clostridium Bacteria in children with autism Spectrum disorder and their anthropometric measurements. J. Mol. Neurosci. 70, 897–907. doi: 10.1007/s12031-020-01482-2
Kovatcheva-Datchary, P., Nilsson, A., Akrami, R., Lee, Y. S., de Vadder, F., Arora, T., et al. (2015). Dietary Fiber-induced improvement in glucose metabolism is associated with increased abundance of Prevotella. Cell Metab. 22, 971–982. doi: 10.1016/j.cmet.2015.10.001
Lacorte, E., Gervasi, G., Bacigalupo, I., Vanacore, N., Raucci, U., and Parisi, P. (2019). A systematic review of the microbiome in children with neurodevelopmental disorders. Front. Neurol. 10:727. doi: 10.3389/fneur.2019.00727
Li, N., Chen, H., Cheng, Y., Xu, F., Ruan, G., Ying, S., et al. (2021). Fecal microbiota transplantation relieves gastrointestinal and autism symptoms by improving the gut microbiota in an open-label study. Front. Cell. Infect. Microbiol. 11:759435. doi: 10.3389/fcimb.2021.759435
Li, N., Yang, J., Zhang, J., Liang, C., Wang, Y., Chen, B., et al. (2019). Correlation of gut microbiome between ASD children and mothers and potential biomarkers for risk assessment. Genomics Proteom. Bioinform. 17, 26–38. doi: 10.1016/j.gpb.2019.01.002
Liu, J., Gao, Z., Liu, C., Liu, T., Gao, J., Cai, Y., et al. (2022). Alteration of gut microbiota: new strategy for treating autism Spectrum disorder. Front. Cell Dev. Biol. 10:792490. doi: 10.3389/fcell.2022.792490
Liu, S., Li, E., Sun, Z., Fu, D., Duan, G., Jiang, M., et al. (2019). Altered gut microbiota and short chain fatty acids in Chinese children with autism spectrum disorder. Sci. Rep. 9:287. doi: 10.1038/s41598-018-36430-z
Lord, C., Elsabbagh, M., Baird, G., and Veenstra-Vanderweele, J. (2018). Autism spectrum disorder. Lancet 392, 508–520. doi: 10.1016/S0140-6736(18)31129-2
Lou, M., Cao, A., Jin, C., Mi, K., Xiong, X., Zeng, Z., et al. (2022). Deviated and early unsustainable stunted development of gut microbiota in children with autism spectrum disorder. Gut 71, 1588–1599. doi: 10.1136/gutjnl-2021-325115
Lozupone, C. A., Stombaugh, J. I., Gordon, J. I., Jansson, J. K., and Knight, R. (2012). Diversity, stability and resilience of the human gut microbiota. Nature 489, 220–230. doi: 10.1038/nature11550
Luna, R. A., Oezguen, N., Balderas, M., Venkatachalam, A., Runge, J. K., Versalovic, J., et al. (2017). Distinct microbiome-neuroimmune signatures correlate with functional abdominal pain in children with autism Spectrum disorder. Cell. Mol. Gastroenterol. Hepatol. 3, 218–230. doi: 10.1016/j.jcmgh.2016.11.008
MacFabe, D. F., Cain, N. E., Boon, F., Ossenkopp, K. P., and Cain, D. P. (2011). Effects of the enteric bacterial metabolic product propionic acid on object-directed behavior, social behavior, cognition, and neuroinflammation in adolescent rats: relevance to autism spectrum disorder. Behav. Brain Res. 217, 47–54. doi: 10.1016/j.bbr.2010.10.005
Maenner, M. J., Shaw, K. A., Baio, J., Washington, A., Patrick, M., DiRienzo, M., et al. (2020). Prevalence of autism Spectrum disorder among children aged 8 years - autism and developmental disabilities monitoring network, 11 sites, United States, 2016. MMWR Surveill. Summ. 69, 1–12. doi: 10.15585/mmwr.ss6904a1
Modabbernia, A., Velthorst, E., and Reichenberg, A. (2017). Environmental risk factors for autism: an evidence-based review of systematic reviews and meta-analyses. Mol. Autism. 8:13. doi: 10.1186/s13229-017-0121-4
Needham, B. D., Kaddurah-Daouk, R., and Mazmanian, S. K. (2020). Gut microbial molecules in behavioural and neurodegenerative conditions. Nat. Rev. Neurosci. 21, 717–731. doi: 10.1038/s41583-020-00381-0
Nogay, N. H., and Nahikian-Nelms, M. (2021). Can we reduce autism-related gastrointestinal and behavior problems by gut microbiota based dietary modulation? A review. Nutr. Neurosci. 24, 327–338. doi: 10.1080/1028415X.2019.1630894
Pulikkan, J., Maji, A., Dhakan, D. B., Saxena, R., Mohan, B., Anto, M. M., et al. (2018). Gut microbial dysbiosis in Indian children with autism Spectrum disorders. Microb. Ecol. 76, 1102–1114. doi: 10.1007/s00248-018-1176-2
Relman, D. A. (2012). The human microbiome: ecosystem resilience and health. Nutr. Rev. 70, S2–S9. doi: 10.1111/j.1753-4887.2012.00489.x
Rose, D. R., Yang, H., Serena, G., Sturgeon, C., Ma, B., Careaga, M., et al. (2018). Differential immune responses and microbiota profiles in children with autism spectrum disorders and co-morbid gastrointestinal symptoms. Brain Behav. Immun. 70, 354–368. doi: 10.1016/j.bbi.2018.03.025
Rubic, T., Lametschwandtner, G., Jost, S., Hinteregger, S., Kund, J., Carballido-Perrig, N., et al. (2008). Triggering the succinate receptor GPR91 on dendritic cells enhances immunity. Nat. Immunol. 9, 1261–1269. doi: 10.1038/ni.1657
Saghazadeh, A., Ataeinia, B., Keynejad, K., Abdolalizadeh, A., Hirbod-Mobarakeh, A., and Rezaei, N. (2019). A meta-analysis of pro-inflammatory cytokines in autism spectrum disorders: effects of age, gender, and latitude. J. Psychiatr. Res. 115, 90–102. doi: 10.1016/j.jpsychires.2019.05.019
Sanlier, N., and Kocabas, Ş. (2021). The effect of probiotic, prebiotic and gut microbiota on ASD: a review and future perspectives. Crit. Rev. Food Sci. Nutr. 63, 2319–2330. doi: 10.1080/10408398.2021.1973957
Sgritta, M., Dooling, S. W., Buffington, S. A., Momin, E. N., Francis, M. B., Britton, R. A., et al. (2019). Mechanisms underlying microbial-mediated changes in social behavior in mouse models of autism Spectrum disorder. Neuron 101, 246–259.e6. doi: 10.1016/j.neuron.2018.11.018
Sharon, G., Cruz, N. J., Kang, D. W., Gandal, M. J., Wang, B., Kim, Y. M., et al. (2019). Human gut microbiota from autism Spectrum disorder promote Behavioral symptoms in mice. Cells 177, 1600–1618. doi: 10.1016/j.cell.2019.05.004
Singer, M. (2010). Pathogen-pathogen interaction: a syndemic model of complex biosocial processes in disease. Virulence 1, 10–18. doi: 10.4161/viru.1.1.9933
Srikantha, P., and Mohajeri, M. H. (2019). The possible role of the microbiota-gut-brain-Axis in autism Spectrum disorder. Int. J. Mol. Sci. 20:2115. doi: 10.3390/ijms20092115
Strati, F., Cavalieri, D., Albanese, D., de Felice, C., Donati, C., Hayek, J., et al. (2017). New evidences on the altered gut microbiota in autism spectrum disorders. Microbiome 5:24. doi: 10.1186/s40168-017-0242-1
Swiss IBD Cohort InvestigatorsYilmaz, B., Juillerat, P., Øyås, O., Ramon, C., Bravo, F. D., et al. (2019). Microbial network disturbances in relapsing refractory Crohn's disease. Nat. Med. 25, 323–336. doi: 10.1038/s41591-018-0308-z
Tomova, A., Husarova, V., Lakatosova, S., Bakos, J., Vlkova, B., Babinska, K., et al. (2015). Gastrointestinal microbiota in children with autism in Slovakia. Physiol. Behav. 138, 179–187. doi: 10.1016/j.physbeh.2014.10.033
Tomova, A., Soltys, K., Repiska, G., Palkova, L., Filcikova, D., Minarik, G., et al. (2020). Specificity of gut microbiota in children with autism spectrum disorder in Slovakia and its correlation with astrocytes activity marker and specific behavioural patterns. Physiol. Behav. 214:112745. doi: 10.1016/j.physbeh.2019.112745
Valentino, F., Bruno, L. P., Doddato, G., Giliberti, A., Tita, R., Resciniti, S., et al. (2021). Exome sequencing in 200 intellectual disability/autistic patients: new candidates and atypical presentations. Brain Sci. 11:936. doi: 10.3390/brainsci11070936
Wan, Y., Zuo, T., Xu, Z., Zhang, F., Zhan, H., Chan, D., et al. (2022). Underdevelopment of the gut microbiota and bacteria species as non-invasive markers of prediction in children with autism spectrum disorder. Gut 71, 910–918. doi: 10.1136/gutjnl-2020-324015
Wang, L., Christophersen, C. T., Sorich, M. J., Gerber, J. P., Angley, M. T., and Conlon, M. A. (2011). Low relative abundances of the mucolytic bacterium Akkermansia muciniphila and Bifidobacterium spp. in feces of children with autism. Appl. Environ. Microbiol. 77, 6718–6721. doi: 10.1128/AEM.05212-11
Wang, L., Christophersen, C. T., Sorich, M. J., Gerber, J. P., Angley, M. T., and Conlon, M. A. (2012). Elevated fecal short chain fatty acid and ammonia concentrations in children with autism spectrum disorder. Dig. Dis. Sci. 57, 2096–2102. doi: 10.1007/s10620-012-2167-7
Wang, T., Guo, H., Xiong, B., Stessman, H. A. F., Wu, H., Coe, B. P., et al. (2016). De novo genic mutations among a Chinese autism spectrum disorder cohort. Nat. Commun. 7:13316. doi: 10.1038/ncomms13316
Wang, F., Lu, L., Wang, S. B., Zhang, L., Ng, C. H., Ungvari, G. S., et al. (2018). The prevalence of autism spectrum disorders in China: a comprehensive meta-analysis. Int. J. Biol. Sci. 14, 717–725. doi: 10.7150/ijbs.24063
Wang, M., Wan, J., Rong, H., He, F., Wang, H., Zhou, J., et al. (2019). Alterations in gut glutamate metabolism associated with changes in gut microbiota composition in children with autism Spectrum disorder. Msystems 4:e00321-18. doi: 10.1128/mSystems.00321-18
Wang, J., Zheng, J., Shi, W., du, N., Xu, X., Zhang, Y., et al. (2018). Dysbiosis of maternal and neonatal microbiota associated with gestational diabetes mellitus. Gut 67, 1614–1625. doi: 10.1136/gutjnl-2018-315988
Xie, X., Li, L., Wu, X., Hou, F., Chen, Y., Shi, L., et al. (2022). Alteration of the fecal microbiota in Chinese children with autism spectrum disorder. Autism Res. 15, 996–1007. doi: 10.1002/aur.2718
Ye, F., Gao, X., Wang, Z., Cao, S., Liang, G., He, D., et al. (2021). Comparison of gut microbiota in autism spectrum disorders and neurotypical boys in China: a case-control study. Synth. Syst. Biotechnol. 6, 120–126. doi: 10.1016/j.synbio.2021.03.003
Zhang, M., Ma, W., Zhang, J., He, Y., and Wang, J. (2018). Analysis of gut microbiota profiles and microbe-disease associations in children with autism spectrum disorders in China. Sci. Rep. 8:13981. doi: 10.1038/s41598-018-32219-2
Zhang, Q., Zou, R., Guo, M., Duan, M., Li, Q., and Zheng, H. (2021). Comparison of gut microbiota between adults with autism spectrum disorder and obese adults. PeerJ 9:e10946. doi: 10.7717/peerj.10946
Zhao, X., Ross, E. J., Wang, Y., and Horwitz, B. H. (2012). Nfkb1 inhibits LPS-induced IFN-β and IL-12 p40 production in macrophages by distinct mechanisms. PLoS One 7:e32811. doi: 10.1371/journal.pone.0032811
Zhao, Y., Wang, Y., Meng, F., Chen, X., Chang, T., Huang, H., et al. (2023). Altered gut microbiota as potential biomarkers for autism Spectrum disorder in early childhood. Neuroscience 523, 118–131. doi: 10.1016/j.neuroscience.2023.04.029
Keywords: gut microbiota, autism spectrum disorder, Chinese children, metabolic pathways, gut-brain axis
Citation: Li H, Guo W, Li S, Sun B, Li N, Xie D, Dong Z, Luo D, Chen W, Fu W, Zheng J and Zhu J (2024) Alteration of the gut microbiota profile in children with autism spectrum disorder in China. Front. Microbiol. 14:1326870. doi: 10.3389/fmicb.2023.1326870
Edited by:
Yu Wang, China Academy of Chinese Medical Sciences, ChinaReviewed by:
Yi Xu, Hefei University of Technology, ChinaYingli Jing, Capital Medical University, China
Copyright © 2024 Li, Guo, Li, Sun, Li, Xie, Dong, Luo, Chen, Fu, Zheng and Zhu. 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.
*Correspondence: Ji Zheng, Jizheng@tmmu.edu.cn; Jingzhen Zhu, zhujz1114@163.com
†These authors have contributed equally to this work