Decomposing Neural Circuit Function into Information Processing Primitives.

It is challenging to measure how specific aspects of coordinated neural dynamics translate into operations of information processing and, ultimately, cognitive functions. An obstacle is that simple circuit mechanisms-such as self-sustained or propagating activity and nonlinear summation of inputs-do not directly give rise to high-level functions. Nevertheless, they already implement simple the information carried by neural activity. Here, we propose that distinct functions, such as stimulus representation, working memory, or selective attention, stem from different combinations and types of low-level manipulations of information or information processing primitives. To test this hypothesis, we combine approaches from information theory with simulations of multi-scale neural circuits involving interacting brain regions that emulate well-defined cognitive functions. Specifically, we track the information dynamics emergent from patterns of neural dynamics, using quantitative metrics to detect where and when information is actively buffered, transferred or nonlinearly merged, as possible modes of low-level processing (storage, transfer and modification). We find that neuronal subsets maintaining representations in working memory or performing attentional gain modulation are signaled by their boosted involvement in operations of information storage or modification, respectively. Thus, information dynamic metrics, beyond detecting which network units participate in cognitive processing, also promise to specify how and when they do it, that is, through which type of primitive computation, a capability that may be exploited for the analysis of experimental recordings.

Different Faces of Medial Beta-Band Activity Reflect Distinct Visuomotor Feedback Signals.

Beta-band (13-35 Hz) modulations following reward, task outcome feedback, and error have been described in cognitive and/or motor adaptation tasks. Observations from different studies are, however, difficult to conciliate. Among the studies that used cognitive response selection tasks, several reported an increase in beta-band activity following reward, whereas others observed increased beta power after negative feedback. Moreover, in motor adaptation tasks, an attenuation of the postmovement beta rebound follows a movement execution error induced by visual or mechanical perturbations. Given that kinematic error typically leads to negative task-outcome feedback (e.g., target missed), one may wonder how contradictory modulations, beta power decrease with movement error versus beta power increase with negative feedback, may coexist. We designed a motor adaptation task in which female and male participants experience varied feedbacks-binary success/failure feedback, kinematic error, and sensory-prediction error-and demonstrate that beta-band modulations in opposite directions coexist at different spatial locations, time windows, and frequency ranges. First, high beta power in the medial frontal cortex showed opposite modulations well separated in time when compared in success and failure trials; that is, power was higher in success trials just after the binary success feedback, whereas it was lower in the postmovement period compared with failure trials. Second, although medial frontal high-beta activity was sensitive to task outcome, low-beta power in the medial parietal cortex was strongly attenuated following movement execution error but was not affected by either the outcome of the task or sensory-prediction error. These findings suggest that medial beta activity in different spatio-temporal-spectral configurations play a multifaceted role in encoding qualitatively distinct feedback signals.SIGNIFICANCE STATEMENT Beta-band activity reflects neural processes well beyond sensorimotor functions, including cognition and motivation. By disentangling alternative spatio-temporal-spectral patterns of possible beta-oscillatory activity, we reconcile a seemingly discrepant literature. First, high-beta power in the medial frontal cortex showed opposite modulations separated in time in success and failure trials; power was higher in success trials just after success feedback and lower in the postmovement period compared with failure trials. Second, although medial frontal high-beta activity was sensitive to task outcome, low-beta power in the medial parietal cortex was strongly attenuated following movement execution error but was not affected by the task outcome or the sensory-prediction error. We propose that medial beta activity reflects distinct feedback signals depending on its anatomic location, time window, and frequency range.

Perturbed Information Processing Complexity in Experimental Epilepsy.

