Articles

Article

Original Article

Exp Neurobiol 2017; 26(5): 307-317

Published online October 31, 2017

https://doi.org/10.5607/en.2017.26.5.307

© The Korean Society for Brain and Neural Sciences

Rapid Assessment of Microbiota Changes in Individuals with Autism Spectrum Disorder Using Bacteria-derived Membrane Vesicles in Urine

Yunjin Lee1, Jin-Young Park1, Eun-Hwa Lee1, Jinho Yang2, Bo-Ri Jeong2,Yoon-Keun Kim2, Ju-Young Seoh3, SoHyun Lee4, Pyung-Lim Han1,5* and Eui-Jung Kim6*

1Departments of Brain and Cognitive Sciences, 4Special Education, and 5Chemistry and Nano Science, Ewha Womans University, Seoul 03760, Korea, 2MD Healthcare Inc., Seoul, Korea; Departments of 3Microbiology and 6Psychiatry, College of Medicine, Ewha Womans University, Seoul 07985, Korea

Correspondence to: *To whom correspondence should be addressed.
Eui-Jung Kim
TEL: 82-2-2650-5163, FAX: 82-2-2650-0984, e-mail: christie@ewha.ac.kr
Pyung-Lim Han
TEL: 82-2-3277-4130, FAX: 82-2-3277-3419, e-mail: plhan@ewha.ac.kr

Received: August 2, 2017; Revised: September 13, 2017; Accepted: September 24, 2017

Individuals with autism spectrum disorder (ASD) have altered gut microbiota, which appears to regulate ASD symptoms via gut microbiota-brain interactions. Rapid assessment of gut microbiota profiles in ASD individuals in varying physiological contexts is important to understanding the role of the microbiota in regulating ASD symptoms. Microbiomes secrete extracellular membrane vesicles (EVs) to communicate with host cells and secreted EVs are widely distributed throughout the body including the blood and urine. In the present study, we investigated whether bacteria-derived EVs in urine are useful for the metagenome analysis of microbiota in ASD individuals. To address this, bacterial DNA was isolated from bacteria-derived EVs in the urine of ASD individuals. Subsequent metagenome analysis indicated markedly altered microbiota profiles at the levels of the phylum, class, order, family, and genus in ASD individuals relative to control subjects. Microbiota identified from urine EVs included gut microbiota reported in previous studies and their up- and down-regulation in ASD individuals were partially consistent with microbiota profiles previously assessed from ASD fecal samples. However, overall microbiota profiles identified in the present study represented a distinctive microbiota landscape for ASD. Particularly, the occupancy of g_Pseudomonas, g_Sphingomonas, g_Agrobacterium, g_Achromobacter, and g_Roseateles decreased in ASD, whereas g_Streptococcus, g_Akkermansia, g_Rhodococcus, and g_Halomonas increased. These results demonstrate distinctively altered gut microbiota profiles in ASD, and validate the utilization of urine EVs for the rapid assessment of microbiota in ASD.

Keywords: Autism spectrum disorder, gut microbiota, Extracellular membrane vesicles, Bacteria-derived EVs, urine marker

Autism spectrum disorder (ASD) is a group of neurodevelopmental disabilities characterized by two domains of core symptoms, persistent social deficits and restricted repetitive patterns of behavior [1]. Most individuals with ASD suffer from various behavioral and physical symptoms, including abnormal preferences regarding specific foods and problems in the digestive system [2,3]. Approximately 70~80% of ASD subjects have food selectivity and restricted food interests due to the texture, smell, or color of specific foods, and food intolerance [3,4]. The limited food intake behavior in ASD subjects leads to health problems including nutrition imbalance and gastrointestinal (GI) symptoms, such as diarrhea and constipation [2,5,6]. Furthermore, studies of a positive correlation between GI symptoms and ASD have been reported; the c-Met promoter variant rs1858830 is associated with ASD and GI symptoms, and the serum level of hepatocyte growth factor (HGF) that binding to the c-Met receptor, is correlated with severity of GI symptom in ASD subjects [7,8,9].

Several lines of evidence indicate that ASD patients have altered microbiota composition in the gut compared to healthy subjects [10,11,12,13,14,15,16,17,18,19,20,21]. The occupancy of the phyla Firmicutes, Fusobacteria, Verrucomicrobia, and Actinobacteria was decreased, whereas Bacteroidetes and Proteobacteria were increased in ASD groups [13,21]. More specifically, in ASD, the genera Bifidobacterum and Akkermansia were found to be decreased in ASD, while Lactobacillus was increased [17,18]. Furthermore, the treatment with a probiotics mix containing Streptococcus (thermophiles), Bifidobacterium breve, B. longum, B. infantis, and Lactobacillus acidophilus, L. plantarum, L. paracasei, and L. delbrueckii (subsp. Bulgaricus) or the transplant of fecal microbiota from healthy subjects to ASD individuals increased overall bacterial diversity and the abundance of Bifidobacterium, Prevotella, and Desulfovibrio among other taxa, and alleviated GI symptoms and ASD core symptoms [22,23]. To date, all available microbiota composition in ASD were mostly assessed from fecal samples [11,13,14,17,18,19,20,21] or directly from the cecum and ileum [15,16]. The fact that some microbiota commonly change in independent studies, and others are not consistently reported (e.g., [13,15,19]), increases the possibility of highly complex dynamics in bodily microbiota composition in ASD individuals under different physiological contexts.

Gram-negative bacteria secrete extracellular membrane vesicles (EVs), also called nanovesicles, to communicate with host cells [24], and are detected in stools, and also in urine and blood serum [25,26,27]. EVs secreted by gram-negative bacteria contain DNA, RNA, proteases, phospholipases, adhesins, toxins, and immunomodulatory compounds. Bacteria-derived EVs are associated with cytotoxicity, bacterial attachment, intercellular DNA transfer, and invasion [24,28,29]. Gram-positive bacteria also produce EVs, which contain peptidoglycan, lipoteichoic acid, virulence proteins, DNA and RNA [30,31]. When bacterial EVs were intraperitoneally injected in mice, EVs were rapidly distributed throughout the body with accumulation in the liver, lung, spleen, and kidney within 3 h [27]. Bacteria-derived EVs in the blood and urine represent the major constituents of microbes in the body, namely the gut microbiota [25,26], and indicate the microbiota that are metabolically or pathologically active [25,27].

