Gene network inference via structural equation modeling in genetical genomics experiments

Our goal is gene network inference in genetical genomics or systems genetics experiments. For species where sequence information is available, we first perform expression quantitative trait locus (eQTL) mapping by jointly utilizing cis-, cis-trans-, and trans-regulation. After using local structural models to identify regulator-target pairs for each eQTL, we construct an encompassing directed network (EDN) by assembling all retained regulator-target relationships. The EDN has nodes corresponding to expressed genes and eQTL and directed edges from eQTL to cis-regulated target genes, from cis-regulated genes to cis-trans-regulated target genes, from trans-regulator genes to target genes, and from trans-eQTL to target genes. For network inference within the strongly constrained search space defined by the EDN, we propose structural equation modeling (SEM), because it can model cyclic networks and the EDN indeed contains feedback relationships. On the basis of a factorization of the likelihood and the constrained search space, our SEM algorithm infers networks involving several hundred genes and eQTL. Structure inference is based on a penalized likelihood ratio and an adaptation of Occam's window model selection. The SEM algorithm was evaluated using data simulated with nonlinear ordinary differential equations and known cyclic network topologies and was applied to a real yeast data set.     

Inferring Regulatory Networks from Expression Data Using Tree-Based Methods

One of the pressing open problems of computational systems biology is the elucidation of the topology of genetic regulatory networks (GRNs) using high throughput genomic data, in particular microarray gene expression data. The Dialogue for Reverse Engineering Assessments and Methods (DREAM) challenge aims to evaluate the success of GRN inference algorithms on benchmarks of simulated data. In this article, we present GENIE3, a new algorithm for the inference of GRNs that was best performer in the DREAM4 In Silico Multifactorial challenge. GENIE3 decomposes the prediction of a regulatory network between p genes into p different regression problems. In each of the regression problems, the expression pattern of one of the genes (target gene) is predicted from the expression patterns of all the other genes (input genes), using tree-based ensemble methods Random Forests or Extra-Trees. The importance of an input gene in the prediction of the target gene expression pattern is taken as an indication of a putative regulatory link. Putative regulatory links are then aggregated over all genes to provide a ranking of interactions from which the whole network is reconstructed. In addition to performing well on the DREAM4 In Silico Multifactorial challenge simulated data, we show that GENIE3 compares favorably with existing algorithms to decipher the genetic regulatory network of Escherichia coli. It doesn''t make any assumption about the nature of gene regulation, can deal with combinatorial and non-linear interactions, produces directed GRNs, and is fast and scalable. In conclusion, we propose a new algorithm for GRN inference that performs well on both synthetic and real gene expression data. The algorithm, based on feature selection with tree-based ensemble methods, is simple and generic, making it adaptable to other types of genomic data and interactions.     

Programmable control of bacterial gene expression with the combined CRISPR and antisense RNA system

A central goal of synthetic biology is to implement diverse cellular functions by predictably controlling gene expression. Though research has focused more on protein regulators than RNA regulators, recent advances in our understanding of RNA folding and functions have motivated the use of RNA regulators. RNA regulators provide an advantage because they are easier to design and engineer than protein regulators, potentially have a lower burden on the cell and are highly orthogonal. Here, we combine the CRISPR system from Streptococcus pyogenes and synthetic antisense RNAs (asRNAs) in Escherichia coli strains to repress or derepress a target gene in a programmable manner. Specifically, we demonstrate for the first time that the gene target repressed by the CRISPR system can be derepressed by expressing an asRNA that sequesters a small guide RNA (sgRNA). Furthermore, we demonstrate that tunable levels of derepression can be achieved (up to 95%) by designing asRNAs that target different regions of a sgRNA and by altering the hybridization free energy of the sgRNA–asRNA complex. This new system, which we call the combined CRISPR and asRNA system, can be used to reversibly repress or derepress multiple target genes simultaneously, allowing for rational reprogramming of cellular functions.