Comorbidities, such as cognitive deficits, which often accompany epilepsies, constitute a basal state, while seizures are rare and transient events. This suggests that neural dynamics, in particular those supporting cognitive function, are altered in a permanent manner in epilepsy. Here, we test the hypothesis that primitive processes of information processing at the core of cognitive function (i.e., storage and sharing of information) are altered in the hippocampus and the entorhinal cortex in experimental epilepsy in adult, male Wistar rats. We find that information storage and sharing are organized into substates across the stereotypic states of slow and theta oscillations in both epilepsy and control conditions. However, their internal composition and organization through time are disrupted in epilepsy, partially losing brain state selectivity compared with controls, and shifting toward a regimen of disorder. We propose that the alteration of information processing at this algorithmic level of computation, the theoretical intermediate level between structure and function, may be a mechanism behind the emergent and widespread comorbidities associated with epilepsy, and perhaps other disorders.SIGNIFICANCE STATEMENT Comorbidities, such as cognitive deficits, which often accompany epilepsies, constitute a basal state, while seizures are rare and transient events. This suggests that neural dynamics, in particular those supporting cognitive function, are altered in a permanent manner in epilepsy. Here, we show that basic processes of information processing at the core of cognitive function (i.e., storage and sharing of information) are altered in the hippocampus and the entorhinal cortex (two regions involved in memory processes) in experimental epilepsy. Such disruption of information processing at the algorithmic level itself could underlie the general performance impairments in epilepsy.

Event-related variability is modulated by task and development.

In carefully designed experimental paradigms, cognitive scientists interpret the mean event-related potentials (ERP) in terms of cognitive operations. However, the huge signal variability from one trial to the next, questions the representability of such mean events. We explored here whether this variability is an unwanted noise, or an informative part of the neural response. We took advantage of the rapid changes in the visual system during human infancy and analyzed the variability of visual responses to central and lateralized faces in 2-to 6-month-old infants compared to adults using high-density electroencephalography (EEG). We observed that neural trajectories of individual trials always remain very far from ERP components, only moderately bending their direction with a substantial temporal jitter across trials. However, single trial trajectories displayed characteristic patterns of acceleration and deceleration when approaching ERP components, as if they were under the active influence of steering forces causing transient attraction and stabilization. These dynamic events could only partly be accounted for by induced microstate transitions or phase reset phenomena. Importantly, these structured modulations of response variability, both between and within trials, had a rich sequential organization, which in infants, was modulated by the task difficulty and age. Our approaches to characterize Event Related Variability (ERV) expand on classic ERP analyses and provide the first evidence for the functional role of ongoing neural variability in human infants.

Repairing Artifacts in Neural Activity Recordings Using Low-Rank Matrix Estimation.

Electrophysiology recordings are frequently affected by artifacts (e.g., subject motion or eye movements), which reduces the number of available trials and affects the statistical power. When artifacts are unavoidable and data are scarce, signal reconstruction algorithms that allow for the retention of sufficient trials become crucial. Here, we present one such algorithm that makes use of large spatiotemporal correlations in neural signals and solves the low-rank matrix completion problem, to fix artifactual entries. The method uses a gradient descent algorithm in lower dimensions to learn the missing entries and provide faithful reconstruction of signals. We carried out numerical simulations to benchmark the method and estimate optimal hyperparameters for actual EEG data. The fidelity of reconstruction was assessed by detecting event-related potentials (ERP) from a highly artifacted EEG time series from human infants. The proposed method significantly improved the standardized error of the mean in ERP group analysis and a between-trial variability analysis compared to a state-of-the-art interpolation technique. This improvement increased the statistical power and revealed significant effects that would have been deemed insignificant without reconstruction. The method can be applied to any time-continuous neural signal where artifacts are sparse and spread out across epochs and channels, increasing data retention and statistical power.

Cortical ignition dynamics is tightly linked to the core organisation of the human connectome.

The capability of cortical regions to flexibly sustain an "ignited" state of activity has been discussed in relation to conscious perception or hierarchical information processing. Here, we investigate how the intrinsic propensity of different regions to get ignited is determined by the specific topological organisation of the structural connectome. More specifically, we simulated the resting-state dynamics of mean-field whole-brain models and assessed how dynamic multistability and ignition differ between a reference model embedding a realistic human connectome, and alternative models based on a variety of randomised connectome ensembles. We found that the strength of global excitation needed to first trigger ignition in a subset of regions is substantially smaller for the model embedding the empirical human connectome. Furthermore, when increasing the strength of excitation, the propagation of ignition outside of this initial core-which is able to self-sustain its high activity-is way more gradual than for any of the randomised connectomes, allowing for graded control of the number of ignited regions. We explain both these assets in terms of the exceptional weighted core-shell organisation of the empirical connectome, speculating that this topology of human structural connectivity may be attuned to support enhanced ignition dynamics.

