Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Comparison of Fecal Microbiota in Children with Autism Spectrum Disorders and Neurotypical Siblings in the Simons Simplex Collection

  • Joshua S. Son,

    Affiliation Department of Medicine, Stony Brook University, Stony Brook, NY, United States of America

  • Ling J. Zheng,

    Affiliation Department of Medicine, Stony Brook University, Stony Brook, NY, United States of America

  • Leahana M. Rowehl,

    Affiliation Department of Medicine, Stony Brook University, Stony Brook, NY, United States of America

  • Xinyu Tian,

    Affiliation Department of Applied Mathematics and Statistics, Stony Brook University, Stony Brook, NY, United States of America

  • Yuanhao Zhang,

    Affiliation Department of Applied Mathematics and Statistics, Stony Brook University, Stony Brook, NY, United States of America

  • Wei Zhu,

    Affiliation Department of Applied Mathematics and Statistics, Stony Brook University, Stony Brook, NY, United States of America

  • Leighann Litcher-Kelly,

    Affiliation Department of Psychiatry, Stony Brook University, Stony Brook, NY, United States of America

  • Kenneth D. Gadow,

    Affiliation Department of Psychiatry, Stony Brook University, Stony Brook, NY, United States of America

  • Grace Gathungu,

    Affiliation Department of Pediatrics, Stony Brook University, Stony Brook, NY, United States of America

  • Charles E. Robertson,

    Affiliation Department of Medicine, University of Colorado Anschutz Medical Campus, Aurora, CO, United States of America

  • Diana Ir,

    Affiliation Department of Medicine, University of Colorado Anschutz Medical Campus, Aurora, CO, United States of America

  • Daniel N. Frank,

    Affiliation Department of Medicine, University of Colorado Anschutz Medical Campus, Aurora, CO, United States of America

  • Ellen Li

    ellen.li@stonybrookmedicine.edu

    Affiliation Department of Medicine, Stony Brook University, Stony Brook, NY, United States of America

Abstract

In order to assess potential associations between autism spectrum disorder (ASD) phenotype, functional GI disorders and fecal microbiota, we recruited simplex families, which had only a single ASD proband and neurotypical (NT) siblings, through the Simons Simplex Community at the Interactive Autism Network (SSC@IAN). Fecal samples and metadata related to functional GI disorders and diet were collected from ASD probands and NT siblings of ASD probands (age 7–14). Functional gastrointestinal disorders (FGID) were assessed using the parent-completed ROME III questionnaire for pediatric FGIDs, and problem behaviors were assessed using the Child Behavior Check List (CBCL). Targeted quantitative polymerase chain reaction (qPCR) assays were conducted on selected taxa implicated in ASD, including Sutterella spp., Bacteroidetes spp. and Prevotella spp. Illumina sequencing of the V1V2 and the V1V3 regions of the bacterial 16S rRNA genes from fecal DNA was performed to an average depth of 208,000 and 107,000 high-quality reads respectively. Twenty-five of 59 ASD children and 13 of 44 NT siblings met ROME III criteria for at least one FGID. Functional constipation was more prevalent in ASD (17 of 59) compared to NT siblings (6 of 44, P = 0.035). The mean CBCL scores in NT siblings with FGID, ASD children with FGID and ASD without FGID were comparably higher (58–62 vs. 44, P < 0.0001) when compared to NT children without FGID. There was no significant difference in macronutrient intake between ASD and NT siblings. There was no significant difference in ASD severity scores between ASD children with and without FGID. No significant difference in diversity or overall microbial composition was detected between ASD children with NT siblings. Exploratory analysis of the 16S rRNA sequencing data, however, identified several low abundance taxa binned at the genus level that were associated with ASD and/or first order ASD*FGID interactions (FDR <0.1).

Introduction

Autism spectrum disorders (ASD) are neurodevelopmental disorders characterized primarily by alterations in reciprocal social interactions and communication as well as repetitive behavior. The reported prevalence of ASD has increased and affects 1 in 68 children in the USA [1]. This may be due in part to heightened awareness for these disorders [2,3]. The diagnosis of ASD is based on specially trained physicians and psychologists administering autism-specific behavioral evaluations. Twin studies provide support for both genetic and environmental contributions to ASD pathogenesis [2,3].

Alterations in the composition of resident gut microbial communities have been reported for an increasing number of human disorders, including ASD and functional GI disorders (FGID) [424]. The term “functional” in FGIDs refers to the absence of structural abnormalities in the gut that can be identified by blood, radiographic or endoscopic tests [25,26]. Consequently, FGIDs are diagnosed based on assessment of GI symptoms, and the exclusion of other GI diseases (e.g. celiac disease, inflammatory bowel diseases). FGIDs are very common in the general population and comprise a major portion of gastroenterology and primary care practices. Abnormalities in the neural regulation of gut function and sensation are implicated in the pathogenesis of FGIDs [20,25,26].

Functional GI symptoms have been frequently reported in ASD children as a non-neurological co-morbidity [2731]. However a major barrier to integrating the ASD studies on GI symptoms with those conducted in the general population has been that the instruments used to assess functional GI symptoms in ASD differed from the ROME III process, which is the current international standard by which functional GI disorders (FGID) are diagnosed and classified [24].

In 2010, the Simons Foundation Autism Research Initiative (SFARI) began assembling the Simons Simplex Collection (SSC) as a resource for ASD research [32]. The SSC includes over 2000 simplex families, each with a single ASD proband but with unaffected parents and siblings. The simplex families were recruited through 12 university affiliated research clinics located in Canada and the United States. Upon recruitment to the SSC, the ASD probands (age 4–18) were evaluated with an extensive battery of diagnostic measures including the Autism Diagnostic Observation Schedule (ADOS) using the same diagnostic criteria across centers. Medical data, including GI symptoms, were collected on the ASD probands, but not on the NT siblings, using an extensive SSC Medical History form. Families were excluded if the proband did not meet criteria for ASD, had medically significant perinatal incidents (e.g. perinatal asphyxia), low mental age, or had primary relatives on the autism spectrum. The SSC thus represents a “unique, well-described sample of able children and adolescents with relatively severe ASD, as indicated by ADI-R and ADOS Calibrated Severity Scores (ADOS-CSS)” [32]. In order to assess the associations between ASD phenotype, FGIDs and fecal microbiota, we recruited a subset of the SSC families with ASD probands and/or NT siblings (age 7–14), to participate in a study collecting fecal samples and metadata related to functional GI disorders and diet.

Materials and Methods

Patient recruitment

Families were recruited via a registry called the Simons Simplex Community through the Interactive Autism Network (SSC@IAN), which is composed of simplex families originally recruited to the Simons Simplex Collection that were willing to be contacted through the Interactive Autism Network for additional studies. The SSC@IAN provided us with contact information for 245 potentially eligible simplex families that expressed their willingness to be contacted about this study. Of the 245 potentially eligible SSC families with 339 children, 107 families with 145 children consented to participate with parental written consents and written child assents for all participating children, and 66 families with 103 children (66 ASD probands and 37 NT siblings) completed the study. Exclusion criteria for ASD probands and NT siblings in this study were age <7 or age >14 years, a past or current history of GI diseases or other serious medical problems, and use of antibiotics within one month of collection of the stool sample. The procedures of the present study were reviewed and approved by the Institutional Review Board of Stony Brook University (IRB 313466)

