Detection of time-, frequency- and direction-resolved communication within brain networks

Electroencephalography (EEG) records fast-changing neuronal signalling and communication and thus can offer a deep understanding of cognitive processes. However, traditional data analyses which employ the Fast-Fourier Transform (FFT) have been of limited use as they do not allow time- and frequency-resolved tracking of brain activity and detection of directional connectivity. Here, we applied advanced qEEG tools using autoregressive (AR) modelling, alongside traditional approaches, to murine data sets from common research scenarios: (a) the effect of age on resting EEG; (b) drug actions on non-rapid eye movement (NREM) sleep EEG (pharmaco-EEG); and (c) dynamic EEG profiles during correct vs incorrect spontaneous alternation responses in the Y-maze. AR analyses of short data strips reliably detected age- and drug-induced spectral EEG changes, while renormalized partial directed coherence (rPDC) reported direction- and time-resolved connectivity dynamics in mice. Our approach allows for the first time inference of behaviour- and stage-dependent data in a time- and frequency-resolved manner, and offers insights into brain networks that underlie working memory processing beyond what can be achieved with traditional methods.

where Y(i) is a vector of the observed time series and ε Y (t) models the observational noise. The second equation describes the AR process. The a(i) are the matrix entries of A(i) rearranged into a vector. The matrix A(i) contains the interactions between the components of the AR process Z(t). The information about the network structure is therefore contained in this matrix. Time-resolved Granger-causality can be derived from the estimated parameters A(i) from which rPDC is calculated. However, while this has been validated previously based on simulated data 9 , its applicability to "real-world" quantitative EEG (qEEG) data remained to be confirmed. Validation of qEEG methods based on actual observed data is inherently more difficult than for simulated data as there is no known ground truth against which to compare results and performance. Therefore, we here first compare the novel approaches with their classical counterparts (i.e., AR to FFT and rPDC to coherence). We demonstrate that AR and FFT are equivalent for a chosen experimental example (age-related qEEG changes in mice) and confirm that the additional directional information provided by rPDC is congruous with results obtained via classical coherence.
We then probe the utility of AR-based methods, first in a preclinical application which does not require temporal resolution (drug effects on NREM sleep) followed by an example requiring high temporal resolution (spatial decision making). In the latter case, we demonstrate that time resolved rPDC can capture brief, behaviourally relevant connectivity patterns which underlie spontaneous alternation behaviour.

