Emergence of large-scale cell death through ferroptotic trigger waves.

Large-scale cell death is commonly observed during organismal development and in human pathologies1-5. These cell death events extend over great distances to eliminate large populations of cells, raising the question of how cell death can be coordinated in space and time. One mechanism that enables long-range signal transmission is trigger waves6, but how this mechanism might be used for death events in cell populations remains unclear. Here we demonstrate that ferroptosis, an iron- and lipid-peroxidation-dependent form of cell death, can propagate across human cells over long distances (≥5 mm) at constant speeds (around 5.5 µm min-1) through trigger waves of reactive oxygen species (ROS). Chemical and genetic perturbations indicate a primary role of ROS feedback loops (Fenton reaction, NADPH oxidase signalling and glutathione synthesis) in controlling the progression of ferroptotic trigger waves. We show that introducing ferroptotic stress through suppression of cystine uptake activates these ROS feedback loops, converting cellular redox systems from being monostable to being bistable and thereby priming cell populations to become bistable media over which ROS propagate. Furthermore, we demonstrate that ferroptosis and its propagation accompany the massive, yet spatially restricted, cell death events during muscle remodelling of the embryonic avian limb, substantiating its use as a tissue-sculpting strategy during embryogenesis. Our findings highlight the role of ferroptosis in coordinating global cell death events, providing a paradigm for investigating large-scale cell death in embryonic development and human pathologies.

Functional buffering via cell-specific gene expression promotes tissue homeostasis and cancer robustness.

Functional buffering that ensures biological robustness is critical for maintaining tissue homeostasis, organismal survival, and evolution of novelty. However, the mechanism underlying functional buffering, particularly in multicellular organisms, remains largely elusive. Here, we proposed that functional buffering can be mediated via expression of buffering genes in specific cells and tissues, by which we named Cell-specific Expression-BUffering (CEBU). We developed an inference index (C-score) for CEBU by computing C-scores across 684 human cell lines using genome-wide CRISPR screens and transcriptomic RNA-seq. We report that C-score-identified putative buffering gene pairs are enriched for members of the same duplicated gene family, pathway, and protein complex. Furthermore, CEBU is especially prevalent in tissues of low regenerative capacity (e.g., bone and neuronal tissues) and is weakest in highly regenerative blood cells, linking functional buffering to tissue regeneration. Clinically, the buffering capacity enabled by CEBU can help predict patient survival for multiple cancers. Our results suggest CEBU as a potential buffering mechanism contributing to tissue homeostasis and cancer robustness in humans.

Multistability maintains redox homeostasis in human cells.

Cells metabolize nutrients through a complex metabolic and signaling network that governs redox homeostasis. At the core of this, redox regulatory network is a mutually inhibitory relationship between reduced glutathione and reactive oxygen species (ROS)-two opposing metabolites that are linked to upstream nutrient metabolic pathways (glucose, cysteine, and glutamine) and downstream feedback loops of signaling pathways (calcium and NADPH oxidase). We developed a nutrient-redox model of human cells to understand system-level properties of this network. Combining in silico modeling and ROS measurements in individual cells, we show that ROS dynamics follow a switch-like, all-or-none response upon glucose deprivation at a threshold that is approximately two orders of magnitude lower than its physiological concentration. We also confirm that this ROS switch can be irreversible and exhibits hysteresis, a hallmark of bistability. Our findings evidence that bistability modulates redox homeostasis in human cells and provide a general framework for quantitative investigations of redox regulation in humans.

Key Immune Events of the Pathomechanisms of Early Cardioembolic Stroke: Multi-Database Mining and Systems Biology Approach.

While inflammation has generally been regarded as a negative factor in stroke recovery, this viewpoint has recently been challenged by demonstrating that inflammation is a necessary and sufficient factor for regeneration in the zebrafish brain injury model. This close relationship with inflammation suggests that a re-examination of the immune system's role in strokes is necessary. We used a systems biology approach to investigate the role of immune-related functions via their interactions with other molecular functions in early cardioembolic stroke. Based on protein interaction models and on microarray data from the blood of stroke subjects and healthy controls, networks were constructed to delineate molecular interactions at four early stages (pre-stroke, 3 h, 5 h and 24 h after stroke onset) of cardioembolic stroke. A comparative analysis of functional networks identified interactions of immune-related functions with other molecular functions, including growth factors, neuro/hormone and housekeeping functions. These provide a potential pathomechanism for early stroke pathophysiology. In addition, several potential targets of miRNA and methylation regulations were derived based on basal level changes observed in the core networks and literature. The results provide a more comprehensive understanding of stroke progression mechanisms from an immune perspective and shed light on acute stroke treatments.

