An approach for extensibly profiling the molecular states of cellular subpopulations

A statistical framework for high-content phenotypic profiling using cellular feature distributions

High-content screening (HCS) uses microscopy images to generate phenotypic profiles of cell morphological data in high-dimensional feature space. While HCS provides detailed cytological information at single-cell resolution, these complex datasets are usually aggregated into summary statistics that do not leverage patterns of biological variability within cell populations. Here we present a broad-spectrum HCS analysis system that measures image-based cell features from 10 cellular compartments across multiple assay panels. We introduce quality control measures and statistical strategies to streamline and harmonize the data analysis workflow, including positional and plate effect detection, biological replicates analysis and feature reduction. We also demonstrate that the Wasserstein distance metric is superior over other measures to detect differences between cell feature distributions. With this workflow, we define per-dose phenotypic fingerprints for 65 mechanistically diverse compounds, provide phenotypic path visualizations for each compound and classify compounds into different activity groups.

GRAN3SAT: Creating Flexible Higher-Order Logic Satisfiability in the Discrete Hopfield Neural Network

One of the main problems in representing information in the form of nonsystematic logic is the lack of flexibility, which leads to potential overfitting. Although nonsystematic logic improves the representation of the conventional k Satisfiability, the formulations of the first, second, and third-order logical structures are very predictable. This paper proposed a novel higher-order logical structure, named G-Type Random k Satisfiability, by capitalizing the new random feature of the first, second, and third-order clauses. The proposed logic was implemented into the Discrete Hopfield Neural Network as a symbolic logical rule. The proposed logic in Discrete Hopfield Neural Networks was evaluated using different parameter settings, such as different orders of clauses, different proportions between positive and negative literals, relaxation, and differing numbers of learning trials. Each evaluation utilized various performance metrics, such as learning error, testing error, weight error, energy analysis, and similarity analysis. In addition, the flexibility of the proposed logic was compared with current state-of-the-art logic rules. Based on the simulation, the proposed logic was reported to be more flexible, and produced higher solution diversity.

Image-based cell profiling enhancement via data cleaning methods

With the advent of high-throughput assays, a large number of biological experiments can be carried out. Image-based assays are among the most accessible and inexpensive technologies for this purpose. Indeed, these assays have proved to be effective in characterizing unknown functions of genes and small molecules. Image analysis pipelines have a pivotal role in translating raw images that are captured in such assays into useful and compact representation, also known as measurements. CellProfiler is a popular and commonly used tool for this purpose through providing readily available modules for the cell/nuclei segmentation, and making various measurements, or features, for each cell/nuclei. Single cell features are then aggregated for each treatment replica to form treatment “profiles”. However, there may be several sources of error in the CellProfiler quantification pipeline that affects the downstream analysis that is performed on the profiles. In this work, we examined various preprocessing approaches to improve the profiles. We consider the identification of drug mechanisms of action as the downstream task to evaluate such preprocessing approaches. Our enhancement steps mainly consist of data cleaning, cell level outlier detection, toxic drug detection, and regressing out the cell area from all other features, as many of them are widely affected by the cell area. Our experiments indicate that by performing these time-efficient preprocessing steps, image-based profiles can preserve more meaningful information compared to raw profiles. In the end, we also suggest possible avenues for future research.

Automated and unbiased discrimination of ALS from control tissue at single cell resolution

Histopathological analysis of tissue sections is invaluable in neurodegeneration research. However, cell-to-cell variation in both the presence and severity of a given phenotype is a key limitation of this approach, reducing the signal to noise ratio and leaving unresolved the potential of single-cell scoring for a given disease attribute. Here, we tested different machine learning methods to analyse high-content microscopy measurements of hundreds of motor neurons (MNs) from amyotrophic lateral sclerosis (ALS) post-mortem tissue sections. Furthermore, we automated the identification of phenotypically distinct MN subpopulations in VCP- and SOD1-mutant transgenic mice, revealing common morphological cellular phenotypes. Additionally we established scoring metrics to rank cells and tissue samples for both disease probability and severity. By adapting this paradigm to human post-mortem tissue, we validated our core finding that morphological descriptors robustly discriminate ALS from control healthy tissue at single cell resolution. Determining disease presence, severity and unbiased phenotypes at single cell resolution might prove transformational in our understanding of ALS and neurodegeneration more broadly.

Single-cell lipidomics with high structural specificity by mass spectrometry

Single-cell analysis is critical to revealing cell-to-cell heterogeneity that would otherwise be lost in ensemble analysis. Detailed lipidome characterization for single cells is still far from mature, especially when considering the highly complex structural diversity of lipids and the limited sample amounts available from a single cell. We report the development of a general strategy enabling single-cell lipidomic analysis with high structural specificity. Cell fixation is applied to retain lipids in the cell during batch treatments prior to single-cell analysis. In addition to tandem mass spectrometry analysis revealing the class and fatty acyl-chain for lipids, batch photochemical derivatization and single-cell droplet treatment are performed to identify the C=C locations and sn-positions of lipids, respectively. Electro-migration combined with droplet-assisted electrospray ionization enables single-cell mass spectrometry analysis with easy operation but high efficiency in sample usage. Four subtypes of human breast cancer cells are correctly classified through quantitative analysis of lipid C=C location or sn-position isomers in ~160 cells. Most importantly, the single-cell deep lipidomics strategy successfully discriminates gefitinib-resistant cells from a population of wild-type human lung cancer cells (HCC827), highlighting its unique capability to promote precision medicine.

Use of artificial intelligence to enhance phenotypic drug discovery

Research and development (R&D) productivity across the pharmaceutical industry has received close scrutiny over the past two decades, especially taking into consideration reports of attrition rates and the colossal cost for drug development. The respective merits of the two main drug discovery approaches, phenotypic and target based, have divided opinion across the research community, because each hold different advantages for identifying novel molecular entities with a successful path to the market. Nevertheless, both have low translatability in the clinic. Artificial intelligence (AI) and adoption of machine learning (ML) tools offer the promise of revolutionising drug development, and overcoming obstacles in the drug discovery pipeline. Here, we assess the potential of target-driven and phenotypic-based approaches and offer a holistic description of the current state of the field, from both a scientific and industry perspective. With the emerging partnerships between AI/ML and pharma still in their relative infancy, we investigate the potential and current limitations with a particular focus on phenotypic drug discovery. Finally, we emphasise the value of public-private partnerships (PPPs) and cross-disciplinary collaborations to foster innovation and facilitate efficient drug discovery programmes.

