Publications

Functional neurosurgery has afforded the opportunity to assess interactions between populations of neurons in the human cerebral cortex and basal ganglia in patients with Parkinson's disease (PD). Interactions occur over a wide range of frequencies, and the functional significance of those >30 Hz is particularly unclear. Do they improve movement, and, if so, in what way? We acquired simultaneously magnetoencephalography and direct recordings from the subthalamic nucleus (STN) in 17 PD patients. We examined the effect of synchronous and sequential finger movements and of the dopamine prodrug levodopa on induced power in the contralateral primary motor cortex (M1) and STN and on the coherence between the two structures. We observed discrete peaks in M1 and STN power at 60-90 Hz and at 300-400 Hz. All these power peaks increased with movement and levodopa treatment. Only STN activity at 60-90 Hz was coherent with activity in M1. Directionality analysis showed that STN gamma activity at 60-90 Hz tended to drive gamma activity in M1. The effects of levodopa on both local and distant synchronization at 60-90 Hz correlated with the degree of improvement in bradykinesia-rigidity as did local STN activity at 300-400 Hz. Despite this, there were no effects of movement type, nor interactions between movement type and levodopa in the STN, nor in the coherence between STN and M1. We conclude that synchronization at 60-90 Hz in the basal ganglia cortical network is prokinetic but likely through a modulatory effect rather than any involvement in explicit motor processing.

Choosing actions that result in advantageous outcomes is a fundamental function of nervous systems. All computational decision-making models contain a mechanism that controls the variability of (or confidence in) action selection, but its neural implementation is unclear-especially in humans. We investigated this mechanism using two influential decision-making frameworks: active inference (AI) and reinforcement learning (RL). In AI, the precision (inverse variance) of beliefs about policies controls action selection variability-similar to decision 'noise' parameters in RL-and is thought to be encoded by striatal dopamine signaling. We tested this hypothesis by administering a 'go/no-go' task to 75 healthy participants, and measuring striatal dopamine 2/3 receptor (D2/3R) availability in a subset (n = 25) using [11C]-(+)-PHNO positron emission tomography. In behavioral model comparison, RL performed best across the whole group but AI performed best in participants performing above chance levels. Limbic striatal D2/3R availability had linear relationships with AI policy precision (P = 0.029) as well as with RL irreducible decision 'noise' (P = 0.020), and this relationship with D2/3R availability was confirmed with a 'decision stochasticity' factor that aggregated across both models (P = 0.0006). These findings are consistent with occupancy of inhibitory striatal D2/3Rs decreasing the variability of action selection in humans.

The default mode network (DMN) is thought to reflect endogenous neural activity, which is considered as one of the most intriguing phenomena in cognitive neuroscience. Previous studies have found that key regions within the DMN are highly interconnected. Here, we characterized the genetic influences on causal or directed information flow within the DMN during the resting state. In this study, we recruited 46 pairs of twins and collected fMRI imaging data using a 3.0 T scanner. Dynamic causal modeling was conducted for each participant, and a structural equation model was used to calculate the heritability of DMN in terms of its effective connectivity. Model comparison favored a full-connected model. Structural equal modeling was used to estimate the additive genetics (A), common environment (C) and unique environment (E) contributions to variance for the DMN effective connectivity. The ACE model was preferred in the comparison of structural equation models. Heritability of DMN effective connectivity was 0.54, suggesting that the genetic made a greater contribution to the effective connectivity within DMN. Establishing the heritability of default-mode effective connectivity endorses the use of resting-state networks as endophenotypes or intermediate phenotypes in the search for the genetic basis of psychiatric or neurological illnesses.

Recently, there has been a lot of interest in characterising the connectivity of resting state brain networks. Most of the literature uses functional connectivity to examine these intrinsic brain networks. Functional connectivity has well documented limitations because of its inherent inability to identify causal interactions. Dynamic causal modelling (DCM) is a framework that allows for the identification of the causal (directed) connections among neuronal systems — known as effective connectivity. This technical note addresses the validity of a recently proposed DCM for resting state fMRI – as measured in terms of their complex cross spectral density – referred to as spectral DCM. Spectral DCM differs from (the alternative) stochastic DCM by parameterising neuronal fluctuations using scale free (i.e., power law) forms, rendering the stochastic model of neuronal activity deterministic. Spectral DCM not only furnishes an efficient estimation of model parameters but also enables the detection of group differences in effective connectivity, the form and amplitude of the neuronal fluctuations or both. We compare and contrast spectral and stochastic DCM models with endogenous fluctuations or state noise on hidden states. We used simulated data to first establish the face validity of both schemes and show that they can recover the model (and its parameters) that generated the data. We then used Monte Carlo simulations to assess the accuracy of both schemes in terms of their root mean square error. We also simulated group differences and compared the ability of spectral and stochastic DCMs to identify these differences. We show that spectral DCM was not only more accurate but also more sensitive to group differences. Finally, we performed a comparative evaluation using real resting state fMRI data (from an open access resource) to study the functional integration within default mode network using spectral and stochastic DCMs.

SummaryThe neural systems involved in hearing and repeating single words were investigated in a series of experiments using PET. Neuropsychological and psycholinguistic studies implicate the involvement of posterior and anterior left perisylvian regions (Wernicke's and Broca's areas). Although previous functional neuroimaging studies have consistently shown activation of Wernicke's area, there has been only variable implication of Broca's area. This study demonstrates that Broca's area is involved in both auditory word perception and repetition but activation is dependent on task (greater during repetition than hearing) and stimulus presentation (greater when hearing words at a slow rate). The peak of frontal activation in response to hearing words is anterior to that associated with repeating words; the former is probably located in Brodmann's area 45, the latter in Brodmann's area 44 and the adjacent precentral sulcus. As Broca's area activation is more subtle and complex than that in Wernicke's area during these tasks, the likelihood of observing it is influenced by both the study design and the image analysis technique employed. As a secondary outcome from the study, the response of bilateral auditory association cortex to 'own voice' during repetition was shown to be the same as when listening to 'other voice' from a prerecorded tape.

