Original Research ARTICLE
Gut Microbiota Dysbiosis Associated With Altered Production of Short Chain Fatty Acids in Children With Neurodevelopmental Disorders
- 1Psychiatry Clinic “Dr Selaković”, Belgrade, Serbia
- 2Department of Biochemistry, Institute for Biological Research “Siniša Stanković”, University of Belgrade, Belgrade, Serbia
- 3Laboratory for Molecular Microbiology, Institute of Molecular Genetics and Genetic Engineering, University of Belgrade, Belgrade, Serbia
While gut microbiota dysbiosis has been linked with autism, its role in the etiology of other neurodevelopmental disorders (NDD) is largely underexplored. To our knowledge this is the first study to evaluate gut microbiota diversity and composition in 36 children from the Republic of Serbia diagnosed with NDD and 28 healthy children. The results revealed an increased incidence of potentially harmful bacteria, closely related to Clostridium species, in the NDD patient group compared to the Control group: Desulfotomaculum guttoideum (P < 0.01), Intestinibacter bartlettii (P < 0.05), and Romboutsia ilealis (P < 0.001). On the other hand, significantly lower diversity of common commensal bacteria in the NDD group of patients was noticed. Enterococcus faecalis (P < 0.05), Enterococcus gallinarum (P < 0.01), Streptococcus pasteurianus (P < 0.05), Lactobacillus rhamnosus (P < 0.01) and Bifidobacteria sp. were detected in lower numbers of patients or were even absent in some NDD patients. In addition, butyrate-producing bacteria Faecalibacterium prausnitzii (P < 0.01), Butyricicoccus pullicaecorum (P < 0.05), and Eubacterium rectale (P = 0.07) were less frequent in the NDD patient group. In line with that, the levels of fecal short chain fatty acids (SCFAs) were determined. Although significant differences in SCFA levels were not detected between NDD patients and the Control group, a positive correlation was noted between number of rDNA amplicons obtained with universal primers and level of propionic acid, as well as a trend for levels of total SCFAs and butyric acid in the Control group. This correlation is lost in the NDD patient group, indicating that NDD patients' microbiota differs from the microbiota of healthy children in the presence or number of strong SCFA-producing bacteria. According to a range-weighted richness index it was observed that microbial diversity was significantly lower in the NDD patient group. Our study reveals that the intestinal microbiota from NDD patients differs from the microbiota of healthy children. It is hypothesized that early life microbiome might have an impact on GI disturbances and accompanied behavioral problems frequently observed in patients with a broad spectrum of NDD.
Neurodevelopmental disorders (NDD) [according to Diagnostic and Statistical Manual of Mental Disorders 5th Edition (DSM-5)] (Swedo et al., 2013) or Disorders of psychological development [according to International Classification of Diseases 10th Edition (ICD10)] (WHO, 2015) are a group of disabilities that occur early in childhood, usually at preschool age. Children with NDDs generally have some degree of speech-language pathology, sensory-motor disorders, specific problems in learning and memory and socio-emotional functioning. Autism, which is the most serious neurodevelopmental condition, has attracted the most attention from the scientific community so far. Although the causes of autism are not yet completely understood, it has been suggested that interactions between some genes and environmental factors are needed for full appearance of the disorder (Muhle et al., 2004). The causes of neurodevelopmental disorders other than autism have not been intensively investigated so far.
Lately, co-morbidities, especially gastrointestinal (GI) disturbances, have been recognized as potential risk factors in the development of autism and other NDDs. It has been observed that individuals with autism and other developmental delays are frequently affected by GI disorders like diarrhea, constipation, bloating and gastro-esophageal reflux (Schieve et al., 2012; Chaidez et al., 2014) and that the prevalence of these GI disturbances is higher in children with some developmental disabilities than in children with typical development (Schieve et al., 2012). While it was noticed that GI complications correlate with the severity of behavioral abnormalities, it has been suggested that these co-morbidities could contribute to the manifestation of autism-related behaviors (Horvath and Perman, 2002; Nikolov et al., 2009; Adams et al., 2011; Hsiao, 2014; Tomova et al., 2015).
Many authors have assumed that GI disturbances detected in patients with autism might be linked to an abnormal composition of the gut microbiota and have proposed a connection between the disturbed composition of gut microbiota and autism (Mulle et al., 2013; Borre et al., 2014; De Angelis et al., 2015; Frye et al., 2015; Reddy and Saier, 2015). Several studies have revealed overgrowth of some enteric bacteria, particularly bacteria belonging to Clostridia clusters, in children with autism spectrum disorder (ASD) (Finegold et al., 2002; Song Y. L. et al., 2004; Parracho et al., 2005; De Angelis et al., 2013). On the other hand, an imbalance in beneficial bacteria, decrease in Bifidobacteria and increase in Lactobacilli, was noticed in fecal samples of children with autism (Adams et al., 2011). Although results differed between studies, making it difficult to state what is the “normal” human microbial diversity, the general finding is that microbiota from patients with autism are less diverse and that the balance in beneficial bacteria was disturbed (Krajmalnik-Brown et al., 2015). The significant role of microbiota in CNS functioning and neurodevelopment due to bidirectional communication in the microbiota-gut-brain axis was revealed by Sudo et al. (2004) and Neufeld K. M. et al. (2011). It is postulated that perturbations of initial microbiota colonization can lead to neurochemical disturbances (Borre et al., 2014) through changes in levels of various neuroactive bacterial metabolites, such as serotonin, dopamine, γ-aminobutyric acid, and short chain fatty acids (SCFAs) (Rogers et al., 2016).
Clinical reports have indicated increased incidence of GI disturbances and GI dysbiosis in NDD patients, as well as the improvement of symptoms and the acceleration of recovery after dysbiosis treatment. These observations together with the lack of knowledge about the causes of NDD inspired us to perform the following study. The aim of this study was to evaluate gut microbiota diversity and to identify bacterial strains whose prevalence were statistically different between NDD patients and the Control group, as well as to determine the levels of SCFAs as an indirect measure of balance among the enteric bacteria (Krajmalnik-Brown et al., 2015). The impact of any exclusion diet, such as a gluten-free, casein-free (GFCF) diet, or any medical treatment was excluded since all patients were selected on the first examination. To our knowledge this is the first study of gut microbiota performed on samples from a wide spectrum of NDD patients from Serbia.
Materials and Methods
Thirty-six pre-adolescent children from the Republic of Serbia diagnosed with NDDs participated in the study. Specifically, 9 patients with expressive and receptive language disorder (ERLD, F80.1 and F80.2), 12 patients with mixed specific developmental disorder (MSDD, F83.0), 10 patients with pervasive developmental disorder non-specified (PDD-NOS, F84.9) and 5 patients with childhood autism (CHA, F84), in total 28 males and eight females between 2 and 9 years of age (median ± SEM, 4.22 ± 0.30 y), participated in the study. The diagnosis was given by a medical doctor specializing in child neuropsychiatry according to the observation of clinical manifestations and anamnesis data given by parents, based on ICD10 criteria, which provides standard criteria for the classification of mental and behavioral disorders. All patients were selected on the first examination, before any medical or nutritional treatment of the developmental disorder was conducted (e.g., didn't follow GFCF diet that could have influence on the gut microbiota composition). From the interview with parents, it was reported that 33% (12/36) of patients suffered from common GI disorders like diarrhea, constipation, gas, bloating, loose stools and undigested food in stool (Supplementary Table 1). Patients hadn't taken probiotics/prebiotics, antibiotics/antimycotics or other supplements for at least 3 months prior to entering the study. All patients were recruited from the Dr. Selaković Psychiatry Clinic, Belgrade, Serbia.
Exclusion criteria for the Control group of 28 unrelated healthy children (19 males and nine females, from 2 to 11 years of age, median ± SEM, 5.96 ± 0.49 y) included the presence of any NDD or GI disorder. Inclusion in the Control group also required that volunteers had not used functional foods (such as probiotics and/or prebiotics or any supplements) or antibiotics/antimycotics for at least 3 months prior to the study.
Written informed consent was obtained from each child's parent prior to inclusion in the study. In addition, a written questionnaire was used during recruitment of volunteers for NDD and Control groups, to assess GI habit and antibiotic usage, and to look for any correlation between the characteristics of individuals and bacterial profiles. This study was carried out under the guidance of the Ethics and Research Committee of the Institute of Molecular Genetics and Genetic Engineering, University of Belgrade (O-EO-002/2015). The study conforms to the World Medical Association's Declaration of Helsinki.
Collection of Stool Specimens
Three aliquots of stool specimens (5–10 g) were collected in standard sterile sample pots with integral spoon after defecation, immediately placed at −20°C in the participant's freezer and kept there until transported to the laboratory (on ice). The samples were stored at −70°C in the laboratory until the beginning of analysis.
Denaturing Gradient Gel Electrophoresis (DGGE) Analysis
Extraction of total bacterial DNA from frozen stool samples was done using the ZR Fecal DNA MiniPrepTM kit (Zymo Research, Irvine, CA, USA) according to the instruction manual. DGGE analysis was performed as described previously (Lukic et al., 2013). Fecal samples were studied by performing a PCR-DGGE analysis with extracted genomic DNA as a template (Heilig et al., 2002). Universal primer set F-968-GC and R-1401 targeting the V6-V8 region of 16S rDNA (Nubel et al., 1996), a 16S rDNA-targeted primer set specific for lactobacilli Lab-0159f and Uni-0515-GCr (Heilig et al., 2002) and a primer set specific for Bifidobacteria Bif 164-f and Bif 662-CG-r (Satokari et al., 2001) were used to obtain an in-depth view of the microbial diversity in patients' and controls' samples (Table 1). The PCR reaction was performed in a thermal cycler (GeneAmp PCR System 2700, Applied Biosystems, Foster City, CA) programmed as follows: initial denaturation of DNA for 5 min at 95°C, 35 cycles of 30 s at 95°C, 20 s at 56°C, and 40 s at 68°C; and extension of incomplete products for 7 min at 68°C. PCR products were quantified by electrophoresis on a 1% (wt/vol) agarose gel containing ethidium bromide, and visualized by CCD camera Biometra BDR2/5/6 (Bio Doc Analyze). The amplification stage was done using KAPA Taq DNA polymerase (KAPA Biosystems, Cape Town, South Africa). PCR amplicons were set on denaturing gradient gel prepared according to Lukic et al. (2013). Optimal separation of the PCR products for the species in the fecal samples was achieved with a 40–80%, 30–60%, and 45–55% urea-formamide denaturant gradient for the Universal primer set, primer set specific for lactobacilli, and primer set specific for Bifidobacteria, respectively, increasing in the direction of electrophoresis. A 100% denaturant corresponds to 7 mol/l urea and 40% (vol/vol) formamide. Electrophoresis was performed at a constant voltage of 85 V at temperature of 60°C for 16 h. The amplicons were visualized by AgNO3 staining and developing basically as described previously (Sanguinetti et al., 1994).
Sequencing of Selected rDNA Amplicons
rDNA amplicons which had incidences that were statistically different between NDD patients and the Control groups (or if a trend existed) were selected for sequencing. The amplicons of interest were excised from the gel and macerated, and the suspension was incubated for 10 min at 98°C (Lukic et al., 2013). After incubation, the suspension was centrifuged and the supernatant (10 μl) was used for PCR amplification with universal primers F-968 and R-1401, Lactobacillus-specific primers (LB) Lab-0159f and Uni-0515r and Bifidobacterium genus-specific primers (BB) Bif164-f and Bif662-r, using the same PCR program as mentioned above (Table 1). The obtained PCR products were purified using the GeneJET PCR Purification Kit and ligated into the pJET1.2/blunt vector system according to the manufacturer's instructions (Thermo Fisher Scientific, MA, USA). Ligated constructs were transformed in Ca2+-induced competent DH5α cells (Hanahan, 1983). Transformants were selected on Luria agar plates containing 100 μg/ml ampicillin. For each excised DNA band two colonies were chosen and plasmids were isolated using the GeneJET Plasmid Miniprep Kit (Thermo Fisher Scientific, MA, USA). The insert-containing pJET1.2/blunt vectors were sequenced with pJET1.2 Forward/Reverse primers (Macrogen Europe Service, Amsterdam, the Netherlands). Sequence annotation and the database searches for sequence similarities were performed using the BLAST tool available online (https://blast.ncbi.nlm.nih.gov/Blast.cgi).
HPLC–UV Quantitative Determination of Fecal SCFAs
All procedures of SCFA isolation from fecal samples were performed as described previously (De Baere et al., 2013). Quantification of SCFAs in fecal samples was performed by high performance liquid chromatography (HPLC).
Chromatographic separation was tested on a Hypersil Gold aQ column (150 × 4.6 mm i.d.) with a particle size of 3 μm (Thermo Scientific, Voltam, MA, SAD). The HPLC–UV system consisted of an UltiMateTM 3000 LPG 3400 HPLC pump, an UltiMateTM 3000 WPS 3000 TSL autosampler, and an UltiMateTM 3000 DAD 3000 detector, all from Thermo Separations Products (Thermo Scientific, Breda, The Netherlands). The HPLC column was protected by a guard column of the same type. The column was thermostatized at 30°C. The mobile phase consisted of 20 mM of sodium dihydrogen phosphate (NaH2PO4) (Merck, Darmstadt, Germany) in HPLC water [pH adjusted to 2.2 using phosphoric acid (H3PO4)] (Merck) (A) and HPLC grade acetonitrile (C2H3N) (Sigma-Aldrich, St. Louis, Missouri, USA) (B). A gradient elution was performed as shown in Supplementary Table 2. The UV detector was set at a wavelength of 210 nm. Data processing was performed using Chromeleon version 6.8 software (Thermo Fisher Scientific, MA, USA).
For quantification of SCFAs the external calibration standards curve method was used. Calibration standards were prepared with concentrations ranging from 0.5 mM to 50 mM for acetic acid (AA), propionic acid (PA), butyric acid (BA), and succinic acid (SA) as internal standards.
Each analysis day, standards for calibration curves were freshly prepared in HPLC grade water. Samples were spiked at an individual organic acid concentration of 25 mM. All measurements were run in triplicate. Standards: acetic, propionic, butyric and succinic acid were purchased from Sigma–Aldrich, USA.
Fecal concentrations of SCFA were calculated by using the following equation: SCFA (AA, PA, BA) = (organic acid in fecal sample × 6 × 10−3)/(succinic acid in fecal sample × mass of fecal sample) × 1000 [mmol/kg].
Data are presented as mean ± SEM for continuous variables and as percentages for categorical variables. Results were considered statistically significant at P < 0.05, P value between 0.05 and 0.1 was considered a trend toward significant difference. Statistical analyses were performed using GraphPad Prism v5 software (San Diego, California, USA).
Comparisons between the two groups (NDD patient vs. Control group) were performed with a Student's t-test and Pearson's Chi-square (with Yates‘ correction for continuity) test for quantitative and categorical variables, respectively. One-way ANOVA followed by a Tukey post-hoc test were used for the analyses of differences between groups with a particular NDD diagnosis and the Control group.
Similarities between samples (NDD patient vs. Control group) were determined by calculating similarity indices based on the Dice similarity coefficient. Two identical profiles create a value of 100%, whereas two profiles without common bands result in a value of 0%. Dice similarity coefficients were used to compare DGGE band profiles between individuals within the same group and between members of the NDD patient and Control group, for every possible pair of samples. The Dice similarity coefficient was calculated using the formula: Dsc=[2j/(a+b)]*100 (%), where “a” is defined as the number of bits set to “1” in sample A, “b” as the numbers of bits set to “1” in sample B and “j” as the numbers of bits that are “1” in both samples A and B. Dsc=100% presents identical DGGE profiles, while Dsc=0% refers to completely different DGGE profiles.
Gender Bias and Increased Incidence of GI Disturbances in NDD Patients
In our study a gender bias was noticed, with the gender ratio being 3.5:1 (28 males: 8 females) in the patient group and 2.1:1 (19 males: 9 females) in the selected corresponding Control group. According to the parents' reports, 33% (12/36) of NDD patients in this study frequently had one or more GI disorders like diarrhea, constipation, gas, bloating, loose stools and undigested food in the stool. The relevant symptoms are given in Supplementary Table 1.
Dice Similarity Coefficient Analysis of Gut Bacteria in NDD Patients and Healthy Children
In order to compare the fecal microbiota of children with NDD and healthy children, DGGE analyses of rDNA amplicons obtained with three sets of primers (universal, Lactobacillus (LB) and Bifidobacteria (BB) specific primers) were performed (Figure 1).
Figure 1. The identity of the rDNA clones obtained from DGGE rDNA amplicons related to the individual children. The numbers on the y-axis correlate to the numbers of the rDNA amplicons observed. The rDNA amplicons whose incidences were statistically different between NDD patient and Control groups (or if trend existed) were excised, cloned, and sequenced and identities of those bacterial species are displayed. The numbers on the x-axis correspond to the numbers of lanes. Each lane represents sample of an individual child. Total of 64 subjects (36 NDD patients and 28 healthy controls) was analyzed. Each symbol represents a corresponding band on DGGE gel. MSDD, mixed specific developmental disorder; PDD-NOS, pervasive developmental disorder-non-specified; ERLD, expressive and receptive language disorder; CHA, childhood autism; NDD, neurodevelopmental disorder.
DGGE profiles obtained with universal primers showed large interpersonal diversity in all groups. Therefore, Dice similarity coefficients didn't show statistically significant differences concerning homogeneity of the groups or similarity between groups (Figures 2A,E). When LB-specific and BB-specific primer sets were used a high similarity of DGGE profiles was observed within the Control (CTRL) group (59.30 ± 15.21 and 71.88 ± 16.16, respectively), while the lowest similarity was scored within the ERLD group (31.12 ± 22.3 and 30.13 ± 37.50, respectively) (Figures 2B,C). When amplicons obtained with LB-specific primers were analyzed, the results revealed that among the NDD patients' groups the most homogenous groups were the MSDD (51.77 ± 18.74) and CHA groups (50.17 ± 23.56) (Figure 2B). With amplicons obtained using BB-specific primers, the most homogeneous group among NDD patients was CHA (52.00 ± 45.41) (Figure 2C). When the amplicons obtained using all three sets of primers were analyzed, a high similarity of DGGE profiles was observed within the CTRL group (64.57 ± 9.35), and the lowest similarity was observed within the ERLD group (49.90 ± 13.64). A lower level of homogeneity than CTRL, but similar to each other, was noticed among other patient groups (MSDD, 57.80 ± 12.54; PDD-NOS, 55.65 ± 12.47, and CHA 58.53 ± 12.59) (Figure 2D).
Figure 2. Box-plot diagram based on Dice similarity coefficient. Homogeneities of the groups and similarities between the groups were evaluated by Dice similarity coefficient (%). Homogeneities of the groups were evaluated by comparing DGGE profiles obtained by universal, Lactobacillus (LB)-, Bifidobacteria (BB)-specific primers and all three sets of primers together (A–D) within Control (CTRL) and particular NDD groups. Similarities between particular NDD groups and similarities of NDD groups with the CTRL group were evaluated by comparing DGGE profiles of samples from particular NDD groups, obtained by universal, LB-, BB-specific primers and all three sets of primers together, between NDD groups (E–H) and with DGGE profiles of the CTRL group (I–L), respectively. Relevant statistical values of Dice analysis are presented in Supplementary Table 3. *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001.
A Dice similarity coefficient, using LB-specific primers and BB-specific primers, revealed that the DGGE profiles of MSDD and CHA groups were similar (51.55 ± 19.89 and 52.85 ± 11.55, respectively) (Figures 2F,G). Analysis with all three sets of primers revealed that MSDD and CHA were similar (59.58 ± 12.16), as well as PDD-NOS and CHA (58.29 ± 11.55), and MSDD and PDD-NOS (56.84 ± 14.17), while the similarities of MSDD, PDD-NOS and CHA with ERLD were lower (52.74 ± 14.50, 51.90 ± 11.20, and 54.15 ± 10.58, respectively) (Figure 2H).
According to analyses of DGGE profiles obtained with universal, LB-specific, and BB-specific primers, so all three sets of primers together, the Dice similarity coefficient revealed that the ERLD group was least similar to the CTRL (59.29 ± 15.76; 40.64 ± 20.99; 46.16 ± 35.38; 50.88 ± 14.32, respectively) (Figures 2I–L). When universal primers were used, the Dice similarity coefficient showed that the PDD-NOS group was the most similar to the CTRL (64.89 ± 13.03) (Figure 2I), LB-specific primers showed that the MSDD group is the most similar to the CTRL (51.18 ± 15.97) (Figure 2J), while DGGE profiles obtained with BB-specific primers revealed that the CHA group is the most similar to the CTRL (57.35 ± 32.66) (Figure 2K). Finally, when all three sets of primers were used, Dice similarity coefficients revealed that MSDD (56.95 ± 10.59) and PDD-NOS (57.00 ± 10.75) are the most similar to the CTRL (Figure 2L). Values of dice analyses are presented in Supplementary Table 3.
Microbial Diversity Comparisons
Microbial diversity was evaluated according to a range-weighted richness index (Rr) which represents the total number of rDNA amplicons obtained with PCR-DGGE analysis, multiplied by the percentage of denaturing gradient, in order to describe the total diversity of the sample analyzed, according to the following formula:
where N represents the total number of rDNA amplicons in the pattern, and Dg the denaturing gradient comprised between the first and the last band of the pattern. A higher Rr value represents higher diversity. Comparisons between Control and NDD patient groups were performed.
It was observed that in the analyzed groups a significant difference in Rr exists between DGGE patterns obtained with universal primers (P < 0.05) (Figure 3A), LB-specific primers (P < 0.001; Figure 3B) and BB-specific primers (P < 0.01; Figure 3C).
Figure 3. Microbial diversity comparison. Student's unpaired t-test was used to compare mean values of range-weighted richness index (Rr) derived from DGGE patterns obtained with universal (Uni) primers (A) and Lactobacillus (LB)-specific primers (B) between NDD patient and Control group. Due to absence of normal distribution of range-weighted richness indexes derived from DGGE patterns obtained with Bifidobacteria (BB)-specific primers, difference between groups was assessed by non-parametric Mann–Whitney U test (C). Comparison of mean values of Rr derived from DGGE patterns obtained with universal (D) and LB-specific primers (E) between particular NDD diagnoses and Control group was performed with one-way ANOVA followed by Tukey post-hoc test, while for comparison of Rr derived from DGGE patterns obtained with BB-specific primers Kruskal-Wallis test was used (F). MSDD- mixed specific developmental disorder, PDD-NOS- pervasive developmental disorder-non-specified, ERLD, expressive and receptive language disorder; CHA, childhood autism; NDD, neurodevelopmental disorder. *p < 0.05; **p < 0.01; ***p < 0.001.
One-way ANOVA revealed statistically significant between-group differences in average Rr values for DGGE patterns obtained with LB-specific primers (P < 0.001). In addition, post hoc analyses of Rr values derived from DGGE patterns obtained with LB-specific primers revealed significant differences in Rr between the Control (131.89 ± 10.66) and each of the NDD patient groups: MSDD (54.60 ± 8.26, P < 0.001), PDD-NOS group (52.98 ± 8.92, P < 0.001), ERLD (45.93 ± 12.25, P < 0.001), and CHA (50.52 ± 15.15, P < 0.01; Figure 3E), indicating significantly lower diversity in NDD patient groups compared to the Control group. Significant differences in Rr derived from DGGE patterns obtained with universal primers between NDD patient groups and the Control group were not detected, while for Rr derived from DGGE patterns obtained with BB-specific primers a trend was noted (P = 0.08) (Figures 3D,F, respectively).
Identification of Key Bacterial Species Linked to NDD
rDNA amplicons which had incidences that were statistically different between NDD patients and the Control group (or if trends existed) were selected for sequencing. The results of sequencing and P values are presented in Table 2. The results revealed lower diversity of regular constituents of gut microbiota and higher incidence of several potentially harmful bacteria in the gut microbiota of NDD patients.
Higher Incidence of Potentially Harmful Bacteria in Fecal Microbiota of NDD Patients
Several bacterial species were found more frequently in the patients' groups compared to the control samples, specifically Desulfotomaculum guttoideum (P < 0.01), Intestinibacter bartlettii (previously known as Clostridium bartlettii) (P < 0.05) and Romboutsia ilealis (P < 0.001; Table 2). D. guttoideum was present in the highest percentage in ERLD samples (55.56%, P < 0.01), in 40% (P = 0.08) of CHA and in 30% (P = 0.08) of the PDD-NOS patients, in comparison to very low occurrence in the control samples (3.6%) (Supplementary Table 4). Intestinibacter bartlettii was present in all (100%) of the MSDD (P < 0.05), ERLD (P = 0.09) and CHA (ns) patients' samples, while it was noticed in 64.3% of the control samples. Romboutsia ilealis was present in 91.7% (P < 0.01) of MSDD, 88.9% (P = 0.01) of ERLD, 80% (P < 0.05) of the PDD-NOS and 60% (ns) of CHA patients, but in only 32% of the control samples. In addition, a high incidence of Erysipelatoclostridium ramosum (previously known as Clostridium ramosum) was noted in the ERLD (55.6%, ns), MSDD (41.7%, ns), and CHA (40%, ns) patients' groups compared to 25% in the control samples.
It was noted that I. bartlettii (100%, P < 0.05) and R. ilealis (91.7%, P < 0.01) were present in a higher number of MSDD patients compared to the Control group, while a high incidence of R. ilealis (80%, P < 0.05) and D. guttoideum (30%, P = 0.08) was observed in the PDD-NOS group. High incidences of all identified Clostridium-related species [D. guttoideum (55.6%, P < 0.01), R. ilealis (88.9%, P < 0.05), I. bartlettii (100%, P = 0.09) and E. ramosum (55.6%, ns)] were detected in the ERLD group of patients. Although in the CHA group a significant difference in comparison to the Control group wasn't achieved, probably due to the small number of samples in this group, high incidences of D. guttoideum (40%, P = 0.08), R. ilealis (60%, ns), I. bartlettii (100%, ns), and E. ramosum (40%, ns) were also noticed.
Lower Incidence of Lactic Acid Bacteria (LAB), Bifidobacteria, and Butyrate-Producing Bacteria in Fecal Microbiota of NDD Patients
Enterococcus faecalis was present in 89.3% of the control samples, while it was found in only 58.3% (P < 0.05) of NDD patients' samples, particularly in 58.3% (P = 0.07) of MSDD and 55.6% (P = 0.08) of ERLD patients' samples. Enterococcus gallinarum wasn't detected in CHA patients' samples (P = 0.08) and was present in only 11.1% of the ERLD patients (P = 0.06) and 16.7% of the MSDD samples (P = 0.07), compared to 53.6% in the control samples. Moreover, Streptococcus pasteurianus was less present in the patients' groups (P < 0.05), being completely absent in the MSDD and CHA group (P < 0.05 and ns, respectively). Among the Lactobacillus genus, Lactobacillus rhamnosus (P < 0.01) and L. sakei (P < 0.01) were found to be less frequently present in the NDD patient group. L. rhamnosus was completely absent in MSDD (P < 0.01) and present in a low percent of ERLD (11.1%, ns), CHA and PDD-NOS (20%, ns) patients, while it was present in 42.9% of the control samples.
A complete absence of Bifidobacterium sp. was noticed in 12% of NDD patients, while such absolute deficit wasn't observed in any of the control samples. Particularly, B. faecale (P < 0.05) was found in all of the control samples while it was absent in 22.2% of patients' samples. B. adolescentis was present in 28.6% of children in the Control group (P < 0.01), but was not detected in any of the patients' samples. B. stercoris (P < 0.01) and B. pseudocatenulatum (P = 0.08) were also less frequently present in the NDD patient group. B. pseudocatenulatum was absent in all of the CHA patients and was present in only 10% of PDD-NOS patients, while among the Control group it was present in 35.7% of samples. Although not significant, a lower incidence of B. longum subsp. infantis was observed in the patients' groups. In contrast, B. animalis subsp. lactis was observed to be present in all patients' groups except CHA and it was present in the highest percentage in PDD-NOS samples (40%, P < 0.01), while not detected in any of the control samples (Supplementary Table 4).
Similarly, the butyrate producing bacteria Faecalibacterium prausnitzii (P < 0.01), Butyricicoccus pullicaecorum (P < 0.05) and Eubacterium rectale (P = 0.07) were less frequent in the patients' groups (Table 2). Faecalibacterium prausnitzii was present in only 20% (ns) of CHA and 22.2% (P < 0.05) of ERLD patients' samples while it was present in 67.9% of the control samples. Butyricicoccus pullicaecorum was completely absent in CHA (ns) and MSDD patients' samples (P < 0.05), while it was present in 42.9% of the control samples. Eubacterium rectale was not detected in any of the CHA patient's samples (ns) and was detected in only 20% (ns) of PDD-NOS samples, while it was present in 46.4% of the control samples.
In addition, Ruminococcus champanellensis (P < 0.01) and Dialister invisus (P = 0.001) were also less frequent in the patient group. R. champanellensis was absent in CHA (P = 0.08) and detected in only 11.1% (P = 0.06) of ERLD samples, while it was present in 53.6% of the control samples. D. invisus was present in 41.7% (P < 0.05) of the MSDD samples, 40% of the PDD-NOS and CHA samples (P < 0.05 and ns, respectively) and 33.3% (P < 0.05) of ERLD samples, while it was present in 82.1% of the control samples (Supplementary Table 4).
Although significant differences between CHA and the Control group weren't detected, deficiency of all of the identified constituents of normal gut microbiota was observed in this group, the most prominent being complete absence of: E. rectale, S. pasteurianus, B. pullicaecorum, R. champanellensis, E. gallinarum, B. pseudocatenulatum and B. adolescentis.
Furthermore, the complete absence of S. pasteurianus, B. pullicaecorum, L. rhamnosus and B. adolescentis was noted in MSDD patients.
Report of Prominent Case
The most striking example of microbiota diversity disturbance associated with NDD in our study was the case of a 2-year-old boy, diagnosed as MSDD (Figure 1: patient number 7) whose DGGE pattern obtained with universal primers showed extremely decreased diversity compared to samples of other patients and controls. The two strains that were most present in this patient's sample were the Clostridium-like strains I. bartlettii and E. ramosum. PCR-DGGE analyses have shown a complete absence of Bifidobacterium sp., L. rhamnosus, E. rectale, B. pullicaecorum, and D. invisus in this individual's sample. The patient's parents reported periodic constipation, while developmental delay, hypotony, and anxiety were clinically observed.
Analysis of SCFAs in Fecal Samples of NDD Patients and Healthy Children
The amount of SCFAs in fecal samples, including acetate, propionate and butyrate, were analyzed. The results revealed that the SCFA content largely varied within groups. Interestingly, the values of total SCFAs, as well as levels of acetic and butyric acid, were the lowest in the CHA patient group, although the differences between patient and control groups were not significant, presumably due to the small number of CHA patients (Figure 4).
Figure 4. Amounts of SCFAs in fecal samples of NDD patients and healthy children. The amounts of acetic acid (A), propionic acid (B), butyric acid (C), and total SCFAs (D) in NDD, Control group and particular NDD diagnoses groups were presented in mmol/kg of feces. Quantification of SCFAs in fecal samples was performed by HPLC analysis. Results are presented as Mean ± SEM. MSDD, mixed specific developmental disorder; PDD-NOS, pervasive developmental disorder-non-specified; ERLD, expressive and receptive language disorder; CHA, childhood autism; NDD, neurodevelopmental disorder.
Correlation analysis revealed that a significant positive correlation exists between number of rDNA amplicons obtained with universal primers and the level of propionic acid (r = 0.40, P < 0.05), as well as a trend for level of total SCFAs (r = 0.38, P = 0.06) and butyric acid (r = 0.37, P = 0.06) in the Control group. This correlation is lost in the NDD patient groups, indicating that the microbiota of NDD patients differs from microbiota of healthy children in the presence and/or number of strong SCFA-producing bacteria (Figure 5).
Figure 5. Correlation between average number of rDNA amplicons and amounts of SCFAs in Control and NDD patient group. Average number of rDNA amplicons obtained with universal primers was correlated with amount of propionic acid (A), butyric acid (B), and total SCFAs (C) in Control group and in NDD group (D–F). Each dot corresponds to one subject. The line denotes linear regression of individual data point; r stands for Pearson's correlation coefficient. P < 0.05 was considered statistically significant; NDD, neurodevelopmental disorder.
Recently, differences in composition of gut microbiota from children with autism and healthy children were intensively studied (Hughes et al., 2018). However, studies evaluating diversity and composition of gut microbiota in other NDDs are lacking. Hence, according to our knowledge, this is the first study to evaluate the specificities of gut microbiota in a broad spectrum of NDD patients. The impact of any exclusion diet, such as a GFCF diet or of any medical treatment was excluded here because all patients were selected on the first examination.
A gender bias toward boys was noticed in the patients' groups (7:2) in accordance with observed diagnosis sex bias in NDDs such as autism (4:1) (Scott et al., 2002) and intellectual disability/developmental delay 2:1 (Gillberg et al., 2006; Ropers, 2008; Polyak et al., 2015). As reported by parents, 1/3 (33%) of NDD patients in this study frequently had one or more GI disorders including diarrhea, constipation, gas, bloating, loose stools and undigested food in stools, in accordance with clinical medical records describing the occurrence of GI symptoms in 24% of ASD patients (Molloy and Manning-Courtney, 2003) and in 22% of PDD patients (Nikolov et al., 2009).
According to the Dice similarity coefficient acquired when DGGE profiles obtained with LB- and BB- specific primers were compared, as well as when DGGE profiles obtained with all three sets of primers were compared together, the microbiota of the Control and CHA groups were the most homogenous, while microbiota of the ERLD group was the most heterogeneous one, indicating large interpersonal diversity in this patient group. DGGE profiles of the MSDD group also showed high homogeneity when LB-specific primers were used.
When DGGE profiles obtained with all three sets of primers were compared between groups, a high similarity was observed between CHA and PDD-NOS (also known as “atypical autism”) (Volkmar and Reichow, 2013) groups, which are also the most similar groups of NDD according to clinical/diagnostic assessments. High similarity was noticed between MSDD and PDD-NOS, and MSDD and CHA groups. The similarities of all these three groups (CHA, PDD-NOS, and MSDD) with the ERLD group was low. This is in concordance to the clinical similarity of MDD, CHA and PDD-NOS, which are more complex and severe disorders than ERLD, having similar levels of functionality and more complex and severe autism-like symptoms compared to ERLD (WHO, 2015). Similarity of microbiota composition in patient groups which have similar neurological and behavioral symptoms might lead to the assumption that intestinal microbiota have a significant role in the pathophysiology of NDDs.
The microbiota of ERLD patients showed the smallest degree of similarity with the microbiota of the Control group. This result indicates that modification of intestinal microbiota, in children with less serious NDD, such as ERLD, might be important. This assumption is in agreement with the observed efficacy of rehabilitation and decline of autism-like symptoms in this group of patients after medical treatment of dysbiosis in clinical practice.
Analysis of DGGE profiles obtained with universal primers revealed that the microbiota of PDD-NOS patients is the most similar to the microbiota of the Control group, which might indicate that disturbances in composition of these bacteria might have lower significance in this group of patients. However, R. ilealis (80%, P < 0.05) was significantly more present in this group of patients compared to the Control group. Further, DGGE profiles obtained with LB-specific primers showed that the MSDD group is the most similar to the Control group, although complete absence of S. pasteurianus, B. pullicaecorum, L. rhamnosus and lower incidence of D. invisus was noted in this group. Finally, even though analysis of the DGGE profile obtained with BB-specific primers showed that the CHA group is the most similar to the Control group, the lowest number of Bifidobacteria and complete absence of B. pseudocatenulatum and B. adolescentis were noted in this group. We assume that this similarity is obtained according to the absence of B. animalis, which is also absent in the Control group, while present in all other NDD groups.
The results of this study reveal that diversity of common commensal bacteria in the patients' groups was significantly lower compared to the Control group, while the incidence of several potentially harmful bacteria was higher. Microbial species such as Desulfotomaculum guttoideum, Intestinibacter bartlettii and Romboutsia ilealis, all closely related to Clostridium clusters (Stackebrandt et al., 1997; Gerritsen et al., 2014), were found to be more frequently present in the NDD patients group. D. guttoideum (reclassification to the Clostridium genus proposed two decades ago based on 16S rDNA analysis) (Stackebrandt et al., 1997; Castro et al., 2000) is a sulfate-reducing bacterium similar to Desulfovibrio, which has been found in increased incidence in children with autism (Finegold et al., 2010, 2012; Tomova et al., 2015). It was suggested that sulfate-reducing organisms could account for much of the abnormality in sulfur metabolism seen in autistic subjects, such as low blood levels of sulfur and high urinary excretion, reduced methylation capacity, decrease in trans-sulfuration and chronic oxidative stress (Finegold et al., 2012). According to increased Desulfovibrio and Clostridia species in the microbiota of children with autism it was speculated that these bacteria could be involved in development of autism and that these bacteria could contribute in some way to autism pathogenesis (Tomova et al., 2015).
Several research groups that have analyzed microbiota in children with autism have detected overgrowth of Clostridium spp. (Finegold et al., 2002; Song Y. et al., 2004; Parracho et al., 2005; De Angelis et al., 2013). It was hypothesized that regressive autism could be a consequence of gut microbiota dysbiosis caused by antimicrobial therapy. Overgrowth of antimicrobial resistant microorganisms, such as Clostridium tetani, leads to elevated production of the neurotoxin p-cresol or toxic metabolites (phenols and indole derivatives) which cause anxiety-like behavior (Bolte, 1998; Finegold et al., 2012; Hsiao et al., 2013). In line with that, it was reported that cognitive and social functions in ASD patients were improved to some extent after vancomycin treatment (Sandler et al., 2000).
In addition, it was revealed that I. bartlettii, previously known as Clostridium bartlettii isolated from the feces of two autistic male children (Song Y. L. et al., 2004), is a strong producer of trans-3-indolylacrylic acid (IAA) (Russell et al., 2013), which can be transformed by glycine conjugation (Clayton, 2012) to Indolyl-3-acryloylglycine (IAG), a putative ASD diagnostic urinary marker (Shattock and Whiteley, 2002; Bull et al., 2003). Moreover, Erysipelatoclostridium ramosum (previously known as Clostridium ramosum), which was noted more frequently in the ERLD, MSDD, and CHA patients compared to the Control group in this study, was formerly isolated from the stool specimens of children with autism (Finegold et al., 2002).
Even though the interpersonal diversity of both Control and NDD patient groups was large, a significant decrease in diversity of LAB was detected. The identified strains that were present in control samples with greater incidence compared to patients' samples were common commensal bacteria and butyrate producing bacteria. En. faecalis, En. gallinarum, and S. pasteurianus, which are part of the commensal microbiota (Layton et al., 2010; Alex et al., 2013; Rajilic-Stojanovic and De Vos, 2014), were detected in low numbers or even absent in some NDD patients' groups. This is in accordance with literature data reporting decreased levels of LAB (Enterococcus sp., Lactobacillus sp., Streptococcus sp., and Lactococcus sp.) in the fecal samples of children with autism and PDD-NOS (De Angelis et al., 2013). Among Lactobacillus sp., Lactobacillus rhamnosus was found to be less frequently present in the NDD patient group. It was shown that ingestion of L. rhamnosus JB-1 has an impact on the regulation of emotional behavior associated with anxiety and depression in normal, healthy mice, through modulation of central GABA receptor expression (Bravo et al., 2011).
Similarly, decreased frequency of Dialister invisus was noticed in the NDD patient group. Although D. invisus species are typically isolated from the oral cavity, they have also been detected in samples of the normal gastrointestinal microbiota (Rajilic-Stojanovic et al., 2007). Decreased populations of D. invisus have been reported in autistic children and its importance in contributing to protective flora was suggested (Finegold et al., 2010).
Moreover, a significant decrease of Bifidobacteria was noticed, including the complete absence of any detectable Bifidobacteria in 12% of NDD patients, while such an absolute deficit wasn't observed in any of the control samples. This is in accordance with previous reports where lower levels of Bifidobacterium spp. were found in samples from children with autism (Adams et al., 2011; Wang et al., 2011; De Angelis et al., 2013). It is important to emphasize that Bifidobacteria are among the first microbial colonizers of the gut immediately after birth (Turroni et al., 2012), a dominant bacterial genus in the infant gut microbiota (Harmsen et al., 2000; Matamoros et al., 2013), and significant member of commensal gut bacteria in adults (Ringel-Kulka et al., 2013; Chaplin et al., 2015). Their beneficial role for the host is assigned to digestion of indigestible polysaccharides, production of B-group vitamins, competitive pathogen exclusion, production of antimicrobial agents, and suppression of inflammatory responses (Deguchi et al., 1985; Sonnenburg et al., 2006; Lee and O'sullivan, 2010; Miyauchi et al., 2013; Mulle et al., 2013; Turroni et al., 2014). Knowing that children with autism have digestion problems, deficiency in vitamin B12, and immune dysregulation and gut inflammation (Rossignol and Frye, 2012; Wasilewska and Klukowski, 2015; Zhang et al., 2016) it could be hypothesized that these problems could be associated with deficiency of intestinal Bifidobacteria.
Additionally, it was suggested that particular Bifidobacterium strains have the potential to treat some mental disorders, like depression, and anxiety. Treatment with B. infantis was shown to be beneficial in the treatment of depression, through elevation of the serotonergic precursor-tryptophan (Desbonnet et al., 2008), and that it can improve hypothalamic-pituitary axis mediated stress response (Sudo et al., 2004), while treatment with B. longum was shown to have an anxiolytic effect (Bercik et al., 2011).
It could be assumed that low Bifidobacteria diversity could be associated with antibiotic treatment in the first years of life, in line with results revealing the high incidence of opportunistic pathogens from Clostridium cluster. This assumption is in accordance with literature data describing the greater diversity of Bifidobacterium spp. in healthy children compared to antibiotic-treated infants (Hussey et al., 2011) and significantly higher oral antibiotic use in children with autism vs. typical children (Konstantareas and Homatidis, 1987; Niehus and Lord, 2006). Interestingly, our results revealed frequent occurrence of Bifidobacterium animalis subsp. lactis in all NDD groups, except in CHA, while it was not detected in any of the control samples, similar to the previously published data (Finegold et al., 2010). The higher frequency of B. animalis might be a consequence of frequent probiotic use in the patient group due to GI problems, since various strains of this species are commonly used as probiotics (Merenstein et al., 2015).
Three butyrate producing bacteria, Faecalibacterium prausnitzii, Butyricicoccus pullicaecorum, and Eubacterium rectale, were less frequent in the patients' groups. Butyric acid is the preferred source of energy for colonocytes and affects their proliferation, differentiation and apoptosis, has anti-inflammatory effects, reinforces the colonic defense barrier by increasing production of mucins and antimicrobial peptides and decreases intestinal epithelial permeability (Hamer et al., 2008; Van Immerseel et al., 2010). Faecalibacterium prausnitzii is one of the most abundant bacteria in the human gut ecosystem (Khan et al., 2012). Faecalibacterium prausnitzii and Butyricicoccus pullicaecorum are important suppliers of butyrate to the colonic epithelium and depletion of these bacteria was associated with irritable bowel syndrome (IBS) (Steppe et al., 2014), while it was shown that B. pullicaecorum strengthens epithelial barrier function (Eeckhaut et al., 2013). Another non-amylolytic member of the Ruminococcacea related to F. prausnitzii, Ruminococcus champanellensis (Ze et al., 2015), was less frequently found in the patients' groups, particularly in CHA and ERLD samples. R. champanellensis is the only human colonic bacterium reported to be capable of degrading crystalline cellulose and fermentation of dietary non-digestible carbohydrates. It has important consequences for health via modulation of microbiota composition by supplying energy and supporting microbial growth (Chassard et al., 2010; Ben David et al., 2015). A study of fecal microbiota of children with autism and PDD-NOS revealed that Faecalibacterium, Ruminococcus, and Eubacteriaceae were found to be present at the lowest level in fecal samples of children with autism (De Angelis et al., 2013), which is in accordance with our results.
Several molecular mechanisms that could explain the influence of microbiota on the brain and that could be involved in the pathogenesis of autism and other NDDs have been revealed so far. SCFAs, potent bioactive molecules, produced by bacterial fermentation of dietary carbohydrates within the gut, of which acetate, propionate, and butyrate are the most abundant (95%), can induce widespread effects on gut, brain, and behavior (Den Besten et al., 2013). In this study it was observed that a positive correlation between number of rDNA amplicons obtained with universal primers and level of propionic acid, as well as a trend for level of total SCFAs and butyric acid in the Control group, was lost in the patients' groups, leading to the conclusion that lower numbers of strong SCFA-producing bacteria are present in the microbiota of NDD patients.
It was found that physiological concentrations of SCFA promote epithelial barrier function in the large intestine (Suzuki et al., 2008). Acetate, propionate and butyrate alone or in combination stimulated the formation of tight junctions (TJ) (Feng et al., 2018) by increasing expression of TJ proteins (Hamer et al., 2008) through their histone deacetylase (HDAC) inhibitor activity (Ohata et al., 2005). Therefore, lower levels of SCFAs might disturb epithelial barrier function in the gut. Since decreased intestinal epithelial barrier function was found in some children with autism (D'eufemia et al., 1996; Boukthir et al., 2010; De Magistris et al., 2010), absence of strong SCFA-producing bacteria and lower levels of SCFAs in the gut could be the cause of this defect.
Along with the intestinal epithelial barrier dysfunction, increased blood-brain barrier permeability was detected in ASD patients (Fiorentino et al., 2016). It was discovered that normal gut microbiota have beneficial effects on stability of the blood-brain barrier (BBB) (Braniste et al., 2014). Increased BBB permeability in germ-free mice was associated with decreased levels of TJ proteins in several brain regions. These changes were restored to normal levels by colonization of germ-free animals with intestinal microbiota from pathogen-free mice (conventionalization), colonization with a single butyrate producing bacteria or oral gavage with sodium butyrate. SCFAs seem to play an important role in maintaining BBB integrity, which has a central role in brain development and the preservation of CNS homeostasis (Silva et al., 2020).
Besides influencing the intestinal epithelial barrier and BBB, it has been shown that microbial products, propionic and butyric acid, have strong epigenetic potential and can induce alterations in expression of genes involved in neurotransmitter systems, neuronal cell adhesion molecules, inflammation, oxidative stress, lipid metabolism and mitochondrial function, all of which have been associated with ASD development (Nankova et al., 2014). In addition, treatment with sodium-butyrate induced an increase in neurotrophin levels (BDNF, NGF, and GDNF) in the hippocampus and frontal cortex, which was associated with improvement in recognition memory (Valvassori et al., 2014), and reduction in hyperlocomotion (Varela et al., 2015) in rats. Sodium-butyrate treatments have also enhanced neurite outgrowth (Suzuki-Mizushima et al., 2002) and shown neuroprotective capacity (Wu et al., 2008). These findings imply that butyrate up-regulates expression of pro-survival, pro-regenerative and pro-plasticity genes by increasing histone acetylation around promoters of these genes (Bourassa et al., 2016).
Although in this study significant differences in the levels of SCFAs weren't detected, it was observed that the values for total SCFAs, as well as butyric acid, were lowest in the CHA group. Literature data revealed that total SCFA levels were significantly lower in autism disorder (AD) and PDD-NOS patients (Adams et al., 2011; De Angelis et al., 2013). Our assumption is that a significant difference in SCFA levels between the CHA subgroup and Control group in our study was not detected due to the small number of CHA patients (ratio CHA:NDD patients is 1: 7.2). This is to some extent in accordance with an estimation by the National Survey of Children's Health in the US for 2007, according to which 1.1% of children ages 3–17 years were reported to have ASDs, while 15% were affected by NDD (ratio 1:15). This estimation included other developmental delays, such as intellectual disability, cerebral palsy, seizures, stuttering or stammering, moderate to profound hearing loss, and blindness (Kogan et al., 2009; Boyle et al., 2011) that weren't included in this study.
As stated by other authors, a positive correlation was found between the level of total SCFAs and the levels of Faecalibacterium, Ruminococcus, and Bifidobacterium (De Angelis et al., 2013), which is in accordance with our observation that F. prausnitzii and R. champanellensis are present with lower frequency in the patient group, as were several strains from the Bifidobacterium genus. Although bacteria belonging to Ruminococcus and Bifidobacterium genera do not produce butyrate, they can promote the production of SCFAs thanks to their ability to process non-digestible carbohydrates and can stimulate butyrate production by cross-feeding of colonic butyrate-producing bacteria (Belenguer et al., 2006).
According to current knowledge, it can be assumed that SCFAs are important molecular mediators in the gut-brain axis and whose effects might be involved in development of autism-like symptoms (Rogers et al., 2016). Therefore, lower numbers of strong SCFA-producing bacteria in the gut microbiota of NDD patients, especially a lower incidence of butyrate producing bacteria, might have a significant impact on the etiology of these disorders.
Recently, it has been recognized that microbiota have a significant role in neurodevelopment due to the intense and dynamic bidirectional communication in the gut-brain axis, and that early-life perturbations of initial microbiota colonization can lead to mental disorders later in life (Diaz Heijtz et al., 2011; Borre et al., 2014; De Theije et al., 2014). Development of the microbiota occurs in parallel with neurodevelopment and they have similar critical developmental windows sensitive to disruption (O'mahony et al., 2017). It has been suggested that colonization by gut microbiota impacts mammalian brain development and subsequent adult behavior in part by altered expression of several genes (Diaz Heijtz et al., 2011). It has also been emphasized that a critical period might exist after which reconstitution of microbiota does not normalize the behavioral phenotype (Sudo et al., 2004; Diaz Heijtz et al., 2011; Neufeld K. A. et al., 2011; Desbonnet et al., 2014) or neurochemical disturbances that have occurred in early phases of brain development (Clarke et al., 2013).
Our study reveals that the intestinal microbiota from NDD patients differs from the microbiota of healthy children. Lower bacterial diversity of beneficial commensal bacteria and increased incidence of resident microbes with pathogenic potential were detected. The alterations reported in the present study are partly consistent with other studies performed on children with autism and support the hypothesis that microbiota dysbiosis could be associated with GI disorders frequently observed in NDD patients and appearance of autism-like behavior (Finegold et al., 2002, 2010; Song Y. et al., 2004; Parracho et al., 2005; Adams et al., 2011; De Angelis et al., 2013; Tomova et al., 2015). The characterization of gut microbiota in our study, as in aforementioned studies, was performed according to analyses of fecal samples. Some authors have pointed out that there are differences in composition of fecal and mucosa associated microbiota in healthy individuals (Carstens et al., 2018), such that mucosal-associated microbiota better discriminate patients from controls (Altomare et al., 2019) and therefore deserve more attention (Ouwehand et al., 2004). Although this must be kept in mind, fecal samples are often used as representative of gut microbiota whenever colonic biopsies are not available.
Our results provide a rationale for further studies of gut microbiota dysbiosis in NDD patients and elucidation of possible pathogenic modes of action of specific bacteria. According to these results and literature data, we suggest that early life microbiome screening could be a potentially useful tool for prevention and treatment of GI disturbances and accompanied behavioral problems frequently observed in autism and other patients from a broad spectrum of NDD. In addition, correction of microbiota composition and supplementation with several health promoting strains, found to be deficient in these patients, could be a safe adjuvant therapy in the treatment of NDD accompanied by GI disturbances.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.
The study involving human participants was reviewed and approved by Ethics and Research Committee of Institute of Molecular Genetics and Genetic Engineering, University of Belgrade (O-EO-002/2015). Written informed consent to participate in this study was provided by the participants' parent.
KB contributed to taking interviews with parents, collecting anamnesis' data and the samples from patients, carrying out the interpretation of data, and preparing the first draft of the manuscript. ÐI contributed to making conception and design of the study, collecting samples from healthy children, carrying out the experiments, analysis and interpretation of data, preparing the first draft of the manuscript, performing checking, and editing of all revisions of the manuscript. SS contributed to carrying out the experiments, preparing the first draft of the manuscript. DV contributed to analysis and interpretation of the data, performing checking and editing of all revisions of the manuscript. MTom contributed to editing final revision of the manuscript. NG contributed to making a conception and design of the study, doing analysis and interpretation of data, performing checking, and editing of all revisions of the manuscript. MTol contributed to making conception, design and supervision of the study, carrying out the experiments, analysis and interpretation of data, preparing the first draft of the manuscript, performing checking, and editing of all revisions of the manuscript.
This work was funded by the Ministry of Education, Science and Technological Development of the Republic of Serbia, grants OI173019 and III 41030, and contract on the implementation and financing of research work in 2020, registration numbers: 451-03-68/2020-14/200042 and 451-03-68/2020-14/200007.
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.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2020.00223/full#supplementary-material
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
Altomare, A., Putignani, L., Del Chierico, F., Cocca, S., Angeletti, S., Ciccozzi, M., et al. (2019). Gut mucosal-associated microbiota better discloses inflammatory bowel disease differential patterns than faecal microbiota. Dig. Liver Dis. 51, 648–656. doi: 10.1016/j.dld.2018.11.021
Belenguer, A., Duncan, S. H., Calder, A. G., Holtrop, G., Louis, P., Lobley, G. E., et al. (2006). Two routes of metabolic cross-feeding between bifidobacterium adolescentis and butyrate-producing anaerobes from the human gut. Appl. Environ. Microbiol. 72, 3593–3599. doi: 10.1128/AEM.72.5.3593-3599.2006
Ben David, Y., Dassa, B., Borovok, I., Lamed, R., Koropatkin, N. M., Martens, E. C., et al. (2015). Ruminococcal cellulosome systems from rumen to human. Environ. Microbiol. 17, 3407–3426. doi: 10.1111/1462-2920.12868
Bercik, P., Park, A. J., Sinclair, D., Khoshdel, A., Lu, J., Huang, X., et al. (2011). The anxiolytic effect of Bifidobacterium longum NCC3001 involves vagal pathways for gut-brain communication. Neurogastroenterol. Motil. 23, 1132–1139. doi: 10.1111/j.1365-2982.2011.01796.x
Borre, Y. E., O'keeffe, G. W., Clarke, G., Stanton, C., Dinan, T. G., and Cryan, J. F. (2014). Microbiota and neurodevelopmental windows: implications for brain disorders. Trends Mol. Med. 20, 509–518. doi: 10.1016/j.molmed.2014.05.002
Bourassa, M. W., Alim, I., Bultman, S. J., and Ratan, R. R. (2016). Butyrate, neuroepigenetics and the gut microbiome: can a high fiber diet improve brain health? Neurosci. Lett. 625, 56–63. doi: 10.1016/j.neulet.2016.02.009
Boyle, C. A., Boulet, S., Schieve, L. A., Cohen, R. A., Blumberg, S. J., Yeargin-Allsopp, M., et al. (2011). Trends in the prevalence of developmental disabilities in US children, 1997-2008. Pediatrics 127, 1034–1042. doi: 10.1542/peds.2010-2989
Braniste, V., Al-Asmakh, M., Kowal, C., Anuar, F., Abbaspour, A., Toth, M., et al. (2014). The gut microbiota influences blood-brain barrier permeability in mice. Sci. Transl Med. 6:263ra158. doi: 10.1126/scitranslmed.3009759
Bravo, J. A., Forsythe, P., Chew, M. V., Escaravage, E., Savignac, H. M., Dinan, T. G., et al. (2011). Ingestion of Lactobacillus strain regulates emotional behavior and central GABA receptor expression in a mouse via the vagus nerve. Proc. Natl. Acad. Sci. U.S.A. 108, 16050–16055. doi: 10.1073/pnas.1102999108
Bull, G., Shattock, P., Whiteley, P., Anderson, R., Groundwater, P. W., Lough, J. W., et al. (2003). Indolyl-3-acryloylglycine (IAG) is a putative diagnostic urinary marker for autism spectrum disorders. Med. Sci. Monit. 9, CR422–425.
Carstens, A., Roos, A., Andreasson, A., Magnuson, A., Agreus, L., Halfvarson, J., et al. (2018). Differential clustering of fecal and mucosa-associated microbiota in 'healthy' individuals. J. Dig. Dis. 19, 745–752. doi: 10.1111/1751-2980.12688
Chaidez, V., Hansen, R. L., and Hertz-Picciotto, I. (2014). Gastrointestinal problems in children with autism, developmental delays or typical development. J. Autism. Dev. Disord. 44, 1117–1127. doi: 10.1007/s10803-013-1973-x
Chaplin, A. V., Efimov, B. A., Smeianov, V. V., Kafarskaia, L. I., Pikina, A. P., and Shkoporov, A. N. (2015). Intraspecies genomic diversity and long-term persistence of bifidobacterium longum. PLoS ONE 10:e0135658. doi: 10.1371/journal.pone.0135658
Chassard, C., Delmas, E., Robert, C., and Bernalier-Donadille, A. (2010). The cellulose-degrading microbial community of the human gut varies according to the presence or absence of methanogens. FEMS Microbiol. Ecol. 74, 205–213. doi: 10.1111/j.1574-6941.2010.00941.x
Clarke, G., Grenham, S., Scully, P., Fitzgerald, P., Moloney, R. D., Shanahan, F., et al. (2013). The microbiome-gut-brain axis during early life regulates the hippocampal serotonergic system in a sex-dependent manner. Mol. Psychiatry 18, 666–673. doi: 10.1038/mp.2012.77
Clayton, T. A. (2012). Metabolic differences underlying two distinct rat urinary phenotypes, a suggested role for gut microbial metabolism of phenylalanine and a possible connection to autism. FEBS Lett. 586, 956–961. doi: 10.1016/j.febslet.2012.01.049
De Angelis, M., Francavilla, R., Piccolo, M., De Giacomo, A., and Gobbetti, M. (2015). Autism spectrum disorders and intestinal microbiota. Gut. Microbes 6, 207–213. doi: 10.1080/19490976.2015.1035855
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
De Baere, S., Eeckhaut, V., Steppe, M., De Maesschalck, C., De Backer, P., Van Immerseel, F., et al. (2013). Development of a HPLC-UV method for the quantitative determination of four short-chain fatty acids and lactic acid produced by intestinal bacteria during in vitro fermentation. J. Pharm. Biomed. Anal. 80, 107–115. doi: 10.1016/j.jpba.2013.02.032
De Magistris, L., Familiari, V., Pascotto, A., Sapone, A., Frolli, A., Iardino, P., et al. (2010). Alterations of the intestinal barrier in patients with autism spectrum disorders and in their first-degree relatives. J. Pediatr. Gastroenterol. Nutr. 51, 418–424. doi: 10.1097/MPG.0b013e3181dcc4a5
De Theije, C. G., Wopereis, H., Ramadan, M., Van Eijndthoven, T., Lambert, J., Knol, J., et al. (2014). Altered gut microbiota and activity in a murine model of autism spectrum disorders. Brain Behav. Immun. 37, 197–206. doi: 10.1016/j.bbi.2013.12.005
Den Besten, G., Van Eunen, K., Groen, A. K., Venema, K., Reijngoud, D. J., and Bakker, B. M. (2013). The role of short-chain fatty acids in the interplay between diet, gut microbiota, and host energy metabolism. J. Lipid. Res. 54, 2325–2340. doi: 10.1194/jlr.R036012
Desbonnet, L., Garrett, L., Clarke, G., Bienenstock, J., and Dinan, T. G. (2008). The probiotic bifidobacteria infantis: an assessment of potential antidepressant properties in the rat. J. Psychiatr. Res. 43, 164–174. doi: 10.1016/j.jpsychires.2008.03.009
Diaz Heijtz, R., Wang, S., Anuar, F., Qian, Y., Bjorkholm, B., Samuelsson, A., et al. (2011). Normal gut microbiota modulates brain development and behavior. Proc. Natl. Acad. Sci. U.S.A. 108, 3047–3052. doi: 10.1073/pnas.1010529108
Eeckhaut, V., Machiels, K., Perrier, C., Romero, C., Maes, S., Flahou, B., et al. (2013). Butyricicoccus pullicaecorum in inflammatory bowel disease. Gut 62, 1745–1752. doi: 10.1136/gutjnl-2012-303611
Feng, Y., Wang, Y., Wang, P., Huang, Y., and Wang, F. (2018). Short-chain fatty acids manifest stimulative and protective effects on intestinal barrier function through the inhibition of NLRP3 inflammasome and autophagy. Cell. Physiol. Biochem. 49, 190–205. doi: 10.1159/000492853
Finegold, S. M., Dowd, S. E., Gontcharova, V., Liu, C., Henley, K. E., Wolcott, R. D., et al. (2010). Pyrosequencing study of fecal microflora of autistic and control children. Anaerobe 16, 444–453. doi: 10.1016/j.anaerobe.2010.06.008
Finegold, S. M., Molitoris, D., Song, Y., Liu, C., Vaisanen, M. L., Bolte, E., et al. (2002). Gastrointestinal microflora studies in late-onset autism. Clin. Infect. Dis. 35(Suppl. 1), S6–S16. doi: 10.1086/341914
Fiorentino, M., Sapone, A., Senger, S., Camhi, S. S., Kadzielski, S. M., Buie, T. M., et al. (2016). Blood-brain barrier and intestinal epithelial barrier alterations in autism spectrum disorders. Mol. Autism. 7:49. doi: 10.1186/s13229-016-0110-z
Frye, R. E., Slattery, J., Macfabe, D. F., Allen-Vercoe, E., Parker, W., Rodakis, J., et al. (2015). Approaches to studying and manipulating the enteric microbiome to improve autism symptoms. Microb. Ecol. Health Dis. 26:26878. doi: 10.3402/mehd.v26.26878
Gerritsen, J., Fuentes, S., Grievink, W., Van Niftrik, L., Tindall, B. J., Timmerman, H. M., et al. (2014). Characterization of Romboutsia ilealis gen. nov., sp. nov., isolated from the gastro-intestinal tract of a rat, and proposal for the reclassification of five closely related members of the genus Clostridium into the genera Romboutsia gen. nov., Intestinibacter gen. nov., Terrisporobacter gen. nov. and Asaccharospora gen. nov. Int. J. Syst. Evol. Microbiol. 64(Pt 5), 1600–1616. doi: 10.1099/ijs.0.059543-0
Gillberg, C., Cederlund, M., Lamberg, K., and Zeijlon, L. (2006). Brief report: “the autism epidemic”. The registered prevalence of autism in a Swedish urban area. J. Autism. Dev. Disord. 36, 429–435. doi: 10.1007/s10803-006-0081-6
Hamer, H. M., Jonkers, D., Venema, K., Vanhoutvin, S., Troost, F. J., and Brummer, R. J. (2008). Review article: the role of butyrate on colonic function. Aliment. Pharmacol. Ther. 27, 104–119. doi: 10.1111/j.1365-2036.2007.03562.x
Harmsen, H. J., Wildeboer-Veloo, A. C., Raangs, G. C., Wagendorp, A. A., Klijn, N., Bindels, J. G., et al. (2000). Analysis of intestinal flora development in breast-fed and formula-fed infants by using molecular identification and detection methods. J. Pediatr. Gastroenterol. Nutr. 30, 61–67. doi: 10.1097/00005176-200001000-00019
Heilig, H. G., Zoetendal, E. G., Vaughan, E. E., Marteau, P., Akkermans, A. D., and De Vos, W. M. (2002). Molecular diversity of Lactobacillus spp. and other lactic acid bacteria in the human intestine as determined by specific amplification of 16S ribosomal DNA. Appl. Environ. Microbiol. 68, 114–123. doi: 10.1128/aem.68.1.114-123.2002
Hsiao, E. Y., Mcbride, S. W., Hsien, S., Sharon, G., Hyde, E. R., Mccue, T., et al. (2013). Microbiota modulate behavioral and physiological abnormalities associated with neurodevelopmental disorders. Cell 155, 1451–1463. doi: 10.1016/j.cell.2013.11.024
Hussey, S., Wall, R., Gruffman, E., O'sullivan, L., Ryan, C. A., Murphy, B., et al. (2011). Parenteral antibiotics reduce bifidobacteria colonization and diversity in neonates. Int. J. Microbiol. 2011:1–6. doi: 10.1155/2011/130574
Khan, M. T., Duncan, S. H., Stams, A. J., Van Dijl, J. M., Flint, H. J., and Harmsen, H. J. (2012). The gut anaerobe Faecalibacterium prausnitzii uses an extracellular electron shuttle to grow at oxic-anoxic interphases. ISME J. 6, 1578–1585. doi: 10.1038/ismej.2012.5
Kogan, M. D., Blumberg, S. J., Schieve, L. A., Boyle, C. A., Perrin, J. M., Ghandour, R. M., et al. (2009). Prevalence of parent-reported diagnosis of autism spectrum disorder among children in the US, 2007. Pediatrics 124, 1395–1403. doi: 10.1542/peds.2009-1522
Krajmalnik-Brown, R., Lozupone, C., Kang, D. W., and Adams, J. B. (2015). Gut bacteria in children with autism spectrum disorders: challenges and promise of studying how a complex community influences a complex disease. Microb. Ecol. Health Dis. 26:26914. doi: 10.3402/mehd.v26.26914
Layton, B. A., Walters, S. P., Lam, L. H., and Boehm, A. B. (2010). Enterococcus species distribution among human and animal hosts using multiplex PCR. J. Appl. Microbiol. 109, 539–547. doi: 10.1111/j.1365-2672.2010.04675.x
Lukic, J., Strahinic, I., Milenkovic, M., Golic, N., Kojic, M., Topisirovic, L., et al. (2013). Interaction of Lactobacillus fermentum BGHI14 with rat colonic mucosa: implications for colitis induction. Appl. Environ. Microbiol. 79, 5735–5744. doi: 10.1128/AEM.0180713
Matamoros, S., Gras-Leguen, C., Le Vacon, F., Potel, G., and De La Cochetiere, M. F. (2013). Development of intestinal microbiota in infants and its impact on health. Trends Microbiol. 21, 167–173. doi: 10.1016/j.tim.2012.12.001
Merenstein, D. J., Tan, T. P., Molokin, A., Smith, K. H., Roberts, R. F., Shara, N. M., et al. (2015). Safety of Bifidobacterium animalis subsp. lactis (B. lactis) strain BB-12-supplemented yogurt in healthy adults on antibiotics: a phase I safety study. Gut. Microbes 6, 66–77. doi: 10.1080/19490976.2015.1005484
Miyauchi, E., Ogita, T., Miyamoto, J., Kawamoto, S., Morita, H., Ohno, H., et al. (2013). Bifidobacterium longum alleviates dextran sulfate sodium-induced colitis by suppressing IL-17A response: involvement of intestinal epithelial costimulatory molecules. PLoS ONE 8:e79735. doi: 10.1371/journal.pone.0079735
Molloy, C. A., and Manning-Courtney, P. (2003). Prevalence of chronic gastrointestinal symptoms in children with autism and autistic spectrum disorders. Autism 7, 165–171. doi: 10.1177/1362361303007002004
Nankova, B. B., Agarwal, R., Macfabe, D. F., and La Gamma, E. F. (2014). Enteric bacterial metabolites propionic and butyric acid modulate gene expression, including CREB-dependent catecholaminergic neurotransmission, in PC12 cells–possible relevance to autism spectrum disorders. PLoS ONE 9:e103740. doi: 10.1371/journal.pone.0103740
Neufeld, K. M., Kang, N., Bienenstock, J., and Foster, J. A. (2011). Reduced anxiety-like behavior and central neurochemical change in germ-free mice. Neurogastroenterol. Motil. 23, 255–264.e119. doi: 10.1111/j.1365-2982.2010.01620.x
Nikolov, R. N., Bearss, K. E., Lettinga, J., Erickson, C., Rodowski, M., Aman, M. G., et al. (2009). Gastrointestinal symptoms in a sample of children with pervasive developmental disorders. J. Autism. Dev. Disord. 39, 405–413. doi: 10.1007/s10803-008-0637-8
Nubel, U., Engelen, B., Felske, A., Snaidr, J., Wieshuber, A., Amann, R. I., et al. (1996). Sequence heterogeneities of genes encoding 16S rRNAs in Paenibacillus polymyxa detected by temperature gradient gel electrophoresis. J. Bacteriol. 178, 5636–5643.
Ohata, A., Usami, M., and Miyoshi, M. (2005). Short-chain fatty acids alter tight junction permeability in intestinal monolayer cells via lipoxygenase activation. Nutrition 21, 838–847. doi: 10.1016/j.nut.2004.12.004
O'mahony, S. M., Clarke, G., Dinan, T. G., and Cryan, J. F. (2017). Early-life adversity and brain development: is the microbiome a missing piece of the puzzle? Neuroscience 342, 37–54. doi: 10.1016/j.neuroscience.2015.09.068
Ouwehand, A. C., Salminen, S., Arvola, T., Ruuska, T., and Isolauri, E. (2004). Microbiota composition of the intestinal mucosa: association with fecal microbiota? Microbiol. Immunol. 48, 497–500. doi: 10.1111/j.1348-0421.2004.tb03544.x
Parracho, H. M., Bingham, M. O., Gibson, G. R., and Mccartney, A. L. (2005). Differences between the gut microflora of children with autistic spectrum disorders and that of healthy children. J. Med. Microbiol. 54(Pt 10), 987–991. doi: 10.1099/jmm.0.46101-0
Rajilic-Stojanovic, M., Smidt, H., and De Vos, W. M. (2007). Diversity of the human gastrointestinal tract microbiota revisited. Environ. Microbiol. 9, 2125–2136. doi: 10.1111/j.1462-2920.2007.01369.x
Ringel-Kulka, T., Cheng, J., Ringel, Y., Salojarvi, J., Carroll, I., Palva, A., et al. (2013). Intestinal microbiota in healthy U.S. young children and adults–a high throughput microarray analysis. PLoS ONE 8:e64315. doi: 10.1371/journal.pone.0064315
Rogers, G. B., Keating, D. J., Young, R. L., Wong, M. L., Licinio, J., and Wesselingh, S. (2016). From gut dysbiosis to altered brain function and mental illness: mechanisms and pathways. Mol. Psychiatry 21, 738–748. doi: 10.1038/mp.2016.50
Rossignol, D. A., and Frye, R. E. (2012). A review of research trends in physiological abnormalities in autism spectrum disorders: immune dysregulation, inflammation, oxidative stress, mitochondrial dysfunction and environmental toxicant exposures. Mol. Psychiatry 17, 389–401. doi: 10.1038/mp.2011.165
Russell, W. R., Duncan, S. H., Scobbie, L., Duncan, G., Cantlay, L., Calder, A. G., et al. (2013). Major phenylpropanoid-derived metabolites in the human gut can arise from microbial fermentation of protein. Mol. Nutr. Food Res. 57, 523–535. doi: 10.1002/mnfr.201200594
Sandler, R. H., Finegold, S. M., Bolte, E. R., Buchanan, C. P., Maxwell, A. P., Vaisanen, M. L., et al. (2000). Short-term benefit from oral vancomycin treatment of regressive-onset autism. J. Child Neurol. 15, 429–435. doi: 10.1177/088307380001500701
Satokari, R. M., Vaughan, E. E., Akkermans, A. D., Saarela, M., and De Vos, W. M. (2001). Bifidobacterial diversity in human feces detected by genus-specific PCR and denaturing gradient gel electrophoresis. Appl. Environ. Microbiol. 67, 504–513. doi: 10.1128/AEM.67.2.504-513.2001
Schieve, L. A., Gonzalez, V., Boulet, S. L., Visser, S. N., Rice, C. E., Van Naarden Braun, K., et al. (2012). Concurrent medical conditions and health care use and needs among children with learning and behavioral developmental disabilities, national health interview survey, 2006-2010. Res. Dev. Disabil. 33, 467–476. doi: 10.1016/j.ridd.2011.10.008
Scott, F. J., Baron-Cohen, S., Bolton, P., and Brayne, C. (2002). Brief report: prevalence of autism spectrum conditions in children aged 5-11 years in Cambridgeshire, UK. Autism 6, 231–237. doi: 10.1177/1362361302006003002
Shattock, P., and Whiteley, P. (2002). Biochemical aspects in autism spectrum disorders: updating the opioid-excess theory and presenting new opportunities for biomedical intervention. Expert. Opin. Ther. Targets 6, 175–183. doi: 10.1517/14728184.108.40.206
Sonnenburg, J. L., Chen, C. T., and Gordon, J. I. (2006). Genomic and metabolic studies of the impact of probiotics on a model gut symbiont and host. PLoS Biol. 4:e413. doi: 10.1371/journal.pbio.0040413
Stackebrandt, E., Sproer, C., Rainey, F. A., Burghardt, J., Pauker, O., and Hippe, H. (1997). Phylogenetic analysis of the genus Desulfotomaculum: evidence for the misclassification of Desulfotomaculum guttoideum and description of Desulfotomaculum orientis as Desulfosporosinus orientis gen. nov., comb. nov. Int. J. Syst. Bacteriol. 47, 1134–1139. doi: 10.1099/00207713-47-4-1134
Steppe, M., Van Nieuwerburgh, F., Vercauteren, G., Boyen, F., Eeckhaut, V., Deforce, D., et al. (2014). Safety assessment of the butyrate-producing Butyricicoccus pullicaecorum strain 25-3(T), a potential probiotic for patients with inflammatory bowel disease, based on oral toxicity tests and whole genome sequencing. Food Chem. Toxicol. 72, 129–137. doi: 10.1016/j.fct.2014.06.024
Sudo, N., Chida, Y., Aiba, Y., Sonoda, J., Oyama, N., Yu, X. N., et al. (2004). Postnatal microbial colonization programs the hypothalamic-pituitary-adrenal system for stress response in mice. J. Physiol. 558(Pt 1), 263–275. doi: 10.1113/jphysiol.2004.063388
Suzuki, T., Yoshida, S., and Hara, H. (2008). Physiological concentrations of short-chain fatty acids immediately suppress colonic epithelial permeability. Br. J. Nutr. 100, 297–305. doi: 10.1017/S0007114508888733
Suzuki-Mizushima, Y., Gohda, E., Okamura, T., Kanasaki, K., and Yamamoto, I. (2002). Enhancement of NGF- and cholera toxin-induced neurite outgrowth by butyrate in PC12 cells. Brain Res. 951, 209–217. doi: 10.1016/s0006-8993(02)03163-3
Swedo, S. E., Baird, G., Cook, E. H., Happe, F. G., Harris, J. C., Kaufmann, W. E., et al. (2013). “Neurodevelopmental disorders,” in American Psychiatric Association: Diagnostic and Statistical Manual of Mental Disorders, eds. S. K. Schultz and E. A. Kuhl (Washington, DC: American Psychiatric Association), 31–87.
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
Turroni, F., Peano, C., Pass, D. A., Foroni, E., Severgnini, M., Claesson, M. J., et al. (2012). Diversity of bifidobacteria within the infant gut microbiota. PLoS ONE 7:e36957. doi: 10.1371/journal.pone.0036957
Turroni, F., Taverniti, V., Ruas-Madiedo, P., Duranti, S., Guglielmetti, S., Lugli, G. A., et al. (2014). Bifidobacterium bifidum PRL2010 modulates the host innate immune response. Appl. Environ. Microbiol. 80, 730–740. doi: 10.1128/AEM.03313-13
Valvassori, S. S., Varela, R. B., Arent, C. O., Dal-Pont, G. C., Bobsin, T. S., Budni, J., et al. (2014). Sodium butyrate functions as an antidepressant and improves cognition with enhanced neurotrophic expression in models of maternal deprivation and chronic mild stress. Curr. Neurovasc. Res. 11, 359–366. doi: 10.2174/1567202611666140829162158
Van Immerseel, F., Ducatelle, R., De Vos, M., Boon, N., Van De Wiele, T., Verbeke, K., et al. (2010). Butyric acid-producing anaerobic bacteria as a novel probiotic treatment approach for inflammatory bowel disease. J. Med. Microbiol. 59(Pt 2), 141–143. doi: 10.1099/jmm.0.017541-0
Varela, R. B., Valvassori, S. S., Lopes-Borges, J., Mariot, E., Dal-Pont, G. C., Amboni, R. T., et al. (2015). Sodium butyrate and mood stabilizers block ouabain-induced hyperlocomotion and increase BDNF, NGF and GDNF levels in brain of Wistar rats. J. Psychiatr. Res. 61, 114–121. doi: 10.1016/j.jpsychires.2014.11.003
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
Wasilewska, J., and Klukowski, M. (2015). Gastrointestinal symptoms and autism spectrum disorder: links and risks - a possible new overlap syndrome. Pediatric Health Med. Ther. 6, 153–166. doi: 10.2147/PHMT.S85717
WHO (2015). “Mental and behavioral disorders: clinical descriptions and diagnostic guidelines,” in International Statistical Classification of Diseases and Related Health Problems (ICD10) (Geneva: WHO Library Cataloguing in Publication Data), 232–259.
Wu, X., Chen, P. S., Dallas, S., Wilson, B., Block, M. L., Wang, C. C., et al. (2008). Histone deacetylase inhibitors up-regulate astrocyte GDNF and BDNF gene transcription and protect dopaminergic neurons. Int. J. Neuropsychopharmacol. 11, 1123–1134. doi: 10.1017/S1461145708009024
Ze, X., Ben David, Y., Laverde-Gomez, J. A., Dassa, B., Sheridan, P. O., Duncan, S. H., et al. (2015). Unique organization of extracellular amylases into amylosomes in the resistant starch-utilizing human colonic firmicutes bacterium Ruminococcus bromii. MBio 6, e01058–e01015. doi: 10.1128/mBio.01058-15
Keywords: autism, gut-brain axis, Lactobacillus, Bifidobacterium, Clostridium like species, microbial diversity
Citation: Bojović K, Ignjatović Ð, Soković Bajić S, Vojnović Milutinović D, Tomić M, Golić N and Tolinački M (2020) Gut Microbiota Dysbiosis Associated With Altered Production of Short Chain Fatty Acids in Children With Neurodevelopmental Disorders. Front. Cell. Infect. Microbiol. 10:223. doi: 10.3389/fcimb.2020.00223
Received: 02 October 2019; Accepted: 22 April 2020;
Published: 19 May 2020.
Edited by:Andrew T. Gewirtz, Georgia State University, United States
Reviewed by:Kun Zhang, Virginia Commonwealth University, United States
Junguk Hur, University of North Dakota, United States
Copyright © 2020 Bojović, Ignjatović, Soković Bajić, Vojnović Milutinović, Tomić, Golić and Tolinački. 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.
†These authors have contributed equally to this work