Dynamics of The Γδtcr Repertoires During The Dedifferentiation Process and Pilot Implications for Immunotherapy of Thyroid Cancer.

γδ T cells are evolutionarily conserved T lymphocytes that manifest unique antitumor efficacy independent of tumor mutation burden (TMB) and conventional human leukocyte antigen (HLA) recognition. However, the dynamic changes in their T cell receptor (TCR) repertoire during cancer progression and treatment courses remain unclear. Here, a comprehensive characterization of γδTCR repertoires are performed in thyroid cancers with divergent differentiation states through cross-sectional studies. The findings revealed a significant correlation between the differentiation states and TCR repertoire diversity. Notably, highly expanded clones are prominently enriched in γδ T cell compartment of dedifferentiated patients. Moreover, by longitudinal investigations of the γδ T cell response to various antitumor therapies, it is found that the emergence and expansion of the Vδ2neg subset may be potentially associated with favorable clinical outcomes after post-radiotherapeutic immunotherapy. These findings are further validated at single-cell resolution in both advanced thyroid cancer patients and a murine model, underlining the importance of further investigations into the role of γδTCR in cancer immunity and therapeutic strategies.

Liver cancer stem cell dissemination and metastasis: uncovering the role of NRCAM in hepatocellular carcinoma.

BACKGROUND: Liver cancer stem cells (LCSCs) play an important role in hepatocellular carcinoma (HCC), but the mechanisms that link LCSCs to HCC metastasis remain largely unknown. This study aims to reveal the contributions of NRCAM to LCSC function and HCC metastasis, and further explore its mechanism in detail. METHODS: 117 HCC and 29 non-HCC patients with focal liver lesions were collected and analyzed to assess the association between NRCAM and HCC metastasis. Single-cell RNA sequencing (scRNA-seq) was used to explore the biological characteristics of cells with high NRCAM expression in metastatic HCC. The role and mechanism of NRCAM in LCSC dissemination and metastasis was explored in vitro and in vivo using MYC-driven LCSC organoids from murine liver cells. RESULTS: Serum NRCAM is associated with HCC metastasis and poor prognosis. A scRNA-seq analysis identified that NRCAM was highly expressed in LCSCs with MYC activation in metastatic HCC. Moreover, NRCAM facilitated LCSC migration and invasion, which was confirmed in MYC-driven LCSC organoids. The in vivo tumor allografts demonstrated that NRCAM mediated intra-hepatic/lung HCC metastasis by enhancing the ability of LCSCs to escape from tumors into the bloodstream. Nrcam expression inhibition in LCSCs blocked HCC metastasis. Mechanistically, NRCAM activated epithelial-mesenchymal transition (EMT) and metastasis-related matrix metalloproteinases (MMPs) through the MACF1 mediated ß-catenin signaling pathway in LCSCs. CONCLUSIONS: LCSCs typified by high NRCAM expression have a strong ability to invade and migrate, which is an important factor leading to HCC metastasis.

Genome-scale functional genomics identify genes preferentially essential for multiple myeloma cells compared to other neoplasias.

Clinical progress in multiple myeloma (MM), an incurable plasma cell (PC) neoplasia, has been driven by therapies that have limited applications beyond MM/PC neoplasias and do not target specific oncogenic mutations in MM. Instead, these agents target pathways critical for PC biology yet largely dispensable for malignant or normal cells of most other lineages. Here we systematically characterized the lineage-preferential molecular dependencies of MM through genome-scale clustered regularly interspaced short palindromic repeats (CRISPR) studies in 19 MM versus hundreds of non-MM lines and identified 116 genes whose disruption more significantly affects MM cell fitness compared with other malignancies. These genes, some known, others not previously linked to MM, encode transcription factors, chromatin modifiers, endoplasmic reticulum components, metabolic regulators or signaling molecules. Most of these genes are not among the top amplified, overexpressed or mutated in MM. Functional genomics approaches thus define new therapeutic targets in MM not readily identifiable by standard genomic, transcriptional or epigenetic profiling analyses.

Generation of a Mouse Spontaneous Autoimmune Thyroiditis Model.

