Search

Abstract Non-invasive neuromodulation represents a major opportunity for brain interventions, and transcranial focused ultrasound (FUS) is one of the most promising approaches. However, some challenges prevent the community from fully understanding its outcomes. We aimed to address one of them and unravel the temporal dynamics of FUS effects in humans. Twenty-two healthy volunteers participated in the study. Eleven received FUS in the right inferior frontal cortex while the other 11 were stimulated in the right thalamus. Using a temporal dynamic approach, we compared resting-state fMRI seed-based functional connectivity obtained before and after FUS. We also assessed behavioural changes as measured with a task of reactive motor inhibition. Our findings reveal that the effects of FUS are predominantly time-constrained and spatially distributed in brain regions functionally connected with the directly stimulated area. In addition, mediation analysis highlighted that FUS applied in the right inferior cortex was associated with behavioural alterations which was directly explained by the applied acoustic pressure and the brain functional connectivity change we observed. Our study underscored that the biological effects of FUS are indicative of behavioural changes observed more than an hour following stimulation and are directly related to the applied acoustic pressure.

Abstract Deep brain stimulation (DBS) can ameliorate motor symptoms in Parkinson’s disease (PD), but its mechanism remains unclear. This work constructs a multi-scale brain model using the fMRI data from 27 PD patients with subthalamic DBS and 30 healthy controls. The model fits microscopic coupling parameters in the cortico-basal ganglia-thalamic neural loop to match individual connectivity, finding the “push-pull” effect of basal ganglia network. Specifically, increased GABAergic projection into the thalamus from basal ganglia worsens rigidity, while reduced GABAergic projection within the cortex exacerbates bradykinesia, suggesting that the dopamine deficiency induces the chain coupling variations to “push” the network to an abnormal state. Conversely, DBS can alleviate rigidity by enhancing GABAergic projections within the basal ganglia, and improve bradykinesia by reducing cortical projections to basal ganglia, exhibiting that DBS “pulls” the network to a healthy state. This work combines the microscopic and macroscopic neural information for understanding PD and its treatment.

Abstract Low-frequency transcranial ultrasound stimulation (TUS) allows to alter brain functioning with a high spatial resolution and to reach deep targets. However, the time-course of TUS effects remains largely unknown. We applied TUS on three brain targets for three different monkeys: the anterior medial prefrontal cortex, the supplementary motor area and the perigenual anterior cingulate cortex. For each, one resting-state fMRI was acquired between 30 and 150 min after TUS as well as one without stimulation (control). We captured seed-based brain connectivity changes dynamically and on an individual basis. We also assessed between individuals and between targets homogeneity and brain features that predicted TUS changes. We found that TUS prompts heterogenous functional connectivity alterations yet retain certain consistent changes; we identified 6 time-courses of changes including transient and long duration alterations; with a notable degree of accuracy we found that brain alterations could partially be predicted. Altogether, our results highlight that TUS induces heterogeneous functional connectivity alterations. On a more technical point, we also emphasize the need to consider brain changes over-time rather than just observed during a snapshot; to consider inter-individual variability since changes could be highly different from one individual to another.

Flexible cognitive functions, such as working memory (WM), usually require a balance between localized and distributed information processing. However, it is challenging to uncover how local and distributed processing specifically contributes to task-induced activity in a region. Although the recently proposed activity flow mapping approach revealed the relative contribution of distributed processing, few studies have explored the adaptive and plastic changes that underlie cognitive manipulation. In this study, we recruited 51 healthy volunteers (31 females) and investigated how the activity flow and brain activation of the frontoparietal systems was modulated by WM load and training. While the activation of both executive control network (ECN) and dorsal attention network (DAN) increased linearly with memory load at baseline, the relative contribution of distributed processing showed a linear response only in the DAN, which was prominently attributed to within-network activity flow. Importantly, adaptive training selectively induced an increase in the relative contribution of distributed processing in the ECN and also a linear response to memory load, which were predominantly due to between-network activity flow. Furthermore, we demonstrated a causal effect of activity flow prediction through training manipulation on connectivity and activity. In contrast with classic brain activation estimation, our findings suggest that the relative contribution of distributed processing revealed by activity flow prediction provides unique insights into neural processing of frontoparietal systems under the manipulation of cognitive load and training. This study offers a new methodological framework for exploring information integration versus segregation underlying cognitive processing.

