Search

Single-gene missense mutations remain challenging to interpret. Here, we deploy scalable functional screening by sequencing (SEUSS), a Perturb-seq method, to generate mutations at protein interfaces of RUNX1 and quantify their effect on activities of downstream cellular programs. We evaluate single-cell RNA profiles of 115 mutations in myelogenous leukemia cells and categorize them into three functionally distinct groups, wild-type (WT)-like, loss-of-function (LoF)-like, and hypomorphic, that we validate in orthogonal assays. LoF-like variants dominate the DNA-binding site and are recurrent in cancer; however, recurrence alone does not predict functional impact. Hypomorphic variants share characteristics with LoF-like but favor protein interactions, promoting gene expression indicative of nerve growth factor (NGF) response and cytokine recruitment of neutrophils. Accessible DNA near differentially expressed genes frequently contains RUNX1-binding motifs. Finally, we reclassify 16 variants of uncertain significance and train a classifier to predict 103 more. Our work demonstrates the potential of targeting protein interactions to better define the landscape of phenotypes reachable by missense mutations.

Human bone marrow (BM) plasma cells are heterogeneous, ranging from newly arrived antibody-secreting cells (ASCs) to long-lived plasma cells (LLPCs). We provide single-cell transcriptional resolution of 17,347 BM ASCs from five healthy adults. Fifteen clusters are identified ranging from newly minted ASCs (cluster 1) expressing MKI67 and high major histocompatibility complex (MHC) class II that progress to late clusters 5-8 through intermediate clusters 2-4. Additional ASC clusters include the following: immunoglobulin (Ig) M predominant (likely of extra-follicular origin), interferon responsive, and high mitochondrial activity. Late ASCs are distinguished by G2M checkpoints, mammalian target of rapamycin (mTOR) signaling, distinct metabolic pathways, CD38 expression, utilization of tumor necrosis factor (TNF)-receptor superfamily members, and two distinct maturation pathways involving TNF signaling through nuclear factor κB (NF-κB). This study provides a single-cell atlas and molecular roadmap of LLPC maturation trajectories essential in the BM microniche. Altogether, understanding BM ASC heterogeneity in health and disease enables development of new strategies to enhance protective ASCs and to deplete pathogenic ones.

The enteric nervous system (ENS) consists of glial cells (EGCs) and neurons derived from neural crest precursors. EGCs retain capacity for large-scale neurogenesis in culture, and in vivo lineage tracing has identified neurons derived from glial cells in response to inflammation. We thus hypothesize that EGCs possess a chromatin structure poised for neurogenesis. We use single-cell multiome sequencing to simultaneously assess transcription and chromatin accessibility in EGCs undergoing spontaneous neurogenesis in culture, as well as small intestine myenteric plexus EGCs. Cultured EGCs maintain open chromatin at genomic loci accessible in neurons, and neurogenesis from EGCs involves dynamic chromatin rearrangements with a net decrease in accessible chromatin. A subset of in vivo EGCs, highly enriched within the myenteric ganglia and that persist into adulthood, have a gene expression program and chromatin state consistent with neurogenic potential. These results clarify the mechanisms underlying EGC potential for neuronal fate transition.

Silencing of transposable elements (TEs) drives the evolution of numerous redundant mechanisms of transcriptional regulation. Arabidopsis MBD5, MBD6, and SILENZIO act as TE repressors downstream of DNA methylation. Here, we show, via single-nucleus RNA-seq of developing male gametophytes, that these repressors are critical for TE silencing in the pollen vegetative cell, a companion cell important for fertilization that undergoes chromatin decompaction. Instead, other silencing mutants (met1, ddm1, mom1, morc) show loss of silencing in all pollen nucleus types and somatic cells. We show that TEs repressed by MBD5/6 gain chromatin accessibility in wild-type vegetative nuclei despite remaining silent, suggesting that loss of DNA compaction makes them sensitive to loss of MBD5/6. Consistently, crossing mbd5/6 to histone 1 mutants, which have decondensed chromatin in leaves, reveals derepression of MBD5/6-dependent TEs in leaves. MBD5/6 and SILENZIO thus act as a silencing system particularly important when chromatin compaction is compromised.

Behavior arises from concerted activity throughout the brain. Consequently, a major focus of modern neuroscience is defining the physiology and behavioral roles of projection neurons linking different brain areas. Single-cell RNA sequencing has facilitated these efforts by revealing molecular determinants of cellular physiology and markers that enable genetically targeted perturbations such as optogenetics, but existing methods for sequencing defined projection populations are low throughput, painstaking, and costly. We developed a straightforward, multiplexed approach, virally encoded connectivity transgenic overlay RNA sequencing (VECTORseq). VECTORseq repurposes commercial retrogradely infecting viruses typically used to express functional transgenes (e.g., recombinases and fluorescent proteins) by treating viral transgene mRNA as barcodes within single-cell datasets. VECTORseq is compatible with different viral families, resolves multiple populations with different projection targets in one sequencing run, and identifies cortical and subcortical excitatory and inhibitory projection populations. Our study provides a roadmap for high-throughput identification of neuronal subtypes based on connectivity.