Optimum concentration-response curve metrics for supervised selection of discriminative cellular phenotypic endpoints for chemical hazard assessment

High-content imaging (HCI) provides quantitative and information-rich measurements of chemical effects on human in vitro cell models. Identification of discriminative phenotypic endpoints from cellular features obtained from HCI is required for accurate assessments of potential chemical hazards. However, the use of suboptimal metrics to quantify the concentration-response curves (CRC) of chemicals based on these features may obscure discriminative features, and lead to non-predictive endpoints and poor chemical classifications or hazard assessments. Here, we present a systematic and data-driven study on the performances of different CRC metrics in identifying image-based phenotypic features that can accurately classify the effects of reference chemicals with known in vivo toxicities. We studied four previous HCI in vitro nephro- or pulmono-toxicity datasets, which contain phenotypic feature measurements from different cell and feature types. Within a feature type, we found that efficacy metrics at higher chemical concentrations tend to give higher classification accuracy, whereas potency metrics do not have obvious trends across different response levels. Across different cell and feature types, efficacy metrics generally gave higher classification accuracy than potency metrics and area under the curve (AUC). Our results suggest that efficacy metrics, especially at higher concentrations, are more likely to help us to identify discriminative phenotypic endpoints. Therefore, HCI experiments for toxicological applications should include measurements at sufficiently high chemical concentrations, and efficacy metrics should always be analyzed. The identified features may be used as specific toxicity endpoints for further chemical hazard assessment.

Improving drug discovery using image-based multiparametric analysis of the epigenetic landscape

High-content phenotypic screening has become the approach of choice for drug discovery due to its ability to extract drug-specific multi-layered data. In the field of epigenetics, such screening methods have suffered from a lack of tools sensitive to selective epigenetic perturbations. Here we describe a novel approach, Microscopic Imaging of Epigenetic Landscapes (MIEL), which captures the nuclear staining patterns of epigenetic marks and employs machine learning to accurately distinguish between such patterns. We validated the MIEL platform across multiple cells lines and using dose-response curves, to insure the fidelity and robustness of this approach for high content high throughput drug discovery. Focusing on noncytotoxic glioblastoma treatments, we demonstrated that MIEL can identify and classify epigenetically active drugs. Furthermore, we show MIEL was able to accurately rank candidate drugs by their ability to produce desired epigenetic alterations consistent with increased sensitivity to chemotherapeutic agents or with induction of glioblastoma differentiation.

Capturing single-cell heterogeneity via data fusion improves image-based profiling

Single-cell resolution technologies warrant computational methods that capture cell heterogeneity while allowing efficient comparisons of populations. Here, we summarize cell populations by adding features’ dispersion and covariances to population averages, in the context of image-based profiling. We find that data fusion is critical for these metrics to improve results over the prior alternatives, providing at least ~20% better performance in predicting a compound’s mechanism of action (MoA) and a gene’s pathway.

A challenge with single-cell resolution methods is that cell heterogeneity should be captured while allowing for comparisons between populations. Here the authors fuse information from the dispersion profiles with the average profiles at the level of profiles’ similarity matrices for single cell imaging data.

Chemical biology approaches targeting the actin cytoskeleton through phenotypic screening

The actin cytoskeleton is dysregulated in cancer, yet this critical cellular machinery has not translated as a druggable clinical target due to cardio-toxic side-effects. Many actin regulators are also considered undruggable, being structural proteins lacking clear functional sites suitable for targeted drug design. In this review, we discuss opportunities and challenges associated with drugging the actin cytoskeleton through its structural regulators, taking tropomyosins as a target example. In particular, we highlight emerging data acquisition and analysis trends driving phenotypic, imaging-based compound screening. Finally, we consider how the confluence of these trends is now bringing functionally integral machineries such as the actin cytoskeleton, and associated structural regulatory proteins, into an expanded repertoire of druggable targets with previously unexploited clinical potential.

Integrating Analysis of Cellular Heterogeneity in High-Content Dose-Response Studies

Heterogeneity is a complex property of cellular systems and therefore presents challenges to the reliable identification and characterization. Large-scale biology projects may span many months, requiring a systematic approach to quality control to track reproducibility and correct for instrumental variation and assay drift that could mask biological heterogeneity and preclude comparisons of heterogeneity between runs or even between plates. However, presently there is no standard approach to the tracking and analysis of heterogeneity. Previously, we demonstrated the use of the Kolmogorov-Smirnov statistic as a metric for monitoring the reproducibility of heterogeneity in a screen and described the use of three heterogeneity indices as a means to characterize, filter, and browse cellular heterogeneity in big data sets (Gough et al., Methods 96:12-26, 2016). In this chapter, we present a detailed method for integrating the analysis of cellular heterogeneity in assay development, validation, screening, and post screen. Importantly, we provide a detailed method for quality control, to normalize cellular data, track heterogeneity over time, and analyze heterogeneity in big data sets, along with software tools to assist in that process. The example screen for this method is from an HCS project, but the approach applies equally to other experimental methods that measure populations of cells.

Compound Functional Prediction Using Multiple Unrelated Morphological Profiling Assays

Phenotypic cell-based assays have proven to be efficient at discovering first-in-class therapeutic drugs mainly because they allow for scanning a wide spectrum of possible targets at once. However, despite compelling methodological advances, posterior identification of a compound's mechanism of action (MOA) has remained difficult and highly refractory to automated analyses. Methods such as the cell painting assay and multiplexing fluorescent dyes to reveal broadly relevant cellular components were recently suggested for MOA prediction. We demonstrated that adding fluorescent dyes to a single assay has limited impact on MOA prediction accuracy, as monitoring only the nuclei stain could reach compelling levels of accuracy. This observation suggested that multiplexed measurements are highly correlated and nuclei stain could possibly reflect the general state of the cell. We then hypothesized that combining unrelated and possibly simple cell-based assays could bring a solution that would be biologically and technically more relevant to predict a drug target than using a single assay multiplexing dyes. We show that such a combination of past screen data could rationally be reused in screening facilities to train an ensemble classifier to predict drug targets and prioritize a possibly large list of unknown compound hits at once.

