fastER: a user-friendly tool for ultrafast and robust cell segmentation in large-scale microscopy.

MOTIVATION: Quantitative large-scale cell microscopy is widely used in biological and medical research. Such experiments produce huge amounts of image data and thus require automated analysis. However, automated detection of cell outlines (cell segmentation) is typically challenging due to, e.g. high cell densities, cell-to-cell variability and low signal-to-noise ratios. RESULTS: Here, we evaluate accuracy and speed of various state-of-the-art approaches for cell segmentation in light microscopy images using challenging real and synthetic image data. The results vary between datasets and show that the tested tools are either not robust enough or computationally expensive, thus limiting their application to large-scale experiments. We therefore developed fastER, a trainable tool that is orders of magnitude faster while producing state-of-the-art segmentation quality. It supports various cell types and image acquisition modalities, but is easy-to-use even for non-experts: it has no parameters and can be adapted to specific image sets by interactively labelling cells for training. As a proof of concept, we segment and count cells in over 200 000 brightfield images (1388 × 1040 pixels each) from a six day time-lapse microscopy experiment; identification of over 46 000 000 single cells requires only about two and a half hours on a desktop computer. AVAILABILITY AND IMPLEMENTATION: C ++ code, binaries and data at https://www.bsse.ethz.ch/csd/software/faster.html . CONTACT: oliver.hilsenbeck@bsse.ethz.ch or timm.schroeder@bsse.ethz.ch. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Inferring causal metabolic signals that regulate the dynamic TORC1-dependent transcriptome.

Cells react to nutritional cues in changing environments via the integrated action of signaling, transcriptional, and metabolic networks. Mechanistic insight into signaling processes is often complicated because ubiquitous feedback loops obscure causal relationships. Consequently, the endogenous inputs of many nutrient signaling pathways remain unknown. Recent advances for system-wide experimental data generation have facilitated the quantification of signaling systems, but the integration of multi-level dynamic data remains challenging. Here, we co-designed dynamic experiments and a probabilistic, model-based method to infer causal relationships between metabolism, signaling, and gene regulation. We analyzed the dynamic regulation of nitrogen metabolism by the target of rapamycin complex 1 (TORC1) pathway in budding yeast. Dynamic transcriptomic, proteomic, and metabolomic measurements along shifts in nitrogen quality yielded a consistent dataset that demonstrated extensive re-wiring of cellular networks during adaptation. Our inference method identified putative downstream targets of TORC1 and putative metabolic inputs of TORC1, including the hypothesized glutamine signal. The work provides a basis for further mechanistic studies of nitrogen metabolism and a general computational framework to study cellular processes.

Accurate cell segmentation in microscopy images using membrane patterns.

MOTIVATION: Identifying cells in an image (cell segmentation) is essential for quantitative single-cell biology via optical microscopy. Although a plethora of segmentation methods exists, accurate segmentation is challenging and usually requires problem-specific tailoring of algorithms. In addition, most current segmentation algorithms rely on a few basic approaches that use the gradient field of the image to detect cell boundaries. However, many microscopy protocols can generate images with characteristic intensity profiles at the cell membrane. This has not yet been algorithmically exploited to establish more general segmentation methods. RESULTS: We present an automatic cell segmentation method that decodes the information across the cell membrane and guarantees optimal detection of the cell boundaries on a per-cell basis. Graph cuts account for the information of the cell boundaries through directional cross-correlations, and they automatically incorporate spatial constraints. The method accurately segments images of various cell types grown in dense cultures that are acquired with different microscopy techniques. In quantitative benchmarks and comparisons with established methods on synthetic and real images, we demonstrate significantly improved segmentation performance despite cell-shape irregularity, cell-to-cell variability and image noise. As a proof of concept, we monitor the internalization of green fluorescent protein-tagged plasma membrane transporters in single yeast cells. AVAILABILITY AND IMPLEMENTATION: Matlab code and examples are available at http://www.csb.ethz.ch/tools/cellSegmPackage.zip.

Near-optimal experimental design for model selection in systems biology.

