Skip to content


  • Research article
  • Open Access

Biological clock function is linked to proactive and reactive personality types

BMC Biology201816:148

  • Received: 16 October 2018
  • Accepted: 4 December 2018
  • Published:



Many physiological processes in our body are controlled by the biological clock and show circadian rhythmicity. It is generally accepted that a robust rhythm is a prerequisite for optimal functioning and that a lack of rhythmicity can contribute to the pathogenesis of various diseases. Here, we tested in a heterogeneous laboratory zebrafish population whether and how variation in the rhythmicity of the biological clock is associated with the coping styles of individual animals, as assessed in a behavioural assay to reliably measure this along a continuum between proactive and reactive extremes.


Using RNA sequencing on brain samples, we demonstrated a prominent difference in the expression level of genes involved in the biological clock between proactive and reactive individuals. Subsequently, we tested whether this correlation between gene expression and coping style was due to a consistent change in the level of clock gene expression or to a phase shift or to altered amplitude of the circadian rhythm of gene expression. Our data show a remarkable individual variation in amplitude of the clock gene expression rhythms, which was also reflected in the fluctuating concentrations of melatonin and cortisol, and locomotor activity. This variation in rhythmicity showed a strong correlation with the coping style of the individual, ranging from robust rhythms with large amplitudes in proactive fish to a complete absence of rhythmicity in reactive fish. The rhythmicity of the proactive fish decreased when challenged with constant light conditions whereas the rhythmicity of reactive individuals was not altered.


These results shed new light on the role of the biological clock by demonstrating that large variation in circadian rhythmicity of individuals may occur within populations. The observed correlation between coping style and circadian rhythmicity suggests that the level of rhythmicity forms an integral part of proactive or reactive coping styles.


  • Circadian rhythm
  • Clock genes
  • Behaviour
  • Coping style
  • Endocrine regulation
  • Cortisol
  • Melatonin
  • RNA sequencing
  • Zebrafish


In humans, personality comprises consistent individual variation of correlated behavioural and habitual traits [1], such as extraversiveness, impulsivity or novelty-seeking [2, 3]. These personality traits are also linked to the physiological and genetic makeup of an individual [48]. Similarly, in animals, behavioural traits such as risk-taking or aggressiveness [9] can be consistently correlated with physiological traits such as metabolic rate [10] or endocrine stress response [11]. These correlated sets of traits, termed coping styles (animal personalities), are recognised in a number of taxa [12] and vary along a proactive-reactive continuum, with proactive individuals being more risk-taker and aggressive and having higher baseline metabolic rates than reactive individuals [9, 10]. Further, characterisation of coping styles is crucial to a fundamental understanding of the occurrence of behavioural polymorphisms [1315] and may provide new perspectives on the application of personalised care and personality disorders in humans [6, 1621].

In experimental animals, coping styles can be assessed using the degree of risk-taking behaviour as a proxy, which is strongly correlated with a number of other behavioural and physiological traits [912]. In various fish species, risk-taking can be evaluated by an emergence test, a well-established test with high repeatability [911, 22]. During this test, a group of individual fish is allowed to emerge from a familiar shelter into a novel and potentially dangerous environment. The order of emergence is considered a measure for the individual tendency of risk-taking [911, 22].

In the present study, we have used this assay to assess the coping styles of individual zebrafish from a heterogeneous laboratory population (AB/TL strain). Subsequently, RNA sequencing was performed on brain samples from proactive and reactive zebrafish to further investigate molecular determinants of coping styles. Our results showed an overrepresentation of genes involved in the biological clock among the genes differentially expressed in the brains from proactive and reactive individuals. The biological clock is an endogenous timing mechanism which allows an organism to anticipate regular changes in the environment that result from the day/night cycle [2327]. In vertebrates, the molecular basis of the biological clock is formed by feedback loops of gene expression. In the core loop, the transcription factors BMAL and CLOCK stimulate the expression of the genes encoding the corepressor proteins PER and CRY, which inhibit the activity of BMAL and CLOCK. Thereby, PER and CRY repress their own expression and create an oscillating pattern of gene expression. This oscillator is stabilised by a second feedback loop which involves the nuclear receptors REV-ERB and ROR. In addition, the biological clock can be synchronised with the actual day-night cycle by environmental cues like light and temperature. The suprachiasmatic nucleus (SCN) of the hypothalamus in mammals and the pineal gland in other vertebrates appear to act as ‘master clocks’. They coordinate the clock of peripheral cells and tissues through a set of systemic signals, including the cycling secretion patterns of the hormones melatonin and cortisol. This way, the biological clock creates diurnal rhythmicity which can be observed at several functional levels, from the oscillating expression of clock genes over endocrine secretion patterns to rhythmic behavioural activity [2327].

The association between clock gene expression levels and coping styles observed in our RNAseq experiments led to three alternative hypotheses. The different gene expression levels may be generally higher or lower between individuals with different coping styles; the observed differences in gene expression may be due to a shift in phase of the biological clock (i.e. chronotype [4, 2325, 28]) or to a different amplitude of the diurnal rhythm. We tested these hypotheses on three levels of biological function: on the molecular level by measuring gene expression, on the endocrine level by investigating hormone production and on the behavioural level by studying locomotion activity. Our results show that the amplitude of the circadian rhythm correlates with the coping style of the individual, ranging between high amplitudes in proactive fish and an absence of rhythmicity in reactive fish.


The emergence test as a proxy for coping style

We assessed different coping styles using a group emergence test, in which each zebrafish of a group of ten was allowed to emerge from a familiar (slightly darkened) shelter into a novel and potentially dangerous environment. The individual fish were ranked according to their latency to emerge (Fig. 1a). We evaluated the temporal consistency of this emergence test and tested whether difference in light intensity was a confounding factor. First, the test was repeated in a single emergence setting and emergence times were measured. A positive correlation between single emergence times and group emergence rank was observed (Fig. 1b; Spearman rank, N = 187, ρ = 0.82, p < 0.0001). Second, single emergence tests were performed on two consecutive days at different times of the day. The results of this experiment show strong correlations between single emergence times of the same individual during two consecutive days at different times of the day (Additional file 1: Figure S1A; Spearman rank, N = 19–24, ρ = 0.38–0.91, p < 0.001 or p < 0.0001). Third, to test the effect of light conditions on emergence behaviour, the emergence test was performed with similar light conditions in both the shelter and the novel environment compartment. A strong correlation was observed between the emergence time measured using this approach and the emergence rank in the initial test (Additional file 1: Figure S1B; Spearman rank, N = 72, ρ = 0.77, p < 0.0001). Fourth, when plotting all pooled emergence times from single emergence experiments over the group emergence times, a strong positive correlation confirmed the robustness of emergence time regardless the social and environmental settings (Additional file 1: Figure S1C; Spearman rank, N = 190, ρ = 0.87, p < 0.0001). In order to further evaluate risk-taking as a proxy for coping style, we tested and confirmed a correlation between emergence rank and frequency of aggressive behaviour during a mirror-biting test (Fig. 1c; Spearman rank, N = 59, ρ = − 0.77, p < 0.0001).
Fig. 1
Fig. 1