EEG surgeries. Animals were housed and tested in accordance with European (FELASA) and UK Home
Office regulations, experiments were approved by the university's Ethics Board and carried out in accordance with the Animal (Scientific Procedures) Act 1986. Electrode implantation surgery was conducted as described previously 15 . Briefly, general anaesthesia was induced using inhaled isoflurane and adequate anaesthetic depth was confirmed by lack of limb withdrawal response to interdigital pinch. In all experiments, gold screw electrodes were skull implanted into burr holes drilled in the parietal plate bilaterally above the dorsal hippocampi (2 mm posterior/1.5 mm lateral (right + left) of Bregma) and into the frontal plate above the right prefrontal cortex (2 mm anterior to Bregma immediately lateral to midline). In Experiment 2 an additional recording electrode was placed above the left somatosensory cortex (SSC; 0.2 mm posterior to Bregma, 3.3 mm lateral to midline). Experiment 3 utilized four recording sites, i.e. bilateral for both prefrontal cortex and parietal/ hippocampal sites (coordinates as above). Reference and ground electrodes were placed at neutral locations superficial to parietal and occipital cortices. Subcutaneous buprenorphine hydrochloride was administered for pain management prior to withdrawal of anaesthesia. Animals were then allowed 2 weeks convalescence prior to testing.
EEG recordings were obtained under freely moving conditions using wireless data logging devices, i.e. either the Neurologger (New Behaviour, CH; refs 15,16 ) or a neural activity tracker (NAT-1) device (Cybula, UK; ref. 17 ). The Neurologger (Experiment 1) has an AC input range of ±750 μV, and input is band-pass limited (1-70 Hz); data were sampled at 200 Hz. The NAT-1 device has an AC input range (configurable) of ±1.75 mV, the input is band pass limited (0.3-500 Hz), and data were sampled at 1000 Hz in Experiment 2 and 200 Hz in Experiment 3. Attachment of the device was conducted 1 hour prior to testing to mitigate handling stress and allow the animal to appropriately habituate to the device. 1) The "traditional" Fast Fourier Transform (FFT) based estimation was run in MatLab 2015a and used to determine the spectrum via smoothing of the periodogram as described in 18 . A triangular window with a smoothing width of h = 1 Hz was chosen to obtain a consistent estimator. 2) The second spectral estimator was based on autoregressive (AR) processes 19 . Briefly, an AR model is fitted to the data and from its coefficients, the respective spectrum is calculated. The model order used here was set as one quarter of sampling rate i.e. p = 50 samples (for data sampled at 200 Hz) and p = 250 samples (for data sampled at 1 kHz). In both cases, the model order therefore corresponds to 250 ms.
Coherence was also determined based on: 1) A "classical" coherence approach, i.e. cross-periodograms were calculated and smoothed (smoothing width = 1 Hz), and 2) A multivariate autoregressive (AR) process to estimate directed Granger causality between data, yielding renormalized partial directed coherence (rPDC) as described in 12 . For rPDC estimation, we again used a model order of sampling rate/4. A combination of rPDC and state space modelling allows a time and frequency-resolved measure for Granger causality 14 .
Statistical analysis. EEG Frequency band definitions were Delta (1-5 Hz), Theta (5-9 Hz), Alpha (9-14 Hz), Beta (14)(15)(16)(17)(18)(19)(20) and Gamma . Statistical comparison of FFT and AR based power spectra between groups/conditions were conducted via non-matching or repeated measures 2-way ANOVAs as appropriate. For within subject comparisons, mean band power (within subject average of power measured across frequency points in band) was also compared using paired 2-tailed t-tests. Differences are reported here for both methods (α = 0.95). Group/condition differences for coherence and rPDC were analysed for the same bands based on the respective distributions under the alternative. For rPDC this is a chi-square distribution with two degrees of freedom. For a group of n independent measurements the sum of the rPDCs is thus chi-square distributed with 2*n degrees of freedom. We used this distribution to generate confidence intervals for each group. u.], means + SEM) illustrated over the 1-50 Hz range for mice aged 5 m (black) and 13 m (red). Left: FFT-based power spectra for each recording site, right: corresponding spectra calculated using autoregressive spectral estimation (AR). Asterisks on plots provide a p-value summary for the effect of age (*p < 0.05, **p < 0.01, ***p < 0.001). (b) Mean "classical" coherence strength [a.u., left column] observed between all possible channel pairs over the 1-50 Hz range for mice aged 5 m (black) and 13 m (red). Corresponding rPDC strength [a.u.] over the same frequency range for all possible channel pairs and directions of coherence are depicted in the middle and right columns. Pictograms on the left illustrate which channel pairs are analysed in each row. Plot titles indicate compared channels and, where applicable, direction of coherence. R.PFC = right prefrontal, R.H = right parietal (hippocampus), L.H = left parietal (hippocampus). Blue lines = 95% confidence limit. Lightly coloured lines indicate the confidence interval of the group means. Significant differences between age groups are illustrated for individual frequency bands (*p < 0.05). In the case of coherence, the distribution is more complex. However, it was shown 20 that tanh −1 [Coh(ω)] is approximately normally distributed. The mean coherence for a group of n independent measurements is therefore also approximately normally distributed. From this distribution confidence intervals for the mean coherence of each group can be derived. For coherence as well as rPDC, based on the confidence intervals for two groups, we can then check for a significant difference by counting at how many frequencies the results for the two groups overlap. If this number is smaller than what is expected at random, we conclude that a significant difference exists. Thresholds for significance were defined as α = 0.95 both for classical coherence analysis and for rPDC.
For time-resolved rPDC, differences between groups were analysed based on the mean of each group. Standard error of the mean (SEM) defines a confidence interval for each group. Differences between groups of more than one sigma are investigated.

