Cell-based high content screening (HCS) is becoming an important and increasingly favored approach in therapeutic drug discovery and functional genomics. In HCS, changes in cellular morphology and biomarker distributions provide an information-rich profile of cellular responses to experimental treatments such as small molecules or gene knockdown probes. One obstacle that currently exists with such cell-based assays is the availability of image processing algorithms that are capable of reliably and automatically analyzing large HCS image sets. HCS images of primary neuronal cell cultures are particularly challenging to analyze due to complex cellular morphology. Here we present a robust method for quantifying and statistically analyzing the morphology of neuronal cells in HCS images. The major advantages of our method over existing software lie in its capability to correct non-uniform illumination using the contrast-limited adaptive histogram equalization method; segment neuromeres using Gabor-wavelet texture analysis; and detect faint neurites by a novel phase-based neurite extraction algorithm that is invariant to changes in illumination and contrast and can accurately localize neurites. Our method was successfully applied to analyze a large HCS image set generated in a morphology screen for polyglutamine-mediated neuronal toxicity using primary neuronal cell cultures derived from embryos of a Drosophila Huntington's Disease (HD) model.

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.

Gene silencing through sequence-specific targeting of mRNAs by RNAi has enabled genome-wide functional screens in cultured cells and in vivo in model organisms. These screens have resulted in the identification of new cellular pathways and potential drug targets. Considerable progress has been made to improve the quality of RNAi screen data through the development of new experimental and bioinformatics approaches. The recent availability of genome-editing strategies, such as the CRISPR (clustered regularly interspaced short palindromic repeats)-Cas9 system, when combined with RNAi, could lead to further improvements in screen data quality and follow-up experiments, thus promoting our understanding of gene function and gene regulatory networks.

Certain pathogens, such as Mycobacterium tuberculosis, survive within the hostile intracellular environment of a macrophage. To identify host factors required for mycobacterial entry and survival within macrophages, we performed a genomewide RNA interference screen in Drosophila macrophage-like cells, using Mycobacterium fortuitum. We identified factors required for general phagocytosis, as well as those needed specifically for mycobacterial infection. One specific factor, Peste (Pes), is a CD36 family member required for uptake of mycobacteria, but not Escherichia coli or Staphylococcus aureus. Moreover, mammalian class B scavenger receptors (SRs) conferred uptake of bacteria into nonphagocytic cells, with SR-BI and SR-BII uniquely mediating uptake of M. fortuitum, which suggests a conserved role for class B SRs in pattern recognition and innate immunity.

Although classical genetic and biochemical approaches have identified hundreds of proteins that function in the dynamic remodeling of cell shape in response to upstream signals, there is currently little systems-level understanding of the organization and composition of signaling networks that regulate cell morphology. We have developed quantitative morphological profiling methods to systematically investigate the role of individual genes in the regulation of cell morphology in a fast, robust, and cost-efficient manner. We analyzed a compendium of quantitative morphological signatures and described the existence of local signaling networks that act to regulate cell protrusion, adhesion, and tension.

Cellular signaling networks have evolved to enable swift and accurate responses, even in the face of genetic or environmental perturbation. Thus, genetic screens may not identify all the genes that regulate different biological processes. Moreover, although classical screening approaches have succeeded in providing parts lists of the essential components of signaling networks, they typically do not provide much insight into the hierarchical and functional relations that exist among these components. We describe a high-throughput screen in which we used RNA interference to systematically inhibit two genes simultaneously in 17,724 combinations to identify regulators of Drosophila JUN NH(2)-terminal kinase (JNK). Using both genetic and phosphoproteomics data, we then implemented an integrative network algorithm to construct a JNK phosphorylation network, which provides structural and mechanistic insights into the systems architecture of JNK signaling.

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.

A large fraction of our genome consists of mobile genetic elements. Governing transposons in germ cells is critically important, and failure to do so compromises genome integrity, leading to sterility. In animals, the piRNA pathway is the key to transposon constraint, yet the precise molecular details of how piRNAs are formed and how the pathway represses mobile elements remain poorly understood. In an effort to identify general requirements for transposon control and components of the piRNA pathway, we carried out a genome-wide RNAi screen in Drosophila ovarian somatic sheet cells. We identified and validated 87 genes necessary for transposon silencing. Among these were several piRNA biogenesis factors. We also found CG3893 (asterix) to be essential for transposon silencing, most likely by contributing to the effector step of transcriptional repression. Asterix loss leads to decreases in H3K9me3 marks on certain transposons but has no effect on piRNA levels.

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.

