Possible association of Firmicutes in the gut microbiota of patients with major depressive disorder

Abstract

Background

Gut microbiota can affect human behavior and mood in many ways. Several studies have shown that patients with depression were also accompanied with gut microbiota disorder, in which Firmicutes are related to the protective function of intestinal barrier. In this study, we explore the changes and effects of Firmicutes in the patients with major depressive disorder (MDD).

Method

We recruited 54 subjects, including 27 patients with MDD. Fecal samples were collected for identification by 16S rRNA sequencing and bioinformatics analysis.

Results

The study shows that the alpha diversity indices of MDD patients are lower than those of the healthy controls. Firmicutes is the most significantly decreased phylum in the MDD samples. There are totally 13 taxonomic biomarkers with P-value <0.01 from Firmicutes. There are differences in 17 KEGG pathways between the two groups.

Conclusion

This study found that there is a significant disorder of gut microbiota in the patients with depression, in which the Firmicutes decreased significantly. Defects of the Firmicutes may lead to the depression in short-chain fatty acids, which could account for the physiological basis of low-level inflammation of depression.

Limitations

This is a cross-sectional study and the sample size is comparatively small. Though several diet-related factors were controlled in the study, there is no quantified assessment of it.

Keywords:

• gut microbiota
• brain–gut axis
• depression
• Firmicutes
• short-chain fatty acids

Introduction

Depression is a serious disabling health problem with high incidence worldwide;¹ however, the mechanism of its occurrence and development remains unclear. Recent studies suggested that the microbial–gut–brain axis may affect people’s mood and behavior in a variety of ways. By interacting with the vagus nerve, directly changing the function of the central nervous system, affecting the intestinal nervous system, changing the plasticity of the brain,² activating the immune system, and even in more ways,^3,4 these conditional pathogenic bacteria can cause the disease. Increasing evidences were found which lead to the belief that the association between gut microbiota and depression is significant.

In the case of mice depression model, the changes of gut microbiota and fecal metabolic phenotype were found to correlate with depression through 16S rRNA sequencing and research methods based on liquid chromatography–mass spectrometry metabolomics.⁵ Moreover, three studies have shown that germ-free mice displayed more depression-like behavior after being transplanted gut microbiota from depressed people.^6–8 These animal experiments suggest that the disorder of gut microbiota may cause depression. Furthermore, growing evidences of continuous low-level immune inflammatory reaction also cannot be ignored in the pathological process of depression development,^5,9 since the source of this immune inflammation reaction is probably related to the disorder of gut microbiota. First, the bacteria of Firmicutes in gut microbiota can ferment carbohydrates to a variety of short-chain fatty acids (SCFAs),¹⁰ and the lack of these SCFAs can lead to decreased intestinal barrier function.¹¹ Then, when many conditional pathogens and their metabolites in the intestinal tract cross the barrier, and stimulate the immune response, the “gut leaky” is formed, which can affect the occurrence and development of the disease.¹² This can be supported by the study of Yu et al, which showed that there was a significant decrease in Firmicutes in depressed mice.¹³ Another study also found a significant correlation between stress-induced behavioral changes in mice and the disorder of Firmicutes in gut microbiota.¹⁴ In patients with inflammatory bowel disease (IBD), the amount of Faecalibacterium prausnitzii in Firmicutes is minimal, and the decrease in the proportion of bacteria was associated with the decrease in intestinal mucosal protection function.¹⁵ These studies suggest that Firmicutes, as a protective factor of the intestine, deserve further exploration.

It can be observed that Firmicutes and Bacteroidetes are still two major focuses in the human studies related to the gut microbiota and depression. At different levels, certain difference in gut microbiota has been shown between patients and the healthy control (HC) group, but the results of studies on Firmicutes are inconsistent. In the study by Jiang et al,¹⁶ it was found that there was a significant decrease in Firmicutes. However, in three other studies, there was no obvious difference of Firmicutes at phylum level.⁶ Furthermore, some bacteria associated with Firmicutes were decreased greatly at lower levels while others displayed certain increment. The inconsistency of these findings may be due to the following factors: 1) The reference HC group was not entirely normal. 2) The individual health status of the recruited patients was different. 3) The patients’ age span varied in those studies. 4) The effects of related treatment. 5) Differences in diet between typical symptoms and atypical symptoms of depression. Although the results of various studies are inconsistent, the disorder of Firmicutes can still be considered as one of the characteristics of patients with depression.