Neural crest (NC) cells migrate throughout vertebrate embryos to give rise to a huge variety of cell types, but when and where lineages emerge and their regulation remain unclear. We have performed single-cell RNA sequencing (RNA-seq) of cranial NC cells from the first pharyngeal arch in zebrafish over several stages during migration. Computational analysis combining pseudotime and real-time data reveals that these NC cells first adopt a transitional state, becoming specified mid-migration, with the first lineage decisions being skeletal and pigment, followed by neural and glial progenitors. In addition, by computationally integrating these data with RNA-seq data from a transgenic Wnt reporter line, we identify gene cohorts with similar temporal responses to Wnts during migration and show that one, Atp6ap2, is required for melanocyte differentiation. Together, our results show that cranial NC cell lineages arise progressively and uncover a series of spatially restricted cell interactions likely to regulate such cell-fate decisions.

The mechanisms of Myc-driven liver tumorigenesis are inadequately understood. Herein we show that Myc-driven hepatocellular carcinoma (HCC) is dramatically aggravated in mice with hepatocyte-specific Ptpn11/Shp2 deletion. However, Myc-induced tumors develop selectively from the rare Shp2-positive hepatocytes in Shp2-deficent liver, and Myc-driven oncogenesis depends on an intact Ras-Erk signaling promoted by Shp2 to sustain Myc stability. Despite a stringent requirement of Shp2 cell autonomously, Shp2 deletion induces an immunosuppressive environment, resulting in defective clearance of tumor-initiating cells and aggressive tumor progression. The basal Wnt/β-catenin signaling is upregulated in Shp2-deficient liver, which is further augmented by Myc transfection. Ablating Ctnnb1 suppresses Myc-induced HCC in Shp2-deficient livers, revealing an essential role of β-catenin. Consistently, Myc overexpression and CTNNB1 mutations are frequently co-detected in HCC patients with poor prognosis. These data elucidate complex mechanisms of liver tumorigenesis driven by cell-intrinsic oncogenic signaling in cooperation with a tumor-promoting microenvironment generated by disrupting the specific oncogenic pathway.

The human fallopian tube harbors the cell of origin for the majority of high-grade serous "ovarian" cancers (HGSCs), but its cellular composition, particularly the epithelial component, is poorly characterized. We perform single-cell transcriptomic profiling of around 53,000 individual cells from 12 primary fallopian specimens to map their major cell types. We identify 10 epithelial subpopulations with diverse transcriptional programs. Based on transcriptional signatures, we reconstruct a trajectory whereby secretory cells differentiate into ciliated cells via a RUNX3high intermediate. Computational deconvolution of advanced HGSCs identifies the "early secretory" population as a likely precursor state for the majority of HGSCs. Its signature comprises both epithelial and mesenchymal features and is enriched in mesenchymal-type HGSCs (p = 6.7 × 10-27), a group known to have particularly poor prognoses. This cellular and molecular compendium of the human fallopian tube in cancer-free women is expected to advance our understanding of the earliest stages of fallopian epithelial neoplasia.

The mammary epithelial cell (MEC) system is a bilayered ductal epithelium of luminal and basal cells, maintained by a lineage of stem and progenitor populations. Here, we used integrated single-cell transcriptomics and chromatin accessibility analysis to reconstruct the cell types of the mouse MEC system and their underlying gene regulatory features in an unbiased manner. We define differentiation states within the secretory type of luminal cells, which forms a continuous spectrum of general luminal progenitor and lactation-committed progenitor cells. By integrating single-cell transcriptomics and chromatin accessibility landscapes, we identify cis- and trans-regulatory elements that are differentially activated in the specific epithelial cell types and our newly defined luminal differentiation states. Our work provides a resource to reveal cis/trans-regulatory elements associated with MEC identity and differentiation that will serve as a reference to determine how the chromatin accessibility landscape changes during breast cancer.

The type 2 cytokine-high asthma endotype (T2H) is characterized by IL-13-driven mucus obstruction of the airways. To further investigate this incompletely understood pathobiology, we characterize IL-13 effects on human airway epithelial cell cultures using single-cell RNA sequencing, finding that IL-13 generates a distinctive transcriptional state for each cell type. Specifically, we discover a mucus secretory program induced by IL-13 in all cell types which converts both mucus and defense secretory cells into a metaplastic state with emergent mucin production and secretion, while leading to ER stress and cell death in ciliated cells. The IL-13-remodeled epithelium secretes a pathologic, mucin-imbalanced, and innate immunity-depleted proteome that arrests mucociliary motion. Signatures of IL-13-induced cellular remodeling are mirrored by transcriptional signatures characteristic of the nasal airway epithelium within T2H versus T2-low asthmatic children. Our results reveal the epithelium-wide scope of T2H asthma and present candidate therapeutic targets for restoring normal epithelial function.

