A multimodal zebrafish developmental atlas reveals the state-transition dynamics of late-vertebrate pluripotent axial progenitors.

Elucidating organismal developmental processes requires a comprehensive understanding of cellular lineages in the spatial, temporal, and molecular domains. In this study, we introduce Zebrahub, a dynamic atlas of zebrafish embryonic development that integrates single-cell sequencing time course data with lineage reconstructions facilitated by light-sheet microscopy. This atlas offers high-resolution and in-depth molecular insights into zebrafish development, achieved through the sequencing of individual embryos across ten developmental stages, complemented by reconstructions of cellular trajectories. Zebrahub also incorporates an interactive tool to navigate the complex cellular flows and lineages derived from light-sheet microscopy data, enabling in silico fate-mapping experiments. To demonstrate the versatility of our multimodal resource, we utilize Zebrahub to provide fresh insights into the pluripotency of neuro-mesodermal progenitors (NMPs) and the origins of a joint kidney-hemangioblast progenitor population.

Dissecting the treatment-naive ecosystem of human melanoma brain metastasis.

Melanoma brain metastasis (MBM) frequently occurs in patients with advanced melanoma; yet, our understanding of the underlying salient biology is rudimentary. Here, we performed single-cell/nucleus RNA-seq in 22 treatment-naive MBMs and 10 extracranial melanoma metastases (ECMs) and matched spatial single-cell transcriptomics and T cell receptor (TCR)-seq. Cancer cells from MBM were more chromosomally unstable, adopted a neuronal-like cell state, and enriched for spatially variably expressed metabolic pathways. Key observations were validated in independent patient cohorts, patient-derived MBM/ECM xenograft models, RNA/ATAC-seq, proteomics, and multiplexed imaging. Integrated spatial analyses revealed distinct geography of putative cancer immune evasion and evidence for more abundant intra-tumoral B to plasma cell differentiation in lymphoid aggregates in MBM. MBM harbored larger fractions of monocyte-derived macrophages and dysfunctional TOX+CD8+ T cells with distinct expression of immune checkpoints. This work provides comprehensive insights into MBM biology and serves as a foundational resource for further discovery and therapeutic exploration.

Progressive Pulmonary Fibrosis Is Caused by Elevated Mechanical Tension on Alveolar Stem Cells.

Fibrosis can develop in most organs and causes organ failure. The most common type of lung fibrosis is known as idiopathic pulmonary fibrosis, in which fibrosis starts at the lung periphery and then progresses toward the lung center, eventually causing respiratory failure. Little is known about the mechanisms underlying the pathogenesis and periphery-to-center progression of the disease. Here we discovered that loss of Cdc42 function in alveolar stem cells (AT2 cells) causes periphery-to-center progressive lung fibrosis. We further show that Cdc42-null AT2 cells in both post-pneumonectomy and untreated aged mice cannot regenerate new alveoli, resulting in sustained exposure of AT2 cells to elevated mechanical tension. We demonstrate that elevated mechanical tension activates a TGF-ß signaling loop in AT2 cells, which drives the periphery-to-center progression of lung fibrosis. Our study establishes a direct mechanistic link between impaired alveolar regeneration, mechanical tension, and progressive lung fibrosis.

Self-Reporting Transposons Enable Simultaneous Readout of Gene Expression and Transcription Factor Binding in Single Cells.

Cellular heterogeneity confounds in situ assays of transcription factor (TF) binding. Single-cell RNA sequencing (scRNA-seq) deconvolves cell types from gene expression, but no technology links cell identity to TF binding sites (TFBS) in those cell types. We present self-reporting transposons (SRTs) and use them in single-cell calling cards (scCC), a novel assay for simultaneously measuring gene expression and mapping TFBS in single cells. The genomic locations of SRTs are recovered from mRNA, and SRTs deposited by exogenous, TF-transposase fusions can be used to map TFBS. We then present scCC, which map SRTs from scRNA-seq libraries, simultaneously identifying cell types and TFBS in those same cells. We benchmark multiple TFs with this technique. Next, we use scCC to discover BRD4-mediated cell-state transitions in K562 cells. Finally, we map BRD4 binding sites in the mouse cortex at single-cell resolution, establishing a new method for studying TF biology in situ.

Tumor-reactive T cell clonotype dynamics underlying clinical response to TIL therapy in melanoma.

