LGR5 expressing skin fibroblasts define a major cellular hub perturbed in scleroderma.

Systemic sclerosis (scleroderma, SSc) is an incurable autoimmune disease with high morbidity and mortality rates. Here, we conducted a population-scale single-cell genomic analysis of skin and blood samples of 56 healthy controls and 97 SSc patients at different stages of the disease. We found immune compartment dysfunction only in a specific subtype of diffuse SSc patients but global dysregulation of the stromal compartment, particularly in a previously undefined subset of LGR5+-scleroderma-associated fibroblasts (ScAFs). ScAFs are perturbed morphologically and molecularly in SSc patients. Single-cell multiome profiling of stromal cells revealed ScAF-specific markers, pathways, regulatory elements, and transcription factors underlining disease development. Systematic analysis of these molecular features with clinical metadata associates specific ScAF targets with disease pathogenesis and SSc clinical traits. Our high-resolution atlas of the sclerodermatous skin spectrum will enable a paradigm shift in the understanding of SSc disease and facilitate the development of biomarkers and therapeutic strategies.

Spatially organized multicellular immune hubs in human colorectal cancer.

Immune responses to cancer are highly variable, with mismatch repair-deficient (MMRd) tumors exhibiting more anti-tumor immunity than mismatch repair-proficient (MMRp) tumors. To understand the rules governing these varied responses, we transcriptionally profiled 371,223 cells from colorectal tumors and adjacent normal tissues of 28 MMRp and 34 MMRd individuals. Analysis of 88 cell subsets and their 204 associated gene expression programs revealed extensive transcriptional and spatial remodeling across tumors. To discover hubs of interacting malignant and immune cells, we identified expression programs in different cell types that co-varied across tumors from affected individuals and used spatial profiling to localize coordinated programs. We discovered a myeloid cell-attracting hub at the tumor-luminal interface associated with tissue damage and an MMRd-enriched immune hub within the tumor, with activated T cells together with malignant and myeloid cells expressing T cell-attracting chemokines. By identifying interacting cellular programs, we reveal the logic underlying spatially organized immune-malignant cell networks.

Ketone Body Signaling Mediates Intestinal Stem Cell Homeostasis and Adaptation to Diet.

Little is known about how metabolites couple tissue-specific stem cell function with physiology. Here we show that, in the mammalian small intestine, the expression of Hmgcs2 (3-hydroxy-3-methylglutaryl-CoA synthetase 2), the gene encoding the rate-limiting enzyme in the production of ketone bodies, including beta-hydroxybutyrate (ßOHB), distinguishes self-renewing Lgr5+ stem cells (ISCs) from differentiated cell types. Hmgcs2 loss depletes ßOHB levels in Lgr5+ ISCs and skews their differentiation toward secretory cell fates, which can be rescued by exogenous ßOHB and class I histone deacetylase (HDAC) inhibitor treatment. Mechanistically, ßOHB acts by inhibiting HDACs to reinforce Notch signaling, instructing ISC self-renewal and lineage decisions. Notably, although a high-fat ketogenic diet elevates ISC function and post-injury regeneration through ßOHB-mediated Notch signaling, a glucose-supplemented diet has the opposite effects. These findings reveal how control of ßOHB-activated signaling in ISCs by diet helps to fine-tune stem cell adaptation in homeostasis and injury.

Intra- and Inter-cellular Rewiring of the Human Colon during Ulcerative Colitis.

Genome-wide association studies (GWAS) have revealed risk alleles for ulcerative colitis (UC). To understand their cell type specificities and pathways of action, we generate an atlas of 366,650 cells from the colon mucosa of 18 UC patients and 12 healthy individuals, revealing 51 epithelial, stromal, and immune cell subsets, including BEST4+ enterocytes, microfold-like cells, and IL13RA2+IL11+ inflammatory fibroblasts, which we associate with resistance to anti-TNF treatment. Inflammatory fibroblasts, inflammatory monocytes, microfold-like cells, and T cells that co-express CD8 and IL-17 expand with disease, forming intercellular interaction hubs. Many UC risk genes are cell type specific and co-regulated within relatively few gene modules, suggesting convergence onto limited sets of cell types and pathways. Using this observation, we nominate and infer functions for specific risk genes across GWAS loci. Our work provides a framework for interrogating complex human diseases and mapping risk variants to cell types and pathways.

T Helper Cell Cytokines Modulate Intestinal Stem Cell Renewal and Differentiation.

