Replicating extensive brain structural heterogeneity in individuals with schizophrenia and bipolar disorder.

Identifying brain processes involved in the risk and development of mental disorders is a major aim. We recently reported substantial interindividual heterogeneity in brain structural aberrations among patients with schizophrenia and bipolar disorder. Estimating the normative range of voxel-based morphometry (VBM) data among healthy individuals using a Gaussian process regression (GPR) enables us to map individual deviations from the healthy range in unseen datasets. Here, we aim to replicate our previous results in two independent samples of patients with schizophrenia (n1 = 94; n2 = 105), bipolar disorder (n1 = 116; n2 = 61), and healthy individuals (n1 = 400; n2 = 312). In line with previous findings with exception of the cerebellum our results revealed robust group level differences between patients and healthy individuals, yet only a small proportion of patients with schizophrenia or bipolar disorder exhibited extreme negative deviations from normality in the same brain regions. These direct replications support that group level-differences in brain structure disguise considerable individual differences in brain aberrations, with important implications for the interpretation and generalization of group-level brain imaging findings to the individual with a mental disorder.

Population-wide cerebellar growth models of children and adolescents.

In the past, the cerebellum has been best known for its crucial role in motor function. However, increasingly more findings highlight the importance of cerebellar contributions in cognitive functions and neurodevelopment. Using a total of 7240 neuroimaging scans from 4862 individuals, we describe and provide detailed, openly available models of cerebellar development in childhood and adolescence (age range: 6-17 years), an important time period for brain development and onset of neuropsychiatric disorders. Next to a traditionally used anatomical parcellation of the cerebellum, we generated growth models based on a recently proposed functional parcellation. In both, we find an anterior-posterior growth gradient mirroring the age-related improvements of underlying behavior and function, which is analogous to cerebral maturation patterns and offers evidence for directly related cerebello-cortical developmental trajectories. Finally, we illustrate how the current approach can be used to detect cerebellar abnormalities in clinical samples.

A 10-Year Longitudinal Study of Brain Cortical Thickness in People with First-Episode Psychosis Using Normative Models.

BACKGROUND: Clinical forecasting models have potential to optimize treatment and improve outcomes in psychosis, but predicting long-term outcomes is challenging and long-term follow-up data are scarce. In this 10-year longitudinal study, we aimed to characterize the temporal evolution of cortical correlates of psychosis and their associations with symptoms. DESIGN: Structural magnetic resonance imaging (MRI) from people with first-episode psychosis and controls (n = 79 and 218) were obtained at enrollment, after 12 months (n = 67 and 197), and 10 years (n = 23 and 77), within the Thematically Organized Psychosis (TOP) study. Normative models for cortical thickness estimated on public MRI datasets (n = 42 983) were applied to TOP data to obtain deviation scores for each region and timepoint. Positive and Negative Syndrome Scale (PANSS) scores were acquired at each timepoint along with registry data. Linear mixed effects models assessed effects of diagnosis, time, and their interactions on cortical deviations plus associations with symptoms. RESULTS: LMEs revealed conditional main effects of diagnosis and time × diagnosis interactions in a distributed cortical network, where negative deviations in patients attenuate over time. In patients, symptoms also attenuate over time. LMEs revealed effects of anterior cingulate on PANSS total, and insular and orbitofrontal regions on PANSS negative scores. CONCLUSIONS: This long-term longitudinal study revealed a distributed pattern of cortical differences which attenuated over time together with a reduction in symptoms. These findings are not in line with a simple neurodegenerative account of schizophrenia, and deviations from normative models offer a promising avenue to develop biomarkers to track clinical trajectories over time.

Computational Modeling of the n-Back Task in the ABCD Study: Associations of Drift Diffusion Model Parameters to Polygenic Scores of Mental Disorders and Cardiometabolic Diseases.

