Spatial motifs reveal patterns in cellular architecture of complex tissues

Spatial organization of cells is crucial to both proper physiological function of tissues and pathological conditions like cancer. Recent advances in spatial transcriptomics have enabled joint profiling of gene expression and spatial context of the cells. The outcome is an information rich map of the tissue where individual cells, or small regions, can be labeled based on their gene expression state. While spatial transcriptomics excels in its capacity to profile numerous genes within the same sample, most existing methods for analysis of spatial data only examine distribution of one or two labels at a time. These approaches overlook the potential for identifying higher-order associations between cell types - associations that can play a pivotal role in understanding development and function of complex tissues. In this context, we introduce a novel method for detecting motifs in spatial neighborhood graphs. Each motif represents a spatial arrangement of cell types that occurs in the tissue more frequently than expected by chance. To identify spatial motifs, we developed an algorithm for uniform sampling of paths from neighborhood graphs and combined it with a motif finding algorithm on graphs inspired by previous methods for finding motifs in DNA sequences. Using synthetic data with known ground truth, we show that our method can identify spatial motifs with high accuracy and sensitivity. Applied to spatial maps of mouse retinal bipolar cells and hypothalamic preoptic region, our method reveals previously unrecognized patterns in cell type arrangements. In some cases, cells within these spatial patterns differ in their gene expression from other cells of the same type, providing insights into the functional significance of the spatial motifs. These results suggest that our method can illuminate the substantial complexity of neural tissues, provide novel insight even in well studied models, and generate experimentally testable hypotheses.

The Spatial Organization of Bacterial Transcriptional Regulatory Networks

The transcriptional regulatory network (TRN) is the central pivot of a prokaryotic organism to receive, process and respond to internal and external environmental information. However, little is known about its spatial organization so far. In recent years, chromatin interaction data of bacteria such as Escherichia coli and Bacillus subtilis have been published, making it possible to study the spatial organization of bacterial transcriptional regulatory networks. By combining TRNs and chromatin interaction data of E. coli and B. subtilis, we explored the spatial organization characteristics of bacterial TRNs in many aspects such as regulation directions (positive and negative), central nodes (hubs, bottlenecks), hierarchical levels (top, middle, bottom) and network motifs (feed-forward loops and single input modules) of the TRNs and found that the bacterial TRNs have a variety of stable spatial organization features under different physiological conditions that may be closely related with biological functions. Our findings provided new insights into the connection between transcriptional regulation and the spatial organization of chromosome in bacteria and might serve as a factual foundation for trying spatial-distance-based gene circuit design in synthetic biology.

Simplicial Persistence of Financial Markets: Filtering, Generative Processes and Structural Risk

We introduce simplicial persistence, a measure of time evolution of motifs in networks obtained from correlation filtering. We observe long memory in the evolution of structures, with a two power law decay regimes in the number of persistent simplicial complexes. Null models of the underlying time series are tested to investigate properties of the generative process and its evolutional constraints. Networks are generated with both a topological embedding network filtering technique called TMFG and by thresholding, showing that the TMFG method identifies high order structures throughout the market sample, where thresholding methods fail. The decay exponents of these long memory processes are used to characterise financial markets based on their efficiency and liquidity. We find that more liquid markets tend to have a slower persistence decay. This appears to be in contrast with the common understanding that efficient markets are more random. We argue that they are indeed less predictable for what concerns the dynamics of each single variable but they are more predictable for what concerns the collective evolution of the variables. This could imply higher fragility to systemic shocks.

Behavioural variation among workers promotes feed-forward loops in a simulated insect colony

Coordinated responses in eusocial insect colonies arise from worker interaction networks that enable collective processing of ecologically relevant information. Previous studies have detected a structural motif in these networks known as the feed-forward loop, which functions to process information in other biological regulatory networks (e.g. transcriptional networks). However, the processes that generate feed-forward loops among workers and the consequences for information flow within the colony remain largely unexplored. We constructed an agent-based model to investigate how individual variation in activity and movement shaped the production of feed-forward loops in a simulated insect colony. We hypothesized that individual variation along these axes would generate feed-forward loops by driving variation in interaction frequency among workers. We found that among-individual variation in activity drove over-representation of feed-forward loops in the interaction networks by determining the directionality of interactions. However, despite previous work linking feed-forward loops with efficient information transfer, activity variation did not promote faster or more efficient information flow, thus providing no support for the hypothesis that feed-forward loops reflect selection for enhanced collective functioning. Conversely, individual variation in movement trajectory, despite playing no role in generating feed-forward loops, promoted fast and efficient information flow by linking together otherwise unconnected regions of the nest.

