Humanized mouse liver reveals endothelial control of essential hepatic metabolic functions.

Hepatocytes, the major metabolic hub of the body, execute functions that are human-specific, altered in human disease, and currently thought to be regulated through endocrine and cell-autonomous mechanisms. Here, we show that key metabolic functions of human hepatocytes are controlled by non-parenchymal cells (NPCs) in their microenvironment. We developed mice bearing human hepatic tissue composed of human hepatocytes and NPCs, including human immune, endothelial, and stellate cells. Humanized livers reproduce human liver architecture, perform vital human-specific metabolic/homeostatic processes, and model human pathologies, including fibrosis and non-alcoholic fatty liver disease (NAFLD). Leveraging species mismatch and lipidomics, we demonstrate that human NPCs control metabolic functions of human hepatocytes in a paracrine manner. Mechanistically, we uncover a species-specific interaction whereby WNT2 secreted by sinusoidal endothelial cells controls cholesterol uptake and bile acid conjugation in hepatocytes through receptor FZD5. These results reveal the essential microenvironmental regulation of hepatic metabolism and its human-specific aspects.

Enteric Nervous System-Derived IL-18 Orchestrates Mucosal Barrier Immunity.

Mucosal barrier immunity is essential for the maintenance of the commensal microflora and combating invasive bacterial infection. Although immune and epithelial cells are thought to be the canonical orchestrators of this complex equilibrium, here, we show that the enteric nervous system (ENS) plays an essential and non-redundant role in governing the antimicrobial protein (AMP) response. Using confocal microscopy and single-molecule fluorescence in situ mRNA hybridization (smFISH) studies, we observed that intestinal neurons produce the pleiotropic cytokine IL-18. Strikingly, deletion of IL-18 from the enteric neurons alone, but not immune or epithelial cells, rendered mice susceptible to invasive Salmonella typhimurium (S.t.) infection. Mechanistically, unbiased RNA sequencing and single-cell sequencing revealed that enteric neuronal IL-18 is specifically required for homeostatic goblet cell AMP production. Together, we show that neuron-derived IL-18 signaling controls tissue-wide intestinal immunity and has profound consequences on the mucosal barrier and invasive bacterial killing.

IL-10 constrains sphingolipid metabolism to limit inflammation.

Interleukin-10 (IL-10) is a key anti-inflammatory cytokine that can limit immune cell activation and cytokine production in innate immune cell types1. Loss of IL-10 signalling results in life-threatening inflammatory bowel disease in humans and mice-however, the exact mechanism by which IL-10 signalling subdues inflammation remains unclear2-5. Here we find that increased saturated very long chain (VLC) ceramides are critical for the heightened inflammatory gene expression that is a hallmark of IL-10 deficiency. Accordingly, genetic deletion of ceramide synthase 2 (encoded by Cers2), the enzyme responsible for VLC ceramide production, limited the exacerbated inflammatory gene expression programme associated with IL-10 deficiency both in vitro and in vivo. The accumulation of saturated VLC ceramides was regulated by a decrease in metabolic flux through the de novo mono-unsaturated fatty acid synthesis pathway. Restoring mono-unsaturated fatty acid availability to cells deficient in IL-10 signalling limited saturated VLC ceramide production and the associated inflammation. Mechanistically, we find that persistent inflammation mediated by VLC ceramides is largely dependent on sustained activity of REL, an immuno-modulatory transcription factor. Together, these data indicate that an IL-10-driven fatty acid desaturation programme rewires VLC ceramide accumulation and aberrant activation of REL. These studies support the idea that fatty acid homeostasis in innate immune cells serves as a key regulatory node to control pathologic inflammation and suggests that 'metabolic correction' of VLC homeostasis could be an important strategy to normalize dysregulated inflammation caused by the absence of IL-10.

Inflammasome activation in infected macrophages drives COVID-19 pathology.

