The emergence of modularity in biological systems

Modular Control of Biological Networks

The concept of control is central to understanding and applications of biological network models. Some of their key structural features relate to control functions, through gene regulation, signaling, or metabolic mechanisms, and computational models need to encode these. Applications of models often focus on model-based control, such as in biomedicine or metabolic engineering. This paper presents an approach to model-based control that exploits two common features of biological networks, namely their modular structure and canalizing features of their regulatory mechanisms. The paper focuses on intracellular regulatory networks, represented by Boolean network models. A main result of this paper is that control strategies can be identified by focusing on one module at a time. This paper also presents a criterion based on canalizing features of the regulatory rules to identify modules that do not contribute to network control and can be excluded. For even moderately sized networks, finding global control inputs is computationally very challenging. The modular approach presented here leads to a highly efficient approach to solving this problem. This approach is applied to a published Boolean network model of blood cancer large granular lymphocyte (T-LGL) leukemia to identify a minimal control set that achieves a desired control objective.

A novel bioinformatics strategy to uncover the active ingredients and molecular mechanisms of Bai Shao in the treatment of non-alcoholic fatty liver disease

Introduction: As a new discipline, network pharmacology has been widely used to disclose the material basis and mechanism of Traditional Chinese Medicine in recent years. However, numerous researches indicated that the material basis of TCMs identified based on network pharmacology was the mixtures of beneficial and harmful substances rather than the real material basis. In this work, taking the anti-NAFLD (non-alcoholic fatty liver disease) effect of Bai Shao (BS) as a case, we attempted to propose a novel bioinformatics strategy to uncover the material basis and mechanism of TCMs in a precise manner. Methods: In our previous studies, we have done a lot work to explore TCM-induced hepatoprotection. Here, by integrating our previous studies, we developed a novel computational pharmacology method to identify hepatoprotective ingredients from TCMs. Then the developed method was used to discover the material basis and mechanism of Bai Shao against Non-alcoholic fatty liver disease by combining with the techniques of molecular network, microarray data analysis, molecular docking, and molecular dynamics simulation. Finally, literature verification method was utilized to validate the findings. Results: A total of 12 ingredients were found to be associated with the anti-NAFLD effect of BS, including monoterpene glucosides, flavonoids, triterpenes, and phenolic acids. Further analysis found that IL1-β, IL6, and JUN would be the key targets. Interestingly, molecular docking and molecular dynamics simulation analysis showed that there indeed existed strong and stable binding affinity between the active ingredients and the key targets. In addition, a total of 23 NAFLD-related KEGG pathways were enriched. The major biological processes involved by these pathways including inflammation, apoptosis, lipid metabolism, and glucose metabolism. Of note, there was a great deal of evidence available in the literature to support the findings mentioned above, indicating that our method was reliable. Discussion: In summary, the contributions of this work can be summarized as two aspects as follows. Firstly, we systematically elucidated the material basis and mechanism of BS against NAFLD from multiple perspectives. These findings further enhanced the theoretical foundation of BS on NAFLD. Secondly, a novel computational pharmacology research strategy was proposed, which would assist network pharmacology to uncover the scientific connotation TCMs in a more precise manner.

Community detection in the human connectome: Method types, differences and their impact on inference

Community structure is a fundamental topological characteristic of optimally organized brain networks. Currently, there is no clear standard or systematic approach for selecting the most appropriate community detection method. Furthermore, the impact of method choice on the accuracy and robustness of estimated communities (and network modularity), as well as method-dependent relationships between network communities and cognitive and other individual measures, are not well understood. This study analyzed large datasets of real brain networks (estimated from resting-state fMRI fromn $$ n $$ = 5251 pre/early adolescents in the adolescent brain cognitive development [ABCD] study), andn $$ n $$ = 5338 synthetic networks with heterogeneous, data-inspired topologies, with the goal to investigate and compare three classes of community detection methods: (i) modularity maximization-based (Newman and Louvain), (ii) probabilistic (Bayesian inference within the framework of stochastic block modeling (SBM)), and (iii) geometric (based on graph Ricci flow). Extensive comparisons between methods and their individual accuracy (relative to the ground truth in synthetic networks), and reliability (when applied to multiple fMRI runs from the same brains) suggest that the underlying brain network topology plays a critical role in the accuracy, reliability and agreement of community detection methods. Consistent method (dis)similarities, and their correlations with topological properties, were estimated across fMRI runs. Based on synthetic graphs, most methods performed similarly and had comparable high accuracy only in some topological regimes, specifically those corresponding to developed connectomes with at least quasi-optimal community organization. In contrast, in densely and/or weakly connected networks with difficult to detect communities, the methods yielded highly dissimilar results, with Bayesian inference within SBM having significantly higher accuracy compared to all others. Associations between method-specific modularity and demographic, anthropometric, physiological and cognitive parameters showed mostly method invariance but some method dependence as well. Although method sensitivity to different levels of community structure may in part explain method-dependent associations between modularity estimates and parameters of interest, method dependence also highlights potential issues of reliability and reproducibility. These findings suggest that a probabilistic approach, such as Bayesian inference in the framework of SBM, may provide consistently reliable estimates of community structure across network topologies. In addition, to maximize robustness of biological inferences, identified network communities and their cognitive, behavioral and other correlates should be confirmed with multiple reliable detection methods.

Identifying Qingkailing ingredients-dependent mesenchymal-epithelial transition factor-axiation "π" structuring module with angiogenesis and neurogenesis effects

OBJECTIVE

To explore the functional role of the drug-dependent mesenchymal-epithelial transition (Met)-axiation "π" structural module of neurogenesis after processing by three components of Qingkailing injection in neurogenesis and angiogenesis in cerebral ischemia.

METHODS

We used a Glutathione S-transferase (GST)-pull down assay, isothermal titration calorimetry assay, and other related methods to identify the relationships among Met, inositol polyphosphate phosphatase like 1 (Inppl1), and death associated protein kinase 3 (Dapk3) in this allosteric module. The biological effects of the modules of neurons generation composed of Met, Inppl1, and Dapk3 were measured through Western blot, apoptosis analysis, and double immunofluorescence labeling.

RESULTS

The GST-pull down assay revealed that proline-serine-threonine rich domain of Met binds to the Src homology domain of Inppl1 to form a protein-protein complex; Dapk3 with a C-terminal domain interacts weakly with the protein kinase C domain of Met in the intracellular region. Thus, we obtained a "π" structuring module considered a neural regeneration module. The biological effects of angiogenesis and neurogenesis modules composed of Met, Inppl1, and Dapk3 were also verified.

CONCLUSION

The study suggested that understanding the functional modules that contribute to pharmaceutics might provide novel signatures that can be used as endpoints to define disease processes under stroke or cerebral ischemia conditions.

Applying precision medicine principles to the management of multimorbidity: the utility of comorbidity networks, graph machine learning, and knowledge graphs

The current management of patients with multimorbidity is suboptimal, with either a single-disease approach to care or treatment guideline adaptations that result in poor adherence due to their complexity. Although this has resulted in calls for more holistic and personalized approaches to prescribing, progress toward these goals has remained slow. With the rapid advancement of machine learning (ML) methods, promising approaches now also exist to accelerate the advance of precision medicine in multimorbidity. These include analyzing disease comorbidity networks, using knowledge graphs that integrate knowledge from different medical domains, and applying network analysis and graph ML. Multimorbidity disease networks have been used to improve disease diagnosis, treatment recommendations, and patient prognosis. Knowledge graphs that combine different medical entities connected by multiple relationship types integrate data from different sources, allowing for complex interactions and creating a continuous flow of information. Network analysis and graph ML can then extract the topology and structure of networks and reveal hidden properties, including disease phenotypes, network hubs, and pathways; predict drugs for repurposing; and determine safe and more holistic treatments. In this article, we describe the basic concepts of creating bipartite and unipartite disease and patient networks and review the use of knowledge graphs, graph algorithms, graph embedding methods, and graph ML within the context of multimorbidity. Specifically, we provide an overview of the application of graph theory for studying multimorbidity, the methods employed to extract knowledge from graphs, and examples of the application of disease networks for determining the structure and pathways of multimorbidity, identifying disease phenotypes, predicting health outcomes, and selecting safe and effective treatments. In today's modern data-hungry, ML-focused world, such network-based techniques are likely to be at the forefront of developing robust clinical decision support tools for safer and more holistic approaches to treating older patients with multimorbidity.

