GUItars: a GUI tool for analysis of high-throughput RNA interference screening data

Genome-wide screens identify specific drivers of mutant hTERT promoters

Mutations in hTERT promoter are seen in over 19% of human cancers, irrespective of the cancer type. Understanding the molecular players that regulate Mut-hTERT promoters may help the design of effective targeting strategies to inhibit telomerase reactivation specifically in cancer cells. Our work uses genome-wide functional screens to identify 30 specific regulators of Mut-hTERT promoters. These candidates identified from the screening serve as an excellent resource to understand how telomerase is reactivated and as targets for making inhibitors to telomerase, a key driver of cancer.

Cancer-specific hTERT promoter mutations reported in 19% of cancers result in enhanced telomerase activity. Understanding the distinctions between transcriptional regulation of wild-type (WT) and mutant (Mut) hTERT promoters may open up avenues for development of inhibitors which specially block hTERT expression in cancer cells. To comprehensively identify physiological regulators of WT- or Mut-hTERT promoters, we generated several isogenic reporter cells driven by endogenous hTERT loci. Genome-wide CRISPR-Cas9 and small interfering RNA screens using these isogenic reporter lines identified specific regulators of Mut-hTERT promoters. We validate and characterize one of these hits, namely, MED12, a kinase subunit of mediator complex. We demonstrate that MED12 specifically drives expression of hTERT from the Mut-hTERT promoter by mediating long-range chromatin interaction between the proximal Mut-hTERT promoter and T-INT1 distal regulatory region 260 kb upstream. Several hits identified in our screens could serve as potential therapeutic targets, inhibition of which may specifically block Mut-hTERT promoter driven telomerase reactivation in cancers.

OpenStats: A robust and scalable software package for reproducible analysis of high-throughput phenotypic data

Reproducibility in the statistical analyses of data from high-throughput phenotyping screens requires a robust and reliable analysis foundation that allows modelling of different possible statistical scenarios. Regular challenges are scalability and extensibility of the analysis software. In this manuscript, we describe OpenStats, a freely available software package that addresses these challenges. We show the performance of the software in a high-throughput phenomic pipeline in the International Mouse Phenotyping Consortium (IMPC) and compare the agreement of the results with the most similar implementation in the literature. OpenStats has significant improvements in speed and scalability compared to existing software packages including a 13-fold improvement in computational time to the current production analysis pipeline in the IMPC. Reduced complexity also promotes FAIR data analysis by providing transparency and benefiting other groups in reproducing and re-usability of the statistical methods and results. OpenStats is freely available under a Creative Commons license at www.bioconductor.org/packages/OpenStats.

Issues Of Z-factor and an approach to avoid them for quality control in high-throughput screening studies

MOTIVATION

High throughput screening (HTS) is a vital automation technology in biomedical research in both industry and academia. The well-known z-factor has been widely used as a gatekeeper to assure assay quality in an HTS study. However, many researchers and users may not have realized that z-factor has major issues.

RESULTS

In this article, the following four major issues are explored and demonstrated so that researchers may use the z-factor appropriately. First, the z-factor violates the Pythagorean Theorem of Statistics. Second, there is no adjustment of sampling error in the application of the z-factor for quality control (QC) in HTS studies. Third, the expectation of the sample-based z-factor does not exist. Fourth, the thresholds in the z-factor based criterion lack a theoretical basis. Here, an approach to avoid these issues was proposed and new QC criteria under homoscedasticity were constructed so that researchers can choose a statistically grounded criterion for QC in the HTS studies. We implemented this approach in an R package and demonstrated its utility in multiple CRISPR/CAS9 or siRNA HTS studies.

AVAILABILITY

The R package qcSSMDhomo is freely available from GitHub: https://github.com/Karena6688/qcSSMDhomo. The file qcSSMDhomo_1.0.0.tar.gz (for Windows) containing qcSSMDhomo is also available at Bioinformatics online. qcSSMDhomo is distributed under the GNU General Public License.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

