Locus coeruleus noradrenergic neurons phase-lock to prefrontal and hippocampal infra-slow rhythms that synchronize to behavioral events

The locus coeruleus (LC) is the primary source of noradrenergic projections to the forebrain, and, in prefrontal cortex, is implicated in decision-making and executive function. LC neurons phase-lock to cortical infra-slow wave oscillations during sleep. Such infra-slow rhythms are rarely reported in awake states, despite their interest, since they correspond to the time scale of behavior. Thus, we investigated LC neuronal synchrony with infra-slow rhythms in awake rats performing an attentional set-shifting task. Local field potential (LFP) oscillation cycles in prefrontal cortex and hippocampus on the order of 0.4 Hz phase-locked to task events at crucial maze locations. Indeed, successive cycles of the infra-slow rhythms showed different wavelengths, as if they are periodic oscillations that can reset phase relative to salient events. Simultaneously recorded infra-slow rhythms in prefrontal cortex and hippocampus could show different cycle durations as well, suggesting independent control. Most LC neurons (including optogenetically identified noradrenergic neurons) recorded here were phase-locked to these infra-slow rhythms, as were hippocampal and prefrontal units recorded on the LFP probes. The infra-slow oscillations also phase-modulated gamma amplitude, linking these rhythms at the time scale of behavior to those coordinating neuronal synchrony. This would provide a potential mechanism where noradrenaline, released by LC neurons in concert with the infra-slow rhythm, would facilitate synchronization or reset of these brain networks, underlying behavioral adaptation.


Introduction
The brain coordinates activity among interconnected regions via coherent oscillatory cycles of excitation and inhibition (Womelsdorf et al., 2007). This can facilitate communication among selected subsets of neurons, groups of neurons, and brain regions. Sensory stimuli or behavioral events can reset the phase of these oscillations (Canavier, 2015;Voloh and Womelsdorf, 2016), linking activity of multiple neurons to process information in concert. However, the principal brain rhythms studied in behaving animals are at the time scale of cell neurophysiological processes, which are much faster (on the order of tens and hundreds of milliseconds) than real life behavioral events, which typically occur at second and supra-second time scales. The brain has several mechanisms linking these two time scales, some of which involve the hippocampus (reviewed in Banquet et al., 2021) and associated networks, including the prefrontal cortex and striatum.
Little is known about brain rhythms that operate in this crucial behavioral time scale during awake behavior. The brain is indeed capable of generating rhythms on the order of 0.1-1.0 Hz, although these have been principally characterized during sleep (Steriade et al., 1993). Furthermore, during sleep or under anesthesia, rat noradrenergic locus coeruleus (LC) and prefrontal cortical (Pfc) neurons are phase-locked to infra-slow rhythms (Lestienne et al., 1997;Eschenko et al., 2012;Totah et al., 2018). LC stimulation exerts powerful influence on neurophysiological activity in Pfc and hippocampus (Hip; Berridge and Foote, 1991). LC actions in prefrontal cortex are implicated in vigilance, decisionmaking, and executive function, while in Hip they are associated with learning and processing contextual information (e.g., Sara, 2009;Wagatsuma et al., 2018 for review). Since oscillations can coordinate activity in brain networks, we reasoned that there might also be rhythmicity on this behavioral time scale in awake animals, and investigated this possibility in rats performing a task engaging Pfc, Hip, and LC (Xiang et al., 2019;Oberto et al., 2022). Such coordinated activity could provide a possible link between neuromodulation and oscillatory coordination of brain areas on the time scale of behavior.