Adoptive cell therapy (ACT) using in vitro expanded tumor-infiltrating lymphocytes (TILs) has inconsistent clinical responses. To better understand determinants of therapeutic success, we tracked TIL clonotypes from baseline tumors to ACT products and post-ACT blood and tumor samples in melanoma patients using single-cell RNA and T cell receptor (TCR) sequencing. Patients with clinical responses had baseline tumors enriched in tumor-reactive TILs, and these were more effectively mobilized upon in vitro expansion, yielding products enriched in tumor-specific CD8+ cells that preferentially infiltrated tumors post-ACT. Conversely, lack of clinical responses was associated with tumors devoid of tumor-reactive resident clonotypes and with cell products mostly composed of blood-borne clonotypes that persisted in blood but not in tumors post-ACT. Upon expansion, tumor-specific TILs lost tumor-associated transcriptional signatures, including exhaustion, and responders exhibited an intermediate exhausted effector state after TIL engraftment in the tumor, suggesting functional reinvigoration. Our findings provide insight into the nature and dynamics of tumor-specific clonotypes associated with clinical response to TIL-ACT, with implications for treatment optimization.

Haploinsufficiency of phosphodiesterase 10A activates PI3K/AKT signaling independent of PTEN to induce an aggressive glioma phenotype.

Glioblastoma is universally fatal and characterized by frequent chromosomal copy number alterations harboring oncogenes and tumor suppressors. In this study, we analyzed exome-wide human glioblastoma copy number data and found that cytoband 6q27 is an independent poor prognostic marker in multiple data sets. We then combined CRISPR-Cas9 data, human spatial transcriptomic data, and human and mouse RNA sequencing data to nominate PDE10A as a potential haploinsufficient tumor suppressor in the 6q27 region. Mouse glioblastoma modeling using the RCAS/tv-a system confirmed that Pde10a suppression induced an aggressive glioma phenotype in vivo and resistance to temozolomide and radiation therapy in vitro. Cell culture analysis showed that decreased Pde10a expression led to increased PI3K/AKT signaling in a Pten-independent manner, a response blocked by selective PI3K inhibitors. Single-nucleus RNA sequencing from our mouse gliomas in vivo, in combination with cell culture validation, further showed that Pde10a suppression was associated with a proneural-to-mesenchymal transition that exhibited increased cell adhesion and decreased cell migration. Our results indicate that glioblastoma patients harboring PDE10A loss have worse outcomes and potentially increased sensitivity to PI3K inhibition.

Transcription factor antagonism regulates heterogeneity in embryonic stem cell states.

Gene expression heterogeneity underlies cell states and contributes to developmental robustness. While heterogeneity can arise from stochastic transcriptional processes, the extent to which it is regulated is unclear. Here, we characterize the regulatory program underlying heterogeneity in murine embryonic stem cell (mESC) states. We identify differentially active and transcribed enhancers (DATEs) across states. DATEs regulate differentially expressed genes and are distinguished by co-binding of transcription factors Klf4 and Zfp281. In contrast to other factors that interact in a positive feedback network stabilizing mESC cell-type identity, Klf4 and Zfp281 drive opposing transcriptional and chromatin programs. Abrogation of factor binding to DATEs dampens variation in gene expression, and factor loss alters kinetics of switching between states. These results show antagonism between factors at enhancers results in gene expression heterogeneity and formation of cell states, with implications for the generation of diverse cell types during development.

Discovering mechanisms of human genetic variation and controlling cell states at scale.

Population-scale sequencing efforts have catalogued substantial genetic variation in humans such that variant discovery dramatically outpaces interpretation. We discuss how single-cell sequencing is poised to reveal genetic mechanisms at a rate that may soon approach that of variant discovery. The functional genomics toolkit is sufficiently modular to systematically profile almost any type of variation within increasingly diverse contexts and with molecularly comprehensive and unbiased readouts. As a result, we can construct deep phenotypic atlases of variant effects that span the entire regulatory cascade. The same conceptual approach to interpreting genetic variation should be applied to engineering therapeutic cell states. In this way, variant mechanism discovery and cell state engineering will become reciprocating and iterative processes towards genomic medicine.

Hippo pathway deletion in adult resting cardiac fibroblasts initiates a cell state transition with spontaneous and self-sustaining fibrosis.

