Skip to main content

Early life and adult stress promote sex dependent changes in hypothalamic miRNAs and environmental enrichment prevents stress-induced miRNA and gene expression changes in rats



The hypothalamus plays a key role in the stress response. While early life stress (ELS) increases susceptibility to psychiatric disorders including major depressive disorder (MDD), acute stress during adulthood can also precipitate MDD after ELS.


Here, we tested the expression of miRNAs following ELS and susceptibility to depression-like behavior and whether sex or acute stress exacerbates this response. We also tested whether environmental enrichment (Enr) promotes early life and adult behavioral stress resilience and its effect on hypothalamic miRNA and gene expression. Following rat maternal separation (MS) as an ELS model, Enr from weaning through adulthood, and restraint (RS) as acute adult stress, we tested both animal behavior and miRNA expression in the hypothalamus. Target genes and their enrichment and ontology were analyzed using bioinformatic tools. Target gene expression changes were tested using qPCR, and miRNA promoter methylation was studied using methylated-DNA immunoprecipitation qPCR.


MS, Enr, RS, and sex altered hypothalamic miRNAs, including several previously reported in MS literature: miRs-29, − 124, − 132, − 144, − 504. Sex had a significant effect on the greatest number of miRNAs. Also, Enr reversed downregulation of miR-29b-1-5p and -301b-3p in MS. qPCR showed that MAPK6 and MMP19, targets of miR-301b-3p, were upregulated in MS and reversed by Enr. Additionally, miR-219a was hypermethylated in MS coinciding with decreased miR-219a expression.


This study found that sex plays a critical role in the hypothalamic miRNA response to both ELS and acute stress, with males expressing greater changes following postnatal stress. Moreover, enrichment significantly altered behavior as well as hypothalamic miRNA expression and their gene targets. Because of its role as the initiator of the autonomic stress response and connection to hedonic and motivational behavior, the hypothalamic miRNA landscape may significantly alter both the short and long-term behavioral response to stress.

Peer Review reports


Early life stress (ELS), namely abuse, neglect, and household dysfunction, is estimated to affect 64% of the US population [1]. This poses a major public health risk as ELS increases susceptibility to psychiatric disorders, including major depressive disorder (MDD). Not only does MDD carry a sizeable economic burden due to lost productivity [2], its symptoms are particularly debilitating. MDD onset has also been strongly associated with a precipitating stressful event, especially within the previous month [3]. Moreover, ELS and stress during adulthood can interact producing more severe or long-lasting symptoms than individually [4], also described as the two/three-hit hypothesis of stress susceptibility [5]. There are also significant sex differences in the response to ELS [6]. In the US, men and women experience similar levels of ELS [7], yet women are twice as likely to develop MDD [8]. There are only a limited number of studies that have examined sex-mediated depressive response to ELS [9].

Previous preclinical and clinical studies have implicated hypothalamic-pituitary-adrenal (HPA)-axis responsiveness in the depressive effects of ELS [10,11,12]. The HPA axis is the primary neurochemical stress response system whereby corticotropin releasing hormone (CRH) is secreted from the hypothalamus triggering a cascade of hormone release by the pituitary (adrenocorticotropin-releasing hormone-ACTH) and adrenal glands (corticosterone-CORT). Not only is the hypothalamus the initial site of this neuroendocrine response, but it also receives inputs from regions central to emotion processing like the amygdala [13] and frontal cortex [14]. Moreover, it plays a particularly important role in the symptomatology of MDD [15]. The hypothalamus receives sensory cues predicting reward (or lack thereof) and stimulates dopamine downstream [16]. In rodents, disruption of dopaminergic projections to the hypothalamus leads to reduced sucrose preference (i.e. more anhedonia-like behavior) [17]. Similarly, antagonizing hypothalamic GABA receptors increases anxiety-related measures in the elevated plus maze (EPM) such as time spent and number of entries into the closed arms [18] and directly injecting calcium channel inhibitors into the hypothalamus increases swimming behavior in the forced swim test (FST) indicated reduced behavioral despair [19]. These findings point to several neurochemical signaling mechanisms in the hypothalamus which might contribute to increased stress and depression-related behaviors. ELS has also been associated with disturbed HPA function including increased cortisol response to ACTH administration and increased cerebrospinal fluid CRH levels in MDD patients [20, 21]. While our understanding of the HPA stress response and its interaction with ELS has improved over the decades, currently available treatments for MDD are only modestly effective.

Recently, microRNAs (miRNAs) have been proposed as potential therapeutic or biomarker targets for MDD [22,23,24,25]. These small non-coding RNAs (~ 22 nucleotides) are synthesized in the cell nucleus as a hairpin loop and then exported into the cytoplasm where the loop is cleaved. This mature, single-stranded miRNA is loaded into the RNA-induced silencing complex (RISC) and targets the 3’UTR of genes with partially matching nucleotide sequences [26]. Canonically, miRNAs block the translation of mRNAs into proteins. Furthermore, because miRNAs are relatively short and can bind to mRNAs without perfect base pairing, they are able to target many different genes concurrently. Only a few studies to date have explored the effect of ELS on miRNAs (reviewed in [27]), one of which employed short duration maternal separation (MS) as a form of enrichment and found increases in miRs-488, − 144, and − 542-5p and decreases in miRs-421 and -376b-5 in hypothalamus [28]. MS has been well established as a rodent model of ELS associated with depression- and anxiety-like behavior [29]. MS protocols vary, but typically rodent pups are separated from their dam for 180 min daily for 12–21 days. Another study combined postnatal MS and restraint stress (RS) during adulthood to test if RS precipitated depressive behavior after MS [30]. Animals who experienced both MS and RS showed increased immobility in the forced swim test and decreased sucrose preference along with increased medial prefrontal cortex expression of miR-124 [30]. Increased miR-124 in the dentate gyrus after 90 min MS has also been reported [31]. While miRNAs serve a primary function to regulate gene expression, they are also regulated by epigenetic modifications such as methylation [32]. Only a few studies have elucidated miRNA methylation changes in psychiatric disorders. In adolescent patients with high risk for MDD, significant hypomethylation was found in miR-4646-3p promoter region [33]. Using chronic CORT administration to induce depression-like behavior in rats, our group found promoter hypomethylation of miR-124-3p [34]. No studies have explored miRNA promoter methylation in an ELS model.

There has been a growing shift toward non-pharmaceutical treatments in mental health. Clinical trials have increasingly pursued behavioral interventions while pharmaceutical-based trials have decreased from 43% in 2007 to 27% by 2018 [35]. Depressed individuals with a history of ELS may respond less to antidepressant drug therapy than patients with no ELS history [36]; however, very little is yet understood about how non-pharmaceutical methods contribute to healthy brain function. Environmental enrichment (Enr) has been found to improve behaviors in rodents related to anxiety [37, 38] and MDD [39] and lower blood CORT levels after RS [38]. Earlier work proposed that Enr reverses HPA dysregulation after ELS [40, 41].

The purpose of the study was: 1) to show whether ELS-induced depressive behavior is associated with miRNA expression changes in the hypothalamus and whether there is an interaction between MS, RS, and sex on miRNA expression; 2) to examine the effect of the interaction between RS, sex and Enr in MS animals on hypothalamic miRNA expression; 3) explore if these MS-induced miRNAs changes in the hypothalamus relate to gene expression changes and whether these miRNAs themselves might be regulated via methylation. To examine these, we assessed the behavioral and physiological consequences of MS and Enr as well as genome-wide changes in miRNA expression using next-generation sequencing. We also used bioinformatic tools to understand the potential functions associated with altered miRNAs and their targets. Finally, we tested the expression of stress-related miRNA gene targets and potential regulation of miRNAs by promoter region methylation.


Differential miRNA expression in MS, RS, and enriched animals

Animals were assigned to MS, RS, and Enr groups or relative controls. As shown in Fig. 1a MS animals were separated for 14 days and enrichment began immediately post-weaning; control animals were handled daily but not separated and were housed in standard conditions. Animals in RS groups were restrained for 7 days prior to behavioral testing and tissue collection; animals assigned to non-RS control groups were briefly handled daily. Finally, RNA isolated from the hypothalamus was sequenced to determine miRNA expression differences across these groups; animal groups for each set of comparisons are shown in Fig. 1b.

