BioMed Research International

BioMed Research International / 2015 / Article
Special Issue

From Micro- to Macroscopic Brain Connectivity Using Multiple Modalities

View this Special Issue

Review Article | Open Access

Volume 2015 |Article ID 915606 | 12 pages | https://doi.org/10.1155/2015/915606

Modeling the Generation of Phase-Amplitude Coupling in Cortical Circuits: From Detailed Networks to Neural Mass Models

Academic Editor: Vincenzo Romei
Received27 Mar 2015
Revised28 Jul 2015
Accepted06 Aug 2015
Published11 Oct 2015

Abstract

Phase-amplitude coupling (PAC), the phenomenon where the amplitude of a high frequency oscillation is modulated by the phase of a lower frequency oscillation, is attracting an increasing interest in the neuroscience community due to its potential relevance for understanding healthy and pathological information processing in the brain. PAC is a diverse phenomenon, having been experimentally detected in at least ten combinations of rhythms: delta-theta, delta-alpha, delta-beta, delta-gamma, theta-alpha, theta-beta, theta-gamma, alpha-beta, alpha-gamma, and beta-gamma. However, a complete understanding of the biophysical mechanisms generating this diversity is lacking. Here we review computational models of PAC generation that range from detailed models of neuronal networks, where each cell is described by Hodgkin-Huxley-type equations, to neural mass models (NMMs) where only the average activities of neuronal populations are considered. We argue that NMMs are an appropriate mathematical framework (due to the small number of parameters and variables involved and the richness of the dynamics they can generate) to study the PAC phenomenon.

1. Introduction

From the theory of signal processing we know that if an input-state-output system is linear its output will have the same frequency content as its inputs. Conversely, in nonlinear systems, the energy at one frequency in the inputs appears at different frequencies in the outputs. This induces cross-frequency coupling (CFC) between any two sources, when the output of one serves as the input to the other [1]. It has been shown that pyramidal cells produce a varied set of intrinsic dynamics based only on the type and compartmental localization of intrinsic conductances [2]. A combination of sodium, potassium, and calcium conductances produces coexistent gamma (~40 Hz) and theta (~6 Hz) rhythms on tonic depolarization. In contrast, combinations of persistent sodium and potassium channels in the soma produce a use-dependent transition between regular spiking at ~10 Hz and a repetitive, brief burst generation at ~20 Hz [2]. Since cortical columns and brain areas generating different brain rhythms are interconnected, the presence of CFC should not be surprising, even if the exact mechanisms responsible for its generation remain imprecise. The question is then whether CFC is only a mechanistic result of the way the brain is constructed or if it also has a role in brain computations. At least six types of CFC have been documented [3, 4]: amplitude-amplitude coupling (AAC), phase-phase coupling (PPC), frequency-frequency coupling (FFC), phase-amplitude coupling (PAC), phase-frequency coupling (PFC), and frequency-amplitude coupling (FAC). PAC, the type of CFC that occurs when the phase of a low frequency oscillation modulates the amplitude of a higher frequency oscillation, has received a lot of attention in the last decade due to its potential relevance for understanding healthy and pathological brain function [511]. PAC has been hypothesized to be the carrier mechanism for the interaction of local and global processes and therefore being directly linked to the integration of distributed information in the brain [12]. For instance, it has been suggested that theta-gamma PAC is used as a coding scheme for multi-item short-term memory in the hippocampus, where different spatial information is represented in different gamma subcycles of a theta cycle [13, 14]. Recent experimental evidence also suggests that PAC links local blood-oxygen-level-dependent (BOLD) signals to BOLD correlation across distributed networks [15].

In parallel to the experimental study of the PAC phenomenon, computational models have been proposed in order to clarify the neurobiological mechanism underlying its generation [1623]. Here we review these models, going from the detailed description of each cell (via the Hodgkin-Huxley formalism) in neuronal networks to neural mass models (NMMs), which are a type of mean field description that focuses on the dynamics of the average activity in a neuronal population while neglecting the second-order statistics (variance and covariances) and from models only focusing on generation of the theta-gamma PAC in the hippocampus to the most recent models capable of simultaneously generating several PAC combinations.

This review is structured as follows. First, in Section 2, we show that there is evidence for at least ten different PAC combinations (of a low and a higher frequency oscillation). Computational models of the PAC phenomenon can be divided into two types: detailed and NMMs. The main characteristics of these two types are briefly discussed in Section 3, followed by two sections describing specific models of both types.

2. Experimental Evidence of the Diversity of the PAC Phenomenon

The classic example of PAC was demonstrated in the CA1 region of the hippocampus [24] where the phase of the theta rhythm was found to modulate the power of gamma oscillations. Later studies found that PAC is neither restricted to theta-gamma coupling nor to the hippocampus. For instance, PAC has also been reported in the frontal, posterior, and parietal human cortices during auditory, visual, linguistic, and memory tasks [2527], in monkey auditory and visual cortices [15, 28, 29] and rodent olfactory bulb [30]. In addition to Bragin et al.’s study [24], other studies have confirmed the existence of theta-gamma coupling in the hippocampus [3134] and other brain areas [3544]. Other PAC combinations of low and high frequency rhythms have also been detected: delta-theta [37, 45], delta-alpha [46, 47], delta-beta [44, 46], delta-gamma [34, 35, 38, 41, 44], theta-alpha [46], theta-beta [44, 46], alpha-beta [45], alpha-gamma [15, 26, 27, 35, 46, 48, 49], and beta-gamma [7, 15].

