Growth-restricting effects of siRNA transfections: a largely deterministic combination of off-target binding and hybridization-independent competition.

Perturbation of gene expression by means of synthetic small interfering RNAs (siRNAs) is a powerful way to uncover gene function. However, siRNA technology suffers from sequence-specific off-target effects and from limitations in knock-down efficiency. In this study, we assess a further problem: unintended effects of siRNA transfections on cellular fitness/proliferation. We show that the nucleotide compositions of siRNAs at specific positions have reproducible growth-restricting effects on mammalian cells in culture. This is likely distinct from hybridization-dependent off-target effects, since each nucleotide residue is seen to be acting independently and additively. The effect is robust and reproducible across different siRNA libraries and also across various cell lines, including human and mouse cells. Analyzing the growth inhibition patterns in correlation to the nucleotide sequence of the siRNAs allowed us to build a predictor that can estimate growth-restricting effects for any arbitrary siRNA sequence. Competition experiments with co-transfected siRNAs further suggest that the growth-restricting effects might be linked to an oversaturation of the cellular miRNA machinery, thus disrupting endogenous miRNA functions at large. We caution that competition between siRNA molecules could complicate the interpretation of double-knockdown or epistasis experiments, and potential interactions with endogenous miRNAs can be a factor when assaying cell growth or viability phenotypes.

PyBDA: a command line tool for automated analysis of big biological data sets.

BACKGROUND: Analysing large and high-dimensional biological data sets poses significant computational difficulties for bioinformaticians due to lack of accessible tools that scale to hundreds of millions of data points. RESULTS: We developed a novel machine learning command line tool called PyBDA for automated, distributed analysis of big biological data sets. By using Apache Spark in the backend, PyBDA scales to data sets beyond the size of current applications. It uses Snakemake in order to automatically schedule jobs to a high-performance computing cluster. We demonstrate the utility of the software by analyzing image-based RNA interference data of 150 million single cells. CONCLUSION: PyBDA allows automated, easy-to-use data analysis using common statistical methods and machine learning algorithms. It can be used with simple command line calls entirely making it accessible to a broad user base. PyBDA is available at https://pybda.rtfd.io.

Improved pathway reconstruction from RNA interference screens by exploiting off-target effects.

Motivation: Pathway reconstruction has proven to be an indispensable tool for analyzing the molecular mechanisms of signal transduction underlying cell function. Nested effects models (NEMs) are a class of probabilistic graphical models designed to reconstruct signalling pathways from high-dimensional observations resulting from perturbation experiments, such as RNA interference (RNAi). NEMs assume that the short interfering RNAs (siRNAs) designed to knockdown specific genes are always on-target. However, it has been shown that most siRNAs exhibit strong off-target effects, which further confound the data, resulting in unreliable reconstruction of networks by NEMs. Results: Here, we present an extension of NEMs called probabilistic combinatorial nested effects models (pc-NEMs), which capitalize on the ancillary siRNA off-target effects for network reconstruction from combinatorial gene knockdown data. Our model employs an adaptive simulated annealing search algorithm for simultaneous inference of network structure and error rates inherent to the data. Evaluation of pc-NEMs on simulated data with varying number of phenotypic effects and noise levels as well as real data demonstrates improved reconstruction compared to classical NEMs. Application to Bartonella henselae infection RNAi screening data yielded an eight node network largely in agreement with previous works, and revealed novel binary interactions of direct impact between established components. Availability and implementation: The software used for the analysis is freely available as an R package at https://github.com/cbg-ethz/pcNEM.git. Supplementary information: Supplementary data are available at Bioinformatics online.

ALPK1 controls TIFA/TRAF6-dependent innate immunity against heptose-1,7-bisphosphate of gram-negative bacteria.

