Search

HIGH CD36 gene signature generated from upregulated genes in AML patients that express the most CD36 mRNA compared to that of the lowest expression of CD36 in TCGA cohort.

Supplementary Methods. Supplementary Figure S1. Clinical distributions of AML patients from TUH and 3 independent patient cohorts. Supplementary Figure S2. In vivo treatment with 60 mg/kg/d of cytarabine (AraC) given daily for 5 days induces a significant reduction of the total cell tumor burden in AML-engrafted mice. Supplementary Figure S3. In vivo cytarabine (AraC) treatment induces significant but heterogeneous response in bone marrow and spleen of AML-xenografted NSG mice. Supplementary Figure S4. Comparative analysis of the in vivo response to cytarabine (AraC) with clinicobiological data of AML patients. Supplementary Figure S5. In vivo cytarabine (AraC) treatment induces changes in CD34+/-CD38+/- phenotypes in AML-engrafted mice. Supplementary Figure S6. In vivo cytarabine (AraC) chemotherapy treatment leads to reduction of the absolute number of human CD34+CD38+/- populations in AML. Supplementary figure S7. No enrichment in G0 quiescent cells was observed in mice treated with sublethal dose of cytarabine (AraC) for 5 days in vivo. Supplementary figure S8. Mitochondrial and energetic features of LOW (KG1, KG1a, U937) and HIGH (MOLM14, MV4-11, HL60) OXPHOS AML cell lines. Supplementary figure S9. Functional analysis of the transcriptomes of LOW (KG1a, U937) versus HIGH (MOLM14, HL60) OXPHOS AML cell lines untreated or treated with metformin. Supplementary figure S10. AML cells surviving after cytarabine (AraC) treatment are resistant to chemotherapies and are pre-existing CD36+CD44+ phenotype with an increased oxidative metabolism. Supplementary figure S11. Culture in galactose induces energetic shift of LOW OXPHOS AML U937 cells toward HIGH OXPHOS state, leading to cytarabine (AraC) resistance in AML. Supplementary figure S12. Energetic shift of mtDNA-depleted Rho0 MOLM14 cells toward LOW OXPHOS state induces AraC sensitivity. Supplementary figure S13. Electron Transfer Chain Complex I inhibition by Phenformin (Phenf) induces energetic shift toward LOW OXPHOS state and increases cytarabine (AraC) sensitivity in MOLM14 cells. Supplementary figure S14. Electron Transfer Chain Complex I inhibition by Metformin (Met) induces energetic shift toward LOW OXPHOS state and increases cytarabine (AraC) sensitivity in MOLM14 cells. Supplementary figure S15. Electron Transfer Chain Complex I inhibition by Rotenone (Rot) induces energetic shift toward LOW OXPHOS state and increases cytarabine (AraC) sensitivity in MOLM14 cells. Supplementary figure S16. Electron Transfer Chain Complex III inhibition by Antimycin A (AntiA) induces energetic shift toward LOW OXPHOS state and increases cytarabine (AraC) sensitivity in MOLM14 cells. Supplementary figure S17. Electron Transfer Chain Complex III inhibition by Atovaquone (ATQ) induces energetic shift toward LOW OXPHOS state and increases cytarabine (AraC) sensitivity in MOLM14 cells. Supplementary figure S18. Working model of the resistance to AraC in vivo.

This document details the preliminary differential gene expression analysis carried out using the Köln dataset

This document contains Supplementary Table S1, which provides the expression data pertaining to the 27-gene signature and relevant clinical data of the Heidelberg cohort.

Inevitable tumor recurrence and a poor median survival are frustrating reminders of the inefficacy of our current standard of care for patients with newly diagnosed glioblastoma (GBM), which includes surgery followed by radiotherapy and chemotherapy with the DNA alkylating agent temozolomide. Because resistance to genotoxic damage is achieved mainly through execution of the DNA damage response (DDR) and DNA repair pathways, knowledge of the changes in DNA repair and cell-cycle gene expression that occur during tumor development might help identify new targets and improve treatment. Here, we performed a gene expression analysis targeting components of the DNA repair and cell-cycle machineries in cohorts of paired tumor samples (i.e., biopsies from the same patient obtained at the time of primary tumor operation and at recurrence) from patients treated with radiotherapy or radiotherapy plus temozolomide. We identified and validated a 27-gene signature that resulted in the classification of GBM specimens into three groups, two of which displayed inverse expression profiles. Each group contained primary and recurrent samples, and the tumor at relapse frequently displayed a gene expression profile different from that of the matched primary biopsy. Within the groups that exhibited opposing gene expression profiles, the expression pattern of the gene signature at relapse was linked to progression-free survival. We provide experimental evidence that our signature exposes group-specific vulnerabilities against genotoxicants and inhibitors of the cell cycle and DDR, with the prospect of personalized therapeutic strategies.Significance: These findings suggest that classification of GBM tumors based on a DNA repair and cell-cycle gene expression signature exposes vulnerabilities to standard-of-care therapies and offers the potential for personalized therapeutic strategies.

This document contains Supplementary Figures S1-S6 as well as their Legends. Specifically, Fig. S1: characterization of the 52 gene signature; Fig. S2: Validation of the 27 DNA repair and cell cycle gene signature; Fig. S3: Analysis of the Murat cohort; Fig. S4: Analysis of the TCGA dataset; Fig. S5: Identification of differentially expressed genes and gene ontology analysis of the G1 and G3 groups; Fig. S6: DNA repair and cell-cycle-based strategies for personalized management of GBM".

This document contains Supplementary Table S4, which regroups the differentially expressed genes identified between the G1 and G3 groups of the Wang dataset and the gene ontology terms associated.

