Transcriptional profiling of retinal astrocytes identifies a specific marker and points to functional specialization.

Astrocyte heterogeneity is an increasingly prominent research topic, and studies in the brain have demonstrated substantial variation in astrocyte form and function, both between and within regions. In contrast, retinal astrocytes are not well understood and remain incompletely characterized. Along with optic nerve astrocytes, they are responsible for supporting retinal ganglion cell axons and an improved understanding of their role is required. We have used a combination of microdissection and Ribotag immunoprecipitation to isolate ribosome-associated mRNA from retinal astrocytes and investigate their transcriptome, which we also compared to astrocyte populations in the optic nerve. Astrocytes from these regions are transcriptionally distinct, and we identified retina-specific astrocyte genes and pathways. Moreover, although they share much of the "classical" gene expression patterns of astrocytes, we uncovered unexpected variation, including in genes related to core astrocyte functions. We additionally identified the transcription factor Pax8 as a highly specific marker of retinal astrocytes and demonstrated that these astrocytes populate not only the retinal surface, but also the prelaminar region at the optic nerve head. These findings are likely to contribute to a revised understanding of the role of astrocytes in the retina.

Human regulatory memory B cells defined by expression of TIM-1 and TIGIT are dysfunctional in multiple sclerosis.

Background: Regulatory B cells (Bregs) play a pivotal role in suppressing immune responses, yet there is still a lack of cell surface markers that can rigorously identify them. In mouse models for multiple sclerosis (MS), TIM-1 or TIGIT expression on B cells is required for maintaining self-tolerance and regulating autoimmunity to the central nervous system. Here we investigated the activities of human memory B cells that differentially express TIM-1 and TIGIT to determine their potential regulatory function in healthy donors and patients with relapsing-remitting (RR) MS. Methods: FACS-sorted TIM-1+/-TIGIT+/- memory B (memB) cells co-cultured with allogenic CD4+ T cells were analyzed for proliferation and induction of inflammatory markers using flow cytometry and cytokine quantification, to determine Th1/Th17 cell differentiation. Transcriptional differences were assessed by SMARTSeq2 RNA sequencing analysis. Results: TIM-1-TIGIT- double negative (DN) memB cells strongly induce T cell proliferation and pro-inflammatory cytokine expression. The TIM-1+ memB cells enabled low levels of CD4+ T cell activation and gave rise to T cells that co-express IL-10 with IFNγ and IL-17A or FoxP3. T cells cultured with the TIM-1+TIGIT+ double positive (DP) memB cells exhibited reduced proliferation and IFNγ, IL-17A, TNFα, and GM-CSF expression, and exhibited strong regulation in Breg suppression assays. The functional activity suggests the DP memB cells are a bonafide Breg population. However, MS DP memB cells were less inhibitory than HC DP memB cells. A retrospective longitudinal study of anti-CD20 treated patients found that post-treatment DP memB cell frequency and absolute number were associated with response to therapy. Transcriptomic analyses indicated that the dysfunctional MS-derived DP memB/Breg population exhibited increased expression of genes associated with T cell activation and survival (CD80, ZNF10, PIK3CA), and had distinct gene expression compared to the TIGIT+ or TIM-1+ memB cells. Conclusion: These findings demonstrate that TIM-1/TIGIT expressing memory B cell subsets have distinct functionalities. Co-expression of TIM-1 and TIGIT defines a regulatory memory B cell subset that is functionally impaired in MS.

WNT signalling control by KDM5C during development affects cognition.

Although KDM5C is one of the most frequently mutated genes in X-linked intellectual disability1, the exact mechanisms that lead to cognitive impairment remain unknown. Here we use human patient-derived induced pluripotent stem cells and Kdm5c knockout mice to conduct cellular, transcriptomic, chromatin and behavioural studies. KDM5C is identified as a safeguard to ensure that neurodevelopment occurs at an appropriate timescale, the disruption of which leads to intellectual disability. Specifically, there is a developmental window during which KDM5C directly controls WNT output to regulate the timely transition of primary to intermediate progenitor cells and consequently neurogenesis. Treatment with WNT signalling modulators at specific times reveal that only a transient alteration of the canonical WNT signalling pathway is sufficient to rescue the transcriptomic and chromatin landscapes in patient-derived cells and to induce these changes in wild-type cells. Notably, WNT inhibition during this developmental period also rescues behavioural changes of Kdm5c knockout mice. Conversely, a single injection of WNT3A into the brains of wild-type embryonic mice cause anxiety and memory alterations. Our work identifies KDM5C as a crucial sentinel for neurodevelopment and sheds new light on KDM5C mutation-associated intellectual disability. The results also increase our general understanding of memory and anxiety formation, with the identification of WNT functioning in a transient nature to affect long-lasting cognitive function.