A certain number of epithelial cells in intestinal crypts are DNA damage resistant and contribute to regeneration. However, the cellular mechanism underlying intestinal regeneration remains unclear. Using lineage tracing, we show that cells marked by an Msi1 reporter (Msi1+) are right above Lgr5high cells in intestinal crypts and exhibit DNA damage resistance. Single-cell RNA sequencing reveals that the Msi1+ cells are heterogeneous with the majority being intestinal stem cells (ISCs). The DNA damage-resistant subpopulation of Msi1+ cells is characterized by low-to-negative Lgr5 expression and is more rapidly cycling than Lgr5high radiosensitive crypt base columnar stem cells (CBCs). This enables an efficient repopulation of the intestinal epithelium at early stage when Lgr5high cells are not emerging. Furthermore, relative to CBCs, Msi1+ cells preferentially produce Paneth cells during homeostasis and upon radiation repair. Together, we demonstrate that the DNA damage-resistant Msi1+ cells are cycling ISCs that maintain and regenerate the intestinal epithelium.

Our knowledge of transcriptional heterogeneities in epithelial stem and progenitor cell compartments is limited. Epidermal basal cells sustain cutaneous tissue maintenance and drive wound healing. Previous studies have probed basal cell heterogeneity in stem and progenitor potential, but a comprehensive dissection of basal cell dynamics during differentiation is lacking. Using single-cell RNA sequencing coupled with RNAScope and fluorescence lifetime imaging, we identify three non-proliferative and one proliferative basal cell state in homeostatic skin that differ in metabolic preference and become spatially partitioned during wound re-epithelialization. Pseudotemporal trajectory and RNA velocity analyses predict a quasi-linear differentiation hierarchy where basal cells progress from Col17a1Hi/Trp63Hi state to early-response state, proliferate at the juncture of these two states, or become growth arrested before differentiating into spinous cells. Wound healing induces plasticity manifested by dynamic basal-spinous interconversions at multiple basal transcriptional states. Our study provides a systematic view of epidermal cellular dynamics, supporting a revised "hierarchical-lineage" model of homeostasis.

Fallopian tube secretory epithelial cells (FTSECs) are likely the main precursor cell type of high-grade serous ovarian cancers (HGSOCs), but these tumors may also arise from ovarian surface epithelial cells (OSECs). We profiled global landscapes of gene expression and active chromatin to characterize molecular similarities between OSECs (n = 114), FTSECs (n = 74), and HGSOCs (n = 394). A one-class machine learning algorithm predicts that most HGSOCs derive from FTSECs, with particularly high FTSEC scores in mesenchymal-type HGSOCs (padj < 8 × 10-4). However, a subset of HGSOCs likely derive from OSECs, particularly HGSOCs of the proliferative type (padj < 2 × 10-4), suggesting a dualistic model for HGSOC origins. Super-enhancer (SE) landscapes were also more similar between FTSECs and HGSOCs than between OSECs and HGSOCs (p < 2.2 × 10-16). The SOX18 transcription factor (TF) coincided with a HGSOC-specific SE, and ectopic overexpression of SOX18 in FTSECs caused epithelial-to-mesenchymal transition, indicating that SOX18 plays a role in establishing the mesenchymal signature of fallopian-derived HGSOCs.

Single-cell transcriptome profiling of heterogeneous tissues can provide high-resolution windows into developmental dynamics and environmental responses, but its application to plants has been limited. Here, we used the high-throughput Drop-seq approach to profile >12,000 cells from Arabidopsis roots. This identified numerous distinct cell types, covering all major root tissues and developmental stages, and illuminated specific marker genes for these populations. In addition, we demonstrate the utility of this approach to study the impact of environmental conditions on developmental processes. Analysis of roots grown with or without sucrose supplementation uncovers changes in the relative frequencies of cell types in response to sucrose. Finally, we characterize the transcriptome changes that occur across endodermis development and identify nearly 800 genes with dynamic expression as this tissue differentiates. Collectively, we demonstrate that single-cell RNA-seq can be used to profile developmental processes in plants and show how they can be altered by external stimuli.

Spermatogenesis has been intensely studied in rodents but remains poorly understood in humans. Here, we used single-cell RNA sequencing to analyze human testes. Clustering analysis of neonatal testes reveals several cell subsets, including cell populations with characteristics of primordial germ cells (PGCs) and spermatogonial stem cells (SSCs). In adult testes, we identify four undifferentiated spermatogonia (SPG) clusters, each of which expresses specific marker genes. We identify protein markers for the most primitive SPG state, allowing us to purify this likely SSC-enriched cell subset. We map the timeline of male germ cell development from PGCs through fetal germ cells to differentiating adult SPG stages. We also define somatic cell subsets in both neonatal and adult testes and trace their developmental trajectories. Our data provide a blueprint of the developing human male germline and supporting somatic cells. The PGC-like and SSC markers are candidates to be used for SSC therapy to treat infertility.

