High-throughput data analysis

Jianrong Lu, Marie-Laure Ruhf, Norbert Perrimon, and Philip Leder. 2007. “A genome-wide RNA interference screen identifies putative chromatin regulators essential for E2F repression.” Proc Natl Acad Sci U S A, 104, 22, Pp. 9381-6.Abstract

Regulation of chromatin structure is critical in many fundamental cellular processes. Previous studies have suggested that the Rb tumor suppressor may recruit multiple chromatin regulatory proteins to repress E2F, a key regulator of cell proliferation and differentiation. Taking advantage of the evolutionary conservation of the E2F pathway, we have conducted a genome-wide RNAi screen in cultured Drosophila cells for genes required for repression of E2F activity. Among the genes identified are components of the putative Domino chromatin remodeling complex, as well as the Polycomb Group (PcG) protein-like fly tumor suppressor, L3mbt, and the related CG16975/dSfmbt. These factors are recruited to E2F-responsive promoters through physical association with E2F and are required for repression of endogenous E2F target genes. Surprisingly, their inhibitory activities on E2F appear to be independent of Rb. In Drosophila, domino mutation enhances cell proliferation induced by E2F overexpression and suppresses a loss-of-function cyclin E mutation. These findings suggest that potential chromatin regulation mediated by Domino and PcG-like factors plays an important role in controlling E2F activity and cell growth.

Dawei Jiang, Linlin Zhao, and David E Clapham. 2009. “Genome-wide RNAi screen identifies Letm1 as a mitochondrial Ca2+/H+ antiporter.” Science, 326, 5949, Pp. 144-7.Abstract

Mitochondria are integral components of cellular calcium (Ca2+) signaling. Calcium stimulates mitochondrial adenosine 5'-triphosphate production, but can also initiate apoptosis. In turn, cytoplasmic Ca2+ concentrations are regulated by mitochondria. Although several transporter and ion-channel mechanisms have been measured in mitochondria, the molecules that govern Ca2+ movement across the inner mitochondrial membrane are unknown. We searched for genes that regulate mitochondrial Ca2+ and H+ concentrations using a genome-wide Drosophila RNA interference (RNAi) screen. The mammalian homolog of one Drosophila gene identified in the screen, Letm1, was found to specifically mediate coupled Ca2+/H+ exchange. RNAi knockdown, overexpression, and liposome reconstitution of the purified Letm1 protein demonstrate that Letm1 is a mitochondrial Ca2+/H+ antiporter.

Marcelo Perez-Pepe, Victoria Slomiansky, Mariela Loschi, Luciana Luchelli, Maximiliano Neme, María Gabriela Thomas, and Graciela Lidia Boccaccio. 2012. “BUHO: a MATLAB script for the study of stress granules and processing bodies by high-throughput image analysis.” PLoS One, 7, 12, Pp. e51495.Abstract

The spontaneous and reversible formation of foci and filaments that contain proteins involved in different metabolic processes is common in both the nucleus and the cytoplasm. Stress granules (SGs) and processing bodies (PBs) belong to a novel family of cellular structures collectively known as mRNA silencing foci that harbour repressed mRNAs and their associated proteins. SGs and PBs are highly dynamic and they form upon stress and dissolve thus releasing the repressed mRNAs according to changes in cell physiology. In addition, aggregates containing abnormal proteins are frequent in neurodegenerative disorders. In spite of the growing relevance of these supramolecular aggregates to diverse cellular functions a reliable automated tool for their systematic analysis is lacking. Here we report a MATLAB Script termed BUHO for the high-throughput image analysis of cellular foci. We used BUHO to assess the number, size and distribution of distinct objects with minimal deviation from manually obtained parameters. BUHO successfully addressed the induction of both SGs and PBs in mammalian and insect cells exposed to different stress stimuli. We also used BUHO to assess the dynamics of specific mRNA-silencing foci termed Smaug 1 foci (S-foci) in primary neurons upon synaptic stimulation. Finally, we used BUHO to analyze the role of candidate genes on SG formation in an RNAi-based experiment. We found that FAK56D, GCN2 and PP1 govern SG formation. The role of PP1 is conserved in mammalian cells as judged by the effect of the PP1 inhibitor salubrinal, and involves dephosphorylation of the translation factor eIF2α. All these experiments were analyzed manually and by BUHO and the results differed in less than 5% of the average value. The automated analysis by this user-friendly method will allow high-throughput image processing in short times by providing a robust, flexible and reliable alternative to the laborious and sometimes unfeasible visual scrutiny.

