article: Alternative cell differentiation pathways are believed to arise from the concerted action of signalling pathways and transcriptional regulatory networks. However, the prediction of mammalian cell differentiation from the knowledge of the presence of specific signals and transcriptional factors is still a daunting challenge. In this respect, the vertebrate hematopoietic system, with its many branching differentiation pathways and cell types, is a compelling case study. In this paper, we propose an integrated, comprehensive model of the regulatory network and signalling pathways controlling Th cell differentiation. As most available data are qualitative, we rely on a logical formalism to perform extensive dynamical analyses. To cope with the size and complexity of the resulting network, we use an original model reduction approach together with a stable state identification algorithm. To assess the effects of heterogeneous environments on Th cell differentiation, we have performed a systematic series of simulations considering various prototypic environments. Consequently, we have identified stable states corresponding to canonical Th1, Th2, Th17 and Treg subtypes, but these were found to coexist with other transient hybrid cell types that co-express combinations of Th1, Th2, Treg and Th17 markers in an environment-dependent fashion. In the process, our logical analysis highlights the nature of these cell types and their relationships with canonical Th subtypes. Finally, our logical model can be used to explore novel differentiation pathways in silico.
  Alternative cell differentiation pathways are believed to arise from the concerted action of signalling pathways and transcriptional regulatory networks. However, the prediction of mammalian cell differentiation from the knowledge of the presence of specific signals and transcriptional factors is currently a daunting challenge. In this respect, the vertebrate hematopoietic system, with its many branching differentiation pathways and cell types, is a compelling case study. In particular, numerous publications describe molecular and genetic interactions involved in the control of the late stages of (TCR+, CD4+) T helper (Th) cell differentiation, and yet our current knowledge on cell lineage branching is clearly incomplete, as demonstrated by recent reports on novel Th cell types [1]–[4]. The existence of some of these cell types was revealed following the identification of novel transcription factors [2] or cytokines [1], [3]. However, several Th cell phenotypes recently described likely result from reassorted expression of already known genes [4], [5]. Beyond the expression of diverse and cell-specific antigen receptor genes, the appreciation of the heterogeneity of late Th cell lineages emerged from the characterisation of Th1 and Th2 cell types [6]. Both cell types arise following the sustained activation of uncommitted Th0 cell precursors and can be characterised by the expression of mutually exclusive sets of cytokines: Th1 cells produce IFN-, whereas Th2 cells express IL-4, IL-5 and IL-6. These cytokine profiles have a critical influence on the selection of a specific immune response, driving pro-inflammatory or allergic responses, and promoting alternative antibody classes. The cellular dichotomy has been mechanistically explained by the mutual inhibition between the master transcription factors T-bet and GATA-3 at the single cell level, as well as by cross-regulatory mechanisms at the cell population level [7]. Indeed, the cytokines produced by each Th subtype drives the differentiation of precursors into the same pathway while inhibiting alternative pathways. Additional T-helper subtypes have been recently identified. Regulatory T cells that depend on the transcription factor Foxp3 are capable of preventing (auto) immunity by inhibiting the activation and proliferation of other cells [8]. Furthermore, pro-inflammatory Th17 cells expressing IL-17 and dependent on RORt have been characterised. Current evidences indicate that late Th cell differentiation pathways are more complex and likely comprise further, non-canonical cell types [5], [9]–[12], whose mechanistic underpinnings and functional roles remain to be established. Effective immunity to many fungi and bacteria requires that the T cell response is dominated by pro-inflammatory effector Th1 or Th17 cells. Allergic reactions, whether beneficial or deleterious, are strictly dependent on Th2 cells. Avoiding spontaneous autoimmunity or controlling the collateral damage of effective immune responses to infection involves a fine balance between regulatory T cells and other Th cells. Genetic defects or accidental failures affecting this delicate balance can lead to irreversible immunopathology. The issue is not only how heterogeneous are these cell types, but, perhaps more importantly, how is the heterogeneity of Th cells sustained and controlled? How does antigen-specific memory correlate with the predominance of an appropriate Th cell branch? How plastic or resilient are the differentiated Th cells? During their lifetime, naive and memory Th cells face a changing environment in time and space. The lymphoid tissues are heteregeneous and provide variable local cytokine contexts to circulating Th cells. How robust are the Th subtypes with respect to these heterogeneous environments? Will a Th cell that differentiated into a Th1 phenotype in a lymph node always remain in this state or can it switch to another cell type if it faces a different environment? Evidences for substantial plasticity have been recently reported. For example, in vitro stimulation of Th2 cells in the presence of TGF- generates a non-canonical cell type expressing IL-10 and IL-9 in the absence of Foxp3 [5], [10]. Furthermore, Foxp3+ regulatory T cells loose the expression of this key transcription factor in the absence of effector T cells in vivo and can then drive inflammatory responses [13]. How many instances of such conversions should be expected? Are some Th cells irreversibly committed and others more plastic? It is difficult to address these questions directly due to the current impossibility to follow a single cell as it circulates in the body, during the rather long time scale of cell differentiation. Instead, studies are usually made on cell populations and therefore measure the predominance of one or another cell type. However, using mathematical modelling, these questions can be addressed in terms of stability and robustness of differentiation states of the molecular network underpinning the cellular phenotypes. Mathematical modelling has been used recurrently to formalise hypothetic regulatory schemes in immunology. Early phenomenological models represented the development of Th1 vs Th2 responses from naive Th0 cells using ordinary differential equations [14]–[17]. In these models, the subcellular molecular network controlling the state of the cell was implicit behind the transitions between cell types. These models accounted for the role of cell interactions in driving population commitment and sustained polarised responses. In general, they featured mutual inhibitions among cell populations, thereby enabling multistability. Alternative population stationary states were interpreted as polarised cell responses. However, such cell population models are unable to predict novel cell types or to question their plasticity because cell properties are hardwired in the model structure. More recently, models of the cellular networks driving Th cell differentiation and polarisation have been formulated using logical [18] (for an earlier logical model of T-cell regulation, see [19]) or ordinary differential equations [20], [21]. These models assume cross-inhibiting master transcription factors to generate canonical Th subtypes, thereby precluding cell plasticity. In this paper, we propose an integrated, comprehensive model of the regulatory network and signalling pathways accounting for the core control of Th cell differentiation. As most available data are qualitative, we rely on a qualitative, logical formalism to perform extensive dynamical analyses. To cope with the size and complexity of the resulting network, we use an original model reduction approach described in detail elsewhere [22]. To assess the effects of heterogeneous environments on Th cell differentiation, we have performed systematic simulations, considering various prototypic environments. As we shall see, stable states corresponding to canonical Th1, Th2, Th17 and Treg subtypes are readily identified, but they are found to coexist with other hybrid cell types that co-express combinations of Th1, Th2, Treg and Th17 markers in an environment-dependent fashion. In the process, our logical analysis highlights the nature of these cell types and their relationships to canonical Th subtypes.
  The precise roles of the different molecular species involved in the regulation of T cell differentiation are sparsely known. Even in the cases where direct regulatory interactions have been documented, little or no quantitative information is available on the relative strengths or rates of these processes. The extended logical formalism [23], [24] is a discrete modelling framework well adapted to biological systems where available information is qualitative. In this framework, a regulatory network is modelled in terms of a regulatory graph, where nodes represent regulatory components (proteins, complexes, genes, etc.), whereas arcs represent interactions between these components (i. e transcriptional activations or inhibitions, phosphorylations, etc.). In addition, each regulatory component is associated with a logical variable denoting its qualitative concentration or level of activity. In many cases, Boolean variables capture the most relevant situations (i. e. a Boolean variable takes the value 1 if the component is present or active, 0 otherwise). It is worth noting that components may represent phenomenological features besides specific molecular species (e. g. cell proliferation, see Table 1). Whenever needed (i. e when different levels of a component have distinct functional consequences), multi-valued variables are introduced. In our Th model, ternary variables have been associated with several interleukin receptor components (IL4RA, IL4R, IL2R, IL12RB1), as well as with STAT5, which can be up-regulated depending on signalling. Next, logical rules are defined for each regulatory component to specify its target activity level according to the levels of its regulators. Whereas some authors consider standard logical functions for all components, we do not impose such a restriction. For example, although STAT1 and IL21 both receive three activatory interactions, their logical rules differ: the rule assigned to STAT1 stipulates that STAT1 = 1 if IFNBR = 1 or IFNGR = 1 or IL27R = 1 (otherwise STAT1 = 0), whereas the rule assigned to IL21 stipulates that IL21 = 1 if NFAT = 1 and proliferation = 1 and STAT3 = 1. Given a regulatory graph and starting from a (set of) initial state (s), successor states can be recursively computed. This results in a new graph, called state transition graph, which describes the dynamical behaviour of the system. In this graph, a node denotes a state of the system (i. e. , a tuple giving the levels of the regulatory components), whereas an arc linking two states denotes a possible state transition. Considering a state, if the activity level of a component differs from the target level defined by its logical rule, there is an updating call on the corresponding variable. We generally assume that a state has as many successors as updated components (fully asynchronous dynamics), potentially leading to alternative behaviours. The computation of the state transition graph can be restricted by considering a set of priority classes. Each regulatory component is then associated to a priority class and is updated only in the absence of concurrent updating call with a higher priority [25]. Considering a state transition graph, it is particularly interesting to identify the attractors, which correspond to potential asymptotical behaviours, being either stable states (states without any successor) or cyclical components (denoting oscillatory behaviours). The asynchronous dynamical analysis of increasingly large regulatory graphs can be very challenging due to the exponential growth of the state transition graphs. A solution consists in reducing the model by removing (intermediate) components. We have recently proposed an algorithm to automatically compute the logical rules for a user-defined reduced model [22]. Starting with a detailed model, the computation of a reduced model is performed by iteratively removing regulatory components. Auto-regulated components are not entitled for removal, thereby avoiding the loss of dynamical properties associated with regulatory circuits (see below). Removing a regulatory component G implies a revised wiring where all targets of G are directly regulated by the regulators of G. The logical rules of the targets of G are modified accordingly to conserve the (indirect) effects of its regulators. This reduction method ensures the preservation of a number of dynamical properties of the original model. In particular, stable states and more complex attractors are conserved. However, additional cyclical attractors may arise from the isolation of transient cycles of the original system. An attractor which is reachable in the reduced model is also reachable in the full model, but the reverse is not always true, as the reduction may lead to the loss of some trajectories (see [22] for further details). One asset of the logical framework is the possibility to analyse the dynamical roles of the regulatory circuits (circular chains of interactions) embedded in a regulatory graph. The dynamical roles of isolated regulatory circuits depend on their signs. A circuit is positive if it encompasses an even number of inhibitions, whereas a negative circuit involves an odd number of inhibitions. During the eighties, R. Thomas conjectured that a positive circuit is required to generate multiple attractors, whereas a negative circuit is necessary for sustained oscillations [26]. Since then, these rules have inspired several theorems in different mathematical frameworks (see [27] for a brief survey). In this work, we focus on positive circuits, which enable the existence of alternative differentiated cell lineages, each corresponding to one attractor. While the presence of a positive circuit is necessary for the existence of multiple stable states, this condition is not sufficient. Indeed, when regulatory circuits are embedded in large networks, external regulators may prevent some circuits to generate the expected dynamical behaviour. Circuit functionality contexts can then be defined in terms of constraints on the values of the regulators of circuit members (see [28] for further details). In practice, in large regulatory graphs, a tiny fraction of the regulatory circuits are usually functional (i. e. have non-empty functionality contexts). Such circuit functionality analysis pinpoints crucial regulatory structures (usually with a predominance of short circuits) that are responsible for salient dynamical features. Developed in Java, the software GINsim provides a user interface to set up logical models, as well as to compute state transition graphs under diverse updating schemes, including the use of priority classes [29]. GINsim further implements an efficient algorithm for the determination of all stable states of a logical model, as well as the reduction method and the regulatory circuit analysis briefly mentioned above (see also [28]). We used GINsim to build a comprehensive logical model by integrating a large set of documented molecular actors and interactions. The resulting model is available in the model repository linked to the GINsim website. As GINsim allows the user to document a model by associating annotations (textual comments, links and references) with each component, extensive documentation has been integrated in the model file (see Dataset S1). In the case of the model presented here, we have to deal with many input nodes that collectively handle local environmental conditions. Varying the values of these inputs, all possible fates of a given cell can be computed in terms of attractors reachable in these diverse local environments. Our simulations aim at determining how naive cells can differentiate into specialised cells, depending on local environments, and how these differentiated cells respond to environmental changes. In this respect, we have defined a set of prototypic local environments (i. e. in terms of specific combinations of values for APC and external cytokines), that are biologically meaningful. For each of these environments, a simulation is performed starting from a state corresponding to naive Th cells. All the resulting attractors consist in stable states (i. e. attractors reduced to a unique state, each corresponding to a cell lineage). These stable states are in turn taken as initial states for a new round of simulations, considering different local environments. We have chained such simulations until no new stable state could be found, thereby allowing us to assess the behaviour of each cell lineage depending on the local environment. These simulations have been performed using priority classes assigning a higher rank to the receptors (TCR and cytokine receptors) and their downstream signal transducers and/or transcription factors (NFAT and the STATs) over other components. As we wish to take into account the roles of various external cytokines on the activation and differentiation of Th cells, we expect to find many possible stable states, some differing from each other only by minor component activities, and thus corresponding to the same Th subtype. Similar stable states can be grouped and denoted by the value of a vector encompassing the major transcription factors (GATA-3, T-bet, RORt, Foxp3) and cytokines expressed by specific Th subtypes (IFN-, TGF-, IL-2, etc.). The resulting patterns corresponding to different (quiescent, activated, or anergic) modalities of the same Th subtype can be further grouped into cell constellations, where simulated cellular state transitions are denoted by arcs linking the corresponding cell types.
  The logical formalism enables a modular approach for the construction of large regulatory network models, where parts of the network are defined and studied separately, before merging them to generate a comprehensive model. The Th cell network has been constructed by combining several “modules” illustrated in Figure 1 and described below. All the modules must then be integrated to build a comprehensive model. The TCR signalling module functions as an input, since it is not regulated by any other component. Several cytokine receptors share subchains and targets. For example, the common gamma chain (CGC) is shared by IL2, IL4, IL7, IL9 and IL15 receptors that lead to the activation of STAT5 (cf. Table 2). Furthermore, the cytokine modules are connected through a number of cross-regulatory interactions. Figure 2 displays the whole regulatory graph assembled for this study, while Table 1 provides a brief description of all the components included in the model (Dataset S1). The logical rules associated with each component are given in Text S1, along with biological justifications and bibliographical references. Which cell types and differentiation pathways can be predicted from the logical model just described? The different cell types correspond to attractors in the state transition graph, i. e. regions of the state transition graph from which the system cannot escape. Among these attractors, stable states can be readily determined [28]. Other attractors may consist in (intertwined) terminal cycles. Their identification requires a thorough analysis of the state transition graph. As our regulatory network encompasses too many components to enable a direct analysis of the full state transition graph, we have applied the reduction method described in Section “Model reduction”. Regarding the selection of the components for reduction, we face a compromise between computational performance and biological readability. Selecting TCR, CD28, cytokines receptors and their subchains, along with the intermediate components RUNX3, IRF1, SMAD3, IKB and NFKB for reduction, we obtain a regulatory graph containing 34 components (see Figure 3 and Dataset S2). Although the resulting logical model is still too large to compute the whole state transition graph (encompassing over states), it is now amenable to simulations starting from selected initial states. Indeed, in this regulatory graph, 13 out of 34 components represent environmental conditions (inputs), which are fixed for each simulation. With 21 internal components, the system has “only” about 4. 5 millions of possible states, most of which are not reachable when starting from a specific initial state. A further reduction leads to a graph conserving the 13 inputs, but only 12 internal nodes: TBET, GATA3, RORGT, FOXP3, STAT3, STAT5, IL2, IL4, IFNG, IL17, TGFB and proliferation. Using this compact model, we could compute the full state transition graphs for relevant input combinations. This led us to identify 28 context-dependent stable states (corresponding to the greyed-out cells in Figure 4) and to verify the absence of cyclic attractors for the polarising environments considered. Since our reduction preserves attractors [22], we can conclude that, for these environments, all attractors of the original model are indeed stable states. Furthermore, whenever a stable state can be reached in the compact model from specific initial conditions, we can conclude that it is also the case for the full model.
  We have reconstructed the molecular network controlling the activation and differentiation of Th cells and asked how many stable states to expect, considering that these cells face a changing local environment during their life span. Components and cross-regulatory links were extracted from the literature and, in some cases, from previous logical models [18], [34]. It is likely that components or links were missed. To this adds the problem of the definition of the logical functions driving the behaviour of the components, particularly those involving complex regulatory mechanisms. This cautionary remark notwithstanding, our current model recapitulates the differentiation of naive cells into Th1, Th2, Th17 and Treg subtypes. Strikingly, our model also gives rise to hybrid states expressing markers characteristic of two or more canonical cell types. Cyclic attractors (potentially corresponding to oscillatory behaviour) have been found only in restricted environmental situations, whose biological relevance remains to be assessed. Furthermore, our model analysis emphasises an unexpected plasticity of the canonical cell types. Indeed, according to the model, both Foxp3+ regulatory T cells and Th17 cells are highly plastic and labile, whereas Th1 and Th2 subtypes are more readily maintained across different environmental conditions. Based on our results, canonical Foxp3+ regulatory T cells would not be truly lineage committed, but would rather correspond to a context-dependent stable state of the underlying regulatory network. This is surprising considering the fundamental role of these cells in avoiding autoimmunity. According to our model analysis, the maintenance of a regulatory T cell phenotype would require sustained TCR/CD28 and IL-2 signals. Indeed, the absence of TCR stimulation can lead to the loss of Foxp3 expression and a conversion into a Th0 phenotype. Depending on inputs, these cells can regain Foxp3 expression but also differentiate into other cell types. Our model further predicts that Treg cells may differentiate into Th1 or Th2 subtypes in proper polarising environments. This plasticity of regulatory Foxp3+ T cells is supported by several recent reports [13], [35], [36] showing that fluorescence sorted Foxp3+ T cells loose Foxp3 expression and their suppressive capacity, under experimental conditions consistent with our simulations (e. g. loss of Foxp3 expression after adoptive transfer of regulatory T cells in the absence of conventional Foxp3- CD4 T cells, which provide IL-2 to the former [13]). However, in addition to cells with labile Foxp3 expression, the regulatory Foxp3+ T pool apparently also contains irreversibly committed cells, perhaps of thymic origin [35], which cannot be easily explained with our model. In our model analysis, Treg lability clearly depends on the assumption that Foxp3 expression requires Stat5 and AP-1/NFAT transcriptional activities. This might turn out to be an oversimplification as we might have overlooked some mechanisms preventing regulatory cells to become conventional T cells and vice-versa. Some data suggest that locus epigenetic regulation is needed for sustainable Foxp3 expression [37], and that regulatory miRNAs are likely involved [38]. Alternatively, the irreversible regulatory T cell state might be explained by an additional positive loop involving a transcription factor upstream of Foxp3. According to this scenario, Foxp3 could merely control anergy and suppressive functions of regulatory T cells, but would not be the Treg master regulator. The recent characterisation of Th cells with partial Treg phenotype in Foxp3-deficient mice supports this scenario [39]. Yet another alternative is that higher order regulatory mechanisms stemming from intercellular interactions and population dynamics may stabilise the Treg phenotype. Indeed, differential equation models [40] suggest that stable regulatory T cells pools could be maintained despite a continuous interconversion of cells. Further evidence is provided by loss or gain of Treg phenotype depending on the remaining populations [13]. Now, if our prediction of pervasive T cell plasticity is correct, how can the classes of immune responses in which they predominate be so robust? This point is particularly intriguing in the context of natural tolerance, presumably depending on labile Tregs, or yet regarding persistent memory responses to bacterial infections. Here again, the stability of the responses could stem from regulatory feedbacks at the level of the T cell populations, as intercellular interactions and cell population dynamics might sustain and stabilise specific combinations of Th subtypes. Regarding Treg mediated tolerance, stable and robust coexistence of conventional Th cells and regulatory T cells has been obtained using a cross-regulatory model encompassing positive and negative feedback circuits at the population level in [41]. In stationary populations, the production of IL-2 by conventional T cells sustains the expression of Foxp3 by Tregs. Disruption of this balance can lead to the conversion of Tregs into conventional Th phenotypes, as observed in adoptive transfers [13]. Regarding Th cell memory, the recent finding that memory cells relocate to the bone marrow, where they remain quiescent in special niches [42], raises the possibility that stromal cells may provide the environment necessary to sustain a context-dependent memory state. Among the model stable states, those corresponding to hybrid cell types expressing two or more Th master regulators are particularly striking. Most of these hybrid cell types co-express RORt with another Th master regulator. The remaining hybrids combine the expression of Foxp3 with other master regulators. The generation of these hybrid stable states might reveal missing cross-inhibitory mechanisms between RORt and Foxp3, or yet among these and the two other master regulators. The introduction of such cross-inhibitory circuits would grant mutually exclusive expression of master regulators (as assumed for GATA-3 and T-bet and generalised in [21]). Some indications of additional mutual inhibitions among master regulators can be found in the recent literature (e. g. between Foxp3 and RORt [37], [43]). Alternatively, since predicted hybrid cell types tend to be anergic (i. e. do not produce cytokines), they could have been overlooked in routine experimental assays characterising cell populations based on cytokine expression. In fact, the existence of some of these hybrid phenotypes is supported by recent studies relying on more sophisticated quantitative assays, thereby suggesting that mutually exclusive master regulator expression might not be a general principle of the Th regulatory network architecture. For example, hybrid cells expressing both IFN- and IL-17, dependent on T-bet and RORt respectively, have been observed by flow cytometry [44] or clonal analysis [45], thus resembling our in silico Th1 RORt+ subtype (row 7 in Figure 6). Other reports point to the generation of hybrid cells similar to our Treg RORt+ subtype (row 15 in Figure 6), following the induction of RORt and IL-17 expression in regulatory Foxp3+ T cells, both in mice [46], [47] and humans [48]. Furthermore, in the presence of IFN-, Foxp3+ regulatory T cells have been shown to up-regulate T-bet [12], leading to hybrid cells akin our Th1 Foxp3+ subtype (row 8 in Figure 6). Finally, genetically enforcing GATA-3 expression does not prevent differentiation of Th17 cells, leading to GATA-3+ IL-17 expressing hybrid cells [49]. Whether these cells are an artefact or reflect naturally occurring activated Th2 RORt+ cells (row 13 in Figure 6) remains to be established. Altogether, these observations support the existence of several of our predicted environment-dependent Th subtypes, thereby warranting the investigation of the other predicted hybrids (e. g. Th1 Foxp3+ RORt+ or Th2 Foxp3+ RORt+). Indeed, we have identified several environment-dependent Foxp3+ subtypes expressing other master transcription factors characteristic of canonical Th1, Th2 or Th17 cells. What are the functional implications of these chimeric regulatory cells? Since Foxp3 inhibits the expression of the cytokines downstream of GATA-3, T-bet, or RORt, these hybrid cells are expected to be anergic. In an environment dominated by other effector Th cells, up-regulation of a subset of genes not suppressed by Foxp3 could allow chimeric cells to mingle within effector T cells, profiting from local growth or survival factors, as suggested by population dynamical modelling [41]. This hypothesis is supported by the recent observation that essential Th1 or Th2 regulatory factors are required in Foxp3+ cells to enable proper control of Th1 or Th2 responses, respectively [12], [50]. The number of reported cytokines grows rapidly. Tentatively, some of these cytokines could be predominantly expressed by specific, yet undiscovered, cell types. One might thus wonder what part of our conclusions could be retained as our knowledge on the Th cell differentiation network grows. In this respect, we have extended earlier models for Th1/Th2 bipolarisation [18], [34], with the aim to account for two additional cell fates (Treg and Th17 subtypes). Interestingly, although our model accounts for many more cell fates, Th0, Th1 and Th2 subtypes are maintained as reference subtypes, around which are organised the cellular constellations displayed in Figure 8. This organisation can be related to the definition of a robust GATA-3 and T-bet cross-regulatory module [51] progressively embedded in a more comprehensive regulatory network. However, even this well established paradigm has been recently challenged by reports demonstrating simultaneous and sustained GATA-3 and T-bet co-expression in reprogrammed Th2 cells [52]. Although mechanistic details still need to be worked out, this finding suggests that Th cell plasticity may be even greater than fostered by our model. Further refinements of the model (e. g. the introduction of an additional level for T-bet, and potentially of additional components) will be required to properly model the generation of such hybrid Th2+1 cells. More generally, the characterisation of the roles of additional cytokines, transcription factors, and regulatory RNAs, or yet the delineation of epigenetic regulatory mechanisms, should enable proper model extensions and refinements, thereby leading to refined predictions, while preserving the functional roles of the most salient regulatory modules. In conclusion, our results indicate that the pool of CD4 Th cells is highly heterogeneous and that the structure of the gene network promotes this heterogeneity. Diversity of antigen receptors and their crossreactivity are the most salient features of the adaptive immune system of the vertebrates, allowing the host to recognize and potentially react to numerous antigens. The heterogeneity and plasticity of Th cell types emphasised here might further contribute to the capacity of the immune system to deal with wide contingencies. Arguably, such plasticity does not fit well with the classical depiction of T cell differentiation potential in terms of a branching tree. Instead, our computational study points to a reticulate network of alternative, environment-dependent, differentiation and reprogramming events.