Statistical Identification of Important Nodes in Biological Systems

Biological systems can be modeled and described by biological networks. Biological networks are typical complex networks with widely real-world applications. Many problems arising in biological systems can be boiled down to the identification of important nodes. For example, biomedical researchers frequently need to identify important genes that potentially leaded to disease phenotypes in animal and explore crucial genes that were responsible for stress responsiveness in plants. To facilitate the identification of important nodes in biological systems, one needs to know network structures or behavioral data of nodes (such as gene expression data). If network topology was known, various centrality measures can be developed to solve the problem; while if only behavioral data of nodes were given, some sophisticated statistical methods can be employed. This paper reviewed some of the recent works on statistical identification of important nodes in biological systems from three aspects, that is, 1) in general complex networks based on complex networks theory and epidemic dynamic models; 2) in biological networks based on network motifs; and 3) in plants based on RNA-seq data. The identification of important nodes in a complex system can be seen as a mapping from the system to the ranking score vector of nodes, such mapping is not necessarily with explicit form. The three aspects reflected three typical approaches on ranking nodes in biological systems and can be integrated into one general framework. This paper also proposed some challenges and future works on the related topics. The associated investigations have potential real-world applications in the control of biological systems, network medicine and new variety cultivation of crops.

Exact Probability Landscapes of Stochastic Phenotype Switching in Feed-Forward Loops: Phase Diagrams of Multimodality

Feed-forward loops (FFLs) are among the most ubiquitously found motifs of reaction networks in nature. However, little is known about their stochastic behavior and the variety of network phenotypes they can exhibit. In this study, we provide full characterizations of the properties of stochastic multimodality of FFLs, and how switching between different network phenotypes are controlled. We have computed the exact steady-state probability landscapes of all eight types of coherent and incoherent FFLs using the finite-butter Accurate Chemical Master Equation (ACME) algorithm, and quantified the exact topological features of their high-dimensional probability landscapes using persistent homology. Through analysis of the degree of multimodality for each of a set of 10,812 probability landscapes, where each landscape resides over 10⁵–10⁶ microstates, we have constructed comprehensive phase diagrams of all relevant behavior of FFL multimodality over broad ranges of input and regulation intensities, as well as different regimes of promoter binding dynamics. In addition, we have quantified the topological sensitivity of the multimodality of the landscapes to regulation intensities. Our results show that with slow binding and unbinding dynamics of transcription factor to promoter, FFLs exhibit strong stochastic behavior that is very different from what would be inferred from deterministic models. In addition, input intensity play major roles in the phenotypes of FFLs: At weak input intensity, FFL exhibit monomodality, but strong input intensity may result in up to 6 stable phenotypes. Furthermore, we found that gene duplication can enlarge stable regions of specific multimodalities and enrich the phenotypic diversity of FFL networks, providing means for cells toward better adaptation to changing environment. Our results are directly applicable to analysis of behavior of FFLs in biological processes such as stem cell differentiation and for design of synthetic networks when certain phenotypic behavior is desired.

Statistical Identification of Important Nodes in Biological Systems

Biological systems can be modeled and described by biological networks. Biological networks are typical complex networks with widely real-world applications. Many problems arising in biological systems can be boiled down to the identification of important nodes. For example, biomedical researchers frequently need to identify important genes that potentially leaded to disease phenotypes in animal and explore crucial genes that were responsible for stress responsiveness in plants. To facilitate the identification of important nodes in biological systems, one needs to know network structures or behavioral data of nodes (such as gene expression data). If network topology was known, various centrality measures can be developed to solve the problem; while if only behavioral data of nodes were given, some sophisticated statistical methods can be employed. This paper reviewed some of the recent works on statistical identification of important nodes in biological systems from three aspects, that is, 1) in general complex networks based on complex networks theory and epidemic dynamic models; 2) in biological networks based on network motifs; and 3) in plants based on RNA-seq data. The identification of important nodes in a complex system can be seen as a mapping from the system to the ranking score vector of nodes, such mapping is not necessarily with explicit form. The three aspects reflected three typical approaches on ranking nodes in biological systems and can be integrated into one general framework. This paper also proposed some challenges and future works on the related topics. The associated investigations have potential real-world applications in the control of biological systems, network medicine and new variety cultivation of crops.

