Identifying single-cell molecular programs by stochastic profiling

Patient-derived response estimates from zero-passage organoids of luminal breast cancer

Background

Primary luminal breast cancer cells lose their identity rapidly in standard tissue culture, which is problematic for testing hormone interventions and molecular pathways specific to the luminal subtype. Breast cancer organoids are thought to retain tumor characteristics better, but long-term viability of luminal-subtype cases is a persistent challenge. Our goal was to adapt short-term organoids of luminal breast cancer for parallel testing of genetic and pharmacologic perturbations.

Methods

We freshly isolated patient-derived cells from luminal tumor scrapes, miniaturized the organoid format into 5 μl replicates for increased throughput, and set an endpoint of 14 days to minimize drift. Therapeutic hormone targeting was mimicked in these "zero-passage" organoids by withdrawing β-estradiol and adding 4-hydroxytamoxifen. We also examined sulforaphane as an electrophilic stress and commercial neutraceutical with reported anti-cancer properties. Downstream mechanisms were tested genetically by lentiviral transduction of two complementary sgRNAs and Cas9 stabilization for the first week of organoid culture. Transcriptional changes were measured by RT-qPCR or RNA sequencing, and organoid phenotypes were quantified by serial brightfield imaging, digital image segmentation, and regression modeling of cellular doubling times.

Results

We achieved >50% success in initiating luminal breast cancer organoids from tumor scrapes and maintaining them to the 14-day zero-passage endpoint. Success was mostly independent of clinical parameters, supporting general applicability of the approach. Abundance of ESR1 and PGR in zero-passage organoids consistently remained within the range of patient variability at the endpoint. However, responsiveness to hormone withdrawal and blockade was highly variable among luminal breast cancer cases tested. Combining sulforaphane with knockout of NQO1 (a phase II antioxidant response gene and downstream effector of sulforaphane) also yielded a breadth of organoid growth phenotypes, including growth inhibition with sulforaphane, growth promotion with NQO1 knockout, and growth antagonism when combined.

Conclusions

Zero-passage organoids are a rapid and scalable way to interrogate properties of luminal breast cancer cells from patient-derived material. This includes testing drug mechanisms of action in different clinical cohorts. A future goal is to relate inter-patient variability of zero-passage organoids to long-term outcomes.

Deciphering the History of ERK Activity from Fixed-Cell Immunofluorescence Measurements

The Ras/ERK pathway drives cell proliferation and other oncogenic behaviors, and quantifying its activity in situ is of high interest in cancer diagnosis and therapy. Pathway activation is often assayed by measuring phosphorylated ERK. However, this form of measurement overlooks dynamic aspects of signaling that can only be observed over time. In this study, we combine a live, single-cell ERK biosensor approach with multiplexed immunofluorescence staining of downstream target proteins to ask how well immunostaining captures the dynamic history of ERK activity. Combining linear regression, machine learning, and differential equation models, we develop an interpretive framework for immunostains, in which Fra-1 and pRb levels imply long term activation of ERK signaling, while Egr-1 and c-Myc indicate recent activation. We show that this framework can distinguish different classes of ERK dynamics within a heterogeneous population, providing a tool for annotating ERK dynamics within fixed tissues.

Continuous sensing of nutrients and growth factors by the mTORC1-TFEB axis

mTORC1 senses nutrients and growth factors and phosphorylates downstream targets, including the transcription factor TFEB, to coordinate metabolic supply and demand. These functions position mTORC1 as a central controller of cellular homeostasis, but the behavior of this system in individual cells has not been well characterized. Here, we provide measurements necessary to refine quantitative models for mTORC1 as a metabolic controller. We developed a series of fluorescent protein-TFEB fusions and a multiplexed immunofluorescence approach to investigate how combinations of stimuli jointly regulate mTORC1 signaling at the single-cell level. Live imaging of individual MCF10A cells confirmed that mTORC1-TFEB signaling responds continuously to individual, sequential, or simultaneous treatment with amino acids and the growth factor insulin. Under physiologically relevant concentrations of amino acids, we observe correlated fluctuations in TFEB, AMPK, and AKT signaling that indicate continuous activity adjustments to nutrient availability. Using partial least squares regression modeling, we show that these continuous gradations are connected to protein synthesis rate via a distributed network of mTORC1 effectors, providing quantitative support for the qualitative model of mTORC1 as a homeostatic controller and clarifying its functional behavior within individual cells.

pYtags enable spatiotemporal measurements of receptor tyrosine kinase signaling in living cells

Receptor tyrosine kinases (RTKs) are major signaling hubs in metazoans, playing crucial roles in cell proliferation, migration, and differentiation. However, few tools are available to measure the activity of a specific RTK in individual living cells. Here, we present pYtags, a modular approach for monitoring the activity of a user-defined RTK by live-cell microscopy. pYtags consist of an RTK modified with a tyrosine activation motif that, when phosphorylated, recruits a fluorescently labeled tandem SH2 domain with high specificity. We show that pYtags enable the monitoring of a specific RTK on seconds-to-minutes time scales and across subcellular and multicellular length scales. Using a pYtag biosensor for epidermal growth factor receptor (EGFR), we quantitatively characterize how signaling dynamics vary with the identity and dose of activating ligand. We show that orthogonal pYtags can be used to monitor the dynamics of EGFR and ErbB2 activity in the same cell, revealing distinct phases of activation for each RTK. The specificity and modularity of pYtags open the door to robust biosensors of multiple tyrosine kinases and may enable engineering of synthetic receptors with orthogonal response programs.

Nucleocytoplasmic transport of active HER2 causes fractional escape from the DCIS-like state

Activation of HER2/ErbB2 coincides with escape from ductal carcinoma in situ (DCIS) premalignancy and disrupts 3D organization of cultured breast-epithelial spheroids. The 3D phenotype is infrequent, however, and mechanisms for its incomplete penetrance have been elusive. Using inducible HER2/ErbB2-EGFR/ErbB1 heterodimers, we match phenotype penetrance to the frequency of co-occurring transcriptomic changes and uncover a reconfiguration in the karyopherin network regulating ErbB nucleocytoplasmic transport. Induction of the exportin CSE1L inhibits nuclear accumulation of ErbBs, whereas nuclear ErbBs silence the importin KPNA1 by inducing miR-205. When these negative feedbacks are incorporated into a validated systems model of nucleocytoplasmic transport, steady-state localization of ErbB cargo becomes ultrasensitive to initial CSE1L abundance. Erbb2-driven carcinomas with Cse1l deficiency outgrow less irregularly from mammary ducts, and NLS-attenuating mutants or variants of HER2 favor escape in 3D culture. We conclude here that adaptive nucleocytoplasmic relocalization of HER2 creates a systems-level molecular switch at the premalignant-to-malignant transition.

ANKLE1 cleaves mitochondrial DNA and contributes to cancer risk by promoting apoptosis resistance and metabolic dysregulation

Alleles within the chr19p13.1 locus are associated with increased risk of both ovarian and breast cancer and increased expression of the ANKLE1 gene. ANKLE1 is molecularly characterized as an endonuclease that efficiently cuts branched DNA and shuttles between the nucleus and cytoplasm. However, the role of ANKLE1 in mammalian development and homeostasis remains unknown. In normal development ANKLE1 expression is limited to the erythroblast lineage and we found that ANKLE1's role is to cleave the mitochondrial genome during erythropoiesis. We show that ectopic expression of ANKLE1 in breast epithelial-derived cells leads to genome instability and mitochondrial DNA (mtDNA) cleavage. mtDNA degradation then leads to mitophagy and causes a shift from oxidative phosphorylation to glycolysis (Warburg effect). Moreover, mtDNA degradation activates STAT1 and expression of epithelial-mesenchymal transition (EMT) genes. Reduction in mitochondrial content contributes to apoptosis resistance, which may allow precancerous cells to avoid apoptotic checkpoints and proliferate. These findings provide evidence that ANKLE1 is the causal cancer susceptibility gene in the chr19p13.1 locus and describe mechanisms by which higher ANKLE1 expression promotes cancer risk.

