General theory of genotype to phenotype mapping: derivation of epigenetic landscapes from N-node complex gene regulatory networks

Irreversibility in bacterial regulatory networks

Irreversibility, in which a transient perturbation leaves a system in a new state, is an emergent property in systems of interacting entities. This property has well-established implications in statistical physics but remains underexplored in biological networks, especially for bacteria and other prokaryotes whose regulation of gene expression occurs predominantly at the transcriptional level. Focusing on the reconstructed regulatory network of Escherichia coli, we examine network responses to transient single-gene perturbations. We predict irreversibility in numerous cases and find that the incidence of irreversibility increases with the proximity of the perturbed gene to positive circuits in the network. Comparison with experimental data suggests a connection between the predicted irreversibility to transient perturbations and the evolutionary response to permanent perturbations.

Stochastic Turing patterns of trichomes in Arabidopsis leaves

Biological patterns that emerge during the morphogenesis of multicellular organisms can display high precision at large scales, while at cellular scales, cells exhibit large fluctuations stemming from cell-cell differences in molecular copy numbers also called demographic noise. We study the conflicting interplay between high precision and demographic noise in trichome patterns on the epidermis of wild-type Arabidopsis thaliana leaves, as a two-dimensional model system. We carry out a statistical characterization of these patterns and show that their power spectra display fat tails-a signature compatible with noise-driven stochastic Turing patterns-which are absent in power spectra of patterns driven by deterministic instabilities. We then present a theoretical model that includes demographic noise stemming from birth-death processes of genetic regulators which we study analytically and by stochastic simulations. The model captures the observed experimental features of trichome patterns.

Uncoding the interdependency of tumor microenvironment and macrophage polarization: insights from a continuous network approach

The balance between pro- and anti-inflammatory immune system responses is crucial to preventing complex diseases like cancer. Macrophages are essential immune cells that contribute to this balance constrained by the local signaling profile of the tumor microenvironment. To understand how pro- and anti-inflammatory unbalance emerges in cancer, we developed a theoretical analysis of macrophage differentiation that is derived from activated monocytes circulating in the blood. Once recruited to the site of inflammation, monocytes can be polarized based on the specific interleukins and chemokines in the microenvironment. To quantify this process, we used a previous regulatory network reconstructed by our group and transformed Boolean Network attractors of macrophage polarization to an ODE scheme, it enables us to quantify the activation of their genes in a continuous fashion. The transformation was developed using the interaction rules with a fuzzy logic approach. By implementing this approach, we analyzed different aspects that cannot be visualized in the Boolean setting. For example, this approach allows us to explore the dynamic behavior at different concentrations of cytokines and transcription factors in the microenvironment. One important aspect to assess is the evaluation of the transitions between phenotypes, some of them characterized by an abrupt or a gradual transition depending on specific concentrations of exogenous cytokines in the tumor microenvironment. For instance, IL-10 can induce a hybrid state that transits between an M2c and an M2b macrophage. Interferon- γ can induce a hybrid between M1 and M1a macrophage. We further demonstrated the plasticity of macrophages based on a combination of cytokines and the existence of hybrid phenotypes or partial polarization. This mathematical model allows us to unravel the patterns of macrophage differentiation based on the competition of expression of transcriptional factors. Finally, we survey how macrophages may respond to a continuously changing immunological response in a tumor microenvironment.

Leveraging developmental landscapes for model selection in Boolean gene regulatory networks

Boolean models are a well-established framework to model developmental gene regulatory networks (DGRNs) for acquisition of cellular identities. During the reconstruction of Boolean DGRNs, even if the network structure is given, there is generally a large number of combinations of Boolean functions that will reproduce the different cell fates (biological attractors). Here we leverage the developmental landscape to enable model selection on such ensembles using the relative stability of the attractors. First we show that previously proposed measures of relative stability are strongly correlated and we stress the usefulness of the one that captures best the cell state transitions via the mean first passage time (MFPT) as it also allows the construction of a cellular lineage tree. A property of great computational importance is the insensitivity of the different stability measures to changes in noise intensities. That allows us to use stochastic approaches to estimate the MFPT and thereby scale up the computations to large networks. Given this methodology, we revisit different Boolean models of Arabidopsis thaliana root development, showing that a most recent one does not respect the biologically expected hierarchy of cell states based on relative stabilities. We therefore developed an iterative greedy algorithm that searches for models which satisfy the expected hierarchy of cell states and found that its application to the root development model yields many models that meet this expectation. Our methodology thus provides new tools that can enable reconstruction of more realistic and accurate Boolean models of DGRNs.

