Behavioural characterisation of chronic unpredictable stress based on ethologically relevant paradigms in rats

The chronic unpredictable stress (CUS) paradigm is extensively used in preclinical research. However, CUS exhibits translational inconsistencies, some of them resulting from the use of adult rodents, despite the evidence that vulnerability for many psychiatric disorders accumulates during early life. Here, we assessed the validity of the CUS model by including ethologically-relevant paradigms in juvenile rats. Thus, socially-isolated (SI) rats were submitted to CUS and compared with SI (experiment 1) and group-housed controls (experiment 1 and 2). We found that lower body-weight gain and hyperlocomotion, instead of sucrose consumption and preference, were the best parameters to monitor the progression of CUS, which also affected gene expression and neurotransmitter contents associated with that CUS-related phenotype. The behavioural characterisation after CUS placed locomotion and exploratory activity as the best stress predictors. By employing the exploratory factor analysis, we reduced each behavioural paradigm to few latent variables which clustered into two general domains that strongly predicted the CUS condition: (1) hyper-responsivity to novelty and mild threats, and (2) anxiety/depressive-like response. Altogether, the analyses of observable and latent variables indicate that early-life stress impairs the arousal-inhibition system leading to augmented and persistent responses towards novel, rewarding, and mildly-threatening stimuli, accompanied by lower body-weight gain.