During infection by invasive bacteria, epithelial cells contribute to innate immunity via the local secretion of inflammatory cytokines. These are directly produced by infected cells or by uninfected bystanders via connexin-dependent cell-cell communication. However, the cellular pathways underlying this process remain largely unknown. Here we perform a genome-wide RNA interference screen and identify TIFA and TRAF6 as central players of Shigella flexneri and Salmonella typhimurium-induced interleukin-8 expression. We show that threonine 9 and the forkhead-associated domain of TIFA are necessary for the oligomerization of TIFA in both infected and bystander cells. Subsequently, this process triggers TRAF6 oligomerization and NF-κB activation. We demonstrate that TIFA/TRAF6-dependent cytokine expression is induced by the bacterial metabolite heptose-1,7-bisphosphate (HBP). In addition, we identify alpha-kinase 1 (ALPK1) as the critical kinase responsible for TIFA oligomerization and IL-8 expression in response to infection with S. flexneri and S. typhimurium but also to Neisseria meningitidis. Altogether, these results clearly show that ALPK1 is a master regulator of innate immunity against both invasive and extracellular gram-negative bacteria.

Specific inhibition of diverse pathogens in human cells by synthetic microRNA-like oligonucleotides inferred from RNAi screens.

Systematic genetic perturbation screening in human cells remains technically challenging. Typically, large libraries of chemically synthesized siRNA oligonucleotides are used, each designed to degrade a specific cellular mRNA via the RNA interference (RNAi) mechanism. Here, we report on data from three genome-wide siRNA screens, conducted to uncover host factors required for infection of human cells by two bacterial and one viral pathogen. We find that the majority of phenotypic effects of siRNAs are unrelated to the intended "on-target" mechanism, defined by full complementarity of the 21-nt siRNA sequence to a target mRNA. Instead, phenotypes are largely dictated by "off-target" effects resulting from partial complementarity of siRNAs to multiple mRNAs via the "seed" region (i.e., nucleotides 2-8), reminiscent of the way specificity is determined for endogenous microRNAs. Quantitative analysis enabled the prediction of seeds that strongly and specifically block infection, independent of the intended on-target effect. This prediction was confirmed experimentally by designing oligos that do not have any on-target sequence match at all, yet can strongly reproduce the predicted phenotypes. Our results suggest that published RNAi screens have primarily, and unintentionally, screened the sequence space of microRNA seeds instead of the intended on-target space of protein-coding genes. This helps to explain why previously published RNAi screens have exhibited relatively little overlap. Our analysis suggests a possible way of identifying "seed reagents" for controlling phenotypes of interest and establishes a general strategy for extracting valuable untapped information from past and future RNAi screens.

NEMix: single-cell nested effects models for probabilistic pathway stimulation.

Nested effects models have been used successfully for learning subcellular networks from high-dimensional perturbation effects that result from RNA interference (RNAi) experiments. Here, we further develop the basic nested effects model using high-content single-cell imaging data from RNAi screens of cultured cells infected with human rhinovirus. RNAi screens with single-cell readouts are becoming increasingly common, and they often reveal high cell-to-cell variation. As a consequence of this cellular heterogeneity, knock-downs result in variable effects among cells and lead to weak average phenotypes on the cell population level. To address this confounding factor in network inference, we explicitly model the stimulation status of a signaling pathway in individual cells. We extend the framework of nested effects models to probabilistic combinatorial knock-downs and propose NEMix, a nested effects mixture model that accounts for unobserved pathway activation. We analyzed the identifiability of NEMix and developed a parameter inference scheme based on the Expectation Maximization algorithm. In an extensive simulation study, we show that NEMix improves learning of pathway structures over classical NEMs significantly in the presence of hidden pathway stimulation. We applied our model to single-cell imaging data from RNAi screens monitoring human rhinovirus infection, where limited infection efficiency of the assay results in uncertain pathway stimulation. Using a subset of genes with known interactions, we show that the inferred NEMix network has high accuracy and outperforms the classical nested effects model without hidden pathway activity. NEMix is implemented as part of the R/Bioconductor package 'nem' and available at www.cbg.ethz.ch/software/NEMix.

Simultaneous analysis of large-scale RNAi screens for pathogen entry.