The multiple meanings of omnivory influence empirical, modular theory and whole food web stability relationships

The persistence of whole communities hinges on the presence of select interactions which act to stabilize communities making the identification of these keystone interactions critical. One potential candidate is omnivory, yet theoretical research on omnivory thus far has been dominated by a modular theory approach whereby an omnivore and consumer compete for a shared resource. Empirical research, however, has highlighted the presence of a broader suite of omnivory modules. Here, we integrate empirical data analysis and mathematical models to explore the influence of both omnivory module (including classic, multi-resource, higher level, mutual predation and cannibalism) and omnivore-resource interaction type on food web stability. We use six classic empirical food webs to examine the prevalence of the different types of omnivory, a multi-species consumer-resource model to determine the stability of these different kinds of omnivory within a module context, and finally extend these models to a 50 species, whole food web model to examine the influence of omnivory on whole food web persistence. Our results challenge the concept that omnivory is broadly stabilizing. In particular, we demonstrate that the impact of omnivory depends on the type of omnivory being examined with multi-resource omnivory having the largest correlation with whole food web persistence. Moreover, our results highlight that we need to exercise caution when scaling modular theory to whole food web theory. Cannibalism, for example, was the most persistent and stable omnivory module in the modular theory analysis, but only demonstrated a weak correlation with whole food web persistence. Lastly, our results demonstrate that the frequency of omnivory modules are more important for whole food web persistence than the frequency of omnivore-resource interactions. Together, these results demonstrate that the role of omnivory often depends both on the type of omnivory being examined and the food web within which it is nested. In whole food web models, omnivory acts less as a keystone interaction, rather, specific types of omnivory, particularly multi-resource omnivory, act as keystone modules. Future work integrating module and whole food web theory is critical for resolving the role of key interactions in food webs.

Mechanistic Computational Models of MicroRNA-Mediated Signaling Networks in Human Diseases

MicroRNAs (miRs) are endogenous non-coding RNA molecules that play important roles in human health and disease by regulating gene expression and cellular processes. In recent years, with the increasing scientific knowledge and new discovery of miRs and their gene targets, as well as the plentiful experimental evidence that shows dysregulation of miRs in a wide variety of human diseases, the computational modeling approach has emerged as an effective tool to help researchers identify novel functional associations between differential miR expression and diseases, dissect the phenotypic expression patterns of miRs in gene regulatory networks, and elucidate the critical roles of miRs in the modulation of disease pathways from mechanistic and quantitative perspectives. Here we will review the recent systems biology studies that employed different kinetic modeling techniques to provide mechanistic insights relating to the regulatory function and therapeutic potential of miRs in human diseases. Some of the key computational aspects to be discussed in detail in this review include (i) models of miR-mediated network motifs in the regulation of gene expression, (ii) models of miR biogenesis and miR–target interactions, and (iii) the incorporation of such models into complex disease pathways in order to generate mechanistic, molecular- and systems-level understanding of pathophysiology. Other related bioinformatics tools such as computational platforms that predict miR-disease associations will also be discussed, and we will provide perspectives on the challenges and opportunities in the future development and translational application of data-driven systems biology models that involve miRs and their regulatory pathways in human diseases.

Motif discovery in biological network using expansion tree