In order to explore a more certain correlation between the disorder of Firmicutes and the occurrence and development of depression, we adjusted inclusion criteria to better limit the possible interference of the abovementioned factors on gut microbiota so as to avoid the inconsistency appeared in previous studies. We are aiming to clarify the changes of Firmicutes in patients with depression and their related effects.

Materials and methods

Participants

This study was approved by the ethics committee of Peking University Sixth Hospital and Beijing Hospital of Chinese Traditional and Western Medicine. Clinical information was collected at Beijing Hospital of Chinese Traditional and Western Medicine. All subjects signed their written informed consent before participation. The clinical information and sample collection were conducted after obtaining all subjects’ informed consent, and the whole procedures conformed to the directives of the Declaration of Helsinki.

We recruited subjects following a redesigned inclusion criteria with emendation to previous studies.^6,8,16–18 Some adjustments were made according to the specific medical standards in Beijing region. From March 30 to June 30, 2018, we recruited 30 patients with depression, in which 27 met the study criteria and formed the major depressive disorder (MDD) group; then 27 healthy subjects were selected as the HC group according to the MDD group’s age and sex. Both groups are Han Chinese residents living in Beijing for a long time, with no special eating habits and their BMI range from 18 to 30 kg/m². The MDD group met the diagnostic criteria of ICD-10 MDD; they were at the first episode and without systemic antidepressive treatment. The depressive episodes caused by organic and substance abuse and those with atypical characteristics were excluded. The HC group was assessed by two certified attending physician to exclude any other mental illness.

In addition, by carefully examining the approaches applied in former researches,¹⁹ we put stricter limitation on the exclusion criteria as well. We reviewed the previous medical data provided by the subjects and none of the following subjects were included in this study: 1) suffering from any other chronic disease that may affect the stability of the gut microbiota, such as hypertension, diabetes mellitus, metabolic syndrome, immune deficiency, autoimmune disease, cancer, IBD, diarrhea in the last 3 months; 2) drugs that affect gut microbiota have been used in the last 6 months, including antibiotics, glucocorticoids, cytokines, large doses of probiotics and biological agents, and so on; 3) gastroscopy, colonoscopy, or barium meal in the digestive tract were performed in the last 6 months; 4) people who had major gastrointestinal surgery (cholecystectomy, appendicectomy, and intestinal tract resection) in the past 5 years; 5) people with restricted movement due to a major physical or mental illness; 6) people who had experienced significant dietary changes in the last 6 months; and 7) gestating women.

Clinical information collection

We collected general information of all subjects through questionnaires. The general information includes age, gender, race, height, weight, past medical history, drug history, smoking history, and drinking history.

16S rRNA amplification and sequencing

Participants’ fecal samples were put in sterile containers by themselves and collected at the fecal center by specialist. All 54 fresh fecal samples were stored at −80°C prior to DNA extraction. DNA was extracted from 200 mg fecal sample using PowerSoil DNA Kit (Missouri Biotechnology Association, Jefferson, MO, USA) and operating as per manufacturer’s instructions. The V3–V4 region of 16S rRNA was amplified and observed under universal primer pairs 341F (5′-GGACTACHVGGGTWTCTAAT-3′) and 805R (5′-ACTCCTACGGGAGGCAGCAG-3′) by KAPA HiFi HotStart ReadyMix (KAPA Biosystems, Inc., Wilmington, MA, USA). Unique 8 nt barcodes were added to primers in different samples. PCR was conducted at cycling conditions: 95°C for 5 minutes, 20 cycles of 98°C for 20 seconds, 58°C for 30 seconds, 72°C for 30 seconds, and 72°C for 5 minutes. We added 10 pmol primers and 100 ng templates to the 50 µL PCR reactions, then the PCR was performed in triplicate and the PCR products were pooled. QIAquick Gel Extraction Kit (QIAGEN, Hilden, Germany) was used to select DNA segments in proper size. All selected DNA segments were sequenced in paired-end mode using Illumina HiSeq2500 in Novogene Bioinformatics Institute, Beijing, China.

Statistical analyses

Demographic analysis

The SPSS 23.0 statistical package for Windows was used for data analysis. Demographic data and clinical characteristics were compared between groups. Continuous variables were performed by independent samples t-test. Level of significance was set at 0.05 (two-tailed).

