A developmental constraint model of cancer cell states and tumor heterogeneity.

Cancer is a disease that stems from a fundamental liability inherent to multicellular life forms in which an individual cell is capable of reneging on the interests of the collective organism. Although cancer is commonly described as an evolutionary process, a less appreciated aspect of tumorigenesis may be the constraints imposed by the organism's developmental programs. Recent work from single-cell transcriptomic analyses across a range of cancer types has revealed the recurrence, plasticity, and co-option of distinct cellular states among cancer cell populations. Here, we note that across diverse cancer types, the observed cell states are proximate within the developmental hierarchy of the cell of origin. We thus posit a model by which cancer cell states are directly constrained by the organism's "developmental map." According to this model, a population of cancer cells traverses the developmental map, thereby generating a heterogeneous set of states whose interactions underpin emergent tumor behavior.

Human microglial state dynamics in Alzheimer's disease progression.

Altered microglial states affect neuroinflammation, neurodegeneration, and disease but remain poorly understood. Here, we report 194,000 single-nucleus microglial transcriptomes and epigenomes across 443 human subjects and diverse Alzheimer's disease (AD) pathological phenotypes. We annotate 12 microglial transcriptional states, including AD-dysregulated homeostatic, inflammatory, and lipid-processing states. We identify 1,542 AD-differentially-expressed genes, including both microglia-state-specific and disease-stage-specific alterations. By integrating epigenomic, transcriptomic, and motif information, we infer upstream regulators of microglial cell states, gene-regulatory networks, enhancer-gene links, and transcription-factor-driven microglial state transitions. We demonstrate that ectopic expression of our predicted homeostatic-state activators induces homeostatic features in human iPSC-derived microglia-like cells, while inhibiting activators of inflammation can block inflammatory progression. Lastly, we pinpoint the expression of AD-risk genes in microglial states and differential expression of AD-risk genes and their regulators during AD progression. Overall, we provide insights underlying microglial states, including state-specific and AD-stage-specific microglial alterations at unprecedented resolution.

Atlas of clinically distinct cell states and ecosystems across human solid tumors.

Determining how cells vary with their local signaling environment and organize into distinct cellular communities is critical for understanding processes as diverse as development, aging, and cancer. Here we introduce EcoTyper, a machine learning framework for large-scale identification and validation of cell states and multicellular communities from bulk, single-cell, and spatially resolved gene expression data. When applied to 12 major cell lineages across 16 types of human carcinoma, EcoTyper identified 69 transcriptionally defined cell states. Most states were specific to neoplastic tissue, ubiquitous across tumor types, and significantly prognostic. By analyzing cell-state co-occurrence patterns, we discovered ten clinically distinct multicellular communities with unexpectedly strong conservation, including three with myeloid and stromal elements linked to adverse survival, one enriched in normal tissue, and two associated with early cancer development. This study elucidates fundamental units of cellular organization in human carcinoma and provides a framework for large-scale profiling of cellular ecosystems in any tissue.

Multi-omics analysis of innate and adaptive responses to BCG vaccination reveals epigenetic cell states that predict trained immunity.

Immune responses are tightly regulated yet highly variable between individuals. To investigate human population variation of trained immunity, we immunized healthy individuals with Bacillus Calmette-Guérin (BCG). This live-attenuated vaccine induces not only an adaptive immune response against tuberculosis but also triggers innate immune activation and memory that are indicative of trained immunity. We established personal immune profiles and chromatin accessibility maps over a 90-day time course of BCG vaccination in 323 individuals. Our analysis uncovered genetic and epigenetic predictors of baseline immunity and immune response. BCG vaccination enhanced the innate immune response specifically in individuals with a dormant immune state at baseline, rather than providing a general boost of innate immunity. This study advances our understanding of BCG's heterologous immune-stimulatory effects and trained immunity in humans. Furthermore, it highlights the value of epigenetic cell states for connecting immune function with genotype and the environment.

CRISPR Tools for Physiology and Cell State Changes: Potential of Transcriptional Engineering and Epigenome Editing.

Given the large amount of genome-wide data that have been collected during the last decades, a good understanding of how and why cells change during development, homeostasis, and disease might be expected. Unfortunately, the opposite is true; triggers that cause cellular state changes remain elusive, and the underlying molecular mechanisms are poorly understood. Although genes with the potential to influence cell states are known, the historic dependency on methods that manipulate gene expression outside the endogenous chromatin context has prevented us from understanding how cells organize, interpret, and protect cellular programs. Fortunately, recent methodological innovations are now providing options to answer these outstanding questions, by allowing to target and manipulate individual genomic and epigenomic loci. In particular, three experimental approaches are now feasible due to DNA targeting tools, namely, activation and/or repression of master transcription factors in their endogenous chromatin context; targeting transcription factors to endogenous, alternative, or inaccessible sites; and finally, functional manipulation of the chromatin context. In this article, we discuss the molecular basis of DNA targeting tools and review the potential of these new technologies before we summarize how these have already been used for the manipulation of cellular states and hypothesize about future applications.