Arc regulates a second-guessing cognitive bias during naturalistic foraging through effects on discrete behavior modules

Foraging in animals relies on innate decision-making heuristics that can result in suboptimal cognitive biases in some contexts. The mechanisms underlying these biases are not well understood, but likely involve strong genetic effects. To explore this, we studied fasted mice using a naturalistic foraging paradigm and discovered an innate cognitive bias called "second-guessing." This involves repeatedly investigating an empty former food patch instead of consuming available food, which hinders the mice from maximizing feeding benefits. The synaptic plasticity gene Arc is revealed to play a role in this bias, as Arc-deficient mice did not exhibit second-guessing and consumed more food. In addition, unsupervised machine learning decompositions of foraging identified specific behavior sequences, or "modules", that are affected by Arc. These findings highlight the genetic basis of cognitive biases in decision making, show links between behavior modules and cognitive bias, and provide insight into the ethological roles of Arc in naturalistic foraging.

Problem of Domain/Building Block Preservation in the Evolution of Biological Macromolecules and Evolutionary Computation

Structurally and functionally isolated domains in biological macromolecular evolution, both natural and artificial, are largely similar to "schemata", building blocks (BBs), in evolutionary computation (EC). The problem of preserving in subsequent evolutionary searches the already found domains / BBs is well known and quite relevant in biology as well as in EC. Both biology and EC are seeing parallel and independent development of several approaches to identifying and preserving previously identified domains / BBs. First, we notice the similarity of DNA shuffling methods in synthetic biology and multi-parent recombination algorithms in EC. Furthermore, approaches to computer identification of domains in proteins that are being developed in biology can be aligned with BB identification methods in EC. Finally, approaches to chimeric protein libraries optimization in biology can be compared to evolutionary search methods based on probabilistic models in EC. We propose to validate the prospects of mutual exchange of ideas and transfer of algorithms and approaches between evolutionary systems biology and EC in these three principal directions. A crucial aim of this transfer is the design of new advanced experimental techniques capable of solving more complex problems of in vitro evolution.

Cellular and Natural Viral Engineering in Cognition-Based Evolution

Neo-Darwinism conceptualizes evolution as the continuous succession of predominately random genetic variations disciplined by natural selection. In that frame, the primary interaction between cells and the virome is relegated to host-parasite dynamics governed by selective influences. Cognition-Based Evolution regards biological and evolutionary development as a reciprocating cognition-based informational interactome for the protection of self-referential cells. To sustain cellular homeorhesis, cognitive cells collaborate to assess the validity of ambiguous biological information. That collective interaction involves coordinate measurement, communication, and active deployment of resources as Natural Cellular Engineering. These coordinated activities drive multicellularity, biological development, and evolutionary change. The virome participates as the vital intercessory among the cellular domains to ensure their shared permanent perpetuation. The interactions between the virome and the cellular domains represent active virocellular cross-communications for the continual exchange of resources. Modular genetic transfers between viruses and cells carry bioactive potentials. Those exchanges are deployed as nonrandom flexible tools among the domains in their continuous confrontation with environmental stresses. This alternative framework fundamentally shifts our perspective on viral-cellular interactions, strengthening established principles of viral symbiogenesis. Pathogenesis can now be properly appraised as one expression of a range of outcomes between cells and viruses within a larger conceptual framework of Natural Viral Engineering as a co-engineering participant with cells. It is proposed that Natural Viral Engineering should be viewed as a co-existent facet of Natural Cellular Engineering within Cognition-Based Evolution.

Modular Screening Reveals Driver Induced Additive Mechanisms of Baicalin and Jasminoidin on Cerebral Ischemia Therapy

Combination therapy with increased efficacy and reduced toxicity plays a crucial role in treating complex diseases, such as stroke, but it remains an insurmountable barrier to elucidate the mechanisms of synergistic effects. Here, we present a Driver-induced Modular Screening (DiMS) strategy integrated synergistic module and driver gene identification to elucidate the additive mechanisms of Baicalin (BA) and Jasminoidin (JA) on cerebral ischemia (CI) therapy. Based on anti-ischemia genomic networks BA, JA, and their combination (BJ), we obtained 4, 3, and 9 On-modules of BA, JA, and BJ by modular similarity analysis. Compared with the monotherapy groups, four additive modules (Add-module, BJ_Mod-4, 7, 9, and 13), 15 driver genes of BJ were identified by modular similarity and network control methods, and seven driver proteins (PAQR8, RhoA, EMC10, GGA2, VIPR1, FAM120A, and SEMA3F) were validated by animal experiments. The functional analysis found neuroprotective roles of the Add-modules and driver genes, such as the Neurotrophin signaling pathway and FoxO signaling pathway, which may reflect the additive mechanisms of BJ. Moreover, such a DiMS paradigm provides a new angle to explore the synergistic mechanisms of combination therapy and screen multi-targeted drugs for complex diseases.

A weighted network analysis framework for the hourglass effect-And its application in the C. elegans connectome

Understanding hierarchy and modularity in natural as well as technological networks is of utmost importance. A major aspect of such analysis involves identifying the nodes that are crucial to the overall processing structure of the network. More recently, the approach of hourglass analysis has been developed for the purpose of quantitatively analyzing whether only a few intermediate nodes mediate the information processing between a large number of inputs and outputs of a network. We develop a new framework for hourglass analysis that takes network weights into account while identifying the core nodes and the extent of hourglass effect in a given weighted network. We use this framework to study the structural connectome of the C. elegans and identify intermediate neurons that form the core of sensori-motor pathways in the organism. Our results show that the neurons forming the core of the connectome show significant differences across the male and hermaphrodite sexes, with most core nodes in the male concentrated in sex-organs while they are located in the head for the hermaphrodite. Our work demonstrates that taking weights into account for network analysis framework leads to emergence of different network patterns in terms of identification of core nodes and hourglass structure in the network, which otherwise would be missed by unweighted approaches.

Modularity in Biological Networks

Network modeling, from the ecological to the molecular scale has become an essential tool for studying the structure, dynamics and complex behavior of living systems. Graph representations of the relationships between biological components open up a wide variety of methods for discovering the mechanistic and functional properties of biological systems. Many biological networks are organized into a modular structure, so methods to discover such modules are essential if we are to understand the biological system as a whole. However, most of the methods used in biology to this end, have a limited applicability, as they are very specific to the system they were developed for. Conversely, from the statistical physics and network science perspective, graph modularity has been theoretically studied and several methods of a very general nature have been developed. It is our perspective that in particular for the modularity detection problem, biology and theoretical physics/network science are less connected than they should. The central goal of this review is to provide the necessary background and present the most applicable and pertinent methods for community detection in a way that motivates their further usage in biological research.

Machine Learning for Causal Inference in Biological Networks: Perspectives of This Challenge