Temporal Genetic Modifications after Controlled Cortical Impact--Understanding Traumatic Brain Injury through a Systematic Network Approach.

Traumatic brain injury (TBI) is a primary injury caused by external physical force and also a secondary injury caused by biological processes such as metabolic, cellular, and other molecular events that eventually lead to brain cell death, tissue and nerve damage, and atrophy. It is a common disease process (as opposed to an event) that causes disabilities and high death rates. In order to treat all the repercussions of this injury, treatment becomes increasingly complex and difficult throughout the evolution of a TBI. Using high-throughput microarray data, we developed a systems biology approach to explore potential molecular mechanisms at four time points post-TBI (4, 8, 24, and 72 h), using a controlled cortical impact (CCI) model. We identified 27, 50, 48, and 59 significant proteins as network biomarkers at these four time points, respectively. We present their network structures to illustrate the protein-protein interactions (PPIs). We also identified UBC (Ubiquitin C), SUMO1, CDKN1A (cyclindependent kinase inhibitor 1A), and MYC as the core network biomarkers at the four time points, respectively. Using the functional analytical tool MetaCore™, we explored regulatory mechanisms and biological processes and conducted a statistical analysis of the four networks. The analytical results support some recent findings regarding TBI and provide additional guidance and directions for future research.

A Systems Biology Approach to the Coordination of Defensive and Offensive Molecular Mechanisms in the Innate and Adaptive Host-Pathogen Interaction Networks.

Infected zebrafish coordinates defensive and offensive molecular mechanisms in response to Candida albicans infections, and invasive C. albicans coordinates corresponding molecular mechanisms to interact with the host. However, knowledge of the ensuing infection-activated signaling networks in both host and pathogen and their interspecific crosstalk during the innate and adaptive phases of the infection processes remains incomplete. In the present study, dynamic network modeling, protein interaction databases, and dual transcriptome data from zebrafish and C. albicans during infection were used to infer infection-activated host-pathogen dynamic interaction networks. The consideration of host-pathogen dynamic interaction systems as innate and adaptive loops and subsequent comparisons of inferred innate and adaptive networks indicated previously unrecognized crosstalk between known pathways and suggested roles of immunological memory in the coordination of host defensive and offensive molecular mechanisms to achieve specific and powerful defense against pathogens. Moreover, pathogens enhance intraspecific crosstalk and abrogate host apoptosis to accommodate enhanced host defense mechanisms during the adaptive phase. Accordingly, links between physiological phenomena and changes in the coordination of defensive and offensive molecular mechanisms highlight the importance of host-pathogen molecular interaction networks, and consequent inferences of the host-pathogen relationship could be translated into biomedical applications.

Identification of network-based biomarkers of cardioembolic stroke using a systems biology approach with time series data.

BACKGROUND: Molecular signaling of angiogenesis begins within hours after initiation of a stroke and the following regulation of endothelial integrity mediated by growth factor receptors and vascular growth factors. Recent studies further provided insights into the coordinated patterns of post-stroke gene expressions and the relationships between neurodegenerative diseases and neural function recovery processes after a stroke. RESULTS: Differential protein-protein interaction networks (PPINs) were constructed at 3 post-stroke time points, and proteins with a significant stroke relevance value (SRV) were discovered. Genes, including UBC, CUL3, APP, NEDD8, JUP, and SIRT7, showed high associations with time after a stroke, and Ingenuity Pathway Analysis results showed that these post-stroke time series-associated genes were related to molecular and cellular functions of cell death, cell survival, the cell cycle, cellular development, cellular movement, and cell-to-cell signaling and interactions. These biomarkers may be helpful for the early detection, diagnosis, and prognosis of ischemic stroke. CONCLUSIONS: This is our first attempt to use our theory of a systems biology framework on strokes. We focused on 3 key post-stroke time points. We identified the network and corresponding network biomarkers for the 3 time points, further studies are needed to experimentally confirm the findings and compare them with the causes of ischemic stroke. Our findings showed that stroke-associated biomarker genes at different time points were significantly involved in cell cycle processing, including G2-M, G1-S and meiosis, which contributes to the current understanding of the etiology of stroke. We hope this work helps scientists reveal more hidden cellular mechanisms of stroke etiology and repair processes.