G3BP1 Is a Tunable Switch that Triggers Phase Separation to Assemble Stress Granules

The mechanisms underlying ribonucleoprotein (RNP) granule assembly, including the basis for establishing and maintaining RNP granules with distinct composition, are unknown. One prominent type of RNP granule is the stress granule (SG), a dynamic and reversible cytoplasmic assembly formed in eukaryotic cells in response to stress. Here, we show that SGs assemble through liquid-liquid phase separation (LLPS) arising from interactions distributed unevenly across a core protein-RNA interaction network. The central node of this network is G3BP1, which functions as a molecular switch that triggers RNA-dependent LLPS in response to a rise in intracellular free RNA concentrations. Moreover, we show that interplay between three distinct intrinsically disordered regions (IDRs) in G3BP1 regulates its intrinsic propensity for LLPS, and this is fine-tuned by phosphorylation within the IDRs. Further regulation of SG assembly arises through positive or negative cooperativity by extrinsic G3BP1-binding factors that strengthen or weaken, respectively, the core SG network.

Discovery of Druggable Host Factors Critical to Plasmodium Liver-Stage Infection

Plasmodium parasites undergo an obligatory and asymptomatic developmental stage within the liver before infecting red blood cells to cause malaria. The hijacked host pathways critical to parasite infection during this hepatic phase remain poorly understood. Here, we implemented a forward genetic screen to identify over 100 host factors within the human druggable genome that are critical to P. berghei infection in hepatoma cells. Notably, we found knockdown of genes involved in protein trafficking pathways to be detrimental to parasite infection. The disruption of protein trafficking modulators, including COPB2 and GGA1, decreases P. berghei parasite size, and an immunofluorescence study suggests that these proteins are recruited to the Plasmodium parasitophorous vacuole in infected hepatocytes. These findings reveal that various host intracellular protein trafficking pathways are subverted by Plasmodium parasites during the liver stage and provide new insights into their manipulation for growth and development.

HTS-Net: An integrated regulome-interactome approach for establishing network regulation models in high-throughput screenings

High-throughput RNAi screenings (HTS) allow quantifying the impact of the deletion of each gene in any particular function, from virus-host interactions to cell differentiation. However, there has been less development for functional analysis tools dedicated to RNAi analyses. HTS-Net, a network-based analysis program, was developed to identify gene regulatory modules impacted in high-throughput screenings, by integrating transcription factors-target genes interaction data (regulome) and protein-protein interaction networks (interactome) on top of screening z-scores. HTS-Net produces exhaustive HTML reports for results navigation and exploration. HTS-Net is a new pipeline for RNA interference screening analyses that proves better performance than simple gene rankings by z-scores, by re-prioritizing genes and replacing them in their biological context, as shown by the three studies that we reanalyzed. Formatted input data for the three studied datasets, source code and web site for testing the system are available from the companion web site at http://htsnet.marseille.inserm.fr/. We also compared our program with existing algorithms (CARD and hotnet2).

Transcription factor ZNF148 is a negative regulator of human muscle differentiation

Muscle differentiation is a complex process in which muscle progenitor cells undergo determination and eventually cellular fusion. This process is heavily regulated by such master transcription factors as MYOD and members of the MEF2 family. Here, we show that the transcription factor ZNF148 plays a direct role in human muscle cell differentiation. Downregulation of ZNF148 drives the formation of a muscle phenotype with rapid expression of myosin heavy chain, even in proliferative conditions. This phenotype was most likely mediated by the robust and swift upregulation of MYOD and MEF2C.

Data Analysis for High-Throughput RNAi Screening

High-throughput RNA interference (HT-RNAi) screening is an effective technology to help identify important genes and pathways involved in a biological process. Analysis of high-throughput RNAi screening data is a critical part of this technology, and many analysis methods have been described. Here, we summarize the workflow and types of analyses commonly used in high-throughput RNAi screening.