Sensory Input, Sex, and Function Shape Hypothalamic Cell Type Development.

Mammalian behavior and physiology undergo dramatic changes in early life. Young animals rely on conspecifics to meet their homeostatic needs, until weaning and puberty initiate nutritional independence and sex-specific social interactions, respectively. How neuronal populations regulating homeostatic functions and social behaviors develop and mature during these transitions remains unclear. We used paired transcriptomic and chromatin accessibility profiling to examine the developmental trajectories of neuronal populations in the hypothalamic preoptic region, where cell types with key roles in physiological and behavioral control have been identified1-6. These data reveal a remarkable diversity of developmental trajectories shaped by the sex of the animal, and the location and behavioral or physiological function of the corresponding cell types. We identify key stages of preoptic development, including the perinatal emergence of sex differences, postnatal maturation and subsequent refinement of signaling networks, and nonlinear transcriptional changes accelerating at the time of weaning and puberty. We assessed preoptic development in various sensory mutants and find a major role for vomeronasal sensing in the timing of preoptic cell type maturation. These results provide novel insights into the development of neurons controlling homeostatic functions and social behaviors and lay ground for examining the dynamics of these functions in early life.

The ORFIUS complex regulates ORC2 localization at replication origins.

High-grade serous ovarian cancer (HGSC) is a lethal malignancy with elevated replication stress (RS) levels and defective RS and RS-associated DNA damage responses. Here we demonstrate that the bromodomain-containing protein BRD1 is a RS suppressing protein that forms a replication origin regulatory complex with the histone acetyltransferase HBO1, the BRCA1 tumor suppressor, and BARD1, ORigin FIring Under Stress (ORFIUS). BRD1 and HBO1 promote eventual origin firing by supporting localization of the origin licensing protein ORC2 at origins. In the absence of BRD1 and/or HBO1, both origin firing and nuclei with ORC2 foci are reduced. BRCA1 regulates BRD1, HBO1, and ORC2 localization at replication origins. In the absence of BRCA1, both origin firing and nuclei with BRD1, HBO1, and ORC2 foci are increased. In normal and non-HGSC ovarian cancer cells, the ORFIUS complex responds to ATR and CDC7 origin regulatory signaling and disengages from origins during RS. In BRCA1-mutant and sporadic HGSC cells, BRD1, HBO1, and ORC2 remain associated with replication origins, and unresponsive to RS, DNA damage, or origin regulatory kinase inhibition. ORFIUS complex dysregulation may promote HGSC cell survival by allowing for upregulated origin firing and cell cycle progression despite accumulating DNA damage, and may be a RS target.

Age-associated H3K9me2 loss alters the regenerative equilibrium between murine lung alveolar and bronchiolar progenitors.

The lung contains multiple progenitor cell types, but how their responses are choreographed during injury repair and whether this changes with age is poorly understood. We report that histone H3 lysine 9 di-methylation (H3K9me2), mediated by the methyltransferase G9a, regulates the dynamics of distal lung epithelial progenitor cells and that this regulation deteriorates with age. In aged mouse lungs, H3K9me2 loss coincided with fewer alveolar type 2 (AT2) cell progenitors and reduced alveolar regeneration but increased the frequency and activity of multipotent bronchioalveolar stem cells (BASCs) and bronchiolar progenitor club cells. H3K9me2 depletion in young mice decreased AT2 progenitor activity and impaired alveolar injury repair. Conversely, H3K9me2 depletion increased chromatin accessibility of bronchiolar cell genes, increased BASC frequency, and accelerated bronchiolar cell injury repair. These findings indicate that during aging, the epigenetic regulation that coordinates lung progenitor cells' regenerative responses becomes dysregulated, aiding our understanding of age-related susceptibility to lung disease.

TP63 fusions drive multicomplex enhancer rewiring, lymphomagenesis, and EZH2 dependence.

