# Transcriptomic Studies of Antidepressant Action in Rodent Models of Depression: A First Meta-Analysis

El Chérif Ibrahim, Victor Gorgievski, Pilar Ortiz-Teba, Raoul Belzeaux, Gustavo Turecki, Etienne Sibille, Guillaume Charbonnier, Eleni T Tzavara

► To cite this version:

El Chérif Ibrahim, Victor Gorgievski, Pilar Ortiz-Teba, Raoul Belzeaux, Gustavo Turecki, et al.. Transcriptomic Studies of Antidepressant Action in Rodent Models of Depression: A First Meta-Analysis. International Journal of Molecular Sciences, 2022, 23, 10.3390/ijms232113543. [https://doi.org/10.3390/ijms232113543](https://doi.org/10.3390/ijms232113543)

HAL Id: hal-03851878

https://hal.science/hal-03851878

Submitted on 14 Nov 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Transcriptomic Studies of Antidepressant Action in Rodent Models of Depression: A First Meta-Analysis

El Chérif Ibrahim<sup>1,2,*</sup>, Victor Gorgievski<sup>3,4</sup>, Pilar Ortiz-Teba<sup>5 ID</sup>, Raoul Belzeaux<sup>1,2,6</sup>, Gustavo Turecki<sup>7 ID</sup>, Etienne Sibille<sup>8</sup>, Guillaume Charbonnier<sup>9 ID</sup> and Eleni T. Tzavara<sup>2,3,4,10,*</sup>

1 Aix-Marseille Univ, CNRS, INT, Inst Neurosci Timone, 13005 Marseille, France

2 Fondation FondaMental, 94000 Créteil, France

3 Université Paris Cité, CNRS, Integrative Neuroscience and Cognition Center, 75006 Paris, France

4 Sorbonne Université, INSERM, CNRS, Neuroscience Paris Seine, Institut de Biologie Paris Seine, 75005 Paris, France

5 Laboratory of Neuropharmacology, Department of Medicine and Life Sciences, Universitat Pompeu Fabra, 08002 Barcelona, Spain

6 Pôle Universitaire de Psychiatrie, CHU de Montpellier, 34295 Montpellier, France

7 Douglas Mental Health University Institute, McGill University, Montreal, QC H4H 1R3, Canada

8 Department of Psychiatry, University of Toronto, Toronto, ON M5T 1R8, Canada

9 Aix-Marseille Univ, INSERM, TAGC, 13009 Marseille, France

10 Hôpital Sainte Marguerite AP-HM, Pôle de Psychiatrie, 13274 Marseille, France

* Correspondence: el-cherif.ibrahim@univ-amu.fr (E.C.I.); eleni.tzavara@inserm.fr (E.T.T.)

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Ibrahim, E.C.; Gorgievski, V.; Ortiz-Teba, P.; Belzeaux, R.; Turecki, G.; Sibille, E.; Charbonnier, G.; Tzavara, E.T. Transcriptomic Studies of Antidepressant Action in Rodent Models of Depression: A First Meta-Analysis. *Int. J. Mol. Sci.* **2022**, *23*, 13543. https://doi.org/10.3390/ijms232113543

Academic Editor: Larisa
Bobrovskaya

Received: 4 October 2022
Accepted: 31 October 2022
Published: 4 November 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Antidepressants (ADs) are, for now, the best everyday treatment we have for moderate to severe major depressive episodes (MDEs). ADs are among the most prescribed drugs in the Western Hemisphere; however, the trial-and-error prescription strategy and side-effects leave a lot to be desired. More than 60% of patients suffering from major depression fail to respond to the first AD they are prescribed. For those who respond, full response is only observed after several weeks of treatment. In addition, there are no biomarkers that could help with therapeutic decisions; meanwhile, this is already true in cancer and other fields of medicine. For years, many investigators have been working to decipher the underlying mechanisms of AD response. Here, we provide the first systematic review of animal models. We thoroughly searched all the studies involving rodents, profiling transcriptomic alterations consecutive to AD treatment in naïve animals or in animals subjected to stress-induced models of depression. We have been confronted by an important heterogeneity regarding the drugs and the experimental settings. Thus, we perform a meta-analysis of the AD signature of fluoxetine (FLX) in the hippocampus, the most studied target. Among genes and pathways consistently modulated across species, we identify both old players of AD action and novel transcriptional biomarker candidates that warrant further investigation. We discuss the most prominent transcripts (immediate early genes and activity-dependent synaptic plasticity pathways). We also stress the need for systematic studies of AD action in animal models that span across sex, peripheral and central tissues, and pharmacological classes.

**Keywords:** major depression; antidepressant; fluoxetine; hippocampus; transcriptomics

## 1. Introduction

Major depression, which affects nearly 5% of the world's adult population, has become the leading cause of disability, and thus represents a heavy burden for our societies [1–3]. However, the clinical management of major depression remains "artisanal", involving many "trial-and-error" attempts based on the intuition of the medical practitioner [4].

Moderate to severe depressive episodes are treated with antidepressant (AD) drugs including selective serotonin reuptake inhibitors (SSRIs), serotonin and norepinephrine reuptake inhibitors (SNRIs), and older, less selective compounds. More than 30 different
ADs are available today at an ambulatory setting, and there is no clear difference in overall efficacy between the different AD treatments [5]. Despite their known activity toward targets involved in the pathophysiology of major depression [6] (all of these ADs but the newly approved ketamine target the monoaminergic system, acting on a more or less selective combination of monoamine transporters and/or receptors), they do not provide an effective response for many patients. More than 60% of patients suffering from a major depressive episode (MDE) fail to achieve remission following the AD they are prescribed. Even more, one-third of patients will not reach remission after up to four AD trials [7]. For those who respond, full response is only observed after several weeks of treatment. This delay negatively affects patient compliance, increases overall pharmacological burden, and potentially dangerous side effects.

The discovery of novel targets for better and safer Ads, as well as the identification of biomarkers related to therapeutic efficacy, is urgently needed [8]. For the moment, despite efforts in this direction, there is no reliable biomarker that can predict the therapeutic response in a patient suffering from an MDE, and there is no absolute predictor to guide the choice of the therapeutic approach [9]. For this, we believe it is necessary to know exactly what the effects of AD drugs are, in particular the changes in the gene expression program that they induce over time [10], in the parts of our body directly involved in the processing of emotions and mood.

Indeed, RNA, as an immediate product of gene expression and epigenetic programs, is a perfect reflection of the functional status of an individual. For this reason, many investigators have compared peripheral RNA expression profiles between subjects suffering from major depression and control subjects or between patients at different times during their AD treatment [11,12]. Unfortunately, the complexity of the conditions of human subjects in terms of personal history (stress at different times of life, age, sex, obesity, cardiovascular problems, infectious history, tobacco and alcohol consumption, hereditary factors, and dietary habits) and the different undergoing pharmacological treatments at the time of study require very large cohorts to achieve sufficient statistical power and allow the identification of universal signatures. Furthermore, it would be necessary to study the signatures induced by AD drugs in a context of healthy subjects to compare them to those of subjects suffering from an MDE if we want to be able to effectively select new therapeutic targets. However, from this point of view, the data available in humans are very limited.

To better understand and manage affective disorders, major efforts have been made in recent years to characterize the transcriptional signatures of major depression and AD action that distinguish patients from healthy subjects or their equivalent in animal models, particularly in mice and rats [13]. Animal models provide a relatively homogeneous study population in which the events triggering the onset of depressive-type symptoms can be precisely controlled, as well as the timing of the administration of a particular drug at a specific dose and for a defined duration [14]. In addition, animal models allow access to all tissues that may be relevant for both pathophysiology and biomarker development, which is not always the case in humans for evident ethical and technical reasons. Furthermore, it should be noted that in humans, for the time being, most of the transcriptional data have been obtained from post-mortem brain tissue and that it is complicated in these circumstances to disentangle the signature of treatments from that of the long disease process and the terminal effect of death itself. This is something that can be much better controlled in an animal model where, also, new technologies may rapidly advance our understanding of psychiatric disorders at the cellular, molecular, and brain circuit levels [15].

Much progress has been made on the mechanisms leading to the phenomena of vulnerability or, on the contrary, resilience to stress using animal models and the modalities of the response to AD treatment, in particular with regard to transcriptional signatures [16–20]. As in humans, animal models reveal absence of homogeneous response, which offers the possibility of selecting animals that respond well to AD treatment to study the biological pathways involved in the recovery of normal social, emotional, hedonic, affective, mnemonic, exploratory, and feeding behaviors [21]. Going a step beyond, recent conver-
gence studies have been undertaken from human and animal transcriptomes and should advance our knowledge of the pathophysiology of major depression [16,22–25]. Concerning the identification of signatures of AD treatments, several teams have tried to gather common features of genome-wide transcriptional variations induced by the action of ADs in patients suffering from affective disorders [11,12,26], but to our knowledge, studies carried out in rodent models are very scarce and fragmented.

Despite the privileged use of SSRIs in the pharmacological armamentarium to treat major depression, and the widespread use of this pharmacological class in omics studies in humans and animal models [27], we found only one review regarding studies conducted in animals to evaluate the transcriptional effects of these drugs [28]. However, this work is necessary to establish the targets common to humans and animals that would deserve special attention to improve therapeutic efficacy, which remains too limited today.

Since, to our knowledge, no study has established common patterns of gene expression program variation in rodent models of AD drug exposure, we searched all available transcriptomic data in rats and mice with AD drug exposure and sought to define recurrent elements of a biological signature with already known biological processes involved in mood disorder pathophysiology, as well as other processes that should merit further investigations.

## 2. Results

### 2.1. Preponderance of Studies Investigating the Effect of Fluoxetine

We first listed all studies reporting effects of a confirmed or putative antidepressant intervention on the variation of transcriptional expression in the mouse or the rat in (i) naïve animals and (ii) in animals subjected to a paradigm modeling negative affects and depression-like behaviors. Focusing on studies using sufficiently mature technologies and considering the whole genome, we identified 108 published articles and 6 publicly available but unpublished transcriptomic datasets (Figure 1, Table S1). We noticed, however, that in some cases the same animal cohorts have been used in multiple publications. We inventoried 107 independent animal cohorts.

**Sex:** Strikingly, most of the studies were conducted on male animals since out of 107 cohorts only 8 involved females [29–38], and, in fact, RNA profiling on males and females subjected to the same paradigm was performed only four times [32–34,36–38]. In contrast, the proportions are much more balanced between studies in mice and rats or between studies conducted in naïve or stressed animals (Figure 2A–C).

**Nature of AD interventions:** When we look at the class of the AD treatments tested (Figure 2D), we see that almost half of the cohorts used SSRIs, in particular fluoxetine (FLX), which was used in more than a third of the studies. The second most represented class is tricyclics, especially imipramine, followed by the NMDA receptor (NMDAR) antagonists, essentially ketamine. Only after that do we find the SNRIs, followed by mood stabilizers (lithium), atypical antipsychotics, electrical stimulation of the brain (essentially by electroconvulsive therapy, ECT), monoamine oxidase inhibitors (MAOIs), and histone deacetylase (HDAC) inhibitors. Finally, while there is a wide variety of AD interventions (both drugs and non-drugs), of the 65 ADs we have listed (Table S2), more than half (N = 40) have only been tested for pan-genomic activity once.

**Animal models:** Numerous paradigms have been developed in rodents to model symptoms related to human mood disorders and this diversity is reflected in the transcriptomic studies that we surveyed (Figure 2E). Thus, while no fewer than 13 different paradigms have been used (Table S3), the most popular is unpredictable chronic moderate stress (UCMS) used in 40% of studies, while electric shocks, restraint stress, social defeat, and maternal separation each concerned 8–10% of the studies.

**Region of interest:** From an anatomical point of view (Figure 2F), the brain is the most studied organ, overwhelmingly with 105 out of 107 cohorts, while for the periphery, just three studies profiled blood [13,21,39,40], and only one study examined the adrenal gland [41], liver and kidney [42,43], or mammary glands [31]. Concerning the brain, there
is a strong heterogeneity in the size of the cerebral areas or the type of cells profiled ranging from the whole brain to the paraventricular nucleus of the hypothalamus [44] to a particular type of cortical neurons [37,45,46], but, in fact, the most studied area is the hippocampus (as a whole or more precisely the dentate gyri) that we found in almost half of the studies, far ahead of the prefrontal cortex, the frontal cortex, the amygdala and the nucleus accumbens. Finally, although animal models offer the possibility of easily recovering several different anatomical areas and thus determine the possible regionalization of expression profiles, most studies only sampled one anatomical area and only four studies distinguished more than three anatomical areas from the same animals [30,47–49].

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Flow diagram outlining the selection procedure of the studies kept in the meta-analysis.</figcaption></figure>

Given the great variability in terms of AD intervention, paradigms used, species, sex, brain regions, and assay methods, it is difficult to make meaningful comparison across treatments. This is even more difficult considering the fact that a little more than half of the expression data are not publicly available (Figure 1). Therefore, we decided to focus on the effect of FLX on the hippocampus, which nevertheless represents 12 studies (Table 1).

## 2.2. *Signature of FLX Response in Stressed Rodents*

After retrieving the complete data for each of the 12 selected studies, we established for each dataset the lists of differentially expressed genes between the control condition and the FLX treatment condition. First, we grouped studies in mice and rats where FLX treatment restored behavior in stressed animals, which represents four studies in mice and three in rats (Table 1). We felt that it would not be appropriate to merge all the raw data (after different
calibration procedures) into one large dataset because of the great heterogeneity in the specifics of each study. In addition to the differences between species and strains, sample sizes, transcriptomic platforms, and behavioral protocols also differed between studies. Instead, we kept from each individual study the list of the top 300 genes with nominal *p* values below the arbitrary threshold of 0.05. Then, we combined the individual lists into a single graph to allow identification of genes with repeated dysregulated expression in multiple experiments (Figure S1). In this way, we extracted 1973 genes that were modulated by FLX in at least one experiment in stressed animals. Of these genes, 98 were dysregulated in at least four different experiments (Table S4). We also defined for each gene the global trend of the direction of variation of the transcriptional expression with a positive consensus score for overexpression and negative for underexpression. Of note, although nine genes (*Arc*, *Ddah1*, *Egr1*, *Hmgcs1*, *Kcng2*, *Klhl5*, *Nr4a1*, *Oxtr*, and *ZfHX2*) show significant expression variation in at least five of the experiments, only four show a clear directional pattern (*Ddah1* and *Oxtr* are predominantly upregulated, whereas *Klhl5* and *ZfHX2* are predominantly downregulated). Interestingly, at best, only 15 genes show a significant direction of variation that is consistent across most experiments (Table 2).

In the literature there is no consensus way to establish a signature based on a collection of different datasets. As an alternative approach, we considered defining the FLX action portrait based on the methodology recently proposed by Stephen Gammie to establish a signature of major depression in the human brain [23]. Applying this methodology, we ranked 39,629 genes and obtained scores ranging from +4.56 for the most overexpressed consensually in the seven studies (immediate early response gene, *Ier5*) to -4.46 for the most underexpressed (*Sh3d19*, Table S5). We listed in Table 3 the signature of effect of FLX with genes showing absolute values of portrait scores equal or above 4. When we compared the integration (Table 2, N = 22) with the portrait score (Table 3, N = 12) signature, we identified two genes (*Hmgcs1* and *Prkar1b*) that are common to both signatures. The advantage of having a sign indicating the overall direction of expression variation, whether with the consensus or the portrait score, allows us to assess whether there is any distortion in the directionality of FLX-induced transcriptional expression variation. We can see that this is not the case because with both scoring methods a little less than 15% of the genes show a directionality of their expression variation in one direction or another and thus for nearly 70% of the genes no trend is apparent. To extract the most salient variations with the portrait score, we separated the genes with scores whose absolute value is equal or greater than 3 and thus obtained 82 overexpressed genes and 80 underexpressed genes (Table S5). As a more general and consensual signature of FLX in a stress paradigm, we combined both scoring procedures with absolute values equal or above 2, resulting in 412 upregulated genes and 411 downregulated genes in an almost perfect balance (Table S5).

To identify the biological processes targeted by the signatures described above, we conducted ontological analysis with the DAVID algorithm. Table S6 shows ontological results for the list of the 98 most commonly deregulated genes from Table S4 (integration method). Two overrepresented processes emerged: the activity of transcription factors and the synapse. Then, we performed ontological analysis on the 823 genes with consensus and portrait scores above 2 after separating up- (Table 4) and downregulated genes (Table 5). Overexpressed genes were strongly associated with the synapse, postsynaptic density, and cell junction, as well as with glutamatergic and oxytocin signaling pathways. We also found recurrent modules classically associated with AD function. On the contrary, underexpressed genes were largely enriched for the translation machinery, the ribosome.

## 2.3. *Signature of FLX Response in Naïve Rodents*

Regarding the signature of the effect of FLX on naïve animals, we surveyed seven different datasets, five in mice and two in rats; 1763 genes showed a variation in expression in at least one experiment (Figure S2) and 117 in at least four experiments (Table S7). Sixteen genes were significantly affected by FLX in five out of seven experiments, and in fact, 36 genes showed significant variation in the same direction in most experiments (Table 6).
The most significantly upregulated gene is *Myo1e* (five experiments), whereas the most significantly downregulated genes are *Isoc1*, *Map1a*, *Scn3b*, and *Zfhx2* (five experiments). The application of the portrait scoring method allowed us to rank 36,483 genes with a score ranging from +5.56, for the most overexpressed consensually in the seven studies (vasoactive intestinal peptide, *Vip*), to -5.58 for the most underexpressed (intersectin 1, *Itsn1*, Table S8). We listed in Table 7 the portrait signature of effect of FLX with 60 genes showing absolute values of portrait scores above 4. Eight upregulated genes (*Cfh*, *Ddr1*, *Gsn*, *Homer1*, *Igfbp6*, *Knstrn*, *Sel1l3*, and *Sema3a*) and nine downregulated genes (*Doc2b*, *Fat4*, *Itga4*, *Itsn1*, *Pcdh19*, *Rasgrf1*, *Scn3b*, *Tnxb*, and *Zfp316*) are common to the signatures we obtained with both methods (Tables 6 and 7).

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Main features of transcriptomic studies conducted in rodents to investigate the effect of AD treatments. Distribution of the sex (A), rodent species (B), and use of a stimulation to induce depression-like behaviors before testing the effect of ADs (C). (D) Diversity of the class of ADs used with detailed proportions for SSRIs and tricyclics. (E) Variety of paradigms used to induce depression-like symptoms. (F) Anatomical distribution of collected tissues. Numbers in parentheses indicate the number of studies where an AD (D), a type of stress (E), or a specific piece of tissue (F) has been examined.</figcaption></figure>
<table><caption>**Table 1.** Studies selected for meta-analysis.</caption><thead><tr><th>Meta-Analysis</th><th>ID</th><th>Ref.</th><th>Rodent</th><th>Stress</th><th>Fluoxetine</th><th>Hippocampus Dissection</th><th>Sample Number</th><th>Platform</th></tr></thead><tbody><tr><td rowspan="4">mouse S-FLX vs. S-Veh</td><td>GSE84185</td><td>[21]</td><td>8-week BALB/c</td><td>8-week UCMS</td><td>In drinking water (10-20 mg/kg/day) for the last 6 weeks</td><td>dentate gyrus</td><td>N = 21: 8 NS-Veh; 8 S-Veh; 5 S-FLX-R</td><td>Agilent SurePrint G3 Mouse Gene Expression 8 × 60 K</td></tr><tr><td>private</td><td>[50]</td><td>8-week BALB/c</td><td>7-week UCMS</td><td>ip administration (20 mg/kg/day) from week 2 to the end of UCMS</td><td>dentate gyrus 5 h after the last injection</td><td>N = 18: 6 NS-Sal; 6 S-Sal; 6 S-FLX</td><td>Affymetrix GeneChip Mouse Genome 430 2.0</td></tr><tr><td>SRP057486</td><td>[48]</td><td>7-8-week C57BL/6NCrl</td><td>fear conditioning by electric foot shock</td><td>in drinking water (20 mg/kg/day) 12 h after conditioning and for 28 days, followed by 28-day wash out</td><td>dorsal cornu ammonis 1 (CA1) 9 weeks after stress</td><td>N = 9: 3 NS-Veh; 3 S-Veh; 3 S-FLX</td><td>Illumina Genome Analyzer II system</td></tr><tr><td>GSE43261</td><td>[51]</td><td>7-8-week C57BL/6J</td><td>corticosterone (35 µg/mL) in drinking water for 21 days</td><td>in drinking water (160 µg/mL) for 21 days</td><td>dorsal or ventral dentate gyrus</td><td>N = 30: (8 NS-Veh; 7 S-FLX-R) per area</td><td>Affymetrix Mouse Genome 430 2.0</td></tr><tr><td rowspan="3">rat S-FLX vs. S-Veh</td><td>GSE56028</td><td>[52]</td><td>2-month Wistar</td><td>6-week UCMS</td><td>ip administration (10 mg/kg/day in ultra-pure water), for the last 2 weeks</td><td>dentate gyrus</td><td>N = 9: 3 NS-Veh; 3 S-Veh; 3 S-FLX</td><td>Affymetrix Rat Gene 1.0 ST</td></tr><tr><td>SRP131063</td><td>NA</td><td>6-7-week Sprague-Dawley</td><td>5-week UCMS</td><td>ip administration (7 mg/kg/day) the last week of UCMS.</td><td>hippocampus</td><td>N = 9: 3 NS-Sal; 3 S-Sal; 3 S-FLX</td><td>Illumina HiSeq 4000</td></tr><tr><td>SRP084288</td><td>[53]</td><td>6-week Sprague-Dawley</td><td>restraint stress (6 h/day) for 28 days</td><td>administration by gavage (10 mg/kg/day) one hour before stress for 28 days</td><td>hippocampus</td><td>N = 3: 1 CTL; 1 S; 1 S-FLX</td><td>Illumina Nextseq 500/151 nt</td></tr><tr><td rowspan="4">mouse FLX vs. Veh</td><td>GSE84185</td><td>[21]</td><td>8-week BALB/c</td><td>No</td><td>In drinking water (10-20 mg/kg/day) for 6 weeks</td><td>dentate gyrus</td><td>N= 16: 8 NS-Veh; 8 NS-FLX</td><td>Agilent SurePrint G3 Mouse Gene Expression 8 × 60 K</td></tr><tr><td>SRP057486</td><td>[48]</td><td>7-8-week C57BL/6NCrl</td><td>No</td><td>in drinking water (20 mg/kg/day) for 28 days, followed by 5-week wash out</td><td>dorsal cornu ammonis 1 (CA1)</td><td>N = 6: 3 Veh; 3 FLX</td><td>Illumina Genome Analyzer II system</td></tr><tr><td>GSE118669</td><td>[54]</td><td>9-week C57BL/6J</td><td>No</td><td>subcutaneous releasing pellet (15 mg/kg/day) in the dorsal interscapular region for 3 weeks</td><td>dentate gyrus</td><td>N = 16: 8 Veh; 8 FLX</td><td>Affymetrix Mouse Genome 430 2.0</td></tr><tr><td>GSE54307</td><td>[55]</td><td>8-week C57BL/6J</td><td>No</td><td>in drinking water (22 mg/kg/day) for 4 weeks</td><td>dentate gyrus 24 h after the end of FLX treatment</td><td>N = 2: 1 Veh; 1 FLX</td><td>Affymetrix Mouse Genome 430 2.0</td></tr><tr><td rowspan="2">rat FLX vs. Veh</td><td>GSE6476</td><td>[56]</td><td>3-5-week DBA/2J</td><td>No</td><td>in drinking water (18 mg/kg/day) for 3 weeks</td><td>hippocampus</td><td>N = 4: 2 CTL; 2 FLX</td><td>Affymetrix Mouse Genome 430 2.0</td></tr><tr><td>SRP056480</td><td>[57]</td><td>Wistar</td><td>No</td><td>oral gavage (12 mg/kg/day) or Veh (1% methylcellulose) from PND 67 to 88</td><td>hippocampus at PND 128</td><td>N = 4: 2 pools of 2 rats per group (Veh; FLX)</td><td>Illumina Genome Analyzer IIx</td></tr><tr><td></td><td>GSE42940</td><td>[58]</td><td>2-month Sprague-Dawley</td><td>No</td><td>po (10 mg/kg) or Veh (5% sucrose) from PND 2 to 21 and wash out for 5 weeks</td><td>hippocampus</td><td>N = 8: 4 Veh; 4 FLX</td><td>Agilent TIFR-Rat-8 × 15 K v1.0</td></tr></tbody></table>

CTL, control; FLX, fluoxetine; ip, intraperitoneal; NS, non-stressed; po, per os; S, stressed; Sal, saline; UCMS, unpredictable chronic mild stress, Veh, vehicle.
<table><caption>**Table 2.** Signature of FLX effect in stressed animals (integration method).</caption><thead><tr><th rowspan="2">Dataset ID<br>Gene</th><th colspan="2">GSE43261</th><th colspan="2">GSE56028</th><th colspan="2">GSE84185</th><th colspan="2">SRP057486</th><th colspan="2">SRP084288</th><th colspan="2">SRP131063</th><th colspan="2">Surget 2009</th><th rowspan="2">Consensus<br>Score</th></tr><tr><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th></tr></thead><tbody><tr><td>Adprm</td><td>-0.096</td><td>4.25 × 10<sup>-2</sup></td><td>-0.251</td><td>5.43 × 10<sup>-3</sup></td><td>-0.105</td><td>4.01 × 10<sup>-2</sup></td><td>0.035</td><td>6.92 × 10<sup>-1</sup></td><td>-0.287</td><td>7.68 × 10<sup>-3</sup></td><td>0.054</td><td>6.83 × 10<sup>-1</sup></td><td>-0.012</td><td>9.11 × 10<sup>-1</sup></td><td>-4</td></tr><tr><td>Arc</td><td>1.796</td><td>9.70 × 10<sup>-3</sup></td><td>-0.379</td><td>3.20 × 10<sup>-2</sup></td><td>1.946</td><td>5.65 × 10<sup>-6</sup></td><td>-0.011</td><td>7.98 × 10<sup>-1</sup></td><td>-0.001</td><td>9.61 × 10<sup>-1</sup></td><td>0.282</td><td>3.87 × 10<sup>-4</sup></td><td>-0.390</td><td>2.18 × 10<sup>-2</sup></td><td>1</td></tr><tr><td>Ddah1</td><td>0.370</td><td>1.03 × 10<sup>-4</sup></td><td>-0.051</td><td>6.22 × 10<sup>-1</sup></td><td>0.180</td><td>1.58 × 10<sup>-2</sup></td><td>0.104</td><td>2.28 × 10<sup>-1</sup></td><td>0.192</td><td>1.96 × 10<sup>-2</sup></td><td>-0.205</td><td>1.75 × 10<sup>-3</sup></td><td>0.174</td><td>2.51 × 10<sup>-2</sup></td><td>3</td></tr><tr><td>Egr1</td><td>1.434</td><td>5.36 × 10<sup>-3</sup></td><td>-0.091</td><td>4.48 × 10<sup>-1</sup></td><td>1.883</td><td>1.13 × 10<sup>-5</sup></td><td>-0.308</td><td>5.97 × 10<sup>-4</sup></td><td>-0.117</td><td>4.20 × 10<sup>-1</sup></td><td>0.255</td><td>1.48 × 10<sup>-3</sup></td><td>-0.286</td><td>2.65 × 10<sup>-2</sup></td><td>1</td></tr><tr><td>Ephb6</td><td>1.084</td><td>8.54 × 10<sup>-12</sup></td><td>-0.176</td><td>7.84 × 10<sup>-2</sup></td><td>1.328</td><td>3.22 × 10<sup>-8</sup></td><td>-0.034</td><td>7.02 × 10<sup>-1</sup></td><td>0.242</td><td>4.84 × 10<sup>-3</sup></td><td>0.269</td><td>1.71 × 10<sup>-4</sup></td><td>-0.037</td><td>7.49 × 10<sup>-1</sup></td><td>4</td></tr><tr><td>Hmgs1</td><td>0.766</td><td>7.36 × 10<sup>-5</sup></td><td>-0.218</td><td>3.04 × 10<sup>-2</sup></td><td>0.336</td><td>4.46 × 10<sup>-2</sup></td><td>-0.145</td><td>7.06 × 10<sup>-2</sup></td><td>0.089</td><td>1.75 × 10<sup>-1</sup></td><td>-0.410</td><td>3.43 × 10<sup>-11</sup></td><td>-0.217</td><td>1.73 × 10<sup>-2</sup></td><td>-1</td></tr><tr><td>Krtg2</td><td>-0.750</td><td>2.29 × 10<sup>-14</sup></td><td>0.146</td><td>4.95 × 10<sup>-2</sup></td><td>-1.648</td><td>1.07 × 10<sup>-5</sup></td><td>0.094</td><td>2.29 × 10<sup>-1</sup></td><td>0.238</td><td>2.84 × 10<sup>-2</sup></td><td>0.127</td><td>4.95 × 10<sup>-2</sup></td><td>0.075</td><td>5.17 × 10<sup>-1</sup></td><td>1</td></tr><tr><td>Klhl5</td><td>-0.394</td><td>2.91 × 10<sup>-4</sup></td><td>-0.073</td><td>3.34 × 10<sup>-1</sup></td><td>-0.562</td><td>1.95 × 10<sup>-5</sup></td><td>-0.013</td><td>8.71 × 10<sup>-1</sup></td><td>0.242</td><td>4.04 × 10<sup>-2</sup></td><td>-0.172</td><td>8.79 × 10<sup>-3</sup></td><td>-0.177</td><td>2.34 × 10<sup>-2</sup></td><td>-3</td></tr><tr><td>Lzts1</td><td>0.301</td><td>1.35 × 10<sup>-3</sup></td><td>0.014</td><td>8.95 × 10<sup>-1</sup></td><td>0.723</td><td>6.19 × 10<sup>-4</sup></td><td>-0.020</td><td>8.23 × 10<sup>-1</sup></td><td>0.345</td><td>7.62 × 10<sup>-5</sup></td><td>0.165</td><td>3.34 × 10<sup>-2</sup></td><td>-0.096</td><td>5.85 × 10<sup>-1</sup></td><td>4</td></tr><tr><td>Mef2d</td><td>0.202</td><td>5.24 × 10<sup>-4</sup></td><td>0.038</td><td>7.29 × 10<sup>-1</sup></td><td>0.646</td><td>1.82 × 10<sup>-7</sup></td><td>0.042</td><td>6.25 × 10<sup>-1</sup></td><td>0.227</td><td>2.63 × 10<sup>-2</sup></td><td>0.263</td><td>1.14 × 10<sup>-5</sup></td><td>-0.101</td><td>1.49 × 10<sup>-1</sup></td><td>4</td></tr><tr><td>Nfib</td><td>-0.224</td><td>4.91 × 10<sup>-3</sup></td><td>0.139</td><td>6.01 × 10<sup>-2</sup></td><td>-0.929</td><td>4.53 × 10<sup>-7</sup></td><td>-0.052</td><td>5.40 × 10<sup>-1</sup></td><td>-0.329</td><td>4.12 × 10<sup>-4</sup></td><td>-0.062</td><td>4.02 × 10<sup>-1</sup></td><td>-0.310</td><td>1.20 × 10<sup>-2</sup></td><td>-4</td></tr><tr><td>Nr4a1</td><td>1.390</td><td>2.01 × 10<sup>-2</sup></td><td>-0.395</td><td>4.02 × 10<sup>-2</sup></td><td>1.906</td><td>1.79 × 10<sup>-4</sup></td><td>-0.201</td><td>2.30 × 10<sup>-2</sup></td><td>-0.072</td><td>9.00 × 10<sup>-1</sup></td><td>0.352</td><td>1.08 × 10<sup>-5</sup></td><td>-0.303</td><td>2.63 × 10<sup>-2</sup></td><td>0</td></tr><tr><td>Oxtr</td><td>1.420</td><td>1.52 × 10<sup>-5</sup></td><td>0.214</td><td>4.93 × 10<sup>-2</sup></td><td>-1.121</td><td>5.44 × 10<sup>-4</sup></td><td>0.138</td><td>4.23 × 10<sup>-2</sup></td><td>-0.012</td><td>8.13 × 10<sup>-1</sup></td><td>0.091</td><td>4.18 × 10<sup>-2</sup></td><td>-0.051</td><td>7.59 × 10<sup>-1</sup></td><td>3</td></tr><tr><td>Ppara</td><td>-0.165</td><td>8.29 × 10<sup>-3</sup></td><td>-0.187</td><td>3.29 × 10<sup>-2</sup></td><td>-0.160</td><td>3.02 × 10<sup>-2</sup></td><td>-0.024</td><td>7.45 × 10<sup>-1</sup></td><td>-0.045</td><td>5.08 × 10<sup>-1</sup></td><td>-0.101</td><td>4.22 × 10<sup>-1</sup></td><td>-0.175</td><td>2.86 × 10<sup>-2</sup></td><td>-4</td></tr><tr><td>Prkar1b</td><td>0.850</td><td>6.71 × 10<sup>-12</sup></td><td>-0.018</td><td>8.51 × 10<sup>-1</sup></td><td>0.774</td><td>3.31 × 10<sup>-7</sup></td><td>-0.005</td><td>9.52 × 10<sup>-1</sup></td><td>0.284</td><td>3.64 × 10<sup>-4</sup></td><td>0.181</td><td>9.95 × 10<sup>-3</sup></td><td>-0.148</td><td>1.14 × 10<sup>-1</sup></td><td>4</td></tr><tr><td>Rimbp2</td><td>0.622</td><td>1.40 × 10<sup>-6</sup></td><td>-0.121</td><td>2.24 × 10<sup>-1</sup></td><td>0.768</td><td>9.56 × 10<sup>-3</sup></td><td>0.179</td><td>3.87 × 10<sup>-2</sup></td><td>0.048</td><td>4.99 × 10<sup>-1</sup></td><td>0.158</td><td>7.54 × 10<sup>-3</sup></td><td>-0.001</td><td>9.95 × 10<sup>-1</sup></td><td>4</td></tr><tr><td>Rps10</td><td>-0.179</td><td>7.62 × 10<sup>-5</sup></td><td>0.000</td><td>1.00</td><td>-0.196</td><td>4.56 × 10<sup>-3</sup></td><td>-0.008</td><td>9.29 × 10<sup>-1</sup></td><td>-0.227</td><td>3.07 × 10<sup>-3</sup></td><td>-0.253</td><td>3.39 × 10<sup>-4</sup></td><td>-0.194</td><td>8.19 × 10<sup>-2</sup></td><td>-4</td></tr><tr><td>Sgsm2</td><td>0.230</td><td>4.56 × 10<sup>-4</sup></td><td>-0.077</td><td>3.07 × 10<sup>-1</sup></td><td>0.180</td><td>2.16 × 10<sup>-2</sup></td><td>0.201</td><td>2.30 × 10<sup>-2</sup></td><td>0.164</td><td>6.88 × 10<sup>-2</sup></td><td>0.123</td><td>4.94 × 10<sup>-2</sup></td><td>-0.048</td><td>4.99 × 10<sup>-1</sup></td><td>4</td></tr><tr><td>Ttkb1</td><td>0.209</td><td>2.86 × 10<sup>-4</sup></td><td>-0.093</td><td>3.12 × 10<sup>-1</sup></td><td>0.386</td><td>1.58 × 10<sup>-6</sup></td><td>-0.019</td><td>8.15 × 10<sup>-1</sup></td><td>0.233</td><td>9.14 × 10<sup>-3</sup></td><td>0.272</td><td>1.71 × 10<sup>-4</sup></td><td>-0.151</td><td>1.97 × 10<sup>-1</sup></td><td>4</td></tr><tr><td>Tlyh3</td><td>0.275</td><td>1.01 × 10<sup>-4</sup></td><td>-0.065</td><td>3.24 × 10<sup>-1</sup></td><td>0.400</td><td>1.68 × 10<sup>-5</sup></td><td>0.035</td><td>5.02 × 10<sup>-2</sup></td><td>0.189</td><td>3.61 × 10<sup>-3</sup></td><td>0.231</td><td>1.14 × 10<sup>-3</sup></td><td>0.061</td><td>4.37 × 10<sup>-1</sup></td><td>4</td></tr><tr><td>Wak1</td><td>-0.158</td><td>1.17 × 10<sup>-3</sup></td><td>0.159</td><td>2.04 × 10<sup>-1</sup></td><td>-0.300</td><td>2.07 × 10<sup>-2</sup></td><td>-0.183</td><td>4.11 × 10<sup>-2</sup></td><td>-0.223</td><td>2.17 × 10<sup>-3</sup></td><td>-0.032</td><td>5.92 × 10<sup>-1</sup></td><td>-0.089</td><td>2.90 × 10<sup>-1</sup></td><td>-4</td></tr><tr><td>Wtbp</td><td>-0.137</td><td>3.81 × 10<sup>-3</sup></td><td>0.004</td><td>9.58 × 10<sup>-1</sup></td><td>-0.271</td><td>4.47 × 10<sup>-3</sup></td><td>0.105</td><td>2.29 × 10<sup>-1</sup></td><td>-0.026</td><td>8.53 × 10<sup>-1</sup></td><td>-0.126</td><td>4.91 × 10<sup>-2</sup></td><td>-0.254</td><td>2.83 × 10<sup>-2</sup></td><td>-4</td></tr><tr><td>Zfxt2</td><td>-0.350</td><td>2.05 × 10<sup>-7</sup></td><td>-0.165</td><td>1.51 × 10<sup>-1</sup></td><td>-0.421</td><td>4.88 × 10<sup>-3</sup></td><td>-0.207</td><td>1.69 × 10<sup>-2</sup></td><td>-0.200</td><td>1.95 × 10<sup>-2</sup></td><td>0.347</td><td>7.94 × 10<sup>-7</sup></td><td>-0.270</td><td>3.20 × 10<sup>-3</sup></td><td>-4</td></tr></tbody></table>

$|log_2FC| > 0.263$ and $p$-values < 0.05 are in bold. Consensus score indicates whether FLX induces mainly upregulation (positive score) or downregulation (negative score) of gene expression.

<table><caption>**Table 3.** Signature of FLX effect in stressed animals (portrait score method).</caption><thead><tr><th rowspan="2">Dataset ID<br>Gene</th><th colspan="2">GSE43261</th><th colspan="2">GSE56028</th><th colspan="2">GSE84185</th><th colspan="2">SRP057486</th><th colspan="2">SRP084288</th><th colspan="2">SRP131063</th><th colspan="2">Surget 2009</th><th rowspan="2">Consensus<br>Score</th><th rowspan="2">Portrait Score</th></tr><tr><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th></tr></thead><tbody><tr><td>Ier5</td><td>1.143</td><td>1.69 × 10<sup>-11</sup></td><td>-0.009</td><td>9.29 × 10<sup>-1</sup></td><td>1.786</td><td>1.39 × 10<sup>-5</sup></td><td>0.107</td><td>2.33 × 10<sup>-1</sup></td><td>0.119</td><td>3.36 × 10<sup>-1</sup></td><td>0.157</td><td>2.00 × 10<sup>-2</sup></td><td>0.238</td><td>8.61 × 10<sup>-2</sup></td><td>3</td><td>4.556666</td></tr><tr><td>Mapk4</td><td>0.928</td><td>9.14 × 10<sup>-7</sup></td><td>-0.040</td><td>6.26 × 10<sup>-1</sup></td><td>1.021</td><td>3.20 × 10<sup>-5</sup></td><td>-0.015</td><td>8.55 × 10<sup>-1</sup></td><td>0.169</td><td>7.03 × 10<sup>-2</sup></td><td>0.176</td><td>1.66 × 10<sup>-3</sup></td><td>0.150</td><td>2.30 × 10<sup>-1</sup></td><td>3</td><td>4.555544</td></tr><tr><td>Dinas2</td><td>0.310</td><td>6.90 × 10<sup>-7</sup></td><td>-0.046</td><td>5.63 × 10<sup>-1</sup></td><td>0.387</td><td>4.61 × 10<sup>-7</sup></td><td>0.110</td><td>2.15 × 10<sup>-1</sup></td><td>-0.028</td><td>4.61 × 10<sup>-1</sup></td><td>0.212</td><td>8.37 × 10<sup>-4</sup></td><td>0.121</td><td>3.74 × 10<sup>-1</sup></td><td>3</td><td>4.454433</td></tr><tr><td>Kcnh3</td><td>0.162</td><td>3.22 × 10<sup>-2</sup></td><td>-0.001</td><td>9.90 × 10<sup>-1</sup></td><td>1.004</td><td>2.57 × 10<sup>-7</sup></td><td>0.169</td><td>5.02 × 10<sup>-2</sup></td><td>0.189</td><td>3.49 × 10<sup>-2</sup></td><td>0.284</td><td>1.83 × 10<sup>-4</sup></td><td>0.048</td><td>5.57 × 10<sup>-1</sup></td><td>4</td><td>4.445666</td></tr><tr><td>Sema7a</td><td>0.786</td><td>7.29 × 10<sup>-6</sup></td><td>0.056</td><td>5.72 × 10<sup>-1</sup></td><td>1.083</td><td>1.89 × 10<sup>-6</sup></td><td>0.037</td><td>6.38 × 10<sup>-1</sup></td><td>0.008</td><td>8.35 × 10<sup>-1</sup></td><td>0.168</td><td>1.43 × 10<sup>-2</sup></td><td>0.158</td><td>1.86 × 10<sup>-1</sup></td><td>3</td><td>4.444466</td></tr><tr><td>Mapk9</td><td>0.350</td><td>1.45 × 10<sup>-6</sup></td><td>0.026</td><td>6.68 × 10<sup>-1</sup></td><td>0.248</td><td>1.13 × 10<sup>-4</sup></td><td>-0.028</td><td>7.22 × 10<sup>-1</sup></td><td>0.218</td><td>5.68 × 10<sup>-3</sup></td><td>-0.031</td><td>5.78 × 10<sup>-1</sup></td><td>0.132</td><td>1.41 × 10<sup>-1</sup></td><td>3</td><td>4.444443</td></tr><tr><td>Nptx2</td><td>2.410</td><td>2.74 × 10<sup>-7</sup></td><td>-0.132</td><td>5.64 × 10<sup>-1</sup></td><td>1.804</td><td>3.80 × 10<sup>-5</sup></td><td>-0.093</td><td>2.98 × 10<sup>-1</sup></td><td>0.209</td><td>3.74 × 10<sup>-2</sup></td><td>0.237</td><td>8.39 × 10<sup>-4</sup></td><td>0.069</td><td>5.27 × 10<sup>-1</sup></td><td>4</td><td>4.333433</td></tr><tr><td>Prkar1b</td><td>0.850</td><td>6.71 × 10<sup>-12</sup></td><td>-0.018</td><td>8.51 × 10<sup>-1</sup></td><td>0.774</td><td>3.31 × 10<sup>-7</sup></td><td>-0.005</td><td>9.52 × 10<sup>-1</sup></td><td>0.284</td><td>3.64 × 10<sup>-4</sup></td><td>0.181</td><td>9.95 × 10<sup>-3</sup></td><td>-0.148</td><td>1.14 × 10<sup>-1</sup></td><td>4</td><td>4.333333</td></tr><tr><td>Hmgs1</td><td>0.766</td><td>7.36 × 10<sup>-5</sup></td><td>-0.218</td><td>3.04 × 10<sup>-2</sup></td><td>0.336</td><td>4.46 × 10<sup>-2</sup></td><td>-0.145</td><td>7.06 × 10<sup>-2</sup></td><td>0.089</td><td>1.75 × 10<sup>-1</sup></td><td>-0.410</td><td>3.43 × 10<sup>-11</sup></td><td>-0.217</td><td>1.73 × 10<sup>-2</sup></td><td>-1</td><td>-4.321111</td></tr><tr><td>Rspo3</td><td>-0.630</td><td>1.34 × 10<sup>-6</sup></td><td>-0.163</td><td>8.89 × 10<sup>-2</sup></td><td>-0.966</td><td>3.40 × 10<sup>-6</sup></td><td>0.063</td><td>4.81 × 10<sup>-1</sup></td><td>-0.269</td><td>2.10 × 10<sup>-2</sup></td><td>0.072</td><td>4.25 × 10<sup>-1</sup></td><td>-0.099</td><td>4.70 × 10<sup>-1</sup></td><td>-3</td><td>-4.443233</td></tr><tr><td>Cnn3</td><td>0.155</td><td>2.09 × 10<sup>-2</sup></td><td>-0.029</td><td>7.31 × 10<sup>-1</sup></td><td>-0.772</td><td>1.08 × 10<sup>-6</sup></td><td>-0.055</td><td>5.13 × 10<sup>-1</sup></td><td>-0.228</td><td>9.77 × 10<sup>-3</sup></td><td>-0.175</td><td>3.09 × 10<sup>-3</sup></td><td>-0.183</td><td>8.63 × 10<sup>-2</sup></td><td>-2</td><td>-4.443445</td></tr><tr><td>Sh3d19</td><td>-0.677</td><td>1.80 × 10<sup>-11</sup></td><td>-0.020</td><td>8.19 × 10<sup>-1</sup></td><td>-0.643</td><td>2.80 × 10<sup>-6</sup></td><td>-0.112</td><td>1.86 × 10<sup>-1</sup></td><td>-0.010</td><td>9.62 × 10<sup>-1</sup></td><td>-0.131</td><td>1.16 × 10<sup>-1</sup></td><td>-0.313</td><td>5.57 × 10<sup>-2</sup></td><td>-2</td><td>-4.455555</td></tr></tbody></table>

$|log_2FC| > 0.263$ and $p$-values < 0.05 are in bold. Genes in bold are also in Table 2. Consensus score indicates whether FLX induces mainly upregulation (positive score) or downregulation (negative score) of gene expression.
The ontological analysis of the processes enriched from the list of 117 genes (Table S7) reveals again a synaptic location with a glutamatergic signaling (Table S9). When performing ontological analysis on the 528 upregulated genes with consensus and portrait scores above or equal to 2, no real specific process was significantly enriched (Table S10). For the 613 downregulated genes, cell junction and synapse were the most significant locations and regulation of NMDA receptor activity, insulin secretion, and circadian entrainment were also affected (Table 8).

## 2.4. Shared Signature of FLX in Stressed and Naïve Rodents

To establish the core of FLX-induced transcriptional activity in the hippocampus, regardless of the paradigm employed, we first examined the convergence between the stress and non-stress signatures previously described and obtained by either integration or scoring methods. Figure 3 shows that only one gene is common to all comparisons, *Zfhx2*. Next, we combined in the same dataset, the 14 comparisons made previously (7 in a stress context and 7 in a naïve context, Figure S3). With the integration method, among 3360 genes modulated in at least one experiment following FLX treatment, 78 were altered in at least 7 experiments, 15 in at least 8 experiments, while *Pdlim5* and *Zfhx2* were altered in 9 and 10 comparisons, respectively (Table S11). Second, by applying the portrait method to the 14 comparisons, we could rank 40,113 genes with top scorer genes, *Sel113* (+7.88) and *Nfia* (−8.13) (Table S12). When we examine the best scoring genes, with absolute values of both integration and portrait scores equal or above 6, we obtain nine upregulated genes: *Ddr1*, *Ier5*, *Igfbp6*, *Nptx2*, *Prkar1b*, *Ptpn5*, *Sel113*, *Tyro3*, and *Zfp703*, and seven downregulated genes: *Akt3*, *Fat4*, *Nfia*, *Pcdh19*, *Rab27a*, *Scn3b*, and *Zfhx2*.

By evaluating lists of equal size obtained by integration method (Table S11) versus consensus and portrait score (absolute values equal or above 5, Table S12), we defined 15 FLX-modulated genes: *Baalc*, *Igfbp6*, *Itga4*, *Nptx2*, *Prkar1b*, *Rasgrf1*, *S100a6*, *Sel113*, *Slc4a4*, *Sorcs1*, *Tmem47*, *Trpm3*, *Zfhx2*, and *Zfp316*. Among these genes, five are also listed in Figure 3 (*Baalc*, *Igfbp6*, *Prkar1b*, *Sel113*, and *Zfhx2*).

Without surprise, the ontological analysis of the processes enriched from the list of 78 genes most affected by FLX action (Table S11 integration method) reveals mainly a synaptic activity (Table S13). When performing ontological analysis on the 113 upregulated genes with consensus and portrait scores above or equal to 4, the axon and the MAPK signaling pathways were significantly enriched (Table 9). For the 151 downregulated genes, the MAPK signaling pathway was only process significantly affected (Table 10). A secondary analysis of the genes involved in the MAPK signaling pathway identified distinct components for up and downregulated genes. Upregulated genes were associated with neurotrophic pathways while downregulated genes were restricted to protein phosphorylation.

The ontological analysis of the processes enriched from the list of 78 genes deregulated after FLX action confirms with both ontological tools that neuronal structure, development, and signaling are the main processes concerned (Tables 8 and 9).
<table><caption>**Table 4.** Ontological analysis of upregulated genes (consensus and portrait scores ≥ 2) after FLX treatment in stressed animals.</caption><thead><tr><th>Category</th><th>ID</th><th>Name</th><th>Genes</th><th>Count</th><th>Fold enrichment</th><th>Padj</th></tr></thead><tbody><tr><td rowspan="4">Cellular component</td><td>GO:0045202</td><td>synapse</td><td>Add1, Adgrb1, Adgrl1, Arhgap44, Arhgdia, Atcay, Baalc, Btbd8, Cacna1g, Cacnb3, Cacng2, Cdkl5, Clstn1, Cyp46a1, Dlg1, Dlg2, Dlgap3, Dnajc6, Drd1, Egr3, Gabbr1, Gabra1, Gabrg2, Git1, Grin1, Grk2, Iqsec3, Kcnb1, Kcnk1, Lrrc4b, Magi2, Mink1, Mthfr, Ncs1, Nlgn2, Ntrk2, Pak5, Palm, Ppm1h, Prkar1b, Psd3, Rims3, Sh2d5, Shank1, Shank3, Slc30a3, Slc4a10, Slc5a7, Slc6a17, Strn4, Stx1a, Sv2a, Syn1, Syndig1, Syngap1, Syngr1, Vgf, Znrf1</td><td>58</td><td>3.61</td><td>2.05 × 10<sup>-14</sup></td></tr><tr><td>GO:0014069</td><td>postsynaptic density</td><td>Add1, Adgrb1, Adgrl1, Arhgap44, Arhgef2, Atp1a1, Baalc, Baiap2, Cacng2, Cdk5r1, Clstn1, Dclk1, Dlg1, Dlg2, Dlgap3, Dnajc6, Git1, Grin1, Iqsec3, Lrp8, Magi2, Mink1, Ncs1, Ntrk2, Palm, Psd3, Sh2d5, Shank1, Shank3, Shisa8, Stx1a, Syn1, Syndig1, Syngap1</td><td>34</td><td>5.15</td><td>1.31 × 10<sup>-11</sup></td></tr><tr><td>GO:0030054</td><td>cell junction</td><td>Adgrb1, Adgrl1, Arhgap44, Arhgef2, Atcay, Baalc, Basp1, Btbd8, Cacng2, Clstn1, Cyp46a1, Dlg1, Dlg2, Dlgap3, Drd1, Egr3, Gabbr1, Gabra1, Gabrg2, Git1, Grin1, Iqsec3, Itgav, Cnb1, Ccnk1, Lrrc4b, Magi2, Mink1, Mpp7, Ncs1, Nlgn2, Ntrk2, Palm, Psd3, Ptk2, Rimbp2, Rims3, Sh2d5, Shank1, Shank3, Slc30a3, Slc4a10, Slc5a7, Slc6a17, Sptbn2, Strn4, Stx1a, Sv2a, Syndig1, Syngap1, Syngr1, Syt5, Vasp, Znrf1</td><td>53</td><td>3.13</td><td>2.12 × 10<sup>-10</sup></td></tr><tr><td>GO:0005891</td><td>voltage-gated calcium channel complex</td><td>Cacna1g, Cacna2d2, Cacna2d4, Cacnb1, Cacnb3, Cacng2,</td><td>6</td><td>13.1</td><td>3.40 × 10<sup>-2</sup></td></tr><tr><td rowspan="2">Biological process</td><td>GO:0035556</td><td>positive regulation of synaptic transmission, glutamatergic</td><td>Cacng2, Cacng3, Drd1, Grin1, Iqsec2, Nlgn2, Ntrk2, Oxtr, Ptgs2, Shank3, Tnr</td><td>11</td><td>15.6</td><td>2.72 × 10<sup>-6</sup></td></tr><tr><td>GO:0010807</td><td>regulation of synaptic vesicle riming</td><td>Nabp, Stx1a, Stx1b, Stxbp1, Stxbp5</td><td>5</td><td>29.9</td><td>2.90 × 10<sup>-2</sup></td></tr><tr><td rowspan="4">KEGG pathway</td><td>mmu04921</td><td>oxytocin signaling pathway</td><td>Adcy2, Cacna2d2, Cacna2d4, Cacnb1, Cacnb3, Cacng2, Cacng3, Camk1g, Camk2d, Kras, Mylk3, Oxtr, Pik3r6, Prkacb, Ptgs2, Rock2</td><td>16</td><td>5.76</td><td>2.30 × 10<sup>-5</sup></td></tr><tr><td>mmu05414</td><td>dilated cardiomyopathy</td><td>Adcy2, Cacna2d2, Cacna2d4, Cacnb1, Cacnb3, Cacng2, Cacng3, Itgav, Lmna, Prkacb, Tgfb3</td><td>11</td><td>6.44</td><td>1.66 × 10<sup>-3</sup></td></tr><tr><td>mmu04261</td><td>adrenergic signaling in cardiomyocytes</td><td>Adcy2, Atp1a1, Cacna2d2, Cacna2d4, Cacnb1, Cacnb3, Cacng2, Cacng3, Camk2d, Mapk11, Pik3r6, Ppp2r2c, Prkacb</td><td>13</td><td>4.71</td><td>4.43 × 10<sup>-3</sup></td></tr><tr><td>mmu04722</td><td>neurotrophin signaling pathway</td><td>Arhgdia, Arhgdig, Bdnf, Camk2d, Grb2, Kras, Mapk9, Mapk11, Ntrk2, Pik3r2, Sort1</td><td>11</td><td>5.01</td><td>1.52 × 10<sup>-2</sup></td></tr><tr><td rowspan="3">Molecular function</td><td>GO:0005245</td><td>voltage-gated calcium channel activity</td><td>Cacna1g, Cacna2d2, Cacna2d4, Cacnb1, Cacnb3, Cacng2, Cacng3, Itgav, Ncs1</td><td>9</td><td>12.0</td><td>3.93 × 10<sup>-4</sup></td></tr><tr><td>GO:0016301</td><td>kinase activity</td><td>Acvr1c, Brsk2, Camk1g, Camk2d, Cdk5r1, Cdkl5, Dclk1, Grk2, Hkdc1, Ikbkg, Itpkc, Lmtk2, Mapk11, Mapk4, Mapk9, Mark4, Melk, Mink1, Mylk3, Ntrk2, Pak5, Pjkm, Pim1, Pip4k2b, Prkacb, Prkar1b, Ptk2, Rock2, Tesk1, Tlbk1, Tyro3</td><td>31</td><td>2.49</td><td>4.97 × 10<sup>-3</sup></td></tr><tr><td>GO:0004672</td><td>protein kinase activity</td><td>Acvr1c, Brsk2, Camk1g, Camk2d, Cdkl5, Dclk1, Ephb6, Grk2, Lmtk2, Mapk11, Mapk4, Mapk9, Mark4, Melk, Mink1, Mylk3, Ntrk2, Pak5, Pim1, Prkacb, Ptk2, Rock2, Tesk1, Tlbk1, Tyro3</td><td>25</td><td>2.56</td><td>3.01 × 10<sup>-2</sup></td></tr></tbody></table>

Genes belonging to the core signatures (Tables 2 and 3) are in bold.
<table><caption>**Table 5.** Ontological analysis (DAVID) of downregulated genes (consensus and portrait scores < −2) after FLX treatment in stressed animals.</caption><thead><tr><th>Category</th><th>ID</th><th>Name</th><th>Genes</th><th>Count</th><th>Fold Enrichment</th><th>Padj</th></tr></thead><tbody><tr><td rowspan="3">Cellular component</td><td>GO:0022626</td><td>cytosolic ribosome</td><td>Apod, Eif2ak4, Rpl4, Rpl11, Rpl17, Rpl21, Rpl24, Rpl26, Rpl31, Rpl34, Rps6, Rps7, Rps10, Rps11, Rps12, Rps13, Rps15a, Rps24, Rps27a, Uba52</td><td>20</td><td>13.9</td><td>4.75 × 10<sup>−14</sup></td></tr><tr><td>GO:0005840</td><td>ribosome</td><td>Mrpl3, Mrpl50, Mrps35, Rpl4, Rpl11, Rpl17, Rpl21, Rpl24, Rpl31, Rpl34, Rps6, Rps7, Rps10, Rps11, Rps12, Rps13, Rps15a, Rps24, Rps27a, Uba52</td><td>20</td><td>5.76</td><td>9.32 × 10<sup>−7</sup></td></tr><tr><td>GO:0005912</td><td>adherens junction</td><td>Ahi1, Cdh2, Dlg5, Ctnna1, Dll1, Fermt2, Frmd4a, Frmd4b, Pdlim5, Pgm5, Pkp2, Pkp5, Tspan33</td><td>13</td><td>4.61</td><td>1.14 × 10<sup>−2</sup></td></tr><tr><td>Molecular function</td><td>G:0003735</td><td>structural constituent of ribosome</td><td>Mrpl3, Mrps35, Rpl4, Rpl11, Rpl17, Rpl21, Rpl24, Rpl26, Rpl31, Rpl34, Rps6, Rps7, Rps10, Rps11, Rps12, Rps13, Rps15a, Rps24, Rps27a, Uba52</td><td>20</td><td>5.85</td><td>9.53 × 10<sup>−7</sup></td></tr><tr><td>Biological process</td><td>GO:0006412</td><td>translation</td><td>Eif2b2, Eif2s2, Eif3m, Mrpl3, Rpl4, Rpl11, Rpl17, Rpl21, Rpl24, Rpl26, Rpl31, Rpl34, Rps6, Rps7, Rps11, Rps12, Rps13, Rps15a, Rps24, Rps27a, Tars2, Uba52</td><td>22</td><td>3.91</td><td>5.58 × 10<sup>−4</sup></td></tr></tbody></table>

Genes belonging to the core signatures (Tables 2 and 3) are in bold.

<table><caption>**Table 6.** Signature of FLX effect in naïve animals (integration method).</caption><thead><tr><th>Dataset ID<br>Gene</th><th>GSE118669<br>log<sub>2</sub>FC</th><th>p</th><th>GSE42940<br>log<sub>2</sub>FC</th><th>p</th><th>GSE4307<br>log<sub>2</sub>FC</th><th>p</th><th>GSE6476<br>log<sub>2</sub>FC</th><th>p</th><th>GSE84185<br>log<sub>2</sub>FC</th><th>p</th><th>SRP056481<br>log<sub>2</sub>FC</th><th>p</th><th>SRP057486<br>log<sub>2</sub>FC</th><th>p</th><th>Consensus<br>Score</th></tr></thead><tbody><tr><td>Arhgef28</td><td>-0.354</td><td>1.14 × 10<sup>-2</sup></td><td>-0.774</td><td>1.60 × 10<sup>-2</sup></td><td>-0.748</td><td>8.14 × 10<sup>-2</sup></td><td>-0.191</td><td>3.42 × 10<sup>-2</sup></td><td>-1.251</td><td>7.89 × 10<sup>-11</sup></td><td>0.140</td><td>2.15 × 10<sup>-1</sup></td><td>-0.122</td><td>2.92 × 10<sup>-1</sup></td><td>-4</td></tr><tr><td>Arrb2</td><td>-0.179</td><td>1.94 × 10<sup>-2</sup></td><td>-1.311</td><td>1.48 × 10<sup>-3</sup></td><td>-0.184</td><td>5.83 × 10<sup>-1</sup></td><td>-0.156</td><td>5.25 × 10<sup>-2</sup></td><td>-0.250</td><td>2.34 × 10<sup>-3</sup></td><td>-0.013</td><td>9.13 × 10<sup>-1</sup></td><td>-0.227</td><td>4.81 × 10<sup>-2</sup></td><td>-4</td></tr><tr><td>Cd68</td><td>0.377</td><td>3.33 × 10<sup>-4</sup></td><td>-0.257</td><td>1.77 × 10<sup>-1</sup></td><td>0.9468</td><td>1.94 × 10<sup>-1</sup></td><td>0.842</td><td>3.67 × 10<sup>-4</sup></td><td>0.193</td><td>4.49 × 10<sup>-2</sup></td><td>0.0034</td><td>9.49 × 10<sup>-1</sup></td><td>0.352</td><td>1.03 × 10<sup>-2</sup></td><td>4</td></tr><tr><td>Cfh</td><td>0.077</td><td>6.62 × 10<sup>-1</sup></td><td>0.461</td><td>2.91 × 10<sup>-2</sup></td><td>1.3316</td><td>2.28 × 10<sup>-2</sup></td><td>0.430</td><td>8.38 × 10<sup>-4</sup></td><td>0.027</td><td>3.68 × 10<sup>-1</sup></td><td>0.054</td><td>5.75 × 10<sup>-1</sup></td><td>0.392</td><td>1.96 × 10<sup>-3</sup></td><td>4</td></tr><tr><td>Cdon</td><td>-0.274</td><td>1.76 × 10<sup>-2</sup></td><td>-0.847</td><td>4.12 × 10<sup>-2</sup></td><td>-3.139</td><td>3.14 × 10<sup>-3</sup></td><td>-0.359</td><td>7.20 × 10<sup>-2</sup></td><td>-0.607</td><td>1.71 × 10<sup>-4</sup></td><td>-0.103</td><td>3.78 × 10<sup>-1</sup></td><td>-0.031</td><td>8.19 × 10<sup>-1</sup></td><td>-4</td></tr><tr><td>Cldb</td><td>0.227</td><td>3.03 × 10<sup>-2</sup></td><td>-0.751</td><td>4.05 × 10<sup>-2</sup></td><td>0.507</td><td>3.48 × 10<sup>-1</sup></td><td>0.396</td><td>1.27 × 10<sup>-3</sup></td><td>-0.376</td><td>1.84 × 10<sup>-5</sup></td><td>-0.249</td><td>3.94 × 10<sup>-2</sup></td><td>0.111</td><td>3.34 × 10<sup>-1</sup></td><td>-1</td></tr><tr><td>Ddr1</td><td>0.245</td><td>9.08 × 10<sup>-4</sup></td><td>-0.276</td><td>2.83 × 10<sup>-1</sup></td><td>1.145</td><td>5.00 × 10<sup>-2</sup></td><td>0.261</td><td>8.30 × 10<sup>-4</sup></td><td>-0.493</td><td>2.22 × 10<sup>-5</sup></td><td>0.097</td><td>4.09 × 10<sup>-1</sup></td><td>0.342</td><td>2.84 × 10<sup>-3</sup></td><td>3</td></tr><tr><td>Doc2b</td><td>-0.506</td><td>1.51 × 10<sup>-3</sup></td><td>0.371</td><td>1.83 × 10<sup>-1</sup></td><td>-0.798</td><td>7.63 × 10<sup>-2</sup></td><td>-0.233</td><td>1.07 × 10<sup>-2</sup></td><td>-3.232</td><td>3.42 × 10<sup>-11</sup></td><td>0.028</td><td>7.90 × 10<sup>-1</sup></td><td>-0.416</td><td>4.12 × 10<sup>-4</sup></td><td>-4</td></tr><tr><td>Fat4</td><td>-0.890</td><td>4.54 × 10<sup>-4</sup></td><td>0.000</td><td>1.00</td><td>0.000</td><td>1.00</td><td>-1.352</td><td>1.80 × 10<sup>-3</sup></td><td>-1.748</td><td>2.48 × 10<sup>-11</sup></td><td>-0.210</td><td>6.03 × 10<sup>-2</sup></td><td>-0.342</td><td>1.25 × 10<sup>-2</sup></td><td>-4</td></tr><tr><td>Gpr12</td><td>-0.402</td><td>3.06 × 10<sup>-2</sup></td><td>-0.708</td><td>2.38 × 10<sup>-2</sup></td><td>-1.981</td><td>1.71 × 10<sup>-1</sup></td><td>-0.837</td><td>2.38 × 10<sup>-4</sup></td><td>-0.491</td><td>1.25 × 10<sup>-6</sup></td><td>-0.115</td><td>3.16 × 10<sup>-1</sup></td><td>-0.001</td><td>9.95 × 10<sup>-1</sup></td><td>-4</td></tr><tr><td>Gst1</td><td>0.303</td><td>2.82 × 10<sup>-2</sup></td><td>0.783</td><td>1.49 × 10<sup>-1</sup></td><td>1.488</td><td>7.55 × 10<sup>-2</sup></td><td>0.549</td><td>4.24 × 10<sup>-3</sup></td><td>0.397</td><td>4.04 × 10<sup>-6</sup></td><td>0.106</td><td>3.13 × 10<sup>-1</sup></td><td>0.495</td><td>1.16 × 10<sup>-4</sup></td><td>4</td></tr><tr><td>H2-D1</td><td>0.210</td><td>9.67 × 10<sup>-2</sup></td><td>0.000</td><td>1.00</td><td>1.173</td><td>3.25 × 10<sup>-2</sup></td><td>0.673</td><td>1.07 × 10<sup>-2</sup></td><td>0.424</td><td>3.22 × 10<sup>-4</sup></td><td>0.000</td><td>1.00</td><td>0.425</td><td>3.57 × 10<sup>-4</sup></td><td>4</td></tr><tr><td>Homer1</td><td>0.174</td><td>4.03 × 10<sup>-2</sup></td><td>0.688</td><td>5.84 × 10<sup>-3</sup></td><td>0.873</td><td>6.39 × 10<sup>-2</sup></td><td>2.919</td><td>1.41 × 10<sup>-3</sup></td><td>1.133</td><td>7.05 × 10<sup>-12</sup></td><td>0.057</td><td>6.17 × 10<sup>-1</sup></td><td>-0.307</td><td>2.55 × 10<sup>-2</sup></td><td>3</td></tr><tr><td>Htr1b</td><td>0.121</td><td>3.19 × 10<sup>-2</sup></td><td>0.544</td><td>9.07 × 10<sup>-3</sup></td><td>1.466</td><td>1.88 × 10<sup>-2</sup></td><td>-1.189</td><td>1.03 × 10<sup>-3</sup></td><td>1.084</td><td>2.35 × 10<sup>-3</sup></td><td>0.076</td><td>3.90 × 10<sup>-1</sup></td><td>-0.262</td><td>5.54 × 10<sup>-1</sup></td><td>3</td></tr><tr><td>Htr5b</td><td>0.920</td><td>3.29 × 10<sup>-3</sup></td><td>0.987</td><td>7.01 × 10<sup>-2</sup></td><td>0.847</td><td>7.39 × 10<sup>-1</sup></td><td>-1.112</td><td>3.30 × 10<sup>-2</sup></td><td>1.252</td><td>7.36 × 10<sup>-3</sup></td><td>0.224</td><td>4.37 × 10<sup>-2</sup></td><td>-0.477</td><td>4.97 × 10<sup>-4</sup></td><td>1</td></tr><tr><td>Igfbp6</td><td>1.110</td><td>2.04 × 10<sup>-4</sup></td><td>0.081</td><td>6.67 × 10<sup>-1</sup></td><td>2.796</td><td>2.01 × 10<sup>-1</sup></td><td>0.821</td><td>1.64 × 10<sup>-2</sup></td><td>2.720</td><td>3.41 × 10<sup>-3</sup></td><td>-0.079</td><td>4.68 × 10<sup>-1</sup></td><td>0.325</td><td>6.21 × 10<sup>-3</sup></td><td>4</td></tr><tr><td>Ints10</td><td>-0.129</td><td>2.54 × 10<sup>-2</sup></td><td>0.000</td><td>1.00</td><td>0.000</td><td>1.00</td><td>-0.562</td><td>6.79 × 10<sup>-4</sup></td><td>-0.259</td><td>2.41 × 10<sup>-3</sup></td><td>0.016</td><td>8.92 × 10<sup>-1</sup></td><td>-0.254</td><td>1.42 × 10<sup>-2</sup></td><td>-4</td></tr><tr><td>Isoc1</td><td>-0.460</td><td>6.82 × 10<sup>-5</sup></td><td>-0.337</td><td>4.55 × 10<sup>-2</sup></td><td>-1.803</td><td>1.54 × 10<sup>-2</sup></td><td>-0.360</td><td>2.98 × 10<sup>-2</sup></td><td>-0.761</td><td>1.26 × 10<sup>-7</sup></td><td>-0.051</td><td>4.85 × 10<sup>-1</sup></td><td>0.035</td><td>7.98 × 10<sup>-1</sup></td><td>-5</td></tr><tr><td>Ilg4</td><td>-0.437</td><td>1.33 × 10<sup>-2</sup></td><td>-0.007</td><td>9.83 × 10<sup>-1</sup></td><td>-3.476</td><td>6.77 × 10<sup>-2</sup></td><td>-0.601</td><td>2.59 × 10<sup>-2</sup></td><td>-1.608</td><td>1.14 × 10<sup>-10</sup></td><td>-0.012</td><td>8.44 × 10<sup>-1</sup></td><td>-0.490</td><td>2.28 × 10<sup>-5</sup></td><td>-4</td></tr><tr><td>Itsn1</td><td>-0.238</td><td>3.80 × 10<sup>-3</sup></td><td>-2.471</td><td>1.91 × 10<sup>-3</sup></td><td>-0.915</td><td>5.91 × 10<sup>-2</sup></td><td>-0.232</td><td>2.93 × 10<sup>-3</sup></td><td>-0.785</td><td>1.52 × 10<sup>-8</sup></td><td>-0.068</td><td>5.71 × 10<sup>-1</sup></td><td>-0.164</td><td>1.84 × 10<sup>-1</sup></td><td>-4</td></tr><tr><td>Kirrel3</td><td>-0.264</td><td>1.02 × 10<sup>-3</sup></td><td>-0.603</td><td>4.93 × 10<sup>-2</sup></td><td>0.000</td><td>1.00</td><td>-0.447</td><td>5.77 × 10<sup>-3</sup></td><td>-0.617</td><td>9.80 × 10<sup>-8</sup></td><td>-0.105</td><td>3.65 × 10<sup>-1</sup></td><td>0.030</td><td>8.29 × 10<sup>-1</sup></td><td>-4</td></tr><tr><td>Knstrn</td><td>0.182</td><td>2.56 × 10<sup>-2</sup></td><td>0.000</td><td>1.00</td><td>3.322</td><td>2.93 × 10<sup>-3</sup></td><td>0.153</td><td>1.00 × 10<sup>-2</sup></td><td>1.627</td><td>1.10 × 10<sup>-11</sup></td><td>0.002</td><td>9.38 × 10<sup>-1</sup></td><td>0.026</td><td>8.34 × 10<sup>-1</sup></td><td>4</td></tr><tr><td>Map1a</td><td>-0.146</td><td>1.24 × 10<sup>-2</sup></td><td>-3.255</td><td>4.42 × 10<sup>-3</sup></td><td>0.000</td><td>1.00</td><td>-0.116</td><td>2.02 × 10<sup>-2</sup></td><td>-0.199</td><td>2.12 × 10<sup>-3</sup></td><td>-0.059</td><td>6.26 × 10<sup>-1</sup></td><td>-0.273</td><td>1.40 × 10<sup>-2</sup></td><td>-5</td></tr><tr><td>Mpdz</td><td>-0.355</td><td>9.19 × 10<sup>-5</sup></td><td>-0.074</td><td>7.18 × 10<sup>-1</sup></td><td>-0.226</td><td>5.95 × 10<sup>-1</sup></td><td>-0.323</td><td>2.64 × 10<sup>-2</sup></td><td>-0.540</td><td>6.63 × 10<sup>-5</sup></td><td>-0.071</td><td>5.55 × 10<sup>-1</sup></td><td>-0.250</td><td>4.02 × 10<sup>-2</sup></td><td>-4</td></tr><tr><td>Mat2a</td><td>0.195</td><td>6.33 × 10<sup>-3</sup></td><td>0.689</td><td>3.65 × 10<sup>-3</sup></td><td>1.323</td><td>2.40 × 10<sup>-2</sup></td><td>-0.210</td><td>1.34 × 10<sup>-2</sup></td><td>-0.398</td><td>8.29 × 10<sup>-5</sup></td><td>0.046</td><td>6.98 × 10<sup>-1</sup></td><td>0.014</td><td>8.81 × 10<sup>-1</sup></td><td>1</td></tr><tr><td>Mylk</td><td>0.333</td><td>3.40 × 10<sup>-2</sup></td><td>-0.616</td><td>2.80 × 10<sup>-3</sup></td><td>1.305</td><td>2.24 × 10<sup>-1</sup></td><td>0.792</td><td>9.15 × 10<sup>-3</sup></td><td>-1.215</td><td>1.07 × 10<sup>-7</sup></td><td>0.034</td><td>7.68 × 10<sup>-1</sup></td><td>0.455</td><td>8.35 × 10<sup>-4</sup></td><td>1</td></tr><tr><td>Myo1e</td><td>0.203</td><td>9.58 × 10<sup>-3</sup></td><td>0.540</td><td>1.79 × 10<sup>-1</sup></td><td>1.855</td><td>1.08 × 10<sup>-2</sup></td><td>0.054</td><td>3.40 × 10<sup>-2</sup></td><td>0.837</td><td>2.74 × 10<sup>-6</sup></td><td>-0.005</td><td>9.59 × 10<sup>-1</sup></td><td>0.290</td><td>2.78 × 10<sup>-4</sup></td><td>5</td></tr><tr><td>Negr1</td><td>-0.255</td><td>9.34 × 10<sup>-3</sup></td><td>0.472</td><td>5.96 × 10<sup>-1</sup></td><td>0.000</td><td>1.00</td><td>-0.326</td><td>3.76 × 10<sup>-3</sup></td><td>-0.588</td><td>3.19 × 10<sup>-3</sup></td><td>-0.246</td><td>3.52 × 10<sup>-2</sup></td><td>0.023</td><td>8.60 × 10<sup>-1</sup></td><td>-4</td></tr></tbody></table>
<table><caption>Table 6. *Cont.*</caption><thead><tr><th rowspan="2">Dataset ID<br>Gene</th><th colspan="2">GSE118669</th><th colspan="2">GSE42940</th><th colspan="2">GSE4307</th><th colspan="2">GSE6476</th><th colspan="2">GSE84185</th><th colspan="2">SRP056481</th><th colspan="2">SRP057486</th><th rowspan="2">Consensus<br>Score</th></tr><tr><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th></tr></thead><tbody><tr><td>Nhs12</td><td>-0.427</td><td>5.59 × 10<sup>-4</sup></td><td>0.000</td><td>1.00</td><td>-0.671</td><td>3.28 × 10<sup>-1</sup></td><td>-0.356</td><td>1.84 × 10<sup>-2</sup></td><td>-0.868</td><td>5.61 × 10<sup>-10</sup></td><td>0.000</td><td>1.00</td><td>-0.302</td><td>6.20 × 10<sup>-3</sup></td><td>-4</td></tr><tr><td>Ntrk3</td><td>-0.253</td><td>3.34 × 10<sup>-3</sup></td><td>0.179</td><td>2.81 × 10<sup>-1</sup></td><td>1.932</td><td>5.89 × 10<sup>-2</sup></td><td>-0.156</td><td>5.25 × 10<sup>-3</sup></td><td>-0.414</td><td>2.52 × 10<sup>-6</sup></td><td>-0.059</td><td>6.23 × 10<sup>-1</sup></td><td>-0.409</td><td>1.31 × 10<sup>-4</sup></td><td>-4</td></tr><tr><td>Pdsh7</td><td>0.355</td><td>2.30 × 10<sup>-4</sup></td><td>0.498</td><td>3.95 × 10<sup>-1</sup></td><td>1.132</td><td>4.86 × 10<sup>-2</sup></td><td>0.222</td><td>4.79 × 10<sup>-2</sup></td><td>1.167</td><td>1.04 × 10<sup>-7</sup></td><td>-0.114</td><td>3.09 × 10<sup>-1</sup></td><td>0.204</td><td>9.10 × 10<sup>-2</sup></td><td>4</td></tr><tr><td>Pdth19</td><td>-0.572</td><td>4.44 × 10<sup>-3</sup></td><td>0.000</td><td>1.00</td><td>0.000</td><td>1.00</td><td>-0.252</td><td>1.31 × 10<sup>-2</sup></td><td>-1.419</td><td>9.69 × 10<sup>-15</sup></td><td>-0.114</td><td>3.35 × 10<sup>-1</sup></td><td>-0.413</td><td>1.08 × 10<sup>-4</sup></td><td>-4</td></tr><tr><td>Pde4b</td><td>0.427</td><td>8.61 × 10<sup>-6</sup></td><td>-0.064</td><td>7.55 × 10<sup>-1</sup></td><td>0.763</td><td>3.36 × 10<sup>-1</sup></td><td>0.180</td><td>2.91 × 10<sup>-2</sup></td><td>0.485</td><td>7.66 × 10<sup>-7</sup></td><td>0.036</td><td>7.31 × 10<sup>-1</sup></td><td>0.517</td><td>6.00 × 10<sup>-7</sup></td><td>4</td></tr><tr><td>Pde7b</td><td>-0.720</td><td>6.86 × 10<sup>-4</sup></td><td>0.063</td><td>8.82 × 10<sup>-1</sup></td><td>-1.251</td><td>3.96 × 10<sup>-2</sup></td><td>-0.587</td><td>1.96 × 10<sup>-2</sup></td><td>0.511</td><td>1.55 × 10<sup>-3</sup></td><td>-0.287</td><td>8.38 × 10<sup>-3</sup></td><td>-0.054</td><td>6.10 × 10<sup>-1</sup></td><td>-3</td></tr><tr><td>Pdlim5</td><td>-0.124</td><td>2.95 × 10<sup>-2</sup></td><td>-0.722</td><td>8.51 × 10<sup>-3</sup></td><td>0.900</td><td>7.04 × 10<sup>-2</sup></td><td>0.459</td><td>9.00 × 10<sup>-3</sup></td><td>-2.857</td><td>2.37 × 10<sup>-12</sup></td><td>-0.165</td><td>9.09 × 10<sup>-2</sup></td><td>0.267</td><td>3.40 × 10<sup>-2</sup></td><td>-1</td></tr><tr><td>Rab27a</td><td>-0.474</td><td>1.50 × 10<sup>-3</sup></td><td>0.967</td><td>1.82 × 10<sup>-1</sup></td><td>-2.928</td><td>5.78 × 10<sup>-3</sup></td><td>-0.680</td><td>2.90 × 10<sup>-2</sup></td><td>-1.773</td><td>8.63 × 10<sup>-10</sup></td><td>-0.047</td><td>6.45 × 10<sup>-1</sup></td><td>-0.096</td><td>4.85 × 10<sup>-1</sup></td><td>-4</td></tr><tr><td>Rasgrf1</td><td>-0.394</td><td>5.36 × 10<sup>-5</sup></td><td>-0.200</td><td>2.02 × 10<sup>-1</sup></td><td>-0.443</td><td>3.55 × 10<sup>-1</sup></td><td>-0.519</td><td>4.75 × 10<sup>-4</sup></td><td>-1.051</td><td>2.38 × 10<sup>-10</sup></td><td>-0.006</td><td>9.59 × 10<sup>-1</sup></td><td>-0.335</td><td>5.16 × 10<sup>-4</sup></td><td>-4</td></tr><tr><td>Rassf5</td><td>0.281</td><td>8.67 × 10<sup>-5</sup></td><td>-1.470</td><td>2.91 × 10<sup>-1</sup></td><td>1.934</td><td>9.78 × 10<sup>-3</sup></td><td>0.156</td><td>8.56 × 10<sup>-3</sup></td><td>0.398</td><td>6.61 × 10<sup>-4</sup></td><td>-0.004</td><td>9.66 × 10<sup>-1</sup></td><td>0.225</td><td>8.59 × 10<sup>-2</sup></td><td>4</td></tr><tr><td>S100a6</td><td>0.473</td><td>4.52 × 10<sup>-3</sup></td><td>-0.668</td><td>1.93 × 10<sup>-2</sup></td><td>2.020</td><td>1.08 × 10<sup>-1</sup></td><td>2.247</td><td>3.72 × 10<sup>-4</sup></td><td>1.012</td><td>4.01 × 10<sup>-11</sup></td><td>-0.027</td><td>7.80 × 10<sup>-1</sup></td><td>0.368</td><td>7.22 × 10<sup>-3</sup></td><td>3</td></tr><tr><td>Scn3b</td><td>-0.245</td><td>7.02 × 10<sup>-3</sup></td><td>-0.479</td><td>1.05 × 10<sup>-2</sup></td><td>-0.678</td><td>1.15 × 10<sup>-1</sup></td><td>-0.594</td><td>5.20 × 10<sup>-5</sup></td><td>-0.971</td><td>4.31 × 10<sup>-8</sup></td><td>0.178</td><td>1.41 × 10<sup>-1</sup></td><td>-0.304</td><td>2.13 × 10<sup>-2</sup></td><td>-5</td></tr><tr><td>Sell13</td><td>0.281</td><td>8.20 × 10<sup>-4</sup></td><td>-1.470</td><td>1.00</td><td>1.934</td><td>3.86 × 10<sup>-2</sup></td><td>0.156</td><td>1.18 × 10<sup>-2</sup></td><td>0.398</td><td>9.19 × 10<sup>-7</sup></td><td>-0.004</td><td>2.82 × 10<sup>-1</sup></td><td>0.225</td><td>8.62 × 10<sup>-1</sup></td><td>4</td></tr><tr><td>Sema3a</td><td>0.948</td><td>6.47 × 10<sup>-4</sup></td><td>-0.388</td><td>4.91 × 10<sup>-1</sup></td><td>2.727</td><td>3.16 × 10<sup>-2</sup></td><td>1.152</td><td>4.04 × 10<sup>-3</sup></td><td>4.176</td><td>2.62 × 10<sup>-11</sup></td><td>-0.004</td><td>9.07 × 10<sup>-1</sup></td><td>-0.084</td><td>5.41 × 10<sup>-1</sup></td><td>4</td></tr><tr><td>Sorc1</td><td>1.183</td><td>1.18 × 10<sup>-3</sup></td><td>-0.366</td><td>3.72 × 10<sup>-1</sup></td><td>1.372</td><td>9.98 × 10<sup>-2</sup></td><td>1.300</td><td>1.04 × 10<sup>-2</sup></td><td>2.102</td><td>8.85 × 10<sup>-8</sup></td><td>0.059</td><td>5.70 × 10<sup>-1</sup></td><td>0.495</td><td>2.66 × 10<sup>-4</sup></td><td>4</td></tr><tr><td>Tfrc</td><td>0.269</td><td>6.62 × 10<sup>-4</sup></td><td>0.610</td><td>3.00 × 10<sup>-1</sup></td><td>0.795</td><td>2.65 × 10<sup>-1</sup></td><td>0.725</td><td>7.75 × 10<sup>-3</sup></td><td>0.435</td><td>2.85 × 10<sup>-6</sup></td><td>-0.103</td><td>3.92 × 10<sup>-1</sup></td><td>0.298</td><td>8.19 × 10<sup>-3</sup></td><td>4</td></tr><tr><td>Tnxb</td><td>-0.323</td><td>4.79 × 10<sup>-4</sup></td><td>0.000</td><td>1.00</td><td>-0.695</td><td>3.40 × 10<sup>-1</sup></td><td>-0.574</td><td>7.37 × 10<sup>-3</sup></td><td>-1.462</td><td>6.94 × 10<sup>-9</sup></td><td>0.000</td><td>1.00</td><td>-0.398</td><td>3.57 × 10<sup>-3</sup></td><td>-4</td></tr><tr><td>Zfxt2</td><td>-0.145</td><td>1.02 × 10<sup>-2</sup></td><td>0.000</td><td>1.00</td><td>0.000</td><td>1.00</td><td>-0.051</td><td>3.74 × 10<sup>-2</sup></td><td>-0.653</td><td>2.51 × 10<sup>-9</sup></td><td>0.035</td><td>7.70 × 10<sup>-1</sup></td><td>-0.479</td><td>9.56 × 10<sup>-5</sup></td><td>-4</td></tr><tr><td>Zfp316</td><td>-0.218</td><td>4.28 × 10<sup>-3</sup></td><td>-1.324</td><td>3.00 × 10<sup>-4</sup></td><td>-1.531</td><td>1.80 × 10<sup>-2</sup></td><td>-0.912</td><td>8.63 × 10<sup>-3</sup></td><td>-0.448</td><td>1.30 × 10<sup>-3</sup></td><td>-0.107</td><td>3.35 × 10<sup>-1</sup></td><td>0.019</td><td>8.72 × 10<sup>-1</sup></td><td>-5</td></tr></tbody></table>

$|\log_2 \text{FC}| > 0.263$ and $p$-values $\le 0.05$ are in bold. Consensus score indicates whether FLX induces mainly upregulation (positive score) or downregulation (negative score) of gene expression. Genes with absolute value of consensus score equal to 5 are in bold.

<table><caption>Table 7. Signature of FLX effect in naïve animals (portrait score method).</caption><thead><tr><th rowspan="2">Dataset ID<br>Gene</th><th colspan="2">GSE118669</th><th colspan="2">GSE42940</th><th colspan="2">GSE4307</th><th colspan="2">GSE6476</th><th colspan="2">GSE84185</th><th colspan="2">SRP056481</th><th colspan="2">SRP057486</th><th rowspan="2">Consensus Score</th><th rowspan="2">Portrait Score</th></tr><tr><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th></tr></thead><tbody><tr><td>Vip</td><td>0.306</td><td>3.19 × 10<sup>-2</sup></td><td>0.141</td><td>5.63 × 10<sup>-1</sup></td><td>0.934</td><td>9.13 × 10<sup>-2</sup></td><td>0.002</td><td>9.84 × 10<sup>-1</sup></td><td>2.129</td><td>6.92 × 10<sup>-9</sup></td><td>0.122</td><td>2.30 × 10<sup>-1</sup></td><td>0.318</td><td>1.61 × 10<sup>-2</sup></td><td>3</td><td>55,566,665</td></tr><tr><td>Gsn</td><td>0.303</td><td>2.82 × 10<sup>-2</sup></td><td>0.783</td><td>1.49 × 10<sup>-1</sup></td><td>1.49</td><td>7.55 × 10<sup>-2</sup></td><td>0.549</td><td>4.24 × 10<sup>-3</sup></td><td>0.397</td><td>4.04 × 10<sup>-6</sup></td><td>0.106</td><td>3.13 × 10<sup>-1</sup></td><td>0.495</td><td>1.16 × 10<sup>-4</sup></td><td>4</td><td>46,777,777</td></tr><tr><td>Cd9</td><td>0.209</td><td>8.41 × 10<sup>-2</sup></td><td>0.458</td><td>2.96 × 10<sup>-2</sup></td><td>1.346</td><td>1.53 × 10<sup>-1</sup></td><td>0.764</td><td>3.18 × 10<sup>-3</sup></td><td>-0.168</td><td>5.47 × 10<sup>-2</sup></td><td>0.205</td><td>8.59 × 10<sup>-2</sup></td><td>0.426</td><td>1.46 × 10<sup>-3</sup></td><td>3</td><td>46,666,655</td></tr><tr><td>S100a10</td><td>0.203</td><td>1.55 × 10<sup>-1</sup></td><td>0.308</td><td>8.39 × 10<sup>-2</sup></td><td>1.600</td><td>5.32 × 10<sup>-2</sup></td><td>0.619</td><td>1.09 × 10<sup>-2</sup></td><td>1.107</td><td>9.75 × 10<sup>-9</sup></td><td>0.068</td><td>5.22 × 10<sup>-1</sup></td><td>0.258</td><td>5.77 × 10<sup>-2</sup></td><td>2</td><td>45,777,777</td></tr><tr><td>Fcgr2b</td><td>0.142</td><td>5.86 × 10<sup>-2</sup></td><td>0.688</td><td>3.57 × 10<sup>-2</sup></td><td>1.566</td><td>3.08 × 10<sup>-2</sup></td><td>0.920</td><td>1.59 × 10<sup>-2</sup></td><td>0.617</td><td>7.66 × 10<sup>-6</sup></td><td>0.004</td><td>9.53 × 10<sup>-1</sup></td><td>0.390</td><td>3.44 × 10<sup>-3</sup></td><td>5</td><td>45,666,666</td></tr><tr><td>Vsll1</td><td>0.212</td><td>4.23 × 10<sup>-2</sup></td><td>0.596</td><td>3.87 × 10<sup>-2</sup></td><td>0.408</td><td>2.37 × 10<sup>-1</sup></td><td>-0.084</td><td>4.32 × 10<sup>-1</sup></td><td>1.372</td><td>4.34 × 10<sup>-12</sup></td><td>0.220</td><td>6.71 × 10<sup>-2</sup></td><td>0.205</td><td>1.23 × 10<sup>-1</sup></td><td>3</td><td>45,666,555</td></tr><tr><td>Bgn</td><td>0.144</td><td>4.83 × 10<sup>-1</sup></td><td>0.495</td><td>4.41 × 10<sup>-2</sup></td><td>1.739</td><td>6.86 × 10<sup>-2</sup></td><td>1.368</td><td>9.10 × 10<sup>-3</sup></td><td>0.377</td><td>3.47 × 10<sup>-2</sup></td><td>0.100</td><td>3.76 × 10<sup>-1</sup></td><td>0.302</td><td>1.76 × 10<sup>-2</sup></td><td>4</td><td>45,555,777</td></tr><tr><td>Sell13</td><td>0.281</td><td>8.20 × 10<sup>-4</sup></td><td>-1.470</td><td>1.00</td><td>1.934</td><td>3.86 × 10<sup>-2</sup></td><td>0.156</td><td>1.18 × 10<sup>-2</sup></td><td>0.398</td><td>9.19 × 10<sup>-7</sup></td><td>-0.004</td><td>2.82 × 10<sup>-1</sup></td><td>0.225</td><td>8.62 × 10<sup>-1</sup></td><td>4</td><td>45,555,555</td></tr><tr><td>Ppp2r5c</td><td>0.392</td><td>8.26 × 10<sup>-4</sup></td><td>-0.289</td><td>2.68 × 10<sup>-1</sup></td><td>1.23</td><td>5.69 × 10<sup>-2</sup></td><td>0.216</td><td>7.36 × 10<sup>-3</sup></td><td>1.25</td><td>3.61 × 10<sup>-10</sup></td><td>0.109</td><td>3.63 × 10<sup>-1</sup></td><td>0.040</td><td>7.18 × 10<sup>-1</sup></td><td>3</td><td>45,444,444</td></tr><tr><td>Ddr1</td><td>0.245</td><td>9.08 × 10<sup>-4</sup></td><td>-0.276</td><td>2.83 × 10<sup>-1</sup></td><td>1.145</td><td>5.00 × 10<sup>-2</sup></td><td>0.261</td><td>8.30 × 10<sup>-4</sup></td><td>-0.493</td><td>2.22 × 10<sup>-5</sup></td><td>0.097</td><td>4.09 × 10<sup>-1</sup></td><td>0.342</td><td>2.84 × 10<sup>-3</sup></td><td>3</td><td>45,333,333</td></tr><tr><td>Anxa5</td><td>0.109</td><td>1.86 × 10<sup>-1</sup></td><td>0.678</td><td>5.12 × 10<sup>-2</sup></td><td>0.981</td><td>7.99 × 10<sup>-2</sup></td><td>0.641</td><td>2.84 × 10<sup>-3</sup></td><td>0.259</td><td>3.55 × 10<sup>-3</sup></td><td>0.083</td><td>4.81 × 10<sup>-1</sup></td><td>0.314</td><td>9.22 × 10<sup>-3</sup></td><td>3</td><td>44,667,777</td></tr><tr><td>Spo23</td><td>0.375</td><td>1.77 × 10<sup>-2</sup></td><td>0.287</td><td>3.05 × 10<sup>-1</sup></td><td>1.462</td><td>2.78 × 10<sup>-2</sup></td><td>0.419</td><td>3.12 × 10<sup>-3</sup></td><td>0.362</td><td>1.12 × 10<sup>-2</sup></td><td>0.053</td><td>6.58 × 10<sup>-1</sup></td><td>0.317</td><td>1.33 × 10<sup>-2</sup></td><td>5</td><td>44,557,777</td></tr><tr><td>Lcam1</td><td>0.027</td><td>7.66 × 10<sup>-1</sup></td><td>0.640</td><td>3.06 × 10<sup>-1</sup></td><td>0.869</td><td>5.82 × 10<sup>-2</sup></td><td>0.164</td><td>6.84 × 10<sup>-3</sup></td><td>-0.008</td><td>5.96 × 10<sup>-1</sup></td><td>0.110</td><td>2.43 × 10<sup>-1</sup></td><td>0.311</td><td>1.19 × 10<sup>-2</sup></td><td>2</td><td>44,555,556</td></tr><tr><td>Rassf8</td><td>0.382</td><td>1.38 × 10<sup>-2</sup></td><td>-0.501</td><td>6.59 × 10<sup>-1</sup></td><td>1.424</td><td>7.20 × 10<sup>-2</sup></td><td>0.577</td><td>2.20 × 10<sup>-3</sup></td><td>0.971</td><td>1.39 × 10<sup>-5</sup></td><td>0.026</td><td>8.17 × 10<sup>-1</sup></td><td>0.279</td><td>3.82 × 10<sup>-2</sup></td><td>4</td><td>44,554,466</td></tr><tr><td>Tspan5</td><td>0.231</td><td>1.90 × 10<sup>-1</sup></td><td>1.609</td><td>1.56 × 10<sup>-2</sup></td><td>1.041</td><td>8.04 × 10<sup>-2</sup></td><td>0.378</td><td>3.79 × 10<sup>-3</sup></td><td>0.787</td><td>2.00 × 10<sup>-8</sup></td><td>0.008</td><td>9.47 × 10<sup>-1</sup></td><td>-0.108</td><td>3.79 × 10<sup>-1</sup></td><td>3</td><td>44,554,444</td></tr><tr><td>Kif5b</td><td>0.218</td><td>3.85 × 10<sup>-3</sup></td><td>0.699</td><td>4.29 × 10<sup>-2</sup></td><td>0.350</td><td>2.86 × 10<sup>-1</sup></td><td>0.055</td><td>4.41 × 10<sup>-1</sup></td><td>0.586</td><td>1.36 × 10<sup>-8</sup></td><td>0.012</td><td>9.17 × 10<sup>-1</sup></td><td>0.257</td><td>1.03 × 10<sup>-2</sup></td><td>4</td><td>44,455,557</td></tr><tr><td>Cft</td><td>0.077</td><td>6.62 × 10<sup>-1</sup></td><td>0.461</td><td>2.91 × 10<sup>-2</sup></td><td>1.3316</td><td>2.28 × 10<sup>-2</sup></td><td>0.430</td><td>8.38 × 10<sup>-4</sup></td><td>0.027</td><td>3.68 × 10<sup>-1</sup></td><td>0.054</td><td>5.75 × 10<sup>-1</sup></td><td>0.392</td><td>1.96 × 10<sup>-3</sup></td><td>4</td><td>44,455,556</td></tr><tr><td>Homer1</td><td>0.174</td><td>4.03 × 10<sup>-2</sup></td><td>0.688</td><td>5.84 × 10<sup>-3</sup></td><td>0.873</td><td>6.39 × 10<sup>-2</sup></td><td>2.919</td><td>1.41 × 10<sup>-3</sup></td><td>1.133</td><td>7.05 × 10<sup>-12</sup></td><td>0.057</td><td>6.17 × 10<sup>-1</sup></td><td>-0.307</td><td>2.55 × 10<sup>-2</sup></td><td>3</td><td>44,455,555</td></tr><tr><td>Tmem47</td><td>0.346</td><td>2.20 × 10<sup>-3</sup></td><td>0.000</td><td>1.00</td><td>1.247</td><td>4.33 × 10<sup>-2</sup></td><td>0.546</td><td>1.13 × 10<sup>-2</sup></td><td>0.000</td><td>1.00</td><td>0.050</td><td>6.67 × 10<sup>-1</sup></td><td>0.250</td><td>2.10 × 10<sup>-2</sup></td><td>4</td><td>44,445,555</td></tr><tr><td>Tmem98</td><td>0.190</td><td>2.81 × 10<sup>-2</sup></td><td>0.000</td><td>1.00</td><td>0.753</td><td>7.81 × 10<sup>-2</sup></td><td>0.312</td><td>2.47 × 10<sup>-1</sup></td><td>-0.222</td><td>2.87 × 10<sup>-2</sup></td><td>0.184</td><td>6.86 × 10<sup>-2</sup></td><td>0.655</td><td>1.19 × 10<sup>-6</sup></td><td>1</td><td>44,445,444</td></tr><tr><td>Knstrn</td><td>0.182</td><td>2.56 × 10<sup>-2</sup></td><td>0.000</td><td>1.00</td><td>3.323</td><td>2.93 × 10<sup>-3</sup></td><td>0.153</td><td>1.00 × 10<sup>-2</sup></td><td>1.627</td><td>1.10 × 10<sup>-11</sup></td><td>0.002</td><td>9.38 × 10<sup>-1</sup></td><td>0.026</td><td>8.34 × 10<sup>-1</sup></td><td>4</td><td>44,444,444</td></tr><tr><td>Sema3a</td><td>0.948</td><td>6.47 × 10<sup>-4</sup></td><td>-0.388</td><td>4.91 × 10<sup>-1</sup></td><td>2.727</td><td>3.16 × 10<sup>-2</sup></td><td>1.152</td><td>4.04 × 10<sup>-3</sup></td><td>4.176</td><td>2.62 × 10<sup>-11</sup></td><td>-0.004</td><td>9.07 × 10<sup>-1</sup></td><td>-0.084</td><td>5.41 × 10<sup>-1</sup></td><td>4</td><td>44,443,312</td></tr><tr><td>C1qb</td><td>0.240</td><td>5.04 × 10<sup>-4</sup></td><td>-0.180</td><td>3.60 × 10<sup>-1</sup></td><td>1.00</td><td>4.37 × 10<sup>-2</sup></td><td>0.873</td><td>1.59 × 10<sup>-4</sup></td><td>-0.011</td><td>8.42 × 10<sup>-1</sup></td><td>0.030</td><td>7.94 × 10<sup>-1</sup></td><td>0.473</td><td>2.69 × 10<sup>-5</sup></td><td>4</td><td>44,433,344</td></tr><tr><td>Dpp4</td><td>-0.111</td><td>1.13 × 10<sup>-1</sup></td><td>0.329</td><td>9.32 × 10<sup>-2</sup></td><td>3.517</td><td>2.42 × 10<sup>-2</sup></td><td>0.028</td><td>3.80 × 10<sup>-1</sup></td><td>1.553</td><td>3.13 × 10<sup>-11</sup></td><td>0.017</td><td>7.28 × 10<sup>-1</sup></td><td>0.211</td><td>5.37 × 10<sup>-2</sup></td><td>2</td><td>44,333,455</td></tr><tr><td>Igfbp6</td><td>1.110</td><td>2.04 × 10<sup>-4</sup></td><td>0.081</td><td>6.67 × 10<sup>-1</sup></td><td>2.796</td><td>2.01 × 10<sup>-1</sup></td><td>0.821</td><td>1.64 × 10<sup>-2</sup></td><td>2.720</td><td>3.41 × 10<sup>-10</sup></td><td>-0.079</td><td>4.68 × 10<sup>-1</sup></td><td>0.325</td><td>6.21 × 10<sup>-3</sup></td><td>3</td><td>43,445,554</td></tr><tr><td>Adk</td><td>-0.161</td><td>5.78 × 10<sup>-2</sup></td><td>0.338</td><td>8.46 × 10<sup>-2</sup></td><td>1.000</td><td>4.86 × 10<sup>-2</sup></td><td>0.117</td><td>1.76 × 10<sup>-1</sup></td><td>2.226</td><td>2.74 × 10<sup>-8</sup></td><td>0.053</td><td>6.57 × 10<sup>-1</sup></td><td>0.261</td><td>4.21 × 10<sup>-2</sup></td><td>0</td><td>43,345,555</td></tr><tr><td>Actr10</td><td>-0.120</td><td>1.99 × 10<sup>-2</sup></td><td>0.554</td><td>4.28 × 10<sup>-3</sup></td><td>0.816</td><td>7.05 × 10<sup>-2</sup></td><td>0.170</td><td>2.88 × 10<sup>-3</sup></td><td>-0.207</td><td>4.49 × 10<sup>-3</sup></td><td>0.129</td><td>2.78 × 10<sup>-1</sup></td><td>0.062</td><td>5.24 × 10<sup>-1</sup></td><td>0</td><td>43,332,233</td></tr></tbody></table>

$|\log_2 \text{FC}| > 0.263$ and $p$-values $\le 0.05$ are in bold. Consensus score indicates whether FLX induces mainly upregulation (positive score) or downregulation (negative score) of gene expression. Genes with absolute value of consensus score equal to 5 are in bold.
<table><caption>Table 7. *Cont.*</caption><thead><tr><th rowspan="2">Dataset ID<br>Gene</th><th colspan="2">GSE118669</th><th colspan="2">GSE42940</th><th colspan="2">GSE4307</th><th colspan="2">GSE6476</th><th colspan="2">GSE84185</th><th colspan="2">SRP056481</th><th colspan="2">SRP057486</th><th rowspan="2">Consensus Score</th><th rowspan="2">Portrait Score</th></tr><tr><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th><th>log<sub>2</sub>FC</th><th>p</th></tr></thead><tbody><tr><td>Drd1</td><td>1.070</td><td>3.46 × 10<sup>-3</sup></td><td>-0.279</td><td>5.07 × 10<sup>-1</sup></td><td>0.000</td><td>1.00</td><td>1.734</td><td>1.68 × 10<sup>-3</sup></td><td>1.472</td><td>1.73 × 10<sup>-8</sup></td><td>-0.077</td><td>4.06 × 10<sup>-1</sup></td><td>0.251</td><td>3.38 × 10<sup>-2</sup></td><td>4</td><td>43,332,223</td></tr><tr><td>Tyrp3</td><td>0.267</td><td>4.15 × 10<sup>-2</sup></td><td>-0.202</td><td>4.39 × 10<sup>-1</sup></td><td>1.572</td><td>2.11 × 10<sup>-2</sup></td><td>-0.438</td><td>4.39 × 10<sup>-2</sup></td><td>1.818</td><td>1.16 × 10<sup>-9</sup></td><td>0.128</td><td>2.77 × 10<sup>-1</sup></td><td>-0.131</td><td>3.40 × 10<sup>-1</sup></td><td>1</td><td>43,321,112</td></tr><tr><td>Sox11</td><td>0.614</td><td>6.75 × 10<sup>-3</sup></td><td>1.267</td><td>1.63 × 10<sup>-3</sup></td><td>1.617</td><td>5.52 × 10<sup>-2</sup></td><td>1.215</td><td>9.06 × 10<sup>-6</sup></td><td>-0.476</td><td>1.15 × 10<sup>-4</sup></td><td>-0.033</td><td>7.03 × 10<sup>-1</sup></td><td>-0.205</td><td>1.32 × 10<sup>-1</sup></td><td>2</td><td>43,321,111</td></tr><tr><td>Trpm3</td><td>-0.437</td><td>1.49 × 10<sup>-4</sup></td><td>0.107</td><td>7.64 × 10<sup>-1</sup></td><td>0.000</td><td>1.00</td><td>-0.946</td><td>1.39 × 10<sup>-2</sup></td><td>-1.568</td><td>9.87 × 10<sup>-9</sup></td><td>-0.188</td><td>1.06 × 10<sup>-1</sup></td><td>-0.016</td><td>8.63 × 10<sup>-1</sup></td><td>-3</td><td>-44,443,334</td></tr><tr><td>Jun</td><td>-0.439</td><td>1.30 × 10<sup>-2</sup></td><td>-0.609</td><td>3.06 × 10<sup>-1</sup></td><td>-1.000</td><td>1.07 × 10<sup>-1</sup></td><td>0.267</td><td>8.55 × 10<sup>-2</sup></td><td>-2.150</td><td>7.95 × 10<sup>-11</sup></td><td>-0.143</td><td>2.20 × 10<sup>-1</sup></td><td>-0.003</td><td>9.71 × 10<sup>-1</sup></td><td>-2</td><td>-44,444,444</td></tr><tr><td>Efnb3</td><td>-0.383</td><td>1.12 × 10<sup>-2</sup></td><td>0.000</td><td>1.00</td><td>-0.694</td><td>2.27 × 10<sup>-1</sup></td><td>-0.686</td><td>1.56 × 10<sup>-3</sup></td><td>-2.080</td><td>2.00 × 10<sup>-12</sup></td><td>-0.128</td><td>2.74 × 10<sup>-1</sup></td><td>0.197</td><td>8.16 × 10<sup>-2</sup></td><td>-3</td><td>-44,444,444</td></tr><tr><td>Lct</td><td>-0.600</td><td>1.75 × 10<sup>-4</sup></td><td>-0.567</td><td>1.15 × 10<sup>-1</sup></td><td>-1.100</td><td>1.16 × 10<sup>-1</sup></td><td>0.266</td><td>5.55 × 10<sup>-1</sup></td><td>-4.916</td><td>4.72 × 10<sup>-12</sup></td><td>0.011</td><td>3.89 × 10<sup>-1</sup></td><td>-0.421</td><td>3.72 × 10<sup>-4</sup></td><td>-3</td><td>-44,444,444</td></tr><tr><td>Pdla6</td><td>-0.390</td><td>1.76 × 10<sup>-4</sup></td><td>0.440</td><td>2.97 × 10<sup>-1</sup></td><td>-1.320</td><td>1.05 × 10<sup>-1</sup></td><td>-0.247</td><td>7.76 × 10<sup>-4</sup></td><td>-0.112</td><td>2.60 × 10<sup>-1</sup></td><td>-0.070</td><td>5.61 × 10<sup>-1</sup></td><td>-0.382</td><td>2.62 × 10<sup>-4</sup></td><td>-3</td><td>-44,444,444</td></tr><tr><td>Kcniq3</td><td>-0.297</td><td>9.57 × 10<sup>-3</sup></td><td>0.000</td><td>1.00</td><td>0.000</td><td>1.00</td><td>-0.452</td><td>1.47 × 10<sup>-2</sup></td><td>-0.176</td><td>6.74 × 10<sup>-2</sup></td><td>-0.176</td><td>1.33 × 10<sup>-1</sup></td><td>-0.367</td><td>3.62 × 10<sup>-3</sup></td><td>-3</td><td>-44,444,455</td></tr><tr><td>Mcm6</td><td>-0.593</td><td>3.24 × 10<sup>-5</sup></td><td>-0.526</td><td>5.41 × 10<sup>-2</sup></td><td>-1.105</td><td>7.72 × 10<sup>-2</sup></td><td>-0.419</td><td>4.23 × 10<sup>-1</sup></td><td>-2.518</td><td>1.14 × 10<sup>-10</sup></td><td>-0.028</td><td>7.34 × 10<sup>-1</sup></td><td>-0.008</td><td>9.38 × 10<sup>-1</sup></td><td>-2</td><td>-44,445,666</td></tr><tr><td>Cacna1d</td><td>-0.153</td><td>1.65 × 10<sup>-1</sup></td><td>-0.613</td><td>1.86 × 10<sup>-2</sup></td><td>-0.830</td><td>9.70 × 10<sup>-2</sup></td><td>-0.309</td><td>1.29 × 10<sup>-2</sup></td><td>-0.159</td><td>6.13 × 10<sup>-2</sup></td><td>-0.082</td><td>4.98 × 10<sup>-1</sup></td><td>-0.240</td><td>4.85 × 10<sup>-2</sup></td><td>-3</td><td>-44,566,677</td></tr><tr><td>Nfia</td><td>-0.282</td><td>1.95 × 10<sup>-3</sup></td><td>0.077</td><td>6.30 × 10<sup>-1</sup></td><td>-0.408</td><td>4.77 × 10<sup>-1</sup></td><td>-0.373</td><td>6.99 × 10<sup>-3</sup></td><td>-1.022</td><td>2.60 × 10<sup>-9</sup></td><td>-0.075</td><td>5.19 × 10<sup>-1</sup></td><td>-0.369</td><td>5.85 × 10<sup>-3</sup></td><td>-4</td><td>-44,655,556</td></tr><tr><td>Nedd4l</td><td>-0.276</td><td>4.08 × 10<sup>-3</sup></td><td>-0.653</td><td>1.53 × 10<sup>-2</sup></td><td>-0.283</td><td>3.90 × 10<sup>-1</sup></td><td>-0.377</td><td>2.13 × 10<sup>-3</sup></td><td>-0.326</td><td>4.03 × 10<sup>-4</sup></td><td>-0.079</td><td>5.10 × 10<sup>-1</sup></td><td>-0.258</td><td>2.86 × 10<sup>-2</sup></td><td>-5</td><td>-44,677,777</td></tr><tr><td>Fndc1</td><td>-0.377</td><td>3.74 × 10<sup>-3</sup></td><td>0.000</td><td>1.00</td><td>-1.356</td><td>1.14 × 10<sup>-1</sup></td><td>0.361</td><td>1.25 × 10<sup>-1</sup></td><td>-1.817</td><td>4.60 × 10<sup>-9</sup></td><td>-0.113</td><td>1.43 × 10<sup>-1</sup></td><td>-0.233</td><td>8.75 × 10<sup>-2</sup></td><td>-2</td><td>-45,444,444</td></tr><tr><td>Ntf3</td><td>-0.938</td><td>2.00 × 10<sup>-3</sup></td><td>0.208</td><td>3.41 × 10<sup>-1</sup></td><td>-3.919</td><td>2.35 × 10<sup>-1</sup></td><td>-1.793</td><td>5.70 × 10<sup>-4</sup></td><td>-3.831</td><td>1.97 × 10<sup>-12</sup></td><td>-0.156</td><td>1.64 × 10<sup>-1</sup></td><td>0.015</td><td>8.04 × 10<sup>-1</sup></td><td>-3</td><td>-45,444,444</td></tr><tr><td>Itga4</td><td>-0.437</td><td>1.33 × 10<sup>-2</sup></td><td>-0.007</td><td>9.83 × 10<sup>-1</sup></td><td>-3.476</td><td>6.77 × 10<sup>-2</sup></td><td>-0.601</td><td>2.59 × 10<sup>-2</sup></td><td>-1.608</td><td>1.14 × 10<sup>-10</sup></td><td>-0.012</td><td>8.44 × 10<sup>-1</sup></td><td>-0.490</td><td>2.28 × 10<sup>-5</sup></td><td>-4</td><td>-45,555,466</td></tr><tr><td>Kbtbd11</td><td>-0.191</td><td>1.58 × 10<sup>-2</sup></td><td>0.000</td><td>1.00</td><td>0.000</td><td>1.00</td><td>-0.998</td><td>1.18 × 10<sup>-2</sup></td><td>-1.353</td><td>2.14 × 10<sup>-10</sup></td><td>-0.158</td><td>1.82 × 10<sup>-1</sup></td><td>-0.231</td><td>7.09 × 10<sup>-2</sup></td><td>-3</td><td>-45,555,555</td></tr><tr><td>Kcnd4</td><td>-0.238</td><td>2.69 × 10<sup>-3</sup></td><td>0.000</td><td>1.00</td><td>-0.544</td><td>3.16 × 10<sup>-1</sup></td><td>-0.516</td><td>2.34 × 10<sup>-3</sup></td><td>-1.775</td><td>2.59 × 10<sup>-10</sup></td><td>-0.223</td><td>6.47 × 10<sup>-2</sup></td><td>0.021</td><td>7.70 × 10<sup>-1</sup></td><td>-3</td><td>-45,555,555</td></tr><tr><td>Tuxb</td><td>-0.323</td><td>4.79 × 10<sup>-4</sup></td><td>0.000</td><td>1.00</td><td>-0.695</td><td>3.40 × 10<sup>-1</sup></td><td>-0.574</td><td>7.37 × 10<sup>-3</sup></td><td>-1.462</td><td>6.94 × 10<sup>-9</sup></td><td>0.000</td><td>1.00</td><td>-0.398</td><td>3.57 × 10<sup>-3</sup></td><td>-4</td><td>-45,555,555</td></tr><tr><td>Scn3b</td><td>-0.245</td><td>7.02 × 10<sup>-3</sup></td><td>-0.479</td><td>1.05 × 10<sup>-2</sup></td><td>-0.678</td><td>1.15 × 10<sup>-1</sup></td><td>-0.594</td><td>5.20 × 10<sup>-5</sup></td><td>-0.971</td><td>4.31 × 10<sup>-8</sup></td><td>0.178</td><td>1.41 × 10<sup>-1</sup></td><td>-0.304</td><td>2.13 × 10<sup>-2</sup></td><td>-5</td><td>-45,555,555</td></tr><tr><td>Foxo1</td><td>-0.500</td><td>2.14 × 10<sup>-3</sup></td><td>-0.208</td><td>3.19 × 10<sup>-1</sup></td><td>-0.840</td><td>1.99 × 10<sup>-1</sup></td><td>-0.587</td><td>1.09 × 10<sup>-3</sup></td><td>-1.836</td><td>1.47 × 10<sup>-12</sup></td><td>-0.137</td><td>2.26 × 10<sup>-1</sup></td><td>-0.037</td><td>7.27 × 10<sup>-1</sup></td><td>-3</td><td>-45,566,666</td></tr><tr><td>Rasgrf1</td><td>-0.394</td><td>5.36 × 10<sup>-5</sup></td><td>-0.200</td><td>2.02 × 10<sup>-1</sup></td><td>-0.443</td><td>3.55 × 10<sup>-1</sup></td><td>-0.519</td><td>4.75 × 10<sup>-4</sup></td><td>-1.051</td><td>2.38 × 10<sup>-10</sup></td><td>-0.006</td><td>9.59 × 10<sup>-1</sup></td><td>-0.335</td><td>5.16 × 10<sup>-4</sup></td><td>-4</td><td>-45,666,667</td></tr><tr><td>Auts2</td><td>-0.223</td><td>1.52 × 10<sup>-3</sup></td><td>0.000</td><td>1.00</td><td>-0.491</td><td>3.48 × 10<sup>-1</sup></td><td>-0.489</td><td>5.37 × 10<sup>-3</sup></td><td>-1.047</td><td>2.47 × 10<sup>-8</sup></td><td>-0.117</td><td>3.18 × 10<sup>-1</sup></td><td>-0.245</td><td>7.51 × 10<sup>-2</sup></td><td>-3</td><td>-46,666,666</td></tr><tr><td>Doc2b</td><td>-0.506</td><td>1.51 × 10<sup>-3</sup></td><td>0.371</td><td>1.83 × 10<sup>-1</sup></td><td>-0.798</td><td>7.63 × 10<sup>-2</sup></td><td>-0.233</td><td>1.07 × 10<sup>-2</sup></td><td>-3.232</td><td>3.42 × 10<sup>-11</sup></td><td>0.028</td><td>7.90 × 10<sup>-1</sup></td><td>-0.416</td><td>4.12 × 10<sup>-4</sup></td><td>-4</td><td>-54,444,433</td></tr><tr><td>Dsp</td><td>-0.749</td><td>2.12 × 10<sup>-6</sup></td><td>0.198</td><td>4.54 × 10<sup>-1</sup></td><td>-1.479</td><td>9.31 × 10<sup>-2</sup></td><td>-1.688</td><td>9.19 × 10<sup>-3</sup></td><td>-5.642</td><td>5.44 × 10<sup>-14</sup></td><td>-0.125</td><td>1.86 × 10<sup>-1</sup></td><td>-0.023</td><td>7.13 × 10<sup>-1</sup></td><td>-3</td><td>-55,544,444</td></tr><tr><td>Slc4a4</td><td>-0.537</td><td>5.83 × 10<sup>-4</sup></td><td>0.052</td><td>8.62 × 10<sup>-1</sup></td><td>-1.062</td><td>1.43 × 10<sup>-1</sup></td><td>-0.912</td><td>1.93 × 10<sup>-3</sup></td><td>-1.025</td><td>2.36 × 10<sup>-12</sup></td><td>-0.166</td><td>1.61 × 10<sup>-1</sup></td><td>0.043</td><td>7.39 × 10<sup>-1</sup></td><td>-3</td><td>-55,554,455</td></tr><tr><td>Pcdh19</td><td>-0.572</td><td>4.43 × 10<sup>-3</sup></td><td>0.000</td><td>1.00</td><td>0.000</td><td>1.00</td><td>-0.252</td><td>1.31 × 10<sup>-2</sup></td><td>-1.419</td><td>9.69 × 10<sup>-15</sup></td><td>-0.114</td><td>3.35 × 10<sup>-1</sup></td><td>-0.413</td><td>1.08 × 10<sup>-4</sup></td><td>-4</td><td>-55,555,555</td></tr><tr><td>Fat4</td><td>-0.890</td><td>4.54 × 10<sup>-4</sup></td><td>0.000</td><td>1.00</td><td>0.000</td><td>1.00</td><td>-1.352</td><td>1.80 × 10<sup>-3</sup></td><td>-1.748</td><td>2.48 × 10<sup>-11</sup></td><td>-0.210</td><td>6.03 × 10<sup>-2</sup></td><td>-0.342</td><td>1.25 × 10<sup>-2</sup></td><td>-4</td><td>-55,555,555</td></tr><tr><td>Akt3</td><td>-0.224</td><td>5.02 × 10<sup>-3</sup></td><td>-0.482</td><td>1.34 × 10<sup>-2</sup></td><td>-0.819</td><td>1.04 × 10<sup>-1</sup></td><td>-0.372</td><td>3.84 × 10<sup>-3</sup></td><td>-0.602</td><td>8.66 × 10<sup>-9</sup></td><td>0.081</td><td>4.84 × 10<sup>-1</sup></td><td>-0.128</td><td>1.81 × 10<sup>-1</sup></td><td>-4</td><td>-55,555,555</td></tr><tr><td>Sipa1l2</td><td>-0.553</td><td>2.89 × 10<sup>-5</sup></td><td>0.000</td><td>1.00</td><td>-0.804</td><td>1.28 × 10<sup>-1</sup></td><td>-0.469</td><td>6.07 × 10<sup>-3</sup></td><td>-0.918</td><td>1.36 × 10<sup>-8</sup></td><td>-0.200</td><td>9.62 × 10<sup>-2</sup></td><td>-0.050</td><td>6.56 × 10<sup>-1</sup></td><td>-3</td><td>-55,555,556</td></tr><tr><td>Zfp316</td><td>-0.218</td><td>4.28 × 10<sup>-3</sup></td><td>-1.324</td><td>3.00 × 10<sup>-4</sup></td><td>-1.531</td><td>1.80 × 10<sup>-2</sup></td><td>-0.912</td><td>8.63 × 10<sup>-3</sup></td><td>-0.448</td><td>1.30 × 10<sup>-10</sup></td><td>-0.107</td><td>3.35 × 10<sup>-1</sup></td><td>0.019</td><td>8.72 × 10<sup>-1</sup></td><td>-5</td><td>-55,556,666</td></tr><tr><td>Sllt1</td><td>-0.284</td><td>1.93 × 10<sup>-4</sup></td><td>-0.198</td><td>4.17 × 10<sup>-1</sup></td><td>-0.807</td><td>1.18 × 10<sup>-1</sup></td><td>0.051</td><td>4.40 × 10<sup>-1</sup></td><td>-2.120</td><td>6.68 × 10<sup>-12</sup></td><td>-0.119</td><td>3.23 × 10<sup>-1</sup></td><td>-0.407</td><td>2.67 × 10<sup>-4</sup></td><td>-3</td><td>-55,566,665</td></tr><tr><td>Itsn1</td><td>-0.238</td><td>3.80 × 10<sup>-3</sup></td><td>-2.471</td><td>1.91 × 10<sup>-3</sup></td><td>-0.915</td><td>5.91 × 10<sup>-2</sup></td><td>-0.232</td><td>2.93 × 10<sup>-3</sup></td><td>-0.785</td><td>1.52 × 10<sup>-8</sup></td><td>-0.068</td><td>5.71 × 10<sup>-1</sup></td><td>-0.164</td><td>1.84 × 10<sup>-1</sup></td><td>-4</td><td>-55,777,777</td></tr></tbody></table>

$|\log_2 \text{FC}| > 0.263$ and $p$-values $\le 0.05$ are in bold. Consensus score indicates whether FLX induces mainly upregulation (positive score) or downregulation (negative score) of gene expression. Genes in bold are also in Table 6.
<table><caption>**Table 8.** Ontological analysis of downregulated genes (consensus and portrait scores< −2) after FLX treatment in naive animals.</caption><thead><tr><th>Category</th><th>ID</th><th>Name</th><th>Genes</th><th>Count</th><th>Fold Enrichment</th><th>Padj</th></tr></thead><tbody><tr><td rowspan="2">Cellular component</td><td>GO:0030054</td><td>cell junction</td><td>Actn2, Ahi1, Arhgap21, Cacna1c, Camk2b, Ccd85a, Chrm1, Chrna7, Clstn2, Ctbl2, Ddn, Dlg5, Dsp, Enah, Eps8, Faim2, Flrt3, Frmd4b, Gabra5, Gabrd, Gria1, Grid1, Hap1, Itsn1, Kctd12, Lgi1, Limd1, Lrrc7, Lrrtm1, Map1b, Mpdz, Mpp2, Nbea, Ncdn, Nfia, Nlgn3, Nos1, Nphlp1, Nrxn1, Palm, Palmd, Pcdh1, Plpr4, Prkcg, Prrt1, Ptk2b, Rapgef2, Rgs12, Sarm1, Scn8a, Sema4f, Shisa6, Slc1a1, Sorbs1, Sptb, Sptbn4, Syt7, Tenm2, Traf4, Tspan33, Wasf1</td><td>61</td><td>2.43</td><td>1.70 × 10<sup>-7</sup></td></tr><tr><td>GO:0045202</td><td>synapse</td><td>Ablim3, Bcan, Cacna1c, Calb1, Camk2b, Capn5, Chrm1, Chrna7, Clstn2, Ctbl2, Ddn, Dlg5, Doc2b, Enah, Eps8, Faim2, Gabra5, Gabrd, Gria1, Grid1, Grip1, Hap1, Itsn1, Kctd12, Lgi1, Lrrc7, Lrrtm1, Lypd1, Map1a, Map1b, Mpdz, Mpp2, Nbea, Ncdn, Nlgn3, Nos1, Nrxn1, Palm, Palmd, Plpr4, Ppfia4, Prkcg, Prrt1, Rapgef2, Rgs12, Sarm1, Sema4f, Shisa6, Slc1a1, Sorbs1, Syt7, Tenm2, Wasf1</td><td>53</td><td>2.23</td><td>5.62 × 10<sup>-5</sup></td></tr><tr><td>Biological process</td><td>GO:2000310</td><td>regulation of NMDA receptor activity</td><td>Dapk1, Mapk8ip2, Nlgn3, Nrxn1, Ptk2b, Rasgrf1, Rasgrf2</td><td>7</td><td>13.9</td><td>1.83 × 10<sup>-2</sup></td></tr><tr><td rowspan="2">KEGG pathway</td><td>mmu04713</td><td>circadian entrainment</td><td>Adcy1, Adcy9, Adcyap1r1, Cacna1c, Cacna1d, Camk2b, Gng7, Gria1, Nos1, Per3, Prkcg, Ryr1, Ryr2</td><td>13</td><td>4.82</td><td>3.85 × 10<sup>-3</sup></td></tr><tr><td>mmu04911</td><td>insulin secretion</td><td>Abcc8, Adcy1, Adcy9, Adcyap1r1, Atf6b, Cacna1c, Cacna1d, Camk2b, Gck, Prkcg, Ryr2</td><td>11</td><td>4.65</td><td>3.12 × 10<sup>-2</sup></td></tr></tbody></table>

Genes belonging to the core signatures (Tables 6 and 7) are in bold.
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Venn diagram of hippocampal transcriptional signatures. The signatures were induced by either FLX in a stress-inducing depressive-like behavior paradigm (S-FLX vs. S-Veh) or FLX in naïve rodents (FLX vs. Veh) and the list of genes were obtained by either integration (significant in at least four experiments) or scoring methods (both consensus and portrait). Genes in red are generally overexpressed, genes in blue are downregulated whereas genes in black show no consensual trend of variation direction.</figcaption></figure>

<table><caption>**Table 9.** Ontological analysis of upregulated genes (consensus and portrait scores ≥4) after FLX treatment.</caption><thead><tr><th>Category</th><th>ID</th><th>Name</th><th>Genes</th><th>Count</th><th>Fold Enrichment</th><th>Padj</th></tr></thead><tbody><tr><td rowspan="3">Cellular component</td><td>GO:0030424</td><td>axon</td><td>Bdnf, Cdk5r1, Cst3, Dclk1, Dcx, Drd1, Ntrk2, Pink1, Ptpn5, Sema3a, Stx1b, Stxbp1, Tubb3</td><td>13</td><td>4.97</td><td>2.26 × 10<sup>-3</sup></td></tr><tr><td>GO:0009986</td><td>cell surface</td><td>Bgn, Cd14, Cd151, Drd1, Ephb6, Fcgr4, H2-K1, H2-Q10, Ifitm3, Itgav, Ntrk2, Prlr, Tfrc, Tyro3</td><td>14</td><td>3.48</td><td>3.44 × 10<sup>-2</sup></td></tr><tr><td>GO:0005576</td><td>extracellular region</td><td>Bdnf, Bgn, Brinp2, C1qb, Ccdc3, Cd14, Ccn4, Cst3, Ephb6, Gsn, Igfbp6, Lyzl4, Nptx2, Olfm3, Pcsk2, Sema3a, Sema3c, Serpina3n, Spock3, Tfrc, Vgf, Vip, Wnt10a</td><td>23</td><td>2.37</td><td>3.62 × 10<sup>-2</sup></td></tr><tr><td>KEGG pathway</td><td>mmu04010</td><td>MAPK signaling pathway</td><td>Bdnf, Cacna2d2, Cd14, Gadd45a, Kras, Mapk9, Mapk11, Ntrk2, Ptpn5, Rasgrp1</td><td>10</td><td>5.65</td><td>9.09 × 10<sup>-3</sup></td></tr><tr><td>Molecular function</td><td>GO:0016301</td><td>kinase activity</td><td>Cdk5r1, Cdkl5, Cmpk1, Dclk1, Mapk4, Mapk9, Mapk11, Melk, Ntrk2, Pak5, Pink1, Prkar1b, Tyro3</td><td>13</td><td>3.73</td><td>4.60 × 10<sup>-2</sup></td></tr></tbody></table>

Genes belonging to the core signature (Tables S11 and S12: genes with absolute scores ≥6, Figure 3) are in bold.
<table><caption>**Table 10.** Ontological analysis of downregulated genes (consensus and portrait scores ≤ −4) after FLX treatment.</caption><thead><tr><th>Category</th><th>ID</th><th>Name</th><th>Genes</th><th>Count</th><th>Fold Enrichment</th><th>Padj</th></tr></thead><tbody><tr><td>KEGG pathway</td><td>mmu04010</td><td>MAPK signaling pathway</td><td>Akt3, Cacna1d, Flt3, Hspa2, Il1r1, Map3k4, Ntf3, Rasgrf1, Rps6kas, Stk4, Tgfb2</td><td>11</td><td>5.17</td><td>7.67 × 10<sup>-3</sup></td></tr></tbody></table>

Genes belonging to the core signature (Table S11 and S12: genes with absolute scores ≥6, Figure 3) are in bold.

### 3. Discussion

In this study we have voluntarily limited our field of research to post-2006 publications on transcriptional signatures of the effect of AD treatments in animal models and have identified 114 studies. Nevertheless, we were confronted with the difficulty of accessing complete transcriptional data for more than 60 studies for which data were not available in a public manner. This unfortunately reflects, at the minimum, a lack of reflexes on the part of investigators, or perhaps also a lack of sufficient incentive on the part of the authors and publishers, to share data [59]. Indeed, omics studies come at a substantial cost, and it is conceivable that some authors would prefer to fully capitalize on future analyses before sharing the data publicly.

#### 3.1. Methodological Constraints

##### 3.1.1. Current Picture: A Fragmented Landscape

A detailed review of the available studies shows a fragmented landscape regarding the design of the studies and an important imbalance in many experimental factors regarding sex, region, and/or ADs studied (Figure 1).

Although women are twice as affected by major depression, preclinical study models of this disease and its treatment are still too often focused on male animals. Thus, among all the studies studying the transcriptome following the action of AD molecules, less than 10% concern females (Figure 1). In addition, recent studies suggest that the transcriptional signatures of depression are different in men and women [60]. It is therefore essential that more preclinical studies seriously address the changes induced in females by the major therapeutic agents that are or will be offered to patients with mood disorders.

A similar imbalance exists for brain regions targeted in transcriptomic studies of AD effects (Figure 1). This is an important caveat because, while stress (acute and prolonged) causes molecular and functional changes in several brain areas such as the medial prefrontal cortex, amygdala, nucleus accumbens, and hippocampus, among the most studied, the stress-induced molecular signatures differ across regions. For instance, activation of the CREB-BDNF axis in the hippocampus is antidepressant, but the same activation in the nucleus accumbens is pro-depressant [61]. Region-dependent effects of stress on neuronal activation and BDNF signaling could also be sex-dependent [62,63]. Comparing AD transcriptomic signatures across regions in male and female mice could provide potential insights into how brain pharmacology can be modulated in a region-dependent and sex-specific manner; however, we could not undertake this analysis because of insufficient power.

Another big discrepancy that we observed across studies concerns pharmacological classes tested and compounds inside the same class. Although SSRIs remain the most prescribed ADs, it is important to revisit older, less selective classes of ADs. Indeed, many open questions remain in terms of both the delay of onset [64] and of the mechanism of action [65]. Again, there was not enough power (numbers) in the available studies to consider for our analysis. Therefore, our review of the literature underlines the need for systematic studies designed to include male and female mice, different regions, and more than one class of ADs.

##### 3.1.2. Methodological Choices

It is not a surprise that the majority of the available studies concerns the most prescribed AD in the world, Prozac, approved by the Food and Drug Administration (FDA)
in 1987, whose active ingredient is FLX [66]. Moreover, the preferred study area for the transcriptional effect induced by ADs is the hippocampus, a region consistently implicated in major depression and AD action [67]. Thus, we decided to focus on the transcriptomic signature produced by FLX in the hippocampus. In the different studies analyzed for some is the whole hippocampus that has been collected; in other cases, it is the dentate gyrus and/or the cornu ammonis 1 (CA1) and sometimes it is specified whether it is the ventral or dorsal part (Table 1). The hippocampus is functionally divided into a dorsal region that is primarily engaged in cognitive functions and a ventral region that regulates emotion [68,69]. Although it has been observed for one of the experiments we included that similar gene expression responses to FLX were found in dorsal and ventral dentate gyrus [51], more recent investigations show that the ventral hippocampus was particularly sensitive to the effects of stress. Therefore, it was proposed to consider the dorsal and ventral hippocampus separately when conducting high-throughput molecular analyses [24,70]. We are aware of this caveat, but we had no choice if we wanted to keep enough power in our analysis. Stephen Gammie has recently drawn a picture of depression based on publicly available transcriptomic data from human brains [23]. His comparison of the human signature with the signatures of nearly 200 data sets, mainly from mice or rats treated with ADs and other drug classes, showed that the effect of FLX assessed at the hippocampus level provided the closest inverted signature to that of human depression [23].

Another choice we made, for the sake of homogeneity, was to select, at least for one experiment, samples from mice responding to FLX (GSE84185). Indeed, if it seems that most investigators retained for transcriptomic analysis only samples from animals showing an improvement in behavior following FLX administration, in at least one case the overall response was profiled [21]. It may be noted that at least two studies had highlighted interindividual variability in behavioral responses to FLX [21,51]. It was even proposed that it would be interesting to evaluate ambiguous responders, as they could be useful in dissociating anti-depression-like effects from anti-anxiety-like effects [51]. Later, such behavioral heterogeneity, especially concerning anxiety, following chronic FLX treatment was confirmed and may reflect a specific gene expression profile [46].

A limitation of our work is that we almost only processed protein-coding mRNAs. However, emerging literature, especially in terms of biomarker identification for mood disorders, concerns microRNAs (miRNAs) [71,72]. Thus, among the studies that we identified as potentially relevant for our analysis, nine concerned miRNAs, two of which were interested in the effect of FLX in the hippocampus [38,73]. It will therefore be important in the future to compare miRNA and mRNA signatures to study whether they are related [74,75], which will help to better understand the mechanism of action of ADs and undoubtedly improve their effectiveness [76]. In addition to miRNAs, a large family of non-coding RNAs could not be considered in our analysis. These are the myriad of long non-coding RNAs (lncRNAs) [77]. However, not only are they very present for a while on microarrays, and a fortiori are widely detected by high-throughput sequencing, but it is often in this class of RNA that we find the strongest expression variations [78,79]. Efforts to standardize annotations and a better knowledge of their exact role should in the future help us to better integrate them into a more global signature of the effect of ADs [80,81]. Circular RNAs are another category that is nowadays emerging in the field of psychiatry transcriptomics [82]. However, the available data in animal models of depression and antidepressant action are still limited, albeit promising [83].

## 3.2. Biological Pathways and Genes

### 3.2.1. Serotonin System and BDNF

The primary target of FLX is the serotonin transporter 5-HTT, encoded by the *SLC6A4* gene. Although there may be transient effects on variations in Slc6a4 mRNA expression levels in rodent models [84,85], there is little evidence to suggest a transcriptional effect of FLX on either *SLC6A4* [86] or the limiting enzyme for serotonin synthesis in brain, tryptophan hydroxylase 2 (TPH2) [28], and we did not observe robust modifications for
these two genes. On the other hand, in the large family of serotonin receptors, several studies have implicated transcriptional variations in the gene encoding the 1B receptor (HTR1B). We confirm in our meta-analysis this effect for *Htr1b* but also for *Htr5b* (Table 6), for which there is no functional equivalent in humans [87]. One target related to serotonin is the overexpressed *Vip* messenger RNA (Tables 7, S8 and S12) encoding vasoactive intestinal peptide, as it has been known for a long time that the neurotransmitter effect of serotonin in the hippocampus can be modulated by *Vip* through regulation of cyclic adenosine 3':5'-monophosphate (cAMP) levels and serotonin receptors [88–90]. It was also found that FLX, while improving depressed behavior in a rat model of chronic stress-induced depression, increased *Vip* expression [91]. Other expected targets for FLX include brain-derived neurotrophic factor (BDNF), which is widely implicated in depression and the mechanisms of its treatment and would be expected to be increased [28], which is confirmed both in a stress paradigm (Tables 4 and S5), or in naive animals (Table S8) and thus also in the overall meta-analysis (Tables 9 and S12).

### 3.2.2. Immediate Early Genes

In the present meta-analysis, among the genes significantly modulated by FLX appear the immediate early genes *Nr4a1*, *Nr4a2*, *Arc*, *Egr1*, *Fosb*, *Fos12*, and *Junb* (Tables 2, S4 and S6). These are genes that have long been identified as particularly important in the program of gene expression modification during the induction of depression in different animal models [92–94]. EGR1, in particular, is the major downstream partner of the ERK/Elk-1 cascade that we have recently proposed as a novel target for AD development [67]. Interestingly, such immediate early genes are regulated not only in the hippocampus but also in other brain areas such as frontal and prefrontal cortex, lateral amygdala, and nucleus accumbens [95–99], and they are the indirect target of other AD drugs than FLX such as agomelatine, duloxetine, imipramine, ketamine, paroxetine, or vortioxetine [100–105] but also of other AD-related interventions such as sleep deprivation and deep brain stimulation [106–108]. In addition to this set of genes, we also noted the immediate early response 5 gene (*Ier5*), which appears upregulated in our meta-analysis (Tables 3, S5, S8 and S12), was previously found dysregulated in peripheral blood mononuclear cells of unmedicated mood-disorder patients compared to healthy controls [109].

### 3.2.3. Signal Transduction Pathways

It is not a surprise that our analysis highlights universal signal transduction pathways. Following receptor activation, second messengers fine-tune neuronal activity, synaptic remodeling, long-term potentiation, and neurogenesis. They thus regulate synaptic plasticity and cellular resilience, processes by which the brain perceives, adapts, and responds to a variety of internal and external stimuli (including stress and other depressogenic factors [110]). Over the last decades, evidence from different research groups showed that kinase-phosphatase pathways are important mediators of AD action, and have used translational models, namely genetically modified animals to phenocopy AD-sensitive behaviors and to highlight AD-sensitive brain circuits [111–114].

#### a MAP kinases

In line with the neurotrophic hypothesis of depression, genes coding for members of the neurotrophic tyrosine receptor kinase family such as *Ntrk2* and *Ntrk3* were normalized by FLX action (Tables 4, 6, 9 and S12) in a reverse direction compared to observation in rodent stress models of depression [115]. These genes encode kinases that, upon neurotrophin binding, phosphorylate members of the MAPK pathway. Thus, it makes sense to find the MAPK signaling pathway as one of the pathways most consistently affected by FLX (Tables 9 and 10). Among the multiple connections between BDNF, MAP kinases, and AD action, we found an upregulation of the *Ptpn5* gene expression (Tables S5, S8, S12 and 9 and Figure 3). It encodes a striatal-enriched tyrosine protein phosphatase (STEP) that inactivates key neuronal signaling proteins such as MAP kinases, tyrosine kinases NMDA, and AMPA
receptors and whose inactivation decreases the expression of BDNF [116]. Therefore, STEP inhibitors have been proposed as a novel target for AD drugs of the new generation [117].

b WNT/catenin

Another potential player in this activation of the MAPK pathway is insulin-like growth factor binding protein 6 encoded by *Igfbp6* mRNA, which is overexpressed in the different paradigms in response to FLX (Tables 6, 7 and 9 and Figure 3) including studies not included in the current meta-analysis [118]. There is even another member of this gene family with *Igfbp4* (Figure 3), which can be linked to the activation of another biological pathway important for the resolution of depressive symptoms in connection with hippocampal neurogenesis, the Wnt/β-catenin pathway [119,120]. That same pathway is also activated by the neuronal pentraxin 2 [121], encoded by the *Nptx2* gene that is upregulated by FLX (Tables 3, 9, S5, S8 and S12) as already demonstrated in a previous study focused on a shared mechanism of AD effect of chronic FLX and exercise [118].

c PI3K/AKT

The response to ADs involves several overlapping mechanisms and, in addition to the signaling mentioned above, the PI3K/AKT pathway can also be mentioned [122]. Thus, a dysregulation of *Akt3* expression is observed (Tables 7 and 10), knowing that the invalidation of this gene in mice revealed an endophenotype reminiscent of psychiatric manifestations such as schizophrenia, anxiety, and depression [123]. The gene *Tyro3*, which we found upregulated after FLX action (Tables 4, 7, 9 and S12), encodes the most widely expressed receptor protein tyrosine kinase and is linked to the PI3K/AKT pathway [124].

d Cyclic AMP

It is noteworthy that among the genes common to the signatures of FLX activity in a stress or naive context is the *Pde4b* gene (Tables 6 and S6, Figure 3). It encodes a cyclic phosphodiesterase, regulating concentrations of cyclic nucleotides, and thereby plays a role in signal transduction. Not only have several studies specifically shown variations in the expression of transcripts and protein encoded by *Pde4b* in preclinical models and humans [125-129], but this gene is in fact the target of a drug with AD properties, rolipram. A link between neurotrophin synthesis and the production of second messengers like cAMP is protein kinase A. It turns out that one of the genes most reliably affected by FLX codes for a subunit of this kinase, *Prkar1b* (Tables 2-4 and 9 and Figure 3).

e Glucorticoids

Perhaps a pivot of the effect of antidepressants on major depression through several previously mentioned signaling pathways, cAMP/PKA, phosphodiesterase, and neurotrophins, materializes through the activation of nuclear receptors, the glucocorticoid receptor (GR) and the mineralocorticoid (MR) [130–134]. Two genes share the task of regulating, among others, inflammatory responses, proliferation and differentiation processes, *NR3C1* and *NR3C2*. In our data, if it seems that the tendency is rather to underexpress these two transcription factors, *Nr3c2* is more consensually deregulated in the different signatures that we obtain (Tables S5, S8 and S12), without being one of the key markers either. On the other hand, considering the partners of the GR, we notice that AHI1, which regulates the nuclear translocation of the GR, and which has already been associated with stress-induced depressive behavior in mice [135], has its mRNA significantly underexpressed in the different signatures we have generated (Tables S5, S8 and S12).

f Synaptic plasticity

In a way, in fine, the beneficial action of SSRIs can be linked to a capacity to achieve synaptic plasticity and from this point of view it seems important to us to mention the decrease in expression of two mRNAs coding for transcription factors involved in brain development [136], *Nfia* and *Nfib* (Tables 2, 7, S5, S6, S8 and S11 and Figure 3). For the latter, a study in rats had shown that its expression variation was during the molecular signature at the level of the frontal cortex characterizing the AD action of quetiapine in a chronic
stress model in rats [137]. Recently, a study on lymphoblastoid cell lines from depression patients treated with citalopram reported a significant association of *NFIB* expression with improvement in depression scale [138]. The NMDA subtype of glutamate receptors (NMDAR) plays a key role in synaptic plasticity in the context of depression [139], and it has been shown that a $Ca^{2+}$/calmodulin-dependent Ras-guanine-nucleotide-releasing factor (RasGRF1) served as an NMDAR-dependent regulator of the ERK kinase pathway. In our meta-analysis, we noted that *Rasgrf1* was one of the most consistent downregulated genes by FLX activity (Tables 6–8 and 10). Consistently, it has been shown that *Rasgrf1* knockdown reversed the effect of UCMS on mice [140], and it has been shown in humans that RASGRF1 may be a potential specific biomarker of treatment response for bipolar disorder [141]. Concerning the NMDAR-positive modulation underlying AD action, one identified initial trigger of such an effect is Drd1-pyramidal cell signaling [142,143]. We uncovered a *Drd1* gene upregulation induced by FLX in our meta-analysis (Tables 4, 7 and 9) in agreement with these observations. Another gene common to the different signatures of FLX activity is *Pdlim5* (Table 6, Figure 3), which encodes a protein that tethers kinases to the Z-discs in striated muscles but also restrains the postsynaptic growth of excitatory synapses. Several investigations suggested that the human gene plays a key role in the pathophysiology of mood disorders but is also likely involved in AD response [144–146].

Overall, whatever the paradigm explored (stress vs. naïve), it appears that FLX has a marked effect on synaptic activity in the hippocampus (Tables 4, 8, S9 and S13). Thus, many genes characterizing the transcriptional signature of FLX code for proteins playing an important role at the synaptic level. This is the case for the *Sorcs1* gene [147], encoding the sortilin-related VPS10 domain containing receptor 1, upregulated by FLX (Tables 6 and S11) as previously observed [148]. Another example is the *Itsn1* gene (Tables 6, 7, 8, S9 and S12), encoding the intersectin 1 that regulates synaptic vesicle recycling [149]. Interestingly, we can even point out a link, perhaps not the most expected, between synapse and transcriptional activators. For example, the *Pcdh19* gene (Tables 6, 7, S7 and S12), which codes for a protocadherin, bridges neuronal activity with gene expression by regulating immediate-early genes expression to favor maintenance of neuronal homeostasis [150]. In addition, *Sox11* (Tables 7, S4, S6, S7 and S11) encodes a transcriptional factor playing a major role in brain development and activated in an activity-dependent fashion specific to the dentate gyrus of the hippocampus [151], which had already been shown to be dysregulated by FLX [118], as well as by paroxetine, another SSRI [100].

### 3.2.4. Emerging Pathways

A novelty in this meta-analysis is the zinc finger homeobox 2 gene, *Zfhx2* (Tables 2 and 6, Figure 3), for which there is generally very little literature. This putative transcriptional regulator is highly expressed in the developing brain and, very interestingly, its global knockout has been reported to show several behavioral abnormalities, namely hyperactivity, enhanced depression-like behaviors, and an aberrantly altered anxiety-like phenotype [152], as well as hyposensitivity to noxious mechanical stimuli [153]. A human variant of *ZFHX2* has been associated with a pain-insensitive phenotype [153]. Another gene related to pain and nociception is *Trpm3* (Tables 7 and S12, Figure 3), encoding a tetrameric cation channel. It had been shown that FLX but also imipramine and the antipsychotic drug chlorpromazine inhibit the channel activity [154]. Among the genes whose relevance in terms of belonging to a biological pathway is not at first obvious to relate to the pathophysiology of depression or the action of ADs is *Slc4a4* (Tables 7, S8 and S12). It encodes a sodium bicarbonate cotransporter (NBC) involved in the regulation of bicarbonate secretion and absorption and intracellular pH, but, interestingly, it has been proposed as one of the best biomarkers for predicting suicidal ideation [155,156]. Conversely, our analysis pinpoints to the oxytocin pathway altered in animals under stress exposure (Table 3) that has received renewed attention in regard to psychiatric disorders [157].

Last but not least, components of the ribosome that we saw as particularly affected in stressed animals responding to FLX (Table 5) are also affected in the study of Hori et al.,
who explored the blood transcriptome of three homogeneous groups, namely “resilient”, “vulnerable”, and “resistant” [158]. Remarkably, among genes dysregulated between the different groups are *RPL17*, *RPL21*, *RPL34*, *RPS15A*, and *RPS27*, which were all dysregulated in rodent models (Table 5). In another study, Guilloux et al. were able to identify a group of six genes with predictive value for response to citalopram treatment for major depression [159]. Half of the genes were ribosomal units including again *RPL17*, as well as *RPL24*, that we found dysregulated (Table 5).

# 4. Methods and Materials

## 4.1. Systematic Search of the Literature

A systematic search was conducted on PubMed on 19 April 2022. The search was limited to papers published after 2006 to exclude studies based on immature DNA microarrays technology. It used combinations of the following terms: (antidepressant OR fluoxetine) AND (microarray OR transcriptome data OR transcriptomic changes OR RNA-sequencing OR transcriptome profiling OR differential molecular signature OR large-scale gene expression OR bacTRAP) AND (mouse OR rat OR rodent OR depression-like) NOT (focused microarray OR cell line OR cancer OR tumor OR traditional medicine) NOT (review [publication type]). The filter “English language” was also activated. In addition, references from eligible studies were examined to include additional studies that would not be retrieved by the PubMed search. To compile a complete list of eligible studies (Table S1), we excluded the studies that did not use an AD or profile the action of an AD genome-wide. Because we focused on the in vivo action of AD in rodents as a model of major depression in humans, we also excluded studies performed in systems other than rodent. Namely, studies on other species, cells in culture, and on bacterial infection or brain injury models were excluded. We searched all eligible studies for either an ID number linked to the public download of full transcriptomic data or a table containing the complete data as supplementary information. We also searched NCBI resources, the Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/, accessed on 19 April 2022), the Sequence Read Archive (SRA, https://www.ncbi.nlm.nih.gov/sra, accessed on 19 April 2022), and the European repository ArrayExpress (https://www.ebi.ac.uk/arrayexpress/, accessed on 19 April 2022) for microarray and high-throughput sequencing data related to published and unpublished studies by typing “antidepressant AND (mouse OR rat)”. Finally, when full transcriptomic data were not publicly available, we emailed the corresponding authors to request raw data.

## 4.2. Bioinformatic Analyses

The analyses were produced using Snakemake (https://doi.org/10.12688/f1000research.29032.1, accessed on 19 April 2022) as workflow manager and R Bookdown for reporting.

### 4.2.1. Individual Reprocessing of Published Experiments—Microarrays

Normalized microarrays from seven studies (GSE84185, GSE43261, GSE56028, GSE118669, GSE54307, GSE6476, and GSE42940) were retrieved from GEO using GEOquery [160]. One last microarray was retrieved directly from the authors [50] and was normalized using RMA from the oligo R package [161]. The microarrays' latest annotations were retrieved from Ensembl Biomart for Affymetrix GeneChip Mouse Genome 430 2.0, Affymetrix GeneChip Mouse 430A 2.0, Affymetrix Rat Gene 1.0 ST arrays, and Agilent SurePrint G3 Mouse Gene Expression 8 × 60 K microarrays. When unavailable, annotations from the original authors were completed using megablast queries (GSE42940). Probes with mapping to multiple genes were discarded. Differential gene expression analyses were produced with Limma [162].

### 4.2.2. Individual Reprocessing of Published Experiments—RNA-Seq

Raw reads from four studies (SRP056481, SRP057486, SRP084288, and SRP131063) were retrieved from SRA using SRA tools. Reads were trimmed using Sickle (-q 20) [163]
and aligned on GRCm38 (for mouse samples) or Rnor6 (for rat samples) with STAR [164] using default settings. Gene-level counts were produced using subread featureCounts [165] with Ensembl release 102 annotations. Differential gene expression analyses were produced with Deseq2 [166] using apeglm shrinkage [167].

### 4.3. *Integration Analyses*

To compare differential expression analyses between different experiments, we kept only genes with *p* values below 0.05 and in the top 300 of any compared designs. For RNA-Seq datasets, the log2 fold change values have been calculated using the “normal” shrinkage type and not the apeglm from Deseq2 because the latter was considered too stringent for this kind of integration analysis. Plots were produced using the R ggplot2 package [168]. To be able to specify the general trend of the direction of expression variation of a gene, we were inspired by the procedure defined in the MetaVolcano tool, the vote counting approach [169]. This consists of counting the number of comparisons with *p* values below 0.05 showing overexpression and subtracting the number of comparisons with *p* values below 0.05 corresponding to underexpression. The resulting integer, whether positive, zero, or negative, is called the consensus score.

### 4.4. *Portraits of FLX Action*

As an alternative procedure to extract the transcriptional signature of FLX treatment, we applied the method proposed by Stephen Gammie to generate a portrait of depression [23]. The *p* value for each gene and for each individual microarray or RNA-Seq dataset used in integration analyses was -log10 transformed and then multiplied by the sign of the direction of change. Then, for each dataset, both upregulated and downregulated genes were ranked from low to high. With the 1000-increment cutoff, for each gene we counted how many times it received a ranking of 1000 or less among the upregulated genes. The same approach was applied to count 2000 or less, but that time the count number was multiplied by 0.1. We used this method up to 8000 by multiplying the sum by $10^{-(k-1)}$, where *k* is the number of the thousands of ranks considered. The same calculation was applied to downregulated genes. Finally, we calculated the sum of the above values for a given gene. Genes were sorted by the absolute value so that the final list contains all genes ordered by magnitude of the difference from controls with also information on direction of expression change.

### 4.5. *Ontological Analysis*

Gene lists were uploaded on DAVID (database for annotation, visualization and integrated discovery) Bioinformatics Resources (2021 Update) [170] for identifying statistically relevant biological processes. The default settings were kept unchanged and only functional annotations related to the gene ontology (GO) biological process, molecular process, and cellular component, as well as the KEGG pathway, with corrected *p*-value (Bonferroni) < 0.05, were retrieved. We limited the search to *Mus musculus* species annotation. To simplify the lists of biological processes obtained, we only kept, for identical gene lists, the process associated with the lowest *p*-value.

## 5. Conclusions

We provided the first meta-analysis of a transcriptional signature of FLX in animal models and highlighted the impact and pitfalls of gene expression variation studies. Future studies aiming at extending such signature to other widely prescribed ADs and considering other tissues and brain regions in both males and females are warranted to provide better insights into optimal targets for AD response in human.

Despite all the limitations mentioned before, our study shows how murine models are a valuable source for understanding the mechanisms involved in the response to antidepressants. Indeed, with the benefit of hindsight on the main processes regulated by FLX, we realize that there are solid points of convergence with what is described in the
context of the follow-up of patients suffering from major depression. It would obviously be
important to know whether the changes induced by FLX, in the context of an animal model
recapitulating the response to AD treatment, are similar to the biological program that is set
up in a depressed patient who responds favorably to AD treatment. The major concern is
that all human studies based on post-mortem tissues cannot distinguish the death-induced
signature from those induced by the treatment and those corresponding to the depressive
illness. In contrast, studies that allow for clinical interviews at the same time the biological
samples are taken, such as studies based on blood samples, can define the AD treatment-
related signature in relation to response or non-response to treatment. In fact, studies
on human cohorts have already underlined the importance of some biological processes
significantly enriched from ontological analyses of genome-wide transcriptional variations
in FLX-treated rodent models of major depression. That includes synaptic plasticity [171],
even when the variations in the corticolimbic regions of the brain are of opposite directions
between men and women [60,172], and in particular at the blood level with glutamatergic
signaling, the intervention of neurotrophin [173,174], and MAPK signaling [74].

Overall, our results highlight activity-dependent gene transcription as the major AD
target. This is a major finding with an important potential for future biomarker clinical
research. At the center of cellular plasticity and resilience, in the CNS and the periphery,
such cell signaling cascades are likely early-stage mediators of AD response body-wide.
The findings of this meta-analysis allow us to therefore hypothesize that ADs do induce
early changes in cellular biology and processes at the interface of CNS and periphery that
could constitute clinically relevant biomarkers and suggest that more focus is needed in
the earliest point of biomarker variations.

**Supplementary Materials:** The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms232113543/s1.

**Author Contributions:** E.C.I. and E.T.T. designed the study. E.C.I. administered the project and performed the ontological analyses. E.S. provided data. G.C. performed the bioinformatic analyses. E.C.I., R.B, E.T.T. and G.T. provided funding for the study. R.B. provided guidance on the use of AD modalities in clinical practice. V.G. and P.O.-T. provided input on animal model protocols. E.C.I. wrote the first draft of the manuscript. E.C.I. and E.T.T. wrote the paper. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by ERA-NET NEURON ANTaRES (https://www.neuron-eranet.eu/projects/ANTaRES/) with reference ANR-19-NEUR-0005 to E.T.T., R.B., G.T., E.C.I., V.G. and P.O.-T., and by Gebra (Imaging GEne and BRAin networks in stress susceptibility: novel pathways and biomarkers in mice and men) with reference ANR-19-CE37-0017 to E.T.T., R.B., G.T., E.C.I. and V.G.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The source code for all the analyses is available at https://github.com/guillaumecharbonnier/mw-ibrahim2022 and a compiled version of the Bookdown report is available at https://guillaumecharbonnier.github.io/mw-ibrahim2022.

**Acknowledgments:** We would like to thank all the researchers who provided access to raw transcriptional data through public portals.

**Conflicts of Interest:** E.S. is founder and acting chief officer of Damona Pharmaceuticals, a drug development company with small molecules in the pipeline for treatment of cognitive deficits across brain disorders and aging. E.T.T. is co-founder of Melkin Pharmaceuticals, a biotech company. Neither company had any role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. All other authors report no biomedical financial interests or potential conflicts of interest.
## References

1. Friedrich, M.J. Depression Is the Leading Cause of Disability around the World. *JAMA* **2017**, *317*, 1517. [CrossRef] [PubMed]

2. Malhi, G.S.; Mann, J.J. Depression. *Lancet* **2018**, *392*, 2299-2312. [CrossRef]

3. Hasin, D.S.; Sarvet, A.L.; Meyers, J.L.; Saha, T.D.; Ruan, W.J.; Stohl, M.; Grant, B.F. Epidemiology of Adult DSM-5 Major Depressive Disorder and Its Specifiers in the United States. *JAMA Psychiatry* **2018**, *75*, 336-346. [CrossRef]

4. Perlis, R.H. Abandoning personalization to get to precision in the pharmacotherapy of depression. *World Psychiatry* **2016**, *15*, 228-235. [CrossRef]

5. Cipriani, A.; Furukawa, T.A.; Salanti, G.; Chaimani, A.; Atkinson, L.Z.; Ogawa, Y.; Leucht, S.; Ruhe, H.G.; Turner, E.H.; Higgins, J.P.T.; et al. Comparative efficacy and acceptability of 21 antidepressant drugs for the acute treatment of adults with major depressive disorder: A systematic review and network meta-analysis. *Lancet* **2018**, *391*, 1357-1366. [CrossRef]

6. Casarotto, P.C.; Girych, M.; Fred, S.M.; Kovaleva, V.; Moliner, R.; Enkavi, G.; Biojone, C.; Cannarozzo, C.; Sahu, M.P.; Kaurinkoski, K.; et al. Antidepressant drugs act by directly binding to TRKB neurotrophin receptors. *Cell* **2021**, *184*, 1299-1313.e19. [CrossRef]

7. Trivedi, M.H.; Daly, E.J. Treatment strategies to improve and sustain remission in major depressive disorder. *Dialogues Clin. Neurosci.* **2008**, *10*, 377-384.

8. Kraus, C.; Kadriu, B.; Zarate, C.A., Jr.; Kasper, S. Prognosis and improved outcomes in major depression: A review. *Transl. Psychiatry* **2019**, *9*, 12-17. [CrossRef]

9. Le-Niculescu, H.; Roseberry, K.; Gill, S.S.; Levey, D.F.; Phalen, P.L.; Mullen, J.; Williams, A.; Bhoiro, S.; Voegtline, T.; Davis, H.; et al. Precision medicine for mood disorders: Objective assessment, risk prediction, pharmacogenomics, and repurposed drugs. *Mol. Psychiatry* **2021**, *26*, 2776-2804. [CrossRef]

10. Cook, I.A.; Congdon, E.; Krantz, D.E.; Hunter, A.M.; Coppola, G.; Hamilton, S.P.; Leuchter, A.F. Time Course of Changes in Peripheral Blood Gene Expression During Medication Treatment for Major Depressive Disorder. *Front. Genet.* **2019**, *10*, 870. [CrossRef]

11. Pisanu, C.; Severino, G.; De Toma, I.; Dierssen, M.; Fusar-Poli, P.; Gennarelli, M.; Lio, P.; Maffioletti, E.; Maron, E.; Mehta, D.; et al. Transcriptional biomarkers of response to pharmacological treatments in severe mental disorders: A systematic review. *Eur. Neuropsychopharmacol.* **2022**, *55*, 112-157. [CrossRef] [PubMed]

12. Mariani, N.; Cattane, N.; Pariante, C.; Cattaneo, A. Gene expression studies in Depression development and treatment: An overview of the underlying molecular mechanisms and biological processes to identify biomarkers. *Transl. Psychiatry* **2021**, *11*, 354. [CrossRef] [PubMed]

13. Carrillo-Roa, T.; Labermaier, C.; Weber, P.; Herzog, D.P.; Lareau, C.; Santarelli, S.; Wagner, K.V.; Rex-Haffner, M.; Harbich, D.; Scharf, S.H.; et al. Common genes associated with antidepressant response in mouse and man identify key role of glucocorticoid receptor sensitivity. *PLoS Biol.* **2017**, *15*, e2002690. [CrossRef]

14. Wang, Q.; Timberlake, M.A., 2nd; Prall, K.; Dwivedi, Y. The recent progress in animal models of depression. *Prog. Neuropsychopharmacol. Biol. Psychiatry* **2017**, *77*, 99-109. [CrossRef] [PubMed]

15. Baker, M.; Hong, S.-I.; Kang, S.; Choi, D.-S. Rodent models for psychiatric disorders: Problems and promises. *Lab. Anim. Res.* **2020**, *36*, 9. [CrossRef] [PubMed]

16. Scarpa, J.R.; Fatma, M.; Loh, Y.-H.E.; Traore, S.R.; Stefan, T.; Chen, T.H.; Nestler, E.J.; Labonté, B. Shared Transcriptional Signatures in Major Depressive Disorder and Mouse Chronic Stress Models. *Biol. Psychiatry* **2020**, *88*, 159-168. [CrossRef]

17. Cathomas, F.; Murrough, J.W.; Nestler, E.J.; Han, M.-H.; Russo, S.J. Neurobiology of Resilience: Interface Between Mind and Body. *Biol. Psychiatry* **2019**, *86*, 410-420. [CrossRef]

18. Lorsch, Z.S.; Hamilton, P.; Ramakrishnan, A.; Parise, E.M.; Salery, M.; Wright, W.; Lepack, A.E.; Mews, P.; Issler, O.; McKenzie, A.; et al. Stress resilience is promoted by a Zfp189-driven transcriptional network in prefrontal cortex. *Nat. Neurosci.* **2019**, *22*, 1413-1423. [CrossRef]

19. Dudek, K.A.; Dion-Albert, L.; Kaufmann, F.N.; Tuck, E.; Lebel, M.; Menard, C. Neurobiology of resilience in depression: Immune and vascular insights from human and animal studies. *Eur. J. Neurosci.* **2021**, *53*, 183-221. [CrossRef]

20. Krishnan, V.; Nestler, E.J. Animal Models of Depression: Molecular Perspectives. *Mol. Funct. Model. Neuropsychiatry* **2011**, *7*, 121-147. [CrossRef]

21. Hervé, M.; Bergon, A.; Le Guisquet, A.-M.; Leman, S.; Consoloni, J.-L.; Fernandez-Nunez, N.; Lefebvre, M.-N.; El-Hage, W.; Belzeaux, R.; Belzung, C.; et al. Translational Identification of Transcriptional Signatures of Major Depression and Antidepressant Response. *Front. Mol. Neurosci.* **2017**, *10*, 248. [CrossRef] [PubMed]

22. Gammie, S.C. Evaluation of animal model congruence to human depression based on large-scale gene expression patterns of the CNS. *Sci. Rep.* **2022**, *12*, 1-14. [CrossRef] [PubMed]

23. Gammie, S.C. Creation of a gene expression portrait of depression and its application for identifying potential treatments. *Sci. Rep.* **2021**, *11*, 3829. [CrossRef]

24. Rayan, N.A.; Kumar, V.; Aow, J.; Rastegar, N.; Lim, M.G.L.; O'Toole, N.; Aliwarga, E.; Arcego, D.M.; Yeo, H.T.G.; Wong, J.Y.; et al. Integrative multi-omics landscape of fluoxetine action across 27 brain regions reveals global increase in energy metabolism and region-specific chromatin remodelling. *Mol. Psychiatry* **2022**, 1-16. [CrossRef]

25. Carboni, L.; Marchetti, L.; Lauria, M.; Gass, P.; Vollmayr, B.; Redfern, A.; Jones, L.; Razzoli, M.; Malki, K.; Begni, V.; et al. Cross-species evidence from human and rat brain transcriptome for growth factor signaling pathway dysregulation in major depression. *Neuropsychopharmacology* **2018**, *43*, 2134-2145. [CrossRef] [PubMed]
26. Lin, E.; Tsai, S.-J. Genome-wide microarray analysis of gene expression profiling in major depression and antidepressant therapy. *Prog. Neuro-Psychopharmacology Biol. Psychiatry* **2016**, *64*, 334–340. [CrossRef]

27. Ramos-Da-Silva, L.; Carlson, P.T.; Silva-Costa, L.C.; Martins-De-Souza, D.; de Almeida, V. Molecular Mechanisms Associated with Antidepressant Treatment on Major Depression. *Complex Psychiatry* **2021**, *7*, 49–59. [CrossRef]

28. Kumar, S.; Gál, Z.; Gonda, X.; Huse, R.J.; Juhasz, G.; Bagdy, G.; Petschner, P. Transcriptomic changes following chronic administration of selective serotonin reuptake inhibitors: A review of animal studies. *Neuropsychopharmacol. Hung.* **2019**, *21*, 26–35.

29. Ramsteijn, A.S.; Verkaik-Schakel, R.N.; Houwing, D.J.; Plösch, T.; Olivier, J.D.A. Perinatal exposure to fluoxetine and maternal adversity affect myelin-related gene expression and epigenetic regulation in the corticolimbic circuit of juvenile rats. *Neuropsychopharmacology* **2022**, *47*, 1620–1632. [CrossRef]

30. Rajkumar, A.P.; Qvist, P.; Donskov, J.G.; Lazarus, R.; Pallesen, J.; Nava, N.; Winther, G.; Liebenberg, N.; la Cour, S.H.; Paternoster, V.; et al. Reduced Brd1 expression leads to reversible depression-like behaviors and gene-expression changes in female mice. *Transl. Psychiatry* **2020**, *10*, 239. [CrossRef]

31. Sheftel, C.M.; Liu, L.; Field, S.L.; Weaver, S.R.; Vezina, C.M.; Peñagaricano, F.; Hernandez, L.L. Impact of Fluoxetine Treatment and Folic Acid Supplementation on the Mammary Gland Transcriptome During Peak Lactation. *Front. Pharmacol.* **2022**, *13*, 828735. [CrossRef] [PubMed]

32. Malki, K.; Lourdusamy, A.; Binder, E.; Payá-Cano, J.; Sluyter, F.; Craig, I.; Keers, R.; McGuffin, P.; Uher, R.; Schalkwyk, L.C. Antidepressant-dependent mRNA changes in mouse associated with hippocampal neurogenesis in a mouse model of depression. *Pharmacogenet. Genom.* **2012**, *22*, 765–776. [CrossRef] [PubMed]

33. Malki, K.; Uher, R.; Paya-Cano, J.; Binder, E.; Rietschel, M.; Zobel, A.; Mors, O.; Hauser, J.; Henigsberg, N.; Jerman, B.; et al. Convergent Animal and Human Evidence Suggests a Role of PPM1A Gene in Response to Antidepressants. *Biol. Psychiatry* **2011**, *69*, 360–365. [CrossRef] [PubMed]

34. Malki, K.; Tosto, M.G.; Mouriño-Talín, H.; Rodríguez-Lorenzo, S.; Pain, O.; Jumhaboy, I.; Liu, T.; Parpas, P.; Newman, S.; Malykh, A.; et al. Highly polygenic architecture of antidepressant treatment response: Comparative analysis of SSRI and NRI treatment in an animal model of depression. *Am. J. Med. Genet. Part B Neuropsychiatr. Genet.* **2017**, *174*, 235–250. [CrossRef] [PubMed]

35. Akkouh, I.A.; Skrede, S.; Holmgren, A.; Ersland, K.M.; Hansson, L.; Bahrami, S.; Andreassen, O.A.; Steen, V.M.; Djurovic, S.; Hughes, T. Exploring lithium's transcriptional mechanisms of action in bipolar disorder: A multi-step study. *Neuropsychopharmacology* **2020**, *45*, 947–955. [CrossRef] [PubMed]

36. Park, J.-Y.; Lee, E.-H.; Kwon, H.-J.; Baek, I.-S.; Seo, J.-S.; Kim, K.-S.; Han, P.-L. Behavioral Engagement With Playable Objects Resolves Stress-Induced Adaptive Changes by Reshaping the Reward System. *Biol. Psychiatry* **2022**, *91*, 676–689. [CrossRef] [PubMed]

37. Frazer, S.; Otomo, K.; Dayer, A. Early-life serotonin dysregulation affects the migration and positioning of cortical interneuron subtypes. *Transl. Psychiatry* **2015**, *5*, e644. [CrossRef]

38. Miao, N.; Jin, J.; Kim, S.-N.; Sun, T. Hippocampal MicroRNAs Respond to Administration of Antidepressant Fluoxetine in Adult Mice. *Int. J. Mol. Sci.* **2018**, *19*, 671. [CrossRef]

39. Ulrich-Merzenich, G.; Kelber, O.; Koptina, A.; Freischmidt, A.; Heilmann, J.; Müller, J.; Zeitler, H.; Seidel, M.; Ludwig, M.; Heinrich, E.; et al. Novel neurological and immunological targets for salicylate-based phytopharmaceuticals and for the anti-depressant imipramine. *Phytomedicine* **2012**, *19*, 930–939. [CrossRef]

40. Ulrich-Merzenich, G.; Koptina, A.; Kelber, O.; Freischmidt, A.; Heilmann, J.; Müller, J.; Sadeghlar, F.; Zeitler, H.; Wagner, H. Prediction of adverse events by in vivo gene expression profiling exemplified for phytopharmaceuticals containing salicylates and the antidepressant imipramine. *Phytomedicine* **2012**, *19*, 322–329. [CrossRef]

41. Iegaki, N.; Narita, Y.; Hattori, N.; Hirata, Y.; Ichihara, K. Royal jelly reduces depression-like behavior through possible effects on adrenal steroidogenesis in a murine model of unpredictable chronic mild stress. *Biosci. Biotechnol. Biochem.* **2020**, *84*, 606–612. [CrossRef] [PubMed]

42. Igarashi, Y.; Nakatsu, N.; Yamashita, T.; Ono, A.; Ohno, Y.; Urushidani, T.; Yamada, H. Open TG-GATEs: A large-scale toxicogenomics database. *Nucleic Acids Res.* **2015**, *43*, D921–D927. [CrossRef] [PubMed]

43. Kawamoto, T.; Ito, Y.; Morita, O.; Honda, H. Mechanism-based risk assessment strategy for drug-induced cholestasis using the transcriptional benchmark dose derived by toxicogenomics. *J. Toxicol. Sci.* **2017**, *42*, 427–436. [CrossRef] [PubMed]

44. Sakaida, M.; Sukeno, M.; Imoto, Y.; Tsuchiya, S.; Sugimoto, Y.; Okuno, Y.; Segi-Nishida, E. Electroconvulsive seizure-induced changes in gene expression in the mouse hypothalamic paraventricular nucleus. *J. Psychopharmacol.* **2013**, *27*, 1058–1069. [CrossRef]

45. Schmidt, E.F.; Warner-Schmidt, J.L.; Otopalik, B.G.; Pickett, S.B.; Greengard, P.; Heintz, N. Identification of the Cortical Neurons that Mediate Antidepressant Responses. *Cell* **2012**, *149*, 1152–1163. [CrossRef]

46. Sargin, D.; Chottekalapanda, R.U.; Perit, K.E.; Yao, V.; Chu, D.; Sparks, D.W.; Kalik, S.; Power, S.K.; Troyanskaya, O.G.; Schmidt, E.F.; et al. Mapping the physiological and molecular markers of stress and SSRI antidepressant treatment in S100a10 corticostriatal neurons. *Mol. Psychiatry* **2020**, *25*, 1112–1129. [CrossRef]

47. Bagot, R.C.; Cates, H.M.; Purushothaman, I.; Vialou, V.; Heller, E.A.; Yieh, L.; LaBonté, B.; Peña, C.J.; Shen, L.; Wittenberg, G.M.; et al. Ketamine and Imipramine Reverse Transcriptional Signatures of Susceptibility and Induce Resilience-Specific Gene Expression Profiles. *Biol. Psychiatry* **2017**, *81*, 285–295. [CrossRef]
48. Kao, C.-Y.; He, Z.; Zannas, A.S.; Hahn, O.; Kühne, C.; Reichel, J.M.; Binder, E.B.; Wotjak, C.T.; Khaitovich, P.; Turck, C.W. Fluoxetine treatment prevents the inflammatory response in a mouse model of posttraumatic stress disorder. *J. Psychiatr. Res.* **2016**, *76*, 74–83. [CrossRef]

49. Conti, B.; Maier, R.; Barr, A.; Morale, M.C.; Lu, X.; Sanna, P.P.; Bilbe, G.; Hoyer, D.; Bartfai, T. Region-specific transcriptional changes following the three antidepressant treatments electro convulsive therapy, sleep deprivation and fluoxetine. *Mol. Psychiatry* **2007**, *12*, 167–189. [CrossRef]

50. Surget, A.; Wang, Y.; Leman, S.; Ibarguen-Vargas, Y.; Edgar, N.; Griebel, G.; Belzung, C.; Sibille, E. Corticolimbic Transcriptional Changes are State-Dependent and Region-Specific in a Rodent Model of Depression and of Antidepressant Reversal. *Neuropsychopharmacology* **2009**, *34*, 1363–1380. [CrossRef]

51. Samuels, B.A.; Leonardo, E.D.; Dranovsky, A.; Williams, A.; Wong, E.; Nesbitt, A.M.; McCurdy, R.D.; Hen, R.; Alter, M. Global state measures of the dentate gyrus gene expression system predict antidepressant-sensitive behaviors. *PLoS ONE* **2014**, *9*, e85136. [CrossRef] [PubMed]

52. Patrício, P.; Mateus-Pinheiro, A.; Irmler, M.; Alves, N.D.; Machado-Santos, A.R.; Morais, M.; Correia, J.S.; Korostynski, M.; Piechota, M.; Stoffel, R.; et al. Differential and Converging Molecular Mechanisms of Antidepressants' Action in the Hippocampal Dentate Gyrus. *Neuropsychopharmacology* **2015**, *40*, 338–349. [CrossRef]

53. Wang, Y.; Jiang, H.; Meng, H.; Lu, J.; Li, J.; Zhang, X.; Yang, X.; Zhao, B.; Sun, Y.; Bao, T. Genome-wide transcriptome analysis of hippocampus in rats indicated that TLR/NLR signaling pathway was involved in the pathogenesis of depressive disorder induced by chronic restraint stress. *Brain Res. Bull.* **2017**, *134*, 195–204. [CrossRef]

54. Hagihara, H.; Ohira, K.; Miyakawa, T. Transcriptomic evidence for immaturity induced by antidepressant fluoxetine in the hippocampus and prefrontal cortex. *Neuropsychopharmacol. Rep.* **2019**, *39*, 78–89. [CrossRef] [PubMed]

55. Imoto, Y.; Segi-Nishida, E.; Suzuki, H.; Kobayashi, K. Rapid and stable changes in maturation-related phenotypes of the adult hippocampal neurons by electroconvulsive treatment. *Mol. Brain* **2017**, *10*, 8. [CrossRef] [PubMed]

56. Miller, B.H.; Schultz, L.E.; Gulati, A.; Cameron, M.D.; Pletcher, M.T. Genetic Regulation of Behavioral and Neuronal Responses to Fluoxetine. *Neuropsychopharmacology* **2008**, *33*, 1312–1322. [CrossRef] [PubMed]

57. Kroeze, Y.; Peeters, D.; Boulle, F.; Hove, D.L.A.V.D.; van Bokhoven, H.; Zhou, H.; Homberg, J.R. Long-term consequences of chronic fluoxetine exposure on the expression of myelination-related genes in the rat hippocampus. *Transl. Psychiatry* **2015**, *5*, e642. [CrossRef]

58. Sarkar, A.; Chachra, P.; Kennedy, P.; Pena, C.J.; Desouza, L.A.; Nestler, E.J.; Vaidya, V.A. Hippocampal HDAC4 Contributes to Postnatal Fluoxetine-Evoked Depression-Like Behavior. *Neuropsychopharmacology* **2014**, *39*, 2221–2232. [CrossRef]

59. Tellam, R.L.; Rushton, P.; Schuerman, P.; Pala, I.; Anane, D. The primary reasons behind data sharing, its wider benefits and how to cope with the realities of commercial data. *BMC Genom.* **2015**, *16*, 626. [CrossRef]

60. Seney, M.L.; Glausier, J.; Sibille, E. Large-Scale Transcriptomics Studies Provide Insight Into Sex Differences in Depression. *Biol. Psychiatry* **2022**, *91*, 14–24. [CrossRef]

61. Malberg, J.E.; Blendy, J.A. Antidepressant action: To the nucleus and beyond. *Trends Pharmacol. Sci.* **2005**, *26*, 631–638. [CrossRef] [PubMed]

62. Kim, W.; Chung, C. Brain-wide cellular mapping of acute stress-induced activation in male and female mice. *FASEB J.* **2021**, *35*, e22041. [CrossRef] [PubMed]

63. Monteggia, L.M.; Luikart, B.; Barrot, M.; Theobold, D.; Malkovska, I.; Nef, S.; Parada, L.F.; Nestler, E.J. Brain-Derived Neurotrophic Factor Conditional Knockouts Show Gender Differences in Depression-Related Behaviors. *Biol. Psychiatry* **2007**, *61*, 187–197. [CrossRef]

64. Tylee, A.; Walters, P. Onset of action of antidepressants. *BMJ* **2007**, *334*, 911–912. [CrossRef] [PubMed]

65. Moncrieff, J.; Cooper, R.E.; Stockmann, T.; Amendola, S.; Hengartner, M.P.; Horowitz, M.A. The serotonin theory of depression: A systematic umbrella review of the evidence. *Mol. Psychiatry* **2022**, *1–14*. [CrossRef]

66. Benton, C.; Wiltshire, T. Biological alterations in depression. In *Psychiatric Disorders—Trends and Developments*; Uehara, T., Ed.; Intechopen: Rijeka, Croatia, 2011; pp. 223–266.

67. Apazoglou, K.; Farley, S.; Gorgievski, V.; Belzeaux, R.; López, J.P.; Grenier, J.; Ibrahim, E.C.; El Khoury, M.-A.; Tse, Y.C.; Mongredien, R.; et al. Antidepressive effects of targeting ELK-1 signal transduction. *Nat. Med.* **2018**, *24*, 591–597. [CrossRef]

68. Fanselow, M.S.; Dong, H.W. Are the dorsal and ventral hippocampus functionally distinct structures? *Neuron* **2010**, *65*, 7–19. [CrossRef]

69. Ding, S.-L.; Yao, Z.; Hirokawa, K.E.; Nguyen, T.N.; Graybuck, L.T.; Fong, O.; Bohn, P.; Ngo, K.; Smith, K.A.; Koch, C.; et al. Distinct Transcriptomic Cell Types and Neural Circuits of the Subiculum and Prosubiculum along the Dorsal-Ventral Axis. *Cell Rep.* **2020**, *31*, 107648. [CrossRef]

70. Floriou-Servou, A.; von Ziegler, L.; Stalder, L.; Sturman, O.; Privitera, M.; Rassi, A.; Cremonesi, A.; Thöny, B.; Bohacek, J. Distinct Proteomic, Transcriptomic, and Epigenetic Stress Responses in Dorsal and Ventral Hippocampus. *Biol. Psychiatry* **2018**, *84*, 531–541. [CrossRef]

71. Žurawek, D.; Turecki, G. The miRNome of Depression. *Int. J. Mol. Sci.* **2021**, *22*, 11312. [CrossRef]

72. Li, Q.S.; Galbraith, D.; Morrison, R.L.; Trivedi, M.H.; Drevets, W.C. Circulating microRNA associated with future relapse status in major depressive disorder. *Front. Psychiatry* **2022**, *13*, 937360. [CrossRef] [PubMed]
73. O'Connor, R.M.; Grenham, S.; Dinan, T.G.; Cryan, J.F. microRNAs as novel antidepressant targets: Converging effects of ketamine and electroconvulsive shock therapy in the rat hippocampus. *Int. J. Neuropsychopharmacol.* **2013**, *16*, 1885–1892. [CrossRef] [PubMed]

74. Lopez, J.P.; Fiori, L.M.; Cruceanu, C.; Lin, R.; LaBonte, B.; Cates, H.M.; Heller, E.A.; Vialou, V.; Ku, S.M.; Gerald, C.; et al. MicroRNAs 146a/b-5 and 425-3p and 24-3p are markers of antidepressant response and regulate MAPK/Wnt-system genes. *Nat. Commun.* **2017**, *8*, 15497. [CrossRef] [PubMed]

75. Belzeaux, R.; Fiori, L.M.; Lopez, J.P.; Boucekine, M.; Boyer, L.; Blier, P.; Farzan, F.; Frey, B.N.; Giacobbe, P.; Lam, R.W.; et al. Predicting Worsening Suicidal Ideation With Clinical Features and Peripheral Expression of Messenger RNA and MicroRNA During Antidepressant Treatment. *J. Clin. Psychiatry* **2019**, *80*, 22559. [CrossRef]

76. Yrondi, A.; Fiori, L.M.; Frey, B.N.; Lam, R.W.; MacQueen, G.M.; Milev, R.; Müller, D.J.; Foster, J.A.; Kennedy, S.H.; Turecki, G. Association Between Side Effects and Blood microRNA Expression Levels and Their Targeted Pathways in Patients With Major Depressive Disorder Treated by a Selective Serotonin Reuptake Inhibitor, Escitalopram: A CAN-BIND-1 Report. *Int. J. Neuropsychopharmacol.* **2020**, *23*, 88–95. [CrossRef]

77. Nojima, T.; Proudfoot, N.J. Mechanisms of lncRNA biogenesis as revealed by nascent transcriptomics. *Nat. Rev. Mol. Cell Biol.* **2022**, *23*, 389–406. [CrossRef]

78. Roy, B.; Wang, Q.; Dwivedi, Y. Long Noncoding RNA-Associated Transcriptomic Changes in Resiliency or Susceptibility to Depression and Response to Antidepressant Treatment. *Int. J. Neuropsychopharmacol.* **2018**, *21*, 461–472. [CrossRef]

79. Zhang, C.-L.; Li, Y.-J.; Lu, S.; Zhang, T.; Xiao, R.; Luo, H.-R. Fluoxetine ameliorates depressive symptoms by regulating lncRNA expression in the mouse hippocampus. *Zool. Res.* **2021**, *42*, 28–42. [CrossRef]

80. Liu, L.; Wang, H.; Chen, X.; Zhang, Y.; Li, W.; Rao, X.; Liu, Y.; Zhao, L.; Pu, J.; Gui, S.; et al. Integrative Analysis of Long Non-coding RNAs, Messenger RNAs, and MicroRNAs Indicates the Neurodevelopmental Dysfunction in the Hippocampus of Gut Microbiota-Dysbiosis Mice. *Front. Mol. Neurosci.* **2021**, *14*, 745437. [CrossRef]

81. Hao, W.-Z.; Chen, Q.; Wang, L.; Tao, G.; Gan, H.; Deng, L.-J.; Huang, J.-Q.; Chen, J.-X. Emerging roles of long non-coding RNA in depression. *Prog. Neuro-Psychopharmacol. Biol. Psychiatry* **2022**, *115*, 110515. [CrossRef]

82. Zimmerman, A.J.; Hafez, A.K.; Amoah, S.K.; Rodriguez, B.A.; Dell'Orco, M.; Lozano, E.; Hartley, B.J.; Alural, B.; Lalonde, J.; Chander, P.; et al. A psychiatric disease-related circular RNA controls synaptic gene expression and cognition. *Mol. Psychiatry* **2020**, *25*, 2712–2727. [CrossRef] [PubMed]

83. Mao, J.; Li, T.; Fan, D.; Zhou, H.; Feng, J.; Liu, L.; Zhang, C.; Wang, X. Abnormal expression of rno_circRNA_014900 and rno_circRNA_005442 induced by ketamine in the rat hippocampus. *BMC Psychiatry* **2020**, *20*, 1. [CrossRef]

84. Lesch, K.; Aulakh, C.S.; Wolozin, B.L.; Tolliver, T.J.; Hill, J.L.; Murphy, D.L. Regional brain expression of serotonin transporter mRNA and its regulation by reuptake inhibiting antidepressants. *Mol. Brain Res.* **1993**, *17*, 31–35. [CrossRef]

85. Neumaier, J.F.; Root, D.C.; Hamblin, M.W. Chronic Fluoxetine Reduces Serotonin Transporter mRNA and 5-HT1B mRNA in a Sequential Manner in the Rat Dorsal Raphe Nucleus. *Neuropsychopharmacology* **1996**, *15*, 515–522. [CrossRef]

86. Koed, K.; Linnet, K. The serotonin transporter messenger RNA level in rat brain is not regulated by antidepressants. *Biol. Psychiatry* **1997**, *42*, 1177–1180. [CrossRef]

87. Maekawa, T.; Kim, S.; Nakai, D.; Makino, C.; Takagi, T.; Ogura, H.; Yamada, K.; Chatton, B.; Ishii, S. Social isolation stress induces ATF-7 phosphorylation and impairs silencing of the 5-HT 5B receptor gene. *EMBO J.* **2010**, *29*, 196–208. [CrossRef]

88. Rostene, W.; Fischette, C.; Rainbow, T.; McEwen, B. Modulation by vasoactive intestinal peptide of serotonin1 receptors in the dorsal hippocampus of the rat brain: An autoradiographic study. *Neurosci. Lett.* **1983**, *37*, 143–148. [CrossRef]

89. Etgen, A.; Browning, E. Activators of cyclic adenosine 3':5'-monophosphate accumulation in rat hippocampal slices: Action of vasoactive intestinal peptide (VIP). *J. Neurosci.* **1983**, *3*, 2487–2493. [CrossRef]

90. Rostene, W.; Fischette, C.; McEwen, B. Modulation by vasoactive intestinal peptide (VIP) of serotonin receptors in membranes from rat hippocampus. *J. Neurosci.* **1983**, *3*, 2414–2419. [CrossRef]

91. Huang, Y.-L.; Yu, J.-P.; Wang, G.-H.; Chen, Z.-H.; Wang, Q.; Xiao, L. Effect of fluoxetine on depression-induced changes in the expression of vasoactive intestinal polypeptide and corticotrophin releasing factor in rat duodenum. *World J. Gastroenterol.* **2007**, *13*, 6060–6065. [CrossRef]

92. Ibi, D.; Takuma, K.; Koike, H.; Mizoguchi, H.; Tsuritani, K.; Kuwahara, Y.; Kamei, H.; Nagai, T.; Yoneda, Y.; Nabeshima, T.; et al. Social isolation rearing-induced impairment of the hippocampal neurogenesis is associated with deficits in spatial memory and emotion-related behaviors in juvenile mice. *J. Neurochem.* **2008**, *105*, 921–932. [CrossRef] [PubMed]

93. Ieraci, A.; Mallei, A.; Popoli, M. Social Isolation Stress Induces Anxious-Depressive-Like Behavior and Alterations of Neuroplasticity-Related Genes in Adult Male Mice. *Neural Plast.* **2016**, *2016*, 1–13. [CrossRef] [PubMed]

94. Anderson, S.T.; Commins, S.; Moynagh, P.N.; Coogan, A.N. Lipopolysaccharide-induced sepsis induces long-lasting affective changes in the mouse. *Brain Behav. Immun.* **2015**, *43*, 98–109. [CrossRef]

95. Navailles, S.; Zimnisky, R.; Schmauss, C. Expression of Glucocorticoid Receptor and Early Growth Response Gene 1 during Postnatal Development of Two Inbred Strains of Mice Exposed to Early Life Stress. *Dev. Neurosci.* **2010**, *32*, 139–148. [CrossRef] [PubMed]

96. Monsey, M.S.; Boyle, L.M.; Zhang, M.L.; Nguyen, C.P.; Kronman, H.G.; Ota, K.T.; Duman, R.S.; Taylor, J.R.; Schafe, G.E. Chronic Corticosterone Exposure Persistently Elevates the Expression of Memory-Related Genes in the Lateral Amygdala and Enhances the Consolidation of a Pavlovian Fear Memory. *PLoS ONE* **2014**, *9*, e91530. [CrossRef] [PubMed]
97. Jeanneteau, F.; Barrère, C.; Vos, M.; De Vries, C.J.; Rouillard, C.; Levesque, D.; Dromard, Y.; Moisan, M.-P.; Duric, V.; Franklin, T.C.; et al. The Stress-Induced Transcription Factor NR4A1 Adjusts Mitochondrial Function and Synapse Number in Prefrontal Cortex. *J. Neurosci.* **2018**, *38*, 1335–1350. [CrossRef] [PubMed]

98. Vialou, V.; Thibault, M.; Kaska, S.; Cooper, S.; Gajewski, P.; Eagle, A.; Mazei-Robison, M.; Nestler, E.J.; Robison, A. Differential induction of FosB isoforms throughout the brain by fluoxetine and chronic stress. *Neuropharmacology* **2015**, *99*, 28–37. [CrossRef]

99. Sibille, E.; Wang, Y.; Joeyen-Waldorf, J.; Gaiteri, C.; Surget, A.; Oh, S.; Belzung, C.; Tseng, G.C.; Lewis, D. A Molecular Signature of Depression in the Amygdala. *Am. J. Psychiatry* **2009**, *166*, 1011–1024. [CrossRef]

100. Sillaber, I.; Panhuysen, M.; Henniger, M.S.H.; Ohl, F.; Kühne, C.; Pütz, B.; Pohl, T.; Deussing, J.M.; Paez-Pereda, M.; Holsboer, F. Profiling of behavioral changes and hippocampal gene expression in mice chronically treated with the SSRI paroxetine. *Psychopharmacology* **2008**, *200*, 557–572. [CrossRef]

101. Molteni, R.; Calabrese, F.; Mancini, M.; Racagni, G.; Riva, M.A. Basal and stress-induced modulation of activity-regulated cytoskeletal associated protein (Arc) in the rat brain following duloxetine treatment. *Psychopharmacology* **2008**, *201*, 285–292. [CrossRef]

102. Brivio, P.; Corsini, G.; Riva, M.A.; Calabrese, F. Chronic vortioxetine treatment improves the responsiveness to an acute stress acting through the ventral hippocampus in a glucocorticoid-dependent way. *Pharmacol. Res.* **2019**, *142*, 14–21. [CrossRef] [PubMed]

103. Zhang, W.-J.; Wang, H.-H.; Lv, Y.-D.; Liu, C.-C.; Sun, W.-Y.; Tian, L.-J. Downregulation of Egr-1 Expression Level via GluN2B Underlies the Antidepressant Effects of Ketamine in a Chronic Unpredictable Stress Animal Model of Depression. *Neuroscience* **2018**, *372*, 38–45. [CrossRef] [PubMed]

104. Erburu, M.; Cajaleon, L.; Guruceaga, E.; Venzala, E.; Muñoz-Cobo, I.; Beltrán, E.; Puerta, E.; Tordera, R. Chronic mild stress and imipramine treatment elicit opposite changes in behavior and in gene expression in the mouse prefrontal cortex. *Pharmacol. Biochem. Behav.* **2015**, *135*, 227–236. [CrossRef] [PubMed]

105. Waller, J.A.; Tamm, J.A.; Abdourahman, A.; Pehrson, A.L.; Li, Y.; Cajina, M.; Sánchez, C. Chronic vortioxetine treatment in rodents modulates gene expression of neurodevelopmental and plasticity markers. *Eur. Neuropsychopharmacol.* **2017**, *27*, 192–203. [CrossRef]

106. Scarpa, J.R.; Jiang, P.; Gao, V.D.; Fitzpatrick, K.; Millstein, J.; Olker, C.; Gotter, A.; Winrow, C.J.; Renger, J.J.; Kasarskis, A.; et al. Cross-species systems analysis identifies gene networks differentially altered by sleep loss and depression. *Sci. Adv.* **2018**, *4*, eaat1294. [CrossRef]

107. Muzio, L.; Brambilla, V.; Calcaterra, L.; D'Adamo, P.; Martino, G.; Benedetti, F. Increased neuroplasticity and hippocampal microglia activation in a mice model of rapid antidepressant treatment. *Behav. Brain Res.* **2016**, *311*, 392–402. [CrossRef]

108. Papp, M.; Gruca, P.; Faron-Górecka, A.; Kusmider, M.; Willner, P. Genomic Screening of Wistar and Wistar-Kyoto Rats Exposed to Chronic Mild Stress and Deep Brain Stimulation of Prefrontal Cortex. *Neuroscience* **2019**, *423*, 66–75. [CrossRef] [PubMed]

109. Savitz, J.; Frank, M.B.; Victor, T.; Bebak, M.; Marino, J.H.; Bellgowan, P.S.; McKinney, B.A.; Bodurka, J.; Teague, K.; Drevets, W.C. Inflammation and neurological disease-related genes are differentially expressed in depressed patients with mood disorders and correlate with morphometric and functional imaging abnormalities. *Brain Behav. Immun.* **2013**, *31*, 161–171. [CrossRef]

110. Manji, H.K.; Quiroz, J.A.; Sporn, J.; Payne, J.L.; Denicoff, K.; Gray, N.A.; Zarate, C.A., Jr.; Charney, D.S. Enhancing neuronal plasticity and cellular resilience to develop novel, improved therapeutics for difficult-to-treat depression. *Biol. Psychiatry* **2003**, *53*, 707–742. [CrossRef]

111. Svenningsson, P.; Tzavara, E.T.; Witkin, J.M.; Fienberg, A.A.; Nomikos, G.G.; Greengard, P. Involvement of striatal and extrastriatal DARPP-32 in biochemical and behavioral effects of fluoxetine (Prozac). *Proc. Natl. Acad. Sci. USA* **2002**, *99*, 3182–3187. [CrossRef]

112. Crozatier, C.; Farley, S.; Mansuy, I.; Dumas, S.; Giros, B.; Tzavara, E. Calcineurin (protein phosphatase 2B) is involved in the mechanisms of action of antidepressants. *Neuroscience* **2007**, *144*, 1470–1476. [CrossRef] [PubMed]

113. Nibuya, M.; Nestler, E.; Duman, R. Chronic antidepressant administration increases the expression of cAMP response element binding protein (CREB) in rat hippocampus. *J. Neurosci.* **1996**, *16*, 2365–2372. [CrossRef] [PubMed]

114. Gould, T.D.; O'Donnell, K.C.; Picchini, A.M.; Dow, E.R.; Chen, G.; Manji, H.K. Generation and behavioral characterization of β-catenin forebrain-specific conditional knock-out mice. *Behav. Brain Res.* **2008**, *189*, 117–125. [CrossRef]

115. Urigüen, L.; Arteta, D.; Díez-Alarcia, R.; Diaz, A.; Pazos, A.; Meana, J.J.; Ferrer-Alcón, M. Gene expression patterns in brain cortex of three different animal models of depression. *Genes Brain Behav.* **2008**, *7*, 649–658. [CrossRef] [PubMed]

116. Kulikova, E.; Kulikov, A.; Kulikov, E.K.A.A. Striatal-enriched Tyrosine Protein Phosphatase (STEP) in the Mechanisms of Depressive Disorders. *Curr. Protein Pept. Sci.* **2017**, *18*, 1152–1162. [CrossRef]

117. Kulikova, E.; Khotskin, N.; Illarionova, N.; Sorokin, I.; Bazhenova, E.; Kondaurova, E.; Volcho, K.; Khomenko, T.; Salakhutdinov, N.; Ponimaskin, E.; et al. Inhibitor of Striatal-Enriched Protein Tyrosine Phosphatase, 8-(Trifluoromethyl)-1,2,3,4,5-Benzopentathiepin-6-Amine hydrochloride (TC-2153), Produces Antidepressant-Like Effect and Decreases Functional Activity and Protein Level of 5-HT2A Receptor in the Brain. *Neuroscience* **2018**, *394*, 220–231. [CrossRef]

118. Huang, G.-J.; Ben-David, E.; Piella, A.T.; Edwards, A.; Flint, J.; Shifman, S. Neurogenomic Evidence for a Shared Mechanism of the Antidepressant Effects of Exercise and Chronic Fluoxetine in Mice. *PLoS ONE* **2012**, *7*, e35901. [CrossRef]

119. Niu, H.; Gou, R.; Xu, Q.; Duan, D. Recombinant insulin-like growth factor binding protein-4 inhibits proliferation and promotes differentiation of neural progenitor cells. *Neurosci. Lett.* **2017**, *642*, 71–76. [CrossRef]
120. Pilar-Cuellar, F.; Vidal, R.; Diaz, A.; Castro, E.; Anjos, S.; Vargas, V.; Romero, B.; Valdizan, E. Signaling Pathways Involved in Antidepressant-Induced Cell Proliferation and Synaptic Plasticity. *Curr. Pharm. Des.* **2014**, *20*, 3776-3794. [CrossRef]

121. Xu, C.; Tian, G.; Jiang, C.; Xue, H.; Kuerbanjiang, M.; Sun, L.; Gu, L.; Zhou, H.; Liu, Y.; Zhang, Z.; et al. NPTX2 promotes colorectal cancer growth and liver metastasis by the activation of the canonical Wnt/β-catenin pathway via FZD6. *Cell Death Dis.* **2019**, *10*, 217. [CrossRef]

122. Huang, W.; Zhao, Y.; Zhu, X.; Cai, Z.; Wang, S.; Yao, S.; Qi, Z.; Xie, P. Fluoxetine upregulates phosphorylated-AKT and phosphorylated-ERK1/2 proteins in neural stem cells: Evidence for a crosstalk between AKT and ERK1/2 pathways. *J. Mol. Neurosci.* **2013**, *49*, 244-249. [CrossRef] [PubMed]

123. Bergeron, Y.; Bureau, G.; Laurier-Laurin, M.É.; Asselin, E.; Massicotte, G.; Cyr, M. Genetic Deletion of Akt3 Induces an Endophenotype Reminiscent of Psychiatric Manifestations in Mice. *Front. Mol. Neurosci.* **2017**, *10*, 102. [CrossRef] [PubMed]

124. Prieto, A.; O'Dell, S.; Varnum, B.; Lai, C. Localization and signaling of the receptor protein tyrosine kinase Tyro3 in cortical and hippocampal neurons. *Neuroscience* **2007**, *150*, 319-334. [CrossRef]

125. Takahashi, M.; Terwilliger, R.; Lane, C.; Mezes, P.S.; Conti, M.; Duman, R.S. Chronic antidepressant administration increases the expression of cAMP-specific phosphodiesterase 4A and 4B isoforms. *J. Neurosci.* **1999**, *19*, 610-618. [CrossRef] [PubMed]

126. Miró, X.; Pérez-Torres, S.; Artigas, F.; Puigdomènech, P.; Palacios, J.M.; Mengod, G. Regulation of cAMP phosphodiesterase mRNAs expression in rat brain by acute and chronic fluoxetine treatment. An in situ hybridization study. *Neuropharmacology* **2002**, *43*, 1148-1157. [CrossRef]

127. Dlaboga, D.; Hajjhussein, H.; O'Donnell, J.M. Regulation of phosphodiesterase-4 (PDE4) expression in mouse brain by repeated antidepressant treatment: Comparison with rolipram. *Brain Res.* **2006**, *1096*, 104-112. [CrossRef]

128. Fatemi, S.H.; Folsom, T.D.; Reutiman, T.J.; Braun, N.N.; Lavergne, L.G. Levels of phosphodiesterase 4A and 4B are altered by chronic treatment with psychotropic medications in rat frontal cortex. *Synapse* **2010**, *64*, 550-555. [CrossRef]

129. Fujita, M.; Richards, E.M.; Niciu, M.; Ionescu, D.F.; Zoghbi, S.S.; Hong, J.; Telu, S.; Hines, C.S.; Pike, V.W.; Zarate, C.A.; et al. cAMP signaling in brain is decreased in unmedicated depressed patients and increased by treatment with a selective serotonin reuptake inhibitor. *Mol. Psychiatry* **2016**, *22*, 754-759. [CrossRef]

130. Anacker, C.; Zunszain, P.A.; Carvalho, L.A.; Pariante, C.M. The glucocorticoid receptor: Pivot of depression and of antidepressant treatment? *Psychoneuroendocrinology* **2011**, *36*, 415-425. [CrossRef]

131. Anacker, C.; Zunszain, P.A.; Cattaneo, A.; Carvalho, L.A.; Garabedian, M.J.; Thuret, S.; Price, J.; Pariante, C.M. Antidepressants increase human hippocampal neurogenesis by activating the glucocorticoid receptor. *Mol. Psychiatry* **2011**, *16*, 738-750. [CrossRef]

132. de Kloet, E.R.; Otte, C.; Kumsta, R.; Kok, L.; Hillegers, M.H.J.; Hasselmann, H.; Kliegel, D.; Joëls, M. Stress and Depression: A Crucial Role of the Mineralocorticoid Receptor. *J. Neuroendocr.* **2016**, *28*. [CrossRef] [PubMed]

133. Herzog, D.P.; Cuadrado, D.P.; Treccani, G.; Jene, T.; Opitz, V.; Hasch, A.; Lutz, B.; Lieb, K.; Sillaber, I.; van der Kooij, M.A.; et al. A distinct transcriptional signature of antidepressant response in hippocampal dentate gyrus granule cells. *Transl. Psychiatry* **2021**, *11*, 4. [CrossRef] [PubMed]

134. Daskalakis, N.P.; Meijer, O.C.; de Kloet, E.R. Mineralocorticoid receptor and glucocorticoid receptor work alone and together in cell-type-specific manner: Implications for resilience prediction and targeted therapy. *Neurobiol. Stress* **2022**, *18*, 100455. [CrossRef] [PubMed]

135. Bin Wang, B.; Xin, N.; Qian, X.; Zhai, L.; Miao, Z.; Yang, Y.; Li, S.; Sun, M.; Xu, X.; Li, X.-J. Ah1 regulates the nuclear translocation of glucocorticoid receptor to modulate stress response. *Transl. Psychiatry* **2021**, *11*, 188. [CrossRef]

136. Sagner, A.; Zhang, I.; Watson, T.; Lazaro, J.; Melchionda, M.; Briscoe, J. A shared transcriptional code orchestrates temporal patterning of the central nervous system. *PLoS Biol.* **2021**, *19*, e3001450. [CrossRef]

137. Orsetti, M.; Di Brisco, F.; Rinaldi, M.; Dallorto, D.; Ghi, P. Some molecular effectors of antidepressant action of quetiapine revealed by DNA microarray in the frontal cortex of anhedonic rats. *Pharmacogenet. Genom.* **2009**, *19*, 600-612. [CrossRef]

138. Barakat, A.K.; Scholl, C.; Steffens, M.; Brandenburg, K.; Ising, M.; Lucae, S.; Holsboer, F.; Laje, G.; Kalayda, G.V.; Jaehde, U.; et al. Citalopram-induced pathways regulation and tentative treatment-outcome-predicting biomarkers in lymphoblastoid cell lines from depression patients. *Transl. Psychiatry* **2020**, *10*, 210. [CrossRef]

139. Machado-Vieira, R.; Zarate, C.A., Jr.; Manji, H.K. Emerging Novel Treatments for Severe Mood Disorders Involving Cellular Plasticity Cascades. *Curr. Psychos. Ther. Rep.* **2006**, *4*, 181-190. [CrossRef]

140. Cheng, Y.; An, Q.; Wang, J.; Wang, Y.; Dong, J.; Yin, J. RasGRF1 participates in the protective effect of tanshinone IIA on depressive like behaviors of a chronic unpredictable mild stress induced mouse model. *Gene* **2020**, *754*, 144817. [CrossRef]

141. Dmitrzak-Weglarz, M.; Szczepankiewicz, A.; Rybakowski, J.; Kapelski, P.; Bilska, K.; Skibinska, M.; Reszka, E.; Lesicka, M.; Jablonska, E.; Wieczorek, E.; et al. Expression Biomarkers of Pharmacological Treatment Outcomes in Women with Unipolar and Bipolar Depression. *Pharmacopsychiatry* **2021**, *54*, 261-268. [CrossRef]

142. Pothula, S.; Kato, T.; Liu, R.-J.; Wu, M.; Gerhard, D.; Shinohara, R.; Sliby, A.-N.; Chowdhury, G.M.I.; Behar, K.L.; Sanacora, G.; et al. Cell-type specific modulation of NMDA receptors triggers antidepressant actions. *Mol. Psychiatry* **2020**, *26*, 5097-5111. [CrossRef] [PubMed]

143. Hare, B.D.; Shinohara, R.; Liu, R.J.; Pothula, S.; Dileone, R.J.; Duman, R.S. Optogenetic stimulation of medial prefrontal cortex Drd1 neurons produces rapid and long-lasting antidepressant effects. *Nat. Commun.* **2019**, *10*, 223. [CrossRef] [PubMed]

144. Belzeaux, R.; Azorin, J.-M.; Ibrahim, E.C. Monitoring candidate gene expression variations before, during and after a first major depressive episode in a 51-year-old man. *BMC Psychiatry* **2014**, *14*, 73. [CrossRef] [PubMed]
145. Liu, Z.; Zhu, F.; Yao, L.; Yang, C.; Xiao, L.; Zhu, J.; Wang, H.; Wang, G.; Liu, W.; Xiao, Z. PDLIM5 gene polymorphisms and short term antidepressant response in Chinese major depressive disorders. *Int. J. Clin. Exp. Med.* **2013**, *6*, 677–682. [PubMed]

146. Horiuchi, Y.; Ishikawa, M.; Kaito, N.; Iijima, Y.; Tanabe, Y.; Ishiguro, H.; Arinami, T. Experimental Evidence for the Involvement of PDLIM5 in Mood Disorders in Hetero Knockout Mice. *PLoS ONE* **2013**, *8*, e59320. [CrossRef]

147. Ribeiro, L.F.; Verpoort, B.; Nys, J.; Vennekens, K.M.; Wierda, K.D.; De Wit, J. SorCS1-mediated sorting in dendrites maintains neurexin axonal surface polarization required for synaptic function. *PLOS Biol.* **2019**, *17*, e3000466. [CrossRef]

148. Lauterbach, E.C. Psychotropic Drug Effects on Gene Transcriptomics Relevant to Alzheimer Disease. *Alzheimer Dis. Assoc. Disord.* **2012**, *26*, 1–7. [CrossRef]

149. Pechstein, A.; Shupliakov, O.; Haucke, V. Intersectin 1: A versatile actor in the synaptic vesicle cycle. *Biochem. Soc. Trans.* **2010**, *38*, 181–186. [CrossRef]

150. Gerosa, L.; Mazzoleni, S.; Rusconi, F.; Longaretti, A.; Lewerissa, E.; Pelucchi, S.; Murru, L.; Giannelli, S.G.; Broccoli, V.; Marcello, E.; et al. The epilepsy-associated protein PCDH19 undergoes NMDA receptor-dependent proteolytic cleavage and regulates the expression of immediate-early genes. *Cell Rep.* **2022**, *39*, 110857. [CrossRef]

151. von Wittgenstein, J.; Zheng, F.; Wittmann, M.-T.; Balta, E.-A.; Ferrazzi, F.; Schäffner, I.; Häberle, B.M.; Valero-Aracama, M.J.; Koehl, M.; Miranda, C.J.; et al. Sox11 is an Activity-Regulated Gene with Dentate-Gyrus-Specific Expression Upon General Neural Activation. *Cereb. Cortex* **2020**, *30*, 3731–3743. [CrossRef]

152. Komine, Y.; Takao, K.; Miyakawa, T.; Yamamori, T. Behavioral Abnormalities Observed in Zfhx2-Deficient Mice. *PLoS ONE* **2012**, *7*, e53114. [CrossRef] [PubMed]

153. Habib, A.M.; Matsuyama, A.; Okorokov, A.L.; Santana-Varela, S.; Bras, J.T.; Aloisi, A.M.; Emery, E.C.; Bogdanov, Y.D.; Follenfant, M.; Gossage, S.J.; et al. A novel human pain insensitivity disorder caused by a point mutation in ZFHX2. *Brain* **2018**, *141*, 365–376. [CrossRef]

154. Thippeswamy, H.; Davies, W. A new molecular risk pathway for postpartum mood disorders: Clues from steroid sulfatase-deficient individuals. *Arch. Women's Ment. Health* **2020**, *24*, 391–401. [CrossRef] [PubMed]

155. Niculescu, A.B.; Levey, D.; Phalen, P.; Le-Niculescu, H.; Dainton-Howard, H.; Jain, N.; Belanger, E.; James, A.; George, S.; Weber, H.; et al. Understanding and predicting suicidality using a combined genomic and clinical risk assessment approach. *Mol. Psychiatry* **2015**, *20*, 1266–1285. [CrossRef] [PubMed]

156. Niculescu, A.B.; Le-Niculescu, H.; Levey, D.F.; Phalen, P.L.; Dainton, H.L.; Roseberry, K.; Niculescu, E.M.; Niezer, J.O.; Williams, A.; Graham, D.L.; et al. Precision medicine for suicidality: From universality to subtypes and personalization. *Mol. Psychiatry* **2017**, *22*, 1250–1273. [CrossRef]

157. Ferreira, A.C.; Osório, F.D.L. Peripheral oxytocin concentrations in psychiatric disorders—A systematic review and methanalysis: Further evidence. *Prog. Neuro-Psychopharmacol. Biol. Psychiatry* **2022**, *117*, 110561. [CrossRef]

158. Hori, H.; Nakamura, S.; Yoshida, F.; Teraishi, T.; Sasayama, D.; Ota, M.; Hattori, K.; Kim, Y.; Higuchi, T.; Kunugi, H. Integrated profiling of phenotype and blood transcriptome for stress vulnerability and depression. *J. Psychiatr. Res.* **2018**, *104*, 202–210. [CrossRef]

159. Guilloux, J.-P.; Bassi, S.; Ding, Y.; Walsh, C.; Turecki, G.; Tseng, G.; Cyranowski, J.M.; Sibille, E. Testing the Predictive Value of Peripheral Gene Expression for Nonremission Following Citalopram Treatment for Major Depression. *Neuropsychopharmacology* **2015**, *40*, 701–710. [CrossRef]

160. Davis, S.; Meltzer, P.S. GEOquery: A bridge between the Gene Expression Omnibus (GEO) and BioConductor. *Bioinformatics* **2007**, *23*, 1846–1847. [CrossRef]

161. Carvalho, B.S.; Irizarry, R.A. A framework for oligonucleotide microarray preprocessing. *Bioinformatics* **2010**, *26*, 2363–2367. [CrossRef]

162. Ritchie, M.E.; Belinda, P.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. limma powers differential expression analyses for RNA-sequencing and microarray studies. *Nucleic Acids Res.* **2015**, *43*, e47. [CrossRef] [PubMed]

163. Joshi, N.A.; Fass, J.N. *Sickle: A Sliding-Window, Adaptive, Quality-Based Trimming Tool for FastQ Files*; Version 1.33; Github: San Francisco, CA, USA, 2011.

164. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. *Bioinformatics* **2013**, *29*, 15–21. [CrossRef] [PubMed]

165. Liao, Y.; Smyth, G.K.; Shi, W. feature Counts: An efficient general purpose program for assigning sequence reads to genomic features. *Bioinformatics* **2014**, *30*, 923–930. [CrossRef]

166. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. *Genome Biol.* **2014**, *15*, 550. [CrossRef] [PubMed]

167. Zhu, A.; Ibrahim, J.G.; Love, M.I. Heavy-tailed prior distributions for sequence count data: Removing the noise and preserving large differences. *Bioinformatics* **2018**, *35*, 2084–2092. [CrossRef]

168. Wickham, H. *ggplot2*, 1st ed.; Use R; Springer: New York, NY, USA, 2009; p. 213.

169. Prada, C.; Lima, D.; Nakaya, H. *MetaVolcanoR: Gene Expression Meta-Analysis Visualization Tool*; R Package Version 1.10.0; GitHub: San Francisco, CA, USA, 2022.

170. Sherman, B.T.; Hao, M.; Qiu, J.; Jiao, X.; Baseler, M.W.; Lane, H.C.; Imamichi, T.; Chang, W. DAVID: A web server for functional enrichment analysis and functional annotation of gene lists (2021 update). *Nucleic Acids Res.* **2022**, *50*, W216–W221. [CrossRef]
171. Kang, H.J.; Voleti, B.; Hajszan, T.; Rajkowska, G.; Stockmeier, C.A.; Licznerski, P.; Lepack, A.; Majik, M.; Jeong, L.S.; Banasr, M.; et al. Decreased expression of synapse-related genes and loss of synapses in major depressive disorder. *Nat. Med.* **2012**, *18*, 1413–1417. [CrossRef]
172. Seney, M.L.; Huo, Z.; Cahill, K.; French, L.; Puralewski, R.; Zhang, J.; Logan, R.W.; Tseng, G.; Lewis, D.A.; Sibille, E. Opposite Molecular Signatures of Depression in Men and Women. *Biol. Psychiatry* **2018**, *84*, 18–27. [CrossRef]
173. Pettai, K.; Milani, L.; Tammiste, A.; Võsa, U.; Kolde, R.; Eller, T.; Nutt, D.; Metspalu, A.; Maron, E. Whole-genome expression analysis reveals genes associated with treatment response to escitalopram in major depression. *Eur. Neuropsychopharmacol.* **2016**, *26*, 1475–1483. [CrossRef]
174. Duric, V.; Banasr, M.; Stockmeier, C.A.; Simen, A.A.; Newton, S.S.; Overholser, J.C.; Jurjus, G.J.; Dieter, L.; Duman, R.S. Altered expression of synapse and glutamate related genes in post-mortem hippocampus of depressed subjects. *Int. J. Neuropsychopharmacol.* **2012**, *16*, 69–82. [CrossRef]