Relations between large-scale brain connectivity and effects of regional stimulation depend on collective dynamical state.

At the macroscale, the brain operates as a network of interconnected neuronal populations, which display coordinated rhythmic dynamics that support interareal communication. Understanding how stimulation of different brain areas impacts such activity is important for gaining basic insights into brain function and for further developing therapeutic neurmodulation. However, the complexity of brain structure and dynamics hinders predictions regarding the downstream effects of focal stimulation. More specifically, little is known about how the collective oscillatory regime of brain network activity-in concert with network structure-affects the outcomes of perturbations. Here, we combine human connectome data and biophysical modeling to begin filling these gaps. By tuning parameters that control collective system dynamics, we identify distinct states of simulated brain activity and investigate how the distributed effects of stimulation manifest at different dynamical working points. When baseline oscillations are weak, the stimulated area exhibits enhanced power and frequency, and due to network interactions, activity in this excited frequency band propagates to nearby regions. Notably, beyond these linear effects, we further find that focal stimulation causes more distributed modifications to interareal coherence in a band containing regions' baseline oscillation frequencies. Importantly, depending on the dynamical state of the system, these broadband effects can be better predicted by functional rather than structural connectivity, emphasizing a complex interplay between anatomical organization, dynamics, and response to perturbation. In contrast, when the network operates in a regime of strong regional oscillations, stimulation causes only slight shifts in power and frequency, and structural connectivity becomes most predictive of stimulation-induced changes in network activity patterns. In sum, this work builds upon and extends previous computational studies investigating the impacts of stimulation, and underscores the fact that both the stimulation site, and, crucially, the regime of brain network dynamics, can influence the network-wide responses to local perturbations.

Dynamic Functional Connectivity between order and randomness and its evolution across the human adult lifespan.

Functional Connectivity (FC) during resting-state or task conditions is not static but inherently dynamic. Yet, there is no consensus on whether fluctuations in FC may resemble isolated transitions between discrete FC states rather than continuous changes. This quarrel hampers advancing the study of dynamic FC. This is unfortunate as the structure of fluctuations in FC can certainly provide more information about developmental changes, aging, and progression of pathologies. We merge the two perspectives and consider dynamic FC as an ongoing network reconfiguration, including a stochastic exploration of the space of possible steady FC states. The statistical properties of this random walk deviate both from a purely "order-driven" dynamics, in which the mean FC is preserved, and from a purely "randomness-driven" scenario, in which fluctuations of FC remain uncorrelated over time. Instead, dynamic FC has a complex structure endowed with long-range sequential correlations that give rise to transient slowing and acceleration epochs in the continuous flow of reconfiguration. Our analysis for fMRI data in healthy elderly revealed that dynamic FC tends to slow down and becomes less complex as well as more random with increasing age. These effects appear to be strongly associated with age-related changes in behavioural and cognitive performance.

Modular slowing of resting-state dynamic functional connectivity as a marker of cognitive dysfunction induced by sleep deprivation.

Dynamic Functional Connectivity (dFC) in the resting state (rs) is considered as a correlate of cognitive processing. Describing dFC as a flow across morphing connectivity configurations, our notion of dFC speed quantifies the rate at which FC networks evolve in time. Here we probe the hypothesis that variations of rs dFC speed and cognitive performance are selectively interrelated within specific functional subnetworks. In particular, we focus on Sleep Deprivation (SD) as a reversible model of cognitive dysfunction. We found that whole-brain level (global) dFC speed significantly slows down after 24h of SD. However, the reduction in global dFC speed does not correlate with variations of cognitive performance in individual tasks, which are subtle and highly heterogeneous. On the contrary, we found strong correlations between performance variations in individual tasks -including Rapid Visual Processing (RVP, assessing sustained visual attention)- and dFC speed quantified at the level of functional sub-networks of interest. Providing a compromise between classic static FC (no time) and global dFC (no space), modular dFC speed analyses allow quantifying a different speed of dFC reconfiguration independently for sub-networks overseeing different tasks. Importantly, we found that RVP performance robustly correlates with the modular dFC speed of a characteristic frontoparietal module.

