Neuropharmacological computational analysis of longitudinal electroencephalograms in clozapine-treated patients with schizophrenia using hierarchical dynamic causal modeling.

The hierarchical characteristics of the brain are prominent in the pharmacological treatment of psychiatric diseases, primarily targeting cellular receptors that extend upward to intrinsic connectivity within a region, interregional connectivity, and, consequently, clinical observations such as an electroencephalogram (EEG). To understand the long-term effects of neuropharmacological intervention on neurobiological properties at different hierarchical levels, we explored long-term changes in neurobiological parameters of an N-methyl-D-aspartate canonical microcircuit model (CMM-NMDA) in the default mode network (DMN) and auditory hallucination network (AHN) using dynamic causal modeling of longitudinal EEG in clozapine-treated patients with schizophrenia. The neurobiological properties of the CMM-NMDA model associated with symptom improvement in schizophrenia were found across hierarchical levels, from a reduced membrane capacity of the deep pyramidal cell and intrinsic connectivity with the inhibitory population in DMN and intrinsic and extrinsic connectivity in AHN. The medication duration mainly affects the intrinsic connectivity and NMDA time constant in DMN. Virtual perturbation analysis specified the contribution of each parameter to the cross-spectral density (CSD) of the EEG, particularly intrinsic connectivity and membrane capacitances for CSD frequency shifts and progression. It further reveals that excitatory and inhibitory connectivity complements frequency-specific CSD changes, notably the alpha frequency band in DMN. Positive and negative synergistic interactions exist between neurobiological properties primarily within the same region in patients treated with clozapine. The current study shows how computational neuropharmacology helps explore the multiscale link between neurobiological properties and clinical observations and understand the long-term mechanism of neuropharmacological intervention reflected in clinical EEG.

A computational framework for optimal control of a self-adjustive neural system with activity-dependent and homeostatic plasticity.

The control of the brain system has received increasing attention in the domain of brain science. Most brain control studies have been conducted to explore the brain network's graph-theoretic properties or to produce the desired state based on neural state dynamics, regarding the brain as a passively responding system. However, the self-adjusting nature of neural system after treatment has not been fully considered in the brain control. In the present study, we propose a computational framework for optimal control of the brain with a self-adjustment process in the effective connectivity after treatment. The neural system is modeled to adjust its outgoing effective connectivity as activity-dependent plasticity after treatment, followed by synaptic rescaling of incoming effective connectivity. To control this neural system to induce the desired function, the system's self-adjustment parameter is first estimated, based on which the treatment is optimized. Utilizing this framework, we conducted simulations of optimal control over a functional hippocampal circuitry, estimated using dynamic causal modeling of voltage-sensitive dye imaging from the wild type and mutant mice, responding to consecutive electrical stimuli. Simulation results for optimal control of the abnormal circuit toward a healthy circuit using a single node treatment, neural-type specific treatment as an analogy of medication, and combined treatments of medication and nodal treatment suggest the plausibility of the current framework in controlling the self-adjusting neural system within a restricted treatment setting. We believe the proposed computational framework of the self-adjustment system would help optimal control of the dynamic brain after treatment.

Empirical Bayes estimation of pairwise maximum entropy model for nonlinear brain state dynamics.

The pairwise maximum entropy model (pMEM) has recently gained widespread attention to exploring the nonlinear characteristics of brain state dynamics observed in resting-state functional magnetic resonance imaging (rsfMRI). Despite its unique advantageous features, the practical application of pMEM for individuals is limited as it requires a much larger sample than conventional rsfMRI scans. Thus, this study proposes an empirical Bayes estimation of individual pMEM using the variational expectation-maximization algorithm (VEM-MEM). The performance of the VEM-MEM is evaluated for several simulation setups with various sample sizes and network sizes. Unlike conventional maximum likelihood estimation procedures, the VEM-MEM can reliably estimate the individual model parameters, even with small samples, by effectively incorporating the group information as the prior. As a test case, the individual rsfMRI of children with attention deficit hyperactivity disorder (ADHD) is analyzed compared to that of typically developed children using the default mode network, executive control network, and salient network, obtained from the Healthy Brain Network database. We found that the nonlinear dynamic properties uniquely established on the pMEM differ for each group. Furthermore, pMEM parameters are more sensitive to group differences and are better associated with the behavior scores of ADHD compared to the Pearson correlation-based functional connectivity. The simulation and experimental results suggest that the proposed method can reliably estimate the individual pMEM and characterize the dynamic properties of individuals by utilizing empirical information of the group brain state dynamics.