Severe COVID-19 is characterized by persistent lung inflammation, inflammatory cytokine production, viral RNA and a sustained interferon (IFN) response, all of which are recapitulated and required for pathology in the SARS-CoV-2-infected MISTRG6-hACE2 humanized mouse model of COVID-19, which has a human immune system1-20. Blocking either viral replication with remdesivir21-23 or the downstream IFN-stimulated cascade with anti-IFNAR2 antibodies in vivo in the chronic stages of disease attenuates the overactive immune inflammatory response, especially inflammatory macrophages. Here we show that SARS-CoV-2 infection and replication in lung-resident human macrophages is a critical driver of disease. In response to infection mediated by CD16 and ACE2 receptors, human macrophages activate inflammasomes, release interleukin 1 (IL-1) and IL-18, and undergo pyroptosis, thereby contributing to the hyperinflammatory state of the lungs. Inflammasome activation and the accompanying inflammatory response are necessary for lung inflammation, as inhibition of the NLRP3 inflammasome pathway reverses chronic lung pathology. Notably, this blockade of inflammasome activation leads to the release of infectious virus by the infected macrophages. Thus, inflammasomes oppose host infection by SARS-CoV-2 through the production of inflammatory cytokines and suicide by pyroptosis to prevent a productive viral cycle.

Paracrine orchestration of intestinal tumorigenesis by a mesenchymal niche.

The initiation of an intestinal tumour is a probabilistic process that depends on the competition between mutant and normal epithelial stem cells in crypts1. Intestinal stem cells are closely associated with a diverse but poorly characterized network of mesenchymal cell types2,3. However, whether the physiological mesenchymal microenvironment of mutant stem cells affects tumour initiation remains unknown. Here we provide in vivo evidence that the mesenchymal niche controls tumour initiation in trans. By characterizing the heterogeneity of the intestinal mesenchyme using single-cell RNA-sequencing analysis, we identified a population of rare pericryptal Ptgs2-expressing fibroblasts that constitutively process arachidonic acid into highly labile prostaglandin E2 (PGE2). Specific ablation of Ptgs2 in fibroblasts was sufficient to prevent tumour initiation in two different models of sporadic, autochthonous tumorigenesis. Mechanistically, single-cell RNA-sequencing analyses of a mesenchymal niche model showed that fibroblast-derived PGE2 drives the expansion οf a population of Sca-1+ reserve-like stem cells. These express a strong regenerative/tumorigenic program, driven by the Hippo pathway effector Yap. In vivo, Yap is indispensable for Sca-1+ cell expansion and early tumour initiation and displays a nuclear localization in both mouse and human adenomas. Using organoid experiments, we identified a molecular mechanism whereby PGE2 promotes Yap dephosphorylation, nuclear translocation and transcriptional activity by signalling through the receptor Ptger4. Epithelial-specific ablation of Ptger4 misdirected the regenerative reprogramming of stem cells and prevented Sca-1+ cell expansion and sporadic tumour initiation in mutant mice, thereby demonstrating the robust paracrine control of tumour-initiating stem cells by PGE2-Ptger4. Analyses of patient-derived organoids established that PGE2-PTGER4 also regulates stem-cell function in humans. Our study demonstrates that initiation of colorectal cancer is orchestrated by the mesenchymal niche and reveals a mechanism by which rare pericryptal Ptgs2-expressing fibroblasts exert paracrine control over tumour-initiating stem cells via the druggable PGE2-Ptger4-Yap signalling axis.

Author Correction: Distinct modes of mitochondrial metabolism uncouple T cell differentiation and function.

An Amendment to this paper has been published and can be accessed via a link at the top of the paper.

Distinct modes of mitochondrial metabolism uncouple T cell differentiation and function.

Activated CD4 T cells proliferate rapidly and remodel epigenetically before exiting the cell cycle and engaging acquired effector functions. Metabolic reprogramming from the naive state is required throughout these phases of activation1. In CD4 T cells, T-cell-receptor ligation-along with co-stimulatory and cytokine signals-induces a glycolytic anabolic program that is required for biomass generation, rapid proliferation and effector function2. CD4 T cell differentiation (proliferation and epigenetic remodelling) and function are orchestrated coordinately by signal transduction and transcriptional remodelling. However, it remains unclear whether these processes are regulated independently of one another by cellular biochemical composition. Here we demonstrate that distinct modes of mitochondrial metabolism support differentiation and effector functions of mouse T helper 1 (TH1) cells by biochemically uncoupling these two processes. We find that the tricarboxylic acid cycle is required for the terminal effector function of TH1 cells through succinate dehydrogenase (complex II), but that the activity of succinate dehydrogenase suppresses TH1 cell proliferation and histone acetylation. By contrast, we show that complex I of the electron transport chain, the malate-aspartate shuttle and mitochondrial citrate export are required to maintain synthesis of aspartate, which is necessary for the proliferation of T helper cells. Furthermore, we find that mitochondrial citrate export and the malate-aspartate shuttle promote histone acetylation, and specifically regulate the expression of genes involved in T cell activation. Combining genetic, pharmacological and metabolomics approaches, we demonstrate that the differentiation and terminal effector functions of T helper cells are biochemically uncoupled. These findings support a model in which the malate-aspartate shuttle, mitochondrial citrate export and complex I supply the substrates needed for proliferation and epigenetic remodelling early during T cell activation, whereas complex II consumes the substrates of these pathways, which antagonizes differentiation and enforces terminal effector function. Our data suggest that transcriptional programming acts together with a parallel biochemical network to enforce cell state.