Experiment 1: Comparison of FFT and AR based methods.
Male wild type (PLB1 WT ) mice 16 were surgically implanted at ~3 months of age. Briefly, 24 hour EEG recordings were conducted at 5 and 13 months of age with animals individually housed in Pheno Typer home cage observation systems (Noldus IT, NL) for the duration of recording. Automated sleep staging of recorded EEG with post hoc manual correction was then carried out using SleepSign (Kissei Comtec Co. Ltd, JP) to identify periods of wakefulness, as well as REM and NREM sleep. Secondary confirmation of wakefulness was provided through detection of movement by the device's on-board accelerometer. For this study, 6 mice from each age group were selected for analysis based on optimal quality of EEG recordings in all channels. A 100 second artefact-free EEG strip (during wakefulness) was extracted for each animal. Power spectra were then calculated using both FFT and AR analysis and pooled by age group. Cross-channel coherence was calculated using both classical and rPDC methods. Experiment 2: Drug effects on AR spectra and rPDC. Ten male PLB1 WT mice were implanted at ~4 months of age to investigate the effects of ketamine and diazepam on sleep and EEG profiles. A 4-week Latin square based drug administration regime was constructed where animals were allocated to 4 groups. Once per week each group received an intraperitoneal injection of either 10 mg/kg ketamine hydrochloride (Tocris Bioscience, UK) in vehicle (saline), saline only (ketamine vehicle), diazepam (Sigma Aldrich, UK) in a vehicle composed of a 1:1:18 ratio of ethanol, Cremaphor and saline, or the diazepam vehicle only. Treatment group allocations were rotated such that each animal received each of the 4 treatments with a 1 week washout between treatments. Following administration, recording devices were attached and animals were placed in individual holding cages where EEG activity was recorded for 6 hours. Five animals were selected for data analysis where high quality recordings had been obtained from recordings conducted both after administration of a drug treatment and administration of its respective vehicle. Sleep staging of EEG was conducted as in experiment 1. Here, a 60 second artefact free strip of EEG from periods of NREM sleep was extracted for each animal and condition. Power spectra and cross-channel coherence were then assessed using AR and rPDC analyses, respectively, and results were pooled by treatment group.

Experiment 3: Time resolved rPDC analysis of spontaneous alternation in the Y-Maze. Five
female PLB1 WT mice were implanted at 6 months of age for Y-maze behavioural testing. As in our previous investigations 21 the Y-maze was constructed from white acrylic and consisted of 3 open topped corridors (arms), each 60 cm long, meeting at a triangular intersection termed 'the central zone' (CZ). Each arm (labelled A, B and C) was subdivided into an inner (I) zone (30 cm proximal to CZ) and an outer (O) zone (30 cm distal to CZ). The demarcation of these zones is presented in Suppl. Fig. 1a. At the start of each trial the animal was placed at the distal end of arm B (BO) facing toward the CZ. Upon release, the motion tracking software was activated and the animal was allowed to freely explore the maze for 10 minutes. Transitions between maze arms during the trial were recorded both automatically by the software and manually by an observer viewing the trial on a monitor in a segregated cubicle. The animal was considered to have entered an arm when all 4 paws had moved from the CZ to one of the inner zones. Transitions between arms (with the exception of the first two) were scored as either correct or incorrect alternations by the observer. A correct alternation was defined by a triplet of arm entries where the animal visited all 3 arms in sequence without repetition (ABC, BCA etc.). An incorrect alternation (or alternate arm re-entries) was defined by a triplet where the first and third entry were the same arm (ABA, BCB etc.). A third category, termed direct repeats, was defined as a triplet which contained one or more consecutive entries into the same arm (AAB, BAA etc.).
Consistent with expectations, correct alternations were the most common triplet type accounting for 58.69% of triplets (SEM = 2.95), followed by incorrect alternations (34.49%, SEM = 4.44) and direct repeats (6.81% SEM = 3.30). Due to the rarity of direct repeats (absent in 2 animals) they were not considered for EEG analysis.
Motion tracking and event time stamping. Animal movement was recorded by a ceiling mounted camera offering a planar view of the maze. Camera output was digitized using a RTV-24 frame grabber PCI card (ADLINK Technology Inc.). Motion was tracked by a PC running the ANYMAZE 4.9 (Ugo Basile Srl, IT) video tracking package. ANYMAZE was configured such that transitions between zones resulted in the generation of TTL level output trigger pulses unique for the 7 zones of the maze. These pulses were sent to an Arduino Uno rev.2 micro-controller which triggered a unique sequence of on-off pulses to drive flashes from an infra-red (IR) LED array via a MOSFET driver. These IR flashes were detected by an IR receiver module on the NAT-1 device and used to time stamp zone transition events onto the EEG signal.
Extraction of behaviourally relevant EEG segments. Data from the NAT-1 devices (4 EEG channels and IR channel) were downloaded as CSV files. Time resolved coordinates representing the animal's location were exported from ANYMAZE (also in CSV format) along with TTL time stamps to spatially and temporally synchronize data, i.e. the animal's location and the EEG time series. This procedure was accomplished using an in-house written MATLAB (MathWorks Inc, USA) script. Six second segments of EEG data corresponding to either correct or incorrect transitions (see above) were identified and extracted for analysis.
To control for gross differences in movement speed and location, strips were taken from arm transitions where the animal had moved from the outer extremity of one arm to the outer extremity of another during the extraction time frame. This step was taken to ensure that the motion paths of animals were comparable for all EEG   Figure 1c. AR power spectra for each alternation type were calculated, averaged within subjects and subsequently pooled. An rPDC analysis for the above data set was conducted using both time stationary and time resolved techniques. For the time resolved analysis, rPDC strength was calculated, resolved in both time and frequency for all 30 EEG samples. Results were pooled to calculate the population median coherence strength observed over all channel pairs. The data were then partitioned into correct and incorrect alternations and grouped heat plots created to provide a visual comparison of changes in coherence strength between alternation types. The population median coherence value was used as a colour scale normalization factor to ensure that grouped heat plots of correct and incorrect alternations were directly comparable.

