Reliability of energy landscape analysis of resting-state functional MRI data.

Energy landscape analysis is a data-driven method to analyse multidimensional time series, including functional magnetic resonance imaging (fMRI) data. It has been shown to be a useful characterization of fMRI data in health and disease. It fits an Ising model to the data and captures the dynamics of the data as movement of a noisy ball constrained on the energy landscape derived from the estimated Ising model. In the present study, we examine test-retest reliability of the energy landscape analysis. To this end, we construct a permutation test that assesses whether or not indices characterizing the energy landscape are more consistent across different sets of scanning sessions from the same participant (i.e. within-participant reliability) than across different sets of sessions from different participants (i.e. between-participant reliability). We show that the energy landscape analysis has significantly higher within-participant than between-participant test-retest reliability with respect to four commonly used indices. We also show that a variational Bayesian method, which enables us to estimate energy landscapes tailored to each participant, displays comparable test-retest reliability to that using the conventional likelihood maximization method. The proposed methodology paves the way to perform individual-level energy landscape analysis for given data sets with a statistically controlled reliability.

State-transition dynamics of resting-state functional magnetic resonance imaging data: model comparison and test-to-retest analysis.

Electroencephalogram (EEG) microstate analysis entails finding dynamics of quasi-stable and generally recurrent discrete states in multichannel EEG time series data and relating properties of the estimated state-transition dynamics to observables such as cognition and behavior. While microstate analysis has been widely employed to analyze EEG data, its use remains less prevalent in functional magnetic resonance imaging (fMRI) data, largely due to the slower timescale of such data. In the present study, we extend various data clustering methods used in EEG microstate analysis to resting-state fMRI data from healthy humans to extract their state-transition dynamics. We show that the quality of clustering is on par with that for various microstate analyses of EEG data. We then develop a method for examining test-retest reliability of the discrete-state transition dynamics between fMRI sessions and show that the within-participant test-retest reliability is higher than between-participant test-retest reliability for different indices of state-transition dynamics, different networks, and different data sets. This result suggests that state-transition dynamics analysis of fMRI data could discriminate between different individuals and is a promising tool for performing fingerprinting analysis of individuals.

Fixation dynamics on hypergraphs.

Hypergraphs have been a useful tool for analyzing population dynamics such as opinion formation and the public goods game occurring in overlapping groups of individuals. In the present study, we propose and analyze evolutionary dynamics on hypergraphs, in which each node takes one of the two types of different but constant fitness values. For the corresponding dynamics on conventional networks, under the birth-death process and uniform initial conditions, most networks are known to be amplifiers of natural selection; amplifiers by definition enhance the difference in the strength of the two competing types in terms of the probability that the mutant type fixates in the population. In contrast, we provide strong computational evidence that a majority of hypergraphs are suppressors of selection under the same conditions by combining theoretical and numerical analyses. We also show that this suppressing effect is not explained by one-mode projection, which is a standard method for expressing hypergraph data as a conventional network. Our results suggest that the modeling framework for structured populations in addition to the specific network structure is an important determinant of evolutionary dynamics, paving a way to studying fixation dynamics on higher-order networks including hypergraphs.

Gene communities in co-expression networks across different tissues.

With the recent availability of tissue-specific gene expression data, e.g., provided by the GTEx Consortium, there is interest in comparing gene co-expression patterns across tissues. One promising approach to this problem is to use a multilayer network analysis framework and perform multilayer community detection. Communities in gene co-expression networks reveal groups of genes similarly expressed across individuals, potentially involved in related biological processes responding to specific environmental stimuli or sharing common regulatory variations. We construct a multilayer network in which each of the four layers is an exocrine gland tissue-specific gene co-expression network. We develop methods for multilayer community detection with correlation matrix input and an appropriate null model. Our correlation matrix input method identifies five groups of genes that are similarly co-expressed in multiple tissues (a community that spans multiple layers, which we call a generalist community) and two groups of genes that are co-expressed in just one tissue (a community that lies primarily within just one layer, which we call a specialist community). We further found gene co-expression communities where the genes physically cluster across the genome significantly more than expected by chance (on chromosomes 1 and 11). This clustering hints at underlying regulatory elements determining similar expression patterns across individuals and cell types. We suggest that KRTAP3-1, KRTAP3-3, and KRTAP3-5 share regulatory elements in skin and pancreas. Furthermore, we find that CELA3A and CELA3B share associated expression quantitative trait loci in the pancreas. The results indicate that our multilayer community detection method for correlation matrix input extracts biologically interesting communities of genes.

Evolutionary Game Dynamics with Environmental Feedback in a Network with Two Communities.