In recent years, Hashimoto's thyroiditis (HT) has become the most common autoimmune thyroid disease. It is characterized by lymphocyte infiltration and the detection of specific serum autoantibodies. Although the potential mechanism is still not clear, the risk of Hashimoto's thyroiditis is related to genetic and environmental factors. At present, there are several types of models of autoimmune thyroiditis, including experimental autoimmune thyroiditis (EAT) and spontaneous autoimmune thyroiditis (SAT). EAT in mice is a common model for HT, which is immunized with lipopolysaccharide (LPS) combined with thyroglobulin (Tg) or supplemented with complete Freund's adjuvant (CFA). The EAT mouse model is widely established in many types of mice. However, the disease progression is more likely associated with the Tg antibody response, which may vary in different experiments. SAT is also widely used in the study of HT in the NOD.H-2h4 mouse. The NOD.H2h4 mouse is a new strain obtained from the cross of the nonobese diabetic (NOD) mouse with the B10.A(4R), which is significantly induced for HT with or without feeding iodine. During the induction, the NOD.H-2h4 mouse has a high level of TgAb accompanied by lymphocyte infiltration in the thyroid follicular tissue. However, for this type of mouse model, there are few studies to comprehensively evaluate the pathological process during the induction of iodine. A SAT mouse model for HT research is established in this study, and the pathologic changing process is evaluated after a long period of iodine induction. Through this model, researchers can better understand the pathological development of HT and screen new treatment methods for HT.

Atp8a1 deletion increases the proliferative activity of hematopoietic stem cells by impairing PTEN function.

PURPOSE: The eukaryotic cell plasma membrane contains several asymmetrically distributed phospholipids, which is maintained by the P4-ATPase flippase complex. Herein, we demonstrated the biological effects and mechanisms of asymmetrical loss in hematopoietic stem cells (HSCs). METHODS: An Atp8a1 knockout mouse model was employed, from which the HSC (long-term HSCs and short-term HSCs) population was analyzed to assess their abundance and function. Additionally, competitive bone marrow transplantation and 5-FU stress assays were performed. RNA sequencing was performed on Hematopoietic Stem and Progenitor Cells, and DNA damage was assayed using immunofluorescence staining and comet electrophoresis. The protein abundance for members of key signaling pathways was confirmed using western blotting. RESULTS: Atp8a1 deletion resulted in slight hyperleukocytosis, associated with the high proliferation of HSCs and BCR/ABL1 transformed leukemia stem cells (LSCs). Atp8a1 deletion increased the repopulation capability of HSCs with a competitive advantage in reconstitution assay. HSCs without Atp8a1 were more sensitive to 5-FU-induced apoptosis. Moreover, Atp8a1 deletion prevented HSC DNA damage and facilitated DNA repair processes. Genes involved in PI3K-AKT-mTORC1, DNA repair, and AP-1 complex signaling were enriched and elevated in HSCs with Atp8a1 deletion. Furthermore, Atp8a1 deletion caused decreased PTEN protein levels, resulting in the activation of PI3K-AKT-mTORC1 signaling, further increasing the activity of JNK/AP-1 signaling and YAP1 phosphorylation. CONCLUSION: We identified the role of Atp8a1 on hematopoiesis and HSCs. Atp8a1 deletion resulted in the loss of phosphatidylserine asymmetry and intracellular signal transduction chaos.

TSPAN32 suppresses chronic myeloid leukemia pathogenesis and progression by stabilizing PTEN.

We report herein that TSPAN32 is a key node factor for Philadelphia (Ph+) leukemia pathogenesis. We found that TSPAN32 expression was repressed by BCR-ABL and ectopic TSPAN32 expression upon Imatinib treatment inhibited the proliferation of Ph+ cell lines. Tspan32 overexpression significantly prevented BCR-ABL induced leukemia progression in a murine model and impaired leukemia stem cell (LSC) proliferation. LSCs represent an obstacle for chronic myeloid leukemia (CML) elimination, which continually replenish leukemia cells and are associated with disease relapse. Therefore, the identification of essential targets that contribute to the survival and self-renewal of LSCs is important for novel curative CML. Mechanistically, TSPAN32 was shown to interact with PTEN, increased its protein level and caused a reduction in PI3K-AKT signaling activity. We also found that TSPAN32 was repressed by BCR-ABL via the suppression of an important transcription factor, TAL1. Ectopic expression of TAL1 significantly increased TSPAN32 mRNA and protein level, which indicated that BCR-ABL repressed TSPAN32 transcription by decreasing TAL1 expression. Overall, we identified a new signaling axis composed of "BCR-ABL-TAL1-TSPAN32-PTEN-PI3K-AKT". Our findings further complement the known mechanisms underlying the transformation potential of BCR-ABL in CML pathogenesis. This new signaling axis also provides a potential means to target PI3K-AKT for CML treatment.

