Denis Thieffry

GOToolBox: functional analysis of gene datasets based on Gene Ontology

We have developed methods and tools based on the Gene Ontology (GO) resource allowing the identification of statistically over- or under-represented terms in a gene dataset; the clustering of functionally related genes within a set; and the retrieval of genes sharing annotations with a query gene. GO annotations can also be constrained to a slim hierarchy or a given level of the ontology. The source codes are available upon request, and distributed under the GPL license.

RSAT 2011: regulatory sequence analysis tools

RSAT (Regulatory Sequence Analysis Tools) comprises a wide collection of modular tools for the detection of cis-regulatory elements in genome sequences. Thirteen new programs have been added to the 30 described in the 2008 NAR Web Software Issue, including an automated sequence retrieval from EnsEMBL (retrieve-ensembl-seq), two novel motif discovery algorithms (oligo-diff and info-gibbs), a 100-times faster version of matrix-scan enabling the scanning of genome-scale sequence sets, and a series of facilities for random model generation and statistical evaluation (random-genome-fragments, random-motifs, random-sites, implant-sites, sequence-probability, permute-matrix). Our most recent work also focused on motif comparison (compare-matrices) and evaluation of motif quality (matrix-quality) by combining theoretical and empirical measures to assess the predictive capability of position-specific scoring matrices. To process large collections of peak sequences obtained from ChIP-seq or related technologies, RSAT provides a new program (peak-motifs) that combines several efficient motif discovery algorithms to predict transcription factor binding motifs, match them against motif databases and predict their binding sites. Availability (web site, stand-alone programs and SOAP/WSDL (Simple Object Access Protocol/Web Services Description Language) web services): http://rsat.ulb.ac.be/rsat/.

Modeling ERBB receptor-regulated G1/S transition to find novel targets for de novo trastuzumab resistance

BackgroundIn breast cancer, overexpression of the transmembrane tyrosine kinase ERBB2 is an adverse prognostic marker, and occurs in almost 30% of the patients. For therapeutic intervention, ERBB2 is targeted by monoclonal antibody trastuzumab in adjuvant settings; however, de novo resistance to this antibody is still a serious issue, requiring the identification of additional targets to overcome resistance. In this study, we have combined computational simulations, experimental testing of simulation results, and finally reverse engineering of a protein interaction network to define potential therapeutic strategies for de novo trastuzumab resistant breast cancer.ResultsFirst, we employed Boolean logic to model regulatory interactions and simulated single and multiple protein loss-of-functions. Then, our simulation results were tested experimentally by producing single and double knockdowns of the network components and measuring their effects on G1/S transition during cell cycle progression. Combinatorial targeting of ERBB2 and EGFR did not affect the response to trastuzumab in de novo resistant cells, which might be due to decoupling of receptor activation and cell cycle progression. Furthermore, examination of c-MYC in resistant as well as in sensitive cell lines, using a specific chemical inhibitor of c-MYC (alone or in combination with trastuzumab), demonstrated that both trastuzumab sensitive and resistant cells responded to c-MYC perturbation.ConclusionIn this study, we connected ERBB signaling with G1/S transition of the cell cycle via two major cell signaling pathways and two key transcription factors, to model an interaction network that allows for the identification of novel targets in the treatment of trastuzumab resistant breast cancer. Applying this new strategy, we found that, in contrast to trastuzumab sensitive breast cancer cells, combinatorial targeting of ERBB receptors or of key signaling intermediates does not have potential for treatment of de novo trastuzumab resistant cells. Instead, c-MYC was identified as a novel potential target protein in breast cancer cells.

Dynamical behaviour of biological regulatory networks—II. Immunity control in bacteriophage lambda*

A number of bacterial and viral genes take part in the decision between lysis and lysogenization in temperate bacteriophages. In the lambda case, at least five viral genes (cI, cro, cII, N and cIII) and several bacterial genes are involved. Several attempts have been made to model this complex regulatory network. Our approach is based on a logical method described in the first paper of the series which formalizes the interactions between the elements of a regulatory network in terms of discrete variables, functions and parameters. In this paper two models are described and discussed, the first (two-variable model) focused on cI and cro interactions, the second (four-variable model) considering, in addition, genes cII and N. The treatment presented emphasizes the roles of positive and negative feedback loops and their interactions in the development of the phage. The role of the loops between cI and cro, and of cI on itself (which both have to be positive loops) was discovered earlier; this groups contribution to this aspect mainly deals with the possibility of treating these loops as parts of a more extended network. In contrast, the role of the negative loop of cro on itself had apparently remained unexplained. We realized that this loop buffers the expression of genes cro itself. cII, O and P against the inflation due to the rapid replication of the phage. More generally, negative auto-control of a gene appears an efficient way to render its expression insensitive (or less sensitive) to gene dosage, whereas a simple negative control would not provide this result.

RSAT peak-motifs: motif analysis in full-size ChIP-seq datasets

