1 General Aspects of Biological Timing
In order to adapt to a changing environment, most living organisms have established the ability to measure time by the development of biological timing systems. There are different time-scales in which these temporal systems operate, thus covering more than 12 orders of magnitude, from microseconds to days or even seasons. All of these biological ‘clocks’ use time-keeping mechanisms that provide organisms with temporal information crucial to fulfill the required biological processes (Buhusi & Meck, 2005; Golombek, Bussi, & Agostino, 2014). Most chapters in this book are focused on temporal processing in the millisecond-to-minutes range, which is crucial to many forms of cognitive processing. In the present chapter, we will cover different aspects of biological timing in the circadian (i.e., ‘about a day’, 24-h) range. Since circadian timing is ubiquitous and regulates many physiological and behavioral functions – including cognitive functions – the study of circadian rhythms has a deep impact in neuroscience, medicine and health research.
Despite early classical experiments such as De Mairan’s 1729 demonstration of an endogenous rhythm in leaf movement of Mimosa plants, among many others, the notion of an inner biological clock remained quite elusive until the mid-twentieth century (Golombek & Rosenstein, 2010). Several lines of evidence revealed the presence of an endogenous daily biological clock, and the term circadian was introduced in the fifties to identify self-sustained rhythms under constant conditions.
At present, it is very well known that the circadian clock is an endogenous timing system that organizes the environmental oscillations every 24-h. The periodic changes in the light/dark cycle produced by the rotation of the Earth influenced the development of circadian rhythms in almost all living organisms. Hence, organisms are able to adapt to and anticipate the changes in lighting conditions, thus settling their physiological processes to specific times of the day. Some physiological and behavioral functions regulated by the circadian system include the sleep-wake cycle, body temperature, hormone release, and gene expression. For example, nearly half of all mammalian genes are rhythmically expressed in one or more tissues (Yan, Wang, Liu, & Shao, 2008).
The circadian system have three fundamental characteristics: (1) the presence of an endogenous oscillator, which generates rhythms with a period close to 24-h in the absence of external cues; (2) the capacity to be synchronized by rhythmic environmental stimuli; and (3) temperature compensation, which means that the biochemical reactions that take place in the circadian clock do not change their speed with temperature variations (Buhr & Takahashi, 2013). Additionally, the circadian system is composed by three main components that allow the regulation of physiology and behavior: (1) an input pathway sensing the environment, (2) an endogenous core oscillator, and (3) an output pathway that couples the core oscillator with the rest of the organism. In mammals, the main circadian oscillator (called the core or master oscillator) is located in the suprachiasmatic nuclei (scn) of the ventral hypothalamus. The scn synchronizes subordinate organ and tissue clocks (called peripheral circadian oscillators) using mainly autonomic pathways from the hypothalamus, driving neurochemical, endocrine, and metabolic signaling pathways (Albrecht, 2012).
2 Molecular Components of the Circadian System
The scn neurons have the ability to generate autonomous circadian rhythms (Welsh, Logothetis, Meister, & Reppert, 1995). The intrinsic rhythmicity of each neuron is based on a molecular oscillation that involves negative feedback loops of transcriptional and translational processes (Figure 1.1). In the mammalian circadian system, a primary feedback loop is composed by the positive elements CLOCK and BMAL1, which heterodimerize (i.e., join themselves to form a complex of non-identical monomers) and promote the transcription of Per (Period) and Cry (Cryptochrome) genes (negative elements) by acting on specific biding sites (E-box sequences) located in their promoter region (Gekakis et al., 1998; Munoz, Brewer, & Baler, 2006; Yoo et al., 2005). After transcription and translation, PER and CRY proteins heterodimerize and translocate to the nucleus to inhibit their own transcription by acting on the dimer CLOCK/BMAL1 (Ye et al., 2014).
Each transcriptional/translational loop takes near to 24-h to complete. In addition to this primary loop, another negative feedback loop contributes to clock precision and robustness. In the second loop, the protein REV-ERBα moves to the nucleus to repress Bmal1 expression; inversely, Rorα can stimulate its transcription, both binding to rore (Retinoic acid-related Orphan receptor Response Element) sites on the Bmal1 promoter (Guillaumond, Dardente, Giguere, & Cermakian, 2005; Preitner et al., 2002; Sato et al., 2004). These oscillations of negative and positive elements generate circadian rhythmicity,
3 Synchronization to the Environment
The circadian system in nature is normally exposed to a rhythmic 24-h environment, so that appropriate signals (called Zeitgebers, from German Zeit, ‘time’; geben, ‘to give’), such as light, temperature, or food, synchronize its oscillation (Golombek & Rosenstein, 2010). This effect in Chronobiology is called entrainment. Therefore, under entrained conditions, circadian rhythms adjust their endogenous period (called τ) to that of the Zeitgeber (called T). It is worth to mention that the concept of entrainment – the phenomenon in which two or more independent rhythmic processes synchronize with each other – has been described in a wide variety of systems and over a wide range of timescales: some examples include fireflies illuminating in synchrony, human individuals adjusting their speech rhythms to match each other in conversation, foot tapping to the beat of a song, and attention to rhythmic patterns in the milliseconds scale (Jones, 1976; Clayton, Sager, & Will, 2005). Circadian rhythms of living organisms are entrained to rhythmic environmental cues, such as the 24-h light/dark alternation or daily temperature cycles. Locomotor activity rhythms are one of the most studied circadian outputs, and are frequently represented as actograms (Figure 1.2), where horizontal lines represent consecutive days and black bars describe the locomotor activity episode. Each horizontal line in Figure 1.2 shows one experimental day along the vertical axis, and the hours of the day are represented in the x axis. Further information about actograms is provided in Section 6.3.
In mammals, the most powerful synchronizer is the daily light/dark cycle. Under normal 24-h light/dark conditions, circadian rhythms are entrained with a period equal to 24-h (Figure 1.2A).
Light stimulates a group of photosensitive retinal ganglion cells (pRGC) that express the photopigment melanopsin (Panda et al., 2002) and project to the scn through the retinohypothalamic tract (Morin, 2013). The effect of light on scn neurons leads to an increase in the intracellular concentrations of Ca2+, which initiates a signal transduction cascade that ultimately results in a phase shift of the circadian clock (Golombek, Agostino, Plano, & Ferreyra, 2004; Golombek et al., 2003; Golombek & Rosenstein, 2010; Morin & Allen, 2006). The circadian clock differs in its temporal responsiveness to external stimuli. For example, in nocturnal rodents, exposure to light synchronizes circadian rhythms by inducing phase delays during the early night and phase advances
The effect of light in synchronizing the circadian clock depends on the time of day, thus defining a phase-response curve (prc) that indicates the times at which light can induce phase delays, phase advances, or no change at all (Figure 1.2B, Dunlap et al., 2004; Golombek & Rosenstein, 2010; Johnson, Elliott, & Foster, 2003). The prc is one of the main tools to study sensitivity of the circadian system to either light (photic entrainment) or many other stimuli potentially capable of resetting the clock (called non-photic entrainment). The prc, its shape and amplitude, is an intrinsic property of the circadian oscillator. Indeed, the prc is characteristic for each species. In humans, it is also the basis for the design of treatments in abnormal entrainment situations, such as jet-lag, shiftwork or circadian-related sleep disruption (Gooley, 2008; Lewy, Emens, Jackman, & Yuhas, 2006).
Other entrainment cues like food intake, environmental changes in temperature, drugs or even social interactions can reset circadian rhythms, eliciting phase shifts mainly during the subjective day (the resting period in rodents; Figure 1.2B).
When photic or non-photic Zeitgebers directly affect a rhythmic output without affecting the circadian oscillator, for instance directly suppressing or stimulating locomotor activity, this mechanism is called masking (Mrosovsky, 1999). Different from entrainment, masking influences a circadian output (such as behavior) without the participation of the endogenous circadian clock.
3.1 Zeitgeber Time and Circadian Time
Since a strong Zeitgeber defines the rhythm of the circadian oscillator, under entrained conditions time is expressed as Zeitgeber time (zt). In the case of photic entrainment, within a light/dark schedule of 12 hours of light and 12 hours of darkness (ld 12:12), ZT0 is defined as the time of ‘lights ON’ (i.e., the beginning of the light phase). Similarly, ZT12 corresponds to the time of ‘lights OFF’, the end of the light phase. Analogous terminology is applied for non-photic Zeitgebers, such as food or temperature.
When organisms are kept under constant conditions by avoiding contact with external time cues, they display so-called freerunning or circadian rhythms, no longer equal to 24-h. Therefore, time cannot be expressed in zt but it is expressed in circadian time (ct) units. One circadian cycle is divided
3.2 In the Limit of Synchronization: T-cycles, Forced Desynchronization and Relative Coordination
Under natural conditions of entrainment to a Zeitgeber of T=24-h, variations of the subject’s endogenous period (τ) may lead to different phases of entrainment, that are associated with chronotypes. In this sense, short endogenous periods τ often lead to early phases (‘morning larks’) and long periods τ correspond to late phases (‘night owls’). Mutations that affect the intrinsic period τ result in large shifts of the entrainment phase. For example, patients with familial advanced sleep phase syndrome have a short endogenous period of τ = 23.3-h leading to a phase advance of more than 3-h (Jones et al., 1999). Thus, the phase of entrainment is a function of the period mismatch τ – T (Bordyugov et al., 2015).
Most organisms can only entrain to Zeitgebers cycling with a period close to 24-h. If the entraining period is too short or too long – thus, exceeding the range of entrainment – the circadian system cannot follow the Zeitgeber anymore. In the laboratory, circadian rhythms can be entrained to variable Zeitgebers with different T periods (called T-cycles) in order to determine the range of entrainment for each species (Aschoff, Daan, & Honma, 1982; Jud, Schmutz, Hampp, Oster, & Albrecht, 2005). To determine this range of entrainment, T-cycles of different periods T≠24-h are applied. For example, 22-h or 26-h T-cycles may be used to define the limits of stable entrainment. These limits differ from species to species and depend on the nature of the applied Zeitgeber. Notably, under T periods close to the limits of circadian entrainment, nocturnal rodents present two patterns of behavior, called forced desynchronization and relative coordination (Campuzano, Vilaplana, Cambras, & Diez-Noguera, 1998; Casiraghi, Oda, Chiesa, Friesen, & Golombek, 2012; de la Iglesia, Cambras, Schwartz, & Diez-Noguera, 2004; Golombek et al., 2013). Forced desynchronization (Figure 1.3A) occurs when a system of dual, weakly coupled circadian oscillators is subjected to a Zeitgeber that is outside the range of entrainment of one of these oscillators. For example, it has been described that when rats were maintained under light/dark cycles with period of 22- and 23-h, two circadian rhythms were simultaneously present in their motor activity, temperature, and feeding activity: one rhythm followed the external light/dark cycle, whereas the other rhythm was under freerunning with a period of around 25-h (Madrid, Lax, Vilaplana, Cambras, & Díez-Noguera, 1992;
3.3 Circadian Desynchronization
When the timing of the light/dark cycle is altered – for example, as a result of transmeridian travel or work schedules – it leads to a misalignment of the circadian system with the environment.
A well-known consequence of circadian alteration is jet-lag disorder, which results from rapidly crossing times zones – with the severity of symptoms typically depending on the number of times zones crossed, the travel frequency, and the direction of travel (east or west). Symptoms vary from insomnia and excessive sleepiness to more severe complications, including increased cancer risk, metabolic disorder, cognitive deficits, reduced fertility, and deregulated immune responses (Castanon-Cervantes et al., 2010; Golombek et al., 2013; Karatsoreos, Bhagat, Bloss, Morrison, & McEwen, 2011; Mahoney, 2010; Scheer, Hilton, Mantzoros, & Shea, 2009). In animal models, interference with the cGMP-related photic entrainment pathway has been proven to be a useful tool for accelerating synchronization to advances of the light/dark cycle, which may reduce jet-lag symptoms due to eastbound flights (Agostino et al., 2007; Plano, Agostino, de la Iglesia, & Golombek, 2012).
Another circadian alteration is produced by shift work disorder, which affects people who frequently rotate shifts or work at night. Under these conditions, the alternation between a normal sleep/wake cycle and being awake at night and sleeping during the day not only alters the light/dark schedule in relation to the endogenous circadian clock, but also modifies feeding patterns, social habits, etc. (Reid & Abbott, 2015; Zee & Goldstein, 2010). Sleep loss, in addition to circadian misalignment, contributes to decreased alertness during night work. Current treatment of shift work disorder includes strategies to achieve and maintain some degree of circadian alignment (e.g., photic phase shift protocols, avoidance of light during the day, melatonin administration, etc.), improve sleep (using hypnotics, melatonin, and behavioral
Rodent models based on disruption of circadian rhythms provide useful information to study desynchronization due to modifications in the environment, such as abrupt changes in the light/dark cycle or nocturnal exposure to light. A commonly used model of circadian disruption is constant light (ll) exposure (Figure 1.3C), which induces period lengthening followed by
4 Circadian Influence on Cognition
As mentioned, the circadian system regulates several physiological and behavioral functions, and disruption of the circadian clock may have severe consequences for physical and mental health. In this section, we will overview an important example of circadian modulation: the influence of the circadian system on cognitive functions.
There is evidence that cognitive performance and learning is influenced by circadian oscillators (Eckel-Mahan & Storm, 2009; Gerstner & Yin, 2010). Also, regularly-timed cognitive processes impact circadian rhythms (Gritton et al., 2012), indicating a bidirectional interaction between cognitive performance and circadian processes. Furthermore, at least for some tasks in nocturnal rodents, night-phase performance can never be equaled by light-phase performance regardless of the strength of entrainment to the schedule, which may have deep implications for shiftwork therapies (Gritton et al., 2012). In addition, desynchrony between internal and environmental time has been associated with impaired cognitive function in animals (Devan et al., 2001; Gibson, Wang, Tjho, Khattar, & Kriegsfeld, 2010; Loh et al., 2010) and humans (Cho, 2001; Cho, Ennaceur, Cole, & Suh, 2000; Folkard, 1996; Santhi, Horowitz, Duffy, & Czeisler, 2007). Although it is well established that cognitive abilities vary as a function of daytime, there is still a widespread view that the circadian clock has a singular role in cognition related to sleep timing, and that the main factor for cognitive maintenance is the quality and duration of sleep. However, there is also evidence that indicates a circadian control of cognition beyond sleep timing, as revealed by forced desynchronization protocols (Kyriacou & Hastings, 2010; Wright, Hull, & Czeisler, 2002).
A fundamental component of cognition is the perception of the passage of time. In particular, temporal processing within the hundreds of milliseconds and the seconds-to-minutes range – known as interval timing – is crucial for many complex behaviors, such as speech comprehension, working memory,
5 Different Approaches to Measuring Circadian Rhythms
Several outputs are used to assess the overall status of the circadian system. These variables (called circadian marker rhythms) can be easily measurable over long periods of time, preferably using non-invasive methods. The most widely used circadian marker rhythms in several organisms are those of activity, temperature, hormones, and clock gene expression (LeSauter & Silver, 1998). In humans, the main variables used are actigraphy, core body temperature (cbt), and plasma or salivary melatonin (Bonmati-Carrion et al., 2014; Klerman, Gershengorn, Duffy, & Kronauer, 2002). In nocturnal rodents, wheel-running activity is the most commonly used and reliable method for assessing the output of the circadian clock. Other measures that can be used include general activity, food-bin approaches, drinking, sleep, and body temperature.
Bioluminiscence has become a powerful tool for studying circadian rhythms. Exogenous luciferase genes, under the control of promoters conferring circadian regulation, have been introduced into a wide variety of organisms, including cyanobacteria (Synechococcus), plants (Arabidopsis), insects (Drosophila), nematodes (C. elegans) and rodents (Welsh & Kay, 2005; Goya, Romanowski, Caldart, Bénard, & Golombek, 2016). These transgenic organisms exhibit robust circadian rhythms of luminescence, useful as a longitudinal assay of rhythmicity for a wide range of recent genetic and biochemical studies of circadian clock mechanisms.
6 Circadian Data Analysis
6.1 Period, Amplitude and Phase
Like other periodic patterns, circadian rhythms can be described using mainly three parameters: a) period (τ), defined as the time in which a cycle is completed, b) amplitude, defined as the difference between the peak (or trough) and the mean value of the variable under study, and c) phase or phase angle (φ), the timing of a reference point in the cycle relative to a fixed event (Figure 1.4). There are several conventions to define the phase of a rhythm; however, the most commonly used is when the amplitude reaches its maximum, which is called acrophase in the Cosinor analysis (see below; Dunlap et al., 2004). Consequently, every time that we mention the phase of a rhythmic signal, it will be referring to its acrophase. Indeed, the algorithms that will be described in the following sections will also compute the acrophase (time of peak) of the signal.
These parameters (period, amplitude, and phase) can describe any periodic pattern, but are normally used to characterize a sinusoidal wave. In many cases, however, this is not the real situation (i.e., locomotor activity may behave as a square-like rhythm). For this reason, most of the algorithms that are usually used to compute period and phase have been tested for different waveforms and are still valid for data analysis. Even though in this chapter we will only analyze the mentioned parameters, it is worth noting that according to the type of experimental data there are many other parameters that can be useful to analyze (e.g., sinusoidal dampening, fragmentation, coupling of the activity and Zeitgeiber, phase angle, etc.).
6.2 Sampling and Signal Conditioning
Before a complete overview of the methods used for analyzing circadian rhythms, it is important to refresh some basic concepts of signal analysis, such as sampling and signal conditioning (filtering and detrending). This information will allow us to know beforehand some of the limitations of our analytic methods.
A fundamental aspect of data analysis in temporal series is sampling. The first step is to decide between uniform sampling and non-uniform sampling (i.e., when data cannot be collected uniformly spaced). When possible, a uniform sample is clearly the best option, because it allows a wider and stronger set of statistical tools. Even if some data points were lost, the missing data could still be interpolated in a reliable and easier way. Another key feature of the experimental design is to define the sampling rate – the time interval between two consecutive samples, or bin – according to the needs and objectives of the undergoing experiment. In this age of fast computers and large hard drives, the trend is to choose a very high frequency of sampling and disregard the storage size; however, sample size does matter and it is a key factor for further data analysis. The sampling rate will define the shortest cycle that can be measured. According to the Nyquist theorem, the sample rate must be at least twice the higher frequency (i.e., 1/ τ) that it is aimed to study (Levine, Funes, Dowse, & Hall, 2002). For example, if we aim only to find the circadian frequency (24-h), we could just sample every 12 hours. However, in this case we cannot compute a phase or an accurate period. On the other hand, if the sampling rate were once per hour, then it would be possible to evaluate periodicities down to two hours (but not shorter), since a minimum of two points is required to describe a cycle. For this reason, the precision of the experimental approach is limited by the sampling procedure. However, a very high sampling rate would implicate low sample size – defining sample size as the length of the integrated time of the measure. In this case, when a very low sample size is used, the time for data processing and visualization greatly increases. According to the experimental output (e.g., locomotor activity, gene or protein oscillation, etc.), a compromise is reached between accuracy and speed. The usual sample size for circadian analysis of locomotor activity data varies between 5 to 30 minutes. This sample size will allow an accurate determination of circadian period and phase, as well as the study of the transient behavior of these parameters (including period and phase shift).
Sometimes there are some elements in the signal that may interfere with the calculation of circadian parameters. In this case, the raw data often need to be filtered for further study; this processing of the raw signal before the analysis is known as signal conditioning. Two of these interfering elements will be addressed in this section: the presence of a shifting temporal baseline (i.e., linear trend) and the presence of a high frequency noise.
In the first case (the presence of a shifting temporal baseline), the two types of linear effects that can alter the subsequent data analyses involve:
In the second case (the presence of high frequency noise), the approach is to use digital filters. These filters diminish the spectral contribution of a certain wavelength to the signal while keeping the others intact. The filters that reduce the high frequency components are called low-pass or smooth. We will mention a simple filter, moving-average smoothing, but there are other types of filters more complex, with different advantages and disadvantages. A filter can be defined as a function that transforms a time series into another. The moving-average filter is one of the most simple and used in the circadian field. It is used to analyze a set of data by creating a series of averages of different subsets of the full data set. Specifically, this filter averages a fixed number of points – it is better to have an odd number of points – in order to find the mean value of the initial subset of data. Then the subset is modified by moving one point forward in time, and the process is repeated until we reach the end of the data series. The resulting plot line connecting all the averages is the moving average: a set of numbers, each of which is the average of the corresponding subset of a larger set of data. In this case we give each point the same weight, but sometimes it is recommended weighting differently each point mainly with a Gaussian distribution (Levine et al., 2002) to give the central point of the window a higher weight. An example of moving average can be found in the online Appendix section (smoothma.m file; see book’s GitHub repository).
As previously shown in Figures 1.2 and 1.3, actograms are a useful tool to represent circadian locomotor activity rhythms in several species. The actogram is usually performed as a ‘double plot’ graph, where two cycles are represented in the same horizontal line. Specifically, double-plotted actograms show the second cycle on the right side of each line, as well as at the start of the following horizontal line, and so on. If the actogram is usually set at a fixed period of 24-h (for a solar day) in the x-axis, then the double plot representation will have 48-h duration (Figure 1.5, left panel). The cycle is divided into a time fraction with main activity episodes (called alpha) and other with major rest episodes (called rho). Sometimes, scattered activity can be observed in the rho time because the animal interrupts its sleep for a short time (Jud et al., 2005). Double-plotting is especially useful to visualize non-24-h rhythms or very long series.
One of the reasons that the actogram became a key representation in Chronobiology is that it allows a fast and reliable visual analysis. With a simple observation we can realize the changes of period and phase. The actogram was first introduced by Colin Pittendrigh at the beginning of the 1960s (Pittendrigh, 1965), well before any computational method for circadian data analysis was established. By that time, data were plot and measured by hand. Currently, there are many computational tools (e.g., a circadian plugin for Imagej and the software El Temps) that still maintain the concept of visual analysis (‘paper and ruler’ analysis, in which the actogram is printed and analyzed with common measuring instruments).
The main downfall of this approach is the variability of the observation between users and the lack of rigorous quantitative analysis. Nevertheless, this kind of analysis is still heavily used. A quick example of this visual inspection in the actogram can be seen in Figure 1.6, where locomotor records with different periods, as well as phase shifts, can be easily identified.
An example of a double-plotted actogram containing real data can be found in the online Appendix, by running the function Actogram (see the files SCRIPT.m and ScriptFilter.m for further details; see book’s GitHub repository).
Periodic data of known period could be divided into sections, with the length of these sections matching the underlying period. Overlaying these sections (cycles) will produce a clear waveform, with peak and trough (Figure 1.5, right panel). In other words, a waveform can be thought as an average of all the lines of the actogram. The length of the waveform is given by the length of the cycle, which can be empirically fixed (e.g., the length of a light/dark cycle) or previously calculated by any method. Even though it can be used in many
Once the limits of alpha and rho phases are defined, their length (tα and tρ) is easily measured by taking into account that T = tα + tρ. The most usual way to represent these times is relative to the length of T (for alpha tα% = 100 tα / T). For example, for T = 24-h and tα = 12-h, then tα% = 50, meaning that half the total period corresponds to the activity phase. Along with the length of alpha and rho phases, it can be also calculated the area under the curve for both time lapses (Aα and Aρ) and to express them as percentages of the total area (Aα% and Aρ%). Finally, with both the length and amount of activity in both states, it can be computed the degree of agglomeration (Aα% / tα%). When this degree is near 1, the distribution of activity is uniform (i.e., we have a weak variation in the rhythm), while when this degree is higher than 1, most of the activity is concentrated in the alpha phase (Madrid & de Lama, 2006).
A key step in the analysis of circadian data is to accurately calculate the underlying period. There are several methods for determining period, which range from graphic approaches, such as the actogram – mentioned in Section 6.3 – to more mathematically complex approaches, such as Fourier-transform and non-Fourier-transform based algorithms (Zielinski, Moore, Troup, Halliday, & Millar, 2014). The latter methods allow a level of significance for period calculation, as well as an unbiased estimation.
A periodogram is a function relating periodic components of a time series to their spectral power. Many algorithms for period estimation had been described in the literature, all of them with a varying degree of complexity and different assumptions. Nevertheless the calculation method, all periodograms possess the same graphic representation, which makes them quite easy to understand (Figure 1.6). Several authors (Refinetti, 1993; Zielinski et al., 2014) compare the different methods for periodogram analysis (e.g., Fourier, Enright, Lomb-Scargle, etc.) and describe the best option for different cases. Each method has its own strength and weakness. In this section, we will mention two widely used methods: Lomb-Scargle (ls, very useful in short series and non-uniform sampling) and the Chi-square method, also known as Sokolove-Bushell (sb).
The most popular and conceptually simpler method is the Chi-square periodogram – its name corresponds to the statistical distribution that it uses. This method consists in splitting the data in sections of a fixed length. When this fixed length matches the period of the series, the overlay of all sections will give a clear waveform (with a strong variation) as the repeating patterns of the series coincide. However, if the length of sections does not coincide, then the
Another general approach is based in the curve-fitting criterion. This method assumes an a-priori periodic signal (mainly sinusoidal) and optimizes the signal parameters to best fit the data. In this case, we plot the correlation coefficient after the fit with different periods. This approach allows using shorter series and even non-uniform sampling. The Fourier method is the most basic approach of this set of algorithms. As we use a sinusoidal signal to fit, the method is especially sensitive to this waveform; however, it is not as good to fit other types of rhythmic signals. The Lomb-Scargle periodogram consists of a spectral analysis (a modification of the Fourier periodogram). This method creates a spectrum representing the significance of each frequency with a normalized power, and is usually the best choice in short series and non-uniform sampling. This method combines both cosine and sine fit to find periodicity regardless of the data phase. Another interesting feature is that it reduces the subsequent sub-harmonics of the main period (if the period is T, the peaks in 2T, 3T, etc.) and, therefore, the periodogram analysis and representation are ‘clean’.
There are many other methods to calculate the underlying period, such as mesa (Maximum Entropy Spectral Analysis) or fft-nlls (Fast Fourier Transform Non-linear Least Square algorithm). mesa is very useful because it can resolve multiple periods that are relatively close (Dowse, 2013). It is beyond the objective of the present chapter to describe all of them, but it is important to recognize which method is the best choice according to the specific experimental data. A comparison of the most used algorithms for determining period can be found in Zielinski et al. (2014).
An example containing the periodogram calculation of real data can be found in the online Appendix, by running the functions chisquare (Chi-square or Sokolov-Bushel periodogram) and lsssigx (a simplify application of the fastlomb.m that gives the complete spectrum of the Lomb-Scargle
As it was already shown in Figure 1.4, a simple rhythm can be described as a sinusoidal wave. By taking into account all the rhythm parameters, the signal can be mathematically expressed as:(1)
where M is the mean of the adjusted data (Midline Estimating Statistic Of Rhythm, mesor), φ is the acrophase, A is the amplitude (difference between φ and M) and t corresponds to time. Finally, ω = 2π/T is the angular frequency, being T the rhythm period. The mathematical procedure that fits the best curve – best fit according to the less-square method – is known as Cosinor. The expression of y leads to a nonlinear optimization to find the values of M, A and φ. The equation (1), by using the Euler’s equation through the formula cos(α – β) = cos(α)cos(β) + sin(α)sin(β), is transformed to:(2)
where a [A cos (φ)] and b [A sin (φ)] are known as the Euler coefficients, and are used to calculate the values of A and φ in the original equation (1) as follows:
If the signal was uniform sampled, the calculation is even simpler:(4)
The Cosinor method is mainly used to analyze a rhythmic time series with a known period (T), which can be obtained from the periodogram analysis (Cornelissen, 2014; Diez-Noguera, 2013). The Cosinor analysis will allow to obtain the phase (φ) and the amplitude (A) of the rhythm. A polar graph is usually performed to plot the results (see Figure 1.7 and online Appendix). In this graph, a circle represents the period T, with a vector that origins in the center of the circle. This vector has a size that depends on the amplitude (A) and its direction depends on the time of the acrophase (φ). Additionally, it is
The Cosinor method is widely used to find the amplitude (as criterion for significant rhythmicity), but more importantly, the acrophase of the time series. However, as many other curve-fitting methods, is strongly dependent on the signal’s waveform, and responds much better (both in accuracy and significance) with smooth sinusoidal data than signals with abrupt changes or nonsymmetrical shapes (Cornelissen, 2014).
The online Appendix includes examples of Cosinor analysis of locomotor activity series. To analyze this kind of data, first select the time period and then run the function cosinor (see the files SCRIPT.m and ScriptFilter.m for further details in book’s GitHub repository).
6.7 Circulars Statistics and Rayleigh
Circular statistics refers to the development of statistical techniques for the use with data represented on an angular scale (Berens, 2009). Due to their circular nature, many circadian data can be analyzed with circular statistics. The phase of time series, when we study a rhythmic signal, is an example of circular data. In this case, there is no designated zero and, in contrast to a linear scale, high/low values are arbitrary. This idea of ‘circular time’ is especially strong in the case of analyzing repeating patterns with a fixed time range (i.e., a cycle of period T and how this pattern changes). For this reason, the rhythmic signal is naturally represented as a circular graph (as shown in the Cosinor’s section). The rhythm’s phase can be studied as a fraction of the complete cycle (in time, angle or radians).
There are many other tests for circular statistics as an alternative for Rayleigh, such as the Omnibus test (Zar, 1999). This test works well with unimodal, bimodal and multimodal data, but its use is more selective and is not as extensive as others (Berens, 2009).
In case that there is a set of signals (e.g., many individuals in the same experiment), each one will have its own phase. We can study phase similarities by representing all the phases in the same circular graph and test them for a significant clustering by means of the Rayleigh test (Figure 1.7). This test asks how large the vectorial sum of the individual vectors must be to indicate a non-uniform distribution, according to the Fisher circular distribution (Fisher, 1995). The Rayleigh test is especially strong for detecting a unimodal deviation from uniformity (Berens, 2009). The test estimates the module (r) of the sum vector of the unit vectors corresponding to the phases of each individual. The higher the value of r, the greater degree of phase homogeneity. The p-value of
Examples of circular statistics can be found in the online Appendix section (Rayleigh.m file; see book’s GitHub repository).
The circadian system is ubiquitous, being essential for most living organisms’ physiology. The interaction of several circadian clock genes and proteins generates oscillations in expression of output target genes which temporally regulate numerous molecular and cellular processes. A precise regulation of circadian rhythms is crucial to keep track with the environment, and perturbations in the circadian system are linked to numerous molecular dysfunctions and may result in severe pathologies.
Hence, a comprehensive knowledge regarding the mechanistic of the circadian system is necessary to develop new procedures to investigate pathologies associated with a deregulated clock.
In the online appendix section (see book’s GitHub repository), we incorporate important tools for circadian data analysis, including codes and sample data. We include the codes for ls, sb, Cosinor, and Rayleigh analysis as well as a code to plot actograms. There are many other software – some of them freeware – that possess a wider range of circadian tools. We include links to some software too.
Devan B.D. E.H. Goad H.L. Petri E.A. Antoniadis N.S. Hong C.H. Ko L. Leblanc S.S. Lebovic Q. Lo M.R. Ralph & R.J. McDonald (2001). Circadian phase-shifted rats show normal acquisition but impaired long-term retention of place information in the water task. Neurobiology of learning and memory75(1) 51–62.
Hampp G. J.A. Ripperger T. Houben I. Schmutz C. Blex S. Perreau-Lenz I. Brunk R. Spanagel G. Ahnert-Hilger J.H. Meijer & U. Albrecht (2008). Regulation of monoamine oxidase A by circadian-clock components implies clock influence on mood. Current Biology18(9) 678–683.
Meng Q.J. E.S. Maywood D.A. Bechtold W.Q. Lu J. Li J.E. Gibbs S.M. Dupré J.E. Chesham F. Rajamohan J. Knafels B. Sneed L.E. Zawadzke J.F. Ohren K.M. Walton T.T. Wager M.H. Hastings & A.S. Loudon (2010). Entrainment of disrupted circadian behavior through inhibition of casein kinase 1 (CK1) enzymes. Proceedings of the National Academy of Sciences USA 107(34) 15240–15245.
Moriya T. Y. Yoshinobu Y. Kouzu A. Katoh H. Gomi M. Ikeda T. Yoshioka S. Itohara & S. Shibata (2000). Involvement of glial fibrillary acidic protein (gfap) expressed in astroglial cells in circadian rhythm under constant lighting conditions in mice. Journal of Neuroscience Research60(2) 212–218.
Yoo S.H. C.H. Ko P.L. Lowrey E.D. Buhr E.J. Song S. Chang O.J. Yoo S. Yamazaki C. Lee & J.S. Takahashi (2005). A noncanonical E-box enhancer drives mouse
Period2 circadian oscillations in vivo. Proceedings of the National Academy of Sciences USA, 102, 2608–2613.