Dynamic causal modeling of hippocampal activity measured via mesoscopic voltage-sensitive dye imaging.

The aim of this paper is to present a dynamic causal modeling (DCM) framework for hippocampal activity measured via voltage-sensitive dye imaging (VSDI). We propose a DCM model of the hippocampus that summarizes interactions between the hilus, CA3 and CA1 regions. The activity of each region is governed via a neuronal mass model with two inhibitory and one/two excitatory neuronal populations, which can be linked to measurement VSDI by scaling neuronal activity. To optimize the model structure for the hippocampus, we propose two Bayesian schemes: Bayesian hyperparameter optimization to estimate the unknown electrophysiological properties necessary for constructing a mesoscopic hippocampus model; and Bayesian model reduction to determine the parameterization of neural properties, and to test and include potential connections (morphologically inferred without direct evidence yet) in the model by evaluating group-level model evidence. The proposed method was applied to model spatiotemporal patterns of accumulative responses to consecutive stimuli in separate groups of wild-type mice and epileptic aristaless-related homeobox gene (Arx) conditional knock-out mutant mice (Arx-/+;Dlx5/6CRE-IRES-GFP) in order to identify group differences in the effective connectivity within the hippocampus. The causal role of each group-differing connectivity in generating mutant-like responses was further tested. The group-level analysis identified altered intra- and inter-regional effective connectivity, some of which are crucial for explaining mutant-like responses. The modelling results for the hippocampal activity suggest the plausibility of the proposed mesoscopic hippocampus model and the usefulness of utilizing the Bayesian framework for model construction in the mesoscale modeling of neural interactions using DCM.

Differences in intrinsic functional networks in patients with essential tremor who had good and poor long-term responses after thalamotomy performed using MR-guided ultrasound.

OBJECTIVE: Thalamotomy at the nucleus ventralis intermedius using MR-guided focused ultrasound has been an effective treatment method for essential tremor (ET). However, this is not true for all cases, even for successful ablation. How the brain differs in patients with ET between those with long-term good and poor outcomes is not clear. To analyze the functional connectivity difference between patients in whom thalamotomy was effective and those in whom thalamotomy was ineffective and its prognostic role in ET treatment, the authors evaluated preoperative resting-state functional MRI in thalamotomy-treated patients. METHODS: Preoperative resting-state functional MRI data in 85 patients with ET, who were experiencing tremor relief at the time of treatment and were followed up for a minimum of 6 months after the procedure, were collected for the study. The authors conducted a graph independent component analysis of the functional connectivity matrices of tremor-related networks. The patients were divided into thalamotomy-effective and thalamotomy-ineffective groups (thalamotomy-effective group, ≥ 50% motor symptom reduction; thalamotomy-ineffective group, < 50% motor symptom reduction at 6 months after treatment) and the authors compared network components between groups. RESULTS: Seventy-two (84.7%) of the 85 patients showed ≥ 50% tremor reduction from baseline at 6 months after thalamotomy. The network analysis shows significant suppression of functional network components with connections between the areas of the cerebellum and the basal ganglia and thalamus, but enhancement of those between the premotor cortex and supplementary motor area in the noneffective group compared to the effective group. CONCLUSIONS: The present study demonstrates that patients in the noneffective group have suppressed functional subnetworks in the cerebellum and subcortex regions and have enhanced functional subnetworks among motor-sensory cortical networks compared to the thalamotomy-effective group. Therefore, the authors suggest that the functional connectivity pattern might be a possible predictive factor for outcomes of MR-guided focused ultrasound thalamotomy.

Inter-species cortical registration between macaques and humans using a functional network property under a spherical demons framework.

Systematic evaluation of cortical differences between humans and macaques calls for inter-species registration of the cortex that matches homologous regions across species. For establishing homology across brains, structural landmarks and biological features have been used without paying sufficient attention to functional homology. The present study aimed to determine functional homology between the human and macaque cortices, defined in terms of functional network properties, by proposing an iterative functional network-based registration scheme using surface-based spherical demons. The functional connectivity matrix of resting-state functional magnetic resonance imaging (rs-fMRI) among cortical parcellations was iteratively calculated for humans and macaques. From the functional connectivity matrix, the functional network properties such as principal network components were derived to estimate a deformation field between the human and macaque cortices. The iterative registration procedure updates the parcellation map of macaques, corresponding to the human connectome project's multimodal parcellation atlas, which was used to derive the macaque's functional connectivity matrix. To test the plausibility of the functional network-based registration, we compared cortical registration using structural versus functional features in terms of cortical regional areal change. We also evaluated the interhemispheric asymmetry of regional area and its inter-subject variability in humans and macaques as an indirect validation of the proposed method. Higher inter-subject variability and interhemispheric asymmetry were found in functional homology than in structural homology, and the assessed asymmetry and variations were higher in humans than in macaques. The results emphasize the significance of functional network-based cortical registration across individuals within a species and across species.