Modeling the effect of environmental cytokines, nutrient conditions and hypoxia on CD4+ T cell differentiation

Upon antigen stimulation and co-stimulation, CD4⁺ T lymphocytes produce soluble factors that promote the activity of other immune cells against pathogens or modified tissues; this task must be performed in presence of a variety of environmental cytokines, nutrient, and oxygen conditions, which necessarily impact T cell function. The complexity of the early intracellular processes taking place upon lymphocyte stimulation is addressed by means of a mathematical model based on a network that integrates variable microenvironmental conditions with intracellular activating, regulatory, and metabolic signals. Besides the phenotype subsets considered in previous works (Th1, Th2, Th17, and Treg) the model includes the main early events in differentiation to the T_FH phenotype. The model describes how cytokines, nutrients and oxygen availability regulate the differentiation of naïve CD4⁺ T cells into distinct subsets. Particularly, it shows that elevated amounts of an all-type mixture of effector cytokines under optimal nutrient and oxygen availability conduces the system towards a highly-polarized Th1 or Th2 state, while reduced cytokine levels allow the expression of the Th17, Treg or T_FH subsets, or even hybrid phenotypes. On the other hand, optimal levels of an all-type cytokine mixture in combination with glutamine or tryptophan restriction implies a shift from Th1 to Th2 expression, while decreased levels of the Th2-inducing cytokine IL-4 leads to the rupture of the Th1-Th2 axis, allowing the manifestation of different (or hybrid) subsets. Modeling proposes that, even under reduced levels of pro-inflammatory cytokines, the sole action of hypoxia boost Th17 expression.

The flowering transition pathways converge into a complex gene regulatory network that underlies the phase changes of the shoot apical meristem in Arabidopsis thaliana

Post-embryonic plant development is characterized by a period of vegetative growth during which a combination of intrinsic and extrinsic signals triggers the transition to the reproductive phase. To understand how different flowering inducing and repressing signals are associated with phase transitions of the Shoot Apical Meristem (SAM), we incorporated available data into a dynamic gene regulatory network model for Arabidopsis thaliana. This Flowering Transition Gene Regulatory Network (FT-GRN) formally constitutes a dynamic system-level mechanism based on more than three decades of experimental data on flowering. We provide novel experimental data on the regulatory interactions of one of its twenty-three components: a MADS-box transcription factor XAANTAL2 (XAL2). These data complement the information regarding flowering transition under short days and provides an example of the type of questions that can be addressed by the FT-GRN. The resulting FT-GRN is highly connected and integrates developmental, hormonal, and environmental signals that affect developmental transitions at the SAM. The FT-GRN is a dynamic multi-stable Boolean system, with 2²³ possible initial states, yet it converges into only 32 attractors. The latter are coherent with the expression profiles of the FT-GRN components that have been experimentally described for the developmental stages of the SAM. Furthermore, the attractors are also highly robust to initial states and to simulated perturbations of the interaction functions. The model recovered the meristem phenotypes of previously described single mutants. We also analyzed the attractors landscape that emerges from the postulated FT-GRN, uncovering which set of signals or components are critical for reproductive competence and the time-order transitions observed in the SAM. Finally, in the context of such GRN, the role of XAL2 under short-day conditions could be understood. Therefore, this model constitutes a robust biological module and the first multi-stable, dynamical systems biology mechanism that integrates the genetic flowering pathways to explain SAM phase transitions.

Continuous Modeling of T CD4 Lymphocyte Activation and Function

