Virtual brain twins: from basic neuroscience to clinical use.

Virtual brain twins are personalized, generative and adaptive brain models based on data from an individual's brain for scientific and clinical use. After a description of the key elements of virtual brain twins, we present the standard model for personalized whole-brain network models. The personalization is accomplished using a subject's brain imaging data by three means: (1) assemble cortical and subcortical areas in the subject-specific brain space; (2) directly map connectivity into the brain models, which can be generalized to other parameters; and (3) estimate relevant parameters through model inversion, typically using probabilistic machine learning. We present the use of personalized whole-brain network models in healthy ageing and five clinical diseases: epilepsy, Alzheimer's disease, multiple sclerosis, Parkinson's disease and psychiatric disorders. Specifically, we introduce spatial masks for relevant parameters and demonstrate their use based on the physiological and pathophysiological hypotheses. Finally, we pinpoint the key challenges and future directions.

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.

The virtual aging brain: Causal inference supports interhemispheric dedifferentiation in healthy aging.

The mechanisms of cognitive decline and its variability during healthy aging are not fully understood, but have been associated with reorganization of white matter tracts and functional brain networks. Here, we built a brain network modeling framework to infer the causal link between structural connectivity and functional architecture and the consequent cognitive decline in aging. By applying in-silico interhemispheric degradation of structural connectivity, we reproduced the process of functional dedifferentiation during aging. Thereby, we found the global modulation of brain dynamics by structural connectivity to increase with age, which was steeper in older adults with poor cognitive performance. We validated our causal hypothesis via a deep-learning Bayesian approach. Our results might be the first mechanistic demonstration of dedifferentiation during aging leading to cognitive decline.

Delta-alpha cross-frequency coupling for different brain regions.

Neural interactions occur on different levels and scales. It is of particular importance to understand how they are distributed among different neuroanatomical and physiological relevant brain regions. We investigated neural cross-frequency couplings between different brain regions according to the Desikan-Killiany brain parcellation. The adaptive dynamic Bayesian inference method was applied to EEG measurements of healthy resting subjects in order to reconstruct the coupling functions. It was found that even after averaging over all subjects, the mean coupling function showed a characteristic waveform, confirming the direct influence of the delta-phase on the alpha-phase dynamics in certain brain regions and that the shape of the coupling function changes for different regions. While the averaged coupling function within a region was of similar form, the region-averaged coupling function was averaged out, which implies that there is a common dependence within separate regions across the subjects. It was also found that for certain regions the influence of delta on alpha oscillations is more pronounced and that oscillations that influence other are more evenly distributed across brain regions than the influenced oscillations. When presenting the information on brain lobes, it was shown that the influence of delta emanating from the brain as a whole is greatest on the alpha oscillations of the cingulate frontal lobe, and at the same time the influence of delta from the cingulate parietal brain lobe is greatest on the alpha oscillations of the whole brain.

Characterization of regional differences in resting-state fMRI with a data-driven network model of brain dynamics.

Model-based data analysis of whole-brain dynamics links the observed data to model parameters in a network of neural masses. Recently, studies focused on the role of regional variance of model parameters. Such analyses however necessarily depend on the properties of preselected neural mass model. We introduce a method to infer from the functional data both the neural mass model representing the regional dynamics and the region- and subject-specific parameters while respecting the known network structure. We apply the method to human resting-state fMRI. We find that the underlying dynamics can be described as noisy fluctuations around a single fixed point. The method reliably discovers three regional parameters with clear and distinct role in the dynamics, one of which is strongly correlated with the first principal component of the gene expression spatial map. The present approach opens a novel way to the analysis of resting-state fMRI with possible applications for understanding the brain dynamics during aging or neurodegeneration.

White-matter degradation and dynamical compensation support age-related functional alterations in human brain.