BACKGROUND: Cognitive dysfunction is common in mental disorders and represents a potential risk factor in childhood. The nature and extent of associations between childhood cognitive function and polygenic risk for mental disorders is unclear. We applied computational modeling to gain insight into mechanistic processes underlying decision making and working memory in childhood and their associations with polygenic risk scores (PRSs) for mental disorders and comorbid cardiometabolic diseases. METHODS: We used the drift diffusion model to infer latent computational processes underlying decision making and working memory during the n-back task in 3707 children ages 9 to 10 years from the Adolescent Brain Cognitive Development (ABCD) Study. Single nucleotide polymorphism-based heritability was estimated for cognitive phenotypes, including computational parameters, aggregated n-back task performance, and neurocognitive assessments. PRSs were calculated for Alzheimer's disease, bipolar disorder, coronary artery disease (CAD), major depressive disorder, obsessive-compulsive disorder, schizophrenia, and type 2 diabetes. RESULTS: Heritability estimates of cognitive phenotypes ranged from 12% to 38%. Bayesian mixed models revealed that slower accumulation of evidence was associated with higher PRSs for CAD and schizophrenia. Longer nondecision time was associated with higher PRSs for Alzheimer's disease and lower PRSs for CAD. Narrower decision threshold was associated with higher PRSs for CAD. Load-dependent effects on nondecision time and decision threshold were associated with PRSs for Alzheimer's disease and CAD, respectively. Aggregated neurocognitive test scores were not associated with PRSs for any of the mental or cardiometabolic phenotypes. CONCLUSIONS: We identified distinct associations between computational cognitive processes and genetic risk for mental illness and cardiometabolic disease, which could represent childhood cognitive risk factors.

Estimating multimodal brain variability in schizophrenia spectrum disorders: A worldwide ENIGMA study.

Objective: Schizophrenia is a multifaceted disorder associated with structural brain heterogeneity. Despite its relevance for identifying illness subtypes and informative biomarkers, structural brain heterogeneity in schizophrenia remains incompletely understood. Therefore, the objective of this study was to provide a comprehensive insight into the structural brain heterogeneity associated with schizophrenia. Methods: This meta- and mega-analysis investigated the variability of multimodal structural brain measures of white and gray matter in individuals with schizophrenia versus healthy controls. Using the ENIGMA dataset of MRI-based brain measures from 22 international sites with up to 6139 individuals for a given brain measure, we examined variability in cortical thickness, surface area, folding index, subcortical volume and fractional anisotropy. Results: We found that individuals with schizophrenia are distinguished by higher heterogeneity in the frontotemporal network with regard to multimodal structural measures. Moreover, individuals with schizophrenia showed higher homogeneity of the folding index, especially in the left parahippocampal region. Conclusions: Higher multimodal heterogeneity in frontotemporal regions potentially implies different subtypes of schizophrenia that converge on impaired frontotemporal interaction as a core feature of the disorder. Conversely, more homogeneous folding patterns in the left parahippocampal region might signify a consistent characteristic of schizophrenia shared across subtypes. These findings underscore the importance of structural brain variability in advancing our neurobiological understanding of schizophrenia, and aid in identifying illness subtypes as well as informative biomarkers.

Deviations from normative brain white and gray matter structure are associated with psychopathology in youth.

Combining imaging modalities and metrics that are sensitive to various aspects of brain structure and maturation may help identify individuals that show deviations in relation to same-aged peers, and thus benefit early-risk-assessment for mental disorders. We used one timepoint multimodal brain imaging, cognitive, and questionnaire data from 1280 eight- to twenty-one-year-olds from the Philadelphia Neurodevelopmental Cohort. We estimated age-related gray and white matter properties and estimated individual deviation scores using normative modeling. Next, we tested for associations between the estimated deviation scores, and with psychopathology domain scores and cognition. More negative deviations in DTI-based fractional anisotropy (FA) and the first principal eigenvalue of the diffusion tensor (L1) were associated with higher scores on psychosis positive and prodromal symptoms and general psychopathology. A more negative deviation in cortical thickness (CT) was associated with a higher general psychopathology score. Negative deviations in global FA, surface area, L1 and CT were also associated with poorer cognitive performance. No robust associations were found between the deviation scores based on CT and DTI. The low correlations between the different multimodal magnetic resonance imaging-based deviation scores suggest that psychopathological burden in adolescence can be mapped onto partly distinct neurobiological features.

Closing the life-cycle of normative modeling using federated hierarchical Bayesian regression.

