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

Systematic data analysis pipeline for quantitative morphological cell phenotyping

Quantitative morphological phenotyping (QMP) is an image-based method used to capture morphological features at both the cellular and population level. Its interdisciplinary nature, spanning from data collection to result analysis and interpretation, can lead to uncertainties, particularly among those new to this actively growing field. High analytical specificity for a typical QMP is achieved through sophisticated approaches that can leverage subtle cellular morphological changes. Here, we outline a systematic workflow to refine the QMP methodology. For a practical review, we describe the main steps of a typical QMP; in each step, we discuss the available methods, their applications, advantages, and disadvantages, along with the R functions and packages for easy implementation. This review does not cover theoretical backgrounds, but provides several references for interested researchers. It aims to broaden the horizons for future phenome studies and demonstrate how to exploit years of endeavors to achieve more with less.

Targeting therapy-persistent residual disease

Disease relapse driven by acquired drug resistance limits the effectiveness of most systemic anti-cancer agents. Targeting persistent cancer cells in residual disease before relapse has emerged as a potential strategy for enhancing the efficacy and the durability of current therapies. However, barriers remain to implementing persister-directed approaches in the clinic. This Perspective discusses current preclinical and clinical complexities and outlines key steps toward the development of clinical strategies that target therapy-persistent residual disease.

A microfluidic array enabling generation of identical biochemical stimulating signals to trapped biological cells for single-cell dynamics

Microfluidic-based analyses of single-cell dynamics in response to dynamic biochemical signals are emerging as pivotal approaches for investigating the effects of extracellular microenvironmental biochemical factors on cellular structure, function, and behavior. However, current devices often fail to consistently apply identical dynamic biochemical signals to trapped cells. In this study, we introduce a novel radially distributed single-cell trapping microfluidic array, designed to quantitatively and consistently apply identical biochemical stimulating signals to each trapped cell. Numerical simulations were employed to optimize microchannel geometry, enhancing trapping efficiency while minimizing signal distortion. Experimental validation demonstrated the trapping success rate and the single-cell trapping efficiency exceeding 99% and 85%, respectively. The microarray's capability to deliver identical dynamic biochemical stimulating signals, with various waveforms, to each unit was confirmed through fluorescein transport tests. Furthermore, we examined the intracellular calcium dynamics of U-2 OS human osteosarcoma cells in response to dynamic ATP signals, observing both single-peak calcium responses and calcium oscillations, which were modelled by a second-order system with a natural frequency of 1.6 mHz. Overall, our proposed microfluidic array offers a robust and valuable framework for advancing the understanding of single-cell dynamics.

Extended methods for spatial cell classification with DBSCAN-CellX

Local cell densities and positioning within cellular monolayers and stratified epithelia have important implications for cell interactions and the functionality of various biological processes. To analyze the relationship between cell localization and tissue physiology, density-based clustering algorithms, such as DBSCAN, allow for a detailed characterization of the spatial distribution and positioning of individual cells. However, these methods rely on predefined parameters that influence the outcome of the analysis. With varying cell densities in cell cultures or tissues impacting cell sizes and, thus, cellular proximities, these parameters need to be carefully chosen. In addition, standard DBSCAN approaches generally come short in appropriately identifying individual cell positions. We therefore developed three extensions to the standard DBSCAN-algorithm that provide: (i) an automated parameter identification to reliably identify cell clusters, (ii) an improved identification of cluster edges; and (iii) an improved characterization of the relative positioning of cells within clusters. We apply our novel methods, which are provided as a user-friendly OpenSource-software package (DBSCAN-CellX), to cellular monolayers of different cell lines. Thereby, we show the importance of the developed extensions for the appropriate analysis of cell culture experiments to determine the relationship between cell localization and tissue physiology.

High throughput microscopy and single cell phenotypic image-based analysis in toxicology and drug discovery

Measuring single cell responses to the universe of chemicals (drugs, natural products, environmental toxicants etc.) is of paramount importance to human health as phenotypic variability in sensing stimuli is a hallmark of biology that is considered during high throughput screening. One of the ways to approach this problem is via high throughput, microscopy-based assays coupled with multi-dimensional single cell analysis methods. Here, we will summarize some of the efforts in this vast and growing field, focusing on phenotypic screens (e.g., Cell Painting), single cell analytics and quality control, with particular attention to environmental toxicology and drug screening. We will discuss advantages and limitations of high throughput assays with various end points and levels of complexity.

Embryonic stem cell ERK, AKT, plus STAT3 response dynamics combinatorics are heterogeneous but NANOG state independent

Signaling is central in cell fate regulation, and relevant information is encoded in its activity over time (i.e., dynamics). However, simultaneous dynamics quantification of several pathways in single mammalian stem cells has not yet been accomplished. Here we generate mouse embryonic stem cell (ESC) lines simultaneously expressing fluorescent reporters for ERK, AKT, and STAT3 signaling activity, which all control pluripotency. We quantify their single-cell dynamics combinations in response to different self-renewal stimuli and find striking heterogeneity for all pathways, some dependent on cell cycle but not pluripotency states, even in ESC populations currently assumed to be highly homogeneous. Pathways are mostly independently regulated, but some context-dependent correlations exist. These quantifications reveal surprising single-cell heterogeneity in the important cell fate control layer of signaling dynamics combinations and raise fundamental questions about the role of signaling in (stem) cell fate control.

Quantification of long-term doxorubicin response dynamics in breast cancer cell lines to direct treatment schedules

While acquired chemoresistance is recognized as a key challenge to treating many types of cancer, the dynamics with which drug sensitivity changes after exposure are poorly characterized. Most chemotherapeutic regimens call for repeated dosing at regular intervals, and if drug sensitivity changes on a similar time scale then the treatment interval could be optimized to improve treatment performance. Theoretical work suggests that such optimal schedules exist, but experimental confirmation has been obstructed by the difficulty of deconvolving the simultaneous processes of death, adaptation, and regrowth taking place in cancer cell populations. Here we present a method of optimizing drug schedules in vitro through iterative application of experimentally calibrated models, and demonstrate its ability to characterize dynamic changes in sensitivity to the chemotherapeutic doxorubicin in three breast cancer cell lines subjected to treatment schedules varying in concentration, interval between pulse treatments, and number of sequential pulse treatments. Cell populations are monitored longitudinally through automated imaging for 600–800 hours, and this data is used to calibrate a family of cancer growth models, each consisting of a system of ordinary differential equations, derived from the bi-exponential model which characterizes resistant and sensitive subpopulations. We identify a model incorporating both a period of growth arrest in surviving cells and a delay in the death of chemosensitive cells which outperforms the original bi-exponential growth model in Akaike Information Criterion based model selection, and use the calibrated model to quantify the performance of each drug schedule. We find that the inter-treatment interval is a key variable in determining the performance of sequential dosing schedules and identify an optimal retreatment time for each cell line which extends regrowth time by 40%-239%, demonstrating that the time scale of changes in chemosensitivity following doxorubicin exposure allows optimization of drug scheduling by varying this inter-treatment interval.

