Background Although the link between gut microbiota and depression has been suggested, changes of gut microbiota vary largely among individuals with depression.
Aims Explore the heterogeneity of microbiota–gut–brain axis and new pathogenic characteristics in murine models of depression.
Methods Adolescent female mice were randomly divided into control (CON) group (n=10), chronic unexpected mild stress (CUMS) group (n=15) and dexamethasone (DEX) group (n=15). Mice in the DEX group were gavaged twice a day with 0.2 mg/kg of DEX for 5 weeks, whereas CON mice were given the same amount of solvent. Mice in the CUMS group were exposed to stressors. After behavioural evaluations, all mice were sacrificed for harvesting tissues and blood samples. Enzyme-linked immunosorbent assay (ELISA) was conducted for measuring levels of corticosterone (CORT) and interleukin-1β (IL-1β) in sera, whereas levels of protein expression in colon and hippocampal tissues were examined by western blot. Faecal microbial communities were analysed by sequencing 16S rDNAs.
Results Mice in CUMS and DEX groups exhibited severe depression-like behaviours. Compared with CON mice, CUMS-exposed mice showed a significant increase in both α and β diversity. Prevotellaceae and Desulfovibrio were enriched, whereas Bacilli were decreased in the faeces of mice in the CUMS group. DEX-treated mice had a decrease in the abundance of Clostridium XVIII. Levels of occludin in colon tissue of DEX-treated mice were reduced. Relative to mice in the CON and CUMS groups, DEX-treated mice contained higher serum levels of CORT and IL-1β. Compared with CON mice, mice in the DEX and CUMS groups had higher levels of IL-1β in sera and lower levels of glial fibrillary acidic protein (GFAP), Nestin, Synapsin-1 and P2Y12 receptor in the hippocampus.
Conclusions Changes of gut microbiota diversity, intestinal integrity and neuroinflammation in the brain contribute to CUMS-induced depression, whereas pathobionts and excessive immunosuppression with damaged neuronal synapses is a basis of the DEX-induced depression.
- depressive disorder
Data availability statement
All data relevant to the study are included in the article or uploaded as supplemental information.
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.
Depression is characterised by depressive mood and loss of interest in daily activities.1 According to the World Health Organization (WHO), there are more than 350 million people of all ages suffering from depression, making it a leading cause of disability over any other diseases.2 Up to now, pharmacological, psychological and physical interventions are the mainstay of depression treatment.3 4 However, a significant population of persons with depression do not respond to these treatments. Depression in children and younger individuals can devastatingly impact their development, ability to adapt and performance at school.5 It is documented that most antidepressants are unable to produce a satisfactory result for young patients. Therefore, further elucidation of the mechanisms that underlie depression is needed.
The microbiota–gut–brain axis, through which the gut and the brain communicate to each other, has been regarded as a great example of biopsychosocial interplay in the pathogenesis of neurological diseases.6 Accumulating evidence suggests that changes in the metabolic states of the host influence the composition of gut microbiota. While varied composition of gut microbiota in patients with depression or animal models of depression has been widely reported, there exists inconsistencies regarding the characteristics of gut microbiota in people with depression compared with healthy adults. Even conflicting results are reported on microbial diversity as well as the abundance of bacteria at phyla, family and genus taxonomic levels in people with depression. Hence, it is necessary to further investigate the role of microbiota–gut–brain axis in the development of depression.
Women are more susceptible to developing depression and anxiety disorders in response to stress than men. Dysfunction of the hypothalamic–pituitary–adrenal (HPA) axis plays a central role in the formation of depression in women.7 There is evidence that depression is more prevalent among females with abdominal obesity and over-activation of the HPA axis.8 Additionally, approximately 50% to 75% of individuals with excessive weight loss due to anorexia nervosa experience lifetime major depressive disorder (MDD).9 We speculate that these subtypes of major depression may arise from alterations in the microbiota–gut–brain axis. In this study, we used adolescent female mice to model depression with two different inducers, specifically dexamethasone (DEX) and chronic unexpected mild stress (CUMS), and then compared the heterogeneity and homogeneity of gut microbiota. We found that there were significant changes in gut microbiota, expression levels of tight junction proteins in the colonic mucosa, as well as signs of neuroinflammation in the hippocampus in mice with depressive behaviours induced by CUMS or DEX. Our study suggests that analysis of the heterogeneity of gut microbiota is helpful to our understanding of depression pathogenesis and to the discovery of targets for treating depression in women.
Materials and methods
Antibodies for glial fibrillary acidic protein (GFAP), Nestin, and Caspase-1 were from Proteintech (Chicago, Illinois, USA). Antibodies against interleukin-1β (IL-1β), glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and horseradish peroxidase conjugated antirabbit and mouse IgG were acquired from Cell Signaling Technology (Boston, Massachusetts, USA). Anti-Claudin-1 and anti-Synapsin-1 were purchased from Bioworld (Dublin, Ohio, USA), whereas anti-Occludin was from Thermofisher (USA). Anti-NLRP3 (nucleotide-binding domain-like receptor protein3) and anti-P2Y12 receptors were purchased from NOVUS (Cambridge, UK) and Abcam (Cambridge, UK), respectively.
All experiments were conducted in 4 to 6-week-old female C57/BL6 mice (14–16 g), which were purchased from Shanghai SLAC Laboratory Animal, Shanghai, China. Mice were maintained in a controlled environment at 22℃±1℃ with a 12/12 hour light–dark cycle. Mice were randomly divided into control (CON) group (n=10), DEX group (n=15) and CUMS group (n=15).
Establishment of depression models
Dysfunction of the HPA axis in female mice was induced by gavage feeding 2 mg/kg DEX twice a day for 5 weeks (figure 1). The CUMS-induced model was developed as previously described.10 In brief, mice were individually caged and subjected to CUMS for 5 weeks with the following stressors: 4℃ ice-water swimming for 5 min, cage tilted by 30° for 24 hours, foot electric shocked for twice (0.8 mA, 5 s), fasting 24 hours, continuous illumination, behavioural restrictions for 2 hours, wet pad for 24 hours, tail-clamp for 90 s, water deprivation for 24 hours and cage shaking for 15 min. To prevent habituation and to provide an unpredictable feature to the stressors, stressors were randomly administered every day for 35 days without the same kind of stimulus. Control mice were fed normally and fed by gavage with the same volume of water.
Behaviours of all mice in each group were evaluated once CUMS and/or DEX exposure was completed. Mice were allowed to adapt to their environments for 1 week before being evaluated. All testing apparatuses were cleaned with 70% alcohol and distilled water to eliminate residual odours before and after each animal was assessed. The experimenters who analysed behaviour data were blinded to the treatment conditions.
Sucrose preference test (SPT) was conducted after 1-week adaption. The SPT, open field test (OFT), elevated plus maze (EPM), tail suspension test (TST) and forced swimming test (FST) were conducted after CUMS and DEX administration were completed. The schedule is shown in figure 1.
Sucrose preference test
Mice were allowed to choose from two bottles of liquids over 24 hours; one bottle contained 1% sucrose solution, whereas the other bottle contained the same volume of pure water. The bottle positions of the two bottles were switched once 12 hours after SPT was started. The sucrose preference was determined as the percentage of the total amount of liquid consumed.
Open field test
The open field was a box (100 cm × 100 cm × 30 cm) with its floor divided into 16 squares; the central four squares were defined as the centre, whereas the rest 12 squares were considered as periphery. Each mouse in each group was placed in the centre of the open arena and allowed to explore for 5 min. Movements of each mouse were recorded with a video-tracking system (Shanghai Mobile Data Information Technology). The total distance moved in the arena, travel speed and rearing frequency were measured.
Elevated plus maze
The apparatus (50 cm height from floor) consisted of two open arms (30 cm×5 cm) and two enclosed arms (30 cm×5 cm×25 cm) and was used for measuring anxiety-like behaviours of animals. Each mouse was placed in the central section facing an open arm and was allowed to explore the arms for 5 min. Activities of each mouse were recorded by video-tracking system (Shanghai Mobile Data Information Technology). The time spent in open arms was calculated.
Tail suspension test
Each mouse was suspended upside down by fixing the tip of its tail with adhesive tape on a metal rod 40 cm above the floor. During the 6 min test period of testing, activities of the mice tested were recorded with a video camera. The immobility time within the last 4 min was calculated.
Forced swimming test
Each mouse in each group was placed in a vertical transparent cylinder (30 cm in height and 12 cm in diameter) containing tap water (25℃±1℃ and 20 cm in depth). Activities of the mice tested over 6 min period of testing were recorded with a video camera. The immobility time within the last 4 min was calculated.
After all behaviour tests were completed, mice were sacrificed by exsanguination. Before being sacrificed, faecal samples were collected in sterilised tubes. Blood samples were collected by puncturing the abdominal aorta. Sera from each mouse were extracted with a centrifugation of 5000 g for 15 min and stored at −80℃. Colon and hippocampal tissues were dissected on ice immediately after euthanisation. Colon tissues were cut open longitudinally and washed with ice-cold phosphate buffer saline (PBS) to remove gut contents. All tissues were all frozen and stored at −80℃ until being used for studies.
Serum levels of corticosterone (CORT) and IL-1β were measured using mouse CORT and IL-1β ELISA KIT (XLPCC company), respectively, according to the supplier’s instructions. In brief, after reagents were restored to room temperature, 50 µL of each sample or standard as well as reagents were added into wells of a microplate sequentially and gently mixed. After incubation for 1 hour in 37℃ incubators, samples in each well were washed five times using the washing buffer. To each sample (well), 50 µL reagent A and B were added and incubated until the reaction was stopped by adding stop buffer. The absorbance at 450 nm was measured in a plate reader (Biotek, USA).
Western blot analysis
Tissues were lysed with ice-cold radioimmunoprecipitation assay (RIPA) lysis buffers (Beyotime Institute of Biotechnology, Jiangsu) supplemented with protease and phosphatase inhibitors (Roche, Indianapolis, Indiana, USA). Lysates were centrifuged at 4℃ 12 000 g for 10 min. The resulting supernatants were collected. Protein concentrations were determined using Bradford reagents (Beyotime Institute of Biotechnology). An aliquot of proteins was mixed with sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) sample buffers, incubated at 95℃ for 8 min and separated on 10% SDS-polyacrylamide gels and transferred onto polyvinylidene difluoride membranes (Millipore, Bedford, Massachusetts, USA). The membranes blots were blocked with 5% bovine serum albumin (BSA) in Tris-buffered saline containing 1% Tween-20 (Tris-buffered salin with Tween-20, TBST) for 2 hours at room temperature and then incubated with primary antibodies overnight at 4℃. Primary antibodies used for western blot analysis included: anti-claudin-1, anti-Occludin, anti-IL-1β, anti-NLRP3, anti-Caspase-1, anti-Nestin, anti-GFAP, anti-synapsin-1 and anti-GAPDH. After washing in TBST, the blots were incubated with corresponding horseradish peroxidase (HRP)-linked secondary antibodies. Protein bands on the blots were detected with chemiluminescence detection reagents (Beyotime Institute of Biotechnology). Optical densities of interested protein bands were obtained using a software (Bio-Rad Laboratory, USA) and protein levels were normalised to GAPDH.
Fecal microbiota analysis
The 16S rDNA identification technology was used for detecting and identifying the composition of intestinal flora of mice through services at the Shanghai Ruiyi Biological Technology Company. The quantity and quality of DNA samples isolated from faeces were examined with Thermo Nanodrop 2000 UV microspectrophotometer and gel electrophoresis, respectively. The V3–V4 region of 16S rDNA was chosen for amplification using the universal primers 341F: 5′ CCTACGGGRSGCAGCAG 3′ and 806R: 5′ GGACTACVVGGGTATCTAATC 3′. The index sequence and linker sequence for Illumina Miseq PE250 sequencing were added to the 5′ end of the universal primer. PCR was performed using the high-fidelity enzyme of Kapa Hifi Hotstart ReadyMix PCR Kit. PCR products were detected by agarose gel electrophoresis and recovered using the gel cutting with AxyPrep DNA Gel Recovery Kit. The quality of the library was determined by using a Thermo Nanodrop 2000 UV spectrophotometer and gel electrophoresis. Qubit was used for library quantitation. According to the data volume requirements of each sample, the corresponding proportion was mixed, and Illumina Miseq PE250 was used for computer sequencing.
Alpha diversity analysis
Alpha diversity was calculated using QIIME (V.1.9.1) software, and the corresponding dilution curve was drawn. The dilution curve uses the relative ratios of various operational taxonomic units (OTUs) known in the measured 16S rDNA sequence to calculate the expected values of each Alpha diversity index when n Reads are extracted (n is less than the total number of measured Reads sequences), and then the curve is made according to the expected values of a set of n values (generally a set of isograms less than the total number of sequences) and the corresponding alpha diversity index. Rank sum test analysis was carried out on each index of alpha diversity (Wilcox. test function in R.V3.5.1 was used for comparison between two groups of samples, and Kruskal. test function in RV.3.5.1 was used for comparison between two groups of samples previously), and Alpha diversity index with significant difference under different conditions was screened by rank sum test.
Beta diversity analysis
Principal component analysis
In order to further show the difference of species diversity among samples, principal coordinate analysis (PCoA) was used to show the difference between samples implemented with the Dudi.pco function in the R (V.3.5.1) language ade4 package.
Anosim (similarity) analysis is a non-parametric test to test whether the difference between groups (two or more groups) is significantly greater than the difference within the group, so as to determine whether the grouping is meaningful. First, the distance between two samples was calculated by Unifrac algorithm, and then all distances were sorted from small to large. The calculation process is realised by anosim function in the vegan package of R (V.3.5.1) language.
Linear discriminant analysis (LDA) effect size (LEfSe analysis V.1.0): LEfSe uses LDA to estimate the impact of each component (species) abundance on the difference effect and to identify communities or species that have significant differences in sample partitioning. The analysis was performed using LEfSe Tools.
Data were expressed using mean (standard deviation, SD). The SPSS statistical software V.16.0 was used for statistical analysis. One-way analysis of variance followed by the Tukey’s test was conducted for comparisons. To reveal differentially abundant features, relative abundance comparisons at levels of genus, family and phylum were performed on non-normalised data, employing non-parametric Mann-Whitney U tests, as well as, Kruskal-Wallis tests with multiple comparison correction. A value of p<0.05 was considered as having statistical significance.
CUMS and DEX treatments induce depression-like behaviours in female mice
We evaluated behaviours of mice with SPT, OFT, EPM, TST and FST to determine whether CUMS and DEX evoked depression-like or anxious behaviours. Figure 2 showed that mice in both CUMS and DEX groups displayed a significant decrease in the total travel distance, movement speed, frequency of rearing and crossing in open arena, as well as the time spent in open arms of the EPM and sucrose preference compared with the mice in the CON group (p<0.01). Mice exposed to CUMS or DEX had a significant increase of immobility time in TST and FST relative to mice in the CON group (p<0.01). While there was no significant difference in body weight between groups at the beginning of experiments (figure 2B), the body weight of mice (p<0.01) was significantly reduced and increased, respectively, relative to that of mice in the CON group (p<0.01). Collectively, these data suggest that we have succeeded in establishing CUMS and or DEX induced models of depression.
Differences in the gut microbiota diversity in two mouse models of depression
To determine whether gut microbiota was altered in CUMS and or DEX treated mice, we collected faeces for analysing bacterial diversity and composition. The dilution curve in online supplemental figure S1 showed that the sequencing depth covered all species in the samples. The α diversity indices including chao1, observed species, shannon and simpson were used as indicators of species richness and uniformity in community ecology. As shown in online supplemental figure S2, there was a significant difference in the abundance and diversity between the CON group and CUMS group (p<0.05), whereas there was no significant difference between the CON group and the DEX group (p>0.05). These data suggest that CUMS but not DEX alters the abundance and richness of intestinal flora. β diversity analysis was used for comparing the differences in species diversity of a pair of intergroup samples. Anosim similarity revealed the differences between intergroups and intragroups. The differences of intergroups were significantly greater than those of intragroups (online supplemental figure S3, R=0.892, p=0.007; R=0.524, p=0.01). To further compare the species diversity, PCoA was conducted (online supplemental figure S4).11 12 As shown in figure 3, the relative abundance in the CUMS group was higher than that in the CON group, but the relative abundance of the DEX group was similar to that in the CON group. These results indicate that gut microbiota of mice subjected to CUMS are significantly different from the CON group, whereas the compositions of intestinal flora between DEX and CON groups are similar.
The species abundance and difference of gut microbiota in CUMS and DEX group
According to the results of species annotation, the corresponding histogram of each sample was analysed at phylum, class, order, family and genus classification levels respectively, so as to intuitively discern the abundance and proportion of species in different classification levels of each sample. At the phylum level, the abundance of Bacteroidetes (72.86 (2.43)%) and of Firmicutes (24.81 (4.24)%) in the CUMS group was higher and lower, respectively, when compared with the CON group (43.25 (14.17)%). There was no significant difference in the abundances of Bacteroidetes and Firmicutes between the DEX and CON groups. At the genus level, the abundances of Bacteroides (p<0.01) and Parabacteroides (p<0.01) were decreased significantly in the DEX group relative to that in the CON group. The abundance of Parabacteroides (p<0.01) as well as Lactobacillus and Staphylococcus (p<0.05) was decreased in the CUMS group. Escherichia/Shigella was found in the DEX group but not the CUMS or the CON group. Moreover, the relative abundances of Alloprevotella, Desulfovibrio, Parasutterella and Flavonifractor were significantly increased in the CUMS group relative to their corresponding abundances in the CON group. LDA effect size (LEfSe) was used to evaluate the impact of the abundance of each component and to identify the communities or species that are different between groups. Compared with the control group, mice in the CUMS group had an increase of the genera in their intestinal microflora, which included:Prevotellaceae, Alloprevotella, Oscillibacter, Desulfovibrionaceae, Mucispirillum, Olsenella, Allobaculum, Candidatus Saccharibacteria and Saccharibacteria genera incertae sedis, and a decrease of Bacilli, Lachnospiracea incertae sedis, Neisseriaceae, Neisseriales, Neisseria, Eubacteriaceae, Eubacterium, Butyricicoccus, Staphylococcus, Bacillales and Staphylococcaceae; Compared with the control group, mice in the DEX group had an increase of Enterobacteriaceae, Escherichia Shigella, Gammaproteobacteria, Enterobacteriales, Anaerofustis and Erysipelotrichaceae_incertae_sedis and a decrease of Clostridium XVIII. Relative to those in the CUMS group, Prevotellaceae, Alloprevotella, Ruminococcaceae, Oscillibacter, Desulfovibrionaceae, Deltaproteobacteria, Desulfovibrionales, Desulfovibrio, Deferribacteres, Deferribacteraceae, Deferribacterales, Burkholderiales, Mucispirillum, Olsenella, Allobaculum, Sutterellaceae, Candidatus Saccharibacteria and Saccharibacteria genera incertae sedis were significantly decreased, whereas Eubacteriaceae, Eubacterium, Enterobacteriaceae, Escherichia/Shigella, Gammaproteobacteria, Streptococcus, Enterobacteriales, Anaerofustis, Erysipelotrichaceae_incertae_sedis, Streptococcaceae, Peptostreptococcaceae, Lachnospiracea_incertae_sedis, Clostridium XI and Firmicutes were increased in the DEX group.
The intestinal barrier integrity and proinflammatory cytokine in the colon of depression mouse models
The intestinal epithelial barrier is the first line of defence of the host. Serious inflammation or many other intestinal diseases are associated with the integrity damage of the intestinal epithelial barrier. Accumulating evidence suggests that patients with chronic colon inflammation frequently have serious depression and anxiety morbidities.13 Therefore, we evaluated the colon epithelial barrier integrity by examining expression levels of claudin-1 and occludin, which are protein markers of the tight junction. While levels of claudin-1 were comparable in colon tissues of mice between CON and DEX groups, levels of occludin were significantly decreased in colon tissues of mice in the DEX group relative to those in colon lysates of mice in the CON group (p<0.05). Levels of IL-1β, which reflects inflammation, were significantly increased in colon tissues of mice in the CUMS group (figure 4A, p<0.01). These results suggest that DEX influences the colon permeability, whereas CUMS evokes inflammation in the gut.
The heterogeneity of serum levels of CORT and IL-1β in CUMS and DEX groups
Glucocorticoids (GCs) are stress hormones released by the adrenal cortex and regulate body metabolism and function. The serum level of corticosterone (CORT) is considered to be an indicator of stress and the severity of depression. The proinflammatory cytokine IL-1β is a cytokine released by macrophages on lipopolysaccharide (LPS)-mediated stimulation of toll-like receptor 4 (TLR4). Compared with that in the control group, the serum level of CORT was significantly increased in the CUMS group (p<0.01). Noticeably, serum levels of both CORT and IL-1β were significantly elevated in the DEX group relative to their levels in CON and CUMS groups (p-values <0.01). The rise of CORT in the serum of mice in the DEX group may partly result from metabolites of dexamethasone (figure 4B).
CUMS and DEX groups mainly affected glia cells instead of neuron in the hippocampus
The hippocampus has been regarded as the main target in the formation of depression.14 Compared with those of the CON group, levels of GFAP and Nestin in the hippocampus of mice in both CUMS and DEX groups were significantly decreased (p-values <0.05). No significant difference in levels of neuron-specific neuronal nuclei (NeuN) among the three groups was noted (p>0.05). However, we found that levels of synapsin-1 and P2Y12 were decreased in the hippocampi of mice in both CUMS and DEX groups when compared with those in the CON group (figure 5, p<0.05). Compared with their expression in the control group, NLRP3 and IL-1β were upregulated in the CUMS group, but only IL-1β was found to increase in the hippocampus of mice in the DEX group (p<0.05). Collectively, these results suggest that the activation of glial cells occurs in the brain of CUMS or DEX induced mouse model of depression.
The microbiota–gut–brain axis is emerging as a particular area of interest and a potential new target for the treatment of central nervous system disorders. Dynamic homeostasis of the intestinal environment is maintained by the inter-relation and interaction between beneficial and pathogenic microbes in the intestinal tract. The dysbiosis could cause long-lasting harmful effects to health by releasing enteric toxins, invading the gut mucosa, triggering inflammation and even interfering with nutrient absorption. Considerable studies have shown that chronic low-grade inflammation increases the susceptibility to depression. We found that compared with the CON group, the CUMS group had a significant change in both α and β diversities of intestinal microbiota, suggestive of dysbiosis in the gut. However, DEX exhibited little effect on the microbiota community. An elevation in the Firmicutes-to-Bacteroidetes ratio signs a greater energy-harvesting capacity from undigested carbohydrates with more lipids production. Obesity is associated with a high proportion of Firmicutes, whereas lean favours the Bacteroidetes. Our results showed that the body weight of CUMS-treated mice significantly declined and coincidently had a significant decrease in the abundance of Firmicutes-to-Bacteroidetes on the phylum level compared with the control group.
On the genus level, we found a high abundance of Alloprevotella, Desulfovibrio, Parasutterella and Flavonifractor, and a low abundance of Parabacteroides, Lactobacillus and Staphylococcus in the CUMS group. In the DEX group, Escherichia/Shigella were highly abundant, and Bacteroides and Parabacteroides were relatively scant. Parabacteroides and Lactobacillus are probiotics and contribute to the intestinal microflora balance and nutrient digestion by regulating bile acid profiles. The abundance of Lactobacillus species in the intestinal tracts of patients with major depression are significantly decreased.15 16 Lactobacillus competes with harmful bacteria in the intestinal mucosa space and nutrition and also promotes the conversion of glutamic acid to γ-aminobutyric acid (GABA).17 GABA is one of the most important inhibitory neurotransmitters in the central nervous system and is involved in the pathogenesis of anxiety and depression associated with functional bowel disorders.15 Increased levels of Desulfovibrio and Parasutterella have been found in models of gout and in patients with inflammatory bowel disease.18 19 Overgrowth of Desulfovibrio in the gut has been suggested to play a role in regressive autism.20–22 Flavonifractor, one of the Bacteroidales order Bacteroidia class Bacteroidetes phylum, displays a high invasion potential in the gut of inflammatory bowel diseases in patients with major depression.23 Alloprevotella is a dominant bacterium in healthy people and produces the short-chain fatty acid. As the abundance of Alloprevotella in the gut and the fasting state are correlated with each other, we speculate that the loss of body weight in the CUMS group might be associated with a higher abundance of Alloprevotella. The most isolated Enterobacteriaceae is Escherichia/Shigella, which was significantly enriched in the DEX group.
Allobacullum and Olsenella are related with intestinal protection effect by producing butyrate. Both are found to decrease in the gut of obese mice with a high fat diet.24–26 We found that Allobacullum and Olsenella were increased in the gut of mice in the CUMS group. The Staphylococcus is included in the Staphylococcaceae family, Bacilli class, Bacillales order and Firmicutes phylum. It is positively associated with energy intake by using amino acids as a major energy source. The decrease of it may further reduce energy utilisation in the CUMS group.15 Additionally, the increased Erysipelotrichaceae_incertae_sedis in gut microbiota composition has also been reported in people with depression. High fat diet-induced obesity is associated with high abundance of Anaerofustis.27 The abundance of Clostridium XVIII in the gut is related to food allergy.28 We found that DEX treatment significantly decreased Clostridium XVIII and increased Anaerofustis, such changes might contribute to the altered body weight of mice in the DEX group. Compared with the DEX group, the CUMS group had higher abundances of Ruminococcaceae, Burkholderiales and Sutterellaceae. Ruminococcaceae belongs to the Clostridia class, which are anaerobic and gram-positive microbes. It has been reported that the abundance of Ruminococcaceae and the abundances of Burkholderiales and Sutterellaceae are increased and decreased, respectively, in the adrenocorticotropic hormone induced rat model of depression.29 In our CUMS-induced mouse model of depression, the abundances of Burkholderiales and Sutterellaceae were increased, likely because of the high levels of oestrogen in female mice. Compared with the CUMS group, the DEX group had higher abundances of Peptostreptococcaceae, Lachnospiraceae and Clostridium XI, which are usually associated with gut inflammation and blood lipids.30
Since most microbiota are located in the colon, we evaluated the integrity of colonic mucosa of depression models via two tight junction associated proteins. DEX significantly decreased the occludin expression in the colon, which may weaken the gut mucosal barrier and facilitate the translocation of gut microbiota. DEX was orally administered in the DEX group; it may cause severe immune inhibition in gastrointestinal tract of mice, which contributed to the low level of local IL-1β in the colon tissue and impaired colon mucosa regeneration in the DEX group. A rise in serum CORT level is a risk factor for major depression or anxiety disorder.31 Compared with the two other groups, the mice treated with DEX had the highest serum level of CORT and IL-1β, which may result from DEX-induced exacerbation of pathobiont infection and bacterial endotoxins with a reduced phagocytic capacity of blood leucocytes. The level of IL-1β was significantly higher in the colon of the CUMS group than that of the control group, which may due to the high relative abundance of opportunistic pathogens.
Considering that high levels of GCs strongly inhibit hippocampal plasticity and functions,32 we further investigated the changes of neural biomarkers in the hippocampus in each group. There was no significant difference in levels of NeuN. We observed that synapsin-1 expression was significantly decreased in the CUMS and DEX groups. Therefore, we conclude that the synaptic plasticity of hippocampus is more prone to being damaged in depression models. The low expression of GFAP and Nestin in DEX and CUMS groups suggest the dysfunction of astrocytes and neural stem/progenitor cells (NSPCs) in the depression-like behaviour mice.33 34 Astroglia cells provide neurotrophic factors for the survival of neuron and proliferation of NSPCs during the course of neurogenesis35 and maintain the balance of neurons energy metabolism via lactate shuttle.36 Dysfunctions of astrocytes and synaptic plasticity in the hippocampus are common characteristics in the two different depression models.
Neuroinflammation contributes to the pathophysiology of MDD.37 P2Y12 purinergic receptor expression is necessary for neuron–microglia interactions; its receptors on microglia are thought to be a key player in the surveillance of the local environment and neurogenesis.38 The decrease in P2Y12 receptors may aggravate the damage of hippocampal neurogenesis in the CUMS and DEX groups compared with the control group. We also found a significant increase in NLRP3 in the CUMS group and a remarkable decrease of it in the DEX group. NLRP3 inflammasome plays a key role in shaping immune/inflammatory responses in brain diseases.39 These results suggest that inflammasome plays an important role in the hippocampal neuroinflammation in the CUMS group. However, we found a significant decrease in levels of caspase-1 and NLRP3 in the mice from the DEX group compared with CON mice. Since DEX is one of the most commonly used anti-inflammatory agents, the microglia activation and increased hippocampal IL-1β level in the DEX group may have originated from the damage of the neuron–microglia crosstalk in the brain and or elevated peripheral cytokines. Our results suggest the heterogeneous constitution of neuroinflammation in the formation of depression.
Although we revealed the differences in the microbiota–gut–brain axis in two common models of depression, we only focused on the neural markers in the hippocampus. It will be interesting to know if there are changes of prefrontal cortex to understand the precise mechanism of neuroinflammation in our models. The gender differences in the microbiota–gut–brain in our models also deserve additional investigation.
Up to now, studies on depression in animal models mainly involve CUMS, which simulates the pathophysiology of MDD. Repeated administration of corticosterone has been widely used to simulate behavioural changes in patients with MDD, including anhedonia, immobility and anxiety behaviour.40 Based on our findings, chronic low-grade inflammation in the peripheral and central nervous system with neural synapses loss is a common contributor to the formation of two different depression models. The increase in gut conditional pathogens with nutrient deficiency may play an important role in the formation of CUMS-induced depression. Excessive immune inhibition, obesity and the increase of gut pathobionts are likely to be key triggers of DEX-induced depression.
Data availability statement
All data relevant to the study are included in the article or uploaded as supplemental information.
Patient consent for publication
Ethics approval Obtained (20181112). The study was approved by the Institutional Ethics Committee of the Shanghai Jiao Tong University animal experiment Center (IORG Number: 20181112). The Institutional Ethics Committee of the Shanghai Jiao Tong University also approved the consent procedure.
Dr Zejian Wang obtained a PHD degree from Shanghai Jiao Tong University, China in 2003. He is currently an associate professor at Shanghai Jiao Tong University, Shanghai, China. He has participated in and presided over two National Natural Science Foundation of China. He is also a member of Shanghai Pharmacological society in China. He has attended and presented in various national and international conferences. His main research interests include pathogenesis of depression, Alzheimer's Disease and cerebral neural activity.
Contributors ZW designed the entire study, monitored data collection for the whole experiment, wrote the statistical analysis plan and analysed the data, and drafted and revised the paper. He is a guarantor. JW implemented the experiments and analysed the data, and drafted and revised the paper. XL revised the paper and thoroughly provided technical supports for the experimental design. JC initiated the collaborative project, co-supervised data collection for the whole trial and revised the paper. CG, UM, YC and JL provided the support to edit the manuscript and design the trial. All authors contributed to the manuscript and read and approved the final manuscript.
Funding This project was supported by the Shanghai Key Laboratory of Psychotic Disorders (NO.13dz2260500), National Natural Science Foundation of China (No. 81871122) and Interdisciplinary Program of Shanghai Jiao Tong University (ZH2018QNA59).
Competing interests None declared.
Provenance and peer review Commissioned; externally peer reviewed.
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.