TGFBR3 supports anoikis through suppressing ATF4 signaling

Epithelial morphogenesis and oncogenic transformation can cause loss of cell adhesion, and detached cells are eliminated by anoikis. Here, we reveal that transforming growth factor beta receptor 3 (TGFBR3) acts as an anoikis mediator through the coordination of activating transcription factor 4 (ATF4). In breast cancer, TGFBR3 is progressively lost, but elevated TGFBR3 is associated with a histologic subtype characterized by cellular adhesion defects. Dissecting the impact of extracellular matrix (ECM) deprivation, we demonstrate that ECM loss promotes TGFBR3 expression, which in turn differentiates cell aggregates to a prosurvival phenotype and drives the intrinsic apoptotic pathway. We demonstrate that inhibition of TGFBR3 impairs epithelial anoikis by activating ATF4 signaling. These preclinical findings provide a rationale for therapeutic inhibition of ATF4 in the subgroup of breast cancer patients with low TGFBR3 expression.

Substratum stiffness regulates Erk signaling dynamics through receptor-level control

The EGFR/Erk pathway is triggered by extracellular ligand stimulation, leading to stimulus-dependent dynamics of pathway activity. Although mechanical properties of the microenvironment also affect Erk activity, their effects on Erk signaling dynamics are poorly understood. Here, we characterize how the stiffness of the underlying substratum affects Erk signaling dynamics in mammary epithelial cells. We find that soft microenvironments attenuate Erk signaling, both at steady state and in response to epidermal growth factor (EGF) stimulation. Optogenetic manipulation at multiple signaling nodes reveals that intracellular signal transmission is largely unaffected by substratum stiffness. Instead, we find that soft microenvironments decrease EGF receptor (EGFR) expression and alter the amount and spatial distribution of EGF binding at cell membranes. Our data demonstrate that the mechanical microenvironment tunes Erk signaling dynamics via receptor-ligand interactions, underscoring how multiple microenvironmental signals are jointly processed through a highly conserved pathway that regulates tissue development, homeostasis, and disease progression.

Metabolic Stress Adaptations Underlie Mammary Gland Morphogenesis and Breast Cancer Progression

Breast cancers display dynamic reprogrammed metabolic activities as cancers develop from premalignant lesions to primary tumors, and then metastasize. Numerous advances focus on how tumors develop pro-proliferative metabolic signaling that differs them from adjacent, non-transformed epithelial tissues. This leads to targetable oncogene-driven liabilities among breast cancer subtypes. Other advances demonstrate how microenvironments trigger stress-response at single-cell resolution. Microenvironmental heterogeneities give rise to cell regulatory states in cancer cell spheroids in three-dimensional cultures and at stratified terminal end buds during mammary gland morphogenesis, where stress and survival signaling juxtapose. The cell-state specificity in stress signaling networks recapture metabolic evolution during cancer progression. Understanding lineage-specific metabolic phenotypes in experimental models is useful for gaining a deeper understanding of subtype-selective breast cancer metabolism.

Laser Capture Proteomics: spatial tissue molecular profiling from the bench to personalized medicine

INTRODUCTION

Laser Capture Microdissection (LCM) uses a laser to isolate, or capture, specific cells of interest in a complex heterogeneous tissue section, under direct microscopic visualization. Recently, there has been a surge of publications using LCM for tissue spatial molecular profiling relevant to a wide range of research topics.

AREAS COVERED

We summarize the many advances in tissue Laser Capture Proteomics (LCP) using mass spectrometry for discovery, and protein arrays for signal pathway network mapping. This review emphasizes: a) transition of LCM phosphoproteomics from the lab to the clinic for individualized cancer therapy, and b) the emerging frontier of LCM single cell molecular analysis combining proteomics with genomic, and transcriptomic analysis. The search strategy was based on the combination of MeSH terms with expert refinement.

EXPERT OPINION

LCM is complemented by a rich set of instruments, methodology protocols, and analytical A.I. (artificial intelligence) software for basic and translational research. Resolution is advancing to the tissue single cell level. A vision for the future evolution of LCM is presented. Emerging LCM technology is combining digital and AI guided remote imaging with automation, and telepathology, to a achieve multi-omic profiling that was not previously possible.

Fragmentation of Small-Cell Lung Cancer Regulatory States in Heterotypic Microenvironments

Small-cell lung cancers derive from pulmonary neuroendocrine cells, which have stem-like properties to reprogram into other cell types upon lung injury. It is difficult to uncouple transcriptional plasticity of these transformed cells from genetic changes that evolve in primary tumors or secondary metastases. Profiling of single cells is also problematic if the required sample dissociation activates injury-like signaling and reprogramming. Here we defined cell-state heterogeneities in situ through laser capture microdissection-based 10-cell transcriptomics coupled with stochastic-profiling fluctuation analysis. In labeled cells from a small-cell lung cancer mouse model initiated by neuroendocrine deletion of Rb1-Trp53, variations in transcript abundance revealed cell-to-cell differences in regulatory state in vitro and in vivo. Fluctuating transcripts in spheroid culture were partly shared among Rb1-Trp53-null models, and heterogeneities increased considerably when cells were delivered intravenously to colonize the liver. Colonization of immunocompromised animals drove a fractional appearance of alveolar type II-like markers and poised cells for paracrine stimulation from immune cells and hepatocytes. Immunocompetency further exaggerated the fragmentation of tumor states in the liver, yielding mixed stromal signatures evident in bulk sequencing from autochthonous tumors and metastases. Dozens of transcript heterogeneities recurred irrespective of biological context; their mapped orthologs brought together observations of murine and human small-cell lung cancer. Candidate heterogeneities recurrent in the liver also stratified primary human tumors into discrete groups not readily explained by molecular subtype but with prognostic relevance. These data suggest that heterotypic interactions in the liver and lung are an accelerant for intratumor heterogeneity in small-cell lung cancer. SIGNIFICANCE: These findings demonstrate that the single-cell regulatory heterogeneity of small-cell lung cancer becomes increasingly elaborate in the liver, a common metastatic site for the disease.See related articles by Singh and colleagues, p. 1840 and Sutcliffe and colleagues, p. 1868.

Pan-Cancer Drivers Are Recurrent Transcriptional Regulatory Heterogeneities in Early-Stage Luminal Breast Cancer

The heterogeneous composition of solid tumors is known to impact disease progression and response to therapy. Malignant cells coexist in different regulatory states that can be accessed transcriptomically by single-cell RNA sequencing, but these methods have many caveats related to sensitivity, noise, and sample handling. We revised a statistical fluctuation analysis called stochastic profiling to combine with 10-cell RNA sequencing, which was designed for laser-capture microdissection (LCM) and extended here for immuno-LCM. When applied to a cohort of late-onset, early-stage luminal breast cancers, the integrated approach identified thousands of candidate regulatory heterogeneities. Intersecting the candidates from different tumors yielded a relatively stable set of 710 recurrent heterogeneously expressed genes (RHEG), which were significantly variable in >50% of patients. RHEGs were not strongly confounded by dissociation artifacts, cell-cycle oscillations, or driving mutations for breast cancer. Rather, RHEGs were enriched for epithelial-to-mesenchymal transition genes and, unexpectedly, the latest pan-cancer assembly of driver genes across cancer types other than breast. These findings indicate that heterogeneous transcriptional regulation conceivably provides a faster, reversible mechanism for malignant cells to evaluate the effects of potential oncogenes or tumor suppressors on cancer hallmarks. SIGNIFICANCE: Profiling intratumor heterogeneity of luminal breast carcinoma cells identifies a recurrent set of genes, suggesting sporadic activation of pathways known to drive other types of cancer.See related articles by Schaff and colleagues, p. 1853 and Sutcliffe and colleagues, p. 1868.