State-Dependent Effective Connectivity in Resting-State fMRI.

The human brain at rest exhibits intrinsic dynamics transitioning among the multiple metastable states of the inter-regional functional connectivity. Accordingly, the demand for exploring the state-specific functional connectivity increases for a deeper understanding of mental diseases. Functional connectivity, however, lacks information about the directed causal influences among the brain regions, called effective connectivity. This study presents the dynamic causal modeling (DCM) framework to explore the state-dependent effective connectivity using spectral DCM for the resting-state functional MRI (rsfMRI). We established the sequence of brain states using the hidden Markov model with the multivariate autoregressive coefficients of rsfMRI, summarizing the functional connectivity. We decomposed the state-dependent effective connectivity using a parametric empirical Bayes scheme that models the effective connectivity of consecutive windows with the time course of the discrete states as regressors. We showed the plausibility of the state-dependent effective connectivity analysis in a simulation setting. To test the clinical applicability, we applied the proposed method to characterize the state- and subtype-dependent effective connectivity of the default mode network in children with combined-type attention deficit hyperactivity disorder (ADHD-C) compared with age-matched, typically developed children (TDC). All 88 children were subtyped according to the occupation times (i.e., dwell times) of the three dominant functional connectivity states, independently of clinical diagnosis. The state-dependent effective connectivity differences between ADHD-C and TDC according to the subtypes and those between the subtypes of ADHD-C were expressed mainly in self-inhibition, magnifying the importance of excitation inhibition balance in the subtyping. These findings provide a clear motivation for decomposing the state-dependent dynamic effective connectivity and state-dependent analysis of the directed coupling in exploring mental diseases.

Structural Brain Connectivity Constrains within-a-Day Variability of Direct Functional Connectivity.

The idea that structural white matter connectivity constrains functional connectivity (interactions among brain regions) has widely been explored in studies of brain networks; studies have mostly focused on the "average" strength of functional connectivity. The question of how structural connectivity constrains the "variability" of functional connectivity remains unresolved. In this study, we investigated the variability of resting state functional connectivity that was acquired every 3 h within a single day from 12 participants (eight time sessions within a 24-h period, 165 scans per session). Three different types of functional connectivity (functional connectivity based on Pearson correlation, direct functional connectivity based on partial correlation, and the pseudo functional connectivity produced by their difference) were estimated from resting state functional magnetic resonance imaging data along with structural connectivity defined using fiber tractography of diffusion tensor imaging. Those types of functional connectivity were evaluated with regard to properties of structural connectivity (fiber streamline counts and lengths) and types of structural connectivity such as intra-/inter-hemispheric edges and topological edge types in the rich club organization. We observed that the structural connectivity constrained the variability of direct functional connectivity more than pseudo-functional connectivity and that the constraints depended strongly on structural connectivity types. The structural constraints were greater for intra-hemispheric and heterologous inter-hemispheric edges than homologous inter-hemispheric edges, and feeder and local edges than rich club edges in the rich club architecture. While each edge was highly variable, the multivariate patterns of edge involvement, especially the direct functional connectivity patterns among the rich club brain regions, showed low variability over time. This study suggests that structural connectivity not only constrains the strength of functional connectivity, but also the within-a-day variability of functional connectivity and connectivity patterns, particularly the direct functional connectivity among brain regions.

Comparative evaluation of the white matter fiber integrity in patients with prelingual and postlingual deafness.

The aim of this study was to perform a comparative investigation of white matter integrity in patients with prelingual and postlingual deafness; we carried out a tract-based statistical analysis of diffusion tensor anisotropy in eight and ten adults with prelingual and postlingual deafness, respectively. Patients with deafness showed significant decreases in diffusion anisotropy at the right internal capsule, the right thalamus, and the splenium of the corpus callosum as well as within the bilateral superior temporal gyrus (including Heschl gyrus) and right temporal white matter. Furthermore, relative to patients with postlingual deafness, those with prelingual deafness showed lower anisotropy in the right superior temporal gyrus, bilateral temporal white matter, and the genu and anterior body of the corpus callosum. We believe that, in patients with deafness, reception of early auditory stimuli before language acquisition might be more critical to white matter maturation and brain reorganization than the nature of auditory stimuli itself or the duration of disuse. These findings provide the theoretical background for early auditory rehabilitation.