Search

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.”

Background & Aims: Noninvasive modalities for assessing active endoscopic and histologic inflammation in Crohn’s disease and ulcerative colitis patients are critically needed. Fecal wash host shed-cell transcriptomics has been shown to be a robust classifier of endoscopic and histologic inflammation in inflammatory bowel disease patients with distal colitis. Whether such fecal washes can inform on inflammatory processes occurring in more proximal intestinal segments is currently unknown. Methods: Fifty-nine inflammatory bowel disease patients and 50 controls were prospectively enrolled. Biopsy specimens and fecal washes from the distal colon, proximal colon, and terminal ileum were compared. Host transcriptomics were performed on the biopsy specimens and fecal washes obtained during colonoscopy at predefined locations throughout the colon and terminal ileum and results were associated with concurrent clinical, endoscopic, and histologic parameters. Results: We found that host transcriptomics of distal fecal washes robustly classify histologic inflammation in ileal and proximal colonic Crohn’s disease, even without distal colonic involvement (area under the receiver operating characteristic curve, 0.94 ± 0.09). We further found that fecal washes consist of modules of co-expressed genes of immune, stromal, and epithelial origin that are indicative of endoscopic disease severity. Fecal wash host transcriptomics also captures expression of gene modules previously associated with a lack of response to biological therapies. Conclusions: Our study establishes the accuracy of distal colonic fecal washes for identifying and scoring inflammatory processes throughout the entire ileal–colonic axis.

Necrotizing enterocolitis (NEC) is a gastrointestinal complication of premature infants with high rates of morbidity and mortality. A comprehensive view of the cellular changes and aberrant interactions that underlie NEC is lacking. This study aimed at filling in this gap. We combine single-cell RNA sequencing (scRNAseq), T-cell receptor beta (TCRβ) analysis, bulk transcriptomics, and imaging to characterize cell identities, interactions, and zonal changes in NEC. We find an abundance of proinflammatory macrophages, fibroblasts, endothelial cells as well as T cells that exhibit increased TCRβ clonal expansion. Villus tip epithelial cells are reduced in NEC and the remaining epithelial cells up-regulate proinflammatory genes. We establish a detailed map of aberrant epithelial-mesenchymal-immune interactions that are associated with inflammation in NEC mucosa. Our analyses highlight the cellular dysregulations of NEC-associated intestinal tissue and identify potential targets for biomarker discovery and therapeutics.

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.

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.

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.

Lactate is a by-product of glycolysis that can function via its G protein receptor GPR81. Here the authors show that LPS or Salmonella infection enhances glycolytic metabolism in bone marrow neutrophils, resulting in lactate production, which increases endothelial barrier permeability and mobilization of these neutrophils by targeting endothelial GPR81.

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 intestine is lined with isolated lymphoid follicles (ILFs) that facilitate sampling of luminal antigens to elicit immune responses. Technical challenges related to the scarcity and small sizes of ILFs and their follicle-associated epithelium (FAE) impeded the characterization of their spatial gene expression programs. Here, we combined RNA sequencing of laser capture microdissected tissues with single-molecule transcript imaging to obtain a spatial gene expression map of the ILF and its associated FAE in the mouse small intestine. We identified zonated expression programs in both follicles and FAEs, with a decrease in enterocyte antimicrobial and absorption programs and a partial induction of expression programs normally observed at the villus tip. We further identified Lepr+ subepithelial telocytes at the FAE top, which are distinct from villus tip Lgr5+ telocytes. Our analysis exposes the epithelial and mesenchymal cell states associated with ILFs.

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: We describe an optimized smFISH protocol for the intact pancreas. The protocol is adapted from Lyubimova et al. (2013), a generic tissue smFISH protocol that works for most tissues but not the pancreas. The main changes implemented include increasing the period of mRNA denaturation from 5 min to at least 3 h and increasing formamide concentrations from 10% to 30%. These modifications yield sensitive single mRNA visualization that is comparable to those achieved in other tissues using the standard protocol.For complete details on the use and execution of this protocol, please refer to Farack et al. (2018, 2019).