Small noncoding microRNAs (miRNA) regulate the expression of target mRNAs by repressing their translation or orchestrating their sequence-specific degradation. In this study, we investigated miRNA and miRNA target gene expression patterns in melanoma to identify candidate biomarkers for early and progressive disease. Because data presently available on miRNA expression in melanoma are inconsistent thus far, we applied several different miRNA detection and profiling techniques on a panel of 10 cell lines and 20 patient samples representing nevi and primary or metastatic melanoma. Expression of selected miRNAs was inconsistent when comparing cell line–derived and patient-derived data. Moreover, as expected, some discrepancies were also detected when miRNA microarray data were correlated with qPCR-measured expression levels. Nevertheless, we identified miRNA-200c to be consistently downregulated in melanocytes, melanoma cell lines, and patient samples, whereas miRNA-205 and miRNA-23b were markedly reduced only in patient samples. In contrast, miR-146a and miR-155 were upregulated in all analyzed patients but none of the cell lines. Whole-genome microarrays were performed for analysis of selected melanoma cell lines to identify potential transcriptionally regulated miRNA target genes. Using Ingenuity pathway analysis, we identified a deregulated gene network centered around microphthalmia-associated transcription factor, a transcription factor known to play a key role in melanoma development. Our findings define miRNAs and miRNA target genes that offer candidate biomarkers in human melanoma. Cancer Res; 70(10); 4163–73. ©2010 AACR.

Supplementary Figure 1 from Signatures of MicroRNAs and Selected MicroRNA Target Genes in Human Melanoma

Supplementary Figure 2 from Signatures of MicroRNAs and Selected MicroRNA Target Genes in Human Melanoma

Supplementary Table 1 from Signatures of MicroRNAs and Selected MicroRNA Target Genes in Human Melanoma

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 broad transcriptional changes in AML cells before cell death induction. 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, which limits both the positive and negative effects of DNR on transcription. Quantitative proteomics shows that proteins that are 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. Its induction is preceded by a SUMO-dependent reconfiguration of chromatin loops engaging its CTCF- and SUMO-bound promoter with distal cis-regulatory regions. Altogether, our work suggests that one of the earliest effects of DNR in AML cells is a SUMO-dependent transcriptional reprogramming.

Inevitable tumor recurrence and a poor median survival are frustrating reminders of the inefficacy of our current standard of care for patients with newly diagnosed glioblastoma (GBM), which includes surgery followed by radiotherapy and chemotherapy with the DNA alkylating agent temozolomide. Because resistance to genotoxic damage is achieved mainly through execution of the DNA damage response (DDR) and DNA repair pathways, knowledge of the changes in DNA repair and cell-cycle gene expression that occur during tumor development might help identify new targets and improve treatment. Here, we performed a gene expression analysis targeting components of the DNA repair and cell-cycle machineries in cohorts of paired tumor samples (i.e., biopsies from the same patient obtained at the time of primary tumor operation and at recurrence) from patients treated with radiotherapy or radiotherapy plus temozolomide. We identified and validated a 27-gene signature that resulted in the classification of GBM specimens into three groups, two of which displayed inverse expression profiles. Each group contained primary and recurrent samples, and the tumor at relapse frequently displayed a gene expression profile different from that of the matched primary biopsy. Within the groups that exhibited opposing gene expression profiles, the expression pattern of the gene signature at relapse was linked to progression-free survival. We provide experimental evidence that our signature exposes group-specific vulnerabilities against genotoxicants and inhibitors of the cell cycle and DDR, with the prospect of personalized therapeutic strategies. Significance: These findings suggest that classification of GBM tumors based on a DNA repair and cell-cycle gene expression signature exposes vulnerabilities to standard-of-care therapies and offers the potential for personalized therapeutic strategies.

Alternative splicing is an essential characteristic of living cells that usually infers a various exon-exon junction governed by different splice sites. The traditional classification based on the mode of use designates splice site to one of the two groups, constitutive or alternative. Here, we considered another criterion and reorganized splice sites into “unisplice” and “multisplice” groups according to the number of undertaken splicing events. This approach provided us with a new insight in the organization and functionality of leukemia cells. We determined features associated with uni- and multisplice sites and found that combinatorics of these sites follows strict rules of the power-law in the t(8;21)-positive leukemia cells. We also found that system splicing characteristics of the transcriptome of leukemia cells remained persistent after drastic changes in the transcript composition caused by knockdown of the RUNX1-RUNX1T1 oncogene. In this work, we show for the first time that leukemia cells possess a sub-set of unisplice sites with a hidden multisplice potential. These findings reveal a new side in organization and functioning of the leukemic cells and open up new perspectives in the study of the t(8;21)-positive leukemia.

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.

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. Electronic Supplementary Material The online version of this article (10.1007/s12035-020-02085-z) contains supplementary material, which is available to authorized users.

Interleukin-6 (IL-6)-type cytokines play important roles in liver (patho-)biology. For instance, they regulate the acute phase response to inflammatory signals and are involved in hepatocarcinogenesis. Much is known about the regulation of protein-coding genes by cytokines whereas their effects on the miRNome is less well understood. We performed a microarray screen to identify microRNAs (miRNAs) in human hepatocytes which are modulated by IL-6-type cytokines. Using samples of 2 donors, 27 and 68 miRNAs (out of 1,733) were found to be differentially expressed upon stimulation with hyper-IL-6 (HIL-6) for up to 72 h, with an overlap of 15 commonly regulated miRNAs. qPCR validation revealed that miR-146b-5p was also consistently up-regulated in hepatocytes derived from 2 other donors. Interestingly, miR-146b-5p (but not miR-146a-5p) was induced by IL-6-type cytokines (HIL-6 and OSM) in non-transformed liver-derived PH5CH8 and THLE2 cells and in Huh-7 hepatoma cells, but not in HepG2 or Hep3B hepatoma cells. We did not find evidence for a differential regulation of miR-146b-5p expression by promoter methylation, also when analyzing the TCGA data set on liver cancer samples. Inducible overexpression of miR-146b-5p in PH5CH8 cells followed by RNA-Seq analysis revealed effects on multiple mRNAs, including those encoding IRAK1 and TRAF6 crucial for Toll-like receptor signaling. Indeed, LPS-mediated signaling was attenuated upon overexpression of miR-146b-5p, suggesting a regulatory loop to modulate inflammatory signaling in hepatocytes. Further validation experiments suggest DNAJC6, MAGEE1, MPHOSPH6, PPP2R1B, SLC10A3, SNRNP27, and TIMM17B to be novel targets for miR-146b-5p (and miR-146a-5p). IL-6-type cytokines up-regulate miR-146b-5p expression in liver cells, an anti-inflammatory miRNA with multiple targets shared with miR-146a-5p.