A Genome-Wide Screen for Machinery Involved in Downregulation of MHC Class I by HIV-1 Nef

The HIV-1-encoded protein, Nef, plays a key role in the development of AIDS. One of Nef’s functions is to keep MHC class I off the surface of infected cells, a process that requires the host proteins clathrin and AP-1. To identify other proteins involved in this pathway, we carried out a genome-wide siRNA library screen on HeLa cells co-expressing HLA-A2 and an inducible form of Nef. Out of 21,121 siRNA pools, 100 were selected for further analysis, based on their ability to either inhibit or enhance downregulation of MHC-I by Nef. When cells were treated with the same siRNA pools as those used in the screen, 79% produced a similar phenotype. However, when the cells were treated with different siRNA reagents targeting the same genes, only 16% produced a similar phenotype. This indicates that most of the hits found in the original screen are likely to have been off-target, an important concern that is often not taken into account in siRNA screening studies. Nevertheless, we identified novel host factors involved in Nef-induced downregulation of MHC-I, including four genes, MIIP, CAMSAP3, SLC6A3, and KCTD19, where multiple reagents produced a strong inhibitory effect on Nef activity. Other hits slightly below our very high stringency cutoff point may also deserve further study. Thus, our dataset is a valuable resource for scientists investigating the pathogenesis of HIV.

PI4K-beta and MKNK1 are regulators of hepatitis C virus IRES-dependent translation

Cellular translation is down-regulated by host antiviral responses. Picornaviridae and Flaviviridae including hepatitis C virus (HCV) evade this process using internal ribosomal entry sequences (IRESs). Although HCV IRES translation is a prerequisite for HCV replication, only few host factors critical for IRES activity are known and the global regulator network remains largely unknown. Since signal transduction is an import regulator of viral infections and the host antiviral response we combined a functional RNAi screen targeting the human signaling network with a HCV IRES-specific reporter mRNA assay. We demonstrate that the HCV host cell cofactors PI4K and MKNK1 are positive regulators of HCV IRES translation representing a novel pathway with a functional relevance for the HCV life cycle and IRES-mediated translation of viral RNA.

Viability of developmental stages of Schistosoma mansoni quantified with xCELLigence worm real-time motility assay (xWORM)

Infection with helminth parasites causes morbidity and mortality in billions of people and livestock worldwide. Where anthelmintic drugs are available, drug resistance is a major problem in livestock parasites, and a looming threat to public health. Monitoring the efficacy of these medicines and screening for new drugs has been hindered by the lack of objective, high-throughput approaches. Several cell monitoring technologies have been adapted for parasitic worms, including video-, fluorescence-, metabolism enzyme- and impedance-based tools that minimize the screening bottleneck. Using the xCELLigence impedance-based system we previously developed a motility-viability assay that is applicable for a range of helminth parasites. Here we have improved substantially the assay by using diverse frequency settings, and have named it the xCELLigence worm real-time motility assay (xWORM). By utilizing strictly standardized mean difference analysis we compared the xWORM output measured with 10, 25 and 50 kHz frequencies to quantify the motility of schistosome adults (human blood flukes) and hatching of schistosome eggs. Furthermore, we have described a novel application of xWORM to monitor movement of schistosome cercariae, the developmental stage that is infectious to humans. For all three stages, 25 kHz was either optimal or near-optimal for monitoring and quantifying schistosome motility. These improvements in methodology sensitivity should enhance the capacity to screen small compound libraries for new drugs both for schistosomes and other helminth pathogens at large.

•

25 kHz on the xCELLigence system dramatically improves the schistosome xWORM assay.

•

xWORM assay can efficiently determine viability of Schistome adults or eggs.

•

First time cercariae have been incorporated into an automated viability assay.

•

Other helminth monitoring may benefit from alternate xCELLigence frequency options.

