The metabolic heterogeneity and metabolic interplay between cells are known as significant contributors to disease treatment resistance. However, with the lack of a mature high-throughput single-cell metabolomics technology, we are yet to establish systematic understanding of the intra-tissue metabolic heterogeneity and cooperative mechanisms. To mitigate this knowledge gap, we developed a novel computational method, namely, single-cell flux estimation analysis (scFEA), to infer the cell-wise fluxome from single-cell RNA-sequencing (scRNA-seq) data. scFEA is empowered by a systematically reconstructed human metabolic map as a factor graph, a novel probabilistic model to leverage the flux balance constraints on scRNA-seq data, and a novel graph neural network–based optimization solver. The intricate information cascade from transcriptome to metabolome was captured using multilayer neural networks to capitulate the nonlinear dependency between enzymatic gene expressions and reaction rates. We experimentally validated scFEA by generating an scRNA-seq data set with matched metabolomics data on cells of perturbed oxygen and genetic conditions. Application of scFEA on this data set showed the consistency between predicted flux and the observed variation of metabolite abundance in the matched metabolomics data. We also applied scFEA on five publicly available scRNA-seq and spatial transcriptomics data sets and identified context- and cell group–specific metabolic variations. The cell-wise fluxome predicted by scFEA empowers a series of downstream analyses including identification of metabolic modules or cell groups that share common metabolic variations, sensitivity evaluation of enzymes with regards to their impact on the whole metabolic flux, and inference of cell–tissue and cell–cell metabolic communications.
Single-cell sequencing studies have characterized the transcriptomic signature of cell types within the kidney. However, the spatial distribution of acute kidney injury (AKI) is regional and affects cells heterogeneously. We first optimized coordination of spatial transcriptomics and single-nuclear sequencing data sets, mapping 30 dominant cell types to a human nephrectomy. The predicted cell-type spots corresponded with the underlying histopathology. To study the implications of AKI on transcript expression, we then characterized the spatial transcriptomic signature of 2 murine AKI models: ischemia/reperfusion injury (IRI) and cecal ligation puncture (CLP). Localized regions of reduced overall expression were associated with injury pathways. Using single-cell sequencing, we deconvoluted the signature of each spatial transcriptomic spot, identifying patterns of colocalization between immune and epithelial cells. Neutrophils infiltrated the renal medulla in the ischemia model. Atf3 was identified as a chemotactic factor in S3 proximal tubules. In the CLP model, infiltrating macrophages dominated the outer cortical signature, and Mdk was identified as a corresponding chemotactic factor. The regional distribution of these immune cells was validated with multiplexed CO-Detection by indEXing (CODEX) immunofluorescence. Spatial transcriptomic sequencing complemented single-cell sequencing by uncovering mechanisms driving immune cell infiltration and detection of relevant cell subpopulations.
Retinal organoids are three-dimensional structures derived from human pluripotent stem cells (hPSCs) which recapitulate the spatial and temporal differentiation of the retina, serving as effective in vitro models of retinal development. However, a lack of emphasis has been placed upon the development and organization of retinal ganglion cells (RGCs) within retinal organoids. Thus, initial efforts were made to characterize RGC differentiation throughout early stages of organoid development, with a clearly defined RGC layer developing in a temporally-appropriate manner expressing a complement of RGC-associated markers. Beyond studies of RGC development, retinal organoids may also prove useful for cellular replacement in which extensive axonal outgrowth is necessary to reach post-synaptic targets. Organoid-derived RGCs could help to elucidate factors promoting axonal outgrowth, thereby identifying approaches to circumvent a formidable obstacle to RGC replacement. As such, additional efforts demonstrated significant enhancement of neurite outgrowth through modulation of both substrate composition and growth factor signaling. Additionally, organoid-derived RGCs exhibited diverse phenotypes, extending elaborate growth cones and expressing numerous guidance receptors. Collectively, these results establish retinal organoids as a valuable tool for studies of RGC development, and demonstrate the utility of organoid-derived RGCs as an effective platform to study factors influencing neurite outgrowth from organoid-derived RGCs.
One of the primary mechanisms of tumor cell immune evasion is the loss of antigenicity, which arises due to lack of immunogenic tumor antigens as well as dysregulation of the antigen processing machinery. In a screen for small-molecule compounds from herbal medicine that potentiate T cell–mediated cytotoxicity, we identified atractylenolide I (ATT-I), which substantially promotes tumor antigen presentation of both human and mouse colorectal cancer (CRC) cells and thereby enhances the cytotoxic response of CD8+ T cells. Cellular thermal shift assay (CETSA) with multiplexed quantitative mass spectrometry identified the proteasome 26S subunit non–ATPase 4 (PSMD4), an essential component of the immunoproteasome complex, as a primary target protein of ATT-I. Binding of ATT-I with PSMD4 augments the antigen-processing activity of immunoproteasome, leading to enhanced MHC-I–mediated antigen presentation on cancer cells. In syngeneic mouse CRC models and human patient–derived CRC organoid models, ATT-I treatment promotes the cytotoxicity of CD8+ T cells and thus profoundly enhances the efficacy of immune checkpoint blockade therapy. Collectively, we show here that targeting the function of immunoproteasome with ATT-I promotes tumor antigen presentation and empowers T cell cytotoxicity, thus elevating the tumor response to immunotherapy.
Characterized by the accumulation of somatic mutations in blood cell lineages, clonal hematopoiesis of indeterminate potential (CHIP) is frequent in aging and involves the expansion of mutated hematopoietic stem and progenitor cells (HSC/Ps) that leads to an increased risk of hematologic malignancy. However, the risk factors that contribute to CHIP-associated clonal hematopoiesis (CH) are poorly understood. Obesity induces a proinflammatory state and fatty bone marrow (FBM), which may influence CHIP-associated pathologies. We analyzed exome sequencing and clinical data for 47,466 individuals with validated CHIP in the UK Biobank. CHIP was present in 5.8% of the study population and was associated with a significant increase in the waist-to-hip ratio (WHR). Mouse models of obesity and CHIP driven by heterozygosity of Tet2, Dnmt3a, Asxl1, and Jak2 resulted in exacerbated expansion of mutant HSC/Ps due in part to excessive inflammation. Our results show that obesity is highly associated with CHIP and that a proinflammatory state could potentiate the progression of CHIP to more significant hematologic neoplasia. The calcium channel blockers nifedipine and SKF-96365, either alone or in combination with metformin, MCC950, or anakinra (IL-1 receptor antagonist), suppressed the growth of mutant CHIP cells and partially restored normal hematopoiesis. Targeting CHIP-mutant cells with these drugs could be a potential therapeutic approach to treat CH and its associated abnormalities in individuals with obesity.