Fig. 1
figure 1

Behavior and physiology associated with maternal separation, restraint stress, sex, and enrichment. A schematic diagram (a) shows the timeline of behavioral testing relative to MS, Enr, and RS. MS was conducted for 180 min daily for the first 14 postnatal days. Enr was provided to animals following weaning until tissues were harvested at PND 95. b) Statistical comparisons of behavior and miRNA expression were conducted using 3-way ANOVA across MS, RS, and sex (n = 6 animals per group; first 4 columns of Fig. 1b). Another 3-way analysis was conducted across RS, Enr, and sex, only in MS animals (last 4 columns of Fig. 1b). Analysis of Enr as a treatment following MS was conducted using one-way ANOVA with control, MS, and MS + Enr as groups (indicated across the top of Fig. 1b; n = 12 animals per group including n = 6 RS and n = 6 Non-RS animals per group). c-k) Behavioral changes based on the interaction between MS, RS, and sex are shown. Importantly, MS reduced sucrose preference (e; F(1,48)=9.24, p = 0.004) regardless of RS or sex. Sex was a significant predictor for escape latency (g; F(1,48)=7.70, p = 0.008), EPM movement (h; F(1,48)=9.20, p = 0.004), and FST swim and climb scores (i; F(1,48)=4.19, p = 0.047; j;F(1,48)=8.61, p = 0.006). l-t) In MS animals, a 3-way ANOVA revealed increased EPM movement (n; F(1,48)=22.91, p < 0.000) and sucrose preference (p:(1,48)=16.05, p < 0.000) and decreased FST climbing (o; F(1,48)=11.92, p = 0.001) and adrenal weight (s; F(1,48)=5.23, p = 0.028) following Enr. RS significantly increased plasma Cort levels (t; F(1,48)=4.39, p < 0.05).One-way ANOVA comparison of control, MS, and MS + Enr groups (u-y) revealed the Enr recovered sucrose preference changes induced by MS (u; F(2,69)=12.049, p < .001) as well as recovering FST climbing changes in male animals (v; F(2, 69)=6.841, p < 0.005). Bars represent Mean ± SEM with males in blue and females in red. Dotted linear trend lines in f indicate an interaction. Significance: *p < 0.05. Abbreviations: PND-postnatal day, MS-maternal separation, RS-restraint stress, Enr-environmental enrichment

Behavior differences resulting from MS, RS, and sex

We tested for main and interaction effects on physiological and behavioral measures of stress (animals included in this analysis are noted in Fig. 1b by the label “MS x RS x sex”) and detected a significant 3-way interaction between MS, RS and sex on EPM closed arm time (F(1,48)=5.30, p = 0.027; Fig. 1c), however there was no significant effect on the EPM index score (p > 0.05; Fig. 1d). Follow-up pairwise t-tests showed that RS-only males spent more time in the closed arms compared to MS, MS + RS, and controls (t(22)= − 2.423, p = 0.056). There was a significant decrease in sucrose preference for all MS animals compared to all non-MS (F(1,48)=9.24, p = 0.004;Fig. 1e). RS significantly altered the effect of sex on body weight (F(1,48)=5.79, p = 0.021;Fig. 1f). There was also a significant main effect of sex on escape latency (F(1,48)=7.70, p = 0.008;Fig. 1g) with all males exhibiting higher escape latency than females. We found significantly decreased total movement in the EPM in all female animals compared to males (F(1,48)=9.20, p = 0.004;Fig. 1h). We also found that forced swim test (FST) swim scores were higher in female than male animals (F(1,48)=4.19, p = 0.047;Fig. 1i) but climb scores were higher in males than females (F(1,48)=8.61, p = 0.006;Fig. 1j); there were no differences in FST float score.

The interaction of MS, RS, and sex on corticosterone levels

We compared plasma CORT level across groups using an independent samples t-tests (Fig. 1k) but there were no significant differences (p > 0.05). In males, RS significantly increased CORT level compared to controls and MS-only males (t(22)= − 2.511, p < .05).

The effect of estrus on behavior and miRNAs

Because of the significant sex differences in behavior and concerns for the effect of female sex hormones on stress-related changes, we tested female animals for estrus phase via vaginal lavage. Representative microscope images of cell cytology and group behavior differences are shown in Supplementary Fig. 1. Because there were no significant behavior differences between the estrus phase groups further analyses did not include this variable. Detailed results are described in Supplementary Table 5.

Behavior differences resulting from RS, enrichment, and sex in MS rats

A RS x Enr x Sex ANOVA was conducted to identify the effect of these factors and their interaction on behavior. For this analysis, only MS animals were considered because we were primarily interested in exploring enrichment as a preventative measure as opposed to its use simply as enrichment in normal controls (animals included in this analysis are noted in Fig. 1b by the label “RS x Enr x sex”).

There were no significant effects of MS, RS, or Enr on EPM closed arm frequency (Fig. 1l) or the anxiety index (Fig. 1m). Enr significantly increased total movement in the EPM, regardless of RS (F(1,48)=22.91, p < 0.000;Fig. 1n). All Enr animals regardless of RS showed a significant reduction in climb scores compared to non-enriched animals (F(1,48)=11.92, p = 0.001;Fig. 1o). Enrichment significantly increased sucrose preference in MS animals regardless of whether they had received RS (F(1,48)=16.05, p < 0.000;Fig. 1p). There were no significant effects on escape latency (p > 0.05;Fig. 1q). As expected, there was a significant main effect of sex on body weight (F(1,48)=495.06, p < 0.000), but we also found that RS decreased body weight (F(1,48)=4.82, p = 0.034;Fig. 1r) compared to all animals without RS. Adrenal weight was significantly reduced by Enr regardless of RS (F(1,48)=5.23, p = 0.028;Fig. 1s). Finally, we found that RS increased plasma CORT levels (Fig. 1t; F(1,48)=4.39, p < 0.05).

miRNAs altered by the interaction of MS, RS, and sex

3-way ANOVA was used to detect significant miRNA expression differences associated with MS, RS, and sex (Table 1). The fold change column headers indicate the two comparison groups and direction of change (e.g. a fold change value of 1.33 for “MS:control” indicates 1.33-times fold change increase in expression in MS relative to control animals). Specific p values are listed alongside fold change in Table 1. Only miR-26a-3p was significantly affected by the 3-way interaction (p < 0.05). 5 miRNAs were affected by the interaction of MS and RS, 11 by sex and MS, and 46 by sex and RS (p < 0.05). Alone, sex significantly altered 24 miRNAs while MS affected 28 and RS affected 43 (p < 0.05). The miRNAs with the greatest fold change were miR-547-3p which was 2.45-fold downregulated in all females compared to all males (p < 0.05), miR-144-3p which was 2.3-fold downregulated in MS males compared to male controls (p < 0.05), and miR-23a-5p which was increased 1.97-fold in all females as compared to males (p < 0.05).

Table 1 miRNAs altered by the interaction between maternal separation, restraint stress, and sex

miRNAs altered by the interaction of RS, enrichment, and sex in maternally separated animals

miRNAs significantly altered by RS, Enr, sex or their interaction (based on 3-way ANOVA) in MS animals are reported in Table 2 along with their fold change and specific significance (p value). The 3-way interaction (RS, Enr, and sex) had a significant effect on 32 miRNAs. The largest number of miRNAs were affected by the interaction between Enr and sex, with 56 miRNAs significantly altered. There was a significant interaction between RS and Enr on miR-135b-3p expression (p < 0.05). miR-666-3p was significantly affected by the interaction between RS and sex (p < 0.05). We detected 2 miRNAs associated with RS main effect: miR-338-5p and − 341 (p < 0.05). There was a main effect of Enr on 39 miRNAs and sex affected expression of 14 miRNAs. The greatest fold change in expression was found for miR-547-3p with females showing a 2.51-fold downregulation compared to males (p < 0.05) regardless of RS or Enr. Next, miR-539-5p showed the next greatest fold change at 1.84-fold increased expression in Enr females compared to non-Enr females (p < 0.05).

Table 2 miRNAs altered by the interaction between restraint stress, environmental enrichment, and sex in maternally separated rats

miRNA-mediated hypothalamic signaling changes after MS and enrichment

Enrichment as prevention of depression- and anxiety-related behavior after MS