Structural connectivity of the brain at different ages is analyzed using diffusion-weighted magnetic resonance imaging (MRI) data. The largest decrease of streamlines is found in frontal regions and for long inter-hemispheric links. The average length of the tracts also decreases, but the clustering is unaffected. From functional MRI we identify age-related changes of dynamic functional connectivity (dFC) and spatial covariation features of functional connectivity (FC) links captured by metaconnectivity. They indicate more stable dFC, but wider range and variance of MC, whereas static features of FC did not show any significant differences with age. We implement individual connectivity in whole-brain models and test several hypotheses for the mechanisms of operation among underlying neural system. We demonstrate that age-related functional fingerprints are only supported if the model accounts for: (i) compensation of the individual brains for the overall loss of structural connectivity and (ii) decrease of propagation velocity due to the loss of myelination. We also show that with these 2 conditions, it is sufficient to decompose the time-delays as bimodal distribution that only distinguishes between intra- and inter-hemispheric delays, and that the same working point also captures the static FC the best, and produces the largest variability at slow time-scales.

Brain simulation as a cloud service: The Virtual Brain on EBRAINS.

The Virtual Brain (TVB) is now available as open-source services on the cloud research platform EBRAINS (ebrains.eu). It offers software for constructing, simulating and analysing brain network models including the TVB simulator; magnetic resonance imaging (MRI) processing pipelines to extract structural and functional brain networks; combined simulation of large-scale brain networks with small-scale spiking networks; automatic conversion of user-specified model equations into fast simulation code; simulation-ready brain models of patients and healthy volunteers; Bayesian parameter optimization in epilepsy patient models; data and software for mouse brain simulation; and extensive educational material. TVB cloud services facilitate reproducible online collaboration and discovery of data assets, models, and software embedded in scalable and secure workflows, a precondition for research on large cohort data sets, better generalizability, and clinical translation.

Normalizing the brain connectome for communication through synchronization.

Networks in neuroscience determine how brain function unfolds, and their perturbations lead to psychiatric disorders and brain disease. Brain networks are characterized by their connectomes, which comprise the totality of all connections, and are commonly described by graph theory. This approach is deeply rooted in a particle view of information processing, based on the quantification of informational bits such as firing rates. Oscillations and brain rhythms demand, however, a wave perspective of information processing based on synchronization. We extend traditional graph theory to a dual, particle-wave, perspective, integrate time delays due to finite transmission speeds, and derive a normalization of the connectome. When applied to the database of the Human Connectome Project, it explains the emergence of frequency-specific network cores including the visual and default mode networks. These findings are robust across human subjects (N = 100) and are a fundamental network property within the wave picture. The normalized connectome comprises the particle view in the limit of infinite transmission speeds and opens the applicability of graph theory to a wide range of novel network phenomena, including physiological and pathological brain rhythms. These two perspectives are orthogonal, but not incommensurable, when understood within the novel, here-proposed, generalized framework of structural connectivity.

Integrating psychosocial variables and societal diversity in epidemic models for predicting COVID-19 transmission dynamics.

During the current COVID-19 pandemic, governments must make decisions based on a variety of information including estimations of infection spread, health care capacity, economic and psychosocial considerations. The disparate validity of current short-term forecasts of these factors is a major challenge to governments. By causally linking an established epidemiological spread model with dynamically evolving psychosocial variables, using Bayesian inference we estimate the strength and direction of these interactions for German and Danish data of disease spread, human mobility, and psychosocial factors based on the serial cross-sectional COVID-19 Snapshot Monitoring (COSMO; N = 16,981). We demonstrate that the strength of cumulative influence of psychosocial variables on infection rates is of a similar magnitude as the influence of physical distancing. We further show that the efficacy of political interventions to contain the disease strongly depends on societal diversity, in particular group-specific sensitivity to affective risk perception. As a consequence, the model may assist in quantifying the effect and timing of interventions, forecasting future scenarios, and differentiating the impact on diverse groups as a function of their societal organization. Importantly, the careful handling of societal factors, including support to the more vulnerable groups, adds another direct instrument to the battery of political interventions fighting epidemic spread.

