High-resolution mapping of cancer cell networks using co-functional interactions.

Powerful new technologies for perturbing genetic elements have recently expanded the study of genetic interactions in model systems ranging from yeast to human cell lines. However, technical artifacts can confound signal across genetic screens and limit the immense potential of parallel screening approaches. To address this problem, we devised a novel PCA-based method for correcting genome-wide screening data, bolstering the sensitivity and specificity of detection for genetic interactions. Applying this strategy to a set of 436 whole genome CRISPR screens, we report more than 1.5 million pairs of correlated "co-functional" genes that provide finer-scale information about cell compartments, biological pathways, and protein complexes than traditional gene sets. Lastly, we employed a gene community detection approach to implicate core genes for cancer growth and compress signal from functionally related genes in the same community into a single score. This work establishes new algorithms for probing cancer cell networks and motivates the acquisition of further CRISPR screen data across diverse genotypes and cell types to further resolve complex cellular processes

Improved Network Performance via Antagonism: From Synthetic Rescues to Multi-drug Combinations

Recent research shows that a faulty or sub-optimally operating metabolic network can often be rescued by the targeted removal of enzyme-coding genes--the exact opposite of what traditional gene therapy would suggest. Predictions go as far as to assert that certain gene knockouts can restore the growth of otherwise nonviable gene-deficient cells. Many questions follow from this discovery: What are the underlying mechanisms? How generalizable is this effect? What are the potential applications? Here, I will approach these questions from the perspective of compensatory perturbations on networks. Relations will be drawn between such synthetic rescues and naturally occurring cascades of reaction inactivation, as well as their analogues in physical and other biological networks. I will specially discuss how rescue interactions can lead to the rational design of antagonistic drug combinations that select against resistance and how they can illuminate medical research on cancer, antibiotics, and metabolic diseases.Comment: Online Open "Problems and Paradigms" articl

Dispensability of Escherichia coli's latent pathways

Gene-knockout experiments on single-cell organisms have established that expression of a substantial fraction of genes is not needed for optimal growth. This problem acquired a new dimension with the recent discovery that environmental and genetic perturbations of the bacterium Escherichia coli are followed by the temporary activation of a large number of latent metabolic pathways, which suggests the hypothesis that temporarily activated reactions impact growth and hence facilitate adaptation in the presence of perturbations. Here we test this hypothesis computationally and find, surprisingly, that the availability of latent pathways consistently offers no growth advantage, and tends in fact to inhibit growth after genetic perturbations. This is shown to be true even for latent pathways with a known function in alternate conditions, thus extending the significance of this adverse effect beyond apparently nonessential genes. These findings raise the possibility that latent pathway activation is in fact derivative of another, potentially suboptimal, adaptive response

Flux imbalance analysis and the sensitivity of cellular growth to changes in metabolite pools

Stoichiometric models of metabolism, such as flux balance analysis (FBA), are classically applied to predicting steady state rates - or fluxes - of metabolic reactions in genome-scale metabolic networks. Here we revisit the central assumption of FBA, i.e. that intracellular metabolites are at steady state, and show that deviations from flux balance (i.e. flux imbalances) are informative of some features of in vivo metabolite concentrations. Mathematically, the sensitivity of FBA to these flux imbalances is captured by a native feature of linear optimization, the dual problem, and its corresponding variables, known as shadow prices. First, using recently published data on chemostat growth of Saccharomyces cerevisae under different nutrient limitations, we show that shadow prices anticorrelate with experimentally measured degrees of growth limitation of intracellular metabolites. We next hypothesize that metabolites which are limiting for growth (and thus have very negative shadow price) cannot vary dramatically in an uncontrolled way, and must respond rapidly to perturbations. Using a collection of published datasets monitoring the time-dependent metabolomic response of Escherichia coli to carbon and nitrogen perturbations, we test this hypothesis and find that metabolites with negative shadow price indeed show lower temporal variation following a perturbation than metabolites with zero shadow price. Finally, we illustrate the broader applicability of flux imbalance analysis to other constraint-based methods. In particular, we explore the biological significance of shadow prices in a constraint-based method for integrating gene expression data with a stoichiometric model. In this case, shadow prices point to metabolites that should rise or drop in concentration in order to increase consistency between flux predictions and gene expression data. In general, these results suggest that the sensitivity of metabolic optima to violations of the steady state constraints carries biologically significant information on the processes that control intracellular metabolites in the cell.Published versio