Most machine learning-based methods predict outcomes rather than understanding causality. Machine learning methods have been proved to be efficient in finding correlations in data, but unskilful to determine causation. This issue severely limits the applicability of machine learning methods to infer the causal relationships between the entities of a biological network, and more in general of any dynamical system, such as medical intervention strategies and clinical outcomes system, that is representable as a network. From the perspective of those who want to use the results of network inference not only to understand the mechanisms underlying the dynamics, but also to understand how the network reacts to external stimuli (e. g. environmental factors, therapeutic treatments), tools that can understand the causal relationships between data are highly demanded. Given the increasing popularity of machine learning techniques in computational biology and the recent literature proposing the use of machine learning techniques for the inference of biological networks, we would like to present the challenges that mathematics and computer science research faces in generalising machine learning to an approach capable of understanding causal relationships, and the prospects that achieving this will open up for the medical application domains of systems biology, the main paradigm of which is precisely network biology at any physical scale.

Revealing the Modular Similarities and Differences Among Alzheimer's Disease, Vascular Dementia, and Parkinson's Disease in Genomic Networks

Alzheimer's disease (AD), vascular dementia (VD), and Parkinson's disease (PD) exert increasingly lethal or disabling effects on humans, but the associations among these diseases at the molecular level remain unclear. In our research, lists of genes related to these three diseases were acquired from public databases. We constructed gene-gene networks of the lists of disease-related genes using the STRING database and selected the plug-in MCODE as the most suitable method to divide the three disease-associated networks into modules through an entropy calculation. Notably, 1173 AD-related, 203 VD-related, and 722 PD-related genes as well as 72 overlapping genes were observed among the three diseases. By dividing the modules from the gene network, we divided the AD-related gene network into 27 modules, the VD-related gene network into 8 modules, and the PD-related gene network into 17 modules. After the enrichment analysis of each disease-related gene, 146 overlapping biological processes and 32 overlapping pathways were identified. Ultimately, through similarity analysis of the genes, biological processes, and pathways, we found that AD and VD were the most closely related at the biological process and pathway levels, with similarity coefficients of 0.2784 and 0.3626, respectively. After analyzing the overlapping gene network, we found that INS might play an important role in the network and that insulin and its signaling pathways may play a key role in these neurodegenerative diseases. Our research illustrates a new method for in-depth research on the three diseases, which may accelerate the progress of developing new therapeutics and may be applied to prevent neurodegenerative diseases.

Gene coexpression networks reveal molecular interactions underlying cichlid jaw modularity

Background

The oral and pharyngeal jaw of cichlid fishes are a classic example of evolutionary modularity as their functional decoupling boosted trophic diversification and contributed to the success of cichlid adaptive radiations. Most studies until now have focused on the functional, morphological, or genetic aspects of cichlid jaw modularity. Here we extend this concept to include transcriptional modularity by sequencing whole transcriptomes of the two jaws and comparing their gene coexpression networks.

Results

We show that transcriptional decoupling of gene expression underlies the functional decoupling of cichlid oral and pharyngeal jaw apparatus and the two units are evolving independently in recently diverged cichlid species from Lake Tanganyika. Oral and pharyngeal jaw coexpression networks reflect the common origin of the jaw regulatory program as there is high preservation of gene coexpression modules between the two sets of jaws. However, there is substantial rewiring of genetic architecture within those modules. We define a global jaw coexpression network and highlight jaw-specific and species-specific modules within it. Furthermore, we annotate a comprehensive in silico gene regulatory network linking the Wnt and AHR signalling pathways to jaw morphogenesis and response to environmental cues, respectively. Components of these pathways are significantly differentially expressed between the oral and pharyngeal jaw apparatus.

Conclusion

This study describes the concerted expression of many genes in cichlid oral and pharyngeal jaw apparatus at the onset of the independent life of cichlid fishes. Our findings suggest that – on the basis of an ancestral gill arch network—transcriptional rewiring may have driven the modular evolution of the oral and pharyngeal jaws, highlighting the evolutionary significance of gene network reuse. The gene coexpression and in silico regulatory networks presented here are intended as resource for future studies on the genetics of vertebrate jaw morphogenesis and trophic adaptation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01787-9.

Decoupling transcription factor expression and activity enables dimmer switch gene regulation

Gene-regulatory networks achieve complex mappings of inputs to outputs through mechanisms that are poorly understood. We found that in the galactose-responsive pathway in Saccharomyces cerevisiae, the decision to activate the transcription of genes encoding pathway components is controlled independently from the expression level, resulting in behavior resembling that of a mechanical dimmer switch. This was not a direct result of chromatin regulation or combinatorial control at galactose-responsive promoters; rather, this behavior was achieved by hierarchical regulation of the expression and activity of a single transcription factor. Hierarchical regulation is ubiquitous, and thus dimmer switch regulation is likely a key feature of many biological systems. Dimmer switch gene regulation may allow cells to fine-tune their responses to multi-input environments on both physiological and evolutionary time scales.

Computational animal welfare: towards cognitive architecture models of animal sentience, emotion and wellbeing

To understand animal wellbeing, we need to consider subjective phenomena and sentience. This is challenging, since these properties are private and cannot be observed directly. Certain motivations, emotions and related internal states can be inferred in animals through experiments that involve choice, learning, generalization and decision-making. Yet, even though there is significant progress in elucidating the neurobiology of human consciousness, animal consciousness is still a mystery. We propose that computational animal welfare science emerges at the intersection of animal behaviour, welfare and computational cognition. By using ideas from cognitive science, we develop a functional and generic definition of subjective phenomena as any process or state of the organism that exists from the first-person perspective and cannot be isolated from the animal subject. We then outline a general cognitive architecture to model simple forms of subjective processes and sentience. This includes evolutionary adaptation which contains top-down attention modulation, predictive processing and subjective simulation by re-entrant (recursive) computations. Thereafter, we show how this approach uses major characteristics of the subjective experience: elementary self-awareness, global workspace and qualia with unity and continuity. This provides a formal framework for process-based modelling of animal needs, subjective states, sentience and wellbeing.

A Middle-Out Modeling Strategy to Extend a Colon Cancer Logical Model Improves Drug Synergy Predictions in Epithelial-Derived Cancer Cell Lines

Cancer is a heterogeneous and complex disease and one of the leading causes of death worldwide. The high tumor heterogeneity between individuals affected by the same cancer type is accompanied by distinct molecular and phenotypic tumor profiles and variation in drug treatment response. In silico modeling of cancer as an aberrantly regulated system of interacting signaling molecules provides a basis to enhance our biological understanding of disease progression, and it offers the means to use computer simulations to test and optimize drug therapy designs on particular cancer types and subtypes. This sets the stage for precision medicine: the design of treatments tailored to individuals or groups of patients based on their tumor-specific molecular cancer profiles. Here, we show how a relatively large manually curated logical model can be efficiently enhanced further by including components highlighted by a multi-omics data analysis of data from Consensus Molecular Subtypes covering colorectal cancer. The model expansion was performed in a pathway-centric manner, following a partitioning of the model into functional subsystems, named modules. The resulting approach constitutes a middle-out modeling strategy enabling a data-driven expansion of a model from a generic and intermediate level of molecular detail to a model better covering relevant processes that are affected in specific cancer subtypes, comprising 183 biological entities and 603 interactions between them, partitioned in 25 functional modules of varying size and structure. We tested this model for its ability to correctly predict drug combination synergies, against a dataset of experimentally determined cell growth responses with 18 drugs in all combinations, on eight cancer cell lines. The results indicate that the extended model had an improved accuracy for drug synergy prediction for the majority of the experimentally tested cancer cell lines, although significant improvements of the model's predictive performance are still needed. Our study demonstrates how a tumor-data driven middle-out approach toward refining a logical model of a biological system can further customize a computer model to represent specific cancer cell lines and provide a basis for identifying synergistic effects of drugs targeting specific regulatory proteins. This approach bridges between preclinical cancer model data and clinical patient data and may thereby ultimately be of help to develop patient-specific in silico models that can steer treatment decisions in the clinic.