Premalignant Oligodendrocyte Precursor Cells Stall in a Heterogeneous State of Replication Stress Prior to Gliomagenesis

Cancer evolves from premalignant clones that adopt unusual cell states to achieve transformation. We previously pinpointed the oligodendrocyte precursor cell (OPC) as a cell of origin for glioma, but the early changes of mutant OPCs during premalignancy remained unknown. Using mice engineered for inducible Nf1-Trp53 loss in OPCs, we acutely isolated labeled mutant OPCs by laser-capture microdissection, determined global gene-expression changes by bulk RNA sequencing, and compared with cell-state fluctuations at the single-cell level by stochastic profiling, which uses RNA-sequencing measurements from random pools of 10 mutant cells. At 12 days after Nf1-Trp53 deletion, bulk differences were mostly limited to mitotic hallmarks and genes for ribosome biosynthesis, and stochastic profiling revealed a spectrum of stem-progenitor (Axl, Aldh1a1), proneural, and mesenchymal states as potential starting points for gliomagenesis. At 90 days, bulk sequencing detected few differentially expressed transcripts, whereas stochastic profiling revealed cell states for neurons and mural cells that do not give rise to glial tumors, suggesting cellular dead-ends for gliomagenesis. Importantly, mutant OPCs that strongly expressed key effectors of nonsense-mediated decay (Upf3b) and homology-dependent DNA repair (Rad51c, Slx1b, Ercc4) were identified along with DNA-damage markers, suggesting transcription-associated replication stress. Analysis of 10-cell transcriptomes at 90 days identified a locus of elevated gene expression containing an additional repair endonuclease (Mus81) and Rin1, a Ras-Raf antagonist and possible counterbalance to Nf1 loss, which was microdeleted or downregulated in gliomas at 150 days. These hidden cell-state variations uncover replication stress as a potential bottleneck that must be resolved for glioma initiation. SIGNIFICANCE: Profiling premalignant cell states in a mouse model of glioma uncovers regulatory heterogeneity in glioma cells-of-origin and defines a state of replication stress that precedes tumor initiation.See related articles by Singh and colleagues, p. 1840 and Schaff and colleagues, p. 1853.

stochprofML: stochastic profiling using maximum likelihood estimation in R

BACKGROUND

Tissues are often heterogeneous in their single-cell molecular expression, and this can govern the regulation of cell fate. For the understanding of development and disease, it is important to quantify heterogeneity in a given tissue.

RESULTS

We present the R package stochprofML which uses the maximum likelihood principle to parameterize heterogeneity from the cumulative expression of small random pools of cells. We evaluate the algorithm's performance in simulation studies and present further application opportunities.

CONCLUSION

Stochastic profiling outweighs the necessary demixing of mixed samples with a saving in experimental cost and effort and less measurement error. It offers possibilities for parameterizing heterogeneity, estimating underlying pool compositions and detecting differences between cell populations between samples.

Transient phases of OXPHOS inhibitor resistance reveal underlying metabolic heterogeneity in single cells

Cell-to-cell heterogeneity in metabolism plays an unknown role in physiology and pharmacology. To functionally characterize cellular variability in metabolism, we treated cells with inhibitors of oxidative phosphorylation (OXPHOS) and monitored their responses with live-cell reporters for ATP, ADP/ATP, or activity of the energy-sensing kinase AMPK. Across multiple OXPHOS inhibitors and cell types, we identified a subpopulation of cells resistant to activation of AMPK and reduction of ADP/ATP ratio. This resistant state persists transiently for at least several hours and can be inherited during cell divisions. OXPHOS inhibition suppresses the mTORC1 and ERK growth signaling pathways in sensitive cells, but not in resistant cells. Resistance is linked to a multi-factorial combination of increased glucose uptake, reduced protein biosynthesis, and G0/G1 cell-cycle status. Our results reveal dynamic fluctuations in cellular energetic balance and provide a basis for measuring and predicting the distribution of cellular responses to OXPHOS inhibition.

Systems-Level Properties of EGFR-RAS-ERK Signaling Amplify Local Signals to Generate Dynamic Gene Expression Heterogeneity

Intratumoral heterogeneity is associated with aggressive tumor behavior, therapy resistance, and poor patient outcomes. Such heterogeneity is thought to be dynamic, shifting over periods of minutes to hours in response to signaling inputs from the tumor microenvironment. However, models of this process have been inferred from indirect or post-hoc measurements of cell state, leaving the temporal details of signaling-driven heterogeneity undefined. Here, we developed a live-cell model system in which microenvironment-driven signaling dynamics can be directly observed and linked to variation in gene expression. Our analysis reveals that paracrine signaling between two cell types is sufficient to drive continual diversification of gene expression programs. This diversification emerges from systems-level properties of the EGFR-RAS-ERK signaling cascade, including intracellular amplification of amphiregulin-mediated paracrine signals and differential kinetic filtering by target genes including Fra-1, c-Myc, and Egr1. Our data enable more precise modeling of paracrine-driven transcriptional variation as a generator of gene expression heterogeneity. A record of this paper's transparent peer review process is included in the Supplemental Information.

Sporadic activation of an oxidative stress-dependent NRF2-p53 signaling network in breast epithelial spheroids and premalignancies

Breast and mammary epithelial cells experience different local environments during tissue development and tumorigenesis. Microenvironmental heterogeneity gives rise to distinct cell regulatory states whose identity and importance are just beginning to be appreciated. Cellular states diversify when clonal three-dimensional (3D) spheroids are cultured in basement membrane, and one such state is associated with stress tolerance and poor response to anticancer therapeutics. Here, we found that this state was jointly coordinated by the NRF2 and p53 pathways, which were costabilized by spontaneous oxidative stress within 3D cultures. Inhibition of NRF2 or p53 individually disrupted some of the transcripts defining the regulatory state but did not yield a notable phenotype in nontransformed breast epithelial cells. In contrast, combined perturbation prevented 3D growth in an oxidative stress-dependent manner. By integrating systems models of NRF2 and p53 signaling in a single oxidative stress network, we recapitulated these observations and made predictions about oxidative stress profiles during 3D growth. NRF2 and p53 signaling were similarly coordinated in normal breast epithelial tissue and hormone-negative ductal carcinoma in situ lesions but were uncoupled in triple-negative breast cancer (TNBC), a subtype in which p53 is usually mutated. Using the integrated model, we correlated the extent of this uncoupling in TNBC cell lines with the importance of NRF2 in the 3D growth of these cell lines and their predicted handling of oxidative stress. Our results point to an oxidative stress tolerance network that is important for single cells during glandular development and the early stages of breast cancer.

A Live-Cell Screen for Altered Erk Dynamics Reveals Principles of Proliferative Control