Biologically Relevant Heterogeneity: Metrics and Practical Insights

Heterogeneity is a fundamental property of biological systems at all scales that must be addressed in a wide range of biomedical applications, including basic biomedical research, drug discovery, diagnostics, and the implementation of precision medicine. There are a number of published approaches to characterizing heterogeneity in cells in vitro and in tissue sections. However, there are no generally accepted approaches for the detection and quantitation of heterogeneity that can be applied in a relatively high-throughput workflow. This review and perspective emphasizes the experimental methods that capture multiplexed cell-level data, as well as the need for standard metrics of the spatial, temporal, and population components of heterogeneity. A recommendation is made for the adoption of a set of three heterogeneity indices that can be implemented in any high-throughput workflow to optimize the decision-making process. In addition, a pairwise mutual information method is suggested as an approach to characterizing the spatial features of heterogeneity, especially in tissue-based imaging. Furthermore, metrics for temporal heterogeneity are in the early stages of development. Example studies indicate that the analysis of functional phenotypic heterogeneity can be exploited to guide decisions in the interpretation of biomedical experiments, drug discovery, diagnostics, and the design of optimal therapeutic strategies for individual patients.

Chemoattractant concentration-dependent tuning of ERK signaling dynamics in migrating neutrophils

The directed migration (chemotaxis) of neutrophils toward the bacterial peptide N-formyl-Met-Leu-Phe (fMLP) is a crucial process in immune defense against invading bacteria. While navigating through a gradient of increasing concentrations of fMLP, neutrophils and neutrophil-like HL-60 cells switch from exhibiting directional migration at low fMLP concentrations to exhibiting circuitous migration at high fMLP concentrations. The extracellular signal-regulated kinase (ERK) pathway is implicated in balancing this fMLP concentration-dependent switch in migration modes. We investigated the role and regulation of ERK signaling through single-cell analysis of neutrophil migration in response to different fMLP concentrations over time. We found that ERK exhibited gradated, rather than all-or-none, responses to fMLP concentration. Maximal ERK activation occurred in response to about 100 nM fMLP, and ERK inactivation was promoted by p38. Furthermore, we found that directional migration of neutrophils reached a maximal extent at about 100 nM fMLP and that ERK, but not p38, was required for neutrophil migration. Thus, our data suggest that, in chemotactic neutrophils responding to fMLP, ERK displays gradated activation and p38-dependent inhibition and that these ERK dynamics promote neutrophil migration.

Cell Painting, a high-content image-based assay for morphological profiling using multiplexed fluorescent dyes

In morphological profiling, quantitative data are extracted from microscopy images of cells to identify biologically relevant similarities and differences among samples based on these profiles. This protocol describes the design and execution of experiments using Cell Painting, which is a morphological profiling assay that multiplexes six fluorescent dyes, imaged in five channels, to reveal eight broadly relevant cellular components or organelles. Cells are plated in multiwell plates, perturbed with the treatments to be tested, stained, fixed, and imaged on a high-throughput microscope. Next, an automated image analysis software identifies individual cells and measures ∼1,500 morphological features (various measures of size, shape, texture, intensity, and so on) to produce a rich profile that is suitable for the detection of subtle phenotypes. Profiles of cell populations treated with different experimental perturbations can be compared to suit many goals, such as identifying the phenotypic impact of chemical or genetic perturbations, grouping compounds and/or genes into functional pathways, and identifying signatures of disease. Cell culture and image acquisition takes 2 weeks; feature extraction and data analysis take an additional 1-2 weeks.

Applications in image-based profiling of perturbations

A dramatic shift has occurred in how biologists use microscopy images. Whether experiments are small-scale or high-throughput, automatically quantifying biological properties in images is now widespread. We see yet another revolution under way: a transition towards using automated image analysis to not only identify phenotypes a biologist specifically seeks to measure ('screening') but also as an unbiased and sensitive tool to capture a wide variety of subtle features of cell (or organism) state ('profiling'). Mapping similarities among samples using image-based (morphological) profiling has tremendous potential to transform drug discovery, functional genomics, and basic biological research. Applications include target identification, lead hopping, library enrichment, functionally annotating genes/alleles, and identifying small molecule modulators of gene activity and disease-specific phenotypes.

A metric and workflow for quality control in the analysis of heterogeneity in phenotypic profiles and screens

Heterogeneity is well recognized as a common property of cellular systems that impacts biomedical research and the development of therapeutics and diagnostics. Several studies have shown that analysis of heterogeneity: gives insight into mechanisms of action of perturbagens; can be used to predict optimal combination therapies; and can be applied to tumors where heterogeneity is believed to be associated with adaptation and resistance. Cytometry methods including high content screening (HCS), high throughput microscopy, flow cytometry, mass spec imaging and digital pathology capture cell level data for populations of cells. However it is often assumed that the population response is normally distributed and therefore that the average adequately describes the results. A deeper understanding of the results of the measurements and more effective comparison of perturbagen effects requires analysis that takes into account the distribution of the measurements, i.e. the heterogeneity. However, the reproducibility of heterogeneous data collected on different days, and in different plates/slides has not previously been evaluated. Here we show that conventional assay quality metrics alone are not adequate for quality control of the heterogeneity in the data. To address this need, we demonstrate the use of the Kolmogorov-Smirnov statistic as a metric for monitoring the reproducibility of heterogeneity in an SAR screen, describe a workflow for quality control in heterogeneity analysis. One major challenge in high throughput biology is the evaluation and interpretation of heterogeneity in thousands of samples, such as compounds in a cell-based screen. In this study we also demonstrate that three heterogeneity indices previously reported, capture the shapes of the distributions and provide a means to filter and browse big data sets of cellular distributions in order to compare and identify distributions of interest. These metrics and methods are presented as a workflow for analysis of heterogeneity in large scale biology projects.

A Versatile Cell Death Screening Assay Using Dye-Stained Cells and Multivariate Image Analysis

