The Gut-Brain Axis in Healthy Females: Lack of Significant Association between Microbial Composition and Diversity with Psychiatric Measures

Abstract

Objective

This study examined associations between the composition and diversity of the intestinal microbiota and measures of depression, anxiety, eating disorder psychopathology, stress, and personality in a group of healthy adult females.

Methods

Female participants (n = 91) ages 19–50 years with BMI 18.5–25 kg/m2 were recruited from central North Carolina between July 2014 and March 2015. Participants provided a single fecal sample and completed an online psychiatric questionnaire that included five measures: (i) Beck Anxiety Inventory; (ii) Beck Depression Inventory-II; (iii) Eating Disorder Examination-Questionnaire; (iv) Perceived Stress Scale; and (v) Mini International Personality Item Pool. Bacterial composition and diversity were characterized by Illumina sequencing of the 16S rRNA gene, and associations were examined using Kendall’s tau-b correlation coefficient, in conjunction with Benjamini and Hochberg’s False Discovery Rate procedure.

Results

We found no significant associations between microbial markers of gut composition and diversity and scores on psychiatric measures of anxiety, depression, eating-related thoughts and behaviors, stress, or personality in a large cohort of healthy adult females.

Discussion

This study was the first specifically to examine associations between the intestinal microbiota and psychiatric measures in healthy females, and based on 16S rRNA taxonomic abundances and diversity measures, our results do not suggest a strong role for the enteric microbe-gut-brain axis in normal variation on responses to psychiatric measures in this population. However, the role of the intestinal microbiota in the pathophysiology of psychiatric illness may be limited to more severe psychopathology.

About the Authors:

Susan C. Kleiman

Affiliation: Department of Nutrition, The University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America

Emily C. Bulik-Sullivan

Affiliation: Department of Medicine, The University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America

Elaine M. Glenny

Affiliation: Department of Nutrition, The University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America

Stephanie C. Zerwas

Affiliation: Department of Psychiatry, The University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America

Eun Young Huh

Affiliation: Department of Psychiatry, The University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America

Matthew C. B. Tsilimigras

Affiliation: Department of Bioinformatics and Genomics, The University of North Carolina at Charlotte, Charlotte, North Carolina, United States of America

Anthony A. Fodor

Affiliation: Department of Bioinformatics and Genomics, The University of North Carolina at Charlotte, Charlotte, North Carolina, United States of America

Cynthia M. Bulik

Affiliations Department of Nutrition, The University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America, Department of Psychiatry, The University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America, Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Stockholm, Sweden

ORCID http://orcid.org/0000-0001-7772-3264

Ian M. Carroll

* E-mail: [email protected]

Affiliations Department of Medicine, The University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America, Center for Gastrointestinal Biology and Disease, The University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America

Introduction

Investigations conducted over the last decade have generated consensus among researchers that the intestinal microbiota plays a vital role in a range of physiologic processes, especially those related to immunologic and metabolic function. Emerging evidence also suggests that a healthy intestinal microbiota is important for normal brain development [1]. The enteric microbe-gut-brain axis has garnered increasing attention as a key, bidirectional communication pathway that influences mood, cognition, and behavior [2–4]. In addition to a direct connection via the vagus nerve, it may be possible that gut bacteria interact with the brain through production of neurotransmitters, hormones, and other metabolites [5].

Whether a dysbiosis in the intestinal microbiota, an unhealthy change in the normal microbial ecology of the gut, plays a direct role in the pathophysiology of psychiatric disorders remains to be determined; however, both preclinical animal studies and clinical human studies are actively investigating this question. Numerous studies in animal models have documented behavioral changes following manipulation of the intestinal microbiota, including effects on behavior associated with stress [6], anxiety [1, 7–9], and depression [5, 10]. Validating animal models of the enteric microbe-gut-brain axis in human populations has reported modest associations and has been limited by small sample sizes and lack of consistency in assessment of psychiatric and microbial outcome measures. Although underpowered, these studies in human cohorts suggest a potential role for the intestinal microbiota in anxiety, depression, stress, cognitive reactivity, and eating disorders [11–15].

The majority of the emerging data reporting a relationship between the intestinal microbiota and psychological/behavioral traits has been demonstrated in either rodents or in humans who suffer from threshold psychiatric disorders [16]. It is unknown whether these associations hold across the entire spectrum of severity of psychopathology, or if they are confined to the more pathological extremes. Our goal, therefore, was to determine whether the associations between the intestinal microbiota and psychological measures hold across the spectrum of severity of our target psychological dimensions in healthy individuals. We examined associations between the composition and diversity of the intestinal microbiota and measures of depression, anxiety, eating disorder psychopathology, stress, and personality in a group of healthy adult females.

Materials and Methods

The study was approved by the Biomedical Institutional Review Board at the University of North Carolina at Chapel Hill. All participants provided written consent before study participation.

Study Population

