Pancreatic beta cells have been shown to be heterogeneous at multiple levels. However, spatially interrogating transcriptional heterogeneity in the intact tissue has been challenging. Here, we developed an optimized protocol for single-molecule transcript imaging in the intact pancreas and used it to identify a sub-population of "extreme" beta cells with elevated mRNA levels of insulin and other secretory genes. Extreme beta cells contain higher ribosomal and proinsulin content but lower levels of insulin protein in fasted states, suggesting they may be tuned for basal insulin secretion. They exhibit a distinctive intra-cellular polarization pattern, with elevated mRNA concentrations in an apical ER-enriched compartment, distinct from the localization of nascent and mature proteins. The proportion of extreme cells increases in db/db diabetic mice, potentially facilitating the required increase in basal insulin. Our results thus highlight a sub-population of beta cells that may carry distinct functional roles along physiological and pathological timescales.

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.

Understanding the dynamics and variability of protein circuitry requires accurate measurements in living cells as well as theoretical models. To address this, we employed one of the best-studied protein circuits in human cells, the negative feedback loop between the tumor suppressor p53 and the oncogene Mdm2. We measured the dynamics of fluorescently tagged p53 and Mdm2 over several days in individual living cells. We found that isogenic cells in the same environment behaved in highly variable ways following DNA-damaging gamma irradiation: some cells showed undamped oscillations for at least 3 days (more than 10 peaks). The amplitude of the oscillations was much more variable than the period. Sister cells continued to oscillate in a correlated way after cell division, but lost correlation after about 11 h on average. Other cells showed low-frequency fluctuations that did not resemble oscillations. We also analyzed different families of mathematical models of the system, including a novel checkpoint mechanism. The models point to the possible source of the variability in the oscillations: low-frequency noise in protein production rates, rather than noise in other parameters such as degradation rates. This study provides a view of the extensive variability of the behavior of a protein circuit in living human cells, both from cell to cell and in the same cell over time.

Neutrophils provide first line of host defense against bacterial infections utilizing glycolysis for their effector functions. How glycolysis and its major byproduct lactate are triggered in bone marrow (BM) neutrophils and their contribution to neutrophil mobilization in acute inflammation is not clear. Here we report that bacterial lipopolysaccharides (LPS) or Salmonella Typhimurium triggers lactate release by increasing glycolysis, NADPH-oxidase-mediated reactive oxygen species and HIF-1α levels in BM neutrophils. Increased release of BM lactate preferentially promotes neutrophil mobilization by reducing endothelial VE-Cadherin expression, increasing BM vascular permeability via endothelial lactate-receptor GPR81 signaling. GPR81-/- mice mobilize reduced levels of neutrophils in response to LPS, unless rescued by VE-Cadherin disrupting antibodies. Lactate administration also induces release of the BM neutrophil mobilizers G-CSF, CXCL1 and CXCL2, indicating that this metabolite drives neutrophil mobilization via multiple pathways. Our study reveals a metabolic crosstalk between lactate-producing neutrophils and BM endothelium, which controls neutrophil mobilization under bacterial infection.

The islets of Langerhans are dynamic structures that can change in size, number of cells, and molecular function in response to physiological and pathological stress. Molecular cues originating from the surrounding "peri-islet" acinar cells that could facilitate this plasticity have not been explored. Here, we combine single-molecule transcript imaging in the intact pancreas and transcriptomics to identify spatial heterogeneity of acinar cell gene expression. We find that peri-islet acinar cells exhibit a distinct molecular signature in db/db diabetic mice that includes upregulation of trypsin family genes and elevated mTOR activity. This zonated expression program seems to be induced by CCK that is secreted from islet cells. Elevated peri-islet trypsin secretion could facilitate the islet expansion observed in this model via modulation of the islet capsule matrix components. Our study highlights a molecular axis of communication between the pancreatic exocrine and endocrine compartments that may be relevant to islet expansion.

Single-cell transcriptomics (scRNA-seq) has revolutionized the understanding of the spatial architecture of tissue structure and function. Advancing the "transcript-centric" view of scRNA-seq analyses is presently restricted by the limited resolution of proteomics and genome-wide techniques to analyze post-translational modifications. Here, by combining spatial cell sorting with transcriptomics and quantitative proteomics/phosphoproteomics, we established the spatially resolved proteome landscape of the liver endothelium, yielding deep mechanistic insight into zonated vascular signaling mechanisms. Phosphorylation of receptor tyrosine kinases was detected preferentially in the central vein area, resulting in an atypical enrichment of tyrosine phosphorylation. Prototypic biological validation identified Tie receptor signaling as a selective and specific regulator of vascular Wnt activity orchestrating angiocrine signaling, thereby controlling hepatocyte function during liver regeneration. Taken together, the study has yielded fundamental insight into the spatial organization of liver endothelial cell signaling. Spatial sorting may be employed as a universally adaptable strategy for multiomic analyses of scRNA-seq-defined cellular (sub)-populations.

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.