Expression of long noncoding RNAs in cancer-associated fibroblasts linked to patient survival in ovarian cancer

Cancer-associated fibroblasts (CAFs) are the most abundant cell type in the tumor microenvironment and are responsible for producing the desmoplastic reaction that is a poor prognostic factor in ovarian cancer. Long non-coding RNAs (lncRNAs) have been shown to play important roles in cancer. However, very little is known about the role of lncRNAs in the tumor microenvironment. We aimed to identify lncRNAs expressed in ovarian CAFs that were associated with patient survival and used computational approaches to predict their function. Increased expression of 9 lncRNAs and decreased expression of 1 lncRNA in ovarian CAFs were found to be associated with poorer overall survival. A "guilt-by-association" approach was used to predict the function of these lncRNAs. In particular, MIR155HG was predicted to play a role in immune response. Further investigation revealed high MIR155HG expression to be associated with higher infiltrates of immune cell subsets. In conclusion, these data indicate expression on several lncRNAs in CAFs are associated with patient survival and are likely to play an important role in regulating CAF function.

A model for the interplay between plastic tradeoffs and evolution in changing environments

Performance tradeoffs are ubiquitous in both ecological and evolutionary modeling, yet they are usually postulated and built into fitness and ecological landscapes. However, tradeoffs depend on genetic background and evolutionary history and can themselves evolve. We present a simple model capable of capturing the key feedback loop: evolutionary history shapes tradeoff strength, which, in turn, shapes evolutionary future. One consequence of this feedback is that genomes with identical fitness can have different evolutionary properties shaped by prior environmental exposure. Another is that, generically, the best adaptations to one environment may evolve in another. Our simple framework bridges the gap between the phenotypic Fisher's Geometric Model and the genotypic properties, such as modularity and evolvability, and can serve as a rich playground for investigating evolution in multiple or changing environments.

On the emergence of cognition: from catalytic closure to neuroglial closure

In an analogous manner as occurred during the development of a connected metabolism that at some point reached characteristics associated with what is called "life"-due mainly to a catalytic closure phenomenon when chemicals started to autocatalyze themselves forming a closed web of chemical reactions-it is here proposed that cognition and consciousness (or features associated with them) arose as a consequence of another type of closure within the nervous system, the brain especially. Proper brain function requires an efficient web of connections and once certain complexity is attained due to the number and coordinated activities of the brain cell networks, the emergent properties of cognition and consciousness take place. Seeking to identify main features of the nervous system organization for optimal function, it is here proposed that while catalytic closure yielded life, neuroglial closure produced cognition/consciousness.

Exploring modularity in biological networks

Network theoretical approaches have shaped our understanding of many different kinds of biological modularity. This essay makes the case that to capture these contributions, it is useful to think about the role of network models in exploratory research. The overall point is that it is possible to provide a systematic analysis of the exploratory functions of network models in bioscientific research. Using two examples from molecular and developmental biology, I argue that often the same modelling approach can perform one or more exploratory functions, such as introducing new directions of research, offering a complementary set of concepts, methods and algorithms for individuating important features of natural phenomena, generating proofs of principle demonstrations and potential explanations for phenomena of interest and enlarging the scope of certain research agendas. This article is part of the theme issue 'Unifying the essential concepts of biological networks: biological insights and philosophical foundations'.

The hourglass organization of the Caenorhabditis elegans connectome

We approach the C. elegans connectome as an information processing network that receives input from about 90 sensory neurons, processes that information through a highly recurrent network of about 80 interneurons, and it produces a coordinated output from about 120 motor neurons that control the nematode's muscles. We focus on the feedforward flow of information from sensory neurons to motor neurons, and apply a recently developed network analysis framework referred to as the "hourglass effect". The analysis reveals that this feedforward flow traverses a small core ("hourglass waist") that consists of 10-15 interneurons. These are mostly the same interneurons that were previously shown (using a different analytical approach) to constitute the "rich-club" of the C. elegans connectome. This result is robust to the methodology that separates the feedforward from the feedback flow of information. The set of core interneurons remains mostly the same when we consider only chemical synapses or the combination of chemical synapses and gap junctions. The hourglass organization of the connectome suggests that C. elegans has some similarities with encoder-decoder artificial neural networks in which the input is first compressed and integrated in a low-dimensional latent space that encodes the given data in a more efficient manner, followed by a decoding network through which intermediate-level sub-functions are combined in different ways to compute the correlated outputs of the network. The core neurons at the hourglass waist represent the information bottleneck of the system, balancing the representation accuracy and compactness (complexity) of the given sensory information.

Stability of spontaneous, correlated activity in mouse auditory cortex

Neural systems can be modeled as complex networks in which neural elements are represented as nodes linked to one another through structural or functional connections. The resulting network can be analyzed using mathematical tools from network science and graph theory to quantify the system’s topological organization and to better understand its function. Here, we used two-photon calcium imaging to record spontaneous activity from the same set of cells in mouse auditory cortex over the course of several weeks. We reconstruct functional networks in which cells are linked to one another by edges weighted according to the correlation of their fluorescence traces. We show that the networks exhibit modular structure across multiple topological scales and that these multi-scale modules unfold as part of a hierarchy. We also show that, on average, network architecture becomes increasingly dissimilar over time, with similarity decaying monotonically with the distance (in time) between sessions. Finally, we show that a small fraction of cells maintain strongly-correlated activity over multiple days, forming a stable temporal core surrounded by a fluctuating and variable periphery. Our work indicates a framework for studying spontaneous activity measured by two-photon calcium imaging using computational methods and graphical models from network science. The methods are flexible and easily extended to additional datasets, opening the possibility of studying cellular level network organization of neural systems and how that organization is modulated by stimuli or altered in models of disease.

Neurons coordinate their activity with one another, forming networks that help support adaptive, flexible behavior. Still, little is known about the organization of these networks at the cellular scale and their stability over time. Here, we reconstruct networks from calcium imaging data recorded in mouse primary auditory cortex. We show that these networks exhibit spatially constrained, hierarchical modular structure, which may facilitate specialized information processing. However, we show that connection weights and modular structure are also variable over time, changing on a timescale of days and adopting novel network configurations. Despite this, a small subset of neurons maintain their connections to one another and preserve their modular organization across time, forming a stable temporal core surrounded by a flexible periphery. These findings represent a conceptual bridge linking network analyses of macroscale and cellular-level neuroimaging data. They also represent a complementary approach to existing circuits- and systems-based interrogation of nervous system function, opening the door for deeper and more targeted analysis in the future.

The Eminence of Co-Expressed Ties in Schizophrenia Network Communities

Exploring gene networks is crucial for identifying significant biological interactions occurring in a disease condition. These interactions can be acknowledged by modeling the tie structure of networks. Such tie orientations are often detected within embedded community structures. However, most of the prevailing community detection modules are intended to capture information from nodes and its attributes, usually ignoring the ties. In this study, a modularity maximization algorithm is proposed based on nonlinear representation of local tangent space alignment (LTSA). Initially, the tangent coordinates are computed locally to identify k-nearest neighbors across the genes. These local neighbors are further optimized by generating a nonlinear network embedding function for detecting gene communities based on eigenvector decomposition. Experimental results suggest that this algorithm detects gene modules with a better modularity index of 0.9256, compared to other traditional community detection algorithms. Furthermore, co-expressed genes across these communities are identified by discovering the characteristic tie structures. These detected ties are known to have substantial biological influence in the progression of schizophrenia, thereby signifying the influence of tie patterns in biological networks. This technique can be extended logically on other diseases networks for detecting substantial gene “hotspots”.

Reusable building blocks in biological systems