Recent developments of eco-evolutionary models have shown that evolving feedbacks between behavioral strategies and the environment of game interactions, leading to changes in the underlying payoff matrix, can impact the underlying population dynamics in various manners. We propose and analyze an eco-evolutionary game dynamics model on a network with two communities such that players interact with other players in the same community and those in the opposite community at different rates. In our model, we consider two-person matrix games with pairwise interactions occurring on individual edges and assume that the environmental state depends on edges rather than on nodes or being globally shared in the population. We analytically determine the equilibria and their stability under a symmetric population structure assumption, and we also numerically study the replicator dynamics of the general model. The model shows rich dynamical behavior, such as multiple transcritical bifurcations, multistability, and anti-synchronous oscillations. Our work offers insights into understanding how the presence of community structure impacts the eco-evolutionary dynamics within and between niches.

Associations of conservatism and jumping to conclusions biases with aberrant salience and default mode network.

AIM: While conservatism bias refers to the human need for more evidence for decision-making than rational thinking expects, the jumping to conclusions (JTC) bias refers to the need for less evidence among individuals with schizophrenia/delusion compared to healthy people. Although the hippocampus-midbrain-striatal aberrant salience system and the salience, default mode (DMN), and frontoparietal networks ("triple networks") are implicated in delusion/schizophrenia pathophysiology, the associations between conservatism/JTC and these systems/networks are unclear. METHODS: Thirty-seven patients with schizophrenia and 33 healthy controls performed the beads task, with large and small numbers of bead draws to decision (DTD) indicating conservatism and JTC, respectively. We performed independent component analysis (ICA) of resting functional magnetic resonance imaging (fMRI) data. For systems/networks above, we investigated interactions between diagnosis and DTD, and main effects of DTD. We similarly applied ICA to structural and diffusion MRI to explore the associations between DTD and gray/white matter. RESULTS: We identified a significant main effect of DTD with functional connectivity between the striatum and DMN, which was negatively correlated with delusion severity in patients, indicating that the greater the anti-correlation between these networks, the stronger the JTC and delusion. We further observed the main effects of DTD on a gray matter network resembling the DMN, and a white matter network connecting the functional and gray matter networks (all P < 0.05, family-wise error [FWE] correction). Function and gray/white matter showed no significant interactions. CONCLUSION: Our results support the novel association of conservatism and JTC biases with aberrant salience and default brain mode.

Establishing Monoclonal Gammopathy of Undetermined Significance as an Independent Pre-Disease State of Multiple Myeloma Using Raman Spectroscopy, Dynamical Network Biomarker Theory, and Energy Landscape Analysis.

Multiple myeloma (MM) is a cancer of plasma cells. Normal (NL) cells are considered to pass through a precancerous state, such as monoclonal gammopathy of undetermined significance (MGUS), before transitioning to MM. In the present study, we acquired Raman spectra at three stages-834 NL, 711 MGUS, and 970 MM spectra-and applied the dynamical network biomarker (DNB) theory to these spectra. The DNB analysis identified MGUS as the unstable pre-disease state of MM and extracted Raman shifts at 1149 and 1527-1530 cm-1 as DNB variables. The distribution of DNB scores for each patient showed a significant difference between the mean values for MGUS and MM patients. Furthermore, an energy landscape (EL) analysis showed that the NL and MM stages were likely to become stable states. Raman spectroscopy, the DNB theory, and, complementarily, the EL analysis will be applicable to the identification of the pre-disease state in clinical samples.

Similarity-Based Analysis of Allele Frequency Distribution among Multiple Populations Identifies Adaptive Genomic Structural Variants.

Structural variants have a considerable impact on human genomic diversity. However, their evolutionary history remains mostly unexplored. Here, we developed a new method to identify potentially adaptive structural variants based on a similarity-based analysis that incorporates genotype frequency data from 26 populations simultaneously. Using this method, we analyzed 57,629 structural variants and identified 576 structural variants that show unusual population differentiation. Of these putatively adaptive structural variants, we further showed that 24 variants are multiallelic and overlap with coding sequences, and 20 variants are significantly associated with GWAS traits. Closer inspection of the haplotypic variation associated with these putatively adaptive and functional structural variants reveals deviations from neutral expectations due to: 1) population differentiation of rapidly evolving multiallelic variants, 2) incomplete sweeps, and 3) recent population-specific negative selection. Overall, our study provides new methodological insights, documents hundreds of putatively adaptive variants, and introduces evolutionary models that may better explain the complex evolution of structural variants.

Perturbation theory for evolution of cooperation on networks.

