A multidimensional investigation of sleep and biopsychosocialprofiles with associated neural signatures.

Sleep is essential for optimal functioning and health. Interconnected to multiple biological, psychological and socio-environmental factors (i.e., biopsychosocial factors), the multidimensional nature of sleep is rarely capitalized on in research. Here, we deployed a data-driven approach to identify sleep-biopsychosocial profiles that linked self-reported sleep patterns to inter-individual variability in health, cognition, and lifestyle factors in 770 healthy young adults. We uncovered five profiles, including two profiles reflecting general psychopathology associated with either reports of general poor sleep or an absence of sleep complaints (i.e., sleep resilience) respectively. The three other profiles were driven by sedative-hypnotics-use and social satisfaction, sleep duration and cognitive performance, and sleep disturbance linked to cognition and mental health. Furthermore, identified sleep-biopsychosocial profiles displayed unique patterns of brain network organization. In particular, somatomotor network connectivity alterations were involved in the relationships between sleep and biopsychosocial factors. These profiles can potentially untangle the interplay between individuals' variability in sleep, health, cognition and lifestyle - equipping research and clinical settings to better support individual's well-being.

DeepResBat: deep residual batch harmonization accounting for covariate distribution differences.

Pooling MRI data from multiple datasets requires harmonization to reduce undesired inter-site variabilities, while preserving effects of biological variables (or covariates). The popular harmonization approach ComBat uses a mixed effect regression framework that explicitly accounts for covariate distribution differences across datasets. There is also significant interest in developing harmonization approaches based on deep neural networks (DNNs), such as conditional variational autoencoder (cVAE). However, current DNN approaches do not explicitly account for covariate distribution differences across datasets. Here, we provide mathematical results, suggesting that not accounting for covariates can lead to suboptimal harmonization outcomes. We propose two DNN-based harmonization approaches that explicitly account for covariate distribution differences across datasets: covariate VAE (coVAE) and DeepResBat. The coVAE approach is a natural extension of cVAE by concatenating covariates and site information with site- and covariate-invariant latent representations. DeepResBat adopts a residual framework inspired by ComBat. DeepResBat first removes the effects of covariates with nonlinear regression trees, followed by eliminating site differences with cVAE. Finally, covariate effects are added back to the harmonized residuals. Using three datasets from three different continents with a total of 2787 participants and 10085 anatomical T1 scans, we find that DeepResBat and coVAE outperformed ComBat, CovBat and cVAE in terms of removing dataset differences, while enhancing biological effects of interest. However, coVAE hallucinates spurious associations between anatomical MRI and covariates even when no association exists. Therefore, future studies proposing DNN-based harmonization approaches should be aware of this false positive pitfall. Overall, our results suggest that DeepResBat is an effective deep learning alternative to ComBat.

Comparison between gradients and parcellations for functional connectivity prediction of behavior.

Resting-state functional connectivity (RSFC) is widely used to predict behavioral measures. To predict behavioral measures, representing RSFC with parcellations and gradients are the two most popular approaches. Here, we compare parcellation and gradient approaches for RSFC-based prediction of a broad range of behavioral measures in the Human Connectome Project (HCP) and Adolescent Brain Cognitive Development (ABCD) datasets. Among the parcellation approaches, we consider group-average "hard" parcellations (Schaefer et al., 2018), individual-specific "hard" parcellations (Kong et al., 2021a), and an individual-specific "soft" parcellation (spatial independent component analysis with dual regression; Beckmann et al., 2009). For gradient approaches, we consider the well-known principal gradients (Margulies et al., 2016) and the local gradient approach that detects local RSFC changes (Laumann et al., 2015). Across two regression algorithms, individual-specific hard-parcellation performs the best in the HCP dataset, while the principal gradients, spatial independent component analysis and group-average "hard" parcellations exhibit similar performance. On the other hand, principal gradients and all parcellation approaches perform similarly in the ABCD dataset. Across both datasets, local gradients perform the worst. Finally, we find that the principal gradient approach requires at least 40 to 60 gradients to perform as well as parcellation approaches. While most principal gradient studies utilize a single gradient, our results suggest that incorporating higher order gradients can provide significant behaviorally relevant information. Future work will consider the inclusion of additional parcellation and gradient approaches for comparison.

A parsimonious description of global functional brain organization in three spatiotemporal patterns.