Results
Experiment 1: Effect of age in wild-type mice. Regional analysis: FFT vs AR Spectra. Power spectra from 100 second EEG strips were calculated using both FFT and AR based methods and compared between young (5 m) and aged (13 m) wild-type mice (PLB1 WT ) (Fig. 1a & Table 1). Visual inspection indicates an improved smoothness (and hence identification of peaks) of AR vs. FFT spectra, likely due the reduced impact of noise. However, for the purpose of quantitative analysis, both methods were statistically equivalent, and detected global reductions in power with age for most frequency bands. These findings confirmed that AR modelling results in accurate reconstruction of the underlying power spectrum and detection of age related changes.
Connectivity: Classical coherence vs rPDC. Connectivity measures were compared between young and aged mice for all possible channel pairs using both 'classical' coherence and rPDC (Fig. 1b). The former confirmed robust connectivity between all channel pairings for both groups and an effect of age regarding a reduction intra-hippocampal connectivity in all frequency bands except gamma. By contrast, a number of additional age-related effects were detected using rPDC: left to right hippocampal connectivity was reduced in all frequency bands except Delta, while Theta coherence between these channels was enhanced in the opposite direction.  Enhanced Theta coherence was also observed in aged mice for more complex connections (i.e. right pre-frontal cortex to the contralateral hippocampus) as well as enhanced Delta coherence from the right prefrontal cortex to both hippocampi.   Figure 2). In the 6 hrs following drug injection, diazepam increased NREM sleep and reduced waking activity, while opposing trends were seen in ketamine treated mice.
Diazepam reduced NREM Delta/Theta power and increased anterior Gamma power. NREM AR power spectra obtained from diazepam/ketamine treated mice alongside respective vehicle controls (Fig. 2 & Tables 2-3) suggest that diazepam significantly reduced NREM Delta and Theta power relative to vehicle controls in all 4 recording sites, while Gamma power in the right prefrontal cortex as well as Beta and Gamma power in the left somatosensory cortex increased significantly. Bilaterally increased Gamma band power was also indicated in the parietal recording sites, this was not confirmed for mean band power. By contrast, despite multiple significant effects of ketamine on frequency-resolved bands, no differences were detected for mean band power.
Diazepam decreased, while ketamine increased, prefrontal -hippocampal functional connectivity. Drug effects on connectivity (rPDC; Fig. 3) suggested a significant reduction in coherence strength for diazepam vs vehicle treatment. In direct contrast, ketamine exclusively increased interregional coherence. This difference was most striking in communication from the right prefrontal cortex to right hippocampus where diazepam reduced coherence across all frequency bands while ketamine increased coherence in all bands except Theta. A similar contrast was observed in Gamma coherence between the same channels in the opposing direction. Ketamine also significantly increased Beta coherence from the left hippocampus to the right prefrontal cortex, and Alpha coherence from the left-to-right hippocampus. Additional band specific reductions in connectivity were also observed under diazepam. Diazepam decreased Delta coherence from the left-to-right hippocampus, Delta and Theta coherence from the somatosensory to the right prefrontal cortex (and Beta coherence in the opposing direction), as well as Beta coherence from the left somatosensory cortex to the right hippocampus.
For methodological comparison, classical coherence analysis was also performed (Suppl. Figure 3). Only one significant drug effect was noted, an elevated theta coherence between the left hippocampus and somatosensory cortex under diazepam. Interestingly, no significant drug effects emerged in this range between these channels based on rPDC. A trend toward elevated theta coherence under diazepam was however noted in the left hippocampus to somatosensory cortex direction (Fig. 3). Experiment 3: AR and rPDC analysis of brain activity during spontaneous alternation. AR spectra during correct vs. incorrect alternations. AR power spectra were calculated for 15 samples of EEG capturing correct and 15 capturing incorrect Y-maze alternations then averaged within subjects (n = 5; Fig. 4a & Tables 4, 5). A significantly higher prefrontal Gamma power was observed in incorrect vs correct alternations. A similar trend was observed in the right hippocampal recording site. Significant differences at lower frequencies were indicated for all channels, however these were not corroborated by mean band power analysis.
Directional connectivity (rPDC) during correct and incorrect alternations. Time invariant functional connectivity (rPDC) over the 1-50 Hz range during correct and incorrect alternations in the Y-maze (Fig. 4b) suggest coherence to be highest (particularly in the Delta-Theta range) between left and right hippocampi, and from the hippocampi to the ipsilateral prefrontal cortices. Differences detected between correct and incorrect alternations were entirely confined to the Gamma band, and detected from the right hippocampus to the left hippocampus and left prefrontal cortex, from the left hippocampus to the left prefrontal cortex and from the left prefrontal cortex to the right hippocampus. However, in all cases, Gamma coherence was affected in a bi-directional manner and consequently, lines intersected repeatedly. Additionally, the magnitude of coherence observed in the Gamma range in these analyses was close to (and occasionally below) the 95% confidence limit at multiple frequency points. Therefore, results relating to the Gamma band should be considered with caution.
Time resolved effective connectivity during correct vs. incorrect alternations. Cognitive processes in a navigating animal must be expected to be highly dynamic. This was indeed demonstrated using heat plots of time-resolved rPDC over the 6 seconds analysed for correct and incorrect alternations (Suppl. Figures 4 and 5). From these plots, the greatest distinction between correct and incorrect alternations emerged in contralateral homotypic (left-to-right hippocampus) and ipsilateral heterotypic (right prefrontal to right hippocampus) directions and in a time window between 1.5 and 4 seconds. This corresponds to the animals' passage through the central zone of the maze (Suppl. Figure 1c). In order to illustrate these differences, plots of coherence (for the frequency range 1-30 Hz) over the identified time period are presented for interhemispheric homotypic connections in Fig. 5a and for ipsilateral heterotypic connections in Fig. 5b.
Qualitative examination of interhemispheric homotypic coherence during central zone crossings revealed an interesting contrast in the coherence patterns of correct vs. incorrect transitions: Between the hippocampi, directed coherence in the 5-15 Hz range was markedly stronger in the left-to-right direction than right-to-left during correct alternations. During incorrect alternations the same coherence was stronger in the opposing direction. Furthermore, the inverse pattern was seen in the prefrontal cortex: 5-15 Hz coherence was stronger in the right-to-left direction than in the left-to-right during correct alternations. Again, during incorrect alternations the opposite was true. This effect was also visible in plots of ipsilateral heterotypic coherence, i.e. 5-10 Hz coherence from the hippocampus to the prefrontal cortex was strongest in the left hemisphere during correct alternations but in the right hemisphere during incorrect alternations.
Subsequent statistical comparison of time resolved coherence between correct and incorrect alternations revealed multiple time-frequency regions where groups differed in rPDC strength by more than one sigma (see Fig. 5e). Most salient among these was a specific temporal pattern of rPDC (highlighted in Fig. 5c) clearly apparent in correct and absent from incorrect alternations: 1. A burst of right prefrontal to right hippocampal Theta coherence, strongest at ~6 Hz at 1.6-1.8 seconds. 2. Simultaneous bursts of high Theta and Beta coherence from the left to right hippocampus, strongest at ~9 Hz and ~16 Hz, respectively, observed from 1.9 seconds onwards. Similar bursts were also observed ~ 3.1 s in both correct and incorrect alternations but remained stronger in the correct case. 3. A sustained period of Theta coherence from the right hippocampus to right prefrontal cortex at ~6 Hz from 2-3 s. This coincided with the period between intra-hippocampal Theta/Beta coherence bursts described in 2.
The rPDC strength (mean ± SEM) at each of the key frequencies (6, 9 or 16 Hz, indicated by black lines in Fig. 5c and 5e) is illustrated over time for both correct and incorrect alternations to visualize both magnitude of difference and variability (Fig. 5d).