In the present study, we investigated bodily microbiota represented by bacterial EVs in the urine of ASD individuals. The results of the present study identify markedly altered microbiota profiles in ASD relative to non-ASD healthy controls and suggest that bacterial EVs in urine can be served as a useful tool for the evaluation of microbiota composition in ASD.

Subjects and urine sample preparation

Individuals who were enrolled at the Ewha Special Education Research Institute (Seoul, Republic of Korea) or Ewha Womans University MokDong Hospital (Seoul, Republic of Korea) were diagnosed according to the DSM-5 diagnostic criteria by a child and adolescent psychiatrist followed by characterization using the Korean Childhood Autism Rating Scale (K-CARS) as described previously [32]. The K-CARS is a well-established scale for the diagnosis of ASD with good agreement with the DSM-5 diagnostic criteria [33]. This questionnaire contained 15 items, each with 4 symptom scales, and all individual scores on each of the questions were summed to obtain the total score. When the total score was higher than 30 points, the subject was classified as autistic. Individuals who had any associated additional psychiatric and neurological diagnoses, or individuals who were on any antipsychotic medications were excluded from the present study.

Among the characterized ASD individuals, 18 male and 2 female ASD individuals (22.4+/−4.9 years) (Table 1) were joined to this study and their urine was collected during the day. The collected urine samples were frozen and stored at −20℃ until use. Age-matched normal healthy subjects (24 males and 4 females, 21.1+/−9.5 years) (Table 1) were selected from the Inje University Haeundae Paik Hospital (IRB No. 1297992-2015-064) and Seoul National University Hospital Healthcare System Gangnam Center (IRB No. 1502-034-647). The control subjects were not related to ASD and had no clinical findings suggestive of gastrointestinal problems or neuropsychiatric disorders. The control subjects of this study had not taken antibiotics, probiotics or prebiotics in the 3 months prior to the sample collection.

The experimental protocol of human subjects was reviewed and approved by the Institutional Review Board of Ewha Womans University Hospital (IRB No. 2015-08-005-002). All eligible participants had been told about the purpose, procedures, risks and benefits of the present study and informed consent was obtained from all ASD subjects.

Isolation of bacteria-derived EVs and DNA extraction from human urine samples

Bacteria EVs were isolated from the urine of ASD individuals following the procedure described previously [25,26]. Briefly, each urine sample was centrifuged at 10,000 × g for 10 min at 4℃. The supernatant was taken and passed through a 0.22-µm membrane filter to eliminate foreign particles. Isolated EVs were dissolved in 100 µl PBS, and quantified on the basis of protein.

Bacterial DNA extraction from prepared EVs was performed as described previously [25,26]. Briefly, isolated EVs (1 µg by protein, each sample) were boiled at 100℃ for 40 min, centrifuged at 13,000 g for 30 min, and the supernatants were collected. Collected samples were then subjected to bacterial DNA extraction using a DNA extraction kit (PowerSoil DNA Isolation Kit, MO BIO, USA) following the manufacturer's instructions, Isolated DNA was quantified by using the QIAxpert system (QIAGEN, Germany).

PCR amplification, library construction, and sequencing of 16S rRNA gene variable regions

Prepared bacterial DNA was used for PCR amplification of the V3-V4 hypervariable regions of the 16S ribosomal RNA genes using the primer set of 16S_V3_F (5′-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGCCTACGGGNGGCWGCAG-3′) and 16S_V4_R (5′-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGGACTACHVGGGTATCTAATCC-3′). The PCR products were used for the construction of 16S rDNA gene libraries following the MiSeq System guidelines (Illumina Inc., San Diego, CA, USA). The 16S rRNA gene libraries for each sample were quantified using QIAxpert (QIAGEN, Germany), pooled at the equimolar ratio, and used for pyrosequencing with the the MiSeq System (Illumina, USA) according to manufacturer's recommendations.

Taxonomic assignments by 16S rRNA genes sequence reads

Obtained raw pyrosequencing reads were filtered on the basis of the barcode and primer sequences using MiSeq Control Software version 1.1.1 (Illumina, USA). Sequence reads were taxonomically assigned using the MDx-Pro ver.1 profiling program (MD Healthcare Inc., Seoul, Korea). Briefly, the quality of sequence reads was retained by controlling an average PHRED score higher than 20 and read length of more than 300 bp. Operational taxonomic units (OTUs) were clustered using CD-HIT sequence clustering algorithms and were assigned using UCLUST [34] and QIIME [35] on the basis of the GreenGenes 8.15.13 16S rRNA sequence database [36]. Based on the sequence similarities, taxonomic assignments were achieved at the following levels: genus, >94% similarity; family, >90% similarity; order, >85% similarity; class, >80% similarity; and phylum, >75% similarity. In cases where clustering was not possible at the genus level due to a lack of sequence information at the database or redundant sequences, the taxon was named based on the higher-level taxonomy with parentheses.

Visualization and principal component analysis (PCA)

Data were normalized to have a mean of 0 and a standard deviation of 1 by linear normalization. PCA and two-dimensional scatter plots with axis of the first and second principal component were calculated and drawn using Matlab 2011a.

Statistical analysis

Two-sample comparisons were performed using Student's t-test. Data clustering of pyrosequencing reads were compared using the χ2 test or t-test, while the comparisons between phylum compositions were tested by Fisher's exact test using GraphPad Prism 6 (San Diego, CA, USA). All data are presented as the mean±SEM, and a statistical difference was accepted at the 5% level.

Metagenome analysis of bodily microbiota in ASD individuals using bacteria-derived EVs in urine