It should be noted that the studies mentioned above do not always use the same frequency values for the boundaries of the different brain rhythms [50] and that sometimes the gamma band is divided into different subbands such as low-gamma, middle-gamma, and fast-gamma, with boundaries that can differ between different studies. Thus, subdivisions of classical bands can potentially increase the number of PAC combinations to be studied. Additionally, a high number of mathematical methods for detecting PAC have been proposed [3, 12, 5157], each with advantages and caveats, and no gold standard has emerged. Furthermore, those methods are not exempted of spurious results, that is, identifying PAC that is not related to true modulations between neuronal subsystems. These issues (reviewed recently in [58]) are out of the scope of this review, but we mention them here to highlight the fact that the experimental study of the PAC phenomenon is far from being complete and new methods and results in the upcoming years will be necessary to complement, inform, and refine past and future computational models of the phenomenon.

3. Detailed Mathematical Models versus Neural Mass Models

There are two main approaches to modeling the dynamics of neuronal populations. One approach is to realistically model each cell in the network, using multiple compartments for the soma, axon, and dendrites. The most prominent example of this approach is the Blue Brain Project [59], which aims to achieve in the next decade a full simulation of human brain dynamics (a network of approximately 86 billion neurons) in a supercomputer. A practical disadvantage of such realistic modeling is that it requires high computational power. For this reason, simplified versions of such models in which only one compartment is taken into account have been used [16, 60]. However, even in this case, the use of such detailed models makes it difficult to determine the influence of each model parameter on the generated average network characteristics. The second approach is based on the use of NMMs, which constitute a class of biophysical models that captures the average activity of neuronal ensembles, rather than modeling each neuron in the network individually [61, 62]. NMMs are described by nonlinear differential equations and can be rigorously obtained from mean field approaches [6365] after neglecting the second-order moments. For instance, the Wilson-Cowan neural mass model [61] can be obtained from a mean field approximation of two coupled networks of FitzHugh-Nagumo neurons [63]. An alternative way of constructing the NMM formalism is to consider that each neuronal population performs two mathematical operations [62]. The first is the conversion of postsynaptic potentials (PSP) into an average density of action potentials (AP) which is modeled using a sigmoid function. The second operation is the conversion of AP into PSP, which is done by means of a linear convolution with an impulse response function. The Wilson-Cowan model is obtained when the impulse response function is , which produces a system of first-order differential equations describing the activity in each population. A more recent neural mass model, the Jansen-Rit model [62], is obtained when the impulse response function has the form . This results in a system of second-order differential equations describing the dynamics of PSPs in each population. Computational models based on Wilson-Cowan and Jansen-Rit models have provided the mathematical framework for simulating the generation of electrical activity in the brain during resting state [62, 6671], stimulation [62, 7274], and disease [67, 7577].

4. Detailed Mathematical Models

Detailed mathematical models of PAC generation [16, 18] have focused on the theta-gamma interaction observed in the hippocampus [24]. These models consist of either purely inhibitory networks [16] or networks with both excitatory and inhibitory cells [1820] and are based on models previously developed to study the generation of theta and gamma rhythms separately [23].

4.1. Inhibitory-Inhibitory () Network

A simulated inhibitory network in the hippocampus containing fast and slow interneurons was shown to generate theta-gamma coupling under restricted conditions [16]. The network comprised single compartment neurons modeled with the Hodgkin-Huxley formalism:where index , counts the cells in the network, is the applied current, and is a normally distributed noise. The sodium , potassium , leak , and synaptic currents are given byThe cell population was divided into half on the basis of fast and slow synaptic dynamics. Synaptic conductances had one of four possible values depending on the types of the cells connected: fast cell to fast cell, fast cell to slow cell, slow cell to slow cell, and slow cell to fast cell. Connectivity was all to all. Equations for the gating variables , , and , as well as parameter values can be found in [16]. The numerical simulations performed in [16] showed that the model can generate PAC under restricted conditions that included strong connections within the same populations, weaker connections between populations (especially from fast to slow populations), and carefully tuned inputs.

4.2. Excitatory-Inhibitory () Networks

Hippocampal networks producing theta-gamma PAC also have pyramidal cells. To consider this situation, a model comprising three neuronal populations was proposed in [18] and was shown to produce theta-gamma PAC [23]. The three populations are pyramidal cells, fast-spiking basket cells, and oriens lacunosum-moleculare (O-LM) interneurons. The outputs of the O-LM cells are projected as slow inhibitory postsynaptic potentials (IPSP) onto the distal apical dendrites of pyramidal cells [18].

Basket cells were modeled with a single compartment, using the fast-spiking interneuron model proposed in [78], similar to (1) and (2). O-LM cells were also modeled with a single compartment. In addition to sodium, potassium, leak, and synaptic current, two other currents were considered: the h-current and the A current [17, 18]. Pyramidal cells were modeled by 5 compartments: 1 for basal dendrites, 1 for soma, and 3 for apical dendrites. The equation for each compartment () iswhere is the current due to electrical coupling between compartments. The expressions for the ionic and synaptic currents as well as the parameter values to simulate the model can be found in the supplementary information section in [18]. Different simulations were performed in [18], but the one with the highest number of cells comprised a total of 180 cells. Their results showed that O-LM cells alone can coordinate cell assemblies and that the same theta rhythm can coordinate different cell assemblies with different frequencies in the gamma range [18, 23].

