Identification of feedback loops in neural networks based on multi-step Granger causality

LncRNAs driving feedback loops to boost drug resistance: sinuous pathways in cancer

Feedback loops mediate signaling pathways to maintain cellular homeostasis. There are two types, positive and negative feedback loops. Both are subject to alterations, and consequently can become pathogenic in the development of diseases such as cancer. Long noncoding RNAs (lncRNAs) are regulators of signaling pathways through feedback loops hidden as the dark regulatory elements yet to be described with great impact on cancer tumorigenesis, development, and drug resistance. Several feedback loops have been studied in cancer, however, how they are regulated by lncRNAs is hardly evident, setting a trending topic in oncological research. In this review, we recapitulate and discuss the feedback loops that are regulated by lncRNAs to promote drug resistance. Furthermore, we propose additional strategies that allow us to identify, analyze and comprehend feedback loops regulated by lncRNAs to induce drug resistance or even to gain insight into novel feedback loops that are stimulated under the pressure of treatment and consequently increase its efficacy. This knowledge will be useful to optimize the therapeutic use of oncological drugs.

Prediction of blood test values under different lifestyle scenarios using time-series electronic health record

Owing to increasing medical expenses, researchers have attempted to detect clinical signs and preventive measures of diseases using electronic health record (EHR). In particular, time-series EHRs collected by periodic medical check-up enable us to clarify the relevance among check-up results and individual environmental factors such as lifestyle. However, usually such time-series data have many missing observations and some results are strongly correlated to each other. These problems make the analysis difficult and there exists strong demand to detect clinical findings beyond them. We focus on blood test values in medical check-up results and apply a time-series analysis methodology using a state space model. It can infer the internal medical states emerged in blood test values and handle missing observations. The estimated models enable us to predict one's blood test values under specified condition and predict the effect of intervention, such as changes of body composition and lifestyle. We use time-series data of EHRs periodically collected in the Hirosaki cohort study in Japan and elucidate the effect of 17 environmental factors to 38 blood test values in elderly people. Using the estimated model, we then simulate and compare time-transitions of participant's blood test values under several lifestyle scenarios. It visualizes the impact of lifestyle changes for the prevention of diseases. Finally, we exemplify that prediction errors under participant's actual lifestyle can be partially explained by genetic variations, and some of their effects have not been investigated by traditional association studies.

Identifying the pulsed neuron networks' structures by a nonlinear Granger causality method

Background

It is a crucial task of brain science researches to explore functional connective maps of Biological Neural Networks (BNN). The maps help to deeply study the dominant relationship between the structures of the BNNs and their network functions.

Results

In this study, the ideas of linear Granger causality modeling and causality identification are extended to those of nonlinear Granger causality modeling and network structure identification. We employed Radial Basis Functions to fit the nonlinear multivariate dynamical responses of BNNs with neuronal pulse firing. By introducing the contributions from presynaptic neurons and detecting whether the predictions for postsynaptic neurons’ pulse firing signals are improved or not, we can reveal the information flows distribution of BNNs. Thus, the functional connections from presynaptic neurons can be identified from the obtained network information flows. To verify the effectiveness of the proposed method, the Nonlinear Granger Causality Identification Method (NGCIM) is applied to the network structure discovery processes of Spiking Neural Networks (SNN). SNN is a simulation model based on an Integrate-and-Fire mechanism. By network simulations, the multi-channel neuronal pulse sequence data of the SNNs can be used to reversely identify the synaptic connections and strengths of the SNNs.

Conclusions

The identification results show: for 2–6 nodes small-scale neural networks, 20 nodes medium-scale neural networks, and 100 nodes large-scale neural networks, the identification accuracy of NGCIM with the Gaussian kernel function was 100%, 99.64%, 98.64%, 98.37%, 98.31%, 84.87% and 80.56%, respectively. The identification accuracies were significantly higher than those of a traditional Linear Granger Causality Identification Method with the same network sizes. Thus, with an accumulation of the data obtained by the existing measurement methods, such as Electroencephalography, functional Magnetic Resonance Imaging, and Multi-Electrode Array, the NGCIM can be a promising network modeling method to infer the functional connective maps of BNNs.