Behavioural testing in zebrafish. Behavioural testing in zebrafish. a Schematic overview of the group emergence test used in this study. A fish tank was used that consisted of two compartments separated by a wall with a closable hatch. A batch of ten fish was introduced into the darkened holding compartment of this tank. After a 10-min acclimation period, the hatch was opened and fish were allowed to emerge into the uncovered, well lit, second compartment, and the emergence time was recorded. After emergence, individual fish were collected and grouped according to emergence rank (1–10), with number 1 designated early emerger (EE) and number 10 late emerger (LE). The emergence order is considered a measure for risk-taking behaviour, which is a widely used proxy for coping style. b Behavioural traits are correlated across situation and time: single emergence time plotted against group emergence ranks. Data were collected from 24 batches of 10 adult zebrafish (N = 187) in the group emergence test (rank) and subsequent single emergence test (time), with a significant correlation between emergence rank in the group test and time in the single test (Spearman rank test). c Aggressiveness correlates with risk-taking behaviour within a coping style: individuals previously ranked during a group emergence test were subjected to a mirror-image stimulation. Upon exposure to their own mirror image, the number of aggressive behaviours (AGR; bites to the image, parallel swimming, circles and strikes) were counted and divided by the measuring period minus the duration of freezing bouts and approach latency. The resulting AGR frequency (s−1) was significantly correlated with emergence rank (Spearman rank test, N = 59)

Since fish from the hybrid AB/TL line were used in these experiments, it could be hypothesised that the observed behavioural variation in this line originated from genetic differences between the original AB and TL lines. Since these two lines have previously been shown to perform differently in other behavioural tests [29], proactivity may have been inherited from one line and reactivity from the other. However, no difference was observed in the behaviour of fish from the AB and TL lines in the emergence test (Additional file 2: Figure S2; X2 test, N = 6, p > 0.05). We therefore conclude that genetic variation in the hybrid AB/TL line which may underlie the observed phenotypical variation most likely originates from genetic variation within the original AB and TL strains, [30, 31]. However, the contribution of genetic variation between these lines (e.g. involving different genes affecting the phenotype in the two lines) can still not be excluded [30, 32].

Diurnal rhythmicity of gene expression is associated with coping style

Subsequently, group emergence tests were performed with ten individuals and we collected individuals of rank 1 (early emergers, EE, highly risk-taking) and rank 10 (late emergers, LE, highly risk-avoiding). RNA was isolated from their brains and processed for transcriptome analysis by RNA sequencing. This yielded RNA sequences derived from 31,398 different genes. The expression of 1478 of these genes differed significantly between EE and LE fish (Additional file 3: Table S1). Gene ontology analysis revealed that 20 of the top 100 most differentially expressed genes were involved in the control of the biological clock (and 43 of the total of 1478 differentially expressed genes, enrichment score 3.66; Fig. 2 and Additional file 3: Table S1), with a fold change ranging from 4.34 to 20.07. These data indicate a strong association between coping styles and the regulation of the biological clock. Other enriched ontology groups to be found in the top 100 most differentially expressed genes were ‘oxygen binding and transport’ (enrichment score 5.81), ‘transcription factor activity’ (3.38) and ‘haematopoiesis’ (3.30).
Fig. 2
Fig. 2

Transcriptome analysis of brain samples from proactive (early emerging, EE) and reactive (late emerging, LE) zebrafish using RNA sequencing. Brains were collected from EE and LE fish (n = 4) at 7 hCT, and subsequently, RNA was isolated and used for RNA sequencing. Analysis of the EE and LE transcriptomes showed that differences in gene expression were found for 3% of genes involved in the regulation of the biological clock (see Additional file 3: Table S1). This overview shows differential gene expression levels of 44 differently (p < 0.05) regulated genes, involved in the core loop, the stabilising loop and the entrainment of the biological clock. Rectangles represent genes, and differential gene expression between EE and LE fish is indicated by red (higher expression level in EE) and blue rectangles (higher level in LE fish). Circles represent regulating protein complexes. Inhibition is indicated by a red line with round head, and promotion by a black line with arrow head

Subsequently, we further investigated the relationship between coping style and the biological clock at different functional levels, i.e. gene expression, endocrinology and locomotion behaviour. First, we measured, by quantitative PCR, the expression patterns of 5 genes, i.e. one member from each of the four core clock genes and one stabilising loop gene that showed diurnal rhythmicity in its expression (bmal1a, clock1a, per1a, cry1a and cipca) and that were among the top 20 most differentially expressed clock genes. The observed expression patterns show that EE fish had higher peak values and a wider range of expression levels of these genes than LE fish (Fig. 3; Two-way ANOVA, Bonferroni post hoc, N = 6, p < 0.05). When fish of all emergence ranks were taken into account (i.e. 1 to 10, Additional file 4: Figure S3), both the area under the curve (AUC, as a measure for net transcriptional activity) and the amplitude of expression (used as an indicator for rhythmicity) were negatively correlated with the emergence rank, for bmal1a and cipca regarding AUC (Spearman rank test, N = 72, ρ = − 0.4619, p = 0.0102 and ρ = − 0.4793, p = 0.0074, respectively) and all genes regarding amplitude (Spearman rank test, N = 72, ρ = − 0.67 to − 0.45, p < 0.0001 to p = 0.0138). These results indicate a more robust circadian rhythmicity in individuals with proactive coping styles compared to individuals with reactive coping styles.
Fig. 3
Fig. 3

The expression patterns of five clock-related genes over time are correlated with the risk-taking behaviour. The brains were collected from EE and LE fish, at six different time points, i.e. 1, 7, 13, 15, 19 and 23 hCT, in a 14–10-h light to dark regime (n = 3), and subsequently, RNA was isolated and used for qPCR analysis. The relative expression of bmal1a, clock1a, per1a, cry1a and cipca was double-plotted over time from EE (red) and LE (blue) fish. Letters indicate significant difference (two-way ANOVA, with time and coping style as variables, Bonferroni post hoc test, p < 0.05, N = 6. Data shown are means ± SEM)

Diurnal rhythmicity of hormone secretion is associated with coping style