Our primary focus for this study was to elucidate hypothalamic miRNA changes resulting from MS and environmental enrichment. One-way ANOVA was used to compare controls, MS, and MS + Enr animals. Each group (controls, MS, and MS + Enr) consisted of 24 animals (12 male and 12 female) including RS animals; animals included in this analysis are noted in Fig. 1b as “one-way ANOVA”. There were significant group differences in sucrose preference (F(2,69)=12.049, p < .001); MS significantly decreased sucrose preference (p < 0.005) which was returned to normal levels by Enr (p < 0.001; Fig. 1u). Climb score (Fig. 1v) in the FST was significantly decreased by Enr as compared to both control (p < 0.05) and MS (p < 0.005) animals (F(2,69)=6.841, p < 0.005). In the EPM, Enr significantly increased total movement (Fig. 1w) compared to controls (p < 0.001) and MS (p < 0.001) animals (F(2,69)=17.077, p < 0.001). One-way ANOVA revealed significant differences (F(2,69)=3.445, p < 0.05) in EPM closed arm time (Fig. 1x), but follow-up pairwise t-tests only showed a nonsignificant reduction in MS (p = 0.077) and MS + Enr (p = 0.08) groups compared to controls. Enrichment significantly reduced adrenal weight (Fig. 1y) compared to both control (p < 0.005) and MS (p < 0.05) groups (F(2,69)=6.275, p < 0.005). Plasma CORT levels showed no significant group differences (F(2,69)=0.485, p > 0.05).

miRNA expression changes after MS and MS + Enrichment

A one-way ANOVA comparing control, MS, and MS + Enr miRNA counts per million (CPM) expression values showed that 29 miRNAs were significantly upregulated by MS and 21 were significantly downregulated. These changes include miRNAs which were only affected by MS as well as those only affected by MS + Enr or by both MS and MS + Enr. CPM-fold change and p values are shown in Table 3 and describe specific expression changes across each of the three groups. miRs-144-3p and − 206-5p exhibited the greatest fold change at 1.783- and 1.766-fold downregulation in MS animals compared to controls (p < 0.05). miR-29b-1-5p, −301b-3p, and 3065-5p showed a treatment response to Enr in that MS significantly downregulated their expression (p < 0.05) and Enr, at least partially, reversed this effect. Changes specific to male and female animals are shown in Supplementary Tables 1 and 2.

Table 3 miRNAs significantly altered by MS and Enrichment

We visualized the chromosomal location of each significantly altered miRNA in Fig. 2a. 12 miRNAs grouped together on chromosome 6 and each of these were upregulated by MS except for miRs-673-5p and − 341. 7 miRNAs grouped on chromosome 10.

Fig. 2
figure 2

miRNA, gene targets, and ontologies associated with maternal separation and enrichment. a two sizeable groups of miRNAs altered by MS and Enr localized onto chromosomes 6 and 10. Direction of fold change in MS relative to controls, MS + Enr relative to controls, and MS + Enr relatve to MS are represented as circles, diamonds, and triangles, respectively (blue indicates upregulation and green, downregulation). b miRNA gene target network shows that 8 significantly altered miRNAs target many stress-related genes and share many targets. c MetascapeGO revealed functions most associated with significant miRNAs’ gene targets (Top 3: rno05200-“pathways in cancer”, rno04014-“Ras signaling pathway”, and GO:0010720-“postive regulation of cell development”; more significant bars are darker colored). d Metascape also mapped similar ontologies based on common member genes (top map) and significance (bottom map). There is a high level of overlap in membership of these gene targets in Metascape ontologies. e Changes in miRNA gene target expression was tested using qPCR. Bars represent M ± SEM and significance were determined by independent samples t-test (n = 24 animals per group) and one-way ANOVA (for GRIN2B). GABRA1, MAPK6, and MMP19 were significantly upregulated by MS (t(46)> 2.035, p < 0.048) and MAPK6 and MMP19 were partially recovered by Enr. GRIN2B was the only gene affected by the interaction between sex and group (F(2, 72)=3.703, p = 0.030). f CPG islands were identified near miR-promoter region using the UCSC genome browser; the blue arrow shows miR-212 to the left and miR-132 on the right. The row marked in green shows a CPG island. g methylation in the 3’UTR of miR-219-1 was significantly increased by MS (t(30)=2.29, p = 0.029); due to limited tissue availability, methylation was tested in n = 8 animals per group (compared to n = 12 for gene and miRNA expression). For miR-207, methylation in two distinct regions of one CPG island was tested. Significance: *p < 0.05. Abbreviations: MS-maternal separation, C-controls, Enr-environmental enrichment, UTR-untranslated region, SEM-standard error of the mean

miRNA gene targets

We used Ingenuity Pathway Analysis (IPA; Hilden, Germany) to determine which miRNAs associated with MS and Enr shared common stress and depression-related gene targets and had the greatest number of these targets. Thus, miRNA-gene target maps were created based on significantly altered miRNAs and their gene targets. To determine hub miRNAs (i.e. those with the most targets and the most shared common targets), we filtered out miRNAs with 10 or fewer gene connections. 3418 gene targets were identified. miRs-132-5p, − 132-3p, −449a-5p, 30e-5p, − 338-3p, −301b-3p, − 144-3p, and let-7 g-3p were hub miRNAs forming a highly integrated network with many shared gene targets (Fig. 2b). Gene ontology (GO) analysis (Fig. 2c) confirmed the involvement of these miRNA gene targets in MAPK signaling as well as Ras signaling, ion binding, and neuron part ontologies. The top ontologies for GO analysis are further described in Supplementary File 1. Hierarchical clustering of the GO terms in Metascape revealed that these pathways and functions were highly overlapped (Fig. 2d). Stress-related gene targets were selected for expression based on a literature search for anxiety, depression, stress, and related molecular pathways (Supplementary Table 3); we also included those genes with the greatest number of targeting miRNAs. 2-way ANOVA revealed that only GRIN2B expression was significantly altered (F(2,72)=3.703, p = 0.030) by the interaction between sex and group (control, MS, and MS + Enr) with only females showing increased expression following Enr, though not in controls or MS animals (Fig. 2e). Using pairwise t-tests, we found significant upregulation of GABRA1 (t(46)=2.035, p = 0.048), MAPK6 (t(46)=2.307, p = 0.026), and MMP19 (t(46)=2.156, p = 0.036) in MS (n = 24) animals compared to controls (n = 24) (Fig. 2e). Several genes were upregulated in MS + Enr (n = 24) animals compared to controls (n = 24), including CALM1 (t(46)=4.307, p < 0.000), ESR1 (t(46)=2.520, p = 0.015), FOXO1 (t(46)=2.137, p = 0.038), GABRA1 (t(46)=2.384, p = 0.021), IL10RB (t(46)=2.271, p = 0.028), KALRN (t(46)=3.634, p = 0.001), PIK3CB (t(46)=3.077, p = 0.004), TBP (t(46)=2.307, p = 0.026), TRAF6 (t(46)=2.320, p = 0.025), TSC22D3 (t(46)=2.313, p = 0.025), and WNT2B (t(46)=2.697, p = 0.010). Only two genes showed expression changes following MS that returned to the level of controls following Enr: MAPK6 and MMP19.

Exploring methylation as a miRNA regulatory mechanism

Of the 50 miRNAs significantly altered by MS, we found 5 miRNAs (miRs-219a, − 207, − 132, − 212, and let-7i) with CPG islands within 1000 kb of their promoter region. miRs-132 and -212 localized closely on chromosome 10 and shared a CPG island in their promoter regions (Fig. 2f). Pairwise t-test revealed increased methylation at the promoter region of miR-219a in MS (2.61-fold increase, t(30)=2.29, p = 0.029; n = 4 male and n = 4 female per group) but not MS + Enr animals (Fig. 2g). There were no significant group differences in methylation for let-7i, miR-207, or miR-132/212 (p > 0.05).


