Linking imaging to omics utilizing image-guided tissue extraction.

Phenotypic heterogeneity is commonly observed in diseased tissue, specifically in tumors. Multimodal imaging technologies can reveal tissue heterogeneity noninvasively in vivo, enabling imaging-based profiling of receptors, metabolism, morphology, or function on a macroscopic scale. In contrast, in vitro multiomics, immunohistochemistry, or histology techniques accurately characterize these heterogeneities in the cellular and subcellular scales in a more comprehensive but ex vivo manner. The complementary in vivo and ex vivo information would provide an enormous potential to better characterize a disease. However, this requires spatially accurate coregistration of these data by image-driven sampling as well as fast sample-preparation methods. Here, a unique image-guided milling machine and workflow for precise extraction of tissue samples from small laboratory animals or excised organs has been developed and evaluated. The samples can be delineated on tomographic images as volumes of interest and can be extracted with a spatial accuracy better than 0.25 mm. The samples remain cooled throughout the procedure to ensure metabolic stability, a precondition for accurate in vitro analysis.

SINCERITIES: inferring gene regulatory networks from time-stamped single cell transcriptional expression profiles.

MOTIVATION: Single cell transcriptional profiling opens up a new avenue in studying the functional role of cell-to-cell variability in physiological processes. The analysis of single cell expression profiles creates new challenges due to the distributive nature of the data and the stochastic dynamics of gene transcription process. The reconstruction of gene regulatory networks (GRNs) using single cell transcriptional profiles is particularly challenging, especially when directed gene-gene relationships are desired. RESULTS: We developed SINCERITIES (SINgle CEll Regularized Inference using TIme-stamped Expression profileS) for the inference of GRNs from single cell transcriptional profiles. We focused on time-stamped cross-sectional expression data, commonly generated from transcriptional profiling of single cells collected at multiple time points after cell stimulation. SINCERITIES recovers directed regulatory relationships among genes by employing regularized linear regression (ridge regression), using temporal changes in the distributions of gene expressions. Meanwhile, the modes of the gene regulations (activation and repression) come from partial correlation analyses between pairs of genes. We demonstrated the efficacy of SINCERITIES in inferring GRNs using in silico time-stamped single cell expression data and single cell transcriptional profiles of THP-1 monocytic human leukemia cells. The case studies showed that SINCERITIES could provide accurate GRN predictions, significantly better than other GRN inference algorithms such as TSNI, GENIE3 and JUMP3. Moreover, SINCERITIES has a low computational complexity and is amenable to problems of extremely large dimensionality. Finally, an application of SINCERITIES to single cell expression data of T2EC chicken erythrocytes pointed to BATF as a candidate novel regulator of erythroid development. AVAILABILITY AND IMPLEMENTATION: MATLAB and R version of SINCERITIES are freely available from the following websites: http://www.cabsel.ethz.ch/tools/sincerities.html and https://github.com/CABSEL/SINCERITIES. The single cell THP-1 and T2EC transcriptional profiles are available from the original publications (Kouno et al., 2013; Richard et al., 2016). The in silico single cell data are available on SINCERITIES websites. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

TRaCE+: Ensemble inference of gene regulatory networks from transcriptional expression profiles of gene knock-out experiments.

BACKGROUND: The inference of gene regulatory networks (GRNs) from transcriptional expression profiles is challenging, predominantly due to its underdetermined nature. One important consequence of underdetermination is the existence of many possible solutions to this inference. Our previously proposed ensemble inference algorithm TRaCE addressed this issue by inferring an ensemble of network directed graphs (digraphs) using differential gene expressions from gene knock-out (KO) experiments. However, TRaCE could not deal with the mode of the transcriptional regulations (activation or repression), an important feature of GRNs. RESULTS: In this work, we developed a new algorithm called TRaCE+ for the inference of an ensemble of signed GRN digraphs from transcriptional expression data of gene KO experiments. The sign of the edges indicates whether the regulation is an activation (positive) or a repression (negative). TRaCE+ generates the upper and lower bounds of the ensemble, which define uncertain regulatory interactions that could not be verified by the data. As demonstrated in the case studies using Escherichia coli GRN and 100-gene gold-standard GRNs from DREAM 4 network inference challenge, by accounting for regulatory signs, TRaCE+ could extract more information from the KO data than TRaCE, leading to fewer uncertain edges. Importantly, iterating TRaCE+ with an optimal design of gene KOs could resolve the underdetermined issue of GRN inference in much fewer KO experiments than using TRaCE. CONCLUSIONS: TRaCE+ expands the applications of ensemble GRN inference strategy by accounting for the mode of the gene regulatory interactions. In comparison to TRaCE, TRaCE+ enables a better utilization of gene KO data, thereby reducing the cost of tackling underdetermined GRN inference. TRaCE+ subroutines for MATLAB are freely available at the following website: http://www.cabsel.ethz.ch/tools/trace.html .