Networks are powerful representation of topological features in biological systems like protein interaction and gene regulation. In order to understand the design principles of such complex networks, the concept of network motifs emerged. Network motifs are recurrent patterns with statistical significance that can be seen as basic building blocks of complex networks. Identification of network motifs leads to many important applications, such as understanding the modularity and the large-scale structure of biological networks, classification of networks into super-families, protein function annotation, etc. However, identification of network motifs is challenging as it involves graph isomorphism which is computationally hard. Though this problem has been studied extensively in the literature using different computational approaches, we are far from satisfactory results. Motivated by the challenges involved in this field, an efficient and scalable network Motif Discovery algorithm based on Expansion Tree (MODET) is proposed. Pattern growth approach is used in this proposed motif-centric algorithm. Each node of the expansion tree represents a non-isomorphic pattern. The embeddings corresponding to a child node of the expansion tree are obtained from the embeddings of the parent node through vertex addition and edge addition. Further, the proposed algorithm does not involve any graph isomorphism check and the time complexities of these processes are <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>O</mml:mi><mml:mo>(</mml:mo><mml:mi>n</mml:mi><mml:mo>)</mml:mo></mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>O</mml:mi><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:math> , respectively. The proposed algorithm has been tested on Protein-Protein Interaction (PPI) network obtained from the MINT database. The computational efficiency of the proposed algorithm outperforms most of the existing network motif discovery algorithms.

Burkitt lymphoma-associated network construction and important network motif analysis

Biological and medical researchers have discovered numerous transcription factors (TFs), microRNAs (miRNAs) and genes associated with Burkitt lymphoma (BL) through individual experiments; however, their regulatory mechanisms remain unclear. In the present study, BL-dysregulated and BL-associated networks were constructed to investigate these mechanisms. All data and regulatory associations were from known data resources and literature. The dysregulated network consisted of dysregulated TFs, miRNAs and genes, and partially determined the pathogenesis mechanisms underlying BL. The BL-associated network consisted of BL-associated TFs, miRNAs and genes. It has been indicated that the network motif consisted of TFs, miRNAs and genes serve potential functions in numerous biological processes within cancer. Two of the most studied network motifs are feedback loop (FBL) and feed-forward loop (FFL). The important network motifs were extracted, including the FBL motif, 3-nodes FFL motif and 4-nodes motif, from BL-dysregulated and BL-associated networks, and 10 types of motifs were identified from BL-associated network. Finally, 26/31 FBL motifs, 45/75 3-nodes FFL motifs and 54/94 4-nodes motifs were obtained from the dysregulated/associated networks. A total of four TFs (E2F1, NFKB1, E2F4 and TCF3) exhibit complicated regulation associations in BL-associated networks. The biological network does not demonstrate the dysregulated status for healthy people. When the individual becomes unwell, their biological network exhibits a dysregulated status. If the dysregulated status is regulated to a normal status by a number of medical methods, the diseases may be treated successfully. BL-dysregulated networks serve important roles in pathogenesis mechanisms underlying BL regulation of the dysregulated network, which may be an effective strategy that contributes to gene therapy for BL.

A Robust Molecular Network Motif for Period-Doubling Devices

Life is sustained by a variety of cyclic processes such as cell division, muscle contraction, and neuron firing. The periodic signals powering these processes often direct a variety of other downstream systems, which operate at different time scales and must have the capacity to divide or multiply the period of the master clock. Period modulation is also an important challenge in synthetic molecular systems, where slow and fast components may have to be coordinated simultaneously by a single oscillator whose frequency is often difficult to tune. Circuits that can multiply the period of a clock signal (frequency dividers), such as binary counters and flip-flops, are commonly encountered in electronic systems, but design principles to obtain similar devices in biological systems are still unclear. We take inspiration from the architecture of electronic flip-flops, and we propose to build biomolecular period-doubling networks by combining a bistable switch with negative feedback modules that preprocess the circuit inputs. We identify a network motif and we show it can be "realized" using different biomolecular components; two of the realizations we propose rely on transcriptional gene networks and one on nucleic acid strand displacement systems. We examine the capacity of each realization to perform period-doubling by studying how bistability of the motif is affected by the presence of the input; for this purpose, we employ mathematical tools from algebraic geometry that provide us with valuable insights on the input/output behavior as a function of the realization parameters. We show that transcriptional network realizations operate correctly also in a stochastic regime when processing oscillations from the repressilator, a canonical synthetic in vivo oscillator. Finally, we compare the performance of different realizations in a range of realistic parameters via numerical sensitivity analysis of the period-doubling region, computed with respect to the input period and amplitude. Our mathematical and computational analysis suggests that the motif we propose is generally robust with respect to specific implementation details: functionally equivalent circuits can be built as long as the species-interaction topology is respected. This indicates that experimental construction of the circuit is possible with a variety of components within the rapidly expanding libraries available in synthetic biology.