In this study, we present the first genome-wide profiling of miRNAs in the hypothalamus following MS and RS. We also parsed sex and estrus phase differences in miRNA expression resulting from these stress paradigms. Finally, we show miRNA changes underlying the use of environmental enrichment as a preventative for ELS-induced depressive-behavior. We found that 29 miRNAs were upregulated and 21 were downregulated by MS. Out of 21 downregulated miRNAs, 3 miRNAs (miR-29b-1-5p, −301b-3p, and − 3065-5p) showed expression levels similar to controls in MS + Enr animals. Chromosomal localization revealed two large groups of MS-induced miRNAs on chromosome 6 and chromosome 10. Based on significantly altered miRNAs in MS, we detected over 3418 miRNA gene targets. qPCR-based expression analysis showed that GABRA1, MAPK6, and MMP19 were significantly increased in MS animals compared to controls. Interestingly, MAPK6 and MMP19 expression were reversed by Enr in MS animals. miRNA-gene target mapping revealed several miRNA regulatory hubs including miR-301b-3p, − 132-3p, − 132-5p, −449a-5p, −30e-5p, − 338-3p, − 144-3p, and let-7 g-3p. GO analysis confirmed that these miRNA gene targets were significantly involved in MAPK signaling. Finally, we found that a CPG island near the promoter of miR-219 was hypermethylated in MS but not significantly reversed by Enr.

MDD onset is usually preceded by a recent stressful event [3]. We did not find that RS precipitated a depression-like phenotype in MS animals. In female MS animals, RS increased escape latency but still this did not reach levels associated with learned helplessness (> 20 s latency [42];). Alternatively, each group variable (MS and RS) was associated with unique behavioral profiles. MS significantly decreased sucrose preference indicating increased anhedonia whereas RS decreased total movement in the EPM. However, RS-only males spent increased time in the closed arms of the EPM which has been mostly associated with increased anxiety [43]. RS animals also exhibited increased CORT levels compared to controls and MS without RS. Studies on ELS often report altered stress reactivity meaning that CORT may be more elevated following an acute stressor in individuals who experience ELS compared to those who have no ELS history [44]. Our data potentially support this, but it also shows that RS was sufficient to increase CORT levels whereas MS was not. The majority of behavioral differences were across sex. Males showed increased escape latency, increased anxiety index, increased EPM total movement, and increased climbing in the FST compared to females. Females showed increased swimming in the FST. There have been mixed reports regarding behavioral profiles following MS [45, 46], but sex differences are well documented [47], especially in the FST [48].

Only one miRNA, miR-26a-3p was significantly associated with the interaction between MS, RS, and sex. In this interaction, MS males showed increased expression of miR-26a compared to control males, while MS did not affect its expression in females. However, in RS animals, MS only increased miR-26a expression in females. miR-26a targets HTR1A, one of the serotonin receptors and is upregulated by antidepressants, fluoxetine and reboxetine [49]; however, this was only tested in males. A few other studies have shown significant changes in miR-26 following stress including increased expression in mouse prefrontal cortex [50] as well as in male rats prefrontal cortex who experienced repeated and ancestral stress [51]. Sex differences following ELS, especially stress sensitivity in males prior to puberty, may result from interactions between stress neurobiology and peripubertal sex hormone changes [52]. Future studies should systematically elucidate sex differences in how timing of stress affect miR-26, among others, and its gene targets such as HTR1A. miR-3593 showed the greatest fold change in RS vs. non-RS animals. Cattaneo, Cattane [53] found that hippocampal miR-3593 was downregulated in a rodent prenatal ELS model, but there are no reports in RS or other acute stress. In MS animals, we tested for an interaction between RS, Enr, and sex. In contrast to the 3-way interaction between MS, RS, and sex, there were several miRNAs significantly affected by the 3-way interaction between RS, Enr, and sex in MS animals. Interestingly, Enr in RS animals caused the opposite fold change direction as compared to Enr in non-RS animals. For example, Enr alone was associated with a 1.54-fold decrease in miR-702-3p expression compared to controls; Enr + RS males showed a 1.24-fold increase in miR-702-3p. This consistently occurred for almost every miRNA. This finding suggests that Enr and RS exhibit opposing miRNA profiles across many miRNAs.

Two large groups of miRNAs associated with MS were colocalized on chromosome 6 and 10. Several significant miRNAs within close proximity could indicate regulation by a common mechanism such as competitive endogenous RNAs [54]. Additionally, when miRNA are closely colocalized they may be transcribed in tandem as a cluster which may indicate similar or dependent functions as well as common regulatory elements [55]. Of the tested gene targets, CALM1 is most closely localized near (~ 9500 kilobases) the miRNA cluster on chromosome 6. Another gene, Delta-like noncanonical notch ligand 1 (DLK1), was identified extremely close (~ 100 kilobases) to the chromosome 6 cluster. DLK1 interacts with NOTCH signaling pathways and is especially involved in pituitary gland development [56]. Furthermore, after adolescence, DLK1 is expressed solely in neuroendocrine tissues in adults. DLK1 and NOTCH signaling changes as mediated by miRNAs may an important avenue for future stress and depression-related studies.

There are previous reports which suggest that estrus phase significantly alters the behavioral response to stress [57]. Our findings showed that only the FST was affected by estrus timing. Proestrus/diestrus (p/d) animals showed decreased swimming, yet increased climbing compared to estrus/metestrus (e/m) animals. This could be attributed to the increased ratio of MS to control animals in p/d group; however, no other behaviors showed significant differences. Kokras, Antoniou [48] also found that animals in the proestrus and diestrus phase exhibited greater climbing duration following treatment with sertraline compared to animals in estrus. A total of 6 miRNAs showed differential expression between the p/d and e/m groups. Of these, miRs-362-3p, − 374-5p, and − 384-5p are localized on the X chromosome. MiR-384-5p has been shown to regulate SOX9, a member of the NFkB pathway [58]. SOX 9 is also an embryonic transcription factor which responds to estrogen during sex determination [59]. Our findings, in general, do not support a particularly robust contribution of estrus timing to stress susceptibility.

In order to identify miRNAs and their regulatory pathways specific to MS and Enr, we compared miRNA expression between control (non-MS) animals, MS animals without Enr, and MS + Enr. While there were a greater number of upregulated miRNAs in MS compared to controls, downregulated miRNAs exhibited greater fold change differences between MS and controls. miRs-144-3p, − 206-5p, −29b-1-5p, and -301b-3p showed the greatest decrease in expression. In an augmented maternal care model using MS, Vogel Ciernia, Laufer [28] reported increased hypothalamic miR-144 expression. Since this model produces increased maternal care behaviors, this could be considered a form of enrichment. However, in comparison, our MS + Enr animals showed similar expression of miR-144 as the MS group. Several miR-29 family miRNAs have been reported in both rodent and human ELS studies. Uchida, Hara [30] found an increase in miR-29a following MS. On the other hand, Cattane, Mora [60] found increased miR-29b-3p and -29c-3p in healthy individuals with a history of ELS. In humans, at least, resilient individuals show opposite expression changes compared to our MS animals. Moreover, our findings show that Enr dampened the effect of MS on miR-29b-1-5p expression. miRNA gene target network mapping revealed several miRNA hubs for potential gene regulation. In another study of 180-min MS, miR-132 was upregulated compared to control [30]. We found similar increases in expression in both MS and MS + Enr groups. In MDD patients, plasma miR-132 expression is elevated and is also reduced by escitalopram treatment [61]. Let-7 g-3p was another hub miRNA that was significantly downregulated by MS. The let-7 family of miRNAs has been strongly implicated in various aspects of MDD pathophysiology [62]. Between various hub miRNAs, several common gene targets were found including PTEN, MAPK6, CALM1, WNT2B, and GRIN2B.

From the list of most commonly target genes, we selected several genes for qCPR-based validation. GABRA1 (GABA receptor subunit alpha-1) showed increased expression in both MS and MS + Enr groups. MAPK6 (mitogen-activated protein kinase 6) and MMP19 (matrix metalloproteinase-19) showed increased expression in MS, but MS + Enr groups showed expression reversal. One study found that GABRA1 was downregulated in postmortem brain of MDD patients who died by suicide [63]. Our findings only show increased expression of WNT2B in MS + Enr animals compared to controls, though not in MS animals without Enr. Wnt family proteins act on the ß-catenin pathway and may play a role in mood disorder pathophysiology [64]. Our group recently found decreased WNT5B expression associated with increased miR-128 in rodents with learned helplessness [65]. Another study found that neonatal MS reduced myelination in the prefrontal cortex via Wnt signaling changes [66].. MAPK6 encodes ERK3 and is widely expressed in the brain. ERK3 plays an important role in neonatal growth [67]. MAPK pathways have also been found to regulate some MMPs [68]. Wojtowicz and Mozrzymas [69] found that FN-439, a synthetic compound which blocks active MMP sites via chelation of Zn2+, led to a loss of long-term potentiation in CA1 of the hippocampus. However, there is currently little research specific to MMP19.