In this study, structural images of 1048 healthy subjects from the Human Connectome Project Young Adult study and 94 from ADNI-3 study were processed by an in-house tractography pipeline and analyzed together with pre-processed data of the same subjects from braingraph.org. Whole brain structural connectome features were used to build a simple correlation-based regression machine learning model to predict intelligence and age of healthy subjects. Our results showed that different forms of intelligence as well as age are predictable to a certain degree from diffusion tensor imaging detecting anatomical fiber tracts in the living human brain. Though we did not identify significant differences in the prediction capability for the investigated features depending on the imaging feature extraction method, we did find that crystallized intelligence was consistently better predictable than fluid intelligence from structural connectivity data through all datasets. Our findings suggest a practical and scalable processing and analysis framework to explore broader research topics employing brain MR imaging.

AbstractNeural systems are shaped by multiple constraints, balancing region communication with the cost of establishing and maintaining physical connections. It has been suggested that the lengths of neural projections be minimized, reducing their spatial and metabolic impact on the organism. However, long-range connections are prevalent in the connectomes across various species, and thus, rather than rewiring connections to reduce length, an alternative theory proposes that the brain minimizes total wiring length through a suitable positioning of regions, termed component placement optimization. Previous studies in nonhuman primates have refuted this idea by identifying a nonoptimal component placement, where a spatial rearrangement of brain regions in silico leads to a reduced total wiring length. Here, for the first time in humans, we test for component placement optimization. We show a nonoptimal component placement for all subjects in our sample from the Human Connectome Project (N = 280; aged 22–30 years; 138 females), suggesting the presence of constraints—such as the reduction of processing steps between regions—that compete with the elevated spatial and metabolic costs. Additionally, by simulating communication between brain regions, we argue that this suboptimal component placement supports dynamics that benefit cognition.

Significant research has investigated synchronisation in brain networks, but the bulk of this work has explored the contribution of brain networks at the macroscale. Here we explore the effects of changing network topology on functional dynamics in spatially constrained random networks representing mesoscale neocortex. We use the Kuramoto model to simulate network dynamics and explore synchronisation and critical dynamics of the system as a function of topology in randomly generated networks with a distance-related wiring probability and no preferential attachment term. We show networks which predominantly make short-distance connections smooth out the critical coupling point and show much greater metastability, resulting in a wider range of coupling strengths demonstrating critical dynamics and metastability. We show the emergence of cluster synchronisation in these geometrically-constrained networks with functional organisation occurring along structural connections that minimise the participation coefficient of the cluster. We show that these cohorts of internally synchronised nodes also behave en masse as weakly coupled nodes and show intra-cluster desynchronisation and resynchronisation events related to inter-cluster interaction. While cluster synchronisation appears crucial to healthy brain function, it may also be pathological if it leads to unbreakable local synchronisation which may happen at extreme topologies, with implications for epilepsy research, wider brain function and other domains such as social networks.

The analysis of whole brain networks started in the 1980s when only a handful of connectomes were available. In these early days, information about the human connectome was absent and one could only dream about having information about connectivity in a single human subject. Thanks to non-invasive methods such as diffusion imaging, we now know about connectivity in many species and, for some species, in many individuals. To illustrate the rapid change in availability of connectome data, the UK Biobank is on track to record structural and functional connectivity in 100,000 human subjects. Moreover, connectome data from a range of species is now available: from Caenorhabditis elegans and the fruit fly to pigeons, rodents, cats, non-human primates, and humans. This review will give a brief overview of what structural connectivity data is now available, how connectomes are organized, and how their organization shows common features across species. Finally, I will outline some of the current challenges and potential future work in making use of connectome information.

Chen at al used fMRI to study the brain connectivity in two ageing cohorts and a mixed cohort of patients with neurodegenerative diseases. They demonstrated that with ageing the brain moves towards a more segregated modular structure, however this is disrupted by neurodegenerative diseases.

Abstract Late-life depression (LLD) may increase the risk of Alzheimer’s dementia (AD). While amyloidopathy accelerates AD progression, its role in such patients has not yet been elucidated. We hypothesized that cerebral amyloidopathy distinctly affects the alteration of brain network topology and may be associated with distinct cognitive symptoms. We recruited 26 and 27 depressed mild cognitive impairment (MCI) patients with (LLD-MCI-A(+)) and without amyloid accumulation (LLD-MCI-A(−)), respectively, and 21 normal controls. We extracted structural brain networks using their diffusion-weighted images. We aimed to compare the distinct network deterioration in LLD-MCI with and without amyloid accumulation and the relationship with their distinct cognitive decline. Thus, we performed a group comparison of the network topological measures and investigated any correlations with neurocognitive testing scores. Topological features of brain networks were different according to the presence of amyloid accumulation. Disrupted network connectivity was highly associated with impaired recall and recognition in LLD-MCI-A(+) patients. Inattention and dysexecutive function were more influenced by the altered networks involved in fronto-limbic circuitry dysfunction in LLD-MCI-A(−) patients. Our results show that alterations in brain network topology may reflect different cognitive dysfunction depending on amyloid accumulation in depressed older adults with MCI.