5. Neural Mass Models

In this section we review three NMMs that are able to generate PAC. The first two studies [21, 22] are based on the works of Wilson-Cowan and Jansen-Rit and only focus on the generation of one PAC combination. The last study [79] is also based on the Jansen-Rit model but is able to simultaneously generate different PAC combinations.

5.1. PAC Generation Using the Wilson-Cowan Model

Onslow et al. [21] used the Wilson-Cowan model to study the generation of theta-gamma PAC in a brain region not necessarily restricted to the hippocampus. The model comprises two coupled populations (Figure 1(a)), one excitatory and one inhibitory. The system of first-order differential equations describing the model iswhere and are the average activity levels of excitatory and inhibitory populations, respectively [61] and and are the external inputs to the two populations. The weight of the connection from the excitatory population to the inhibitory population is and from the inhibitory to the excitatory population is , and the self-connections are and . and are the time constants for each population. The nonlinearity in the model is introduced by means of a sigmoid function:where parameter determines the steepness of the sigmoid curve, determines the position of the sigmoid function, and determines the amplitude of the response.

System (4) is capable of producing oscillations due to the reciprocal connections between the two populations. Numerical simulations showed [21] that this model generates gamma oscillations that are locked to a certain phase of theta oscillations when considering oscillatory inputs.

Figure 2(a) shows a realization of the model where the phase of a 4 Hz oscillation modulates the amplitude of a 55 Hz oscillation. The parameter values used in this simulation were , , , , , , , , , , and . Additional simulations showed [21] that the amplitude, frequency, and phase-locking characteristics of the PAC activity generated were dependent on the strength of the connectivity parameters and on the amplitude and mean value of the low frequency input signal. It was possible to tune the parameters of the model to produce different frequencies of activity phase-locked to different phases of the theta rhythm [21].


Parameter (units)InterpretationValue

(mV)Gain, , , , , , , , , , , , , and

(s−1)Reciprocal of time constant, , , , , , , , , ,  , , , and

External input for , , and

Damping coefficient for all populations

(s−1)Maximum firing rate for all populations

(mV)Position of the sigmoid function for all populations

(mV−1)Steepness of the sigmoid function for all populations


L2/3L4L5L6
RSIBLTSFSRSLTSFSRSIBLTSFSRSLTSFS

L2/3RS25101015025300000000
IB1025550000000000
LTS00000000
FS00000000

L4RS10001530251800000
LTS0000000000
FS0000000000

L5RS000000000
IB000000000
LTS0000000
FS0000000

L6RS0000450100000
LTS00000000000
FS00000000000

5.2. The Jansen-Rit Model of a Cortical Column

The Jansen-Rit model of a cortical column [62] comprises three neuronal populations (Figure 1(b)): pyramidal cells, excitatory interneurons, and inhibitory interneurons. The model is mathematically described by a system of second-order differential equations:where is the excitatory postsynaptic potential (EPSP) that feeds into the two populations of interneurons and and are EPSP and inhibitory postsynaptic potentials (IPSP) that enter into the pyramidal cell population, respectively. , , , and are the connection strengths between the populations. In this model, the electroencephalography (EEG) signal is considered to be proportional to .

Figure 2(b) shows a realization of model (6) where delta (3 Hz)-alpha (11 Hz) PAC is produced when considering an oscillatory input . The parameter values used in this simulation were , , , , , , , , , , , and .

Alternatively, EEG signals presenting PAC can be obtained by coupling multiple Jansen-Rit models (see Figures 5, 8, and 9 in [67]). In a more recent work [22], several Jansen-Rit models were also coupled and the cross-frequency transfer was studied in a setting where oscillators (generating the different rhythms) were coupled unidirectionally and thus the driving between them was passive. This study showed that this passive driving can also account for CFC in the brain as a result of the complex nonlinear dynamics of the underlying neuronal activity.

5.3. Cortical Column Model Comprising 4 Layers and 14 Neuronal Populations

A more complex neural mass model of the cortical column was recently proposed [79] in which 4 cortical layers and 14 neuronal populations are considered. Figure 1(c) shows the model obtained by distributing four cell classes in four cortical layers (L2/3, L4, L5, and L6). This produced 14 different neuronal populations, since not all cell classes are present in every layer [80]. Excitatory neurons were either regular spiking (RS) or intrinsically bursting (IB) ones, and inhibitory neurons were either fast-spiking (FS) or low-threshold spiking (LTS) neurons.

The model is based on the Jansen-Rit model and the dynamics of the average PSP in each neuronal population is obtained by solving a system of 14 second-order differential equations:where , and . The populations are numbered from 1 to 14 following the order: [L2RS, L2IB, L2LTS, L2FS, L4RS, L4LTS, L4FS, L5RS, L5IB, L5LTS, L5FS, L6RS, L6LTS, L6FS]. Layer 2/3 was labelled as 2. As can be seen in (7), neuronal populations interact via the connectivity matrix (Figure 1(d)). External inputs are accounted for via which can be any arbitrary function including white noise [62]. The “damping” parameter critically determines the behavior of the system. For (which corresponds to the Jansen-Rit model) an individual population is not capable of oscillating, and it is the presence of interpopulation connections (nonzero , ) that produces oscillatory behavior that mimics observed EEG signals. To account for the possibility of an oscillatory population [78, 81] a nonzero value for was used.