Functional connectivity dynamics: modeling the switching behavior of the resting state.

Functional connectivity (FC) sheds light on the interactions between different brain regions. Besides basic research, it is clinically relevant for applications in Alzheimer's disease, schizophrenia, presurgical planning, epilepsy, and traumatic brain injury. Simulations of whole-brain mean-field computational models with realistic connectivity determined by tractography studies enable us to reproduce with accuracy aspects of average FC in the resting state. Most computational studies, however, did not address the prominent non-stationarity in resting state FC, which may result in large intra- and inter-subject variability and thus preclude an accurate individual predictability. Here we show that this non-stationarity reveals a rich structure, characterized by rapid transitions switching between a few discrete FC states. We also show that computational models optimized to fit time-averaged FC do not reproduce these spontaneous state transitions and, thus, are not qualitatively superior to simplified linear stochastic models, which account for the effects of structure alone. We then demonstrate that a slight enhancement of the non-linearity of the network nodes is sufficient to broaden the repertoire of possible network behaviors, leading to modes of fluctuations, reminiscent of some of the most frequently observed Resting State Networks. Because of the noise-driven exploration of this repertoire, the dynamics of FC qualitatively change now and display non-stationary switching similar to empirical resting state recordings (Functional Connectivity Dynamics (FCD)). Thus FCD bear promise to serve as a better biomarker of resting state neural activity and of its pathologic alterations.

Gamma oscillatory complexity conveys behavioral information in hippocampal networks.

The hippocampus and entorhinal cortex exhibit rich oscillatory patterns critical for cognitive functions. In the hippocampal region CA1, specific gamma-frequency oscillations, timed at different phases of the ongoing theta rhythm, are hypothesized to facilitate the integration of information from varied sources and contribute to distinct cognitive processes. Here, we show that gamma elements -a multidimensional characterization of transient gamma oscillatory episodes- occur at any frequency or phase relative to the ongoing theta rhythm across all CA1 layers in male mice. Despite their low power and stochastic-like nature, individual gamma elements still carry behavior-related information and computational modeling suggests that they reflect neuronal firing. Our findings challenge the idea of rigid gamma sub-bands, showing that behavior shapes ensembles of irregular gamma elements that evolve with learning and depend on hippocampal layers. Widespread gamma diversity, beyond randomness, may thus reflect complexity, likely functional but invisible to classic average-based analyses.

Dynamic effective connectivity of inter-areal brain circuits.

Anatomic connections between brain areas affect information flow between neuronal circuits and the synchronization of neuronal activity. However, such structural connectivity does not coincide with effective connectivity (or, more precisely, causal connectivity), related to the elusive question "Which areas cause the present activity of which others?". Effective connectivity is directed and depends flexibly on contexts and tasks. Here we show that dynamic effective connectivity can emerge from transitions in the collective organization of coherent neural activity. Integrating simulation and semi-analytic approaches, we study mesoscale network motifs of interacting cortical areas, modeled as large random networks of spiking neurons or as simple rate units. Through a causal analysis of time-series of model neural activity, we show that different dynamical states generated by a same structural connectivity motif correspond to distinct effective connectivity motifs. Such effective motifs can display a dominant directionality, due to spontaneous symmetry breaking and effective entrainment between local brain rhythms, although all connections in the considered structural motifs are reciprocal. We show then that transitions between effective connectivity configurations (like, for instance, reversal in the direction of inter-areal interactions) can be triggered reliably by brief perturbation inputs, properly timed with respect to an ongoing local oscillation, without the need for plastic synaptic changes. Finally, we analyze how the information encoded in spiking patterns of a local neuronal population is propagated across a fixed structural connectivity motif, demonstrating that changes in the active effective connectivity regulate both the efficiency and the directionality of information transfer. Previous studies stressed the role played by coherent oscillations in establishing efficient communication between distant areas. Going beyond these early proposals, we advance here that dynamic interactions between brain rhythms provide as well the basis for the self-organized control of this "communication-through-coherence", making thus possible a fast "on-demand" reconfiguration of global information routing modalities.