Gene fusions involving tumor protein p63 gene (TP63) occur in multiple T and B cell lymphomas and portend a dismal prognosis for patients. The function and mechanisms of TP63 fusions remain unclear, and there is no target therapy for patients with lymphoma harboring TP63 fusions. Here, we show that TP63 fusions act as bona fide oncogenes and are essential for fusion-positive lymphomas. Transgenic mice expressing TBL1XR1::TP63, the most common TP63 fusion, develop diverse lymphomas that recapitulate multiple human T and B cell lymphomas. Here, we identify that TP63 fusions coordinate the recruitment of two epigenetic modifying complexes, the nuclear receptor corepressor (NCoR)-histone deacetylase 3 (HDAC3) by the N-terminal TP63 fusion partner and the lysine methyltransferase 2D (KMT2D) by the C-terminal TP63 component, which are both required for fusion-dependent survival. TBL1XR1::TP63 localization at enhancers drives a unique cell state that involves up-regulation of MYC and the polycomb repressor complex 2 (PRC2) components EED and EZH2. Inhibiting EZH2 with the therapeutic agent valemetostat is highly effective at treating transgenic lymphoma murine models, xenografts, and patient-derived xenografts harboring TP63 fusions. One patient with TP63-rearranged lymphoma showed a rapid response to valemetostat treatment. In summary, TP63 fusions link partner components that, together, coordinate multiple epigenetic complexes, resulting in therapeutic vulnerability to EZH2 inhibition.

Low-dose interleukin-2 promotes immune regulation in face transplantation: A pilot study.

Face transplantation is a life-changing procedure for patients with severe composite facial defects. However, it is hampered by high acute rejection rates due to the immunogenicity of skin allograft and toxicity linked to high doses of immunosuppression. To reduce immunosuppression-associated complications, we, for the first time in face transplant recipients, used low-dose interleukin 2 (IL-2) therapy to expand regulatory T cells (Tregs) in vivo and to enhance immune modulation, under close immunological monitoring of peripheral blood and skin allograft. Low-dose IL-2 achieved a sustained expansion (∼4-fold to 5-fold) of circulating Tregs and a reduction (∼3.5-fold) of B cells. Post-IL-2 Tregs exhibited greater suppressive function, characterized by higher expression of TIM-3 and LAG3co-inhibitory molecules. In the skin allograft, Tregs increased after low-dose IL-2 therapy. IL-2 induced a distinct molecular signature in the allograft with reduced cytotoxicity-associated genes (granzyme B and perforin). Two complications were observed during the trial: one rejection event and an episode of autoimmune hemolytic anemia. In summary, this initial experience demonstrated that low-dose IL-2 therapy was not only able to promote immune regulation in face transplant recipients but also highlighted challenges related to its narrow therapeutic window. More specific targeted Treg expansion strategies are needed to translate this approach to the clinic.

Transcription factor SOX15 regulates stem cell pluripotency and promotes neural fate during differentiation by activating the neurogenic gene Hes5.

SOX2 and SOX15 are Sox family transcription factors enriched in embryonic stem cells (ESCs). The role of SOX2 in activating gene expression programs essential for stem cell self-renewal and acquisition of pluripotency during somatic cell reprogramming is well-documented. However, the contribution of SOX15 to these processes is unclear and often presumed redundant with SOX2 largely because overexpression of SOX15 can partially restore self-renewal in SOX2-deficient ESCs. Here, we show that SOX15 contributes to stem cell maintenance by cooperating with ESC-enriched transcriptional coactivators to ensure optimal expression of pluripotency-associated genes. We demonstrate that SOX15 depletion compromises reprogramming of fibroblasts to pluripotency which cannot be compensated by SOX2. Ectopic expression of SOX15 promotes the reversion of a postimplantation, epiblast stem cell state back to a preimplantation, ESC-like identity even though SOX2 is expressed in both cell states. We also uncover a role of SOX15 in lineage specification, by showing that loss of SOX15 leads to defects in commitment of ESCs to neural fates. SOX15 promotes neural differentiation by binding to and activating a previously uncharacterized distal enhancer of a key neurogenic regulator, Hes5. Together, these findings identify a multifaceted role of SOX15 in induction and maintenance of pluripotency and neural differentiation.

Ehf and Fezf2 regulate late medullary thymic epithelial cell and thymic tuft cell development.