Complex, time-varying responses have been observed widely in cell signaling, but how specific dynamics are generated or regulated is largely unknown. One major obstacle has been that high-throughput screens are typically incompatible with the live-cell assays used to monitor dynamics. Here, we address this challenge by screening a library of 429 kinase inhibitors and monitoring extracellular-regulated kinase (Erk) activity over 5 h in more than 80,000 single primary mouse keratinocytes. Our screen reveals both known and uncharacterized modulators of Erk dynamics, including inhibitors of non-epidermal growth factor receptor (EGFR) receptor tyrosine kinases (RTKs) that increase Erk pulse frequency and overall activity. Using drug treatment and direct optogenetic control, we demonstrate that drug-induced changes to Erk dynamics alter the conditions under which cells proliferate. Our work opens the door to high-throughput screens using live-cell biosensors and reveals that cell proliferation integrates information from Erk dynamics as well as additional permissive cues.

The receptor tyrosine kinase Ror is required for dendrite regeneration in Drosophila neurons

While many regulators of axon regeneration have been identified, very little is known about mechanisms that allow dendrites to regenerate after injury. Using a Drosophila model of dendrite regeneration, we performed a candidate screen of receptor tyrosine kinases (RTKs) and found a requirement for RTK-like orphan receptor (Ror). We confirmed that Ror was required for regeneration in two different neuron types using RNA interference (RNAi) and mutants. Ror was not required for axon regeneration or normal dendrite development, suggesting a specific role in dendrite regeneration. Ror can act as a Wnt coreceptor with frizzleds (fzs) in other contexts, so we tested the involvement of Wnt signaling proteins in dendrite regeneration. We found that knockdown of fz, dishevelled (dsh), Axin, and gilgamesh (gish) also reduced dendrite regeneration. Moreover, Ror was required to position dsh and Axin in dendrites. We recently found that Wnt signaling proteins, including dsh and Axin, localize microtubule nucleation machinery in dendrites. We therefore hypothesized that Ror may act by regulating microtubule nucleation at baseline and during dendrite regeneration. Consistent with this hypothesis, localization of the core nucleation protein γTubulin was reduced in Ror RNAi neurons, and this effect was strongest during dendrite regeneration. In addition, dendrite regeneration was sensitive to partial reduction of γTubulin. We conclude that Ror promotes dendrite regeneration as part of a Wnt signaling pathway that regulates dendritic microtubule nucleation.

SIN-KNO: A method of gene regulatory network inference using single-cell transcription and gene knockout data

As a tool of interpreting and analyzing genetic data, gene regulatory network (GRN) could reveal regulatory relationships between genes, proteins, and small molecules, as well as understand physiological activities and functions within biological cells, interact in pathways, and how to make changes in the organism. Traditional GRN research focuses on the analysis of the regulatory relationships through the average of cellular gene expressions. These methods are difficult to identify the cell heterogeneity of gene expression. Existing methods for inferring GRN using single-cell transcriptional data lack expression information when genes reach steady state, and the high dimensionality of single-cell data leads to high temporal and spatial complexity of the algorithm. In order to solve the problem in traditional GRN inference methods, including the lack of cellular heterogeneity information, single-cell data complexity and lack of steady-state information, we propose a method for GRN inference using single-cell transcription and gene knockout data, called SINgle-cell transcription data-KNOckout data (SIN-KNO), which focuses on combining dynamic and steady-state information of regulatory relationship contained in gene expression. Capturing cell heterogeneity information could help understand the gene expression difference in different cells. So, we could observe gene expression changes more accurately. Gene knockout data could observe the gene expression levels at steady-state of all other genes when one gene is knockout. Classifying the genes before analyzing the single-cell data could determine a large number of non-existent regulation, greatly reducing the number of regulation required for inference. In order to show the efficiency, the proposed method has been compared with several typical methods in this area including GENIE3, JUMP3, and SINCERITIES. The results of the evaluation indicate that the proposed method can analyze the diversified information contained in the two types of data, establish a more accurate gene regulation network, and improve the computational efficiency. The method provides a new thinking for dealing with large datasets and high computational complexity of single-cell data in the GRN inference.

Astrocytic trans-Differentiation Completes a Multicellular Paracrine Feedback Loop Required for Medulloblastoma Tumor Growth

The tumor microenvironment (TME) is critical for tumor progression. However, the establishment and function of the TME remain obscure because of its complex cellular composition. Using a mouse genetic system called mosaic analysis with double markers (MADMs), we delineated TME evolution at single-cell resolution in sonic hedgehog (SHH)-activated medulloblastomas that originate from unipotent granule neuron progenitors in the brain. First, we found that astrocytes within the TME (TuAstrocytes) were trans-differentiated from tumor granule neuron precursors (GNPs), which normally never differentiate into astrocytes. Second, we identified that TME-derived IGF1 promotes tumor progression. Third, we uncovered that insulin-like growth factor 1 (IGF1) is produced by tumor-associated microglia in response to interleukin-4 (IL-4) stimulation. Finally, we found that IL-4 is secreted by TuAstrocytes. Collectively, our studies reveal an evolutionary process that produces a multi-lateral network within the TME of medulloblastoma: a fraction of tumor cells trans-differentiate into TuAstrocytes, which, in turn, produce IL-4 that stimulates microglia to produce IGF1 to promote tumor progression.

The new technologies of high-throughput single-cell RNA sequencing

A wealth of genome and transcriptome data obtained using new generation sequencing (NGS) technologies for whole organisms could not answer many questions in oncology, immunology, physiology, neurobiology, zoology and other fields of science and medicine. Since the cell is the basis for the living of all unicellular and multicellular organisms, it is necessary to study the biological processes at its level. This understanding gave impetus to the development of a new direction – the creation of technologies that allow working with individual cells (single-cell technology). The rapid development of not only instruments, but also various advanced protocols for working with single cells is due to the relevance of these studies in many fields of science and medicine. Studying the features of various stages of ontogenesis, identifying patterns of cell differentiation and subsequent tissue development, conducting genomic and transcriptome analyses in various areas of medicine (especially in demand in immunology and oncology), identifying cell types and states, patterns of biochemical and physiological processes using single cell technologies, allows the comprehensive research to be conducted at a new level. The first RNA-sequencing technologies of individual cell transcriptomes (scRNA-seq) captured no more than one hundred cells at a time, which was insufficient due to the detection of high cell heterogeneity, existence of the minor cell types (which were not detected by morphology) and complex regulatory pathways. The unique techniques for isolating, capturing and sequencing transcripts of tens of thousands of cells at a time are evolving now. However, new technologies have certain differences both at the sample preparation stage and during the bioinformatics analysis. In the paper we consider the most effective methods of multiple parallel scRNA-seq using the example of 10XGenomics, as well as the specifics of such an experiment, further bioinformatics analysis of the data, future outlook and applications of new high-performance technologies.

Computational deconvolution of transcriptomics data from mixed cell populations

Summary

Gene expression analyses of bulk tissues often ignore cell type composition as an important confounding factor, resulting in a loss of signal from lowly abundant cell types. In this review, we highlight the importance and value of computational deconvolution methods to infer the abundance of different cell types and/or cell type-specific expression profiles in heterogeneous samples without performing physical cell sorting. We also explain the various deconvolution scenarios, the mathematical approaches used to solve them and the effect of data processing and different confounding factors on the accuracy of the deconvolution results.

Contact

katleen.depreter@ugent.be.

Supplementary information

Supplementary data are available at Bioinformatics online.

Clonal Isolation of Human Pluripotent Stem Cells on Nanofibrous Substrates Reveals an Advanced Subclone for Cardiomyocyte Differentiation