summary: T lymphocytes play a key role in the regulation of the immune response in mammals. Various T-helper subtypes (Th1, Th2, Th17, Treg, …) have been identified over the years, characterised by the expression of specific transcription factors and cytokines, which have a critical influence on the selection of different immune responses, driving pro-inflammatory or allergic responses, promoting alternative antibody classes, or preventing (auto) immunity by inhibiting the activation and proliferation of other cells. To gain insight into the heterogeneity and the plasticity of late T-helper lineages, we have built an integrated model of the regulatory network and signalling pathways controlling Th cell differentiation. Relying on a logical modelling framework, we have performed a systematic series of simulations to assess the effects of heterogeneous environments on Th cell differentiation. We have identified stable states corresponding to canonical Th1, Th2, Th17 and Treg subtypes, but also to hybrid cell types co-expressing combinations of Th1, Th2, Treg and Th17 markers in an environment-dependent fashion. Our analysis highlights the nature of these cell types and their relationships with canonical Th subtypes.
keywords: immunology/leukocyte development
  computational biology/transcriptional regulation
  immunology/immune response
  computational biology/signaling networks
  immunology/leukocyte signaling and gene expression
  computational biology/systems biology
title: Diversity and Plasticity of Th Cell Types Predicted from Regulatory Network Modelling