Resting-state functional magnetic resonance imaging (MRI) has yielded seemingly disparate insights into large-scale organization of the human brain. The brain's large-scale organization can be divided into two broad categories: zero-lag representations of functional connectivity structure and time-lag representations of traveling wave or propagation structure. In this study, we sought to unify observed phenomena across these two categories in the form of three low-frequency spatiotemporal patterns composed of a mixture of standing and traveling wave dynamics. We showed that a range of empirical phenomena, including functional connectivity gradients, the task-positive/task-negative anti-correlation pattern, the global signal, time-lag propagation patterns, the quasiperiodic pattern and the functional connectome network structure, are manifestations of these three spatiotemporal patterns. These patterns account for much of the global spatial structure that underlies functional connectivity analyses and unifies phenomena in resting-state functional MRI previously thought distinct.

Control theory illustrates the energy efficiency in the dynamic reconfiguration of functional connectivity.

The brain's functional connectivity fluctuates over time instead of remaining steady in a stationary mode even during the resting state. This fluctuation establishes the dynamical functional connectivity that transitions in a non-random order between multiple modes. Yet it remains unexplored how the transition facilitates the entire brain network as a dynamical system and what utility this mechanism for dynamic reconfiguration can bring over the widely used graph theoretical measurements. To address these questions, we propose to conduct an energetic analysis of functional brain networks using resting-state fMRI and behavioral measurements from the Human Connectome Project. Through comparing the state transition energy under distinct adjacent matrices, we justify that dynamic functional connectivity leads to 60% less energy cost to support the resting state dynamics than static connectivity when driving the transition through default mode network. Moreover, we demonstrate that combining graph theoretical measurements and our energy-based control measurements as the feature vector can provide complementary prediction power for the behavioral scores (Combination vs. Control: t = 9.41, p = 1.64e-13; Combination vs. Graph: t = 4.92, p = 3.81e-6). Our approach integrates statistical inference and dynamical system inspection towards understanding brain networks.

Identification of community structure-based brain states and transitions using functional MRI.

Brain function relies on a precisely coordinated and dynamic balance between the functional integration and segregation of distinct networks. Characterizing the way in which brain regions reconfigure their interactions to give rise to distinct but hidden brain states remains an open challenge. In this paper, we propose a Bayesian method for characterizing community structure-based latent brain states and showcase a novel strategy based on posterior predictive discrepancy using the latent block model to detect transitions between community structures in blood oxygen level-dependent (BOLD) time series. The set of estimated parameters in the model includes a latent label vector that assigns network nodes to communities, and also block model parameters that reflect the weighted connectivity within and between communities. Besides extensive in-silico model evaluation, we also provide empirical validation (and replication) using the Human Connectome Project (HCP) dataset of 100 healthy adults. Our results obtained through an analysis of task-fMRI data during working memory performance show appropriate lags between external task demands and change-points between brain states, with distinctive community patterns distinguishing fixation, low-demand and high-demand task conditions.

Macroscale and microcircuit dissociation of focal and generalized human epilepsies.

Thalamo-cortical pathology plays key roles in both generalized and focal epilepsies, but there is little work directly comparing these syndromes at the level of whole-brain mechanisms. Using multimodal imaging, connectomics, and computational simulations, we examined thalamo-cortical and cortico-cortical signatures and underlying microcircuits in 96 genetic generalized (GE) and 107 temporal lobe epilepsy (TLE) patients, along with 65 healthy controls. Structural and functional network profiling highlighted extensive atrophy, microstructural disruptions and decreased thalamo-cortical connectivity in TLE, while GE showed only subtle structural anomalies paralleled by enhanced thalamo-cortical connectivity. Connectome-informed biophysical simulations indicated modest increases in subcortical drive contributing to cortical dynamics in GE, while TLE presented with reduced subcortical drive and imbalanced excitation-inhibition within limbic and somatomotor microcircuits. Multiple sensitivity analyses supported robustness. Our multiscale analyses differentiate human focal and generalized epilepsy at the systems-level, showing paradoxically more severe microcircuit and macroscale imbalances in the former.

Charting brain growth in tandem with brain templates at school age.