One of the most widely recognized features of biological systems is their modularity. The modules that constitute biological systems are said to be redeployed and combined across several conditions, thus acting as building blocks. In this work, we analyse to what extent are these building blocks reusable as compared with those found in randomized versions of a system. We develop a notion of decompositions of systems into phenotypic building blocks, which allows them to overlap while maximizing the number of times a building block is reused across several conditions. Different biological systems present building blocks whose reusability ranges from single use (e.g. condition specific) to constitutive, although their average reusability is not always higher than random equivalents of the system. These decompositions reveal a distinct distribution of building block sizes in real biological systems. This distribution stems, in part, from the peculiar usage pattern of the elements of biological systems, and constitutes a new angle to study the evolution of modularity.

Divergent Allometric Trajectories in Gene Expression and Coexpression Produce Species Differences in Sympatrically Speciating Midas Cichlid Fish

The mechanisms of speciation without geographic isolation (i.e., sympatric speciation) remain debated. This is due in part to the fact that the genomic landscape that could promote or hinder species divergence in the presence of gene flow is still largely unknown. However, intensive research is now centered on understanding the genetic architecture of adaptive traits associated with this process as well as how gene expression might affect these traits. Here, using RNA-Seq data, we investigated gene expression of sympatrically speciating benthic and limnetic Neotropical cichlid fishes at two developmental stages. First, we identified groups of coexpressed genes (modules) at each stage. Although there are a few large and well-preserved modules, most of the other modules are not preserved across life stages. Second, we show that later in development more and larger coexpression modules are associated with divergence between benthic and limnetic fish compared with the earlier life stage. This divergence between benthic and limnetic fish in coexpression mirrors divergence in overall expression between benthic and limnetic fish, which is more pronounced later in life. Our results reveal that already at 1-day posthatch benthic and limnetic fish diverge in (co)expression, and that this divergence becomes more substantial when fish are free-swimming but still unlikely to have divergent swimming and feeding habits. More importantly, our study describes how the coexpression of several genes through development, as opposed to individual genes, is associated with benthic–limnetic species differences, and how two morphogenetic trajectories diverge as fish grow older.

A Network Pharmacology Approach to Explore the Mechanisms of Qishen Granules in Heart Failure

This study aimed to investigate the intrinsic mechanisms of Qishen granules (QSG) in the treatment of HF, and to provide new evidence and insights for its clinical application. Information on QSG ingredients was collected from Traditional Chinese medicine systems pharmacology (TCMSP), TCM@Taiwan, TCMID, and Batman, and input into SwissTargetPrediction to identify the compound targets. HF-related targets were detected from Therapeutic Target Database (TTD), Disgenet-Gene, Drugbank database, and Online Mendelian Inheritance in Man (OMIM) database. The overlap targets of QSG and HF were identified for pathway enrichment analysis by utilizing the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. The protein-protein interaction (PPI) network of QSG-HF was constructed, following by the generation of core targets, construction of core modules, and KEGG analysis of the core functional modules. There were 1909 potential targets predicted from the 243 bioactive compounds in QSG which shared 129 common targets with HF-related targets. KEGG pathway analysis of common targets indicated that QSG could regulated 23 representative pathways. In the QSG-HF PPI network analysis, 10 key targets were identified, including EDN1, AGT, CREB1, ACE, CXCR4, ADRBK1, AGTR1, BDKRB1, ADRB2, and F2. Further cluster and enrichment analysis suggested that neuroactive ligand-receptor interaction, cGMP-PKG signaling pathway, renin secretion, vascular smooth muscle contraction, and the renin-angiotensin system might be core pathways of QSG for HF. Our study elucidated the possible mechanisms of QSG from a systemic and holistic perspective. The key targets and pathways will provide new insights for further research on the pharmacological mechanism of QSG.

On a Simple General Principle of Brain Organization

A possible framework to characterize nervous system dynamics and its organization in conscious and unconscious states is introduced, derived from a high level perspective on the coordinated activity of brain cell ensembles. Some questions are best addressable in a global framework and here we build on past observations about the structure of configurations of brain networks in conscious and unconscious states and about neurophysiological results. Aiming to bind some results together into some sort of coherence with a central theme, the scenario that emerges underscores the crucial importance of the creation and dissipation of energy gradients in brain cellular ensembles resulting in maximization of the configurations in the functional connectivity among those networks that favor conscious awareness and healthy conditions. These considerations are then applied to indicate approaches that can be used to improve neuropathological syndromes.

Emergence of Hierarchical Modularity in Evolving Networks Uncovered by Phylogenomic Analysis

Networks describe how parts associate with each other to form integrated systems which often have modular and hierarchical structure. In biology, network growth involves two processes, one that unifies and the other that diversifies. Here, we propose a biphasic (bow-tie) theory of module emergence. In the first phase, parts are at first weakly linked and associate variously. As they diversify, they compete with each other and are often selected for performance. The emerging interactions constrain their structure and associations. This causes parts to self-organize into modules with tight linkage. In the second phase, variants of the modules diversify and become new parts for a new generative cycle of higher level organization. The paradigm predicts the rise of hierarchical modularity in evolving networks at different timescales and complexity levels. Remarkably, phylogenomic analyses uncover this emergence in the rewiring of metabolomic and transcriptome-informed metabolic networks, the nanosecond dynamics of proteins, and evolving networks of metabolism, elementary functionomes, and protein domain organization.

Penicillin Trunk Injection Affects Bacterial Community Structure in Citrus Trees

Huanglongbing (HLB), caused by Candidatus Liberibacter asiaticus (CLas), an uncultured α-proteobacterium, is the most destructive disease of citrus trees worldwide. In previous studies, trunk injections of penicillin reduced CLas titers and HLB symptoms in citrus. However, antibiotic effects on the whole plant microbial community, which include effects on taxa that interact with CLas, have not yet been addressed. In this study, we investigated the effects of penicillin injection (0, 1000, and 6000 mg L^-1) on rhizospheric and endophytic bacterial communities of grapefruit trees in field and greenhouse experiments through culture-independent high-throughput sequencing. DNA extractions from petioles and roots were subjected to 16S rRNA high-throughput sequencing, and reads were clustered by sequence similarity into operational taxonomic units (OTUs). Principal coordinates analysis based on weighted-UniFrac distances did not reveal differences in bacterial communities among treatments in any of the sample sources. However, pairwise linear discriminant analysis indicated significant differences in relative abundance of some taxa (including CLas) among treatments. Network analysis showed that penicillin produced major changes in root bacterial community structure by affecting interspecific microbial associations. This study provides new knowledge of the effect of antimicrobial treatments on interspecific relationships in citrus microbial communities.

The Energy Homeostasis Principle: Neuronal Energy Regulation Drives Local Network Dynamics Generating Behavior

A major goal of neuroscience is understanding how neurons arrange themselves into neural networks that result in behavior. Most theoretical and experimental efforts have focused on a top-down approach which seeks to identify neuronal correlates of behaviors. This has been accomplished by effectively mapping specific behaviors to distinct neural patterns, or by creating computational models that produce a desired behavioral outcome. Nonetheless, these approaches have only implicitly considered the fact that neural tissue, like any other physical system, is subjected to several restrictions and boundaries of operations. Here, we proposed a new, bottom-up conceptual paradigm: The Energy Homeostasis Principle, where the balance between energy income, expenditure, and availability are the key parameters in determining the dynamics of neuronal phenomena found from molecular to behavioral levels. Neurons display high energy consumption relative to other cells, with metabolic consumption of the brain representing 20% of the whole-body oxygen uptake, contrasting with this organ representing only 2% of the body weight. Also, neurons have specialized surrounding tissue providing the necessary energy which, in the case of the brain, is provided by astrocytes. Moreover, and unlike other cell types with high energy demands such as muscle cells, neurons have strict aerobic metabolism. These facts indicate that neurons are highly sensitive to energy limitations, with Gibb's free energy dictating the direction of all cellular metabolic processes. From this activity, the largest energy, by far, is expended by action potentials and post-synaptic potentials; therefore, plasticity can be reinterpreted in terms of their energy context. Consequently, neurons, through their synapses, impose energy demands over post-synaptic neurons in a close loop-manner, modulating the dynamics of local circuits. Subsequently, the energy dynamics end up impacting the homeostatic mechanisms of neuronal networks. Furthermore, local energy management also emerges as a neural population property, where most of the energy expenses are triggered by sensory or other modulatory inputs. Local energy management in neurons may be sufficient to explain the emergence of behavior, enabling the assessment of which properties arise in neural circuits and how. Essentially, the proposal of the Energy Homeostasis Principle is also readily testable for simple neuronal networks.