Network structure is a mechanism for promoting cooperation in social dilemma games. In the present study, we explore graph surgery, i.e., to slightly perturb the given network, towards a network that better fosters cooperation. To this end, we develop a perturbation theory to assess the change in the propensity of cooperation when we add or remove a single edge to/from the given network. Our perturbation theory is for a previously proposed random-walk-based theory that provides the threshold benefit-to-cost ratio, [Formula: see text], which is the value of the benefit-to-cost ratio in the donation game above which the cooperator is more likely to fixate than in a control case, for any finite networks. We find that [Formula: see text] decreases when we remove a single edge in a majority of cases and that our perturbation theory captures at a reasonable accuracy which edge removal makes [Formula: see text] small to facilitate cooperation. In contrast, [Formula: see text] tends to increase when we add an edge, and the perturbation theory is not good at predicting the edge addition that changes [Formula: see text] by a large amount. Our perturbation theory significantly reduces the computational complexity for calculating the outcome of graph surgery.

Fixation probability in evolutionary dynamics on switching temporal networks.

Population structure has been known to substantially affect evolutionary dynamics. Networks that promote the spreading of fitter mutants are called amplifiers of selection, and those that suppress the spreading of fitter mutants are called suppressors of selection. Research in the past two decades has found various families of amplifiers while suppressors still remain somewhat elusive. It has also been discovered that most networks are amplifiers of selection under the birth-death updating combined with uniform initialization, which is a standard condition assumed widely in the literature. In the present study, we extend the birth-death processes to temporal (i.e., time-varying) networks. For the sake of tractability, we restrict ourselves to switching temporal networks, in which the network structure deterministically alternates between two static networks at constant time intervals or stochastically in a Markovian manner. We show that, in a majority of cases, switching networks are less amplifying than both of the two static networks constituting the switching networks. Furthermore, most small switching networks, i.e., networks on six nodes or less, are suppressors, which contrasts to the case of static networks.

Web-Based Social Networks of Individuals With Adverse Childhood Experiences: Quantitative Study.

BACKGROUND: Adverse childhood experiences (ACEs), which include abuse and neglect and various household challenges such as exposure to intimate partner violence and substance use in the home, can have negative impacts on the lifelong health of affected individuals. Among various strategies for mitigating the adverse effects of ACEs is to enhance connectedness and social support for those who have experienced them. However, how the social networks of those who experienced ACEs differ from the social networks of those who did not is poorly understood. OBJECTIVE: In this study, we used Reddit and Twitter data to investigate and compare social networks between individuals with and without ACE exposure. METHODS: We first used a neural network classifier to identify the presence or absence of public ACE disclosures in social media posts. We then analyzed egocentric social networks comparing individuals with self-reported ACEs with those with no reported history. RESULTS: We found that, although individuals reporting ACEs had fewer total followers in web-based social networks, they had higher reciprocity in following behavior (ie, mutual following with other users), a higher tendency to follow and be followed by other individuals with ACEs, and a higher tendency to follow back individuals with ACEs rather than individuals without ACEs. CONCLUSIONS: These results imply that individuals with ACEs may try to actively connect with others who have similar previous traumatic experiences as a positive connection and coping strategy. Supportive interpersonal connections on the web for individuals with ACEs appear to be a prevalent behavior and may be a way to enhance social connectedness and resilience in those who have experienced ACEs.

Clinical Efficacy of Pre-Hospital Electrocardiogram Transmission in Patients Undergoing Primary Percutaneous Coronary Intervention for ST-Segment Elevation Myocardial Infarction.

Rapid reperfusion by primary percutaneous coronary intervention (pPCI) is an established strategy for the treatment of patients with ST-segment elevation myocardial infarction (STEMI). Pre-hospital electrocardiogram (PH-ECG) transmission by the emergency medical services (EMS) facilitates timely reperfusion in these patients. However, evidence regarding the clinical benefits of PH-ECG in individual hospitals is limited.This retrospective, observational study investigated the clinical efficacy of PH-ECG in STEMI patients who underwent pPCI. Of a total of 382 consecutive STEMI patients, 237 were enrolled in the study and divided into 2 groups: a PH-ECG group (n = 77) and non-PH-ECG group (n = 160). Door-to-balloon time (D2BT) was significantly shorter in the PH-ECG group (66 [52-80] min), compared to the non-PH-ECG group (70 [57-88] minutes, P = 0.01). The 30-day all-cause mortality rate was 6% in the PH-ECG group, which was significantly lower than that in the non-PH-ECG group (16%) (P = 0.037, hazard ratio [HR]: 0.38, 95% CI: 0.15-0.98). This trend was particularly evident in severely ill patients when stratified by GRACE score.The use of PH-ECG improved the survival rate of STEMI patients undergoing pPCI due to the improved pre-arrival preparation based on the EMS information. Coordination between EMS and PCI-capable institutes is essential for the management of PH-ECG.