A novel dye-based method for measuring cell death in image-based screens is presented. Unlike conventional high- and medium-throughput cell death assays that measure only one form of cell death accurately, using multivariate analysis of micrographs of cells stained with the inexpensive mix, red dye nonyl acridine orange, and a nuclear stain, it was possible to quantify cell death induced by a variety of different agonists even without a positive control. Surprisingly, using a single known cytotoxic agent as a positive control for training a multivariate classifier allowed accurate quantification of cytotoxicity for mechanistically unrelated compounds enabling generation of dose–response curves. Comparison with low throughput biochemical methods suggested that cell death was accurately distinguished from cell stress induced by low concentrations of the bioactive compounds Tunicamycin and Brefeldin A. High-throughput image-based format analyses of more than 300 kinase inhibitors correctly identified 11 as cytotoxic with only 1 false positive. The simplicity and robustness of this dye-based assay makes it particularly suited to live cell screening for toxic compounds.

Morphological Profiles of RNAi-Induced Gene Knockdown Are Highly Reproducible but Dominated by Seed Effects

RNA interference and morphological profiling-the measurement of thousands of phenotypes from individual cells by microscopy and image analysis-are a potentially powerful combination. We show that morphological profiles of RNAi-induced knockdown using the Cell Painting assay are in fact highly sensitive and reproducible. However, we find that the magnitude and prevalence of off-target effects via the RNAi seed-based mechanism make morphological profiles of RNAi reagents targeting the same gene look no more similar than reagents targeting different genes. Pairs of RNAi reagents that share the same seed sequence produce image-based profiles that are much more similar to each other than profiles from pairs designed to target the same gene, a phenomenon previously observed in small-scale gene-expression profiling experiments. Various strategies have been used to enrich on-target versus off-target effects in the context of RNAi screening where a narrow set of phenotypes are measured, mostly based on comparing multiple sequences targeting the same gene; however, new approaches will be needed to make RNAi morphological profiling (that is, comparing multi-dimensional phenotypes) viable. We have shared our raw data and computational pipelines to facilitate research.

[High content screening in chemical biology: overview and main challenges]

The last two decades have seen the development of high content screening (HCS) methodology and its adaptation for the evaluation of small molecules as drug candidates or their use as chemical tools for research purpose. HCS was initially set-up for the understanding of the mechanism of action of compounds by testing them on cell based-assays for pharmacological and toxicological studies. Since the last decade, the use of HCS has been extended to academic research laboratories and this technology has become the starting point for numerous projects aiming at the identification of molecular targets and cellular pathways for a given disease on which novel type of drugs could act. This screening approach relies on image capture of fluorescently labeled cells therefore generating a large amount of data that must be handled by appropriate automated image analysis methods and storage instrumentation. These latter in addition to the integration and data sharing are current challenges that HCS must still tackle.

Detection of cell aggregation and altered cell viability by automated label-free video microscopy: a promising alternative to endpoint viability assays in high-throughput screening

Automated phase-contrast video microscopy now makes it feasible to monitor a high-throughput (HT) screening experiment in a 384-well microtiter plate format by collecting one time-lapse video per well. Being a very cost-effective and label-free monitoring method, its potential as an alternative to cell viability assays was evaluated. Three simple morphology feature extraction and comparison algorithms were developed and implemented for analysis of differentially time-evolving morphologies (DTEMs) monitored in phase-contrast microscopy videos. The most promising layout, pixel histogram hierarchy comparison (PHHC), was able to detect several compounds that did not induce any significant change in cell viability, but made the cell population appear as spheroidal cell aggregates. According to recent reports, all these compounds seem to be involved in inhibition of platelet-derived growth factor receptor (PDGFR) signaling. Thus, automated quantification of DTEM (AQDTEM) holds strong promise as an alternative or complement to viability assays in HT in vitro screening of chemical compounds.

Single-cell and multivariate approaches in genetic perturbation screens

Large-scale genetic perturbation screens are a classical approach in biology and have been crucial for many discoveries. New technologies can now provide unbiased quantification of multiple molecular and phenotypic changes across tens of thousands of individual cells from large numbers of perturbed cell populations simultaneously. In this Review, we describe how these developments have enabled the discovery of new principles of intracellular and intercellular organization, novel interpretations of genetic perturbation effects and the inference of novel functional genetic interactions. These advances now allow more accurate and comprehensive analyses of gene function in cells using genetic perturbation screens.

On comparing heterogeneity across biomarkers

Microscopy reveals complex patterns of cellular heterogeneity that can be biologically informative. However, a limitation of microscopy is that only a small number of biomarkers can typically be monitored simultaneously. Thus, a natural question is whether additional biomarkers provide a deeper characterization of the distribution of cellular states in a population. How much information about a cell's phenotypic state in one biomarker is gained by knowing its state in another biomarker? Here, we describe a framework for comparing phenotypic states across biomarkers. Our approach overcomes the current limitation of microscopy by not requiring costaining biomarkers on the same cells; instead, we require staining of biomarkers (possibly separately) on a common collection of phenotypically diverse cell lines. We evaluate our approach on two image datasets: 33 oncogenically diverse lung cancer cell lines stained with 7 biomarkers, and 49 less diverse subclones of one lung cancer cell line stained with 12 biomarkers. We first validate our method by comparing it to the "gold standard" of costaining. We then apply our approach to all pairs of biomarkers and use it to identify biomarkers that yield similar patterns of heterogeneity. The results presented in this work suggest that many biomarkers provide redundant information about heterogeneity. Thus, our approach provides a practical guide for selecting independently informative biomarkers and, more generally, will yield insights into both the connectivity of biological networks and the complexity of the state space of biological systems.

Identifying and quantifying heterogeneity in high content analysis: application of heterogeneity indices to drug discovery

One of the greatest challenges in biomedical research, drug discovery and diagnostics is understanding how seemingly identical cells can respond differently to perturbagens including drugs for disease treatment. Although heterogeneity has become an accepted characteristic of a population of cells, in drug discovery it is not routinely evaluated or reported. The standard practice for cell-based, high content assays has been to assume a normal distribution and to report a well-to-well average value with a standard deviation. To address this important issue we sought to define a method that could be readily implemented to identify, quantify and characterize heterogeneity in cellular and small organism assays to guide decisions during drug discovery and experimental cell/tissue profiling. Our study revealed that heterogeneity can be effectively identified and quantified with three indices that indicate diversity, non-normality and percent outliers. The indices were evaluated using the induction and inhibition of STAT3 activation in five cell lines where the systems response including sample preparation and instrument performance were well characterized and controlled. These heterogeneity indices provide a standardized method that can easily be integrated into small and large scale screening or profiling projects to guide interpretation of the biology, as well as the development of therapeutics and diagnostics. Understanding the heterogeneity in the response to perturbagens will become a critical factor in designing strategies for the development of therapeutics including targeted polypharmacology.

