BioSig: an imaging bioinformatics system for phenotypic analysis

BioSig3D: High Content Screening of Three-Dimensional Cell Culture Models

BioSig3D is a computational platform for high-content screening of three-dimensional (3D) cell culture models that are imaged in full 3D volume. It provides an end-to-end solution for designing high content screening assays, based on colony organization that is derived from segmentation of nuclei in each colony. BioSig3D also enables visualization of raw and processed 3D volumetric data for quality control, and integrates advanced bioinformatics analysis. The system consists of multiple computational and annotation modules that are coupled together with a strong use of controlled vocabularies to reduce ambiguities between different users. It is a web-based system that allows users to: design an experiment by defining experimental variables, upload a large set of volumetric images into the system, analyze and visualize the dataset, and either display computed indices as a heatmap, or phenotypic subtypes for heterogeneity analysis, or download computed indices for statistical analysis or integrative biology. BioSig3D has been used to profile baseline colony formations with two experiments: (i) morphogenesis of a panel of human mammary epithelial cell lines (HMEC), and (ii) heterogeneity in colony formation using an immortalized non-transformed cell line. These experiments reveal intrinsic growth properties of well-characterized cell lines that are routinely used for biological studies. BioSig3D is being released with seed datasets and video-based documentation.

Identification of fluorescent compounds with non-specific binding property via high throughput live cell microscopy

INTRODUCTION

Compounds exhibiting low non-specific intracellular binding or non-stickiness are concomitant with rapid clearing and in high demand for live-cell imaging assays because they allow for intracellular receptor localization with a high signal/noise ratio. The non-stickiness property is particularly important for imaging intracellular receptors due to the equilibria involved.

METHOD

Three mammalian cell lines with diverse genetic backgrounds were used to screen a combinatorial fluorescence library via high throughput live cell microscopy for potential ligands with high in- and out-flux properties. The binding properties of ligands identified from the first screen were subsequently validated on plant root hair. A correlative analysis was then performed between each ligand and its corresponding physiochemical and structural properties.

RESULTS

The non-stickiness property of each ligand was quantified as a function of the temporal uptake and retention on a cell-by-cell basis. Our data shows that (i) mammalian systems can serve as a pre-screening tool for complex plant species that are not amenable to high-throughput imaging; (ii) retention and spatial localization of chemical compounds vary within and between each cell line; and (iii) the structural similarities of compounds can infer their non-specific binding properties.

CONCLUSION

We have validated a protocol for identifying chemical compounds with non-specific binding properties that is testable across diverse species. Further analysis reveals an overlap between the non-stickiness property and the structural similarity of compounds. The net result is a more robust screening assay for identifying desirable ligands that can be used to monitor intracellular localization. Several new applications of the screening protocol and results are also presented.

Morphometic analysis of TCGA glioblastoma multiforme

Background

Our goals are to develop a computational histopathology pipeline for characterizing tumor types that are being generated by The Cancer Genome Atlas (TCGA) for genomic association. TCGA is a national collaborative program where different tumor types are being collected, and each tumor is being characterized using a variety of genome-wide platforms. Here, we have developed a tumor-centric analytical pipeline to process tissue sections stained with hematoxylin and eosin (H&E) for visualization and cell-by-cell quantitative analysis. Thus far, analysis is limited to Glioblastoma Multiforme (GBM) and kidney renal clear cell carcinoma tissue sections. The final results are being distributed for subtyping and linking the histology sections to the genomic data.

Results

A computational pipeline has been designed to continuously update a local image database, with limited clinical information, from an NIH repository. Each image is partitioned into blocks, where each cell in the block is characterized through a multidimensional representation (e.g., nuclear size, cellularity). A subset of morphometric indices, representing potential underlying biological processes, can then be selected for subtyping and genomic association. Simultaneously, these subtypes can also be predictive of the outcome as a result of clinical treatments. Using the cellularity index and nuclear size, the computational pipeline has revealed five subtypes, and one subtype, corresponding to the extreme high cellularity, has shown to be a predictor of survival as a result of a more aggressive therapeutic regime. Further association of this subtype with the corresponding gene expression data has identified enrichment of (i) the immune response and AP-1 signaling pathways, and (ii) IFNG, TGFB1, PKC, Cytokine, and MAPK14 hubs.

Conclusion

While subtyping is often performed with genome-wide molecular data, we have shown that it can also be applied to categorizing histology sections. Accordingly, we have identified a subtype that is a predictor of the outcome as a result of a therapeutic regime. Computed representation has become publicly available through our Web site.

Persistence of γ-H2AX and 53BP1 foci in proliferating and non-proliferating human mammary epithelial cells after exposure to γ-rays or iron ions