HDAC I/IIb selective inhibitor Purinostat Mesylate combined with GLS1 inhibition effectively eliminates CML stem cells.

Purinostat Mesylate (PM) is a novel highly selective and active HDAC I/IIb inhibitor, and the injectable formulation of PM (PMF) based on the compound prescription containing cyclodextrin completely can overcome PM's poor solubility and improves its stability and pharmacokinetic properties. Here, we showed that PM effectively repressed the survival of Ph+ leukemia cells and CD34+ leukemia cells from CML patients in vitro. In vivo studies demonstrated that PMF significantly prevented BCR-ABL(T315I) induced CML progression by restraining leukemia stem cells (LSCs), which are insensitive to chemotherapy and responsible for CML relapse. Mechanism studies revealed that targeting HDAC I/IIb repressed several important factors for LSCs survival including c-Myc, ß-Catenin, E2f, Ezh2, Alox5, and mTOR, as well as interrupted some critical biologic processes. Additionally, PMF increased glutamate metabolism in LSCs by increasing GLS1. The combination of PMF and glutaminase inhibitor BPTES synergistically eradicated LSCs by altering multiple key proteins and signaling pathways which are critical for LSC survival and self-renewal. Overall, our findings represent a new therapeutic strategy for eliminating LSCs by targeting HDAC I/IIb and glutaminolysis, which potentially provides a guidance for PMF clinical trials in the future for TKI resistance CML patients.

Pan-cancer single-cell analysis reveals the heterogeneity and plasticity of cancer-associated fibroblasts in the tumor microenvironment.

Cancer-associated fibroblasts (CAFs) are the predominant components of the tumor microenvironment (TME) and influence cancer hallmarks, but without systematic investigation on their ubiquitous characteristics across different cancer types. Here, we perform pan-cancer analysis on 226 samples across 10 solid cancer types to profile the TME at single-cell resolution, illustrating the commonalities/plasticity of heterogenous CAFs. Activation trajectory of the major CAF types is divided into three states, exhibiting distinct interactions with other cell components, and relating to prognosis of immunotherapy. Moreover, minor CAF components represent the alternative origin from other TME components (e.g., endothelia and macrophages). Particularly, the ubiquitous presentation of endothelial-to-mesenchymal transition CAF, which may interact with proximal SPP1+ tumor-associated macrophages, is implicated in endothelial-to-mesenchymal transition and survival stratifications. Our study comprehensively profiles the shared characteristics and dynamics of CAFs, and highlight their heterogeneity and plasticity across different cancer types. Browser of integrated pan-cancer single-cell information is available at https://gist-fgl.github.io/sc-caf-atlas/ .

Pancreatic tumor eradication via selective Pin1 inhibition in cancer-associated fibroblasts and T lymphocytes engagement.

Cancer associated fibroblasts (CAFs) support tumors via multiple mechanisms, including maintaining the immunosuppressive tumor microenvironment and limiting infiltration of immune cells. The prolyl isomerase Pin1, whose overexpression in CAFs has not been fully profiled yet, plays critical roles in tumor initiation and progression. To decipher effects of selective Pin1 inhibition in CAFs on pancreatic cancer, here we formulate a DNA-barcoded micellular system (DMS) encapsulating the Pin1 inhibitor AG17724. DMS functionalized with CAF-targeting anti-FAP-α antibodies (antiCAFs-DMS) can selectively inhibit Pin1 in CAFs, leading to efficacious but transient tumor growth inhibition. We further integrate DNA aptamers (AptT), which can engage CD8+ T lymphocytes, to obtain a bispecific antiCAFs-DMS-AptT system. AntiCAFs-DMS-AptT inhibits tumor growth in subcutaneous and orthotopic pancreatic cancer models.

Dasatinib and interferon alpha synergistically induce pyroptosis-like cell death in philadelphia chromosome positive acute lymphoblastic leukemia.