MOTIVATION: Biological systems are understood through iterations of modeling and experimentation. Not all experiments, however, are equally valuable for predictive modeling. This study introduces an efficient method for experimental design aimed at selecting dynamical models from data. Motivated by biological applications, the method enables the design of crucial experiments: it determines a highly informative selection of measurement readouts and time points. RESULTS: We demonstrate formal guarantees of design efficiency on the basis of previous results. By reducing our task to the setting of graphical models, we prove that the method finds a near-optimal design selection with a polynomial number of evaluations. Moreover, the method exhibits the best polynomial-complexity constant approximation factor, unless P = NP. We measure the performance of the method in comparison with established alternatives, such as ensemble non-centrality, on example models of different complexity. Efficient design accelerates the loop between modeling and experimentation: it enables the inference of complex mechanisms, such as those controlling central metabolic operation. AVAILABILITY: Toolbox 'NearOED' available with source code under GPL on the Machine Learning Open Source Software Web site (mloss.org).

Using CellX to quantify intracellular events.

Methods to quantify features of individual cells using light microscopy have become widely used in biology. A multitude of computational tools has been developed for image analysis; however, they are often only for specific cell types and microscopy techniques. This unit describes CellX, an open-source software package for cell segmentation, intensity quantification, and cell tracking on a variety of microscopy images. CellX can perform cell segmentation largely independently of cell shapes, and can also cope with images that are crowded with cells. The basic protocol describes how to use CellX for cell segmentation and quantification. This protocol remains the same whether there is a collection of images to be analyzed or whether cell tracking on a sequence of images is to be performed. The CellX output comprises control images for visual validation, text files for post-processing statistics, and MATLAB objects for advanced subsequent analysis.

A specialized ODE integrator for the efficient computation of parameter sensitivities.

BACKGROUND: Dynamic mathematical models in the form of systems of ordinary differential equations (ODEs) play an important role in systems biology. For any sufficiently complex model, the speed and accuracy of solving the ODEs by numerical integration is critical. This applies especially to systems identification problems where the parameter sensitivities must be integrated alongside the system variables. Although several very good general purpose ODE solvers exist, few of them compute the parameter sensitivities automatically. RESULTS: We present a novel integration algorithm that is based on second derivatives and contains other unique features such as improved error estimates. These features allow the integrator to take larger time steps than other methods. In practical applications, i.e. systems biology models of different sizes and behaviors, the method competes well with established integrators in solving the system equations, and it outperforms them significantly when local parameter sensitivities are evaluated. For ease-of-use, the solver is embedded in a framework that automatically generates the integrator input from an SBML description of the system of interest. CONCLUSIONS: For future applications, comparatively 'cheap' parameter sensitivities will enable advances in solving large, otherwise computationally expensive parameter estimation and optimization problems. More generally, we argue that substantially better computational performance can be achieved by exploiting characteristics specific to the problem domain; elements of our methods such as the error estimation could find broader use in other, more general numerical algorithms.

Systems analysis of cellular networks under uncertainty.

Besides the often-quoted complexity of cellular networks, the prevalence of uncertainties about components, interactions, and their quantitative features provides a largely underestimated hallmark of current systems biology. This uncertainty impedes the development of mechanistic mathematical models to achieve a true systems-level understanding. However, there is increasing evidence that theoretical approaches from diverse scientific domains can extract relevant biological knowledge efficiently, even from poorly characterized biological systems. As a common denominator, the methods focus on structural, rather than more detailed, kinetic network properties. A deeper understanding, better scaling, and the ability to combine the approaches pose formidable challenges for future theory developments.

DNA replication in the fission yeast: robustness in the face of uncertainty.

DNA replication, the process of duplication of a cell's genetic content, must be carried out with great precision every time the cell divides, so that genetic information is preserved. Control mechanisms must ensure that every base of the genome is replicated within the allocated time (S-phase) and only once per cell cycle, thereby safeguarding genomic integrity. In eukaryotes, replication starts from many points along the chromosome, termed origins of replication, and then proceeds continuously bidirectionally until an opposing moving fork is encountered. In contrast to bacteria, where a specific site on the genome serves as an origin in every cell division, in most eukaryotes origin selection appears highly stochastic: many potential origins exist, of which only a subset is selected to fire in any given cell, giving rise to an apparently random distribution of initiation events across the genome. Origin states change throughout the cell cycle, through the ordered formation and modification of origin-associated multisubunit protein complexes. State transitions are governed by fluctuations of cyclin-dependent kinase (CDK) activity and guards in these transitions ensure system memory. We present here DNA replication dynamics, emphasizing recent data from the fission yeast Schizosaccharomyces pombe, and discuss how robustness may be ensured in spite of (or even assisted by) system randomness.