PURPOSE

To investigate γ-H2AX (phosphorylated histone H2AX) and 53BP1 (tumour protein 53 binding protein No. 1) foci formation and removal in proliferating and non-proliferating human mammary epithelial cells (HMEC) after exposure to sparsely and densely ionising radiation under different cell culture conditions.

MATERIAL AND METHODS

HMEC cells were grown either as monolayers (2D) or in extracellular matrix to allow the formation of acinar structures in vitro (3D). Foci numbers were quantified by image analysis at various time points after exposure.

RESULTS

Our results reveal that in non-proliferating cells under 2D and 3D cell culture conditions, iron-ion induced γ-H2AX foci were still present at 72 h after exposure, although 53BP1 foci returned to control levels at 48 h. In contrast in proliferating HMEC, both γ-H2AX and 53BP1 foci decreased to control levels during the 24-48 h time interval after irradiation under 2D conditions. Foci numbers decreased faster after γ-ray irradiation and returned to control levels by 12 h regardless of marker, cell proliferation status, and cell culture condition.

CONCLUSIONS

The disappearance of radiation-induced γ-H2AX and 53BP1 foci in HMEC has different dynamics that depend on radiation quality and proliferation status. Notably, the general patterns do not depend on the cell culture condition (2D versus 3D). We speculate that the persistent γ-H2AX foci in iron-ion irradiated non-proliferating cells could be due to limited availability of double-strand break (DSB) repair pathways in G0/G1-phase, or that repair of complex DSB requires replication or chromatin remodelling.

Sparse multitask regression for identifying common mechanism of response to therapeutic targets

Motivation: Molecular association of phenotypic responses is an important step in hypothesis generation and for initiating design of new experiments. Current practices for associating gene expression data with multidimensional phenotypic data are typically (i) performed one-to-one, i.e. each gene is examined independently with a phenotypic index and (ii) tested with one stress condition at a time, i.e. different perturbations are analyzed separately. As a result, the complex coordination among the genes responsible for a phenotypic profile is potentially lost. More importantly, univariate analysis can potentially hide new insights into common mechanism of response.

Results: In this article, we propose a sparse, multitask regression model together with co-clustering analysis to explore the intrinsic grouping in associating the gene expression with phenotypic signatures. The global structure of association is captured by learning an intrinsic template that is shared among experimental conditions, with local perturbations introduced to integrate effects of therapeutic agents. We demonstrate the performance of our approach on both synthetic and experimental data. Synthetic data reveal that the multi-task regression has a superior reduction in the regression error when compared with traditional L₁-and L₂-regularized regression. On the other hand, experiments with cell cycle inhibitors over a panel of 14 breast cancer cell lines demonstrate the relevance of the computed molecular predictors with the cell cycle machinery, as well as the identification of hidden variables that are not captured by the baseline regression analysis. Accordingly, the system has identified CLCA2 as a hidden transcript and as a common mechanism of response for two therapeutic agents of CI-1040 and Iressa, which are currently in clinical use.

Contact:b_parvin@lbl.gov

Linking changes in epithelial morphogenesis to cancer mutations using computational modeling

Most tumors arise from epithelial tissues, such as mammary glands and lobules, and their initiation is associated with the disruption of a finely defined epithelial architecture. Progression from intraductal to invasive tumors is related to genetic mutations that occur at a subcellular level but manifest themselves as functional and morphological changes at the cellular and tissue scales, respectively. Elevated proliferation and loss of epithelial polarization are the two most noticeable changes in cell phenotypes during this process. As a result, many three-dimensional cultures of tumorigenic clones show highly aberrant morphologies when compared to regular epithelial monolayers enclosing the hollow lumen (acini). In order to shed light on phenotypic changes associated with tumor cells, we applied the bio-mechanical IBCell model of normal epithelial morphogenesis quantitatively matched to data acquired from the non-tumorigenic human mammary cell line, MCF10A. We then used a high-throughput simulation study to reveal how modifications in model parameters influence changes in the simulated architecture. Three parameters have been considered in our study, which define cell sensitivity to proliferative, apoptotic and cell-ECM adhesive cues. By mapping experimental morphologies of four MCF10A-derived cell lines carrying different oncogenic mutations onto the model parameter space, we identified changes in cellular processes potentially underlying structural modifications of these mutants. As a case study, we focused on MCF10A cells expressing an oncogenic mutant HER2-YVMA to quantitatively assess changes in cell doubling time, cell apoptotic rate, and cell sensitivity to ECM accumulation when compared to the parental non-tumorigenic cell line. By mapping in vitro mutant morphologies onto in silico ones we have generated a means of linking the morphological and molecular scales via computational modeling. Thus, IBCell in combination with 3D acini cultures can form a computational/experimental platform for suggesting the relationship between the histopathology of neoplastic lesions and their underlying molecular defects.