HuR regulates cytoplasmic mRNA stability and translatability, and the HuR expression level has been shown to correlate with poor disease outcome in several cancer types; however, the prognostic value and potential pro-oncogenic properties of HuR in meningioma remain unclear. Thus, in the present study, we analysed 85 meningioma tissue samples to establish the relationship between HuR expression, tumour cell proliferation, and/or patient survival. In addition, we examined the anti-proliferative effects of HuR knockdown in two meningioma cell lines (IOMM-Lee and Ben-Men-1) and conducted transcriptome-wide analyses (IOMM-Lee cells) to elucidate the molecular consequences of HuR knockdown. The results of the present study showed HuR cytoplasmic expression to correlate positively with tumour grade (p = 1.2 × 10-8 ) and negatively with progression-free and overall survival (p = 0.01) time in human meningioma tissues. In vitro, siHuR-induced HuR knockdown was shown to reduce the growth of both Ben-Men-1 (p = 2 × 10-8 ) and IOMM-Lee (p = 4 × 10-9 ) cells. Transcriptome analyses revealed HuR knockdown in IOMM-Lee cells to deregulate the HIF1A signalling pathway (p = 1.5 × 10-6 ) and to up-regulate the expression of genes essential for the assembly of the cytoplasmic mRNA processing body, global genome nucleotide-excision repair, poly(A)-specific ribonuclease activity, the positive regulation of apoptosis and of cell cycle arrest, and the negative regulation of RNA splicing [p(FDR)

Glioblastoma (GBM), the most frequent and aggressive malignant primary brain tumor, is characterized by a highly invasive growth. In our previous study we showed that overexpression of Carboxypeptidase E (CPE) mitigated glioma cell migration. In the present study we aimed at deciphering the regulatory mechanisms of the secreted form of CPE (sCPE). By transcriptome analysis and inhibition of signaling pathways involved in the regulation of cell growth and motility, we discovered that overexpression of sCPE was accompanied by differential regulation of mRNAs connected to the motility-associated networks, among others FAK, PAK, Cdc42, integrin, STAT3 as well as TGF-β. Especially SLUG was downregulated in sCPE-overexpressing glioma cells, paralleled by reduced expression of matrix-metalloproteinases (MMP) and, in consequence, by decreased cell migration. Expression of SLUG was regulated by ERK since inhibition of ERK reverted sCPE-mediated SLUG downregulation and enhanced cell motility. In a mouse glioma model, overexpression of sCPE significantly prolonged survival. Our results implicate a novel role for sCPE that mainly affects the expression of motility-associated genes via several signal pathways.

Interleukin-27 (IL27) is a type-I-cytokine of the IL6/IL12 family predominantly secreted by activated macrophages and dendritic cells. In the liver, IL27 expression was observed to be upregulated in patients with hepatitis B, and sera of hepatocellular carcinoma (HCC) patients contain significantly elevated levels of IL27 compared to healthy controls or patients with hepatitis and/or liver cirrhosis. In this study, we show that IL27 induces STAT1 and STAT3 phosphorylation in 5 HCC lines and 3 different types of non-transformed liver cells. We were especially interested in the relevance of the IL27-induced STAT3 activation in liver cells. Thus, we compared the IL27 responses with those induced by IFNγ (STAT1-dominated response) or IL6-type cytokines (IL6, hyper-IL6 (hy-IL6) or OSM) (STAT3-dominated response) by microarray analysis and find that in HCC cells, IL27 induces an IFNγ-like, STAT1-dependent transcriptional response, but we do not find an effective STAT3-dependent response. Validation experiments corroborate the finding from the microarray evaluation. Interestingly, the availability of STAT1 seems critical in the shaping of the IL27 response, as the siRNA knock-down of STAT1 revealed the ability of IL27 to induce the acute-phase protein γ-fibrinogen, a typical IL6 family characteristic. Moreover, we describe a crosstalk between the signaling of IL6-type cytokines and IL27: responses to the gp130-engaging cytokine IL27 (but not those to IFNs) can be inhibited by IL6-type cytokine pre-stimulation, likely by a SOCS3-mediated mechanism. Thus, IL27 recapitulates IFNγ responses in liver cells, but differs from IFNγ by its sensitivity to SOCS3 inhibition.

Background Invasion and angiogenesis are major hallmarks of glioblastoma (GBM) growth. While invasive tumor cells grow adjacent to blood vessels in normal brain tissue, tumor cells within neovascularized regions exhibit hypoxic stress and promote angiogenesis. The distinct microenvironments likely differentially affect metabolic processes within the tumor cells. Methods In the present study, we analyzed gene expression and metabolic changes in a human GBM xenograft model that displayed invasive and angiogenic phenotypes. In addition, we used glioma patient biopsies to confirm the results from the xenograft model. Results We demonstrate that the angiogenic switch in our xenograft model is linked to a proneural-to-mesenchymal transition that is associated with upregulation of the transcription factors BHLHE40, CEBPB, and STAT3. Metabolic analyses revealed that angiogenic xenografts employed higher rates of glycolysis compared with invasive xenografts. Likewise, patient biopsies exhibited higher expression of the glycolytic enzyme lactate dehydrogenase A and glucose transporter 1 in hypoxic areas compared with the invasive edge and lower-grade tumors. Analysis of the mitochondrial respiratory chain showed reduction of complex I in angiogenic xenografts and hypoxic regions of GBM samples compared with invasive xenografts, nonhypoxic GBM regions, and lower-grade tumors. In vitro hypoxia experiments additionally revealed metabolic adaptation of invasive tumor cells, which increased lactate production under long-term hypoxia. Conclusions The use of glycolysis versus mitochondrial respiration for energy production within human GBM tumors is highly dependent on the specific microenvironment. The metabolic adaptability of GBM cells highlights the difficulty of targeting one specific metabolic pathway for effective therapeutic intervention.

Background and purpose Mineralocorticoid receptor (MR) activation contributes to heart failure (HF) progression. Its overactivity in obesity is thought to accelerate cardiac remodelling and HF development. Given that MR antagonists (MRA) are beneficial in chronic HF patients, we hypothesized that early MRA treatment may target obesity-related disorders and consequently delay the development of HF. Experimental approach Twenty spontaneously hypertensive HF dyslipidaemic obese SHHF(cp/cp) rats and 18 non-dyslipidaemic lean SHHF(+/+) controls underwent regular monitoring for their metabolic and cardiovascular phenotypes with or without MRA treatment [eplerenone (eple), 100 mg∙kg(-1) ∙day(-1) ] from 1.5 to 12.5 months of age. Key results Eleven months of eple treatment in obese rats (SHHF(cp/cp) eple) reduced the obesity-related metabolic disorders observed in untreated SHHF(cp/cp) rats by reducing weight gain, triglycerides and total cholesterol levels and by preserving adiponectinaemia. The MRA treatment predominantly preserved diastolic and systolic functions in obese rats by alleviating the eccentric cardiac hypertrophy observed in untreated SHHF(cp/cp) animals and preserving ejection fraction (70 ± 1 vs. 59 ± 1%). The MRA also improved survival independently of these pressure effects. Conclusion and implications Early chronic eple treatment resulted in a delay in cardiac remodelling and HF onset in both SHHF(+/+) and SHHF(cp/cp) rats, whereas SHHF(cp/cp) rats further benefited from the MRA treatment through a reduction in their obesity and dyslipidaemia. These findings suggest that preventive MRA therapy may provide greater benefits in obese patients with additional risk factors of developing cardiovascular complications.

Boutzen et al. show that the IDH1 mutation and its oncometabolite, (R)-2-hydroxyglutarate, dysregulate downstream target pathways of myeloid-specific TFs, especially CEBPα, priming mutant IDH1-R132H AML blasts to the granulomonocytic lineage., Acute myeloid leukemia (AML) is characterized by the accumulation of malignant blasts with impaired differentiation programs caused by recurrent mutations, such as the isocitrate dehydrogenase (IDH) mutations found in 15% of AML patients. These mutations result in the production of the oncometabolite (R)-2-hydroxyglutarate (2-HG), leading to a hypermethylation phenotype that dysregulates hematopoietic differentiation. In this study, we identified mutant R132H IDH1-specific gene signatures regulated by key transcription factors, particularly CEBPα, involved in myeloid differentiation and retinoid responsiveness. We show that treatment with all-trans retinoic acid (ATRA) at clinically achievable doses markedly enhanced terminal granulocytic differentiation in AML cell lines, primary patient samples, and a xenograft mouse model carrying mutant IDH1. Moreover, treatment with a cell-permeable form of 2-HG sensitized wild-type IDH1 AML cells to ATRA-induced myeloid differentiation, whereas inhibition of 2-HG production significantly reduced ATRA effects in mutant IDH1 cells. ATRA treatment specifically decreased cell viability and induced apoptosis of mutant IDH1 blasts in vitro. ATRA also reduced tumor burden of mutant IDH1 AML cells xenografted in NOD–Scid–IL2rγnull mice and markedly increased overall survival, revealing a potent antileukemic effect of ATRA in the presence of IDH1 mutation. This therapeutic strategy holds promise for this AML patient subgroup in future clinical studies.

Deregulated cell migration and invasion are hallmarks of metastatic cancer cells. Phosphorylation on residue Ser5 of the actin-bundling protein L-plastin activates L-plastin and has been reported to be crucial for invasion and metastasis. Here, we investigate signal transduction leading to L-plastin Ser5 phosphorylation using 4 human breast cancer cell lines. Whole-genome microarray analysis comparing cell lines with different invasive capacities and corresponding variations in L-plastin Ser5 phosphorylation level revealed that genes of the ERK/MAPK pathway are differentially expressed. It is noteworthy that in vitro kinase assays showed that ERK/MAPK pathway downstream ribosomal protein S6 kinases α-1 (RSK1) and α-3 (RSK2) are able to directly phosphorylate L-plastin on Ser5. Small interfering RNA- or short hairpin RNA-mediated knockdown and activation/inhibition studies followed by immunoblot analysis and computational modeling confirmed that ribosomal S6 kinase (RSK) is an essential activator of L-plastin. Migration and invasion assays showed that RSK knockdown led to a decrease of up to 30% of migration and invasion of MDA-MB-435S cells. Although the presence of L-plastin was not necessary for migration/invasion of these cells, immunofluorescence assays illustrated RSK-dependent recruitment of Ser5-phosphorylated L-plastin to migratory structures. Altogether, we provide evidence that the ERK/MAPK pathway is involved in L-plastin Ser5 phosphorylation in breast cancer cells with RSK1 and RSK2 kinases able to directly phosphorylate L-plastin residue Ser5.

The t(8;21) translocation is the most widespread genetic defect found in human acute myeloid leukemia. This translocation results in the RUNX1-RUNX1T1 fusion gene that produces a wide variety of alternative transcripts and influences the course of the disease. The rules of combinatorics and splicing of exons in the RUNX1-RUNX1T1 transcripts are not known. To address this issue, we developed an exon graph model of the fusion gene organization and evaluated its local exon combinatorics by the exon combinatorial index (ECI). Here we show that the local exon combinatorics of the RUNX1-RUNX1T1 gene follows a power-law behavior and (i) the vast majority of exons has a low ECI, (ii) only a small part is represented by "exons-hubs" of splicing with very high ECI values, and (iii) it is scale-free and very sensitive to targeted skipping of "exons-hubs". Stochasticity of the splicing machinery and preferred usage of exons in alternative splicing can explain such behavior of the system. Stochasticity may explain up to 12% of the ECI variance and results in a number of non-coding and unproductive transcripts that can be considered as a noise. Half-life of these transcripts is increased due to the deregulation of some key genes of the nonsense-mediated decay system in leukemia cells. On the other hand, preferred usage of exons may explain up to 75% of the ECI variability. Our analysis revealed a set of splicing-related cis-regulatory motifs that can explain "attractiveness" of exons in alternative splicing but only when they are considered together. Cis-regulatory motifs are guides for splicing trans-factors and we observed a leukemia-specific profile of expression of the splicing genes in t(8;21)-positive blasts. Altogether, our results show that alternative splicing of the RUNX1-RUNX1T1 transcripts follows strict rules and that the power-law component of the fusion gene organization confers a high flexibility to this process.

Chemotherapy-resistant human acute myeloid leukemia (AML) cells are thought to be enriched in quiescent immature leukemic stem cells (LSC). To validate this hypothesis in vivo, we developed a clinically relevant chemotherapeutic approach treating patient-derived xenografts (PDX) with cytarabine (AraC). AraC residual AML cells are enriched in neither immature, quiescent cells nor LSCs. Strikingly, AraC-resistant preexisting and persisting cells displayed high levels of reactive oxygen species, showed increased mitochondrial mass, and retained active polarized mitochondria, consistent with a high oxidative phosphorylation (OXPHOS) status. AraC residual cells exhibited increased fatty-acid oxidation, upregulated CD36 expression, and a high OXPHOS gene signature predictive for treatment response in PDX and patients with AML. High OXPHOS but not low OXPHOS human AML cell lines were chemoresistant in vivo. Targeting mitochondrial protein synthesis, electron transfer, or fatty-acid oxidation induced an energetic shift toward low OXPHOS and markedly enhanced antileukemic effects of AraC. Together, this study demonstrates that essential mitochondrial functions contribute to AraC resistance in AML and are a robust hallmark of AraC sensitivity and a promising therapeutic avenue to treat AML residual disease. Significance: AraC-resistant AML cells exhibit metabolic features and gene signatures consistent with a high OXPHOS status. In these cells, targeting mitochondrial metabolism through the CD36–FAO–OXPHOS axis induces an energetic shift toward low OXPHOS and strongly enhanced antileukemic effects of AraC, offering a promising avenue to design new therapeutic strategies and fight AraC resistance in AML. Cancer Discov; 7(7); 716–35. ©2017 AACR. See related commentary by Schimmer, p. 670. This article is highlighted in the In This Issue feature, p. 653

Background RNA sequencing (RNA-seq) and microarrays are two transcriptomics techniques aimed at the quantification of transcribed genes and their isoforms. Here we compare the latest Affymetrix HTA 2.0 microarray with Illumina 2000 RNA-seq for the analysis of patient samples - normal lung epithelium tissue and squamous cell carcinoma lung tumours. Protein coding mRNAs and long non-coding RNAs (lncRNAs) were included in the study. Results Both platforms performed equally well for protein-coding RNAs, however the stochastic variability was higher for the sequencing data than for microarrays. This reduced the number of differentially expressed genes and genes with predictive potential for RNA-seq compared to microarray data. Analysis of this variability revealed a lack of reads for short and low abundant genes; lncRNAs, being shorter and less abundant RNAs, were found especially susceptible to this issue. A major difference between the two platforms was uncovered by analysis of alternatively spliced genes. Investigation of differential exon abundance showed insufficient reads for many exons and exon junctions in RNA-seq while the detection on the array platform was more stable. Nevertheless, we identified 207 genes which undergo alternative splicing and were consistently detected by both techniques. Conclusions Despite the fact that the results of gene expression analysis were highly consistent between Human Transcriptome Arrays and RNA-seq platforms, the analysis of alternative splicing produced discordant results. We concluded that modern microarrays can still outperform sequencing for standard analysis of gene expression in terms of reproducibility and cost. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3819-y) contains supplementary material, which is available to authorized users.