Families who consented to participate in the study were asked to complete the Pediatric ROME III Version (QPGS-RIII), the GI symptoms form in the Simon Simplex Collection Medical History (v 4.0) [31], the Child Behavior Checklist (CBCL/6-18) [3337], a one-week food diary for the ASD proband and/or the NT sibling prior to collecting a fecal sample, and provide current information on the child’s height and weight. The families were asked to submit the stool samples after the subjects were off probiotics and antibiotics for at least one month.

Pediatric ROME III process of diagnosing FGID and the GI symptoms section of the SSC medical history form

Functional GI disorders (FGID) were assessed using the ROME III process [25]. Information on GI symptoms in both ASD and NT siblings was also collected using the GI symptoms section of the SSC Medical History form [32]. The GI symptoms section of the SSC Medical History form included the seven symptom categories: constipation, diarrhea, abnormal stool smell, bloating, excessive gas, severe abdominal pain, and vomiting. The parents were asked to rank whether these categories were problematic as “yes”, “not sure” and “no”. Many of these symptom categories in the SSC form overlap with the following symptoms included in the calculation of the GI severity index score [27]: constipation, diarrhea, stool smell, flatulence and abdominal pain, but did not include “unexplained daytime irritability”, “nighttime awakening” or “abdominal tenderness during exam”.

Child Behavior Checklist (CBCL)

The CBCL/6-18 is a parent report form to screen for emotional, behavioral, and social problems in children ages 6–18 [33]. This instrument as well as the earlier version of this instrument has been previously used by several research groups studying school-aged children with ASD [3437]. The data from ASD and NT siblings was entered and analyzed using the AGES 6–18 ASEBA-PC MODULE (ASEBA) to yield norm-referenced T-scores (M = 50, SD = 10). The CBCL contains eight syndrome scales: Anxious/Depressed, Withdrawn/Depressed, Somatic Complaints, Rule Breaking Behavior, Aggressive Behavior. The Internalizing Domain score is a measure of emotional problems and contains three syndrome scales: Anxious/Depressed, Withdrawn/Depressed, and Somatic Complaints. The Externalizing Domain score is a measure of behavioral problems and contains the Rule Breaking Behavior and Aggressive Behavior syndrome scales.

Dietary information

The parents were asked to record the daily dietary intake of ASD and NT siblings on the Stony Brook University Medical Center Department of Pediatrics Food Diary/Calorie Count sheet for 7 days prior to the collection of the stool sample. The data were entered and analyzed using Nutritionist Pro nutrient analysis software version 5.2.0 (Axxya Systems, Nutritionist Pro, Stafford, TX). We focused on analyzing the following six macronutrient categories: total kcal/day, total protein g/day, total carbohydrates g/day, total fat g/day, total dietary fiber/day, total sugar g /day.

Fecal sample collection and DNA extraction

Fecal samples were collected in a stool “hat” placed in the toilet and following evacuation 2 ml of stool were immediately transferred, using the spoon built into the lid into a Para-Pak Clean vial (Meridian Bioscience, Cincinnati, OH) that contained 10 ml of RNAlater for metagenomics studies and a Para-Pak Clean vial with no preservative. The fecal samples were shipped in cold packs overnight to our laboratory. Fecal DNA was extracted from stool samples immediately upon arrival using ZR Fecal DNA MiniPre (Zymo Research Corporation, Irvine, CA) and following manufacturer’s protocol. Fecal DNA and stool samples were stored in -80°C.

Quantitative PCR (qPCR) for targeted bacterial subgroups

QPCR assays were performed using established primers for total bacteria (forward, 5’- GTG STG CAY GGY TGT CGT CA-3’ and reverse 5’- ACG TCR TCC MCA CCT TCC TC-3’) [38], Sutterella subgroup (forward, 5’-CGC GAA AAA CCT TAC CTA GCC-3’ and reverse 5’- GAC GTG TGA GGC CCT AGC C-3’) [10], Bacteroidetes subgroup (forward 5'-AACGCTAGCTACAGGCTT-3’ and reverse 5'-CCAATGTGGGGGACCTTC-3’) [39]; Prevotella (forward, 5’-CACCAAGGCGACGATCA-3’ and reverse 5’-GGATAACGCYGGACCT-3’) [40], C. coccoides-E. rectales subgroup (forward, 5’—CGGTACCTGACTAAGAAGC-3’and reverse 5’-AGTTT(C/T)ATTCTTGCGAACG-3’) [41], Faecalibacterium prausnitzii (forward 5’-CCC TTC AGT GCC GCA GT-3’ and reverse 5’-GTC GCA GGA TGT CAA GAC-3’) [42], Escherichia coli subgroup (forward, 5’- GTT AAT ACC TTT GCT CAT TGA-3’ and reverse, 5’- ACC AGG GTA TCT AAT CCT GTT-3’) [42]. The log2 transformation of the relative abundance of each bacterial subgroup was measured by ΔCt = Ct (threshold cycle) total bacteria−Ct subgroup. All assays were carried out in triplicate. Plasmid quantification standards were prepared from representative clones of the target organisms to insure that the assays were conducted within the linear range and with similar slopes.

16S rRNA amplicon library construction and Illumina V1V2 and V1V3 sequencing analysis

Fecal bacterial profiles were determined using broad-range amplification of the ~300 base pair (bp) V1V2 variable region using primers 27FYM (AGAGTTTGATYMTGGCTCAG) and 338R (TGCTGCCTCCCGTAGGAGT) and the ~500 bp V1V3 region using primers 27FYM (AGAGTTTGATYMTGGCTCAG) and 534R (ATTACCGCGGCKGCTGG) [43,44]. Sequencing was performed on the Illumina MiSeq platform using 600-cycle version 3 reagent kits. Demultiplexed, paired end sequence data were deposited in the NCBI Short Read Archive under BioProject Accession Number: PRJNA282013 (www.ncbi.nlm.nih.gov/bioproject/PRJNA2282203). The sorted paired reads were assembled using phrap [45] and paired reads that did not assemble were discarded. Assembled sequence ends were trimmed over a moving window of 5 nucleotides until average quality met or exceeded 20. Trimmed contigs with more than 1 ambiguity or shorter than 200 nucleotides were discarded. Potential chimeras identified with Uchime (usearch6.0.203_i86linux32) [46] using the Schloss [47] Silva reference sequences were removed from subsequent analyses. Assembled sequences were aligned and classified with SINA (1.2.11) using the 418,497 bacterial sequences in Silva 115NR99 as reference configured to yield the Silva taxonomy [48,49]. Operational taxonomic units (OTUs) were produced by clustering sequences with identical taxonomic assignments. OTU counts were normalized between samples by dividing sequence counts by the total number of sequences generated per sample to calculate the relative abundance. Phylum-level OTU tables were generated by collapsing lower level OTUs into higher-level categories. OTUs with a maximal relative abundance <0.0001 and with a prevalence <0.01 were culled. The remaining OTU relative abundances were then transformed using the square root function.

Statistical analysis

Comparisons of gender, FGIDs, and SSC GI symptoms between ASD and NT siblings were conducted using the unconditional Fisher exact test with one tail [50]. The threshold significance was set at P < 0.05.