All 4 of these genes are targeted by miR-301b-3p. Following the canonical pattern of expression between miRNAs and their targets, miR-301b-3p was significantly downregulated while all 4 genes were upregulated. Furthermore, MS + Enr animals showed miR-301b expression more similar controls—the opposite as seen in MAPK6 and MMP19 expression. miR-144 was also a targeting miRNA for both MMP19 and MAPK6. This is one of the first studies to support coregulation of these pathways by miRNAs. Furthermore, they may respond to Enr as a preventative for depressive behavior, particularly anhedonia. In our study, Enr was associated with significant improvement in several behavior. Sucrose preference levels in MS + Enr animals were closer to that of control animals, showing the potential for Enr to protect against the depressogenic effects of MS. Enr animals also showed significantly increased total movement in the EPM indicating a reduced anxiety phenotype. Reduced adrenal weight in Enr animals suggests that Enr reduced engagement of the HPA stress axis, although CORT levels were unchanged.

We identified CPG islands near 4 miRNA promoter regions. miR-219a showed a significant increase in relative methylation following MS. Following the canonical relationship between methylation and decreased gene transcription [70], miR-219 was significantly downregulated in both MS and MS + Enr cohorts. This is the first study to implicate methylation as a potential mechanism of miRNA regulation following MS. TGFBR2 and ESR1 (Estrogen receptor 1) are both gene targets of miR-219a; however, neither showed significant upregulation following MS. A study of environmental enrichment in rats found increased expression of miR-219 in serum exosomes [71]. This study also transfused the serum of young, enriched animals into aging animals and found evidence of increased oligodendrocyte development. Further research is needed to parse the interaction between miRNA expression, miRNA methylation, estrogens, and environmental enrichment.


Both ELS and acute stress have widespread effects on miRNA expression which may mediate changes in stress-related behaviors. We found that sex significantly altered miRNA expression in response to stress, often with male animals showing more robust and extensive changes. This effect may be due to pre-pubertal timing of MS and highlights the need for more detailed characterization of the effect of ELS timing on behavior and neurobiology. Furthermore, this study identified several hypothalamic miRNAs of interest with regards to ELS. It is important for future studies to test whether direct manipulation of these candidate miRNAs lead to or even prevent later stress susceptibility and whether these behavioral changes are immediate and/or long-lasting. Enrichment had a particularly strong effect on animal behavior and miRNA expression and even reversed some of the effects of MS especially in the MAPK signaling pathway. These findings highlight enrichment as a noninvasive means of altering miRNAs which could prove useful in treating or preventing MDD onset. Finally, increased methylation may mediate some of the changes in miRNA expression resulting from ELS. Future studies utilizing next-generation sequencing will be crucial to detecting methylation of miRNAs with better resolution.



The study was approved by the Institutional Animal Care and Use Committee at the University of Alabama at Birmingham and was performed in accordance with relevant regulations, including ARRIVE guidelines. Animals were housed under standard care conditions (ad libitum food and water, 27 °C, 12-h light-dark cycle) for the duration of the experiment. Each group included 6 animals. The expanded methods are described in Supplementary File 1.

Maternal separation

An overview of the experiment is shown in Fig. 1a. Each litter of Holtzman rats (Envigo, Indianapolis, IN, USA) was randomly assigned to either the control (non-MS) or MS group. Non-MS controls were handled for 5 min from PND 1–14. MS pups were separated from the dam and housed individually on a 33 °C heating pad for 180 min each morning until PND 14. Pups were weaned from the dam on PND 21.

Restraint stress

Animals from both control and MS groups were randomly assigned to RS or control groups (non-RS). At PND 80, RS animals were placed in 20 cm restrainer tubes for 120 min each day for 7 days. Non-RS controls were handled but not restrained.

Environmental enrichment

A subset of MS and MS + RS animals was randomly assigned to receive Enr from PND 21 to 90. Enrichment included colored toys, tubes, shreddable cotton and paper objects, and manzanita wood which were rotated weekly to maintain novelty. Non-enriched animals were housed conventionally.

Animal behavior

Sucrose preference test

Sucrose preference was tested immediately following RS to assess anhedonia as previously described [72]. Briefly, on day 1, the animals were given sucrose1% (w/v). After acclimation on day 2, the animals were given access to both regular water and sucrose. Then the animals fasted for 24 h. On day 4, the animals were housed individually and given a premeasured 500 mL bottle of 1% sucrose and regular water. After 8 h, each of the bottles was measured for consumption and sucrose preference was calculated as previously described [72].

Elevated plus maze

Animals were placed individually in the center of a raised plus-shaped platform (50 × 50 cm) with two open and two walled arms (15 cm tall, open roof). For 5 min, the animals were recorded using Noldus Ethovision XT 11.5. Open and closed arm time and frequency were recorded and anxiety index was calculated as published earlier [73].

Forced swim test

Animals were acclimated for 15-min to an acrylic cylinder (28 cm diameter × 46 cm tall) filled to 25 cm with room temperature water. 24 h later, each animal was recorded in the same swim condition for 6 min and was returned to their home cage. A rater blinded to the conditions scored the videos using Kinoscope [74].

Shuttle escape test

As reported previously [72], escape latency was tested using a two-chamber shuttlebox with an electrified grid floor (70 cm × 20 cm × 20 cm, Med Associates, IN, USA). For the initial 5 trials, foot shocks were delivered on a variable interval schedule (0.6 mA, 60 s average interval) and terminated when the animal crossed into the opposite chamber. For the next 25 trials, the foot shocks terminated after the animal crossed into the opposite chamber and back. Escape latency was recorded by a PC connected to the shuttle box and shock generator.

Vaginal cytology, tissue collection, and RNA isolation

The morning following escape testing, prior to tissue collection, female animals were tested for estrous phase by vaginal lavage. A blinded-rater assessed the samples according to cell morphology as previously described [75, 76]. The animals were anesthetized with isofluorane and blood was collected via cardiac puncture. The brain and adrenal glands were dissected, flash frozen, and stored at − 80 °C until further use. RNA was extracted from hypothalamus using TRIzol (Invitrogen, NY, USA) as described in Roy, Dunbar [34]. RNA concentration and quality were tested using Nanodrop 2000C Spectrophotometer; 260/280 nm > 1.7 was considered pure.

Corticosterone, estradiol, and progesterone ELISAs

CORT, Estradiol, and Progesterone were quantified in platelet-free plasma using enzyme-linked immunosorbent assay (ELISA) (Enzo Life Sciences, NY, USA).

miRNA sequencing

RNA sequencing libraries were prepared using the Qiaseq miRNA library kit (Qiagen, Hilden, Germany) and were sequenced on a NextSeq 500 (Illumina, CA, USA). Data were extracted using Qiagen Gene Globe and CPM was calculated using edgeR.

Statistical approach

A 2 × 2 × 2 ANOVA was conducted in R to examine behavior and miRNA expression differences between MS and RS groups as well as sex. We also conducted a 2 × 2 × 2 ANOVA to examine the interaction of RS, sex, and Enr within MS animals on behavior and miRNA expression. Figure 1b is a visual representation of the animals/groups which were used in each analysis. We used a 1-way ANOVA to compare control, MS, and MS + Enr group behavior and miRNA expression. Subsequent bioinformatic analysis and follow-up gene expression and methylation studies were based on miRNA expression changes found when comparing control, MS, and MS + enrichment groups. To determine if sex played a role in MS and Enr-related gene target expression, we conducted a 2-way ANOVA. Only 1 gene showed a significant interaction, so the remaining gene expression comparisons were based on pairwise t-tests across group (control, MS, or MS + Enr). Group differences in qPCR-based methylation were tested using pairwise t-tests.

Bioinformatic analysis

Chromosomal localization was used to identify the chromosomal loci for significantly altered miRNAs. The Phenogram application (Ritchie Lab, University of Pennsylvania, USA) was used to plot the location of miRNAs across the rat karyotype.

