BVAR-Connect: A Variational Bayes Approach to Multi-Subject Vector Autoregressive Models for Inference on Brain Connectivity Networks

Bayesian varying-effects vector autoregressive models for inference of brain connectivity networks and covariate effects in pediatric traumatic brain injury

In this article, we develop an analytical approach for estimating brain connectivity networks that accounts for subject heterogeneity. More specifically, we consider a novel extension of a multi-subject Bayesian vector autoregressive model that estimates group-specific directed brain connectivity networks and accounts for the effects of covariates on the network edges. We adopt a flexible approach, allowing for (possibly) nonlinear effects of the covariates on edge strength via a novel Bayesian nonparametric prior that employs a weighted mixture of Gaussian processes. For posterior inference, we achieve computational scalability by implementing a variational Bayes scheme. Our approach enables simultaneous estimation of group-specific networks and selection of relevant covariate effects. We show improved performance over competing two-stage approaches on simulated data. We apply our method on resting-state functional magnetic resonance imaging data from children with a history of traumatic brain injury (TBI) and healthy controls to estimate the effects of age and sex on the group-level connectivities. Our results highlight differences in the distribution of parent nodes. They also suggest alteration in the relation of age, with peak edge strength in children with TBI, and differences in effective connectivity strength between males and females.

Effective connectivity between resting-state networks in depression

RATIONALE

Although depression has been widely researched, findings characterizing how brain regions influence each other remains scarce, yet this is critical for research on antidepressant treatments and individual responses to particular treatments.

OBJECTIVES

To identify pre-treatment resting state effective connectivity (rsEC) patterns in patients with major depressive disorder (MDD) and explore their relationship with treatment response.

METHODS

Thirty-four drug-free MDD patients had an MRI scan and were subsequently treated for 6 weeks with an SSRI escitalopram 10 mg daily; the response was defined as ≥50% decrease in Hamilton Depression Rating Scale (HAMD) score.

RESULTS

rsEC networks in default mode, central executive, and salience networks were identified for patients with depression. Exploratory analyses indicated higher connectivity strength related to baseline depression severity and response to treatment.

CONCLUSIONS

Preliminary analyses revealed widespread dysfunction of rsEC in depression. Functional rsEC may be useful as a predictive tool for antidepressant treatment response. A primary limitation of the current study was the small size; however, the group was carefully chosen, well-characterized, and included only medication-free patients. Further research in large samples of placebo-controlled studies would be required to confirm the results.

Recommendations and publication guidelines for studies using frequency domain and time-frequency domain analyses of neural time series

Since its beginnings in the early 20th century, the psychophysiological study of human brain function has included research into the spectral properties of electrical and magnetic brain signals. Now, dramatic advances in digital signal processing, biophysics, and computer science have enabled increasingly sophisticated methodology for neural time series analysis. Innovations in hardware and recording techniques have further expanded the range of tools available to researchers interested in measuring, quantifying, modeling, and altering the spectral properties of neural time series. These tools are increasingly used in the field, by a growing number of researchers who vary in their training, background, and research interests. Implementation and reporting standards also vary greatly in the published literature, causing challenges for authors, readers, reviewers, and editors alike. The present report addresses this issue by providing recommendations for the use of these methods, with a focus on foundational aspects of frequency domain and time-frequency analyses. It also provides publication guidelines, which aim to (1) foster replication and scientific rigor, (2) assist new researchers who wish to enter the field of brain oscillations, and (3) facilitate communication among authors, reviewers, and editors.

Effective connectivity in the default mode network after paediatric traumatic brain injury

Children who experience a traumatic brain injury (TBI) are at elevated risk for a range of negative cognitive and neuropsychological outcomes. Identifying which children are at greatest risk for negative outcomes can be difficult due to the heterogeneity of TBI. To address this barrier, the current study applied a novel method of characterizing brain connectivity networks, Bayesian multi-subject vector autoregressive modelling (BVAR-connect), which used white matter integrity as priors to evaluate effective connectivity-the time-dependent relationship in functional magnetic resonance imaging (fMRI) activity between two brain regions-within the default mode network (DMN). In a prospective longitudinal study, children ages 8-15 years with mild to severe TBI underwent diffusion tensor imaging and resting state fMRI 7 weeks after injury; post-concussion and anxiety symptoms were assessed 7 months after injury. The goals of this study were to (1) characterize differences in positive effective connectivity of resting-state DMN circuitry between healthy controls and children with TBI, (2) determine if severity of TBI was associated with differences in DMN connectivity and (3) evaluate whether patterns of DMN effective connectivity predicted persistent post-concussion symptoms and anxiety. Healthy controls had unique positive connectivity that mostly emerged from the inferior temporal lobes. In contrast, children with TBI had unique effective connectivity among orbitofrontal and parietal regions. These positive orbitofrontal-parietal DMN effective connectivity patterns also differed by TBI severity and were associated with persisting behavioural outcomes. Effective connectivity may be a sensitive neuroimaging marker of TBI severity as well as a predictor of chronic post-concussion symptoms and anxiety.

Latent Network Estimation and Variable Selection for Compositional Data Via Variational EM

Network estimation and variable selection have been extensively studied in the statistical literature, but only recently have those two challenges been addressed simultaneously. In this article, we seek to develop a novel method to simultaneously estimate network interactions and associations to relevant covariates for count data, and specifically for compositional data, which have a fixed sum constraint. We use a hierarchical Bayesian model with latent layers and employ spike-and-slab priors for both edge and covariate selection. For posterior inference, we develop a novel variational inference scheme with an expectation-maximization step, to enable efficient estimation. Through simulation studies, we demonstrate that the proposed model outperforms existing methods in its accuracy of network recovery. We show the practical utility of our model via an application to microbiome data. The human microbiome has been shown to contribute too many of the functions of the human body, and also to be linked with a number of diseases. In our application, we seek to better understand the interaction between microbes and relevant covariates, as well as the interaction of microbes with each other. We call our algorithm simultaneous inference for networks and covariates and provide a Python implementation, which is available online.