BACKGROUND: Large-scale RNAi screening has become an important technology for identifying genes involved in biological processes of interest. However, the quality of large-scale RNAi screening is often deteriorated by off-targets effects. In order to find statistically significant effector genes for pathogen entry, we systematically analyzed entry pathways in human host cells for eight pathogens using image-based kinome-wide siRNA screens with siRNAs from three vendors. We propose a Parallel Mixed Model (PMM) approach that simultaneously analyzes several non-identical screens performed with the same RNAi libraries. RESULTS: We show that PMM gains statistical power for hit detection due to parallel screening. PMM allows incorporating siRNA weights that can be assigned according to available information on RNAi quality. Moreover, PMM is able to estimate a sharedness score that can be used to focus follow-up efforts on generic or specific gene regulators. By fitting a PMM model to our data, we found several novel hit genes for most of the pathogens studied. CONCLUSIONS: Our results show parallel RNAi screening can improve the results of individual screens. This is currently particularly interesting when large-scale parallel datasets are becoming more and more publicly available. Our comprehensive siRNA dataset provides a public, freely available resource for further statistical and biological analyses in the high-content, high-throughput siRNA screening field.

Multi-scale Gaussian representation and outline-learning based cell image segmentation.

BACKGROUND: High-throughput genome-wide screening to study gene-specific functions, e.g. for drug discovery, demands fast automated image analysis methods to assist in unraveling the full potential of such studies. Image segmentation is typically at the forefront of such analysis as the performance of the subsequent steps, for example, cell classification, cell tracking etc., often relies on the results of segmentation. METHODS: We present a cell cytoplasm segmentation framework which first separates cell cytoplasm from image background using novel approach of image enhancement and coefficient of variation of multi-scale Gaussian scale-space representation. A novel outline-learning based classification method is developed using regularized logistic regression with embedded feature selection which classifies image pixels as outline/non-outline to give cytoplasm outlines. Refinement of the detected outlines to separate cells from each other is performed in a post-processing step where the nuclei segmentation is used as contextual information. RESULTS AND CONCLUSIONS: We evaluate the proposed segmentation methodology using two challenging test cases, presenting images with completely different characteristics, with cells of varying size, shape, texture and degrees of overlap. The feature selection and classification framework for outline detection produces very simple sparse models which use only a small subset of the large, generic feature set, that is, only 7 and 5 features for the two cases. Quantitative comparison of the results for the two test cases against state-of-the-art methods show that our methodology outperforms them with an increase of 4-9% in segmentation accuracy with maximum accuracy of 93%. Finally, the results obtained for diverse datasets demonstrate that our framework not only produces accurate segmentation but also generalizes well to different segmentation tasks.

The E3 Ubiquitin Ligase Mind Bomb 1 Controls Adenovirus Genome Release at the Nuclear Pore Complex.

Adenoviruses (AdVs) cause respiratory, ocular, and gastrointestinal tract infection and inflammation in immunocompetent people and life-threatening disease upon immunosuppression. AdV vectors are widely used in gene therapy and vaccination. Incoming particles attach to nuclear pore complexes (NPCs) of post-mitotic cells, then rupture and deliver viral DNA (vDNA) to the nucleus or misdeliver to the cytosol. Our genome-wide RNAi screen in AdV-infected cells identified the RING-type E3 ubiquitin ligase Mind bomb 1 (Mib1) as a proviral host factor for AdV infection. Mib1 is implicated in Notch-Delta signaling, ciliary biogenesis, and RNA innate immunity. Mib1 depletion arrested incoming AdVs at NPCs. Induced expression of full-length but not ligase-defective Mib1 in knockout cells triggered vDNA uncoating from NPC-tethered virions, nuclear import, misdelivery of vDNA, and vDNA expression. Mib1 is an essential host factor for AdV uncoating in human cells, and it provides a new concept for licensing virion DNA delivery through the NPC.

A Role for the VPS Retromer in Brucella Intracellular Replication Revealed by Genomewide siRNA Screening.