Construction of a lncRNA-mediated feed-forward loop network reveals global topological features and prognostic motifs in human cancers

Long non-coding RNAs (lncRNAs), transcription factors and microRNAs can form lncRNA-mediated feed-forward loops (L-FFLs), which are functional network motifs that regulate a wide range of biological processes, such as development and carcinogenesis. However, L-FFL network motifs have not been systematically identified, and their roles in human cancers are largely unknown. In this study, we computationally integrated data from multiple sources to construct a global L-FFL network for six types of human cancer and characterized the topological features of the network. Our approach revealed several dysregulated L-FFL motifs common across different cancers or specific to particular cancers. We also found that L-FFL motifs can take part in other types of regulatory networks, such as mRNA-mediated FFLs and ceRNA networks, and form the more complex networks in human cancers. In addition, survival analyses further indicated that L-FFL motifs could potentially serve as prognostic biomarkers. Collectively, this study elucidated the roles of L-FFL motifs in human cancers, which could be beneficial for understanding cancer pathogenesis and treatment.

Detection of Network Motif Based on a Novel Graph Canonization Algorithm from Transcriptional Regulation Networks

Network motifs are patterns of complex networks occurring significantly more frequently than those in random networks. They have been considered as fundamental building blocks of complex networks. Therefore, the detection of network motifs in transcriptional regulation networks is a crucial step in understanding the mechanism of transcriptional regulation and network evolution. The search for network motifs is similar to solving subgraph searching problems, which has proven to be NP-complete. To quickly and effectively count subgraphs of a large biological network, we propose a novel graph canonization algorithm based on resolving sets. This method has been implemented in a command line interface (CLI) program sgip using the SeqAn library. Comparing to Babai's algorithm, this approach has a tighter complexity bound, o ( exp ( n log 2 n + 4 log n ) ) , on strongly regular graphs. Results on several simulated datasets and transcriptional regulation networks indicate that sgip outperforms nauty on many graph cases. The source code of sgip is freely accessible in https://github.com/seqan/seqan/tree/master/apps/sgip and the binary code in http://packages.seqan.de/sgip/.

Form and function in gene regulatory networks: the structure of network motifs determines fundamental properties of their dynamical state space

Network motifs have been studied extensively over the past decade, and certain motifs, such as the feed-forward loop, play an important role in regulatory networks. Recent studies have used Boolean network motifs to explore the link between form and function in gene regulatory networks and have found that the structure of a motif does not strongly determine its function, if this is defined in terms of the gene expression patterns the motif can produce. Here, we offer a different, higher-level definition of the ‘function’ of a motif, in terms of two fundamental properties of its dynamical state space as a Boolean network. One is the basin entropy, which is a complexity measure of the dynamics of Boolean networks. The other is the diversity of cyclic attractor lengths that a given motif can produce. Using these two measures, we examine all 104 topologically distinct three-node motifs and show that the structural properties of a motif, such as the presence of feedback loops and feed-forward loops, predict fundamental characteristics of its dynamical state space, which in turn determine aspects of its functional versatility. We also show that these higher-level properties have a direct bearing on real regulatory networks, as both basin entropy and cycle length diversity show a close correspondence with the prevalence, in neural and genetic regulatory networks, of the 13 connected motifs without self-interactions that have been studied extensively in the literature.

Modular Organization of Residue-Level Contacts Shapes the Selection Pressure on Individual Amino Acid Sites of Ribosomal Proteins