Early-stage lung cancer is driven by a transitional cell state dependent on a KRAS-ITGA3-SRC axis.

Glycine-12 mutations in the GTPase KRAS (KRASG12) are an initiating event for development of lung adenocarcinoma (LUAD). KRASG12 mutations promote cell-intrinsic rewiring of alveolar type-II progenitor (AT2) cells, but to what extent such changes interplay with lung homeostasis and cell fate pathways is unclear. Here, we generated single-cell RNA-seq (scRNA-seq) profiles from AT2-mesenchyme organoid co-cultures, mice, and stage-IA LUAD patients, identifying conserved regulators of AT2 transcriptional dynamics and defining the impact of KRASG12D mutation with temporal resolution. In AT2WT organoids, we found a transient injury/plasticity state preceding AT2 self-renewal and AT1 differentiation. Early-stage AT2KRAS cells exhibited perturbed gene expression dynamics, most notably retention of the injury/plasticity state. The injury state in AT2KRAS cells of patients, mice, and organoids was distinguishable from AT2WT states via altered receptor expression, including co-expression of ITGA3 and SRC. The combination of clinically relevant KRASG12D and SRC inhibitors impaired AT2KRAS organoid growth. Together, our data show that an injury/plasticity state essential for lung repair is co-opted during AT2 self-renewal and LUAD initiation, suggesting that early-stage LUAD may be susceptible to interventions that target specifically the oncogenic nature of this cell state.

Single-cell analysis of treatment-resistant prostate cancer: Implications of cell state changes for cell surface antigen-targeted therapies.

Targeting cell surface molecules using radioligand and antibody-based therapies has yielded considerable success across cancers. However, it remains unclear how the expression of putative lineage markers, particularly cell surface molecules, varies in the process of lineage plasticity, wherein tumor cells alter their identity and acquire new oncogenic properties. A notable example of lineage plasticity is the transformation of prostate adenocarcinoma (PRAD) to neuroendocrine prostate cancer (NEPC)-a growing resistance mechanism that results in the loss of responsiveness to androgen blockade and portends dismal patient survival. To understand how lineage markers vary across the evolution of lineage plasticity in prostate cancer, we applied single-cell analyses to 21 human prostate tumor biopsies and two genetically engineered mouse models, together with tissue microarray analysis on 131 tumor samples. Not only did we observe a higher degree of phenotypic heterogeneity in castrate-resistant PRAD and NEPC than previously anticipated but also found that the expression of molecules targeted therapeutically, namely PSMA, STEAP1, STEAP2, TROP2, CEACAM5, and DLL3, varied within a subset of gene-regulatory networks (GRNs). We also noted that NEPC and small cell lung cancer subtypes shared a set of GRNs, indicative of conserved biologic pathways that may be exploited therapeutically across tumor types. While this extreme level of transcriptional heterogeneity, particularly in cell surface marker expression, may mitigate the durability of clinical responses to current and future antigen-directed therapies, its delineation may yield signatures for patient selection in clinical trials, potentially across distinct cancer types.

Cell state transitions: definitions and challenges.

A fundamental challenge when studying biological systems is the description of cell state dynamics. During transitions between cell states, a multitude of parameters may change - from the promoters that are active, to the RNAs and proteins that are expressed and modified. Cells can also adopt different shapes, alter their motility and change their reliance on cell-cell junctions or adhesion. These parameters are integral to how a cell behaves and collectively define the state a cell is in. Yet, technical challenges prevent us from measuring all of these parameters simultaneously and dynamically. How, then, can we comprehend cell state transitions using finite descriptions? The recent virtual workshop organised by The Company of Biologists entitled 'Cell State Transitions: Approaches, Experimental Systems and Models' attempted to address this question. Here, we summarise some of the main points that emerged during the workshop's themed discussions. We also present examples of cell state transitions and describe models and systems that are pushing forward our understanding of how cells rewire their state.

The Heterogeneous Cellular States of Glioblastoma Stem Cells Revealed by Single Cell Analysis.