Natural vision engages a wide range of higher-level regions that integrate visual information over the large-scale brain network. How interareal connectivity reconfigures during the processing of ongoing natural visual scenes and how these dynamic functional changes relate to the underlaying anatomical links between regions is not well understood. Here, we hypothesized that macaque visual brain regions are poly-functional sharing the capacity to change their configuration state depending on the nature of visual input. To address this hypothesis, we reconstructed networks from in-vivo diffusion-weighted imaging (DWI) and functional magnetic resonance imaging (fMRI) data obtained in four alert macaque monkeys viewing naturalistic movie scenes. At first, we characterized network properties and found greater interhemispheric density and greater inter-subject variability in free-viewing networks as compared to structural networks. From the structural connectivity, we then captured modules on which we identified hubs during free-viewing that formed a widespread visuo-saccadic network across frontal (FEF, 46v), parietal (LIP, Tpt), and occipitotemporal modules (MT, V4, TEm), and that excluded primary visual cortex. Inter-subject variability of well-connected hubs reflected subject-specific configurations that largely recruited occipito-parietal and frontal modules. Across the cerebral hemispheres, free-viewing networks showed higher correlations among long-distance brain regions as compared to structural networks. From these findings, we hypothesized that long-distance interareal connectivity could reconfigure depending on the ongoing changes in visual scenes. Testing this hypothesis by applying temporally resolved functional connectivity we observed that many structurally defined areas (such as areas V4, MT/MST and LIP) were poly-functional as they were recruited as hub members of multiple network states that changed during the presentation of scenes containing objects, motion, faces, and actions. We suggest that functional flexibility in macaque macroscale brain networks is required for the efficient interareal communication during active natural vision. To further promote the use of naturalistic free-viewing paradigms and increase the development of macaque neuroimaging resources, we share our datasets in the PRIME-DE consortium.

Neuromodulation is an established treatment for numerous neurological conditions, but to expand the therapeutic scope there is a need to improve the spatial, temporal and cell-type specificity of stimulation. Optogenetics is a promising area of current research, enabling optical stimulation of genetically-defined cell types without interfering with concurrent electrical recording for closed-loop control of neural activity. We are developing an open-source system to provide a platform for closed-loop optogenetic neuromodulation, incorporating custom integrated circuitry for recording and stimulation, real-time closed-loop algorithms running on a microcontroller and experimental control via a PC interface. We include commercial components to validate performance, with the ultimate aim of translating this approach to humans. In the meantime our system is flexible and expandable for use in a variety of preclinical neuroscientific applications. The platform consists of a Controlling Abnormal Network Dynamics using Optogenetics (CANDO) Control System (CS) that interfaces with up to four CANDO headstages responsible for electrical recording and optical stimulation through custom CANDO LED optrodes. Control of the hardware, inbuilt algorithms and data acquisition is enabled via the CANDO GUI (Graphical User Interface). Here we describe the design and implementation of this system, and demonstrate how it can be used to modulate neuronal oscillations in vitro and in vivo.

Flexibility is a hallmark of human intelligence. Emerging studies have proposed several flexibility measurements at the level of individual regions, to produce a brain map of neural flexibility. However, flexibility is usually inferred from separate components of brain activity (i.e., intrinsic/task-evoked), and different definitions are used. Moreover, recent studies have argued that neural processing may be more than a task-driven and intrinsic dichotomy. Therefore, the understanding to neural flexibility is still incomplete. To address this issue, we propose a multifaceted definition of neural flexibility according to three key features: broad cognitive engagement, distributed connectivity, and adaptive connectome dynamics. For these three features, we first review the advances in computational approaches, their functional relevance, and their potential pitfalls. We then suggest a set of metrics that can help us assign a flexibility rating to each region. Subsequently, we present an emergent probabilistic view for further understanding the functional operation of individual regions in the unified framework of intrinsic and task-driven states. Finally, we highlight several areas related to the multifaceted definition of neural flexibility for future research. This review not only strengthens our understanding of flexible human brain, but also suggests that the measure of neural flexibility could bridge the gap between understanding intrinsic and task-driven brain function dynamics.