Human pluripotent stem cells (hPSCs) have been widely used for various applications including disease modeling and regenerative medicine, among others. Recently, an increasing number of studies has focused on heterogeneity among hPSCs, which could affect cell quality and subsequent applications. In this study, a nanofibrous platform is developed for single human induced pluripotent stem cell isolation and culture. One type of single cell-derived subclone is established and found to have a distinct morphology compared to other subclones. When used for differentiation toward cardiomyocytes, this type of subclone demonstrates higher differentiation efficiency, increased maturation, and stronger beating compared to those derived from the other subclones. The findings provide a convenient method for single-cell isolation and culture, and demonstrate that variations in differentiation tendencies exist among subclones from the same cell line. This substrate adhesion-based selection process could be used to obtain cell lines with improved differentiation efficiency toward cardiomyocytes and other cell types, which would be advantageous for studies in various fields.

GRAPHICAL MODELS FOR ZERO-INFLATED SINGLE CELL GENE EXPRESSION

Bulk gene expression experiments relied on aggregations of thousands of cells to measure the average expression in an organism. Advances in microfluidic and droplet sequencing now permit expression profiling in single cells. This study of cell-to-cell variation reveals that individual cells lack detectable expression of transcripts that appear abundant on a population level, giving rise to zero-inflated expression patterns. To infer gene co-regulatory networks from such data, we propose a multivariate Hurdle model. It is comprised of a mixture of singular Gaussian distributions. We employ neighborhood selection with the pseudo-likelihood and a group lasso penalty to select and fit undirected graphical models that capture conditional independences between genes. The proposed method is more sensitive than existing approaches in simulations, even under departures from our Hurdle model. The method is applied to data for T follicular helper cells, and a high-dimensional profile of mouse dendritic cells. It infers network structure not revealed by other methods; or in bulk data sets. An R implementation is available at https://github.com/amcdavid/HurdleNormal.

In situ 10-cell RNA sequencing in tissue and tumor biopsy samples

Single-cell transcriptomic methods classify new and existing cell types very effectively, but alternative approaches are needed to quantify the individual regulatory states of cells in their native tissue context. We combined the tissue preservation and single-cell resolution of laser capture with an improved preamplification procedure enabling RNA sequencing of 10 microdissected cells. This in situ 10-cell RNA sequencing (10cRNA-seq) can exploit fluorescent reporters of cell type in genetically engineered mice and is compatible with freshly cryoembedded clinical biopsies from patients. Through recombinant RNA spike-ins, we estimate dropout-free technical reliability as low as ~250 copies and a 50% detection sensitivity of ~45 copies per 10-cell reaction. By using small pools of microdissected cells, 10cRNA-seq improves technical per-cell reliability and sensitivity beyond existing approaches for single-cell RNA sequencing (scRNA-seq). Detection of low-abundance transcripts by 10cRNA-seq is comparable to random 10-cell groups of scRNA-seq data, suggesting no loss of gene recovery when cells are isolated in situ. Combined with existing approaches to deconvolve small pools of cells, 10cRNA-seq offers a reliable, unbiased, and sensitive way to measure cell-state heterogeneity in tissues and tumors.

A single-cell translocation and secretion assay (TransSeA)

Understanding biological heterogeneity at the single cell level is required for advancing insights into the complexity of human physiology and diseases. While advances in technological and analytical methods have afforded unprecedented glimpses of this heterogeneity, the information captured to date largely represents one-time "snap" shots of single cell physiology. To address the limits of existing methods and to accelerate discoveries from single cell studies, we developed a single-cell translocation and secretion assay (TransSeA) that supports time lapse analysis, enables molecular cargo analysis of secretions such as extracellular vesicles (EVs) from single cells, allows massively parallel single cell transfer according to user-defined cell selection criteria, and supports tracking of phenotypes between parental and progeny cells derived from single cells. To demonstrate the unique capabilities and efficiencies of the assay, we present unprecedented single cell studies related to cell secretions, EV cargos and cell intrinsic properties. Although used as examples to demonstrate the feasibility and versatility of the technology, the studies already provided insights into key unanswered questions such as the microRNAs carried by EVs, the relationships between EV secretion rates and gene expressions, and the spontaneous, trans-generational phenotypic changes in EV secretion between parental and progeny cells.

Harnessing Tumor Evolution to Circumvent Resistance

High-throughput sequencing can be used to measure changes in tumor composition across space and time. Specifically, comparisons of pre- and post-treatment samples can reveal the underlying clonal dynamics and resistance mechanisms. Here, we discuss evidence for distinct modes of tumor evolution and their implications for therapeutic strategies. In addition, we consider the utility of spatial tissue sampling schemes, single-cell analysis, and circulating tumor DNA to track tumor evolution and the emergence of resistance, as well as approaches that seek to forestall resistance by targeting tumor evolution. Ultimately, characterization of the (epi)genomic, transcriptomic, and phenotypic changes that occur during tumor progression coupled with computational and mathematical modeling of tumor evolutionary dynamics may inform personalized treatment strategies.

Automated brightfield morphometry of 3D organoid populations by OrganoSeg

Spheroid and organoid cultures are powerful in vitro models for biology, but size and shape diversity within the culture is largely ignored. To streamline morphometric profiling, we developed OrganoSeg, an open-source software that integrates segmentation, filtering, and analysis for archived brightfield images of 3D culture. OrganoSeg is more accurate and flexible than existing platforms, and we illustrate its potential by stratifying 5167 breast-cancer spheroid and 5743 colon and colorectal-cancer organoid morphologies. Organoid transcripts grouped by morphometric signature heterogeneity were enriched for biological processes not prominent in the original RNA sequencing data. OrganoSeg enables complete, objective quantification of brightfield phenotypes, which may give insight into the molecular and multicellular mechanisms of organoid regulation.

Tumor-Suppressor Inactivation of GDF11 Occurs by Precursor Sequestration in Triple-Negative Breast Cancer

Triple-negative breast cancer (TNBC) is an aggressive and heterogeneous carcinoma in which various tumor-suppressor genes are lost by mutation, deletion, or silencing. Here we report a tumor-suppressive mode of action for growth-differentiation factor 11 (GDF11) and an unusual mechanism of its inactivation in TNBC. GDF11 promotes an epithelial, anti-invasive phenotype in 3D triple-negative cultures and intraductal xenografts by sustaining expression of E-cadherin and inhibitor of differentiation 2 (ID2). Surprisingly, clinical TNBCs retain the GDF11 locus and expression of the protein itself. GDF11 bioactivity is instead lost because of deficiencies in its convertase, proprotein convertase subtilisin/kexin type 5 (PCSK5), causing inactive GDF11 precursor to accumulate intracellularly. PCSK5 reconstitution mobilizes the latent TNBC reservoir of GDF11 in vitro and suppresses triple-negative mammary cancer metastasis to the lung of syngeneic hosts. Intracellular GDF11 retention adds to the concept of tumor-suppressor inactivation and reveals a cell-biological vulnerability for TNBCs lacking therapeutically actionable mutations.

Akt regulation of glycolysis mediates bioenergetic stability in epithelial cells

Cells use multiple feedback controls to regulate metabolism in response to nutrient and signaling inputs. However, feedback creates the potential for unstable network responses. We examined how concentrations of key metabolites and signaling pathways interact to maintain homeostasis in proliferating human cells, using fluorescent reporters for AMPK activity, Akt activity, and cytosolic NADH/NAD⁺ redox. Across various conditions, including glycolytic or mitochondrial inhibition or cell proliferation, we observed distinct patterns of AMPK activity, including both stable adaptation and highly dynamic behaviors such as periodic oscillations and irregular fluctuations that indicate a failure to reach a steady state. Fluctuations in AMPK activity, Akt activity, and cytosolic NADH/NAD⁺ redox state were temporally linked in individual cells adapting to metabolic perturbations. By monitoring single-cell dynamics in each of these contexts, we identified PI3K/Akt regulation of glycolysis as a multifaceted modulator of single-cell metabolic dynamics that is required to maintain metabolic stability in proliferating cells.

