key: cord-0883965-yduf3niz
authors: Guo, Yinrui; Zhu, Xiangxiang; Zeng, Miao; Qi, Longkai; Tang, Xiaocui; Wang, Dongdong; Zhang, Mei; Xie, Yizhen; Li, Hongye; Yang, Xin; Chen, Diling
title: A diet high in sugar and fat influences neurotransmitter metabolism and then affects brain function by altering the gut microbiota
date: 2021-05-27
journal: Transl Psychiatry
DOI: 10.1038/s41398-021-01443-2
sha: 0f4f1ef0a9bf4009aedc78b4169925a8261a4d24
doc_id: 883965
cord_uid: yduf3niz

Gut microbiota (GM) metabolites can modulate the physiology of the host brain through the gut–brain axis. We wished to discover connections between the GM, neurotransmitters, and brain function using direct and indirect methods. A diet with increased amounts of sugar and fat (high-sugar and high-fat (HSHF) diet) was employed to disturb the host GM. Then, we monitored the effect on pathology, neurotransmitter metabolism, transcription, and brain circularRNAs (circRNAs) profiles in mice. Administration of a HSHF diet-induced dysbacteriosis, damaged the intestinal tract, changed the neurotransmitter metabolism in the intestine and brain, and then caused changes in brain function and circRNA profiles. The GM byproduct trimethylamine-n-oxide could degrade some circRNAs. The basal level of the GM decided the conversion rate of choline to trimethylamine-n-oxide. A change in the abundance of a single bacterial strain could influence neurotransmitter secretion. These findings suggest that a new link between metabolism, brain circRNAs, and GM. Our data could enlarge the “microbiome–transcriptome” linkage library and provide more information on the gut–brain axis. Hence, our findings could provide more information on the interplay between the gut and brain to aid the identification of potential therapeutic markers and mechanistic solutions to complex problems encountered in studies of pathology, toxicology, diet, and nutrition development.

Gut microbiota (GM) metabolites can potentially modulate nearly all aspects of host physiology 1 , from regulating immunity 2 and metabolism 3 in the gut to shaping mood and behavior 4 . These metabolites can act locally in the intestine or can accumulate up to millimolar concentrations in the serum and organs 5 . Studies have shown that formation of a gut-brain neural circuit for sensory transduction of nutrients enables the gut to inform the brain of all occurrences, and make sense of what has been eaten 6 . Recent studies have revealed that the GM is important in neurodegenerative diseases, including Alzheimer's disease (AD) 7, 8 and Parkinson's disease (PD), and that targeting the GM or/GM metabolites could be used to treat neurodegenerative diseases 9, 10 .

An imbalanced diet that includes a high intake of sugar and fat and insufficient dietary fiber over a long time can cause enteric dysbacteriosis. The latter increases the permeability of the intestinal mucosa, and results in abnormalities in intestinal immunity and glucolipid metabolism. At this time, the dominant bacteria can change readily. For example, Fujisaka and colleagues showed that the relative abundance of Bifidobacterium species and Bacillus species decreased in mice fed a highfat diet, whereas the abundance of Gram-negative bacteria increased 11 . Recent studies have shown that obesity is not necessary for dysfunction of the intestinal barrier. That is, hyperglycemia is more likely to drive intestinal-barrier dysfunction and the risk of enteric infection. Hyperglycemia increases the permeability of the intestinal barrier, which provides microbes with more chances to enter the body and cause proliferation of pathogenic bacteria and focal shifts 12 . Studies have demonstrated that some bacteria can produce bioactive neurotransmitters 13, 14 , and these neurotransmitters are thought to regulate the nervous system activity and behavior of the host 15, 16 . Recently, O Donnell and coworkers revealed that the neuromodulator tyramine produced by commensal bacteria of Providencia species (which colonize the gut) bypassed the requirement for host tyramine biosynthesis and manipulated a host sensory decision in Caenorhabditis elegans 17 . However, how these bacteria release signals to activate the brain is not known.

Using animal models, several pathways of communication have been identified along the "gut-brain-axis", including those driven by the immune system, vagus nerve, or by modulation of neuroactive compounds by the microbiota 18, 19 . In recent years, microbiota have been shown to produce and/or consume a wide range of mammalian neurotransmitters, including dopamine, norepinephrine, serotonin, or gamma-aminobutyric acid 13, 20 . Accumulating evidence in animals suggests that manipulation of these neurotransmitters by bacteria may have an impact in host physiology. Preliminary clinical studies have revealed that microbiota-based interventions can also alter neurotransmitter levels 13 . Nonetheless, substantially more work is required to determine if microbiota-mediated manipulation of human neurotransmission has physiological implications and, if so, how it may be exploited therapeutically.

We chose a diet with increased amounts of sugar and fat (i.e., a high-sugar and high-fat (HSHF) diet) to disturb the GM, then monitored the effect on pathology, neurotransmitters, metabolism, and transcription of circu-larRNAs (circRNAs) in mice. We aimed to identify some associations between the functions of the GM, neurotransmitters, and the brain. We also aimed to enlarge the "microbiome-transcriptome" linkage library. This approach would provide more information on the interplay between the gut and brain to aid identification of potential therapeutic markers and mechanistic solutions to complex problems encountered in studies on pathology, toxicology, diets and nutrition development.

Ethical approval of the study protocol All experimental protocols were approved by the Center of Laboratory Animals of the Guangdong Institute of Microbiology (Guangzhou, China). All efforts were made to minimize the number of animals used.

Adult male KM mice (18-22 g, 6 weeks) were obtained from the Center of Laboratory Animals of Guangdong Province (certificate number: SCXK [Yue] 2008-0020, SYXK [Yue] 2008-0085). They were pair-housed in plastic cages in a temperature-controlled (25 ± 2°C) colony room at a 12-h light-dark cycle. Food and water were available ad libitum.

Mice were allocated randomly into two groups of 12: control and model. The mice in the control group were fed a standard diet. The mice in the model groups were fed a HSHF diet. Water was available freely. These treatments lasted for 3 months.

The components of the HSHF diet were 20% sucrose, 15% fat, 1.2% cholesterol, 0.2% of bile acid sodium, 10% casein, 0.6% calcium hydrogen phosphate, 0.4% stone powder, 0.4% premix, and 52.2% basic feed. Heat ratio: protein 17%, fat 17%, carbohydrate 46%.

Male SAMP8 mice (5 months; mean bodyweight, 20 ± 5 g) were purchased from Beijing HFK Bioscience All mice were pair-housed in plastic cages in a temperature-controlled (25 ± 2°C) colony room with a 12-h light-dark cycle. Mice had free access to food and water. All animals were allowed to acclimatize to their surroundings for ≥1 week before initiation of experimentation.

Twenty male Sprague-Dawley rats (180-220 g) obtained from the Center of Laboratory Animal of Guangdong Province (SCXK [Yue] 2008-0020, SYXK [Yue] 2008-0085). They were pair-housed in plastic cages in a temperature-controlled (25 ± 2°C) colony room at a 12-h light-dark cycle. Food and water were available ad libitum.

Male Sprague-Dawley rats were divided into two groups of 10: normal group, trimethylamine (TMA)induced group (2 mL/kg/d of 2.5% TMA purchased from Shanghai Aladdin Biochemical Technology, Shanghai China).

The abundance of C. albicans (a symbiotic opportunistic pathogen in humans) and K. pneumoniae was found to be high in AD patients in our other study (data not shown). We chose these pathogens to ascertain the potential routes of communication/interaction between the host and its resident bacteria on neurotransmitter metabolism and brain function. C. albicans and K. pneumoniae were administered (i.g.) to normal C57 mice by monotherapy or in combination.

