Cell Type Purification by Single-Cell Transcriptome-Trained Sorting.

Much of current molecular and cell biology research relies on the ability to purify cell types by fluorescence-activated cell sorting (FACS). FACS typically relies on the ability to label cell types of interest with antibodies or fluorescent transgenic constructs. However, antibody availability is often limited, and genetic manipulation is labor intensive or impossible in the case of primary human tissue. To date, no systematic method exists to enrich for cell types without a priori knowledge of cell-type markers. Here, we propose GateID, a computational method that combines single-cell transcriptomics with FACS index sorting to purify cell types of choice using only native cellular properties such as cell size, granularity, and mitochondrial content. We validate GateID by purifying various cell types from zebrafish kidney marrow and the human pancreas to high purity without resorting to specific antibodies or transgenes.

Differentiated Troy+ chief cells act as reserve stem cells to generate all lineages of the stomach epithelium.

Proliferation of the self-renewing epithelium of the gastric corpus occurs almost exclusively in the isthmus of the glands, from where cells migrate bidirectionally toward pit and base. The isthmus is therefore generally viewed as the stem cell zone. We find that the stem cell marker Troy is expressed at the gland base by a small subpopulation of fully differentiated chief cells. By lineage tracing with a Troy-eGFP-ires-CreERT2 allele, single marked chief cells are shown to generate entirely labeled gastric units over periods of months. This phenomenon accelerates upon tissue damage. Troy(+) chief cells can be cultured to generate long-lived gastric organoids. Troy marks a specific subset of chief cells that display plasticity in that they are capable of replenishing entire gastric units, essentially serving as quiescent "reserve" stem cells. These observations challenge the notion that stem cell hierarchies represent a "one-way street."

Long-term expanding human airway organoids for disease modeling.

Organoids are self-organizing 3D structures grown from stem cells that recapitulate essential aspects of organ structure and function. Here, we describe a method to establish long-term-expanding human airway organoids from broncho-alveolar resections or lavage material. The pseudostratified airway organoids consist of basal cells, functional multi-ciliated cells, mucus-producing secretory cells, and CC10-secreting club cells. Airway organoids derived from cystic fibrosis (CF) patients allow assessment of CFTR function in an organoid swelling assay. Organoids established from lung cancer resections and metastasis biopsies retain tumor histopathology as well as cancer gene mutations and are amenable to drug screening. Respiratory syncytial virus (RSV) infection recapitulates central disease features, dramatically increases organoid cell motility via the non-structural viral NS2 protein, and preferentially recruits neutrophils upon co-culturing. We conclude that human airway organoids represent versatile models for the in vitro study of hereditary, malignant, and infectious pulmonary disease.

Mapping the physical network of cellular interactions.

A cell's function is influenced by the environment, or niche, in which it resides. Studies of niches usually require assumptions about the cell types present, which impedes the discovery of new cell types or interactions. Here we describe ProximID, an approach for building a cellular network based on physical cell interaction and single-cell mRNA sequencing, and show that it can be used to discover new preferential cellular interactions without prior knowledge of component cell types. ProximID found specific interactions between megakaryocytes and mature neutrophils and between plasma cells and myeloblasts and/or promyelocytes (precursors of neutrophils) in mouse bone marrow, and it identified a Tac1+ enteroendocrine cell-Lgr5+ stem cell interaction in small intestine crypts. This strategy can be used to discover new niches or preferential interactions in a variety of organs.

Single-cell messenger RNA sequencing reveals rare intestinal cell types.

Understanding the development and function of an organ requires the characterization of all of its cell types. Traditional methods for visualizing and isolating subpopulations of cells are based on messenger RNA or protein expression of only a few known marker genes. The unequivocal identification of a specific marker gene, however, poses a major challenge, particularly if this cell type is rare. Identifying rare cell types, such as stem cells, short-lived progenitors, cancer stem cells, or circulating tumour cells, is crucial to acquire a better understanding of normal or diseased tissue biology. To address this challenge we first sequenced the transcriptome of hundreds of randomly selected cells from mouse intestinal organoids, cultured self-organizing epithelial structures that contain all cell lineages of the mammalian intestine. Organoid buds, like intestinal crypts, harbour stem cells that continuously differentiate into a variety of cell types, occurring at widely different abundances. Since available computational methods can only resolve more abundant cell types, we developed RaceID, an algorithm for rare cell type identification in complex populations of single cells. We demonstrate that this algorithm can resolve cell types represented by only a single cell in a population of randomly sampled organoid cells. We use this algorithm to identify Reg4 as a novel marker for enteroendocrine cells, a rare population of hormone-producing intestinal cells. Next, we use Reg4 expression to enrich for these rare cells and investigate the heterogeneity within this population. RaceID confirmed the existence of known enteroendocrine lineages, and moreover discovered novel subtypes, which we subsequently validated in vivo. Having validated RaceID we then applied the algorithm to ex vivo-isolated Lgr5-positive stem cells and their direct progeny. We find that Lgr5-positive cells represent a homogenous abundant population of stem cells mixed with a rare population of Lgr5-positive secretory cells. We envision broad applicability of our method for discovering rare cell types and the corresponding marker genes in healthy and diseased organs.