Materials and methods
All experiments were carried out in accordance with local (Comité d'éthique en matière d'expérimentation animale no. 59), institutional (Scientific Committee of the animal facilities of the Collège de France) and international (US National Institutes of Health guidelines; Declaration of Helsinki) standards, legal regulations (Certificat no. B751756), and European/national requirements (European Directive 2010/63/EU; French Ministère de l'Enseignement Supérieur et de la Recherche 2016061613071167) regarding the use and care of animals. The data analyzed here were recorded in experiments described by Xiang et al. (2019) and further details can be found there.

Animals
Nine male Long-Evans rats (Janvier Labs, Le Genest-Saint Isle France; weight, 280-400 g) were maintained on a 12 h:12 h lightdark cycle (lights on at 7 a.m.). The rats were handled on each workday. To motivate animals for behavioral training on the Tmaze, water was restricted except for a 10-30 min period daily to maintain body weight at 85% of normal values according to age. Rats were rehydrated during weekends. Food was restricted to 14 g of rat chow daily (the normal daily requirement) to prevent the animals from becoming overweight. Recording properties and behavioral correlates of the LC neurons of four of the rats were reported in Xiang et al. (2019), but no LFP data are presented there. Recordings of LC neurons during task performance was not possible in the remaining five rats included here.

The automated T-maze with return arms
The behavioral task took place in an elevated automated T-maze (see Figure 3A) consisting of a start area, a central arm, two reward arms and two return arms, which connected the reward arms to the start area. Small wells at the end of each reward arm delivered liquid reward (30 µl of 0.25% saccharin solution in water) via solenoid valves controlled by a CED Power1401 system (Cambridge Electronic Design, Cambridge, UK) with a custom-written script. As the rats crossed a central photo-detector, visual cues (VCs) were displayed in pseudo-random sequence on video monitors positioned behind, and parallel to the two reward arms. This is the "VC, " or "central arm PD" event. The VCs were either lit or dim uniform fields. The rat then selected the left or right arm and crossed the Reward arm (Rew) photodetector (PD), triggering reward release from an audible solenoid valve. Crossing the photodetector in the middle of the return arm of the T-maze (Return arm PD; Ret) triggered the visual cue to be turned off. Photodetectors detected task events and triggered cues and rewards via the CED Spike2 script.

Viral vector preparation, injection, and immunohistochemistry
The Canine Adenoviral vector (CAV2-PRS-ChR2-mCherry) was produced at the University of Bristol. Details about it and injections procedures appear in Xiang et al. (2019). Also see Xiang et al. (2019) for histological and immunohistochemical procedures and neuron characterization.

Electrode and optrode implants
Following VC task pre-training, at least 1 day before surgery, rats were returned to ad libitum water and food. Surgical procedures and electrode construction are the same as in Xiang et al. (2019). Moveable tungsten microelectrodes (insulated with epoxylite R , impedance = 2-4 M , FHC Inc., Bowdoin, ME, USA) were used for LC recordings. A single microelectrode, or two or three such electrodes glued together, was implanted at AP −3.8 to -4 mm relative to lambda, and ML 1.1-1.2 mm, with a 15 • rostral tilt. A stainless steel wire (Teflon coated, diameter = 178 µm, A-M Systems Inc., Sequim, WA, USA) implanted in the midbrain area about 1-2 mm anterior to the LC electrode tip served as a fixed LC reference electrode, permitting differential recording. The rat with the virus injection (R328) was implanted with an optrode composed of a tungsten microelectrode (insulated with epoxylite, impedance = 2-4 M , FHC Inc., Bowdoin, ME, USA) glued to a 200 µm optic fiber implant with a ferrule (0.37 numerical aperture, hard polymer clad, silica core, multimode, Thorlabs, Bergkirchen, Germany), with tip distances 1 mm apart (the electrode was deeper). The optic fiber implant and optic fiber cables were constructed at the NeuroFabLab (CPN, Ste. Anne Hospital, Paris). Two screws (diameter = 1 mm, Phymep, Paris) with wire leads were placed in the skull above the cerebellum to serve as ground. LC electrodes were progressively lowered under electrophysiological control until characteristic LC spikes were identified (located ∼ 5-6 mm below the cerebellar surface, see Xiang et al., 2019 for details). For the virus-injected rat, LC spikes could also be identified by responses to laser stimulations (described below). Following implantation, the microelectrode was fixed to a micro-drive allowing for adjustments along the dorsal-ventral axis. The headstage was fixed to the skull with dental cement, and surrounded by wire mesh stabilized with dental cement for protection and shielding. After the surgery, animals were returned to their home cages for at least 1week recovery with ad libitum water and food, and regular observation.

Electrophysiological recordings
Rats were then returned to dietary restriction. The movable electrodes were gradually advanced until a well-discriminated LC unit was encountered, and then all channels were recorded simultaneously while the rat performed in the T-maze. If no cells could be discriminated, the electrodes were advanced and there was at least a 2 h delay before the next recording session.
For daily online monitoring of LC spikes, pre-amplified signals were filtered between 300-3,000 Hz for verification on the computer screen (Lynx-8, Neuralynx, Bozeman, MT, USA) and also transmitted to an audio monitor (audio analyzer, FHC, Bowdoin, ME). For recordings, brain signals were pre-amplified at unity gain (Preamp 32, Noted Bt, Pecs, Hungary) and then led through a flexible cable to amplifiers (x500, Lynx-8, Neuralynx, Bozeman, MT, USA) and filters (0.1-9 kHz, Lynx-8, Neuralynx, Bozeman, MT, USA). Brain signals were digitized at ∼20 kHz using CED Power1401 converter and Spike2 data acquisition software (CED, Cambridge, UK). The LC unit activity was identified as in Xiang et al. (2019) by: (1) spike waveform durations ≥0.6 ms; (2) low average firing rate (1-2 Hz) during quiet immobility; (3) brief responses to unexpected acoustic stimuli followed by prolonged (around 1 s) inhibition; (4) for the virus-injected rat (R328), LC units were verified by responses to laser stimulation. A laser driver (Laserglow Technologies, North York, ON, Canada, wavelength 473 nm) was controlled by signals from a stimulator (Grass Technologies, West Warwick, RI, USA, Model SD9). Light intensity from the tip of optic fiber was measured by a power meter (Thorlabs, Bergkirchen, Germany, Model PM100D). If unit firing was entrained to the pulses with an increased rate (to at least twice the baseline firing rate) averaged over all the stimulations, they were considered to be noradrenergic LC units.
A light emitting diode (LED) was mounted on the cable that was plugged into the headstage. This was detected by a video camera mounted above the T-maze and transmitted to the data acquisition system at a sampling rate of ∼30 Hz for the purpose of position tracking.
2.6. Signal processing, spike sorting, and data analyses For off-line spike detection of LC activity in three of the rats, the wide-band signals were converted and digitally high-pass filtered (non-linear median-based filter). Waveforms with amplitudes passing a threshold were extracted, and then subjected to principal component analysis (PCA). All of these processes were performed with NDManager (Hazan et al., 2006). Spikes were sorted with a semi-automatic cluster cutting procedure combining KlustaKwik (KD Harris) and Klusters (Hazan et al., 2006). 1 Spikes with durations less than 0.6 ms were rejected. In one rat (R311), the LC signal was filtered from 300-3,000 Hz during recording, and the spike sorting was performed with Spike2 software (which employs a waveform template-matching algorithm). Most data analyses were performed using Matlab (R2010a) with the statistical toolbox FMAToolbox (developed by M. Zugaro) 2 and scripts developed in the laboratory, as well as some statistical analyses performed with Microsoft© Excel©. Phase was computed with the "Phase" function of the FMAToolbox, employing the Hilbert transform. To characterize periods with infra-slow rhythms, a criterion for salient phase-locking to task events was established as when the SEM range of LFP phase was less than 0.75 * π radians (cf., Figure 4, middle column). This is termed "regular phaselocking." Sessions tallied for phase-locking of LC neurons to infraslow rhythms were included only if they had at least 1,000 LC neuron spikes.

LC neuron phase-locking to prefrontal and hippocampal infra-slow rhythms
Infra-slow rhythms were readily apparent in visual inspections of hippocampal (Hip) and prefrontal cortical (Pfc) local field potentials (LFPs) (Figure 1A). These were rendered more salient by filtering the signal in a 0.1-1.0 Hz window ( Figure 1B). We applied an amplitude threshold to examine data from those periods when the infra-slow rhythm amplitude was elevated ( Figure 1C). This was intended to limit analyses to those periods when the infra-slow oscillation was sufficiently robust, thus avoiding possibly spurious computations of phase angle from low amplitude oscillations. In each of four rats, LC neurons were phase-locked to Pfc (e.g., Figures 1D, 2), as well as Hip infra-slow LFP rhythms. The incidence of phase-locking of the LC neurons in sessions with regular phase-locking of the infra-slow rhythms to task events was 20 out 23 for Hip LFP and 15/20 for Pfc LFP (Rayleigh test, p < 0.05; for histology, see Figure 2 of Xiang et al., 2019). In the animal where noradrenergic LC neurons were identified optogenetically (see Materials and methods), all were phase-locked to the infra-slow rhythms (n = 8 for both Pfc and Hip; Rayleigh test, p < 0.05).
The modal preferred infra-slow phase among these neurons was 0.35 * π radians for Hip infra-slow and 0.15 * π radians for Pfc infra-slow oscillations (p < 0.05, Rayleigh test; not shown).

Prefrontal and hippocampal infra-slow rhythms are synchronized to maze events
The infra-slow rhythms were phase-locked to task-relevant positions on the maze (Figure 3B and Supplementary Figure 1B). To quantify this phase-locking, the mean (± SEM) phase of the rhythm was plotted in peri-event time color plots (see Figure 4 and Supplementary Figure 1A for examples) over all trials in 57 sessions from nine rats (including the four with LC recordings). Infra-slow rhythms were phase-locked to the reward arm photodetector crossing (Rwd) in 51 of the recording sessions for Hip, and 46 sessions for Pfc (see Table 1). The other maze events had fewer incidences of regular phase-locking (Pfc return arm photodetector crossing, or Rtn: 11; Pfc central arm visual cue onset PD, or VC: 18; Hip Rtn: 18; Hip VC: 20). The mean phases at the respective PD crossings (in those cases when SEM ≤ 0.75 * π radians there) were 0.70 * π and 0.24 * π radians for Pfc and Hip Rtn, 0.25 * π and 0.22 * π radians for Pfc and Hip Rwd, and 0.19 * π and 0.01 * π radians for Pfc and Hip VC. The root-mean-square differences between Pfc and Hip mean phase (calculated pairwise by session) at the respective PD crossings were 0.14 * π, 0.13 * π, and 0.12 * π rad. The regular phase-locking could last from less than one to over 2.5 successive rhythmic cycles (Supplementary  Figures 1, 2 and Table 1) and could continue from one event to the next (Figure 4 and Supplementary Figure 1). For PL Rwd and Hip Rwd, 30 and 37 sessions had durations of regular phaselocking lasting one or more cycles, respectively. These permitted quantification of the temporal duration of the cycles, which ranged from 2.0 to 2.6 s, the equivalent of 0.4-0.5 Hz. In the six cases of Rwd PD phase-locking which had a second complete cycle, the mean of the first was 2.3 s, while the second was lower, 2.0 s (pairwise t-test, p = 0.0009, df = 5). Thus, these are not regular periodic oscillations, but, rather, this is consistent with phaselocking to task events. Pfc and Hip infra-slow rhythms sometimes resembled one another (e.g., Figure 2). To compare them, sessions were classified as having Pfc and Hip regular phase-locking in several ranges of cycles (see Table 1). In 17 of the 57 sessions, these numbers of cycles were different between Pfc and Hip for VC, Rwd and/or Rtn (e.g., Supplementary Figure 2). This indicates that it is unlikely that Pfc and Hip infra-slow rhythms are related by volume conduction, and suggests that they could be independently generated.
The infra-slow rhythms were regularly phase-locked to two (in 24 sessions), or even all three (in 8 sessions) different task events. Thus, they were not linked to any specific task-related behavior. To test whether infra-slow rhythms were triggered by rapid head movements, regression analysis compared the onset of regular phase-locking and times of peak acceleration, or deceleration around the Rwd PD crossings, and were not significant (r 2 = 0.034, p = 0.49 and r 2 = 0.0056, p = 0.80, respectively; df = 15; e.g., see Supplementary Figure 3). In Xiang et al. (2019), we showed that LC neurons fire more during accelerations. Indeed, the periods with the greatest increase in LC activity were not those most frequent for the start of regular phase-locking (i.e., reset) of the infra-slow rhythm; rather phase-locking occurred most frequently to Rwd PD crossing (see above), where no consistent accelerations occurred (see Supplementary Figures 3, 4). These results indicate A B

FIGURE 2
Spike phase-locking to infra-slow rhythms from two example LC neurons (A,B). Radius values are spike counts. Red arrows represent resultant vectors. it is unlikely that Pfc and Hip infra-slow rhythms are due to a biomechanical artifact, for example, from locomotion or head rocking.

Coordination of neuronal activity across time scales
In the four sessions where Hip and Pfc neurons could be discriminated from the LFP electrodes, most were also modulated by infra-slow rhythms (Pfc LFP modulated 6/12 Pfc units and 8/11 Hip units; Hip LFP modulated 8/12 Pfc units and 8/11 Hip units; Rayleigh test p < 0.05). The LC neurons could have relatively consistent phases with respect to the two infra-slow rhythms (not shown). LC neurons could be phase-locked to oscillations in the delta frequency range (1-4 Hz) in Pfc (n = 15/37) and Hip (11/37) as well as theta (5-10 Hz; 7/37 and 5/37, respectively) for Pfc and Hip (see Supplementary Figure 5). While phase-locking of LC neurons to gamma (40-80 Hz) was rare (n = 2 for both structures' LFPs), the infra-slow rhythm did modulate the amplitude of gamma oscillations at 35-45 Hz in hippocampus and prefrontal cortex (Figure 5).

Discussion
LFP oscillation cycles on the order of 0.4 Hz in prefrontal cortex and hippocampus were phase-locked to task events at crucial points on the maze. Successive cycles had different cycle lengths, indicating that, if they are indeed periodic oscillations, their phase can reset to salient events. Simultaneous recordings in prefrontal cortex and hippocampus could have different cycle lengths as well, while still phase-locking to task events. This would seem to exclude any single structure from entraining these independent rhythms simultaneously. This also would rule out a contribution of volume conduction. The intriguing issue of the origin of these rhythms merits further investigation. Most of the LC neurons were phaselocked to these infra-slow prefrontal cortical and hippocampal An example of simultaneous recordings of Pfc and Hip LFP infra-slow oscillations phase-locked to principal maze events, the PD crossings (at time zero). Each row of the color plots corresponds to a single trial. The phase of the infra-slow LFP is color-coded. Black rings correspond to the PD crossing prior to (left) or after (right) the event at zero for each plot. Note that the time scales vary among the events, in order to display prior and subsequent PD's. The dark traces in the 2nd and 4th rows show mean phase while the lighter shaded upper and lower lines are ±SEM. In the middle column, the blue vertical bars and blue double-headed arrow illustrate the calculation of the range of regular phase-locking (defined here as the period with the criterion of SEM range <0.75*π radians; pink double-headed arrows). Here, desynchronization (zones with wider SEM ranges) and discontinuities in the mean phase result from inter-trial variability in speed and distance from the synchronization point (PD-photodetector crossing). This is from the same session as the recording in Figure 3B. In the six cases of two or more cycles, only data from the first cycle were counted for mean cycle period and frequency. Cycles are only counted in the period from the previous trial event to the next one, even though successive cycles could extend before or after (cf., Figure 4 and Supplementary Figure 1).
LFPs, including all of the optogenetically identified noradrenergic neurons. Hippocampal and prefrontal units were also phase-locked to the infra-slow oscillations. While the number of LC neurons recorded may appear low, this is typical for the rare chronic recording studies of this structure in behaving rats, likely because its diminutive dimensions and deep location render accurate electrode placements challenging. This is consistent with previous work showing neuronal activity adapting to the time scale of behavioral events. For example, in behavioral tasks with delays, several brain structures show "time cell" activity: neurons with sequential "tiling" activity lasting on the order of several seconds. These periods can expand or contract depending upon the duration of task-imposed intervals (MacDonald et al., 2011). We speculate that this infra-slow rhythm may originate in the hippocampal-prefrontal system since neurophysiological activity there tracks time intervals on the order of several seconds based upon regularities in temporal structure of behavioral or environmental events. Steriade et al. (1993) observed infra-slow (0.3-1.0 Hz) rhythms in neocortical activity in anesthetized and naturally sleeping cats. Eschenko et al. (2012) showed that LC neuronal activity in sleeping rats is synchronized with the sleep infra-slow wave cycle (1 Hz) and is out of phase with Pfc neuronal activity. Similarly, in rats under ketamine anesthesia, there is a negative correlation between activity of LC NE neurons and prefrontal neurons, when neuron activation oscillates at ∼1 Hz (Sara and Hervé-Minvielle, 1995;FIGURE 5 Example of infra-slow modulation of gamma rhythm LFP in Pfc (top) and Hip (bottom). The black trace schematically illustrates the phase of the infra-slow rhythms. Gamma power was elevated between the phases of zero and pi radians of the infra-slow rhythm, while this was not apparent at other frequencies. a.u.-arbitrary units. Lestienne et al., 1997). Furthermore, when the latter authors pooled their LC recording data, they were significantly phase-locked to cortical LFP delta oscillations. While these infra-slow cycles of UP-DOWN state transitions are not generally observed in awake animals, this does demonstrate that LC can fire rhythmically, and that these structures can coordinate their activity at this time scale. Furthermore, in rats under urethane anesthesia, Totah et al. (2018) found that the firing rate of locus coeruleus neurons oscillates at 0.4-0.5 Hz. And, in head-fixed awake mice, cortical noradrenergic axons exhibited rhythmic Ca 2+ activity at 0.5-0.6 Hz (Oe et al., 2020). Thus, the LC could also be associated with the Pfc-Hip in the origin, maintenance and communication of behaviorally relevant infra-slow rhythms in the brain. Further work is required to elucidate the respective roles of these structures in these processes.
In the awake state, there is evidence for infra-slow neural processing although this was not observed as rhythms per se. Molter et al. (2012) observed a 0.7 Hz modulation of the power of theta rhythm recorded in rat hippocampus. This 0.7 Hz rhythm modulated hippocampal neuronal activity during sleep, as well as during behavior in a maze, a running wheel, and an open field. Positions on a figure-8 maze corresponded to specific phases of this modulatory rhythm, similar to the infra-slow rhythm recorded here. (Their filter settings excluded 0.7 Hz rhythms and thus this could not be directly measured in that work.) In Molter et al. (2012), the 0.7 Hz modulation of the power of the theta infra-slow modulation was locked at π radians to junction points in the maze (their Figure 7B), where accelerations might be expected. However, they found no overall correlation between phase and acceleration. Halgren et al. (2018) observed a rhythm at less than 3 Hz generated in the superficial layers of the cerebral cortex in awake humans. The phase of this rhythm reset to infrequent tones in their oddball task, similar to the reset of the infra-slow rhythm here in relation to salient task events. Villette et al. (2015) used calcium imaging to observe CA1 pyramidal cells in head-fixed mice moving in the dark on a non-motorized treadmill. They found that different neurons fired sequentially in cycles at the same time scale as the infra-slow oscillations observed here. Furthermore, the cycles could occur singly, or consecutively in groups of two or three. The authors interpreted this as representing an intrinsic metric for representing distance walked. This resembles time cell activity (Pastalkova et al., 2008;MacDonald et al., 2011) evoked above, where the length of the cycle extends to the time scale of the ongoing task (Kraus et al., 2013;Ravassard et al., 2013). The 2-5 s durations of the cycles in the Villette et al. (2015) study may represent a default value since their task had no temporal structure. This is on the order of the time scale of the infra-slow rhythm recorded here, and the variable numbers of cycles they observed might flexibly adapt to the positions of task-relevant events to lead to the results found here.
The present observations of phase-locking of LC neurons to infra-slow rhythms in the hippocampus could ostensibly be due to independent synchrony to task events by the infra-slow rhythms and the LC neurons. However, the LC neurons showed phase preferences in the infra-slow rhythms in data pooled over multiple task events. Furthermore, we did not observe any simple relation between infra-slow rhythms and motor events (e.g., as we showed for LC neurons with acceleration or deceleration in Xiang et al., 2019), since regular phase-locking could start before (Supplementary Figure 1) or after the same task events in different sessions (not shown), and continue over periods including a variety of behaviors.
The phase-locking of LC neurons to infra-slow rhythms in Hip and Pfc, as well as to oscillations in the delta, theta and gamma frequency bands could reveal coordinated neuronal processing within a unified temporal framework (cf., Totah et al., 2019). The time scale of the infra-slow rhythms corresponded to the temporal and spatial regularities characterizing the current behavioral patterns. Cross-frequency coupling could serve as a mechanism to link processing at different time scales. This could facilitate both "Communication through coherence" (CTC, Fries, 2005;Bosman et al., 2012) and "Binding by synchrony" (Eckhorn et al., 1990;Engel et al., 1999;Buehlmann and Deco, 2010). Thus, infra-slow rhythms would serve as a scaffold to link the time scales of dynamics of neuronal processes to those of behavior and cognitive processes. Noradrenaline, released by LC neurons in concert with the infra-slow rhythm, would participate in synchronizing or resetting those brain networks underlying behavioral adaptation to these events (Bouret and Sara, 2005;Sara and Bouret, 2012).

Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.