Thymic epithelial cells are indispensable for T cell maturation and selection and the induction of central immune tolerance. The self-peptide repertoire expressed by medullary thymic epithelial cells is in part regulated by the transcriptional regulator Aire (Autoimmune regulator) and the transcription factor Fezf2. Due to the high complexity of mTEC maturation stages (i.e., post-Aire, Krt10+ mTECs, and Dclk1+ Tuft mTECs) and the heterogeneity in their gene expression profiles (i.e., mosaic expression patterns), it has been challenging to identify the additional factors complementing the transcriptional regulation. We aimed to identify the transcriptional regulators involved in the regulation of mTEC development and self-peptide expression in an unbiased and genome-wide manner. We used ATAC footprinting analysis as an indirect approach to identify transcription factors involved in the gene expression regulation in mTECs, which we validated by ChIP sequencing. This study identifies Fezf2 as a regulator of the recently described thymic Tuft cells (i.e., Tuft mTECs). Furthermore, we identify that transcriptional regulators of the ELF, ESE, ERF, and PEA3 subfamily of the ETS transcription factor family and members of the Krüppel-like family of transcription factors play a role in the transcriptional regulation of genes involved in late mTEC development and promiscuous gene expression.

Whole-genome CRISPR screening identifies genetic manipulations to reduce immune rejection of stem cell-derived islets.

Human embryonic stem cells (hESCs) provide opportunities for cell replacement therapy of insulin-dependent diabetes. Therapeutic quantities of human stem cell-derived islets (SC-islets) can be produced by directed differentiation. However, preventing allo-rejection and recurring autoimmunity, without the use of encapsulation or systemic immunosuppressants, remains a challenge. An attractive approach is to transplant SC-islets, genetically modified to reduce the impact of immune rejection. To determine the underlying forces that drive immunogenicity of SC-islets in inflammatory environments, we performed single-cell RNA sequencing (scRNA-seq) and whole-genome CRISPR screen of SC-islets under immune interaction with allogeneic peripheral blood mononuclear cells (PBMCs). Data analysis points to "alarmed" populations of SC-islets that upregulate genes in the interferon (IFN) pathway. The CRISPR screen in vivo confirms that targeting IFNγ-induced mediators has beneficial effects on SC-islet survival under immune attack. Manipulating the IFN response by depleting chemokine ligand 10 (CXCL10) in SC-islet grafts confers improved survival against allo-rejection compared with wild-type grafts in humanized mice. These results offer insights into the nature of immune destruction of SC-islets during allogeneic responses and provide targets for gene editing.

Nanchangmycin regulates FYN, PTK2, and MAPK1/3 to control the fibrotic activity of human hepatic stellate cells.

Chronic liver injury causes fibrosis, characterized by the formation of scar tissue resulting from excessive accumulation of extracellular matrix (ECM) proteins. Hepatic stellate cell (HSC) myofibroblasts are the primary cell type responsible for liver fibrosis, yet there are currently no therapies directed at inhibiting the activity of HSC myofibroblasts. To search for potential anti-fibrotic compounds, we performed a high-throughput compound screen in primary human HSC myofibroblasts and identified 19 small molecules that induce HSC inactivation, including the polyether ionophore nanchangmycin (NCMC). NCMC induces lipid re-accumulation while reducing collagen expression, deposition of collagen in the extracellular matrix, cell proliferation, and migration. We find that NCMC increases cytosolic Ca2+ and reduces the phosphorylated protein levels of FYN, PTK2 (FAK), MAPK1/3 (ERK2/1), HSPB1 (HSP27), and STAT5B. Further, depletion of each of these kinases suppress COL1A1 expression. These studies reveal a signaling network triggered by NCMC to inactivate HSC myofibroblasts and reduce expression of proteins that compose the fibrotic scar. Identification of the antifibrotic effects of NCMC and the elucidation of pathways by which NCMC inhibits fibrosis provide new tools and therapeutic targets that could potentially be utilized to combat the development and progression of liver fibrosis.

Molecular characterization and cell type composition deconvolution of fibrosis in NAFLD.