Reg4+ deep crypt secretory cells function as epithelial niche for Lgr5+ stem cells in colon.

Leucine-rich repeat-containing G-protein coupled receptor 5-positive (Lgr5(+)) stem cells reside at crypt bottoms of the small and large intestine. Small intestinal Paneth cells supply Wnt3, EGF, and Notch signals to neighboring Lgr5(+) stem cells. Whereas the colon lacks Paneth cells, deep crypt secretory (DCS) cells are intermingled with Lgr5(+) stem cells at crypt bottoms. Here, we report regenerating islet-derived family member 4 (Reg4) as a marker of DCS cells. To investigate a niche function, we eliminated DCS cells by using the diphtheria-toxin receptor gene knocked into the murine Reg4 locus. Ablation of DCS cells results in loss of stem cells from colonic crypts and disrupts gut homeostasis and colon organoid growth. In agreement, sorted Reg4(+) DCS cells promote organoid formation of single Lgr5(+) colon stem cells. DCS cells can be massively produced from Lgr5(+) colon stem cells in vitro by combined Notch inhibition and Wnt activation. We conclude that Reg4(+) DCS cells serve as Paneth cell equivalents in the colon crypt niche.

Self-organizing models of human trunk organogenesis recapitulate spinal cord and spine co-morphogenesis.

Integrated in vitro models of human organogenesis are needed to elucidate the multi-systemic events underlying development and disease. Here we report the generation of human trunk-like structures that model the co-morphogenesis, patterning and differentiation of the human spine and spinal cord. We identified differentiation conditions for human pluripotent stem cells favoring the formation of an embryo-like extending antero-posterior (AP) axis. Single-cell and spatial transcriptomics show that somitic and spinal cord differentiation trajectories organize along this axis and can self-assemble into a neural tube surrounded by somites upon extracellular matrix addition. Morphogenesis is coupled with AP patterning mechanisms, which results, at later stages of organogenesis, in in vivo-like arrays of neural subtypes along a neural tube surrounded by spine and muscle progenitors contacted by neuronal projections. This integrated system of trunk development indicates that in vivo-like multi-tissue co-morphogenesis and topographic organization of terminal cell types can be achieved in human organoids, opening windows for the development of more complex models of organogenesis.

Enteropathogenic Escherichia coli and vaccinia virus do not require the family of WASP-interacting proteins for pathogen-induced actin assembly.

The human pathogens enteropathogenic Escherichia coli (EPEC) and vaccinia virus trigger actin assembly in host cells by activating the host adaptor Nck and the actin nucleation promoter neural Wiskott-Aldrich syndrome protein (N-WASP). EPEC translocates effector molecules into host cells via type III secretion, and the interaction between the translocated intimin receptor (Tir) and the bacterial membrane protein intimin stimulates Nck and N-WASP recruitment, leading to the formation of actin pedestals beneath adherent bacteria. Vaccinia virus also recruits Nck and N-WASP to generate actin tails that promote cell-to-cell spread of the virus. In addition to Nck and N-WASP, WASP-interacting protein (WIP) localizes to vaccinia virus tails, and inhibition of actin tail formation upon ectopic expression of WIP mutants led to the suggestion that WIP is required for this process. Similar studies of WIP mutants, however, did not affect the ability of EPEC to form actin pedestals, arguing against an essential role for WIP in EPEC-induced actin assembly. In this study, we demonstrate that Nck and N-WASP are normally recruited by vaccinia virus and EPEC in the absence of WIP, and neither WIP nor the WIP family members CR16 and WIRE/WICH are essential for pathogen induced actin assembly. In addition, although Nck binds EPEC Tir directly, N-WASP is required for its localization during pedestal formation. Overall, these data highlight similar pathogenic strategies shared by EPEC and vaccinia virus by demonstrating a requirement for both Nck and N-WASP, but not WIP or WIP family members in pathogen-induced actin assembly.

Integration of multiple lineage measurements from the same cell reconstructs parallel tumor evolution.