To limit the natural variability in the composition of the intestinal microbiota between healthy individuals, we recruited healthy adult females (n = 100) ages 19–50 years with BMI 18.5–25.0 kg/m2 from central North Carolina via listserv announcements, targeted emails, and social media to serve as controls for ongoing research (clinicaltrials.gov—NCT01916538). Participants were recruited between July 2014 and March 2015. Due to possible impact on the intestinal microbiota, potential participants were excluded for the following reasons: (i) history of gastrointestinal tract surgery (other than appendectomy or cholecystectomy); (ii) history of inflammatory bowel diseases, irritable bowel syndrome, or celiac disease; (iii) history of eating disorders (anorexia nervosa, bulimia nervosa, binge-eating disorder); (iv) treatment in the last two months with antibiotics or steroids; (v) intentional use of probiotics during the last two months (via food or supplement); and/or (vi) abuse of laxatives within the last month.

Body Composition and Assessments

Participants self-reported current height and weight during the screening process. Participants completed an online psychiatric questionnaire that included five widely used and validated measures: (i) Beck Anxiety Inventory (BAI) [17, 18]; (ii) Beck Depression Inventory-II (BDI) [19, 20]; (iii) Eating Disorder Examination-Questionnaire (EDE-Q) [21, 22]; (iv) Perceived Stress Scale-10 (PSS) [23, 24]; and (v) Mini International Personality Item Pool (Mini IPIP) [25].

Sample Collection, Processing, and Storage

During the consent process, participants were provided with an at-home stool collection kit and trained in sample collection procedures. Each kit included: Styrofoam container, disposable collection hat, stool collection tube, biohazard bag, pair of non-latex gloves, two ice packs, and stool collection record sheet. Participants were instructed to return the sample (in the biohazard bag, with ice packs, in the Styrofoam box) to the research office within 24 hours of collection and to keep the sample refrigerated during any interim period. Samples were then immediately transferred to the laboratory, where they were mechanically homogenized with a sterile spatula, aliquoted into sterile 2 ml cryotubes, and stored in a −80°C freezer for future DNA isolation and molecular microbiological analysis.

DNA Isolation

Bacterial DNA was isolated from collected samples using a phenol/chloroform extraction method combined with physical disruption of bacterial cells and a DNA clean-up kit (QIAmp DNA Stool Mini Kit [Qiagen, Valencia, CA]), as previously described [26, 27].

Sequencing of 16S rRNA Genes

Bacterial community composition in isolated DNA samples was characterized by amplification of the V4 variable region of the 16S rRNA gene by polymerase chain reaction (PCR) (forward primer 515, 5'-GA GTG CCA GCM GCC GCG GTA A-3'; reverse primer 806, 5'-ACG GAC TAC HVG GGT WTC TAA T-3'). Forward and reverse primers incorporated single-nucleotide phase shifts of six different lengths each, to improve the quality of the sequence data generated. Generation of 16S rRNA sequences consisted of two separate amplifications: (1) 95°C for three minutes, then 10 cycles of 95°C for 30 seconds, 50°C for 30 seconds, and 72°C for 30 seconds, followed by one cycle of 72°C for five minutes using 120 ng of fecal DNA as template, mixes of the 6 forward and 6 reverse 16S V4 primers at a final concentration of 10 μM, and the KAPA2G Robust PCR kit (Kapa Biosystems, Wilmington, MA); and (2) 95°C for three minutes, then 22 cycles of 95°C for 30 seconds, 50°C for 30 seconds, and 72°C for 30 seconds, followed by one cycle of 72°C for five minutes using 5 μL of purified PCR product from the first amplification as template, 10 μM of forward and reverse primers that contain Illumina MiSeq adapter sequences with a 12-base error-correcting Golay barcode incorporated in the reverse primer, and the KAPA HiFi HotStart ReadyMix PCR kit [28]. Purification of PCR products was carried out after each amplification using the HighPrep PCR clean-up kit (MagBio, Lausanne, Switzerland) with a DynaMag-96 side magnet (Life Technologies, Carlsbad, CA). Equimolar 16S rRNA PCR products were then quantified and pooled for sequencing. Sequencing was performed on an Illumina MiSeq desktop sequencer (Illumina, San Diego, CA) by the High-Throughput Sequencing Facility in the Carolina Center for Genome Sciences at the UNC School of Medicine. Demultiplexed MiSeq data for each of the 91 participants can be accessed on the MG-RAST online metagenomics server [29]: http://metagenomics.anl.gov/linkin.cgi?project=mgp20265 (Project ID: mgp20265).

Analysis of 16S rRNA Sequences

16S rRNA sequencing data were processed by the Quantitative Insights Into Microbial Ecology (QIIME) pipeline [30], with quality filtering as previously described [27]. Forward sequence reads (250 bp) were clustered into Operational Taxonomic Units (OTUs) based on their sequence similarity at a 97% threshold using BLAST and assigned taxonomy using the Greengenes database [31]. Principal coordinates were generated using unweighted and weighted UniFrac distances [32–34].

Results were validated using an alternate pipeline, in which forward reads from the 16S rRNA sequencing data were classified with version 2.10.1 of the RDP classifier with a threshold of a 50% RDP score [35].

Statistical significance was determined using Kendall’s tau-b correlation coefficient in R [36]. R scripts are available at: https://github.com/mcbtBINF/healthyCohort/.

The diversity of the intestinal microbiota was characterized by the Shannon diversity index [37, 38].

Statistical Analysis