RNA m6A demethylase ALKBH5 regulates the development of γδ T cells.

γδ T cells are an abundant T cell population at the mucosa and are important in providing immune surveillance as well as maintaining tissue homeostasis. However, despite γδ T cells' origin in the thymus, detailed mechanisms regulating γδ T cell development remain poorly understood. N6-methyladenosine (m6A) represents one of the most common posttranscriptional modifications of messenger RNA (mRNA) in mammalian cells, but whether it plays a role in γδ T cell biology is still unclear. Here, we show that depletion of the m6A demethylase ALKBH5 in lymphocytes specifically induces an expansion of γδ T cells, which confers enhanced protection against gastrointestinal Salmonella typhimurium infection. Mechanistically, loss of ALKBH5 favors the development of γδ T cell precursors by increasing the abundance of m6A RNA modification in thymocytes, which further reduces the expression of several target genes including Notch signaling components Jagged1 and Notch2. As a result, impairment of Jagged1/Notch2 signaling contributes to enhanced proliferation and differentiation of γδ T cell precursors, leading to an expanded mature γδ T cell repertoire. Taken together, our results indicate a checkpoint role of ALKBH5 and m6A modification in the regulation of γδ T cell early development.

An in vivo screen of noncoding loci reveals that Daedalus is a gatekeeper of an Ikaros-dependent checkpoint during haematopoiesis.

Haematopoiesis relies on tightly controlled gene expression patterns as development proceeds through a series of progenitors. While the regulation of hematopoietic development has been well studied, the role of noncoding elements in this critical process is a developing field. In particular, the discovery of new regulators of lymphopoiesis could have important implications for our understanding of the adaptive immune system and disease. Here we elucidate how a noncoding element is capable of regulating a broadly expressed transcription factor, Ikaros, in a lymphoid lineage-specific manner, such that it imbues Ikaros with the ability to specify the lymphoid lineage over alternate fates. Deletion of the Daedalus locus, which is proximal to Ikaros, led to a severe reduction in early lymphoid progenitors, exerting control over the earliest fate decisions during lymphoid lineage commitment. Daedalus locus deletion led to alterations in Ikaros isoform expression and a significant reduction in Ikaros protein. The Daedalus locus may function through direct DNA interaction as Hi-C analysis demonstrated an interaction between the two loci. Finally, we identify an Ikaros-regulated erythroid-lymphoid checkpoint that is governed by Daedalus in a lymphoid-lineage-specific manner. Daedalus appears to act as a gatekeeper of Ikaros's broad lineage-specifying functions, selectively stabilizing Ikaros activity in the lymphoid lineage and permitting diversion to the erythroid fate in its absence. These findings represent a key illustration of how a transcription factor with broad lineage expression must work in concert with noncoding elements to orchestrate hematopoietic lineage commitment.

A bedside to bench study of anti-PD-1, anti-CD40, and anti-CSF1R indicates that more is not necessarily better.