Applying NGS Data to Find Evolutionary Network Biomarkers from the Early and Late Stages of Hepatocellular Carcinoma.

Hepatocellular carcinoma (HCC) is a major liver tumor (~80%), besides hepatoblastomas, angiosarcomas, and cholangiocarcinomas. In this study, we used a systems biology approach to construct protein-protein interaction networks (PPINs) for early-stage and late-stage liver cancer. By comparing the networks of these two stages, we found that the two networks showed some common mechanisms and some significantly different mechanisms. To obtain differential network structures between cancer and noncancer PPINs, we constructed cancer PPIN and noncancer PPIN network structures for the two stages of liver cancer by systems biology method using NGS data from cancer cells and adjacent noncancer cells. Using carcinogenesis relevance values (CRVs), we identified 43 and 80 significant proteins and their PPINs (network markers) for early-stage and late-stage liver cancer. To investigate the evolution of network biomarkers in the carcinogenesis process, a primary pathway analysis showed that common pathways of the early and late stages were those related to ordinary cancer mechanisms. A pathway specific to the early stage was the mismatch repair pathway, while pathways specific to the late stage were the spliceosome pathway, lysine degradation pathway, and progesterone-mediated oocyte maturation pathway. This study provides a new direction for cancer-targeted therapies at different stages.

Dynamic network-based relevance score reveals essential proteins and functional modules in directed differentiation.

The induction of stem cells toward a desired differentiation direction is required for the advancement of stem cell-based therapies. Despite successful demonstrations of the control of differentiation direction, the effective use of stem cell-based therapies suffers from a lack of systematic knowledge regarding the mechanisms underlying directed differentiation. Using dynamic modeling and the temporal microarray data of three differentiation stages, three dynamic protein-protein interaction networks were constructed. The interaction difference networks derived from the constructed networks systematically delineated the evolution of interaction variations and the underlying mechanisms. A proposed relevance score identified the essential components in the directed differentiation. Inspection of well-known proteins and functional modules in the directed differentiation showed the plausibility of the proposed relevance score, with the higher scores of several proteins and function modules indicating their essential roles in the directed differentiation. During the differentiation process, the proteins and functional modules with higher relevance scores also became more specific to the neuronal identity. Ultimately, the essential components revealed by the relevance scores may play a role in controlling the direction of differentiation. In addition, these components may serve as a starting point for understanding the systematic mechanisms of directed differentiation and for increasing the efficiency of stem cell-based therapies.

A systems biology approach to study systemic inflammation.

Systemic inflammation needs a precise control on the sequence and magnitude of occurring events. The high throughput data on the host-pathogen interactions gives us an opportunity to have a glimpse on the systemic inflammation. In this article, a dynamic Candida albicans-zebrafish interactive infectious network is built as an example to demonstrate how systems biology approach can be used to study systematic inflammation. In particular, based on microarray data of C. albicans and zebrafish during infection, the hyphal growth, zebrafish, and host-pathogen intercellular PPI networks were combined to form an integrated infectious PPI network that helps us understand the systematic mechanisms underlying the pathogenicity of C. albicans and the immune response of the host. The signaling pathways for morphogenesis and hyphal growth of C. albicans were 2 significant interactions found in the intercellular PPI network. Two cellular networks were also developed corresponding to the different infection stages (adhesion and invasion), and then compared with each other to identify proteins to gain more insight into the pathogenic role of hyphal development in the C. albicans infection process. Important defense-related proteins in zebrafish were predicted using the same approach. This integrated network consisting of intercellular invasion and cellular defense processes during infection can improve medical therapies and facilitate development of new antifungal drugs.

On the crucial cerebellar wound healing-related pathways and their cross-talks after traumatic brain injury in Danio rerio.