It has been recently discussed that the CUS effects and the reliability of the model depend on the severity of the stress protocol 20 , which suggests that the addition of multiple and more severe stressors would proportionally strengthen the CUS-related phenotype. With this in mind, and considering that we failed to reproduce some of the typical CUS effects, we made our protocol more intense by including two additional stressors (Fig. 1D,H). Also, the stress exposures were now more frequent (i.e., from 15 to 24 days) and longer (i.e., from 20 to 22 hours) than before. This new protocol was assessed in detail in experiment 2. There, the progression of CUS effects was monitored after 10, 20, and 30 days of CUS, which corresponded to PNDs 42/43, 52/53, and 62/63, respectively (Fig. 1A).
We measured body weight, sucrose consumption and preference, and included locomotor activity in the OFT to have an easy but informative behavioural parameter known to be responsive to stress 14,17 . Testing took place 10-days apart, to minimise habituation and cross-over effects between and within tests, and to enable the comparison of stress chronicity between studies ranging from days to weeks 20 . Data correspond to 40 rats, that is, 20 unstressed (CON) and 20 stressed (CUS) rats that started the experiment (Fig. 1A). Before CUS, no significant differences in body weight, locomotion, water and sucrose consumption, and sucrose preference were observed (Fig. 3). As a consequence of CUS, the stressed animals had significantly less weight than the ones in the control group (P = 0.0001, η 2 = 0.45). Although both groups showed a significant increase in body weight throughout the study (P = 0.0001, η 2 = 0.99), weight gain was less pronounced in CUS rats (P = 0.0001, η 2 = 0.25) (Fig. 3A). The locomotor activity of all animals decreased from day 10 to day 20 and remained unchanged thereafter (P = 0.004, η 2 = 0.14). CUS rats, however, displayed higher levels of locomotion throughout the three-time points measured (P = 0.0001, η 2 = 0.35) (Fig. 3B). Sucrose consumption (P = 0.0001, η 2 = 0.59) (Fig. 3C) and preference (P = 0.0001, η 2 = 0.28) (Fig. 3D)  28 rats were weighted and screened behaviourally using the Sucrose Preference Test (SPT) and the Open Field Test (OFT). Based on the total distance travelled during the OFT, the sucrose preference, and body weight, the animals were balanced and semi-randomly allocated to either the control (CON) or chronic unpredictable stress (CUS) group (n = 20 each). No animals were excluded from the study. The CUS protocol started on PND 33 and ended on PND 63. Body weight, sucrose consumption and preference, and locomotor activity were measured at PND 42/43, PND 52/53, and PND 62/63 to monitor the progression of CUS. Afterwards, animals from each group were balanced based on the level of those variables at PNDs 62 and 63 and then allocated semi-randomly to subgroups for either neurochemical or behavioural analyses (20 animals each, 10 per group). Animals allocated into the gene expression/neurochemical analysis were euthanised at PND 65, whereas the rest of the animals were behaviourally tested from PND 65 to PND 69 (see methods for details). The CUS protocol consisted of the following stressors: (B) Wet bedding: 300 mL of water was poured on and mixed with 1 L of sawdust bedding. (C) Sleep deprivation: a cylindrical, wooden pedestal (6 cm diameter and 5 cm height) was placed on the floor of the cage opposite to the food/water compartment. The cage was flooded with tap water 3 cm deep, allowing the animal to stand on the bottom of the pedestal but denying the possibility of sleeping. (D) Electric Footshocks: rats were placed in an Ugo Basile Automatic Reflex Conditioner (Global Biotech, USA; 40 cm length × 20 cm width × 22 cm height) for the administration of 0.8 mA shocks through the grid floor of the chamber. Each 10-minutes session consisted of two trials of 5 minutes with two stages each: the delivery stage with six presentations of 5 shocks interspaced by 30-seconds intervals and the resting stage of 2 minutes in the chamber without receiving shocks. (F) Cage tilting: the cage was tilted up to 45 degrees with food and water located at the higher top. (G) Water deprivation. (F) Food deprivation. (H) Confinement: rats were individually housed in small cages (20 cm length × 10 cm width × 13 cm height). Rats continued on social isolation in standard cages during the remaining time between stressors and from PND 63 onwards, in which the behavioural testing took place. The CON group animals were group-housed (4-6 animals per cage) throughout the entire experiment, except during the SPT. Stressors followed a semi-random order to reduce their predictability. All stressors lasted 22 hours, except for the foot-shocks sessions that lasted 10 minutes each. Animals were exposed to every stressor 3 to 4 times throughout the protocol.
increased with time in all rats with no differences observed between groups. When correcting sucrose consumption by body weight (Fig. 3E), the opposite pattern appeared, that is, intake slightly decreased over days in both groups (P = 0.003, η 2 = 0.14). This pattern was more noticeable in the CUS group and also when comparing the second with the third and fourth SPT measurements, which did not differ to each other. The reason behind the gradual reduction in sucrose consumption was that body weight increased at a higher rate than sucrose and water intake, which made the ratio (ml/g) to remain almost the same over the 10-days blocks. When compared by groups, the corrected consumption of sucrose was significantly higher in CUS animals throughout the protocol (P = 0.002, η 2 = 0.22) (Fig. 3E). To identify which combination of variables has the largest explicative power for predicting the CUS condition, a linear multiple regression analysis was conducted with the dependent variable groups categorised with dummy codes of "0" for controls and "1" for CUS -the condition to be predicted in the models. According to the stepwise method of analysis, body weight was the best predictor of CUS (adjusted R 2 = 0.44, P = 0.001), followed by the combination of body weight and locomotion (adjusted R 2 = 0.56, P = 0.001). Sucrose consumption and preference variables yielded no significant contributions to the whole prediction and were excluded from the analysis. When including the variable by blocks, adding sucrose consumption and preference in a second block and in different combinations over body weight, made no significant improvements to the prediction of CUS. Although the body weight alone showed the highest adjusted regression coefficient, a model including locomotion activity had the lowest Akaike information criterion (AIC) values with a ΔAIC of 9 below a model with only body weight (for details see statistical analysis section).
Discussion. In experiment 1, rats were exposed to a milder CUS protocol in order to determine the contribution of SI to the CUS effects on body weight and sucrose consumption and preference. According to the eta square coefficients, the unique parameter strong enough to discriminate the effects of CUS from the other conditions was Figure 2. Effects of chronic unpredictable stress (CUS) and social isolation (SI) before and after a 30-days protocol (experiment 1). Body weight (A), water (B) and sucrose consumption (C), sucrose preference (D), and sucrose consumption corrected by body weight (E). Sucrose and water consumption was expressed as percentages [(intake ml/200 ml)x100]. Similarly, preference was calculated as percentages [(sucrose consumption/(sucrose + water consumption)) X 100]. Sucrose corrected by body weight was computed as follows: consumption ml/body weight in grams. Narrower bars within the bold-line bars corresponded to baseline levels. Control group (CON). Between-groups planned comparisons: CUS differed from the other groups, ++ P < 0.01. SI differed from the other groups, # P < 0.05. CON differed from the other groups,*P < 0.05. body weight, with a coefficient of 36%. Surprisingly, our 30-days CUS protocol was insufficient to reduce sucrose consumption and preference, contrary to what is expected with this test. When correcting the sucrose consumption by body weight, a weak effect of CUS was revealed (10%). Contrary to the main trend in the field, we found an increase, instead of a decrease, in sucrose consumption in CUS rats. As SI and CUS animals showed almost the same intake after the correction, the other CUS stressors seem to contribute little or nothing to the whole effect, suggesting that SI was responsible for increasing sucrose intake and preference relative to group-housed controls. In addition, CUS significantly decreased sucrose preference relative to SI, an effect that coincides with that reported in early CUS studies in which control animals were housed in SI 14,15,17 . In experiment 2, a more intense protocol was implemented to strengthen the weak CUS-related profile obtained in experiment 1 and to better characterise our CUS protocol on a wider range of neurobehavioral readouts. According to the eta square coefficients, the reduction in body weight was, again, the most important alteration induced by stress (45%), followed by the increase in locomotor activity (35%), and sucrose consumption (18%). When comparing all variables with each other in order to assess their individual and additive prediction capacities, the regression analyses showed that SPT variables are poor predictors which did not sum up to the effects already produced by body weight and locomotion. The lower body weight of stressed animals (Figs 2A and 3A) replicated previous findings and could have resulted from a combination of factors such as an increased caloric requirement to maintain the body temperature after sleep deprivation and wet bedding 33 , with SI impeding communal thermoregulation via curling. Hypophagia may be an alternative, non-exclusive explanation 34 ; however, in stressed animals, sucrose consumption was consistently higher throughout the experiment 2. As SI alone failed to reduce body weight (experiment 1, Fig. 2A), we attributed this effect to the CUS stressors other than SI.
The second most important parameter affected by stress was locomotion, which was consistently higher in CUS animals (Fig. 3B), supporting evidence obtained from SI, CUS, and other stress models [35][36][37] . It is unlikely to consider such an effect on locomotion as a by-product of SI. Many studies conducted in different rat strains showed that hyperlocomotion is an expected outcome of CUS, regardless of the type of housing and the strain used [36][37][38][39] . It is worth noting that many studies have reported a reduction of locomotion in the OFT. However, most of these studies have in common the use of larger OFT arenas, shorter periods of observation (e.g., 3-5 minutes), and more severe CUS protocols [40][41][42][43][44][45][46][47] . Thus, as stressed rats may be more sensitive to the anxiogenic properties associated with larger arenas, they would freeze soon after being placed in the OF, leading to rather Figure 3. Assessment of the CUS progression within the 30-days protocol in Experiment 2. Body weight (A), locomotor activity (B), sucrose and water consumption (C), sucrose preference (D), and sucrose consumption corrected by body weight (E). Sucrose and water consumption was expressed as percentages [(intake ml/200 ml)x100]. Similarly, preference was calculated as percentages [(sucrose consumption/(sucrose + water consumption)) X 100]. Sucrose corrected by body weight was computed as follows: consumption ml/ body weight in grams. Dashed lines corresponded to water consumption. Control group (CON). Chronic unpredictable stress group (CUS). Between-groups comparison: ***P < 0.0001, *P < 0.05. low levels of activity. Such an effect is even more noticeable if only the first minutes of the test are analysed. When freezing behaviour was measured in a 5-minutes OFT, it was observed that CUS rats spent about 66% of the time immobile. If the anxiolytic drug diazepam is previously administered, such an effect could be prevented 48 . In a very similar CUS experiment, the use of an arena 33% smaller (40 cm 2 ), produced a shift in the activity levels, with stressed animals now exhibiting significantly higher locomotion than controls 49 . Finally, intense CUS protocols coincide in producing very low levels of exploratory activity regardless of the arena dimensions and experimental settings 40,42,47,50 . Overall, this evidence supports our opinion that CUS-induced hypolocomotion is quite uncommon and may result from the combination of different methodological factors.
The SPT was the weakest test in assessing the effects of CUS. In both experiments, sucrose consumption was consistently higher rather than lower after stress (Figs 2C, 3D), but only when correcting by body weight (Figs 2E, 3E). The latter highlights the importance of performing the correction to unmask the effects of CUS, especially when recognising the well-known relationship between body weight and fluid intake and that CUS affects body weight substantially 51 . For instance, significant reductions in sucrose consumption and preference can disappear after the correction 52,53 or, in contrast, significant differences emerge when comparing the corrected data, as it was in our case. Our results suggest that SI alone or within the CUS protocol is capable of increasing sucrose consumption and preference in agreement with SI studies 25,26,28,36,39 , but in disagreement with the CUS evidence [14][15][16]39,46,[54][55][56][57][58] . As we did not include a condition of CUS rats housed in groups, we could not clear up how much of our particular CUS effect was due to SI. We decided to maintain our control animals grouped, as the practice of maintaining controls in SI has been discontinued because of misunderstandings as regards the interpretations given to the effects, or the lack thereof when comparing SI versus CUS. Although having our control rats in groups affects comparability with a bulk of literature in the field, it was still necessary in order to better characterise our CUS protocol, especially in experiment 2. The downside was, however, the impossibility to separate out the effects of SI per se from the full CUS procedure.
We acknowledged that many CUS studies including SI as one stressor found a reduction in sucrose consumption or preference when compared with SI alone 15,45,49,54,58,59 . However, other CUS protocols also including SI were unable to replicate the reduction in sucrose consumption or preference when compared with SI controls 52,53,60 or found differences only when measuring the intake during the night-time but not throughout the day 61 . In other experiments, marginal and significant increases in sucrose or saccharin intake and preference after CUS were observed on different rat strains when compared with SI 36,39,52,59,62 . It is worth noting that, except in a few cases, most CUS studies have been conducted in adult rats, so little is known about the effects of CUS at earlier stages. The few available studies in young animals indicate that juvenile rats tended to consume and prefer sucrose more than adults 39,62,63 . After CUS, juvenile rats showed significantly higher sucrose preference than older counterparts 63 , which somehow support our findings obtained in juvenile rats. All this evidence indicates that the SPT may lead to contradictory results depending on the variations in stress intensity and duration, the type of parameter used (e.g., intake or preference), the strain, and the age of the subjects. In our hands, however, neither the sucrose concentration nor the rat strain accounted for those discrepancies. In fact, either with CUS or SI, we have obtained almost the same results when using sucrose at 1% (current data) or 32% 24,28 or when employing Wistar or Sprague-Dawley rats 24,28 . In conclusion, to assess the progressive effects of CUS, the best variables are the reduction in body weight gain followed by hyperlocomotion, and to a lesser extent, the increase in sucrose consumption.