BACKGROUND: Stimulating inflammatory tumor associated macrophages can overcome resistance to PD-(L)1 blockade. We previously conducted a phase I trial of cabiralizumab (anti-CSF1R), sotigalimab (CD40-agonist) and nivolumab. Our current purpose was to study the activity and cellular effects of this three-drug regimen in anti-PD-1-resistant melanoma. METHODS: We employed a Simon's two-stage design and analyzed circulating immune cells from patients treated with this regimen for treatment-related changes. We assessed various dose levels of anti-CSF1R in murine melanoma models and studied the cellular and molecular effects. RESULTS: Thirteen patients were enrolled in the first stage. We observed one (7.7%) confirmed and one (7.7%) unconfirmed partial response, 5 patients had stable disease (38.5%) and 6 disease progression (42.6%). We elected not to proceed to the second stage. CyTOF analysis revealed a reduction in non-classical monocytes. Patients with prolonged stable disease or partial response who remained on study for longer had increased markers of antigen presentation after treatment compared to patients whose disease progressed rapidly. In a murine model, higher anti-CSF1R doses resulted in increased tumor growth and worse survival. Using single-cell RNA-sequencing, we identified a suppressive monocyte/macrophage population in murine tumors exposed to higher doses. CONCLUSIONS: Higher anti-CSF1R doses are inferior to lower doses in a preclinical model, inducing a suppressive macrophage population, and potentially explaining the disappointing results observed in patients. While it is impossible to directly infer human doses from murine studies, careful intra-species evaluation can provide important insight. Cabiralizumab dose optimization is necessary for this patient population with limited treatment options. TRIAL REGISTRATION: ClinicalTrials.gov Identifier: NCT03502330.

Longitudinal single-cell analysis of a patient receiving adoptive cell therapy reveals potential mechanisms of treatment failure.

Adoptive cell therapy (ACT) using tumor infiltrating lymphocytes (TIL) is being studied in multiple tumor types. However, little is known about clonal cell expansion in vitro and persistence of the ACT product in vivo. We performed single-cell RNA and T-Cell Receptor (TCR) sequencing on serial blood and tumor samples from a patient undergoing ACT, who did not respond. We found that clonal expansion varied during preparation of the ACT product, and only one expanded clone was preserved in the ACT product. The TCR of the preserved clone which persisted and remained activated for five months was previously reported as specific for cytomegalovirus and had upregulation of granzyme family genes and genes associated with effector functions (HLA-DQB1, LAT, HLA-DQA1, and KLRD1). Clones that contracted during TIL preparation had features of exhaustion and apoptosis. At disease progression, all previously detected clonotypes were detected. New clonotypes appearing in blood or tumor at disease progression were enriched for genes associated with cytotoxicity or stemness (FGFBP2, GNLY, GZMH, GZMK, IL7R, SELL and KLF2), and these might be harnessed for alternative cellular therapy or cytokine therapy. In-depth single-cell analyses of serial samples from additional ACT-treated patients is warranted, and viral- versus tumor-specificity should be carefully analyzed.

RNAex: an RNA secondary structure prediction server enhanced by high-throughput structure-probing data.

Several high-throughput technologies have been developed to probe RNA base pairs and loops at the transcriptome level in multiple species. However, to obtain the final RNA secondary structure, extensive effort and considerable expertise is required to statistically process the probing data and combine them with free energy models. Therefore, we developed an RNA secondary structure prediction server that is enhanced by experimental data (RNAex). RNAex is a web interface that enables non-specialists to easily access cutting-edge structure-probing data and predict RNA secondary structures enhanced by in vivo and in vitro data. RNAex annotates the RNA editing, RNA modification and SNP sites on the predicted structures. It provides four structure-folding methods, restrained MaxExpect, SeqFold, RNAstructure (Fold) and RNAfold that can be selected by the user. The performance of these four folding methods has been verified by previous publications on known structures. We re-mapped the raw sequencing data of the probing experiments to the whole genome for each species. RNAex thus enables users to predict secondary structures for both known and novel RNA transcripts in human, mouse, yeast and Arabidopsis The RNAex web server is available at http://RNAex.ncrnalab.org/.

Generation and optimization of off-the-shelf immunotherapeutics targeting TCR-Vß2+ T cell malignancy.