Bacteria-derived EVs were isolated from the urine of 20 ASD individuals and 28 normal healthy subjects. The average age of the control and ASD subjects was 21.1+/−9.5 years and 22.4+/−4.9 years, respectively (Table 1). ASD subjects showed mild impairment of social interaction and stereotypies. The average K-CARS values of these ASD individuals was in the range between 31.5 and 33.5 and IQ values were in the range between 65 and 86. The control group was composed of healthy volunteers who had no medical problems including those related to ASD.

After the extraction of bacterial genomic DNA from the isolated EVs, variable regions of the 16S rRNA genes were amplified by PCR, and the libraries were constructed, as described previously [25,26]. Subsequent DNA sequencing analyses led us to identify over 2,000 operational taxonomic units (OTUs) for ASD and normal individuals. There was no significant difference in the alpha diversity between the two groups (Fig. 1A). Among the identified OTUs, we assigned 30 OTUs at the phylum level, 75 OTUs at the class level, 141 OTUs at the order level, 279 OTUs at the family level, and 619 OTUs at the genus level. Among these OTUs, we primarily focused on OTUs that occupied more than 0.1% of the identified taxons in the following analyses.

Altered microbiota profiles between ASD individuals and healthy subjects

Sequence readings of EVs-based 16S rDNA indicated that the top five members of the phyla p_Proteobacteria, p_Firmicutes, p_Actinobacteria, p_Bacteroidetes, and p_Cyanobacteria comprised 95.2% of the identified OTUs in healthy subjects, whereas these members covered 90.65% of the total OTUs in ASD individuals, suggesting that ASD individuals have altered phyla composition. More specifically, the occupancy of p_Proteobacteria decreased from 49.12 to 35.30%, p_Cyanobacteria decreased from 4.36 to 1.92%, and p_Armatimonadetes decreased from 0.38 to 0.00% in ASD individuals (Fig. 1B~D). In contrast, the occupancy of p_Firmicutes increased from 24.96 to 33.07% and p_Verrucomicrobia increased from 0.58 to 2.37% in ASD.

The microbiota whose occupancy decreased or increased in ASD individuals were further analyzed at the class, order and family levels (Table 2, Supplemental Fig. S1). The decrease of f_Sphingomonadaceae and f_Rhizobiaceae accounted for the major decrease in p_Proteobacteria in ASD. In p_Cyanobacteria, o_Streptophyta decreased from 3.8 to 1.68%. In p_Verrucomicrobia, f_Verrucomicrobiaceae was the major increase (0.52 to 2.35%). In p_Firmicutes, f_Streptococcaceae, f_Clostridiaceae, an unclassified member of o_Clostridiales and f_Eubacteriaceae increased from 3.43 to 8.09% (Table 2).

The members of the genus occupied by more than 0.1% in either control or ASD individuals are summarized in Table 3 and Supplemental Table 1. Overall, 14 members at the genus level were down-regulated in ASD and their total occupancy in ASD dropped from 34.77 to 14.06%. On the contrary, 17 genus members were up-regulated in ASD and their total occupancy in ASD increased from 6.47 to 22.58%.

More specifically, an unclassified member of f_Enterobacteriaceae decreased from 8.85 to 6.08%, g_Pseudomonas decreased from 7.48 to 5.10%, g_Sphingomonas decreased from 4.17 to 0.71%, g_Agrobacterium decreased from 3.83 to 0.11%, an unclassified member of o_Streptophyta decreased from 3.80 to 1.68%, g_Achromobacter decreased from 2.42 to 0.05%, g_Roseateles decreased from 1.12 to 0.02%, and an unclassified member of f_mitochondria decreased from 0.84 to 0.10% (Table 3, Fig. 2).

On the other hand, g_Streptococcus increased from 1.58 to 4.77%, an unclassified member of o_Clostridiales increased from 1.06 to 1.87%, an unclassified member of f_Comamonadaceae increased from 0.92 to 3.79%, an unclassified member of f_S24-7 increased from 0.84 to 2.02%, g_Akkermansia increased from 0.52 to 2.35%, g_Rhodococcus increased from 0.40% to 1.56%, and g_Halomonas increased from 0.12 to 1.72% (Table 3, Fig. 2).

Metagenome analysis of bacterial EVs in urine identifies altered microbiota profiles in ASD

In the present study, we demonstrated that bacteria-derived EVs in urine were useful for the rapid assessment of microbiota profiles in ASD. The metagenome analysis of urine EVs indicated that p_Verrucomicrobia (0.58 to 2.37%, p=0.02) and p_Firmicutes (24.96 to 33.07%, p=0.03) increased in ASD, whereas p_Cyanobacteria (4.36 to 1.92%, p=0.01) and p_Proteobacterium (49.12 to 35.3%, p=0.01) decreased. There was no significant change in p_Bacteroidetes (5.85 to 8.62%, p=0.11) and p_Actinobacterium (10.91 to 11.74%, p=0.56). The altered microbiota compositions identified from urine EVs of ASD were partially consistent with microbiota compositions assessed from fecal samples reported in recent studies. The analyses of fecal microbiota compositions in previous studies reported that p_Actinobacteria, p_Verrucomicrobia and p_Cyanobacteria decreased or tended to decrease in ASD, but there were conflicting results for p_Firmicutes, p_Bacteroidetes and p_Proteobacteria (Table 4) [13,15,16,19].

In p_Firmicutes, g_Streptococcus (1.58 to 4.77%, p=0.02), g_Jeotgalicoccus (0.03 to 0.5%, p<0.01), g_Oscillospira (0.1 to 0.47%, p=0.02), and an unclassified member of f_Clostridiaceae (0.27 to 0.7%, p=0.04) significantly increased in ASD by more than two-fold. The genera g_Akkermansia occupied the greatest proportion of the increased phylum p_Verrucomicrobia (0.52 to 2.35%, p=0.02), whereas in p_Proteobacterium, 12 genera decreased and 5 genera increased in ASD. The genera g_Pseudomonas, g_Citrobacter, an unclassified member of f_Enterobacteriaceae, g_Sphingomonas, g_Agrobacterium, unclassified members of f_Rhizobiaceae, an unclassified member of f_Bradyrhizobiaceae, and an unclassified member of f_mitochondria, g_Roseateles, g_Delftia and g_Achromobacter were decreased (total 30.58 to 12.37%), whereas g_Halomonas, g_Erwinia, g_Rhodobacter, g_Comamonas, and an unclassified member of f_Comamonadaceae were increased (total 1.39 to 6.74%).