Brucella, the agent causing brucellosis, is a major zoonotic pathogen with worldwide distribution. Brucella resides and replicates inside infected host cells in membrane-bound compartments called Brucella-containing vacuoles (BCVs). Following uptake, Brucella resides in endosomal BCVs (eBCVs) that gradually mature from early to late endosomal features. Through a poorly understood process that is key to the intracellular lifestyle of Brucella, the eBCV escapes fusion with lysosomes by transitioning to the replicative BCV (rBCV), a replicative niche directly connected to the endoplasmic reticulum (ER). Despite the notion that this complex intracellular lifestyle must depend on a multitude of host factors, a holistic view on which of these components control Brucella cell entry, trafficking, and replication is still missing. Here we used a systematic cell-based small interfering RNA (siRNA) knockdown screen in HeLa cells infected with Brucella abortus and identified 425 components of the human infectome for Brucella infection. These include multiple components of pathways involved in central processes such as the cell cycle, actin cytoskeleton dynamics, or vesicular trafficking. Using assays for pathogen entry, knockdown complementation, and colocalization at single-cell resolution, we identified the requirement of the VPS retromer for Brucella to escape the lysosomal degradative pathway and to establish its intracellular replicative niche. We thus validated the VPS retromer as a novel host factor critical for Brucella intracellular trafficking. Further, our genomewide data shed light on the interplay between central host processes and the biogenesis of the Brucella replicative niche.IMPORTANCE With >300,000 new cases of human brucellosis annually, Brucella is regarded as one of the most important zoonotic bacterial pathogens worldwide. The agent causing brucellosis resides inside host cells within vacuoles termed Brucella-containing vacuoles (BCVs). Although a few host components required to escape the degradative lysosomal pathway and to establish the ER-derived replicative BCV (rBCV) have already been identified, the global understanding of this highly coordinated process is still partial, and many factors remain unknown. To gain deeper insight into these fundamental questions, we performed a genomewide RNA interference (RNAi) screen aiming at discovering novel host factors involved in the Brucella intracellular cycle. We identified 425 host proteins that contribute to Brucella cellular entry, intracellular trafficking, and replication. Together, this study sheds light on previously unknown host pathways required for the Brucella infection cycle and highlights the VPS retromer components as critical factors for the establishment of the Brucella intracellular replicative niche.

Multi-omic measurements of heterogeneity in HeLa cells across laboratories.

Reproducibility in research can be compromised by both biological and technical variation, but most of the focus is on removing the latter. Here we investigate the effects of biological variation in HeLa cell lines using a systems-wide approach. We determine the degree of molecular and phenotypic variability across 14 stock HeLa samples from 13 international laboratories. We cultured cells in uniform conditions and profiled genome-wide copy numbers, mRNAs, proteins and protein turnover rates in each cell line. We discovered substantial heterogeneity between HeLa variants, especially between lines of the CCL2 and Kyoto varieties, and observed progressive divergence within a specific cell line over 50 successive passages. Genomic variability has a complex, nonlinear effect on transcriptome, proteome and protein turnover profiles, and proteotype patterns explain the varying phenotypic response of different cell lines to Salmonella infection. These findings have implications for the interpretation and reproducibility of research results obtained from human cultured cells.

Explicit Modeling of siRNA-Dependent On- and Off-Target Repression Improves the Interpretation of Screening Results.