Analysis of Hierarchical Organization in Gene Expression Networks Reveals Underlying Principles of Collective Tumor Cell Dissemination and Metastatic Aggressiveness of Inflammatory Breast Cancer

Clusters of circulating tumor cells (CTCs), despite being rare, may account for more than 90% of metastases. Cells in these clusters do not undergo a complete epithelial-to-mesenchymal transition (EMT), but retain some epithelial traits as compared to individually disseminating tumor cells. Determinants of single cell dissemination versus collective dissemination remain elusive. Inflammatory breast cancer (IBC), a highly aggressive breast cancer subtype that chiefly metastasizes via CTC clusters, is a promising model for studying mechanisms of collective tumor cell dissemination. Previous studies, motivated by a theory that suggests physical systems with hierarchical organization tend to be more adaptable, have found that the expression of metastasis-associated genes is more hierarchically organized in cases of successful metastases. Here, we used the cophenetic correlation coefficient (CCC) to quantify the hierarchical organization in the expression of two distinct gene sets, collective dissemination-associated genes and IBC-associated genes, in cancer cell lines and in tumor samples from breast cancer patients. Hypothesizing that a higher CCC for collective dissemination-associated genes and for IBC-associated genes would be associated with retention of epithelial traits enabling collective dissemination and with worse disease progression in breast cancer patients, we evaluated the correlation of CCC with different phenotypic groups. The CCC of both the abovementioned gene sets, the collective dissemination-associated genes and the IBC-associated genes, was higher in (a) epithelial cell lines as compared to mesenchymal cell lines and (b) tumor samples from IBC patients as compared to samples from non-IBC breast cancer patients. A higher CCC of both gene sets was also correlated with a higher rate of metastatic relapse in breast cancer patients. In contrast, neither the levels of CDH1 gene expression nor gene set enrichment analysis (GSEA) of the abovementioned gene sets could provide similar insights. These results suggest that retention of some epithelial traits in disseminating tumor cells as IBC progresses promotes successful breast cancer metastasis. The CCC provides additional information regarding the organizational complexity of gene expression in comparison to GSEA. We have shown that the CCC may be a useful metric for investigating the collective dissemination phenotype and a prognostic factor for IBC.

How Organisms Gained Causal Independence and How It Might Be Quantified

Two broad features are jointly necessary for autonomous agency: organisational closure and the embodiment of an objective-function providing a ‘goal’: so far only organisms demonstrate both. Organisational closure has been studied (mostly in abstract), especially as cell autopoiesis and the cybernetic principles of autonomy, but the role of an internalised ‘goal’ and how it is instantiated by cell signalling and the functioning of nervous systems has received less attention. Here I add some biological ‘flesh’ to the cybernetic theory and trace the evolutionary development of step-changes in autonomy: (1) homeostasis of organisationally closed systems; (2) perception-action systems; (3) action selection systems; (4) cognitive systems; (5) memory supporting a self-model able to anticipate and evaluate actions and consequences. Each stage is characterised by the number of nested goal-directed control-loops embodied by the organism, summarised as will-nestedness N. Organism tegument, receptor/transducer system, mechanisms of cellular and whole-organism re-programming and organisational integration, all contribute to causal independence.

CONCLUSION

organisms are cybernetic phenomena whose identity is created by the information structure of the highest level of causal closure (maximum N), which has increased through evolution, leading to increased causal independence, which might be quantifiable by ‘Integrated Information Theory’ measures.

On the role of sparseness in the evolution of modularity in gene regulatory networks

Modularity is a widespread property in biological systems. It implies that interactions occur mainly within groups of system elements. A modular arrangement facilitates adjustment of one module without perturbing the rest of the system. Therefore, modularity of developmental mechanisms is a major factor for evolvability, the potential to produce beneficial variation from random genetic change. Understanding how modularity evolves in gene regulatory networks, that create the distinct gene activity patterns that characterize different parts of an organism, is key to developmental and evolutionary biology. One hypothesis for the evolution of modules suggests that interactions between some sets of genes become maladaptive when selection favours additional gene activity patterns. The removal of such interactions by selection would result in the formation of modules. A second hypothesis suggests that modularity evolves in response to sparseness, the scarcity of interactions within a system. Here I simulate the evolution of gene regulatory networks and analyse diverse experimentally sustained networks to study the relationship between sparseness and modularity. My results suggest that sparseness alone is neither sufficient nor necessary to explain modularity in gene regulatory networks. However, sparseness amplifies the effects of forms of selection that, like selection for additional gene activity patterns, already produce an increase in modularity. That evolution of new gene activity patterns is frequent across evolution also supports that it is a major factor in the evolution of modularity. That sparseness is widespread across gene regulatory networks indicates that it may have facilitated the evolution of modules in a wide variety of cases.

Modular systems have performance and design advantages over non-modular systems. Thus, modularity is very important for the development of a wide range of new technological or clinical applications. Moreover, modularity is paramount to evolutionary biology since it allows adjusting one organismal function without disturbing other previously evolved functions. But how does modularity itself evolve? Here I analyse the structure of regulatory networks and follow simulations of network evolution to study two hypotheses for the origin of modules in gene regulatory networks. The first hypothesis considers that sparseness, a low number of interactions among the network genes, could be responsible for the evolution of modular networks. The second, that modules evolve when selection favours the production of additional gene activity patterns. I found that sparseness alone is neither sufficient nor necessary to explain modularity in gene regulatory networks. However, it enhances the effects of selection for multiple gene activity patterns. While selection for multiple patterns may be decisive in the evolution of modularity, that sparseness is widespread across gene regulatory networks suggests that its contributions should not be neglected.

Ecological and evolutionary dynamics of interconnectedness and modularity

In this contribution, we develop a theoretical framework for linking microprocesses (i.e., population dynamics and evolution through natural selection) with macrophenomena (such as interconnectedness and modularity within an ecological system). This is achieved by developing a measure of interconnectedness for population distributions defined on a trait space (generalizing the notion of modularity on graphs), in combination with an evolution equation for the population distribution. With this contribution, we provide a platform for understanding under what environmental, ecological, and evolutionary conditions ecosystems evolve toward being more or less modular. A major contribution of this work is that we are able to decompose the overall driver of changes at the macro level (such as interconnectedness) into three components: (i) ecologically driven change, (ii) evolutionarily driven change, and (iii) environmentally driven change.

Network-Wide Screen Identifies Variation of Novel Precise On-Module Targets Using Conformational Modudaoism

Modular targeting is promising in drug research at the network level, but it is challenging to quantificationally identify the precise on-modules. Based on a proposed Modudaoism (MD), we defined conserved MD (MDc) and varied MD (MDv) to quantitatively evaluate the conformational and energy variations of modules, and thereby identify the conserved and discrepant allosteric modules (AMs). Compared to the Z_summary , MDc/MDv got an optimized result of module preserved ratio and modular structure. In the mice anti-ischemic networks, 3, 5, and 1 conserved AMs as well as 4, 1, and 3 on-modules of baicalin (BA), jasminoidin (JA), and ursodeoxycholic acid (UA) were identified by MDc and MDv, 5 unique AMs and their characteristic actions were revealed. Besides, co-immunoprecipitation (Co-IP) experiments validated the representative modular structure. MDc/MDv method can quantitatively define the conformational variations of modules and screen the precise on-modules network-wide, which may provide a promising strategy for drug discovery.