Linear Integration of ERK Activity Predominates over Persistence Detection in Fra-1 Regulation

ERK signaling regulates the expression of target genes, but it is unclear how ERK activity dynamics are interpreted. Here, we investigate this question using simultaneous, live, single-cell imaging of two ERK activity reporters and expression of Fra-1, a target gene controlling epithelial cell identity. We find that Fra-1 is expressed in proportion to the amplitude and duration of ERK activity. In contrast to previous "persistence detector" and "selective filter" models in which Fra-1 expression only occurs when ERK activity persists beyond a threshold duration, our observations demonstrate that the network regulating Fra-1 expression integrates total ERK activity and responds to it linearly. However, exploration of a generalized mathematical model of the Fra-1 coherent feedforward loop demonstrates that it can perform either linear integration or persistence detection, depending on the basal mRNA production rate and protein production delays. Our data indicate that significant basal expression and short delays cause Fra-1 to respond linearly to integrated ERK activity.

Translation initiation events on structured eukaryotic mRNAs generate gene expression noise

Gene expression stochasticity plays a major role in biology, creating non-genetic cellular individuality and influencing multiple processes, including differentiation and stress responses. We have addressed the lack of knowledge about posttranscriptional contributions to noise by determining cell-to-cell variations in the abundance of mRNA and reporter protein in yeast. Two types of structural element, a stem–loop and a poly(G) motif, not only inhibit translation initiation when inserted into an mRNA 5΄ untranslated region, but also generate noise. The noise-enhancing effect of the stem–loop structure also remains operational when combined with an upstream open reading frame. This has broad significance, since these elements are known to modulate the expression of a diversity of eukaryotic genes. Our findings suggest a mechanism for posttranscriptional noise generation that will contribute to understanding of the generally poor correlation between protein-level stochasticity and transcriptional bursting. We propose that posttranscriptional stochasticity can be linked to cycles of folding/unfolding of a stem–loop structure, or to interconversion between higher-order structural conformations of a G-rich motif, and have created a correspondingly configured computational model that generates fits to the experimental data. Stochastic events occurring during the ribosomal scanning process can therefore feature alongside transcriptional bursting as a source of noise.

Single-Cell Sequencing Technologies for Cardiac Stem Cell Studies

Today with the rapid advancements in stem cell studies and the promising potential of using stem cells in clinical therapy, there is an increasing demand for in-depth comprehensive analysis on individual cell transcriptome and epigenome, as they play critical roles in a number of cell functions such as cell differentiation, growth, and reprogramming. The development of single-cell sequencing technologies has helped in revealing some exciting new perspectives in stem cells and regenerative medicine research. Among the various potential applications, single-cell analysis for cardiac stem cells (CSCs) holds tremendous promises in understanding the mechanisms of heart development and regeneration, which might light up the path toward cell therapy for cardiovascular diseases. This review briefly highlights the recent progresses in single-cell sequencing analysis technologies and their applications in CSC research.

Fibroblasts: Diverse Cells Critical to Biomaterials Integration

Fibroblasts are key participants in wound healing and inflammation, and are capable of driving the progression of tissue repair to fully functional tissue or pathologic scar, or fibrosis, depending on the specific mechanical and biochemical cues with which they are presented. Thus, understanding and modulating the fibroblastic response to implanted materials is paramount to achieving desirable outcomes, such as long-term implant function or tissue regeneration. However, fibroblasts are remarkably heterogeneous and can differ vastly in their contributions to regeneration and fibrosis. This heterogeneity exists between tissues and within tissues, down to the level of individual cells. This review will discuss the role of fibroblasts, the pitfalls of describing them as a collective, the specifics of their function, and potential future directions to better understand and organize their highly variable biology.

Inhibition of extracellular matrix mediated TGF-β signalling suppresses endometrial cancer metastasis

Although aggressive invasion and distant metastases are an important cause of morbidity and mortality in patients with endometrial cancer (EC), the requisite events determining this propensity are currently unknown. Using organotypic three-dimensional culture of endometrial cancer cell lines, we demonstrated anti-correlated TGF-β signalling gene expression patterns that arise among extracellular matrix (ECM)-attached cells. TGF-β pathway seemed to be active in EC cells forming non-glandular colonies in 3D-matrix but weaker in glandular colonies. Functionally we found that out of several ECM proteins, fibronectin relatively promotes Smad phosphorylation suggesting a potential role in regulating TGF-β signalling in non-glandular colonies. Importantly, alteration of TGF-β pathway induced EMT and MET in both type of colonies through slug protein. The results exemplify a crucial role of TGF-β pathway during EC metastasis in human patients and inhibition of the pathway in a murine model impaired tumour cell invasion and metastasis depicting an attractive target for therapeutic intervention of malignant tumour progression. These findings provide key insights into the role of ECM-derived TGF-β signalling to promote endometrial cancer metastasis and offer an avenue for therapeutic targeting of microenvironment derived signals along with tumour cells.

Environment Tunes Propagation of Cell-to-Cell Variation in the Human Macrophage Gene Network

Cell-to-cell variation in gene expression and the propagation of such variation (PoV or "noise propagation") from one gene to another in the gene network, as reflected by gene-gene correlation across single cells, are commonly observed in single-cell transcriptomic studies and can shape the phenotypic diversity of cell populations. While gene network "rewiring" is known to accompany cellular adaptation to different environments, how PoV changes between environments and its underlying regulatory mechanisms are less understood. Here, we systematically explored context-dependent PoV among genes in human macrophages, utilizing different cytokines as natural perturbations of multiple molecular parameters that may influence PoV. Our single-cell, epigenomic, computational, and stochastic simulation analyses reveal that environmental adaptation can tune PoV to potentially shape cellular heterogeneity by changing parameters such as the degree of phosphorylation and transcription factor-chromatin interactions. This quantitative tuning of PoV may be a widespread, yet underexplored, property of cellular adaptation to distinct environments.

Profiling Subcellular Protein Phosphatase Responses to Coxsackievirus B3 Infection of Cardiomyocytes

Cellular responses to stimuli involve dynamic and localized changes in protein kinases and phosphatases. Here, we report a generalized functional assay for high-throughput profiling of multiple protein phosphatases with subcellular resolution and apply it to analyze coxsackievirus B3 (CVB3) infection counteracted by interferon signaling. Using on-plate cell fractionation optimized for adherent cells, we isolate protein extracts containing active endogenous phosphatases from cell membranes, the cytoplasm, and the nucleus. The extracts contain all major classes of protein phosphatases and catalyze dephosphorylation of plate-bound phosphosubstrates in a microtiter format, with cellular activity quantified at the end point by phosphospecific ELISA. The platform is optimized for six phosphosubstrates (ERK2, JNK1, p38α, MK2, CREB, and STAT1) and measures specific activities from extracts of fewer than 50,000 cells. The assay was exploited to examine viral and antiviral signaling in AC16 cardiomyocytes, which we show can be engineered to serve as susceptible and permissive hosts for CVB3. Phosphatase responses were profiled in these cells by completing a full-factorial experiment for CVB3 infection and type I/II interferon signaling. Over 850 functional measurements revealed several independent, subcellular changes in specific phosphatase activities. During CVB3 infection, we found that type I interferon signaling increases subcellular JNK1 phosphatase activity, inhibiting nuclear JNK1 activity that otherwise promotes viral protein synthesis in the infected host cell. Our assay provides a high-throughput way to capture perturbations in important negative regulators of intracellular signal-transduction networks.