Philadelphia chromosome-positive acute lymphoblastic leukemia (Ph+ALL) is a high-risk disease subtype with a dismal prognosis. Inhibiting BCR-ABL kinase alone is insufficient to eradicate Ph+ALL clones, and alternative BCR-ABL-dependent and -independent pathways need to be targeted as an effective strategy. Our study revealed that the combination of dasatinib and interferon-α showed synergistic activity against Ph+ALL, inducing mitochondrial dysfunction and causing necrosis-like cell lysis. Mechanistic studies showed that the induced cell death was caspase-3-independent. Canonical necroptosis signals, such as RIP1 and MLKL, were not activated; instead, the pyroptosis executor Gasdermin D was upregulated expression and activated. The expression levels of extracellular ATP and IL-1ß were also upregulated, both of which are markers of pyroptotic cell death. In a murine Ph+ALL model, the dual drug treatment prolonged the survival of tumor-bearing mice. More importantly, we incorporated the dual drugs to maintenance therapy in 39 patients who were unfit for allogeneic stem cell transplantation (allo-HSCT). The median follow-up was 28.5 months, the 4-year disease-free survival and overall survival rates were 52.2% and 65.2%, respectively. Our data suggest that the combination of dasatinib and interferon-α has potential synergistic activity against Ph+ALL and shows promise as a maintenance therapy for Ph+ALL patients who are unfit for allo-HSCT.

The suppressive functions of Rora in B lineage cell proliferation and BCR/ABL1-induced B-ALL pathogenesis.

RORA plays an important role in regulating circadian rhythms, inflammation, metabolism and cellular development. Herein, we explore the roles of Rora in B cell proliferation and differentiation, as well as in Ph+ B-ALL. By using Roraloxp/loxp Mx-1-Cre mice, Rora was deleted in hematopoietic cells post Pipc induction. Rora deficiency mice were associated with an obvious accumulation of B cells in the peripheral blood, bone marrow, and spleen. On the other hand, activation of Rora with Cholesterol sulfate (CS) was associated with decreased B cell numbers. RNA-seq analysis revealed that the transcription level of Lmo1 was decreased in Rora deficient B cells. Moreover, the expression of RORA was shown to be decreased in Ph+ B-ALL cells compared to peripheral blood derived B cells from healthy donors. The overexpression of Rora in BaF3 cells with BCR/ABL1 was also associated with impeded the cell growth and an increased apoptotic rate compared to cells transduced with BCR/ABL1 alone. The co-expression of BCR/ABL1 and Rora induced B-ALL mouse model was associated with the significant inhibition of BCR/ABL1-transformed cell growth and prolonged the survival of the diseased mice. These results suggest a novel role for Rora in B cell development and Ph+ leukemogenesis.

High-dose vitamin D metabolite delivery inhibits breast cancer metastasis.

Besides its well-known benefits on human health, calcitriol, the hormonally active form of vitamin D3, has been being evaluated in clinical trials as an anticancer agent. However, currently available results are contradictory and not fundamentally deciphered. To the best of our knowledge, hypercalcemia caused by high-dose calcitriol administration and its low bioavailability limit its anticancer investigations and translations. Here, we show that the one-step self-assembly of calcitriol and amphiphilic cholesterol-based conjugates leads to the formation of a stable minimalist micellar nanosystem. When administered to mice, this nanosystem demonstrates high calcitriol doses in breast tumor cells, significant tumor growth inhibition and antimetastasis capability, as well as good biocompatibility. We further reveal that the underlying molecular antimetastatic mechanisms involve downregulation of proteins facilitating metastasis and upregulation of paxillin, the key protein of focal adhesion, in primary tumors.

Differential m6A RNA landscapes across hematopoiesis reveal a role for IGF2BP2 in preserving hematopoietic stem cell function.

N6-methyladenosine (m6A) on mRNA plays critical roles in various cellular processes. However, the landscape and dynamics of m6A modification in hematopoietic system remain unknown. Here, we delineate a comprehensive m6A landscape across hematopoietic hierarchy and uncover that IGF2BP2 is required for preserving the function of hematopoietic stem cells (HSCs). Our data reveal a cell-type-specific m6A landscape in hematopoiesis. m6A modifications arise mostly in the early stage of hematopoiesis and prefer to play distinct roles for determining mRNA fates in HSCs and committed progenitors. Mechanistically, increased m6A-IGF2BP2 expression controls transcriptional state and maintenance of HSCs. IGF2BP2 deficiency induces quiescence loss and impairs HSC function. Moreover, IGF2BP2 loss increases mitochondrial activity of HSCs by accelerating Bmi1 mRNA decay, leading to de-repression of mitochondria-related genes. Collectively, our results present a fascinating portrait of m6A modification of hematopoietic hierarchy and reveal a key role of IGF2BP2 in maintaining HSC function by restraining mitochondrial activity.