Associations between psychiatric and microbial measures were examined using Kendall’s tau-b correlation coefficient, in conjunction with Benjamini and Hochberg’s False Discovery Rate (FDR) procedure to correct for multiple comparisons [39]. Psychiatric measures included: BAI (anxiety), BDI (depression), EDE-Q (total + subscales for dietary restraint, eating concern, shape concern, and weight concern), PSS (stress), and Mini-IPIP (scales for extraversion, agreeableness, conscientiousness, neuroticism, and imagination). Microbial measures included: alpha diversity (Shannon diversity index) and taxa abundance of bacterial groups at the phylum, class, order, family, and genus levels. Linear models were additionally constructed wherein the first two principal coordinates were regressed against each of the psychiatric measures and other participant metadata. The FDR procedure was applied to the number of comparisons per outcome and per taxonomic rank. The α level used was 0.05, but for FDR correction, a more lenient criterion of 0.1 was used [40]. The minimal effect sizes for >80% and >95% power, respectively, are 0.08 and 0.14 at the genus level based on Kendall’s tau-b using simulated data with various degrees of correlation at 10% FDR. All analyses were conducted in R [36]. R scripts are available at: https://github.com/mcbtBINF/healthyCohort/.

Results

Of 100 participants who consented to participate in the study, 94 completed the psychiatric questionnaires and submitted a fecal sample, of which sequencing results from 91 samples met minimal sequencing depth standards for analysis. Demographic and clinical characteristics of the final participant sample (n = 91) are shown in Table 1. In brief, the participants had a mean (SD) age of 29.0 (7.9) years and were within the normal or healthy weight range for adults [41]. On average, their scores indicate normal or minimal levels of anxiety (BAI), depression (BDI), and stress (PSS) and are in line with, or lower than, those of similar non-clinical samples [18, 42–44]. Total scores on the EDE-Q and its four subscales (dietary restraint and eating, weight, and shape concerns) are lower than norms for U.S. college students and young adult females in Sweden and Australia [45–47], which is likely a reflection of the participant recruitment and screening process, which eliminated individuals with a lifetime eating disorder history.

[Figure omitted. See PDF.]

Table 1. Demographic and clinical characteristics of participants in this study (n = 91) as compared to clinical and normative values.

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

Following sequencing of 16S rRNA genes, we had 91 samples with complete data, after excluding those samples with insufficient depth of sequence reads for our downstream analysis.

The total number of 16S rRNA sequence reads was 15,391,194, and the mean number of reads was 169,134 per sample (range: 47,492–317,380 sequence reads).

When examining associations between psychiatric measures and the composition and diversity of the intestinal microbiota, there were no associations that met established significance thresholds. We considered 17 different measures from our human participants—the 15 measures in Table 1 plus participant height and weight. The RDP classifier reported 232 non-rare taxa (12 phyla, 19 classes, 22 orders, 46 families, and 133 genera) that were present in at least 25% of our samples. At each taxonomic level, we also calculated the Shannon diversity index. We therefore evaluated 3,944 hypotheses [17 measures * (232 taxa + 5 Shannon diversity metrics)] using the non-parametric Kendall’s tau-b test for association. Histograms of generated p-values across all possible associations (Fig 1) are largely uniform, suggesting that the null hypothesis of no association is generally supported across all taxonomic levels. Using FDR correction for all 3,944 hypotheses, there were no significant hits even if the threshold were set to 93% FDR (S1 Table).

[Figure omitted. See PDF.]

Fig 1. Histograms of p-values for associations with psychiatric measures by taxonomic level.

Associations between psychiatric and microbial measures were examined using Kendall’s tau-b correlation coefficient, in conjunction with Benjamini and Hochberg’s False Discovery Rate procedure, using data generated by the RDP classifier. Psychiatric measures included: Beck Anxiety Inventory, Beck Depression Inventory-II, Eating Disorder Examination-Questionnaire, Perceived Stress Scale, and Mini-International Personality Item Pool. P-value frequencies were examined at each taxonomic level: (a) phylum; (b) class; (c) order; (d) family; and (e) genus.

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

We also used a less conservative correction, in which associations between each of the 17 human measurements and each taxonomic level were corrected independently (for example, the comparisons of BDI and the 12 phyla were corrected only for the 12 phyla independent of all the other tests that we ran). Even using this much less stringent threshold, where we might expect some spurious correlations, there were no significant hits at a 5% FDR. We conclude that there is a striking lack of correlation between microbial community composition and the measurements we have gathered from this human cohort.

To further visualize the associations in our data set, we generated principal coordinate plots using unweighted UniFrac distances and colored these plots by quartiles of the main psychiatric measures of interest (BAI, BDI, EDE-Q total, PSS) (Fig 2). These plots are based on the first three principal coordinates, which explain 11.5% (PC1), 5.16% (PC2), and 3.98% (PC3) of the variance in microbial composition. The plots do not show evidence of clustering or segregation based on extreme values on psychiatric measures, which further supports a lack of microbial markers for these psychiatric outcomes in this population. The regression of PC1 and PC2 against each of the psychiatric measures and other participant variables did not indicate any significant linear relationships after FDR correction.

[Figure omitted. See PDF.]

Fig 2. Principal coordinate plots of psychiatric measures by quartile.

