1
|
Abramian D, Larsson M, Eklund A, Aganj I, Westin CF, Behjat H. Diffusion-informed spatial smoothing of fMRI data in white matter using spectral graph filters. Neuroimage 2021; 237:118095. [PMID: 34000402 PMCID: PMC8356807 DOI: 10.1016/j.neuroimage.2021.118095] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2020] [Revised: 03/07/2021] [Accepted: 04/13/2021] [Indexed: 12/15/2022] Open
Abstract
Brain activation mapping using functional magnetic resonance imaging (fMRI) has been extensively studied in brain gray matter (GM), whereas in large disregarded for probing white matter (WM). This unbalanced treatment has been in part due to controversies in relation to the nature of the blood oxygenation level-dependent (BOLD) contrast in WM and its detectability. However, an accumulating body of studies has provided solid evidence of the functional significance of the BOLD signal in WM and has revealed that it exhibits anisotropic spatio-temporal correlations and structure-specific fluctuations concomitant with those of the cortical BOLD signal. In this work, we present an anisotropic spatial filtering scheme for smoothing fMRI data in WM that accounts for known spatial constraints on the BOLD signal in WM. In particular, the spatial correlation structure of the BOLD signal in WM is highly anisotropic and closely linked to local axonal structure in terms of shape and orientation, suggesting that isotropic Gaussian filters conventionally used for smoothing fMRI data are inadequate for denoising the BOLD signal in WM. The fundamental element in the proposed method is a graph-based description of WM that encodes the underlying anisotropy observed across WM, derived from diffusion-weighted MRI data. Based on this representation, and leveraging graph signal processing principles, we design subject-specific spatial filters that adapt to a subject's unique WM structure at each position in the WM that they are applied at. We use the proposed filters to spatially smooth fMRI data in WM, as an alternative to the conventional practice of using isotropic Gaussian filters. We test the proposed filtering approach on two sets of simulated phantoms, showcasing its greater sensitivity and specificity for the detection of slender anisotropic activations, compared to that achieved with isotropic Gaussian filters. We also present WM activation mapping results on the Human Connectome Project's 100-unrelated subject dataset, across seven functional tasks, showing that the proposed method enables the detection of streamline-like activations within axonal bundles.
Collapse
Affiliation(s)
- David Abramian
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden; Center for Medical Image Science and Visualization, Linköping University, Linköping, Sweden.
| | - Martin Larsson
- Centre of Mathematical Sciences, Lund University, Lund, Sweden
| | - Anders Eklund
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden; Center for Medical Image Science and Visualization, Linköping University, Linköping, Sweden; Department of Computer and Information Science, Linköping University, Linköping, Sweden
| | - Iman Aganj
- Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Boston, USA; Computer Science and Artificial Intelligence Lab, Massachusetts Institute of Technology, Cambridge, USA
| | - Carl-Fredrik Westin
- Department of Radiology, Brigham and Women's Hospital, Harvard Medical School, Boston, USA
| | - Hamid Behjat
- Department of Biomedical Engineering, Lund University, Lund, Sweden; Department of Radiology, Brigham and Women's Hospital, Harvard Medical School, Boston, USA; Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Boston, USA.
| |
Collapse
|
2
|
Friston KJ, Fagerholm ED, Zarghami TS, Parr T, Hipólito I, Magrou L, Razi A. Parcels and particles: Markov blankets in the brain. Netw Neurosci 2021; 5:211-251. [PMID: 33688613 PMCID: PMC7935044 DOI: 10.1162/netn_a_00175] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2020] [Accepted: 11/24/2020] [Indexed: 11/04/2022] Open
Abstract
At the inception of human brain mapping, two principles of functional anatomy underwrote most conceptions-and analyses-of distributed brain responses: namely, functional segregation and integration. There are currently two main approaches to characterizing functional integration. The first is a mechanistic modeling of connectomics in terms of directed effective connectivity that mediates neuronal message passing and dynamics on neuronal circuits. The second phenomenological approach usually characterizes undirected functional connectivity (i.e., measurable correlations), in terms of intrinsic brain networks, self-organized criticality, dynamical instability, and so on. This paper describes a treatment of effective connectivity that speaks to the emergence of intrinsic brain networks and critical dynamics. It is predicated on the notion of Markov blankets that play a fundamental role in the self-organization of far from equilibrium systems. Using the apparatus of the renormalization group, we show that much of the phenomenology found in network neuroscience is an emergent property of a particular partition of neuronal states, over progressively coarser scales. As such, it offers a way of linking dynamics on directed graphs to the phenomenology of intrinsic brain networks.
Collapse
Affiliation(s)
- Karl J. Friston
- Wellcome Centre for Human Neuroimaging, University College London, London, United Kingdom
| | - Erik D. Fagerholm
- Department of Neuroimaging, King’s College London, London, United Kingdom
| | - Tahereh S. Zarghami
- Bio-Electric Department, School of Electrical and Computer Engineering, University of Tehran, Amirabad, Tehran, Iran
| | - Thomas Parr
- Wellcome Centre for Human Neuroimaging, University College London, London, United Kingdom
| | - Inês Hipólito
- Berlin School of Mind and Brain, and Institut für Philosophie, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Loïc Magrou
- Univ Lyon, Université Claude Bernard Lyon 1, Inserm, Stem Cell and Brain Research Institute U1208, Bron, France
| | - Adeel Razi
- Wellcome Centre for Human Neuroimaging, University College London, London, United Kingdom
- Turner Institute for Brain and Mental Health, Monash University, Clayton, Australia
| |
Collapse
|
3
|
Miller RL, Erhardt EB, Agcaoglu O, Allen EA, Michael AM, Turner JA, Bustillo J, Ford JM, Mathalon DH, Van Erp TGM, Potkin S, Preda A, Pearlson G, Calhoun VD. Multidimensional frequency domain analysis of full-volume fMRI reveals significant effects of age, gender, and mental illness on the spatiotemporal organization of resting-state brain activity. Front Neurosci 2015; 9:203. [PMID: 26136646 PMCID: PMC4468831 DOI: 10.3389/fnins.2015.00203] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2014] [Accepted: 05/21/2015] [Indexed: 12/30/2022] Open
Abstract
Clinical research employing functional magnetic resonance imaging (fMRI) is often conducted within the connectionist paradigm, focusing on patterns of connectivity between voxels, regions of interest (ROIs) or spatially distributed functional networks. Connectivity-based analyses are concerned with pairwise correlations of the temporal activation associated with restrictions of the whole-brain hemodynamic signal to locations of a priori interest. There is a more abstract question however that such spatially granular correlation-based approaches do not elucidate: Are the broad spatiotemporal organizing principles of brains in certain populations distinguishable from those of others? Global patterns (in space and time) of hemodynamic activation are rarely scrutinized for features that might characterize complex psychiatric conditions, aging effects or gender—among other variables of potential interest to researchers. We introduce a canonical, transparent technique for characterizing the role in overall brain activation of spatially scaled periodic patterns with given temporal recurrence rates. A core feature of our technique is the spatiotemporal spectral profile (STSP), a readily interpretable 2D reduction of the native four-dimensional brain × time frequency domain that is still “big enough” to capture important group differences in globally patterned brain activation. Its power to distinguish populations of interest is demonstrated on a large balanced multi-site resting fMRI dataset with nearly equal numbers of schizophrenia patients and healthy controls. Our analysis reveals striking differences in the spatiotemporal organization of brain activity that correlate with the presence of diagnosed schizophrenia, as well as with gender and age. To the best of our knowledge, this is the first demonstration that a 4D frequency domain analysis of full volume fMRI data exposes clinically or demographically relevant differences in resting-state brain function.
Collapse
Affiliation(s)
| | - Erik B Erhardt
- Department of Mathematics and Statistics, University of New Mexico Albuquerque, NM, USA
| | - Oktay Agcaoglu
- The Mind Research Network Albuquerque, NM, USA ; Department of Electrical and Computer Engineering, University of New Mexico Albuquerque, NM, USA
| | - Elena A Allen
- Bergen Center for Neuropsychiatric Research Bergen, Norway
| | - Andrew M Michael
- Geisinger Autism and Developmental Medicine Institute Lewisburg, PA, USA
| | - Jessica A Turner
- Department of Psychology and Neuroscience, Georgia State University Atlanta, GA, USA
| | - Juan Bustillo
- Department of Psychiatry, University of New Mexico School of Medicine Albuquerque, NM, USA
| | - Judith M Ford
- Department of Psychiatry, University of California San Francisco School of Medicine San Francisco, CA, USA
| | - Daniel H Mathalon
- Department of Psychiatry, University of California San Francisco School of Medicine San Francisco, CA, USA
| | - Theo G M Van Erp
- Department of Psychiatry and Human Behavior, School of Medicine, University of California Irvine Irvine, CA, USA
| | - Steven Potkin
- Department of Psychiatry and Human Behavior, School of Medicine, University of California Irvine Irvine, CA, USA
| | - Adrian Preda
- Department of Psychiatry and Human Behavior, School of Medicine, University of California Irvine Irvine, CA, USA
| | - Godfrey Pearlson
- Department of Psychiatry, Yale University School of Medicine New Haven, CT, USA
| | - Vince D Calhoun
- The Mind Research Network Albuquerque, NM, USA ; Department of Psychiatry, University of New Mexico School of Medicine Albuquerque, NM, USA ; Department of Psychiatry, Yale University School of Medicine New Haven, CT, USA ; Department of Electrical and Computer Engineering, University of New Mexico Albuquerque, NM, USA
| |
Collapse
|
4
|
Zalesky A, Breakspear M. Towards a statistical test for functional connectivity dynamics. Neuroimage 2015; 114:466-70. [PMID: 25818688 DOI: 10.1016/j.neuroimage.2015.03.047] [Citation(s) in RCA: 211] [Impact Index Per Article: 23.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2014] [Revised: 02/09/2015] [Accepted: 03/18/2015] [Indexed: 12/15/2022] Open
Abstract
Sliding-window correlation is an emerging method for mapping time-resolved, resting-state functional connectivity. To avoid mapping spurious connectivity fluctuations (false positives), Leonardi and Van De Ville recently recommended choosing a window length exceeding the longest wavelength composing the BOLD signal, usually assumed to be ~100s. Here, we provide further statistical support for this rule of thumb. However, we demonstrate that non-stationary fluctuations in functional connectivity can in theory be detected with much shorter window lengths (e.g. 40s), while maintaining nominal control of false positives. We find that statistical power is near-maximal for window lengths chosen according to Leonardi and Van De Ville's rule of thumb. Furthermore, we lay some foundations for a parametric test to identify non-stationary fluctuations in functional connectivity, also noting limitations of the sinusoidal model upon which our work, and the work of Leonardi and Van De Ville, is based. Most notably, our analytical results pertain to covariances, as does our statistical test, whereas functional connectivity is more commonly measured using correlations.
Collapse
Affiliation(s)
- Andrew Zalesky
- Melbourne Neuropsychiatry Centre and Melbourne Health, The University of Melbourne, Victoria 3010, Australia; Melbourne School of Engineering, The University of Melbourne, Victoria 3010, Australia.
| | - Michael Breakspear
- QIMR Berghofer Medical Research Institute, Brisbane, Queensland 4029, Australia; The Royal Brisbane and Womens Hospital, Brisbane, Queensland 4029, Australia.
| |
Collapse
|
5
|
Spatiotemporal hemodynamic response functions derived from physiology. J Theor Biol 2014; 347:118-36. [PMID: 24398024 DOI: 10.1016/j.jtbi.2013.12.027] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2012] [Revised: 11/28/2013] [Accepted: 12/27/2013] [Indexed: 12/12/2022]
Abstract
Probing neural activity with functional magnetic resonance imaging (fMRI) relies upon understanding the hemodynamic response to changes in neural activity. Although existing studies have extensively characterized the temporal hemodynamic response, less is understood about the spatial and spatiotemporal hemodynamic responses. This study systematically characterizes the spatiotemporal response by deriving the hemodynamic response due to a short localized neural drive, i.e., the spatiotemporal hemodynamic response function (stHRF) from a physiological model of hemodynamics based on a poroelastic model of cortical tissue. In this study, the model's boundary conditions are clarified and a resulting nonlinear hemodynamic wave equation is derived. From this wave equation, damped linear hemodynamic waves are predicted from the stHRF. The main features of these waves depend on two physiological parameters: wave propagation speed, which depends on mean cortical stiffness, and damping which depends on effective viscosity. Some of these predictions were applied and validated in a companion study (Aquino et al., 2012). The advantages of having such a theory for the stHRF include improving the interpretation of spatiotemporal dynamics in fMRI data; improving estimates of neural activity with fMRI spatiotemporal deconvolution; and enabling wave interactions between hemodynamic waves to be predicted and exploited to improve the signal to noise ratio of fMRI.
Collapse
|
6
|
Aquino KM, Schira MM, Robinson PA, Drysdale PM, Breakspear M. Hemodynamic traveling waves in human visual cortex. PLoS Comput Biol 2012; 8:e1002435. [PMID: 22457612 PMCID: PMC3310706 DOI: 10.1371/journal.pcbi.1002435] [Citation(s) in RCA: 68] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2011] [Accepted: 02/06/2012] [Indexed: 11/18/2022] Open
Abstract
Functional MRI (fMRI) experiments rely on precise characterization of the blood oxygen level dependent (BOLD) signal. As the spatial resolution of fMRI reaches the sub-millimeter range, the need for quantitative modelling of spatiotemporal properties of this hemodynamic signal has become pressing. Here, we find that a detailed physiologically-based model of spatiotemporal BOLD responses predicts traveling waves with velocities and spatial ranges in empirically observable ranges. Two measurable parameters, related to physiology, characterize these waves: wave velocity and damping rate. To test these predictions, high-resolution fMRI data are acquired from subjects viewing discrete visual stimuli. Predictions and experiment show strong agreement, in particular confirming BOLD waves propagating for at least 5–10 mm across the cortical surface at speeds of 2–12 mm s-1. These observations enable fundamentally new approaches to fMRI analysis, crucial for fMRI data acquired at high spatial resolution. Functional magnetic resonance imaging (fMRI) experiments have advanced our understanding of the structure and function of the human brain. Dynamic changes in the flow and concentration of oxygen in blood are observed experimentally in fMRI data via the blood oxygen level dependent (BOLD) signal. Since neuronal activity induces this hemodynamic response, the BOLD signal provides a noninvasive measure of neuronal activity. Understanding the mechanisms that drive this BOLD response is fundamental for accurately inferring the underlying neuronal activity. The goal of this study is to systematically predict spatiotemporal hemodynamics from a biophysical model, then test these in a high resolution fMRI study of the visual cortex. Using this theory, we predict and empirically confirm the existence of hemodynamic waves in cortex – a striking and novel finding.
Collapse
Affiliation(s)
- Kevin M Aquino
- School of Physics, University of Sydney, New South Wales, Australia.
| | | | | | | | | |
Collapse
|
7
|
Exploiting the potential of three dimensional spatial wavelet analysis to explore nesting of temporal oscillations and spatial variance in simultaneous EEG-fMRI data. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2010; 105:67-79. [PMID: 21094179 DOI: 10.1016/j.pbiomolbio.2010.11.003] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/13/2010] [Revised: 11/10/2010] [Accepted: 11/11/2010] [Indexed: 11/20/2022]
Abstract
Synchronization of the activity in neural networks is a fundamental mechanism of brain function, putatively serving the integration of computations on multiple spatial and temporal scales. Time scales are thought to be nested within distinct spatial scales, so that whereas fast oscillations may integrate local networks, slow oscillations might integrate computations across distributed brain areas. We here describe a newly developed approach that provides potential for the further substantiation of this hypothesis in future studies. We demonstrate the feasibility and important caveats of a novel wavelet-based means of relating time series of three-dimensional spatial variance (energy) of fMRI data to time series of temporal variance of EEG. The spatial variance of fMRI data was determined by employing the three-dimensional dual-tree complex wavelet transform. The temporal variance of EEG data was estimated by using traditional continuous complex wavelets. We tested our algorithm on artificial signals with known signal-to-noise ratios and on empirical resting state EEG-fMRI data obtained from four healthy human subjects. By employing the human posterior alpha rhythm as an exemplar, we demonstrated face validity of the approach. We believe that the proposed method can serve as a suitable tool for future research on the spatiotemporal properties of brain dynamics, hence moving beyond analyses based exclusively in one domain or the other.
Collapse
|
8
|
Robinson PA, Henderson JA, Matar E, Riley P, Gray RT. Dynamical reconnection and stability constraints on cortical network architecture. PHYSICAL REVIEW LETTERS 2009; 103:108104. [PMID: 19792345 DOI: 10.1103/physrevlett.103.108104] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/05/2009] [Indexed: 05/26/2023]
Abstract
Stability under dynamical changes to network connectivity is invoked alongside previous criteria to constrain brain network architecture. A new hierarchical network is introduced that satisfies all these constraints, unlike more commonly studied regular, random, and small-world networks. It is shown that hierarchical networks can simultaneously have high clustering, short path lengths, and low wiring costs, while being robustly stable under large scale reconnection of substructures.
Collapse
Affiliation(s)
- P A Robinson
- School of Physics, University of Sydney, Sydney, New South Wales 2006, Australia
| | | | | | | | | |
Collapse
|
9
|
Stephan KE, Riera JJ, Deco G, Horwitz B. The Brain Connectivity Workshops: moving the frontiers of computational systems neuroscience. Neuroimage 2008; 42:1-9. [PMID: 18511300 DOI: 10.1016/j.neuroimage.2008.04.167] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2008] [Revised: 04/03/2008] [Accepted: 04/11/2008] [Indexed: 11/30/2022] Open
Abstract
Understanding the link between neurobiology and cognition requires that neuroscience moves beyond mere structure-function correlations. An explicit systems perspective is needed in which putative mechanisms of how brain function is constrained by brain structure are mathematically formalized and made accessible for experimental investigation. Such a systems approach critically rests on a better understanding of brain connectivity in its various forms. Since 2002, frontier topics of connectivity and neural system analysis have been discussed in a multidisciplinary annual meeting, the Brain Connectivity Workshop (BCW), bringing together experimentalists and theorists from various fields. This article summarizes some of the main discussions at the two most recent workshops, 2006 at Sendai, Japan, and 2007 at Barcelona, Spain: (i) investigation of cortical micro- and macrocircuits, (ii) models of neural dynamics at multiple scales, (iii) analysis of "resting state" networks, and (iv) linking anatomical to functional connectivity. Finally, we outline some central challenges and research trajectories in computational systems neuroscience for the next years.
Collapse
Affiliation(s)
- Klaas Enno Stephan
- Wellcome Trust Centre for Neuroimaging, University College London, 12 Queen Square, London WC1N3BG, UK.
| | | | | | | |
Collapse
|
10
|
Bathellier B, Van De Ville D, Blu T, Unser M, Carleton A. Wavelet-based multi-resolution statistics for optical imaging signals: Application to automated detection of odour activated glomeruli in the mouse olfactory bulb. Neuroimage 2006; 34:1020-35. [PMID: 17185002 DOI: 10.1016/j.neuroimage.2006.10.038] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2006] [Revised: 10/04/2006] [Accepted: 10/26/2006] [Indexed: 11/29/2022] Open
Abstract
Optical imaging techniques offer powerful solutions to capture brain networks processing in animals, especially when activity is distributed in functionally distinct spatial domains. Despite the progress in imaging techniques, the standard analysis procedures and statistical assessments for this type of data are still limited. In this paper, we perform two in vivo non-invasive optical recording techniques in the mouse olfactory bulb, using a genetically expressed activity reporter fluorescent protein (synaptopHfluorin) and intrinsic signals of the brain. For both imaging techniques, we show that the odour-triggered signals can be accurately parameterized using linear models. Fitting the models allows us to extract odour specific signals with a reduced level of noise compared to standard methods. In addition, the models serve to evaluate statistical significance, using a wavelet-based framework that exploits spatial correlation at different scales. We propose an extension of this framework to extract activation patterns at specific wavelet scales. This method is especially interesting to detect the odour inputs that segregate on the olfactory bulb in small spherical structures called glomeruli. Interestingly, with proper selection of wavelet scales, we can isolate significantly activated glomeruli and thus determine the odour map in an automated manner. Comparison against manual detection of glomeruli shows the high accuracy of the proposed method. Therefore, beyond the advantageous alternative to the existing treatments of optical imaging signals in general, our framework propose an interesting procedure to dissect brain activation patterns on multiple scales with statistical control.
Collapse
Affiliation(s)
- Brice Bathellier
- Flavour Perception Group, Brain Mind Institute, Ecole Polytechnique Fédérale de Lausanne, (EPFL), CH-1015, Switzerland
| | | | | | | | | |
Collapse
|