Epithelial zonation along the mouse and human small intestine defines five discrete metabolic domains.

A key aspect of nutrient absorption is the exquisite division of labour across the length of the small intestine, with individual nutrients taken up at different proximal:distal positions. For millennia, the small intestine was thought to comprise three segments with indefinite borders: the duodenum, jejunum and ileum. By examining the fine-scale longitudinal transcriptional patterns that span the mouse and human small intestine, we instead identified five domains of nutrient absorption that mount distinct responses to dietary changes, and three regional stem cell populations. Molecular domain identity can be detected with machine learning, which provides a systematic method to computationally identify intestinal domains in mice. We generated a predictive model of transcriptional control of domain identity and validated the roles of Ppar-δ and Cdx1 in patterning lipid metabolism-associated genes. These findings represent a foundational framework for the zonation of absorption across the mammalian small intestine.

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.

Epithelial zonation along the mouse and human small intestine defines five discrete metabolic domains.

A key aspect of nutrient absorption is the exquisite division of labor across the length of the small intestine, with individual classes of micronutrients taken up at different positions. For millennia, the small intestine was thought to comprise three segments with indefinite borders: the duodenum, jejunum, and ileum. By examining fine-scale longitudinal segmentation of the mouse and human small intestines, we identified transcriptional signatures and upstream regulatory factors that define five domains of nutrient absorption, distinct from the three traditional sections. Spatially restricted expression programs were most prominent in nutrient-absorbing enterocytes but initially arose in intestinal stem cells residing in three regional populations. While a core signature was maintained across mice and humans with different diets and environments, domain properties were influenced by dietary changes. We established the functions of Ppar-ẟ and Cdx1 in patterning lipid metabolism in distal domains and generated a predictive model of additional transcription factors that direct domain identity. Molecular domain identity can be detected with machine learning, representing the first systematic method to computationally identify specific intestinal regions in mice. These findings provide a foundational framework for the identity and control of longitudinal zonation of absorption along the proximal:distal small intestinal axis.

Single-cell biology: what does the future hold?

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

Single-cell atlas of the human neonatal small intestine affected by necrotizing enterocolitis.

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.

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.

Physically interacting beta-delta pairs in the regenerating pancreas revealed by single-cell sequencing.

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.

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 discordances between mRNAs and proteins in the intestinal epithelium.

The use of transcriptomes as reliable proxies for cellular proteomes is controversial. In the small intestine, enterocytes operate for 4 days as they migrate along villi, which are highly graded microenvironments. Spatial transcriptomics have demonstrated profound zonation in enterocyte gene expression, but how this variability translates to protein content is unclear. Here we show that enterocyte proteins and messenger RNAs along the villus axis are zonated, yet often spatially discordant. Using spatial sorting with zonated surface markers, together with a Bayesian approach to infer protein translation and degradation rates from the combined spatial profiles, we find that, while many genes exhibit proteins zonated toward the villus tip, mRNA is zonated toward the villus bottom. Finally, we demonstrate that space-independent protein synthesis delays can explain many of the mRNA-protein discordances. Our work provides a proteomic spatial blueprint of the intestinal epithelium, highlighting the importance of protein measurements for inferring cell states in tissues that operate outside of steady state.

Single cell biology-a Keystone Symposia report.

Single cell biology has the potential to elucidate many critical biological processes and diseases, from development and regeneration to cancer. Single cell analyses are uncovering the molecular diversity of cells, revealing a clearer picture of the variation among and between different cell types. New techniques are beginning to unravel how differences in cell state-transcriptional, epigenetic, and other characteristics-can lead to different cell fates among genetically identical cells, which underlies complex processes such as embryonic development, drug resistance, response to injury, and cellular reprogramming. Single cell technologies also pose significant challenges relating to processing and analyzing vast amounts of data collected. To realize the potential of single cell technologies, new computational approaches are needed. On March 17-19, 2021, experts in single cell biology met virtually for the Keystone eSymposium "Single Cell Biology" to discuss advances both in single cell applications and technologies.

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.

A spatial vascular transcriptomic, proteomic, and phosphoproteomic atlas unveils an angiocrine Tie-Wnt signaling axis in the liver.

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.

Interleukin-6 produced by enteric neurons regulates the number and phenotype of microbe-responsive regulatory T cells in the gut.

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.

Bi-fated tendon-to-bone attachment cells are regulated by shared enhancers and KLF transcription factors.

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.

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.