Abstract Mammalian glucose homeostasis is controlled by the antagonistic hormones insulin and glucagon, secreted by pancreatic beta and alpha cells respectively. These two cell types are adjacently located in the islets of Langerhans and affect each others’ secretions in a paradoxical manner: while insulin inhibits glucagon secretion from alpha cells, glucagon seems to stimulate insulin secretion from beta cells. Here we ask what are the design principles of this negative feedback loop. We systematically simulate the dynamics of all possible islet inter-cellular connectivity patterns and analyze different performance criteria. We find that the observed circuit dampens overshoots of blood glucose levels after reversion of glucose drops. This feature is related to the temporal delay in the rise of insulin concentrations in peripheral tissues, compared to the immediate hormone action on the liver. In addition, we find that the circuit facilitates coordinate secretion of both hormones in response to protein meals. Our study highlights the advantages of a paradoxical paracrine feedback loop in maintaining metabolic homeostasis.

Abstract Tissue stem cells produce a constant flux of differentiated cells with distinct proportions. Here, we show that stem cells in colonic crypts differentiate early to form precisely 1:3 ratio of secretory to absorptive cells. This precision is surprising, as there are only eight stem cells making irreversible fate decisions, and so large stochastic effects of this small pool should have yielded much larger noise in cell proportions. We use single molecule FISH, lineage‐tracing mice and simulations to identify the homeostatic mechanisms facilitating robust proportions. We find that Delta‐Notch lateral inhibition operates in a restricted spatial zone to reduce initial noise in cell proportions. Increased dwell time and dispersive migration of secretory cells further averages additional variability added during progenitor divisions and breaks up continuous patches of same‐fate cells. These noise‐reducing mechanisms resolve the trade‐off between early commitment and robust differentiation and ensure spatially uniform spread of secretory cells. Our findings may apply to other cases where small progenitor pools expand to give rise to precise tissue cell proportions.

The human genome contains thousands of long noncoding RNAs which have been preserved by evolution, through their functions are poorly described. Here the authors show that NORAD binds the Pumilo homologues PUM1 and PUM2 to regulate mRNA levels of genes involved in chromosome segregation.

mRNA is thought to predominantly reside in the cytoplasm, where it is translated and eventually degraded. Although nuclear retention of mRNA has a regulatory potential, it is considered extremely rare in mammals. Here, to explore the extent of mRNA retention in metabolic tissues, we combine deep sequencing of nuclear and cytoplasmic RNA fractions with single-molecule transcript imaging in mouse beta cells, liver, and gut. We identify a wide range of protein-coding genes for which the levels of spliced polyadenylated mRNA are higher in the nucleus than in the cytoplasm. These include genes such as the transcription factor ChREBP, Nlrp6, Glucokinase, and Glucagon receptor. We demonstrate that nuclear retention of mRNA can efficiently buffer cytoplasmic transcript levels from noise that emanates from transcriptional bursts. Our study challenges the view that transcripts predominantly reside in the cytoplasm and reveals a role of the nucleus in dampening gene expression noise.

Cancer therapy exerts a strong selection pressure that shapes tumor evolution, yet our knowledge of how tumors change during treatment is limited. Here, we report the analysis of cellular heterogeneity for genetic and phenotypic features and their spatial distribution in breast tumors pre- and post-neoadjuvant chemotherapy. We found that intratumor genetic diversity was tumor-subtype specific, and it did not change during treatment in tumors with partial or no response. However, lower pretreatment genetic diversity was significantly associated with pathologic complete response. In contrast, phenotypic diversity was different between pre- and posttreatment samples. We also observed significant changes in the spatial distribution of cells with distinct genetic and phenotypic features. We used these experimental data to develop a stochastic computational model to infer tumor growth patterns and evolutionary dynamics. Our results highlight the importance of integrated analysis of genotypes and phenotypes of single cells in intact tissues to predict tumor evolution.