Current treatments for T cell malignancies encounter issues of disease relapse and off-target toxicity. Using T cell receptor (TCR)Vß2 as a model, here we demonstrate the rapid generation of an off-the-shelf allogeneic chimeric antigen receptor (CAR)-T platform targeting the clone-specific TCR Vß chain for malignant T cell killing while limiting normal cell destruction. Healthy donor T cells undergo CRISPR-induced TRAC, B2M and CIITA knockout to eliminate T cell-dependent graft-versus-host and host-versus-graft reactivity. Second generation 4-1BB/CD3zeta CAR containing high affinity humanized anti-Vß scFv is expressed efficiently on donor T cells via both lentivirus and adeno-associated virus transduction with limited detectable pre-existing immunoreactivity. Our optimized CAR-T cells demonstrate specific and persistent killing of Vß2+ Jurkat cells and Vß2+ patient derived malignant T cells, in vitro and in vivo, without affecting normal T cells. In parallel, we generate humanized anti-Vß2 antibody with enhanced antibody-dependent cellular cytotoxicity (ADCC) by Fc-engineering for NK cell ADCC therapy.

Gene trajectory inference for single-cell data by optimal transport metrics.

Single-cell RNA sequencing has been widely used to investigate cell state transitions and gene dynamics of biological processes. Current strategies to infer the sequential dynamics of genes in a process typically rely on constructing cell pseudotime through cell trajectory inference. However, the presence of concurrent gene processes in the same group of cells and technical noise can obscure the true progression of the processes studied. To address this challenge, we present GeneTrajectory, an approach that identifies trajectories of genes rather than trajectories of cells. Specifically, optimal transport distances are calculated between gene distributions across the cell-cell graph to extract gene programs and define their gene pseudotemporal order. Here we demonstrate that GeneTrajectory accurately extracts progressive gene dynamics in myeloid lineage maturation. Moreover, we show that GeneTrajectory deconvolves key gene programs underlying mouse skin hair follicle dermal condensate differentiation that could not be resolved by cell trajectory approaches. GeneTrajectory facilitates the discovery of gene programs that control the changes and activities of biological processes.

Wnt5 controls splenic myelopoiesis and neutrophil functional ambivalency during DSS-induced colitis.

Neutrophils are important innate immune cells with plasticity, heterogenicity, and functional ambivalency. While bone marrow is often regarded as the primary source of neutrophil production, the roles of extramedullary production in regulating neutrophil plasticity and heterogenicity in autoimmune diseases remain poorly understood. Here, we report that the lack of wingless-type MMTV integration site family member 5 (WNT5) unleashes anti-inflammatory protection against colitis in mice, accompanied by reduced colonic CD8+ T cell activation and enhanced splenic extramedullary myelopoiesis. In addition, colitis upregulates WNT5 expression in splenic stromal cells. The ablation of WNT5 leads to increased splenic production of hematopoietic niche factors, as well as elevated numbers of splenic neutrophils with heightened CD8+ T cell suppressive capability, in part due to elevated CD101 expression and attenuated pro-inflammatory activities. Thus, our study reveals a mechanism by which neutrophil plasticity and heterogenicity are regulated in colitis through WNT5 and highlights the role of splenic neutrophil production in shaping inflammatory outcomes.

IL-4 Licenses B Cell Activation Through Cholesterol Synthesis.

Lymphocyte activation involves a transition from quiescence and associated catabolic metabolism to a metabolic state with noted similarities to cancer cells such as heavy reliance on aerobic glycolysis for energy demands and increased nutrient requirements for biomass accumulation and cell division 1-3 . Following antigen receptor ligation, lymphocytes require spatiotemporally distinct "second signals". These include costimulatory receptor or cytokine signaling, which engage discrete programs that often involve remodeling of organelles and increased nutrient uptake or synthesis to meet changing biochemical demands 4-6 . One such signaling molecule, IL-4, is a highly pleiotropic cytokine that was first identified as a B cell co-mitogen over 30 years ago 7 . However, how IL-4 signaling mechanistically supports B cell proliferation is incompletely understood. Here, using single cell RNA sequencing we find that the cholesterol biosynthetic program is transcriptionally upregulated following IL-4 signaling during the early B cell response to influenza virus infection, and is required for B cell activation in vivo . By limiting lipid availability in vitro , we determine cholesterol to be essential for B cells to expand their endoplasmic reticulum, progress through cell cycle, and proliferate. In sum, we demonstrate that the well-known ability of IL-4 to act as a B cell growth factor is through a previously unknown rewiring of specific lipid anabolic programs, relieving sensitivity of cells to environmental nutrient availability.