Figure 2(c) presents the temporal evolution of the average PSP in each neuronal population. Time series colored in red correspond to excitatory PSP (EPSP) whereas inhibitory PSP (IPSP) are presented in blue. As seen in the figure, both EPSP and IPSP present the characteristic “waxing and waning” (i.e., amplitude modulation) observed in real EEG signals. Parameters values are presented in Tables 1 and 2. To quantify the PAC phenomenon, a causality measure between time series, the information flow [82], was computed using phases and amplitudes of the signals shown in Figure 2(c). Figure 3 shows the information flow from the phase to the amplitude for nine different combinations of phases and amplitudes: delta-theta, delta-alpha, delta-beta, delta-gamma, theta-alpha, theta-beta, theta-gamma, alpha-beta, and alpha-gamma. A negative value of the information flow means that the phase tends to stabilize the amplitude whereas a positive value means that the phase tends to make the amplitude more uncertain. An exploratory analysis of the influence of the parameters on PAC showed that changes in external input and time constants produced theta-gamma PAC values higher than alpha-gamma PAC, whereas changes in connectivity produced higher alpha-gamma PAC values. Additional information can be found in [79].

6. Conclusions

In conclusion, we have shown that PAC is a diverse phenomenon, not restricted to the theta-gamma coupling in the hippocampus. In order to model the complexity of the PAC phenomenon, which is hypothesized to bridge local and global scales in the brain [12, 15], reduced models of neuronal activity such as NMMs are needed, since detailed models are computationally expensive and their results are difficult to interpret due to the high number of variables and parameters involved. An open problem to be explored with NMMs is how the different PAC combinations are related.

While both types of models reviewed here, detailed models and NMMs, are capable of generating signals reflecting PAC, only in a few studies a quantitative measure of the phenomenon has been provided. This is probably related to the lack of a gold standard for PAC detection, which has resulted in the development of numerous methods.

The computational models summarized here focused on the mechanistic generation of the PAC phenomenon. NMMs are simple (in the sense of the few variables and parameters involved) but complex (in the sense of the richness of the dynamics they can generate) enough to approach important questions related to the functional role of the PAC phenomenon.

Conflict of Interests

The author declares that no conflict of interests exists.

