ABSTRACT
Temporal lobe epilepsy (TLE) is the most common pharmaco-resistant epilepsy in adults. While primarily associated with mesiotemporal pathology, recent evidence suggests that brain alterations in TLE extend beyond the paralimbic epicenter and impact macroscale function and cognitive functions, particularly memory. Using connectome-wide manifold learning and generative models of effective connectivity, we examined functional topography and directional signal flow patterns between large-scale neural circuits in TLE at rest. Studying a multisite cohort of 95 patients with TLE and 95 healthy controls, we observed atypical functional topographies in the former group, characterized by reduced differentiation between sensory and transmodal association cortices, with most marked effects in bilateral temporo-limbic and ventromedial prefrontal cortices. These findings were consistent across all study sites, present in left and right lateralized patients, and validated in a subgroup of patients with histopathological validation of mesiotemporal sclerosis and post-surgical seizure freedom. Moreover, they were replicated in an independent cohort of 30 TLE patients and 40 healthy controls. Further analyses demonstrated that reduced differentiation related to decreased functional signal flow into and out of temporolimbic cortical systems and other brain networks. Parallel analyses of structural and diffusion-weighted MRI data revealed that topographic alterations were independent of TLE-related cortical thinning but partially mediated by white matter microstructural changes that radiated away from paralimbic circuits. Finally, we found a strong association between the degree of functional alterations and behavioral markers of memory dysfunction. Our work illustrates the complex landscape of macroscale functional imbalances in TLE, which can serve as intermediate markers bridging microstructural changes and cognitive impairment.
Subject(s)
Connectome , Epilepsy, Temporal Lobe , Humans , Epilepsy, Temporal Lobe/physiopathology , Epilepsy, Temporal Lobe/diagnostic imaging , Epilepsy, Temporal Lobe/pathology , Female , Male , Adult , Middle Aged , Magnetic Resonance Imaging , Young Adult , Brain/diagnostic imaging , Brain/physiopathology , Brain/pathology , Cohort Studies , Nerve Net/diagnostic imaging , Nerve Net/physiopathology , Nerve Net/pathologyABSTRACT
Temporal lobe epilepsy (TLE) is one of the most common pharmaco-resistant epilepsies in adults. While hippocampal pathology is the hallmark of this condition, emerging evidence indicates that brain alterations extend beyond the mesiotemporal epicenter and affect macroscale brain function and cognition. We studied macroscale functional reorganization in TLE, explored structural substrates, and examined cognitive associations. We investigated a multisite cohort of 95 patients with pharmaco-resistant TLE and 95 healthy controls using state-of-the-art multimodal 3T magnetic resonance imaging (MRI). We quantified macroscale functional topographic organization using connectome dimensionality reduction techniques and estimated directional functional flow using generative models of effective connectivity. We observed atypical functional topographies in patients with TLE relative to controls, manifesting as reduced functional differentiation between sensory/motor networks and transmodal systems such as the default mode network, with peak alterations in bilateral temporal and ventromedial prefrontal cortices. TLE-related topographic changes were consistent in all three included sites and reflected reductions in hierarchical flow patterns between cortical systems. Integration of parallel multimodal MRI data indicated that these findings were independent of TLE-related cortical grey matter atrophy, but mediated by microstructural alterations in the superficial white matter immediately beneath the cortex. The magnitude of functional perturbations was robustly associated with behavioral markers of memory function. Overall, this work provides converging evidence for macroscale functional imbalances, contributing microstructural alterations, and their associations with cognitive dysfunction in TLE.
ABSTRACT
After a first episode of major depressive disorder (MDD), there is substantial risk for a long-term remitting-relapsing course. Prevention and early interventions are thus critically important. Various studies have examined the feasibility of detecting at-risk individuals based on out-of-sample predictions about the future occurrence of depression. However, functional magnetic resonance imaging (fMRI) has received very little attention for this purpose so far. Here, we explored the utility of generative models (i.e. different dynamic causal models, DCMs) as well as functional connectivity (FC) for predicting future episodes of depression in never-depressed adults, using a large dataset (N = 906) of task-free ("resting state") fMRI data from the UK Biobank (UKB). Connectivity analyses were conducted using timeseries from pre-computed spatially independent components of different dimensionalities. Over a three-year period, 50% of selected participants showed indications of at least one depressive episode, while the other 50% did not. Using nested cross-validation for training and a held-out test set (80/20 split), we systematically examined the combination of 8 connectivity feature sets and 17 classifiers. We found that a generative embedding procedure based on combining regression DCM (rDCM) with a support vector machine (SVM) enabled the best predictions, both on the training set (0.63 accuracy, 0.66 area under the curve, AUC) and the test set (0.62 accuracy, 0.64 AUC; p < 0.001). However, on the test set, rDCM was only slightly superior to predictions based on FC (0.59 accuracy, 0.61 AUC). Interpreting model predictions based on SHAP (SHapley Additive exPlanations) values suggested that the most predictive connections were widely distributed and not confined to specific networks. Overall, our analyses suggest (i) ways of improving future fMRI-based generative embedding approaches for the early detection of individuals at-risk for depression and that (ii) achieving accuracies of clinical utility may require combination of fMRI with other data modalities.
Subject(s)
Brain , Depressive Disorder, Major , Adult , Humans , Magnetic Resonance Imaging/methods , Support Vector Machine , Models, NeurologicalABSTRACT
Regression dynamic causal modeling (rDCM) is a novel and computationally highly efficient method for inferring effective connectivity at the whole-brain level. While face and construct validity of rDCM have already been demonstrated, here we assessed its test-retest reliability-a test-theoretical property of particular importance for clinical applications-together with group-level consistency of connection-specific estimates and consistency of whole-brain connectivity patterns over sessions. Using the Human Connectome Project dataset for eight different paradigms (tasks and rest) and two different parcellation schemes, we found that rDCM provided highly consistent connectivity estimates at the group level across sessions. Second, while test-retest reliability was limited when averaging over all connections (range of mean intraclass correlation coefficient 0.24-0.42 over tasks), reliability increased with connection strength, with stronger connections showing good to excellent test-retest reliability. Third, whole-brain connectivity patterns by rDCM allowed for identifying individual participants with high (and in some cases perfect) accuracy. Comparing the test-retest reliability of rDCM connectivity estimates with measures of functional connectivity, rDCM performed favorably-particularly when focusing on strong connections. Generally, for all methods and metrics, task-based connectivity estimates showed greater reliability than those from the resting state. Our results underscore the potential of rDCM for human connectomics and clinical applications.
ABSTRACT
In generative modeling of neuroimaging data, such as dynamic causal modeling (DCM), one typically considers several alternative models, either to determine the most plausible explanation for observed data (Bayesian model selection) or to account for model uncertainty (Bayesian model averaging). Both procedures rest on estimates of the model evidence, a principled trade-off between model accuracy and complexity. In the context of DCM, the log evidence is usually approximated using variational Bayes. Although this approach is highly efficient, it makes distributional assumptions and is vulnerable to local extrema. This paper introduces the use of thermodynamic integration (TI) for Bayesian model selection and averaging in the context of DCM. TI is based on Markov chain Monte Carlo sampling which is asymptotically exact but orders of magnitude slower than variational Bayes. In this paper, we explain the theoretical foundations of TI, covering key concepts such as the free energy and its origins in statistical physics. Our aim is to convey an in-depth understanding of the method starting from its historical origin in statistical physics. In addition, we demonstrate the practical application of TI via a series of examples which serve to guide the user in applying this method. Furthermore, these examples demonstrate that, given an efficient implementation and hardware capable of parallel processing, the challenge of high computational demand can be overcome successfully. The TI implementation presented in this paper is freely available as part of the open source software TAPAS. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s11571-021-09696-9.
ABSTRACT
The basal ganglia (BG) are a group of subcortical nuclei responsible for motor and executive function. Central to BG function are striatal cells expressing D1 (D1R) and D2 (D2R) dopamine receptors. D1R and D2R cells are considered functional antagonists that facilitate voluntary movements and inhibit competing motor patterns, respectively. However, whether they maintain a uniform function across the striatum and what influence they exert outside the BG is unclear. Here, we address these questions by combining optogenetic activation of D1R and D2R cells in the mouse ventrolateral caudoputamen with fMRI. Striatal D1R/D2R stimulation evokes distinct activity within the BG-thalamocortical network and differentially engages cerebellar and prefrontal regions. Computational modeling of effective connectivity confirms that changes in D1R/D2R output drive functional relationships between these regions. Our results suggest a complex functional organization of striatal D1R/D2R cells and hint toward an interconnected fronto-BG-cerebellar network modulated by striatal D1R and D2R cells.
Subject(s)
Basal Ganglia/metabolism , Corpus Striatum/metabolism , Neostriatum/metabolism , Neurons/metabolism , Optogenetics , Receptors, Dopamine D1/metabolism , Receptors, Dopamine D2/metabolism , Animals , Female , Male , MiceABSTRACT
Dynamic Causal Modeling (DCM) is a Bayesian framework for inferring on hidden (latent) neuronal states, based on measurements of brain activity. Since its introduction in 2003 for functional magnetic resonance imaging data, DCM has been extended to electrophysiological data, and several variants have been developed. Their biophysically motivated formulations make these models promising candidates for providing a mechanistic understanding of human brain dynamics, both in health and disease. However, due to their complexity and reliance on concepts from several fields, fully understanding the mathematical and conceptual basis behind certain variants of DCM can be challenging. At the same time, a solid theoretical knowledge of the models is crucial to avoid pitfalls in the application of these models and interpretation of their results. In this paper, we focus on one of the most advanced formulations of DCM, i.e. conductance-based DCM for cross-spectral densities, whose components are described across multiple technical papers. The aim of the present article is to provide an accessible exposition of the mathematical background, together with an illustration of the model's behavior. To this end, we include step-by-step derivations of the model equations, point to important aspects in the software implementation of those models, and use simulations to provide an intuitive understanding of the type of responses that can be generated and the role that specific parameters play in the model. Furthermore, all code utilized for our simulations is made publicly available alongside the manuscript to allow readers an easy hands-on experience with conductance-based DCM.
Subject(s)
Brain Mapping/methods , Brain/physiology , Computer Simulation , Bayes Theorem , Electrophysiological Phenomena , Humans , Magnetic Resonance Imaging/methods , Models, Neurological , Nerve Net/physiology , Neurons , SoftwareABSTRACT
Interoception, the perception of internal bodily states, is thought to be inextricably linked to affective qualities such as anxiety. Although interoception spans sensory to metacognitive processing, it is not clear whether anxiety is differentially related to these processing levels. Here we investigated this question in the domain of breathing, using computational modeling and high-field (7 T) fMRI to assess brain activity relating to dynamic changes in inspiratory resistance of varying predictability. Notably, the anterior insula was associated with both breathing-related prediction certainty and prediction errors, suggesting an important role in representing and updating models of the body. Individuals with low versus moderate anxiety traits showed differential anterior insula activity for prediction certainty. Multi-modal analyses of data from fMRI, computational assessments of breathing-related metacognition, and questionnaires demonstrated that anxiety-interoception links span all levels from perceptual sensitivity to metacognition, with strong effects seen at higher levels of interoceptive processes.
Subject(s)
Interoception , Anxiety , Anxiety Disorders , Heart Rate , Humans , RespirationABSTRACT
Dynamic causal models (DCMs) of electrophysiological data allow, in principle, for inference on hidden, bulk synaptic function in neural circuits. The directed influences between the neuronal elements of modeled circuits are subject to delays due to the finite transmission speed of axonal connections. Ordinary differential equations are therefore not adequate to capture the ensuing circuit dynamics, and delay differential equations (DDEs) are required instead. Previous work has illustrated that the integration of DDEs in DCMs benefits from sophisticated integration schemes in order to ensure rigorous parameter estimation and correct model identification. However, integration schemes that have been proposed for DCMs either emphasize speed (at the possible expense of accuracy) or robustness (but with computational costs that are problematic in practice). In this technical note, we propose an alternative integration scheme that overcomes these shortcomings and offers high computational efficiency while correctly preserving the nature of delayed effects. This integration scheme is available as open-source code in the Translational Algorithms for Psychiatry-Advancing Science (TAPAS) toolbox and can be easily integrated into existing software (SPM) for the analysis of DCMs for electrophysiological data. While this paper focuses on its application to the convolution-based formalism of DCMs, the new integration scheme can be equally applied to more advanced formulations of DCMs (e.g. conductance based models). Our method provides a new option for electrophysiological DCMs that offers the speed required for scientific projects, but also the accuracy required for rigorous translational applications, e.g. in computational psychiatry.
Subject(s)
Brain Mapping/methods , Electrophysiological Phenomena/physiology , Models, Statistical , Algorithms , Brain/physiology , Computer Simulation , Humans , Magnetic Resonance Imaging/methods , Models, Neurological , SoftwareABSTRACT
Psychiatry faces fundamental challenges with regard to mechanistically guided differential diagnosis, as well as prediction of clinical trajectories and treatment response of individual patients. This has motivated the genesis of two closely intertwined fields: (i) Translational Neuromodeling (TN), which develops "computational assays" for inferring patient-specific disease processes from neuroimaging, electrophysiological, and behavioral data; and (ii) Computational Psychiatry (CP), with the goal of incorporating computational assays into clinical decision making in everyday practice. In order to serve as objective and reliable tools for clinical routine, computational assays require end-to-end pipelines from raw data (input) to clinically useful information (output). While these are yet to be established in clinical practice, individual components of this general end-to-end pipeline are being developed and made openly available for community use. In this paper, we present the Translational Algorithms for Psychiatry-Advancing Science (TAPAS) software package, an open-source collection of building blocks for computational assays in psychiatry. Collectively, the tools in TAPAS presently cover several important aspects of the desired end-to-end pipeline, including: (i) tailored experimental designs and optimization of measurement strategy prior to data acquisition, (ii) quality control during data acquisition, and (iii) artifact correction, statistical inference, and clinical application after data acquisition. Here, we review the different tools within TAPAS and illustrate how these may help provide a deeper understanding of neural and cognitive mechanisms of disease, with the ultimate goal of establishing automatized pipelines for predictions about individual patients. We hope that the openly available tools in TAPAS will contribute to the further development of TN/CP and facilitate the translation of advances in computational neuroscience into clinically relevant computational assays.
ABSTRACT
Drugs affecting neuromodulation, for example by dopamine or acetylcholine, take centre stage among therapeutic strategies in psychiatry. These neuromodulators can change both neuronal gain and synaptic plasticity and therefore affect electrophysiological measures. An important goal for clinical diagnostics is to exploit this effect in the reverse direction, i.e., to infer the status of specific neuromodulatory systems from electrophysiological measures. In this study, we provide proof-of-concept that the functional status of cholinergic (specifically muscarinic) receptors can be inferred from electrophysiological data using generative (dynamic causal) models. To this end, we used epidural EEG recordings over two auditory cortical regions during a mismatch negativity (MMN) paradigm in rats. All animals were treated, across sessions, with muscarinic receptor agonists and antagonists at different doses. Together with a placebo condition, this resulted in five levels of muscarinic receptor status. Using a dynamic causal model - embodying a small network of coupled cortical microcircuits - we estimated synaptic parameters and their change across pharmacological conditions. The ensuing parameter estimates associated with (the neuromodulation of) synaptic efficacy showed both graded muscarinic effects and predictive validity between agonistic and antagonistic pharmacological conditions. This finding illustrates the potential utility of generative models of electrophysiological data as computational assays of muscarinic function. In application to EEG data of patients from heterogeneous spectrum diseases, e.g. schizophrenia, such models might help identify subgroups of patients that respond differentially to cholinergic treatments. SIGNIFICANCE STATEMENT: In psychiatry, the vast majority of pharmacological treatments affect actions of neuromodulatory transmitters, e.g. dopamine or acetylcholine. As treatment is largely trial-and-error based, one of the goals for computational psychiatry is to construct mathematical models that can serve as "computational assays" and infer the status of specific neuromodulatory systems in individual patients. This translational neuromodeling strategy has great promise for electrophysiological data in particular but requires careful validation. The present study demonstrates that the functional status of cholinergic (muscarinic) receptors can be inferred from electrophysiological data using dynamic causal models of neural circuits. While accuracy needs to be enhanced and our results must be replicated in larger samples, our current results provide proof-of-concept for computational assays of muscarinic function using EEG.
Subject(s)
Auditory Cortex/physiology , Auditory Perception/physiology , Electrocorticography/methods , Evoked Potentials, Auditory/physiology , Muscarinic Agonists/pharmacology , Muscarinic Antagonists/pharmacology , Receptors, Muscarinic/physiology , Animals , Auditory Cortex/drug effects , Auditory Perception/drug effects , Behavior, Animal/physiology , Electrocorticography/drug effects , Evoked Potentials, Auditory/drug effects , Muscarinic Agonists/administration & dosage , Muscarinic Antagonists/administration & dosage , Pilocarpine/pharmacology , Proof of Concept Study , Rats , Scopolamine/pharmacology , Support Vector MachineABSTRACT
"Resting-state" functional magnetic resonance imaging (rs-fMRI) is widely used to study brain connectivity. So far, researchers have been restricted to measures of functional connectivity that are computationally efficient but undirected, or to effective connectivity estimates that are directed but limited to small networks. Here, we show that a method recently developed for task-fMRI-regression dynamic causal modeling (rDCM)-extends to rs-fMRI and offers both directional estimates and scalability to whole-brain networks. First, simulations demonstrate that rDCM faithfully recovers parameter values over a wide range of signal-to-noise ratios and repetition times. Second, we test construct validity of rDCM in relation to an established model of effective connectivity, spectral DCM. Using rs-fMRI data from nearly 200 healthy participants, rDCM produces biologically plausible results consistent with estimates by spectral DCM. Importantly, rDCM is computationally highly efficient, reconstructing whole-brain networks (>200 areas) within minutes on standard hardware. This opens promising new avenues for connectomics.
Subject(s)
Brain/physiology , Connectome/methods , Magnetic Resonance Imaging/methods , Nerve Net/physiology , Adolescent , Adult , Brain/diagnostic imaging , Connectome/standards , Humans , Magnetic Resonance Imaging/standards , Middle Aged , Models, Theoretical , Nerve Net/diagnostic imaging , Regression Analysis , Young AdultABSTRACT
Connectomics is essential for understanding large-scale brain networks but requires that individual connection estimates are neurobiologically interpretable. In particular, a principle of brain organization is that reciprocal connections between cortical areas are functionally asymmetric. This is a challenge for fMRI-based connectomics in humans where only undirected functional connectivity estimates are routinely available. By contrast, whole-brain estimates of effective (directed) connectivity are computationally challenging, and emerging methods require empirical validation. Here, using a motor task at 7T, we demonstrate that a novel generative model can infer known connectivity features in a whole-brain network (>200 regions, >40,000 connections) highly efficiently. Furthermore, graph-theoretical analyses of directed connectivity estimates identify functional roles of motor areas more accurately than undirected functional connectivity estimates. These results, which can be achieved in an entirely unsupervised manner, demonstrate the feasibility of inferring directed connections in whole-brain networks and open new avenues for human connectomics.
Subject(s)
Brain/diagnostic imaging , Connectome/methods , Motor Cortex/diagnostic imaging , Adult , Aged , Brain/physiology , Female , Functional Neuroimaging/methods , Humans , Magnetic Resonance Imaging , Male , Middle Aged , Models, Neurological , Models, Statistical , Motor Cortex/physiology , Neural Pathways/diagnostic imaging , Neural Pathways/physiology , Regression AnalysisABSTRACT
Aspirin is considered a potential confound for functional magnetic resonance imaging (fMRI) studies. This is because aspirin affects the synthesis of prostaglandin, a vasoactive mediator centrally involved in neurovascular coupling, a process underlying blood oxygenated level dependent (BOLD) responses. Aspirin-induced changes in BOLD signal are a potential confound for fMRI studies of at-risk individuals or patients (e.g. with cardiovascular conditions or stroke) who receive low-dose aspirin prophylactically and are compared to healthy controls without aspirin. To examine the severity of this potential confound, we combined high field (7 Tesla) MRI during a simple hand movement task with a biophysically informed hemodynamic model. We compared elderly individuals receiving aspirin for primary or secondary prophylactic purposes versus age-matched volunteers without aspirin medication, testing for putative differences in BOLD responses. Specifically, we fitted hemodynamic models to BOLD responses from 14 regions activated by the task and examined whether model parameter estimates were significantly altered by aspirin. While our analyses indicate that hemodynamics differed across regions, consistent with the known regional variability of BOLD responses, we neither found a significant main effect of aspirin (i.e., an average effect across brain regions) nor an expected drug × region interaction. While our sample size is not sufficiently large to rule out small-to-medium global effects of aspirin, we had adequate statistical power for detecting the expected interaction. Altogether, our analysis suggests that patients with cardiovascular risk receiving low-dose aspirin for primary or secondary prophylactic purposes do not show strongly altered BOLD signals when compared to healthy controls without aspirin.
Subject(s)
Aspirin , Cardiovascular Diseases , Aged , Brain/diagnostic imaging , Brain Mapping , Heart Disease Risk Factors , Hemodynamics , Humans , Magnetic Resonance Imaging , Oxygen , Risk FactorsABSTRACT
The mesiotemporal lobe (MTL) is implicated in many cognitive processes, is compromised in numerous brain disorders, and exhibits a gradual cytoarchitectural transition from six-layered parahippocampal isocortex to three-layered hippocampal allocortex. Leveraging an ultra-high-resolution histological reconstruction of a human brain, our study showed that the dominant axis of MTL cytoarchitectural differentiation follows the iso-to-allocortical transition and depth-specific variations in neuronal density. Projecting the histology-derived MTL model to in-vivo functional MRI, we furthermore determined how its cytoarchitecture underpins its intrinsic effective connectivity and association to large-scale networks. Here, the cytoarchitectural gradient was found to underpin intrinsic effective connectivity of the MTL, but patterns differed along the anterior-posterior axis. Moreover, while the iso-to-allocortical gradient parametrically represented the multiple-demand relative to task-negative networks, anterior-posterior gradients represented transmodal versus unimodal networks. Our findings establish that the combination of micro- and macrostructural features allow the MTL to represent dominant motifs of whole-brain functional organisation.
Subject(s)
Cognition/physiology , Hippocampus/physiology , Models, Biological , Parahippocampal Gyrus/physiology , Temporal Lobe/physiology , Adult , Aged , Brain Mapping , Female , Humans , Image Processing, Computer-Assisted , Magnetic Resonance Imaging , Male , Neural Pathways/physiologyABSTRACT
Patients with major depressive disorder (MDD) show heterogeneous treatment response and highly variable clinical trajectories: while some patients experience swift recovery, others show relapsing-remitting or chronic courses. Predicting individual clinical trajectories at an early stage is a key challenge for psychiatry and might facilitate individually tailored interventions. So far, however, reliable predictors at the single-patient level are absent. Here, we evaluated the utility of a machine learning strategy - generative embedding (GE) - which combines interpretable generative models with discriminative classifiers. Specifically, we used functional magnetic resonance imaging (fMRI) data of emotional face perception in 85 MDD patients from the NEtherlands Study of Depression and Anxiety (NESDA) who had been followed up over two years and classified into three subgroups with distinct clinical trajectories. Combining a generative model of effective (directed) connectivity with support vector machines (SVMs), we could predict whether a given patient would experience chronic depression vs. fast remission with a balanced accuracy of 79%. Gradual improvement vs. fast remission could still be predicted above-chance, but less convincingly, with a balanced accuracy of 61%. Generative embedding outperformed classification based on conventional (descriptive) features, such as functional connectivity or local activation estimates, which were obtained from the same data and did not allow for above-chance classification accuracy. Furthermore, predictive performance of GE could be assigned to a specific network property: the trial-by-trial modulation of connections by emotional content. Given the limited sample size of our study, the present results are preliminary but may serve as proof-of-concept, illustrating the potential of GE for obtaining clinical predictions that are interpretable in terms of network mechanisms. Our findings suggest that abnormal dynamic changes of connections involved in emotional face processing might be associated with higher risk of developing a less favorable clinical course.
Subject(s)
Brain Mapping/methods , Depressive Disorder, Major/diagnostic imaging , Image Interpretation, Computer-Assisted/methods , Support Vector Machine , Adult , Disease Progression , Female , Humans , Magnetic Resonance Imaging/methods , Male , Middle Aged , PrognosisABSTRACT
A recently introduced hierarchical generative model unified the inference of effective connectivity in individual subjects and the unsupervised identification of subgroups defined by connectivity patterns. This hierarchical unsupervised generative embedding (HUGE) approach combined a hierarchical formulation of dynamic causal modelling (DCM) for fMRI with Gaussian mixture models and relied on Markov chain Monte Carlo (MCMC) sampling for inference. While well suited for the inversion of complex hierarchical models, MCMC-based sampling suffers from a computational burden that is prohibitive for many applications. To address this problem, this paper derives an efficient variational Bayesian (VB) inversion scheme for HUGE that simultaneously provides approximations to the posterior distribution over model parameters and to the log model evidence. The face validity of the VB scheme was tested using two synthetic fMRI datasets with known ground truth. Additionally, an empirical fMRI dataset of stroke patients and healthy controls was used to evaluate the practical utility of the method in application to real-world problems. Our analyses demonstrate good performance of our VB scheme, with a marked speed-up of model inversion by two orders of magnitude compared to MCMC, while maintaining a similar level of accuracy. Notably, additional acceleration would be possible if parallel computing techniques were applied. Generally, our VB implementation of HUGE is fast enough to support multi-start procedures for whole-group analyses, a useful strategy to ameliorate problems with local extrema. HUGE thus represents a potentially useful practical solution for an important problem in clinical neuromodeling and computational psychiatry, i.e., the unsupervised detection of subgroups in heterogeneous populations that are defined by effective connectivity.
Subject(s)
Algorithms , Brain Mapping/methods , Image Processing, Computer-Assisted/methods , Models, Neurological , Adult , Aged , Bayes Theorem , Datasets as Topic , Female , Humans , Magnetic Resonance Imaging/methods , Male , Middle AgedABSTRACT
The development of whole-brain models that can infer effective (directed) connection strengths from fMRI data represents a central challenge for computational neuroimaging. A recently introduced generative model of fMRI data, regression dynamic causal modeling (rDCM), moves towards this goal as it scales gracefully to very large networks. However, large-scale networks with thousands of connections are difficult to interpret; additionally, one typically lacks information (data points per free parameter) for precise estimation of all model parameters. This paper introduces sparsity constraints to the variational Bayesian framework of rDCM as a solution to these problems in the domain of task-based fMRI. This sparse rDCM approach enables highly efficient effective connectivity analyses in whole-brain networks and does not require a priori assumptions about the network's connectivity structure but prunes fully (all-to-all) connected networks as part of model inversion. Following the derivation of the variational Bayesian update equations for sparse rDCM, we use both simulated and empirical data to assess the face validity of the model. In particular, we show that it is feasible to infer effective connection strengths from fMRI data using a network with more than 100 regions and 10,000 connections. This demonstrates the feasibility of whole-brain inference on effective connectivity from fMRI data - in single subjects and with a run-time below 1â¯min when using parallelized code. We anticipate that sparse rDCM may find useful application in connectomics and clinical neuromodeling - for example, for phenotyping individual patients in terms of whole-brain network structure.
Subject(s)
Brain/physiology , Connectome/methods , Models, Neurological , Models, Theoretical , Nerve Net/physiology , Bayes Theorem , Humans , Magnetic Resonance Imaging/methodsABSTRACT
Despite the success of modern neuroimaging techniques in furthering our understanding of cognitive and pathophysiological processes, translation of these advances into clinically relevant tools has been virtually absent until now. Neuromodeling represents a powerful framework for overcoming this translational deadlock, and the development of computational models to solve clinical problems has become a major scientific goal over the last decade, as reflected by the emergence of clinically oriented neuromodeling fields like Computational Psychiatry, Computational Neurology, and Computational Psychosomatics. Generative models of brain physiology and connectivity in the human brain play a key role in this endeavor, striving for computational assays that can be applied to neuroimaging data from individual patients for differential diagnosis and treatment prediction. In this review, we focus on dynamic causal modeling (DCM) and its use for Computational Psychiatry. DCM is a widely used generative modeling framework for functional magnetic resonance imaging (fMRI) and magneto-/electroencephalography (M/EEG) data. This article reviews the basic concepts of DCM, revisits examples where it has proven valuable for addressing clinically relevant questions, and critically discusses methodological challenges and recent methodological advances. We conclude this review with a more general discussion of the promises and pitfalls of generative models in Computational Psychiatry and highlight the path that lies ahead of us. This article is categorized under: Neuroscience > Computation Neuroscience > Clinical Neuroscience.
Subject(s)
Computational Biology/methods , Models, Neurological , Psychiatry/methods , Translational Research, Biomedical/methods , Bayes Theorem , Humans , Neuroimaging/methodsABSTRACT
The development of brain imaging techniques, in particular functional magnetic resonance imaging (fMRI), made it possible to non-invasively study the hemispheric lateralization of cognitive brain functions in large cohorts. Comprehensive models of hemispheric lateralization are, however, still missing and should not only account for the hemispheric specialization of individual brain functions, but also for the interactions among different lateralized cognitive processes (e.g., language and visuospatial processing). This calls for robust and reliable paradigms to study hemispheric lateralization for various cognitive functions. While numerous reliable imaging paradigms have been developed for language, which represents the most prominent left-lateralized brain function, the reliability of imaging paradigms investigating typically right-lateralized brain functions, such as visuospatial processing, has received comparatively less attention. In the present study, we aimed to establish an fMRI paradigm that robustly and reliably identifies right-hemispheric activation evoked by visuospatial processing in individual subjects. In a first study, we therefore compared three frequently used paradigms for assessing visuospatial processing and evaluated their utility to robustly detect right-lateralized brain activity on a single-subject level. In a second study, we then assessed the test-retest reliability of the so-called Landmark task-the paradigm that yielded the most robust results in study 1. At the single-voxel level, we found poor reliability of the brain activation underlying visuospatial attention. This suggests that poor signal-to-noise ratios can become a limiting factor for test-retest reliability. This represents a common detriment of fMRI paradigms investigating visuospatial attention in general and therefore highlights the need for careful considerations of both the possibilities and limitations of the respective fMRI paradigm-in particular, when being interested in effects at the single-voxel level. Notably, however, when focusing on the reliability of measures of hemispheric lateralization (which was the main goal of study 2), we show that hemispheric dominance (quantified by the lateralization index, LI, with |LI| >0.4) of the evoked activation could be robustly determined in more than 62% and, if considering only two categories (i.e., left, right), in more than 93% of our subjects. Furthermore, the reliability of the lateralization strength (LI) was "fair" to "good". In conclusion, our results suggest that the degree of right-hemispheric dominance during visuospatial processing can be reliably determined using the Landmark task, both at the group and single-subject level, while at the same time stressing the need for future refinements of experimental paradigms and more sophisticated fMRI data acquisition techniques.