Understanding the molecular evolution of macromolecular complexes in the light of their structure, assembly, and stability is of central importance. Here, we address how the modular organization of native molecular contacts shapes the selection pressure on individual residue sites of ribosomal complexes. The bacterial ribosomal complex is represented as a residue contact network where nodes represent amino acid/nucleotide residues and edges represent their van der Waals interactions. We find statistically overrepresented native amino acid-nucleotide contacts (OaantC, one amino acid contacts one or multiple nucleotides, internucleotide contacts are disregarded). Contact number is defined as the number of nucleotides contacted. Involvement of individual amino acids in OaantCs with smaller contact numbers is more random, whereas only a few amino acids significantly contribute to OaantCs with higher contact numbers. An investigation of structure, stability, and assembly of bacterial ribosome depicts the involvement of these OaantCs in diverse biophysical interactions stabilizing the complex, including high-affinity protein-RNA contacts, interprotein cooperativity, intersubunit bridge, packing of multiple ribosomal RNA domains, etc. Amino acid-nucleotide constituents of OaantCs with higher contact numbers are generally associated with significantly slower substitution rates compared with that of OaantCs with smaller contact numbers. This evolutionary rate heterogeneity emerges from the strong purifying selection pressure that conserves the respective amino acid physicochemical properties relevant to the stabilizing interaction with OaantC nucleotides. An analysis of relative molecular orientations of OaantC residues and their interaction energetics provides the biophysical ground of purifying selection conserving OaantC amino acid physicochemical properties.

The free energy cost of accurate biochemical oscillations

Oscillation is an important cellular process that regulates timing of different vital life cycles. However, in the noisy cellular environment, oscillations can be highly inaccurate due to phase fluctuations. It remains poorly understood how biochemical circuits suppress phase fluctuations and what is the incurred thermodynamic cost. Here, we study three different types of biochemical oscillations representing three basic oscillation motifs shared by all known oscillatory systems. In all the systems studied, we find that the phase diffusion constant depends on the free energy dissipation per period following the same inverse relation parameterized by system specific constants. This relationship and its range of validity are shown analytically in a model of noisy oscillation. Microscopically, we find that the oscillation is driven by multiple irreversible cycles that hydrolyze the fuel molecules such as ATP; the number of phase coherent periods is proportional to the free energy consumed per period. Experimental evidence in support of this general relationship and testable predictions are also presented.

Probabilistic generation of random networks taking into account information on motifs occurrence

Because of the huge number of graphs possible even with a small number of nodes, inference on network structure is known to be a challenging problem. Generating large random directed graphs with prescribed probabilities of occurrences of some meaningful patterns (motifs) is also difficult. We show how to generate such random graphs according to a formal probabilistic representation, using fast Markov chain Monte Carlo methods to sample them. As an illustration, we generate realistic graphs with several hundred nodes mimicking a gene transcription interaction network in Escherichia coli.

Network Analysis of Cancer-focused Association Network Reveals Distinct Network Association Patterns

Cancer is a complex and heterogeneous disease. Genetic methods have uncovered thousands of complex tissue-specific mutation-induced effects and identified multiple disease gene targets. Important associations between cancer and other biological entities (eg, genes and drugs) in cancer network, however, are usually scattered in biomedical publications. Systematic analyses of these cancer-specific associations can help highlight the hidden associations between different cancer types and related genes/drugs. In this paper, we proposed a novel network-based computational framework to identify statistically over-expressed subnetwork patterns called network motifs (NMs) in an integrated cancer-specific drug-disease-gene network extracted from Semantic MEDLINE, a database containing extracted associations from MEDLINE abstracts. Eight significant NMs were identified and considered as the backbone of the cancer association network. Each NM corresponds to specific biological meanings. We demonstrated that such approaches will facilitate the formulization of novel cancer research hypotheses, which is critical for translational medicine research and personalized medicine in cancer.

Validation of Aintegumenta as a gene to modify floral size in ornamental plants