T CD4+ cells are central to the adaptive immune response against pathogens. Their activation is induced by the engagement of the T-cell receptor by antigens, and of co-stimulatory receptors by molecules also expressed on antigen presenting cells. Then, a complex network of intracellular events reinforce, diversify and regulate the initial signals, including dynamic metabolic processes that strongly influence both the activation state and the differentiation to effector cell phenotypes. The regulation of cell metabolism is controlled by the nutrient sensor adenosine monophosphate-activated protein kinase (AMPK), which drives the balance between oxidative phosphorylation (OXPHOS) and glycolysis. Herein, we put forward a 51-node continuous mathematical model that describes the temporal evolution of the early events of activation, integrating a circuit of metabolic regulation into the main routes of signaling. The model simulates the induction of anergy due to defective co-stimulation, the CTLA-4 checkpoint blockade, and the differentiation to effector phenotypes induced by external cytokines. It also describes the adjustment of the OXPHOS-glycolysis equilibrium by the action of AMPK as the effector function of the T cell develops. The development of a transient phase of increased OXPHOS before induction of a sustained glycolytic phase during differentiation to the Th1, Th2 and Th17 phenotypes is shown. In contrast, during Treg differentiation, glycolysis is subsequently reduced as cell metabolism is predominantly polarized towards OXPHOS. These observations are in agreement with experimental data suggesting that OXPHOS produces an ATP reservoir before glycolysis boosts the production of metabolites needed for protein synthesis, cell function, and growth.

Spatio-Temporal Dynamics of the Patterning of Arabidopsis Flower Meristem

The qualitative model presented in this work recovers the onset of the four fields that correspond to those of each floral organ whorl of Arabidopsis flower, suggesting a mechanism for the generation of the positional information required for the differential expression of the A, B, and C identity genes according to the ABC model for organ determination during early stages of flower development. Our model integrates a previous model for the emergence of WUS pattern in the floral meristem, and shows that this pre-pattern is a necessary but not sufficient condition for the posterior information of the four fields predicted by the ABC model. Furthermore, our model predicts that LFY diffusion along the L1 layer of cells is not a necessary condition for the patterning of the floral meristem.

Assessing the efficacy of eDNA metabarcoding for measuring microbial biodiversity within forest ecosystems

We investigated the nascent application and efficacy of sampling and sequencing environmental DNA (eDNA) in terrestrial environments using rainwater that filters through the forest canopy and understory vegetation (i.e., throughfall). We demonstrate the utility and potential of this method for measuring microbial communities and forest biodiversity. We collected pure rainwater (open sky) and throughfall, successfully extracted DNA, and generated over 5000 unique amplicon sequence variants. We found that several taxa including Mycoplasma sp., Spirosoma sp., Roseomonas sp., and Lactococcus sp. were present only in throughfall samples. Spiroplasma sp., Methylobacterium sp., Massilia sp., Pantoea sp., and Sphingomonas sp. were found in both types of samples, but more abundantly in throughfall than in rainwater. Throughfall samples contained Gammaproteobacteria that have been previously found to be plant-associated, and may contribute to important functional roles. We illustrate how this novel method can be used for measuring microbial biodiversity in forest ecosystems, foreshadowing the utility for quantifying both prokaryotic and eukaryotic lifeforms. Leveraging these methods will enhance our ability to detect extant species, describe new species, and improve our overall understanding of ecological community dynamics in forest ecosystems.

Dynamical modeling predicts an inflammation-inducible CXCR7+ B cell precursor with potential implications in lymphoid blockage pathologies

Background

The blockage at the early B lymphoid cell development pathway within the bone marrow is tightly associated with hematopoietic and immune diseases, where the disruption of basal regulatory networks prevents the continuous replenishment of functional B cells. Dynamic computational models may be instrumental for the comprehensive understanding of mechanisms underlying complex differentiation processes and provide novel prediction/intervention platforms to reinvigorate the system.

Methods

By reconstructing a three-module regulatory network including genetic transcription, intracellular transduction, and microenvironment communication, we have investigated the early B lineage cell fate decisions in normal and pathological settings. The early B cell differentiation network was simulated as a Boolean model and then transformed, using fuzzy logic, to a continuous model. We tested null and overexpression mutants to analyze the emergent behavior of the network. Due to its importance in inflammation, we investigated the effect of NFkB induction at different early B cell differentiation stages.

Results

While the exhaustive synchronous and asynchronous simulation of the early B cell regulatory network (eBCRN) reproduced the configurations of the hematopoietic progenitors and early B lymphoid precursors of the pathway, its simulation as a continuous model with fuzzy logics suggested a transient IL-7R⁺ ProB-to-Pre-B subset expressing pre-BCR and a series of dominant B-cell transcriptional factors. This conspicuous differentiating cell population up-regulated CXCR7 and reduced CXCR4 and FoxO1 expression levels. Strikingly, constant but intermediate NFkB signaling at specific B cell differentiation stages allowed stabilization of an aberrant CXCR7⁺ pre-B like phenotype with apparent affinity to proliferative signals, while under constitutive overactivation of NFkB, such cell phenotype was aberrantly exacerbated from the earliest stage of common lymphoid progenitors. Our mutant models revealed an abnormal delay in the BCR assembly upon NFkB activation, concomitant to sustained Flt3 signaling, down-regulation of Ebf1, Irf4 and Pax5 genes transcription, and reduced Ig recombination, pointing to a potential lineage commitment blockage.