Host shutoff is a common strategy used by viruses to repress cellular mRNA translation and concomitantly allow the efficient translation of viral mRNAs. Here we use RNA-sequencing and ribosome profiling to explore the mechanisms that are being utilized by the Influenza A virus (IAV) to induce host shutoff. We show that viral transcripts are not preferentially translated and instead the decline in cellular protein synthesis is mediated by viral takeover on the mRNA pool. Our measurements also uncover strong variability in the levels of cellular transcripts reduction, revealing that short transcripts are less affected by IAV. Interestingly, these mRNAs that are refractory to IAV infection are enriched in cell maintenance processes such as oxidative phosphorylation. Furthermore, we show that the continuous oxidative phosphorylation activity is important for viral propagation. Our results advance our understanding of IAV-induced shutoff, and suggest a mechanism that facilitates the translation of genes with important housekeeping functions.

How are cells able to maintain constant levels of mRNA when the number of genes in a cell doubles ahead of cell division?

BACKGROUND: Necdin, a MAGE family protein expressed primarily in the nervous system, has been shown to interact with both nuclear and cytoplasmic proteins, but the mechanism of its nucleocytoplasmic transport are unknown. METHODOLOGY/PRINCIPAL FINDINGS: We carried out a large-scale interaction screen using necdin as a bait in the yeast RRS system, and found a wide range of potential interactors with different subcellular localizations, including over 60 new candidates for direct binding to necdin. Integration of these interactions into a comprehensive network revealed a number of coherent interaction modules, including a cytoplasmic module connecting to necdin through huntingtin-associated protein 1 (Hap1), dynactin and hip-1 protein interactor (Hippi); a nuclear P53 and Creb-binding-protein (Crebbp) module, connecting through Crebbp and WW domain-containing transcription regulator protein 1 (Wwtr1); and a nucleocytoplasmic transport module, connecting through transportins 1 and 2. We validated the necdin-transportin1 interaction and characterized a sequence motif in necdin that modulates karyopherin interaction. Surprisingly, a D234P necdin mutant showed enhanced binding to both transportin1 and importin β1. Finally, exclusion of necdin from the nucleus triggered extensive cell death. CONCLUSIONS/SIGNIFICANCE: These data suggest that necdin has multiple roles within protein complexes in different subcellular compartments, and indicate that it can utilize multiple karyopherin-dependent pathways to modulate its localization.

Fundamental aspects of embryonic and post-natal development, including maintenance of the mammalian female germline, are largely unknown. Here we employ a retrospective, phylogenetic-based method for reconstructing cell lineage trees utilizing somatic mutations accumulated in microsatellites, to study female germline dynamics in mice. Reconstructed cell lineage trees can be used to estimate lineage relationships between different cell types, as well as cell depth (number of cell divisions since the zygote). We show that, in the reconstructed mouse cell lineage trees, oocytes form clusters that are separate from hematopoietic and mesenchymal stem cells, both in young and old mice, indicating that these populations belong to distinct lineages. Furthermore, while cumulus cells sampled from different ovarian follicles are distinctly clustered on the reconstructed trees, oocytes from the left and right ovaries are not, suggesting a mixing of their progenitor pools. We also observed an increase in oocyte depth with mouse age, which can be explained either by depth-guided selection of oocytes for ovulation or by post-natal renewal. Overall, our study sheds light on substantial novel aspects of female germline preservation and development.