Cardiac fibroblasts (CFs) respond to injury by transitioning through multiple cell states, including resting CFs, activated CFs, and myofibroblasts. We report here that Hippo signaling cell-autonomously regulates CF fate transitions and proliferation, and non-cell-autonomously regulates both myeloid and CF activation in the heart. Conditional deletion of Hippo pathway kinases, Lats1 and Lats2, in uninjured CFs initiated a self-perpetuating fibrotic response in the adult heart that was exacerbated by myocardial infarction (MI). Single cell transcriptomics showed that uninjured Lats1/2 mutant CFs spontaneously transitioned to a myofibroblast cell state. Through gene regulatory network reconstruction, we found that Hippo-deficient myofibroblasts deployed a network of transcriptional regulators of endoplasmic reticulum (ER) stress, and the unfolded protein response (UPR) consistent with elevated secretory activity. We observed an expansion of myeloid cell heterogeneity in uninjured Lats1/2 CKO hearts with similarity to cells recovered from control hearts post-MI. Integrated genome-wide analysis of Yap chromatin occupancy revealed that Yap directly activates myofibroblast cell identity genes, the proto-oncogene Myc, and an array of genes encoding pro-inflammatory factors through enhancer-promoter looping. Our data indicate that Lats1/2 maintain the resting CF cell state through restricting the Yap-induced injury response.

Hippo signaling does it again: arbitrating cardiac fibroblast identity and activation.

The Hippo pathway is an evolutionarily conserved kinase cascade that is fundamental for tissue development, homeostasis, and regeneration. In the developing mammalian heart, Hippo signaling regulates cardiomyocyte numbers and organ size. While cardiomyocytes in the adult heart are largely postmitotic, Hippo deficiency can increase proliferation of these cells and affect cardiac regenerative capacity. Recent studies have also shown that resident cardiac fibroblasts play a critical role in disease responsiveness and healing, and in this issue of Genes and Development, Xiao and colleagues (pp. 1491-1505) demonstrate that Hippo signaling also integrates the activity of fibroblasts in the heart. They show that Hippo signaling normally maintains the cardiac fibroblast in a resting state and, conversely, its inactivation during disease-related stress results in a spontaneous transition toward a myofibroblast state that underlies fibrosis and ventricular remodeling. This phenotypic switch is associated with increased cytokine signaling that promotes nonautonomous resident fibroblast and myeloid cell activation.

Methods and Insights from Single-Cell Expression Quantitative Trait Loci.

Recent advancements in single-cell technologies have enabled expression quantitative trait locus (eQTL) analysis across many individuals at single-cell resolution. Compared with bulk RNA sequencing, which averages gene expression across cell types and cell states, single-cell assays capture the transcriptional states of individual cells, including fine-grained, transient, and difficult-to-isolate populations at unprecedented scale and resolution. Single-cell eQTL (sc-eQTL) mapping can identify context-dependent eQTLs that vary with cell states, including some that colocalize with disease variants identified in genome-wide association studies. By uncovering the precise contexts in which these eQTLs act, single-cell approaches can unveil previously hidden regulatory effects and pinpoint important cell states underlying molecular mechanisms of disease. Here, we present an overview of recently deployed experimental designs in sc-eQTL studies. In the process, we consider the influence of study design choices such as cohort, cell states, and ex vivo perturbations. We then discuss current methodologies, modeling approaches, and technical challenges as well as future opportunities and applications.

Dissecting the polygenic basis of atherosclerosis via disease-associated cell state signatures.

Coronary artery disease (CAD) is a pandemic disease where up to half of the risk is explained by genetic factors. Advanced insights into the genetic basis of CAD require deeper understanding of the contributions of different cell types, molecular pathways, and genes to disease heritability. Here, we investigate the biological diversity of atherosclerosis-associated cell states and interrogate their contribution to the genetic risk of CAD by using single-cell and bulk RNA sequencing (RNA-seq) of mouse and human lesions. We identified 12 disease-associated cell states that we characterized further by gene set functional profiling, ligand-receptor prediction, and transcription factor inference. Importantly, Vcam1+ smooth muscle cell state genes contributed most to SNP-based heritability of CAD. In line with this, genetic variants near smooth muscle cell state genes and regulatory elements explained the largest fraction of CAD-risk variance between individuals. Using this information for variant prioritization, we derived a hybrid polygenic risk score (PRS) that demonstrated improved performance over a classical PRS. Our results provide insights into the biological mechanisms associated with CAD risk, which could make a promising contribution to precision medicine and tailored therapeutic interventions in the future.

Clearing the slate: RNA turnover to enable cell state switching?

The distribution of mRNA in tissue is determined by the balance between transcription and decay. Understanding the control of RNA decay during development has been somewhat neglected compared with transcriptional control. Here, we explore the potential for mRNA decay to trigger rapid cell state transitions during development, comparing a bistable switch model of cell state conversion with experimental evidence from different developmental systems. We also consider another potential role for large-scale RNA decay that has emerged from studies of stress-induced cell state transitions, in which removal of mRNA unblocks the translation machinery to prioritise the synthesis of proteins that establish the new cell state.