Abstract Introduction The pathophysiological process of Alzheimer's disease is thought to begin years before clinical decline, with evidence suggesting prion‐like spreading processes of neurofibrillary tangles and amyloid plaques. Methods Using diffusion magnetic resonance imaging data from the Alzheimer's Disease Neuroimaging Initiative database, we first identified relevant features for dementia diagnosis. We then created dynamic models with the Nathan Kline Institute‐Rockland Sample database to estimate the earliest detectable stage associated with dementia in the simulated disease progression. Results A classifier based on centrality measures provides informative predictions. Strength and closeness centralities are the most discriminative features, which are associated with the medial temporal lobe and subcortical regions, together with posterior and occipital brain regions. Our model simulations suggest that changes associated with dementia begin to manifest structurally at early stages. Discussion Our analyses suggest that diffusion magnetic resonance imaging–based centrality measures can offer a tool for early disease detection before clinical dementia onset.

Resting-state networks usually show time-limited self-sustaining oscillatory patterns (TLSOPs) with the characteristic features of multiscaled rhythms and frequent switching between different rhythms, but the underlying mechanisms remain unclear. To reveal the mechanisms of multiscaled rhythms, we present a simplified reaction-diffusion model of activation propagation to reproduce TLSOPs in real brain networks. We find that the reproduced TLSOPs do show multiscaled rhythms, depending on the activating threshold and initially chosen activating nodes. To understand the frequent switching between different rhythms, we present an approach of dominant activation paths and find that the multiscaled rhythms can be separated into individual rhythms denoted by different core networks, and the switching between them can be implemented by a time-dependent activating threshold. Further, based on the microstates of TLSOPs, we introduce the concept of a return loop to study the distribution of the return times of microstates in TLSOPs and find that it satisfies the Weibull distribution. Then, to check it for real data, we present a method of a shifting window to transform a continuous time series into a discrete two-state time series and interestingly find that the Weibull distribution also exists in resting-state EEG and fMRI data. Finally, we show that the TLSOP lifetime depends exponentially on the core network size and can be explained by a theory of the complete graphs.

Structural covariance analysis is a widely used structural MRI analysis method which characterises the co-relations of morphology between brain regions over a group of subjects. To our knowledge, little has been investigated in terms of the comparability of results between different data sets of healthy human subjects, as well as the reliability of results over the same subjects in different rescan sessions, image resolutions, or FreeSurfer versions.In terms of comparability, our results show substantial differences in the structural covariance matrix between data sets of age- and sex-matched healthy human adults. These differences persist after univariate site correction, they are exacerbated by low sample sizes, and they are most pronounced when using average cortical thickness as a morphological measure. Down-stream graph theoretic analyses further show statistically significant differences.In terms of reliability, substantial differences were also found when comparing repeated scan sessions of the same subjects, image resolutions, and even FreeSurfer versions of the same image. We could further estimate the relative measurement error and showed that it is largest when using cortical thickness as a morphological measure. Using simulated data, we argue that cortical thickness is least reliable because of larger relative measurement errors.Practically, we make the following recommendations (1) combining subjects across sites into one group should be avoided, particularly if sites differ in image resolutions, subject demographics, or preprocessing steps; (2) surface area and volume should be preferred as morphological measures over cortical thickness; (3) a large number of subjects (n≫30 for the Desikan-Killiany parcellation) should be used to estimate structural covariance; (4) measurement error should be assessed where repeated measurements are available; (5) if combining sites is critical, univariate (per ROI) site-correction is insufficient, but error covariance (between ROIs) should be explicitly measured and modelled.

Computational studies of the influence of different network parameters on the dynamic and topological network effects of brain stimulation can enhance our understanding of different outcomes between individuals. In this study, a brain stimulation session along with the subsequent post-stimulation brain activity is simulated for a period of one day using a network of modified Wilson-Cowan oscillators coupled according to diffusion imaging based structural connectivity. We use this computational model to examine how differences in the inter-region connectivity and the excitability of stimulated regions at the time of stimulation can affect post-stimulation behaviours. Our findings indicate that the initial inter-region connectivity can heavily affect the changes that stimulation induces in the connectivity of the network. Moreover, differences in the excitability of the stimulated regions seem to lead to different post-stimulation connectivity changes across the model network, including on the internal connectivity of non-stimulated regions.

Neuronal circuits can be modelled in detail allowing us to predict the effects of stimulation on individual neurons. Electrical stimulation of neuronal circuits in vitro and in vivo excites a range of neurons within the tissue and measurements of neural activity, e.g the local field potential (LFP), are again an aggregate of a large pool of cells. The previous version of our Virtual Electrode Recording Tool for EXtracellular Potentials (VERTEX) allowed for the simulation of the LFP generated by a patch of brain tissue. Here, we extend VERTEX to simulate the effect of electrical stimulation through a focal electric field. We observe both direct changes in neural activity and changes in synaptic plasticity. Testing our software in a model of a rat neocortical slice, we determine the currents contributing to the LFP, the effects of paired pulse stimulation to induce short term plasticity (STP), and the effect of theta burst stimulation (TBS) to induce long term potentiation (LTP).