A Statistical Description of Plant Shoot Architecture

Plant architectures can be characterized statistically by their spatial density function, which specifies the probability of finding a branch at each location in the territory occupied by a plant. Using high-precision 3D scanning, we analyzed 557 plant shoot architectures, representing three species, grown across three to five environmental conditions, and through 20-30 developmental time points. We found two elegant properties in the spatial density functions of these architectures: all functions could be nearly modified in one direction without affecting the density in orthogonal directions (called "separability"), and all functions shared the same underlying shape, aside from stretching and compression (called "self-similarity"). Surprisingly, despite their striking visual diversity, we discovered that all architectures could be described as variations on a single underlying function: a Gaussian density function truncated at roughly two SDs. We also observed systematic variation in the spatial density functions across species, growth conditions, and time, which suggests functional specialization despite following the same general design form.

Plasmid Replicons from Pseudomonas Are Natural Chimeras of Functional, Exchangeable Modules

Plasmids are a main factor for the evolution of bacteria through horizontal gene exchange, including the dissemination of pathogenicity genes, resistance to antibiotics and degradation of pollutants. Their capacity to duplicate is dependent on their replication determinants (replicon), which also define their bacterial host range and the inability to coexist with related replicons. We characterize a second replicon from the virulence plasmid pPsv48C, from Pseudomonas syringae pv. savastanoi, which appears to be a natural chimera between the gene encoding a newly described replication protein and a putative replication control region present in the widespread family of PFP virulence plasmids. We present extensive evidence of this type of chimerism in structurally similar replicons from species of Pseudomonas, including environmental bacteria as well as plant, animal and human pathogens. We establish that these replicons consist of two functional modules corresponding to putative control (REx-C module) and replication (REx-R module) regions. These modules are functionally separable, do not show specificity for each other, and are dynamically exchanged among replicons of four distinct plasmid families. Only the REx-C module displays strong incompatibility, which is overcome by a few nucleotide changes clustered in a stem-and-loop structure of a putative antisense RNA. Additionally, a REx-C module from pPsv48C conferred replication ability to a non-replicative chromosomal DNA region containing features associated to replicons. Thus, the organization of plasmid replicons as independent and exchangeable functional modules is likely facilitating rapid replicon evolution, fostering their diversification and survival, besides allowing the potential co-option of appropriate genes into novel replicons and the artificial construction of new replicon specificities.

Nestedness across biological scales

Biological networks pervade nature. They describe systems throughout all levels of biological organization, from molecules regulating metabolism to species interactions that shape ecosystem dynamics. The network thinking revealed recurrent organizational patterns in complex biological systems, such as the formation of semi-independent groups of connected elements (modularity) and non-random distributions of interactions among elements. Other structural patterns, such as nestedness, have been primarily assessed in ecological networks formed by two non-overlapping sets of elements; information on its occurrence on other levels of organization is lacking. Nestedness occurs when interactions of less connected elements form proper subsets of the interactions of more connected elements. Only recently these properties began to be appreciated in one-mode networks (where all elements can interact) which describe a much wider variety of biological phenomena. Here, we compute nestedness in a diverse collection of one-mode networked systems from six different levels of biological organization depicting gene and protein interactions, complex phenotypes, animal societies, metapopulations, food webs and vertebrate metacommunities. Our findings suggest that nestedness emerge independently of interaction type or biological scale and reveal that disparate systems can share nested organization features characterized by inclusive subsets of interacting elements with decreasing connectedness. We primarily explore the implications of a nested structure for each of these studied systems, then theorize on how nested networks are assembled. We hypothesize that nestedness emerges across scales due to processes that, although system-dependent, may share a general compromise between two features: specificity (the number of interactions the elements of the system can have) and affinity (how these elements can be connected to each other). Our findings suggesting occurrence of nestedness throughout biological scales can stimulate the debate on how pervasive nestedness may be in nature, while the theoretical emergent principles can aid further research on commonalities of biological networks.

Identify bilayer modules via pseudo-3D clustering: applications to miRNA-gene bilayer networks

Module identification is a frequently used approach for mining local structures with more significance in global networks. Recently, a wide variety of bilayer networks are emerging to characterize the more complex biological processes. In the light of special topological properties of bilayer networks and the accompanying challenges, there is yet no effective method aiming at bilayer module identification to probe the modular organizations from the more inspiring bilayer networks. To this end, we proposed the pseudo-3D clustering algorithm, which starts from extracting initial non-hierarchically organized modules and then iteratively deciphers the hierarchical organization of modules according to a bottom-up strategy. Specifically, a modularity function for bilayer modules was proposed to facilitate the algorithm reporting the optimal partition that gives the most accurate characterization of the bilayer network. Simulation studies demonstrated its robustness and outperformance against alternative competing methods. Specific applications to both the soybean and human miRNA-gene bilayer networks demonstrated that the pseudo-3D clustering algorithm successfully identified the overlapping, hierarchically organized and highly cohesive bilayer modules. The analyses on topology, functional and human disease enrichment and the bilayer subnetwork involved in soybean fat biosynthesis provided both the theoretical and biological evidence supporting the effectiveness and robustness of pseudo-3D clustering algorithm.

Maintaining Genome Stability in Defiance of Mitotic DNA Damage

The implementation of decisions affecting cell viability and proliferation is based on prompt detection of the issue to be addressed, formulation and transmission of a correct set of instructions and fidelity in the execution of orders. While the first and the last are purely mechanical processes relying on the faithful functioning of single proteins or macromolecular complexes (sensors and effectors), information is the real cue, with signal amplitude, duration, and frequency ultimately determining the type of response. The cellular response to DNA damage is no exception to the rule. In this review article we focus on DNA damage responses in G2 and Mitosis. First, we set the stage describing mitosis and the machineries in charge of assembling the apparatus responsible for chromosome alignment and segregation as well as the inputs that control its function (checkpoints). Next, we examine the type of issues that a cell approaching mitosis might face, presenting the impact of post-translational modifications (PTMs) on the correct and timely functioning of pathways correcting errors or damage before chromosome segregation. We conclude this essay with a perspective on the current status of mitotic signaling pathway inhibitors and their potential use in cancer therapy.

Quantitative Determination of Flexible Pharmacological Mechanisms Based On Topological Variation in Mice Anti-Ischemic Modular Networks

Targeting modules or signalings may open a new path to understanding the complex pharmacological mechanisms of reversing disease processes. However, determining how to quantify the structural alteration of these signalings or modules in pharmacological networks poses a great challenge towards realizing rational drug use in clinical medicine. Here, we explore a novel approach for dynamic comparative and quantitative analysis of the topological structural variation of modules in molecular networks, proposing the concept of allosteric modules (AMs). Based on the ischemic brain of mice, we optimize module distribution in different compound-dependent modular networks by using the minimum entropy criterion and then calculate the variation in similarity values of AMs under various conditions using a novel method of SimiNEF. The diverse pharmacological dynamic stereo-scrolls of AMs with functional gradient alteration, which consist of five types of AMs, may robustly deconstruct modular networks under the same ischemic conditions. The concept of AMs can not only integrate the responsive mechanisms of different compounds based on topological cascading variation but also obtain valuable structural information about disease and pharmacological networks beyond pathway analysis. We thereby provide a new systemic quantitative strategy for rationally determining pharmacological mechanisms of altered modular networks based on topological variation.