Acquired chemoresistance is a common cause of treatment failure in cancer. The scheduling of a multi-dose course of chemotherapeutic treatment may influence the dynamics of acquired chemoresistance, and drug schedule optimization may increase the duration of effectiveness of a particular chemotherapeutic agent for a particular patient. Here we present a method for experimentally optimizing an in vitro drug schedule through iterative rounds of experimentation and computational analysis, and demonstrate the method’s ability to improve the performance of doxorubicin treatment in three breast carcinoma cell lines. Specifically, we find that the interval between drug exposures can be optimized while holding drug concentration and number of treatments constant, suggesting that this may be a key variable to explore in future drug schedule optimization efforts. We further use this method’s model calibration and selection process to extract information about the underlying biology of the doxorubicin response, and find that the incorporation of delays on both cell death and regrowth are necessary for accurate parameterization of cell growth data.

Quality Control for Single Cell Imaging Analytics Using Endocrine Disruptor-Induced Changes in Estrogen Receptor Expression

BACKGROUND

Diverse toxicants and mixtures that affect hormone responsive cells [endocrine disrupting chemicals (EDCs)] are highly pervasive in the environment and are directly linked to human disease. They often target the nuclear receptor family of transcription factors modulating their levels and activity. Many high-throughput assays have been developed to query such toxicants; however, single-cell analysis of EDC effects on endogenous receptors has been missing, in part due to the lack of quality control metrics to reproducibly measure cell-to-cell variability in responses.

OBJECTIVE

We began by developing single-cell imaging and informatic workflows to query whether the single cell distribution of the estrogen receptor-α (ER), used as a model system, can be used to measure effects of EDCs in a sensitive and reproducible manner.

METHODS

We used high-throughput microscopy, coupled with image analytics to measure changes in single cell ER nuclear levels on treatment with ∼100 toxicants, over a large number of biological and technical replicates.

RESULTS

We developed a two-tiered quality control pipeline for single cell analysis and tested it against a large set of biological replicates, and toxicants from the EPA and Agency for Toxic Substances and Disease Registry lists. We also identified a subset of potentially novel EDCs that were active only on the endogenous ER level and activity as measured by single molecule RNA fluorescence in situ hybridization (RNA FISH).

DISCUSSION

We demonstrated that the distribution of ER levels per cell, and the changes upon chemical challenges were remarkably stable features; and importantly, these features could be used for quality control and identification of endocrine disruptor toxicants with high sensitivity. When coupled with orthogonal assays, ER single cell distribution is a valuable resource for high-throughput screening of environmental toxicants. https://doi.org/10.1289/EHP9297.

Live-Cell Imaging Shows Uneven Segregation of Extrachromosomal DNA Elements and Transcriptionally Active Extrachromosomal DNA Hubs in Cancer

Oncogenic extrachromosomal DNA elements (ecDNA) play an important role in tumor evolution, but our understanding of ecDNA biology is limited. We determined the distribution of single-cell ecDNA copy number across patient tissues and cell line models and observed how cell-to-cell ecDNA frequency varies greatly. The exceptional intratumoral heterogeneity of ecDNA suggested ecDNA-specific replication and propagation mechanisms. To evaluate the transfer of ecDNA genetic material from parental to offspring cells during mitosis, we established the CRISPR-based ecTag method. ecTag leverages ecDNA-specific breakpoint sequences to tag ecDNA with fluorescent markers in living cells. Applying ecTag during mitosis revealed disjointed ecDNA inheritance patterns, enabling rapid ecDNA accumulation in individual cells. After mitosis, ecDNAs clustered into ecDNA hubs, and ecDNA hubs colocalized with RNA polymerase II, promoting transcription of cargo oncogenes. Our observations provide direct evidence for uneven segregation of ecDNA and shed new light on mechanisms through which ecDNAs contribute to oncogenesis. SIGNIFICANCE: ecDNAs are vehicles for oncogene amplification. The circular nature of ecDNA affords unique properties, such as mobility and ecDNA-specific replication and segregation behavior. We uncovered fundamental ecDNA properties by tracking ecDNAs in live cells, highlighting uneven and random segregation and ecDNA hubs that drive cargo gene transcription.See related commentary by Henssen, p. 293.This article is highlighted in the In This Issue feature, p. 275.

Multi-omics analysis defines 5-fluorouracil drug resistance in 3D HeLa carcinoma cell model

Cervical cancer is a serious health problem in women around the globe. However, the use of clinical drug is seriously dampened by the development of drug resistance. Efficient in vitro tumor model is essential to improve the efficiency of drug screening and the accuracy of clinical application. Multicellular tumor spheroids (MTSs) can in a way recapitulates tumor traits in vivo, thereby representing a powerful transitional model between 2D monolayer culture and xenograft. In this study, based on the liquid overlay method, a protocol for rapid generation of the MTSs with uniform size and high reproducibility in a high-throughput manner was established. As expected, the cytotoxicity results showed that there was enhanced 5-fluorouracil (5-FU) resistance of HeLa carcinoma cells in 3D MTSs than 2D monolayer culture with a resistance index of 5.72. In order to obtain a holistic view of the molecular mechanisms that drive 5-FU resistance in 3D HeLa carcinoma cells, a multi-omics study was applied to discover hidden biological regularities. It was observed that in the 3D MTSs mitochondrial function-related proteins and the metabolites of the tricarboxylic acid cycle (TCA cycle) were significantly decreased, and the cellular metabolism was shifted towards glycolysis. The differences in the protein synthesis, processing, and transportation between 2D monolayer cultures and 3D MTSs were significant, mainly in the heat shock protein family, with the up-regulation of protein folding function in endoplasmic reticulum (ER), which promoted the maintenance of ER homeostasis in the 3D MTSs. In addition, at the transcript and protein level, the expression of extracellular matrix (ECM) proteins (e.g., laminin and collagen) were up-regulated in the 3D MTSs, which enhanced the physical barrier of drug penetration. Summarizing, this study formulates a rapid, scalable and reproducible in vitro model of 3D MTS for drug screening purposes, and the findings establish a critical role of glycolytic metabolism, ER hemostasis and ECM proteins expression profiling in tumor chemoresistance of HeLa carcinoma cells towards 5-FU.