Clinical neuroimaging data availability has grown substantially in the last decade, providing the potential for studying heterogeneity in clinical cohorts on a previously unprecedented scale. Normative modeling is an emerging statistical tool for dissecting heterogeneity in complex brain disorders. However, its application remains technically challenging due to medical data privacy issues and difficulties in dealing with nuisance variation, such as the variability in the image acquisition process. Here, we approach the problem of estimating a reference normative model across a massive population using a massive multi-center neuroimaging dataset. To this end, we introduce a federated probabilistic framework using hierarchical Bayesian regression (HBR) to complete the life-cycle of normative modeling. The proposed model provides the possibilities to learn, update, and adapt the model parameters on decentralized neuroimaging data. Our experimental results confirm the superiority of HBR in deriving more accurate normative ranges on large multi-site neuroimaging datasets compared to the current standard methods. In addition, our approach provides the possibility to recalibrate and reuse the learned model on local datasets and even on datasets with very small sample sizes. The proposed method will facilitate applications of normative modeling as a medical tool for screening the biological deviations in individuals affected by complex illnesses such as mental disorders.

The normative modeling framework for computational psychiatry.

Normative modeling is an emerging and innovative framework for mapping individual differences at the level of a single subject or observation in relation to a reference model. It involves charting centiles of variation across a population in terms of mappings between biology and behavior, which can then be used to make statistical inferences at the level of the individual. The fields of computational psychiatry and clinical neuroscience have been slow to transition away from patient versus 'healthy' control analytic approaches, probably owing to a lack of tools designed to properly model biological heterogeneity of mental disorders. Normative modeling provides a solution to address this issue and moves analysis away from case-control comparisons that rely on potentially noisy clinical labels. Here we define a standardized protocol to guide users through, from start to finish, normative modeling analysis using the Predictive Clinical Neuroscience toolkit (PCNtoolkit). We describe the input data selection process, provide intuition behind the various modeling choices and conclude by demonstrating several examples of downstream analyses that the normative model may facilitate, such as stratification of high-risk individuals, subtyping and behavioral predictive modeling. The protocol takes ~1-3 h to complete.

Charting brain growth and aging at high spatial precision.

Defining reference models for population variation, and the ability to study individual deviations is essential for understanding inter-individual variability and its relation to the onset and progression of medical conditions. In this work, we assembled a reference cohort of neuroimaging data from 82 sites (N=58,836; ages 2-100) and used normative modeling to characterize lifespan trajectories of cortical thickness and subcortical volume. Models are validated against a manually quality checked subset (N=24,354) and we provide an interface for transferring to new data sources. We showcase the clinical value by applying the models to a transdiagnostic psychiatric sample (N=1985), showing they can be used to quantify variability underlying multiple disorders whilst also refining case-control inferences. These models will be augmented with additional samples and imaging modalities as they become available. This provides a common reference platform to bind results from different studies and ultimately paves the way for personalized clinical decision-making.

How does the presence of neural probes affect extracellular potentials?

OBJECTIVE: Mechanistic modeling of neurons is an essential component of computational neuroscience that enables scientists to simulate, explain, and explore neural activity. The conventional approach to simulation of extracellular neural recordings first computes transmembrane currents using the cable equation and then sums their contribution to model the extracellular potential. This two-step approach relies on the assumption that the extracellular space is an infinite and homogeneous conductive medium, while measurements are performed using neural probes. The main purpose of this paper is to assess to what extent the presence of the neural probes of varying shape and size impacts the extracellular field and how to correct for them. APPROACH: We apply a detailed modeling framework allowing explicit representation of the neuron and the probe to study the effect of the probes and thereby estimate the effect of ignoring it. We use meshes with simplified neurons and different types of probe and compare the extracellular action potentials with and without the probe in the extracellular space. We then compare various solutions to account for the probes' presence and introduce an efficient probe correction method to include the probe effect in modeling of extracellular potentials. MAIN RESULTS: Our computations show that microwires hardly influence the extracellular electric field and their effect can therefore be ignored. In contrast, multi-electrode arrays (MEAs) significantly affect the extracellular field by magnifying the recorded potential. While MEAs behave similarly to infinite insulated planes, we find that their effect strongly depends on the neuron-probe alignment and probe orientation. SIGNIFICANCE: Ignoring the probe effect might be deleterious in some applications, such as neural localization and parameterization of neural models from extracellular recordings. Moreover, the presence of the probe can improve the interpretation of extracellular recordings, by providing a more accurate estimation of the extracellular potential generated by neuronal models.

Functional Relevance of Different Basal Ganglia Pathways Investigated in a Spiking Model with Reward Dependent Plasticity.