Multiple-protein detections of single-cells reveal cell-cell heterogeneity in human cells

Cell population represents an intrinsically heterogeneous and stochastic system, in which individual cells often behave very differently in molecular contents, functions and even genotypes from the population average in response to uniform physiological stimuli. The traditional bulk cellular analysis often overlooks cellular heterogeneity and does not provide information on cell-cell variations. Single-cell measurements can reveal information obscured in population averages, and enable us to determine distributions rather than averaged properties within a cell population. The level of complexity, with numerous variables acting at the same time, requires multiparametric and dynamic investigation of a large number of single cells. Multiplexed study can provide quantitative correlations or inter-relationships among multiple cellular components and molecular markers within a protein network or family in biological processes. In this paper, we applied multiple fluorophore-conjugated primary antibodies to detect multiple proteins expressed on the same singe cells from a clonal population. To reveal cell-cell heterogeneity, we quantified the histograms of six proteins within a cell population as functions of TNF-α stimulation time. Then, we quantified noise and noise strength of these protein histograms as functions of TNF-α stimulation time. Thirdly, we quantified correlation coefficients of multiple proteins expressed on same single-cells as functions of TNF-α stimulation time. Above parameters demonstrated nonlinear relationships with TNF-α stimulation. Quantification of above parameters on independent cell subpopulations further reveals the cell-cell heterogeneity when exposed to identical environmental conditions. Such cellular heterogeneity will be useful to characterize the disease progression and disease diagnoses.

Quantitative protein localization signatures reveal an association between spatial and functional divergences of proteins

Protein subcellular localization is a major determinant of protein function. However, this important protein feature is often described in terms of discrete and qualitative categories of subcellular compartments, and therefore it has limited applications in quantitative protein function analyses. Here, we present Protein Localization Analysis and Search Tools (PLAST), an automated analysis framework for constructing and comparing quantitative signatures of protein subcellular localization patterns based on microscopy images. PLAST produces human-interpretable protein localization maps that quantitatively describe the similarities in the localization patterns of proteins and major subcellular compartments, without requiring manual assignment or supervised learning of these compartments. Using the budding yeast Saccharomyces cerevisiae as a model system, we show that PLAST is more accurate than existing, qualitative protein localization annotations in identifying known co-localized proteins. Furthermore, we demonstrate that PLAST can reveal protein localization-function relationships that are not obvious from these annotations. First, we identified proteins that have similar localization patterns and participate in closely-related biological processes, but do not necessarily form stable complexes with each other or localize at the same organelles. Second, we found an association between spatial and functional divergences of proteins during evolution. Surprisingly, as proteins with common ancestors evolve, they tend to develop more diverged subcellular localization patterns, but still occupy similar numbers of compartments. This suggests that divergence of protein localization might be more frequently due to the development of more specific localization patterns over ancestral compartments than the occupation of new compartments. PLAST enables systematic and quantitative analyses of protein localization-function relationships, and will be useful to elucidate protein functions and how these functions were acquired in cells from different organisms or species. A public web interface of PLAST is available at http://plast.bii.a-star.edu.sg.

Proteins are fundamental building blocks of cells. They perform a variety of biological functions, many of which are essential to the vitality and normal functioning of cells. Proteins have to be located at the appropriate regions inside a cell to perform their functions. Therefore, when proteins change their locations, they may acquire new or different functions. However, the relationships between the locations and functions of proteins are difficult to analyze because protein locations are often represented in distinct and manually-defined categories of subcellular regions. Many proteins have complex or subtle differences in their localization patterns that can be hardly represented by these categories. Here, we present an automated analysis tool for generating quantitative signatures of protein localization patterns without requiring manual or automated assignments of proteins into distinct categories. We show that our tool can identify proteins located at the same subcellular regions more accurately than existing categorization-based methods. Our tool allows comprehensive and more accurate studies of the relationships between protein localizations and functions. By knowing where proteins are located and how their locations were changed, we may discover their functions and better understand how they acquire these functions.

Parameterizing cell-to-cell regulatory heterogeneities via stochastic transcriptional profiles

Regulated changes in gene expression underlie many biological processes, but globally profiling cell-to-cell variations in transcriptional regulation is problematic when measuring single cells. Transcriptome-wide identification of regulatory heterogeneities can be robustly achieved by randomly collecting small numbers of cells followed by statistical analysis. However, this stochastic-profiling approach blurs out the expression states of the individual cells in each pooled sample. Here, we show that the underlying distribution of single-cell regulatory states can be deconvolved from stochastic-profiling data through maximum-likelihood inference. Guided by the mechanisms of transcriptional regulation, we formulated plausible mixture models for cell-to-cell regulatory heterogeneity and maximized the resulting likelihood functions to infer model parameters. Inferences were validated both computationally and experimentally for different mixture models, which included regulatory states for multicellular function that were occupied by as few as 1 in 40 cells of the population. Importantly, when the method was extended to programs of heterogeneously coexpressed transcripts, we found that population-level inferences were much more accurate with pooled samples than with one-cell samples when the extent of sampling was limited. Our deconvolution method provides a means to quantify the heterogeneous regulation of molecular states efficiently and gain a deeper understanding of the heterogeneous execution of cell decisions.

A novel phenotypic dissimilarity method for image-based high-throughput screens

Background

Discovering functional relationships of genes through cell-based phenotyping has become an important approach in functional genomics. High-throughput imaging offers the ability to quantitatively assess complex phenotypes after perturbation by RNA interference (RNAi). Such image-based high-throughput RNAi screening studies have facilitated the discovery of novel components of gene networks and their interactions. Images generated by automated microscopy are typically analyzed by extracting quantitative features of individual cells, resulting in large multidimensional data sets. Robust and sensitive methods to interpret these data sets and to derive biologically relevant information in a high-throughput and unbiased manner remain to be developed.

Results

