MOTIVATION: The reconstruction of signaling pathways from gene knockdown data is a novel research field enabled by developments in RNAi screening technology. However, while RNA interference is a powerful technique to identify genes related to a phenotype of interest, their placement in the corresponding pathways remains a challenging problem. Difficulties are aggravated if not all pathway components can be observed after each knockdown, but readouts are only available for a small subset. We are then facing the problem of reconstructing a network from incomplete data. RESULTS: We infer pathway topologies from gene knockdown data using Bayesian networks with probabilistic Boolean threshold functions. To deal with the problem of under-determined network parameters, we employ a Bayesian learning approach, in which we can integrate arbitrary prior information on the network under consideration. Missing observations are integrated out. We compute the exact likelihood function for smaller networks, and use an approximation to evaluate the likelihood for larger networks. The posterior distribution is evaluated using mode hopping Markov chain Monte Carlo, distributions over topologies and parameters can then be used to design additional experiments. We evaluate our approach on a small artificial dataset, and present inference results on RNAi data from the Jak/Stat pathway in a human hepatoma cell line.

RNA interference (RNAi) is a powerful tool to study gene function in cultured cells. Transfected cell microarrays in principle allow high-throughput phenotypic analysis after gene knockdown by microscopy. But bottlenecks in imaging and data analysis have limited such high-content screens to endpoint assays in fixed cells and determination of global parameters such as viability. Here we have overcome these limitations and developed an automated platform for high-content RNAi screening by time-lapse fluorescence microscopy of live HeLa cells expressing histone-GFP to report on chromosome segregation and structure. We automated all steps, including printing transfection-ready small interfering RNA (siRNA) microarrays, fluorescence imaging and computational phenotyping of digital images, in a high-throughput workflow. We validated this method in a pilot screen assaying cell division and delivered a sensitive, time-resolved phenoprint for each of the 49 endogenous genes we suppressed. This modular platform is scalable and makes the power of time-lapse microscopy available for genome-wide RNAi screens.

West Nile virus (WNV), and related flaviviruses such as tick-borne encephalitis, Japanese encephalitis, yellow fever and dengue viruses, constitute a significant global human health problem1. However, our understanding of the molecular interaction of such flaviviruses with mammalian host cells is limited1. WNV encodes only 10 proteins, implying that it may use many cellular proteins for infection1. WNV enters the cytoplasm through pH-dependent endocytosis, undergoes cycles of translation and replication, assembles progeny virions in association with endoplasmic reticulum, and exits along the secretory pathway1, 2, 3. RNA interference (RNAi) presents a powerful forward genetics approach to dissect virus–host cell interactions4, 5, 6. Here we report the identification of 305 host proteins that affect WNV infection, using a human-genome-wide RNAi screen. Functional clustering of the genes revealed a complex dependence of this virus on host cell physiology, requiring a wide variety of molecules and cellular pathways for successful infection. We further demonstrate a requirement for the ubiquitin ligase CBLL1 in WNV internalization, a post-entry role for the endoplasmic-reticulum-associated degradation pathway in viral infection, and the monocarboxylic acid transporter MCT4 as a viral replication resistance factor. By extending this study to dengue virus, we show that flaviviruses have both overlapping and unique interaction strategies with host cells. This study provides a comprehensive molecular portrait of WNV–human cell interactions that forms a model for understanding single plus-stranded RNA virus infection, and reveals potential antiviral targets.

Insulin resistance is one of the dominant symptoms of type 2 diabetes (T2D). Although the molecular mechanisms leading to this resistance are largely unknown, experimental data support that the insulin signaling pathway is impaired in patients who are insulin resistant. To identify novel components/modulators of the insulin signaling pathway, we designed siRNAs targeting over 300 genes and tested the effects of knocking down these genes in an insulin-dependent, anti-lipolysis assay in 3T3-L1 adipocytes. For 126 genes, significant changes in free fatty acid release were observed. However, due to off-target effects (in addition to other limitations), high-throughput RNAi-based screens in cell-based systems generate significant amounts of noise. Therefore, to obtain a more reliable set of genes from the siRNA hits in our screen, we developed and applied a novel network-based approach that elucidates the mechanisms of action for the true positive siRNA hits. Our analysis results in the identification of a core network underlying the insulin signaling pathway that is more significantly enriched for genes previously associated with insulin resistance than the set of genes annotated in the KEGG database as belonging to the insulin signaling pathway. We experimentally validated one of the predictions, S1pr2, as a novel candidate gene for T2D.