Model-free reconstruction of excitatory neuronal connectivity from calcium imaging signals.

A systematic assessment of global neural network connectivity through direct electrophysiological assays has remained technically infeasible, even in simpler systems like dissociated neuronal cultures. We introduce an improved algorithmic approach based on Transfer Entropy to reconstruct structural connectivity from network activity monitored through calcium imaging. We focus in this study on the inference of excitatory synaptic links. Based on information theory, our method requires no prior assumptions on the statistics of neuronal firing and neuronal connections. The performance of our algorithm is benchmarked on surrogate time series of calcium fluorescence generated by the simulated dynamics of a network with known ground-truth topology. We find that the functional network topology revealed by Transfer Entropy depends qualitatively on the time-dependent dynamic state of the network (bursting or non-bursting). Thus by conditioning with respect to the global mean activity, we improve the performance of our method. This allows us to focus the analysis to specific dynamical regimes of the network in which the inferred functional connectivity is shaped by monosynaptic excitatory connections, rather than by collective synchrony. Our method can discriminate between actual causal influences between neurons and spurious non-causal correlations due to light scattering artifacts, which inherently affect the quality of fluorescence imaging. Compared to other reconstruction strategies such as cross-correlation or Granger Causality methods, our method based on improved Transfer Entropy is remarkably more accurate. In particular, it provides a good estimation of the excitatory network clustering coefficient, allowing for discrimination between weakly and strongly clustered topologies. Finally, we demonstrate the applicability of our method to analyses of real recordings of in vitro disinhibited cortical cultures where we suggest that excitatory connections are characterized by an elevated level of clustering compared to a random graph (although not extreme) and can be markedly non-local.

State-switching and high-order spatiotemporal organization of dynamic functional connectivity are disrupted by Alzheimer's disease.

Spontaneous activity during the resting state, tracked by BOLD fMRI imaging, or shortly rsfMRI, gives rise to brain-wide dynamic patterns of interregional correlations, whose structured flexibility relates to cognitive performance. Here, we analyze resting-state dynamic functional connectivity (dFC) in a cohort of older adults, including amnesic mild cognitive impairment (aMCI, N = 34) and Alzheimer's disease (AD, N = 13) patients, as well as normal control (NC, N = 16) and cognitively "supernormal" controls (SNC, N = 10) subjects. Using complementary state-based and state-free approaches, we find that resting-state fluctuations of different functional links are not independent but are constrained by high-order correlations between triplets or quadruplets of functionally connected regions. When contrasting patients with healthy subjects, we find that dFC between cingulate and other limbic regions is increasingly bursty and intermittent when ranking the four groups from SNC to NC, aMCI and AD. Furthermore, regions affected at early stages of AD pathology are less involved in higher order interactions in patient than in control groups, while pairwise interactions are not significantly reduced. Our analyses thus suggest that the spatiotemporal complexity of dFC organization is precociously degraded in AD and provides a richer window into the underlying neurobiology than time-averaged FC connections.

Synchronous chaos and broad band gamma rhythm in a minimal multi-layer model of primary visual cortex.