Here we propose a new analysis method, PhenoDissim, which computes the phenotypic dissimilarity between cell populations via Support Vector Machine classification and cross validation. Applying this method to a kinome RNAi screening data set, we demonstrate that the proposed method shows a good replicate reproducibility, separation of controls and clustering quality, and we are able to identify siRNA phenotypes and discover potential functional links between genes.

Conclusions

PhenoDissim is a novel analysis method for image-based high-throughput screen, relying on two parameters which can be automatically optimized without a priori knowledge. PhenoDissim is freely available as an R package.

cellXpress: a fast and user-friendly software platform for profiling cellular phenotypes

BACKGROUND

High-throughput, image-based screens of cellular responses to genetic or chemical perturbations generate huge numbers of cell images. Automated analysis is required to quantify and compare the effects of these perturbations. However, few of the current freely-available bioimage analysis software tools are optimized for efficient handling of these images. Even fewer of them are designed to transform the phenotypic features measured from these images into discriminative profiles that can reveal biologically meaningful associations among the tested perturbations.

RESULTS

We present a fast and user-friendly software platform called "cellXpress" to segment cells, measure quantitative features of cellular phenotypes, construct discriminative profiles, and visualize the resulting cell masks and feature values. We have also developed a suite of library functions to load the extracted features for further customizable analysis and visualization under the R computing environment. We systematically compared the processing speed, cell segmentation accuracy, and phenotypic-profile clustering performance of cellXpress to other existing bioimage analysis software packages or algorithms. We found that cellXpress outperforms these existing tools on three different bioimage datasets. We estimate that cellXpress could finish processing a genome-wide gene knockdown image dataset in less than a day on a modern personal desktop computer.

CONCLUSIONS

The cellXpress platform is designed to make fast and efficient high-throughput phenotypic profiling more accessible to the wider biological research community. The cellXpress installation packages for 64-bit Windows and Linux, user manual, installation guide, and datasets used in this analysis can be downloaded freely from http://www.cellXpress.org.

Benchmarking of Multivariate Similarity Measures for High-Content Screening Fingerprints in Phenotypic Drug Discovery

High-content screening (HCS) is a powerful tool for drug discovery being capable of measuring cellular responses to chemical disturbance in a high-throughput manner. HCS provides an image-based readout of cellular phenotypes, including features such as shape, intensity, or texture in a highly multiplexed and quantitative manner. The corresponding feature vectors can be used to characterize phenotypes and are thus defined as HCS fingerprints. Systematic analyses of HCS fingerprints allow for objective computational comparisons of cellular responses. Such comparisons therefore facilitate the detection of different compounds with different phenotypic outcomes from high-throughput HCS campaigns. Feature selection methods and similarity measures, as a basis for phenotype identification and clustering, are critical for the quality of such computational analyses. We systematically evaluated 16 different similarity measures in combination with linear and nonlinear feature selection methods for their potential to capture biologically relevant image features. Nonlinear correlation-based similarity measures such as Kendall’s τ and Spearman’s ρ perform well in most evaluation scenarios, outperforming other frequently used metrics (such as the Euclidian distance). We also present four novel modifications of the connectivity map similarity that surpass the original version, in our experiments. This study provides a basis for generic phenotypic analysis in future HCS campaigns.

A Novel Organelle Map Framework for High-Content Cell Morphology Analysis in High Throughput

A screening procedure was developed that takes advantage of the cellular normalization by micropatterning and a novel quantitative organelle mapping approach that allows unbiased and automated cell morphology comparison using black-box statistical testing. Micropatterns of extracellular matrix proteins force cells to adopt a reproducible shape and distribution of intracellular compartments avoiding strong cell-to-cell variation that is a major limitation of classical culture conditions. To detect changes in cell morphology induced by compound treatment, fluorescently labeled intracellular structures from several tens of micropatterned cells were transformed into probabilistic density maps. Then, the similarity or difference between two given density maps was quantified using statistical testing that evaluates differences directly from the data without additional analysis or any subjective decision. The versatility of this organelle mapping approach for different magnifications and its performance for different cell shapes has been assessed. Density-based analysis detected changes in cell morphology due to compound treatment in a small-scale proof-of-principle screen demonstrating its compatibility with high-throughput screening. This novel tool for high-content and high-throughput cellular phenotyping can potentially be used for a wide range of applications from drug screening to careful characterization of cellular processes.

A simple high-content cell cycle assay reveals frequent discrepancies between cell number and ATP and MTS proliferation assays

In order to efficiently characterize both antiproliferative potency and mechanism of action of small molecules targeting the cell cycle, we developed a high-throughput image-based assay to determine cell number and cell cycle phase distribution. Using this we profiled the effects of experimental and approved anti-cancer agents with a range mechanisms of action on a set of cell lines, comparing direct cell counting versus two metabolism-based cell viability/proliferation assay formats, ATP-dependent bioluminescence, MTS (3-(4,5-dimethylthiazol-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium) reduction, and a whole-well DNA-binding dye fluorescence assay. We show that, depending on compound mechanisms of action, the metabolism-based proxy assays are frequently prone to 1) significant underestimation of compound potency and efficacy, and 2) non-monotonic dose-response curves due to concentration-dependent phenotypic ‘switching’. In particular, potency and efficacy of DNA synthesis-targeting agents such as gemcitabine and etoposide could be profoundly underestimated by ATP and MTS-reduction assays. In the same image-based assay we showed that drug-induced increases in ATP content were associated with increased cell size and proportionate increases in mitochondrial content and respiratory flux concomitant with cell cycle arrest. Therefore, differences in compound mechanism of action and cell line-specific responses can yield significantly misleading results when using ATP or tetrazolium-reduction assays as a proxy for cell number when screening compounds for antiproliferative activity or profiling panels of cell lines for drug sensitivity.

A computational platform for robotized fluorescence microscopy (I): high-content image-based cell-cycle analysis