Glioblastoma stem cells (GSCs) contributed to the progression, treatment resistance, and relapse of glioblastoma (GBM). However, current researches on GSCs were performed usually outside the human tumor microenvironment, ignoring the importance of the cellular states of primary GSCs. In this study, we leveraged single-cell transcriptome sequencing data of 6 independent GBM cohorts from public databases, and combined lineage and stemness features to identify primary GSCs. We dissected the cell states of GSCs and correlated them with the clinical outcomes of patients. As a result, we constructed a cellular hierarchy where GSCs resided at the center. In addition, we identified and characterized 2 different and recurrent GSCs subpopulations: proliferative GSCs (pGSCs) and quiescent GSCs (qGSCs). The pGSCs showed high cell cycle activity, indicating rapid cell division, while qGSCs showed a quiescent state. Then we traced the processes of tumor development by pseudo-time analysis and tumor phylogeny, and found that GSCs accumulated throughout the whole tumor development period. During the process, pGSCs mainly contributed to the early stage and qGSCs were enriched in the later stage. Finally, we constructed an 8-gene prognostic signature reflecting pGSCs activity and found that patients whose tumors were enriched for the pGSC signature had poor clinical outcomes. Our study highlights the primary GSCs heterogeneity and its correlation to tumor development and clinical outcomes, providing the potential targets for GBM treatment.

scEvoNet: a gradient boosting-based method for prediction of cell state evolution.

BACKGROUND: Exploring the function or the developmental history of cells in various organisms provides insights into a given cell type's core molecular characteristics and putative evolutionary mechanisms. Numerous computational methods now exist for analyzing single-cell data and identifying cell states. These methods mostly rely on the expression of genes considered as markers for a given cell state. Yet, there is a lack of scRNA-seq computational tools to study the evolution of cell states, particularly how cell states change their molecular profiles. This can include novel gene activation or the novel deployment of programs already existing in other cell types, known as co-option. RESULTS: Here we present scEvoNet, a Python tool for predicting cell type evolution in cross-species or cancer-related scRNA-seq datasets. ScEvoNet builds the confusion matrix of cell states and a bipartite network connecting genes and cell states. It allows a user to obtain a set of genes shared by the characteristic signature of two cell states even between distantly-related datasets. These genes can be used as indicators of either evolutionary divergence or co-option occurring during organism or tumor evolution. Our results on cancer and developmental datasets indicate that scEvoNet is a helpful tool for the initial screening of such genes as well as for measuring cell state similarities. CONCLUSION: The scEvoNet package is implemented in Python and is freely available from https://github.com/monsoro/scEvoNet . Utilizing this framework and exploring the continuum of transcriptome states between developmental stages and species will help explain cell state dynamics.

Taming human brain organoids one cell at a time.

Human brain organoids are self-organizing three-dimensional structures that emerge from human pluripotent stem cells and mimic aspects of the cellular composition and functionality of the developing human brain. Despite their impressive self-organizing capacity, organoids lack the stereotypic structural anatomy of their in vivo counterpart, making conventional analysis techniques underpowered to assess cellular composition and gene network regulation in organoids. Advances in single cell transcriptomics have recently allowed characterization and improvement of organoid protocols, as they continue to evolve, by enabling identification of cell types and states along with their developmental origins. In this review, we summarize recent approaches, progresses and challenges in resolving brain organoid's complexity through single-cell transcriptomics. We then discuss emerging technologies that may complement single-cell RNA sequencing by providing additional readouts of cellular states to generate an organ-level view of developmental processes. Altogether, these integrative technologies will allow monitoring of global gene regulation in thousands of individual cells and will offer an unprecedented opportunity to investigate features of human brain development and disease across multiple cellular modalities and with cell-type resolution.

Dissecting cellular states of infiltrating microenvironment cells in melanoma by integrating single-cell and bulk transcriptome analysis.

BACKGROUND: Cellular states of different immune cells can affect the activity of the whole immune microenvironment. METHODS: Here, leveraging reference profiles of microenvironment cell states that were constructed based on single-cell RNA-seq data of melanoma, we dissected the composition of microenvironment cell states across 463 skin cutaneous melanoma (SKCM) bulk samples through CIBERSORT-based deconvolution of gene expression profiles and revealed high heterogeneity of their distribution. Correspondence analysis on the estimated cellular fractions of melanoma bulk samples was performed to identify immune phenotypes. Based on the publicly available clinical survival and therapy data, we analyzed the relationship between immune phenotypes and clinical outcomes of melanoma. RESULTS: By analysis of the relationships among those cell states, we further identified three distinct tumor microenvironment immune phenotypes: "immune hot/active", "immune cold-suppressive" and "immune cold-exhausted". They were characterized by markedly different patterns of cell states: most notably the CD8 T Cytotoxic state, CD8 T Mixed state, B non-regulatory state and cancer-associated fibroblasts (CAFs), depicting distinct types of antitumor immune response (or immune activity). These phenotypes had prognostic significance for progression-free survival and implications in response to immune therapy in an independent cohort of anti-PD1 treated melanoma patients. CONCLUSIONS: The proposed strategy of leveraging single-cell data to dissect the composition of microenvironment cell states in individual bulk tumors can also extend to other cancer types, and our results highlight the importance of microenvironment cell states for the understanding of tumor immunity.