Behavioural characterisation of the CUS effects.
Results. This phase comprised the remaining 20 rats (10 per group) that did not undergo neurochemical analyses. To provide further information about the CUS effects on spontaneous exploratory activity, the third OFT carried out at PND 63 ( Fig. 1A) was fully examined by breaking down the timeframe of analysis (i.e., minute by minute) and by measuring new behaviours (i.e., rearing and grooming). We found that locomotor activity decreased over minutes in both groups (P = 0.0001, η 2 = 0.55), but stressed rats showed significantly higher locomotion than CON rats at almost every single minute of the test (P = 0.01, η 2 = 0.31) (Fig. 4A), with no interaction detected between Minutes and Treatment. Time spent rearing showed an irregular tendency to increase during the first minutes of the test (Fig. 4B), which was significantly higher in CUS rats (P = 0.047, η 2 = 0.10) (Fig. 4B). No main effects for Session and Treatment were observed. With regards to grooming (Fig. 4C), the time spent on this behaviour increased over minutes in both groups with a rather irregular pattern (P = 0.024, η 2 = 0.13) characterised by three bouts of pronounced activity, which were significantly lower in CUS than in CON group (P = 0.047, η 2 = 0.12) (other OFT parameters are shown in Supplementary Table 1). At PND 65, animals were exposed to the EPM for 5-minutes. In this test, CUS neither affected the traditional anxiety-like parameters (i.e., the time spent in and the number of entries to the closed and open arms, Fig. 4D,E) nor the risk-assessment behaviours -like the stretch attempt posture (SAP)-, except head-dipping (HD) time (P = 0.03, η 2 = 0.23), which was lower in CUS animals than in CON rats ( Fig. 4F) (other EPM parameters are shown in Supplementary Table 2). The next test was the ORT, which consisted of three phases. In the first one, carried out at PND 66 and 24 hours before the sample trial, all animals were habituated to the arena for 5 minutes without any objects. We found that locomotion decreased over minutes in both groups (P = 0.0001, η 2 = 0.44), with CON rats showing a more pronounced activity decline when compared to CUS animals (P = 0.022, η 2 = 0.15) (data not shown), but with a marginal effect of Treatment (Fig. 4G). Other behaviours, such as rearing and grooming, did not differ between groups in the habituation phase (for details see Supplementary Table 3). During the sample trial, rats were exposed to two identical objects for 5 minutes. The groups did not differ in the exploration time of the objects and showed no significant preferences for neither of them (Fig. 4H left). Thirty minutes later, the test trial took place just like the sample trial, except that one of the objects was changed. There, both groups explored more the novel object than the familiar one (P = 0.009, η 2 = 0.32), with CUS animals exploring significantly more than CON rats irrespective of the object (P = 0.004, η 2 = 0.38) (Fig. 4H right). When comparing total exploration time between sample and test trials, (2019) 9:17403 | https://doi.org/10.1038/s41598-019-53624-1 www.nature.com/scientificreports www.nature.com/scientificreports/ no significant differences were found. A detailed within-group analysis showed that exploration of CON rats decreased 23% from sample to test trial (P = 0.03, η 2 = 0.44), whereas in CUS rats it rather increased 10%, without being significant (Fig. 4H right). As the perception of novelty during the test trial depends on how much the objects were explored during the sample trial, the exploration time of each object was normalised accordingly. Relative to the sample trial, stressed rats spent significantly more time exploring the familiar object during the test trial (P = 0.016, η 2 = 0.28), with no group differences observed for the novel object (Fig. 4I). A discrimination index was calculated, using the normalised exploration times for each object to obtain an additional recognition memory parameter, which was unbiased by the previous differences in object exploration. We found no group differences in the discrimination index (Supplementary Table 4). Finally, the analysis of the spontaneous activity www.nature.com/scientificreports www.nature.com/scientificreports/ revealed that, in all animals, locomotion decreased between trials (P = 0.02, η 2 = 0.15), with stress rats displaying higher locomotion in all sessions as compared with CON animals (P = 0.004, η 2 = 0.21) (Fig. 4G). The last test was the FST, carried out at PND 68 and 69. In the 15-minutes pre-test session (Fig. 4J), stress increased the immobility time and decreased struggling activity as compared with the CON group (P = 0.026, η 2 = 0.25). A detailed analysis of the active behaviours showed that climbing, but not swimming, was significantly reduced by stress (P = 0.01, η 2 = 0.34) (Fig. 4J). Twenty-four hours later, the test session (5 minutes) was performed. There, no significant group differences for any FST behaviour were observed (Fig. 4K). When comparing pre-test and test, immobility and struggling activity (P = 0.028, η 2 = 0.24) increased and decreased in all animals, respectively. Out of these active behaviours, swimming, but not climbing, was significantly different between sessions (P = 0.001, η 2 = 0.48). A detailed within-group comparison showed that immobility time in CON rats increased significantly from pre-test to test (13%) until levelling off with the immobility time of CUS rats (P = 0.027, η 2 = 0.43), which showed a non-significant increase of 4% over the already high levels of immobility seen in the pre-test.
The linear regression analysis was conducted only with the variables that were significantly different between groups. As CUS had the largest effect on body weight on previous analyses, the initial model tested whether body weight -measured when the behavioural characterisation took place-was still the best predictor of CUS in the subset of rats that completed this phase. According to both the stepwise method and the individual analysis by blocks, body weight appeared as the best predictor of CUS (adjusted R 2 = 0.36, P = 0.003, AIC = 25.2). Thus, subsequent block models were built up by adding predictors to body weight to identify the minimal combination of variables with a greater predictive capacity than a model with only body weight (reference model) and a model including all variables (saturated model). After testing all possible combinations, the best fit was obtained (adjusted R 2 = 0.73, P = 0.001, AIC = 11.5) with a model including body weight plus locomotion in the OFT and locomotion and object exploration (i.e., time exploring both the familiar and novel objects and the normalised exploration of the familiar object) in the test phase of the ORT. This model was also better than the saturated model (adjusted R 2 = 0.67, P = 0.08, AIC = 12.1).
Discussion. According to the eta squared coefficients, the behavioural characterisation of CUS after completing the protocol revealed that locomotor activity during the third OFT (31%) and throughout all ORT sessions (21%) had, again, one of the largest size effects, with CUS rats displaying the highest activity levels. Rearing (10%) and grooming (12%) were weakly affected by CUS and only during the OFT, with no effects detected on the ORT. Object exploration during the test session of the ORT offered the most considerable differentiation between groups (38%), with CUS rats not only displaying higher overall exploration but also spending more time exploring the familiar object (21%). In fact, when all variables were allowed to compete with each other to determine their individual and additive predictive capacity, the hyperlocomotion and object exploration in mild-stressing contexts were the best predictors of CUS -after body weight. Hyperlocomotion and sustained exploration after repeated testing on similar environments may result from an impairment in the arousal-inhibition system leading to hypervigilance, which constitutes a key factor of anticipatory anxiety induced by unpredictable stress 64 . Such a failure to habituate may be attributable to an augmented response to mild threats (e.g., OFT and ORT arenas) and a compromised ability for recognising cues that should have already signalled safety 64,65 . Although no differences in recognition memory indexes were observed, the fact that CUS rats explored more the familiar object may constitute an index of perseveration, which refers to an impairment in memory inhibition of previously learned pieces of information compromising object discrimination 66 . Certainly, perseveration has been proved to be augmented by stress 67 . It has been repeatedly observed that CUS decreases the discrimination index either marginally 68 or significantly [69][70][71][72][73][74] , although others have reported that only a subset of CUS responders shows memory impairments in the ORT 75 . The lack of significant differences in the discrimination index may be related to the features of the objects used. As our objects were quite plain, smooth, and unclimbable, the discrimination between novel and familiar objects could have become too hard. In this sense, it has been reported that objects that can be climbed over are explored longer than objects that can only be touched, leading to greater discrimination indexes 76 . Nevertheless, in other tasks not related to recognising objects, CUS has induced impairments on different memory domains 62,65,77 . On the EPM, most parameters were unaffected by CUS, except for the time spent on HD (23%). As this behaviour is part of the risk-assessment repertoire of the rat, a reduced motivation to explore putative dangerous zones (e.g., the borders of the open arm) can be interpreted as a state of heightened, negative emotionality, as suggested elsewhere [78][79][80] . However, the lack of effects of CUS in the EPM is not surprising at all. For instance, CUS protocols in rats extended from 2 to 6 weeks increased -instead of decreasing-the time spent on the open arms 10,57,71 , whereas a shorter protocol (i.e., 1.5 weeks) had no effects on this parameter 12 . Interestingly, a lighter, but considerably longer CUS protocol of 11 weeks was not able to produce anxiogenic effects on the EPM 81 . In mice, a reduction in the time spent on the open arms was observed after six weeks 69 , but not after nine weeks of CUS 82 . In this latter study, CUS mice stayed longer in the open arms of the EPM, spent more time in the centre of an OFT, and the bright side of the dark/light test 82 , suggesting a consistent anxiolytic profile induced by CUS, contrary to what can be expected. Although all these discrepancies could be attributed to differences in CUS protocols, the EPM seems to be an unstable test that is more likely to yield false positives when treatments increase locomotion and exploratory activity, as may be the case for CUS.
Regarding the FST, we found that CUS increased immobility and diminished climbing times only in the pre-test session. Others have also found that CUS induced high immobility levels only at the first FST exposure 45,48,63,83 . These results contrast with evidence showing that CUS lasting three to six weeks is sufficient to increase immobility in test session 39,50,57,67,75 . However, in a six-weeks CUS protocol immobility increased neither in the pre-test nor in the test session in male rats 49 ; and in a mild but more extended CUS protocol (11 weeks), immobility in the test session was unaffected by stress 49,81 , suggesting that differences in CUS intensity may lead to conflicting results. Immobility in the FST can be initially interpreted as an adaptive response towards inescapability 84 ; nevertheless, longer immobility periods after stress may suppose a failure in displaying active, www.nature.com/scientificreports www.nature.com/scientificreports/ escape-oriented behaviours 85 . In this sense, our results indicate that CUS increased passive-coping response towards acute and intense stress, which can be taken as a depressogenic effect of CUS. In the FST, immobility increases in the test session as animals learn that the situation is inescapable. We found that CON rats increased this behaviour, but CUS animals almost did not. It can be argued that CUS affected learning processes or the ability to display adaptive postures in the test session, as suggested elsewhere 84 . Interestingly, this lack of change in the immobility time from the pre-test to the test session has also been reported in at least one additional CUS study with male rats 49 . Altogether, the behavioural characterisation after CUS showed that its best predictor was the impairment in the arousal-inhibition system evidenced by the excessive exploratory and locomotor activity seen throughout all ORT and OFT sessions, which positioned these tests as the most informative paradigms in this section. It must be considered, nevertheless, that the behavioural pattern described here might not generalise to other studies in which CUS rats were housed in groups or control rats were housed in SI.
Discussion. CUS slightly reduced the BDNF levels, especially in the hippocampus, suggesting that the BNDF-dependent signalling related to neuronal differentiation, survival, and structural and synaptic plasticity was down-regulated after stress, in agreement with previous reports 46,77,86 . Such a decrease has been associated with the development of anxiety-and depression-related behaviours in animals and humans 87 . Indeed, we found www.nature.com/scientificreports www.nature.com/scientificreports/ that CUS animals with lower accumbal BDNF levels showed hyperlocomotion (r = −0.601, P = 0.033) and higher sucrose consumption (r = −0.592, P = 0.036). The reduced accumbal BDNF levels of CUS rats were also related with higher sucrose consumption (r = −0.566, P = 0.044) and lower body weight (r = 0.624, P = 0.027) -two variables used to monitor the progression of CUS. Under normal physiological conditions, CRF orchestrates both hypothalamic and extra-hypothalamic stress responses. In that regard, CON rats with higher hippocampal and accumbal CRF expression showed hyperlocomotion (r = 0.656, P = 0.02) and higher sucrose consumption (r = −0.568, P = 0.043), respectively. These data coincide with the fact that under normal physiological conditions, CRF transcriptional activity is involved in psychomotor arousal and reward motivation, as suggested elsewhere 88 . After chronic stress, CRF has been linked to the induction of anxiety-and depression-like behaviours 89 . We found lower CRF expression after CUS, which may be inconsistent with a recent report showing higher CRF expression in cortical and subcortical regions associated with the extra-hypothalamic response to stress 90 . Such discordance may be dependent on the duration of the stress protocols: De Andrade et al. 90 found an increase after two weeks; we found a reduction at four weeks and others observed no significant differences after five weeks 58 . In this line, our data would correspond to a later compensatory phase after sustained CRF signalling, with CRF transcription returning to baseline levels over time in a region-dependent manner.
On the other hand, CUS induced no effects on serotoninergic and glutamatergic contents (Supplementary Table 6). Even though the 5-HT system plays an essential role in the neurobiology of neuropsychiatric disorders such as depression, at the preclinical level, there are almost as many reports showing no effects of CUS as those showing a decrease on the 5-HT contents or turnover 13,49 . As compared with monoamines, the role of glutamate concentrations in depression has been less investigated both at preclinical and clinical levels. For example, one study found higher glutamate concentration in the frontal cortex and the hippocampus at 24 hours, but not four weeks after CUS 91,92 . In general, discrepancies can be due to differences in the time elapsed from the completion of the CUS protocol and the neurochemical measurements and to variations in the anatomical specificity of tissue samples among studies 13 .
Regarding GABA, CUS reduced its brain concentrations, specifically in the nucleus accumbens. Our findings coincide with the evidence that lower GABA concentrations have been consistently associated with depression and anxiety symptoms in adolescents and adults [93][94][95] . Depressive individuals resistant to conventional antidepressant drugs, as well as subjects exposed to traumatic stress experiences, have the lowest GABA levels in the brain 93,95 . In this regard, sleep deprivation and chronic stress have been found to induce hyper-arousal and insomnia via reduction of GABA concentrations 96 . Regarding norepinephrine and dopamine, we found that CUS reduced both the norepinephrine contents in the nucleus accumbens and the dopamine turnover in the hippocampus. These data agree with several clinical, experimental, and pre-clinical studies, which have associated the reduction of catecholamines contents and signalling with symptoms of depression, anxiety, and other stress-related disorders 21,[97][98][99] . The neurochemical changes produced by stress were somewhat related to the behavioural alterations it induced. For instance, CUS animals with a reduced dopamine turnover in hippocampus (r = 0.757, P = 0.006) and lower concentrations of norepinephrine (r = 0.640, P = 0.023) and GABA (r = 0.863, P = 0.001) in nucleus accumbens and hippocampus (r = 0.760, P = 0.005) spent less time self-grooming in the OFT, and high levels of this behaviour have been repeatedly associated with stress de-arousal and OFT habituation 27,100,101 . In addition, CUS animals with lower accumbal dopamine turnover drank an even higher amount of sucrose (r = −0.626, P = 0.026) and had the lowest body weights (r = 0.660, P = 0.019) -two parameters affected by CUS in the present study. In general, the significant correlations suggest that the alterations induced by CUS were consistent and integrated well between the behavioural and neurochemical parameters (e.g., gene expression and neurotransmitter concentrations) measured here. We are aware that other preclinical studies have found contradictory results for the same neurotransmitters and in the same brain regions reported here 13 , so these data should be interpreted with caution. These correlations may not generalise to other studies where the pattern of behavioural results differs from that described here or to CUS studies in which stressed rats were housed in groups or control rats were housed in SI. exploratory factor analysis (efA) of behavioural paradigms. When judging the relevance of our behavioural results, it becomes evident that only one or two variables per test served well to discriminate between CON and CUS conditions, with the rest of the variables measured being weakly affected by stress and, thus, being barely informative if analysed separately. The same impression is obtained by examining the literature discussed in previous sections. Such inconsistencies may result from the combination of the following factors: a) implementing the CUS model at ages with less translational potential, b) not analysing the same variables over the studies, c) not considering the interrelationships among variables or analysing them separately, d) and not reporting indexes of size effects for judging their particular contributions to the whole effect. This situation often impedes having a clear take-a-home-message about which variables are more responsive to capture the impact of stress. EFA was then implemented to summarise the effect of CUS by reducing the data set to a more manageable number of variables while retaining the intrinsic behavioural meaning of each test. Once the latent variables were extracted, we compared them between groups to determine if they differentiate by CUS. Subsequently, those variables that were significantly affected by CUS were further reduced to obtain general behavioural domains. Finally, the odds ratios for belonging to the CUS group were estimated according to the predictive capacity of all significant latent variables. Please note that for all these analyses having completely orthogonal conditions is rather ideal, namely, comparing totally unstressed rats (i.e., group-housed controls) versus stressed rats (i.e., SI rats submitted to CUS stressors). Table 1, EFA computed with all SPT variables extracted two factors explaining 52% and 30% of the total variance, respectively. The KMO was in the limit of adequacy, but the BST www.nature.com/scientificreports www.nature.com/scientificreports/ indicated that the EFA was suitable for the SPT. The first factor retained the three sucrose preference measurements with factorial loads ranging from 0.95 to 0.79. As the preference of sucrose depends on the animal's capability to respond to natural rewards, this factor was then named as 'Reward sensitivity' . The second factor retained the two last sucrose consumption measurements with factorial loads ranging from 0.94 to 0.91, and it was consequently called 'Sucrose Consumption' . For each factor, all observable variables were positively associated with the latent variable. Loading coefficients of water consumption variables were lower than the cut-off criterion (i.e., 0.40) and, therefore, were not retained on any factor. When comparing the factor scores between groups, the reward sensitivity was significantly higher in CUS than in CON rats (P = 0.009, η 2 = 0.325). The sucrose consumption, on the contrary, did not differ significantly between groups but was descriptively higher in CUS rats.