Sequencing data analysis

Raw reads were demultiplexed using seqtk (https://github.com/lh3/seqtk). Paired-end reads were merged using FLASH and quality filtered using Trimmomatic: pairs with >15 nt overlap were merged;^20,21 merged sequences were trimmed where average quality score over 4 bases window was <20, and sequences containing ambiguous bases or <400 bp were removed. All qualified sequences were pooled and operational taxonomic units (OTUs) were picked using pick_open_reference_otus.py script from QIIME 1.9.1,²² and chimeras were removed by aligning sequences to the “gold” reference database using UCHIME. OTU sequences were assigned to taxonomy using assign_taxonomy.py from QIIME. All representative sequences (OTUs) were mapped against Greengenes database²³ using UCLUST algorithm with 97% identity.²⁴ Representative sequences were aligned by mafft,²⁵ and phylogenetic tree was generated by FastTree using QIIME.²⁶ Singletons and OTUs only occur in one sample were removed, and OTU table was rarefied using QIIME.

ACE, Chao1, and Shannon diversity values were calculated using vegan,²⁷ and statistical tests were performed using R.²⁸ Faith’s phylogenetic diversity was analyzed using alpha_diversity.py and compare_alpha_diversity.py from QIIME. Unweighted and weighted Unifrac distances were calculated using beta_diversity.py from QIIME, and principal coordinates analysis (PCoA) was performed using R. Significance test of Faith’s phylogenetic diversity was performed using Monte Carlo permutation test in QIIME, and all other significance tests were performed using Wilcoxon test in R.

Taxonomic biomarkers of HC and MDD groups were analyzed using LEfSe (linear discriminant analysis [LDA] Effect Size), and the taxa with P-value <0.01 and LDA score >2.0 were picked out as biomarkers. The metagenome functional profiling was predicted using PICRUSt,²⁹ and the de novo OTUs were removed before the functional profiling prediction according to the manual of PICRUSt. Predicted KOs (KEGG orthology) and pathways were analyzed using STAMP,³⁰ and P-value <0.01 was used to pick out differential KOs and pathways between HC and MDD samples.

Results

Demographic data and clinical characteristics of the subjects

We totally recruited 54 subjects, including 27 patients with MDD and 27 HCs; both groups had the same male to female ratio of 7:20. The mean age of the patient group was 48.7±12.8 and HCs was 42.3±14.1. As the table shows, there was no significant difference in age, height, weight, and BMI between the two groups (Table 1).

Table 1 Demographic and clinical characteristics

Characteristics	MDD (mean ± SD)	HC (mean ± SD)	t-value	P-value
Age	48.7±12.8	42.3±14.1	1.749	0.086
Height	1.65±0.1	1.65±0.1	0.390	0.698
Weight	65.4±11.0	63.6±8.5	0.693	0.492
Body mass index	23.8±2.8	23.4±2.9	0.488	0.628

Abbreviations: HC, healthly control; MDD, major depressive disorder.

OTU picking

Raw sequencing read pairs of samples range from 11,015 to 1,035,838, and the length of reads is 250 bp. After merging pair-end reads, quality filtering, and OTU clustering, the available sequences of samples are between 3,505 and 662,238. The utilization rate of total reads is 52.26%. After OTUs picking and taxonomy assignment, 2,888 OTUs were picked from all sequences, and 183 known taxa were identified.

Lower gut microbiota diversity in MDD patients

Our results show that alpha diversity indices of HC are higher than those of MDD patients (Figure 1). Chao1 and ACE diversity indices can be used to assess the species richness of samples. These two indices are both significantly higher in HC than in MDD (P<0.0008, Wilcox test), indicating that there are richer species in healthy people. The Shannon index can be used to estimate the species evenness and richness of sample, which is significantly higher in HC than in MDD samples (P=0.003, Wilcox test), indicating that healthy people have higher species diversity. The Faith’s phylogenetic diversity index can be used to estimate the phylogenetic diversity of species within a sample, and this index is also significantly higher in HC than in MDD samples (P=0.04, Monte Carlo permutation test). These all indicate that there is a significant decrease in the diversity of gut microbiota in MDD than in HC people. The PCoA plot based on weighted Unifrac distance also shows that the samples from MDD and HC are obviously different in community profile (Figure 2). After taxonomy assignment, the relative abundance of Bacteroides and Firmicutes are the highest two phyla in both HC and MDD samples, which together add up to 92% relative abundance in HC samples and 90% in MDD samples (Figure 3A). One other major difference between HC and MDD samples is the percentage of Firmicutes phylum (Figure 3B). The average relative abundance of Firmicutes in HC samples is 43.46%, while in MDD samples it is only 28.72% (P=0.00016, Wilcox test).

Figure 1 Alpha diversity of HC and MDD samples.

Notes: (A–D) The ACE, Chao1, Shannon, and Faith’s phylogenetic diversity of HC and MDD calculated using rarefied OTU matrix. All four diversity indices in HC are significantly higher than MDD. ^*P<0.05, ^**P<0.01, ^***P<0.001.
Abbreviations: HC, healthly control; MDD, major depressive disorder; OTU, operational taxonomic unit.

Figure 2 Beta diversity of HC and MDD.

Notes: (A and B) PCoA plot of unweighted Unifrac and weighted Unifrac distance of HC and MDD samples using rarefied OTU matrix. Green dots represent HC samples and red dots represent MDD samples.
Abbreviations: HC, healthly control; MDD, major depressive disorder; OTU, operational taxonomic unit; PCoA, principal coordinates analysis.

Figure 3 Taxa in phylum level of HC and MDD.

Notes: (A) The chart plot of taxa in HC and MDD samples. (B) The statistical test of Firmicutes relative abundance between HC and MDD samples. ^***P<0.001.
Abbreviations: HC, healthly control; MDD, major depressive disorder.

Taxonomic biomarkers in HC are all from Firmicutes

In total, there are 13 taxonomic biomarkers found with P-value <0.01 (Kruskal–Wallis test) and LDA score (log 10) >2.0, and among which seven are enriched in HC and six are enriched in MDD (Figure 4). The six biomarkers in HC are all from Firmicutes, including Lachnospiraceae, Ruminococcaceae, Coprococcus, Blautia, Clostridiaceae, and Dorea. The six biomarkers enriched in MDD are from Proteobacteria (Oxalobacter and Pseudomonas) and Firmicutes (Parvimonas, Bulleidia, Peptostreptococcus, and Gemella).¹ This suggests that Firmicutes is the most important phylum that is correlated to depression.

Figure 4 Taxonomic biomarkers in HC and MDD.

Notes: (A) Taxonomic biomarkers found by LEfSe in HC (green) and MDD (red). All taxa enriched in HC are from Firmicutes. (B) Cladogram plot of biomarkers. The size of node represents the abundance of taxa. Only taxa with P-value <0.01 and LDA score (log 10) >2 are shown.
Abbreviations: HC, healthly control; LEfSe, LDA effect size; MDD, major depressive disorder.

Function profiling prediction

There are eleven KEGG pathways enriched in MDD (P<0.01, Welch’s t-test), including lipopolysaccharide biosynthesis, ubiquinone and other terpenoidquinone biosynthesis, glycosaminoglycan degradation, glycosphingolipid biosynthesis, toluene degradation, cellular antigens, protein digestion and absorption, steroid hormone biosynthesis, lipoic acid metabolism. Six pathways are enriched in HC, including sporulation, bacterial motility proteins, bacterial chemotaxis, nitrotoluene degradation, germination, synthesis and degradation of ketone bodies (Figure 5). These micro-biota change in MDD patients and the effects in metabolite could be further explored in future studies.

Figure 5 Predicted differential KEGG pathways in HC and MDD.

Notes: The extended error bar plot of significantly differential KEGG pathways predicted using PICRUSt. Only P-value <0.01 are shown.
Abbreviations: HC, healthly control; MDD, major depressive disorder.

Discussion

Firmicutes account for about 40%–65% of the colon or fecal microbiota. According to the results of previous 16S rRNA sequencing, the dominant abundance flora includes three major Clostridium clusters (IV, IX, and XIV), while other clusters have lower abundance.¹⁰ In our study, it shows that the overall content of Firmicutes in patients with depression was significantly lower than that in the healthy group; this is consistent with the findings of Jiang et al.¹⁶ At the genus level, the significantly reduced genus of Firmicutes mainly fell into three families, which are the Faecalibacterium of the Ruminococcaceae and the Dorea, while Coprococcus of Lachnospiraceae have the most significant difference (P<0.001). These genera belong to Clostridium cluster IV and XIVa, respectively, and can metabolize various carbohydrate substrates to form various SCFAs such as acetate, butyrate, and lactate.¹⁰ The reduction of these fermentation-related bacteria leads to a decrease in the production of SCFA, which in turn leads to intestinal barrier dysfunction.¹¹ This natural barrier function is weakened, multiple antigenic substances are exposed, and the weak intestinal tract becomes the source of inflammation.

Previous studies emphasized that SCFAs produced in the gut play an important role in improving chronic inflammatory diseases and promoting colonic epithelial cells. It has been reported that SCFAs can inhibit the production of proinflammatory cytokines, enhance IL-10 expression, activate regulatory T cells (Tregs), and alleviate colonic inflammation.^31,32 SCFAs mainly include acetate, propionate and butyric acid, which have significant effects on intestinal epithelial cell proliferation, differentiation, and metabolism. Among them, butyrate not only can provide energy for the long epithelium but also strengthen the colon defense barrier. In addition, butyric acid can also play a role in immunoregulating cell cycle inhibition, inducing programmed cell death, and cell differentiation in various cell types. Recent evidence suggests that butyrate and propionate are keys to the regulation of Foxp3+ production of Tregs, whereas Tregs play an important role in suppressing inflammatory responses.³³ Since most bacteria that produce butyric acid belong to the Firmicutes,¹¹ with the decrease of Firmicutes, these protective factors are weakened, and the body is further exposed to the risk of inflammation.

A number of studies have indicated that cytokines and inflammation are closely related to depressive symptoms in patients with depression. It has been suggested that depression can be seen as a group of symptoms caused by peripheral inflammation and a response to inflammation.^34,35 A meta-analysis suggests that the concentration of IL-6 and TNF-α in the blood is significantly elevated in depression without any physical illness.³⁶ A large number of longitudinal studies have shown that exogenous cytokines can aggravate depressive symptoms.^37–41 Similarly, injection of lipopolysaccharide endotoxin or related vaccine can increase both proinflammatory cytokine concentrations and depressive symptoms.^42–44 In the mice study of Zhang et al, it was found that there was a significant decrease of Firmicutes in the social defeat stress model, while the change of proteobacteria was not significant. They also found that intravenous injection of MR16-1 induce antidepressant effects by normalizing the altered composition of the gut microbiome.⁴⁵ This is consistent with our conclusion.

In addition, some cell-mediated immune processes may also be involved in the development of depression.⁴⁶ Two meta-analyses have indicated that there are multiple activations of cellular immune pathways in patients with depression.^36,47 Although there is currently no direct evidence that low-grade inflammation in depressive patients is derived from the intestine, there is increasing evidence that gut microbiota is important in causing this inflammatory process. Although some preclinical and clinical studies have confirmed the positive effects of probiotic supplementation on depressive symptoms, a meta-analysis showed that probiotic supplementation has an overall insignificant effect on mood.⁴⁸ Therefore, it is still necessary to further clarify the changes of gut microbiota in depression, which might help the targeted supplementary to achieve a better effect.

In conclusion, our study found that there is a significant disorder of gut microbiota in patients with depression, in which the Firmicutes decreased significantly. Defects of the Firmicutes may lead to the depression in SCFA, which may be the physiological basis for low-level inflammation of depression. In the future, we can further explore the role of Firmicutes in depression through the method of multiomics.

Limitation

This study still has some limitations. First, the sample size we used was comparatively small due to financial limit. Second, although the results of this study support that gut microbiota plays a role in the development of depression, we are currently unable to investigate how exactly gut microbiota changed along with this process. In future study, the changes of gut microbiota should be further observed in high-risk groups throughout their possible development of symptom. Third, there is a lack of relevant inflammatory indicators in this study. Finally, even though we carefully selected subjects aiming to reduce the influence of related factors on gut microbiota, some confounding factors, such as diet, still need more control or detailed assessments. Moreover, atypical symptoms of depression, such as greediness and somnolence, could also have potential impact on gut microbiota, which requests a more detailed classification of depression in future research.

Acknowledgments

We thank all the subjects who took part in this study. All the costs of this study are self-financed.

Disclosure

The authors declare no conflicts of interest in this work.