Principal coordinates were generated using unweighted UniFrac distances from the QIIME pipeline and allocated to quartiles (red: top quartile; orange: middle two quartiles; blue: bottom quartile) based on scores from the (a) Beck Anxiety Inventory; (b) Beck Depression Inventory-II; (c) Eating Disorder Examination-Questionnaire; and (d) Perceived Stress Scale. Plots are based on the first three principal coordinates, which explain 11.5% (PC1), 5.16% (PC2), and 3.98% (PC3) of the variance in microbial composition, and do not cluster by quartile—supporting a lack of association between microbial markers and these psychiatric measures in healthy individuals.

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

Discussion

Our results provide evidence for a lack of association in physically and psychologically healthy adult females between microbial markers of gut composition and diversity and a collection of psychiatric measures, including anxiety, depression, eating-related thoughts and behaviors, stress, and personality. No associations between these measures met established significance thresholds in our analysis. Consistent with our results, recent work in large (>1000) Dutch and Flemish cohorts suggests that effect sizes for a wide variety of clinical and lifestyle variables associated with the microbiome in healthy individuals are on the order of 0.01 or smaller, likely below what would be detectable with our sample size of ~100 [56, 57].

Animal models suggest a role for the intestinal microbiota in anxiety, depression, and stress, and many animal studies have documented behavioral changes following manipulation of the intestinal microbiota using prebiotics, probiotics, antibiotics, infection with pathogenic bacteria, or microbial transfer to germ-free (GF) mice (i.e., mice raised in a sterile environment and lacking an intestinal microbiota). Seminal work by Sudo et al. (2004) on hypothalamic-pituitary-adrenal (HPA) axis activity showed that GF mice have exaggerated stress response when compared to conventionally raised mice [6]. GF mice also have reduced anxiety-like behavior compared to conventional mice [1, 7–9], which can be reversed via early-life colonization with intestinal bacteria [1, 7]. Anxiety-like behavior can also be increased in mice with pathogenic infection [58–60] or transferred between mice with a characteristic anxiety phenotype and non-anxious GF mice using microbial transfer [61]. Probiotic formulations, such as Lactobacillus rhamnosus and Bifidobacterium infantis, have been shown in animal models to reduce depressive and anxiety-like behavior at effect sizes similar to antidepressant treatment [5, 10], and prebiotic human milk oligosaccharides may reduce stress-induced anxiety-like behavior and stimulate changes in microbial diversity [62]. Furthermore, using GF mice, it has recently been discovered that the intestinal microbiota is necessary for appropriate and dynamic regulation of myelin-related genes [63]. Altogether, these findings, together with ours, suggest that changes within the intestinal microbiota may be of central importance to the development or maintenance of depression and anxiety—but that the effects are only observed in more extreme expressions of the traits.

The majority of studies investigating the enteric microbe-gut-brain axis has been in animal models with few attempts to translate these findings to a psychologically healthy human population [16]. Significant results generally lack replication. Mixed evidence has emerged from investigations comparing the intestinal microbiotas of individuals with major depressive disorder to healthy controls, with one study failing to find significant between-group differences in microbial diversity or taxonomic composition [64], while the other found increased diversity and significant taxonomic differences at the phylum, family, and genus levels [11]. In patients with acute anorexia nervosa, which is frequently comorbid with depression, work from our laboratory has shown that microbial diversity was both associated with depression and significantly lower than in healthy controls [12]. Composition and diversity of the intestinal microbiota may also be associated with temperament in young children, but how such links may evolve during the development of adult personality is unclear [65].

Prebiotic and probiotic supplementation has emerged in human clinical studies as potential means for altering mood, but connecting post-intervention changes in mood to differences in microbial composition or diversity is lacking. Studies have reported improvement in measures of depression, anxiety, cognitive reactivity, and stress levels in healthy volunteers after placebo-controlled supplementation trials of prebiotic or probiotic formulas [13–15], but these supplements may not be associated with observable compositional changes to the intestinal microbiota. As probiotic supplementation has been found to have little effect on the composition of the intestinal microbiota, changes in mood or behavior may be mediated by the metatranscriptome (i.e., functional activity of enteric microbes) rather than the intestinal microbiota community composition [66, 67].

Nevertheless, clinical trials of novel treatments for depression based on manipulation of the intestinal microbiota are underway. Seminal work that used Lactobacillus rhamnosus to demonstrate the central importance of the enteric microbe-gut-brain axis [5] is currently being tested in human trials, which are investigating the effects of probiotic supplementation in healthy volunteers and as an augmentation to antidepressant treatment for individuals with treatment-resistant depression [3]. Another ongoing randomized controlled trial is examining the potential benefit of minocycline (a tetracycline antibiotic) treatment as an adjunctive treatment for individuals with moderate to severe depression [68]. Yet because these intervention studies use microbial strains that are not permanent members of the enteric microbiota, it is possible that any evidence of behavioral change is not due to changes to the intestinal microbiota but rather more direct mechanisms of action.