Aging is associated with tissue-level changes in cellular composition that are correlated with increased susceptibility to disease. Aging human mammary tissue shows skewed progenitor cell potency, resulting in diminished tumor-suppressive cell types and the accumulation of defective epithelial progenitors. Quantitative characterization of these age-emergent human cell subpopulations is lacking, impeding our understanding of the relationship between age and cancer susceptibility. We conducted single-cell resolution proteomic phenotyping of healthy breast epithelia from 57 women, aged 16-91 years, using mass cytometry. Remarkable heterogeneity was quantified within the two mammary epithelial lineages. Population partitioning identified a subset of aberrant basal-like luminal cells that accumulate with age and originate from age-altered progenitors. Quantification of age-emergent phenotypes enabled robust classification of breast tissues by age in healthy women. This high-resolution mapping highlighted specific epithelial subpopulations that change with age in a manner consistent with increased susceptibility to breast cancer.

Aging is associated with tissue-level changes in cellular composition that are correlated with increased susceptibility to disease. Aging human mammary tissue shows skewed progenitor cell potency, resulting in diminished tumor-suppressive cell types and the accumulation of defective epithelial progenitors. Quantitative characterization of these age-emergent human cell subpopulations is lacking, impeding our understanding of the relationship between age and cancer susceptibility. We conducted single-cell resolution proteomic phenotyping of healthy breast epithelia from 57 women, aged 16-91 years, using mass cytometry. Remarkable heterogeneity was quantified within the two mammary epithelial lineages. Population partitioning identified a subset of aberrant basal-like luminal cells that accumulate with age and originate from age-altered progenitors. Quantification of age-emergent phenotypes enabled robust classification of breast tissues by age in healthy women. This high-resolution mapping highlighted specific epithelial subpopulations that change with age in a manner consistent with increased susceptibility to breast cancer.

The developmental origins of most adult stem cells are poorly understood. Here, we report the identification of a transcription factor-RHOX10-critical for the initial establishment of spermatogonial stem cells (SSCs). Conditional loss of the entire 33-gene X-linked homeobox gene cluster that includes Rhox10 causes progressive spermatogenic decline, a phenotype indistinguishable from that caused by loss of only Rhox10. We demonstrate that this phenotype results from dramatically reduced SSC generation. By using a battery of approaches, including single-cell-RNA sequencing (scRNA-seq) analysis, we show that Rhox10 drives SSC generation by promoting pro-spermatogonia differentiation. Rhox10 also regulates batteries of migration genes and promotes the migration of pro-spermatogonia into the SSC niche. The identification of an X-linked homeobox gene that drives the initial generation of SSCs has implications for the evolution of X-linked gene clusters and sheds light on regulatory mechanisms influencing adult stem cell generation in general.

Age-related alterations in the immune system are starting to emerge as key contributors to impairments found in aged organs. A decline in regenerative capacity is a hallmark of tissue aging; however, the contribution of immune aging to regenerative failure is just starting to be explored. Here, we apply a strategy combining single-cell RNA sequencing with flow cytometry, histological analysis, and functional assays to perform a complete analysis of the immune environment of the aged regenerating skeletal muscle on a time course following injury with single-cell resolution. Our results reveal an unanticipated complexity and functional heterogeneity in immune populations within the skeletal muscle that have been regarded as homogeneous. Furthermore, we uncover a profound remodeling of both myeloid and lymphoid compartments in aging. These discoveries challenge established notions on immune regulation of skeletal muscle regeneration, providing a set of potential targets to improve skeletal muscle health and regenerative capacity in aging., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

The composition and cellular heterogeneity of the corpus cavernosum (CC) microenvironment have been characterized, but the spatial heterogeneity at the molecular level remains unexplored. In this study, we integrate single-cell RNA sequencing (scRNA-seq) and spatial transcriptome sequencing to comprehensively chart the spatial cellular landscape of the human and rat CC under normal and disease conditions. We observe differences in the proportions of cell subtypes and marker genes between humans and rats. Based on the analysis of the fibroblast (FB) niche, we also find that the enrichment scores of mechanical force signaling vary across different regions and correlate with the spatial distribution of FB subtypes. In vitro, the soft and hard extracellular matrix (ECM) induces the differentiation of FBs into apolipoprotein (APO)+ FBs and cartilage oligomeric matrix protein (COMP)+ FBs, respectively. In summary, our study provides a cross-species and physiopathological transcriptomic atlas of the CC, contributing to a further understanding of the molecular anatomy and regulation of penile erection., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