The identification and significance of cancer stem-like cells in malignant gliomas remains controversial. It has been proposed that cancer stem-like cells display increased drug resistance, through the expression of ATP-binding cassette transporters that detoxify cells by effluxing exogenous compounds. Here, we investigated the 'side population' phenotype based on efflux properties of ATP-binding cassette transporters in freshly isolated human glioblastoma samples and intracranial xenografts derived thereof. Using fluorescence in situ hybridization analysis on sorted cells obtained from glioblastoma biopsies, as well as human tumour xenografts developed in immunodeficient enhanced green fluorescence protein-expressing mice that allow an unequivocal tumour-stroma discrimination, we show that side population cells in human glioblastoma are non-neoplastic and exclusively stroma-derived. Tumour cells were consistently devoid of efflux properties regardless of their genetic background, tumour ploidy or stem cell associated marker expression. Using multi-parameter flow cytometry we identified the stromal side population in human glioblastoma to be brain-derived endothelial cells with a minor contribution of astrocytes. In contrast with their foetal counterpart, neural stem/progenitor cells in the adult brain did not display the side population phenotype. Of note, we show that CD133-positive cells often associated with cancer stem-like cells in glioblastoma biopsies, do not represent a homogenous cell population and include CD31-positive endothelial cells. Interestingly, treatment of brain tumours with the anti-angiogenic agent bevacizumab reduced total vessel density, but did not affect the efflux properties of endothelial cells. In conclusion our findings contribute to an unbiased identification of cancer stem-like cells and stromal cells in brain neoplasms, and provide novel insight into the complex issue of drug delivery to the brain. Since efflux properties of endothelial cells are likely to compromise drug availability, transiently targeting ATP-binding cassette transporters may be a valuable therapeutic strategy to improve treatment effects in brain tumours.