High-Content Single-Cell Förster Resonance Energy Transfer Imaging of Cultured Striatal Neurons Reveals Novel Cross-Talk in the Regulation of Nuclear Signaling by Protein Kinase A and Extracellular Signal-Regulated Kinase 1/2

Genetically encoded biosensors can be used to track signaling events in living cells by measuring changes in fluorescence emitted by one or more fluorescent proteins. Here, we describe the use of genetically encoded biosensors based on Förster resonance energy transfer (FRET), combined with high-content microscopy, to image dynamic signaling events simultaneously in thousands of neurons in response to drug treatments. We first applied this approach to examine intercellular variation in signaling responses among cultured striatal neurons stimulated with multiple drugs. Using high-content FRET imaging and immunofluorescence, we identified neuronal subpopulations with unique responses to pharmacological manipulation and used nuclear morphology to identify medium spiny neurons within these heterogeneous striatal cultures. Focusing on protein kinase A (PKA) and extracellular signal-regulated kinase 1/2 (ERK1/2) signaling in the cytoplasm and nucleus, we noted pronounced intercellular differences among putative medium spiny neurons, in both the magnitude and kinetics of signaling responses to drug application. Importantly, a conventional "bulk" analysis that pooled all cells in culture yielded a different rank order of drug potency than that revealed by single-cell analysis. Using a single-cell analytical approach, we dissected the relative contributions of PKA and ERK1/2 signaling in striatal neurons and unexpectedly identified a novel role for ERK1/2 in promoting nuclear activation of PKA in striatal neurons. This finding adds a new dimension of signaling crosstalk between PKA and ERK1/2 with relevance to dopamine D1 receptor signaling in striatal neurons. In conclusion, high-content single-cell imaging can complement and extend traditional population-level analyses and provides a novel vantage point from which to study cellular signaling. SIGNIFICANCE STATEMENT: High-content imaging revealed substantial intercellular variation in the magnitude and pattern of intracellular signaling events driven by receptor stimulation. Since individual neurons within the same population can respond differently to a given agonist, interpreting measures of intracellular signaling derived from the averaged response of entire neuronal populations may not always reflect what happened at the single-cell level. This study uses this approach to identify a new form of cross-talk between PKA and ERK1/2 signaling in the nucleus of striatal neurons.

A Microfluidic Array Device for Single Cell Capture and Intracellular Ca2+ Response Analysis Induced by Dynamic Biochemical Stimulus

A microfluidic array was constructed for trapping single cell and loading identical dynamic biochemical stimulation for gain a better understanding of Ca2+ signalling in single cells by applying extracellular dynamic biochemical stimulus. This microfluidic array consists of multiple radially aligned flow channels with equal intersection angles, which was designed by a combination of stagnation point flow and physical barrier. Numerical simulation results and trajectory analysis shown the effectiveness of this single cell trapping device. Fluorescent experiment results demonstrated the effects of flow rate and frequency of dynamic stimulus on the profiles of biochemical concentration which exposed on captured cells. In this array chip, the captured single cells in each trapping channels were able to receive identical extracellular dynamic biochemical stimuli which being transmitted from the entrance at the middle of the microfluidic array. Besides, after loading dynamic Adenosine Triphosphate (ATP) stimulation on captured cells by this device, consistent average intracellular Ca2+ dynamics phase and cellular heterogeneity were observed in captured single K562 cells. Furthermore, this device is able to be used for investigating cellular respond in single cells to temporally varying environments by modulating the stimulation signal in terms of concentration, pattern, and duration of exposure.

Cell-to-cell variation in gene expression and the aging process

There is tremendous variation in biological traits, and much of it is not accounted for by variation in DNA sequence, including human diseases and lifespan. Emerging evidence points to differences in the execution of the genetic program as a key source of variation, be it stochastic variation or programmed variation. Here we discuss variation in gene expression as an intrinsic property and how it could contribute to variation in traits, including the rate of aging. The review is divided into sections describing the historical context and evidence to date for nongenetic variation, the different approaches that may be used to detect nongenetic variation, and recent findings showing that the amount of variation in gene expression can be both genetically programmed and epigenetically controlled. Finally, we present evidence that changes in cell-to-cell variation in gene expression emerge as part of the aging process and may be linked to disease vulnerability as a function of age. These emerging concepts are likely to be important across the spectrum of biomedical research and may well underpin what we understand as biological aging.

Lactic Acidosis Interferes With Toxicity of Perifosine to Colorectal Cancer Spheroids: Multimodal Imaging Analysis

Colorectal cancer (CRC) is a disease with constantly increasing incidence and high mortality. The treatment efficacy could be curtailed by drug resistance resulting from poor drug penetration into tumor tissue and the tumor-specific microenvironment, such as hypoxia and acidosis. Furthermore, CRC tumors can be exposed to different pH depending on the position in the intestinal tract. CRC tumors often share upregulation of the Akt signaling pathway. In this study, we investigated the role of external pH in control of cytotoxicity of perifosine, the Akt signaling pathway inhibitor, to CRC cells using 2D and 3D tumor models. In 3D settings, we employed an innovative strategy for simultaneous detection of spatial drug distribution and biological markers of proliferation/apoptosis using a combination of mass spectrometry imaging and immunohistochemistry. In 3D conditions, low and heterogeneous penetration of perifosine into the inner parts of the spheroids was observed. The depth of penetration depended on the treatment duration but not on the external pH. However, pH alteration in the tumor microenvironment affected the distribution of proliferation- and apoptosis-specific markers in the perifosine-treated spheroid. Accurate co-registration of perifosine distribution and biological response in the same spheroid section revealed dynamic changes in apoptotic and proliferative markers occurring not only in the perifosine-exposed cells, but also in the perifosine-free regions. Cytotoxicity of perifosine to both 2D and 3D cultures decreased in an acidic environment below pH 6.7. External pH affects cytotoxicity of the other Akt inhibitor, MK-2206, in a similar way. Our innovative approach for accurate determination of drug efficiency in 3D tumor tissue revealed that cytotoxicity of Akt inhibitors to CRC cells is strongly dependent on pH of the tumor microenvironment. Therefore, the effect of pH should be considered during the design and pre-clinical/clinical testing of the Akt-targeted cancer therapy.