Stem cell dynamics in vivo are often being studied by lineage tracing methods. Our laboratory has previously developed a retrospective method for reconstructing cell lineage trees from somatic mutations accumulated in microsatellites. This method was applied here to explore different aspects of stem cell dynamics in the mouse colon without the use of stem cell markers. We first demonstrated the reliability of our method for the study of stem cells by confirming previously established facts, and then we addressed open questions. Our findings confirmed that colon crypts are monoclonal and that, throughout adulthood, the process of monoclonal conversion plays a major role in the maintenance of crypts. The absence of immortal strand mechanism in crypts stem cells was validated by the age-dependent accumulation of microsatellite mutations. In addition, we confirmed the positive correlation between physical and lineage proximity of crypts, by showing that the colon is separated into small domains that share a common ancestor. We gained new data demonstrating that colon epithelium is clustered separately from hematopoietic and other cell types, indicating that the colon is constituted of few progenitors and ruling out significant renewal of colonic epithelium from hematopoietic cells during adulthood. Overall, our study demonstrates the reliability of cell lineage reconstruction for the study of stem cell dynamics, and it further addresses open questions in colon stem cells. In addition, this method can be applied to study stem cell dynamics in other systems.

Myofiber cultures give rise to myogenic as well as to non-myogenic cells. Whether these myofiber-associated non-myogenic cells develop from resident stem cells that possess mesenchymal plasticity or from other stem cells such as mesenchymal stem cells (MSCs) remain unsolved. To address this question, we applied a method for reconstructing cell lineage trees from somatic mutations to MSCs and myogenic and non-myogenic cells from individual myofibers that were cultured at clonal density.Our analyses show that (i) in addition to myogenic progenitors, myofibers also harbor non-myogenic progenitors of a distinct, yet close, lineage; (ii) myofiber-associated non-myogenic and myogenic cells share the same muscle-bound primordial stem cells of a lineage distinct from bone marrow MSCs; (iii) these muscle-bound primordial stem-cells first part to individual muscles and then differentiate into myogenic and non-myogenic stem cells.

Cells need to allocate their limited resources to express a wide range of genes. To understand how Escherichia coli partitions its transcriptional resources between its different promoters, we employ a robotic assay using a comprehensive reporter strain library for E. coli to measure promoter activity on a genomic scale at high-temporal resolution and accuracy. This allows continuous tracking of promoter activity as cells change their growth rate from exponential to stationary phase in different media. We find a heavy-tailed distribution of promoter activities, with promoter activities spanning several orders of magnitude. While the shape of the distribution is almost completely independent of the growth conditions, the identity of the promoters expressed at different levels does depend on them. Translation machinery genes, however, keep the same relative expression levels in the distribution across conditions, and their fractional promoter activity tracks growth rate tightly. We present a simple optimization model for resource allocation which suggests that the observed invariant distributions might maximize growth rate. These invariant features of the distribution of promoter activities may suggest design constraints that shape the allocation of transcriptional resources.

The depth of a cell of a multicellular organism is the number of cell divisions it underwent since the zygote, and knowing this basic cell property would help address fundamental problems in several areas of biology. At present, the depths of the vast majority of human and mouse cell types are unknown. Here, we show a method for estimating the depth of a cell by analyzing somatic mutations in its microsatellites, and provide to our knowledge for the first time reliable depth estimates for several cells types in mice. According to our estimates, the average depth of oocytes is 29, consistent with previous estimates. The average depth of B cells ranges from 34 to 79, linearly related to the mouse age, suggesting a rate of one cell division per day. In contrast, various types of adult stem cells underwent on average fewer cell divisions, supporting the notion that adult stem cells are relatively quiescent. Our method for depth estimation opens a window for revealing tissue turnover rates in animals, including humans, which has important implications for our knowledge of the body under physiological and pathological conditions.

BackgroundPhylogenetic mitochondrial DNA haplogroups are highly partitioned across global geographic regions. A unique exception is the X haplogroup, which has a widespread global distribution without major regions of distinct localization.Principal findingsWe have examined mitochondrial DNA sequence variation together with Y-chromosome-based haplogroup structure among the Druze, a religious minority with a unique socio-demographic history residing in the Near East. We observed a striking overall pattern of heterogeneous parental origins, consistent with Druze oral tradition, together with both a high frequency and a high diversity of the mitochondrial DNA (mtDNA) X haplogroup within a confined regional subpopulation. Furthermore demographic modeling indicated low migration rates with nearby populations.ConclusionsThese findings were enabled through the use of a paternal kindred based sampling approach, and suggest that the Galilee Druze represent a population isolate, and that the combination of a high frequency and diversity of the mtDNA X haplogroup signifies a phylogenetic refugium, providing a sample snapshot of the genetic landscape of the Near East prior to the modern age.