Gene target prediction, miRNA-gene target networks, and gene ontology

A primary aim was to determine how MS leads to depression-like behavior via miRNA-based gene regulation. Thus, IPA was used to identify validated and highly predicted gene targets of miRNAs which were significantly altered by MS. This gene-target list was filtered based on each gene’s membership in canonical stress-related pathways. Futhermore, in order to identify a set of miRNAs with the greatest potential impact across these pathways only those pathways with > 16 gene targets were included. Because miRNAs can target many different genes and each gene may also be targeted by many miRNAs [26], pathways targeted by fewer miRNAs and with fewer compatible gene targets are less likely to be functionally affected by changes in an individual miRNA’s expression; it is possible that changes in a single miRNA may be negated by the activity of other miRNAs or even other regulatory mechanism. However, when several miRNAs targeting many genes within the same molecular pathway are altered, it is less likely that this cumulative effect will be countered by another miRNA or other mechanism. The identified miRNA targets were narrowed to 99 genes involved in Axonal Guidance Signaling (27 genes), Glucocorticoid Receptor Signaling (27 genes), Neuroinflammation Signaling (26 genes), Synaptogenesis Signaling (23 genes), Estrogen Receptor Signaling (23 genes), Protein Kinase A Signaling (19 genes), AMPK Signaling (18 genes), and ERK/MAPK Signaling (17 genes). miRNA-gene target networks were created using these 99 genes and significantly altered miRNAs based on one-way ANOVA. Among these miRNA-gene targeting relationships, miRNAs with more than 10 gene targets were identified as gene regulatory hubs and were also visualized in IPA as a miRNA-gene network. ShinyGO v0.61 [77] was used to explore GO associated with miRNA target genes. Metascape [78] was also used to confirm our GO analysis and to cluster ontologies with similar gene components.

Gene targets for qPCR follow-up were selected from the miRNA hub-gene target network based on their number of miRNA targets. Gene targets with the highest fold change as well as those which showed reversal of expression changes after Enr were also selected for qPCR-based expression testing.

CDNA synthesis and gene target expression by qPCR

1 μg of RNA was reverse transcribed using an oligo dT priming method to synthesize first strand complimentary DNA (cDNA). Relative gene abundance was tested using qPCR and BrightGreen chemistry (Applied Biological Material, Canada). Primer sequences for the selected genes are listed in Supplementary Table 4. Livak’s ΔΔ CT [79] was method to quantify group differences. The data was normalized to the geometric mean of GAPDH, ß-actin, and 18srRNA. Because there was little variability in the housekeeping genes, each sample was tested individually for each gene of interest.

Methylated DNA immunoprecipitation (MeDIP) and miRNA promoter region methylation by qPCR

Using the UCSC rodent genome browser tool, we searched for CPG islands near the promoter region of significantly altered miRNAs. We only found CPG islands upstream of 5 miRNAs and designed primers (Supplementary Table 4) to target these regions using Primer3 (v0.04.0; Untergasser, Cutcutache [80]). Genomic DNA (gDNA) was isolated from hypothalamus using phenol: chloroform: isoamylalcohol (25:24:1 V/V) chemistry. gDNA was sheared and immunoprecipitated with 5-methyl cytosine antibody (Zymo Research, CA). BrighGreen-based qPCR was used to quantify relative methylation enrichment in the immunoprecipitated DNA using the designed primers. Because of limited tissue availability, miRNA promoter methylation was tested in 48 of 72 original samples (n = 4 per group, including males and females).

Availability of data and materials

The datasets generated and analyzed during the current study are available in the NCBI Sequence Read Archive (SRA) repository, at and accession #: PRJNA728805.



early life stress


major depressive disorder


maternal separation


restraint stress


hypothalamic pituitary adrenal


corticotropin releasing hormone


adrenocorticotropin releasing hormone



miRNA, miR:



RNA-induced silencing complex


environmental enrichment


elevated plus maze


forced swim test






counts per million


gene ontology


complimentary DNA