The ASD and NT siblings were further subdivided with respect to whether they also carried a diagnosis of an FGID based on the ROME III process. Study characteristics including body mass index (BMI), CBCL scores, dietary scores and the qPCR ΔCt scores were compared between the following four groups: 1) ASD with FGID, 2) ASD without FGID, 3) NT siblings with FGID and 4) NT siblings without FGID. Outcome variables including CBCL, dietary macronutrients, qPCR ΔCt were also modeled using mixed effect models with ASD, FGID and ASD*FGID as the three fixed effects, and family as random effects. Each continuous variable was fitted to a linear model while each categorical variable was fitted to a logistic model. Comparison of continuous variables between family matched ASD and NT siblings was conducted using two-tailed Student’s T test or the Wilcoxon signed rank test with the threshold of significance set as P < 0.05.

Alpha diversity indices (e.g. Chao1, Shannon complexity H, Shannon Evenness H/Ho) were calculated using Explicet [51] for each of the four groups with 10,000 replicate resamplings: 1) ASD with FGID, 2) ASD without FGID, 3) NT siblings with FGID and 4) NT siblings without FGID.

Beta diversity was compared between ASD and NT siblings using the adonis function in the R vegan package as previously described [52,53] at the phyla, family and genus level. This function uses a non-parametric multivariate analysis of variance test (PERMANOVA) [54], and was applied using Bray-Curtis, Jaccard and Morisita-Horn indices as distance measurements [55]. Family (i.e. whether ASD or NT siblings belonged to the same family) was selected as a strata variable to restrict permutations as previously described [56]. Overall microbial composition binned as phyla, was also compared between family matched ASD and NT siblings using the permutation Hotelling T test [57].

Paired two class comparisons of the transformed (square root) relative abundances of individual OTUs (phyla and genera level) between ASD probands and NT siblings from the same family were conducted using the Wilcoxon signed rank test with the threshold set as FDR < 0.05.

Because over-dispersion is often observed in microbiome sequence count data [58], a negative binomial regression model was also used to assess the effect of the ASD phenotype, FGID phenotype and first order interactions between ASD and FGID (ASD*FGID) on each individual OTU, in both paired and unpaired ASD and NT samples as follows: where Yijk denotes the raw sequence count of OTU k for child j in family i. In this model we treat family, which may reflect genetic and possibly dietary factors, as a zero mean random coefficient. Coefficients b0, β1, β2 and β3 are fixed constants representing grand mean, the effect of ASD, FGID and the interaction of ASD and FGID respectively. The log total sequence count () for each individual is considered as an offset. Because the distribution of the raw sequence count is often skewed by a large proportion of zero counts, a zero-inflated version of the negative binomial model was also fitted to each individual OTU: The final model was selected based on the Akaike Information Criterion (AIC). The p-values for each OTU were adjusted for multiple comparisons to calculate a false discovery rate (FDR) using the Benjamini-Hochberg procedure [59]. The threshold of significance for ASD, FGID or ASD*FGID first order interactions was set at FDR < 0.1. The logistic and linear mixed models were built with R package glmmADMB [60] and lme4-R [61] respectively.

Comparative analysis of the parallel square root transformed relative abundance [62] measurements of Bacteroidetes phylum, the Sutterella genus and the Prevotella genus by MiSeq V1V2, MiSeq V1V3 and targeted qPCR was conducted using latent structural equation modeling [63]. Because of the large proportion of zero counts for Sutterella and Provetella genera, samples with zero counts were excluded from the analysis. To avoid the complication caused by correlated samples, SEM models were fitted with grouping variable ASD and the constraint that loadings are equal for different groups. The corresponding model for each group is: The corresponding model for each group is: Here ξij is the underlying true bacteria relative abundance in the square root scale, while δ1, δ2 and δ3 are the corresponding measurement errors. In order to render the model identifiable, we chose qPCR as the reference level with a1 = 0 and b1 = 1. The model was fitted with the R package lavaan [64].

Results

Demographics and functional GI disorders (FGID) in ASD and NT siblings

A total of 59 ASD children and 44 NT siblings from 66 families completed the study (Table 1). Matched ASD and NT sibling pairs were from 37 families. In some instances NT siblings of ASD probands were recruited, but not the ASD proband s because their age exceeded 14 years. The ASD children were predominantly male. The gender distribution for the NT siblings was roughly equal. The mean age in both ASD children and NT siblings was 10 years, and they were predominantly of Caucasian descent. There are three major groups of FGIDs (see Table 1) [24]: H1) Vomiting and aerophagia, H2). Abdominal-pain related FGIDs, and H3) Constipation and fecal incontinence. Group H1 is further subdivided into three FGIDs: H1a) Adolescent rumination syndrome; H1b) Cyclic vomiting syndrome; H1c) Aerophagia. Group H2 is further subdivided into four FGIDS: H2a) Functional dyspepsia, H2b) Irritable bowel syndrome (IBS), H2c) abdominal migraine, H2d) Childhood functional abdominal pain with a subset termed H2d1. Childhood functional abdominal pain syndrome. Subjects classified within the H2d1 subset have additional somatic symptoms such as headache, limb pain or difficulty sleeping.

thumbnail
Table 1. Demographics and functional GI disorders (FGID) in ASD and NT siblings.

The study characteristics of the ASD and NT siblings are listed below. For categorical variables, the percent are shown in parenthesis (). For continuous variables the mean ± standard deviations are listed.

https://doi.org/10.1371/journal.pone.0137725.t001

Twenty-five of 59 (42%) ASD children and 13 of 44 (30%) NT siblings were diagnosed with at least one FGID. Fifteen ASD children and 13 NT siblings in 37 families met ROME III criteria for at least one FGID, with concordance for FGIDs in 6 (18%) ASD/NT sibling pairs. Some ASD and NT siblings were diagnosed with more than one FGID. ASD children exhibited a higher prevalence of functional constipation in compared to NT siblings (P = 0.035) when both matched and unmatched individuals were included in the analysis. We detected no significant difference in the prevalence of children with at least one FGID in family matched ASD/NT sibling pairs. We also detected no significant difference in the prevalence of specific FGIDs such as functional constipation in the family matched ASD/NT sibling pairs.

When the SSC registry was established, medical history forms including one on GI symptoms were completed by the parents on the ASD children but not the NT siblings. The same SSC GI symptoms form was then completed (~ 2–3 years later) by the parents for the ASD and NT children for the current study. As shown in Table 2, there was no significant change in the prevalence of the GI symptoms occurred over time. Furthermore, the SSC results indicate that ASD children exhibited a higher prevalence of constipation compared to NT siblings (P = 0.05).

thumbnail
Table 2. SSC problematic GI symptoms results.

The previous and current percent of ASD children and NT siblings having a problematic GI symptom is listed below. The previous percent were based on parental responses to the SSC medical form for ASD children only on first entering the SSC registry. The current percent were based on parental responses to the same forms for both ASD and NT siblings during the current study ~2–3 years later. The P-values reflect comparison of the current responses of the ASD and NT siblings. There was no significant change in the prevalence of any of the seven problematic GI symptoms in the ASD children at the two time points.

https://doi.org/10.1371/journal.pone.0137725.t002

Comparison of CBCL and macronutrient intake