Single-cell transcriptome and epigenomic reprogramming of cardiomyocyte-derived cardiac progenitor cells

The molecular basis underlying the dedifferentiation of mammalian adult cardiomyocytes (ACMs) into myocyte-derived cardiac progenitor cells (mCPCs) during cardiac tissue regeneration is poorly understood. We present data integrating single-cell transcriptome and whole-genome DNA methylome analyses of mouse mCPCs to understand the epigenomic reprogramming governing their intrinsic cellular plasticity. Compared to parental cardiomyocytes, mCPCs display epigenomic reprogramming with many differentially-methylated regions, both hypermethylated and hypomethylated, across the entire genome. Correlating well with the methylome, our single-cell transcriptomic data show that the genes encoding cardiac structure and function proteins are remarkably down-regulated in mCPCs, while those for cell cycle, proliferation, and stemness are significantly up-regulated. In addition, implanting mCPCs into infarcted mouse myocardium improves cardiac function with augmented left ventricular ejection fraction. This dataset suggests that the cellular plasticity of mammalian cardiomyocytes is the result of a well-orchestrated epigenomic reprogramming and a subsequent global transcriptomic alteration. Understanding cardiomyocyte epigenomic reprogramming may enable the design of future clinical therapies that induce cardiac regeneration, and prevent heart failure.

End Sequence Analysis Toolkit (ESAT) expands the extractable information from single-cell RNA-seq data

RNA-seq protocols that focus on transcript termini are well suited for applications in which template quantity is limiting. Here we show that, when applied to end-sequencing data, analytical methods designed for global RNA-seq produce computational artifacts. To remedy this, we created the End Sequence Analysis Toolkit (ESAT). As a test, we first compared end-sequencing and bulk RNA-seq using RNA from dendritic cells stimulated with lipopolysaccharide (LPS). As predicted by the telescripting model for transcriptional bursts, ESAT detected an LPS-stimulated shift to shorter 3′-isoforms that was not evident by conventional computational methods. Then, droplet-based microfluidics was used to generate 1000 cDNA libraries, each from an individual pancreatic islet cell. ESAT identified nine distinct cell types, three distinct β-cell types, and a complex interplay between hormone secretion and vascularization. ESAT, then, offers a much-needed and generally applicable computational pipeline for either bulk or single-cell RNA end-sequencing.

Robust Inference of Cell-to-Cell Expression Variations from Single- and K-Cell Profiling

Quantifying heterogeneity in gene expression among single cells can reveal information inaccessible to cell-population averaged measurements. However, the expression level of many genes in single cells fall below the detection limit of even the most sensitive technologies currently available. One proposed approach to overcome this challenge is to measure random pools of k cells (e.g., 10) to increase sensitivity, followed by computational “deconvolution” of cellular heterogeneity parameters (CHPs), such as the biological variance of single-cell expression levels. Existing approaches infer CHPs using either single-cell or k-cell data alone, and typically within a single population of cells. However, integrating both single- and k-cell data may reap additional benefits, and quantifying differences in CHPs across cell populations or conditions could reveal novel biological information. Here we present a Bayesian approach that can utilize single-cell, k-cell, or both simultaneously to infer CHPs within a single condition or their differences across two conditions. Using simulated as well as experimentally generated single- and k-cell data, we found situations where each data type would offer advantages, but using both together can improve precision and better reconcile CHP information contained in single- and k-cell data. We illustrate the utility of our approach by applying it to jointly generated single- and k-cell data to reveal CHP differences in several key inflammatory genes between resting and inflammatory cytokine-activated human macrophages, delineating differences in the distribution of ‘ON’ versus ‘OFF’ cells and in continuous variation of expression level among cells. Our approach thus offers a practical and robust framework to assess and compare cellular heterogeneity within and across biological conditions using modern multiplexed technologies.

Different cells can make different amounts of biomolecules such as RNA transcripts of genes. New technologies are emerging to measure the transcript level of many genes in single cells. However, accurate quantification of the biological variation from cell to cell can be challenging due to the low transcript level of many genes and the presence of substantial measurement noise. Here we present a flexible, novel computational approach to quantify biological cell-to-cell variation that can use different types of data, namely measurements directly obtained from single cells, and/or those from random pools of k-cells (e.g., k = 10). Assessment of these different inputs using simulated and real data revealed that each data type can offer advantages under different scenarios, but combining both single- and k-cell measurements tend to offer the best of both. Application of our approach to single- and k-cell data obtained from resting and inflammatory macrophages, an important type of immune cells implicated in diverse diseases, revealed interesting changes in cell-to-cell variation in transcript levels upon inflammatory stimulation, thus suggesting that inflammation can shape not only the average expression level of a gene but also the gene’s degree of expression variation among single cells.

TNF-insulin crosstalk at the transcription factor GATA6 is revealed by a model that links signaling and transcriptomic data tensors

Signal transduction networks coordinate transcriptional programs activated by diverse extracellular stimuli, such as growth factors and cytokines. Cells receive multiple stimuli simultaneously, and mapping how activation of the integrated signaling network affects gene expression is a challenge. We stimulated colon adenocarcinoma cells with various combinations of the cytokine tumor necrosis factor (TNF) and the growth factors insulin and epidermal growth factor (EGF) to investigate signal integration and transcriptional crosstalk. We quantitatively linked the proteomic and transcriptomic data sets by implementing a structured computational approach called tensor partial least squares regression. This statistical model accurately predicted transcriptional signatures from signaling arising from single and combined stimuli and also predicted time-dependent contributions of signaling events. Specifically, the model predicted that an early-phase, AKT-associated signal downstream of insulin repressed a set of transcripts induced by TNF. Through bioinformatics and cell-based experiments, we identified the AKT-repressed signal as glycogen synthase kinase 3 (GSK3)-catalyzed phosphorylation of Ser(37) on the long form of the transcription factor GATA6. Phosphorylation of GATA6 on Ser(37) promoted its degradation, thereby preventing GATA6 from repressing transcripts that are induced by TNF and attenuated by insulin. Our analysis showed that predictive tensor modeling of proteomic and transcriptomic data sets can uncover pathway crosstalk that produces specific patterns of gene expression in cells receiving multiple stimuli.

Cell type-specific transcriptome profiling in mammalian brains

A mammalian brain contains numerous types of cells. Advances in neuroscience in the past decade allow us to identify and isolate neural cells of interest from mammalian brains. Recent developments in high-throughput technologies, such as microarrays and next-generation sequencing (NGS), provide detailed information on gene expression in pooled cells on a genomic scale. As a result, many novel genes have been found critical in cell type-specific transcriptional regulation. These differentially expressed genes can be used as molecular signatures, unique to a particular class of neural cells. Use of this gene expression-based approach can further differentiate neural cell types into subtypes, potentially linking some of them with neurological diseases. In this article, experimental techniques used to purify neural cells are described, followed by a review on recent microarray- or NGS-based transcriptomic studies of common neural cell types. The future prospects of cell type-specific research are also discussed.

Small Molecule Inhibitor of CBFβ-RUNX Binding for RUNX Transcription Factor Driven Cancers