Discussion
This paper sought to validate improved qEEG algorithms in key neuroscientific applications. Deciphering directional changes as well as direct vs indirect routes of communication between channels was found to be critical for our understanding of network phenomena. Previous attempts to provide time resolution of directional connectivity typically employed non-linear auto-regression based approaches [22][23][24] . However, these approaches were limited in their temporal resolution and only easily applied to bivariate analyses. We therefore propose that the time-resolved combination of rPDC with state space modelling offers a superior solution, whilst also dealing with non-linear and noisy properties of biopotentials.
Comparison of classical functional connectivity (coherence) with directed functional connectivity (or effective connectivity; rPDC) indicated logical agreements between the methods. For example, a generally high level of coherence was observed between contralateral homotypic sites densely anatomically interconnected by callosal/commissural fibres 25 . Yet, rPDC produced a wealth of additional details and was able to elaborate on hippocampal-prefrontal network communications and their modulation by physiological states. Some connectivity patterns are particularly difficult to resolve with traditional methods, for example, ipsilateral prefrontal-hippocampal coherence, where rPDC indicated significant coherence in both directions. As the prefrontal -hippocampal connection is anatomically indirect [26][27][28] communication along this tract can be assumed to be underestimated by classical coherence analysis.
Effect of age. AR spectra confirmed a global age-dependent reduction of EEG power for all frequencies at prefrontal channels, and at frequencies above Delta in the hippocampal channel. The first salient finding from the rPDC analysis was an increased Delta coherence from the prefrontal cortex to the ipsilateral and contralateral hippocampus. Increased Delta coherence is a frequent observation in human AD and has been associated with cholinergic dysfunction 29,30 . Reports on the effect of normal physiological aging specifically on Delta coherence are however absent from the literature. More generally, interhemispheric coherence between homotypic regions reliably decreases in resting awake subjects across multiple frequency bands with normal aging 31,32 and is thought to result from loss of callosal projections 31,33 . This is consistent with our observation of reduced left -right hippocampal connectivity in aged mice and with results of classical coherence analysis. However, Theta coherence was found to be enhanced in aged mice in the right-left direction under rPDC. As classical coherence essentially sums coherence in both directions, the opposing Theta changes observed here were likely rendered invisible thus far. Additionally, right-left coherence was markedly stronger in the opposing direction for both age groups. rPDC results therefore indicate that a symmetrical loss of callosal projection fibres between the hippocampi is an incomplete explanation. Rather, age-related degradation of hippocampus dependent memory in rodents 34,35 may involve both loss of function alongside compensatory increase in specific frequency bands. rPDC can therefore offer new avenues to study hemisphere-specific changes and functional lateralization (e.g., refs 36,37 ). Critically, findings provided by rPDC do not violate the predictions of classical coherence nor those of existing literature, i.e. summation of changes in rPDC over both directions would result in the cancellation of opposing changes, and thus similar connectivity estimates to those provided by classical methods.