As shown in Table 3 and S1 Table, the ASD severity score, ADOS-CSS, did not differ significantly between the ASD w FGID group and the ASD w/o FGID group. Mean total, internalizing, externalizing and subscale CBCL scores were significantly increased (p < 0.05) in the ASD children in both the ASD/NT sibling pairs and in comparison of all the ASD with NT children (matched and unmatched). These CBCL scores were also significantly increased (p < 0.05) in NT siblings with FGID compared to NT siblings without FGID (Table 3). There was no significant difference in the mean CBCL scores between ASD children with FGID compared with ASD children without FGID.

thumbnail
Table 3. Comparison of ADOS-CSS, CBCL scores and macronutrient intake between ASD children with FGID, ASD children without FGID, NT siblings with FGID and NT siblings without FGID.

For continuous variables the mean values ± standard deviations are listed. N.A., not available.

https://doi.org/10.1371/journal.pone.0137725.t003

No significant differences (p < 0.05) between the four study groups with respect to daily intake of macronutrients (calories, protein, fats, carbohydrates, sugars or dietary fiber) were noted (Table 3). Only 4 children were reported to adhere to a gluten and/or casein free diet. The study participants had BMIs that were within the normal range and were not obese.

No associations (FDR < 0.1) detected in targeted bacterial subgroups between ASD, FGID and/or first order ASD*FGID by QPCR

Targeted qPCR assays of fecal samples were conducted for Sutterella, Prevotella and total Bacteroidetes subgroups, because of previous reports that the relative abundance of these subgroups differed between ASD and NT children [7,8,1014]. QPCR assays targeting the C. coccoides-E-rectales group, Faecalibacterium prausnitzii and Escherichia coli were also performed because altered relative abundances have been reported previously in GI disorders such as inflammatory bowel diseases [23]. The mean ΔCt values of these subgroups, which represent a log2 transformation of the relative abundance, for 1) ASD w FGID, 2) ASD w/o FGID, 3) NT w FGID and 4) NT w/o FGID are shown in Table 4. No significant difference was observed for any of these bacterial subgroups in the 37 ASD/NT sibling pairs (see Methods). Furthermore, in mixed models, no significant effect of ASD, FGID or ASD*FGID was observed in an analysis of all 59 ASD and 44 NT siblings.

thumbnail
Table 4. QPCR comparisons of selected bacterial subgroups in ASD children with FGID, ASD children without FGID, NT siblings with FGID and NT siblings without FGID.

The qPCR assays were conducted using established primers as described in Methods. The mean ΔCt values (~Log2 relative abundance of targeted bacterial subgroups) ± standard deviations are listed. The data from family matched and unmatched ASD and NT siblings are included.

https://doi.org/10.1371/journal.pone.0137725.t004

16S rRNA sequence analysis of ASD children and NT siblings

We performed Illumina sequencing of the V1V3 hypervariable region of the bacterial 16S rRNA gene to facilitate comparisons of this data with previous 16S rRNA 454 pyrosequencing datasets of the V1V3 hypervariable region [16,23]. A total of 11,050,770 high-quality sequences were generated (median 1.1E+05 sequences/sample, IQR 8.0e+04 to 1.2E+05). All libraries had a Good’s coverage score ≥ 99.9% at the rarefaction point of 15,000 sequences, indicating deep sequence coverage of the intestinal microbiome was achieved for each sample. To corroborate results generated through V1V3 sequencing, we also performed Illumina sequencing of the V1V2 hypervariable region. A total of 14,783,994 high-quality sequences were generated (median 1.4E+05 sequences/sample, IQR 1.3e+05 to 1.6E+05; minimum Good’s coverage = 99.9%). No significant effect of ASD, FGID or ASD*FGID was observed on OTU complexity or richness, two measures of α-diversity, Shannon diversity and Chao 1 (see Table 5 for mean values). No significant difference (p < 0.05) was noted in β-diversity between ASD and NT sibling bacterial communities analyzed at the phyla, family or genus level and using Bray-Curtis, Jaccard or Morisita-Horn indices as dissimilarity measurements.

thumbnail
Table 5. Shannon H and SChao1 scores in ASD children with FGID, ASD children without FGID, NT siblings with FGID and NT siblings without FGID.

The mean values for each of these scores ± standard deviation are shown below for the V1V2 and V1V3 datasets. The data from family matched and unmatched ASD children and NT siblings are included.

https://doi.org/10.1371/journal.pone.0137725.t005

The relative abundances of the major four phyla, (Firmicutes, Bacteroidetes, Actinobacteria, Proteobacteria) for the ASD children with FGID, ASD children without FGID, NT siblings with FGID and NT siblings without FGID groups for the V1V2 and V1V3 datasets are shown in Fig 1 (see S2 Table). Phyla of lower abundance, including RF3, Tenericutes, Cyanobacteria, Verrucomicrobia, Lentisphaerae, and Fusobacterium were also detected (see S3 Table). Permutational Hotelling T test revealed no significant difference in the overall microbial composition when binned at the phyla level, between ASD children and NT siblings. With the exception of the Verrucomicrobia phylum, the relative abundances of the phyla were similar between the V1V2 and V1V3 data sets. Exploratory analyses of individual phyla detected no significant effect (FDR <0.1) of ASD, FGID or ASD*FGID first order interactions for any of the individual phyla, including Bacteroidetes, using a negative binomial model. At the genus level, 359 operational taxonomic units (OTU) were detected for the V1V2 dataset and 189 OTUs for the V1-V3 datasets, respectively. No significant associations were detected between the relative abundance of the Sutterella or Prevotella genera and ASD, FGID or ASD*FGID first order interactions (see Table 6) in either the V1V2 or V1V3 datasets. Exploratory analyses of individual genera revealed a significant effect (FDR < 0.1) of ASD and ASD*FGID first order interactions on the Cyanobacteria/Chloroplast genus for both the V1V2 and/or V1V3 datasets using a zero-inflated negative binomial model (see Table 6). The mean relative abundance of this genus was increased in the ASD-with FGID group compared to those of the other three categories (ASD without FGID, NT siblings with FGID and NT siblings without FGID, see S4 Table for mean relative abundances). 16S RNA sequence analyses of the human gut microbiota have previously detected sequences binned within the Cyanobacteria phylum. These sequences could represent chloroplast 16S rRNA genomic sequences present in dietary plant foods or novel uncultivated bacterial species. Metagenomic studies have detected a novel nonphotosynthetic taxa termed “Melainabacteria” present in human gut flora that represents either a sibling phylum or a class within the Cyanobacteria phylum [65,66]. The Silva database used in the current study, however, binned the Melainabacteria sequences in a taxon that was distinct from the Chloroplast taxon. Inspection of the daily food diaries of all the subjects revealed that two ASD children with FGID (functional constipation), who had high relative abundances of this genus and consequently dominated the signal associated with ASD with FGID category, were fed chia seeds added to smoothies. None of the other subjects consumed chia seeds except one ASD child who ingested a serving of mixed grain bread that included chia grain. Of note, chia seeds have been touted as an autism cure on the internet [67,68].

thumbnail
Fig 1. Phyla comparison between ASD children with FGID, ASD children without FGID, NT siblings with FGID and NT siblings without FGID, using Illumina V1-V2 and V1-V3 datasets.

The average relative abundance of each phyla (relative abundance ≥ 0.02) is shown for ASD w FGID, ASD w/o FGID, NT w FGID and NT w/o FGID for each of the two sequencing datasets (See S1 and S2 Tables for means ± standard deviations).