Patients with idiopathic generalised epilepsy (IGE) typically have normal conventional magnetic resonance imaging (MRI), hence diagnosis based on MRI is challenging. Anatomical abnormalities underlying brain dysfunctions in IGE are unclear and their relation to the pathomechanisms of epileptogenesis is poorly understood.In this study, we applied connectometry, an advanced quantitative neuroimaging technique for investigating localised changes in white-matter tissues in vivo. Analysing white matter structures of 32 subjects we incorporated our in vivo findings in a computational model of seizure dynamics to suggest a plausible mechanism of epileptogenesis.Patients with IGE have significant bilateral alterations in major white-matter fascicles. In the cingulum, fornix, and superior longitudinal fasciculus, tract integrity is compromised, whereas in specific parts of tracts between thalamus and the precentral gyrus, tract integrity is enhanced in patients. Combining these alterations in a logistic regression model, we computed the decision boundary that discriminated patients and controls. The computational model, informed with the findings on the tract abnormalities, specifically highlighted the importance of enhanced cortico-reticular connections along with impaired cortico-cortical connections in inducing pathological seizure-like dynamics.We emphasise taking directionality of brain connectivity into consideration towards understanding the pathological mechanisms; this is possible by combining neuroimaging and computational modelling. Our imaging evidence of structural alterations suggest the loss of cortico-cortical and enhancement of cortico-thalamic fibre integrity in IGE. We further suggest that impaired connectivity from cortical regions to the thalamic reticular nucleus offers a therapeutic target for selectively modifying the brain circuit for reversing the mechanisms leading to epileptogenesis. Keywords: Computational model, Diagnosis, Diffusion MRI, Epilepsy mechanism, Generalised epilepsy

We studied the dynamic functional connectivity profile of dementia with Lewy bodies (DLB) and Alzheimer's disease (AD) compared to controls, how it differs between the two dementia subtypes, and a possible relation between dynamic connectivity alterations and temporally transient clinical symptoms in DLB.Resting state fMRI data from 31 DLB, 29 AD, and 31 healthy control participants were analyzed using dual regression to determine between-network functional connectivity. Subsequently, we used a sliding window approach followed by k-means clustering and dynamic network analyses to study dynamic functional connectivity. Dynamic connectivity measures that showed significant group differences were tested for correlations with clinical symptom severity.Our results show that AD and DLB patients spent more time than controls in sparse connectivity configurations with absence of strong positive and negative connections and a relative isolation of motor networks from other networks. Additionally, DLB patients spent less time in a more strongly connected state and the variability of global brain network efficiency was reduced in DLB compared to controls. There were no significant correlations between dynamic connectivity measures and clinical symptom severity.An inability to switch out of states of low inter-network connectivity into more highly and specifically connected network configurations might be related to the presence of dementia in general as it was observed in both AD and DLB. In contrast, the loss of global efficiency variability in DLB might indicate the presence of an abnormally rigid brain network and the lack of economical dynamics, factors which could contribute to cognitive slowing and an inability to respond appropriately to situational demands. Keywords: Resting state fMRI, Dual regression, Sliding window, Dynamic network analysis, Neurodegeneration, Cognitive fluctuations

Lewy body dementia (LBD) and Alzheimer’s disease (AD) are common forms of dementia that have different clinical profiles but are both commonly associated with attentional deficits. The aim of this study was to investigate efficiency of different attentional systems in LBD and AD and its association with brain structural abnormalities. We studied reaction time (RT) data from 45 LBD, 31 AD patients and 22 healthy controls (HCs) using the Attention Network Test (ANT) to assess the efficiency of three different attentional systems: alerting, orienting and executive conflict. Voxel-based morphometry (VBM) was used to investigate relations between different attention components and cortical volume. Both dementia groups showed slower overall RTs than controls, with additional slowing in LBD relative to AD. There was a significant alerting effect in controls which was absent in the dementia groups, the executive conflict effect was greater in both dementia groups compared to controls, but the orienting effect did not differ between groups. Mean RT in AD was negatively correlated with occipital gray matter (GM) volume and in LBD orienting efficiency was negatively related to occipital white matter (WM) volume. Given that previous studies in less impaired patients suggest a maintenance of the alerting effect, the absent alerting effect in our study suggests a loss of alerting efficiency with dementia progression. While orienting was largely preserved, it might be related to occipital structural abnormalities in LBD. Executive function was markedly impaired in both dementia groups, however, the absence of relations to brain volume suggests that it might be more related to functional rather than macrostructural pathophysiological changes.

