# Perinatal Stress Programs Sex Differences in the Behavioral and Molecular Chronobiological Profile of Rats Maintained Under a 12-h Light-Dark Cycle

Sara Morley-Fletcher, Jerome Mairesse, Gilles van Camp, Marie-Line Reynaert, Eleonora Gatta, Jordan Marrocco, Hammou Bouwalerh,

Ferdinando Nicoletti, Stefania Maccari

► To cite this version:

Sara Morley-Fletcher, Jerome Mairesse, Gilles van Camp, Marie-Line Reynaert, Eleonora Gatta, et al. Perinatal Stress Programs Sex Differences in the Behavioral and Molecular Chronobiological Profile of Rats Maintained Under a 12-h Light-Dark Cycle. Frontiers in Molecular Neuroscience, 2019, 12, 10.3389/fnmol.2019.00089. [hal-03094662](https://hal-univ-lille.fr/hal-03094662)

HAL Id: hal-03094662

https://hal.univ-lille.fr/hal-03094662

Submitted on 7 Jan 2021

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

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

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

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

Check for
updates

# Perinatal Stress Programs Sex Differences in the Behavioral and Molecular Chronobiological Profile of Rats Maintained Under a 12-h Light-Dark Cycle

Sara Morley-Fletcher<fnref n="1" />,²†, Jerome Mairesse³†, Gilles Van Camp<fnref n="1" />,², Marie-Line Reynaert<fnref n="1" />, Eleonora Gatta⁴, Jordan Marrocco⁵, Hammou Bouwalerh<fnref n="1" />,², Ferdinando Nicoletti⁶,⁷† and Stefania Maccari<fnref n="1" />,⁸*‡

**OPEN ACCESS**

**Edited by:**
Gi Hoon Son,
Korea University, South Korea

**Reviewed by:**
Henrik Oster,
Universität zu Lübeck, Germany
Han Kyoung Choe,
Daegu Gyeongbuk Institute
of Science and Technology (DGIST),
South Korea

*Correspondence:
Stefania Maccari
stefania.maccari@univ-lille.fr

† These authors have contributed
equally to this work

‡ Present address:
Ferdinando Nicoletti and
Stefania Maccari,
Co-director of International
Associated Laboratory (LIA),
University of Lille 1 - CNRS
and Sapienza University of Rome -
IRCCS Neuromed, Rome, Italy

Received: 28 November 2018

Accepted: 21 March 2019

Published: 01 May 2019

**Citation:**
Morley-Fletcher S, Mairesse J,
Van Camp G, Reynaert M-L, Gatta E,
Marrocco J, Bouwalerh H, Nicoletti F
and Maccari S (2019) Perinatal Stress
Programs Sex Differences
in the Behavioral and Molecular
Chronobiological Profile of Rats
Maintained Under a 12-h Light-Dark
Cycle. *Front. Mol. Neurosci.* 12:89.
doi: 10.3389/fnmol.2019.00089

Stress and the circadian systems play a major role in an organism's adaptation to environmental changes. The adaptive value of the stress system is reactive while that of the circadian system is predictive. Dysfunctions in these two systems may account for many clinically relevant disorders. Despite the evidence that interindividual differences in stress sensitivity and in the functioning of the circadian system are related, there is limited integrated research on these topics. Moreover, sex differences in these systems are poorly investigated. We used the perinatal stress (PRS) rat model, a well-characterized model of maladaptive programming of reactive and predictive adaptation, to monitor the running wheel behavior in male and female adult PRS rats, under a normal light/dark cycle as well as in response to a chronobiological stressor (6-h phase advance/shift). We then analyzed across different time points the expression of genes involved in circadian clocks, stress response, signaling, and glucose metabolism regulation in the suprachiasmatic nucleus (SCN). In the unstressed control group, we found a sex-specific profile that was either enhanced or inverted by PRS. Also, PRS disrupted circadian wheel-running behavior by inducing a phase advance in the activity of males and hypoactivity in females and increased vulnerability to chronobiological stress in both sexes. We also observed oscillations of several genes in the SCN of the unstressed group in both sexes. PRS affected males to greater extent than females, with PRS males displaying a pattern similar to unstressed females. Altogether, our findings provide evidence for a specific profile of dysmasculinization induced by PRS at the behavioral and molecular level, thus advocating the necessity to include sex as a biological variable to study the set-up of circadian system in animal models.

**Keywords:** predictive adaptation, reactive adaptation, circadian rhythms, locomotor activity, chronobiological stressor, mRNA expression

<footnote n="1">¹ UMR 8576, Unité de Glycobiologie Structurale et Fonctionnelle, Campus Cité Scientifique, CNRS, University of Lille, Lille, France, ² University Lille - CNRS-UMR 8576, International Associated Laboratory (LIA) "Prenatal Stress and Neurodegenerative Diseases," Sapienza University of Rome - IRCCS Neuromed, Rome, Italy, ³ Division of Neonatology, Department of Pediatrics, University of Geneva, Geneva, Switzerland, ⁴ Department of Psychiatry, College of Medicine, Psychiatric Institute, Center for Alcohol Research in Epigenetics, University of Illinois at Chicago, Chicago, IL, United States, ⁵ Laboratory of Neuroendocrinology, The Rockefeller University, New York, NY, United States, ⁶ Department of Physiology and Pharmacology "V. Erspamer," Sapienza University of Rome, Rome, Italy, ⁷ Istituto di Ricovero e Cura a Carattere Scientifico (IRCCS), NEUROMED, Pozzilli, Italy, ⁸ Department of Science and Medical - Surgical Biotechnology, Sapienza University of Rome, Rome, Italy</footnote>
INTRODUCTION

The optimal functioning of an organism is based on its ability to adapt to changes in the environment. The adaptation can be reactive to unexpected environmental changes, or predictive, as it prepares the organism to anticipate the daily changes in the environment (Moore-Ede, 1986). The glucocorticoid stress response is an essential mediator of allostatic processes (McEwen, 1998), which are engaged to allow the organism to find a new equilibrium in response to the specific stressor. This stress response is adaptive when it is transitory and shuts off quickly with the help of feedback loops (McEwen and Stellar, 1993; Maccari et al., 2003, 2017). However, in the long run, prolonged and/or uncontrollable sources of stress lead to a loss of organism's resistance and entering a phase of exhaustion that might account for many clinically relevant disorders at the behavioral, endocrine, metabolic, cardiovascular, and immune level (de Kloet, 1999; Russell et al., 2018). Glucocorticoids also have a circadian rhythm, and the circadian system is another major physiological system involved in the adaptation of the organism to environmental changes. In contrast to the stress system, whose adaptive value is reactive, the circadian system has a predictive adaptive role (Turek, 2016). As with the stress response, disruption of circadian rhythms leads to clinically relevant disorders including obesity, cardiovascular disease, inflammation, and cognitive impairments (Wirz-Justice and Van den Hoofdakker, 1999; Bunney and Potkin, 2008; Gerstner and Yin, 2010; Wulff et al., 2010; Medic et al., 2017). Circadian rhythms are under control of an internal clock present in nearly every cell, with a core mechanism governed by a genetic network. These clocks are organized hierarchically, with a central pacemaker located in the suprachiasmatic nuclei (SCN) of the hypothalamus at the top, that receive photic cues and in turn coordinate local circadian clocks in the periphery (Moore and Eichler, 1972; Klein et al., 1991; Buijs and Kalsbeek, 2001; Cirelli, 2009). Glucocorticoid hormones and body temperature, among other signals, act as internal systemic synchronizers to entrain peripheral clocks through the SCN (Dibner et al., 2010).

Stress-induced pathologies are often characterized by alterations in circadian rhythms (Weibel et al., 2002; Spiga et al., 2014; Liyanarachchi et al., 2017; Nicolaides et al., 2017), and the glucocorticoid stress response and many clock-controlled processes display sexual dimorphism (for a review, see Darnaudéry and Maccari, 2008; Mong et al., 2011; Heck and Handa, 2018). The corticosterone stress response is markedly greater in adult female rodents than in males (Heck and Handa, 2018), and in both rodents and humans the circadian rhythm of glucocorticoids secretion has higher amplitude in females (Mong et al., 2011; Vestergaard et al., 2017, respectively). In humans, women also display higher rates of sleep disorders than men and a differential ability to adapt to changing shift-work schedules and jet lag (Bailey and Silver, 2014). A sexual dimorphism and first demonstration of gonad-dependence of circadian rhythm concerns circadian locomotor activity in rodents (Albers et al., 1981; Wollnik and Turek, 1988). Gonadal steroids in the SCN seem to act

through different mechanisms in males and females. Indeed, androgen receptors are densely located in the core SCN and therefore act on the SCN via a direct pathway, unlike the estrogen receptors that are located mostly extra-SCN and act on SCN afferents (Mong et al., 2011). Nevertheless, the number of studies that associate reactive and predictive adaptations is limited, with the females being underrepresented (Zucker and Beery, 2010; Beery and Zucker, 2011; Kiryanova et al., 2017; McCarthy et al., 2017).

An individual's capabilities to cope with environmental challenges are shaped during development by the interplay of genetic and epigenetic determinants (McEwen, 2016). Exposure to stressful events in early life strongly programs an individual's phenotype and adaptive capabilities by modifying reactive (Weaver et al., 2004; Seckl, 2008; Brunton et al., 2011; Maccari et al., 2017) as well as predictive adaptation (Van Reeth et al., 2000; Maccari and Morley-Fletcher, 2007). A well-characterized animal model of early programming of stress response is the PRS model in the rat (Maccari et al., 1995; for a review see Maccari et al., 2017), in which exposure to gestational stress and altered maternal behavior programs a life-long disruption in the reactive adaptation such as a hyperactive response to stress and a defective feedback of the hypothalamus-pituitary-adrenal (HPA) axis (Vallée et al., 1999; Gatta et al., 2018) together with long-lasting modifications in stress/anti-stress gene expression balance in the hippocampus (Mairesse et al., 2015; Gatta et al., 2018). The changes induced by PRS in predictive adaptation include a profound alteration of the sleep-wake cycle architecture together with increased sleep fragmentation when the animals reach adulthood (Dugovic et al., 1999; Mairesse et al., 2013) and alterations in the rhythm of plasma corticosterone secretion (Koehl et al., 1999). Nevertheless, there are no studies on changes in circadian rhythms at the molecular level in the PRS model.

Perinatal stress induces sex-specific effects on reactive and predictive adaptations as well as on sexual maturation induced by gonadal hormones (Reynaert et al., 2016). Indeed, PRS reduces the ratio of testosterone/dihydrotestosterone in males and estradiol levels in females. Such an effect is causally related to sex-dependent alterations in gene expression in the hypothalamus and in the hedonic sensitivity to palatable food (Reynaert et al., 2016). However, studies have generally been conducted in PRS males and have focused on reactive adaptation to a greater extent.

The aim of this study was to explore the influence of PRS on the circadian oscillations of gene expression in the SCN and on circadian locomotor behavior, in a sex-dependent manner. Of note, research on transcriptional rhythms has shown that more than half of all genes in the human and rodent genome follow a circadian pattern (Panda et al., 2002; Hughes et al., 2009; Summa and Turek, 2011; Zhang et al., 2014). Thus, we focused on genes belonging to four functional classes, namely the circadian clock, HPA axis stress response regulation, signaling and glucose metabolism in male and female adult PRS rats. We also monitored the running wheel behavior, first under a regular light/dark (L/D) cycle and then after an abrupt 6-h advance light shift (chronobiological stress).
# MATERIALS AND METHODS

## Ethics Statement

All the experiments performed in this study followed the guideline of the European Communities Council Directive 2010/63/EU. The protocol of perinatal stress, running wheel behavior and sample collection was approved by the Local Committee CEEA-75 (Comité d'Ethique en Experimentation Animale Nord-Pas de Calais, 75).

## Animals and Perinatal Stress Procedure

### Animals

In total, 48 nulliparous female Sprague-Dawley rats, weighing approximately 250 g, were purchased from a commercial breeder (Harlan, France). Animals were housed at constant temperature (22 ± 2°C) and under a regular 12 h light/dark (L/D) cycle (light on at 8:00 am). A vaginal smear using endocrine serum (NaCl 0.9%) was performed on the morning following mating with an experienced male. The day on which the smear was sperm positive was considered to be embryonic day 0 (E0). After mating, pregnant females were individually housed with *ad libitum* access to food and water at constant temperature (22°C ± 2°C), and under a regular 12 h light/dark cycle (light on at 08:00). On E11, pregnant females were randomly assigned to stress or control groups (n = 24 in each group). Control females were left undisturbed, with an exception made for weighing one time per week in order to follow gestation.

### Perinatal Stress Procedure

The pregnant female group was subjected to a restraint stress procedure according to a standard protocol (Maccari et al., 1995; Morley-Fletcher et al., 2003). From day 11 of pregnancy until delivery, pregnant female rats were subjected to three stress sessions daily (45 min each), during which they were placed in plastic transparent cylinders with a conical extremity and exposed to bright light or were left undisturbed (control dams). Stress sessions were conducted during the light phase (between 09:00 and 15:00) with a minimum interval of 2 h between each stress session. The local ethical committee approved the gestational restraint procedure. Maternal behavior was monitored for 24 h every day during the first seven post-partum days. Constant monitoring was performed with small infrared cameras placed on the animal cage rack where cages containing lactating females were placed. Within each observation period, the behavior of each mother was scored every minute from post-partum day 1 to day 7 (60 observations/h with 2 h of observation per day, 1 h before lights off and 1 h after lights on). The active behavior of the mother (nursing behavior, grooming, licking, and carrying pups) was scored and the data obtained were expressed as percentages with respect to the total number of observations. Since gestational stress induces a reduction of maternal behavior (Gatta et al., 2018), we refer to the whole procedure as PRS having a prenatal and postnatal effect). In the present study, only male and female offspring from dams presenting a stress-reduced maternal behavior (with a cutoff below 40% of maternal care in the PRS group vs. a cutoff above 60% of maternal care in the control group), and from litters of 10-14 rats with a

similar number of males and females, were used. Weaning of the offspring occurred at 21 days after birth.

## Experimental Design

Separate sets of rats were used for behavioral studies (n = 40, 10 rat/sex/group) and transcriptomic analysis (n = 60, n = 5/group/time point). Within each litter, we took only two males and two females to minimize any litter effect (Chapman and Stern, 1979). Therefore, there were no more than two sibling within the same time point for the transcriptomic study. Housing conditions were two same-sex animals per cage from each perinatal and experimental group (CONT vs. PRS and behavioral or transcriptomic study) after weaning and until experiments started (at 4 months of age).

## Running Wheel Activity

A first group of rats (n = 10/sex/group) was used to analyze the rhythm of circadian activity under a regular 12/12 L/D cycle (08:00-20:00). Rats were housed in light, air tight chambers equipped with continuously operating ventilating fans and single-housed in individual cages equipped with a running wheel that allowed continuous recording of locomotor activity via an on-line computer (Chronobiology kit, Stanford Software System, CA, United States), with light intensity set at 30-40 lux at cage floor level. Each chamber contained 10 single-cages, and males and females were housed separately. During the course of the experiments, food and water were provided *ad libitum*; room temperature (22°C) and humidity (60%) were kept constant. After 10-15 days of adaptation to the running wheels, the rhythm of activity was individually analyzed over 10 consecutive days. The onset of activity was identified with a 5 min resolution and was defined as the first time point at which the mean intensity of activity was above 10% of the maximum and remained above that point for at least 50% of the time during the following 30 min. The reversed procedure was used for the cessation (offset) of activity (first time point below 10% of maximum and activity remained below that point for at least 50% of the time during the following 30 min). The time elapsed between the onset and offset of activity was defined as the total time of nocturnal activity, the peak value of activity and peak hour of activity were directly determined on the actogram for each animal. The mean 24 h integrated activity was determined by adding the mean number of revolutions in the wheel, every 5 min over 10 consecutive days for each animal. The data were then plotted with a 30 min or 5 min resolution; this represented the mean number of wheel revolutions performed by the animals.

## Exposure to a Chronobiological Stressor

At the end of the recording of normal activity, rats (n = 10/sex/group) were subjected to an abrupt 6-h advance shift in the L/D cycle. On the day of the shift, lights were turned-off 6 h before the current time, and the new 12 h L/D cycle (light on, 02:00) was maintained thereafter. The time to resynchronize to the new L/D cycle was defined as the number of days for the animal to exhibit a regular activity for at least three consecutive days. The time taken for re-entrainment to the new L/D cycle was individually assessed using the onset of activity criterion
determined as previously described. This criterion was defined as the smallest number of days required for the shifted activity onset to occur within 30 min of lights-off and to be stable for 3 days under the new L/D cycle.

## Gene Expression Analysis in the Suprachiasmatic Nucleus (SCN)

A separate set of male and female CONT and PRS rats was used for gene expression analysis in the brain (n = 5 per group/time point) at three selected time points (15:00, 19:00 and 02:00) under a regular 12/12 L/D cycle (with light on at 08:00) in order to extract SCN during light phase (15:00), before the onset of dark phase (19:00) and during it (02:00).

The SCN of male and female CONT and PRS rats was rapidly dissected and kept frozen at -80°C. Special attention was taken for tissue collection at 2 am: in order to avoid any influence of light on gene expression, animals were killed under red light before tissue dissection, which did occur under white light. RNA extraction was performed using the RNeasy Plus mini kit (Qiagen, France) following manufacturer's instructions. RNA concentration was determined using Nanodrop (ND-1000, Labtech, Germany), and quality verified by RIN (RNA Integrity Number; Bioanalyzer 2100, Agilent Technologies, France). All RNA samples had a quality score above 7.5 RIN. Retrotranscription was performed with the High-Capacity cDNA Reverse Transcription kit (Applied Biosystems, France). Transcript levels were measured by Custom Taqman qRT-PCR (Applied Biosystems, France). The following Taqman real-time PCR probes were obtained from Applied Biosystems:

### Clock-Related Genes

Clock (Rn00573120_m1), aryl hydrocarbon receptor nuclear translocator-like (Bmal, Arntl, Rn00577590_m1), Cry1 (Rn01503063), period circadian regulator 1 (Per1, Rn01496757_m1), period circadian regulator 2 (Per2, Rn01426757_m1), period circadian regulator 3 (Per3, Rn00709499_m1), nuclear receptor subfamily 1-group D-member 1 (Rev-erbA-α, Nrd1, Rn00595671_m1), melatonin-1A receptor (Mel1AR, Rn01488022_m1), melatonin-1B receptor (Mel1BR, Rn01447987_m1), neuronal PAS domain protein 2 (Npas2, Rn01438223_m1).

### HPA Axis Stress-Related Genes

Corticotropin releasing hormone receptor 1 (CrhR1, Rn00578611_m1), corticotropin releasing hormone receptor 2 (CrhR2, Rn00575617_m1), glucocorticoid receptor (GR, Nrc1, Rn00565562_m1), mineralocorticoid receptor (MR, Nrc2, Rn00561369_m1), Arginin vasopressin receptor 1a (Avpr1a, Rn00583910_m1), oxytocin receptor (OxtR, Rn00563503_m1).

### Signaling-Related Genes

Metabotropic glutamate receptor 2 (mGluR2, Grm2, Rn01447672_m1), metabotropic glutamate receptor 3 (mGluR3, Grm3, Rn01755352_m1), metabotropic glutamate receptor 5 (mGluR5, Grm5, Rn00566628_m1), 5-hydroxytryptamine receptor 2C (Htr2c, Rn00562748_m1), mechanistic target of rapamycin (mTor, Rn00693900_m1), v-akt murine thymoma

viral oncogene homolog 1 (Akt1, Rn00583646_m1), disrupted in schizophrenia 1 (Disc1, Rn00598264_m1), dickkopf WNT signaling pathway inhibitor 1 (Dkk1, Rn01501537_m1), early-growth transcription factor 1 (Ngf1A, Egr1, Rn00561139_m1), erb-b2 receptor tyrosine kinase 2 (Erbb2, Rn00566561_m1), erb-b2 receptor tyrosine kinase 3 (Erbb3, Rn00568107_m1), erb-b2 receptor tyrosine kinase 4 (Erbb4, Rn00572447_m1), solute carrier family 17-member 7, vesicle bound sodium-dependent inorganic phosphate cotransporter (vGlut1, Slc17a7, Rn01462431_m1), solute carrier family 1 member 3, glial high affinity glutamate transporter (Glast, Slc1a3, Rn00570130_m1).

### Glucose Metabolism-Related Genes

Glycogen synthase kinase 3 beta (Gsk3β, (Rn00583429_m1), insulin-like growth factor 1 (Igf1, Rn00710306_m1), insulin-like growth factor 1 receptor (Igf1R, Rn00583837_m1), insulin receptor (InsR, Rn00690703_m1), phosphoinositide-3-kinase regulatory subunit 1 (Pik3R1, Rn00564547_m1), solute carrier family 2, facilitated glucose transporter, member 1 (Glut1, Slc2a1, Rn01417099_m1), solute carrier family 2, facilitated glucose transporter, member 3 (Glut3, Slc2a3, Rn00567331_m1), solute carrier family 2, facilitated glucose transporter, member 4 (Glut4, Slc2a4, Rn00562597_m1), O-GlcNAcase (Oga, Mgea5, Rn00590870_m1), O-GlcNAc-transferase (Ogt, Rn00820779_m1), glutamine fructose-6-phosphate transaminase 1 (Gfpt1, Rn01765495_m1) and glutamine fructose-6-phosphate transaminase 2 (Gfpt2, Rn01456720_m1).

Transcript levels were normalized based on endogenous glyceraldehyde-3-phosphate dehydrogenase (Gapdh, Rn01462662_g1) and tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein (Ywhaz, Rn00755072_m1) expression.

Assay validity was assessed by using appropriate negative controls, in which cDNA was omitted. These negative controls were not read by the software. Acquisition of data (threshold cycle, Ct) was performed by StepOnePlus™ software. A ΔCt (Ct of the considered gene - mean Ct (Gapdh and Ywhaz genes) was calculated.

## Statistical Analysis

All data were analyzed using two-way Analysis Of Variance (ANOVA, with group (CONT vs. PRS) and sex (Males vs. Females) as independent variables, and time (hours) as repeated measure. Fisher's LSD test was used to analyze the differences. Significance was set at a *p*-value of 0.05. In mRNA analysis, in order to subdivide data accounting for significant lower-order effects when F was larger than, trends with 0.05 < *p* < 0.08 were also considered as significant (Snedecor and Cochran, 1967).

## RESULTS

### PRS Induces a Sex-Dependent Disruption of the Circadian Rhythm of Running Wheel Activity

The circadian rhythms of running wheel activity under a regular 12/12 LD cycle were individually analyzed in male and female
CONT unstressed and PRS rats over 10 consecutive days of continuous registration (Figure 1A). PRS induced changes in the circadian rhythm (Group x Sex x Time effect, F(1,47) = 13.241, p < 0.05) with PRS males anticipating the light being switched off by starting the wheel running activity at 19:10 (Figure 1B); in contrast, all the other groups started around 20:00 (which is the precise time of light switch off). PRS females presented reduced number of wheel revolutions with respect to all other groups during the dark phase (Figures 1A,C). Indeed, there was a clear-cut effect of sex on the effect of PRS on the levels of activity defined by the number of wheel revolutions (Figure 1C); first, during the period of lower activity (light phase), both CONT and PRS females were more active than males. During the light phase, PRS increased activity in males, which reached levels of CONT females (Group x Sex effect, F(1,36) = 7.68, p < 0.01). More interestingly, during the period of activity (dark phase), male PRS rats were more active than male CONT rats, in contrast, female PRS rats were less active than CONT females (Group x Sex effect, F(1,36) = 39.48, p < 0.01). In addition, during the dark phase, CONT female rats were less active than CONT male rats. Overall, PRS males presented a phase advance and increased wheel activity, while PRS females were hypoactive and presented a more fragmented profile of activity than did CONT females (Figure 1D).

## PRS Increases Vulnerability to Chronobiological Stress

We measured the number of days needed for the circadian rhythm of locomotor activity to become resynchronized after an abrupt 6-h advance shift in the L/D cycle (Figure 2A). Both PRS and sex influenced the time necessary for resynchronization (ANOVA, Group x Sex effect, F(1,36) = 6.316, p < 0.05). In the CONT group, females required a longer time to resynchronize to the new light shift than the males did (p < 0.05). PRS increased the time to become resynchronized to the new L/D cycle in both sexes with PRS males reaching levels of CONT females. Moreover, the PRS effect occurred to a larger extent in female than in males (p < 0.05 vs. PRS males). Further, PRS females presented a more fragmented activity than PRS males (Figure 2B).

## Effect of PRS and Sex on Circadian Gene Expression in the SCN

The expression of the 42 genes analyzed in the SCN during discrete time points (15:00, 19:00 and 02:00) is reported in Table 1, with genes such as *Gsk3β*, *Igf1*, *mGluR2*, *mGluR3*, *Dkk1*, and *Glast* that were only influenced by PRS or sex (Supplementary Figure S1) and the other genes that followed a circadian rhythm in a sex-dependent manner and were also affected by PRS. Seven genes, did not show changes associated with PRS, sex, or time (*CrhR2*, *mGluR5*, *Htr2c*, *Erbb2*, *Erbb3*, *Erbb4*, *Gfpt1*, *Gfpt2*).

### Signaling- and Glucose Metabolism-Related Genes

The expression of genes involved in signaling such as *mTor*, *Akt1*, *Disc1*, *Dkk1*, *mGluR2*, *mGluR3*, and *Glast* or in glucose

metabolism, such as *Igf1*, *InsR*, *Gsk3β* and *Glut1*, displayed a marked effect of both sex and PRS variables or sex only, whereas circadian oscillations were evident for *mTor*, *Glut1*, *Akt1*, *Disc1*, and *InsR*. As shown in Supplementary Figure S1, most of the sex-biased genes were down-regulated in females, with only *Dkk1* and *Igf1* being upregulated. Interestingly, PRS in males downregulated levels of gene expression compared to CONT males, with levels identical to CONT females (*mTor*, *Gsk3β*, and *Glut1*).

### Clock- and Glucose Metabolism-Related Genes

Interestingly, we observed circadian variations in gene expression that were influenced by both PRS and sex, in the genes that were measured and were involved in clock and glucose metabolism regulation, such as *Clock*, *Bmal*, *Mell1aR*, *Rev-erbA-α*, *Ngf1A*, and *Oga* (Figure 3). Among these genes, sex differences observed in the unstressed group for *Clock* and *Bmal* disappeared in the PRS group. mRNA levels of *Ngf1A* and *Oga* in PRS males showed opposing oscillations when compared to the CONT male group, and this occurred specifically before the onset of the dark phase (19:00) or during it (02:00). Additionally, PRS flipped the sex-specific profile observed in *Rev-erbA-α* and *Mell1aR* when compared to CONT, with PRS males showing the same profile that CONT females, and PRS females displaying an expression profile mimicking CONT males (ANOVA, *Group x Sex x Time effect*, *Clock*, $F_{(2,32)} = 3.229$, *p* = 0.05; *Bmal*, $F_{(2,32)} = 3.901$, *p* < 0.05; *Ngf1A*, $F_{(2,32)} = 3.595$, *p* < 0.05; *Oga*, $F_{(2,32)} = 5.355$, *p* < 0.01; *Mell1aR*, $F_{(2,32)} = 2.811$, *p* = 0.07; *Rev-erbA-α*, $F_{(2,32)} = 7.901$, *p* < 0.01).

### HPA Axis Stress-Related Genes

The effects of PRS and sex were also observed on most of the genes involved in HPA axis activity and regulation that were included in the present study (Figure 4), such as *CrhR1*, *Nr3c2* (*MR*), *AvpR1a*. Here, sex differences observed in unstressed CONT rats in the levels of *CrhR1* and *MR*, which were upregulated in males but not in females, disappeared in the PRS group (ANOVA, *CrhR1* (*Group x Sex effect*, $F_{(1,16)} = 18.012$, *p* < 0.001, *Time effect*, $F_{(2,32)} = 4.434$, *p* < 0.05); *MR* (*Group x Sex effect*, $F_{(1,16)} = 4.520$, *p* < 0.05, *Time effect*, $F_{(2,32)} = 3.140$, *p* = 0.06; *AvpR1a*, *Sex effect*, $F_{(1,16)} = 5.701$, *p* < 0.05, *Time effect*, $F_{(2,32)} = 11.228$, *p* < 0.001). On the other hand, *Nr3c1* (*GR*) and *OxtR* displayed oscillations that were independent of the main variables effect of group or sex (ANOVA, *Time effect only*: *GR*, $F_{(2,32)} = 6.093$, *p* < 0.01; *OxtR*, $F_{(2,32)} = 2.965$, *p* = 0.07) and *CrhR2* expression was unchanged.

## DISCUSSION

We showed a sex- specific profile in the unstressed group that was either exacerbated or flipped by PRS. In unstressed animals, circadian wheel running behavior was greater during the day in females and higher during the night in males. Interestingly, PRS increased and decreased total locomotor activity in males and females, respectively, and it induced a significant phase advance in the rhythm of circadian activity
<figure><img src="image_3.png" /></figure>

exclusively in males. The advance phase observed in PRS males is responsible at least in part for their hyperactivity in the dark cycle. Since in the light period PRS males behaved similarly to unstressed females, PRS induced a dysmasculinization of males' profile. In the dark phase, PRS increased sex differences, with PRS males being hyperactive compared to PRS females, which were hypoactive with respect to all groups. The phase advance observed in PRS males is consistent with the altered circadian functioning of HPA axis activity observed in the PRS model (Koehl et al., 1999). Similar results have also been found using other models of early life stress (i.e., prenatal hypoxia, Joseph et al., 2002). The reduced locomotor activity observed in PRS females could also be related to the changes in the gonadal status that characterize the PRS model. Indeed, the reduced levels in estradiol and impaired estrous cycle observed in PRS females (Reynaert et al., 2016; Van Camp et al., 2018) is associated with the time and intensity of locomotor activity

(Wollnik and Turek, 1988). Consistent with this, we have recently shown that a correction of estradiol levels in PRS females increased their rhythm of locomotor activity (Van Camp et al., 2018). The hyperactivity shown by PRS males may be related to the reduced testosterone levels observed in PRS males (Reynaert et al., 2016). In this line, in male mice, reduced testosterone levels induce a phase advance and increase levels of locomotor activity, while testosterone treatment corrects this profile (Butler et al., 2012).

Sex differences in the unstressed group were also observed in response to a chronobiological stress such as a 6-h advance phase shift, with females being more vulnerable than males since they required three more days to resynchronize to the new L/D cycle. The locomotor response to the chronobiological stress in unstressed females is consistent with the greater activation of the glucocorticoid response to stress that is generally observed in female rodents (Heck and Handa, 2018). Sex differences
<figure><img src="image_4.png" /></figure>

**FIGURE 2** | Locomotor response to a chronobiological stressor in PRS and CONT animals of both sexes. (A) Circadian running wheel activity from CONT-unstressed and PRS rats of both sexes that were subjected to abrupt 6-h phase advance of the L/D cycle. PRS increased the time taken for resynchronization to the new L/D cycle in both sexes, showing increased sensitivity to the chronobiological stressor in PRS animals. Moreover, this occurred to a larger extent in females than in male PRS rats, with 7 days required by PRS females and 3 days by PRS males. (B) Representative activity records of the circadian running wheel activity on day 7. The first L/D cycle and the new L/D cycle, after a 6-h advance in the light period are shown at the top of each panel. Values are expressed as mean ± S.E.M. *n* = 10 rats per group. \*p < 0.05 PRS vs. CONT within the same sex group; \#p < 0.05 females vs. males within the CONT or PRS group.

were observed also in PRS rats, with a greater number of days required by PRS females to re-entrain with respect to PRS males. However, we found that male and female PRS rats resynchronized their activity rhythm to the new L/D cycle more slowly than unstressed rats. Therefore, PRS males behaved like unstressed females. This specific pattern of increased time of resynchronization in the PRS animals is consistent with the

increased stress corticosterone response and rhythm of secretion that has been reported in PRS animals and to a greater extent, in females (Koehl et al., 1997, 1999). Of note, significant increase in the number of days animals required in order to re-entrain to the new LD cycle following a phase shift, have been also observed in other prenatal stress paradigms in both mice and rats, following an 8-h advance/phase shift (chronic mild gestational
<table><caption>TABLE 1 | Effect of perinatal stress and sex on the temporal gene expression in the SCN.</caption><thead><tr><th rowspan="3">Gene name</th><th rowspan="3">Statistical effect</th><th colspan="12">Males</th><th colspan="4">Females</th></tr><tr><th colspan="3">CONT</th><th colspan="3">PRS</th><th colspan="3">CONT</th><th colspan="3">PRS</th><th rowspan="2">02</th><th rowspan="2">Taqman ID (FtX_m1)</th></tr><tr><th>15</th><th>19</th><th>02</th><th>15</th><th>19</th><th>02</th><th>15</th><th>19</th><th>02</th><th>15</th><th>19</th><th>02</th></tr></thead><tbody><tr><td>Clock</td><td>GST</td><td>1.00 ± 0.01</td><td>0.91 ± 0.03</td><td>1.16 ± 0.11</td><td>0.89 ± 0.01</td><td>0.94 ± 0.01</td><td>1.06 ± 0.06</td><td>0.98 ± 0.03</td><td>0.93 ± 0.05</td><td><strong>0.90# ± 0.04</strong></td><td>0.93 ± 0.02</td><td>0.91 ± 0.05</td><td><strong>1.05# ± 0.04</strong></td><td>00573120</td></tr><tr><td>Bmal</td><td>GST</td><td>1.00 ± 0.04</td><td>0.79 ± 0.06</td><td>1.12 ± 0.13</td><td>0.97 ± 0.04</td><td>0.81 ± 0.04</td><td>1.13 ± 0.11</td><td>1.11 ± 0.03</td><td><strong>1.07# ± 0.09</strong></td><td>0.96 ± 0.06</td><td><strong>0.89# ± 0.03</strong></td><td>0.88 ± 0.05</td><td><strong>1.28# ± 0.13</strong></td><td>00577590</td></tr><tr><td>Oga</td><td>GST</td><td>1.00 ± 0.03</td><td>0.95 ± 0.07</td><td>1.08 ± 0.11</td><td>1.05 ± 0.04</td><td>0.90 ± 0.02</td><td>1.00 ± 0.05</td><td>1.13 ± 0.04</td><td>1.01 ± 0.07</td><td><strong>0.89# ± 0.04</strong></td><td><strong>0.94# ± 0.01</strong></td><td>0.96 ± 0.09</td><td><strong>1.09# ± 0.07</strong></td><td>00590870</td></tr><tr><td>Ngrf1A</td><td>GST</td><td>1.00 ± 0.24</td><td>0.56 ± 0.08</td><td>0.52 ± 0.14</td><td><strong>0.66# ± 0.12</strong></td><td>0.73 ± 0.05</td><td><strong>0.91# ± 0.08</strong></td><td>0.76 ± 0.07</td><td>0.61 ± 0.09</td><td>0.56 ± 0.09</td><td>0.86 ± 0.11</td><td>0.80 ± 0.04</td><td><strong>0.62# ± 0.04</strong></td><td>00561138</td></tr><tr><td>Rev-erbA-α</td><td>GST</td><td>1.00 ± 0.06</td><td>1.01 ± 0.09</td><td>0.55 ± 0.07</td><td>1.12 ± 0.07</td><td><strong>0.85# ± 0.05</strong></td><td>0.61 ± 0.06</td><td>1.12 ± 0.02</td><td><strong>0.85# ± 0.04</strong></td><td>0.58 ± 0.02</td><td>1.04 ± 0.03</td><td>0.95 ± 0.08</td><td>0.52 ± 0.02</td><td>00595671</td></tr><tr><td>Mel1AR</td><td>GST</td><td>1.00 ± 0.37</td><td>0.77 ± 0.24</td><td>0.47 ± 0.09</td><td>1.38 ± 0.46</td><td><strong>1.74# ± 0.37</strong></td><td>0.59 ± 0.22</td><td>1.09 ± 0.22</td><td><strong>2.19# ± 0.39</strong></td><td>0.68 ± 0.21</td><td>0.87 ± 0.28</td><td><strong>0.90# ± 0.26</strong></td><td>0.49 ± 0.02</td><td>01488022</td></tr><tr><td>CnrR1</td><td>GS/T</td><td>1.00 ± 0.06</td><td>0.91 ± 0.03</td><td>1.08 ± 0.06</td><td><strong>0.80# ± 0.03</strong></td><td>0.81 ± 0.04</td><td><strong>0.94# ± 0.09</strong></td><td><strong>0.84# ± 0.02</strong></td><td><strong>0.75# ± 0.03</strong></td><td><strong>0.86# ± 0.04</strong></td><td>0.84 ± 0.03</td><td>0.83 ± 0.04</td><td>0.85 ± 0.03</td><td>00578611</td></tr><tr><td>MR</td><td>GS/T</td><td>1.00 ± 0.02</td><td>0.94 ± 0.03</td><td>1.07 ± 0.07</td><td><strong>0.86# ± 0.01</strong></td><td>0.90 ± 0.02</td><td>0.99 ± 0.07</td><td>0.97 ± 0.03</td><td>0.95 ± 0.04</td><td>0.96 ± 0.01</td><td>0.96 ± 0.02</td><td>0.95 ± 0.08</td><td>1.01 ± 0.04</td><td>00561369</td></tr><tr><td>Mtor</td><td>GS/T</td><td>1.00 ± 0.01</td><td>0.91 ± 0.03</td><td>1.00 ± 0.04</td><td><strong>0.87# ± 0.02</strong></td><td>0.92 ± 0.01</td><td>0.99 ± 0.05</td><td>0.93 ± 0.02</td><td>0.86 ± 0.03</td><td><strong>0.92# ± 0.02</strong></td><td>0.95 ± 0.01</td><td>0.91 ± 0.03</td><td>0.92 ± 0.04</td><td>00693900</td></tr><tr><td>Glut1</td><td>GS/T</td><td>1.00 ± 0.06</td><td>1.13 ± 0.09</td><td>1.27 ± 0.17</td><td>0.84 ± 0.03</td><td>0.99 ± 0.04</td><td>1.15 ± 0.09</td><td>0.80 ± 0.02</td><td><strong>0.88# ± 0.03</strong></td><td><strong>0.93# ± 0.02</strong></td><td>0.84 ± 0.03</td><td>1.08 ± 0.17</td><td>0.93 ± 0.03</td><td>01417099</td></tr><tr><td>vGlut</td><td>GT/ST</td><td>1.00 ± 0.28</td><td>0.76 ± 0.11</td><td>2.40 ± 0.38</td><td>1.23 ± 0.50</td><td>0.94 ± 0.19</td><td><strong>3.62# ± 0.56</strong></td><td>1.66 ± 0.45</td><td>0.80 ± 0.13</td><td><strong>0.73# ± 0.13</strong></td><td>1.77 ± 0.14</td><td>1.29 ± 0.17</td><td><strong>2.02# ± 0.23</strong></td><td>01462431</td></tr><tr><td>Mel1BR</td><td>ST</td><td>1.00 ± 0.16</td><td>1.08 ± 0.22</td><td>1.39 ± 0.48</td><td>1.32 ± 0.14</td><td>1.03 ± 0.15</td><td>1.30 ± 0.22</td><td>1.14 ± 0.18</td><td>1.14 ± 0.17</td><td><strong>0.82# ± 0.13</strong></td><td>0.89 ± 0.09</td><td>1.18 ± 0.10</td><td><strong>0.73# ± 0.09</strong></td><td>01447987</td></tr><tr><td>Npas2</td><td>ST</td><td>1.00 ± 0.02</td><td>0.83 ± 0.04</td><td>1.23 ± 0.13</td><td>0.88 ± 0.03</td><td>0.90 ± 0.04</td><td>1.23 ± 0.04</td><td>0.90 ± 0.01</td><td>0.86 ± 0.05</td><td><strong>0.97# ± 0.03</strong></td><td>0.89 ± 0.01</td><td>0.86 ± 0.06</td><td><strong>1.04# ± 0.06</strong></td><td>01438223</td></tr><tr><td>Disc1</td><td>ST</td><td>1.00 ± 0.09</td><td>0.95 ± 0.03</td><td>1.45 ± 0.25</td><td>0.88 ± 0.04</td><td>0.96 ± 0.09</td><td>1.48 ± 0.17</td><td>0.82 ± 0.03</td><td>1.01 ± 0.09</td><td><strong>0.98# ± 0.02</strong></td><td>0.91 ± 0.05</td><td>1.04 ± 0.12</td><td><strong>0.92# ± 0.09</strong></td><td>00598264</td></tr><tr><td>Glut4</td><td>ST</td><td>1.00 ± 0.15</td><td>1.61 ± 0.39</td><td>1.92 ± 0.34</td><td>1.06 ± 0.10</td><td>1.16 ± 0.09</td><td>2.15 ± 0.35</td><td>1.01 ± 0.08</td><td>1.34 ± 0.04</td><td><strong>1.27# ± 0.13</strong></td><td>1.00 ± 0.12</td><td>1.31 ± 0.23</td><td><strong>1.08# ± 0.11</strong></td><td>00562597</td></tr><tr><td>Igf1R</td><td>ST</td><td>1.00 ± 0.03</td><td>1.04 ± 0.05</td><td>1.26 ± 0.13</td><td>0.94 ± 0.02</td><td>0.95 ± 0.02</td><td>1.38 ± 0.17</td><td>1.01 ± 0.01</td><td>1.11 ± 0.03</td><td><strong>1.11# ± 0.03</strong></td><td>1.03 ± 0.02</td><td>1.23 ± 0.17</td><td><strong>1.17# ± 0.14</strong></td><td>00583837</td></tr><tr><td>Pik3R1</td><td>ST</td><td>1.00 ± 0.07</td><td>1.17 ± 0.07</td><td>1.32 ± 0.09</td><td>1.03 ± 0.04</td><td>1.11 ± 0.04</td><td>1.27 ± 0.07</td><td><strong>1.24# ± 0.05</strong></td><td>1.14 ± 0.05</td><td><strong>1.11# ± 0.03</strong></td><td>1.11 ± 0.05</td><td>1.16 ± 0.10</td><td>1.33 ± 0.03</td><td>00564547</td></tr><tr><td>Cry1</td><td>S/T</td><td>1.00 ± 0.04</td><td>1.23 ± 0.07</td><td>1.47 ± 0.22</td><td>1.00 ± 0.04</td><td>1.14 ± 0.07</td><td>1.39 ± 0.10</td><td>0.90 ± 0.03</td><td>1.05 ± 0.06</td><td>1.32 ± 0.05</td><td>0.89 ± 0.01</td><td>1.15 ± 0.16</td><td>0.79 ± 0.06</td><td>01503063</td></tr><tr><td>Pgr1</td><td>S/T</td><td>1.00 ± 0.08</td><td>1.13 ± 0.12</td><td>0.89 ± 0.08</td><td>0.85 ± 0.11</td><td>1.04 ± 0.05</td><td>1.03 ± 0.11</td><td>0.84 ± 0.06</td><td>0.92 ± 0.05</td><td>0.76 ± 0.01</td><td>0.86 ± 0.02</td><td>1.18 ± 0.11</td><td>0.79 ± 0.06</td><td>01496757</td></tr><tr><td>Avpr1a</td><td>S/T</td><td>1.00 ± 0.09</td><td>1.13 ± 0.10</td><td>1.32 ± 0.22</td><td>1.04 ± 0.05</td><td>1.05 ± 0.04</td><td>1.40 ± 0.10</td><td>1.14 ± 0.04</td><td>1.10 ± 0.02</td><td>1.51 ± 0.09</td><td>1.18 ± 0.07</td><td>1.18 ± 0.11</td><td>1.43 ± 0.08</td><td>00583910</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 1 | Continued</caption><thead><tr><th rowspan="3">Gene name</th><th rowspan="3">Statistical effect</th><th colspan="10">Males</th><th colspan="4">Females</th></tr><tr><th colspan="3">CONT</th><th colspan="3">PRS</th><th colspan="3">CONT</th><th colspan="3">PRS</th><th rowspan="2">Taqman ID (RnX_m1)</th></tr><tr><th>15</th><th>19</th><th>02</th><th>15</th><th>19</th><th>02</th><th>15</th><th>19</th><th>02</th><th>15</th><th>19</th><th>02</th></tr></thead><tbody><tr><td>Akt1</td><td>S/T</td><td>1.00 ± 0.03</td><td>0.95 ± 0.01</td><td>1.07 ± 0.04</td><td>0.91 ± 0.03</td><td>0.92 ± 0.03</td><td>1.13 ± 0.13</td><td><strong>0.86# ± 0.03</strong></td><td>0.89 ± 0.03</td><td><strong>0.93# ± 0.02</strong></td><td>0.96 ± 0.01</td><td>0.91 ± 0.05</td><td><strong>0.93# ± 0.02</strong></td><td>00583646</td></tr><tr><td>nsR</td><td>S/T</td><td>1.00 ± 0.02</td><td>0.90 ± 0.04</td><td>1.08 ± 0.12</td><td>0.91 ± 0.02</td><td>0.88 ± 0.02</td><td>0.97 ± 0.07</td><td>0.92 ± 0.03</td><td>0.85 ± 0.02</td><td>0.88 ± 0.02</td><td>0.88 ± 0.02</td><td>0.90 ± 0.05</td><td>0.95 ± 0.03</td><td>00690703</td></tr><tr><td>Gsk3b</td><td>GS</td><td>1.00 ± 0.02</td><td>1.01 ± 0.03</td><td>1.01 ± 0.02</td><td>0.95 ± 0.03</td><td>0.95 ± 0.04</td><td>0.92 ± 0.02</td><td>0.96 ± 0.01</td><td>0.88 ± 0.04</td><td>0.94 ± 0.03</td><td>0.93 ± 0.02</td><td>0.93 ± 0.06</td><td>0.91 ± 0.04</td><td>00583429</td></tr><tr><td>Igfr1</td><td>S</td><td>1.00 ± 0.13</td><td>1.16 ± 0.16</td><td>1.08 ± 0.11</td><td>1.13 ± 0.08</td><td>1.18 ± 0.05</td><td>1.28 ± 0.12</td><td>1.80 ± 0.46</td><td>1.74 ± 0.19</td><td>1.59 ± 0.12</td><td>1.68 ± 0.25</td><td>1.75 ± 0.14</td><td>1.91 ± 0.25</td><td>00710306</td></tr><tr><td>mGlu2</td><td>S</td><td>1.00 ± 0.15</td><td>0.90 ± 0.12</td><td>1.04 ± 0.13</td><td>1.02 ± 0.07</td><td>0.70 ± 0.10</td><td>1.19 ± 0.39</td><td>0.92 ± 0.09</td><td>0.60 ± 0.11</td><td>0.52 ± 0.08</td><td>0.72 ± 0.10</td><td>0.79 ± 0.09</td><td>0.68 ± 0.07</td><td>01447672</td></tr><tr><td>mGlu3</td><td>S</td><td>1.00 ± 0.07</td><td>0.98 ± 0.08</td><td>1.11 ± 0.10</td><td>0.82 ± 0.04</td><td>0.81 ± 0.06</td><td>1.20 ± 0.26</td><td>0.89 ± 0.05</td><td>0.78 ± 0.10</td><td>0.82 ± 0.04</td><td>0.82 ± 0.04</td><td>0.95 ± 0.12</td><td>0.89 ± 0.06</td><td>01755352</td></tr><tr><td>Dkk1</td><td>S</td><td>1.00 ± 0.12</td><td>0.89 ± 0.09</td><td>1.46 ± 0.20</td><td>1.03 ± 0.11</td><td>0.87 ± 0.09</td><td>1.63 ± 0.38</td><td>1.04 ± 0.15</td><td>1.74 ± 0.39</td><td>1.17 ± 0.22</td><td>1.41 ± 0.19</td><td>1.22 ± 0.27</td><td>1.56 ± 0.09</td><td>01501537</td></tr><tr><td>Glast</td><td>S</td><td>1.00 ± 0.03</td><td>1.06 ± 0.03</td><td>1.05 ± 0.08</td><td>1.02 ± 0.03</td><td>0.99 ± 0.08</td><td>1.22 ± 0.14</td><td>0.98 ± 0.02</td><td>0.93 ± 0.03</td><td>0.97 ± 0.03</td><td>0.98 ± 0.02</td><td>0.99 ± 0.03</td><td>0.94 ± 0.02</td><td>00570130</td></tr><tr><td>Per2</td><td>T</td><td>1.00 ± 0.05</td><td>1.43 ± 0.13</td><td>1.12 ± 0.11</td><td>0.98 ± 0.08</td><td>1.27 ± 0.09</td><td>1.19 ± 0.11</td><td>1.03 ± 0.01</td><td>1.22 ± 0.05</td><td>1.00 ± 0.04</td><td>1.02 ± 0.02</td><td>1.38 ± 0.20</td><td>0.91 ± 0.05</td><td>01427704</td></tr><tr><td>Per3</td><td>T</td><td>1.00 ± 0.02</td><td>1.20 ± 0.07</td><td>1.09 ± 0.11</td><td>0.93 ± 0.05</td><td>1.11 ± 0.02</td><td>1.04 ± 0.07</td><td>0.97 ± 0.02</td><td>1.10 ± 0.05</td><td>1.07 ± 0.04</td><td>0.94 ± 0.01</td><td>1.25 ± 0.18</td><td>0.96 ± 0.03</td><td>00709499</td></tr><tr><td>GR</td><td>T</td><td>1.00 ± 0.03</td><td>1.04 ± 0.07</td><td>1.29 ± 0.12</td><td>0.95 ± 0.04</td><td>1.04 ± 0.07</td><td>1.16 ± 0.07</td><td>1.01 ± 0.05</td><td>1.08 ± 0.06</td><td>1.09 ± 0.06</td><td>0.95 ± 0.02</td><td>1.07 ± 0.10</td><td>1.04 ± 0.07</td><td>00565562</td></tr><tr><td>Oxfr</td><td>T</td><td>1.00 ± 0.16</td><td>0.97 ± 0.16</td><td>0.76 ± 0.12</td><td>0.89 ± 0.12</td><td>1.23 ± 0.19</td><td>0.96 ± 0.20</td><td>0.87 ± 0.08</td><td>1.05 ± 0.11</td><td>0.75 ± 0.03</td><td>0.87 ± 0.07</td><td>1.22 ± 0.26</td><td>1.01 ± 0.11</td><td>00563503</td></tr><tr><td>Glut3</td><td>T</td><td>1.00 ± 0.03</td><td>0.88 ± 0.02</td><td>0.99 ± 0.04</td><td>0.87 ± 0.02</td><td>0.91 ± 0.03</td><td>0.99 ± 0.04</td><td>0.95 ± 0.02</td><td>0.87 ± 0.04</td><td>0.95 ± 0.03</td><td>0.93 ± 0.00</td><td>0.89 ± 0.03</td><td>0.91 ± 0.05</td><td>00567331</td></tr><tr><td>Ogt</td><td>T</td><td>1.00 ± 0.06</td><td>1.02 ± 0.04</td><td>1.24 ± 0.14</td><td>1.09 ± 0.06</td><td>1.26 ± 0.10</td><td>1.31 ± 0.18</td><td>1.28 ± 0.05</td><td>1.07 ± 0.09</td><td>1.17 ± 0.08</td><td>1.00 ± 0.03</td><td>1.18 ± 0.14</td><td>1.26 ± 0.11</td><td>00820779</td></tr><tr><td>mGlu5</td><td>–</td><td>1.00 ± 0.08</td><td>0.82 ± 0.05</td><td>0.88 ± 0.06</td><td>0.90 ± 0.06</td><td>0.91 ± 0.04</td><td>1.10 ± 0.16</td><td>0.94 ± 0.04</td><td>0.83 ± 0.07</td><td>0.88 ± 0.04</td><td>0.95 ± 0.04</td><td>0.82 ± 0.06</td><td>0.84 ± 0.03</td><td>00566628</td></tr><tr><td>Htr2c</td><td>–</td><td>1.00 ± 0.06</td><td>1.05 ± 0.06</td><td>1.10 ± 0.17</td><td>1.07 ± 0.09</td><td>1.03 ± 0.11</td><td>1.00 ± 0.22</td><td>0.93 ± 0.06</td><td>1.17 ± 0.17</td><td>0.92 ± 0.05</td><td>0.85 ± 0.01</td><td>1.23 ± 0.08</td><td>1.00 ± 0.07</td><td>00562748</td></tr><tr><td>CnrR2</td><td>–</td><td>1.00 ± 0.12</td><td>0.99 ± 0.08</td><td>0.89 ± 0.07</td><td>0.92 ± 0.11</td><td>1.05 ± 0.11</td><td>0.86 ± 0.02</td><td>1.00 ± 0.11</td><td>1.00 ± 0.12</td><td>1.00 ± 0.05</td><td>1.07 ± 0.07</td><td>1.00 ± 0.06</td><td>0.94 ± 0.03</td><td>00575617</td></tr><tr><td>ErbB2</td><td>–</td><td>1.00 ± 0.09</td><td>1.07 ± 0.07</td><td>1.00 ± 0.12</td><td>1.02 ± 0.03</td><td>1.02 ± 0.04</td><td>1.12 ± 0.11</td><td>0.97 ± 0.06</td><td>1.37 ± 0.15</td><td>1.07 ± 0.07</td><td>0.95 ± 0.07</td><td>1.15 ± 0.11</td><td>1.04 ± 0.05</td><td>00566661</td></tr><tr><td>ErbB3</td><td>–</td><td>1.00 ± 0.07</td><td>1.01 ± 0.09</td><td>1.27 ± 0.13</td><td>0.87 ± 0.05</td><td>0.93 ± 0.07</td><td>0.96 ± 0.11</td><td>1.01 ± 0.06</td><td>0.97 ± 0.07</td><td>1.09 ± 0.06</td><td>0.90 ± 0.04</td><td>1.07 ± 0.10</td><td>1.01 ± 0.09</td><td>00568107</td></tr><tr><td>ErbB4</td><td>–</td><td>1.00 ± 0.03</td><td>0.97 ± 0.03</td><td>1.12 ± 0.12</td><td>0.93 ± 0.03</td><td>0.91 ± 0.03</td><td>1.07 ± 0.10</td><td>1.08 ± 0.06</td><td>0.93 ± 0.02</td><td>0.92 ± 0.03</td><td>0.96 ± 0.02</td><td>0.99 ± 0.07</td><td>0.98 ± 0.06</td><td>00572447</td></tr><tr><td>Gpt1</td><td>–</td><td>1.00 ± 0.03</td><td>0.99 ± 0.04</td><td>1.11 ± 0.08</td><td>0.96 ± 0.04</td><td>0.97 ± 0.03</td><td>1.06 ± 0.04</td><td>1.01 ± 0.04</td><td>0.97 ± 0.02</td><td>0.97 ± 0.02</td><td>0.99 ± 0.02</td><td>1.00 ± 0.06</td><td>1.00 ± 0.04</td><td>01765496</td></tr><tr><td>Gpt2</td><td>–</td><td>1.00 ± 0.06</td><td>0.97 ± 0.04</td><td>1.10 ± 0.09</td><td>0.93 ± 0.03</td><td>0.90 ± 0.02</td><td>1.12 ± 0.11</td><td>0.90 ± 0.04</td><td>0.93 ± 0.03</td><td>0.95 ± 0.02</td><td>0.92 ± 0.04</td><td>1.00 ± 0.10</td><td>0.95 ± 0.06</td><td>01456720</td></tr></tbody></table>

Expression of selected genes in the SCN of unstressed or PRS rats of both sexes was assayed by Custom Taqman qRT-PCR. For each sample, expression levels were normalized to the average levels of unstressed male rats obtained at 15:00 and expressed as Relative Quantitation (RQ). Genes are listed according to the main statistical effect: GST, Group × Stress × Time effect; GS, Group × Sex effect; GT, group × time; ST, sex × time effect; G, S, T for Group, Sex, or Time effect only. Data are means ± S.E.M of five determinations for each group. *p < 0.05 for PRS vs. CONT of the same-sex group; #p < 0.05 for females vs. males within the PRS or CONT group. Bolded values are those with greater significance.
Expression of clock, stress and glucose-related genes in the SCN during light/dark phase

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

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

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

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

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

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

**FIGURE 3** | Circadian expression of clock and glucose metabolism-related genes in the SCN of PRS and CONT animals of both sexes. Expression of selected genes in the SCN of CONT or PRS rats of both sexes was assayed by Custom Taqman qRT-PCR. Expression levels of each sample were normalized to the average levels in CONT male rats obtained at 15:00 and expressed as Relative Quantitation (RQ). Effect of PRS on gene expression was observed in both sexes for all the genes considered. PRS induced a dysmasculinization in the profile of expression for *Clock* and *Mel1aR* in males, which showed the same expression profile as CONT females. In *Bmal*, *Rev-erb-α* and *Mel1aR*, PRS females behaved as CONT males. Data are depicted as mean ± S.E.M of five determinations for each group.
Expression of HPA axis stress-related genes in the SCN during light/dark phase

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

**FIGURE 4** | Circadian expression of HPA axis stress-related genes in the SCN of PRS and CONT animals of both sexes. Expression of selected genes in the SCN of CONT or PRS rats of both sexes was assayed by Custom Taqman qRT-PCR. Expression levels of each sample were normalized to the average levels in CONT male rats obtained at 15:00 and expressed as Relative Quantitation (RQ). All genes presented oscillation with the time exception made for *CrhR2*. PRS induced a dysmasculinization in the profile of expression for *CrhR1* in males, which showed the same expression profile as CONT females. All values are represented as mean ± S.E.M of five animals per group.
stress, Kiryanova et al., 2017) or an abrupt 6-h delay/phase shift (prenatal hypoxia Joseph et al., 2002).

Remarkably, the pattern of locomotor activity in PRS rats under a 12:12 h L/D cycle as well as in response to the chronobiological stress was erratic and fragmented, and this was particularly evident in female PRS rats. This specific fragmentation profile has already been shown for another circadian rhythm – the sleep-wake cycle – with PRS males showing discontinuous and shorter sleep when compared to unstressed animals (Dugovic et al., 1999; Mairesse et al., 2013). Moreover, fragmented circadian activity rhythms reflect age-related alterations in the biological clock, because circadian rhythms and sleep patterns change with age (Turek et al., 1995; Van Someren, 2000). Along these lines, the characteristic pattern of fragmentation in PRS rats would lean toward an accelerated aging of the circadian clock induced by PRS. This would broaden the concept of PRS-mediated accelerated aging processes to predictive adaptation, which has already been put forward at the level of reactive adaptation. Indeed, adult PRS males display prolonged secretion of plasma corticosterone in response to stress compared to unstressed aged animals (Vallée et al., 1999), and PRS females display reduced levels of estradiol already in adulthood (Reynaert et al., 2016) similar to unstressed females during aging (Van Camp et al., 2018).

The main actor in the regulation of central and peripheral circadian rhythms is the SCN. Interestingly, we observed patterns of oscillations for the genes analyzed in the SCN that were changed by PRS and in a sex-dependent manner. Indeed, the sex-specific profile of gene expression that was reported in almost all of the genes studied in the SCN in unstressed animals, was flipped by PRS, with PRS males displaying a pattern similar to that of unstressed females. As observed in the behavioral phenotype, this pattern of dysmasculinization also affected gene expression changes.

We observed that females had downregulated expression of the genes involved in signaling regulation, with only *Dkk1* and *Igf1* being upregulated. In particular, for glucose transporters, marked sex differences exist in the gene expression profiles in rodents' tissues (Nagai et al., 2014). Here PRS downregulated *mTor*, *Glut1*, and *Gsk3b* compared to unstressed males. Moreover, PRS males displayed a profile identical to unstressed females, again indicating a pattern of dysmasculinization. Downregulation of *Glut1* in the SCN of PRS males suggests an impairment of glucose metabolism in the brain. Indeed, we previously reported downregulation of *GLUT1* in the placenta of male PRS fetuses (Mairesse et al., 2007).

Among the genes involved in clock and glucose regulation, we showed upregulation of *Clock* in unstressed males compared to unstressed females and upregulation of *Bmal*, *Oga*, and to a greater extent for *Mel1aR* in unstressed females. Levels of *Clock* in PRS males, were downregulated at night (02:00) with respect to unstressed males, thus displaying an unstressed female-like profile. Recent studies in rats indicate that *Clock* is expressed across brain tissues in a sex-specific pattern (Chun et al., 2015). PRS in males also downregulated *Ngf1A* during

the day and higher levels during the night when compared to both unstressed animals and PRS females. Interestingly, *Bmal* and *Rev-erb-α* in PRS females displayed a masculinized profile, in that it mimicked the one observed in unstressed males. Remarkably, we observed that PRS completely flipped the sex-dependent expression of *Mel1aR* transcript levels, with PRS males showing the same profile as unstressed females with higher levels of *Mel1aR* expression, and PRS females displaying an expression profile similar to that of unstressed males, with reduced levels of the transcript. Within the stress response genes, we observed minor effects with upregulation of *CrhR1* mRNA levels in the SCN in males compared to females in the unstressed group. Overall, among all the genes studied, the unique profile in *Mel1aR* expression indicates that this gene is a key regulator in the programming of the circadian system induced by PRS.

## CONCLUSION

Our results provide a first characterization of sex differences in the behavioral and molecular chronobiological profile in unstressed rats. PRS affected males to greater extent than females on the set-up of the circadian system, and lead to the dysmasculinization of PRS males related to locomotor activity phenotype, resynchronization to the new L/D cycle, and in gene expression in the SCN by suggesting a critical role for *Mel1aR* in both sexes. To our knowledge, this is the first demonstration of sex-dependent studies on predictive adaptation programmed by early life stress. Moreover, our results provide additional information, this time at the gene expression level, on the specific pattern of dysmasculinization, which has already been reported at the sexual, physiological, and behavioral level in animal studies of prenatal stress in the past decades (Ward, 1972, 1984; Becker and Kowall, 1977; Dahlöf et al., 1977; Dörner et al., 1983; Ward and Weisz, 1984; Koehl et al., 1997, 1999; Weinstock, 2001; Reznikov and Tarasenko, 2007). The disrupting effect induced by PRS has already been shown in reactive adaptation with its effects on stress response and glutamatergic transmission in the hippocampus. The results obtained in the present study extend this disruption to the main biological clock - and therefore to predictive adaptation. Overall, perturbation of the early life environment has profound consequences on the ability of the organism to react to changes (stress response) or anticipate those changes (circadian rhythms). Therefore, it is necessary to consider sex differences in translational research, in order to gain knowledge of this under-explored field. Animal models of stress- and sleep-disorders as well as sex differences satisfy translational criteria of validity and thus provide a great opportunity of investigation.

## ETHICS STATEMENT

This study was carried out in accordance with the recommendations of the committee on animal experimentation CEEA75, Lille France. All experiments followed the rules of the European Communities Council Directive 2010/63/EU.
## AUTHOR CONTRIBUTIONS

SM-F, JeM, GVC, JoM, FN, and SM designed the study and wrote the manuscript. SM-F, JeM, GVC, M-LR, EG, JoM, and HB conducted the experiments and analyzed the data.

## FUNDING

This study was supported by the University of Lille, Sapienza University of Rome (Framework Agreement signed on February 15th, 2007), and the CNRS within the framework of the Prenatal Stress and Neurodegenerative Diseases International Associated Laboratory (LIA-PSND), co-directed by Professors SM (France) and FN (Italy). Behavioral assessments were carried out using the platform PHExMAR of the University of Lille. We are indebted to the Research Federation FRABio (University of Lille, CNRS, FR 3688, FRABio, Biochimie Structurale et Fonctionnelle des Assemblages Biomoléculaires) for providing the scientific and technical environment conducive to completing this body of work. Transcriptomic analysis was performed using the platform QPCR Imagif-CNRS, at Gif-sur-Yvette. English editing was provided by Editing Elsevier services.

## REFERENCES

Albers, H., Gerall, A., and Axelson, J. (1981). Effect of reproductive state on circadian periodicity in the rat. *Physiol. Behav.* 26, 21–25. doi: 10.1016/0031-9384(81)90073-1

Bailey, M., and Silver, R. (2014). Sex differences in circadian timing systems: implications for disease. *Front. Neuroendocrinol.* 35, 111–139. doi: 10.1016/j. yfrne.2013.11.003

Becker, G., and Kowall, M. (1977). Crucial role of the postnatal maternal environment in the expression of prenatal stress effects in the male rats. *J. Comp. Physiol. Psychol.* 91, 1432–1446. doi: 10.1037/h0077401

Beery, A. K., and Zucker, I. (2011). Sex bias in neuroscience and biomedical research. *Neurosci. Biobehav. Rev.* 35, 565–572. doi: 10.1016/j.neubiorev.2010.07.002

Brunton, P. J., Donadio, M. V. F., and Russell, J. A. (2011). Sex differences in prenatally programmed anxiety behaviour in rats: differential corticotropin-releasing hormone receptor mRNA expression in the amygdaloid complex. *Stress* 14, 634–643. doi: 10.3109/10253890.2011.604750

Buijs, R. M., and Kalsbeek, A. (2001). Hypothalamic integration of central and peripheral clocks. *Nat. Rev. Neurosci.* 2, 521–526. doi: 10.1038/35081582

Bunney, J. N., and Potkin, S. G. (2008). Circadian abnormalities, molecular clock genes and chronobiological treatments in depression. *Br. Med. Bull.* 86, 23–32. doi: 10.1093/bmb/ldn019

Butler, M. P., Karatsoreos, I. N., LeSauter, J., and Silver, R. (2012). Dose-dependent effects of androgens on the circadian timing system and its response to light. *Endocrinology* 153, 2344–2352. doi: 10.1210/en.2011-1842

Chapman, R. H., and Stern, J. M. (1979). Failure of severe maternal stress or ACTH during pregnancy to affect emotionality of male rat offspring: implications of litter effects for prenatal studies. *Dev. Psychobiol.* 12, 255–267. doi: 10.1002/dev. 420120309

Chun, L. E., Woodruff, E. R., Morton, S., Hinds, L. R., and Spencer, R. L. (2015). Variations in phase and amplitude of rhythmic clock gene expression across prefrontal cortex, hippocampus, amygdala, and hypothalamic paraventricular and suprachiasmatic nuclei of male and female rats. *J. Biol. Rhythm.* 30, 417–436. doi: 10.1177/0748730415598608

Cirelli, C. (2009). The genetic and molecular regulation of sleep: from fruit flies to humans. *Nat. Rev. Neurosci.* 10, 549–560. doi: 10.1038/nrn2683

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnmol. 2019.00089/full#supplementary-material

**FIGURE S1** | Effect of PRS and sex differences on the expression levels of signaling- and glucose metabolism-related genes in the SCN. Expression of selected genes in the SCN of CONT or PRS rats of both sexes was assayed by Custom Taqman qRT-PCR. Expression levels of each sample were normalized to the average levels in CONT male rats obtained at 15:00 and expressed as Relative Quantitation (RQ). All genes displayed a marked effect of both sex and PRS variables or sex only. Here, gene expression is pulled over time variable. Most of the sex differences were in the direction of reduced expression levels for females, with only *Dkk1* and *Igf1* being upregulated in females. PRS males displayed reduced levels of expression for *mTor*, *Glut1*, and *Gsk3b* and therefore, displayed a profile identical to CONT females. Data are depicted as mean ± S.E.M of five determinations for each group (ANOVA, *Group x Sex effect*, *mTor F*(1,16) = 5.750, *p* < 0.05; *Glut1 F*(1,16) = 5.051, *p* < 0.05; *Gsk3b F*(1,16) = 3.288, *p* = 0.07; *Sex effect only*, *mGlu2*, *F*(1,16) = 9.268, *p* < 0.01; *mGlu3 F*(1,16) = 5.894, *p* < 0.05; *Akt1*, *F*(1,16) = 9.747, *p* < 0.01; *Disc*, *F*(1,16) = 5.991 *p* < 0.05; *Dkk1*, *F*(1,16) = 4.855 *p* < 0.05; *Glast*, *F*(1,16) = 7.617, *p* < 0.01; *InsR*, *F*(1,16) = 3.893, *p* = 0.06; *Time effect*, *mTor F*(2,32) = 4.180, *p* < 0.05; *Glut1*, *F*(2,32) = 6.046, *p* < 0.01; *Akt1*, *F*(2,32) = 5.151, *p* < 0.05; *Disc1*, *F*(2,32) = 8.906, *p* < 0.001; *InsR*, *F*(2,32) = 4.050, *p* < 0.05). \* *p* < 0.05 PRS vs. CONT within the same sex group; \*\* *p* < 0.05 and \*\*\* *p* < 0.01 females vs. males within the CONT or PRS group.

Dahlöf, L. G., Härd, E., and Larsson, K. (1977). Influence of maternal stress on offspring sexual behavior. *Anim. Behav.* 25, 958–968. doi: 10.1016/0003-3472(77)90047-1

Darnaudéry, M., and Maccari, S. (2008). Epigenetic programming of the stress response in male and female rats by prenatal restraint stress. *Brain Res. Rev.* 57, 571–585. doi: 10.1016/j.brainresrev.2007.11.004

de Kloet, E. R. (1999). Neurobiology of adaptation disorders. *Acta Neuropsychiatr.* 11, 11–14. doi: 10.1017/S0924270800036292

Dibner, C., Schibler, U., and Albrecht, U. (2010). The mammalian circadian timing system: organization and coordination of central and peripheral clocks. *Annu. Rev. Physiol.* 72, 517–549. doi: 10.1146/annurev-physiol-021909-135821

Dörner, G., Götz, F., and Döcke, W. D. (1983). Prevention of demasculinization and feminization of the brain in prenatally stressed male rats by perinatal androgen treatment. *Exp. Clin. Endocrinol.* 81, 88–90. doi: 10.1055/s-0029-1210211

Dugovic, C., Maccari, S., Weibel, L., Turek, F. W., and Van Reeth, O. (1999). High corticosterone levels in prenatally stressed rats predict persistent paradoxical sleep alterations. *J. Neurosci.* 19, 8656–8664. doi: 10.1523/JNEUROSCI.19-19-08656.1999

Gatta, E., Mairesse, J., Deruyter, L., Marrocco, J., Van Camp, G., Bouwalerh, H., et al. (2018). Reduced maternal behavior caused by gestational stress is predictive of life span changes in risk-taking behavior and gene expression due to altering of the stress/anti-stress balance. *Neurotoxicology* 66, 138–149. doi: 10.1016/j.neuro.2018.04.005

Gerstner, J. R., and Yin, J. C. P. (2010). Circadian rhythms and memory formation. *Nat. Rev. Neurosci.* 11, 577–588. doi: 10.1038/nrn2881

Heck, A. L., and Handa, R. J. (2018). Sex differences in the hypothalamic-pituitary-adrenal axis response to stress: an important role for gonadal hormones. *Neuropsychopharmacology* 44, 45–58. doi: 10.1038/s41386-018-0167-9

Hughes, M. E., DiTacchio, L., Hayes, K. R., Vollmers, C., Pulivarthy, S., Baggs, J. E., et al. (2009). Harmonics of circadian gene transcription in mammals. *PLoS Genet.* 5:e1000442. doi: 10.1371/journal.pgen.1000442

Joseph, V., Mamet, J., Lee, F., Dalmaz, Y., and Van Reeth, O. (2002). Prenatal hypoxia impairs circadian synchronisation and response of the biological clock to light in adult rats. *J. Physiol.* 543, 387–395. doi: 10.1113/jphysiol.2002.022236

Kiryanova, V., Smith, V. M., Dyck, R. H., and Antle, M. C. (2017). Circadian behavior of adult mice exposed to stress and fluoxetine during development. *Psychopharmacology* 234, 793–804. doi: 10.1007/s00213-016-4515-3
Klein, D. C., Moore, R. Y., and Reppert, S. M. (1991). *Suprachiasmatic Nucleus: The Mind's Clock*. New York, NY: Oxford University Press.

Koehl, M., Barbazanges, A., Le Moal, M., and Maccari, S. (1997). Prenatal stress induces a phase advance of circadian corticosterone rhythm in adult rats which is prevented by postnatal stress. *Brain Res.* 759, 317-320. doi: 10.1016/S0006-8993(97)00394-6

Koehl, M., Darnaudéry, M., Dulluc, J., Van Reeth, O., Le Moal, M., and Maccari, S. (1999). Prenatal stress alters circadian activity of hypothalamo-pituitary-adrenal axis and hippocampal corticosteroid receptors in adult rats of both gender. *J. Neurobiol.* 40, 302-315. doi: 10.1002/(SICI)1097-4695(19990905)40: 3<302::AID-NEU3>3.0.CO;2-7

Liyanarachchi, K., Ross, R., and Debono, M. (2017). Human studies on hypothalamo-pituitary-adrenal (HPA) axis. *Best Pract. Res. Clin. Endocrinol. Metab.* 31, 459-473. doi: 10.1016/j.beem.2017.10.011

Maccari, S., Darnaudery, M., Morley-Fletcher, S., Zuena, A. R., Cinque, C., and Van Reeth, O. (2003). Prenatal stress and long-term consequences: implications of glucocorticoid hormones. *Neurosci. Biobehav. Rev.* 27, 119-127. doi: 10.1016/S0149-7634(03)00014-9

Maccari, S., and Morley-Fletcher, S. (2007). Effects of prenatal restraint stress on the hypothalamus-pituitary-adrenal axis and related behavioural and neurobiological alterations. *Psychoneuroendocrinology* 32(Suppl. 1), S10-S15. doi: 10.1016/j.psyneuen.2007.06.005

Maccari, S., Piazza, P. V., Kabbaj, M., Barbazanges, A., Simon, H., and Le Moal, M. (1995). Adoption reverses the long-term impairment in glucocorticoid feedback induced by prenatal stress. *J. Neurosci.* 15, 110-116. doi: 10.1523/JNEUROSCI. 15-01-00110.1995

Maccari, S., Polese, D., Reynaert, M. L., Amici, T., Morley-Fletcher, S., and Fagioli, F. (2017). Early-life experiences and the development of adult diseases with a focus on mental illness: the human birth theory. *Neuroscience* 342, 232-251. doi: 10.1016/j.neuroscience.2016.05.042

Mairesse, J., Gatta, E., Reynaert, M. L., Marrocco, J., Morley-Fletcher, S., Soichot, M., et al. (2015). Activation of presynaptic oxytocin receptors enhances glutamate release in the ventral hippocampus of prenatally restraint stressed rats. *Psychoneuroendocrinology* 62, 36-46. doi: 10.1016/j.psyneuen.2015.07.005

Mairesse, J., Lesage, J., Breton, C., Bréant, B., Hahn, T., Darnaudéry, M., et al. (2007). Maternal stress alters endocrine function of the feto-placental unit in rats. *Am. J. Physiol. Endocrinol. Metab.* 292, E1526-E1533. doi: 10.1152/ ajpendo.00574.2006

Mairesse, J., Silletti, V., Laloux, C., Zuena, A. R., Giovine, A., Consolazione, M., et al. (2013). Chronic agomelatine treatment corrects the abnormalities in the circadian rhythm of motor activity and sleep/wake cycle induced by prenatal restraint stress in adult rats. *Int. J. Neuropsychopharmacol.* 16, 323-338. doi: 10.1017/S1461145711001970

McCarthy, M. M., Woolley, C. S., and Arnold, A. P. (2017). Incorporating sex as a biological variable in neuroscience: what do we gain? *Nat. Rev. Neurosci.* 12, 707-708. doi: 10.1038/nrn.2017.137

McEwen, B. S. (1998). Stress, adaptation, and disease. *allostasis and allostatic load*. *Ann. N.Y. Acad. Sci.* 1, 33-44. doi: 10.1111/j.1749-6632.1998.tb09546.x

McEwen, B. S. (2016). In pursuit of resilience: stress, epigenetics, and brain plasticity. *Ann. N.Y. Acad. Sci.* 1373, 56-64. doi: 10.1111/nyas.13020

McEwen, B. S., and Stellar, E. (1993). Stress and the individual mechanisms leading to disease. *Arch. Intern. Med.* 153, 2093-2101. doi: 10.1001/archinte. 1993.00410180039004

Medic, G., Wille, M., and Hemels, M. E. (2017). Short- and long-term health consequences of sleep disruption. *Nat. Sci. Sleep* 9, 151-161. doi: 10.2147/NSS. S134864

Mong, J. A., Baker, F. C., Mahoney, M. M., Paul, K. N., Schwartz, M. D., Semba, K., et al. (2011). Sleep, rhythms, and the endocrine brain: influence of sex and gonadal hormones. *J. Neurosci.* 31, 16107-16116. doi: 10.1523/JNEUROSCI. 4175-11.2011

Moore, R. Y., and Eichler, V. B. (1972). Loss of a circadian adrenal corticosterone rhythm following suprachiasmatic lesions in the rat. *Brain Res.* 42, 201-206. doi: 10.1016/0006-8993(72)90054-6

Moore-Ede, M. C. (1986). Physiology of the circadian timing system: predictive versus reactive homeostasis. *Am. J. Physiol.* 250, R737-R752. doi: 10.1152/ ajpregu.1986.250.5.R737

Morley-Fletcher, S., Darnaudery, M., Koehl, M., Casolini, P., Van Reeth, O., and Maccari, S. (2003). Prenatal stress in rats predicts immobility behavior in the

forced swim test. effects of a chronic treatment with tianeptine. *Brain Res.* 989, 246-251. doi: 10.1016/S0006-8993(03)03293-1

Nagai, K., Yoshida, S., and Konishi, H. (2014). Gender differences in the gene expression profiles of glucose transporter GLUT class I and SGLT in mouse tissues. *Pharmazie* 69, 856-859.

Nicolaides, N. C., Charmandari, E., Kino, T., and Chrousos, G. P. (2017). Stress-related and circadian secretion and target tissue actions of glucocorticoids: impact on health. *Front. Endocrinol.* 8:70. doi: 10.3389/fendo.2017.00070

Panda, S., Antoch, M. P., Miller, B. H., Su, A. I., Schook, A. B., Straume, M., et al. (2002). Coordinated transcription of key pathways in the mouse by the circadian clock. *Cell* 109, 307-320. doi: 10.1016/S0092-8674(02)00722-5

Reynaert, M. L., Marrocco, J., Mairesse, J., Lionetto, L., Simmaco, M., Deruyter, L., et al. (2016). Hedonic sensitivity to natural rewards is affected by prenatal stress in a sex-dependent manner. *Addict. Biol.* 21, 1072-1085. doi: 10.1111/adb.12270

Reznikov, A. G., and Tarasenko, L. V. (2007). Hormonal protection of gender-related peculiarities of testosterone metabolism in the brain of prenatally stressed rats. *Neuro Endocrinol Lett.* 28, 671-674.

Russell, A. L., Tasker, J. G., Lucion, A. B., Fiedler, J., Munhoz, C. D., Wu, T. J., et al. (2018). Factors promoting vulnerability to dysregulated stress reactivity and stress-related disease. *J. Neuroendocrinol.* 24:e12641. doi: 10.1111/jne.12641

Seckl, J. R. (2008). Glucocorticoids, developmental 'programming' and the risk of affective dysfunction. *Prog. Brain Res.* 167, 17-34. doi: 10.1016/S0079-6123(07) 67002-2

Snedecor, G. W., and Cochran, W. G. (1967). *Statistical Methods*, 6th Edn. Ames Iowa: Iowa State Univ. Press.

Spiga, F., Walker, J. J., Terry, J. R., and Lightman, S. L. (2014). HPA axis-rhythms. *Compr. Physiol.* 4, 1273-1298. doi: 10.1002/cphy.c140003

Summa, K. C., and Turek, F. W. (2011). The genetics of sleep: insight from rodent models. *Sleep Med. Clin.* 6, 141-154. doi: 10.1016/j.jsmc.2011.04.004

Turek, F. W. (2016). Circadian clocks: not your grandfather's clock. *Science* 354, 992-993. doi: 10.1126/science.aal2613

Turek, F. W., Penev, P., Zhang, Y., Van Reeth, O., and Zee, P. (1995). Effects of age on the circadian system. *Neurosci. Biobehav. Rev.* 19, 53-58. doi: 10.1016/0149- 7634(94)00030-5

Vallée, M., Maccari, S., Dellu, F., Simon, H., Le Moal, M., and Mayo, W. (1999). Long-term effects of prenatal stress and postnatal handling on age-related glucocorticoid secretion and cognitive performance: a longitudinal study in the rat. *Eur. J. Neurosci.* 11, 2906-2916. doi: 10.1046/j.1460-9568.1999.00705.x

Van Camp, G., Cigalotti, J., Bouwalerh, H., Mairesse, J., Gatta, E., Palanza, P., et al. (2018). Consequences of a double hit of stress during the perinatal period and midlife in female rats: mismatch or cumulative effect? *Psychoneuroendocrinology* 93, 45-55. doi: 10.1016/j.psyneuen.2018.04.004

Van Reeth, O., Weibel, L., Spiegel, K., Leproult, R., Dugovic, C., and Maccari, S. (2000). Interactions between stress and sleep: from basic research to clinical situations. *Sleep Med. Rev.* 4, 201-219. doi: 10.1053/smrv.1999.0097

Van Someren, E. J. (2000). Circadian and sleep disturbances in the elderly. *Exp. Gerontol.* 35, 1229-1237. doi: 10.1016/S0531-5565(00)00191-1

Vestergaard, M., Holm, S. K., Uldall, P., Siebner, H. R., Paulson, O. B., Baaré, W. F. C., et al. (2017). Glucocorticoid treatment earlier in childhood and adolescence show dose-response associations with diurnal cortisol levels. *Dev. Psychobiol.* 59, 1010-1020. doi: 10.1002/dev.21559

Ward, I. L. (1972). Prenatal stress feminizes and demasculinizes the behavior of males. *Science* 175, 82-84. doi: 10.1126/science.175.4017.82

Ward, I. L. (1984). The prenatal stress syndrome: current status. *Psychoneuroendocrinology* 9, 3-11. doi: 10.1016/0306-4530(84)90016-7

Ward, I. L., and Weisz, J. (1984). Differential effects of maternal stress on circulating levels of corticosterone, progesterone, and testosterone in male and female rat fetuses and their mothers. *Endocrinology* 114, 1635-1644. doi: 10. 1210/endo-114-5-1635

Weaver, I. C., Cervoni, N., Champagne, F. A., D'Alessio, A. C., Sharma, S., Seckl, J. R., et al. (2004). Epigenetic programming by maternal behavior. *Nat. Neurosci.* 7, 847-854. doi: 10.1038/nn1276

Weibel, L., Maccari, S., and Van Reeth, O. (2002). Circadian clock functioning is linked to acute stress reactivity in rats. *J. Biol. Rhythms* 17, 438-446. doi: 10.1177/074873002237138

Weinstock, M. (2001). Alterations induced by gestational stress in brain morphology and behaviour of the offspring. *Prog. Neurobiol.* 65, 427-451. doi: 10.1016/S0301-0082(01)00018-1
Wirz-Justice, A., and Van den Hoofdakker, R. H. (1999). Sleep deprivation in depression: what do we know, where do we go? *Biol. Psychiatry* 46, 445-453.

Wollnik, F., and Turek, F. W. (1988). Estrous correlated modulations of circadian and ultradian wheel-running activity rhythms in LEW/Ztm rats. *Physiol. Behav.* 43, 389-396. doi: 10.1016/0031-9384(88)90204-1

Wulff, K., Gatti, S., Wettstein, J. G., and Foster, R. G. (2010). Sleep and circadian rhythm disruption in psychiatric and neurodegenerative disease. *Nat. Rev. Neurosci.* 11, 589-599. doi: 10.1038/nrn2868

Zhang, R., Lahens, N. F., Ballance, H. I., Hughes, M. E., and Hogenesch, J. B. (2014). A circadian gene expression atlas in mammals: implications for biology and medicine. *Proc. Natl. Acad. Sci. U.S.A.* 111, 16219-16224. doi: 10.1073/pnas.1408886111

Zucker, I., and Beery, A. K. (2010). Males still dominate animal studies. *Nature* 465:690. doi: 10.1038/465690a

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2019 Morley-Fletcher, Mairesse, Van Camp, Reynaert, Gatta, Marrocco, Bouwalerh, Nicoletti and Maccari. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.