This technical note presents a dynamic causal modelling (DCM) procedure for evaluating different models of neurovascular coupling in the human brain - using combined electromagnetic (M/EEG) and functional magnetic resonance imaging (fMRI) data. This procedure compares the evidence for biologically informed models of neurovascular coupling using Bayesian model comparison. First, fMRI data are used to localise regionally specific neuronal responses. The coordinates of these responses are then used as the location priors in a DCM of electrophysiological responses elicited by the same paradigm. The ensuing estimates of model parameters are then used to generate neuronal drive functions, which model pre- or post-synaptic activity for each experimental condition. These functions form the input to a model of neurovascular coupling, whose parameters are estimated from the fMRI data. Crucially, this enables one to evaluate different models of neurovascular coupling, using Bayesian model comparison - asking, for example, whether instantaneous or delayed, pre- or post-synaptic signals mediate haemodynamic responses. We provide an illustrative application of the procedure using a single-subject auditory fMRI and MEG dataset. The code and exemplar data accompanying this technical note are available through the statistical parametric mapping (SPM) software.

Abstract As with the Spanish Flu a century ago, authorities have responded to the current COVID-19 pandemic with extraordinary public health measures. In particular, lockdown and related social distancing policies are motivated in some countries by the need to slow virus propagation—so that the primary wave of patients suffering from severe forms of COVID infection do not exceed the capacity of intensive care units. But unlocking poses a critical issue because relaxing social distancing may, in principle, generate secondary waves. Ironically however, the dynamic repertoire of established epidemiological models that support this kind of reasoning is limited to single epidemic outbreaks. In turn, predictions regarding secondary waves are tautologically derived from imposing assumptions about changes in the so-called “effective reproduction number”. In this work, we depart from this approach and extend the LIST (Location-Infection-Symptom-Testing) model of the COVID pandemic with realistic nonlinear feedback mechanisms that under certain conditions, cause lockdown-induced secondary outbreaks. The original LIST model captures adaptive social distancing, i . e . the transient reduction of the number of person-to-person contacts (and hence the rate of virus transmission), as a societal response to salient public health risks. Here, we consider the possibility that such pruning of socio-geographical networks may also temporarily isolate subsets of local populations from the virus. Crucially however, such unreachable people will become susceptible again when adaptive social distancing relaxes and the density of contacts within socio-geographical networks increases again. Taken together, adaptive social distancing and network unreachability thus close a nonlinear feedback loop that endows the LIST model with a mechanism that can generate autonomous (lockdown-induced) secondary waves. However, whether and how secondary waves arise depend upon the interaction with other nonlinear mechanisms that capture other forms of transmission heterogeneity. We apply the ensuing LIST model to numerical simulations and exhaustive analyses of regional French epidemiological data. In brief, we find evidence for this kind of nonlinear feedback mechanism in the empirical dynamics of the pandemic in France. However, rather than generating catastrophic secondary outbreaks (as is typically assumed), the model predicts that the impact of lockdown-induced variations in population susceptibility and transmission may eventually reduce to a steady-state endemic equilibrium with a low but stable infection rate.

A bstract 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. S ignificance S tatement 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 neuro-modelling 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.

Prominent accounts of sentient behaviour depict brains as generative models of organismic interaction with the world, raising points of contact with current work in Generative AI. However, because they contend with the control of purposive, life-maintaining sensorimotor interactions, the generative models of living organisms are inextricably anchored to the body and world. Unlike the passive models learnt by Generative AIs, they must capture and control the sensory consequences of action. This allows embodied agents to intervene upon their worlds in ways that constantly put their best models to the test; providing a solid bedrock that is—we argue—essential to the development of genuine understanding. Here, we review the resulting implications, and consider future directions of travel for Generative AI.

Huntington's disease is an incurable neurodegenerative disease caused by inheritance of an expanded cytosine-adenine-guanine (CAG) trinucleotide repeat within the Huntingtin gene. Extensive volume loss and altered diffusion metrics in the basal ganglia, cortex and white matter are seen when patients with Huntington's disease (HD) undergo structural imaging, suggesting that changes in basal ganglia-cortical structural connectivity occur. The aims of this study were to characterise altered patterns of basal ganglia-cortical structural connectivity with high anatomical precision in premanifest and early manifest HD, and to identify associations between structural connectivity and genetic or clinical markers of HD. 3-Tesla diffusion tensor magnetic resonance images were acquired from 14 early manifest HD subjects, 17 premanifest HD subjects and 18 controls. Voxel-based analyses of probabilistic tractography were used to quantify basal ganglia-cortical structural connections. Canonical variate analysis was used to demonstrate disease-associated patterns of altered connectivity and to test for associations between connectivity and genetic and clinical markers of HD; this is the first study in which such analyses have been used. Widespread changes were seen in basal ganglia-cortical structural connectivity in early manifest HD subjects; this has relevance for development of therapies targeting the striatum. Premanifest HD subjects had a pattern of connectivity more similar to that of controls, suggesting progressive change in connections over time. Associations between structural connectivity patterns and motor and cognitive markers of disease severity were present in early manifest subjects. Our data suggest the clinical phenotype in manifest HD may be at least partly a result of altered connectivity.