Transcription factors have traditionally been viewed with skepticism as viable drug targets, but they offer the potential for completely novel mechanisms of action that could more effectively address the stem cell like properties, such as self-renewal and chemo-resistance, that lead to the failure of traditional chemotherapy approaches. Core binding factor is a heterodimeric transcription factor comprised of one of 3 RUNX proteins (RUNX1-3) and a CBFβ binding partner. CBFβ enhances DNA binding of RUNX subunits by relieving auto-inhibition. Both RUNX1 and CBFβ are frequently mutated in human leukemia. More recently, RUNX proteins have been shown to be key players in epithelial cancers, suggesting the targeting of this pathway could have broad utility. In order to test this, we developed small molecules which bind to CBFβ and inhibit its binding to RUNX. Treatment with these inhibitors reduces binding of RUNX1 to target genes, alters the expression of RUNX1 target genes, and impacts cell survival and differentiation. These inhibitors show efficacy against leukemia cells as well as basal-like (triple-negative) breast cancer cells. These inhibitors provide effective tools to probe the utility of targeting RUNX transcription factor function in other cancers.

Multiscale analysis of the murine intestine for modeling human diseases

When functioning properly, the intestine is one of the key interfaces between the human body and its environment. It is responsible for extracting nutrients from our food and excreting our waste products. It provides an environment for a host of healthful microbes and serves as a first defense against pathogenic ones. These processes require tight homeostatic controls, which are provided by the interactions of a complex mix of epithelial, stromal, neural and immune cells, as well as the resident microflora. This homeostasis can be disrupted by invasive microbes, genetic lesions, and carcinogens, resulting in diseases such Clostridium difficile infection, inflammatory bowel disease (IBD) and cancer. Enormous strides have been made in understanding how this important organ functions in health and disease using everything from cell culture systems to animal models to human tissue samples. This has resulted in better therapies for all of these diseases, but there is still significant room for improvement. In the United States alone, 14,000 people per year die of C. difficile, up to 1.6 million people suffer from IBD, and more than 50,000 people die every year from colon cancer. Because these and other intestinal diseases arise from complex interactions between the different components of the gut ecosystem, we propose that systems approaches that address this complexity in an integrative manner may eventually lead to improved therapeutics that deliver lasting cures. This review will discuss the use of systems biology for studying intestinal diseases in vivo with particular emphasis on mouse models. Additionally, it will focus on established experimental techniques that have been used to drive this systems-level analysis, and emerging techniques that will push this field forward in the future.

Single-cell states versus single-cell atlases - two classes of heterogeneity that differ in meaning and method

Recent advances have created new opportunities to dissect cellular heterogeneity at the omics level. The enthusiasm for deep single-cell profiling has obscured a discussion of different types of heterogeneity and the most-appropriate techniques for studying each type. Here, I distinguish heterogeneity in regulation from heterogeneity in lineage. Snapshots of lineage heterogeneity provide a cell atlas that catalogs cellular diversity within complex tissues. Profiles of regulatory heterogeneity seek to interrogate one lineage deeply to capture an ensemble of single-cell states. Single-cell atlases require molecular signatures from many cells at a throughput afforded by mass cytometry-based, microfluidic-based, and microencapsulation-based methods. Single-cell states are more dependent on time, microenvironment, and low-abundance transcripts, emphasizing in situ methods that stress depth of profiling and quantitative accuracy.

Single-cell triple omics sequencing reveals genetic, epigenetic, and transcriptomic heterogeneity in hepatocellular carcinomas

Single-cell genome, DNA methylome, and transcriptome sequencing methods have been separately developed. However, to accurately analyze the mechanism by which transcriptome, genome and DNA methylome regulate each other, these omic methods need to be performed in the same single cell. Here we demonstrate a single-cell triple omics sequencing technique, scTrio-seq, that can be used to simultaneously analyze the genomic copy-number variations (CNVs), DNA methylome, and transcriptome of an individual mammalian cell. We show that large-scale CNVs cause proportional changes in RNA expression of genes within the gained or lost genomic regions, whereas these CNVs generally do not affect DNA methylation in these regions. Furthermore, we applied scTrio-seq to 25 single cancer cells derived from a human hepatocellular carcinoma tissue sample. We identified two subpopulations within these cells based on CNVs, DNA methylome, or transcriptome of individual cells. Our work offers a new avenue of dissecting the complex contribution of genomic and epigenomic heterogeneities to the transcriptomic heterogeneity within a population of cells.

Differentiation of human gingival mesenchymal stem cells into neuronal lineages in 3D bioconjugated injectable protein hydrogel construct for the management of neuronal disorder

The success of regeneration attempt is based on an ideal combination of stem cells, scaffolding and growth factors. Tissue constructs help to maintain stem cells in a required area for a desired time. There is a need for easily obtainable cells, potentially autologous stem cells and a biologically acceptable scaffold for use in humans in different difficult situations. This study aims to address these issues utilizing a unique combination of stem cells from gingiva and a hydrogel scaffold, based on a natural product for regenerative application. Human gingival mesenchymal stem cells (HGMSCs) were, with due induction, differentiated to neuronal lineages to overcome the problems associated with birth tissue-related stem cells. The differentiation potential of neuronal lineages was confirmed with suitable specific markers. The properties of mesenchymal stem cells in encapsulated form were observed to be similar to free cells. The encapsulated cells (3D) were then subjected to differentiation into neuronal lineages with suitable inducers, and the morphology and gene expression of transient cells were analyzed. HGMSCs was differentiated into neuronal lineages as both free and encapsulated forms without any significant differences. The presence of Nissl bodies and the neurite outgrowth confirm the differentiation. The advantages of this new combination appear to make it a promising tissue construct for translational application.

Epigenomic Reprogramming of Adult Cardiomyocyte-Derived Cardiac Progenitor Cells

It has been believed that mammalian adult cardiomyocytes (ACMs) are terminally-differentiated and are unable to proliferate. Recently, using a bi-transgenic ACM fate mapping mouse model and an in vitro culture system, we demonstrated that adult mouse cardiomyocytes were able to dedifferentiate into cardiac progenitor-like cells (CPCs). However, little is known about the molecular basis of their intrinsic cellular plasticity. Here we integrate single-cell transcriptome and whole-genome DNA methylation analyses to unravel the molecular mechanisms underlying the dedifferentiation and cell cycle reentry of mouse ACMs. Compared to parental cardiomyocytes, dedifferentiated mouse cardiomyocyte-derived CPCs (mCPCs) display epigenomic reprogramming with many differentially-methylated regions, both hypermethylated and hypomethylated, across the entire genome. Correlated well with the methylome, our transcriptomic data showed that the genes encoding cardiac structure and function proteins are remarkably down-regulated in mCPCs, while those for cell cycle, proliferation, and stemness are significantly up-regulated. In addition, implantation of mCPCs into infarcted mouse myocardium improves cardiac function with augmented left ventricular ejection fraction. Our study demonstrates that the cellular plasticity of mammalian cardiomyocytes is the result of a well-orchestrated epigenomic reprogramming and a subsequent global transcriptomic alteration.

Receptor Level Mechanisms Are Required for Epidermal Growth Factor (EGF)-stimulated Extracellular Signal-regulated Kinase (ERK) Activity Pulses

In both physiological and cell culture systems, EGF-stimulated ERK activity occurs in discrete pulses within individual cells. Many feedback loops are present in the EGF receptor (EGFR)-ERK network, but the mechanisms driving pulsatile ERK kinetics are unknown. Here, we find that in cells that respond to EGF with frequency-modulated pulsatile ERK activity, stimulation through a heterologous TrkA receptor system results in non-pulsatile, amplitude-modulated activation of ERK. We further dissect the kinetics of pulse activity using a combination of FRET- and translocation-based reporters and find that EGFR activity is required to maintain ERK activity throughout the 10-20-minute lifetime of pulses. Together, these data indicate that feedbacks operating within the core Ras-Raf-MEK-ERK cascade are insufficient to drive discrete pulses of ERK activity and instead implicate mechanisms acting at the level of EGFR.