The cellular states and fates of shed intestinal cells.

The intestinal epithelium is replaced every few days1. Enterocytes are shed into the gut lumen predominantly from the tips of villi2,3 and have been believed to rapidly die upon their dissociation from the tissue4,5. However, technical limitations prohibited studying the cellular states and fates of shed intestinal cells. Here we show that shed epithelial cells remain viable and upregulate distinct anti-microbial programmes upon shedding, using bulk and single-cell RNA sequencing of male mouse intestinal faecal washes. We further identify abundant shedding of immune cells, which is elevated in mice with dextran sulfate sodium-induced colitis. We find that faecal host transcriptomics reflect changes in the intestinal tissue following perturbations. Our study suggests potential functions of shed cells in the intestinal lumen and demonstrates that host cell transcriptomes in intestinal washes can be used to probe tissue states.

Distal Fecal Wash Host Transcriptomics Identifies Inflammation Throughout the Colon and Terminal Ileum.

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.

A spatiotemporally resolved single-cell atlas of the Plasmodium liver stage.

Malaria infection involves an obligatory, yet clinically silent liver stage1,2. Hepatocytes operate in repeating units termed lobules, exhibiting heterogeneous gene expression patterns along the lobule axis3, but the effects of hepatocyte zonation on parasite development at the molecular level remain unknown. Here we combine single-cell RNA sequencing4 and single-molecule transcript imaging5 to characterize the host and parasite temporal expression programmes in a zonally controlled manner for the rodent malaria parasite Plasmodium berghei ANKA. We identify differences in parasite gene expression in distinct zones, including potentially co-adaptive programmes related to iron and fatty acid metabolism. We find that parasites develop more rapidly in the pericentral lobule zones and identify a subpopulation of periportally biased hepatocytes that harbour abortive infections, reduced levels of Plasmodium transcripts and parasitophorous vacuole breakdown. These 'abortive hepatocytes', which appear predominantly with high parasite inoculum, upregulate immune recruitment and key signalling programmes. Our study provides a resource for understanding the liver stage of Plasmodium infection at high spatial resolution and highlights the heterogeneous behaviour of both the parasite and the host hepatocyte.

The spatiotemporal program of zonal liver regeneration following acute injury.

The liver carries a remarkable ability to regenerate rapidly after acute zonal damage. Single-cell approaches are necessary to study this process, given the spatial heterogeneity of liver cell types. Here, we use spatially resolved single-cell RNA sequencing (scRNA-seq) to study the dynamics of mouse liver regeneration after acute acetaminophen (APAP) intoxication. We find that hepatocytes proliferate throughout the liver lobule, creating the mitotic pressure required to repopulate the necrotic pericentral zone rapidly. A subset of hepatocytes located at the regenerating front transiently upregulate fetal-specific genes, including Afp and Cdh17, as they reprogram to a pericentral state. Zonated endothelial, hepatic stellate cell (HSC), and macrophage populations are differentially involved in immune recruitment, proliferation, and matrix remodeling. We observe massive transient infiltration of myeloid cells, yet stability of lymphoid cell abundance, in accordance with a global decline in antigen presentation. Our study provides a resource for understanding the coordinated programs of zonal liver regeneration.

Host transcriptome signatures in human faecal-washes predict histological remission in patients with IBD.

BACKGROUND: Colonoscopy is the gold standard for evaluation of inflammation in inflammatory bowel diseases (IBDs), yet entails cumbersome preparations and risks of injury. Existing non-invasive prognostic tools are limited in their diagnostic power. Moreover, transcriptomics of colonic biopsies have been inconclusive in their association with clinical features. AIMS: To assess the utility of host transcriptomics of faecal wash samples of patients with IBD compared with controls. METHODS: In this prospective cohort study, we obtained biopsies and faecal-wash samples from patients with IBD and controls undergoing lower endoscopy. We performed RNAseq of biopsies and matching faecal-washes, and associated them with endoscopic and histological inflammation status. We also performed faecal mass-spectrometry proteomics on a subset of samples. We inferred cell compositions using computational deconvolution and used classification algorithms to identify informative genes. RESULTS: We analysed biopsies and faecal washes from 39 patients (20 IBD, 19 controls). Host faecal-transcriptome carried information that was distinct from biopsy RNAseq and faecal proteomics. Transcriptomics of faecal washes, yet not of biopsies, from patients with histological inflammation were significantly correlated to one another (p=5.3×10-12). Faecal-transcriptome had significantly higher statistical power in identifying histological inflammation compared with transctiptome of intestinal biopsies (150 genes with area under the curve >0.9 in faecal samples vs 10 genes in biopsy RNAseq). These results were replicated in a validation cohort of 22 patients (10 IBD, 12 controls). Faecal samples were enriched in inflammatory monocytes, regulatory T cells, natural killer-cells and innate lymphoid cells. CONCLUSIONS: Faecal wash host transcriptome is a statistically powerful biomarker reflecting histological inflammation. Furthermore, it opens the way to identifying important correlates and therapeutic targets that may be obscured using biopsy transcriptomics.

