A comprehensive modular map of molecular interactions in RB/E2F pathway
2008; Springer Nature; Volume: 4; Issue: 1 Linguagem: Inglês
10.1038/msb.2008.7
ISSN1744-4292
AutoresLaurence Calzone, Amélie Gelay, Andreï Zinovyev, François Radvanyi, Emmanuel Barillot,
Tópico(s)Genomics and Chromatin Dynamics
ResumoReview4 March 2008Open Access A comprehensive modular map of molecular interactions in RB/E2F pathway Laurence Calzone Laurence Calzone Institut Curie, Service Bioinformatique, Paris INSERM, U900, Paris, F-75248 France Ecole des Mines de Paris, ParisTech, Fontainebleau, F-77300 France Search for more papers by this author Amélie Gelay Amélie Gelay Institut Curie, Service Bioinformatique, Paris Institut Curie, CNRS UMR144, Paris, France Search for more papers by this author Andrei Zinovyev Corresponding Author Andrei Zinovyev Institut Curie, Service Bioinformatique, Paris INSERM, U900, Paris, F-75248 France Ecole des Mines de Paris, ParisTech, Fontainebleau, F-77300 France Search for more papers by this author François Radvanyi François Radvanyi Institut Curie, Centre de Recherche, Paris, F-75248 France Institut Curie, CNRS UMR144, Paris, FranceThese two authors are the joint senior authors of this paper. Search for more papers by this author Emmanuel Barillot Emmanuel Barillot Institut Curie, Service Bioinformatique, Paris INSERM, U900, Paris, F-75248 France Ecole des Mines de Paris, ParisTech, Fontainebleau, F-77300 FranceThese two authors are the joint senior authors of this paper. Search for more papers by this author Laurence Calzone Laurence Calzone Institut Curie, Service Bioinformatique, Paris INSERM, U900, Paris, F-75248 France Ecole des Mines de Paris, ParisTech, Fontainebleau, F-77300 France Search for more papers by this author Amélie Gelay Amélie Gelay Institut Curie, Service Bioinformatique, Paris Institut Curie, CNRS UMR144, Paris, France Search for more papers by this author Andrei Zinovyev Corresponding Author Andrei Zinovyev Institut Curie, Service Bioinformatique, Paris INSERM, U900, Paris, F-75248 France Ecole des Mines de Paris, ParisTech, Fontainebleau, F-77300 France Search for more papers by this author François Radvanyi François Radvanyi Institut Curie, Centre de Recherche, Paris, F-75248 France Institut Curie, CNRS UMR144, Paris, FranceThese two authors are the joint senior authors of this paper. Search for more papers by this author Emmanuel Barillot Emmanuel Barillot Institut Curie, Service Bioinformatique, Paris INSERM, U900, Paris, F-75248 France Ecole des Mines de Paris, ParisTech, Fontainebleau, F-77300 FranceThese two authors are the joint senior authors of this paper. Search for more papers by this author Author Information Laurence Calzone1,2,3,‡, Amélie Gelay1,5,‡, Andrei Zinovyev 1,2,3,‡, François Radvanyi4,5 and Emmanuel Barillot1,2,3 1Institut Curie, Service Bioinformatique, Paris 2INSERM, U900, Paris, F-75248 France 3Ecole des Mines de Paris, ParisTech, Fontainebleau, F-77300 France 4Institut Curie, Centre de Recherche, Paris, F-75248 France 5Institut Curie, CNRS UMR144, Paris, France ‡These authors contributed equally to this work. *Corresponding author. Institut Curie, Service Bioinformatique, 26 rue d'Ulm, Paris 75248, France. Tel. +33 1 53 10 70 50; Fax: +33 1 44 41 09 08; E-mail: [email protected] Molecular Systems Biology (2008)4:0174https://doi.org/10.1038/msb.2008.7 PDFDownload PDF of article text and main figures. ToolsAdd to favoritesDownload CitationsTrack CitationsPermissions Figures & Info We present, here, a detailed and curated map of molecular interactions taking place in the regulation of the cell cycle by the retinoblastoma protein (RB/RB1). Deregulations and/or mutations in this pathway are observed in most human cancers. The map was created using Systems Biology Graphical Notation language with the help of CellDesigner 3.5 software and converted into BioPAX 2.0 pathway description format. In the current state the map contains 78 proteins, 176 genes, 99 protein complexes, 208 distinct chemical species and 165 chemical reactions. Overall, the map recapitulates biological facts from approximately 350 publications annotated in the diagram. The network contains more details about RB/E2F interaction network than existing large-scale pathway databases. Structural analysis of the interaction network revealed a modular organization of the network, which was used to elaborate a more summarized, higher-level representation of RB/E2F network. The simplification of complex networks opens the road for creating realistic computational models of this regulatory pathway. Introduction The cell cycle is the succession of four phases called G1, S, G2 and M. In dividing cells, DNA replication (S phase) and mitosis (M phase) alternate (Alberts et al, 1994), and are separated by two gap phases, G1 and G2 phases. In quiescent cells, the cells are considered to be in G0 phase. When they receive external signals, such as growth factors, a series of activations push the cell from a G0 to a G1 state and enters the cell cycle. The whole process of cell division is mainly orchestrated by complexes composed of two subunits, a kinase and a cyclin partner. These complexes phosphorylate a certain number of proteins, either activating or inhibiting them. Among them, the retinoblastoma tumour suppressor protein RB (RB1) is a key regulator in cell-cycle entry (transition G1/S). It sequesters a family of transcription factors, the E2Fs, responsible for the transcription of many genes involved in cell-cycle regulation, DNA replication and other functions like the activation of the apoptotic pathway (Muller et al, 2001). RB functions as a brake in the cell cycle, which is released when external signals trigger S-phase entry. The main targets of the external signals are the G1 cyclin/CDK complexes. Once active, the complexes, among them CycD1/CDK4,6, act as starters of the cell cycle (Novak et al, 2007) and phosphorylate RB, which then releases E2F (DeGregori, 2004). RB is a member of a family of proteins called the pocket proteins (Knudsen and Wang, 1997). These proteins RB, p107 and p130, share sequence similarities, especially in the 'pocket domain' (Stevaux and Dyson, 2002), which is responsible for their repressor function. RB protein contains domains where the binding sites for co-repressors (E2F proteins and viral oncoproteins) are situated. These sites are subjected to most mutations. RB is a tumour suppressor gene. Because of its implication in so many, if not all, cancers (Sherr and McCormick, 2002), the study of RB regulation requires a special attention. More specifically, the RB/E2F pathway is commonly deregulated in cancer through genetic or epigenetic mechanisms, resulting in E2F activation. Several common oncogenes (involved in many cancer types) are the activators of the pathway, whereas several common tumour suppressor genes are inhibitors of the pathway. For example, cyclin D1 (CCND1), E2F3 and the two cyclin-dependent kinases CDK4 and CDK6 can be activated by translocation, amplification or mutation, whereas RB (RB1) and the cyclin-dependent kinase inhibitors p16INK4a (CNKN2A) and p15INK4b (CDKN2B) can be inactivated by point mutation, homozygous deletion or DNA methylation. In addition, RB can be inactivated by several oncogenic viral proteins including E7 from human papillomavirus, which is responsible for more than 90% of cervical carcinomas (Munger et al, 2001). Tumour suppressor gene inactivation is found not only in sporadic tumours but also in tumour-prone families. Germline mutations of RB1 results in retinoblastoma with a high penetrance early in young individuals and late in life in sarcomas and lung and bladder carcinomas (Knudson, 1971; Nevins, 2001; Giacinti and Giordano, 2006). Germinal mutations of p16INK4a results in melanomas (Hussussian et al, 1994). Finally, the importance of the RB/E2F pathway in tumour formation based on genetic or epigenetic arguments, has been confirmed by in vivo experiments using transgenic mice (Classon and Harlow, 2002). Often, textbooks present a simplified picture of RB regulation (Figure 1): mitogens activate CycD1/CDK complexes, which phosphorylate RB, and which, in turn, releases the hold on E2F transcription factors. E2F activity controls cell-cycle progression. Figure 1.Textbook view of the RB/E2F pathway. When CycD1/CDK4,6 is activated, it phosphorylates RB, which in turn is inactivated and releases E2F transcription factors leading to cell proliferation. Download figure Download PowerPoint However, the real picture of interactions around RB is much more complex due to combinatorial complexity and the number of important regulators, which are not taken into account in the simple pathway view. There are not only several E2F transcription factors but also several proteins besides RB with which the E2F proteins bind to be inactivated. Moreover, the RB-related proteins not only can sequester the E2Fs but also can actively repress transcription. There are 10 known E2F proteins—7 of which bind to dimerization partners (DP1 and DP2) to become active—that can associate with the three pocket proteins RB, p107 and p130, and proteins from the polycomb group. These E2F proteins can be (1) activators of transcription (E2F1, E2F2 and E2F3a) or (2) repressors of transcription (E2F3b, E2F4, E2F5 and E2F6). The last E2F family proteins, E2F7a, E2F7b and E2F8, form homodimers, do not bind to pocket proteins and seem to repress transcription. We believe that this complexity has to be taken into account in any realistic consideration of the RB network (such as the creation of a computational model). Thus, we stated our problem as a very careful and focused consideration of the whole corpus of available biological facts published in high-level biological literature to reconstruct a detailed network of interactions around RB protein. We started from the simplified picture (Figure 1) and added several tens of proteins involved in cell-cycle regulation, carefully connecting them to the rest of the network by means of known biochemical transformations. There were several questions that had to be answered during this reconstruction: What should the language be for such a detailed description? It should be standard and expressive enough to readily incorporate available biological facts. A choice was made to use Systems Biology Graphical Notation (SBGN) standard partially implemented in CellDesigner software (Kitano et al, 2005). This choice also defined for us the level of details that we should keep in the description. How to deliver the collected information to the public? Together with CellDesigner diagram, we chose to use the BioPAX format for storing pathway information. An automatic conversion tool from CellDesigner (Funahashi et al, 2003) to BioPAX was developed (Zinovyev et al, 2007). How to define the 'borders' of RB pathway? RB interacts with more than 100 proteins (Morris and Dyson, 2001). We chose to concentrate mainly on genes involved in cell-cycle regulation and in apoptosis entry. How to make the resulting detailed network usable and understandable? We performed a structural analysis of the resulting network, and it was found that the network has a modular structure. A module was defined as a cluster of relevant cycles in the reaction graph. A special software, BiNoM, plugin of Cytoscape (Zinovyev et al, 2007), was developed to perform such analysis (see Materials and methods). How to maintain and expand the pathway? A web-page is associated with the pathway (http://bioinfo-out.curie.fr/projects/rbpathway/) and will be regularly updated as some new papers appear and with the input we expect from the RB experts worldwide (listed in 'Contributors' section on our web-page). One purpose behind the construction of this diagram, once validated, is to provide a map of RB/E2F pathway that can become not only a reference while studying different cancers and mutations but also a tool to analyse formally the pathway and predict its behaviour in response to different types of deregulations (Nevins, 2001). Table I recapitulates all current names for genes used in this article and their corresponding HUGO names. RB, p107 and p130 are the common names for the proteins referenced as RB1, RBL1 and RBL2 (HUGO names), respectively. Table 1. List of protein names used in the diagram and the corresponding HUGO names Cell designer name Used name HUGO name AKT1 AKT1 AKT1 APC Anaphase-promoting complex ATM ATM ATM ATR ATR ATR BAT8 G9a BAT8 BMI1* BMI1 PCGF4 CDC14A CDC14A CDC14A CDC2 CDK1 CDC2 CDC20 CDC20 CDC20 CDC25C CDC25C CDC25C CDC37 CDC37 CDC37 CDH1* CDH1 FZR1 CDK2 CDK2 CDK2 CDK3 CDK3 CDK3 CDK4 CDK4 CDK4 CDK6 CDK6 CDK6 CDK7 CDK7 CDK7 CDKN3 KAP1 CDKN3 CHEK1 Chk1 CHEK1 CHEK2 Chk2 CHEK2 CREBBP CBP CREBBP cyclin A2* cyclin A2, CycA2 CCNA2 cyclin B1* cyclin B1, CycB1 CCNB1 cyclin C* cyclin C, CycC CCNC cyclin D1* cyclin D1, CycD1 CCND1 cyclin E1* cyclin E1, CycE1 CCNE1 cyclin H* cyclin H, CycH CCNH DP1* DP1 TFDP1 DP2* DP2 TFDP2 E2F1 E2F1 E2F1 E2F2 E2F2 E2F2 E2F3 E2F3 E2F3 E2F4 E2F4 E2F4 E2F5 E2F5 E2F5 E2F6 E2F6 E2F6 E2F7 E2F7 E2F7 E2F8 E2F8 E2F8 EED EED EED EHMT1 Eu-HMT1 EHMT1 EP300 p300 EP300 EPC1 EPC1 EPC1 EZH2 EZH2 EZH2 GSK3B GSK3B GSK3B HDAC1 HDAC1 HDAC1 HP1gamma* HP1gamma CBX3 HSP90 HSP90 family MAX MAX MAX MDM2 MDM2 MDM2 MEL-18* MEL-18 PCGF2 MGA MGA MGA MNAT1 MAT1 MNAT1 NBS1 NBS1 NBS1 p107* p107 RBL1 p130* p130 RBL2 p14ARF* p14ARF CDKN2A p15INK4b* p15INK4b CDKN2B p16INK4a* p16INK4a CDKN2A p18INK4c* p18INK4c CDKN2C p19INK4d* p19INK4d CDKN2D p21Cip1* p21Cip1 CDKN1A p27Kip1* p27Kip1 CDKN1B p53* p53 TP53 p57Kip2* p57Kip2 CDKN1C PCAF P/CAF PCAF PCNA PCNA PCNA PHC1 Mph1 PHC1 PKMYT1 Myt1 PKMYT1 PP1 PP1 PP1 pRB* RB, pRB RB1 RING1 RING1 RING1 RYBP RYBP RYBP SERTAD1 SERTAD1 SERTAD1 SIN3B SIN3B SIN3B SUV39H1 SUV39H1 SUV39H1 SWI/SNF SWI/SNF SWI/SNF TFIIH* TFIIH TOPBP1 TOPBP1 TOPBP1 WEE1 WEE1 WEE1 Asterisks are added to the CellDesigner names each time the name used does not correspond to the HUGO name. A comprehensive map of RB/E2F pathway Map of the pathway in CellDesigner and BioPAX representations Figure 2 shows the resulting comprehensive map of RB pathway. The diagram utilizes SBGN system (http://www.sbgn.org) to represent proteins and their specific modifications, protein complexes and genes, as well as various protein transformations (binding, unbinding, phosphorylation, acetylation, transport and so on) and their effects of activation or inhibition of chemical reactions, including transcriptional activation or inhibition. Figure 2.The textbook pathway of RB has been expanded by integrating data from the literature. The E2F transcription factors (represented here by single proteins in the nuclear compartment) are connected by activation and inhibition arrows to their gene targets. (A) Map of target genes of E2F transcription factors. Each E2F associates with different cofactors to activate or inhibit the transcription of many genes; pointed arrows mean activation and flat arrows mean inhibitions (B) Map of protein–protein interaction network. Each icon on the diagram represents distinct chemical species. See Kitano and co-workers' description of CellDesigner's standard notation (Kitano et al, 2005) for a detailed meaning of shapes. When the information is available (from Atlas Oncology web-page: http://www.atlasgeneticsoncology.org/), tumour suppressor genes and the corresponding proteins are coloured in blue and oncogenes in red, the other proteins are in green. To read and navigate through the map, visit our webpage: http://bioinfo-out.curie.fr/projects/rbpathway/. The map is clickable and allows easy access to all included information (such as literature references or standard protein ids) and hyperlinked to other databases. Download figure Download PowerPoint In the general organization of the diagram, there is a part representing the set of genes regulated by the E2F family of transcription factors (2A) and another one representing the location of the different proteins in the three different cellular compartments (the nucleus, nucleolus and cytoplasm) (2B). The resulting map has a total of 78 proteins, which are represented in 208 distinct chemical species (158 of them are located in the nucleus, 47 in the cytoplasm and 3 in the nucleolus), 526 reactions and regulations (among them, 57 protein associations, 13 protein dissociations, 68 post-translational modifications, 361 transcriptional regulations and 27 transport pseudo-reactions), 176 genes, and compiles experimental results from more than 350 publications (301 research papers and 56 reviews). Once compiled in CellDesigner software, the pathway information was translated into BioPAX format using BiNoM Cytoscape plugin (see Materials and methods). All the files in SBML and BioPAX formats are available in the Supplementary information. An interactive version of the map can also be found on the web site. We compared our RB/E2F map to existing pathway databases such as Reactome (Joshi-Tope et al, 2005) or Transpath (Krull et al, 2003) and concluded that we provide a more systematic description of RB/E2F network than that contained in the general purpose pathway databases (details at http://bioinfo-out.curie.fr/projects/rbpathway/Comparison.html). Modular decomposition of the pathway Using clustering of relevant cycles in the reaction graph (see Materials and methods for details) and consequent manual curation, the detailed RB pathway was decomposed into 16 network modules. As a result, almost every module can be thought as a detailed sequence of events that occur with a particular protein or protein complex whose name designates the whole module. Thus, we identified RB, E2F1-3, E2F4-5, E2F6-8, CycC/CDK3, CycH/CDK7, CycD1/CDK4,6, CycE1/CDK2, CycA2/CDK2, CycB1/CDC2, p16INK4a/p15INK4b, p27KIP1/p21CIP1, CDC25C and Wee1, APC and Apoptosis entry modules. Eight additional modules were considered representing the transcriptional targets of each of the E2Fs described in the RB/E2F map. Identification of these modules allowed us to compile a modular RB pathway view (Figure 3), where network modules are connected by 'activation' and 'inhibition' relations. The information about these relations is derived from the detailed diagram. For example, in the detailed map, E2F1 is phosphorylated by CycA2/CDK2 and is subsequently recognized for degradation, which is translated in the modular map by CycA2/CDK2 module inhibiting E2F1-3 module. Figure 3.Modular view of the comprehensive map presented in Figure 2. The comprehensive map has been simplified using curated structural analysis techniques to divide the graph in modules. The links between the nodes represent the influence that a module has on the others. There are 16 protein modules (green rounded rectangles) and 8 E2F target gene modules (yellow rectangles). Download figure Download PowerPoint Because of the space limitation, the module descriptions and figures are provided below only for E2Fs and RB modules (Figure 4). However, the 'clickable' modular view leading to a description and the corresponding figures of each of the modules can be found in the Supplementary information (http://bioinfo-out.curie.fr/projects/rbpathway/Modules.html). Figure 4.Cytoscape view of (A) E2F1-3 module and (B) RB module. Download figure Download PowerPoint The RB module Two cycles appear in this module revealing two roles of RB. (1) A transcriptional repressor: HDAC1 is a marker of transcriptional repression and seems to require the chromatin remodeler SWI-SNF to form a repressor complex with DP1/E2F1/RB (Frolov and Dyson, 2004). The complex recruits both HP1-γ and SUV39H1 and continues to prevent transcription (Nielsen et al, 2001). (2) A repressor of E2F1/DP1: RB binds to E2F1/DP1 and blocks its transcriptional activity. The affinity between E2F1/DP1 and RB is decreased through sequential phosphorylations by the Cyclin/CDK complexes. The first phosphorylations by CycC/CDK3 favours the passage from G0 to G1 phase, then CycD1/CDK4,6 modifies its conformation and releases HDAC1 (Zhang et al, 2000), revealing a new site of phosphorylation targeted by CycE1/CDK2 (Vidal and Koff, 2000; Muchardt and Yaniv, 2001). The phosphorylation by CycD1/CDK4,6 already allows some genes to be transcribed (such as Cyclin E). The complexes CycD1/CDK4 and CycD1/CDK6 act as sensors of growth factors. When activated, they precipitate the cells into S phase by further phosphorylating RB (Weinberg, 1995; Planas-Silva and Weinberg, 1997). When RB is hyperphosphorylated, the complex dissociates and E2F1/DP1 is released from the inactive complex. Later, RB is dephosphorylated by the phosphatase PP1 towards the end of M phase and able to repress E2Fs again (Vidal and Koff, 2000). The E2F transcription factors modules The pocket proteins RB, p130 and p107 inhibit a family of transcription factors, the E2F, through association and can also act as active repressors by recruiting other partners. We already mentioned three subgroups of E2F transcription factors. The first group contains activators of transcription: E2F1, E2F2 and E2F3a, which bind to RB; the second one regroups inhibitors of transcription: E2F3b, E2F4 and E2F5, which bind to either p107 or p130 and to some extent RB; and the third group includes E2F6, E2F7 and E2F8, which do not need to bind to pocket proteins to become active repressors of transcription. These three groups correspond to three different modules that are described below. Activator E2F1-3 module The E2F1-3 module includes the E2F activator family of transcription E2F1, E2F2 and E2F3a. Even though the proteins slightly differ in their cell-cycle role, we chose to describe the three E2F transcription factors as one, as they share a lot of functional similarities. In a future version of the pathway, we plan to differentiate the activity of the three transcription activators. Experiments show that dimerization between E2F1 and its partner DP1 is stable and that E2F1 stimulates nuclear localization of DP1 (Magae et al, 1996). E2F1/DP1 is acetylated by three acetyltransferases—PCAF, CREB binding protein and p300—to stabilize the E2F1 protein (Frolov and Dyson, 2004). The acetylated complex is capable of binding to PCAF to form an active dimer. The complex ability of binding to DNA on the promoter sites of its target genes along with its transcriptional activity is increased during the G1'S transition. At G2, the complex is phosphorylated by CycA2/CDK2 (He and Cress, 2002). The affinity between E2F1 and DP1 is then diminished leading to the dissociation of the complex (Tsantoulis and Gorgoulis, 2005) and the release of PCAF. The proteins undergo further modifications before degradation: E2F1 is deacetylated by HDAC1 (Martinez-Balbas et al, 2000), dephosphorylated and phosphorylated de novo during S phase by TFIIH kinase for rapid degradation (Ianari et al, 2004). On DNA damage, the complex PCAF/E2F1/DP1 can be phosphorylated and stabilized either by CHEK1 and CHEK2 through phosphorylation at Ser-364 or by ATM and ATR (Dimova and Dyson, 2005; Powers et al, 2004), preventing E2F1 ubiquitination (Wang et al, 2004). E2F1 mediates the transcription of many genes involved in apoptosis. However, E2F1 transcriptional activity can also be inhibited when bound to the topoisomerase TopBP1 to give time to the cell to repair the damage (Liu et al, 2003). Repressor E2F4–5 module E2F4 associates successively with two different pocket proteins: p130 in quiescent cells and p107 in proliferating cells. E2F4 is initially found in the cytoplasm. When in complex, E2F4 is translocated in the nucleus, where it acts as a repressor of transcription during G0 and G1 phases (Verona et al, 1997). Once in the nucleus, the complex binds to DNA and recruits some co-repressors of transcription: the chromatin remodeler Sin3B, the deacetylase histone HDAC1 and the methyltransferase histones SUV39H1 (Rayman et al, 2002; Liu et al, 2005). At the G0–G1 transition, when the quiescent cells receive signals from growth factors, p130 starts to be phosphorylated and gets dissociated from the complex it was forming with E2F4/DP2. Later, p130 can also be phosphorylated by CycD1/CDK4,6 and CycE1/CDK2 when present and degraded in late G1. It is then replaced by p107 (Farkas et al, 2002). When E2F4/DP2 is in complex with p107, it continues to repress transcription of target genes until it is phosphorylated by CycD1/CDK4,6. E2F4 is then translocated to the cytoplasm, where it can no longer repress transcription, whereas both p130 and p107 are able to inhibit CycE1/CDK2 and CycA2/CDK2 activities (Litovchick et al, 2004). In this module, E2F3b and E2F5 should also be considered. Their repressive role in the pathway is not described yet but will be detailed in future versions of the pathway. Repressor E2F6–8 module E2F6 seems to play a role in S-phase entry. It binds to both DP partners, DP1 and DP2, in the cytoplasm, and the complexes are then translocated in the nucleus. As opposed to other E2F family members, E2F6 does not associate with pocket proteins but rather recruits some proteins from the polycomb group (PcG) to repress transcription. The first polycomb group with which E2F6 is involved includes BMI1, Mel-18, PHC1, RING1 and RYBP and acts as a repressor of transcription (Trimarchi et al, 2001; Sánchez-Beato et al, 2004). Similarly, E2F6, Max and HP1-γ bind to form a complex that has also revealed a repressive transcription role in quiescent cells (Ogawa et al, 2002). However, E2F6 seems to intervene in other transitions of the cell cycle than only that of G0–G1. Another complex involving EPC1 and EZH2 has been found in proliferating cells: E2F6/DP1 is repressing transcription when in complex with EPC1 alone or with both EPC1 and EZH2 (Attwooll et al, 2005). E2F6/DP1 binds consecutively to EPC1 and then to both EZH2 and EED. The roles of the different complexes are not clearly established yet and more experiments will be needed to confidently describe their specific actions in proliferating and quiescent cells. E2F7 and E2F8 regulations have not been carefully defined in our pathway yet. Both are known to inhibit transcription although (de Bruin et al, 2003; Logan et al, 2005). The details will be added as more publications on their role in RB pathway appear. RB pathway transcriptional activity modules The eight 'E2F*_targets' modules correspond to the gene targets of the eight E2F transcription factors. The three transcription factors E2F1, E2F2 and E2F3, although not having well-documented differences at the level of protein–protein interactions and represented by a generic entity in our diagram (Figure 2B), show their specificities at the transcription level. For this reason, the generic entity is decomposed into the three individual components in the upper part of the diagram (Figure 2A). According to our diagram, RB/E2F pathway is a self-regulating molecular mechanism, as there exist multiple positive and negative feedbacks through transcription: among the 78 proteins described in the diagram, 23 are targets of the E2Fs transcription factors, both activators and inhibitors. The structure of these feedbacks is detailed in Supplementary Figure S1 and Supplementary Table S1. E2F1 is a target gene for all E2Fs with the exception of E2F5. In turn, E2F1 regulates expression of E2F1, E2F2, E2F6, E2F7 and E2F8 genes resulting in several negative and positive feedback control circuits. To analyse the differences between the E2Fs and their roles in other contexts, we calculated the significance of the overlaps of all E2F1-8 transcriptional targets with other known pathways such as MSigDB database (BROAD, MIT) (Subramanian et al, 2005). The results are provided in Supplementary Table S2. As expected, these lists have highly significant overlaps with cell-cycle-, G1'S transition-, Rb/E2F- or p21-related pathways. Some other overlaps highlight the involvement of RB/E2F in some differentiation-related pathways and reciprocally. However, only E2F1 target list is linked to various apoptosis-related pathways (APOPTOSIS_GENMAP, APOPTOSIS_KEGG, DEATH_PATHWAY and others). In particular, E2F1, but not the other E2Fs, targets p53 tumour suppressor gene in RB/E2F pathway. This confirms the recent findings that E2F1 is the only specific inducer of apoptosis among the E2F transcription factors, even though its level can be regulated by other E2Fs (Lazzerini Denchi and Helin, 2005). Case study of bladder tumour data As an example of the potentialities of the map we have assembled, we performed a study on 55 bladder tumours (Stransky et al, 2006). This study exemplified how the map and its modular decomposition can be used to explain the differences between two different tumour progression pathways and/or different stages of cancer progression. More specifically, we analysed transcriptome and comparative genomics hybridization data collected for 55 patients with bladder cancers. We verified which groups of genes (modules) are up- or downregulated in invasive cancers and identified two different paths that both lead to invasive aggressive cancers. This study confirmed known facts about bladder cancer, for example overexpression of CCND1 in low-grade tumours, and led to new observations, for example downregulation of E2F4-5, E2F6-8, Wee1, APC modules in invasive cancers. The approach and our map can similarly be used for comparing other biological contexts. Indeed, further analyses will be developed to obtain more insights into molecular mechanisms of cancer progression. The details and results of this preliminary study can be found on our webpage: http://bioinfo-out.curie.fr/projects/rbpathway/case_study.html. Discussion In this paper, we present a comp
Referência(s)