A 12‐genus Bacterial Signature Identifies a Group of Severe Autistic Children with Differential Sensory Behavior and Brain Structures
Kai Shi,Lingli Zhang,Juehua Yu,Zilin Chen,Senying Lai,Xingzhong Zhao,Wei‐Guang Li,Qiang Luo,Wei Lin,Jianfeng Feng,Peer Bork,Xing‐Ming Zhao,Fēi Li
DOI: https://doi.org/10.1002/ctm2.314
IF: 8.554
2021-01-01
Clinical and Translational Medicine
Abstract:Autism spectrum disorder (ASD) refers to a group of heterogeneous neurodevelopmental disorders with neuropsychological and behavioral deficits,1 whose accurate subtyping is important but no efficient method is available.2, 3 Emerging studies suggest a possible mechanism that involves the microbiota–gut–brain axis4, 5; however, the associations underlying the gut microbiota community, brain structure,6 and behavioral symptoms are poorly defined in ASD. Herein, we aimed to investigate and identify a potential microbial signature that is linked to brain structure variations and severity of symptoms in patients with ASD by using a machine learning framework.7 All methods are fully detailed in Materials and Methods in the Supporting Information. Among 128 eligible patients recruited, the initially collected cohort was used as the discovery set (n = 78), and the subsequently collected cohort was used as the test set (n = 50). The clinical characteristics including age, gender, and CARS8/ADOS9 total scores were well balanced (Table 1). Two other public sets were used as the validation sets (SRP093968, PRJEB15418). First, we developed a computational pipeline utilizing a greedy search strategy and identified a microbial signature of 12 bacterial genera in the discovery set, by which the subjects were classified into two subpopulations, denoted as ASDmp1 and ASDmp2 (Figure 1A and B). We subsequently trained a k-nearest neighbors (kNN) classifier using the abundance profile of the 12 bacterial genera, which performed very well with the classification accuracy of 0.87 in leave-one-out cross-validation (Figure S1). Distinct microbial patterns in subpopulations of children with ASD. (A) A signature of 12 genera that was most associated with the behavioral characteristics using Spearman correlation. (B) Two ASD subpopulations with a 12-genus signature in the discovery set. Samples were categorized by hierarchical clustering using average linkage with the similarity of 12 genera on the Bray-Curtis metric. The first color bar indicates the severity of the disorder (CARS total score: red ≥ 36, blue < 36), and the second color bar shows the CARS scores. (C) Stacked horizontal bar charts depict the variability in genus-level composition for individuals by groups for 12 microbial markers. (D and E) Principal component analysis (PCoA) of a 12-genus signature based on Bray-Curtis dissimilarity for the discovery and test sets. (F and G) The comparison of the relative abundance ratio between Parabacteriodes and Blautia in the discovery and test sets. (H) α-Diversity in ASDmp1 and ASDmp2 as measured by observed amplicon sequence variants. (I) Differential MetaCyc pathway analysis as calculated by edgeR. Fold change (ASDmp1/ASDmp2) as a factor for Benjamini–Hochberg corrected p-values. CARS, Children Autism Rating Scale; ASDmp1, ASD subpopulation 1; ASDmp2, ASD subpopulation 2 The diversity analyses of the gut microbiome indicated that the two subpopulations ASDmp1 and ASDmp2 were significantly distinguishable (Figure 1C–E). The alpha diversity of the gut microbiome in the subjects of ASDmp1 was much greater than that of ASDmp2 (p = 0.0017, Figure 1H) with respect to Pielou's evenness diversity index. Similar results were observed for the other diversity measures: Shannon's entropy (p = 0.016), Faith's phylogenetic diversity (p = 0.04), and the number of observed OTUs in the discovery set (p = 0.05, Figure S2). Among the 12 genera, three displayed a significant differential relative abundance between two subpopulations in both discovery and test sets (Table S1): Parabacteroides (p = 8.03 × 10 − 9 ), Streptococcus (p = 0.03), and Granulicatella (p = 0.021). Interestingly, the opposite abundance pattern of Parabacteriodes and Blautia was observed ( χ 2 = 29.6, p = 5.12 × 10 − 8 , Figure 1F). We then performed MetaCyc Pathway analysis and identified a list of nine differential metabolic pathways (fdr.p < 0.05), in which seven were enriched in ASDmp1, such as nicotinate degradation (fdr.p = 3.44 × 10 − 9 ), glucose degradation (fdr.p = 0.012), and mevalonate (fdr.p = 0.012) pathways. On the contrary, heparin (fdr.p = 0.00039) and l-glutamate degradation (fdr.p = 0.04) were enriched in ASDmp2 (Figure 1I). (the FDR adjusted p value: fdr.p; the permuttion adjusted p valure: perm.p) In parallel, the two subpopulations that were identified by the 12-genus signature also displayed significant differences in the CARS total score (perm. p = 0.001), ADOS total score (perm. p = 0.012), and CARS item scores ≥3 (perm. p = 0.0058) in the discovery set (Table 1). The proportion of subjects with severe symptoms in the ASDmp2 was significantly higher than that of the ASDmp1 (42% versus 77%, odds-ratio = 0.2, p = 0.005). Using the kNN classifier previously trained, the subjects in the test set were also classified into two subpopulations, ASDmp1 (n = 32) and ASDmp2 (n = 18). Similarly, the differences between the two subpopulations in the test set were observed in CARS total score (perm. p = 0.032), ADOS total score (perm.p = 0.043), CARS item score ≥3 (perm.p = 0.012), and microbial profiles (Figure 1C–E). An opposite abundance of Parabacteriodes and Blautia was also observed ( χ 2 = 27.3, p = 1.7 × 10 − 7 , Figure 1G). The proportion of severe subjects in the ASDmp2 was much higher than that in the ASDmp1 (41% versus 78%, odds ratio = 0.2, p = 0.018). Furthermore, the kNN classifier was applied to the two geographically separate sets. Although the microbiota diversity from individuals with a geographical difference was large, the 12-genus signature was able to categorize the subjects into the ASDmp1 and ASDmp2. The Parabacteriodes genus displayed a higher average abundance level in ASDmp1 as compared with that of ASDmp2 (Figure S3), which was consistent with the previous results. We further explored the behavioral measurements in the three primary CARS domains8 representing diverse behavioral and emotional aspects between these two subpopulations. Compared to ASDmp1, the subjects in ASDmp2 displayed markedly higher scores on the CARS domains of social impairment (perm. p = 0.0026) and distorted sensory response (perm. p = 0.048), and also indicated a trend in negative emotionality in the discovery set (perm. p = 0.06) (Table 1). Consistently, the comparison in the test set revealed significant differences in social impairment (perm. p = 0.03) and distorted sensory response (perm. p = 0.009), suggesting that the unique microbial profile is correlated with severity of the disorder at specific behavioral and emotional domains in human. Meanwhile, we sought to determine whether dysbiosis, gastrointestinal dysfunction, and feeding behavior assessments are differentiated between the two subpopulations. By combining the discovery and test sets, we explored dietary and gastrointestinal assessments in 75 subjects (ASDmp1: n = 47 versus ASDmp2: n = 28). Although we observed more constipation in severe ASD patients (ASDmp1:29/47, ASDmp2:20/28, p = 0.45), there was no significant difference in gastrointestinal symptoms between the two subpopulations. Notably, significant differences in feeding behavior “narrowed food spectrum” (p = 0.044) and “resistance to accept new food” (p = 0.031) were identified between the two subpopulations, suggesting that there were indeed more eating problems in patients with severe ASD and perhaps the severity of autistic feeding behavior may lead to altered microbial composition in ASD patients. To better understand how the signature of 12 genera affects the autistic behavior through the microbiota–gut–brain axis, we compared the brain structures between the subjects from ASDmp1 (n = 62) and ASDmp2 (n = 37) by combining the discovery and test sets. Brain volumetric differences between the subjects from the two subpopulations were observed in 12 regions (Figure 2A, Table S2), which spread across primary sensory regions, such as superior occipital gyrus, middle occipital gyrus, superior temporal gyrus, and postcentral gyrus. Significant volumetric differences were also observed in PCC, cuneus cortex, calcarine gyrus, supramarginal gyrus, and superior parietal gyrus. Interestingly, the PCC volume was negatively correlated with the autistic symptoms related to the distorted sensory response (r = −0.268, perm. p = 0.033, Figure 2B) in the ASDmp1. Similarly, a negative correlation was observed in the left superior temporal gyrus (r = −0.287, perm. p = 0.030) and right middle occipital gyrus (r = −0.260, perm. p = 0.049) with autistic sensory scores (Figure 2B). The differential brain image regions between ASD patients from subpopulations. (A) Volume differences between subjects in ASDmp1 and ASDmp2, defined by the 12-genus signature. (B) Association between left PCC, left STG and right MOG volume and CARS sensory score of subjects in ASDmp1 and ASDmp2. Abbreviations: lPCC, - left posterior cingulate cortex; rPCC - right PCC; CUNE-cuneus cortex; SOG, superior occipital gyrus; MOG, middle occipital gyrus; PoCG, postcentral gyrus; STG, superior temporal gyrus; SPG, superior parietal gyrus; SMarg, supramarginal gyrus; CALC, calcarine gyrus; CARS, Childhood Autism Rating Scale Taken together, to our best knowledge this is the first report demonstrating a potential microbial signature associated with brain structure variations and behavioral characteristics in ASD, which would further strengthen our understanding of this severe and heterogeneous neurodevelopmental disorder. However, more research with a large cohort study to clarify the role of the influence of dysbiosis on brain function and autistic behavior is needed in the future. We thank Luis Pedro Coelho (Fudan University) for helpful discussions. This study was funded by the Shanghai Municipal Commission of Health and Family Planning (2017ZZ02026, 2018BR33, 2017EKHWYX-02, and GWV-10.1-XK07), Shanghai Shenkang Hospital Development Center (16CR2025B), Shanghai Municipal Education Commission (20152234), National Natural Science Foundation of China (81571031, 81761128035, 81930095, 81873909, 81703249, 61932008, 61772368, 11925103, and 31860306), Shanghai Committee of Science and Technology (17ZR1444400 and 19410713500), Xinhua Hospital of Shanghai Jiao Tong University School of Medicine (2018YJRC03), National Human Genetic Resources Sharing Service Platform (2005DKA21300), the National Key Research and Development Program of China (2018YFC0910503, 2018YFC0910500, 2020YFA0712403, and 2018YFC0116600), 111 Project (B18015), the Shanghai Municipal Science and Technology Major Project (2018SHZDZX01), Key Scientific and Technological Projects of Guangdong Province (2018B030335001), Collaborative Innovation Program of Shanghai Municipal Health Commission (2020CXJQ01), Shanghai Clinical Key Subject Construction Project (shslczdzk02902), and Science and Technology Department of Yunnan Province (202001AV070010). The authors declare that they have no conflict of interest. This study was approved by the Institutional Review Board of the Xinhua Hospital Shanghai Jiao Tong University (XHEC-C-2019-076). The study was performed in accordance with the Helsinki Declaration and Rules of Good Clinical Practice (GCP). Written informed consents were signed by the caregivers of all applicants after the study protocol was fully explained. Kai Shi, Lingli Zhang, and Juehua Yu designed the study, analyzed the data, and wrote the manuscript. Zilin Chen, Shenying Lai, Wei-Guang Li, and Qiang Luo analyzed the data and helped design experiments. Fei Li and Xing-Ming Zhao designed the study, supervised all work, and helped write the manuscript. Data are available on a reasonable request from the corresponding author. Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.