FunRes: resolving tissue-specific functional cell states based on a cell-cell communication network model.

The functional specialization of cell types arises during development and is shaped by cell-cell communication networks determining a distribution of functional cell states that are collectively important for tissue functioning. However, the identification of these tissue-specific functional cell states remains challenging. Although a plethora of computational approaches have been successful in detecting cell types and subtypes, they fail in resolving tissue-specific functional cell states. To address this issue, we present FunRes, a computational method designed for the identification of functional cell states. FunRes relies on scRNA-seq data of a tissue to initially reconstruct the functional cell-cell communication network, which is leveraged for partitioning each cell type into functional cell states. We applied FunRes to 177 cell types in 10 different tissues and demonstrated that the detected states correspond to known functional cell states of various cell types, which cannot be recapitulated by existing computational tools. Finally, we characterize emerging and vanishing functional cell states in aging and disease, and demonstrate their involvement in key tissue functions. Thus, we believe that FunRes will be of great utility in the characterization of the functional landscape of cell types and the identification of dysfunctional cell states in aging and disease.

MicroRNAs organize intrinsic variation into stem cell states.

Pluripotent embryonic stem cells (ESCs) contain the potential to form a diverse array of cells with distinct gene expression states, namely the cells of the adult vertebrate. Classically, diversity has been attributed to cells sensing their position with respect to external morphogen gradients. However, an alternative is that diversity arises in part from cooption of fluctuations in the gene regulatory network. Here we find ESCs exhibit intrinsic heterogeneity in the absence of external gradients by forming interconverting cell states. States vary in developmental gene expression programs and display distinct activity of microRNAs (miRNAs). Notably, miRNAs act on neighborhoods of pluripotency genes to increase variation of target genes and cell states. Loss of miRNAs that vary across states reduces target variation and delays state transitions, suggesting variable miRNAs organize and propagate variation to promote state transitions. Together these findings provide insight into how a gene regulatory network can coopt variation intrinsic to cell systems to form robust gene expression states. Interactions between intrinsic heterogeneity and environmental signals may help achieve developmental outcomes.

A periodic table of cell types.

Single cell biology is currently revolutionizing developmental and evolutionary biology, revealing new cell types and states in an impressive range of biological systems. With the accumulation of data, however, the field is grappling with a central unanswered question: what exactly is a cell type? This question is further complicated by the inherently dynamic nature of developmental processes. In this Hypothesis article, we propose that a 'periodic table of cell types' can be used as a framework for distinguishing cell types from cell states, in which the periods and groups correspond to developmental trajectories and stages along differentiation, respectively. The different states of the same cell type are further analogous to 'isotopes'. We also highlight how the concept of a periodic table of cell types could be useful for predicting new cell types and states, and for recognizing relationships between cell types throughout development and evolution.

Turning Data to Knowledge: Online Tools, Databases, and Resources in microRNA Research.

MicroRNAs (miRNAs) provide a fundamental layer of regulation in cells. miRNAs act posttranscriptionally through complementary base-pairing with the 3'-UTR of a target mRNA, leading to mRNA degradation and translation arrest. The likelihood of forming a valid miRNA-target duplex within cells was computationally predicted and experimentally monitored. In human cells, the miRNA profiles determine their identity and physiology. Therefore, alterations in the composition of miRNAs signify many cancer types and chronic diseases. In this chapter, we introduce online functional tools and resources to facilitate miRNA research. We start by introducing currently available miRNA catalogs and miRNA-gateway portals for navigating among different miRNA-centric online resources. We then sketch several realistic challenges that may occur while investigating miRNA regulation in living cells. As a showcase, we demonstrate the utility of miRNAs and mRNAs expression databases that cover diverse human cells and tissues, including resources that report on genetic alterations affecting miRNA expression levels and alteration in binding capacity. Introducing tools linking miRNAs with transcription factor (TF) networks reveals miRNA regulation complexity within living cells. Finally, we concentrate on online resources that analyze miRNAs in human diseases and specifically in cancer. Altogether, we introduce contemporary, selected resources and online tools for studying miRNA regulation in cells and tissues and their utility in health and disease.

