Frontiers in network physiology最新文献

The study of specific physiological processes from the perspective of network physiology has gained recent attention. Modeling the global information integration among the separated functionalized modules in structural and functional brain networks is a central problem. In this article, the preferentially cutting-rewiring operation (PCRO) is introduced to approximatively describe the above physiological process, which consists of the cutting procedure and the rewiring procedure with specific preferential constraints. By applying the PCRO on the classical Erdös-Rényi random network (ERRN), three types of isolated nodes are generated, based on which the common leaves (CLs) are formed between the two hubs. This makes the initially homogeneous ERRN experience drastic changes and become heterogeneous. Importantly, a statistical analysis method is proposed to theoretically analyze the statistical properties of an ERRN with a PCRO. Specifically, the probability distributions of these three types of isolated nodes are derived, based on which the probability distribution of the CLs can be obtained easily. Furthermore, the validity and universality of our statistical analysis method have been confirmed in numerical experiments. Our contributions may shed light on a new perspective in the interdisciplinary field of complexity science and biological science and would be of great and general interest to network physiology.

Introduction: Statistical shape analysis (SSA) with clustering is often used to objectively define and categorise anatomical shape variations. However, studies until now have often focused on simplified anatomical reconstructions, despite the complexity of studied anatomies. This work aims to provide insights on the anatomical detail preservation required for SSA of highly diverse and complex anatomies, with particular focus on the left atrial appendage (LAA). This anatomical region is clinically relevant as the location of almost all left atrial thrombi forming during atrial fibrillation (AF). Moreover, its highly patient-specific complex architecture makes its clinical classification especially subjective.

Methods: Preliminary LAA meshes were automatically detected after robust image selection and wider left atrial segmentation. Following registration, four additional LAA mesh datasets were created as reductions of the preliminary dataset, with surface reconstruction based on reduced sample point densities. Utilising SSA model parameters determined to optimally represent the preliminary dataset, SSA model performance for the four simplified datasets was calculated. A representative simplified dataset was selected, and clustering analysis and performance were evaluated (compared to clinical labels) between the original trabeculated LAA anatomy and the representative simplification.

Results: As expected, simplified anatomies have better SSA evaluation scores (compactness, specificity and generalisation), corresponding to simpler LAA shape representation. However, oversimplification of shapes may noticeably affect 3D model output due to differences in geometric correspondence. Furthermore, even minor simplification may affect LAA shape clustering, where the adjusted mutual information (AMI) score of the clustered trabeculated dataset was 0.67, in comparison to 0.12 for the simplified dataset.

Discussion: This study suggests that greater anatomical preservation for complex and diverse LAA morphologies, currently neglected, may be more useful for shape categorisation via clustering analyses.

The concept of multifunctionality has enabled reservoir computers (RCs), a type of dynamical system that is typically realized as an artificial neural network, to reconstruct multiple attractors simultaneously using the same set of trained weights. However, there are many additional phenomena that arise when training a RC to reconstruct more than one attractor. Previous studies have found that in certain cases, if the RC fails to reconstruct a coexistence of attractors, then it exhibits a form of metastability, whereby, without any external input, the state of the RC switches between different modes of behavior that resemble the properties of the attractors it failed to reconstruct. In this paper, we explore the origins of these switching dynamics in a paradigmatic setting via the "seeing double" problem.

The stability of wave conduction in the heart is strongly related to the proper interplay between the electrophysiological activation and mechanical contraction of myocytes and extracellular matrix (ECM) properties. In this study, we statistically compare bioengineered cardiac tissues cultured on soft hydrogels ( $E\simeq 12$ kPa) and rigid glass substrates by focusing on the critical threshold of alternans, network-physiological tissue properties, and the formation of stable spiral waves that manifest after wave breakups. For the classification of wave dynamics, we use an improved signal oversampling technique and introduce simple probability maps to identify and visualize spatially concordant and discordant alternans as V- and X-shaped probability distributions. We found that cardiac tissues cultured on ECM-mimicking soft hydrogels show a lower variability of the calcium transient durations among cells in the tissue. This lowers the likelihood of forming stable spiral waves because of the larger dynamical range that tissues can be stably entrained with to form alternans and larger spatial spiral tip movement that increases the chance of self-termination on the tissue boundary. Conclusively, we show that a dysfunction in the excitation-contraction coupling dynamics facilitates life-threatening arrhythmic states such as spiral waves and, thus, highlights the importance of the network-physiological interplay between contractile myocytes and the ECM.