Discussion

For the first time, an inducible CXCR7^hi B cell precursor endowed with the potential capability of shifting central lymphoid niches, is inferred by computational modeling. Its phenotype is compatible with that of leukemia-initiating cells and might be the foundation that bridges inflammation with blockage-related malignancies and a wide range of immunological diseases. Besides the predicted differentiation impairment, inflammation-inducible phenotypes open the possibility of newly formed niches colonized by the reported precursor. Thus, emergent bone marrow ecosystems are predicted following a pro-inflammatory induction, that may lead to hematopoietic instability associated to blockage pathologies.

Type 2 diabetes progression: A regulatory network approach

In order to elucidate central elements underlying type 2 diabetes, we constructed a regulatory network model involving 37 components (molecules, receptors, processes, etc.) associated to signaling pathways of pancreatic beta-cells. In a first approximation, the network topology was described by Boolean rules whose interacting dynamics predicted stationary patterns broadly classified as health, metabolic syndrome, and diabetes stages. A subsequent approximation based on a continuous logic analysis allowed us to characterize the progression of the disease as transitions between these states associated to alterations of cell homeostasis due to exhaustion or exacerbation of specific regulatory signals. The method allowed the identification of key transcription factors involved in metabolic stress as essential for the progression of the disease. Integration of the present analysis with existent mathematical models designed to yield accurate account of experimental data in human or animal essays leads to reliable predictions for beta-cell mass, insulinemia, glycemia, and glycosylated hemoglobin in diabetic fatty rats.

A system-level mechanistic explanation for asymmetric stem cell fates: Arabidopsis thaliana root niche as a study system

Asymmetric divisions maintain long-term stem cell populations while producing new cells that proliferate and then differentiate. Recent reports in animal systems show that divisions of stem cells can be uncoupled from their progeny differentiation, and the outcome of a division could be influenced by microenvironmental signals. But the underlying system-level mechanisms, and whether this dynamics also occur in plant stem cell niches (SCN), remain elusive. This article presents a cell fate regulatory network model that contributes to understanding such mechanism and identify critical cues for cell fate transitions in the root SCN. Novel computational and experimental results show that the transcriptional regulator SHR is critical for the most frequent asymmetric division previously described for quiescent centre stem cells. A multi-scale model of the root tip that simulated each cell's intracellular regulatory network, and the dynamics of SHR intercellular transport as a cell-cell coupling mechanism, was developed. It revealed that quiescent centre cell divisions produce two identical cells, that may acquire different fates depending on the feedback between SHR's availability and the state of the regulatory network. Novel experimental data presented here validates our model, which in turn, constitutes the first proposed systemic mechanism for uncoupled SCN cell division and differentiation.

From Discrete to Continuous Modeling of Lymphocyte Development and Plasticity in Chronic Diseases

The molecular events leading to differentiation, development, and plasticity of lymphoid cells have been subject of intense research due to their key roles in multiple pathologies, such as lymphoproliferative disorders, tumor growth maintenance and chronic diseases. The emergent roles of lymphoid cells and the use of high-throughput technologies have led to an extensive accumulation of experimental data allowing the reconstruction of gene regulatory networks (GRN) by integrating biochemical signals provided by the microenvironment with transcriptional modules of lineage-specific genes. Computational modeling of GRN has been useful for the identification of molecular switches involved in lymphoid specification, prediction of microenvironment-dependent cell plasticity, and analyses of signaling events occurring downstream the activation of antigen recognition receptors. Among most common modeling strategies to analyze the dynamical behavior of GRN, discrete dynamic models are widely used for their capacity to capture molecular interactions when a limited knowledge of kinetic parameters is present. However, they are less powerful when modeling complex systems sensitive to biochemical gradients. To compensate it, discrete models may be transformed into regulatory networks that includes state variables and parameters varying within a continuous range. This approach is based on a system of differential equations dynamics with regulatory interactions described by fuzzy logic propositions. Here, we discuss the applicability of this method on modeling of development and plasticity processes of adaptive lymphocytes, and its potential implications in the study of pathological landscapes associated to chronic diseases.