Adult male C57 mice (18-22 g, 6 weeks) obtained from Center of Laboratory Animal of Guangdong Province, SCXK [Yue] 2008-0020, SYXK [Yue] 2008-0085 were pairhoused in plastic cages in a temperature-controlled (25 ± 2°C) colony room at a 12-h light-dark cycle. Twenty C57 mice were divided into four groups: normal, C. albicanstreated, K. pneumoniae-treated, and C. albicans + K. pneumoniae-treated (Ca + Kp).

The appearance, behavior, and fur color of animals were documented every day. Bodyweight was measured every 3 days during the period of drug administration. Blood samples were drawn by removing the eye under anesthesia (isoflurane). Serum was acquired by centrifugation and stored at −80°C until measurement. Levels of triglycerides (TG), total cholesterol, (T-CHO) and highdensity lipoprotein-cholesterol (HDL-C) were measured with commercially available kits from Nanjing Jiancheng Bioengineering Institute (Jiangsu, China). Serum and brain-tissue levels of trimethylamine-n-oxide (TMAO), and neurotransmitter levels were quantified using liquid chromatography-mass spectrometry (LC-MS).

Brain, liver, renal, spinal marrow, spleen, and adipose tissues were removed and fixed in 4% paraformaldehyde at pH 7.4 for further pathologic observation. These tissue samples were made into paraffin sections after drawing materials, fixation, washing, dehydration, transparency, dipping in wax, and embedding. Obesityrelated parameters or other related pathologic changes were measured.

The brains of animals were dissected. Four brains from each group were fixed in 4% paraformaldehyde solution and prepared as paraffin sections. Sections were stained with hematoxylin and eosin (H&E). silver, or underwent Nissl staining and TUNEL staining. Immunostaining using paraffin-embedded sections (thickness = 3 μm) and a two-step method involving a peroxidaseconjugated polymer kit (Envision ® ; DAKO, Carpinteria, CA, USA) was also done. Slides were observed under light microscopy.

Fresh samples of intestinal content were collected 12 h before the fasting of rats and stored at −80°C. Microbial DNA isolated from these samples, with a total mass ranging from 1.2 ng to 20.0 ng, was stored at −20°C. Microbial 16S rRNA genes were amplified using a forward primer (5′-ACTCCTACGGGAGGCAGCA-3′) and a reverse primer (5′-GGACTACHVGGGTWTCTAAT-3′). Polymerase chain reaction (PCR) amplification was done and the products were purified, quantified, and homogenized to form a sequencing library. The established library was firstly inspected by the library, and the qualified library was sequenced by Illumina HiSeq 2500.

Aliquots of each standard solution were mixed to generate a stock standard mixture of 4 μg/mL in 50% acetonitrile. This 4 μg/mL standard mixture (100 μL) was mixed with 50 μL of sodium carbonate (100 mM) and 50 μL of 2% benzoyl chloride (BZ) or 2% benzoyl chloride-(phenyl-13 C 6 ) ( 13 C 6 BZ). The reaction mixture was vortexmixed and diluted to 2400 ng/mL for BZ derivatives and 50 ng/mL with 50% acetonitrile for 13 C 6 BZ derivatives as the internal standard (IS) stock solution. BZ derivatives were serially diluted to 240, 120, 60, 30, 12, 3, 1.2, and 0.02 ng/mL. To prepare the standard curve, the standard BZ derivative solutions stated above were mixed isometrically with 50 ng/mL of IS stock solution to generate calibration levels covering a range of 0.01-1200 ng/mL for all analytes.

Samples of brain tissue (20 ± 1 mg) were homogenized in a 4-volume (vol/wt) precooled aqueous solution of ascorbic acid (20 mM) with TissueLyser™ JX-24 (Jingxin, Shanghai, China) beads at 30 Hz for 3 min. The homogenates were supplemented with 15-volume prechilled acetonitrile (−40°C) and vortex-mixed for 1 min before centrifugation at 4°C and 14,000 × g for 15 min. Derivatization was started by addition of 40 μL of 2% BZ to a mixture of 40 μL of supernatant solution and 20 μL of sodium carbonate (100 mM). Derivatized samples were mixed with IS stock solution to analyze high concentrations of neurotransmitters by ultra-high pressure-tandem mass spectrometry (UPLC-MS/MS) with an injection volume of 1 μL. The derivatized sample was diluted 10-fold, and the dilution was mixed isometricallywith ISs to analyze low concentrations of neurotransmitters by UPLC-MS/MS with an injection volume of 5 μL.

Frozen samples of intestinal content (20 ± 1 mg) were homogenized in 5-volume (vol/wt) of a precooled aqueous solution of ascorbic acid (20 mM). After sonication for 3 min, the homogenates were supplemented with 15-volume prechilled acetonitrile (−40°C) and vortexmixed for 1 min before centrifugation at 4°C and 14,000 × g for 15 min. Derivatization was started by addition of 40 μL of 2% BZ to a mixture of 40 μL of supernatant and 20 μL of sodium carbonate (100 mM). Derivatized samples were mixed with IS stock solution to analyze low concentrations of neurotransmitters by UPLC-MS/MS with an injection volume of 1 μL. The derivatized sample was diluted 10-fold, and the dilution was mixed isometrically with ISs to analyze high concentrations of neurotransmitters with UPLC-MS/MS with an injection volume of 10 μL.

UPLC-MS/MS was done on a Acquity UPLC system (Waters, Milford, MA, USA) coupled to a Triple Quad™ 5500 tandem mass spectrometer (AB Sciex, Framingham, MA, USA). Each sample or standard mixture was injected onto a UPLC BEH C18 column (100 mm × 2.1 mm, 1.7 μm; Waters) at a flow rate of 0.4 mL/min. The mobile phase consisted of 2 mM ammonium acetate with 0.1% formic acid in water (A) and acetonitrile (B). Chromatographic separation was conducted by a gradient elution program: 0 min, 1% B; 0.5 min, 1% B; 1 min, 45% B; 4 min, 65% B; 4.2 min, 100% B; 5.2 min, 100% B; 5.3 min, 1% B; 7 min, 1% B. The column temperature was 40°C.