Collective signalling drives rapid jumping between cell states.

Development can proceed in 'fits and starts', with rapid transitions between cell states involving concerted transcriptome-wide changes in gene expression. However, it is not clear how these transitions are regulated in complex cell populations, in which cells receive multiple inputs. We address this issue using Dictyostelium cells undergoing development in their physiological niche. A continuous single cell transcriptomics time series identifies a sharp 'jump' in global gene expression marking functionally different cell states. By simultaneously imaging the physiological dynamics of transcription and signalling, we show the jump coincides with the onset of collective oscillations of cAMP. Optogenetic control of cAMP pulses shows that different jump genes respond to distinct dynamic features of signalling. Late jump gene expression changes are almost completely dependent on cAMP, whereas transcript changes at the onset of the jump require additional input. The coupling of collective signalling with gene expression is a potentially powerful strategy to drive robust cell state transitions in heterogeneous signalling environments. Based on the context of the jump, we also conclude that sharp gene expression transitions may not be sufficient for commitment.

Topological and geometric analysis of cell states in single-cell transcriptomic data.

Single-cell RNA sequencing (scRNA-seq) enables dissecting cellular heterogeneity in tissues, resulting in numerous biological discoveries. Various computational methods have been devised to delineate cell types by clustering scRNA-seq data, where clusters are often annotated using prior knowledge of marker genes. In addition to identifying pure cell types, several methods have been developed to identify cells undergoing state transitions, which often rely on prior clustering results. The present computational approaches predominantly investigate the local and first-order structures of scRNA-seq data using graph representations, while scRNA-seq data frequently display complex high-dimensional structures. Here, we introduce scGeom, a tool that exploits the multiscale and multidimensional structures in scRNA-seq data by analyzing the geometry and topology through curvature and persistent homology of both cell and gene networks. We demonstrate the utility of these structural features to reflect biological properties and functions in several applications, where we show that curvatures and topological signatures of cell and gene networks can help indicate transition cells and the differentiation potential of cells. We also illustrate that structural characteristics can improve the classification of cell types.

scSTAR reveals hidden heterogeneity with a real-virtual cell pair structure across conditions in single-cell RNA sequencing data.

Cell-state transition can reveal additional information from single-cell ribonucleic acid (RNA)-sequencing data in time-resolved biological phenomena. However, most of the current methods are based on the time derivative of the gene expression state, which restricts them to the short-term evolution of cell states. Here, we present single-cell State Transition Across-samples of RNA-seq data (scSTAR), which overcomes this limitation by constructing a paired-cell projection between biological conditions with an arbitrary time span by maximizing the covariance between two feature spaces using partial least square and minimum squared error methods. In mouse ageing data, the response to stress in CD4+ memory T cell subtypes was found to be associated with ageing. A novel Treg subtype characterized by mTORC activation was identified to be associated with antitumour immune suppression, which was confirmed by immunofluorescence microscopy and survival analysis in 11 cancers from The Cancer Genome Atlas Program. On melanoma data, scSTAR improved immunotherapy-response prediction accuracy from 0.8 to 0.96.

Dynamical hallmarks of cancer: Phenotypic switching in melanoma and epithelial-mesenchymal plasticity.

Phenotypic plasticity was recently incorporated as a hallmark of cancer. This plasticity can manifest along many interconnected axes, such as stemness and differentiation, drug-sensitive and drug-resistant states, and between epithelial and mesenchymal cell-states. Despite growing acceptance for phenotypic plasticity as a hallmark of cancer, the dynamics of this process remains poorly understood. In particular, the knowledge necessary for a predictive understanding of how individual cancer cells and populations of cells dynamically switch their phenotypes in response to the intensity and/or duration of their current and past environmental stimuli remains far from complete. Here, we present recent investigations of phenotypic plasticity from a systems-level perspective using two exemplars: epithelial-mesenchymal plasticity in carcinomas and phenotypic switching in melanoma. We highlight how an integrated computational-experimental approach has helped unravel insights into specific dynamical hallmarks of phenotypic plasticity in different cancers to address the following questions: a) how many distinct cell-states or phenotypes exist?; b) how reversible are transitions among these cell-states, and what factors control the extent of reversibility?; and c) how might cell-cell communication be able to alter rates of cell-state switching and enable diverse patterns of phenotypic heterogeneity? Understanding these dynamic features of phenotypic plasticity may be a key component in shifting the paradigm of cancer treatment from reactionary to a more predictive, proactive approach.