methylated DNA immunoprecipitation


  1. Middlebrooks JS, Audage NC. The effects of childhood stress on health across the lifespan. Atlanta: Centers for Disease Control and Prevention, National Center for Injury Prevention and Control. 2008. p. 3-16.

  2. Greenberg PE, Fournier A-A, Sisitsky T, Pike CT, Kessler RC. The economic burden of adults with major depressive disorder in the United States (2005 and 2010). J Clin Psychiatry. 2015;76(2):155–62.

    Article  PubMed  Google Scholar 

  3. Kendler KS, Karkowski LM, Prescott CA. Causal relationship between stressful life events and the onset of major depression. Am J Psychiatry. 1999;156(6):837–41.

    Article  CAS  PubMed  Google Scholar 

  4. Seok BJ, Jeon S, Lee J, Cho SJ, Lee YJ, Kim SJ. Effects of early trauma and recent stressors on depression, anxiety, and anger. Front Psychiatry. 2020;11:744.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Daskalakis NP, Bagot RC, Parker KJ, Vinkers CH, de Kloet ER. The three-hit concept of vulnerability and resilience: toward understanding adaptation to early-life adversity outcome. Psychoneuroendocrinology. 2013;38(9):1858–73.

    Article  PubMed  PubMed Central  Google Scholar 

  6. McEwen BS, Gray JD, Nasca C. 60 YEARS OF NEUROENDOCRINOLOGY: redefining neuroendocrinology: stress, sex and cognitive and emotional regulation. J Endocrinol. 2015;226(2):T67–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. U.S. Department of Health & Human Services, Administration for Children and Families, Administration on Children, Youth and Families, Children’s Bureau. Child maltreatment 2016. 2018. Available from

  8. CfBHSa Q. 2016 National Survey on drug use and health: methodological summary and definitions. Substance Abuse and Mental Health Services Administration: Rockville, MD; 2017.

    Google Scholar 

  9. Coiro P, Pollak DD. Sex and gender bias in the experimental neurosciences: the case of the maternal immune activation model. Transl Psychiatry. 2019;9(1):90.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Heim C, Newport DJ, Heit S, Graham YP, Wilcox M, Bonsall R, et al. Pituitary-adrenal and autonomic responses to stress in women after sexual and physical abuse in childhood. Jama. 2000;284(5):592–7.

    Article  CAS  PubMed  Google Scholar 

  11. Vazquez DM, Lopez JF, Van Hoers H, Watson SJ, Levine S. Maternal deprivation regulates serotonin 1A and 2A receptors in the infant rat. Brain Res. 2000;855(1):76–82.

    Article  CAS  PubMed  Google Scholar 

  12. Horii-Hayashi N, Sasagawa T, Matsunaga W, Matsusue Y, Azuma C, Nishi M. Developmental changes in desensitisation of c-Fos expression induced by repeated maternal separation in pre-weaned mice. J Neuroendocrinol. 2013;25(2):158–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Price JL. Comparative aspects of amygdala connectivity. Ann N Y Acad Sci. 2003;985(1):50–8.

    Article  PubMed  Google Scholar 

  14. Öngür D, Price JL. The Organization of Networks within the orbital and medial prefrontal cortex of rats. Cereb Cortex: Monkeys and Humans; 2000.

    Google Scholar 

  15. Price JL, Drevets WC. Neurocircuitry of mood disorders. Neuropsychopharmacology : Official Publication of the American College of Neuropsychopharmacology. 2010;35(1):192–216.

    Article  Google Scholar 

  16. Sakurai T. The role of orexin in motivated behaviours. Nat Rev Neurosci. 2014;15(11):719–31.

    Article  CAS  PubMed  Google Scholar 

  17. Martinez-Hernandez J, Lanuza E, Martinez-Garcia F. Selective dopaminergic lesions of the ventral tegmental area impair preference for sucrose but not for male sexual pheromones in female mice. Eur J Neurosci. 2006;24(3):885–93.

    Article  PubMed  Google Scholar 

  18. Shekhar A. GABA receptors in the region of the dorsomedial hypothalamus of rats regulate anxiety in the elevated plus-maze test. I Behav Meas Brain Res. 1993;627(1):9–16.

    Article  CAS  Google Scholar 

  19. Bidzinski A, Jankowska E, Pucilowski O. Antidepressant-like action of nicardipine, verapamil and hemicholinium-3 injected into the anterior hypothalamus in the rat forced swim test. Pharmacol Biochem Behav. 1990;36(4):795–8.

    Article  CAS  PubMed  Google Scholar 

  20. Heim C, Nemeroff CB. The role of childhood trauma in the neurobiology of mood and anxiety disorders: preclinical and clinical studies. Biol Psychiatry. 2001;49(12):1023–39.

    Article  CAS  PubMed  Google Scholar 

  21. Heim C, Newport DJ, Mletzko T, Miller AH, Nemeroff CB. The link between childhood trauma and depression: insights from HPA axis studies in humans. Psychoneuroendocrinology. 2008;33(6):693–710.

    Article  CAS  PubMed  Google Scholar 

  22. Dwivedi Y. microRNA-124: a putative therapeutic target and biomarker for major depression. Expert Opin Ther Targets. 2017;21(7):653–6.

    Article  CAS  PubMed  Google Scholar 

  23. O'Connor RM, Dinan TG, Cryan JF. Little things on which happiness depends: microRNAs as novel therapeutic targets for the treatment of anxiety and depression. Mol Psychiatry. 2012;17(4):359–76.

    Article  CAS  PubMed  Google Scholar 

  24. Dwivedi Y. Emerging role of microRNAs in major depressive disorder: diagnosis and therapeutic implications. Dialogues Clin Neurosci. 2014;16(1):43–61.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Roy B, Yoshino Y, Allen L, Prall K, Schell G, Dwivedi Y. Exploiting circulating MicroRNAs as biomarkers in psychiatric disorders. Mol Diagn Ther. 2020;24(3):279–98.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Bartel DP. Metazoan MicroRNAs. Cell. 2018;173(1):20–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Allen L, Dwivedi Y. MicroRNA mediators of early life stress vulnerability to depression and suicidal behavior. Mol Psychiatry. 2020;25(2):308–20.

    Article  CAS  PubMed  Google Scholar 

  28. Vogel Ciernia A, Laufer BI, Dunaway KW, Mordaunt CE, Coulson RL, Totah TS, et al. Experience-dependent neuroplasticity of the developing hypothalamus: integrative epigenomic approaches. Epigenetics. 2018;13(3):318–30.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Vetulani J. Early maternal separation: a rodent model of depression and a prevailing human condition. Pharmacol Rep. 2013;65(6):1451–61.

    Article  PubMed  Google Scholar 

  30. Uchida S, Hara K, Kobayashi A, Funato H, Hobara T, Otsuki K, et al. Early life stress enhances behavioral vulnerability to stress through the activation of REST4-mediated gene transcription in the medial prefrontal cortex of rodents. J Neurosci. 2010;30(45):15007–18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Bahi A. Sustained lentiviral-mediated overexpression of microRNA124a in the dentate gyrus exacerbates anxiety- and autism-like behaviors associated with neonatal isolation in rats. Behav Brain Res. 2016;311:298–308.

    Article  CAS  PubMed  Google Scholar 

  32. Parodi F, Carosio R, Ragusa M, Di Pietro C, Maugeri M, Barbagallo D, et al. Epigenetic dysregulation in neuroblastoma: a tale of miRNAs and DNA methylation. Biochim Biophys Acta. 2016;1859(12):1502–14.

    Article  CAS  PubMed  Google Scholar 

  33. Bostrom AE, Ciuculete DM, Attwood M, Krattinger R, Nikontovic L, Titova OE, et al. A MIR4646 associated methylation locus is hypomethylated in adolescent depression. J Affect Disord. 2017;220:117–28.

    Article  CAS  PubMed  Google Scholar 

  34. Roy B, Dunbar M, Shelton RC, Dwivedi Y. Identification of microRNA-124-3p as a putative epigenetic signature of major depressive disorder. Neuropsychopharmacology : Official Publication of the American College of Neuropsychopharmacology. 2017;42(4):864–75.

    Article  CAS  Google Scholar 

  35. Wortzel JR, Turner BE, Weeks BT, Fragassi C, Ramos V, Truong T, et al. Trends in mental health clinical research: Characterizing the registry from 2007–2018. PLoS ONE. 2020;15(6):e0233996.

  36. Kaplan MJ, Klinetob NA. Childhood emotional trauma and chronic posttraumatic stress disorder in adult outpatients with treatment-resistant depression. J Nerv Ment Dis. 2000;188(9):596–601.

    Article  CAS  PubMed  Google Scholar 

  37. Durairaj RV, Koilmani ER. Environmental enrichment modulates glucocorticoid receptor expression and reduces anxiety in Indian field male mouse Mus booduga through up-regulation of microRNA-124a. Gen Comp Endocrinol. 2014;199:26–32.

    Article  CAS  PubMed  Google Scholar 

  38. Sztainberg Y, Kuperman Y, Tsoory M, Lebow M, Chen A. The anxiolytic effect of environmental enrichment is mediated via amygdalar CRF receptor type 1. Mol Psychiatry. 2010;15(9):905–17.

    Article  CAS  PubMed  Google Scholar 

  39. Brenes Saenz JC, Villagra OR, Fornaguera TJ. Factor analysis of forced swimming test, sucrose preference test and open field test on enriched, social and isolated reared rats. Behav Brain Res. 2006;169(1):57–65.

    Article  PubMed  Google Scholar 

  40. Morley-Fletcher S, Rea M, Maccari S, Laviola G. Environmental enrichment during adolescence reverses the effects of prenatal stress on play behaviour and HPA axis reactivity in rats. Eur J Neurosci. 2003;18(12):3367–74.

    Article  PubMed  Google Scholar 

  41. Weinstock M. Does Prenata stress impair coping and regulation of hypothalamic-pituitary-adrenal Axis? Neurosci Biobehav Rev. 1997;21(1):1–10.

    Article  CAS  PubMed  Google Scholar 

  42. Smalheiser NR, Lugli G, Rizavi HS, Zhang H, Torvik VI, Pandey GN, et al. MicroRNA expression in rat brain exposed to repeated inescapable shock: differential alterations in learned helplessness vs. non-learned helplessness. Int J Neuropsychopharmacol. 2011;14(10):1315–25.

    Article  CAS  PubMed  Google Scholar 

  43. Walf AA, Frye CA. The use of the elevated plus maze as an assay of anxiety-related behavior in rodents. Nat Protoc. 2007;2(2):322–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Barha CK, Pawluski JL, Galea LA. Maternal care affects male and female offspring working memory and stress reactivity. Physiol Behav. 2007;92(5):939–50.

    Article  CAS  PubMed  Google Scholar 

  45. Murthy S, Gould E. Early life stress in rodents: animal models of illness or resilience? Front Behav Neurosci. 2018;12:157.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Wang D, Levine JLS, Avila-Quintero V, Bloch M, Kaffman A. Systematic review and meta-analysis: effects of maternal separation on anxiety-like behavior in rodents. Transl Psychiatry. 2020;10(1):174.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Kokras N, Dalla C. Sex differences in animal models of psychiatric disorders. Br J Pharmacol. 2014;171(20):4595–619.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Kokras N, Antoniou K, Mikail HG, Kafetzopoulos V, Papadopoulou-Daifoti Z, Dalla C. Forced swim test: what about females? Neuropharmacology. 2015;99:408–21.

    Article  CAS  PubMed  Google Scholar 

  49. Xie L, Chen J, Ding YM, Gui XW, Wu LX, Tian S, et al. MicroRNA-26a-2 maintains stress resiliency and antidepressant efficacy by targeting the serotonergic autoreceptor HTR1A. Biochem Biophys Res Commun. 2019;511(2):440–6.

    Article  CAS  PubMed  Google Scholar 

  50. Rinaldi A, Vincenti S, De Vito F, Bozzoni I, Oliverio A, Presutti C, et al. Stress induces region specific alterations in microRNAs expression in mice. Behav Brain Res. 2010;208(1):265–9.

    Article  CAS  PubMed  Google Scholar 

  51. Ambeskovic M, Babenko O, Ilnytskyy Y, Kovalchuk I, Kolb B, Metz GAS. Ancestral stress alters lifetime mental health trajectories and cortical Neuromorphology via epigenetic regulation. Sci Rep. 2019;9(1):6389.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Bale TL, Epperson CN. Sex differences and stress across the lifespan. Nat Neurosci. 2015;18(10):1413–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Cattaneo A, Cattane N, Malpighi C, Czamara D, Suarez A, Mariani N, et al. FoxO1, A2M, and TGF-beta1: three novel genes predicting depression in gene X environment interactions are identified using cross-species and cross-tissues transcriptomic and miRNomic analyses. Mol Psychiatry. 2018;23(11):2192–208.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Bak RO, Mikkelsen JG. miRNA sponges: soaking up miRNAs for regulation of gene expression. Wiley Interdiscip Rev RNA. 2014;5(3):317–33.

    Article  CAS  PubMed  Google Scholar 

  55. Lee Y, Jeon K, Lee JT, Kim S, Kim VN. MicroRNA maturation: stepwise processing and subcellular localization. EMBO J. 2002;21(17):4663–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Falix FA, Aronson DC, Lamers WH, Gaemers IC. Possible roles of DLK1 in the notch pathway during development and disease. Biochim Biophys Acta. 2012;1822(6):988–95.

    Article  CAS  PubMed  Google Scholar 

  57. Jaric I, Rocks D, Cham H, Herchek A, Kundakovic M. Sex and estrous cycle effects on anxiety- and depression-related phenotypes in a two-hit developmental stress model. Front Mol Neurosci. 2019;12:74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Zhang W, Cheng P, Hu W, Yin W, Guo F, Chen A, et al. Inhibition of microRNA-384-5p alleviates osteoarthritis through its effects on inhibiting apoptosis of cartilage cells via the NF-kappaB signaling pathway by targeting SOX9. Cancer Gene Ther. 2018;25(11–12):326–38.

    Article  CAS  PubMed  Google Scholar 

  59. Mork L, Capel B. Oestrogen shuts the door on SOX9. BMC Biol. 2010;8(1):110.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Cattane N, Mora C, Lopizzo N, Borsini A, Maj C, Pedrini L, et al. Identification of a miRNAs signature associated with exposure to stress early in life and enhanced vulnerability for schizophrenia: new insights for the key role of miR-125b-1-3p in neurodevelopmental processes. Schizophr Res. 2018;205:63–75.

    Article  PubMed  Google Scholar 

  61. Fang Y, Qiu Q, Zhang S, Sun L, Li G, Xiao S, et al. Changes in miRNA-132 and miR-124 levels in non-treated and citalopram-treated patients with depression. J Affect Disord. 2018;227:745–51.

    Article  CAS  PubMed  Google Scholar 

  62. Gururajan A, Naughton ME, Scott KA, O'Connor RM, Moloney G, Clarke G, et al. MicroRNAs as biomarkers for major depression: a role for let-7b and let-7c. Transl Psychiatry. 2016;6(8):e862.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Sequeira A, Klempan T, Canetti L, ffrench-Mullen J, Benkelfat C, Rouleau GA, et al. Patterns of gene expression in the limbic system of suicides with and without major depression. Mol Psychiatry. 2007;12(7):640–55.

    Article  CAS  PubMed  Google Scholar 

  64. Sani G, Napoletano F, Forte AM, Kotzalidis GD, Panaccione I, Porfiri GM, et al. The wnt pathway in mood disorders. Curr Neuropharmacol. 2012;10(3):239–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Roy B, Dunbar M, Agrawal J, Allen L, Dwivedi Y. Amygdala-based altered miRNome and epigenetic contribution of miR-128-3p in conferring susceptibility to depression-like behavior via Wnt signaling. Int J Neuropsychopharmacol. 2020;23(3):165–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Yang Y, Cheng Z, Tang H, Jiao H, Sun X, Cui Q, et al. Neonatal maternal separation impairs prefrontal cortical myelination and cognitive functions in rats through activation of Wnt signaling. Cereb Cortex. 2017;27(5):2871–84.

    Article  PubMed  Google Scholar 

  67. Klinger S, Turgeon B, Levesque K, Wood GA, Aagaard-Tillery KM, Meloche S. Loss of Erk3 function in mice leads to intrauterine growth restriction, pulmonary immaturity, and neonatal lethality. Proc Natl Acad Sci U S A. 2009;106(39):16710–5.

    Article  PubMed  PubMed Central  Google Scholar 

  68. Reunanen N, Li SP, Ahonen M, Foschi M, Han J, Kahari VM. Activation of p38 alpha MAPK enhances collagenase-1 (matrix metalloproteinase (MMP)-1) and stromelysin-1 (MMP-3) expression by mRNA stabilization. J Biol Chem. 2002;277(35):32360–8.

    Article  CAS  PubMed  Google Scholar 

  69. Wojtowicz T, Mozrzymas JW. Late phase of long-term potentiation in the mossy fiber-CA3 hippocampal pathway is critically dependent on metalloproteinases activity. Hippocampus. 2010;20(8):917–21.

    Article  PubMed  Google Scholar 

  70. Kass SU, Pruss D, Wolffe A. How does DNA methylation repress transcription? Trends Genet. 1997;13(11):444–9.

    Article  CAS  PubMed  Google Scholar 

  71. Pusic AD, Kraig RP. Youth and environmental enrichment generate serum exosomes containing miR-219 that promote CNS myelination. Glia. 2014;62(2):284–99.

    Article  PubMed  Google Scholar 

  72. Timberlake Ii M, Roy B, Dwivedi Y. A novel animal model for studying depression featuring the induction of the unfolded protein response in hippocampus. Mol Neurobiol. 2019;56(12):8524–36.

    Article  CAS  PubMed  Google Scholar 

  73. Cohen H, Matar MA, Joseph Z. Animal models of post-traumatic stress disorder. Curr Protocols Neurosci. 2013;Chapter 9:Unit 9 45.

  74. Kokras N, Baltas D, Theocharis F, Dalla C. Kinoscope: an open-source computer program for behavioral pharmacologists. Front Behav Neurosci. 2017;11:88.

    Article  PubMed  PubMed Central  Google Scholar 

  75. McLean AC, Valenzuela N, Fai S, Bennett SA. Performing vaginal lavage, crystal violet staining, and vaginal cytological evaluation for mouse estrous cycle staging identification. J Vis Exp. 2012;67(67):e4389.

    Article  CAS  Google Scholar 

  76. Marcondes FK, Bianchi FJ, Tanno AP. Determination of the estrous cycle phases of rats: some helpful considerations. Braz J Biol. 2002;62(4A):609–14.

    Article  CAS  PubMed  Google Scholar 

  77. Ge SX, Jung D, Yao R. ShinyGO: a graphical gene-set enrichment tool for animals and plants. Bioinformatics. 2020;36(8):2628–9.

    Article  CAS  PubMed  Google Scholar 

  78. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10(1):1523.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  79. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25(4):402–8.

    Article  CAS  PubMed  Google Scholar 

  80. Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. Primer3--new capabilities and interfaces. Nucleic Acids Res. 2012;40(15):e115.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