The analytes eluted from the column were ionized in an electrospray ionization source in positive mode. The conditions were: source temperature, 600°C; curtain gas, 30 psi; ion source gas 1, 50 psi; ion source gas 2, 50 psi; collision gas, 8 psi; ion spray voltage, 5500 V; entrance potential, 10 V; collision cell exit potential, 14 V. Scheduled multiple reaction monitoring (sMRM) was used to acquire data in optimized MRM transition (precursor > product), declustering potential, and collision energy ( Table 1 ). The total scan time per cycle was 0.25 s. The samples and standard-curve samples were analyzed simultaneously. AB Sciex Analyst v1.5.2 (https://sciex. com/products/software/analyst-software) was used to control instruments and acquire data using the default parameters and assist manual inspection to ensure the qualitative and quantitative accuracies of each compound. The peak areas of target compounds were integrated and outputted for quantitative calculation.

Total RNA was isolated using QIAzol™ (Qiagen, Germany) and miRNeasy™ kits (Qiagen, Germany), including additional DNase I digestion. Then, ribosomal RNA was removed using the Ribo-Zero™ Magnetic Gold kit (Qiagen, Germany), and enzymatic degradation of linear RNA was done using Rnase R enzyme. Fragmentation buffer was added to fragments. The first chain of complementary (c)DNA was synthesized with six-base random hexamers. Then, the buffer, dNTPs, RNase H, and DNA polymerase I were added to synthesize the second chain of cDNA. Purification using the QiaQuick™ PCR kit and elution using EB buffer after terminal repair, processing of base A, and addition of sequencing joints were carried out. For next-generation sequencing, 0.5 μg of ribosomal RNAdepleted RNA was fragmented and primed. Sequencing libraries were constructed using TruSeq™ RNA Sample Preparation kits (Illumina, San Diego, CA, USA) and were sequenced by Illumina HiSeq™ 2500 flow cells.

First, reads were mapped to the latest University of California Santa Cruz transcript set using Bowtie2 v2.1.0 (http://bowtie-bio.sourceforge.net/bowtie2/index.shtml). Gene expression was estimated using RSEM v1.2.15. For analyses of lincRNA expression, we used the transcripts set from Lncipedia (www.lncipedia.org/). Trimmed mean of M-values (TMM) was used to normalize gene expression. Differentially expressed genes were identified using the edgeR program. Genes showing altered expression with p < 0.05 and >1.5-fold changes were considered to have differential expression. Pathway analysis and network analysis were undertaken using Ingenuity (IPA). IPA computes a score for each network according to the fit of a set of supplied focus genes. These scores indicate the likelihood of focus genes belonging to a network versus those obtained by chance. Score >2 indicates a ≤99% confidence that a focus gene network was not generated by chance alone. The canonical pathways generated by IPA are the most significant for the uploaded dataset. Fischer's exact test with a false discovery rate (FDR) option was used to calculate the significance of the canonical pathway.

For analyses of circRNA expression, reads were to mapped a genome using STAR. DCC was employed to identify and estimate expression of a particular circRNA. TMM was used to normalize gene expression. Differentially expressed genes were identified using edgeR. miRanda was employed to predict the miRNA target of the circRNA. R (R Institute for Statistical Computing, Vienna, Austria) was used to generate figures.

circRNA verification by real-time reverse transcriptionquantitative polymerase chain reaction (RT-qPCR)

We wished to validate the reliability of high-throughput RNA-sequencing (RNA-seq) and to explore expression of circRNAs during aging. Hence, expression of circRNAs was measured by RT-qPCR. With reference to the method described by Memczak, two sets of primers for each cir-cRNA were designed using Primer Express v5.0 ( Table 2) : an outward-facing set to amplify only the circRNA, and an opposite-directed set to amplify the linear form. Total RNA was extracted (TRIzol® Reagent, Life Technologies, Carlsbad, CA, USA), digested using RNase R, and purified. cDNA was synthesized using the Geneseed ® II First Strand cDNA Synthesis kit (Agilent Technologies, Santa Clara, CA, USA). Outward-facing primers were designed to amplify the fragment across the junction from cDNA, then the fragment was sequenced by Sangon Biological Engineering (Shanghai China). RT-qPCR was undertaken using Geneseed ® qPCR SYBR ® Green Master Mix (Agilent Technologies). PCR-specific amplification was conducted with an ABI 7500 system (Applied Biosystems, Foster City, CA, USA). Expression of cir-cRNAs was defined based on the threshold cycle (Ct), and relative expression was calculated via the 2 −ΔΔCt method. Glyceraldehyde 3-phosphate dehydrogenase served as an IS control with all reactions done in triplicate.

Briefly, global brain tissue was dissected from treated mice (purchased from Beijing HFK Bioscience; SCXK (Jing) 2014-0004) and proteins extracted with a radioimmunoprecipitation assay lysis buffer (T-PER ™ Tissue 

Data are the mean ± SD of at least three independent experiments. Significant differences between treatments were analyzed by one-way analysis of variance (ANOVA) at p < 0.05 using SPSS (IBM, Armonk, NY, USA) and Prism 5 (GraphPad, San Diego, CA, USA).

The HSHF diet disrupts the GM The bodyweights of mice fed a HSHF diet were higher than those fed a standard diet (control, p < 0.05) (Fig. 1A) ; levels of blood glucose (Fig. 1B) and TMAO (Fig. 1C) were also higher (p < 0.05). These data suggested that the HSHF diet utilized caused increases in bodyweight, blood glucose level, and TMAO level. Sequencing of the 16 S rRNA gene showed that the HSHF diet significantly decreased the operational taxonomic units (OTUs) of bacteria of the phyla Acidobacteria, Verrucomicrobia, Tenericutes, and Firmicutes, while increasing the OTUs of bacteria of the phyla Bacteroidetes, Proteobacteria, Deferribacteres, Cyanobacteria and Actinobacteria (p < 0.05) ( Fig. 1D and (Fig. 1E and Fig. S1B) following administration of the HSHF diet. The whole GM was different between the two groups of mice: control mice and those fed a HSHF diet. Principal component analysis (Fig. S1C) , Venn diagram ( Fig. S1D) , and linear discriminant analysis effect size (LEfSe) analysis ( Fig. S1E and S1F) could be used to distinguish between these two groups. Hence, the HSHF diet could induce dysbacteriosis in mice. These results are almost identical to data from studies showing that intestinal dysbacteriosis interacts with fat 21 .

The HSHF diet changes the steady state of the intestine Intestinal dysbacteriosis were considered first affects the intestinal physiology 22 . The pathology of the small intestine and colon of mice fed the HSHF diet was different from that of mice fed the standard diet. Changes induced by the HSHF diet included: cell shrinkage; reduction of cell turnover; cytoplasmic vacuolar changes; blurred and indistinct cell boundaries; severe shedding of intestinal villi; scattering of damaged tissue blocks (Fig. 1F, G) .

Immunohistochemical analyses showed that expression of some proinflammatory markers was affected. Expression of tumor necrosis factor (TNF)-α and nuclear factorkappa B (NF-κB p65) was activated (p < 0.05) (Fig. 1H-K) , whereas expression of interleukin (IL)−2 and peroxisome proliferator-activated receptor (PPAR)-γ was inhibited (p < 0.05) (Fig. 1H-K) . RNA-seq of small-intestinal tissues revealed that expression of 42 RNAs was upregulated and that of 68 RNAs was downregulated (|log2FC| > 1.0, FDR < 0.05 vs. control) (Table S1, Fig. S2A ). For example, expression of Hspa1a, Hspa1b, P2ry4, Enpp7, Ano3, Hsph1, Nos1ap, Slc5a12, Slc5a4a, Gm5286, Paqr9, G0s2, Trim50, Lep, Apoa2, Retnlb, Gsdmc, TCONS_00019039, Fam205a4, Gm20708, Exosc6, Pde2a and Gm10184 mRNAs was changed, and these were all related to intestinal injury/leakage, inflammation, and immunity. Pathway analysis using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database showed that the important enrichment pathways were the: "PPAR signaling pathway" (mRNA of APOA2, Acaa1b, Cyp4a10, Fabp1, Hmgcs2, Me1); "biosynthesis of unsaturated fatty acids"; "fructose and mannose metabolism"; "glycolysis/gluconeogenesis"; "fatty acid elongation"; "carbon metabolism"; "renin secretion"; "renin-angiotensin system"; "AMPK signaling pathway" (mRNA of G6pc, Fbp1, Srebf1, Lep). Together these results demonstrated that dysbacteriosis influenced intestinal injury/leakage, the inflammatory response, and energy metabolism. Furthermore, histopathology of the liver (Fig. S2B) , kidney (Fig. S2E) , spinal marrow (Fig. S2E) , spleen (Fig. S2D) , adipose (Fig. S2C) and heart (Fig. S2F) tissues showed that the HSHF diet not only induced dysbacteriosis, but also activated inflammation and damage to multiple organs.

The HSHF diet influences the brain cholinergic system and inflammation

The HSHF diet was administered to mice to disturb the GM. Brain sections from mice were analyzed using E Influences of a diet high in sugar and fat on the gut microbiota at the genus level, see Fig. S2 . F, G Influences of a diet high in sugar and fat on intestine (F) and colon (G) histopathology using H&E staining, see the pathologic observations in other organs (livers, renal, spinal marrow, spleen, adipose and heart) in Fig. S2 . H-K Influences of a diet high in sugar and fat on intestine and colon histopathology using immunohistochemical analyses. Data are the mean ± SD of more than 8 independent experiments (the histopathology data were the mean ± SD of more than 3 independent experiments). *p < 0.05 and **p < 0.01 vs. the model group by one-way ANOVA, followed by the Holm-Sidak test. various stains. H&E (Fig. 2Aa, Fig. S3A-A) , Nissl (Fig. 2Ab, D and Fig. S3A-B) , silver (Fig. 2Ac , E and Fig. S3A -C) and TUNEL (Fig. 2B-F and Fig. S3A-D) staining showed obvious changes in pathology in mice with HSHF diet-induced dysbacteriosis: reduction in the size and turnover of neurons; cytoplasmic vacuolar changes; nerve-fiber reduction. In particular, the apoptotic percentage in the hypothalamus was increased in mice suffering from dysbacteriosis, which suggested that the appetite of HSHF-diet mice increased and bodyweight soared, as seen previously 23 . Immunofluorescent staining of microglia and astrocytes using antibodies to IBA-1 (Fig. 2C, G, Fig. S3B -B) and GFAP (Fig. 3A, Fig. S3B-A) showed that the number of astrocytes was reduced significantly and microglial cells were activated in mice suffering from dysbacteriosis. Expression of the cholinergic neuron of AChE (Acetylcholinesterase) (Fig. 3B , p < 0.05), Amphiphysin (AMP) (Fig. 3B , p < 0.05), Acetylcholine receptor of CHRNB1 and CHRNA1 measurement showed that the cholinergic system was affected.

Expression of TNF-α (p < 0.05) (Fig. 3C) and NF-κB p65 (p < 0.05) (Fig. 3C) was activated, whereas expression of PPAR-γ and IL-2 showed no obvious differences (p > 0.05) (Fig. 3C) . These data demonstrated that atrophy, inflammation, or the immune response were imbalanced in the brains of mice with HSHF diet-induced dysbacteriosis.

A total of 7362 circRNAs were identified in the brains of mice fed a HSHF diet (Fig. 3) . Among them, 287 differentially expressed circRNAs were in the control group vs. model group (fold change >1.50, p < 0.05) (Fig. 3D , and Table 1 , circBase): expression of 88 cir-cRNAs was upregulated, and expression of 199 circRNAs was downregulated. This finding indicated that the cir-cRNAs had different expression profiles in slender and obese mice. Bioinformatics analysis of differentially expressed circRNAs showed that dysbacteriosis could influence the brain in terms of metabolism, synaptic Table 1 . Data are the mean ± SD of 3 independent experiments. *p < 0.05 and **p < 0.01 vs. the model group by one-way ANOVA, followed by the Holm-Sidak test. transmission and plasticity, and endogenous hormone levels. These actions would affect the activities and functions of the brain, but also affect heart function (e.g., circadian rhythm, dilated cardiomyopathy, and arrhythmogenic right ventricular cardiomyopathy). These data demonstrated that dysbacteriosis could affect brain circRNA-seq directly or indirectly.

To identify the influences of the GM on neurotransmitters in the intestine, combined data analyses of the GM and neurotransmitters in fimo were done (R 2 > 0.85). Some members of Aerococcaceae, Corynebacteriaceae, Brucellaceae, F16, Paraprevotellaceae, Veillonellaceae, and Dehalobacteriaceae were the key bacteria related to the production and/or secretion of neurotransmitters (Fig. 4A) . They influenced proline, glycine, leucine, and serine (Fig.  4A) . Phenylalanine, tryptophan, and tyrosine were influenced by bacterial in Veillonellaceae, Dehalobacteriaceae, Paraprevotellaceae, Coriobacteriaceae, Brucellaceae, and F16, but how these bacteria produce them, and then participate or intervene in host metabolic pathways, merits further study.

The norepinephrine level was influenced by bacteria from Streptophyta, RF23, Flavobacteriaceae, Lactobacillaceae, Ruminococcaceae, Turicibacteraceae, Clostridiales, Mogibacteriaceae and Peptostreptococcaceae (Fig. 4A) . The 5-hydroxy tryptamine (5-HT) level was influenced by bacteria from Streptococcaceae, Helicobacteraceae, Clostridiales, Turicibacteraceae, Clostridiaceae, and Brucellaceae (Fig. 4A) . The level of 5-Hydroxyindoleacetic acid (5-HIAA), the primary metabolite of 5-HT, was influenced by bacteria from YS2, Lactobacillaceae, Turicibacteraceae, Christensenellaceae, Peptostreptococcaceae, Paraprevotellaceae, and Bacteroidales (Fig. 4A ).

Upon combination with data for intestinal mRNAs, we found that intestinal neurotransmitters also influenced intestinal genes, and then bacteria. As shown in Fig. 4B , leucine was positive to the mRNA of glutathione S-transferase Mu 3 (Gstm3), integrin subunit alpha L (Itgal), Krt12, small proline rich protein 2 A 2 (Sprr2a2), small proline rich protein 2 A 1 (Sprr2a1), tropomyosin 3 (Tpm3), dual specificity phosphatase 6 (Dusp6), and then to bacteria from Coprobacillus, Enterococcus, Enterobacter, Enterobacteriaceae, and Peptostreptococcaceae. The gene of cytochrome P450 family 2 subfamily E member 1 (Cyp2e1) influenced taurine, tyrosine, and 4-Aminobutyric acid (Fig. 4C) , and then became involved with the bacteria of Akkermansia, Allobaculum, Coriobacteriaceae, Alcaligenaceae, Corynebacterium Alistipes, and Facklamia (Fig. 4C) . Tenascin R was involved with several types of bacteria (Fig. 4C) .

The bacteria of Flavobacteriaceae, Porphyromonadaceae, Bifidobacteriaceae, Bacteroidaceae, Desulfovibrionaceae Enterococcaceae, Prevotellaceae, Lactobacillaceae and Peptostreptococcaceae could interact with most brain neurotransmitters: isoleucine, valine, histamine, serine, proline, leucine, phenylalanine, tryptophan, tyrosine, glutamic acid, 3,4-Dihydroxyphenylalanine, kynurenine, 5-Hydroxyindoleacetic acid, tryptamine, glycine, and ornithine (Fig. 5A) . The brain neurotransmitters homovanillic acid and 3,4-Dihydroxyphenylacetic acid were influenced by bacteria from Bacteroidales, YS2, Paraprevotellaceae, Enterobacteriaceae and Bacteroidaceae (Fig. 5A) .

The "cholinergic hypothesis" has an important role in AD therapy 24 . The degenerative dysfunction of cholinergic neurons is a partial cause of memory deficit in dementia patients. The neurotransmitter acetylcholine (ACh) in the brain was implicated with several intestinal genes (Fig. 4D , Table S1 and enrichment analysis using the KEGG database), and might be intervened and/or promoted by bacteria from Paraprevotella, [Ruminococcus], Streptophyta and Comamonadaceae (Fig. 4D) . These data indicated that the GM has complex interactions with the host cholinergic system. We discovered that the norepinephrine level in the brain could: (i) influence bacteria of Peptostreptococcaceae, Sutterella, Flavobacteriaceae and Coprobacillus; interact with bacteria of Enterococcus, Enterobacteriaceae, and Enterobacter, and then interact with 3,4-dihydroxyphenylacetic acid; interact with bacteria from Lactobacillus, and then influence 5hydroxyindoleacetic acid (Fig. 5B) .

Glutamine, isoleucine, phenylalanine, tyrosine, leucine, tryptophan, proline and valine did not interact with a single type of bacteria (Fig. 5B) . Hence, the interactions between GM, neurotransmitters, and host gene are complicated. More other complex interactions between the GM, mRNA-seq of intestinal tissue, and neurotransmitters of brain tissues are shown in Figs. 4 and 5. The neurotransmitters ACh, 4-aminobutyric acid, 5-hydroxyindoleacetic acid (5-HIAA), glutamine, serotonin, leucine, and kynurenine may also have important roles in the brain-GM axis, or the GM may be indispensable for the synthesis and secretion of neurotransmitters: the exact mechanism needs much more work.

Network interaction of the GM, brain neurotransmitters, and circRNAs Next, the interactions of the GM, neurotransmitters, and circRNAs in brain tissues were analyzed. The ACh level in brain tissues was important for the interactions between the brain circRNA of chr5_21278044/21281602_+, chrX_162934286/162941844_+ and chr16_37121798/ 37129964_−, and bacteria from Paraprevotella, Helicobacteraceae, Erysipelotrichaceae, [Ruminococcus] and Streptophyta (Fig. 5C) . The neurotransmitters 5-HIAA, leucine, 3,4-dihydroxyphenylacetic acid, kynurenine, ornithine, tryptamine, ACh, and tyrosine were the main metabolites influencing the GM and circRNAs (Fig. 5D ), or they could influence multiple microorganisms and circRNAs (Fig. 5) . Much more experimental evidence is needed to explain these associations.

5-HIAA is the main metabolite of 5-HT. A recent study showed that 5-HT directly stimulates and inhibits the growth of commensal bacteria in vitro, and exhibits concentration-dependent and species-specific effects 25 . Some bacteria may have important roles in the conversion of 5-HT to 5-HIAA. Bacteria from the family Peptostreptococcaceae could affect the mRNA of Tpm3 and Dusp6 (Figs. 4 and 5) . Research has revealed a genetic proof-of-function for choline acetyltransferase (ChAT) in T cells during viral infection and identified a pathway of T-cell migration that sustains antiviral immunity 26 . Levels of ACh and ChAT were influenced primarily by bacteria from Helicobacteraceae, Paraprevotella, [Ruminococcus], Erysipelotrichaceae, Comamonadaceae, and Streptophyta (Figs. 4 and 5). The role of the GM in TMAO byproducts from choline has been demonstrated 27 , and TMAO has been shown to be harmful in several diseases 28 . Therefore, the GM acts as a bridge to contribute to host nutrition/health by regulating the metabolism. The above-mentioned data, together with data from other studies, suggest that gut dysbacteriosis could impact circRNA sequences in the brain.

In this additional experiment, expression of some cir-cRNAs in a murine microglial cell line (BV2) was monitored after incubation with TMAO (18.45 mg/mL) (p < 0.05) (Fig. 6A) . Expression of MAOA (Amine oxidase [flavincontaining] A), MAOB, COMT (Catechol O-methyltransferase), AChE, AMP, CHRNA1 and CHRNB1 was influenced by this concentration of TMAO (p < 0.05) (Fig. 6B) , but a low concentration of TMAO did not cause obvious damage to BV2 cells. Percentage inhibition in BV2 cells did not differ at a TMAO concentration <16 mg/mL (Fig. 6A) . Expression of some circRNAs, including circNF1-419 and circ_0001239, was influenced by a low concentration of TMAO (p < 0.05) (Fig. 6C) , indicating that some circRNAs might be sensitive responsive signaling molecules to TMAO. These data suggested that microbial metabolites might influence the formation and degradation of host circRNAs.

Levels of MAOA, MAOB, and COMT were affected after TMAO administration to BV2 cells (Fig. 6) . Levels of some monoamine neurotransmitter were changed in mice fed the HSHF diet (Figs. 4 and 5) . Taken together, these data suggested that behavior and stress would be influenced and, in these processes, the GM had an accelerant role. Therefore, the GM decides the transformation efficiency from choline to TMAO. The latter influences the levels of MAOA, MAOB, and COMT, which participate in the regulation of metabolism of neuroactive and vasoactive amines in the central nervous system (CNS) and peripheral tissues. circRNAs spliced from Rims1/Syt1/Unc13b could be enriched in the ACh release cycle (Fig. 7A) . Levels of some of these circRNAs were measured in the brain tissues from different mice after choline treatment. Choline and/or its metabolite could change circRNA expression (Fig. 7B) , which suggested that the GM also influences the absorption and conversion of choline. Pecoraro and colleagues suggested that a high-affinity PAC1 (Proteasome assembly chaperone 1) receptor presynaptically modulates hippocampal glutamatergic transmission acting through AChE in hippocampal CA1 29 , indicating that choline plays an important part in neurotransmitter homeostasis.

The GM plays an important part in lipid homeostasis in blood, and in the causes and development of nervoussystem diseases. In this study, the levels of total cholesterol (T-CHO, Fig. 7Ca) , triglyceride (TG, Fig. 7Cb ) and lowdensity lipoprotein (LDL, Fig. 7Cc) were also increased Fig. 5 Dysbacteriosis affects neurotransmitter metabolism by targeting the gut-brain axis in mice fed a diet high in sugar and fat for 4 months. A Combined data analyses between the gut microbiota (genera) and neurotransmitters in brain tissues. B Combined data analyses between the gut microbiota (genera), transcriptome in the intestine, and neurotransmitters in brain tissues. C, D Combined data analyses between the gut microbiota (genera), neurotransmitters, and circRNAs in brain tissues. Data are the mean ± SD of more than eight independent experiments.

(p < 0.05), and the levels of serum TMAO (Fig. 3C) , when compared to the HSHF-diet group (p < 0.05), and the highdensity lipoprotein (HDL, Fig. 7Cd) were decreased (p < 0.05).

Rats treated with TMA could put on weight continuously (Fig. 8Aa) , had increased levels of TMAO in serum (Fig. 8Ab) and increased levels of total T-CHO (Fig. 8Ac) , TG (Fig. 8Ac) and LDL (Fig. 8Ac) (p < 0.05) when compared with those in the model group (p < 0.05), whereas the HDL content was reduced (p < 0.05) (Fig. 8Ac) . Data for sequencing of the 16 S rRNA gene showed that TMA decreased the OTUs of bacteria from Verrucomicrobia, TM7 and Acidobacteria (p < 0.05) (Fig. 8B) , while increasing the OTUs for bacteria from Tenericutes, Deferribacteres, Actinobacteria, and cyanobacteria (p < 0.05) (Fig. 8B) . At the genus level, the OTUs of bacteria from Bifidobacterium, Aerococcus, Facklamia, Blautia and Enterobacteriace were increased by TMA Fig. 6 Effects of TMAO on BV2 cells. A Effects of TMAO on proliferation of BV2 cells, measured using the MTT assay. B Expression of MAOA, MAOB, COMT, AChE, AMP, CHRNA1, and CHRNB1 in TMAO-damaged BV2 cells. C Expression of circNF1-419 and circ_0001239 in TMAO-damaged BV2 cells. D Expression of circNF1-419 and circ_0001239 in Kombucha tea-treated mice. E Ratio of abundance of Bifidobacterium species and Lactobacillus species was imbalanced in mice fed a HSHF diet and in AD-like mice. Data are the mean ± SD of more than five independent experiments (and three independent experiments when using western blotting). *p < 0.05 and **p < 0.01 vs. the model group by one-way ANOVA, followed by the Holm-Sidak test. , and protein expression in liver tissues after intragastric administration of TMA at 0.2 mL/100 g/d of 2.5% for 6 weeks. B Effects on the gut microbiota at the phylum level, and (C) at genus level. D Pathologic changes in hippocampus, thalamus, brainstem, and midbrain after TMA administration. Data are the mean or mean ± SD of more than 8 independent experiments. *p < 0.05 and **p < 0.01 vs. the model group by one-way ANOVA, followed by the Holm-Sidak test. administration (p < 0.05) (Fig. 8C) , whereas those of bacteria from Prevotella, Lactobacillus, Helicobacteraceae and Akkermansia were decreased significantly (p < 0.05) (Fig. 8C) , which indicated that excessive TMA in the diet can induce dysbacteriosis. Histopathologic changes in the hippocampus, thalamus, brainstem, and midbrain of the TMA-induced group compared with those in the control group (Fig. 8D) included cell shrinkage, cytoplasmic vacuolar changes, blurred and indistinct cell boundaries, and scattered damaged tissue blocks. Also, the differential expression of proteins in the liver from TMA-treated rats revealed that TMA induced damage in the host. For example, the level of mitochondrial brown fat uncoupling protein 1 was decreased after treatment with excess TMA, together with a decrease in Ndufb2 expression (Fig. 8Ad) , which indicated that energy metabolism was imbalanced due to dysfunction of NADH dehydrogenase activity and oxidoreductase. Meanwhile, Ncam1 expression was upregulated (Fig. 8Ad) , suggesting that neurite outgrowth, synaptic plasticity, as well as learning and memory might have been influenced. Increased expression of heat shock 70 kDa protein 1B and B-cell receptor-associated protein 31 (Fig. 8Ad) indicated that the host stress response was activated and immune-system function might have been compromised. In agreement with results from other studies, TMA resulted in negative consequences for the host 29, 30 and, in the present study, TMA led to negative consequences for the nervous system.

We wished to ascertain the potential routes of communication/interaction between the host and its resident bacteria with regard to neurotransmitter metabolism and brain function. C. albicans and K. pneumoniae were administered (i.g.) to normal C57 mice as monotherapy or in combination. Expression of AChE, AMP, CHRNA1, CHRNB1, and GAD65 was changed in the intestine (p < 0.05) (Fig. 9A ) and brain (Fig. 9B) , especially in the C. albicans-treated group in the intestine. These data indicated that changes in two core microorganisms influenced GM-mediated compounds and neurotransmitters.

Increasing numbers of studies have indicated that the GM is closely related to disorders of the nervous system, including autism 31 , depression 32 , schizophrenia 33 , anorexia nervosa 34 , multiple sclerosis 35 , epilepsy 36 , PD 37 , and AD 18 . In recent years, evidence has shown that circRNAs may be important for brain function: CDR1 is found at a high level in mammalian neurons 38 , and the circRNA transcripts from neuron genes can accumulate during aging in Drosophila 39 . Hentze and colleagues showed that circRNAs are dynamic during the activities of daily life 40 .

Also, circRNAs are expected to be biomarkers or drug targets for neurodegenerative diseases, and studies have shown that the CDR1 level is reduced in AD 38 . Here, we demonstrated that gut dysbacteriosis can implicate the brain circRNA-sequencing directly or indirectly.

Previously, using in an overexpressed circNF1-419 adeno-associated virus system, we showed that cir-cRNAs in the brain influenced the cholinergic system of the brain, and changed the GM composition, intestinal homeostasis/physiology, and even the GM trajectory in newborn mice 41 . Another system, circZCCHC11 (mmu_circ_0001239), showed a "sponge" miRNA function and mediated a series of chain reactions in the brain, which then influenced gut function and GM engraftment from their parent 41 . Those data demonstrated a link between circRNAs and the GM, enlarged the microbiome-transcriptome linkage library, and provided additional information on the gut-brain axis. Those data and the observations mentioned above indicated that the GM influences the conversion rate of choline to TMAO, and that the TMAO level influences the cholinergic system. In an additional experiment, expression of some circRNAs in the BV2 cell line was measured after incubation with TMAO, which indicated that some circRNAs might be sensitive responsive signaling molecules to TMAO. We concluded that microbial metabolites influenced the formation and degradation of host circRNAs.

Studies have demonstrated that catalyzes the oxidative deamination of biogenic and xenobiotic amines, and has important functions in the metabolism of neuroactive and vasoactive amines in the CNS and peripheral tissues. For example, MAOA preferentially oxidizes biogenic amines such as 5-HT, norepinephrine, and epinephrine. MAOB preferentially degrades benzylamine and phenylethylamine. COMT catalyzes the transfer of a methyl group from S-adenosylmethionine to catechol substrates such as the neurotransmitters dopamine, epinephrine, and norepinephrine, and degrades catecholamine transmitters 42 . All of the above-mentioned monoamine neurotransmitters have indispensable roles in the regulation of CNS functions.

We found that the levels of MAOA, MAOB, and COMT in BV2 cells were affected after TMAO administration, and that the levels of some monoamine neurotransmitters were changed in mice fed the HSHF diet. These data suggested that behavior and stress would be influenced and that, in these processes, the GM had an accelerant role. Therefore, we concluded that the GM decided the transformation efficiency from choline to TMAO. The latter then influenced the levels of MAOA, MAOB, and COMT, which participate in the regulation of metabolism of neuroactive and vasoactive amines in the CNS and peripheral tissues. How circRNAs are involved in this process merits further investigation. A major goal of any microbiome study is to move beyond correlation, and parse out potential routes of communication/interaction between the host and its resident bacteria. To identify their communication/ interaction, the species level was identified using metagenomics sequencing. Twelve species of Archaea were detected, including Methanobrevibacter sp. AbM4 (781.63-fold changes compared with control, Fig. S4A ), Methanosarcina sp. MTP4 (1.21-fold, Fig. S4A) , and especially Methanosalsum zhilinae, Methanomethylovorans hollandica, Candidatus Methanomethylophilus alvus, and Thermococcus cleftensis, which were new residents (Fig. S4A ). Most are methanogens (anaerobic prokaryotes from the domain Archaea that utilize hydrogen to reduce carbon dioxide, acetate, and various methyl compounds to methane 43 Methanogens are emerging pathogens associated with abscesses in the brain and muscles 44 . They have been implicated in dysbiosis of the oral microbiota, periodontitis, and peri-implantitis 43 . They have also been associated with dysbiosis of the digestive-tract microbiota linked to metabolic disorders (anorexia, malnutrition, and obesity) and with lesions of the digestive tract (colon cancer) 43 . One clinical investigation showed that a negative association between methane concentrations in breath and anthropometric biomarkers of obesity 45 , and that methane significantly decreased the neurological deficit induced by cerebral ischemia and reperfusion via the antioxidant pathway of PI3K/Akt/HO-1 46 . Special diets have been used to change anaerobic prokaryotes to involve the digestive and nervous systems 47 . ACh is a neurotransmitter in mammalian central and peripheral cholinergic nervous systems. However, it is also widely expressed in non-neuronal animal tissues as well as in plants, fungi, and bacteria, where it is likely involved in the transport of water, electrolytes and nutrients 48 . With the changes in levels of methanogens, based on histopathology, we found that the colon and brain was damaged; neurotransmitter levels were also changed.

Mycobiota are crucial for human health 49 . Surprisingly, a small number of species can trigger huge changes in the human body 50 . Dysbiosis of and invasion by mycobiota can cause disease in different parts of the body 51, 52 . Meanwhile, the body also produces corresponding immune changes upon mycobiota infection 53 . Several recent studies have made a connection between intestinal mycobiota and the human immune system 52, 53 . In HSHFdiet mice, 10 species of Eukaryota were detected. For example, the relative abundance of Kluyveromyces lactis, Leishmania major, Saccharomycetales species and Theileria orientalis was different (p < 0.05) (Fig. S4B ).

C. albicans is a very common yeast found in the gastrointestinal tract and oral cavity, and can occasionally cause oral-cavity ulcers. C. albicans can pass readily through the blood-brain barrier, where it can cause asymptomatic infection in the cortex, form fungal-induced glial granulomas, and cause short-term memory disorders 54 . Candida dubliniensis is an opportunistic fungal pathogen. Bacher and colleagues revealed that human immunity based on T-helper 17 cells against fungi is reliant on cross-reactivity against C. albicans 53 . Witchley and coworkers showed that programs of C. albicans morphogenesis control the balance between gut commensalism and invasive infection 55 . Neuronal ACh and non-neuronal ACh have been demonstrated to modulate the inflammatory response: ACh protects against C. albicans infection by inhibiting biofilm formation and promoting hemocyte function in a model of Galleria mellonella infection 56 . Hence, C. albicans and C. dubliniensis seem to have important roles in host immunity. In comparison with normal groups, the relative abundance of C. dubliniensis was increased (Fig. S4B) , and levels of proinflammatory mediators were upregulated, upon consumption of the HSHF diet (Fig. 3) . Also, in C. albicans-and K. pneumoniae-treated C57 mice, expression of AChE, AMP, CHRNA1, CHRNB1, and GAD65 were changed in the intestine (Fig. 9A ) and brain (Fig. 9A) , especially in the C. albicans-treated group in the intestine. These data indicated that changes in a few microorganisms could influence microbiota-mediated compounds (including neurotransmitters).

Cooperation within the GM is complex and important for human health, such as mood and bodyweight control. In most cases, the intestinal flora cooperate and influence each other. Cooperative phenotypes are thought to be at the core of microbial-community functions, including through quorum sensing, biofilm formation, and antibiotic resistance 57, 58 . We found that with an increase in some Saccharomyces species who converted sugars to carbon dioxide and ethanol, the abundance of some methanogens was increased. The reason may have been because the metabolism of alcohol imbalances the proportion of NADH to NAD, galactose tolerance, TG synthesis, and lipid peroxidation (Figs. 3 and 7) . Xu and colleagues revealed that chronic exposure to alcoholinduced GM dysbiosis, and was correlated with neuropsychic behaviors 59 . Furthermore, a metabolite from Saccharomyces, acetic acid, is the growth substrate of methanogens. The cooperation of these two microorganisms promotes food digestion/absorption, energy storage, and the accumulation of some harmful products 60 . Studies have demonstrated that the PI3K-I/ Akt-mTOR signaling network can regulate anabolic processes such as the synthesis of lipids, fatty acids, and nucleotides, and requires an abundant supply of reducing power in the form of NADPH, and the growth factorstimulated PI3K-Akt-mTORC1 signaling network 61, 62 , and the NADPH oxidase nox can be derived by gut microbiome 63 .

Twenty-two species of viruses were detected (Fig. S4C) 64 . Furthermore, expression of retinoic acid-inducible gene I supported the notion that too much of a HSHF diet affects immunity and increases the risk of virus invasion (p < 0.05) (Fig. S4D) . Studies have revealed that Herpesvirus species are associated with AD 65 . Therefore, gut dysbacteriosis appears to be a critical factor in inducing changes in the circRNA-expression profile in the brain. Our results also suggest that the stability of circRNAs in terms of structure and quantity is needed for health (i.e., circRNAs may have non-negligible roles in physiological function/regulation in the organism). Also, studies have revealed that the bacteria, diet, and host genes are related to the invasion, infection degree, and drug resistance of viruses 66, 67 .

A total of 622 species of bacteria were detected. The relative abundance of species was changed obviously in the model group (Fig. S5) . Changes were observed in the abundance of Bifidobacterium, Lactobacillus, Bacteroides, Prevotella, Streptococcus and Clostridium species ( Figure S5A-D) . With an increase of Lactobacillus prophage Lj771 abundance (Fig. S4C) , the abundance of most of the species of Lactobacillus was reduced significantly, including that of L. acidipiscis, L. agilis, L. allii, L. amylophilus, L. amylovorus, L. brevis, L. coryniformis, L. fermentum, L. gasseri, L. ginsenosidimutans, L. helveticus, L. jensenii, L. johnsonii, L. kefiranofaciens, L. paraplantarum, L. pentosus, L. plantarum, L. reuteri, L. rhamnosus, L. Ruminis and L. salivarius, which were found only in the control group (Fig. S4A-D) . Simultaneously, the abundance of Bifidobacterium species, including that of B. adolescentis, B. asteroids, B. bifidum, B. breve, and B. pseudocatenulatum, was also inhibited. In contrast, the abundance of most Bacteroides species was increased, including B. caccae, B. caecimuris, B. cellulosilyticus, B. dorei, B. fragilis, B. helcogenes, B. ovatus, B. salanitronis, B. Thetaiotaomicron and B. vulgatus. Phages may have been the reason why the ratio of Bifidobacterium species and Lactobacillus species was imbalanced in mice fed the HSHF diet (Fig. 6E) . The abundance of Bifidobacteriaceae species and Lactobacillaceae species is crucial for homeostatic balance in the intestine. The relative abundance of Lactobacillaceae species was increased and that of Bifidobacteriaceae species was decreased in APP/PS1 mice (Fig. 6E) . These data provide evidence that probiotics are beneficial for health under certain physiological conditions.

The abundance of members of the family Lachnospiraceae (Fig. S5) , including Butyrivibrio hungatei, Anaerostipes hadrus, Butyrivibrio proteoclasticus, Blautia sp. YL58, Herbinix luporum, Blautia hansenii, Roseburia hominis, Lachnoclostridium phocaeense, [Clostridium] bolteae, phytofermentans, Lachnoclostridium Lachnoclostridium sp. YL32 and Pediococcus pentosaceus, was reduced in the intestinal contents of mice fed the HSHF diet ( Fig. S5A-D) . These Lachnospiraceae members encode a composite inositol catabolism-butyrate biosynthesis pathway, the presence of which is associated with a lower risk of host metabolic disease 68 . Members of the Lachnospiraceae family are among the main producers of short-chain fatty acids. Different taxa of Lachnospiraceae are also associated with different intra-and extraintestinal diseases 69 . Supplement the butyrate-producing Lachnospiraceae is beneficial for the intestinal barrier 70 , probiotic for treating stress-induced visceral hypersensitivity 71 , the butyrate-producing species R. inulinivorans includes strains able to grow on inulin and FOS in pure culture 72 . Furuya and colleagues revealed that Blautia hansenii can hydrolyze glucosylceramide to ceramide in plants 73 . Previously, we showed that circNF1-419 could regulate the synthesis and metabolism of ceramides, and the changed ceramides were enriched in three signaling pathways (neurotrophin, sphingolipid, and adipocytokine) 41 . We could conclude that the couples of Blautia and ceramide, butyrate biosynthesis bacteria, and butyrate might be one of interactions between circRNA and microbiome-gut-brain axis, but we still need much more evidence.

Peptostreptococcaceae, a family within the order Clostridiales (Fig. S5) , includes the genera Peptostreptococcus, Acetoanaerobium, Proteocatella, Sporacetigenium, Filifactor, and Tepidibacter. The genera Acetoanaerobium, Sporacetigenium and Proteocatella are monospecific. Representatives of the family have different cell morphology, which varies from cocci to rods and filaments.

Species of Filifactor, Proteocatella, Sporacetigenium, and Tepidibacter form endospores. All members of the family are anaerobes with have a fermentative type of metabolism. The genus Tepidibacter contains moderately thermophilic species. Members of Peptostreptococcaceae are found in different habitats, including the human body, manure, soil, and sediments. Species of Peptostreptococcus and Filifactor are components of the human oral microbiome 74 . Metagenomics sequencing showed that the relative abundance of Peptoclostridium acidaminophilum and Clostridioides difficile was increased in HSHF diet-fed mice. P. acidaminophilum can ferment amino acids 75 . C. difficile is a spore-forming, anaerobic, intestinal pathogen that causes severe diarrhea that can lead to death 76 . P. acidaminophilum and C. difficile had a negative correlation with expression of Tpm3 and Dusp6, and a positive correlation with expression of 5-HT and 5-HIAA, but the mechanism of action needs further analysis. Bacteria can produce a range of major neurotransmitters, and substantial evidence has accumulated around the microbiotamediated influence of those compounds 13 . However, the microbiota can also influence levels of neurotransmitters, including histamine, gasotransmitters, neuropeptides, steroids, and endocannabinoids.

We demonstrated again that consumption of a HSHF diet-induced dysbacteriosis, damaged the intestinal tract, and changed the neurotransmitter metabolism in the intestine and brain. Our new findings were that consumption of a HSHF diet caused changes in brain function and circRNA profiles. Additional experiments found that the GM byproduct of TMAO could degrade some circRNAs, and the basal level of the GM decided the conversion rate of choline to TMAO. A change in the abundance of C. albicans and/or K. pneumoniae could influence the cholinergic system. These findings demonstrate a new link between metabolism, brain circRNAs, and the GM, enlarge the microbiome-transcriptome linkage library and provide more information on the gut-brain axis.

There are several strengths with this study, but there are also limitations. For example, the study lacks a rescue test with large number of single bacteria, so we have not fully identified the one-on-one communication/interaction mechanism between neurotransmitter, circRNA and single bacteria, and therefore could not find the signaliccgeted strains in germ-free mice with multiomics studies to reveal the interaction of neurotransmitter and circRNAs on the microbiome-gut-brain axis, or knock-out mice are needed in the future for positive validations. 

A forward chemical genetic screen reveals gut microbiota metabolites that modulate host physiology

Gut microbiota, metabolites and host immunity

The microbiotagut-brain axis in obesity

Microbiome, probiotics and neurodegenerative diseases: deciphering the gut brain axis

The neuropharmacology of butyrate: The bread and butter of the microbiota-gut-brain axis?

Neuropods übermitteln Informationen über Nahrungsmittel im Darm über vagale Neuronen in Millisekunden and as Gehirn

Altered microbiomes distinguish Alzheimer's disease from amnestic mild cognitive impairment and health in a Chinese cohort

Mild cognitive impairment has similar alterations as Alzheimer's disease in gut microbiota

Gut bacterial tyrosine decarboxylases restrict levels of levodopa in the treatment of Parkinson's disease

Gut microbes metabolize Parkinson's disease drug

Diet, genetics, and the gut microbiome drive dynamic changes in plasma metabolites

Hyperglycemia drives intestinal barrier dysfunction and risk for enteric infection

Neurotransmitter modulation by the gut microbiota

The gut microbiota mediates the anti-seizure effects of the ketogenic diet

GABA-modulating bacteria of the human gut microbiota

Identification of attractive odorants released by preferred bacterial food found in the natural habitats of C. elegans

A neurotransmitter produced by gut bacteria modulates host sensory behaviour

The gut microbiota and Alzheimer's disease

A gut-brain neural circuit for nutrient sensory transduction

The gut microbiome from patients with schizophrenia modulates the glutamate-glutamine-GABA cycle and schizophrenia-relevant behaviors in mice

Gut microbiome and serum metabolome alterations in obesity and after weight-loss intervention

Small intestinal microbial dysbiosis underlies symptoms associated with functional gastrointestinal disorders

Cholinergic neurons in the dorsomedial hypothalamus regulate food intake

The cholinergic system in the pathophysiology and treatment of Alzheimer's disease

Modulation of gut microbiota composition by serotonin signaling influences intestinal immune response and susceptibility to colitis

Choline acetyltransferase-expressing T cells are required to control chronic viral infection

Gut microbial metabolite TMAO enhances platelet hyperreactivity and thrombosis risk

The gut microbiota-derived metabolite trimethylamine Noxide is elevated in Alzheimer's disease

A subnanomolar concentration of Pituitary Adenylate Cyclase-Activating Polypeptide (PACAP) pre-synaptically modulates glutamatergic transmission in the rat hippocampus acting through acetylcholine

Trimethylamine N-oxide: breathe new life

Gut microbiota, inflammation, and probiotics on neural development in autism spectrum disorder

Anxiety, depression, and the microbiome: a role for gut peptides

Altered gut microbiota and mucosal immunity in patients with schizophrenia

The impact of starvation on the microbiome and gut-brain interaction in anorexia nervosa

The intestinal barrier in multiple sclerosis: implications for pathophysiology and therapeutics

The gut microbiome and epilepsy

Gut microbiota regulate motor deficits and neuroinflammation in a model of Parkinson's disease

Circular RNAs are a large class of animal RNAs with regulatory potency

An insulin-sensitive circular RNA that regulates lifespan in drosophila

Circular RNAs: splicing's enigma variations

CircNF1-419 improves the gut microbiome structure and function in AD-like mice

Brain monoamine oxidase B and A in human parkinsonian dopamine deficiency disorders

Methanogens as emerging pathogens in anaerobic abscesses

Evidence of archaeal methanogens in brain abscess

Breath methane concentrations and markers of obesity in patients with functional gastrointestinal disorders

Inhaled methane protects rats against neurological dysfunction induced by cerebral ischemia and reperfusion injury: PI3K/Akt/HO-1 pathway involved

Microbiome and its role in irritable bowel syndrome

Expression of acetylcholine (ACh) and ACh-synthesizing activity in Archaea

The mycobiota of the human body: a spark can start a prairie fire

Early gut mycobiota and mother-offspring transfer

Gut mycobiota in immunity and inflammatory disease

The gut mycobiota: insights into analysis, environmental interactions and role in gastrointestinal diseases

Human anti-fungal Th17 immunity and pathology rely on cross-reactivity against Candida albicans

Binary interactions between the yeast Candida albicans and two gut-associated Bacteroides species

Candida albicans morphogenesis programs control the balance between gut commensalism and invasive infection

Repurposing pilocarpine hydrochloride for treatment of Candida albicans infections

Impacts of microbiome metabolites on immune regulation and autoimmunity

Microbiota and the social brain

Chronic alcohol exposure induced gut microbiota dysbiosis and its correlations with neuropsychic behaviors and brain BDNF/Gabra1 changes in mice

The acetate uptake transporter family motif "NPAPLGL(M/S)" is essential for substrate uptake

Sestrin2 is a leucine sensor for the mTORC1 pathway

Mechanisms of mTORC1 activation by RHEB and inhibition by PRAS40

Microbiota-derived lactate activates production of reactive oxygen species by the intestinal NADPH oxidase nox and shortens Drosophila lifespan

Structure and degradation of circular RNAs regulate PKR activation in innate immunity

Alzheimer's disease-associated beta-amyloid is rapidly seeded by herpesviridae to protect against brain infection

Gut microbiota and Covid-19-possible link and implications

The human gut virome is highly diverse, stable, and individual specific

Higher risk of stroke is correlated with increased opportunistic pathogen load and reduced levels of butyrate-producing bacteria in the gut

The controversial role of human gut Lachnospiraceae. Microorganisms

Construction of a model culture system of human colonic microbiota to detect decreasedlachnospiraceae abundance and butyrogenesis in the feces of ulcerative colitis patients

Beneficial effect of butyrate-producing Lachnospiraceae on stress-induced visceral hypersensitivity in rats

Diversity of human colonic butyrate-producing bacteria revealed by analysis of the butyryl-CoA:acetate CoA-transferase gene

Isolation of a novel bacterium, Blautia glucerasei sp. nov., hydrolyzing plant glucosylceramide to ceramide

Phylogenomic analysis of the family Peptostreptococcaceae (Clostridium cluster XI) and proposal for reclassification of Clostridium litorale (Fendrich et al. 1991) and Eubacterium acidaminophilum (Zindel et al. 1989) as Peptoclostridium litorale gen. nov. comb. nov. and Peptoclostridium acidaminophilum comb. nov

Identifying anaerobic amino acids degraders through the comparison of short-term and long-term enrichments

We would like to thank Zhang Maolei for helpful advises in the preparation of this study. Sequencing services were provided by Guangzhou Geneseed Biotech Co., Ltd. Guangzhou, and Biomarker Technologies Co., Ltd. Beijing, China. The present work was supported by the financial support from the National Natural Science Foundation of China (81701086) 

Data collection and sharing for this project was funded by the National Natural Science Foundation of China.

The animal protocols used in this work were approved by the Institutional Animal Care and Use committee of the Center of Laboratory Animals of the Guangdong Institute of Microbiology.

The authors declare no competing interests.Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

The online version contains supplementary material available at https://doi.org/10.1038/s41398-021-01443-2.