Insulin is expressed by enteroendocrine cells during human fetal development.

Generation of beta cells via transdifferentiation of other cell types is a promising avenue for the treatment of diabetes. Here we reconstruct a single-cell atlas of the human fetal and neonatal small intestine. We identify a subset of fetal enteroendocrine K/L cells that express high levels of insulin and other beta cell genes. Our findings highlight a potential extra-pancreatic source of beta cells and expose its molecular blueprint.

Spatial gene expression maps of the intestinal lymphoid follicle and associated epithelium identify zonated expression programs.

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.

Clump sequencing exposes the spatial expression programs of intestinal secretory cells.

Single-cell RNA sequencing combined with spatial information on landmark genes enables reconstruction of spatially-resolved tissue cell atlases. However, such approaches are challenging for rare cell types, since their mRNA contents are diluted in the spatial transcriptomics bulk measurements used for landmark gene detection. In the small intestine, enterocytes, the most common cell type, exhibit zonated expression programs along the crypt-villus axis, but zonation patterns of rare cell types such as goblet and tuft cells remain uncharacterized. Here, we present ClumpSeq, an approach for sequencing small clumps of attached cells. By inferring the crypt-villus location of each clump from enterocyte landmark genes, we establish spatial atlases for all epithelial cell types in the small intestine. We identify elevated expression of immune-modulatory genes in villus tip goblet and tuft cells and heterogeneous migration patterns of enteroendocrine cells. ClumpSeq can be applied for reconstructing spatial atlases of rare cell types in other tissues and tumors.

Space-time logic of liver gene expression at sub-lobular scale.

The mammalian liver is a central hub for systemic metabolic homeostasis. Liver tissue is spatially structured, with hepatocytes operating in repeating lobules, and sub-lobule zones performing distinct functions. The liver is also subject to extensive temporal regulation, orchestrated by the interplay of the circadian clock, systemic signals and feeding rhythms. However, liver zonation has previously been analysed as a static phenomenon, and liver chronobiology has been analysed at tissue-level resolution. Here, we use single-cell RNA-seq to investigate the interplay between gene regulation in space and time. Using mixed-effect models of messenger RNA expression and smFISH validations, we find that many genes in the liver are both zonated and rhythmic, and most of them show multiplicative space-time effects. Such dually regulated genes cover not only key hepatic functions such as lipid, carbohydrate and amino acid metabolism, but also previously unassociated processes involving protein chaperones. Our data also suggest that rhythmic and localized expression of Wnt targets could be explained by rhythmically expressed Wnt ligands from non-parenchymal cells near the central vein. Core circadian clock genes are expressed in a non-zonated manner, indicating that the liver clock is robust to zonation. Together, our scRNA-seq analysis reveals how liver function is compartmentalized spatio-temporally at the sub-lobular scale.

The pro-oncogenic effect of the lncRNA H19 in the development of chronic inflammation-mediated hepatocellular carcinoma.

The oncofetal long noncoding RNA (lncRNA) H19 is postnatally repressed in most tissues, and re-expressed in many cancers, including hepatocellular carcinoma (HCC). The role of H19 in carcinogenesis is a subject of controversy. We aimed to examine the role of H19 in chronic inflammation-mediated hepatocarcinogenesis using the Mdr2/Abcb4 knockout (Mdr2-KO) mouse, a well-established HCC model. For this goal, we have generated Mdr2-KO/H19-KO double knockout (dKO) mice and followed spontaneous tumor development in the dKO and control Mdr2-KO mice. Cellular localization of H19 and effects of H19 loss in the liver were determined in young and old Mdr2-KO mice. Tumor incidence and tumor load were both significantly decreased in the liver of dKO versus Mdr2-KO females. The expression levels of H19 and Igf2 were variable in nontumor liver tissues of Mdr2-KO females and were significantly downregulated in most matched tumors. In nontumor liver tissue of aged Mdr2-KO females, H19 was expressed mainly in hepatocytes, and hepatocyte proliferation was increased compared to dKO females. At an early age, dKO females displayed lower levels of liver injury and B-cell infiltration, with higher percentage of binuclear hepatocytes. In human samples, H19 expression was higher in females, positively correlated with cirrhosis (in nontumor liver samples) and negatively correlated with CTNNB1 (beta-catenin) mutations and patients' survival (in tumors). Our data demonstrate that the lncRNA H19 is pro-oncogenic during the development of chronic inflammation-mediated HCC in the Mdr2-KO mouse model, mainly by increasing liver injury and decreasing hepatocyte polyploidy in young mice.

A single cell atlas of the human liver tumor microenvironment.

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.

Zonation of Pancreatic Acinar Cells in Diabetic Mice.

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.

Lgr5+ telocytes are a signaling source at the intestinal villus tip.