Subsequently, we studied the diurnal rhythmicity of levels of cortisol and melatonin, important endocrine regulators of diurnal activity patterns [23]. In EE fish, whole body melatonin and cortisol concentrations fluctuated significantly between light and dark phase, whereas LE fish had constantly high values for cortisol and intermediate values for melatonin, with no detectable rhythmicity (Fig. 4; two-way ANOVA, with time and coping style as variables, Bonferroni post hoc test, N = 6, p < 0.05). Additionally, taking the data of all emergence ranks into account, a positive correlation between emergence rank and net secretion activity (AUC), and a negative correlation between emergence rank and circadian rhythmicity (amplitude), was observed for both hormones (Additional file 5: Figure S4; Spearman rank test, N = 72; cortisol: ρ = 0.81 and p < 0.0001 (AUC), ρ = − 0.64 and p = 0.0002 (amplitude); melatonin: ρ = 0.61 and p = 0.0003 (AUC), ρ = − 0.37 and p = 0.0429 (amplitude)). This indicates a stronger rhythmicity in the secretion of melatonin and cortisol in proactive individuals compared to reactive individuals.
Fig. 4
Fig. 4

Diurnal patterns in cortisol and melatonin levels are associated with the risk-taking behaviour. Full body melatonin (a) and cortisol (b) concentrations (pg g−1 and ng g−1 body weight, respectively) were determined by ELISA in early (EE, red) and late emerging fish (LE, blue) over 24 h during a 14:10 light to dark cycle (indicated by white and grey bars, double plotted). Samples were taken at 1, 7, 13, 15, 19 and 23 h circadian time (hCT). Data shown are means ± SE (N = 3). Letters indicate significant difference (two-way ANOVA, with time and coping style as variables, Bonferroni post hoc test, p < 0.05, N = 6. Data shown are means ± SEM)

Diurnal rhythmicity of behavioural activity is associated with coping style

Finally, we monitored the locomotor activity during three consecutive diurnal cycles. Plotting average swimming velocity data of EE and LE fish revealed that EE fish reached higher peak activity levels (especially in the first hours of the light phase) and a wider range of activity levels than LE fish (Fig. 5a; two-way ANOVA with time and coping style as variables, Sidack’s post hoc test, N = 6, p < 0.05). Further analysis (Spearman rank) of the activity data using fish of all emergence ranks showed that emergence rank was negatively correlated with rhythmicity (Additional file 6: Figure S5A, N = 72, ρ = − 0.34 and p < 0.0037 (AUC)), total behavioural activity (Additional file 6: Figure S5B, N = 72, ρ = − 0.39 and p < 0.0008 (amplitude)) and rhythm strength (Additional file 6: Figure S5C, N = 72, ρ = − 0.53 and p < 0.0001). However, maximum swimming speeds did not vary significantly with emergence ranks, suggesting that there was no difference in physiological capacity for locomotion between fish of different emergence ranks (Additional file 6: Figure S5D). Again, these data demonstrate a more pronounced rhythmicity in the more proactive individuals.
Fig. 5
Fig. 5

Diurnal activity pattern is associated with emergence behaviour. Actogram of early (EE, emergence rank 1, red) and late emerging fish (LE, rank 10, blue) during 3 days under light-dark regime (a) and during three subsequent days in constant light (b). Locomotion activity is represented as average swimming velocities over 30 min (V, mm s−1) as a function of time (hours circadian time, hCT). Activity recording started at 7 hCT of day 1 and ended at 7 hCT of day 6. The light regime switched at day 3 from light-dark (14:10) cycles to constant light, indicated by white and grey bars. Asterisks and horizontal bars indicate significant different time intervals between coping styles (two-way ANOVA with time and coping style as variables, Sidack’s post hoc test, N = 6, p < 0.05. Data shown are means ± SEM)

Finally, we subjected the fish to a constant light challenge and monitored their behaviour over three consecutive diurnal cycles. Under these constant light conditions, the differences between EE and LE fish in maximum values and range of activity level gradually decreased, with EE fish approaching the activity pattern observed for LE fish (Fig. 5b; two-way ANOVA with time and coping style as variables, Sidack’s post hoc test, N = 6, p < 0.05). The activity data of fish of all emergence ranks showed that rhythm strength was correlated (Spearman rank) with emergence rank (Additional file 6: Figure S5G; N = 72, ρ = − 0.61, p < 0.0001), but not total activity or amplitude (Additional file 6: Figure S5E-F, N = 72, p > 0.05). These results indicate that activity patterns of proactive individuals gradually lose rhythmicity in the absence of time-related cues and become similar to the arrhythmic behaviour of reactive individuals. This demonstrates that the strong diurnal rhythmicity of proactive individuals is dependent on external light cues.


In the present study, we have classified individual zebrafish from the same batch by their coping style, along a proactive-reactive continuum, using risk-taking behaviour as a proxy for coping style. Subsequently, we have demonstrated a strong association between coping style and the amplitude of the circadian rhythmicity at the molecular, endocrine and behavioural level: proactive individuals showed a robust diurnal rhythm with a large amplitude, whereas reactive individuals virtually lacked rhythmicity.

Transcriptome analysis using RNA sequencing on brain samples from proactive and reactive individuals showed an enrichment of genes involved in the regulation of the biological clock among the set of genes that was differentially expressed. In several recent studies, transcriptome analysis has also been performed to study differential gene expression related to coping styles in zebrafish [9, 33]. In one of these studies [9], an approach to assess coping styles was used which was very similar to ours, and the RNA sequencing analysis on hindbrain samples demonstrated that genes involved in the biological clock were overrepresented among the differentially expressed genes. However, in other brain regions, no correlation between clock-related genes and coping style was observed. In another study [34], a microarray analysis was performed on whole brain samples from zebrafish strains with different behavioural phenotypes, which also revealed differential expression of biological clock-related genes.

The results from our RNA sequencing raised two issues, which we answered in additional control experiments. First, the time of the day may have affected the separation of behavioural phenotypes during the emergence test, and the observed differences in coping style could have been a result of phase-shifted diurnal rhythmicities (i.e. chronotypes [4, 2325, 28]). However, the performance in the emergence test of individual fish at different times on two consecutive days appeared to be highly consistent, indicating that the observed behaviour in the emergence test was independent of the time of the day. Second, many of the differentially expressed clock genes are known to be light-sensitive [35]. Therefore, the issue arose whether light responsiveness was a confounding factor in the emergence test, since in the initial setup there was a difference in light intensity between the two compartments. When the test was repeated with identical light conditions on both sides of the hatch, the emergence order was not altered, which demonstrated that the emergence order was not affected by differences in light responsiveness between individual fish.