The majority of the emerging data relating to associations between mood and the composition of the intestinal microbiota has been demonstrated in either rodents or in humans with diagnosable psychiatric disorders such as major depressive disorder [16]. Our goal was to see if this phenomenon translates to subthreshold psychological variation in healthy individuals. This would inform us as to whether the link between enteric microbial communities and behavior is found in all individuals—across the full spectrum of a trait (e.g., depression/anxiety)—or is only detectable in individuals at the extreme end of the distribution such as those suffer from a psychiatric illness. To this end, we designed this study explicitly to investigate the relationship between the composition of the intestinal microbiota and continuous measures of psychological traits in healthy individuals. Our findings are in line with an emerging trend in microbiome research: little taxonomic association with healthy variation. Associations between microbial composition and diversity and a wide range of variables reflecting anthropometry, lifestyle, diet, disease, and medication were recently examined in two large European cohorts (Flemish Gut Flora Project: n = 1106; Dutch Lifelines-DEEP: n = 1135). In addition to the high levels of inter-individual variation reported, findings suggest that, while many factors may be associated with microbial variation, any individual factor would have a very modest effect size [56, 57]. These papers did not measure psychiatric variables as broadly or deeply as we did, in which we have demonstrated that effect sizes for healthy cohorts with regard to psychological measures are likely also very small.

These results should be considered in connection with several limitations. With respect to psychopathology, our sample was, on average, healthier than other non-clinical samples of young adults and had less variability on psychiatric measures than would have been expected. As such, restriction of range on psychiatric measures may have played a role in the lack of significant associations with microbial markers. In addition, our analysis focused on taxonomic and diversity measures of 16S rRNA sequencing data, which describes microbial composition but does not account for metabolic activity or functional impact of intestinal bacteria. We may have seen different results with RNA-seq or whole-genome metagenomic shotgun sequencing, longitudinal sampling, or had we analyzed microbial communities from intestinal biopsies. Our sample comprised adult females (age range: 19–50 years), which may differ with respect to these outcome measures from adult males or individuals in younger or older age brackets. Additionally, we excluded participants who had undergone any GI surgery with the exception of appendectomy or cholecystectomy. Recent studies reporting the influence of the appendix and bile acids and salts on the composition of the intestinal microbiota highlight the limitation of including of these surgeries in our study [69, 70]. As we are unable to determine what proportion of our study population had undergone an appendectomy or cholecystectomy, we cannot predict the effect of these surgeries on our results. However, given the relatively low population base rate of these surgeries, and given that our study does not report any significant associations, we believe that excluding subjects with these surgeries would not impact our results or conclusions. Lastly, our study did not take into account that diet or food allergies may influence our observations. Indeed, associations between the enteric microbiota and food intolerance have been reported [71] and changes in diet can alter the composition of this complex community of microorganisms within 24 hours [72].

The enteric microbe-gut-brain axis has attracted considerable attention in recent years, with much focus on the potential role of enteric microorganisms in the development or maintenance of psychiatric illness. Studies involving GF mouse models or clinical populations present extreme cases of psychopathology, which may not reflect microbial mechanisms in a healthy human population.

This study was the first specifically to examine associations between composition and diversity of the intestinal microbiota and psychiatric measures in healthy females, and our results do not reveal associations between the intestinal microbiota and low levels of symptomatology in a healthy population. However, the role of the intestinal microbiota in the pathophysiology of psychiatric illness and evidence of the enteric microbe-gut-brain axis may only be observable in the presence of wider variability of symptom measures and more severe psychopathology.

Supporting Information

[Figure omitted. See PDF.]

S1 Table. Associations between clinical measures and composition and diversity of the intestinal microbiota.

We considered associations between 17 different clinical and psychiatric measures from our human cohort (column B) and 232 bacterial taxa (12 phyla, 19 classes, 22 orders, 46 families, and 133 genera) (column A) present in at least 25% of our samples, as well as the Shannon diversity index. We evaluated 3,944 hypotheses [17 measures * (232 taxa + 5 Shannon diversity metrics)] using the non-parametric Kendall’s tau-b test for association (column D), and there were no associations that met established significance thresholds (column C), even after FDR correction (column E).

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

(TXT)

Acknowledgments

We would like to thank Natalie Bulik-Sullivan for assistance with participant recruitment and sample collection and Amy Perou at the UNC High-Throughput Sequencing Facility for advice and assistance with sequencing preparation.

Author Contributions

1. Conceptualization: SCK SCZ CMB IMC.

2. Data curation: SCK MCBT AAF.

3. Formal analysis: SCK MCBT AAF.

4. Funding acquisition: CMB IMC.

5. Methodology: SCK ECBS EMG EYH IMC.

6. Supervision: CMB IMC.

7. Writing – original draft: SCK ECBS CMB IMC.

8. Writing – review & editing: SCK ECBS EMG SCZ EYH MCBT AAF CMB IMC.

Citation: Kleiman SC, Bulik-Sullivan EC, Glenny EM, Zerwas SC, Huh EY, Tsilimigras MCB, et al. (2017) The Gut-Brain Axis in Healthy Females: Lack of Significant Association between Microbial Composition and Diversity with Psychiatric Measures. PLoS ONE 12(1): e0170208. https://doi.org/10.1371/journal.pone.0170208

References

1. Diaz Heijtz R, Wang S, Anuar F, Qian Y, Bjorkholm B, Samuelsson A, et al. Normal gut microbiota modulates brain development and behavior. Proc Natl Acad Sci U S A. 2011;108(7):3047–52. Epub 2011/02/02. pmid:21282636

2. Dinan TG, Cryan JF. The impact of gut microbiota on brain and behaviour: implications for psychiatry. Curr Opin Clin Nutr Metab Care. 2015;18(6):552–8. pmid:26372511