RNAi is broadly used to map gene regulatory networks, but the identification of genes that are responsible for the observed phenotypes is challenging, as small interfering RNAs (siRNAs) simultaneously downregulate the intended on targets and many partially complementary off targets. Additionally, the scarcity of publicly available control datasets hinders the development and comparative evaluation of computational methods for analyzing the data. Here, we introduce PheLiM (https://github.com/andreariba/PheLiM), a method that uses predictions of siRNA on- and off-target downregulation to infer gene-specific contributions to phenotypes. To assess the performance of PheLiM, we carried out siRNA- and CRISPR/Cas9-based genome-wide screening of two well-characterized pathways, bone morphogenetic protein (BMP) and nuclear factor κB (NF-κB), and we reanalyzed publicly available siRNA screens. We demonstrate that PheLiM has the overall highest accuracy and most reproducible results compared to other available methods. PheLiM can accommodate various methods for predicting siRNA off targets and is broadly applicable to the identification of genes underlying complex phenotypes.

Microscopy-based Assays for High-throughput Screening of Host Factors Involved in Brucella Infection of Hela Cells.

Brucella species are facultative intracellular pathogens that infect animals as their natural hosts. Transmission to humans is most commonly caused by direct contact with infected animals or by ingestion of contaminated food and can lead to severe chronic infections. Brucella can invade professional and non-professional phagocytic cells and replicates within endoplasmic reticulum (ER)-derived vacuoles. The host factors required for Brucella entry into host cells, avoidance of lysosomal degradation, and replication in the ER-like compartment remain largely unknown. Here we describe two assays to identify host factors involved in Brucella entry and replication in HeLa cells. The protocols describe the use of RNA interference, while alternative screening methods could be applied. The assays are based on the detection of fluorescently labeled bacteria in fluorescently labeled host cells using automated wide-field microscopy. The fluorescent images are analyzed using a standardized image analysis pipeline in CellProfiler which allows single cell-based infection scoring. In the endpoint assay, intracellular replication is measured two days after infection. This allows bacteria to traffic to their replicative niche where proliferation is initiated around 12 hr after bacterial entry. Brucella which have successfully established an intracellular niche will thus have strongly proliferated inside host cells. Since intracellular bacteria will greatly outnumber individual extracellular or intracellular non-replicative bacteria, a strain constitutively expressing GFP can be used. The strong GFP signal is then used to identify infected cells. In contrast, for the entry assay it is essential to differentiate between intracellular and extracellular bacteria. Here, a strain encoding for a tetracycline-inducible GFP is used. Induction of GFP with simultaneous inactivation of extracellular bacteria by gentamicin enables the differentiation between intracellular and extracellular bacteria based on the GFP signal, with only intracellular bacteria being able to express GFP. This allows the robust detection of single intracellular bacteria before intracellular proliferation is initiated.

A Genome-Wide siRNA Screen Implicates Spire1/2 in SipA-Driven Salmonella Typhimurium Host Cell Invasion.

Salmonella Typhimurium (S. Tm) is a leading cause of diarrhea. The disease is triggered by pathogen invasion into the gut epithelium. Invasion is attributed to the SPI-1 type 3 secretion system (T1). T1 injects effector proteins into epithelial cells and thereby elicits rearrangements of the host cellular actin cytoskeleton and pathogen invasion. The T1 effector proteins SopE, SopB, SopE2 and SipA are contributing to this. However, the host cell factors contributing to invasion are still not completely understood. To address this question comprehensively, we used Hela tissue culture cells, a genome-wide siRNA library, a modified gentamicin protection assay and S. TmSipA, a sopBsopE2sopE mutant which strongly relies on the T1 effector protein SipA to invade host cells. We found that S. TmSipA invasion does not elicit membrane ruffles, nor promote the entry of non-invasive bacteria "in trans". However, SipA-mediated infection involved the SPIRE family of actin nucleators, besides well-established host cell factors (WRC, ARP2/3, RhoGTPases, COPI). Stage-specific follow-up assays and knockout fibroblasts indicated that SPIRE1 and SPIRE2 are involved in different steps of the S. Tm infection process. Whereas SPIRE1 interferes with bacterial binding, SPIRE2 influences intracellular replication of S. Tm. Hence, these two proteins might fulfill non-redundant functions in the pathogen-host interaction. The lack of co-localization hints to a short, direct interaction between S. Tm and SPIRE proteins or to an indirect effect.

Genome-Wide siRNA Screen Identifies Complementary Signaling Pathways Involved in Listeria Infection and Reveals Different Actin Nucleation Mechanisms during Listeria Cell Invasion and Actin Comet Tail Formation.

UNLABELLED: Listeria monocytogenes enters nonphagocytic cells by a receptor-mediated mechanism that is dependent on a clathrin-based molecular machinery and actin rearrangements. Bacterial intra- and intercellular movements are also actin dependent and rely on the actin nucleating Arp2/3 complex, which is activated by host-derived nucleation-promoting factors downstream of the cell receptor Met during entry and by the bacterial nucleation-promoting factor ActA during comet tail formation. By genome-wide small interfering RNA (siRNA) screening for host factors involved in bacterial infection, we identified diverse cellular signaling networks and protein complexes that support or limit these processes. In addition, we could precise previously described molecular pathways involved in Listeria invasion. In particular our results show that the requirements for actin nucleators during Listeria entry and actin comet tail formation are different. Knockdown of several actin nucleators, including SPIRE2, reduced bacterial invasion while not affecting the generation of comet tails. Most interestingly, we observed that in contrast to our expectations, not all of the seven subunits of the Arp2/3 complex are required for Listeria entry into cells or actin tail formation and that the subunit requirements for each of these processes differ, highlighting a previously unsuspected versatility in Arp2/3 complex composition and function. IMPORTANCE: Listeria is a bacterial pathogen that induces its internalization within the cytoplasm of human cells and has been used for decades as a major molecular tool to manipulate cells in order to explore and discover cellular functions. We have inactivated individually, for the first time in epithelial cells, all the genes of the human genome to investigate whether each gene modifies positively or negatively the Listeria infectious process. We identified novel signaling cascades that have never been associated with Listeria infection. We have also revisited the role of the molecular complex Arp2/3 involved in the polymerization of the actin cytoskeleton, which was shown previously to be required for Listeria entry and movement inside host cells, and we demonstrate that contrary to the general dogma, some subunits of the complex are dispensable for both Listeria entry and bacterial movement.

Autophagy Proteins Promote Repair of Endosomal Membranes Damaged by the Salmonella Type Three Secretion System 1.

Salmonella Typhimurium (S.Tm) is an enteropathogen requiring multiple virulence factors, including two type three secretion systems (T1 and T2). T1 triggers epithelium invasion in which the bacteria are taken up into endosomes that mature into Salmonella-containing vacuoles (SCV) and trigger T2 induction upon acidification. Mechanisms controlling endosome membrane integrity or pathogen egress into the cytosol are incompletely understood. We screened for host factors affecting invasion and SCV maturation and identified a role for autophagy in sealing endosomal membranes damaged by T1 during host cell invasion. S.Tm-infected autophagy-deficient (atg5(-/-)) cells exhibit reduced SCV dye retention and lower T2 expression but no effects on steps preceding SCV maturation. However, in the absence of T1, autophagy is dispensable for T2 induction. These findings establish a role of autophagy at early stages of S.Tm infection and suggest that autophagy-mediated membrane repair might be generally important for invasive pathogens and endosomal membrane function.

gespeR: a statistical model for deconvoluting off-target-confounded RNA interference screens.

Small interfering RNAs (siRNAs) exhibit strong off-target effects, which confound the gene-level interpretation of RNA interference screens and thus limit their utility for functional genomics studies. Here, we present gespeR, a statistical model for reconstructing individual, gene-specific phenotypes. Using 115,878 siRNAs, single and pooled, from three companies in three pathogen infection screens, we demonstrate that deconvolution of image-based phenotypes substantially improves the reproducibility between independent siRNA sets targeting the same genes. Genes selected and prioritized by gespeR are validated and shown to constitute biologically relevant components of pathogen entry mechanisms and TGF-ß signaling. gespeR is available as a Bioconductor R-package.

Imaging InlC secretion to investigate cellular infection by the bacterial pathogen Listeria monocytogenes.

Bacterial intracellular pathogens can be conceived as molecular tools to dissect cellular signaling cascades due to their capacity to exquisitely manipulate and subvert cell functions which are required for the infection of host target tissues. Among these bacterial pathogens, Listeria monocytogenes is a Gram positive microorganism that has been used as a paradigm for intracellular parasitism in the characterization of cellular immune responses, and which has played instrumental roles in the discovery of molecular pathways controlling cytoskeletal and membrane trafficking dynamics. In this article, we describe a robust microscopical assay for the detection of late cellular infection stages of L. monocytogenes based on the fluorescent labeling of InlC, a secreted bacterial protein which accumulates in the cytoplasm of infected cells; this assay can be coupled to automated high-throughput small interfering RNA screens in order to characterize cellular signaling pathways involved in the up- or down-regulation of infection.