For patients with refractory epilepsy, the seizure onset zone (SOZ) plays an essential role in determining the specific regions of the brain that will be surgically resected. High-frequency oscillations (HFOs) and connectivity-based approaches have been identified among the potential biomarkers to localize the SOZ. However, there is no consensus on how connectivity between HFO events should be estimated, nor on its subject-specific short-term reliability. Therefore, we propose the channel-level connectivity dispersion (CLCD) as a metric to quantify the variability in synchronization between individual electrodes and to identify clusters of electrodes with abnormal synchronization, which we hypothesize to be associated with the SOZ. In addition, we developed a specialized filtering method that reduces oscillatory components caused by filtering broadband artifacts, such as sharp transients, spikes, or direct current shifts. Our connectivity estimates are therefore robust to the presence of these waveforms. To calculate our metric, we start by creating binary signals indicating the presence of high-frequency bursts in each channel, from which we calculate the pairwise connectivity between channels. Then, the CLCD is calculated by combining the connectivity matrices and measuring the variability in each electrode's combined connectivity values. We test our method using two independent open-access datasets comprising intracranial electroencephalography signals from 89 to 15 patients with refractory epilepsy, respectively. Recordings in these datasets were sampled at approximately 1000 Hz, and our proposed CLCDs were estimated in the ripple band (80-200 Hz). Across all patients in the first dataset, the average ROC-AUC was 0.73, and the average Cohen's d was 1.05, while in the second dataset, the average ROC-AUC was 0.78 and Cohen's d was 1.07. On average, SOZ channels had lower CLCD values than non-SOZ channels. Furthermore, based on the second dataset, which includes surgical outcomes (Engel I-IV), our analysis suggested that higher CLCD interquartile (as a measure of CLCD distribution spread) is associated with favorable outcomes (Engel I). This suggests that CLCD could significantly assist in identifying SOZ clusters and, therefore, provide an additional tool in surgical planning for epilepsy patients.

This study investigates the dynamic interactions between physiological systems during competitive gaming, utilizing a Network Physiology approach. By examining the physiological responses of a gamer with attention-deficit/hyperactivity disorder playing a real-time strategy game, we explore the relationships and temporal lag effects between pupil dilation, skin temperature, and heart rate. Our findings highlight the interconnectedness of these physiological systems and demonstrate how different physiological states are associated with unique patterns of network interactions. The study employs the concept of Time Delay Stability towards a deeper understanding of the complex dynamics involved. This research contributes to the growing field of Network Physiology by offering new insights into the physiological underpinnings of competitive gaming, potentially informing targeted training and recovery protocols for eSports athletes.

Epilepsy is a common neurological disorder, affecting over 65 million people worldwide. Unfortunately, despite resective surgery, over 30 $\%$ of patients with drug-resistant epilepsy continue to experience seizures. Retrospective studies considering connectivity using intracranial electrocorticography (ECoG) obtained during neuromonitoring have shown that treatment failure is likely driven by failure to consider critical components of the seizure network, an idea first formally introduced in 2002. However, current studies only capture snapshots in time, precluding the ability to consider seizure network development. Over the past few years, multiwell microelectrode arrays have been increasingly used to study neuronal networks in vitro. As such, we sought to develop a novel in vitro MEA seizure model to allow for study of seizure networks. Specifically, we used 4-aminopyridine (4-AP) to capture hyperexcitable activity, and then show increased network changes after 2 days of chronic treatment. We characterize network changes using functional connectivity measures and a novel technique using dimensionality reduction. We find that 4-AP successfully captures persistently elevated mean firing rate and significant changes in underlying connectivity patterns. We believe this affords a robust in vitro seizure model from which longitudinal network changes can be studied, laying groundwork for future studies exploring seizure network development.

