Search

Abstract The Human Cell Atlas (HCA) is a global partnership “to create comprehensive reference maps of all human cells—the fundamental units of life – as a basis for both understanding human health and diagnosing, monitoring, and treating disease.” ( https://www.humancellatlas.org/ ) The atlas shall characterize cells from diverse individuals across the globe to better understand human biology. HCA proactively considers the priorities of, and benefits accrued to, contributing communities. Here, we lay out principles and action items that have been adopted to affirm HCA’s commitment to equity so that the atlas is beneficial to all of humanity.

Graphical Abstract In this Editorial, our Chief Editor and members of our Advisory Editorial Board discuss recent breakthroughs, current challenges, and emerging opportunities in single‐cell biology and share their vision of “where the field is headed.”

Mono-allelic germline disruptions of the transcription factor GATA2 result in a propensity for developing myelodysplastic syndrome (MDS) and acute myeloid leukemia (AML), affecting more than 85% of carriers. How a partial loss of GATA2 functionality enables leukemic transformation years later is unclear. This question has remained unsolved mainly due to the lack of informative models, as Gata2 heterozygote mice do not develop hematologic malignancies. Here we show that two different germline Gata2 mutations (TgErg/Gata2het and TgErg/Gata2L359V) accelerate AML in mice expressing the human hematopoietic stem cell regulator ERG. Analysis of Erg/Gata2het fetal liver and bone marrow-derived hematopoietic cells revealed a distinct pre-leukemic phenotype. This was characterized by enhanced transition from stem to progenitor state, increased proliferation, and a striking mitochondrial phenotype, consisting of highly expressed oxidative-phosphorylation-related gene sets, elevated oxygen consumption rates, and notably, markedly distorted mitochondrial morphology. Importantly, the same mitochondrial gene-expression signature was observed in human AML harboring GATA2 aberrations. Similar to the observations in mice, non-leukemic bone marrows from children with germline GATA2 mutation demonstrated marked mitochondrial abnormalities. Thus, we observed the tumor suppressive effects of GATA2 in two germline Gata2 genetic mouse models. As oncogenic mutations often accumulate with age, GATA2 deficiency-mediated priming of hematopoietic cells for oncogenic transformation may explain the earlier occurrence of MDS/AML in patients with GATA2 germline mutation. The mitochondrial phenotype is a potential therapeutic opportunity for the prevention of leukemic transformation in these patients.

Objectives: Until recently, communication between neighboring cells in islets of Langerhans was overlooked by genomic technologies, which require rigorous tissue dissociation into single cells. Methods: We utilize sorting of physically interacting cells (PICs) with single-cell RNA-sequencing to systematically map cellular interactions in the endocrine pancreas after pancreatectomy. Results: The pancreas cellular landscape features pancreatectomy associated heterogeneity of beta-cells, including an interaction-specific program between paired beta and delta-cells. Conclusions: Our analysis suggests that the particular cluster of beta-cells that pairs with delta-cells benefits from stress protection, implying that the interaction between beta- and delta-cells might safeguard against pancreatectomy associated challenges. The work encourages testing the potential relevance of physically-interacting beta-delta-cells also in diabetes mellitus.

From single-cell transcriptome analyses to defining culture media for spheroids, the authors provide a census of information to understand the development of human pancreatic progenitors. This approach identifies signalling pathways (EGF and FGF) regulating progenitor proliferation.

The control of off-target activity is a challenge for adapting CRISPR to therapeutic use. Here the authors present CRISPECTOR, a software tool to detect, evaluate and quantify editing activity, including translocations, from NGS data.

Combining scRNA-seq with spatial information to enable the reconstruction of spatially-resolved cell atlases is challenging for rare cell types. Here the authors present ClumpSeq, an approach for sequencing small clumps of tissue attached cells, and apply it to establish spatial atlases for all secretory cell types in the small intestine.

Integrated studies of matched tissue sites and cell types in COVID-19 patients are important to define the immune mechanisms of pathology. Here, the authors describe an immune signature in fatal COVID-19 patients harmonizing single-cell RNA sequencing of blood and matched BAL cells with deep clinical, immunological and functional data.

Abstract NF-κB is a well-characterized transcription factor, widely known for its roles in inflammation and immune responses, as well as in control of cell division and apoptosis. However, its function in β-cells is still being debated, as it appears to depend on the timing and kinetics of its activation. To elucidate the temporal role of NF-κB in vivo, we have generated two transgenic mouse models, the ToIβ and NOD/ToIβ mice, in which NF-κB activation is specifically and conditionally inhibited in β-cells. In this study, we present a novel function of the canonical NF-κB pathway during murine islet β-cell development. Interestingly, inhibiting the NF-κB pathway in β-cells during embryogenesis, but not after birth, in both ToIβ and NOD/ToIβ mice, increased β-cell turnover, ultimately resulting in a reduced β-cell mass. On the NOD background, this was associated with a marked increase in insulitis and diabetes incidence. While a robust nuclear immunoreactivity of the NF-κB p65-subunit was found in neonatal β-cells, significant activation was not detected in β-cells of either adult NOD/ToIβ mice or in the pancreata of recently diagnosed adult T1D patients. Moreover, in NOD/ToIβ mice, inhibiting NF-κB post-weaning had no effect on the development of diabetes or β-cell dysfunction. In conclusion, our data point to NF-κB as an important component of the physiological regulatory circuit that controls the balance of β-cell proliferation and apoptosis in the early developmental stages of insulin-producing cells, thus modulating β-cell mass and the development of diabetes in the mouse model of T1D.

Abstract Malignant cell growth is fueled by interactions between tumor cells and the stromal cells composing the tumor microenvironment. The human liver is a major site of tumors and metastases, but molecular identities and intercellular interactions of different cell types have not been resolved in these pathologies. Here, we apply single cell RNA‐sequencing and spatial analysis of malignant and adjacent non‐malignant liver tissues from five patients with cholangiocarcinoma or liver metastases. We find that stromal cells exhibit recurring, patient‐independent expression programs, and reconstruct a ligand–receptor map that highlights recurring tumor–stroma interactions. By combining transcriptomics of laser‐capture microdissected regions, we reconstruct a zonation atlas of hepatocytes in the non‐malignant sites and characterize the spatial distribution of each cell type across the tumor microenvironment. Our analysis provides a resource for understanding human liver malignancies and may expose potential points of interventions.

Genome editing of human cluster of differentiation 34+ (CD34+) hematopoietic stem and progenitor cells (HSPCs) holds great therapeutic potential. This study aimed to optimize on-target, ex vivo genome editing using the CRISPR-Cas9 system in CD34+ HSPCs and to create a clear workflow for precise identification of off-target effects. Modified synthetic guide RNAs (gRNAs), either 2-part gRNA or single-guide RNA (sgRNA), were delivered to CD34+ HSPCs as part of ribonucleoprotein (RNP) complexes, targeting therapeutically relevant genes. The addition of an Alt-R electroporation enhancer (EE), a short, single-stranded oligodeoxynucleotide (ssODN), significantly increased editing efficiency in CD34+ HSPCs. Notably, similar editing improvement was observed when excess gRNA over Cas9 protein was used, providing a DNA-free alternative suitable for therapeutic applications. Furthermore, we demonstrated that sgRNA may be preferable over 2-part gRNA in a locus-specific manner. Finally, we present a clear experimental framework suitable for the unbiased identification of bona fide off-target sites by Genome-Wide, Unbiased Identification of Double-Strand Breaks (DSBs) Enabled by Sequencing (GUIDE-seq), as well as subsequent editing quantification in CD34+ HSPCs using rhAmpSeq. These findings may facilitate the implementation of genome editing in CD34+ HSPCs for research and therapy and can be adapted for other hematopoietic cells.

Epithelial gene expression has been shown to be zonated along the crypt-villus axis, but mechanisms shaping this spatial variability were unknown. Here, Bahar Halpern et al. uncover zonation of mesenchymal cells, including Lgr5+ telocytes, which regulate epithelial gene expression at the villus tip.

The lineage-specific receptor CSF1R controls macrophage development and homeostasis. Here the authors show that deletion of a conserved Csf1r enhancer (FIRE) selectively depletes brain microglia and resident macrophages in the epidermis, kidney, heart and peritoneum of otherwise healthy mice.

Blocking the PD-1 pathway was shown to be effective in amyloid beta mouse models, yet little is known about its therapeutic potential in models of tauopathy. The authors show here that blocking PD-L1, a PD-1 ligand, is similarly effective, and that both treatments reversed cognitive deficiencies, and modified disease pathology not only in an animal model of AD, but also in the DM-hTAU mouse tauopathy model, through a mechanism that involves monocyte-derived macrophages.

Senescent cells accumulate with aging contributing to age-related disease and the role of the immune system in removing senescent cells is not completely understood. Here, the authors show that perforin deficient mice accumulate more senescent cells and have a shorter lifespan, and that this phenotype can be reversed with administration of a senolytic drug.

The mechanical challenge of attaching elastic tendons to stiff bones is solved by the formation of a unique transitional tissue. Here, we show that murine tendon-to-bone attachment cells are bi-fated, activating a mixture of chondrocyte and tenocyte transcriptomes, under regulation of shared regulatory elements and Krüppel-like factors (KLFs) transcription factors. High-throughput bulk and single-cell RNA sequencing of humeral attachment cells revealed expression of hundreds of chondrogenic and tenogenic genes, which was validated by in situ hybridization and single-molecule ISH. ATAC sequencing showed that attachment cells share accessible intergenic chromatin areas with either tenocytes or chondrocytes. Epigenomic analysis revealed enhancer signatures for most of these regions. Transgenic mouse enhancer reporter assays verified the shared activity of some of these enhancers. Finally, integrative chromatin and motif analyses and transcriptomic data implicated KLFs as regulators of attachment cells. Indeed, blocking expression of both Klf2 and Klf4 in developing limb mesenchyme impaired their differentiation.

Summary: Hematopoietic stem cells (HSCs) are considered a heterogeneous cell population. To further resolve the HSC compartment, we characterized a retinoic acid (RA) reporter mouse line. Sub-fractionation of the HSC compartment in RA-CFP reporter mice demonstrated that RA-CFP-dim HSCs were largely non-proliferative and displayed superior engraftment potential in comparison with RA-CFP-bright HSCs. Gene expression analysis demonstrated higher expression of RA-target genes in RA-CFP-dim HSCs, in contrast to the RA-CFP reporter expression, but both RA-CFP-dim and RA-CFP-bright HSCs responded efficiently to RA in vitro. Single-cell RNA sequencing (RNA-seq) of >1,200 HSCs showed that differences in cell cycle activity constituted the main driver of transcriptional heterogeneity in HSCs. Moreover, further analysis of the single-cell RNA-seq data revealed that stochastic low-level expression of distinct lineage-affiliated transcriptional programs is a common feature of HSCs. Collectively, this work demonstrates the utility of the RA-CFP reporter line as a tool for the isolation of superior HSCs. : HSCs are considered a functional heterogeneous population. Lauridsen et al. use scRNA-seq to demonstrate that most transcriptional heterogeneity within the HSC compartment is associated with differences in cell cycle status. They further use an RA-CFP reporter mouse line to isolate slow-cycling HSCs characterized by superior engraftment potential. Keywords: hematopoietic stem cells, single-cell RNA-sequencing, retinoic acid, transcriptional heterogeneity, hematopoiesis

Microglia cells in the brain regulate immune responses, but in ageing can negatively affect brain function. Here the authors show that the chronic presence of type I interferon in aged mouse brain impedes cognitive ability by altering microglia transcriptome and limiting Mef2C, a microglia ‘off’ signal.

The extent and nature of epigenomic changes associated with melanoma progression is poorly understood. Through systematic epigenomic profiling of 35 epigenetic modifications and transcriptomic analysis, we define chromatin state changes associated with melanomagenesis by using a cell phenotypic model of non-tumorigenic and tumorigenic states. Computation of specific chromatin state transitions showed loss of histone acetylations and H3K4me2/3 on regulatory regions proximal to specific cancer-regulatory genes in important melanoma-driving cell signaling pathways. Importantly, such acetylation changes were also observed between benign nevi and malignant melanoma human tissues. Intriguingly, only a small fraction of chromatin state transitions correlated with expected changes in gene expression patterns. Restoration of acetylation levels on deacetylated loci by histone deacetylase (HDAC) inhibitors selectively blocked excessive proliferation in tumorigenic cells and human melanoma cells, suggesting functional roles of observed chromatin state transitions in driving hyperproliferative phenotype. Through these results, we define functionally relevant chromatin states associated with melanoma progression.

Summary: The glucocorticoid receptor (GR) acts as a ubiquitous cortisol-dependent transcription factor (TF). To identify co-factors, we used protein-fragment complementation assays and found that GR recognizes FLI1 and additional ETS family proteins, TFs relaying proliferation and/or migration signals. Following steroid-dependent translocation of FLI1 and GR to the nucleus, the FLI1-specific domain (FLS) binds with GR and strongly enhances GR’s transcriptional activity. This interaction has functional consequences in Ewing sarcoma (ES), childhood and adolescence bone malignancies driven by fusions between EWSR1 and FLI1. In vitro, GR knockdown inhibited the migration and proliferation of ES cells, and in animal models, antagonizing GR (or lowering cortisol) retarded both tumor growth and metastasis from bone to lung. Taken together, our findings offer mechanistic rationale for repurposing GR-targeting drugs for the treatment of patients with ES. : The single oncogene of Ewing sarcoma (ES), a childhood cancer, encodes a FLI1 fusion protein. Srivastava et al. report physical interactions between the fusion protein and the glucocorticoid receptor. Drug-induced inhibition of these interactions retards the progression of ES in mouse models. Keywords: cancer therapy, Ewing sarcoma, glucocorticoid receptor, metastasis, protein-fragment complementation assay

Dendritic cells (DC) are unrivaled in their potential to prime naive T cells by presenting antigen and providing costimulation. DC are furthermore believed to decode antigen context by virtue of pattern recognition receptors and to polarize T cells through cytokine secretion toward distinct effector functions. Diverse polarized T helper (TH) cells have been explored in great detail. In contrast, studies of instructing DC have to date largely been restricted to in vitro settings or adoptively transferred DC. Here we report efforts to unravel the DC response to cognate T cell encounter in antigen-challenged lymph nodes (LN). Mice engrafted with antigen-specific T cells were immunized with nanoparticles (NP) entrapping adjuvants and absorbed with antigen to study the immediate DC response to T cell encounter using bulk and single cell RNA-seq profiling. NP induced robust antigen-specific TH1 cell responses with minimal bystander activation. Fluorescent-labeled NP allowed identification of antigen-carrying DC and focus on transcriptional changes in DC that encounter T cells. Our results support the existence of a bi-directional crosstalk between DC and T cells that promotes TH1 responses, including involvement of the ubiquitin-like molecule Isg15 that merits further study.

It is well established that inducible transcription is essential for the consolidation of salient experiences into long-term memory. However, whether inducible transcription relays information about the identity and affective attributes of the experience being encoded, has not been explored. To this end, we analyzed transcription induced by a variety of rewarding and aversive experiences, across multiple brain regions. Our results describe the existence of robust transcriptional signatures uniquely representing distinct experiences, enabling near-perfect decoding of recent experiences. Furthermore, experiences with shared attributes display commonalities in their transcriptional signatures, exemplified in the representation of valence, habituation and reinforcement. This study introduces the concept of a neural transcriptional code, which represents the encoding of experiences in the mouse brain. This code is comprised of distinct transcriptional signatures that correlate to attributes of the experiences that are being committed to long-term memory.

The recent advent of methods for high-throughput single-cell molecular profiling has catalyzed a growing sense in the scientific community that the time is ripe to complete the 150-year-old effort to identify all cell types in the human body. The Human Cell Atlas Project is an international collaborative effort that aims to define all human cell types in terms of distinctive molecular profiles (such as gene expression profiles) and to connect this information with classical cellular descriptions (such as location and morphology). An open comprehensive reference map of the molecular state of cells in healthy human tissues would propel the systematic study of physiological states, developmental trajectories, regulatory circuitry and interactions of cells, and also provide a framework for understanding cellular dysregulation in human disease. Here we describe the idea, its potential utility, early proofs-of-concept, and some design considerations for the Human Cell Atlas, including a commitment to open data, code, and community.

Depletion of essential nutrients triggers regulatory programs that prolong cell growth and survival. Starvation-induced processes increase nutrient transport, mobilize nutrient storage, and recycle nutrients between cellular components. This leads to an effective increase in intracellular nutrients, which may act as a negative feedback that downregulates the starvation program. To examine how cells overcome this potential instability, we followed the transcription response of budding yeast transferred to medium lacking phosphate. Genes were induced in two temporal waves. The first wave was stably maintained and persisted even upon phosphate replenishment, indicating a positive feedback loop. This commitment was abolished after 2 hr with the induction of the second expression wave, coinciding with the reduction in cell growth rate. We show that the overall temporal stability of the expression response depends on the sequential pattern of gene induction. Our results emphasize the key role of gene expression dynamics in optimizing cellular adaptation.

Abstract Hundreds of immune cell types work in coordination to maintain tissue homeostasis. Upon infection, dramatic changes occur with the localization, migration, and proliferation of the immune cells to first alert the body of the danger, confine it to limit spreading, and finally extinguish the threat and bring the tissue back to homeostasis. Since current technologies can follow the dynamics of only a limited number of cell types, we have yet to grasp the full complexity of global in vivo cell dynamics in normal developmental processes and disease. Here, we devise a computational method, digital cell quantification (DCQ), which combines genome‐wide gene expression data with an immune cell compendium to infer in vivo changes in the quantities of 213 immune cell subpopulations. DCQ was applied to study global immune cell dynamics in mice lungs at ten time points during 7 days of flu infection. We find dramatic changes in quantities of 70 immune cell types, including various innate, adaptive, and progenitor immune cells. We focus on the previously unreported dynamics of four immune dendritic cell subtypes and suggest a specific role for CD103+ CD11b− DCs in early stages of disease and CD8+ pDC in late stages of flu infection.

Abstract Recognition of microbial danger signals by toll‐like receptors (TLR) causes re‐programming of macrophages. To investigate kinase cascades triggered by the TLR4 ligand lipopolysaccharide (LPS) on systems level, we performed a global, quantitative and kinetic analysis of the phosphoproteome of primary macrophages using stable isotope labelling with amino acids in cell culture, phosphopeptide enrichment and high‐resolution mass spectrometry. In parallel, nascent RNA was profiled to link transcription factor (TF) phosphorylation to TLR4‐induced transcriptional activation. We reproducibly identified 1850 phosphoproteins with 6956 phosphorylation sites, two thirds of which were not reported earlier. LPS caused major dynamic changes in the phosphoproteome (24% up‐regulation and 9% down‐regulation). Functional bioinformatic analyses confirmed canonical players of the TLR pathway and highlighted other signalling modules (e.g. mTOR, ATM/ATR kinases) and the cytoskeleton as hotspots of LPS‐regulated phosphorylation. Finally, weaving together phosphoproteome and nascent transcriptome data by in silico promoter analysis, we implicated several phosphorylated TFs in primary LPS‐controlled gene expression.

Abstract Robust biological signaling networks evolved, through gene duplications, from simple, relatively fragile cascades. Architectural features such as layered configuration, branching and modularity, as well as functional characteristics (e.g., feedback control circuits), enable fail‐safe performance in the face of internal and external perturbations. These universal features are exemplified here using the receptor tyrosine kinase (RTK) family. The RTK module is richly mutated and overexpressed in human malignancies, and pharmaceutical interception of its signaling effectively retards growth of specific tumors. Therapy‐induced interception of RTK‐signaling pathways and the common evolvement of drug resistance are respectively considered here as manifestations of fragility and plasticity of robust networks. The systems perspective we present views pathologies as hijackers of biological robustness and offers ways for identifying fragile hubs, as well as strategies to overcome drug resistance.