Sucrose preference test (SPT). As shown in
Open field test (OFT). EFA was computed with all variables measured for the third OFT at PND 63 (see Table 2). Two factors explained 51% and 19% of the total variance extracted. The OFT was suitable for EFA as indicated by an acceptable KMO value with a highly significant BST. The first factor retained four variables with factorial loads ranging from 0.95 to 0.52. All variables related to the level and spatial distribution of the locomotor activity in the OFT (e.g., distance travelled and entries to and time spent in the central area) and were positively associated with the latent variable. Consequently, we named this factor as ' Ambulatory activity' . The second factor retained four variables, including vertical exploration (i.e., time and frequency of rearing), grooming time, and distance travelled. The factorial loads ranged from 0.95 to −0.40. As expected, grooming time was negatively associated with the latent variable 100 , which was named as 'Risk-assessment' in agreement with previous factorisations including almost the same behaviours 102 . When comparing the factor scores between groups, we found significantly high levels of ambulatory activity in CUS rats (P = 0.048, η 2 = 0.199), without differences on the risk assessment factor, which was also descriptively higher in stressed animals. Table 3, EFA computed with all EPM variables extracted three factors explaining 40%, 25%, and 9% of the total variance, respectively. The KMO was in the limit of adequacy, but with a highly significant BST. The first factor retained both traditional (i.e., time spent in and entries to the arms) and non-traditional (i.e., SAP and HD) anxiety-like parameters with factorial loads ranging from 0.84 to 0.46. Distance, time, and entries to the open arms and the central area were positively associated with the latent variable, whereas the time in the closed arms and SAP were negatively associated with it. Consequently, we called this factor ' Anxiolytic-like response' . All these parameters are related to the avoidance-approach conflict between novelty salience and the aversiveness induced by the different maze compartments 103 . The second factor retained all variables related to vertical and horizontal exploration within the maze, with factorial loads ranging from 0.84 to 0.51. As all observable variables were positively associated with the latent variable, we named this factor as 'General exploratory activity' , in agreement with a previous EFA description of a latent factor clustering these variables 102 . The third factor retained four variables, with factorial loads ranging from 0. 84  www.nature.com/scientificreports www.nature.com/scientificreports/ included grooming frequency and time, which were negatively associated with the latent variable. Also, the factor contained the time spent rearing and in the central area, which showed positive coefficients. The appearance of grooming in adverse situations may denote the activation of a de-arousal system as the stress situation is being overcome 100,101 . In contrast, rearing and time spent in the centre are activities aimed at vigilance and risk assessment 102 . Hence, we named this factor as 'Emotional distress' . When comparing the factor scores between groups, we found significantly lower levels of anxiolytic-like response (P = 0.043, η 2 = 0.209) and significantly higher levels of emotional distress in stressed rats (P = 0.001, η 2 = 0.266). The general exploratory activity was descriptively higher in CUS rats without reaching the significance level.

