Changlin Wan

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.

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.

Immune evasion is a pivotal event in tumor progression. To eliminate human cancer cells, current immune checkpoint therapy is set to boost CD8+ T cell-mediated cytotoxicity. However, this action is eventually dependent on the efficient recognition of tumor-specific antigens via T cell receptors. One primary mechanism by which tumor cells evade immune surveillance is to downregulate their antigen presentation. Little progress has been made toward harnessing potential therapeutic targets for enhancing antigen presentation on the tumor cell. Here, we identified MAL2 as a key player that determines the turnover of the antigen-loaded MHC-I complex and reduces the antigen presentation on tumor cells. MAL2 promotes the endocytosis of tumor antigens via direct interaction with the MHC-I complex and endosome-associated RAB proteins. In preclinical models, depletion of MAL2 in breast tumor cells profoundly enhanced the cytotoxicity of tumor-infiltrating CD8+ T cells and suppressed breast tumor growth, suggesting that MAL2 is a potential therapeutic target for breast cancer immunotherapy.

A key challenge in modeling single-cell RNA-seq data is to capture the diversity of gene expression states regulated by different transcriptional regulatory inputs across individual cells, which is further complicated by largely observed zero and low expressions. We developed a left truncated mixture Gaussian (LTMG) model, from the kinetic relationships of the transcriptional regulatory inputs, mRNA metabolism and abundance in single cells. LTMG infers the expression multi-modalities across single cells, meanwhile, the dropouts and low expressions are treated as left truncated. We demonstrated that LTMG has significantly better goodness of fitting on an extensive number of scRNA-seq data, comparing to three other state-of-the-art models. Our biological assumption of the low non-zero expressions, rationality of the multimodality setting, and the capability of LTMG in extracting expression states specific to cell types or functions, are validated on independent experimental data sets. A differential gene expression test and a co-regulation module identification method are further developed. We experimentally validated that our differential expression test has higher sensitivity and specificity, compared with other five popular methods. The co-regulation analysis is capable of retrieving gene co-regulation modules corresponding to perturbed transcriptional regulations. A user-friendly R package with all the analysis power is available at https://github.com/zy26/LTMGSCA.