Medical Systems Biology

The aim of this volume is to encourage the use of systems-level methodologies to contribute to the improvement of human-health . We intend to motivate biomedical researchers to complement their current theoretical and empirical practice with up-to-date systems biology conceptual approaches. Our perspective is based on the deep understanding of the key biomolecular regulatory mechanisms that underlie health, as well as the emergence and progression of human-disease . We strongly believe that the contemporary systems biology perspective opens the door to the effective development of novel methodologies to the improvement of prevention . This requires a deeper and integrative understanding of the involved underlying systems-level mechanisms. In order to explain our proposal in a simple way, in this chapter we privilege the conceptual exposition of our chosen framework over formal considerations. The formal exposition of our proposal will be expanded and discussed later in the next chapters.

Role of Cytokine Combinations on CD4+ T Cell Differentiation, Partial Polarization, and Plasticity: Continuous Network Modeling Approach

Purpose: We put forward a theoretical and dynamical approach for the semi-quantitative analysis of CD4+ T cell differentiation, the process by which cells with different functions are derived from activated CD4+ T naïve lymphocytes in the presence of particular cytokine microenvironments. We explore the system-level mechanisms that underlie CD4+ T plasticity-the conversion of polarized cells to phenotypes different from those originally induced. Methods: In this paper, we extend a previous study based on a Boolean network to a continuous framework. The network includes transcription factors, signaling pathways, as well as autocrine and exogenous cytokines, with interaction rules derived using fuzzy logic. Results: This approach allows us to assess the effect of relative differences in the concentrations and combinations of exogenous and endogenous cytokines, as well as of the expression levels of diverse transcription factors. We found either abrupt or gradual differentiation patterns between observed phenotypes depending on critical concentrations of single or multiple environmental cytokines. Plastic changes induced by environmental cytokines were observed in conditions of partial phenotype polarization in the T helper 1 to T helper 2 transition. On the other hand, the T helper 17 to induced regulatory T-cells transition was highly dependent on cytokine concentrations, with TGFβ playing a prime role. Conclusion: The present approach is useful to further understand the system-level mechanisms underlying observed patterns of CD4+ T differentiation and response to changing immunological challenges.

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.

A physico-genetic module for the polarisation of auxin efflux carriers PIN-FORMED (PIN)

Intracellular polarisation of auxin efflux carriers is crucial for understanding how auxin gradients form in plants. The polarisation dynamics of auxin efflux carriers PIN-FORMED (PIN) depends on both biomechanical forces as well as chemical, molecular and genetic factors. Biomechanical forces have shown to affect the localisation of PIN transporters to the plasma membrane. We propose a physico-genetic module of PIN polarisation that integrates biomechanical, molecular, and cellular processes as well as their non-linear interactions. The module was implemented as a discrete Boolean model and then approximated to a continuous dynamic system, in order to explore the relative contribution of the factors mediating PIN polarisation at the scale of single cell. Our models recovered qualitative behaviours that have been experimentally observed and enable us to predict that, in the context of PIN polarisation, the effects of the mechanical forces can predominate over the activity of molecular factors such as the GTPase ROP6 and the ROP-INTERACTIVE CRIB MOTIF-CONTAINING PROTEIN RIC1.

Dynamic network modelling to understand flowering transition and floral patterning

Differentiation and morphogenetic processes during plant development are particularly robust. At the cellular level, however, plants also show great plasticity in response to environmental conditions, and can even reverse apparently terminal differentiated states with remarkable ease. Can we understand and predict both robust and plastic systemic responses as a general consequence of the non-trivial interplay between intracellular regulatory networks, extrinsic environmental signalling, and tissue-level mechanical constraints? Flower development has become an ideal model system to study these general questions of developmental biology, which are especially relevant to understanding stem cell patterning in plants, animals, and human disease. Decades of detailed study of molecular developmental genetics, as well as novel experimental techniques for in vivo assays in both wild-type and mutant plants, enable the postulation and testing of experimentally grounded mathematical and computational network dynamical models. Research in our group aims to explain the emergence of robust transitions that occur at the shoot apical meristem, as well as flower development, as the result of the collective action of key molecular components in regulatory networks subjected to intra-organismal signalling and extracellular constraints. Here we present a brief overview of recent work from our group, and that of others, focusing on the use of simple dynamical models to address cell-fate specification and cell-state stochastic dynamics during flowering transition and cell-state transitions at the shoot apical meristem of Arabidopsis thaliana. We also focus on how our work fits within the general field of plant developmental modelling, which is being developed by many others.