Pharmaco-EEG. Spectral analysis indicated that diazepam reduced Delta-Theta activity and enhanced
Gamma activity in the prefrontal and somatosensory cortices. This is a known effect of diazepam often termed "the benzodiazepine fingerprint" 38 . By contrast, the anticipated increase in NREM Delta band power under ketamine 39 was not detected, likely attributable to dosage and timing differences.
Striking differences between the actions of ketamine and diazepam were confirmed via rPDC analysis particually for the ipsilateral prefrontal-hippocampal coherence: Diazepam reduced this coherence in all frequency bands and Gamma coherence in the opposing direction, while ketamine significantly increased coherence strength for these directions and frequencies. This implies a weakening of prefrontal -hippocampal network coupling under diazepam and strengthening of the same under ketamine. Reductions in intra-hemispheric Theta and Alpha coherence have been observed in awake humans in response to diazepam 40 , while decreased inter-parietal coherence has been reported in rats 41 . The descent into NREM sleep is associated with a reduction in cortical connectivity and eventually a functional decoupling of the prefrontal cortex [42][43][44] . Our findings therefore indicate that the enhanced induction and maintenance of NREM sleep under diazepam may result from facilitated decoupling of the prefrontal cortex. Conversely, the increased NREM onset latency observed under ketamine may result from a stimulating effect of enhanced prefrontal coupling. These results are particularly exciting as they demonstrate altered communication between specific brain structures in response to a systemically administered drug.  Furthermore, the opposing action of these drugs on prefrontal-hippocampal network coupling as predicted by rPDC are completely in accordance with their opposing actions at the behavioural level.