Background: Previous studies reported reduced volumes of many brain regions for temporal lobe epilepsy (TLE). It has also been suggested that there may be widespread changes in network features of TLE patients. It is not fully understood, however, how these two observations are related. Methods: Using magnetic resonance imaging data, we perform parcellation of the brains of 22 patients with left TLE and 39 non-epileptic controls. In each parcellated region of interest (ROI) we computed the surface area and, using diffusion tensor imaging and deterministic tractography, infer the number of streamlines and their average length between each pair of connected ROIs. For comparison to previous studies, we use a connectivity ‘weight’ and investigate how ROI surface area, number of streamlines & mean streamline length contribute to such weight. Results: We find that although there are widespread significant changes in surface area and position of ROIs in patients compared to controls, the changes in connectivity are much more subtle. Significant changes in connectivity weight can be accounted for by decreased surface area and increased streamline count. Conclusion: Changes in the surface area of ROIs can be a reliable biomarker for TLE with a large influence on connectivity. However, changes in structural connectivity via white matter streamlines are more subtle with a relatively lower influence on connection weights.

Focal seizures are episodes of pathological brain activity that appear to arise from a localised area of the brain. The onset patterns of focal seizure activity have been studied intensively, and they have largely been distinguished into two types-low amplitude fast oscillations (LAF), or high amplitude spikes (HAS). Here we explore whether these two patterns arise from fundamentally different mechanisms. Here, we use a previously established computational model of neocortical tissue, and validate it as an adequate model using clinical recordings of focal seizures. We then reproduce the two onset patterns in their most defining properties and investigate the possible mechanisms underlying the different focal seizure onset patterns in the model. We show that the two patterns are associated with different mechanisms at the spatial scale of a single ECoG electrode. The LAF onset is initiated by independent patches of localised activity, which slowly invade the surrounding tissue and coalesce over time. In contrast, the HAS onset is a global, systemic transition to a coexisting seizure state triggered by a local event. We find that such a global transition is enabled by an increase in the excitability of the "healthy" surrounding tissue, which by itself does not generate seizures, but can support seizure activity when incited. In our simulations, the difference in surrounding tissue excitability also offers a simple explanation of the clinically reported difference in surgical outcomes. Finally, we demonstrate in the model how changes in tissue excitability could be elucidated, in principle, using active stimulation. Taken together, our modelling results suggest that the excitability of the tissue surrounding the seizure core may play a determining role in the seizure onset pattern, as well as in the surgical outcome.

Many real-world networks contain highly connected nodes called hubs. Hubs are often crucial for network function and spreading dynamics. However, classical models of how hubs originate during network development unrealistically assume that new nodes attain information about the connectivity (for example the degree) of existing nodes. Here, we introduce hub formation through nonlinear growth where the number of nodes generated at each stage increases over time and new nodes form connections independent of target node features. Our model reproduces variation in number of connections, hub occurrence time, and rich-club organization of networks ranging from protein–protein, neuronal and fibre tract brain networks to airline networks. Moreover, nonlinear growth gives a more generic representation of these networks compared with previous preferential attachment or duplication–divergence models. Overall, hub creation through nonlinear network expansion can serve as a benchmark model for studying the development of many real-world networks.

Cognitive fluctuations are a core symptom in dementia with Lewy bodies (DLB) and may relate to pathological alterations in distributed brain networks. To test this we analysed resting state fMRI changes in a cohort of fluctuating DLB patients (n = 16) compared with age matched controls (n = 17) with the aim of finding functional connectivity (FC) differences between these two groups and whether these associate with cognitive fluctuations in DLB. Resting state networks (RSNs) were estimated using independent component analysis and FC between the RSN maps and the entirety of the brain was assessed using dual regression. The default mode network (DMN) appeared unaffected in DLB compared to controls but significant cluster differences between DLB and controls were found for the left fronto-parietal, temporal, and sensory–motor networks. Desynchronization of a number of cortical and subcortical areas related to the left fronto-parietal network was associated with the severity and frequency of cognitive fluctuations. Our findings provide empirical evidence for the potential role of attention–executive networks in the aetiology of this core symptom in DLB.