Distinct propagation patterns of right pulmonary veins through multiple epicardial connections during right atrial pacing and sinus rhythm.

A 47-year-old man with symptomatic paroxysmal atrial fibrillation (AF) underwent AF ablation. Activation maps during right atrial pacing and sinus rhythm before the ablation revealed distinctive left atrial (LA) propagations with multiple LA breakthrough sites via epicardial connections. A wide area circumferential ablation was not able to achieve a right pulmonary vein (RPV) isolation and required an inner PV ablation to isolate the RPV. Activation maps during different rhythms before the ablation may be helpful to unmask multiple epicardial connections between the RPV and right atrium.

Effect of proximal optimization technique on coronary bifurcation stent failure: Insights from the multicenter randomized PROPOT trial.

OBJECTIVE: We investigated the effect of proximal optimization technique (POT) on coronary bifurcation stent failure (BSF) in cross-over stenting by comparing with the kissing balloon technique (KBT) in a multicenter randomized PROPOT trial. BACKGROUND: POT is recommended due to increased certainty for optimal stent expansion and side branch (SB) wiring. METHODS: We randomized 120 patients treated with crossover stenting into the POT group, which was followed by SB dilation (SBD), and the KBT group. Finally, 52 and 57 patients were analyzed by optical coherence tomography before SBD and at the final procedure, respectively. Composite BSF was defined as a maximal malapposition distance of >400 µm, or malapposed and SB-jailed strut rates of >5.95% and >21.4%, respectively. RESULTS: Composite BSF before SBD in the POT and KBT groups was observed in 29% and 26% of patients, respectively. In the POT group, differences in stent volumetric index between the proximal and distal bifurcation (odds ratio [OR] 60.35, 95% confidential interval [CI] 0.13-0.93, p = 0.036) and between the proximal bifurcation and bifurcation core (OR: 3.68, 95% CI: 1.01-13.40, p = 0.048) were identified as independent risk factors. Composite BSF at final in 27% and 32%, and unplanned additional procedures in 38% and 25% were observed, respectively. Composite BSF before SBD was a risk factor for the former (OR: 6.33, 95% CI: 1.10-36.50, p = 0.039) and the latter (OR: 6.43, 95% CI: 1.25-33.10, p = 0.026) in the POT group. CONCLUSION: POT did not result in a favorable trend in BSF. Insufficient expansion of the bifurcation core after POT was associated with BSF.

Epidemic dynamics on metapopulation networks with node2vec mobility.

Metapopulation models have been a powerful tool for both theorizing and simulating epidemic dynamics. In a metapopulation model, one considers a network composed of subpopulations and their pairwise connections, and individuals are assumed to migrate from one subpopulation to another obeying a given mobility rule. While how different mobility rules affect epidemic dynamics in metapopulation models has been studied, there have been relatively few efforts on comparison of the effects of simple (i.e., unbiased) random walks and more complex mobility rules. Here we study a susceptible-infectious-susceptible (SIS) dynamics in a metapopulation model in which individuals obey a parametric second-order random-walk mobility rule called the node2vec. We map the second-order mobility rule of the node2vec to a first-order random walk in a network whose each node is a directed edge connecting a pair of subpopulations and then derive the epidemic threshold. For various networks, we find that the epidemic threshold is large (therefore, epidemic spreading tends to be suppressed) when the individuals infrequently backtrack or infrequently visit the common neighbors of the currently visited and the last visited subpopulations than when the individuals obey the simple random walk. The amount of change in the epidemic threshold induced by the node2vec mobility is in general not as large as, but is sometimes comparable with, the one induced by the change in the diffusion rate for individuals.

Modelling state-transition dynamics in resting-state brain signals by the hidden Markov and Gaussian mixture models.

Recent studies have proposed that one can summarize brain activity into dynamics among a relatively small number of hidden states and that such an approach is a promising tool for revealing brain function. Hidden Markov models (HMMs) are a prevalent approach to inferring such neural dynamics among discrete brain states. However, the impact of assuming Markovian structure in neural time series data has not been sufficiently examined. Here, to address this situation and examine the performance of the HMM, we compare the model with the Gaussian mixture model (GMM), which is with no temporal regularization and thus a statistically simpler model than the HMM, by applying both models to synthetic time series generated from empirical resting-state functional magnetic resonance imaging (fMRI) data. We compared the GMM and HMM for various sampling frequencies, lengths of recording per participant, numbers of participants and numbers of independent component signals. We find that the HMM attains a better accuracy of estimating the hidden state than the GMM in a majority of cases. However, we also find that the accuracy of the GMM is comparable to that of the HMM under the condition that the sampling frequency is reasonably low (e.g., TR = 2.88 or 3.60 s) or the data are relatively short. These results suggest that the GMM can be a viable alternative to the HMM for investigating hidden-state dynamics under this condition.