Non-alcoholic fatty liver disease (NAFLD) is the most common cause of liver disease worldwide. In adults with NAFLD, fibrosis can develop and progress to liver cirrhosis and liver failure. However, the underlying molecular mechanisms of fibrosis progression are not fully understood. Using total RNA-Seq, we investigated the molecular mechanisms of NAFLD and fibrosis. We sequenced liver tissue from 143 adults across the full spectrum of fibrosis stage including those with stage 4 fibrosis (cirrhosis). We identified gene expression clusters that strongly correlate with fibrosis stage including four genes that have been found consistently across previously published transcriptomic studies on NASH i.e. COL1A2, EFEMP2, FBLN5 and THBS2. Using cell type deconvolution, we estimated the loss of hepatocytes versus gain of hepatic stellate cells, macrophages and cholangiocytes with advancing fibrosis stage. Hepatocyte-specific functional analysis indicated increase of pro-apoptotic pathways and markers of bipotent hepatocyte/cholangiocyte precursors. Regression modelling was used to derive predictors of fibrosis stage. This study elucidated molecular and cell composition changes associated with increasing fibrosis stage in NAFLD and defined informative gene signatures for the disease.

Immunoregulatory and lipid presentation pathways are upregulated in human face transplant rejection.

BACKGROUNDRejection is the primary barrier to broader implementation of vascularized composite allografts (VCAs), including face and limb transplants. The immunologic pathways activated in face transplant rejection have not been fully characterized.METHODSUsing skin biopsies prospectively collected over 9 years from 7 face transplant patients, we studied rejection by gene expression profiling, histology, immunostaining, and T cell receptor sequencing.RESULTSGrade 1 rejection did not differ significantly from nonrejection, suggesting that it does not represent a pathologic state. In grade 2, there was a balanced upregulation of both proinflammatory T cell activation pathways and antiinflammatory checkpoint and immunomodulatory pathways, with a net result of no tissue injury. In grade 3, IFN-γ-driven inflammation, antigen-presenting cell activation, and infiltration of the skin by proliferative T cells bearing markers of antigen-specific activation and cytotoxicity tipped the balance toward tissue injury. Rejection of VCAs and solid organ transplants had both distinct and common features. VCA rejection was uniquely associated with upregulation of immunoregulatory genes, including SOCS1; induction of lipid antigen-presenting CD1 proteins; and infiltration by T cells predicted to recognize CD1b and CD1c.CONCLUSIONOur findings suggest that the distinct features of VCA rejection reflect the unique immunobiology of skin and that enhancing cutaneous immunoregulatory networks may be a useful strategy in combatting rejection.Trial registrationClinicalTrials.gov NCT01281267.FUNDINGAssistant Secretary of Defense and Health Affairs, through Reconstructive Transplant Research (W81XWH-17-1-0278, W81XWH-16-1-0647, W81XWH-16-1-0689, W81XWH-18-1-0784, W81XWH-1-810798); American Society of Transplantation's Transplantation and Immunology Research Network Fellowship Research Grant; Plastic Surgery Foundation Fellowship from the American Society of Plastic Surgeons; Novo Nordisk Foundation (NNF15OC0014092); Lundbeck Foundation; Aage Bangs Foundation; A.P. Moller Foundation for the Advancement of Medical Science; NIH UL1 RR025758.

Drosophila PDGF/VEGF signaling from muscles to hepatocyte-like cells protects against obesity.

PDGF/VEGF ligands regulate a plethora of biological processes in multicellular organisms via autocrine, paracrine, and endocrine mechanisms. We investigated organ-specific metabolic roles of Drosophila PDGF/VEGF-like factors (Pvfs). We combine genetic approaches and single-nuclei sequencing to demonstrate that muscle-derived Pvf1 signals to the Drosophila hepatocyte-like cells/oenocytes to suppress lipid synthesis by activating the Pi3K/Akt1/TOR signaling cascade in the oenocytes. Functionally, this signaling axis regulates expansion of adipose tissue lipid stores in newly eclosed flies. Flies emerge after pupation with limited adipose tissue lipid stores and lipid level is progressively accumulated via lipid synthesis. We find that adult muscle-specific expression of pvf1 increases rapidly during this stage and that muscle-to-oenocyte Pvf1 signaling inhibits expansion of adipose tissue lipid stores as the process reaches completion. Our findings provide the first evidence in a metazoan of a PDGF/VEGF ligand acting as a myokine that regulates systemic lipid homeostasis by activating TOR in hepatocyte-like cells.

