Microbiota Transfer Therapy alters gut ecosystem and improves gastrointestinal and autism symptoms: an open-label study
- Dae-Wook Kang†1,
- James B. Adams†2,
- Ann C. Gregory†3, 15,
- Thomas Borody4,
- Lauren Chittick5, 15,
- Alessio Fasano6,
- Alexander Khoruts7, 8, 9,
- Elizabeth Geis2,
- Juan Maldonado1,
- Sharon McDonough-Means10,
- Elena L. Pollard2,
- Simon Roux5, 15,
- Michael J. Sadowsky8, 11,
- Karen Schwarzberg Lipson12,
- Matthew B. Sullivan3, 5, 15, 16Email author,
- J. Gregory Caporaso12, 13Email author and
- Rosa Krajmalnik-Brown1, 14Email authorView ORCID ID profile
© The Author(s). 2017
Received: 7 October 2016
Accepted: 21 December 2016
Published: 23 January 2017
Autism spectrum disorders (ASD) are complex neurobiological disorders that impair social interactions and communication and lead to restricted, repetitive, and stereotyped patterns of behavior, interests, and activities. The causes of these disorders remain poorly understood, but gut microbiota, the 1013 bacteria in the human intestines, have been implicated because children with ASD often suffer gastrointestinal (GI) problems that correlate with ASD severity. Several previous studies have reported abnormal gut bacteria in children with ASD. The gut microbiome-ASD connection has been tested in a mouse model of ASD, where the microbiome was mechanistically linked to abnormal metabolites and behavior. Similarly, a study of children with ASD found that oral non-absorbable antibiotic treatment improved GI and ASD symptoms, albeit temporarily. Here, a small open-label clinical trial evaluated the impact of Microbiota Transfer Therapy (MTT) on gut microbiota composition and GI and ASD symptoms of 18 ASD-diagnosed children.
MTT involved a 2-week antibiotic treatment, a bowel cleanse, and then an extended fecal microbiota transplant (FMT) using a high initial dose followed by daily and lower maintenance doses for 7–8 weeks. The Gastrointestinal Symptom Rating Scale revealed an approximately 80% reduction of GI symptoms at the end of treatment, including significant improvements in symptoms of constipation, diarrhea, indigestion, and abdominal pain. Improvements persisted 8 weeks after treatment. Similarly, clinical assessments showed that behavioral ASD symptoms improved significantly and remained improved 8 weeks after treatment ended. Bacterial and phagedeep sequencing analyses revealed successful partial engraftment of donor microbiota and beneficial changes in the gut environment. Specifically, overall bacterial diversity and the abundance of Bifidobacterium, Prevotella, and Desulfovibrio among other taxa increased following MTT, and these changes persisted after treatment stopped (followed for 8 weeks).
This exploratory, extended-duration treatment protocol thus appears to be a promising approach to alter the gut microbiome and virome and improve GI and behavioral symptoms of ASD. Improvements in GI symptoms, ASD symptoms, and the microbiome all persisted for at least 8 weeks after treatment ended, suggesting a long-term impact.
This trial was registered on the ClinicalTrials.gov, with the registration number NCT02504554
KeywordsAutism spectrum disorders (ASD) Fecal microbiota transplant (FMT) Clinical trial Gut bacteria Gut bacteriophage Microbiome Virome
Autism spectrum disorders (ASDs) are complex neurobiological disorders that impair social interactions and communication and lead to restricted, repetitive, and stereotyped patterns of behavior, interests, and activities . While ASD diagnoses are increasing, with ~1–2% of children currently diagnosed worldwide , the causes of this disorder remain poorly understood and appear to involve a complex interplay of genetic and environmental factors, of which the microbiome is an environmental factor that is partially inherited from the mother . Despite increased ASD diagnoses, there remains no US Food and Drug Administration (FDA)-approved pharmaceutical treatment to alleviate core ASD symptoms . Coincident with ASD, many children and adults also experience significant gastrointestinal (GI) symptoms, such as constipation, diarrhea, and alternating constipation/diarrhea , which correlate with ASD severity [6, 7]. Such GI symptoms appear to be due, in part, to dysbiotic gut microbiota  and perhaps their missing roles on modulating metabolites (e.g., 4-ethylphenylsulfate, indolepyruvate, and corticosterone) that affect GI function and neurobiological conditions, such as ASD and anxiety [9, 10]. Many children with ASD often undergo increased oral antibiotic treatment during the first 3 years of life , which is thought to destabilize their gut microbiota  and open opportunities for competitive potential pathogens to contribute to ASD severity [13, 14]. A number of studies reported that children with ASD have altered gut bacteria profiles compared with neurotypcial children [13–18], although in certain cohorts, no significant difference has been reported [19, 20]. Because children with ASD have lower abundances of fermentative bacteria (e.g., Prevotella copri), and lower overall bacterial diversity, it has also been hypothesized that lack of beneficial gut microbiota impairs neurological health . Consistent with this, experiments done in an ASD mouse model demonstrated that augmentation with Bacteroides fragilis alone could alter gut microbiota and blood metabolite profiles, correct increased gut permeability (gaps in cell-to-cell junctions), and improve ASD-associated behaviors . In children with ASD, a small open-label study found that 8 weeks of treatment with oral vancomycin (a non-absorbable antibiotic which acts only in the gut) led to major improvements in both GI symptoms and ASD symptoms, although the benefits were lost within a few weeks after treatment was stopped . Thus, gut microbiota appears strongly associated with ASD. Viruses are also abundant in the gut  and may also impact ASD symptoms by modulating the abundance, evolutionary trajectories, and metabolic outputs of gut microbiota like they do in other environments .
Interest in rebalancing human gut microbiota to treat disease is growing . Diet, antibiotics, probiotics, prebiotics, and fecal microbiota transplants are treatments with reported potential [26–30]. For ASD, however, only temporary symptom improvements have been reported from vancomycin treatment , and probiotics have had mixed clinical results with minimal microbiota analysis or long-term follow-up . Contrasting to probiotics which contain a few bacterial species from milk cultures, fecal microbiota transplant (FMT) contains approximately a thousand bacterial species native to the gut and has helped treat recurrent Clostridium difficile infection  and is promising for the treatment of chronic inflammatory diseases such as inflammatory bowel disease  and insulin sensitivity . Therefore, ASD’s GI and behavioral symptoms may derive, at least in part, from gut microbiota dysbiosis and FMT may effectively rebalance the gut microbiota and alleviate some GI and ASD symptoms.
FMT therapy usually involves only a single dose for recurrent C. difficile infection  and other GI conditions, although there is a growing interest in the use of several doses . For this study, a prolonged, daily treatment regimen was implemented based on the clinical experiences of team member Thomas Borody who found that only C. difficile infection is responsive to one or two FMT infusions. All other GI problems—originally described in ulcerative colitis —require multiple infusions of donor microbiota to achieve measurable and long-lasting benefits, including those associated with ASD. An open-label trial was designed to investigate the safety, tolerability, and efficacy of FMT for GI and behavior symptoms in children with ASD. Long-term FMT treatment was administered to 18 children with GI problems and ASD. Clinical responses, gut bacteria, and phage double-stranded DNA profiles were monitored for 18 weeks. Briefly, a modified FMT protocol, termed Microbiota Transfer Therapy (MTT), involved 14 days of oral vancomycin treatment followed by 12–24 h fasting with bowel cleansing, then repopulating gut microbiota by administering a high initial dose of Standardized Human Gut Microbiota (SHGM)  either orally or rectally followed by daily, lower maintenance oral doses with a stomach acid suppressant for 7–8 weeks. A stomach-acid suppressant was used to increase the survival of SHGM through the stomach. Participants were followed for an additional 8 weeks after treatment ended, to determine if treatment effects were temporary or long-lasting. This report focuses on the safety and tolerability of MTT and its effects on microbiota, GI symptoms, and other ASD-related symptoms.
The goals of the study were to follow gut microbiota in healthy and treated children with ASD longitudinally as well as to evaluate an investigational new treatment, MTT, for its effectiveness in children with ASD in treating both GI symptoms (primary outcome) and ASD-related symptoms (secondary outcomes), and to determine the effect of MTT on the gut microbiome.
The study physician first assessed inclusion-exclusion criteria through an extensive review of the participants’ last 2 years of medical records and height/weight/growth charts. Once qualified, autism spectrum diagnosis was verified using the ADI-R, which involved a phone interview of the parents by an ADI-R evaluator. Once qualified and enrolled, participants engaged in an initial 30-min meeting which included a general physical health examination by the study physician and discussion with a project staff member. Participant exclusion criteria included antibiotics use in the prior 6 months or probiotics use in the prior 3 months; dependence on tube feeding; severe GI problems that require immediate treatment (life-threatening); recent/scheduled surgeries; diagnosed as severely malnourished or underweight; and diagnosed with a single-gene disorder, major brain malformations, ulcerative colitis, Crohn’s disease, celiac disease, or eosinophilic esophagitis. None of the neurotypical children had been diagnosed with mental disorders including ASD, attention-deficit hyperactivity disorder (ADHD), depression, or anxiety. None of the neurotypical children had first-degree relatives (i.e., parents and siblings) with ASD. From participants, initial blood and stool samples were collected. Parents were asked to complete a 1-week diet assessment on behalf of their child at the beginning of the study. Participants were recruited primarily from the greater Phoenix, Arizona area; three were from outside that area. Neurotypical families were recruited from friends of the ASD families and professionals who work with ASD families.
The MTT treatment protocol consisted of four key parts: (1) oral vancomycin, (2) MoviPrep, (3) SHGM, and (4) Prilosec. As summarized in Fig. 1, the treatment began with 14 days of oral vancomycin, a non-absorbable broad spectrum antibiotic that stays in the GI tract. A 14-day course of vancomycin was used to ensure that pathogenic bacteria were profoundly suppressed. Prilosec (an acid pump inhibitor) was administered starting on the 12th day of vancomycin, and continued until the end of the lower dosage of SHGM in order to reduce stomach acidity and increase the survival rate of SHGM through the stomach. On day 15, parents administered MoviPrep, a drink that flushes the bowels, to remove most remaining gut bacteria and vancomycin. To enhance its effectiveness, a fasting period of 1 day was implemented during which participants were only allowed to consume clear liquids (children under 12 years were allowed a light breakfast), and then at 4 pm and 8 pm, parents administered the two doses of MoviPrep. On day 16, the participants began either oral administration of SHGM (2.5 × 1012 cells/day) mixed in a chocolate milk, milk substitute, or juice for 2 days (divided into three daily doses), or a single rectal dose of SHGM (2.5 × 1012 cells), given similar to an enema. The rectal dose was administered slowly over 1 h, and participants remained prone for at least several hours, and delayed defecation for at least several hours. The rectal dose was administered under the direct supervision of the study physician, and the first oral dose was similarly administered in the presence of the physician. Participants were randomly assigned to either the oral or rectal route of administration. If one administration route was not tolerated, or if the family preferred the other route, then participants had the option of trying the other route. For participants who received the major initial rectal dose, they waited for 1 week (so the effect of the rectal dose could be evaluated by itself) and then received a lower oral maintenance dose (2.5 × 109 cells) for 7 weeks. In contrast, for participants who received major initial oral doses, they received a lower oral maintenance dose (2.5 × 109 cells) for 8 weeks, directly after the major initial oral dose. The lower maintenance SHGM doses were self-administered orally every day up to the end of week 10. After treatment was stopped, participants were monitored for another 8 weeks.
Standardized human gut microbiota
Instead of pure stool, this study involved the use of standardized human gut microbiota that is > 99% bacteria and prepared as previously described using stool from healthy individuals as starting material . Briefly, donors underwent rigorous screening that involved regular questionnaires, review of medical history, and physical examinations to rule out infectious disease, metabolic syndrome, gastrointestinal disorders, and neurologic or neurodevelopmental problems. Serologic testing was performed to rule out infection with human immunodeficiency virus-1 and -2; hepatitis A, B, and C; and syphilis. The stool used in preparation was tested for potential bacterial pathogens (C. difficile toxin B, Campylobacter, Salmonella, toxin-producing Escherichia coli, Vibrio, Yersinia, Listeria, methicillin-resistant Staphylococcus aureus, and vancomycin-resistant Enterococcus), potential parasites (Giardia, Cryptosporidium, Cyclospora, and Isospora), and potential viral infections (Rotavirus A, Adenovirus, and Norovirus). Metabolic health of donor individuals was assessed with physical examinations and serologic testing (fasting glucose, lipid panel, liver function tests, and high sensitivity C-reactive protein). In addition, the fluorescent antinuclear antibody was employed as a screen for autoimmunity risk. Any single abnormality resulted in disqualification of the donor and prevents material release. The donated material was then extensively filtered and standardized under anaerobic conditions, following FDA good manufacturing processes (GMP), resulting in > 99% microbiota. The final product was in liquid form which can be frozen and was proven to be highly effective for treating C. difficile . The SHGM was stored in −80 °C freezers at Arizona State University (ASU), and then delivered to families on dry ice every week during the study. Families were instructed to keep the SHGM in a container with dry ice and thaw it shortly before use.
Participants received two different doses of SHGM; the high major dose and a lower maintenance dose. The high-dose SHGM was at a daily dosage of 2.5 × 1012 cells, with 2 days for oral and 1 day for rectal administration. The rationale for the high dose was that after the MoviPrep, 1-day fast is presumably the most critical time in which to provide new beneficial bacteria. The maintenance dose of SHGM for the following 7–8 weeks was 2.5 × 109 cells/day.
Evaluation and sample collection
Parents were asked to collect stool samples from their child on approximately 0, 21, 70, and 126 days and to collect fecal swabs bi-weekly on 0, 14, 21, 28, 42, 56, 70, 84, 98, 112, and 126 days. The stool samples were analyzed to determine the types and amounts of gut microbiota present. For safety tests, blood samples were collected on approximately 0, 19, 33, and 74 days. During the study, the participants met with the physician for an initial physical evaluation (including review of medical history) and following evaluations on 16, 30, and 74 days. The physician had a phone consult with families on 7, 21, 42, and 130 days, and more frequently if adverse symptoms occurred, or if families had any questions. Neurotypical participants did not receive any treatment. They simply provided stool samples (at weeks 0 and 19) and swab samples every 2 weeks for 4 months.
Assessments of gastrointestinal symptoms
Parents/guardians were asked to fill in the Gastrointestinal Symptom Rating Scale (GSRS) and the daily stool records (DSR). The GSRS is an assessment of GI symptoms during the previous week, based on 15 questions, which are then scored in five domains: abdominal pain, reflux, indigestion, diarrhea, and constipation. A score for each domain was reported based on the average within the questions in that domain. The original GSRS used a 4-point scale, but this study employed a revised version which included 7-point Likert scale which also has simpler language . The GSRS were assessed on 0, 7, 14, 21, 28, 35, 42, 56, 74, and 130 days, and the children with ASD were defined as non-responders when they achieved less than 50% reduction in the average GSRS. The baseline DSR was collected daily, for 2 weeks, during the treatment phase, and the last 2 weeks of the observation period. The DSR primarily included a rating of the stool using the Bristol Stool Form scale (1 = very hard, 7 = liquid).
Assessments of autism and related symptoms
The ADI-R is a 2-h structured interview and is one of the primary tools used for clinical diagnosis of autism and autism spectrum disorders. It is not designed to be a measure of autism severity but higher scores are generally consistent with more severe symptoms . The ADI-R was used to verify the diagnosis of ASD for admission into the study. The Parent Global Impressions-III (PGI-III) was introduced here as an expanded version of PGI-R  by using a 7-point scale ranging from “much worse” to “much better.” An “average change” is calculated by computing the average in all 18 scores of the PGI-III-final. This tool was chosen because it was found to be more reliable to ask parents directly about observed changes than to have them estimate symptom severity at beginning and end and then compute a difference . Also, the use of a 7-point scale to detect changes seems to yield a high sensitivity to changes. The Childhood Autism Rating Scale (CARS) is a 15-item scale that can be used to both diagnose autism and ASD and assess the overall severity of the symptoms. The Aberrant Behavior Checklist (ABC) assesses problem behaviors in five areas common in children with ASD, including irritability, lethargy, stereotypy, hyperactivity, and inappropriate speech. The Social Responsiveness Scale (SRS) is a 65-item scale that assesses social impairments, a core issue in autism, including social awareness, social information processing, capacity for reciprocal social communication, social anxiety/avoidance, and autistic preoccupations and traits. The Vineland Adaptive Behavior Scale II (VABS-II) is a measure of the functioning level in four different domains: communication, daily living skills, socialization, and motor skills, and 11 sub-domains. The raw scores were converted into an age equivalent score. Its assessment of adaptive skills complements the ABC, which assesses problem behaviors.
PGI-III on 0, 7, 14, 21, 28, 35, 42, 56, 74, and 130 days and the CARS, ABC, and SRS at baseline, at the end of treatment, and at the end of the observation period were assessed, whereas the VABS-II was assessed at baseline and at the end of the observation period only, because it is lengthy and likely less sensitive to short time periods since it assesses changes in specific adaptive skills. The same professional evaluator assessed the ADI-R and the CARS, and parents assessed the PGI-III, ABC, SRS, and VABS-II.
Microbial DNA extraction and next-generation sequencing
Microbial DNA was extracted from feces, swabs, and donor samples using the PowerSoil® DNA Isolation Kit (Mobio Carlsbard, CA). A 16S rRNA library for MiSeq Illumina platform was prepared according to the protocol from Earth Microbiome Project (http://www.earthmicrobiome.org/emp-standard-protocols/). The barcoded primer set 515f-806r were used for pair-ended sequencing to target the 16S rRNA V4 region . Library preparation and sequencing work were performed at the Microbiome Analysis Laboratory in the Swette Center for Environmental Biotechnology (http://krajmalnik.environmentalbiotechnology.org/microbiome-lab.html). These primers amplify both bacterial and archaeal 16S rRNA genes. Archaea-specific changes were not observed and are not discussed in this manuscript.
Microbiome sequencing data were analyzed using Quantitative Insights Into Microbial Ecology (QIIME) 1.9.1 , biom-format version 2.1.5 , VSEARCH version 1.7.0 (https://github.com/torognes/vsearch), SSU-ALIGN 0.1 , and FastTree , as well as custom analytic software (source code at https://github.com/caporaso-lab/autism-fmt1) being prepared for release in QIIME 2. Sequence quality control and demultiplexing using QIIME’s split_libraries_fastq.py with default parameters was performed as described in Bokulich et al.  on a per-run basis. The sequences were combined across runs by merging the resulting files using the cat Unix command, and sequences were clustered into operational taxonomic units (OTUs) at sequence similarities of 100 and 97%. One-hundred percent OTUs were computed using a pipeline designed for this study. First, sequences were clustered into 100% OTUs with VSEARCH, and the resulting data were loaded into a BIOM table using the biom from-uc command. OTUs that occurred in only one sample were filtered from the table for computational efficiency. OTU representative sequences were aligned with ssu-align, and high entropy positions were filtered with ssu-mask. A phylogenetic tree of representative sequences was built using FastTreeMP for use in phylogenetic diversity analyses, and representative sequences were taxonomically annotated using QIIME’s RDP Classifier wrapper against the Greengenes 13_5 reference database. After filtering OTUs that were observed in only a single sample, a median of 28,486 sequences per sample was observed. Alpha and beta diversity analyses were performed using QIIME’s core_diversity_analyses.py, at rarefaction depths of 5721 (to retain as many samples as possible) and 10,000 to confirm that the results were similar with more sequences per sample. In a parallel analysis, OTUs were clustered at 97% similarity using QIIME’s pick_open_reference_otus.py with the Greengenes 13_5 reference database and default parameters. Engraftment analyses were performed by using custom software that is provided in the GitHub repository referenced above. Statistics were performed using scipy 0.17.0, visualizations were created with seaborn 0.6.0, and all analyses were performed using Project Jupyter (notebook version 4.0.6).
Isolation and sequencing of viral DNA
Viral DNA was isolated from stool samples as previously described by Minot et al.  with slight modifications. Briefly, 0.5 g of stool was resuspended into 40 mL of SM buffer, spun down at 4000 rpm for 30 min, and filtered the supernatant at 0.2 μm. The filtrate was ultra-centrifuged through a CsCl step gradient as detailed in Thurber et al. . To target dsDNA bacteriophages, the 1.35–1.5 g/mL fraction was collected from the CsCl column and was treated with chloroform and then with DNase I (100 U/mL) followed by the addition of 0.1 M EDTA and 0.1 M EGTA to halt enzyme activity as described . Viral DNA was then extracted using the DNeasy Blood & Tissue Kit. Following DNA extraction, the sequencing libraries were prepared using the NexteraXT kit with two minor changes. During the library preparation, input DNA was PCR amplified with 18–25 cycles. When input DNA concentrations were low, the buffer ATM was added at a 1:10 dilution. Sequencing was carried out on a MiSeq v3 2 × 300 at one sixth of a lane per sequencing library.
The quality control was performed on sequence reads using Trimmomatic  to remove adaptors, trim low-quality ends of reads (reads were cut as soon as the base quality dropped below 20 on a 4 bp window), and discard short reads (< 50 bp). Then, the reads were assembled from each sample using Idba_ud  with kmer size varying from 20 to 100 by increment of 10. The assembled contigs were screened with VirSorter  to identify and remove all microbial genomes sequences (i.e., all contigs >10 kb and not detected as viral by VirSorter in “virome decontamination” mode). Then, a non-redundant dataset of viral contigs was generated by clustering all viral contigs with Cd-hit  using the thresholds previously established (95% ANI on 80% of the shortest sequence) [54, 55]. This resulted in 4759 non-redundant viral sequences longer than 10 kb.
Analyses of viral populations
To determine the viral population relative abundances in the initial samples, the QC reads were mapped back to this non-redundant contigs database with bowtie2 (option—non-deterministic and non-sensitive, default otherwise) . A contig was considered as detected in a sample if covered by reads on more than 75% of its length, and its abundance was computed as the contig average coverage (number of base pairs mapped to the contig divided by contig length) normalized by the total number of base pairs sequenced in the metagenome . The diversity indices, Shannon’s H′ and Peilou’s J, and Bray-Curtis distances were calculated by using the vegan package  in R version 3.2.3 . Bray-Curtis distances were statistically ordinated using the nonmetric multidimensional scaling (NMDS) and then evaluated the influence of the metadata on sample ordination using the “envfit” function with a total of 9999 permutations in the vegan package. Engraftment analyses were performed by using custom Perl scripts. The scripts can be found in the project’s GitHub repository. Viral genes for each viral population were predicted using Prodigal (https://github.com/hyattpd/prodigal/releases/). A blastx for all identified viral genes was performed against the Viral Protein RefSeq to obtain the top three hits with a bitscore of > 50. The familial taxonomy was then obtained for the three hits for each protein. If more than two of the hits had the same familial taxonomy, the viral protein was then assigned that taxonomy. To assign viral taxonomy to the whole viral contig, > 50% of the genes within the contig had to have the same familial taxonomy. To determine if a viral population was similar to the core viral dataset in Manrique et al. , the core contigs genomes were obtained from Manrique et al. and used as a blast database. A blastn of the 1651 viral populations in this dataset was performed against the core 23 phage contigs. If a populations had a blastn alignment length of > 500 bp to one the 23 core gut phage contigs at a percent identify greater than 75%, it was considered related to the core 23 phage contigs.
All commands that were applied for the microbiome analyses are provided in the GitHub repository available at http://github.com/caporaso-lab/autism-fmt1 to facilitate reproducibility of these bioinformatics methods.
Statistical analysis was not utilized to predetermine sample size, since the effect size was unknown. Instead, the study was designed based on our previous research in which statistically significant differences within a similar sample size were detected . The previous study was a case-control comparison that did not include an intervention, and so similar or larger differences were assumed to appear as a result of treatment. Since the sample size is still relatively small, and the data are assumed as non-normally distributed, nonparametric analyses were performed, including the Mann-Whitney U test, Wilcoxon signed-rank test, and Spearman’s correlation test. All p values reported in the study were from two-tailed tests, except the hypothesis on low fiber consumption and low microbial diversity in children with ASD at baseline. p values lower than 0.05 were accepted as significant in clinical data analysis. All p values for bacterial microbiome analyses were corrected using the Benjamini-Hochberg false discovery rate correction, and the resulting corrected values were referred to as q values. q values lower than 0.05 were accepted as significant. For some previously hypothesized beneficial bacteria (Bifidobacterium and Prevotella), q values were not significant, but they were considered to be suggestive of statistical significance (q values less than 0.1 but greater than 0.05). Statistical significance of variance is reported as indicated per experiment in figure legends. All center values in the box plots are median. The top and bottom edges of the box are of the 75th and 25th percentiles of the sample. p values for the phageome analyses are permutation p values calculated from 9999 randomized permutations, with p values lower than 0.05 accepted as significant.
Results and discussion
Eighteen children with ASD each from a different family and 20 neurotypical children from 13 families (6 families had 1 neurotypical participant and 7 families had 2 neurotypical participants) were enrolled in the study reported here. All ASD participants completed the 18-week treatment study (neurotypical children were not treated). Neurotypical children had no first-degree relatives of individuals with ASD. Participants in both groups were of similar age, gender distribution, and body mass index (BMI), but the ASD group had more individuals that were delivered by C-section, used non-standard formula during infancy, and had food allergies and eczema (Additional file 1: Table S1). Children with ASD had marginally lower fiber consumption (one-tailed Mann-Whitney U test, p = 0.07), and their mothers also had significantly lower fiber consumption compared with mothers of neurotypical children (two-tailed Mann-Whitney U test, p < 0.01). Children with ASD were breastfed significantly shorter time than neurotypical children (two-tailed Mann-Whitney U test, p < 0.05). Consumptions on carbohydrate, fat, protein, and calorie were comparable between children with ASD and neurotypical children (Additional file 1: Table S1). Other larger studies reported that children with ASD had more antibiotics administered during the first few years of life , but this ASD group reported a comparable number of antibiotic administrations to the control group during the first 4 years of life (Additional file 1: Table S1). Children with ASD who had moderate or severe GI problems were recruited, which reflected higher GSRS scores in the ASD group than the control group. A summary on participants’ characteristics and their medical and diet history is listed in Additional file 2: Dataset S1.
GI and ASD evaluations
Beyond these GI improvements, ASD-related behavior also improved following MTT. The PGI-II assessment, which evaluates 17 ASD-related symptoms, revealed significant improvement during treatment and no reversion 8 weeks after treatment ended (Fig. 2b). Further, a significant negative correlation between change in GSRS and PGI-III (Spearman’s correlation test showed r = −0.59 and p < 0.001, Additional file 3: Figure S2) suggests that GI symptoms worsen directly with ASD behaviors, and that these can be altered via MTT. The scores on CARS, which rates core ASD symptoms, decreased by 22% from beginning to end of the treatment and 24% (relative to baseline) after 8 weeks of no treatment (Wilcoxon signed-rank test, p < 0.001, Fig. 2c). Children with ASD saw improvement in their scores in the SRS, which assesses social skill deficits, and the ABC, which evaluates irritability, hyperactivity, lethargy, stereotypy, and aberrant speech (Fig. 2d, e). The VABS-II scoring, which evaluates adaptive behaviors such as communication, daily living skills, and socialization, found that the average developmental age increased by 1.4 years (p < 0.001) and across all sub-domain areas (Additional file 3: Figure S3) during MTT, though the final VABS-II age equivalent was still lower than their chronological age. Finally, MTT appears to be beneficial for children ages 7–16 years old (no significant correlations between age and GSRS or CARS improvement), and there was no significant difference in clinical outcomes between those who received the initial SHGM dose orally or rectally.
The MTT treatments were generally well-tolerated, with only temporary adverse effects (primarily mild to moderate hyperactivity and tantrums/aggression) at the beginning of vancomycin treatment, no major changes in blood chemistry or long-term adverse effects were noted. Detailed information is provided in Additional file 4. The improvements in GI and ASD symptoms are consistent with a previous 8-week trial of the use of vancomycin for treating children with ASD , but a key difference is that in the previous study, benefits were lost within a few weeks of stopping vancomycin therapy (despite the use of standard probiotics in some children), whereas in this study, the benefits continued for at least 8 weeks. It is also relevant to note that GI and ASD symptoms slowly improved over the 10-week MTT treatment and 8-week observation period, since this observation is very different from FMT treatment for C. difficile, where a single dose generally leads to recovery within a few days . Thus, it appears likely that extended treatment with FMT over many weeks, as done in this study, is necessary to observe these benefits.
Bacterial changes after MTT
Importantly, the donor bacterial community was at least partially engrafted in the recipient gut, consistent with earlier work  and a recent study of the efficacy of FMT . Specifically, the unweighted UniFrac distance (i.e., a qualitative measure of the dissimilarity of a pair of microbial communities based on shared OTUs) between the host gut and their most recent donor sample significantly decreased over time (Fig. 3c; two-tailed Mann-Whitney U test p < 0.01 at 3 weeks and p < 0.001 at 10 and 18 weeks) and remained more similar to the donor’s bacterial community 8 weeks after treatment stopped. By the end of the treatment (week 10) and 8 weeks after the treatment stopped (week 18), the distance between the recipient and the donor bacterial community was less than normal interpersonal bacterial community variation (in this case, defined by variation between the neurotypical controls) (Fig. 3c). These signatures of engraftment suggest that MTT overcame “colonization resistance” .
Specific genera that significantly changed in their relative abundances with treatment included Bifidobacterium, Prevotella, and Desulfovibrio (Fig. 3e–g and Additional file 5: Dataset S2). Bifidobacterium was reported to be underrepresented in children with ASD [7, 21, 66], also observed in this study at baseline (two-tailed Mann-Whitney U test p < 0.05), but following MTT, the relative abundance of Bifidobacterium significantly increased fourfold and became comparable to its relative abundance in neurotypical children (Fig. 3e). This suggests strong engraftment by these microbes in particular. Additionally, relative abundances of Prevotella and Desulfovibrio significantly increased after MTT from baseline to 8 weeks following treatment (Fig. 3f, g). Initially, the relative abundance of Prevotella was comparable between neurotypical children and children with ASD at baseline, which was not consistent with our previous cohort study with 20 neurotypical children and 19 children with ASD . However, the increase in the relative abundance of Prevotella after MTT is consistent with their potentially beneficial role in the gut of children with ASD. The increased relative abundance of Desulfovibrio is intriguing, since their role in the human gut has been controversially proposed as either commensal  or detrimental [18, 67]. Both Prevotella and Desulfovibrio were on average more abundant in MTT recipients following treatment than in the donor samples, illustrating that the transferred microbiota changes the gut environment in a way that is more hospitable to recruit new commensal bacteria. Taken together, these data suggest that MTT successfully shifts the ASD bacterial community toward that of age/gender-matched healthy controls and to that of their donors.
Phage community changes after MTT
While the role of phages in the gut is largely understudied in comparison to the role of bacteria, this study and recent research has begun to uncover the potential role of phages in the gut. In healthy individuals, the gut virome is highly stable over time [23, 47, 69], with some phage populations hypothesized to provide a non-host-derived protective barrier to invading bacterial pathogens [70, 71]. While there is high inter-individual variation , recent analyses have identified a distinct subset of phages that are found across the majority of healthy individuals . These “healthy” phage populations represent < 5% of the phage population identified in our study (Additional file 3: Figure S5). In individuals with gastrointestinal disease (i.e., ulcerative colitis and Crohn’s disease), these phage populations represent a significantly smaller percentage of the gut phage community . This shift in viral community structure is hypothesized to allow potentially harmful bacteria and viruses to proliferate through phage-mediated dysbiosis, whereby perturbations to the healthy gut phage community leads to increased abundances of phages and selected reduction in bacterial species [70–72].
Studies looking at the effects of perturbations in the gut from Crohn’s disease in humans  and from diet in mice  show this response with increased diversity of phage communities paired with decreased diversity of the bacterial community. In this study that looks at reversing the negative responses to gut perturbations caused by ASD, no significant changes were observed in the diversity of the phage community, but an altered phage community paired with an increase in the diversity of the bacterial community. This suggests that MTT may be able to reverse phage-mediated dysbiosis of the ASD gut, though further study is necessary to test this assertion.
Study limitations and recommendations
Although study observations are noteworthy, the current open-label trial is not placebo controlled, blinded, or randomized. Here, we list some limitations and how they should be addressed in a follow-up blinded trial with a placebo control arm. First, this exploratory study only looked at the consequences of the combined treatment of MTT. Follow-up studies are needed to determine whether MTT efficacy in our study results solely from vancomycin, MoviPrep, SHGM, Prilosec, or a combination of these four factors. Sandler et al.  reported a temporal efficacy of vancomycin treatment in GI and ASD symptoms, but this study involved a small number of participants who were younger than 7 years old. Previous studies have shown that vancomycin  and proton pump inhibitors  significantly alter gut microbiota. Further studies are, however, also essential in order to clarify how each factor in MTT contributes on changes in gut microbiota in the context of ASD. Second, in this study, participants had a range of GI issues, including constipation, diarrhea, and alternating diarrhea/constipation. Larger studies in future would allow us to look at those groups separately. Alternatively, a more homogeneous cohort (e.g., children with shared GI issues and ASD etiologies and similar ages) would allow for better disentanglement of the signal from inter-individual variation from FMT efficacy, since inter-individual gut bacterial  and viral  community variation is high. Third, since this was an open-label study, the effect on GI and ASD symptoms are likely to be subject to placebo effects and should be cautiously interpreted and viewed as preliminary. Fourth, a clinical trial with extended longer observation period after treatment would help determining long-term safety and possible benefits. Lastly, a larger sample size will be essential to clarify associations with other variables, such as the efficacy of oral versus rectal administration of the SHGM.
Further, in follow-up studies, continued use of GI and behavior assessments to carefully track changes in ASD severities, with some additional modifications, is recommended. In this study, the GSRS, SRS, ABC, PGI-III, and VABS-II assessments are reported by parents/guardians, consulting with subjects verbally if the subjects were adolescents. Previous ASD clinical studies have reported disagreement between parent report and that of a pediatric gastroenterologist in terms of the specific GI symptoms and diagnoses, for some metrics (e.g., GI assessment instrument—QPGS-Rome III ). As a result, clinical expertise, in addition to parent/subject reports, could provide more reliable and independent assessments.
Together, these findings suggest that MTT is safe and well-tolerated in children with ASD ages 7–16 years. MTT led to significant improvements in both GI- and ASD-related symptoms, and the improvements were sustained at least 8 weeks after treatment. Coincident with these clinical improvements, both microbiota and phage from the donors appear to have engrafted, at least partially, in the recipients. This shifted gut microbiota of children with ASD toward that of neurotypical children is consistent with the hypothesis that gut microbiota may be at least partially responsible for GI and ASD symptoms. While this study was an open-label trial that is subject to placebo effects, these results are promising and provide a crucial step for understanding the connection between the microbiome and ASD. A randomized, double-blind, placebo-controlled study is the next step to investigate the value of MTT in treating children with ASD and GI problems.
Aberrant Behavior Checklist
Attention-deficit hyperactivity disorder
Autism Diagnostic Interview-Revised
Autism spectrum disorder
Body mass index
Childhood Autism Rating Scale
Daily stool records
Food and Drug Administration
Fecal microbiota transplant
Good manufacturing processes
Gastrointestinal Symptom Rating Scale
Microbiota transfer therapy
Nonmetric multidimensional scaling
Operational taxonomic units
Parent Global Impressions–III
Quantitative Insights into Microbial Ecology
Standardized human gut microbiota
Social Responsiveness Scale
Vineland Adaptive Behavior Scale II
We gratefully thank all the children with ASD, neurotypical children, and their families for participating in the study. We would like to thank S. Dessoy for help with study coordination; G. Harrington for experimental support; A. Hanratty and J. Patel for sample deliveries; M. Moreno for technical support; and G. Ackermann and C. Staley for help with data deposition. We also thank N. Tkacenko, B. McCall, and J. Huerta for phlebotomy/blood processing; D. Zenner, C. Minkner, and A. Millagracia for data entry; Crestovo for FMT products; LabCorp for doing the ChemPanel and complete blood count tests; and E. Linden and Walgreens Pharmacy for providing the vancomycin, MoviPrep, and Prilosec for the study.
This work was supported mainly by the Arizona Board of Regents. The Autism Research Institute also provided supplemental funding. Partial funding was provided by a Gordon and Betty Moore Foundation grant (GBMF #3790) to MBS.
Availability of data and materials
The datasets supporting the conclusions of this article are available in the open-source microbiome database “Qiita” with the study ID number 10532 ( https://qiita.microbio.me) for 16S rRNA gene sequence reads. The virome sequence reads are available in the public repository “iVirus" at: http://mirrors.iplantcollaborative.org/browse/iplant/home/shared/iVirus/ABOR. The datasets supporting the conclusions of this article are also included within the article as Additional files 2 and 5.
D-WK, JBA, MBS, JGC, and RKB designed the study and wrote funding proposal. TB, AF, AK, SM-M, and MJS contributed to the clinical protocol. D-WK, JBA, EG, JM, ELP, and RKB coordinated the study. AK and MJS produced the FMT material. JBA, SM-M, and RKB supervised clinical trial. JBA, EG, SM-M, and ELP gathered the clinical data. D-WK, ACG, LC, JM, and JGC processed the biosamples and performed the data acquisition. D-WK, JBA, ACG, JM, SR, KSL, MBS, JGC, and RKB performed the analysis and interpreted the data. D-WK, JBA, ACG, MBS, JGC, and RKB wrote the main paper, and TB, LC, AF, AK, JM, SM-M, ELP, SR, MJS, and KSL revised the manuscript. All authors read and approved the final manuscript.
Several authors (JBA, D-WK, RKB, TB, AK, JGC, and MJS) have pending/approved patents related to the use of FMT and/or probiotics for various conditions including autism. MJS, AK, JBA, RKB, and D-WK have received research funding from Crestovo for FMT research. JBA, JGC, RKB, and MJS are part-time consultants for Crestovo. The other authors do not have any financial conflicts of interest.
Consent for publication
All authors have read and approved the paper for submission.
Ethics approval and consent to participate
The protocol was approved by US FDA (Investigational new drug number 15886) and the Institutional Review Board of Arizona State University (ASU IRB Protocol #: 00001053). The study was advertised by e-mail to approximately 2500 ASD families in Arizona, USA, using the contact list of the Autism Society of Greater Phoenix and the Autism/Asperger’s Research Program at Arizona State University. Families with children who met the study inclusion and exclusion criteria had a 1-h individual phone call to discuss the study. After the phone call, families who signed the parent permission and child assent forms were provided with initial questionnaires to complete. They were also sent a letter so that their personal physicians would double-check their medications and prepare for the delivery of the vancomycin, Prilosec, and the fecal transplant. All participants provided informed consent. The trial is registered on the ClinicalTrials.gov (NCT02504554).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- World Health Organization. International statistical classification of diseases and health related problems, Tenth Revision (ICD-10). Geneva: World Health Organization; 2004.
- Lai MC, Lombardo MV, Baron-Cohen S. Autism. Lancet. 2014;383(9920):896–910.View ArticlePubMedGoogle Scholar
- Dominguez-Bello MG, Costello EK, Contreras M, Magris M, Hidalgo G, Fierer N, Knight R. Delivery mode shapes the acquisition and structure of the initial microbiota across multiple body habitats in newborns. Proc Natl Acad Sci U S A. 2010;107(26):11971–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Neul JL, Sahin M. Therapeutic advances in autism and other neurodevelopmental disorders. Neurotherapeutics. 2015;12(3):519–20.View ArticlePubMedPubMed CentralGoogle Scholar
- McElhanon BO, McCracken C, Karpen S, Sharp WG. Gastrointestinal symptoms in autism spectrum disorder: a meta-analysis. Pediatrics. 2014;133(5):872–83.View ArticlePubMedGoogle Scholar
- Chaidez V, Hansen RL, Hertz-Picciotto I. Gastrointestinal problems in children with autism, developmental delays or typical development. J Autism Dev Disord. 2014;44(5):1117–27.View ArticlePubMedPubMed CentralGoogle Scholar
- Adams JB, Johansen LJ, Powell LD, Quig D, Rubin RA. Gastrointestinal flora and gastrointestinal status in children with autism-comparisons to typical children and correlation with autism severity. BMC Gastroenterol. 2011;11:22.View ArticlePubMedPubMed CentralGoogle Scholar
- Krajmalnik-Brown R, Lozupone C, Kang DW, Adams JB. 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. 2015;26:26914.
- Hsiao EY, McBride SW, Hsien S, Sharon G, Hyde ER, McCue T, Codelli JA, Chow J, Reisman SE, Petrosino JF, et al. Microbiota modulate behavioral and physiological abnormalities associated with neurodevelopment disorders. Cell. 2013;155(7):1451–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Bravo JA, Forsythe P, Chew MV, Escaravage E, Savignac HM, Dinan TG, Bienenstock J, Cryan JF. 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. 2011;108(38):16050–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Niehus R, Lord C. Early medical history of children with autism spectrum disorders. J Dev Behav Pediatr. 2006;27(2):S120–7.View ArticlePubMedGoogle Scholar
- Willing BP, Russell SL, Finlay BB. Shifting the balance: antibiotic effects on host-microbiota mutualism. Nat Rev Microbiol. 2011;9(4):233–43.View ArticlePubMedGoogle Scholar
- Finegold SM, Dowd SE, Gontcharova V, Liu C, Henley KE, Wolcott RD, Youn E, Summanen PH, Granpeesheh D, Dixon D, et al. Pyrosequencing study of fecal microflora of autistic and control children. Anaerobe. 2010;16(4):444–53.View ArticlePubMedGoogle Scholar
- Williams BL, Hornig M, Parekh T, Lipkin WI. Application of novel PCR-based methods for detection, quantification, and phylogenetic characterization of Sutterella species in intestinal biopsy samples from children with autism and gastrointestinal disturbances. MBio. 2012;3(1):e00261-11.View ArticlePubMedPubMed CentralGoogle Scholar
- Song YL, Liu CX, Finegold SA. Real-time PCR quantitation of Clostridia in feces of autistic children. Appl Environ Microbiol. 2004;70(11):6459–65.View ArticlePubMedPubMed CentralGoogle Scholar
- De Angelis M, Piccolo M, Vannini L, Siragusa S, De Giacomo A, Serrazzanetti DI, Cristofori F, Guerzoni ME, Gobbetti M, Francavilla R. Fecal microbiota and metabolome of children with autism and pervasive developmental disorder not otherwise specified. Plos One. 2013;8(10):18.View ArticleGoogle Scholar
- Wang L, Christophersen CT, Sorich MJ, Gerber JP, Angley MT, Conlon MA. Increased abundance of Sutterella spp. and Ruminococcus torques in feces of children with autism spectrum disorder. Molecular Autism. 2013;4:42.
- Tomova A, Husarova V, Lakatosova S, Bakos J, Vlkova B, Babinska K, Ostatnikova D. Gastrointestinal microbiota in children with autism in Slovakia. Physiol Behav. 2015;138:179–87.View ArticlePubMedGoogle Scholar
- Gondalia SV, Palombo EA, Knowles SR, Cox SB, Meyer D, Austin DW. Molecular characterisation of gastrointestinal microbiota of children with autism (with and without gastrointestinal dysfunction) and their neurotypical siblings. Autism Res. 2012;5(6):419–27.View ArticlePubMedGoogle Scholar
- Son JS, Zheng LJ, Rowehl LM, Tian XY, Zhang YH, Zhu W, Litcher-Kelly L, Gadow KD, Gathungu G, Robertson CE, et al. Comparison of fecal microbiota in children with autism spectrum disorders and neurotypical siblings in the Simons simplex sollection. Plos One. 2015;10(10):19.Google Scholar
- Kang DW, Park JG, Ilhan ZE, Wallstrom G, Labaer J, Adams JB, Krajmalnik-Brown R. Reduced incidence of Prevotella and other fermenters in intestinal microflora of autistic children. PLoS One. 2013;8(7):e68322.View ArticlePubMedPubMed CentralGoogle Scholar
- Sandler RH, Finegold SM, Bolte ER, Buchanan CP, Maxwell AP, Vaisanen ML, Nelson MN, Wexler HM. Short-term benefit from oral vancomycin treatment of regressive-onset autism. J Child Neurol. 2000;15(7):429–35.View ArticlePubMedGoogle Scholar
- Minot S, Bryson A, Chehoud C, Wu GD, Lewis JD, Bushman FD. Rapid evolution of the human gut virome. Proc Natl Acad Sci U S A. 2013;110(30):12450–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Brum JR, Sullivan MB. Rising to the challenge: accelerated pace of discovery transforms marine virology. Nat Rev Microbiol. 2015;13(3):147–59.View ArticlePubMedGoogle Scholar
- Gilbert JA, Krajmalnik-Brown R, Porazinska DL, Weiss SJ, Knight R. Toward effective probiotics for autism and other neurodevelopmental disorders. Cell. 2013;155(7):1446–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Tillisch K, Labus J, Kilpatrick L, Jiang Z, Stains J, Ebrat B, Guyonnet D, Legrain-Raspaud S, Trotin B, Naliboff B, et al. Consumption of fermented milk product with probiotic modulates brain activity. Gastroenterology. 2013;144(7):1394–U1136.View ArticlePubMedGoogle Scholar
- Youngster I, Sauk J, Pindar C, Wilson RG, Kaplan JL, Smith MB, Alm EJ, Gevers D, Russell GH, Hohmann EL. Fecal microbiota transplant for relapsing Clostridium difficile infection using a frozen inoculum from unrelated donors: a randomized, open-label, controlled pilot study. Clin Infect Dis. 2014;58(11):1515–22.View ArticlePubMedPubMed CentralGoogle Scholar
- Prantera C, Lochs H, Grimaldi M, Danese S, Scribano ML, Gionchetti P, Retic Study Grp R-E. Rifaximin-extended intestinal release induces remission in patients with moderately active Crohn's disease. Gastroenterology. 2012;142(3):473–U124.View ArticlePubMedGoogle Scholar
- Brandt LJ, Aroniadis OC, Mellow M, Kanatzar A, Kelly C, Park T, Stollman N, Rohlke F, Surawicz C. Long-term follow-up of colonoscopic fecal microbiota transplant for recurrent Clostridium difficile infection. Am J Gastroenterol. 2012;107(7):1079–87.View ArticlePubMedGoogle Scholar
- Moayyedi P, Quigley EMM, Lacy BE, Lembo AJ, Saito YA, Schiller LR, Soffer EE, Spiegel BMR, Ford AC. The effect of fiber supplementation on irritable bowel syndrome: a systematic review and meta-analysis. Am J Gastroenterol. 2014;109(9):1367–74.View ArticlePubMedGoogle Scholar
- Parracho HM, Gibson GR, Knott F, Bosscher D, Kleerebezem M, McCartney AL. A double-blind, placebo-controlled, crossover-designed probiotic feeding study in children diagnosed with autistic spectrum disorders. Int J Probiot Prebiot. 2010;5(2):69.Google Scholar
- Bagdasarian N, Rao K, Malani PN. Diagnosis and treatment of Clostridium difficile in adults a systematic review. Jama-Journal of the American Medical Association. 2015;313(4):398–408.View ArticleGoogle Scholar
- Moayyedi P, Surette MG, Kim PT, Libertucci J, Wolfe M, Onischi C, Armstrong D, Marshall JK, Kassam Z, Reinisch W, et al. Fecal microbiota transplantation induces remission in patients with active ulcerative colitis in a randomized controlled trial. Gastroenterology. 2015;149(1):102.View ArticlePubMedGoogle Scholar
- Vrieze A, Van Nood E, Holleman F, Salojarvi J, Kootte RS, Bartelsman J, Dallinga-Thie GM, Ackermans MT, Serlie MJ, Oozeer R, et al. Transfer of intestinal microbiota from lean donors increases insulin sensitivity in individuals with metabolic syndrome. Gastroenterology. 2012;143(4):913.View ArticlePubMedGoogle Scholar
- Kunde S, Pham A, Bonczyk S, Crumb T, Duba M, Conrad H, Cloney D, Kugathasan S. Safety, tolerability, and clinical response after fecal transplantation in children and young adults with ulcerative colitis. J Pediatr Gastroenterol Nutr. 2013;56(6):597–601.View ArticlePubMedGoogle Scholar
- Borody TJ, Warren EF, Leis S, Surace R, Ashman O. Treatment of ulcerative colitis using fecal bacteriotherapy. J Clin Gastroenterol. 2003;37(1):42–7.View ArticlePubMedGoogle Scholar
- Hamilton MJ, Weingarden AR, Sadowsky MJ, Khoruts A. Standardized frozen preparation for transplantation of fecal microbiota for recurrent Clostridium difficile infection. Am J Gastroenterol. 2012;107(5):761–7.View ArticlePubMedGoogle Scholar
- Revicki DA, Wood M, Wiklund I, Crawley J. Reliability and validity of the gastrointestinal symptom rating scale in patients with gastroesophageal reflux disease. Qual Life Res. 1998;7(1):75–83.View ArticlePubMedGoogle Scholar
- Lord C, Rutter M, Lecouteur A. Autism Diagnostic Interview-Revised: a revised version of a diagnostic interview for caregivers of individuals with possible pervasive developmental disorders. J Autism Dev Disord. 1994;24(5):659–85.View ArticlePubMedGoogle Scholar
- Adams JB, Audhya T, McDonough-Means S, Rubin RA, Quig D, Geis E, Gehn E, Loresto M, Mitchell J, Atwood S, et al. Effect of a vitamin/mineral supplement on children and adults with autism. BMC Pediatr. 2011;11:30.View ArticleGoogle Scholar
- Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, Owens SM, Betley J, Fraser L, Bauer M, et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. Isme Journal. 2012;6(8):1621–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335–6.View ArticlePubMedPubMed CentralGoogle Scholar
- McDonald D, Clemente JC, Kuczynski J, Rideout JR, Stombaugh J, Wendel D, Wilke A, Huse S, Hufnagle J, Meyer F, et al. The biological observation matrix (BIOM) format or: how I learned to stop worrying and love the ome-ome. Gigascience. 2012;1:6.View ArticleGoogle Scholar
- Nawrocki EP, Kolbe DL, Eddy SR. Infernal 1.0: inference of RNA alignments. Bioinformatics. 2009;25(10):1335–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Price MN, Dehal PS, Arkin AP. FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009;26(7):1641–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Bokulich NA, Subramanian S, Faith JJ, Gevers D, Gordon JI, Knight R, Mills DA, Caporaso JG. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat Methods. 2013;10(1):57–U11.View ArticlePubMedGoogle Scholar
- Minot S, Sinha R, Chen J, Li H, Keilbaugh SA, Wu GD, Lewis JD, Bushman FD. The human gut virome: inter-individual variation and dynamic response to diet. Genome Res. 2011;21(10):1616–25.View ArticlePubMedPubMed CentralGoogle Scholar
- Thurber RV, Willner-Hall D, Rodriguez-Mueller B, Desnues C, Edwards RA, Angly F, Dinsdale E, Kelly L, Rohwer F. Metagenomic analysis of stressed coral holobionts. Environ Microbiol. 2009;11(8):2148–63.View ArticleGoogle Scholar
- Hurwitz BL, Deng L, Poulos BT, Sullivan MB. Evaluation of methods to concentrate and purify ocean virus communities through comparative, replicated metagenomics. Environ Microbiol. 2013;15(5):1428–40.View ArticlePubMedPubMed CentralGoogle Scholar
- Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Peng Y, Leung HCM, Yiu SM, Chin FYL. IDBA-UD: a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth. Bioinformatics. 2012;28(11):1420–8.View ArticlePubMedGoogle Scholar
- Roux S, Enault F, Hurwitz BL, Sullivan MB. VirSorter: mining viral signal from microbial genomic data. Peerj. 2015;3:e985.
- Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22(13):1658–9.View ArticlePubMedGoogle Scholar
- Brum JR, Ignacio-Espinoza JC, Roux S, Doulcier G, Acinas SG, Alberti A, Chaffron S, Cruaud C, de Vargas C, Gasol JM, et al. Patterns and ecological drivers of ocean viral communities. Science. 2015;348(6237):1261498.
- Deng L, Ignacio-Espinoza JC, Gregory AC, Poulos BT, Weitz JS, Hugenholtz P, Sullivan MB. Viral tagging reveals discrete populations in Synechococcus viral genome sequence space. Nature. 2014;513(7517):242.View ArticlePubMedGoogle Scholar
- Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–U354.View ArticlePubMedPubMed CentralGoogle Scholar
- Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, Minchin PR, O'Hara RB, Simpson GL, Solymos P, Stevens MHH, Szoecs E, Wagner H. Vegan: Community Ecology Package. R package version. 2016;2:4-1. https://CRAN.R-project.org/package=vergan.
- R Core Team R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria; 2015. https://www.R-project.org/.
- Manrique P, Bolduc B, Walk ST, van der Oost J, de Vos WM, Young MJ. Healthy human gut phageome. Proc Natl Acad Sci U S A. 2016;113(37):10400–5.View ArticlePubMedGoogle Scholar
- Aroniadis OC, Brandt LJ. Fecal microbiota transplantation: past, present and future. Curr Opin Gastroenterol. 2013;29(1):79–84.View ArticlePubMedGoogle Scholar
- Faith DP. Conservation evaluation and phylogenetic diversity. Biol Conserv. 1992;61(1):1–10.View ArticleGoogle Scholar
- Qin JJ, Li RQ, Raes J, Arumugam M, Burgdorf KS, Manichanh C, Nielsen T, Pons N, Levenez F, Yamada T, et al. A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010;464(7285):59–U70.View ArticlePubMedPubMed CentralGoogle Scholar
- Grehan MJ, Borody TJ, Leis SM, Campbell J, Mitchell H, Wettstein A. Durable alteration of the colonic microbiota by the administration of donor fecal flora. J Clin Gastroenterol. 2010;44(8):551–61.View ArticlePubMedGoogle Scholar
- Li SS, Zhu A, Benes V, Costea PI, Hercog R, Hildebrand F, Huerta-Cepas J, Nieuwdorp M, Salojarvi J, Voigt AY, et al. Durable coexistence of donor and recipient strains after fecal microbiota transplantation. Science. 2016;352(6285):586–9.View ArticlePubMedGoogle Scholar
- Vollaard EJ, Clasener HAL. Colonization resistance. Antimicrob Agents Chemother. 1994;38(3):409–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang L, Christophersen CT, Sorich MJ, Gerber JP, Angley MT, Conlon MA. Low relative abundances of the mucolytic bacterium Akkermansia muciniphila and Bifidobacterium spp. in feces of children with autism. Appl Environ Microbiol. 2011;77(18):6718–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Finegold SM. Desulfovibrio species are potentially important in regressive autism. Med Hypotheses. 2011;77(2):270–4.View ArticlePubMedGoogle Scholar
- Rodriguez-Brito B, Li LL, Wegley L, Furlan M, Angly F, Breitbart M, Buchanan J, Desnues C, Dinsdale E, Edwards R, et al. Viral and microbial community dynamics in four aquatic environments. Isme Journal. 2010;4(6):739–51.View ArticlePubMedGoogle Scholar
- Reyes A, Haynes M, Hanson N, Angly FE, Heath AC, Rohwer F, Gordon JI. Viruses in the faecal microbiota of monozygotic twins and their mothers. Nature. 2010;466(7304):334–U381.View ArticlePubMedPubMed CentralGoogle Scholar
- Barr JJ, Auro R, Furlan M, Whiteson KL, Erb ML, Pogliano J, Stotland A, Wolkowicz R, Cutting AS, Doran KS, et al. Bacteriophage adhering to mucus provide a non-host-derived immunity. Proc Natl Acad Sci U S A. 2013;110(26):10771–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Barr JJ, Youle M, Rohwer F. Innate and acquired bacteriophage-mediated immunity. Bacteriophage. 2013;3(3):10771–6.View ArticleGoogle Scholar
- Ogilvie LA, Jones BV. The human gut virome: a multifaceted majority. Frontiers in Microbiology. 2015;6.
- Norman JM, Handley SA, Baldridge MT, Droit L, Liu CY, Keller BC, Kambal A, Monaco CL, Zhao G, Fleshner P, et al. Disease-specific alterations in the enteric virome in inflammatory bowel disease. Cell. 2015;160(3):447–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Howe A, Ringus DL, Williams RJ, Choo ZN, Greenwald SM, Owens SM, Coleman ML, Meyer F, Chang EB. Divergent responses of viral and bacterial communities in the gut microbiome to dietary disturbances in mice. Isme Journal. 2016;10(5):1217–27.View ArticlePubMedGoogle Scholar
- Vrieze A, Out C, Fuentes S, Jonker L, Reuling I, Kootte RS, van Nood E, Holleman F, Knaapen M, Romijn JA, et al. Impact of oral vancomycin on gut microbiota, bile acid metabolism, and insulin sensitivity. J Hepatol. 2014;60(4):824–31.View ArticlePubMedGoogle Scholar
- Freedberg DE, Toussaint NC, Chen SP, Ratner AJ, Whittier S, Wang TC, Wang HH, Abrams JA. Proton pump inhibitors alter specific taxa in the human gastrointestinal microbiome: a crossover trial. Gastroenterology. 2015;149(4):883–U531.View ArticlePubMedPubMed CentralGoogle Scholar
- Huttenhower C, Gevers D, Knight R, Abubucker S, Badger JH, Chinwalla AT, Creasy HH, Earl AM, FitzGerald MG, Fulton RS, et al. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486(7402):207–14.View ArticleGoogle Scholar
- Gorrindo P, Williams KC, Lee EB, Walker LS, McGrew SG, Levitt P. Gastrointestinal dysfunction in autism: parental report, clinical evaluation, and associated factors. Autism Res. 2012;5(2):101–8.View ArticlePubMedPubMed CentralGoogle Scholar