MicroRNAs (miRNAs) are ubiquitously expressed small non-coding RNAs that, in most cases, negatively regulate gene expression at the post-transcriptional level. miRNAs are involved in fine-tuning fundamental cellular processes such as proliferation, cell death and cell cycle control and are believed to confer robustness to biological responses. Here, we investigated simultaneously the transcriptional changes of miRNA and mRNA expression levels over time after activation of the Janus kinase/Signal transducer and activator of transcription (Jak/STAT) pathway by interferon-γ stimulation of melanoma cells. To examine global miRNA and mRNA expression patterns, time-series microarray data were analysed. We observed delayed responses of miRNAs (after 24-48 h) with respect to mRNAs (12-24 h) and identified biological functions involved at each step of the cellular response. Inference of the upstream regulators allowed for identification of transcriptional regulators involved in cellular reactions to interferon-γ stimulation. Linking expression profiles of transcriptional regulators and miRNAs with their annotated functions, we demonstrate the dynamic interplay of miRNAs and upstream regulators with biological functions. Finally, our data revealed network motifs in the form of feed-forward loops involving transcriptional regulators, mRNAs and miRNAs. Additional information obtained from integrating time-series mRNA and miRNA data may represent an important step towards understanding the regulatory principles of gene expression.

Transcriptional profiling was performed on 452 RNA preparations isolated from various types of pancreatic tissue from tumour patients and healthy donors, with a particular focus on peritumoral samples. Pancreatic ductal adenocarcinomas (PDAC) and cystic tumours were most different in these non‐tumorous tissues surrounding them, whereas the actual tumours exhibited rather similar transcript patterns. The environment of cystic tumours was transcriptionally nearly identical to normal pancreas tissue. In contrast, the tissue around PDAC behaved a lot like the tumour, indicating some kind of field defect, while showing far less molecular resemblance to both chronic pancreatitis and healthy tissue. This suggests that the major pathogenic difference between cystic and ductal tumours may be due to their cellular environment rather than the few variations between the tumours. Lack of correlation between DNA methylation and transcript levels makes it unlikely that the observed field defect in the peritumoral tissue of PDAC is controlled to a large extent by such epigenetic regulation. Functionally, a strikingly large number of autophagy‐related transcripts was changed in both PDAC and its peritumoral tissue, but not in other pancreatic tumours. A transcription signature of 15 autophagy‐related genes was established that permits a prognosis of survival with high accuracy and indicates the role of autophagy in tumour biology., What's new? To date, relatively few RNA expression data are available on pancreatic tumour types other than pancreatic ductal adenocarcinoma (PDAC) or peritumoral tissues. Here, the authors found that transcriptionally, peritumoral tissues of PDAC behave much like the tumours, although resembling normal tissue in cell composition. In contrast, transcription in the environment of cystic pancreatic tumours is basically identical to that of normal tissue, suggesting an effect of peritumoral tissues on the strong pathogenic difference between PDAC and cystic tumours. At the functional level, autophagy‐related genes are uniquely changed in PDAC and its peritumoral tissue, with their expression variation predicting patient survival.