Relative stability of network states in Boolean network models of gene regulation in development

Progress in cell type reprogramming has revived the interest in Waddington's concept of the epigenetic landscape. Recently researchers developed the quasi-potential theory to represent the Waddington's landscape. The Quasi-potential U(x), derived from interactions in the gene regulatory network (GRN) of a cell, quantifies the relative stability of network states, which determine the effort required for state transitions in a multi-stable dynamical system. However, quasi-potential landscapes, originally developed for continuous systems, are not suitable for discrete-valued networks which are important tools to study complex systems. In this paper, we provide a framework to quantify the landscape for discrete Boolean networks (BNs). We apply our framework to study pancreas cell differentiation where an ensemble of BN models is considered based on the structure of a minimal GRN for pancreas development. We impose biologically motivated structural constraints (corresponding to specific type of Boolean functions) and dynamical constraints (corresponding to stable attractor states) to limit the space of BN models for pancreas development. In addition, we enforce a novel functional constraint corresponding to the relative ordering of attractor states in BN models to restrict the space of BN models to the biological relevant class. We find that BNs with canalyzing/sign-compatible Boolean functions best capture the dynamics of pancreas cell differentiation. This framework can also determine the genes' influence on cell state transitions, and thus can facilitate the rational design of cell reprogramming protocols.

A Minimal Regulatory Network of Extrinsic and Intrinsic Factors Recovers Observed Patterns of CD4+ T Cell Differentiation and Plasticity

CD4+ T cells orchestrate the adaptive immune response in vertebrates. While both experimental and modeling work has been conducted to understand the molecular genetic mechanisms involved in CD4+ T cell responses and fate attainment, the dynamic role of intrinsic (produced by CD4+ T lymphocytes) versus extrinsic (produced by other cells) components remains unclear, and the mechanistic and dynamic understanding of the plastic responses of these cells remains incomplete. In this work, we studied a regulatory network for the core transcription factors involved in CD4+ T cell-fate attainment. We first show that this core is not sufficient to recover common CD4+ T phenotypes. We thus postulate a minimal Boolean regulatory network model derived from a larger and more comprehensive network that is based on experimental data. The minimal network integrates transcriptional regulation, signaling pathways and the micro-environment. This network model recovers reported configurations of most of the characterized cell types (Th0, Th1, Th2, Th17, Tfh, Th9, iTreg, and Foxp3-independent T regulatory cells). This transcriptional-signaling regulatory network is robust and recovers mutant configurations that have been reported experimentally. Additionally, this model recovers many of the plasticity patterns documented for different T CD4+ cell types, as summarized in a cell-fate map. We tested the effects of various micro-environments and transient perturbations on such transitions among CD4+ T cell types. Interestingly, most cell-fate transitions were induced by transient activations, with the opposite behavior associated with transient inhibitions. Finally, we used a novel methodology was used to establish that T-bet, TGF-β and suppressors of cytokine signaling proteins are keys to recovering observed CD4+ T cell plastic responses. In conclusion, the observed CD4+ T cell-types and transition patterns emerge from the feedback between the intrinsic or intracellular regulatory core and the micro-environment. We discuss the broader use of this approach for other plastic systems and possible therapeutic interventions.

Reshaping the epigenetic landscape during early flower development: induction of attractor transitions by relative differences in gene decay rates

BACKGROUND

Gene regulatory network (GRN) dynamical models are standard systems biology tools for the mechanistic understanding of developmental processes and are enabling the formalization of the epigenetic landscape (EL) model.

METHODS

In this work we propose a modeling framework which integrates standard mathematical analyses to extend the simple GRN Boolean model in order to address questions regarding the impact of gene specific perturbations in cell-fate decisions during development.

RESULTS

We systematically tested the propensity of individual genes to produce qualitative changes to the EL induced by modification of gene characteristic decay rates reflecting the temporal dynamics of differentiation stimuli. By applying this approach to the flower specification GRN (FOS-GRN) we uncovered differences in the functional (dynamical) role of their genes. The observed dynamical behavior correlates with biological observables. We found a relationship between the propensity of undergoing attractor transitions between attraction basins in the EL and the direction of differentiation during early flower development - being less likely to induce up-stream attractor transitions as the course of development progresses. Our model also uncovered a potential mechanism at play during the transition from EL basins defining inflorescence meristem to those associated to flower organs meristem. Additionally, our analysis provided a mechanistic interpretation of the homeotic property of the ABC genes, being more likely to produce both an induced inter-attractor transition and to specify a novel attractor. Finally, we found that there is a close relationship between a gene's topological features and its propensity to produce attractor transitions.