Visually induced neuronal activity in V1 displays a marked gamma-band component which is modulated by stimulus properties. It has been argued that synchronized oscillations contribute to these gamma-band activity. However, analysis of Local Field Potentials (LFPs) across different experiments reveals considerable diversity in the degree of oscillatory behavior of this induced activity. Contrast-dependent power enhancements can indeed occur over a broad band in the gamma frequency range and spectral peaks may not arise at all. Furthermore, even when oscillations are observed, they undergo temporal decorrelation over very few cycles. This is not easily accounted for in previous network modeling of gamma oscillations. We argue here that interactions between cortical layers can be responsible for this fast decorrelation. We study a model of a V1 hypercolumn, embedding a simplified description of the multi-layered structure of the cortex. When the stimulus contrast is low, the induced activity is only weakly synchronous and the network resonates transiently without developing collective oscillations. When the contrast is high, on the other hand, the induced activity undergoes synchronous oscillations with an irregular spatiotemporal structure expressing a synchronous chaotic state. As a consequence the population activity undergoes fast temporal decorrelation, with concomitant rapid damping of the oscillations in LFPs autocorrelograms and peak broadening in LFPs power spectra. We show that the strength of the inter-layer coupling crucially affects this spatiotemporal structure. We predict that layer VI inactivation should induce global changes in the spectral properties of induced LFPs, reflecting their slower temporal decorrelation in the absence of inter-layer feedback. Finally, we argue that the mechanism underlying the emergence of synchronous chaos in our model is in fact very general. It stems from the fact that gamma oscillations induced by local delayed inhibition tend to develop chaos when coupled by sufficiently strong excitation.

How Many Gammas? Redefining Hippocampal Theta-Gamma Dynamic During Spatial Learning.

The hippocampal formation is one of the brain systems in which the functional roles of coordinated oscillations in information representation and communication are better studied. Within this circuit, neuronal oscillations are conceived as a mechanism to precisely coordinate upstream and downstream neuronal ensembles, underlying dynamic exchange of information. Within a global reference framework provided by theta (θ) oscillations, different gamma-frequency (γ) carriers would temporally segregate information originating from different sources, thereby allowing networks to disambiguate convergent inputs. Two γ sub-bands were thus defined according to their frequency (slow γ, 30-80 Hz; medium γ, 60-120 Hz) and differential power distribution across CA1 dendritic layers. According to this prevalent model, layer-specific γ oscillations in CA1 would reliably identify the temporal dynamics of afferent inputs and may therefore aid in identifying specific memory processes (encoding for medium γ vs. retrieval for slow γ). However, this influential view, derived from time-averages of either specific γ sub-bands or different projection methods, might not capture the complexity of CA1 θ-γ interactions. Recent studies investigating γ oscillations at the θ cycle timescale have revealed a more dynamic and diverse landscape of θ-γ motifs, with many θ cycles containing multiple γ bouts of various frequencies. To properly capture the hippocampal oscillatory complexity, we have argued in this review that we should consider the entirety of the data and its multidimensional complexity. This will call for a revision of the actual model and will require the use of new tools allowing the description of individual γ bouts in their full complexity.

Cerebellar connectivity maps embody individual adaptive behavior in mice.

The cerebellar cortex encodes sensorimotor adaptation during skilled locomotor behaviors, however the precise relationship between synaptic connectivity and behavior is unclear. We studied synaptic connectivity between granule cells (GCs) and Purkinje cells (PCs) in murine acute cerebellar slices using photostimulation of caged glutamate combined with patch-clamp in developing or after mice adapted to different locomotor contexts. By translating individual maps into graph network entities, we found that synaptic maps in juvenile animals undergo critical period characterized by dissolution of their structure followed by the re-establishment of a patchy functional organization in adults. Although, in adapted mice, subdivisions in anatomical microzones do not fully account for the observed spatial map organization in relation to behavior, we can discriminate locomotor contexts with high accuracy. We also demonstrate that the variability observed in connectivity maps directly accounts for motor behavior traits at the individual level. Our findings suggest that, beyond general motor contexts, GC-PC networks also encode internal models underlying individual-specific motor adaptation.

Early memory deficits and extensive brain network disorganization in the AppNL-F/MAPT double knock-in mouse model of familial Alzheimer's disease.