Synaptic wiring of neurons in Caenorhabditis elegans is largely invariable between animals. It has been suggested that this feature stems from genetically encoded molecular markers that guide the neurons in the final stage of synaptic formation. Identifying these markers and unraveling the logic by which they direct synapse formation is a key challenge. Here, we address this task by constructing a probabilistic model that attempts to explain the neuronal connectivity diagram of C. elegans as a function of the expression patterns of its neurons. By only considering neuron pairs that are known to be connected by chemical or electrical synapses, we focus on the final stage of synapse formation, in which neurons identify their designated partners. Our results show that for many neurons the neuronal expression map of C. elegans can be used to accurately predict the subset of adjacent neurons that will be chosen as its postsynaptic partners. Notably, these predictions can be achieved using the expression patterns of only a small number of specific genes that interact in a combinatorial fashion.

The cell lineage tree of a multicellular organism represents its history of cell divisions from the very first cell, the zygote. A new method for high-resolution reconstruction of parts of such cell lineage trees was recently developed based on phylogenetic analysis of somatic mutations accumulated during normal development of an organism. In this study we apply this method in mice to reconstruct the lineage trees of distinct cell types. We address for the first time basic questions in developmental biology of higher organisms, namely what is the correlation between the lineage relation among cells and their (1) function, (2) physical proximity and (3) anatomical proximity. We analyzed B-cells, kidney-, mesenchymal- and hematopoietic-stem cells, as well as satellite cells, which are adult skeletal muscle stem cells isolated from their niche on the muscle fibers (myofibers) from various skeletal muscles. Our results demonstrate that all analyzed cell types are intermingled in the lineage tree, indicating that none of these cell types are single exclusive clones. We also show a significant correlation between the physical proximity of satellite cells within muscles and their lineage. Furthermore, we show that satellite cells obtained from a single myofiber are significantly clustered in the lineage tree, reflecting their common developmental origin. Lineage analysis based on somatic mutations enables performing high resolution reconstruction of lineage trees in mice and humans, which can provide fundamental insights to many aspects of their development and tissue maintenance.

Trinucleotide hereditary diseases such as Huntington disease and Friedreich ataxia are cureless diseases associated with inheriting an abnormally large number of DNA trinucleotide repeats in a gene. The genes associated with different diseases are unrelated and harbor a trinucleotide repeat in different functional regions; therefore, it is striking that many of these diseases have similar correlations between their genotype, namely the number of inherited repeats and age of onset and progression phenotype. These correlations remain unexplained despite more than a decade of research. Although mechanisms have been proposed for several trinucleotide diseases, none of the proposals, being disease-specific, can account for the commonalities among these diseases. Here, we propose a universal mechanism in which length-dependent somatic repeat expansion occurs during the patient's lifetime toward a pathological threshold. Our mechanism uniformly explains for the first time to our knowledge the genotype-phenotype correlations common to trinucleotide disease and is well-supported by both experimental and clinical data. In addition, mathematical analysis of the mechanism provides simple explanations to a wide range of phenomena such as the exponential decrease of the age-of-onset curve, similar onset but faster progression in patients with Huntington disease with homozygous versus heterozygous mutation, and correlation of age of onset with length of the short allele but not with the long allele in Friedreich ataxia. If our proposed universal mechanism proves to be the core component of the actual mechanisms of specific trinucleotide diseases, it would open the search for a uniform treatment for all these diseases, possibly by delaying the somatic expansion process.