Elevated plus maze (EPM). As shown in
Spontaneous activity in object recognition arena. The ORT was carried out in an arena 38% larger than the OFT. This size difference was enough to elicit exploratory and defensive behaviours, although rats were already familiarised to the OFT. Thus, we measured the same OFT parameters and in each of the three ORT phases (i.e., habituation, trial, and test). To capture the common variance related to the spontaneous activity in the arena, we first computed one EFA only with these variables (see Table 4). Two factors explaining 43% and 17% of the total variance were extracted. The ORT was in the limit of adequacy with a moderate KMO coefficient and with a rather significant BST. The first factor retained six variables with factorial loads ranging from 0.91 to −0.70. The variables were time and frequency of rearing and grooming displayed during the habituation and trial sessions of the ORT. As in the OFT, the rearing and grooming were positively and negatively associated with the latent variable, respectively. We then named this factor as 'Risk-assessment' . The second factor included the distance travelled during the habituation, trial, and test phases of the ORT and was, therefore, named as ' Ambulatory activity' . The factorial loads ranged from 0.81 to 0.71. When comparing the factor scores between groups, CUS rats showed significantly higher levels of ambulatory activity (P = 0.002, η 2 = 0.441), without differences in the risk-assessment scores.
Object recognition test. The second EFA included all variables traditionally associated with the ORT (see Table 5). As with the previous EFA, the KMO was in the limit of adequacy with a rather significant BST. EFA extracted two factors explaining 40% and 26% of the total variance. The first factor retained six variables with factorial loads ranging from 0.81 to 0.61. The factor included the time and frequency of exploration of both objects during the trial and test phases, with all variables being positively associated with the latent variable. For this reason, we named the factor as 'Object exploration' . The second factor retained four variables with factorial loads ranging from 0.97 to −0.67. All variables related to the exploration and memory of the changed objects. The variables of object exploration corresponded to the trial session and were inversely associated with the latent variable, which included two memory discrimination variables. Accordingly, we named the factor as 'Object recognition' . When comparing the factor scores between groups, object exploration was significantly higher in CUS than in CON rats (P = 0.001, η 2 = 0.448), whereas object recognition was unaffected by stress.
Forced swimming test. The last test analysed was the FST (see Table 6). The KMO was in the limit adequacy, but with a highly significant BST. The EFA extracted two factors explaining 59% and 18% of the total variance. The  Table 3. www.nature.com/scientificreports www.nature.com/scientificreports/ factors retained the same behaviours (i.e., immobility, swimming, and climbing) but separated them by session, with factorial loads ranging from −0.91 to 0.68. The first factor retained the behaviours of the test session. As discussed above, the response to the test involves the so-called behavioural despair learning, in which animals learned from the pre-test that the stress is unavoidable and that no active coping responses will be useful, leading to an immobile posture. As the immobility behaviour correlated negatively with the latent variable, we named this factor as ' Antidepressive-like response' .
The second factor corresponded to the pre-test, which captures the unconditioned response to a strong uncontrollable and inescapable stress 85,104 . There, animals display active behaviours like swimming and climbing for long periods trying to escape. As long as time passes, immobility gradually displaces struggling activity. Again, the immobility behaviour correlated negatively with the latent variable and, therefore, we called this factor 'Stress-coping response' . When comparing the factor scores between groups, no differences in the antidepressant-like response were found, whereas the stress-coping response was significantly lower in stressed rats (P = 0.044, η 2 = 0.206).
Discussion. When assessing the validity of a given preclinical model, behavioural paradigms are chosen and implemented assuming they resemble -as a whole and not any single parameter-some characteristics of a given psychological or psychiatric construct. However, when it comes to presenting the effects of treatments and experimental conditions, the statistical analyses are usually applied to discrete variables separately, without considering the total variance derived from the interaction of all parameters measured. Here, we first used EFA to reduce a large number of behavioural parameters per test to obtain a more meaningful and manageable number of constructs. We were aware that we have a small sample size with many parameters. However, when using simulations or real-life data to investigate the effect of different participant-to-variable ratios, it has been concluded that changes in this ratio made little difference to the stability of EFA solutions 105 . The latter is particularly true when the EFA has extracted few latent variables explaining a considerable proportion of the variance with factors containing high loading coefficients. If factors have four or more loading coefficients greater than 0.6, then they are reliable regardless of the sample size and, therefore, increasing the number of subjects is very unlikely to worsen the outcome 106 . In all our EFA analyses those criteria were met, especially for the latent variables that were   significantly different between treatments. Despite the small sample size, the KMO values were over the limit of acceptance. This lower participant-to-variable ratio neither interfered with obtaining rotated factorial solutions, which were theoretically coherent for all behavioural paradigms. In fact, the analyses of the common variance shared by all observable variables were effectively reduced to few factors within each test. It is worth noting that EFA was computed without restricting the number of factors to be retained. We have already employed EFA to reduce behavioural paradigms in small sample sizes with quite successful results 27 .
The analysis of the object exploration and ambulatory activity displayed in the ORT included variables from three independent sessions and provided the two greatest size effects (e.g., 44-45%) when comparing CON and CUS groups. The other related factor was the ambulatory activity in the OFT, with an eta squared coefficient of 20%. All these latent variables were substantially increased by CUS, despite the repeated experience animals had with similar contexts throughout the experiment. Ambulation and exploratory activity are driven by the necessity of gathering information about the likely threat sources and the opportunities to escape or hide 66,101 . Thus, the persistence of these behavioural processes relates to a possible impairment in the arousal-inhibition and information-processing systems, leading to a heightened emotionality 100,107 . The second most important factor was the reward sensitivity from the SPT, increased by 33% by CUS (Fig. 6). Here, the shared variance between sucrose preference measurements was higher and more consistent than that for the sucrose consumption, even though consumption was also significantly different between groups. The latter is not surprising because preference already comprised both sucrose and water intake and was also corrected by body weight. In any case, these data suggest that hedonic or motivational components of reward were increased after CUS. It has been already shown that CUS substantially increased the intracerebral stimulation threshold of ventral tegmental area 108 . This evidence is consistent with the idea of CUS rats drinking significantly more sucrose to barely experience the same rewarding properties of the sweetened solution that would otherwise experience the unstressed counterparts. A link between sucrose consumption and rewarding/motivational processes has been previously established by the study of individual differences in sucrose feeding 109,110 . The EPM provided two factors with moderate effect sizes. The EPM parameters analysed separately did not differ between groups (Fig. 4D-F). However, as these observable variables were highly interdependent (e.g., time spent in closed and open arms), the latent variable did capture the effect of CUS that otherwise could not be detected. Finally, the FST contributed with one factor (i.e., stress-coping response) with an eta squared coefficient of 21%. Figure 6 summarised EFA factors that were significantly different between groups and ranked their effect sizes (eta squared coefficients) from the smallest to the largest one (i.e., 20% to 45%). In general, there were three latent variables related with one behavioural domain, including locomotor activity, exploration, and novelty habituation and other three latent variables belonging to a second behavioural domain comprising unconditioned anxiety, stress-coping response, and reward sensitivity. To summarise these two general domains, an additional factor analysis with all these factors was performed. As expected, EFA extracted two factors explaining 45% and 20% of the total variance ( Table 7). The first and most important factor retained the object exploration and ambulatory activity both in the ORT and OFT, with factorial loads ranging from 0.89 to 0.85. In a previous EFA, variables related to locomotion and exploratory activity also loaded into the same factor 27 . Thus, due to the nature of these variables and considering that all were positively interrelated, we named the first factor as 'Hyper-responsivity to novelty and mild threats' (Fig. 6B). The second factor retained the anxiolytic-like response, the stress-coping response, and the reward sensitivity factors from the EPM, the FST, and the SPT, respectively. The factorial loads ranged from −0.81 to 0.57. Similar to previous findings 27 , the anxiety-related parameters loaded together with the SPT and the FST into the same factor, supporting the consistency and the interpretation of our current factors. The reward sensitivity was positively associated with the latent variable, and considering that CUS increased sucrose preference, we named this factor as ' Anxiety/depressive-like response' (Fig. 6B). When comparing the factors between groups, CUS rats showed significantly higher scores on the first (P = 0.008, η 2 = 0.328) and second (P = 0.003, η 2 = 0.388) factor (Fig. 6B). These data agree with our previous study in which chronic SI differed from grouped-housing when comparing their scores on two main behavioural domains related to depression and unconditioned anxiety 27 . Finally, the eta squared coefficients shown in Fig. 6A provided us with information on how much those latent variables were affected by CUS separately, but do not allow judging their relative contribution to predicting the CUS condition when comparing them to each other. The binomial logistic regression showed that the best predictor of CUS was the latent variable 'object exploration' with an odds ratio of 26.68 (P = 0.03, AIC = 18.8). Based on the forward likelihood ratio method, only this factor was retained in  Table 6. www.nature.com/scientificreports www.nature.com/scientificreports/ the equation. From the excluded factors, the emotional distress (odd ratio = 5.17, P = 0.04, AIC = 25.1) and the reward sensitivity predictors (odd ratio = 5.13, P = 0.03, AIC = 24.1) showed much lower yet significant odds ratios. Regarding the two general behavioural domains, they showed rather similar predictive capacities (hyper-responsivity to novelty and mild threats: odd ratio = 11.42, P = 0.046, AIC = 22.7; anxiety/depressive-like response: odd ratio = 12.43, P = 0.042, AIC = 21.4).

