Frontiers in network physiology最新文献

In this mini review, we propose the use of the Julia programming language and its software as a strong candidate for reproducible, efficient, and sustainable physiological signal analysis. First, we highlight available software and Julia communities that provide top-of-the-class algorithms for all aspects of physiological signal processing despite the language's relatively young age. Julia can significantly accelerate both research and software development due to its high-level interactive language and high-performance code generation. It is also particularly suited for open and reproducible science. Openness is supported and welcomed because the overwhelming majority of Julia software programs are open source and developed openly on public platforms, primarily through individual contributions. Such an environment increases the likelihood that an individual not (originally) associated with a software program would still be willing to contribute their code, further promoting code sharing and reuse. On the other hand, Julia's exceptionally strong package manager and surrounding ecosystem make it easy to create self-contained, reproducible projects that can be instantly installed and run, irrespective of processor architecture or operating system.

The nervous system, especially the human brain, is characterized by its highly complex network topology. The neurodevelopment of some of its features has been described in terms of dynamic optimization rules. We discuss the principle of adaptive rewiring, i.e., the dynamic reorganization of a network according to the intensity of internal signal communication as measured by synchronization or diffusion, and its recent generalization for applications in directed networks. These have extended the principle of adaptive rewiring from highly oversimplified networks to more neurally plausible ones. Adaptive rewiring captures all the key features of the complex brain topology: it transforms initially random or regular networks into networks with a modular small-world structure and a rich-club core. This effect is specific in the sense that it can be tailored to computational needs, robust in the sense that it does not depend on a critical regime, and flexible in the sense that parametric variation generates a range of variant network configurations. Extreme variant networks can be associated at macroscopic level with disorders such as schizophrenia, autism, and dyslexia, and suggest a relationship between dyslexia and creativity. Adaptive rewiring cooperates with network growth and interacts constructively with spatial organization principles in the formation of topographically distinct modules and structures such as ganglia and chains. At the mesoscopic level, adaptive rewiring enables the development of functional architectures, such as convergent-divergent units, and sheds light on the early development of divergence and convergence in, for example, the visual system. Finally, we discuss future prospects for the principle of adaptive rewiring.

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.