The Evolutionary Origins of Hierarchy

Hierarchical organization-the recursive composition of sub-modules-is ubiquitous in biological networks, including neural, metabolic, ecological, and genetic regulatory networks, and in human-made systems, such as large organizations and the Internet. To date, most research on hierarchy in networks has been limited to quantifying this property. However, an open, important question in evolutionary biology is why hierarchical organization evolves in the first place. It has recently been shown that modularity evolves because of the presence of a cost for network connections. Here we investigate whether such connection costs also tend to cause a hierarchical organization of such modules. In computational simulations, we find that networks without a connection cost do not evolve to be hierarchical, even when the task has a hierarchical structure. However, with a connection cost, networks evolve to be both modular and hierarchical, and these networks exhibit higher overall performance and evolvability (i.e. faster adaptation to new environments). Additional analyses confirm that hierarchy independently improves adaptability after controlling for modularity. Overall, our results suggest that the same force-the cost of connections-promotes the evolution of both hierarchy and modularity, and that these properties are important drivers of network performance and adaptability. In addition to shedding light on the emergence of hierarchy across the many domains in which it appears, these findings will also accelerate future research into evolving more complex, intelligent computational brains in the fields of artificial intelligence and robotics.

Quantitative assessment of gene expression network module-validation methods

Validation of pluripotent modules in diverse networks holds enormous potential for systems biology and network pharmacology. An arising challenge is how to assess the accuracy of discovering all potential modules from multi-omic networks and validating their architectural characteristics based on innovative computational methods beyond function enrichment and biological validation. To display the framework progress in this domain, we systematically divided the existing Computational Validation Approaches based on Modular Architecture (CVAMA) into topology-based approaches (TBA) and statistics-based approaches (SBA). We compared the available module validation methods based on 11 gene expression datasets, and partially consistent results in the form of homogeneous models were obtained with each individual approach, whereas discrepant contradictory results were found between TBA and SBA. The TBA of the Zsummary value had a higher Validation Success Ratio (VSR) (51%) and a higher Fluctuation Ratio (FR) (80.92%), whereas the SBA of the approximately unbiased (AU) p-value had a lower VSR (12.3%) and a lower FR (45.84%). The Gray area simulated study revealed a consistent result for these two models and indicated a lower Variation Ratio (VR) (8.10%) of TBA at 6 simulated levels. Despite facing many novel challenges and evidence limitations, CVAMA may offer novel insights into modular networks.

Retroactivity in the Context of Modularly Structured Biomolecular Systems

Synthetic biology has intensively promoted the technical implementation of modular strategies in the fabrication of biological devices. Modules are considered as networks of reactions. The behavior displayed by biomolecular systems results from the information processes carried out by the interconnection of the involved modules. However, in natural systems, module wiring is not a free-of-charge process; as a consequence of interconnection, a reactive phenomenon called retroactivity emerges. This phenomenon is characterized by signals that propagate from downstream modules (the modules that receive the incoming signals upon interconnection) to upstream ones (the modules that send the signals upon interconnection). Such retroactivity signals, depending of their strength, may change and sometimes even disrupt the behavior of modular biomolecular systems. Thus, analysis of retroactivity effects in natural biological and biosynthetic systems is crucial to achieve a deeper understanding of how this interconnection between functionally characterized modules takes place and how it impacts the overall behavior of the involved cell. By discussing the modules interconnection in natural and synthetic biomolecular systems, we propose that such systems should be considered as quasi-modular.

Diverse array-designed modes of combination therapies in Fangjiomics

In line with the complexity of disease networks, diverse combination therapies have been demonstrated potential in the treatment of different patients with complex diseases in a personal combination profile. However, the identification of rational, compatible and effective drug combinations remains an ongoing challenge. Based on a holistic theory integrated with reductionism, Fangjiomics systematically develops multiple modes of array-designed combination therapies. We define diverse "magic shotgun" vertical, horizontal, focusing, siege and dynamic arrays according to different spatiotemporal distributions of hits on targets, pathways and networks. Through these multiple adaptive modes for treating complex diseases, Fangjiomics may help to identify rational drug combinations with synergistic or additive efficacy but reduced adverse side effects that reverse complex diseases by reconstructing or rewiring multiple targets, pathways and networks. Such a novel paradigm for combination therapies may allow us to achieve more precise treatments by developing phenotype-driven quantitative multi-scale modeling for rational drug combinations.

Evolution of bow-tie architectures in biology

Bow-tie or hourglass structure is a common architectural feature found in many biological systems. A bow-tie in a multi-layered structure occurs when intermediate layers have much fewer components than the input and output layers. Examples include metabolism where a handful of building blocks mediate between multiple input nutrients and multiple output biomass components, and signaling networks where information from numerous receptor types passes through a small set of signaling pathways to regulate multiple output genes. Little is known, however, about how bow-tie architectures evolve. Here, we address the evolution of bow-tie architectures using simulations of multi-layered systems evolving to fulfill a given input-output goal. We find that bow-ties spontaneously evolve when the information in the evolutionary goal can be compressed. Mathematically speaking, bow-ties evolve when the rank of the input-output matrix describing the evolutionary goal is deficient. The maximal compression possible (the rank of the goal) determines the size of the narrowest part of the network—that is the bow-tie. A further requirement is that a process is active to reduce the number of links in the network, such as product-rule mutations, otherwise a non-bow-tie solution is found in the evolutionary simulations. This offers a mechanism to understand a common architectural principle of biological systems, and a way to quantitate the effective rank of the goals under which they evolved.

Many biological systems show bow-tie (also called hourglass) architecture. A bow-tie means that a large number of inputs are converted to a small number of intermediates, which then fan out to generate a large number of outputs. For example, cells use a wide variety of nutrients; process them into 12 metabolic precursors, which are then used to make all of the cells biomass. Similar principles exist in biological signaling and in the information processing in the visual system. Despite the ubiquity of bow-tie structures in biology, there is no explanation of how they evolved. Here, we find that bow-ties spontaneously evolve when the information in the evolutionary goal they evolved to satisfy can be compressed. Mathematically, this means that the matrix representing the goal has deficient rank. The maximal compression possible determines the width of the bow-tie—the narrowest part in the network (equal to the rank of the goal matrix). This offers a mechanism to understand a common architectural principle of biological systems, and a way to quantitate the rank of the goals under which they evolved.

Convergence and divergence of genetic and modular networks between diabetes and breast cancer

Diabetes mellitus (DM) and breast cancer (BC) can simultaneously occur in the same patient populations, but the molecular relationship between them remains unknown. In this study, we constructed genetic networks and used modularized analysis approaches to investigate the multi-dimensional characteristics of two diseases and one disease subtype. A text search engine (Agilent Literature Search 2.71) and MCODE software were applied to validate potential subnetworks and to divide the modules, respectively. A total of 793 DM-related genes, 386 type 2 diabetes (T2DM) genes and 873 BC-related genes were identified from the Online Mendelian Inheritance in Man database. For DM and BC, a total of 99 overlapping genes, 9 modules, 29 biological processes and 7 pathways were identified. Meanwhile, for T2DM and BC, 56 overlapping genes, 5 modules, 20 biological processes and 12 pathways were identified. Based on the Gene Ontology functional enrichment analysis of the top 10 non-overlapping modules of the two diseases, 10 biological functions and 5 pathways overlapped between them. The glycosphingolipid and lysosome pathways verified molecular mechanisms of cell death related to both DM and BC. We also identified new biological functions of dopamine receptors and four signalling pathways (Parkinson's disease, Alzheimer's disease, Huntington's disease and long-term depression) related to both diseases; these warrant further investigation. Our results illustrate the landscape of the novel molecular substructures between DM and BC, which may support a new model for complex disease classification and rational therapies for multiple diseases.