The intestinal epithelium is a structured organ composed of crypts harboring Lgr5+ stem cells, and villi harboring differentiated cells. Spatial transcriptomics have demonstrated profound zonation of epithelial gene expression along the villus axis, but the mechanisms shaping this spatial variability are unknown. Here, we combine laser capture micro-dissection and single cell RNA sequencing to uncover spatially zonated populations of mesenchymal cells along the crypt-villus axis. These include villus tip telocytes (VTTs) that express Lgr5, a gene previously considered a specific crypt epithelial stem cell marker. VTTs are elongated cells that line the villus tip epithelium and signal through Bmp morphogens and the non-canonical Wnt5a ligand. Their ablation is associated with perturbed zonation of enterocyte genes induced at the villus tip. Our study provides a spatially-resolved cell atlas of the small intestinal stroma and exposes Lgr5+ villus tip telocytes as regulators of the epithelial spatial expression programs along the villus axis.

Spatial sorting enables comprehensive characterization of liver zonation.

The mammalian liver is composed of repeating hexagonal units termed lobules. Spatially resolved single-cell transcriptomics revealed that about half of hepatocyte genes are differentially expressed across the lobule, yet technical limitations impeded reconstructing similar global spatial maps of other hepatocyte features. Here, we show how zonated surface markers can be used to sort hepatocytes from defined lobule zones with high spatial resolution. We apply transcriptomics, miRNA array measurements and mass spectrometry proteomics to reconstruct spatial atlases of multiple zonated features. We demonstrate that protein zonation largely overlaps with mRNA zonation, with the periportal HNF4α as an exception. We identify zonation of miRNAs such as miR-122, and inverse zonation of miRNAs and their hepatocyte target genes, highlighting potential regulation of protein levels through zonated mRNA degradation. Among the targets we find the pericentral Wnt receptors Fzd7 and Fzd8 and the periportal Wnt inhibitors Tcf7l1 and Ctnnbip1. Our approach facilitates reconstructing spatial atlases of multiple cellular features in the liver and other structured tissues.

Single-Cell Analysis of Diverse Pathogen Responses Defines a Molecular Roadmap for Generating Antigen-Specific Immunity.

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.

Transcriptional Heterogeneity of Beta Cells in the Intact Pancreas.

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.

Spatial Reconstruction of Single Enterocytes Uncovers Broad Zonation along the Intestinal Villus Axis.

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.

Induction of Nitric-Oxide Metabolism in Enterocytes Alleviates Colitis and Inflammation-Associated Colon Cancer.

Nitric oxide (NO) plays an established role in numerous physiological and pathological processes, but the specific cellular sources of NO in disease pathogenesis remain unclear, preventing the implementation of NO-related therapy. Argininosuccinate lyase (ASL) is the only enzyme able to produce arginine, the substrate for NO generation by nitric oxide synthase (NOS) isoforms. Here, we generated cell-specific conditional ASL knockout mice in combination with genetic and chemical colitis models. We demonstrate that NO derived from enterocytes alleviates colitis by decreasing macrophage infiltration and tissue damage, whereas immune cell-derived NO is associated with macrophage activation, resulting in increased severity of inflammation. We find that induction of endogenous NO production by enterocytes with supplements that upregulate ASL expression and complement its substrates results in improved epithelial integrity and alleviation of colitis and of inflammation-associated colon cancer.

Dynamic zonation of liver polyploidy.

The liver is a polyploid organ, consisting of hepatocytes with one or two nuclei each containing 2, 4, 8 or more haploid chromosome sets. The dynamic changes in the spatial distributions of polyploid classes across the liver lobule, its repeating anatomical unit, have not been characterized. Identifying these spatial patterns is important for understanding liver homeostatic and regenerative turnover, as well as potential division of labor among ploidy classes. Here, we use single molecule-based tissue imaging to reconstruct the spatial zonation profiles of liver polyploid classes in mice of different ages. We find that liver polyploidy proceeds in spatial waves, advancing more rapidly in the mid-lobule zone compared to the periportal and perivenous zones. We also measure the spatial zonation profiles of S-phase entry at different ages and identify more rapid S-phase entry in the mid-lobule zone at older ages. Our findings reveal fundamental features of liver spatial heterogeneity and highlight their dynamic changes during development and aging.

Microglia development follows a stepwise program to regulate brain homeostasis.

Microglia, the resident myeloid cells of the central nervous system, play important roles in life-long brain maintenance and in pathology. Despite their importance, their regulatory dynamics during brain development have not been fully elucidated. Using genome-wide chromatin and expression profiling coupled with single-cell transcriptomic analysis throughout development, we found that microglia undergo three temporal stages of development in synchrony with the brain--early, pre-, and adult microglia--which are under distinct regulatory circuits. Knockout of the gene encoding the adult microglia transcription factor MAFB and environmental perturbations, such as those affecting the microbiome or prenatal immune activation, led to disruption of developmental genes and immune response pathways. Together, our work identifies a stepwise microglia developmental program integrating immune response pathways that may be associated with several neurodevelopmental disorders.