The gene AINTEGUMENTA (AtANT) is an APETALA2 transcription factor in Arabidopsis activating growth downstream of auxin signalling. Lateral organ size is positively correlated with ANT expression in Arabidopsis. We tested the use of AtANT as a tool to modify floral size in two different plants used as model organisms and ornamental crops, Petunia × hybrida and Antirrhinum majus. Petunia plants expressing PhANT RNAi showed a decrease in PhANT expression correlated with smaller petal limbs. In contrast Petunia plants overexpressing AtANT had larger petal limbs. Petal tube length was less affected in down-regulation of PhANT or overexpression of AtANT. Overexpression of AtANT in Antirrhinum caused increased flower size via increased petal limb width and tube length. Down-regulation of PhANT showed an effect on cell size while overexpression of AtANT in Petunia and Antirrhinum caused significant increases in cell expansion that could explain the differences in floral organ size. The endogenous expression levels of PhANT and AmANT tended to be higher in the limb than in the tube in both Antirrhinum and Petunia. AtANT overexpression caused significant AmANT up-regulation in Antirrhinum limbs but not of PhANT in Petunia, indicating differences in the regulatory network. The differential effect of AtANT on limb and tube in Petunia and Antirrhinum correspond to phenotypic differences observed in natural variation in the corresponding genus indicating a relation between the phenotypic space of a genus and the effect of modified ANT levels, validating ANT as a gene to modify floral size.

Robust development of synfire chains from multiple plasticity mechanisms

Biological neural networks are shaped by a large number of plasticity mechanisms operating at different time scales. How these mechanisms work together to sculpt such networks into effective information processing circuits is still poorly understood. Here we study the spontaneous development of synfire chains in a self-organizing recurrent neural network (SORN) model that combines a number of different plasticity mechanisms including spike-timing-dependent plasticity, structural plasticity, as well as homeostatic forms of plasticity. We find that the network develops an abundance of feed-forward motifs giving rise to synfire chains. The chains develop into ring-like structures, which we refer to as "synfire rings." These rings emerge spontaneously in the SORN network and allow for stable propagation of activity on a fast time scale. A single network can contain multiple non-overlapping rings suppressing each other. On a slower time scale activity switches from one synfire ring to another maintaining firing rate homeostasis. Overall, our results show how the interaction of multiple plasticity mechanisms might give rise to the robust formation of synfire chains in biological neural networks.

Current innovations and future challenges of network motif detection

Network motif detection is the search for statistically overrepresented subgraphs present in a larger target network. They are thought to represent key structure and control mechanisms. Although the problem is exponential in nature, several algorithms and tools have been developed for efficiently detecting network motifs. This work analyzes 11 network motif detection tools and algorithms. Detailed comparisons and insightful directions for using these tools and algorithms are discussed. Key aspects of network motif detection are investigated. Network motif types and common network motifs as well as their biological functions are discussed. Applications of network motifs are also presented. Finally, the challenges, future improvements and future research directions for network motif detection are also discussed.

Regulatory network features in Listeria monocytogenes-changing the way we talk

Our understanding of how pathogens shape their gene expression profiles in response to environmental changes is ever growing. Advances in Bioinformatics have made it possible to model complex systems and integrate data from variable sources into one large regulatory network. In these analyses, regulatory networks are typically broken down into regulatory motifs such as feed-forward loops (FFL) or auto-regulatory feedbacks, which serves to simplify the structure, while the functional implications of different regulatory motifs allow to make informed assumptions about the function of a specific regulatory pathway. Here we review the basic concepts of network features and use this language to break down the regulatory networks that govern the interactions between the main regulators of stress response, virulence, and transmission in Listeria monocytogenes. We point out the advantage that taking a “systems approach” could have for our understanding of gene functions, the detection of distant regulatory inputs, interspecies comparisons, and co-expression.

Noise attenuation in the ON and OFF states of biological switches

Biological switches must sense changes in signal concentration and at the same time buffer against signal noise. While many studies have focused on the response of switching systems to noise in the ON state, how systems buffer noise at both ON and OFF states is poorly understood. Through analytical and computational approaches, we find that switching systems require different dynamics at the OFF state than at the ON state in order to have good noise buffering capability. Specifically, we introduce a quantity called the input-associated Signed Activation Time (iSAT) that concisely captures an intrinsic temporal property at either the ON or OFF state. We discover a trade-off between achieving good noise buffering in the ON versus the OFF states: a large iSAT corresponds to noise amplification in the OFF state in contrast to noise buffering in the ON state. To search for biological circuits that can buffer noise in both ON and OFF states, we systematically analyze all three-node circuits and identify mutual activation as a central motif. We also study connections among signal sensitivity, iSAT, and noise amplification. We find that a large iSAT at the ON state maintains signaling sensitivity while minimizing noise propagation. Taken together, the analysis of iSATs helps reveal the noise properties of biological networks and should aid in the design of robust switches that can both repress noise at the OFF state and maintain a reliable ON state.

