1
|
Assadzadeh S, Annen J, Sanz L, Barra A, Bonin E, Thibaut A, Boly M, Laureys S, Gosseries O, Robinson PA. Method for quantifying arousal and consciousness in healthy states and severe brain injury via EEG-based measures of corticothalamic physiology. J Neurosci Methods 2023; 398:109958. [PMID: 37661056 DOI: 10.1016/j.jneumeth.2023.109958] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2023] [Revised: 08/09/2023] [Accepted: 08/27/2023] [Indexed: 09/05/2023]
Abstract
BACKGROUND Characterization of normal arousal states has been achieved by fitting predictions of corticothalamic neural field theory (NFT) to electroencephalographic (EEG) spectra to yield relevant physiological parameters. NEW METHOD A prior fitting method is extended to distinguish conscious and unconscious states in healthy and brain injured subjects by identifying additional parameters and clusters in parameter space. RESULTS Fits of NFT predictions to EEG spectra are used to estimate neurophysiological parameters in healthy and brain injured subjects. Spectra are used from healthy subjects in wake and sleep and from patients with unresponsive wakefulness syndrome, in a minimally conscious state (MCS), and emerged from MCS. Subjects cluster into three groups in parameter space: conscious healthy (wake and REM), sleep, and brain injured. These are distinguished by the difference X-Y between corticocortical (X) and corticothalamic (Y) feedbacks, and by mean neural response rates α and β to incoming spikes. X-Y tracks consciousness in healthy individuals, with smaller values in wake/REM than sleep, but cannot distinguish between brain injuries. Parameters α and β differentiate deep sleep from wake/REM and brain injury. COMPARISON WITH EXISTING METHODS Other methods typically rely on laborious clinical assessment, manual EEG scoring, or evaluation of measures like Φ from integrated information theory, for which no efficient method exists. In contrast, the present method can be automated on a personal computer. CONCLUSION The method provides a means to quantify consciousness and arousal in healthy and brain injured subjects, but does not distinguish subtypes of brain injury.
Collapse
Affiliation(s)
- S Assadzadeh
- School of Physics, The University of Sydney, NSW 2006, Australia; Center for Integrative Brain Function, The University of Sydney, NSW 2006, Australia
| | - J Annen
- Coma Science Group, GIGA-Consciousness, University of Liège, Belgium; Centre du Cerveau, University Hospital of Liège, Belgium
| | - L Sanz
- Coma Science Group, GIGA-Consciousness, University of Liège, Belgium; Centre du Cerveau, University Hospital of Liège, Belgium
| | - A Barra
- Coma Science Group, GIGA-Consciousness, University of Liège, Belgium; Centre du Cerveau, University Hospital of Liège, Belgium
| | - E Bonin
- Coma Science Group, GIGA-Consciousness, University of Liège, Belgium; Centre du Cerveau, University Hospital of Liège, Belgium
| | - A Thibaut
- Coma Science Group, GIGA-Consciousness, University of Liège, Belgium; Centre du Cerveau, University Hospital of Liège, Belgium
| | - M Boly
- Department of Psychiatry, University of Wisconsin-Madison, Madison, WI, USA; Department of Neurology, School of Medicine and Public Health, University of Wisconsin, Madison, WI, USA
| | - S Laureys
- Coma Science Group, GIGA-Consciousness, University of Liège, Belgium; Centre du Cerveau, University Hospital of Liège, Belgium; Joint International Research Unit on Consciousness, CERVO Brain Research Centre, U Laval, Canada; International Consciousness Science Institute, Hangzhou Normal University, Hangzhou, China
| | - O Gosseries
- Coma Science Group, GIGA-Consciousness, University of Liège, Belgium; Centre du Cerveau, University Hospital of Liège, Belgium
| | - P A Robinson
- School of Physics, The University of Sydney, NSW 2006, Australia; Center for Integrative Brain Function, The University of Sydney, NSW 2006, Australia.
| |
Collapse
|
2
|
Liu X, Robinson PA. Mutual consistency of multiple visual feature maps constrains combined map topology. Phys Rev E 2023; 107:064401. [PMID: 37464602 DOI: 10.1103/physreve.107.064401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2021] [Accepted: 05/09/2023] [Indexed: 07/20/2023]
Abstract
The primary visual cortex (V1) is the first cortical area that processes visual information relayed from the thalamus. The topologies permitted in joint ocular dominance (OD), orientation preference (OP), and direction preference (DP) maps in V1 are considered, with the aim of finding a maximally symmetric periodic case that can serve as the basis for perturbations toward natural realizations. It is shown that mutual consistency of the maps selects just two possible such lattice structures, and that one of these is much closer to experiment than the other. This comprises a hexagonal lattice of alternating positive and negative OP singularities, with each unit cell or hypercolumn containing four such singularities, each of which radiates three DP discontinuities that follow OP contours and end at OP singularities of opposite sign. Other DP discontinuities emanate at 90 degrees to the midpoints of the ones that link OP singularities, and cross OP contours perpendicularly. These features explain experimentally observed relationships between DP discontinuities and OP contours, including sudden approximately 90-degree changes of direction in the former.
Collapse
Affiliation(s)
- X Liu
- School of Physics, The University of Sydney, NSW 2006, Australia and Center for Integrative Brain Function, The University of Sydney, NSW 2006, Australia
| | - P A Robinson
- School of Physics, The University of Sydney, NSW 2006, Australia and Center for Integrative Brain Function, The University of Sydney, NSW 2006, Australia
| |
Collapse
|
3
|
Robinson PA, McInnes A, Sarkar S. Spatiotemporal evolution of urban populations and housing: A dynamic utility-driven market-mediated model. PLoS One 2023; 18:e0282583. [PMID: 37027371 PMCID: PMC10081807 DOI: 10.1371/journal.pone.0282583] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2022] [Accepted: 01/29/2023] [Indexed: 04/08/2023] Open
Abstract
A model of the spatiotemporal evolution of urban areas is developed that simultaneously includes the effects on household utility of geography, population density, income distribution, and household preference for characteristics of dwellings and neighbors. The result is a utility function whose structure is similar to that of the energy of interacting spin systems in external fields. Spatiotemporal housing market evolution then results via transactions driven by increases in utility and changes in numbers of households and dwellings. It is shown that the model successfully predicts formation of monocentric and polycentric urban areas, stratification by wealth, segregation due to preferences for housing or neighbors, and the balance of supply and demand. These results go well beyond those of prior models that each dealt with subsets of these phenomena, and do so within a single, unified framework. Potential generalizations are discussed and further applications are suggested.
Collapse
Affiliation(s)
- P A Robinson
- School of Physics, University of Sydney, Sydney, New South Wales, Australia
| | - A McInnes
- School of Physics, University of Sydney, Sydney, New South Wales, Australia
| | - Somwrita Sarkar
- School of Architecture, Design, and Planning, The University of Sydney, Sydney, New South Wales, Australia
| |
Collapse
|
4
|
Robinson PA. Ten rules for effective modeling. Neuroimage 2022; 263:119622. [PMID: 36096279 DOI: 10.1016/j.neuroimage.2022.119622] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2022] [Revised: 09/02/2022] [Accepted: 09/08/2022] [Indexed: 10/31/2022] Open
Abstract
Modeling a natural system such as the brain aims to deepen understanding and to help to explain and link multiple phenomena into a coherent picture. In any specific case, this requires a clear view of the aims of each modeling project, followed by coordinated selection of the model's style and components; theoretical, numerical, and statistical analysis methods; distillation and presentation of results; and resulting well supported conclusions. The ten rules presented here apply to modeling of the brain and other systems and are designed to assist in carrying out integrated modeling with valid and well-supported outcomes that effectively achieve the modeling aims; referees can also use them when assessing the validity of modeling in submitted manuscripts.
Collapse
Affiliation(s)
- P A Robinson
- School of Physics, The University of Sydney, Sydney, NSW 2006, Australia.
| |
Collapse
|
5
|
Aghili Yajadda MM, Robinson PA, Henderson JA. Generalized neural field theory of cortical plasticity illustrated by an application to the linear phase of ocular dominance column formation in primary visual cortex. Biol Cybern 2022; 116:33-52. [PMID: 34773503 DOI: 10.1007/s00422-021-00901-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/09/2021] [Accepted: 09/30/2021] [Indexed: 06/13/2023]
Abstract
Physiologically based neural field theory (NFT) is extended to encompass cortical plasticity dynamics. An illustrative application is provided which treats the evolution of the connectivity of left- and right-eye visual stimuli to neuronal populations in the primary visual cortex (V1), and the initial, linear phase of formation of approximately one-dimensional (1D) ocular dominance columns (ODCs) that sets their transverse spatial scale. This links V1 activity, structure, and physiology within a single theory that already accounts for a range of other brain activity and connectivity phenomena, thereby enabling ODC formation and many other phenomena to be interrelated and cortical parameters to be constrained across multiple domains. The results accord with experimental ODC widths for realistic cortical parameters and are based directly on a unified description of the neuronal populations involved, their connection strengths, and the neuronal activity they support. Other key results include simple analytic approximations for ODC widths and the parameters of maximum growth rate, constraints on cortical excitatory and inhibitory gains, elucidation of the roles of specific poles of the V1 response function, and the fact that ODCs are not formed when input stimuli are fully correlated between eyes. This work provides a basis for further generalization of NFT to model other plasticity phenomena, thereby linking them to the range multiscale phenomena accounted for by NFT.
Collapse
Affiliation(s)
- M M Aghili Yajadda
- School of Physics, University of Sydney, Sydney, NSW, 2006, Australia
- Center for Integrative Brain Function, University of Sydney, Sydney, NSW, 2006, Australia
| | - P A Robinson
- School of Physics, University of Sydney, Sydney, NSW, 2006, Australia
- Center for Integrative Brain Function, University of Sydney, Sydney, NSW, 2006, Australia
| | - J A Henderson
- School of Physics, University of Sydney, Sydney, NSW, 2006, Australia.
- Center for Integrative Brain Function, University of Sydney, Sydney, NSW, 2006, Australia.
| |
Collapse
|
6
|
Robinson PA. Integrals and series related to propagators of neural and haemodynamic waves. R Soc Open Sci 2021; 8:211562. [PMID: 34966557 PMCID: PMC8633804 DOI: 10.1098/rsos.211562] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/09/2021] [Accepted: 11/08/2021] [Indexed: 06/14/2023]
Abstract
The propagator, or Green function, of a class of neural activity fields and of haemodynamic waves is evaluated exactly. The results enable a number of related integrals to be evaluated, along with series expansions of key results in terms of Bessel functions of the second kind. Connections to other related equations are also noted.
Collapse
Affiliation(s)
- P. A. Robinson
- School of Physics, University of Sydney, Sydney, NSW 2006, Australia
- Center for Integrative Brain Function, University of Sydney, Sydney, NSW 2006, Australia
| |
Collapse
|
7
|
Abstract
The problem of finding a compact natural representation of brain dynamics and connectivity is addressed using an expansion in terms of physical spatial eigenmodes and their frequency resonances. It is demonstrated that this discrete expansion via the system transfer function enables linear and nonlinear dynamics to be analyzed in compact form in terms of natural dynamic "atoms," each of which is a frequency resonance of an eigenmode. Because these modal resonances are determined by the system dynamics, not the investigator, they are privileged over widely used phenomenological patterns, and obviate the need for artificial discretizations and thresholding in coordinate space. It is shown that modal resonances participate as nodes of a discrete spectral network, are noninteracting in the linear regime, but are linked nonlinearly by wave-wave coalescence and decay processes. The modal resonance formulation is shown to be capable of speeding numerical calculations of strongly nonlinear interactions. Recent work in brain dynamics, especially based on neural field theory (NFT) approaches, allows eigenmodes and their resonances to be estimated from data without assuming a specific brain model. This means that dynamic equations can be inferred using system identification methods from control theory, rather than being assumed, and resonances can be interpreted as control-systems data filters. The results link brain activity and connectivity with control-systems functions such as prediction and attention via gain control and can also be linked to specific NFT predictions if desired, thereby providing a convenient bridge between physiologically based theories and experiment. Amplitudes of modes and resonances can also be tracked to provide a more direct and temporally localized representation of the dynamics than correlations and covariances, which are widely used in the field. By synthesizing many different lines of research, this work provides a way to link quantitative electrophysiological and imaging measurements, connectivity, brain dynamics, and function. This underlines the need to move between coordinate and spectral representations as required. Moreover, standard theoretical-physics approaches and mathematical methods can be used in place of ad hoc statistical measures such as those based on graph theory of artificially discretized and decimated networks, which are highly prone to selection effects and artifacts.
Collapse
Affiliation(s)
- P A Robinson
- School of Physics, University of Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| |
Collapse
|
8
|
Mukta KN, Robinson PA, Pagès JC, Gabay NC, Gao X. Evoked response activity eigenmode analysis in a convoluted cortex via neural field theory. Phys Rev E 2021; 102:062303. [PMID: 33466049 DOI: 10.1103/physreve.102.062303] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2019] [Accepted: 07/15/2020] [Indexed: 11/07/2022]
Abstract
Neural field theory of the corticothalamic system is used to explore evoked response potentials (ERPs) caused by spatially localized impulse stimuli on the convoluted cortex and on a spherical cortex. Eigenfunctions are calculated analytically on the spherical cortex and numerically on the convoluted cortex via eigenfunction expansions. Eigenmodes on a convoluted cortex are similar to those of the spherical cortex, and a few such modes are found to be sufficient to reproduce the main ERP features. It is found that the ERP peak is stronger in spherical cortex than convoluted cortex, but in both cases the peak decreases monotonically with increasing distance from the stimulus point. In the convoluted case, cortical folding causes ERPs to differ between locations at the same distance from the stimulus point and spherical symmetries are only approximately preserved.
Collapse
Affiliation(s)
- K N Mukta
- School of Physics, University of Sydney, New South Wales 2006, Australia.,Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - P A Robinson
- School of Physics, University of Sydney, New South Wales 2006, Australia.,Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - J C Pagès
- School of Physics, University of Sydney, New South Wales 2006, Australia.,Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia.,School of Physics, University of Zurich, Zürich, Canton of Zürich, Switzerland
| | - N C Gabay
- School of Physics, University of Sydney, New South Wales 2006, Australia.,Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - Xiao Gao
- School of Physics, University of Sydney, New South Wales 2006, Australia.,Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| |
Collapse
|
9
|
Robinson PA. Neural field theory of neural avalanche exponents. Biol Cybern 2021; 115:237-243. [PMID: 33939016 DOI: 10.1007/s00422-021-00875-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/18/2020] [Accepted: 04/10/2021] [Indexed: 06/12/2023]
Abstract
The power-law exponents of observed size and lifetime distributions of near-critical neural avalanches are calculated from neural field theory using diagrammatic methods. This brings neural avalanches within the ambit of neural field theory, which has also previously explained near-critical 1/f spectra and many other observed features of neural activity. This strengthens the case for near-criticality of the brain and opens the way for these other phenomena to be interrelated with avalanches and their dynamics.
Collapse
Affiliation(s)
- P A Robinson
- School of Physics, The University of Sydney, Sydney, New South Wales, 2006, Australia.
- Center of Excellence for Integrative Brain Function, The University of Sydney, Sydney, New South Wales, 2006, Australia.
| |
Collapse
|
10
|
Robinson PA, Gao X, Han Y. Relationships between lognormal distributions of neural properties, activity, criticality, and connectivity. Biol Cybern 2021; 115:121-130. [PMID: 33825983 DOI: 10.1007/s00422-021-00871-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/27/2020] [Accepted: 03/15/2021] [Indexed: 06/12/2023]
Abstract
Relationships between convergence of inputs onto neurons, divergence of outputs from them, synaptic strengths, nonlinear firing response properties, and randomness of axonal ranges are systematically explored by interrelating means and variances of synaptic strengths, firing rates, and soma voltages. When self-consistency is imposed, it is found that broad distributions of synaptic strength are a necessary concomitant of the known massive convergence of inputs to individual neurons, and observed widths of lognormal distributions of synaptic strength and firing rate are explained provided the brain is in a near-critical state, consistent with independent observations. The strongest individual synapses are shown to have an effect on soma voltage comparable to the effect of all others combined, which supports suggestions that they may have a key role in neural communication. Remarkably, inclusion of moderate randomness in characteristic axonal ranges is shown to account for the observed [Formula: see text]-fold variability in two-point connectivity at a given separation and [Formula: see text]-fold overall when the known mean exponential fall-off is included, consistent with observed near-lognormal distributions. Inferred axonal deviations from straight-line paths are also consistent with independent estimates.
Collapse
Affiliation(s)
- P A Robinson
- School of Physics, University of Sydney, New South Wales 2006, Sydney, Australia.
- Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Sydney, Australia.
| | - Xiao Gao
- School of Physics, University of Sydney, New South Wales 2006, Sydney, Australia
- Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Sydney, Australia
| | - Y Han
- School of Physics, University of Sydney, New South Wales 2006, Sydney, Australia
- Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Sydney, Australia
| |
Collapse
|
11
|
Abstract
Spectral analysis and neural field theory are used to investigate the role of local connections in brain connectivity matrices (CMs) that quantify connectivity between pairs of discretized brain regions. This work investigates how the common procedure of omitting such self-connections (i.e., the diagonal elements of CMs) in published studies of brain connectivity affects the properties of functional CMs (fCMs) and the mutually consistent effective CMs (eCMs) that correspond to them. It is shown that retention of self-connections in the fCM calculated from two-point activity covariances is essential for the fCM to be a true covariance matrix, to enable correct inference of the direct total eCMs from the fCM, and to ensure their compatibility with it; the deCM and teCM represent the strengths of direct connections and all connections between points, respectively. When self-connections are retained, inferred eCMs are found to have net inhibitory self-connections that represent the local inhibition needed to balance excitation via white matter fibers at longer ranges. This inference of spatially unresolved connectivity exemplifies the power of spectral connectivity methods, which also enable transformation of CMs to compact diagonal forms that allow accurate approximation of the fCM and total eCM in terms of just a few modes, rather than the full [Formula: see text] CM entries for connections between N brain regions. It is found that omission of fCM self-connections affects both local and long-range connections in eCMs, so they cannot be omitted even when studying the large-scale. Moreover, retention of local connections enables inference of subgrid short-range inhibitory connectivity. The results are verified and illustrated using the NKI-Rockland dataset from the University of Southern California Multimodal Connectivity Database. Deletion of self-connections is common in the field; this does not affect case-control studies but the present results imply that such fCMs must have self-connections restored before eCMs can be inferred from them.
Collapse
Affiliation(s)
- Xiao Gao
- Department of Biomedical Engineering, University of Melbourne, Parkville, VIC 3052 Australia
- School of Physics, The University of Sydney, Sydney, NSW 2006 Australia
- Center of Excellence for Integrative Brain Function, The University of Sydney, Sydney, NSW 2006 Australia
| | - P. A. Robinson
- School of Physics, The University of Sydney, Sydney, NSW 2006 Australia
- Center of Excellence for Integrative Brain Function, The University of Sydney, Sydney, NSW 2006 Australia
| |
Collapse
|
12
|
Abstract
It is shown that the statistical properties of connections between regions of the brain and their dependence on coarse-graining and thresholding in published data can be reproduced by a simple distance-based physical connectivity model. This allows studies with differing parcellation and thresholding to be interrelated objectively, and for the results of future studies on more finely grained or differently thresholded networks to be predicted. As examples of the implications, it is shown that the dependences of network measures on thresholding and parcellation imply that chosen brain regions can appear to form a small world network, even though the network at finer scales, or ultimately of individual neurons, may not be small world networks themselves.
Collapse
Affiliation(s)
- T. C. Lacy
- School of Physics, University of Sydney, Sydney, NSW, Australia
- Center for Integrative Brain Function, University of Sydney, Sydney, NSW, Australia
| | - P. A. Robinson
- School of Physics, University of Sydney, Sydney, NSW, Australia
- Center for Integrative Brain Function, University of Sydney, Sydney, NSW, Australia
- * E-mail:
| |
Collapse
|
13
|
Ferdousi M, Babaie-Janvier T, Robinson PA. Nonlinear wave-wave interactions in the brain. J Theor Biol 2020; 500:110308. [PMID: 32389568 DOI: 10.1016/j.jtbi.2020.110308] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2019] [Revised: 03/30/2020] [Accepted: 04/27/2020] [Indexed: 11/24/2022]
Abstract
Neural field theory of the corticothalamic system is used to analyze nonlinear wave-wave interactions in steady state visual evoked potential responses. The nonlinear power spectrum is analytically calculated by convolving the linear power spectrum with itself and other factors. Periodic sine and square wave stimuli are used to generate steady state visual evoked potential responses and to study stimulus-driven nonlinear corticothalamic dynamic interactions. Moreover, we use dual sine drives to analyze the driven dynamics. Numerical analysis shows that the nonlinear power spectrum embodies key nonlinear features, including harmonic and subharmonic generation, entrainment of the alpha rhythm to periodic stimuli at the drive frequency, sum and difference frequencies due to wave-wave coalescence and decay. Further, the scaling properties of the key phenomena observed in nonlinear interactions are studied, verifying some of the theoretical predictions for these being generated by three-wave processes.
Collapse
Affiliation(s)
- M Ferdousi
- School of Physics, University of Sydney, New South Wales 2006, Australia; Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia.
| | - T Babaie-Janvier
- School of Physics, University of Sydney, New South Wales 2006, Australia; Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - P A Robinson
- School of Physics, University of Sydney, New South Wales 2006, Australia; Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| |
Collapse
|
14
|
Abstract
This paper generalizes and extends previous work on using neural field theory to quantitatively analyze the two-dimensional (2D) spatiotemporal correlation properties of gamma-band (30-70 Hz) oscillations evoked by stimuli arriving at the primary visual cortex, and modulated by patchy connectivities that depend on orientation preference (OP). Correlation functions are derived analytically for general stimulus and measurement conditions. The theoretical results reproduce a range of published experimental results. These include (i) the existence of two-point oscillatory temporal cross correlations with zero time lag between neurons with similar OP; (ii) the influence of spatial separation of neurons on the strength of the correlations; and (iii) the effects of differing stimulus orientations. They go beyond prior work by incorporating experimentally observed patchy projection patterns to predict the 2D correlation structure including both OP and ocular dominance effects, thereby relaxing assumptions of translational invariance implicit in prior one-dimensional analysis.
Collapse
Affiliation(s)
- X Liu
- School of Physics, University of Sydney, Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, Sydney, New South Wales 2006, Australia
| | - P Sanz-Leon
- School of Physics, University of Sydney, Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, Sydney, New South Wales 2006, Australia
| | - P A Robinson
- School of Physics, University of Sydney, Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, Sydney, New South Wales 2006, Australia
| |
Collapse
|
15
|
Deeba F, Sanz-Leon P, Robinson PA. Effects of physiological parameter evolution on the dynamics of tonic-clonic seizures. PLoS One 2020; 15:e0230510. [PMID: 32240175 PMCID: PMC7117716 DOI: 10.1371/journal.pone.0230510] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2020] [Accepted: 03/03/2020] [Indexed: 02/02/2023] Open
Abstract
The temporal and spectral characteristics of tonic-clonic seizures are investigated using a neural field model of the corticothalamic system in the presence of a temporally varying connection strength between the cerebral cortex and thalamus. Increasing connection strength drives the system into ∼ 10 Hz seizure oscillations once a threshold is passed and a subcritical Hopf bifurcation occurs. In this study, the spectral and temporal characteristics of tonic-clonic seizures are explored as functions of the relevant properties of physiological connection strengths, such as maximum strength, time above threshold, and the ramp rate at which the strength increases or decreases. Analysis shows that the seizure onset time decreases with the maximum connection strength and time above threshold, but increases with the ramp rate. Seizure duration and offset time increase with maximum connection strength, time above threshold, and rate of change. Spectral analysis reveals that the power of nonlinear harmonics and the duration of the oscillations increase as the maximum connection strength and the time above threshold increase. A secondary limit cycle at ∼ 18 Hz, termed a saddle-cycle, is also seen during seizure onset and becomes more prominent and robust with increasing ramp rate. If the time above the threshold is too small, the system does not reach the 10 Hz limit cycle, and only exhibits 18 Hz saddle-cycle oscillations. It is also seen that the time to reach the saturated large amplitude limit-cycle seizure oscillation from both the instability threshold and from the end of the saddle-cycle oscillations is inversely proportional to the square root of the ramp rate.
Collapse
Affiliation(s)
- F. Deeba
- School of Physics, University of Sydney, Sydney, NSW, Australia
- Center for Integrative Brain Function, University of Sydney, Sydney, NSW, Australia
- * E-mail: ,
| | - P. Sanz-Leon
- School of Physics, University of Sydney, Sydney, NSW, Australia
- Center for Integrative Brain Function, University of Sydney, Sydney, NSW, Australia
| | - P. A. Robinson
- School of Physics, University of Sydney, Sydney, NSW, Australia
- Center for Integrative Brain Function, University of Sydney, Sydney, NSW, Australia
| |
Collapse
|
16
|
Abstract
Hemodynamic modeling is used to explore the origin, predict, and analyze the power spectrum of the resting-state blood-oxygen-level-dependent (BOLD) signal measured by functional magnetic resonance imaging (fMRI), which has been reported to have a power-law form, i.e., P(f)∝f^{-s}, where P(f) is the power, f is the frequency, and s>0 is the power-law exponent. However, current fMRI experimental paradigms have limited acquisition durations, affecting the spectral resolution of fMRI data at the low-frequency regime. Here, the claimed power-law spectrum is investigated by using a recent hemodynamic model to analytically derive the BOLD power spectrum, with parameters that are related to neurophysiology. The theoretical results show that, for all realistic parameter combinations, the BOLD power spectrum is flat at f≲0.01Hz, has a weak resonance originating from intrinsic oscillations of vasodilatory response, and becomes a power law for high frequencies, all of which is in agreement with an empirical data set that describes the spectrum of one subject and brain region. However, the results are contrary to studies reporting a pure power-law spectrum at f≲0.2Hz. The discrepancy is attributed largely to data averaging employed by current approaches that averages together important properties of the BOLD power spectrum, such as its resonance, that biases the spectrum to only show a power law. Data averaging also reduces the high-frequency power-law exponent relative to individual cases. Overall, this work demonstrates how the model can reproduce BOLD dynamics and further analyze its low-frequency behavior. Moreover, it also uses the model to explain the impact of procedures, such as data averaging, on the reported features of the BOLD power spectrum.
Collapse
Affiliation(s)
- J C Pang
- School of Physics, University of Sydney, Sydney, NSW 2006, Australia.,Center for Integrative Brain Function, University of Sydney, Sydney, NSW 2006, Australia.,QIMR Berghofer Medical Research Institute, Herston, QLD 4006, Australia
| | - P A Robinson
- School of Physics, University of Sydney, Sydney, NSW 2006, Australia.,Center for Integrative Brain Function, University of Sydney, Sydney, NSW 2006, Australia
| |
Collapse
|
17
|
Palczynski LJ, Bleach ECL, Brennan ML, Robinson PA. Giving calves 'the best start': Perceptions of colostrum management on dairy farms in England. Anim Welf 2020. [DOI: 10.7120/09627286.29.1.045] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Good colostrum management can confer protective immunity to newborn calves, making calves less susceptible to infectious disease, and fundamentally improving both their short- and long-term health, welfare and productivity. Industry recommendations commonly refer to 'The Three 'Q's'
of colostrum management: the need for calves to receive sufficient 'Quantity' of high 'Quality' colostrum 'Quickly' after birth; some also include 'sQueaky clean' and 'Quantification of passive transfer'. However, research to date suggests that the failure of passive transfer of colostral
antibodies is common on commercial dairy farms, contributing to sub-optimal calf health and mortality. This paper explores why this may be the case by investigating stakeholder perceptions of colostrum management and how these perceptions might affect the practice of ensuring adequate colostrum
administration to newborn calves. Calf rearing and youngstock management practices on English dairy farms were investigated using 40 in-depth semi-structured interviews: 26 with dairy farmers and 14 with advisors (including veterinarians, feed and pharmaceutical company representatives). Interviews
were audio-recorded, transcribed and thematically coded for analysis. 'The Three 'Q's' were found to act as useful reminders about the goals of colostrum management, and a case can be made for further publicising the inclusion of 'sQueaky clean' and 'Quantification of passive transfer' as
there remains a lack of focus on colostrum hygiene and measurement of successful antibody transfer. Knowledge of the 'Q's' did not guarantee implementation, and time and labour constraints alongside farmer misconceptions must be addressed when offering professional advice on improving calf
health. Further research to encourage on-farm collection and analysis of monitoring data including rates of passive transfer is particularly needed. Advisors must not overlook the importance of colostrum management when assessing farm practices and ensure that they promote evidence-based recommendations
if dairy calf morbidity and mortality is to be reduced.
Collapse
|
18
|
Abstract
The dynamics of interictal events between absence seizures and their relationship to seizures themselves are investigated by employing a neural field model of the corticothalamic system. Interictal events are modeled as being due to transient parameter excursions beyond the seizure threshold, in the present case by sufficiently temporally varying the connection strength between the cerebral cortex and the thalamus. Increasing connection strength drives the system into ∼3-Hz seizure oscillations via a supercritical Hopf bifurcation once the linear instability threshold is passed. Depending on the time course of the excursion above threshold, different interictal activity event dynamics are seen in the time series of corticothalamic fields. These resemble experimental interictal time series observed via electroencephalography. It is found that the morphology of these events depends on the magnitude and duration of the excursion above threshold. For a large-amplitude excursion of short duration, events resemble interictal spikes, where one large spike is seen, followed by small damped oscillations. For a short excursion with long duration, events like observed interictal periodic sharp waves are seen. When both amplitude and duration above threshold are large, seizure oscillations are seen. Using these outcomes, proximity to seizure can be estimated and tracked.
Collapse
Affiliation(s)
- F Deeba
- Department of Physics, Dhaka University of Engineering and Technology, Gazipur 1700, Bangladesh; School of Physics, University of Sydney, New South Wales 2006, Australia; and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - P Sanz-Leon
- School of Physics, University of Sydney, New South Wales 2006, Australia, and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - P A Robinson
- School of Physics, University of Sydney, New South Wales 2006, Australia, and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| |
Collapse
|
19
|
Yang DP, Robinson PA. Unified analysis of global and focal aspects of absence epilepsy via neural field theory of the corticothalamic system. Phys Rev E 2019; 100:032405. [PMID: 31639915 DOI: 10.1103/physreve.100.032405] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2019] [Indexed: 06/10/2023]
Abstract
Absence epilepsy is characterized by a sudden paroxysmal loss of consciousness accompanied by oscillatory activity propagating over many brain areas. Although primary generalized absence seizures are supported by the global corticothalamic system, converging experimental evidence supports a focal theory of absence epilepsy. Here a physiology-based corticothalamic model is investigated with spatial heterogeneity due to focal epilepsy to unify global and focal aspects of absence epilepsy. Numeric and analytic calculations are employed to investigate the emergent spatiotemporal dynamics as well as their underlying dynamical mechanisms. They can be categorized into three scenarios: suppressed epilepsy, focal seizures, or generalized seizures, as summarized from a phase diagram vs focal width and characteristic axon range. The corresponding temporal frequencies and spatial extents of cortical waves in generalized seizures and focal seizures agree well with experimental observations of global and focal aspects of absence epilepsy, respectively. The emergence of the spatiotemporal dynamics corresponding to focal seizures provides a biophysical explanation of the temporally higher frequency but spatially more localized cortical waves observed in genetic rat models that display characteristics of human absence epilepsy. Predictions are also presented for further experimental test.
Collapse
Affiliation(s)
- Dong-Ping Yang
- School of Physics, University of Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - P A Robinson
- School of Physics, University of Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| |
Collapse
|
20
|
Abstract
Effective connectivity embodied in transfer functions is derived from symmetric-network activity correlations under task-free conditions via a recent causal spectral factorization method. This generalizes previous covariance-based analyses to include frequency dependencies and time delays. Results are verified against analytic solutions of equations that have reproduced many aspects of experimental brain dynamics and against cases of more complex connectivity. Robustness to noise is also demonstrated.
Collapse
Affiliation(s)
- J N MacLaurin
- School of Physics, University of Sydney, New South Wales 2006, Australia, and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - P A Robinson
- School of Physics, University of Sydney, New South Wales 2006, Australia, and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| |
Collapse
|
21
|
Abstract
Evoked response potentials (ERPs) are calculated in spherical and planar geometries using neural field theory of the corticothalamic system. The ERP is modeled as an impulse response and the resulting modal effects of spherical corticothalamic dynamics are explored, showing that results for spherical and planar geometries converge in the limit of large brain size. Cortical modal effects can lead to a double-peak structure in the ERP time series. It is found that the main difference between infinite planar geometry and spherical geometry is that the ERP peak is sharper and stronger in the spherical geometry. It is also found that the magnitude of the response decreases with increasing spatial width of the stimulus at the cortex. The peak is slightly delayed at large angles from the stimulus point, corresponding to group velocities of 6-10 m s^{-1}. Strong modal effects are found in the spherical geometry, with the lowest few modes sufficing to describe the main features of ERPs, except very near to spatially narrow stimuli.
Collapse
Affiliation(s)
- K N Mukta
- School of Physics, University of Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - Xiao Gao
- School of Physics, University of Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - P A Robinson
- School of Physics, University of Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| |
Collapse
|
22
|
Robinson PA. Neural field theory of effects of brain modifications and lesions on functional connectivity: Acute effects, short-term homeostasis, and long-term plasticity. Phys Rev E 2019; 99:042407. [PMID: 31108595 DOI: 10.1103/physreve.99.042407] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2018] [Indexed: 11/07/2022]
Abstract
Neural field theory is used to predict the functional connectivity effects of lesions or other modifications to effective connectivity. Widespread initial changes are predicted after localized or diffuse changes to white or gray matter, consistent with observations, and enabling lesion severity indexes to be defined. It is shown how short-term homeostasis and longer-term plasticity can reduce perturbations while maintaining brain criticality under conditions where some connections remain fixed because of damage in the lesion core. The extent to which such effects can compensate for initial connectivity changes is then explored, showing that the strongest corrective changes are concentrated toward the edges of the perturbation if it is localized and its core is fixed. The results are applicable to inferring underlying connectivity changes and to interpreting and monitoring functional connectivity modifications after lesions, injury, surgery, drugs, or brain stimulation.
Collapse
Affiliation(s)
- P A Robinson
- School of Physics, University of Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| |
Collapse
|
23
|
Abstract
Brain connectivity and structure-function relationships are analyzed from a physical perspective in place of common graph-theoretic and statistical approaches that overwhelmingly ignore the brain's physical structure and geometry. Field theory is used to define connectivity tensors in terms of bare and dressed propagators, and discretized representations are implemented that respect the physical nature and dimensionality of the quantities involved, retain the correct continuum limit, and enable diagrammatic analysis. Eigenfunction analysis is used to simultaneously characterize and probe patterns of brain connectivity and activity, in place of statistical or phenomenological patterns. Physically based measures that characterize the connectivity are then developed in coordinate and spectral domains; some of which generalize or rectify graph-theoretic measures to implement correct dimensionality and continuum limits, and some replace graph-theoretic quantities. Traditional graph-based measures are shown to be highly prone to artifacts introduced by discretization and threshold, often because essential physical constraints have not been imposed, dimensionality has not been included, and/or distinctions between scalar, vector, and tensor quantities have not been considered. The results can replace them in ways that converge correctly and measure properties of brain structure, rather than of its discretization, and thus potentially enable physical interpretation of the many phenomenological results in the literature. Geometric effects are shown to dominate in determining many brain properties and care must be taken not to interpret geometric differences as differences in intrinsic neural connectivity. The results demonstrate the need to use systematic physical methods to analyze the brain and the potential of such methods to obtain new insights from data, make new predictions for experimental test, and go beyond phenomenological classification to dynamics and mechanisms.
Collapse
Affiliation(s)
- P A Robinson
- School of Physics, University of Sydney, Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, Sydney, New South Wales 2006, Australia
| |
Collapse
|
24
|
Assadzadeh S, Robinson PA. Necessity of the sleep-wake cycle for synaptic homeostasis: system-level analysis of plasticity in the corticothalamic system. R Soc Open Sci 2018; 5:171952. [PMID: 30473798 PMCID: PMC6227995 DOI: 10.1098/rsos.171952] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/03/2018] [Accepted: 09/13/2018] [Indexed: 06/09/2023]
Abstract
Neural field theory is used to study the system-level effects of plasticity in the corticothalamic system, where arousal states are represented parametrically by the connection strengths of the system, among other physiologically based parameters. It is found that the plasticity dynamics have no fixed points or closed cycles in the parameter space of the connection strengths, but parameter subregions exist where flows have opposite signs. Remarkably, these subregions coincide with previously identified regions that correspond to wake and slow-wave sleep, thus demonstrating state dependence of the sign of synaptic modification. We then show that a closed cycle in the parameter space is possible when the plasticity dynamics are driven by the ascending arousal system, which cycles the brain between sleep and wake to complete a closed loop that includes arcs through the opposite-flow subregions. Thus, it is concluded that both wake and sleep are necessary, and together are able to stabilize connection weights in the brain over the daily cycle, thereby providing quantitative realization of the synaptic homeostasis hypothesis.
Collapse
Affiliation(s)
- S. Assadzadeh
- School of Physics, The University of Sydney, New South Wales 2006, Australia
- Center for Integrative Brain Function, The University of Sydney, New South Wales 2006, Australia
| | - P. A. Robinson
- School of Physics, The University of Sydney, New South Wales 2006, Australia
- Center for Integrative Brain Function, The University of Sydney, New South Wales 2006, Australia
| |
Collapse
|
25
|
Abstract
Evoked response potentials (ERPs) and other transients are modeled as impulse responses using physiology-based neural field theory (NFT) of the corticothalamic system of neural activity in the human brain that incorporates synaptic and dendritic dynamics, firing response, axonal propagation, and corticocortical and corticothalamic pathways. The properties of model-predicted ERPs are explored throughout the stability zone of the corticothalamic system, and predicted time series and wavelet spectra are also analyzed. This provides a unified treatment of predicted ERPs for both normal and abnormal states within the brain's stability zone, including likely parameters to represent abnormal states of reduced arousal.
Collapse
Affiliation(s)
- M S Zobaer
- School of Physics, The University of Sydney, Sydney, NSW, 2006, Australia.
- Center for Integrative Brain Function, The University of Sydney, Sydney, NSW, 2006, Australia.
- Center for Research Excellence, Neurosleep, 431 Glebe Point Rd, Glebe, NSW, 2037, Australia.
- Department of Physics, Bangladesh University of Textiles, Dhaka, 1208, Bangladesh.
| | - P A Robinson
- School of Physics, The University of Sydney, Sydney, NSW, 2006, Australia
- Center for Integrative Brain Function, The University of Sydney, Sydney, NSW, 2006, Australia
- Center for Research Excellence, Neurosleep, 431 Glebe Point Rd, Glebe, NSW, 2037, Australia
| | - C C Kerr
- School of Physics, The University of Sydney, Sydney, NSW, 2006, Australia
- Center for Integrative Brain Function, The University of Sydney, Sydney, NSW, 2006, Australia
- Department of Physiology and Pharmacology, State University of New York Downstate Medical Center, 450 Clarkson Ave, Brooklyn, NY, USA
| |
Collapse
|
26
|
Roy N, Sanz-Leon P, Robinson PA. Spectrum of connectivity fluctuations including the effect of activity-dependent feedback. Phys Rev E 2018; 98:022319. [PMID: 30253627 DOI: 10.1103/physreve.98.022319] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2018] [Indexed: 11/07/2022]
Abstract
The spatiotemporal spectrum of feedback-driven fluctuations of brain connectivity is investigated using nonlinear neural field theory of the corticothalamic system. Weakly nonlinear dynamics of neural feedbacks are expanded in terms of first order perturbations of neural activity relative to a fixed point. Susceptibilities are used to quantify the change in connectivity per unit change in presynaptic or postsynaptic activity caused by nonlinear feedbacks such as facilitation, depression, sensitization, potentiation, and the effects of discrete eigenmode structure are included for a spherical brain geometry. Spectral signatures such as resonances are identified that allow the presence of particular presynaptic and postsynaptic feedback effects to be inferred. These include additional resonances at high frequencies and shifts of existing spectral peaks, mostly visible in the lowest spatial modes of the response.
Collapse
Affiliation(s)
- N Roy
- School of Physics, University of Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - P Sanz-Leon
- School of Physics, University of Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - P A Robinson
- School of Physics, University of Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| |
Collapse
|
27
|
Deeba F, Sanz-Leon P, Robinson PA. Dependence of absence seizure dynamics on physiological parameter evolution. J Theor Biol 2018; 454:11-21. [PMID: 29807025 DOI: 10.1016/j.jtbi.2018.05.029] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2018] [Revised: 05/22/2018] [Accepted: 05/24/2018] [Indexed: 12/30/2022]
Abstract
A neural field model of the corticothalamic system is applied to investigate the temporal and spectral characteristics of absence seizures in the presence of a temporally varying connection strength between the cerebral cortex and thalamus. Increasing connection strength drives the system into an absence seizure-like state once a threshold is passed and a supercritical Hopf bifurcation occurs. The dynamics and spectral characteristics of the resulting model seizures are explored as functions of maximum connection strength, time above threshold, and the rate at which the connection strength increases (ramp rate). Our results enable spectral and temporal characteristics of seizures to be related to changes in the underlying physiological evolution of connections via nonlinear dynamics and neural field theory. Spectral analysis reveals that the power of the harmonics and the duration of the oscillations increase as the maximum connection strength and the time above threshold increase. It is also found that the time to reach the stable limit-cycle seizure oscillation from the instability threshold decreases with the square root of the ramp rate.
Collapse
Affiliation(s)
- F Deeba
- School of Physics, University of Sydney, NSW 2006, Australia; Center for Integrative Brain Function, University of Sydney, NSW 2006, Australia.
| | - Paula Sanz-Leon
- School of Physics, University of Sydney, NSW 2006, Australia; Center for Integrative Brain Function, University of Sydney, NSW 2006, Australia
| | - P A Robinson
- School of Physics, University of Sydney, NSW 2006, Australia; Center for Integrative Brain Function, University of Sydney, NSW 2006, Australia
| |
Collapse
|
28
|
Robinson PA, Pagès JC, Gabay NC, Babaie T, Mukta KN. Neural field theory of perceptual echo and implications for estimating brain connectivity. Phys Rev E 2018; 97:042418. [PMID: 29758729 DOI: 10.1103/physreve.97.042418] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2018] [Indexed: 06/08/2023]
Abstract
Neural field theory is used to predict and analyze the phenomenon of perceptual echo in which random input stimuli at one location are correlated with electroencephalographic responses at other locations. It is shown that this echo correlation (EC) yields an estimate of the transfer function from the stimulated point to other locations. Modal analysis then explains the observed spatiotemporal structure of visually driven EC and the dominance of the alpha frequency; two eigenmodes of similar amplitude dominate the response, leading to temporal beating and a line of low correlation that runs from the crown of the head toward the ears. These effects result from mode splitting and symmetry breaking caused by interhemispheric coupling and cortical folding. It is shown how eigenmodes obtained from functional magnetic resonance imaging experiments can be combined with temporal dynamics from EC or other evoked responses to estimate the spatiotemporal transfer function between any two points and hence their effective connectivity.
Collapse
Affiliation(s)
- P A Robinson
- School of Physics, University of Sydney, NSW 2006, Australia and Center for Integrative Brain Function, University of Sydney, NSW 2006, Australia
| | - J C Pagès
- School of Physics, University of Sydney, NSW 2006, Australia and Center for Integrative Brain Function, University of Sydney, NSW 2006, Australia
| | - N C Gabay
- School of Physics, University of Sydney, NSW 2006, Australia and Center for Integrative Brain Function, University of Sydney, NSW 2006, Australia
| | - T Babaie
- School of Physics, University of Sydney, NSW 2006, Australia and Center for Integrative Brain Function, University of Sydney, NSW 2006, Australia
| | - K N Mukta
- School of Physics, University of Sydney, NSW 2006, Australia and Center for Integrative Brain Function, University of Sydney, NSW 2006, Australia
| |
Collapse
|
29
|
Abstract
Transfers of large-scale neural activity into, within and between corticothalamic neural populations and brain hemispheres are analysed using time-integrated transfer functions and state parameters obtained from neural field theory for a variety of arousal states. It is shown that the great majority of activity results from feedbacks within the corticothalamic system, including significant transfer between hemispheres, but only a small minority arises via net input from the external world, with the brain thus in a near-critical, highly introspective state. Notably, the total excitatory and inhibitory influences on cortical neurons are balanced to within a few per cent across arousal states. Strong negative intrahemispheric feedforward exists to the cortex, and even larger interhemispheric positive feedforward, but these are modified by feedback loops to yield near-critical positive overall gain. The results underline the utility of transfer functions for the analysis of brain activity.
Collapse
Affiliation(s)
- P A Robinson
- School of Physics, University of Sydney, Sydney, New South Wales 2006, Australia .,Center for Integrative Brain Function, University of Sydney, Sydney, New South Wales 2006, Australia
| |
Collapse
|
30
|
Mukta KN, MacLaurin JN, Robinson PA. Theory of corticothalamic brain activity in a spherical geometry: Spectra, coherence, and correlation. Phys Rev E 2017; 96:052410. [PMID: 29347754 DOI: 10.1103/physreve.96.052410] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2017] [Indexed: 11/07/2022]
Abstract
Corticothalamic neural field theory is applied to a spherical geometry to better model neural activity in the human brain and is also compared with planar approximations. The frequency power spectrum, correlation, and coherence functions are computed analytically and numerically. The effects of cortical boundary conditions and resulting modal aspects of spherical corticothalamic dynamics are explored, showing that the results of spherical and finite planar geometries converge to those for the infinite planar geometry in the limit of large brain size. Estimates are made of the point at which modal series can be truncated and it is found that for physiologically plausible parameters only the lowest few spatial eigenmodes are needed for an accurate representation of macroscopic brain activity. A difference between the geometries is that there is a low-frequency 1/f spectrum in the infinite planar geometry, whereas in the spherical geometry it is 1/f^{2}. Another difference is that the alpha peak in the spherical geometry is sharper and stronger than in the planar geometry. Cortical modal effects can lead to a double alpha peak structure in the power spectrum, although the main determinant of the alpha peak is corticothalamic feedback. In the spherical geometry, the cross spectrum between two points is found to only depend on their relative distance apart. At small spatial separations the low-frequency cross spectrum is stronger than for an infinite planar geometry and the alpha peak is sharper and stronger due to the partitioning of the energy into discrete modes. In the spherical geometry, the coherence function between points decays monotonically as their separation increases at a fixed frequency, but persists further at resonant frequencies. The correlation between two points is found to be positive, regardless of the time lag and spatial separation, but decays monotonically as the separation increases at fixed time lag. At fixed distance the correlation has peaks at multiples of the period of the dominant frequency of system activity.
Collapse
Affiliation(s)
- K N Mukta
- School of Physics, University of Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - J N MacLaurin
- School of Physics, University of Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - P A Robinson
- School of Physics, University of Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| |
Collapse
|
31
|
Roy N, Sanz-Leon P, Robinson PA. Spectral signatures of activity-dependent neural feedback in the corticothalamic system. Phys Rev E 2017; 96:052310. [PMID: 29347805 DOI: 10.1103/physreve.96.052310] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2017] [Indexed: 11/07/2022]
Abstract
The modulation of neural quantities by presynaptic and postsynaptic activities via local feedback processes is investigated by incorporating nonlinear phenomena such as relative refractory period, synaptic enhancement, synaptic depression, and habituation. This is done by introducing susceptibilities, which quantify the response in either firing threshold or synaptic strength to unit change in either presynaptic or postsynaptic activity. Effects on the power spectra are then analyzed for a realistic corticothalamic model to determine the spectral signatures of various nonlinear processes and to what extent these are distinct. Depending on the feedback processes, there can be enhancements or reductions in low-frequency and/or alpha power, splitting of the alpha resonance, and/or appearance of new resonances at high frequencies. These features in the power spectra allow processes to be fully distinguished where they are unique, or partly distinguished if they are common to only a subset of feedbacks, and can potentially be used to constrain the types, strengths, and dynamics of feedbacks present.
Collapse
Affiliation(s)
- N Roy
- School of Physics, University of Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - P Sanz-Leon
- School of Physics, University of Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - P A Robinson
- School of Physics, University of Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| |
Collapse
|
32
|
Lacy TC, Aquino KM, Robinson PA, Schira MM. Shock-like haemodynamic responses induced in the primary visual cortex by moving visual stimuli. J R Soc Interface 2017; 13:rsif.2016.0576. [PMID: 27974572 DOI: 10.1098/rsif.2016.0576] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2016] [Accepted: 11/17/2016] [Indexed: 11/12/2022] Open
Abstract
It is shown that recently discovered haemodynamic waves can form shock-like fronts when driven by stimuli that excite the cortex in a patch that moves faster than the haemodynamic wave velocity. If stimuli are chosen in order to induce shock-like behaviour, the resulting blood oxygen level-dependent (BOLD) response is enhanced, thereby improving the signal to noise ratio of measurements made with functional magnetic resonance imaging. A spatio-temporal haemodynamic model is extended to calculate the BOLD response and determine the main properties of waves induced by moving stimuli. From this, the optimal conditions for stimulating shock-like responses are determined, and ways of inducing these responses in experiments are demonstrated in a pilot study.
Collapse
Affiliation(s)
- T C Lacy
- School of Physics, University of Sydney, New South Wales, Australia .,Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - K M Aquino
- School of Physics, University of Sydney, New South Wales, Australia.,Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia.,Sir Peter Mansfield Imaging Centre, School of Physics and Astronomy, University of Nottingham, University Park, Nottingham NG7 2RD, UK
| | - P A Robinson
- School of Physics, University of Sydney, New South Wales, Australia.,Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - M M Schira
- School of Psychology, University of Wollongong, Wollongong, New South Wales 2522, Australia.,Neuroscience Research Australia, Royal Hospital for Women, Randwick, New South Wales 2031, Australia
| |
Collapse
|
33
|
Pang JC, Robinson PA, Aquino KM. Response-mode decomposition of spatio-temporal haemodynamics. J R Soc Interface 2017; 13:rsif.2016.0253. [PMID: 27170653 DOI: 10.1098/rsif.2016.0253] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2016] [Accepted: 04/12/2016] [Indexed: 12/12/2022] Open
Abstract
The blood oxygen-level dependent (BOLD) response to a neural stimulus is analysed using the transfer function derived from a physiologically based poroelastic model of cortical tissue. The transfer function is decomposed into components that correspond to distinct poles, each related to a response mode with a natural frequency and dispersion relation; together these yield the total BOLD response. The properties of the decomposed components provide a deeper understanding of the nature of the BOLD response, via the components' frequency dependences, spatial and temporal power spectra, and resonances. The transfer function components are then used to separate the BOLD response to a localized impulse stimulus, termed the Green function or spatio-temporal haemodynamic response function, into component responses that are explicitly related to underlying physiological quantities. The analytical results also provide a quantitative tool to calculate the linear BOLD response to an arbitrary neural drive, which is faster to implement than direct Fourier transform methods. The results of this study can be used to interpret functional magnetic resonance imaging data in new ways based on physiology, to enhance deconvolution methods and to design experimental protocols that can selectively enhance or suppress particular responses, to probe specific physiological phenomena.
Collapse
Affiliation(s)
- J C Pang
- School of Physics, University of Sydney, Sydney, New South Wales 2006, Australia Center for Integrative Brain Function, University of Sydney, Sydney, New South Wales 2006, Australia
| | - P A Robinson
- School of Physics, University of Sydney, Sydney, New South Wales 2006, Australia Center for Integrative Brain Function, University of Sydney, Sydney, New South Wales 2006, Australia
| | - K M Aquino
- School of Physics, University of Sydney, Sydney, New South Wales 2006, Australia Sir Peter Mansfield Imaging Centre, School of Physics and Astronomy, University of Nottingham, University Park, Nottingham NG7 2RD, UK
| |
Collapse
|
34
|
Abstract
Perturbation analysis of neural field theory is used to derive eigenmodes of neural activity on a cortical hemisphere, which have previously been calculated numerically and found to be close analogs of spherical harmonics, despite heavy cortical folding. The present perturbation method treats cortical folding as a first-order perturbation from a spherical geometry. The first nine spatial eigenmodes on a population-averaged cortical hemisphere are derived and compared with previous numerical solutions. These eigenmodes contribute most to brain activity patterns such as those seen in electroencephalography and functional magnetic resonance imaging. The eigenvalues of these eigenmodes are found to agree with the previous numerical solutions to within their uncertainties. Also in agreement with the previous numerics, all eigenmodes are found to closely resemble spherical harmonics. The first seven eigenmodes exhibit a one-to-one correspondence with their numerical counterparts, with overlaps that are close to unity. The next two eigenmodes overlap the corresponding pair of numerical eigenmodes, having been rotated within the subspace spanned by that pair, likely due to second-order effects. The spatial orientations of the eigenmodes are found to be fixed by gross cortical shape rather than finer-scale cortical properties, which is consistent with the observed intersubject consistency of functional connectivity patterns. However, the eigenvalues depend more sensitively on finer-scale cortical structure, implying that the eigenfrequencies and consequent dynamical properties of functional connectivity depend more strongly on details of individual cortical folding. Overall, these results imply that well-established tools from perturbation theory and spherical harmonic analysis can be used to calculate the main properties and dynamics of low-order brain eigenmodes.
Collapse
Affiliation(s)
- Natasha C Gabay
- School of Physics, University of Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - P A Robinson
- School of Physics, University of Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| |
Collapse
|
35
|
Abstract
Neural field theory of the corticothalamic system is used to analyze the properties of its steady-state solutions, including their linear stability, in the parameter space of synaptic couplings for physiological parameter ranges representing normal arousal waking states in adult humans. The independent connections of the corticothalamic model define an eight-dimensional parameter space, while specific combinations of these connections parameterize intracortical, corticothalamic, and intrathalamic loops. Multistable regions are systematically identified and the existence of up to five steady-state solutions is confirmed, up to three of which are linearly stable. A key determinant for the existence of five steady states is found to be the number of nonzero connections. This finding had not been previously proposed as the determining factor of high multiplicities of multistability in mesoscopic models of the brain. In the corticothalamic model presented here, multistability occurs when the intrathalamic loop is present (i.e., the reticular nucleus inhibits the relay nuclei), and when the net synaptic effect of the intracortical loop is inhibitory. The signature of these additional waking states is an overall increased level of thalamic activity. It is argued that the additional steady states found may represent hyperarousal states which occur when the corticothalamic projections do not attenuate the activity of the cortex.
Collapse
Affiliation(s)
- Paula Sanz-Leon
- School of Physics, University of Sydney, NSW 2006, Australia; Center for Integrative Brain Function, University of Sydney, NSW 2006, Australia.
| | - P A Robinson
- School of Physics, University of Sydney, NSW 2006, Australia; Center for Integrative Brain Function, University of Sydney, NSW 2006, Australia
| |
Collapse
|
36
|
Rubchinsky LL, Ahn S, Klijn W, Cumming B, Yates S, Karakasis V, Peyser A, Woodman M, Diaz-Pier S, Deraeve J, Vassena E, Alexander W, Beeman D, Kudela P, Boatman-Reich D, Anderson WS, Luque NR, Naveros F, Carrillo RR, Ros E, Arleo A, Huth J, Ichinose K, Park J, Kawai Y, Suzuki J, Mori H, Asada M, Oprisan SA, Dave AI, Babaie T, Robinson P, Tabas A, Andermann M, Rupp A, Balaguer-Ballester E, Lindén H, Christensen RK, Nakamura M, Barkat TR, Tosi Z, Beggs J, Lonardoni D, Boi F, Di Marco S, Maccione A, Berdondini L, Jędrzejewska-Szmek J, Dorman DB, Blackwell KT, Bauermeister C, Keren H, Braun J, Dornas JV, Mavritsaki E, Aldrovandi S, Bridger E, Lim S, Brunel N, Buchin A, Kerr CC, Chizhov A, Huberfeld G, Miles R, Gutkin B, Spencer MJ, Meffin H, Grayden DB, Burkitt AN, Davey CE, Tao L, Tiruvadi V, Ali R, Mayberg H, Butera R, Gunay C, Lamb D, Calabrese RL, Doloc-Mihu A, López-Madrona VJ, Matias FS, Pereda E, Mirasso CR, Canals S, Geminiani A, Pedrocchi A, D’Angelo E, Casellato C, Chauhan A, Soman K, Srinivasa Chakravarthy V, Muddapu VR, Chuang CC, Chen NY, Bayati M, Melchior J, Wiskott L, Azizi AH, Diba K, Cheng S, Smirnova EY, Yakimova EG, Chizhov AV, Chen NY, Shih CT, Florescu D, Coca D, Courtiol J, Jirsa VK, Covolan RJM, Teleńczuk B, Kempter R, Curio G, Destexhe A, Parker J, Klishko AN, Prilutsky BI, Cymbalyuk G, Franke F, Hierlemann A, da Silveira RA, Casali S, Masoli S, Rizza M, Rizza MF, Masoli S, Sun Y, Wong W, Farzan F, Blumberger DM, Daskalakis ZJ, Popovych S, Viswanathan S, Rosjat N, Grefkes C, Daun S, Gentiletti D, Suffczynski P, Gnatkovski V, De Curtis M, Lee H, Paik SB, Choi W, Jang J, Park Y, Song JH, Song M, Pallarés V, Gilson M, Kühn S, Insabato A, Deco G, Glomb K, Ponce-Alvarez A, Ritter P, Gilson M, Campo AT, Thiele A, Deeba F, Robinson PA, van Albada SJ, Rowley A, Hopkins M, Schmidt M, Stokes AB, Lester DR, Furber S, Diesmann M, Barri A, Wiechert MT, DiGregorio DA, Dimitrov AG, Vich C, Berg RW, Guillamon A, Ditlevsen S, Cazé RD, Girard B, Doncieux S, Doyon N, Boahen F, Desrosiers P, Laurence E, Doyon N, Dubé LJ, Eleonora R, Durstewitz D, Schmidt D, Mäki-Marttunen T, Krull F, Bettella F, Metzner C, Devor A, Djurovic S, Dale AM, Andreassen OA, Einevoll GT, Næss S, Ness TV, Halnes G, Halgren E, Halnes G, Mäki-Marttunen T, Pettersen KH, Andreassen OA, Sætra MJ, Hagen E, Schiffer A, Grzymisch A, Persike M, Ernst U, Harnack D, Ernst UA, Tomen N, Zucca S, Pasquale V, Pica G, Molano-Mazón M, Chiappalone M, Panzeri S, Fellin T, Oie KS, Boothe DL, Crone JC, Yu AB, Felton MA, Zulfiqar I, Moerel M, De Weerd P, Formisano E, Boothe DL, Crone JC, Felton MA, Oie K, Franaszczuk P, Diggelmann R, Fiscella M, Hierlemann A, Franke F, Guarino D, Antolík J, Davison AP, Frègnac Y, Etienne BX, Frohlich F, Lefebvre J, Marcos E, Mattia M, Genovesio A, Fedorov LA, Dijkstra TM, Sting L, Hock H, Giese MA, Buhry L, Langlet C, Giovannini F, Verbist C, Salvadé S, Giugliano M, Henderson JA, Wernecke H, Sándor B, Gros C, Voges N, Dabrovska P, Riehle A, Brochier T, Grün S, Gu Y, Gong P, Dumont G, Novikov NA, Gutkin BS, Tewatia P, Eriksson O, Kramer A, Santos J, Jauhiainen A, Kotaleski JH, Belić JJ, Kumar A, Kotaleski JH, Shimono M, Hatano N, Ahmad S, Cui Y, Hawkins J, Senk J, Korvasová K, Tetzlaff T, Helias M, Kühn T, Denker M, Mana P, Grün S, Dahmen D, Schuecker J, Goedeke S, Keup C, Goedeke S, Heuer K, Bakker R, Tiesinga P, Toro R, Qin W, Hadjinicolaou A, Grayden DB, Ibbotson MR, Kameneva T, Lytton WW, Mulugeta L, Drach A, Myers JG, Horner M, Vadigepalli R, Morrison T, Walton M, Steele M, Anthony Hunt C, Tam N, Amaducci R, Muñiz C, Reyes-Sánchez M, Rodríguez FB, Varona P, Cronin JT, Hennig MH, Iavarone E, Yi J, Shi Y, Zandt BJ, Van Geit W, Rössert C, Markram H, Hill S, O’Reilly C, Iavarone E, Shi Y, Perin R, Lu H, Zandt BJ, Bryson A, Rössert C, Hadrava M, Hlinka J, Hosaka R, Olenik M, Houghton C, Iannella N, Launey T, Kameneva T, Kotsakidis R, Meffin H, Soriano J, Kubo T, Inoue T, Kida H, Yamakawa T, Suzuki M, Ikeda K, Abbasi S, Hudson AE, Heck DH, Jaeger D, Lee J, Abbasi S, Janušonis S, Saggio ML, Spiegler A, Stacey WC, Bernard C, Lillo D, Bernard C, Petkoski S, Spiegler A, Drakesmith M, Jones DK, Zadeh AS, Kambhampati C, Karbowski J, Kaya ZG, Lakretz Y, Treves A, Li LW, Lizier J, Kerr CC, Masquelier T, Kheradpisheh SR, Kim H, Kim CS, Marakshina JA, Vartanov AV, Neklyudova AA, Kozlovskiy SA, Kiselnikov AA, Taniguchi K, Kitano K, Schmitt O, Lessmann F, Schwanke S, Eipert P, Meinhardt J, Beier J, Kadir K, Karnitzki A, Sellner L, Klünker AC, Kuch L, Ruß F, Jenssen J, Wree A, Sanz-Leon P, Knock SA, Chien SC, Maess B, Knösche TR, Cohen CC, Popovic MA, Klooster J, Kole MH, Roberts EA, Kopell NJ, Kepple D, Giaffar H, Rinberg D, Koulakov A, Forlim CG, Klock L, Bächle J, Stoll L, Giemsa P, Fuchs M, Schoofs N, Montag C, Gallinat J, Lee RX, Stephens GJ, Kuhn B, Tauffer L, Isope P, Inoue K, Ohmura Y, Yonekura S, Kuniyoshi Y, Jang HJ, Kwag J, de Kamps M, Lai YM, dos Santos F, Lam KP, Andras P, Imperatore J, Helms J, Tompa T, Lavin A, Inkpen FH, Ashby MC, Lepora NF, Shifman AR, Lewis JE, Zhang Z, Feng Y, Tetzlaff C, Kulvicius T, Li Y, Pena RFO, Bernardi D, Roque AC, Lindner B, Bernardi D, Vellmer S, Saudargiene A, Maninen T, Havela R, Linne ML, Powanwe A, Longtin A, Naveros F, Garrido JA, Graham JW, Dura-Bernal S, Angulo SL, Neymotin SA, Antic SD. 26th Annual Computational Neuroscience Meeting (CNS*2017): Part 2. BMC Neurosci 2017. [PMCID: PMC5592442 DOI: 10.1186/s12868-017-0371-2] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
37
|
Müller EJ, van Albada SJ, Kim JW, Robinson PA. Unified neural field theory of brain dynamics underlying oscillations in Parkinson's disease and generalized epilepsies. J Theor Biol 2017. [PMID: 28633970 DOI: 10.1016/j.jtbi.2017.06.016] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The mechanisms underlying pathologically synchronized neural oscillations in Parkinson's disease (PD) and generalized epilepsies are explored in parallel via a physiologically-based neural field model of the corticothalamic-basal ganglia (CTBG) system. The basal ganglia (BG) are approximated as a single effective population and their roles in the modulation of oscillatory dynamics of the corticothalamic (CT) system and vice versa are analyzed. In addition to normal EEG rhythms, enhanced activity around 4 Hz and 20 Hz exists in the model, consistent with the characteristic frequencies observed in PD. These rhythms result from resonances in loops formed between the BG and CT populations, analogous to those that underlie epileptic oscillations in a previous CT model, and which are still present in the combined CTBG system. Dopamine depletion is argued to weaken the dampening of these loop resonances in PD, and network connections then explain the significant coherence observed between BG, thalamic, and cortical population activity around 4-8 Hz and 20 Hz. Parallels between the afferent and efferent connection sites of the thalamic reticular nucleus (TRN) and BG predict low dopamine to correspond to a reduced likelihood of tonic-clonic (grand mal) seizures, which agrees with experimental findings. Furthermore, the model predicts an increased likelihood of absence (petit mal) seizure resulting from pathologically low dopamine levels in accordance with experimental observations. Suppression of absence seizure activity is demonstrated when afferent and efferent BG connections to the CT system are strengthened, which is consistent with other CTBG modeling studies. The BG are demonstrated to have a suppressive effect on activity of the CTBG system near tonic-clonic seizure states, which provides insight into the reported efficacy of current treatments in BG circuits. Sleep states of the TRN are also found to suppress pathological PD activity in accordance with observations. Overall, the findings demonstrate strong parallels between coherent oscillations in generalized epilepsies and PD, and provide insights into possible comorbidities.
Collapse
Affiliation(s)
- E J Müller
- School of Physics, The University of Sydney, Sydney, NSW 2006, Australia; Center for Integrative Brain Function, The University of Sydney, NSW 2006, Australia.
| | - S J van Albada
- School of Physics, The University of Sydney, Sydney, NSW 2006, Australia; Institute of Neuroscience and Medicine (INM-6) and Institute for Advanced Simulation (IAS-6) and JARA BRAIN Institute I, Jülich Research Center, Jülich, Germany
| | - J W Kim
- School of Physics, The University of Sydney, Sydney, NSW 2006, Australia; Center for Integrative Brain Function, The University of Sydney, NSW 2006, Australia
| | - P A Robinson
- School of Physics, The University of Sydney, Sydney, NSW 2006, Australia; Center for Integrative Brain Function, The University of Sydney, NSW 2006, Australia
| |
Collapse
|
38
|
Yang DP, Robinson PA. Critical dynamics of Hopf bifurcations in the corticothalamic system: Transitions from normal arousal states to epileptic seizures. Phys Rev E 2017; 95:042410. [PMID: 28505725 DOI: 10.1103/physreve.95.042410] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2016] [Indexed: 06/07/2023]
Abstract
A physiologically based corticothalamic model of large-scale brain activity is used to analyze critical dynamics of transitions from normal arousal states to epileptic seizures, which correspond to Hopf bifurcations. This relates an abstract normal form quantitatively to underlying physiology that includes neural dynamics, axonal propagation, and time delays. Thus, a bridge is constructed that enables normal forms to be used to interpret quantitative data. The normal form of the Hopf bifurcations with delays is derived using Hale's theory, the center manifold theorem, and normal form analysis, and it is found to be explicitly expressed in terms of transfer functions and the sensitivity matrix of a reduced open-loop system. It can be applied to understand the effect of each physiological parameter on the critical dynamics and determine whether the Hopf bifurcation is supercritical or subcritical in instabilities that lead to absence and tonic-clonic seizures. Furthermore, the effects of thalamic and cortical nonlinearities on the bifurcation type are investigated, with implications for the roles of underlying physiology. The theoretical predictions about the bifurcation type and the onset dynamics are confirmed by numerical simulations and provide physiologically based criteria for determining bifurcation types from first principles. The results are consistent with experimental data from previous studies, imply that new regimes of seizure transitions may exist in clinical settings, and provide a simplified basis for control-systems interventions. Using the normal form, and the full equations from which it is derived, more complex dynamics, such as quasiperiodic cycles and saddle cycles, are discovered near the critical points of the subcritical Hopf bifurcations.
Collapse
Affiliation(s)
- Dong-Ping Yang
- School of Physics, University of Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - P A Robinson
- School of Physics, University of Sydney, New South Wales 2006, Australia and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| |
Collapse
|
39
|
Zobaer MS, Anderson RM, Kerr CC, Robinson PA, Wong KKH, D'Rozario AL. K-complexes, spindles, and ERPs as impulse responses: unification via neural field theory. Biol Cybern 2017; 111:149-164. [PMID: 28251306 DOI: 10.1007/s00422-017-0713-2] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/04/2016] [Accepted: 02/12/2017] [Indexed: 06/06/2023]
Abstract
To interrelate K-complexes, spindles, evoked response potentials (ERPs), and spontaneous electroencephalography (EEG) using neural field theory (NFT), physiology-based NFT of the corticothalamic system is used to model cortical excitatory and inhibitory populations and thalamic relay and reticular nuclei. The impulse response function of the model is used to predict the responses to impulses, which are compared with transient waveforms in sleep studies. Fits to empirical data then allow underlying brain physiology to be inferred and compared with other waves. Spontaneous K-complexes, spindles, and other transient waveforms can be reproduced using NFT by treating them as evoked responses to impulsive stimuli with brain parameters appropriate to spontaneous EEG in sleep stage 2. Using this approach, spontaneous K-complexes and sleep spindles can be analyzed using the same single theory as previously been used to account for waking ERPs and other EEG phenomena. As a result, NFT can explain a wide variety of transient waveforms that have only been phenomenologically classified to date. This enables noninvasive fitting to be used to infer underlying physiological parameters. This physiology-based model reproduces the time series of different transient EEG waveforms; it has previously reproduced experimental EEG spectra, and waking ERPs, and many other observations, thereby unifying transient sleep waveforms with these phenomena.
Collapse
Affiliation(s)
- M S Zobaer
- School of Physics, The University of Sydney, Sydney, NSW, 2006, Australia.
- Center for Integrative Brain Function, University of Sydney, Sydney, NSW, 2006, Australia.
- Center for Research Excellence, Neurosleep, 431 Glebe Point Rd, Glebe, NSW, 2037, Australia.
- Department of Physics, Bangladesh University of Textiles, Dhaka, 1208, Bangladesh.
| | - R M Anderson
- School of Physics, The University of Sydney, Sydney, NSW, 2006, Australia
| | - C C Kerr
- School of Physics, The University of Sydney, Sydney, NSW, 2006, Australia
- Center for Integrative Brain Function, University of Sydney, Sydney, NSW, 2006, Australia
- Department of Physiology and Pharmacology, State University of New York Downstate Medical Center, 450 Clarkson Ave, Brooklyn, NY, USA
| | - P A Robinson
- School of Physics, The University of Sydney, Sydney, NSW, 2006, Australia
- Center for Integrative Brain Function, University of Sydney, Sydney, NSW, 2006, Australia
- Center for Research Excellence, Neurosleep, 431 Glebe Point Rd, Glebe, NSW, 2037, Australia
| | - K K H Wong
- CIRUS, Woolcock Institute of Medical Research, The University of Sydney, Sydney, NSW, Australia
- Respiratory and Sleep Disorders Department, Royal Prince Alfred Hospital and Sydney Local Health District, Sydney, NSW, Australia
- Sydney Medical School, The University of Sydney, Sydney, NSW, Australia
| | - A L D'Rozario
- CIRUS, Woolcock Institute of Medical Research, The University of Sydney, Sydney, NSW, Australia
- Respiratory and Sleep Disorders Department, Royal Prince Alfred Hospital and Sydney Local Health District, Sydney, NSW, Australia
- School of Psychology, Faculty of Science, Brain and Mind Centre and Charles Perkins Centre, The University of Sydney, Sydney, NSW, Australia
| |
Collapse
|
40
|
Mehta-Pandejee G, Robinson PA, Henderson JA, Aquino KM, Sarkar S. Inference of direct and multistep effective connectivities from functional connectivity of the brain and of relationships to cortical geometry. J Neurosci Methods 2017; 283:42-54. [PMID: 28342831 DOI: 10.1016/j.jneumeth.2017.03.014] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2016] [Revised: 03/15/2017] [Accepted: 03/18/2017] [Indexed: 01/26/2023]
Abstract
BACKGROUND The problem of inferring effective brain connectivity from functional connectivity is under active investigation, and connectivity via multistep paths is poorly understood. NEW METHOD A method is presented to calculate the direct effective connection matrix (deCM), which embodies direct connection strengths between brain regions, from functional CMs (fCMs) by minimizing the difference between an experimental fCM and one calculated via neural field theory from an ansatz deCM based on an experimental anatomical CM. RESULTS The best match between fCMs occurs close to a critical point, consistent with independent published stability estimates. Residual mismatch between fCMs is identified to be largely due to interhemispheric connections that are poorly estimated in an initial ansatz deCM due to experimental limitations; improved ansatzes substantially reduce the mismatch and enable interhemispheric connections to be estimated. Various levels of significant multistep connections are then imaged via the neural field theory (NFT) result that these correspond to powers of the deCM; these are shown to be predictable from geometric distances between regions. COMPARISON WITH EXISTING METHODS This method gives insight into direct and multistep effective connectivity from fCMs and relating to physiology and brain geometry. This contrasts with other methods, which progressively adjust connections without an overarching physiologically based framework to deal with multistep or poorly estimated connections. CONCLUSIONS deCMs can be usefully estimated using this method and the results enable multistep connections to be investigated systematically.
Collapse
Affiliation(s)
- Grishma Mehta-Pandejee
- School of Physics, The University of Sydney, Sydney, New South Wales 2006, Australia; Center of Excellence for Integrative Brain Function, The University of Sydney, New South Wales 2006, Australia.
| | - P A Robinson
- School of Physics, The University of Sydney, Sydney, New South Wales 2006, Australia; Center of Excellence for Integrative Brain Function, The University of Sydney, New South Wales 2006, Australia
| | - James A Henderson
- School of Physics, The University of Sydney, Sydney, New South Wales 2006, Australia; Center of Excellence for Integrative Brain Function, The University of Sydney, New South Wales 2006, Australia; School of Information Technology and Electrical Engineering, University of Queensland, St Lucia, Queensland 4072, Australia
| | - K M Aquino
- School of Physics, The University of Sydney, Sydney, New South Wales 2006, Australia; Center of Excellence for Integrative Brain Function, The University of Sydney, New South Wales 2006, Australia; Sir Peter Mansfield Imaging Center, University of Nottingham, Nottingham NG7 2RD, United Kingdom
| | - Somwrita Sarkar
- Center of Excellence for Integrative Brain Function, The University of Sydney, New South Wales 2006, Australia; Design Lab, University of Sydney, Sydney, New South Wales 2006, Australia
| |
Collapse
|
41
|
Abstract
A recent physiologically based model of the ascending arousal system is used to analyze the dynamics near the transition from wake to sleep, which corresponds to a saddle-node bifurcation at a critical point. A normal form is derived by approximating the dynamics by those of a particle in a parabolic potential well with dissipation. This mechanical analog is used to calculate the power spectrum of fluctuations in response to a white noise drive, and the scalings of fluctuation variance and spectral width are derived versus distance from the critical point. The predicted scalings are quantitatively confirmed by numerical simulations, which show that the variance increases and the spectrum undergoes critical slowing, both in accord with theory. These signals can thus serve as potential precursors to indicate imminent wake-sleep transition, with potential application to safety-critical occupations in transport, air-traffic control, medicine, and heavy industry.
Collapse
Affiliation(s)
- Dong-Ping Yang
- School of Physics, University of Sydney, New South Wales 2006, Australia.,Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | | | - Angela Karanjai
- School of Physics, University of Sydney, New South Wales 2006, Australia
| | - P A Robinson
- School of Physics, University of Sydney, New South Wales 2006, Australia.,Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| |
Collapse
|
42
|
Sarkar S, Chawla S, Robinson PA, Fortunato S. Publisher's Note: Eigenvector dynamics under perturbation of modular networks [Phys. Rev. E 93, 062312 (2016)]. Phys Rev E 2016; 94:019903. [PMID: 27575238 DOI: 10.1103/physreve.94.019903] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2016] [Indexed: 06/06/2023]
Abstract
This corrects the article DOI: 10.1103/PhysRevE.93.062312.
Collapse
|
43
|
Abstract
A quantitative, physiology-based model of the ascending arousal system is developed, using continuum neuronal population modeling, which involves averaging properties such as firing rates across neurons in each population. The model includes the ventrolateral preoptic area (VLPO), where circadian and homeostatic drives enter the system, the monoaminergic and cholinergic nuclei of the ascending arousal system, and their interconnections. The human sleep-wake cycle is governed by the activities of these nuclei, which modulate the behavioral state of the brain via diffuse neuromodulatory projections. The model parameters are not free since they correspond to physiological observables. Approximate parameter bounds are obtained by requiring consistency with physiological and behavioral measures, and the model replicates the human sleep-wake cycle, with physiologically reasonable voltages and firing rates. Mutual inhibition between the wake-promoting monoaminergic group and sleep-promoting VLPO causes ``flip-flop'' behavior, with most time spent in 2 stable steady states corresponding to wake and sleep, with transitions between them on a timescale of a few minutes. The model predicts hysteresis in the sleep-wake cycle, with a region of bistability of the wake and sleep states. Reducing the monoaminergic-VLPO mutual inhibition results in a smaller hysteresis loop. This makes the model more prone to wake-sleep transitions in both directions and makes the states less distinguishable, as in narcolepsy. The model behavior is robust across the constrained parameter ranges, but with sufficient flexibility to describe a wide range of observed phenomena.
Collapse
Affiliation(s)
- A J K Phillips
- School of Physics, University of Sydney, New South Wales 2006, Australia.
| | | |
Collapse
|
44
|
Wilson MT, Fung PK, Robinson PA, Shemmell J, Reynolds JNJ. Calcium dependent plasticity applied to repetitive transcranial magnetic stimulation with a neural field model. J Comput Neurosci 2016; 41:107-25. [DOI: 10.1007/s10827-016-0607-7] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2015] [Revised: 05/05/2016] [Accepted: 05/12/2016] [Indexed: 10/21/2022]
|
45
|
Sarkar S, Chawla S, Robinson PA, Fortunato S. Eigenvector dynamics under perturbation of modular networks. Phys Rev E 2016; 93:062312. [PMID: 27415285 DOI: 10.1103/physreve.93.062312] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2015] [Indexed: 06/06/2023]
Abstract
Rotation dynamics of eigenvectors of modular network adjacency matrices under random perturbations are presented. In the presence of q communities, the number of eigenvectors corresponding to the q largest eigenvalues form a "community" eigenspace and rotate together, but separately from that of the "bulk" eigenspace spanned by all the other eigenvectors. Using this property, the number of modules or clusters in a network can be estimated in an algorithm-independent way. A general argument and derivation for the theoretical detectability limit for sparse modular networks with q communities is presented, beyond which modularity persists in the system but cannot be detected. It is shown that for detecting the clusters or modules using the adjacency matrix, there is a "band" in which it is hard to detect the clusters even before the theoretical detectability limit is reached, and for which the theoretically predicted detectability limit forms the sufficient upper bound. Analytic estimations of these bounds are presented and empirically demonstrated.
Collapse
Affiliation(s)
- Somwrita Sarkar
- Design Lab, Faculty of Architecture, Design and Planning, University of Sydney, Australia NSW 2006 and ARC Centre of Excellence for Integrative Brain Function
| | - Sanjay Chawla
- Qatar Computing Research Institute, Hamad Bin Khalifa University (HBKU), Qatar Foundation, Doha, Qatar and Faculty of Engineering and IT, University of Sydney, Australia NSW 2006
| | - P A Robinson
- School of Physics, University of Sydney, Australia NSW 2006 and ARC Centre of Excellence for Integrative Brain Function
| | | |
Collapse
|
46
|
Robinson PA, Zhao X, Aquino KM, Griffiths JD, Sarkar S, Mehta-Pandejee G. Eigenmodes of brain activity: Neural field theory predictions and comparison with experiment. Neuroimage 2016; 142:79-98. [PMID: 27157788 DOI: 10.1016/j.neuroimage.2016.04.050] [Citation(s) in RCA: 70] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2015] [Revised: 03/13/2016] [Accepted: 04/21/2016] [Indexed: 12/20/2022] Open
Abstract
Neural field theory of the corticothalamic system is applied to predict and analyze the activity eigenmodes of the bihemispheric brain, focusing particularly on their spatial structure. The eigenmodes of a single brain hemisphere are found to be close analogs of spherical harmonics, which are the natural modes of the sphere. Instead of multiple eigenvalues being equal, as in the spherical case, cortical folding splits them to have distinct values. Inclusion of interhemispheric connections between homologous regions via the corpus callosum leads to further splitting that depends on symmetry or antisymmetry of activity between brain hemispheres, and the strength and sign of the interhemispheric connections. Symmetry properties of the lowest observed eigenmodes strongly constrain the interhemispheric connectivity strengths and unihemispheric mode spectra, and it is predicted that most spontaneous brain activity will be symmetric between hemispheres, consistent with observations. Comparison with the eigenmodes of an experimental anatomical connectivity matrix confirms these results, permits the relative strengths of intrahemispheric and interhemispheric connectivities to be approximately inferred from their eigenvalues, and lays the foundation for further experimental tests. The results are consistent with brain activity being in corticothalamic eigenmodes, rather than discrete "networks" and open the way to new approaches to brain analysis.
Collapse
Affiliation(s)
- P A Robinson
- School of Physics, University of Sydney, New South Wales 2006, Australia; Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia.
| | - X Zhao
- School of Physics, University of Sydney, New South Wales 2006, Australia; Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - K M Aquino
- School of Physics, University of Sydney, New South Wales 2006, Australia; Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia; Sir Peter Mansfield Imaging Center, University of Nottingham, Nottingham NG7 2RD, UK, EU
| | - J D Griffiths
- School of Physics, University of Sydney, New South Wales 2006, Australia; Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia; Rotman Research Institute at Baycrest, 3560 Bathurst St, Toronto, Ontario, M6A 2E1, Canada
| | - S Sarkar
- Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia; Design Lab, School of Architecture, Design, and Planning, University of Sydney, New South Wales 2006, Australia
| | - Grishma Mehta-Pandejee
- School of Physics, University of Sydney, New South Wales 2006, Australia; Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| |
Collapse
|
47
|
Abeysuriya RG, Robinson PA. Real-time automated EEG tracking of brain states using neural field theory. J Neurosci Methods 2015; 258:28-45. [PMID: 26523766 DOI: 10.1016/j.jneumeth.2015.09.026] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2015] [Revised: 09/13/2015] [Accepted: 09/16/2015] [Indexed: 12/01/2022]
Abstract
A real-time fitting system is developed and used to fit the predictions of an established physiologically-based neural field model to electroencephalographic spectra, yielding a trajectory in a physiological parameter space that parametrizes intracortical, intrathalamic, and corticothalamic feedbacks as the arousal state evolves continuously over time. This avoids traditional sleep/wake staging (e.g., using Rechtschaffen-Kales stages), which is fundamentally limited because it forces classification of continuous dynamics into a few discrete categories that are neither physiologically informative nor individualized. The classification is also subject to substantial interobserver disagreement because traditional staging relies in part on subjective evaluations. The fitting routine objectively and robustly tracks arousal parameters over the course of a full night of sleep, and runs in real-time on a desktop computer. The system developed here supersedes discrete staging systems by representing arousal states in terms of physiology, and provides an objective measure of arousal state which solves the problem of interobserver disagreement. Discrete stages from traditional schemes can be expressed in terms of model parameters for backward compatibility with prior studies.
Collapse
Affiliation(s)
- R G Abeysuriya
- School of Physics, University of Sydney, New South Wales 2006, Australia; Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia; Neurosleep, 431 Glebe Point Rd, Glebe, New South Wales 2037, Australia; Brain Dynamics Center, Sydney Medical School - Western, University of Sydney, Westmead, New South Wales 2145, Australia.
| | - P A Robinson
- School of Physics, University of Sydney, New South Wales 2006, Australia; Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia; Neurosleep, 431 Glebe Point Rd, Glebe, New South Wales 2037, Australia; Brain Dynamics Center, Sydney Medical School - Western, University of Sydney, Westmead, New South Wales 2145, Australia
| |
Collapse
|
48
|
Zhao X, Robinson PA. Generalized seizures in a neural field model with bursting dynamics. J Comput Neurosci 2015; 39:197-216. [PMID: 26282528 DOI: 10.1007/s10827-015-0571-7] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2014] [Revised: 07/02/2015] [Accepted: 07/26/2015] [Indexed: 11/27/2022]
Abstract
The mechanisms underlying generalized seizures are explored with neural field theory. A corticothalamic neural field model that has accounted for multiple brain activity phenomena and states is used to explore changes leading to pathological seizure states. It is found that absence seizures arise from instabilities in the system and replicate experimental studies in numerous animal models and clinical studies.
Collapse
Affiliation(s)
- X Zhao
- School of Physics, The University of Sydney, Sydney, New South Wales, 2006, Australia.
- Center for Integrative Brain Function, University of Sydney, NSW, 2006, Australia.
- Neurosleep, 431 Glebe Point Rd, Glebe, New South Wales, 2037, Australia.
- Cooperative Research Center for Alertness, Safety, and Productivity, University of Sydney, NSW, 2006, Australia.
| | - P A Robinson
- School of Physics, The University of Sydney, Sydney, New South Wales, 2006, Australia
- Center for Integrative Brain Function, University of Sydney, NSW, 2006, Australia
- Neurosleep, 431 Glebe Point Rd, Glebe, New South Wales, 2037, Australia
- Cooperative Research Center for Alertness, Safety, and Productivity, University of Sydney, NSW, 2006, Australia
| |
Collapse
|
49
|
Robinson PA, Roy N. Neural field theory of nonlinear wave-wave and wave-neuron processes. Phys Rev E Stat Nonlin Soft Matter Phys 2015; 91:062719. [PMID: 26172747 DOI: 10.1103/physreve.91.062719] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/23/2015] [Indexed: 06/04/2023]
Abstract
Systematic expansion of neural field theory equations in terms of nonlinear response functions is carried out to enable a wide variety of nonlinear wave-wave and wave-neuron processes to be treated systematically in systems involving multiple neural populations. The results are illustrated by analyzing second-harmonic generation, and they can also be applied to wave-wave coalescence, multiharmonic generation, facilitation, depression, refractoriness, and other nonlinear processes.
Collapse
Affiliation(s)
- P A Robinson
- School of Physics, University of Sydney, New South Wales 2006, Australia
- Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
- Neurosleep, 431 Glebe Point Road, Glebe, New South Wales 2037, Australia
| | - N Roy
- School of Physics, University of Sydney, New South Wales 2006, Australia
- Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| |
Collapse
|
50
|
Zhao X, Kim JW, Robinson PA. Slow-wave oscillations in a corticothalamic model of sleep and wake. J Theor Biol 2015; 370:93-102. [PMID: 25659479 DOI: 10.1016/j.jtbi.2015.01.028] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2014] [Revised: 01/21/2015] [Accepted: 01/24/2015] [Indexed: 11/27/2022]
Abstract
A physiologically-based corticothalamic neural field model is used to study slow wave oscillations including cortical UP and DOWN states in deep sleep by extending it to incorporate bursting dynamics of neurons in the thalamic reticular nucleus. The interplay of local bursting dynamics and network interactions produces the cortical UP and DOWN states of slow wave sleep while preserving previously verified model predictions in the wake state. Results show that EEG spectral features in wake and sleep are reproduced. The bursting is subthreshold but acts to intensify the amplitude of oscillations in slow wave sleep with deep UP/DOWN oscillations on the cortex emerging naturally. Furthermore, there is a continuous cycle between the two regimes, rather than a flip-flop between discrete states.
Collapse
Affiliation(s)
- X Zhao
- School of Physics, The University of Sydney, Sydney, New South Wales 2006, Australia; Center of Research Excellence, Neurosleep, 431 Glebe Point Rd, Glebe, New South Wales 2037, Australia; Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia; Cooperative Research Center for Alertness, Safety, and Productivity, University of Sydney, New South Wales 2006, Australia.
| | - J W Kim
- School of Physics, The University of Sydney, Sydney, New South Wales 2006, Australia; Center of Research Excellence, Neurosleep, 431 Glebe Point Rd, Glebe, New South Wales 2037, Australia
| | - P A Robinson
- School of Physics, The University of Sydney, Sydney, New South Wales 2006, Australia; Center of Research Excellence, Neurosleep, 431 Glebe Point Rd, Glebe, New South Wales 2037, Australia; Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia; Cooperative Research Center for Alertness, Safety, and Productivity, University of Sydney, New South Wales 2006, Australia
| |
Collapse
|