We would like to thank Dr. Bhaskar Roy for his help in molecular studies. Kevin Prall, Emma Jones, and Grant Schell helped with behavior data and animal tissue collection. We would also like to acknowledge Subhan Ahmed and Genelle Samson for their blind ratings of behavior data. miRNA sequencing for this study was conducted by UAB’s Heflin Genomics Core (Birmingham, AL, USA).


The research was partly supported by grants from National Institute of Mental Health (MH082802; MH101890; MH100616; MH107183; MH112014; MH118884) and American Foundation for Suicide Prevention (DIG-0-041-18) to Dr. Dwivedi. The funding agency had no input on the study design, collection, analysis, or interpretation of the data.

Author information

Authors and Affiliations



LAM contributed to the study design, data collection, formal analysis, data curation and wrote the original draft. YD conceptualized the project and methodology, provided resources, and reviewed and edited the final draft.

Corresponding author

Correspondence to Yogesh Dwivedi.

Ethics declarations

Ethics approval and consent to participate

Animal experiments as a part of this study were approved by the Institutional Animal Care and Use Committee at the University of Alabama at Birmingham (#20086). All experiments were carried out in accordance with relevant regulations and guidelines. The ARRIVE guidelines were followed throughout the study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

McKibben, L.A., Dwivedi, Y. Early life and adult stress promote sex dependent changes in hypothalamic miRNAs and environmental enrichment prevents stress-induced miRNA and gene expression changes in rats. BMC Genomics 22, 701 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: