Original Research ARTICLE
Perinatal Stress Programs Sex Differences in the Behavioral and Molecular Chronobiological Profile of Rats Maintained Under a 12-h Light-Dark Cycle
- 1UMR 8576, Unité de Glycobiologie Structurale et Fonctionnelle, Campus Cité Scientifique, CNRS, University of Lille, Lille, France
- 2University Lille – CNRS-UMR 8576, International Associated Laboratory (LIA) “Prenatal Stress and Neurodegenerative Diseases,” Sapienza University of Rome – IRCCS Neuromed, Rome, Italy
- 3Division of Neonatology, Department of Pediatrics, University of Geneva, Geneva, Switzerland
- 4Department of Psychiatry, College of Medicine, Psychiatric Institute, Center for Alcohol Research in Epigenetics, University of Illinois at Chicago, Chicago, IL, United States
- 5Laboratory of Neuroendocrinology, The Rockefeller University, New York, NY, United States
- 6Department of Physiology and Pharmacology “V. Erspamer,” Sapienza University of Rome, Rome, Italy
- 7Istituto di Ricovero e Cura a Carattere Scientifico (IRCCS), NEUROMED, Pozzilli, Italy
- 8Department of Science and Medical – Surgical Biotechnology, Sapienza University of Rome, Rome, Italy
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.
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
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
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.
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 occurr 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 (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).
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 StepOnePlusTM software. A ΔCt (Ct of the considered gene – mean Ct (Gapdh and Ywhaz genes) was calculated.
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).
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).
Figure 1. Circadian locomotor activity in a wheel running paradigm in PRS and CONT animals of both sexes. (A) The 24 h representation of the mean number of wheel revolutions calculated over 10 consecutive days and plotted with a 30 min resolution, obtained from adult PRS and CONT rats of both sexes showing hyperactivity in PRS males and hypoactivity in PRS females. (B) Details of the light switch off period with data plotted with a 5 min resolution allowing a clear visualization of the advance in activity phase in male PRS animals. (C) Total number of wheel revolutions over 10 consecutive days during the period of inactivity (light phase) and activity (dark phase) showing the sex effect in the CONT group with females being more active in the light and less active in the dark phase. This effect was observed in the PRS group only in the dark phase, with females being hypoactive, and PRS males being always hyperactive. (D) Representative activity records of the circadian running wheel activity of male and female PRS and CONT rats over 10 consecutive days, where the phase advance in PRS males and fragmented activity in both PRS males and females is clearly observed. The first L/D cycle is shown at the top of each panel. Each horizontal line represents 24 h of the animal’s life composed from the vertical bar, each representing a wheel revolution. Values are expressed as means ± 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.
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).
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.
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, Mel1aR, 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 Mel1aR 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; Mel1aR, F(2,32) = 2.811, p = 0.07; Rev-erbA-α, F(2,32) = 7.901, p < 0.01).
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.
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.
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.
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 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 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 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. Down-regulation 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.
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.
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.
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.
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.
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.
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; Gsk3ß 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.
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
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
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. Rhyth. 30, 417–436. doi: 10.1177/0748730415598608
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
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
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
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
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
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
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
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
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
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., 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
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
Keywords: predictive adaptation, reactive adaptation, circadian rhythms, locomotor activity, chronobiological stressor, mRNA expression
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
Received: 28 November 2018; Accepted: 21 March 2019;
Published: 01 May 2019.
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
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.
*Correspondence: Stefania Maccari, firstname.lastname@example.org
†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