Upon injury, the direct damage and the subsequent secondary injury in the brain often result in chronic neurological disorders. Due to multifactorial nature of secondary injury and subsequent complex cellular responses, much of the underlying mechanisms are unclear. This study used an adult zebrafish cerebellum injury model to investigate the phenotypes and the secondary injury responses for recovery mechanisms of injured brain. Using the time course microarray analysis, a candidate protein-protein interaction (PPI) network was refined as cerebellar wound healing PPI network by dynamic modeling and big data mining. Pathway enrichment and ontological analysis were incorporated into the refined network to highlight the main molecular scheme of cerebellar wound healing. Several significant pathways, including chemokine, Phosphatidylinositide 3-kinases, and axon guidance signaling pathway and their cross-talks through PI3K, PAK2, and PLXNA3 were identified to coordinate for neurogenesis and angiogenesis, which are essential for the restoration of the injured brain. Our finding provides an insight into the molecular restoration mechanisms after traumatic brain injury, and open up new opportunity to devise the treatment for traumatic brain injury in human.

Systems biology as an integrated platform for bioinformatics, systems synthetic biology, and systems metabolic engineering.

Systems biology aims at achieving a system-level understanding of living organisms and applying this knowledge to various fields such as synthetic biology, metabolic engineering, and medicine. System-level understanding of living organisms can be derived from insight into: (i) system structure and the mechanism of biological networks such as gene regulation, protein interactions, signaling, and metabolic pathways; (ii) system dynamics of biological networks, which provides an understanding of stability, robustness, and transduction ability through system identification, and through system analysis methods; (iii) system control methods at different levels of biological networks, which provide an understanding of systematic mechanisms to robustly control system states, minimize malfunctions, and provide potential therapeutic targets in disease treatment; (iv) systematic design methods for the modification and construction of biological networks with desired behaviors, which provide system design principles and system simulations for synthetic biology designs and systems metabolic engineering. This review describes current developments in systems biology, systems synthetic biology, and systems metabolic engineering for engineering and biology researchers. We also discuss challenges and future prospects for systems biology and the concept of systems biology as an integrated platform for bioinformatics, systems synthetic biology, and systems metabolic engineering.

On the calculation of signal transduction ability of signaling transduction pathways in intracellular communication: systematic approach.

MOTIVATION: The major function of signal transduction pathways in cells is to sense signals from the environment and process the information through signaling molecules in order to regulate the activity of transcription factors. On the molecular level, the information transmitted by a small number of signal molecules is amplified in the internal signaling pathway through enzyme catalysis, molecular modification and via the activation or inhibition of interactions. However, the dynamic system behavior of a signaling pathway can be complex and, despite knowledge of the pathway components and interactions, it is still a challenge to interpret the pathways behavior. Therefore, a systematic method is proposed in this study to quantify the signal transduction ability. RESULTS: Based on the non-linear signal transduction system, signal transduction ability can be investigated by solving a Hamilton-Jacobi inequality (HJI)-constrained optimization problem. To avoid difficulties associated with solving a complex HJI-constrained optimization problem for signal transduction ability, the Takagi-Sugeno fuzzy model is introduced to approximate the non-linear signal transduction system by interpolating several local linear systems so that the HJI-constrained optimization problem can be replaced by a linear matrix inequality (LMI)-constrained optimization problem. The LMI problem can then be efficiently solved for measuring signal transduction ability. Finally, the signal transduction ability of two important signal transduction pathways was measured by the proposed method and confirmed using experimental data, which is useful for biotechnological and therapeutic application and drug design.

Three-dimensional reconstruction of brain-wide wiring networks in Drosophila at single-cell resolution.

BACKGROUND: Animal behavior is governed by the activity of interconnected brain circuits. Comprehensive brain wiring maps are thus needed in order to formulate hypotheses about information flow and also to guide genetic manipulations aimed at understanding how genes and circuits orchestrate complex behaviors. RESULTS: To assemble this map, we deconstructed the adult Drosophila brain into approximately 16,000 single neurons and reconstructed them into a common standardized framework to produce a virtual fly brain. We have constructed a mesoscopic map and found that it consists of 41 local processing units (LPUs), six hubs, and 58 tracts covering the whole Drosophila brain. Despite individual local variation, the architecture of the Drosophila brain shows invariance for both the aggregation of local neurons (LNs) within specific LPUs and for the connectivity of projection neurons (PNs) between the same set of LPUs. An open-access image database, named FlyCircuit, has been constructed for online data archiving, mining, analysis, and three-dimensional visualization of all single neurons, brain-wide LPUs, their wiring diagrams, and neural tracts. CONCLUSION: We found that the Drosophila brain is assembled from families of multiple LPUs and their interconnections. This provides an essential first step in the analysis of information processing within and between neurons in a complete brain.