consICA: an R package for robust reference-free deconvolution of multi-omics data.

Motivation: Deciphering molecular signals from omics data helps understanding cellular processes and disease progression. Effective algorithms for extracting these signals are essential, with a strong emphasis on robustness and reproducibility. Results: R/Bioconductor package consICA implements consensus independent component analysis (ICA)-a data-driven deconvolution method to decompose heterogeneous omics data and extract features suitable for patient stratification and multimodal data integration. The method separates biologically relevant molecular signals from technical effects and provides information about the cellular composition and biological processes. Build-in annotation, survival analysis, and report generation provide useful tools for the interpretation of extracted signals. The implementation of parallel computing in the package ensures efficient analysis using modern multicore systems. The package offers a reproducible and efficient data-driven solution for the analysis of complex molecular profiles, with significant implications for cancer research. Availability and implementation: The package is implemented in R and available under MIT license at Bioconductor (https://bioconductor.org/packages/consICA) or at GitHub (https://github.com/biomod-lih/consICA).

PARK7/DJ-1 deficiency impairs microglial activation in response to LPS-induced inflammation.

BACKGROUND: Specific microglia responses are thought to contribute to the development and progression of neurodegenerative diseases, including Parkinson's disease (PD). However, the phenotypic acquisition of microglial cells and their role during the underlying neuroinflammatory processes remain largely elusive. Here, according to the multiple-hit hypothesis, which stipulates that PD etiology is determined by a combination of genetics and various environmental risk factors, we investigate microglial transcriptional programs and morphological adaptations under PARK7/DJ-1 deficiency, a genetic cause of PD, during lipopolysaccharide (LPS)-induced inflammation. METHODS: Using a combination of single-cell RNA-sequencing, bulk RNA-sequencing, multicolor flow cytometry and immunofluorescence analyses, we comprehensively compared microglial cell phenotypic characteristics in PARK7/DJ-1 knock-out (KO) with wildtype littermate mice following 6- or 24-h intraperitoneal injection with LPS. For translational perspectives, we conducted corresponding analyses in human PARK7/DJ-1 mutant induced pluripotent stem cell (iPSC)-derived microglia and murine bone marrow-derived macrophages (BMDMs). RESULTS: By excluding the contribution of other immune brain resident and peripheral cells, we show that microglia acutely isolated from PARK7/DJ-1 KO mice display a distinct phenotype, specially related to type II interferon and DNA damage response signaling, when compared with wildtype microglia, in response to LPS. We also detected discrete signatures in human PARK7/DJ-1 mutant iPSC-derived microglia and BMDMs from PARK7/DJ-1 KO mice. These specific transcriptional signatures were reflected at the morphological level, with microglia in LPS-treated PARK7/DJ-1 KO mice showing a less amoeboid cell shape compared to wildtype mice, both at 6 and 24 h after acute inflammation, as also observed in BMDMs. CONCLUSIONS: Taken together, our results show that, under inflammatory conditions, PARK7/DJ-1 deficiency skews microglia towards a distinct phenotype characterized by downregulation of genes involved in type II interferon signaling and a less prominent amoeboid morphology compared to wildtype microglia. These findings suggest that the underlying oxidative stress associated with the lack of PARK7/DJ-1 affects microglia neuroinflammatory responses, which may play a causative role in PD onset and progression.

Glioblastoma-instructed microglia transition to heterogeneous phenotypic states with phagocytic and dendritic cell-like features in patient tumors and patient-derived orthotopic xenografts.

BACKGROUND: A major contributing factor to glioblastoma (GBM) development and progression is its ability to evade the immune system by creating an immune-suppressive environment, where GBM-associated myeloid cells, including resident microglia and peripheral monocyte-derived macrophages, play critical pro-tumoral roles. However, it is unclear whether recruited myeloid cells are phenotypically and functionally identical in GBM patients and whether this heterogeneity is recapitulated in patient-derived orthotopic xenografts (PDOXs). A thorough understanding of the GBM ecosystem and its recapitulation in preclinical models is currently missing, leading to inaccurate results and failures of clinical trials. METHODS: Here, we report systematic characterization of the tumor microenvironment (TME) in GBM PDOXs and patient tumors at the single-cell and spatial levels. We applied single-cell RNA sequencing, spatial transcriptomics, multicolor flow cytometry, immunohistochemistry, and functional studies to examine the heterogeneous TME instructed by GBM cells. GBM PDOXs representing different tumor phenotypes were compared to glioma mouse GL261 syngeneic model and patient tumors. RESULTS: We show that GBM tumor cells reciprocally interact with host cells to create a GBM patient-specific TME in PDOXs. We detected the most prominent transcriptomic adaptations in myeloid cells, with brain-resident microglia representing the main population in the cellular tumor, while peripheral-derived myeloid cells infiltrated the brain at sites of blood-brain barrier disruption. More specifically, we show that GBM-educated microglia undergo transition to diverse phenotypic states across distinct GBM landscapes and tumor niches. GBM-educated microglia subsets display phagocytic and dendritic cell-like gene expression programs. Additionally, we found novel microglial states expressing cell cycle programs, astrocytic or endothelial markers. Lastly, we show that temozolomide treatment leads to transcriptomic plasticity and altered crosstalk between GBM tumor cells and adjacent TME components. CONCLUSIONS: Our data provide novel insights into the phenotypic adaptation of the heterogeneous TME instructed by GBM tumors. We show the key role of microglial phenotypic states in supporting GBM tumor growth and response to treatment. Our data place PDOXs as relevant models to assess the functionality of the TME and changes in the GBM ecosystem upon treatment.

Early-life influenza A (H1N1) infection independently programs brain connectivity, HPA AXIS and tissue-specific gene expression profiles.

Early-life adversity covers a range of physical, social and environmental stressors. Acute viral infections in early life are a major source of such adversity and have been associated with a broad spectrum of later-life effects outside the immune system or "off-target". These include an altered hypothalamus-pituitary-adrenal (HPA) axis and metabolic reactions. Here, we used a murine post-natal day 14 (PND 14) Influenza A (H1N1) infection model and applied a semi-holistic approach including phenotypic measurements, gene expression arrays and diffusion neuroimaging techniques to investigate HPA axis dysregulation, energy metabolism and brain connectivity. By PND 56 the H1N1 infection had been resolved, and there was no residual gene expression signature of immune cell infiltration into the liver, adrenal gland or brain tissues examined nor of immune-related signalling. A resolved early-life H1N1 infection had sex-specific effects. We observed retarded growth of males and altered pre-stress (baseline) blood glucose and corticosterone levels at PND42 after the infection was resolved. Cerebral MRI scans identified reduced connectivity in the cortex, midbrain and cerebellum that were accompanied by tissue-specific gene expression signatures. Gene set enrichment analysis confirmed that these were tissue-specific changes with few common pathways. Early-life infection independently affected each of the systems and this was independent of HPA axis or immune perturbations.

Cross-Talk between miRNAs from the Dlk1-Dio3 Locus and Histone Methylation to Protect Male Cerebellum from Methyl Donor Deficiency.

SCOPE: Disruption of the one carbon metabolism during development, i.e., following a gestational vitamin B9 and B12 deficiencies, is involved in birth defects and brain development delay. Using a rat nutritional model, consisting of pups born to dams fed a vitamin B9 and B12 deficient diet (MDD), the study previously reports molecular and cellular alterations in the brain, in a sex dependent manner, with females being more affected than males. The study hypothesizes that epigenetic modifications could participate in the sex differences is observed. METHODS AND RESULTS: The study investigates lysine methylation of histones and expression of microRNAs in the cerebellum of MDD male and female pups. The study reports a differential regulation of H3K36Me2 and H4K20Me3 between males and females, in response to MDD. Moreover, distinct regulation of Kmt5b and Kdm2a expression by miR-134-5p and miR-369-5p from the Dlk1-Dio3 locus, contributes to the maintenance of expression of genes involved in synaptic plasticity. CONCLUSION: These results could explain the neuroprotection to MDD that male pups display. The work will contribute to the understanding of the consequences of vitamin starvation on brain development, as well as how the epigenome is affected by one carbon metabolism disruption.

DeSUMOylation of chromatin-bound proteins limits the rapid transcriptional reprogramming induced by daunorubicin in acute myeloid leukemias.

Genotoxicants have been used for decades as front-line therapies against cancer on the basis of their DNA-damaging actions. However, some of their non-DNA-damaging effects are also instrumental for killing dividing cells. We report here that the anthracycline Daunorubicin (DNR), one of the main drugs used to treat Acute Myeloid Leukemia (AML), induces rapid (3 h) and broad transcriptional changes in AML cells. The regulated genes are particularly enriched in genes controlling cell proliferation and death, as well as inflammation and immunity. These transcriptional changes are preceded by DNR-dependent deSUMOylation of chromatin proteins, in particular at active promoters and enhancers. Surprisingly, inhibition of SUMOylation with ML-792 (SUMO E1 inhibitor), dampens DNR-induced transcriptional reprogramming. Quantitative proteomics shows that the proteins deSUMOylated in response to DNR are mostly transcription factors, transcriptional co-regulators and chromatin organizers. Among them, the CCCTC-binding factor CTCF is highly enriched at SUMO-binding sites found in cis-regulatory regions. This is notably the case at the promoter of the DNR-induced NFKB2 gene. DNR leads to a reconfiguration of chromatin loops engaging CTCF- and SUMO-bound NFKB2 promoter with a distal cis-regulatory region and inhibition of SUMOylation with ML-792 prevents these changes.

Impact of IDH Mutations, the 1p/19q Co-Deletion and the G-CIMP Status on Alternative Splicing in Diffuse Gliomas.

By generating protein diversity, alternative splicing provides an important oncogenic pathway. Isocitrate dehydrogenase (IDH) 1 and 2 mutations and 1p/19q co-deletion have become crucial for the novel molecular classification of diffuse gliomas, which also incorporates DNA methylation profiling. In this study, we have carried out a bioinformatics analysis to examine the impact of the IDH mutation, as well as the 1p/19q co-deletion and the glioma CpG island methylator phenotype (G-CIMP) status on alternative splicing in a cohort of 662 diffuse gliomas from The Cancer Genome Atlas (TCGA). We identify the biological processes and molecular functions affected by alternative splicing in the various glioma subgroups and provide evidence supporting the important contribution of alternative splicing in modulating epigenetic regulation in diffuse gliomas. Targeting the genes and pathways affected by alternative splicing might provide novel therapeutic opportunities against gliomas.

Actin cytoskeleton depolymerization increases matrix metalloproteinase gene expression in breast cancer cells by promoting translocation of cysteine-rich protein 2 to the nucleus.

The actin cytoskeleton plays a critical role in cancer cell invasion and metastasis; however, the coordination of its multiple functions remains unclear. Actin dynamics in the cytoplasm control the formation of invadopodia, which are membrane protrusions that facilitate cancer cell invasion by focusing the secretion of extracellular matrix-degrading enzymes, including matrix metalloproteinases (MMPs). In this study, we investigated the nuclear role of cysteine-rich protein 2 (CRP2), a two LIM domain-containing F-actin-binding protein that we previously identified as a cytoskeletal component of invadopodia, in breast cancer cells. We found that F-actin depolymerization stimulates the translocation of CRP2 into the nucleus, resulting in an increase in the transcript levels of pro-invasive and pro-metastatic genes, including several members of the MMP gene family. We demonstrate that in the nucleus, CRP2 interacts with the transcription factor serum response factor (SRF), which is crucial for the expression of MMP-9 and MMP-13. Our data suggest that CRP2 and SRF cooperate to modulate of MMP expression levels. Furthermore, Kaplan-Meier analysis revealed a significant association between high-level expression of SRF and shorter overall survival and distant metastasis-free survival in breast cancer patients with a high CRP2 expression profile. Our findings suggest a model in which CRP2 mediates the coordination of cytoplasmic and nuclear processes driven by actin dynamics, ultimately resulting in the induction of invasive and metastatic behavior in breast cancer cells.

Glioblastoma-instructed microglia transition to heterogeneous phenotypic states with phagocytic and dendritic cell-like features in patient tumors and patient-derived orthotopic xenografts.

Background: A major contributing factor to glioblastoma (GBM) development and progression is its ability to evade the immune system by creating an immune-suppressive environment, where GBM-associated myeloid cells, including resident microglia and peripheral monocyte-derived macrophages, play critical pro-tumoral roles. However, it is unclear whether recruited myeloid cells are phenotypically and functionally identical in GBM patients and whether this heterogeneity is recapitulated in patient-derived orthotopic xenografts (PDOXs). A thorough understanding of the GBM ecosystem and its recapitulation in preclinical models is currently missing, leading to inaccurate results and failures of clinical trials. Methods: Here, we report systematic characterization of the tumor microenvironment (TME) in GBM PDOXs and patient tumors at the single-cell and spatial levels. We applied single-cell RNA-sequencing, spatial transcriptomics, multicolor flow cytometry, immunohistochemistry and functional studies to examine the heterogeneous TME instructed by GBM cells. GBM PDOXs representing different tumor phenotypes were compared to glioma mouse GL261 syngeneic model and patient tumors. Results: We show that GBM tumor cells reciprocally interact with host cells to create a GBM patient-specific TME in PDOXs. We detected the most prominent transcriptomic adaptations in myeloid cells, with brain-resident microglia representing the main population in the cellular tumor, while peripheral-derived myeloid cells infiltrated the brain at sites of blood-brain barrier disruption. More specifically, we show that GBM-educated microglia undergo transition to diverse phenotypic states across distinct GBM landscapes and tumor niches. GBM-educated microglia subsets display phagocytic and dendritic cell-like gene expression programs. Additionally, we found novel microglial states expressing cell cycle programs, astrocytic or endothelial markers. Lastly, we show that temozolomide treatment leads to transcriptomic plasticity and altered crosstalk between GBM tumor cells and adjacent TME components. Conclusions: Our data provide novel insights into the phenotypic adaptation of the heterogeneous TME instructed by GBM tumors. We show the key role of microglial phenotypic states in supporting GBM tumor growth and response to treatment. Our data place PDOXs as relevant models to assess the functionality of the TME and changes in the GBM ecosystem upon treatment.

Neurodegeneration and neuroinflammation are linked, but independent of alpha-synuclein inclusions, in a seeding/spreading mouse model of Parkinson's disease.

A key pathological process in Parkinson's disease (PD) is the transneuronal spreading of α-synuclein. Alpha-synuclein (α-syn) is a presynaptic protein that, in PD, forms pathological inclusions. Other hallmarks of PD include neurodegeneration and microgliosis in susceptible brain regions. Whether it is primarily transneuronal spreading of α-syn particles, inclusion formation, or other mechanisms, such as inflammation, that cause neurodegeneration in PD is unclear. We used a model of spreading of α-syn induced by striatal injection of α-syn preformed fibrils into the mouse striatum to address this question. We performed quantitative analysis for α-syn inclusions, neurodegeneration, and microgliosis in different brain regions, and generated gene expression profiles of the ventral midbrain, at two different timepoints after disease induction. We observed significant neurodegeneration and microgliosis in brain regions not only with, but also without α-syn inclusions. We also observed prominent microgliosis in injured brain regions that did not correlate with neurodegeneration nor with inclusion load. Using longitudinal gene expression profiling, we observed early gene expression changes, linked to neuroinflammation, that preceded neurodegeneration, indicating an active role of microglia in this process. Altered gene pathways overlapped with those typical of PD. Our observations indicate that α-syn inclusion formation is not the major driver in the early phases of PD-like neurodegeneration, but that microglia, activated by diffusible, oligomeric α-syn, may play a key role in this process. Our findings uncover new features of α-syn induced pathologies, in particular microgliosis, and point to the necessity for a broader view of the process of α-syn spreading.

Mitochondrial metabolism supports resistance to IDH mutant inhibitors in acute myeloid leukemia.

Mutations in IDH induce epigenetic and transcriptional reprogramming, differentiation bias, and susceptibility to mitochondrial inhibitors in cancer cells. Here, we first show that cell lines, PDXs, and patients with acute myeloid leukemia (AML) harboring an IDH mutation displayed an enhanced mitochondrial oxidative metabolism. Along with an increase in TCA cycle intermediates, this AML-specific metabolic behavior mechanistically occurred through the increase in electron transport chain complex I activity, mitochondrial respiration, and methylation-driven CEBPα-induced fatty acid ß-oxidation of IDH1 mutant cells. While IDH1 mutant inhibitor reduced 2-HG oncometabolite and CEBPα methylation, it failed to reverse FAO and OxPHOS. These mitochondrial activities were maintained through the inhibition of Akt and enhanced activation of peroxisome proliferator-activated receptor-γ coactivator-1 PGC1α upon IDH1 mutant inhibitor. Accordingly, OxPHOS inhibitors improved anti-AML efficacy of IDH mutant inhibitors in vivo. This work provides a scientific rationale for combinatory mitochondrial-targeted therapies to treat IDH mutant AML patients, especially those unresponsive to or relapsing from IDH mutant inhibitors.

Fra-1 regulates its target genes via binding to remote enhancers without exerting major control on chromatin architecture in triple negative breast cancers.

The ubiquitous family of dimeric transcription factors AP-1 is made up of Fos and Jun family proteins. It has long been thought to operate principally at gene promoters and how it controls transcription is still ill-understood. The Fos family protein Fra-1 is overexpressed in triple negative breast cancers (TNBCs) where it contributes to tumor aggressiveness. To address its transcriptional actions in TNBCs, we combined transcriptomics, ChIP-seqs, machine learning and NG Capture-C. Additionally, we studied its Fos family kin Fra-2 also expressed in TNBCs, albeit much less. Consistently with their pleiotropic effects, Fra-1 and Fra-2 up- and downregulate individually, together or redundantly many genes associated with a wide range of biological processes. Target gene regulation is principally due to binding of Fra-1 and Fra-2 at regulatory elements located distantly from cognate promoters where Fra-1 modulates the recruitment of the transcriptional co-regulator p300/CBP and where differences in AP-1 variant motif recognition can underlie preferential Fra-1- or Fra-2 bindings. Our work also shows no major role for Fra-1 in chromatin architecture control at target gene loci, but suggests collaboration between Fra-1-bound and -unbound enhancers within chromatin hubs sometimes including promoters for other Fra-1-regulated genes. Our work impacts our view of AP-1.

Adrenomedullin-CALCRL axis controls relapse-initiating drug tolerant acute myeloid leukemia cells.

Drug tolerant/resistant leukemic stem cell (LSC) subpopulations may explain frequent relapses in acute myeloid leukemia (AML), suggesting that these relapse-initiating cells (RICs) persistent after chemotherapy represent bona fide targets to prevent drug resistance and relapse. We uncover that calcitonin receptor-like receptor (CALCRL) is expressed in RICs, and that the overexpression of CALCRL and/or of its ligand adrenomedullin (ADM), and not CGRP, correlates to adverse outcome in AML. CALCRL knockdown impairs leukemic growth, decreases LSC frequency, and sensitizes to cytarabine in patient-derived xenograft models. Mechanistically, the ADM-CALCRL axis drives cell cycle, DNA repair, and mitochondrial OxPHOS function of AML blasts dependent on E2F1 and BCL2. Finally, CALCRL depletion reduces LSC frequency of RICs post-chemotherapy in vivo. In summary, our data highlight a critical role of ADM-CALCRL in post-chemotherapy persistence of these cells, and disclose a promising therapeutic target to prevent relapse in AML.

Mitochondrial inhibitors circumvent adaptive resistance to venetoclax and cytarabine combination therapy in acute myeloid leukemia.

Therapy resistance represents a major clinical challenge in acute myeloid leukemia (AML). Here we define a 'MitoScore' signature, which identifies high mitochondrial oxidative phosphorylation in vivo and in patients with AML. Primary AML cells with cytarabine (AraC) resistance and a high MitoScore relied on mitochondrial Bcl2 and were highly sensitive to venetoclax (VEN) + AraC (but not to VEN + azacytidine). Single-cell transcriptomics of VEN + AraC-residual cell populations revealed adaptive resistance associated with changes in oxidative phosphorylation, electron transport chain complex and the TP53 pathway. Accordingly, treatment of VEN + AraC-resistant AML cells with electron transport chain complex inhibitors, pyruvate dehydrogenase inhibitors or mitochondrial ClpP protease agonists substantially delayed relapse following VEN + AraC. These findings highlight the central role of mitochondrial adaptation during AML therapy and provide a scientific rationale for alternating VEN + azacytidine with VEN + AraC in patients with a high MitoScore and to target mitochondrial metabolism to enhance the sensitivity of AML cells to currently approved therapies.

A New Synuclein-Transgenic Mouse Model for Early Parkinson's Reveals Molecular Features of Preclinical Disease.

Understanding Parkinson's disease (PD), in particular in its earliest phases, is important for diagnosis and treatment. However, human brain samples are collected post-mortem, reflecting mainly end-stage disease. Because brain samples of mouse models can be collected at any stage of the disease process, they are useful in investigating PD progression. Here, we compare ventral midbrain transcriptomics profiles from α-synuclein transgenic mice with a progressive, early PD-like striatal neurodegeneration across different ages using pathway, gene set, and network analysis methods. Our study uncovers statistically significant altered genes across ages and between genotypes with known, suspected, or unknown function in PD pathogenesis and key pathways associated with disease progression. Among those are genotype-dependent alterations associated with synaptic plasticity and neurotransmission, as well as mitochondria-related genes and dysregulation of lipid metabolism. Age-dependent changes were among others observed in neuronal and synaptic activity, calcium homeostasis, and membrane receptor signaling pathways, many of which linked to G-protein coupled receptors. Most importantly, most changes occurred before neurodegeneration was detected in this model, which points to a sequence of gene expression events that may be relevant for disease initiation and progression. It is tempting to speculate that molecular changes similar to those changes observed in our model happen in midbrain dopaminergic neurons before they start to degenerate. In other words, we believe we have uncovered molecular changes that accompany the progression from preclinical to early PD.

Patient-derived organoids and orthotopic xenografts of primary and recurrent gliomas represent relevant patient avatars for precision oncology.

Patient-based cancer models are essential tools for studying tumor biology and for the assessment of drug responses in a translational context. We report the establishment a large cohort of unique organoids and patient-derived orthotopic xenografts (PDOX) of various glioma subtypes, including gliomas with mutations in IDH1, and paired longitudinal PDOX from primary and recurrent tumors of the same patient. We show that glioma PDOXs enable long-term propagation of patient tumors and represent clinically relevant patient avatars that retain histopathological, genetic, epigenetic, and transcriptomic features of parental tumors. We find no evidence of mouse-specific clonal evolution in glioma PDOXs. Our cohort captures individual molecular genotypes for precision medicine including mutations in IDH1, ATRX, TP53, MDM2/4, amplification of EGFR, PDGFRA, MET, CDK4/6, MDM2/4, and deletion of CDKN2A/B, PTCH, and PTEN. Matched longitudinal PDOX recapitulate the limited genetic evolution of gliomas observed in patients following treatment. At the histological level, we observe increased vascularization in the rat host as compared to mice. PDOX-derived standardized glioma organoids are amenable to high-throughput drug screens that can be validated in mice. We show clinically relevant responses to temozolomide (TMZ) and to targeted treatments, such as EGFR and CDK4/6 inhibitors in (epi)genetically defined subgroups, according to MGMT promoter and EGFR/CDK status, respectively. Dianhydrogalactitol (VAL-083), a promising bifunctional alkylating agent in the current clinical trial, displayed high therapeutic efficacy, and was able to overcome TMZ resistance in glioblastoma. Our work underscores the clinical relevance of glioma organoids and PDOX models for translational research and personalized treatment studies and represents a unique publicly available resource for precision oncology.

Reference-free deconvolution, visualization and interpretation of complex DNA methylation data using DecompPipeline, MeDeCom and FactorViz.

DNA methylation profiling offers unique insights into human development and diseases. Often the analysis of complex tissues and cell mixtures is the only feasible option to study methylation changes across large patient cohorts. Since DNA methylomes are highly cell type specific, deconvolution methods can be used to recover cell type-specific information in the form of latent methylation components (LMCs) from such 'bulk' samples. Reference-free deconvolution methods retrieve these components without the need for DNA methylation profiles of purified cell types. Currently no integrated and guided procedure is available for data preparation and subsequent interpretation of deconvolution results. Here, we describe a three-stage protocol for reference-free deconvolution of DNA methylation data comprising: (i) data preprocessing, confounder adjustment using independent component analysis (ICA) and feature selection using DecompPipeline, (ii) deconvolution with multiple parameters using MeDeCom, RefFreeCellMix or EDec and (iii) guided biological inference and validation of deconvolution results with the R/Shiny graphical user interface FactorViz. Our protocol simplifies the analysis and guides the initial interpretation of DNA methylation data derived from complex samples. The harmonized approach is particularly useful to dissect and evaluate cell heterogeneity in complex systems such as tumors. We apply the protocol to lung cancer methylomes from The Cancer Genome Atlas (TCGA) and show that our approach identifies the proportions of stromal cells and tumor-infiltrating immune cells, as well as associations of the detected components with clinical parameters. The protocol takes slightly >3 d to complete and requires basic R skills.

Extracellular ATP and CD39 Activate cAMP-Mediated Mitochondrial Stress Response to Promote Cytarabine Resistance in Acute Myeloid Leukemia.

Relapses driven by chemoresistant leukemic cell populations are the main cause of mortality for patients with acute myeloid leukemia (AML). Here, we show that the ectonucleotidase CD39 (ENTPD1) is upregulated in cytarabine-resistant leukemic cells from both AML cell lines and patient samples in vivo and in vitro. CD39 cell-surface expression and activity is increased in patients with AML upon chemotherapy compared with diagnosis, and enrichment in CD39-expressing blasts is a marker of adverse prognosis in the clinics. High CD39 activity promotes cytarabine resistance by enhancing mitochondrial activity and biogenesis through activation of a cAMP-mediated adaptive mitochondrial stress response. Finally, genetic and pharmacologic inhibition of CD39 ecto-ATPase activity blocks the mitochondrial reprogramming triggered by cytarabine treatment and markedly enhances its cytotoxicity in AML cells in vitro and in vivo. Together, these results reveal CD39 as a new residual disease marker and a promising therapeutic target to improve chemotherapy response in AML. SIGNIFICANCE: Extracellular ATP and CD39-P2RY13-cAMP-OxPHOS axis are key regulators of cytarabine resistance, offering a new promising therapeutic strategy in AML.This article is highlighted in the In This Issue feature, p. 1426.

Machine learning-assisted neurotoxicity prediction in human midbrain organoids.

INTRODUCTION: Brain organoids are highly complex multi-cellular tissue proxies, which have recently risen as novel tools to study neurodegenerative diseases such as Parkinson's disease (PD). However, with increasing complexity of the system, usage of quantitative tools becomes challenging. OBJECTIVES: The primary objective of this study was to develop a neurotoxin-induced PD organoid model and to assess the neurotoxic effect on dopaminergic neurons using microscopy-based phenotyping in a high-content fashion. METHODS: We describe a pipeline for a machine learning-based analytical method, allowing for detailed image-based cell profiling and toxicity prediction in brain organoids treated with the neurotoxic compound 6-hydroxydopamine (6-OHDA). RESULTS: We quantified features such as dopaminergic neuron count and neuronal complexity and built a machine learning classifier with the data to optimize data processing strategies and to discriminate between different treatment conditions. We validated the approach with high content imaging data from PD patient derived midbrain organoids. CONCLUSIONS: The here described model is a valuable tool for advanced in vitro PD modeling and to test putative neurotoxic compounds.

Oncogenic KIT mutations induce STAT3-dependent autophagy to support cell proliferation in acute myeloid leukemia.

Autophagy is associated with both survival and cell death in myeloid malignancies. Therefore, deciphering its role in different genetically defined subtypes of acute myeloid leukemia (AML) is critical. Activating mutations of the KIT receptor tyrosine kinase are frequently detected in core-binding factor AML and are associated with a greater risk of relapse. Herein, we report that basal autophagy was significantly increased by the KITD816V mutation in AML cells and contributed to support their cell proliferation and survival. Invalidation of the key autophagy protein Atg12 strongly reduced tumor burden and improved survival of immunocompromised NSG mice engrafted with KITD816V TF-1 cells. Downstream of KITD816V, STAT3, but not AKT or ERK pathways, was identified as a major regulator of autophagy. Accordingly, STAT3 pharmacological inhibition or downregulation inhibited autophagy and reduced tumor growth both in vitro and in vivo. Taken together, our results support the notion that targeting autophagy or STAT3 opens up an exploratory pathway for finding new therapeutic opportunities for patients with CBF-AML or others malignancies with KITD816V mutations.