Translocation of Viable Gut Microbiota to Mesenteric Adipose Drives Formation of Creeping Fat in Humans.

A mysterious feature of Crohn's disease (CD) is the extra-intestinal manifestation of "creeping fat" (CrF), defined as expansion of mesenteric adipose tissue around the inflamed and fibrotic intestine. In the current study, we explore whether microbial translocation in CD serves as a central cue for CrF development. We discovered a subset of mucosal-associated gut bacteria that consistently translocated and remained viable in CrF in CD ileal surgical resections, and identified Clostridium innocuum as a signature of this consortium with strain variation between mucosal and adipose isolates, suggesting preference for lipid-rich environments. Single-cell RNA sequencing characterized CrF as both pro-fibrotic and pro-adipogenic with a rich milieu of activated immune cells responding to microbial stimuli, which we confirm in gnotobiotic mice colonized with C. innocuum. Ex vivo validation of expression patterns suggests C. innocuum stimulates tissue remodeling via M2 macrophages, leading to an adipose tissue barrier that serves to prevent systemic dissemination of bacteria.

Effects of tesamorelin on hepatic transcriptomic signatures in HIV-associated NAFLD.

Nonalcoholic fatty liver disease (NAFLD) is a common comorbidity among people living with HIV that has a more aggressive course than NAFLD among the general population. In a recent randomized placebo-controlled trial, we demonstrated that the growth hormone-releasing hormone analog tesamorelin reduced liver fat and prevented fibrosis progression in HIV-associated NAFLD over 1 year. As such, tesamorelin is the first strategy that has shown to be effective against NAFLD among the population with HIV. The current study leveraged paired liver biopsy specimens from this trial to identify hepatic gene pathways that are differentially modulated by tesamorelin versus placebo. Using gene set enrichment analysis, we found that tesamorelin increased hepatic expression of hallmark gene sets involved in oxidative phosphorylation and decreased hepatic expression of gene sets contributing to inflammation, tissue repair, and cell division. Tesamorelin also reciprocally up- and downregulated curated gene sets associated with favorable and poor hepatocellular carcinoma prognosis, respectively. Notably, among tesamorelin-treated participants, these changes in hepatic expression correlated with improved fibrosis-related gene score. Our findings inform our knowledge of the biology of pulsatile growth hormone action and provide a mechanistic basis for the observed clinical effects of tesamorelin on the liver.

A single-cell survey of Drosophila blood.

Drosophila blood cells, called hemocytes, are classified into plasmatocytes, crystal cells, and lamellocytes based on the expression of a few marker genes and cell morphologies, which are inadequate to classify the complete hemocyte repertoire. Here, we used single-cell RNA sequencing (scRNA-seq) to map hemocytes across different inflammatory conditions in larvae. We resolved plasmatocytes into different states based on the expression of genes involved in cell cycle, antimicrobial response, and metabolism together with the identification of intermediate states. Further, we discovered rare subsets within crystal cells and lamellocytes that express fibroblast growth factor (FGF) ligand branchless and receptor breathless, respectively. We demonstrate that these FGF components are required for mediating effective immune responses against parasitoid wasp eggs, highlighting a novel role for FGF signaling in inter-hemocyte crosstalk. Our scRNA-seq analysis reveals the diversity of hemocytes and provides a rich resource of gene expression profiles for a systems-level understanding of their functions.

A cell atlas of the adult Drosophila midgut.

Studies of the adult Drosophila midgut have led to many insights in our understanding of cell-type diversity, stem cell regeneration, tissue homeostasis, and cell fate decision. Advances in single-cell RNA sequencing provide opportunities to identify new cell types and molecular features. We used single-cell RNA sequencing to characterize the transcriptome of midgut epithelial cells and identified 22 distinct clusters representing intestinal stem cells, enteroblasts, enteroendocrine cells (EEs), and enterocytes. This unbiased approach recovered most of the known intestinal stem cells/enteroblast and EE markers, highlighting the high quality of the dataset, and led to insights on intestinal stem cell biology, cell type-specific organelle features, the roles of new transcription factors in progenitors and regional variation along the gut, 5 additional EE gut hormones, EE hormonal expression diversity, and paracrine function of EEs. To facilitate mining of this rich dataset, we provide a web-based resource for visualization of gene expression in single cells. Altogether, our study provides a comprehensive resource for addressing functions of genes in the midgut epithelium.