Poikilosis - pervasive biological variation

Biological systems are dynamic and display heterogeneity at all levels. Ubiquitous heterogeneity, here called for poikilosis, is an integral and important property of organisms and in molecules, systems and processes within them. Traditionally, heterogeneity in biology and experiments has been considered as unwanted noise, here poikilosis is shown to be the normal state. Acceptable variation ranges are called as lagom. Non-lagom, variations that are too extensive, have negative effects, which influence interconnected levels and once the variation is large enough cause a disease and can lead even to death. Poikilosis has numerous applications and consequences e.g. for how to design, analyze and report experiments, how to develop and apply prediction and modelling methods, and in diagnosis and treatment of diseases. Poikilosis-aware new and practical definitions are provided for life, death, senescence, disease, and lagom. Poikilosis is the first new unifying theory in biology since evolution and should be considered in every scientific study.

Poikilosis - pervasive biological variation

Biological systems are dynamic and display heterogeneity at all levels. Ubiquitous heterogeneity, here called for poikilosis, is an integral and important property of organisms and in molecules, systems and processes within them. Traditionally, heterogeneity in biology and experiments has been considered as unwanted noise, here poikilosis is shown to be the normal state. Acceptable variation ranges are called as lagom. Non-lagom, variations that are too extensive, have negative effects, which influence interconnected levels and once the variation is large enough cause a disease and can lead even to death. Poikilosis has numerous applications and consequences e.g. for how to design, analyze and report experiments, how to develop and apply prediction and modelling methods, and in diagnosis and treatment of diseases. Poikilosis-aware new and practical definitions are provided for life, death, senescence, disease, and lagom. Poikilosis is the first new unifying theory in biology since evolution and should be considered in every scientific study.

Identification of an early cell fate regulator by detecting dynamics in transcriptional heterogeneity and co-regulation during astrocyte differentiation

There are an increasing number of reports that characterize the temporal behavior of gene expression at the single-cell level during cell differentiation. Despite accumulation of data describing the heterogeneity of biological responses, the dynamics of gene expression heterogeneity and its regulation during the differentiation process have not been studied systematically. To understand transcriptional heterogeneity during astrocyte differentiation, we analyzed single-cell transcriptional data from cells representing the different stages of astrocyte differentiation. When we compared the transcriptional variability of co-expressed genes between the undifferentiated and differentiated states, we found that there was significant increase in transcriptional variability in the undifferentiated state. The genes showing large changes in both "variability" and "correlation" between neural stem cells (NSCs) and astrocytes were found to be functionally involved in astrocyte differentiation. We determined that these genes are potentially regulated by Ascl1, a previously known oscillatory gene in NSCs. Pharmacological blockade of Ntsr2, which is transcriptionally co-regulated with Ascl1, showed that Ntsr2 may play an important role in the differentiation from NSCs to astrocytes. This study shows the importance of characterizing transcriptional heterogeneity and rearrangement of the co-regulation network between different cell states. It also highlights the potential for identifying novel regulators of cell differentiation that will further increase our understanding of the molecular mechanisms underlying the differentiation process.

Intratumor heterogeneity index of breast carcinomas based on DNA methylation profiles

Background

Cancer cells evolve and constitute heterogeneous populations that fluctuate in space and time and are subjected to selection generating intratumor heterogeneity. This phenomenon is determined by the acquisition of genetic/epigenetic alterations and their selection over time which has clinical implications on drug resistance.

Methods

DNA extracted from different tumor cell populations (breast carcinomas, cancer cell lines and cellular clones) were analyzed by MS-MLPA. Methylation profiles were used to generate a heterogeneity index to quantify the magnitude of epigenetic heterogeneity in these populations. Cellular clones were obtained from single cells derived of MDA-MB 231 cancer cell lines applying serial limiting dilution method and morphology was analyzed by optical microscopy and flow cytometry. Clones characteristics were examined through cellular proliferation, migration capacity and apoptosis. Heterogeneity index was also calculated from beta values derived from methylation profiles of TCGA tumors.

Results

The study of methylation profiles of 23 fresh breast carcinomas revealed heterogeneous allele populations in these tumor pieces. With the purpose to measure the magnitude of epigenetic heterogeneity, we developed an heterogeneity index based on methylation information and observed that all tumors present their own heterogeneity level. Applying the index calculation in pure cancer cell populations such as cancer cell lines (MDA-MB 231, MCF-7, T47D, HeLa and K-562), we also observed epigenetic heterogeneity. In addition, we detected that clones obtained from the MDA-MB 231 cancer cell line generated their own new heterogeneity over time. Using TCGA tumors, we determined that the heterogeneity index correlated with prognostic and predictive factors like tumor size (p = 0.0088), number of affected axillary nodes (p = 0.007), estrogen receptor expression (p < 0.0001) and HER2 positivity (p = 0.0007). When we analyzed molecular subtypes we found that they presented different heterogeneity levels. Interestingly, we also observed that all mentioned tumor cell populations shared a similar Heterogeneity index (HI) mean.

Conclusions

Our results show that each tumor presents a unique epigenetic heterogeneity level, which is associated with prognostic and predictive factors. We also observe that breast tumor subtypes differ in terms of epigenetic heterogeneity, which could serve as a new contribution to understand the different prognosis of these groups.

High-throughput microfluidic single-cell trapping arrays for biomolecular and imaging analysis

Single-cell analysis is of critical importance in revealing population heterogeneity, identifying minority sub-populations of interest, as well as discovering unique characteristics of individual cells. Microfluidic platforms work at the scale comparable to cell diameter and is suitable for single-cell manipulation. Here we present a microfluidic trapping array which is able to rapidly and deterministically trap single-cells in highly-packed microwells. This chapter first describes the design and fabrication protocols of the trapping array, and then presents its two representative applications: single-cell mRNA probing when integrated with a dielectrophoretic nanotweezer (DENT), and live-cell real-time imaging when combined with fluorescence lifetime imaging microscopy (FLIM). As the single-cell trapping efficiency is determined by the channel design instead of the flow rate, this trapping array can be coupled with different microfluidic sample processing units with different flow rates for various single-cell analyses.

Shift from stochastic to spatially-ordered expression of serine-glycine synthesis enzymes in 3D microtumors

