key: cord-0958384-mubw907e authors: Abedpoor, Navid; Taghian, Farzaneh; Hajibabaie, Fatemeh title: Cross Brain–Gut Analysis Highlighted Hub Genes and LncRNA Networks Differentially Modified During Leucine Consumption and Endurance Exercise in Mice with Depression-Like Behaviors date: 2022-04-27 journal: Mol Neurobiol DOI: 10.1007/s12035-022-02835-1 sha: e019abde4224f7aef6a9c8cde0097c19509732aa doc_id: 958384 cord_uid: mubw907e Depression is a frequent mood disorder that might impair the brain-gut axis. In this study, we divided 30 mice into five groups: untreated mice, mice with depression-like behaviors, mice with depression-like behaviors treated with consumed leucine, mice with depression-like behaviors treated with exercise training, mice with depression-like behaviors treated with exercise training along with consumed leucine. According to artificial intelligence biological analysis, we found some mediators such as lncRNAs profile and Kdr/Vegfα/Pten/Bdnf interactions network in the hippocampus region and ileum tissue which could be decisive molecules in the brain-gut axis. Moreover, KDR as a principal cutpoint protein in the network was identified as the pharmaceutical approach for major depressive ameliorating based on pharmacophore modeling and molecular docking outcomes. Furthermore, we indicated that the mRNA and protein level of the Pten enhanced and Vegfα/Kdr/Bdnf mRNAs, as well as the protein level of KDR, decreased in mice with depression-like behaviors. Moreover, exercise and leucine ameliorated the brain-gut axis in mice with depression-like behaviors. Exercise and leucine regulated the lncRNAs network in the hippocampus and ileum of mice with depression-like behaviors. We suggest that the lncRNAs profiles could be considered as diagnosis and prognosis biomarkers, and exercise + leucine might be a practical approach to improve depression. In developing and developed countries, psychological disorders have imposed heavy economic and social burdens on health care systems [1, 2] . Depression is a frequent mood disorder that could lead to increased cerebrovascular disorders, cardiac diseases, gastrointestinal diseases, disruption of microbiota in the ileum, other medical causes of mortality, and the risk of anxiety and depression speculation during the coronavirus pandemic [3] [4] [5] . Although the pathomechanism of depression is not entirely elucidated, an artificial intelligence survey indicated that several signaling pathways such as neuroinflammation, neurotransmitter functions, phosphatidylinositol signaling system, gap junction, muscle contraction, insulin signaling pathway, and dopaminergic synapse could be involved in this pathogenesis [5, 6] . Growing evidence has indicated various neurohumoral transmissions enabling the communication between the gut and central nervous system (CNS). Also, there are mutual interactions between the mucosal immune system and enteric nervous systems, and this cross-talk is called the gut-brain axis [7] [8] [9] . The gut-brain axis could modulate behavior, mucosal immune system, and emotion [10] . Hence, the gut function is essential for healthy brain function [11] . Based on the evidence, the gut microbiota is a critical environmental compound that affects gastrointestinal system function [12] . Notably, microflora could regulate human health and illness by bidirectional communication between the microbiota and the host [13] . Interestingly, the consumption of antibiotics could deplete the gut microbiota, disarrange intestinal epithelium, disrupt the expression of neuromodulators, and lead to cognitive disability in mice [14] . The gut microbiota might interact with genetic variables to co-regulate illness symptoms in genetically vulnerable people [7] . Furthermore, immense evidence has revealed that depression impairs the intestinal bacterial community, increasing inflammation in the intestinal epithelium and disrupting the gut-brain axis [15] . Thereby, the gut function is essential for healthy brain function [16, 17] . However, the mechanistic and molecular pathway is not elucidated. Recent evidence has indicated that the gut-brain axis had a critical role in digestion and absorbance via regulating microbiota in the intestinal epithelium [18] . On the other hand, Darren W. Roddy et al. revealed that the hippocampus area size in depressive conditions was reduced compared with healthy individuals. Moreover, they have shown a direct association between hippocampus size and depression stages [19] . Hence, the hippocampus might play a decisive role in the brain-gut axis. There is no substantial treatment strategy to halt the progress of depression. Notably, only 60% of depression treatments reported improved symptoms [18] . Therefore, considering effective treatment accompanied without any adverse effect could provide insight into diminishing this disorder and ameliorating the ecosystem in the gastrointestinal system. In addition, based on the epidemiologic evidence, environmental agents might play crucial roles in depression [20] . Therefore, current studies focus on physical activity and dietary interventions in managing depressive symptoms [21] . Moreover, it was well-established that sufficient physical activity and nutrients could be requirements for producing neurotransmitters [22] . Physical activity has been reported to have a relationship with reduced risk of depression [23, 24] . Exercise is considered a non-pharmacological intervention that could be a valuable tool to halt and manage the onset of depression [25] . In addition, physical activity, possibly through positive effects on self-efficacy, coping strategies, quality of life, and body image, caused decreasing mental disorders [26] . Interestingly, aromatic amino acids, including tyrosine, phenylalanine, and tryptophan, could benefit by producing neurotransmitters in therapeutic depression [27, 28] . Based on recent evidence, leucine competed with aromatic amino acids for transport across the blood-brain barrier (BBB) and also could trigger a delay in serotonin secretion, suppression of central fatigue, and decreased level concentration of aromatic amino acids, which led to the possibility of reduced neurotransmitters related to aromatic amino acids. However, although the concentration of leucine significantly declined in depressive conditions, the pathomechanism of this phenomenon is not clear [29] . Recently, evidence has suggested that compared with protein-coding genes in body tissues and fluids, lncRNAs (long non-coding RNAs), roughly more than 200 nucleotides, are feasible to be identified via various techniques, have higher stability in body tissues and fluids, and also have specific expression patterns in tissue [7] . Furthermore, growing evidence has emerged showing that lncRNAs could modulate the expression level of genes through different molecular mechanisms [30] . Therefore, there are several possibilities in lncRNA function associated with the various diseases such as lncRNAs that might influence the neighboring genes that regulate the expression level of networks and could be affected by SNPs [31] . In this artificial system biology study, the principal aim was to evaluate the decisive lncRNAs and hub genes in depression and also demonstrate that the route network between the brain and intestine might improve with regular exercise and consumption of leucine. To predict the master genes involved in the brain-gut axis, we analyzed the gene expressions between two microarray datasets (GSE151807 and GSE171275) with major depressive disorders in various brain regions [32] . Moreover, we engaged GSE64004 in ileum tissue suffering from depression and stress using the Bioconductor package in R statistical programming language. These datasets were analyzed using the Mouse430_2 and Mouse Gene 2.1 Affymetrix Array platform, collaborating with Oligo and Limma packages. Significant differential gene expressions obtained from R statistical analysis of microarray datasets (brain and ileum) separately were illustrated on the HeatMap diagram (P < 0.001). We limited these genes based on P < 0.01 and logFC ± 0.2 cut off to determine essential hub genes among the list of the genes collected from the dataset analysis. Hence, we engaged 1308 and 1003 genes as significant differential gene expressions in the hippocampus and ileum in major depressive conditions. The protein-protein interaction networks constituted were drawn based on the medium confidence threshold in the STRING 11.0 database [33] . Furthermore, calculating essential network parameters such as degree, betweenness, and closeness was utilized by Cytoscape 3.6.0 software [34] . By applying network parameters, we obtained 120 hub genes from ileum tissue analysis based on betweenness: 0.006, closeness: 0.259, and degree: 9. Moreover, 112 hub nodes were detected from analysis of significant gene expressions of hippocampus tissue by considering betweenness:0.012, closeness: 0.175, and degree: 6. To remarkable pivotal molecular and cellular signaling pathways involved in the mood disorder pathomechanism, molecular networks were investigated in specific biomedical servers, a KEGG pathway analysis as a sub-browser of the KOBAS-i database, by P < 0.05 [35] . To determine the common genes between the two hub node lists, we used the Venny 2.1.0 tool [36] . Based on the interactive tool results, 46 overlapped hub genes were concluded. A protein-protein interaction network including 46 overlapped hub genes based on network parameters and eigenvector centrality revealed that four genes, VEGF-α (vascular endothelial growth factor A), KDR (vascular endothelial growth factor receptor 2 or kinase insert domain receptor), PTEN (phosphatase and tensin homolog), and BDNF (brain-derived neurotrophic factor), having the highest degree and most betweenness, could be affected by other hub gene expression patterns under the influencing of endurance exercise and leucine consumption. Enrichment of 46 hub proteins based on KEGG, Panther [36] , and KOBAS-i algorithms highlighted significant molecular signaling pathways while considering P < 0.05. To determine significant lncRNAs with differential transcription profiles in depression, we browsed the gene expression omnibus database and selected GSE189233 for analysis. On the other hand, in this in silico study, we predicted long non-coding RNAs related to VEGFα and PTEN based on diseases and target genes by LncTard [37] , LncRNADisease, mammalian ncRNA-disease repository (MNDR), and LncBook servers. Considering the expression levels of predicted lncRNAs in the ileum tissue and hippocampus region in the LncSEA database, GAS5 (growth arrest specific 5), HOTAIR (HOX transcript antisense RNA), MEG3 (maternally expressed 3), and TUG1 (taurine upregulated 1) were selected. Next, we showed that four long non-coding RNAs had an overlap between PTEN and VEGF-α. In Fig. 1 , we depicted the workflow of the study. Based on evidence and visualization of networks between selected hub genes for experimental assays, VEGFR2 (KDR) as a pivotal receptor has been involved in initiating molecular signaling pathways. Also, we realized that VEGFR2, as a potential cut point with a pivotal role in molecular and cellular pathways, could be effective in other expression profiles of genes involved in the network. Therefore, we obtained VEGFR2's ligands from Binding Database for pharmacophore model design. Moreover, the three-dimensional structure of hub proteins was downloaded from the Protein Data Bank server [38] , and the structure-data file (SDF) of leucine was stored from the PubChem database [39] . We investigated leucine's pharmacokinetic and pharmacodynamic features based on the Swiss ADME database [40] . We indicated that leucine was a necessary complement with high absorbance in the gastrointestinal system, and with In-silico analysis Interaction between LncRNA and mRNA in depression condition Fig. 1 In silico machine flowchart of the study lipophilicity feature Log P = 1.14 and topological polar surface area (TPSA) = 63.32 Å 2 , it could cross the blood-brain barrier (BBB). According to these results, we hypothesized that leucine consumption could affect the gene expression profile of the ileum and the brain. Thus, we designed a suitable pharmacophore model based on active ligands obtaining the binding database library. First, the pharmacophore model design was performed based on the OPLS-2005 force field option by Ligprep and Phase-Pharma applications with 4 to 5 features and 80% of active ligands (IC50 < 1) in the Schrödinger server [41] . Next, we tried matching the leucine structure as a ligand on the pharmacophore model built. Moreover, we computed the leucine's binding affinity and stability on VEGFR2, VEGF-α, PTEN, and BDNF by molecular docking method in PyRx software [42] . This calculation made a contract based on binding affinity < − 5 and root-mean-square deviation of atomic positions (RMSD) < 2 as suitable binding energy between macromolecule and small molecule. Then, optimization and preparation of three-dimensional structures of macromolecules were performed following extra chains ablation and removing ligands and non-complex compounds, based on dock prep tools in Chimera 1.8.1 software [43] . Finally, molecular docking of VEGFR2, VEGF-α, BDNF, and PTEN was applied in the search space box with central dimensions: X: 65 All procedures were conducted and approved following the Research Ethics Committees of Islamic Azad University Isfahan (Khorasgan) Branch (IR.IAU.KHUISF. REC.1400.160). The 6-week-old male C57BL/6 mice were provided and housed at the animal house of the Isfahan (Khorasgan) Branch, Islamic Azad University, Isfahan, Iran. Mice were kept under a standard condition of 24 ± 3 °C temperature, 55 to 60% humidity, and 12-h light and dark cycles and fed normal diets and ad libitum water. Moreover, all animals fed a normal diet (contained 15% (w/w) fat, 58% (w/w) carbohydrate, 27% (w/w) protein, ad libitum, 3.2 kcal/1 g). Mice were adapted for 1 week prior to the start of experiments. Mice with an approximate weight of 23 ± 2 g were randomly assigned into five groups (n = 6 per group): (control) untreated mice, (depression) mice with depression-like behaviors, (Leu) mice with depression and consumed the leucine complement, (EXr) mice with depression with exercise training, (Leu + EXr) mice with depression with exercise training along with consumed leucine complement. The calorie intake, weight, and water drinking were monitored weekly. Mice fasted for 6 h before euthanasia, after 8 weeks of consuming the leucine complement and exercise training. Mice were euthanized by administration of ketamine (80 mg/kg body weight per mouse) and xylazine (10 mg/kg body weight per mouse). This study induced depression in all groups except untreated mice (control). Based on previous evidence, mice were depressed following protocols for 14 days [44, 45] . Subsequently, we divided depression mice into four groups, as mentioned above. First, mice were exposed to the electric tail and foot shock for 1 s with repetition ten times (intensity, 0.5 mA), deprivation of food, and fixation randomly conducted for 14 days [46, 47] . After that, the behavior test, including the open-field test, elevated plus maze test, tail suspension test, and social interaction test, was conducted to validate depression features. Then, we assessed the effect of the exercise training and consumed the leucine complement on the mice with depression-like behaviors. The open-field test was conducted on days 16 and 17. Moreover, the resting times and locomotor activity were measured. The open-field area was (40 × 40 × 40 cm), divided into 16 squares. Each mouse was placed in the central area and monitored for 1-h sessions. Distance moved (cm), duration time (s), and rest times (s) were scored and calculated [44, 45] . For evaluating the levels of depression and anxiety, this test was used. This device contained four arms with two opposite closed arms (35 × 10 × 1 cm) and open arms (35 × 10 × 1 cm). Each mouse was placed in the central area. In this study, the behavior was monitored for 1 h. In addition, the time spent in open and closed arms was recorded, and the number of entries into each arm was recorded. Enhancing the number of entrances into the closed arm and time spent in the closed arm as described was considered to measure depression and anxiety [45] . We evaluated the social interaction using an open-field apparatus. Moreover, we used two plastic chambers placed on both sides of the box. In one of them, we placed a wild-type mouse (as a target) with no prior contact with the depressed mouse, and another chamber was empty (no target). Then, we placed the depressed mouse on the middle open-field apparatus and allowed it to freely explore in the box for 15 min (we considered 10 min a habituation trial). We calculated the amount of time spent around each chamber (wildtype mouse or empty) in this test. Total distance moved (cm) with no target, time spent in interaction zone (s) with no target, total distance traveled (cm) with the target, and time spent in interaction zone (s) with the target were measured and scored [45, 48] . The mouse was stuck approximately 1 cm by an attached adhesive tape on its tail, and the behavior was assessed for 6 min. In this test, mice were allowed to adapt for 2 min. Total freezing was recorded as a depression-like behavior [45] . Mice were fed with free access to standard food (contained 15% (w/w) fat, 58% (w/w) carbohydrate, 27% (w/w) protein, ad libitum, 3.2 kcal/1 g) and tap water. In addition, l-leucine (Catalogue Number: 105360, Merck) was injected (50 mg/ kg, intraperitoneal injection) once per day for 8 weeks [45] . Endurance exercise training (EXr) was conducted on a motorized treadmill. The intense exercise was moderatehigh on the treadmill for 2 months (6 days/week). After that, exercise duration and running speed gradually increased to reach ~ 75% VO 2 max (32 m/min) for 45 min. Moreover, the treadmill slope was considered to be 0% [49, 50] . Total RNA was extracted from the hippocampus region and ileum tissue using TRIzol reagent (Sigma, USA). According to the manufacturer's instructions, cDNA synthesis was conducted with 1 μg of total RNA using a cDNA synthesis kit (TaKaRa, Japan). The qRT-PCR was conducted with CYBR Green (TaKaRa, Japan) using Corbet rotor gene 6000 (Qiagen, Australia). Detection of gene expression was evaluated according to the 2 −ΔΔCT method. The relative expression of genes was calculated based on glyceraldehyde-3-phosphate dehydrogenase (Gapdh) expression levels. Primer sequences were ordered through the Micro-Gen company (South Korea), and their sequences are listed in Table 1 . Hippocampus was lysed via TRI reagent (Thermo Scientific, 15,596-018, USA). Moreover, we loaded 30 μg of sample protein for the SDS-PAGE, and then protein bands were transferred to the PVDF membrane (Bio-Rad, 162-0176, USA). The primary antibodies used were rabbit KDR antibody (Abcam, ab11939, USA), anti-PTEN (ab32199, Abcam), and mouse anti-GAPDH antibody (Santa Cruz, sc-32233). Membranes were incubated with primary TTG CTC TCC TTA TCA TCC ATT AGT GCC TCT CAG TCC 132 51 GAS5 AGT GTG GAC CTC TGT GAT CTC GTA AGC CTT CAT CCT CCT 78 51 TUG1 CTC TGG AGG TGG ACG TTT TGT GTG AGG TCG TGT CTC TCT TTT CTC 72 51 GAPDH TGC CGC CTG GAG AAACC TGA AGT CGC AGG AGA CAA CC 121 58..6 antibodies for 2 h at room temperature, then for 1 h at room temperature with an appropriate secondary antibody. We detected each protein bandwidth Amersham ECL Advance Western Blotting Detection Kit (GE Healthcare, USA). Furthermore, the intensity of bands was calculated by ImageJ software. The sample size was estimated based on 80% power and an alpha level of 0.05. Statistical analysis was conducted using GraphPad Prism (Version 9; GraphPad Software). The Shapiro-Wilk test was used for normalizing distribution, and variables were normally distributed. Data were analyzed by one-way analysis of variance (ANOVA) with Tukey's post hoc test due to multiple comparisons. Differences at P < 0.01 were considered to be significant. Moreover, data were indicated as the mean ± SD. Based on R statistical analysis, we collected 1308 genes with differential expression in mice with major depressive disorders in the hippocampus region compared to normal tissue. Moreover, 1003 differential gene expressions were identified in the ileum tissue of depressive mice compared to normal mice. We showed these differential gene expressions via the HeatMap diagram in Fig. 2A and B, with P < 0.001. According to the optical network's parameters, such as betweenness, degree, and closeness, we calculated the centrality parameters of hub genes. Next, we designed the protein-protein interaction network to consider the nodes' degree, modularity, and betweenness centrality (Fig. 2C, D) Hippocampus, significant genes with P-value<.001 sub-browser of the KOBAS-i server for determining significantly molecular signaling pathways with P < 0.05 marked several signaling and molecular mechanisms such as dopaminergic synapse, platelet activation, focal adhesion, Rap1 signaling pathway, Notch signaling pathway, cell cycle, and AGE-RAGE signaling pathway involved in the hippocampus region with depressive disorders (Fig. 3A) . On the other hand, the bar graph of ileum's hub genes enrichment manifested cytokine-cytokine receptor interaction signaling pathway, immunodeficiency, TNF signaling pathway, IL-17 pathway, Toll-Like receptor, and cell differentiation as crucial molecular signaling pathways involved in the pathomechanism of the ileum tissue in depressive conditions (Fig. 3B) . Moreover, based on the Venn diagram results, 46 common hub genes were distributed between the hippocampus region and ileum tissue (Fig. 4A) . Based on the design protein-protein interaction network construction, we indicated that four genes could be pivotal in the brain-gut axis, including Pten, Vegf-α, Kdr, and Bdnf. We selected these common genes via the critical network's highest betweenness and maximum degree. Significant molecular signaling pathways associated with common genes based on the enrichment analysis in the KEGG pathway, Wiki pathway, Reactome, and Panther have marked several pathways such as PI3K-AKT signaling pathway, AGE-RAGE signaling pathway, P53 signaling pathway, TGF-B signaling pathway, MAPK signaling pathway, RAS signaling pathway, and focal adhesion as vital pathways involved in the gut-brain axis (Fig. 4B-D) . Integrative analysis of non-coding RNAs (LncRNAs) and coding RNAs (mRNAs) was performed in neurodevelopmental impairments in the gut dysbacteriosis mice based on GEO2R analysis. To specify significant lncRNAs, P < 0.05 was applied, which highlighted 358 lncRNAs (195 overexpressed and 163 down expressed). Based on non-coding microarray analysis, GAS5 with downregulation pattern expression was confirmed in depression-like behavior in mice compared to healthy mice. Moreover, interaction types of lncRNAs-mRNAs were browsed in the lncRNAs databases and admitted by text mining ( Table 2) . On the other hand, the prediction of non-coding RNAs, especially long non-coding RNAs (lncRNAs) as regulators of the gene's expression, was performed based on genes and diseases. First, we prepared a list of lncRNAs related to selected genes (Table 3) . Then, we confirmed the expression profiles of lncRNAs in the hippocampus region and ileum tissue based on the database survey and in silico analysis. Eventually, we specified four lncRNAs for expression measurement in two tissues (hippocampus and ileum) after regular exercise and consumption of leucine in depressive mice. Furthermore, we drew a comprehensive network between hub genes and selected lncRNAs that showed that four lncRNAs, including MEG3, HOTAIR, GAS5, and TUG1, are common between Pten and Vegf-α (Fig. 5A, B) . Visualizing the network of selected genes involved in the PI3K-Akt signaling pathway in STRING and data mining led to the selection of KDR as a strong candidate for preventing and treating major depressive disorders and gut mucosal epithelium damage. Drug design with the pharmacophore modeling method was performed based on these data. The AHRRR pharmacophore model, including Acceptor, hydrophobic, and aromatic ring features, with the highest survival score of 5.762731 based on 251 active ligands (IC50 < 1), was designed (Fig. 6A) . In contrast to the designed pharmacophore model, the 3D structure of leucine did not match the model built, but the molecular docking technique with binding affinity − 5.1 kcal/mol and RMSD < 2 confirmed the effect of leucine complement on KDR (Fig. 6B) . Furthermore, the drug design results based on the molecular docking method revealed that leucine with binding ability > − 5 kcal/mol and RMSD threshold < 2 on BDNF, PTEN, and VEGFα proteins, due to the non-suitable binding affinity between leucine and the mentioned proteins, could not be influenced by the function of these proteins (Fig. 6C-E) . This study found that mice with depression-like behaviors significantly gained weight more than did the control group (Fig. 7A , B: P < 0.01). Moreover, we indicated that drinking water consumption and food intake were enhanced in mice with depression-like behaviors (Fig. 7C , D: P < 0.01). Based on these physical features, we could conclude that the excessive food and water intake significantly increased compared with the control group (Fig. 7C , D: P < 0.01). In addition, we exhibited polydipsia and hyperphagia-like symptoms following an increase in weight and drinking water consumption. Here, we evaluated depression levels and the social defeat of the mice. Thus, we found that the distance moved (unite) and movement time (s) of mice with depression-like behaviors were reduced compared with the control group (Fig. 7E , F: P < 0.01). In addition, the rest time of the mice with depression-like behaviors significantly increased with the control group ( Fig. 7G: P < 0.01) . Furthermore, the elevated plus-maze test revealed that the number of entries into the opened and closed arms had decreased and increased, respectively (Fig. 7H , I: P < 0.01). Moreover, the total movement of mice with depressionlike behaviors significantly declined compared with the control group (Fig. 7J: P < 0.01) . Also, the freezing time of the mice with depression-like behaviors was increased compared with the control group ( Fig. 7K: P < 0.01) . Furthermore, the social interaction test demonstrated that the total distance moved (cm) with no target, time spent in interaction zone (s) with no target, total distance traveled (cm) with a target, and time spent in interaction zone (s) with the target were predominantly decreased compared with the control group (Table 4 : P < 0.01). Therefore, depression diminished mobility and induced social deficit-associated behavior based on these data. Based on the in silico machine, we found that Kdr/Vegfα/ Pten/Bdnf predicted a vital role in the brain and gut axis in mice with depression-like behaviors. Moreover, these molecules were common between the hippocampus region and ileum tissue. Therefore, we assessed the expression levels of Kdr/Vegfα/Pten/Bdnf in the hippocampus region and ileum tissue of mice with depression-like behaviors ( Fig. 8A -H: P < 0.01). Based on our results, we indicated that the relative expression of Pten in mice with depression-like behaviors significantly increased compared with the control group in the hippocampus region and ileum tissue (Fig. 8A , B: P < 0.01). In addition, the data indicated that the expression level of the VEGFα/ Kdr/Bdnf in the hippocampus region and ileum tissue was significantly reduced in mice with depression-like behaviors compared with the control group (Fig. 8C -H: P < 0.01). In addition, exercise and leucine consumption decreased the relative expression of the Pten (Fig. 8A , B: P < 0.01) and significantly enhanced the expression level of VEGFα/Kdr/Bdnf in the hippocampus region and ileum tissue (Fig. 8C -H: P < 0.01). On the other hand, based on pharmacophore modeling and in silico analysis, KDR as a principal cutpoint protein in the network was identified. Our data found that KDR was a vital molecule in the Kdr/ Pten/Vegfα/Bdnf network. Hence, we assessed the protein expression of the KDR in the hippocampus region ( Fig. 8I : P < 0.01). In this study, we indicated that the expression protein of KDR was significantly decreased in mice with depression-like behaviors compared with the control group ( Fig. 8I : P < 0.01). Also, exercise and the leucine complement enhanced the protein level of KDR ( Fig. 8I : P < 0.01). Notably, we observed the same pattern between protein and mRNA expression of KDR in the hippocampus region (Fig. 8F , I: P < 0.01). Moreover, the protein level of PTEN significantly increased in mice with depression-like behaviors vs. the control group in the hippocampus region ( Fig. 8J : P < 0.01). Our data revealed that the level of PTEN protein was reduced by exercise and the leucine complement ( Fig. 8J : P < 0.01). Interestingly, the interaction of the leucine complement and exercise significantly decreased the level of PTEN protein ( Fig. 8J : P < 0.01). Based on this data, the leucine complement and physical activity could suppress the expression level of PTEN in the hippocampus region ( Fig. 8J: P < 0.01) . Interestingly, the interaction between exercise and the leucine complement significantly affected the Kdr/Vegfα/ Pten/Bdnf signaling pathway (Fig. 8A -J: P < 0.01). Furthermore, based on these results, we discovered that the interaction between exercise and the leucine complement had a synergistic effect on depression (Fig. 8A -J: P < 0.01). As a result, physical activity and the leucine complement ameliorated the brain and gut axis in mice with depression-like behaviors. To confirm that physical activity and the leucine complement could improve the physical features of mice with depression-like behaviors, we evaluated the social deficit-associated behavior Moreover, the plus-maze test demonstrated that the number of closed and opened arms entries decreased and increased in the Leu, EXr, and Leu + EXr groups, respectively (Fig. 8N, O: P < 0.01) . Also, the whole movement of mice significantly increased in Leu, EXr, and Leu + EXr groups compared with the depression group ( Fig. 8P : P < 0.01). Therefore, we found that the physical activity and leucine complement had a synergetic effect on the activity and sociability of the mice with depression-like behaviors. To validate the function of lncRNA prediction, we evaluated the relative expression of the four lncRNA associated with the Kdr/Vegfα/Pten/Bdnf in the hippocampus region and ileum tissue (Fig. 9A -H: P < 0.01). Our data indicated that the expression level of HOTAIR was increased in the mice with depression-like behaviors in the hippocampus region and ileum tissue (Fig. 9A , B: P < 0.01). Moreover, endurance exercise (EXr group) and leucine consumption (Leu group) were lower than in mice with depression-like behaviors in the hippocampus region and ileum tissue (Fig. 9A , B: P < 0.01). In addition, we found the synergetic effect on the interaction between exercise and the leucine complement compared with other groups (Leu + EXr group). In addition, we demonstrated that the relative expression of the MEG3 and It should be noted that for reducing the error in this study, the data of mice behavioral tests were obtained on triplicate datasets for each sample Based on these data, we could conclude that HOTAIR had positive regulation as well as MEG3 and GAS5 lncR-NAs had negative regulation with the Pten gene. Moreover, MEG3 and GAS5 lncRNAs had positive regulation with the Vegfα gene, and HOTAIR lncRNA had negative regulation with Vegfα in the hippocampus region and ileum tissue. Notably, TUG1 lncRNA was significantly changed in the mice's hippocampus region and ileum tissue with depression-like behaviors compared with the control group (Fig. 9G, H: P < 0.01) . Furthermore, exercise and the leucine complement significantly enhanced the relative expression of the TUG1 in the hippocampus region ( Fig. 9G: P < 0.01) . However, we found different patterns of the TUG1 level in the ileum tissue ( Fig. 9H : P < 0.01). Our data indicated that the expression level of the TUG1 was enhanced in mice with depression-like behaviors in ileum tissue ( Fig. 9H : P < 0.01). Moreover, expression of the TUG1 predominantly declined with exercise and consumption of leucine ( Fig. 9H : P < 0.01). Hence, TUG1 lncRNA had different functions in the hippocampus region and ileum tissue. Moreover, based on these data, we might discover the lncRNA profile, which could play a vital role in depression and dysregulate these lncRNAs causing depression and exercise, and the leucine complement could regulate these lncRNAs and improve the depression pathomechanism. In this study, we addressed the role of the brain-gut axis in mice with depression-like behaviors. Based on the in silico machine, we explored the pivotal hub genes and precise pathomechanism Exercise and leucine complement modulated the lncRNA network. A-H The relative lncRNA expression was measured by qRT-PCR for HOTAIR, MEG, TUG1, and GAS5 transcripts relative to GAPDH in the hippocampus region and ileum tissue (mean ± SD; n = 6). **P < .01, ***P < .001. Data were analyzed by one-way analysis of variance (ANOVA) and Tukey's post hoc test in depression conditions. Moreover, we found communication between network mapping hub genes in the hippocampus region and ileum tissue via artificial intelligence. In addition, systems biology analysis indicated that various signaling and molecular mechanisms, including dopaminergic synapse, platelet activation, focal adhesion, Rap1 signaling pathway, Notch signaling pathway, cell cycle, AGE-RAGE signaling pathway, cytokine-cytokine receptor interaction signaling pathway, immunodeficiency, TNF signaling pathway, IL-17 pathway, Toll-Like receptor, and cell differentiation, are involved in the pathomechanism of the hippocampus region and ileum tissue in depression. Text mining of molecular signaling pathways indicated the role of these pathways in gut-brain pathogenesis based on in silico surveys and bioinformatics outcomes. Dopaminergic neurons, defined as the especially enteric nervous system, with the change of norepinephrine precursor (dopamine) to the active form of norepinephrine, could be regulated by both plexuses of bowel functions in the absence of central nervous system (CNS) input [51] . Moreover, many neurotransmitters launched in the CNS have also been found in the enteric nervous system [51] . On the other hand, platelet activation as the critical process has been associated with brain aging and neurodegenerative pathogenesis. The platelet activation process might trigger brain vascular abnormalities such as microvascular integrity disruption in depression conditions to preserve vascular integrity [52] . In a recent study, Ge Zhang and colleagues illustrated a disposition to platelet hyperactivity and microthrombosis risk in intestinal tissue in inflammatory bowel diseases, especially Crohn's [32] . Furthermore, growing evidence reported that AGE-receptor expression levels in neurodegenerative disorders were upregulated compared to normal subjects' samples [53] . These outcomes obtained from previous studies verified that RAGE is a facilitator for neuropathologic alteration [54, 55] . In the other study, Cristina Luceri et al. mentioned that the expression profile of advanced glycation end-products (AGEs) could increase the inflammation condition of intestinal tissue [56] . These results suggested that the AGEs-RAGE signaling pathway might be critical in inflammatory bowel pathogenesis. In fact, due to the multi-ligand receptor conditions for RAGE implicated in immune activation and inflammatory agent secretion, this signaling is identified as one of the potential therapeutic approaches in enteric inflammation [56] . Since the proliferation-differentiation and cell cycle programming are crucial in body homeostasis, Notch signaling pathways could be an effective cell fate determination in various tissues [57] . Therefore, Notch signaling pathways were disrupted in mice with depression-like behaviors in the hippocampus and gut tissues. As new results, we demonstrated that the Kdr/Vegfα/Pten/Bdnf interaction network could be decisive molecules in the brain-gut axis of mice with depression-like behaviors through bioinformatic analysis. In this present study, we found that the relative expression of Kdr/Vegfα/ Pten/Bdnf was dysregulated in major depressive conditions. Our data indicated that the expression level of the Pten was significantly increased in the brain and ileum tissues. Wang et al. have reported that Pten could play an indispensable role in regulating depression in mice [58] . Furthermore, they indicated that Pten might be a mediator for neuronal atrophy and depression. This data was in agreement with our results. Moreover, Kdr, Vegfα, and Bdnf were reduced in the brain and ileum tissues in mice with depression-like behaviors. Furthermore, we evaluated the profile of lncRNAs that could regulate this interaction network. Based on system biology analysis, we explored that HOTAIR, MEG3, GAS5, and TUG1 lncRNAs were associated with Kdr/Vegfα/Pten/Bdnf interaction network and could be a principal modulator hippocampus region and ileum tissue of the depression condition. Interestingly, these lncRNAs profile could be a precise candidate and strategy to predict the prognosis and diagnosis of depressive disorder. There is currently no effective treatment method for depression. Notably, only 60% of depression therapies effectively alleviated depressive symptoms [18] . Meanwhile, psychotherapy and antidepressant pills have been indicated to be effective in the therapeutic treatment of depression; not all are subjected to these treatments, and the effect sizes of treatments are generally modest [59] . Hence, several strategies are required to target prognosis and diagnosis or even halt psychiatric disorders [60] . One beneficial and convenient approach might consider exercise and physical activity [61] . Accordingly, physical activity and exercise training have been identified as helpful treatment techniques and mind-body practice in complementary and alternative health approaches [62] . Regular exercise training may significantly influence physical performance, cognitive improvement, metabolism, microbiota balance, neurohumoral adaptation, and hormone function on a long-term basis [55, 63] . Moreover, as a physiological and natural complement, leucine is a beneficial substance with various applications in cellular and molecular processes. This natural complement may be complementary and alternative medicine in preventing and alleviating illness hallmarks synchronized with medical prescription. Exercise is considered a nonpharmacological intervention that might be a valuable strategy to halt and manage the onset of depression [25, 26] . Evidence has been reported that exercise reduces the risk of depression [23, 24] . Hu et al. performed a merged systematic review and meta-analyses and indicated that exercise might benefit depression [61] . These results agreed with our data; we showed that physical activity improved the behavior tests related to depression and mobility in mice with depressionlike behaviors. Interestingly, leucine consumption could improve essential neurotransmitters and halt depression hallmarks. Moreover, it was well-established that sufficient physical activity and essential nutrients might be requirements for producing neurotransmitters [22, [27] [28] [29] . In the cross-sectional study on 3175 individuals, consumption of branchedchain amino acids significantly reduced depression [64] . In addition, leucine is a beneficial substance with various applications in cellular and molecular processes. This natural complement may be complementary and alternative medicine in preventing and alleviating illness hallmarks synchronized with prescription medication. According to this evidence, leucine might be an effective and safe complement in treating several neurodegenerative disorders and psychological disorders. In addition, the relative expression of the mRNA-lncRNA network correlation with depression condition was modulated by endurance exercise. Interestingly, based on pharmacophore modeling and molecular docking surveys, we indicated that KDR as the initiator element among four selected hub proteins could be a strong candidate for the depressive therapeutic approach and improve the gut-brain axis function. Although the leucine complement was not aligned with KDR's pharmacophore model, computing binding power illustrated that supplementation with suitable and stable binding affinity on the KDR protein surface positively affected brain and intestine functions. Furthermore, we demonstrated that consumption of the leucine complement regulated the Kdr/Vegfα/Pten/Bdnf interaction network and HOTAIR, MEG3, GAS5, and TUG1 lncRNAs. Sideromenos and colleagues indicated that VEGF could modulate depression-related behavior and be the candidate antidepressant agent [65] . In addition, they found that treatment based on VEGF stimulated the phosphorylation of ERK in the hippocampus. Notably, our study revealed that exercise and leucine complement increased the relative expression of Vegfα mRNA and MEG3 and GAS5 lncRNAs. Based on our data, MEG3 and GAS5 had positive regulation related to the Vegfα gene, and HOTAIR lncRNA negatively correlated with Vegfα in the brain and ileum tissues. Thus, we could conclude that consumption of leucine and endurance exercise might have a synergetic effect on regulating the Kdr/Vegfα/Pten/Bdnf interaction network and lncRNAs. One of the limitations of our study was the lack of data on protein expression of the Vegfα and Bdnf. Moreover, we did not consider the exercise training with various intensities, repetitions, and duration for selecting the most efficient exercise program in depression therapy. We suggest that physical activity and leucine consumption in both genders can investigate in the following studies. Here, according to artificial intelligence, biological analysis, and pharmaceutical approaches for major depressive amelioration, we could conclude that synchronizing leucine consumption and regular endurance exercise might have a synergetic effect on the predicted regulatory network containing hub proteins, Kdr/Vegfα/Pten/Bdnf interaction network, and network long non-coding RNAs in the improved gut-brain axis. Prevalence of common mental health disorders in adults who are high or costly users of healthcare services: protocol for a systematic review and meta-analysis The molecular neurobiology of depression Increases in depression and anxiety symptoms in adolescents and young adults during the COVID-19 pandemic Depression in patient with heart failure: a scoping review Influence of stressor-induced nervous system activation on the intestinal microbiota and the importance for immunomodulation. Microbial endocrinology: the microbiotagut-brain axis in health and disease The gutbrain axis: interactions between enteric microbiota, central and enteric nervous systems Integrative analysis of long non-coding RNAs, messenger RNAs, and microRNAs indicates the neurodevelopmental dysfunction in the hippocampus of gut microbiotadysbiosis mice Genetic studies of major depressive disorder: why are there no GWAS findings, and what can we do about it? Molecular adaptations underlying susceptibility and resistance to social defeat in brain reward regions Omics studies of the murine intestinal ecosystem exposed to subchronic and mild social defeat stress Microbiota and neurodevelopmental windows: implications for brain disorders Gene-environment interaction and psychiatric disorders: review and future directions, Seminars in cell & developmental biology The human intestinal microbiome in health and disease Gut microbiota depletion from early adolescence in mice: implications for brain and behaviour Gut-brain axis: how the microbiome influences anxiety and depression Impact of stressor exposure on the interplay between commensal microbiota and host inflammation Structural & functional consequences of chronic psychosocial stress on the microbiome & host Twelve-month use of mental health services in the United States: results from the National Comorbidity Survey Replication The hippocampus in depression: more than the sum of its parts? Advanced hippocampal substructure segmentation in depression Psychiatric genome-wide association study analyses implicate neuronal, immune and histone pathways The effects of dietary improvement on symptoms of depression and anxiety: a metaanalysis of randomized controlled trials Effects of nutrients on brain function Depression and physical activity in a sample of Nigerian adolescents: levels, relationships and predictors Associations between physical activity and reduced rates of hopelessness, depression, and suicidal behavior among college students Physical activity and the prevention of depression: a systematic review of prospective studies Exercise therapy improves both mental and physical health in patients with major depression The effects of tyrosine depletion in normal healthy volunteers: implications for unipolar depression Serotonin by stress interaction: a susceptibility factor for the development of depression? Branched-chain amino acids and brain function Diagnostic developments involving cellfree (circulating) nucleic acids RNA in unexpected places: long noncoding RNA functions in diverse cellular contexts Activation of Platelet NLRP3 inflammasome in Crohn's disease The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets Cytoscape: a software environment for integrated models of biomolecular interaction networks KOBAS-i: intelligent prioritization and exploratory visualization of biological functions for gene enrichment analysis An interactive tool for comparing lists with Venn Diagrams LncTarD: a manually-curated database of experimentallysupported functional lncRNA-target regulations in human diseases RCSB Protein Data Bank: powerful new tools for exploring 3D structures of biological macromolecules for basic and applied research and education in fundamental biology, biomedicine, biotechnology, bioengineering and energy sciences PubChem in 2021: new data content and improved web interfaces SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules Investigating protein-peptide interactions using the Schrödinger computational suite, Modeling peptide-protein interactions Small-molecule library screening by docking with PyRx UCSF Chimera-a visualization system for exploratory research and analysis Antidepressant effect of fraxinus rhynchophylla hance extract in a mouse model of chronic stress-induced depression Subchronic and mild social defeat stress accelerates food intake and body weight gain with polydipsia-like features in mice The recent progress in animal models of depression Towards a psychobiotic therapy for depression: Bifidobacterium breve CCFM1025 reverses chronic stressinduced depressive symptoms and gut microbial abnormalities in mice GluD1 knockout mice with a pure C57BL/6N background show impaired fear memory, social interaction, and enhanced depressive-like behavior Oxygen consumption in adult and AGED C57BL/6J mice during acute treadmill exercise of different intensity PPARγ/Pgc-1α-Fndc5 pathway up-regulation in gastrocnemius and heart muscle of exercised, branched chain amino acid diet fed mice Enteric dopaminergic neurons: definition, developmental lineage, and effects of extrinsic denervation Platelets in amyloidogenic mice are activated and invade the brain The RAGE signaling pathway is involved in intestinal inflammation and represents a promising therapeutic target for Inflammatory Bowel Diseases Increased expression of the receptor for advanced glycation end products in neurons and astrocytes in a triple transgenic mouse model of Alzheimer's disease Physical activity ameliorates the function of organs via adipose tissue in metabolic diseases Analysis of oxidative stress-related markers in Crohn's disease patients at surgery and correlations with clinical findings Expression of notch receptors and ligands in the adult gut PTEN in prefrontal cortex is essential in regulating depressionlike behaviors in mice Effects of mindfulness exercises as stand-alone intervention on symptoms of anxiety and depression: systematic review and meta-analysis Moderating effects of exercise duration and intensity in neuromuscular vs. endurance exercise interventions for the treatment of depression: a meta-analytical review Exercise interventions for the prevention of depression: a systematic review of meta-analyses Complementary medicine, exercise, meditation, diet, and lifestyle modification for anxiety disorders: a review of current evidence Exercise training-induced changes in metabolic syndrome parameters, carotid wall thickness, and thyroid function in middle-aged women with subclinical hypothyroidism, Pflügers Archiv-European Dietary intake of branched-chain amino acids in relation to depression, anxiety and psychological distress VEGF treatment ameliorates depression-like behavior in adult offspring after maternal immune activation Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations Acknowledgements We thank our colleagues for their association and helpful discussions in this study.