Organoid evolution models complemented with integrated single-cell sequencing technology provide a powerful platform to characterize intra-tumor heterogeneity (ITH) and tumor evolution. Here, we conduct a parallel evolution experiment to mimic the tumor evolution process by evolving a colon cancer organoid model over 100 generations, spanning 6 months in time. We use single-cell whole-genome sequencing (WGS) in combination with viral lineage tracing at 12 time points to simultaneously monitor clone size, CNV states, SNV states, and viral lineage barcodes for 1,641 single cells. We integrate these measurements to construct clonal evolution trees with high resolution. We characterize the order of events in which chromosomal aberrations occur and identify aberrations that recur multiple times within the same tumor sub-population. We observe recurrent sequential loss of chromosome 4 after loss of chromosome 18 in four unique tumor clones. SNVs and CNVs identified in our organoid experiments are also frequently reported in colorectal carcinoma samples, and out of 334 patients with chromosome 18 loss in a Memorial Sloan Kettering colorectal cancer cohort, 99 (29.6%) also harbor chromosome 4 loss. Our study reconstructs tumor evolution in a colon cancer organoid model at high resolution, demonstrating an approach to identify potentially clinically relevant genomic aberrations in tumor evolution.

Patched Receptors Sense, Interpret, and Establish an Epidermal Hedgehog Signaling Gradient.

By using the sensitivity of single-molecule fluorescent in situ hybridization, we have precisely quantified the levels and defined the temporal and spatial distribution of Hedgehog signaling activity during embryonic skin development and discovered that there is a Hedgehog signaling gradient along the proximal-distal axis of developing hair follicles. To explore the contribution of Hedgehog receptors Ptch1 and Ptch2 in establishing the epidermal signaling gradient, we quantitated the level of pathway activity generated in Ptch1- and Ptch1;Ptch2-deficient skin and defined the contribution of each receptor to regulation of the levels of Hedgehog signaling identified in wild-type skin. Moreover, we show that both the cellular phenotype and level of pathway activity featured in Ptch1;Ptch2-deficient cells faithfully recapitulates the Peak level of endogenous Hedgehog signaling detected at the base of developing follicles, where the concentration of endogenous Shh is predicted to be highest. Taken together, these data show that both Ptch1 and Ptch2 play a crucial role in sensing the concentration of Hedgehog ligand and regulating the appropriate dose-dependent response.

A Single-Cell RNA Sequencing Study Reveals Cellular and Molecular Dynamics of the Hippocampal Neurogenic Niche.

Adult neurogenesis in the murine dentate gyrus occurs in a specialized microenvironment that sustains the generation of neurons during life. To fully understand adult neurogenesis, it is essential to determine the neural stem cell (NSC) and progenitor developmental stages, their molecular determinants, and the niche cellular and molecular composition. We report on a single-cell RNA sequencing study of the hippocampal niche, performed by isolating all the non-neuronal cell populations. Our analysis provides a comprehensive description of the dentate gyrus cells, and it allows the identification of exclusive cell-type-specific markers. We define the developmental stages and transcriptional dynamics of NSCs and progenitors, and we find that, while NSCs represent a heterogeneous cellular continuum, progenitors can be grouped into distinct subtypes. We determine the oligodendrocyte lineage and transcriptional dynamics, and we describe the microglia transcriptional profile and activation state. The combined data constitute a valuable resource to understand regulatory mechanisms of adult neurogenesis.

De Novo Prediction of Stem Cell Identity using Single-Cell Transcriptome Data.

Adult mitotic tissues like the intestine, skin, and blood undergo constant turnover throughout the life of an organism. Knowing the identity of the stem cell is crucial to understanding tissue homeostasis and its aberrations upon disease. Here we present a computational method for the derivation of a lineage tree from single-cell transcriptome data. By exploiting the tree topology and the transcriptome composition, we establish StemID, an algorithm for identifying stem cells among all detectable cell types within a population. We demonstrate that StemID recovers two known adult stem cell populations, Lgr5+ cells in the small intestine and hematopoietic stem cells in the bone marrow. We apply StemID to predict candidate multipotent cell populations in the human pancreas, a tissue with largely uncharacterized turnover dynamics. We hope that StemID will accelerate the search for novel stem cells by providing concrete markers for biological follow-up and validation.

Single-molecule mRNA detection and counting in mammalian tissue.

We present a protocol for visualizing and quantifying single mRNA molecules in mammalian (mouse and human) tissues. In the approach described here, sets of about 50 short oligonucleotides, each labeled with a single fluorophore, are hybridized to target mRNAs in tissue sections. Each set binds to a single mRNA molecule and can be detected by fluorescence microscopy as a diffraction-limited spot. Tissue architecture is then assessed by counterstaining the sections with DNA dye (DAPI), and cell borders can be visualized with a dye-coupled antibody. Spots are detected automatically with custom-made software, which we make freely available. The mRNA molecules thus detected are assigned to single cells within a tissue semiautomatically by using a graphical user interface developed in our laboratory. In this protocol, we describe an example of quantitative analysis of mRNA levels and localization in mouse small intestine. The procedure (from tissue dissection to obtaining data sets) takes 3 d. Data analysis will require an additional 3-7 d, depending on the type of analysis.