Recurrence quantification analysis of dynamic brain networks.

Evidence suggests that brain network dynamics are a key determinant of brain function and dysfunction. Here we propose a new framework to assess the dynamics of brain networks based on recurrence analysis. Our framework uses recurrence plots and recurrence quantification analysis to characterize dynamic networks. For resting-state magnetoencephalographic dynamic functional networks (dFNs), we have found that functional networks recur more quickly in people with epilepsy than in healthy controls. This suggests that recurrence of dFNs may be used as a biomarker of epilepsy. For stereo electroencephalography data, we have found that dFNs involved in epileptic seizures emerge before seizure onset, and recurrence analysis allows us to detect seizures. We further observe distinct dFNs before and after seizures, which may inform neurostimulation strategies to prevent seizures. Our framework can also be used for understanding dFNs in healthy brain function and in other neurological disorders besides epilepsy.

Division of labour promotes the spread of information in colony emigrations by the ant Temnothorax rugatulus.

The fitness of group-living animals often depends on how well members share information needed for collective decision-making. Theoretical studies have shown that collective choices can emerge in a homogeneous group of individuals following identical rules, but real animals show much evidence for heterogeneity in the degree and nature of their contribution to group decisions. In social insects, for example, the transmission and processing of information is influenced by a well-organized division of labour. Studies that accurately quantify how this behavioural heterogeneity affects the spread of information among group members are still lacking. In this paper, we look at nest choices during colony emigrations of the ant Temnothorax rugatulus and quantify the degree of behavioural heterogeneity of workers. Using clustering methods and network analysis, we identify and characterize four behavioural castes of workers-primary, secondary, passive and wandering-covering distinct roles in the spread of information during an emigration. This detailed characterization of the contribution of each worker can improve models of collective decision-making in this species and promises a deeper understanding of behavioural variation at the colony level.

Age-related changes in the ease of dynamical transitions in human brain activity.

Executive functions, a set of cognitive processes that enable flexible behavioral control, are known to decay with aging. Because such complex mental functions are considered to rely on the dynamic coordination of functionally different neural systems, the age-related decline in executive functions should be underpinned by alteration of large-scale neural dynamics. However, the effects of age on brain dynamics have not been firmly formulated. Here, we investigate such age-related changes in brain dynamics by applying "energy landscape analysis" to publicly available functional magnetic resonance imaging data from healthy younger and older human adults. We quantified the ease of dynamical transitions between different major patterns of brain activity, and estimated it for the default mode network (DMN) and the cingulo-opercular network (CON) separately. We found that the two age groups shared qualitatively the same trajectories of brain dynamics in both the DMN and CON. However, in both of networks, the ease of transitions was significantly smaller in the older than the younger group. Moreover, the ease of transitions was associated with the performance in executive function tasks in a doubly dissociated manner: for the younger adults, the ability of executive functions was mainly correlated with the ease of transitions in the CON, whereas that for the older adults was specifically associated with the ease of transitions in the DMN. These results provide direct biological evidence for age-related changes in macroscopic brain dynamics and suggest that such neural dynamics play key roles when individuals carry out cognitively demanding tasks.

Zero-determinant strategies in finitely repeated games.

Direct reciprocity is a mechanism for sustaining mutual cooperation in repeated social dilemma games, where a player would keep cooperation to avoid being retaliated by a co-player in the future. So-called zero-determinant (ZD) strategies enable a player to unilaterally set a linear relationship between the player's own payoff and the co-player's payoff regardless of the strategy of the co-player. In the present study, we analytically study zero-determinant strategies in finitely repeated (two-person) prisoner's dilemma games with a general payoff matrix. Our results are as follows. First, we present the forms of solutions that extend the known results for infinitely repeated games (with a discount factor w of unity) to the case of finitely repeated games (0 < w < 1). Second, for the three most prominent ZD strategies, the equalizers, extortioners, and generous strategies, we derive the threshold value of w above which the ZD strategies exist. Third, we show that the only strategies that enforce a linear relationship between the two players' payoffs are either the ZD strategies or unconditional strategies, where the latter independently cooperates with a fixed probability in each round of the game, proving a conjecture previously made for infinitely repeated games.