Animal models are essential tools for basic research and for preclinical therapeutic interventions. Although numerous clinical trials have been conducted, most fail due to inappropriate selection of the compounds at the preclinical stage. Therefore better preclinical models are crucial for achieving important clinical impact. Here we show that glioblastoma patient-derived xenografts (PDXs) based on organotypic spheroids represent a reliable and clinically-relevant animal model. Our PDX models are characterized by a tumor take rate close to 100% and a very reproducible tumor development time. We observe three distinct histological tumor phenotypes: a highly “invasive”, a highly “angiogenic” and an “intermediate” phenotype which combines invasion and vascular abnormalities. Typical glioblastoma characteristics such as pseudopalisading necrosis, invasion or microvascular proliferation are maintained. Interestingly, angiogenic features appear stronger in nude rats compared to immunodeficient mice.

The zebrafish has the capacity to regenerate its heart after severe injury. While the function of a few genes during this process has been studied, we are far from fully understanding how genes interact to coordinate heart regeneration. To enable systematic insights into this phenomenon, we generated and integrated a dynamic co-expression network of heart regeneration in the zebrafish and linked systems-level properties to the underlying molecular events. Across multiple post-injury time points, the network displays topological attributes of biological relevance. We show that regeneration steps are mediated by modules of transcriptionally coordinated genes, and by genes acting as network hubs. We also established direct associations between hubs and validated drivers of heart regeneration with murine and human orthologs. The resulting models and interactive analysis tools are available at http://infused.vital-it.ch. Using a worked example, we demonstrate the usefulness of this unique open resource for hypothesis generation and in silico screening for genes involved in heart regeneration. We thank C. Hoffmann and J. Esposito for technical support, A. Feenstra (VUA) for co-supervision of N.K., and R. Schneider (LCSB) for co-supervision of G.A. during her work at LIH and current support at LCSB. This research was funded by Luxembourg's National Research Fund (FNR) and the Swiss National Research Foundation (SNF), INFUSED project (www.infused-project.eu). N.M. was supported by Comunidad de Madrid (Fibroteam P2010/BMD-2321) and the ERC (Starting Grant 2013 337703 zebraHeart). Zebrafish drawing in Figure 1 was adapted from Wikimedia Commons (http://bit.ly/1PL92vj). Sí