Dll1+ secretory progenitor cells revert to stem cells upon crypt damage.

Lgr5+ intestinal stem cells generate enterocytes and secretory cells. Secretory lineage commitment requires Notch silencing. The Notch ligand Dll1 is expressed by a subset of immediate stem cell daughters. Lineage tracing in Dll1(GFP-ires-CreERT2) knock-in mice reveals that single Dll1(high) cells generate small, short-lived clones containing all four secretory cell types. Lineage specification thus occurs in immediate stem cell daughters through Notch lateral inhibition. Cultured Dll1(high) cells form long-lived organoids (mini-guts) on brief Wnt3A exposure. When Dll1(high) cells are genetically marked before tissue damage, stem cell tracing events occur. Thus, secretory progenitors exhibit plasticity by regaining stemness on damage.

Single-molecule transcript counting of stem-cell markers in the mouse intestine.

Determining the molecular identities of adult stem cells requires technologies for sensitive transcript detection in tissues. In mouse intestinal crypts, lineage-tracing studies indicated that different genes uniquely mark spatially distinct stem-cell populations, residing either at crypt bases or at position +4, but a detailed analysis of their spatial co-expression has not been feasible. Here we apply three-colour single-molecule fluorescent in situ hybridization to study a comprehensive panel of intestinal stem-cell markers during homeostasis, ageing and regeneration. We find that the expression of all markers overlaps at crypt-base cells. This co-expression includes Lgr5, Bmi1 and mTert, genes previously suggested to mark distinct stem cells. Strikingly, Dcamkl1 tuft cells, distributed throughout the crypt axis, co-express Lgr5 and other stem-cell markers that are otherwise confined to crypt bases. We also detect significant changes in the expression of some of the markers following irradiation, indicating their potential role in the regeneration process. Our approach can enable the sensitive detection of putative stem cells in other tissues and in tumours, guiding complementary functional studies to evaluate their stem-cell properties.

Neural Wiskott-Aldrich syndrome protein modulates Wnt signaling and is required for hair follicle cycling in mice.

The Rho family GTPases Cdc42 and Rac1 are critical regulators of the actin cytoskeleton and are essential for skin and hair function. Wiskott-Aldrich syndrome family proteins act downstream of these GTPases, controlling actin assembly and cytoskeletal reorganization, but their role in epithelial cells has not been characterized in vivo. Here, we used a conditional knockout approach to assess the role of neural Wiskott-Aldrich syndrome protein (N-WASP), the ubiquitously expressed Wiskott-Aldrich syndrome-like (WASL) protein, in mouse skin. We found that N-WASP deficiency in mouse skin led to severe alopecia, epidermal hyperproliferation, and ulceration, without obvious effects on epidermal differentiation and wound healing. Further analysis revealed that the observed alopecia was likely the result of a progressive and ultimately nearly complete block in hair follicle (HF) cycling by 5 months of age. N-WASP deficiency also led to abnormal proliferation of skin progenitor cells, resulting in their depletion over time. Furthermore, N-WASP deficiency in vitro and in vivo correlated with decreased GSK-3beta phosphorylation, decreased nuclear localization of beta-catenin in follicular keratinocytes, and decreased Wnt-dependent transcription. Our results indicate a critical role for N-WASP in skin function and HF cycling and identify a link between N-WASP and Wnt signaling. We therefore propose that N-WASP acts as a positive regulator of beta-catenin-dependent transcription, modulating differentiation of HF progenitor cells.

The NF2 tumor suppressor Merlin and the ERM proteins interact with N-WASP and regulate its actin polymerization function.

The function of the NF2 tumor suppressor merlin has remained elusive despite increasing evidence for its role in actin cytoskeleton reorganization. The closely related ERM proteins (ezrin, radixin, and moesin) act as linkers between the cell membrane and cytoskeleton, and have also been implicated as active actin reorganizers. We report here that merlin and the ERMs can interact with and regulate N-WASP, a critical regulator of actin dynamics. Merlin and moesin were found to inhibit N-WASP-mediated actin assembly in vitro, a function that appears independent of their ability to bind actin. Furthermore, exogenous expression of a constitutively active ERM inhibits N-WASP-dependent Shigella tail formation, suggesting that the ERMs may function as inhibitors of N-WASP function in vivo. This novel function of merlin and the ERMs illustrates a mechanism by which these proteins directly exert their effects on actin reorganization and also provides new insight into N-WASP regulation.