It is increasingly understood that the epilepsies are characterized by network pathology that can span multiple spatial and temporal scales. Recent work indicates that infraslow (<0.2 Hz) envelope correlations may form a basis for distant spatial coupling in the brain. We speculated that infraslow correlation structure may be preserved even with some time lag between signals. To this end, we studied intracranial EEG (icEEG) data collected from 22 medically refractory epilepsy patients. For each patient, we selected hour-long background, awake icEEG epochs before and after antiseizure medication (ASM) taper. For each epoch, we selected 5,000 random electrode contact pairs and estimated magnitude-squared coherence (MSC) below 0.15 Hz of band power time-series in the traditional EEG frequency bands. Using these same contact pairs, we shifted one signal of the pair by random durations in 15-s increments between 0 and 300 s. We aggregated these data across all patients to determine how infraslow MSC varies with duration of lag. We further examined the effect of ASM taper on infraslow correlation structure. We also used surrogate data to empirically characterize MSC estimator and to set optimal parameters for estimation specifically for the study of infraslow activity. Our empirical analysis of the MSC estimator showed that hour-long segments with MSC computed using 3-min windows with 50% overlap was sufficient to capture infraslow envelope correlations while minimizing estimator bias and variance. The mean MSC decreased monotonically with increasing time lag until 105 s of lag, then plateaued between 106 and 300 s. Significantly nonzero infraslow envelope MSC was preserved in all frequency bands until about 1 min of time lag, both pre- and post-ASM taper. We also saw a slight, but significant increase in infraslow MSC post-ASM taper, consistent with prior work. These results provide evidence for the feasibility of examining infraslow activity via its modulation of higher-frequency activity in the absence of DC-coupled recordings. The use of surrogate data also provides a general methodology for benchmarking measures used in network neuroscience studies. Finally, our study points to the clinical relevance of infraslow activity in assessing seizure risk.

Epilepsy is a complex, multifaceted disease that affects patients in several ways in addition to seizures, including psychological, social, and quality of life issues, but epilepsy is also known to interact with sleep. Seizures often occur at the boundary between sleep and wake, patients with epilepsy often experience disrupted sleep, and the rate of inter-ictal epileptiform discharges increases during non-REM sleep. The Network Theory of Epilepsy did not address a role for sleep, but recent emphasis on the interaction between epilepsy and sleep suggests that post-seizure sleep may also be involved in the process by which seizures arise and become more severe with time ("epileptogenesis") by co-opting processes related to the formation of long-term memories. While it is generally acknowledged that recurrent seizures arise from the aberrant function of neural circuits, it is possible that the progression of epilepsy is aided by normal, physiological function of neural circuits during sleep that are driven by pathological signals. Studies recording multiple, single neurons prior to spontaneous seizures have shown that neural assemblies activated prior to the start of seizures were reactivated during post-seizure sleep, similar to the reactivation of behavioral neural assemblies, which is thought to be involved in the formation of long-term memories, a process known as Memory Consolidation. The reactivation of seizure-related neural assemblies during sleep was thus described as being a component of Seizure-Related Consolidation (SRC). These results further suggest that SRC may viewed as a network-related aspect of epilepsy, even in those seizures that have anatomically restricted neuroanatomical origins. As suggested by the Network Theory of Epilepsy as a means of interfering with ictogenesis, therapies that interfered with SRC may provide some anti-epileptogenic therapeutic benefit, even if the interference targeted structures that were not involved originally in the seizure. Here, we show how the Network Theory of Epilepsy can be expanded to include neural plasticity mechanisms associated with learning by providing an overview of Memory Consolidation, the mechanisms thought to underlie MC, their relation to Seizure-Related Consolidation, and suggesting novel, anti-epileptogenic therapies targeting interference with network activation in epilepsy following seizures during post-seizure sleep.

Oscillatory complex networks in the metastable regime have been used to study the emergence of integrated and segregated activity in the brain, which are hypothesised to be fundamental for cognition. Yet, the parameters and the underlying mechanisms necessary to achieve the metastable regime are hard to identify, often relying on maximising the correlation with empirical functional connectivity dynamics. Here, we propose and show that the brain's hierarchically modular mesoscale structure alone can give rise to robust metastable dynamics and (metastable) chimera states in the presence of phase frustration. We construct unweighted 3-layer hierarchical networks of identical Kuramoto-Sakaguchi oscillators, parameterized by the average degree of the network and a structural parameter determining the ratio of connections between and within blocks in the upper two layers. Together, these parameters affect the characteristic timescales of the system. Away from the critical synchronization point, we detect the emergence of metastable states in the lowest hierarchical layer coexisting with chimera and metastable states in the upper layers. Using the Laplacian renormalization group flow approach, we uncover two distinct pathways towards achieving the metastable regimes detected in these distinct layers. In the upper layers, we show how the symmetry-breaking states depend on the slow eigenmodes of the system. In the lowest layer instead, metastable dynamics can be achieved as the separation of timescales between layers reaches a critical threshold. Our results show an explicit relationship between metastability, chimera states, and the eigenmodes of the system, bridging the gap between harmonic based studies of empirical data and oscillatory models.