At the genus level, the present study identified g_Sphingomonas, g_Agrobacterium, an unclassified member of o_Streptophyta, g_Achromobacter, and g_Roseateles as being decreased in ASD by more than two-fold among the total OTUs having occupancy 0.1% or higher in either healthy control or ASD subjects, and g_Streptococcus, an unclassified member of f_Comamonadaceae, an unclassified member of f_S24-7, g_Akkermansia, g_Rhodococcus, and g_Halomonas were increased in ASD. Particularly, g_Desulfovibrio increased in ASD (0.04 to 0.48%, p<0.01), g_Lactobacillus tended to increase (2.56 to 5.45%, p=0.08), g_Bifidobacterium tended to decrease (1.9 to 0.8%, p=0.06), and g_Turicibacter tended to decrease (0.19 to 0.02%, p=0.07) in ASD, whereas g_Enterococcus (0.4 to 0.78%, p=0.06), g_Enterobacter (0.03 to 0.37, p=0.08), and g_Clostridium (0.11 to 0.31%, p=0.09) tended to increase in ASD. The changes in these genus members are broadly consistent with the results of previous reports assessed for fecal microbiota compositions (Table 4) [11,13,17,18,21].

The EV levels of g_Oscillospira, unclassified members of f_Clostridiaceae and f_Eubacteriaceae, and an unclassified member of o_Clostridiales were increased in ASD subjects. Previous studies have reported that several species of c_Clostridia produced 4-ethylphenyl sulfate (4-EPS) and p-cresol, which were found at high concentrations in the urine of ASD children. Administration of 4-EPS in healthy mice produced myelination deficits in the prefrontal cortex and sociability defects [37,38,39]. It was also reported that c_Clostridia produced propionic acid (PPA) and its related short-chain fatty acids (SCFAs) as fermentation products, and PPA infusions in rats induced ASD-linked neurochemical and behavioral changes [40]. These results suggest that bacteria-derived metabolites induce neurochemical and structural changes and shape behavioural abnormalities.

Oral treatment with Bifidobacteria fragilis ameliorated ASD-related gastrointestinal deficits and associated behavioural abnormalities behavioral abnormalities in the poly (I:C)-injection model [41]. Bifidobacteria infantis attenuated pro-inflammatory immune responses and production of serotonergic precursor, tryptophan, and has potential anti-depressant properties [42,43]. Considering these results, ASD groups with decreased EV levels of g_Bifidobacterium might have benefits by probiotic treatment with g_Bifidobacterium.

Bacterial EVs in urine are useful for rapid assessment of bodily microbiota profiles in ASD

The microbiota profile assessed from urine EVs might reflect a large part of the gut microbiota. Nonetheless, we do believe that the microbiota profile assessed from urine EVs is not likely a simple alternative for microbiota profile assessed from stool. Possible sources for metagenome analysis of bodily microbiota may include stool bacteria, stool EVs, gut (ex, stomach and/or specific regions of the small and large intestines) bacteria, respiratory exhale EVs, oral/nasal bacteria and EVs, urinary system bacteria and EVs, and blood EVs. Generally speaking, microbiota in stool represents the intestinal compartment, whereas microbiota in urine or blood reflects the whole body including the intestinal compartments, oral system, respiratory system, and urinary system. Nonetheless, among the body parts, the gut is the major source of bodily microbiota. It was reported that the metabolites of intestinal microbiota activities, including phenyllactate, p-cresol sulfate, concentrations, and serotonin in urine, plasma, and stool of mouse pups undernourished by timed separation from lactating dams, then resumed ad libitum nursing, were different from each other, although they had some correlations [44,45]. Similar to metabolite profiles of intestinal microbiota activities, available information suggests that metagenome analysis assessed from these sources might be closely related, but represent some distinct landscapes. For an example, metagenome analysis of bacteria and bacteria-derived EV in stool of inflammatory bowel disease model mice indicated that the EV composition in stool was more drastically altered compared to that of bacterial composition in stool [25]. Considering that bacteria-derived EV indicates the metabolically or pathologically activated microbiota [25,27], urine EV may be more representative of the host's microbiota activities than stool bacteria.

To the best of our knowledge, this is the first report characterizing microbiota in ASD individuals on the basis of urine EVs. Compared to blood and feces, urine is easily obtained in large volume and is readily available via a non-invasive method. Considering the general difficulties in repeated sampling microbiota sources from ASD individuals, particularly low functioning individuals with ASD or toddlers with ASD, using urine as a sample source would be a great advantage for rapid and repeated assessments of microbiota changes under varying physiological contexts compared to the use of blood and feces. Comparative analyses of EV profiles from urine, blood and stool of ASD individuals will be valuable. Also it will be worth to understand EV profiles of ASD with diverse factors including age, sex, familial history, genetics, and ethnics.

Overall, the present study assessed urine EVs from individuals with mildly autistic subjects. We believe that further systematic and unbiased analyses of male and female subjects with broad ASD spectrums are necessary. This study focused on young adult subjects. Considering that ASD should be diagnosed in young children as early as 1.5~3 yr of age, this analysis should be expanded to toddlers and infants.

Supplemental Figure S1

The diversity of microbiota at the class, order, family, and genus levels in control vs. ASD subjects. (A~D) Principal component analysis (PCA) of microbiota diversity at the class (A), order (B), family (C), and genus (D) levels based on the weighted UniFrac distance and Bray-Curtis dissimilarity. Data were normalized to have a mean of 0 and a standard deviation of 1. Control (blue) and ASD (red).

en-26-307-s001.pdf

Supplemental Table 1

The percent composition of microbiota at the genus level in control and ASD subjects