Identification of Selective Inhibitors of the Plasmodium falciparum Hexose Transporter PfHT by Screening Focused Libraries of Anti-Malarial Compounds

Development of resistance against current antimalarial drugs necessitates the search for novel drugs that interact with different targets and have distinct mechanisms of action. Malaria parasites depend upon high levels of glucose uptake followed by inefficient metabolic utilization via the glycolytic pathway, and the Plasmodium falciparum hexose transporter PfHT, which mediates uptake of glucose, has thus been recognized as a promising drug target. This transporter is highly divergent from mammalian hexose transporters, and it appears to be a permease that is essential for parasite viability in intra-erythrocytic, mosquito, and liver stages of the parasite life cycle. An assay was developed that is appropriate for high throughput screening against PfHT based upon heterologous expression of PfHT in Leishmania mexicana parasites that are null mutants for their endogenous hexose transporters. Screening of two focused libraries of antimalarial compounds identified two such compounds that are high potency selective inhibitors of PfHT compared to human GLUT1. Additionally, 7 other compounds were identified that are lower potency and lower specificity PfHT inhibitors but might nonetheless serve as starting points for identification of analogs with more selective properties. These results further support the potential of PfHT as a novel drug target.

Regulators of autophagosome formation in Drosophila muscles

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.

Since the identification of the core autophagy genes in yeast, tissue culture cell lines have been the primary tool to evaluate the role and regulation of autophagy in higher organisms. However, since autophagy is a tissue-specific, context dependent process, stable cell lines can only give a limited view of the autophagic process. Here, we focus on the role of putative autophagy regulators in the specific context of the skeletal muscle, which has one of the most robust autophagy responses in mammals. We describe a fruitfly model of autophagy for skeletal muscles that integrates rapid genetic screening in primary cultured cells with robust in vivo validation in the larval muscle. We screened a set of genes previously linked to the autophagy pathway in humans, and identified both positive and negative regulators of autophagy. Our observation that genes involved in sugar metabolism impact muscle autophagy has important implications for both skeletal and cardiac myopathies associated with aberrant sugar storage. Surprisingly, we found that the proteasome is required to maintain autophagy in the muscle, suggesting that therapeutic treatments aiming to induce autophagy by proteasome inhibition must be carefully calibrated to ensure that the opposite phenotype does not occur.

Stability of the human pregnane X receptor is regulated by E3 ligase UBR5 and serine/threonine kinase DYRK2

The hPXR (human pregnane X receptor), a major chemical toxin sensor, is a ligand-induced transcription factor activated by various xenobiotics and toxins, resulting in the transcriptional up-regulation of detoxifying enzymes. To date, little is known about the upstream regulation of hPXR. Using MS analysis and a kinome-wide siRNA screen, we report that the E3 ligase UBR5 (ubiquitin protein ligase E3 component n-recognin 5) and DYRK2 (dual-specificity tyrosine-phosphorylation-regulated kinase 2) regulate hPXR stability. UBR5 knockdown resulted in accumulation of cellular hPXR and a concomitant increase in hPXR activity, whereas the rescue of UBR5 knockdown decreased the cellular hPXR level and activity. Importantly, UBR5 exerted its effect in concert with the serine/threonine kinase DYRK2, as the knockdown of DYRK2 phenocopied UBR5 knockdown. hPXR was shown to be a substrate for DYRK2, and DYRK2-dependent phosphorylation of hPXR facilitated its subsequent ubiquitination by UBR5. This is the first report of the post-translational regulation of hPXR via phosphorylation-facilitated ubiquitination by DYRK2 and UBR5. The results of the present study reveal the role of the ubiquitin-proteasomal pathway in modulating hPXR activity and indicate that pharmacological inhibitors of the ubiquitin-proteasomal pathway that regulate hPXR stability may negatively affect treatment outcome from unintended hPXR-mediated drug-drug interactions.