Temporal lobe epilepsy (TLE) is a prevalent neurological disorder resulting in disruptive seizures. In the case of drug resistant epilepsy resective surgery is often considered. This is a procedure hampered by unpredictable success rates, with many patients continuing to have seizures even after surgery. In this study we apply a computational model of epilepsy to patient specific structural connectivity derived from diffusion tensor imaging (DTI) of 22 individuals with left TLE and 39 healthy controls. We validate the model by examining patient-control differences in simulated seizure onset time and network location. We then investigate the potential of the model for surgery prediction by performing in silico surgical resections, removing nodes from patient networks and comparing seizure likelihood post-surgery to pre-surgery simulations. We find that, first, patients tend to transit from non-epileptic to epileptic states more often than controls in the model. Second, regions in the left hemisphere (particularly within temporal and subcortical regions) that are known to be involved in TLE are the most frequent starting points for seizures in patients in the model. In addition, our analysis also implicates regions in the contralateral and frontal locations which may play a role in seizure spreading or surgery resistance. Finally, the model predicts that patient-specific surgery (resection areas chosen on an individual, model-prompted, basis and not following a predefined procedure) may lead to better outcomes than the currently used routine clinical procedure. Taken together this work provides a first step towards patient specific computational modelling of epilepsy surgery in order to inform treatment strategies in individuals.

We describe some general results that constrain the dynamical fluctuations that can occur in non-equilibrium steady states, with a focus on molecular dynamics. That is, we consider Hamiltonian systems, coupled to external heat baths, and driven out of equilibrium by non-conservative forces. We focus on the probabilities of rare events (large deviations). First, we discuss a PT (parity-time) symmetry that appears in ensembles of trajectories where a current is constrained to have a large (non-typical) value. We analyse the heat flow in such ensembles, and compare it with non-equilibrium steady states. Second, we consider pathwise large deviations that are defined by considering many copies of a system. We show how the probability currents in such systems can be decomposed into orthogonal contributions that are related to convergence to equilibrium and to dissipation. We discuss the implications of these results for modelling non-equilibrium steady states.

Evidence from anatomical and functional imaging studies have highlighted major modifications of cortical circuits during adolescence. These include reductions of gray matter (GM), increases in the myelination of cortico-cortical connections and changes in the architecture of large-scale cortical networks. It is currently unclear, however, how the ongoing developmental processes impact upon the folding of the cerebral cortex and how changes in gyrification relate to maturation of GM/WM-volume, thickness and surface area. In the current study, we acquired high-resolution (3 Tesla) magnetic resonance imaging (MRI) data from 79 healthy subjects (34 males and 45 females) between the ages of 12 and 23 years and performed whole brain analysis of cortical folding patterns with the gyrification index (GI). In addition to GI-values, we obtained estimates of cortical thickness, surface area, GM and white matter (WM) volume which permitted correlations with changes in gyrification. Our data show pronounced and widespread reductions in GI-values during adolescence in several cortical regions which include precentral, temporal and frontal areas. Decreases in gyrification overlap only partially with changes in the thickness, volume and surface of GM and were characterized overall by a linear developmental trajectory. Our data suggest that the observed reductions in GI-values represent an additional, important modification of the cerebral cortex during late brain maturation which may be related to cognitive development.

Glioma is the most common form of primary brain tumor. Demographically, the risk of occurrence increases until old age. Here we present a novel computational model to reproduce the probability of glioma incidence across the lifespan. Previous mathematical models explaining glioma incidence are framed in a rather abstract way, and do not directly relate to empirical findings. To decrease this gap between theory and experimental observations, we incorporate recent data on cellular and molecular factors underlying gliomagenesis. Since evidence implicates the adult neural stem cell as the likely cell-of-origin of glioma, we have incorporated empirically-determined estimates of neural stem cell number, cell division rate, mutation rate and oncogenic potential into our model. We demonstrate that our model yields results which match actual demographic data in the human population. In particular, this model accounts for the observed peak incidence of glioma at approximately 80 years of age, without the need to assert differential susceptibility throughout the population. Overall, our model supports the hypothesis that glioma is caused by randomly-occurring oncogenic mutations within the neural stem cell population. Based on this model, we assess the influence of the (experimentally indicated) decrease in the number of neural stem cells and increase of cell division rate during aging. Our model provides multiple testable predictions, and suggests that different temporal sequences of oncogenic mutations can lead to tumorigenesis. Finally, we conclude that four or five oncogenic mutations are sufficient for the formation of glioma.

Complex networks have been characterised by their specific connectivity patterns (network motifs), but their building blocks can also be identified and described by node-motifs-a combination of local network features. One technique to identify single node-motifs has been presented by Costa et al. (L. D. F. Costa, F. A. Rodrigues, C. C. Hilgetag, and M. Kaiser, Europhys. Lett., 87, 1, 2009). Here, we first suggest improvements to the method including how its parameters can be determined automatically. Such automatic routines make high-throughput studies of many networks feasible. Second, the new routines are validated in different network-series. Third, we provide an example of how the method can be used to analyse network time-series. In conclusion, we provide a robust method for systematically discovering and classifying characteristic nodes of a network. In contrast to classical motif analysis, our approach can identify individual components (here: nodes) that are specific to a network. Such special nodes, as hubs before, might be found to play critical roles in real-world networks.