In addition, we tested the validity of risk-taking behaviour in the emergence test as a proxy for coping style by correlating it with aggressiveness measured using a mirror test. We found a strong correlation between these two behavioural traits, which is in line with previous studies using an evolutionary model [13] and an experimental approach [36] and has been reviewed elaborately [14, 37, 38]. The results of our extensive validation of the emergence test support the robustness of the emergence test and explain its wide use in coping style research [911, 22].

In our population of zebrafish, we found large individual variation in diurnal rhythmicity within a continuum that ranges from very robust diurnal rhythms to a complete absence of rhythmicity. The spontaneous occurrence of individuals with a low level of rhythmicity within a healthy laboratory population has been shown before in Djungarian hamsters [39], and it will be interesting to study whether such individuals can be found under natural conditions as well. Under certain conditions, a lack of an internal rhythmicity can be advantageous for individuals, for example when external cues are absent or fluctuate unpredictably [2325, 28] or when competition with other individuals can be avoided by a larger flexibility in the timing of certain activities [40]. In humans, the enormous variation observed in the level and rhythm of the melatonin secretion [4143] suggests that variability in circadian rhythmicity may be a more general phenomenon than previously assumed. Recognition of this diversity in individual phenotypes [1215] will be of great relevance for our understanding of the regulation and function of the biological clock, but also for the personalization of diagnosis and treatment of psychiatric and metabolic disorders [2325, 28, 44] and degenerative diseases [45] that are attributed to dysregulation of the biological clock.

The observed variation of circadian rhythmicity was associated with coping style, with proactive individuals showing robust rhythms and reactive fish displaying very low rhythmicity. We therefore suggest the diurnal rhythmicity to be an integral part of the coping style of an individual. In general, proactive individuals have been shown to be mainly intrinsically organised, i.e. routine-based and relatively resistant to the influence of external stimuli [1214]. The biological clock itself is a proactive system, anticipating rhythmic changes in external conditions, so the stronger rhythmicity found in proactive zebrafish is in line with a more intrinsic regulation of diurnal rhythmicity. In contrast, for reactive individuals, a robust diurnal rhythm may be constraining the flexibility that characterises their reactive phenotype and would therefore be maladaptive in combination with this behavioural phenotype.

Although future research should indicate how generally this phenomenon occurs in nature, our study demonstrates that there are mechanisms that under certain circumstances link coping style and circadian rhythmicity. Possibly, this correlation between coping style and circadian rhythmicity may be a result of genetic correlation, i.e. one trait responding as a consequence of selection acting upon another trait with which it shares genes [46]. The main mechanism underlying genetic correlation is pleiotropy, in which one genetic variant affects more than one phenotype [47]. Genes related to the biological clock are often pleiotropic genes, which can influence both behaviour and physiology [2325, 27, 28, 48]. Other possible pleiotropic genes that are involved in the observed correlation between coping style and the biological clock are those involved in the biosynthesis and function of neurotransmitters and hormones. For example, cortisol has been found to be a ‘permissive cue’ for diurnal rhythmicity through a variety of molecular mechanism [49]. Additionally, glucocorticoids are responsible for many changes in behaviour related to coping style in both humans [50] and non-human vertebrates [51]. Besides this direct form of pleiotropy, an indirect form of pleiotropy exists, in which a genetic variant affects one phenotype, which lies on the causal path to a second phenotype. In our model, changes in circadian activity patterns may have dramatic effects on certain endocrine systems [52, 53], which in turn may influence other behavioural traits.


It is generally accepted that a disturbed or weak diurnal rhythm compromises optimal biological functioning and is linked to psychiatric and metabolic disorders and degenerative diseases. However, here, we demonstrate large variation in diurnal rhythmicity within a single laboratory population of zebrafish, ranging from robust rhythms to a complete absence of rhythmicity. This variation is correlated with the coping style of the individual fish. This association suggests that the biological clock may be part of the complex phenotypes that represent individual survival strategies in nature as well.

Materials and methods

Zebrafish husbandry conditions