ChIP-seq is increasingly used to characterize transcription factor binding and chromatin marks at a genomic scale. Various tools are now available to extract binding motifs from peak data sets. However, most approaches are only available as command-line programs, or via a website but with size restrictions. We present peak-motifs, a computational pipeline that discovers motifs in peak sequences, compares them with databases, exports putative binding sites for visualization in the UCSC genome browser and generates an extensive report suited for both naive and expert users. It relies on time- and memory-efficient algorithms enabling the treatment of several thousand peaks within minutes. Regarding time efficiency, peak-motifs outperforms all comparable tools by several orders of magnitude. We demonstrate its accuracy by analyzing data sets ranging from 4000 to 1 28 000 peaks for 12 embryonic stem cell-specific transcription factors. In all cases, the program finds the expected motifs and returns additional motifs potentially bound by cofactors. We further apply peak-motifs to discover tissue-specific motifs in peak collections for the p300 transcriptional co-activator. To our knowledge, peak-motifs is the only tool that performs a complete motif analysis and offers a user-friendly web interface without any restriction on sequence size or number of peaks.

Mathematical modelling of cell-fate decision in response to death receptor engagement.

Cytokines such as TNF and FASL can trigger death or survival depending on cell lines and cellular conditions. The mechanistic details of how a cell chooses among these cell fates are still unclear. The understanding of these processes is important since they are altered in many diseases, including cancer and AIDS. Using a discrete modelling formalism, we present a mathematical model of cell fate decision recapitulating and integrating the most consistent facts extracted from the literature. This model provides a generic high-level view of the interplays between NFκB pro-survival pathway, RIP1-dependent necrosis, and the apoptosis pathway in response to death receptor-mediated signals. Wild type simulations demonstrate robust segregation of cellular responses to receptor engagement. Model simulations recapitulate documented phenotypes of protein knockdowns and enable the prediction of the effects of novel knockdowns. In silico experiments simulate the outcomes following ligand removal at different stages, and suggest experimental approaches to further validate and specialise the model for particular cell types. We also propose a reduced conceptual model implementing the logic of the decision process. This analysis gives specific predictions regarding cross-talks between the three pathways, as well as the transient role of RIP1 protein in necrosis, and confirms the phenotypes of novel perturbations. Our wild type and mutant simulations provide novel insights to restore apoptosis in defective cells. The model analysis expands our understanding of how cell fate decision is made. Moreover, our current model can be used to assess contradictory or controversial data from the literature. Ultimately, it constitutes a valuable reasoning tool to delineate novel experiments.

Logical modelling of regulatory networks with GINsim 2.3

Many important problems in cell biology require the consideration of dense nonlinear interactions between functional modules. The requirement of computer simulation for the understanding of cellular processes is now widely accepted, and a variety of modelling frameworks have been designed to meet this need. Here, we present a novel public release of the Gene Interaction Network simulation suite (GINsim), a software designed for the qualitative modelling and analysis of regulatory networks. The main functionalities of GINsim are illustrated through the analysis of a logical model for the core network controlling the fission yeast cell cycle. The last public release of GINsim (version 2.3), as well as development versions, can be downloaded from the dedicated website (http://gin.univ-mrs.fr/GINsim/), which further includes a model library, along with detailed tutorial and user manual.

Diversity and Plasticity of Th Cell Types Predicted from Regulatory Network Modelling

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.

Segmenting the fly embryo: a logical analysis of the pair-rule cross-regulatory module

This manuscript reports a dynamical analysis of the pair-rule cross-regulatory module controlling segmentation in Drosophila melanogaster. We propose a logical model accounting for the ability of the pair-rule module to determine the formation of alternate juxtaposed Engrailed- and Wingless-expressing cells that form the (para)segmental boundaries. This module has the intrinsic capacity to generate four distinct expression states, each characterized by the expression of a particular combination of pair-rule genes or expression mode. The selection of one of these expression modes depends on the maternal and gap inputs, but also crucially on cross-regulations among pair-rule genes. The latter are instrumental in the interpretation of the maternal-gap pre-pattern. Our logical model allows the qualitative reproduction of the patterns of pair-rule gene expressions corresponding to the wild type situation, to loss-of-function and cis-regulatory mutations, and to ectopic pair-rule expressions. Furthermore, this model provides a formal explanation for the morphogenetic role of the initial bell-shaped expression of the gene even-skipped, i.e. for the distinct effects of different levels of the Even-skipped protein on its target pair-rule genes. It also accounts for the requirement of Even-skipped for the formation of all Engrailed-stripes. Finally, it provides new insights into the roles and evolutionary origins of the apparent redundancies in the regulatory architecture of the pair-rule module.

Functional organisation of Escherichia coli transcriptional regulatory network

Taking advantage of available functional data associated with 115 transcription and 7 sigma factors, we have performed a structural analysis of the regulatory network of Escherichia coli. While the mode of regulatory interaction between transcription factors (TFs) is predominantly positive, TFs are frequently negatively autoregulated. Furthermore, feedback loops, regulatory motifs and regulatory pathways are unevenly distributed in this network. Short pathways, multiple feed-forward loops and negative autoregulatory interactions are particularly predominant in the subnetwork controlling metabolic functions such as the use of alternative carbon sources. In contrast, long hierarchical cascades and positive autoregulatory loops are overrepresented in the subnetworks controlling developmental processes for biofilm and chemotaxis. We propose that these long transcriptional cascades coupled with regulatory switches (positive loops) for external sensing enable the coexistence of multiple bacterial phenotypes. In contrast, short regulatory pathways and negative autoregulatory loops enable an efficient homeostatic control of crucial metabolites despite external variations. TFs at the core of the network coordinate the most basic endogenous processes by passing information onto multi-element circuits. Transcriptional expression data support broader and higher transcription of global TFs compared to specific ones. Global regulators are also more broadly conserved than specific regulators in bacteria, pointing to varying functional constraints.