The brain enables animals to behaviorally adapt in order to survive in a complex and dynamic environment, but how reward-oriented behaviors are achieved and computed by its underlying neural circuitry is an open question. To address this concern, we have developed a spiking model of the basal ganglia (BG) that learns to dis-inhibit the action leading to a reward despite ongoing changes in the reward schedule. The architecture of the network features the two pathways commonly described in BG, the direct (denoted D1) and the indirect (denoted D2) pathway, as well as a loop involving striatum and the dopaminergic system. The activity of these dopaminergic neurons conveys the reward prediction error (RPE), which determines the magnitude of synaptic plasticity within the different pathways. All plastic connections implement a versatile four-factor learning rule derived from Bayesian inference that depends upon pre- and post-synaptic activity, receptor type, and dopamine level. Synaptic weight updates occur in the D1 or D2 pathways depending on the sign of the RPE, and an efference copy informs upstream nuclei about the action selected. We demonstrate successful performance of the system in a multiple-choice learning task with a transiently changing reward schedule. We simulate lesioning of the various pathways and show that a condition without the D2 pathway fares worse than one without D1. Additionally, we simulate the degeneration observed in Parkinson's disease (PD) by decreasing the number of dopaminergic neurons during learning. The results suggest that the D1 pathway impairment in PD might have been overlooked. Furthermore, an analysis of the alterations in the synaptic weights shows that using the absolute reward value instead of the RPE leads to a larger change in D1.

Optogenetic stimulation in a computational model of the basal ganglia biases action selection and reward prediction error.

Optogenetic stimulation of specific types of medium spiny neurons (MSNs) in the striatum has been shown to bias the selection of mice in a two choices task. This shift is dependent on the localisation and on the intensity of the stimulation but also on the recent reward history. We have implemented a way to simulate this increased activity produced by the optical flash in our computational model of the basal ganglia (BG). This abstract model features the direct and indirect pathways commonly described in biology, and a reward prediction pathway (RP). The framework is similar to Actor-Critic methods and to the ventral/dorsal distinction in the striatum. We thus investigated the impact on the selection caused by an added stimulation in each of the three pathways. We were able to reproduce in our model the bias in action selection observed in mice. Our results also showed that biasing the reward prediction is sufficient to create a modification in the action selection. However, we had to increase the percentage of trials with stimulation relative to that in experiments in order to impact the selection. We found that increasing only the reward prediction had a different effect if the stimulation in RP was action dependent (only for a specific action) or not. We further looked at the evolution of the change in the weights depending on the stage of learning within a block. A bias in RP impacts the plasticity differently depending on that stage but also on the outcome. It remains to experimentally test how the dopaminergic neurons are affected by specific stimulations of neurons in the striatum and to relate data to predictions of our model.

Action selection performance of a reconfigurable basal ganglia inspired model with Hebbian-Bayesian Go-NoGo connectivity.

Several studies have shown a strong involvement of the basal ganglia (BG) in action selection and dopamine dependent learning. The dopaminergic signal to striatum, the input stage of the BG, has been commonly described as coding a reward prediction error (RPE), i.e., the difference between the predicted and actual reward. The RPE has been hypothesized to be critical in the modulation of the synaptic plasticity in cortico-striatal synapses in the direct and indirect pathway. We developed an abstract computational model of the BG, with a dual pathway structure functionally corresponding to the direct and indirect pathways, and compared its behavior to biological data as well as other reinforcement learning models. The computations in our model are inspired by Bayesian inference, and the synaptic plasticity changes depend on a three factor Hebbian-Bayesian learning rule based on co-activation of pre- and post-synaptic units and on the value of the RPE. The model builds on a modified Actor-Critic architecture and implements the direct (Go) and the indirect (NoGo) pathway, as well as the reward prediction (RP) system, acting in a complementary fashion. We investigated the performance of the model system when different configurations of the Go, NoGo, and RP system were utilized, e.g., using only the Go, NoGo, or RP system, or combinations of those. Learning performance was investigated in several types of learning paradigms, such as learning-relearning, successive learning, stochastic learning, reversal learning and a two-choice task. The RPE and the activity of the model during learning were similar to monkey electrophysiological and behavioral data. Our results, however, show that there is not a unique best way to configure this BG model to handle well all the learning paradigms tested. We thus suggest that an agent might dynamically configure its action selection mode, possibly depending on task characteristics and also on how much time is available.