While the autophagic process is mainly regulated at the post-translational level, a growing body of evidence suggests that autophagy might also be regulated at the transcriptional level. The identification of transcription factors involved in the regulation of autophagy genes has provided compelling evidence for such regulation. In this context, a powerful high throughput analysis tool to simultaneously monitor the expression level of autophagy genes is urgently needed. Here we describe setting up the first comprehensive human autophagy database (HADb, available at www.autophagy.lu) and the development of a companion Human Autophagy-dedicated cDNA Microarray which comprises 234 genes involved in or related to autophagy. The autophagy microarray tool used on breast adenocarcinoma MCF-7 cell line allowed the identification of 47 differentially expressed autophagy genes associated with the acquisition of resistance to the cytotoxic effect of TNFα. The autophagy-core machinery genes DRAM (Damage-Regulated Autophagy Modulator), BNIP3L (BCL2/adenovirus E1B 19 kDa interacting protein 3-like), BECN1 (Beclin 1), GABARAP (Gamma-AminoButyric Acid Receptor-Associated Protein) and UVRAG (UV radiation resistance associated gene) were found upregulated in TNF-resistant cells, suggesting a constitutive activation of the autophagy machinery in these cells. More interestingly, we identified NPC1 as the most upregulated genes in TNF-resistant compared to TNF-sensitive MCF-7 cells, suggesting a relation between the intracellular transport of cholesterol, the regulation of autophagy and NPC1 expression in TNF-resistant tumor cells. In conclusion, we describe here new tools that may help investigating autophagy gene regulation in various cellular models and diseases.

4; International audience; Epithelial to mesenchymal transition (EMT) is a key step toward metastasis. MCF7 breast cancer cells conditionally expressing the EMT master regulator SNAI1 were used to identify early expressed microRNAs (miRNAs) and their targets that may contribute to the EMT process. Potential targets of miRNAs were identified by matching lists of in silico predicted targets and of inversely expressed mRNAs. MiRNAs were ranked based on the number of predicted hits, highlighting miR-661, a miRNA with so far no reported role in EMT. MiR-661 was found required for efficient invasion of breast cancer cells by destabilizing two of its predicted mRNA targets, the cell-cell adhesion protein Nectin-1 and the lipid transferase StarD10, resulting, in turn, in the downregulation of epithelial markers. Reexpression of Nectin-1 or StarD10 lacking the 3'-untranslated region counteracted SNAI1-induced invasion. Importantly, analysis of public transcriptomic data from a cohort of 295 well-characterized breast tumor specimen revealed that expression of StarD10 is highly associated with markers of luminal subtypes whereas its loss negatively correlated with the EMT-related, basal-like subtype. Collectively, our non-a priori approach revealed a nonpredicted link between SNAI1-triggered EMT and the down-regulation of Nectin-1 and StarD10 through the up-regulation of miR-661, which may contribute to the invasion of breast cancer cells and poor disease outcome.

Small noncoding microRNAs (miRNA) regulate the expression of target mRNAs by repressing their translation or orchestrating their sequence-specific degradation. In this study, we investigated miRNA and miRNA target gene expression patterns in melanoma to identify candidate biomarkers for early and progressive disease. Because data presently available on miRNA expression in melanoma are inconsistent thus far, we applied several different miRNA detection and profiling techniques on a panel of 10 cell lines and 20 patient samples representing nevi and primary or metastatic melanoma. Expression of selected miRNAs was inconsistent when comparing cell line–derived and patient-derived data. Moreover, as expected, some discrepancies were also detected when miRNA microarray data were correlated with qPCR-measured expression levels. Nevertheless, we identified miRNA-200c to be consistently downregulated in melanocytes, melanoma cell lines, and patient samples, whereas miRNA-205 and miRNA-23b were markedly reduced only in patient samples. In contrast, miR-146a and miR-155 were upregulated in all analyzed patients but none of the cell lines. Whole-genome microarrays were performed for analysis of selected melanoma cell lines to identify potential transcriptionally regulated miRNA target genes. Using Ingenuity pathway analysis, we identified a deregulated gene network centered around microphthalmia-associated transcription factor, a transcription factor known to play a key role in melanoma development. Our findings define miRNAs and miRNA target genes that offer candidate biomarkers in human melanoma. Cancer Res; 70(10); 4163–73. ©2010 AACR.

International audience; The transcription regulator SNAI1 triggers a transcriptional program leading to epithelial to mesenchymal transition (EMT), providing epithelial cells with mesenchymal features and invasive properties during embryonic development and tumor progression. To identify early transcriptional changes occurring during SNAI1-induced EMT, we performed a time-resolved genome-scale study using human breast carcinoma cells conditionally expressing SNAI1. The approach we developed for microarray data analysis, allowed identifying three distinct EMT stages and the temporal classification of genes. Importantly, we identified unexpected, biphasic expression profiles of EMT-associated genes, supporting their pivotal role during this process. Finally, we established early EMT gene networks by identifying transcription factors and their potential targets which may orchestrate early events of EMT. Collectively, our work provides a framework for the identification and future systematic analysis of novel genes which contribute to SNAI1-triggered EMT.

Micro(mi)RNAs are small noncoding RNAs that orchestrate many key aspects of cell physiology and their deregulation is often linked to distinct diseases including cancer. Here, we studied the contribution of miRNAs in a well-characterized human myeloid leukemia, acute promyelocytic leukemia (APL), targeted by retinoic acid and trioxide arsenic therapy. We identified several miRNAs transcriptionally repressed by the APL-associated PML-RAR oncogene which are released after treatment with all-trans retinoic acid. These coregulated miRNAs were found to control, in a coordinated manner, crucial pathways linked to leukemogenesis, such as HOX proteins and cell adhesion molecules whose expressions are thereby repressed by the chemotherapy. Thus, APL appears linked to transcriptional perturbation of miRNA genes, and clinical protocols able to successfully eradicate cancer cells may do so by restoring miRNA expression. The identification of abnormal miRNA biogenesis in cancer may therefore provide novel biomarkers and therapeutic targets in myeloid leukemias.