BioFNet: biological functional network database for analysis and synthesis of biological systems

In synthetic biology and systems biology, a bottom-up approach can be used to construct a complex, modular, hierarchical structure of biological networks. To analyze or design such networks, it is critical to understand the relationship between network structure and function, the mechanism through which biological parts or biomolecules are assembled into building blocks or functional networks. A functional network is defined as a subnetwork of biomolecules that performs a particular function. Understanding the mechanism of building functional networks would help develop a methodology for analyzing the structure of large-scale networks and design a robust biological circuit to perform a target function. We propose a biological functional network database, named BioFNet, which can cover the whole cell at the level of molecular interactions. The BioFNet takes an advantage in implementing the simulation program for the mathematical models of the functional networks, visualizing the simulated results. It presents a sound basis for rational design of biochemical networks and for understanding how functional networks are assembled to create complex high-level functions, which would reveal design principles underlying molecular architectures.

The interplay between transcription factors and microRNAs in genome-scale regulatory networks

Metazoan genomes contain thousands of protein-coding and non-coding RNA genes, most of which are differentially expressed, i.e., at different locations, at different times during development, or in response to environmental signals. Differential gene expression is achieved through complex regulatory networks that are controlled in part by two types of trans-regulators: transcription factors (TFs) and microRNAs (miRNAs). TFs bind to cis-regulatory DNA elements that are often located in or near their target genes, while miRNAs hybridize to cis-regulatory RNA elements mostly located in the 3' untranslated region of their target mRNAs. Here, we describe how these trans-regulators interact with each other in the context of gene regulatory networks to coordinate gene expression at the genome-scale level, and discuss future challenges of integrating these networks with other types of functional networks.

Building trophic modules into a persistent food web

Understanding what maintains species and perpetuates their coexistence in a network of feeding relationships (the food web) is of great importance for biodiversity conservation. A food web can be viewed as consisting of a number of simple subunits called trophic modules. Intraguild predation (IGP), in which a prey and its predator compete for the same resource, is one of the best-studied trophic modules. According to theory, there are two ways to yield a large persistent system from such modules: (i) to use persistent subunits as building blocks or (ii) to arrange the subunits in a way that externally supports the nonpersistent subunits. Here, I show that the complex food web of the Caribbean marine ecosystem is constructed in both ways. I show that IGP modules, which convey internal persistence because of the fact that prey are superior competitors for the resources, are overrepresented in the Caribbean ecosystem. The other modules, consisting of competitively inferior prey, are not persistent in isolation. However, competitively inferior prey in these modules tend to receive more advantage from extra-module interactions, which allows persistence of the IGP module. In addition, those exterior interactions tend to be provided by intrinsically persistent IGP modules to prevent cascading extinction of interacting IGP modules. The food web can be viewed as a set of interacting modules, nonrandomly arranged to enhance the maintenance of biodiversity.

Centrality analysis methods for biological networks and their application to gene regulatory networks

The structural analysis of biological networks includes the ranking of the vertices based on the connection structure of a network. To support this analysis we discuss centrality measures which indicate the importance of vertices, and demonstrate their applicability on a gene regulatory network. We show that common centrality measures result in different valuations of the vertices and that novel measures tailored to specific biological investigations are useful for the analysis of biological networks, in particular gene regulatory networks.

Combinatorics of feedback in cellular uptake and metabolism of small molecules

We analyze the connection between structure and function for regulatory motifs associated with cellular uptake and usage of small molecules. Based on the boolean logic of the feedback we suggest four classes: the socialist, consumer, fashion, and collector motifs. We find that the socialist motif is good for homeostasis of a useful but potentially poisonous molecule, whereas the consumer motif is optimal for nutrition molecules. Accordingly, examples of these motifs are found in, respectively, the iron homeostasis system in various organisms and in the uptake of sugar molecules in bacteria. The remaining two motifs have no obvious analogs in small molecule regulation, but we illustrate their behavior using analogies to fashion and obesity. These extreme motifs could inspire construction of synthetic systems that exhibit bistable, history-dependent states, and homeostasis of flux (rather than concentration).