conclusions
The CUS paradigm is one of the most used preclinical models in behavioural neuroscience. Surprisingly, the majority of CUS studies has been conducted on adult rats, despite the well-known evidence that risk factors for mood disorders accumulate during the first two decades of life. Most of these studies report few variables from only one or two behavioural tests -a situation that has led to a vast array of inconsistencies. The purpose of the present study was, therefore, to assess the validity of the CUS model by including ethologically relevant behavioural paradigms to characterise the effects of chronic stress in juvenile rats. Traditionally, the reduction of sucrose consumption/preference after CUS has been considered a marker of anhedonia -a core symptom of depressive disorders. The ability of this measure to discriminate the effects of CUS was compared to other parameters, such as body weight and locomotor activity. For this purpose, two independent CUS experiments were carried out in post-weaning male rats. Experiment 1 included a milder CUS protocol, which was compared with SI and group-housed controls. There, CUS reduced body weight and SI increased sucrose consumption and preference. When correcting the intake by body weight, the sucrose consumption in CUS rats was now similar to that in SI and, in consequence, higher than in controls. As we were unable to reproduce some of the expected CUS alterations, we made the protocol more severe to better characterise the effects of CUS. In experiment 2, body weight was, again, the most responsive parameters to monitor the progression of CUS, followed by the locomotor activity. The weakest parameters belonged to the SPT, with sucrose consumption being increased instead of decreased by CUS after correcting by body weight. The behavioural characterisation of CUS showed that only one or two variables per test served well to discriminate between groups. Again, locomotor and exploratory activity were the most consistent parameters. Although we acknowledge that our neurochemical data may not generalise to other studies where the pattern of behavioural results differs from that described here, our CUS   www.nature.com/scientificreports www.nature.com/scientificreports/ protocol was found to affect brain parameters typically associated with neural plasticity, anxiety, and incentive motivation. By employing the EFA, we were able to reduce each test to a few but meaningful numbers of latent variables, which could be successfully clustered into two general behavioural domains. The first one relates to the hyper-responsivity to novelty and mild threats (including locomotor activity, exploration, and novelty habituation factors) and the second one deals with the anxiety/depressive-like response (including unconditioned anxiety, stress-coping response, and reward sensitivity). The differences obtained when comparing observable and latent variables support our point that in most cases no individual parameter is robust enough to capture the complexity of a paradigm, especially when revisiting the consistency and construct validity of a model. The latter was best illustrated with the SPT, which was the third most important latent factor but the worst observable variables when analysed separately. Altogether, the analyses applied to both observable and latent variables suggest that stress during adolescence impairs the arousal-inhibition system leading to an augmented and persistent response towards novel, rewarding, and mildly threatening stimuli, accompanied by lower body-weight gain. Finally, we recommend the inclusion of relevant behavioural parameters to increase the robustness and reliability of a given model in combination with appropriate analytic techniques to uncovering the complex underlying structure of the data. As we did not include a condition of CUS rats housed in groups or control rats housed in SI (experiment 2), we were unable to separate out the specific contribution of SI from the other CUS stressors. However, to characterise the neurobehavioural effects of CUS, we still find of value comparing totally unstressed rats (i.e., group-housed controls) versus stressed rats (i.e., SI rats submitted to CUS stressors).

