Metabolomic signatures and microbial community profiling of depressive rat model induced by adrenocorticotrophic hormone
Journal of Translational Medicine volume 17, Article number: 224 (2019)
Adrenocorticotrophic hormone (ACTH)-treatment rat model has been utilized as a widely accepted model of treatment-resistant depression. Metabolomic signatures represent the pathophysiological phenotype of diseases. Recent studies in gut microbiota and metabolomics analysis revealed the dramatic role of microbiome in psychoneurological system diseases, but still, the mechanisms underlying gut microbiome–host interaction remain unclear.
Male Wistar rats were s.c. injection of ACTH fragment 1–24 for 14 days to induce treatment-resistant depression. Depression-related behavioral tests, analysis of serum monoamine neurotransmitters and hypothalamic–pituitary–adrenal (HPA) axis-related hormones were determined for assessment of ACTH-induced depression rat model. A gas chromatography-time-of-flight mass spectrometer based urinary metabolomic signatures integrated 16S rRNA sequence analysis based gut microbial profiling was performed, as well as Spearman’s correlation coefficient analysis was used to manifest the covariation between the differential urinary metabolites and gut microbiota of genus level.
Chronic injection of ACTH-induced depression-like phenotype (increased immobility time in forced swimming test and tail suspension test) was accompanied by peripheral serotonin down-regulation and HPA axis overactivation (ACTH and corticosterone up-regulation). Urinary metabolomics analysis indicated that pyruvic acid, l-threonine, mannitol, d-gluconic acid, 4-hydroxybenzoic acid, d-arabitol, myo-inositol and ascorbic acid levels were reduced in ACTH-treated rats’ urine, while hippurate level was elevated. In addition, microbial community profiling revealed bacterial enrichment (e.g. Ruminococcus, Klebsiella) and reduction (e.g. Akkermansia, Lactobacillus) in the ACTH-induced depression rat model. Correlation analysis showed that Akkermansia and Lactobacillus were closely relevant to metabolites myo-inositol and hippurate, which were included in host inositol phosphate metabolism, and phenylalanine, tyrosine and tryptophan biosynthesis.
Depression rat model induced by ACTH is associated with disturbance of pyruvate metabolism, ascorbate and aldarate metabolism, inositol phosphate metabolism, glycine, serine and threonine metabolism, and glycolysis or gluconeogenesis, as well as changes in microbial community structure. Gut microbiota may participate in the mediation of systemic metabolomic changes in ACTH-induced depression model. Therefore, integrated metabolomic signatures and gut microbial community profiling would provide a basis for further studies on the pathogenesis of depression.
As the latest World Health Organization (WHO) report showed, depression is a high-risk psychiatric disease that affects people of all races around the world, and ranks among the world’s largest contributors of years lived with disability. It is also one of the major causes of the global economic burden of disease . However, up to 60% of patients show no respond to the initial antidepressant treatment, which is called treatment resistant depression (TRD), which may be a subtype of depression patients with unique pathophysiological features [2, 3]. Due to lack of sufficient relief from current antidepressant therapy, TRD patients have an increasing risk of suffering from chronic psychosocial disorders, relapses and suicide .
The abnormality of hypothalamic–pituitary–adrenal (HPA) axis is one of the most concerned research areas in depression [5,6,7]. Stress can activate the HPA axis, including increased release of corticotropin-releasing hormone (CRH) and adrenocorticotropic hormone (ACTH) [8, 9]. In clinical studies, the dysregulation of HPA axis is also a common feature found in depressive patients . Numerous CRH receptors are resident in extra-hypothalamic area, and most cells of the 5-hydroxytryptamine (5-HT) and noradrenaline (NA) systems are also included herein . CRH can activate 5-HT neurons in the dorsal rape nucleus, and regulate NA activity. Glucocorticoids can stimulate 5-HT and NA systems via mineralocorticoid and glucocorticoid receptors. Therefore, the disorder of HPA axis and its interaction with central monoaminergic system might be one of the main causes of TRD [12, 13].
As an interesting animal model of TRD, chronic administration of ACTH to rodent, mimicking the stress induced the HPA axis activation, results in resistance to the treatment of imipramine in the forced swimming test (FST) , as well as resistance to other antidepressants [15, 16]. Though, ACTH-induced depression rat model has been widely accepted as an animal model for TRD, the underlying mechanism of this model is not completely understood.
Deregulation of the gut microbial community has been extensively associated with the diseases of the central nervous system (CNS). The crosstalk between the gut microbiota and the CNS, so-called gut microbiota–brain axis, occurs through neuroendocrine, enteric, autonomic, and immune system pathways, and may impact the host physiological and pathological mechanisms, including the HPA axis activation, and the communications between neurotransmitters and immune system . Therefore, the perturbed gut microbiota could be associated with neurological diseases via the gut microbiota–brain axis . The alterations of specific species of gut microbiota may contribute to the occurrence and the development of depression, while depressive states may cause the changes of specific species of gut microbiota, which eventually result in more severe depression. These specific species of gut microbiota are mainly included in the bacteria phyla of Firmicutes, Actinobacteria, Bacteroidetes and Proteobacteria . Meanwhile, a number of studies have shown that probiotics, such as Lactobacillus and Bifidobacterium, have beneficial psychoneurological effects [20, 21]. It is reported that CBM588, a specific phenotype of the strain Clostridium butyricum, in combination with antidepressants is effective and well tolerated in the treatment of TRD . Recently, a 16S ribosomal RNA gene sequencing of feces based study showed that the antidepressant effects of ketamine, an N-methyl-d-aspartate receptor (NMDAR) antagonist showed rapid and sustained antidepressant effects in TRD patients, might be partly mediated by the restoration of altered gut microbiota .
The crosstalk between gut microbiota and brain might be indirect. The metabolites of gut microbiota, which circulate in the systematic circulation and metabolism, may play a vital role in the gut microbiota–brain axis. Metabolomics is a systematic biological method for measuring the changes of endogenous small molecule metabolites. It has been increasingly applied to the diagnosis of diseases, the discovery of disease biomarkers, and the pathogenesis study of diseases . In metabolomics studies, urinary sample is non-invasive, and more convenient for collection compared to other samples (blood or specific organ). Furthermore, the greatest advantage of urinary sample for metabolomics studies is that it contains a combination of metabolites in a relatively longer window period (12 or 24 h), while the metabolites in blood can only represent the transient state of the research object. At present, urinary samples of animal or human have been widely used for metabolomics research to reveal the pathophysiological mechanism of diseases systemically . Therefore, to characterize the global characteristics of ACTH-induced TRD, it is necessary to combine the gut microbiota study with the metabolomic signatures.
In the present study, an integrative gas chromatography-time-of-flight mass spectrometer (GC-TOFMS)-based urinary metabolomics and 16S rRNA gene sequencing approach was performed on chronic ACTH-treated rat model with depression phenotype to reveal the differentiated gut microbiota and urinary metabolites. Further, the relevancies of above two were discovered to reveal the interactions between the gut microbiota and the host metabolism in ACTH-induced TRD rats. The results will help to enhance comprehensive understanding of the pathogenesis of ACTH-induced TRD.
Materials and methods
Animals and treatments
A total of 20 male Wistar rats (10 weeks, 230–250 g) were obtained from Shanghai Sippr-BK Laboratory Animal Co. Ltd., China. The animals were under a constant light–dark cycle (07:00–19:00 at 40 W light condition), with free access to food and water ad libitum in a standard laboratory environment (20 ± 5 °C and 55 ± 15% humidity). All the study was conducted under the Guidelines for Animal Experimentation of Shanghai University of Traditional Chinese Medicine.
The rats were randomly divided into two groups: normal control group (Control, n = 10) and ACTH-treated TRD model group (ACTH, n = 10). ACTH group received s.c. injection of adrenocorticotropic hormone fragment 1–24 (100 μg/rat, 14 days, CAS No. 16960-16-0, purchased from Sigma Chemical Co., MO, USA) , and at the same time, Control group received the same volume of saline injection.
Forced swimming test
Forced swimming test (FST) , a representative behavioral test for depression, is frequently used to evaluate the depressive state in rodent models. Forced immersion of rats in water for an extended period produces a characteristic behavior of immobility, and rats of depression are more inclined to exhibit a longer immobility time in a test period. After 14 days of ACTH-treatment, each rat was placed in a transparent glass container (50 cm height × 18 cm diameter) filled to 30 cm with 25 ± 1 °C water. Rats were forced to swim for 6 min and the duration of immobility within the last 4 min was measured. All behavioral tests were performed between 14:00 and 18:00 p.m. Immobility was defined as floating in the water without struggling and only keeping the head above the water without other motions. The immobility time was determined by an observer blinded to the purpose of the experiment and the water must be changed after every test.
Tail suspension test
Tail suspension test (TST)  is based on the fact that animals subjected to the short-term, inescapable stress of being suspended by their tail, will develop an immobile posture. Each rat was suspended 50 cm above the floor by tape. Every test was performed for 6 min and the duration of immobility within the last 4 min was measured. Immobility means that the animal gives up its active struggle and has not moved at all. Immobility time was determined by an observer who was unaware the purpose of the test.
Sample collection and preparation
After the last behavior test, all rats were fasted and placed in metabolic cages individually 24 h for urinary sample collection. The urinary samples were centrifuged at 1500×g for 5 min at 4 °C. The resultant supernatants were stored at − 80 °C pending GC-TOFMS analysis. As for fecal samples, approximately 2 pellets of feces were attained from each rat, then transferred into EP tubes, and stored at − 80 °C until further analysis. Blood was collected from the abdominal aorta after anesthesia by chloral hydrate (300 mg/kg body weight, i.p.). Blood samples were centrifuged at 3000×g for 15 min at 4 °C, and the supernatants were stored at − 80 °C.
Serum parameters analysis
Serum samples were conducted for monoamine neurotransmitter and hormone analysis with ELISA kits, such as 5-HT, norepinephrine (NE), CRH, ACTH, and corticosterone (CORT), according to the manufacturer’s protocols (Cusabio Biotechnology Co., Ltd, Wuhan, China).
Urinary metabolomic signatures
According to our previous study , metabolomics analysis of urinary samples was based on a GC-TOFMS metabolomics method. Rat urinary samples were thawed at room temperature 2 h before the experiment and 2.0 mL urinary samples were taken in an EP tube and centrifuged at 12,000 rpm for 10 min. 100 μL of supernatant was subjected to urea degradation by adding urease (70 IU, 37 °C, 15 min), and then pretreated with l-2-chlorophenyl alanine as an internal standard and ethyl chloroformate as a derivatizing agent. The sample analysis was performed using GC-TOFMS (Pegasus HT, Leco Corp., St. Joseph, MI; electron ionization mode).
Microbial community profiling
To analyze the taxonomic composition of the microbial community, 16S rRNA sequencing of fecal sample was performed following the method we published previously with minor modifications . Total DNA in fecal samples was extracted using E.Z.N.A.® soil DNA Kit (Omega Bio-Tek, Norcross, GA, USA). The quality of DNA extraction was detected by 0.8% agarose gel electrophoresis and DNA was quantified using an ultraviolet spectrophotometer. The V3–V4 region of the bacteria 16S rRNA gene was amplified by PCR using Q5 High-Fidelity DNA polymerase (NEB, USA) and primers 520 F (5′-AYTGGGYDTAAAGNG-3′) and 802 R (5′-TACNVGGGTATCTAATCC-3′). The amplification products of PCR were qualified by 2% agarose gel electrophoresis, and purified by AxyPrep DNA Gel Extraction Kit (AXYGEN Biosciences, Union City, CA, USA). Then the amplification products were quantified on Microplate reader (BioTek, FLx800) using Quant-iT PicoGreen dsDNA Assay Kit (Invitrogen, USA). The purified amplicons were finally combined at equimolar concentrations and paired-end sequencing was performed on an Illumina MiSeq platform (Majorbio Bio-Pharm Technology Co. Ltd., Shanghai, China).
The raw data from GC-TOFMS analysis were exported in NetCDF format using the Chroma TOF software (v4.44, Leco Co., Los Angeles, CA), and further preprocessed by R 2.13.2 (Lucent Technologies). Then the peak area, retention time and compounds name were included in output data for multivariate analysis using SIMCA-P 11 software (Umetrics, Umea, Sweden) to perform principal component analysis (PCA) and orthogonal partial least squares discriminant analysis (OPLS-DA). Furthermore, OPLS-DA was used to select the differential metabolites with variable importance in the projection (VIP) value > 1.0. The significance identification was then performed by Student’s t test (P < 0.05) and fold change (FC, compared to Control group, a range of > 1.2 or < 0.8). Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.genome.jp/kegg/) and Human Metabolome Database (HMDB; http://www.hmdb.ca/) were applied to confirmed potential depression biomarkers. The metabolic pathway analysis module of MetaboAnalyst 3.0 (http://www.metaboanalyst.ca/) was used to identify the metabolic pathways of differential metabolites.
For sequencing data, Trimmomatic, a flexible read trimming tool for Illumina sequence data, was used to quality-filter the raw fastq files. The reads were truncated in a sliding window of 10 bp with an average quality score > 20. The paired-end sequences passed through the quality-filter were merged using FLASH (v1.2.7, http://ccb.jhu.edu/software/FLASH/). Bases mismatch was not allowed. The sequence of each sample was isolated based on the barcode (exactly matching) and the primer (exactly matching); readings containing ambiguous bases were discarded. Operational Taxonomic Units (OTUs) were clustered with 97% similarity cut-off using UPARSE (v7.1 http://drive5.com/uparse/). The taxonomy of each 16S rRNA gene sequence was classified against the Greengenes database using a confidence threshold of 70%. Partial least squares discriminant analysis (PLS-DA) was used to find out whether ACTH-treated rats could be separated from the control rats. Furtherly, the linear discriminant analysis (LDA) effect size (LEfSe) method was used to identify the dominant phylotypes responsible for differences between ACTH group and Control group. LEfSe aims to provide biomarkers that explain the effects of biomarker discovery comparisons and hypothesis-driven researches on phenotypes of interest (two or more). The visualization of the biomarkers found on taxonomic trees provides an effective tool of concluding the results in a biologically meaningful manner, because it captures statistically and visually the hierarchical relationships inherent in 16S-based taxonomies/phylogenies or in pathways and biomolecular functions.
The data was analyzed using SPSS 21.0. Statistical significance among the groups was conducted by two-tailed Student’s t-test. All values were expressed as the mean ± standard deviation (SD). Probability values of less than 0.05 were considered to show a statistical significance.
Correlation analysis between metabolomic signatures and microbial community profiling
To explore the relationship between altered gut microbiota and perturbed urinary metabolites, correlation coefficient of Spearman’s between metabolomic signatures and gut microbial alterations of genus level were exhibited in the form of a heat map. Coefficient r values were distinguished by different color blocks. Additionally, the metabolic associations of each well-correlated members of gut microbe (|r| > 0.3) were visualized as a cross-correlation diagram showing positive (in red line) and negative (in blue line) connections.
Quality assessment of ACTH-induced depression model
Combined FST and TST results were used to access the depressive state of ACTH model. As shown in Additional file 1: Fig. S1, compared to Control group, immobility time of rats treated with ACTH in FST and TST was significantly increased (P < 0.01).
NE and 5-HT are monoamine neurotransmitters closely associated with depression. Serum 5-HT and NE levels of two groups were shown in Additional file 1: Fig. S2. As the results of s.c. injection of ACTH, serum concentration of 5-HT was significantly decreased (P < 0.01) in ACTH group, compared with Control group. At the same time, the level of NE in ACTH group showed a slight decrease trend with no significance.
CRH, ACTH and CORT are the major hormones in the HPA axis. In our study, ACTH and CORT levels were markedly increased in the serum of ACTH-treated rats compared to Control group (P < 0.001). The CRH level of ACTH group was lower than that in Control group, however there was no statistical difference (P > 0.05) (Additional file 1: Fig. S2). This may be the result of negative feedback from ACTH injection.
Urinary metabolomic signatures of ACTH-treated rats
Urinary metabolomics is a typical metabonomic research approach, which is characterized by convenient collection, non-injury, and continuity of metabolites. As shown in Fig. 1a, the metabolomic signatures of ACTH-treated rats were remarkably separated from Control group in PCA (R2X = 0.540, Q2 = 0.43) and OPLS-DA score plots (R2X = 0.646, R2Y = 0.988, Q2 = 0.955, Additional file 1: Fig. S3), suggesting that ACTH injection caused significant metabolic changes. In addition, the variables with VIP > 1 (surrounded by red boxes) in V-plot (Fig. 1b) were further analyzed. We finally identified 9 metabolites. Compared with Control group, the levels of pyruvic acid, l-threonine, mannitol, d-gluconic acid, 4-hydroxybenzoic acid, d-arabitol, myo-inositol, ascorbic acid were lower in ACTH group, and the level of hippurate was increased (Fig. 1c, Additional file 2: Table S1).
Metabolic pathway and network analysis
KEGG and HMDB were used to link the 9 urinary differential metabolites to potential relevant pathway, and MetaboAnalyst 3.0 was applied to identify the impact value of relevant pathways. The results showed that five main metabolic pathways were affected, including pyruvate metabolism, ascorbate and aldarate metabolism, inositol phosphate metabolism, glycine, serine and threonine metabolism, and glycolysis or gluconeogenesis (Fig. 2, impact factor ≥ 0.1). Based on the relationship between metabolomic signatures, the disturbed metabolic pathways were represented in diagram (Fig. 3).
Microbial community profiling of ACTH-induced depression
In order to analyze whether there was a difference in the structure of the microbial community in the gut microbiota between ACTH group and Control group, PLS-DA analysis was performed as shown in Fig. 4a. The results showed ACTH group had a distinct microbiota composition that clustered differently from Control groups. Moreover, the microbial community structures at the phylum level of ACTH group and Control group were shown in Fig. 4b. The abundance of each sample is presented as the percentage of the total number of sequences. In ACTH-treated rats, the relative abundances of Firmicutes (P < 0.05) and Bacteroidetes (P < 0.01) were increased significantly, and the relative abundance of Verrucomicrobia (P < 0.01) was decreased significantly. In addition, the ratio of Firmicutes-to-Bacteroidetes in ACTH group was higher than that of Control group (Additional file 1: Fig. S4). At the genus level (Additional file 1: Fig. S5), gut microbiota of two groups were both dominated by Ruminococcus, Klebsiella, Lactobacillus, and Akkermansia. The relative abundances of Ruminococcus (P < 0.05) and Klebsiella (P < 0.01) in ACTH group were significantly higher than those in Control group, but the relative abundances of Lactobacillus (P < 0.01) and Akkermansia (P < 0.01) were decreased significantly. The metagenome analysis of LEfSe approach was applied to explore the differential gut microbiota (from phylum to genus) in two groups of rats (Fig. 4c, LDA score > 3.0). Anaeroplasma, Anaeroplasmataceae, Anaeroplasmatales, Rumincoccus and Oscillospira were the most abundant differential microbiota in ACTH group. In Control group, Lactobacillus, Lactobacillaceae, Coprococcus, Actinobacteria, Actinomycetales, Lachnobacterium and Burkholderiales were significantly increased compared to ACTH group.
Correlation analysis between metabolomic signatures and microbial community profiling
The covariation of urinary differential metabolites and the gut microbiota of genus level was presented as a heat map diagram (Fig. 5a). The metabolic connections of well-predicted bacteria in gut microbiota were shown in Fig. 5b (|r| > 0.3). There were multiple correlations between gut microbiota of genus level and metabolites, such as negative correlation between Ruminococcus and pyruvic acid (r = − 0.338), mannitol (r = − 0.313), and positive correlation with hippurate (r = 0.313). Akkermansia has a positive correlation with pyruvic acid (r = 0.5), l-threonine (r = 0.335), 4-hydroxybenzoic acid (r = 0.484), myo-inositol (r = 0.637, P < 0.05), d-arabitol (r = 0.370) and ascorbic acid (r = 0.336), but negatively related to hippurate (r = − 0.415). Lactobacillus was negatively correlated with hippurate (r = − 0.601, P < 0.05), but positively correlated with pyruvic acid (r = 0.334) and d-arabitol (r = 0.695, P < 0.05). In contrast, Klebsiella was positively associated with hippurate (r = 0.411) and negatively correlated with l-threonine (r = − 0.340) and mannitol (r = − 0.412).
Depression is a multifactorial disease which is regulated by genetic and environmental factors. As a systematic disease, it is a complex phenotype associated abnormalities ranging from central nervous system to many peripheral systems [30,31,32]. As an important component of the neuroendocrine system, HPA axis is a complex collection of direct influences and feedback interactions, and controls the response to stress and regulates many body processes including digestion, the immune system, mood and emotion, sexuality, energy storage and expenditure [33, 34]. In humans, gut microbiota contains the largest number of bacteria and the greatest number of species compared to other areas of the body. It not only plays a barrier to pathogenic organisms, but also serves as an endocrine organ by providing short-chain fatty acids (SCFAs), vitamin, and xenobiotics to the host . Although the vital role of HPA axis dysfunction or gut microbiota dysregulation in pathogenesis of TRD has been widely accept, the interactions between the structural features of gut microbiota and host metabolic phenotype in depression accompanied by HPA axis dysfunction have rarely been reported until now. In the present study, an integrative metabolomic signatures and microbial community profiling was employed on ACTH-induced depression rats compared with normal control rats.
For FST and TST (Additional file 1: Fig. S1), striking helpless and depression-like behaviors observed in ACTH group were consistent with other reports [25, 36]. Elevated levels of ACTH and CORT in serum and negative feedback of CRH level (Additional file 1: Fig. S2) confirmed that ACTH injection induced hyperactivity of HPA axis, which was reported to be closely related to depression [37, 38].
The metabonomic results indicated that the differential metabolites were mainly related to disorders of energy metabolism, amino acid metabolism, ascorbate and aldarate metabolism, and inositol phosphate metabolism (Fig. 3). The disorder of energy metabolism was mainly related to pyruvate metabolism and glycolysis or gluconeogenesis. Pyruvic acid is an important intermediate in the glucose metabolism and its carboxylate anion is called pyruvate. Pyruvate is the end product of glycolysis and can be converted to acetyl coenzyme A (acetyl-CoA) by oxidative decarboxylation, which further participates in tricarboxylic acid (TCA) cycle to generate ATP . It was reported that perturbed energy metabolism is ubiquitous in animal models of depression [40, 41]. In ACTH group, pyruvic acid concentration was significantly reduced (Additional file 2: Table S1), interfering with glycolysis or gluconeogenesis pathway, which may further lead to TCA circulatory dysfunction and depressive symptoms. Threonine, the essential dietary amino acid, is the precursor to isoleucine. Since Threonine is mainly found in the central nervous system, it can be of great help in the treatment of different types of depression. Threonine contributes to the synthesis of glycine and serine, and plays a key role in maintaining the integrity and barrier function of the intestinal mucosal . Previous studies showed that abnormal amino acid metabolism might usually happen in depression and disturbed amino acid level might be served as clinical trait-biomarker for depression [43, 44]. Our research found that l-threonine level was significantly reduced in ACTH group and played a key role in the glycine, serine and threonine metabolic pathways, which may have great impact in the occurrence of depression (Additional file 2: Table S1). As the major metabolite of human, hippurate is formed by the reaction of many aromatic compounds such as benzoic acid and toluene with glycine and is involved in phenylalanine, tyrosine and tryptophan biosynthesis  (Fig. 3). Glycine is one of the excitatory amino acids closely related to major depression . The increase of hippurate concentration is closely related to the change of glycine concentration. It has been reported that gut microbiota can signal to the brain through the regulation of tryptophan metabolism, thus affecting the development, function and behavior of the brain, leading to the occurrence of neurodevelopmental diseases derived from gut microbiota . In our study, the hippurate concentration in ACTH group was significantly higher than that in Control group (Additional file 2: Table S1), and a critical signal molecule in the brain-gut axis, 5-HT, was significantly decreased (Additional file 1: Fig. S2), possibly leading to the occurrence of brain-gut axis disorders and depression.
In addition, ascorbate and aldarate metabolism was also disrupted in ACTH-induced depression rats compared to Control group. Ascorbic acid, an important antioxidant in the CNS, promotes the synthesis of neurotransmitters (5-HT and NE) . Compared with Control group, the level of ascorbic acid in ACTH group was significantly reduced (Additional file 2: Table S1), resulting in lower levels of 5-HT and NE in serum (Additional file 1: Fig. S2). The abnormality of these monoamines is widely accepted as the main biological mechanism of depression, and targeted by most of antidepressants.
Myo-inositol, a stereoisomer of inositol, is a carbocyclic sugar, which is widely present in brain and other mammalian tissues. It mediates cell signal transduction through various hormones, neurotransmitters and growth factors and is involved in osmoregulation. Myo-inositol has the highest concentration in the brain and plays an important role here, combining other neurotransmitters and some steroid hormones with their receptors . It’s also an essential participant and regulator of inositol phosphate metabolism . Chiappelli et al.  found that the occurrence of depression in schizophrenic patients was closely related to the reduction of myo-inositol levels by proton magnetic resonance spectroscopy analysis. This finding is consistent with the decrease in myo-inositol levels in ACTH group (Additional file 2: Table S1).
In the microbial community profiling, the PLS-DA analysis and the compositional structure both presented the significant differentiation between the groups (Fig. 4a, b). According to the reports, the increase of the ratio of Firmicutes to Bacteroidetes of phylum level was commonly observed in both depression patients  and patients with irritable bowel syndrome . Therefore, the increased ratio of Firmicutes to Bacteroidetes is considered to be a marker of the gut–brain axis disorders, which also occurred in ACTH group (Additional file 1: Fig. S4).
The LEfSe difference analysis showed that the differential microbiota was mainly concentrated on the cumulative number of Oscillospira and Ruminococcus (Fig. 4c). Oscillospira is a mysterious bacterial genus that has never been cultured but is continuously detected by sequencing of 16S rRNA gene in the human microbiome . In terms of the metabolites excreted by Oscillospira, it is inversely correlated with inflammatory diseases and BMI . BMI represented the best predictive value for gut dysbiosis and metabolic alterations. It is speculated that Oscillospira may produce short-chain fatty acid butyrate, similar to trends of other genera (e.g., Roseburia and Faecalibacterium) . Butyrate concentrations were negatively correlated with anxiety levels. In addition, Ruminococcus, a genus of class Clostridia which pointing to an energy metabolism dysfunction, was identified as one of characteristic genus of individuals with inflammatory bowel disease (IBD) .
By integrating gut microbiota data with differential metabolites, our results highlighted the tight crosstalk between the gut microbiota and host metabolism in ACTH-induced depression rats. As shown in Fig. 5, the relevant genus of gut microbiota predicted from each of differential urinary metabolites were Ruminococcus, Akkermansia, Lactobacillus and Klebsiella. Akkermansia muciniphila, the only currently known species within genus Akkermansia, has the ability to continuity regeneration of the intestinal epithelium and degrade mucin, and also producing acetic acid, propionic acid and oligosaccharides. These products become the substrate for Faecalibacterium prausnitzii, one of the main producers of butyrate in the intestine. Anaerobic bacteria that produce butyrate help to inhibit inflammation in the gastrointestinal tract, and prevent an increase in intestinal permeability . In our study, 6 metabolites were statistically relevant to Akkermansia, including pyruvic acid, 4-hydroxybenzoic acid, myo-inositol, hippurate, d-arabitol and ascorbic acid, many of which are closely related to depression, as we discussed above.
Lactobacillus is a genus of Gram-positive, facultative anaerobic or microaerophilic, is the main component the lactic acid bacteria group (i.e. they convert sugars to lactic acid). It was reported ingestion of Lactobacillus strain could regulate emotional behavior and central GABA receptor expression in mouse via the vagus nerve , as well as improve behavioral, cognitive, and biochemical aberrations caused by chronic restraint stress . As shown in Fig. 5, there are 3 metabolites relevant to Lactobacillus, namely pyruvic acid, hippurate and d-arabitol. Arabitol is a sugar alcohol, which can be formed by reducing arabinose or lyxose. Some organic acid tests examined the presence of d-arabitol, which may indicate overgrowth of intestinal microbes. Arabitol showed very similar inhibition effects to its isomer xylitol on oral Lactobacilli . As we discussed above, the correlation analysis between differential metabolites and bacteria revealed that active metabolites of Akkermansia and Lactobacillus were closely related to the host inositol phosphate metabolism, phenylalanine, tyrosine and tryptophan biosynthesis, which indicated there was a significant crosstalk between the gut microbiota and the host metabolism in ACTH-induced TRD rats.
In addition, it has been reported that the abundance of Klebsiella and Ruminococcus was changed in patients with MDD compared with healthy controls [62, 63]. In our study, these two intestinal microbes were observed to have varying degrees of correlation with l-threonine, pyruvic acid, mannitol and hippurate, which are closely related to the development of depression as we discussed earlier, but there was no statistical significance (P > 0.05). These results suggest that the changes in gut microbiota and urinary metabolic phenotype are associated with the development of depression. The relevance between the gut microbiota of genus level and the differential urinary metabolites further confirmed the correlation between disturbances of gut microbiota and the disturbed microorganism-host metabolic balance in ACTH-treated rats.
In this study, 16S rRNA gene sequencing was combined with GC-TOFMS based metabolomics to assess effects of ACTH-treatment on gut microbial community profiling and urinary metabolomic signatures. Our results showed that ACTH-induced depression not only disturbs the composition of gut microbiota, such as Oscillospira, Ruminococcus, Akkermansia, Lactobacillus and Klebsiella, but also alters the homeostasis of metabolic pathway that including pyruvate metabolism, ascorbate and aldarate metabolism, inositol phosphate metabolism, glycine, serine and threonine metabolism, and glycolysis or gluconeogenesis. In addition, the correlation analysis revealed that altered gut microbes were closely related to differential urinary metabolites. These gut microbiota with strong host metabolic connectivity can provide a basis for the disorder of gut microbiota community in ACTH-induced depression.
Availability of data and materials
The data are included in the article as figures and its additional files.
gas chromatography-time-of-flight mass spectrometer
World Health Organization
treatment resistant depression
forced swimming test
central nervous system
tail suspension test
principal component analysis
orthogonal partial least squares discriminant analysis
variable importance in the projection
Kyoto Encyclopedia of Genes and Genomes
Human Metabolome Database
operational taxonomic units
partial least squares discriminant analysis
linear discriminant analysis
linear discriminant analysis effect size
short-chain fatty acids
inflammatory bowel disease
World Health Organization. Depression and other common mental disorders. 2017; WHO/MSD/MER/2017.2.
Kola I, Landis J. Can the pharmaceutical industry reduce attrition rates? Nat Rev Drug Discov. 2004;3:711–5.
Srikumar BN, Paschapur M, Kalidindi N, Adepu B, Das ML, Sreedhara MV, et al. Characterization of the adrenocorticotrophic hormone-induced mouse model of resistance to antidepressant drug treatment. Pharmacol Biochem Behav. 2017;161:53–61.
Miyaoka T, Kanayama M, Wake R, Hashioka S, Hayashida M, Nagahama M, et al. Clostridium butyricum MIYAIRI 588 as adjunctive therapy for treatment-resistant major depressive disorder. Clin Neuropharmacol. 2018;41:151–5.
Wang S, Wang C, Yu Z, Wu C, Peng D, Liu X, et al. Agarwood essential oil ameliorates restrain stress-induced anxiety and depression by inhibiting HPA axis hyperactivity. Int J Mol Sci. 2018;19:3468.
Zhao J, Niu C, Wang J, Yang H, Du Y, Wei L, et al. The depressive-like behaviors of chronic unpredictable mild stress-treated mice, ameliorated by Tibetan medicine Zuotai: involvement in the hypothalamic–pituitary–adrenal (HPA) axis pathway. Neuropsychiatr Dis Treat. 2018;14:129–41.
Hamada M, Nishigawa T, Maesono S, Aso K, Ikeda H, Furuse M. Decreased stress-induced depression-like behavior in lactating rats is associated with changes in the hypothalamic–pituitary–adrenal axis, brain monoamines, and brain amino acid metabolism. Stress. 2019;2019:1–10.
Nicolaides NC, Kyratzi E, Lamprokostopoulou A, Chrousos GP, Charmandari E. Stress, the stress system and the role of glucocorticoids. Neuroimmunomodulation. 2015;22:6–19.
Clow A, Hucklebridge F, Stalder T, Evans P, Thorn L. The cortisol awakening response: more than a measure of HPA axis function. Neurosci Biobehav Rev. 2010;35:97–103.
Farrell C, Doolin K, O’Leary N, Jairaj C, Roddy D, Tozzi L, et al. DNA methylation differences at the glucocorticoid receptor gene in depression are related to functional alterations in hypothalamic–pituitary–adrenal axis activity and to early life emotional abuse. Psychiatry Res. 2018;265:341–8.
Su J, Tanaka Y, Muratsubaki T, Kano M, Kanazawa M, Fukudo S. Injection of corticotropin-releasing hormone into the amygdala aggravates visceral nociception and induces noradrenaline release in rats. Neurogastroenterol Motil. 2015;27:30–9.
Fischer S, Macare C, Cleare AJ. Hypothalamic–pituitary–adrenal (HPA) axis functioning as predictor of antidepressant response-meta-analysis. Neurosci Biobehav Rev. 2017;83:200–11.
Joels M, Baram TZ. The neuro-symphony of stress. Nat Rev Neurosci. 2009;10:459–66.
Kitamura Y, Araki H, Gomita Y. Influence of ACTH on the effects of imipramine, desipramine and lithium on duration of immobility of rats in the forced swim test. Pharmacol Biochem Behav. 2002;71:63–9.
Iwai T, Ohnuki T, Sasaki-Hamada S, Saitoh A, Sugiyama A, Oka J. Glucagon-like peptide-2 but not imipramine exhibits antidepressant-like effects in ACTH-treated mice. Behav Brain Res. 2013;243:153–7.
Walker AJ, Burnett SA, Hasebe K, McGillivray JA, Gray LJ, McGee SL, et al. Chronic adrenocorticotrophic hormone treatment alters tricyclic antidepressant efficacy and prefrontal monoamine tissue levels. Behav Brain Res. 2013;242:76–83.
Rieder R, Wisniewski PJ, Alderman BL, Campbell SC. Microbes and mental health: a review. Brain Behav Immun. 2017;66:9–17.
Liang S, Wu X, Hu X, Wang T, Jin F. Recognizing depression from the microbiota–gut–brain axis. Int J Mol Sci. 2018;19:1592. https://doi.org/10.3390/ijms19061592.
Chen JJ, Zheng P, Liu YY, Zhong XG, Wang HY, Guo YJ, et al. Sex differences in gut microbiota in patients with major depressive disorder. Neuropsychiatr Dis Treat. 2018;14:647–55.
Messaoudi M, Violle N, Bisson JF, Desor D, Javelot H, Rougeot C. Beneficial psychological effects of a probiotic formulation (Lactobacillus helveticus R0052 and Bifidobacterium longum R0175) in healthy human volunteers. Gut Microbes. 2011;2:256–61.
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:755–64.
Qu Y, Yang C, Ren Q, Ma M, Dong C, et al. Comparison of (r)-ketamine and lanicemine on depression-like phenotype and abnormal composition of gut microbiota in a social defeat stress model. Sci Rep. 2017;7:15725.
McCartney A, Vignoli A, Biganzoli L, Love R, Tenori L, Luchinat C, et al. Metabolomics in breast cancer: a decade in review. Cancer Treat Rev. 2018;67:88–96.
Zhang A, Sun H, Wu X, Wang X. Urine metabolomics. Clin Chim Acta. 2012;414:65–9.
Tokita K, Fujita Y, Yamaji T, Hashimoto K. Depressive-like behavior in adrenocorticotropic hormone-treated rats blocked by memantine. Pharmacol Biochem Behav. 2012;102:329–34.
Armario A, Gavalda A, Marti O. Forced swimming test in rats: effect of desipramine administration and the period of exposure to the test on struggling behavior, swimming, immobility and defecation rate. Eur J Pharmacol. 1988;158:207–12.
Steru L, Chermat R, Thierry B, Simon P. The tail suspension test: a new method for screening antidepressants in mice. Psychopharmacology. 1985;85:367–70.
Wang X, Zhao T, Qiu Y, Su M, Jiang T, Zhou M, et al. Metabonomics approach to understanding acute and chronic stress in rat models. J Proteome Res. 2009;8:2511–8.
Zhao L, Zhang Q, Ma W, Tian F, Shen H, Zhou M. A combination of quercetin and resveratrol reduces obesity in high-fat diet-fed rats by modulation of gut microbiota. Food Funct. 2017;8:4644–56.
Dantzer R, O’Connor JC, Freund GG, Johnson RW, Kelley KW. From inflammation to sickness and depression: when the immune system subjugates the brain. Nat Rev Neurosci. 2008;9:46–56.
Dinan TG, Cryan JF. Melancholic microbes: a link between gut microbiota and depression? Neurogastroenterol Motil. 2013;25:713–9.
Albertini G, Deneyer L, Ottestad-Hansen S, Zhou Y, Ates G, Walrave L, et al. Genetic deletion of xCT attenuates peripheral and central inflammation and mitigates LPS-induced sickness and depressive-like behavior in mice. Glia. 2018;66:1845–61.
Mitani H, Shirayama Y, Yamada T, Kawahara R. Plasma levels of homovanillic acid, 5-hydroxyindoleacetic acid and cortisol, and serotonin turnover in depressed patients. Prog Neuropsychopharmacol Biol Psychiatry. 2006;30:531–4.
Fumagalli F, Molteni R, Racagni G, Riva MA. Stress during development: impact on neuroplasticity and relevance to psychopathology. Prog Neurobiol. 2007;81:197–217.
Clarke G, Stilling RM, Kennedy PJ, Stanton C, Cryan JF, Dinan TG. Minireview: gut microbiota: the neglected endocrine organ. Mol Endocrinol. 2014;28:1221–38.
Kim Y, McGee S, Czeczor JK, Walker AJ, Kale RP, Kouzani AZ, et al. Nucleus accumbens deep-brain stimulation efficacy in ACTH-pretreated rats: alterations in mitochondrial function relate to antidepressant-like effects. Transl Psychiatry. 2016;6:e842.
Schuhmacher A, Lennertz L, Wagner M, Hofels S, Pfeiffer U, Guttenthaler V, et al. A variant of the neuronal amino acid transporter SLC6A15 is associated with ACTH and cortisol responses and cognitive performance in unipolar depression. Int J Neuropsychopharmacol. 2013;16:83–90.
Cai L, Li R, Tang WJ, Meng G, Hu XY, Wu TN. Antidepressant-like effect of geniposide on chronic unpredictable mild stress-induced depressive rats by regulating the hypothalamus–pituitary–adrenal axis. Eur Neuropsychopharmacol. 2015;25:1332–41.
Kladna A, Marchlewicz M, Piechowska T, Kruk I, Aboul-Enein HY. Reactivity of pyruvic acid and its derivatives towards reactive oxygen species. Luminescene. 2015;30:1153–8.
Mallei A, Failler M, Corna S, Racagni G, Mathe AA, Popoli M. Synaptoproteomic analysis of a rat gene-environment model of depression reveals involvement of energy metabolism and cellular remodeling pathways. Int J Neuropsychopharmacol. 2014;18:pyu067.
Reus GZ, Nacif MP, Abelaira HM, Tomaz DB, Dos SM, Carlessi AS, et al. Ketamine treatment partly reverses alterations in brain derived-neurotrophic factor, oxidative stress and energy metabolism parameters induced by an animal model of depression. Curr Neurovasc Res. 2015;12:73–84.
Mao X, Zeng X, Qiao S, Wu G, Li D. Specific roles of threonine in intestinal mucosal integrity and barrier function. Front Biosci. 2011;3:1192–200.
Lu YR, Fu XY, Shi LG, Jiang Y, Wu JL, Weng XJ, et al. Decreased plasma neuroactive amino acids and increased nitric oxide levels in melancholic major depressive disorder. BMC Psychiatry. 2014;14:123.
Baranyi A, Amouzadeh-Ghadikolai O, von Lewinski D, Rothenhausler HB, Theokas S, Robier C, et al. Branched-chain amino acids as new biomarkers of major depression—a novel neurobiology of mood disorder. PLoS ONE. 2016;11:e160542.
Pero RW. Health consequences of catabolic synthesis of hippuric acid in humans. Curr Clin Pharmacol. 2010;5:67–73.
Altamura C, Maes M, Dai J, Meltzer HY. Plasma concentrations of excitatory amino acids, serine, glycine, taurine and histidine in major depression. Eur Neuropsychopharmacol. 1995;5(Suppl):71–5.
Kelly JR, Borre Y, O’Brien C, Patterson E, El Aidy S, Deane J, et al. Transferring the blues: depression-associated gut microbiota induces neurobehavioural changes in the rat. J Psychiatr Res. 2016;82:109–18.
Covarrubias-Pinto A, Acuna AI, Beltran FA, Torres-Diaz L, Castro MA. Old things new view: ascorbic acid protects the brain in neurodegenerative disorders. Int J Mol Sci. 2015;16:28194–217.
Croze ML, Soulage CO. Potential role and therapeutic interests of myo-inositol in metabolic diseases. Biochimie. 2013;95:1811–27.
Njau S, Joshi SH, Leaver AM, Vasavada M, Van Fleet J, Espinoza R, et al. Variations in myo-inositol in fronto-limbic regions and clinical response to electroconvulsive therapy in major depression. J Psychiatr Res. 2016;80:45–51.
Chiappelli J, Rowland LM, Wijtenburg SA, Muellerklein F, Tagamets M, McMahon RP, et al. Evaluation of myo-inositol as a potential biomarker for depression in schizophrenia. Neuropsychopharmacology. 2015;40:2157–64.
Zheng P, Zeng B, Zhou C, Liu M, Fang Z, Xu X, et al. Gut microbiome remodeling induces depressive-like behaviors through a pathway mediated by the host’s metabolism. Mol Psychiatry. 2016;21:786–96.
Wong ML, Inserra A, Lewis MD, Mastronardi CA, Leong L, Choo J, et al. Inflammasome signaling affects anxiety- and depressive-like behavior and gut microbiome composition. Mol Psychiatry. 2016;21:797–805.
Goodrich JK, Waters JL, Poole AC, Sutter JL, Koren O, Blekhman R, et al. Human genetics shape the gut microbiome. Cell. 2014;159:789–99.
Verdam FJ, Fuentes S, de Jonge C, Zoetendal EG, Erbil R, Greve JW, et al. Human intestinal microbiota composition is associated with local and systemic inflammation in obesity. Obesity. 2013;21:E607–15.
Konikoff T, Gophna U. Oscillospira: a central, enigmatic component of the human gut microbiota. Trends Micobiol. 2016;24:523–4.
Berry D, Reinisch W. Intestinal microbiota: a source of novel biomarkers in inflammatory bowel diseases? Best Pract Res Clin Gastroenterol. 2013;27:47–58.
Szachta P, Bartnicka A, Galecka M. Microbiota—a key to healing the gastrointestinal tract? Pomeranian J Life Sci. 2016;62:21–4.
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. 2011;108:16050–5.
Liang S, Wang T, Hu X, Luo J, Li W, Wu X, et al. Administration of Lactobacillus helveticus NS8 improves behavioral, cognitive, and biochemical aberrations caused by chronic restraint stress. Neuroscience. 2015;310:561–77.
Loman AA, Ju LK. Inhibitory effects of arabitol on caries-associated microbiologic parameters of oral Streptococci and Lactobacilli. Arch Oral Biol. 2015;60:1721–8.
Cheung SG, Goldenthal AR, Uhlemann AC, Mann JJ, Miller JM, Sublette ME. Systematic review of gut microbiota and major depression. Front Psychiatry. 2019;10:34.
Xiao S, Fei N, Pang X, Shen J, Wang L, Zhang B, et al. A gut microbiota-targeted dietary intervention for amelioration of chronic inflammation underlying metabolic syndrome. FEMS Microbiol Ecol. 2014;87:357–67.
Our sincere gratitude is extended to Xianghui Yun, Lina Xing, and Yun Li for their efforts in animal experiments and sample collection. We also acknowledge the support from the Majorbio Company (Shanghai, China).
This work was supported by Program of Shanghai Academic/Technology Research Leader (Grant Number: 17XD1403500), Innovation Project for Undergraduates of Shanghai University of Traditional Chinese Medicine (2019SHUTCM182), the National Natural Science Foundation of China (Grant Number: 81773927, 81573814) and Xinglin Young Scholar (2017.10-2020.9 to Jiaye Jiang).
Ethics approval and consent to participate
All experimental protocols were approved by the animal committee of Shanghai University of Traditional Chinese Medicine.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Neurobehavioral alternations in ACTH-treated rats. (A) Immobile time of forced swimming test (FST). (B) Immobile time of tail suspension test (TST). Data are expressed as mean ± SD (n = 10). **P < 0.01 compared with Control group. Fig. S2. Changes of serum biochemical parameters in ACTH-induced depression rats. The serum concentrations of 5-HT (A), NE (B), CRH (C), ACTH (D), and CORT (E). Data are expressed as mean ± SD (n = 10). ***P < 0.001, **P < 0.01 compared with Control group. Fig. S3. OPLS-DA score plots showing the distributions of urinary metabolites in ACTH group and Control group (n = 6). Fig. S4. Ratio of Firmicutes to Bacteroidetes in ACTH group and Control group. Data are expressed as mean ± SD (n = 6). **P < 0.05 compared with Control group. Fig. S5. (A) Genus-level taxonomic distributions of the microbial communities in fecal contents. (B) The relative abundance of bacterial genus detected in fecal samples. Data are expressed as mean ± SD (n = 6). *P < 0.05, **P < 0.01 compared with Control group.
Differential metabolites between ACTH group and Control group.
About this article
Cite this article
Song, J., Ma, W., Gu, X. et al. Metabolomic signatures and microbial community profiling of depressive rat model induced by adrenocorticotrophic hormone. J Transl Med 17, 224 (2019). https://doi.org/10.1186/s12967-019-1970-8