In the small intestine, a niche of accessory cell types supports the generation of mature epithelial cell types from intestinal stem cells (ISCs). It is unclear, however, if and how immune cells in the niche affect ISC fate or the balance between self-renewal and differentiation. Here, we use single-cell RNA sequencing (scRNA-seq) to identify MHC class II (MHCII) machinery enrichment in two subsets of Lgr5+ ISCs. We show that MHCII+ Lgr5+ ISCs are non-conventional antigen-presenting cells in co-cultures with CD4+ T helper (Th) cells. Stimulation of intestinal organoids with key Th cytokines affects Lgr5+ ISC renewal and differentiation in opposing ways: pro-inflammatory signals promote differentiation, while regulatory cells and cytokines reduce it. In vivo genetic perturbation of Th cells or MHCII expression on Lgr5+ ISCs impacts epithelial cell differentiation and IEC fate during infection. These interactions between Th cells and Lgr5+ ISCs, thus, orchestrate tissue-wide responses to external signals.

Bacterial infection disrupts established germinal center reactions through monocyte recruitment and impaired metabolic adaptation.

Consecutive exposures to different pathogens are highly prevalent and often alter the host immune response. However, it remains unknown how a secondary bacterial infection affects an ongoing adaptive immune response elicited against primary invading pathogens. We demonstrated that recruitment of Sca-1+ monocytes into lymphoid organs during Salmonella Typhimurium (STm) infection disrupted pre-existing germinal center (GC) reactions. GC responses induced by influenza, plasmodium, or commensals deteriorated following STm infection. GC disruption was independent of the direct bacterial interactions with B cells and instead was induced through recruitment of CCR2-dependent Sca-1+ monocytes into the lymphoid organs. GC collapse was associated with impaired cellular respiration and was dependent on TNFα and IFNγ, the latter of which was essential for Sca-1+ monocyte differentiation. Monocyte recruitment and GC disruption also occurred during LPS-supplemented vaccination and Listeria monocytogenes infection. Thus, systemic activation of the innate immune response upon severe bacterial infection is induced at the expense of antibody-mediated immunity.

A revised airway epithelial hierarchy includes CFTR-expressing ionocytes.

The airways of the lung are the primary sites of disease in asthma and cystic fibrosis. Here we study the cellular composition and hierarchy of the mouse tracheal epithelium by single-cell RNA-sequencing (scRNA-seq) and in vivo lineage tracing. We identify a rare cell type, the Foxi1+ pulmonary ionocyte; functional variations in club cells based on their location; a distinct cell type in high turnover squamous epithelial structures that we term 'hillocks'; and disease-relevant subsets of tuft and goblet cells. We developed 'pulse-seq', combining scRNA-seq and lineage tracing, to show that tuft, neuroendocrine and ionocyte cells are continually and directly replenished by basal progenitor cells. Ionocytes are the major source of transcripts of the cystic fibrosis transmembrane conductance regulator in both mouse (Cftr) and human (CFTR). Knockout of Foxi1 in mouse ionocytes causes loss of Cftr expression and disrupts airway fluid and mucus physiology, phenotypes that are characteristic of cystic fibrosis. By associating cell-type-specific expression programs with key disease genes, we establish a new cellular narrative for airways disease.

Epithelial microRNAs regulate gut mucosal immunity via epithelium-T cell crosstalk.

Colonic homeostasis entails epithelium-lymphocyte cooperation, yet many participants in this process are unknown. We show here that epithelial microRNAs mediate the mucosa-immune system crosstalk necessary for mounting protective T helper type 2 (T(H)2) responses. Abolishing the induction of microRNA by gut-specific deletion of Dicer1 (Dicer1(Δgut)), which encodes an enzyme involved in microRNA biogenesis, deprived goblet cells of RELMß, a key T(H)2 antiparasitic cytokine; this predisposed the host to parasite infection. Infection of Dicer1(Δgut) mice with helminths favored a futile T(H)1 response with hallmarks of inflammatory bowel disease. Interleukin 13 (IL-13) induced the microRNA miR-375, which regulates the expression of TSLP, a T(H)2-facilitating epithelial cytokine; this indicated a T(H)2-amplification loop. We found that miR-375 was required for RELMß expression in vivo; miR-375-deficient mice had significantly less intestinal RELMß, which possibly explains the greater susceptibility of Dicer1(Δgut) mice to parasites. Our findings indicate that epithelial microRNAs are key regulators of gut homeostasis and mucosal immunity.

A single-cell survey of the small intestinal epithelium.