Brain growth charts and age-normed brain templates are essential resources for researchers to eventually contribute to the care of individuals with atypical developmental trajectories. The present work generates age-normed brain templates for children and adolescents at one-year intervals and the corresponding growth charts to investigate the influences of age and ethnicity using a common pediatric neuroimaging protocol. Two accelerated longitudinal cohorts with the identical experimental design were implemented in the United States and China. Anatomical magnetic resonance imaging (MRI) of typically developing school-age children (TDC) was obtained up to three times at nominal intervals of 1.25 years. The protocol generated and compared population- and age-specific brain templates and growth charts, respectively. A total of 674 Chinese pediatric MRI scans were obtained from 457 Chinese TDC and 190 American pediatric MRI scans were obtained from 133 American TDC. Population- and age-specific brain templates were used to quantify warp cost, the differences between individual brains and brain templates. Volumetric growth charts for labeled brain network areas were generated. Shape analyses of cost functions supported the necessity of age-specific and ethnicity-matched brain templates, which was confirmed by growth chart analyses. These analyses revealed volumetric growth differences between the two ethnicities primarily in lateral frontal and parietal areas, regions which are most variable across individuals in regard to their structure and function. Age- and ethnicity-specific brain templates facilitate establishing unbiased pediatric brain growth charts, indicating the necessity of the brain charts and brain templates generated in tandem. These templates and growth charts as well as related codes have been made freely available to the public for open neuroscience (https://github.com/zuoxinian/CCS/tree/master/H3/GrowthCharts).

Individual-specific fMRI-Subspaces improve functional connectivity prediction of behavior.

There is significant interest in using resting-state functional connectivity (RSFC) to predict human behavior. Good behavioral prediction should in theory require RSFC to be sufficiently distinct across participants; if RSFC were the same across participants, then behavioral prediction would obviously be poor. Therefore, we hypothesize that removing common resting-state functional magnetic resonance imaging (rs-fMRI) signals that are shared across participants would improve behavioral prediction. Here, we considered 803 participants from the human connectome project (HCP) with four rs-fMRI runs. We applied the common and orthogonal basis extraction (COBE) technique to decompose each HCP run into two subspaces: a common (group-level) subspace shared across all participants and a subject-specific subspace. We found that the first common COBE component of the first HCP run was localized to the visual cortex and was unique to the run. On the other hand, the second common COBE component of the first HCP run and the first common COBE component of the remaining HCP runs were highly similar and localized to regions within the default network, including the posterior cingulate cortex and precuneus. Overall, this suggests the presence of run-specific (state-specific) effects that were shared across participants. By removing the first and second common COBE components from the first HCP run, and the first common COBE component from the remaining HCP runs, the resulting RSFC improves behavioral prediction by an average of 11.7% across 58 behavioral measures spanning cognition, emotion and personality.

Inversion of a large-scale circuit model reveals a cortical hierarchy in the dynamic resting human brain.

We considered a large-scale dynamical circuit model of human cerebral cortex with region-specific microscale properties. The model was inverted using a stochastic optimization approach, yielding markedly better fit to new, out-of-sample resting functional magnetic resonance imaging (fMRI) data. Without assuming the existence of a hierarchy, the estimated model parameters revealed a large-scale cortical gradient. At one end, sensorimotor regions had strong recurrent connections and excitatory subcortical inputs, consistent with localized processing of external stimuli. At the opposing end, default network regions had weak recurrent connections and excitatory subcortical inputs, consistent with their role in internal thought. Furthermore, recurrent connection strength and subcortical inputs provided complementary information for differentiating the levels of the hierarchy, with only the former showing strong associations with other macroscale and microscale proxies of cortical hierarchies (meta-analysis of cognitive functions, principal resting fMRI gradient, myelin, and laminar-specific neuronal density). Overall, this study provides microscale insights into a macroscale cortical hierarchy in the dynamic resting brain.

Co-activated yet disconnected-Neural correlates of eye closures when trying to stay awake.

Spontaneous eye-closures that herald sleep onset become more frequent when we are sleep deprived. Although these are typically associated with decreased responsiveness to external stimuli, it is less clear what occurs in the brain at these transitions to drowsiness and light sleep. To investigate this, task-free fMRI of sleep-deprived participants was acquired. BOLD activity associated with periods of spontaneously occurring eye closures were marked and analyzed. We observed concurrent and extensive hypnagogic co-activation of the extrastriate visual, auditory, and somatosensory cortices as well as the default mode network, consistent with internal sensory activity without external stimulation. Co-activation of fronto-parietal areas known to mediate attentional control could correspond with participants resisting sleep or additional engagement of mental imagery. This constellation of signal changes differed from those elicited by cued eye closures of similar duration and distribution in the same, rested participants. They also differ from signal changes associated with mind-wandering and consolidated light sleep. Concurrent with the observed event-related changes, eye closures elicited additional reduction in functional connectivity within nodes of the DMN and DAN, superposed on already reduced connectivity associated with sleep deprivation. There was concurrent deactivation of the thalamus during eye-closure during the sleep-deprived state but almost similar changes occurred in the well-rested state that may also be relevant. These findings highlight the dynamic shifts in brain activity and connectivity at border between wakefulness and sleep.