Bioimage informatics for experimental biology

Over the past twenty years there have been great advances in light microscopy with the result that multidimensional imaging has driven a revolution in modern biology. The development of new approaches of data acquisition is reported frequently, and yet the significant data management and analysis challenges presented by these new complex datasets remain largely unsolved. As in the well-developed field of genome bioinformatics, central repositories are and will be key resources, but there is a critical need for informatics tools in individual laboratories to help manage, share, visualize, and analyze image data. In this article we present the recent efforts by the bioimage informatics community to tackle these challenges, and discuss our own vision for future development of bioimage informatics solutions.

Screening the cellular microenvironment: a role for microfluidics

The cellular microenvironment is an increasingly discussed topic in cell biology as it has been implicated in the progression of cancer and the maintenance of stem cells. The microenvironment of a cell is an organized combination of extracellular matrix (ECM), cells, and interstitial fluid that influence cellular phenotype through physical, mechanical, and biochemical mechanisms. Screening can be used to map combinations of cells and microenvironments to phenotypic outcomes in a way that can help develop more predictive in vitro models and to better understand phenotypic mechanisms from a systems biology perspective. This paper examines microenvironmental screening in terms of outcomes and benefits, key elements of the screening process, challenges for implementation, and a possible role for microfluidics as the screening platform. To assess microfluidics for use in microenvironmental screening, examples and categories of micro-scale and microfluidic technology are highlighted. Microfluidic technology shows promise for simultaneous control of multiple parameters of the microenvironment and can provide a base for scaling advanced cell-based experiments into automated high-throughput formats.

Geometric approach to segmentation and protein localization in cell culture assays

Cell-based fluorescence imaging assays are heterogeneous and require the collection of a large number of images for detailed quantitative analysis. Complexities arise as a result of variation in spatial nonuniformity, shape, overlapping compartments and scale (size). A new technique and methodology has been developed and tested for delineating subcellular morphology and partitioning overlapping compartments at multiple scales. This system is packaged as an integrated software platform for quantifying images that are obtained through fluorescence microscopy. Proposed methods are model based, leveraging geometric shape properties of subcellular compartments and corresponding protein localization. From the morphological perspective, convexity constraint is imposed to delineate and partition nuclear compartments. From the protein localization perspective, radial symmetry is imposed to localize punctate protein events at submicron resolution. Convexity constraint is imposed against boundary information, which are extracted through a combination of zero-crossing and gradient operator. If the convexity constraint fails for the boundary then positive curvature maxima are localized along the contour and the entire blob is partitioned into disjointed convex objects representing individual nuclear compartment, by enforcing geometric constraints. Nuclear compartments provide the context for protein localization, which may be diffuse or punctate. Punctate signal are localized through iterative voting and radial symmetries for improved reliability and robustness. The technique has been tested against 196 images that were generated to study centrosome abnormalities. Corresponding computed representations are compared against manual counts for validation.

Iterative voting for inference of structural saliency and characterization of subcellular events

Saliency is an important perceptual cue that occurs at different levels of resolution. Important attributes of saliency are symmetry, continuity, and closure. Detection of these attributes is often hindered by noise, variation in scale, and incomplete information. This paper introduces the iterative voting method, which uses oriented kernels for inferring saliency as it relates to symmetry. A unique aspect of the technique is the kernel topography, which is refined and reoriented iteratively. The technique can cluster and group nonconvex perceptual circular symmetries along the radial line of an object's shape. It has an excellent noise immunity and is shown to be tolerant to perturbation in scale. The application of this technique to images obtained through various modes of microscopy is demonstrated. Furthermore, as a case example, the method has been applied to quantify kinetics of nuclear foci formation that are formed by phosphorylation of histone gammaH2AX following ionizing radiation. Iterative voting has been implemented in both 2-D and 3-D for multi image analysis.

Systems cell biology based on high-content screening

A new discipline of biology has emerged since 2004, which we call "systems cell biology" (SCB). Systems cell biology is the study of the living cell, the basic unit of life, an integrated and interacting network of genes, proteins, and myriad metabolic reactions that give rise to function. SCB takes advantage of high-content screening platforms, but delivers more detailed profiles of cellular systemic function, including the application of advanced reagents and informatics tools to sophisticated cellular models. Therefore, an SCB profile is a cellular systemic response as measured by a panel of reagents that quantify a specific set of biomarkers.