Intestinal epithelial cells absorb nutrients, respond to microbes, function as a barrier and help to coordinate immune responses. Here we report profiling of 53,193 individual epithelial cells from the small intestine and organoids of mice, which enabled the identification and characterization of previously unknown subtypes of intestinal epithelial cell and their gene signatures. We found unexpected diversity in hormone-secreting enteroendocrine cells and constructed the taxonomy of newly identified subtypes, and distinguished between two subtypes of tuft cell, one of which expresses the epithelial cytokine Tslp and the pan-immune marker CD45, which was not previously associated with non-haematopoietic cells. We also characterized the ways in which cell-intrinsic states and the proportions of different cell types respond to bacterial and helminth infections: Salmonella infection caused an increase in the abundance of Paneth cells and enterocytes, and broad activation of an antimicrobial program; Heligmosomoides polygyrus caused an increase in the abundance of goblet and tuft cells. Our survey highlights previously unidentified markers and programs, associates sensory molecules with cell types, and uncovers principles of gut homeostasis and response to pathogens.

Distinct Tissue-Specific Roles for the Disease-Associated Autophagy Genes ATG16L2 and ATG16L1.

The clear role of autophagy in human inflammatory diseases such as Crohn disease was first identified by genome-wide association studies and subsequently dissected in multiple mechanistic studies. ATG16L1 has been particularly well studied in knockout and hypomorph settings as well as models recapitulating the Crohn disease-associated T300A polymorphism. Interestingly, ATG16L1 has a single homolog, ATG16L2, which is independently implicated in diseases, including Crohn disease and systemic lupus erythematosus. However, the contribution of ATG16L2 to canonical autophagy pathways and other cellular functions is poorly understood. To better understand its role, we generated and analyzed the first, to our knowledge, ATG16L2 knockout mouse. Our results show that ATG16L1 and ATG16L2 contribute very distinctly to autophagy and cellular ontogeny in myeloid, lymphoid, and epithelial lineages. Dysregulation of any of these lineages could contribute to complex diseases like Crohn disease and systemic lupus erythematosus, highlighting the value of examining cell-specific effects. We also identify a novel genetic interaction between ATG16L2 and epithelial ATG16L1. These findings are discussed in the context of how these genes may contribute distinctly to human disease.

Human microRNAs regulate stress-induced immune responses mediated by the receptor NKG2D.

MICA and MICB are stress-induced ligands recognized by the activating receptor NKG2D. A microRNA encoded by human cytomegalovirus downregulates MICB expression by targeting a specific site in the MICB 3' untranslated region. As this site is conserved among different MICB alleles and a similar site exists in the MICA 3' untranslated region, we speculated that these sites are targeted by cellular microRNAs. Here we identified microRNAs that bound to these MICA and MICB 3' untranslated region sequences and obtained data suggesting that these microRNAs maintain expression of MICA and MICB protein under a certain threshold and facilitate acute upregulation of MICA and MICB during cellular stress. These microRNAs were overexpressed in various tumors and we demonstrate here that they aided tumor avoidance of immune recognition.

Forefront: MiR-34a-Knockout Mice with Wild Type Hematopoietic Cells, Retain Persistent Fibrosis Following Lung Injury.

MicroRNAs (miRs) are known to limit gene expression at the post-transcriptional level and have important roles in the pathogenesis of various conditions, including acute lung injury (ALI) and fibrotic diseases such as idiopathic pulmonary fibrosis (IPF). In this study, we found increased levels of miR-34 at times of fibrosis resolution following injury, in myofibroblasts from Bleomycin-treated mouse lungs, which correlates with susceptibility to cell death induced by immune cells. On the contrary, a substantial downregulation of miR-34 was detected at stages of evolution, when fibroblasts resist cell death. Concomitantly, we found an inverse correlation between miR-34 levels with that of the survival molecule FLICE-like inhibitory protein (FLIP) in lung myofibroblasts from humans with IPF and the experimental model. Forced upregulation of miR-34 with miR-34 mimic in human IPF fibrotic-lung myofibroblasts led to decreased cell survival through downregulation of FLIP. Using chimeric miR-34 knock-out (KO)-C57BL/6 mice with miR34KO myofibroblasts but wild-type (WT) hematopoietic cells, we found, in contrast to WT mice, increased and persistent FLIP levels with a more severe fibrosis and with no signs of resolution as detected in pathology and collagen accumulation. Moreover, a mimic of miR-34a decreased FLIP expression and susceptibility to cell death was regained in miR-34KO fibroblasts. Through this study, we show for the first time an inverse correlation between miR-34a and FLIP expression in myofibroblasts, which affects survival, and accumulation in lung fibrosis. Reprogramming fibrotic-lung myofibroblasts to regain susceptibility to cell-death by specifically increasing their miR34a and downregulating FLIP, may be a useful strategy, enabling tissue regeneration following lung injury.