A critical challenge in current research on Alzheimer's disease (AD) is to clarify the relationship between network dysfunction and the emergence of subtle memory deficits in itspreclinical stage. The AppNL-F/MAPT double knock-in (dKI) model with humanized ß-amyloid peptide (Aß) and tau was used to investigate both memory and network dysfunctions at an early stage. Young male dKI mice (2 to 6 months) were tested in three tasks taxing different aspects of recognition memory affected in preclinical AD. An early deficit first appeared in the object-place association task at the age of 4 months, when increased levels of ß-CTF and Aß were detected in both the hippocampus and the medial temporal cortex, and tau pathology was found only in the medial temporal cortex. Object-place task-dependent c-Fos activation was then analyzed in 22 subregions across the medial prefrontal cortex, claustrum, retrosplenial cortex, and medial temporal lobe. Increased c-Fos activation was detected in the entorhinal cortex and the claustrum of dKI mice. During recall, network efficiency was reduced across cingulate regions with a major disruption of information flow through the retrosplenial cortex. Our findings suggest that early perirhinal-entorhinal pathology is associated with abnormal activity which may spread to downstream regions such as the claustrum, the medial prefrontal cortex and ultimately the key retrosplenial hub which relays information from frontal to temporal lobes. The similarity between our findings and those reported in preclinical stages of AD suggests that the AppNL-F/MAPT dKI model has a high potential for providing key insights into preclinical AD.

Virtual Connectomic Datasets in Alzheimer's Disease and Aging Using Whole-Brain Network Dynamics Modelling.

Large neuroimaging datasets, including information about structural connectivity (SC) and functional connectivity (FC), play an increasingly important role in clinical research, where they guide the design of algorithms for automated stratification, diagnosis or prediction. A major obstacle is, however, the problem of missing features [e.g., lack of concurrent DTI SC and resting-state functional magnetic resonance imaging (rsfMRI) FC measurements for many of the subjects]. We propose here to address the missing connectivity features problem by introducing strategies based on computational whole-brain network modeling. Using two datasets, the Alzheimer's Disease Neuroimaging Initiative (ADNI) dataset and a healthy aging dataset, for proof-of-concept, we demonstrate the feasibility of virtual data completion (i.e., inferring "virtual FC" from empirical SC or "virtual SC" from empirical FC), by using self-consistent simulations of linear and nonlinear brain network models. Furthermore, by performing machine learning classification (to separate age classes or control from patient subjects), we show that algorithms trained on virtual connectomes achieve discrimination performance comparable to when trained on actual empirical data; similarly, algorithms trained on virtual connectomes can be used to successfully classify novel empirical connectomes. Completion algorithms can be combined and reiterated to generate realistic surrogate connectivity matrices in arbitrarily large number, opening the way to the generation of virtual connectomic datasets with network connectivity information comparable to the one of the original data.

Pathologically reduced neural flexibility recovers during psychotherapy of OCD patients.

Flexibility is a key feature of psychological health, allowing the individual to dynamically adapt to changing environmental demands, which is impaired in many psychiatric disorders like obsessive-compulsive disorder (OCD). Adequately responding to varying demands requires the brain to switch between different patterns of neural activity, which are represented by different brain network configurations (functional connectivity patterns). Here, we operationalize neural flexibility as the dissimilarity between consecutive connectivity matrices of brain regions (jump length). In total, 132 fMRI scans were obtained from 17 patients that were scanned four to five times during inpatient psychotherapy, and from 17 controls that were scanned at comparable time intervals. Significant negative correlations were found between the jump lengths and the symptom severity scores of OCD, depression, anxiety, and stress, suggesting that high symptom severity corresponds to inflexible brain functioning. Further analyses revealed that impaired reconfiguration (pattern stability) of the brain seems to be more related to general psychiatric impairment rather than to specific symptoms, e.g., of OCD or depression. Importantly, the group × time interaction of a repeated measures ANOVA was significant, as well as the post-hoc paired t-tests of the patients (first vs. last scan). The results suggest that psychotherapy is able to significantly increase the neural flexibility of patients. We conclude that psychiatric symptoms like anxiety, stress, depression, and OCD are associated with an impaired adaptivity of the brain. In general, our results add to the growing evidence that dynamic functional connectivity captures meaningful properties of brain functioning.