Experimental and Computational Study on Motor Control and Recovery After Stroke: Toward a Constructive Loop Between Experimental and Virtual Embodied Neuroscience.

Being able to replicate real experiments with computational simulations is a unique opportunity to refine and validate models with experimental data and redesign the experiments based on simulations. However, since it is technically demanding to model all components of an experiment, traditional approaches to modeling reduce the experimental setups as much as possible. In this study, our goal is to replicate all the relevant features of an experiment on motor control and motor rehabilitation after stroke. To this aim, we propose an approach that allows continuous integration of new experimental data into a computational modeling framework. First, results show that we could reproduce experimental object displacement with high accuracy via the simulated embodiment in the virtual world by feeding a spinal cord model with experimental registration of the cortical activity. Second, by using computational models of multiple granularities, our preliminary results show the possibility of simulating several features of the brain after stroke, from the local alteration in neuronal activity to long-range connectivity remodeling. Finally, strategies are proposed to merge the two pipelines. We further suggest that additional models could be integrated into the framework thanks to the versatility of the proposed approach, thus allowing many researchers to achieve continuously improved experimental design.

Dynamical Mechanisms of Interictal Resting-State Functional Connectivity in Epilepsy.

Drug-resistant focal epilepsy is a large-scale brain networks disorder characterized by altered spatiotemporal patterns of functional connectivity (FC), even during interictal resting state (RS). Although RS-FC-based metrics can detect these changes, results from RS functional magnetic resonance imaging (RS-fMRI) studies are unclear and difficult to interpret, and the underlying dynamical mechanisms are still largely unknown. To better capture the RS dynamics, we phenomenologically extended the neural mass model of partial seizures, the Epileptor, by including two neuron subpopulations of epileptogenic and nonepileptogenic type, making it capable of producing physiological oscillations in addition to the epileptiform activity. Using the neuroinformatics platform The Virtual Brain, we reconstructed 14 epileptic and 5 healthy human (of either sex) brain network models (BNMs), based on individual anatomical connectivity and clinically defined epileptogenic heatmaps. Through systematic parameter exploration and fitting to neuroimaging data, we demonstrated that epileptic brains during interictal RS are associated with lower global excitability induced by a shift in the working point of the model, indicating that epileptic brains operate closer to a stable equilibrium point than healthy brains. Moreover, we showed that functional networks are unaffected by interictal spikes, corroborating previous experimental findings; additionally, we observed higher excitability in epileptogenic regions, in agreement with the data. We shed light on new dynamical mechanisms responsible for altered RS-FC in epilepsy, involving the following two key factors: (1) a shift of excitability of the whole brain leading to increased stability; and (2) a locally increased excitability in the epileptogenic regions supporting the mixture of hyperconnectivity and hypoconnectivity in these areas.SIGNIFICANCE STATEMENT Advances in functional neuroimaging provide compelling evidence for epilepsy-related brain network alterations, even during the interictal resting state (RS). However, the dynamical mechanisms underlying these changes are still elusive. To identify local and network processes behind the RS-functional connectivity (FC) spatiotemporal patterns, we systematically manipulated the local excitability and the global coupling in the virtual human epileptic patient brain network models (BNMs), complemented by the analysis of the impact of interictal spikes and fitting to the neuroimaging data. Our results suggest that a global shift of the dynamic working point of the brain model, coupled with locally hyperexcitable node dynamics of the epileptogenic networks, provides a mechanistic explanation of the epileptic processes during the interictal RS period. These, in turn, are associated with the changes in FC.

Natural rhythms of periodic temporal attention.

That attention is a fundamentally rhythmic process has recently received abundant empirical evidence. The essence of temporal attention, however, is to flexibly focus in time. Whether this function is constrained by an underlying rhythmic neural mechanism is unknown. In six interrelated experiments, we behaviourally quantify the sampling capacities of periodic temporal attention during auditory or visual perception. We reveal the presence of limited attentional capacities, with an optimal sampling rate of ~1.4 Hz in audition and ~0.7 Hz in vision. Investigating the motor contribution to temporal attention, we show that it scales with motor rhythmic precision, maximal at ~1.7 Hz. Critically, motor modulation is beneficial to auditory but detrimental to visual temporal attention. These results are captured by a computational model of coupled oscillators, that reveals the underlying structural constraints governing the temporal alignment between motor and attention fluctuations.