Slow-Wave Recordings From Micro-Sized Neural Clusters Using Multiwell Type Microelectrode Arrays

OBJECTIVE

The use of microelectrode array (MEA) recordings is a very effective neurophysiological method because it is able to continuously and noninvasively obtain the spatiotemporal information of electrical activity from many neurons constituting a neural network. Very recently, studies have been published that used MEAs for the measurement of a low-frequency component of electrical activity as an indicator of diverse activity of cultured neurons. The occurrence of low-frequency activities has electrophysiological information that does not include the information from fast spikes. However, there is no in vitro experimental model suitable for measuring the low-frequency activities (slow-waves) for further study.

METHODS

Neural clusters consisting of dozens of neurons were placed directly onto each electrode of an MEA from which fast spikes and slow-waves were measured.

RESULTS

We obtained sufficient data on the early development patterns of the slow-waves and the spikes measured from many independent neural clusters confirming that the slow-waves occurred first before the emergence of the spikes in the neural clusters. We also showed that changes in the occurrence frequency of the slow-waves for synaptic blockers were measured from a large number of independent cultures.

CONCLUSION

Microsized neural cluster arrays, which can be combined with conventional MEAs, are suitable for multiple simultaneous recordings of slow-waves.

SIGNIFICANCE

Our technology provides a simple but useful method to study the generation of a low-frequency component of the electrical activity in cultured neural networks that are not yet well known as well as to expand the use of conventional MEAs.

Inferring Broad Regulatory Biology from Time Course Data: Have We Reached an Upper Bound under Constraints Typical of In Vivo Studies?

There is a growing appreciation for the network biology that regulates the coordinated expression of molecular and cellular markers however questions persist regarding the identifiability of these networks. Here we explore some of the issues relevant to recovering directed regulatory networks from time course data collected under experimental constraints typical of in vivo studies. NetSim simulations of sparsely connected biological networks were used to evaluate two simple feature selection techniques used in the construction of linear Ordinary Differential Equation (ODE) models, namely truncation of terms versus latent vector projection. Performance was compared with ODE-based Time Series Network Identification (TSNI) integral, and the information-theoretic Time-Delay ARACNE (TD-ARACNE). Projection-based techniques and TSNI integral outperformed truncation-based selection and TD-ARACNE on aggregate networks with edge densities of 10-30%, i.e. transcription factor, protein-protein cliques and immune signaling networks. All were more robust to noise than truncation-based feature selection. Performance was comparable on the in silico 10-node DREAM 3 network, a 5-node Yeast synthetic network designed for In vivo Reverse-engineering and Modeling Assessment (IRMA) and a 9-node human HeLa cell cycle network of similar size and edge density. Performance was more sensitive to the number of time courses than to sample frequency and extrapolated better to larger networks by grouping experiments. In all cases performance declined rapidly in larger networks with lower edge density. Limited recovery and high false positive rates obtained overall bring into question our ability to generate informative time course data rather than the design of any particular reverse engineering algorithm.

Axonal conduction slowing induced by spontaneous bursting activity in cortical neurons cultured in a microtunnel device

We demonstrated that spontaneous bursting activity can decrease the axonal conduction velocity of cortical neurons cultured in a microtunnel device.

Short-term dynamics of causal information transfer in thalamocortical networks during natural inputs and microstimulation for somatosensory neuroprosthesis

