- Open Access
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 volume 22, Article number: 701 (2021)
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.
Early life stress (ELS), namely abuse, neglect, and household dysfunction, is estimated to affect 64% of the US population . 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 , its symptoms are particularly debilitating. MDD onset has also been strongly associated with a precipitating stressful event, especially within the previous month . Moreover, ELS and stress during adulthood can interact producing more severe or long-lasting symptoms than individually , also described as the two/three-hit hypothesis of stress susceptibility . There are also significant sex differences in the response to ELS . In the US, men and women experience similar levels of ELS , yet women are twice as likely to develop MDD . There are only a limited number of studies that have examined sex-mediated depressive response to ELS .
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  and frontal cortex . Moreover, it plays a particularly important role in the symptomatology of MDD . The hypothalamus receives sensory cues predicting reward (or lack thereof) and stimulates dopamine downstream . In rodents, disruption of dopaminergic projections to the hypothalamus leads to reduced sucrose preference (i.e. more anhedonia-like behavior) . 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  and directly injecting calcium channel inhibitors into the hypothalamus increases swimming behavior in the forced swim test (FST) indicated reduced behavioral despair . 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 . 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 ), 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 . MS has been well established as a rodent model of ELS associated with depression- and anxiety-like behavior . 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 . 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 . Increased miR-124 in the dentate gyrus after 90 min MS has also been reported . While miRNAs serve a primary function to regulate gene expression, they are also regulated by epigenetic modifications such as methylation . 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 . Using chronic CORT administration to induce depression-like behavior in rats, our group found promoter hypomethylation of miR-124-3p . 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 . Depressed individuals with a history of ELS may respond less to antidepressant drug therapy than patients with no ELS history ; 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  and lower blood CORT levels after RS . 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.
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).
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).
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.
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.
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 . 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 ;). 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 . 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 . 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 , especially in the FST .
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 ; 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  as well as in male rats prefrontal cortex who experienced repeated and ancestral stress . Sex differences following ELS, especially stress sensitivity in males prior to puberty, may result from interactions between stress neurobiology and peripubertal sex hormone changes . 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  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 . 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 . 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 . 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 . 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  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 . SOX 9 is also an embryonic transcription factor which responds to estrogen during sex determination . 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  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  found an increase in miR-29a following MS. On the other hand, Cattane, Mora  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 . 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 . 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 . 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 . 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 . Our group recently found decreased WNT5B expression associated with increased miR-128 in rodents with learned helplessness . Another study found that neonatal MS reduced myelination in the prefrontal cortex via Wnt signaling changes .. MAPK6 encodes ERK3 and is widely expressed in the brain. ERK3 plays an important role in neonatal growth . MAPK pathways have also been found to regulate some MMPs . Wojtowicz and Mozrzymas  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 , 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 . 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.
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.
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.
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.
Sucrose preference test
Sucrose preference was tested immediately following RS to assess anhedonia as previously described . 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 .
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 .
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 .
Shuttle escape test
As reported previously , 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 . 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).
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.
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.
miRbase.org 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 , 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  was used to explore GO associated with miRNA target genes. Metascape  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  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 ). 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 http://www.ncbi.nlm.nih.gov/bioproject/728805 and accession #: PRJNA728805.
early life stress
major depressive disorder
- HPA :
hypothalamic pituitary adrenal
corticotropin releasing hormone
adrenocorticotropin releasing hormone
- miRNA, miR:
RNA-induced silencing complex
- EPM :
elevated plus maze
forced swim test
counts per million
methylated DNA immunoprecipitation
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.
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. https://doi.org/10.4088/JCP.14m09298.
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. https://doi.org/10.1176/ajp.156.6.837.
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. https://doi.org/10.3389/fpsyt.2020.00744.
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. https://doi.org/10.1016/j.psyneuen.2013.06.008.
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. https://doi.org/10.1530/JOE-15-0121.
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 https://www.acf.hhs.gov/cb/research-data-technology/statistics-research/child-maltreatment.
CfBHSa Q. 2016 National Survey on drug use and health: methodological summary and definitions. Substance Abuse and Mental Health Services Administration: Rockville, MD; 2017.
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. https://doi.org/10.1038/s41398-019-0423-8.
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. https://doi.org/10.1001/jama.284.5.592.
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. https://doi.org/10.1016/S0006-8993(99)02307-0.
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. https://doi.org/10.1111/j.1365-2826.2012.02377.x.
Price JL. Comparative aspects of amygdala connectivity. Ann N Y Acad Sci. 2003;985(1):50–8. https://doi.org/10.1111/j.1749-6632.2003.tb07070.x.
Öngür D, Price JL. The Organization of Networks within the orbital and medial prefrontal cortex of rats. Cereb Cortex: Monkeys and Humans; 2000.
Price JL, Drevets WC. Neurocircuitry of mood disorders. Neuropsychopharmacology : Official Publication of the American College of Neuropsychopharmacology. 2010;35(1):192–216. https://doi.org/10.1038/npp.2009.104.
Sakurai T. The role of orexin in motivated behaviours. Nat Rev Neurosci. 2014;15(11):719–31. https://doi.org/10.1038/nrn3837.
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. https://doi.org/10.1111/j.1460-9568.2006.04944.x.
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. https://doi.org/10.1016/0006-8993(93)90742-6.
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. https://doi.org/10.1016/0091-3057(90)90079-W.
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. https://doi.org/10.1016/S0006-3223(01)01157-X.
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. https://doi.org/10.1016/j.psyneuen.2008.03.008.
Dwivedi Y. microRNA-124: a putative therapeutic target and biomarker for major depression. Expert Opin Ther Targets. 2017;21(7):653–6. https://doi.org/10.1080/14728222.2017.1328501.
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. https://doi.org/10.1038/mp.2011.162.
Dwivedi Y. Emerging role of microRNAs in major depressive disorder: diagnosis and therapeutic implications. Dialogues Clin Neurosci. 2014;16(1):43–61. https://doi.org/10.31887/DCNS.2014.16.1/ydwivedi.
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. https://doi.org/10.1007/s40291-020-00464-9.
Bartel DP. Metazoan MicroRNAs. Cell. 2018;173(1):20–51. https://doi.org/10.1016/j.cell.2018.03.006.
Allen L, Dwivedi Y. MicroRNA mediators of early life stress vulnerability to depression and suicidal behavior. Mol Psychiatry. 2020;25(2):308–20. https://doi.org/10.1038/s41380-019-0597-8.
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. https://doi.org/10.1080/15592294.2018.1451720.
Vetulani J. Early maternal separation: a rodent model of depression and a prevailing human condition. Pharmacol Rep. 2013;65(6):1451–61. https://doi.org/10.1016/S1734-1140(13)71505-6.
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. https://doi.org/10.1523/JNEUROSCI.1436-10.2010.
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. https://doi.org/10.1016/j.bbr.2016.05.033.
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. https://doi.org/10.1016/j.bbagrm.2016.10.006.
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. https://doi.org/10.1016/j.jad.2017.05.017.
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. https://doi.org/10.1038/npp.2016.175.
Wortzel JR, Turner BE, Weeks BT, Fragassi C, Ramos V, Truong T, et al. Trends in mental health clinical research: Characterizing the ClinicalTrials.gov registry from 2007–2018. PLoS ONE. 2020;15(6):e0233996.
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. https://doi.org/10.1097/00005053-200009000-00006.
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. https://doi.org/10.1016/j.ygcen.2014.01.005.
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. https://doi.org/10.1038/mp.2009.151.
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. https://doi.org/10.1016/j.bbr.2005.12.001.
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. https://doi.org/10.1111/j.1460-9568.2003.03070.x.
Weinstock M. Does Prenata stress impair coping and regulation of hypothalamic-pituitary-adrenal Axis? Neurosci Biobehav Rev. 1997;21(1):1–10. https://doi.org/10.1016/S0149-7634(96)00014-0.
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. https://doi.org/10.1017/S1461145710001628.
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. https://doi.org/10.1038/nprot.2007.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. https://doi.org/10.1016/j.physbeh.2007.06.022.
Murthy S, Gould E. Early life stress in rodents: animal models of illness or resilience? Front Behav Neurosci. 2018;12:157. https://doi.org/10.3389/fnbeh.2018.00157.
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. https://doi.org/10.1038/s41398-020-0856-0.
Kokras N, Dalla C. Sex differences in animal models of psychiatric disorders. Br J Pharmacol. 2014;171(20):4595–619. https://doi.org/10.1111/bph.12710.
Kokras N, Antoniou K, Mikail HG, Kafetzopoulos V, Papadopoulou-Daifoti Z, Dalla C. Forced swim test: what about females? Neuropharmacology. 2015;99:408–21. https://doi.org/10.1016/j.neuropharm.2015.03.016.
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. https://doi.org/10.1016/j.bbrc.2019.02.078.
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. https://doi.org/10.1016/j.bbr.2009.11.012.
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. https://doi.org/10.1038/s41598-019-42691-z.
Bale TL, Epperson CN. Sex differences and stress across the lifespan. Nat Neurosci. 2015;18(10):1413–20. https://doi.org/10.1038/nn.4112.
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. https://doi.org/10.1038/s41380-017-0002-4.
Bak RO, Mikkelsen JG. miRNA sponges: soaking up miRNAs for regulation of gene expression. Wiley Interdiscip Rev RNA. 2014;5(3):317–33. https://doi.org/10.1002/wrna.1213.
Lee Y, Jeon K, Lee JT, Kim S, Kim VN. MicroRNA maturation: stepwise processing and subcellular localization. EMBO J. 2002;21(17):4663–70. https://doi.org/10.1093/emboj/cdf476.
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. https://doi.org/10.1016/j.bbadis.2012.02.003.
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. https://doi.org/10.3389/fnmol.2019.00074.
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. https://doi.org/10.1038/s41417-018-0029-y.
Mork L, Capel B. Oestrogen shuts the door on SOX9. BMC Biol. 2010;8(1):110. https://doi.org/10.1186/1741-7007-8-110.
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. https://doi.org/10.1016/j.schres.2018.07.030.
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. https://doi.org/10.1016/j.jad.2017.11.090.
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. https://doi.org/10.1038/tp.2016.131.
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. https://doi.org/10.1038/sj.mp.4001969.
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. https://doi.org/10.2174/157015912803217279.
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. https://doi.org/10.1093/ijnp/pyz071.
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. https://doi.org/10.1093/cercor/bhw121.
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. https://doi.org/10.1073/pnas.0900919106.
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. https://doi.org/10.1074/jbc.M204296200.
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. https://doi.org/10.1002/hipo.20787.
Kass SU, Pruss D, Wolffe A. How does DNA methylation repress transcription? Trends Genet. 1997;13(11):444–9. https://doi.org/10.1016/S0168-9525(97)01268-7.
Pusic AD, Kraig RP. Youth and environmental enrichment generate serum exosomes containing miR-219 that promote CNS myelination. Glia. 2014;62(2):284–99. https://doi.org/10.1002/glia.22606.
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. https://doi.org/10.1007/s12035-019-01687-6.
Cohen H, Matar MA, Joseph Z. Animal models of post-traumatic stress disorder. Curr Protocols Neurosci. 2013;Chapter 9:Unit 9 45.
Kokras N, Baltas D, Theocharis F, Dalla C. Kinoscope: an open-source computer program for behavioral pharmacologists. Front Behav Neurosci. 2017;11:88. https://doi.org/10.3389/fnbeh.2017.00088.
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. https://doi.org/10.3791/4389.
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. https://doi.org/10.1590/S1519-69842002000400008.
Ge SX, Jung D, Yao R. ShinyGO: a graphical gene-set enrichment tool for animals and plants. Bioinformatics. 2020;36(8):2628–9. https://doi.org/10.1093/bioinformatics/btz931.
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. https://doi.org/10.1038/s41467-019-09234-6.
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. https://doi.org/10.1006/meth.2001.1262.
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. https://doi.org/10.1093/nar/gks596.
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.
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
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
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). https://doi.org/10.1186/s12864-021-08003-4
- Early life stress
- Maternal separation
- Restraint stress
- microRNA sequencing