Spatial transcriptomic brain imaging reveals the effects of immunomodulation therapy on specific regional brain cells in a mouse dementia model.

Increasing evidence of brain-immune crosstalk raises expectations for the efficacy of novel immunotherapies in Alzheimer's disease (AD), but the lack of methods to examine brain tissues makes it difficult to evaluate therapeutics. Here, we investigated the changes in spatial transcriptomic signatures and brain cell types using the 10x Genomics Visium platform in immune-modulated AD models after various treatments. To proceed with an analysis suitable for barcode-based spatial transcriptomics, we first organized a workflow for segmentation of neuroanatomical regions, establishment of appropriate gene combinations, and comprehensive review of altered brain cell signatures. Ultimately, we investigated spatial transcriptomic changes following administration of immunomodulators, NK cell supplements and an anti-CD4 antibody, which ameliorated behavior impairment, and designated brain cells and regions showing probable associations with behavior changes. We provided the customized analytic pipeline into an application named STquantool. Thus, we anticipate that our approach can help researchers interpret the real action of drug candidates by simultaneously investigating the dynamics of all transcripts for the development of novel AD therapeutics.

Activation of alveolar epithelial ER stress by ß-coronavirus infection disrupts surfactant homeostasis in mice: implications for COVID-19 respiratory failure.

COVID-19 syndrome is characterized by acute lung injury, hypoxemic respiratory failure, and high mortality. Alveolar type 2 (AT2) cells are essential for gas exchange, repair, and regeneration of distal lung epithelium. We have shown that the causative agent, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), and other members of the ß-coronavirus genus induce an endoplasmic reticulum (ER) stress response in vitro; however, the consequences for host AT2 cell function in vivo are less understood. To study this, two murine models of coronavirus infection were used-mouse hepatitis virus-1 (MHV-1) in A/J mice and a mouse-adapted SARS-CoV-2 strain. MHV-1-infected mice exhibited dose-dependent weight loss with histological evidence of distal lung injury accompanied by elevated bronchoalveolar lavage fluid (BALF) cell counts and total protein. AT2 cells showed evidence of both viral infection and increased BIP/GRP78 expression, consistent with activation of the unfolded protein response (UPR). The AT2 UPR included increased inositol-requiring enzyme 1α (IRE1α) signaling and a biphasic response in PKR-like ER kinase (PERK) signaling accompanied by marked reductions in AT2 and BALF surfactant protein (SP-B and SP-C) content, increases in surfactant surface tension, and emergence of a reprogrammed epithelial cell population (Krt8+ and Cldn4+). The loss of a homeostatic AT2 cell state was attenuated by treatment with the IRE1α inhibitor OPK-711. As a proof-of-concept, C57BL6 mice infected with mouse-adapted SARS-CoV-2 demonstrated similar lung injury and evidence of disrupted surfactant homeostasis. We conclude that lung injury from ß-coronavirus infection results from an aberrant host response, activating multiple AT2 UPR stress pathways, altering surfactant metabolism/function, and changing AT2 cell state, offering a mechanistic link between SARS-CoV-2 infection, AT2 cell biology, and acute respiratory failure.NEW & NOTEWORTHY COVID-19 syndrome is characterized by hypoxemic respiratory failure and high mortality. In this report, we use two murine models to show that ß-coronavirus infection produces acute lung injury, which results from an aberrant host response, activating multiple epithelial endoplasmic reticular stress pathways, disrupting pulmonary surfactant metabolism and function, and forcing emergence of an aberrant epithelial transition state. Our results offer a mechanistic link between SARS-CoV-2 infection, AT2 cell biology, and respiratory failure.

Single cell spatial biology over developmental time can decipher pediatric brain pathologies.

Pediatric low grade brain tumors and neurodevelopmental disorders share proteins, signaling pathways, and networks. They also share germline mutations and an impaired prenatal differentiation origin. They may differ in the timing of the events and proliferation. We suggest that their pivotal distinct, albeit partially overlapping, outcomes relate to the cell states, which depend on their spatial location, and timing of gene expression during brain development. These attributes are crucial as the brain develops sequentially, and single-cell spatial organization influences cell state, thus function. Our underlying premise is that the root cause in neurodevelopmental disorders and pediatric tumors is impaired prenatal differentiation. Data related to pediatric brain tumors, neurodevelopmental disorders, brain cell (sub)types, locations, and timing of expression in the developing brain are scant. However, emerging single cell technologies, including transcriptomic, spatial biology, spatial high-resolution imaging performed over the brain developmental time, could be transformational in deciphering brain pathologies thereby pharmacology.