The "innate-like" T cell compartment, known as T inn , represents a diverse group of T cells that straddle the boundary between innate and adaptive immunity. We explore the transcriptional landscape of T inn compared to conventional T cells (T conv ) in the human thymus and blood using single-cell RNA sequencing (scRNA-seq) and flow cytometry. In human blood, the majority of T inn cells share an effector program driven by specific transcription factors, distinct from those governing T conv cells. Conversely, only a fraction of thymic T inn cells displays an effector phenotype, while others share transcriptional features with developing T conv cells, indicating potential divergent developmental pathways. Unlike the mouse, human T inn cells do not differentiate into multiple effector subsets but develop a mixed type 1/type 17 effector potential. Cross-species analysis uncovers species-specific distinctions, including the absence of type 2 T inn cells in humans, which implies distinct immune regulatory mechanisms across species., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Through the use of bulk measurements in metabolic organs, the circadian clock was shown to play roles in organismal energy homeostasis. However, the relationship between metabolic and circadian oscillations has not been studied in vivo at a single-cell level. Also, it is unknown whether the circadian clock controls metabolism in stem cells. We used a sensitive, noninvasive method to detect metabolic oscillations and circadian phase within epidermal stem cells in live mice at the single-cell level. We observe a higher NADH/NAD+ ratio, reflecting an increased glycolysis/oxidative phosphorylation ratio during the night compared to the day. Furthermore, we demonstrate that single-cell metabolic heterogeneity within the basal cell layer correlates with the circadian clock and that diurnal fluctuations in NADH/NAD+ ratio are Bmal1 dependent. Our data show that, in proliferating stem cells, the circadian clock coordinates activities of oxidative phosphorylation and glycolysis with DNA synthesis, perhaps as a protective mechanism against genotoxicity.

Natural silks are renewable proteins with impressive mechanical properties and biocompatibility that are useful in various fields. However, the cellular and spatial organization of silk-secreting organs remains unclear. Here, we combined single-nucleus and spatially resolved transcriptomics to systematically map the cellular and spatial composition of the silk glands (SGs) of mulberry silkworms late in larval development. This approach allowed us to profile SG cell types and cell state dynamics and identify regulatory networks and cell-cell communication related to efficient silk protein synthesis; key markers were validated via transgenic approaches. Notably, we demonstrated the indispensable role of the ecdysone receptor (ultraspiracle) in regulating endoreplication in SG cells. Our atlas presents the results of spatiotemporal analysis of silk-secreting organ architecture late in larval development; this atlas provides a valuable reference for elucidating the mechanism of efficient silk protein synthesis and developing sustainable products made from natural silk. [Display omitted] • Cell-type distribution and spatial organization in silkworm silk glands were mapped • Spatiotemporal gene expression patterns in the silk glands of late larvae were profiled • The characteristics and roles of types of cells in silk protein synthesis were analyzed • Resources were applied to explore the cell biology of silk-secreting organs Ma et al. present a spatiotemporal transcriptomic map of silkworm silk-secreting organs. This atlas provides insights into cellular heterogeneity, dynamic changes in the cell state, regulatory networks, and endoreplication regulation in silk glands, as well as a framework for high-resolution exploration of the mechanism of efficient silk protein synthesis. [ABSTRACT FROM AUTHOR]

The transition from two-dimensional (2D) to 3D growth likely facilitated plants to colonize land, but its heterogeneity is not well understood. In this study, we utilized single-cell RNA sequencing to analyze the moss Physcomitrium patens, whose morphogenesis involves a transition from 2D to 3D growth. We profiled over 17,000 single cells covering all major vegetative tissues, including 2D filaments (chloronema and caulonema) and 3D structures (bud and gametophore). Pseudotime analyses revealed larger numbers of candidate genes that determine cell fates for 2D tip elongation or 3D bud differentiation. Using weighted gene co-expression network analysis, we identified a module that connects β-type carbonic anhydrases (βCAs) with auxin. We further validated the cellular expression patterns of βCAs and demonstrated their roles in 3D gametophore development. Overall, our study provides insights into cellular heterogeneity in a moss and identifies molecular signatures that underpin the 2D-to-3D growth transition at single-cell resolution., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024. Published by Elsevier Inc.)