LMD: Multiscale Marker Identification in Single-cell RNA-seq Data.

Accurate cell marker identification in single-cell RNA-seq data is crucial for understanding cellular diversity and function. An ideal marker is highly specific in identifying cells that are similar in terms of function and state. Current marker identification methods, commonly based on clustering and differential expression, capture general cell-type markers but often miss markers for subtypes or functional cell subsets, with their performance largely dependent on clustering quality. Moreover, cluster-independent approaches tend to favor genes that lack the specificity required to characterize regions within the transcriptomic space at multiple scales. Here we introduce Localized Marker Detector (LMD), a novel tool to identify "localized genes" - genes with expression profiles specific to certain groups of highly similar cells - thereby characterizing cellular diversity in a multi-resolution and fine-grained manner. LMD's strategy involves building a cell-cell affinity graph, diffusing the gene expression value across the cell graph, and assigning a score to each gene based on its diffusion dynamics. We show that LMD exhibits superior accuracy in recovering known cell-type markers in the Tabula Muris bone marrow dataset relative to other methods for marker identification. Notably, markers favored by LMD exhibit localized expression, whereas markers prioritized by other clustering-free algorithms are often dispersed in the transcriptomic space. We further group the markers suggested by LMD into functional gene modules to improve the separation of cell types and subtypes in a more fine-grained manner. These modules also identify other sources of variation, such as cell cycle status. In conclusion, LMD is a novel algorithm that can identify fine-grained markers for cell subtypes or functional states without relying on clustering or differential expression analysis. LMD exploits the complex interactions among cells and reveals cellular diversity at high resolution.

Evaluation of zero counts to better understand the discrepancies between bulk and single-cell RNA-Seq platforms.

Recent advances in sample preparation and sequencing technology have made it possible to profile the transcriptomes of individual cells using single-cell RNA sequencing (scRNA-Seq). Compared to bulk RNA-Seq data, single-cell data often contain a higher percentage of zero reads, mainly due to lower sequencing depth per cell, which affects mostly measurements of low-expression genes. However, discrepancies between platforms are observed regardless of expression level. Using four paired datasets with multiple samples each, we investigated technical and biological factors that can contribute to this expression shift. Using two separate machine learning models we found that, in addition to expression level, RNA integrity, gene or UTR3 length, and the number of transcripts potentially also influence the occurrence of zeros. These findings could enable the development of novel analytical methods for cross-platform expression shift correction. We also identified genes and biological pathways in our diverse datasets that consistently showed differences when assessed at the single cell versus bulk level to assist in interpreting analysis across transcriptomic platforms. At the gene level, 25 genes (0.12%) were found in all datasets as discordant, but at the pathway level, 7 pathways (2.02%) showed shared enrichment in discordant genes.

IL-10 constrains sphingolipid metabolism via fatty acid desaturation to limit inflammation.

Unchecked chronic inflammation is the underlying cause of many diseases, ranging from inflammatory bowel disease to obesity and neurodegeneration. Given the deleterious nature of unregulated inflammation, it is not surprising that cells have acquired a diverse arsenal of tactics to limit inflammation. IL-10 is a key anti-inflammatory cytokine that can limit immune cell activation and cytokine production in innate immune cell types; however, the exact mechanism by which IL-10 signaling subdues inflammation remains unclear. Here, we find that IL-10 signaling constrains sphingolipid metabolism. Specifically, we find increased saturated very long chain (VLC) ceramides are critical for the heightened inflammatory gene expression that is a hallmark of IL-10-deficient macrophages. Genetic deletion of CerS2, the enzyme responsible for VLC ceramide production, limited exacerbated inflammatory gene expression associated with IL-10 deficiency both in vitro and in vivo , indicating that "metabolic correction" is able to reduce inflammation in the absence of IL-10. Surprisingly, accumulation of saturated VLC ceramides was regulated by flux through the de novo mono-unsaturated fatty acid (MUFA) synthesis pathway, where addition of exogenous MUFAs could limit both saturated VLC ceramide production and inflammatory gene expression in the absence of IL-10 signaling. Together, these studies mechanistically define how IL-10 signaling manipulates fatty acid metabolism as part of its molecular anti-inflammatory strategy and could lead to novel and inexpensive approaches to regulate aberrant inflammation.