Single-Cell Transcriptome of Wet AMD Patient-Derived Endothelial Cells in Angiogenic Sprouting.

Age-related macular degeneration (AMD) is a global leading cause of visual impairment in older populations. 'Wet' AMD, the most common subtype of this disease, occurs when pathological angiogenesis infiltrates the subretinal space (choroidal neovascularization), causing hemorrhage and retinal damage. Gold standard anti-vascular endothelial growth factor (VEGF) treatment is an effective therapy, but the long-term prevention of visual decline has not been as successful. This warrants the need to elucidate potential VEGF-independent pathways. We generated blood out-growth endothelial cells (BOECs) from wet AMD and normal control subjects, then induced angiogenic sprouting of BOECs using a fibrin gel bead assay. To deconvolute endothelial heterogeneity, we performed single-cell transcriptomic analysis on the sprouting BOECs, revealing a spectrum of cell states. Our wet AMD BOECs share common pathways with choroidal neovascularization such as extracellular matrix remodeling that promoted proangiogenic phenotype, and our 'activated' BOEC subpopulation demonstrated proinflammatory hallmarks, resembling the tip-like cells in vivo. We uncovered new molecular insights that pathological angiogenesis in wet AMD BOECs could also be driven by interleukin signaling and amino acid metabolism. A web-based visualization of the sprouting BOEC single-cell transcriptome has been created to facilitate further discovery research.

What we know about cardiomyocyte dedifferentiation.

Cardiomyocytes (CMs) lost during cardiac injury and heart failure (HF) cannot be replaced due to their limited proliferative capacity. Regenerating the failing heart by promoting CM cell-cycle re-entry is an ambitious solution, currently vigorously pursued. Some genes have been proven to promote endogenous CM proliferation, believed to be preceded by CM dedifferentiation, wherein terminally differentiated CMs are initially reversed back to the less mature state which precedes cell division. However, very little else is known about CM dedifferentiation which remains poorly defined. We lack robust molecular markers and proper understanding of the mechanisms driving dedifferentiation. Even the term dedifferentiation is debated because there is no objective evidence of pluripotency, and could rather reflect CM plasticity instead. Nonetheless, the significance of CM transition states on cardiac function, and whether they necessarily lead to CM proliferation, remains unclear. This review summarises the current state of knowledge of both natural and experimentally induced CM dedifferentiation in non-mammalian vertebrates (primarily the zebrafish) and mammals, as well as the phenotypes and molecular mechanisms involved. The significance and potential challenges of studying CM dedifferentiation are also discussed. In summary, CM dedifferentiation, essential for CM plasticity, may have an important role in heart regeneration, thereby contributing to the prevention and treatment of heart disease. More attention is needed in this field to overcome the technical limitations and knowledge gaps.

Distinct cellular states determine calcium signaling response.

The heterogeneity in mammalian cells signaling response is largely a result of pre-existing cell-to-cell variability. It is unknown whether cell-to-cell variability rises from biochemical stochastic fluctuations or distinct cellular states. Here, we utilize calcium response to adenosine trisphosphate as a model for investigating the structure of heterogeneity within a population of cells and analyze whether distinct cellular response states coexist. We use a functional definition of cellular state that is based on a mechanistic dynamical systems model of calcium signaling. Using Bayesian parameter inference, we obtain high confidence parameter value distributions for several hundred cells, each fitted individually. Clustering the inferred parameter distributions revealed three major distinct cellular states within the population. The existence of distinct cellular states raises the possibility that the observed variability in response is a result of structured heterogeneity between cells. The inferred parameter distribution predicts, and experiments confirm that variability in IP3R response explains the majority of calcium heterogeneity. Our work shows how mechanistic models and single-cell parameter fitting can uncover hidden population structure and demonstrate the need for parameter inference at the single-cell level.

Control of cell state transitions by post-transcriptional regulation.

Cell state transitions are prevalent in biology, playing a fundamental role in development, homeostasis and repair. Dysregulation of cell state transitions can lead to or occur in a wide range of diseases. In this letter, I explore and highlight the role of post-transcriptional regulatory mechanisms in determining the dynamics of cell state transitions. I propose that regulation of protein levels after transcription provides an under-appreciated regulatory route to obtain fast and sharp transitions between distinct cell states. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.