Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 1 of 38 1 Diverse phase relations among neuronal rhythms 2 and their potential function 3 4 Eric Maris*, Pascal Fries#*, Freek van Ede*^ 5 *Radboud 6 EZ, Nijmegen, The Netherlands 7 #Ernst 8 Planck Society, 60528, Frankfurt, Germany 9 ^Oxford Centre for Human Brain Activity, Department of Psychiatry, University of 10 University, Donders Institute for Brain, Cognition and Behaviour, 6525 Strüngmann Institute (ESI) for Neuroscience in Cooperation with Max Oxford, OX3 7JX, Oxford, United Kingdom 11 12 Corresponding Author: 13 Maris E. ([email protected]) 14 15 Keywords: neuronal oscillations, correlated neuronal activity, phase relations, 16 travelling waves, selective neuronal communication. 17 18 The total number of words of the manuscript and the original three boxes, 19 excluding abstract, glossary, new box (nr 1) and figure legends: 4773 20 The number of words of the abstract: 106 21 The number of boxes (excluding Glossary): 4 22 The number of figures: 3 1 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity 23 Page nr. 2 of 38 Key figure: Number 3 24 25 26 Acknowledgements: and very useful input. 27 28 29 30 The authors would like to thank the two expert reviewers for their sharp FvE was supported by The British Academy and The Royal Society (Newton International Fellowship). PF was supported by the Human Connectome Project (WU-Minn 31 Consortium, NIH grant 1U54MH091657), a European Young Investigator 32 Award, the European Union (HEALTH F2 2008 200728) and the LOEWE 33 program (“Neuronale Koordination Forschungsschwerpunkt Frankfurt”) 34 2 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 3 of 38 35 Abstract 36 Neuronal oscillations at nearby sites in the brain often show phase relations that 37 are consistent across time, yet diverse across space. We discuss recent 38 demonstrations of this phase-relation diversity, and show that, contrary to 39 earlier beliefs, this diversity is a general property of oscillations that is neither 40 restricted to low-frequency oscillations, nor to periods outside of stimulus 41 processing. Arguing for the computational relevance of phase-relation diversity, 42 we discuss that it can be modulated by sensory and motor events, and put 43 forward the idea that phase-relation diversity may support effective neuronal 44 communication by (1) enhancing selectivity and (2) allowing for the concurrent 45 segregation of multiple information streams. 46 3 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 4 of 38 47 Existence and potential relevance of phase-relation 48 diversity 49 A central question in neuroscience is how the billions of neurons in the human 50 brain are coordinated such that they perform useful computations. Looking for 51 an answer, it is obvious to consider the fact that neuronal activity is usually 52 coordinated across time and space. This holds both for sub-threshold membrane 53 potentials (reflecting the inputs to a neuron) and action potentials (spikes, the 54 neuronal output). Correlations among neuronal output, such as spike synchrony 55 and relative spike timing, have a substantial impact on neuronal function [1-7]. 56 Correlations across time often occur within a limited frequency band, as typically 57 identified by a rhythmic pattern in the auto-correlogram (see Glossary). This 58 rhythmic neuronal activity is typically denoted as a neuronal oscillation. Neuronal 59 oscillations are involved in a whole range of sensory [8], motor [9, 10], and 60 cognitive processes [11-13]. They have been described in terms of their 61 frequency, amplitude, synchronization and between-site phase relations. Here, 62 we focus on the surprising diversity in their between-site phase relations, even 63 between nearby sites (see Box 1, for a description of how these phase relations 64 can be assessed empirically). Historically, this diversity has mainly been studied 65 in the context of traveling waves, and was considered typical for ongoing 66 oscillations in the absence of stimulus processing. In this paper, we discuss 67 recent demonstrations of phase-relation diversity, showing that it is a general 68 property of oscillations that is neither restricted to low-frequency oscillations, 69 nor to periods outside of stimulus processing. 4 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 5 of 38 70 71 To argue for the functional relevance of phase-relation diversity, it is important 72 to demonstrate (1) that it can be modulated by sensory and motor events, and 73 (2) that it plays a role in the flow of neuronal information. With respect to the 74 role of neuronal oscillations in the flow of neuronal information, much attention 75 has gone to interactions between communicating sites. Specifically, it has been 76 proposed that effective communication between two oscillating neuronal 77 populations depends, at least in part, on their phase relation [14-17]. 78 Communication is most effective when the neuronal output of the sending 79 population arrives at the receiving population at its most excitable phase. 80 Because oscillations in both the sending and the receiving populations may be 81 characterized by local diversity in their phase relations (as we will review), it 82 makes sense to incorporate this type of diversity into models of how neuronal 83 oscillations contribute to selective routing of information. As a first attempt to 84 this, toward the end of this review, we will propose two ways by which such 85 diversity may support effective neuronal communication: (1) by enhancing 86 selectivity and (2) by allowing for the concurrent segregation of multiple 87 information streams. 88 89 In the following, we will first discuss recent observations demonstrating that 90 phase-relation diversity is a general property of oscillations. Next, we highlight 91 several recent demonstrations of phase-relation diversity being modulated by 92 sensory and motor events. Finally, we review recent studies that investigated the 93 potential role of phase relations in selective neuronal communication. 5 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 6 of 38 94 95 Phase-relation diversity is a general property of neuronal 96 oscillations 97 Neuronal activity is often correlated across space. Strong evidence for correlated 98 neuronal input comes from optical recordings using voltage-sensitive dyes [VSDs; 99 18] and from dual patch-clamp recordings [19], two methods which are only 100 weakly or not affected by volume conduction and therefore ideal to assess 101 correlation. VSD images show orderly spatio-temporal patterns in subthreshold 102 membrane potentials, of which travelling waves are the most intensively studied 103 [20]. 104 105 When neuronal activity is correlated across space, this correlation can include 106 systematic time- and phase relations. Importantly, neuronal activity can be cross- 107 correlated (“cross” denotes “across locations”) with zero or non-zero time 108 relations; in the latter case, the cross-correlations are said to be time-lagged. 109 Depending on their frequency content, these patterns can be better described in 110 either the time or the frequency domain. A single-waveform (non-rhythmic) 111 travelling wave can be better described by its trajectory in the time domain, 112 whereas rhythmic standing and travelling waves can be better described in the 113 frequency domain. We will use the terms time and phase relations 114 interchangeably; when using the latter, their dependence on frequency will be 115 implicit (see Box 2). In the current review, we will consider traveling waves as a 116 special case of the more general phenomenon of phase-relation diversity, in 6 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 7 of 38 117 which phase relations between sites are diverse, regardless of whether they vary 118 systematically with space (see Box 3). Systematic phase relations can be ongoing 119 or can be evoked by a sensory or motor event. In the latter case, the associated 120 neuronal activity is often denoted as an evoked response, which is to be 121 distinguished from ongoing neuronal activity of which the phases are not locked 122 to some event. We will not distinguish between systematic phase relations in 123 evoked or ongoing neuronal activity, although their functions may be different. 124 125 Reliable phase-relation diversity (schematically depicted in Fig. 1A) has been 126 shown on multiple spatial scales, ranging from a few millimeters (using VSD 127 imaging and extracellular wire recordings) to over ten centimeters (using electro- 128 corticography, electro-encephalography, and magneto-encephalography). In 129 investigations of travelling waves, the focus has been on the low-frequency 130 oscillations that are prominent in the absence of task-relevant sensory 131 processing. For example, studies based on extracranial recordings typically 132 focused on travelling alpha waves [21-26], which are prominent in the absence of 133 visual stimulation. These waves appear to be widespread, travelling from 134 occipital to frontal sites over more than ten centimeters, although they are likely 135 to also contain volume conduction effects from much more local travelling waves 136 of primary currents in occipito-parietal cortex [21]. A prominent review article on 137 the subject expresses the view that this type of phase-relation diversity is 138 “typically present during periods outside of stimulation, while synchronous 139 activity dominates in the presence of a strong stimulus” [20]. In contrast to this 7 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 8 of 38 140 view, we now review recent studies that demonstrate robust phase-relation 141 diversity during sensory stimulation and in periods of cognitive engagement. 142 143 First, phase-relation diversity was demonstrated for band-limited oscillations in 144 the gamma-frequency range during sensory processing of task-relevant visual 145 stimuli [27, 28]. Specifically, under similar visual stimulation conditions, gamma- 146 band phase-relation diversity was observed both on a very small spatial scale 147 (inter-electrode distance <900 𝜇m) in macaque V4 [27], and on a large spatial 148 scale (using magneto-encephalography) in humans [28; see Fig. 1A and 1B]. 149 Importantly, this gamma-band phase-relation diversity was not only observed 150 between the site-specific local field potentials (LFPs), but also between the multi- 151 unit spiking activity observed at these same sites [27], implying that the phase- 152 relation diversity is also visible in the synaptic input to the receiving areas. 153 154 Second, gamma-band phase-relation diversity was also observed across the 155 layers of a V1 cortical column [29], and this diversity was moreover shown to 156 have a spatial pattern that is consistent with feed-forward processing: current 157 sinks occurred first in input layer 4 and propagated to deep and superficial layers 158 of cortex [30]. In addition, the latter study showed that the temporal progression 159 across the layers was the opposite for alpha-band oscillations: current sinks 160 occurred first in the feedback-recipient layers 1, 2 and 5 and propagated to the 161 input layer 4. Thus, the spatial pattern in the phase-relation diversity across 162 layers provided evidence for the gamma-band rhythm indexing feed-forward and 163 the alpha-band rhythm indexing feedback processing. 8 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 9 of 38 164 165 Third, phase-relation diversity was demonstrated in one of the most prominent 166 mammalian rhythms: the hippocampal theta rhythm. More precisely, it was 167 demonstrated that the hippocampal theta rhythm in the CA1 layer is a travelling 168 wave with dominant movement direction along the septo-temporal axis [31], 169 with the maximum spatial phase difference (between the septal and the 170 temporal pole) ranging from 0 to 𝜋 [32]. Crucially, it was also shown that, due to 171 the locking of the CA1 neurons’ spikes to their local theta phase, exactly the 172 same travelling wave was also observed in the spiking activity [31]. Because the 173 different hippocampal subregions along the septo-temporal axis project to 174 different targets (retrosplenial and perirhinal cortex for the septal subregion and 175 hypothalamus, lateral septum, amygdala and medial prefrontal cortex for the 176 temporal one), travelling waves ensure that these distinct hippocampal targets 177 receive peak CA1 input in a particular order. This can be relevant for how these 178 target areas integrate hippocampal input (see also Phase-relation Diversity and 179 the Flow of Neuronal Information). 180 181 Fourth and last, there is not only diversity in the between-site phase relations at 182 a particular frequency, but also in the between-site phase relations involved in 183 phase-amplitude coupling (see Fig. 1C). Using ECoG, recorded in humans during 184 working-memory encoding and maintenance, strong diversity was observed in 185 the between-site phase relations between low-frequency phases and high- 186 frequency amplitude envelopes [33]. The diversity in these cross-frequency 187 coupling phases was mainly due to diversity in the low-frequency phase 9 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 10 of 38 188 relations; bursts of high-frequency amplitudes tended to be synchronized [33]. 189 Some of the diversity in these low-frequency phase relations has a travelling 190 wave pattern [34]. 191 192 Interactions with sensory and motor events 193 Evidence for the functional relevance of phase-relation diversity comes from 194 studies that show how phase-relation diversity is affected by sensory input and 195 motor output. Related to this, we consider the question how well a stimulus 196 configuration can be decoded from the pattern of between-site phase relations. 197 A relevant variable in this respect is stimulus strength, and we therefore also 198 discuss the results of recent modeling work on the relation between stimulus 199 strength, neuronal oscillation frequency, and phase-relation diversity. Although 200 all reviewed studies demonstrate that phase-relation diversity exists, not all of 201 them also quantified its strength [specifically, this was not the case for 35, 36, 202 37]. 203 204 In one study, we investigated the effect of a visual grating stimulus (versus a 205 blank screen) on phase-relation diversity in macaque area V4 [27]. It is well 206 known that, in visual cortical areas, the presentation of a visual stimulus reduces 207 the alpha band and increases the gamma band power [38]. Using LFPs recorded 208 from four electrodes at the corners of a 650 by 650 µm square, we observed a 209 decrease in the alpha-band phase-relation diversity and an increase in the 210 gamma-band phase-relation diversity. That is, phase-relation diversity increased 10 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 11 of 38 211 with power, opposite to what would be expected if, at this spatial scale, power 212 increases would only result from increases in lag-zero synchronization. 213 214 In a related study, gamma-band phase-relation diversity was demonstrated at 215 the level of isolated single units recorded in V1 [39]. The results of this study 216 pertain to spike-field coherence, the phenomenon that neuronal spiking is 217 concentrated at a particular phase of the LFP. For the gamma-band LFP, spikes 218 preferably occur at the trough of the gamma phase [40-43]. This study 219 investigated the preferred phases of isolated single units relative to an aggregate 220 gamma band LFP obtained from the other electrodes (i.e., excluding the 221 electrode with which the unit was recorded). Subsequently, it demonstrated that 222 the diversity in these spike-LFP phase relations could be predicted by variables 223 that index neuronal activation relative to the local population: stimulus 224 orientation (relative to the unit’s preferred orientation), time after stimulus 225 onset, and spike density (calculated over a 250 ms window centered at the unit 226 spike) [39]. Specifically, for all these indices, with increasing neuronal activation, 227 the preferred phase shifted forward in the gamma cycle. 228 229 Stimulus-induced shifts in spike-LFP phase are not always observed. For instance, 230 they were not observed in a study that evaluated the effects of stimulus contrast 231 on different properties of gamma band synchronization [44]. They were also not 232 observed in a study that investigated the spiking of putatitve inhibitory neurons 233 in relation to a pooled gamma phase (pooled across multiple recording sites in 234 V4) [45]. On the positive side, besides the two studies that were already 11 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 12 of 38 235 mentioned [27, 39], two more studies reported shifts in spike-LFP phase as a 236 function of visual stimulation [46, 47]. 237 238 Phase-relation diversity is also triggered by motor events. This was demonstrated 239 in a recent study that investigated neuronal activity in macaque V4 as a function 240 of its time relative to saccadic eye movements [35]. This study showed that 241 saccadic eye movements evoke a wave with a single wavefront, travelling across 242 the V4 retinotopic map, originating from the foveal representation and travelling 243 towards the periphery (Fig. 2A). Importantly, these waves co-occur with a 244 reorganization of the post-saccadic firing rates, which follow the same 245 retinotopic pattern as the LFPs, and therefore may contribute to a temporal 246 pattern in the prioritizing of visual stimuli. 247 248 Highly relevant for a possible functional role of phase-relation diversity would be 249 a quantification of the degree to which a stimulus configuration can be decoded 250 from the pattern of between-site phase relations. Such a quantification of the 251 decoding performance was recently performed using recordings from the CA1 252 region in rat hippocampus [36]. The hippocampus receives input from all sensory 253 modalities and uses this information, amongst others, to determine the animal’s 254 position relative to landmarks in the environment. The firing rate of the 255 hippocampal neurons varies as a function of the animal’s position, forming so- 256 called place fields [48], allowing a decoder to determine the animal’s position 257 from the joint activity of multiple hippocampal neurons [49]. Crucially, it has now 12 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 13 of 38 258 been demonstrated that the animal’s position could be equally well determined 259 from the diversity in the theta-band phase relations [36; see Fig. 2B]. 260 261 Finally, there is recent modeling work on the relation between stimulus strength, 262 neuronal oscillation frequency, and phase-relation diversity. Here, the starting 263 point is the observation that the frequency of gamma-band oscillations in visual 264 cortex increases with stimulus strength, as implemented experimentally by 265 stimulus contrast [44, 50, 51]. Modeling studies have shown that differences in 266 excitatory drive not only determine gamma peak frequency but also between- 267 site phase relations [52-54]. Specifically, varying a whole range of model details 268 (involving conductance-based, Izhikevitch, and phase oscillator network models), 269 it was demonstrated that, within a population of model neurons oscillating at the 270 same frequency, neurons that received the stronger excitatory drive were phase- 271 advanced relative to those that received a weaker drive [53]. This demonstrates 272 a possible computational mechanism via which differences in excitatory drive are 273 translated into phase-relation diversity. 274 275 The relation between neurons’ intrinsic frequencies and their between-neuron 276 phase relations is well understood in the theory of weakly coupled oscillators [for 277 a review, see 55]. The basic ingredients of the theory are (1) a drive to each of 278 the network nodes (model neurons) which manifests itself as a node-specific 279 intrinsic frequency, and (2) a coupling strength between all model neurons that 280 co-determines whether a given neuron pair will synchronize (possibly at a non- 281 zero phase, as determined by their relative intrinsic frequencies). Thus, via 13 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 14 of 38 282 network-level synchronization, different input strengths may be translated into 283 non-zero phase relations. The quantitative relation between input strength, 284 synchronization frequency, and between-neuron phase relations, is 285 characterized by so-called Arnold tongues [55; see Fig. 2C]. An important 286 question is whether Arnold tongues can also be identified in vivo and, if so, 287 whether synchronization frequency and between-neuron phase relations are 288 modulated by input strength as specified by the theory of weakly coupled 289 oscillators. And further, because selective attention increases both firing rates 290 [56] and gamma peak frequency [57], the question is whether selective attention 291 modulates between-neuron phase relations as prescribed by the Arnold tongues. 292 Both questions require further empirical research. 293 294 Phase-relation diversity and the flow of neuronal 295 information 296 The influence of sensory and motor events on phase-relation diversity calls for an 297 investigation of the role such diversity may have in the routing of information in 298 the brain. A recent study investigated how the activity of a local neuronal group 299 depends on its preceding phase relation to another group [58]. The authors 300 investigated gamma band (50-80Hz) phase relations in macaque area V1 (with 301 electrodes sampling within a 16 mm2 cortical patch) while the animal watched 302 short video clips. They used a direction-specific information-theoretic measure 303 [transfer entropy; 59] with which they quantified how much a between- 304 population gamma band phase relation (the gamma phase of a sending 14 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 15 of 38 305 population relative to the one of a receiving population) influences the spiking 306 activity of the receiving population above and beyond what can be predicted by 307 the past firing rate of the receiving population itself. The study demonstrates 308 positive transfer entropy, indicating a causal effect (in the Wiener-Granger 309 sense) of the gamma band phase relation (sending-versus-receiving) on the firing 310 rate in the receiving location. Moreover, the study showed that Granger- 311 causation predominantly flows from the neuronal group that leads the phase 312 relation to the one that lags. These effects are likely functionally relevant, 313 because visually induced changes in gamma band phase relations resulted in an 314 increase in transfer entropy when they were in same direction as the average 315 phase relation. This study thus suggests that transient shifts in gamma band 316 phase relations mediate a dynamic reconfiguration of the pattern of causal 317 interactions (but see Box 4). 318 319 This study [58] also raises the question of how a neuronal network with fixed 320 anatomical connections (and thus also fixed axonal conduction delays) is able to 321 produce phase relation shifts on such a short time scale (i.e., the time scale at 322 which different movie segments trigger different phase relations). A recent 323 computational study showed that such stimulus-induced shifts in phase relations 324 can be produced by a network of connected neuronal populations if the intra- 325 population inhibition is larger than the inter-population inhibition [60]. This 326 study showed further, that the sign of the phase relation indicates the dominant 327 direction of information flow: information predominantly flows from the location 15 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 16 of 38 328 with the earlier phase to the location with the later phase, as confirmed 329 experimentally in [58]. 330 331 Phase relations cannot only be affected by bottom-up stimulation, but also by 332 top-down cognitive control. We will briefly consider work that has revealed an 333 important role for the alignment (and misalignment) of oscillatory phase in the 334 filtering of relevant and irrelevant information [37, 61]. In this work, the authors 335 investigated neuronal entrainment to stimuli repeating quasi-rhythmically at 1.5- 336 2 Hz, a form of temporal attention that optimizes stimulus processing for the 337 time window in which a relevant stimulus will appear [37]. They demonstrated 338 phase-relation diversity across different frequency-tuned sites in primary 339 auditory cortex (A1) and – crucially – its modulation by top-down attention. This 340 work built on previous work that demonstrated a frequency-specific phase reset 341 in A1: pure tones that correspond to a site’s preferred or non-preferred 342 frequency reset local oscillatory activity to its high, respectively, low excitability 343 phase [61]. Now, using dual rhythmic auditory streams, of which only one was 344 attended, it was demonstrated that the temporal structure of the LFP phase 345 reflected the attended frequency: differently tuned A1 regions were entrained 346 with opposite phases, with the phase of the attended frequency corresponding 347 to high neural excitability (as reflected by a firing rate increase) [37]. Thus, this 348 pattern of phase-relation diversity results in both stronger and more selective 349 neural responses at the attended time points, effectively functioning as a 350 spectro-temporal filter. 351 16 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 17 of 38 352 The study described in the previous paragraph has revealed that phase-relation 353 diversity of neural oscillations can separate relevant from irrelevant external 354 stimuli by virtue of the temporal predictability of the input streams. Although 355 this remains largely untested, the same concept may also apply in the context of 356 a neuronal population that receives input from both a relevant and an irrelevant 357 lower-area population. This requires two ingredients: (1) the higher-area 358 population must be able to predict the timing of inputs from the lower-area 359 population, and (2) temporally segregated outputs of the lower-area 360 populations. Predictability in the higher-area population follows from the lower- 361 area populations sending their output rhythmically. And temporal segregation in 362 the lower-area populations’ output follows from phase-relation diversity across 363 these populations. This configuration may allow a higher-area population to filter 364 the incoming information from multiple phase-diverse lower-area populations by 365 means of a phase alignment to the input from the relevant lower-area 366 population. 367 368 Phase-relation diversity may thus provide a mechanism for achieving selectivity 369 in inter-areal communication. This may be particularly relevant for oscillations 370 that are coherent over a wide area. For such widely coherent oscillations, if there 371 would be no diversity in the phase relations across space, it would be difficult to 372 conceive how selectivity could be achieved. This is also depicted schematically in 373 Figure 3A. However, in the presence of phase-relation diversity – either in the 374 lower- (as depicted in Fig. 3B) or the higher-area population – the information 375 could potentially be filtered for a smaller population within the larger population 17 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 18 of 38 376 that oscillates coherently. As such, phase-relation diversity may introduce 377 selectivity in inter-areal communication. Finally, if we consider that phase- 378 relation diversity may be prominent in both the lower- and the higher-area 379 populations, then this leads to the possibility that such diversity may support the 380 concurrent segregation of multiple communication channels, as schematically 381 depicted in Figure 3C. 382 383 This hypothesis is untested and therefore speculative. However, two recent 384 studies have provided support for the notion that the alignment – as well as the 385 misalignment – of phases in distant populations may indeed shape inter-areal 386 communication and vary dynamically as a function of task-demands. 387 Interestingly, both of these studies have demonstrated this for oscillations in the 388 beta frequency band that are coherent across wide parts of the cortex and in 389 particular between frontal and parietal brain areas. One study employed a visual 390 working memory task requiring an occulo-motor response [62]. This task had 391 previously been shown to engage coherent beta oscillations between prefrontal 392 and posterior parietal areas [63]. Interestingly, this study showed that this beta- 393 band inter-areal coherence exhibits phase-relation diversity, with the phase 394 relations forming a bimodal distribution, with modes at 0 and ±𝜋 [62]. 395 Moreover, it was demonstrated that these inter-areal phase relations vary with 396 task-demands. Another study also recorded from parietal and frontal areas (i.e., 397 the parietal reach region and the dorsal premotor area), but now in the context 398 of a planned movement task that had to be performed with either a hand or the 399 eyes [64]. In this task, monkeys maintained a movement plan for one second 18 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 19 of 38 400 until a go signal instructed movement execution. These authors also observed 401 pronounced beta-band inter-areal coherence, but in their dataset this occurred 402 predominantly at phase relations close to 0. Critically, however, when 403 statistically controlling for the contribution of the dominant parietal LFP 404 component (using a partial coherence analysis), they found that there are 405 periods of the task in which the spikes of parietal neurons actually show a phase- 406 opposition with the premotor cortex LFP. This pattern was particularly 407 pronounced during the movement planning phase and was much stronger when 408 a saccade was planned as compared to when a hand movement was planned. 409 This phase-opposition may contribute to the down-regulation of the 410 communication between these parietal and frontal areas at the stage of the 411 motor planning (when the movement plan should not yet be communicated to 412 output areas). This study therefore suggests that the dynamic regulation of inter- 413 areal phase-relation diversity may not only serve to enhance selectivity and 414 segregate information streams, but also to temporarily and selectively inhibit 415 effective inter-areal interactions. 416 417 An outstanding question is to what extent these observations generalize to 418 oscillations in other frequency bands and/or brain areas. A crucial requirement 419 for the proposed mechanisms (selection and segregation, depicted in Fig. 3B and 420 3C) is coherence between neuronal populations, both the interacting and the 421 non-interacting ones. Because of this requirement, it is uncertain whether these 422 mechanisms can hold for some oscillations, such as gamma-band oscillations 423 over cortical areas that span more than a few millimeters (e.g., V1, V2 and V4 19 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 20 of 38 424 combined). For the latter oscillations, part of the available data suggests that 425 selectivity is not implemented through phase-relation diversity, but through the 426 selective presence of coherence (see Box 4). 427 428 Conclusion 429 We have reviewed a series of recent studies that demonstrated that phase- 430 relation diversity is a general property of oscillations that is neither restricted to 431 low-frequency oscillations, nor to periods outside of stimulus processing, as 432 suggested by the earlier studies. Further, arguing for the functional role of 433 phase-relation diversity, we have discussed that (1) this phase-relation diversity 434 is modulated by sensory and motor events, to the extent that a stimulus 435 configuration can be decoded from it, and (2) that it is a potentially relevant 436 parameter in determining the flow of neuronal information, allowing for flexible 437 information selection and segregation mechanisms. One core target for future 438 research will be the identification of the precise conditions (in terms of 439 anatomical connectivity motif and network input) under which phase-relation 440 diversity is relevant for effective inter-areal communication (see Outstanding 441 Questions Box). 442 443 1. Branco, T., B.A. Clark, and M. Hausser, Dendritic discrimination of 444 temporal input sequences in cortical neurons. Science, 2010. 329(5999): p. 445 1671-5. 20 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity 446 2. 3. 4. 451 452 Markram, H., W. Gerstner, and P.J. Sjöström, A history of spike-timingdependent plasticity. Frontiers in synaptic neuroscience, 2011. 3. 5. 453 454 London, M. and M. Häusser, Dendritic computation. Annu. Rev. Neurosci., 2005. 28: p. 503-532. 449 450 Branco, T. and M. Häusser, Synaptic integration gradients in single cortical pyramidal cell dendrites. Neuron, 2011. 69(5): p. 885-892. 447 448 Page nr. 21 of 38 Caporale, N. and Y. Dan, Spike timing-dependent plasticity: a Hebbian learning rule. Annu. Rev. Neurosci., 2008. 31: p. 25-46. 6. Azouz, R. and C.M. Gray, Dynamic spike threshold reveals a mechanism 455 for synaptic coincidence detection in cortical neurons in vivo. Proceedings 456 of the National Academy of Sciences, 2000. 97(14): p. 8110-8115. 457 7. 458 459 neural information. Nature reviews neuroscience, 2001. 2(8): p. 539-550. 8. 9. Engel, A.K. and P. Fries, Beta-band oscillations - signalling the status quo? Current Opinion in Neurobiology, 2010. 20(2): p. 156-165. 462 463 Koepsell, K., et al., Exploring the function of neural oscillations in early sensory systems. Frontiers in Neuroscience, 2010. 4: p. 53. 460 461 Salinas, E. and T.J. Sejnowski, Correlated neuronal activity and the flow of 10. Davis, N.J., S.P. Tomlinson, and H.M. Morgan, The role of beta-frequency 464 neural oscillations in motor control. Journal of Neuroscience, 2012. 32(2): 465 p. 403-404. 466 11. 467 468 469 Ward, L.M., Synchronous neural oscillations and cognitive processes. Trends in Cognitive Sciences, 2003. 7(12): p. 553-559. 12. Düzel, E., W.D. Penny, and N. Burgess, Brain oscillations and memory. Current Opinion in Neurobiology, 2010. 20(2): p. 143-149. 21 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity 470 13. Page nr. 22 of 38 Bosman, C.A., C.S. Lansink, and C.M.A. Pennartz, Functions of gamma- 471 band synchronization in cognition: from single circuits to functional 472 diversity across cortical and subcortical systems. European Journal of 473 Neuroscience, 2014. 39(11): p. 1982-1999. 474 14. Börgers, C., S. Epstein, and N.J. Kopell, Gamma oscillations mediate 475 stimulus competition and attentional selection in a cortical network 476 model. Proceedings of the National Academy of Sciences, 2008. 105(46): 477 p. 18023-18028. 478 15. Fries, P., A mechanism for cognitive dynamics: neuronal communication 479 through neuronal coherence. Trends in Cognitive Sciences, 2005. 9(10): p. 480 474-480. 481 16. Tiesinga, P., J.M. Fellous, and T.J. Sejnowski, Regulation of spike timing in 482 visual cortical circuits. Nature Reviews Neuroscience, 2008. 9(2): p. 97- 483 109. 484 17. 485 486 Fries, P., Rhythms for Cognition: Communication through Coherence. Neuron, 2015. 88(1): p. 220-35. 18. Wu, J.Y., H. Xiaoying, and Z. Chuan, Propagating waves of activity in the 487 neocortex: what they are, what they do. Neuroscientist, 2008. 14(5): p. 488 487-502. 489 19. 490 491 Yu, J. and D. Ferster, Membrane potential synchrony in primary visual cortex during sensory stimulation. Neuron, 2010. 68(6): p. 1187-1201. 20. Ermentrout, G.B. and D. Kleinfeld, Traveling electrical waves in cortex: 492 insights from phase dynamics and speculation on a computational role. 493 Neuron, 2001. 29(1): p. 33-44. 22 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity 494 21. Hindriks, R., M.J. van Putten, and G. Deco, Intra-cortical propagation of EEG alpha oscillations. Neuroimage, 2014. 103: p. 444-53. 495 496 Page nr. 23 of 38 22. Ito, J., A.R. Nikolaev, and C. van Leeuwen, Spatial and temporal structure 497 of phase synchronization of spontaneous alpha EEG activity. Biological 498 cybernetics, 2005. 92(1): p. 54-60. 499 23. Nunez, P.L., B.M. Wingeier, and R.B. Silberstein, Spatial‐temporal 500 structures of human alpha rhythms: theory, microcurrent sources, 501 multiscale measurements, and global binding of local networks. Human 502 brain mapping, 2001. 13(3): p. 125-164. 503 24. 504 505 and correlations with responses. PloS one, 2012. 7(6): p. e38392. 25. 506 507 Fellinger, R., et al., Evoked traveling alpha waves predict visual-semantic categorization-speed. Neuroimage, 2012. 59(4): p. 3379-3388. 26. Burkitt, G.R., et al., Steady-state visual evoked potentials and travelling waves. Clinical Neurophysiology, 2000. 111(2): p. 246-258. 508 509 Patten, T.M., et al., Human cortical traveling waves: dynamical properties 27. Maris, E., et al., Rhythmic neuronal synchronization in visual cortex entails 510 spatial phase relation diversity that is modulated by stimulation and 511 attention. Neuroimage, 2013. 74: p. 99-116. 512 28. van Ede, F., et al., Both ongoing alpha and visually-induced gamma 513 oscillations show reliable diversity in their across-site phase relations 514 Journal of Neurophysiology, 2014. 515 29. Livingstone, M., Oscillatory firing and interneuronal correlations in 516 squirrel monkey striate cortex. Journal of Neurophysiology, 1996. 517 75(cb5e1d62-641d-b3bc-81b4-6adfa592e067): p. 2467-2552. 23 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity 518 30. Page nr. 24 of 38 van Kerkoerle, T., et al., Alpha and gamma oscillations characterize 519 feedback and feedforward processing in monkey visual cortex. 520 Proceedings of the National Academy of Sciences, 2014. 111(40): p. 521 14332-14341. 522 31. 523 524 travelling waves. Nature, 2009. 459(7246): p. 534-539. 32. 525 526 Lubenov, E.V. and A.G. Siapas, Hippocampal theta oscillations are Patel, J., et al., Traveling theta waves along the entire septotemporal axis of the hippocampus. Neuron, 2012. 75(3): p. 410-7. 33. van der Meij, R., M. Kahana, and E. Maris, Phase-Amplitude Coupling in 527 Human Electrocorticography Is Spatially Distributed and Phase Diverse. 528 Journal of Neuroscience, 2012. 32(1): p. 111-123. 529 34. Bahramisharif, A., et al., Propagating neocortical gamma bursts are 530 coordinated by traveling alpha waves. J Neurosci, 2013. 33(48): p. 18849- 531 54. 532 35. 533 534 visual cortex. Neuron, 2015. 85(3): p. 615-27. 36. 535 536 Agarwal, G., et al., Spatially distributed local fields in the hippocampus encode rat position. Science, 2014. 344(6184): p. 626-30. 37. 537 538 Zanos, T.P., et al., A sensorimotor role for traveling waves in primate Lakatos, P., et al., The spectrotemporal filter mechanism of auditory selective attention. Neuron, 2013. 77(4): p. 750-61. 38. Fries, P., et al., The effects of visual stimulation and selective visual 539 attention on rhythmic neuronal synchronization in macaque area v4. 540 Journal of Neuroscience, 2008. 28(18): p. 4823-4835. 24 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity 541 39. 40. 41. 546 547 Fries, P., et al., Modulation of oscillatory neuronal synchronization by selective visual attention. Science, 2001. 291(5508): p. 1560-1563. 544 545 Vinck, M., et al., Gamma-Phase Shifting in Awake Monkey Visual Cortex. Journal of Neuroscience, 2010. 30(4): p. 1250-1257. 542 543 Page nr. 25 of 38 Jia, X., S. Tanabe, and A. Kohn, Gamma and the coordination of spiking activity in early visual cortex. Neuron, 2013. 77(4): p. 762-774. 42. Ray, S. and J.H. Maunsell, Network rhythms influence the relationship 548 between spike-triggered local field potential and functional connectivity. 549 The Journal of Neuroscience, 2011. 31(35): p. 12674-12682. 550 43. Gregoriou, G.G., et al., High-frequency, long-range coupling between 551 prefrontal and visual cortex during attention. Science, 2009. 324(5931): p. 552 1207-10. 553 44. Ray, S. and J.H.R. Maunsell, Differences in Gamma Frequencies across 554 Visual Cortex Restrict Their Possible Use in Computation. Neuron, 2010. 555 67(5): p. 885-896. 556 45. Vinck, M., et al., Attentional modulation of cell-class-specific gamma- 557 band synchronization in awake monkey area v4. Neuron, 2013. 80(4): p. 558 1077-1089. 559 46. Havenith, M.N., et al., Synchrony Makes Neurons Fire in Sequence, and 560 Stimulus Properties Determine Who Is Ahead. Journal of Neuroscience, 561 2011. 31(23): p. 8570-8584. 562 563 47. Wang, P., et al., Time delays in the β/γ cycle operate on the level of individual neurons. NeuroReport, 2010. 21(11): p. 746-750. 25 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity 564 48. Page nr. 26 of 38 O'Keefe, J. and J. Dostrovsky, The hippocampus as a spatial map. 565 Preliminary evidence from unit activity in the freely-moving rat. Brain 566 research, 1971. 34(1): p. 171-175. 567 49. Brown, E.N., et al., A statistical paradigm for neural spike train decoding 568 applied to position prediction from ensemble firing patterns of rat 569 hippocampal place cells. Journal of Neuroscience, 1998. 18(18): p. 7411- 570 7425. 571 50. V2 by dynamic frequency matching. Neuron, 2013. 78(3): p. 523-36. 572 573 Roberts, M.J., et al., Robust gamma coherence between macaque V1 and 51. Jia, X., D. Xing, and A. Kohn, No consistent relationship between gamma 574 power and peak frequency in macaque primary visual cortex. The Journal 575 of Neuroscience, 2013. 33(1): p. 17-25. 576 52. Tiesinga, P.H. and T.J. Sejnowski, Mechanisms for phase shifting in cortical 577 networks and their role in communication through coherence. Frontiers in 578 human neuroscience, 2010. 4. 579 53. Lowet, E., et al., Input-Dependent Frequency Modulation of Cortical 580 Gamma Oscillations Shapes Spatial Synchronization and Enables Phase 581 Coding. PLoS computational biology, 2015. 11(2): p. e1004072-e1004072. 582 54. 583 584 Eur J Neurosci, 2014. 39(5): p. 705-19. 55. 585 586 587 Cannon, J., et al., Neurosystems: brain rhythms and cognitive processing. Pikovsky, A., et al., Synchronization: A universal concept in nonlinear sciences. Vol. 2. 2002: Cambridge University Press Cambridge. 56. Reynolds, J.H., T. Pasternak, and R. Desimone, Attention increases sensitivity of V4 neurons. Neuron, 2000. 26(3): p. 703-714. 26 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity 588 57. Page nr. 27 of 38 Bosman, C.A., et al., Attentional stimulus selection through selective 589 synchronization between monkey visual areas. Neuron, 2012. 75(5): p. 590 875-88. 591 58. Besserve, M., et al., Shifts of Gamma Phase across Primary Visual Cortical 592 Sites Reflect Dynamic Stimulus-Modulated Information Transfer. PLoS 593 Biol, 2015. 13(9): p. e1002257. 594 59. 595 596 2000. 85(2): p. 461. 60. 597 598 Battaglia, D., et al., Dynamic effective connectivity of inter-areal brain circuits. PLoS Comput Biol, 2012. 8(3): p. e1002438-e1002438. 61. O'Connell, M.N., et al., Dual mechanism of neuronal ensemble inhibition in primary auditory cortex. Neuron, 2011. 69(4): p. 805-17. 599 600 Schreiber, T., Measuring information transfer. Physical review letters, 62. Dotson, N.M., R.F. Salazar, and C.M. Gray, Frontoparietal correlation 601 dynamics reveal interplay between integration and segregation during 602 visual working memory. The Journal of Neuroscience, 2014. 34(41): p. 603 13600-13613. 604 63. 605 606 Salazar, R.F., et al., Content-specific fronto-parietal synchronization during visual working memory. Science, 2012. 338(6110): p. 1097-100. 64. Stetson, C. and R.A. Andersen, The parietal reach region selectively anti- 607 synchronizes with dorsal premotor cortex during planning. The Journal of 608 Neuroscience, 2014. 34(36): p. 11948-11958. 609 65. Nolte, G., et al., Robustly estimating the flow direction of information in 610 complex physical systems. Physical review letters, 2008. 100(23): p. 611 234101. 27 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity 612 66. Page nr. 28 of 38 Schoffelen, J.M., R. Oostenveld, and P. Fries, Neuronal coherence as a 613 mechanism of effective corticospinal interaction. Science, 2005. 308 614 (5718): p. 111-113. 615 67. 616 617 Baldauf, D. and R. Desimone, Neural mechanisms of object-based attention. Science, 2014. 344(6182): p. 424-7. 68. Maldonado, P.E., S. Friedman-Hill, and C.M. Gray, Dynamics of striate 618 cortical activity in the alert macaque: II. Fast time scale synchronization. 619 Cerebral Cortex, 2000. 10(11): p. 1117-1131. 620 69. Grothe, I., et al., Switching neuronal inputs by differential modulations of 621 gamma-band phase-coherence. The Journal of Neuroscience, 2012. 622 32(46): p. 16172-16180. 623 70. Singer, W. and C.M. Gray, Visual feature integration and the temporal 624 correlation hypothesis. Annual review of neuroscience, 1995. 18(1): p. 625 555-586. 626 627 628 28 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity 629 Boxes 630 Glossary Box Page nr. 29 of 38 631 632 Auto-correlogram: correlation across time between a signal and a time-shifted 633 copy of that same signal. 634 Coherence: degree to which the phase relations between two signals are 635 consistent across non-overlapping time windows. 636 Cross-correlogram: correlation across time between a signal and a time-shifted 637 copy of another signal. 638 Excitable phase: phase at which a neuron with an oscillating membrane 639 potential has the highest probability of generating an action potential in 640 response to excitatory synaptic input. 641 Decoding: data-analytic technique that estimates a variable (e.g., an animal’s 642 current position) on the basis of a set of other variables (e.g., neuronal signals). 643 Membrane potential: the difference in electric potential between the interior 644 and the exterior of a cell. The membrane potential is especially relevant for cells 645 with voltage-sensitive ion channels, such as neurons: the more positive the 646 membrane potential, the higher the probability that an excitatory synaptic input 647 results in an action potential (produced by a massive opening of voltage sensitive 648 sodium channels). 649 Neuronal oscillation: neuronal signal with a repeating waveform. Such a signal is 650 also said to exhibit periodicity or rhythmicity, and it is reflected in the fact that 651 the auto-correlogram has multiple peaks. 29 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 30 of 38 652 Phase-amplitude coupling: particular relation between a low- and a high- 653 frequency signal component, involving that the time-varying amplitude of the 654 high-frequency component predominantly occurs a particular phase of the low- 655 frequency component. 656 Phase relation: time relation between two periodic signal components, 657 expressed as the fraction of the repeating waveform by which the signals are 658 shifted. 659 Phase-relation diversity: variability across signal pairs with respect to their phase 660 relations. 661 Traveling wave: waveform over (as a function of) space that changes its position 662 over time. 663 Voltage-sensitive dye: dye that changes its spectral properties in response to 664 voltage changes, allowing it to be used for optical imaging of potential 665 distributions. 666 667 Box 1. Calculating phase relations, their diversity, and their coherence 668 In a first step, the phases of single-sites are obtained from a Fourier (or wavelet) 669 transform of some epoch of the recorded signal (see Box Fig. IA). The Fourier 670 coefficients are typically handled as real and imaginary components of 671 frequency-indexed complex numbers whose phase angle and magnitude are the 672 ingredients for several derived quantities. First, from the single-site phases 673 obtained from multiple sites, one can calculate the between-site phase relations 674 (see Box Fig. IB). In Box Figure IC, we consider an example situation involving 675 three sites (color-coded blue, green and red), resulting in two site pairs (blue30 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 31 of 38 676 green and green-red), for which we obtained recordings in a series of 677 trials/epochs. The phase relations are shown in the grey middle columns. They 678 are obtained as the product of the immediate left and the conjugate of the 679 immediate right column. This operation produces complex numbers whose 680 phase angles equal the phase differences (phase relations) between the two 681 sites, and magnitudes that equal the product between the magnitudes of the 682 single-site coefficients. Importantly, one must distinguish the epoch-specific 683 phase relations from their average across all epochs, which is also called the 684 average phase relation. Phase-relation diversity pertains to the diversity in the 685 average phase relations across the site pairs (here illustrated for 2 pairs – the 686 minimum number). Crucially, this across-site-pair diversity only makes sense if 687 the average phase relations are also reliable, implying that the corresponding 688 epoch-wise phase relations are similar to the average phase relation. This 689 reliability can be quantified using a split-half procedure (see Fig. 1A), as well as 690 by measures of coherence, which is obtained by normalizing the cross-spectrum 691 (the average of the complex numbers in the middle columns in Fig. IB) by the 692 square-roots of the sites’ power (the average of the squared magnitudes of the 693 numbers in the outer columns). Thus, phase-relation diversity is only relevant for 694 sites that are coherent, since otherwise the diverse average phase relations are 695 not representative of the epoch-wise phase relations. 696 697 Box 2. Phase-relation diversity across frequencies 698 The main focus of this review is on phase-relation diversity across space. 699 However, phase-relation diversity across frequencies also turns out to be an 31 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 32 of 38 700 important concept, and it is important to distinguish between the two. Diversity 701 across frequencies in the phase relations between two populations corresponds 702 to the conduction time delay between these populations, which may be positive 703 or negative. This time delay can be used to infer the directionality of the coupling 704 between the populations. Consider two populations that are coherently 705 oscillating in a range of frequencies, say between 15 and 30 Hz. Provided there is 706 a fixed time delay between these populations, the phase relations between 707 these populations will vary linearly with frequency (at least, within the frequency 708 band in which they are coherent). Critically, with increasing frequency, the phase 709 lags will increase when the flow of information is from population A to B (a 710 positive time delay), whereas they will decrease when the flow is in the opposite 711 direction (a negative time delay). As such, directionality of the coupling between 712 two oscillating populations can be inferred from the pattern of the across- 713 frequency phase-relation diversity. This rationale form the basis of the phase 714 slope index [65] that has proven to be highly informative in studying the roles of 715 neuronal oscillations in the flow of neuronal communication [43, 66, 67]. 716 717 Box 3. Phase-relation diversity and traveling waves 718 The defining feature of a traveling wave is that its phase varies monotonically as 719 a function of anatomical space – such as along the septo-temporal axis of the 720 hippocampus [31, 32]. As such, traveling waves are always characterized by 721 phase-relation diversity. However, the reverse is not necessarily true: phase- 722 relation diversity need not always reflect a traveling wave. For example, phase- 723 relation diversity may also emerge when (1) populations with different (non32 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 33 of 38 724 spatial) tuning profiles (e.g., V1 populations that differ with respect to their 725 orientation tuning), (2) nearby cortical areas (e.g., V1, V2, and V3) or (3) different 726 layers, oscillate at different phases. In all these cases, phase relations will be 727 diverse across space, but will not necessarily vary monotonically as a function of 728 anatomical space. Importantly, in such cases, oscillatory phase may still vary 729 systematically as a function of a latent coding space (e.g., orientation tuning) and 730 may thus serve important functional roles in neuronal information processing 731 (see Phase-relation diversity and the flow of neuronal information). Thus, if the 732 underlying connectivity is mostly patchy rather than continuous, then wave-like 733 propagation across a latent coding space will produce non-wave-like phase 734 diversity across anatomical space. 735 736 Box 4. Coherence in gamma-band oscillations over visual cortex 737 Gamma-band oscillations in visual cortex are strongly dependent on the 738 presence of a visual stimulus in the neuronal population’s receptive field [RF; see 739 57, supplementary material, 68]. Gamma-band oscillations are coherent 740 between neuronal populations with overlapping RFs: two studies showed this for 741 populations in V1 and V4 [57, 69], and one study for populations in V4 and FEF 742 [43]. The situation without RF overlap (but with visual input in both RFs) is 743 equally relevant for the proposed selection and segregation mechanisms (see 744 Fig. 3B and 3C), but for this configuration, gamma-band coherence is less well 745 established. Specifically, early studies in cat visual cortex [reviewed in 70] 746 demonstrated gamma-band coherence between the spiking activity of neurons 747 with non-overlapping receptive fields (RFs), but also showed that this coherence 33 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 34 of 38 748 was more likely when the RFs overlapped. This pattern was replicated in 749 macaque V1 [68]. The latter study also showed that gamma-band coherence is a 750 rare phenomenon overall (8% of the single unit, and 15% of the multi-unit pairs) 751 but quite common among the single/multi-unit pairs whose firing is synchronized 752 (60% of the single unit, and 69% of the multi-unit pairs). Further, in a recent 753 study in macaques, no gamma-band coherence was found between FEF multi- 754 units and V4 LFPs when the underlying populations had non-overlapping RFs 755 [43]. However, it is unclear whether RF overlap is the essential factor that 756 determines whether gamma-band coherence is and is not observed. In fact, in 757 one of the studies we discussed, coherent gamma-band oscillation were shown 758 to exhibit non-zero phase relations that had a Granger-causal effect on the firing 759 rate [58]. Unfortunately, this study did not document RF overlap. Thus, part of 760 the results show that gamma-band coherence is limited to neuronal populations 761 that are communicating, and suggest that selectivity is implemented through the 762 selective presence of coherence. Under which conditions gamma-band 763 oscillations can also implement selectivity through phase-relation diversity 764 requires more empirical work. 765 766 767 34 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 35 of 38 768 Figure Legends 769 Figure 1. Reliable phase-relation diversity occurs during sensory processing, in 770 high (>20Hz) frequency ranges, and in spatially distributed phase amplitude 771 coupling. (A) Schematic of reliable phase-relation diversity. Left panel: 772 hypothetical phase distribution across a patch of cortex, together with nine 773 hypothetical recording sites. Polar plots (see also Box 1) represent the phases in 774 each site relative to an arbitrary reference site and depict diversity in phase 775 relations across space. Right panel: a split-half procedure can be used to 776 separate systematic from unsystematic phase-relation diversity. Phase-relation 777 diversity discussed in this review refers to phase relations that are both diverse 778 (spread out across the x and y axis in the plot) as well as reliable (similar in the 779 two random partitions of the data and therefore close to the diagonal in the 780 plot). Dots represent individual site-pairs. (B) Empirical data demonstrating 781 reliable phase-relation diversity of gamma oscillations during visual stimulation 782 in monkey LFP recordings in V4 (left panel) as well as in human MEG recordings 783 from posterior sites (right panel). Scatterplots show the outcomes of the split- 784 half approach explained in A and form the basis for the depicted frequency 785 spectra of reliable diversity. (C) Phase-relation diversity is also prominent in 786 patterns of spatially-distributed phase-amplitude coupling (PAC). A 787 decomposition technique was used to extract spatial and spectral profiles of 788 distributed PAC components. For one representative component, the phase- 789 providing spatial map (depicted to the right) shows diversity in the preferred 790 coupling phase across space. In the lower panels, we zoomed in on two 791 representative site pairs. Whereas site 43 shows similar coupling strength with 35 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 36 of 38 792 both site 29 and 57, the high frequency oscillations in site 43 are nested at 793 opposite phases of the low-frequency oscillations in sites 29 (left) and 57 (right). 794 Panel B adapted from [27] and [28]; panel C adapted from [33]. 795 796 Figure 2. Phase-relation diversity is modulated by sensorimotor events, as well 797 as by the strength of the excitatory drive. (A) Beta oscillations following a 798 saccade are characterized by phase-relation diversity in the form of a traveling 799 wave from sites representing foveal visual input to sites representing more 800 peripheral input. (B) Upper plot: phase-relation diversity of hippocampal theta 801 oscillations is dynamic during spatial navigation. Lower plot: reconstructions of 802 the rat’s location based on spikes and LFP. Phase-relation diversity in the LFP 803 allows for a reconstruction (decoding) of the rat’s location and this 804 reconstruction is approximately as accurate as the reconstruction that is based 805 on spikes. (C) Modelling study demonstrating that the excitatory drive can 806 modulate phase-relation diversity across a neuronal population. Panel A adapted 807 from Zanos, Mineault [35]; panel B adapted from Agarwal, Stevenson [36]; panel 808 C adapted from Lowet, Roberts [53]. 809 810 Figure 3. Schematic of how phase-relation diversity may contribute to the flow 811 of neuronal information. For this schematic we make two key assumptions. First, 812 oscillatory phase (here represented in polar plots; see also box 1) reflects 813 excitability fluctuations of some neuronal population. As a result, the 814 effectiveness of the communication between a lower (sending) and a higher 815 (receiving) area will be determined by their phase relation. For simplicity, we 36 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 37 of 38 816 assume a conduction delay of 0 ms, such that communication is most effective 817 when the lower- and the higher-area populations have a phase relation of 0 (i.e. 818 when output that is concentrated at the most excitable phase of the sending 819 population will arrive in the most excitable phase of the receiving population). In 820 reality, the most effective phase relation will depend on the precise conduction 821 delay. Second, all depicted sites are equally coherent with each other; only the 822 average phase relations differ between the three depicted scenarios. (A) 823 Scenario without diversity in the phase relations across the lower-area 824 populations. As a consequence, output from the three depicted lower-area 825 populations will arrive in the same phase of the oscillation in the higher-area 826 population. Because coherence is assumed to be equal for all three populations, 827 their impact is equally effective. The graph on the right-hand side shows the 828 impact onto the higher-area population (Y-axis) as function of the location of the 829 lower-area population (X-axis). The flat line in this graph indicates that 830 communication is not selective. (B) Scenario in which diversity is introduced 831 across the lower-area populations. As a consequence, across the lower-area, 832 some populations will have more effective phase relations with the higher-area 833 population than others. This allows for selective communication between the 834 lower- and higher-area populations, as depicted in the right panel. This 835 selectivity can be modulated by changing the degree and the spatial extent of 836 the phase-relation diversity. (C) When diversity is also introduced across the 837 higher-area receiving populations, this may further allow for the concurrent 838 segregation of multiple information streams. This is depicted in the right-hand 839 panel, in which the colours of the lines denote the three higher-area populations. 37 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Phase-relation diversity Page nr. 38 of 38 840 841 842 Figure I. Calculating phase relations, their diversity, and their coherence. See 843 Box 1 text. 844 38 Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Outstanding Questions What are the precise conditions (in terms of anatomical connectivity motif, network input, …) under which phase-relation diversity is relevant for effective inter-areal communication? Do the modulations of phase-relation diversity by sensory and motor events imply modulations of the flow of neuronal information? How is phase-relation diversity generated and/or controlled at the level of the underlying neuronal network? What is the relation between oscillatory amplitude and phase-relation diversity with respect to (1) how they are modulated by task- and stimulus variables, and (2) how these parameters are determined by the underlying neuronal network? Does phase-relation diversity observed in macroscopic EEG/MEG measurements reflect the same neurophysiological processes as phase-relation diversity in local LFP recordings? How does phase-relation diversity change with cohort variables, such as ageing and disease, and could it serve as a useful biomarker for changes in neuronal information processing? Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 Trends box Neuronal oscillations are a prominent feature of neuronal population activity. Across populations, oscillations exhibit stable phase-relations that can be highly diverse across nearby sites. Phase-relation diversity is prominent even for high frequency (>40 Hz) oscillations during sustained visual stimulation, and is therefore neither restricted to low frequencies nor to periods outside of sensory processing. Arguing for its computation relevance, this diversity is modulated by both sensory and motor events and, for hippocampal theta oscillations, it even allows for reconstructing a rat’s location. Phase-relation diversity may enhance selectivity of neuronal communication and allow for the concurrent segregation of multiple information streams. This may be particularly relevant for beta oscillations that are coherent across frontal and parietal brain areas. A Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99magnitude Fourier doi: 10.1016/j.neuron.2015.12.018 Transform phase B real site 2 ... ... site 2 phase relation site 3 ... phase relation ... site 1 ... trials / epochs C ... phase of site x phase relation phase of site y between site x and site y imaginary average Σ power crossspectrum normalize coherence [1,2] power Σ power phase-relation diversity [1,2] vs [2,3] crossspectrum normalize coherence [2,3] power 1 beta LFP A Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 foveal foveal peripheral peripheral 0 -1 0 time rel. to saccade end (ms) B theta phase hippocamp. sites coupling strength C 40 20 rat’s location phase-locking 0.2 0.7 0 -0.02 1 0 ∆ frequency 10 −π/2 -10 0.02 -0.02 0 ∆ phase 0.02 -0.02 excitatory drive 0 π/2 0.02 impact higher area lower area no local diversity impact on higher area A Published in final edited form in: Neurosciences 39 (2), February 2016, Pages 86-99 doi: 10.1016/j.neuron.2015.12.018 location in lower area local diversity impact on higher area B C local diversity local diversity impact on higher area location in lower area location in lower area
© Copyright 2026 Paperzz