https://doi.org/10.1371/journal.pone.0137725.g001

thumbnail
Table 6. Low abundance taxa binned at the genera level that are significantly associated with ASD, FGID or ASD*FGID first order interactions.

The FDRs for each low abundance taxa are listed below for the V1V2 and V1V3 datasets, with the threshold of significant <0.1. N.D. is not detected based on filtering criteria. The FDRs are also listed for two genera, Sutterella and Provetella, which have been previously reported to be increased and decreased, respectively, in ASD children compared to unrelated control children [10,11].

https://doi.org/10.1371/journal.pone.0137725.t006

For the V1V2 dataset only, significant associations were observed for 1) ASD and ASD*FGID first order interactions with Firmicutes/Asteroleplasma, 2) ASD, FGID and ASD*FGID first order interactions with Proteobacteria/Thalassospira, and 3) FGID and ASD*FGID first order interactions with Proteobacteria/Burkholderia (see Table 5 and S4 Table). The Firmicutes/Asteroleplasma genus appears to represent a group of uncultured wall-less mycobacterium like organisms within the class Erysipelotrichia [69]. The Proteobacteria/Thalassospira genus appears to represent alpha-proteobacteria present in bacterioplankton and a number of uncultivated organisms [70]. The sequences binned in this Proteobacteria/Burkholderia taxon represent unclassified sequences within the Burkholderia order that did not bin within the Comamonadaceae family (see below).

For the V1V3 dataset only, significant associations were observed for 1) ASD and ASD*FGID first order interactions with the Proteobacteria/Comamonadaceae, 2) ASD, FGID and ASD*FGID first order interactions with Fusobacteria.Fusobacteriales, 3) ASD*FGID first order interactions with Bacteroidetes/Prevotellaceae, and 4) ASD and ASD*FGID interactions with Actinobacteria/Mobiluncus. These associations, however, did not reach significance in the V1V2 dataset. Proteobacteria/Comamonadaceae is a large and diverse bacterial family within the Burkholderiales order (see above). The sequences binned within this Fusobacteria/Fusobacteriales taxon represent unclassified sequences within the Fusobacteriales family that did not bin within the Fusobacterium genus. The sequences binned within the Bacteroidetes/Prevotellaceae taxon represent unclassified sequences within the Prevotellaceae family that did not bin within the Prevotella genus. The Actinobacteria/Mobiluncus genus has previously been associated with bacterial vaginosis [71].

Comparison of MiSeq V1V2, V1V3 and qPCR measurements of selected bacterial categories

Comparisons of the relative abundances of the Bacteroidetes phylum, and the Sutterella and Prevotella genera reveal high correlations between MiSeq V1V2 and V1V3 (ρ≈0.9) and small proportional biases (S5 Table and S1 Fig). The different reverse primers, 338R and 534R, used for amplifying respectively the V1V2 and V1V3 regions, might have subtle differences in specificity. While there were generally positively correlated measurements between the MiSeq and qPCR platforms, we observed large systematic biases particularly with the two relative low abundance genera, Sutterella and Provetella. This may reflect at least in part that the depth of sequencing was not sufficient to accurately measure relatively low abundance bacterial categories.

Discussion

This study of 59 ASD children and 44 NT siblings in the Simons Simplex Community through the Interactive Autism Network (SSC@IAN) registry, was motivated by previous reports linking ASD, GI symptoms with alterations in the gut microbiota (see Fig 2) and some simply linking ASD with the gut microbiota (see Fig 3). These previous studies of gut microbial communities in ASD children used both unrelated control children and NT siblings of ASD children. One study also included related disorders, such as Aspergers or Pervasive developmental disorder not otherwise specified (PDD-NOS) as the non-ASD control population. Our current study compared microbial composition and diversity between ASD children with only NT siblings. Of the 44 NT siblings, 37 were family matched to the ASD children in the study and the remaining seven NT children were siblings of ASD children who were not in the current study, in some cases because the ASD child was older than 14 years of age. This study did not replicate previously reported differences in the relative abundances of bacterial categories, including the Bacteroidetes phylum [1214], the Sutterella genus [7,8,10] and the Prevotella genus [11] between ASD and unrelated control children. This may reflect the smaller number of subjects enrolled in the previous studies, possibly resulting in Type 1 errors. One study used intestinal mucosal [10] rather than fecal samples. Of note, previous studies [68,12,14], with the exception of one study [13], failed to detect differences in microbial composition between ASD and NT siblings [11]. Our study also detected no differences in either α or β diversity between ASD and NT siblings.

thumbnail
Fig 2. Summary of studies comparing gut microbiota in ASD and NT children with and without GI symptoms.

ASD w GI, ASD children with GI symptoms, ASD w/o GI, ASD children without GI symptoms, NT Sib w GI, NT siblings with GI symptoms; NT sib w/o GI, NT siblings without GI symptoms, NT Un Con w GI, unrelated controls with GI symptoms; Un Con w/o GI, unrelated controls without GI symptoms; PDD-NoS, pervasive developmental disorder not otherwise specified.

https://doi.org/10.1371/journal.pone.0137725.g002

thumbnail
Fig 3. Summary of studies comparing gut microbiota in ASD and NT children with no information on GI symptoms.

NT Sib, neurotypical sibling, Un Con, unrelated control.

https://doi.org/10.1371/journal.pone.0137725.g003

In this study, 16S rRNA sequencing using MiSeq was carried out to a greater depth than previous studies utilizing pyrosequencing. Both platforms detected an increased relative abundance of the low abundance Chloroplast genus in ASD compared to NT siblings. Inspection of the daily dietary logs of the ASD children with high counts of Chloroplast sequences revealed that these particular children were fed shakes containing chia seeds. We speculate that the observed increase in the relative abundance of Chloroplast sequences could reflect ingestion of chia seeds by these children.

An alternative explanation for why a study comparing ASD children with NT siblings did not replicate differences reported between ASD children with unrelated controls, is that NT siblings of ASD children have altered microbiomes compared to that of unrelated children [24]. Studies tracking siblings of ASD children have detected vulnerabilities in some neurocognitive domains in the absence of an ASD diagnosis [72] suggesting that some of the NT siblings may exhibit a broader subclinical ASD phenotype.

Review of previous studies suggest that GI symptoms are more common in ASD children than in NT children, but these studies suffer from methodological differences in assessing GI symptoms [2831]. Our assessment of GI symptoms differs from previous studies in that we used the international standard for assessing FGIDs in the general pediatric population. Because of potential overlaps between ROME III FGID diagnoses, we subdivided the ASD children and NT siblings each into two categories 1) those with at least one FGID and 2) those without any FGID. A high prevalence of children diagnosed with at least one FGID was observed in both the ASD (~40%) and NT siblings (~30%) in the current study. Previous epidemiological studies using the ROME III process to survey the prevalence of a child having at least one FGID in general pediatric clinics reported a prevalence of 2% and 16%, respectively [73,74], suggesting that NT siblings may also have a higher prevalence of FGID than the general pediatric population.

The current study detected an increased prevalence of functional constipation FGID in ASD children compared to NT siblings reached significance (see Tables 1 and 2). Constipation has been previously recognized as a significant comorbidity in ASD children [2831]. Further studies on a larger sample of ASD and NT siblings (preferably family matched) as well as unrelated NT children will need to be conducted to determine 1) whether FGIDs, particularly functional constipation, are more prevalent in ASD children than their NT siblings and 2) whether FGIDs, particularly functional constipation, are more prevalent in NT siblings of ASD children than unrelated control children.