Cell-to-cell differences in protein expression in normal tissues and tumors are a common phenomenon, but the underlying principles that govern this heterogeneity are largely unknown. Here, we show that in monolayer cancer cell-line cultures, the expression of the five metabolic enzymes of serine-glycine synthesis (SGS), including its rate-limiting enzyme, phosphoglycerate dehydrogenase (PHGDH), displays stochastic cell-to-cell variation. By contrast, in cancer cell line-derived three-dimensional (3D) microtumors PHGDH expression is restricted to the outermost part of the microtumors' outer proliferative cell layer, while the four other SGS enzymes display near uniform expression throughout the microtumor. A mathematical model suggests that metabolic stress in the microtumor core activates factors that restrict PHGDH expression. Thus, intracellular enzyme expression in growing cell ecosystems can shift to spatially ordered patterns in 3D structured environments due to emergent cell-cell communication, with potential implications for the design of effective anti-metabolic cancer therapies.

Phenotypic Plasticity, Bet-Hedging, and Androgen Independence in Prostate Cancer: Role of Non-Genetic Heterogeneity

It is well known that genetic mutations can drive drug resistance and lead to tumor relapse. Here, we focus on alternate mechanisms-those without mutations, such as phenotypic plasticity and stochastic cell-to-cell variability that can also evade drug attacks by giving rise to drug-tolerant persisters. The phenomenon of persistence has been well-studied in bacteria and has also recently garnered attention in cancer. We draw a parallel between bacterial persistence and resistance against androgen deprivation therapy in prostate cancer (PCa), the primary standard care for metastatic disease. We illustrate how phenotypic plasticity and consequent mutation-independent or non-genetic heterogeneity possibly driven by protein conformational dynamics can stochastically give rise to androgen independence in PCa, and suggest that dynamic phenotypic plasticity should be considered in devising therapeutic dosing strategies designed to treat and manage PCa.

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.

Beyond molecular tumor heterogeneity: protein synthesis takes control

One of the daunting challenges facing modern medicine lies in the understanding and treatment of tumor heterogeneity. Most tumors show intra-tumor heterogeneity at both genomic and proteomic levels, with marked impacts on the responses of therapeutic targets. Therapeutic target-related gene expression pathways are affected by hypoxia and cellular stress. However, the finding that targets such as eukaryotic initiation factor (eIF) 4E (and its phosphorylated form, p-eIF4E) are generally homogenously expressed throughout tumors, regardless of the presence of hypoxia or other cellular stress conditions, opens the exciting possibility that malignancies could be treated with therapies that combine targeting of eIF4E phosphorylation with immune checkpoint inhibitors or chemotherapy.

Digging deep into Golgi phenotypic diversity with unsupervised machine learning

Structural alterations of the Golgi apparatus may lead to phenotypes that human vision cannot easily discriminate. In this work, we present a high-content analysis framework including an unsupervised clustering step to automatically uncover Golgi phenotypic diversity. We use this deep phenotyping to quantitatively compare the effects of gene depletion.

The synthesis of glycans and the sorting of proteins are critical functions of the Golgi apparatus and depend on its highly complex and compartmentalized architecture. High-content image analysis coupled to RNA interference screening offers opportunities to explore this organelle organization and the gene network underlying it. To date, image-based Golgi screens have based on a single parameter or supervised analysis with predefined Golgi structural classes. Here, we report the use of multiparametric data extracted from a single marker and a computational unsupervised analysis framework to explore Golgi phenotypic diversity more extensively. In contrast with the three visually definable phenotypes, our framework reproducibly identified 10 Golgi phenotypes. They were used to quantify and stratify phenotypic similarities among genetic perturbations. The derived phenotypic network partially overlaps previously reported protein–protein interactions as well as suggesting novel functional interactions. Our workflow suggests the existence of multiple stable Golgi organizational states and provides a proof of concept for the classification of drugs and genes using fine-grained phenotypic information.

Phyllotactic regularity requires the Paf1 complex in Arabidopsis

In plants, aerial organs are initiated at stereotyped intervals, both spatially (every 137° in a pattern called phyllotaxis) and temporally (at prescribed time intervals called plastochrons). To investigate the molecular basis of such regularity, mutants with altered architecture have been isolated. However, most of them only exhibit plastochron defects and/or produce a new, albeit equally reproducible, phyllotactic pattern. This leaves open the question of a molecular control of phyllotaxis regularity. Here, we show that phyllotaxis regularity depends on the function of VIP proteins, components of the RNA polymerase II-associated factor 1 complex (Paf1c). Divergence angles between successive organs along the stem exhibited increased variance in vip3-1 and vip3-2 compared with the wild type, in two different growth conditions. Similar results were obtained with the weak vip3-6 allele and in vip6, a mutant for another Paf1c subunit. Mathematical analysis confirmed that these defects could not be explained solely by plastochron defects. Instead, increased variance in phyllotaxis in vip3 was observed at the meristem and related to defects in spatial patterns of auxin activity. Thus, the regularity of spatial, auxin-dependent, patterning at the meristem requires Paf1c.

Data-analysis strategies for image-based cell profiling

Image-based cell profiling is a high-throughput strategy for the quantification of phenotypic differences among a variety of cell populations. It paves the way to studying biological systems on a large scale by using chemical and genetic perturbations. The general workflow for this technology involves image acquisition with high-throughput microscopy systems and subsequent image processing and analysis. Here, we introduce the steps required to create high-quality image-based (i.e., morphological) profiles from a collection of microscopy images. We recommend techniques that have proven useful in each stage of the data analysis process, on the basis of the experience of 20 laboratories worldwide that are refining their image-based cell-profiling methodologies in pursuit of biological discovery. The recommended techniques cover alternatives that may suit various biological goals, experimental designs, and laboratories' preferences.

Combination Therapy Targeting BCL6 and Phospho-STAT3 Defeats Intratumor Heterogeneity in a Subset of Non-Small Cell Lung Cancers