en-26-307-s002.pdf
Fig. 1. The diversity and percent composition of microbiota at the phylum level in control vs. ASD subjects. (A) Rarefication curves representing the mean OTUs over the identified sequences of variable regions of 16S rRNA gene in control (blue) and ASD (red) subjects. Data are the mean +/− SEM (n=5, each). (B) Principal component analysis of microbiota diversity based on the weighted UniFrac distance and Bray-Curtis dissimilarity. Data were normalized to have a mean of 0 and a standard deviation of 1. Control (blue) and ASD (red). (C) Overall composition of microbiota at the phylum level in control (blue) and ASD (red) subjects. Those with occupancy 0.1% or higher in control and/or ASD subjects are presented. (D) The percent composition of microbiota at the phylum level in control and ASD subjects. ↑ and ↓ denote an increase and decrease in the percent composition, respectively. Data are the mean +/− SEM (n=5, each). * and ** denote the differences between the indicated groups at p<0.05 and p<0.01, respectively (Student's t-test).
Fig. 2. Genus members that were downregulated or upregulated in ASD subjects. The percent composition of microbiota whose occupancy was significantly changed at the genus level among those occupying 0.1% or higher in either control (blue) or ASD (red) subjects. Data are presented as the mean percentage±SEM. * and ** denote the differences between the indicated groups at p<0.05 and p<0.01, respectively (Student's t-test).
Table. 1. The number, age, and sex of control and ASD subjects
ControlASDp-value
Age (years)21.1±9.522.4±4.90.57
N2820-
Sex (Male/Female)24/418/20.66

Table. 2. The percent composition of microbiota at the class, order, and family taxonomic levels in control and ASD subjects
ClassOrderFamily
TaxonMean (%)Fold changep-valueTaxonMean (%)Fold changep-valueTaxonMean (%)Fold changep-value
ControlASDControlASDControlASD
Gammaproteobacteria23.5123.4710.99Oceanospirillales0.212.6212.680.02*Halomonadaceae0.22.61↑ 13.160.02*
Alphaproteobacteria15.294.840.320**Sphingomonadales6.291.980.310**Sphingomonadaceae6.271.98↓ 0.320**
Rhizobiales5.691.050.190**Rhizobiaceae4.580.18↓ 0.040**
Bradyrhizobiaceae0.460.08↓ 0.170.01*
Rickettsiales0.940.10.110.04*mitochondria0.840.1↓ 0.120.02**
Betaproteobacteria10.096.260.620.1Burkholderiales9.55.730.60.11Comamonadaceae2.554.37↑ 1.710.01*
Alcaligenaceae2.530.12↓ 0.050**
Deltaproteobacteria0.20.643.230.01*Desulfovibrionales0.060.539.660**Desulfovibrionaceae0.060.53↑ 9.660**
Bacilli12.8417.551.370.07Lactobacillales7.9211.941.510.07Streptococcaceae1.934.88↑ 2.520.04*
Clostridia11.6315.41.320.18Clostridiales11.5915.381.330.18Unclassified1.061.87↑ 1.750.04*
Clostridiaceae0.441.19↑ 2.700.02*
Eubacteriaceae00.15↑ 146.230.04*
Actinobacteria9.911.091.120.38Actinomycetales7.8410.291.310.09Nocardiaceae0.41.56↑ 3.910.02*
Bifidobacteriales2.060.80.390.03*Bifidobacteriaceae2.060.8↓ 0.390.03*
Bacteroidia5.057.281.440.21Bacteroidales5.057.281.440.21S24-70.842.02↑ 2.400.04*
Flavobacteriia0.561.192.110.03*Flavobacteriales0.561.192.110.03*[Weeksellaceae]0.431.12↑ 2.630.01*
Chloroplast3.811.760.460.03*Streptophyta3.81.680.440.03*Unclassified3.81.68↓ 0.440.03*
Verrucomicrobiae0.522.354.490.02*Verrucomicrobiales0.522.354.490.02*Verrucomicrobiaceae0.522.35↑ 4.490.02*
[Fimbriimonadia]0.38000**[Fimbriimonadales]0.38000**[Fimbriimonadaceae]0.380↓ 0.000**
Deinococci0.070.476.470**Thermales0.030.216.30.02*Thermaceae0.030.21↓ 6.300.02*
TM7-30.050.376.960.05*

Microbiota at the family level whose occupancy was significantly different in ASD subjects are presented with associated higher taxonomy levels. Those occupying 0.1% or higher in either normal healthy or ASD subjects were included. ↑ and ↓ denote an increase and decrease in the percent composition, respectively. Data are the mean percentages. * and ** denote significant differences between the indicated groups at p<0.05 and p<0.01, respectively (Student's t-test)..


Table. 3. The percent composition of microbiota at the genus level in control and ASD subjects
ClassOrderFamilyTaxonMean±SEM (%)Fold changep-value
ControlASD
GammaproteobacteriaOceanospirillalesHalomonadaceaeHalomonas0.12±0.061.72±0.51↑ 14.610.01*
PseudomonadalesPseudomonadaceaePseudomonas7.48±0.865.10±0.6↓ 0.680.03*
EnterobacterialesEnterobacteriaceaeErwinia0.26±0.10.64±0.15↑ 2.410.04*
Citrobacter0.66±0.240.08±0.05↓ 0.120.02*
Unclassified8.85±1.016.08±0.58↓ 0.690.02*
AlphaproteobacteriaSphingomonadalesSphingomonadaceaeSphingomonas4.17±0.830.71±0.2↓ 0.170.00**
RhizobialesRhizobiaceaeAgrobacterium3.83±0.960.11±0.05↓ 0.030.00**
Unclassified 10.63±0.170.07±0.03↓ 0.110.00**
Unclassified 20.11±0.030.00±0.00↓ 0.000.00*
BradyrhizobiaceaeUnclassified0.24±0.070.05±0.02↓ 0.210.02*
RickettsialesMitochondriaUnclassified0.84±0.30.1±0.06↓ 0.120.02*
BetaproteobacteriaRhodobacteralesRhodobacteraceaeRhodobacter0.01±0.010.23±0.08↑ 20.340.02*
BurkholderialesComamonadaceaeRoseateles1.12±0.350.02±0.01↓ 0.020.00**
Delftia0.22±0.080.01±0.00↓ 0.050.02*
Comamonas0.08±0.050.36±0.12↑ 4.480.04*
Unclassified0.92±0.153.79±0.49↑ 4.140.00**
AlcaligenaceaeAchromobacter2.42±0.750.05±0.03↓ 0.020.00**
DeltaproteobacteriaDesulfovibrionalesDesulfovibrionaceaeDesulfovibrio0.04±0.020.48±0.14↑ 10.880.00*
BacilliLactobacillalesStreptococcaceaeStreptococcus1.58±0.264.77±1.28↑ 3.030.02*
BacillalesStaphylococcaceaeJeotgalicoccus0.03±0.020.50±0.11↑ 14.700.00*
ClostridiaClostridialesUnclassifiedUnclassified1.06±0.261.87±0.26↑ 1.750.04*
ClostridiaceaeUnclassified0.27±0.080.70±0.18↑ 2.580.04*
RuminococcaceaeOscillospira0.10±0.040.47±0.14↑ 4.530.02*
ActinobacteriaActinomycetalesNocardiaceaeRhodococcus0.40±0.131.56±0.43↑ 3.910.02*
MicrococcaceaeKocuria0.06±0.050.30±0.08↑ 4.750.01*
BacteroidiaBacteroidalesS24-7Unclassified0.84±0.362.02±0.44↑ 2.400.04*
FlavobacteriiaFlavobacteriales[WeeksellaceaeCloacibacterium0.13±0.060.62±0.2↑ 4.920.03*
ChloroplastStreptophytaUnclassifiedUnclassified3.80±0.841.68±0.38↓ 0.440.03*
VerrucomicrobiaeVerrucomicrobialesVerrucomicrobiaceaeAkkermansia0.52±0.182.35±0.68↑ 4.520.02*
[Fimbriimonadia][Fimbriimonadales][Fimbriimonadaceae]Fimbriimonas0.38±0.120.00±0.00↓ 0.000.00**
DeinococciThermalesThermaceaeThermus0.03±0.020.21±0.07↑ 6.300.02*

Microbiota at the genus level whose occupancy was significantly different in ASD subjects are presented with associated higher taxonomy levels. Microbiota with occupancy 0.1% or higher in either normal healthy or ASD subjects were considered. ↑ and ↓ denote an increase and decrease in the percent composition, respectively..

Data are the mean percentage±SEM. * and ** denote significant differences between the indicated groups at p<0.05 and p<0.01, respectively (Student's t-test)..


Table. 4. Summary and comparison of microbiota characterized in the present study with those identified from fecal samples in previous studies
TaxonsMean (%)Fold changep-valueLiteratures
ControlASD
PhylumProteobacteria49.1235.30↓ 0.720.01**↓ [19]; ↑ [13]
Firmicutes24.9633.07↑ 1.330.03*↓ [13]; ↑ [15]
Actinobacteria10.9111.74↑ 1.080.56↓ [13]
Bacteroidetes5.858.62↑ 1.470.11↑ [13]; ↓ [15]; ↓ [19]
Verrucomicrobia0.582.37↑ 4.120.02*↓ [19]
ClassBetaproteobacteria10.096.26↓ 0.620.10↑ [15]
OrderClostridiales11.5915.38↑ 1.330.18↑ [15]
FamilyRuminococcaceae5.466.11↑ 1.120.68↑ [15]
Lachnospiraceae3.122.57↓ 0.820.54↑ [15]
Corynebacterium2.603.38↑ 1.30.34↑ [11]
Alcaligenaceae2.530.12↓ 0.050.00**↑ [16]
Pseudomonas7.485.10↓ 0.680.03*↑ [21]
Lactobacillus2.565.45↑ 2.130.08↑ [18]; ↑ [11]
Bacteroides2.482.93↑ 1.180.65↑ [13]; ↑ [21]
Staphylococcus2.232.57↑ 1.150.66↓ [21]
Faecalibacterium2.142.05↓ 0.960.90↑ [20]
Bifidobacterium1.900.80↓ 0.420.06↓ [13]; ↓ [18]; ↓ [21]; ↓ [17]
Streptococcus1.584.77↑ 3.030.02*↓ [13]; ↓ [21]
Akkermansia0.522.35↑ 4.520.02*↓ [19]; ↓ [17]
Blautia0.470.28↓ 0.590.38↓ [20]
Enterococcus0.400.78↑ 1.920.06↑ [21]; ↓ [18]
Collinsella0.370.11↓ 0.290.16↓ [13]; ↑ [11]
Veillonella0.370.50↑ 1.340.43↓ [11]
Lactococcus0.350.11↓ 0.30.16↓ [13]; ↓ [21]
[Ruminococcus]0.340.15↓ 0.450.16↓ [13]
Coprococcus0.280.27↓ 0.940.90↓ [19]
Leuconostoc0.270.01↓ 0.020.14↓ [13]
Dialister0.231.39↑ 6.160.18↓ [13]; ↓ [11]
Parabacteroides0.210.18↓ 0.830.76↓ [11]; ↑ [13]
Weissella0.210.07↓ 0.310.33↓ [13]
Turicibacter0.190.02↓ 0.120.07↓ [13]
Dorea0.170.10↓ 0.60.44↑ [11]
Clostridium0.110.31↑ 2.720.09↑ [21]; ↓ [13]
[Prevotella]0.090.12↑ 1.340.74↑ [21]; ↓ [19]
Desulfovibrio0.040.48↑ 10.880.00**↑ [13]; ↓ [19]
Genus0.030.37↑ 11.530.08↑ [21]

The microbiota whose percent composition were significantly different in ASD subjects as characterized in the present study were compared with those identified from fecal samples in previous studies. ↑ and ↓ denote an increase and decrease in the percent composition, respectively. * and ** denote significant differences between the indicated groups at p<0.05 and p<0.01, respectively (Student's t-test)..


  1. American Psychiatric Association. Diagnostic and statistical manual of mental disorders: DSM-5. 5th ed. Washington, D.C.: American Psychiatric Association, 2013.
  2. Kral TV, Eriksen WT, Souders MC, Pinto-Martin JA. Eating behaviors, diet quality, and gastrointestinal symptoms in children with autism spectrum disorders: a brief review. J Pediatr Nurs 2013;28:548-556.
    Pubmed
  3. Marí-Bauset S, Zazpe I, Mari-Sanchis A, Llopis-González A, Morales-Suárez-Varela M. Food selectivity in autism spectrum disorders: a systematic review. J Child Neurol 2014;29:1554-1561.
    Pubmed
  4. Bryant-Waugh R, Markham L, Kreipe RE, Walsh BT. Feeding and eating disorders in childhood. Int J Eat Disord 2010;43:98-111.
    Pubmed
  5. Hyman SL, Stewart PA, Schmidt B, Cain U, Lemcke N, Foley JT, Peck R, Clemons T, Reynolds A, Johnson C, Handen B, James SJ, Courtney PM, Molloy C, Ng PK. Nutrient intake from food in children with autism. Pediatrics 2012;130:S145-S153.
    Pubmed
  6. McElhanon BO, McCracken C, Karpen S, Sharp WG. Gastrointestinal symptoms in autism spectrum disorder: a meta-analysis. Pediatrics 2014;133:872-883.
    Pubmed
  7. Campbell DB, Buie TM, Winter H, Bauman M, Sutcliffe JS, Perrin JM, Levitt P. Distinct genetic risk based on association of MET in families with co-occurring autism and gastrointestinal conditions. Pediatrics 2009;123:1018-1024.
    Pubmed
  8. Russo AJ, Krigsman A, Jepson B, Wakefield A. Decreased serum hepatocyte growth factor (HGF) in autistic children with severe gastrointestinal disease. Biomark Insights 2009;4:181-190.
    Pubmed
  9. Hsiao EY. Gastrointestinal issues in autism spectrum disorder. Harv Rev Psychiatry 2014;22:104-111.
    Pubmed
  10. Mulle JG, Sharp WG, Cubells JF. The gut microbiome: a new frontier in autism research. Curr Psychiatry Rep 2013;15:337.
    Pubmed
  11. Strati F, Cavalieri D, Albanese D, De Felice C, Donati C, Hayek J, Jousson O, Leoncini S, Renzi D, Calabrò A, De Filippo C. New evidences on the altered gut microbiota in autism spectrum disorders. Microbiome 2017;5:24.
    Pubmed
  12. Rosenfeld CS. Microbiome disturbances and autism spectrum disorders. Drug Metab Dispos 2015;43:1557-1571.
    Pubmed
  13. Finegold SM, Dowd SE, Gontcharova V, Liu C, Henley KE, Wolcott RD, Youn E, Summanen PH, Granpeesheh D, Dixon D, Liu M, Molitoris DR, Green JA. Pyrosequencing study of fecal microflora of autistic and control children. Anaerobe 2010;16:444-453.
    Pubmed
  14. Tomova A, Husarova V, Lakatosova S, Bakos J, Vlkova B, Babinska K, Ostatnikova D. Gastrointestinal microbiota in children with autism in Slovakia. Physiol Behav 2015;138:179-187.
    Pubmed
  15. Williams BL, Hornig M, Buie T, Bauman ML, Cho Paik M, Wick I, Bennett A, Jabado O, Hirschberg DL, Lipkin WI. Impaired carbohydrate digestion and transport and mucosal dysbiosis in the intestines of children with autism and gastrointestinal disturbances. PLoS One 2011;6:e24585.
    Pubmed
  16. Williams BL, Hornig M, Parekh T, Lipkin WI.. Application of novel PCR-based methods for detection, quantitation, and phylogenetic characterization of Sutterella species in intestinal biopsy samples from children with autism and gastrointestinal disturbances. MBio 2012;3:e00261-e00211.
    Pubmed
  17. Wang L, Christophersen CT, Sorich MJ, Gerber JP, Angley MT, Conlon MA. Low relative abundances of the mucolytic bacterium Akkermansia muciniphila and Bifidobacterium spp. in feces of children with autism. Appl Environ Microbiol 2011;77:6718-6721.
    Pubmed
  18. Adams JB, Johansen LJ, Powell LD, Quig D, Rubin RA. Gastrointestinal flora and gastrointestinal status in children with autism--comparisons to typical children and correlation with autism severity. BMC Gastroenterol 2011;11:22.
    Pubmed
  19. Kang DW, Park JG, Ilhan ZE, Wallstrom G, Labaer J, Adams JB, Krajmalnik-Brown R. Reduced incidence of Prevotella and other fermenters in intestinal microflora of autistic children. PLoS One 2013;8:e68322.
    Pubmed
  20. Inoue R, Sakaue Y, Sawai C, Sawai T, Ozeki M, Romero-Pérez GA, Tsukahara T. A preliminary investigation on the relationship between gut microbiota and gene expressions in peripheral mononuclear cells of infants with autism spectrum disorders. Biosci Biotechnol Biochem 2016;80:2450-2458.
    Pubmed
  21. De Angelis M, Piccolo M, Vannini L, Siragusa S, De Giacomo A, Serrazzanetti DI, Cristofori F, Guerzoni ME, Gobbetti M, Francavilla R. Fecal microbiota and metabolome of children with autism and pervasive developmental disorder not otherwise specified. PLoS One 2013;8:e76993.
    Pubmed
  22. Santocchi E, Guiducci L, Fulceri F, Billeci L, Buzzigoli E, Apicella F, Calderoni S, Grossi E, Morales MA, Muratori F. Gut to brain interaction in Autism Spectrum Disorders: a randomized controlled trial on the role of probiotics on clinical, biochemical and neurophysiological parameters. BMC Psychiatry 2016;16:183.
    Pubmed
  23. Kang DW, Adams JB, Gregory AC, Borody T, Chittick L, Fasano A, Khoruts A, Geis E, Maldonado J, McDonough-Means S, Pollard EL, Roux S, Sadowsky MJ, Lipson KS, Sullivan MB, Caporaso JG, Krajmalnik-Brown R. Microbiota Transfer Therapy alters gut ecosystem and improves gastrointestinal and autism symptoms: an open-label study. Microbiome 2017;5:10.
    Pubmed
  24. Kuehn MJ, Kesty NC. Bacterial outer membrane vesicles and the host-pathogen interaction. Genes Dev 2005;19:2645-2655.
    Pubmed
  25. Kang CS, Ban M, Choi EJ, Moon HG, Jeon JS, Kim DK, Park SK, Jeon SG, Roh TY, Myung SJ, Gho YS, Kim JG, Kim YK. Extracellular vesicles derived from gut microbiota, especially Akkermansia muciniphila, protect the progression of dextran sulfate sodium-induced colitis. PLoS One 2013;8:e76520.
    Pubmed
  26. Yoo JY, Rho M, You YA, Kwon EJ, Kim MH, Kym S, Jee YK, Kim YK, Kim YJ. 16S rRNA gene-based metagenomic analysis reveals differences in bacteria-derived extracellular vesicles in the urine of pregnant and non-pregnant women. Exp Mol Med 2016;48:e208.
    Pubmed
  27. Jang SC, Kim SR, Yoon YJ, Park KS, Kim JH, Lee J, Kim OY, Choi EJ, Kim DK, Choi DS, Kim YK, Park J, Di Vizio D, Gho YS. In vivo kinetic biodistribution of nano-sized outer membrane vesicles derived from bacteria. Small 2015;11:456-461.
    Pubmed
  28. Ellis TN, Kuehn MJ. Virulence and immunomodulatory roles of bacterial outer membrane vesicles. Microbiol Mol Biol Rev 2010;74:81-94.
    Pubmed
  29. Maldonado R, Wei R, Kachlany SC, Kazi M, Balashova NV. Cytotoxic effects of Kingella kingae outer membrane vesicles on human cells. Microb Pathog 2011;51:22-30.
    Pubmed
  30. Lee EY, Choi DY, Kim DK, Kim JW, Park JO, Kim S, Kim SH, Desiderio DM, Kim YK, Kim KP, Gho YS. Grampositive bacteria produce membrane vesicles: proteomics-based characterization of Staphylococcus aureus-derived membrane vesicles. Proteomics 2009;9:5425-5436.
    Pubmed
  31. Brown L, Kessler A, Cabezas-Sanchez P, Luque-Garcia JL, Casadevall A. Extracellular vesicles produced by the Gram-positive bacterium Bacillus subtilis are disrupted by the lipopeptide surfactin. Mol Microbiol 2014;93:183-198.
    Pubmed
  32. Shim MS, Kim YH. Standardization study for the Korean version of childhood autism rating scale: reliability, validity and cut-off score. Korean J Clin Psychol 1998;17:1-15.
  33. Park EY, Kim J. Factor structure of the Childhood Autism Rating Scale as per DSM-5. Pediatr Int 2016;58:139-145.
    Pubmed
  34. Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010;26:2460-2461.
    Pubmed
  35. Lozupone C, Hamady M, Knight R. UniFrac--an online tool for comparing microbial community diversity in a phylogenetic context. BMC Bioinformatics 2006;7:371.
    Pubmed
  36. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Peña AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R. QIIME allows analysis of high-throughput community sequencing data. Nat Methods 2010;7:335-336.
    Pubmed
  37. Altieri L, Neri C, Sacco R, Curatolo P, Benvenuto A, Muratori F, Santocchi E, Bravaccio C, Lenti C, Saccani M, Rigardetto R, Gandione M, Urbani A, Persico AM. Urinary p-cresol is elevated in small children with severe autism spectrum disorder. Biomarkers 2011;16:252-260.
    Pubmed
  38. Gabriele S, Sacco R, Cerullo S, Neri C, Urbani A, Tripi G, Malvy J, Barthelemy C, Bonnet-Brihault F, Persico AM. Urinary p-cresol is elevated in young French children with autism spectrum disorder: a replication study. Biomarkers 2014;19:463-470.
    Pubmed
  39. Gacias M, Gaspari S, Santos PM, Tamburini S, Andrade M, Zhang F, Shen N, Tolstikov V, Kiebish MA, Dupree JL, Zachariou V, Clemente JC, Casaccia P. Microbiota-driven transcriptional changes in prefrontal cortex override genetic differences in social behavior. Elife 2016;5:e13442.
    Pubmed
  40. Macfabe DF. Short-chain fatty acid fermentation products of the gut microbiome: implications in autism spectrum disorders. Microb Ecol Health Dis 2012;23:19260.
  41. Hsiao EY, McBride SW, Hsien S, Sharon G, Hyde ER, McCue T, Codelli JA, Chow J, Reisman SE, Petrosino JF, Patterson PH, Mazmanian SK. Microbiota modulate behavioral and physiological abnormalities associated with neurodevelopmental disorders. Cell 2013;155:1451-1463.
    Pubmed
  42. Desbonnet L, Garrett L, Clarke G, Bienenstock J, Dinan TG. The probiotic Bifidobacteria infantis: an assessment of potential antidepressant properties in the rat. J Psychiatr Res 2008;43:164-174.
    Pubmed
  43. Desbonnet L, Garrett L, Clarke G, Kiely B, Cryan JF, Dinan TG. Effects of the probiotic Bifidobacterium infantis in the maternal separation model of depression. Neuroscience 2010;170:1179-1188.
    Pubmed
  44. Preidis GA, Ajami NJ, Wong MC, Bessard BC, Conner ME, Petrosino JF. Microbial-derived metabolites reflect an altered intestinal microbiota during catch-up growth in undernourished neonatal mice. J Nutr 2016;146:940-948.
    Pubmed
  45. Armstrong CW, McGregor NR, Lewis DP, Butt HL, Gooley PR. The association of fecal microbiota and fecal, blood serum and urine metabolites in myalgic encephalomyelitis/chronic fatigue syndrome. Metabolomics 2017;13:8.