CONCLUSIONS

The study of how the state-space associated with a dynamical model of a GRN can be restructured by modulation of genes' characteristic expression times is an important aid for understanding underlying mechanisms occurring during development. Our contribution offers a simple framework to approach such problem, as exemplified here by the case of flower development. Different GRN models and the effect of diverse inductive signals can be explored within the same framework. We speculate that the dynamical role of specific genes within a GRN, as uncovered here, might give information about which genes are more likely to link a module to other regulatory circuits and signaling transduction pathways.

XAANTAL2 (AGL14) Is an Important Component of the Complex Gene Regulatory Network that Underlies Arabidopsis Shoot Apical Meristem Transitions

In Arabidopsis thaliana, multiple genes involved in shoot apical meristem (SAM) transitions have been characterized, but the mechanisms required for the dynamic attainment of vegetative, inflorescence, and floral meristem (VM, IM, FM) cell fates during SAM transitions are not well understood. Here we show that a MADS-box gene, XAANTAL2 (XAL2/AGL14), is necessary and sufficient to induce flowering, and its regulation is important in FM maintenance and determinacy. xal2 mutants are late flowering, particularly under short-day (SD) condition, while XAL2 overexpressing plants are early flowering, but their flowers have vegetative traits. Interestingly, inflorescences of the latter plants have higher expression levels of LFY, AP1, and TFL1 than wild-type plants. In addition we found that XAL2 is able to bind the TFL1 regulatory regions. On the other hand, the basipetal carpels of the 35S::XAL2 lines lose determinacy and maintain high levels of WUS expression under SD condition. To provide a mechanistic explanation for the complex roles of XAL2 in SAM transitions and the apparently paradoxical phenotypes of XAL2 and other MADS-box (SOC1, AGL24) overexpressors, we conducted dynamic gene regulatory network (GRN) and epigenetic landscape modeling. We uncovered a GRN module that underlies VM, IM, and FM gene configurations and transition patterns in wild-type plants as well as loss and gain of function lines characterized here and previously. Our approach thus provides a novel mechanistic framework for understanding the complex basis of SAM development.

Modeling the epigenetic attractors landscape: toward a post-genomic mechanistic understanding of development

Robust temporal and spatial patterns of cell types emerge in the course of normal development in multicellular organisms. The onset of degenerative diseases may result from altered cell fate decisions that give rise to pathological phenotypes. Complex networks of genetic and non-genetic components underlie such normal and altered morphogenetic patterns. Here we focus on the networks of regulatory interactions involved in cell-fate decisions. Such networks modeled as dynamical non-linear systems attain particular stable configurations on gene activity that have been interpreted as cell-fate states. The network structure also restricts the most probable transition patterns among such states. The so-called Epigenetic Landscape (EL), originally proposed by C. H. Waddington, was an early attempt to conceptually explain the emergence of developmental choices as the result of intrinsic constraints (regulatory interactions) shaped during evolution. Thanks to the wealth of molecular genetic and genomic studies, we are now able to postulate gene regulatory networks (GRN) grounded on experimental data, and to derive EL models for specific cases. This, in turn, has motivated several mathematical and computational modeling approaches inspired by the EL concept, that may be useful tools to understand and predict cell-fate decisions and emerging patterns. In order to distinguish between the classical metaphorical EL proposal of Waddington, we refer to the Epigenetic Attractors Landscape (EAL), a proposal that is formally framed in the context of GRNs and dynamical systems theory. In this review we discuss recent EAL modeling strategies, their conceptual basis and their application in studying the emergence of both normal and pathological developmental processes. In addition, we discuss how model predictions can shed light into rational strategies for cell fate regulation, and we point to challenges ahead.

An extended gene protein/products Boolean network model including post-transcriptional regulation

Background