Ian Flockhart, Matthew Booker, Amy Kiger, Michael Boutros, Susan Armknecht, Nadire Ramadan, Kris Richardson, Andrew Xu, Norbert Perrimon, and Bernard Mathey-Prevot. 2006. “FlyRNAi: the Drosophila RNAi screening center database.” Nucleic Acids Res, 34, Database issue, Pp. D489-94.Abstract

RNA interference (RNAi) has become a powerful tool for genetic screening in Drosophila. At the Drosophila RNAi Screening Center (DRSC), we are using a library of over 21,000 double-stranded RNAs targeting known and predicted genes in Drosophila. This library is available for the use of visiting scientists wishing to perform full-genome RNAi screens. The data generated from these screens are collected in the DRSC database (http://flyRNAi.org/cgi-bin/RNAi_screens.pl) in a flexible format for the convenience of the scientist and for archiving data. The long-term goal of this database is to provide annotations for as many of the uncharacterized genes in Drosophila as possible. Data from published screens are available to the public through a highly configurable interface that allows detailed examination of the data and provides access to a number of other databases and bioinformatics tools.

Amy M Wiles, Dashnamoorthy Ravi, Selvaraj Bhavani, and Alexander JR Bishop. 2008. “An analysis of normalization methods for Drosophila RNAi genomic screens and development of a robust validation scheme.” J Biomol Screen, 13, 8, Pp. 777-84.Abstract

Genome-wide RNA interference (RNAi) screening allows investigation of the role of individual genes in a process of choice. Most RNAi screens identify a large number of genes with a continuous gradient in the assessed phenotype. Screeners must decide whether to examine genes with the most robust phenotype or the full gradient of genes that cause an effect and how to identify candidate genes. The authors have used RNAi in Drosophila cells to examine viability in a 384-well plate format and compare 2 screens, untreated control and treatment. They compare multiple normalization methods, which take advantage of different features within the data, including quantile normalization, background subtraction, scaling, cellHTS2 (Boutros et al. 2006), and interquartile range measurement. Considering the false-positive potential that arises from RNAi technology, a robust validation method was designed for the purpose of gene selection for future investigations. In a retrospective analysis, the authors describe the use of validation data to evaluate each normalization method. Although no method worked ideally, a combination of 2 methods, background subtraction followed by quantile normalization and cellHTS2, at different thresholds, captures the most dependable and diverse candidate genes. Thresholds are suggested depending on whether a few candidate genes are desired or a more extensive systems-level analysis is sought. The normalization approaches and experimental design to perform validation experiments are likely to apply to those high-throughput screening systems attempting to identify genes for systems-level analysis.

Ralph A Neumüller and Norbert Perrimon. 2011. “Where gene discovery turns into systems biology: genome-scale RNAi screens in Drosophila.” Wiley Interdiscip Rev Syst Biol Med, 3, 4, Pp. 471-8.Abstract

Systems biology aims to describe the complex interplays between cellular building blocks which, in their concurrence, give rise to the emergent properties observed in cellular behaviors and responses. This approach tries to determine the molecular players and the architectural principles of their interactions within the genetic networks that control certain biological processes. Large-scale loss-of-function screens, applicable in various different model systems, have begun to systematically interrogate entire genomes to identify the genes that contribute to a certain cellular response. In particular, RNA interference (RNAi)-based high-throughput screens have been instrumental in determining the composition of regulatory systems and paired with integrative data analyses have begun to delineate the genetic networks that control cell biological and developmental processes. Through the creation of tools for both, in vitro and in vivo genome-wide RNAi screens, Drosophila melanogaster has emerged as one of the key model organisms in systems biology research and over the last years has massively contributed to and hence shaped this discipline. WIREs Syst Biol Med 2011 3 471-478 DOI: 10.1002/wsbm.127

AA Kiger, B Baum, S Jones, MR Jones, A Coulson, C Echeverri, and N Perrimon. 2003. “A functional genomic analysis of cell morphology using RNA interference.” J Biol, 2, 4, Pp. 27.Abstract

BACKGROUND: The diversity of metazoan cell shapes is influenced by the dynamic cytoskeletal network. With the advent of RNA-interference (RNAi) technology, it is now possible to screen systematically for genes controlling specific cell-biological processes, including those required to generate distinct morphologies. RESULTS: We adapted existing RNAi technology in Drosophila cell culture for use in high-throughput screens to enable a comprehensive genetic dissection of cell morphogenesis. To identify genes responsible for the characteristic shape of two morphologically distinct cell lines, we performed RNAi screens in each line with a set of double-stranded RNAs (dsRNAs) targeting 994 predicted cell shape regulators. Using automated fluorescence microscopy to visualize actin filaments, microtubules and DNA, we detected morphological phenotypes for 160 genes, one-third of which have not been previously characterized in vivo. Genes with similar phenotypes corresponded to known components of pathways controlling cytoskeletal organization and cell shape, leading us to propose similar functions for previously uncharacterized genes. Furthermore, we were able to uncover genes acting within a specific pathway using a co-RNAi screen to identify dsRNA suppressors of a cell shape change induced by Pten dsRNA. CONCLUSIONS: Using RNAi, we identified genes that influence cytoskeletal organization and morphology in two distinct cell types. Some genes exhibited similar RNAi phenotypes in both cell types, while others appeared to have cell-type-specific functions, in part reflecting the different mechanisms used to generate a round or a flat cell morphology.

Ramanuj DasGupta, Kent Nybakken, Matthew Booker, Bernard Mathey-Prevot, Foster Gonsalves, Binita Changkakoty, and Norbert Perrimon. 2007. “A case study of the reproducibility of transcriptional reporter cell-based RNAi screens in Drosophila.” Genome Biol, 8, 9, Pp. R203.Abstract

Off-target effects have been demonstrated to be a major source of false-positives in RNA interference (RNAi) high-throughput screens. In this study, we re-assess the previously published transcriptional reporter-based whole-genome RNAi screens for the Wingless and Hedgehog signaling pathways using second generation double-stranded RNA libraries. Furthermore, we investigate other factors that may influence the outcome of such screens, including cell-type specificity, robustness of reporters, and assay normalization, which determine the efficacy of RNAi-knockdown of target genes.

Amy M Wiles, Mark Doderer, Jianhua Ruan, Ting-Ting Gu, Dashnamoorthy Ravi, Barron Blackman, and Alexander JR Bishop. 2010. “Building and analyzing protein interactome networks by cross-species comparisons.” BMC Syst Biol, 4, Pp. 36.Abstract

BACKGROUND: A genomic catalogue of protein-protein interactions is a rich source of information, particularly for exploring the relationships between proteins. Numerous systems-wide and small-scale experiments have been conducted to identify interactions; however, our knowledge of all interactions for any one species is incomplete, and alternative means to expand these network maps is needed. We therefore took a comparative biology approach to predict protein-protein interactions across five species (human, mouse, fly, worm, and yeast) and developed InterologFinder for research biologists to easily navigate this data. We also developed a confidence score for interactions based on available experimental evidence and conservation across species. RESULTS: The connectivity of the resultant networks was determined to have scale-free distribution, small-world properties, and increased local modularity, indicating that the added interactions do not disrupt our current understanding of protein network structures. We show examples of how these improved interactomes can be used to analyze a genome-scale dataset (RNAi screen) and to assign new function to proteins. Predicted interactions within this dataset were tested by co-immunoprecipitation, resulting in a high rate of validation, suggesting the high quality of networks produced. CONCLUSIONS: Protein-protein interactions were predicted in five species, based on orthology. An InteroScore, a score accounting for homology, number of orthologues with evidence of interactions, and number of unique observations of interactions, is given to each known and predicted interaction. Our website http://www.interologfinder.org provides research biologists intuitive access to this data.

Clemens Bergwitz, Mark J Wee, Sumi Sinha, Joanne Huang, Charles DeRobertis, Lawrence B Mensah, Jonathan Cohen, Adam Friedman, Meghana Kulkarni, Yanhui Hu, Arunachalam Vinayagam, Michael Schnall-Levin, Bonnie Berger, Lizabeth A Perkins, Stephanie E Mohr, and Norbert Perrimon. 2013. “Genetic determinants of phosphate response in Drosophila.” PLoS One, 8, 3, Pp. e56753.Abstract

Phosphate is required for many important cellular processes and having too little phosphate or too much can cause disease and reduce life span in humans. However, the mechanisms underlying homeostatic control of extracellular phosphate levels and cellular effects of phosphate are poorly understood. Here, we establish Drosophila melanogaster as a model system for the study of phosphate effects. We found that Drosophila larval development depends on the availability of phosphate in the medium. Conversely, life span is reduced when adult flies are cultured on high phosphate medium or when hemolymph phosphate is increased in flies with impaired malpighian tubules. In addition, RNAi-mediated inhibition of MAPK-signaling by knockdown of Ras85D, phl/D-Raf or Dsor1/MEK affects larval development, adult life span and hemolymph phosphate, suggesting that some in vivo effects involve activation of this signaling pathway by phosphate. To identify novel genetic determinants of phosphate responses, we used Drosophila hemocyte-like cultured cells (S2R+) to perform a genome-wide RNAi screen using MAPK activation as the readout. We identified a number of candidate genes potentially important for the cellular response to phosphate. Evaluation of 51 genes in live flies revealed some that affect larval development, adult life span and hemolymph phosphate levels.

Yousang Gwack, Sonia Sharma, Julie Nardone, Bogdan Tanasa, Alina Iuga, Sonal Srikanth, Heidi Okamura, Diana Bolton, Stefan Feske, Patrick G Hogan, and Anjana Rao. 2006. “A genome-wide Drosophila RNAi screen identifies DYRK-family kinases as regulators of NFAT.” Nature, 441, 7093, Pp. 646-50.Abstract

Precise regulation of the NFAT (nuclear factor of activated T cells) family of transcription factors (NFAT1-4) is essential for vertebrate development and function. In resting cells, NFAT proteins are heavily phosphorylated and reside in the cytoplasm; in cells exposed to stimuli that raise intracellular free Ca2+ levels, they are dephosphorylated by the calmodulin-dependent phosphatase calcineurin and translocate to the nucleus. NFAT dephosphorylation by calcineurin is countered by distinct NFAT kinases, among them casein kinase 1 (CK1) and glycogen synthase kinase 3 (GSK3). Here we have used a genome-wide RNA interference (RNAi) screen in Drosophila to identify additional regulators of the signalling pathway leading from Ca2+-calcineurin to NFAT. This screen was successful because the pathways regulating NFAT subcellular localization (Ca2+ influx, Ca2+-calmodulin-calcineurin signalling and NFAT kinases) are conserved across species, even though Ca2+-regulated NFAT proteins are not themselves represented in invertebrates. Using the screen, we have identified DYRKs (dual-specificity tyrosine-phosphorylation regulated kinases) as novel regulators of NFAT. DYRK1A and DYRK2 counter calcineurin-mediated dephosphorylation of NFAT1 by directly phosphorylating the conserved serine-proline repeat 3 (SP-3) motif of the NFAT regulatory domain, thus priming further phosphorylation of the SP-2 and serine-rich region 1 (SRR-1) motifs by GSK3 and CK1, respectively. Thus, genetic screening in Drosophila can be successfully applied to cross evolutionary boundaries and identify new regulators of a transcription factor that is expressed only in vertebrates.

Sylvia Erhardt, Barbara G Mellone, Craig M Betts, Weiguo Zhang, Gary H Karpen, and Aaron F Straight. 2008. “Genome-wide analysis reveals a cell cycle-dependent mechanism controlling centromere propagation.” J Cell Biol, 183, 5, Pp. 805-18.Abstract

Centromeres are the structural and functional foundation for kinetochore formation, spindle attachment, and chromosome segregation. In this study, we isolated factors required for centromere propagation using genome-wide RNA interference screening for defects in centromere protein A (CENP-A; centromere identifier [CID]) localization in Drosophila melanogaster. We identified the proteins CAL1 and CENP-C as essential factors for CID assembly at the centromere. CID, CAL1, and CENP-C coimmunoprecipitate and are mutually dependent for centromere localization and function. We also identified the mitotic cyclin A (CYCA) and the anaphase-promoting complex (APC) inhibitor RCA1/Emi1 as regulators of centromere propagation. We show that CYCA is centromere localized and that CYCA and RCA1/Emi1 couple centromere assembly to the cell cycle through regulation of the fizzy-related/CDH1 subunit of the APC. Our findings identify essential components of the epigenetic machinery that ensures proper specification and propagation of the centromere and suggest a mechanism for coordinating centromere inheritance with cell division.

Jennifer L Rohn, David Sims, Tao Liu, Marina Fedorova, Frieder Schöck, Joseph Dopie, Maria K Vartiainen, Amy A Kiger, Norbert Perrimon, and Buzz Baum. 2011. “Comparative RNAi screening identifies a conserved core metazoan actinome by phenotype.” J Cell Biol, 194, 5, Pp. 789-805.Abstract

Although a large number of actin-binding proteins and their regulators have been identified through classical approaches, gaps in our knowledge remain. Here, we used genome-wide RNA interference as a systematic method to define metazoan actin regulators based on visual phenotype. Using comparative screens in cultured Drosophila and human cells, we generated phenotypic profiles for annotated actin regulators together with proteins bearing predicted actin-binding domains. These phenotypic clusters for the known metazoan "actinome" were used to identify putative new core actin regulators, together with a number of genes with conserved but poorly studied roles in the regulation of the actin cytoskeleton, several of which we studied in detail. This work suggests that although our search for new components of the core actin machinery is nearing saturation, regulation at the level of nuclear actin export, RNA splicing, ubiquitination, and other upstream processes remains an important but unexplored frontier of actin biology.

Susan Armknecht, Michael Boutros, Amy Kiger, Kent Nybakken, Bernard Mathey-Prevot, and Norbert Perrimon. 2005. “High-throughput RNA interference screens in Drosophila tissue culture cells.” Methods Enzymol, 392, Pp. 55-73.Abstract

This chapter describes the method used to conduct high-throughput screening (HTs) by RNA interference in Drosophila tissue culture cells. It covers four main topics: (1) a brief description of the existing platforms to conduct RNAi-screens in cell-based assays; (2) a table of the Drosophila cell lines available for these screens and a brief mention of the need to establish other cell lines as well as cultures of primary cells; (3) a discussion of the considerations and protocols involved in establishing assays suitable for HTS in a 384-well format; and (A) a summary of the various ways of handling raw data from an ongoing screen, with special emphasis on how to apply normalization for experimental variation and statistical filters to sort out noise from signals.

Eric J Wagner, Brandon D Burch, Ashley C Godfrey, Harmony R Salzler, Robert J Duronio, and William F Marzluff. 2007. “A genome-wide RNA interference screen reveals that variant histones are necessary for replication-dependent histone pre-mRNA processing.” Mol Cell, 28, 4, Pp. 692-9.Abstract

Metazoan replication-dependent histone mRNAs are not polyadenylated and instead end in a conserved stem loop that is the cis element responsible for coordinate posttranscriptional regulation of these mRNAs. Using biochemical approaches, only a limited number of factors required for cleavage of histone pre-mRNA have been identified. We therefore performed a genome-wide RNA interference screen in Drosophila cells using a GFP reporter that is expressed only when histone pre-mRNA processing is disrupted. Four of the 24 genes identified encode proteins also necessary for cleavage/polyadenylation, indicating mechanistic conservation in formation of different mRNA 3' ends. We also unexpectedly identified the histone variants H2Av and H3.3A/B. In H2Av mutant cells, U7 snRNP remains active but fails to accumulate at the histone locus, suggesting there is a regulatory pathway that coordinates the production of variant and canonical histones that acts via localization of essential histone pre-mRNA processing factors.