Oncogene-specific changes in cellular signaling have been widely observed in lung cancer. Here, we investigated how these alterations could affect signaling heterogeneity and suggest novel therapeutic strategies. We compared signaling changes across six human bronchial epithelial cell (HBEC) strains that were systematically transformed with various combinations of TP53, KRAS, and MYC-oncogenic alterations commonly found in non-small cell lung cancer (NSCLC). We interrogated at single-cell resolution how these alterations could affect classic readouts (β-CATENIN, SMAD2/3, phospho-STAT3, P65, FOXO1, and phospho-ERK1/2) of key pathways commonly affected in NSCLC. All three oncogenic alterations were required concurrently to observe significant signaling changes, and significant heterogeneity arose in this condition. Unexpectedly, we found two mutually exclusive altered subpopulations: one with STAT3 upregulation and another with SMAD2/3 downregulation. Treatment with a STAT3 inhibitor eliminated the upregulated STAT3 subpopulation, but left a large surviving subpopulation with downregulated SMAD2/3. A bioinformatics search identified BCL6, a gene downstream of SMAD2/3, as a novel pharmacologically accessible target of our transformed HBECs. Combination treatment with STAT3 and BCL6 inhibitors across a panel of NSCLC cell lines and in xenografted tumors significantly reduced tumor cell growth. We conclude that BCL6 is a new therapeutic target in NSCLC and combination therapy that targets multiple vulnerabilities (STAT3 and BCL6) downstream of common oncogenes, and tumor suppressors may provide a potent way to defeat intratumor heterogeneity. Cancer Res; 77(11); 3070-81. ©2017 AACR.

An end-to-end software solution for the analysis of high-throughput single-cell migration data

The systematic study of single-cell migration requires the availability of software for assisting data inspection, quality control and analysis. This is especially important for high-throughput experiments, where multiple biological conditions are tested in parallel. Although the field of cell migration can count on different computational tools for cell segmentation and tracking, downstream data visualization, parameter extraction and statistical analysis are still left to the user and are currently not possible within a single tool. This article presents a completely new module for the open-source, cross-platform CellMissy software for cell migration data management. This module is the first tool to focus specifically on single-cell migration data downstream of image processing. It allows fast comparison across all tested conditions, providing automated data visualization, assisted data filtering and quality control, extraction of various commonly used cell migration parameters, and non-parametric statistical analysis. Importantly, the module enables parameters computation both at the trajectory- and at the step-level. Moreover, this single-cell analysis module is complemented by a new data import module that accommodates multiwell plate data obtained from high-throughput experiments, and is easily extensible through a plugin architecture. In conclusion, the end-to-end software solution presented here tackles a key bioinformatics challenge in the cell migration field, assisting researchers in their high-throughput data processing.

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.

Single cell dynamic phenotyping

Live cell imaging has improved our ability to measure phenotypic heterogeneity. However, bottlenecks in imaging and image processing often make it difficult to differentiate interesting biological behavior from technical artifact. Thus there is a need for new methods that improve data quality without sacrificing throughput. Here we present a 3-step workflow to improve dynamic phenotype measurements of heterogeneous cell populations. We provide guidelines for image acquisition, phenotype tracking, and data filtering to remove erroneous cell tracks using the novel Tracking Aberration Measure (TrAM). Our workflow is broadly applicable across imaging platforms and analysis software. By applying this workflow to cancer cell assays, we reduced aberrant cell track prevalence from 17% to 2%. The cost of this improvement was removing 15% of the well-tracked cells. This enabled detection of significant motility differences between cell lines. Similarly, we avoided detecting a false change in translocation kinetics by eliminating the true cause: varied proportions of unresponsive cells. Finally, by systematically seeking heterogeneous behaviors, we detected subpopulations that otherwise could have been missed, including early apoptotic events and pre-mitotic cells. We provide optimized protocols for specific applications and step-by-step guidelines for adapting them to a variety of biological systems.

High-Sensitivity Mass Spectrometry for Probing Gene Translation in Single Embryonic Cells in the Early Frog (Xenopus) Embryo

Direct measurement of protein expression with single-cell resolution promises to deepen the understanding of the basic molecular processes during normal and impaired development. High-resolution mass spectrometry provides detailed coverage of the proteomic composition of large numbers of cells. Here we discuss recent mass spectrometry developments based on single-cell capillary electrophoresis that extend discovery proteomics to sufficient sensitivity to enable the measurement of proteins in single cells. The single-cell mass spectrometry system is used to detect a large number of proteins in single embryonic cells in the 16-cell embryo of the South African clawed frog (Xenopus laevis) that give rise to distinct tissue types. Single-cell measurements of protein expression provide complementary information on gene transcription during early development of the vertebrate embryo, raising a potential to understand how differential gene expression coordinates normal cell heterogeneity during development.

The impact of mechanical compression on cortical microtubules in Arabidopsis: a quantitative pipeline

Exogenous mechanical perturbations on living tissues are commonly used to investigate whether cell effectors can respond to mechanical cues. However, in most of these experiments, the applied mechanical stress and/or the biological response are described only qualitatively. We developed a quantitative pipeline based on microindentation and image analysis to investigate the impact of a controlled and prolonged compression on microtubule behaviour in the Arabidopsis shoot apical meristem, using microtubule fluorescent marker lines. We found that a compressive stress, in the order of magnitude of turgor pressure, induced apparent microtubule bundling. Importantly, that response could be reversed several hours after the release of compression. Next, we tested the contribution of microtubule severing to compression-induced bundling: microtubule bundling seemed less pronounced in the katanin mutant, in which microtubule severing is dramatically reduced. Conversely, some microtubule bundles could still be observed 16 h after the release of compression in the spiral2 mutant, in which severing rate is instead increased. To quantify the impact of mechanical stress on anisotropy and orientation of microtubule arrays, we used the nematic tensor based FibrilTool ImageJ/Fiji plugin. To assess the degree of apparent bundling of the network, we developed several methods, some of which were borrowed from geostatistics. The final microtubule bundling response could notably be related to tissue growth velocity that was recorded by the indenter during compression. Because both input and output are quantified, this pipeline is an initial step towards correlating more precisely the cytoskeleton response to mechanical stress in living tissues.

Signal Transduction at the Single-Cell Level: Approaches to Study the Dynamic Nature of Signaling Networks

Signal transduction, or how cells interpret and react to external events, is a fundamental aspect of cellular function. Traditional study of signal transduction pathways involves mapping cellular signaling pathways at the population level. However, population-averaged readouts do not adequately illuminate the complex dynamics and heterogeneous responses found at the single-cell level. Recent technological advances that observe cellular response, computationally model signaling pathways, and experimentally manipulate cells now enable studying signal transduction at the single-cell level. These studies will enable deeper insights into the dynamic nature of signaling networks.

Exploiting Single-Cell Quantitative Data to Map Genetic Variants Having Probabilistic Effects

