RESUMO
This review summarizes current advances in the role of transcriptional stochasticity in HIV-1 latency, which were possible in a large part due to the development of single-cell approaches. HIV-1 transcription proceeds in bursts of RNA production, which stem from the stochastic switching of the viral promoter between ON and OFF states. This switching is caused by random binding dynamics of transcription factors and nucleosomes to the viral promoter and occurs at several time scales from minutes to hours. Transcriptional bursts are mainly controlled by the core transcription factors TBP, SP1 and NF-κb, the chromatin status of the viral promoter and RNA polymerase II pausing. In particular, spontaneous variability in the promoter chromatin creates heterogeneity in the response to activators such as TNF-α, which is then amplified by the Tat feedback loop to generate high and low viral transcriptional states. This phenomenon is likely at the basis of the partial and stochastic response of latent T cells from HIV-1 patients to latency-reversing agents, which is a barrier for the development of shock-and-kill strategies of viral eradication. A detailed understanding of the transcriptional stochasticity of HIV-1 and the possibility to precisely model this phenomenon will be important assets to develop more effective therapeutic strategies.
Assuntos
Soropositividade para HIV , HIV-1 , Humanos , HIV-1/genética , Cromatina , NF-kappa B , NucleossomosRESUMO
Monitoring transcription in living cells gives access to the dynamics of this complex fundamental process. It reveals that transcription is discontinuous, whereby active periods (bursts) are separated by one or several types of inactive periods of distinct lifetimes. However, decoding temporal fluctuations arising from live imaging and inferring the distinct transcriptional steps eliciting them is a challenge. We present BurstDECONV, a novel statistical inference method that deconvolves signal traces into individual transcription initiation events. We use the distribution of waiting times between successive polymerase initiation events to identify mechanistic features of transcription such as the number of rate-limiting steps and their kinetics. Comparison of our method to alternative methods emphasizes its advantages in terms of precision and flexibility. Unique features such as the direct determination of the number of promoter states and the simultaneous analysis of several potential transcription models make BurstDECONV an ideal analytic framework for live cell transcription imaging experiments. Using simulated realistic data, we found that our method is robust with regards to noise or suboptimal experimental designs. To show its generality, we applied it to different biological contexts such as Drosophila embryos or human cells.
Assuntos
Drosophila , Transcrição Gênica , Animais , Humanos , Drosophila/genética , Regiões Promotoras GenéticasRESUMO
The MAPK/ERK pathway is an essential intracellular signaling pathway. Its deregulation is involved in tumor transformation and progression. The discovery of activating mutations of BRAF in various cancers has opened new therapeutic avenues with BRAF protein kinase inhibitors. Depending on the type of cancers, these inhibitors have shown either insufficient efficacy due to primary resistance of tumor cells or transient efficacy due to the development of acquired resistance. In this review, we revisit the discoveries that led to the development of BRAF inhibitors and detail the molecular and cellular mechanisms of resistance in cancers treated with these inhibitors. Understanding these mechanisms is crucial for developing more efficient therapeutic strategies.
Title: La résistance aux inhibiteurs de BRAF - Les leçons de la clinique. Abstract: La voie de signalisation MAPK/ERK est une voie centrale de la signalisation intracellulaire. Sa dérégulation participe à la transformation et la progression tumorales. Dans plusieurs cancers, la découverte de mutations activatrices de BRAF, à l'origine de l'activation de cette voie, a ouvert de nouvelles perspectives thérapeutiques avec le développement d'inhibiteurs spécifiques de la protéine. Selon les cancers, ces inhibiteurs ont cependant montré soit une efficacité insuffisante, due à la résistance primaire des cellules tumorales, soit une efficacité transitoire, due à l'apparition d'une résistance acquise. Dans cette revue, nous revenons sur les découvertes qui ont conduit au développement de ces inhibiteurs de BRAF. Nous détaillons également les mécanismes moléculaires et cellulaires de la résistance à ces inhibiteurs observée dans différents types de cancers. Comprendre ces mécanismes est en effet primordial pour développer des stratégies thérapeutiques qui soient plus efficaces.
Assuntos
Neoplasias , Proteínas Proto-Oncogênicas B-raf , Linhagem Celular Tumoral , Resistencia a Medicamentos Antineoplásicos/genética , Humanos , Sistema de Sinalização das MAP Quinases , Mutação , Neoplasias/tratamento farmacológico , Neoplasias/genética , Neoplasias/metabolismo , Inibidores de Proteínas Quinases/farmacologia , Inibidores de Proteínas Quinases/uso terapêutico , Proteínas Proto-Oncogênicas B-raf/genética , Proteínas Proto-Oncogênicas B-raf/metabolismoRESUMO
Cutaneous melanoma is a highly invasive tumor and, despite the development of recent therapies, most patients with advanced metastatic melanoma have a poor clinical outcome. The most frequent mutations in melanoma affect the BRAF oncogene, a protein kinase of the MAPK signaling pathway. Therapies targeting both BRAF and MEK are effective for only 50% of patients and, almost systematically, generate drug resistance. Genetic and non-genetic mechanisms associated with the strong heterogeneity and plasticity of melanoma cells have been suggested to favor drug resistance but are still poorly understood. Recently, we have introduced a novel mathematical formalism allowing the representation of the relation between tumor heterogeneity and drug resistance and proposed several models for the development of resistance of melanoma treated with BRAF/MEK inhibitors. In this paper, we further investigate this relationship by using a new computational model that copes with multiple cell states identified by single cell mRNA sequencing data in melanoma treated with BRAF/MEK inhibitors. We use this model to predict the outcome of different therapeutic strategies. The reference therapy, referred to as "continuous" consists in applying one or several drugs without disruption. In "combination therapy", several drugs are used sequentially. In "adaptive therapy" drug application is interrupted when the tumor size is below a lower threshold and resumed when the size goes over an upper threshold. We show that, counter-intuitively, the optimal protocol in combination therapy of BRAF/MEK inhibitors with a hypothetical drug targeting cell states that develop later during the tumor response to kinase inhibitors, is to treat first with this hypothetical drug. Also, even though there is little difference in the timing of emergence of the resistance between continuous and adaptive therapies, the spatial distribution of the different melanoma subpopulations is more zonated in the case of adaptive therapy.
RESUMO
To maintain cellular identities during development, gene expression profiles must be faithfully propagated through cell generations. The reestablishment of gene expression patterns upon mitotic exit is mediated, in part, by transcription factors (TF) mitotic bookmarking. However, the mechanisms and functions of TF mitotic bookmarking during early embryogenesis remain poorly understood. In this study, taking advantage of the naturally synchronized mitoses of Drosophila early embryos, we provide evidence that GAGA pioneer factor (GAF) acts as a stable mitotic bookmarker during zygotic genome activation. We show that, during mitosis, GAF remains associated to a large fraction of its interphase targets, including at cis-regulatory sequences of key developmental genes with both active and repressive chromatin signatures. GAF mitotic targets are globally accessible during mitosis and are bookmarked via histone acetylation (H4K8ac). By monitoring the kinetics of transcriptional activation in living embryos, we report that GAF binding establishes competence for rapid activation upon mitotic exit.
Assuntos
Cromatina , Histonas , Acetilação , Animais , Cromatina/genética , Drosophila/genética , Histonas/genética , Histonas/metabolismo , Mitose/genética , Fatores de Transcrição/genética , Fatores de Transcrição/metabolismoRESUMO
Summary: Recently, symbolic computation and computer algebra systems have been successfully applied in systems biology, especially in chemical reaction network theory. One advantage of symbolic computation is its potential for qualitative answers to biological questions. Qualitative methods analyze dynamical input systems as formal objects, in contrast to investigating only part of the state space, as is the case with numerical simulation. However, corresponding tools and libraries have a different set of requirements for their input data than their numerical counterparts. A common format used in mathematical modeling of biological processes is Systems Biology Markup Language (SBML). We illustrate that the use of SBML data in symbolic computation requires significant pre-processing, incorporating external biological and mathematical expertise. ODEbase provides suitable input data derived from established existing biomodels, covering in particular the BioModels database. Availability and implementation: ODEbase is available free of charge at https://odebase.org.
RESUMO
Promoter-proximal pausing of RNA polymerase II is a key process regulating gene expression. In latent HIV-1 cells, it prevents viral transcription and is essential for latency maintenance, while in acutely infected cells the viral factor Tat releases paused polymerase to induce viral expression. Pausing is fundamental for HIV-1, but how it contributes to bursting and stochastic viral reactivation is unclear. Here, we performed single molecule imaging of HIV-1 transcription. We developed a quantitative analysis method that manages multiple time scales from seconds to days and that rapidly fits many models of promoter dynamics. We found that RNA polymerases enter a long-lived pause at latent HIV-1 promoters (>20 minutes), thereby effectively limiting viral transcription. Surprisingly and in contrast to current models, pausing appears stochastic and not obligatory, with only a small fraction of the polymerases undergoing long-lived pausing in absence of Tat. One consequence of stochastic pausing is that HIV-1 transcription occurs in bursts in latent cells, thereby facilitating latency exit and providing a rationale for the stochasticity of viral rebounds.
Assuntos
Regulação Viral da Expressão Gênica , Infecções por HIV/genética , HIV-1/genética , Regiões Promotoras Genéticas/genética , Latência Viral/genética , Algoritmos , RNA Polimerases Dirigidas por DNA/metabolismo , Infecções por HIV/metabolismo , Infecções por HIV/virologia , HIV-1/fisiologia , Células HeLa , Humanos , Modelos Genéticos , Processos Estocásticos , Fatores de Tempo , Ativação Viral/genética , Produtos do Gene tat do Vírus da Imunodeficiência Humana/genéticaRESUMO
Genes are expressed in stochastic transcriptional bursts linked to alternating active and inactive promoter states. A major challenge in transcription is understanding how promoter composition dictates bursting, particularly in multicellular organisms. We investigate two key Drosophila developmental promoter motifs, the TATA box (TATA) and the Initiator (INR). Using live imaging in Drosophila embryos and new computational methods, we demonstrate that bursting occurs on multiple timescales ranging from seconds to minutes. TATA-containing promoters and INR-containing promoters exhibit distinct dynamics, with one or two separate rate-limiting steps respectively. A TATA box is associated with long active states, high rates of polymerase initiation, and short-lived, infrequent inactive states. In contrast, the INR motif leads to two inactive states, one of which relates to promoter-proximal polymerase pausing. Surprisingly, the model suggests pausing is not obligatory, but occurs stochastically for a subset of polymerases. Overall, our results provide a rationale for promoter switching during zygotic genome activation.
Assuntos
Drosophila melanogaster/genética , Embrião não Mamífero/metabolismo , Regiões Promotoras Genéticas/genética , TATA Box/genética , Imagem com Lapso de Tempo/métodos , Transcrição Gênica/genética , Algoritmos , Animais , Animais Geneticamente Modificados , Drosophila melanogaster/embriologia , Drosophila melanogaster/metabolismo , Embrião não Mamífero/embriologia , Proteínas de Fluorescência Verde/genética , Proteínas de Fluorescência Verde/metabolismo , Cinética , Proteínas Luminescentes/genética , Proteínas Luminescentes/metabolismo , Microscopia Confocal , Modelos Teóricos , Proteína Vermelha FluorescenteRESUMO
Spleen tyrosine kinase (SYK) can behave as an oncogene or a tumor suppressor, depending on the cell and tissue type. As pharmacological SYK inhibitors are currently evaluated in clinical trials, it is important to gain more information on the molecular mechanisms underpinning these opposite roles. To this aim, we reconstructed and compared its signaling networks using phosphoproteomic data from breast cancer and Burkitt lymphoma cell lines where SYK behaves as a tumor suppressor and promoter. Bioinformatic analyses allowed for unveiling the main differences in signaling pathways, network topology and signal propagation from SYK to its potential effectors. In breast cancer cells, the SYK target-enriched signaling pathways included intercellular adhesion and Hippo signaling components that are often linked to tumor suppression. In Burkitt lymphoma cells, the SYK target-enriched signaling pathways included molecules that could play a role in SYK pro-oncogenic function in B-cell lymphomas. Several protein interactions were profoundly rewired in the breast cancer network compared with the Burkitt lymphoma network. These data demonstrate that proteomic profiling combined with mathematical network modeling allows untangling complex pathway interplays and revealing difficult to discern interactions among the SYK pathways that positively and negatively affect tumor formation and progression.
Assuntos
Neoplasias da Mama/metabolismo , Neoplasias da Mama/patologia , Linfoma de Burkitt/metabolismo , Linfoma de Burkitt/patologia , Quinase Syk/metabolismo , Neoplasias da Mama/genética , Linfoma de Burkitt/genética , Linhagem Celular Tumoral , Feminino , Humanos , Células MCF-7 , Modelos Teóricos , Fosfoproteínas/metabolismo , Proteômica , Transdução de Sinais/genética , Transdução de Sinais/fisiologia , Quinase Syk/genéticaAssuntos
Proteínas de Drosophila/fisiologia , Desenvolvimento Embrionário/fisiologia , Genoma , Proteínas Nucleares/fisiologia , Ativação Transcricional/fisiologia , Zigoto/crescimento & desenvolvimento , Animais , Ciclo Celular/fisiologia , Cromatina/metabolismo , Elementos Facilitadores Genéticos/fisiologia , Mitose/fisiologia , Organogênese/fisiologia , RNA Mensageiro/fisiologiaRESUMO
Organisms must ensure that expression of genes is directed to the appropriate tissues at the correct times, while simultaneously ensuring that these gene regulatory systems are robust to perturbation. This idea is captured by a mathematical concept called r-robustness, which says that a system is robust to a perturbation in up to r - 1 randomly chosen parameters. r-robustness implies that the biological system has a small number of sensitive parameters and that this number can be used as a robustness measure. In this work we use this idea to investigate the robustness of gene regulation using a sequence level model of the Drosophila melanogaster gene even-skipped. We consider robustness with respect to mutations of the enhancer sequence and with respect to changes of the transcription factor concentrations. We find that gene regulation is r-robust with respect to mutations in the enhancer sequence and identify a number of sensitive nucleotides. In both natural and in silico predicted enhancers, the number of nucleotides that are sensitive to mutation correlates negatively with the length of the sequence, meaning that longer sequences are more robust. The exact degree of robustness obtained is dependent not only on DNA sequence, but also on the local concentration of regulatory factors. We find that gene regulation can be remarkably sensitive to changes in transcription factor concentrations at the boundaries of expression features, while it is robust to perturbation elsewhere.
Assuntos
Elementos Facilitadores Genéticos/genética , Regulação da Expressão Gênica no Desenvolvimento/genética , Análise de Sequência de DNA/métodos , Animais , Sítios de Ligação/genética , Padronização Corporal/genética , Simulação por Computador , Proteínas de Drosophila/genética , Drosophila melanogaster/genética , Evolução Molecular , Proteínas de Homeodomínio/genética , Proteínas de Homeodomínio/metabolismo , Modelos Teóricos , Fatores de Transcrição/genética , Fatores de Transcrição/metabolismoRESUMO
Protein phosphorylation acts as an efficient switch controlling deregulated key signaling pathway in cancer. Computational biology aims to address the complexity of reconstructed networks but overrepresents well-known proteins and lacks information on less-studied proteins. A bioinformatic tool to reconstruct and select relatively small networks that connect signaling proteins to their targets in specific contexts is developed. It enables to propose and validate new signaling axes of the Syk kinase. To validate the potency of the tool, it is applied to two phosphoproteomic studies on oncogenic mutants of the well-known phosphatidyl-inositol 3-kinase (PIK3CA) and the unfamiliar Src-related tyrosine kinase lacking C-terminal regulatory tyrosine and N-terminal myristoylation sites (SRMS) kinase. By combining network reconstruction and signal propagation, comprehensive signaling networks from large-scale experimental data are built and multiple molecular paths from these kinases to their targets are extracted. Specific paths from two distinct PIK3CA mutants are retrieved, and their differential impact on the HER3 receptor kinase is explained. In addition, to address the missing connectivities of the SRMS kinase to its targets in interaction pathway databases, phospho-tyrosine and phospho-serine/threonine proteomic data are integrated. The resulting SRMS-signaling network comprises casein kinase 2, thereby validating its currently suggested role downstream of SRMS. The computational pipeline is publicly available, and contains a user-friendly graphical interface (http://doi.org/10.5281/zenodo.3333687).
Assuntos
Neoplasias/metabolismo , Proteômica , Transdução de Sinais , Linhagem Celular Tumoral , Humanos , Mutação/genética , Proteínas de Neoplasias/metabolismo , Fosforilação , Interface Usuário-ComputadorRESUMO
The original version of this Article contained an error in Fig. 4a, in which the "=" sign of the equation was inadvertently replaced with a "-" sign. This has been corrected in the PDF and HTML versions of the Article.
RESUMO
Although novel targeted therapies have significantly improved the overall survival of patients with advanced melanoma, understanding and combatting drug resistance remains a major clinical challenge. Using partial differential equations, we describe the evolution of a cellular population through time, space, and phenotype dimensions, in the presence of various drug species. We then use this framework to explore models in which resistance is attained by either mutations (irreversible) or plasticity (reversible). Numerical results suggest that punctuated evolutionary assumptions are more consistent with results obtained from murine melanoma models than gradual evolution. Furthermore, in the context of an evolving tumour cell population, sequencing the treatment, for instance applying immunotherapy before BRAF inhibitors, can increase treatment effectiveness. However, drug strategies which showed success within a spatially homogeneous tumour environment were unsuccessful under heterogeneous conditions, suggesting that spatio-environmental heterogeneity may be the greatest challenge to tumour therapies. Plastic metabolic models are additionally capable of reproducing the characteristic resistant tumour volume curves and predicting re-sensitisation to secondary waves of treatment observed in patient derived xenograft (PDX) melanomas treated with MEK and BRAF inhibitors. Nevertheless, secondary relapse due to a pre-adapted subpopulation, remaining after the first wave of treatment, results in a more rapid development of resistance. Our model provides a framework through which tumour resistance can be understood and would suggest that carefully phased treatments may be able to overcome the development of long-term resistance in melanoma.
Assuntos
Imunoterapia , Melanoma , Modelos Biológicos , Mutação , Recidiva Local de Neoplasia , Inibidores de Proteínas Quinases/uso terapêutico , Animais , Humanos , Melanoma/genética , Melanoma/imunologia , Melanoma/metabolismo , Melanoma/terapia , Camundongos , Recidiva Local de Neoplasia/genética , Recidiva Local de Neoplasia/imunologia , Recidiva Local de Neoplasia/metabolismo , Recidiva Local de Neoplasia/terapiaRESUMO
Pioneer transcription factors can engage nucleosomal DNA, which leads to local chromatin remodeling and to the establishment of transcriptional competence. However, the impact of enhancer priming by pioneer factors on the temporal control of gene expression and on mitotic memory remains unclear. Here we employ quantitative live imaging methods and mathematical modeling to test the effect of the pioneer factor Zelda on transcriptional dynamics and memory in Drosophila embryos. We demonstrate that increasing the number of Zelda binding sites accelerates the kinetics of nuclei transcriptional activation regardless of their transcriptional past. Despite its known pioneering activities, we show that Zelda does not remain detectably associated with mitotic chromosomes and is neither necessary nor sufficient to foster memory. We further reveal that Zelda forms sub-nuclear dynamic hubs where Zelda binding events are transient. We propose that Zelda facilitates transcriptional activation by accumulating in microenvironments where it could accelerate the duration of multiple pre-initiation steps.
Assuntos
Proteínas de Drosophila/metabolismo , Drosophila/metabolismo , Regulação da Expressão Gênica no Desenvolvimento , Fatores de Transcrição/metabolismo , Animais , Núcleo Celular/genética , Núcleo Celular/metabolismo , Drosophila/citologia , Drosophila/enzimologia , Drosophila/genética , Proteínas de Drosophila/genética , Cinética , Mitose , Proteínas Nucleares , Fatores de Transcrição/genética , Transcrição Gênica , Ativação TranscricionalRESUMO
During development, transcriptional properties of progenitor cells are stably propagated across multiple cellular divisions. Yet, at each division, chromatin faces structural constraints imposed by the important nuclear re-organization operating during mitosis. It is now clear that not all transcriptional regulators are ejected during mitosis, but rather that a subset of transcription factors, chromatin regulators and epigenetic histone marks are able to 'bookmark' specific loci, thereby providing a mitotic memory. Here we review mechanisms of mitotic bookmarking and discuss their impact on transcriptional dynamics in the context of multicellular developing embryos. We document recent discoveries and technological advances, and present current mathematical models of short-term transcriptional memory.
RESUMO
Mitosis is induced by the activation of the cyclin B/cdk1 feedback loop that creates a bistable state. The triggering factor promoting active cyclin B/cdk1 switch has been assigned to cyclin B/cdk1 accumulation during G2. However, this complex is rapidly inactivated by Wee1/Myt1-dependent phosphorylation of cdk1 making unlikely a triggering role of this kinase in mitotic commitment. Here we show that cyclin A/cdk1 kinase is the factor triggering mitosis. Cyclin A/cdk1 phosphorylates Bora to promote Aurora A-dependent Plk1 phosphorylation and activation and mitotic entry. We demonstrate that Bora phosphorylation by cyclin A/cdk1 is both necessary and sufficient for mitotic commitment. Finally, we identify a site in Bora whose phosphorylation by cyclin A/cdk1 is required for mitotic entry. We constructed a mathematical model confirming the essential role of this kinase in mitotic commitment. Overall, our results uncover the molecular mechanism by which cyclin A/cdk1 triggers mitotic entry.
Assuntos
Proteína Quinase CDC2/metabolismo , Proteínas de Ciclo Celular/metabolismo , Ciclina A/metabolismo , Mitose/fisiologia , Proteínas de Xenopus/metabolismo , Xenopus laevis/metabolismo , Animais , Proteína Quinase CDC2/genética , Proteínas de Ciclo Celular/genética , Ciclina A/genética , Ativação Enzimática , Feminino , Modelos Teóricos , Fosforilação , Proteínas de Xenopus/genética , Xenopus laevis/genéticaRESUMO
Sensing and reciprocating cellular systems (SARs) are important for the operation of many biological systems. Production in interferon (IFN) SARs is achieved through activation of the Jak-Stat pathway, and downstream upregulation of IFN regulatory factor (IRF)-7 and IFN transcription, but the role that high- and low-affinity IFNs play in this process remains unclear. We present a comparative between a minimal spatio-temporal partial differential equation model and a novel spatio-structural-temporal (SST) model for the consideration of receptor, binding, and metabolic aspects of SAR behaviour. Using the SST framework, we simulate single- and multi-cluster paradigms of IFN communication. Simulations reveal a cyclic process between the binding of IFN to the receptor, and the consequent increase in metabolism, decreasing the propensity for binding due to the internal feedback mechanism. One observes the effect of heterogeneity between cellular clusters, allowing them to individualise and increase local production, and within clusters, where we observe 'subpopular quiescence'; a process whereby intra-cluster subpopulations reduce their binding and metabolism such that other such subpopulations may augment their production. Finally, we observe the ability for low-affinity IFN to communicate a long range signal, where high affinity cannot, and the breakdown of this relationship through the introduction of cell motility. Biological systems may utilise cell motility where environments are unrestrictive and may use fixed system, with low-affinity communication, where a localised response is desirable.
Assuntos
Modelos Biológicos , Transdução de Sinais/fisiologia , Animais , Comunicação Celular , Simulação por Computador , Humanos , Interferons/metabolismo , Ligantes , Conceitos Matemáticos , Redes e Vias Metabólicas , Receptores de Interferon/metabolismo , Análise Espaço-TemporalRESUMO
We consider continuous-time recurrent neural networks as dynamical models for the simulation of human body motions. These networks consist of a few centers and many satellites connected to them. The centers evolve in time as periodical oscillators with different frequencies. The center states define the satellite neurons' states by a radial basis function (RBF) network. To simulate different motions, we adjust the parameters of the RBF networks. Our network includes a switching module that allows for turning from one motion to another. Simulations show that this model allows us to simulate complicated motions consisting of many different dynamical primitives. We also use the model for learning human body motion from markers' trajectories. We find that center frequencies can be learned from a small number of markers and can be transferred to other markers, such that our technique seems to be capable of correcting for missing information resulting from sparse control marker settings.
Assuntos
Movimento (Física) , Algoritmos , Corpo Humano , Humanos , Aprendizagem , Redes Neurais de Computação , NeurôniosRESUMO
MOTIVATION: Integration of metabolic networks with '-omics' data has been a subject of recent research in order to better understand the behaviour of such networks with respect to differences between biological and clinical phenotypes. Under the conditions of steady state of the reaction network and the non-negativity of fluxes, metabolic networks can be algebraically decomposed into a set of sub-pathways often referred to as extreme currents (ECs). Our objective is to find the statistical association of such sub-pathways with given clinical outcomes, resulting in a particular instance of a self-contained gene set analysis method. In this direction, we propose a method based on sparse group lasso (SGL) to identify phenotype associated ECs based on gene expression data. SGL selects a sparse set of feature groups and also introduces sparsity within each group. Features in our model are clusters of ECs, and feature groups are defined based on correlations among these features. RESULTS: We apply our method to metabolic networks from KEGG database and study the association of network features to prostate cancer (where the outcome is tumor and normal, respectively) as well as glioblastoma multiforme (where the outcome is survival time). In addition, simulations show the superior performance of our method compared to global test, which is an existing self-contained gene set analysis method. AVAILABILITY AND IMPLEMENTATION: R code (compatible with version 3.2.5) is available from http://www.abi.bit.uni-bonn.de/index.php?id=17. CONTACT: samal@combine.rwth-aachen.de or frohlich@bit.uni-bonn.de. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.