3. Friedrich MJ. Unraveling the influence of gut microbes on the mind. JAMA. 2015;313(17):1699–701. pmid:25942709

4. Sampson TR, Mazmanian SK. Control of brain development, function, and behavior by the microbiome. Cell Host Microbe. 2015;17(5):565–76. pmid:25974299

5. Bravo JA, Forsythe P, Chew MV, Escaravage E, Savignac HM, Dinan TG, et al. 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. Epub 2011/08/31. pmid:21876150

6. Sudo N, Chida Y, Aiba Y, Sonoda J, Oyama N, Yu XN, et al. Postnatal microbial colonization programs the hypothalamic-pituitary-adrenal system for stress response in mice. J Physiol. 2004;558(Pt 1):263–75. Epub 2004/05/11. pmid:15133062

7. Clarke G, Grenham S, Scully P, Fitzgerald P, Moloney RD, Shanahan F, et al. The microbiome-gut-brain axis during early life regulates the hippocampal serotonergic system in a sex-dependent manner. Mol Psychiatry. 2013;18(6):666–73. Epub 2012/06/13. pmid:22688187

8. Neufeld KA, Kang N, Bienenstock J, Foster JA. Effects of intestinal microbiota on anxiety-like behavior. Commun Integr Biol. 2011;4(4):492–4. Epub 2011/10/04. pmid:21966581

9. Neufeld KM, Kang N, Bienenstock J, Foster JA. Reduced anxiety-like behavior and central neurochemical change in germ-free mice. Neurogastroenterol Motil. 2011;23(3):255–64, e119. Epub 2010/11/09. pmid:21054680

10. Desbonnet L, Garrett L, Clarke G, Kiely B, Cryan JF, Dinan TG. Effects of the probiotic Bifidobacterium infantis in the maternal separation model of depression. Neuroscience. 2010;170(4):1179–88. Epub 2010/08/11. pmid:20696216

11. Jiang H, Ling Z, Zhang Y, Mao H, Ma Z, Yin Y, et al. Altered fecal microbiota composition in patients with major depressive disorder. Brain Behav Immun. 2015;48:186–94. pmid:25882912

12. Kleiman SC, Watson HJ, Bulik-Sullivan EC, Huh EY, Tarantino LM, Bulik CM, et al. The intestinal microbiota in acute anorexia nervosa and during renourishment: relationship to depression, anxiety, and eating disorder psychopathology. Psychosom Med. 2015;77(9):969–81. pmid:26428446

13. Messaoudi M, Lalonde R, Violle N, Javelot H, Desor D, Nejdi A, et al. Assessment of psychotropic-like properties of a probiotic formulation (Lactobacillus helveticus R0052 and Bifidobacterium longum R0175) in rats and human subjects. Br J Nutr. 2011;105(5):755–64. pmid:20974015

14. Schmidt K, Cowen PJ, Harmer CJ, Tzortzis G, Errington S, Burnet PW. Prebiotic intake reduces the waking cortisol response and alters emotional bias in healthy volunteers. Psychopharmacology (Berl). 2015;232(10):1793–801.

15. Steenbergen L, Sellaro R, van Hemert S, Bosch JA, Colzato LS. A randomized controlled trial to test the effect of multispecies probiotics on cognitive reactivity to sad mood. Brain Behav Immun. 2015;48:258–64. pmid:25862297

16. Cryan JF, Dinan TG. Mind-altering microorganisms: the impact of the gut microbiota on brain and behaviour. Nat Rev Neurosci. 2012;13(10):701–12. pmid:22968153

17. Beck AT, Steer RA. Manual for the Beck Anxiety Inventory. San Antonio, TX: Psychological Corporation; 1993.

18. Creamer M, Foran J, Bell R. The Beck Anxiety Inventory in a non-clinical sample. Behav Res Ther. 1995;33(4):477–85. pmid:7755538

19. Beck AT, Steer RA, Brown GK. Manual for the Beck Depression Inventory-II. San Antonio, TX: Psychological Corporation; 1996.

20. Carmody DP. Psychometric characteristics of the Beck Depression Inventory-II with college students of diverse ethnicity. Int J Psych Clin Pract. 2005;9(1):22–8.

21. Fairburn CG, Beglin SJ. Assessment of eating disorders: interview or self-report questionnaire? Int J Eat Disord. 1994;16(4):363–70. pmid:7866415

22. Mond JM, Hay PJ, Rodgers B, Owen C, Beumont PJ. Validity of the Eating Disorder Examination Questionnaire (EDE-Q) in screening for eating disorders in community samples. Behav Res Ther. 2004;42(5):551–67. pmid:15033501

23. Cohen S, Williamson G. Perceived stress in a probability sample of the United States. In: Spacapan S, Oskamp S, editors. The Social Psychology of Health: Claremont Symposium on Applied Social Psychology. Newbury Park, CA: Sage; 1988.

24. Taylor JM. Psychometric analysis of the Ten-Item Perceived Stress Scale. Psychol Assess. 2015;27(1):90–101. pmid:25346996

25. Donnellan MB, Oswald FL, Baird BM, Lucas RE. The mini-IPIP scales: tiny-yet-effective measures of the Big Five factors of personality. Psychol Assess. 2006;18(2):192–203. pmid:16768595