Despite the recent progress in sequencing technologies, genome-wide association studies (GWAS) remain limited by a statistical-power issue: many polymorphisms contribute little to common trait variation and therefore escape detection. The small contribution sometimes corresponds to incomplete penetrance, which may result from probabilistic effects on molecular regulations. In such cases, genetic mapping may benefit from the wealth of data produced by single-cell technologies. We present here the development of a novel genetic mapping method that allows to scan genomes for single-cell Probabilistic Trait Loci that modify the statistical properties of cellular-level quantitative traits. Phenotypic values are acquired on thousands of individual cells, and genetic association is obtained from a multivariate analysis of a matrix of Kantorovich distances. No prior assumption is required on the mode of action of the genetic loci involved and, by exploiting all single-cell values, the method can reveal non-deterministic effects. Using both simulations and yeast experimental datasets, we show that it can detect linkages that are missed by classical genetic mapping. A probabilistic effect of a single SNP on cell shape was detected and validated. The method also detected a novel locus associated with elevated gene expression noise of the yeast galactose regulon. Our results illustrate how single-cell technologies can be exploited to improve the genetic dissection of certain common traits. The method is available as an open source R package called ptlmapper.

Label-free Quantification of Proteins in Single Embryonic Cells with Neural Fate in the Cleavage-Stage Frog (Xenopus laevis) Embryo using Capillary Electrophoresis Electrospray Ionization High-Resolution Mass Spectrometry (CE-ESI-HRMS)

Quantification of protein expression in single cells promises to advance a systems-level understanding of normal development. Using a bottom-up proteomic workflow and multiplexing quantification by tandem mass tags, we recently demonstrated relative quantification between single embryonic cells (blastomeres) in the frog (Xenopus laevis) embryo. In this study, we minimize derivatization steps to enhance analytical sensitivity and use label-free quantification (LFQ) for single Xenopus cells. The technology builds on a custom-designed capillary electrophoresis microflow-electrospray ionization high-resolution mass spectrometry platform and LFQ by MaxLFQ (MaxQuant). By judiciously tailoring performance to peptide separation, ionization, and data-dependent acquisition, we demonstrate an ∼75-amol (∼11 nm) lower limit of detection and quantification for proteins in complex cell digests. The platform enabled the identification of 438 nonredundant protein groups by measuring 16 ng of protein digest, or <0.2% of the total protein contained in a blastomere in the 16-cell embryo. LFQ intensity was validated as a quantitative proxy for protein abundance. Correlation analysis was performed to compare protein quantities between the embryo and n = 3 different single D11 blastomeres, which are fated to develop into the nervous system. A total of 335 nonredundant protein groups were quantified in union between the single D11 cells spanning a 4 log-order concentration range. LFQ and correlation analysis detected expected proteomic differences between the whole embryo and blastomeres, and also found translational differences between individual D11 cells. LFQ on single cells raises exciting possibilities to study gene expression in other cells and models to help better understand cell processes on a systems biology level.

A Perspective on Implementing a Quantitative Systems Pharmacology Platform for Drug Discovery and the Advancement of Personalized Medicine

Drug candidates exhibiting well-defined pharmacokinetic and pharmacodynamic profiles that are otherwise safe often fail to demonstrate proof-of-concept in phase II and III trials. Innovation in drug discovery and development has been identified as a critical need for improving the efficiency of drug discovery, especially through collaborations between academia, government agencies, and industry. To address the innovation challenge, we describe a comprehensive, unbiased, integrated, and iterative quantitative systems pharmacology (QSP)-driven drug discovery and development strategy and platform that we have implemented at the University of Pittsburgh Drug Discovery Institute. Intrinsic to QSP is its integrated use of multiscale experimental and computational methods to identify mechanisms of disease progression and to test predicted therapeutic strategies likely to achieve clinical validation for appropriate subpopulations of patients. The QSP platform can address biological heterogeneity and anticipate the evolution of resistance mechanisms, which are major challenges for drug development. The implementation of this platform is dedicated to gaining an understanding of mechanism(s) of disease progression to enable the identification of novel therapeutic strategies as well as repurposing drugs. The QSP platform will help promote the paradigm shift from reactive population-based medicine to proactive personalized medicine by focusing on the patient as the starting and the end point.

Leveraging Large-scale Behavioral Profiling in Zebrafish to Explore Neuroactive Polypharmacology

Many psychiatric drugs modulate the nervous system through multitarget mechanisms. However, systematic identification of multitarget compounds has been difficult using traditional in vitro screening assays. New approaches to phenotypic profiling in zebrafish can help researchers identify novel compounds with complex polypharmacology. For example, large-scale behavior-based chemical screens can rapidly identify large numbers of structurally diverse and phenotype-related compounds. Once these compounds have been identified, a systems-level analysis of their structures may help to identify statistically enriched target pathways. Together, systematic behavioral profiling and multitarget predictions may help researchers identify new behavior-modifying pathways and CNS therapeutics.

Visualizing quantitative microscopy data: History and challenges

Data visualization is a fundamental aspect of science. In the context of microscopy-based studies, visualization typically involves presentation of the images themselves. However, data visualization is challenging when microscopy experiments entail imaging of millions of cells, and complex cellular phenotypes are quantified in a high-content manner. Most well-established visualization tools are inappropriate for displaying high-content data, which has driven the development of new visualization methodology. In this review, we discuss how data has been visualized in both classical and high-content microscopy studies; as well as the advantages, and disadvantages, of different visualization methods.

Variation is function: Are single cell differences functionally important?: Testing the hypothesis that single cell variation is required for aggregate function

There is a growing appreciation of the extent of transcriptome variation across individual cells of the same cell type. While expression variation may be a byproduct of, for example, dynamic or homeostatic processes, here we consider whether single-cell molecular variation per se might be crucial for population-level function. Under this hypothesis, molecular variation indicates a diversity of hidden functional capacities within an ensemble of identical cells, and this functional diversity facilitates collective behavior that would be inaccessible to a homogenous population. In reviewing this topic, we explore possible functions that might be carried by a heterogeneous ensemble of cells; however, this question has proven difficult to test, both because methods to manipulate molecular variation are limited and because it is complicated to define, and measure, population-level function. We consider several possible methods to further pursue the hypothesis that variation is function through the use of comparative analysis and novel experimental techniques.

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.

Intra-tumor heterogeneity of cancer cells and its implications for cancer treatment

Recent studies have revealed extensive genetic and non-genetic variation across different geographical regions of a tumor or throughout different stages of tumor progression, which is referred to as intra-tumor heterogeneity. Several causes contribute to this phenomenon, including genomic instability, epigenetic alteration, plastic gene expression, signal transduction, and microenvironmental differences. These variables may affect key signaling pathways that regulate cancer cell growth, drive phenotypic diversity, and pose challenges to cancer treatment. Understanding the mechanisms underlying this heterogeneity will support the development of effective therapeutic strategies.