Optimal design of gene knockout experiments for gene regulatory network inference.

MOTIVATION: We addressed the problem of inferring gene regulatory network (GRN) from gene expression data of knockout (KO) experiments. This inference is known to be underdetermined and the GRN is not identifiable from data. Past studies have shown that suboptimal design of experiments (DOE) contributes significantly to the identifiability issue of biological networks, including GRNs. However, optimizing DOE has received much less attention than developing methods for GRN inference. RESULTS: We developed REDuction of UnCertain Edges (REDUCE) algorithm for finding the optimal gene KO experiment for inferring directed graphs (digraphs) of GRNs. REDUCE employed ensemble inference to define uncertain gene interactions that could not be verified by prior data. The optimal experiment corresponds to the maximum number of uncertain interactions that could be verified by the resulting data. For this purpose, we introduced the concept of edge separatoid which gave a list of nodes (genes) that upon their removal would allow the verification of a particular gene interaction. Finally, we proposed a procedure that iterates over performing KO experiments, ensemble update and optimal DOE. The case studies including the inference of Escherichia coli GRN and DREAM 4 100-gene GRNs, demonstrated the efficacy of the iterative GRN inference. In comparison to systematic KOs, REDUCE could provide much higher information return per gene KO experiment and consequently more accurate GRN estimates. CONCLUSIONS: REDUCE represents an enabling tool for tackling the underdetermined GRN inference. Along with advances in gene deletion and automation technology, the iterative procedure brings an efficient and fully automated GRN inference closer to reality. AVAILABILITY AND IMPLEMENTATION: MATLAB and Python scripts of REDUCE are available on www.cabsel.ethz.ch/tools/REDUCE CONTACT: rudi.gunawan@chem.ethz.ch SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Ensemble inference and inferability of gene regulatory networks.

The inference of gene regulatory network (GRN) from gene expression data is an unsolved problem of great importance. This inference has been stated, though not proven, to be underdetermined implying that there could be many equivalent (indistinguishable) solutions. Motivated by this fundamental limitation, we have developed new framework and algorithm, called TRaCE, for the ensemble inference of GRNs. The ensemble corresponds to the inherent uncertainty associated with discriminating direct and indirect gene regulations from steady-state data of gene knock-out (KO) experiments. We applied TRaCE to analyze the inferability of random GRNs and the GRNs of E. coli and yeast from single- and double-gene KO experiments. The results showed that, with the exception of networks with very few edges, GRNs are typically not inferable even when the data are ideal (unbiased and noise-free). Finally, we compared the performance of TRaCE with top performing methods of DREAM4 in silico network inference challenge.

Kinetics of optical control of enzyme activity with photoswitchable inhibitors.

This paper presents theoretical and simulation studies on controlling enzymatic reactions with photoswitchable inhibitors. It is found that the maximum attainable switching ratio (ratio of the steady state rates of product formation in the "on" and the "off" state) of a photoswitchable inhibitor is dependent on its photoswitching factor (ratio of the equilibrium constants of the photostationary states under the "off" and the "on" illuminations). Attachment of multiple photoswitchable groups to an inhibitor molecule increases the theoretically attainable switching ratio. The affinity of the enzyme for the substrate and the inhibitor is the rate-limiting factor of the switching between active and inactive states. Use of inhibitors with high enzyme affinity and photoswitchable groups with high photoswitching factor would provide high switching ratio. These results may help to design better systems for optical control of biochemical processes.