An integrative, multi-scale, genome-wide model reveals the phenotypic landscape of Escherichia coli.

Given the vast behavioral repertoire and biological complexity of even the simplest organisms, accurately predicting phenotypes in novel environments and unveiling their biological organization is a challenging endeavor. Here, we present an integrative modeling methodology that unifies under a common framework the various biological processes and their interactions across multiple layers. We trained this methodology on an extensive normalized compendium for the gram-negative bacterium Escherichia coli, which incorporates gene expression data for genetic and environmental perturbations, transcriptional regulation, signal transduction, and metabolic pathways, as well as growth measurements. Comparison with measured growth and high-throughput data demonstrates the enhanced ability of the integrative model to predict phenotypic outcomes in various environmental and genetic conditions, even in cases where their underlying functions are under-represented in the training set. This work paves the way toward integrative techniques that extract knowledge from a variety of biological data to achieve more than the sum of their parts in the context of prediction, analysis, and redesign of biological systems

Essential plasticity and redundancy of metabolism unveiled by synthetic lethality analysis

We unravel how functional plasticity and redundancy are essential mechanisms underlying the ability to survive of metabolic networks. We perform an exhaustive computational screening of synthetic lethal reaction pairs in Escherichia coli in a minimal medium and we find that synthetic lethal pairs divide in two different groups depending on whether the synthetic lethal interaction works as a backup or as a parallel use mechanism, the first corresponding to essential plasticity and the second to essential redundancy. In E. coli, the analysis of pathways entanglement through essential redundancy supports the view that synthetic lethality affects preferentially a single function or pathway. In contrast, essential plasticity, the dominant class, tends to be inter-pathway but strongly localized and unveils Cell Envelope Biosynthesis as an essential backup for Membrane Lipid Metabolism. When comparing E. coli and Mycoplasma pneumoniae, we find that the metabolic networks of the two organisms exhibit a large difference in the relative importance of plasticity and redundancy which is consistent with the conjecture that plasticity is a sophisticated mechanism that requires a complex organization. Finally, coessential reaction pairs are explored in different environmental conditions to uncover the interplay between the two mechanisms. We find that synthetic lethal interactions and their classification in plasticity and redundancy are basically insensitive to medium composition, and are highly conserved even when the environment is enriched with nonessential compounds or overconstrained to decrease maximum biomass formation.Comment: 22 pages, 4 figure

Mapping genetic interactions in cancer: a road to rational combination therapies.

The discovery of synthetic lethal interactions between poly (ADP-ribose) polymerase (PARP) inhibitors and BRCA genes, which are involved in homologous recombination, led to the approval of PARP inhibition as a monotherapy for patients with BRCA1/2-mutated breast or ovarian cancer. Studies following the initial observation of synthetic lethality demonstrated that the reach of PARP inhibitors is well beyond just BRCA1/2 mutants. Insights into the mechanisms of action of anticancer drugs are fundamental for the development of targeted monotherapies or rational combination treatments that will synergize to promote cancer cell death and overcome mechanisms of resistance. The development of targeted therapeutic agents is premised on mapping the physical and functional dependencies of mutated genes in cancer. An important part of this effort is the systematic screening of genetic interactions in a variety of cancer types. Until recently, genetic-interaction screens have relied either on the pairwise perturbations of two genes or on the perturbation of genes of interest combined with inhibition by commonly used anticancer drugs. Here, we summarize recent advances in mapping genetic interactions using targeted, genome-wide, and high-throughput genetic screens, and we discuss the therapeutic insights obtained through such screens. We further focus on factors that should be considered in order to develop a robust analysis pipeline. Finally, we discuss the integration of functional interaction data with orthogonal methods and suggest that such approaches will increase the reach of genetic-interaction screens for the development of rational combination therapies

Coupling metabolic footprinting and flux balance analysis to predict how single gene knockouts perturb microbial metabolism