Recording the activity of large populations of neurons requires new methods to analyze and use the large volumes of time series data thus created. Fast and clear methods for finding functional connectivity are an important step toward the goal of understanding neural processing. This problem presents itself readily in somatosensory neuroprosthesis (SSNP) research, which uses microstimulation (MiSt) to activate neural tissue to mimic natural stimuli, and has the capacity to potentiate, depotentiate, or even destroy functional connections. As the aim of SSNP engineering is artificially creating neural responses that resemble those observed during natural inputs, a central goal is describing the influence of MiSt on activity structure among groups of neurons, and how this structure may be altered to affect perception or behavior. In this paper, we demonstrate the concept of Granger causality, combined with maximum likelihood methods, applied to neural signals recorded before, during, and after natural and electrical stimulation. We show how these analyses can be used to evaluate the changing interactions in the thalamocortical somatosensory system in response to repeated perturbation. Using LFPs recorded from the ventral posterolateral thalamus (VPL) and somatosensory cortex (S1) in anesthetized rats, we estimated pair-wise functional interactions between functional microdomains. The preliminary results demonstrate input-dependent modulations in the direction and strength of information flow during and after application of MiSt. Cortico-cortical interactions during cortical MiSt and baseline conditions showed the largest causal influence differences, while there was no statistically significant difference between pre- and post-stimulation baseline causal activities. These functional connectivity changes agree with physiologically accepted communication patterns through the network, and their particular parameters have implications for both rehabilitation and brain—machine interface SSNP applications.

Inference of gene regulatory networks incorporating multi-source biological knowledge via a state space model with L1 regularization

Comprehensive understanding of gene regulatory networks (GRNs) is a major challenge in the field of systems biology. Currently, there are two main approaches in GRN analysis using time-course observation data, namely an ordinary differential equation (ODE)-based approach and a statistical model-based approach. The ODE-based approach can generate complex dynamics of GRNs according to biologically validated nonlinear models. However, it cannot be applied to ten or more genes to simultaneously estimate system dynamics and regulatory relationships due to the computational difficulties. The statistical model-based approach uses highly abstract models to simply describe biological systems and to infer relationships among several hundreds of genes from the data. However, the high abstraction generates false regulations that are not permitted biologically. Thus, when dealing with several tens of genes of which the relationships are partially known, a method that can infer regulatory relationships based on a model with low abstraction and that can emulate the dynamics of ODE-based models while incorporating prior knowledge is urgently required. To accomplish this, we propose a method for inference of GRNs using a state space representation of a vector auto-regressive (VAR) model with L1 regularization. This method can estimate the dynamic behavior of genes based on linear time-series modeling constructed from an ODE-based model and can infer the regulatory structure among several tens of genes maximizing prediction ability for the observational data. Furthermore, the method is capable of incorporating various types of existing biological knowledge, e.g., drug kinetics and literature-recorded pathways. The effectiveness of the proposed method is shown through a comparison of simulation studies with several previous methods. For an application example, we evaluated mRNA expression profiles over time upon corticosteroid stimulation in rats, thus incorporating corticosteroid kinetics/dynamics, literature-recorded pathways and transcription factor (TF) information.

Heat shock protein 70 (Hsp70) inhibits oxidative phosphorylation and compensates ATP balance through enhanced glycolytic activity

To address possible effects of heat shock protein 70 (Hsp70) on energy metabolism, we established a cell line expressing different levels of Hsp70 and evaluated changes in glucose and lactate metabolites, as well as ATP levels accordingly. In addition, activities of enzymes involved in glycolysis [phosphofructokinase (PFK) and lactate dehydrogenase (LDH)], Krebs cycle [citric synthase (CS)], and oxidative phosphorylation {NADH dehydrogenase [complex I (CI)] and ubiquinol:cytochrome-c reductase [complex III (CIII)]} were analyzed. The results show that both glucose consumption and lactate excretion were elevated significantly in cells expressing increased levels of Hsp70. Simultaneously, the activities of glycolytic enzymes PFK and LDH were increased markedly in cells overexpressing Hsp70. Activities of enzymes CI and CIII, both involved in oxidative phosphorylation, decreased upon increased expression of Hsp70. These findings were supported by nonsignificant reductions of CS activities in cells that overexpressed Hsp70, whereas intracellular ATP levels remained constant over a wide range of Hsp70 expression. In conclusion, overexpression of Hsp70 in HeLa cells results in downregulation of oxidative phosphorylation, in particular, multiprotein CIII, the main source of reactive oxygen species. In exchange, upregulation of the glycolytic pathway compensates for the homeostasis of cellular ATP supply.