Networks Biology allows the study of complex interactions between biological systems using formal, well structured, and computationally friendly models. Several different network models can be created, depending on the type of interactions that need to be investigated. Gene Regulatory Networks (GRN) are an effective model commonly used to study the complex regulatory mechanisms of a cell. Unfortunately, given their intrinsic complexity and non discrete nature, the computational study of realistic-sized complex GRNs requires some abstractions. Boolean Networks (BNs), for example, are a reliable model that can be used to represent networks where the possible state of a node is a boolean value (0 or 1). Despite this strong simplification, BNs have been used to study both structural and dynamic properties of real as well as randomly generated GRNs.

Results

In this paper we show how it is possible to include the post-transcriptional regulation mechanism (a key process mediated by small non-coding RNA molecules like the miRNAs) into the BN model of a GRN. The enhanced BN model is implemented in a software toolkit (EBNT) that allows to analyze boolean GRNs from both a structural and a dynamic point of view. The open-source toolkit is compatible with available visualization tools like Cytoscape and allows to run detailed analysis of the network topology as well as of its attractors, trajectories, and state-space. In the paper, a small GRN built around the mTOR gene is used to demonstrate the main capabilities of the toolkit.

Conclusions

The extended model proposed in this paper opens new opportunities in the study of gene regulation. Several of the successful researches done with the support of BN to understand high-level characteristics of regulatory networks, can now be improved to better understand the role of post-transcriptional regulation for example as a network-wide noise-reduction or stabilization mechanisms.

The topology of robustness and evolvability in evolutionary systems with genotype-phenotype map

In this paper we formulate a topological definition of the concepts of robustness and evolvability. We start our investigation by formulating a multiscale model of the evolutionary dynamics of a population of cells. Our cells are characterised by a genotype-phenotype map: their chances of survival under selective pressure are determined by their phenotypes, whereas the latter are determined their genotypes. According to our multiscale dynamics, the population dynamics generates the evolution of a genotype-phenotype network. Our representation of the genotype-phenotype network is similar to previously described ones, but has a novel element, namely, our network contains two types of nodes: genotype and phenotype nodes. This network representation allows us to characterise robustness and evolvability in terms of its topological properties: phenotypic robustness by means of the clustering coefficient of the phenotype nodes, and evolvability as the emergence of giant connected component which allows navigation between phenotypes. This topological definition of evolvability allows us to characterise the so-called robustness of evolvability, which is defined in terms of the robustness against attack (i.e. edge removal) of the giant connected component. An investigation of the factors that affect the robustness of evolvability shows that phenotypic robustness and the cryptic genetic variation are key to the integrity of the ability to innovate. These results fit within the framework of a number of models which point out that robustness favours rather than hindering evolvability. We further show that the corresponding phenotype network, defined as the one-component projection of the whole genotype-phenotype network, exhibits the small-world phenomenon, which implies that in this type of evolutionary system the rate of adaptability is enhanced.

Principle for performing attractor transits with single control in Boolean networks

We present an algebraic approach to reveal attractor transitions in Boolean networks under single control based on the recently developed matrix semitensor product theory. In this setting, the reachability of attractors is estimated by the state transition matrices. We then propose procedures that compute the shortest control sequence and the result of each step of input (control) exactly. The general derivation is exemplified by numerical simulations for two kinds of gene regulation networks, the protein-nucleic acid interactions network and the cAMP receptor of Dictyostelium discoideum network.

Molecular evolution constraints in the floral organ specification gene regulatory network module across 18 angiosperm genomes

The gene regulatory network of floral organ cell fate specification of Arabidopsis thaliana is a robust developmental regulatory module. Although such finding was proposed to explain the overall conservation of floral organ types and organization among angiosperms, it has not been confirmed that the network components are conserved at the molecular level among flowering plants. Using the genomic data that have accumulated, we address the conservation of the genes involved in this network and the forces that have shaped its evolution during the divergence of angiosperms. We recovered the network gene homologs for 18 species of flowering plants spanning nine families. We found that all the genes are highly conserved with no evidence of positive selection. We studied the sequence conservation features of the genes in the context of their known biological function and the strength of the purifying selection acting upon them in relation to their placement within the network. Our results suggest an association between protein length and sequence conservation, evolutionary rates, and functional category. On the other hand, we found no significant correlation between the strength of purifying selection and gene placement. Our results confirm that the studied robust developmental regulatory module has been subjected to strong functional constraints. However, unlike previous studies, our results do not support the notion that network topology plays a major role in constraining evolutionary rates. We speculate that the dynamical functional role of genes within the network and not just its connectivity could play an important role in constraining evolution.