Analysis of high-throughput data increasingly relies on pathway annotation and functional information derived from Gene Ontology. This approach has limitations, in particular for the analysis of network dynamics over time or under different experimental conditions, in which modules within a network rather than complete pathways might respond and change. We report an analysis framework based on protein complexes, which are at the core of network reorganization. We generated a protein complex resource for human, Drosophila, and yeast from the literature and databases of protein-protein interaction networks, with each species having thousands of complexes. We developed COMPLEAT (http://www.flyrnai.org/compleat), a tool for data mining and visualization for complex-based analysis of high-throughput data sets, as well as analysis and integration of heterogeneous proteomics and gene expression data sets. With COMPLEAT, we identified dynamically regulated protein complexes among genome-wide RNA interference data sets that used the abundance of phosphorylated extracellular signal-regulated kinase in cells stimulated with either insulin or epidermal growth factor as the output. The analysis predicted that the Brahma complex participated in the insulin response.
The way in which cells adopt different morphologies is not fully understood. Cell shape could be a continuous variable or restricted to a set of discrete forms. We developed quantitative methods to describe cell shape and show that Drosophila haemocytes in culture are a heterogeneous mixture of five discrete morphologies. In an RNAi screen of genes affecting the morphological complexity of heterogeneous cell populations, we found that most genes regulate the transition between discrete shapes rather than generating new morphologies. In particular, we identified a subset of genes, including the tumour suppressor PTEN, that decrease the heterogeneity of the population, leading to populations enriched in rounded or elongated forms. We show that these genes have a highly conserved function as regulators of cell shape in both mouse and human metastatic melanoma cells.
RNA interference (RNAi) is a widely adopted tool for loss-of-function studies but RNAi results only have biological relevance if the reagents are appropriately mapped to genes. Several groups have designed and generated RNAi reagent libraries for studies in cells or in vivo for Drosophila and other species. At first glance, matching RNAi reagents to genes appears to be a simple problem, as each reagent is typically designed to target a single gene. In practice, however, the reagent-gene relationship is complex. Although the sequences of oligonucleotides used to generate most types of RNAi reagents are static, the reference genome and gene annotations are regularly updated. Thus, at the time a researcher chooses an RNAi reagent or analyzes RNAi data, the most current interpretation of the RNAi reagent-gene relationship, as well as related information regarding specificity (e.g., predicted off-target effects), can be different from the original interpretation. Here, we describe a set of strategies and an accompanying online tool, UP-TORR (for Updated Targets of RNAi Reagents; www.flyrnai.org/up-torr), useful for accurate and up-to-date annotation of cell-based and in vivo RNAi reagents. Importantly, UP-TORR automatically synchronizes with gene annotations daily, retrieving the most current information available, and for Drosophila, also synchronizes with the major reagent collections. Thus, UP-TORR allows users to choose the most appropriate RNAi reagents at the onset of a study, as well as to perform the most appropriate analyses of results of RNAi-based studies.
Besides its essential and well established role as a component of the cytoskeleton, actin is also present in the cell nucleus, where it has been linked to many processes that control gene expression. For example, nuclear actin regulates the activity of specific transcription factors, associates with all three RNA polymerases, and is a component of many chromatin remodelling complexes. Despite the fact that two export receptors, Crm1 and exportin 6, have been linked to nuclear export of actin, the mechanism by which actin enters the nucleus to elicit these essential functions has not been determined. It is also unclear whether actin is actively exchanged between the nucleus and the cytoplasm, and whether this connection has any functional significance for the cell. By applying a variety of live-cell imaging techniques we revealed that actin constantly shuttles in and out of the nucleus. The fast transport rates, which depend on the availability of actin monomers, suggest an active transport mechanism in both directions. Importantly, we identified importin 9 as the nuclear import factor for actin. Furthermore, our RNAi experiments showed that the active maintenance of nuclear actin levels by importin 9 is required for maximal transcriptional activity. Measurements of nuclear export rates and depletion studies also clarified that nuclear export of actin is mediated by exportin 6, and not by Crm1. These results demonstrate that cytoplasmic and nuclear actin pools are dynamically connected and identify the nuclear import and export mechanisms of actin.
When cells swell in hypo-osmotic solutions, chloride-selective ion channels (Cl(swell)) activate to reduce intracellular osmolality and prevent catastrophic cell rupture. Despite intensive efforts to assign a molecular identity to the mammalian Cl(swell) channel, it remains unknown. In an unbiased genome-wide RNA interference (RNAi) screen of Drosophila cells stably expressing an anion-sensitive fluorescent indicator, we identify Bestrophin 1 (dBest1) as the Drosophila Cl(swell) channel. Of the 23 screen hits with mammalian homologs and predicted transmembrane domains, only RNAi specifically targeting dBest1 eliminated the Cl(swell) current (I(Clswell)). We further demonstrate the essential contribution of dBest1 to Drosophila I(Clswell) with the introduction of a human Bestrophin disease-associated mutation (W94C). Overexpression of the W94C construct in Drosophila cells significantly reduced the endogenous I(Clswell). We confirm that exogenous expression of dBest1 alone in human embryonic kidney (HEK293) cells creates a clearly identifiable Drosophila-like I(Clswell). In contrast, activation of mouse Bestrophin 2 (mBest2), the closest mammalian ortholog of dBest1, is swell-insensitive. The first 64 residues of dBest1 conferred swell activation to mBest2. The chimera, however, maintains mBest2-like pore properties, strongly indicating that the Bestrophin protein forms the Cl(swell) channel itself rather than functioning as an essential auxiliary subunit. dBest1 is an anion channel clearly responsive to swell; this activation depends upon its N-terminus.
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.
Regulation of genes that initiate and amplify inflammatory programs of gene expression is achieved by signal-dependent exchange of coregulator complexes that function to read, write, and erase specific histone modifications linked to transcriptional activation or repression. Here, we provide evidence for the role of trimethylated histone H4 lysine 20 (H4K20me3) as a repression checkpoint that restricts expression of toll-like receptor 4 (TLR4) target genes in macrophages. H4K20me3 is deposited at the promoters of a subset of these genes by the SMYD5 histone methyltransferase through its association with NCoR corepressor complexes. Signal-dependent erasure of H4K20me3 is required for effective gene activation and is achieved by NF-κB-dependent delivery of the histone demethylase PHF2. Liver X receptors antagonize TLR4-dependent gene activation by maintaining NCoR/SMYD5-mediated repression. These findings reveal a histone H4K20 trimethylation/demethylation strategy that integrates positive and negative signaling inputs that control immunity and homeostasis.
Reactive Oxygen Species (ROS) are a natural by-product of cellular growth and proliferation, and are required for fundamental processes such as protein-folding and signal transduction. However, ROS accumulation, and the onset of oxidative stress, can negatively impact cellular and genomic integrity. Signalling networks have evolved to respond to oxidative stress by engaging diverse enzymatic and non-enzymatic antioxidant mechanisms to restore redox homeostasis. The architecture of oxidative stress response networks during periods of normal growth, and how increased ROS levels dynamically reconfigure these networks are largely unknown. In order to gain insight into the structure of signalling networks that promote redox homeostasis we first performed genome-scale RNAi screens to identify novel suppressors of superoxide accumulation. We then infer relationships between redox regulators by hierarchical clustering of phenotypic signatures describing how gene inhibition affects superoxide levels, cellular viability, and morphology across different genetic backgrounds. Genes that cluster together are likely to act in the same signalling pathway/complex and thus make "functional interactions". Moreover we also calculate differential phenotypic signatures describing the difference in cellular phenotypes following RNAi between untreated cells and cells submitted to oxidative stress. Using both phenotypic signatures and differential signatures we construct a network model of functional interactions that occur between components of the redox homeostasis network, and how such interactions become rewired in the presence of oxidative stress. This network model predicts a functional interaction between the transcription factor Jun and the IRE1 kinase, which we validate in an orthogonal assay. We thus demonstrate the ability of systems-biology approaches to identify novel signalling events.
FlyRNAi (http://www.flyrnai.org), the database and website of the Drosophila RNAi Screening Center (DRSC) at Harvard Medical School, serves a dual role, tracking both production of reagents for RNA interference (RNAi) screening in Drosophila cells and RNAi screen results. The database and website is used as a platform for community availability of protocols, tools, and other resources useful to researchers planning, conducting, analyzing or interpreting the results of Drosophila RNAi screens. Based on our own experience and user feedback, we have made several changes. Specifically, we have restructured the database to accommodate new types of reagents; added information about new RNAi libraries and other reagents; updated the user interface and website; and added new tools of use to the Drosophila community and others. Overall, the result is a more useful, flexible and comprehensive website and database.
Hedgehog (Hh) proteins are secreted molecules that function as organizers in animal development. In addition to being palmitoylated, Hh is the only metazoan protein known to possess a covalently-linked cholesterol moiety. The absence of either modification severely disrupts the organization of numerous tissues during development. It is currently not known how lipid-modified Hh is secreted and released from producing cells. We have performed a genome-wide RNAi screen in Drosophila melanogaster cells to identify regulators of Hh secretion. We found that cholesterol-modified Hh secretion is strongly dependent on coat protein complex I (COPI) but not COPII vesicles, suggesting that cholesterol modification alters the movement of Hh through the early secretory pathway. We provide evidence that both proteolysis and cholesterol modification are necessary for the efficient trafficking of Hh through the ER and Golgi. Finally, we identified several putative regulators of protein secretion and demonstrate a role for some of these genes in Hh and Wingless (Wg) morphogen secretion in vivo. These data open new perspectives for studying how morphogen secretion is regulated, as well as provide insight into regulation of lipid-modified protein secretion.
Sex chromosome dosage compensation in Drosophila provides a model for understanding how chromatin organization can modulate coordinate gene regulation. Male Drosophila increase the transcript levels of genes on the single male X approximately two-fold to equal the gene expression in females, which have two X-chromosomes. Dosage compensation is mediated by the Male-Specific Lethal (MSL) histone acetyltransferase complex. Five core components of the MSL complex were identified by genetic screens for genes that are specifically required for male viability and are dispensable for females. However, because dosage compensation must interface with the general transcriptional machinery, it is likely that identifying additional regulators that are not strictly male-specific will be key to understanding the process at a mechanistic level. Such regulators would not have been recovered from previous male-specific lethal screening strategies. Therefore, we have performed a cell culture-based, genome-wide RNAi screen to search for factors required for MSL targeting or function. Here we focus on the discovery of proteins that function to promote MSL complex recruitment to "chromatin entry sites," which are proposed to be the initial sites of MSL targeting. We find that components of the NSL (Non-specific lethal) complex, and a previously unstudied zinc-finger protein, facilitate MSL targeting and display a striking enrichment at MSL entry sites. Identification of these factors provides new insight into how MSL complex establishes the specialized hyperactive chromatin required for dosage compensation in Drosophila.
The pairing of homologous chromosomes is a fundamental feature of the meiotic cell. In addition, a number of species exhibit homolog pairing in nonmeiotic, somatic cells as well, with evidence for its impact on both gene regulation and double-strand break (DSB) repair. An extreme example of somatic pairing can be observed in Drosophila melanogaster, where homologous chromosomes remain aligned throughout most of development. However, our understanding of the mechanism of somatic homolog pairing remains unclear, as only a few genes have been implicated in this process. In this study, we introduce a novel high-throughput fluorescent in situ hybridization (FISH) technology that enabled us to conduct a genome-wide RNAi screen for factors involved in the robust somatic pairing observed in Drosophila. We identified both candidate "pairing promoting genes" and candidate "anti-pairing genes," providing evidence that pairing is a dynamic process that can be both enhanced and antagonized. Many of the genes found to be important for promoting pairing are highly enriched for functions associated with mitotic cell division, suggesting a genetic framework for a long-standing link between chromosome dynamics during mitosis and nuclear organization during interphase. In contrast, several of the candidate anti-pairing genes have known interphase functions associated with S-phase progression, DNA replication, and chromatin compaction, including several components of the condensin II complex. In combination with a variety of secondary assays, these results provide insights into the mechanism and dynamics of somatic pairing.
RNA interference (RNAi) leads to sequence-specific knockdown of gene function. The approach can be used in large-scale screens to interrogate function in various model organisms and an increasing number of other species. Genome-scale RNAi screens are routinely performed in cultured or primary cells or in vivo in organisms such as C. elegans. High-throughput RNAi screening is benefitting from the development of sophisticated new instrumentation and software tools for collecting and analyzing data, including high-content image data. The results of large-scale RNAi screens have already proved useful, leading to new understandings of gene function relevant to topics such as infection, cancer, obesity, and aging. Nevertheless, important caveats apply and should be taken into consideration when developing or interpreting RNAi screens. Some level of false discovery is inherent to high-throughput approaches and specific to RNAi screens, false discovery due to off-target effects (OTEs) of RNAi reagents remains a problem. The need to improve our ability to use RNAi to elucidate gene function at large scale and in additional systems continues to be addressed through improved RNAi library design, development of innovative computational and analysis tools and other approaches.
The Drosophila MSL complex mediates dosage compensation by increasing transcription of the single X chromosome in males approximately two-fold. This is accomplished through recognition of the X chromosome and subsequent acetylation of histone H4K16 on X-linked genes. Initial binding to the X is thought to occur at "entry sites" that contain a consensus sequence motif ("MSL recognition element" or MRE). However, this motif is only ∼2 fold enriched on X, and only a fraction of the motifs on X are initially targeted. Here we ask whether chromatin context could distinguish between utilized and non-utilized copies of the motif, by comparing their relative enrichment for histone modifications and chromosomal proteins mapped in the modENCODE project. Through a comparative analysis of the chromatin features in male S2 cells (which contain MSL complex) and female Kc cells (which lack the complex), we find that the presence of active chromatin modifications, together with an elevated local GC content in the surrounding sequences, has strong predictive value for functional MSL entry sites, independent of MSL binding. We tested these sites for function in Kc cells by RNAi knockdown of Sxl, resulting in induction of MSL complex. We show that ectopic MSL expression in Kc cells leads to H4K16 acetylation around these sites and a relative increase in X chromosome transcription. Collectively, our results support a model in which a pre-existing active chromatin environment, coincident with H3K36me3, contributes to MSL entry site selection. The consequences of MSL targeting of the male X chromosome include increase in nucleosome lability, enrichment for H4K16 acetylation and JIL-1 kinase, and depletion of linker histone H1 on active X-linked genes. Our analysis can serve as a model for identifying chromatin and local sequence features that may contribute to selection of functional protein binding sites in the genome.
While the 26S proteasome is a key proteolytic complex, little is known about how proteasome levels are maintained in higher eukaryotic cells. Here we describe an RNA interference (RNAi) screen of Drosophila melanogaster that was used to identify transcription factors that may play a role in maintaining levels of the 26S proteasome. We used an RNAi library against 993 Drosophila transcription factor genes to identify genes whose suppression in Schneider 2 cells stabilized a ubiquitin-green fluorescent protein reporter protein. This screen identified Cnc (cap 'n' collar [CNC]; basic region leucine zipper) as a candidate transcriptional regulator of proteasome component expression. In fact, 20S proteasome activity was reduced in cells depleted of cnc. Immunoblot assays against proteasome components revealed a general decline in both 19S regulatory complex and 20S proteasome subunits after RNAi depletion of this transcription factor. Transcript-specific silencing revealed that the longest of the seven transcripts for the cnc gene, cnc-C, was needed for proteasome and p97 ATPase production. Quantitative reverse transcription-PCR confirmed the role of Cnc-C in activation of transcription of genes encoding proteasome components. Expression of a V5-His-tagged form of Cnc-C revealed that the transcription factor is itself a proteasome substrate that is stabilized when the proteasome is inhibited. We propose that this single cnc gene in Drosophila resembles the ancestral gene family of mammalian nuclear factor erythroid-derived 2-related transcription factors, which are essential in regulating oxidative stress and proteolysis.
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.
BACKGROUND: High-throughput screening using RNAi is a powerful gene discovery method but is often complicated by false positive and false negative results. Whereas false positive results associated with RNAi reagents has been a matter of extensive study, the issue of false negatives has received less attention. RESULTS: We performed a meta-analysis of several genome-wide, cell-based Drosophila RNAi screens, together with a more focused RNAi screen, and conclude that the rate of false negative results is at least 8%. Further, we demonstrate how knowledge of the cell transcriptome can be used to resolve ambiguous results and how the number of false negative results can be reduced by using multiple, independently-tested RNAi reagents per gene. CONCLUSIONS: RNAi reagents that target the same gene do not always yield consistent results due to false positives and weak or ineffective reagents. False positive results can be partially minimized by filtering with transcriptome data. RNAi libraries with multiple reagents per gene also reduce false positive and false negative outcomes when inconsistent results are disambiguated carefully.
Wnt proteins are secreted, lipid-modified glycoproteins that control animal development and adult tissue homeostasis. Secretion of Wnt proteins is at least partly regulated by a dedicated machinery. Here, we report a genome-wide RNA interference screen for genes involved in the secretion of Wingless (Wg), a Drosophila Wnt. We identify three new genes required for Wg secretion. Of these, Emp24 and Eclair are required for proper export of Wg from the endoplasmic reticulum (ER). We propose that Emp24 and Eca act as specific cargo receptors for Wg to concentrate it in forming vesicles at sites of ER export.
The DNA damage checkpoint, the first pathway known to be activated in response to DNA damage, is a mechanism by which the cell cycle is temporarily arrested to allow DNA repair. The checkpoint pathway transmits signals from the sites of DNA damage to the cell cycle machinery through the evolutionarily conserved ATM (ataxia telangiectasia mutated) and ATR (ATM- and Rad3-related) kinase cascades. We conducted a genome-wide RNAi (RNA interference) screen in Drosophila cells to identify previously unknown genes and pathways required for the G₂-M checkpoint induced by DNA double-strand breaks (DSBs). Our large-scale analysis provided a systems-level view of the G₂-M checkpoint and revealed the coordinated actions of particular classes of proteins, which include those involved in DNA repair, DNA replication, cell cycle control, chromatin regulation, and RNA processing. Further, from the screen and in vivo analysis, we identified previously unrecognized roles of two DNA damage response genes, mus101 and mus312. Our results suggest that the DNA replication preinitiation complex, which includes MUS101, and the MUS312-containing nuclease complexes, which are important for DSB repair, also function in the G₂-M checkpoint. Our results provide insight into the diverse mechanisms that link DNA damage and the checkpoint signaling pathway.
To identify Huntington's Disease therapeutics, we conducted high-content small molecule and RNAi suppressor screens using a Drosophila primary neural culture Huntingtin model. Drosophila primary neurons offer a sensitive readout for neurotoxicty, as their neurites develop dysmorphic features in the presence of mutant polyglutamine-expanded Huntingtin compared to nonpathogenic Huntingtin. By tracking the subcellular distribution of mRFP-tagged pathogenic Huntingtin and assaying neurite branch morphology via live-imaging, we identified suppressors that could reduce Huntingtin aggregation and/or prevent the formation of dystrophic neurites. The custom algorithms we used to quantify neurite morphologies in complex cultures provide a useful tool for future high-content screening approaches focused on neurodegenerative disease models. Compounds previously found to be effective aggregation inhibitors in mammalian systems were also effective in Drosophila primary cultures, suggesting translational capacity between these models. However, we did not observe a direct correlation between the ability of a compound or gene knockdown to suppress aggregate formation and its ability to rescue dysmorphic neurites. Only a subset of aggregation inhibitors could revert dysmorphic cellular profiles. We identified lkb1, an upstream kinase in the mTOR/Insulin pathway, and four novel drugs, Camptothecin, OH-Camptothecin, 18β-Glycyrrhetinic acid, and Carbenoxolone, that were strong suppressors of mutant Huntingtin-induced neurotoxicity. Huntingtin neurotoxicity suppressors identified through our screen also restored viability in an in vivo Drosophila Huntington's Disease model, making them attractive candidates for further therapeutic evaluation.