It has also been proposed that ASD children with GI symptoms may differ from ASD children without GI symptoms. Adams et al. [15] reported that GI symptoms, which were assessed by the modified 6-GSI score, correlated with ASD severity that was assessed by the parental completed Autism Treatment Evaluation Checklist (ATEC) [15]. However, in our study, the mean ADOS-CSS severity scores, which are based on physician rather than parental evaluations, were not significantly different between the ASD children with at least one FGID and those without any FGID.

To assess the overall emotional and behavioral status of ASD and NT siblings we analyzed parental completed CBCL6/18 normalized T scores. The results of this study demonstrated that ASD children with or without FGID exhibited significantly higher scores than the NT siblings without FGID, consistent with previous studies reporting higher scores in ASD children compared to control children [3437]. However, this study demonstrated that the NT siblings with at least one FGID had normalized CBCL6/18 T-scores that approximated those of ASD children with or without at least one FGID, and were also significantly higher than NT siblings without any FGIDs.

Although this is one of the larger studies comparing ASD, FGID and the gut microbiota in ASD children and NT siblings, our sample size is still modest given the number of variables included. Recruitment of additional ASD children and their NT siblings (preferably family matched), as well as unrelated control children that use the same phenotyping tools will need to be conducted to confirm the findings reported in this study. Furthermore, it may be necessary to possibly restrict dietary intake of plant-derived supplements, such as chia seeds, to determine whether the Cyanobacteria/Chloroplast signal is of biological significance with respect to ASD pathophysiology.

Supporting Information

S1 Fig. Measurement platform comparison for MiSeq V1V2, MiSeq V1V3 and qPCR.

The square root transformed relative abundance are plotted for each pair of platforms with the fitted line in proposed SEM model. For two perfectly consistent platforms, the fitted line is y = x.

https://doi.org/10.1371/journal.pone.0137725.s001

(TIF)

S1 Table. P-values for pairwise comparisons of CBCL scores (See Table 2) between ASD w FGID, ASD w/o FGID, NT sib w FGID, and NT sib w/o FGID.

https://doi.org/10.1371/journal.pone.0137725.s002

(DOCX)

S2 Table. Relative abundances of major phyla for ASD children with FGID, ASD children without FGID, NT siblings with FGID, and NT siblings without FGID.

The mean relative abundance (sequence count for phylum/total sequence count) ± standard deviation for the V1V2 and V1V3 datasets are listed below. The data from family matched and unmatched ASD children and NT siblings are included.

https://doi.org/10.1371/journal.pone.0137725.s003

(DOCX)

S3 Table. Relative abundances of low abundance phyla for ASD w FGID, ASD w/o FGID, NT w FGID, and NT w/o FGID, based on the V1V2 and V1V3 datasets.

The mean relative abundance (sequence count for phylum/total sequence count) ± standard deviation are listed for each phylum. The data from family matched ASD children and NR siblings are included.

https://doi.org/10.1371/journal.pone.0137725.s004

(DOCX)

S4 Table. Relative abundance of selected taxa in ASD children with FGID, ASD w/o FGID, NT siblings with FGID and NT siblings without FGID.

The mean relative abundances (sequence count for genera/total sequence count) ± standard deviations are listed for the V1V2 and the V1V3 datasets. N.D. is not detected based on filtering criteria. The taxa that exhibited significant effects of ASD, FGID and/or ASD*FGID first order interactions are bolded. The data from family matched and unmatched ASD children are included.

https://doi.org/10.1371/journal.pone.0137725.s005

(DOCX)

S5 Table. Measurement platform comparison for MiSeq V1V2, MiSeq V1V3 and qPCR.

The 95% CI of estimated coefficients in the proposed SEM model are listed in a pairwise manner. For two perfectly consistent platforms, the slope is 1 and the intercept is 0.

https://doi.org/10.1371/journal.pone.0137725.s006

(DOCX)

Acknowledgments

We are grateful to all of the families at the participating SFARI Simplex Collection (SSC) sites, the staff at SFARI and the Interactive Autism Network (Kennedy Kreiger Institute) and Maryann Parrish R.N. for assistance with the study recruitment process. We thank Erum Ahmed for assistance in extracting the fecal DNA samples and carrying out the qPCR assays and Christine Collins for assistance in extracting the fecal DNA, carrying out the qPCR assays and entering data for the food diaries.

Author Contributions

Conceived and designed the experiments: JS LJZ GG LL-K KDG DF EL. Performed the experiments: JS LJZ LMR CER DI DF. Analyzed the data: JS LJZ XT YZ WZ LL-K GG CER DF EL. Contributed reagents/materials/analysis tools: JS LJZ LR XT YZ WZ GG LL-K CER DF EL. Wrote the paper: EL JS LJZ LMR LL-K LMR XT YZ WZ GG KDG CR DI DF.