Recent studies by our group and others demonstrated a required and conserved role of Stim in store-operated Ca(2+) influx and Ca(2+) release-activated Ca(2+) (CRAC) channel activity. By using an unbiased genome-wide RNA interference screen in Drosophila S2 cells, we now identify 75 hits that strongly inhibited Ca(2+) influx upon store emptying by thapsigargin. Among these hits are 11 predicted transmembrane proteins, including Stim, and one, olf186-F, that upon RNA interference-mediated knockdown exhibited a profound reduction of thapsigargin-evoked Ca(2+) entry and CRAC current, and upon overexpression a 3-fold augmentation of CRAC current. CRAC currents were further increased to 8-fold higher than control and developed more rapidly when olf186-F was cotransfected with Stim. olf186-F is a member of a highly conserved family of four-transmembrane spanning proteins with homologs from Caenorhabditis elegans to human. The endoplasmic reticulum (ER) Ca(2+) pump sarco-/ER calcium ATPase (SERCA) and the single transmembrane-soluble N-ethylmaleimide-sensitive (NSF) attachment receptor (SNARE) protein Syntaxin5 also were required for CRAC channel activity, consistent with a signaling pathway in which Stim senses Ca(2+) depletion within the ER, translocates to the plasma membrane, and interacts with olf186-F to trigger CRAC channel activity.

Nearly 1.7 billion people are infected with Mycobacterium tuberculosis. Its ability to survive intracellularly is thought to be central to its success as a pathogen, but how it does this is poorly understood. Using a Drosophila model of infection, we identify three host cell activities, Rab7, CG8743, and the ESCRT machinery, that modulate the mycobacterial phagosome. In the absence of these factors the cell no longer restricts growth of the non-pathogen Mycobacterium smegmatis. Hence, we identify factors that represent unique vulnerabilities of the host cell, because manipulation of any one of them alone is sufficient to allow a nonpathogenic mycobacterial species to proliferate. Furthermore, we demonstrate that, in mammalian cells, the ESCRT machinery plays a conserved role in restricting bacterial growth.

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.

Akt represents a nodal point between the Insulin receptor and TOR signaling, and its activation by phosphorylation controls cell proliferation, cell size, and metabolism. The activity of Akt must be carefully balanced, as increased Akt signaling is frequently associated with cancer and as insufficient Akt signaling is linked to metabolic disease and diabetes mellitus. Using a genome-wide RNAi screen in Drosophila cells in culture, and in vivo analyses in the third instar wing imaginal disc, we studied the regulatory circuitries that define dAkt activation. We provide evidence that negative feedback regulation of dAkt occurs during normal Drosophila development in vivo. Whereas in cell culture dAkt is regulated by S6 Kinase (S6K)-dependent negative feedback, this feedback inhibition only plays a minor role in vivo. In contrast, dAkt activation under wild-type conditions is defined by feedback inhibition that depends on TOR Complex 1 (TORC1), but is S6K-independent. This feedback inhibition is switched from TORC1 to S6K only in the context of enhanced TORC1 activity, as triggered by mutations in tsc2. These results illustrate how the Akt-TOR pathway dynamically adapts the routing of negative feedback in response to the activity load of its signaling circuit in vivo.

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.

Given the diversity of autophagy targets and regulation, it is important to characterize autophagy in various cell types and conditions. We used a primary myocyte cell culture system to assay the role of putative autophagy regulators in the specific context of skeletal muscle. By treating the cultures with rapamycin (Rap) and chloroquine (CQ) we induced an autophagic response, fully suppressible by knockdown of core ATG genes. We screened D. melanogaster orthologs of a previously reported mammalian autophagy protein-protein interaction network, identifying several proteins required for autophagosome formation in muscle cells, including orthologs of the Rab regulators RabGap1 and Rab3Gap1. The screen also highlighted the critical roles of the proteasome and glycogen metabolism in regulating autophagy. Specifically, sustained proteasome inhibition inhibited autophagosome formation both in primary culture and larval skeletal muscle, even though autophagy normally acts to suppress ubiquitin aggregate formation in these tissues. In addition, analyses of glycogen metabolic genes in both primary cultured and larval muscles indicated that glycogen storage enhances the autophagic response to starvation, an important insight given the link between glycogen storage disorders, autophagy, and muscle function.