Topology, Cross-Frequency, and Same-Frequency Band Interactions Shape the Generation of Phase-Amplitude Coupling in a Neural Mass Model of a Cortical Column

Abstract

Phase-amplitude coupling (PAC), a type of cross-frequency coupling (CFC) where the phase of a low-frequency rhythm modulates the amplitude of a higher frequency, is becoming an important indicator of information transmission in the brain. However, the neurobiological mechanisms underlying its generation remain undetermined. A realistic, yet tractable computational model of the phenomenon is thus needed. Here we analyze a neural mass model of a cortical column, comprising fourteen neuronal populations distributed across four layers (L2/3, L4, L5 and L6). A control analysis showed that the conditional transfer entropy (cTE) measure is able to correctly estimate the flow of information between neuronal populations. Then, we computed cTE from the phases to the amplitudes of the oscillations generated in the cortical column. This approach provides information regarding directionality by distinguishing PAC from APC (amplitude-phase coupling), i.e. the information transfer from amplitudes to phases, and can be used to estimate other types of CFC such as amplitude-amplitude coupling (AAC) and phase-phase coupling (PPC). While experiments often only focus on one or two PAC combinations (e.g., theta-gamma or alpha-gamma), we found that a cortical column can simultaneously generate almost all possible PAC combinations, depending on connectivity parameters, time constants, and external inputs. PAC interactions with and without an anatomical equivalent (direct and indirect interactions, respectively) were analyzed. We found that the strength of PAC between two populations was strongly correlated with the strength of the effective connections between the populations and, on average, did not depend on whether the PAC connection was direct or indirect. When considering a cortical column circuit as a complex network, we found that neuronal populations making indirect PAC connections had, on average, higher local clustering coefficient, efficiency, and betweenness centrality than populations making direct connections and populations not involved in PAC connections. This suggests that their interactions were more effective when transmitting information. Since approximately 60% of the obtained interactions represented indirect connections, our results highlight the importance of the topology of cortical circuits for the generation of the PAC phenomenon. Finally, our results demonstrated that indirect PAC interactions can be explained by a cascade of direct CFC and same-frequency band interactions, suggesting that PAC analysis of experimental data should be accompanied by the estimation of other types of frequency interactions for an integrative understanding of the phenomenon.

How to Build a Functional Connectomic Biomarker for Mild Cognitive Impairment From Source Reconstructed MEG Resting-State Activity: The Combination of ROI Representation and Connectivity Estimator Matters

“…A recent study concluded that the network topology, the CFC and the intra-frequency interactions shaped the PAC generation in a cortical column using a novel neural mass model (Sotero, 2016 ). Here, in order to reduce the computational time needed to run the pipeline from the reviewers in order to evaluate the whole analysis, we did not run surrogate analysis.…”

Design of optimal nonlinear network controllers for Alzheimer's disease

“…The inclusion of nonlinearities in our model causes several control tasks to fail for each subject, a fact that, to the best of our knowledge, has not been reported for linear brain dynamics. However, we do not expect that inputs to every neuronal conglomerate in real stimulation experiments are able to steer the (AD) brain to the desired state, given its complexity and nonlinear character [ 18 – 20 ]. The order in which areas were ranked according to the energy used for controlling the network, changed with the strength of the nonlinearity, γ .…”

“…Some studies have focused on identifying the most suitable sites for network controllability from a structural viewpoint only [ 17 ] while simplifying the dynamical interactions occurring on top of the connectivity scaffold. Other studies [ 9 – 12 ] used linear dynamics to model neural processes, which are known to be intrinsically nonlinear [ 13 , 18 – 20 ]. Hence, their predictions on neural network control should be taken with caution.…”

Detecting brain network communities: considering the role of information flow and its different temporal scales

“…Components are different as to conveyance of information (Sotero, 2016). For notation purposes, IMF1 denotes the fastest mode (highest frequency).…”