Tese de mestrado. Biologia (Bioinformática e Biologia Computacional). Universidade de Lisboa, Faculdade de Ciências, 2012The model organisms Caenorhabditis elegans and E. coli form one of the simplest gut microbe host interaction models. Interventions in the microbe that increase the host longevity including inhibition of folate synthesis have been reported previously. To find novel single gene knockouts with an effect on lifespan, a screen of the Keio collection of E. coli was undertaken, and some of the genes found are directly involved in metabolism. The next step in those specific cases is to understand how these mutations perturb metabolism systematically, so that hypotheses can be generated. For that, I employed dynamic Flux Balance Analysis (dFBA), a constraint-based modeling technique capable of simulating the dynamics of metabolism in a batch culture and making predictions about changes in intracellular flux distribution. Since the specificities of the C. elegans lifespan experiments demand us to culture microbes in conditions differing from most of the published literature on E. coli physiology, novel data must be acquired to characterize and make dFBA simulations as realistic as possible. To do this exchange fluxes were measured using quantitative H NMR Time-Resolved Metabolic Footprinting. Furthermore, I also investigate the combination of TReF and dFBA as a tool in microbial metabolism studies. These approaches were tested by comparing wild type E. coli with one of the knockout strains found, ΔmetL, a knockout of the metL gene which encodes a byfunctional enzyme involved in aspartate and threonine metabolism. I found that the strain exhibits a slower growth rate than the wild type. Model simulation results revealed that reduced homoserine and methionine synthesis, as well as impaired sulfur and folate metabolism are the main effects of this knockout and the reasons for the growth deficiency. These results indicate that there are common mechanisms of the lifespan extension between ΔmetL and inhibition of folate biosynthesis and that the flux balance analysis/metabolic footprinting approach can help us understand the nature of these mechanisms.Os organismos modelo Caenorhabditis elegans e E. coli formam um dos modelos mais simples de interacções entre micróbio do tracto digestivo e hospedeiro. Intervenções no micróbio capazes de aumentar a longevidade do hospedeiro, incluindo inibição de síntese de folatos, foram reportadas previamente. Para encontrar novas delecções génicas do micróbio capazes de aumentar a longevidade do hospedeiro, a colecção Keio de deleções génicas de E. coli foi rastreada. Alguns dos genes encontrados participam em processos metabólicos, e nesses casos, o próximpo passo é perceber como as deleções perturbam o metabolismo sistémicamente, para gerar hipóteses. Para isso, utilizo dynamic Flux Balance Analysis (dFBA), uma técnica de modelação metabólica capaz de fazer previsões sobre alterações na distribuição intracelular de fluxos. As especificidades das experiências de tempo de vida em C.elegans obrigam-nos a trabalhar em condições diferentes das usadas na maioria da literatura publicada em fisiologia de E. coli, e para dar o máximo realismo às simulações de dFBA novos dados foram adquiridos, utilizando H NMR Time-Resolved Metabolic Footprinting para medir fluxos de troca de metabolitos entre microorganismo e meio de cultura. A combinação de TReF e dFBA como ferramenta de estudo do metabolism microbiano é também investigada. Estas abordagens foram testadas ao comparar E. coli wild-type com uma das estirpes encontradas no rastreio, ΔmetL, knockout do gene metL, que codifica um enzima bifunctional participante no metabolismo de aspartato e treonina, e que exibe uma taxa de crescimento reduzida comparativamente ao wild-type. Os resultados das simulações revelaram que os principais efeitos da deleção deste gene, e as razões para a menor taxa de crescimento observada, são a produção reduzida de homoserina e metionina e os efeitos que provoca no metabolismo de folatos e enxofre. Estes resultados indicam que há mecanismos comuns na extensão da longevidade causada por esta deleção e inibição de síntese de folatos, e que a combinação metabolic footprinting/flux balance analysis pode ajudar-nos a compreender a natureza desses mecanismos

How to understand the cell by breaking it: network analysis of gene perturbation screens

Modern high-throughput gene perturbation screens are key technologies at the forefront of genetic research. Combined with rich phenotypic descriptors they enable researchers to observe detailed cellular reactions to experimental perturbations on a genome-wide scale. This review surveys the current state-of-the-art in analyzing perturbation screens from a network point of view. We describe approaches to make the step from the parts list to the wiring diagram by using phenotypes for network inference and integrating them with complementary data sources. The first part of the review describes methods to analyze one- or low-dimensional phenotypes like viability or reporter activity; the second part concentrates on high-dimensional phenotypes showing global changes in cell morphology, transcriptome or proteome.Comment: Review based on ISMB 2009 tutorial; after two rounds of revisio