Longitudinal single-cell data informs deterministic modelling of inflammatory bowel disease.

Single-cell-based methods such as flow cytometry or single-cell mRNA sequencing (scRNA-seq) allow deep molecular and cellular profiling of immunological processes. Despite their high throughput, however, these measurements represent only a snapshot in time. Here, we explore how longitudinal single-cell-based datasets can be used for deterministic ordinary differential equation (ODE)-based modelling to mechanistically describe immune dynamics. We derived longitudinal changes in cell numbers of colonic cell types during inflammatory bowel disease (IBD) from flow cytometry and scRNA-seq data of murine colitis using ODE-based models. Our mathematical model generalised well across different protocols and experimental techniques, and we hypothesised that the estimated model parameters reflect biological processes. We validated this prediction of cellular turnover rates with KI-67 staining and with gene expression information from the scRNA-seq data not used for model fitting. Finally, we tested the translational relevance of the mathematical model by deconvolution of longitudinal bulk mRNA-sequencing data from a cohort of human IBD patients treated with olamkicept. We found that neutrophil depletion may contribute to IBD patients entering remission. The predictive power of IBD deterministic modelling highlights its potential to advance our understanding of immune dynamics in health and disease.

Dendritic cell ICAM-1 strengthens synapses with CD8 T cells but is not required for their early differentiation.

Lymphocyte priming in lymph nodes (LNs) was postulated to depend on the formation of stable T cell receptor (TCR)-specific immune synapses (ISs) with antigen (Ag)-presenting dendritic cells (DCs). The high-affinity LFA-1 ligand ICAM-1 was implicated in different ISs studied in vitro. We dissect the in vivo roles of endogenous DC ICAM-1 in Ag-stimulated T cell proliferation and differentiation and find that under type 1 polarizing conditions in vaccinated or vaccinia virus-infected skin-draining LNs, Ag-presenting DCs engage in ICAM-1-dependent stable conjugates with a subset of Ag-specific CD8 blasts. Nevertheless, in the absence of these conjugates, CD8 lymphocyte proliferation and differentiation into functional cytotoxic T cells (CTLs) and skin homing effector lymphocytes takes place normally. Our results suggest that although CD8 T cell blasts engage in tight ICAM-1-dependent DC-T ISs, firm ISs are dispensable for TCR-triggered proliferation and differentiation into productive effector lymphocytes.

Leveraging polygenic enrichments of gene features to predict genes underlying complex traits and diseases.

Genome-wide association studies (GWASs) are a valuable tool for understanding the biology of complex human traits and diseases, but associated variants rarely point directly to causal genes. In the present study, we introduce a new method, polygenic priority score (PoPS), that learns trait-relevant gene features, such as cell-type-specific expression, to prioritize genes at GWAS loci. Using a large evaluation set of genes with fine-mapped coding variants, we show that PoPS and the closest gene individually outperform other gene prioritization methods, but observe the best overall performance by combining PoPS with orthogonal methods. Using this combined approach, we prioritize 10,642 unique gene-trait pairs across 113 complex traits and diseases with high precision, finding not only well-established gene-trait relationships but nominating new genes at unresolved loci, such as LGR4 for estimated glomerular filtration rate and CCR7 for deep vein thrombosis. Overall, we demonstrate that PoPS provides a powerful addition to the gene prioritization toolbox.

γδ T cells regulate the intestinal response to nutrient sensing.

The intestine is a site of direct encounter with the external environment and must consequently balance barrier defense with nutrient uptake. To investigate how nutrient uptake is regulated in the small intestine, we tested the effect of diets with different macronutrient compositions on epithelial gene expression. We found that enzymes and transporters required for carbohydrate digestion and absorption were regulated by carbohydrate availability. The "on-demand" induction of this machinery required γδ T cells, which regulated this program through the suppression of interleukin-22 production by type 3 innate lymphoid cells. Nutrient availability altered the tissue localization and transcriptome of γδ T cells. Additionally, transcriptional responses to diet involved cellular remodeling of the epithelial compartment. Thus, this work identifies a role for γδ T cells in nutrient sensing.

Dietary suppression of MHC class II expression in intestinal epithelial cells enhances intestinal tumorigenesis.