References

  1. C. C. Chen, S. J. Kiebel, and K. J. Friston, “Dynamic causal modelling of induced responses,” NeuroImage, vol. 41, no. 4, pp. 1293–1312, 2008. View at: Publisher Site | Google Scholar
  2. R. D. Traub, E. H. Buhl, T. Gloveli, and M. A. Whittington, “Fast rhythmic bursting can be induced in layer 2/3 cortical neurons by enhancing persistent Na+ conductance or by blocking BK channels,” Journal of Neurophysiology, vol. 89, no. 2, pp. 909–921, 2003. View at: Publisher Site | Google Scholar
  3. V. Jirsa and V. Müller, “Cross-frequency coupling in real and virtual brain networks,” Frontiers in Computational Neuroscience, 2013. View at: Publisher Site | Google Scholar
  4. O. Jensen and L. L. Colgin, “Cross-frequency coupling between neuronal oscillations,” Trends in Cognitive Sciences, vol. 11, no. 7, pp. 267–269, 2007. View at: Publisher Site | Google Scholar
  5. J. López-Azcárate, M. Tainta, M. C. Rodríguez-Oroz et al., “Coupling between beta and high-frequency activity in the human subthalamic nucleus may be a pathophysiological mechanism in Parkinson's disease,” Journal of Neuroscience, vol. 30, no. 19, pp. 6667–6677, 2010. View at: Publisher Site | Google Scholar
  6. S. A. Shimamoto, E. S. Ryapolova-Webb, J. L. Ostrem, N. B. Galifianakis, K. J. Miller, and P. A. Starr, “Subthalamic nucleus neurons are synchronized to primary motor cortex local field potentials in Parkinson's disease,” Journal of Neuroscience, vol. 33, no. 17, pp. 7220–7233, 2013. View at: Publisher Site | Google Scholar
  7. C. de Hemptinne, E. S. Ryapolova-Webb, E. L. Air et al., “Exaggerated phase-amplitude coupling in the primary motor cortex in Parkinson disease,” Proceedings of the National Academy of Sciences of the United States of America, vol. 110, no. 12, pp. 4780–4785, 2013. View at: Publisher Site | Google Scholar
  8. E. A. Allen, J. Liu, K. A. Kiehl et al., “Components of cross-frequency modulation in health and disease,” Frontiers in Systems Neuroscience, vol. 5, article 59, 2011. View at: Publisher Site | Google Scholar
  9. L. V. Moran and L. E. Hong, “High vs low frequency neural oscillations in schizophrenia,” Schizophrenia Bulletin, vol. 37, no. 4, pp. 659–663, 2011. View at: Publisher Site | Google Scholar
  10. K. Kirihara, A. J. Rissling, N. R. Swerdlow, D. L. Braff, and G. A. Light, “Hierarchical organization of gamma and theta oscillatory dynamics in schizophrenia,” Biological Psychiatry, vol. 71, no. 10, pp. 873–880, 2012. View at: Publisher Site | Google Scholar
  11. V. Miskovic, D. A. Moscovitch, D. L. Santesso, R. E. McCabe, M. M. Antony, and L. A. Schmidt, “Changes in EEG cross-frequency coupling during cognitive behavioral therapy for social anxiety disorder,” Psychological Science, vol. 22, no. 4, pp. 507–516, 2011. View at: Publisher Site | Google Scholar
  12. R. T. Canolty and R. T. Knight, “The functional role of cross-frequency coupling,” Trends in Cognitive Sciences, vol. 14, no. 11, pp. 506–515, 2010. View at: Publisher Site | Google Scholar
  13. J. E. Lisman and M. A. P. Idiart, “Storage of 7+/-2 short-term memories in oscillatory subcycles,” Science, vol. 267, no. 5203, pp. 1512–1515, 1995. View at: Google Scholar
  14. J. E. Lisman and O. Jensen, “The theta-gamma neural code,” Neuron, vol. 77, no. 6, pp. 1002–1016, 2013. View at: Publisher Site | Google Scholar
  15. L. Wang, Y. B. Saalmann, M. A. Pinsk, M. J. Arcaro, and S. Kastner, “Electrophysiological low-frequency coherence and cross-frequency coupling contribute to BOLD connectivity,” Neuron, vol. 76, no. 5, pp. 1010–1020, 2012. View at: Publisher Site | Google Scholar
  16. J. A. White, M. I. Banks, R. A. Pearce, and N. J. Kopell, “Networks of interneurons with fast and slow γ-aminobutyric acid type A (GABAA) kinetics provide substrate for mixed gamma-theta rhythm,” Proceedings of the National Academy of Sciences of the United States of America, vol. 97, no. 14, pp. 8128–8133, 2000. View at: Publisher Site | Google Scholar
  17. H. G. Rotstein, D. D. Pervouchine, C. D. Acker et al., “Slow and fast inhibition and an H-current interact to create a theta rhythm in a model of CA1 interneuron network,” Journal of Neurophysiology, vol. 94, no. 2, pp. 1509–1518, 2005. View at: Publisher Site | Google Scholar
  18. A. B. L. Tort, H. G. Rotstein, T. Dugladze, T. Gloveli, and N. J. Kopell, “On the formation of gamma-coherent cell assemblies by oriens lacunosum-moleculare interneurons in the hippocampus,” Proceedings of the National Academy of Sciences of the United States of America, vol. 104, no. 33, pp. 13490–13495, 2007. View at: Publisher Site | Google Scholar
  19. X. Zhang, K. M. Kendrick, H. Zhou, Y. Zhan, and J. Feng, “A computational study on altered theta-gamma coupling during learning and phase coding,” PLoS ONE, vol. 7, no. 6, Article ID e36472, 2012. View at: Publisher Site | Google Scholar
  20. E. Spaak, M. Zeitler, and S. Gielen, “Hippocampal theta modulation of neocortical spike times and gamma rhythm: a biophysical model study,” PLoS ONE, vol. 7, no. 10, Article ID e45688, 2012. View at: Publisher Site | Google Scholar
  21. A. C. Onslow, M. W. Jones, R. Bogacz, and A. B. Tort, “A canonical circuit for generating phase-amplitude coupling,” PLoS ONE, vol. 9, no. 8, Article ID e102591, 2014. View at: Publisher Site | Google Scholar
  22. M. Jedynak, A. J. Pons, and J. Garcia-Ojalvo, “Cross-frequency transfer in a stochastically driven mesoscopic neuronal model,” Frontiers in Computational Neuroscience, vol. 9, article 14, 2015. View at: Publisher Site | Google Scholar
  23. N. Kopell, C. Börgers, D. Pervouchine, P. Malerba, and A. Tort, “Gamma and theta rhythms in biophysical models of hippocampal circuits,” in Hippocampal Microcircuits: A Computational Modeler's Resource Book, V. Cutsuridis, B. P. Graham, S. Cobb, and I. Vida, Eds., Springer, New York, NY, USA, 2010. View at: Google Scholar
  24. A. Bragin, G. Jandó, Z. Nádasdy, J. Hetke, K. Wise, and G. Buzsáki, “Gamma (40–100 Hz) oscillation in the hippocampus of the behaving rat,” The Journal of Neuroscience, vol. 15, no. 1, pp. 47–60, 1995. View at: Google Scholar
  25. R. T. Canolty, E. Edwards, S. S. Dalal et al., “High gamma power is phase-locked to theta oscillations in human neocortex,” Science, vol. 313, no. 5793, pp. 1626–1628, 2006. View at: Publisher Site | Google Scholar
  26. D. Osipova, D. Hermes, and O. Jensen, “Gamma power is phase-locked to posterior alpha activity,” PLoS ONE, vol. 3, no. 12, Article ID e3990, 2008. View at: Publisher Site | Google Scholar
  27. B. Voytek, R. T. Canolty, A. Shestyuk, N. E. Crone, J. Parvizi, and R. T. Knight, “Shifts in gamma phase-amplitude coupling frequency from theta to alpha over posterior cortex during visual tasks,” Frontiers in Human Neuroscience, vol. 4, article 191, 2010. View at: Publisher Site | Google Scholar
  28. P. Lakatos, C.-M. Chen, M. N. O'Connell, A. Mills, and C. E. Schroeder, “Neuronal oscillations and multisensory interaction in primary auditory cortex,” Neuron, vol. 53, no. 2, pp. 279–292, 2007. View at: Publisher Site | Google Scholar
  29. P. Lakatos, G. Karmos, A. D. Mehta, I. Ulbert, and C. E. Schroeder, “Entrainment of neuronal oscillations as a mechanism of attentional selection,” Science, vol. 320, no. 5872, pp. 110–113, 2008. View at: Publisher Site | Google Scholar
  30. D. Rojas-Líbano and L. M. Kay, “Olfactory system gamma oscillations: the physiological dissection of a cognitive neural system,” Cognitive Neurodynamics, vol. 2, no. 3, pp. 179–194, 2008. View at: Publisher Site | Google Scholar
  31. R. Scheffer-Teixeira, H. Belchior, F. V. Caixeta, B. C. Souza, S. Ribeiro, and A. B. L. Tort, “Theta phase modulates multiple layer-specific oscillations in the CA1 region,” Cerebral Cortex, vol. 22, no. 10, pp. 2404–2414, 2012. View at: Publisher Site | Google Scholar
  32. M. J. Gillies, R. D. Traub, F. E. N. LeBeau et al., “A model of atropine-resistant theta oscillations in rat hippocampal area CA1,” Journal of Physiology, vol. 543, no. 3, pp. 779–793, 2002. View at: Publisher Site | Google Scholar
  33. B. Lega, J. Burke, J. Jacobs, and M. J. Kahana, “Slow-theta-to-gamma phase-amplitude coupling in human hippocampus supports the formation of new episodic memories,” Cerebral Cortex, 2014. View at: Publisher Site | Google Scholar
  34. J. Gross, N. Hoogenboom, G. Thut et al., “Speech rhythms and multiplexed oscillatory sensory coding in the human brain,” PLoS Biology, vol. 11, no. 12, Article ID e1001752, 2013. View at: Publisher Site | Google Scholar
  35. E. Florin and S. Baillet, “The brain's resting-state activity is shaped by synchronized cross-frequency coupling of neural oscillations,” NeuroImage, vol. 111, pp. 26–35, 2015. View at: Publisher Site | Google Scholar
  36. S. Dürschmid, T. Zaehle, K. Kopitzki et al., “Phase-amplitude cross-frequency coupling in the human nucleus accumbens tracks action monitoring during cognitive control,” Frontiers in Human Neuroscience, vol. 7, article 635, 2013. View at: Publisher Site | Google Scholar
  37. P. Lakatos, A. S. Shah, K. H. Knuth, I. Ulbert, G. Karmos, and C. E. Schroeder, “An oscillatory hierarchy controlling neuronal excitability and stimulus processing in the auditory cortex,” Journal of Neurophysiology, vol. 94, no. 3, pp. 1904–1911, 2005. View at: Publisher Site | Google Scholar
  38. J. Lee and J. Jeong, “Correlation of risk-taking propensity with cross-frequency phase-amplitude coupling in the resting EEG,” Clinical Neurophysiology, vol. 124, no. 11, pp. 2172–2180, 2013. View at: Publisher Site | Google Scholar
  39. R. J. McGinn and T. A. Valiante, “Phase-amplitude coupling and interlaminar synchrony are correlated in human neocortex,” Journal of Neuroscience, vol. 34, no. 48, pp. 15923–15930, 2014. View at: Publisher Site | Google Scholar
  40. T. Demiralp, Z. Bayraktaroglu, D. Lenz et al., “Gamma amplitudes are coupled to theta phase in human EEG during visual perception,” International Journal of Psychophysiology, vol. 64, no. 1, pp. 24–30, 2007. View at: Publisher Site | Google Scholar
  41. S. M. Szczepanski, N. E. Crone, R. A. Kuperman et al., “Dynamic changes in phase-amplitude coupling facilitate spatial attention control in fronto-parietal cortex,” PLoS Biology, vol. 12, no. 8, Article ID e1001936, 2014. View at: Publisher Site | Google Scholar
  42. M. van Wingerden, R. van der Meij, T. Kalenscher, E. Maris, and C. M. A. Pennartz, “Phase-amplitude coupling in rat orbitofrontal cortex discriminates between correct and incorrect decisions during associative learning,” The Journal of Neuroscience, vol. 34, no. 2, pp. 493–505, 2014. View at: Publisher Site | Google Scholar
  43. J. Wang, D. Li, X. Li et al., “Phase-amplitude coupling between theta and gamma oscillations during nociception in rat electroencephalography,” Neuroscience Letters, vol. 499, no. 2, pp. 84–87, 2011. View at: Publisher Site | Google Scholar
  44. C. Nakatani, A. Raffone, and C. van Leeuwen, “Efficiency of conscious access improves with coupling of slow and fast neural oscillations,” Journal of Cognitive Neuroscience, vol. 26, no. 5, pp. 1168–1179, 2014. View at: Publisher Site | Google Scholar
  45. R. C. Sotero, A. Bortel, S. Naaman et al., “Laminar distribution of cross-frequency coupling during spontaneous activity in rat area S1Fl,” in Proceedings of the 43rd Annual Meeting of the Society-for-Neuroscience, San Diego, Calif, USA, 2013. View at: Google Scholar
  46. M. X. Cohen, C. E. Elger, and J. Fell, “Oscillatory activity and phase-amplitude coupling in the human medial frontal cortex during decision making,” Journal of Cognitive Neuroscience, vol. 21, no. 2, pp. 390–402, 2009. View at: Publisher Site | Google Scholar
  47. J. Ito, P. Maldonado, and S. Grün, “Cross-frequency interaction of the eye-movement related LFP signals in V1 of freely viewing monkeys,” Frontiers in Systems Neuroscience, vol. 7, article 1, 2013. View at: Google Scholar
  48. E. Spaak, M. Bonnefond, A. Maier, D. A. Leopold, and O. Jensen, “Layer-specific entrainment of gamma-band neural activity by the alpha rhythm in monkey visual cortex,” Current Biology, vol. 22, no. 24, pp. 2313–2318, 2012. View at: Publisher Site | Google Scholar
  49. T. Yanagisawa, O. Yamashita, M. Hirata et al., “Regulation of motor representation by phase-amplitude coupling in the sensorimotor cortex,” The Journal of Neuroscience, vol. 32, no. 44, pp. 15467–15475, 2012. View at: Publisher Site | Google Scholar
  50. C. Magri, A. Mazzoni, N. K. Logothetis, and S. Panzeri, “Optimal band separation of extracellular field potentials,” Journal of Neuroscience Methods, vol. 210, no. 1, pp. 66–78, 2012. View at: Publisher Site | Google Scholar
  51. W. D. Penny, E. Duzel, K. J. Miller, and J. G. Ojemann, “Testing for nested oscillation,” Journal of Neuroscience Methods, vol. 174, no. 1, pp. 50–61, 2008. View at: Publisher Site | Google Scholar
  52. A. B. L. Tort, R. Komorowski, H. Eichenbaum, and N. Kopell, “Measuring phase-amplitude coupling between neuronal oscillations of different frequencies,” Journal of Neurophysiology, vol. 104, no. 2, pp. 1195–1210, 2010. View at: Publisher Site | Google Scholar
  53. B. Voytek, M. D'Esposito, N. Crone, and R. T. Knight, “A method for event-related phase/amplitude coupling,” NeuroImage, vol. 64, no. 1, pp. 416–424, 2013. View at: Publisher Site | Google Scholar
  54. D. Dvorak and A. A. Fenton, “Toward a proper estimation of phase-amplitude coupling in neural oscillations,” Journal of Neuroscience Methods, vol. 225, pp. 42–56, 2014. View at: Publisher Site | Google Scholar
  55. B. Pittman-Polletta, W.-H. Hsieh, S. Kaur, M.-T. Lo, and K. Hu, “Detecting phase-amplitude coupling with high frequency resolution using adaptive decompositions,” Journal of Neuroscience Methods, vol. 226, pp. 15–32, 2014. View at: Publisher Site | Google Scholar
  56. T. E. Özkurt, “Statistically reliable and fast direct estimation of phase-amplitude cross-frequency coupling,” IEEE Transactions on Biomedical Engineering, vol. 59, no. 7, pp. 1943–1950, 2012. View at: Publisher Site | Google Scholar
  57. M. A. Kramer and U. T. Eden, “Assessment of cross-frequency coupling with confidence using generalized linear models,” Journal of Neuroscience Methods, vol. 220, no. 1, pp. 64–74, 2013. View at: Publisher Site | Google Scholar
  58. J. Aru, J. Aru, V. Priesemann et al., “Untangling cross-frequency coupling in neuroscience,” Current Opinion in Neurobiology C, vol. 31, pp. 51–61, 2014. View at: Google Scholar
  59. R. D. Traub, R. Duncan, A. J. C. Russell et al., “Spatiotemporal patterns of electrocorticographic very fast oscillations (> 80 Hz) consistent with a network model based on electrical coupling between principal neurons,” Epilepsia, vol. 51, no. 8, pp. 1587–1597, 2010. View at: Publisher Site | Google Scholar
  60. J. Rinzel, D. Terman, X.-J. Wang, and B. Ermentrout, “Propagating activity patterns in large-scale inhibitory neuronal networks,” Science, vol. 279, no. 5355, pp. 1351–1355, 1998. View at: Publisher Site | Google Scholar
  61. H. R. Wilson and J. D. Cowan, “Excitatory and inhibitory interactions in localized populations of model neurons,” Biophysical Journal, vol. 12, no. 1, pp. 1–24, 1972. View at: Publisher Site | Google Scholar
  62. B. H. Jansen and V. G. Rit, “Electroencephalogram and visual-evoked potential generation in a mathematical-model of coupled cortical columns,” Biological Cybernetics, vol. 73, no. 4, pp. 357–366, 1995. View at: Publisher Site | Google Scholar
  63. H. Hasegawa, “Dynamical mean-field theory of spiking neuron ensembles: response to a single spike with independent noises,” Physical Review E, vol. 67, no. 4, Article ID 041903, 2003. View at: Publisher Site | Google Scholar
  64. H. Hasegawa, “Dynamical mean-field theory of noisy spiking neuron ensembles: application to the Hodgkin-Huxley model,” Physical Review E, vol. 68, no. 4, Article ID 041909, 2003. View at: Google Scholar
  65. R. Rodriguez and H. C. Tuckwell, “Noisy spiking neurons and networks: useful approximations for firing probabilities and global behavior,” BioSystems, vol. 48, no. 1–3, pp. 187–194, 1998. View at: Publisher Site | Google Scholar
  66. O. David and K. J. Friston, “A neural mass model for MEG/EEG: coupling and neuronal dynamics,” NeuroImage, vol. 20, no. 3, pp. 1743–1755, 2003. View at: Publisher Site | Google Scholar
  67. R. C. Sotero, N. J. Trujillo-Barreto, Y. Iturria-Medina, F. Carbonell, and J. C. Jimenez, “Realistically coupled neural mass models can generate EEG rhythms,” Neural Computation, vol. 19, no. 2, pp. 478–512, 2007. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  68. M. Zavaglia, L. Astolfi, F. Babiloni, and M. Ursino, “A neural mass model for the simulation of cortical activity estimated from high resolution EEG during cognitive or motor tasks,” Journal of Neuroscience Methods, vol. 157, no. 2, pp. 317–329, 2006. View at: Publisher Site | Google Scholar
  69. P. A. Valdes-Sosa, J. M. Sanchez-Bornot, R. C. Sotero et al., “Model driven EEG/fMRI fusion of brain oscillations,” Human Brain Mapping, vol. 30, no. 9, pp. 2701–2721, 2009. View at: Publisher Site | Google Scholar
  70. F. Cona, M. Zavaglia, and M. Ursino, “Binding and segmentation via a neural mass model trained with hebbian and anti-Hebbian mechanisms,” International Journal of Neural Systems, vol. 22, no. 2, Article ID 1250003, 2012. View at: Publisher Site | Google Scholar
  71. R. C. Sotero and A. Shmuel, “Energy-based stochastic control of neural mass models suggests time-varying effective connectivity in the resting state,” Journal of Computational Neuroscience, vol. 32, no. 3, pp. 563–576, 2012. View at: Publisher Site | Google Scholar
  72. O. David, S. J. Kiebel, L. M. Harrison, J. Mattout, J. M. Kilner, and K. J. Friston, “Dynamic causal modeling of evoked responses in EEG and MEG,” NeuroImage, vol. 30, no. 4, pp. 1255–1272, 2006. View at: Publisher Site | Google Scholar
  73. R. C. Sotero and N. J. Trujillo-Barreto, “Biophysical model for integrating neuronal activity, EEG, fMRI and metabolism,” NeuroImage, vol. 39, no. 1, pp. 290–309, 2008. View at: Publisher Site | Google Scholar
  74. R. C. Sotero, A. Bortel, R. Martínez-Cancino et al., “Anatomically-constrained effective connectivity among layers in a cortical column modeled and estimated from local field potentials,” Journal of Integrative Neuroscience, vol. 9, no. 4, pp. 355–379, 2010. View at: Publisher Site | Google Scholar
  75. F. Wendling, J. J. Bellanger, F. Bartolomei, and P. Chauvel, “Relevance of nonlinear lumped-parameter models in the analysis of depth-EEG epileptic signals,” Biological Cybernetics, vol. 83, no. 4, pp. 367–378, 2000. View at: Publisher Site | Google Scholar
  76. F. Wendling and F. Bartolomei, “Modeling EEG signals and interpreting measures of relationship during temporal-lobe seizures: an approach to the study of epileptogenic networks,” Epileptic Disorders, vol. 3, no. 1, pp. SI67–SI78, 2001. View at: Google Scholar
  77. B. S. Bhattacharya, D. Coyle, and L. P. Maguire, “A thalamo-cortico-thalamic neural mass model to study alpha rhythms in Alzheimer's disease,” Neural Networks, vol. 24, no. 6, pp. 631–645, 2011. View at: Publisher Site | Google Scholar
  78. X.-J. Wang and G. Buzsáki, “Gamma oscillation by synaptic inhibition in a hippocampal interneuronal network model,” Journal of Neuroscience, vol. 16, no. 20, pp. 6402–6413, 1996. View at: Google Scholar
  79. R. C. Sotero, “Generation of phase-amplitude coupling of neurophysiological signals in a neural mass model of a cortical column,” BioRxiv, 2015. View at: Publisher Site | Google Scholar
  80. S. A. Neymotin, K. M. Jacobs, A. A. Fenton, and W. W. Lytton, “Synaptic information transfer in computer models of neocortical columns,” Journal of Computational Neuroscience, vol. 30, no. 1, pp. 69–84, 2011. View at: Publisher Site | Google Scholar
  81. L. Tattini, S. Olmi, and A. Torcini, “Coherent periodic activity in excitatory Erdös-Renyi neural networks: the role of network connectivity,” Chaos, vol. 22, no. 2, Article ID 023133, 2012. View at: Publisher Site | Google Scholar
  82. X. S. Liang, “Unraveling the cause-effect relation between time series,” Physical Review E, vol. 90, no. 5, Article ID 052150, 2014. View at: Publisher Site | Google Scholar

Copyright © 2015 Roberto C. Sotero. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1484 Views | 713 Downloads | 8 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19.