26. Carroll IM, Ringel-Kulka T, Keku TO, Chang YH, Packey CD, Sartor RB, et al. Molecular analysis of the luminal- and mucosal-associated intestinal microbiota in diarrhea-predominant irritable bowel syndrome. Am J Physiol Gastrointest Liver Physiol. 2011;301(5):G799–807. Epub 2011/07/09. pmid:21737778

27. Carroll IM, Ringel-Kulka T, Siddle JP, Ringel Y. Alterations in composition and diversity of the intestinal microbiota in patients with diarrhea-predominant irritable bowel syndrome. Neurogastroenterol Motil. 2012;24(6):521–30, e248. Epub 2012/02/22. pmid:22339879

28. Werner JJ, Zhou D, Caporaso JG, Knight R, Angenent LT. Comparison of Illumina paired-end and single-direction sequencing for microbial 16S rRNA gene amplicon surveys. ISME J. 2012;6(7):1273–6. pmid:22170427

29. Meyer F, Paarmann D, D'Souza M, Olson R, Glass EM, Kubal M, et al. The metagenomics RAST server—a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics. 2008;9(1):386.

30. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335–6. Epub 2010/04/13. pmid:20383131

31. McDonald D, Price MN, Goodrich J, Nawrocki EP, DeSantis TZ, Probst A, et al. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012;6(3):610–8. Epub 2011/12/03. pmid:22134646

32. Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005;71(12):8228–35. Epub 2005/12/08. pmid:16332807

33. Lozupone C, Lladser ME, Knights D, Stombaugh J, Knight R. UniFrac: an effective distance metric for microbial community comparison. ISME J. 2011;5(2):169–72. Epub 2010/09/10. pmid:20827291

34. Lozupone CA, Hamady M, Kelley ST, Knight R. Quantitative and qualitative beta diversity measures lead to different insights into factors that structure microbial communities. Appl Environ Microbiol. 2007;73(5):1576–85. Epub 2007/01/16. pmid:17220268

35. Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73(16):5261–7. Epub 2007/06/26. pmid:17586664

36. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2015.

37. Hill MO. Diversity and evenness: A unifying notation and its consequences. Ecology. 1973;54(2):427–32.

38. Haegeman B, Hamelin J, Moriarty J, Neal P, Dushoff J, Weitz JS. Robust estimation of microbial diversity in theory and in practice. ISME J. 2013;7(6):1092–101. pmid:23407313

39. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol. 1995;57:289–300.

40. Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I. Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001;125(1–2):279–84. pmid:11682119

41. CDC. About Adult BMI: Division of Nutrition, Physical Activity, and Obesity, National Center for Chronic Disease Prevention and Health Promotion; 2015 [updated May 15, 2015; cited 2015 December 8]. http://www.cdc.gov/healthyweight/assessing/bmi/adult_bmi/index.html.

42. Wang YP, Gorenstein C. Psychometric properties of the Beck Depression Inventory-II: a comprehensive review. Rev Bras Psiquiatr. 2013;35(4):416–31. pmid:24402217

43. Osman A, Kopper BA, Barrios FX, Osman JR, Wade T. The Beck Anxiety Inventory: reexamination of factor structure and psychometric properties. J Clin Psychol. 1997;53(1):7–14. pmid:9120035

44. Birks Y, McKendree J, Watt I. Emotional intelligence and perceived stress in healthcare students: a multi-institutional, multi-professional survey. BMC Med Educ. 2009;9:61. pmid:19761603

45. Welch E, Birgegard A, Parling T, Ghaderi A. Eating disorder examination questionnaire and clinical impairment assessment questionnaire: general population and clinical norms for young adult women in Sweden. Behav Res Ther. 2011;49(2):85–91. Epub 2010/12/28. pmid:21185552

46. Mond JM, Hay PJ, Rodgers B, Owen C. Eating Disorder Examination Questionnaire (EDE-Q): Norms for young adult women. Behav Res Ther. 2006;44(1):53–62. pmid:16301014

47. Quick VM, Byrd-Bredbenner C. Eating Disorders Examination Questionnaire (EDE-Q): norms for US college students. Eat Weight Disord. 2013;18(1):29–35. pmid:23757248

48. North Carolina State Center for Health Statistics. 2015 BRFSS Survey Results: North Carolina Derived Variables and Risk Factors-Females: Body Mass Index Grouping-Underweight, Recommended Range, Overweight and Obese Web2015 [cited 2016]. http://www.schs.state.nc.us/data/brfss/2015/nc/female/rf1.html.

49. Julian LJ. Measures of anxiety: State-Trait Anxiety Inventory (STAI), Beck Anxiety Inventory (BAI), and Hospital Anxiety and Depression Scale-Anxiety (HADS-A). Arthritis Care Res. 2011;63(S11):S467–S72.

50. Gillis MM, Haaga DAF, Ford GT. Normative values for the Beck Anxiety Inventory, Fear Questionnaire, Penn State Worry Questionnaire, and Social Phobia and Anxiety Inventory. Psychol Assess. 1995;7(4):450–5.

51. Hunt M, Auriemma J, Cashaw ACA. Self-report bias and underreporting of depression on the BDI-II. J Pers Assess. 2003;80(1):26–30. pmid:12584064