Highly multiplexed imaging of single cells using a high-throughput cyclic immunofluorescence method

Single-cell analysis reveals aspects of cellular physiology not evident from population-based studies, particularly in the case of highly multiplexed methods such as mass cytometry (CyTOF) able to correlate the levels of multiple signalling, differentiation and cell fate markers. Immunofluorescence (IF) microscopy adds information on cell morphology and the microenvironment that are not obtained using flow-based techniques, but the multiplicity of conventional IF is limited. This has motivated development of imaging methods that require specialized instrumentation, exotic reagents or proprietary protocols that are difficult to reproduce in most laboratories. Here we report a public-domain method for achieving high multiplicity single-cell IF using cyclic immunofluorescence (CycIF), a simple and versatile procedure in which four-colour staining alternates with chemical inactivation of fluorophores to progressively build a multichannel image. Because CycIF uses standard reagents and instrumentation and is no more expensive than conventional IF, it is suitable for high-throughput assays and screening applications.

Multiplexed single cell measurements provide insight into connections between cell state and phenotype. Here Lin et al. present CycIF, a high throughput, public domain immunofluorescence method for multiplexed single-cell analysis of adherent cells following live-cell imaging.

Stochastic sensitivity analysis and kernel inference via distributional data

Cellular processes are noisy due to the stochastic nature of biochemical reactions. As such, it is impossible to predict the exact quantity of a molecule or other attributes at the single-cell level. However, the distribution of a molecule over a population is often deterministic and is governed by the underlying regulatory networks relevant to the cellular functionality of interest. Recent studies have started to exploit this property to infer network states. To facilitate the analysis of distributional data in a general experimental setting, we introduce a computational framework to efficiently characterize the sensitivity of distributional output to changes in external stimuli. Further, we establish a probability-divergence-based kernel regression model to accurately infer signal level based on distribution measurements. Our methodology is applicable to any biological system subject to stochastic dynamics and can be used to elucidate how population-based information processing may contribute to organism-level functionality. It also lays the foundation for engineering synthetic biological systems that exploit population decoding to more robustly perform various biocomputation tasks, such as disease diagnostics and environmental-pollutant sensing.

Modeling Signaling Networks to Advance New Cancer Therapies

Cell signaling pathways control cells' responses to their environment through an intricate network of proteins and small molecules partitioned by intracellular structures, such as the cytoskeleton and nucleus. Our understanding of these pathways has been revised recently with the advent of more advanced experimental techniques; no longer are signaling pathways viewed as linear cascades of information flowing from membrane-bound receptors to the nucleus. Instead, such pathways must be understood in the context of networks, and studying such networks requires an integration of computational and experimental approaches. This understanding is becoming more important in designing novel therapies for diseases such as cancer. Using the MAPK (mitogen-activated protein kinase) and PI3K (class I phosphoinositide-3' kinase) pathways as case studies of cellular signaling, we give an overview of these pathways and their functions. We then describe, using a number of case studies, how computational modeling has aided in understanding these pathways' deregulation in cancer, and how such understanding can be used to optimally tailor current therapies or help design new therapies against cancer.

Surviving apoptosis: life-death signaling in single cells

Tissue development and homeostasis are regulated by opposing pro-survival and pro-death signals. An interesting feature of the Tumor Necrosis Factor (TNF) family of ligands is that they simultaneously activate opposing signals within a single cell via the same ligand-receptor complex. The magnitude of pro-death events such as caspase activation and pro-survival events such as Nuclear Factor (NF)-κB activation vary not only from one cell type to the next but also among individual cells of the same type due to intrinsic and extrinsic noise. The molecules involved in these pro-survival and/or pro-death pathways, and the different phenotypes that result from their activities, have been recently reviewed. Here we focus on the impact of cell-to-cell variability in the strength of these opposing signals on shaping cell fate decisions.

Analysis of single-cell cytokine secretion reveals a role for paracrine signaling in coordinating macrophage responses to TLR4 stimulation

Macrophages not only produce multiple cytokines but also respond to multiple cytokines, which likely shapes the ultimate response of the population. To determine the role of paracrine signaling in shaping the profile of inflammatory cytokines secreted by macrophages in response to stimulation of Toll-like receptor 4 (TLR4) with lipopolysaccharide (LPS), we combined multiplexed, microwell-based measurements of cytokine secretion by single cells with analysis of cytokine secretion by cell populations. Loss of paracrine signaling as a result of cell isolation reduced the secretion by macrophage-like U937 cells and human monocyte-derived macrophages (MDMs) of a subset of LPS-stimulated cytokines, including interleukin-6 (IL-6) and IL-10. Graphical Gaussian modeling (GGM) of the single-cell data defined a regulatory network of paracrine signals, which was validated experimentally in the population through antibody-mediated neutralization of individual cytokines. Tumor necrosis factor-α (TNF-α) was the most influential cytokine in the GGM network. Paracrine signaling by TNF-α secreted from a small subpopulation of "high-secreting" cells was necessary, but not sufficient, for the secretion of large amounts of IL-6 and IL-10 by the cell population. Decreased relative IL-10 secretion by isolated MDMs was linked to increased TNF-α secretion, suggesting that inhibition of the inflammatory response also depends on paracrine signaling. Our results reveal a previously uncharacterized role for cell-to-cell communication within a population in coordinating a rapid innate immune response despite underlying cell-to-cell heterogeneity.

Akt/PKB: one kinase, many modifications

Akt/PKB, a serine/threonine kinase member of the AGC family of proteins, is involved in the regulation of a plethora of cellular processes triggered by a wide diversity of extracellular signals and is thus considered a key signalling molecule in higher eukaryotes. Deregulation of Akt signalling is associated with a variety of human diseases, revealing Akt-dependent pathways as an attractive target for therapeutic intervention. Since its discovery in the early 1990s, a large body of work has focused on Akt phosphorylation of two residues, Thr308 and Ser473, and modification of these two sites has been established as being equivalent to Akt activation. More recently, Akt has been identified as a substrate for many different post-translational modifications, including not only phosphorylation of other residues, but also acetylation, glycosylation, oxidation, ubiquitination and SUMOylation. These modifications could provide additional regulatory steps for fine-tuning Akt function, Akt trafficking within the cell and/or for determining the substrate specificity of this signalling molecule. In the present review, we provide an overview of these different post-translational modifications identified for Akt, focusing on their consequences for this kinase activity.