Computational methods for characterizing and learning from heterogeneous cell signaling data

Stimulus-response signaling dynamics characterize macrophage polarization states

The functional state of cells is dependent on their microenvironmental context. Prior studies described how polarizing cytokines alter macrophage transcriptomes and epigenomes. Here, we characterized the functional responses of 6 differentially polarized macrophage populations by measuring the dynamics of transcription factor nuclear factor κB (NF-κB) in response to 8 stimuli. The resulting dataset of single-cell NF-κB trajectories was analyzed by three approaches: (1) machine learning on time-series data revealed losses of stimulus distinguishability with polarization, reflecting canalized effector functions. (2) Informative trajectory features driving stimulus distinguishability ("signaling codons") were identified and used for mapping a cell state landscape that could then locate macrophages conditioned by an unrelated condition. (3) Kinetic parameters, inferred using a mechanistic NF-κB network model, provided an alternative mapping of cell states and correctly predicted biochemical findings. Together, this work demonstrates that a single analyte's dynamic trajectories may distinguish the functional states of single cells and molecular network states underlying them. A record of this paper's transparent peer review process is included in the supplemental information.

Characterizing heterogeneous single-cell dose responses computationally and experimentally using threshold inhibition surfaces and dose-titration assays

Single cancer cells within a tumor exhibit variable levels of resistance to drugs, ultimately leading to treatment failures. While tumor heterogeneity is recognized as a major obstacle to cancer therapy, standard dose-response measurements for the potency of targeted kinase inhibitors aggregate populations of cells, obscuring intercellular variations in responses. In this work, we develop an analytical and experimental framework to quantify and model dose responses of individual cancer cells to drugs. We first explore the connection between population and single-cell dose responses using a computational model, revealing that multiple heterogeneous populations can yield nearly identical population dose responses. We demonstrate that a single-cell analysis method, which we term a threshold inhibition surface, can differentiate among these populations. To demonstrate the applicability of this method, we develop a dose-titration assay to measure dose responses in single cells. We apply this assay to breast cancer cells responding to phosphatidylinositol-3-kinase inhibition (PI3Ki), using clinically relevant PI3Kis on breast cancer cell lines expressing fluorescent biosensors for kinase activity. We demonstrate that MCF-7 breast cancer cells exhibit heterogeneous dose responses with some cells requiring over ten-fold higher concentrations than the population average to achieve inhibition. Our work reimagines dose-response relationships for cancer drugs in an emerging paradigm of single-cell tumor heterogeneity.

Quantifying the phenotypic information in mRNA abundance

Quantifying the dependency between mRNA abundance and downstream cellular phenotypes is a fundamental open problem in biology. Advances in multimodal single‐cell measurement technologies provide an opportunity to apply new computational frameworks to dissect the contribution of individual genes and gene combinations to a given phenotype. Using an information theory approach, we analyzed multimodal data of the expression of 83 genes in the Ca²⁺ signaling network and the dynamic Ca²⁺ response in the same cell. We found that the overall expression levels of these 83 genes explain approximately 60% of Ca²⁺ signal entropy. The average contribution of each single gene was 17%, revealing a large degree of redundancy between genes. Using different heuristics, we estimated the dependency between the size of a gene set and its information content, revealing that on average, a set of 53 genes contains 54% of the information about Ca²⁺ signaling. Our results provide the first direct quantification of information content about complex cellular phenotype that exists in mRNA abundance measurements.

D-Xylose Sensing in Saccharomyces cerevisiae: Insights from D-Glucose Signaling and Native D-Xylose Utilizers

Extension of the substrate range is among one of the metabolic engineering goals for microorganisms used in biotechnological processes because it enables the use of a wide range of raw materials as substrates. One of the most prominent examples is the engineering of baker's yeast Saccharomyces cerevisiae for the utilization of d-xylose, a five-carbon sugar found in high abundance in lignocellulosic biomass and a key substrate to achieve good process economy in chemical production from renewable and non-edible plant feedstocks. Despite many excellent engineering strategies that have allowed recombinant S. cerevisiae to ferment d-xylose to ethanol at high yields, the consumption rate of d-xylose is still significantly lower than that of its preferred sugar d-glucose. In mixed d-glucose/d-xylose cultivations, d-xylose is only utilized after d-glucose depletion, which leads to prolonged process times and added costs. Due to this limitation, the response on d-xylose in the native sugar signaling pathways has emerged as a promising next-level engineering target. Here we review the current status of the knowledge of the response of S. cerevisiae signaling pathways to d-xylose. To do this, we first summarize the response of the native sensing and signaling pathways in S. cerevisiae to d-glucose (the preferred sugar of the yeast). Using the d-glucose case as a point of reference, we then proceed to discuss the known signaling response to d-xylose in S. cerevisiae and current attempts of improving the response by signaling engineering using native targets and synthetic (non-native) regulatory circuits.