The intestinal epithelium is a highly structured tissue composed of repeating crypt-villus units. Enterocytes perform the diverse tasks of absorbing a wide range of nutrients while protecting the body from the harsh bacterium-rich environment. It is unknown whether these tasks are spatially zonated along the villus axis. Here, we extracted a large panel of landmark genes characterized by transcriptomics of laser capture microdissected villus segments and utilized it for single-cell spatial reconstruction, uncovering broad zonation of enterocyte function along the villus. We found that enterocytes at villus bottoms express an anti-bacterial gene program in a microbiome-dependent manner. They next shift to sequential expression of carbohydrates, peptides, and fat absorption machineries in distinct villus compartments. Finally, they induce a Cd73 immune-modulatory program at the villus tips. Our approach can be used to uncover zonation patterns in other organs when prior knowledge of landmark genes is lacking.

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 cells hijack and remodel existing metabolic pathways for their benefit. Argininosuccinate synthase (ASS1) is a urea cycle enzyme that is essential in the conversion of nitrogen from ammonia and aspartate to urea. A decrease in nitrogen flux through ASS1 in the liver causes the urea cycle disorder citrullinaemia. In contrast to the well-studied consequences of loss of ASS1 activity on ureagenesis, the purpose of its somatic silencing in multiple cancers is largely unknown. Here we show that decreased activity of ASS1 in cancers supports proliferation by facilitating pyrimidine synthesis via CAD (carbamoyl-phosphate synthase 2, aspartate transcarbamylase, and dihydroorotase complex) activation. Our studies were initiated by delineating the consequences of loss of ASS1 activity in humans with two types of citrullinaemia. We find that in citrullinaemia type I (CTLN I), which is caused by deficiency of ASS1, there is increased pyrimidine synthesis and proliferation compared with citrullinaemia type II (CTLN II), in which there is decreased substrate availability for ASS1 caused by deficiency of the aspartate transporter citrin. Building on these results, we demonstrate that ASS1 deficiency in cancer increases cytosolic aspartate levels, which increases CAD activation by upregulating its substrate availability and by increasing its phosphorylation by S6K1 through the mammalian target of rapamycin (mTOR) pathway. Decreasing CAD activity by blocking citrin, the mTOR signalling, or pyrimidine synthesis decreases proliferation and thus may serve as a therapeutic strategy in multiple cancers where ASS1 is downregulated. Our results demonstrate that ASS1 downregulation is a novel mechanism supporting cancerous proliferation, and they provide a metabolic link between the urea cycle enzymes and pyrimidine synthesis.

The intestinal epithelium holds an immense regenerative capacity mobilized by intestinal stem cells (ISCs), much of it supported by Wnt pathway activation. Several unique regulatory mechanisms ensuring optimal levels of Wnt signaling have been recognized in ISCs. Here, we identify another Wnt signaling amplifier, CKIε, which is specifically upregulated in ISCs and is essential for ISC maintenance, especially in the absence of its close isoform CKIδ. Co-ablation of CKIδ/ε in the mouse gut epithelium results in rapid ISC elimination, with subsequent growth arrest, crypt-villous shrinking, and rapid mouse death. Unexpectedly, Wnt activation is preserved in all CKIδ/ε-deficient enterocyte populations, with the exception of Lgr5+ ISCs, which exhibit Dvl2-dependent Wnt signaling attenuation. CKIδ/ε-depleted gut organoids cease proliferating and die rapidly, yet survive and resume self-renewal upon reconstitution of Dvl2 expression. Our study underscores a unique regulation mode of the Wnt pathway in ISCs, possibly providing new means of stem cell enrichment for regenerative medicine.

Innate lymphoid cells (ILCs) are critical modulators of mucosal immunity, inflammation, and tissue homeostasis, but their full spectrum of cellular states and regulatory landscapes remains elusive. Here, we combine genome-wide RNA-seq, ChIP-seq, and ATAC-seq to compare the transcriptional and epigenetic identity of small intestinal ILCs, identifying thousands of distinct gene profiles and regulatory elements. Single-cell RNA-seq and flow and mass cytometry analyses reveal compartmentalization of cytokine expression and metabolic activity within the three classical ILC subtypes and highlight transcriptional states beyond the current canonical classification. In addition, using antibiotic intervention and germ-free mice, we characterize the effect of the microbiome on the ILC regulatory landscape and determine the response of ILCs to microbial colonization at the single-cell level. Together, our work characterizes the spectrum of transcriptional identities of small intestinal ILCs and describes how ILCs differentially integrate signals from the microbial microenvironment to generate phenotypic and functional plasticity.

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.