52. Mond JM, Myers TC, Crosby RD, Hay PJ, Rodgers B, Morgan JF, et al. Screening for eating disorders in primary care: EDE-Q versus SCOFF. Behav Res Ther. 2008;46(5):612–22. pmid:18359005

53. Cohen S, Kamarck T, Mermelstein R. A global measure of perceived stress. J Health Soc Behavior. 1983;24(4):385–96.

54. Cole SR. Assessment of differential item functioing in the Perceived Stress Scale-10. J Epidemiol Community Health. 1999;53:319–20. pmid:10396541

55. Cooper AJ, Smillie LD, Corr PJ. A confirmatory factor analysis of the Mini-IPIP five-factor model personality scale. Pers Individ Diff. 2010;48(5):688–91.

56. Falony G, Joossens M, Vieira-Silva S, Wang J, Darzi Y, Faust K, et al. Population-level analysis of gut microbiome variation. Science. 2016;352(6285):560–4. pmid:27126039

57. Zhernakova A, Kurilshikov A, Bonder MJ, Tigchelaar EF, Schirmer M, Vatanen T, et al. Population-based metagenomics analysis reveals markers for gut microbiome composition and diversity. Science. 2016;352(6285):565–9. pmid:27126040

58. Goehler LE, Park SM, Opitz N, Lyte M, Gaykema RP. Campylobacter jejuni infection increases anxiety-like behavior in the holeboard: possible anatomical substrates for viscerosensory modulation of exploratory behavior. Brain Beh Immun. 2008;22(3):354–66. Epub 2007/10/09.

59. Lyte M, Li W, Opitz N, Gaykema RP, Goehler LE. Induction of anxiety-like behavior in mice during the initial stages of infection with the agent of murine colonic hyperplasia Citrobacter rodentium. Physiol Behav. 2006;89(3):350–7. Epub 2006/08/05. pmid:16887154

60. Lyte M, Varcoe JJ, Bailey MT. Anxiogenic effect of subclinical bacterial infection in mice in the absence of overt immune activation. Physiol Behav. 1998;65(1):63–8. Epub 1998/11/12. pmid:9811366

61. Bercik P, Denou E, Collins J, Jackson W, Lu J, Jury J, et al. The intestinal microbiota affect central levels of brain-derived neurotropic factor and behavior in mice. Gastroenterology. 2011;141(2):599–609, e1–3. Epub 2011/06/21. pmid:21683077

62. Tarr AJ, Galley JD, Fisher SE, Chichlowski M, Berg BM, Bailey MT. The prebiotics 3'Sialyllactose and 6'Sialyllactose diminish stressor-induced anxiety-like behavior and colonic microbiota alterations: Evidence for effects on the gut-brain axis. Brain Behav Immun. 2015;50:166–77. pmid:26144888

63. Hoban AE, Stilling RM, Ryan FJ, Shanahan F, Dinan TG, Claesson MJ, et al. Regulation of prefrontal cortex myelination by the microbiota. Transl Psychiatry. 2016;6:e774. pmid:27045844

64. Naseribafrouei A, Hestad K, Avershina E, Sekelja M, Linlokken A, Wilson R, et al. Correlation between the human fecal microbiota and depression. Neurogastroenterol Motil. 2014;26(8):1155–62. Epub 2014/06/04. pmid:24888394

65. Christian LM, Galley JD, Hade EM, Schoppe-Sullivan S, Kamp Dush C, Bailey MT. Gut microbiome composition is associated with temperament during early childhood. Brain Behav Immun. 2015;45:118–27. pmid:25449582

66. Kristensen NB, Bryrup T, Allin KH, Nielsen T, Hansen TH, Pedersen O. Alterations in fecal microbiota composition by probiotic supplementation in healthy adults: a systematic review of randomized controlled trials. Genome Med. 2016;8(1).

67. McNulty NP, Yatsunenko T, Hsiao A, Faith JJ, Muegge BD, Goodman AL, et al. The impact of a consortium of fermented milk strains on the gut microbiome of gnotobiotic mice and monozygotic twins. Sci Transl Med. 2011;3(106):106ra.

68. Dean OM, Maes M, Ashton M, Berk L, Kanchanatawan B, Sughondhabirom A, et al. Protocol and rationale-the efficacy of minocycline as an adjunctive treatment for major depressive disorder: a double blind, randomised, placebo controlled trial. Clin Psychopharmacol Neurosci. 2014;12(3):180–8. pmid:25598820

69. Ridlon JM, Kang D-J, Hylemon PB. Bile salt biotransformations by human intestinal bacteria. J Lipid Res. 2005;47(2):241–59. pmid:16299351

70. Guinane CM, Tadrous A, Fouhy F, Ryan CA, Dempsey EM, Murphy B, et al. Microbial composition of human appendices from patients following appendectomy. MBio. 2013;4(1).

71. Blázquez AB, Berin MC. Microbiome and food allergy. Transl Res. 2016.

72. David LA, Maurice CF, Carmody RN, Gootenberg DB, Button JE, Wolfe BE, et al. Diet rapidly and reproducibly alters the human gut microbiome. Nature. 2013;505(7484):559–63. pmid:24336217