References

  1. 1. Wingate M, Kirby RS, Pettygrove S, Cunniff C, Schulz E, Ghosh T, et al. (2014) Prevalence of autism spectrum disorder among children aged 8 Years—Autism and developmental disabilities monitoring network, 11 Sites, United States, 2010. MMWR Surveill Summ. 63:1–21.
  2. 2. Kim YS, Leventhal BL. (2015) Genetic epidemiology and insights into interactive genetic and environmental effects in autism spectrum disorders. Biol Psychiatry. 77:66–74 pmid:25483344
  3. 3. Colvert E, Tick B, McEwen F, Stewart C, Curran SR, Woodhouse E, et al. (2015) Heritability of autism spectrum disorder in a UK population-based twin sample. JAMA Psychiatry. 72:415–23. pmid:25738232
  4. 4. Frank DN, Zhu W, Sartor RB, Li E. (2011) Investigating the biological and clinical significance of human dysbioses. Trends Microbiol. 19:427–34. pmid:21775143
  5. 5. Krajmalnik-Brown R, Lozupone C, Kang DW, Adams JB. (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 pmid:25769266
  6. 6. Parracho HM, Bingham MO, Gibson GR, McCartney AL. (2005) Differences between the gut microbiota of children with autistic spectrum disorders and that of healthy children. J Med Microbiol. 2005;54:987–91. pmid:16157555
  7. 7. Wang L, Christophersen CT, Sorich MJ, Gerber JP, Angley MT, Conlon MA. (2011) Low relative abundances of the mucolytic bacterium Akkermansia muciniphila and Bifidobacterium spp. in feces of children with autism. Appl Environ Microbiol 77:6718–21. pmid:21784919
  8. 8. Wang L, Christophersen CT, Sorich MJ, Gerber JP, Angley MT, Conlon MA. (2013) Increased abundance of Sutterella spp. and Ruminococcus torques in feces of children with autism spectrum disorder. Mol Autism. 2013;4:42. pmid:24188502
  9. 9. Gondalia SV, Palombo EA, Knowles SR, Cox SB, Meyer D, Austin DW. (2012) Molecular characterization of gastrointestinal microbiota of children with autism (with and without gastrointestinal dysfunction) and their neurotypical siblings. Autism Res 5: 419–27. pmid:22997101
  10. 10. Williams BL, Hornig M, Parekh T, Lipkin WI (2012) 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 3: pii: e00261–11. pmid:22233678
  11. 11. Kang DW, Park JG, Ilhan ZE, Wallstrom G, Labaer J, Adams JB, et al. (2013) Reduced incidence of Prevotella and other fermenters in intestinal microbiota of autistic children. PLoS One 8: e68322. pmid:23844187
  12. 12. 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–87. pmid:25446201
  13. 13. De Angelis M, Piccolo M, Vannini L, Siragusa S, De Giacomo A, Serrazzanetti DI, et al. (2013). Fecal microbiota and metabolome of children with autism and pervasive developmental disorder not otherwise specified. PLoS One, 8(10), e76993. pmid:24130822
  14. 14. Finegold SM, Dowd SE, Gontcharova V, Liu C, Henley KE, Wolcott RD, et al. (2010) Pyrosequencing study of fecal microbiota of autistic and control children. Anaerobe 16: 444–53. pmid:20603222
  15. 15. Adams JB, Johansen LJ, Powell LD, Quig D, Rubin RA. (2011) Gastrointestinal flora and gastrointestinal status in children with autism—comparisons to typical children and correlation with autism severity. BMC Gastroenterol 11: 22. pmid:21410934
  16. 16. Saulnier DM, Riehle K, Mistretta TA, Diaz MA, Mandal D, Raza S, et al. (2011) Gastrointestinal microbiome signatures of pediatric patients with irritable bowel syndrome. Gastroenterology. 141:1782–91. pmid:21741921
  17. 17. Simrén M, Barbara G, Flint HJ, Spiegel BM, Spiller RC, Vanner S, et al. (2013) Intestinal microbiota in functional bowel disorders: a Rome foundation report. Gut 69:159–76
  18. 18. Shankar V, Agans R, Holmes B, Raymer M, Paliy O. (2013) Do gut microbial communities differ in pediatric IBS and health? Gut Microbes. 4:347–52 pmid:23674073
  19. 19. Kim SE, Choi SC, Park KS, Park MI, Shin JE, Lee TH, et al. (2015) Change of fecal flora and effectiveness of the short-term VSL#3 probiotic treatment in patients with functional constipation. J Neurogastroenterol Motil. 2015;21:111–20. pmid:25537674
  20. 20. Ringel Y, Maharshak N. (2013) Intestinal microbiota and immune function in the pathogenesis of irritable bowel syndrome. Am J Physiol Gastrointest Liver Physiol. 305:G529–41. pmid:23886861
  21. 21. Rajilić-Stojanović M, Biagi E, Heilig HG, Kajander K, Kekkonen RA, Tims S, et al. (2011) Global and deep molecular analysis of microbiota signatures in fecal samples from patients with irritable bowel syndrome. Gastroenterology. 141:1792–801. pmid:21820992
  22. 22. Jeffery IB, O'Toole PW, Öhman L, Claesson MJ, Deane J, Quigley EM, et al. (2012) An irritable bowel syndrome subtype defined by species-specific alterations in faecal microbiota. Gut. 61:997–1006. pmid:22180058
  23. 23. Li E, Hamm CM, Gulati AS, Sartor RB, Chen H, Wu X, et al. (2012) Inflammatory bowel diseases phenotype, C. difficile and NOD2 genotype are associated with shifts in human ileum associated microbial composition. PLoS One.7:e26284. pmid:22719818
  24. 24. Krajmalnik-Brown R, Lozupone C, Kang DW, Adams JB. (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. pmid:25769266
  25. 25. Rasquin A, Di Lorenzo C, Forbes D, Guiraldes E, Hyams JS, Staiano A, et al. (2006) Childhood functional gastrointestinal disorders: child/adolescent. Gastroenterology. 130:1527–37. pmid:16678566
  26. 26. Jones MP, Dilley JB, Drossman D, Crowell MD. (2006) Brain-gut connections in functional GI disorders: anatomic and physiologic relationships. Neurogastroenterol Motil. 18:91–103. pmid:16420287
  27. 27. Schneider CK, Melmed RD, Barstow LE, Enriquez FJ, Ranger-Moore J, Ostrem JA: (2006) Oral human immunoglobulin for children with autism and gastrointestinal dysfunction: a prospective, open-label study. J Autism Dev Disord. 36:1053–64. pmid:16845577
  28. 28. McElhanon BO, McCracken C, Karpen S, Sharp WG. (2014) Gastrointestinal symptoms in autism spectrum disorder: a meta-analysis. Pediatrics. 133:872–83. pmid:24777214
  29. 29. Hsiao EY (2014) Gastrointestinal issues in autism spectrum disorder. Harv Rev Psychiatry 22: 104–11. pmid:24614765
  30. 30. Furuta GT, Williams K, Kooros K, Kaul A, Panzer R, Coury DL, et al. (2012) Management of constipation in children and adolescents with autism spectrum disorders. Pediatrics. 130 Suppl 2:S98–105. pmid:23118260
  31. 31. Gorrindo P, Williams KC, Lee EB, Walker LS, McGrew SG, Levitt P. (2012) Gastrointestinal dysfunction in autism: parental report, clinical evaluation, and associated factors. Autism Res. 5:101–8. pmid:22511450
  32. 32. Fischback GD, Lord C. (2010) The Simons Simplex Collection: a resource for identification of autism genetic risk factors. Neuron. 68:192–95. pmid:20955926
  33. 33. Achenbach TM, Rescorla LA. (2001) Manual for the ASEBA school-age forms & profiles: an integrated system of multi-informant assessment. Burlington: University of Vermont, Research Center for Children, Youth & Families.
  34. 34. Duarte CS, Bordin IA, de Oliveira A, Bird H. (2003) The CBCL and the identification of children with autism and related conditions in Brazil: pilot findings. J Autism Dev Disord. 33:703–7. pmid:14714937
  35. 35. Mazefsky CA, Anderson R, Conner CM, Minshew N. (2011) Child behavior checklist scores for school-aged children with autism: preliminary evidence of patterns suggesting the need for referral. J Psychopathol Behav Assess. 33:31–37. pmid:22661827
  36. 36. Pandolfi V, Magyar CI, Dill CA. (2012) An initial psychometric evaluation of the CBCL 6–18 in a sample of youth with autism spectrum disorders. Res Autism Spectr Disord. 6: 96–108. pmid:22059091
  37. 37. Skokauskas N, Gallagher L (2012) Mental health aspects of autistic spectrum disorders in children. J Intellect Disabil Res. 56:248–57. pmid:21554467
  38. 38. Maeda H, Fujimoto C, Haruki Y, Maeda T, Kokequchi S, Petelin M, et al. (2003) Quantitative real-time PCR using TaqMan and SYBR Green for Actinobacillus actinomycetemcomitans, Porphyromonas gingivalis, Prevotella intermedia, tetQ gene and total bacteria. FEMS Immunol Med Microbiol. 39:81–86. pmid:14557000
  39. 39. Frank DN, St Amand AL, Feldman RA, Boedeker EC, Harpaz N, Pace NR. (2007) Molecular-phylogenetic characterization of microbial community imbalances in human inflammatory bowel diseases. Proc Natl Acad Sci U S A. 104:13780–85. pmid:17699621
  40. 40. Larsen N, Vogensen FK, van den Berg FWJ, Nielsen DS, Andreasen AS, Pedersen BK, et al. (2010) Gut microbiota in human adults with type 2 diabetes differs from non-diabetic adults. PLoS One 5: e9085. pmid:20140211
  41. 41. Rinttilä T, Kassinen A, Malinen E, Krogius L, Palva A. (2004) Development of an extensive set of 16S rDNA-targeted primer for quantification of pathogenic and indigenous bacteria in faecal samples by real-time PCR. J. Appl. Microbiol. 97: 1166–77. pmid:15546407
  42. 42. Gulati AS, Shanahan MT, Arthur JC, Grossniklaus E, von Furstenberg RJ, Kreuk L, et al. (2012) Mouse background strain profoundly influences Paneth cell function and intestinal microbial composition. PLoS One. 7:e32403. pmid:22384242
  43. 43. Hara N, Alkanani AK, Ir D, Robertson CE, Wagner BD, Frank DN, et al. (2012) Prevention of virus-induced type 1 diabetes with antibiotic therapy. J Immunol 189: 3805–14. pmid:22988033
  44. 44. Markle JG, Frank DN, Mortin-Toth S, Robertson CE, Feazel LM, Rolle-Kampczyk U, et al. (2013) Sex differences in the gut microbiome drive hormone-dependent regulation of autoimmunity. Science 339: 1084–88. pmid:23328391
  45. 45. Ewing B, Hillier L, Wendl MC, Green P. (1998) Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res 8: 175–85. pmid:9521921
  46. 46. Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R (2011) UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27: 2194–2200. pmid:21700674
  47. 47. Schloss PD, Westcott SL (2011) Assessing and improving methods used in operational taxonomic unit-based approaches for 16S rRNA gene sequence analysis. Appl Environ Microbiol 77: 3219–26. pmid:21421784
  48. 48. Pruesse E, Peplies J, Glockner FO (2012) SINA: accurate high-throughput multiple sequence alignment of ribosomal RNA genes. Bioinformatics 28: 1823–29. pmid:22556368
  49. 49. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P,et al. (2013) The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res 41: D590–96. pmid:23193283
  50. 50. McDonald JH. (2009). Handbook of biological statistics (Vol. 2, pp. 64–68). Baltimore, MD: Sparky House Publishing.
  51. 51. Robertson CD, Harris JK, Wagner BD, Granger D, Browne K, Tatem B, et al. (2013) Explicet: graphical user interface software for metadata driven management, analysis and visualization of microbiome data. Bioinformatics 29:3100–01. pmid:24021386
  52. 52. Frank DN, Bales ES, Monks J, Jackman MJ, MacLean PS, Ir D, et al. (2015) Perilipin-2 modulates lipid absorption and microbiome responses in the mouse intestine. PLoS One 10:e0131944. pmid:26147095
  53. 53. Oksanen J, Kindt R, Legendre P, O'Hara B, Simpson GL, Solymos P, O’Hara RB, et al. (2008) Vegan: Community Ecology Package. R package version 1.15–1. Available: http://cran.r-project.org/, http://vegan.r-forge.r-project.org.
  54. 54. Anderson MJ (2001) A new method for non-parametric multivariate analysis of variance. Austral Ecology 26: 32–46.
  55. 55. Anderson MJ, Crist TO, Chase JM, Vellend M, Inouye BD, Freestone AL, et al. (2011) Navigating the multiple meanings of β diversity: a roadmap for the practicing ecologist. Ecol Lett. 14:19–28. pmid:21070562
  56. 56. Kim RJ, Biswas K, Hoggard M, Taylor MW, Douglas RG. (2015) Paired analysis of the microbiota of surface mucus and whole‐tissue specimens in patients with chronic rhinosinusitis. Int Forum Allergy Rhinol. pmid:26215930
  57. 57. Curran JM (2013). Hotelling: Hotelling's t-squared test and variants. R package version 1.0–2. Available: http://cran.r-project.org/package=hotelling.
  58. 58. McMurdie PJ, Holmes S. (2014) Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput Biol. 10: e1003531. pmid:24699258
  59. 59. Benjamini Y, Hochberg Y. (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B 57: 289–300.
  60. 60. Skaug H, Fournier D, Bolker B, Magnusson A, Nielsen A. (2015). Generalized linear mixed models using AD model builder. R package version 0.8.1.
  61. 61. Bates D, Maechler M, Bolker B, Walker S. (2014). lme4: Linear mixed-effects models using Eigen and S4. R package version 1.1–7, Available: http://CRAN.R-project.org/package=lme4.
  62. 62. Sokal RR, Rohlf FJ. (1994) Biometry: the principles and practice of statistics in biological research. Freeman, New York.
  63. 63. Wu X, Berkow K, Frank DN, Li E, Gulati AS, Zhu W. (2013) Comparative analysis of microbiome measurement platforms using latent variable structural equation modeling. BMC Bioinformatics. 14:79. pmid:23497007
  64. 64. Rosseel Y. (2012) lavaanAn R package for structural equation modeling. JSS 48: 1–36.
  65. 65. Di Rienzi SC, Sharon I, Wrighton KC, Koren O, Hug LA, Thomas BC, et al. (2013) The human gut and groundwater harbor non-photosynthetic bacteria belonging to a new candidate phylum sibling to Cyanobacteria. Elife. 2:e01102. pmid:24137540
  66. 66. Soo RM, Skennerton CT, Sekiguchi Y, Imelfort M, Paech SJ, Dennis PG, et al. (2014) An expanded genomic representation of the phylum cyanobacteria. Genome Biol Evol. 6:1031–45. pmid:24709563
  67. 67. Available: http://freedom-articles.toolsforfreedom.com/autism-cure-omega-3-chia-hemp/.
  68. 68. Available: http://www.newautism.com/ch-ch-ch-chia-how-chia-seeds-may-help-with-autism/2021/.
  69. 69. Davis JJ, Xia F, Overbeek RA, Olsen GJ. (2013) Genomes of the class Erysipelotrichia clarify the firmicute origin of the class Mollicutes. Int J Syst Evol Microbiol. 63:2727–41. pmid:23606477
  70. 70. Lai Q, Liu Y, Yuan J, Du J, Wang L, Sun F, Shao Z. (2014) Multilocus sequence analysis for assessment of phylogenetic diversity and biogeography in Thalassospira bacteria from diverse marine environments. PLoS One. 9:e106353 pmid:25198177
  71. 71. Lamont RF, Sobel JD, Akins RA, Hassan SS, Chaiworapongsa T, Kusanovic JP, et al. (2011) The vaginal microbiome: new information about genital tract flora using molecular based techniques. BJOG. 118:533–49. pmid:21251190
  72. 72. Wade M, Prime H, Madigan S. (2015) Using sibling designs to understand neurodevelopmental disorders: from genes and environments to prevention programming. Biomed Res Int. 2015:672784. pmid:26258141
  73. 73. Miele E, Simeone D, Marino A, Greco L, Auricchio R, Novek SJ, et al. (2004) Functional gastrointestinal disorders in children: an Italian prospective survey. Pediatrics 114:73–8. pmid:15231910
  74. 74. Karabulut GS, Beşer OF, Erginöz E, Kutlu T, Cokuğraş FÇ, Erkan T. (2013) The incidence of irritable bowel syndrome in children using the Rome III criteria and the effect of trimebutine treatment. J Neurogastroenterol Motil. 19:90–3. pmid:23350053