Hardware automation and software development have allowed a dramatic increase of throughput in both acquisition and analysis of images by associating an optimized statistical significance with fluorescence microscopy. Despite the numerous common points between fluorescence microscopy and flow cytometry (FCM), the enormous amount of applications developed for the latter have found relatively low space among the modern high-resolution imaging techniques. With the aim to fulfill this gap, we developed a novel computational platform named A.M.I.CO. (Automated Microscopy for Image-Cytometry) for the quantitative analysis of images from widefield and confocal robotized microscopes. Thanks to the setting up of both staining protocols and analysis procedures, we were able to recapitulate many FCM assays. In particular, we focused on the measurement of DNA content and the reconstruction of cell-cycle profiles with optimal parameters. Standard automated microscopes were employed at the highest optical resolution (200 nm), and white-light sources made it possible to perform an efficient multiparameter analysis. DNA- and protein-content measurements were complemented with image-derived information on their intracellular spatial distribution. Notably, the developed tools create a direct link between image-analysis and acquisition. It is therefore possible to isolate target populations according to a definite quantitative profile, and to relocate physically them for diffraction-limited data acquisition. Thanks to its flexibility and analysis-driven acquisition, A.M.I.CO. can integrate flow, image-stream and laser-scanning cytometry analysis, providing high-resolution intracellular analysis with a previously unreached statistical relevance.

Concise review: a high-content screening approach to stem cell research and drug discovery

High-throughput screening (HTS) is a technology widely used for early stages of drug discovery in pharmaceutical and biotechnology industries. Recent hardware and software improvements have enabled HTS to be used in combination with subcellular resolution microscopy, resulting in cell image-based HTS, called high-content screening (HCS). HCS allows the acquisition of deeper knowledge at a single-cell level such that more complex biological systems can be studied in a high-throughput manner. The technique is particularly well-suited for stem cell research and drug discovery, which almost inevitably require single-cell resolutions for the detection of rare phenotypes in heterogeneous cultures. With growing availability of facilities, instruments, and reagent libraries, small-to-moderate scale HCS can now be carried out in regular academic labs. We envision that the HCS technique will play an increasing role in both basic mechanism study and early-stage drug discovery on stem cells. Here, we review the development of HCS technique and its past application on stem cells and discuss possible future developments.

Intensity quantile estimation and mapping--a novel algorithm for the correction of image non-uniformity bias in HCS data

MOTIVATION

Image non-uniformity (NU) refers to systematic, slowly varying spatial gradients in images that result in a bias that can affect all downstream image processing, quantification and statistical analysis steps. Image NU is poorly modeled in the field of high-content screening (HCS), however, such that current conventional correction algorithms may be either inappropriate for HCS or fail to take advantage of the information available in HCS image data.

RESULTS

A novel image NU bias correction algorithm, termed intensity quantile estimation and mapping (IQEM), is described. The algorithm estimates the full non-linear form of the image NU bias by mapping pixel intensities to a reference intensity quantile function. IQEM accounts for the variation in NU bias over broad cell intensity ranges and data acquisition times, both of which are characteristic of HCS image datasets. Validation of the method, using simulated and HCS microtubule polymerization screen images, is presented. Two requirements of IQEM are that the dataset consists of large numbers of images acquired under identical conditions and that cells are distributed with no within-image spatial preference.

AVAILABILITY AND IMPLEMENTATION

MATLAB function files are available at http://nadon-mugqic.mcgill.ca/.

Shedding light on filovirus infection with high-content imaging

Microscopy has been instrumental in the discovery and characterization of microorganisms. Major advances in high-throughput fluorescence microscopy and automated, high-content image analysis tools are paving the way to the systematic and quantitative study of the molecular properties of cellular systems, both at the population and at the single-cell level. High-Content Imaging (HCI) has been used to characterize host-virus interactions in genome-wide reverse genetic screens and to identify novel cellular factors implicated in the binding, entry, replication and egress of several pathogenic viruses. Here we present an overview of the most significant applications of HCI in the context of the cell biology of filovirus infection. HCI assays have been recently implemented to quantitatively study filoviruses in cell culture, employing either infectious viruses in a BSL-4 environment or surrogate genetic systems in a BSL-2 environment. These assays are becoming instrumental for small molecule and siRNA screens aimed at the discovery of both cellular therapeutic targets and of compounds with anti-viral properties. We discuss the current practical constraints limiting the implementation of high-throughput biology in a BSL-4 environment, and propose possible solutions to safely perform high-content, high-throughput filovirus infection assays. Finally, we discuss possible novel applications of HCI in the context of filovirus research with particular emphasis on the identification of possible cellular biomarkers of virus infection.

Closed-form density-based framework for automatic detection of cellular morphology changes

A primary method for studying cellular function is to examine cell morphology after a given manipulation. Fluorescent markers attached to proteins/intracellular structures of interest in conjunction with 3D fluorescent microscopy are frequently exploited for functional analysis. Despite the central role of morphology comparisons in cell biological approaches, few statistical tools are available that allow biological scientists without a high level of statistical training to quantify the similarity or difference of fluorescent images containing multifactorial information. We transform intracellular structures into kernels and develop a multivariate two-sample test that is nonparametric and asymptotically normal to directly and quantitatively compare cellular morphologies. The asymptotic normality bypasses the computationally intensive calculations used by the usual resampling techniques to compute the P-value. Because all parameters required for the statistical test are estimated directly from the data, it does not require any subjective decisions. Thus, we provide a black-box method for unbiased, automated comparison of cell morphology. We validate the performance of our test statistic for finite synthetic samples and experimental data. Employing our test for the comparison of the morphology of intracellular multivesicular bodies, we detect changes in their distribution after disruption of the cellular microtubule cytoskeleton with high statistical significance in fixed samples and live cell analysis. These results demonstrate that density-based comparison of multivariate image information is a powerful tool for automated detection of cell morphology changes. Moreover, the underlying mathematics of our test statistic is a general technique, which can be applied in situations where two data samples are compared.

High-content screening data management for drug discovery in a small- to medium-size laboratory: results of a collaborative pilot study focused on user expectations as indicators of effectiveness

High-content screening (HCS) technology provides a powerful vantage point to approach biological problems; it allows analysis of cell parameters, including changes in cell or protein movement, shape, or texture. As part of a collaborative pilot research project to improve bioscience research data integration, we identified HCS data management as an area ripe for advancement. A primary goal was to develop an integrated data management and analysis system suitable for small- to medium-size HCS programs that would improve research productivity and increase work satisfaction. A system was developed that uses Labmatrix, a Web-based research data management platform, to integrate and query data derived from a Cellomics STORE database. Focusing on user expectations, several barriers to HCS productivity were identified and reduced or eliminated. The impact of the project on HCS research productivity was tested through a series of 18 lab-requested integrated data queries, 7 of which were fully enabled, 7 partially enabled, and 4 enabled through data export to standalone data analysis tools. The results are limited to one laboratory, but this pilot suggests that through an "implementation research" approach, a network of small- to medium-size laboratories involved in HCS projects could achieve greater productivity and satisfaction in drug discovery research.