Little is known about how interactions of diet, intestinal stem cells (ISCs), and immune cells affect early-stage intestinal tumorigenesis. We show that a high-fat diet (HFD) reduces the expression of the major histocompatibility complex class II (MHC class II) genes in intestinal epithelial cells, including ISCs. This decline in epithelial MHC class II expression in a HFD correlates with reduced intestinal microbiome diversity. Microbial community transfer experiments suggest that epithelial MHC class II expression is regulated by intestinal flora. Mechanistically, pattern recognition receptor (PRR) and interferon-gamma (IFNγ) signaling regulates epithelial MHC class II expression. MHC class II-negative (MHC-II-) ISCs exhibit greater tumor-initiating capacity than their MHC class II-positive (MHC-II+) counterparts upon loss of the tumor suppressor Apc coupled with a HFD, suggesting a role for epithelial MHC class II-mediated immune surveillance in suppressing tumorigenesis. ISC-specific genetic ablation of MHC class II increases tumor burden cell autonomously. Thus, HFD perturbs a microbiome-stem cell-immune cell interaction that contributes to tumor initiation in the intestine.

The Taste Receptor TAS1R3 Regulates Small Intestinal Tuft Cell Homeostasis.

Tuft cells are an epithelial cell type critical for initiating type 2 immune responses to parasites and protozoa in the small intestine. To respond to these stimuli, intestinal tuft cells use taste chemosensory signaling pathways, but the role of taste receptors in type 2 immunity is poorly understood. In this study, we show that the taste receptor TAS1R3, which detects sweet and umami in the tongue, also regulates tuft cell responses in the distal small intestine. BALB/c mice, which have an inactive form of TAS1R3, as well as Tas1r3-deficient C57BL6/J mice both have severely impaired responses to tuft cell-inducing signals in the ileum, including the protozoa Tritrichomonas muris and succinate. In contrast, TAS1R3 is not required to mount an immune response to the helminth Heligmosomoides polygyrus, which infects the proximal small intestine. Examination of uninfected Tas1r3-/- mice revealed a modest reduction in the number of tuft cells in the proximal small intestine but a severe decrease in the distal small intestine at homeostasis. Together, these results suggest that TAS1R3 influences intestinal immunity by shaping the epithelial cell landscape at steady-state.

Regenerative potential of prostate luminal cells revealed by single-cell analysis.

Androgen deprivation is the cornerstone of prostate cancer treatment. It results in involution of the normal gland to ~90% of its original size because of the loss of luminal cells. The prostate regenerates when androgen is restored, a process postulated to involve stem cells. Using single-cell RNA sequencing, we identified a rare luminal population in the mouse prostate that expresses stemlike genes (Sca1 + and Psca +) and a large population of differentiated cells (Nkx3.1 +, Pbsn +). In organoids and in mice, both populations contribute equally to prostate regeneration, partly through androgen-driven expression of growth factors (Nrg2, Rspo3) by mesenchymal cells acting in a paracrine fashion on luminal cells. Analysis of human prostate tissue revealed similar differentiated and stemlike luminal subpopulations that likewise acquire enhanced regenerative potential after androgen ablation. We propose that prostate regeneration is driven by nearly all persisting luminal cells, not just by rare stem cells.

miR-103 inhibits proliferation and sensitizes hemopoietic tumor cells for glucocorticoid-induced apoptosis.

Glucocorticoid (GC) hormones are an important ingredient of leukemia therapy since they are potent inducers of lymphoid cell apoptosis. However, the development of GC resistance remains an obstacle in GC-based treatment. In the present investigation we found that miR-103 is upregulated in GC-sensitive leukemia cells treated by the hormone. Transfection of GC resistant cells with miR-103 sensitized them to GC induced apoptosis (GCIA), while miR-103 sponging of GC sensitive cells rendered them partially resistant. miR-103 reduced the expression of cyclin dependent kinase (CDK2) and its cyclin E1 target, thereby leading to inhibition of cellular proliferation. miR-103 is encoded within the fifth intron of PANK3 gene. We demonstrate that the GC receptor (GR) upregulates miR-103 by direct interaction with GC response element (GRE) in the PANK3 enhancer. Consequently, miR-103 targets the c-Myc activators c-Myb and DVL1, thereby reducing c-Myc expression. Since c-Myc is a transcription factor of the miR-17~92a poly-cistron, all six miRNAs of the latter are also downregulated. Of these, miR-18a and miR-20a are involved in GCIA, as they target GR and BIM, respectively. Consequently, GR and BIM expression are elevated, thus advancing GCIA. Altogether, this study highlights miR-103 as a useful prognostic biomarker and drug for leukemia management in the future.