Zebrafish were maintained and handled according to the guidelines from the Zebrafish Model Organism Database (ZFIN, and in compliance with the directives of the local animal welfare committee of Leiden University (DEC number 14058). The wildtype zebrafish (Danio rerio) strain used in this study was an AB/TL strain, which originates from crossbreeding of the AB and Tüpfel Long Fin (TL) strains. This strain was originally obtained from the Hubrecht Laboratory (Utrecht, The Netherlands) and maintained in our laboratory for at least ten generations at the time of the experiments described here. New generations were generated by mating of fish from the previous generation, using 60–80 individuals (1:1 ratio male to female). The AB/TL line, considered a segregated hybrid line, was chosen because it contains large genetic variation. In addition to the genetic variation existing within the original AB and TL strains [30, 31], the variation between the two original lines [30, 32] also contributes to the genetic variation of the AB/TL line.

The fish were reared in densities of ± 40 individuals (male to female 1:1) per 7.5-l tanks in standardised recirculation systems (Fleuren & Nooijen, Nederweert, The Netherlands); water temperature was maintained at 28 ± 1 °C (n = 5), with a conductivity of 518 ± 12 μS (n = 5) and oxygen concentration of 7.9 ± 0.4 mg l−1 (n = 5). Light cycles were maintained at 14-h light to 10-h darkness cycle, with light periods from 8:00/7:00 (0 h Circadian Time, (hCT)) to 22:00/21:00 (14 hCT) summer time/winter time, with a linearly decreasing/increasing light intensity between 0 and 320 ± 21 ln m−2 (N = 3) over a period of 15 min. Fish were fed twice daily, at 1 ± 1 hCT and at 8 ± 1 hCT, with dry food (DuplaRinM, Gelsdorf, Germany) and frozen artemia (Dutch Select Food, Aquadistri BV, Klundert, The Netherlands). The fish used in the experiments were between 1 and 2 years old and had a standard length of 32.1 ± 2.3 mm and a body weight of 150.61 ± 17.99 mg (mean ± SD). There was no correlation between emergence rank and body weight (Spearman rank, N = 144, p > 0.05).

Behavioural tests

Group emergence

The setup for the emergence test consisted of a Plexiglas tank (33 × 13 × 13 cm) with a volume of 3 l. It was divided into a darkened holding compartment and an uncovered novel area compartment, by a wall with a hatch (2 × 2 cm) at its mid bottom. The hatch was manually closable by means of a trap door. Ten fish of mixed sex were transferred from the housing tank into the holding compartment of the emergence setup where they were acclimated for 10 min. Thereafter, the trap door was opened, enabling the emergence of one fish at a time into the novel area compartment. After each emergence event, the trap door was closed and the emerged fish was manually transferred to holding tanks (33 × 13 × 13 cm), separated by emergence rank (1–10). The time of emergence from the moment of the opening of the trap door until the actual passage through the hatch was recorded for all individuals. The entire test did not last longer than 10 min and was performed between 5:00 to 8:00 hCT, unless otherwise indicated. Fish not emerging after 10 min were excluded from further analysis. Animals were kept overnight in the holding tanks until further experimentation the next day.

Individual emergence

In order to test the consistency of the degree of risk-taking over the diurnal cycle, individual emergence was investigated at three times of the day (1:00 h, 7:00 h and 13:00 hCT). Nine groups (n = 8) were created, each tested at a combination of two different times of the day on consecutive days (1:00/1:00, 1:00/7:00, 1:00/13:00, 7:00/1:00, 7:00/7:00, 7:00/13:00, 13:00/1:00, 13:00/7:00, 13:00/13:00 hCT). The individual emergence setup consisted of three vertical rows of eight Plexiglass emergence tanks, as used for the group emergence tests, with the trap doors connected per row. A flow through system (47.6 l h−1 per tank) provided each tank with fresh system water, and LED strips (Eurolite, Steinigke Showtechnic GmbH, Waldbüttelbrunn, Germany) were placed in front of the uncovered compartment, resulting in similar light intensity as in the husbandry facility (425 ± 11 ln m−2). The light intensity of the holding compartment was 5 ± 1 ln m−2. This setup enabled testing the emergence of eight individuals simultaneously. Fish were placed individually in the holding compartment and acclimated overnight. At testing time, hatches were opened manually and emergence was recorded using a video camera (HDC-SD90, Panansonic Corporation, Kadoma Osaka, Japan) placed 1.2 m in front of the setup. The test was terminated when all fish had emerged. Fish were placed back into the holding compartment, and the test was repeated on the same individual the following day, according to the design described above.

In order to test the consistency of risk-taking behaviour at identical light intensities of the holding compartment and the novel environment compartment of the emergence test, a similar setup was used with light intensities of 5 ± 1 ln m−2 in both compartments. Here too, nine groups with eight fish each (N = 72) were created, each tested at 7:00 h ZT in the individual emergence setup according to the protocol described above.

A group emergence test was performed 1 day prior to the individual emergence tests to investigate the correlation between the behaviour in the two tests. Each row of 8 tanks in the individual emergence setup contained fish from the same group emergence test. Since only 8 out of the 10 individuals could be tested simultaneously, two individuals were randomly eliminated, using the excel randomisation function. The group emergence test was performed nine times, yielding 72 fish (9 groups of 8) to be tested in the individual emergence test. As indicated above, each group was tested at different times on two consecutive days. This combination of individual and group emergence test was performed three times, so a total of 216 fish was tested. Individual emergence times were determined by the observation from the video footage and were correlated with the emergence time of the group emergence test using Spearman rank correlation with significance accepted at p < 0.05.

Emergence of mixed AB and TL groups

In order to test whether the observed behavioural variation in the degree of risk-taking originated from the difference between the AB and TL lines used to create the wildtype (AB/TL) line of this study, the emergence test was repeated with groups created by mixing AB and TL fish in a ratio of 1:1. Six groups of ten fish were created before the group emergence test was performed, and housed under the same conditions as previously described. The group emergence test was subsequently performed for each group according to the protocol described above. Data analysis was performed using a Χ2 test (N = 6, significance accepted at p < 0.05).

Mirror-image stimulation

For the mirror-image stimulation (15), three rows of eight Plexiglas tanks (33 × 13 × 13 cm, capacity 3 l) were placed on a white table surface (1.0 × 1.0 m). A HD video camera (Panasonic, HDC-SD90, Panasonic Inc., Japan) was mounted 1.7 m above the table. Fish previously tested in the group emergence test (ca. 2 hCT) in sets of ten (with random elimination of two) were placed individually in a circular holding chamber, consisting of a grey tube (10 cm diameter, 15 cm height) placed vertically in the middle of the tanks. A vertical mirror (13 × 15 cm) was placed at one of the short ends of the tank. Each row contained fish from the same group emergence test, i.e. 8 fish in three rows, resulting in 24 fish per experimental run, with three experimental runs. Fish were acclimated 30 min prior to recording, after which the holding chambers were removed simultaneously from all 24 tanks in a swift vertical movement, exposing them to the open field and the mirror, and their behaviour was recorded for 10 min (ca. 4 hCT). Measured variables were the duration of freezing bouts (FRZ; complete lack of movement only for eyes and gills), the latency to realise the first approach to the mirror until a distance of 1 body length, i.e. ca. 4 cm, i.e. 1 body length (LFA), and the number of aggressive behaviours (AGR; bites to the image, parallel swimming, circles and strikes). The variables were determined by the observation from the video footage. During subsequent analysis, AGR were counted and divided by the measuring period minus the FRZ and LFA, therefore regarding only the period dedicated to potentially aggressive interactions with the mirror image. Non-responsive fish, i.e. fish with FRZ of 10 min or longer, were excluded from the analysis. The correlation between the resulting AGR frequency (s−1) and the emergence time of the group emergence test was evaluated using Spearman rank correlation (n = 59, N = 9, significance accepted at p < 0.05).

Transcriptional analysis

RNA sequencing

In order to establish a first estimation of the range of RNA expression levels over emergence rank variation, four individuals of emergence rank 1 (EE) and four individuals of emergence rank 10 (LE) were sacrificed in ice water, at 7 ± 0.5 hCT. Their brains were dissected out and preserved in RNALater (Ambion, Austin, TX). Brains were homogenised in 500 μl of Trizol reagent (Qiagen). Total RNA was extracted and column-purified using the RNeasy MinElute Cleanup Kit (Qiagen), according to the manufacturer’s instructions. RNA samples were treated with DNaseI (Life Technologies) to remove residual genomic DNA. RNA integrity was analysed by Lab-on-a-chip analysis (Agilent, Amstelveen, The Netherlands). A total of 2 μg of RNA was used to make RNAseq libraries using the Illumina TruSeq RNA Sample Preparation Kit v2 (Illumina, Inc., San Diego, CA, USA) according to the manufacturer’s instructions (with two minor modifications: in the adapter ligation step, 1 μl, instead of 2.5 μl, adaptor was used, and in the library size selection step, the library fragments were isolated with a double Ampure XP purification with a 0.7× beads to library ratio (Beckman Coulter, Woerden, The Netherlands)). The resulting mRNAseq library was sequenced using an Illumina HiSeq2500 Instrument (Illumina, Inc.) with a read length of 2 × 50 nucleotides. Image analysis and base-calling were done using the Illumina HCS version 2.0.12. The raw data has been submitted to the GEO database (accession number GSE64570). The data was analysed using the GeneTiles software ( with a cutoff p value of 0.05. In brief, Genetiles used fastq files as input for the program Bowtie2 ( to align the reads to the zebrafish genome (obtained from Ensembl version Zv9). Subsequently, the programs SAMtools (, DESeq2 and DESeq ( were used for data processing. Gene ontology of the significantly (p < 0.05) different regulated 1478 genes was analysed using the online functional classification tool DAVID ( In addition, for genes not classified by DAVID, information was gathered on their function, using the websites GeneCards (, NCBI ( and Genetics Home Reference ( Using this information, all genes involved in the biological clock were identified. In the “Results” section, we mention the enrichment scores of the three most enriched functional annotation clusters for the top hundred most differentially expressed genes, assigned by DAVID.

Quantitative PCR

For the estimation of transcriptomic regulation, one member from each core clock gene family was selected (bmal1a, clock1a, per1a, cry1a), and one gene involved in the stabilising loop that showed diurnal rhythmicity in its expression (cipca). The transcriptional regulation of these genes was measured in individuals from all emergence ranks (1–10). After the group emergence test, random elimination of two fish per group and overnight group housing separated by emergence order, triplicates of 3 × 8 fish (n = 72) were housed individually in separate containers (10 × 10 × 10 cm). Sampling occurred at six time points, i.e. in the middle of the light phase (7:00 hCT), 1 h before the dark phase (13:00 hCT), 1 h into the dark phase (15:00 hCT), in the middle of the dark phase (19:00 hCT), 1 h before the light phase (23:00 hCT) and 1 h into the light phase (1:00 hCT). Fish were sacrificed in ice water, and brains were dissected and preserved in RNALater (Ambion, Austin, TX). The bodies were preserved in liquid nitrogen for further analysis of hormones (see below). Triplicate samples were taken. The tissue was homogenised in 1 ml TRIzol reagent (Invitrogen), and total RNA was extracted according to the manufacturer’s instructions. cDNA synthesis was performed using the iScript kit (Bio-Rad Laboratories) according to the manufacturer’s instructions with forward and reverse primers for the aforementioned target genes and ippA (inducing PCN production A) as a reference (housekeeping) gene. Real-time quantitative PCR was performed using the Chromo4 Real-time PCR detection system (Bio-Rad Laboratories). All reactions were performed as technical triplicates. PCR analysis was performed using the following protocol: 95 °C 3 min, 40 cycles of 95 °C 15 s and 60 °C 45 s, and final melting curve of 81 cycles from 95 °C 1 min to 55 °C 10 s. Results were analysed using the ΔΔCt method. In short, ΔCt was calculated for each gene by subtracting the Ct values determined for the housekeeping gene from Ct values determined for the target genes. Then, the ΔCt values for EE were subtracted from the ΔCt values for LE fish. The values of the resulting ΔΔCt were then taken as the potentiator of 2 (2(ΔΔCt)), resulting in average fold change of emergence rank per time point.

Data analysis was performed using standard procedures [54]: A LOWESS curve (5-point window) was fitted over the data plotted over time and a sine wave was fitted subsequently over the resulting curve. The area under the curve (AUC) and the amplitude of the sine wave were further plotted against emergence rank (GraphPad Prism 6.0), and the correlation was analysed by a Spearman rank test (N = 72, significance accepted at p < 0.05).

Determination of cortisol and melatonin levels

Whole body cortisol and melatonin concentrations were measured in the bodies of the fish used for the qPCR analysis of the brains. After weighing and pulverising them in liquid nitrogen, 500 μl ice-cold 1 × PBS buffer was added and the samples were vortexed for 1 min. Subsequently, 500 μl chloroform was added and vortexed for 1 min and then centrifuged at 1500g for 5 min. Following centrifugation, the organic layer of each sample containing melatonin and cortisol was transferred to a separate test tube. Extraction was repeated three times. Samples were kept overnight in the fume hood for evaporation of chloroform. The next day, 200 μl ice-cold 1 × PBS buffer was added, and 100 μl was used for the determination of melatonin, and 50 μl for cortisol, using ELISA kits (melatonin: Abelissa, Alachua, FL, 32615 USA; cortisol: Demetic Diagnostics GmbH, Kiel, Germany). Two-way ANOVA, with concentrations and time points as variables, indicated overall effects, and Sidaks’s multiple comparison test was used for the post hoc determination of differences between values per time point. Significance was accepted at p < 0.05.

A LOWESS curve (5-point window) was fitted over the data plotted over time, and a sine wave was fitted subsequently over the resulting curve. The AUC and the amplitude of the sine wave were further plotted against emergence rank (GrapPad Prism 6.0), and the correlation was analysed by a Spearman rank test (N = 72, significance accepted at p < 0.05).

Activity patterns over diurnal/circadian cycle

In order to determine activity levels over diurnal and circadian cycles, individuals were automatically tracked over a period of 3 days under the usual light-dark regime (14:10 LD). Three rows of eight Plexiglas tanks (33 × 13 × 13 cm, capacity 3 l) were placed on a table surface (1.0 × 1.0 m). A flow through system (47.6 l h−1) provided each tank with fresh system water. Underneath, a retro-reflecting tape surface (1.2 × 1.2 m, Scotchlite 3 M, St. Paul, USA) reflected infrared light from three spots (U48R Univivi, Shenzhen, China), placed around an infrared firewire camera (Dragonfly, Point Grey Research Inc., Richmond, Canada) in a distance of 7 cm. The camera light unit was placed 1.7 m above the table.

Fish previously tested in the group emergence test in sets of ten, with random elimination of two, were placed individually in the tanks. Each row contained fish from the same group emergence test, i.e. 8 fish in three rows, resulting in 24 fish per experimental run, with three experimental runs (N = 3), resulting in n = 72. Food was provided ad libitum in the form of negatively buoyant pallets (Weekend, Tetra, USA). Fish were acclimated 30–48 h prior to recording. During the exposure to diurnal cycles (L:D), filming started at the middle of the light phase (7:00 hCT) at day 0 and lasted until 7:00 hCT of day 3 (L:D). The light was then kept on for another 3 days (72 h) in order to create constant light conditions (L:L), and filming continued until 7:00 hCT of day 3 L:L. Locomotion patterns were quantified as average velocity (V, mm s−1) over periods of 30 min, using EthoVision XT 6 (Noldus Information Technology b.v., Wageningen, The Netherlands).

For the analysis of the locomotion activity, a LOWESS curve (5-point window) was fitted over the activity data plotted over time and subsequently fitted sine wave yielded area under the curve (AUC, dimensionless) as an indicator for general locomotion activity and amplitude (mm s−1) as an indicator for activity fluctuation.

The rhythm strength of the 24-h component in activity data was gained as follows: after testing for multicollinearity in the explanatory variables by using variance inflation factors, a multivariate model predicting emergence rank was made, based on the remaining variables. Cook’s distance was calculated for all data points, but no outliers were found. Using likelihood ratio tests and Akaike’s information criterion, the optimal model was determined. The only explanatory variable of this optimal model was rhythm strength of the 24-h component. Visual inspection of the residuals revealed no violation by this model of statistical assumptions.

Finally, the maximum velocity (Vmax, mm s−1), i.e. the average of the maximum values for V during the first three 24-h periods (L:D), gave an indication for absolute swimming capacity. AUC, amplitude, rhythm strength and Vmax were plotted against emergence rank (GraphPad Prism 6.0), and the correlation was analysed by a Spearman rank test (N = 72, significance accepted at p < 0.05).



We would like to thank Peter Snelderwaard for his invaluable help with the setup and animal care.


No funding source to report

Availability of data and materials

All data generated during this study are included in this published article and its additional information files, and all materials generated during this study are available upon request.

Authors’ contributions

CT designed the study, conducted the experiments, analysed the results and wrote the paper. HS designed the study and wrote the paper. YR analysed the results and wrote the paper. EH conducted the experiments. JM wrote the paper. HS analysed the results and wrote the paper. MS designed the study and wrote the paper. All authors read and approved the final manuscript

Ethics approval and consent to participate

Zebrafish were maintained and handled according to the guidelines from the Zebrafish Model Organism Database (ZFIN, and in compliance with the directives of the local animal welfare committee of Leiden University (DEC number 14058).

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

Institute of Biology, Leiden University, Leiden, The Netherlands
Molecular Cell Biology, Leiden University Medical Center, Leiden, The Netherlands


  1. McAdams DP, Olson BD. Personality development: continuity and change over the life course. Ann Rev Psychol. 2010;61:517–42.View ArticleGoogle Scholar
  2. Bardo MT, Donohew RL, Harrington NG. Psychobiology of novelty seeking and drug seeking behavior. Behav Brain Res. 1996;77:23–43.View ArticleGoogle Scholar
  3. Caci H, Robert P, Boyer P. Novelty seekers and impulsive subjects are low in morningness. Europ Psych. 2004;19:79–84.View ArticleGoogle Scholar
  4. Maestripieri D. Night owl women are similar to men in their relationship orientation, risk-taking propensities, and cortisol levels: implications for the adaptive significance and evolution of eveningness. Evol Psychol. 2014;12:130–47.View ArticleGoogle Scholar
  5. Maestripieri DM, Baran N, Sapienza P, Zingales L. Between- and within-sex variation in hormonal responses to psychological stress in a large sample of college students. Stress. 2010;13:413–24.View ArticleGoogle Scholar
  6. Ursin H. The psychology in psychoneuroendocrinology. Psychoneuroendocrinol. 1998;23;6:555–70.Google Scholar
  7. Ebstein RP. The molecular genetic architecture of human personality: beyond self-report questionnaires. Mol Psychiatry. 2006;11:427–45.View ArticleGoogle Scholar
  8. Benjamin J, Li L, Patterson C, Greenberg BD, Murphy DL, Hamer DH. Population and familial association between the D4 dopamine receptor gene and measures of novelty seeking. Nature Genet. 1996;12:81–4.View ArticleGoogle Scholar
  9. Rey S, Boltana S, Vargas R, Roher N, MacKenzie S. Combining animal personalities with transcriptomics resolves individual variation within a wild-type zebrafish population and identifies underpinning molecular differences in brain function. Mol Ecol. 2013;22:6100–15.View ArticleGoogle Scholar
  10. Huntingford FA, Andrew G, MacKenzie S, Morera D, Coyle SM, Pilarczyk M, et al. Coping strategies in a strongly schooling fish, the common carp Cyprinus carpio. J Fish Biol. 2010;76:1576–91.View ArticleGoogle Scholar
  11. Tudorache C, Schaaf M, Slabbekoorn H. Covariation between behaviour and physiology indicators of coping style in zebrafish (Danio rerio). J Endocrinol. 2013;219:251–8.View ArticleGoogle Scholar
  12. Koolhaas JM, Korte SM, de Boer SF, Van Der Vegt BJ, Van Reenen CG, Hopster H, et al. Coping styles in animals: current status in behavior and stress-physiology. Neurosci Biobeh Rev. 1999;23:925–35.View ArticleGoogle Scholar
  13. Wolf M, van Doorn GS, Leimar O, Weissing FJ. Life-history trade-offs favour the evolution of animal personalities. Nature. 2007;447:581–4.View ArticleGoogle Scholar
  14. Conrad JL, Weinersmith KL, Brodin T, Saltz JB, Sih A. Behavioural syndromes in fishes: a review with implications for ecology and fisheries management. J Fish Biol. 2011;78:395–435.View ArticleGoogle Scholar
  15. Øverli Ø, Sørensen C, Pulman KGT, Pottinger TG, Korzan WJ, Summers CH, et al. Evolutionary background for stress-coping styles: relationships between physiological, behavioral, and cognitive traits in non-mammalian vertebrates. Neurosci BiobehavRev. 2007;31:396–412.View ArticleGoogle Scholar
  16. Weiss A. Personality traits: a view from the animal kingdom. J Pers. 2018;86:12–22.View ArticleGoogle Scholar
  17. Friedman HS. The multiple linkages of personality and disease. Brain Beh Imm. 2008;22:668–75.View ArticleGoogle Scholar
  18. van der Greef J, Hankemeier T, McBurney RN. Metabolomics-based systems biology and personalized medicine: moving towards n = 1 clinical trials? Pharmacogenomics. 2006;7:1087–94.View ArticleGoogle Scholar
  19. Schork NJ. Personalized medicine: time for one-person trials. Nature. 2015;520:609–11.View ArticleGoogle Scholar
  20. Kumar D. The personalised medicine. A paradigm of evidence-based medicine. Ann Ist Super Sanita. 2011;47:31–40.PubMedGoogle Scholar
  21. Kaushik A, Vasudev A, Arya SK, Pasha SK, Bhansali S. Recent advances in cortisol sensing technologies for point-of-care application. Biosens Bioelectron. 2014;53:499–512.View ArticleGoogle Scholar
  22. Tudorache C, Braake ter A, Tromp M, Slabbekoorn H, Schaaf M. Behavioural and physiological indicators of stress coping styles in larval zebrafish. Stress. 2014;18:121–8.View ArticleGoogle Scholar
  23. Roenneberg T, Merrow M. The circadian clock and human health. Curr Biol. 2016;26:R432–43.View ArticleGoogle Scholar
  24. McKenna HT, Reiss IK, Martin DS. The significance of circadian rhythms and dysrhythmias in critical illness. J Intens Care Soc. 2017;18:121–9.View ArticleGoogle Scholar
  25. Bollinger T, Schibler U. Circadian rhythms – from genes to physiology and disease. Swiss Med Wkly. 2014;144:1–11.Google Scholar
  26. Hardin PE. From biological clock to biological rhythms. Genome Biol. 2000;1(4):1023.1–5.View ArticleGoogle Scholar
  27. Foulkes NS, Whitmore D, Vallone D, Bertolucci C. Studying the evolution of the vertebrate circadian clock: the power of fish as comparative models. Adv Genet. 2016;95:1–30.View ArticleGoogle Scholar
  28. Adan A, Archer SN, Hidalgo MP, Di Milia L, Natale V, Randler C. Circadian typology: a comprehensive review. Chronobiol Int. 2012;29:1153–75.View ArticleGoogle Scholar
  29. Gorissen M, Manuel R, Pelgrim TNM, Mes W, de Wolf MJS, Zethof J, et al. Differences in inhibitory avoidance, cortisol and brain gene expression in TL and AB zebrafish. Genes Brain Behav. 2015;14:428–38.View ArticleGoogle Scholar
  30. Guryev V. Genetic variation in the zebrafish. Genome Res. 2006;16:491–7.View ArticleGoogle Scholar
  31. Brown KH, Dobrinski KP, Lee AS, Gokcumen O, Mills RE, Shi X, et al. Extensive genetic diversity and substructuring among zebrafish strains revealed through copy number variant analysis. Proc Natl Acad Sci. 2012;109:529–34.View ArticleGoogle Scholar
  32. Stickney HL. Rapid mapping of zebrafish mutations with SNPs and oligonucleotide microarrays. Genome Res. 2002;12:1929–34.View ArticleGoogle Scholar
  33. Wong RY, Oxendine SE, Godwin J. Behavioral and neurogenomic transcriptome changes in wild-derived zebrafish with fluoxetine treatment. BMC Genomics. 2013;14:348.View ArticleGoogle Scholar
  34. Drew RE, Settles ML, Churchill EJ, Williams SM, Balli S, Robison BD. Brain transcriptome variation among behaviorally distinct strains of zebrafish (Danio rerio). BMC Genomics. 2012;13:323.View ArticleGoogle Scholar
  35. Weger BD, Sahinbas M, Otto GW, Mracek P, Armant O, Dolle D, et al. The light responsive transcriptome of the zebrafish: function and regulation. Mueller F, editor. PLoS One 2011; 6: e17080–e17015.Google Scholar
  36. Bell AM, Sih A. Exposure to predation generates personality in threespined sticklebacks (Gasterosteus aculeatus). Ecol Lett. 2007;10:828–34.View ArticleGoogle Scholar
  37. Sih A, Bell A, Johnson JC. Behavioral syndromes: an ecological and evolutionary overview. Trends Ecol Evol. 2004;19:372–8.View ArticleGoogle Scholar
  38. Bell AM. Evolutionary biology: animal personalities. Nature. 2007;447:539–40.View ArticleGoogle Scholar
  39. Weinert D, Schottner K. An inbred lineage of Djungarian hamsters with a strongly attenuated ability to synchronize. Chronobiol Int. 2007;24:1065–79.View ArticleGoogle Scholar
  40. Smith JM. Evolution and the theory of games. Cambridge: Cambridge University Press; 2009.Google Scholar
  41. Griefahn B, Bröde P, Remer T, Blaszkewicz M. Excretion of 6-hydroxymelatonin sulfate (6-OHMS) in siblings during childhood and adolescence. Neuroendocrinology. 2003;78:241–3.View ArticleGoogle Scholar
  42. Waldhauser F, Weiszenbacher G, Tatzer E, Gisinger B, Waldhauser M, Schemper M, et al. Alterations in nocturnal serum melatonin levels in humans with growth and aging. J Clin Endocrin Metabol. 1988;66:648–52.View ArticleGoogle Scholar
  43. Zeitzer JM, Daniels JE, Duffy JF, Klerman EB, Shanahan TL, Dijk D-J, et al. Do plasma melatonin concentrations decline with age? Am J Med. 1999;107:432–6.View ArticleGoogle Scholar
  44. Srinivasan V, Pandi-Perumal SR, Trakht I, Spence DW, Hardeland R, Poeggeler B, et al. Pathophysiology of depression: role of sleep and the melatonergic system. PsychiatrRes. 2009;165:201–14.Google Scholar
  45. Michel S, Block GD, Meijer JH. The aging clock. In: Colwell CS, editor. Circadian medicine. Hoboken: Wiley, Blackwell; 2015. p. 321–35.Google Scholar
  46. Roff DA. The evolution of genetic correlations: an analysis of patterns. Evolution. 1996;50:1392.View ArticleGoogle Scholar
  47. Solovieff N, Cotsapas C, Lee PH, Purcell SM, Smoller JW. Pleiotropy in complex traits: challenges and strategies. Nature Rev Gen. 2013;14:483–95.View ArticleGoogle Scholar
  48. Ben-Moshe Livne Z, Alon S, Vallone D, Bayleyen Y, Tovin A, Shainer I, et al. Genetically blocking the zebrafish pineal clock affects circadian behavior. PLoS Gen. 2016;12:e1006445–26.View ArticleGoogle Scholar
  49. Dickmeis T, Weger BD, Weger M. The circadian clock and glucocorticoids--interactions across many time scales. Mol Cell Endocrinol. 2013;380:2–15.View ArticleGoogle Scholar
  50. de Kloet RE, Joëls M, Holsboer F. Stress and the brain: from adaptation to disease. Nat Rev Neurosci. 2005;6:463–75.View ArticleGoogle Scholar
  51. Gosling SD, John O. Personality dimensions in nonhuman animals: a cross-species review. Curr Dir Psychol Sci. 2016;8:69–75.Google Scholar
  52. van der Vinne V, Riede SJ, Gorter JA, Eijer WG, Sellix MT, Menaker M, et al. Cold and hunger induce diurnality in a nocturnal mammal. Proc Natl Acad Sci. 2014;111:15256–60.View ArticleGoogle Scholar
  53. van Oosterhout F, Lucassen EA, Houben T, HT vL, Antle MC, Meijer JH. Amplitude of the SCN clock enhanced by the behavioral activity rhythm. PLoS One. 2012;7:e39693.View ArticleGoogle Scholar
  54. Zuur AF, Ieno EN, Walker N, Saveliev AA, Smith GM. Mixed effects models and extensions in ecology with R. New York: Springer New York; 2009.View ArticleGoogle Scholar


© The Author(s). 2018