The immune and enteric nervous (ENS) systems monitor the frontier with commensal and pathogenic microbes in the colon. We investigated whether FoxP3+ regulatory T (Treg) cells functionally interact with the ENS. Indeed, microbe-responsive RORγ+ and Helios+ subsets localized in close apposition to nitrergic and peptidergic nerve fibers in the colon lamina propria (LP). Enteric neurons inhibited in vitro Treg (iTreg) differentiation in a cell-contact-independent manner. A screen of neuron-secreted factors revealed a role for interleukin-6 (IL-6) in modulating iTreg formation and their RORγ+ proportion. Colonization of germfree mice with commensals, especially RORγ+ Treg inducers, broadly diminished colon neuronal density. Closing the triangle, conditional ablation of IL-6 in neurons increased total Treg cells but decreased the RORγ+ subset, as did depletion of two ENS neurotransmitters. Our findings suggest a regulatory circuit wherein microbial signals condition neuronal density and activation, thus tuning Treg cell generation and immunological tolerance in the gut.

sBLISS (in-suspension breaks labeling in situ and sequencing) is a versatile and widely applicable method for identification of endogenous and induced DNA double-strand breaks (DSBs) in any cell type that can be brought into suspension. sBLISS provides genome-wide profiles of the most consequential DNA lesion implicated in a variety of pathological, but also physiological, processes. In sBLISS, after in situ labeling, DSB ends are linearly amplified, followed by next-generation sequencing and DSB landscape analysis. Here, we present a step-by-step experimental protocol for sBLISS, as well as a basic computational analysis. The main advantages of sBLISS are (i) the suspension setup, which renders the protocol user-friendly and easily scalable; (ii) the possibility of adapting it to a high-throughput or single-cell workflow; and (iii) its flexibility and its applicability to virtually every cell type, including patient-derived cells, organoids, and isolated nuclei. The wet-lab protocol can be completed in 1.5 weeks and is suitable for researchers with intermediate expertise in molecular biology and genomics. For the computational analyses, basic-to-intermediate bioinformatics expertise is required.

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.

Lung development and function arises from the interactions between diverse cell types and lineages. Using single-cell RNA sequencing (RNA-seq), we characterize the cellular composition of the lung during development and identify vast dynamics in cell composition and their molecular characteristics. Analyzing 818 ligand-receptor interaction pairs within and between cell lineages, we identify broadly interacting cells, including AT2, innate lymphocytes (ILCs), and basophils. Using interleukin (IL)-33 receptor knockout mice and in vitro experiments, we show that basophils establish a lung-specific function imprinted by IL-33 and granulocyte-macrophage colony-stimulating factor (GM-CSF), characterized by unique signaling of cytokines and growth factors important for stromal, epithelial, and myeloid cell fates. Antibody-depletion strategies, diphtheria toxin-mediated selective depletion of basophils, and co-culture studies show that lung resident basophils are important regulators of alveolar macrophage development and function. Together, our study demonstrates how whole-tissue signaling interaction map on the single-cell level can broaden our understanding of cellular networks in health and disease.

The immune system generates pathogen-tailored responses. The precise innate immune cell types and pathways that direct robust adaptive immune responses have not been fully characterized. By using fluorescent pathogens combined with massively parallel single-cell RNA-seq, we comprehensively characterized the initial 48 h of the innate immune response to diverse pathogens. We found that across all pathogens tested, most of the lymph node cell types and states showed little pathogen specificity. In contrast, the rare antigen-positive cells displayed pathogen-specific transcriptional programs as early as 24 h after immunization. In addition, mycobacteria activated a specific NK-driven IFNγ response. Depletion of NK cells and IFNγ showed that IFNγ initiated a monocyte-specific signaling cascade, leading to the production of major chemokines and cytokines that promote Th1 development. Our systems immunology approach sheds light on early events in innate immune responses and may help further development of safe and efficient vaccines.

Thousands of long noncoding RNA (lncRNA) genes are encoded in the human genome, and hundreds of them are evolutionarily conserved, but their functions and modes of action remain largely obscure. Particularly enigmatic lncRNAs are those that are exported to the cytoplasm, including NORAD-an abundant and highly conserved cytoplasmic lncRNA. Here we show that most of the sequence of NORAD is comprised of repetitive units that together contain at least 17 functional binding sites for the two mammalian Pumilio homologues. Through binding to PUM1 and PUM2, NORAD modulates the mRNA levels of their targets, which are enriched for genes involved in chromosome segregation during cell division. Our results suggest that some cytoplasmic lncRNAs function by modulating the activities of RNA-binding proteins, an activity which positions them at key junctions of cellular signalling pathways.

Intestinal crypts in mammals are comprised of long-lived stem cells and shorter-lived progenies. These two populations are maintained in specific proportions during adult life. Here, we investigate the design principles governing the dynamics of these proportions during crypt morphogenesis. Using optimal control theory, we show that a proliferation strategy known as a "bang-bang" control minimizes the time to obtain a mature crypt. This strategy consists of a surge of symmetric stem cell divisions, establishing the entire stem cell pool first, followed by a sharp transition to strictly asymmetric stem cell divisions, producing nonstem cells with a delay. We validate these predictions using lineage tracing and single-molecule fluorescence in situ hybridization of intestinal crypts in infant mice, uncovering small crypts that are entirely composed of Lgr5-labeled stem cells, which become a minority as crypts continue to grow. Our approach can be used to uncover similar design principles in other developmental systems.