Materials and Methods
Animals. Behavioural tests. The Sucrose Preference Test (SPT) was carried out at PNDs 42, 52, and 62 ( Fig. 1A) similar as previously reported 24 . Rats were housed individually in standard polycarbonate cages for 24 hours with ad libitum access to food and two bottles: one containing 200 ml of 1% sucrose solution (w/v) and another filled with 200 ml of tap water. Sucrose and water consumption was first computed as percentages [(intake ml/200 ml)x100] and then corrected by body weight (consumption%/body weight in grams). Similarly, preference was calculated as percentages [(sucrose consumption/(sucrose + water consumption)) X 100] and corrected by body weight. The Open Field Test (OFT) took place at PNDs 43, 53 and 63 (Fig. 1A) as previously reported 101 . During the 10-minutes session the following variables were measured: the total distance travelled, the total distance travelled in the central area, time (in seconds) spent in and the number of entries to the centre, the time and frequency of rearing (posture sustained with hind paws on the floor) and grooming (including washing or mouthing of forelimbs, hind paws, face, body and genitals). The Elevated Plus Maze (EPM) was carried out at PND 65, similarly as previously reported 100 . During the 5-minutes session the following parameters were scored: distance travelled in, time spent on, and the number of entries (i.e., four paws being placed on the area) to the open arms, closed arms and central area; as well as the total distance travelled in the maze, time and frequency of rearing, grooming, head-dipping (HD), and stretch-attempt-posture (SAP). Object Recognition Test (ORT) was performed similarly as recently reported 111 . The test consisted of three 5-minutes sessions: habituation, sample trial, and test trial. On the habituation session, animals were allowed to explore and habituate to the empty open field. Twenty-four hours later, animals were exposed to the arena containing two identical objects placed in the back corners of the box, situated 15 cm away from the walls: either two red iron cylinders (5 cm in diameter, 8 cm high) or two blue iron pillars (6 cm in diameter, 8 cm high). During the inter-trial interval, animals were returned to their home cages for 30 minutes. Afterwards, rats were again allowed to explore the objects in the test trial, but now one of the familiar objects was changed by a novel one. Objects and spatial location were counterbalanced across groups. Object exploration was scored whenever the rat's nose touched the object or when it was directed toward the object within a distance of 2 cm. Climbing onto an object was not recorded as exploration unless the snout of the rat was directed towards it by less than 2 cm. For all sessions, the following variables were measured when corresponded: familiar and novel object exploration time and frequency, total distance travelled, rearing (time and frequency), and grooming (time and frequency). The objects' exploration time during test trial was normalised by time spent exploring the same object/position during sample trial (NET) as follows: (object exploration time in the test trial in seconds)/(object exploration time in the sample trial in seconds). The Discrimination Index was also computed as [(NET novel object-NET familiar object)/(NET novel object + NET familiar object)]. Such parameter was interpreted as a memory index. Forced Swimming Test (FST) was carried out as previously reported 20 and consisted of two sessions on two consecutive days: a 15-minutes pre-test and a 5-minutes test (at PNDs 68 and 69, respectively). The following parameters were observed: duration of immobility, swimming and climbing, and total time spent in active behaviours (swimming time + climbing time). Only the first 5 minutes of the pre-test were scored to get a fair comparison between the test and pre-test sessions scored. Etholog 2.25 112 and ANY-maze (version 4.72, Stoelting Co., USA) software packages were used for manual and automatic behavioural analysis, respectively. All apparatuses were cleaned with 70% ethanol between animals and sessions. Other details regarding behavioural testing and apparatuses are included in the Supplementary methods. www.nature.com/scientificreports www.nature.com/scientificreports/ Analysis of gene expression and ex-vivo neurotransmitters contents. Animals were euthanised by decapitation once the CUS protocol finished (PND 65). Brains were quickly dissected on ice, and three different areas were collected: The medial prefrontal cortex, the hippocampus, and the nucleus accumbens. Both hemispheres were pooled in the case of the medial prefrontal Cortex, whereas for the hippocampus, and the nucleus accumbens only one sample per hemisphere was used following a right-and-left alternating method. For RT-qPCR, all samples were run in duplicates, and their mean values were used for further calculations. Each run included both CON and CUS group samples. Furthermore, each gene was run individually according to the sample maximisation method 113 . Non-template controls and minus RT controls were included in each run. The absence of amplification in the non-template and the minus RT controls excluded the possibility of genomic DNA contamination. Fluorescence data were collected, and the threshold cycle (Ct) was calculated using the Rotor-Gene Q Series Software (QIAgen, Germany). The remaining samples from the hippocampus and the nucleus accumbens were used for neurochemical analyses. High-performance liquid chromatography coupled with electrochemical detection (HPLC-EC) was used, and all procedures were carried out as previously reported 24 . All samples were analysed for their contents on norepinephrine (NE), dopamine (DA) and its metabolite 3,4-dihydroxyphenylacetic acid (DOPAC), serotonin (5-HT) and its metabolite 5-hydroxyindoleacetic acid (5-HIAA), using the internal standard method. Furthermore, the DA and 5-HT turnover were also calculated (DOPAC/DA and 5-HIAA/5-HT, respectively). Glutamate (Glu) and gamma-aminobutyric acid (GABA) were analysed by reverse-phase HPLC with fluorescence detection (HPLC-FD) (Agilent Technologies, USA). The amino acid concentration was determined using the peak area and the external standard method. Data for both monoamines and amino acids concentrations were expressed as nanograms per milligram of wet tissue weight.
Statistical analysis. Data were presented as mean ± standard error of the mean (SEM). Only significance (P) and eta squared (η 2 ) values were shown in the main text (all statistical details are shown in Supplementary results). In experiment 1, the between-groups differences were estimated following planned contrasts. The planned contrasts were specifically designed for each parameter after a detailed, visual inspection of the data and following a priori hypotheses regarding the likely effects of the treatments. Similarly, as with posthoc test, planned contrasts were corrected to avoid committing a type I error. However, instead of comparing all possible pairwise comparisons, planned contrasts were limited to a restrict number of comparisons (i.e., k-1), with "k" being the number of groups. The t-tests and p-values from 1000 bootstrap samples were reported for each contrast following the default function displayed by SPSS package (compare means/one-way ANOVA/contrasts/ bootstrap) (for details see Supplementary results). In experiment 2, the body weight, SPT, and OFT were analysed using mixed analysis of variance (ANOVA) with Session as within-subjects factor with three levels (i.e., PND 42,PND 52,and PND 62) and Treatment as between-subjects factor with two levels (groups CON and CUS). Additionally, the OFT carried at PND 63 had Minutes as within-subjects factor with ten levels (i.e., from minute 1 to minute 10). For the EPM and the OFT carried out as the habituation session of the ORT, two factors were used: Minutes as within-subjects factor (five levels: from minute 1 to minute 5), and Treatment as between-subjects factor. In the case of the ORT, a comparison of Trials (i.e., sample vs test) as within-subject factor and Objects (i.e., familiar vs novel) and Treatments as between-subject factors were computed. In the FST, the within-subject factor was Session with two levels (i.e., pre-test and test) and Treatment as between-subjects factor. For analysis of gene expression in the brain samples, the Region was included as a between-subject factor with three levels (i.e., mPFC, HPC, and NAc), in addition to Treatment. Likewise, the comparison of the ex-vivo concentration of neurotransmitters included the Region with two levels (i.e., HPC and NAc) and the Treatment, both as between-subject factors. For all analyses of brain samples, the Mother was included as an additional between-subject factor. Based on our previous studies 114, 115 and preliminary results (data not shown), we identify that for genetic analysis the variance shared by littermates may sometimes exceed the variance between groups, which may eventually increase the likelihood of committing error type I or II. The total sample consisted of 40 rats from 9 mothers, which were balanced between groups to the largest extent possible. In the subgroup of 20 rats selected for brain analysis, the number of littermates per group was not as equally balanced as with the 40 rats, because other variables (i.e., the sucrose preference, locomotion, and body weight) were also incorporated to split the samples. The latter justified, even more, the use of the mother as a variable in the analysis. For all repeated-measures analyses, the Greenhouse-Geisser correction was used when appropriate. A p-value < 0.05 was considered statistically significant.
A multiple linear regression analysis was conducted with the dependent variable groups categorised with dummy codes of "0" for controls and "1" for CUS -the condition to be predicted in the model. The first analysis included the variables measured within the 30 days of CUS. We averaged the data of body weight, sucrose consumption and preference, and locomotor activity of PND 42, PND 52, and PND 62 to obtain only one variable (i.e., a general mean) per parameter. Thus, the regression analysis was run using groups as a dependent variable and body weight, sucrose consumption and preference, and locomotor activity as predictors. The second analysis included as predictors only the behaviours that were significantly different between groups during the behavioural characterisation phase. We performed the analysis, first, using the stepwise mode of the SPSS package. Once a solution was obtained, this stepwise model was manually confirmed by including and removing all the variables in all possible combinations on successive blocks of analysis using the statistical package Jamovi (Jamovi project 2018, Version 0.9.5.12, retrieved from https://www.jamovi.org). For all analyses, the adjusted regression coefficient was reported (adjusted R 2 ). The Jamovi package also displayed the Akaike information criterion (AIC) parameter, which was used to compare and select the best model. The best model was the one with the lowest AIC value, and a delta (Δ) of AIC values between models was computed to choose the model. The simplest model consisted of the dependent variable "groups" and one predictor, and on subsequent models, more predictors were added as long as they produced significant contributions to the whole prediction.
www.nature.com/scientificreports www.nature.com/scientificreports/ The exploratory factor analysis (EFA) was carried out using the observable variables (e.g., behavioural parameters) of each behavioural test. The EFA was conducted as previously reported 27 . Briefly, the factorial solution was made with Principal Components with orthogonal Varimax rotation. Kaiser-Meyer-Olkin test (KMO) is a measure of sample adequacy and was calculated to determine the sample adjustment. KMO values from 0 to 0.49 indicate that the data are unsuitable for EFA. As the KMO may be sensitive to rather small sample sizes, it would yield low coefficients even when the EFA is still appropriate. Thus, to complement the KMO test, Bartlett's Sphericity test (BST) was computed. The BST compares the observed correlation matrix to the identity matrix. If the null hypothesis -i.e. that both matrixes are equal-is rejected, then it is assumed that there are enough commonalities among variables that can be well summarised with a few numbers of factors 116 . Regarding the factorial loads, a value of ±1 indicates a perfect correlation of the variable with the component (factor). Loading values ranging from ±0.60 to ±0.80 indicate a strong correlation, values ranging from ±0.40 to ±0.60 indicate a moderate correlation and values lower than ±0.40, a poor correlation. Factorial loadings were therefore restricted to >0.40 to retain variables with at least moderate correlation to the factor. Positive component values indicated that a variable is directly related to the behavioural meaning of the factor, whereas a negative component value means the opposite. The factor scores for each latent variable extracted were computed using the regression method included by default in the SPSS package. These scores are standardised regression coefficients with mean 0 and variance 1. In this method, the factor loadings are adjusted to take account of the initial correlations between variables, with the advantage that differences in scaling between variables are stabilised 116 . Then, each animal received an individual score indicating how it performed on each factor extracted, which allowed us to further compare the group differences per construct or latent variable by employing an ANOVA analysis. The variables that were significantly different between groups were further reduced to obtain general behavioural domains. Finally, a binomial logistic regression analysis was conducted to estimate the odds ratios to belong to the CUS group according to the predictive capacity of all significant latent variables. To estimate the particular contribution of each predictor, the forward method with the likelihood ratio statistic was employed using the SPSS. With the Jamovi package, the variables were tested on different combinations in successive blocks. For each predictor, the significance of the Wald test, the odds ratios, and the AIC scores were reported.

Data availability
All datasets analysed in the current study will be made available upon reasonable request. All data generated from this work other than those presented in figures tables are available in the Supplemental Information file.