Article, Background The reconstruction of context-specific metabolic models from easily and reliably measurable features such as transcriptomics data will be increasingly important in research and medicine. Current reconstruction methods suffer from high computational effort and arbitrary threshold setting. Moreover, understanding the underlying epigenetic regulation might allow the identification of putative intervention points within metabolic networks. Genes under high regulatory load from multiple enhancers or super-enhancers are known key genes for disease and cell identity. However, their role in regulation of metabolism and their placement within the metabolic networks has not been studied. Methods Here we present FASTCORMICS, a fast and robust workflow for the creation of high-quality metabolic models from transcriptomics data. FASTCORMICS is devoid of arbitrary parameter settings and due to its low computational demand allows cross-validation assays. Applying FASTCORMICS, we have generated models for 63 primary human cell types from microarray data, revealing significant differences in their metabolic networks. Results To understand the cell type-specific regulation of the alternative metabolic pathways we built multiple models during differentiation of primary human monocytes to macrophages and performed ChIP-Seq experiments for histone H3 K27 acetylation (H3K27ac) to map the active enhancers in macrophages. Focusing on the metabolic genes under high regulatory load from multiple enhancers or super-enhancers, we found these genes to show the most cell type-restricted and abundant expression profiles within their respective pathways. Importantly, the high regulatory load genes are associated to reactions enriched for transport reactions and other pathway entry points, suggesting that they are critical regulatory control points for cell type-specific metabolism. Conclusions By integrating metabolic modelling and epigenomic analysis we have identified high regulatory load as a common feature of metabolic genes at pathway entry points such as transporters within the macrophage metabolic network. Analysis of these control points through further integration of metabolic and gene regulatory networks in various contexts could be beneficial in multiple fields from identification of disease intervention strategies to cellular reprogramming., published version

The histopathological and molecular heterogeneity of glioblastomas represents a major obstacle for effective therapies. Glioblastomas do not develop autonomously, but evolve in a unique environment that adapts to the growing tumour mass and contributes to the malignancy of these neoplasms. Here, we show that patient-derived glioblastoma xenografts generated in the mouse brain from organotypic spheroids reproducibly give rise to three different histological phenotypes: (i) a highly invasive phenotype with an apparent normal brain vasculature, (ii) a highly angiogenic phenotype displaying microvascular proliferation and necrosis and (iii) an intermediate phenotype combining features of invasion and vessel abnormalities. These phenotypic differences were visible during early phases of tumour development suggesting an early instructive role of tumour cells on the brain parenchyma. Conversely, we found that tumour-instructed stromal cells differentially influenced tumour cell proliferation and migration in vitro, indicating a reciprocal crosstalk between neoplastic and non-neoplastic cells. We did not detect any transdifferentiation of tumour cells into endothelial cells. Cell type-specific transcriptomic analysis of tumour and endothelial cells revealed a strong phenotype-specific molecular conversion between the two cell types, suggesting co-evolution of tumour and endothelial cells. Integrative bioinformatic analysis confirmed the reciprocal crosstalk between tumour and microenvironment and suggested a key role for TGFβ1 and extracellular matrix proteins as major interaction modules that shape glioblastoma progression. These data provide novel insight into tumour-host interactions and identify novel stroma-specific targets that may play a role in combinatorial treatment strategies against glioblastoma.

Aberrant activation of oncogenic kinases is frequently observed in human cancers, but the underlying mechanism and resulting effects on global signaling are incompletely understood. Here, we demonstrate that the oncogenic FIP1L1-PDGFRα kinase exhibits a significantly different signaling pattern compared to its PDGFRα wild type counterpart. Interestingly, the activation of primarily membrane-based signal transduction processes (such as PI3-kinase- and MAP-kinase- pathways) is remarkably shifted toward a prominent activation of STAT factors. This diverging signaling pattern compared to classical PDGF-receptor signaling is partially coupled to the aberrant cytoplasmic localization of the oncogene, since membrane targeting of FIP1L1-PDGFRα restores activation of MAPK- and PI3K-pathways. In stark contrast to the classical cytokine-induced STAT activation process, STAT activation by FIP1L1-PDGFRα does neither require Janus kinase activity nor Src kinase activity. Furthermore, we investigated the mechanism of STAT5 activation via FIP1L1-PDGFRα in more detail and found that STAT5 activation does not involve an SH2-domain-mediated binding mechanism. We thus demonstrate that STAT5 activation occurs via a non-canonical activation mechanism in which STAT5 may be subject to a direct phosphorylation by FIP1L1-PDGFRα.

Major efforts have been put in anti-angiogenic treatment for glioblastoma (GBM), an aggressive and highly vascularized brain tumor with dismal prognosis. However clinical outcome with anti-angiogenic agents has been disappointing and tumors quickly develop escape mechanisms. In preclinical GBM models we have recently shown that bevacizumab, a blocking antibody against vascular endothelial growth factor, induces hypoxia in treated tumors, which is accompanied by increased glycolytic activity and tumor invasiveness. Genome-wide transcriptomic analysis of patient derived GBM cells including stem cell lines revealed a strong up-regulation of glycolysis-related genes in response to severe hypoxia. We therefore investigated the importance of glycolytic enzymes in GBM adaptation and survival under hypoxia, both in vitro and in vivo. We found that shRNA-mediated attenuation of glycolytic enzyme expression interfered with GBM growth under normoxic and hypoxic conditions in all cellular models. Using intracranial GBM xenografts we identified seven glycolytic genes whose knockdown led to a dramatic survival benefit in mice. The most drastic effect was observed for PFKP (PFK1, +21.8%) and PDK1 (+20.9%), followed by PGAM1 and ENO1 (+14.5% each), HK2 (+11.8%), ALDOA (+10.9%) and ENO2 (+7.2%). The increase in mouse survival after genetic interference was confirmed using chemical inhibition of PFK1 with clotrimazole. We thus provide a comprehensive analysis on the importance of the glycolytic pathway for GBM growth in vivo and propose PFK1 and PDK1 as the most promising therapeutic targets to address the metabolic escape mechanisms of GBM.