Heterogeneous resistance to immunotherapy remains a major challenge in cancer treatment, often leading to disease progression and death. Using CITE-seq and matched 40-plex PhenoCycler tissue imaging, we performed longitudinal multimodal single-cell analysis of tumors from metastatic melanoma patients with innate resistance, acquired resistance, or response to immunotherapy. We established the multimodal integration toolkit to align transcriptomic features, cellular epitopes, and spatial information to provide deeper insights into the tumors. With longitudinal analysis, we identified an "immune-striving" tumor microenvironment marked by peri-tumor lymphoid aggregates and low infiltration of T cells in the tumor and the emergence of MITF + SPARCL1 + and CENPF + melanoma subclones after therapy. The enrichment of B cell-associated signatures in the molecular composition of lymphoid aggregates was associated with better survival. These findings provide further insights into the establishment of microenvironmental cell interactions and molecular composition of spatial structures that could inform therapeutic intervention., Competing Interests: Declaration of interests G.V.L. is a consultant advisor for Agenus, Amgen, Array Biopharma, AstraZeneca, Boehringer Ingelheim, Bristol Myers Squibb, Evaxion, Hexal AG (Sandoz), Highlight Therapeutics S.L., Innovent Biologics USA, Merck Sharpe & Dohme, Novartis, OncoSec, PHMR Ltd., Pierre Fabre, Provectus, Qbiotics, and Regeneron. R.A.S. has received fees for professional services from SkylineDx BV, IO Biotech ApS, MetaOptima Technology, Roche, Evaxion, Provectus Biopharmaceuticals Australia, Qbiotics, Novartis, Merck Sharp & Dohme, NeraCare, AMGEN., Bristol Myers Squibb, Myriad Genetics, and GlaxoSmithKline. I.P.d.S. is on the advisory board of Merck Sharp & Dohme; has received fees for professional services from Roche, Bristol Myers Squibb, Merck Sharp & Dohme, Novartis, and Pierre Fabre; and has had travel support from Bristol Myers Squibb and Merck Sharp & Dohme. R.P.M.S. has received honoraria for advisory board participation from Merck Sharp & Dohme, Novartis, and Qbiotics, and speaking honoraria from Bristol Myers Squibb and Novartis. A.S. has received honoraria for advisory board participation from Nanostring and Phenomic. AI and research support was provided by Nanostring and 10X Genomics., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

A key issue for research on COVID-19 pathogenesis is the lack of biopsies from patients and of samples at the onset of infection. To overcome these hurdles, hamsters were shown to be useful models for studying this disease. Here, we further leverage the model to molecularly survey the disease progression from time-resolved single-cell RNA sequencing data collected from healthy and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-infected Syrian and Roborovski hamster lungs. We compare our data to human COVID-19 studies, including bronchoalveolar lavage, nasal swab, and postmortem lung tissue, and identify a shared axis of inflammation dominated by macrophages, neutrophils, and endothelial cells, which we show to be transient in Syrian and terminal in Roborovski hamsters. Our data suggest that, following SARS-CoV-2 infection, commitment to a type 1- or type 3-biased immunity determines moderate versus severe COVID-19 outcomes, respectively., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Quiescent adult neural stem cells (NSCs) in the mammalian brain arise from proliferating NSCs during development. Beyond acquisition of quiescence, an adult NSC hallmark, little is known about the process, milestones, and mechanisms underlying the transition of developmental NSCs to an adult NSC state. Here, we performed targeted single-cell RNA-seq analysis to reveal the molecular cascade underlying NSC development in the early postnatal mouse dentate gyrus. We identified two sequential steps, first a transition to quiescence followed by further maturation, each of which involved distinct changes in metabolic gene expression. Direct metabolic analysis uncovered distinct milestones, including an autophagy burst before NSC quiescence acquisition and cellular reactive oxygen species level elevation along NSC maturation. Functionally, autophagy is important for the NSC transition to quiescence during early postnatal development. Together, our study reveals a multi-step process with defined milestones underlying establishment of the adult NSC pool in the mammalian brain., Competing Interests: Declaration of interests The authors declare no other competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Mammary epithelial stem cells are vital to tissue expansion and remodeling during various phases of postnatal mammary development. Basal mammary epithelial cells are enriched in Wnt-responsive cells and can reconstitute cleared mammary fat pads upon transplantation into mice. Lgr5 is a Wnt-regulated target gene and was identified as a major stem cell marker in the small intestine, colon, stomach, and hair follicle, as well as in kidney nephrons. Here, we demonstrate the outstanding regenerative potential of a rare population of Lgr5-expressing (Lgr5(+)) mammary epithelial cells (MECs). We found that Lgr5(+) cells reside within the basal population, are superior to other basal cells in regenerating functional mammary glands (MGs), are exceptionally efficient in reconstituting MGs from single cells, and exhibit regenerative capacity in serial transplantations. Loss-of-function and depletion experiments of Lgr5(+) cells from transplanted MECs or from pubertal MGs revealed that these cells are not only sufficient but also necessary for postnatal mammary organogenesis.

Summary: High-dimensional data, such as those generated by single-cell RNA sequencing (scRNA-seq), present challenges in interpretation and visualization. Numerical and computational methods for dimensionality reduction allow for low-dimensional representation of genome-scale expression data for downstream clustering, trajectory reconstruction, and biological interpretation. However, a comprehensive and quantitative evaluation of the performance of these techniques has not been established. We present an unbiased framework that defines metrics of global and local structure preservation in dimensionality reduction transformations. Using discrete and continuous real-world and synthetic scRNA-seq datasets, we show how input cell distribution and method parameters are largely determinant of global, local, and organizational data structure preservation by 11 common dimensionality reduction methods.