Extracting fluorescent reporter time courses of cell lineages from high-throughput microscopy at low temporal resolution

The extraction of fluorescence time course data is a major bottleneck in high-throughput live-cell microscopy. Here we present an extendible framework based on the open-source image analysis software ImageJ, which aims in particular at analyzing the expression of fluorescent reporters through cell divisions. The ability to track individual cell lineages is essential for the analysis of gene regulatory factors involved in the control of cell fate and identity decisions. In our approach, cell nuclei are identified using Hoechst, and a characteristic drop in Hoechst fluorescence helps to detect dividing cells. We first compare the efficiency and accuracy of different segmentation methods and then present a statistical scoring algorithm for cell tracking, which draws on the combination of various features, such as nuclear intensity, area or shape, and importantly, dynamic changes thereof. Principal component analysis is used to determine the most significant features, and a global parameter search is performed to determine the weighting of individual features. Our algorithm has been optimized to cope with large cell movements, and we were able to semi-automatically extract cell trajectories across three cell generations. Based on the MTrackJ plugin for ImageJ, we have developed tools to efficiently validate tracks and manually correct them by connecting broken trajectories and reassigning falsely connected cell positions. A gold standard consisting of two time-series with 15,000 validated positions will be released as a valuable resource for benchmarking. We demonstrate how our method can be applied to analyze fluorescence distributions generated from mouse stem cells transfected with reporter constructs containing transcriptional control elements of the Msx1 gene, a regulator of pluripotency, in mother and daughter cells. Furthermore, we show by tracking zebrafish PAC2 cells expressing FUCCI cell cycle markers, our framework can be easily adapted to different cell types and fluorescent markers.

Automated microscopy and image analysis for androgen receptor function

Systems-level approaches have emerged that rely on analytical, microscopy-based technology for the discovery of novel drug targets and the mechanisms driving AR signaling, transcriptional activity, and ligand independence. Single cell behavior can be quantified by high-throughput microscopy methods through analysis of endogenous protein levels and localization or creation of biosensor cell lines that can simultaneously detect both acute and latent responses to known and unknown androgenic stimuli. The cell imaging and analytical protocols can be automated to discover agonist/antagonist response windows for nuclear translocation, reporter gene activity, nuclear export, and subnuclear transcription events, facilitating access to a multiplex model system that is inherently unavailable through classic biochemical approaches. In this chapter, we highlight the key steps needed for developing, conducting, and analyzing high-throughput screens to identify effectors of AR signaling.

Sequential array cytometry: multi-parameter imaging with a single fluorescent channel

Heterogeneity within the human population and within diseased tissues necessitates a personalized medicine approach to diagnostics and the treatment of diseases. Functional assays at the single-cell level can contribute to uncovering heterogeneity and ultimately assist in improved treatment decisions based on the presence of outlier cells. We aim to develop a platform for high-throughput, single-cell-based assays using well-characterized hydrodynamic cell isolation arrays which allow for precise cell and fluid handling. Here, we demonstrate the ability to extract spatial and temporal information about several intracellular components using a single fluorescent channel, eliminating the problem of overlapping fluorescence emission spectra. Integrated with imaging technologies such as wide field-of-view lens-free fluorescent imaging, fiber-optic array scanning technology, and microlens arrays, use of a single fluorescent channel will reduce the cost of reagents and optical components. Specifically, we sequentially stain hydrodynamically trapped cells with three biochemical labels all sharing the same fluorescence excitation and emission spectrum. These markers allow us to analyze the amount of DNA, and compare nucleus-to-cytoplasm ratio, as well as glycosylation of surface proteins. By imaging cells in real-time we enable measurements of temporal localization of cellular components and intracellular reaction kinetics, the latter is used as a measurement of multi-drug resistance. Demonstrating the efficacy of this single-cell analysis platform is the first step in designing and implementing more complete assays, aimed toward improving diagnosis and personalized treatments to complex diseases.

Patterns of basal signaling heterogeneity can distinguish cellular populations with different drug sensitivities

Phenotypic heterogeneity has been widely observed in cellular populations. However, the extent to which heterogeneity contains biologically or clinically important information is not well understood. Here, we investigated whether patterns of basal signaling heterogeneity, in untreated cancer cell populations, could distinguish cellular populations with different drug sensitivities. We modeled cellular heterogeneity as a mixture of stereotyped signaling states, identified based on colocalization patterns of activated signaling molecules from microscopy images. We found that patterns of heterogeneity could be used to separate the most sensitive and resistant populations to paclitaxel within a set of H460 lung cancer clones and within the NCI-60 panel of cancer cell lines, but not for a set of less heterogeneous, immortalized noncancer human bronchial epithelial cell (HBEC) clones. Our results suggest that patterns of signaling heterogeneity, characterized as ensembles of a small number of distinct phenotypic states, can reveal functional differences among cellular populations.

A Personal Perspective on High-Content Screening (HCS)

High-content screening (HCS) was introduced in 1997 based on light microscope imaging technologies to address the need for an automated platform that could analyze large numbers of individual cells with subcellular resolution using standard microplates. Molecular specificity based on fluorescence was a central element of the platform taking advantage of the growing list of reagent classes and the ability to multiplex. In addition, image analysis coupled to data management, data mining, and data visualization created a tool that focused on biological information and knowledge to begin exploring the functions of genes identified in the genomics revolution. This overview looks at the development of HCS, the evolution of the technologies, and the market up to the present day. In addition, the options for adopting uniform definitions is suggested along with a perspective on what advances are needed to continue building the value of HCS in biomedical research, drug discovery, and development and diagnostics.

Cellular heterogeneity: do differences make a difference?

A central challenge of biology is to understand how individual cells process information and respond to perturbations. Much of our knowledge is based on ensemble measurements. However, cell-to-cell differences are always present to some degree in any cell population, and the ensemble behaviors of a population may not represent the behaviors of any individual cell. Here, we discuss examples of when heterogeneity cannot be ignored and describe practical strategies for analyzing and interpreting cellular heterogeneity.