ALCAM-EGFR interaction regulates myelomagenesis.

Multiple myeloma, a plasma cell malignancy in the bone marrow, remains largely incurable with currently available therapeutics. In this study, we discovered that the activated leukocyte cell adhesion molecule (ALCAM) interacted with epidermal growth factor receptor (EGFR), and regulated myelomagenesis. ALCAM was a negative regulator of myeloma clonogenicity. ALCAM expression was positively correlated with patients' survival. ALCAM-knockdown myeloma cells displayed enhanced colony formation in the presence of bone marrow stromal cells (BMSCs). BMSCs supported myeloma colony formation by secreted epidermal growth factor (EGF), which bound with its receptor (EGFR) on myeloma cells and activated Mek/Erk cell signaling, PI3K/Akt cell signaling, and hedgehog pathway. ALCAM could also bind with EGFR, block EGF from binding to EGFR, and abolish EGFR-initiated cell signaling. Hence, our study identifies ALCAM as a novel negative regulator of myeloma pathogenesis.

Characterizing dedifferentiation of thyroid cancer by integrated analysis.

Understanding of dedifferentiation, an indicator of poo prognosis for patients with thyroid cancer, has been hampered by imprecise and incomplete characterization of its heterogeneity and its attributes. Using single-cell RNA sequencing, we explored the landscape of thyroid cancer at single-cell resolution with 46,205 cells and delineated its dedifferentiation process and suppressive immune microenvironment. The developmental trajectory indicated that anaplastic thyroid cancer (ATC) cells were derived from a small subset of papillary thyroid cancer (PTC) cells. Moreover, a potential functional role of CREB3L1 on ATC development was revealed by integrated analyses of copy number alteration and transcriptional regulatory network. Multiple genes in differentiation-related pathways (e.g., EMT) were involved as the downstream targets of CREB3L1, increased expression of which can thus predict higher relapse risk of PTC. Collectively, our study provided insights into the heterogeneity and molecular evolution of thyroid cancer and highlighted the potential driver role of CREB3L1 in its dedifferentiation process.

Genome-scale screens identify factors regulating tumor cell responses to natural killer cells.

To systematically define molecular features in human tumor cells that determine their degree of sensitivity to human allogeneic natural killer (NK) cells, we quantified the NK cell responsiveness of hundreds of molecularly annotated 'DNA-barcoded' solid tumor cell lines in multiplexed format and applied genome-scale CRISPR-based gene-editing screens in several solid tumor cell lines, to functionally interrogate which genes in tumor cells regulate the response to NK cells. In these orthogonal studies, NK cell-sensitive tumor cells tend to exhibit 'mesenchymal-like' transcriptional programs; high transcriptional signature for chromatin remodeling complexes; high levels of B7-H6 (NCR3LG1); and low levels of HLA-E/antigen presentation genes. Importantly, transcriptional signatures of NK cell-sensitive tumor cells correlate with immune checkpoint inhibitor (ICI) resistance in clinical samples. This study provides a comprehensive map of mechanisms regulating tumor cell responses to NK cells, with implications for future biomarker-driven applications of NK cell immunotherapies.

Functional Genomics Identify Distinct and Overlapping Genes Mediating Resistance to Different Classes of Heterobifunctional Degraders of Oncoproteins.

Heterobifunctional proteolysis-targeting chimeric compounds leverage the activity of E3 ligases to induce degradation of target oncoproteins and exhibit potent preclinical antitumor activity. To dissect the mechanisms regulating tumor cell sensitivity to different classes of pharmacological "degraders" of oncoproteins, we performed genome-scale CRISPR-Cas9-based gene editing studies. We observed that myeloma cell resistance to degraders of different targets (BET bromodomain proteins, CDK9) and operating through CRBN (degronimids) or VHL is primarily mediated by prevention of, rather than adaptation to, breakdown of the target oncoprotein; and this involves loss of function of the cognate E3 ligase or interactors/regulators of the respective cullin-RING ligase (CRL) complex. The substantial gene-level differences for resistance mechanisms to CRBN- versus VHL-based degraders explains mechanistically the lack of cross-resistance with sequential administration of these two degrader classes. Development of degraders leveraging more diverse E3 ligases/CRLs may facilitate sequential/alternating versus combined uses of these agents toward potentially delaying or preventing resistance.