Elucidating the complex relationships between mRNA and protein expression at high spatiotemporal resolution is critical for unraveling multilevel gene regulation and enhancing mRNA-based developmental analyses. In this study, we conduct a single-cell analysis of mRNA and protein expression of transcription factors throughout C. elegans embryogenesis. Initially, cellular co-presence of mRNA and protein is low, increasing to a medium-high level (73%) upon factoring in delayed protein synthesis and long-term protein persistence. These factors substantially affect mRNA-protein concordance, leading to potential inaccuracies in mRNA-reliant gene detection and specificity characterization. Building on the learned relationship, we infer protein presence from mRNA expression and demonstrate its utility in identifying tissue-specific genes and elucidating relationships between genes and cells. This approach facilitates identifying the role of sptf-1/SP7 in neuronal lineage development. Collectively, this study provides insights into gene expression dynamics during rapid embryogenesis and approaches for improving the efficacy of transcriptome-based developmental analyses. [Display omitted] • Analyzing single-cell mRNA and protein concordance during C. elegans embryogenesis • Delayed protein synthesis and prolonged persistence diminish mRNA-protein concordance • Identifying limitations of mRNA expression in developmental gene characterization • Inferring protein expression enhances mRNA-based developmental analyses Fan et al. conduct a cell-by-cell comparison of mRNA and protein expression throughout embryogenesis, characterizing that delayed protein synthesis and long-term protein persistence compromise developmental analyses that rely solely on mRNA. They further show that inferring protein expression from mRNA data enhances developmental analyses at both gene and cell levels. [ABSTRACT FROM AUTHOR]

Glioblastoma (GBM) is a deadly brain tumor, and the kinesin motor KIF11 is an attractive therapeutic target with roles in proliferation and invasion. Resistance to KIF11 inhibitors, which has mainly been studied in animal models, presents significant challenges. We use lineage-tracing barcodes and single-cell RNA sequencing to analyze resistance in patient-derived GBM neurospheres treated with ispinesib, a potent KIF11 inhibitor. Similar to GBM progression in patients, untreated cells lose their neural lineage identity and become mesenchymal, which is associated with poor prognosis. Conversely, cells subjected to long-term ispinesib treatment exhibit a proneural phenotype. We generate patient-derived xenografts and show that ispinesib-resistant cells form less aggressive tumors in vivo, even in the absence of drug. Moreover, treatment of human ex vivo GBM slices with ispinesib demonstrates phenotypic alignment with in vitro responses, underscoring the clinical relevance of our findings. Finally, using retrospective lineage tracing, we identify drugs that are synergistic with ispinesib., Competing Interests: Declaration of interests P.A.S. is listed as an inventor on patent applications and issued patents filed by Columbia University related to the microwell technology described here. P.A.S. receives patent royalties from Guardant Health., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

During cell fate transitions, cells remodel their transcriptome, chromatin, and epigenome; however, it has been difficult to determine the temporal dynamics and cause-effect relationship between these changes at the single-cell level. Here, we employ the heterokaryon-mediated reprogramming system as a single-cell model to dissect key temporal events during early stages of pluripotency conversion using super-resolution imaging. We reveal that, following heterokaryon formation, the somatic nucleus undergoes global chromatin decompaction and removal of repressive histone modifications H3K9me3 and H3K27me3 without acquisition of active modifications H3K4me3 and H3K9ac. The pluripotency gene OCT4 (POU5F1) shows nascent and mature RNA transcription within the first 24 h after cell fusion without requiring an initial open chromatin configuration at its locus. NANOG, conversely, has significant nascent RNA transcription only at 48 h after cell fusion but, strikingly, exhibits genomic reopening early on. These findings suggest that the temporal relationship between chromatin compaction and gene activation during cellular reprogramming is gene context dependent., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