The nematode Caenorhabditis elegans, with information on neural connectivity, three-dimensional position and cell linage, provides a unique system for understanding the development of neural networks. Although C. elegans has been widely studied in the past, we present the first statistical study from a developmental perspective, with findings that raise interesting suggestions on the establishment of long-distance connections and network hubs. Here, we analyze the neuro-development for temporal and spatial features, using birth times of neurons and their three-dimensional positions. Comparisons of growth in C. elegans with random spatial network growth highlight two findings relevant to neural network development. First, most neurons which are linked by long-distance connections are born around the same time and early on, suggesting the possibility of early contact or interaction between connected neurons during development. Second, early-born neurons are more highly connected (tendency to form hubs) than later-born neurons. This indicates that the longer time frame available to them might underlie high connectivity. Both outcomes are not observed for random connection formation. The study finds that around one-third of electrically coupled long-range connections are late forming, raising the question of what mechanisms are involved in ensuring their accuracy, particularly in light of the extremely invariant connectivity observed in C. elegans. In conclusion, the sequence of neural network development highlights the possibility of early contact or interaction in securing long-distance and high-degree connectivity.

An essential requirement for the representation of functional patterns in complex neural networks, such as the mammalian cerebral cortex, is the existence of stable regimes of network activation, typically arising from a limited parameter range. In this range of limited sustained activity (LSA), the activity of neural populations in the network persists between the extremes of either quickly dying out or activating the whole network. Hierarchical modular networks were previously found to show a wider parameter range for LSA than random or small-world networks not possessing hierarchical organization or multiple modules. Here we explored how variation in the number of hierarchical levels and modules per level influenced network dynamics and occurrence of LSA. We tested hierarchical configurations of different network sizes, approximating the large-scale networks linking cortical columns in one hemisphere of the rat, cat, or macaque monkey brain. Scaling of the network size affected the number of hierarchical levels and modules in the optimal networks, also depending on whether global edge density or the numbers of connections per node were kept constant. For constant edge density, only few network configurations, possessing an intermediate number of levels and a large number of modules, led to a large range of LSA independent of brain size. For a constant number of node connections, there was a trend for optimal configurations in larger-size networks to possess a larger number of hierarchical levels or more modules. These results may help to explain the trend to greater network complexity apparent in larger brains and may indicate that this complexity is required for maintaining stable levels of neural activation.

Multiple local neuronal circuits support different, discrete frequencies of network rhythm in neocortex. Relationships between different frequencies correspond to mechanisms designed to minimise interference, couple activity via stable phase interactions, and control the amplitude of one frequency relative to the phase of another. These mechanisms are proposed to form a framework for spectral information processing. Individual local circuits can also transform their frequency through changes in intrinsic neuronal properties and interactions with other oscillating microcircuits. Here we discuss a frequency transformation in which activity in two coactive local circuits may combine sequentially to generate a third frequency whose period is the concatenation sum of the original two. With such an interaction, the intrinsic periodicity in each component local circuit is preserved – alternate, single periods of each original rhythm form one period of a new frequency - suggesting a robust mechanism for combining information processed on multiple concurrent spatiotemporal scales.

It has been suggested that neural systems across several scales of organization show optimal component placement, in which any spatial rearrangement of the components would lead to an increase of total wiring. Using extensive connectivity datasets for diverse neural networks combined with spatial coordinates for network nodes, we applied an optimization algorithm to the network layouts, in order to search for wire-saving component rearrangements. We found that optimized component rearrangements could substantially reduce total wiring length in all tested neural networks. Specifically, total wiring among 95 primate (Macaque) cortical areas could be decreased by 32%, and wiring of neuronal networks in the nematode Caenorhabditis elegans could be reduced by 48% on the global level, and by 49% for neurons within frontal ganglia. Wiring length reductions were possible due to the existence of long-distance projections in neural networks. We explored the role of these projections by comparing the original networks with minimally rewired networks of the same size, which possessed only the shortest possible connections. In the minimally rewired networks, the number of processing steps along the shortest paths between components was significantly increased compared to the original networks. Additional benchmark comparisons also indicated that neural networks are more similar to network layouts that minimize the length of processing paths, rather than wiring length. These findings suggest that neural systems are not exclusively optimized for minimal global wiring, but for a variety of factors including the minimization of processing steps.