Novel Recurrent Altered Genes in Chinese Patients With Anaplastic Thyroid Cancer.

BACKGROUND: Anaplastic thyroid cancer (ATC) is a rare but lethal malignancy, and few systematic investigations on genomic profiles of ATC have been performed in Chinese patients. METHODS: Fifty-four ATC patients in West China Hospital between 2010 to 2020 were retrospectively analyzed, while 29 patients with available samples were sequenced by whole-exome sequencing (WES). The associations between genomic alterations and clinical characteristics were statistically evaluated. RESULTS: The median overall survival was 3.0 months in the entire cohort, which was impacted by multiple clinical features, including age, tumor size, and different treatment strategies. In the WES cohort, totally 797 nonsilent mutations were detected; the most frequently altered genes were TP53 (48%), BRAF (24%), PIK3CA (24%), and TERT promoter (21%). Although these mutations have been well-reported in previous studies, ethnic specificity was exhibited in terms of mutation frequency. Moreover, several novel significantly mutated genes were identified including RBM15 (17%), NOTCH2NL (14%), CTNNA3 (10%), and KATNAL2 (10%). WES-based copy number alteration analysis also revealed a high frequent gain of NOTCH2NL (41%), which induced its increased expression. Gene mutations and copy number alterations were enriched in phosphatidylinositol 3-kinase/AKT/mechanistic target of rapamycin (mTOR), NOTCH, and WNT pathways. CONCLUSIONS: This study reveals shared and ethnicity-specific genomic profiles of ATC in Chinese patients and suggests NOTCH2NL may act as a novel candidate driver gene for ATC tumorigenesis.

A deep learning diagnostic platform for diffuse large B-cell lymphoma with high accuracy across multiple hospitals.

Diagnostic histopathology is a gold standard for diagnosing hematopoietic malignancies. Pathologic diagnosis requires labor-intensive reading of a large number of tissue slides with high diagnostic accuracy equal or close to 100 percent to guide treatment options, but this requirement is difficult to meet. Although artificial intelligence (AI) helps to reduce the labor of reading pathologic slides, diagnostic accuracy has not reached a clinically usable level. Establishment of an AI model often demands big datasets and an ability to handle large variations in sample preparation and image collection. Here, we establish a highly accurate deep learning platform, consisting of multiple convolutional neural networks, to classify pathologic images by using smaller datasets. We analyze human diffuse large B-cell lymphoma (DLBCL) and non-DLBCL pathologic images from three hospitals separately using AI models, and obtain a diagnostic rate of close to 100 percent (100% for hospital A, 99.71% for hospital B and 100% for hospital C). The technical variability introduced by slide preparation and image collection reduces AI model performance in cross-hospital tests, but the 100% diagnostic accuracy is maintained after its elimination. It is now clinically practical to utilize deep learning models for diagnosis of DLBCL and ultimately other human hematopoietic malignancies.

Targeting glutaminase1 and synergizing with clinical drugs achieved more promising antitumor activity on multiple myeloma.

Multiple myeloma (MM) pathogenesis remains incompletely understood and biomarkers predicting treatment response still remain lacking. Here we describe the rational mechanisms of combining targeting glautaminase1 (GLS1) with other chemo-reagents for MM treatment. Gls1 is highly expressed cMYC/KRAS12V-drived plasmacytoma (PCT) cells. Down-regulation of Gls1 with miRNAi in cMYC/KRAS12V-expressing BaF3 cells prevented them from growing independence of interleukin 3 (IL3). By using our cMYC/KRAS12V-transduced adoptive plasmacytoma mouse model, we found that Gls1 is involved in PCT pathogenesis. Down-regulation of Gls1 significantly prolonged the survival of PCT recipients. Knockdown of Gls1 increased the expression of Cdkn1a and Cdkn1b and decreased the expression of some critical oncogenes for cancer cell survival, such as c-Myc, Cdk4, and NfκB, as well as some genes which are essential for MM cell survival, such as Irf4, Prdm1, Csnk1α1, and Rassf5. Combination of Gls1 inhibition with LBH589, Bortezomib, or Lenalidomide significantly impaired tumor growth in a MM xenograft mouse model. Our data strongly suggest that Gls1 plays an important role for MM pathogenesis and that combination of GLS1 inhibitor with other MM therapy agents could benefit to MM patients.