Transmission time delays organize the brain network synchronization.

The timing of activity across brain regions can be described by its phases for oscillatory processes, and is of crucial importance for brain functioning. The structure of the brain constrains its dynamics through the delays due to propagation and the strengths of the white matter tracts. We use self-sustained delay-coupled, non-isochronous, nonlinearly damped and chaotic oscillators to study how spatio-temporal organization of the brain governs phase lags between the coherent activity of its regions. In silico results for the brain network model demonstrate a robust switching from in- to anti-phase synchronization by increasing the frequency, with a consistent lagging of the stronger connected regions. Relative phases are well predicted by an earlier analysis of Kuramoto oscillators, confirming the spatial heterogeneity of time delays as a crucial mechanism in shaping the functional brain architecture. Increased frequency and coupling are also shown to distort the oscillators by decreasing their amplitude, and stronger regions have lower, but more synchronized activity. These results indicate specific features in the phase relationships within the brain that need to hold for a wide range of local oscillatory dynamics, given that the time delays of the connectome are proportional to the lengths of the structural pathways. This article is part of the theme issue 'Nonlinear dynamics of delay systems'.

Controlling seizure propagation in large-scale brain networks.

Information transmission in the human brain is a fundamentally dynamic network process. In partial epilepsy, this process is perturbed and highly synchronous seizures originate in a local network, the so-called epileptogenic zone (EZ), before recruiting other close or distant brain regions. We studied patient-specific brain network models of 15 drug-resistant epilepsy patients with implanted stereotactic electroencephalography (SEEG) electrodes. Each personalized brain model was derived from structural data of magnetic resonance imaging (MRI) and diffusion tensor weighted imaging (DTI), comprising 88 nodes equipped with region specific neural mass models capable of demonstrating a range of epileptiform discharges. Each patient's virtual brain was further personalized through the integration of the clinically hypothesized EZ. Subsequent simulations and connectivity modulations were performed and uncovered a finite repertoire of seizure propagation patterns. Across patients, we found that (i) patient-specific network connectivity is predictive for the subsequent seizure propagation pattern; (ii) seizure propagation is characterized by a systematic sequence of brain states; (iii) propagation can be controlled by an optimal intervention on the connectivity matrix; (iv) the degree of invasiveness can be significantly reduced via the proposed seizure control as compared to traditional resective surgery. To stop seizures, neurosurgeons typically resect the EZ completely. We showed that stability analysis of the network dynamics, employing structural and dynamical information, estimates reliably the spatiotemporal properties of seizure propagation. This suggests novel less invasive paradigms of surgical interventions to treat and manage partial epilepsy.

Phase-lags in large scale brain synchronization: Methodological considerations and in-silico analysis.