The preoptic area of the hypothalamus (POA) is essential for sleep regulation. However, the cellular makeup of the POA is heterogeneous, and the molecular identities of the sleep-promoting cells remain elusive. To address this question, this study compares mice during recovery sleep following sleep deprivation to mice allowed extended sleep. Single-nucleus RNA sequencing (single-nucleus RNA-seq) identifies one galanin inhibitory neuronal subtype that shows upregulation of rapid and delayed activity-regulated genes during recovery sleep. This cell type expresses higher levels of growth hormone receptor and lower levels of estrogen receptor compared to other galanin subtypes. single-nucleus RNA-seq also reveals cell-type-specific upregulation of purinergic receptor (P2ry14) and serotonin receptor (Htr2a) during recovery sleep in this neuronal subtype, suggesting possible mechanisms for sleep regulation. Studies with RNAscope validate the single-nucleus RNA-seq findings. Thus, the combined use of single-nucleus RNA-seq and activity-regulated genes identifies a neuronal subtype functionally involved in sleep regulation., Competing Interests: Declaration of interests B.C.R. receives funding from Novo Nordisk and Boehringer Ingelheim that was not utilized for this project., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Ovarian endometriosis is characterized by the growth of endometrial tissue within the ovary, causing infertility and chronic pain. However, its pathophysiology remains unclear. Utilizing high-precision single-cell RNA sequencing, we profile the normal, eutopic, and ectopic endometrium from 34 individuals across proliferative and secretory phases. We observe an increased proportion of ciliated cells in both eutopic and ectopic endometrium, characterized by a diminished expression of estrogen sulfotransferase, which likely confers apoptosis resistance. After translocating to ectopic lesions, endometrial epithelium upregulates nicotinamide N-methyltransferase expression that inhibits apoptosis by promoting deacetylation and subsequent nuclear exclusion of transcription factor forkhead box protein O1, thereby leading to the downregulation of the apoptotic gene BIM. Moreover, epithelial cells in ectopic lesions elevate HLA class II complex expression, which stimulates CD4 + T cells and consequently contributes to chronic inflammation. Altogether, our study provides a comprehensive atlas of ovarian endometriosis and highlights potential therapeutic targets for modulating apoptosis and inflammation., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Cells receive time-varying signals from the environment and generate functional responses by secreting their own signaling molecules. Characterizing dynamic input-output relationships in single cells is crucial for understanding and modeling cellular systems. We developed an automated microfluidic system that delivers precisely defined dynamical inputs to individual living cells and simultaneously measures key immune parameters dynamically. Our system combines nanoliter immunoassays, microfluidic input generation, and time-lapse microscopy, enabling study of previously untestable aspects of immunity by measuring time-dependent cytokine secretion and transcription factor activity from single cells stimulated with dynamic inflammatory inputs. Employing this system to analyze macrophage signal processing under pathogen inputs, we found that the dynamics of TNF secretion are highly heterogeneous and surprisingly uncorrelated with the dynamics of NF-κB, the transcription factor controlling TNF production. Computational modeling of the LPS/TLR4 pathway shows that post-transcriptional regulation by TRIF is a key determinant of noisy and uncorrelated TNF secretion dynamics in single macrophages.

Summary We have developed a high-throughput, microfluidics-based platform to perform kinetic analysis of viral infections in individual cells. We have analyzed thousands of individual poliovirus infections while varying experimental parameters, including multiplicity of infection, cell cycle, viral genotype, and presence of a drug. We make several unexpected observations masked by population-based experiments: (1) viral and cellular factors contribute uniquely and independently to viral infection kinetics; (2) cellular factors cause wide variation in replication start times; and (3) infections frequently begin later and replication occurs faster than predicted by population measurements. We show that mutational load impairs interaction of the viral population with the host, delaying replication start times and explaining the attenuated phenotype of a mutator virus. We show that an antiviral drug can selectively extinguish the most-fit members of the viral population. Single-cell virology facilitates discovery and characterization of virulence determinants and elucidation of mechanisms of drug action eluded by population methods. [ABSTRACT FROM AUTHOR]

The stomach is an important digestive organ with various biological functions. However, because of the complexity of its cellular and glandular composition, its precise cellular biology has yet to be elucidated. In this study, we conducted single-cell RNA sequencing (scRNA-seq) and subcellular-level spatial transcriptomics analysis of the human stomach and constructed the largest dataset to date: a stomach encyclopedia. This dataset consists of approximately 380,000 cells from scRNA-seq and the spatial transcriptome, enabling integrated analyses of transcriptional and spatial information of gastric and metaplastic cells. This analysis identified LEFTY1 as an uncharacterized stem cell marker, which was confirmed through lineage tracing analysis. A wide variety of cell-cell interactions between epithelial and stromal cells, including PDGFRA + BMP4 + WNT5A + fibroblasts, was highlighted in the developmental switch of intestinal metaplasia. Our extensive dataset will function as a fundamental resource in investigations of the stomach, including studies of development, aging, and carcinogenesis., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 The Author(s). Published by Elsevier Inc. All rights reserved.)

Summary Cells receive time-varying signals from the environment and generate functional responses by secreting their own signaling molecules. Characterizing dynamic input-output relationships in single cells is crucial for understanding and modeling cellular systems. We developed an automated microfluidic system that delivers precisely defined dynamical inputs to individual living cells and simultaneously measures key immune parameters dynamically. Our system combines nanoliter immunoassays, microfluidic input generation, and time-lapse microscopy, enabling study of previously untestable aspects of immunity by measuring time-dependent cytokine secretion and transcription factor activity from single cells stimulated with dynamic inflammatory inputs. Employing this system to analyze macrophage signal processing under pathogen inputs, we found that the dynamics of TNF secretion are highly heterogeneous and surprisingly uncorrelated with the dynamics of NF-κB, the transcription factor controlling TNF production. Computational modeling of the LPS/TLR4 pathway shows that post-transcriptional regulation by TRIF is a key determinant of noisy and uncorrelated TNF secretion dynamics in single macrophages. [ABSTRACT FROM AUTHOR]