Time-resolved qEEG tracks networks during spontaneous alternation.
A key challenge in cognitive neuroscience is the assessment of transient cognitive processes and behavioural states. EEG offers high temporal resolution but limitations imposed by analytic methods have squandered this advantage. Here, we demonstrate that time-stationary rPDC analyses during spontaneous alternation in the Y-Maze yielded no overt effect of alternation type on inter-regional connectivity and only a subtle elevation of hippocampal Gamma power. Time-resolved rPDC however demonstrated an astounding level of dynamism in network connectivity within a short time frame. Furthermore, we discovered that correct choices are uniquely accompanied by a specific temporal pattern of Theta coherence on approach to the maze intersection. The so far undetected pattern of a prefrontal-hippocampal initiation signal expands on the previously suggested role of Theta communication between the two regions, thought to be critical for cognitive function 45 . Synchronization during decision making in working memory tasks [46][47][48] was assumed to be primarily hippocampally-driven as suggested by phase locking analysis 45 . However, this assumption is not definitive, particularly when there is an indirect anatomical route between recording sites. Enhanced hippocampal-prefrontal Theta coherence was previously proposed in a reward-based Y-Maze paradigm and a role in learning suggested 49 . Our data now confirms Theta coherence in the same maze location in the absence of a reward, and stipulates that it uniquely precedes 'correct' spontaneous alternations, in line with previous observations 48 . Interestingly, a recent study has indicated modulation of hippocampal encoding during route planning via the prefrontal-thalamic-hippocampal pathway 50 . Our data suggest that -regardless of whether a decision is made based on anticipation of reward 49 or spontaneously -similar neural mechanisms are responsible for executing this choice.
Overall, we propose that our analyses provide compelling evidence that rPDC provides exciting new opportunities for a wide range of neuroscience applications. For example, it is now possible to determine the time course and anatomical conduits through which representations of space encoded by hippocampal place and grid cells 51 interact with prefrontal executive control functions. Within the realm of neurodegenerative research, hippocampal-prefrontal directional connectivity can also be probed in e.g. AD models and MCI patients. This will ultimately provide valuable information regarding the impact of disease processes on network properties, and may offer a wealth of new biomarkers for diagnosis.  Table 5. Effect of alternation type on EEG band power. Summary of paired 2-tailed t-test statistics comparing mean band power of mice during correct and incorrect Y-Maze alternations (within subjects design). Degrees of freedom (DF) as well as T and F values are provided for each EEG band.