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.

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.

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.

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

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

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.

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.

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.

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.