Architecture of phase relationships among neural oscillations is central for their functional significance but has remained theoretically poorly understood. We use phenomenological model of delay-coupled oscillators with increasing degree of topological complexity to identify underlying principles by which the spatio-temporal structure of the brain governs the phase lags between oscillatory activity at distant regions. Phase relations and their regions of stability are derived and numerically confirmed for two oscillators and for networks with randomly distributed or clustered bimodal delays, as a first approximation for the brain structural connectivity. Besides in-phase, clustered delays can induce anti-phase synchronization for certain frequencies, while the sign of the lags is determined by the natural frequencies and by the inhomogeneous network interactions. For in-phase synchronization faster oscillators always phase lead, while stronger connected nodes lag behind the weaker during frequency depression, which consistently arises for in-silico results. If nodes are in anti-phase regime, then a distance π is added to the in-phase trends. The statistics of the phases is calculated from the phase locking values (PLV), as in many empirical studies, and we scrutinize the method's impact. The choice of surrogates do not affects the mean of the observed phase lags, but higher significance levels that are generated by some surrogates, cause decreased variance and might fail to detect the generally weaker coherence of the interhemispheric links. These links are also affected by the non-stationary and intermittent synchronization, which causes multimodal phase lags that can be misleading if averaged. Taken together, the results describe quantitatively the impact of the spatio-temporal connectivity of the brain to the synchronization patterns between brain regions, and to uncover mechanisms through which the spatio-temporal structure of the brain renders phases to be distributed around 0 and π. TRIAL REGISTRATION: South African Clinical Trials Register: http://www.sanctr.gov.za/SAClinicalbrnbspTrials/tabid/169/Default.aspx, then link to respiratory tract then link to tuberculosis, pulmonary; and TASK Applied Sciences Clinical Trials, AP-TB-201-16 (ALOPEXX): https://task.org.za/clinical-trials/.

The multiscale entropy: Guidelines for use and interpretation in brain signal analysis.

BACKGROUND: Multiscale entropy (MSE) estimates the predictability of a signal over multiple temporal scales. It has been recently applied to study brain signal variability, notably during aging. The grounds of its application and interpretation remain unclear and subject to debate. METHOD: We used both simulated and experimental data to provide an intuitive explanation of MSE and to explore how it relates to the frequency content of the signal, depending on the amount of (non)linearity and stochasticity in the underlying dynamics. RESULTS: The scaling and peak-structure of MSE curves relate to the scaling and peaks of the power spectrum in the presence of linear autocorrelations. MSE also captures nonlinear autocorrelations and their interactions with stochastic dynamical components. The previously reported crossing of young and old adults' MSE curves for EEG data appears to be mainly due to linear stochastic processes, and relates to young adults' EEG dynamics exhibiting a slower time constant. COMPARISON WITH EXISTING METHODS: We make the relationship between MSE curve and power spectrum as well as with a linear autocorrelation measure, namely multiscale root-mean-square-successive-difference, more explicit. MSE allows gaining insight into the time-structure of brain activity fluctuations. Its combined use with other metrics could prevent any misleading interpretations with regard to underlying stochastic processes. CONCLUSIONS: Although not straightforward, when applied to brain signals, the features of MSE curves can be linked to their power content and provide information about both linear and nonlinear autocorrelations that are present therein.

Heterogeneity of time delays determines synchronization of coupled oscillators.

Network couplings of oscillatory large-scale systems, such as the brain, have a space-time structure composed of connection strengths and signal transmission delays. We provide a theoretical framework, which allows treating the spatial distribution of time delays with regard to synchronization, by decomposing it into patterns and therefore reducing the stability analysis into the tractable problem of a finite set of delay-coupled differential equations. We analyze delay-structured networks of phase oscillators and we find that, depending on the heterogeneity of the delays, the oscillators group in phase-shifted, anti-phase, steady, and non-stationary clusters, and analytically compute their stability boundaries. These results find direct application in the study of brain oscillations.

Alterations in the coupling functions between cortical and cardio-respiratory oscillations due to anaesthesia with propofol and sevoflurane.

The precise mechanisms underlying general anaesthesia pose important and still open questions. To address them, we have studied anaesthesia induced by the widely used (intravenous) propofol and (inhalational) sevoflurane anaesthetics, computing cross-frequency coupling functions between neuronal, cardiac and respiratory oscillations in order to determine their mutual interactions. The phase domain coupling function reveals the form of the function defining the mechanism of an interaction, as well as its coupling strength. Using a method based on dynamical Bayesian inference, we have thus identified and analysed the coupling functions for six relationships. By quantitative assessment of the forms and strengths of the couplings, we have revealed how these relationships are altered by anaesthesia, also showing that some of them are differently affected by propofol and sevoflurane. These findings, together with the novel coupling function analysis, offer a new direction in the assessment of general anaesthesia and neurophysiological interactions, in general.